{ "cells": [ { "cell_type": "code", "execution_count": 4, "id": "fc68c9e1", "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "from imports import *\n", "from matplotlib.patches import ConnectionPatch\n", "np.set_printoptions(suppress=True)" ] }, { "cell_type": "code", "execution_count": 5, "id": "deaf42ab", "metadata": {}, "outputs": [ { "data": { "application/javascript": [ "IPython.OutputArea.prototype._should_scroll = function(lines) {\n", " return false;\n", "}\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%%javascript\n", "IPython.OutputArea.prototype._should_scroll = function(lines) {\n", " return false;\n", "}" ] }, { "cell_type": "code", "execution_count": 6, "id": "a74c43da", "metadata": {}, "outputs": [], "source": [ "source = '/home/sobolev/nevermind/Andrey/data'\n", "session = '008229_hippoSIT_2022-05-17_21-44-43'\n", "COLORS = ['tab:blue', 'tab:orange', 'tab:green', 'tab:red', 'tab:gray'] # for diff neurons\n", "\n", "animal = session.split('_')[0]\n", "sessionpath = os.path.join(source, animal, session)\n", "h5name = os.path.join(source, animal, session, session + '.h5')" ] }, { "cell_type": "markdown", "id": "54597948", "metadata": {}, "source": [ "## Head direction / unit" ] }, { "cell_type": "code", "execution_count": 9, "id": "254504be", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0.0, 0.5)" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA8wAAAFlCAYAAAA+pGFhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACrzUlEQVR4nOzdd3hUVfrA8e+ZSe+9J6RQQ4dQBSsqVuwFuyjWVX+6u9ZV17K77tq7KIodey9gQekl9A4hkEAgvfdk5vz+uAkghtSZzEzyfp4nD8nNnXvfkJM7973nnPcorTVCCCGEEEIIIYT4I5OjAxBCCCGEEEIIIZyRJMxCCCGEEEIIIUQLJGEWQgghhBBCCCFaIAmzEEIIIYQQQgjRAkmYhRBCCCGEEEKIFkjCLIQQQgghhBBCtMDN0QF0p7CwMJ2YmOjoMIQQQgghhBBC2MHq1asLtdbhtjper0qYExMTSU9Pd3QYQgghhBBCCCHsQCmVZcvjyZBsIYQQQgghhBCiBZIwCyGEEEIIIYQQLZCEWQghhBBCCCGEaIEkzEIIIYQQQgghRAskYRZCCCGEEEIIIVogCbMQQgghhBBCCNECl06YlVJeSqmVSqn1SqnNSql/OjomIYQQQgghhBA9g6uvw1wHnKi1rlRKuQOLlVI/aK2XOzowIYQQQgghhBCuzaUTZq21BiqbvnRv+tCOi0gIIYQQQgghRE/h0kOyAZRSZqXUOiAf+ElrvcLBIQkhhBBCCCGE6AFcPmHWWlu01iOAOGCsUmrI4d9XSs1USqUrpdILCgocEqMQQgghhBBCCNfj8glzM611KbAAmHrE9lla6zStdVp4eLhDYhNCCCGEEEII4XpcOmFWSoUrpYKaPvcGTga2OTQoIYQQQgghhBA9gksX/QKigbeVUmaM5P9jrfW3Do5JCCGEEEIIIUQP4NIJs9Z6AzDS0XEIIYQQQgghhOh5XHpIthBCCCGEEEIIYS+SMAshhBBCCCGEEC2QhFkIIYQQQgghhGiBJMxCCCGEEEIIIUQLJGEWQgghhBBCCCFaIAmzEEIIIYQQQgjRAkmYhRBCCCGEEEKIFkjCLIQQQgghhBBCtEASZiGEEEIIIYQQogWSMAshhBBCCCGEEC2QhFkIIYQQQgghhGiBJMxCCCGEEEIIIUQL3BwdgBBCCCGEEEKIP6ptsLA2u5RlmUVs2V9GVlE1pTUN1DZY8HI3E+rrQVKYL4OiAxifHMrIhCDczdIfamuSMAshhBBCCCGEE2iwWPl1Wz6fpO9l0c5C6hqtmBT0jfAjKcyXEF8PvNzN1DVayC+vY3tuBT9uzkVrCPZx5/Sh0Uwfl8DgmEBH/yg9hiTMQgghhBBCCOFA9Y1WPl29j5cWZJBTWkOEvyfTxyVwTEoYY5JCCPR2P+pry6obWJZZyPcbc/lszT7eX5HN+OQQbj6+L5P7haGU6safpOdRWmtHx9Bt0tLSdHp6uqPDEEIIIYQQQggsVs1Hq/YeTJSHxwdxy/EpnDgwArdODK8uq27go/Rs5izZw/6yWsYlhXDf6YMYHh9k++CdlFJqtdY6zWbHk4RZCCGEEEIIIbrX1gPl3PP5RtbvLWVUQhC3T+nPsTbqEa5rtDB35V5e+HUnhZX1XJQWx9+nDiTMz9MGkTs3SZi7QBJmIYQQwvVkF1WzfHcRRZX1WLUmJdyPkQlBRAZ4OTo0IYTosLpGC8/9vJNZCzMJ9HbnwbNSOXt4jF2GTlfUNvDCrxm8uXg33h5m7jq5P5eP79Op3mtXIQlzF0jCLIQQQrgGq1Uzf0seLy3IYGNO2Z++b1IwZVAkNx2fwsiEYAdEKIQQHZdXXsuN761mbXYpF46O477TBxHs62H382bkV/Lw15tZnFHIoOgAHp02mLTEELuf1xEkYT6MUioeeAeIBDQwS2v93NH2l4RZCCGEcH6ZBZX8/dMNpGeVkBjqw+Xj+3D8gHBigrwB2JFXyfzNucxdtZeS6nqunpjI304dgI+H1DIVQjivNdkl3PjuaipqG3nywuGcMSy6W8+vteaHTbk8+u0WDpTVcv6oOO45bSDh/j1rmLYkzIdRSkUD0VrrNUopf2A1cI7WektL+0vCLIQQQji3r9fv5+5PN+BuVjxwRirnjYo96tDByrpG/vvjNt5ZlsXAKH/evHrMwaRaCCGcyUersvnHl5uJDPTk9SvTGBgV4LBYqusbeeHXDN5YlImXe88bpi0JcyuUUl8BL2qtf2rp+5IwCyGEEM5Ja80zP+/k+V92ktYnmBenjyIqsH1zlBfuKOCW99fg42nmravHkhrjuBtRIYQ4XH2jlUe/3cK7y7OY3C+MFy4dSZCP/Ydgt8euAmOY9qKdPWuYtiTMR6GUSgQWAkO01uUt7SMJsxBCCOF8rFbNQ19v5t3lWVyUFsdj5wzFw61jPR3bcsu55q1V1DZY+OiGCfSP9LdTtELYTnFVPdtyy8kpqaGu0Yqnm4nkcD8GRPnj5ylTDFxdYWUdN7+3hpV7ipl5bDJ/P3WA0/XiHjlM+8SBEVw/OZnxySEuu36zJMwtUEr5Ab8Dj2utPz/iezOBmQAJCQmjs7Ky7BpLfaOVnNIaahsseLiZiA3yxsvdbNdzCiGEEK5Ka839X27igxXZ3HBcMvdMHdjpm7Q9hVVc9NoyNPDJDRNIDPO1bbBC2MDe4mo+Xb2PeZtz2ZZb0eI+ZpNiQnIoZw+PYdrIGDzd5F7S1WzcV8bMd9MprqrnvxcMY9qIWEeH1Krq+kbeWLSbt5fuoaiqniGxAVwzMYkpqZEEers7OrwOkYT5CEopd+BbYJ7W+unW9rVXD3NOaQ1frs1h3uZcth4op8Hyx//TlHBfJvUN49QhUYxPCsVkcs2nNUIIIYQtaa159NutvLlkNzcfn8Lfpw7s8jF35lVw0WvLCPLx4IubJzrN0EchduRV8OzPO/hxUy4aGJsYwrH9wxkRH0RcsDfeHmaq6yzszK9kTXYJP27KZXdhFRH+nsw8NpkrJyR2eOSFcIzP1+zj3s83EubnyWtXjGZIbKCjQ2q32gYLX6zN4fVFmWQWVOFmUkxICeXk1EjS+oTQP9LP6XrJjyQJ82GU8Qj6baBYa31HW/vbOmHeX1rDC7/u5JP0fTRaNaMSghiTFEK/CH98PczUNlrYW1zD6qwSVuwuorbBSp9QHy4f14fp4xLwlaE2QggherEn523nxQUZXD0xkYfOSrXZ8L/0PcVMf30Fo/sE886Msbg7+c2d6NmKKut45ucdfLAiG19PNy4b14crJ/Rps0Cd1pqlu4p4+bcMlmQUkRzuy8NnDebY/uHdFLnoqEaLlX//sI3Zi3czPjmEl6aPItTPNStQW62atXtLmb8ll/mb89hdWAWAl7uJxFBfUsL9iA32JjrQi5ggb2ICvYkO8iLU18PhQ7klYT6MUmoSsAjYCFibNt+ntf6+pf1tlTBrrXlnWRb//XEbDRbNJWPjuX5yMvEhPkd9TW2DhR835fLBymxW7i4m2MedGZOSuPqYJJmjIoQQotd5+bcM/vvjdi4ZE8+/zxtq8xusz9fs486P13Pp2AT+fd5Qmx5biPaob7TyzrI9PPfLTqrrLVw+LoE7pvTv1Jq7C7bl889vNrOnqJqL0+J54MxB+Hu51jDZnq64qp5bP1jD0l1FXD0xkfvPGNRjHtZprckurmZtdikbc8rILKhkd2EV+8tqqW+0/mFfTzcT0YFeRAd6G4l0kPH5gCh/hsYGdssoCUmYu8AWCXNFbQN/+2QDP27O5dj+4Tx+zpBWE+WWrMku4YVfdrJgewHBPu7ccFwKV07oI+tHCiGE6BXeXZ7FP77cxLQRMTx90QjMdpqq9MSP23jlt1389/xhXDQm3i7nEOJIWmt+2pLHv77fyp6iao7rH84DZwyiXxcL0dU1Wnj255289vsuYoK8eerC4YxLDrVR1KIrNu8v44Z3V5NfUcfj5wzhwrTecb3RWlNUVc+B0lpySms4UFbD/tIa9pfVcqC0hv2lteRX1GJtSje93c2M6hPEpL7hnDYkym51JiRh7oKuJswZ+ZXc8G46e4qqufe0gcyYlNSlJ+Lr9pby9E87WLijgDA/D246vi+XjUuwa5Ewq1Wzv6yGXQVVFFTUUVXXiMmk8HE3ExXoRXywD3HB3jLPWgghhF18tS6HOz5ax0kDI3jl8tF27YGxWDVXvrmC9D0lfHbTRJeaRyhc09YD5Tz23RaWZBTRN8KP+88YxAkDImx6jtVZxdz58Xqyi6u58bgU/m9Kf5nb7ECfrt7H/V9sJMjHndeuSGNEfJCjQ3IqDRYruWW1bMopY8XuYpZnFh0sdjcoOoDThkRx+tAo+kbYbmUDSZi7oCsJ8+b9ZVz+xgrMJsULl45iQortnuil7ynmqfk7WJZZRFSAF7ec2JeL0+JtcvHTWrP1QAVLMgpZnFHIqj3FVNdbWn2Nv5cbw+ICGZcUyjF9QxkWF9RjhpQIIUR3s1o1SuHwOV3O4Ndtecx8ZzWj+wTz9rVju2UViaLKOs58YTHuZhPf/GWSy1V7Fa4hv6KWp+fv4KP0vQR6u/N/U/ozfVyC3e6fquoaefTbLcxdtZchsQE8e/FI+kb42eVcomU19RYe/W4LH6zIZkJyKM9fOpJwf9ecr9zd9pVU8+OmXH7YlMvqrBIABkb5c9bwGM4cFk2f0K71PEvC3AWdTZg35ZRx+ewV+Lib+XDm+C7/Eo9m6a5Cnp6/g/SsEmKDvLntpL6cOzKuw4lzXaOF5ZnF/LQll5+35JNbXgsY1bonpoQxMNqfvuF+RAV64efphkVrquosHCirIauomk05ZazNLmVrbjlag5+nGxNSQjm2XxiT+4XTJ9RHbvyEEKIFtQ0WFmzLZ+HOQjbmlLK3uIby2gbMSuHv5UZ8iA/9I/0ZlRDMxJTQXnU9/X1HATPfSad/pD8fXD+uW+dfrs4q4eLXlnH8gHBmXZEmo6iEzdQ2WHhjUSYv/7aLBouVqyYk8pcT+xHo0z3te97mXO75bAM1DRbuP30Ql4/v02uuKY60KaeM2+euZVdBFTcel8JfT+nv9JWjnVVeeS0/bDzAtxsOkN6UPA+LC+SsYTGcMSy6zeJ4LZGEuQs6kzBv3FfGZW8sx9/LnQ+vH09CaMfmK3eU1pqFOwt5ev521u8rI9DbndOGRHHCwAjGJoa0WCiivtFKRn4lq/YYwxwW7Syksq4Rb3czx/YP46RBkUzuF0Z0YMcaXElVPcsyi1icUcjCHQXsK6kBID7Em8n9wpmYEsqw2CDiQ7zl4iwAo+DF0l2FbMwpY39pLTX1FrzcTcSH+DAwyp8JyaFEBHg5OkwhbK6oso7Zi3fz7rIsKuoa8fdyY3hcEMnhvgR5u9No1ZTVNLCnqIrtuRUUVtYDxvV06uAopg6JZmR8UI9N5BZsy+eGd1fTN8KP964bR0gnih511VtLdvPPb7Zw99SB3HR8SrefX/QsNfUW3l+RxayFmeRX1DF1cBT3nDbQIWt/55fX8rdPN/D7jgJOGBDOY+cOJbYTSYYzsFo123IrWJNdwp7CKvYUVbGvpIa6RisNFisWq8bHw0ywjwdxwd4khfkxKNqftMSQbrmuVNU18vwvO3lj8W7C/Dx46sIRTOoXZvfz9hY5pTV8t2E/36w/wMacMgDS+gRz5rBoTh8WTYR/++4hJWHugo4mzOv2lnLF7BUEehvJckeLe3VFc+L8xZp9zN+Sd3AYdZifJ5EBnvh4mKlvtFJa08C+khosTbPpY4O8ObZ/GCenRjIxJcxmw9201mQVVbNoZwELdxaybFcRlXWNAAR4uZEQ6kNUgBdRgV6E+nri4WbCbFK4NX1YNVi1bvow5pVprbFYje1aayxa4+1uJsTXkxBfd0J8PQ+WqrdXQRhhG0t3FfL20j38vDUfi1XjYTYRE+SFt4cbNfWN5JTW0GAxhqWOSgjmwtFxnDMytluGYwphT7UNFl74dSdvLt5DbaOF04dEM31cAuOSQo7a26C1JrOwiqUZhfy6LZ/FGYU0WDQR/p6cNCiSKYMiOKav7a7fjvbthv3830frGBgVwLszxjpsXWStNbd+uJYfNh7g/evG23RqlbCN3LJa1mSXsLuwityyWmobLFg1BPm4ExXgRd9IPwbHBLT7ptkeymoa+HBlNq8vzKSoqp4JyaH838n9GZsU4rCYwGjf7y7P4l/fbwXgxuNSuPG4FJe4juwtrj44dXDZriKKqowHip5uxvJF8SE++HiYcTMrzEpRXW+hqKqOvcU17C+roTmVSQ73Ja1PMOOSQpmQEtqpnsmjabBY+XT1Pp7/ZScHymq5OC2ee08fKOu829Gewiq+bUqet+dVoBQMigpgQkooY5NCGBwTQGxQy512kjB3QUcS5jXZJVw1eyVBvkayHBfcfcnykeoaLazfW8ba7BJ2FVRSVFlPVX0jnm5mArzdSQz1oW+EH6P7BHdbnA0WK1sPlLMpp5zN+8vIKa0ht6yW3PJaSqsbOnw8pcCk1MHE/3DuZkV8sA+JYb4MivYnNTqQwTEBJIT49NjeGFexI6+CR7/dwqKdhYT4enBhWhynDo5iWGzgH5KF+kYrO/Iq+HVbPt9tMC58wT7u3HR8CldOSHSJN3QhjrQ2u4S7PllPZkEVZw+P4baT+naqaEl5bQMLtuUzb3Muv28voKppZMbElDBGxgcxNC6QobGBLreWp9aa1xdl8q/vt5HWJ5jZV49x+PzhyrpGpr24mLKaRr6/bZKMeHEC2UXVfLZmH99u2M+ugqqD24N93PFuem8orWn4Q/2V5DBfJqSEcsrgKCYkh3ZLwautB8p5Z1kWX67NoabBwrH9w7ntxL6kJTo2UT7SvpJq/v39Nr7beIDoQC+un5zMxWPi8XWiJUyLKutYllnEkoxClmQUkV1cDUCEvyeT+oYxsW8Y45JCiA1quwhtTb2FjTllpGcVs3pPCauzSw7ehyaG+jAhJZTxyUYC3ZkHLcVV9Xy6ei/vLMtiX0kNIxOCeOCMQYzu41y/955uR14FP27KZemuQtZklx5cyirAy41B0QEkhfkSGeBFdKAXEQGenDQoShLmzmpvwrw6q5ir3lxFqJ8HH14/3qZPqHoDq1XTYDWGzTRYNFarxqQUygRmpYzPFZhNxuemw4rh1DZYKKmup6iynuKqenJKa9hTVEV2UTW7CirZVVB1MKn29TAzOCaQ4fGBDI8PYnhcEHHBMjy8O1itmpd/y+C5X3bi6+nGX07sx+XjE/B0azvx1VqzYncxL/+2i4U7CogO9OKe0wZy9vAY+d0Jl6C15r3lWfzzmy1EBnjxn/OHMrlfuE2OXddoYUVmMb9szWNRRiGZhyUQAV5uhPl5EurnQZCPB24t3EgqBW4mE+5mE+5mhZe7mRBfD0L9PIgO9KJfhH+7bkK7qqK2gfu+2MQ36/dzxrBonrpwuNM8GNuRV8G0F5cwNDaQD64fJ/MOHWRTThnP/bKTn7fmATA+KZSTBkUwNimElHC/PyV4JVX17MyvZN3eElZkFrMss4jqegv+nm6cOCiCqYOjOG5AuE2X6CyuqueHTQf4cm0Oq/aU4OlmYtqIGK6ckOj0FdeX7SrimZ93sHJ3MT4eZs4YGs3pQ6OZkBLa7X+LhZV1rGqqjrw8s5jteUaFZH9PN8anhHJMSiiT+oWREu7X5fuA5iHdyzKLWLariBW7i6ioNUZE9o3wY0Ky0TuZFOZLXLA3gd7ufzhnWU0DGfkVrM4qYeGOQpZlFmGxasYmhnDDccmcODBC7lUcrLbBwub95Ww9UM6WA8a/+0pqKKysOzjaIOuJMyVh7qz2JMyr9hRz9ZsriQjw4oPrx3V43q+wr9oGCzvzKtlyoIzN+8vZmGP82/ykKdTX42DyPDw+kOFxQS3O+xadV17bwP/NXccv2/I5Y1g0j5w9uNM9X8t2FfHvH7ayYV8Zk/qG8a9zh9q9ToAQXVHbYOGBLzfx6ep9nDgwgmcuHmHXXtOK2gbjWrfPGMlTWFlHUWU9JdX1WFt4/7ZqaLRYabBoGixWahosB28Wm/l6mOkX6W/M++sTwtikEJs+bFy6q5D7Pt9IdnE1d57cn5uP7+t0o4G+XGssbXXDscnce/ogR4fTq+RX1PLYt1v5ev1+Ar3duWJ8Hy4bn9Dh+63aBgtLMgqZtzmXn7fmU1xVj6ebiWP7hzN1cBRTBkV2qvBWfnktv+0o4PuNB1i8s5BGqyY53JdLxsRzUVq8yw3BXZNdwkcr9/Lthv1U1VvwdDMxOCaA1JgA+oT4Eh3kRZC3BwHebphNyvhQCqUOfa4xptPppml1HPy6eVod1FusFFXWUVhZR2FlPXnltezMq2Rn/qGaDT4eZtISQxifHML45NA/jUazB4tVs3l/GUt3GQn0kavF+HqY8fZwQymorG2kpuGwkQzhvpw6OIqzh8cwKDrArnGKrqtvtJJfUUtBRR2j+oRIwtxZbSXMKzKLuGbOKqICvfjw+vFEylAtl1DfaGV7bgXr9pWyfq/xkVFQefApU59Qn6YEOogR8UEMjglwmp4OV5ORX8nMd9PJLqrmwbNSucIG1TgtVs0HK7L474/bsWrNQ2cP5sLRcfIEVzidspoGrn87nZV7irn9pH7cflI/p0sEW1LfaG0asVPNjrxKtudWsCOvgk05ZZQ3JdNRAV6MSTKS57GJIfSL8Ovwz7ZlfzkvLcjgu40HiA/x5qkLRzh8XmdrHvhyI+8tz+a1K0Zz6uAoR4fT42mt+WxNDo9+u4WaBgszJycz87hkAmxQLb3RYmXVnhLmbc7lx0255JbX4mZSDI0LJK1PMENiA4kJ8iYmyJsQHw/MJkVto4WyaqMOzO7CKjbmGFPfmteHjQ3y5szh0Zw9PIbU6ACXf0+qbbCwYncxi3YUsGFfGVtzy//0MM2W/D3dSI7wY0Ckn7EyQJ9ghsYGOnyZ0gaLcc+4r6SGfSXV7C+tpbbRgtYaP083Qv08SQn3Y3h8oEPnyouukTnMXdBawrw0o5AZb6cTE2QkyzKvybVV1DawMaeM9XvLjCR6XykHyozltdxMioHR/geT6FEJwaSE+7r8m6G9rc4q4eq3VuLpZuKl6aMYl2zbgjk5pTXc9fE6lmcWM3VwFE+cP6zbluUQoi155bVc9eZKdhVU8vRFIzhreIyjQ+oyq1WzI7+ClbuLWbG7mFW7i8mvqAOMIktG73Mww+KC6BvhR6ivxx+uk80rNCzLLGLeplxW7jGGfl4/OZmbjnf+YkN1jRYuenUZmQVVfPOXSQ6pbtxbFFXWcdcn6/ltewFpfYJ54oJhpITbZ81gq1WzIaeM+ZtzWbm7mA37yqi3WNt8XYCXG8Pjg5iYEsZx/cMZFO3f4+8LyqobyC2vpaymgfKaBhqt+mCBVou1+V+aps4ZtWaAg1PrTEqhMKbVuZsVIb4ehPl5EubnibeHc//9i55NEuYuOFrCvHBHAde/k06fUB/ev268LDreQ+WV1x5MnpsT6YqmSt9xwd4cPyCcEwb0rOq0ttI8VSHc35P3rhtnt+JyVqvmjcWZ/G/edqICvXjlstFOP0/MVhosVtZml7JubwmZBVWUVjdQ12ghwNuoDpsaE8CohOBurdbfHjmlNWzdX052cfXBoWzhfp4kh/syJDawR/wt7S6s4orZKyiuque1K0bbbL6ys9Fak11czcrdxazcXcyqPcXsKao++H0PNxM+HmZ83M3UW4xe6+Y6jQOj/Dl7RAyXje3jUg+69pVUc8bzi4kJ8uaLmyf2iPbqbDbsK+Wm99ZQUFnHfacN5MoJid06MqO2wcLe4mpySms4UFZLSXU9FovGy91MoLc7MUHe9An1kRooQvQgkjB3QUsJ86/b8rjx3TWkRPjx3oyxLleFVHSe1arZXVTFsl1F/La9gKW7Cqmut+Dn6cYpqZGcNTyGSf3CHD58yNGWZxZxbdNUhQ+uG09UoP1HX6zJLuGW99dQVFXPI2cP5uIx8T3yRkZrzcrdxXyUvpf5m/MOLtUW2lSkycPNRHlNI7nltQfn6Q+I9Of0odFcNCbOYTUWMvIr+Dh9H79uyycjv/Ko+3mYTYxJCuasYTGcOTwGPyeq0tpem3LKuOrNlWhgzjVjGBYX5OiQulV+eS1bcyvIyK8kv8JYW7263oKHm4kwXw/6RfozIj7I6R7kdMSCbflcM2cV00bE8OzFI3rktcZRPk7fywNfbiLcz5NXLx/N0Lje8QBUCOFYkjB3wZEJ8yfpe7nn842kRjt2bUjhHJqr03634QA/bDpAeW0jEf6eXDwmnkvGJhDbC6ulL8koZMbbq4gP9uH968d163yeoso6bp+7jsUZhVw4Oo7Hzh3SrircrmLRzgKenLed9fvK8PN04/ShUZw4MJKxSSGEHFGorsHSNPR1VxE/bs5l1Z5iTEpx6uBIbj6+b7f1wm89UM6Lv2bw/aYDuJkU45JCOWFgBCPig0gM9cHPyw2toaCijq0HyknPKmH+5lz2FFXj7+nGpeMSuG5yksvMC1uRWcSMt9MJ9Hbn3RljSbbTEFLheC8tyOB/87Zz+0n9+L+T+zs6HJenteb5XzJ45ucdHNM3lBcuHfWn65oQQtiLJMxd0JwwW62a537ZyXO/7GRyvzBevmwU/jYoOiF6jrpGC79tL+CjVXtZsD0fBRw/IILLxydwXP8IzC5Q6Kerft9RwMx30kkK8+W968YR5oDRFxar5tmfd/DCrxmk9Qnm1StGOyQOW9pbXM2j325h/pY84kO8ufG4FM4bGdeh+V7ZRdW8vyKLD1ZmU1HbyEkDI7h9Sj+79X5uyinj+V92Mn9LHn6eblw1sQ8zJiW36wZYa83avaW8uXg3P2zKxcNsYsakJJsV+7GXBdvyufG91cQFe/PedbJiQk+nteZvn27g09X7eObi4Zw7Ms7RIbksq1Xz6HdbeGvJHs4fFccT5w+VpbuEEN1KEuYuSEtL0wsWL+Ouj9czf0se54+K49/nDe2WBe+F69pXUs3clXv5KH0vBRV1xAV7M31cAhenxffYIfwLtuVzw7ur6Rvhx3vXjXN4z8C3G/Zz18frCfPzZPbVaQyMcr3lHSxWzSu/ZfD8rxmYleLWE/ty3eSkLvWal9U08M7SPcxespvS6gamDIrgjin9bdbjvDa7hBd+zeDXbfkEeLlx7aQkrpmY1Ok5qnsKq3jqpx18s34/YX6e/OPMQU65/vY36/fzfx+tY2C0P29fI1N1eov6RitXvbmS9Kxi3psxzuaFDXuDRouVuz/byGdr9nHNMYn844xUl6gkL4ToWSRh7oLUYSN1xBVPs7ekhgfOGMTVExOd7kZNOK8Gi5X5m/N4b3kWyzKL8DCbOH1oFFdM6MOohOAe05Z+2pLHze+vZmCUc01VWL+3lOvfSaeqrpHnLx3JSYMiHR1Su+0rqeb/PlrHqj0lnD40in+cmWrTHsuK2gbmLNnD64syKa9t5OTUSO6Y0o/BMR1PnLXWLNtVxCu/72LRzkKCfdy5bnIyV0zoY7Me4Q37Snngy01s2FfG5H5hPHbOEPqEOkeF4g9XZnPfFxsZ0yeEN65Oc+pecGF7ZdUNnPvKEoqr6vlo5gQGRPk7OiSXUdtg4bYP1zJ/Sx53ntyfv5zYt8e8LwohXIskzF3gFd1Pj7rtVZ65eIQ8ORZdsjOvgvdXZPPZ6n1U1DUyKDqAy8cncM6IWHxdsLBRs+82HOD2uWsZHBvIO9eOJdDbuZKF3LJarntnFZv3l3P/6YOYMSnJ6W/Ilu0q4qb3V9No0Tx6zmDOGRFrt5jLaxt4a/Ee3licSUVtIycMCOeitHhOHBTRZk92dX0jP2zM5Z1le1i/r4wwP0+um5zEFeP72KVNW6ya95Zn8b9522m0WrljSn9mTEpyaJG9WQt38a/vt3H8gHBeuWy0LIvSS2UXVXPBq0uxavjkxgkkyXJTbaqsa2TmO+ks3VXEw2elcvUxSY4OSQjRi0nC3AWRyal6x+b1TpcECNdVVdfI1+v38+6yLLYcKMfP043zR8Vy+fg+9It0rZ6JL9bu466P1zMqIZg3rxnjtD1r1fWN3PnRen7cnMulY+N5ZNoQp61k/sGKbB78ahOJYb68cWVat63zWlbTwFtLdvPhymzyyusI8nHn2H7hjEkMJjHMlyBvY9RAQWUtmQVVrNhdzNKMQqrqLSSG+nD9scmcPyquW5bYyS2r5aGvNzFvcx6DogN44vyh3V6JutFi5dFvt/D2sizOGBbNMxeNkKk6vVxGfgUXv7YcTzcTH984wW5L6fUEJVX1XP3WSjbtL+fJC4fJ/G8hhMNJwtwFR1uHWYiu0lqzJruU95Zn8d2GA9RbrKRGB3DK4EhOHRzFwCh/p+4JfXvpHh7+ZjMTkkN546o0fDycu5fcatU8OX87L/+2i4kpobx82SinGToORgL22HdbmbN0D8cPCOf5S0c65AGExapZnFHIl2tzWJJRSH5FXYv79Qn14Zi+YZwzIpYxiY6ZXvDjplwe/GoThZV1XHNMEnee3L9bRmuU1TRw+9y1/La9gOsnJ3HPaYN6RVE/0bYt+8u5ZNYygnw8+OTGCUQGuEZ19+6UW1bLFbNXkFVczcvTRzEl1XWmygghei5JmA+jlHoTOBPI11oPaWt/SZhFdyiqrOOLtTn8uCmX1dklaA3h/p4ckxLKxL5hHNM3zGmWqLJaNf/5cRuzFmZycmokL1w6slt6FW3ls9X7uOfzDcQH+zD76jFOMXSyrLqBWz9cw6KdhcyYlMR9pztHAqa15kBZLdnF1VTUNqK1JtTPkz6hPk5Teby8toEnftjG+yuyiQ3y5vFzh3D8gAi7nW/jvjJu/mA1B0preWTaEKaPS7DbuYRrWptdwuVvrCDc35N3Z4xz6fWmbW1PYRWXz15BaXUDr1+ZxoQUmeomhHAOkjAfRil1LFAJvCMJs3BG+RW1LNiWz5KMIpbuKqSwsh6AhBAfxieHMC4plPEpoQ5JoEur67njo3X8tr2AK8b34eGzBztFYtdRK3cXc8O76Vg1vHr5aIfetGUWVHLd2+nsLanmsXOGcPEYScA6Y9WeYu79fCMZ+ZVMGxHDP85MtWlSb7FqZi/O5Ml5Owj18+DF6aMY3SfYZscXPcua7BKueWsVXu4m3psxzuWm29jDlv3lXPnmSixWK29fO7bbp1EIIURrJGE+glIqEfhWEmbh7LTWbM+rYElGEcszi1i5u5iymgYA4oK9jeQ5OYTxyaHEBXvbdVjs6qxibp+7jrzyWh48M5XLx/dx6iHjbckqqmLG2+nsKazi8XMdk6gu3lnIze+vxs1s4pXLRklhwS6qa7Twym+7eGlBBj4ebsw8NpmrJibi18Vh2muzS3j4my2s31vKyamRPHH+MIcvmyac3/bcCq6YvYJ6i5W3rh7DyITe+4Bl5e5iZry9Cj9PN96dMY6+EX6ODkkIIf5AEuYjtJUwK6VmAjMBEhISRmdlZXVjdEIcndVqJNDLM4tYkVnMit1FlFQbCXRskDfjkkIY15RAJ4T42CShLa9t4On5O3h72R5iAr15cfrIHnPjV17bwC3vG0OhZx6bzN1TB3ZLj7nWRrXnh7/ZQkq4L7OvGiPDNm0oI7+Cx7/byoLtBQeXuLp4THyHe5w37Cvltd8z+W7jAcL8PHnwrFTOGhbt0g+KRPfKLqrm8tkrKKys46XLRnGCHacLOKv5m3P5y4driQ325t0Z45xmepEQQhxOEuYjSA+z6CmsVs3O/Eojgd5tJNFFVcYQ7qgAL2MId3Io45JCSArz7dCNfklVPR+uymbWwkxKqxu4akIf/jZ1YJd765xNo8XKP7/ZwrvLs5gyKJLnLhlh18JR9Y1WHvp6Mx+uzOakgRE8e8kI/J20urirW5tdwrM/7+T3HQW4mxXHD4hgcr8wJqaEkhLu96e/h0aLlYyCShbuKOC7jbms31uKr4eZGZOTmXlsco9r+6J75JfXcvVbq9iWW869pw3iusnOv7SdLWiteX1RJv/+YRvDYgN58+oxhDpJ7QMhhDiSJMxHkIRZ9FRaazLyK1m+u5gVmUUszyymsNKochzm58nIhCBGJgQxKDqAxFBfYoO88XAzobWmsq6RrKJq1u0tZdHOAhZsK6DeYuX4AeH89ZQBDIkNdPBPZ19vL93DP7/ZzICoAGZflUaMHXpBCirquOm91aRnlXDz8SncdcoAl5wD7mp25lXw4cq9zNucS05pDQD+nm7EBnvj7+WGUorS6nqyi6upbbACMDgmgPNGxXFhWpzTLpcmXEd1fSN3fbyeHzblcsHoOB4/d0ib65y7sqq6Rh74chNfrM3h9KFRPHXhCFmjXAjh1CRhPoIkzKK30FqTWVjF8swiVu8pYe3eUnYXVv1hH3ezwmLVWA/7s44O9OLUwVFcPCaeQdEB3Ry14/y2PZ9bP1iLt4eZVy8fbdOiTot2FvB/H62nsq6B/10wnLOGx9js2KJ9tNbsLa5h6a5Cth4oZ19JDdX1FixaE+zjTmyQD0PjAkjrEyJD5IXNWa2a53/dybM/72RkQhAvTh/VI4cnbz1Qzq0frGF3YRV3TOnPrSf0xSQPBoUQTk4S5sMopT4EjgfCgDzgIa317KPtLwmz6GmKq+rZVVBJVlE1+0uNhMHdrPDzdCMhxIcBUf4dHr7dk+zIq2DG26vIKanh+mOT+b8p/bu0bFZ1fSNPz9/BG4t30y/Cjxenj2JAlFTMFaK3+mHjAf726QZMCv57wTCmDol2dEg2Uddo4eUFRtG9YF8PnrtkBBNTwhwdlhBCtIskzF0gCbMQvU95bQP/+m4rc1ftJTnclwfPTOW4/uEdeoigtWbe5lwe/XYrOaU1XDYugQfOSJVhiUIIsoqquO3DtazfV+by1waLVfPthv38b9529pXUcO7IWB48M5VgqSQvhHAhkjB3gSTMQvRei3YWcP8Xm8gurmZ0n2BmTEpiyqBIPNxMR31NbYOF+VvymLVwF5tyyukX4ce/zhvKmMSQboxcCOHs6hutPDV/O68tzCQu2JtHzxniUlW0Cyvr+Grdft5euofs4mpSowO47/RBTOonvcpCCNcjCXMXSMIsRO9W32hl7qpsXl+Uyd7iGvy93Dh+QAQj4oOIDfLC19ON6noLOSU1rM4uYdGOAsprG0kI8eH2k/oxbUQMbuajJ9hCiN5tRWYR932xkV0FVZwxNJoHz0olMsDL0WEdZLFqSqvrKaqqJ6e0hlW7i1mxu5h1e0uxWDWjEoK4bnIyUwdHyVxlIYTLkoS5CyRhFkKAcdNoLDd0gEU7C8grr/vTPjGBXoxPDuW8UXFMTAmVm0chRLvUNVqY9XsmLyzIwKRg+tg+3HhcMhF2Spy11hRU1pFfXkdeeS35FXUUVdZRWGkkxkWVdRRV1lNUVUdxVf0fikKaTYqhsYFM6hvGtBEx9IuUmgxCCNcnCXMXSMIshGhJQUUdBRV1VNY14uNhJiLAkwh/5+kVEkK4nqyiKp7/JYMv1+VgNinOHxXLRWnxjIgP6lQhRq01OaU1bD1QwY68CjLyK8nIr2RXQSXV9ZY/7e/v6Uaonwehfp6E+hr/hvl5HPw83N+TobGBdl2rXgghHEES5i6QhFkIIYQQ3SmrqIpXftvFV+v2U9NgIT7EmxMHRDA6MYRBUf7EBfv8oUiY1Wr0GOeU1pCRV8mWA+Vsbfoor208uF9UgBf9Iv1ICfcjOdyXyAAvIvw9iQjwItTXo0srAgghhCuThLkLJGEWQgghhCNU1Dbw3YYD/LQlj6W7iqhpONQr7OlmwsvdTKPFStURvcU+HmYGRPkzKDqAQdEBpEb70z/SH38v9+7+EYQQwiXYOmGWcThCCCGEEHbm7+XOJWMTuGRsAg0WK9tzK9hVUMm+khrKaxqobbDgbjbh4+lGuL8nsUFeJIb6khjqKzUUhBDCgSRhFkIIIYToRu5mE0NiAxkSG+joUIQQQrRB1kcRQgghhBBCCCFaIAmzEEIIIYQQQgjRAkmYhRBCCCGEEEKIFkjCLIQQQgghhBBCtEASZiGEEEIIIYQQogWSMAshhBBCCCGEEC2QhFkIIYQQQgghhGiBJMxCCCGEEEIIIUQLJGEWQgghhBBCCCFaIAmzEEIIIYQQQgjRAkmYhRBCCCGEEEKIFrh8wqyUmqqU2q6UylBK3ePoeIQQQgghhBBC9AwunTArpczAS8BpQCpwqVIq1bFRCSGEEEIIIYToCVw6YQbGAhla60ytdT0wF5jm4JiEEEIIIYQQQvQArp4wxwJ7D/t6X9O2g5RSM5VS6Uqp9IKCgm4NTgghhBBCCCGE63JzdAD2prWeBcwCUEpVKKW2OzgkIWwhDCh0dBBC2IC0ZdETSDsWPYW0ZdETDLDlwVw9Yc4B4g/7Oq5p29Fs11qn2TckIexPKZUubVn0BNKWRU8g7Vj0FNKWRU+glEq35fFcfUj2KqCfUipJKeUBXAJ87eCYhBBCCCGEEEL0AC7dw6y1blRK3QrMA8zAm1rrzQ4OSwghhBBCCCFED+DSCTOA1vp74Pt27j7LnrEI0Y2kLYueQtqy6AmkHYueQtqy6Als2o6V1tqWxxNCCCGEEEIIIXoEV5/DLIQQQgghhBBC2IXLJ8xKqTeVUvlKqU0tfO8upZRWSoU1fa2UUs8rpTKUUhuUUqMO2/cqpdTOpo+ruvNnEAI63JYHKqWWKaXqlFJ/PWLfqUqp7U3t/J7uil8I6HA7vqzpWrxRKbVUKTX8sH2lHQuH6mBbntbUltcppdKVUpMO21fuL4TDdKQdH7Z9jFKqUSl1wWHbpB0Lh+rgNfl4pVRZ0zV5nVLqwcP27fD9hcsnzMAcYOqRG5VS8cApQPZhm08D+jV9zAReado3BHgIGAeMBR5SSgXbNWoh/mwO7W/LxcBtwJNH7GsGXsJo66nApUqpVDvFK0RL5tD+drwbOE5rPRR4lKY5R9KOhZOYQ/vb8i/AcK31COBa4I2mfeX+QjjaHNrfjpuvv08A8w/bJu1YOIM5dKAtA4u01iOaPh5p2rdT9xcunzBrrRdiJA9Hegb4O3D4JO1pwDvasBwIUkpFA6cCP2mti7XWJcBPtPALEcKeOtKWtdb5WutVQMMR+44FMrTWmVrremAuRrsXolt0sB0vbbrmAiwH4po+l3YsHK6DbblSHyoK43vY9+T+QjhUB++TAf4CfAbkH7ZN2rFwuE605ZZ06v7C5RPmliilpgE5Wuv1R3wrFth72Nf7mrYdbbsQDtVKWz4aacvC6bSzHc8Afmj6XNqxcEqttWWl1LlKqW3Adxi9zCBtWTiho7VjpVQscC5NIzAPI+1YOKU27i8mKKXWK6V+UEoNbtrWqbbs8stKHUkp5QPch9E1L4TLkrYseoL2tGOl1AkYCfOko+0jhKO11Za11l8AXyiljsWYYjClG8MTol3aaMfPAndrra1KqW6NS4iOaqMtrwH6aK0rlVKnA19iTMntlJ7Yw5wCJAHrlVJ7MIb4rVFKRQE5QPxh+8Y1bTvadiEcqbW2fDTSloWzabUdK6WGYcz3nKa1Lmp6jbRj4YzadU1uGjaY3FR8RtqycDatteM0YG7T9guAl5VS5yDtWDino7ZlrXW51roSQGv9PeDelWtyj+th1lpvBCKav276D0zTWhcqpb4GblVKzcUoXFCmtT6glJoH/OuwAganAPd2c+hC/EFrbbmVl60C+imlkjAuAJcA0+0ZpxCtaeOanAB8Dlyhtd5x2MukHQun00Zb7gvs0lprZazA4QkUAXJ/IZxKG/cWSYdtnwN8q7X+sqnol7Rj4VTauCZHAXlN1+SxGJ3ERUApnbi/cPkeZqXUh8AyYIBSap9SakYru38PZAIZwOvAzQBa62KM4VOrmj4eadomRLfpSFtWSkUppfYBdwIPNO0foLVuBG7FuEnbCnystd7cHfELAR2+Jj8IhGL0YqxTSqUDSDsWzqCDbfl8YJNSah1GBdaLmwqMyv2FcKgOtuMWSTsWzqCDbfkCjGvyeuB54JKma3Kn7i/UoaKOQgghhBBCCCGEaObyPcxCCCGEEEIIIYQ9SMIshBBCCCGEEEK0QBJmIYQQQgghhBCiBT2uSnZrwsLCdGJioqPDOKrKukb2FlfTJ9QHH4/WfzUVtQ3sKaoGYGhsYHeEJ4QQLqe2wUpZTT3lNY1o2q7ZoZSiT4gPHm7yPFkIIZxJWU0D2cXGva/ZpIgJ9CbIx93BUQlntHr16kKtdbitjterEubExETS09MdHcZRvb8ii/u/2ITJ041ZM8YyKiG4xf025ZRx4avLiG6wEOjtTvpDLa3XLYQQvY/WmqveWsWB0hrKahqorKjDXcEZyaGE+nm2+trS6noW7SzkmavSOGlQZDdFLIQQoj2+3bCfWz9Yy9MXDeedZVms21vKyNRI/nXuUML9W7++i95FKZVly+P1qoTZ2TUXLPdyN3PV7JW8e904RsQH/Wm/v3y4lhBfD/pH+rH1QEX3BimEEE6s3mJl4Y4CAM4bFcuohGBOSY0kIsCrzddu2FfKop2FyOIRQgjhfJqvzUNjA/nsponMXpzJk/N3cNpzi3jm4uFM7mezDkUh/kASZic0+6o0/vLhWqa/vpwHzkjlorQ4zCbFjrxK1maXsLuwirtO7s++khqnSpi11mzKKWdHXgWNVivh/p70i/AnLtgbpZSjwxNC9ALNN1R/O3UAt5zQ17HBCLsrq26gsKqO8poGGq2tP+nw83RjUHRAN0UmhLAns0kx89gUju0fzl8+WMsVs1dyw7HJ3HXKAJlSI2xOEmYn0vxWHx3kxcc3TOCOj9Zy3xcb+ff3W3EzK0qqGw7uOzA6gH0lNe2ak9cd1maXcPdnG9iRV/mn78UFe3Nc/3BOGxLN+OQQ3MxyIRNC2FdnntEpjBc5x1VVtKS6vpFv1u/npy35rN9XSkFFXYdeP++OYxkQ5W+n6IQQ9tR8bT78+j4wKoCvb53Eo99t4bWFmSzLLOL5S0aSGObrkBhFzyQJsxNSKKICPfnguvH8ui2f33bkY7FqhscFMT45FH8vN0L9PPlla56jQwVgdVYxl76+gnA/T/57wTBG9wnG081EblktWw6Us2hnIV+szeH9FdmE+nowdUgUZwyLZlxSKGaT9DwLIYRondaaj9P38sSP2ymuqj/4IHZApD8RAZ74e7nhYTYf9fXbcst57LutFFbWMQBJmIXoSbw9zPzr3KEc2y+Mv3+6gTOeX8Tj5w7lnJGxjg5N9BCSMDuTIybOmUyKKamRTEk9evEZR8+1q2+0cufH64kO9OKLm48hxNfj4Pfign1ISwzhygmJ1DZY+G17Pt9uOMDna4zkOczPk9OakucxiSGSPAshbKa5t7hDr2l6iXb0hVX8gdWqeeCrTXywIpuxSSH87dQBpPUJ7tBUn+ZKulV1jfYKUwhhZ4euzS3/7U8dEs3QuCD+b+467vhoHQt3FPDIOUPw85R0R3SNtCAn1N57AGeYFvz5mn1kFVUz55oxf0iWj+TlbmbqkGimDommpt7Cgu35fLfhAJ+s3su7y7MI9/fk9CFRnDk8htEJwZgkeRZCdILkuj3Pf+dt54MV2dx4XAp/P3VAp94fvD2M3ueaBoutwxNCOJHYIG8+uH4cLy7I4PlfdrImu4TnLx3JsLggR4cmXJgkzC7O0feGX6zNISXcl+P6t78yobeHmdOHRnP60Giq6hr5dZuRPM9dtZe3l2URGeDJlRMSuX5yshRuEEJ0SlceKDr6uioOWbijgFd/38X0cQncPXVApwtIhvh44GZSvLxgF8E+HkzuFybFKIVwUW396bqZTdwxpT8TU8K4Y+5azn9lKX8/dSAzJiVJh4zoFEmYnUjrA01a4tg/+tLqelbuKea2E/t1+sbD19ONs4bHcNbwGCrrGvllax6fr8nhf/O2M39zLrOvHkNYG2unCiFEM2cphCi6rq7Rwv1fbqRvhB8PnpnapQQ32NeDVy4fzYNfbeLKN1cSH+LN5H7GHOij3UB7mk2cNTzmYO+0EMK1jE0K4fvbJ3PPZxt5/PutLNxZwFMXDSfCv+1lBoU4nCTMLs6Rww+Lq+rRGpLDbVOJ0M/TjWkjYpk2IpYfNh7g/z5ex6WzlvP5zRPx93K3yTmE88guqua9FVlsPVCOm0kR4O3e4iOguGAf7jqlv/QGiQ7pTGs5NIfZpqGITnp/eTZ7i2t4+9qxeLl3PWk9OTWSyf3C+G7DAX7YdICv1uZQVd/6EG03s+K8UXFdPrcQouuar80dub4H+XjwyuWj+GBlNo98s4XTn1vMB9ePo3+kFP8T7ScJsxM5eCFoZ2Lg6PyheS6YLW5kjnTa0GgCvN258s2V3PfFJl64dKTNzyEc68GvN7F4ZyGDYwJotGoyC6v+tE95TQMl1Q1cf2wygd7y0ES0TZLdnmP24t1MSA7l2H5hNjuml7uZ80fHcf7oOBotVkprGlpsMw0WKxP/8ytfrM0BoF+EP+H+ng5/3z2Sl5uZQB+5NgrRGqUUl43rQ1qfEK6YvYLpry/noxsmkBLu5+jQhIuQhNnlOe7usLYpYfa2Q8IMcEzfMG47sR/P/LyD80bFcsKACLucRzhGZW0jY5NC+OD68Ufd591le/jHV5tpsFi7MTLRE3RlHWaZxewciqrqOHNYtN1Gl7iZTa1O+ZkyKJLfd+SzaGehXc5vCyZlrC3dT3rLRC/QPOWms9eEAVH+fDhzPBe/tozr3k7nq1uPIUBGMIp2kITZiTSXy2/vZcDRD7rfW54NYNf5XTcen8xX63N4+OvNTLgj1C692cIx6i1W/LxavwQ1F32rb5SEWbTPoVoQjr5Ciq6wWjW1DVaHXvPfuCqNBouVrKIqduRVUlrd4LBYWrKnqIpZCzMpqKiThFmIdkoJ9+OVy0dz6azl3PnROl6/Mk2mfIk2ScLs4hw1/LCituHgULX4YB+7ncfTzcwjZw/h8tkrmPrsQmZMTuaK8X3sdj7RfeobrXi2UQXd3Wx8X3qYRXt1ZQ1lmcPsPGobm0YwObjglrvZRN8If/pGOF9CumpPMbMWZmJxkga7aGcByzOLGBITyOCYwHb97gK83fB0kwfhon06M4e5JWMSQ7jv9EE88u0WPlm9j4vS4rscm+jZJGF2Igd7RlxgHebs4moAXr5sFFGB9q02OKlfGA+cMYgXF2Tw46YDkjA7mf2lNTz67RZmHpvMyITgdr+urtGKRxs3SpIwi86SDgPXZLFqnvlpB/tKjPcYe0356QlMTY3c6hz5Mi/8msHK3cUdek1KuC+/3HW8fQISohVXT0zkx825PPrtFo7rH05kgFTOFkcnCbOLc9T7ZG5ZLQDRdk6Wm103OZlFOwsprXGuIXECNuWU8cOmXOZtzuX6Y5O5ZmJSu9bPrm2w4GFuXw9zfaOT3BEKp9eVlnKwh9kmkYjO+HbDfl5ckEG4vycJIT4MjQt0dEhOq3k1LGs7MubVWcW8vzybXQWVduuR3plXydTBUdx4fArbDpTT0EZcS3YW8uPmXG56bzVKNU2jUMZDEj9PN3w8zPh6uhHh70limC99w/0I9vWwS+zCNRwqjtv1Y5lMiv+eP4xTnl3If37YxjMXj+j6QUWPJQmzEzk01KSdVbIdOEcvt7w5YfbutnP6ebkd7HUQzqP5nmhCSiiv/Z7Ja79ntvu1vp6t9x41D9mWHmYher7q+kb+++N2BkUH8N1fJh11fWRhMDf9/1jaSEwr6xq58NVlBHi7MzQ2sM0HlZ0VFeDFZeMTGBEfxIj4oDb3HxEXRHZxNRn5lQcfUlm1prbeQmVdI1X1lj/9bENiA5g2PJZLxsbLcpOiyxLDfLl+chIvLdjFFRP6MKoDo+RE7+JUCbNS6k3gTCBfaz2kaVsI8BGQCOwBLtJalyhjhv5zwOlANXC11nqNI+J2pK7M1+uK3LJaTArC/LrvaW9UgBc/b8mjpKpenjI7keY2+I8zU7FYNel7Str1OqWMKrStaf49L9lVyPB23IAJ0ZVLYvNDSCeZEtrrPPfzTnJKa3jm4hGSLLfDoSHZrTfYmnoLVg13nTLAqaY0DY0L5PvbJx/1+1obhd9yy2vJKqpiU04Zv27L5/Hvt/Lawl38+7xhnJza+nuI6FnsUdTxpuP78nH6Ph75Zguf3zRRrj2iRU6VMANzgBeBdw7bdg/wi9b6P0qpe5q+vhs4DejX9DEOeKXpX5d18C3PBeYw55bVEuHvhZudnlS35KK0eGYv3s2Hq7K5+fi+3XZe0brmDgCTUgyMCWBwjO2GUA6PC2Tq4Cie/WknCSE+BHk794OSQdH+hLayTI3oPlL11LXM25zLawszuXRsAmOTQhwdjktob8Lc3Evr5mKJgFIKbw8zSWG+JIX5cvyACG49sR/r9pZy/xcbuf6ddB6ZNpgrJyQ6OlThwvw83fj7qQP426cb+Gp9DueOjHN0SMIJOVXCrLVeqJRKPGLzNOD4ps/fBn7DSJinAe9oo3truVIqSCkVrbU+0E3hOgWHzWEuryWym+YvNxsQ5c+kvmG8szSL6ycnH5zfKhyreT6cPe7FlFI8fu4Qpj63iFs/WGv7E9jYKamRzLoyzdFh9G5d6WE+OIdZupi7i9aaT1fv474vNjI8PoiHzkp1dEgu49CQ7Nb3a7Ra/7C/qxsRH8RnN03k1g/W8uBXm4kK8OKUwVGODkt0g4PLr9q4KZ8/Ko53l2fx7++3cXJqFH6e9kuPduRVMG9TLoWVxnJw542KxcfDqdIx0QJX+A1FHpYE5wLN429igb2H7bevadsfEmal1ExgJkBCQoJ9I+2ijg6vduRbX25ZLSnhft1+3muOSWTG2+l8v/EA00bEdvv5xZ/pgwmzfVpkqJ8n8+84loyCSrsc31Ye/GozFbWNjg5DNOlMa+wZ6YRrqG2wsHRXIW8s2s3SXUWMTw7htSvSHLrusqs5WPSrjXuHRotr9jC3xsvdzIvTR3LRa8u46+P1zL8zsFtrqoiexWRS/PPswZz78lJe+HUn9542yObnyC+v5eFvNvP9xlyUMnq2K2obmbUwk7OGR2M2mfB0M3HFhD4EyPx8p9OuhFkp9eBRvqWBWiAD+FFrXWOrwFo8mdZaKdWhrFJrPQuYBZCWluYS3QYdyTscNdeutKbBIfOITxgQQf9IP575aQenD42WXmYnYLVzwgzGXOYxvs49TDPQ281p1kPtzWzROyy/Rvsoqarnl235/Lwlj4U7C6iutxDm58mDZ6Zy1cTEHtMD2l2a51q2mTA3Dcnuaf+/Xu5mXrx0FKc8+zuPfLOFVy4f7eiQhAsbmRDMBaPjeHPxbi5OiyfZhp1Ci3YWcOsHa6lpsHD7Sf24ckIfQnw9WJZZxENfbealBbsO7hsX7C0dQk6ovT3MFwIJgC+wv2lbDFAFFADxQL5S6jitdftL5LZPXvNQa6VUNJDftD2n6bzN4pq29RqOnKNntWqHPK02mRT3nDaQa+ek886yLGZMSur2GMQfNY32s2vC7ArMJkVjg1TzdhadaY69vAnb1bvLs3joq01YNUQGeHLuyFimpEYyITlUepU7qaNzmHviA+aEUB/+cmI//jdvOysyixiXHOrokIQd2ftZ5t+nDmDeplzu/XwjH1w/vssPmbTWvL10D49+t5V+EX68fNmoPyTiE1PC+OnO4wDYX1rDxP/8Sk29pUvnFPbR3qvnU8AqIFFrnaC1TsCoWr0CeAQjed4BPG2HGL8Grmr6/Crgq8O2X6kM44GynjJ/uSN/no6qkm3R2i5zVtvjhAERnDAgnP/+uI3tuRWOCUIcZLHTnCJXY1JKepidwKHl+bpwDJtEIg63cV8p/l7ufH3rMSy/9yQeP3coJwyIkGS5C8yqd85hPtKMSUmE+Xny3C87HR2KcHER/l7848xUVuwu5vVFXev/a7BYuf/LTTz8zRZOGBDBpzdNbLXX2rvpWljTIAmzM2pvwvwQcKfWel/zhqbP/w48orUuAu4HJnQlGKXUh8AyYIBSap9SagbwH+BkpdROYErT1wDfA5kYw8FfB27uyrmdgSvda1us2mGl95VS/PeC4fh7uXPtnFXsLZa1mR2p+aFNT70Zay+zSWFtYz1U4ex6dxu2J4vVmLM3LC5IKpjbiKnpDq6t605PnMN8OC93Mzcel8zSXUWk7yl2dDjCnpofiNqxKV+YFsdpQ6J4ct52lmQUduoYJVX1XDF7BR+syOam41OYdcXoNguJeXtIwuzM2pswRwItlUT2BCKaPs8DfLoSjNb6Uq11tNbaXWsdp7WerbUu0lqfpLXup7WeorUubtpXa61v0VqnaK2Haq3Tu3JuZ9KRmwlH3Z5brdqhQ3DD/T2Zc80YKmobOPflJSzt5EVNdN3hy0r1ZtLD7BwOrtPZhfboqJE7PZlV64MJnrCN9g7J7qlzmA932bg+BPm4d7lXUAilFP85fxgp4X7c+O5qNuWUdej1O/MqmPbSEtZkl/LMxcO5e+rAdnUweboZF0gpHuqc2vv29TPwmlJqjFLK1PQxBmPt45+a9hkK7LZHkL2FKy1lYtWOf/MdEhvIpzdNJNDbnctmr+Df32+lVp7MdTurHZeVciUmpdocGinsryvJbi9/5mNXFqs+OIRY2MbBZaW0ZvHOQr5Yu6/Fj9+2G6VfeuIc5mbeHmYuHZvAT1vyZNRZD9Z8n2zvUSqB3u68dc0YArzduWTWchZsy2/zNVprPlqVzTkvLaG63sLcmeM7tKazUoqxSSF8tGovZdUNXQlf2EF7i35dB7yDMWe5OSMxAfOB65u+rgD+atPoeqn2XgYcee9hzGF2/M1P/0h/vr51Eo9/v5XXFmby2/YCnr54OINjAh0dWq9hPThEyvHtwZHMpraHRoru08ubo9OxaO3wh6w9TXMbf2dpFtvz2q7nEerX/StbdKcrJ/Rh1sJM3lm2h/vPkPW8RdfEBHnz2U0TuXbOKq6Zs4qL0uK4fUp/YoP+vHzZ+r2l/HfeNpZkGEvkPXPxiE4tc/bwWYM584VFPDl/O4+eM8QWP4awkXYlzFrrfGCqUmoAMKBp8zat9Y7D9llgh/h6lU51jDjo/twYku2Ycx/J19ONf507lJNTI/n7pxs456Ul3D114J/K8isFob4eNk/sMgsq+WT1PnblVxLs48GU1EhOGBCOWw9+mn84ay8Y7tceZpMMyXYGXfkN9O4WbF9WqyTMttbcY789r4Kzh8dwx5R+R32Q7eNhJiKgpZl1PUd0oDenDYli7qq93DGlP75tzBkVrscWRR07IirQi89vnsgzP+3gzSW7+WxNDmMTQxgeH4S/lxsFFXUszyxiW24FwT7u/PPswVwxvk+na/ykxgRw5YRE3l62hzOHRUvVdyfSoauJ1nq7UqoUKNBay+BDO2lvPqcceHtndcLeghMGRDD/jmP5+2cbeOy7rTz23dY/7XPHlH7cMaW/zc753vIsHv56M0pBYqgvKyqK+Sh9L0E+7kddeL5/pD9vXJVmsxgcTYZkG0xKin45ky5VyZZfo81ZHFz3oify9XTD38uNaSNieOTsIQ4rxOlMrjkmiW83HODzNfu4YkKio8MRPYCXu5l7Tx/ElRMT+WBFFr9szWf24kwaLBpvdzND4wL559mDOXdU7FHv+zrir6cO4PcdBdw2dy3f3TaZMD9PG/wUoqvalTArpdyBx4GbAG+gP5CplHoCyNJav2y/EHsPV+lg1lpj1c5Z5CnY14PXLh/NvM25FFbV/+F7//txGzklNTY935uLdzMoOoA3rx5DuL8nDRYrv27L5+cteQcLrRxu64Fyft6ah9a6xwxhliHZBulhdg5d+RX09jZsT874kNXVebmbSX9gCp5usjRXs1EJQQyPC+StpXu4bFzne/qEczpU1LH7zx0b5M3fTh3I304diMWqabBY8XQz2fx9w8/TjRenj+S8l5cyY84q3r9+fJsVtoX9tfc38BBwFnA58MFh21cCdwOSMNtQe3uOHXVv5+xVkU0mxWlDo/+0/fWFmTTYuCpTdb2FtMRgwv2NJ4DuZhOnDo7i1MFRLe7//C872ZZbYRRNc87/vg7T0sMMGMMjLdLD7Dy6UiXbhQowugqLDMm2C0mW/0gpxbWTkrh97jp+3prHKUd5LxaiK8wmhdlkv7+9wTGBvDR9FDe8t5oZc1Yx++oxkjQ7WHsnWV4K3Ki1/go4POPYhNHbLGygMz0jjlj+xHpw3d1uP3WXuJsVDRbb/n/VNloOLjbfHs03jI3WnjOjwSJzmAHjQY0MyXa8riS7vbsF25fFSUcliZ7njKHRJIT48NKCDFkiroc5NIe5519LpqRG8vRFw0nPKmH668spPmLUpOhe7X1cEQNkHeX1LvPIY/P+cgb948c/bY8I8OTH2489uGi4o7V/DrNjNCdIrjbUycPNTL0NepjLqhtoaEp4q+steHWg3bg1LwPSgxIrZx9x0F3Msg6zU5E5zM5Fin6J7uJmNnHT8Snc+/lGFmcUMrlfuKNDEqJTpo2IxdfDjVs+WMOFry7l3RnjiGmhSrewv/Ymu5uBY4E9R2y/CFhty4DsKcTXgysm9PnDtq0Hylm0s5DSmnq8PRzbCDvTM+KI+7pDRZ5c6+bHw6y6NCRba80TP27n1d93/WF7R4o8HOph7thvTmvNttwKp1xnOqfUWPPSxZqDzZnNioKKOk59ZiE+nmZ8PdwI9vUgKsCTyAAvogO9iQo0Pvfx+POl18fDjFcHRiuIozg4p77jL+3tbdieZB1m0Z3OGxXLcz/v5OmfdjCpb5jUJ+ghDq3D7OBAutGU1EjeuXYs172dzsWzlvHh9eOJC/ZxdFi9TnsT5n8C7yml4gEzcKFSaiAwHTjDXsHZWnSgF/edPugP2z5cmc2inYUu2aPgqAvGwSG4LnbFcjebOp0wa6353zwjWT53ZCyjEoKApvnSQ/48X/poDvYwd2Bo+NrsEv71/VZW7SnpUMzdydPN5HLtwdYuHB1HTb2FqrpGahosVNY1srekmvlltdQ1tt3u/D3dWHH/SS0m06J7ueL7gbOzaI3JxabxCNfl6WbmzlP68/dPN/D1+v1/WmZSCFcyLjmU964bx+WzV3Dp68slaXaA9q7D/I1S6iLgPow5zA8Ba4CztNY/2zE+u3P1W3xH3NgdHILrYsPr3M0mSqoaOlV85uXfdvHyb7u4dGwCj5/T+eU7zE0Tv9vbw/zvH7by2u+ZhPl58NBZqSSG+XbqvPYWE+jda9adPpqRCcGMTAj+03atNWU1DRwoqyW3vJbcslrqj0ig1+0t5Yu1ORRV1uMTIglzVxysotqJq3tvmBfnKFarxt29d18jRPe6YFQcby/dw39+2MYpqVFOM+1OdF53r8PsTIbHB/HejENJ8+c3HXOw4Kywv3bfmWmt5wHz7BiLQzR3ijlDh4Lu4FDC7h5ipLVm9uLdZBZWAa5XFXlCSihP/7SDs15YzJUT+jAhJZS4YJ8/Jc9aa6rqLVTXNx78nXy34QCj+wTzr3OHdOn/vaNzmD9atZfJ/cJ45fLRUiHRRSmlCPLxIMjHg0HRAS3uE7bhAF+szaG63vmG3PdGzvB+4Eh7i6vJLq6moKKOwso6ymoaqKhtpLymgfLaBsqbPq+obWzX6AmA0up6JvYNs3PkQhxiMikeOmswF722jKfmb+eBM1MdHZIQXTI8Poh3Z4zj0lnLue7tVcydOUEeBHWTXn8H3tyj4KqVFLtz+ZPdhVU89t1WvN3NhPh6MCDKv9vObQt/ObEvfUJ9eO6Xndzz+UbASPpDfD1wM5kwmxR1jVbKaupbrKZ96dj4Lj+k6GiV7LoGKwOj/CVZ7uF8PI03vKr6RgdH4vo6+uDxcL18VgFgXOdPePK3P2xTypgyEODtjr+XOwFebsQF+xDg7Yanm7nd/2+nDZElfkT3GpsUwhXj+zB7yW5OTo1kXHKoo0MSXXDwzqwXX6tHxAfx/KUjmfluOrfPXcurl492uRGfruiod+FKKSvtfNCutXbdxxvNPcxOlC+3ex1mO8dxpKW7igD4/vbJJDnp0ODWKKWYNiKWs4fHsDO/ktVZJewvraGwsh6L1YrFCh5uTb2B3u74eLodnJerFJwwIKLLMXSkh1lrTV2jRQpB9QIhPh4ArMkqYVQLw7pF+9niIaKrPkC1healS+4/fRAnDoogzM8Tf083uSETLuue0waycGcBd368nm//MolgXw9HhyREl5ycGskDZ6Ty6LdbeOX3XdxyQl9Hh9TjtdZtdRGHEuZI4BHgC2BZ07YJwDkY85mFg3Tnfd3vOwqICfQiMdS1Cw0opegf6U//yO7vIW/uYW7PetANFo1VGwW1RM82LC6QEwdG8L952zmufzj9HNA2expJ7zqn+WFeakwAKeF+Do5GiK7z9XTjuUtGctGry7ht7lrmXDNWljhzVU03vVJvAq49JpF1e0t5av52RiYEMTGle6a8aK3ZvL+cspoGRsQH4dtLRkAe9afUWn/a/LlS6mvgXq3164ft8qZSaiVG0vyy3SK0s+Y/OWfoUGju1Wj3sMBuvF5U1jXy+44Cpo9NkOUZuqC5quHdn23gygl9Wq2I3LyElPQw93xKKZ44fxinPruQOz5ax1UTEv/wfT8vN6YOjpJevnawxbXcCd4OHKZ5uogkFKInGREfxCPTBnPP5xv537zt3HPaQEeHJESXKKX4z3lD2XqgnNs+XMu3f5lMVKCXXc+5aGcBj3yzhZ35lYAxatK7hXvUAVH+fHrTRLvG0t3a+1jgRODOFrYvAJ61WTQO0Jz8dedcYFvqrqgXbMunvtEqc9C6aHSfYJ67ZAT//GYLt89d167XhPrJ8LHeINzfk3+dO4RbPljL3z/b8Kfvf3nLMYyID+r+wFyUzGHunOYeZjdJmEUPc8nYBDbmlPHq77uIDfbmivF9HB2SEF3i6+nGq5eP4uwXl3DrB2v4cOZ43O2wYkl9o5V//7CVt5bsoU+oD/89fxiRgV4szyyiruGPNXnWZJewOtt5l0HtrPYmzIXABcB/jth+AVBg04g6SCk1FXgOY33oN7TWR8bY+uub/nWOHmbj3/Z3MHffDc2Pm3IJ8/MkLTGk287ZU00bEcvUIVFkFVX/aXmhI7mZFf0jZHhubzF1SDTp94dS3XCoWvaGvaXc9P4aKmobHBiZ67DJpdwJ3g8cpXnJO+lhFj3RP88eTF55LQ9+tYlQXw9OHxrt6JBEBxxcNlAuTwf1jfDnP+cP47YP1/KfH7bxDxtXg8+vqOWm99awOquEa45J5J7TBuLpZvQqH9c//E/7P/fzTtbtLcVq1T1qVFx7E+YHgbeUUidwaA7zeGAKMMMegbWHUsoMvAScDOwDVimlvtZab2n/MYx/Xfb+qBsCr22w8Ou2fM4bFSs3UTbi6WZ2yBxq4fyCfT04vOxXUWUdALUN7ausLgydWodZ7sKwNNVXsEcvhRCO5mY28cKlo7hi9grumLuOQG93jpHlzoSLO3t4DGuySpi9eDejEoI5Y5htHgRlF1VzxZsryC+v48XpIzlzWEybr2l+62i0ajx6UM7QrndErfU7wESMnuazmz6KgGO01m/bL7w2jQUytNaZWut6YC4wrSMHOJgwO0EX86EnZ+2skt1N7XDhjgJqGixMleHYQnS75jnsdY2yRnN72OJa7qpTdGxB5jCLns7bw8zsq8aQHO7L9e+ks9aBw0cbLVbyy2vl+t5OHR2J2Zvcd/ogRiYE8fdP17OroLLLx9t6oJzzX11KWU0DH1w/rl3JMoDZZKSWVifIq2yp3Y+QtdYrtNaXaa1HNX1cprVeYc/g2iEW2HvY1/uatrXbwXWYbRdTt+qOG7v5W/II8HJjvKxfKES382oa+iQ9zB3UmTnMto/C5TTKHGbRCwT6uPPOtWMJ8/Pkmjmr2J5b0a3nzy2r5d7PNzDsn/MZ+69fGPnIT9z/xUbKamTqjegcDzcTL00fhae7mRvfXU1lXWOnj7VydzEXvbYMN5PikxsmMLIDy102v3c0tmP5VFfS48dcKaVmKqXSlVLpBQV/nm7tTCPwOj6H2f4aLVZ+2ZrHSYMiZYieEA7g52XMnPnvj9t4+OvN/LwlT+Yzt8IWPRA97MF4h1hkDrPoJSICvHhvxjg8zCaumL2CvcXV3XLehTsKOO25hXy2JoezhsXw6DlDOGtYDB+uzObi15YdnIYj/uzQajJyfWpJTJA3z18ykszCKm58d3WnRi78vCWPK2avINzfk09vmtjhZS6b5y1b2rF8qitx9cWzcoD4w76Oa9p2kNZ6FjALIC0t7ai/PVe9QbJ33Cv3FFNS3cApqZH2PZEQokUhvh68cWUac1ft5cOV2cxZugezSTE8LpBj+oYRH+LTanI4PjmU+BDXXju9u8g9GDRamnuY5QGp6PkSQn14d8Y4LnptGZe9sYJPb5xARID9lub5ZWseN763mpRwP16+bBTJh611ftbwGGa8vYpbP1jLe9eNk4dWolMm9QvjifOH8ddP1nPnR+t59pIR7e7w+njVXu79YiODYwJ46+oxhPp5dvj8zT3MFldNrI7C1RPmVUA/pVQSRqJ8CTC9Iwc49JTK8b/Y5uHV7b1p646bu/mb8/B0M3HcgD9XwhNCdI8pqZFMSY2ktsHCmuwSlmYUsWRXIS//tutgj+DRnDYkilcuH91NkTqPrvRAOP7dwLbWZJfwnx+2HVzbvTXFVfUASL4seosBUf7MuWYMl72xgivfXMlnN03E19P2t8cb95Vx0/trGBQdwHvXjSPAy/0P35/UL4xHzxnC3z/dwJyle5gxKcnmMbi6g7V+HBqF87tgdByl1fU89t1WquobefmyUfh4HL1Na6158dcMnvppB5P7hfHK5aPx6+TfgPngkOyeNY3MpRNmrXWjUupWYB7GslJvaq03d+5YNg2t29gzbK01P23JY3K/8Fb/0IQQ3cPL3czElDAmpoTxVwZQUdtAafXRh2df/046Ne1IkoShO5fq6y6/7yjgxndXE+TjzsCotofWhfp6MDEllOhA726ITgjnMDIhmFcuH83Vb63k/i828szFI2w67LesuoGbP1hNmK8Hc64Z+6dkudmFo+P4bsMBnv15B+eMiOlUD58QANdNTsbX0437v9jItBeX8OwlIxgcE/in/Uqq6rn7sw3M35LHuSNjeeL8YXi4df6JaXPC3NbDfFfj8lmQ1vp74PvOvt55+pcPm3vX3irZdr6525RTTk5pDXdM6WfX8wghOsffyx3/o9x4gZFg97Q3rbZ0ZQ7zoVUTbBaOQ32zfj93fryOvhH+vH3tGCL87TfUVAhXd1z/cO6c0p+nftpBWmIIl4/vY7Nj/+fHrRworeXjGycQ4utx1P2UUvzjzFSmPruQF37N4OGzB9sshp7g0H2yY+NwFZeOTSAu2Ju7Pl7PWS8sZtqIWM4eEUPfcD9Kquv5eWs+by3ZTW2DhQfOGMSMSUldflAkCfNRKKUeBOZrrZfbIJ5u5+o3SPZcDmve5lxMCk4aJPOXhXBFZpPqcUs7iPb5Yu0+7vx4PWP6hPD6VWkEeh/9wYoQwnDLCX1ZlVXCY99tYVLfMBLDfLt8zK0Hyvlo1V6uOSaJUe2oNtw3wo9zRsYyd1U2t5zQl3B/6WUWnTe5Xzjz7jiWl3/L4L3l2Xyx9lCpJ6XghAER3D11IAPaMQKpPcyqZybMtpil9DCwWCn1pg2O1e0OLSvl+F+s4yP4o/lbchmbFNLq01AhhPMyK9Xj3rTa0tFaEIc7NOLItf/PvlqXw10fr2dCcihvXztWkmUh2slkUvz3/GG4m038/bMNWLt4/dRa8/h3Wwnwdue2E9s/Wu/m41Ooa7Qye/HuLp2/pzk0h1m6mDsi2NeD+89IZe2DJ/PujLH894JhvHr5KJbdcxJvXj3GZskygJtZlpVqkdbaBPQD0rseTvdz5R5mpeyXZO8urGJHXiWnDo6y0xmEEPZmMvW8p7xtccVruS0t2JbP/320jjGJIcy+agzeHmZHhySES4kK9OIfZ6aycncx7y7P6tKxfttewOKMQm47sR+BPu1/cJUc7seZw2J4d9keylqpUyFER3i5m5ncL5yL0uKZOiSaqEDbT9MJa5p3//mafTY/tiPZpA6m1nq31vplWxyrux3sUXCGmyytnWZexrzNuQCcIgmzEC7LzWTqdQlzs05dS134AWqzL9bmEOLryZtXS7IsRGddODqOyf3C+O+P29hfWtOpYzRYrDz23RaSwnw7NR/6lhNSqKq3MGfpnk6dvyfSXSlSIbrFxJRQLhkTz0sLdvH8LzvtOnW0O7UrYVZKBSilzlNK/VUpdZdS6lyllO367x3oYA+zCw7Bs+f1Yv7mXIbEBhAbJJVShXBVJpPC4nqXti7pZT/un9Q0WAj397TLsjhC9BZKKf517lCsGv7x5aZO3fTPXZnNroIq7j1tYKeqDg+MCmDKoAjeWrqbqrrGDr9eCEdQSvHoOUM4b1QsT/+0g3s/30iDxfWXmGrzL1gpdRGwB/gU+C/wP+AzIEspdaFdo+sWTXOYneAuS9PxJNgeceeX17Imu5RTU6V3WQhXZlZg6WFrIbZXZ+a4Happ4bpqGyx4u8siykJ0VXyID3ed0p9ftuXz3cYDHXptWU0Dz/y8k/HJIZyc2vnCqTef0JfS6gY+WJHd6WP0RM4yGlO0zN1s4qkLh3PrCX2Zu2ovV85eSWl1vaPD6pJW31WVUkOB94AfgNGAN+ADjMFY+/i9pn1clkv/0dkp+Obh2KcOkYRZCFdmNpnoAQ92O6S5J8ilr+1dUFNvwctdhmILYQtXT0xkaGwgD3+9pUNziV9ekEFJdT0PnJHapWV6RiUEMzEllNcXZVLbYOn0cYTobkop/nrqAJ6+aDirs0o49+WlZBZUOjqsTmtrzNZtwA9a68uO2L4auFQp9VXTPtfbI7ju4Ez3VFq3fw1me/p2wwH6RfjRP7JHjLoXotcym+hyldfeRDlVUYs/q2u0sHl/OeuyS1mdXcKKzCLKa/84VLO+0cqUQREOilCInsXNbOI/5w/l7BeX8K/vt/LEBcPafE1WURVvLdnDeSPjGBIb2OUYbj2hL9PfWMGnq/fZdG1oVyRTmF3PeaPiSAjxYea7qznnpSW8cvlojukb5uiwOqythHkycEcr338FeNZWwThCc4LqpPdHrbLHBSO/opaVe4q5/aT2L38ghHBOZpPC4ooXty7oST9tfaOVBdvzWZ5ZxJrsUrbsL6OhaVJ6dKAXx/YLJyLgz1VOTx3c+SGgQog/GhwTyPWTk3n1911MGxnDxJTWb/Yf/XYL7mbF36cOsMn5J6SEMiI+iFd+28UFo+NkBIlwOWmJIXx1yzHMeHsVV765kn+fO5SLxsQ7OqwOaSthjgW2tvL9bUCc7cJxHGco+qXRnUqCtdY265n+YWMuWsMZQ6NtcjwhhOOYeuE6zF1xaB1m5/D20j08/v1WvNxNDI8L4tpJSYyMD2JEfLBdlgMRQrTsjin9+GHTAe77fCM/3nHsUZPWeZtz+XlrPvecNpDIFh5mdYZSir+eMoDLZ6/gzSW7ufn4vjY5ritqvld3htGYomPiQ3z47KaJ3Pz+Gv7+2QZKquu54bgUR4fVbm1VBvEFWqunX4Mxr9llOfkIvFbZ+nqhtebj9L0Mig6gnwzHFsLluZkU1fWNVPaiCqsHh+z1gBuqilpjzuT6h07hoxsmcO9pg+y2dqYQ4ui83M38+9yh7Cmq5p7PNrQ41eVAWQ13f7aB1OgArj0myabnn9QvjJNTI3nx1wzyy2ttemwhuou/lzuzrxrDmcOi+fcP2/j3D1tdZtmp9qw7MUwpVXyU77neIPQjHFpWyvE622aMuc9dP//6fWVs3l/Oo+cM6frBhBAOlxzux5fr9jP8n/NJDvMlyMcdgIpaI4muqbfQYLHiZjbhZlJ4e5gZmxhy1PV7Pd1M3Hx8X4J9Pbrzx+g2zjZFx6rBpMDTTYZgCuFoE/uG8bdTB/C/edvxcjfzyLQhB5eLyi+v5crZK2lotPLC9JGdWkaqLfefPohTnlnIw99s5qXpo3rEQ0HR+3i4mXjukpEEervz2u+ZlFY18Pi5Q3AzO/fqDu1JmOfR+nRZJ7m16JyDCbOT3CF15PrXmWVTWvPO0j34eJg5Z0SMTY8rhHCM207qx+R+Yfy8NY8deZVU1TWiNSSE+ODn6YaPpxk3k4lGq5VGi6agoo5ftuW3eD1stGoqahsZHh/EmcOc+RrRNGTPwVHYgkVrzKae8JMI0TPcfHwK1fWNvLRgF+lZJZw7MpbKukbmrsymtsHKm1ePISXczy7nTgzz5fYp/fjfvO18tW4/54yMtct5nJkU/eoZzCbFY+cMIdTXg+d/zaCspoFnLxnh1PPz20qYbTumxAn1hHU3bRH7nsIqvlq/n6snJuLv5W6DIwohnMHIhGBGJgR3+Ti7C6s44cnfaOjB61QdmqLjHO8IVqvGJL1IQjgNpRR/O3UgI+ODefqnHfxv3naUgmP7hfOPMwfRN8K+09luPC6FX7fl84+vNjEiPojEMF+7nk8Ie1FKcecpAwjy8eCRb7dw7ZxVzLoyDT/P9vTldr9Wo9JaZ3VXIA5zsIfZsWGAkfh2pNfYlvdRz/68AzeT4objkm13UCFEj+HW1NPZXKXZWR2aw+zYOGzBYpUeZiGc0ZTUSKakRlJR24C72dRtPWNmk+KZi0Zw9kuLmfH2Kj6/+RgCvXtPJ0fzu09PuL4Lw7WTkgj2deevn2xg+uvLeevqMYT6eTo6rD9pdcC4UiqhPR/dFaw9HPqbc+6bwNZ0tTdkeWYRX67bz4xJSUT4SzEZIcSfuZmNq2WjsyfMXXitM9W0AGNItvQwC+G8/L3cu30YaUKoD69ePprs4mqufye9VxV1FD3TuSPjmHXFaLbnVnD+K0tZm13i6JD+pK0Z1nuA3a18NH/fZTlTkRdtdDG3my1uoypqG7jnsw3EBXvzlxNl7WUhRMvcTMbbhcXqGkOybV3jwRGMIdmOjkII4WzGJ4fy9EUjWJ1VwuVvrOg1lbMPzWGWC2NPc9KgSN6/bhz1jVbOf2Upj3+3hbKaBkeHdVBbCfMYYGzTxzigFjjvsG3N33dZPeFPrrO5vtWqueezjewtqeHpi0YctTKuEEK4m3v+kOyDNS2c5Ee0amRIthCiRWcNj+GVy0axLbec055bxNfr97e43JUQriItMYQf/+9YLkqL5/VFu5n8xK+88MtOCirqHB1a6wmz1nr1YR/pgBXYeMT21V0NQil1oVJqs1LKqpRKO+J79yqlMpRS25VSpx62fWrTtgyl1D2dP7fxrzNcYjS6Qwl8V0bqWa2a+7/cxHcbD3D31AGMTQrp/MGEED1e85IPjS7Sw9wTSJVsIURrThkcxbd/mURUoBe3fbiWU59dyEsLMli5u5jCyjrqG3vW9Vo3r4Igl8UeK8DLnf+cP4zvb5tMWmIIT/20g4n/+YUb313NZ6v3UVTpmOTZWUqRbcLouX7t8I1KqVTgEmAwEAP8rJTq3/Ttl4CTgX3AKqXU11rrLZ0NwFl6FDqjo7FrrXnw6018uDKbm49P4frJUuhLCNG65qJfjU7eg6G7sqyUEz1ABamSLYRoW98If76+dRLfrN/P28v28L952//wfQ+zCU83E57uJjzdzHi4GV8nhPjQP9KfY/qGMSYx2OnXwRW9S2pMAG9ePYZdBZV8uCKbr9bv58fNuQAkhvowOCaQPqE+hPl5EuDtjlVrGixWymoaKKmqt3k8TpEwa623Ai0twj4NmKu1rgN2K6UyODQEPENrndn0urlN+3Y4YT40BM8JbpF0B9dh7sSNlNaah77ezHvLs7nhuGT+duqATh1HCNG7uDf3MDv5kOyeRKpkCyHaw2xSnDMylnNGxpJfXsvmA+VkFVZRWddIZZ2FukYL9Y1W6po+auobySys4pdt+by4IIMwP0+un5zE1cck4unmvNPznOFWXXSvlHA/HjgzlftOH8SWA+X8vqOATTllbMwp48fNuVhaeIjvbYdCfJ1JmLuzucYCyw/7el/TNoC9R2wf19IBlFIzgZkACQl/LujtTEOyO0u3M3qtNY98u4V3lmVx/eQk7pk6UJJlIUS7mJRxvfx5ax59Qn0YGR9MbLC30yV0XZrDfHCZQed4R5Aq2UKIjooI8CIiwAsGtL1vVV0jC3cU8MHKbP79wzY+X5PDMxePIDUmwP6BCtEBJpNiSGwgQ2IDD26zWjVlNQ1U1DZiNivcTIpAb6NyvXrMtudvNWFWSn19xCYv4HWlVPXhG7XWZ7d1IqXUz0BUC9+6X2v9VVuv7yyt9SxgFkBaWtqf7oLUwf3sFUH7dXQd5o56av4O3lqyh2uPSeK+0wdJsiyEaDelFHdO6c+cpXu4fe46wEiiw/098fV0w8NswsPN1OnriknBHVP6c1z/cBtG7dq0BpOMkhRC2ImvpxunDY3mtKHRLNiWz92fbeCcl5YwODaAUF9PXpw+stuXzWoPuX0VYCTRwb4eBPt62P1cbfUwFx3x9XudPZHWekonXpYDxB/2dVzTNlrZ3jEHe5idIGPupPYk+3OW7ObFBRlcMiaef5wpybIQouP+clI/bjo+hW25FWzYV8aBshryymupqjeG+zVYrJ1++Lg4o5Dftxd0OWE+dP6OX+Oc7aposWrMcq0WQnSDEwZG8MPtk/nvj9vZmFPGz1vz2F1YxaBo6W0WotWEWWt9TXcFchRfAx8opZ7GKPrVD1iJcV/TTymVhJEoXwJM78wJlBNVedFad3AOc/v2S99TzGPfbWXKoEgeP3eoJMtCiE5zM5v+NCzKFkY8Mt9pKnDba8RRo8WYP1jbYDk4l7Cu0cLRfuyS6npMTjbkXQjRc4X6efLEBcNYnlnEJbOWU+igisRH0zxdRtZhFt3NKYp+KaXOBV4AwoHvlFLrtNanaq03K6U+xijm1QjcorW2NL3mVmAeYAbe1Fpv7ty5jX+dIF+2i4raBv7y4Vpig7156qLhTjffUAghANxMJpus8dyVZUds8TCxpt7CN+v3szijkF0FlRRV1lNe20Bdo7XF4iRtGWrjBxNCCNGWvhF++HqYeXLedsYkhjjlsGwhupNTJMxa6y+AL47yvceBx1vY/j3wfVfP7VRzmHXHhgS25wnb0z/tILe8ls9vmkigt3vngxNCCDtyNysaLU7Sw9yFR6hnv7iYnfmVRAZ4khodwOCYAAK8jCIkhy/t4nXYEi+tFfYaGOXf6ViEEKIzwvw8eebiEcx8dzX3fLaBZy4e4RSjE7tS1FGIrnCKhNmRmi8APXEO8+b9Zby9dA+XjUtgZEJw9wYlhBAd4GZWNlnj+eANVSdea4t7sKziai4bl8Bj5wxxihtMIYTojFMGR/HXU/rz5PwdDIwO4MbjUhwdkhAO0+vrbzrT/YzGNkMCmz09fwcB3u787ZSBNjumEELYg7vJRIMNe5g7cy09tKxU58/baLES7OMhybIQwuXdckJfzhwWzRM/buPXbXmODudg15ZcXUV3kx7mpn+PdoNUUdvAroIqMvIrKa6qo7reQpC3O/0j/UlLDMHDzXHPHA7Nv/5z8Bv2lfLLtnz+ekp/An1kKLYQwrm5mRWNNpjD7EhWq8aqkVoRQogeQSnF/y4Yzp6iKm77cB1f3DyRfpEyTUT0Pr0+YW6mgZKqetbvK2X93jI27Ctl8/5ycstrj/qaYB93Zh6bwnWTk3A3dz1x7ugc5tY8/8tOAr3duWpioo2OKIQQ9uNmMtmkSnbXhmQ3T9HpHEvTyd3NkjALIXoGbw8zs65I4+wXl3DdO+l8dcsxBPnYf93blhyawyzXWNG9en3C3Pw3d89nGzhQVntwW0q4HxNSQukb4XfwIzLAC293MyXV9azNLuWDFVk88eM25m/J5ZXLRhMV6NW9sTf9e2Tv+KacMn7ems9dJ/fH30t6l4UQzs/drMgsqOKrdTmkRgeQEOqDp5trVWZt7iE3m3r9bCchRA8SE+TNa1eM5tJZy7nlgzXMuWasTTqKhHAVvT5hjg3yIdTXg7hgby4bl8CoPsEMjQ1sNdEM8/Pk5NRITk6N5NsN+7n70w1c+NpSPrhuPPEhPp2ORWObLubnftlJgJcbVx2T2PWDCSFEN5jYN4zZi3dz+9x1gPHgMibQm6QwXxLDfEgM9SUu2LvNZDSrqOrg6zuqq3OYm3vI3WRIthCihxndJ5h/nTeUv36ynse+3cLDZw9ud0/vvpJqvlq3n437ymi0avpH+nH2iBgGRgXYOWohbKPXJ8xRgV6s/sfJnX79mcNiiA/24co3V3LRa8uYO3M8fUJ9bRjh0bV0ndqUU8ZPW/L4vyn9CZDeZSGEi7h76kDuPLk/O/Iq2JlXye7CKvYUVbGnsIqv1+2nvLaxQ8cLcMAyes3rLLvJkGwhRA90weg4tueW8/qi3Xh7uHH31AGtJs3FVfU8OX87c1dmo4GkUF883Ews2J7Py7/t4rxRsfzz7MHtHg3ZXLNHrrCiu/X6hNkWhscH8eH145n+xnIunbWcuTMnkBDa8Z7mzs5hPrwz5Nmfjd7lq6V3WQjhYtzNJgbHBDI4JvAP27XWlFQ3sL+0pl3H8fYwkxzW+QeXnV1msKFpSLb0MAsheqp7TxtEdb2FV3/fRXltAw+dlfqn6TMWq+b9FVk8OW87VfUWrpyQyHWTk4gLNu6NS6vreW1hJq/9vos1WSXMuWYsiV24Zgthb5Iw20hqTADvXzeO6a+v4NLXlzN3ZteGZ7eHOiK93rCvlJ+35nHXyf0JdEDvihBC2INSihBfD0J8HVNopr0O9TDL3D4hRM9kMikeO2cI/l7uvPr7LtZll/LXU/tzTN8w6hqt/La9gFd/28WWA+Uc0zeUh88a/KfK2kE+Htw9dSAnDoxg5jvpXPDqUuZcM5YhsYFHOavhUNEve/10QrRMEmYbGhwT2JQ0Lz+YNDc/TWuvzlT+001XkGd/3kmQj7v0LgshRCc0X36LKuvZeqCc+kYr9RYr9Y1W6hotFFXWU1BZh9Xacg9087BxWVZKCNGTKaW457SBjEoI4uGvN3PtnPQ/fD8x1IeXpo/i9KFRrd7XjkkM4dObJnLl7JVcMms5s64czcSUMHuHL0SHScJsY0NiA3nvunFc9sYKznlpCf85bxgnDYqwSwn8ww+5NKOQX7fl87dTB0hlbCGE6ASTUriZFLMX72b24t2dOoabSdHHzqOLhBDCGZwyOIrjB0Tw2/Z8thwox91sYmR8EOOTQzG188FhSrgfn900kSvfXMHVb67i+UtHMHVIdIv7Nj+qlGWlRHeThNkOhsUF8dlNE7ntw7Vc9046o/sEM21EDENjA+kT6kuwj3uLf+xa604NM2m0aB76ejPxId7MmJRkg59ACCF6H3ezifeuG0deeS0eZhMebsaHp5sZDzcTIT4ehPt7tlrUSyFDsoUQvYeHm4lTBkdxyuCoTh8jKtCLj2+YwLVzVnHz+2v4+9SBzJyc3O6kWwh7k4TZTvpH+vPVrcfw8aq9vLlkDw9+tfng93w9zMQF+zA4NoCzh8dwXP/wLj0te2vpHnbmVzLritF4ubvWuqVCCOFMxieHOjoEIYTodYJ8PHj/uvHc9ck6/vPDNhZsy+cfZ6b+cV5zZ9f8E6KLJGG2I083M1dMSOSKCYlkFVWxI6+SrKIqckpr2Ftcwy9b8/l8TQ6jEoJ47JyhQOeqZL+8IIMTB0ZwcmqkbX8AIYQQQgghuoG3h5mXpo/ik/R9PP79Vs58YTEnDYzgkrEJnDAg3NHhiV5MEuZu0ifU90/rMzdYrHy2eh//m7ed059fBEBoB6rANvdK+3u58Z/zhsqcDiGEEEII4bKUUlw0Jp6pQ6N4c/Fu3l+RzS/vpBMZ4Emor6dUyBYOIQmzA7mbTVwyNoGpQ6L4cOVequoaGRwT0O7Xe7kb8+T+fd4wIgK87BWmEEIIIYQQ3SbAy507pvTnlhP6smBbPh+uzOa3HQV4y9RD4QBK96L5AEqpCmC7o+MQwgbCgEJHByGEDUhbFj2BtGPRU0hbFj3BAK21f9u7tU9v62HerrVOc3QQQnSVUipd2rLoCaQti55A2rHoKaQti55AKZXe9l7tJ2tfCCGEEEIIIYQQLZCEWQghhBBCCCGEaEFvS5hnOToAIWxE2rLoKaQti55A2rHoKaQti57Apu24VxX9EkIIIYQQQggh2qu39TALIYQQQgghhBDt4vIJs1LqTaVUvlJqUwvfu0sppZVSYU1fK6XU80qpDKXUBqXUqMP2vUoptbPp46ru/BmEgA635YFKqWVKqTql1F+P2HeqUmp7Uzu/p7viFwI63I4va7oWb1RKLVVKDT9sX2nHwqE62JanNbXldUqpdKXUpMP2lfsL4TAdaceHbR+jlGpUSl1w2DZpx8KhOnhNPl4pVdZ0TV6nlHrwsH07fH/h8gkzMAeYeuRGpVQ8cAqQfdjm04B+TR8zgVea9g0BHgLGAWOBh5RSwXaNWog/m0P723IxcBvw5BH7moGXMNp6KnCpUirVTvEK0ZI5tL8d7waO01oPBR6lac6RtGPhJObQ/rb8CzBcaz0CuBZ4o2lfub8QjjaH9rfj5uvvE8D8w7ZJOxbOYA4daMvAIq31iKaPR5r27dT9hcsnzFrrhRjJw5GeAf4OHD5JexrwjjYsB4KUUtHAqcBPWutirXUJ8BMt/EKEsKeOtGWtdb7WehXQcMS+Y4EMrXWm1roemIvR7oXoFh1sx0ubrrkAy4G4ps+lHQuH62BbrtSHisL4HvY9ub8QDtXB+2SAvwCfAfmHbZN2LByuE225JZ26v3D5hLklSqlpQI7Wev0R34oF9h729b6mbUfbLoRDtdKWj0basnA67WzHM4Afmj6XdiycUmttWSl1rlJqG/AdRi8zSFsWTuho7VgpFQucS9MIzMNIOxZOqY37iwlKqfVKqR+UUoObtnWqLbt1PVTnopTyAe7D6JoXwmVJWxY9QXvasVLqBIyEedLR9hHC0dpqy1rrL4AvlFLHYkwxmNKN4QnRLm2042eBu7XWVqVUt8YlREe10ZbXAH201pVKqdOBLzGm5HZKT+xhTgGSgPVKqT0YQ/zWKKWigBwg/rB945q2HW27EI7UWls+GmnLwtm02o6VUsMw5ntO01oXNb1G2rFwRu26JjcNG0xuKj4jbVk4m9bacRowt2n7BcDLSqlzkHYsnNNR27LWulxrXQmgtf4ecO/KNbnH9TBrrTcCEc1fN/0HpmmtC5VSXwO3KqXmYhQuKNNaH1BKzQP+dVgBg1OAe7s5dCH+oLW23MrLVgH9lFJJGBeAS4Dp9oxTiNa0cU1OAD4HrtBa7zjsZdKOhdNpoy33BXZprbUyVuDwBIoAub8QTqWNe4ukw7bPAb7VWn/ZVPRL2rFwKm1ck6OAvKZr8liMTuIioJRO3F+4fA+zUupDYBkwQCm1Tyk1o5XdvwcygQzgdeBmAK11McbwqVVNH480bROi23SkLSulopRS+4A7gQea9g/QWjcCt2LcpG0FPtZab+6O+IWADl+THwRCMXox1iml0gGkHQtn0MG2fD6wSSm1DqMC68VNBUbl/kI4VAfbcYukHQtn0MG2fAHGNXk98DxwSdM1uVP3F+pQUUchhBBCCCGEEEI0c/keZiGEEEIIIYQQwh4kYRZCCCGEEEIIIVogCbMQQgghhBBCCNGCHlcluzVhYWE6MTHR0WEIIYQQQgghhLCD1atXF2qtw211vF6VMCcmJpKenu7oMIQQQgghhBBC2IFSKsuWx5Mh2UIIIYQQQgghRAskYRZCCCGEEEIIIVogCbMQQgghhBBCCNECSZiFEEIIIYQQQogWSMIshBBCCCGEEEK0QBJmIYQQQgghhBCiBZIwCyGEEEIIIYQQLZCEWQghhBBCCCGEaIEkzEIIIYQQQgghRAskYRZCCCGEEEIIIVogCbMQQgghhBBCCNECSZiFEEIIIYQQQogWSMIshBBCCCGEEEK0QBJmIYQQQgghhBCiBU6bMCulpiqltiulMpRS97Sy3/lKKa2USuvO+IQQQgghhBBC9GxOmTArpczAS8BpQCpwqVIqtYX9/IHbgRXdG6EQQgghhBBCiJ7OKRNmYCyQobXO1FrXA3OBaS3s9yjwBFDbncEJIYQQQgghhOj5nDVhjgX2Hvb1vqZtBymlRgHxWuvvWjuQUmqmUipdKZVeUFBg+0iFEEIIIYQQQvRIzpowt0opZQKeBu5qa1+t9SytdZrWOi08PNz+wQkhhBBCCCGE6BGcNWHOAeIP+zquaVszf2AI8JtSag8wHvhaCn8JIYQQQgghhLAVZ02YVwH9lFJJSikP4BLg6+Zvaq3LtNZhWutErXUisBw4W2ud7phwhRBCCCGEEEL0NE6ZMGutG4FbgXnAVuBjrfVmpdQjSqmzHRudEEIIIYQQQojewM3RARyN1vp74Psjtj14lH2P746YhBBCCCGEEEL0Hk7ZwyyEEEIIIYQQQjiaJMxCCCGEEEIIIUQLJGEWQgghhBBCCCFaIAmzEEIIIYQQQgjRAkmYhRBCCCGEEEKIFkjCLIQQQgghhBBCtEASZiGEEEIIIUTPoDVkZ8PSpbB+PdTUODoi0R7V1fCPf0BCAgQGwrRpsG2bo6MCJGEWQgghhBBCuLqGBnjlFRg4EPr0gWOOgREjjORr6lT45BNjH+F8iorg2GPhsceM39n06bB4MYwaBb/84ujoJGEWQgghhBBCuLCMDBg7Fm6+GUJD4fnn4fvv4eOP4fbbYcsWuOgiSEqC//wHiosdHbFoVl8P558PGzfC118bH6+8Aps2Qd++TtHTrLTWDg2gO6Wlpen09HRHhyGEEEIIIYSwhSVL4PTTwWyG2bPhnHNAqT/uY7HADz/As88aPZbe3nDVVXDbbTBoUOvHLyqCzEzYswfq6sDLC/r3h8GDjXM6g7IyeO01+OILI1Z/fxg/Hq6/3ui5PfL/w5n885/w8MPw7rtw+eV//N7+/TB8OMTHw8qV4ObWrkMqpVZrrdNsFaIkzEIIIYQQQgjXs3ChkSzHxMD8+ZCY2PZrNm6E556D994zEuATToBJk4zXenkZCXJ2trHf+vWQm9vyccLC4LLL4G9/g9hYW/5UHfPTT3DllUac48bB0KFQWgoLFhg/y9SpMGuWkXQ6m23bjIT4vPPgww9b3ueTT4zRAS+/DDfd1K7DSsLcBZIwCyGEEEII0QNs22b0okZHw6+/Gv92REEBvPqqMWx7yxawWg99z8MDUlONZG7YMGNocGIi+PhAZSVs3mwMHf7sMyPJfvJJuOGG7u/JffdduOYao5f8zTdhzJhD36upMX6+Bx4wetS//homTuze+NoybRr8/jts3w6RkS3vozWcdBJs2GD08vv5tXlYSZi7QBJmIYQQQgghXFxRkdGbWlFhDNXt06drx6uuNhLomhpjDnRISPuGW2dmwo03Gr28119v9IK2c9hwl82bB2ecAccdZwzFDghoeb8dO4z99u41kuZTTume+NqyZg2MHg2PPmok9a1ZvhwmTICnnoI772zz0JIwd4EkzEIIIYQQwm5qa2HdOmNIbEyMc81z7Snq642kb/lyY9jxhAmOjUdrePBBo8LzZZfBO++Ayc51lXfuNJLN5GRYtMiYs9yawkKYMsV43U8/OUdP87RpRuy7dxuVzNty4onGqILdu8HTs9VdbZ0wS5VsIYQQQgghuiInx+hpDAkxErjTTjOG88bFGT1oVVWOjrBn0NqohP3770aBL0cny2AMw370UXj8cXj/fbj7bvuer7HRKFjm5gbffNN2sgzGfOt584y51meeaSSdjrR9u9Hbffvt7UuWAe65Bw4cMIbBdzNJmIUQQgghhOisuXONnuQ5c4z1Y7/80qjc/O67Ri/ggw/CkCGwapWjI3V9Tz9tJMoPPGD05jqTe+81kvknn4TPP7ffeZ58EpYtg5de6lghr8hIo1K41kaRrepq+8XYlpdfBnd34yFTe02ZAikpRgGzbiZDsoUQQgghhOgorY01fe+7D445Bt5+27ihP9LChXDFFUbv2DvvwCWX2CeeujqjF3HtWuPrkSONocteXvY5X3f7+mtjyajzz4ePPrL/sOfOqK83Km7v2GEMzW9P1e6O2LfPWNLqtNPg0087V2Ts+++NXuYbbjDWO+5ulZVGT/cZZ8AHH3TstU88YfQ0b9sGAwYcdTcZki1c29at8MIL8K9/GQUKHPl0SwghhBCis/77XyNZnj7dqNLcUrIMxjq4a9caw4enT4c33rBtHFarkfj06WPMC334YeNj2jSjuvPHH9v2fI6wbp3xfzd6tPFgwhmTZTCqa3/0kfEw5brrjH9t6R//MNaUfuqpzlfkPv10uOsuo4L2vHm2ja893n8fysvhlls6/tqrrzaGor/+us3Dao2TtjbR4xQVwaWXGiX6b7sN7r/fGA4SF2csIn94KX8hhBBCCGf29ttGT9ellxpDrz08Wt8/JMQYDnvqqUY1ZVv17OXnG8WQbr7ZWFrohx+MzojaWqMnMToaLr7YWCvYVUeV7txpJHnBwUYvs4+PoyNqXVKS0RP6yy/GiAJbWbfOaHe33971nutHHzXay7XXQkmJLaJrvzlzjCkKnSk8FhlptIUPPzQeHHQTSZiF/e3fbzxd/fxz48lYdrZxMf/lFxg7Fv7v/4xF1YuLHR2pEEIIIUTrVq+GmTONtWHnzGl/b6ePjzG/+cwzjQS3q0nzhg3GfdSKFca83l9/Ne6nvL2NKsKnnQZLlxo9eU8+adxvuVrSvHOn8UCgoQF+/LHjay07ysyZxjD9O+80Hmp0ldbw178aDw3uu6/rx/PyMpL53Fxjjn13ycw0qptffnnne8gvvdTILRYvtm1srZCEWdhXba0xJCg7G+bPh0ceMQoUeHsbF8AffjAm7//+O5xwAuTlOTpiIYQQQoiWlZTAhRdCRIRR7KutnuUjeXoac0+7mjTPn28kZA0NxtI8117bcgLi7m5MhbvjDnjuOXjxxc6dzxF++w3GjzfWRv75Z6OwmqswmYz724oK21TN/uEHo6PpoYcgKKjrxwNISzOKbr38svHwpTs0z1nuyjz+s84yHj59+KFtYmoHKfol7Kv5zeDLL43E+Wh++sko5JCQYDwxCg3trgiFEEIIIdqmNVx0kXFPs3Bh15Y0qquDCy6Ab781Epabbmr/a995B2bMMKa5ff+9UUCpLVYrnHuucb5ffoHjj+906O1SWGjM5V282Oglrq837u1SUmDcOGM47qBBLffO5+UZtW5eeMEocPXdd0efH+7s7r3XKAy3cCFMnty5YzQ2GkuU1dfD5s0df0jTmuJi4/84NdXovOpsr297aG2cJzzc+P/oiunTjYdGBw4YD4WOYOuiX5IwC/tZssSoFHjnnUZxgrb8/rsxt2fUKONJorPPURFCCCFE7/Hxx8Z84H//25i/3FUdTZq1NhLJBx4whoN/9ln717AFo7dz9GhjWtyGDca8alsrKjLmx77yipHgJSQYSZKXlzE0eft2Yx+AgAAjnmHDjGS6utqYp7tggfHaG2805gK3Z51hZ1VVZfz8AQGwZk2LyV2bZs0yKlp/9plR/8fWXn/dGEL+6adGBXJ7WbvWuMd/5ZWOLSfVkq+/NjrivvvOmNN8BEmYu0AS5m7U0GD8UZSXw5Yt4Ovbvtd9/rnx5nH22caFwWy2b5xCCCGcW0WFcfMcEmLf3g8hWlNQYCQ+SUnGvGA3N9sc9/Ck+emnjaHTLbXz4mK45hojUbjsMnjzzc71NK5ebfSMn302fPKJbf+m5s+HK680/q9mzIC//AWGDv3jPlpDRobRqbJqlfGxebORLLu5GVW9p041Hh7072+72Bzpq6+MUZRPPWV0InVERQX062d8LFxon2ugxWL8nrSGTZvsd+/9t78ZhX5zc7s+krSuzpgWceGFLVadt3XCjNa613yMHj1ai27y/vtag9afftrx177wgvHa22+3eVhCCCFcxHffaT1hgvF+AFr376/13LmOjkr0VhddpLWHh9YbN9r+2LW1Wk+bZrTzk07SeuFCrS0W43tFRVo/84zWYWFau7lp/dxzWlutXTvfE08Y53r99a5GfsiTTxrHHDxY63XrOv76+vpDP3NPY7VqfcYZWvv5ab13b8dee//9xv/r8uX2ia3Zp58a55kzxz7Ht1i0jo3V+swzbXfMSy7ROjxc68bGP30LSNc2zCGl6JewPa2Np2gDBhjzZTrq1luNSo7PPWfMXxFCCNF7NA/FPOMMYy7jP/8J//ufMVLpkkuMZQl70eg44QQ++8wYjv3QQ8ZyOLbm6WmMsHvpJWNI8rHHGtWQ4+IgLMy4Jxo8GNLTjaU5u9rL+Ne/GkO6b7/dGCLdFVobf5N//avR27dypTHftqPc3Z13beWuUgqef96Yi9yRHuasLKO6+WWXGfO+7em884zh8Q8/bFyDbW3hQsjJMX4WWzn3XGM0w5IltjvmUciQbGF7v/1mVLyeNctYa7AzLBbjwvvVV0ZxjbPOsmWEQgghnFFdnVFU6euvjeF7jz9+aM6fxWIM03z99a69vwjREcXFRnGquDhjOZzOzEHtiKoq+OIL41zV1cYQ8DPPhJEjbXuenBxj7nBiIixb1rnh3VarMez65ZcPrS0tU+mO7rHHjOVV582DU05pe/9LLjGuhdu3GyvM2NuPPxpLkb3+Olx3nW2PPXOmUSE7L6/90zTbUlFhPFC65RZjOsNhZA5zF0jC3E3OPNN4wpj9/+3dd3iUVdoG8PuIlKChinRZmiJFqoIVC7LYQFQQWBFWFHdBZUERcP2woDQFRDfouoq0RcBFMQiIWFhZ0F1C76EaEsGEJMYAIfV8f9wz1wwhmWSSaZncv+uaK8lkMnmTnLzzPuc8z3Pi2OShpM6eZRfHPXs4M9W5s88OUUREQoy13F9z6VJufTNy5IWPyc3lyvP69ay1a9Ei4Icp5cywYcD8+Vzd7dAh2EfjWytWcJVu7Fhg+nTvvjY7Gxg6lEHQ88+zE7R6DHiWmemqFd61y/M18tq1rOV+6SWu+gaCtdxqKj0d2LfPd5MfmZlAvXo8dy9a5JvndLrnHvZKOnLkvPHn64A5THMfJGj27WPHupEjSxcsA+ySvXIl28/fey8DcBERCU9TpzJYnjq14GAZ4AWcs9nRM88oNVv867vvON6efTb8gmWAjaiefJIlD2vXFv/rzp5lCu/ixewYPm2aguXiqFyZk4GHDvF3Vpi0NK7wXn21b7qxF5cxwAsvcBuw5ct997xffgn8+qtv07Gd+vYFjh3z+z7SIRswG2N6GWMOGGMOGWMuGC3GmDHGmL3GmJ3GmG+MMU2CcZySz6xZDJRHjPDN89Wtyz0GMzI4i5SW5pvnFRGR0LF+PesgBw7kapUnDRpwxWXNGm5HKOHHWq7qdurETri9ejF4DaSMDKaRNmvGVb5wNXMmVz3792fWRlFSUphOvGoVU7ADGdCFg549eZ579VWew/LLzWVgefIk8NFHpV988lbfvuxBNGWK7yYkFy9m6nSPHr55Pnf33cdA/7PPfP/cbkIyYDbGVAAQBeAuAK0BDDTGtM73sG0AulhrrwHwLwBe5pKIzyUmAgsWcEuBOnV897ytW3Oma/9+1jVnZ/vuuUVEJLjS07ldTvPmrJ0rzkrVn//MFL9Jk/x/fBJYeXkMVIcO5cd9+rA06/bbWf8ZqKyCSZO4Evj3vzPjLVxVrcrg99JLWb8aG1v4Y2NjgZtv5lZQy5aVfi/d8ur991k/3r//+UFzZibw2GP8e8ye7f9GXwW56CJg3Dg2n/vyy9I/32+/sQ67f3//1P/XrQvceCPLC/woJANmANcBOGStPWKtzQKwBEAf9wdYa7+z1p51fPgjgEYBPkbJb84c/rOPHu37577jDl5IrVvH1Wul4YmIhIdx49gNdt684jeDiYhg3eW337I5koSPqVO5r+oLL7Bu+B//YKA2bBibJhWVgeALO3cyTXnIEP+sioWaxo2ZzZeZCdxwA8vh3GVns4N3p05s2rR2LfeOlpK59FL+vlu0YPZk375cqb/mGi48vfIKJwWD5Q9/4JiYMqX0z/Xpp8C5c8Ajj5T+uQrTty+wYwfrmP0kVAPmhgCOu30c77ivMMMAFJDXABhjhhtjYowxMUlJST48RDlPRgZPpvfeC7Rq5Z/vMXQo8OKLfCEdO1ZBs4hIWbdlC/Dee9ze5sYbvfva4cOBatVYEyjhYfNmriIPHMjg2LnNUEQEA+cRI7jNTlSU/44hM5OZcrVqcYvM8qJ9e2DTJpY89O7N5k+jRrGWtmlTbvnZtStXHm+9NdhHW/bVrw9s2MBA+b//5QRNtWoMpCdODG5NeKVK3CZsw4bSb9m0aBGzh7p1882xFeT++/nWj6vMIdkl2xjzEIBe1trHHR8PBtDVWvtUAY99BMBTALpbazM9Pa+6ZPvR+++zccR33/n3RGotG7387W9MBYqKCt99+0REwpm1TO+MjWWTmerVvX+OZ55hymx8vG9LgSTwrOWkyZEjHBPVql34mNxcriatXs36dW8nWYpj/Hg2ZIqOLp9bWmZn85puwQJ2H46IYLDzxBNcFFFzL9+zlrdQup49cwZo0gS4/voLMw6KKyGBK9UTJ/q/03fHjiwvcAT45aVLdgIA9w3HGjnuO48xpgeAvwLoXVSwLH6Ul8emEZ06Ad27+/d7OTd/HzeOqxKDBnF1W0REypalS3lxM3lyyYJlgBOnWVnsZCxl27Jl3A948uSCg2WAXdIXLuTewQ8/DPg6c3DDBm6v9Pjj5TNYBlhnOnIkVz3T09mfxjl5oGDZP4wJrWAZYHnMqFHAF19wC6yS+PhjTgT4ozt2fg88wAyJEyf88vShusJ8MYBYAHeAgfJmAIOstXvcHtMRbPbVy1p7sDjPqxVmP/niC55I//lPBrCBYC1f1CZMYKC+YgXQqJhl7OnpvEg7fpxNY269FYiM9OfRioiIu5wcbplStSqwdWvp9vu89VauMB88qAv6sspapgTn5vLivKjgYds2rnx1786mSb4INn75hdcTERF8fl0XSHmXksJV5j59SrZ/cocO7PIdiD4Te/cCbdow83TEiPKxwmytzQHTrNcC2AdgmbV2jzHmVWNMb8fD3gBwKYBPjDHbjTHRQTpcmTGDwWq/foH7nsZwlXnFCuDAAf5TLl7sua55yxbWQderx06Qw4ezTqdJE6b0heDkkYhIWFq0iB2IX321dMEywPP64cNq/lWWrVnDQHncuOIFvx07Mtvsq6980yk9O5tdfFNTuSuHgmUR1vE/+SSwZAlw9Kh3XxsTw0Zcgwf759jyu/pqbof16ad+efqQDJgBwFq72lp7pbW2ubX2dcd9E6210Y73e1hr61prOzhuvT0/o/jF1q3cP3PUKP+0iy9K795sEtKyJVM+unVjCkhiIlPF4+PZebV7dzawWL6cj1u3jl1Zv/0W6NyZaX3PPaegWUTE37KzGeR06sRzeGk9+CBXBRcsKP1zSXBMn85ax4EDi/81TzzB5lwvvwx8/nnJv7e17Ej8/fes3W3fvuTPJRJuRo/mJNabb3r3de+9xwwif3bHdmcMXwvWrweSk33+9CEbMEsZMW0aZ2KfeCJ4x9CqFfCf/3CVOCmJaeF163LVonFj7u95/DjrrBMS+ILYowdwxRXAbbdxe4SnnuLnZ80K3s8hIlIeLFzIxk4vv+ybFOrISDaCWrqUHY6lbNm7lw28nn7au4l3Y3hR3qULL8r37vX+e1vLJl8ffsju3IG6uBcpKxo25PZqc+eybKE4fv2VWZ+DBpW8P0VJPPAAyzqifZ90HJI1zP6iGmYfO3iQwerYsdw3MRTk5gL/+x+bVaSmApdfzjqnjh09X5jl5XFPwZUr+bWdOgXumEVEyousLKbNXXYZz9W+qjn+8kuW2ixfzosmKTtGj2bdYUJCyTqdx8czaK5UiatLzZoV7+tyc9llfc4crjBHRakGXqQgsbGu6/1p04p+/KxZwJgxTMvu3Nn/x+dkLbdAa9cO5osvfFrDrIBZSm74cKbAHTvGuuCyLiUFaNcOqFGD+wwGI8VcRCScffghOxCvWgXcfbfvnjcnhxlFXbv6dS9O8bFz57jvb8+erJMsqe3bgTvuAC69lI1I27Xz/PjERNa+r1kDPP88MGVK6HUpFgklgwcDn3zCvkFNmhT+uHPnuO9yy5acwAq0MWOAqCiYrKzwb/olZcCxY8D8+Ux3DodgGWBzg3ffZVrXnDnBPhoRkfCSk8Mtg7p04WqwL118MftTrF4NnDrl2+cW/1m+nNlgpS3r6tAB+Ppr1sdffz3wzjscb/llZvL1vU0b9jB5912umClYFvFs8mT+n4wf7/lxH30E/PwzSxyCoV8/ZjL5mFaYy5u0NBbDX3IJ05VLmn40aBBn8WNji7+dU1lgLdCrF9OyDx4sWXqYiIhcaOFCNmlasYLblPjajh0MnP72N+4jK6HPuSVYbKxvgtYTJ1hvuW4dV6779mUJQFYWu3BHRzNA796d46Rt29J/T5HyYuJENmz8+mtmdOSXmsrU7ZYtuad5MEocrAWaN4c5elQrzOKlxETODLVpw3Tj5s25KlyvHut29u3z7vl++IGdqEePDq9gGeA/91tvca/m4tRpiIhI0XJzgddfB665xjedsQvSvj2ff+FC/zy/+NaBA2z29fjjvlvhrV+fjTxXrmQvknnzWKf83HMsA+jVi8H0d98pWBbx1vjxnIB69NGCO1GPG8cMn3feCV4/AGO4qOfrp9UKcxg7d46dn6dMAU6f5ozqnXcyyE1LAzZu5IvKuXOs5XnjDaB2bc/PmZHBBloZGcDu3eG7V+GQIcCyZezkWr9+sI9GRKRsW7oUGDCA59V+/fz3fWbMYHC0fz8v7CR0jR3LCerjx/1X2pWXxwv4ihW5YKCmXiKls20bt3Dt0oU9AKpV4/3vvcdFuOefD/6CU0ICTKNGavpVUuUqYD58mBcl27Yx9W3KFG7qnV9SEgf27Nl8MZk5k9sqFPSiYi1nlRYt4gxuz55+/zGC5vBhXmyNHMnfjYiIlExuLlOlc3M50erPetETJzgpPGEC8Npr/vs+UjrZ2fw73XAD8NlnwT4aEfHG8uXAww9ze9bHHmOm6uLF7E0RHc2eEkFmjFFKthRh+XKmIh09ylqxFSsKDpYB1ui++SawdStTtR99FLjlFnacdJedzb2KFy1i/UI4B8sAfxdDh3Jv5/j4YB+NiEjZtWgRA+WXXvJ/c6X69ZlJtWgRVxclNK1axXKxYcOCfSQi4q0HH2Q5RY0abO4VHc2V5RUrQiJY9getMIeTrCymOL39NnDddUyB+93viv/1ubnsbjdhAmsT7rgDuO02Pu+SJaw3eu45YPr08pHWdOwYcOWVfEF/991gH42ISNmTkcHzaP36bKYYiNeOxYvZMXv9epYiSei57z5gyxYgLi5sL7BFyoX0dCAiIuT+j7XCLAU7dgy46SYGy6NGsTudN8EyAFSowOYbsbFcCTh6FPjrX4FXXmGNQnQ065zLQ7AM8Pf3+OPcN/To0WAfjYhI2TNzJrN0Avnacf/93I93wYLAfD/xzs8/c/uvIUNC7iJbRLwUGVku/o+1whwOoqP5wpOXB8ydy1QJXzl9mv8IVar47jnLkoQEpmcPGsTfrYhIeXLqFBu7/PQT96q/7bbCS3zy27+ftcv33Qd88olfD/MCf/wjy5NOngSqVg3s9xbPpk5lJltsLLefERHxMa0wi0tmJleT+/QBmjVjHbIvg2WAs/TlNVgGgIYNgREjgPnz+eIuIlIeZGYCL7zAxkyPPso6tZEjgdatuXfuf//r+evPnePXXXIJtxgJtMceY6rg4sWB/95SOGs5+XzLLQqWRaTMUMBcVu3eDVx/vSsFe9MmroSK740fz/qMl14K9pGIiPjfyZPAzTdzd4V+/VhrmpnJ0p833mA/i27dgMGDmV6bX04OA9bNm4EPPvDflkGe3HQT92R+5x0GaRIavvkGOHiQ5U4iImWEAuayJiMDePFF7oV8/Djw+efcx7By5WAfWfi6/HJOSixZAuzcGeyjERHxn6QkNnzcu5cpzQsXcteFSpWAJk3Y+DE2lim1y5axoddrrzHIBpiGfc89wMcfM/W2b9/g/BzGAE8/zXP2hg3BOQa50OzZfE3t3z/YRyIiUmyqYS4ua1nPevQou0ZHRDBVrWFDNsvyt8xMdrCeNIkz+o8+CsyYAVx2mf+/twCpqUx779wZWLeu/DQ+E5HyIzUVuP12Br1r1jD12pMjR4Bnn+VWIgDTr8+c4duZM4Hhw/19xJ6dPcvX6TvuCHwNtVzo0CFOsPzf/7GZqIiIn/i6hjn825qV1r59TOlasQI4ceLCz0dEAO3bM5Dq2JENTtq08U3db24uEBPDF/p587jV0403cub+lltK//xSfDVrAq+/zhq+Zcu4YbtIYTIzmf2xYQPw22/AFVcAd93FMgpNtkgo+u03oFcvrixHRxcdLAOcRPzsM75OrlzJ18hmzZjGHYw07PyqVgWefJJbIe7fD7RqFewjKt+iorjA8Kc/BftIRES8ohXmwqSlMZ1r0SIGv3ffzZn3li358enTTInet4/NtrZtY4MRgC8IrVoBbdvyQrlRI6BBA6B6dTbRcjbSys113c6cYSpcUhLre3btYlOV5GR2qe7TB/jzn3kMuuAOjtxcoGtXrvDv38+ttkTyW7aMq27x8fxfr1mT2Sl5eZxce+MN4M47g32UUlalp/P1p2ZN3zVkPHOGwfKPPzINu3dv3zxvKEhKApo25c+kBmDBk5zMrRr79OF1lYiIH2mFORBiYriC+NNPwLhxvPgtKvU5L4/paTt2ANu387Z5M2ffs7K8+/4VK3LbjnvvBX7/e6BnT6B27ZL+NOIrFSoAc+ZwlXDkSNb2iTjl5XGSbc4coEsXNjvq0YPjJj2dgfRrr/H/ecgQ1vJVrx7so5ZQZy3w1VfMMvr2WyAx0fW5Vq04+TJ0KDOcSjKZeuYMJ4Q3bWKfhnAKlgGgTh3+X06bxrrrdu2CfUTl08yZHGsTJgT7SEREvKYVZnfW8iL2+eeZTrZkCXDDDaX7ptZyZjUhwbUycOYMm3dVqMDbRRex5qtOHd4aNGDQLKHplVeAl1/mBeyQIcE+GgkF1jID5O9/5wTb1KnMDMkvM5N9CKZMYebJ3Lmsr/S33Fymy6ak8HxTuzbPcd4GWHl5wOHDzKjJzATq1weuvVaBvz+cPcuVuLfeYiZTnTqcQG3XDoiMZOC8eTPw9df8W3TtyoaQ99xT/L9raiqbcm3YAPzzn8CAAX79kYImOZm1s23aAP/+t7K0As25unzPPbyuEhHxM1+vMCtgdkpO5jYY0dGcYf/oI6BWrcAeoJQNublcOdy0iSs/3bsH+4gkmKzlClZUFLcgmzy56AvyH3/kZEtsrKvG0pcp/r/9BvznPwwO/v1vlo1kZ5//mIgIlpi0aXP+rVkzTuSdOcPj27WLX791KzNnnKUnThddxAvhv/yFJSPh4OxZvhZ8+y1/7l9+YaZQZCRQowa7/Napw7fO952ZBO4To6dP8/6aNYHGjbmHcZs2njOWjhzhxMsHH3CCo2NHYPRoZj1VqnTh41NTGezOmMFtnzp04CreAw8UPGnjtGsX8NBDbGQ5fz4wcGApf2kh7sMPuZXR3LnAH/8Y7KMpX556Cnj3XXYsb9Mm2EcjIuWAAuZSKDRg3riRM+u//AK8+SYvfjUDLZ6kpLAB24kTwKpVfF/KH2sZzMyeze12pk8v/rkjI4PdYmfOZFbJxIlMrS0oKPIkJ4eBknP7nA0buAKcl8dMleuuY6ZM8+YM1PLyWNd5+DBr8ffsYfmJU8WKDIIzM133Va3KQKxjR24x1LEjg8ejR9k1fv58rnjedRf3hm/RwrufIVScOMHU3blzGfjWqMEV9IYN+Xc5fZoBalISf97ERODcuQufxxhmDV1yCX/fqan8OznVrcuV4rZt+T4AxMXxtWjnTgbZffpwO7ubby7emMrOZo3u5Mmc6GjcmM2V+vZl6rYxHK8HD7JsICqKgfzy5fwe4S4vj5ObO3ey7Kply2AfUfmwfTuboo4YwQaqIiIBoIC5FC4ImLOz2fn4tdeYLrRkCWsPRYojLo71g3FxvMAO9xUaOZ+1DJJnzmTQPGNGySbafvgBGDOGq861awMPPgjcdBNTSKtX54V+ejonaX7+meUdztuxYwyAnH0SqlQBunVjF/3u3fl+1apFH0N6OtN+9+wBDhzgfdWqudJYr7zS8/Z5mZm8GJ40yXVefeaZwGy55wtJSQyU58zh73LgQGDYMP4eL7qo8K+zlivJiYn8O0VG8hYRcf5YcG5LuGcPsHv3+W/PnuVjqlfnZETPnsAjjzBlvyRyczmJN2sWsH6967nr1nUF+xdfzK0Jp00rX1sTxsXxd1y/PieWatQI9hGFt6wsnsuOHuUkTs2awT4iESknFDCXwnkB88aNTCGMieHFSVSUuh6L95KSgPvvZ3r2wIHMUGjQINhHJf5mLdOvp09nRsrs2aXLSnFv7LRyJYMwT+rU4arnFVdw9bBVK6b7dugAVK5c8uMorYQEppivWsVV7Q8+YAPDUJWWxgmPmTMZuD7yCFf9A7VCbi0nG/Lyijex4a24OO6nvHs3g/rISAaMffpw/JRHX3/NJmedOvF3oyDOf8aP56TMv/7FiUARkQBRwFwKXdq2tTFjxrCRynffMbCZNQvo3z/YhyZlWU4O0yAnTWI668iRbADVrFmwj0z8wVo2BnzzTf6do6J8W8KRk8PVmKNHWYtcoQK3p6pRg0FOvXrBDYqLYi3PsaNGceV61Cimmxd3QvL0ae4FnJLCYLJ6df7MLVp4rsn1RmIi64Rnz2b/ioce4v+v9uktH1as4Ot+kybAp5+qc7Y/LFzILIbhw/m/JiISQAqYS6GLMTYGYCAzfDgbUVxySbAPS8LFkSMMDD7+mCtGPXpwVv2++8rvak64yctjZso773Bi5O23PafslmeJicALL7BcITKS9dl9+7KeMTKSv8v4eKYl79zJuutt25hiXtDrUqVKXEXv1ImlM126ANdcU/zJg9RU4MsvudXf558zXfSuuxgod+7s0x9dyoCNGzlRcuoU/6fHjmUDNym9pUuZrXHzzVzFD+UJPhEJSwqYS6FL06Y2ZuVK1uSpqZf4S3w8O7IuXMjGSgAv7G+8kbdu3VgzX1bqO4XS04HBgxlsjRnDFWadR4q2dSvruz/5xNWpu3Ll85uKAVzt69iRt/btGbxUrsxV9uPH2dV5505gyxYGOQAzOtq2ZRDduDFrwCtXZh3vmTNMET9+HNixg4E4wOft148TplpRLt+SkoBx47grRpUqLKsZMAC47TZt7VgS2dmcgJo0ibXLX3yhLedEJCgUMJdCkfswi/iStWyk9PnnLAH44QemmwK8OLvySl6wN27MJjT16zP1tHp1V/OgyEhmQZSFwMy553hKiuv266+sDXXeMjLOf+tseFSxIlcQIyLYnKhePf4+WrZkRkiwL17Xr2dWypEjrHdVJ33vpaUB33/P/4nkZNffuk0bBr3F3cbPWtbmxsTwtmULV6adQbS7qlWZ3dG2LVeknY3QNFkl7vbv5//1xx/zHF2rFvdHv/12vm3RQv/vnljL4PjFFzmpNWQI8N57fJ0TEQkCBcyloIBZgio3l6tkMTG8QNu3jx2JExIK3prGyX2LGuft0ksvvK+wz1erxo67TZr4prHd2bNcrTtwgD+H821srGtCwJMqVRjIVK3q6iacnc0U2TNnGGS7u/hiXrBefTVw1VXn3/y5V3peHjvpzpzJPXmbNWN6sfbdDk05OQzEs7M5ZqpU4eSTAh0prowMYO1apu1/8w3PzQAn77p1461rV06+hFo5l7WcpExM5Lk0O5uvOVWrul4Tatb03eRjdjazR1av5j7ghw9z67o33mDphYhIEJWbgNkY0wvAbAAVAHxgrZ2a7/OVASwA0BlAMoCHrbXHPD2nAmYJSdYy7fTECeDkSb6fns6b8/3TpxlMFnZz/7yn/+kaNdhZuUkTBtENG7re1qzJdNZKlXjhmJbGAOTIEV4MHT7M4Nh9z15j+HxXXcXV8mbNuE1NrVq8Va/OizVncFylStE1v5mZvOhLSGAQvn+/a4Lh8GFXWi/AbtFNmvDtZZe5UnIrVXLdnKuJxrhu7h/n/1xKCrsKf/89/x61anGLpLFj/dPJWERCj7XAoUMMnDdu5LZvhw7xcxUquDIjWrfm+1ddxfNoZKR/Jmmys3lOjIvjOTgu7sL3i+quD/A8Wa+e63b55Txv1q7tOoc6z9PG8Hycmura1i42lhOm27Zx8tQYrsQPGcJ09mBnA4mIoJwEzMaYCgBiAdwJIB7AZgADrbV73R4zAsA11to/GWMGAOhrrX3Y0/MqYJawZy1Xq92D6bQ01nH+9JPrFhfHi6/k5OI9b40aXD1wruy2asW3LVsGNojMyWH36AMHXLfjx1mLeOqUq7Oyc1/ikmrWDLj2WqB3b27BE2qrSSISeKdOAf/7H4PnmBh2c3efQAR4PmzQgBOQ+TOBKlZk9kOFCnzrvDk/zss7/9ydksLz9M8/cxIx//VanTqcsHROgl5xhauLfsWKDHozMlyTqsnJnAR0Ts6ePAn88gsfU1z16rGcqH17NvW65RaWVoiIhJDyEjBfD+Bla+3vHR9PAABr7RS3x6x1POYHY8zFAE4CqGM9/EAKmEXyycjgBVl8PFezs7IYcEZEcHW4Zk02KPNn6rM/WMvgOiuLF6HWum7Oz+d/3/mxM4VRRKQo6emukpQTJxjcnjjB0hL37J+zZ3lOKuzmFBHhCrCdW8k1bMi08MaNXYFx48a+m6zMyGAw7bxlZvJ8mJfHLJ2aNXmrW5cr6CIiIc7XAbOPNrX0uYYAjrt9HA+ga2GPsdbmGGPSANQGUEDnFxEpUEQE64NbtAj2kfiWMVxhUXqgiPhTZCSzUa69tnTPk5fHt8HYpi4igqU5jRoF/nuLiJQBoRow+4wxZjiA4Y4PM40xu4N5PCI+chk0OSThQWNZwoHGsYQLjWUJB1f58slCNWBOANDY7eNGjvsKeky8IyW7Otj86zzW2vcBvA8AxpgYXy7PiwSLxrKEC41lCQcaxxIuNJYlHBhjfFqDG4Tcn2LZDKClMaapMaYSgAEAovM9JhrAEMf7DwH41lP9soiIiIiIiIg3QnKF2VGT/BSAteC2UnOttXuMMa8CiLHWRgP4EMBCY8whAClgUC0iIiIiIiLiEyEZMAOAtXY1gNX57pvo9v45AP28fNr3fXBoIqFAY1nChcayhAONYwkXGssSDnw6jkNyWykRERERERGRYAvVGmYRERERERGRoCrzAbMxZq4xJrGg7aKMMc8aY6wx5jLHx8YY87Yx5pAxZqcxppPbY4cYYw46bkPyP5eIv3k5llsZY34wxmQaY57L99hexpgDjnE+PlDHLwJ4PY7/4DgX7zLGbDLGtHd7rMaxBJWXY7mPYyxvN8bEGGNucnusri8kaLwZx273X2uMyTHGPOR2n8axBJWX5+RbjTFpjnPydmPMRLfHen19UeYDZgDzAPTKf6cxpjGAngDi3O6+C0BLx204gHcdj60F4CUAXQFcB+AlY0xNvx61yIXmofhjOQXAMwDezPfYCgCiwLHeGsBAY0xrPx2vSEHmofjj+CiA7tbadgAmwbUFoMaxhIJ5KP5Y/gZAe2ttBwCPAfjA8VhdX0iwzUPxx7Hz/DsNwFdu92kcSyiYBy/GMoAN1toOjturjseW6PqizAfM1trvweAhv1kAngfgXqTdB8ACSz8CqGGMqQ/g9wDWWWtTrLWpANahgD+IiD95M5attYnW2s0AsvM99joAh6y1R6y1WQCWgONeJCC8HMebHOdcAPgRQCPH+xrHEnRejuXTbltbXuL2OV1fSFB5eZ0MAE8DWA4g0e0+jWMJuhKM5YKU6PqizAfMBTHG9AGQYK3dke9TDQEcd/s43nFfYfeLBJWHsVwYjWUJOcUcx8MArHG8r3EsIcnTWDbG9DXG7AewClxlBjSWJQQVNo6NMQ0B9IUjA9ONxrGEpCKuL643xuwwxqwxxrRx3FeisRyy20qVlDGmKoAXwKV5kTJLY1nCQXHGsTHmNjBgvqmwx4gEW1Fj2Vr7GYDPjDG3gCUGPQJ4eCLFUsQ4fgvAOGttnjEmoMcl4q0ixvJWAE2staeNMXcDWAGW5JZIOK4wNwfQFMAOY8wxMMVvqzGmHoAEAI3dHtvIcV9h94sEk6exXBiNZQk1HsexMeYasN6zj7U22fE1GscSiop1TnakDTZzNJ/RWJZQ42kcdwGwxHH/QwDmGGPuh8axhKZCx7K19jdr7WkAsNauBlCxNOfksFthttbuAnC582PHL7CLtfaUMSYawFPGmCVg44I0a+0JY8xaAJPdGhj0BDAhwIcuch5PY9nDl20G0NIY0xQ8AQwAMMifxyniSRHn5CsAfApgsLU21u3LNI4l5BQxllsAOGyttYY7cFQGkAxA1xcSUoq4tmjqdv88AF9Ya1c4mn5pHEtIKeKcXA/AL45z8nXgInEygF9RguuLMr/CbIz5GMAPAK4yxsQbY4Z5ePhqAEcAHALwDwAjAMBamwKmT2123F513CcSMN6MZWNMPWNMPIAxAF50PL6atTYHwFPgRdo+AMustXsCcfwigNfn5IkAaoOrGNuNMTEAoHEsocDLsfwggN3GmO1gB9aHHQ1GdX0hQeXlOC6QxrGEAi/H8kPgOXkHgLcBDHCck0t0fWFcTR1FRERERERExKnMrzCLiIiIiIiI+IMCZhEREREREZECKGAWERERERERKYACZhEREREREZECKGAWERERERERKYACZhEREREREZECKGAWERERERERKYACZhEREREREZEC/D+XtiEumPJMpgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "tp1, tp2 = 1400, 1450\n", "unit_name = '8-2'\n", "\n", "with h5py.File(h5name, 'r') as f:\n", " tl = np.array(f['processed']['timeline'])\n", " i_rate = np.array(f['units'][unit_name][H5NAMES.inst_rate['name']])\n", " s_times = np.array(f['units'][unit_name][H5NAMES.spike_times['name']])\n", "\n", "i_rate = instantaneous_rate(s_times, tl[:, 0], k_width=50)\n", " \n", "fig = plt.figure(figsize=(16, 6))\n", "\n", "ax = fig.add_subplot(311)\n", "ax.plot(tl[:, 0], i_rate)\n", "ax.set_xlim(tp1, tp2)\n", "\n", "# head direction\n", "hd = np.rad2deg(tl[:, 4])\n", "ax = fig.add_subplot(312)\n", "ax.plot(tl[:, 0], hd)\n", "ax.set_xlim(tp1, tp2)\n", "ax.set_ylim(-180, 180)\n", "ax.set_ylabel('HD, deg', fontsize=14)\n", "\n", "# speed\n", "width = 300 # 100 points ~= 1 sec with at 100Hz\n", "kernel = signal.gaussian(width, std=(width) / 7.2)\n", "s_smooth = np.convolve(tl[:, 3], kernel, 'same') / kernel.sum()\n", "\n", "ax = fig.add_subplot(313)\n", "ax.plot(tl[:, 0], s_smooth, color='red')\n", "ax.set_xlim(tp1, tp2)\n", "ax.set_ylim(0, .5)" ] }, { "cell_type": "code", "execution_count": null, "id": "5cbf9a58", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }