{ "cells": [ { "cell_type": "code", "execution_count": 2, "id": "06db5f6c", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import os, h5py\n", "import matplotlib.pyplot as plt\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 5, "id": "b5a09f6a", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "6095" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "source = '/home/sobolev/nevermind/Andrey/data'\n", "\n", "session = '009266_hippoSIT_2023-04-26_08-20-17' # sample session with some PPC units\n", "#session = '009266_hippoSIT_2023-04-28_09-04-09'\n", "#session = '009266_hippoSIT_2023-05-04_09-11-06' # coma session with baseline AEPs\n", "animal = session.split('_')[0]\n", "sessionpath = os.path.join(source, animal, session)\n", "aeps_file = os.path.join(sessionpath, 'AEPs.h5')\n", "\n", "with h5py.File(aeps_file, 'r') as f:\n", " aeps = np.array(f['aeps'])\n", " aeps_events = np.array(f['aeps_events'])\n", " \n", "t_N1 = 20 # all in ms\n", "N1_l, N1_r = 15, 25\n", "P1_l, P1_r = 30, 75\n", "\n", "len(aeps)" ] }, { "cell_type": "markdown", "id": "dcdd5bc8", "metadata": {}, "source": [ "### N1 / P1 Integrals" ] }, { "cell_type": "code", "execution_count": 55, "id": "358f920f", "metadata": {}, "outputs": [], "source": [ "# normalize to the LFP N1 start\n", "aeps_N1 = np.array([aep - aep[25] for aep in aeps])" ] }, { "cell_type": "code", "execution_count": 56, "id": "08017857", "metadata": {}, "outputs": [], "source": [ "aeps_N1_areas = aeps[:, N1_l:N1_r].sum(axis=1)\n", "aeps_P1_areas = aeps[:, P1_l:P1_r].sum(axis=1)" ] }, { "cell_type": "code", "execution_count": 57, "id": "ba915a39", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWsAAAEDCAYAAADz4SVPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAPtUlEQVR4nO3dbYxc5XnG8euqDbgKpH7bUgtDbFc1bZIqQDbULmnF0jhxHBTaqpVAbZQ2RKvmBREaFOHypc0XElLVJGorWAFpopA0DjENspRQp1maRHJM1gQDNtgxhlBbtF4cUcKXSiZ3P8xjdnYzszOzO2dm79n/Txpx5jlndu+Hc/bax895WUeEAAAL2y/1uwAAQGuENQAkQFgDQAKENQAkQFgDQAKENQAkUFlY277X9knbT3bp671q+7HyerAbXxMAsnBV11nb/n1Jr0j6YkS8uQtf75WIOHf+lQFAPpWNrCPiu5J+Wt9m+9dtf8v2ftvfs/2bVX1/ABgkvZ6zHpN0Q0S8VdLNkv65g88usz1h+we2/7CS6gBggVraq29k+1xJvyvpa7bPNJ9T1v2xpE82+NiJiHhXWX5DRJywvUHSd2w/ERHPVF03ACwEPQtr1UbxL0XEJTNXRMQuSbtm+3BEnCj/PWb7YUmXSiKsASwKPZsGiYiXJT1r+08lyTVvaeeztlfYPjMKXy3pCkmHKisWABaYKi/d+4qkvZIutn3c9vWS/kzS9bYPSDoo6Zo2v9xvSZoonxuX9KmIIKwBLBqVXboHAOge7mAEgAQqOcG4evXqWLduXRVfGgAG0v79+1+MiKFm6ysJ63Xr1mliYqKKLw0AA8n2T2Zb31ZY235O0s8kvSrpdEQMz780AEC7OhlZj0TEi5VVAgBoihOMAJBAu2Edkv69PIBptNEGtkfLszsmJicnu1chAKDtsH57RFwm6d2SPlIefzpNRIxFxHBEDA8NNT2hCQCYg7bCuu65HCclPSDp8iqLAgBM1zKsbb/O9nlnliW9U1JX/voLAKA97VwNcr6kB8pjTZdK+nJEfKvSqgAA07QM64g4Jqmtp+MBAKrRy+dZA/M3flvnnxnZ3v06gB7jOmsASICwBoAEmAbBwjSX6Q5ggDGyBoAECGsASICwBoAECGsASICwBoAEuBoEfbVjz5GG7TdxZALTMLIGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAR4ajAWpL3HTjVs37xhVY8rARYGRtYAkABhDQAJENYAkABz1kil2Vy2xHw2BhsjawBIgLAGgAQIawBIoO2wtr3E9o9s766yIADAL+pkZH2jpKeqKgQA0FxbYW17raT3SLq72nIAAI20e+neHZI+Iem86koB5qfpLeojPS4EqEDLkbXtqyWdjIj9LbYbtT1he2JycrJrBQIA2htZXyHpvba3SVom6fW2vxQRf16/UUSMSRqTpOHh4eh6pUhtx54j/S4BSK3lyDoitkfE2ohYJ+laSd+ZGdQAgGpxuzkG3/htjdtHtve2DmAeOgrriHhY0sOVVAIAaIo7GAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAcIaABIgrAEgAf6sFwbe3mOnGrZvHulxIcA8MLIGgAQIawBIgGkQ9NWm58f6XQKQAiNrAEiAsAaABAhrAEiAsAaABAhrAEiAq0HQVXvvubnxiotGe1sIMGAYWQNAAoQ1ACRAWANAAoQ1ACRAWANAAoQ1ACRAWANAAi3D2vYy24/YPmD7oO2/60VhAIAp7dwU83+SroqIV2yfJen7tr8ZET+ouDYAQNEyrCMiJL1S3p5VXlFlUQCA6dqas7a9xPZjkk5K2hMR+xpsM2p7wvbE5ORkl8sEgMWtrbCOiFcj4hJJayVdbvvNDbYZi4jhiBgeGhrqcpkAsLh1dDVIRLwkaVzS1kqqAQA01M7VIEO2l5flX5a0RdLTFdcFAKjTztUgayR9wfYS1cJ9Z0TsrrYsDBr+MC4wP+1cDfK4pEt7UAsAoAnuYASABAhrAEiAsAaABAhrAEiAsAaABPjr5li0duw50rD9pi0be1wJ0BojawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAS4KQYda3YziSRt6mEdwGLCyBoAEiCsASABwhoAEiCsASABwhoAEuBqEHSMv1QO9B4jawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgLAGgAQIawBIgGeDYNFq/oyTv+9pHUA7GFkDQAItw9r2hbbHbR+yfdD2jb0oDAAwpZ1pkNOSPh4Rj9o+T9J+23si4lDFtQH9MX5b83Uj23tXB1Cn5cg6Il6IiEfL8s8kPSXpgqoLAwBM6WjO2vY6SZdK2tdg3ajtCdsTk5OTXSoPACB1ENa2z5X0dUkfi4iXZ66PiLGIGI6I4aGhoW7WCACLXlthbfss1YL6vojYVW1JAICZ2rkaxJLukfRURPxD9SUBAGZqZ2R9haT3SbrK9mPlta3iugAAdVpeuhcR35fkHtQCAGiCOxgBIAHCGgASIKwBIAHCGgASIKwBIAHCGgASIKwBIAHCGgASIKwBIAH+BiMww95jp5qu2zzSw0KAOoysASABRtZoaseeIw3bN/W4DgCMrAEgBcIaABJgGgRNbXp+rN8lACgYWQNAAoQ1ACRAWANAAoQ1ACRAWANAAoQ1ACRAWANAAoQ1ACRAWANAAoQ1ACRAWANAAjwbBOhAs8fG3rRlY48rwWLDyBoAEiCsASABwhoAEiCsASABwhoAEiCsASCBlmFt+17bJ20/2YuCAAC/qJ2R9b9I2lpxHQCAWbQM64j4rqSf9qAWAEATXbuD0faopFFJuuiii7r1ZdEDe++5ud8lAGihaycYI2IsIoYjYnhoaKhbXxYAIK4GAYAUCGsASKCdS/e+ImmvpIttH7d9ffVlAQDqtTzBGBHX9aIQAEBzTIMAQAKENQAkQFgDQAL8WS+gA5ueH2u8YnxV4/aR7dUVg0WFkTUAJEBYA0ACTIMsJuO39bsCAHPEyBoAEiCsASABwhoAEiCsASABwhoAEiCsASABLt0DumDvsVMN2zeP9LgQDCzCehFpFigAFj6mQQAgAcIaABIgrAEgAeasgSo1ex4Lj05FhxhZA0AChDUAJMA0yCDiUajAwCGsgQpxswy6hWkQAEiAsAaABAhrAEiAOesBxDNAgMHDyBoAEmBkDfTBjj1HGrbftGVjjytBFoysASABwhoAEmAaJCvuUkxt0/NjjVeMr2r+IR7+tKgR1klxxcdgmm2/ctfj4sY0CAAk0NbI2vZWSZ+VtETS3RHxqUqrWoyY1gAwi5ZhbXuJpH+StEXScUk/tP1gRByqujgw3YEpXO63uLUzsr5c0tGIOCZJtv9V0jWSCOsmmv1QSc1/sAhlzFWnIT7b8dmpTr/HXH6xdKve7L/UHBGzb2D/iaStEfHB8v59kn4nIj46Y7tRSaPl7cWSDs+xptWSXpzjZxeaQenLoPRDoi8L0aD0Q5pfX94QEUPNVnbtapCIGJPU5Hqk9tmeiIjhLpTUd4PSl0Hph0RfFqJB6YdUbV/auRrkhKQL696vLW0AgB5pJ6x/KOk3bK+3fbakayU9WG1ZAIB6LadBIuK07Y9Keki1S/fujYiDFdY076mUBWRQ+jIo/ZDoy0I0KP2QKuxLyxOMAID+4w5GAEiAsAaABHoS1rZvsP207YO2b69r3277qO3Dtt9V1761tB21fUtd+3rb+0r7V8sJT9k+p7w/Wtavq7g/H7cdtleX97b9ufL9H7d9Wd2277f94/J6f137W20/UT7zOdsu7Stt7ynb77G9ooL6P1P2x+O2H7C9vG5dyn3SSrP6+832hbbHbR8qPx83lvaGx0E3j7WK+rPE9o9s7y7vOz4+Oj0GK+rHctv3l5+Tp2xv7vs+iYhKX5JGJH1b0jnl/a+W/75R0gFJ50haL+kZ1U5gLinLGySdXbZ5Y/nMTknXluU7JX2oLH9Y0p1l+VpJX62wPxeqdrL1J5JWl7Ztkr4pyZI2SdpX2ldKOlb+u6IsryjrHinbunz23aX9dkm3lOVbJH26gj68U9LSsvzpM98j6z5po79N6+/3S9IaSZeV5fMkHSn7oeFx0M1jraL+/LWkL0vaPZfjYy7HYEX9+IKkD5blsyUt7/c+6cXBuFPSOxq0b5e0ve79Q5I2l9dDM7crnXpRUyHz2nZnPluWl5btXFF/7pf0FknPaSqs75J0Xd02h8sP4XWS7qprv6u0rZH0dF37a9ud+WxZXiPpcMX7548k3Zd5n7TRx4b196OWNmr9hmrP4Wl4HHTzWKug9rWS/kPSVZJ2z+X46PQYrKgfvyLp2ZnHa7/3SS+mQTZK+r3yT53/tP220n6BpP+q2+54aWvWvkrSSxFxekb7tK9V1v9v2b6rbF8j6UREHJixqtO+XFCWZ7ZL0vkR8UJZ/m9J53en+qY+oNpvdinhPmlTs/oXlDIVcKmkfWp+HHTzWOu2OyR9QtLPy/u5HB+d9q8K6yVNSvp8mdK52/br1Od90pXbzW1/W9KvNVh1a/keK1Ub8r9N0k7bG7rxfavQoi9/o9oUQk9ERNie07WVs/UjIr5RtrlV0mlJ9829SnSD7XMlfV3SxyLi5fopzPkcB71i+2pJJyNiv+0r+1zOfC2VdJmkGyJin+3Pqjbt8Zp+7JOuhHVEvKPZOtsfkrQrauP9R2z/XLWHncx2G3uj9lOSltteWn4T129/5msdt71UtX/GzOkxds36Yvu3VfuNe6D8IK2V9Kjty2fpywlJV85of7i0r22wvST9j+01EfGC7TWSTnazH3X9+QtJV0v6g7JvNEs/1KS9J/ukCxb0IxNsn6VaUN8XEbtKc7PjoJvHWjddIem9trdJWibp9ao9A7/T46PTY7AKxyUdj4h95f39qoV1f/dJFXM+M+Z5/krSJ8vyRtX+WWBJb9L0EwnHVDuJsLQsr9fUiYQ3lc9/TdNPVny4LH9E009W7OxBv57T1Jz1ezT9BMMjpX2lanNfK8rrWUkry7qZJxi2lfbPaPpJjNsrqH2rao+4HZrRnnqfzNLfpvX3+1X2/xcl3TGjveFx0M1jrcI+XampE4wdHR9zOQYr6sP3JF1clv+27I++7pNeHIxnS/qSpCclPSrpqrp1t6p2hvew6s6GqnZ29UhZd2td+4bSyaPlIDhzhcmy8v5oWb+hB/16TlNhbdX+QMMzkp6QNFy33QdKXUcl/WVd+3D5f/KMpH/U1N2kq1Q7SfNj1a6iWVlB7UdV+6X5WHndOQj7pEWfG9bf75ekt0sKSY/X7Y9tzY6Dbh5rFfbpSk2FdcfHR6fHYEV9uETSRNkv/6Za2PZ1n3C7OQAkwB2MAJAAYQ0ACRDWAJAAYQ0ACRDWAJAAYQ0ACRDWAJDA/wNHruFu/I902gAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "_ = plt.hist(aeps_N1_areas[aeps_events[:, 1] == 1], bins=np.linspace(-60000, 60000, 50), density=True, alpha=0.5)\n", "_ = plt.hist(aeps_N1_areas[aeps_events[:, 1] == 2], bins=np.linspace(-60000, 60000, 50), density=True, alpha=0.5)" ] }, { "cell_type": "code", "execution_count": 59, "id": "543fe098", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEDCAYAAADOc0QpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAP7klEQVR4nO3df4ylVX3H8fengJhUW5EllgDrgkVb06QKE11qY8AWBdK4tqVxTSNqabax0h+k/CE1UeM/aJvWaLCSLW50jQGrte22oaVY12AT1jIYfht1QAq7pbKAQomNdptv/7jP2us4d+6dnTv3x5n3K7mZ557n3OeeOTv7mTPnee55UlVIktryY9NugCRp/Ax3SWqQ4S5JDTLcJalBhrskNchwl6QGTTXck+xJ8liSe8d0vP9Ncmf32DeOY0rSPMo0r3NP8mrgGWBvVf3cGI73TFU9Z/0tk6T5NtWRe1XdCjzZX5bkRUn+KckdSb6U5Gem1DxJmluzOOe+G/i9qjoXuAr4izW89tlJFpMcSPKGDWmdJM2B46fdgH5JngP8AvCZJEeLT+z2/RrwvhVedqiqXtdtv7CqDiU5C/hCknuq6oGNbrckzZqZCnd6f0l8p6petnxHVX0O+NxqL66qQ93XB5N8EXg5YLhL2nRmalqmqp4GvpnkNwDS8/OjvDbJSUmOjvK3AK8C7t+wxkrSDJv2pZA3ALcBL0lyMMnlwG8Clye5C7gP2DHi4X4WWOxetx94f1UZ7pI2paleCilJ2hgzNS0jSRqPqZ1Q3bJlS23btm1aby9Jc+mOO+54vKpOGVZvauG+bds2FhcXp/X2kjSXkvz7KPWclpGkBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAbN2nrukvrtv2bl8guunmw7NHccuUtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1aGi4Jzkjyf4k9ye5L8kfrFAnST6cZCnJ3UnO2ZjmSpJGcfwIdY4Af1RVX0nyXOCOJLdU1f19dS4Gzu4erwQ+2n2VJE3B0JF7VT1aVV/ptv8L+Cpw2rJqO4C91XMAeF6SU8feWknSSNY0555kG/By4MvLdp0GPNL3/CA/+guAJLuSLCZZPHz48BqbKkka1SjTMgAkeQ7w18AfVtXTx/JmVbUb2A2wsLBQx3IMScD+a1Yuv+DqtdVf7TWaayON3JOcQC/YP1VVn1uhyiHgjL7np3dlkqQpGOVqmQAfA75aVX8+oNo+4LLuqpntwFNV9egY2ylJWoNRpmVeBbwZuCfJnV3ZHwNbAarqOuAm4BJgCfgu8Laxt1SSNLKh4V5V/wpkSJ0C3jGuRkmS1sdPqEpSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lq0MhL/kqaA6st7atNxZG7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapCXQkoT8sFbvj5w35UXvniCLdFmYLhLY7ZaiEuT4rSMJDXIcJekBjktIx0jp180ywx3aQYM+kVxpf9DdYyclpGkBhnuktQg/+iTJmT7w7sH7juwddcEW6LNwJG7JDXIcJekBhnuktQg59ylGTBwPv6skyfbEDXDkbskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0aGu5J9iR5LMm9A/afn+SpJHd2j3ePv5mSpLUY5UNMHweuBfauUudLVfUrY2mRJGndho7cq+pW4MkJtEWSNCbjWn7gvCR3Af8BXFVV961UKckuYBfA1q1bx/TWUrtue/CJFcvPc1kCDTGOcP8K8MKqeibJJcDfAmevVLGqdgO7ARYWFmoM7y1tOO+Vqnm07qtlqurpqnqm274JOCHJlnW3TJJ0zNYd7kl+Kkm67Vd0x1z5b0lJ0kQMnZZJcgNwPrAlyUHgPcAJAFV1HXAp8PYkR4D/BnZWlVMukjRFQ8O9qt40ZP+19C6VlCTNCD+hKkkNMtwlqUGGuyQ1yHCXpAYZ7pLUoHEtPyDNNT+FqtY4cpekBhnuktQgw12SGuScu7TZ7b9m5fILrp5sOzRWjtwlqUGGuyQ1yHCXpAY55y5pbZyjnwuO3CWpQY7cpTnkjbM1jOGuTcVlBrRZOC0jSQ1y5C41xOkaHeXIXZIaZLhLUoOclpG0skHXs2suOHKXpAYZ7pLUIKdlpDHb/vDuaTdBcuQuSS0y3CWpQYa7JDXIOXfpGDm3rlnmyF2SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yEsh1SRvp6fNbujIPcmeJI8luXfA/iT5cJKlJHcnOWf8zZQkrcUo0zIfBy5aZf/FwNndYxfw0fU3S5K0HkPDvapuBZ5cpcoOYG/1HACel+TUcTVQkrR24zihehrwSN/zg13Zj0iyK8liksXDhw+P4a0lSSuZ6NUyVbW7qhaqauGUU06Z5FtL0qYyjnA/BJzR9/z0rkySNCXjCPd9wGXdVTPbgaeq6tExHFeSdIyGXuee5AbgfGBLkoPAe4ATAKrqOuAm4BJgCfgu8LaNaqwkaTRDw72q3jRkfwHvGFuLJEnr5vIDktQgw12SGmS4S1KDXDhMc80FwqSVOXKXpAYZ7pLUIMNdkhpkuEtSgzyhKm0Ctz34xMB955118gRboklx5C5JDTLcJalBhrskNchwl6QGeUJVGmL7w7un3QRpzRy5S1KDHLlLm9ygyyS9RHK+OXKXpAY5cpc0HvuvWbn8gqsn2w4BjtwlqUmGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSg7wUUnPBG2FLa+PIXZIaZLhLUoMMd0lqkOEuSQ3yhKpmhidNpfFx5C5JDTLcJalBTsuoSWu9Nd6Brbs2qCXSdDhyl6QGOXKXtCJvvzffHLlLUoNGGrknuQj4EHAccH1VvX/Z/rcCfwoc6oqurarrx9hOSfNq0O33wFvwbaCh4Z7kOOAjwIXAQeD2JPuq6v5lVT9dVVdsQBslSWs0yrTMK4Clqnqwqr4P3Ajs2NhmSZLWY5RwPw14pO/5wa5suV9PcneSzyY5Y6UDJdmVZDHJ4uHDh4+huZKkUYzrapm/B26oqu8l+R3gE8Brlleqqt3AboCFhYUa03tL67bW6+KlWTfKyP0Q0D8SP53/P3EKQFU9UVXf655eD5w7nuZJko7FKOF+O3B2kjOTPAvYCezrr5Dk1L6nrwe+Or4mSpLWaui0TFUdSXIFcDO9SyH3VNV9Sd4HLFbVPuD3k7weOAI8Cbx1A9ssSRpipDn3qroJuGlZ2bv7tq8GvGBVkmaEn1CVpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDfJmHZLWxJt4zAfDXRP3wVu+Pu0mSM1zWkaSGuTIXRvGEbo0PYa7pNkz6NZ83pZvZE7LSFKDDHdJapDTMppr3kFpk3G6ZmSO3CWpQYa7JDXIcJekBhnuktQgT6hq3fywkjR7DHdJY+GCYrPFcJe0oQaFPsB5F0ywIZuM4a6R3faxq1besXXXZBsiaShPqEpSgwx3SWqQ4S5JDXLOXTNjtXViDjivL62J4a4fsto169sn2A5J62O4S5qaQYOJK02mdXPOXZIa5O/HebPB61m7PromaeDPm59qXTdH7pLUIMNdkhrktMwmNehE1iSuiDmWqR+nizaXiSxC1vgt+wz3xrkcr7Q5Ge5zZuCIxtX1JPUx3CXNDQc3ozPcp2jgBzgufPHYjjVornqcH+efxHtIqxm4HDUTuFnIoLl7mOr8/UjhnuQi4EPAccD1VfX+ZftPBPYC5wJPAG+sqofG29TZN86wljQem/UOUUPDPclxwEeAC4GDwO1J9lXV/X3VLge+XVU/nWQn8AHgjRvR4HFbayAfywnKtb7G9V0krdcoI/dXAEtV9SBAkhuBHUB/uO8A3tttfxa4NkmqqsbY1h84loB19Cyp36AR/YEjGz9dOok8yrD8TXIpcFFV/Xb3/M3AK6vqir4693Z1DnbPH+jqPL7sWLuAoxOxLwG+Nq5vZA22AI8PrbV52T/D2Uers3+GW08fvbCqThlWaaInVKtqNzDVT6MkWayqhWm2YZbZP8PZR6uzf4abRB+NsvzAIeCMvuend2Ur1klyPPCT9E6sSpKmYJRwvx04O8mZSZ4F7AT2LauzD3hLt30p8IWNmm+XJA03dFqmqo4kuQK4md6lkHuq6r4k7wMWq2of8DHgk0mWgCfp/QKYVS5Ssjr7Zzj7aHX2z3Ab3kdDT6hKkuaPS/5KUoMMd0lq0NyHe5L3JjmU5M7ucUnfvquTLCX5WpLX9ZVf1JUtJXlnX/mZSb7clX+6O4FMkhO750vd/m0T/SYnZFC/tCrJQ0nu6X5uFruy5ye5Jck3uq8ndeVJ8uGub+5Ock7fcd7S1f9Gkrf0lZ/bHX+pe20m/12uTZI9SR7rPrtytGzD+2TQe8yaAf0zmxlUVXP9oPfJ2KtWKH8pcBdwInAm8AC9E8LHddtnAc/q6ry0e81fATu77euAt3fbvwtc123vBD497e97A/pxYL+0+gAeArYsK/sT4J3d9juBD3TblwD/CITeKhBf7sqfDzzYfT2p2z6p2/dvXd10r7142t/zCH3yauAc4N5J9smg95i1x4D+mckMmvuR+yp2ADdW1feq6pvAEr2lFH6wnEJVfR+4EdjRjSBeQ2/5BIBPAG/oO9Ynuu3PAr80D6OwNVqxX6bcpmno/7de/jOwt3oOAM9LcirwOuCWqnqyqr4N3AJc1O37iao6UL3/kXv7jjWzqupWele89ZtEnwx6j5kyoH8GmWoGtRLuV3R/Fu7p+3PuNOCRvjoHu7JB5ScD36mqI8vKf+hY3f6nuvotGdQvLSvgn5Pckd7SGAAvqKpHu+3/BF7Qba/15+m0bnt5+TyaRJ8Meo95MXMZNBfhnuTzSe5d4bED+CjwIuBlwKPAn02zrZorv1hV5wAXA+9I8ur+nd3o0muF+0yiT+aw32cyg+biZh1V9cuj1Evyl8A/dE9XWzZhpfIn6P1ZeXz3m7G//tFjHUy7yyuMssxEU6rqUPf1sSR/Q+/P5W8lObWqHu2mER7rqg/qn0PA+cvKv9iVn75C/Xk0iT4Z9B4zr6q+dXR7ljJoLkbuq+l+EI76VeDoWex9wM7uLPOZwNn0TuasuJxCN1rYT2/5BOgtp/B3fcdqfXmFUZaZaEaSH0/y3KPbwGvp/ez0/1sv/xm4rLtCZDvwVDeNcDPw2iQndX+Ovxa4udv3dJLt3dzoZX3HmjeT6JNB7zHzZjaDpn32eb0P4JPAPcDdXQec2rfvXfTOSn+NvisV6J3l/3q371195Wd1nb8EfAY4sSt/dvd8qdt/1rS/7w3qyxX7pcVH9299V/e47+j3S28e81+AbwCfB57flYfeTWse6H7eFvqO9Vvdz8YS8La+8oXuP/oDwLV0nwif5QdwA72phf+hN+d7+ST6ZNB7zNpjQP/MZAa5/IAkNWjup2UkST/KcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkN+j+P8UIYK+GCVAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "_ = plt.hist(aeps_P1_areas[aeps_events[:, 1] == 1], bins=np.linspace(-60000, 150000, 50), density=True, alpha=0.5)\n", "_ = plt.hist(aeps_P1_areas[aeps_events[:, 1] == 2], bins=np.linspace(-60000, 150000, 50), density=True, alpha=0.5)" ] }, { "cell_type": "code", "execution_count": 63, "id": "57803047", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 63, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZsAAAD8CAYAAAChHgmuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACRH0lEQVR4nOz9eYxl2XWfC3577zPeOcaMnCtrZrFYLJJlUZMltQXZlNtoyniCW2rAUtuCacMSHhr9gPbQQOvBQ0Puhp9hP9hqyM+CpcazabXcttR+dNNs2RotSiyyyGJVsYaszMo55jufcQ/9xzkRFZmVmZVZGTXm+YBAROx77rn3xnDWXXv91m8J5xwNDQ0NDQ3vJvL9fgINDQ0NDR99mmDT0NDQ0PCu0wSbhoaGhoZ3nSbYNDQ0NDS86zTBpqGhoaHhXacJNg0NDQ0N7zp3HGyEEL8shNgUQrxwYO2/F0JcEUJ8s/74swdu+1tCiLNCiFeEEH/mwPrn6rWzQoi/eWD9jBDij+r1fyOECOr1sP7+bH37A/f8qhsaGhoa3lPuJrP5l8DnbrL+j5xzT9cfXwIQQjwB/ATw8fo+/0wIoYQQCvinwI8CTwA/WR8L8A/qcz0MDIGfqdd/BhjW6/+oPq6hoaGh4UPEHQcb59zvArt3ePjngS8653Ln3HngLPBd9cdZ59w551wBfBH4vBBCAH8K+PX6/r8C/NiBc/1K/fWvAz9cH9/Q0NDQ8CHBO4Rz/JwQ4qeAZ4H/zjk3BI4DXz1wzOV6DeDSDeufBZaAkXNO3+T443v3cc5pIcS4Pn77xicihPgC8AWAdrv9mccff/zeX13DRw5jHUmhkUIghcA6h3WOVuCh5Lv7PuaVV14B4LHHHntXH6eh4Z3y9a9/fds5t3LY573XYPOLwN8FXP35HwJ/+V6f1DvFOfdLwC8BPPPMM+7ZZ599v55KwweY5y4OKbSlFbz5558UmsCTfOrUwrv62D/0Qz8EwG//9m+/q4/T0PBOEUJceDfOe09qNOfchnPOOOcs8M+ptskArgAnDxx6ol671foOMBBCeDesX3eu+vZ+fXxDwztikpbEvrpuLfYVk7R8n55RQ8NHn3sKNkKIowe+/fPAnlLtN4GfqJVkZ4BHgD8GvgY8UivPAioRwW+6yg30vwA/Xt//p4HfOHCun66//nHgP7vGPbThHujFPmlprltLS0Mv9t+nZ9TQ8NHnjrfRhBD/GvghYFkIcRn4eeCHhBBPU22jvQH8VQDn3ItCiF8DXgI08LPOOVOf5+eALwMK+GXn3Iv1Q/wN4ItCiL8HPAf8i3r9XwD/TyHEWSqBwk+80xfb0ABwZrnNcxeHQJXRpKUhKTSPrb27W2gNDfcz4qOaJDQ1m4bbMUoKzm/PmaQlvdjnzHKbQSt41x+3qdk0fNARQnzdOffMYZ/3MNRoDQ0fOgatgE+deveDy41o68hKw++8svmeBrmGhvebxq6moeE9YpQUzLIS6xwLrYBCW567OGSUFO/3U2toeNdpgk1Dw3vE+e05UgqUEAghaAUercDj/Pb8/X5qDQ3vOk2waWh4j5ikJeoG84tGct1wv9AEm4aG94he7GNuEOQ0kuuG+4Um2DQ0vEecWW5jrcM4h3OVZU5SaM4st9/vp9bQ8K7TBJuGhveIQSugE/lIIRgmxb49TqNGa7gfaKTPDQ3vIZ4UdEKPH3xs9f1+Kg0N7ylNsGlo+IDwfjWaNjS8FzTbaA0NHwBGSbHvRt304DR8FGmCTUPDB4Dz2/P9vpumB6fho0izjdZwX/NB2bqapCULNzxu7CuGTWbT8BGhCTYN9y17W1etwGOhFZCWhucuDt91hZi2jucuDq8LcHtjDw4OdGt6cBo+SjTBpuG+5eDW1TQruTpK2ZkXbE5y/syTa4cScEZJwfOXR5zdnOIQ7M4LnHP7tZm9APfgSodzWzOgGXvQ8NGkCTYN9y1XRylJrtmaFezMco4NYpbbAdvz/FAynFFS8HuvbXF5mBJIwdY8Y2ua4SvJOClpDbz9TGZ3XvCpUwuc354zTAp6sc9ja00PTsNHhybYNNyXjJKCa+MUTwjy0qCEYH2cArDUDrEWvvzCOmv96B3Xcs5vzxmnJYGSrI8zQk8hhcA6xx+d3+aHP7ZGN/L3azPv19iDhob3gkaN1nBfcn57zpmlDg7BONO0QoUAro5TepHPhd05o7S4JxnyJC3R2jFJC0JPEXgSJQXOgXFVZgW3r83s1ZV+55XNRgrd8KHmjoONEOKXhRCbQogXDqz934UQLwshnhdC/DshxKBef0AIkQohvll//D8O3OczQohvCyHOCiH+iRCVDa4QYlEI8RUhxGv154V6XdTHna0f59OH9uob7lsmaclKN+SxtS6+hNc25mxOc4rScm57xtVhwtY059WNKca6dyRD7sU+nicYpwbfq9yelaw+x55impW39Udrem8aPkrcTWbzL4HP3bD2FeBJ59xTwKvA3zpw2+vOuafrj792YP0Xgb8CPFJ/7J3zbwK/5Zx7BPit+nuAHz1w7Bfq+zc03BNCwPNXRnzz0qjOLDxWOgGehD84u80015xealMaxyvrU7Sxdz0K4Mxym37sY7DMMk1eapyDyFcsdHws3NYfrem9afgoccc1G+fc7wohHrhh7T8d+ParwI/f7hxCiKNAzzn31fr7XwV+DPiPwOeBH6oP/RXgt4G/Ua//qnPOAV8VQgyEEEedc9fu9Lk3NBxklBRM0pJZVjJNS2LfY5qXpNqCEAzigFbg0Qnf3No6vzPn48f6+/e/k96cQSvgTz6ygpKC33t1i0BJIl8R+pJHVrtvK0Boem8aPkocZs3mL1MFjT3OCCGeE0L8jhDiT9Zrx4HLB465XK8BHDkQQNaBIwfuc+kW97kOIcQXhBDPCiGe3drauoeX0vBR5vz2nJVuxFMnFiitIysNRWlBQORJnjzeo9CGrNQ455jnBS9eGXF1lPK7r27ye69tve3W1t4W2LcujTjaj/lrP/QwP/bpE0S+wld35va813tzkKb3puHDyqEEGyHE/xnQwP9cL10DTjnnPgX8H4F/JYTo3en56izGve2Bb73fLznnnnHOPbOysnK3d2+4T5ikJbGv6EY+pxZbSCFYbAdM05Jro4yvXxiy3InwlWR9knFhN+XMSoeTCy0u7iRcHqYY6265tXWw1uJJwYtXxvzGc1eYZiWdyKMTenekbDuz3N6v6dw4/6YRDjR82Lhn6bMQ4n8P/Dngh+sggXMuB/L6668LIV4HHgWuACcO3P1EvQawsbc9Vm+3bdbrV4CTt7hPQ8PbsrftdXWUMss0W7OcQezz4EoHgFwbNiYF7dBjuRNwdmvGtXHCd51ZJNMGKSKePDZACIG2jn7kc3WU8thalWHcuLW1V2sx1vHNyyNmuSbNDH9wdptxUtC/Qwl1JYV+a+8N8L44HzQ03Av3FGyEEJ8D/k/ADzrnkgPrK8Cuc84IIR6kKu6fc87tCiEmQojvBv4I+Cngf6zv9pvATwO/UH/+jQPrPyeE+CLwWWDc1Gsa7pS9DMBa2JhkSEA4x9YsZ5aXWAdRIFGZoBP5dEKf731oiaujlJc3pgB84viAblQFlk7oUWjDLH9ze+vGra29ZtFvXRmzPspoBRIhBGbmSEuDl+u7fh2zXLM+zvYD5nIn3G8I3ft8fnve9Ok0fGC542AjhPjXVAX8ZSHEZeDnqdRnIfCVWsH81Vp59gPA3xFClIAF/ppzbrc+1V+nUrbFVDWevTrPLwC/JoT4GeAC8Bfq9S8BfxY4CyTAX3onL7Th/mQvy7i0mxD7Xl2gV5TGIqXg9a0ZaWk5sRDjkAjh8KTkyWN9VntVQ2eh7f75jg1inr88pBP5OOfeYitzsFl0Z5Yzy0oyLVlsB0SBQhvH7A6DzXWBcpwihWBeaCapZp6XxIHaD4IHs6sPirloQ8NB7kaN9pM3Wf4Xtzj23wL/9ha3PQs8eZP1HeCHb7LugJ+90+fZcP9ys4vsfpZxecxSO6Adesxzzc4855MnBpxebPGda1Py0jJoeyS55dJwzKnFFhvTjHbgsTnNOT6IObPcxlOSEwsterG/v7W11u/UjztifZyx0onYnReMk5LQEwgh2JrkHB1EFMZSGseodgy4HdcFysAj8j2y0jCcF0jhXbeVt5ddvV/mog0Nb0djV9PwoeTGwLLYDji3NbvuIvt7r23xxvacbuix1K6Cw8sbU450QpbaAfNc88ZuwtMnB5zdnDJKS0JVaWaujFJOLbYQVFtnO7OCF6+OOL7Q4qkTA546MWDQChglBb//2hajpGSSlnz76phWoHh4tUsUSKaZJlKC0licAwEoyR0FgD3p8yzX9OoMJvQkrdDDAjvz/C3Z1cHeHGi22Bo+ODR2NQ0fOm7WWf+VlzawlusaIMdpSehLHNCLA0ZpgbCOUVrSiwMcMIh9PCn44Y+t8eSxPsY5umF13+1ZwTCpGjEnWcnJxTZZYfnOtQm//vXLXNiZ8/zlEZeHCbk27CYloZTM05JLuwmelPhSMckLlBT7H63wzpoz96TPndAj11WNKNeWlU7I6cU2gzhgmBTXNYbuKe0OEvvqrhtSGxoOmyazafhQMUoKvvzCOqO0YKkdcmwQ0418jLVcHlaqrVlu6ISKSaKJA8ljaz2ujlJCzyNUDoSgH/scG8RcGSbsJgWPH+3z2JrP5jRnY5LiScEg9plmmlfWJzhcvQ1XqdgmacFXXlonKw2LrYCNSUHkKY4utHh9c8prG1MeX+tiTUpeCNqBpBv5hL6i5as7as5cbAd85aV1ZplmkhYsdyKiQLHabSMlNx2D0MzFafig0gSbhg8NexnNtXGGtZZr44zvXBvz2TPLRJ7kxatTPnVqgV7kkWvLtWnKyToY7dU2ZlmJlFWh/T+/vME0KymtQwjwpeTV9SlKwcmFFtuznIu7KZvTHCEcnUAhhCApDb04YGuWszMvWGyFpKWhHSgCIYm8qnEz9j3WBi3OrHRISksn9OiGHkqKtw0Ao6Tg3NaMM0sdduc5uba8eG3CWi+kHXh890NLN92CO7Pc5rmLQwC0sZzfnrM7L3j61OCO6kQNDe8WTbBpeF+4sfelE3kcq4vwt7ognt+eYy1VgNCOwhjGacn/51tXObEQE/qSN3uBHYv1FltSaGJfsdAKePHKmI1pxuYkw5OCyFNYHF9NChbbIbO8JPYkSkle25gihKTlV2aaV8Ypq72Ib18e8eiRLostH20M46xEG8vFYcY8M7yxM+P4QoynJI8uxqz1Yn7n1Q2+c220v+21Nc34/keubzw+WIdaH2csd0JWexFxoJjlmiPdkHbocWqpzbmtGf3Yf8vPaq835/nLI755cchiJ+QzpxfwlGyEAg3vK02waXjPuZWkN/IUo6S45QVxkpbsznMiX/Li1SEIQTdQJEXJi1dLfuxTx7B1faUTenzXA4uM0pLAkwyTAt8TCCnYnhZ4StEKqiFm2jiOL8Qc70f0wqqW8uz5XQrr6IcCX/nsJiXawO6sYJ4btLF85vQircDj4m7CxZ05k6ySJGvjGM0rCbSSgq1pTqZtvd3nmGaar18Y0ov964QGeyoyTwpeujam1IbHjlZ+bJGvCD3FJNNvW/QftAK6kc9nH1y+bjvtdvdpaHi3aYJNw3vOzSS9u/Ocb16q+lduNZa5F/s8d3HIpZ2UpU5AqR3TwqCAlV5AUhg+efLNMcpJoTk2iPnUqTe77rWxrPQClJBsTDIujxICTxL4kuVuBFSKr3FS0os8pplGCDg2aDHPS7LSMIh9CmP52hs7/LmnjiMFlMayeXlELwroxh6TXDPNKrPP7XnBUtvn+EKLfy8Fnqy22i7uJhjr9l0C9lwHXt2Y0Q58ynro2iQt+NjRPrm2dMKq+P92NZ/GxLPhg0YTbBrec26U9M5zzbVxhjaOB5batxzLfGa5zW98syApNSvdEGOhFRp6kU+gJLvzYn/L7GbNln9wdptX16cMk5yktPQiH+dgnJTMMkNpLA+vdulEPue3E6yz7MxLpnlJL/ZYboe0g4BTSx2mWUngCVZ7EZeGKUvtgOMLbQaRx7FBi6TUXN5JuDpKQMDDq4tsz3KkqPpujHVow74qbe9n8urGlMhXHBvEvLEzJy0M3djn/PacEwstTi91Adia5mzPcv6X56/edBtyb4SCNtAJq/MpKRqhQMP7RhNsGt5zbpT0bk1zJIJBS1EYy1I73L8IH9zyGbQCnjjaZ3taMJxXxe7lTghAJ/J5+tRgf8vsoI/Y7766yR++vsMLV0cM04L1cY62lrTQ4CArLUpaLu4m9GOfcarpRx6ldXRCj1mumWeGrEg5NmhxYXvOTpIzaPlcG6V0QsUb24ZB5DErKneA0jhaoWRnbugG1etMCoMQAuccnpJ0wjdVaXs/k1lu6EWVfPtoP2SUVE4Bl3YTTi226IQem5OMF66MeGC5wxtbM4ZJQVIaHljqcHmY8NSJwf4IhV5UOSB86/KIEwsxf/KRxqC24f2hCTYN7zl7iqmFVsCFnRk7s5TQV3SjkKw0nF5q33LL55EjXUJP8fzlEcaCEjAvDPNcA4PrBAZ7dZAro5RJUjJONMaB70mctkxTjXXgKUHkC/LS8NrmlH4rIPQ8ytIyyTVOgLEOXwpGSY62hlI7fCn58kvXODFocWFnxjzT5MawPa0yrNVuyGLLZ7ET8vL6hHbgYWzV3NkJq0xkT5W29zPxFGSlQQiQQvA9Dy2jpODJ4326UdWYuj3LefL4gCujhM1ZTjesnu8oKTDOsT3NeeJYn8V2WAkwck039G4qKGhoeK9ogk3De85BN+NMx2zPqyL+Xu9LN/JJCn3dls9B9drWLOepEwNGSc6LVyeEnuKHHlsh9BTPXRzy4EqH3XnBNy4OCT3JJNVc3E3ohh7OObaNIy8d2lbatdgTWCuqBlArKAvDVmEotEUIaHmSeWGY5hpPSTypWBn4SCkYJwVZYRi0fLZneTUfJ8nRxrEtYKHl8/haj61ZxoXtOdpUMuur4zkXvjlHCsenTy8C8OBKByUF37w4YrEd8OiRLkoKkkLv153Ob8/59uURUsCrG1P6cYBxjuE8Z5xpHlhscXaW8ZnTiwgh9iXfzrmmXtPwvtIEm4b3hSrgVEHn+x5e3ldixb7an9tysN6yd/vJhRaRpzi/MyMrDU+dHLDaiZhkJdfGGYU2vHBlzHedWSIvDduTnHPbMy4NE4yzlNohhMMBe77N9WAMfCkpjGY4tkjh8JViuRNSGkfsCyYmR1uHFNAOPK6MU0pjGM01nzzV58ljPcaJ5tzOjOMLEcvtkIVWyAtXxvypx1c5Pmjx7+KqRqWNq6aEhh6b44zFVrivxHvqxGC/jtMK5VvGChzpRcxzzdY0B+eYFRaBoB95aGOYJCVb05zVXrT/824aOxveb5pg0/C+c6u5LXtbPjf6fa32IjqRxwtXx5xZau8X1XuRz2sbKTvznMfWeuzMi0qp1gn5RjEkLy1KVdmMJ6GozZxL4wg8yTgtURJ8KYgCj1mm2ZimRJ5CSoGx1Iq0gBeujhkmJf2WQilJUVo2y4Je5NGPAzwluTxK2Z2XjLOSCztzIl8yy0o8JYl8n3lh2JhkXBtlDNOS73lwua5TLezXqvYyuucu7hIoxYMrHY4vtHhlfcJKN+KVayNCz8Pg6EcBO/OcTuTze2e3+JMPr7DSDd8ilmhoeD9ogk3DoXEv1vZ7mc7NuJWMV+B48eqYUVKibbW2Mc3JS8Ov/uEbxIHEkxIlYLUbcXl3jtbgK4GUgkDWDaDOkeSmqt84iCMJCDwl0caSYShLR+RLnBBcGibkZSVDzkvLoOUzTEqiQPHa5pR5bsmHFussncBDSSi1od8KKIzFAb9/dgtPCjqhwpOCN7YTcFuM0pJf/9pFQt/jiWNdjHWUxvH65oxu5DPJSh5e7ZJrx9VRwvntlNV+yEo7YJgVhFLyzOkuaWmq7E8bjg3i64J3Q8P7QRNsGg6Fd9PavldLf89uThkmJQstn4dXu6z1Y77ywjrL3Yh2qNid55zdmNGPPQpjWeuFTDLN+jSn5Uv6kU9uHMY5IikJfUgLTaEdpnZkjkOF5wTTtMRTAgBrwAowzhJIRVJajIVpZuhEAuEsG5OMUlduAgvtkEAJshJGWRUoB62QtX6IdZAUBpzDCYFDkJSWVmT5o/O7tHzFZx5YpNSGf//NK4RKcXqpxTyrlGqbk4xz23M6gSLwJEf6Edo61qcFZ5ZanFhoM041D660ObnY2jfpbGh4v2mCTcOh8G5a23tS8F9e3mAQ+7QDybmtKc++scszDyzywEqHTBuSwpJry2o3YJpXPmUCQTvw2BUFrUixQkxeGoxzjOY5qXa1ygywVe3GOdA4BFDoqphvncNaMAKiWhLtsAiq265Ncgpt8aSgHfpY5yhMpTiLAsk8r+olK90I5xzWOkJPkBSGxAmQcG5rhjWO5dUukV81us6yIZeTlEwbVrshSWEZpwXWWj52tM/WNOf0YotJVrKTFMwLg3OOcVpwbLB8S0VfM1yt4f2gGTHQcCi8m9b2L12d8OhalyhQXNhNUVLSDhRfeWmD5y7usjUtONIL6e+5KgeKU0ttpoVGSujGCuEEpxdjVnsh2lgSbTDWoKRACGhHilYgUEKQ5AZdS6I7kUfgCQRQWpgWJVpXmY1xlSRaIAiVRCpBK1QstkIiT5JpS1ZajAFPVOf2pCT0JEJIAr/aQutHPllhKbTj9a0p37i4y8YkZZwVaGvxlcRXCk+CMZZJWjkZHB+06EQ+R/oxoZIkhWGYFjy21qUb+TcVBdxsPMNzF4eMGqVaw7vMXQUbIcQvCyE2hRAvHFhbFEJ8RQjxWv15oV4XQoh/IoQ4K4R4Xgjx6QP3+en6+NeEED99YP0zQohv1/f5J6KeNX2rx2j44LDXlHiQw1JAbU0zjnQjIl9xYiHGWEhLi5DQDn22pxl/fH6Hs9tTxmnJsUHE0ycX+f6HVjjSjZBCcnKxxVInRFvHMKmMPLEOTwoKA1lu0Bbysgo01oK1jrQOAhrQForSkhv2j8lLS641SkE/8ulFPmmpSUtbDTYrDEi4Ok742vkdcm3ItWU4L5jnGodjXhTMCk1qqp/ftVHGNy4OwQl8KZFS4ElB6Ck8Vb2WVuhxfCGmMAZrHScX2hzpBigpeWils6/oO7Pcvu5neTAD3Zv7cyezdRoa7pW7zWz+JfC5G9b+JvBbzrlHgN+qvwf4UeCR+uMLwC9CFTiAnwc+C3wX8PMHgscvAn/lwP0+9zaP0fAB4cxye/8C55y75cXunbDSjZhmJUlhmWWaTBuEhCOdEOsc06xEiCqT0hYu7iZcHlY9LL4niQPFGztzXtuccm5zSuhJeqGHVIpZbjGuDh6uqueEqpJF5wYKbSntm88lN2Cp6juG6n6lAUdlQXN1mLI9zUnKAgTEgaIdKC4NE17fnmGdozSW3BjmmWZjknF+e04v8mn7HmlhmGQFm5OUcVLQb/l0giqQKyFYaAWcXuqwMUp5ZX3C1jTn/PaMWVFyernDx4520NZdN1DtIM1wtYb3i7sKNs653wV2b1j+PPAr9de/AvzYgfVfdRVfBQZCiKPAnwG+4pzbdc4Nga8An6tv6znnvuqcc8Cv3nCumz1GwweEPfnynl3MrS52d8Pelo8Q8J31KUmhmaYl88JgjOPMahdfCmJfcW2cMUw0p5diWr7i2Td2uTrO0MbiCUHLV4zmJdPMoI0DKSlN1TOjAF9C6CuUhKJO0BxVMLkZjuuP2U00m5OcXJs6q5GESvLk8QGPrHZoBx5x4O07FgRKoW1VP4p8Rdv3OL3cRklJYRyh7xEHPp3QIy011lmiQHFmpc3xhYgoVJzbmnFpmGKdY7kT4pzjhz+2xg8+tnrLn/27mYE2NNyOwxAIHHHOXau/XgeO1F8fBy4dOO5yvXa79cs3Wb/dY1yHEOILVFkUp06deievpeEeuJ18+U7ZK16f3Zzy4tUxgzjk+CDikycGvHBlxDAr8YTggeUWbd/jmnGVDYyBpVZQbT15Cm0dLV/Ra/nMc8PRgeLiMEFbwygzxJ6HcRYlIHeV4syTDtybmcst4swtKR3Y0hH50A0VUaC4tDPjwm7KNCtRUmBtNSlUGwtSYB3khSXJ5kyKEingSC9ipRNwdZxxcrGFMZbAk2hjONJrkxSWhTjk0SM98npezySpsrv/8K0rfO/DK7cs+h8crnYzw9KGhneLQ1WjOeecEOJu/0cP7TGcc78E/BLAM888864+j4bD5+Ccm+9cm4CFWa6ZZBopBD/8sTVyYzi3OeO1zRnDZMrOLOfKqHJdlnO4NDSM0wLjLK9vzXnqRJ/IV8yLyghTm6rQb63GHNgeE0Bu3twyeyfBBkBJsE6wOy8YZQXCVS7PuErdZp3DOci1pdCWOBAIHKWttrh6kWJ3XjDLS54+UU0dvbCbcKLVYrUb8er6hI8fW+Dbl8ckpSbJLYXWXMg1f6K3wHBe7hf9b5bdvF0D7QeRRj330eAwgs2GEOKoc+5avRW2Wa9fAU4eOO5EvXYF+KEb1n+7Xj9xk+Nv9xgNHyGevzziyjDljZ05G+Ock0stPFF5m631I3bnOa3QY60f84fnd9kcZbQjSaAEW9OUS8OkHjImcRZmecml3YRjg5hvXBzWI6Drekvti7YXUAzgDgQf+9and0doC74HpQZtQAiHlNV6HXNqabTBk5XIwJiqYdQYR1Za+i1FoSuZs7aWo/2YJ472yUrDuZ0E5wQvXx2TaUscKMrSUDrHi1fGCClACFZ7AUoKfuDR1bc8x3vJQN/rC/+72b/V8N5yGMHmN4GfBn6h/vwbB9Z/TgjxRSoxwLgOFl8G/q8HRAF/GvhbzrldIcRECPHdwB8BPwX8j2/zGA0fES7szPnSt68RepKNac5wlvPa5pTAqwaNfeb0ACEkUsIr6xOG84wwkNS2ZkxTgxMOT0q0NQjA82TtjFyQlxZrHJI6sAC+qLa+9ninAeYgFkiL6qSOqm9H2sqLTbo9H7YqEClZOTyHPoSewo8kSW5I8mq2zktXx5xeatOJPJLCkGtNy1N8Z32KtpXaLSsNuTE4C+vTjDPLbUJPcmWYMpyX+5NA4d4Dxftx4X83+7ca3lvuKtgIIf41VVayLIS4TKUq+wXg14QQPwNcAP5CffiXgD8LnAUS4C8B1EHl7wJfq4/7O865PdHBX6dSvMXAf6w/uM1jNHwEGCUFX3lpnU6gcAh2ZwXbs5JSG7S2yKgSCBhTvZO/spuRG4OvYHeWA4LSVo2Xha78zQDagUI7h6O6UBeuyjrQVcAx9YVf1M/jsPZdbwxae9+bG9aFgG5LYjSkpaYXt9DWYowjUIpZXnJ5mLI+qZwRVrohjx/t8XuvbdOOfGa5Ji8tWWEIPIV0gkEcEnqKUlusdfsX5cMIFO/Hhb+ZOPrR4a6CjXPuJ29x0w/f5FgH/OwtzvPLwC/fZP1Z4MmbrO/c7DEaPhqc355jLJxe7vDsG7t4UoCz+HXDZeBVyjEhBJdHKa3Aw+aWYVIyzasBaM7WF3NXZQ2egKS0WO1IZSUEwFbKMV0/rrvh880Ia8POewlEN9Z/9mxxktzgajdQZ2E3yfCkJDcWUadAk7RAKUGuDVlpSEuDFJW/W+hVsmohIMlKJsZybZzSChW5rgLQNy4OObPcPpRA8X5c+PfUc3vPFxr13IeVxq6m4Y44rL36g3Np9sYZr08yIq/qru/Hlall6Cu0sax2I04utvCV5PlLQ4bzgpnSpLkh1RptqiBz4wVdO9C1blnaKsDsbbndDRaQ4s0s6J1w410dVcOoJ2sHaiEocaS5QUqDEhIhRLXNph1tTyIdTFJNVloW2z6jRFNYSzvwEaLaMtTG8sbOnCujhKeOL3C83yL0JM9dHJIUhpMLreuex90Givfjwt+o5z46NMGm4W05rL36g2qzjUmGpLLsNxYmeUnoV6Oi+5HHcC4r3zJPMpwXWBzDTBOqKssprSU50Id4u1hwL7WY8jAKOTdhVlgkVd3GDyWBVxl0Bp5ECoUUjrwwhEpirUMGEiGqptNpphECcI7zW3MCX1FaiAMPVc/q+cbFXVqB4lOnFmgFHpuTnK1pzjApmOWGTqhYaAUsdu789/d+XPg/jOq5hpvTBJuGt+Ww9ur3znNpNyH2PSJfsTvP2JpmjFNN7CsGsc/WLCc3Bq0N18Yp1jpKU/mCSSWYzUsK+y5FgfeAvRqRpdo+SwtDHEiEqFRyCEduNKUBT4HvKWJfMsk1SlicVbQjj6K0lM4R4DgxiMm1ZWua0Y09QiHYnuX8zqubPHqkgwNeuJLSj326kc80K3l9c8onTg6YvHJn2er7deE/jP6thvefJtg0vC2HtVe/d55ZbuhFHvNcc21cjVD+xPEer23O2JylfPrUAi+vTyrDSW3JjCErNSvdgElm8H2B1YL8Xva23kfcDV/nBnRa1ZUKZ9HGVoEI0NZhrUYCaWHpt33asUfbq3qHBrVSzViLkoKjvRiLRSIprSXyFc++MSIrNU+eGCCkYJprhBQEftXTc2LQuuNstbnwN7xTmmDT8LYc1l793nk6oarfhedIIPQl57bnjJKSrVnGC1cmABzpx0RKEQeSr10YkhaVwiqqPdAE5tAUZO8Fe95Qgsoixxz43lDLok0VZJQAWVvnlAZKa/EUzAtNaSwzT6EU9FoBic7QxqEkRKFiY1xyfODR8j3WRym5sSx3wkrxZj0eW+vx2uaUXBu+dWnMxjgjDjx6kXfL3pyGhnulCTYNb8vd7NXfTkiwd56FVsCF3Tk7sxQnICkq2/zFVsD2JGM7KQi8qgu/MJYjvQhtLLm1dCMf5ySTTH+oAg1cXzsSB752VP+ISlXBRdS3O1cFKAEkuaUVKbCO1FRza0Ln0fYF3/9QNU56nJZ0Q5/VXkDgKTyvEg6cWGwT+ZJxqon8aproq+tTlBDkhSYvfdKyIFDwzYvD63pzGhoOiybYNLwtd7pXf1BI4EnBi1fG/N6rWxwd1K7NuaEVKtb6MUd6Eduzgp15BgIeWO7w8rUxqXYU2lEYC9YicFzcmXOkF2OcZNDyWR/nFKW9Ljv4MOF4U369h4b9rAbAukoBt9eEiqvGHyhRycGFEDgsmbZ8+vQCH1vr8dK1CcOkZCkKWOkEbM8KnIOlboASglmucQ4u7syJfMW1ccqxfkwn8im0ZXeuObYQNQ2TDe8KTbBpuCPuZK9+TwBgrOPVjSmRX/WAfOmFaxzrRTy+1sNYx2vrU37s0yf4voeX+dX/+ga78wJPwsWdhEJrAiVIcouQDl9I0rpx8dRKi5V2gDWVlcs4+zCGmltzsP9n75UdzIZKDVo4fA9kbeYpAticFBTG8MyZBZQUfOvymDe2E0ZpiXWOSzszHl7r8r0PrjAvNUlpePxIl3FaEPoKh9uf8Pl9Dy814wYa3hWaYNNwaOwJAJ67NGRzmmOs4/XNCcZWxpNfvzjk9GKH0JN89fUd/syTa3Qjj+9cmzBOy3oOjsVTEqUEuDrDcbA9z0hKzaXQpzSGXH+0As2dYADlIPKqHqTQUxgLrVAiCsfpxQ6Xh3NK7Zjl1VZnWhimUgJTFlohnzjR50efPEroKRCC9UnGPK982h5b6+IpSStsBvg2HD7NX1XDodGLfbamOa+uT5AC2qFXbZVNM0pT2fcb67g6Svn917b51f96ntI4ji/EzHLNrNRkhSXNNYVxZLrqc1GyeidfGsPWNGV7WpDduA91nyAFZIUhN4600HgKvnVxxCzXXB4lnN2cszPPsc4ipaQXBTgsW/OS17emKCl46sSArWlGrg3jeYGxlnbocXzQOrSBdw0NN9JkNg2HxpnlNl+/MCQKPHCiGr2MoBVUPl792Of87oxXrk2wxrHU8VnuCDYmGVlpKMqqOq5rexhLrdqqA05pHLk+HMPMDyKS652ob4qoRAShNBQIytJyYXte29lYZnmJs5a8tFX/ji+RorIqKI3l3Nac731oGUfVL7Xcibg4nLM+ThH1nKBvXRo1Vv4Nh04TbBreEbdSnR3tR/gSXt2Y0Y8Dzqy0uTxM2J5kAFytL2qduMp6Lg8Tcm3JysrLS2DfnP3Cm8XxonZs/qgGGngzyNzudQZKYKVDSYkSEs9TGGN5YydhJykxxuJJwTgrcdaRAtpUPnPztOTaKOH89pzVbsRSO6TUjqODmO1pxvOXR5zfnvPksS4nFtqMkqKx8m84NJpg03DX3M6+5tggZrkT8siRHldHKYWx+EpwbZwxyUp8JTm9GGEs+KoabaykICsdkScpSot1Dk9V0zP3fM+shY962dpx80Cz5/vmCXDOEQWStHB02hLnKuueSVbSbwXMS4NDMElLhIM4cgRKEPkepRPk2u7X1vZEHMY6XtucAY7VTsjGJMc6wanFFs9fHtGN/GZwWcM90wSbhrvmZvY1s0zz5RfW6UQe57ZmhL4kkIrVXoAn4figxSwrGaYlWWk5Vg9Dm2UahyM3JUoI4kCijaOwlrDeMtIftoaae+BmGc1eoOnGitJYrIHIl1jj2MwyrAUpq9k2vlK1MadAIVBInHUsd0I6kUfoqf3m2j0nh1fWp2zOciIl2ZnnBJ7i5KLi8nDOONV89sGlZnBZwz3TCAQa7ppJWhL7av/7aVZyYWfGKC0ZxD7aOi7tpiSFpuV7nFlu04s9itpeZndW8OK1CddGCZmphoUFUpKVlmlmKIxFSkHkK3xV1W3ud5QAox0tz0M7S5IbduYlWluss4RKMUtLjHOcWWzRDwMCJRm0A1Z6MWEgWe2GHB1EnFluk9Tigt15zutbU6Rw9FseuXb7E0Jf35qz2A5oBR5CiP03GOe35+/3j6PhQ0gTbBrumr13xntcHaVk2jLLSv7jC9c4uzljfZTy8voE6xzLnYiHV7ucWW6xMcnoRApfOq6Mc/JSk5WmbjisLFesg7J0lNbhKXFdt/39ipKC3FhSrQmlh1dH4MJU/TeWyi3AWbg2zWgFkiio5t0oAfNMM2gH++4Anzq1wGI74I/O7TBNS3qBT1paSmNY6YS8sTOnMPYtyrTYV00fTsM7ogk2DXfN3jvjpKgCxOVRyqWdhNCTXB1l7M5yJlnJxjjld1/d5j+/vM76OMMBDyy1sQ6GicZXEk8qtHX4ShKHPkoppBAEHpTakmv3lm77+xEhBc5BoBRhoJBCokQVZDRQGFsJArSu6mO5IfI8EI7COjLt0Mby1InB/jljX/Hx430+eXJASVUPWmoHlMYyzzU/+Ogqnrr+EtEMLmt4p9xzzUYI8Rjwbw4sPQj8X4AB8FeArXr9bzvnvlTf528BP0NV//1vnXNfrtc/B/xjqp2T/8k59wv1+hngi8AS8HXgLzrnmrmw7xM32tcU2nJysUWuHQ6HUoJZogl9xUK7Gkn80rUJ7cDjsSNdppnm0SNdXl6fkBaV4aRSAmMcxlZmm3D9wDJPVBnPR1mNdjvm9UTPnUQTqWqeTV7/MARgjWM3KTC1oqIdSIg8Cu3oRh5L7ZDhgQFAe3W3EwstSuM4tdTh6iihMJYHltucWmzx1IkBv//aFqOkpDQWX1V2Qd//yMp7/wNo+NBzz5mNc+4V59zTzrmngc8ACfDv6pv/0d5tBwLNE8BPAB8HPgf8MyGEEkIo4J8CPwo8AfxkfSzAP6jP9TAwpApUDe8jewHnBx9b5ePHekSBYpyVRPU2i0MQehKcQ1vHIPZxwIWdOcY64qDKYKx9c6ZLXmqscVWBW0mkoPIBoxIJ3K+B5kYyA+mBdM9R9SCZ+gfkqWrbLckNeWkQQDtU+Ers11v26m7HBjFZaVBS8NBKh5VuyPFBzFMnBozTkvPbc17dmHJlmJIUHy6X7YYPFoe9jfbDwOvOuQu3OebzwBedc7lz7jxwFviu+uOsc+5cnbV8Efi8EEIAfwr49fr+vwL82CE/74Z74Ngg5vRim26okEjiwGfQ8vGVxFEp0RZa1fjib10ec2k34fXNKbNck2lH6Eu0BiWrAWJ7EmCAuu2mgUokcCuMqxpfo0Cg6gidlQZtqy2xeabxpODqKAWqIP78lRGvbszqhlnD9rxgEFdvIgC+8tIG3dDjkSNdAk9yYSfh2jjl+cuj9+DVNnzUOOxg8xPAvz7w/c8JIZ4XQvyyEGLPj/44cOnAMZfrtVutLwEj55y+Yf0tCCG+IIR4Vgjx7NbW1s0OaXgXOLPcRkp4+uQCTxzr0o89Or7koZU2nhRMs5LnLo2YZSVr/YhQSTZnBUmuEcLh+4p2KIk8CaK21heC2Fd8SOejHTqSKkDcCgdoA1lRjY4udTV+rTCOeWHotwKOL7SYZZpRUjBJS2ZZSaDAk5J5YTg2iPgzT64xaAWc355jbOVTd3EnQYpqC204K/nmxRGjuxyc19BwaMFGCBEA/xvg/1Uv/SLwEPA0cA34h4f1WLfCOfdLzrlnnHPPrKw0+8rvBnsNnb/zyibPXRwySor9LTXfEzigF3lk1rE+zhgmOTuzgrw0zDJNWmrCUNGLFEIIIk+R5QYnBL6naAceTsA4M8wy86EcIfBuoGRVt7oVDvazwj0HhnlhSAuDNY7z2zOuDBOkrOo1K92Ip04sEHiKUVoymhdc2Jlzfnu+H4wWWwFXximhpwg8SaAkqTYstoNG/txw1xxmU+ePAt9wzm0A7H0GEEL8c+A/1N9eAU4euN+Jeo1brO8AAyGEV2c3B49veA+5nXMAVCabf+KBRX7gkRXOb8/5zW9dYaEdMMtLZpnha2/sIgX79ZykLNC6mkRpS5hl5rq6TFOjeZPSvvXn4XO9q4Kos0BjKzGFL6utt7huvs20ZWOc0Qo8Ti60EEJwbADTTDNYapNrQ6Etz10ckpaGYVLy8tUJi52A1W6EMQ4lqky2kT833C2HGWx+kgNbaEKIo865a/W3fx54of76N4F/JYT4H4BjwCPAH1PVgR+plWdXqLbk/nfOOSeE+C/Aj1PVcX4a+I1DfN4Nd8jNnAP21ve+31srjeVoL6QVemxOMl5en5CVFiUgN5Y5kOgqtfaoxyK/9y/pQ8XBYOMJ8JTA1dLwPZsbA3iuCjJOVL03/biypBmlmnPbc15en9EOJavdmHbo0Y99wNGN/H03iNc3Z/hKcHq5xe6s4NWNKUd7ET9Qy6GbMQQNd8uhBBshRBv4EeCvHlj+vwkhnqa6hryxd5tz7kUhxK8BL1G1CPysc87U5/k54MtU0udfds69WJ/rbwBfFEL8PeA54F8cxvNuuDv2PLUOoo3lhatj0lxzpBdxfKFFN/KZ5YbFdsh31idsT3O0dXhSMMsN9bQB4M0+kYa3x1dvjozWtsokfQ+ErepcfjU1GlkHIiUk2liujjKO9CNW2wEvXRvTiXwGcZtpWvDtS0M+c2aRTuhxeqlq4Nyd5wSe5BPHB3B5zNVRhldbCSkp9keC324EeEPDjRxKsHHOzakK+QfX/uJtjv/7wN+/yfqXgC/dZP0clVqt4X1kzzlgL3uZZiXfvjKiE/l4wuOV9QnfvDTksbUeUlDJoLNqWmQgJaOsAPfWDKbZLrszCgOBhMCrepKolOV4teuCpxTaGIwFbR1SGCJPoiLBNNM8f3mEwyGogv6RXkToS85uzfnfPnOSblQ1a+4mVb0GIPQEnz2zxDgp2JjmnN+Z8SNPrAHccku1CTgNN6Mx4my4Y84st3nu4hCous/Pbc1wCFY7Ed++MubqOMM4R3JxxGovRAJHBy2KYsa6ztG2etfd7Je9c0oLtqzGL/i1Rtw58D2BsYbSgS8g8CSlseTaUhiHsZphUtKLPCIleXVesDHJOLUYMS+qPhvnHGlpUBIW2yGvbU73J656UvDxYz0eW+uxOy/YnRe33FJ9u/HhDfcnzcZrw00VZjdjT3UWeLJyDjCGp473WZ+kTNKCY4OYhThglmuSwvDQaoenjvdYW4hxwuGrJs7cK0pUPTUWKCxoXQWbyPewDjq+RzvyEUoQhR6d2CctNCDohQptKk+6wJOEnuD1rTmhUuTaMEwKAk/yI0+skZaaF6+MEFRWQvNCM8s02lQjCm40Y4XGN63h9jSZzX3O7RRmN9sOqQJOtd6LfQpteW1zSpIbspklKTVKSGZZyR+e2+XzTx+jHfic3ZhybZSgTVWQE+L+Gh1wWGhX1WwUlRw69BVpbhjPS6yDbgiBkoRS0g4UvpIMk5ITiy3SrOT87hxZ/34EjtJYTi+36Ub+vqoQKqeHhXZIUjr6keSR1R6qlk1//Hgf4Lot1b3vG9+0hlvRBJv7nNspzG62HXKwKCxE5fh8ZZQSKckkM+Ta1E2ZAThNWfulHenFDJOCKBBkpSYt3ow0ispipTRN/eZO2MsnqhHSDqWqhk5PQVZaenE1FkAbyyy3nBhEGGPYmhcIJ1BeNclTSMHRXoQn5VsyEufgBx5Z4dWNGZGvCD3JcJ5zdmvGQjugG3lM0pKVbkTsK9LS7AsHGhpuRrONdp9zN9shF3bm/PrXL/PH53fYnGSU2rE1zVntRGzOSpR0dEJFFCjmheHEYosro4RhUvLQSocHlzs4LJl2yAN/eQ5QCg68SW7+MG+D3vuwUJYWXU80dXVnp6cE7UDu13FCX7HSCZFSVM2hDrqhTytQOOfYmKZvyUh6sY+nJI+tdfGVYGOa8cZOwpnlNicXWoSewsF122+NOKDhdjSZzX3OjQozuPl2yCgp+MpL63hCsNiJyLXh4u4c4+BjR7toa3EOro0zpIBZVnJtnPHS1QkPrrTRzhJ4Vc2gGyhKB8oYSl0Fm7Rgf26NoMlw7gRJJRjY85ITgBKC0bxgnhsePdLBV5Ktac76JEdri7aOTGtKHGGgGGUlO7PiLXNr9sQgrcDj0SNdMm04vhDx5PHB/iC11W60H2QaGt6O5g3kfc6Ns2n2vr7x4lN5ZUEvDiqbGd8j8hVaV5nKUycGHOlFDFo+o6Rgd6556dqYc1sz/uvr22yMc04shASeQkpwtprGaakCi6s/N2K1O0NSZTN7PzNZz7YJfEm/FdCPq8zFmMrV4eNHe+TGoI2lFXpIYJqWZIVllJT7NjXw5lZpUhje2J5zaZiQa8snjg/25dHQCAIa7o4ms7lPuFUD3o2zaXqxz2Nrb90OqbyyfHJtieptt9BTeJ6k0AYhHOOkZHeec2F3jkRWHe7AuDZ/bAU+ZxbbvLwxoSgdN7tMNYHmzjiY+Umqeo2vBKWxDJOC00sttHZ0Io9+HDIvLIFS4AuUFPSjgEwbBJZ2oPZtah5c6XBua7ZvabNXi3lktdMMUmu4J5pgcx/wdoqzgwqzW91/fZxxbZwyzTQLbZ9SW7bnBc46HjnSpbSW3Bpe3ZhRGkc7hMhTpNogLVzeTbmym5IZw6xoQsphEYgqq9EaQiUIhGSh7TNoBUwLzfFBTK4Nm5OMdqAQQF5vpxnjmJSWUVbwjYtDuqHH+e0Zpxc7XNpNmOWaTlj9zSgJSVF5PTSCgIZ3QhNs7gPuVnF2kL1AtdwJmRcaYx3fujik3w7pRR79KGCUloSeZKUT0Yt9nLMooei3AmxasDvNKWrJc9D8xR0qhQPpqp9rHHooBJ6S9eC6qrlTSRilOc5Vs24WIp+iNDjnyLRhnhu+9sYOpxZaXB2lpEVlNdSLfHJtuLAz40g/5vseXn7bDLih4VY0//r3ATfzNIt9xfAOZpIcDFRxoPjD17erC4xzrHQjdmYF2hpi38NYRyfyyXJDaS2TtGCWlJTmzeKgbir/h041qROmqeaRI10+ttYlKUouTDLGacnDqx1OLcS8ujGnMJY4UGCrLbe27zGo+6U26uOvjTOODaqaXeR7++Mh3i4Dbmi4HY1A4D5gT3F2kDvdbz8oje5GPnHg0Y8DPCnJCsP6KOWFy2N2ZhmtQHG8FxEGlepsY5KTaYujKmaLA2OeGw4PB/WQOcf6JOVrF4aUxnF6scWxQcS5rRnnthM8JXn4SJvSOiIl8aTgxGJMO/TxhGCUlBwfxGxOMrK9zKesxj50ouZ9acO90fwF3Qfc6Gl2N/vtN0qj00Izz0uS0mIsnFhqVT0Y2ykPrXZphQpfQTvyyLTB2jebEJWsCtuerCS7DYeLc460MAznBVlpWGn7WAdx5HOkH3Fyoc0krWTRx/oxb+zOSUuDnma06trMkV6Mn+T4SjDJNJ1Qsdpts9j5YGY0jfP0h4cms7kPuNHT7G4a8G6URispWB9nFNpyZZSwPcs5NoiIAsGVYcpyJ+DkYpsj3ZBBK2C5ExB6lcQ5M1Wnu20CzbtCqqsBdMbAMMm4OEq5OskoCk1eGC7uJMxzQ8v3SApNN/CrbDVUTLMSbR1Xx2lV12kFfPrUgJOLLaTkLVL4DwJ79cRCWxZawb6irhlZ/cGkyWzuE97pfvtBafSlYcI0K5kWGodACkEsFXHg8fhan8CXpLnm+MDQj6qM6Oo4QyqBbx15PW0yUhBJgXUOa+vZLIf/ku87HG9Kx/MCIl8gcJzbmTNoBax1Bc5JokBRGMtCy6d0kOWGVFvGaUE3brHWj3np2oQHdZtjg/gDKwS4F+FLw3tPE2wa3pZBK+DMcvVO8kg/YnUcIYXAkzBo+ZSm8uf69KkFro1TXrgyZmuWM880vpKUxiGkoy0dQgk6kU8n8Nicpjgr8Jxl3vQGHgoHk8Zca5ZaAYXRaONItGE5UMShx1rk0Q09Ftohr1wb4/uCyJMstQN25wXGOR5abX+g3QHuRfjS8N7TbKM13BF77yID5fHJkwOkEGS6sps/tRCTFZbFdsDrmzPakUfkKTJtcNbR8hX9yGO1FxP7Hm3f45nTCwgEhbYY14gGDhtfgTaOaV4NUAtUZWNzbmeOJyHyJGlpWe2GdOOA04tt4tAjUIqFdkDkSX7v1c0P9JbUvQhfGt57Di3YCCHeEEJ8WwjxTSHEs/XaohDiK0KI1+rPC/W6EEL8EyHEWSHE80KITx84z0/Xx78mhPjpA+ufqc9/tr5vc316D9lTpXVCRS8K+K4zSzx2pFfNRQkUT58asDuv60EnF3jkSAepJEoJAg/CwNsv4Bba8o1aMZVbyHTjHHDoGIiUwlOSeWHIjWWxFbDWD9maFezMCzpR5RxwbZxwdZSy0qnshASCQCkCT3F+e/5+v5JbcqdWSw0fDA47s/lfOeeeds49U3//N4Hfcs49AvxW/T3AjwKP1B9fAH4RquAE/DzwWaox0D+/F6DqY/7Kgft97pCfe8Mt2HMQ+MNz2ySFZpjkKCnoRop24LE7y4Fq3MBiK0BJwSdPLvKnHl1hEPs4J+iHVff65d0UT1VjiY1xqNs/dMNd4IvqA6rttH7Lpx14xIG3by2kLcSepNCWo/2YwliWOyHbs5z1ccbVccIkLdmZ5wRK8AevbX1gi+73InxpeO95t2s2nwd+qP76V4DfBv5Gvf6rzjkHfFUIMRBCHK2P/YpzbhdACPEV4HNCiN8Ges65r9brvwr8GPAf3+Xnf99zo4OAtQ5nHZuTjDe2Z3zseI+PH+3jKcmr4yntwOPs5hTjIPAUndBDCEPgK0ZpyawoafmVW7CvJEJYTKMOOBSsq2Tlg1CilGKlFzFKSlbaPq3QJ/Y95rkm14ZO5PPpU4tMs5LtWc7L6xO2Zxm7SbXl1o18Ti21WOmE+yqvD+KFvGk0/fBwmJmNA/6TEOLrQogv1GtHnHPX6q/XgSP118eBSwfue7leu9365ZusX4cQ4gtCiGeFEM9ubW3d6+tp4M1azWov4vG1Hp3Ix/MUubH8uU8e53seXNkf1rXSifjmpRFLnZBe6LE+zlBSsNqL6EY+kkoxNCs1kVdt8chmM/TQEAKkrCxqVroBi62AQAk8pejHAQ8d6bDY8UlLyyTXTLOSs5tThrOcR1a7KCkpSsv2NEciiH3F8YXWvuLrg7yl1vDB5zAzm+93zl0RQqwCXxFCvHzwRuecE0K8q1vzzrlfAn4J4JlnnmnKAIfAQcVPN/I5NqiaB799eYQU1HNNQo4NYgptODqIWe1GbM1ykKCtIy8MxhPkppp5M89LJklZD996X1/ehx5FJRsPVP3PLKhGQCiJLwXHBhHz3BB6km+8scM8MyR1Uf0/fOsKm9OcbuTRiXw+drTHLNO8tjnl2iThf71wbH+kQKPyarhXDi3YOOeu1J83hRD/jqrmsiGEOOqcu1Zvk23Wh18BTh64+4l67Qpvbrvtrf92vX7iJsc3vMscdBC4Nkr5o/M7jLOS3SSnO1VkUdW0Oc1Kxmlld3JsEDPNSlqhhycFW/OcbujRi3wupwmFdkhRWaw0/Z13j187MAhASIhVNaVTCEngCZ442mOaaZLScKQfMogdf3xxiNGWVqjoKMU802SlZZ4bwFXZaBTQiwJKY5nnJdfGKZO0ZJZrPCk4tdR6v196w4eYQ9lGE0K0hRDdva+BPw28APwmsKco+2ngN+qvfxP4qVqV9t3AuN5u+zLwp4UQC7Uw4E8DX65vmwghvrtWof3UgXM13AF7tZffeWXzLQXf2922p/jZnGR89dw2xlgKbTi50GY30RSlZZJqBDBKSxbbIVdHKZGv6IU+ma7Ge0X1zBQBhL7Ek5Ubcdhso909rhICRJ6gX8uVrQOHxfck3SjgEycGRL7ipSsTfv/sNoEQPHqkQ+QpSu2IA49xVnByMUYbx2he4nBMsoJJVjLLDb/7SiV9DpRkmmvGafmBFAo0fDg4rMzmCPDvajWyB/wr59z/VwjxNeDXhBA/A1wA/kJ9/JeAPwucBRLgLwE453aFEH8X+Fp93N/ZEwsAfx34l0BMJQxoxAF3yO3m2QB3MOtmgS+/sM680Kz1YpwQLLYDOpHPNM3ZmTtOLMQ8cbRHWmqevzwiVJKdWU6uLVluuFymOOcIfQ+0wQpHywvItUZntnEQuEMEUDpo1xZA2jp8T6At5CV4QvLqxhRvW+CcY5pqSmtISsHrW3Oy0iCFIJ9mbM0yTi9bhINzOzOGSUFSGo72IpY7AbPM8vLGlIXYJ/I9Xrgy5sLOnO97eKXxIGu4aw4l2DjnzgGfvMn6DvDDN1l3wM/e4ly/DPzyTdafBZ685yd7H3I7W4+9729n+TFoBaz1I545vYi2DiGg1JZeWI0VePJYj5OLLXJtGNf9OOO0ZJoVaGNZ6gXszDLmRWXMaawlDBSBJ9md3T7QCEDVW25NEe7Az0BAy5dMc4NxEmcqNdrmNGdjkuN50A4UCkEnqpo0r00yPAmtWpUWBBLhqBweXDVMzROC7aRACcfTJxfYnZUYWzVLCuDqMGV3VjBKig+kOq3hg0vjIHAfcHBMwB578+Nvd9tBerHPYjskKw3dyCfThuE8p9SG3XnOH53b4fXNGS3f4wcfW6UdeSzEIacXW8wyTeh7RL7COktuQFiYZ8VbAs2Nu2oO0E2geQuFdiRF5Q4Q+RJDVf+yrhIM5BpGiaG01U8uK6sMRgpJUhqEFJxaaDPLDdbB9zy0QqI1/ZZPP/QotOP3zm6zPk75zvqY9XHGMC0Ig0oo0KjTGu6WJtjcB9zO1uNOLT/OLLeREk4ttunV6qXcWAJP0Ao9PnN6gcI4LuzMAFjuhCAqh+fltk9Yd6RLIYg8SLUh17ylqbMJKrdH1B/aVQPT0twymmt0PQn1xp9faSzGOjwlUAqsrbyBWoHHOCmZZQVLHR8lBSudiLVeXDXshh7GOC7uJozSEikcWWHQxrI1y2/6hqSh4XY0Rpz3AW83z+ZOZt0cdH/2lODhI12mWUnoKVqBxzQrmWUlw3nOMClZqeXQmbZoZ9mcFbQDRVF6eApGiSYOJaEnmBeGopGl3RF7wURQ9dXYOuvbC0KyXt/LBq1zSMBXlUu3AQaxTz/yMdaxO9e0A49xVnJiIWJ9nGFc9fue5iVpoemHAdNc8/jRPkYbktzcsQdZM2+mYY8m2NwHHAwUN5sff7vb3nqeN9d/55VNYl/tS6KTQjPLNJmxRL7CWMs8K7m8m6KtZW4hDhQKSRg4ssLSCRVKSWZ5Sdb03NwRHpXk2VeCVuyhhGNnrtGuCjZ7xqai/rowlhOtmKS09CMPB+zMczqBx8mFymWgFXqsDGImma7qOZ6iFVbzbo4PYq6MMox2GOdQStx2+N5egLk6Srk2zjiz1GalG75FfNJwf9EEm/uE29l6vFPLj17sszXN+aPzOygh6ccBo6RgcydHCcFKN+SRtS6vb81RQlAYQ15WV0JPCFJr6UQRhbbM8mZL5k6QvJnBeJ7kWD9imml251WktgeOawWSbuTTjioT1NOeYLkdsjsrafuKIFA8sNgiLQzf+9Ayr21M6MU+D690aIUevidYiAOSQrPcCQl9gZI+R3vRLQPGQeVjkms8ARd358SB2m8QbebN3J80wabhHXNmuc3XL+ySaUs3FJzbmjMvNYPYp9COK6OMP/eJY8xyzTfeqLbqUmfQxmGsQQrYmeV4niBQlTlks5t2awSV9xkOAk9wtBuTFJpMW/qxIikMSZ0dtkJBP/Z5fK1biQac44GlNpNMs9gJCb2Q0Pc40g2ZFYZWoHj61AKzXHN2Y0boK77rgUXOb8/oxj7f9/AynpIkhb5tZnJQ+TgvLL04INeGq6OUx9b8xongPqYJNg3vmEEr4Gg/ZjQveeHamJ1ZTif0sM7VUmfLf35lncjz8DyBNIJASjwJUiqywlYOxMbVxp1Q6OrduaARCxyk5Vdii3boIyX0Ih+Eo7QOh2OxFbHcE6yPMnwlWOqELLQCji+0maQlTsAgDrg2SulEPmkJG9M5m5OMH/3E0f0eK2Mdn//UcXbnBZO05IljfaDq52mF8rZTO6usZheJpBN5SOHItSX0FJOsylybeTP3L02wuc+51wLusUHVgf4Hr2/Rj3zaoccoLRjNS470Il5dn4IQjBONsxYrBNY62qFPmufMC01uLFKI6wKMApoSToUEjnRjui2f2JNMs5JSw/YkZ6Uf0QocxsLxQcRyy+fapBr5MGh5tANF5EumaclSO+C7Hlji5Y0pw6Qg9hWPHe1xZrkDvNljtTsv7npC5972WaAUUlS9O7PMgNC0fEU79Pbnzdyq1tPw0aYJNvcxt3MWuNOAs7eV1o48ktwwSkuK0jJo+2xMMiaZ5kgvRAGTwuB7Eu0cu3Xm41zl66WkwwKxVxW1dWMpsI8FxlnBmeUO1yYpnpT4kSAIIrLCYq0j04admQ84ktLga0thHIU2+Ery6FqXduCBsHzvQ8scG8S8uj7F3ZA/amN54eqYq6OUWabpRB7HBvHbvgnZ2z57cKXDK+sTIl8waPnMihLtHK3QI/Bunxk1fLRpgs19zO2cBe6mgDtKStJUkxSWlX5Amhs2xhmpNvVwtZLdtMBYh3UWJSWpsVhbXUgDVQUXB6S2qksoVfl/Fbr6urzPg888N1wdJ2S1KWqmLUd7Ea9OZpTGICSc35pSWMNKJ+ShlQ7WOk4ttTk+aLE9yzm11N7/HQN4ngD3ZhvtNCv59pURhXZ8fbhLrh2hEuiTC2/rGLDnDi6E4LG1HldHKdOsREnJf/PpE02AaWiCzf3MwfEBe9xNAXcvMyq144njfTYnObPCUGhLVtrqnXHg0YsV80Ixyyp5bqktUoLvQanB2OtVVMaCFBAGErBIKVFYsvs04AigMLA1yznai9hNNP3IY5yWDGLFVlKNfVZCEHgKbRzOOc6sdKo5Q92QTFf9U/BmP1U/9hFAUmhiX3Fua0ZaWnZmOS1f0e8GzDPD85dHPHViwJdfWGetH910u/WgO3g38nlszScpNIEnm0DTADQOAvc1d+oecDNGScGXX1jn1Y0pSVmCEJxcarHWC6s5NdYSh5LIV1WDoRC0QkWgJKGnCGTVX+N7oGSVzUgBnqo+I2CaWzINRltCX+If+Gu9n/5wXf0xzTUWUbk9twM2pznaOWLfQ0lBK/BYjH20dWxMC6x1bE5znr8y4vz2jPPbM772xi7PXtgh14Y/+cgK3//Iyv5Y5cIYFls+npK0Ix+BoB0qkkLzrUu7jOo3J3uTO2/mDp4UGufc/tdnltvv28+t4YPF/fQ/23AD7/QCsZfRjNKS5XZIJ/SZZZrtWc68MEgpWGoF9KOg9uRylUqt9ukatAPaoUI5MKZSWTlb1W5wVfDRdaqjqGo448xWY48FdMMq8Cje6qX2UcUTlZ/cziyn0BZPgJKCpDQ4a4mUonQWiyBQinaouDLKuDZO2ZrmdW+OoBUoTi92MHu/i7rh9wcfW+VTpxYrGXXkU9a/gNI4jLUU2rLUrrbJbja5c+88e4Er8GTTvNlwHU2wuY95pxeIvVrPUjugMI7FVsA0K3DAg8ttFtsBEljphQzaPklhKEuNEJWVSlZoCmORqs5kZLUee9VIY1MHmlBWt1N3xitROxAIgaMynLxf5NEC8FUleV7uhry6OWWhFVQOzqUlKw3UAfvxtQ6xJ9iZ5ZxcjFlsByy2Q/qtkNj3bmmkeWa5jZKS0JdkpWaWlWTaVHU1X3FsEO8fezNvtOsDVxNoGq6nqdnc59yte8DBXgopYZZXFjWnF9usTzPGacnDKx1WOyHfuDSkG/p0Io9ZIZjkhthXGOsYtEKSvGQ5VoxyTegpfCWYZAWzohIHWFdtr3m+xNqq4XOtHzFNi0pWe4i8n309t3vsqmpVzbBRgLOOeW4oSotzlihQBLnEGouvBIEnKBwsdSMGbcfGOOPyKOX0YotTSx1avmKS6ZvW5gatgB954ghfeWmdTlRtx3lSYIzPJ08s7DsAwK23WxsvtIZb0QSbDyHv9j/0rc5/YWfOV15a59JOQitQtEOf9XHKhd0E4yzd0OextS5rvZivnt/mY0f7DJOC4W5BXhjWuiEOQehJlBJ4yifyFINuiLGOjUnGLDNVjcJVGQ8OstpZwBMQSEE79Bkl5aHIo+uHAKqLuQTeDeOc2wUU74YR2YJq9DOi2mZUQOSB73vMSk2SV4X+K8OUR9a6PHG0x6zQTOtJqUd7IbFfOQpIIekGHq9vzbg0TFloBSy1A7amOYudt/7NnF5q8+OfOXnd73+xHXBua7YvJLiVWethSOkbPro0weZDxrv9D32r8z+40uErL23gicrx+cUrI17dnLLcDim1pRMpVnsRpbFcGSUESnFqoUVeGo4vxOxMCwpj2Z0XtAKFJwT9lo9xjpXYZycp8JWkHShCa5kXVaNiyZsX6cgTbEwz0lpVFdQX6VvFnD0zSsetL/SWN7MHQ7Wd926kOG93Sl9B6CnmucFX0A48Ul2NDvBVtY211o8ZJTkOaPuKzDh2Zjk6dsjaIeCh5TYOx7GFGCUk65MMIQRXhinOwUovYqkd8MKVET/26RP7j/92b2D6sf+2Zq2HJaVv+Ghyz8FGCHES+FWq0dAO+CXn3D8WQvz3wF8BtupD/7Zz7kv1ff4W8DNU/9//rXPuy/X654B/TPVm7n9yzv1CvX4G+CKwBHwd+IvOufvSYOmw/6FvvMhMs/Km5//q6zvMshLnIC1t5Q5caDbGGa3QY6kTEXsek1QT+pLLuwkPr3URAsZJSaYtRa18m6YFgVIstHxaoces0AznBeBY6YVsT3McVYazN6Ml8qkaR0sDVC4E9jaBZg8Bt/VbO5jZ7Nn2i/d4WFvpoOsJ4kCy1A4wzqGNRVtFEDmcg9IYdmc52hpagY/FsdTy0MaxPctphYo//6kThEry8sYUrWG5H2Cs4+Vrk2pbTBs8CUlp+OSJARd25uzOi9qdOeXMUueW7sx3st16r1L6ho82hyEQ0MB/55x7Avhu4GeFEE/Ut/0j59zT9cdeoHkC+Ang48DngH8mhFBCCAX8U+BHgSeAnzxwnn9Qn+thYEgVqO5L7nSy5p2wl8UU2u5LWr95cYQ211+eY19xcbcKSGlRGWhuzwuUqEY7n1poE3oSISxXxxnb04wTiy16kYdzgp15TlFWvSCdwCOOfDxPoDzJdz+4TD9UnF5osxCHCBSeJwkVhAo6oeLUYkwvDJhkJXlZ9ZPEgcK/A4utO9lpOxhYrH1/ajeyVnkpKejFPgvtSlzRiwPW+i16USXCmCSGaVa5aoehqrKcyOfMUhtPSi6NUhyOVzcnXB2lnN2cUWhLL1KcXGrz9MkFHl7psDvP+Wb9u6/cmQUXdxNmub6p2uxOuBcpfcNHn3sONs65a865b9RfT4HvAMdvc5fPA190zuXOufPAWeC76o+zzrlzddbyReDzQggB/Cng1+v7/wrwY/f6vD+sHOY/9MEsaU/SutgOOL9z/UUmLQ0OwXI3xFG9kw6VYpIUDJOS3BgCT6KkZK0XcmqxzcOrHR5f69GPK2POtNQIHBaHdKAQDGKfx492sQg6kaITeSy0PWJfEfqqHgQmKK2rjB1l1aPTjQM6kU95E/O0PSm0pMqKbpwEeiN7pp97Kf6dBKe3O+fd3seToIRgmlUOzsN5gaq3xfpxpebLS0tpqqFoeVm9UZilBk9JAgWdOODquKrJ9EKPy7sJ/+mlda5NElqBZJYZ5kU1KC30FOe25ix2wtqd2dCLAyJf7Xf+X9yZ8wevbb2ln+Z2NL02DbfjUGs2QogHgE8BfwR8H/BzQoifAp6lyn6GVIHoqwfudpk3g9OlG9Y/S7V1NnLO6Zscf+PjfwH4AsCpU6cO4RV98Hi7qZt3w822Pfa8zm4sBp9ZbjFNS9b6MZeHKbO8oLSW4wsx1jjObU4pjENKQeRJPn16kV7k0408jvdjUu3ISsMsL1nqRCy1fZRQfOU7Gyx1ApyF9UlW+X6pSkCgCxClxloLsYdw0ImqOSvTtLxpfcXd8FlJMLfZR9s7xduZfh58qHeiS3i7ulJW2mpejIBpabgyTCmx+EIipcA4h6DuLRIgcYzmOQvtAG0URWnYyDXnt2b0Qp/VbkRpLZeGKS1fIoQkVJLtWeUKXRjLmaU206xkZ5bzxs6cfuRXzaNZiQCO9KL9Bs47qQm+3ZC+hvubQws2QogO8G+B/4NzbiKE+EXg71L9j/5d4B8Cf/mwHu9mOOd+CfglgGeeeeYj2YJxmP/QBy1G9vCU5OkDvTd75z+/PWd3VjBMCnxPcGqxTTfyKbRFCrg2ybDWcma5SyfyeGN7znMXRvieYF5atDF0Qx9PCpQQtHyPwBOESnBysYN1lnlu2JpnxL5COsGg5WGsozSOtLCcWGyR5JpCV5Y45W2CyF7KLkXdGHqLY8QtbrsbvJucY6/epCS3HXktqbbujHC0gyr/sQYKU8/9Ua7aNhXgexIlBEoKAiVIS42vBJ3Q48owoRUqSu1QUmKt48RiTKEd7ahyFrg8SviDszt88kSfz5xeICstF3cTBq2AtNCMk4LtpOCJtR5x4HF8oXXXNcF3Ooiv4aPPoQQbIYRPFWj+Z+fc/xvAObdx4PZ/DvyH+tsrwMkDdz9Rr3GL9R1gIITw6uzm4PH3JYf1D32rLOlm72LPLFdbNycXWzy+1mN9nCEEfOJEn2ujlLV+SD8O+PjxAa+uTxglBTvznMCTLHcCpmlJkmtk3QvSiz2ErFRXl4dzHj3S5eMn+lza8dic5pTWUWqLcZZBHLLQ9lnrRryyMWF9nFf+ady++O+oDDzFDQfuZRO+urm79I2jDiy3r+NoIFJVoNamcjpwDqSqpMu3ErjtSZ73nqyvFEVt8yNLyDHVNqKxaFuNY7YOIiXxPQmi2mY8OgjZmhTkZRWwBDBMC157ecpyJ6QdekxzTeR5PHIk4PRKm9hXvHRtTDf09rPbq+MUr3YlePrUm301TZG/4TC455pNXVP5F8B3nHP/w4H1owcO+/PAC/XXvwn8hBAirFVmjwB/DHwNeEQIcUYIEVCJCH7TOeeA/wL8eH3/nwZ+416fd8PdOQgcPDb0FScWWjy0UtmeFKaayLjajWgHXt2D45GWVff5J04s8MwDy5xcbnNqIebBlS5/+omj9MJg3wbFV4rIU3z2wWUi36MbKE4vtnjm9DIPr3YpjeXZC7uMU81C26ff8pCyusirm3jWeLLKNgxV5hDVgoNQvqlA25NNC6oLvy+qc4UKYvVm1nMnKXJpIPQlR/oRnUjRb/l0fA9PvVkTurFuo2tnhNivTEfTUjPJCtLcUBpHJ/RZ7YV0Qx9jHcI5Sm3IS4OxFiUkSWEYJZpuy8f3JMY5Lu4k+EJirEMJwRvbMyJfcnKxxWdOLxIoxUo3QghBO/SYZCX92OeHHz/C9z20TDfy76iBs6HhbjiMzOb7gL8IfFsI8c167W9TqcmepvpffQP4qwDOuReFEL8GvET1pvBnnXMGQAjxc8CXqf4vf9k592J9vr8BfFEI8feA56iCW8MhcLMs6VY9F3vH7mVEraAq5ntK8Mr6dP+CZJ1grR/XneoSBXQjxZFuCA6ODGLGacFSJ+DybkIvDjDWcnmY8OrGhH4s6YUdJrmhNJatacb2LGeWV89nnJQEnqITVI2Lle/Xm6jaOX/v4u6rSrHlCxhl5X5WFEr2xxy42pOtHSiEEDgHnjIkhduv+9xqu01SBTdtHau9CF9KxllJagymVrftCREO3ifwBN1I1eIAQyos1jgsVXZkPEkgFYuLIVIJxmlB4FUZEw6Q1ZyY7VnO8UFMJ/R48cqIdlB5o51YiPGUJC0lZzdnfPzYgCvjlLVeROwr4kC9ZezAYjtknJVv28B5J38vDQ0Huedg45z7fW7uh/il29zn7wN//ybrX7rZ/Zxz56jUag3vMgebOj0p+NobO/za1y5wfKHFUycGPHViwKAV8OBKh6++vsPWNKMdejy43CbXhrQoUQJG85wjvZC1Xoy2llFSstwJeXi1S1YaNiYZR3oRD692uTZO+c61Md3Y52SnhScFL1ydVBnXvGCSFWSFJlQe3SggK3NybWo7fYvnKdJCU5o3sxzroBNXo4mFEGhjmWuDpwSxc1hbZRaBVx0bKkFpHEoKQl/hSVhoxaxPMzwh2ZnlmFtEGwEEvmS1E7IQ+yy2AnaTnPNbc5SwJLl9iwAh8EAJxzzXBEoR+xIrILUOT8BSO6q2zXCMZjnHFiI6kUc39Ci0ZZoZ5nlJkhukEIzTEucqEUYn8hACHj/aY5yWTNKScaJJ8pJCG7qhx9Y055HVzlvGDkgJP/LEGrvz4o5qgo1rQMOd0jgINFzHnhzaWMe3Lo3ZmGb7rs4vXR0zSUs+cWLAua0ZDyy3+djRHmlp2JxWSrKr45R25NEKPU4vtdieZkih6MU+pxfbSAk/+NgK57fnFNrSCjwmacmTxweAwFeCY4OYea4ZJVWdJ/YUZeBYaPuEnsdYlkzLEqUEQggWY588UIzSHGurLS1fQOxJQikZFbqyx7HVvBchIdUGW9dBlHF4nqIVCdq+h0MwzgrmeUIv8kAoWr6HQJPeosajtaUoDFuzgkHs87G1AbHv89rGBOtK5oXbP9ZR1YqMgEBC6AkssBD5tD2PTBv6LQ9jHN2WhyclHz/WY7kTcWk35fWtKb4nybRmpR1QOmr3Z8diOwQExxZaHF9osTkdIgUY57g8yhi0PFaLiPM7c378Myf2f+c3BpbTS3cmV25cAxrulCbYNFzHnhz61Y0ps1zTDf1KUZZr+nHAKCn56us7PLDcvu4Cs9qNCDzJj36iKtXtba0oKfbHCy92gv0tljPL7IsTpllJ6ClybTi91KUb+XzXmSVe3phyernFxrjg3PaMSVYyzTUOWOoECAEzKfGUZND2iXyFdY7hvCBQEu0sZWkRWIwBTyqEEwjpcK6SE2NqP7TS0g58klJzcqGNwJEUhklmUcqhncP363PeoEowABYmpWG6M2cYVXN7VF3AH8QhxuZo62rZMnSjKkNxzlI6yEtDVrpaXeaz2otoBYoHlzssdwIybTmz1OEbF4csdwLmuUZJSW4dvSig5SuWuyHjJGdzVqCAs+sTksKgpOCZBxZRQjBONee2Z/zAoyv7mce9BIXDdg1otuQ+ujTBpuE69uTQs9ygjaUdepTa0goUoacYpwWjpOBjR3vX3e/GC8xefedTp26+139Qwm0Bi+Oxte5+YdpTkkdWOwBsTXew1jJMqmZSAIQgUJLVnkfkSVqBx6AdME01ha5k1A5b118ExlrasSIvDNZYnLA4C0hXZTuC+uLveH1rTlZqjHU4oBd5eBLGqUWIN1VkNwoHnIV+yyMtNN++OqYdKPKyutgf7UfMcsM1IbDO4XsCiaxsfAqLtY7UVO4Ms6yy/PnsmUU+cXyAkoLzO3M2pxmLrYDSWq6NqzpNoR1JrgmV2Bc9/MCjKxjjeG1rRuBJHjnS4/igBcDRQTVQzd1FY8DNAgBU2cvZrRmhJ3lwubP/u3ungoJmS+6jTRNsGq5jr/jvqWo41zzXOBxrgw65NvhKstKN3tKfc/ACc7t3pxd25vu1npVuxHc/tMR/8+kTPHdxiJIC5xxpadiaZsxyjbaObuSx0o04vz3HGUM7UFWPieex2gsIlWJrXoC2DFo+K92A89tzdmY5Skg6UcAozZDWYbQm01W9xhfQC/2qgdQ5ssJgbOVy4EuJtlVtZ5prIiUJPZBKkJWOQFQD3g6q2UprmWYlSWnwhKUfVgKKRFcD5XxZ9chgq6Fk1jkKDQ5L4EsCLNqAFZDlmiPdkHPbU04stPiRJ47wW99ZpxN6BJ5ipROirePqMOWN3YRpoRlvTLHO8UCueerEgJOLLa6NU2ZZVavxvar2NEpyLu4mPHdx+LaZw80CwO+9toUAVroRjx/p8u0rI751ecRTx/t4Sr7jJuNmS+6jTRNsGq5jL+NQsnIKTrXhgcUWEsc4LTmx0Nqv2cBbXQxuvDhtTXO+fmGXo/0YKQTPXxpypBdxtB8zzUr+/Tcu82OfPrEvOLi4m1RjpbVhmpW0A48k12xMc1S9NeWkoB+HrPSqYWDaWJx1FNrRCWCWGawF31P0WwFHuhG9xOPqKMEg8KSj5SscDs+TOOswxjHLq/2xyBdIIfCVIggE2lmUlMS+QihJWeaYegSCs1XjqHXUGYqoGlcVjHONFJVDsy8EE1edfy+o+qo6LtegjSVQkiiQCFeZjL62OefJE72q3rXU5uHVLto6LmynhF7Vf9Nv+SxlPhuznH7o8cSxAfOs5NefvcT3PLzEKC3ohj5SCNYnGdvTnI+t9Ti5GN+RO8DNAsA4LcHB6aUq83zqxALntma8sjHhU6cW33GTcWPk+dGmCTYNb2HQCviBR1d56sSA5y+POLs5JdWWJ47199Vot7Kc3ws0rcCrPLZ2EzwhSHLN2c05s6zk+EIbKSX9VgjAf/7OBg8st1nuhMzzEikEv/PqJgLBJjnTrMRTApylMJUV/6Nr3dpdwLLUCUhLzaVhyuVRCkIQBpJZLhjNC4azglAJPE/w0HKbYaYJlSQtDJ6CcV5V/feCRlY6rDUICbOsWjfSIQJFUZZYqqK+FIJ5UTkoOFvNoBFV4sJSHICAWaoZtAKUEvSi6qIP0PI9srKSRkde1dBZGovD0Y88wkBxbZSSac041ZxeanNtlHJpJ8FYxzCx5NowiAOO9ENOL7VpRz5pYZgXlY3QPDN8+tQif3B2m6P9iKVWwMmFmNi/c3eAmwUArR3uwCZiN/J56sSAYVLcctv0TriZo0XT4/PRoQk2DbdkL+j8wKOrN73tZheoSVriScEr6xNeXp8S7KnLCsO8KOnX2U47rP70upHPNy7s8sSxPpd2E+LAI/I9cmNxTjCcF2jjCH1J6PsIAb242iZb60f0Yx/jYFYHjMhTbE5zTL2l1QoUgadItSHNDZFXyYeh8iMbJRocRJ5Ee5VXQKGhtBDWY6yVrLKRpNBEvqIbKgrtkBIGLQ+JYF6USCXoBT6z0jBKSgJPVsFLV2Obj/YjhIDSWHqxR+hLSmPItENiEAhCVSnLYl/RjRS+kJzbmPJvvnYRT0oeWO4wSQvGqSYpBccXospUsx2w2ovZmub0Ih9PwfYs58xyh3ZQ9eJsTTNWOiHHF1p37A5wU0sjr25kOsBhBIXD9P1r+ODRBJuGQ0UIeP7KmEEc1N5jgtc2ZxztR2TasbU9I/YVK72Qdp39hH5V25jlml59EWx5Hhd255TWEUpJrit/tdVezMePdXl5fUo3qrzWro0y5kXJOCsZmYLQlySlRgqBlJLSOPLS4rCMkoJTSy2SwnBqMeb89py01GSmylAyXWU3Agh9hW/B96utNqerjAQJnQhagaIV+sSeYpgWzDJNXlqkqy7Iuh5lHfuK3FiujvNq1LWSnFhoo61FAjtJgbGWUCl8JevHloSBz7TO9L5xccTHjnQ4s9xhsR0yLzSvrk/YmOQc7cdsTTMu7CbMc81SOyT0JQt15rjSDfGU4NG17r7cfI+3CxI3CwD92EfAXTV+3gmNkedHmybYNBwKe6KA5y+P2ZhkBEoSB4q0MOTacHmY8MBii+cuZpQIzm9NWeqETNOSpW7AH57bZpZVBfrFdkgUVL0507QSCUhZ9aKkpebaOOPUYpvIU4yyglmpCTzFIPIZpgV5afCEIo4EQlbebM7BlVFCYSpjy4VWwMYko7SWTughRe07VhomuSY3UBjHwyttdmYFVlo85eGEI5BVE6bvVRJnzxMsdQLGqQZROQ7kpaEw1Wvphh6rgWJzlld9Ns7x2uYUoJoHNGjTiqoa0va0wGFRqsqkuqFP4MF3rk157pKurIIW22xPs0qBJgQPrXb49pURsVedY5wUZNry+U9VF/+taV47MOi3DEnbnGb0Y5//5fmr+xL1Y4P4BteI6wPAn3xkBbh5f8690hh5fnRpgk3DPXNQFNAJPKLFVtXcGXgY54g8iUOw1Al54tiArNRMMo0j49EjPWLf41uXhkyykpevFXzs6IDIk5xaaLEb5mzPSkpjCJXAl5JpVnKkF/HgSpvNqUcn9Li0M6coKtPO0BPszorKRSBUDOIAbS0Pr3SrLnkh6EQ+oV9Jj4dJTlZUBfogULRw9JWi4yuy0qKdI1SKqLb6DzzB8UHMvDQESnBpJ6m83ZRgmFvSsnI3GISCKFBYHHlp8GW1hSZF5dRc/exKcmlABRzphFDLtHNdDZoDGKW6tp0xvHh1yryoZOm+kjy82sZYx584s8jFnZRr44RWy+fRoz20sWxOMl64MuLJ4wNWuiGRpzi/MyfThm7kIYBSOzbGaV2D0lUQr+svB22KbqQJCg13QxNsGu6Zg4qlTuRRGsexvmCUlCx1Qr59ecRSx8dXgu95aIlu5OOc49kLO6x0Ii7uJlUdIS1BiMqdYKVFL6oaNb99ZcTWtKDUltCXPHlswGo/5OigqoFM0pydaYnG0Q4ks8JUzZ/OomTlu7Y9K4h9yadPLxAoxceO9vi91zb37VusKUlLyzTX+BJkIHC+wjnHIPKZFZqFVoAvJVLAKC352NEeDjg2aPH85REgaYU+S52QtLQY4zDWoY1jNys40g3xVXX/hZZPJ/KBOduzkkgJ+i2fvLRspxqNZZ6XhL5iZ1oQKoknYDfJENvVMLlHjnR5eLXLqxsz+nHAcsfQiz2WOyE7s4yvnks40o957EiX1V4EwGqvsr0JvMpNrtSOb14aMs01/SigF3sMa3fvRnLccJg0wabhnjmoWDo2iPnmpSEXduZMM8NqN8QBrcCnF/lcHaXM8hmFNlzYnnNlmBJ6ikE7AASLcUDmG5441if21b51TnhSkWu73/jpnOPSMGGSaY4PWrQf8/mtlza5tJvinKMbK8qyus8oKelHHnHg1Vtggqw0XNiZ0/IUTgiE0ChZd2oKQFbHZFpwpBsipMAYB87Sj336sc8jR7qc256z3PERCNb6IbtJQV5YulEl2a62whSnl1sEXjVG25OVD9tL1yYUZaUq256VDFqa1X7IicWYWa6ZF4btaUHkS44OWuAck7zcH4v98GoHJQUCy8tXp3hKcmwh5uooozSWT5zok2vLzjxnpRu+RRQwyzUbk4xpphm0fLRxXBtn5Nry6JFuIzluOFSaYNNwzxxULHVree9oXuB5kn7ssdzp8erGlP/fy+s8fWIBbQyvbUwQQmJzAw6+/sYuJxbbREqgpMe3L4/5xIk+57dnXN7NWOtHPHm8TzfymWYl57ZmnNuestaPSbSj0JYHllsYZ9iZaQZxSLuvMM5xdBBzpBOwOc1JS82gHTHNC4yxHBnEXB2lFNpSaIuQYEpQShOGAQutgEHLZ7UfM0wKssJUFjBK8PWLI4bTnMJWj781zSht1dw5CH0iP2BZCDqBTy/yWOtHxPW8mU7soYeWyIdAVbWW3XnOsYWYQMr9DOQbF4b4ShJ6ksIYzix32J5mTNNqLECuDdPcIJXk1FKbSVoSKomgMiBdaofMspKro5TH1qpgs1fDeXljghSCyJdoUzkpFNqSFrqRHDccOk2wabhnblQsXR2lHOnHPLrWo13XHXZmOeuTnMJYxqnm8aN9AL55cURaWOJAMU4KZMunXxfvd+cFz5xe4vRizrNv7PDClTHaWq6NUhbaIe3AZ3desD0rmGUlpXasdGIGLcNKLybNNcO0YGeWIxA8dXLAw6tdnIP1ccZTJxfAOTYmWdUwai3GVlJnTyhCJehEHhd3E8aZJvAEi3UGN00Kzm5MGbRC0kJTFI7Sg0BJSm1RUnBqscXHjvaQUnBmqUMcKH5FSea5ZnOc142lgjgQeFJSGMsbWzMeXu3QiUJ+5Ik1zm3N2Bhn+O1KXbY5qbKW1V5I6Kl6ZHebT50ccG2ccX57xlI74FS3hbaOY4OYl68V7MwLnHNsTfM3azjzkI1xytw4fCmI/So4K+k3kuOGQ6cJNg33zI2KJU9Kjg2i/UADYKmK6p85vcjXL+zuS5zPLLd5ZWOCJwTOc6z1e1wdp5xebKFNZWQZB4rAV+zMCrQ1BL7CORjOC1qBz5FuiLEOKTSBJ3FOsjnOUErQDRRR4CFwrPVjnKsysVmu+b6Hl/mdVzaRQrDaDUkKj1FSsNKLKLQhNyCEI9eOS7sJD660KzNPDPPSEUgPQRWQdnTJLNf4strKijxJYRyPH+2x0Ar4o/M7mNrAU1vHldG8ciQQEiEkgRJo49iaFXzmgZDPnllid15wbBCjraMdeOzOq9cUeB6x7/GdaxMKbVifpBzpxSy1Qx490sVX1cCbyFd0I5/TSx22ZznDpGB7lvPk8QGrvYhhUhB6iu1ZzjTXhL6HUoKjvajxI2s4dJpg03AoHFQsCQEvXR2TlXrfzdk6x2ovBqATeuS6chU7tdRiuRtwYTdBAP3YpzBVL8heEfvqKKUXeszzqp6w1A6ZZiW+9/9v78yD3MjOw/77+gAa5wBzcUgOhxzuktQe4l60VirrWMs61msnsitKIuewbDlWfMgVV1zlSNl/XHZcsaVKfMSOJcV2IidKSbJ8SKVal7xSrI3syq688kp7c0ku9+DwGs4MgMGNRr/88RogBgSGGM6FGb5fFTiN168b7/Vr9tfvfZeFEgVYJCM2FnCpUMWydIiZpaJPzLHY4zoECh4/s8Ab94+0dROzYwkeeMMeXr5S5HLBxxbt3NnwA+qNAIVwKV8nHXOo+xbJqMOVUo3FYo1qo8loIorrWGFYHQtRevlrNBFl34iHF7E5OJbg5fkis2NJXl8qsVxtUGs0iVg2QaBIRG1cW6jUm1QbAbdMJLllMsGrCyUmUh73HMjSaCoWS3Vc2yIVdShUfaZGPB1SKFchV64zlvAoVhuUaj71QBtDvHF/hsuFKmcXSuwd8dpCdiKl/W/2ZWIsV30OZOPU/CYzownOLpRIeg5nr5SYHWeFwDERmQ3rYd1poQ2Gbo5PZ5jOxmkGinylTjNQ3LlvhL0Zj3LdZ++Ip6NHV+rsHfHYn4mTjjq85fA4R/ekSEUd8tUG+zJaOM0Xa7y6WMa1LcYSetlqLlelGSgOjiZQShF1bY5Pj3Bsb5pGoHAsi/0jHpMprZOpNbRuyA8Ury2WmExFObtQZCTm8sDRPezNxsJkZTrlgWsL43FHpwlA4TcDyrUmrqX1UeW6j+sIQRAQsW0QwXEsElGXlKdzQe/PxHj8zIJOwZD2tBl1oHAdm6irzaiXq3VteCBwYCzGmw6Psliqc25Jh6VJxyLcPzvGsT0pUjGXiGNz+740UUfrtS4VatiWhWVpc27HscnEXG7fN0Ku0uDsQpHZsQQHsnHqfsCFfJX55Rqgozccm0oRoCjV/WvqPvXaErnQSKBl3l73A7LxyDX7DYbrsWNmNiLyIPDb6Gy/f6CU+vVtbpKhD5l4hLcemegblr7a0NZmoB/+o8kIP3zvdDs75MxYnHyl0Q5YmQs97PdndZj8V68UsW1YrvrYlp4pJT2HSiMg6th837EJbBEuFmoslWrkqzqiQD5MCe25OvT/3pEYNb/JpWIVC2EyFSMbj5AvN7SJsq2DpZXrPul4hMvFGq4teFGdxXS56pOJuVgS8PpSlcmUx1gygt8MODtf4s2Hxzh1abmdjuH58wVEdKDOhOeS8iJcLFSYy5WYGUvxpkNjlOsBL1xYoFxvcnq+zHtu38PekRjHpzNkExHOXC6iFLyyUGK5qlMLjMZ1eu27pjNtK7K3H53kqdeWGE9GVwTRnB1LcHahSNLTURtsSy9vzozGiTp234jLJiKzYb3sCGEjIjbwe8C7gXPA34nIl5VSz29vywz96OcIODtOO7inQjgymWwvx3Rmh2wt2SyV66Q8bTZtixB1LG2BVvfJlX2Wqw0UEAQK14KILTT8gFfzVWwR5pdrLFd8CtUmpXTAUvkit00liUVd3nFUe8LfMpFkOhPnmXP5cEaiWChqYwZB8P2AyWSEiudwbrFMNh5hJOZQrDcpVHVq5nTM5dBYvB2OZnYiwfxybUU6hnqYryZQipirTbpnsgnmizVu35filYVSaKYdYSRmUaw1+MaLl3jgDXsYibnsy8RIeQ5/e/oKflORjtkkvQieY+E6FudzFQ6MxttWZL2CaE6kolT9pk653eH9/93Xc8Rce0XdzrhpJiKzYb2IWksWpW1CRN4C/LJS6r3h948BKKX+Y79jUqmUuu+++7aohYZB8AMdSqXmBziWtMuijkUi6uAHipof0AwUtqUFi+faVBtNGs0AP3SStC3BsXXUZR0zLWibB7eOL9f9MF9MQCNQBIHSof/bH4vpbIxqQwfJdCyhHP6OUlD3ddIzEQmND3TMtEYzIBFxdHRnBU6Y88e2pB14M1CKiK2XtiZTUaoNnc9mqVRn7syLKBTZ6aM0lUIphSDaL6ehl+mi4czLsvQsSBDSMSd0AoVL+QqubaGASl0HII1FbJqBIh6x9XKaJRRr+hrYcjVoZlOpFREMWlyv7lrOZdjZPPbYY99WSp3Y6PPuFJ3NfuD1ju/nwrIViMiHReRJEXmy0WhsWeMMg1FtNPFDfYplSfthWvcDlsoNGs2Aht8kCBQNP6DRDChWG+EDVwenTHk6WrIAiahDxLYYibnEI3ZbQIEWYq5tEXHCwJaOhVg6/pn2LbF1W2zRkZgDLciCQL98iWjrMCc8ZzMApfQSXL2pM3amPL2sFI86uI4VzoTQwTQFbBEc2yLpuajwd0EHJw1CQaOU7pc2c9DGFUqp9myt0Qx0/LZQgDihbkah2xOL2G1B0/qtliD3XG2E0AxfKJtK98/rmsEMUnct5zIYerGrXkmUUp8GPg1w4sQJ9Y1vfGN7G2RYwWMnL3N2vkQ65iLhg1cpxXPnc0ymPcYSURrNAC/M9eLawoHROBHHYnY80dMSqqW07oxkfLlQ5esvXkSUMF+ssVCqtWdFe0Y8Hjg2SbOpmJ1IcNeBDN88Nc83T81Tq+tZTqHqMxLXVnFjCZ27p1j1da4cFKWaz9uPTJBNeCSjWsh9/YVLlMIMmemYiyUwM5pgNBlp+yHFIw4/9U//AYVKgzt/6j+zP+tx78EsaS/CUqnOs3M5Kr5POhZhIqk9/peKNbDgF951rG351RmLrjPqci9z5bVYkF2vrrFGuzmQjtnrRrJThM0ccKDj+3RYZthBpGMujqOXvlpvxDW/id9UjMYjK1IMRB2LQtVv6wX664CuDYFvWfDgnXtZKjX42vMXcR2LkbhDzLFxbKiHMdbSMZdMPMLx6Qx//eIl/CDAi1jsz6YYT0WZSns88uwFVLgcd3g8QcSxqdYbLJYb3H94HICTFwt44QwjX/Up1n3unx1jIhXl9aUyL11cJldpMJaIaOOAqMNDx6dYKNVJey5RxyIesTg0oROkjXguyahDqaojR99/cGyFIr5fKH6Ap15bukYYdF63lqDqJTBMxGXDZrJThM3fAUdEZBYtZD4A/LPtbZJhLeTKdZarDeaWKhRrPodG48QiNoVqg4lUlNFEFJE6Nb+pk6f5Acmofd2wKas9eL/47XPctjdNplBlPl8FFGkvwisLJd40O9q2kFss1bl3ZhTbEjy3lVitSbnR5NhUmkv5CuV6k1ikyfRonJgb55lzOZ6dy1Go+uRKNWp1HbctE49wPlfmW2cXSHsury6UyFfqJCIOY8koxZoPwJ6Uh4jg2haFaoNE1OGOtMehsQSvLZa5UqyTjbucODTJ1Ih3jSK+nxBppeOuNJrXpHwepM5q43ejxxoMsEOEjVLKF5GPAF9Fmz7/kVLquW1ulmFAOh9Ub711nOcu5Dk9v8z+bJzj05m242M2HuHVhSK1RpMAmEwleoZN6bWc052OOFeuU200WSrVUEoxkfao+QH5so/rCiMdAqxQaTA7nuClS8uU6z6Fik++UqdQ85lKe3iOxZ2hM+h8ocrUSKy9rBcoxZ6Ux9SIxdkrRVzbYjShZ2knLxYIFIzGI4glXF6uUW0E2Jb25ZlIRjk2ldax3q4UKRSqJCMO9xzItqM0g05Sdr04ZYOYJq/HfNmYPhvWy44QNgBKqUeAR7a7HYa10/2gesvhCe6azhJxrLaQGIm52gfHj7WTeLX0Hd16g0Hf4DNxl5hjcaVUZ365StS1OTSeYiKp44q1jkvHXOp+wP5MPAwroxXfgdLKfhDKtYCEZ1P3m5zPV9iT8rhj3wiea7f1TAulOvlKg3IjwA/AtixEBVxartJoaiu6ut/EsS3ylQYHRhMUKnWenssjKN64P0O1EejYZWTaCc4GiVM2iGnyesyXjemzYb3sGGFj2LkM8qBqLQt1z1C6WcsbfCLi8H9fmscSqNUDHEend751MrXiuJbe53Kxyi0TSUQI0wsEjHguMcdhvlhlqRSQ9mws26La8LljXxoR4eTFAqDNoFOey8HRGEvlBucWiliWUK43ycS1kYEfKEQU77xtD36geOq1RVJRh8MTSe1PFIM7yXClWMOx5ZosmP2U9J2Rt1t0L0EOUqcf6znWYICdY/ps2MG0HlSd3OiDqlBp9HQ+LFQaK+r4zYBnz+eZSkXJxCOU6j6Fco19mXi7buu4lt6n5uvcMq5tcWwqzUQyCiiUKA6OxXFsyFcapD2bu2eyOLYVhnxJ49rax8dzbI5NpVEKsskohapPxLZIRlySnoNtaZ8ZP1DcM5PllokUx6cz7VwzoB0vp0Y83nFssueMrVfImNlxveRYrvsopdrbLb0UMFCdfqznWIMBzMzGsAX0shi70RD2g77BPzeXpxnAWMrTUQCaoc9KM2gr6TuPy8Qj3DuTXWFGvS8DT5xdYLFY47a9aQ6PJ8lXG0ykvLaeCXS4nAOj8bb/j20JsYjO5fO6EuJRm2LdJwj0slw84rQFXnd/Wvqbmh+smLnkynW++uxFcpU6Y4loGE1At13P6rI9DSU6lyD7GVMMouBfz7EGAxhhY1gDN+pnsZEPqkEE1+x4gm++NE/Msag3gnaaghg2V0o1prOx9pt593Gd57ZDJ8oDY3HqTUhGLW4ZT3J5ucrXX7jIrZMpan6Taii03nZEh785e6VE1LU5PJ7CsSwuF7RhQCxq4draobMl5Dp/028GK/Q3rZnL4YkkL88XuZCvEASKC/kaL1wscP/sGFPpq5Zqg5gur8e82ZhGG9aDETaGgViv6etGPagGEVyZeIS7ZzK8eGGZVxdKjMQiHN2TYrnS4HKxRjyq0xf00oWU600uF2okPYd9mRiHJ5IcyMYREZarDU5eLBB1LCysdvKy7mtwz8xVR86pEY//d3qeqKMTkwVKsVCsc+byso7rpvRMqOY3OX15eYX+psXjZxYYT0Z14FGxyMZdFoo1vvzdOW6dTLE37ZELfZEMhmHFCBvDQAyT6esgguv4dIZ8uFR1uVDl1KVlJlJRfuJ7Z68J+Pn0uRzfeS3HaEILifFkdIU+orXMderysk4tXW+S8pwwFpmz4hp0zv5sS7hSrLJnJEap5rNYqBEEOmrCn3z7dabSce6dGWE6m8CyYE861hZsLWKuzfxyFUt0/pmLhWpomt3Q0Q7Kdd64b8T4vBiGHiNsDANxPYuyYQxlIkA8YrM/G8O1LTJxd4V/TWu2NperMJ6MIgIvXVrm2FS6LURaM5Ri1efkhTyJqIsTnuvkxQIjMZdnz+f421PzxKM6RP/seLI9+ytWm3zPoVHOzBdZrvoolBYSFZ+05/PqQoVA6RTSxarfUx81kfKYy1V1JINGwOu5MqJgLBlhPOUxmdZ5ggYV/L3GChi68TPsLow1mmEgVrMoG8bEWmev6EyXdx3IcuLQGHcdyDKR8jh7pbSiTjzi4Dd1oEnPdfBcm/O5yjWWaleKNVzHxrVtDo0lGE14lGo+X3n6PILF3pEYc0tVnnx1iXy5gYgQjziMJiKcXShx6lKR5TCfDqJD1pTrPhfzVTzXZrGkl+56WXzdvi/NhXyFy4UKVd8nX65TrNaZSHlMhlk3uy3y+tFrrP4mjA3Xa/xa9R87eXnbx9SwszHCxjAQq5m+di6xtR6yrZnBdtHPRPp8rtLx8FzEb+qwODU/ACDq2BRr/jWWalMjHu+5fYqJVLSd1O3MfIkgUMxOJLEsC9sSsrEIz8zl2r85O55gsVjjSqmGY1sESiFoCzbHEpZrPlHHYrGsM5PeM5Nt55ppOb36geLuAxkWyg3yoQBMei6nLhXaseQGNSV/+lyOuaUKL1wo8NKlZZqBIheet3v8nj6XG7qXCMPOxQgbw0C03vC7H4SZeGQg35etptdMbH65xoV8tf3wjNg2T8/lSXuuduJsNKk2fJ3bpsuHJB3Ty2fHplK4tlCo6sRtR/YkSbT1WDauDUvlq/12bEv75Fja+bPm63w5xUpDx3+L2BQqdWyLjsCZ2RU+NoWKTrNw/+wYb5hKM52Nk/Ac0jGXfKU+sM9LrlznO68tYQmkPZdGM+DkxQKFagPfX5nXKubanL68fM1LRBDAV5+9aGY6hjVjhI1hYHo9CGFjnTY3il4zsbMLRWbHEu2H5+GJJILilcUilihevJjnmbkco8nINcr21vlsSzi6J8Vte1NMj8bIJqLtOuMpj3ylQSJqr5j9OZaQr+hU046lTZ9zlToCpOMOvoJ33z7VV0eSjrkslrXO7OBYkrsPZLlrOstUOsalQnWF4F+Ns1dKjCajiFiISHvZMF/2cZyVYeUrjSYKWfESsVxt8OpiiVylbmY6hjVjhI1h3Qyjd3mvmdjekRgTqavCIeW5zI4neXm+RNUPuHPfCO+6be81s7R+5/tH9x6g4QfkyzWCIMD3m6Q9l7sOZNp1Dk8k+avnL7E/E2M6G9cJ0Gwhm4gScW1+8Pg+3n/f9AoLuW5mxxPYFhQqdZRSVBs+lsA9M1m+98jEwFZohUqD2bFEexbXSt5mWTo2Xff4HZlMrniJOJ+rYAFjiejQLJcadg7GGs2wbobVu7yXiXS3tdd8scad+0a468BV585+ll29zpeOuTx+ZoEL+QoTKY8P3H9wheB46rUlKvUme0c8MokokTDr56HxBJMpj7cfnRyoH+++fYpHn7/EfLHKaDzCzKg2l16LQG8FHD02leJ8rkKh6uPY8JZbxjg+nembHwf0stpCqY5raRPsFiYYp2FQjLAxbAib4V2+0ebUvaIPLBZr3HdwdEW9Xg/Qfm05OJZYdVZSqOh8PeVak6TnEnVsFBBxbGZG4wO3/eBYgvffN72u69GZMfTonlQ7+sLx6Uzf8et8icjEXMbDDKIttnu51LBzMMLGMJRsRrKuXjOwVkDN5WqD8zmd2M2xhJmxq4JgPW1Jx1yOTKZ44uUrACgUflPhNwPefMvYDbT/xoXtjcxAO3+zdR3KdX/dMe4MNx9G2BhWMCzOmZsVsaD7gZ0r1/mbU/OcWyqT9lwitkW+qk2BWyFgVmvL7PjqzpCz4wly5Tr3Hx7n1OVl6n6AZQn/4s0HGYm5fVM0d7NR47L+2GjDt1xq2Bmsy0BARD4hIi+KyNMi8ucikgnLD4lIRUS+E34+2XHMfSLyjIicFpHfkTA2h4iMisijInIq/JsNyyWsdzr8nXvX02ZDf4bJOXOrzKlbuWCSnku9CRHH4q7pDJMdDqDX89lZ7Xq1HtDZhMtowiUWzo7SoaAZ5FoP07h0WyQCxunTMBDrtUZ7FLhTKXUceAn4WMe+M0qpu8PPT3eU/z7wU8CR8PNgWP5R4OtKqSPA18PvAD/QUffD4fGGTWCYnDO30pxaKTi+P8N9B7Mcm0qT8twVgq1fW4pVf+Dr1QwUd+zLkPYcRODR5y8RBAx07DCNSyfDJAQNw8+6hI1S6q+UUn749XFgerX6IrIXSCulHldKKeCPgR8Od78P+Ey4/Zmu8j9WmseBTHgewwYzTM6ZW2lOfT3B1q8tSc9Z9Xq1HsZ/9vevM7dUoRlox0lbhGYQsFiq9T22k2Eal06GVQgahpON9LP5EPCXHd9nReQpEXlMRN4Wlu0HznXUOReWAexRSl0Ity8CezqOeb3PMYYNZJicM1eLWLDRXE+w9WvLvkxsoHhxFhaWwMmLhbbAGY1HWCw3eh7bzVrHZavimQ2rEDQMJ9c1EBCRrwFTPXY9rJT6UljnYcAHPhvuuwDMKKUWROQ+4C9E5I5BG6WUUiKirl/zmrZ+GL3UxszMzFoPv+nZyIyaG8FWJesaNEdOd1tmx+l7vTrf+pOeQ6Op8Fyh5gfEIzajiSj5amMgy661jMtmWPH1Y5CsqQZDi+vObJRS71JK3dnj0xI0Pw78EPDPw6UxlFI1pdRCuP1t4AxwFJhj5VLbdFgGcKm1PBb+vRyWzwEH+hzT3dZPK6VOKKVOTExMDNB9QydbOZvYTjbqzX/QeHH7MrHQY1/rbppKYVk6RM0g13ot47KVS1vDGDnCMLysy/RZRB4Efgl4h1Kq3FE+ASwqpZoichit3H9ZKbUoIgUReTPwBPBjwH8JD/sy8EHg18O/X+oo/4iIfA64H8h3LLcZNpitmk1sF73e/L95ah4BJlLemmcD/a5X51t/ynM5NpXi5StFFGCJtM+9mkPoIL/TzfXyDt0o/UyvjSm0YVDW62fzu0AUeDS0YH48tDx7O/ArItIAAuCnlVKL4TE/C/wPIIbW8bT0PL8OfEFEfhJ4FfgnYfkjwEPAaaAM/MQ622zYQWy0308vn5l8pQEKDo4l286dC6Ualws13ntn/wCZq9G99GVbwv6Mjs3mWLJpD+TNWNrqJ6BHYi5K6d+860DGCBnDqqxL2Cilbu1T/qfAn/bZ9yRwZ4/yBeD7e5Qr4OfW007DzmQ9+od+QqrXm7/vKxSK5WqDkxcLeK7NeCLKlVL9hvQdrd8u15tcLuikaPsyMY5N6VQDm8lm6N26BXQzUJxbqpAvNzg+ndlUvZBh92CiPhuGlhvVP6zm/9HLsstxBNe2OJ+rtDN21puKsUSk5++tpvPp/O0D2TiHxhPEI/aWRWLYDL1bt9XZ+VyFEc/FD5QxeTYMjAlXYxhablT/sHp4mWvf/EdiLgLM5SrEXZtTiwXylQbHplL4zYBqh3C63mxrs8LsrIWN1rt1L80Va00its422sJEfzZcDzOzMQwtN+r3s5r/R683/7cdmeCtRyZwbeGZuTwAb9ibxrUtnpnLIR0rX9ebbe1G35NuqzPHhkK1sSLVgDF5NlwPM7MxDC2dsxC/GXB2ocRiscbdM9l2kMxeXE9J3u/N/9bJFDVfkYlFiDqWTuHMSh3L9WZbu9H3pNvqbGY0TqHSwLYEpdS2+2MZdgZmZmMYWloPuZrf5NuvLoGC+w6OEnXsVf1jbtT/Q8dIG8G1hULVx7WF4/tHUB3uxTca2man+550BuB8+9FJ3npkYtf7Yxk2FjOzMQw1mXiElOdy/+GxFbMF6K8HycQjHJ5I8viZBeaXq0ykPN58y1jPh2Gn1drFfJXxZJRjU+n2/nLdJx69+k52PWuvtfqeDEtKh7Wy2/2xDBuPETaGoWethgK5cp2X54scGk9w2940lUaTl+eLjMTcFQ/ybmW/31Q8O5fjTjJMpKI9l4duNLRNv3ZuVWgZg2G7McLGMPSsVQ8yqEVYd73JtMedZLhSrOHY0ndWslFv9cNguWYwbBVG2BiGnrU6Kg46E+pVbyIVxbGFdxyb3MAe9GazQssYDMOIMRAwDD1rdVQc1GR6u1MqbPfvGwxbiZnZGHYEa1m6GnQmNDue4G9OzZMrN2g0A1zbIhN3eeuRjY0Y3mkEUKz5eKEfzrCldDAYNhMjbAy7jtWU+J0PfhEo1XwQEP0Pa06idB26jQACpShWG20/IRM12XCzYISNYVfSaybU/eB/ei5Hsdrg+HSWlKeXrsp1f0MV9N1GALYIWFeNAIwJseFmwQgbw01D94Pfb0Laczmfq3BsSgubjVbQ9zICsEU2PHzNTvXXMdw8GAMBw01Dd9yyZNQGhGLNb5dttIK+lxFAUynSMXfDMoauFuXaYBgWjLAx3DR0P/j3ZWLkqw2cMMbXZoSW6Q5f01SKIFCMJiIbJiC2MhW0wXCjGGFjuGnofvDbljCdjTEzFt+0GF/dZtuWCEnPZbFU3zABsRsjTRt2H+sSNiLyyyIyJyLfCT8Pdez7mIicFpGTIvLejvIHw7LTIvLRjvJZEXkiLP+8iETC8mj4/XS4/9B62my4eemXXuDtRyd5x7HJTQsT0xnEMhl1cCzZUAFh/HUMO4GNmNn8plLq7vDzCICI3A58ALgDeBD4ryJii4gN/B7wA8DtwI+GdQF+IzzXrcAS8JNh+U8CS2H5b4b1DIYbovPBv50xyDZSQOzWSNOG3cVmLaO9D/icUqqmlDoLnAbeFH5OK6VeVkrVgc8B7xMRAd4JfDE8/jPAD3ec6zPh9heB7w/rGww7lo0UEJuRCtpg2Gg2wvT5IyLyY8CTwC8qpZaA/cDjHXXOhWUAr3eV3w+MATmllN+j/v7WMUopX0TyYf0r3Q0RkQ8DHwaYmZlZf88Mhk1iox06jb+OYdi5rrARka8BUz12PQz8PvCraMfrXwX+E/ChjWzgWlBKfRr4NMCJEyc22hncsM3sNl8SIyAMNxPXFTZKqXcNciIR+W/AV8Kvc8CBjt3TYRl9yheAjIg44eyms37rXOdExAFGwvqGmwiT+8Vg2Nms1xptb8fXHwGeDbe/DHwgtCSbBY4A3wL+DjgSWp5F0EYEX1ZKKeCvgfeHx38Q+FLHuT4Ybr8f+D9hfcNNRMuXpBkoXrq0zAsXlpnLVXj6XG67m2YwGAZgvTqbj4vI3ehltFeAfw2glHpORL4APA/4wM8ppZoAIvIR4KuADfyRUuq58Fz/DviciPwH4CngD8PyPwT+p4icBhbRAspwk1GoaOfLly4V8VybtOdQbTT5zms5jk9ndsXsZrctExoMnaxL2Cil/uUq+34N+LUe5Y8Aj/QofxltrdZdXgX+8Xraadj5pGMuz53P47l2O0S/CIwmIrsis6VZJjTsdkwEAcOOYHY8wWKxhlIBSimqDZ9qo8nseGJXeMqbkDOG3Y4RNoYdQSYe4e6ZLIGCQrWBa1scm0rj2Nau8JQ3IWcMux2TYsCwYzg+naEZKOIRZ9dltmxFFGilPwATcsawuzAzG8OOYTd7ypuQM4bdjpnZGHYUu9UR0qSINux2ZLe6rIjIMnByu9uxiYzTI2TPLmI392839w1M/3Y6x5RSqY0+6W6e2ZxUSp3Y7kZsFiLypOnfzmQ39w1M/3Y6IvLkZpzX6GwMBoPBsOkYYWMwGAyGTWc3C5tPb3cDNhnTv53Lbu4bmP7tdDalf7vWQMBgMBgMw8NuntkYDAaDYUgwwsZgMBgMm86OEjYi8vMi8qKIPCciH+8o/5iInBaRkyLy3o7yB8Oy0yLy0Y7yWRF5Iiz/fJhbhzD/zufD8idE5NCWdlC34RdFRInIePhdROR3wjY9LSL3dtT9oIicCj8f7Ci/T0SeCY/5HRGRsHxURB4N6z8qIlsW50VEPhGO3dMi8ucikunYt2vG73r069OwISIHROSvReT58P/bvwnLe95DG3mfbnE/bRF5SkS+En5f87211vt3C/uWEZEvhv/vXhCRt2zr+CmldsQH+D7ga0A0/D4Z/r0d+C4QBWaBM+hcOXa4fRiIhHVuD4/5AvCBcPuTwM+E2z8LfDLc/gDw+S3u4wF0rp9XgfGw7CHgLwEB3gw8EZaPAi+Hf7Phdjbc962wroTH/kBY/nHgo+H2R4Hf2MK+vQdwwu3faP32bhq/Aa5B3z4N2wfYC9wbbqeAl8Kx6nkPbeR9usX9/LfA/wa+ciP31o3cv1vYt88A/yrcjgCZ7Ry/bb+p13DhvgC8q0f5x4CPdXz/KvCW8PPV7nrhhbnC1Qdfu17r2HDbCevJFvbxi8Bd6ER0LWHzKeBHO+qcDB8EPwp8qqP8U2HZXuDFjvJ2vdax4fZetOPrdozljwCf3W3jN0C/e/Zpu9s1YNu/BLy73z20kffpFvZpGvg68E50Svs131trvX+3sG8jwNnu+387x28nLaMdBd4WTmEfE5HvCcv3A6931DsXlvUrHwNySim/q3zFucL9+bD+piMi7wPmlFLf7dq11v7tD7e7ywH2KKUuhNsXgT0b0/o18yH0mxDskvEbkH59GmrCJaN7gCfofw9t5H26VfwW8EtAEH6/kXtrrf3eKmaBeeC/h8uEfyAiCbZx/IYqXI2IfA2Y6rHrYXRbR9HTtu8BviAih7eweevmOv379+ilpi1BKaVEZEPt3lfrn1LqS2Gdh9Gpwj+7kb9t2BxEJAn8KfALSqlC57L8ZtxDW4WI/BBwWSn1bRF5YJubsxk4wL3AzyulnhCR30Yvm7XZ6vEbKmGjlHpXv30i8jPAnyk9Z/uWiATogHhzaF1Hi+mwjD7lC0BGRJzwDaWzfutc50TEQU9FF9bdsZB+/RORN6LfRL4b/meeBv5eRN5E//7NAQ90lX8jLJ/uUR/gkojsVUpdEJG9wOV1dmkFq40fgIj8OPBDwPeH4wg7aPw2gNX6OnSIiIsWNJ9VSv1ZWNzvHtrI+3Qr+F7gH4rIQ4AHpIHfZu331lrv363iHHBOKfVE+P2LaGGzfeO3lWuk61yD/GngV8Lto+ipnQB3sFJB9zJaOeeE27NcVdDdER7/J6xUAv5suP1zrFQCfmGb+voKV3U2P8hKxd23wvJR9JpsNvycBUbDfd2Ku4fC8k+wUjn48S3s04PA88BEV/muG79VrkHfPg3bJ7x3/hj4ra7ynvfQRt6n29DXB7hqILCme+tG7t8t7Nc30RGcAX45HLttG79tv6nXcOEiwP8CngX+Hnhnx76H0ZYfJ+mwiEBbWLwU7nu4o/xweKFOhzdXy8LNC7+fDvcf3qa+vsJVYSPA74V9eAY40VHvQ2FbTwM/0VF+IrxOZ4Df5WqkiDG0QvQU2rJvdAv7dBr9gvCd8PPJ3Tp+17kOPfs0bB/grYACnu4Ys4f63UMbeZ9uQ18f4KqwWfO9tdb7dwv7dTfwZDiGf4EWFts2fiZcjcFgMBg2nZ1kjWYwGAyGHYoRNgaDwWDYdIywMRgMBsOmY4SNwWAwGDYdI2wMBoPBsOkYYWMwGAyGTccIG4PBYDBsOv8fSagFxo+2ibUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 1)\n", "ax.scatter(aeps_N1_areas, aeps_P1_areas, alpha=0.2)\n", "ax.set_xlim(-60000, 60000)\n", "ax.set_ylim(-60000, 150000)\n", "ax.axhline(0, color='black')\n", "ax.axvline(0, color='black')" ] }, { "cell_type": "code", "execution_count": null, "id": "e5daed50", "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 }