diff --git a/SiPM_peaks_mejorado.ipynb b/SiPM_peaks_mejorado.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..7209677bb8435ffea826ae74a938a34941330823 --- /dev/null +++ b/SiPM_peaks_mejorado.ipynb @@ -0,0 +1,440 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "id": "1e9e0f71", + "metadata": {}, + "outputs": [], + "source": [ + "import numpy as np\n", + "import pandas as pd\n", + "import matplotlib.pyplot as plt\n", + "import os\n", + "\n", + "%matplotlib inline" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "id": "054e11a4", + "metadata": {}, + "outputs": [], + "source": [ + "plt.rc('axes', labelsize=20)\n", + "plt.rc('xtick', labelsize=20)\n", + "plt.rc('ytick', labelsize=20)" + ] + }, + { + "cell_type": "markdown", + "id": "235c5c85", + "metadata": {}, + "source": [ + "### Lectura de los datos\n", + "\n", + "- Datos obtenidos para la caracterización del SiMP." + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "73c4c6ff", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "'/home/christian/MEGA/PostDoc_UIS/Conga_2022/Notebooks_CONGA/SiPM'" + ] + }, + "execution_count": 3, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "#A partir de aquà python sabe en que carpeta se encuentran los datos\n", + "os.chdir(\"/home/christian/MEGA/PostDoc_UIS/Conga_2022/Notebooks_CONGA/SiPM/\")\n", + "os.getcwd()" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "18a64f20", + "metadata": {}, + "outputs": [], + "source": [ + "#Lectura de los datos\n", + "ds=pd.read_csv(\"picos_285292.tmp.txt\", sep=\"\\t\")" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "id": "4bd27abe", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "<div>\n", + "<style scoped>\n", + " .dataframe tbody tr th:only-of-type {\n", + " vertical-align: middle;\n", + " }\n", + "\n", + " .dataframe tbody tr th {\n", + " vertical-align: top;\n", + " }\n", + "\n", + " .dataframe thead th {\n", + " text-align: right;\n", + " }\n", + "</style>\n", + "<table border=\"1\" class=\"dataframe\">\n", + " <thead>\n", + " <tr style=\"text-align: right;\">\n", + " <th></th>\n", + " <th>ADC Channel - Charge distribution</th>\n", + " <th>Counts - Charge distribution</th>\n", + " </tr>\n", + " </thead>\n", + " <tbody>\n", + " <tr>\n", + " <th>0</th>\n", + " <td>-249</td>\n", + " <td>0</td>\n", + " </tr>\n", + " <tr>\n", + " <th>1</th>\n", + " <td>-247</td>\n", + " <td>0</td>\n", + " </tr>\n", + " <tr>\n", + " <th>2</th>\n", + " <td>-245</td>\n", + " <td>0</td>\n", + " </tr>\n", + " <tr>\n", + " <th>3</th>\n", + " <td>-243</td>\n", + " <td>0</td>\n", + " </tr>\n", + " <tr>\n", + " <th>4</th>\n", + " <td>-241</td>\n", + " <td>1</td>\n", + " </tr>\n", + " </tbody>\n", + "</table>\n", + "</div>" + ], + "text/plain": [ + " ADC Channel - Charge distribution Counts - Charge distribution\n", + "0 -249 0\n", + "1 -247 0\n", + "2 -245 0\n", + "3 -243 0\n", + "4 -241 1" + ] + }, + "execution_count": 5, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "ds.head()" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "id": "57e0f701", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "<class 'pandas.core.frame.DataFrame'>\n", + "RangeIndex: 2936 entries, 0 to 2935\n", + "Data columns (total 2 columns):\n", + " # Column Non-Null Count Dtype\n", + "--- ------ -------------- -----\n", + " 0 ADC Channel - Charge distribution 2936 non-null int64\n", + " 1 Counts - Charge distribution 2936 non-null int64\n", + "dtypes: int64(2)\n", + "memory usage: 46.0 KB\n" + ] + } + ], + "source": [ + "ds.info()" + ] + }, + { + "cell_type": "markdown", + "id": "a6a7ae18", + "metadata": {}, + "source": [ + "### Ajuste gaussiano de cada uno de los picos" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "id": "15d0c568", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAApkAAAHyCAYAAACgdNVWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABt5ElEQVR4nO3dd5gURf4G8PdLzkkkCEpSTIh6YsK0mAPq3emZFfOZz3iGM6Bnzp4YzvATlVNMZ8B8IisqKsGAEVCSIFGUnKnfH9VD1/T0zHTPdJx5P8+zz3T3dKjtHdh3q7qqRCkFIiIiIqIg1Ym7AERERERUeRgyiYiIiChwDJlEREREFDiGTCIiIiIKHEMmEREREQWOIZOIiIiIAlcv7gJUg7Zt26quXbuGeo2lS5eiadOmoV6jGvG+hof3Nhy8r+HgfQ0P7204yrmv48aNm6+U2rDcMjBkRqBr164YO3ZsqNeora1FTU1NqNeoRryv4eG9DQfvazh4X8PDexuOcu6riEwLogxsLiciIiKiwDFkEhEREVHgGDKJiIiIKHAMmUREREQUOIZMIiIiIgocQyYRERERBY4hk4iIiIgCx5BJRERERIFjyCQiIiKiwDFkEhEREVHgGDKJiIiIKHAMmUREREQUOIZMIiIiIgocQyYRERERBY4hk4iIiIgCx5BJRERERIFjyKTUevBBoH9/YPTouEtCRERETvXiLgBRKSZNAs49Vy+/9x6wYkW85SEiIqJsrMmkVDJrL1eujK8cxdx+O9ClC/DQQ3GXhIiIKFoMmZRKa9fGXYLili8HLr8cmD4dOOecuEtDREQULYZMSqU1a+IuQXHLlsVdAiIiovgwZFIqMWQSERElG0MmpVIamsuvuy7uEhAREcWHIZNSKQ01mS+8EHcJiIiI4sOQSank7K2tVDzlKGTJkuz1JJaRiIgoLAyZlErff5+9/vzz8ZQjH7eaVoZMIiKqJgyZVBEeeyzuEmR76qncbevWRV8OIiKiuDBkUkWYODHuEmT7+efcbQyZRERUTRgyqSJMnw6MHRt3KWxuvd8ZMomIqJowZFLF6N8/7hLY3AIlQyYREVUThkyqGHPmxF0Cm1snH3b8ISKiasKQSRSCb7/N3caaTCIiqiYMmUQhePXV3G0MmUREVE0YMokiwpBJRETVhCGTKCIMmUREVE0YMokiwo4/RERUTRgyiQK2erX7dtZkEhFRNWHIJArYiBHu2xkyiYiomjBkEgXs6qvdt69YEW05iIiI4sSQSRSwMWPctz/0ULTlICIiihNDJlFEbr897hIQERFFhyGTiIiIiALHkElEREREgWPIJCIiIqLAMWQSERERUeAYMil1pkyJuwRERERUDEMmpc6ll8ZdAiIiIiqGIZNSZ/HiuEtARERExTBkUuo0aBB3CYiIiKgYhkxKnYYN4y5BaXbbLe4SEBERRYchk1Knfv24S1CaZs3iLgEREVF0GDIpdZIcMpXK/96qVdGVg4iIKG4MmZQ6SQ6Zn3yS/z2GTCIiqiYMmZQ6SQ6Zixblf2/16ujKQUREFDeGTEqdevXiLkF+devmf481mUREVE0YMil1mjSJuwT51SnwL4ohk4iIqglDJqXOpElxlyA/1mQSERFpDJmUOq++mv+9deuiK4ebQiGTz2QSEVE1YcikirJ2bbzXZ00mERGRxpBJFSXukMlnMomIiDSGTKoocYdMNpcTERFpDJmUanffnb0e9zOZIvnfY00mERFVE4ZMSp3ttrOX99gDaNnSXo+7JrNQyGXIJCKiasKQSaljBrn69bObqOMOmYXmLl+3Lv7yERERRYUhk1JnzRp7uV697M42cYc4t5rMBg3sZT6XSURE1YIhk1LHGTKTVJPpDJkDBmTPtc4mcyIiqhYMmZQ6hUJm3B1/nM3ld9/NmkwiIqpODJmUOmmpydx9d6BNm+yQyZpMIiKqFgyZlDpJDpm//movL1igX9lcTkRE1Yghk1InyR1//v53e/m77/Qrm8uJiKgaMWRS6jhDphks4w6ZP/6Yu80MmStXRlcWIiKiODFkUuqYQbJePWDaNHv93XejL08x5mDxc+fGVw4iIqIoMWRS6jhrMk0XXBBtWbzo2tVe/uWX2IpBREQUKYZMSp1CITOJGje2l9lcTkRE1SL2kCkiG4jI6SLysoj8KCLLRWShiHwkIqeJiGsZRaSviLwpIgusY8aLyIUiUtdtf+uY/iJSa51/iYh8JiIDipRvgIiMtvZfaB3fv9zvm0qXtpDZqJG9vGJFfOUgIiKKUuwhE8BfADwKYGcAnwG4F8BLAHoBeAzA8yIi5gEicjiAkQD2BPAygEEAGgC4B8BQt4uIyHkAhlnnHWJdcyMAg0XkzjzH3AlgMICO1v5DAGwDYJh1PoqYUtnPZNZJwie4iIYN7WXWZBIRUbVIQj3QRACHAXhDKbV+KGsRuQrAaABHAPgzdPCEiLSADnxrAdQopcZa268B8D6AI0XkGKXUUONcXQHcCWABgD5KqanW9hsAjAFwiYi8pJT6xDimL4BLAPwEYEel1G/W9jsAjANwp4i8njkXRcPZ6Sf7z49kYk0mERFVo9jrgZRS7yulhpkB09o+G8DD1mqN8daRADYEMDQTMK39VwC42lo923GZUwE0BDDIDIVWcLzZWj3LcUxm/aZMwLSOmQrgAet8pxT/DilIaWsqB1iTSURE1Sn2kFlEZuhqI1pgb+v1bZf9RwJYBqCviBi/2gse85Zjn3KOoZClPWRef3185SAiIopSYkOmiNQDcJK1aga9za3Xic5jlFJrAEyBfgygu8djZgFYCqCziDSxrt0UQCcAS6z3nSZZrz09fTMUGDNk1nXp4tWpU3Rl8coZhpcti6ccREREUUpyXdCt0J103lRKvWNszwxtvTDPcZntrXwe09Tab1mJ18giImcCOBMA2rdvj9ra2ny7BmLJkiWhXyMJ5s1rCGBXAED9+itRW/sJevXaHt98o39kffv+jNranwK7nt/72qjRHlixwk6/tbW1+P77LgC6rd/2/vsfoVmzNS5HV5dq+cxGjfc1HLyv4eG9DUcS7msiQ6aIXADd6eYHACfGXJySKKUeAfAIAPTp00fV1NSEer3a2lqEfY0k+OYbe7ldu4aoqanBKacAl1yit3XqtDFqajYO7Hp+7+uBBwKvvKKX77sPqKmpweuvZ++z0067o127wIqYWtXymY0a72s4eF/Dw3sbjiTc18Q1l1tDA90H4DsA/ZRSCxy7ZGoRW8JdZvvvJRyz0PHq5xoUgd9/t5dbtdKvZnP0mpgrCDMBEwC6dNGvq1dn77NqVWTFISIiik2iQqaIXAjgfgDfQAfM2S67TbBec56HtJ7j7AbdUWiyx2M6QjeVz1BKLQMApdRSADMBNLPed9rMes15xpPC9dtv9rJbyFyY7wGHCHzxRfZ6plzr1mVvZ8gkIqJqkJiQKSKXQw+m/iV0wJybZ9f3rdcDXd7bE0ATAKOUUuZgMYWOOcixTznHUMjcajLNDkBPPx1labKZTfmAXS6GTCIiqkaJCJnWQOq3Qg9yvo9San6B3V8EMB/AMSLSxzhHIwA3WqsPOY55AsBKAOdZA7NnjmkN4Cpr9WHHMZn1f1j7ZY7pCuBc63xPFPveKFiLF9vLzZvr11lu/f8TIF/IdDafExERVaLYO/5Yc4ffAD2Dz4cALpDcaVymKqUGA4BSapGInAEdNmtFZCj0TD6HQQ9V9CKA58yDlVJTROQyAP8CMFZEngOwCnpg984A7jJn+7GOGSUidwO4GMB4EXkReurKowG0AXA+Z/uJnjnjT/368ZXDjfNjy5pMIiKqZrGHTNhju9QFcGGefT6AnkMcAKCUekVE9gLwD+hpJxsB+BE6EP5LKaWcJ1BK3S8iUwFcCj3+Zh3ozkVXK6WedLuoUuoSEfkauubyTADrAHwO4A6l1Otux1C4zMCWtHnLnSEz80ymGYwB1mQSEVF1iD1kKqUGAhhYwnEfAzjY5zHDAAzzecxgGAGX4uUWMpM6f3mmJvP444HHH7e3O0MnERFRJUpYXRBRYWkMmf36ZW9nyCQiomrAkEmpkqaQaQ6ttNde9jJDJhERVQOGTEqVJIfMfB1/nMsMmUREVA0YMilV0tTxhyGTiIiqWcJ+TRMVluSaTCeGTCIiqmYMmZQqbiEzKTWa+YYwAhgyiYio+iTk1zORN2mtyTSDMEMmERFVA4ZMShW3kNmzZzxlKSZfc7lzBiAiIqJKxJBJqeIWMv/0p+x9cud7igefySQiomrGkEmpkq+53Gwyjytk8plMIiIiG0MmpUq+IYySGDJZk0lERNWMIZNSxUvITMozjwyZRERUzRgyKVXyhUxzmTWZRERE8WPIpFRJcnO5E5/JJCKiasaQSanyzDP2shkmkxAyWZNJRERkY8ikVJk9215+9ll7mc9kEhERJQtDJqXWTz/Zy0l4JtOJIZOIiKoZQyallhkmly1z3x6nfNNKJqWmlYiIKEwMmZRaDzzgvn3s2GjLkVHH8a/JbMJnTSYREVUbhkxKldat7eV99nHf59ZboylLIe3aZa8zZBIRUbVhyKRUMQNa48bu+yxYEE1ZnMxm8N13z36PIZOIiKoNQyalyvLl9nK+kBmXfEMqAQyZRERUfRgyKTXWrgVWr9bLIkCDBvGWx8kMmc7nMxkyiYio2jBkUmo4azGdtYVxM5vLWZNJRETVjiGTUiPJTeUAMG+eveysyTTXP/ggmvIQERHFiSGTUiPpIfPSS+3lL7/Mfu/FF+1lhkwiIqoGDJmUGkkPmatW2cs//JD93tdfR1sWIiKiuDFkUmokPWQSERGRjSGTUsNryExahyAAOPLIuEtAREQULYZMSo0012Refrm93KZNfOUgIiKKCkMmpcbKlfZy0sbILMacDrNly/jKQUREFBWGTEoNcxxKc9xJp622Cr8sxTz+ePZ6w4b2shmWiYiIKhVDJqWGOYi5M2Seeaa93KJFNOVx6tnTXu7bN/s9hkwiIqo2DJmUGoVqMnv1ct8vSuYQRvXrZ7+XtJD5+efAllsChxxiT9VJREQUJIZMSg2zJrPQ3OBJCJlmqASynyFNQsjcbz89luebbwKDBsVdGiIiqkQMmZQahWoyzdAZ19zgZsh0dkwy11evBpSKpkz5LFhgL3/4YXzlICKiysWQSalRqCbTXE9CTaYzZNapk92Ebu4bt5dfjrsERERUiRgyKTUK1WQmvbncuS0JTeZERERhYsik1EhzTaZzW5JqMomIiMLAkEmpkeRnMteutctXp477OJ5JCML5xP2MKBERVR6GTEqNQjWZcffeLlaLCSQ7ZL7+etwlICKiSsOQSalRqCazaVN7eenSaMpj8jLlpRkyk1ZzeMUVcZeAiIgqDUMmpUahmswmTezlZcuiKY/JS03m7Nn28owZ4ZanELfHCX79NfpyEBFRZWPIpNRIck2ml5BpOvvs8MpSjFunI/Z2JyKioDFkUmoUqsls1MheXrEimvKY/IbMb74JryzFMGQSEVEUGDIpNcyQWWiczDh6l/sNmXFym6t8+fLoy0FERJWNIZNSw2wud9Zk1qtnL69ZE015TH5DZpwdfzhGJxERRYEhk1IjLTWZbrP9OMUZMr/6yn37qFHRloOIiCobQyalhllDadZcAskKmV5qMuMcJ/Pgg/1tJyIiKgVDJqVGodpCNpeXb+HCuEtARESVhCGTUqPQgOdpq8lMYsgkIiIKEkMmpUahIBd3yPQy4w8REVE1Ycik1GBzORERUXowZFJqVFJzORERUaVjyKTUWLLEXmbIJCIiSjaGTEqNSZPs5Y02yn7PbC5nyMwvzqGTiIioujBkUmqYc5J37Jj9nlmTuXRpNOUxmdMyNm4c/fW9eumluEtARETVgiGTUsOsoXQOxl6/fvb63Lnhl8dkhsxGjaK9th933BF3CYiIqFowZFJqFJrxp149oGVLe33+/GjKlGHWsia5JjOORwmIiKg6MWRSapgh0zl3OQBssom9vHp1+OUxpaW5/PPP4y4BERFVC4ZMSo1CzeVAdpN51GNlpiVkEhERRYUhk1KjUHO5c1vUITMtzeVERERRYcik1CjWXB5nyExDx59Zs+IuARERVROGTEqNJNdkpqG5/LffcreddJK93L17dGUhIqLKx5BJqVHsmcwJE+zlL78MvThZvDSXmx2T4iCSu+3KK+1lt3tKRERUKoZMSo1iNZlz5tjLF14YenGymCEzX3N53CHOLWTG2VmKiIgqG0MmpUaxZzLjZA6Z5BwYPiPuMruFzDgfMSAiosrGkEmpUawmM05mU36+MOmc0zzqecSLhczZs6MrCxERVb5EhEwROVJE7heRD0VkkYgoERmSZ9+u1vv5voYWuM4AERktIktEZKGI1IpI/wL71xWRi0RkvIgsF5EFIvKmiPQN4vsm75TKri10BjanPfYItzxOXgJw06bZ62bzfhSKhcxVq6IPvkREVLmSUh90NYBtASwBMAPAFh6O+QrAKy7bv3HbWUTuBHCJdf5HATQAcAyAYSJyvlJqkGN/ATAUwJEAJgAYBKANgKMBjBSRI5RSr3ooJwVg1arsdbfAdO21wA036OUddgi/TCYvNZnOKR2XLg2vPF45g+/MmcDGG8dTFiIiqixJCZkXQYe/HwHsBWCEh2O+VEoN9HJyq+bxEgA/AdhRKfWbtf0OAOMA3CkiryulphqHHQMdMEcB2EcptcI65mEAHwF4VETeV0ot9lIGKs/w4cX3advWXo76+UIvNZkdO2avr1wZXnm8atYse/377xkyiYgoGIloLldKjVBKTVJKqZAucZb1elMmYFrXnQrgAQANAZziOOZs6/XqTMC0jhkD4DkAG0KHUIrAo48W38esQYw6ZHqpybz//ux1s0d6FLw0hR91VPjlICKi6pCIkFmijUTkryJylfXau8C+e1uvb7u895ZjH4hIIwB9ASwD8KGXYyhc+Xpsm8waRGfTdNi89Hzv2hXYbDN7PeqaTC8hc+HC8MtBRETVISnN5aXYz/paT0RqAQxQSk03tjUF0AnAEqWU28R6k6zXnsa2HgDqApislHKrE3M7hkLkpTd5nMPxFBsoPqNTJ2CS9emJuibT2U7w3HPRXp+IiKpLGkPmMgD/hO70M9na1hvAQAD9AAwXke2UUpluFS2t13x1NJntrYxtpRyTRUTOBHAmALRv3x61tbX5dg3EkiVLQr9GnFq23Bg6+2tu3+vEie0BbAkAmDFjNmprfyj7ul7v64QJNeuXx479DLNnL3fdb+nS3tD9x4AxY8ajTp0FZZfRqylTmgLY0SjLGNTWLgVQk7VfVJ+jSv/MxoX3NRy8r+HhvQ1HEu5r6kKmUmougGsdm0eKyP7QHXJ2BnA6gPuiLptJKfUIgEcAoE+fPqqmpibU69XW1iLsa8Tpl1+Ahx/Wy5tuCtfvdeZMe7lt2w6oqelQ9nW93FdzaCUA2G23ndGtm/u+7drZy1tt1RtR/sjatMle79NnR2yzTe5+UX2OKv0zGxfe13DwvoaH9zYcSbivaX4mM4vVrP2Ytbqn8Vam1rEl3GW2/17mMRQis6l3p53c94nrmcxM83dGoZl9zPeifm503rxor0dERNWtYkKmJfNrdP3of1az+UwAzUSko8sxma4YE41tPwFYC6C7iLjV9rodQyEyO624jZEJxNe73BkWC3WwiTNkXnxx9nrLfH9CERERBaDSQuYu1utkx/b3rdcDXY45yLEPrCGLRgFoAsBt7picYyhcZnCrk+dTG1fHn8WOkVIL9YSPM2SOH5+9vskm0V6fiIiqS+pCpoj8QURyyi0i+0AP6g4Azikpraf58A8RaW0c0xXAuQBWAnjCccxD1uuN1pBGmWN2hJ71Zx6Al0r8Nsgns7ncS8iMMsAtd/Tx6dQp/75x9oAnIiKKUiI6/ojIHwH80VrN9NbYVUQGW8vzlVKXWst3A9hMREZBzxIE6N7lmTErr1FKjTLPr5QaJSJ3A7gYwHgReRF6Wsmjobv6nu+Y7QfQU0r+GXrA9S9EZBiADaxj6gI4Qym1qNTvmfxJcnO5ea399su/H5AdkGfPDqc8RERESZCUmsztAAywvg6wtnU3tpkz6zwN4AvosVjOAHAO9DOSzwPYUyl1o9sFlFKXQM/qMxt6aKGTAHwL4FDnvOXW/grAsdDBdA2A86FD50jrOpy3PEJJbi73MqVkxssv28uXXBJOefy65RZ7+Q9/iK8cRERUWRJRk2nNQT7Q476PA3i8xOsMBjDYx/5rANxjfVGMktxc7nUgdiC3aT0JehtzZZlDLBEREZUjKTWZRAV5aS5PQ01mEpmPGXiZepKIiMgLhkxKBS/N5Ul4JrPQGJlJZd5PhkwiIgoKQyalgt/m8rjGyUxjTSZDJhERhYEhk1LBb3N5lM9kpr253AyZUY/dSURElYshk1LBS01mEprL0xIyzXnMzfv5wQfAypXRl4eIiCoPQyalQpKHMEpjc/kFF9jLzudIH3ww2rIQEVFlCiRkisgWInKRiPxVRDgjMgXOS3N5gwb2cpS1cX46/jz0kL28557hlCcf8/6ccoq97Azt338fTXmIiKiy+QqZInKtiMwSkTbGtn2hB0e/E8CDAD4XkQ2CLSZVOy/N5U2b2stLl4ZbHpOf5vLOne3l5s3DKU8+rVrZyw0b2ssrVmTvZ4ZRIiKiUvmtyTwIwA9KqQXGtlsAKADXQc/33Q3A34IpHpHmpbk8rpC5erW9XCxkxtmTO18YbtQoez8z0BMREZXKb8jsCmB9Y5qIdAKwA4AHlVI3KqXOA/A+7HnIiQLhpbncGTKjCktm07xZQ+gmiSFz112z91u8OJryEBFRZfMbMlsDMGsxd4OuxXzd2DYOwCZllosoi9dxMjMhT6nopnD0EzLjnF0nX42rM7Rz/nIiIgqC35A5D0AnY70fgNUAPjO2NSjhvEQFeWkuB+JpMk97TSYA9O9vLxf7HoiIiLzwGwa/BHCYiPQSkU0BHA3gI6WUWWfUFcCsYIpHpHlpLgfiCZlPP20v169feF8zZA4fHk558ikUMrt2dd+PiIioVH5D5u0AWgL4CsAEa/muzJsiUhe6CX1sUAUkArw1lwPxhMwZM+zl554rvG9cnWrWrbOvLZI71JIZOletiq5cRERUuXyFTKXUhwD6A3gFwMsAjlRKvWXs0hfATOs9osB4HYuycWN72Tk0TxTatSv8vnP8TvM5yTAVG2Zp8mR7+Z57wi8PERFVPt/zkyil3gbwdp73PgSwfbmFInLy+txjXLP+ZFxySeH3naFy8eLsKR7DMm9e/jIAwGuv2cszZ4ZfHiIiqnzsoEOp8Oqr9nKhkGk+ExlVLWHv3vZy27aF93WWadGi4Mvj5rrrCr9/xBHRlIOIiKpHyTMti0hn6J7mrr/ylVIjSz03JcPkybpTy6GHxj+sjTnVYaFpD+OoyfQz448zZEY1zJLZOcnN8ccDL72klw8+OPzyEBFR5fMdMkVkfwD3ANiiyK5FZnGmpDv0UOC774B//hNYtiw50w3+8EP+9+KoyVy71l4uNne5s0zmsWFq1Khwhx7zWdY33wy/PEREVPn8zl2+C/TA660ADAIgAEYCeBTAD9b6MAA3BFpKisV33+nXtWuzO4bErVAwmzPHXv766/DLAvgLmS1a5D82TMXGvnTWwEY9hicREVUev89kXglgBYAdlVKZ+clHKKXOAtALwI0A9gXwYnBFpCRI0nzWhxyS/71vvrGXi3XCCYqf5nJz0HPnsWEq1rnIGXY5jBEREZXLb8jcFcBrSqlfnOdQ2rXQc5tfH1D5KCGSFDJPOSXuEmTzU5NZty6wvTH+QlQ1meajDt265b7vbMaPY/gnIiKqLH5DZksA0431VQCaOvb5GMCe5RSKkifukNm6tb3cqFF85XBjBsViNZlAdtN1VDWZ5qMDd9+d+z5DJhERBc1vyJwLoLVjvYdjn/oAGoMqStwh00+TdNS8DhSfYZY/qppMk1uHKGfYZcgkIqJy+Q2ZE5EdKj8FsJ+I9AQAEekA4AgAk4IpHiVF3CHTDEaF5gc/4AB7Oaphl/w0lzv3iWPAeLeQeeCB2esMmUREVC6/IfNtAHuJSKYbwX3QtZZfiMgY6B7mGwK4N7ASUiLEHTLNMFYoZJ53nr3csWN45TH5bS6PuibT+bPbeOPcfZo3z16PavxOIiKqXH5D5r+hn7dcDQBKqY8B/AXAFOje5bMAnK2UeirIQlL8dtghvh7HSnlvLjcDaByDsfttLo+ijM6QufvuxY/JDMxORERUKl8hUym1SCn1mVJqsbHtZaVUL6VUY6XUlkqpR4IvJkVtwoTs9TVrgH//O56yOEOcSP5942iKNsO3lwHrzTJGUZNpjnlZp07h+5dx003hlYeIiKqD38HY9xSRTYrss7GIsHd5ym3hMp/TTz9FXw7AX6efOGoJ/YbMqMvoDJlERERR8PsrZwSAk4vsc5K1H1WYuHp1e+30A0T/vKMzAHsJcVHXZJrN5QyZREQUFb+/cjw0tEEAJGjobgqKl+cNw+C10w8QfXO531pMIN6aTC9N5UREREEIo16jC4DFRfeixMrXkzyukGnWZCatuXzRInt52TJvxyS1JnPAAHv5xBPDKw8REVWHog2gInKtY1ONuFeH1AWwCYBjAHxUftEoLm7jKALJCJlJay5/6CH/xyS1JnOPPYAnn9TLxe4zERFRMV6eshtoLCsANdZXPjMBXFFyiSh2+YYqSkJzebGazKiby5298L2IuibTDOmFalvNe5vvDw0iIiKvvITMftarAHgfwGAAT7rstxbArwAmKKXWubxPKZEvYMT1PJ9ZnpUrC+8bdS1h797A88/7OybqMnqtbY26XEREVNmKhkyl1AeZZRF5EsAr5jaqPPlqMgcOBK67LtKiAABGjbKX58wpvG+czztedZW3Y6Iu4/vve9svjoHsiYiocvkalEYpdUpYBaHkiGtmn3xee837vlHXxvmdtxyIvoxef56sySQioiCVPPKhiDQB0Bq6w08OpdT0Us9N8Sr0PJ5S0TebH3QQ8N//ets3bSEziprMTz7xth+fySQioiD5DpkiciKAywFsWWA3Vcq5KRkK1XytWxd9ByDzeief7H3fpM6mE3UZvV6DNZlERBQkX0FQRE4G8H/QnXw+BPAzAP46qjCFarHWro0+ZJq1fcWCXNS1hGmoyfSKz2QSEVGQ/NY2XgrgNwC7K6W+D6E8lACFajLXrPE+s01QzNrCYkEuDc3lUddkesXmciIiCpLfGX82BfACA2ZlK1aTGTU/NZlpaC43w9wLLwRbnnKwuZyIiILkN2QuAFBkpEJKu2I1mVErtSYzqc3l5lSUG28cbHnKweZyIiIKkt+Q+Tr0tJIxDctNUSgUMpNek5mG5vLdd3c/Pm5sLicioiD5DZlXAmgI4GERaRZCeSgBxo7N/17cNZmV0FzeooW9nKQxSdlcTkREQfLb8ecFAMsAnA7gOBGZBOB3l/2UUmqfMstGMbnyyvzvxR0yK6HjT8OG9vL//hdsecrBkElEREHyGzJrjOWmALbLs5/Ks51SLunN5c73ly0DmjQJvkwZZk1kKUMYZc4RdY99N3wmk4iIguSruVwpVcfjV8QjKVJUkl6T6TRkSLBlcfq//7OXvTaXO5vIV6wIrjzlMEPmxInJel6UiIjSx+8zmVTlkl6T6bRsWbBlKcRrAF++PHs9KSGzefPs9aFD4ykHERFVBoZM8iVtNZlNmwZbFpNyPBRilrOQLR0TsjpDZ1xatsxeHz8+nnIQEVFl8B0yRaSOiJwvIp+KyEIRWWO8t72IPCgiPYMtJkWpY0d7uVev7PfSVpPZunWwZTE574XXkNm9e/Z6Umoync+FhnnviIio8vn6lS0iDQD8D8C9AHoAWAzAHDNzCoBTARwfUPkoBl262Ms77pj9XhpqMo87zl4OM8A5Q6afWtbeve3lMGsynbWt5r0ppnHjYMtCRETVxW9N5mUA+gG4HkB7AI+ZbyqlfgcwEsABQRSO4mF2TDGH2wHSUZPZqJG9HGaAc94LP7WsUc1M5KxdHTSo8P4XX2wvs4c5ERGVw2/IPB7Ax0qpG5RS6+A+VNEUAJuUXTKKjTnbywknZL+XhppMs5l3wYLgy5NRTsg09/XazF4K8+dVv37xJnCOlUlEREHxGzK7Afi0yD4LALQprTiUBGbIbNMmu9dxGmoy27e3l+fMCb48GeWETDMsh3lPzZpcL83fnFqSiIiC4jdkrgDQqsg+m8B9FiBKCbO5vEEDYLvt7PU4arfmzrWXzbEc8zF7lEfZXO7nmcw4QqaXQelZk0lEREHxGzK/BLC/1QEoh4i0hH4ec3SZ5aIYmTVY9etH9/xgPoMH28tempbNIBpmbZzzXuy3n/djo2ouN8cJ9RIyOesPEREFxW/IfATAxgD+IyItzDdEpBWAwQBaA3g4iMJRPH791V6uXz87TEU5uLmbT4s9rIH4Qmbnzt6Pjaom0/x5+W0uZ8gkIqJy+J1W8lnoIHkEgHkAzgYAERkLYBaAwwE8qJR6M9hiUlTWrMkOJg0aACNH2uuPPZZ7TJS8PPdojvcYVcg0xxb1Io6Q6be5nM9kEhFROXwPxq6UOhV6LMzvAGwIPU7mHwD8COA0pdT5gZaQIvX++9nrzmcgX301urK48RIyzTI75wkPUjkzEZnfR1I7/rAmk4iIylHStJJKqcFKqe0BNAPQGUBzpdQ2SqknAi0dxc45C0zc/IbMqGoy/YbMadPs5S++CKY8bt54w17mM5lERBSlsuYuV0otV0r9opRaGlSBKF4tWmSve+nNHaVKCZmTJtnLV1wRTHnc3Hmnvfz228X3Z00mEREFpayQSZXHGZbq1gVOP91eHzAg2vI4JfWZTL8hM6n4TCYREQWlXvFdbCIy2eOuSinVo4TyUMzM4LT99vp1jz3sDj9xDGFkEim+T1TPZKYtZDZrVnwfM2TOmhVeWYiIqPL5CpnQNZ9uU0m2AtDSWv4FAOtAUsoMTpk5wM25wFeujLY8TuecU3yfNDSXx+Ggg/zt76V5nYiIKB9fIVMp1TXfeyKyKYB/AWgKPSA7pZBbj+mGDe1tK1ZEWx6nnXYqvg9DpjsvZfzuu/DLQURE1SGwZzKVUj8C+DOATgCuC+q8FC23ecLNmswRI6Itjxl6RfhMZjk23DDuEvizfDnwwQd8NpSIKK0C7fijlFoB4H8Ajg3yvBQdt+Bk1mQuWRLuNIheylMMn8l0d9llcZfAn332AWpqgKOPjrskRERUijB6l68B0MHrziJypIjcLyIfisgiEVEiMqTIMX1F5E0RWSAiy0VkvIhcKCJ5f9WLSH8RqRWRhSKyREQ+E5GCfaVFZICIjLb2X2gd39/r95ZGXoJTlM9lmsPo1PP4cEcam8vDGi7IHBuzdWv/xyu3J7AjMG8e8Mknevnll+MpAxERlSfQkCkibQH8CcDPPg67GsB5ALYDMNPDNQ4HMBLAngBeBjAIQAMA9wAYmueY8wAMA9ALwBAAjwLYCMBgEbkzzzF3Qk+h2dHafwiAbQAMs85XkaZOtZfnztWvzucwo2y+NKdFNGtUC0lDyDzjjOz1558vvzxu/JbR2TloyZJgy+PVM8/Ec10iIgqO3yGMri1wno2h5y5vCeBKH6e9CMAM6Gkp9wKQ96k/EWkBHfjWAqhRSo21tl8D4H0AR4rIMUqpocYxXQHcCWABgD5KqanW9hsAjAFwiYi8pJT6xDimL4BLAPwEYEel1G/W9jsAjANwp4i8njlXJfnrX+3l8eP1a6tW2ftEGTLnzLGX27f3dkwaQuaVVwKPPmqvZwJ90MxHG7w8z1pTk73+8MPRNrNPmQJsuy2weHH29tGjvXX6IiKi5PA7hNHAIu8vAnCjUup2rydUSq0PlVJ8EMQjoedLfyoTMK1zrBCRqwEMB3A2sms0TwXQEMBtZihUSv0mIjcDeBzAWQA+MY45y3q9KRMwrWOmisgDAK4BcAqqpIPTzjtnr0cZMn/7zV5u08bbMWbHn6Q+k+kMzF4fBfCr3Cb9H34IrixenHhibsAEdPg1a7WJiCj5/P5q65dn+zoAvwH4QSkV5mR0e1uvbiP4jQSwDEBfEWmolFrp4Zi3HPt4uc5b0CFzb1RJyHRm/y+/BA48MJprm89/mr3cC0lDTaaz6T+sjkN+azKdNtoouLJ48fHH7tuXL4+2HEREVD5fv3aUUh/k+fpQKfVNyAETADa3Xie6lG0NgCnQwbm7x2NmAVgKoLOINAEAEWkKPQzTEut9p8ys0z1L+QYqQZTNp+bzoEl7JvPJJ+1lvyHRuX8YHWycowB4DZl3320vmzXJREREfhT9tSMiDawe1sNFpH6R/YaLyKeF9itTZlahhXnez2xvVcIxLR2vfq5RVbxM7RiUpNZkrl0L/Oc/2etJU2otptkLfdGi4MpDRETVxUtz+QkAdgBwqFIq769spdQqq2PMmwCOh+6ZXbVE5EwAZwJA+/btUVtbG+r1lixZEsg1unTZEdOmNQUA1Ku3DrW1I613atbvs2LFYtTWjiv7Wl588UU7AFsBABYunIva2uJT0uhwVQNAh8wRI2pLDsb57uuKFXWgBzjQhg9HCfe/Zv3S7NnfobY22N4/kyc3BbAjAH1PvJZv2rS20AMxAFOmzENt7beBlivDeW8XLaoHYPe8+7//fm1JTf7VJqj/Cygb72t4eG/DkYj7qpQq+AXgdQCTiu1n7D8BwBte93ccWwM9N/qQPO+Psd7fIc/731jvb2lsm2dt2yDPMUus95tY602t9cV59m9rvT/H6/e1ww47qLCNGDEikPOcdppSuvFWqdtvt7dntgFK9e8fyKU8GTzYvu6JJ3o/rm5d+7hVq0q/fr77+vvv2fcE8H/uY46xj3366dLLmE/fvqWV75137GP22Sf4cmU47+1ZZ+XeU/Nr4cLwylJJgvq/gLLxvoaH9zYc5dxXAGNVCTnO+eWlXmB7ALU+cutI6DEvwzDBes15HlJE6gHoBj0Y/GSPx3SEDpUzlFLLAEAptRR6vM5m1vtOm1mvOc94VgJzUHBzGkKzCfqACGemnz/fXm7WzPtxYTeZB3FO8xnTMMpoDv/kR/Pm9rJbT++wPPxw4ffZdE9ElC5eQmZbAH5+Xc0BsEFpxSnqfevVrW/zngCaABil7J7lxY45yLFPOcdUhHwz7Jx1lr0c5fOHPxvD+nfvnn8/pzSETLOMYfSeLnUWobhCZrHOU1GWhYiIyuclZC4H4KMOCc0ArCi6V2leBDAfwDEi0iezUUQaAbjRWn3IccwTAFYCOM8amD1zTGsAV1mrzjqUzPo/rP0yx3QFcK51vifK+UaSKl/INANAlCHTvJbX3uVA+GNlOkPmEUf4P8fGG9vL5kxLQZk2rbTjzJAZ5Yw/9Yt0FwxzzFMiIgqel44/PwPoU3QvWx8A073uLCJ/BPBHazUz5/muIjLYWp6vlLoUAJRSi0TkDOiwWSsiQ6Fn8jkMeqiiFwE8Z55fKTVFRC4D8C8AY0XkOQCroAd27wzgLmXM9mMdM0pE7gZwMYDxIvIi9NSVRwNoA+B8VYGz/QDJDpl+On2YITOMmkxnLWGxgOSmQwd7OUlDBbVoYS9H2URdbEB6c6QBIiJKPi+/tmuhQ1/RoCkiOwDoiwJTQ7rYDsAA6yvztF93Y9uR5s5KqVegp58cCeAIAOcDWA0dCI+xHliF45j7oYPotwBOgu71PRvAyZkA63LMJdCz+sy29j/JOv5QpdQgH99fqiQtZJY6DM8sY4TTMWOCK0+GM7hecon/czRtai8vXVpeeYLkbC4PYwxPN8VqTe+5J5pyEBFRMLz82h4E3Zv6BRHZMt9OIrIFgBeg5xV/0GsBlFIDlVJS4KuryzEfK6UOVkq1Vko1Vkpto5S6RymVN/4opYYppfZSSjVXSjVVSu2olHoy3/7WMYOt/Zpax+2llHrd6/eWRpUSMk0nnRRMWUzOkNnHT12/xezIFEaztFmbu9VW3o+rV88OwOvWRddkXmza0KFDC79PRETJUvTXtlJqAoAbAHQB8IWIDBGRU0Vkf+vrFBEZAuALAF0BXG8dQymUL2Say3GFzFKnXgyjydcMmdtuW9o5wq7J7N3bXn7gAX/Hmk3mC/NNSxCwww+P5jpERBQNT3OXK6VuEJE10HN1HwfgWMcuAt1k/Q+l1C3BFpGi5KUmM8ypGp2CqMkMg3kPSnkeEwi/JtP8Y8DP8E8A0LixvbwirG58DuYfA1tsAfzwQzTXJSKicHj+ta2Uuhl6rMl/Qj9z+YP1NQK6prMnA2b65QuZv/xiL98S4U+51I4/YTNDZrEOK/mYwS+MmsxyaoHNKTyj6HAzZgzw0kv2+gZhDYJGRESR8fXrUSk1Dbo2kypUvpBpzkyVtmcyw2Dep1JrMs3m8rBrMssJmVHUZB58cPb6ySfrIZhmzAj/2kREFI4E/dqmJMgXMs8+O/qyAME8k7ll3u5qpQu6uTyMmsxyQqY5JmkUNZnmzE6Abq6fMCH3WdIoH9UgIqLyMGRSlnwhc5dd7OWdd46uPKXWZB56qPtyUIIImazJzG/lSqBJE6CmJnv76NHRl4WIiErDkElZ8oXMuEJHqc9k7rGH+zmCEkTIbNwYENHLK1YEX86JE+3lpNdkOn32mX41Z0UCSp+PnYiIoseQSVm++speTkLILLUm0yx7qXN4FxJEyBQJbxgj5zzffkNm2HO/F3PCCfq1eXOgXTt7+/33R18WIiIqDUMmreesSTNDXRJCpp+gFPa4nkH0LgfCG8ZoypTsdb8h0/yeTj657OL4ttlm9vLcufay2QGNiIiSjSGT1lu+PHvdnDEmCSEzSTWZQfQuB8Lr/OMsk9+QaYboOOZVNz97TlGObkBERKVjyKT1nGFs003tZYbMbEE0lwPhdf4555zsdb/zj5v3PQ6FQuZzz0VXDiIiKh1DJq1nhjHnYNhmR5A0dPwxa+6SHDLDqMlcvTq3Wdlv5x1nyIw6dBYKmffcE105iIiodAyZtF6+nuVAbm9jvzVjpQqiJvPDD4MrT8arr9rL5YRMc/rGZctKP4/JbZ7yTTbxd463385eDyOoF1Lonpr3jIiIkoshk9YrFDLr1MmuXVq1Kpoyldrxxzzup5+CK0/Gm2/ay+UMqxPGUEHjx+duK1Qz6EXUITMztBOQO1ZmkyaRFoWIiErEkEnrFesxHcdzmaXWZC5alL0eZs2rOa+7X2bIDCq4h/G9hhkyi537ttuy11u0CK8sREQUHIZMWq9QTSYQT8g0r+OnNs5Z6xlmZxEzKPpl3udZs8ovCxDO85Nhhszp07PXMwOxZzifD95oo/DKQkREwWHIpPWSGDLNGsmWLb0f56z1HD48mPK4Kacp2mx2v+yy8ssChDPET5g/b+fA8c6aSnMwdgBo0ya8shARUXAYMmm9YmM/xhEyFy60l/00kzpr88KsiSsnZJrDFq1YARx7bPnlcYbMUs75pz9lr7/wQunlKebf/85ed84w1Lw5sPnm9jrHySQiSgeGTFovic9kmjWZfkKm87nEMDsqldNc7jR0KDB6dHnncIawf/7T/zmc0zdeeGHJxSnq00+z1906QJ14or0cVaczIiIqD0Mmrfff/9rLX36Z+37UIVOp4EJmmPNvl1OT+e67udvMaRRL4QyZPXr4P0enTuWVwY/jj89e79Ahdx+zt/mtt0Y3hBYREZWOIZPWu+mmwu9HHTKXL7cDU6NG/sKcs7k8yNov57nLmbt8v/1ya0LLrRkNoznZLfgFxdl7vHPn3H2GDs1e/+678MpDRETBYMgkz6IOmaU+jwmEW5PpfL7T77zgTltsEez5zIHiy2EO4L7XXsGc0828efbyHnu47/P119nrQc7zTkRE4WDIJM+iDpml9iwHdA2hKcyQ6Wf8TjdffZW9/tJL5Z0vKPfeay+PGBFNE/Wvv7pv79o1e53N5UREyceQSZ5FPX95qc9jAsDWW2evBxkynecqN2Q6PfhgsOcrlTmzzty5wBtvBH8NZ1jM96iAGXgB9jAnIkoDhkzyLOqaTHMmnVatyjvX9tuXd7xp+fLs9aBDZlI0bZq9fuihwV/j+++z180OPqaePbPXwxhwnoiIglWhvx4pDGbNlnMA7TBMm2YvO0OGF+edZy+7jftZqmXLstfzBSOv+vcv7/iwRDFHuLPWNl94dD6nGvVc6kRE5B9DJnlm9jCeMSP865mBo5Qe192728tB9i53DjFUbk1muSE1LM6azDC0b5+97rWG0m0sTSIiShaGTPJs003t5S++CP965nN3pQQ5c8gj5zN95Xj++ez1ckOmW7CaMKH087VuXfqxpihqMs1HMADvHXoYMomIko8hkzzr1s1eNjvlhMUMX6UEOWcTa1g9ko8+urzjmzXL3TZiROnnMwdSv+WW0s8TRU2msxNVvppMZ8cvhkwiouRjyCTPGje2l6Po+GMGjlLGjnQ2QwfVI9kM2wCw557lnc8tCJYTiM0A5pyD3A9nTWYYHZycITPf9+0cDJ4hk4go+Rgyab22be3l007LfT/q3uXmNUoJOM6QGdRzmWb4veCC8s/nDK0A8MgjpZ/PDGDlTHkZ5Jzs+Tg78BRqoj/1VHuZIZOIKPkYMmm9Ll3s5TPPzH0/6pA5cKC9XMoML2GFzFtvtZfDGr7Ibe54r8zvs5yg6Lx/QQ8btGRJvZzpIh9/PP/+5vfinAGIiIiShyGT1isWTszmcudYkWEbPNj/Mc6QFETt1+zZ+iuj3CkgC3HOBOSV+X1GURtZqiuv3AY//mivH3ss0Lt3/v3NkGseR0REycSQSeuZz8e5jStpbgtyBh0vSnme0lnLGERN5qxZ2ethhkxznE8/fvvNXk5qyFQK+Oab7LlCnQOzO5nPZYYx+xAREQWLIZPWM0OYW8g0A1XU0/qVMpZkGM3lznOGOdtPKbXFM2dmryc1ZLo1vc+bV/gY5/zlRESUbAyZtJ5ZO+nWYcQMVFFP61eNIXPcOP/Pojpn0KlXr7wyHHecvbzlluWdy+T2R0qxP1ycP7+whqQiIqJgMGTSemagMZ+/zIizJjMpvcud53z33fLPWcjtt/vb/+abs9fLnU3oppvs5e+/Dy7Yuf2RctddhY9x1mQGOYsTEREFjyGTAOihZMxn+dq0yd0nbc3l222XvR5GyPz88/LPWcgNN5R+bLt25V/fWRO62WblnxNw//zsvnvhY/bbL3udwxgRESUbQyYB0DP4ZGqpWrZ0b2ZNW8jcdtvs9SBCprMmL8zm8oxp00o77qmnyr+2s2PTTz+Vf07A/fNTbExPkew/fhgyiYiSjSGTAGQ/j5mvs0jaQiYA1NTYy0GEzJdfLv8cfvXp460T0FVXZa87Q3Yp3P7YeOKJ7FrvUrj9LNw6mzmZn002lxMRJRtDJgHInnklX2cRs9ZOqWg7XpQaMs3asSBCiTlAPFD+M49ezJ8PvPZa4X3Gjs2dntLtkQe/3IZoOvVU4LDDyjvv/ffnbvMyO5EZMlmTSUSUbAyZBMBbyBTJDlVR9jBPSsh0iuoeFCv7m2/mbitnSsmMfOOAfvRReX9kvP127rZCU0pmMGQSEaUHQyYB8BYygexgEWWTeVJDZlS1ucWuE1Y5Cn0WygnY33yTu83LwPYMmURE6cGQSQCKz/bjZu7ccMriptQONmGHzE6dgjnPH/9Y+P1CIXLdutKnoCymUPArddan1auBZcuyt221lbdjx4+3l52PLhARUbIwZBIA7zWZpnvvDaUornr0KO24sEPmqacGc56HHwbOPTf/+4VC5nHH5XZI2nnnYMpVKGSanxk/nnsud1spteKvvsoB2YmIkowhkwCUFjKjbK78179KO84MmX4HNvfi738P5jzt2wODBrk/Wwno2r81a3QHoK+/zn7PLbRde20w5QojZC5alLutZ8/SzuWsESUiouRgyCQAwIQJ9rLXptcwa5Hmz89e79y5/HN+913pTbwZ5gDv3bp566ziR77ho5YvB+64Azj8cF2G6dP19nzzfbdvH0x5Cj2mUGrIdHu+9uyzSzuXW2AlIqJkYMgkAMDll/s/JsyQ+eOP2eulPpM5a1b2ut+5wJ3atrWX77ijvHO5cU6dmPG3v9njYK5bZ9dUXn+9+/7lzlnuxd//Dowe7f84t5B50EHejnWGZ4ZMIqLkYsgkAMCSJf6PcZvfPCiNGmWvlxoynTWX5Tavmsd36FDeudx07+6tVu/JJ/XrL7+4v++lp3a5nnhCj5dZbu2wH//4R/b64sXRXZuIiPxhyCQAwHnn+T9mk02CL0dGvmZjv5wdSoIMmWGF7EGDckO2m0ID4kcRMgFgzpzcRxuKcRu+yKsWLbLX//e/0s9FREThYsgkAEDz5vbyhRfm38/sUDJnTmjFyRmDsdTmX2fTbLnTIZohM+jnMTPq1AF++KH4fmPGRDMY/F13FX7f7yMIDzxQelmc4dk5lSYRESUHQyYByA4KhQKdOS7kzTeHVx5nDWTLlqWdp1Wr7PVyh10y5xAPK2QCQJcuwD77FN6n0DBFQQ6Uf/HFhd/386jFZ5/lbjv5ZO/HlzqUFRERRY8hkwAA111nLxfqTBHVNIpmSDJ7dPt1443Z6888U/q5gGhqMjO6dCm+T745zYOejemaa/K/96c/5X821Ondd3O3Pfig93Lsuqv3fYmIKF4MmZTj0Ufzv7diRTRlMENSOc8Xbr55+WUxmTW+YXZ8AoDTTiv92FKn4cxn4EDdKcnN1Km693upvDx/aurXr/RrERFRdBgyKUehoYmcTelhDWMUVMgM0po1dsiuUyf8msy+fUs7rm5dYJttgi1LnTrA228DRx3l/v6LL3o7j1stud9A7Az3Eyf6O56IiKLBkEk5Cg0X5AyVYQ1fYw70nZSQac6bvW5d8LWFQRk3LpyybbYZ8H//l//9YoOzL1oE3Hln+eVwTr/55z+Xf04iIgoeQyblSELINGsyyx1Y/Pzzs9fNzjt+7L57eeWIirOzU5CaNs3/XrGOSvfcE0wZamqy17/9NpjzEhFRsBgyKUehWrATTsheX7UqnDIE2Vx++OHZ64Vq4/L55JPSw2k53nvP/zFhzsRUyMiRwE476fFTH3kk9/0FC3K3DR3q/zpuf3RE1SGNiIi8Y8ikHIVqMtu0yV4PK2SuXGkvlxuanM/wLVzo/xxHHFFeGUq1zz7AFVfEc+18CnUMGzMG+Pln4K9/BfbeG3jlFfs9t5/joYf6v75byBw2zP95iIgoXAyZlKPYFI7mWJlhhcz//tde/uCD8s7lDJmlhFbnHOhRuuIKPU+4V0HNlpTPKacArVsX32/ECD280amn6uc533gjdx+/PcsB98/nlVf6Pw8REYWLIZNyFOs0Yj7zN2NGOGV46qngzuUMmTfdVP45O3Ys/xxetWwJ3Habt3332CP8stWtqwOkV088Afz4IzB5cu57pc5J7/T998Gch4iIgsOQSTmK/eI3x0sMa2rJc86xl3v2LO9cztqy5cu9D7mTz8CB5R0fhlmzyq/19Srs4ZtKMXdu3CUgIiITQyblKBYyzdBmPjsZJPMaJ51U3rncBk3/y1/KO2eLFuUdH4YOHaIbVimIgejbtw92ZP+ffgr0dEREVCaGTMpRLGSaz/yFFTLNoZHq1y/vXPkC0YQJ5Z03SS66KNrrFRrKyKvHHx9T/kkMSR23lIioWjFkUo6rry78fhQh0+xQ1KBBeefKFzK9zrdtzleeEdcwQU6nnKKbyIMY5NyPIGpymzYtfYJ1t44+DJlERMnCkEkAsoPjccd53zcNNZn5jvfa5PvOO+VdPyy77ALcdx+w557BdaDxKu5ZmK6/Pr5rjx4NbL45cNhh4U1GQERUCRgyCUB282exUGfWLIYVMqdP916eUnkJCKtW1UnMtIW9etnL226rB4hv3jy+8sSpfn1g+PDsbZ98AixZEs711q4F7r9f16DuvLOeL33YMD0L1JVXxjvEFRFRUjFkEpYt0/NKZxSr4YuiJtMcUzGskLnnnsDs2YX3mT07dyDHJk1KG0S8XEOH6mbqli2BZ5+N/vpOEyeWfuxDD5V//W7dstcvugjYeutwaheHDgUuuAC49dbs7aNH621nnRX8NYmI0o4hkzBvHrBmjV7u1Kl4pw4zZIYxGLvzGcggnn/MN+TQ5psXPs5tusJPPwWaNSu7SL5tvbV+jnTmTGDLLaO/vtNmmwFffGGv+5kzPYjmdmfIBHQN+Kuvln9up8svL/z+a68Ff00iorRLbcgUkakiovJ8udZPiUhfEXlTRBaIyHIRGS8iF4pI3l95ItJfRGpFZKGILBGRz0RkQHjfWfTMoOhlBpawazKdY28G0aHDbSpCQNfgZgK2m7Vrcy++zTbll6dUTZsG07M7KNttB/z2G/Dww3ru8qj175+7za2jVrmS0tGLiChN8vzqTY2FAO512Z7zZJaIHA7gJQArADwHYAGAQwHcA2A3ADkjJ4rIeQDuB/ArgCEAVgE4EsBgEdlGKXVpIN9FzPz25A47ZC5dGvw5C3n2WeDEE93fW7kytX+HRaZVKz1XuVfNmwNHHZVdC1oqtyk0w+hl7iVkrlwZ/pSeRERpkvaQ+btSamCxnUSkBYBHAawFUKOUGmttvwbA+wCOFJFjlFJDjWO6ArgTOoz2UUpNtbbfAGAMgEtE5CWl1CeBfkcxSFrIDKMJvlBIOOkkYKONgH79gMWL9TOPADBuHHDRRdsFX5gKds89+cfs/Okn3Vln553te1yu99/P3RZGT3svHXvmz9ePmxARkVYt1TRHAtgQwNBMwAQApdQKAJlRIc92HHMqgIYABmUCpnXMbwButlYr4nH/pIXMMDpuFKuJ2ndfPed3+/bAf/+r59nu0wdYtSr7SYqttgq+bJXkb38DPvsMeOut3Pe6dwfOOAPo3Tu46/32W+62Tz/VvcGjtiLYCYyIiFIv7TWZDUXkBACbAFgKYDyAkUop56+Yva3Xt13OMRLAMgB9RaShUmqlh2PecuyTamao8xIywx7CKIyaTC8yc18fcQSw227u+5xxRnTlSSMRYKed9PKMGUDnznp53LjoyjBoELBwIfDUU9FdE2DIJCJySnvI7ADgace2KSJyilLqA2Nbpg9xzqArSqk1IjIFwNYAugP43sMxs0RkKYDOItJEKRVCV4PoJK0mM4yQ6bezzscfu28vd87zatKpU3wdZp5+miGTiChuaQ6ZTwD4EMC3ABZDB8TzAJwJ4C0R2VUp9ZW1b+YJsIV5zpXZ3srY5uWYptZ+OSFTRM60yoL27dujtra28HdTpiVLlpR8jbFj2wDobZ1nAWprxxfcf9KktgD0yOC//DIftbXflHRdL+UBgB9++AG1tUUGtCyiZUugX7+t8N13LTBnjocu9C5at16FSZNGYdKksopClnI+sxmtW/fFb7+5/2X07rsj0aCByxhUPsyd2xBXX90LQPFR70eN+hyLFy8qul/YgrivlIv3NTy8t+FIxH1VSlXUF3RnHQXgZWPbRGvbpnmO+dh6f1dj2yprW708x8y03u9YrEw77LCDCtuIESNKPvbVV5XSdU5KHXpo8f3feMPe/8ADS76sp/IASn39dXDnXrcu+9x+vo46KrhyUHmf2YzPP8//82rTRqkffyzv/Icc4v3z8d57ZX87gQjivlIu3tfw8N6Go5z7CmCsCiCTVWLHn4et1z2NbZnayHx9WjPbfy/hmHw1namRtOZyZ7OjOZ1iucoZ3ibfWJsUn+23zz980oIFwKmnlnd+c+apYvbdV/euJyIirRJD5jzr1RyyeoL12tO5s4jUA9ANwBoAkz0e09E6/wyV8ucxgfJCZhBjHTrNn28vn3568OcvVRhD41D5Ck0g8NVX+d8Lw8UXc+B2IqKMSvy1uYv1agbGzGh6B7rsvyeAJgBGKbtnebFjDnLsk2rlhMzffw9+8HRzxp8OHYI9dzkYMpOp0IxNCxfqcVBLnQWolJ95ofIQEVWTVP7aFJEtRSRncj1rAPVB1uoQ460XAcwHcIyI9DH2bwTgRmv1IcfpngCwEsB51nkzx7QGcJW1+jAqgBky69cvvr+z5ijo6QTNkNm+fbDnLkeXLnGXgNwUG1f16aeBu+4q7dxuc9dnPPOM+/YwHiEhIkqjVIZMAEcDmC0ib4jIgyJym4i8CD380KYA3oTuAAQAUEotAnAGgLoAakXkMRG5HcCXAHaFDqHPmRdQSk0BcBmANgDGisgDInIP9FicPQDcpSpgth/Af03mRhtlrzdrFmx55s2zl9u1C/bcAPCPf5R23MUXB1sOCoaXwfuvvTa4640bB3z9NXDsse7vcygjIiItrSFzBIDXocPecQAuBrAXgI8ADADQXymVNdqiUuoVa5+RAI4AcD6A1daxx1i9qeA45n4Ah0EPk3QS9JBEswGcrCpk3nIAeM6I115qMtu0yV73cowfi4xRYIKaftB0ww3+j7n6aj1HNyWP1xmi3nzT+/OSy5cDU6a4v9epk90Z7bbbct9nyCQi0lLZX1bpgdY/KLpj7nEfAzjY5zHDAAzze600+egje/nrr4vvLwJ062b/Eg568PTFi+3l5sWHJ/StTh39HOkHHwAbbqjn1D7mmMLHlBJMKRpeRww45BBgxAigpqbwflOmAJtumr+pvK4x0+i55wKXX579ft++urZzww29lYuIqFKltSaTQnLYYd7269rVXg56rvGwQyYANGkCHHSQnp/86KML73vDDeUNfUThuu467/v261d8nzPOKPwsphkymzYFXn45+/2ffwb+/nfvZSIiqlQMmYS2be3l/v29HWM+u5m2mkyvRo0CrrrqOwaGhOvRA/jEx9PR+aYMBfQfTMOHFz7e+QfHoYfm7jN4sPfyEBFVKoZMyhpkvHFjb8eYtTmVEDKvvDJ32667AvvtNzdryCZKpl12Kb5Pxu6753/v3/8ufrz52XdbJyIijSGTsjoqFBrY2mT+Yn366eDKolQ8IfOmm4APP7TXN944mutScI46yvu++YYf+vzzwseddlq8tetERGnCkElYvtxe9lqT+fvv9vJLLwVblrVr9XLDht6GVAqCiK7h+vBD4JJLgLffjua6FJxnngE+/dTbvscfD0yenL1t1izgiSfyH3PSScBjj5VePiKiasOQWeXWrcsePNprTeaZZ2avL1kSTHkWLLCXnUMlRWH33YE77wS22ir6a1N56tYFdt7Z+/6XXaZfx40Dzjknd/xXp+OOK71sRETViCGzypkBs2FD772oncOzBDXLya+/2stxhEyqHm++qUdJ6NMHeMg535eL/ffP/97ttwdWLCKiisGQWeVKaSoHcge1DmoAarMmc4MNgjknkZsVK4Bp07zte9RRhf8AO/zwYMpERFRJGDKrXCmdfoDckGmG1XJMn24vsyaTSnHuucGfs3Pnwu/37Jm9vueewZeBiChtGDKrXKk1mebYms7zlMOc4rJjx2DOSdXl/vuBMWMKj4fpR9u2ubP6uHn+eXt55Eg98LuzcxERUTVhyKxy5nBBzZp5P65Pn+z1oJrLzdmDnEGWyAsR/fns2xd48MHyz/f990C7dsX3cw57VVtbfLpSIqJKxpBZ5RYutJdbtvR+XGbIn4ygajLNsLrffsGck6rXmWe6D7Tvh9c/dtzmKh8zprxrExGlGUNmlVu0yF5u0cLfsWbzelAh0zyPn2dEidzUrQvcfHM4z2k6eantJCKqJgyZVc4c39LvTCZmCAyqubzUjkhEhdx/f2nHbbed9339PG5CRFQNGDKr3LJl9nKTJv6ODaMmkyGTwiACPPxw7vZi04f6mapSBNhiC3/lIiKqZAyZVa6ckGnub56nHAyZFJa//lVPHWk6/nhg1SqgU6fs7ZttBhxyCHDRRf6u8coruduc1yQiqhYMmVWunJBpNq+bvdTLYYZMP0MqEXnRoUP2ev36+uvzz7O3f/st8Prr/v/QcesktNFGekgjIqJqw5BZ5ZYutZebNvV3rNlRyOxAVA7WZFLYzGGFBgzQr87ZfOrVK+3c+Wap2nff0s5HRJRmDJlVrpyaTHPIo2uvBdatK788ZuhlyKQw3HsvcN11wGuvAT166G3Oz26hKSRLYY7/SkRULRgyq1w5IbNbt+z1114rryyTJmWvl1qbRFRI+/bAwIHAoYfa29q1A7beWi+XW+uYrzaTiKjaMGRWuYcespf9hswDDshed4ZEv955p7zjiUolAgwfDgwZkj09ZCnMqVFN48aVd14iorRhyKT16vj8NAQ9LuCUKcGej8iP9u11b/PWrcs7zz77AJdemrvdnIp13TpgxAjgl1/KuxYRUZKxQbKKrV2bvb7//v6Odz63Vk4z4bp1wN13l348UZLk6wj3/PPAggXA9OnALbfoP9RmzPA3pSsRUVqwJrOKffdd9rrb3MvFmD11yxmQfebM0o8lSpq+fd23H300cPbZOmACesatQYOiKxcRUZQYMqvYqlXZ66V0tDEHsTZ7hvvl7H176qmln4sobsce633foKZkJSJKGobMKrZmjb28446lncMcW7OckLlyZfZ69+6ln4sobg0a6GZwL4IeLomIKCkYMquY+Uxm3bqlnSOokOmsVXU+L0qUNp06eWsKf+ON8MtCRBQHhswqZtZkljompTnrz4IFpZdl4cLsdbMnLlFanX568Vr5zz8Hrr8+mvIQEUWJIbOKBVGT2bWrvfzBB6WX5a67stcPOqj0cxElRcOGwI8/Ft9v4MDQi0JEFDmGzCpmhsxSazI7dLCXJ08uvTbTOVsQn1OjSuH1s/ztt+VdRyngqaeAN97ogMmT9XUzX99/X965iYhKwZBZxczm8lJrMtu0yV6/557Sy0NUzU46ydt+a9YATzwBPPNM9pzrzz4LDBgA3HnnFuvnZM/YaqvgyklE5BVDZhULoibTOYj0Y49l/+IjIuCTT4rv8/nnujVAqcL7PfWUHuLr+OOBt94Cfv4Z6NFDrxfCf5dEFDWGzCoWRE2mc77z2bOB//7X3zmWLCnt2kRpscsuwGmnFd+vRw89vevcufqPwJ9+yt3HPE///sAmm+hwWoxzBAciorAxZFaxIGoy69fP3faXv/g7x4MPlnZtojS5667sIb8Kad8e2G47YNNNgc6d7drNcob2mjQpdzxaIqIwMWRWsSBqMr2cu5gvvgj22kRJ1LIlMGyY9/2/+Ua/zpxpj9wwe3bp1+/dG+jZs7zpX4mI/GDIrGJmT/CgQ+Z//uN939atg702UVI1a1bacWefrV8zc56Xavp0/dw0EVEUGDKrmDl+X5cuwZ775JOL7zN9uu64MG9e9vaePYMtC1FS9OlT2hSuq1frPwofeKD8MowcWf45iIi8KPFJPKoE8+fby5ttVvp5unQBpk3L3b52LTBuHDBihB6epWNH+72FC4FevYDFi3OPe+ml0stClGQiwKhR+g+8H38EDj3U23E//QRssEEwZXjxxWDOQ0RUDGsyq9ivv9rLbduWfp7bb3fffsIJwG67AVdcAWy0UXanhYcfdg+YI0fq8ElUqerVA7bYQvcMVwq4+OJwrlMowM6ZE841iYhMDJlVzAyZ5dSS5OtNPnRodgeggw4C9t0X2HhjHTzddOpUejmI0uiOO4B33/U2lmYhzueqn3oq/74dOgCPPFLe9YiIimHIrGLmL7VyQqYIMHVq8f3+9z9g+HBgxoz8+wT9bChR0tWpA+y3nx5Lc/Xq0s6xxRb68ZfMrEGHHAK0alW4heKvfy39ekREXjBkVqnRo7PXy33eq3Pn8o7PCLqXO1Ga1KunZ/DxY/x4PTd5q1bA4MHAk0+Oxmuv6feOO67wsRw3k4jCxJBZpcaPz153zkHuF8MhUTA6dwZ++MHbvu3aAdtsY6+LAJtssgx1rP/ZBw4EDj44//Fz5wI33wyceSYwa1bJRSYicsWQWaWczWhuM/f49d135R3fo0f5ZSCqBJtvDnz2We72Cy7QHehuvFEHw6++Knye1q2BN96wB3Z36tED+Mc/gEcfBc45p/xyExGZOIRRlVq3zl4utxYzY8stgUWLgBYtSjv+66+DKQdRJXBOUpCZWhLQwdCPrbcuvs8rr/g7JxFRMazJrFLm1HIHHBDceZs3B3bd1f9xq1YBjRsHVw6itOvaFdhwQ728887RXHPRIuD884G//53PaxJR+ViTWaVWrLCXGzUK9txPPw1suqn3/QcNCqa5nqiS1K8PvPeenu+8WAceLx57DDj99ML7HH44UFurlzt0CG8MTyKqDqzJrFJmTWbQNYg9egCzZ3vf/9xzg70+UaXo3Vs3jXfrVv65Tjut+D6ZgAmUP086ERFDZpUKM2QCQPv2ek7yjh11J6MmTbLf791bv150UfDXJqLymdPOEhGVgiGzSpmzgQTdXJ7Rti0wc6YeGmXePN0s/re/AWPGAJ9/DkyaBNx1VzjXJqJct97qb//DDiveg52IKB8+k1mlzCFNwqyxENEDTNerl9ss7ue5TSIq39/+pv+933mnt/2HDdNf77wDbLedHpWiQ4dQi0hEFYQ1mVXIHAoFCG62HiJKtkaN9Fzp77zj77gDDtCPwHTsqP9wnDDB+7G//w789JP+f+e//wUuu0wff8stemilF1/0VxYiSg/WZFahJUuy14PouUpE6bHffsA11wD//Gdpx2+xBbDvvrqD38kn62GP6tTRLRaA/j9m+HA9I9HBB+ugucMOwLhx+n2zJvUvf8n9w5eIKgNDZhVasMBebtEC6N49vrIQUfREgBtu0H9gDhsGvP46MHKkv3O8955+vfRS/QUARx6ph0d7/fXc/TMB083KlUDDhv6uT0TJx5BZhX791V7u2jW2YhBRzLbYQn9ddhnQpQswfXp55yu16bt3bz3E0sqVwJVX2jWiRJRu/KdchcyazA02iK8cRJQcn38O7Lmn/iN0zpxorz1xInD55Xp58mTgiSeivT4RhYMdf6rQVVfZy87xK4moOm2wAfDtt3rYsdWrgS+/1D3RzeHOojB4sO7FTkTpx5BZhcaMsZffeiu+chBR8tStq5urt91WB89jj41+PNu6dYGlS6O9JhEFjyGzyrFXJxEVUq+ensN88WKgV6/ortu3L3D77cDPP0d3TSIKFkNmlWvTJu4SEFEaNGsGvPFGdNcbP14/p7nJJnpO9dmzgR9/jO76RFQ+hswq45zd5/bb4ykHEaXPJpvo6WHr1NE1jeeeqwd3HzGi+HEXXADcd5/+P2joUODAA71ft18/PRD8ZpsBjz3GFhiitGDv8iozdWr2+sknx1EKIkqrc88Fzj5bB03Tq68Cv/wCbL+9fqbyqKOAKVP0e3feqQddzzj6aD2L0FNP6ekqhwwBHn3U2/XPOEN/ffwx8NBDunb1ww/17EFElCwMmVXm2Wft5V12yf1FQURUjNv/G4cdlr3+6qt6JqBevfQg7U6tWunaTQDYaSfvITNjt93s5V69dFitU0fPYnTwwcAnn+iZhg44QG9//HHg9NN1R6Y//1nXim67ra4VXbECaNzY3/WJqDiGzCpz99328qefxlcOIqps22yjn6X0olEjPWTSdtuVfr0vv9SvhxxSeL9nn9Vf9erpWY7+9Cc9LugJJwBPPmnPPrRiBYd4IyoX67Gq2H77xV0CIiJt222jvd6aNfq50szA80OG6Gb+Jk30a9Om+rGAYcOA0aOjLRtRpWDIrCKzZ2evX3ZZPOUgInKTtHF7H35YPwaw885Av341+Pbb7PeXLWMnJKJC2FxeRTp2zF7v1y+echARuTnwQD0e5803606KDRoANTXAKafEXTLNHCd0m22Ar7/Wy7/+CkyYAEybBvzxj7r5n4gYMqvGDz/kbqvHnz4RJUyzZjpkmrp0Af7zH91xZ+VK3UnoP//JPbZdO2Du3GjKmQmYgJ4ZyWnnnYHPPtPLV1+tnxVduRLYay/v11i2TIfX7bYDRMoqLlEs2FxehIh0FpH/E5FfRGSliEwVkXtFpHXcZfPj/fez1zncBxGlRb9+enzMXXbRIW3IEGDePD00Ukbjxrr2c/jw2IqZJRMwAeDGG4Fdd9W1siLAQQfpdRH9LOo11wCdOgEXXaT3X7tWD8vUtCnwhz/o3vEnnqhrTE3r1gEffaQ7KWW8/bYOtTNmhP4tEhXFuqwCRKQHgFEA2gF4FcAPAHYC8DcAB4rIbkqpXwucIjGGDs1eN/9zJiJKm7Zt9XSXu+6qn4vceWfdOrP33npdKeD114GxY3UnxxUrgOnTgdNO08cffDDw5pvxlP3tt+3l8eP1FwDce6/+cjNkiP568kld2zt3ru6YlLHllsD339vrN90EPPEEsNVWwDPPAAMGAJtvrpvyOXQdRYUhs7AHoQPmBUqp+zMbReRuABcBuAnAWTGVzbNrr9V/FZt22imeshARBUVE9xDP996hh+ov08EH6+btevX0MG6tWwM9ewKrVukm+E6dgEWL9IDxAPDee3rMzWuuCfd78WrAAPftZsDMMJ9lve8+e3nbbfWg+l9/DXzxha4NBYBbb9Wh/brrgL//XY9FOnGiDqevvqqnIf79d2DwYD0gfmaAfaWAJUuA5s3dyzZnjj5u8819frOUegyZeVi1mPsDmArgAcfb1wE4E8CJInKJUmppxMXzZPFi/dftzJnZ2084gXOWE1F16tDBXt51V3u5USO7lhPQzfKtW+vOR/vsA+y+ey26dq1B48Z6Ot6vvgJ+/lmHsLT56ivgzDNzt19xhb1cbIzT//2v9OtvthkwaZK5pQYAcMstOozedpt+NOKOO3St79ZbAy1a6LDarp3+GbZqpYPtl1/qDmNTpwLjxula63btgPr19TOwra0H2xYu1H94LF8OTJ6sfze2aKF/P86bB6xerTtzNW2qH0OYPVt3ljWfhV22TD+W4fZ87MKFQMuW3u/B2rXAmDG6prlFCz93L10YMvPL9L1+Vym1znxDKbVYRD6GDqG7AIjtKaD33tPzBo8bt8X6KdaWFom8Tz0VTdmIiNKqffvcbV276te77sre/u23ujm+Qwcdapo0Ab75Rk97+eCDoRc1dbIDpu3KK+3lTz8F9tgjmvKUqnt3HViddtsN+PxzHWi92H9/4N13s7d16pRdQdSnj65tr1cPuPxyPW1rGjBk5pep2M/3d+ok6JDZEzGGzNdfzzSDdCi2KwD91x57KRIRBcetI2WvXsADD+gvpXQN3Wef6f+DW7bUNWcXX6wHfl+4UNe6AbrJ+YYbgHvu0R19OnXSQyPdd5/+cmsWp3i4BUwA+Phjf+dxBkwgtwVy7Fh7+d//ZsisBJmK74V53s9sb+X2poicCd2kjvbt26PW6/xqPs2YsSmAzp72/fe/x2LKlCWYMiWUolSkJUuWhPazq3a8t+HgfQ1HEPe1USNgiy3s9Rdf1K8LF9ZDgwbr0KjRuvWVAIMH51YIPPggMHNmY8yY0Rhbb70QTZuuXf/e/PkNscEGK/HFF63x7bct0KPHEixZUg/LltXD7NmNMGlSM8yf3xBduy5Fp07L8dZbHbB0aT2sWcNeQGmzdOmvqK39uuh+Sfi/gCEzJEqpRwA8AgB9+vRRNTU1oVxn8WLdRPPMM6uwaFEDrFuXOwPFsGFA//4A0CeUMlSy2tpahPWzq3a8t+HgfQ1HWu7r3nsX26PwhOxKZYdbpYBfftGPD5hjK8+dq2taW7XSz0DWrZt9zMyZ9jFK6Z79n32me9IPGKC3DRyon6vcfPMJUGpzLF2qtzdvrgP5BhsA3brpZueuXYGXX9Zjpr78sv4+331X9y9o316v//ST3rZ8uR5bdNw4/ejCtGn6ucvMo2SdOgGbbKI7dNWvr5uiP/lEf9/77283g7/zTu79+fOfda3i9Om6A9VXX+ntbdsC8+frbd27A7/9Zt+T6dP1+dZafxP066cfc2veXHeeql9fn3PcOP1+48b6mdQRI7KvLaI7c+2++waePotJ+MwyZOaXqanM9yhvZvvv4Rclv0zvySOPHBX7h4mIiNLNWXsqokOZU7t2+ivfOTp3zl7v2lV/ZXrtA8Czz+rX2tpZqKkp3vU801Hr3HP161VX5e5jdl6i+LGePL8J1mvPPO9vZr2msG8hERERUbgYMvPLVFTvLyJZ90lEmgPYDcAyAJ9GXTAiIiKipGPIzEMp9ROAdwF0BXCu4+3rATQF8HRSx8gkIiIiihOfySzsHOhpJf8lIvsA+B7AztBjaE4E8I8Yy0ZERESUWKzJLMCqzewDYDB0uLwEQA8A9wHYJS3zlhMRERFFjTWZRSilfgZwStEdiYiIiGg91mQSERERUeAYMomIiIgocAyZRERERBQ4hkwiIiIiChxDJhEREREFjiGTiIiIiALHkElEREREgWPIJCIiIqLAMWQSERERUeAYMomIiIgocAyZRERERBQ4UUrFXYaKJyLzAEwL+TJtAcwP+RrViPc1PLy34eB9DQfva3h4b8NRzn3topTasNwCMGRWCBEZq5TqE3c5Kg3va3h4b8PB+xoO3tfw8N6GIwn3lc3lRERERBQ4hkwiIiIiChxDZuV4JO4CVCje1/Dw3oaD9zUcvK/h4b0NR+z3lc9kEhEREVHgWJNJRERERIFjyCQiIiKiwDFkppiIdBaR/xORX0RkpYhMFZF7RaR13GVLAut+qDxfs/Mc01dE3hSRBSKyXETGi8iFIlK3wHX6i0itiCwUkSUi8pmIDAjvO4uGiBwpIveLyIcissi6b0OKHBPJ/RORASIy2tp/oXV8/1K/16j5ubci0rXA51iJyNAC1/F1n0SkrohcZP3clls/xzdFpG8Q33eYRGQDETldRF4WkR+t8i8UkY9E5DQRcf19x89scX7vLT+z3onIbSIyXER+Nsr/hYhcJyIb5DkmPZ9ZpRS/UvgFoAeAOQAUgFcA3ArgfWv9BwAbxF3GuL8ATAXwO4CBLl+Xuux/OIA1AJYAeBzAHda9VABeyHON86z35wN4AMA9AH62tt0Z9z0o8/59aX0fiwF8by0PKbB/JPcPwJ3W+z9b+z8A4Fdr23lx37eg7y2Artb7X+b5LB8ZxH0CIABeMP4PucP6OS6xfq6Hx33fitzTs6yy/wLgPwBuAfB/1v8BCsCLsPoh8DMb7r3lZ9bXvV0F4FPrft4K4H4AY6zvaSaAjdP8mY39BvOr5A/mO9YP+3zH9rut7Q/HXca4v6BD5lSP+7YAMBfASgB9jO2NAIyy7ukxjmO6Alhh/cPramxvDeBH65hd474PZdy/fgA2s/4jr0HhIBTJ/QPQ19r+I4DWjnP9ap2vaznfdwLvbVfr/cE+zu/7PgE41jrmYwCNjO07Wj/XuQCax33vCnzPewM4FEAdx/YOAKZb39sR/MxGcm/5mfX+fTfKs/0m63t7MM2fWTaXp5CI9ACwP3SIesDx9nUAlgI4UUSaRly0NDsSwIYAhiqlxmY2KqVWALjaWj3bccypABoCGKSUmmoc8xuAm63Vs8IqcNiUUiOUUpOU9b9LEVHdv8z6TdZ+mWOmQv9baAjgFA/ljZXPe1uKUu5T5udztfVzyxwzBsBz0D/fI0Mqb9mUUu8rpYYppdY5ts8G8LC1WmO8xc+sRyXc21JU3WcWWP95c/O89bqZsS11n1mGzHTqZ72+6/KPfjH0X3VNAOwSdcESqKGInCAiV4nI30SkX57nVva2Xt92eW8kgGUA+opIQ4/HvOXYp9JFdf+q+Z5vJCJ/tT7LfxWR3gX29XWfRKQRdO3FMgAfejkmZVZbr2uMbfzMBsPt3mbwM1u6Q63X8ca21H1m63nZiRJnc+t1Yp73J0HXdPYEMDySEiVXBwBPO7ZNEZFTlFIfGNvy3lOl1BoRmQJgawDdoZ+hK3bMLBFZCqCziDRRSi0r55tIgdDvn1Uz3wnAEqXULJcyTLJee5bxfSTZftbXeiJSC2CAUmq6sa2U+9QDQF0Ak5VSbmEhtfdWROoBOMlaNX9p8jNbpgL3NoOfWY9E5FIAzQC0BNAHwO7QAfNWY7fUfWZZk5lOLa3XhXnez2xvFX5REu0JAPtAB82mALYB8G/o50reEpFtjX1Luadej2mZ5/1KEsX9q9bP/TIA/wSwA/RzVK0B7AVgBHQT5XDHozFh/ixa5Xk/yW4F0AvAm0qpd4zt/MyWL9+95WfWv0uhH3e7EDpgvg1gf6XUPGOf1H1mGTKpYimlrreeJZqjlFqmlPpGKXUWdOeoxtC9HIkSTSk1Vyl1rVLqc6XU79bXSOjWis8AbArg9HhLmUwicgGAS6B7354Yc3EqSqF7y8+sf0qpDkopga4U+TN0beQXIvKHeEtWHobMdCpWQ5bZ/nv4RUmlzIPqexrbSrmnXo/J9xdhJYni/vFzb7CaCB+zVqP6LP+e5/3EEZHzANwH4DsA/ZRSCxy78DNbIg/31hU/s8VZlSIvQwfyDQA8Zbydus8sQ2Y6TbBe8z0TkemNlu+ZzWqXaX4wm2vy3lPruaNu0A+2T/Z4TEfr/DOq4HlMIIL7p5RaCj1uXDPrfadq/NznfJZLvE8/AVgLoLv18/JyTGKJyIXQ4w1+Ax2C3CZf4Ge2BB7vbSH8zHqglJoGHeK3FpG21ubUfWYZMtNphPW6v8tMC80B7Ab9TMynURcsJTK97s1/iO9brwe67L8ndG/9UUqplR6POcixT6WL6v7xnmdz+ywDPu+TNQTKKOif0x5ejkkqEbkcevDoL6FD0Nw8u/Iz65OPe1sIP7PebWS9rrVe0/eZdQ6cya90fIGDsRe7P1sCaOqyvSt07zgF4Cpjewvov7D9DHLbDRU8GLvje61B8cHYQ79/qICBrUu4t3+AYxBsa/s+1verAPQt9z7B28DWLeK+X0Xu5TXW9zAWQJsi+/IzG9695WfW2z3tCaCly/Y6sAdj/zjNn1mxDqSUsQZkHwWgHYBXoYcr2Bl6DM2J0P+Af42vhPESkYHQD6WPBDANevq+HgAOgf4H+SaAPymlVhnH/BF6erQVAIYCWADgMOghIF4EcJRy/IMRkfMB/Av6H95z0FOEHQmgM4C7lFKXhvU9hs26H3+0VjsAOAC69iEzJt188/uL6v6JyF0ALgYwwzpvAwBHQz+/dL5SalBZ33gE/Nxba8iXzaD/vc+w3u8Ne5y6a5RSN7pcw9d9EhGBHgD6SOjOHMOsfY+G/jdzhFLq1dK/63BZ8zAPhq71uR/uz0JPVUoNNo75I/iZLcrvveVn1hvr0YNbAHwEYAr0Z6o9dE/87gBmA9hHKfWdccwfkabPbNxJnl+lfwHYGHqYnlnWh2YagHth/OVRrV/WP9Jnof/j+R16wOB5AP4HPa6b5DluN+gA+huA5QC+BnARgLoFrnUogA+gg+xS6HlnB8R9DwK4hwOh/5rN9zU1rvsH4GRrv6XWcR8A6B/3PQvj3gI4DcDr0DN8LYGuxZgO/ctijyDvE/TYyRdZP7fl1s/xTThqnZL45eGeKgC1/MyGf2/5mfV8X3sBGAT9+MF86OcpF1rf/0DkqTFO02eWNZlEREREFDh2/CEiIiKiwDFkEhEREVHgGDKJiIiIKHAMmUREREQUOIZMIiIiIgocQyYRERERBY4hk4iIiIgCx5BJRJRCIjJYRJSIdI27LEREbhgyiajqiMgWInK/iHwjIgtFZJWI/CIib4jIaSLSMO4yhklENhaRtVZIvbnIvsrxtVJE5onI5yLymIgcJCJ1i5yjqYhcKCLvi8hc637/LiKjReQmEeke7HdIREnAGX+IqKqIyLUAroP+I/sTAGOhp75rD6AGes7gcUqpPnGV0QsRGQxgAIBuSqmpPo+9HsC10NMBzgGwsVJqTZ59M78krrde6wJoBWBr6OntGkDfw+OVUhNdjt8Feu7jTtDzIA8H8AuApgC2B9AXwDoAuyilPvfzfRBRstWLuwBERFERkaugw9LPAP6ilPrMZZ/+AC6JumxRsWodTwWwCMAQAOcAOAzAfwsdp5Qa6HKu9gDuB/AXAO+JSB+l1Fzj/S0AvAOgGYArANzlDLMi0g3AbQBalP5dEVESsbmciKqC9eziQACrARzsFjABQCn1OoADHceeLCIvichkEVkuIotE5GMROSHPtWqtpuV6InKViEyympl/FpHbRKSByzF/FJEhIjJRRJZaX+NE5AIRCfL/6oMAdAbwHICHrG1nlHIipdQcAMcAqAWwMYCrHLvcDx0eb1NK3eZWW6qUmqKUOgq6VpmIKghDJhFVi1MA1AfwklLqm0I7KqVWOjY9BKALgJEA7gUw1Fp/WkT+WeBUzwA4H8CH1jmWA/g7gH+77HsrgD8A+Aw6nD0FXQN4H4AnC5XXpzOt18HWfRgHYH8R6VLKyZRS6wDcaK0eKyICrK+h3BfACgC3eziP854TUcqxuZyIqsXu1uvwEo7tpZT6ydxg1Ua+BeAKEXlYKTXT5bgeALZWSi2wjvkHgK8AnCQiVyqlZhv7HuJyjToAnrD2H5Sv9tUrEekE4GAAE5VSo6zNgwHsAOB0ANeUeOqPAKwB0A5AVwBTYN/vcUqp30s8LxGlGGsyiahadLReZ/g90Bn+rG2rADwA/cf6PnkOvTwTMK1jlgL4D/T/vVkdi/JcYx10TSYAHOC33C5Ohe64M9jY9gyAVQBOLdZLPB+rFvJXa3VD67Xk+01ElYEhk4ioCBHZREQeEJEfRGRZZjgfAC9Zu3TKc+hYl20/W6+tHdfYQERuFZHxIrLEuMa4Itfw+j3UAXAadE/upzLbrRA8DMBGAA4p5xKZU5ZxDiKqIGwuJ6JqMQvAlvAZ1qwxHEdDh8IPAbwLYCGAtdBNwwMAuI6rmaeZONP5ZX2toYi0AjAGQDfrWk8BWGDt2wrA3/Jdw4cDoJ8jfcelaX8wgCOgn9d8ze+JRaQRgDbW6jzrdZb1WlY4JqL0YsgkomrxEYC9oZu2H/dx3MUANgBwilJqsPmGiBwLHTLLdTp0wLzeOVSQiOwKHTLLlenwc4Ax9qXTgSKysVLq5zzv57M79O+TOcaYnR9Zr31EpKVSaqHPcxJRyrG5nIiqxRPQwxcdISJbFdrRMePPptbrSy677hVQ2UK9hoh0ANAfemzMx/N8fQxdu3qqz3PXAfAPa/WZzHal1BQA7wFoBOAyD+ep6FmWiKoRQyYRVQWrhm0g9Aw1b4iI64w+InIgdK/xjKnWa41jvwOgayCDkO8a2wO4MoDznwpd0/gfpdTpbl8AToZ+nvI0r+Nyikg76OGcagBMB+CcovIC6GB7pYhcIiI5rWfW865DAexa4vdGRAnF5nIiqhpKqZutoHMdgDEiMgrZ00ruCWAzZHfYeRB6jM0XRORF6CkRe0EP2P48gKMDKNpT0LV994pIPwCTrHL0h56Jp+RrWONWZsLwY/n2U0r9KCIfQAfGgwC84TjPQGuxDuxpJXeHDu2joaeVnO845/dWGH8RwJ0A/iYi5rSS20JPTamgZ/0hogrCkElEVUUpdYOIvAA9nWI/6ADZCHoIni+hw84QY//xVvC7Ebr3dT3osS7/DOB3BBAylVK/iMge0AOy7w7dSecHq4zvlXmNfaGf9/zCw9zgj0KHzDPhCJnQwRzQwx0tBjANOhy/BOBda7ilHEqpT63pJc+Anr7yEOhOVMsA/AjgLgCPWM3rRFRBRCmONkFEREREweIzmUREREQUOIZMIiIiIgocQyYRERERBY4hk4iIiIgCx5BJRERERIFjyCQiIiKiwDFkEhEREVHgGDKJiIiIKHAMmUREREQUOIZMIiIiIgrc/wNnX9zdvNrHiwAAAABJRU5ErkJggg==\n", + "text/plain": [ + "<Figure size 720x576 with 1 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "x=np.linspace(1,2936,2936)\n", + "\n", + "plt.figure(figsize=(10, 8))\n", + "plt.plot(x, ds[\"Counts - Charge distribution\"], color=\"b\", linewidth=3)\n", + "plt.xlabel(\"Canal ADC\")\n", + "plt.ylabel(\"Cuentas\")\n", + "plt.grid()\n", + "plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "id": "b6692d1c", + "metadata": {}, + "outputs": [], + "source": [ + "from scipy.optimize import curve_fit\n", + "\n", + "# Define the Gaussian function\n", + "def Gauss(x, a, mu, sig, c):\n", + " y = a*np.exp(-(x-mu)**2/(2*sig**2)) +c\n", + " return y\n", + "\n", + "parameters=()\n", + "def fit(i,j):\n", + " xdata = np.asarray(ds[\"ADC Channel - Charge distribution\"][i:j])\n", + " ydata = np.asarray(ds[\"Counts - Charge distribution\"][i:j])\n", + " param, covariance = curve_fit(Gauss, xdata, ydata, p0=[1000, i*2, 200, 0])\n", + " fit_A=param[0]\n", + " fit_B=param[1]\n", + " fit_C=param[2]\n", + " fit_D=param[3]\n", + " print(\"mu=\", fit_B, \"\\nsigma=\", fit_C, \"\\n\")\n", + " fit_y = Gauss(xdata, fit_A, fit_B, fit_C, fit_D)\n", + " plt.plot(xdata, ydata, c=\"darkblue\", label='data')\n", + " plt.plot(xdata, fit_y, '-r', label='fit')\n", + " plt.xlabel(\"Canal ADC\")\n", + " plt.ylabel(\"Cuentas\")\n", + " plt.grid()" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "id": "f15890ec", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "mu= 19.157181333696382 \n", + "sigma= 49.24544222362323 \n", + "\n", + "mu= 479.0099909259773 \n", + "sigma= 51.310704584285276 \n", + "\n", + "mu= 939.6075158534509 \n", + "sigma= 52.7541516809037 \n", + "\n", + "mu= 1396.856500246315 \n", + "sigma= 58.13163757023065 \n", + "\n", + "mu= 1851.4849042741075 \n", + "sigma= 66.47690197710251 \n", + "\n", + "mu= 2292.507518508694 \n", + "sigma= 90.54113833966873 \n", + "\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa4AAAEYCAYAAAAEZhLyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABeAklEQVR4nO2dd3hUVfrHP2dmkkx6J0AghCJFEAVBBUWxCzbsrq4N67rqYt1Vf5ZV13VX3dW166qo2FZcdG3YQEAUkab03hNI723K+f1x74T0zCQzmZnk/TzPfe7Mveec+85Mcr/3nPOe91VaawRBEAQhXLAE2wBBEARB8AURLkEQBCGsEOESBEEQwgoRLkEQBCGsEOESBEEQwgpbsA3oCaSlpens7OxgmyEIghA2LF++vEBrnd7SORGuLiA7O5tly5YF2wxBEISwQSm1s7VzMlQoCIIghBUiXIIgCEJYIcIlCIIghBUiXIIgCEJYIcIlCIIghBUiXIIgCEJYIcIlCIIghBUiXEK3o7i4hvfe2xBsMwRBCBAiXEK348ILP+E3v/mUXbvKgm2KX1i2bB8ffbQ52GYIQsggkTOEbsf69YUAuN3dI0nq+PGzAND6jiBbIgihgfS4hG5HdbUTgLo6V5At6Tz791cG2wRBCDlEuIRuhcvlpqioBoAlS3KDbE3n+eSTrcE2QRBCDhEuoVuxdOm++tdXXPFFEC3pPGvXFnDttV8F2wxBCDlEuIRuRU5ORaP3t946n8rKuiBZ0zl+/TU/2CYIQkgiwiV0K7Ru7JDx1FPLeeKJ8Ewp43C4g22CIIQkIlxCt8fpDE8BCPehTkEIFCJcQrfC5eoeLvAACjdf8jIb+DspiHehIHiQdVxCt8LpdBOJky95BQdWpnI1SgXbKt/Zu7ecKWzkFIyFx1ezFLdbY7GE4YcRBD8jPS6hW+F0upnCBiazjZPZzEmEZ8SJfv1eYirrKSeKlfTlDNY3czwRhJ6KCJfQrXA6NVPZQDU2arEylQ2ocOxyAceyne/JZkPmGI5gNyce82awTRKEkECES+hWOJ1uDiOH7xnIT2RxOHuCbVKHsOFiGPmsoi9RE4/AjpPInduCbZYghAQiXEK3YvvWIg5mP2vJYBV9GU0uSoeXV6HbrRlCAZG4WEsG+5OzABhOXpAtE4TQQIRL6Fa8//jXxFHHIRcdx3p6EUcdrj05wTbLJ6qrHYxkPwDryKCy30DcKEaIcAkCIMIldDM8vZIJV5/EdlIAmP/6wmCa5DPV1c76z3Hvm1eQ1j+VHSSLcAmCibjDC90Gt1vTnxIAVPaAeuHKpiiIVvlOSUktWZSwnzhGH5XNoEGJfHtVKiPtpcE2TRBCAulxCd2G2lon/SjFhcLWry87SQZgYJgJ12mnfUh/SthNEjabwmq1EDF4IL0d4fU5BCFQiHAJ3YbaWhf9KGUf8djskdQSQQ4JZFMcbNN8YuvWElO4ErHZjH/RiuTe9HKVQm1tkK0ThOAjwiV0G/btqySTUvaQiFKK6Ggbu0msHz4MJ7Lqe1zGv2hNeh/jxN69QbRKEEIDES6h23DzzfPoRynWrP4ALFx4MfuIJ4PwijiRQDUJ1LKrgXDV9c409lu2B9M0QQgJRLiEbkNlpYN+lFIYnQrA6NHp7CeeDMqDbJlvZJk9xIY9LndmPwBqNotwCYIIl9BtmHZiLxKpYcJ54wCIjLSSODSLdCrB5Qqydd7Tv4FwWa1GuCprP6PH5dgdXmvSBCEQiHAJ3YaoQmOdU/zwgfXHdK8MrGgoKAiWWT5z5rgYAHJIqD8W0yuJSiJw7RXhEgQRLqHbYC0qBEBlZNQfq04w1nKxf38wTOoQiTXGeq1r7j6R+PhIABIS7ewnHr1PFiELggiX0G2IKDF7Venp9cdqktKMF2EkXHlrtlNOFP/36In1ke0TEiLZRzyWvPD5HIIQKIIuXEqpVKXUNUqpOUqpLUqpaqVUqVLqe6XU1UqpFm1USk1USn2ulCoy6/yqlJqhlLK2ca0zlFLfme1XKKV+Ukpd0Y59VyillprlS836Z3T2cwv+J7KsxHiRllZ/rNYULmeYxCt0uzW9qGA/cY2OJyREsp84bEXS4xKEoAsXcAHwCnAk8BPwFPAhMAr4N/Af1SShklLqbGAhcCwwB3gWiAT+CbzX0kWUUjcBn5jtzjKv2ReYqZR6opU6TwAzgT5m+VnAIcAnZntCCGEvNyNLNBCuxVsMp4y5bywJhkk+U1XloBcV5DURrgEDEii0JWIryA+SZYIQOoRCrMJNwFnAZ1ofyD+hlLoHWAqcB5yLIWYopRIwRMQFTNZaLzOP3wfMA85XSl2stX6vQVvZwBNAETBOa73DPP4Q8DNwu1LqQ631jw3qTARuB7YC47XWxebxx4HlwBNKqU89bQnBJ7qyhCpLFDHR0fXHdhZpqrFRszM8Fu56hGsbqY2OR0XZqE1KI7agDJxOsIXCv64gBIeg97i01vO01p80FC3z+D7gRfPt5AanzgfSgfc8omWWrwH+z3z7uyaXmQ5EAc82FBpTjB41397QpI7n/V88omXW2QE8Z7Z3VfufUOgqYqpKKY1IaHRMWRQFxBJfFx6LkCsrW+5xAVQnpmFBQ770uoSeTdCFqx0c5t7Z4NgJ5n5uC+UXAlXARKVUlJd1vmhSpjN1hCASX1NKeVRj4bJYFIXEEOcIj0XIVRV1pFPZbI4LwJFiDIHuWrqpq80ShJAiZIVLKWUDLjffNhSPYea+2X+v1toJbMcYAh3kZZ1coBLop5SKMa8dC2QCFeb5pmw290O9+jBClxBXU0ZlTGKjYw8/fAxFxNA7MjyC09bm5mFFt9jj2l5tDIE+cddnXW2WIIQUIStcwGMYjhSfa62/bHDcc2dqLTmR53hSB+okNtn7co1GKKWuU0otU0oty5ehnYBTV+civra0kSs8wEknDaBExRJbGx5Dhd+9vxyAS2Yc1+zc1mLj3zUhTHqPghAoQlK4lFK3YDhGbAAuC7I5HUJr/bLWepzWelx6k5up4H+2bCkmjUrsmRnNzhXoaNwFhUGwync+ec3wD4of3K/Zuf1OI6JGiruyS20ShFAj5ITLdDN/GlgHHK+1bpo9r2nvqCme4yUdqFPaZO/LNYQgkreriDjqiGpBuIqIIYUqtNvdQs3QopcZyd7SuwUBdkXhRpHkCo/eoyAEipASLqXUDOAZYA2GaO1rodhGc99sfsmcFxuI4cyxzcs6fYBYYI/WugpAa10J7AXizPNNOcjcyyx5iKCLS4x9cnKzc0VEY8eJqzz0eyoe4bL27tXs3GmnD6GYaJKlxyX0cEJGuJRSf8RYQLwKQ7RaCxEwz9yf1sK5Y4EY4AetdcPZ+LbqTGlSpjN1hGBRUmLsE5t3koswhthc+aEfaDeFKgAiMpoPL7/00smUWONIcsocl9CzCQnhMhcPP4axsPdErXVbd5jZQAFwsVJqXIM27MAj5tsXmtR5HagFbjIXI3vqJAP3mG9fbFLH8/5es5ynTjbwe7O919v7bEIXUWqO7iYmNTvlEa5wmOdKoZpioomOi2p2LirKRk1sInF1IlxCzyboy+/NWIEPYUTCWATc0iTCE8AOrfVMAK11mVLqWgwB+04p9R5GRIyzMNzeZwPvN6ystd6ulLoT+BewTCn1PlCHsZi5H/Bkw6gZZp0flFL/AG4DflVKzcYIK3URkALcLFEzQgdLmSFclpSkZuc8wqULm06Xhh4pVFFENCn2lv81K6ISSKoM/Z6jIASSoAsXxpwUgBWY0UqZBRgxAwHQWn+klDoOuBcjJJQd2IIhMv/SWuumDWitn1FK7QDuwFgfZsFwAPk/rfUbLV1Ua327Umo1Rg/rOsANrAAe11p/6tOnFALKAeFqeY4LwqPHlUoVRcTQ195yrOiK6ESySiQLstCzCbpwaa0fBB7sQL3FwFQf63yCEWjXlzozaSCaQmhiKS8z9slJzc6dfdWR8DroonDpccVgb6XHVR2dQKJTvAqFnk1IzHEJQmexVRg9Lmtq8x5Xv9FZxoswGSosJIYWhssBqI5NJEbXQXV1F1smCKGDCJfQLVClpTiwYEtoHiqJ6BiqsaGKQn+o0NPjao2a2CTjRWHofxZBCBQiXEK34OdvN1KKncio5kNsNpvFEIPi0O5xuZ0ukqluU7hq45KMFyJcQg9GhEvoFiRRQwnR2GzN/6StVkUx0bgKQlu4XEXFWNFtCpcz0RwKLRDPQqHnIsIldAuSqKaEaOLiIpudKy+voxQ7vyzYEgTLvMfj9eiMT2q1jCMxBQAtwiX0YES4hG5BirWW+P7NwyQBVFc7KSGaGEdoh0py5RvCdfx5h7ZeJskQLtd+yTgg9FxEuISwp7raQbyrCp2Q0OL5ujoXpdhJpKaLLfMNzwJpR0JSq2XcScZQoStPelxCz0WESwh78vKqSKIG1UKAXTggXEmEtgv5wo9+AaDO44DRAjtzaygjii1LNrdaRhC6OyJcQthjCFc1thbWcAE4HG5KiDZ6XM2DqoQMn721DICamJZ7jgA//phDITHkrtvdVWYJQsghwiWEPXk5ZUYurl4pLZ739LgicUFN6A4XeiLDO+NbSwMHL7xwEoXEkhXr6CqzBCHkEOESwp6a/cbcUER6y8I1fnwfSsx4hfXpT0KQ3pG1lGLnuBMHtlpm7NgMSlQMUVVlXWiZIIQWIlxC2KNKS4CWA+wCXHzxcOIz04w3nvQnIciQZDeVUfGMHJnWahmlFBVR8dgrQvdzNKW0tJbHH1+K2x26w7RCeCHCJYQ92hQjawspTTwkeFzlQ7jHFV1dRlVUfLvlqqITiK4Jnx7XrbfO5667FvL22+uCbYrQTRDhEsIeSzs9LoAauxnDMIR7XLG15VTHtD6/5aEuNpHYugpwu7vAqs5TXGzMK15++RdBtkToLohwCWGP1UxpYktreY4LwBFj9mRCuMeV4KigNq594aqNTcSKDmkRbkjDSPclJaHrHCOEDyJcQthjqTBS2Uektd7jqnfcCNGbvdaaRHdlm+GePFRGm+7yYZBfDODLLw8kvpw06b0gWiJ0F0S4hLAnwnRUiEhPbbVMcnYGAHUhGnGiqqKWFKpwt7KIuiHzVhmJJH/9bmOgzfILVVXO+tdr1oTm9y+EFyJcQtgTUVmOG4UlqfVhtvjeyTixUJtX3IWWeU/pngKsaFRK68OdHgrN6PG7Vm5vp2Ro4nKFx9ycELqIcAlhj62qnHKiwNL6n3N0TASl2HEXl3SdYT5QsWs/ANY2eo0ePGlPoqtD37Pw559zicTJbSzgdAyvwosu+iTIVgnhTvOse4IQZkRUlFFmiaYttwa73UopdlJCVLiq9xjCZctIb7esR7jslaE5X9eQI454m2f5H7/nRwCO4wY+/BDcbo3FotqpLQgtIz0uIeyxVZZRaYtts0x0dIQRPaMsNG/2tblGmpKoPu0LV7EZBcQeBtEzUqjkOn7i17FTySWBO1kAGBH9BaGj+EW4lFLDlVK3KqWuV0q1788rCH4ksqqC6qi2hcvT41Ih6lXo2G84LURntpxTrCEurJRgx14Vmp+lIWexjgjcrDz6QvRvfsMpbCKBaiorRbiEjuOTcCml7ldK5SqlUhocOwlYCTwBPA+sUEq1P1AvCH5Aa40qK8ER23pEdQC73UYpdizlodlL8SSRjMvK8Kp8ETEhP1SYl1fJGaxjJ0lsihuI88yzicTFqWwS4RI6ha89rinABq11wwUkfwU08ADwAjAQ+IN/zBOEtnE43MS7qrBntP2sFB8fSQl2XIUhuvapyBCuhAG9vSpeSAxRId7jOurIt5nIThYyiNS0GGwTj6KSCI5mRyMXeUHwFV+FKxtY73mjlMoEDgee11o/orW+CZgHTPOXgYLQFtXVThKpwdVGKhCA0aPTKSWaqJrKLrLMN1RxMeVEYU+I8ap8ETHsX78nwFZ1DueOnfShnL39R3HLLWOJT4llKVlMZIf0uIRO4atwJQMNH1mPxuhtfdrg2HIgq5N2CYJXVFXUkkgN7vi2g9NaLApHTDxRdVUhGePPVlZMicU70QKjx5XgqAigRZ1j//5KJrATgH7nn4zNZiE2NoIlZHEYOVSXhK7tQujjq3DlA5kN3h8POICfGhyL7EC7gtAhagpKsKLRiUntlq20xWBBQ1nozXNZS0uoiGx7ns7DrbceThEx9YknQ5GrrprLOPZQg42izIMA4+HhF/oSgRu1MTyifgihia8Cswo4Syk1Sik1BLgI+F5rXd2gTDaQ6x/zBKFt6vKNSBg6KandsjvLzD/3EPQsjKoopaYdBxMPM2YYwpVMdUj2HgGKimo4mP1sIB0iIuqPP/LRDQBEblrfWlVBaBdfhevvQCLwC7DRfP2k56RSyooxfLjMXwYKQls4TG+8tsI9eSjFDoAuDr2wT3GOcmpivBMui0VRiNF7DMXPAkZYpxHksZ6MRgFNrCOGUYcV67q1wTNOCHt8Ei6t9SLgDOAjYA5wvta6YZKdicBe85wgBJyty415lMh2vAoBkgcYa6QcBaF3s09wVlIT7Z1wORyu+ugZeRtC00HDVldDNsWspxdW64HbTGJaHBtIp3rZquAZJ4Q9Ps9Faa3naq3P01qfr7We0+TcIq31GK31bP+ZKAit8+8nFwEQ3bv1dPceTjh3NAB1eYUBtclntCbJXUltXJJXxdPTY+oD7ZbvCM1R+X6VOVjQrCOjUWinlJRo1tCb/mW7WLQoNEVXCH3EiUIIaxIxplcj09uPqq7MeTBHfoit5SovJwK3V0kkAeLiInnomTMBqN6bF0jLOsyAqr0ArKcXl1wyotG5zaSRRQkfvb8mGKYJ3YAOB9lVSvXD8DCMaum81nphR9sWBG8ZkgoUQt+D+7Vb1pJi5LpyFpUE1igfyV27iz5AvtPudR1PaKi6ffkBsqpzZFftwYmFMRccTXx8ZKNzW0jDiqZ3TWjaLoQ+PguXUuoU4J/A8HaKWjtkkSD4wJA0BYVgS20/AaM1NQkAV1FozXF98NISbgE++aGYa7ysY8swhkYtxSHWezRJK81hB8n0G9R87nELxrE+laE5zCmEPr7GKjwKY7FxEvAsoICFwCvABvP9J8BDfrVSEFrBXldJjSUSIiPbLRudEEMVEegQE65kt7EY93f3nOB1HU+2Z0tJaH0WD9kUs4NkHnro6GbntmCIbmrx3q42S+gm+DrHdTdQA4zXWnviEc7XWt8AjAIeAU4CxDlD6BJi6ioot3oXcaI+tUlJaK3jiqgw7Dn1ksO9rmOPjaKYaBZ89GugzOowBQVVDKCYnSQTGdl84OWBf51JKXb6Vu0PgnVCd8BX4ZoA/E9rndO0DW1wP0Yswz/7yT5BaJPYukoqI9pOaeIhOtpGCXZUaUlgjfIRW5nRa1Kp3idVsNutFBJDaghGz3j+Hz/Sh3J20PLw7YUXDWcLqSQV7u5iy4Tugq/ClQjsavC+Dmh611gMHNsZowTBW6IdVVTZvOtxxcTYKCEaFWKpTSLLzR5gcvvzdB7sdlvIhn1KrzHm3Xa2IlyRkVa2kEZCngiX0DF8Fa48aPTXmAcMblImAswUrYIQYOIcvvW4SrFjKw+tocLIylIqVRREteig2yJRUUaPKxSFK6XMcLpoTbgiIixsI4W4gr246iRKvOA7vgrXJhoL1RLgZKXUUAClVG/gPGCzf8wThLaJdVZSFemdcMXFRVJCNNaK0Opx2StLKbV69xk8WK0WikJ0qDCjxsjmfP+rv2nxfGSklZ0kE4GbRf/5uStNE7oJvgrXXOC4BhmQn8boXa1USv2M4VmYDjzlNwsFoQ3iXNVURcZ5VbZXrxhKiSaiMrSEK7q6jLII7z5DQ5IH9yHdUhMAizqHff9enFgYd3bLziY2m6W+N5ZQLC7xgu/4KlwvYcxfOQC01ouBC4DtGF6FucDvtNZv+tNIQWgRrYl3VVET5V1vxWaz4IxLMJJJah1g47wnuqqMKrt3cQobUh2dSLy7ClyuAFjVcWLyc9hLAjEJLS+oVkrVC1fUPnGJF3zH1yC7ZVrrn7TW5Q2OzdFaj9JaR2utR2itX/a/mYLQAjU1RGon1XbveysqORGb2wnV1e0X7gLcbo2tvIQS5X0SSQ/lEXEhGSE+viiHXSqFiIjWYxC89s11AFh272q1jCC0hq8LkI9VSrWZ3Vgp1V8pJV6FQuAx82rV2tvOftwQqxn2iZKSABjkO1u2FJNCFZuLVPuFm/DNSuP58bv//uJvszpFYsk+9tradu1P6JNMATHsWSzxCgXf8XWocD5wZTtlLjfLCUJAqTHj9OkE74fZXPFmINsQEa7aGiepVDHmhPYiqDWn0FyJUrgpdKKs5+WUklBewGZH279JenoMO0nGuW1H1xgmdCt8FS5vHgsVEDoTCILXfPzxFoYPfw2nMzSz6jbl3Rd/AuDrn713b68XuRDJglxXVEoEbuIG9PG5ricnl70qdJxNdv2wHhvuVl3hPXiEawChNcwphAeBSGsyAChvt5QQclx77Zds3FhEYWFozP+0h73WiPGXPqS313VWbq8D4JcFmwJik6+4zAzOpLSflqUpUy41vPYiK0NDhAFqNm4FWl/D1RCPcJUUh8ffmxA6tBsdXil1f5NDk5VqseNlBbKAi4HvO2+a0NV4Ev65XOHRYU631gJw56OneF3nh/XGuqfNP2/n0IBY5Rse4VKpvgvXbY+cAm+DrTR0ei2v//lzjgFm/PPcdsvuJJlYHGzfkkPS+KZxDAShdbxJa/Jgg9camGxurbEX+FOHLRKChifFel1daLlXt4bHmy5xQIbXdX5zw1HwIsS7Q+MpXxcawmVN9z5OoYfI9BTcqAMho0KAvg5j8fGp17Tvn7WTJAAq1m4GES7BB7wZKjze3E7AmL+a2eBYw+1YYCQwQGu9PBDGCoElJ8cYenv66RUsX74vyNa0j6XMuGHHmEkVveGev50KgC1Uomd0RriiIykmmsiKEj8b1XGyKSaXeGxx7bv33/fKxQB8/sJ3AbZK6G602+PSWi/wvFZKvQF81PCY0D2YM+dAlK6nnlrOU08tR+s7gmhR+1jKSqnBht2Lm6QHW3wsdVhDJuyTzjM8I2OyfHfOsFgURcQQVRE6Pa6RsZWUR/bBm0+TMmYYAPuWrgusUUK3w9cFyFdprf8XKGOE4HHuuR8H2wSfiagoo1T5GM9ZKUpVNBEhIlyOHCMnVXy278IFUGKJwV4VOsLVp66Agjjvhm5TD8qknCjxLBR8xps5rhZRSsVgRIpvcXm81lqWxIcRUVFWamtd9KIcC5p9+B6CqKuJrCyl3BqL9zNcBmWWGCKqKwJik6+sXbCe0USTnOR75AyAEkscGaHiDu9208dRxK9x3g3dxidEsYYkDomvDLBhQnfDZ3d4pdRlSqk1GC7vuzDiFDbdtvnTSCHwpKZGM4w8tvIYa3mCTEpCfj1XVU4BxbQcD68tKq3RRFUHf8WG261Jp5J84mjFU7ddSq2xRNeEhnCt+WYNkTj5aV+E13Wq0zPJdBYF0CqhO+JTj0spdSXwGuACFgG7Aaf/zRK6mrKyWu5jEXEY65yuZwkzZszj2WdPCrJlrRPjqKQA33sqFRGxpNYEv8dVV+cinUpiOrD42EN5VDwxISJcqz7+mVHAymLvf5OiuF4cVLQxcEYJ3RJfhwrvAIqBY7TW6wNgjxAEXC43FRV1TGEDszmEZKqZxlpO/9/WkBauZKop69Xf53pVEbH0r80JgEW+YQhXBTXx/TrcRm1MArFlleB0gq3DI/9+4dPn5/NbYIcXi489FCdkkOSqhIoKiPM9tYvQM/F1qHAI8IGIVveiosJBNsUMoIT5DGYBgxjJfmIcoZek0IPWmiSqSczy3hXeQ3VUHNEh8Nlqa40eV12C74uPPTjik4wXIRB70eNk4U3UDA+lSWbUk507A2GS0E3xVbiKgNpAGCIEj7KyWo5mOwCuCcewhCwsaEZWhe5UpaPORRLV1MX67kRSY48lzhl8h4DaagepVFGX2HHhKjUTUNaY3onBJJsiComhkiiv69RkZBovRLgEH/BVuD7FCPnUsZlkISSZN28Xo9iPyxrBFY+dx1KMzDUHV20PsmWtU1dcRgRunHGJvteNTsDuroPa4D6DOfIKseHGkeT74mMPw44aAkDptuBnEh5ACTtJxmbz/rbi7m/8rbm2he7fmhB6+CpcdwNRwItKKRmQ7iZceeVcRrKP6v6DmHBsNqVEs5cEBjtDN3qGw4zx54zzvceVU2XMBdXlB9ebrWyrkY7ElZLW4Tbs/Y2htpce+cYvNnWGARSzg+T6mJfeYOufSR1WajeFbu9eCD18Fa4PgCrgGiBXKbVCKTWvhe1b/5sqBJKR7Mcx7EBOqA30Yjj5QbSobVwFhui4Enzvcf2yy/CcXPdjcG+WN138HgDu1I4L1z5tPD/uWb7FLzZ1FJfTRTZF7CQZl8v7ZRTpGXHsJpF5MyUut+A9vrohTW7wOhY4rJVy4RFeXAAgEifZFOMYOxqA3/3uUDa8kM6lrAStIQRHhp2mcLkTknyumzWqP6wBV2Fwe1zpGPNsOq3jwnXO9UfDQzBhYItxALqM2pw8YnGwk2Ti4iK9rpeVFc9OkkkpD/4cnRA++BryyeLlFtz/IsEnThmmsKCJGm7Ml8TFRbKRXiRRA/mh2etyF5UY+w70uKbfPgmAij15/jTJZ9Ix1pJFdyBOoYe0vkmUWGNJdwd3LZdji9F73W9P44cffuN1vYkTM8mJSGWQNXTCVgmhTyASSQphRnKpOZc1YAAA9903ge2Ynm47dgTHqHZwm1HVdVKSz3X7jTQ82XLW7fWnST7j6XF1ZgEyQHFkIonVwY3355mjuuiuUzn4YN96kPEjD6KXqxTq6gJhmtANEeESSPUM02RnAxAfH0nmhJHGsRB1U14ydwMAEem+u5Lbe6cD4CgM7s0+nUrKiCK9n/frnlqiOCqJxJoS/xjVQarWGtkFYg4+yOe65Sl9sKBh925/myV0UzoSq9CilLpZKbVEKVWqlHI2ODdGKfW8Umqof80UAkV+fhWplXm4lAUyM+uPf7C0BjgwBBRqrF60CYCoXh1YA5VoDC8GO0L80CQnpZEJpKV1LMCuh9KoJJJrS/xjVAdxbd1OCXbSDurrc93qXrKWS/ANn4RLKRUJfA08BQzGCLTbcOZ+OzAduNRP9gkB5umnV5BNEbt0UqOQQUWuKEqwU7Fmc+uVg0gSNZQRRXS8j2lNAOLicKGwVQY30G5qXQkV8R13zPCQ644jsbbEJ28+f6K1Zvv8X9hBMn36xPpcv7a3IVyylkvwFl97XHdiZDv+M5AB/LvhSa11CbAQONUfxgmBJzs7gWyKcWY2jvk3depAI+ZciD4FJ1FNCdG43R1wYLVYqLDGYA9yOpA0RzGlMR2PmuHh13wrCdTy0/ytfrDKd77+eie9a/LYQQrp6b73Hh29M3GjcGwR4RK8w1fhuhRYrLV+SGvtpmW39+1ghl7wAqXU+UqpZ5RSi5RSZUoprZSa1U6diUqpz5VSRUqpaqXUr0qpGUqpVr0ZlVJnKKW+M4c3K5RSPymlrmjnOlcopZaa5UvN+md4+9nCAbvdRjbFZBxxcKPj9903gZ0kY90bmvMOyVRTTHSHnvAByiPjiakOrnClO0spj+t8j2s/8QBEFAXHA1S73WSbi499iZrhwRptJ5d49i2VEKiCd/j6VzYQWNJOmSLAl8fI/wNuwlgT1q6bl1LqbIxe3bHAHOBZIBL4J/BeK3VuAj4BRgGzgFeAvsBMpdQTrdR5ApgJ9DHLzwIOAT4x2+sWOCpr6EsZ2vQo9JCaGs1OkonM3RMky9omK8FNXXQ8Awb47g4PUGlPIKY6iC7YFRXE6xoqEjoe7snDPlO4LHnBWQuV4KggjroDnqg+snt3GTtJZvt3v/rZssCwcWMRWstS1WDiq3DVAEntlMkCSnxo81ZgKJAA/K6tgkqpBAwRcQGTtdZXa63vxBC9H4HzlVIXN6mTDTyBIajjtNa/11rfCowGtgK3K6UmNKkzEbjdPD9aa32r1vr3wOFmO0+Y7YY91rxcLGhU/8apNVJT7ewlEXttJVQGPyBtU+IcFR2KmuGhLj6J2Nog9rhyjdiCVYm+R7dvSqHVCHv18M0fd7qtjmDbYyQ7P+P3J3aovtOp2UkyWTr0E0rOmrWO4cNf44031gbblB6Nr8K1CjjFdNJohlIqEWN+a6m3DWqt52utN2vvHmHOB9KB97TWyxq0UYPRc4Pm4jcdI77is1rrHQ3qFAOPmm9vaFLH8/4vZjlPnR3Ac2Z7V3lhb8jz3P1fAmDJbOwNlpRkJwczDmBu8AO4NsTt1sTWlOFK7LgbuSsxmQRnZdAcGjzfaXVS54cK/7v4WgAyKMfhcHW6PV+x7THmQTOOHNGh+nfdNZ6dJNOfUpx1oZeXNiengrffXseMGfO47LLPAbjqqrkdm18V/IKvwvUy0B942+z91KOUSsIYWksGXvSHcS1wgrmf28K5hRhxFCcqpRrmVWirzhdNynSmTljSB6PXYevXWLgsFkWuOQRVsWlHV5vVJsVF1aToSuyZGR1uw5WUTCpVVFY6/GiZ97j3Goks61I73+PqN9aIeNKb8qB8nsi9Ro/LnZXdofp9+sSxkyQicbFtcejNc5188gf89ref8/TTKxod//rrHcExSPA55NO7GOJ0HpCP2btRSi0DcoGzgee11p/718x6hpn7TS3Y5sRwDLEBg7yskwtUAv2UUjEASqlYIBOoMM83xeMf3i3WqvXGcAmPyMpsdi4HYyhuw7w1XWpTexTvLSIaJ6pXx3srruRU4qmlvLDCj5Z5z19+/18AalM7Lr71RESQRyx9KaOiouuFy7V1O0VEE5HW8R7wbU+dB0DuknX+MstvrF9vRGkZRAFTWY8d4zv++uvQ9LjtCfjsAqS1no4x/LYOY9hOAWOBLcDVWuub/WphYzyTGq3NqnuOJ3WgTmKTvS/XaIZS6jql1DKl1LL8EI33B0aPy43C2qd3s3OeHldybXAjTDTFk3vKlpHe8UZSDaeIqt1dn7pFa01UcR412Jj9rX/+NmrT+9KfUioquj5s0p7vV7ODFCIjOx6iNGmM8Xzp2BjcKPdNueOO79AaLmcZG3mcz3iNn3maXpTz5JPL2m9ACAgdCvmktZ6ptR4DxAH9gHit9SFa69f9al0Yo7V+WWs9Tms9Lj29EzfYAHNQXC3l9sRGi489PDfrfKqxYcsLrbxcG38wOs8ZI/q3U7J1LOmGcNXmdH2g3bo6F30oJ5d4jpvs9cqRNokYPID+lJCfX+2X9nxhIEXsIJnIyI5HkEs+fCRuFJZtwVmL1hpPPrmMg8jnFWazgEFcwiUMppBXmB1s03o0nYpVqLWu1lrnaK27yu2sae+oKZ7jJR2oU9pk78s1wpbkuhLKY1t2Yx5yUDI5JGALkpt1a7jzjF5K2vCOC5fNHGYs2tz1gXarqpxkUkouCTzwwIT2K3hB4sghZFHCjz/m+KU9b1k4fweDKGQzaaSldSCKiYktNpocWwoRO0JLuAAe43NqsXEJl/AuY3iAUziLdUxkO7W1oedM0hMItyC7G819s/klpZQNY52ZE9jmZZ0+GHnF9mitqwBMEd4LxJnnm+KJItpsziwcSXeWUhHX8loiu91GLglEFIRWjyuirAQAa6+O92Qj+xhOEU/d/7U/TPKJqioHWZSQH51GVJSvKfFaJnrYIBKpIX9r1/5Wleu3EoWLTaQTG+t9Hq6WKEjuR2J+aC14H0gh57KGz4eeRZ45dP4cR5NHLHczn5tvnhdkC3smvsYq3OblFqjHJs9fyWktnDsWiAF+0FrXellnSpMynakTdmit6eUuozqxNeGykkMCEfmh1eOKqigxXqR2fPGuxyMxlSo/WOQbleW19KeEISce6r9G+xu9T8e2rnUYiM/dAcAdL3qfg6s1qvsOoF/NvpBZ3FtSUsOD/TfgQnH2p39D6zvIyoqnikhe4iimsIH134THounuhq89LguGM0bTLRnINrfIDrTrLbOBAuBipdQ4z0GllB14xHz7QpM6rwO1wE0NFw0rpZKBe8y3Td33Pe/vNct56mQDvzfbC/v5vNoaJ70ppya55Z5LTEwEuSRg2R9aPS67R7hSOh7nz5ZhDBWm0vWLq2t37SUKF86+HR/qbIYpXPa8rh36tO824gtGjBze6baq+w0khWoqd4XG31tq8r84cfc8liSNxn7QQAD++McjAEi6+VqsaC5yiINGMPBpnEJrnd3aOaXUEOBfGENvXgfZVUpNA6aZbz2ubROUUjPN1wVa6zvM65cppa7FELDvlFLvYUSyOAvD7X028H4Tm7crpe40bVumlHofqMNYzNwPeFJr/WOTOj8opf4B3Ab8qpSajSHIF2GEs7q54WLmcKVqz35ScFGX0vJaoszMOHJIIJEaqKiAuLgutrBl7JWllKpoEiMiOtxGn8G9qMYWlB6X0wwmq7MGtFPSB7IMJ4/Yoq7tHUfv2U45UUT061wyTADnwMEAlK9YS1wnk2t2lqoqB0eyi0zKuKNkJEebx2+8cQw33jgGp9PNkmfu4qg9i9mzp5x+/eKDam9Pw289I631FuBcjDVQD/hQ9TDgCnPzCN6gBsfOb3Kdj4DjMBYcnwfcDDgwRObiliJwaK2fwRC3tcDlwHXAPuBKjyi2UOd2jOgY+8zyl5v1z9RaP+vD5wtZitfvAMDRyiJYpRTJBxtPmqEUPSO6qpRiS8eC63qIiLBSaU9gaErXR87QO4zhPOugbP812qcPbmXBnp/TpREddn27nE2kYY/u+ENEPUOMhdRFP63ufFudpKiohrNYhwMLX9C8N2mzWdg99gTGsYc/X/ufIFjYs/HrkJ4ZeulrwOsBb631g1pr1caW3UKdxVrrqVrrZK11tOmK/0+tdavxbrTWn2itj9Nax2utY7XW47XWb7Rj20yzXKxZ7zit9afefrZQ56s3jHjJ+y0JrZZxppsLZHO61lutLWKqSymxdr73VxmVQFxd1+fkspix/SLN4Se/YLNRldyLjNoCdu7smuDBZWW1DCWfzaQRFdXxNVwe9MCBuFDM/tvn1NQE11uvutrJWaxlAYOYfusxLZY57okbAbDO/ZytW0u60DohEHNRTg4M+QkhTKbViBox5ao2XLL7mEM2+0Jj3gEgtqaMElvnhasqOoEER9cLV2TOHgqJIbZ35yPDN8TRP5vBFLJjR+CDB9fVuUhL/CfZFLOJdL8IV3RSLDtJZjj5QQvF5cG9aRMHk8f/GInVqlosE3/koWwlhTNZx5Ah/26xjBAY/CpcSqk04BwgtHxahRaxmuuzkkdkt1omor8Rw9C1J3R6XHE1ZZRFdH5OoVDFEldbzs8/d90waG2tk70LVrKFVOLjO+c+3hQ9dChDKSAvL/Dzdq++upohFGJFG0OF9s679cfERLCG3oxiH9XVwe1xuT43wpR+xnD69m35ISk6JpK5DGMyW4lA1nN1JT79tSml7m+jnf4YsQoTgbs7aZfQBUQW5VFFJDEJrQ8VxvbPoBYrzh176Nyskp/QmoTaUsriWrfZWzYUWDiHCg4+4m1amer0O8XFxvDaQgYxLs4P80INsA4fRgqVLJ27losu6ryXX1tUVTk4BEPw9yZlo1TLvRJfsNksrKY3U9nA9pIKCKLDw7rnPySWJLaRyi23jG213PKkQ/h9yY8cya4utE7w9THpwXbOlwGPaK3/3jFzhK4ktqKIfFsCA9q46aSlx7KPeBJ3hUhCyYoK7K5a3L06H5y2MCKJVEdVlz4t15WWM4ASNpGG1erfkXr7aCOtyPczF/DhGYeQnGznhBP8E1KqKW635hD24cRCv5OP8EubsbERrKYPNty8cc9/eeR/v/dLuz7jdnM8W/mEEcyZM63N3+mqN27CdfYrXNkvtNY6dnd8Fa7jWznuBoqBDWaUdiHE0VpTs303u4mjLafstLRo9hFPfG5ozHG5cnKxAlEtRLP3laLIJKiCXnRdhHjHeiPgyib8H78y6pCDARhKPuef/z+AgPUktYZDyGUTaTz376l+afOgg5I59JLJ8M7bbPvkB4wlk12Pc+UvpFLFfIZwVVJUm2UnnXUIq+0DGLpnJcuX7+Pww2V6vyvwNa3Jgla2RVrrNSJa4UNJSS29Ka9P+94aaWnR5JKANUQWIeevNqJ5JRzU+cW7pfYkADK6ULj0xsAJFwMHoq1WhhL4bASffLKVQ9jHr/QhIaHtm7svnHfPmdRhrR+GDAaf3fkqAJXjj2GyF0GQP6sZyFHs5LhxrwbaNMGkXeFSSkUqpZYqpb5VSrU6KG+W+1YptaStckJoUFFRZ0Yob3uuyNPjiiwOjdQs+381Fu+mjuq8K/ld/zgTOJCTrCuwbDJCZ26m85mPmxEZiRo4kGENhGvFisAMYa36fiuDKGI1/l0oPHRkBoXpAxjFPq677iu/tu018+ezhVR6jRvWflngW4YQgZtJbA+wYYIHb3pcvwUOx4gw0aqPqta6DngcOAK41D/mCYGioqCMZKq97HHFE11eDI7guigD7PzJ6LEMObrzzgeDJhg3pgzKcbm6ZiFyxNrVbCWFSvzXS2nEiBGM5IBY7doVGNf4QzB64KsDsPKlfMBQRpPLK6/8isPR6tLMgJC7p5Tj2MZ8BuN0ereQe9rfr6AOK8cTWrnEujPeCNe5wDZvshprrediZAi+oLOGCYGleofh3t6ecNntNiL7m/NJ+4M/AR1RlI8LRcpQP8T5yzAcPHpT3mXZbB3LVrCSzs/PtcrYsQwnjxiMhJLnnPNxQILWTk0x8phFTDjS721HTjiCAZTQi/IuX8+18+OFJFHDfIbgdHr3MOOKimYJWRxP6KVk6a54I1xjgO98aHMhRhgnIYRx7DKCsea2I1wA/cYbGWF0CIR9iq0ootASD9bOL3glOpq66Dh6U86UKR92vr32KCtjkCufVfTl3XfPCMw1xozBimY0B9bdFRXV+P0yo6u2kRuZyutzr/J72zEnTALgKHZ2uXDF/PQ9APMZ7HUvfNq0IcxnMGPZCyUlAbRO8OCNcKUBvjxq7wf8GxJA8DuuvcaN7co/nthuWc8i5MotwV+rUrxhJ7lu/60oq01K6zLnjF2fGjfFuoNHc/HFAVpnNdZYczSmgXClpT1HRUWdXy9zaM029g88xK+OGR6iJ46nDitHsYuKiq4TrpUr95OwfDEbSGcfCZxzzkHtVwKyshKInnIyVjSu774LrJEC4J1wVQO+xNeJA/z/iCf4F7P3NPG8w9stGjnAGNqq2hbctVwul9srT0hfcKamd5lzxr8ufR6A7SlDAneRfv0gLY3xTYLXnHfe//x2ifLNOxlAMfkDD/Fbmw2JS09gFX1N4fKv4LbGZ59tZfzYmaSsW8Z8BpOb+zumTfNOuACKhx1GDTZWPvl++4WFTuONcO0GxrVb6gDjQJaRhzqWvDzcKGKzvZhvyTCjxwd5qNDhCIBwpfUio4uE6zi2sZF09jgDGINEKZg0iclN5lu++moHixb558HjnslGfIHaMeP90l5TlFIsIYsj2MWWdV0zr7plSwlj2UsCtcxnCAkJvoXjOv28g/mBAdi+XxAgC4WGeCNc32Hkx2pXvJRShwMTgfmdtEsIMM49ORSoWBJT27+JRsTGUEAMKi+4a7nqap1kUMF+nwYA2iaifyZ9ukK4XC4msZ0FDGLJkgA/AEyezECKGUBRo8PHHvseu3d3zsuwrKyWETkrKSeKqfdf2Km22uLqWbcQi4PnL3+RJ574OWDX8RARYal3rljAIGJifFvRM3p0OvMYwmHkQmFhIEwUGuCNcD0LaOADpdSI1goppYYDHwAu4Hn/mCcEiojC/RRHJWOxtB9jLirKai5CDq5XoWNfHnac7CXRb20mjRxEIjX0ig7s2nm9YgVJ1LCAQQG9DgDHGwFuTmjBPfuiiz5l+vS5HfI0/OabnQwc+AqnsIncYeOx2APk0g/YTjkJBxZOYyN33rmA2lonq1blBWzZwqpV+UxmK2vJ4Njz2x8+b0pCQhSlYyYabxZIryvQtCtcWuuNwEPAAGClUmqWUmq6UuoUc7tKKTULWAlkA3826wghTGJVMSX2ZK/KRkVZ2Uc8tvzgCpd7hzECvZsk/zU6wAh4lV5TGNAEjDXv/AcnFubi3aLWTjFyJLtJZBprm5368cccXn99DaWltT43e/LJH5BUtIchFNLr0rP8YWmrRKYls5hsTsO4ldx332LGjHmTxx5bGpDrvf7KSiaxne8YxLZtHctnljfgYKpUJMyXAadA41XIJ631Q8D/AQq4BHgF+MLc/m0eswD3aq0fCYypgj9Jri2mNMZ74colHkeQU5u4dxrCtcufwmWmvO+niwPnCKA1lv9+yDwGU9QVMfYtFuanH8VpbCSJllOc5OR470m5e3cZNTWGd98F/ApA0mWBXaqplGIuwxhDDv0o4emnVwCwfLn/H55eeGEVh7OHOOqYzxBmzz6zQ+04LREs0tmUfvylny0UmuJ1rEKt9aPAUOBhjDmsDeY2H6NHNlRr/ddAGCn4GbeblLoSymK9CzsUEWFhH/HElRcY0VWDRNVGI6SOX3tcpnANoJgXXljlv3YbsmgRUbu28765vHHWLP8EpW2LC2bfTyQuch5P5Oyzm3sx5uZW1r/WWrNsWcvzl1u3lpCV9TJ//OMiAC7iF/ZkHgzZ2QGxuyEfMBqAS1hJXZ0RQaO1pI6d4cYbv2k0vzVwYFKH2vnvfzczn8Ek7t4MeXl+tFBoiq9BdndqrR/QWp+ktR5pbidprR/UWndN6AGh0+St30UEbhZs9m5ep6LCQS4JROEK6gLLb15fRA02/vDwqf5rtE8f3BYrWZTwpz8t8l+7DXn8cWrik3mHMTz++HFceunBgblOA6InTYAjjyT6+acZPqS5F+ZDD/3I7Nkb2b+/kszMFxk/fhbvvLO+WblJk94F4F//WsGh7GUMOURf9duA2w+wjTQWM4DLWY4xzQ4WiyIvr5Lvv/fv0ozJbGU1vbnloY7/bd188xjmYzwkVHwapDiLPQT/JgQSwoLqrcYan+Io74YKDz449YALehBd4pPK8thDInkFflwmaLVSm96bLEoA/D9c+OWX8OmnLJt4ITVEcO21o/3bfmsoBffcA9u383D6Cl588WSGDEmqP71w4R4uuOATxo+fVd/7uvTSz8jPbzy02LBn9ge+p5IInJde1iUfAeDfHMlI9nMyRozK//xnI4MG/ZtJk97zSyir6moHUTiYxHbmMYSrr+742rR//vN4lpNJGVFseH52p20TWkeEqwdStdWYK7rq3lO8Kp+REcvR5kLlvF+DF4/tsNRadpPEddf5+eafNaBeuB566Ef/tbttG1x5JQwbxqfDzyU62ubz+qBOceaZcOaZRDx4P9ePKGX9+un8/e/HNiqye3fjpQDnnPMxr722mtraxr3xIeRzKSt5nfFE9Q5AZPsWWLHiMt5hDHtI5D6+xdPr8oSBKi/v3ENGTY2TmJinOY5txOBgLsNITOz472O1WnBhZREDiV/+Q0BiRAoGIlw9kL/fZkRRcGd4H9m7PM6I4vXE7Z8ExCZv6FWVzx5LMgcf7N8bp33YIA6KMm7gjz/uhzVDVVXw6qtw5JFQVwcffsiOfbX07h3rlxT3XqOUYUd2Npx6Kra/P8adV7cdtWPx4r1cffWXnHHGHAYNeqX++ON8Ri02HuFEoqN9zT/bMcaMyWDXvpt5Z8j5TGI7l7G80fnExGeYPPm9DrdfWFgNwBQ2UI2N/5U+RWxs5x4sfvnlCuYxhGHkU7NVZk8CRdf8BQohRV+MRai+CFdFXIrxYl+QhgqrqkisyGdHlP+jNaisLHo7i7Hgxo0Fp9ONzdaBZ7q1a2H6dFixApxOGD8e3nyT4oxs3n//C6ZO7XwOMZ9JTzfWFV13Hdx7L9x3H7XDhvPNRjf5xFKLDScWHFhxodCYwvpNw0Y0u0jir5Gn8shz5xMV1XW3jYyMWO5Y/xKuY9Yxc81nrKrsy2r61p9fsGAPeXmVREdHEB/vvejs2lXG9OlzATiNjXzHYKYkdH5h++jR6XyNESrqloPu5KmKtzsthkJzpMfVA+lDGUVE44rwfgHp0VMOpooIDknr2vxI9WwxFtPujgpAavSBA7G4nPXDhTNmzOtYO716QVQU3HEHzJsHP/3EF9sjSUl5FoDFi4O0nKBXL/joI1i5Eu67j8ghA+mtKpjMVs5gPRfwK5eznGtZynX81MK2lA8ZzXe9JnLNNV00R9cAi82K9cMPsCQmsDjiFaaxGs+wIUBGxgskJPyLZcv24XK5+etff2LPnrajofzud1/z7be7yKaI4eRjO91/np5XPvEbtpPM2axl9uxNfmtXOID0uHogfSgjlwRcLu/H4KeePpgdlgQGRbe8LijgbN4MwLIy/8UprGe4Ean94kM1j/0CH3+8lWefPcn3dtLTYeHCRoemTv1v/eu+fbtgDVdbHHaYsQFjgQ8+2MiLL/7CvHnehRY9/iA/5EDrKJmZsHgxcdOmMeeXN1lJX/7LKFaSyV4SKSCW08a/xP8+O5+/3PMNX8xZx5TTh/DPp1aQl3cjAC+9tIrTTx9McnIUtRU12HBxBusAOPkfN/jP1H7xfMxIbmAJaVd+zBVXjPJb24KBCFcPpC9l5BLvdaI8D1WJaahghX0yhSsgKe+HGdEsDnLnAans2VPOwIEvs2rVFSQm+i+s0aefnuu3tvzBBRcM44ILhnHZZZ8za9a6RucWLbqYSZOM+aMffriEuXO3c9BB3nmhBozsbNTPP8Nrr5F8/z95OK8Fl/PTHzSS1PxsbHcDRM4A4HqAm41iDUdCC4gh7SDvI8G3R2JiFC8ykhl8zylsYvTomfz665V+a18Q4eqR9KGcRQwkxscQR46UXqRsD9LQx+bN7FfxXHfr0f5vOy0NUlIYUJ0DGOE4d+woIynpGbS+o8PNNg2rNGhQUieMDBx//eukZsJ19NGZXHTRMD76aAsTJvRlwoS+rdTuYiIi4PrreXLNEN549kdGsY8MykmjCjsOrGhsuLDhxmoOJ/720hHMerv5GjUPv//jJMORxU8MH57CIgZSRDTTWMMVqwOT/qUnI8LV09CavqqcHJ3AjDN8C/hak5JO/63LcLu1V8F5/Ylev4GNOi0giQtRCoYNo98O/81B5eZWcPfdAVrQ7Gc8btuZmXHcfPMYfv21AKUU773XsdBHXUFGRgzl2PmR7HbLPvo2QMvCe8klI+j72Ol+tS07O5GH/nIcc+4dxYX8yu/ompxi3jJr1joOPjiVsWMzgm1KhxHnjJ5GURGR2knS8GyfvcOcqb1IoZry/M6lxvAZtxvnipWsoq9PnmM+MXw4A2uaC9eYMW+ycqXvw6NjxrzJG28cCHKblBS4SOqdJT09BoCHHz6aP/7xSN5+27838kBw662HN/tbWLDgIjIzffMMfO659jOAd4R77jmKi7/4K/HUci6rcTrd3HDD16xZkx+Q67XFV1/toFev59i1q4zExH9x2WWfc/jhb7F/fyXTpn3Egw8u7nKbOosIVw/DscsIlZMy0nfXbEum8dRasGabX21ql61biaitZhV9A7eG6LDDiCwu4ORRjfMwrVqVx9ixb7Fhg285lvbvP+DEUlMzg337fucXMwOB3W5D6zu46qrwGdKKjY2ktPRmyspuYeHCi3n99dM49tj+7NnjvZNFr14xJCXZA2fjqSeQE53BlSzjb39byksv/cI113R9KKg//nEh+fnVDBjwMmVlB3p/Z5/9ER9/vIU//9mPi+67CBGuHkblFsODzNLPi8zHTUgcbqQAyft1u19tapdVq4wdff05FdGYcUae1IMrW/5sxxzzHtu2lbBihdH70lpTUtI89NTSpbksWdK45xYVZevStU89BaUU8fGRTJrUjyuvPOC5N3HigWHBiy8eXv/6rrvGM2vWVJYvv4zHHz+O7767KNAGsufEczmRLbz5f4Z3qd1uDew1W2DVqpYD/v7004E1mZ7F2OGC/Df1MGq3G3EKPb0nX4ge2A8A5+69frWpXVauxIGFdQRwTP6ww8BiYXjlDqAf6enR5Ocf+GcuLKzmkENmUlXlZM6cs1mxYj8PP7yEnJwb6NMnDofDxciRM9m8ubhRs8OGpQTOZqFFFi++pNH7qCgrb7yxlgcemFif2bir5neWHX42oz99hdtZyPWcj8WieOSRHxk7NoOpU7sgqSigcHMc25jEduKpZQupfMxI9pNQX2bv3gpSU6O7xB5/ID2uHoZzpyFcEQP6+VzXPtBcx7O/5RQYAeOHH1hNH2rxLZ26T8TEwMEHc352KYMGJfLDD5eQkRHTqEhVlRG/75xzPubhh5cA8OKLvwBQWFjTTLQAvv02sHmrhPaZOXMKWt9RL1pdSURmH2YyjstZTiYl7NlTwX33Leb00/9LWVktDkfgFvS/8sqv7J39FSt5ivm8xIN8zU0s5iX+y27LY/yNT4nCiPu4b19lO62FFiJcPY1du8knlrheiT5XjRnQFxcKW14XCldNDXrJEuYzGGg89ON3JkwgbcMKtm68iiFDkvnss/bXXT300I84HC7efbe5u/VLL51MZmYAFkwLYcP06aNIfPR+AP7C3EYPN4mJz3DRRZ8G5Lp1dS5+vu4hMi6YQjLVXM5FxPMIMTzKwdxBxBW/5S4WsIjn6UMpxcU1PPfcSnbt6mLHqw4iwtXDsObsZhdJHVpYG5cUTS4JbPzmlwBY1gpLlqBqa/mOwbz44skBnUzn5JOhrAyWGunhvXUEiYz8J7fd9l2z4+E09CIEBqvVwm/unsJTTOIKljOJxo5Nc+Zs5oMPNjJnzma/Xtf1+JO8zId8xVBGcTtvMY6/PzsFUKwnA157DT1nDodEFDKfl5hx8TvcdNO3nH32R361I1CIcPUwIvft7bBwWa0WdpLEAEq6LmXDV1/hVhYWMZCjjuoT2GudeKKxpusrw/OrswuGA26vEDY8wolsIo23eYdUGg/LXXjhJ5x77sfMn+9d6K12ee01ov/vj/yH0ZzFlZRjp2/fOG688TC++eYCFi/+DQBq2jTU3C/IpJRveJkEqpvlYwtVRLh6GLacPewiuX7tjq9UpvYhi+L6nEgBRWuYPZvN/Q7DERPPyJEBzgOVkgITJsCHH4LW2O021q27iscfP47nnz+J/PwbeeONKY281lrioYeOJjf3dzJMKNRTSRSXcAlpVDKXf5NIcy++1asLUOqJZlFMvOV//9vCh3e9hfPa69k7aiK/5Te4MLwYU1PtKKU48cQBTJx4wKM48vhjefaUPzOUfN7hHXL3lvH666tD3stQvAp7EqWlJFDLLpI6vJC317gR9P/yZwpKaoiLC3C6hqVLYfNm5h32OwYmJHYs1YivXH453HADLF8O48YxYkQqI0akNjg9klNOyaZPnxcaVfvii/NYuTKPtWsLmD59FL17BzmgrhBSfP75uezfX8X5V1UwhzdYxtNczsWNIn/84Q9GVoLLLvucykoHV1wxErvdu1v0p59u5aqz32YlT7GbeMauORUHNkaMSGH9+qJWHzSVUvzpyz+x6TY3p//zXv7K50yfbuGss7bw8cfndPpzBwrpcfUgHFuNNUoTLjii42307U8kLj55uQvCGT39NMTH89+ow7tOCC66COx2eOaZVov07h3Lq6+eSk6Osdj1kEPSOO20gdx995HMmnW69LSEZkyZMogrrhjJ54zgzSv/hR0nP/Acn/NvfstysimiYaqWG274mvj4f/Hdd94NH5515ofM4l0yKOd8LqcEY0TlxhsPw2az8NBDbcf4HPqPe3gt8mjuYgHn80ujBfShiPS4ehBla7aQCliyB3S4jcEnHgqvwxsPf8Gld50cuF7XsmXw3ntw112se8vBySd3kRgkJcFNN8GTT8Ltt8PolvNPTZ9uRJkoLPy910/FQs9GKVUftHnK1lgOX/Qu17OEtzCi8FcSQS4J5BNLDRHUOq2UHP8ae4/qR2a/eFAKl1YU/2Y6P9gGc9ZZQ7juuq/46adc7uVbprCR6zmXFRxY6mKsMbzNK/tmT/o9I77dzev8h5v0GP9/AX5E/uN6ENUbDY+mqIM6nok35TAjBUgWJRQXB2i4sKAALr4Y+vSh9MbbyPnbmwwd2oUpNe6+G954w+h9/fADJLd+7ZQU8RwUfCexbwp/4SQe5QQOJZcj2cVQ8ulDOalUYsdJMnXYKce9tozyPZEUFVUzvupa8mevBlbXt3USm/gzX/MWY3mZo+qPP/nkZM491/t0LW+9P43RaetYwVM8uuGfUHYNJCS0XzEIiHD1IBybt1GHlfghnUgImJUFwACKKS2tpb8/cws6nTB3LvzhD7B3L8ybx6PPGeuj2nOI8CspKfCf/8App8BRR8FTTxmvrV0frkfonrz44slMntwfm83C1KmDyMx8sfXC5ebWApmU8A7vsI5e3MB5gOKuu8bz0ku/cs01h6B8iJGWmhrNI69dyEXTC/im7GX2Tb2QxTOe5pxzh3Z5Noj2UF3m1tyDGTdunF62bFmwzWDb2JOoXbma2J2bycrq+JNUTWwir1cdzKiFHzBpku8ROJgzB+bPh+pqqKkx9vv2wdq1UFICgwbBrFkwYQLHH/8+3323G7f7dp/+Cf3CwoWGs8bOnUava9Qo6N0b4uMhOhosFsN9vul2//2Q6PsCb6HnUllZxymnzGbs2AyefXalV3UicbKAFxjJfsZzCxvpBUBFxS1ER0d0WGyOOeZdjlz8Dk/yKXdyOsNeeZRrrml5yDyQKKWWa63HtXROelw9iOg921mv0pjSr3PzRe6Bgxm8tpDc3IqONfDjj/DWW8bN37OlpsKFF8Kpp8IZZ0CkMQRZUFDNtGlDul60AI49FjZsgE8+ga+/Nl6vXg0VFVBVZbjru93GvuF2++0iXIJPxMZG1sdYHDEihd///ttmZaZOHYjVauGTT7YC8A/+x1Hs4gJ+y6EXTmLjfzYCEBMT0an/lz59YvkHx3Iku3iMz3nju5MhCMLVFiJcPQW3m5SiveTGHdPpbn/EqOEMXTuXd7eWdqyBv//d2LwgP78quNl37Xa44AJjE4Qu4MYbx3D00ZlYLIrRo98AYMiQJD777Dy01nz77S4sb7/FCTN/5HGOYzaHot8/k4qKOs49d2inH/KefvoEli3bx9U7LmAU+zjr3f/D8dDZWLMHhMyQobjD9xRyc4ly1VKcltXppiIOHk4WJaxeuhO3O3BDzVprCgqqSUsTBwihZ3Hoob045JB0NmyYDsDNN48FDM/Ek6J2ccK7j6CPP567mVJf57PPzuPqqzufU61v3zjmzJlGBXbO5Qoi3Q5WDD6WaVPe73Tb/kKEq6ew2YiFVtW3467w9QwdigXNrx8tZerUDzvfXivMm7cLl0t3TZQOQQhBhg1LISfnBm6+2XRPX7sWzjoLBg5EffABLqxcc43/E4AedpgxX7aRXlzJRRzJbs766p846px+v1ZHEOHqIeQsMCZ8K/pkd76xoUONHfl8+eWOzrfXArW1Tk466QOgiz0KBSHE6NMnzhj+q6qCKVOMOeG5cyE1lbq6W3nppVMCct2lSy/l7ruPZA6H8Agncg1LeSFqMnt2N48gX1ZW26WR5UW4egg/vrWQWqyc8bvjOt9YA+EaNCgwTggNPasuvHBYQK4hCGFFTAw89hh88QUMMEZOIiKsAZt3Gj++T30aofs4lSc4lltYzN6pF0NdHVde+QXPPWf8n44bN4sBA14OiB0tIcLVQ8iq2M0m0jnuhOzONxYXB337cvoQB7W1gUmEt2LFgXTjQfEoFIRQ5JJL4NBDu+xyo0enM3PmaQwcmMSdnMFDnMSRa75gUdRwvn7jB266yfB+bCmJaiAR4eohDCzbyfaEbP81OGIEA6v2kptbSWlprf/aNVm9Oh+A44/35wpnQRB85YorRrFt27U8/fSJPMCpXMSlHM4e1vIE17IEm/pbfVmlnuDpp5cH3CYRrp5AcTFp1QVsj+u8R2E9hx1Gn4JtKLeTk0/+wH/tAhs2FLJ6dQEAc+ZM82vbgiB0jFtuGYvWd/AfDmM0t7GSTF7mQ9aZAuZJ1TJjxvwWM4L7ExGunsBqI67ZzsRs/7V52GFY62oZTj4//7yPBQt2U1PjH4+j//73QDbYjiS8FAQhcDz66CS2ksaJXMc0rqCSSF7mQ/bzZz7n39zNt/z7kqfY9cNaY0F+ABDh6gn88gsAe1M6Hly3GWMM99zD2AvA5Mnvc/31X3W62aKiau699/tOtyMIQmC4++4jyc+/kcFDUviYUYxlBuO5heeZSBYlPMpcvuVlso4eRWViOtrt9rsNEjmjB6CXLGGfJZGqpHT/NTpsGNjtHF6zl7c5HIA331xHv37x/OUvkzrc7I8/5vjLQkEQAkRaWgybN1/DjTd+zQsv/MIy+rMMYz46hUrGsJeDKCC2vI5HndoTwc1vSI+rB+BcsIjv3QM49jg/OjrYbDB+PKfE7m10+NFHf+Knn3Lr3zscrnaT4b3zznoefvhHtNZ89dVOwPBmuvvuI/1nryAIfufxx4/jpZdOxu2+nVdfPZVbbz2cImL5lqG8yEQuWPIKkZH+z6og0eG7gKBGh9+7F/r1YwZncenSFxk/vo//2r7vPvRf/0qC60EqsDc6dfXVh/Dyy6eQnv4cRUU19O4dy5o1V1JV5aCkpJYNG4pISbFz4okDUOqJZk17Eu4JghB+LFy4m8LCGs45x/t8YE2R6PA9mUWLAFhMNn8dlebfto89FvXII5TPHU/29TvZufPAyvlXX13Nq68eSHa3b18lDzywmNmzNzVKCz537nnNmj30UD8OaQqC0OUce2xgl7HIUGE3R3/6KfnEsoJMoqMj/Nv4hAlG+pEvv2TVqsvbLf7cc6saiRbAaac1j3X4yiuBCWEjCEL3QISrO+N0oj/7nM8Yzu13BmC+KC4OTj4Z5swhKTEKre+gtPTmDjd3553j+emnS/07nCkIQrdDhKs78+23WEqK+R8jOeQQPw8Tejj3XNixA1asACAhIYrS0pvZs+d6+vc/kLBy+vRRjarZ7VZ++unS+vcLF17M3/9+HEccIaIlCELbyBxXd+bFF6mISeKzqhFc3ysmMNeYNg1uuglefBFeeQUwxCshIYpdu64nJ6eCHTtKmTgxk9tuG8dnn21j9+5yrrxyJIcf3psff7wEt1szcWJmYOwTBKHbIcLVXdm4Ef3JJzzrOpY6bGRkxAbmOikpcNll8MYb8PDD0Lt3o9N9+8bRt28cACNHpjFyZOOe31FHScoSQRB8Q4YKuyNaw623QkwM/+BYhgxJCqyn3h13GNe8/fbAXUMQBMFEhKu7oTU88AB88QVbr7iNfOJ4+OFjApsa5KCD4E9/gnfegaefDtx1BEEQEOHqPrhcsHgxnHOOMWR35ZUc9KwxPDhmTK/AX/++++Dss2HGDCNn0LJlhk2CIAh+Rua4QpmLLjLSdbtcxuZ2H3jd8H1pKezcCTU1hov6Y49RfsMfYOazAAwdmhx4W202+OADeOghePJJePddiI2FrCxITTXWe0VEGJvFAmPHGj1DQRAEHxHhCmV27oS6OrBajc1iOfA6KurAsawsOPNMGD8eTj0VEhN5/V+Ge/qJJ2Z1XQbhiAijtzdjBsydC0uXGiGnioqMz1FRAQ6HMZzZr1/X2CQIQrdDYhW2g1KqH/AQcBqQCuQCHwF/1lp7la+6q2MVut0aq/VJ8/XtXSdcgiAIfkJiFXYQpdRg4AegF/AxsAE4AvgDcJpS6mitdWEQTWyG0+lm4sR3AGOIUERLEITuhghX2zyPIVq3aK2f8RxUSv0DuBX4C3BDkGxrkYiIf9S//vnn3wbREkEQhMAgXoWtYPa2TgF2AM81Of0AUAlcppQK0Mpe33niiZ/rXz/22CQSEiTtvSAI3Q/pcbXO8eb+K611o9zTWutypdRiDGE7Cvi2q40DWL++kB07Srnkks8oKamtP37mmYO5664jgmGSIAhCwBHhap1h5n5TK+c3YwjXUAIkXIcf/hYVFXU4HG6cTjcOh5vKSgcWi6K0tLbFOp99di5Tpw4KhDmCIAghgQhX6ySa+9JWznuOJ7V0Uil1HXAdQFZWVocMGDEiBZdLY7NZiIiwYLNZsNkUubmV5OVVsXVrCVarhYsvHsZVV41i1ChJwCgIQvdHhCtAaK1fBl4Gwx2+I23MmnW6X20SBEHoDohzRut4elSJrZz3HC8JvCmCIAiCBxGu1tlo7oe2cv4gc9/aHJggCIIQAES4Wme+uT9FKdXoe1JKxQNHA1XAkq42TBAEoScjwtUKWuutwFdANvD7Jqf/DMQCb2mtK7vYNEEQhB6NOGe0zY0YIZ/+pZQ6EVgPHImxxmsTcG8QbRMEQeiRSI+rDcxe1zhgJoZg3Q4MBp4Gjgq1OIWCIAg9AelxtYPWejdwVbDtEARBEAykxyUIgiCEFZKPqwtQSuUDO4NthxekAQXBNkLwCfnNwgv5vbxngNa6xXBAIlxCPUqpZa0lbhNCE/nNwgv5vfyDDBUKgiAIYYUIlyAIghBWiHAJDXk52AYIPiO/WXghv5cfkDkuQRAEIayQHpcgCIIQVohwCYIgCGGFCJcgCIIQVohw9XCUUv2UUq8ppXKUUrVKqR1KqaeUUsnBtq07oJQ6Xyn1jFJqkVKqTCmllVKz2qkzUSn1uVKqSClVrZT6VSk1QyllbaPOGUqp75RSpUqpCqXUT0qpK9q5zhVKqaVm+VKz/hkd/azdAaVUqlLqGqXUHKXUFvP7L1VKfa+UurppiqMG9eQ360q01rL10A0jYPB+QAMfAY8B88z3G4DUYNsY7huwyvw+yzGyC2hgVhvlzwacQAXwKvC4+Vto4INW6txkni8AngP+Cew2jz3RSp0nzPO7zfLPAYXmsZuC/b0F8fe6wfwOcoC3gb8Cr2FkOtfAbEynNvnNgvg7BdsA2YL448OX5h/9zU2O/8M8/mKwbQz3DSMFzkGAAia3JVxAApAH1ALjGhy3Y6TX0cDFTepkAzXmDSy7wfFkYItZZ0KTOhPN41uA5CZtFZrtZXfmc4frBpwAnAlYmhzvDewyv7fz5DcL7iZDhT0UpdRg4BRgB8aTW0MeACqBy5RSsV1sWrdCaz1fa71Zm3eZdjgfSAfe01ova9BGDfB/5tvfNakzHYgCntVa72hQpxh41Hx7Q5M6nvd/Mct56uzA+FuIoodmRNBaz9Naf6K1djc5vg940Xw7ucEp+c2CgAhXz+V4c/9VC/+k5cBiIAY4qqsN68GcYO7ntnBuIVAFTFRKRXlZ54smZTpTRwCHuXc2OCa/WRAQ4eq5DDP3m1o5v9ncD+0CWwSDVn8TrbUT2I6RQ2+Ql3VyMXrO/ZRSMQBmDzoTqDDPN0V+9xZQStmAy823DcVDfrMgIMLVc0k096WtnPccTwq8KYJJR34Tb+skNtnL7+4bjwGjgM+11l82OC6/WRAQ4RIEQWgDpdQtwO0YnoKXBdkcARGunkzTp7qmeI6XBN4UwaQjv4m3dUqb7OV39wKl1E3A08A64HitdVGTIvKbBQERrp7LRnPf2rj4Qea+tTkwwf+0+puYcywDMRwDtnlZpw8QC+zRWlcBaK0rgb1AnHm+KfK7myilZgDPAGswRGtfC8XkNwsCIlw9l/nm/pSm0QCUUvHA0RgeUUu62rAezDxzf1oL547F8PL8QWtd62WdKU3KdKZOj0Ip9UeMRb6rMEQrr5Wi8psFg2AvJJMteBuyALmrv+/JtL8AOR/fFrMORBaz+vt3us/8fpYBKe2Uld8sCJvk4+rBmIuQfwB6AR9jhCQ6EmON1yZgota6MHgWhj9KqWnANPNtb+BUjGGjReaxAq31HU3Kz8a4Eb0HFAFnYbhQzwYu1E3+aZVSNwP/wriBvQ/UYSyM7Qc82bD9BnWeBG4D9pjtRgIXAakYDzLPduqDhylmrMCZgAtjmLAlT74dWuuZDepMQ36zriXYyilbcDegP/A6kIvxz7MTeIoGT3Wyder7fRDjSbm1bUcLdY4GPgeKgWpgNXArYG3jOmcCCzBiIlYCPwNXtGPblWa5SrPeAuCMYH9nIf57aeA7+c2Cu0mPSxAEQQgrxDlDEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuAShB6OUmqmU0kqp7GDbIgjeIsIlCB1EKTVcKfWMUmqNUqpUKVWnlMpRSn2mlLq6Sbr2bodSqr9SymUK36PtlNVNtlqlVL5SaoVS6t9KqSlKKWs7bcQqpWYopeYppfLM77tEKbVUKfUXpdSgtuoL3QeJnCEIHUApdT/wAMbD348YAVkrgAyMYLqDgOVa63HBstEblFIzgSuAgVrrHT7W/TNwP0YYpP1Af22kq2+prOdG82dzb8XI2DsSI1xSJMZ3eKnWull6DqXUURjx+TIxYvV9C+RgpAAZgxGE1g0cpbVe4cvnEMIPW7ANEIRwQyl1D8YNeDdwgdb6pxbKnIGRNbdbYvaOpgNlwCzgRozAsv9tq57W+sEW2srACGh7AfCNUmqcbpBGRCk1HCOTQRzwJ4wgtM4mbQwE/oYRrV3o5shQoSD4gDkX9CDgAKa2JFoAWutPaZI7SSl1pVLqQ6XUNqVUtVKqTCm1WCn121au9Z05rGZTSt2jlNpsDrHtVkr9TSkV2UKdaUqpWUqpTUqpSnNbrpS6pWnetU4yBSOS+fvAC+axazvSkNZ6P3Ax8B1G0Od7mhR5BkOQ/qa1/ltLvTqt9Xat9YUYvV+hmyPCJQi+cRUQAXyotV7TVkHdOHkgGDf4AcBCjAj875nv31JKPdxGU+8AN2OkQnkBI/r4XcBLLZR9DBgL/IRxw38To6fyNPBGW/b6yHXmfqb5PSzHSEo6oCONaa3dwCPm298opRTU96ROwkgZ8ncv2mn6nQvdEBkqFATfOMbcf9uBuqO01lsbHjB7TV8Af1JKvai13ttCvcHASK11kVnnXuAX4HKl1N26cUr501u4hgUjdc3lSqlnW+sleotSKhOYCmzSWv9gHp4JHA5cg5GIsSN8j5HmvhdGgsTtHPi+l2utSzrYrtDNkB6XIPhGH3O/x9eKTQXFPFYHPIfxEHliK1X/6BEts04l8DbG/28j549WruHG6HGBkciys0zHcK6Y2eDYOxj53Ka35x3YGmZvyZO4NN3cd/j7FrovIlyC0EUopbKUUs8ppTYopao8ruHAh2aRzFaqLmvh2G5zn9zkGqlKqceUUr8qpSoaXGN5O9fw9jNYgKsxPPje9Bw3hfUToC9wemcu4WmyE20I3RwZKhQE38gFRuCjAJhrjJZiCM0i4CuMtPAujGGxK4AW1321MkTmcVCo790opZIwsuMONK/1JkYaeSeG6/kfWruGD5yKMS/3ZQvDmjOB8zDmv/7na8NKKTuQYr7NN/e55r5Tgit0L0S4BME3vgdOwBjWe9WHercBqcBVWuuZDU8opX6DIVyd5RoM0fpzU7dzpdQEDOHqLB6njFMbrM1qymlKqf5a692tnG+NYzDuSfsbrCn73tyPU0olaq1LfWxT6IbIUKEg+MbrGK7w5ymlDm6rYJPIGUPM/YctFD3OT7YF9BpKqd7AGRhrt15tZVuM0Quc7mPbFuBe8+07nuNa6+3AN4AduNOLdrp1tBLBQIRLEHzA7Ak8iBHp4TOlVIuRMZRSp2F4C3rYYe4nNyl3KkZPyR+0do0xwN1+aH86Ro/oba31NS1twJUY81NXe7tuTCnVC2NpwGRgF9A0fNQtGGJ5t1LqdqVUs5Eic/7wPWBCBz+bEEbIUKEg+IjW+lHz5vkA8LNS6gcah3w6FjiIxk4Vz2OsAftAKTUbI1zRKIxFyv8BLvKDaW9i9EqeUkodD2w27TgDI6JFh69hrqvyCOy/Wyuntd6ilFqAIUJTgM+atPOg+dLCgZBPx2A8CCzFCPlU0KTN9abAzwaeAP6glGoY8ulQjLBRGiN6htDNEeEShA6gtX5IKfUBRqij4zFEyY7hzr0K4wY6q0H5X00xeQTD686GsRbrXKAEPwiX1jpHKTUJYxHyMRiOFBtMG7/p5DVOwpg/W+lFLMBXMITrOpoIF4bYg+E6Xw7sxBDcD4GvTNf9Zmitl5ihn67FCC11OoajSxWwBXgSeNkcWhS6ORJkVxAEQQgrZI5LEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCv+HxQ2gk/QavdsAAAAAElFTkSuQmCC\n", + "text/plain": [ + "<Figure size 432x288 with 1 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fit(0,250)\n", + "fit(250,500)\n", + "fit(500,750)\n", + "fit(750,1000)\n", + "fit(1000,1200)\n", + "fit(1200, 1400)" + ] + }, + { + "cell_type": "markdown", + "id": "12ee35ee", + "metadata": {}, + "source": [ + "### Cálculo de la distancia entre picos\n", + "\n", + "- Primero se buscan los máximos locales y después se calcula la distancia entre ellos." + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "id": "8f6d8519", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAApkAAAHyCAYAAACgdNVWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB1NklEQVR4nO3dd5wcdf3H8dfncqkXUmihk9CRXkVqABWkl9AVBAVRQEpEARXBHypKE0FEEIlIlx6qaDiqKXRRhAQSekIgJCHJpd19f398Z7J7e7t7W2Zmd27fz8fjHt/d2ZnZ781tcu/7fuf7/ZpzDhERERGRKDXVugIiIiIi0vMoZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiESuudYVaAQrrriiGz58eKzvMX/+fFpaWmJ9j0ak6xofXdt46LrGR9c2Hrqu8ajmur7wwgufOOdWqrYOCpkJGD58OM8//3ys79Ha2srIkSNjfY9GpOsaH13beOi6xkfXNh66rvGo5rqa2TtR1EHd5SIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5ApIiIiIpFTyBQRERGRyClkioiIiEjkFDJFREREJHIKmSIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5Ap6eIcTJgAhx0GLS3Q1OTLww+HiRP96yIiIlJzCpmSHkuWwNFHwx57wD33wIIFPlQuWAB33+23H320309ERERqSiFT0sE5OPZYeOABHyo7Ojq/3tEB8+fD/ff7/WrZoqnWVhEREYVMSYmJE2HsWB8wi2lr8/tNmpRMvXKptVVERARQyJS0uOwyHyBL0dbm909amlpbRUREYqaQKenw0ENdQts1fBfDMYdBnfft6PD7Jy0tra0iIiIJUMiUdMjTinkFZwIwg2El7R+7Aq2t77N6131r1doqIiKSEIVMSYf+/btsaqcXAM0sLWn/2OVpbb2dI1iT92llt8771qq1VUREJCEKmZIO++7rR2lnWUozAL1o77xvU5PfP2l5WjGfZScAXmXzkvYXERHpKRQyJR1Gj+7SOhm2ZDqs8779+vn9k5an9XQeAwEw8gzyqUVrq4iISEIUMiUdtt8e9t+/UzALWzLDsAn41w84ALbbLuka5m1tHcPxQJ6QWavWVhERkYQoZEo6mMFNN8GBBy6b4LxTyGxqggED/Os33eT3T1qe1tZQEznTGdWqtVVERCQhCpmSHr17w623wrhxcOihzGIFANr7LwejRkFrK9x2m9+vFvK0toY6hcxatraKiIgkRCFT0sXMh7k771y26Y8nvQB33FH70JbT2vqZLZ95CVcfra0iIiIJUciUVFq8OPP4yitrV48uslpb9x46ftlm69OnPlpbRUREEtJc6wqIVGLOnFrXoIigtXXirKxNV18FJ9auSiIiIklTS6akUu4Uk//5T23qUSr1jIuISKNRyJRUcjkzAm26aW3qUaom/UsTEZEGo199kko5qzfWPYVMERFpNPrVJ6mU25IpIiIi9UUhU1IpbS2Z7e3d7yMiItKTKGRKKqklU0REpL4pZEoqpa0lc8iQWtdAREQkWQqZkkppC5lqeRURkUajkCmpVO+h7d13Oz9PWygWERGplkKmpFK9h7aTT+78vN7rKyIiEjWFTEmlem/JfOSRzs8VMkVEpNEoZEoqpS20ffObta6BiIhIshQyJZXSEjLDUeVLltS0GiIiIolTyJRUqvfu8lCfPrWugYiISG0oZEoqpaUlUyFTREQalUKmpFJaWjJ79651DURERGpDIVNSKS0tmf3717oGIiIitaGQKak0YUKta1CagQNrXQMREZHaUMiUVDr9dF8+9lht69EdhUwREWlUCpmSamaZx/XShZ5dj5/+tHb1EBERqSWFTEm1pqxPcL20ara3Zx5vvnnt6iEiIlJLCpmSatktmfUy4fnSpZnHZjB8eM2qIiIiUjMKmZJq9RIss2WHzKYmOOII6Nu3dvURERGpBYVMSbWFCzOPm+rk0zxnTuZxU5OfK3Px4vTM7SkiIhKFOvm1LFKZ7JCZ3XVeS9/6VuaxmV/1x7nO92qKiIj0dAqZkmrZIbNevPFG5nFTU2ZpyXrs2hcREYmLQqakWj22ZGZ3i5tllpZcvLg29REREakFhUxJtVVXzTyux5CZ3ZKpkCkiIo1EIVNSaYcdfHnAAbWtRz6FQqa6y0VEpJEoZEoqjR/fdVu9jC5Xd7mIiIhCpvQg6i4XERGpHwqZkjqF1iivl7XLc1sy1V0uIiKNSCFTUid3vsmzz86/vVbUXS4iIqKQKSkULtu4006+POIIX9ZjyAR1l4uISGNSyJTUCUPmwQf7slcvX9ZLyNx2287P1V0uIiKNSCFTUicMmc3NvgxDZr3ck7nLLp2fq7tcREQakUKmpE5uyAynLqqXlszcsKvuchERaUQKmZI6hVoy6y1k7rGHL9VdLiIijUghU1Kn3kNmOPDnscd8qe5yERFpRAqZkjr1HjLDlsxwcnh1l4uISCNSyJTUqfeQGYbJ8F5RdZeLiEgjUsiU1Kn30eW/+IUvw5ZMdZeLiEgjUsiU1Kn30eW51F0uIiKNSCFTUqdQS+Y119SmPt0ZPNi3as6cWeuaiIiIJEchU1KnUMh89dXa1Kc7/frBsGHw/vu1romIiEhyFDIldQqFzHrWvz8sWlTrWoiIiCSn5iHTzFYws2+b2b1mNsXM2sxsjpk9Y2bfMrO8dTSzHc3sYTObFRzzqpmdYWYFI4eZ7WdmrcH555nZBDM7rpv6HWdmE4P95wTH71ft9y2VS2PI7NsXFi6sdS1ERESSU/OQCRwGXA98EZgA/Ba4G9gU+BNwp1k4TtczswOBp4BdgXuBq4E+wBXA7fnexMxOBcYG5705eM/VgDFmdmmBYy4FxgCrBvvfDGwGjA3OJzVQaOBPPevXTy2ZIiLSWJprXQHgTeAA4CHn3LJJaMzsPGAicChwCD54YmaD8IGvHRjpnHs+2P5TYBwwysyOdM7dnnWu4cClwCxgW+fctGD7z4FJwGgzu9s596+sY3YERgNvAds55z4Ltl8CvABcamYPhueS5KglU0REpP7VvA3IOTfOOTc2O2AG26cD1wZPR2a9NApYCbg9DJjB/guBnwRPv5vzNicAfYGrs0NhEBx/GTw9OeeY8PkvwoAZHDMN+H1wvuO7/w4lamkMmWrJFBGRRlPzkNmNcI2UpVnb9gjKR/Ps/xSwANjRzPqWeMwjOftUc4wkII0hs08feOaZzLrmIiIiPV3dhkwzawaODZ5mB70Ng/LN3GOcc0uBqfjbANYp8ZiPgPnAGmY2IHjvFmB1YF7weq7JQblBSd+MRCoMmWG4TMM9mY8/7ieLv+GGWtdEREQkGfVwT2YhF+MH6TzsnHssa/vgoJxT4Lhw+5Ayj2kJ9ltQ4Xt0YmYnAScBDBs2jNbW1kK7RmLevHmxv0e9eP75YcDGvPbaBGbPbgNglVW+SJ8+HbS2Tor0vSq7riMBco7z25544h3WW29q9RXrARrpM5skXdf46NrGQ9c1HvVwXesyZJrZ9/GDbv4HfKPG1amIc+464DqAbbfd1o0cOTLW92ttbSXu96gXr73my732+iIrreQfb7cdvPcekV+DSq7rkCHwjW/kr8tmm63NyJFrR1K3tGukz2ySdF3jo2sbD13XeNTDda27jsZgaqArgf8CuzvnZuXsErYiDia/cPvsCo6Zk1OW8x6SkNmzfTk466fT3JzpRq+ljg5fv0L3iTbX5Z91IiIi0aurkGlmZwBXAa/hA+b0PLu9EZRd7ocM7uMcgR8o9HaJx6yK7yp/3zm3AMA5Nx/4ABgYvJ5r/aDsco+nxG/2bBgwwA+mCTU3+3sea+388335+ef5X+8846uIiEjPVTch08x+hJ9M/WV8wPy4wK7jgnLvPK/tCgwAnnPOZU8YU+yYr+XsU80xkoDZs2Ho0M7bPv8cXn8dlizJe0hibrvNl3MK3c0rIiLSIOoiZAYTqV+Mn+R8T+fcJ0V2vwv4BDjSzLbNOkc/4KLg6R9yjrkRWAScGkzMHh4zFDgveHptzjHh8x8H+4XHDAdOCc53Y3ffm0Rv7lxYbrnO2x4N5h+4Pe96T8kr1F2ehpHwIiIiUaj5HWLB2uE/x6/g8zTwfevapzjNOTcGwDk318xOxIfNVjO7Hb+SzwH4qYruAu7IPtg5N9XMzgZ+BzxvZncAi/ETu68BXJa92k9wzHNmdjlwFvCqmd2FX7ryCGB54DSt9lMb7e2F722s9TyU4Ue3UJjs6Mi/XUREpKepecjE30MJ0As4o8A+T+LXEAfAOXefme0G/Bi/7GQ/YAo+EP7Oua5Rwzl3lZlNA36An3+zCT+46CfOub/ke1Pn3Ggz+ze+5fIkoAN4EbjEOfdgWd+lRKajo3CIq/XE7GGILFSPWnfni4iIJKXmIdM5dwFwQQXHPQvsU+YxY4GxZR4zhqyAK7WXxpD54x/DL35RH4OTREREkqA7xCR1ioXMWt/zGIbI3JB54YWdXxcREenpFDIldeo5ZOYueRkK66WQKSIijUIhU1Injd3lZr7OCpkiItIoFDIldeo5ZIYhMl/9evVSyBQRkcahkCmpU8/d5cVGlytkiohII1HIlNRJQ0umQqaIiDQ6hUxJnWIhs9a6C5majF1ERBpFnf6qFiksX8gcNMiX9dKS2XXRKg38ERGRxqKQKamTL2TeESwk2rt38vXJpoE/IiIinkKmpE6+kDlggC9rvXa5QqaIiIinkCmpky9kht3TtQ6Z4T2X+brLFTJFRKSRKGRK6hQLmfUysEYhU0REGp1CpqROvpAZPq91S2ZIIVNERBqdQqakTj13l4d0T6aIiDQ6hUxJnRdfhHnzOm+rt5CplkwREWl0CpmSKkuX+rK1tfP2ersnUy2ZIiLS6BQyJVXa2vJvr7d7MhUyRUSk0SlkSqoUCplp6C5vaqqfllYREZG4KWRKqoQhc6+9Om8PQ91f/pJsfXJ95Su+POywrq+pJVNERBqJQqakShgyv/nNztvDkHnXXYlWp4uNNoIhQ2DDDbu+ppApIiKNRCFTUiUMmf37d96e7x7IWsg3vVJIIVNERBpJnfxqFilNGDL79eu8Pd89kLXgXOG6KGSKiEgjUciUVCnUklkvA2qcU0umiIgIKGRKyhQKmfUS3jo61JIpIiIC0FzrCoiUI80tme3t8PHHydZHRESkVhQyJVXqPWRed13h1555Jrl6iIiI1Jq6yyVV6r27vBT1EohFRETipJApqVLvLZmlWLSo1jUQERGJn0KmpEo9h8xSW1MVMkVEpBEoZEqqtLX50dt9+3beXg8hc8mS4q///ve+rGnIdA4mTPDrXra0+FFKLS1w+OEwcWL9LP4uIiKpp5ApqbJoEfTp03WaoHq4J3Px4uKvh8G4ZiFzyRI4+mjYYw+45x5YsMCHygUL4O67/fajj+4+LYuIiJRAIVNSpaPDzzeZa5ttkq9LrjBkHnxw/tdrGjKdg2OPhQce8KEyt+m3owPmz4f77/f7qUVTRESqpJApqdLenn8eyqFD4aCDYPPNE6/SMmED4F575X+9piFz4kQYO9YHzGLa2vx+kyYlUy8REemxFDIlVQq1ZILfXst7M8OWzN69878ehszuutVjcdllmVFTgZmsSBv9uu7b1ub3FxERqYJCpqRKoZZM8NvrIWT26ZP/9Zq2ZD70UJeLszIz2ZWnuu7b0eH3FxERqYJCpqRKsZbMpqbaDgDqLmSG22sSMnNaMUPPs11Z+4uIiJRKIVNSJQ0tmbnTK4Vq2pKZO7FolrxDfIrsLyIiUgqFTEmVNNyTWZfd5fvuWzCdv8IWnTc0Nfn9RUREqqCQKamShpbM7rrLazLwZ/Togq2TbeRs79fP7y8iIlIFhUxJlTTfkxmG45oE4e23h/33zxs028m6oP37wwEHwHYF7tUUEREpkUKmpEqxlsw+fWrUShio65BpBjfdBAceCC0tOMtcxJ9zvq/cgAH+9Ztu6rqkkoiISJkUMiVVirVktrT4RWtqpdSQWbPFdHr3hltvhXHjaD/ksGWbH+erMGoUtLbCbbcVnuhTRESkDM21roBIOYq1ZA4Y4EOmc7VpiCs1ZL71VjL1ycsMtt+eu0bdDndnbb/jjppVSUREeia1ZEqqdNeS2dFRo9HbZN63u5D5058mU59inn++1jUQEZGeTiFTUqVYS2a/YIXEhQuTq0+2Ulsy64FuuRQRkbjV0a89ke61txduyWxuzuxTC92FzHoKdvVUFxER6ZkUMiVVOjoKtwiG4XPp0uTqk627FX/qqSVz0KDOz9V9LiIiUaujX3si3SvWkhlur9eWzHoKmWus0fn5U0/Vph4iItJz1dGvPZHuLV2a6RbPpZBZmrlz4fjjO2/TAj8iIhK1Ovm1J1KaxYsLd0eH4bPW3eWFppmsl5D573/XugYiItII6uTXnkhpFi0q3FJYDy2ZvXsXHlRTL4NtCl0/ERGRKClkSqosXlzfIbNYgKuXlszsltZ6mLNTRER6pjr5tSdSmlK6y2sVMou1skL9hMzseuy5Z+3qISIiPVud/NoTKU0pLZm1vCczDSEzO4T371+7eoiISM9WJ7/2REqzcKG6y6uVHcILDVISERGpVp382hPpXlsbTJsGa6+d//Vad5d3FzLrZeBPdsjUICAREYmLQqakxqJFvlxxxfyv17q7vK2tePdzvbRknnVW5rFaMkVEJC518mtPpHtheCw0GfvAgb6cMiWZ+uRKS8gcPz7zWCFTRETiUie/9kS6113I3GADX86dm0x9ci1cCP36FX69XkJm6A9/6Nxd7lzt6iIiIj1Pnf3aEyksDJmF1i4PW+XS0F1eq/tGs225ZeeWzHqok4iI9BwKmZIaYQgq1JIZbl+yJJn65ConZC5cGH99utPc3Lkls1bhXEREeiaFTEmN7rrLa92SuXBh8ZCZPbq8VnXs6Mg8bm7u3JKpkCkiIlFSyJTU6K67PAyfaegur1Vr62WXZR4rZIqISJwUMiU1umvJrIcpjIoN/MluyVy8OP765PPss5nHvXr57vIdd/TPaxV8RUSkZ1LIlNTo7p7MsKXw5ZcTqU4X3XWXZ6tVyMweQR5ex298w5dqyRQRkSgpZEpqdNeSGbrvvtirkld3Uxhlq6eQWet7WUVEpGdSyJTU6O6ezFpyztev1MnNa9U1nTvwJ7tUyBQRkSgpZEpqlNqSWQtheCs1AH/+eXx1KSa7JTOsq0KmiIjEoS5CppmNMrOrzOxpM5trZs7Mbi6w7/Dg9UJftxd5n+PMbKKZzTOzOWbWamb7Fdm/l5mdaWavmlmbmc0ys4fNbMcovm8pT7h2ed++ta1HPuUG4GnTYqtKUflaMsN7XT/4IPn6iIhIz1UvbUI/AbYA5gHvAxuVcMwrwH15tr+Wb2czuxQYHZz/eqAPcCQw1sxOc85dnbO/AbcDo4A3gKuB5YEjgKfM7FDn3P0l1FMiEt7HmD2BeK4DDoB3302mPtnCoFZqS+b8+fHVpZjskBnWNZwYftw4GDky8SqJiEgPVS8h80x8+JsC7AY8UcIxLzvnLijl5EHL42jgLWA759xnwfZLgBeAS83sQefctKzDjsQHzOeAPZ1zC4NjrgWeAa43s3HOuRp1fDaeCy7wZbGWzObm2nT7ltuSWasVf7JDZjga/2tf8+WqqyZfHxER6bnqorvcOfeEc26yc9l3jEXq5KD8RRgwg/edBvwe6Ascn3PMd4PyJ2HADI6ZBNwBrIQPoZKQ8eN9Wawls1ev2oTMUlsyn3nGl/UQMsN5O8Npl773veTrIyIiPVddhMwKrWZm3zGz84Jy8yL77hGUj+Z57ZGcfTCzfsCOwALg6VKOkeQ0FfnUNjdnAl+SSh35vu22vgzvL01a9p9xLS2+rMd7XEVEJP3qpbu8El8JvpYxs1bgOOfcu1nbWoDVgXnOuY/ynGdyUG6QtW1doBfwtnMuX7tYvmMkIcXmmKxVd3l3E8WHwlbYWrVkhiFzl10y0y2VOreniIhIOdIYMhcA/4cf9PN2sG1z4AJgd+CfZralcy4cWjE4KOcUOF+4fUjWtkqO6cTMTgJOAhg2bBitra2Fdo3EvHnzYn+PWhswYGcWLGjm/fefYc6c/Ely5swNmT9/KK2t4yN5z1Kv68yZfYEvMWXKG7S25vtbJqN3712ZPPl9WlvfLrpfHD77bEtgCHPnzqa19WUgDJ4jARL9DDXCZ7YWdF3jo2sbD13XeNTFdXXO1dUX/redA24u87hmYHxw7OlZ21cLtr1f4LjeweuLsrbtGGx7psAx6wevv1FK3bbZZhsXtyeeeCL296i1E090bpVViu/z7W87t9pq0b1nqde1tdU5cO6qq7rft6XFubPOqq5eldplF1/PkSM7b/dRM9m6NMJnthZ0XeOjaxsPXdd4VHNdgeddBJkuzfdkduJ8t/afgqe7Zr0UtjoOJr9w++wqj5GYOVf8fkyoXXf5WWf58pFHiu8H/r7NWtw3CjBzpi+7u44iIiLV6mm/aoJfobSEG5zvNv8AGGhm+SZpWT8o38za9hbQDqxjZvluKch3jMSsoyMzIrqQWo0uL2dehFqGzP/9z5dbbFGb9xcRkcbR00LmDkGZe7PbuKDcO88xX8vZB+enLHoOGADsUsoxEr+OjvptyQznmNxmm+73rWXIDP3sZ7V9fxER6flSFzLNbGsz61JvM9sTP6k7QO6SlNcG5Y/NbGjWMcOBU4BFwI05x/whKC8KpjQKj9kOv+rPTODuCr8NqUCp3eW1CHBbb+3LUsJbrYJwtnBkuYiISFzqYnS5mR0EHBQ8XSUov2RmY4LHnzjnfhA8vhxY38yew68SBH50eThn5U+dc89ln98595yZXQ6cBbxqZnfhl5U8Ar9U5Gmu82o/4JeUPAQ/4fpLZjYWWCE4phdwonNubqXfs5SvlO7yWq7406dPactK1kNLZnfXUUREpFp1ETKBLYHjcratE3wBvAOEIfOvwMHAdvhu697ADOBO4GrnXL7J03HOjTazf+NbLk8COoAXgUuccw/m2d+Z2VH4bvMTgNOAhcBTwEW5QVbiV0pLZq3uyVy6tPR1yz/4AP78Z7jhhnjrVEzudTzkELjnntrURUREeqa6CJnOr0F+QYn73gBU9OvZOTcGGFPG/kuBK4IvqbFS78lsb/eBNMnWuvb20tctr5X58zOPc6/juutqUnYREYlW6u7JlMZVand5uG+Sli6t/5CZPQI+N2T26pX8NRMRkZ5NIVNSo9Tucki+yzwNITP7muRex6YmhUwREYmWQqakRqnd5ZD8wJr29tLvyayVJUsyj3NbhBUyRUQkagqZkhrldJerJbOrYtdEIVNERKKmkCmpoe7y6nQXMgE++yyZuoiISM+nkCmpUU53edIhMw2jy4tdkzCc77hjMnUREZGeL5KQaWYbmdmZZvYdMxscxTlFcpXSXd6njy8XLYq/PtnKmSfz5JMz9UxSGDJvzl0Pi0x4D9c2FxERqVZZIdPMzjezj8xs+axtXwZeAi4FrgFeNLMVoq2mSGnd5QMH+jJ7TsgklNNd3tJSm2Udw5CZr57PPptsXUREpOcrtyXza8D/nHOzsrb9CnDAz/DrfY8ATo+meiIZpXSXt7T4MumQuWRJ6cGxVoNsioXMV19Nti4iItLzlRsyhwOvh0/MbHVgG+Aa59xFzrlTgXFk1iEXiUwp3eW1CpmLFkHfvqXtW48h849/TLYuIiLS85UbMocC2a2YO+FbMbPX/n4BWKvKeol0UUp3eS1DZqn3WdYqZIbzZOYLmXvs4ctNN02uPiIi0rOVGzJnAqtnPd8dWAJMyNrWp4LzinSrnrvL096S2a8fbLghfOELydZJRER6rnLD4MvAAWa2qZmtBxwBPOOca8vaZzjwUTTVE8kop7t83rz465Nt0qTSVxmqdcgsdO9o797JT/0kIiI9V7kh8zfAYOAV4I3g8WXhi2bWC9+F/nxUFZQacA4mTIDDDvOpranJl4cfDhMn+tdrVK16bMlcssSHsyeeKG3/l17y38t//hNvvXIVa8kMtytkiohIVMoKmc65p4H9gPuAe4FRzrlHsnbZEfggeE3SaMkSOPpof5PePffAggU+ES1YAHff7bcffXTnhbATUspclP37+3LhwvjrE1q8uLz9HwzuYL434X8l3YXMl1+GBx5IrDoiItLDlb1GiXPuUeDRAq89DWxVbaWkRpyDY4/1SWPBgq6vd3T4JsL77/f73Xpr9/3XEVq0CJZbrvg+YVdwki1yYTf5OeeUtn+vXv6YJIMwZG4hqPeViUREpGfQAB3JmDgRxo7NHzCztbX5/SZNSqZegRkzuh/BHQaoJBtaw5A5bFhp+4cDhJIOmYce6kuFTBERSULFv27MbA38SPO8Y2qdc09Vem6pkcsu8wGyFG1tfv877oi3ToFPPoGpU/1XMbVYu7y7buhctVpfPff9c511lubLFBGR6JQdMs3sq8AVwEbd7FriSs5SNx56qNOw5zkMYiSt3MrRbEzOotYdHX7/hMyZU9p+Zr47uhYtmaWuXR6Om+puEFNcCoXM3r2Tn/pJRER6rnLXLt8BP/H6EOBqwICngOuB/wXPxwI/j7SWkoycVsxH2ZuX2YoLuKCk/eNUaoADH/qmTImvLvneD0pvyVxvPV/26xdPfbpTqAX1ySd9+fe/J1cXERHpucptSzkXWAhs55wL1yd/wjl3MrApcBHwZeCu6KooiQmHZgcM3+TmKDC4J2f/OJXbtXznnfHUI5+wbqUG4Xvu8eUKK8RTn+4UmoFq/HhfPvNMcnUREZGeq9yQ+SXgAefch7nncN75+LXNL4yofpKkffft1IdbNGQ2Nfn9ExJ2f998c2JvWbJyu8tXXrnzcUlpaYGhQ2GLLfK/HtY/6XqJiEjPVG7IHAy8m/V8MdCSs8+zwK7VVEpqZPToTq2TYcjMq18/v39CwpBZqy7mYsrtLg/DXJIDf2bM8PdbjhxZeJ+w/gqZIiIShXJD5sfA0Jzn6+bs0xtIrh9VorP99rD//l26wbu0ZPbvDwccANttl1jVyh3BnaRyu8trEeZuucWXxSaAV0umiIhEqdyQ+SadQ+V44CtmtgGAma0CHApMjqZ6kigzuOkmOPBAaGnBgonWO8KPSVMTDBjgX7/ppkQnYg9bMgutu51t4MB465Kr3O7y8I6EJFsyBwzofp+7gjup11473rqIiEhjKDdkPgrsZmbLB8+vxLdavmRmk/AjzFcCfhtZDSVZvXv7lXzGjcN2+CIAjmDt8lGjoLUVbruttLQXoTCQlfK2J52UWcM8CeV2l4fTLCXZYhgOMgpHtufzpS/5UuuXi4hIFMoNmX/E32+5BMA59yxwGDAVP7r8I+C7zrmboqykJMwMtt8eO/sHALy6zkF89t48P/F6gl3k2cKWzFKCXO/etZmMvZxplnr1SraOYci85prC+4T3u44ZE3t1RESkAZR1h5tzbi4wIWfbvUCRO70kre6+25dvvw277Qavvlq7upTTXZ50gAvr1t2Sl9mam5NtyQzn2C82cCpc7vKVV+Kvj4iI9HzlTsa+q5mt1c0+a5qZRpf3AOFgEYB//7t29YDyBv6EAa7QfJBRW7TIl+WEzKSDcCmrDCV4i62IiDSAcrvLnwC+2c0+xwb7SYq9+GKta9BZOS2ZYRDNWiEzVosX+zJsCSxFrVoya7WUpYiINJ5yf+WU0tZhUGyCRUmD6dNrXYPOyhn4E4bMpFoKw5BZzy2ZYchUa6WIiCQljnaNtYHPYzivJCjBFSNLUs7An6QnOw8DeT3fk1lKdzn4EebDhsVfHxER6fm6/ZVtZufnbBpp+ZtDegFrAUcCWv045cJQVy8q6S5PKsR9+9u+LGdWp3ptyVx3Xb86kIiISLVKGV1+QdZjB4wMvgr5ADin4hpJXQi7gOtFuQN/so9JSjkDjZJuyVy40JfdtWQ2N2ueTBERiUYpIXP3oDRgHDAG+Eue/dqBT4E3nHMJDbmQuKS5JTPp7vKdd4ZnnoFNNin9mKRbMg87zJfdtWQ2N9ffz15ERNKp25DpnHsyfGxmfwHuy94mPVO+lsz29vImHI9SOEF4PXaX9+0LO+5Y3jFJt2SGumvJTHoiexER6bnKGvjjnDveOfdAXJWR+pEvZP7pT8nXI/Tcc74sZXBN0t3llYTvpFsyQ+ouFxGRpJS14k82MxsADMUP+OnCOfdupeeW2svXZTpnTvL1CB1yCNxzDwwd2v2+SXeXVxIya9WSqe5yERFJStkh08y+AfwI2LjIbq6Sc0v9yNeSmdQKOvkMHAjDh5e2b9Ld5R0d5Y0sh2RbMmfNyjxWS6aIiCSlrCBoZt8E/owf5PM08B6gX0k9UL7WrFqGzPb20lerSUN3eZItmdOmZR5315KpezJFRCQq5bY2/gD4DNjZOfd6DPWROhG2ZB5+ONx5p39cy5DZ0VF6kEtDd3m935PpnFYHEhGR6pS74s96wN8UMHu+sCWzX7/MNrVk5tfRUf6a4M3NMGFCPPXJlR0WS7knE2pzv6iIiPQs5YbMWcCiOCoi9SVsycxuoUtLS2bSQamSlsxZs5IbYFNOyAzvLVWXuYiIVKvckPkgfllJdaT1cIsX++D0xS9mtqWlJTMN3eVf+1py3dHZ79Pdz7BWqyWJiEjPU27IPBfoC1xrZgNjqI/UienTfXg66aRa18Qrp0s6Dd3lffokt3RnJSFT0xiJiEi1yh348zdgAfBt4GgzmwzMzrOfc87tWWXdpIZuuMGX5QSUOPW07vI+fXyQqySgVqO7n6G6y0VEJCrlhsyRWY9bgC0L7FfDOCJx6ajhivSVDPyZOhV23TW+OoUqCZl33eXLa6+F730v+jplyw7b6i4XEZGklLusZFOJXzVa4VrilJaWzPnzffnNb8ZWnWWWLoXJk8sPZeHclW+9FXmVuignZIYtmeE1FBERqVSCHXWSdmkZ+DNvXrx1yXbddb4MWybrUSUtmfvsE199RESkMShkSrfC+zJr2V1eTkvmwASHpFUaaE8+2ZerrRZdXQrJDpnd/Qzb2nw5eXJ89RERkcZQdsg0syYzO83MxpvZHDNbmvXaVmZ2jZltEG01JUlhEDn8cF/Ww2CQcloy997bl2GQi1O592KGRo/2ZRKBuJyWzI02ircuIiLSOMoKmWbWB3gc+C2wLvA5kD3b31TgBOCYiOonNRBOX7PFFr4MQ2Ytp7UppyXTzLcQJhGKKx0ZHq6kFLYcxikMmeuuC5tuWnzf3XbzdTvooNirJSIiPVy5vyLPBnYHLgSGAX/KftE5Nxt4CtgrispJbYTzN/bp48t6mDuxnJZM8EEpiQBXaUtmktMshe/x5z+XVt8vfEGjy0VEpHrlhsxjgGedcz93znWQf6qiqcBaVddMaiYMk2EL5uOP+3KVVWpTHyivJRNg+eX90o1xqzRkhoE5iftcw8BYzjyjCpkiIlKtckPmCGB8N/vMApavrDpSD8KQGbZkbrmlL9MyTybAsGF+1aK4VTriPgx8SbRkLlzoywEDSttfIVNERKJQbshcCAzpZp+1yL8KkKRE2F0etmQmvRZ4PkuWZOpTigEDMuEqTpXeQpBkyFywwJflhEwtKykiItUqd8Wfl4Gvmlkf51yXlZfNbDD+fsznIqib1Ehud3lTkx9Mk9QyjbmWLoWXXiqv1bB372TWBl9pJV/+4hflHReGzCRah8OQ2b9/afv37p1MQBcRkZ6t3JbM64A1gVvMbFD2C2Y2BBgDDAWujaJyUhu5A3/Ah6JatWT+61++fPnl0o8J1waPWxgWDzmkvOPCrv96bclUd7mIiFSr3GUlb8MHyUOBmcB3AczseeAj4EDgGufcw9FWU5L0/vu+zO6eXroUfvWr2tSnkmmCevdOJmSGLZHl1tEsudbhcJS9ustFRCRJZf/6ds6dgJ8L87/ASvh5MrcGpgDfcs6dFmkNJXFHHunLeln1pZ5DZhgSKxll3qtXst3l4dyc3VFLpoiIRKHcezIBcM6NAcaYWX989/gc59z8KCsmtfPZZ75Mcg3wYioJcEndk1lNyFy6FF57Ldr65LNggQ+YpYb13r0VMkVEpHpVrV3unGtzzn2ogNmzDBvmy759a1uPUCUtmUndk1lNyAR44IHo6pLPBx/AJZeUN5BHLZkiIhKFqkKm9Eynntq5BNhzT79UY1qkobs8CW+8Uf4xNb0n0zmYMAEOOwxaWvxfGC0tcPjhMHFi5ROTiohI4srqLjezt0vc1Tnn1q2gPlIHwpbD7Clv1lqrssAShUoGx4Rdvs75ATZxqfeQmT1DQKlq1l2+ZAkce6xv3l24MHPD6oIFcPfd8PDDsP/+cNNN5U2aKiIiNVFuS2YTfqBP7tdQYHjw1aeC80odyRec+vWDRYtqU58w8Bx8cOnHhBkk7ha5nhgyP/oI3nsv4Z+3c5mAuWBB1xFRHR0wfz7cf7/fTy2aIiJ1r9wpjIY750bk+RoKbAA8CrwFbBxHZSUZ4e/37ODUt2/tJugOg1x29313FDK9SkLmo4/6MlyzPhETJ8LYsZmh8IW0tfn9Jk1Kpl4iIlKxyFocnXNTgEOA1YGfRXVeSV4YnLIH3NRDS2Y5QS4MV40eMqvRXNHcExW67LLMhJ7daWvz+4uISF2LtFvbObcQeBw4KsrzSrLa2zOThYc++MBPCfS//9WmPlBe6FFLphe2Sl96afnHJhoyH3qoSxf5fAZwMT9iKTkXt6PD7y8iInUtjnsnlwKrlLqzmY0ys6vM7Gkzm2tmzsxu7uaYHc3sYTObZWZtZvaqmZ1hZgV/1ZvZfmbWamZzzGyemU0ws+O6eZ/jzGxisP+c4Pj9Sv3e0qq9vWtouvNOX954Y/L1qaQlMwyZcc+VGUXIfPHFaOqST1i/DTcs/9hEg3OeVsyfcz7ncjE38/WS9hcRkfoSacg0sxWBg4H3yjjsJ8CpwJbAByW8x4HAU8CuwL3A1fjBRlcAtxc45lRgLLApcDNwPbAafkL5vG08wfYxwKrB/jcDmwFjg/P1WB9/3HV0cS1b7MI8Uc68nWlqydxjj2jqkk8l9dtpJ1/OmBF9fQrKnsogcAdHAPA5y5W0v4iI1JeyQqaZnV/g6+dmdiPwJplQVqoz8YOGBhGshV7k/QcF524HRjrnvuWcOxsfUP8FjDKzI3OOGQ5cCswCtnXOneKcOxPYHD9IabSZfSnnmB2B0cHrmzvnznTOnQJsE5zn0uC8PdINN3TdtuWWvhw6NNGqAJmws0rJ7ePpuCfznHN8OWdOdPXJVcna6r/8pS+PSvKml333XVbJDowX2Jp3GA7Ae6zZed+mJr+/iIjUtXJbMi8o8PUT4LjgfBc5535T6gmdc0845yY7V9KcJKPw66Xf7px7PuscC4M6QNegegLQF7jaOTct65jPgODXKSfnHBM+/0WwX3jMNOD3wfmOL6G+Pca11/pyrbWSf+9Zs3y5wgqlH5N0S2YlqxKtv360dcmnkhBckwFeo0cva53sRQfb8sKyly7hh5337dfP7y8iInWt3F+Nuxf42g3fMriicy7OkeVhx+KjeV57ClgA7Ghm2R2rxY55JGefao7p0QYN8uV99yX/3mHoKWf+7STvyay0qzyJqR4racncc09fbr119PUpaPvt/UTr3XWD9+8PBxwA222XTL1ERKRi5c6T+WSBr6edc6855+JeJyQcvvBmnrotBabiVzFap8RjPgLmA2uY2QAAM2vBT8M0L3g91+Sg3KCSbyCtwtB2e967XuO1aJFvvCpn5Z6kWjIXLarvkFlJS2ZTk5/4Pu6A3omZX8nnwAMLV2rAAP/6TTfFu4yTiIhEotuQaWZ9ghHW/zSzgm1JwX7/NLPxxfar0uCgLHQXW7h9SAXHDM4py3mPHq+Wq/gtWlTeoB9IJmSOH++na6w0jA0cGG198sk3sX4pBg2K917RvHr3xt1ya96X3KGjoLUVbrtNS0qKiKREKTPhfR0/4GV/51zBX9nOucVmdgnwMHAMfmR2wzKzk4CTAIYNG0Zra2us7zdv3rwI32MkQKfzffppH2DHLtuTMHXq+jQ1rURr63MlH/Pf/w4BtmTixJdYvLjytFTsut5332qEDdqVXJOVVwYYyV57Tae1NZ4JSG+5ZW1gBC+//CJLl84t+bi5c9fjs8+G0dr6bCz1gvzXdvFiw99909kP176QfedP90FTior2/wLJpmsbD13XeNTFdXXOFf0CHgQmd7df1v5vAA+Vun/OsSMBB9xc4PVJwevbFHj9teD1jbO2zQy2rVDgmHnB6wOC5y3B888L7L9i8PqMUr+vbbbZxsXtiSeeiOxc667r3DHHdN72ySfO+Q7eyN6mZN/8pnNrrlneMU8/7ev6979X997FruuNN1Z/TVZbzblvfavy47sT1m/8+PKOO/dc53r1cq6jI556OZf/2ra0ZOqc/fWTn8RXj54myv8LpDNd23jousajmusKPO8qyHG5X6Xck7kV0FpGbn0KP6VQHN4Iyi73Q5pZMzACPxn82yUesyo+VL7vnFsA4Jybj5+vc2Dweq5wTHCXezx7iqVLu672Eo7s/trXkq/PrFmZgUelSqK7PIo5Q3v3jv++USh/9Ptyy/n7OZNer37+/Pzbk66HiIhUr5RfPSsC5UzLPAMoY7KZsowLyr3zvLYrMAB4zjmXPQlLsWO+lrNPNcf0GPlCJsBmm5V/b2QU3nkHRowo75gkQmYl0xbl6tMnmQBVbiBeLpj//PPPo69LJRQyRUTSp5Rfk21AOUMUBgJx/Uq4C/gEONLMtg03mlk/4KLg6R9yjrkRWAScmj2BupkNBc4Lnl6bc0z4/MfBfuExw4FTgvPVYIHFZBQKmb16ZUYrJ6m9PTO5eqnC/eMcIZ2z1HZF1lrLh+i4lXv96i1k5q5AJSIi9a+UgT/vAdt2u1fGtsC7pe5sZgcBBwVPwzVdvmRmY4LHnzjnfgDgnJtrZifiw2armd2OX4HnAPxURXcBd2Sf3zk31czOBn4HPG9mdwCL8RO7rwFc5pz7V84xz5nZ5cBZwKtmdhd+6cojgOWB01zWxO49TT2GzHJbDZNY8SeKc6+8Mrxb8r+W8mSH4CFDyju23kJmLT53IiJSnVJ+dbfiQ1+3QdPMtsEPQX6ijDpsiV8t6Dhgr2DbOlnbRmXv7Jy7Dz/89CngUOA0YAk+EB4Z3LBKzjFX4YPof4Bj8aO+pwPfDANsnmNG41f1mR7sf2xw/P7OuavL+P5Sp95CZkdH5SHzH/+Ivj6hMGReeGHl52hpgXnzoqlPruyfVbn3tIb7zy19QHqs8i11KiIi9a2UX91X40dT/83MNi60k5ltBPwNv674NaVWwDl3gXPOinwNz3PMs865fZxzQ51z/Z1zmznnrnDOFYxAzrmxzrndnHPLOedanHPbOef+0k3dxgT7tQTH7eace7DU7y2t6jFklntPYRgyb4zxpoYwZH43dyHTMgwcWHiwS7Wyf1Zhy2SpBgezxSY+V2Ygd2WpKG5NEBGRZHXbXe6ce8PMfo5fo/yloOt4HPB+sMvqwJ74VsW+wPnOuTfynUvSoVDIbG5OX0tmnMKQWc3c4GFLpnPRL2IT/qwuvrj8c4chM+mWzLXX9mu6H3igvzZxBXAREYlfKfdk4pz7uZktBX4GHA0clbOL4busf+yc+1W0VZSkFWvJrMUAjJ4cMgcO9N9fuHRmlCpd7QcyS4i3tUVXn1I0N4eT1GvlSBGRtCspZAI4535pZrcAJwA7AeEckh8BzwA3OucSGCcrcXLOt4DlC5mvvQaffJJ8naoZ+BOnMHBXGzLBt2ZGHTIrWbc8FNZl0aLi+0XppZfgrbdgR7+wFGuvDf/5T+b1Qn/8iIhIfSrrv+wgRP4sprpIHQiDSb5f5rUImFBZS2YSy1uHLZnVBJ+WFl/Omwcrrlh9nbJVEzLD+VCTnJ9y6619GV7Pxx/3q0gefbR/vnixQqaISJpEMJ209CRh61y+X+bhaj9JD8KoZOBPFKvxdGfJEh9+q5mUPWzJjOPew7S1ZIbC+q66Khx1FFxxhX/+0UfJ10VERCqnkCmdFAuZu+7qy6SDRyUtmQDbbZdpHYvDkiXVt5hmt2RGrZqQ2bu3vyeyFivtTJ/e+flbb/nyrLOSr4uIiFROIVM6KRYya9W6Vck9mQArrRTv4JEoQma9tmSa+S7zWrRkPpgzSdg3v+nLcieUFxGR2lLIlE5KCZlJt25V2pIZ95RLUYbMOFoyx43zZaX30vbuXZvZBMJ7MEPbbOMDb5wT64uISPQUMqWTf//bl5Mnd32tliGzkta45uZ4Q9LSpdF1l8fRkhneyzh+fGXHm8Ff/+pnHIhb9ntstVXX1xctgg8/jL8eIiISHYVM6eSpp3z5Rp7p9NPYkhlnyFyypPrRznG2ZIbBvNLW3LlzYeZMuPnm6OpUSPZgsmLBPYnAKyIi0VDIlE622MKXF1zQ9TWFzM6iHPgTdUumc/DCC/5xtbMBfPpp9fXpTqlB+DvfibceIiISHYVM6SQMZWELW7ZahcxKB/7EvUJRPY8uHzMm87ja+1KrmaKpVNk/p2Kh+Prr46+LiIhEQyFTOqnXgT/1eE/m3LnVh8zevf3XggXR1Cl0wgmZxz+rcvmEJOYczQ7C6hIXEekZtH6GdFKvIbOS1rQXX4T334e334Z11om2TtOmwUMPRXOuPn38ajZx2Wyz6o5PoiUzO2SOGBH/+4mISPzUkimdFAuZtVhqECoPma+84stHHom2PgBvvhndueKej7La1tYkQubbb2ceH3xw19erbY0VEZHkKWRKJ8XW465FS+bSpf6rT5/Kz3HqqdHVJxTl0pp9+tR3yIxzrtHQeef5crvt8r8+c2b8dRARkWgpZEonYUtmvmBSixV/Pv/cl4MHJ/eepYjyvsHp0+Md0FJtS2S+QWBRC+9JLRSIR43yZdS3PYiISHwUMqWTersnc+5cXw4alNx7liLKlszQa6/Fc29mtUtr/uEP0dSjmKef9mWhFuvdd4eRI2GNNeKvi4iIREMhUzpRyCxNHCOgN9sMzjgj2nM++2zlx4arP1W6YlAlit0WEfdsASIiEi2FTOnk9tt9WS8hc84cX/bk7vLf/CbzeOLE6M5breWWS+69wm7wAw4ovM+iRfDcc9FP9yQiIvFQyJROwpBTL6PL67UlM8ru8qOOyjwOr3E9SHK+yvnzYa214HvfK7xP2KV+5pnJ1ElERKqjeTIlr3whs6nJd2cqZEYbwLIH1tRTyBw+3Ldexz0Z+xtvwIwZ/nEp94++9Va89RERkWioJVPyKhQs+vWrTXd5JSHz8MOjrUu2KFsyw6UlAZ54IrrzQvVh+MQTq58CqTuzZ5e3f5KzG4iISOUUMqUsSYfMsCWzknsyi3W9VivKlsy4Q1w1Wlp8CLz77vjeo62ttP3WWsuXSS8GICIilVHIlLIkHTJnzPBd9NmtfaXafntfVnJsd+JcBjJK1U5fNGCAL8N5KuMwdmxp+229tS8VMkVE0kH3ZEpZBgyAefOSe7+pU2HddSsLS/37w957w6xZ0der1Na3WoiylTWOgJ7r8stL2y+8X1Xd5SIi6aCWTCnLqqvChx8m937t7dUtKdmnTzytjuG9ovUoey7JagNn2JIZp1Lvtw1XLlLIFBFJB4VMKctGG8FLLyU3vU17e3XLIvbuDS+/HH19zzor2vNlu+666o7PDmFrr13duZJoydx55/L2V3e5iEg6KGRKWdZYw3cVL1mSzPt1dFQ3hU44YOVvf4umPqGhQ315wQXRnhfgO9+p7vgwZF50UfXLMCbRklnqZ2mrrXy5yirx1UVERKKjkCnLlDItT7jqT1L3JHZ0VNeSGXrjjerPkW333WGTTeBnP4vmfE89Fc15IBMyV1yx+nMl0ZJZ6lKRo0f7crvt4quLiIhERyFTlglblNZfv/A+SS8tWW13eSjq7vJqW1hzldtlXEwYMqOY2L2eWjKbmmDECN2TKSKSFhpdLsuEv+xPPLHwPv37+zKpkPn447DCCtWfJ8qQ6Rx8+mk04TeUO3reucqnHwpDWDUDpkJJtGSGn7t//rP7ffv2hZkz462PiIhEQy2Zskw4CrtYOEmyuzycKunTT6s/V5Qr9Fx/vV9He8qU6M6Z6+GHKz82/DnWe0vmwoVNjB4NEyb4LvA99uj+mKVLM2uYi4hIfVPIlGXCFqViK9CEryU18CcqUbZkPvSQL+OcL3S//So/Nsru8nxr2Efl7rvXWDZH5qRJpR0zYgQsWBDtHw0iIhIPhUxZppSWzPA+xPb2+OsThbDLOcpQUu0qOnG77z5fljqgppjslsyof+bt7eVfyMcf9+Wdd0ZbFxERiZ5CpixTSktmeB9iEi1JUYSasDUvygnZ6z1ktrb6ctttqz/X6qv7kfQAn31W/fmyNTVV3rysuTJFROqfQqYss2CBL8PBPfkk2ZIZRZBN81KEJ51U2XHPPefLaufIDB17rC/nzo3mfKHssD5yZGnHjBjhyyQGJImISHUUMmWZTz7x5UorFd4nyZAZxXuEQW3gwOrPFYqiG7oU119f/jFxLKEZ3pd5zDHR3tua3ZJ5/PGlHXPTTb6M4n5TERGJl0KmLBN2h4ar2eSTtpbMiy7yZZSjpJNsFS13gNWsWb786U+jq0MYMsePz3TFRyG7JbPU6ZYGD/ZlGlumRUQajUKmLFPK1DdpC5nh/aVRtvCFg0+ScMYZ8NZbpe//3//6cpNNoqtD9qTzH3wQ3XmXLs2kzGL3AWdL8+0PIiKNRiFTlgm7gYtNW5O27nIz30oWRzdyEq65Bg45pLR9f/1r2HNP//jzz6OrQ3bI/MY34PTTM/fvVuOGG9ZZ9lghU0Sk51HIlGXKCZlJjC6P6j369El3KCk10J17bubxV74S3fvnLp/5u9/BVVdFd34oPtgsm0KmiEh6KGTKMmHILGUKo7S0ZEK6WzKh9LCdPShn7bWje/98f3RE/fMvZbUfiGdKKhERiYdCpiwTDjIppSXzhRfir0+ULZlpCCXDh+cfdFXr1W3yrdFe7VyhM2Z0fp7bWlpIGDK1tKSISP1TyJRlSukuD8PFeefFX596D5mbbRbt+aZO9aPDjzmm8/bu5qdsb4drr422Lrnnj9q++1Z2XBgy77knmT90RESkcgqZskwpITPJNcvDcHPhhdWdp08feO+96uuT66WXoj8ndG0lXH/94vuPGQPf/W7m+Xe+E2198v3Mq23JfOedyo7LbvG8667q6iAiIvFSyJRlnn/el8VCZpIDLsIQEc6NWKkpU2DcOHjsserrlK3ULt5y5Qa4FVfMPP7tbzMr+oQmTOj8POqWxzhCZvbn6Morqz+HiIjUH4VMWeaWW3xZbOBPkr/Y//hHX778cjTne/31aM7T3Bzv7QLbbNP5+dy5/uutt+DMM2GnnTq/nrsyUNQhM9+tBlGGzP33r/4cIiJSfxQypYtiLXTLLZdcPaKeSL3UuRiLWbLE31YQ5QpCuU47LROwwQ9yGTwY1luvtOOTaMmcORO+9KXyJorPlv0zLffe2yFDup5DRETqj0KmLLPRRr4s1l2+yy6+3Gef+OsTLjVYT2EinLMyzpDZ1JRZc70SUa+tvs46XbddeqlfZjKKAUflhsxBg3zZ1lb9e4uISHwUMmWZ9daDrbfufr/ddoM5c+KvT7g04hZbRHO+KFr4wmBT6uTh1Zg0qfBr2S2duaJuydxtt8KDr4otQVqqcrvewxbpW26BTz6p/v1FRCQeCpmyTHt7/jkRc62+enT3SRaz666+zB45XY3TT6/+HEm0ZIa23bbwayefXPi1OKYc2m+//NsruS/yoosyj1dYAdZdt7zjs1czqrS7XkRE4qeQKcs88gj85z/d73f//TB/PvzrX/HWJwxLxbrvS/HQQ9XXJZRkyOyOme9q3njjztuj7i4H38K9885dt1cSMs8/P/P4nHPKb8n83e8yj8NbKkREpP4oZAqQGdxRyn1u8+f78rXX4qsPZEJmtVMFRXn/aBgyk+guB7jkkuKvH3VU11HzLS3x1OXOO+Ggg/yAn9BVV5U/xVT28peV/AGRPYArisFcIiISD4VMASpbXafaFsbuRBUyo/T5574MB5/E7Qc/KP76nXd23Zbd0helVVeFe+8tf0WiYg48sLo61XrJTRERKUwhU4DKuljjDn9hneopZL76qi+TnMqpHMcfD8svH+97rL125ce+/Xbn5yNGVHaesOs+qkFhIiISPYVMASobLFLthNzdieqeTMiEkmqXlzzrLF8m1ZJZruyu6LgMH9512+9+V9qSo1ttFU0dzj038ziJ71lERMqnkClAeSHzT3/yZVz3/oXCOpUy4r074So5I0dWfo6nn848TrIlM3fZyGKSCFybbtp12+mnw4kn+vd///3Cx2Z3re+554yK65D9h0cp4VZERJKnkClAeSFzxx19Gfcv9/b26LrKwwEiud215QinVIJkWzK33x4OPzz+luNyZF+L0F/+Ag8+CGuu6Wcq6M6mm1Y+2Wp2yMwerS4iIvVDIVOA8kJm1Ms9FvLLX0Y352PUg5SimIS8HHfcAU891f1+a60Vf13Ah8h8c3U+9pgvr7jCh+Lx4/29tTPyNFo2N1fe7Jo9qvzXv674NCIiEqOYxwdLWpQT5sKBJR99FE9dQlF2/cY9Ej4JxW5PuOwyf6/k/vsnU5cBAzLLkGb7/e99+fjjvnzwQX97xQ03dL7dAKBXr8p/wD3h5yki0tOpJVOA8kPm4MHw4Yfx1Qdggw2iO1d2KLnjjurOdfDB1R1fqS237LptjTV8ecgh/ivJeSNLmZC+udkHTMisex9aa60Fkb13vqmcIuecv0H2sMN84m9q8uXhh8PEiRqBJCKSQyFTgPK7pfv2rWy1l3KssAJ8+cvRnCs7fB15ZGZS9Uq89FL19alEvnsyV1nFl7WY5qlfv+73KVavTTapfILNIUM6Pw+DbGyWLIGjj4Y99oB77vEfIOd8effdfvvRR2sUkohIFoVMAcoPmX36xB8ylyyJbtnA3O7Vs88u7/jsSb+Tvh+zmHvvhSuv9INtklbKQKRCk6VXu9Rn7lygsQ6Kcg6OPRYeeMCHytxvqqPDL4N1//1+P7VoiogACpkSCEPmsGGl7Z9ES+bixdF1/+aGzFJGP2f7858zj+slQ6y8su8u//73a/P+pfxsLrgg//ZqQ+Fyy8Hmm2eeRzHNVUETJ8LYsd03f7e1+f0mTYqxMiIi6aGQKUBmdZ0rryxt/yRC5pIl0YXM3POUGxTnVD7bTiw++gjefLO2daimlbmUrvbuvPxy5nGsLZmXXeYDZJaPWYk28nwTbW1+fxERUcgUr9x1wpMKmVF1l+dOnl7OmtevvNL9GuJJW2UVP/iqlvbZp/Jjq5kUP2SWmbN1/vzqz1fQQw91+sA4YBgfsxtPdt23o6P6ewFERHoIhUwB4JNPfFnqSjZ9+sQ7T+bcuTBlSnTh4bDDOj+fPj3/3I35XHxx5vE++/hxH7Xy3//6r3rQt68fBPXVr/rn665b2h8eq68eXctjOCn80qV+qqTp06M5bydBK+aPuJiv81dW5wMAJrE9t3EkW/JS3v1FRBqdZpsTAN5915frrVfa/nG3ZF5+uS/vvz+a8/Xq5bvMw8G/ixf71sBSus2zB0Xdc09tB/5svHHt3jufLbf097eecgqcemqyUygB/N//+T8Chg/3y1rusAP8618Rv0n//rBgAb/hR11eOprbAFhIX/qxKLO/iIiktyXTzKaZmSvwlbc9w8x2NLOHzWyWmbWZ2atmdoaZFewkNrP9zKzVzOaY2Twzm2Bmx8X3ndVGGL5KDVBxh8yZM6M/Z75AObeEWXSyQ2ZU3fc9SVMT/OEPsMkmpbVQRj3J/nLLweef++cffxzduZfZd99uRxbNIbh3oanJ7y8iIqlvyZwD/DbP9nm5G8zsQOBuYCFwBzAL2B+4AtgJOCzPMacCVwGfAjcDi4FRwBgz28w5V2d36lUu7PouNUT17ZvpYo+zPqVM+F2N/fbrfrnGcFAU1Nf64Wly2WV+NPj99/sWx6g98IAvY/n5jB4NDz8MRW7deImt2JvH/Iim0aNjqISISPqkPWTOds5d0N1OZjYIuB5oB0Y6554Ptv8UGAeMMrMjnXO3Zx0zHLgUH0a3dc5NC7b/HJgEjDazu51zUXfO1UQlITPOlsywtSvK7td8LWi5Sx3mE9X66Y3iL3+Bt9+GUaNgs838tpNP9n8wRDW5frawFRNiCpnbb+/X67y98C6fsoLvJj/gANhuuxgqISKSPqntLi/TKGAl4PYwYAI45xYCPwmefjfnmBOAvsDVYcAMjvkM+GXw9OS4Kpy0ckNmEpOxh+8TlULdtGut5cPmzJl+SsRQWxvsvvtIDRYu07HH+vkxN900sy2p2xRjCZlmLLr+pqK7LOwzCA48EG66Sc3dIiKBtLdk9jWzrwNr4TuzXgWecs7ltj3tEZSP5jnHU8ACYEcz6+ucW1TCMY/k7JN6lbRkxjm6PAyESdwD+d57vofzk09g6lT/3m1tflnLbCutFH9depqOjmQz1/z5cMUVcPrp0U7QvmBJ8Sb1haefA79ZO7o3FBHpAdLekrkK8FfgF/h7M8cBk81st5z9NgzKLtNXO+eWAlPxgXudEo/5CB9q1zCzmO8aTMbixT4M1Ms8meFApCjX5D7kEF9ecknX1yZN8gEz9N57XWeiqeXURWmVdKPehx/CWWf5WyiTtHCYAqaISK40t2TeCDwN/Af4HB8QTwVOAh4xsy85514J9g2nrS60bku4fUjWtlKOaQn267LenJmdFNSFYcOG0draWvy7qdK8efOqeo8pU9ahuXkNnnyym1EwgZkz12P+/FVobX2m4vcs5v33vwCszKJFC2ltHR/JOU880Tj88N6stNJiYGTB/VpbW5k+vR+wQ6ftS5e2EvOPsaFU+5nNGNlly4QJ/2XgwGiGmn/+eTNz5/YGvthp+5lnvsEVV/i/RV9//W1aW9+N5P2qFd11lVy6tvHQdY1HXVxX51yP+sIP1nHAvVnb3gy2rVfgmGeD17+UtW1xsK25wDEfBK+v2l2dttlmGxe3J554oqrjzzzTueWWK33/s892rl+/qt6yqBNPdA6c+81v4jm/7xTP//XlLzv35ptdt0u0qv3Mhl57revP6txznZs8OZLTF/yczJ6defzTn0bzXlGI6rpKV7q28dB1jUc11xV43kWQydLeXZ7PtUG5a9a2sDWy0EJ84fbZFRxTZ6taV2bx4vLufwwH/kQ552G23r1h+eXh7LPjOX8x//gHTJuW/PtKZTbZpOu2X/0K1l8/3vcdODDzeOFCP6JeREQyemLIDKfxbsna9kZQbpC7s5k1AyOApcDbJR6zanD+951zXbrK06jckNm3rw+Y2XNIRumTT7oOvElSuFSiSCG9esFdd/nHt9zil9V85JHix4iINJKeGDLDG+myA+O4oNw7z/67AgOA51xmZHl3x3wtZ5/UKzdkNgd38776ajz1mTHDL/soUs8OPdSXH37oy1deKbyviEijSWXINLONzawlz/bhwNXB05uzXroL+AQ40sy2zdq/H3BR8PQPOae7EVgEnBqcNzxmKHBe8PRaeojFi8ub+Px///PlYV3WSYrGjBkwbFg8565Ec5qHyDUwM9hnn+jO9+c/w0svwb33Fn4/ERHx0vqr8wj8ijtPAe/gR5evC+wL9AMexg8AAsA5N9fMTsSHzVYzux2/ks8B+KmK7sIvNUnWMVPN7Gzgd8DzZnYHmWUl1wAucz1ktR8ovyVz0CBfLojpZoGPP4aVV47n3JV4s8tEVpIWUXZhv/UWHH88bLll/tejnJtTRCTt0vpf4hPAg/hgeTRwFrAb8AxwHLCfc67TVOHOufuCfZ4CDgVOA5YExx4ZjKYi55ir8EH0P8Cx+CmJpgPfdD1o3XLwLZPl/II8L2jLPf746OviHMydmwmycfjkEzjjDL9Ay3PPFd/3j3+EESPiq4ukR74W7SOOyDxWS6aISEYqWzKdc08CT1Zw3LNAWZ1nzrmxwNhy3ytN2trgP/8p75iwlTGO5QIXLfIDipZbLvpzh1ZYwa8ME5o8ufBo5JNOiq8ekoz99oPDD/dLXpZiyhT/b+LjnKk2w0UCst1+O9wR9IP07VtdPUVEepK0tmRKhD7/3JfZa013p1cv32qT75duVPWJM2Tmyr0f9Sc/gb/+FY4/fmr+A6SujBsHaxdZdOehh+C440o/3/rrw0EHdf0Do9DnPVxN6vvf91NviYiIQqbg5/gD331cjj594lm/vB5C5nrrwde/Dsce+05ylZCK7b47nHNO9/uVMq/r/PmFXysUMv/2t8zjzz7r/j1ERBqBQqYsW6O73K7v3r17TkvmaqvB9dfDbsGq9+WMtJf6EH5eVlqp8D5NTTB7dvHz7Lln4deGDy98XhER6Uz/NcqylsxyQ+a8eXBtDJM41SJkAnz72/4+zdVX12TsaXTkkXD55Z1bFfOZNKn46xMm5N/+2GNw2mmV1U1EpBEpZMqylsx+/So7tr092vrMnevLOEeXF7LVVvD++7Diism/t1SnVy8488zuB9989avdB81Cx6nFUkSkdPovUypuyVxjDV/+7nfR1mfWLF9qAIVUopQ/erbfPv/2QvNffv/7FVdHRKRhKWRKxS2ZYRh8991o6xOet5Zrl0t6ldOy7hy8844fxHb99YWXhbzyymjqJiLSSBQypeKWzDCULl0abX0+/dRPjzRkSLTnlcZQ6lyVM2f67u/hw/0AtmrnQz355OqOFxHpaRQypeLR5eEv8zCkRmXWLB8we/WK9rzSGLbfHq65Bi65BC67rPB+pS5bWupMA1dfXdp+IiKNQiFTKu4u32EHXw4YEG197rzTj/AWqYQZfPe78IMfwFlnVXeu3/8e/vvf0vbt1Qt++EP/+KOPqntfEZGeQCFTKp4y6LrrfBkOAIrKxx/rfkyJzne/W/mx3/qWn5i/VOFgtdVWg6FD45lHVkQkLRQypeIpg8J7JsOW0CiE93cWmxBbpBzXXANTK1wdtE+f8vbPngh+9uzMIDYRkUakkCnMmQMtLeXfA9nc7L+ivCdz0SJfVjJnp0gha68Np5zSdXu4wlMhZuW9T+5gtagHxYmIpIlCpjB/PgwcWNmx/ftH25IZBlaFTImSGZx7btftt9/edduvflX5++S2fEY9KE5EJE0UMoUFCyofvBN1yKxm9SGRYlZfHZ59tvO2VVbxXdqbbZbZFrbob755+e+xzz6+1TQU5b8NEZG0UcgU5s/33eWV6N/fh9SoVDpnp0gpdtwRXnut87ahQ+Gpp2DXXf3zXr3gs89g/Pjyz9/U1Hmg0aOPVl5XEZG0U8iUqloyBw2CefOiq4u6yyVuG27YdduQIbDFFv5xr17+eaV/6Ky2Wubxxx/DtGmVnUdEJO0UMqXqkBmOTo+CQqbErVcv2HlnuOuuztvD5SirXQTgyCMzjy+5BEaMqO58IiJp1VzrCkjtzZ8Pw4ZVduyQIfDQQ/7esyi6uBUyJW5m8PTTXbc3BX9yl7rCTyHVHi8i0lMoZEpVLZnh5NO77gqTJlVflxtv9GWzPpmSsAsvBOfguOOiP7dz5U+HJCKSduouF15/vfJRsLvs4svnn4+mLn/+sy/D+TJFkrL88n798Sha0Z94ovNKWK+/Xv05RUTSRiGzwU2f7ss5cyo7Pup1y/fd15df+Uq05xVJ0siRsPHGmeebbBLtADkRkTRQyGxw//ufLytd33nUKF/ut1/1dfnkE39/J2TujxNJq8cf7/z8lls6P7/mmsy/PxGRnki/yhvcddf5crnlKju+b1/YaKNoBv1ccUX15xCpF3vs0fn5okVw330werS/R/OUU2DbbWtSNRGRRChkNrittvJlOBF1JVpa/Ah1Ecl4+OHOz8ePh4MPhssvhyVL/Db9uxGRnkwhs8F1dPiymmlXogqZues+i6RZ376+xTJ0222Zx4sXx/jGzsGECXDYYf4fZ1OTLw8/HCZO7FwpEZEYKWQ2uCgmoI4qZGp+QemJDjqo67bYBgEtWcLGF13k++rvucfPT+acL+++228/+uhMU6qISIwUMhtcGDKrmZdy8GC/1nO1wqB72GHVn0ukXtx7L5x6audtq66aeXzxxRG9kXNw7LGs+OyzPlSG3RShjg7/1+D998Oxx6pFU0Rip5DZ4JYu9WU1o7lHjIB33skE1kqdc44v//KX6s4jUm+KTcl17rnw7rsRvMnEiTB2LL26m2S2rQ3Gjo1m9QQRkSIUMhtce7tvQaxmNZJhw3xYffbZaOqkJSWlp+mup+C99yJ4k8su67SqwmJ6M4uhLKRv133b2vz+IiIx0uJ9DW7p0uruxwQYOtSXu+1WeQ9cds+elt+TnmbrrYu/vvPOEfReP/QQV3acSgdNHMMtDOPjZS85cv5RdXRkJqUVEYmJQmaDa2+vfp3wlpbq6/Gf/1R/DpF6tcoqPtdFtcjAL3/p/7jrtIhCWxtncCUAN/CtTvv/m03ZjNc6n6TStWRFREqk7vIGF0VLZvYvzkonVP/gg+rqIFLvzGCddQq/fv31cOutpbVo/vjH8L3v+cfz5/tzm8t0B/yHTTvtfxejup4kihUURESKUMhscFG0ZGaHzLPOquwc++9fXR1E0uCtt2CXXfK/dtJJcMwx/t/TnDm+oTFf6My+f/Oss2DgwO7ft0t3eVMT7LtveZUXESmTQmaDi6IlM4p7KMNR7iI9XSn/XoYMgSOO8KHzlFM6v/Z//5d5XGrPwcPswyPszUes4jf06+fXtxQRiZFCZoNra6s+ZA4e3Pn5L35R/jk23ri6OoikxVprlbbf2LG+/MMfOm9/5pny3/MFtmUfHmE3nvTd5AccANttV/6JRETKoJDZ4F57DTbcsLpz5K57/pOflH+OKCZzF0mDa66Bm26Cp58u/9hJk+D11yt/78lsAAce6CugaRxEJGYKmQ1u5kxYe+3qzlHp76rXX89M1ffZZ/CNb8DUqdXVRaTeLbec/6yH0xbdemv3x7z8Mrzyil8ZslR/3OaPebcv/ettWsNVRBKhkNngPv0UVlih+vPkdul95zu+fOEF2GsvWLy46zE77AA/+IGfvmjRIvjCF2D48OrrIpImRx0F110HN99ceJ+ttoItt4Rf/7rwPuefn3l8yilw0vPfYfbsrvvNnVtpTUVEyqOQ2cAWL4Z586IJmSef3Pn5ddfBG2/4wQt//3smdIJvtXQu88tu02C2lY8/RqQhnXiiH+TjXPcTt+ez5ZZw4YWZoDpsmC9z75cOX/v884qrKiJSMoXMBjZzpi+jCJkA553X+fmJJ/opWwDGjPHLTu61Fyy/PJx+etfjjzkmmnqIpNn48eXt/9xz8NJL/vFRR8HZZ/+Pc84pvP/SpTBoEDz8cOV1FBEphUJmAwtbH6Oak3mzzTo/zx3YsN9+vlUT4Kqruh6/zTbR1EMkzXr3hokT4Utf6n7f8eM779fUBPvsM73TLZerrpr/2CeeqK6eIiLdUchsYA8+6Mt890tW4ogjir+e7/4wEelqu+3gySe7bj/pJN+lfvrpsNpq8MUvdn+uqVNhwYKu2zW4XETippDZwHbbzZdf/nI05zPzS9yJSPV694bDDvOPzz3XD9K79lr//Le/LX0p1r59fW/F5Zd33n7JJb4888zMPZwiIlFSyGxgm2/uBwYUW0+5XAMGlN9iucMOmidTJJ++fX35hS/4e5mraX3Mt3Srcz6watCdiMShylWrJc0WLozufsxsgwb5X4a5ay4Xcu21fhk9Eenswgv9WuX5AmK51lyz67Z8A/BERKKilswG1tYWT8g0g46O4vvssIMvv/Md2GKL6Osg0hOssw60tuafiqhcffvCb37TeVv2ALw33/TlF78IBx1U/fuJiChkNrC4QmYo/KUFcOednV8LR5IfdVR87y8inZ19duHXNtzQ/58wcSLcf39MFXAOJkzwN5u2tPjh8C0tcPjh/o1L7f4QkVRQd3kDiztkrr++v9fyk0/8Sj5nnOFHzI4ZA+uuCzvu2HXdcxGpnRNPjPHkS5bAscfCAw/4e3XC7o4FC/x6mQ8/7O8LuOkmLXsp0kOoJbOBffgh9OsX73sMGQLrrQfNzXDFFfDii37AUUsLHH20plERSdp++/mvfG65JfP4lVcifFPnMgFzwYKu99N0dPipKe6/3++nFk2RHkEhs0FNmQIvv+xX4RGRxjF2rP/qzpZbwk47+Xk2qzZxon/TfBN2Zmtr8/tNmhTBm4pIrSlkNqj33qt1DUSklv785+73ee45P/jotdd87nvyyQrnwr3sMh8gszhgFkO77tvW5vcXkdRTyGxQ8+b58vrra1sPEamN44+Hp57yk7Hvu2/xfTfbDLbfHkaOhIEDob299PeZMgUefKAdOjq4ku+zCh8xkxX5I99hBWbxJut3PqCjAx56qOzvR0Tqjwb+NKhZs3y5xx61rYeI1M4uu/gv8L3UBxxQ2nHNzf5+7oUL4Xvfgx/+0C9Pu36QF19+2Y/lGTUKDj4Ypi66m3dYizO4EoBTuZo5+HmZprAeGzC58xvktHqKSDopZDaoTz/15fLL17YeIlIfsid833//7u/bXLjQl9dc478AjjgChg7NLH950UWZgeJr8+6yY+/kiGWP/8Jx7MLTLMe8zMnjnPZCRBKjkNmgZs2CXr2imeRZRHqGmTP9PZdDhsAvfwn77OMnZy81891xR9dtS5YUP+ZOjqCdXgxiLrNYnvuaDu2+/15EUkEhs0H98Y/+vipNISQioRVX9F8Av/51Zvull8IPfhDf+97NqMyTfv1g9Oj43kxEEqOBPw1owQI/QbqISClGj/ZTGT3+ePzv9adNrsBtu138byQisVPIbEBz59a6BiKSNsOHw5e/nP+155+P7n1OnHQSw0eoi0WkJ1DIbEAKmSJSqYkTfbnXXn6985/8BLbaqut+b72VeTxjBmywAfxgtOOte17h62s9xZP99lr2+saD3u907Lvvwve/7yeEf/rpGL4J0DrqIgnQPZkN6PXXffmFL9S2HiKSPttt55ec3HRTn8tCf/sbvPQS7LqrH3m+zjqZ11ZeGd54A8CALfjrweDcrnz9WJ/xBg9eg5EjO7/PVVf5ctdd/RK4DzwAa6wR0ZggraMukgiFzAY0ZYovn3uutvUQkXTafPOu20aN8l+5jjkm/znM4K9/zTzfYQcYPz7/vqutlnn8wAN+ft+WltLr20nuOuq5ctdRv/VWjZAUqZC6yxvQP/7hy0GDalsPEenZFi3yjYGlGDOmtP0OOMCvOvSPf/iJ5M1g8uTuj1tG66iLJEYhswE9+qgv9ce5iMSpT5/OXerFbLhh97kv21e+As884x9vsIH//8wsM/vRQw/BDTf4hsu33oJ77vGvH3IoXLbgZJ5hJwA+YhX+y8b536QW66jrXlHpQdRd3mC6mxhZRKRW+vWr/hyXX+6/QqeemlmdCODeD77IvXwRAIexGh8B8D1+z2/4Ia+yOavxIW+wIXt2/JNeSa6jrntFpYdRS2aDueceX4brFYuI1Aszf8/47rv7hruBA+E3v4Ejj/Svn3RS+efMDphd3o9Mq+A1nMJA5rMj/2I477AXf+dCfsbvF3wzkflBu9wrGgbMUO69omrRlBRQyGww4X/Wv/1tTashIpLXuuvCuHF+icrPP4ezz4Yf/tC/ds45/rVTT81/7L/+FW1d/o/zOdVdzVe/6gNweztMn555/aWX/BK9kajne0XVhS8VUshsIH/4Q+bxZpvVrh4iIuXYaiufY0aM8K2cV13lJ4B/6aXMPptt5keox/kH9Je/PJJVV/WBc9VVYeutYYUV4LXXfIvreed13n/BApg9u8STX3aZD5ClSPJe0SVL4Oij/ZD+e+7x35RzmS78Pfbwr9fiXqw0ht801rkKCpndMLM1zOzPZvahmS0ys2lm9lszG1rrupXr/PMzj3U7j4ik2Tbb+Mna//pXuPBCeOIJv/300/3v6Vdf9bMPrbEGnHkmnHBCZqAQwDh2L3juPfhnt++f3aK52Wbwox/Br36VGYA0YoTPDkOH+gB6+eVw1lmFGyo/evAFOjo6B4xPWIE/8a2uO3d0+JFNcavnLvx6Dr+FpLHO1XLO6avAF7AuMANwwH3AxcC44Pn/gBVKOc8222zj4vbEE090u4//NDt3/PGxV6fHKOW6SmV0beOh61rcP//R4T488GTn+vd3t3Kk+ztfdv9lI9eOuQs4393DQW4pTcv+v3zh+Q533HGZ/z+j+BoyxLk+fTLP//hHX36f37oOcO2Yu4pTOh2zKa92PklTk5s717kTTnBu1qzM9/ePfzi3447OLV4cwcUaP965lpbSvqmWFucmTKjobcr+zHZ0OHfkkc4NGFC8Tv37+/06OiqqV6RqUOdq/i8AnndR5KgoTtJTv4DHgkB5Ws72y4Pt15ZynnoImZMnZz7Dd90Ve3V6DP3Cjo+ubTx0XUuweLH/Rd7S4lxTU+df8k1Nzg0Y4J7c4wL3/tRMUvvsM+e+/e233DPPODdypHP77dd97orja2P+407iWndwr/s7bV999a77/vnPzk2a5Nzyy/u8+Oqrzn34Yf7LMX16nut02GFdr0+hr6Ym5w4/vKIfR9mf2YTCb6RqUOd6CJnqLi/AzNYFvgpMA36f8/LPgPnAN8ys0nUnEnPDDbD++v7xPvvAoYfWtj4iIjXVu7fvSx83zv+HmH1v3KhR0NrKrv/8GasPz9xXNGQIHHPMu+y0k++aHzvWJ4L2dnjkEbjvPj/D0OTJvgd5wgRYb73oq/46X+A6vsO97Qd02v7BB133PeEEvwzorFn+ftXNN/erJzU1wcUXw1pr+a79Pn1glVX8gKobbvDbfvc7+M/Yt/lVxw+ZzwBO4AaO5hbO5ZcYjpv4Bh+zEtMZRhv9eLJjZzoefJj2dn9Nwjp1dPhlSB9+OMKLUK/3rxaTxjpHQPNkFhbesPN351ynG1Gcc5+b2bP4ELoDlHADT4I6Ovx/Eq+84m+Yz3bffTWpkohIfTGD7beHO++s6jRNTbD33l23b799ZiWiGTPg4499MAVYvNiH1jPPhC228P9XZ/sF5/FjfgnAYGYzhyGcw6+4mHOrqmvIOTg3z6l+n9WccvrpAM8DcB6/6rLvceRZymkBFaaKkYAPuhtuCE8+Cb16wc9/7m/3HDoU3ngD+vaFTTeFvcauRHvHibzHmjzNLoziLhbThwfZjz0Yx1q8yzBm8B5rslfHY3zwwCe8eDUMGOC/99ZW2HVXf1/vuHGw/PI+iB9xhA/Is2bB44/7n09zM7z7rv95PfmkP65fP+jfP1P7JUv8LbI77QQrrug/Ws75Fa/mzfPbeOghlnQ08RGrsybv8SS78Qw7813+wAr4KQocYFDwnlvn0reIikJmYRsG5ZsFXp+MD5kbUOOQ+b3vwaOPbkufPv75G2/k3++99zTgR0QkacOG+a/cWT1OOcWHmBdfhJYBjgE/Oo0hj/+NQQs/5nhupDdLWJFPl+3/K87jb32OZvyIo/n+o/swf4HRvz8MHuzH5my6Kdx2m1/VbautYOZMv3zwI4/41lWAvfaCqVPhzUK/2Wpo+vTMgKr2dvjxj7vu88YbcHdO5+LT7Lrs8ZNBYO1kIXBa500339x1t7PP7vz8nHOK13ellfw1zrXaavDZZ7kNl/PznuOnXMTOPM0z+Mmrm1nCUnozbP50ZgSBsl8/P5vBhx/CJpvACy8Ur1c9UcgsbHBQzinwerh9SL4Xzewk4CSAYcOG0draGmXdOpkyZSM++mhFoJ3lllsCdF424/DD3+Nb35rKlCkdTJkSWzV6pHnz5sX6s2tkurbx0HWNT5zX9vPPwb5/EBvNf52W555j2KIZNPkxAAB0mOH69GHXnaaz8rl9mTrtyU7HDx/uW832399/ZTvlFF8uXmz07u2KtoZNnjyQ9defx4IFvdjq1z+j46n/cQ+Hchh/Ywiz6cciOjAeYy9u5Wh25Sk25TU241U+3XUX/nfeeZjB448PY4UVFtHUBB991I8pUwby+efNvPjiUA4++APWX38eP/3ppgCsvfZcZs4cwIIFzQwatIQBA5ay3HJLGTx4CSNGzOeNN5ZjwIB2XnxxCGuvvYDJk5djbx7hUb4GwIrM5FDu5u98lamsA8BavMO7rM2avMt7rFXdDyfLWmvNZ/bsPsyd25sBA+YDne+YW3nlhfTq1cHQoYZzfWhvN5ZffjHzZixlPgP998s03mE4AFvyEovoy8rM4GOG0YfFDGABa9gHzHCr0K9fOwsX9mLaNIdzxsKFc2htfamkutbF/wVR3NjZE7+A6/Ct198u8PovgtfP7e5c9TDwRyqj6xofXdt46LrGJ5Fr29HhB30cdlhmYFJLix9UM3Fi/O+frZzBKgMGJDe6PKEBSZGqQZ018Ke+hS2Vgwu8Hm6fHX9VRESkIWTfKzpvnu83njfPL4G03XbJ1mX77X2zaPYNiPn07w8HHJBc/UaP7r5OoX79/P61lsY6R0Ahs7DwzsYNCrwejNcueM+miIhIepnBTTfBgQdmRuBna2ryo2kOPNDvl9SolHoNv8Wksc4RUMgs7Img/KqZdbpOZrYcsBN+LN34pCsmIiKSiBKme+K225IdVVqv4beYNNY5AgqZBTjn3gL+DgwHTsl5+UL83b5/dc7lHzImIiLSE9RTF36oHsNvd9JY5yppdHlx3wOeA35nZnsCrwNfxM+h+SaQZ4IFERERiV1Ec50mKo11roJaMosIWjO3Bcbgw+Vo/HrmVwI7OOc+LXy0iIiISONSS2Y3nHPvAcfXuh4iIiIiaaKWTBERERGJnEKmiIiIiEROIVNEREREIqeQKSIiIiKRU8gUERERkcgpZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiEROIVNEREREImfOuVrXocczs5nAOzG/zYrAJzG/RyPSdY2Prm08dF3jo2sbD13XeFRzXdd2zq1UbQUUMnsIM3veObdtrevR0+i6xkfXNh66rvHRtY2Hrms86uG6qrtcRERERCKnkCkiIiIikVPI7Dmuq3UFeihd1/jo2sZD1zU+urbx0HWNR82vq+7JFBEREZHIqSVTRERERCKnkCkiIiIikVPITDEzW8PM/mxmH5rZIjObZma/NbOhta5bvQiuiSvwNb3AMTua2cNmNsvM2szsVTM7w8x6FXmf/cys1czmmNk8M5tgZsfF953Fz8xGmdlVZva0mc0NrtnN3RyTyLUzs+PMbGKw/5zg+P0q/V6TVs61NbPhRT7DzsxuL/I+ZV0nM+tlZmcGP7e24Of4sJntGMX3HSczW8HMvm1m95rZlKD+c8zsGTP7lpnl/X2nz2z3yr22+syWzsx+bWb/NLP3sur/kpn9zMxWKHBMej6zzjl9pfALWBeYATjgPuBiYFzw/H/ACrWuYz18AdOA2cAFeb5+kGf/A4GlwDzgBuCS4Ho64G8F3uPU4PVPgN8DVwDvBdsurfU1qOLavRx8D58DrwePby6yfyLXDrg0eP29YP/fA58G206t9XWL+toCw4PXXy7wOR4VxXUCDPhb1v8hlwQ/x3nBz/XAWl+3bq7pyUHdPwRuAX4F/Dn49++AuwjGIegzG++11We2rGu7GBgfXM+LgauAScH39AGwZpo/szW/wPqq+IP5WPDDPi1n++XB9mtrXcd6+MKHzGkl7jsI+BhYBGybtb0f8FxwXY/MOWY4sDD4xzc8a/tQYEpwzJdqfR0qvHa7A+sH/5GPpHgQSuTaATsG26cAQ3PO9WlwvuHVfN91eG2HB6+PKeP8ZV8n4KjgmGeBflnbtwt+rh8Dy9X62hX5nvcA9geacravArwbfG+H6jObyLXVZ7b077tfge2/CL63a9L8mVV3eQqZ2brAV/EB6vc5L/8MmA98w8xaEq5a2o0CVgJud849H250zi0EfhI8/W7OMScAfYGrnXPTso75DPhl8PTkuCocJ+fcE865yS7436UbSV278Pkvgv3CY6bh/y30BY4vob41Vea1rUQl1yn8+fwk+LmFx0wC7sD/fEfFVN+qOefGOefGOuc6crZPB64Nno7Mekmf2RJVcG0r0XCfWVj2ecvnzqBcP2tb6j6zCpnptHtQ/j3PP/rP8X/VDQB2SLpidaqvmX3dzM4zs9PNbPcC967sEZSP5nntKWABsKOZ9S3xmEdy9unJkrp2jXy9VzOz7wSf4++Y2eZF9i3rOplZP3zrxQLg6VKOSZklQbk0a5s+s9HId21D+sxWbv+gfDVrW+o+s82l7CR1Z8OgfLPA65PxLZ0bAP9MpEb1bRXgrznbpprZ8c65J7O2FbyuzrmlZjYV2ARYB38fXXfHfGRm84E1zGyAc25BNd9EnYv92gUt86sD85xzH+Wpw+Sg3KCK76OefSX4WsbMWoHjnHPvZm2r5DqtC/QC3nbO5QsLqb22ZtYMHBs8zf6lqc9slYpc25A+syUysx8AA4HBwLbAzviAeXHWbqn7zKolM50GB+WcAq+H24fEX5W6dyOwJz5otgCbAX/E31vyiJltkbVvJde11GMGF3i9p0ji2jXq534B8H/ANvj7qIYCuwFP4Lso/5lza0ycP4shBV6vZxcDmwIPO+cey9quz2z1Cl1bfWbL9wP87W5n4APmo8BXnXMzs/ZJ3WdWIVN6NOfchcH9RDOccwucc685507GD5Dqjx/pKFK3nHMfO+fOd8696JybHXw9he+tmACsB3y7trWsT2b2fWA0fvTtN2pcnR6l2LXVZ7Z8zrlVnHOGbxA5BN8a+ZKZbV3bmlVHITOdumsdC7fPjr8qqRXerL5r1rZKrmupxxT6q7CnSOLa6XOfJegi/FPwNKnP8ewCr9cdMzsVuBL4L7C7c25Wzi76zFaohGublz6z3QsaRO7FB/IVgJuyXk7dZ1YhM53eCMpC90SEo9EK3bMpEHZBZHfZFLyuwb1HI/A3t79d4jGrBud/v4ffjwkJXDvn3Hz8vHEDg9dzNeLnvsvnuMLr9BbQDqwT/LxKOaZumdkZ+PkGX8OHoHwLL+gzW4ESr20x+syWwDn3Dj7Eb2JmKwabU/eZVchMpyeC8qt5VlpYDtgJf0/M+KQrliLhyPvsf4zjgnLvPPvvih+x/5xzblGJx3wtZ5+eLKlrp+vdWb7PMZR5nYIpUJ7D/5x2KeWYemVmP8JPHv0yPgR9XGBXfWbLVMa1LUaf2dKtFpTtQZm+z2zuxJn6SscXmoy9lGu0MdCSZ/tw/Ag5B5yXtX0Q/q/scia6HUEPnYw95/scSfeTscd+7egBE1tXcG23JmcS7GD7nsH364Adq71OlDax9aBaX69uruVPg+/heWD5bvbVZza+a6vPbGnXdANgcJ7tTWQmY382zZ9ZCw6UlAkmZH8OWBm4Hz9dwRfxc2i+if8H/Gntalh7ZnYB/sb0p4B38Ev4rQvsi/9H+TBwsHNucdYxB+GXSFsI3A7MAg7ATwNxF3C4y/lHY2anAb/D/+O7A79M2ChgDeAy59wP4voe4xRci4OCp6sAe+FbH8I56T7J/t6SunZmdhlwFvB+cN4+wBH4+5dOc85dXdU3noByrm0w5cv6+H/v7wevb05mnrqfOucuyvMeZV0nMzP8BNCj8IM5xgb7HoH/93Koc+7+yr/reAXrMI/Bt/pcRf77oKc558ZkHXMQ+sx2q9xrq89saYJbD34FPANMxX+mhuFH4q8DTAf2dM79N+uYg0jTZ7bWSV5flX8Ba+Kn6Pko+NC8A/yWrL88Gvkr+Id6G/4/n9n4SYNnAo/j53azAsfthA+gnwFtwL+BM4FeRd5rf+BJfJCdj1979rhaX4Mqr98F+L9mC31Nq9W1A74Z7Dc/OO5JYL9aX7M4ri3wLeBB/Apf8/CtGO/if1nsEuV1ws+dfGbwc2sLfo4Pk9PqVI9fJVxTB7TqMxv/tdVntuTruilwNf72g0/w91POCb7/CyjQYpymz6xaMkVEREQkchr4IyIiIiKRU8gUERERkcgpZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiEROIVNEJIXMbIyZOTMbXuu6iIjko5ApIg3HzDYys6vM7DUzm2Nmi83sQzN7yMy+ZWZ9a13HOJnZmmbWHoTUX3azr8v5WmRmM83sRTP7k5l9zcx6dXOOFjM7w8zGmdnHwfWebWYTzewXZrZOtN+hiNQDrfgjIg3FzM4Hfob/I/tfwPP4pe+GASPxawa/4JzbtlZ1LIWZjQGOA0Y456aVeeyFwPn45QBnAGs655YW2Df8JXFhUPYChgCb4Je364O/hsc4597Mc/wO+LWPV8evg/xP4EOgBdgK2BHoAHZwzr1YzvchIvWtudYVEBFJipmdhw9L7wGHOecm5NlnP2B00nVLStDqeAIwF7gZ+B5wAHBPseOccxfkOdcw4CrgMOAfZratc+7jrNc3Ah4DBgLnAJflhlkzGwH8GhhU+XclIvVI3eUi0hCCexcvAJYA++QLmADOuQeBvXOO/aaZ3W1mb5tZm5nNNbNnzezrBd6rNehabjaz88xsctDN/J6Z/drM+uQ55iAzu9nM3jSz+cHXC2b2fTOL8v/qrwFrAHcAfwi2nVjJiZxzM4AjgVZgTeC8nF2uwofHXzvnfp2vtdQ5N9U5dzi+VVlEehCFTBFpFMcDvYG7nXOvFdvRObcoZ9MfgLWBp4DfArcHz/9qZv9X5FS3AqcBTwfnaAN+CPwxz74XA1sDE/Dh7CZ8C+CVwF+K1bdMJwXlmOA6vAB81czWruRkzrkO4KLg6VFmZrCshfLLwELgNyWcJ/eai0jKqbtcRBrFzkH5zwqO3dQ591b2hqA18hHgHDO71jn3QZ7j1gU2cc7NCo75MfAKcKyZneucm56177553qMJuDHY/+pCra+lMrPVgX2AN51zzwWbxwDbAN8GflrhqZ8BlgIrA8OBqWSu9wvOudkVnldEUkwtmSLSKFYNyvfLPTA3/AXbFgO/x/+xvmeBQ38UBszgmPnALfj/ezsNLCrwHh34lkyAvcqtdx4n4AfujMnadiuwGDihu1HihQStkJ8GT1cKyoqvt4j0DAqZIiLdMLO1zOz3ZvY/M1sQTucD3B3ssnqBQ5/Ps+29oBya8x4rmNnFZvaqmc3Leo8XunmPUr+HJuBb+JHcN4XbgxA8FlgN2LeatwhPWcU5RKQHUXe5iDSKj4CNKTOsBXM4TsSHwqeBvwNzgHZ81/BxQN55NQt0E4eDX5a1GprZEGASMCJ4r5uAWcG+Q4DTC71HGfbC30f6WJ6u/THAofj7NR8o98Rm1g9YPng6Myg/CsqqwrGIpJdCpog0imeAPfBd2zeUcdxZwArA8c65MdkvmNlR+JBZrW/jA+aFuVMFmdmX8CGzWuGAn72y5r7MtbeZremce6/A64XsjP99MiNrzs5ngnJbMxvsnJtT5jlFJOXUXS4ijeJG/PRFh5rZF4rtmLPiz3pBeXeeXXeLqG6xvoeZrQLsh58b84YCX8/iW1dPKPPcTcCPg6e3htudc1OBfwD9gLNLOE+PXmVJpBEpZIpIQwha2C7Ar1DzkJnlXdHHzPbGjxoPTQvKkTn77YVvgYxCoffYCjg3gvOfgG9pvMU59+18X8A38fdTfqvUeTnNbGX8dE4jgXeB3CUqv48Ptuea2Wgz69J7FtzvejvwpQq/NxGpU+ouF5GG4Zz7ZRB0fgZMMrPn6Lys5K7A+nQesHMNfo7Nv5nZXfglETfFT9h+J3BEBFW7Cd/a91sz2x2YHNRjP/xKPBW/RzBvZRiG/1RoP+fcFDN7Eh8YvwY8lHOeC4KHTWSWldwZH9on4peV/CTnnK8HYfwu4FLgdDPLXlZyC/zSlA6/6o+I9CAKmSLSUJxzPzezv+GXU9wdHyD74afgeRkfdm7O2v/VIPhdhB993Yyf6/IQYDYRhEzn3Idmtgt+Qvad8YN0/hfU8R9VvseX8fd7vlTC2uDX40PmSeSETHwwBz/d0efAO/hwfDfw92C6pS6cc+OD5SVPxC9fuS9+ENUCYApwGXBd0L0uIj2IOafZJkREREQkWronU0REREQip5ApIiIiIpFTyBQRERGRyClkioiIiEjkFDJFREREJHIKmSIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5ApIiIiIpH7f3WO0mR4X75JAAAAAElFTkSuQmCC\n", + "text/plain": [ + "<Figure size 720x576 with 1 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "from scipy.signal import find_peaks\n", + "\n", + "x=ds[\"Counts - Charge distribution\"]\n", + "\n", + "peaks, _ = find_peaks(x, distance=200)\n", + "\n", + "plt.figure(figsize=(10,8))\n", + "plt.plot(peaks, x[peaks], \"or\", markersize=12) \n", + "plt.plot(x, \"-b\")\n", + "plt.xlabel(\"Canal ADC\")\n", + "plt.ylabel(\"Cuentas\")\n", + "plt.grid()" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "id": "94385c3f", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([ 131, 361, 599, 821, 1056, 1278, 1520, 1750, 1950, 2172, 2387,\n", + " 2605, 2862])" + ] + }, + "execution_count": 11, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "#Posición de los picos\n", + "peaks" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "id": "b1915eef", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "distancia entre el pico 1 y 2 = 238\n", + "distancia entre el pico 2 y 3 = 222\n", + "distancia entre el pico 3 y 4 = 235\n", + "distancia entre el pico 4 y 5 = 222\n", + "distancia entre el pico 5 y 6 = 242\n", + "distancia entre el pico 6 y 7 = 230\n", + "distancia entre el pico 7 y 8 = 200\n", + "distancia entre el pico 8 y 9 = 222\n", + "distancia entre el pico 9 y 10 = 215\n", + "distancia entre el pico 10 y 11 = 218\n", + "distancia entre el pico 11 y 12 = 257\n" + ] + } + ], + "source": [ + "dist=[]\n", + "for i in range(1, len(peaks)-1):\n", + " d=peaks[i+1]-peaks[i]\n", + " dist.append(d)\n", + " print(\"distancia entre el pico\", i, \"y\", i+1, \"=\", d )" + ] + }, + { + "cell_type": "code", + "execution_count": 13, + "id": "61992f52", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "promedio de la distancia entre picos: 227.36363636363637\n" + ] + } + ], + "source": [ + "print(\"promedio de la distancia entre picos: \", sum(dist)/(len(peaks)-2))" + ] + } + ], + "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 +}