{ "cells": [ { "cell_type": "markdown", "id": "970a0408", "metadata": {}, "source": [ "# Correlating neutrino data from different experiments" ] }, { "cell_type": "code", "execution_count": 1, "id": "a1666102", "metadata": {}, "outputs": [], "source": [ "import pyvo as vo\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "#retrieve the full Antares test dataset\n", "service_km3net = vo.dal.TAPService(\"http://vo.km3net.de/__system__/tap/run/tap\")\n", "dataset_antares = service_km3net.search(\"SELECT * FROM ant20_01.main\")\n", "# retrieve the full icecube dataset\n", "service_icecube = vo.dal.TAPService(\"https://dc.zah.uni-heidelberg.de/__system__/tap/run/tap\")\n", "dataset_icecube = service_icecube.search(\"SELECT * FROM icecube.nucand\")" ] }, { "cell_type": "code", "execution_count": 2, "id": "acd07570", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEGCAYAAACO8lkDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACyiElEQVR4nO39fZhkV3UfjK7d1d2a6kGiiuEjQElqEMYEGxmDQpsEY3GTYANJHBs7QGg99usMsjF5lTgXGRzfEJtrcwN+ndf2GwenEJCJQODI4ITY8VdsNMSyJWskkBCyRzKSGM0wIKlnRkyPuqe7qtb945xV9Tu/s/Y+p7p7PiTOep7zVNWpffZee++11/feJ6iqNNBAAw000ADCzLlGoIEGGmiggfMPGuHQQAMNNNBACRrh0EADDTTQQAka4dBAAw000EAJGuHQQAMNNNBACWbPNQI7AU9/+tN1cXHxXKPRQAMNNPCEgttvv/1RVX2G99+TQjgsLi7KgQMHzjUaDTTQQANPKAghfCX2X+NWaqCBBhpooASNcGiggQYaaKAEjXBooIEGGmigBE+KmEMDDTTwzQubm5ty+PBhWV9fP9eonLewa9cu6fV6Mjc3V/uZRjg00EADT2g4fPiwXHjhhbK4uCghhHONznkHqiorKyty+PBhed7znlf7ucat1EADDTyhYX19Xfbs2dMIhgiEEGTPnj1TW1aNcGiggQae8NAIhjRsZXwa4dBAAw00cJbgERG5K/8836ERDluAD4nIxflnA2cfmvFv4HyE3/7t35YQgvzVX/1VtMxREdkQkV/6lV+Rxx9//KzhthVohMMW4L0icjj/PBvQMMMinO3xb6CBOvCJT3xCXvWqV8knP/nJaJlni8i8iHxyC8JhOBxuD8EpoREOW4D3iEgv/zwbMC0zfLILk7Mx/md6DJ/sc/TNBqurq3LzzTfLhz/84bFwuOmmm+TKK6+UH/qhH5IXvehF8ta3vlWerio3/dqvyde++lV5zWteI695zWtEROTtb3+7XHHFFfJt3/Zt8m//7b8d17u4uCjvfe975VWvepXceOON8od/+Ifyyle+Ul72spfJD//wD8vq6qqIiLz73e+WF7/4xXL55ZfLO9/5zp3plKo+4a+Xv/zleiahr6q9/PNcwLTt9zQbmN4Zw0h1WVVb+SfCuR6rnYIzOYZ9zcauqv4ny1iqntm+3HPPPVOVf1hV78w/dwquv/56/bEf+zFVVX3lK1+pt99+u372s5/Viy66SB966CEdDof6Xd/1Xfq///f/VlXVSy+9VB955JHx8ysrK6qqOhgM9G9/z/fojXfeqQ/n5d7//verquojjzyi3/3d362rq6uqqvrv/t2/05//+Z/XlZUVfeELX6ij0UhVVY8fP+7i6I2TiBzQCF9tLAep1uLOtRvjbSLyUP5ZB6bRrLeqwX5CRIb5J8I0Y7Ud7XmaZ7fSTmoMt6v1/7RkYxci9Ru8S7KxfFeizBPFAjnXawjhqIhc9KEPyVMvvljkQzszcp/4xCfkzW9+s4iIvPnNb5ZPfOITckJEXvyKV8gFvZ7MzMzIS1/6UrnrwQflLhEZ0fP/9b/+V3nZy14m3/md3yl/+aUvyb333CNH8//e9KY3iYjIH9xyi3zxnntk6e/8HXnpS18q+/btk6985Sty0UUXya5du2Tv3r3y6U9/WhYWFnakT+dc69+Ja7uWQ0xLNG1nWc+dBrdTOMQ0t61qyDthOWxHO5/mWSvbqonXTrbtQTd/vpso01fVkJfrbBOXc2WBIO12NevHNDjUpf2tWA6nez1VEdVefOTqWhiPPvqo7tq1Sy+55BK99NJLtdfr6cUXX6wf+pM/0Ve94Q16Z17uHe94h773ox/V21T1OWA53H///XrZZZfpsWPHVFX1TT/yI/rej350bDlYuV/9zGf0tW9+87g+hPX1df3d3/1dveqqq/Q1r3mNi+e0lsM5Z+w7cW1XOOw049wJMJw6OmFsjMtWGDEzyFQdZ5qpePXXbXMa3Dw3znb6tqwZ417YIp512rb5Ei0L4Fj9sXqNhlJCZqcAcUCa28o6qvv8tMIhQ7SfCYZ+fBbuVNXb8s8U/MZv/IZeffXVhXuvfvWr9V3vfa9+9xveMBYu73jHO/TXPvpRvVNV/+a3f7vef//9qqr6hS98QS+//HIdDof6ta99TZ/5zGfqRz/6UVUtCod7Hn5Y/8bFF+st992nD6vqLadO6Z8fPKgnT57Ur3/966qauae6XV/taITDDsJWGAg/s1Xmawujq3HtaRrhhQyyk2jXwyFW/5kQHmdKIDOuO2G11Hl+q+3UjUvUKV/HUtkpsP52dGItbNXqPVOWQ12oazl8z/d8j/7e7/1e4d6v/uqv6ote9CJ9wxveML73jne8Y8z0f+3Xfk2/9Vu/Va+88kpVVf2RH/kRfdGLXqSvf/3r9Qd+4Adc4aCq+qk//mP9tiuu0Be85CX6gpe8RH/1v/93/epXv6p/62/9LX3JS16i3/7t367/+T//ZxfPRjjsAGxHKNhCtAWaYg6p//o6WWAxPKbFM4ZjVflY/Yx/XS15qwJzJ4VRzC3G7XnukDpzg2Wn6U8dS8ADzzKcxiW5U2M7LY3tBGxFOCDj/7JmFsKXdxatMwJmzdyh0wfVG+GwA9DRifYTA2YcqDFVWQ6pReuZ5Wdai96pejx8696rAztpydTBYRoLYStjmhqbOnTE7aPAmtZSnHZO6uBzJl2SCB7Tq9L60WV0G1x1geuva2VsF7bTTiMcdgDqmOHIOExA1DWj61gTptVOaxl0NPOFG052/2ws1Lqa8FbcDDEtHmFaN5uHA+O7oFl8ocrCMEZsbsA6fUuN1zQWqDplkI6wrjr9rgNnSnGpA4yrx/Sq4gXbtRy4/jrxiboCJSUAGuFwloUDE1vVb9WMWRgCJkRwQXoLp4o5GnO3LJWtatZ2xRjLVl0XKYgxuo5OtNmtMBR7rqOTsY3hiq6irWq2jKO120ngiOPe1nQfkQZwbKrwS82ZR1ex8jvF1M+kwlFVN/dhK5bDduFhzVw7d+Tf67SHVsrDGhcoKUFTN0juQSMctgBVCyZl/ns+3qX8/pL6bhd7JuaW6ejWFl5dywH7w8yvDhPy8MI6maEb04wxQxQiLDit3q5WZ2xNw6Rj1iEyWitj5VJad8qF4+Fo5bejBFid02yos34hfdSFOtbbtPV5NOGttxRdGtM7EwLBEwJ3wDUNo35YJ8LB8GR8uT2vjie95SAi3yoiX4DrGyLyL0Xk50TkCNx/fVVdZyKVtUq75nuotcYsCF7IvMhjwU/P2qhj7cQWMuLKTBJxX0j85zFIZoBdnVhBKWsImbo3ZtYPE7rLzrPWttUTtDxnCG0omxLUdvG8pvrP93G/AlsO8/l/CwlcGWLCJvU8jtNWrQeep+0KCKzPEwRI8ymcjeltR6tW9Rkvavt3RH6nGDXWaYz/9vzyBAD3YScF3hNGOBSQEGmJyNdE5NJcOLxzmud3Qjh0tKjBpIjRI15jaq3I//xsXXPfmDczVy4f+21MCdtH/zj3G//nTVhVWqrX79RY2DNdzZi1WQ6xufAEB5dFwZcCtAo8holMHYXPVtOBJe8j99drI1Uvj2dsXGPPxZSdOmC4pwR91XxzWZ4/BM8a9fptTG+7WUeecEFt35i5p9nHmDjWycFvT5BxPeyKYphGeDxRhcNrReTm/PtZFQ644Jk5xIjcY1bIlOoscI8pe890dMJYehp3WXmWgy3irlOf1ekt8jrCLTUuQeMLXtW3shAHj1F7bXrP1mV8fS0Gm2PWYWyeqlI1rdySTsbbgtVIvIbDAjwTm5OOTuZ0qxvLdgJSY+wpMymrKiWwjE6Q1r21N63lEBMizPhR04+5eQzau3ePmfifHzyor3vd6/Syyy7Tb3nRi/Tv//AP6x997Wv65Ry3e/Jyt+d1/tIv/ZKKSGE/w/ve9z5dvOwyvfSFL9Rf+/3fj7qiprGWnqjC4SMi8s91IhwelOydGB8RkW7kmatF5ICIHLjkkktqDI0PVQwNF7Npo1XaUWqB97TIoKsWNy6aZXiuym3Cz9pv1FQxs8a0wk5FnQbWhw61V8eX7llZMesq9btKcHv9wGerGGxKs2ULhX8zE8M4E1tKhk9HVefUz45COvQy2aYRiluxGmL1MB6szFh/TRDGhHtKSUAlx2t32phDKnUVme0d6scIvux8mnD407U1vfgFL9Bf/cxnxnh86E/+RD/xxS+6mU2HDh3S1772tXrJJZeMhcOXvvQlvfzyy/W29XX9b/ffr73nP1/vGAwKQew7CIcnpeUg2fHmj4rIs/Lfz8rdTDMi8osi8pGqOrZrOcS0F2MOuCinqTOm6XrunDqAjLeOeycmtDxBiMynTp2pgG6sb6hNe9YPtxPT0mPC155PjU2MaXuAc8/9ROHY1yJB4twb069ynbScOpCOPIHizV+VslE1doar4b2gfno2W0ABytlzhu8slV3Q8vynlASzcFERQZh2ExxbDihU8DsKB2TM3tXevVvvUNWf+/CH9Q1XXRWNGzysmcVg7b/xjW/UL3zhC+Od0A+r6jXve5/+7PveN37uyte+Vv/sz/6s8OztmnY31R2nlHA4H05lfZ2I3KGqXxcRUdWvq+pQVUeSHTb5ijPZ+NtE5EoRebuI/AuZnBz5ofzeCSh7RaSO2MmYvyPFkyjfIyJdEVmX7ETO9+f37dmqEzax/ZFMTg3l5/AETPzPThp9q4jsproNt44UTwq15+2E0B/Py71QMgn+BqrnbSJyXESOyeQU2Q/lz/x4XscRERnkn4clG2fr/9sBd+vHet7WlXl9V+a/n5vjdhX08W0i8sG8vVUpjuWHRORU3se35GPxQSmfdmt9Pg331un/x/LvQconjR4WkZskO0n3ehF5Sj4m74XneU6YtlZlMuafyJ/fI9m43pTf/2mnnqqTeFPl3isZvR8XkY/n3x/PPz8hZVpuwbOal7sp7+/jIrKW12WvqJnNPx/P7/2FTE4bxpOHbXwkv/erOc4fqOgbwyMicodk7oW/lMwVcb9kY3upiDw/L2dvZzsqIs8Qkcvzz55kmuul+W97Uc/T4NO0WMnLH737bnnRy18uIiJzeZsnAaejMjmR9Xc+8xl57nOfK9/xHd9R+P/okSOycPHFY1ye0+vJnx05Iofh2UB1nhGISY2zdYnIJ0Xk/4Dfz4bvPyUin6yqY7sBaQzAsm8TL3TDILA21kuUZ/cAasVcD2ra6A4SLbod2GzHrB6vTg40p8CeZ+1vGsuF+4yWmZe9ZWXMX49tmxZt8xXDBTV77gtbCzGXEAZf0XLA/izpRFPGOUIXJFsOnvZu+AYt4u5ZoF78IeVewv5xXxHQAprVsuXQ1aIVgRYgWzb428otaJGGRNNWCa5Fz6Kw526NWA6o+fP1r1X1OXldKXcUauv3OOWsjfbu3aqquvxTP6X/6ld+xW0TLZKbT53Sl73iFXrixAlV1YLl8Kaf/En99euvH7f/j37sx/T9v/VbY4thmr0VCE8ot5KILIjIiog8Fe5dLyJflEzofgaFRezarnDwFgyb9B2dLFpktDH/ccyvywu7A2XZpcKuBmMEMT+vxyyr6ky5tvpaZnrIFGN9tMWNzI0ZJI+ftTELZThryvrW0SL+3nh7bi/sP447u5gQV55jHkdOShAt0ojH6Dw6Mca54OCgWpxH7B+3x0IRGX5Ly4oQM30ba4wZ8JgyPXBfLcDeAbwNJ09ApE7M9ZQbHvs/hJgDBpDR9WKXBYOfCfWmmKwnYO50/jfh8P+77jp9fe5WOqBFofLmH/1R/Y7v+A593etepzfddZd2n/EMvfjSS/XSSy/VVqulF198sR49elR/9n3v02vArfRdr32tfvjP/mxbabqqTzDhsFPXThyfEfPTp8qgBom+UbtMC0xpPTG/vC3ylmZBSquTFwz7iVmj8wRD0HIKZVX/vL61tcw0jZGiUPGsC64fL+vTvE6CsyjA2eLwICbwDS8b066qLubfF7XIyDA2FBsrnEe2KFgQcGAVmZyNrzHlqjhLh8qz0OA2JMeBrbGYAsLtc5vWlwBtp5QPpEubT6MhUx5wvq0+HBejL6vb5utz99xTCiBb0JYFhGnd/0Yzy+GX4b/bNW4ZGLOPpbAu5MLh1scf195ll+n//Tu/M372P/7e7+ldd91VKM8b6NBy+OTdd+u3XH65/vn6uv6P++/X5zzveXrLYLDtwwEb4bBFiLkcPPCsCnTX8GVlWzpZlLaoYkwPXVLMtPEZj5mhYLG6vKwZW+TmFmBm3tdi//pa1Nr5woVbJ+C75NRhTBXr9IRolXvEaz8mjNg9g8wZr7YWGaHn7mDmiHXaPZwj7icLiVh7sfaNDrBNbwysPXZZWrKCp7zYmC/SmDBtcn32H68JpDUcn5Zzz8p2nHJ/eM89JWFgaaIcRL5DixALMnNGEjNztjYWdu8el/9vf/mX+prv/V593gteoM//m39T3/CmN+kffe1reodOrIkDUP8dqvrsSy/Vv3zkkXE7b/+FX9DnPv/5eskLX6i/9j//Z3RfxDTQCIctgqeNMuNOmbo9nSwa2/WKGRwe4j0tamDeYvaYgOeqmId6bcHaAscFhZoZugGYmcUsqZQQRHeCaYNd+o5MBJ/FdkxotLTIuK0ejEG0tThHzHSREaIw9+pgS4sFhMdkPUvCc7uwoLW6ea7QvcS42LM8Hz2qs+v855XHftgY2fh5Fpk3tmhRdgl/Xk88J9ifvk6suSUtj78JGVRqzAL5vXvuGQsDY664Ye12Le5KxvRPK4OMOxYnQKsB2/F+I8QEkPefhwNaP1X7LVLQCIdtgDElI05kCkaoC1rWpj3/swFrk3t0IkD6OmFS5mtPCQKsk8shQ0I/u2oxp9/wxHpsgaOWx8zLc011tax5osDyfNTWD3Z3YJ88V0fsQuERswwwFdJcS+K077lKYgKCLRecJ2tvVuPMGQUVCxR023CMhS0x+6+jvqspJuQZn5gQ8oQP4ulZWPYfj9GcU8YTVDgWC1rG39bpYv75e7nloFpk0sh4UXB4bqaHtcyQU2copYQDCxP8bTEPTKNNBc5jsZOtCIhGOGwRWBggozGNFRv1/Mam/aFJztYAllctLm7Dg4UP1+kBauLsarE2ZnXi5425I6wdzk7qaXHxohWwRPc79BuFj9XP5yQhczWmYpo8M3O7FrWsRQbNBK/hb3ERe4YZlLUZEwIpS8nowMbf+oNCkRk00kMXvqcIfIFwQIsU4xiGA2rvWI+NdYymPGsIXZNMI9Y3D+c2PNvRopWM/RCoq+OUsXVh4LX3e/fc4+5ZsP0MvKmNtXMvXmHWhTFou9DdhPdQGKBQqvtSHhMAB7QslKwOrHcr7qUn4j6Hcw4fkmy7tcG8ZLn0PRHpS5Y+9YH8t+VqH5csz3woWa7zB0Xk1ZLlUH9cJjnhV+b/v0WyHO63wG+Bej8gkzz/45LljL9dRD4sWR7zb8okt92DV8sk7/wm+u/x/HMgIrskyz/H3Py3SZY3bfnth2WSm76Q43dl3rcF6P+P579vze9pfj/kl+b335jj9laZ5OjbXofr8zKYf/9qyXLK/5lkee6/ItneBGvb4KG8L8/Nn/lE3uZQRN6ct7mW99lgE74vymT/hML9kNfB+xc8OClZTv7v5s/cJNn+lVZe59tlsl/hvTLZHzCUyf6N2P4Zw2Utr6sL9Z6A/opk+04+nte7mT/3ONX1cSnuWfmYFPfDrOZtXJmXQ/q2NneLyH35cx8Wkc8Rrssy2cszzO89lpfHMTZ6sDm3PRY4VwZXSbanoCs+/Y8kU3Lvyn/bPoXHoMyzZbJvAZleK/8PcRPJxuLZUtzLMZRsn8lG/n8L7tl6sf0Q1obKZA+FwSOSpWI+At87gJuN+VPz8jYmG9SfaSCTA9NBIxykzATWJFt4z5VscV8F/83B98ckI1jbSGUErpJN7nskYxbGND6Uf35QJkwRN/+8VyYLSvLvt8pkwYtkC97bKGfPiog8nONlZZCo1iQjxFNUh5HOLplslOtJxpjfI5NNUWtQVvPfBnvyZ14IZY6LyKdhDGyz3aqUN6/ZOPy0FAWhbax7GrQVZLKR7T2SMTzr/1Cy+bNF9n7JmFag5x+Dsv9AJpsAre9XSqY0sOCYhd+DHNeTgIttxDPmoc5/PcnGeiiZoOxAnUuAS0cm9PQGyRiI5Djg4r0V8NwknBdlIqyNSRvTM+FoGz53y0SAP5a3/0GZKDHvkYkycKtkdGFgwl+g/ZaUGa+IyEWSCXSb8ytlIggN2jIRgqYw4cZEg7/etUseWVmRDVV3Q1hLMqZtm8qM6T5NRL5TsrmzzWUzMpk3cb7b/0+B+wYqIp+XbOyM1oJkc/ZsmQiCwzIRGLYB70T+uSsv/xTJhI5ItlY3oM5pmbaqysrKiuza5YneOIStSJTzDa644go9cODAlp//kGS7o9cl2459QIoT3xKRCyRjzPMy0UhFsgXzUP79Ksm0MZFsgR+RjOh/VzLCMe3cLA3bDfpeKLcuGYEcr8AZ2+W2DbqSEdhumWiRbajf/rcxMM2WdwxfLBlBi0w0rVnJDsDycJqV4vgFyQSt1W31GcPqSLG/iG8f8LlKMmY0J5kV8QH4ryvZAmtLJqSulIkwEskEzWNSZFSzkjF369eFeR3WD+y39e9KyXa+b0hROAaZMJgP5e09LhmjfoVktGD9xzm/SSZzf0ImguQhqMvKmsCrA0EyJvitktFzgL7ukWxz0ayI/EeZ7EQXyej2oExOBrDxRLqYy+uazb+vQdk1KY5rCqyfHu1afevQzrxM1pDBsoj8/uam/JvDh+UF6+tyYf7MUyUTJKcko6enSyYEHpPc0pCMyXZksgYkL3sq/45MuE33UZio+AJwRiZjIhXlRvD7aYSTV9Zwf0yyvl4YKW+wa9cu6fV6Mjc3V7gfQrhdVX3jNeZveiJdO/2aUPPJol/cyyoSzfzbHS36k9Gfiumr7OdWLfqg63aYff6q8UwS64fFGrAv7M+NjQVnjSjhy4FFi0Hs0WJcAevEWAL6yzm+g4F4L2DJfnQvJuONjXdheqfFPOy/BcIhhgeXwbnlOfdiNEZLsYC3d1n981qOI1TRVYf6yRsPMS5ggOsjFY+puiymxHXManHvB7aPuHpzj2OJdGtzyusjRMpg7Mbowu5zAgfPDY9d1cXPd5y67DsmO0yzjmMgTUB6+4BB0o7G0xXrIIzBayPIRS0Tb0eLgUMLdqLg8bJq7K1rHHzFoLItLsveCfAb28QsrRY9Y20ZQ+9qOegaC6DjmNonPmOLYAHqYKFi5WNg5THNF7OEeE+B9VmovM3tMv03p77w43Le8eo4110tMmWmL8wkYrrDOQ5aFpCYFhyjR04X9miOj8ZAHIODl3cxM+ZrUYvzyxltSuPEmVWqRdroVOCDa8Y24FndKGjaWkyzxftd+B/XUJ0kA+474tXTCX0inXkZgJydOA00wqECqtL9sJxNOmo1szpZOB7x8726Ghcza485eJklEinb0aLm1dGyZuJp50b03I613aE6WePj8e07ZWxx7NFy9hCnMyLzXXbqN0C87JrL2zaBYXPnaYC88Hg8Y+/D8CwHFnAsUFmgsLWJzArb5hRlxqcD5ZGOuJ+mUKCV2Vd/oyLTNlt7xiiF7rM1xpfhzbvx7dll9Xd2Y3/5DX+Y9sregFgmWgBcvXXEGXAeDXA/F7SYWadapBPcLGmCGS0X9maYkGIPwrTQCIcEoLbfdv5H7TZm2k/jEkqVD1okqpRGFjvuGhclMxxsG81m7L/1l60OXtgxq4WFCrulYto549lx+m7pkbYwTcPqwjM4Jh4TjF32HC/6DowvM2xM7/RSg9FiYRcE44pC2dOweb5MYDITQlz6VN7Ao+NpaRivRcAbLUrTqmfye3MaT0vuAN6xMkxbOK6qRToyYdjRyfzi/NjYp7R8diGnxqurRcHKtOXhgfia0IjRM+Ph8aqtQCMcEuD50w2YaZl2xQTkMTdbGG0tbnQzC6OjPvOLEaNXv6ctsxbNG7LQXcDts2mNuPOFCzLoZHOgVw4XC7uNWAtb1OKxHt7CDYSXLSJrH60pxskTejHtnV1GvBsYF23QspDwLAhv06Bqcd5iNMBMxyuLzKdHzxrT5XlHxhTT7OtYurhhr6N+P7x6OjRuuMZQc0f3F5/rpZG6zU3nzUtMIcA9MixITCHgeFRqXOzZDuGbmmemPW+udwIa4ZAAJER0neB9Xnwek0DGMpvXY5PpmfvWDjI/W6S4IGa1TLwC92wRxkxLjBkgcdYRQr1IOXTlsL8aj7ww/PBZPJDOY/z8PJbzxsHKoGtqUX1GO6f+xjbTPhEfYxA4vihozN3Ez2HfzPT3Nt7ZxRYHCnAbO8QVrZGe80yf7pmSkmLuyPi8snZiQIf+M4vALFw8urujZXcTuqDmqY0OjEFHJ/RqzzNd8Lrqa7x/7JNHRs3KA1rPXqyIxwjngBM/eBNgh57xhEpLfSvD2ud+8WbSaaERDglAYoxpyUb0HS0uPtRUceJMkxadMPcuPNum/5B58HEJosVgFRIGLzxPQBiuHXoWNSLPXMcdwh2dMIJFLTLEBS0Hy/Bi10qX2vEuz91lOHuMKybosJ9sVTCuKFwMP6MLrMtzZ9gYGZ6IYyxbzepAMJyMmSGT4udShwD2oJ5OxVjHLrYCrH/M1BkvrsOEQZVVZf+zIoVja5atrcdZ+I2Bah7vjjNO1j4rHAtaVlrMJedZgrhGTKDxvOC8ekoLWm6maGJcpMqiNFy2Ao1wqIDUoDMjNUATM2VpMIHWNbWrLu9cG89FgfjGGBsGylgTZdcOHi3C2Ug2BpaK6LmuOlq2vgynRfVdL9g3vI8ZNDHrCvuPmrUtZtNyPd8ztjUL9eCREjZWHQdn72KrCPGz/8y9iK63GD2pFt0kVhcLWDyuo+pC12fXqWsrl807WimmcNg5YzG68lJacR1hHA01dnYp4ThZOeybWUAe/mgZYBID44I4GnjrkV2YRg94n92xtiYxlXYx/76kW4NGOFSAt5i8N3wh0bJ2Ns2i4/tVZj/XgcSMDM6Ihy0aXJSsmXF2h2rc74wZHMiYbexmtCxkOLvCMlmWoN+eoLDvqI3FXAc9LTN2ezYlJK1ewxPn3hMUWIfnC49drMnaGCHgvKaYsWnu6CrxmA/jvwT1zutk/vY4baAVyPRt/v6O+sd9B/VTqNEC60F5xM+zztlatHEzxWNBi0wb1+oi1M3jhJaD0WCV8PTcw0YLHFf0rHgEbww81xnzCVxTPE9bgUY41ABvkGNavi0M3qiDBCmaLZKOlgkA6+vohNl6jICFiZVngkGNyQgVXS4m2Nhk5UWsdM8WOpvWqJ0x3j0tamkGXiaQfXr1MMHHLLwlZ5xQqPcidcxH6vQEUVuLTMtwj7mNFvP2jRFz/4yOWBu08ehp8eh3/o+1UPY7e0pIjAatf/gb+8TKgrXdofuWhtymerFsB8YwFkNCJm+0zbh24P+Yu4dp0rOWTWh1tBgLsbgBWtGLMB6sIHi0a2uVYxY27/YdmX0sUcObP6aZrUAjHGqAEZsRXle3tgMU88+NKDz3il1s1rbpt2rZ9aNa9t236NMubBu1No5t4GJAXBe1GJzjBcP7PTpa1vysD95YWZlO5H/W+mNph2jB8DjgYuZNYx5eQf35RzeEtzEpNbdcv4cn99kbE+8tfzj36Bb08LJ5juFTB8eYMuPVge5JY8ZIE9bObud59unjfzj+PIe2RrzU4x7VZXPtjZVBX8uxM29/CbaD4+e5wmIKlOf+xPUVG+sqSyUGjXBIAC6mmHujy2VGoyRC5m5B7SG1gDhjiVPpVH2/MpqzVt7TstCNgv9jv9Ad5REhulA8rTX1TmnE3xa2fTecbKy4bRQgqdgO9pHHpaVl7VK06KKwPvBccHvsSkKfvBes5XvGHFGD5GwnA4/pjZlOToO9kycLloPhm/KdqxYFl0f36PbzAuKp3zbmKaXI6CS1v8LaTpWxtE8cZxTagcrGxsQL+HoKENcz5/xuUVlWEu0Zs5CW1bf+NH8W3x7p4Y/vXZ8WzlvhINnZbV8UkS8Ykvm5U38k2cnAfyQiVa7AbQkHIwhbVDzoOGGopYXBQOfW1wsLVaCuFv1mIpuBT4sboCnNwoDdQawBYVveIKFWaYTWde51tEiE5gsuLKYjR7S/f3/JvYJ9ZY2f/cnMBETjGpz9Z+U8N80MjI+3yG1honWDgEoCPmsxD1vIi1qMvXh4xggVz+HCNnlXMsIy/TcWeqORtk+e1O7x467gNLwLwmk00uUbb1Tt9wtzEsMX6XGahW20lqrbaI+VAXQpVV1BVUO+9sJw6FpqVYIMaYbjG55Q2q1FWvTcxIW1ks8jb25lGsPf6MKMjcVMpJ1p4XwXDk+nex8QkXfn398tIu+vqmc7wsGYbiGINhxqazDQ+cGgMKFjH/BopP1rr1XtZ0t5+cYbNQwGOru+Xkh/NGIxjc6eX9AyEzGTuw2ftoA4kGzPevnzJULb2NDOsWPaXlvT7vHjunT0aCmIyP55Y0jGzAr7Oo4f16Wbbx4LRLQ+Olpk9HZ/jDsIFV4M9lnHXYFCzdPMDWcMZnfhWQZchCygYhsTPV93Ks0Tr07eLgpl+8/THGNCMQyHhbHGOcU+m3Ytqtra3Mxo1+uXYxEb8+R0a2Za81rOuce59KyPlCCdg7baqrqQ0/Hs5maBFpevv15bm5u6fMMN5fqc/ti6MlrrAB7sUq2iR7bUbdyNBxjtqfMcj0WJxvP16rVr1mfMmp8GnmjC4aCIPDv//mwROVhVz05YDh5R2eLzLpyM3smTWfnRqMBMY1qQaSkLGxsqw2GBiD1fJWcUpXyPuFg7qtq/9lrtHTo0+S/RJzNxeVyQcfb37y/gu6gTzbbj1DlmSoOB7vna18bPtjc2tDUY6PLBg25aYeryfNZ8dbRofeEYoqbIcYyO1WvCb22tpAXbOKEGuZjAg+lqYWOjoLlzv9HySmnt2Hdsn1NhO2trunzjjQVXlI0f4zd/+nRprNG6Npy9jDYM7vbgNwsvHBtkcOa2Rddqob9orQ+Hurxvn3ZWVrS7sqLLN9xQHOfhUJduu61kCXd0AiicTWjE5jF2LQOuc5qtj96RI4X1q6q6eOKEymikiydOlKxTc8li27Pr6xpIOcU+dOn3k1E4PCAid4jI7SJydX7vBJU5Hnn2asmOqj9wySWXbHFo/IDk7vV1bQ0G2YTCfdYMzdfHDMZzz+AzYyI0JgvM1ogUJxwXkFJ9RtieCdxTVe33C4wB24ppuuY+sUWKVkxK2+w647CgOXO64YZC2f6116qKaP/aawsuqYJP3dH80L/qaZ5z0DYKAa8sC41ChsqNN2rv0KGxlt2ncqYQtIZDd1xQ0JbGajgspSQzbq3hUPvXXqvLBw9OvbCYKaOCgP74giUyGunSbbeNhaKo6vLBg6q93tjaq3MIH1q2RpfjTWubmxqGQ108caLgIsXYRioLjscQmWeLlJ4wHGp3ba2UTViV4FDHpcVWVIGmjhzR/t692acC7Rl+jnLG+zMK/XTWgBcffTK6lZ6Tfz5TRO6U7OVQJ6jM8ap6tputxJLcJsTL/OnQPXTDcGDR7lf5MefMVM59wv39+ws7bjlDBZ/txMzP0WhcHjfNWB/Q38yuJW+xGJMtBOTh0/yky1DfWPgMh9pZWRnXuaia+b2vvbawqC0I6s2FtxB4LkSLAsrKedZWR8uWRyHWc/KkLt94o7aGQ5cpzq6vTzRxLboTUJhqPiatwUBbGxuF9qsWSGtzUzvHjpXbrvGsjXNPVZeOHi1YqGypmFuCaWvh1Cld3rdvzOg82g1aPK563F91AqlEN+2NDdVer8zYNZ0FVriGw7Fg6WgmgDxLHJUFq7ujmQBEDwG77ti3j25HTj9dyOvrHDumCxsbtVykAnV5a9j73Tt5sjA2LPSmgfNWOBQQEfk5EXnn2XYrGZQsiOHQZShMsDECwMBxXydaFwuLQHWEwaDEEEwD7OiECHGPRWdlRRdOnSoRES+6WWivr+VjsktCUovZQp7VYK4h+x8FUaEPKyvZotFirnt0Uh3NqaNFXy7i2z51SlvDYXGXd75QgyNg2vn4LD7wgIbhUNunTml///6xW6C/f3+9wOho5CYs2LyNod/Xdm5tSG5tpHzPxrDmQKDwXG51sZnwHo/rsWPjIHUJpxwPT9lpQz0LCTop1IlzMRyqiujSLbeUyrvCKBdws5oJloVTpyaaeG75sGY+zqIDd08gfFBop6wjjAF6bkm2nPny9q2YlwHf6xLLhjRrdvnGGwsbUbcqGFT1/BQOkr2R70L4/mci8n0i8ksUkP5AVV079bKfvqrLlLpra+WUttFI25GFa4sJfd4poildRByogaFmuqyZFmHB8X7+e3nfvjERxfY/xHDmANu4HYXURNYA1deExwwot4jMFSM6CdTHGLBlii1qMUsLrQIb39ZgMLZM2lDemAFevHExDAaqrVa2FHq97Mq/ezEQO3YaNUzOBMPxs7FbgvGa3dwsBPGRyRZoIOJaw/rrCDA7k6hFONnYoXXgMfYePMNp1qiZ27x2aGyW6TfSz9LNNxfcQ6LxXdQ8FskAdP69A2O/eN99riWGdFy1TlhLH+M5Go0FqRtgdsaW5zLFAzo6EXDGdxY2NsaKzFbhfBUOz89dSXeKyJdE5Gfz+3tE5I/zVNY/FpGnVdW1U5ZDVyfEP6OTBbr4wAPlXOfc7RBDyphqTJNOEkbNC4VFRyfm8ti0Ho0KJ5SmzNygmabtMlTTHoG588WbhGyx2QJmbdw0Jq7HmJE919GyRoqbo0wrXLrtttKCmh8MXOaKlsTc+rrq8rL2r702c9GtrRWELY9RzGftpdeKppmNWRYcf0gxRLMaUZi0qQ7bmc+a5diiUohrgUbNbhKMNzEjtvgCPofzZvMbi/d4StgcCU3U1FtgIVj5jhINkUvJxgbbrEvDqfWJ/v2xpZHj1xoMCrSJcS+7t3T0aCZIEskhOJdM/+ym6x05oluF81I47OS1E5vgCgQc8fXh1drY0E6C0Gxx1enAsk78tHtOnarUGGNXizQww90YnZcdhVf3+HFd3rev5LNdOHVKZ9fWNAyHW8JvZmNDOysr2r/22sJCX9jYGAs4XkSdfE66a2tZ34ZDXT54MHMljEZjBmiMr7Zllmey2G+roxAvOnlyHIjljLUwGhUExPLBg9m4R8aD90NYkF+0rH2PaWY00hYxOj6zBy08Gws8+qHwxrNcAHRsLG2+FawaG4+1NVdBMK16D/QLM648vKIb4WgO7FpYXdXO2lrJIrN22lpk4G2qv0rJWnTG27sCtNfRMm21NzbGCgQL8jAYFBRDs9RwrJZvvDGbXxuLBP7tfFxRuSuUGw61f80147T6aaERDgno6GQSCsw8xvxA4xn7+adglPhCHSSKcV7+sWMFMzvmAy2lyeZmcyE4PRrp3Pq69q+5prjgh0Nfqx6NSs934bt9YnC51gUByJ5mi6O7slJgvLjDGt1xmGnVWVmZLCqYt57WTEF0hBrvwJ3TiZA0YRsGgwIzQ83RxrV96pQroDDwb1cHPvG/ksWAWnIeF1BV1X6/MJ+plGtRndApBG9jwV4T5JyFxuV5oxtbfNwnY5L2u3CMyWCg3ZWVSUwGLnMvenUK/MdxlNJ4As4lpY3GD/cjxZSO3qFDUdejqJbiXIWsKbCCZi01F/qLv+0dJFGBlsdttLe1fKVGOCTAiLk0ycOhhtwtE9Rh6p4GHtGI+KpKCWybdq4TTRM1EFyEXcKpd/JkxvCGw0mgLc8JRyZiKXdsbbR1wmw6x475+eynThWJH/aFIL6lBWM4Hzvmak6G1/LBgxNXBTCM7sqKa9lIpD3MNDFNsC5R2biMU0o/97kscL26mgWudeIisKC8qr9jFudoEXCJ4uS5wgaDcVqt9nra37tX26dOaRgOdenmm8fBfrc/uUY7HmNN+L/J3RNGI9VeL8t4qjFuHRiD1nA4FkzmkvOs6eWDB1U7He3m8YAQW1vO2PDBhpyazG1V9lsddxWVa6+ujmNoZvGX1vDqqs5sbLhzuUh0WWXJVCVu4IbcaaERDgmwRR6bIGMSptW6RArf26ur2l1Z0dm1tYmwIALhs1eQaZsQ6BBOCx7eFveg+r2+LOSMxBaT7TjlHdNz1p/cjVPLXQO55/Z76bbbSoJnXDeMoW2M26OT4KhZL+ZG6agWYgGIU0oL7mjRvHdxrxDmFsQda5Orq9o9dmw8ZxiP8TQ8OyqbA8GpsUR6mMv3BojqeMf+0tGjhTTg3pEjY9chul0WiF7Hrqn9+8cxFvfcHqCn5RtvVBXxXZZ02RoZu5aGw8LGx9g6s82jlkTQ37u3NC9hYyMTcJ6Ll/uc3+OEka4W3V54tVdXM7drXgdaNrNkrcpoNG4nlqHV378/OmYtqCcMh7p09GhZaJpwtnTm0WQDZXt1tTAOW7MZMmiEQwJikt8mqH3ypHaOHRv7Pm0SY6ZuZ2VF+3v3TjQ1WJwYOEaG0d+/X7vHj2sn341rhMaajzFLw3nptttKbhZu0+tXwQeObjJ6JqqJWvv2efRoIYvC6jPNaf706bImBtYG/m4NBuN9HxiMQ0gdhOb1NSoYRsU0xvE8euPmjKlp7WEw0PapU7Usk1LAOdYOXLNrayVGg8zIUmk9wROIPsxVZ+nPnYrxWz54ULXbLW1itHd3YAZUNL6Qf59XWj9eUHY00tm1NZ1fX9cwHLqBaKvL5mD54MGJ5XDq1FhQdykzqaU6iRNSfVHBnVsK3th0tJjNZ4Hm9sZG6ZRgvJY0t4gtu2kwKKWdl9ZJHstAXMztZfO/FWiEQwJKggE0W9QOXSGSlysEnI4dK/hOTdtZPHGikHY69qvnWp9q8XTQnvrE5eaMG2PVciYDEpI3eHPg80RGXdg34VxIkNGMFMDFMmSQeZT8xBFh1oN2YuMSY0qdtbWxRhjDbWFjo+j+cRhSyI/74HFfcBiH5/92mRvHQCKBfpsXy0JbOnq0ZKVYLKGqXaPnuourpTpO78V41Oz6emlfSZV7hOnHNqDNDwbus73DhyfxLbBO+9ddN05PXbr5ZtVeb5JhtbaW7alZW8v2T5Dl3nXiVkZjhew0mIs5p7wn5ANs9jReMX/6dGkO5jRTgkyBnF9fL2c+Op6JluqYls3awvWxFWiEQwLMp156W9dgkLlidLJIjCgwJoBb30uMeTQabxAbZ4jAZrD26mqWoXHsWNHcHQ6zhVhhAfBlVgWnyeEpj/xMG/yiM7kgm11f1/7+/RNLIHf74HMdhYCkx9DoN5rftpDHZj4Gu+k51IzGqXz5prXxc5jZQ4tq/vRpXXzggSjjNVwxS8pcbhhXWcg3yfVOntQ9jzxSdKU59ZpbYjHvJ+LXXVsbByJjfmmrY0Eh5pS7eLTXc+cyejl9NwHW0bLlYIIGaaRwlAf1ecwIp8CJ6aNA/7lbs726qks331zM3BpRQsRwWMwaRMEXid0s3Xxz5vrNd7kb/gsOXoVna6SfLpw6pZ382I7u8eN+eaZV4BNeVtbYRa3FTanGdxrLIXGdkYP38smY16I231cdn1cUO1oBr/bq6uQcHt797Cwwl6hrCAd7D29fteTSMould/LkeAfmzGAwDiB77XTyzXQxXGxjYB08Z9T3y3pn9Xtpod18zPvXXDOxbnChejhGYkJ4Yd9N+xrn7ucHt40tEKCZwvxq5PgVWPgF19lwmFkbHu6sYbJL6NFHMy2+3y8yRBRUHi4R16PRhMW4TFh4ikQpiD2F9eGuC67/yJHx4XTsukvW5bil3HKjUUF56xw7pv29e1V7PV8gOAJ1+eBB7V9zje/KdfDprK2lLXCq31sjLS1a5fgaAczq2yo0wiEBVRrYOJagxU0tqbfELaoWmMPiAw9kxwrfeOP4eG8kIu/obU+7GBOfEa7DRD1mu3zDDSWCLmULUV3tkyf9lNXhcJw1ZHsQ8D9vTHDsYgsKLRhmZm3VsWujlO4YGyf7XsVcRpPguxKOYTDI/PNa3JmOz5vWVkVHLXadAL4F1wN+ooDJcTFok+88DAaZMCO69eZWNItXzdr85YLCDthjJtXWiT99Nreqlg8e9F18NRSZhdwSwXcxlDLnUEDo5Oyk6DpJ3cvXIJfrHTmi/euuiwrojpY9CgunTmVjX7Pd3d/4hnvfE3rtjY1MYcvHwRIIzCXboXlt6fZOZFXVRjhUgce0bOAtnTJ2tvqYKFKLIv9vdn19zJS9836i9cVM6vz+gpI/357PU1h7hw+X6qliSN38KGQPv+7KyviUTyRWs16YgZqW2F5dHWfcFE5vHQwyDS4iqDrHjo1TSb1xtlNf7cCzwv+ooafGPBdszOzaeTwC7/EC7e/fr/1rrimPY+5PN2ZaODoCmANmr3CAHMfDjtpWVe3COU3jzxoadyefI11eLpTHfSfGjG082PUXTJCPRuNssXY+v20ItMauzrFjqiJjK7qztqba6xVOoOV06G7eb1fBwHFnGslTkTGwHYZDnd3czBSSiBsqDP3Tdsf/j0bl+F/e/pyHG1xzGtkgS/SPR9+YNcf0iTHLrUAjHBIwPj7A27GpOvYd4qYoUTheI//kjUPuFWF+qaurOvaZL91ySyF1NQwGYz9pLOOks7KiC+vrJT9y4eK0wdzf2853MHP52fX18U7lwnM6yeG2dwIXUk2B+XRgTMLGhnYfe2zczz2Em5faiFZB9/hx7R06VBZmo3imiTsvNa9ZpeNCQNuT0eSEWoQCkyGLZgnngJUCm2stu30WqAyOGfqp92jxlNT26qr2HnponGYahsOyMKmjGeeXpcZ6AXvRcixi6d57NQwGOp8LXks3Fi2e04Un35orqIPWqtPvEj2PRqo68c+3bd8Fp5nmAkNGWTLG0tGj2XlJHM+KKCh4r3fkiGsteGNnO88XsNxoNBbK7LngOs01tlVohEMCxjtNK7QdM0FtERWOroZNYzHXQW3BAD73jiGZH2/dXVuL+1hHo4LGZwsAGWatTWDEbLzNWqWFpWS1wH3TABc2Ngq7QT0XGGtMolkmR+WChHkIo9GYCVrqb8xKGzOkkyd14fHH/dRK6F/0dFG78vZFJ/GpTo7H8sGDE3eERwsxHE+ezOYQfdla8WIkqss0Ty9hIkn79H/g4DkzTRyv0Ujbp05pd22tcDR8m/Dzgthd1ajCFourdNbWshjWYKCL9903ppk9lGK8dPSo9g4dKronR8VECDf9GVxBbYgjjLOWYCwKsboKixU3cXbov7aWhWbh1N/c4mosh8S1Xcuhe/x4/ITVPIDYv+467V9zzTjbyCaypZMD6/rXXaf9a691jw0en8oJvv6l224bv7thnDliqWq5WataYUrT9469h2BzU5f37Su80lNGo/ppjETc6PpYfOABNz00DIc671gUhUWUX7Z7uWUM2QRDnh7Z0Yz4A2mBdXC3/RHmEhz71qnvfZ1sutJWa/LGPBQ2ifFqDQZFawzws3cc2G+zbuoIhMK4DQYlRlaK9RgDStQ33hzGAhAZe84kkdlZ1l5h/GrgbUzNhGSJ4UaE/czGRrbfZ//+cjA31u5wqL0jRybvWajYsLewvl7IEutoIh7GbSOtqh9rwk1rti5wzizJhQVjT8vupq43djk+e772Ne09+mgTc0hd2w1I49b7njpn4aysqLZak001MPFLUE/n2DHtrqwU9g4YEXmatZ3Pok6bMsrS27TfLzIoj2EBwS6sro6ZXJddXaPyWS522dlOsV2opQUQ+y/y7GzkuAy+MPC7RHXOwKJ3j7mGi1ODuW1z/fROntTlG27Q3kMP6fL115fGuiBcEZ/hUBe//OVKZh/yvnQg+Jvqv2XAxRhpyIO5lZYDxyB0IjBZe0fGh2nXovkRKvfeW1RsaigYHZifyv0PI/Df57jYjvgOC0Ln2ZlcUcHjxavww6sq43A8N+RiMqa86Mw54hc2NnTWU6bw+2g02ZcBz+L5VaX5MgGuWw9KN8IhAejHs8ninP72yZPae/TRqHWBJyYWJi925QRjgcpu7NlcI5oHhm7+ZlsELXquY2mrhw+7B5mNGQdkFs2rFl57acypZVlVTn/aTtsx11lLdZLCS4xlUWswD81PZ4V6LXsjNcbjzyqG7B2wF3P9xNwqkblH5SFZfjRy3XVcLgyHYzdC6iiO7spKwf0hoyz+Yhrs2J0xHI7payE/GLEwzxi09fCPufyqLqhrz8MPl3emj7INbmwFcKYf05mleaLw7CgI5+FwcjpxCi+H/gpegdx91L/uumJyRVV/6SoICOrrwuqqaq/nW4pU51ajDo1wSEBp8w1+atEdMva9E5NIaTeFg8DQZ2uLPWJe42fn2LHJgXo62XNh7hfE1fZhFPYEeHUrvBjn4EHtX311mTmZSymSJ78nP6+pSjgsnjiR4RNhfoupyc0vPka6r5re+ezgEWVi6FqpoIfOsWNjGkhtYMM6xofW5ULZXnzfPX58HBRevO++jPlU0ANuaEwdtOe+1EazYHQnf7a9sVHYQIi7nQvKTmJPBjLvWAypcm5i8+Tcn8W422gU3S3foXvzp09n/+XWEf+fwnV8AoJtXqwQ6B11rJeYkrWxkVns6O4DS66zslJIUUarcI9O4jiN5RC5duLI7iiBgMZm2kOXNrKNNXt6nonU3TE5HGr75MliVgRc9vIQVWe7vvN7+eDBDJ8EsZurwyyWQNkt3uWeP1XhUhGovxSkozILTjsmvJa0+OIeO0+ntvZnwjzFrHjxgsAIGxu6AKdsGqOpfVG9doaW9vuF/Sbs3gjqv7fcrIBUm6VAb47D/GAw0URpPLrHjqkuLGQnpHJWkDOmrc1NXTp6NDsd1hvbOlZFQmAv3XZb2Tql3158ob2xof39+6NWn41rqwau5o5DYTnH62s0Gr+HQ1S1HXknRkzZ6PH/9N07ARnXSrMJLnHtxJHdBSIZDscaVWtzUzsrK+N9Du4GGGfSZxVeUWj3NjfLhKXZGTLjvQhY53CoS/feO36LGp6gGTtjKOkGya+O5htrOODr9SUmMPF/TyNSWFCedUH12buc8c1qtuGw4+AX2zA1ZpyO4Ci8KAaOBuHPGTsBFMqzhTQ+1jynI2ZG42NMrP9en5WsSbgWFY7ziMzP+PTfPLd+NqJgdLjtyJy3T55UFdH+3r3aofdt8PPmsky5UZZuuSXzowNORsulsnwk+Gg03r3dU7AS87pm19Yy4XTzzZNTSvNr+YYbsmPGLV4UEUBz1H/JU5K9U1Kxn3zApDc3C3Ryqt1vbW6OEyTCYKD9H//xYhqwU2/BijMlqbEc6l07FZCW0UiXP/c51VarsNmsd/hwpqnGCA2l/HBY6Q+fyxlCazgs5K17qabjF7wfOVKoc1GdQBozbHAHzcVwypnp2Fw27WpjI3NlOa/ftIWUshwsrbZwHEYdVwL8tnnB9zLMrq9r7/DhLFfe0/a4reFwvNDndKJtdRLvlIgKMvQ3f+5zxbeg9XqldNPSUeFV9da5UIixMHDcP2E0Km4Cw3qGw2J6aj7/cxFLeCyYh8M4XRDtdvHgvNHkoLxSgsZwWAxKA75GS13Ns3lw3p304/n1dV/opvDFMRsMSq7i8RsQU/tBqH52QY2tt3wz7PK+fdrfuzd73/u+fRPFkeopvBgpQitPqpiDiFwsIp8Vkb/M3yH9L/L7PyciR0TkC/n1+qq6tiscCibhqVPZfgYIUi7v21dcXJ6QyJkxn73PZQpvvBpl2UOtnJC9dxPMDwZZ7vaJE+XFx+2bG2Qw0Da5iYJTt13dlRXlM/u7a2uqqr7Jq5DVxJfXhmdZsHZfYYnYfzP5EQPGwApvYMuZT1IIqZaswxKuqX5YrvvqatGsz8/bSjGgxQceyN5m5+wXcAO/sT5MMc6LDzxQzlqDeY8K7qo2YoKOro7ttDfBsrmp2ulkzJaen1tf197Jk1HFKnW0fhSnOnRlQflU2nSqv1g3fJ+NZagBDfUefTSjoUcfdY/557aWjh4tb0AdjZ5cloOIPFtEXpZ/v1BE7hWRF+fC4Z3T1LWjB++NRoVdlP2rry5l/ZQOTqPJbIGQwKudb1rpellE6iyAEaSeEqHspnKVBKyTDVnjUx/NSsgFIJrLbdXyQWhVTLeuZqaT45rdchH3iFdfO99k1dXMtaO9nnvEcuGaRlOP9QPwi7rAeP43NqqtFW8fQk08okwR79mca/ElNgHfL15H+KSYLigLS/feW1CGwsZGVFjJaFTeO1KXxkajstCNPRP5PTOIpHGPRqX3YmC51uamS1N2hLcJnrn19ZJb2tKc+9deOxEmPK4OreH92fV13SqkhMOMnCNQ1aOqekf+/aRkFsRzzzYepQZDkLV2W0REdGZGfvpXf1XWFhYmf6vKD37qUyIhjMvjs6Iqb7nhBnnFl788ua8qIiLr8/Ny1UMPic7Oytzp0+P7ojr+PiRchrOzLt6PY915uwVcsO4cTqjKU++/X/YcOSL/5N57pRuCBFXRVks+dtVVsjY3Ny6/JiIPXnRRVif21eqze4A7tlUCwG1ORD72wheKzmTkN8vlZmeLfWGAPq/Nz8sJEfnGaDT++0333y9BRNrr69La3Czh1l5bkzAaSRgOC7jPisiCiLRFpJt/xtpF/NZE5GMiclhgXhxYn53N5tObK/w9A8uyYhxmh8NiWW9+uC5VufSxx2Rg7eQ0IDMz0hoOi/WMRsU6bb6RJgjH7mOPje//xWWXTdZPCKJzc3L8aU8r4gY4DlstOeD00+tDAUKQ4dxceYysLNeB/cphNDMzrn9WRGasTAjZ+GB5wtlbe7tmZiQ87WmyccEFIjMzsjk3J2u7dxfKHhcRWVuT9/6f/6eMWi2/n0RrH6cyF6bW3XYgJaTP1iUiiyJySEQuksxyeFBE7hKRj4hIzBsyvnb8ZT90r/BycC7jlO8dOlTOJnGsi8X77ptI/7zcDPvoTZsjl8O8h7uHT+Q3au3RF+7UNMm9/k11pWIH1L7lsntl3fdio0ae/27bS5dUC7twOzo5cddtJzW+Nfq/+xvf0Nbmpn+4HtTXOXZMF1ZXtb26WntTYgyn0kY9KFMKug8GumBtQj1z+dlcdSzIZc1cH2PLLpXuG7k/q5TeXEXTU15hVP0yq3G8LNZeAofCplejvUj5MKIj3fn/Glbuk3aHtIg8RURuF5EfzH8/S0RaIjIjIr8oIh+JPHe1iBwQkQOXXHLJlgfHNWGRGSu5cKquUSRbxjMTPeZblf4H9Y43onk+zQTBcf+q+hNj1MkyVXVGcucL/lSq19J6vTHprK1l7qQqNwv8nlFwrQyHyV3Hi5oFYmcoa2bxgQcm4w14z25ult8CxoKKxyT/z94NXjiGAVw1VWNcUjC4Dfi0tOf59fVxijaf5VRHKNjV37+/8CweqmfXWLGpoqvYGoGrreUNku7z9H2O6k5uXkvRcx26d/pgae6FzaEwx5aJVOcFSnaG2FbgvBUOknkY/kBE/lXk/0URubuqnjPxsh/3RTepy2PQDmGF4bDMAIfDpPZQuF/DH9/WyZlRU2tdhjO0E9NeYwuqdEBbqj/E8DxNCU+J7WhZo5/d3IwfqT7FHPIRDgXBFKmvJJDyqxRYdQRi6WhnYh5b7k9MgDgCxqyE1ubmOHOmv3dvgR7Hb9KrgUOXzhXClNCpGXCFsAi5oC7hx21Po7ikcPDGkxMzHPosXHlWljH05YMHi+srV5Aqz6FietOtQUo4hOz/NIQQ/k7u7rlUMndcyDxS+vyturNCCEFE9onIMVX9l3D/2ap6NP/+UyKypKpvTtV1xRVX6IEDJU9lLfiur31Nbn3Ws8q+XQU/aw57RGTFq0Rzn+YMhXBGo8k9JZ+navwT2x2Nin5dnfggg4iMa6X7zxWRVRE5Yf8NhyLmG+U2oY6lP/9zOXLxxXLqKU+R491ucky6kvtM4f/OiROyurAggwsukHHNw2EWX/BiIjHfO+C4FIIcEIrHMPA45RCsbW6L2xOR+eFQNgxPGs+F06fl1Px81BeM+M5tbMjuxx+XU92ubHrt5Z9tEdkUkYFki2oA5RZC8OMXSFNct/UxMo4xfEVV5gcDWdjYkH/w6U/LTa95jcxubsqDz3vepKzXpojMS0aDm3m5MBpN/PM8F5Ipo8k4Sj4WQ3Fo2+JKMzPSPXYso0+nv7MbG3Lh5qa84fd/X37n+75PNkLImF0IoiKyq9WS07t2jcc3qMpMCDKUzGXxbMnWz608voaDQ2cuxOYkh95Xvypy0UVy+ClPKbYRe15VwszMZFxy6IrI+0XkbdUYlSCEcLuqXuH+V1M4/JWI/JRk7p/xGlVVl1fWROpVIvK/ReSLImLRxH8tIm8RkZdKRhsPisiPm7CIwXaEw8xoNGEeIu4kiogsxYQIPhMTALF6VTNiw0BUAge83xaRN0oWCE2CCS7DyRMyo1FGcDMz4/ud48fl8YUF2UCGSBBUZVcIsua1FxnTsLkpigFnT0jwc7FxtLbqLtYcZk+floH1C+uLMd4YLpKZvptUtDUcZkFKfh77mRKOwMhcPGLPcTmbhyqGPI0wESmN04KIPG11VY5ccIEoBoWp3vmNDdmYmyspTK3hUIYmlEPIkgk2N+VxqCvkTFlz3Jbvu08+9sIXltsajaR74kQW9AZhNZfPuYYg7VOn5PSuXeMA8LJkQV4c7yXJgp5rXH9qvEWmouPZ06dlJgTZmJuTmeFQRq1WliBiZWCcWoPBOOC+IJPEh7akkyCqICUcqqwVyYXHrdNYYmf72o5bKXo+D7uDIucCmcvDe37p6NHiRpiqdkajdHl0F/BLSFSniiPs+drXxr5y122Q8rWCO8J1FdTtbwWO0fLTusro6tL5PKV+p9qF+8ta9rGLJjYnOnVU+tlj8QlzP3nzFHMn1Rlndvd5Y0Lptm3NIIp/xRj3Dh0qnXjcOXZsHDuZGQxKr6ztafmAzPFVdd4XjSm+pKpyzKal3/x7CzccVs059sHph6Vqh+EwS9/eBsgOpLJ+NoTwSyGEV4YQXmbXNgTWeQM/COmbIjLWuPasrEgHyqlpYAA9EflPMzPyrN27i2l5mmkItz7zmVnqIlsVCKhphDDRaLEuqHP8DGrf/JnDrNXBnyHI45ZSNzMjD156aVnjqdDE3/rbvy291VXZxX9o5lbhe8l6AbeFU6ekc/y4tAa5kwVTDm0cMAXV6tRME4uCZlZaazSS47t2xS0Ez6x3/l8SkZtE5HMPP1z4e0lE/oKfN3fEaCQzg8Gk8Ggk85tsd4gE1DJnZmRWMgulgEOOh4YgexzXYDAXjDpartdf79O7p7k1AjQdTp2Si1dWxGyl1mCQzV+Of6k9qu/wc58rJ8y1kt/7B488Is/5+tdFQpDnfPWr8oF3vUt6X/2qLEnmnjom5OJVlTkbW3Rh4udoJLPr6yXrTUXkIaprbMWMRtKRzLoY9ykGqjK3uTm2crDPw5mZzKqiuQqjUVbeUpJ5jHHccrwsBVlF5F2XXy4fimO0LajrVvqsc1tV9f+18yhND9txK10sWX76mCDsM+XbDUH2iEh7dVWu/P3fl4+98Y1RE3zu9GnZjLlmvIWbWszIbGKMC8vE6kvEMUrlEm6Jroi8QUQ+JVJ2LSFO2IbD1K3M7OnToq2WXPEXfyF3XX65rC8sSGs0kkFsb4CqLM7MyFdEZNfjj8sbf+u35GNXXZV2i1QxS55/zfzRo0jxghsAn099z6GzuiofeMpT5J+LCInTMk42XjWgc/y4nOh0qt0gMUFYhT9/Z9xUJQwGE/dhrK0K19iC5K6ezU0ZjEYymJ+XGdXMBeVAW0TW7RmnztbmZlFZi40nM3dVmR2NZG4wkLXZ2UnsLvLs8l//tXzsBS8o8o/RSGY3NmRwwQXV9Glj49GsahZHAxx6QsJtCth2zOF8h+0Ih6tE5GNVxAu/W8OhXNhqyWOSEyELEZ5IC9B5//M904xjDHlaYcLAxFaXYSYgCPk9KxhiCbwAX6oO+t2WjLEO7T+bD+5XhaArAM1/d309szamGX+cy1i/VGXx/vvlwcsuK9QxIxIVRpUwGkl3YyPDF9uMMJokPVQJjO2CZxHXeSa2VqvmONWvVBt1/8MyxsTZc1Dn2ZTgBQjDobRPn5b5mRn5wK5dOx6QruVWCiE8NYTw70MIB/Lrl0MIT90CLucd3CQSl+DO5AxnZuSEZIKhZD46hDPO3IhpS/TM0m23+UTMz+I9q8cT9LH78OwsujqmBBWRx2F3MhL2/MZGuW2vT9yvOosp/39NIEMihKJWFzK3QF9E5mIZSyIT11V+zW9sFModn5+PWi6lhYx4sysSn8//LwiGvOzIG7PUPALMbW7KG37nd4q45LBw6lSxMAsL77ddKdqK4aUqy/feOzkNwCs3MzNZRziGg8HENYZ98Zi7J1ywrhjNxWgrxcBjY8FlZmbKlhMpBi4wr0iMfWtzU9YuuEBO7Nol701jtCWoG3P4iIicFJF/kl/fEJGPngF8zjpcKZn2Oz7CwSbD0/hokkpTBpM/e/p0kbjteY+Rwu9bX/GKNMIeoSQILwwG5YVFWuBgbk78QzpqQgjZMRWE28YFF5Rx9xYqLmQ4sqDEcK18DLD+/FNDkA+LyGZiMc5tbhYYTKkv6P9F/FiJ8HBJCTdWRFLMihhfGFJib47D5txc5ua058xCefhhebzdLtWzfP31cZrE+TDg3yFI5/Rp6Z44MVEy4P/fvPRS2TRXiqcEiEziecjkW61xZlIBL4WYFuPHQsIbX8RdVfbQ+Lk05wn+GHgCM6Z88nMiRVqLKZT5vcGuXaKtloThUN4Tx2jLUFc4XKaq/1ZV78+vnxeRLe9xOJ/gJsmY/ECkzJg8YsaJo8Dc+H7OcEv59aqyfMMN0j1+fLyQStq1pwkhcVZpNWy5zM4WUwfzYO4sBdee4tUXYah8JpGEMPHlIi74ad89fzGWnZnJAn+8sD2hwrhG2r+V/0NLIQTZtJTJvMza7t1lHI0RxZgc3BsHLiMCLojILFknJYbmASofbO3lbc174xCCfOXpT3fp6dV3313so9OvgoXM5z6piqytyVNDyGJDCKpZvI3rjmnyONeedq8q85ubMmvCMQTpnjjh14vPspvRyg2HshJTnDxlLgZspXAd9n9eV0BlBPGtskgiY9jZ3NySS6kK6gqHtXxfgoiI5Jvi1hLlnzBwJf5ANwAvsIT2MLexURYsngYwMyMf+6f/VNbn5mQwOytBSLuu0l6QWGMaK/4H+IwXeM7E5+HgvyAi30pVhNFowqSpLSXXjUjG7GbX132NyNHoUzC/sVG0uryxJGabhBgzZ2vA6iUloG3PcAaJU3cYDOSK228va6tQn4pkWWlYBzMpBBzDnLYKz0M9G5HxVdp3EYZDkZkZ+en3vtcX5ICPonLhMOwT3W52SCP3J3aAoIcj9j9mkYUgG/PzsrZrl+jMjIThUE4uLPj4e/Xz2HJgmdvFOa8jtBF3/CRrs7TPJ7UmaiiOumvXGclYqisc3i4ivx5CeDCE8BUR+Q8i8hNnAJ+zDjfhDyKgVsxkpcnctI09nnYyGkn7cdimMjMz1kxb2CYCLyQjCk7R86yWvA1mcG/99KcL5TZzAg2jkewS2g0qInOq2aYbb7E6xBxaLRl6WVnAiMNwmLm5KmBtYUEufeCBElMtjU1+X1lb4/HA7954485xLh+CJBJkC+Xbp05JZziUW5eWyozCwBEWUauHf+d1tkYjufQb3/Dx5dNDPYskBGmfPi29hx6S4KR2RoG1e24nxvTp/qytK15fqTFwcNFWK8v+8XD3rBOuNyKAYu0l68d55LZwPVatd77vrTtq+4TIuYs5qOoXVPU7RORyEXmJqn6nqt55BvA56/AeyTVPnNh8skbDoSx/6lPlRciT6WlI+cQuPvigBFw48PyABYmnuaQWroeXV25mRj72gz8oi1BuM9/focHZ4Swim9RuqFh8m7OzxZ3mOWB+fhDJtKZYP6DfD152WaHIUsh2zrZQC8vLRrVtb9FVabAOEykFiL1nQ5C1hYVJllCsb/bdLCNmil47xMQ+ODMjD154Ybx/TDMObTzebsuxPXvk+FOeUha83GZC4EcZOjJGwmVkmrN3pIqnkKUgZjV4zNxT9KpohMFwzMuM9xLVpcMYnp7ANeuL+YP9B3BlHOMtQzKVNYSwrKofCyH8K+9/Vf33ZwCnqWE7qawimYQcjwJNRHswkOFolB0jYfenASaOFOGzOeoJj1i9sbrrloP/wmgkTzt+XFb27KlfH2uoqtIKIX0eUgwnamNJZHK2Uh1rIoUXPmdg9729Laqy5+GHZeWZz/THIdaGhwMCxi8SfR//BtzaIrLmtU/thsFgcpxFrFwKUvQH/QiqxbOzUrjF+srPOPGd0j4GqLO0FyU2jvB7bjgUHQ5lgJvTUmPq1FXa5xIrz3XXmY9cWLREons7DLa612E7qawWmbvQudwY5hMNPiREbCTl1+bmsrgAEioGNEV8DRTrc6yGaFnvE7+j+WqgKrPr61nmCdwrtI842uKLaDnaamWCwdOGeMGpTnZ3Ep61BQM9N0+xi1tHo7JgsGdiDD+msWm2gxjvtQYDCcOhLN1yi2vFrDzzmfG687FPWipmfeF90prnLDEhvzc/HBY01HnIOlvD/njzkd8fW2kpJs20zP+rFncGc7mZmfKhijENGhljldLkWBVvvPfeCcOiNXXR+rose/XEIATZnJ3NYjce/oizs75bm5vZjvuNDWmV/i23lVzbCGzhzMxMBEOCz1xZhcMWoPaprKp6c9W9cwU7skOaYVpNqw5R1q0zoQ3GtOYwHMqutTVZYzcBPsPaKrUxOxjIAIN0VXjEtLzYWPB9qmfu9Gn59QsukKurtD4PUmX4v9gO8ZQGH2kjDAbFAH0Cgogo5vQntMH26qqst9vjfTKFzZRbpTVn/Jc/+Um56R/+w8nJoF491FZLRHaJyCkuV5duvO/cjke7IkU3FQrX06ezlNlU/2MQo0uC3SJyqq7mX9d68p6pgx89fy4sB4P/p+a9JxxcSb+j+f4p7UqkyGxiz7EWj59ct6WLwv1ZfLUoafXaahVex1giTM3SAEuWjKEvkqUixhYv9ofqXrBXS+L/MS06wYRnZ2bk7Zg26S2imDKDY+tZSo5FUOqnp+FxG4SDsjBNgIqUkwUKBSZ479rYkLd+/OPjrK3CZsoQyudXcT0xIA35Y29+sy8YPN9//ns0HGZM0qk3iBQzzRgfby3kGnL0PCETpF5sL/9/E7O3PNqtGpNYuwDr+F+d/nl14frIz9sq9Ynr9vCjslfGe7dlSO59CiG8UkT+tog8g+IOF4lUW1NPBLiJfo9zaZgImMF4C9xZTGN/LDMh/PS+O2mTg/l52bOyMnH5IF6eQCBGvMHpkwXUa2h7EdwfTzFfry4PZ5Hs3JoqrSwmKLz+p/qSC6GCX57LonbvzV/qu1dfVRn7X1Ueu/BCue+Vr8z82dyuqszOzMicFI9rXnj8cTk9P19877jC8d8sEFNj5W0CzctyWizWoyLFjKnU+HhrymPidbRpwoO/73n0UVnBvR6xeqCvYXNTZkTG4znK038La5/H0msf66eybjyrqu8O3p8Qkev9kluGKsthXrLYwqwU4w3fEJEf2mFczglcKZmUK71MXmQy2R5Tj2mo9HwpWFXFHFCbYpiZKQoG1pSZOLG9GHPL/ytkI8UEDGtAdhlT94SaalE7ijENHmePSaRwxLq8Z7x5COSXh2cWv/IV6XjauTcG/LwDpWwvZ34XoT+juTm59QUv8PsTsljYmojM5BZmGI1k7YIL5OIHHyzgNRec90I4wsb9hDbn8MRb6v/iN74hvVgb9t3RxF1cSGFwlQFvfcQsQBFZ/uQn5QSumxwWRYpp5polH1g7Ojs72dxptEKB9/GpuinBR2toHjLVxoKhgn5KQH2cKr5Xt4maMYdLVfUrZ6D9HYHtxBy6kuUJj9/EldKmPPC0MCKK3le/Koef85zx//iyjmQ9qfpTmjn+NjxGIwn5ZymA6LUDz7ZXV2XPiRNy+LnPTbfrPMtlwnCYCcxU9kXV4ve+p56pY03EPvMsobnTp2UwOyvaasVfFpTqD1ohEdza4pxu69GCFzuatq8i5TFEl15ifO303CEI1u7p03IcTyxlzZ0t3qp5SIxB+/HHy7vYE+sj5DTfEpEPisi7JHuDYRApZCmG0ahImzhvzOjzz/ZgIGtseSb60hoMJMzNSWm3D88NK3sVde+W7M2P08JOxBwez9/n8D9DCH9i1xZwOe/A9KEB5p07mlPhO2pBMW0cfq889amF53V9XQTPPIo8h/gEr00HlxZnDgF+/+ld75KOt3XfIzhgPmu7d8tzDx2Kazcp7S+/9jz6aNb3OsHbmIDTzAqZPX26uE+AGUNK4eG5rMOQRiOZ3dyUXWsZ6x7Mz2fvZQCBVwk13h+xVmVdGC3gf7HPGG3xTn6ec67Doe//eOut8pYbbii0cXzXruI7RqjdsWDw1lZMsOefgc66Wtu9O32sBa0LlSzZYjQayef+y38RfewxERHZtbpaoKMSbdL6EVVpr64W/t9gF5uVhWewP8O5uWx/UwxvaxPXUKxugPXoYGwd6gqHj4vIX4nI80Tk5yV7fedtZwCfsw6FEz297wbApGbzy12ATOghZBlERlynTskaLyTP6qC6S3nkQEjLn/zkJLiFqYkY8ApB3v6+98kJznhhvJm55m3c+spXxtNfmeHSsxLCxN/rMeLY81ZHXibkWtxg167CWOy23cK4uBgPFjJUd6GsyKSveSAU51BC5vYRzc74ufSBB4QheOm9HlT13ymjnNpYR8nA+3jmF9MgjM/C2po/33/6p3LT3/27ccUIlIGYxl1qE/FDGsHTTXPoHj8+STvOtfHCGWUcF5iZkUGrJTozIx+76ip5bGFBRFUuv+su6fDZTAklT2ZmsqQPuDeMjbMneCv4Sok2cTy4PArA0Uje4mOxLajrVrpdVV8eQrhLVS/P7+1X1e85AzhNDTv2DumaaYZRTZOJHxaepSK2T52aEBgJhbmNjclRHFYv1VP6zmUcPMPmpszMzGTvNfYWQN1+1oD54VB2P/aYPHbhhZOdsDm0RGSGUw5541mVa8OAmMWJiy7yX2yv2WanN/23/ya/8/f/vpy46KLC/4UX48BYeu6vlsBGPKQRp83Z06czISbgsmSIWS72Xz4eC+vr8niEZiTkbyubmZHjXG8KUrTl/Q/QXl2VEEIRJyub931+OJz87/XXICY0YmMruVtrdnbynu7UOEomqFU1U8gQRiPpPvaYHHfmH3Fe+vrX5da/8Tf8ttTZbIh11E3z9sYnNjdUPoxGMqrYJBeDnXAr2bsMj4YQ3hBC+E6RcgxqJyGE8H0hhIMhhL8OIbz7TLXz1t/8zbL0Zg3dINcUl26/vXjIXP45awFM0hLCaCSvOHBAgmpR8yCN6VkPP5xpa1YnaBOzp09L+/HHJ2fje9qWp3WLiM7Oygd/5meyoyemYPbu6z5ZW80/25IRxH9oteTYpz4lw2c9qxDsa21uygc/9zmZnZkpp/B51gb3A8cDFl73xAk5/tSnljNo4NnB7Kx86od+qBR0b6+uygfwzClrS0S01ZKFtbUsLTN/bojCIxf4sT4MLrhg3EeNuEDGLsAQSof/zZ4+nbUjMnmlq4GNXz4uuzY25IUHDkgYjaS9ujqhw4gVMj5V1xNIqPEyI89hbffuIk6kIQ927Sr/j7hbn/P6x1YAz6+ItPIyeIz6YH5+Ihiwfq89VflPP/ET5ThCXv54t+tr91DmyOamtNfXJ2VIQGBSw8Ljj5esDRzHwgm3bInnuNn7RMJwWH7XBfYhH8tWHbfmVqDOu7JF5B+IyFNF5NtF5LMicruI/KM6z27lkkxJ+7Jkx4LPi8idIvLiWPmXv/zlW3m3dgadjnbp5eWiqjIaaRgMii8i39jQ3smTunT0aOFF4Mv79mnv8OHsxfXwPH628KXm9OL4oKoLp05pZ2VFl2+4ofRsZ2VFe4cOZfXkLxfnl6gv3XzzpE56MfniffepLi1pX1Xb9pJzxNHpe/vkyalesN5dW1Pt9VT7fe1r9oL4hdXVMU6dlRXtHT48Kb+yMukL9Gf3+nrWn6NHdRn7Eunb3Pp6AbdlVe06419oJ//eWVlRbbVKL4MfjwHgj/MdBgNdvuEG7e/fX66f6uqsrOjSbbcV2g/Oy+SZ1mQ41N6hQ8X+Qb1Lt9yivZMns74CXfQOHy7R855TpzQMh7qwsaEdo1Gjn83NMR6zGxvaGgx0fmOjSPcO3YbBoDzGsTHPf7dPnSqMv90f0wmOPbbL8450k/elffJktobW1opzlq9PvDeL69SZswXr92CgnZUV7e/dWxzTvK72yZNFGrZ1A79ncCzzdTq7tqYyGuk8zW0L6u+srGj/2mt1OR+H2Zyu+/v3a+fYsSJ/OHZMtwoicoDJyy735rm+ROSVIvIH8PtnRORnYuW3JRz6fe1fe63PIFR1Dn4HZhT5996hQ9q/+mrt79+v7Y0NDcOhLp44oS1VXXzgAQ2Dgc6ur+uC6pgxdFdWtHvs2KTu/LN38qQua7bYl26+Oat7717t792rvYce0uV9+7S7sjJhJvbckSOFxbWwujrGsXfokGoIWX97Pe0YoRuR5+22T50qLIrSosQL7gVV7R4/rv29e1V7Pe1BuW4u2Pp792r/mmu0s7KSlb3uOu3v3Zv9XlnR/jXXqPb7hanpOe11jh2b4IZ4jEbaPnVK+/v3q+Z9YgYSBgNduvnm8Vi3T57U3qFDuvjlL5f6yMx6bn1dW5uburxvX7Zsej1VpInI2MhwOC7TGgyycbj6am1B/WEw0OXf/d2sTZuTtTXVVkvbwMgWoI2Wka/RzI03amdtTbvHjo1ppHPsWDYevd4Y52XCtXfokHaPHy+NJc5tO2dwnWPHdPG++8YKh81D++RJ7R0+rMsHD2Z0ZLSXM+f59XXtrK1p32nDFTKJq7W5WRxzW0/Hj2vv5Ent0zjZM8v79hXWGq/lMBhod21Nl3PaCaORLt12m/Yeemi8/rorK5mwv/nmrM7rr9f+3r0FfLrHj2v3+PEx/fA6Ha/7fNwKY431HDum3bW1goDsnTyp2usV1nZrc1P7v/iLulXYtnAQkWeIyL8Wkb5kb4X7iIh8pM6zW7kk20NxHfy+SkT+A5W5WrLz2A5ccsklWx4cVVXt93X5xhs15Nri0s03T5juo49q/9prtXfoUGYhHDqUaYI50zQtpHfoULYIYSGqZhM6ntyvfW3M5PvXXKMqMtYC2mtr2jtyJJvoVkt1cTFj6Gjst1qqy8sFK2J5374Mx/37x4uyffKk9n/8x7WbM9L+3r3Zs3lf+9deq71HH9XlG28cLyjtdrW/d+9YE+qurGjvoYeKmiZodyG3QObW18cWU+/IkYnloJk10b/mGtWFBdVOR3VpadwH7feze9bHfLwK06KTRTWzsTEev/7evbqwvq5BNVvAhw6NGVXvyJHxs8Y0UZvuHTo0nk+8t3zDDRpGo0xDW1nR5X37tLOykmnbuQAbz7OIarud0c3Bg74FQIyuNRhkjDrvb3/vXm2fPJkJhn37sv4vL2v/6qvHNIf00cn7tKAZw17O59KsNe33MwUBcbRxhXImcINpp6Z4HDlSsIjDcJjNn+p4vPrXXFNg5P23v30s+MfttFrj+vr5fOC84Hi1NjcLSsjS0aPaWVtzGXf75MlMidi7V5f37SsI0lY+PqKZQmHa9fzGxtjK03Zbl/fty9aMqvavuy6b9//yXyZ9aLUK69U0+d5DDxXWoa2/3qFDqvPzmZJz4kSm9OzfX15jjz46XpudY8e0k7fRWVnR/tvfPhY6y/v2jdfN2DpQsGB+8idV+/3xfyFXNsZrewuwE8Lhz0Tk/ZK9IvSNdtV5diuXiPywIxz+n1j57VoO2mplQwHMuL9379gi0OXlCXGEkD3TbhfLmea7tJSVW1rKqlctMmmrp9NR7fW0f801E+JkQWD4tdtZu8vL2a39+7PFZ8+0Wllde/dq7/DhjLF0OsX6REqaeQHy8v1rrskWmA2PavH3/v1ZG9dcMxmDn/iJrA/XXZceVxQEJkTzcSgIDpye2PiYMMkZX/+667IxyS0HVVXtdseM3BZr/9prx2M6njur2+pcXp7gktfR37s3E05XX12cQ50Iov4112T45tZfVzOtG5nGeAza7ez5+fns98LCBN8Qsv+73Wy8YfwLAgHLLyyU+9PtluZ8PJ+58BnTT17OGPqY4Vub3a5qp6NL9947dmuN+xGC6sxMcY10OhmTNFq1eel0CngWcG61xkpK7/DhsTLWv+aaMj2HMB7rPtMpKmg2XkxfXCan1f7VV2cWmGaCbIy7lcc1jzSJ67XTycZreXniarVxuO66TAGw51EBnJub0DNYKuN2cut/TG9vf3uBL2wFdkI4fKFOuZ26zqpbCSa9dIUwISq7Z4uay3a7WX3ICI0wcyY6vhYWJovWhAlfuXApATIu/I5MQ3VCpLOzZcFmZYF4C3XVBWfBjuvFMZubm4xdrD0UJKANl8YkZ4RJQbe8XFx0ORMvAI5Bt5vhljPkwhjZ3BnuqCjY/W53wqgN906nWB/SGjDj8XNWF9MBzgcKl5xZF8aev5vg9cbKxsjGEunH/pudneCPuGPdJtzwMrzYGkT8yRouzD+OJ45BbFz5P/sf62W8l5cnTNz6b2VsHFCQ23gjXjjeLHS4TcbD+slrBceb/zMBvLAwoRWPtqeAnRAOvyAir69TdicuybL/7pdsX4UFpL8tVn7bloO36FnDNW3IIw6cJGYmTBy8YHDyRSaLzXGzFMpXmZLYL1yQvV6ZqXe7PuPCejwG0+9nhIpCDxcAX4uLRW0O27OFYOOXa5EFxhnDg3FlBrS05D+LfcNFbQsXmQn+5voNP7I2CmOu6gtzFAg4lnjFBEmnkwneELKxZWGMTAqVBu4P/m61ynh0OtllTAm1YsZ1aWmCF/aLx5jxwvXHzNkULwNy3RbuxdauJ8Ttu9UfowfsBwogdomCR6FAdzhvNj4o3LCvzOwNF54TVkK3CDshHE6KyEiy3f3fyH9/o86zW71E5PUicm+etfSzqbJnTDjYIjatwTQMZg4LC5PFgto4+YTd79iuZwEwVGn4vPg9DYf/s74ZUzdhgYukjrBii4Q1H178nsD0XC+oJSMTjvUZF26MGakW++ZZcDZe6Cpg14/HdNgCsPGwBY7CmOchNl6qRaHOAsD622r5fecxt2cMJ2zbo2+r2/732oj1C9s3K40tNBRONkfW3xjDZCWGGSYLYGvT6Am1b492vT6iQOW1znOHQgPXCVuZRn9GU97aZ4UO6TW2NmvAEy5badprW8IhZgngoLN/EgmLF0XKlDfABeSZ9SmoKudZK7yADX/ThFCD6XaLGlFMqKmWBZuHK7q1UINHhstarY0takWeVWCLk7Vge4YFgyeIPKbJNMDzj8qCAQptz93G9SL9oDsOn52dLQpCxJPrZJedN+YxgcRuFlaUYi4QwxHrsvHBMvPzZdcaj7u3bjz3Vko5itEqChluc2GhKCSRXjwlh5m7tYvzEXN7eX1ldx3/j0IL1wsqKefCchCRF+WfL/Ou1LNn89ruPgeXKXiMiTVpdCGxBsRaAS/yGPNgE9o0hMVFP7uHF4Qt1tlZfzEj4/UYM5vACLywcSEZxPylaFkgcSNTw3sxNxwyLm+xYRuYHGCMBjUzxBWZtblRjCFy7MZcf4Y7anHoWoppnjz3nr8a6cITaJ6liv+jlWJzsbQ00dxnZ4tt29ibwsC4o2uM/2NNmoVDCOU4GcfykNZi1m+v58emEGLrkGMGRo8xYesJCKQxVBBYoHtWjdURs5yQ+dvc4TpgizlmVU0J2xEOH8o/P+tcf5J69mxeO+ZWSgW7POmNrg5P+/K0UavXCz56moBXBk1vZNjMXAxPZnxeAAwZrcf0bRyYgG1BsEbDl2mw1iYvRnQdoZnvaUUxocvjgYD4oqBCIcPMWjVtWcYuq5+VAO9ixo/MCv37puFWacxMA2yhoMaJ7cfokZkSWqZs+bH17I2LN+a8NlBJwLjK0tLEZWjzyUoSW5Hs0kOByTSEihjjt7CQtWPrGttHGkH3HveH3W7s5rR+8NpgNx5lQm0HGrdSCljzZsmuOplQ9GmatumZpN49ZgbonkKtFBl2r1cMNPIi8vDnRY6M3AQPat+4kG2RsEWEENNurW/YPi4i7BNaL57fNyUQ8D8UBKbRzs2V3WWq8SwZz8pBxmf48Rx5l2nIMabr+fJ5LtHV4TEffs4L5PL8I6NnhcbqwftekBwtB96rghYWrhsra8wWY2rInJnpIb6qvkBjuo+5p3jt2vN239Ywzr8naG0TKVvDOA4LC0U6MS8AKz2ogHEwHC0bjCtxTAzjgtuA7VgOP5i6Us+ezWtHUlk5O4aZDqaP8eX595kRmP8YFwprUZixYItyYaFIpLOzE9w5HRHve3jaYsP62u2JuYvtxjanoWZjQtKew1iM/YfBZfYBWxl2kzCTUS1nlmD/OQDNKYnmFjLt08uCQUYXy5LytGFOmfXcPja2do8DpFUxAY/ubF6Q7pgpI+Oan5/0n2Mm2HePKaKyhNZBKniO44pM3JtbZJ5Yp7lojEnab6YRxJeTOtgty+4cz52EuBqzt3VivzGrMOYGQ3yZkaPSgc+iuxfn3Wh8bq5oZfN4TwnbEQ4fTVxnbIf0tNe2hAMyUhQC7G4xEzC2Z8EYlV1eSh8zFSNWM23tGdY2kSEjQ2NtR7Xsmuh2i/sMjOhjmhQTNlokKAyYGXpBYU/wcl67leU8cCZ8z8fKQg5x9xiHjR/niLNVwtq3xyywXRtfFNJe360NG3sO5sayrFgAYoDXi6UwYH/QtWPji7giLceErRdTwLn1NHemMRbM9oyn2HASArZvVh4qVogvZzAxk/ZiKLHfrAiwMsHzl4oZMR/B+ILNAwoMHGOrx0uMmBIat1IKPDcJMgnanewyG09r4IUXu3DReEzaFr3HwIzAYkyJNVV0IRlxo1sLNzSZleH1ARkA+5kRV/TnshDGbCkWHCw0WRtEl5zV41kPVjaV2aPqCx40/VmYmhXi7evAeBQyrJR14I29WZnMoG2+EdjSYZ90TLP23GpIg+i3N6GK82VlPOsHs+JQW2dmbYB12Nibtc0Bf84UREFr42h0w3OOaxPbjHkFRIrKFfIDHicMJNt4GW1yyrMpbrg2PRchCs268bUpYCf2ObxPRDrwuysiv1Dn2bNx7UjMgYkjhAmhoV8QfawYWOWsJ1yk6K7xBBEuVi97yltInusEmSUu4lg2j2qcWWEGFmfzeOOBQssLNqcCu+wn9lwbMfcLC1TUspAJMLNHQNxs8WEMw9PmERfUPpl58tig0IgxJGPgnM2DfUTwaAb3rKDbD4UvWjSYReTFETjo6qVC27gxA2fN2ctEQ6ZpjNf6hdlN6BZEwcR9ZwvD23yGihS6Q1P0iQyZ1xK6S2P0PjdXjI94wsa8Cp5FwLyCXXRTwk4Ih8879+6o8+zZuLYlHAz6/XgQmX31RtRI0DHtwxgpxw08JmVlTHNEbZK1KmaaHmPgjBJm5jFXCTM5jh/w/7jIY3559Bt7fcB58NJvMeYQ6wf7k1nYcCohas7YHvaDs4e8bDZsg+fIGxtjHpjUgO7HXq98LAVmwHn7RdgV6THYmHAWqd7oZuOLgtezSJj2PAHiuUQ8YcL18bhgndwfVvhi+1zYBYt05jHfmGuK03PZ2sSLXbE8Z5jEwW3wPJ2rmMO4kMhdInIB/G6LyJfqPHs2rh3LVopp0eibRoJhd4QxUPYlYzspiW9lMHPE3BdeaixaHni/yl3Av1GDQcaDTM5jLpzmiGPErh52J2BGFD6Lfld0K3muHwbc1+DNK48Bjg0LB8815jGGlB97fj7DF7Vdb/xRgBqT8Mbac+WxBRnDL5Zaje3HyqFbFd0bjCfOHe95YOsQhS/SGjJJFPRenMObe4vxeHV6bli2ZDz3jTeueB/jAZ6S4MUXYgpZr1dWpFAo4zx5OE0JOyEcflpE/lRE/pmI/Fj+/afrPHs2rh3ZIc0BaJ40JmqPkHCxx1w5vChUfU0XNQ9mJB5BeUFWXKzmvzbtA7NnWIPBzB70/9t/7LZiFwtqP2ZGc5aQx+zZEouZ8NwWmuQeU0VhjBlemElki9XSGnlOcZxsbDjf3SCGCwvwWMLCwoKfuohj4WmmKMDQDWO/Y5qsPWuAShBuqmSt1usn7yExBQI1XpvD2J6OmJbPwtpbK0g/bKHE0lW9+Az+ju2R4HgHbxqNCbaUgE79tvH0rOUtwo4EpEXk+0Tk/xKRXxaR76373Nm4dmSHtMdUcYLQ1GVCQjeBN2FMcJyVwO4qxsHzQZsmWVczSwmWFOOIpTfi4vMsBG/8PO0QFybOBRI/a0nIpK2PuEi9q2rsMVhvddtcYHA1VbdqmamzJWnzmRIgaF3ZWLECEXPfVbmPPEGUcpHF5pP3Q9ihipzCi/EBVA68LCFUglJ+d4QY440FfE0x8OaeFblU3TYGnuBjYK8B1oVHnfNlm/68dRvL6psSdko4XCoify//viAiF9Z99kxfO+ZWwnvoGvH8jh6zYAbNbcTcPVUL34jOGI+lMnrag6cNqaru2eMTOF5mMXiLll0IsSwtztxBPLAOttRQyKDvlt1faCWZgOQ9El5ygecu8aw/ZOQxYcR14zhz/Z5f2eaJ8dmzp7wHxhY+4xkT+DYHqPFzYDsWLDbwrGOzQNgyxPnEseLsMcy88mjMcPO0bXYD4Xgjs8VsIqZzjm3hmCJz95QVjHnZmMYsAp4v70RgtoJj8+hZhxyLrDq+vgJ2wq30NhG5TUS+nP/+FhH54zrPno1r2wFpJjibWC+3OEYUmNpoG9diLiUOFMfwib3rwbRaZiB4D60Z9Iva5Z3D7y1Y85sjMzA3APbFsxx4saHQw9RCzm5h7Z7H3cYONUvPTYGLDNtC1wSPoRdM9oQN99fqYxeD587AeWZaQhq0frAGj8LWS3lFCxaD3myBepouWywGrHWzBjw352dNTXPZOjM8PKYbE5q4sdAUhJSbCuvy6NdzZaLFYIKS1wweE+Ph74EXY+KMJ4w5eS60LcKOvOwnf6/C5+HeF+s8ezaubQsHJrhY3CGmbdukeZYEB7u8Re4JEtV09hRbDTGBYu2ZVYD32W3D1kVKcCCxsqWF2htbSLjAbBF7ATuvX8icUENFa4LHnxkpWyXsl2brTLXM9PiIdk+746Aqj7UX70GGjfd7Pd9lZsyexw4ZHQtrPkqFAfuKLkWm25jignPO64THn+eM3WexwL63r8eYqeHPQjhm0WOqLzN+/MRU1ypLku/ViQtg6rSVxyC71bW4WB77c+lWEpFb88/P55+zInJXnWfPxrXjloM3KVUX7pQ1ImFfJGoUXrwh5v+2K3b8Mu9DYEZm+HlxjJhLxGsbGRbmmHtHbmA8wks5rQrSxdxW5sZJ5aXbAvasDh5/1Np6vaJWbbSBmigKbNOsPY2ZBVAqZdhL3cXfnqbuMf9YzCU2Rh6z9cbbcMKysbXBjM3qQDrw3DnM7G0NxLKssC68j8fTxNa4p0jF9hygNYp0z3FJ3PjGuHrMO2aVWL9S6wPX/nngVvqAiPxrEfkrEfn7IvLbIvKLdZ49G9eO7HMwqCJ8zidPMTgrw5kVzDBNe+Egni12e7csgncgn9cuM0dvIRqRpgK6liKIdaDWZgwkphV7GntK+BozNJz27CmOd8q0R+2Lff48BvisZWVh6ipq4dyWufc4Ayrm1vAuZnAeTXkJCcYYuKy1y4wf/ecePsbArGzKPRijFRsPpC0WYCgMY6nQTK/I+Hk9ekoUM2QcX3Zh4jygIIhZeqjk2EZAxMfwX16eKIsWS0FLH4UPK12xPtkYIK7bsBpUdUeEw0wed7hRRH4r/x7qPHs2rh0VDrGJYcbnaS0pRucF8yyIZ212u0XNlU1dFBAefsiIuXxV/rbXbw4Oc8A9FqzE2EBsPFCD4/gHCjVk6IgjpsDG2mFGgBaFgXcirkEq2G5jjvNnCxVjTxzIDKGs3WK/zI3iMS9WJLxUWKYpZPp22iu/SwHPB0PB6DEzkTJ+eDgd0wZalWh5oXvGS+tOJW+wkoBCHd9ZYf/VUeLQ2rH+oNXnpbCite/xiBjtsnKEMUibV9w8iPt3bIximXxTwk5lKz1DRJ5Rt/zZvHZUOCDDQ7eI58uMHYnBv6t2XtuCQCZn/kdehGj+ItPyNHlkdGwlIEEhsfHrBzmegFpprH329Xvjg+OI/cY0Q44VsNVlEBNuuPC43ZirJiZ00OeNlgnSjDF7ro/rZbyNIaDrwmP8yCjZ3eExKGYgjBePI7vUeF3EGKsX2EdawJgOjgdr9HyGEsezUKlIxe5YqFu8w9tT4sVnbH15cTruE44zWztoGfX7xSw7Twh47iVvnSIt87hPCds5lTWIyM+JyKMisiIix0TkERF5T+q5s33tSMyBJxg1BMzLRu2ICc0IgXOXY8FsIxIkIA7+pS7W5O151OSwj71e0TTnhYyMkxdaFU7oQmFLwmNUnpluC4UtK44L8NxZfegXR8bFfcI+41x4GVCxZ4wWUFAbbl69HjNIjbv9ZvpCOuM0TtZqq+aO2/VenmSAffeC0azBelap1R3zzXtX7LgI7DsLJxb6tg6QPmNj0mqls67YXcmJGUaL+HImdpXxXHp4Y/qsZT+yixTbPNuWg4j8lIj8kYg8D+49X0T+QER+KvVsRb2/lMcv7srjF538/qKIrOXZUV8Qkd+oU9+2hYMXhEMmaowNFxUvXPyP/dhe0NHqxfs2wZgxgzhhvZbBwUTZ65XPv4+5o5Axe3slWLs2H6u3sNBX6qUFIrBJHYsdWJsxJoT1cFaYJ0xYU8djDHAxsiaKmSRsteFi5/lDHJjGYpphSoAxA8bxRlrAdE7rG+JsG9a4HnZ/xtySrLnbszE6QosnxpTRImUN3+pA4Z0KRDPdspXiWfec7ssCgoUT0yXTeGy9s5IUE1SGE+LL1uU2BUPGDrYuHD4vIk937j/DO4yv7iUirxWR2fz7+0Xk/fn3RRG5e9r6ti0cqvKzWUszzQDLICNGQvRcSrj42I2j6rtnVMuuHHuGmQgTGxIu+y/ROlItB7v4N2vPvADRckEmGNMkvfGJ7RplBo6aseejxn7iuKbGjn3svV55n4IxD/Rr8wJlv7Dn/mFAZs/P2fh6MQF01fC4GG6coopCgI8C8dyS3IaNIz+L9MaWtzenqEB4ZfDlVjgnMTci4ow0l0oKQMsDBSTSI3sNMMDPFnDs7CQvQ9Drh+feSq2jLbqUVHVbwiHKqLfCxCP1/ICIfDz/fm6EA2uU/Lo/j0l6mlOMGFst1d27JxOPDAUXmufDRk0QBYn5vdlfzb5bxM366jEt84kjY8FnMMiHrhI0h63v7A9WLRJyzDLhxeH5e9GPywvUc3eg/xzdChic9fDFOUBGYQzCCxLGtDicY9w8yM8x3eBYsebLAhyFAtIuChqj7Zj2bowRtWJ2o+BvdAGii41TfDH4zIwcLYKYRRpbp14fLPiL84D+fqYztnI8PBEvz8qLWb2ekobJCrgePWHFbihPATmHlkP0WO7Uf9NcIvI/RGQ5/74oIqdyi2W/iHx34rmrReSAiBy45JJLtjw4JcBFyJo7MgTPv86ujVjaHF5eZhFqMcwsUNPwNKjY60xtIWB5dMmwxWNBPC+NEjUVTOPE+tnVgS4HXlg8Ft5JqBZTYMZZlY3i+XxZuHOg23NLeeOGDCOlxXmMBONN2B/GxesjHhLITDCmUXt7bfgywccaPOKWOm0Xx9dbS7w3wPBiep2f9w/YMyHCsRRPsOA81HGbeevf2uD9EFwH0wMqbvgfCxFsn/FH6xD7uU1hwLAd4TAUkW8410kR2ax49n+JyN3O9f1Q5mfzmEPIf18gInvy7y8XkYdE5KJUO6o7fLYSxwA8a8DL5kCzmrXrbje+hwAXLVsOXCcyWtyNnBI+/J4AthyQ0D1NHdPwvFRQFjBeFgwzTtRCmekyw0TgufGyXZgBYfu4uHgxYi4/L2DrV+yVr7hrO2ZdpdwIzLhtXDxmjrvY0frAcfeO+7DxTGUdWRnUWOfnfU3WLozNxRQKzjzyGCdeKEg8oeYxaGO8dmFQPaZhe/e5j7iXgdNvDRfvZARsA+kZxwutMI8/sILkpf5uE3YklXWnLxH5ERH5cxFZSJS5SUSuqKprR47sNsK2xWfuAvavhjDxQeLEpRYbMlckOmRIyAw9ywEZZZUPFZkCBrVsgdo927RleOERB0yssZ2YaE3FXEe8CNF/zdoRMwUEdrchfmjBecFFDpR6zNNzXeFYpyyUqriVp02zUGI3nQUwY6nQnsCxfmLfcPyxLhZ8Nma8NwTply0Hjw4x64+tGxt/7K+3F4UZNGd6efVx0gj+79GAR2voCvOynhj3WHZjrG0cL96Q6I0lKqKo1GAb24DzTjhIdvz3PbxvIg90t/LvzxeRIyLytKr6tiUcYpp3LEDL6WQpwYAL0LRKfB+tl9aJJqz5u9m8riMYkGl5xOwRJG/qihEpAgoyXASetm4QW6ymQZkpnzKhGUfOWrJ2vbOE8FlkhpzdxQwEs3/Q1ePNB2cI2fyxpWT9RsbOKbnoplxcLAdH2RLCtvkdA55VxHTBFpKXrsvuPoyDMP3hnHt1eacNo9BAnKvcsOyWUy32zXtfuvUx5tL04kRsfXtKEOKGyoHNLR9miRlX7Lb2FKsdgPNROPx17jIqpKyKyBtF5EsicqeI3CEi/7BOfTvyPgePqHEBxI4qsDJeMBSJKaV9oiDi+zGrgl0w3oUaES4gLMN+0FimBTIZ9qGjCwwXQatVfA83Puf56D2/PruWbGEibpitZIvbcwHhvgkUuDHhgGPNix9x5dRlj2lj1tD8vB+4xUwgb/wxzsJjGguK8vhxHACtJewP01XK3+25m5j+vDJmuWJ5nuuYKwyFFAfCmbZia9zbb4HzxTRu9XsWiicYPNzwPy/xAunU60vKIt8CnHfCYaevbccccIF42m+v5zPhmD/Q3E8mPGKLh4UIMwZMX7Qdwmj+G/68D8MuJtQYUdWxRpA5oJuB0+5QYCLhI+PDwKJnNaGfu47Pfm6uzDhZyOCYoOvAEzZWjpkCZzCZBsg44pwxjogTWnU2JsvL8diG0aInTHq9olCPxTwQP45pWZ9YOfBcOUhTKFxDmOwHMX+955f35getIQYWIh49WCYW72CPMWKk5xh+aH1ahhHTj6rvUmJaRQGRUvLQKsB+s9XoKVhTQiMcqgC1f1zwKPVjmplNKjIo9gEjMXnPI1PhoFrKbYWv8cRYAgsuzw8fc/UwYzIfsBG/me4cvMbfnh8dhZKNR7dbdgdZfYYTLmr0kXvjgTEWvIfl0ZphBsFMgRcw1smWTUpI2zMs6NEthb5rHAvb84Eb11jQI43h+KmmLTS2HDy3o/eGN2+Mu914bATrtzRUFl6c6YTrEucdmb6tF9sQ6jH+TqcsND3FxLNObG0yDdpawvWHfIDXFNbn9Y3pB/FKrf/YCbRTQCMc6kJMe8YsnVQZAyQ0dj2hZhXL6kECr8ou4cVg+OEC5qBXyg9qjNPe/cwCgH21nu+UF7mqb7nEmLyNGY81ut+QwWBfvR3V/X75OHX7jlqx5/v25tt7c5h3WN7sbFEo47jzXGNfMEGgToaK1VW1aQrp0+jQtH3cr8B0zTQUo//U/7jpL8UYY/56dvHxemUGH1Os0P3JJ6pyWWO6lmG0uFjeexDjA9hH3hPDEMPbnsfUZW9stxGcboRDXYi9KAeJwXzaJrU5u4kZCsYuPO1ItahRspbCRIuuJkyLwxx1JlrUnPbsKbZX94yb2GVuLsOZLQYvEG6ZTyg0OdjKuez2HesyhsLCwVuMnG+O9XAdyDBxPKtiD3bx2VKe5st0Yq6Q1FjH3AdoPcQydpDWvLpjgVqcV48Z8mtIvWdtzg3QckEmjcKfaQgtA4sH8OkC2LeY5WBj5CkmPN82fuz6su9Gr55ASQm2VNYdupmMjpDeY8kPW4RGONQFXujeuSeYo43/G2F7fkabXNxqj+DldNuCQAKoihOge8LS5DxTP8aEFhbib4OzuquEiS2oKt8/MlE+tgHx995b4blAMLbh7ULF9F1boCYAUmmoHEhHRswCjhkTxha8tElmDDzf7KpLBR6R4XqZYIyTXZ5QYh88Z+kgTcXSkXndYDkT/qwAWDs4fxxAZ7ckWmbMxL29PFWuGi9OgzgsLMSD9Zjqjvjj2LDbkMeeFUSs3+iNcd5G5lIjHOoA++xxcphpeURlmj5rZrZYvUlkzQbdAZ6Q8TQh07RjgWEvVuJpILOzPuHxWf+8ixk3nqGbJRVUjV22oFO+e2ZWGMi1BYRlVcuBYGReHE+Yny9v7IoFxJHZsEWBjI/7gtaRx9CQqWP/2PJgBsYpmGwheYFhHjPTQtkHbs96G748IeSlYrJVY/e43zH3ifUbX4KF9XkHVuK8GlS5ar2YGaagx2JhXj3YZ+sju50xboD0ha5GHB8Tnhjn3CI0wqEKYr5lY9ZICEi4rGV7vmRP8/QWC7Zv/3OAzdN4eGFinSgEWKOMpQcys0ZNyK5YULHfjx9I6LnsmAlg9ocJNWvb4jSm6XuvuPSEKtaFQtV7p3HqsmOTMRMN+4R+X5w7Y55V9eMcVNEn7rkwQZlKh+VYETM2K49Cq46mjdowWwUew0LLyHO/eFaYXexiQXw5psP7KXitpCwHj8FjX61+PvnYU2a8M86YD6CrmS0w3FhqOKD1HvNETAGNcEgBM2FmjqhtInHzO3Jt0lLmrGpRi0UtbM+eckaRV6dZDpYtgS4oLINmO/vK2ZrAgDsGK2MLxrNEYsFzdkkg00yNOQqZWOoiM1RPQ4+51kxDTmUXdbt+9pYJC6+fWCcyLbNGjKF5FmgsZTTFzExood8+puwwozWasD7GgvaWlcP1WNzNE7Ax95aXMGEJAOgSxbRtfhshZgt5VisKND5A0/DgOYnRAK691Jh6ackx2kKNH61DpinLtGNBx+61WCyqBjTCIQU4gUaUODHInDjXGheJx6xU/WAUTjI+h4Ewew5dD0j4SKxIOKzBxlwiSIC4eFmwxBi45ztO7QVBXFPvqmahxf1C7ZItB05ZZJx5ganGhZTHEFGgx7RSrBOVC2YkquUcdgScZ3QVWnmzfBAfFvCIt8ckPZpA9yTewzH35ofriAXGrV9eLK8qnoXC1rNoOJ6DDJjT0tGVg7EUs06RXlC58oLZnE6K8++5phlfbw7Ya+Dtq0FBtEXroREOKcCJxswi08C9N3fZQrX/mInhRCFjRC3WmDL7SU0omVZkC4Y3y3CmCzIN3u9g4Gk/iCsSaUzzwYXCedkcfPOYsscoefEzM7N6WNgxA2A88beVS+1dQabGabw219ZXtBzZT29CM5b3j4KOs6h4LhgPsxo9IYUb7ZimPbcPuzbn5spCicfUFKH5eV/I43x4eftIix4tpKwkc6PgeoiNLa4/FEYpoR7jC4gb+/1ZQHhWXGw87TkWujgOLEDQFYX0ukXroREOKUCGgZtKvNhBLA3ONETelBZz46BwYBy8BWnfUbCkNhxxOynfMZrNrC2yewkFHwdmeUy9xYjxBy+4jZpcFePAOWJNa3a2POYc6PcEFsYRFheLcxlbrDGXF+8+5gwcz08dQnEMPYHrZUClLhaKGKvgMUAfdiqDC8cUf2PShEcr7Nbh9OYYA+U2UEHyYlmeGyv2Aim2SL14HCt17OplhQ1dTOxa9ix/jiexlWI42Rrh1O1txB0a4VAF3iKN+W3Zt+65i6wcuqiY8dqzVb5Mj5i9hZvSzpFJdLu+xpdyESCTareL2SKsUbNWj0wa/fIoQGLBXG/BxBhFrK9ouaS0eWaYzKg5HRHHEuNExhixDyhMOEmg1ytmq8UWOrvT2m3fivHo2LMmmDF7Co9XnyVhmPDEuAAzfBMQLCysDNIhZ/bEcLDxwTllYcZp5VWxAixbhzYQZ3TtsJLDiojnGuV07aUlP1iNApa9HduARjhUAS+M2Vk/7c8mnoN17EO1g9UwQ4SZwjQCwcOB3Q2YOVVHu2VGGNPY+OhyT3BULQBk5p7WaEwslmJoPuBOp7iQjFEww8e9JSmfL5ZHDdsYkwH2EYWB1wfUPm1cURP0/NidTpFWOIEg5vpDi81TGLBdjjkYfrFx6XbLGq7RjRcMRbrCdeS5zVBpYgacikcZPbZaE4sWmSP6+tmFxto2KzWedS+SrStbc1XWFCoPqmUXqQkQnGevDhSenJ7MLrstWgwGjXCoAvbHpoKUMQ0DzXZkSGgapvyPMQbtma9WFk/3ZKZg/eJ9AJ6J62nlzIC9MTEhysxD1R8n3EDELqter7xYvCAljq3nzkHcbPxY28Jrbq4cM8F02Vg2jsdYeacvChUve0w1HvjHfnMigheEjMWTUONEppJy3/AOdi8DzmjQo+2UQObMOnStxFyfPKd8zwS2uYW991l7fUdXJdMIv9OElcWqYDzGcyyjjNepN2e2prCNWDxhB4REIxymAdYy8DIC9tLOPKHCi8wLZiGT8SwKJAzPF2qfaO5jX9ivyYSGm5Uw9oIEjAsWXQixFFEMQsZcHnzZe53t9549vmaJCxtdR16dyED7/WJuemwcYwue3QreRits19NEOaipOqGZlMZsY8rMDoVELJbCcbBUzIVPmOVYAdMUuwv5P9PQvbbm58uCx9ZKjPGmxpzny9tvwYJnYaF41AquUy/miL9NEKGCMz9ftlhxfLBOfvc4WqxeokdVQL8JSJ8F4cBEsHt3cUJTmjYyelskmIpqpjo+w8EnDtYic0OC43dFcxzD04Q8Qotpd+iSQg0n5crwMi6QSZqZ7AUHq2Im5pJAYcauNTbbecFgeQsq2mUMnYODjKPn47YL3QDIJFPaHbqfUntLWLAgHZmA8JQWdP2hmwjHD5kTPsc+dIM6Pn6PScYuFnwx96KNP88d7vvBsow7M+JYjMnKIh3ELFi8b8oj9hmtRBt7THxBi4vjDPYfHopo95E2G8vhDAuHmM8RGYxH/HzZBHoMmokLs3NSFgP/752Hw0dqmADB+nDR2YKOpddhm0ioSIR2r8pPnBozb0MYWy+IEwcZPe085XrB/nu+c2MwaN6n5hdx997cxoyVga0ya4P7xRlw/C5irCuWnsqWA/rdYzTK88907e0twPGvyj6yOfMEH1tr6Pf3ArGewtZqlcugZo5zx31gVxriipo+HquuWqQbzHrzBDf3L6b8MA0i3tuIOzTCoQ54E8GBQdQIO52y5uERbcyXj4Hb2KLB3H5vQdl3iz0gQ2Bth7Vi87V7rhLvDP/YZpsYY/H64sVaMBMDF43nLvKYlCfUkXl584HzZnGQpaWy5p5yiVm92H9LaeYAogccf0J3HfqrkRnEEhlYM/UOhPQSBFJ14lwgU+YkBtRaWTO3Z7l+xBED9khrsRcOxfrF1owxXrYccJ2ymwe1dBaMxuDZsmeFJxW/q3PFlB7836P3LUIjHOqAMRk8SdELLBnh49uuTAvBA7GQkDw3imXDMGPs9YqTb0yZ8+ZTBIgaqBGbFx9AQuRFjQsjlofNuHL7VRcGM71n0CVn7hlmRMi0MJPI+mLnY2Hw1GO8KQbJQU5kTjh/Vhfi7Qk0xMH6hWNr9Ihtea4fm1cD739jjp51F0tX9Y4GiaWfInPylBibD1aucJ1UBaKxDvvNAW2PUaOLldtgwcXzzxYiM32bK0+B4fmbdm2gYmd4oFBCZSKmtNWE8044iMjPicgReIf06+G/n8nfMX1QRL63Tn07HpBGDaIqswMJJcVgODCHu2pNYzWNnV8v6bWlWtQUY4QYO4LbNHnUdA0/szx4QTKzsn5VEXvs8vK8sY1+f4KjzQELLDP7OVgd83kzk8UNjLGNUgsL5XmYny/fY1cVxx+QkbGWy8wFhR4zJlNMmMF6TMjoxZt/j1EZbfEY8RgbczJhktKY2UXq/c97h3jXMKaTxuJkKSvWLH3OQPIyp0QmKazoKvQsDayLY4SMKya02JhiBhsKIlNWTblhembl5MkSkM6Fwzud+y8WkTtF5AIReZ6IfFlEWlX17ZhwQOJC89PzPXNAGDMf6lxeKiZOvMeU8Rk0pz0NJnbt2eO7OzxNlmMNOA7oUsCAep3AckqIWt2saZr7xxMk3tywmwIZv3cGk+femuZCTQ6tFA7s45ylkgEwFRdpxOqqwoetF34xVJUVwnELq8ezupDxGVOz4zXYtWPrLITiew9YGHnWBLuCPPynmS+sN5ZkYGNpc8YWIaY+o0WNCkHM5czrGL0XKesN8fIymaaAJ5Jw+BkR+Rn4/Qci8sqq+s7Ia0L5aGbOMjBgzcAjvthl0p5dHbE0TWPInKVkz2LQyzSOVGDX89sbY0OCN0hZR8YYud/MFK0eZljWF94nMg2Dxs1GCJ5G6bmrvPgFMnxMhfWsCc8CwLa9zUw8nuzbR9emMds6jA/dD5yO7CUHeP55xoutAxu7mBsS59FjYjFmjLTO1jbe47EwfFL7Hry5wjVkliRvtrQx8TLPUHuPxQdZmfLqia2tlOsNPQlbgPNVODwoIneJyEdEpJvf/w8isgzlPiwiPxSp42oROSAiBy655JItD07BtYLfPU2NzWNkcPieV5tM9HsbE8YzdnBSkYDZojCzMsUYvDN7EF+vb94iRpeItc2MvQoXXuQsYMzFVmXtcAwE8ZqZ8RcTa6ueYLCxr9LAPV8uWkkmULBeZMg4h4Y7+pINP4/J4StTU0yU++XFFowR8VEfuHGT6YXHDv38Nn+GI84J9pWPXWdmFsv4iwkKK28WBD+Prh7M+OHX4cbiZKwQcRwylr2HzN7Wl7lzPZci0gSOeew9Fkx7nPb8RAtIi8j/EpG7nev7ReRZItISkRkR+UUR+Uj+zK87wuGNVW1ty3KI5QvjwuD9DexXNgLG7AvVYoYI+5+NkbEgQuGADL1qIVlQmQk9toHGFo1nsaBVEXMxsWZpjMLTbm2HKFtIqLF5i9bGE7OlkFHa2Hjj4fna2fKrIyA4EM5+ZmZ6nCWFwpAXP9fH+DMePKaYNRN7pwJaYqxlIr3gfoGUYEUmlxo7tMQwvoXMDJmqMVOjCX4VKVoEtrHPznWyLLuYILV01pi/ny/DD/9HRQrHlAUVp1xjTMsTJlVu4RCK/IoVuyeb5VBAQGRRRO7Ov599txIvAPPre9LdC4TGdoGioMGJjDEz7+U2aK4i0XvtxfLyOYDFfcAF7uXko7BgV4yXCeJlVzGzwvZQE4+NjeGBG6CQ4XnPLC4WFzdqlZxBZMzFhLVXp9EA7oPAOA0f326vXUXBzHPI82GMp90u7wbGOAy6/WwO2Lqy/no0wZo0j2MqiI0H7dkcM64cHEZByYqK5+tHV6cxco+m616osHFQ2RsHni8TmKmYGgsM66tnxfP6SglZszB5PcY2KE4J551wEJFnw/efEpFP5t+/jQLS95+VgDQHBzEzwnufAzObGMFwJpFNpGc+sjY+Pz9p32OIzEjZV+199zR3rMPLHffGB8tj0BeZP2aAYFAZGSq6VDDY6KUVcnaPjYnn+7cyHAOIudbQOsL+Y79QiCANoKuoaoFjOWaQqTgEzo2q7/I03Obny++a4AwtD7/l5YlLKnaGmIcXuptQG44lUnCfPe0d54djYoavh1sqtoBzi8/EPAdeXCAlpDBgjHPlrSW0zHk9YpaWN25Wx5PVchCR60Xki3nM4TMkLH42z1I6KCKvq1Pftt/nwPEAZMDsZ/SIGGMLmIXBxMtMx8vm8AK6SIApvyfHQJiYPU2y1YpnxrB5LzJJnUQXj5ei6WVXmNbp9Q2DbuzLt/851pHSIL1sFLYcYnEVj0ZijBFx41RJvJBhIvNgiCkcLExQaMUsBGyL3aScaYRtoyD13FveuKMvvdfzX4bjuVBM8POpsUgL1jdUksyNafTImr0F3W09oO8e3apmhbMGjn0xQDxMEcJ9ThxvYh7DypnRJq5hHBtOrPCsLMZxSjjvhMNOX9sSDqlMA3YHmTnHk8JBWvz0mDUuOJ5YfNmMtYd+W7Qm+PRQ1rCMGeO7fpm4cB8EpvTF9negi8DKYfuonWNbXs531SYebN8WN7obMNefg7ceI2Ytjvd5oOC2ulKpt2YJIf5YD5ZD2oq5AjjYjYfCIb5eDMzGG11Sth/C6mQmzYwRxwcVBNZSDQ+02LwxNNr0BAgLL7TkmBaq4m2o4PC6QtwZUNOPafj4GzefoSuQ14q1beOEeyW8mBmPOSsEMaGAZ3ltERrhkALWvm2iTRswIkffoec79cxwDrialoD18kJgjQUJOOZ2QWJLuQ+8gC/+ZnMeCdL2R1gZ3CSIjAQFIzIGdDPxhraY1o6WmgXmcAy8FFkM/nr+bbyXUgw8BYF3qSPT9bQ3ZpY43siQ2RLDuhBnL+DM9ILjg7TAmrmnGadwRyXEozGsx8PBXCSoHaPwwrlEoa0aj0XFLlzPXnag9Q1xw/Fhoeu5f1Ltm6BkPmDKBgbF+ZWiSF8GMeGIlu8WoREO0wC7HowQMCDIWhsSNE+gaaCc02zAvk0mStPOcGMQWxCcohgLlFvmRGpRsTaNhIgMyupFrb1KMLHGjhoYunY87RsXLWe98NzFmB4zRO6rNzZ4mq23bwDTcmNuKc9VxK4cGwOmFXb1MG151lCKiVl/sS9Wh3ciMLbP84m/Y5l+KQuOFSf7LzVGZp0aQzcmi9Yr189aPqbEemPk3Y8JWquXU0s5GYJplt1uHn8xqIppeZvrakIjHFLAloMn1ZG5oVaEWn1MK+csExYSMYbGi5GzFgy8IBj3D60UFHhGWLiwuD4viItME+uqOp3ViJgtIAwKsxuDhRfjhIBMwxPGPBYsLHhsbB5Z8CPjROvK5poVBxY8nC6KwsXDKcb8vQQCtqysP5hOzH7tmMJgdGvZWCwM7Du6TXh+PWZn/yMTZveSCQKeb48eY0qCWVbeYYPYN+63974P7K8JVm8ns82jF3tDBYmTJWxO8MVANm6xeJK3KXdKaIRDCjzXQsxsZPMe36XMZe2efcbiDGw54H0kGo4HcDA8pXkwxExmZHaojWB5XsysJXpj5z1jTByFViohAH3R3mLDuYwJVLY6PJeiN+8e80S3kEcnHpP3LB6vXMrSZObJgsfby8F04TH5VmsSH1tcLG+QM8ZlTDFmQWObMZde7HvsWRwjpgcua7SFZ06xAMcYgzfPbG3jOkfLw1MKUYGxseKED/zfhElMuTR3I2Yr2hx5/GRKaIRDCjxCjTFc/uTMmfl5f4u/uXnMBTONpPdcMEjI5o7xMoeq+ozmdWyxqhYXA+d+ozvCY1YWSMc6Y9aS95yXAYKL3XNX4Jzgf7HAKltF+DwG/r0MH1U/ZdbLgvGEEQa+7TenX/Lc4Rx7eHsKiD3Dey1imjgyShRWqJVjqndV3r09Z67NVFDeY/ie5eAJUs/a9OiXM5/QnRYLEtucc+zMhBHOtzcvbOkhXlyW+QgrpnjFEjpqQCMcpgWPQFWLlgX6+jG3355DLYzdKNNIen4W88p7vXI7nu81xpx5V3esPFoTaCkwM0fXXCrGgm+2YibgaU88ZljOs3C8V6Z6rg5sB11OHnPgBez56j3GEHMPcRlPi2Vt2lNSYpsNcT5je2DYjYU4susuZm16AtXrr/3GVNPUuMSsKyyPgg419lhGF1qHXko1r0srzynn7Laq46JEK5efR6XD1jHWaQIa+8Kn1m7RemiEQ13ABeAxCk9yI7NBIeGZvymXT8yPHnvWC4BzVpUnWGK58bjgWVvxNEeRsi+XtSGvHykmy4sX3XWeJWPPGzCDDaE8n547gU9xtfsWu2CrgfeDIJOwueEgNbvr+v1yzv7sbPGICGMCMSaW2tuA88u/eezZVZlyZyENoeVg6bO289tj8h6de0yVBZQXIOfssZgywoB1YmpqShlky5pjeLG1jQIi5qZFvuHta8B547HEsWgshzMgHJAYmZky47XyzCxMq0BCSDF7j7hY80o9oxoPIuK7GFDLYCZu/cAd4LiZL6URIZPiXb9otntClX2nrVZRO0dB5MU/VMv7JwxszPGcHv5E4Yifnr8Z5wR/m3Bg94KNkcdUTKBhbjozaxyLWOYUMstYFotnydrccV88ujeXETJcfjbGyNCaQeEYs8ZRSPJaQQUL5xhxQIWjikky7fI+B6/PWA4VLptrVGY8AcH7lrCMPYtrJYZ3bOxSz9WARjikwDNZ2WTzNGHV8oYW+4xp+ezO8NwJs7PZMygQUFNm8znGICyTAjNN0DxFhmB4WeAMhYNHgHgv5j5AbQcXEAoPFmKsISEjYN89Myp2MyBTZ4GEGhjWxQHY2MWZNR4TZo3Qc5mxFYWBXvtEFxkzcWuXjz5nfLCv3W457TLlMoq5STiN1rMM0Wr24lsGiDviipq9rQMWBNYnL04Tc1nZs6h84DpkHFP0jjSHyiFboR6f4NhmTJCmhMA2BUQjHFLgafKeBhUzG1N+aiNUdutgEM8C2ayRIbHyTlTG3zNR2d/JxM2ZQB4TZDPWG7uY4EAtFxdNlRmOQVoWvtw2p4cy/pi5g0w/liaIrjNk0p4gUq32y2OWSsyVxdYNM6Yq9yQzIMbHc6cZTccCmaxFe0zOlAhWcoy+TZNmSxLp07OeWFHgDDh2a2EmHNMTuzhTTBf7PI0lwvjj80yTnIFn7WE5PH+MLSNctyYwU8kLNaARDlWARFSVJulNAhMdZ1RwWiALDbxQSNkCZq03xZQ5WInmOpvVxvg8xuoFRLE9ZBiey8n7xPpiQUB2sVk5O2XVGAOeCBpLPUaGi3OMcxmbb+x3LBOHaYI1e9xjYIFNY5Q8x7z72CwLfhUnA9IQbphEBst04zGbWFZVjE65Xo/+er3iDmDVsmWBdaFl7dElMkKkE3TTYaoqWgZ1NWy2RDyBzvPf6fgp0qj0eTu0UWDjZfOHbicWLLG1OiU0wiEFqL2h+8GLMXj+SLYUmEHi4ogJDXxREDPudrt4vpJnmaimhReXQS3U02xwwcdMdiRSjLFYP/BoZ9MmzSJANxePpwWPvXiBx1QMZ9becHHOzxcFi5eCmUpj9AQc45z69AKKLCDZ3YYMMaUdIk2xS8yUHbYQkHY9ywS1dhNQyMAxQ47XEisFHL/A9cS7/7kO3qzHmr63Yx3pAoW6FwvjscA1wkoZXkbTyBM84eNlU3m04yWWYBmrp9WajP2ePWmloSY0wiEFnqaCTACJDifYpDq//SuW4meBYqsfNVbvLCfWzhA4i8MYg+EQ05I8s9rcOF62By84zAm3TxZsvKhRA/I0JMOD/ei24I2BzM8Xz03yXBfMSDiFEzXOmFbqaYC4gD0GXSWYPXcQji3SDWq9sVRIby5Vi2cWpVx37DLyGCgLYKSZWAYOCvKYkI2tK2b8SDcxzRjH3QTYzIw/l0yDbF2h9eu5E2MXP4fj5AkQdl3ZM52Of0oz0wILEs/inQIa4ZACDuwis8YJNfeFTRwzHlvY7MP2YgC2gL1J9szMmDvDixkgA0xZESn/K5bxXFt2IQNGhsL7J1ijNYHCgd1Wq5x+ia6YWPaSJ3yYuTPjYZy8Ywu88Yox3KqxVC0GINlyiI0/+tQ9GsAYRSxzDMcrpuzgOJqi4gknfB5jNDgHXhpzzCK3Ojzh6L3HHAEFG1tAxphR8UI3F3oLUi5DW3+4rjFAbtZJ7AgRz7rnebPxYAUE72GQmwVEap1XQCMcUsDpZDFT0CM+XkTGwHifgGk2MU07xtBi2iqW4ZNCYyYuPseuKS9+wGD4YGYMuy3YDcVZIannYq6sWNyA+8XtMhOPMXArx3GF2Bh6WjDX5bmRVItznponr04uw0qNZ2168RR2x1TRhxeM99w5rBBwvSjkPCuBYy0pemRtGteL9YkVKM8ywDFiZs1rAk93RfBcnCz40DKwOnj+PbqKxcM8q2SL0AiHFLB2ZcyeJxCZB2sFVoZT43DSWNPDABpqSN7CSTGPOsFqxoGZJxJ4jBl5baMv1MqzRs7ChzXQlJsmFgMxIexZOXgWDloEMYZtwDn1Mfw8bQ7HAa0ensMqIYxMJbYxisfI2uJxNyblKRjTWJao/LClh/TmMVmP8fPcYz9icQrEx7OuWcgx/qk1FMtQYmaNZ6UhmFBCy4HXOlu2NmfsWmJ80TuxDQGQgkY4pICZKWpMqK3iQuJnOAaQYuZVk4yLhYOJnlXBVoOngRiOsVx/xD/lW/f6xAwxZgnFtO1pxyrm1mAGzZcFpy2Qxxouv1inDn68a9oT8sZU8J3RsXpZ+08xBq6HGRArLixIY/c8ZhpTmLzncb14/eG08JgyxEqZMV8UntOkcXo4s5BDBYHHAS3Zqvrtu+FssQR8i2IdJRKtI6/sDkAjHKYB1qQ93yBPEmudddqoYjy4sEz7Re0mxgg4GMw44p4HxgUFo6ede7ibxs7+aXzLW5U1YnWk/MtYFhky44ynYCI+nsAwnLyMrGmFAwtVFpQeU0kx+IWFarciMiLOkKqaNwYcR7ZYsC+mNKXmki0Ofq+Ix3wZWLAY3aKFUpdmsA8pJS9maamWMw3rMGovIYLb9XDwLDFPSfX6MCU0wqEKvAnCbBiPgaIGj5utMAXVXAlsVbBLIMaAPL+7tetlGXEOOS5A1owNUKtnn7lHjB7unsWA96o0u2nKqhY1L3TVxPzL2J9WKxtPy2rhY1DQGuF5Zn+2WXdoNWBqLLvrYv1GZmrHmHMGmgG7plj7RaaDrhk8GoVdNjhWePSF1eVl7sTmiQWnWT5s2VQF47E879XxmGUVk6z7f+p02Vjb3l4RK4+xAU8gxnhBSuhje57QmwLOO+EgIr8pIl/IrwdF5Av5/UURWYP/fqNOfdt2K3m+Yhz02AL0LAbOHPK0fGZaVQsNU1w5kIbPso8y1gesnwN7dj/24nXG3fBC4kfh5cVSYv2LWSveuFi7vOM2pjUzo+W9BJbZYkInNc/cFxxbL35Tx40Si3nwvCAtpRgFPu+5d5gWGVccC7YqUq6uWKqll37Na43djinLiQVcLJMtBoiPlyrqpc8anhxXQYHHgMzf608VL/CSMLy9KU9Gy0FEfllE3qMT4XD3tHXsyD4HTKP0AqGeGR/T7NhyMI0wpq1NYxZzMJTrQyaLC97TFlMWTJVGkiLqlHXhlbN7HpOItY1uMl48qfJWtzEHc0FhfAnfmqYa1wx5LKxeZupVQtraML9+zDJEhj8/Xx7/mOWAWjhv4OJYjeHFqaXM0D1gCzJm4fC4xTLzYuvDo5NYJlsMvNgUp8XWcSljnz36MIsdX/rl1Rf77gWlrc7Z2S0LBYPzVjiISBCRh0TkW/RcCQeeFI9IPe0spWFw3TF/YeqZWCpmVfYOEqsxhljgLiWc6gqumOD0GFfMmsB0UsQ91v9U4NDDi5kWauJeimVsY1TdhcjlYwufy7OGzxYMMnxv93AMFy9tmPeTpOitKnYUsz5iCpYBauLTxA+8dVo1xgxsOZjli4Kzzrr1rBgEFEIx5YzrRyvI66tnVW4Rzmfh8GpELhcOp0Tk8yKyX0S+O/Hs1SJyQEQOXHLJJdsaoAJ4hOWZcal9CAY2iZgGV2U5sPDx6mahw5oPvvQezXsv/bNKM5pmzKqYKhM1P4fjzHXVEcqxuYkJXE8DZWaK5ba5EKNg9eMufWTcKXeDZzV4dfM4etaiRw8xDTvG5KpoyGPinEZaB1LtxPCJCWabb9wzFOuLd4+ZOf6fSl+O9YGtIKRfdIHNzz8xLQcR+V8icrdzfT+U+aCI/L/h9wUisif//vLcqrioqq0dfxMcQyrglNJ4vJRTJJDYJqWU1oUL1jQW3G2Nvk27cHcxth9jHHVMe4QUw0mNFz+X6nOdzCdPWKbwrhLO2xWaqXb4f4x7xTYI8jPeWFYxoDqWi6eh8uXtZ0jhYRCjsWlSqFPQ7/ub9Ho9v20UCFUxlRj+yMzRGqpjyXiWR5WgR/reJpyXloOIzIrI10Wklyhzk4hcUVXXGRUOVRObImqbTCNWPiTLywryhAMSJAoUDk6jxu1ZFEysMVdCFbOcZoxSZVOMD/HzBJsxRwsk8yawabRQbtNL5/TKp4Lo0zBM1eq4RgxX+8QdvHU2MmJdfC4W7uzt9YoHP6YSIrjfKdyr5jw1njhmhjdv+rP5Q8brCUqkGcs483BI4ZeyhqoULWT4VYLRBIl3BtMW4XwVDt8nIvvp3jNEpJV/f76IHBGRp1XVtSPCITb5KZM5xTg9jYDBfOytVtnPGUv14w00npbCwVo0UVWnC94h86gb7PPAmJ8x8RjzRvy9LDKDmFZr7rOqvH9rKyWAU/3l9jl4WeXamIbOYm3zOKFgrMu4PW0UrVLcM+IxumkEUQoPxofHAftsbfLeCSvDG1gxBdxbJyYYYhYGK4F1lKHUWOM48gbDswznq3D4zyLyE3TvjSLyJRG5U0TuEJF/WKeuHREOvOCQMDBds0rL5/pS2oC3OcsjZIRpNG8WJgZ8r4rYp00T9HCJnSWF5TB47gWbcbzRakC3mm2+q0pNZKHH41ol2PF5PA+IaSOm6aIA5xx4T+jzOLEAxOymGO14NMn9YMsAL6/uFJ17Y1BlbePYc/YgC0N+gZNnGbBlzQLA6CK1rrhsnbXtjUHKwkAecxaFxHkpHHbyOqOWg+pkAvkseIOURhA7iKzXK76W80wQRozJxYg/Ruwp4REbN7tvC9xOv2TLgRdJyqrgHae8WFmLZqaC/cNnPKFXhwGwsKrKtuE9ELwnxupkZQQVhu345usoE3jPLAfTbGPKUx0rwMrWyTSz52Lvc8DU7DqWC9NoLLUbAd18MaFTZ616Y4D0bnOb2nga688O8IpGOGwHULOKCYcqcz3GCKushO1CXQYXO6VzmjZi6XZeSiouTBQeVQyWd82imY/aNNeTYn6x9rbDADzo9/0381nqpPm7vb02rPmacEnBDjKQQn3TxjNwrlMaeEqrNjpCiNFdHagzXzjWdcfQG6OYIubtoUjFR7D+aVLjK6ARDtsBDABX7WhGMGaAAoKFQl3N3eqbdrHXeYZTIlPPpQiWmXJdjdLbbBXDFzVPz+Uybd9jAqJKcGylrdQcoyKB2WY8H9NYDinX1FZhKzQYEwBcn6dhpwL+vON6GjqIWdRcP85JnX7XEVh115ZHL2hR1U1eqIBGOGwHUjnMVYsFF31qoW6XuWwH0JxNvYa0Dg4xEzrV3xSD8KAug6pTF86PtwirmOs0zDJVloVRVdk6NLhdLTNlDU9TV90xqmKOPEZ1cdkqLSNOKaXQo+WtZDux0pOqw9xt20kQyaERDlUQY1ZVGkYVgW3FFI89mzoQbDuAhBnT8BinlLaFOE7DTGJjnUrxTGn4O2U5pJjrtMxyq4wyBXVpsErocPlY+nCMYe2EZYL1xQLfMUWtqr6t0EedstMKqVT8C70Lqfq2kyBC0AiHKsAJjn33YKcILIUTB1pTWvlWhdG0bhTvOXR5xI5kSNUTM8VT5+jHNP+dhKpxniaRYLvabl38tlOvlamzg3paXOtAzPrE/Sdn6uU3WxF0Vc/wmqy7TyJlGe6gQG6EQxWkLIed0oq2omWmMju8+lLayXYgNibInFNHX1TViXjzWON7l706qoTaTmu2qtMpEB4uVULzTOBct95pteQzZTl4dJ5SInYCdlrQYZ1bPRqE+7/DdNEIhzMN2110W3lmJy2HKogxQ7YcphGqdYl+uwviTCz47SoQhtNOBos93HYazmTddds8k4JzK9ZznXqnFWixtXEGaLkRDmcazgQDOh/AEzY7tTjPFqM5E+1st057vuo90VsBo0V8cVGs/bPJ5M8U7FRfzuQajuG4U/e3AY1wONPwZFpsCE9Wobdd2KlxmeYYk7qAmmosRnUuGOGZgp3qy7lYw+fB+mqEQwNbg/NV6J1rvM5366nfLwfKz7DvegxniuGdRW36rMF5gHtKOITs/yc2XHHFFXrgwIFzjUYDZwsuvljk8GGRXk/koYfONTZPDPjQh0Te+16R97xH5G1ve+K108z5GYEQwu2qeoX338zZRqaBBrYN73lPxiTe855zjckTB972toypnknBcCbbaeb8rENjOTTQQAMNfJNCYzk00EADDTQwFTTCoYEGGmiggRI0wqGBBhpooIESNMKhgQYaaKCBEjTCoYEGGmiggRI0wqGBBhpooIESNMKhgQYaaKCBEjwp9jmEEB4Rka9so4qni8ijO4TOmYQGz52FJwqeIk8cXBs8dxbONJ6XquozvD+eFMJhuxBCOBDbCHI+QYPnzsITBU+RJw6uDZ47C+cSz8at1EADDTTQQAka4dBAAw000EAJGuGQQf9cI1ATGjx3Fp4oeIo8cXBt8NxZOGd4NjGHBhpooIEGStBYDg000EADDZSgEQ4NNNBAAw2U4JtaOIQQvi+EcDCE8NchhHefa3wQQggPhhC+GEL4QgjhQH7vaSGEPwoh3Jd/ds8Rbh8JITwcQrgb7kVxCyH8TD7GB0MI33uO8fy5EMKRfFy/EEJ4/XmA58UhhM+GEP4yhPClEMK/yO+fV2OawPO8GtMQwq4Qwl+EEO7M8fz5/P55NZ4VuJ77MY29P/TJfolIS0S+LCLPF5F5EblTRF58rvEC/B4UkafTvQ+IyLvz7+8WkfefI9xeLSIvE5G7q3ATkRfnY3uBiDwvH/PWOcTz50TknU7Zc4nns0XkZfn3C0Xk3hyf82pME3ieV2MqIkFEnpJ/nxORW0Xku8638azA9ZyP6Tez5fAKEflrVb1fVTdE5JMi8v3nGKcq+H4R2Zd/3yci//hcIKGqnxORY3Q7htv3i8gnVfW0qj4gIn8t2difKzxjcC7xPKqqd+TfT4rIX4rIc+U8G9MEnjE4V3iqqq7mP+fyS+U8G88KXGNw1nD9ZhYOzxURfFP5YUkT+tkGFZE/DCHcHkK4Or/3LFU9KpItVBF55jnDrgwx3M7Hcf7nIYS7creTuRbOCzxDCIsi8p2SaZDn7ZgSniLn2ZiGEFohhC+IyMMi8keqet6OZwRXkXM8pt/MwiE4986nvN6/o6ovE5HXicg7QgivPtcIbRHOt3H+oIhcJiIvFZGjIvLL+f1zjmcI4Ski8ikR+Zeq+o1UUefeWcPVwfO8G1NVHarqS0WkJyKvCCF8e6L4OR3PCK7nfEy/mYXDYRG5GH73ROSr5wiXEqjqV/PPh0XktyUzHb8eQni2iEj++fC5w7AEMdzOq3FW1a/ni3EkIh+SiUl+TvEMIcxJxnA/rqqfzm+fd2Pq4Xm+jmmO2wkRuUlEvk/Ow/FEQFzPhzH9ZhYOt4nIt4QQnhdCmBeRN4vIZ84xTiIiEkLYHUK40L6LyGtF5G7J8PuRvNiPiMh/PzcYuhDD7TMi8uYQwgUhhOeJyLeIyF+cA/xEZMwUDH5AsnEVOYd4hhCCiHxYRP5SVf89/HVejWkMz/NtTEMIzwghdPLvbRH5eyLyV3KejWcK1/NiTM9GRP58vUTk9ZJlXHxZRH72XOMDeD1fsoyEO0XkS4abiOwRkT8Wkfvyz6edI/w+IZmpuymZJvPPUriJyM/mY3xQRF53jvG8XkS+KCJ3SbbQnn0e4PkqyVwDd4nIF/Lr9efbmCbwPK/GVEQuF5HP5/jcLSLvye+fV+NZges5H9Pm+IwGGmiggQZK8M3sVmqggQYaaCACjXBooIEGGmigBI1waKCBBhpooASNcGiggQYaaKAEjXBooIEGGmigBI1waOBJASGEYX565d0hhP8BuePPCSH8Vo3nVyP3/3EI4cUVz94ZQvjElhDfIajbzwYaqAuNcGjgyQJrqvpSVf12yQ7be4dIttNcVX9oG/X+Y8lOwnQhhPA3JVtHr843LJ4T2IF+NtBAARrh0MCTEf5c8sPIQgiLIX+fQwhhIYTwX/PDzH4zhHBrCOEKeyiE8Iu5FXBLCOFZIYS/LSL/SER+KbdKLnPa+qeSbVj6w7ys1XVNCOGevK1P5veeEkL4aMje03FXCOGN+f3XhhD+PIRwRwjhxvzsInunx8/n978YQnhRfv974Jz/z4cQLqR+7oJ2Ph9CeE1+/0dDCJ8OIfx+yN5p8IEdHvcGnkTQCIcGnlQQQmiJyN8V/yiUnxSR46p6uYj8f0Xk5fDfbhG5RVW/Q0Q+JyJvU9U/y+u5NrdKvuzU+SYR+U3JdmO/Be6/W0S+M2/rJ/J7/0ZEHlPVl+T3/ySE8HQR+f+IyN/T7KDFAyLyr6CeR/P7HxSRd+b33iki79DssLbvFpE1wsmsppfkOO0LIezK/3tpjvNLRORNIYSLpYEGHGiEQwNPFmjnxx6viMjTROSPnDKvkuy9HaKqd0t2NIHBhoj8Tv79dhFZrGowhPC3ROQRVf2KZMcxvAyOVr5LRD4eQlgWkUF+7++JyK/b86p6XLIXu7xYRG7O8f8REbkUmrFD+BCnm0Xk34cQrhGRjqoOpAivksyaEVX9KxH5ioi8MP/vj1X1MVVdF5F7qK0GGhhDIxwaeLLAWq5JXyrZm/3e4ZTxjjs22NTJWTJDEZmt0eZbRORFIYQHJTvr5iIReWP+3xskEwQvF5HbQwizeft8Xk2Q7Az/l+bXi1X1n8H/pxknVf13IrJXRNoicou5m2r28zR8r9vPBr4JoREODTypQFUfE5FrROSd+fHSCH8qIv9ERCTPQHpJjSpPSvZKzAKEEGZE5IdF5HJVXVTVRcne0vWW/L+LVfWzIvLTItIRkadIFpf451BHV0RuEZG/E0J4QX5vIYTwQklACOEyVf2iqr5fMjcUC4fPichb87IvFJFLJDukrYEGakMjHBp40oGqfl6yE23fTH/9RxF5RgjhLhF5l2Sun8cqqvukiFybB3YxIP1qETmiqkfg3uckcxE9V0Q+FkL4omQnbv7fmp3V/wsi0s3Tbe8Ukdeo6iMi8qMi8okcr1ukzOwZ/iXUsSYiv+f0s5W3/5si8qOqeporaaCBFDSnsjbwTQN5sHpOVddzRv/HIvJCzd4h3kADDQA0/sYGvplgQUQ+m7ubgoi8vREMDTTgQ2M5NNBAAw00UIIm5tBAAw000EAJGuHQQAMNNNBACRrh0EADDTTQQAka4dBAAw000EAJGuHQQAMNNNBACf7/xZ3V/85aOJcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# superimpose the RA vs Declination scatter plots for both Antares & IC data\n", "ra_antares = []\n", "decl_antares = []\n", "i=0\n", "while i < len(dataset_antares): \n", " for row in dataset_antares:\n", " ra_antares.append(dataset_antares['ra', i ])\n", " decl_antares.append(dataset_antares['decl', i ])\n", " i=i+1\n", "\n", "ra_icecube = []\n", "decl_icecube = []\n", "i=0\n", "while i < len(dataset_icecube): \n", " for row in dataset_icecube:\n", " ra_icecube.append(dataset_icecube['nualpha', i ])\n", " decl_icecube.append(dataset_icecube['nudelta', i ])\n", " i=i+1\n", "\n", "plt.scatter( ra_antares , decl_antares , s=3 , label=\"Antares\" , c=\"red\" )\n", "plt.scatter( ra_icecube , decl_icecube , s=3 , label=\"IC-40\" , c=\"cyan\" )\n", "plt.xlabel(\"Right Ascension\")\n", "plt.ylabel(\"Declination\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 3, "id": "47d6b5a0", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAbJklEQVR4nO3dfZBcdb3n8feHkJtESLw8DFSciekBopJYMkmmckOBV9irlxCuBlfBAUviU8WlcMvHqg3XqhUtKcRFoPAht7DABEuSGxWKFIqSDXFVCggTNpgnAoOZMUNiMo5eSYpNlgzf/aN/k+1Mumd6Znq6e3I+r6pTffrbv9Pn22d6vn36d875tSICMzPLhlNqnYCZmVWPi76ZWYa46JuZZYiLvplZhrjom5llyKm1TmAoZ599duRyuVqnYWY2rmzevPnPEdEwMF73RT+Xy9He3l7rNMzMxhVJXcXi7t4xM8sQF30zswxx0Tczy5C679M3s2x6/fXX6e7u5vDhw7VOpa5NnjyZpqYmJk6cWFZ7F30zq0vd3d1MnTqVXC6HpFqnU5cigt7eXrq7u2lubi5rGXfvmFldOnz4MGeddZYL/iAkcdZZZw3r25CLvpnVLRf8oQ13G7nom5lliIu+mY0PuRxIlZvKvNL/4YcfRhIvvPDCkG3vvvtuXnvttdG9zjHmoj8KOUBFplztUjI7eXV1QUTlpq6iF6yeYPXq1Vx66aWsWbNmyLYjKfp9fX3Daj9aQxZ9SZMlbZL0vKTtkr6W4rdIekXSljQtLljmZkkdknZJuqIgPl/S1vTYPRrnHXZdQBSZynsrmVm9O3ToEE8++ST33XffsaL/61//mssuu4wPf/jDvOMd7+CjH/0oEcE999zD3r17ufzyy7n88ssBuPHGG2ltbWXOnDl89atfPfa8uVyOr3/961x66aX85Cc/4fHHH+fiiy9m3rx5XHPNNRw6dAiA5cuXM3v2bN71rnfx5S9/uTIvKiIGncjvvJ6e5icCzwALgVuALxdpPxt4HpgENAMvAxPSY5uAi9NzPgZcOdT658+fH/WKYcbNrHw7duw4PkCF/7PKeL4f/ehH8clPfjIiIi6++OLYvHlzbNy4MaZNmxZ79uyJvr6+WLhwYfz2t7+NiIiZM2dGT0/PseV7e3sjIuLo0aPxnve8J55//vlj7W6//faIiOjp6Yl3v/vdcejQoYiI+OY3vxlf+9rXore3N972trfFG2+8ERERf/3rX0vmecK2igigPYrU1CH39NPyh9LdiWka7Id1lwBrIuJIROwGOoAFkqYD0yLiqZTQA8DVZX0ymZnVwOrVq2lrawOgra2N1atXA7BgwQKampo45ZRTaGlpobOzs+jya9euZd68ecydO5ft27ezY8eOY4995CMfAeDpp59mx44dXHLJJbS0tLBq1Sq6urqYNm0akydP5tOf/jQPPfQQb3rTmyrymsq6OEvSBGAzcAHwvYh4RtKVwGcl3QC0A1+KiL8CjcDTBYt3p9jraX5gvNj6lgHLAN761rcO6wWZmVVCb28vTzzxBNu2bUMSfX19SGLx4sVMmjTpWLsJEyZw9OjRE5bfvXs3d9xxB88++yxnnHEGH//4x487n/60004D8r0t73vf+459oBTatGkTGzZsYM2aNXz3u9/liSeeGPXrKutAbkT0RUQL0ER+r/2dwArgfKAF2Ad8OzUv1k8fg8SLre/eiGiNiNaGhhOGgzYzG3M//elPueGGG+jq6qKzs5M9e/bQ3NzM7373u5LLTJ06lYMHDwLw6quvctppp/HmN7+Z/fv389hjjxVdZuHChTz55JN0dHQA8Nprr/Hiiy9y6NAh/va3v7F48WLuvvtutmzZUpHXNaxhGCLiPyT9GlgUEXf0xyX9AHg03e0GZhQs1gTsTfGmInEzs6HNnJk/1bKSzzeI1atXs3z58uNiH/rQh1ixYgXnn39+0WWWLVvGlVdeyfTp09m4cSNz585lzpw5nHfeeVxyySVFl2loaGDlypVcd911HDlyBIBvfOMbTJ06lSVLlnD48GEigrvuumsEL/JEynevD9JAagBeTwV/CvA4cDuwOSL2pTZfAP4hItokzQEeBBYAbwE2ALMiok/Ss8B/JX8w+BfAdyLiF4Otv7W1Ner1R1RE8a8qpeJmVr6dO3dy4YUX1jqNcaHYtpK0OSJaB7YtZ09/OrAq9eufAqyNiEcl/UhSC/n61gl8BiAitktaC+wAjgI3RUT/iag3AiuBKeTP3in+fcfMzMbEkEU/In4PzC0S/9ggy9wK3Fok3g68c5g5mplZhfiKXDOzDHHRNzPLEBd9M7MMcdE3M8sQF30zGxdyFB/VdqRTrox1nn766cfmX3zxRRYvXswFF1zAhRdeyLXXXsv+/ftLLnvHHXcgiT//+c/HYrfddhsXXHABb3/72/nVr35V3guvMP9GrpmNC/2j2lbKcC7zOnz4MFdddRV33nkn73//+wHYuHEjPT09nHvuuSe037NnD+vXrz9uGJkdO3awZs0atm/fzt69e3nve9/Liy++yIQJE0b7UobFe/pmZkN48MEHufjii48VfIDLL7+cd76z+BnoX/jCF/jWt7513E8ZPvLII7S1tTFp0iSam5u54IIL2LRp05jnPpCLvpnZELZt28b8+fPLartu3ToaGxu56KKLjou/8sorzJjx/0eoaWpq4pVXXqlonuVw946ZWYW89tpr3HrrrTz++OMnPFZsyJta/I6U9/TNzIYwZ84cNm/eXPSxT3ziE7S0tLB48WJefvlldu/ezUUXXUQul6O7u5t58+bxpz/9iaamJvbs2XNsue7ubt7ylrdU6yUc4z19M7MhXH/99dx22238/Oc/56qrrgLgl7/8JY2Njfzwhz88ru2BAweOzedyOdrb2zn77LP5wAc+wPXXX88Xv/hF9u7dy0svvcSCBQuq+jrAe/pmNk7MpLKnbA4+sPLxpkyZwqOPPsp3vvMdZs2axezZs1m5ciXnnHNO2c8xZ84crr32WmbPns2iRYv43ve+V/Uzd6CMoZVrzUMrm2WTh1Yu33CGVvaevplZhrjom5lliIu+mdWteu9+rgfD3UYu+mZWlyZPnkxvb68L/yAigt7eXiZPnlz2Mj5l08zqUlNTE93d3fT09NQ6lbo2efJkmpqaym7vom9mdWnixIk0NzfXOo2Tjrt3zMwyZMiiL2mypE2Snpe0XdLXUvxMSeslvZRuzyhY5mZJHZJ2SbqiID5f0tb02D2qxcATZmYZVs6e/hHgP0XERUALsEjSQmA5sCEiZgEb0n0kzQbagDnAIuD7kvovO1sBLANmpWlR5V6KmZkNZciiH3mH0t2JaQpgCbAqxVcBV6f5JcCaiDgSEbuBDmCBpOnAtIh4KvKH4x8oWMbMzKqgrD59SRMkbQEOAOsj4hng3IjYB5Bu+wehaAT2FCzenWKNaX5gvNj6lklql9TuI/dmZpVTVtGPiL6IaAGayO+1F/+5mLxi/fQxSLzY+u6NiNaIaG1oaCgnRTMzK8Owzt6JiP8Afk2+L35/6rIh3faPJ9oNzChYrAnYm+JNReJmZlYl5Zy90yDp79P8FOC9wAvAOmBparYUeCTNrwPaJE2S1Ez+gO2m1AV0UNLCdNbODQXLmJlZFZRzcdZ0YFU6A+cUYG1EPCrpKWCtpE8BfwSuAYiI7ZLWAjuAo8BNEdGXnutGYCUwBXgsTWZmViUeT38UPJ6+mdUrj6dvZmYu+mZmWeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGTJk0Zc0Q9JGSTslbZf0uRS/RdIrkrakaXHBMjdL6pC0S9IVBfH5kramx+6RpLF5WWZmVsypZbQ5CnwpIp6TNBXYLGl9euyuiLijsLGk2UAbMAd4C/A/Jb0tIvqAFcAy4GngF8Ai4LHKvBQzMxvKkHv6EbEvIp5L8weBnUDjIIssAdZExJGI2A10AAskTQemRcRTERHAA8DVo30BZmZWvmH16UvKAXOBZ1Los5J+L+l+SWekWCOwp2Cx7hRrTPMD48XWs0xSu6T2np6e4aRoZmaDKLvoSzod+Bnw+Yh4lXxXzflAC7AP+HZ/0yKLxyDxE4MR90ZEa0S0NjQ0lJuimZkNoayiL2ki+YL/44h4CCAi9kdEX0S8AfwAWJCadwMzChZvAvameFOR+ElnJvlPuIFTroY5mZlBeWfvCLgP2BkRdxbEpxc0+yCwLc2vA9okTZLUDMwCNkXEPuCgpIXpOW8AHqnQ66grneS/wgycumqYk5kZlHf2ziXAx4Ctkrak2L8C10lqIV/POoHPAETEdklrgR3kz/y5KZ25A3AjsBKYQv6sHZ+5Y2ZWRcqfSFO/Wltbo729vdZpFCVKHJSoUHszs5GStDkiWgfGfUVuGXIU76OfWcOczMxGopzunczrwnvoZnZy8J6+mVmGuOibmWWIi76ZWYa46JuZZYiLvplZhrjom5lliIu+mVmGuOibmWWIi76ZWYa46JuZZYiLvplZhrjom5lliIu+mVmGuOibmWWIi76ZWYa46JuZZYiLvplZhrjom5llyJBFX9IMSRsl7ZS0XdLnUvxMSeslvZRuzyhY5mZJHZJ2SbqiID5f0tb02D2SNDYvy8zMiilnT/8o8KWIuBBYCNwkaTawHNgQEbOADek+6bE2YA6wCPi+pAnpuVYAy4BZaVpUwddiZmZDGLLoR8S+iHguzR8EdgKNwBJgVWq2Crg6zS8B1kTEkYjYDXQACyRNB6ZFxFMREcADBcuYmVkVDKtPX1IOmAs8A5wbEfsg/8EAnJOaNQJ7ChbrTrHGND8wXmw9yyS1S2rv6ekZTopmZjaIsou+pNOBnwGfj4hXB2taJBaDxE8MRtwbEa0R0drQ0FBuimZmNoSyir6kieQL/o8j4qEU3p+6bEi3B1K8G5hRsHgTsDfFm4rEzcysSso5e0fAfcDOiLiz4KF1wNI0vxR4pCDeJmmSpGbyB2w3pS6gg5IWpue8oWAZMzOrglPLaHMJ8DFgq6QtKfavwDeBtZI+BfwRuAYgIrZLWgvsIH/mz00R0ZeWuxFYCUwBHkuTmZlVifIn0tSv1tbWaG9vr2kOosTBhxo9j5nZUCRtjojWgXFfkWtmliEu+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGTJk0Zd0v6QDkrYVxG6R9IqkLWlaXPDYzZI6JO2SdEVBfL6kremxeySp8i/HzMwGU86e/kpgUZH4XRHRkqZfAEiaDbQBc9Iy35c0IbVfASwDZqWp2HOamdkYGrLoR8RvgL+U+XxLgDURcSQidgMdwAJJ04FpEfFURATwAHD1CHMet2YCKjHlapeWmWXIaPr0Pyvp96n754wUawT2FLTpTrHGND8wXpSkZZLaJbX39PSMIsX60glEiamrdmmZWYaMtOivAM4HWoB9wLdTvFg/fQwSLyoi7o2I1ohobWhoGGGKZmY20IiKfkTsj4i+iHgD+AGwID3UDcwoaNoE7E3xpiJxMzOrohEV/dRH3++DQP+ZPeuANkmTJDWTP2C7KSL2AQclLUxn7dwAPDKKvM3MbAROHaqBpNXAZcDZkrqBrwKXSWoh30XTCXwGICK2S1oL7ACOAjdFRF96qhvJnwk0BXgsTWZmVkXKn0xTv1pbW6O9vb2mOYhBDkCMo3WYWXZI2hwRrQPjviLXzCxDXPTNzDLERd/MLENc9M3MMsRF38wsQ1z0zcwyxEXfzCxDXPTNzDLERd/MLENc9M3MMsRF38wsQ1z0zcwyxEXfzCxDXPTNzDLERd/MLENc9M3MMsRF38wsQ1z0610uB9KJUy5X68zMbBwa8jdyrca6uqDYT1pK1c/FzMY97+mPVzNn+huAmQ3bkEVf0v2SDkjaVhA7U9J6SS+l2zMKHrtZUoekXZKuKIjPl7Q1PXaP5F3VUenszH8DGDh1ddU6MzOrY+Xs6a8EFg2ILQc2RMQsYEO6j6TZQBswJy3zfUkT0jIrgGXArDQNfE4zMxtjQxb9iPgN8JcB4SXAqjS/Cri6IL4mIo5ExG6gA1ggaTowLSKeiogAHihYxszMqmSkffrnRsQ+gHR7Too3AnsK2nWnWGOaHxgvStIySe2S2nt6ekaYYka5r9/MBlHpA7nF+uljkHhREXFvRLRGRGtDQ0PFkssE9/Wb2SBGWvT3py4b0u2BFO8GZhS0awL2pnhTkbglM8l/Mg6ccrt31zItMzvJjLTorwOWpvmlwCMF8TZJkyQ1kz9guyl1AR2UtDCdtXNDwTIGdJL/6jNw6nK3jJlV0JAXZ0laDVwGnC2pG/gq8E1graRPAX8ErgGIiO2S1gI7gKPATRHRl57qRvJnAk0BHkuTmZlVkaLY1Z51pLW1Ndrb22uagxjkAESl5HJF+90VUZl1S8Wv7DWzk5KkzRHROjDuYRjqRanhFszMKsjDMJiZZYiLflb4/H0zw9072dHZWTzuIZDMMsV7+mZmGeKib2aWIS76ZmYZ4qJvZpYhLvpmZhniom9mliEu+mZmGeKin3W+aMssU3xxVtb5oi2zTPGevplZhrjom5lliIu+mVmGuOibmWWIi76ZWYa46Ne5meR/rnHglBvzFftUTrOTkU/ZrKYSv4ML5ItsEZ0lnmrMT6j0qZxmJ6VR7elL6pS0VdIWSe0pdqak9ZJeSrdnFLS/WVKHpF2Srhht8uNO/+/gFptKFVkzswqqRPfO5RHRUvCr68uBDRExC9iQ7iNpNtAGzAEWAd+XNKEC6zczszKNRZ/+EmBVml8FXF0QXxMRRyJiN9ABLBiD9ZuZWQmjLfoBPC5ps6RlKXZuROwDSLfnpHgjsKdg2e4UO4GkZZLaJbX39PSMMkUzM+s32gO5l0TEXknnAOslvTBI22JHAKNYw4i4F7gXoLW1tWgbMzMbvlHt6UfE3nR7AHiYfHfNfknTAdLtgdS8G5hRsHgTsHc06zczs+EZcdGXdJqkqf3zwD8D24B1wNLUbCnwSJpfB7RJmiSpGZgFbBrp+s3MbPhG071zLvCw8udtnwo8GBG/lPQssFbSp4A/AtcARMR2SWuBHcBR4KaI6BtV9mZmNiwjLvoR8QfgoiLxXuCfSixzK3DrSNdpZmaj42EYbHg8PIPZuOZhGMap/jF5isU7x3LFHp7BbFxz0R+nOkvEXXrNbDDu3jEzyxAXfTOzDHHRt8rwAV6zccF9+lYZPsBrNi54T38s5HLF93pL/FCKmVm1eE9/LPT/WIqZWZ3xnv5Jpma/qWtm44KL/kmmk/x41QOnEr/MO/Z8gNesrrh7x8aWD/Ca1RXv6WdE3XX7+BuAWU14Tz8jOkvEa7a/7W8AZjXhPX0zswxx0TczyxAX/YwbN3397u83qwj36WdcZ4l4jjobrx/c329WAd7Tt6I6KX6+P9Twm4HP+DEbNRd9G5ZOanjxV2dnfniLgRP4w8CsTFUv+pIWSdolqUPS8mqvv6I8sNoxpY4NlJpylVz5cD8MfMzAMqyqRV/SBOB7wJXAbOA6SbOrmUNF9Q+sNnAarF/6JNVJ8W8ApSYY3ofEiD4oSn0YlJrAHxJ20qv2nv4CoCMi/hAR/xdYAyypcg5WBzoZ3ofESD8ohjV1dqKI8qfOzuIfTp2dw/vwqMbkDyhLqn32TiOwp+B+N/APAxtJWgYsS3cPSdo1wvWdDfx5hMsen1PJB0o+MpiK5VVhzmt4iubVxSDvl+o4Ma+urpG+VytpXP0d68Bo8yraz1ztol/sXXfCwPMRcS9w76hXJrVHROton6fSnNfwOK/hcV7Dk7W8qt290w3MKLjfBOytcg5mZplV7aL/LDBLUrOkvwPagHVVzsHMLLOq2r0TEUclfRb4FTABuD8ito/hKkfdRTRGnNfwOK/hcV7Dk6m8FP4tVzOzzPAVuWZmGeKib2aWISdl0a+XoR4kzZC0UdJOSdslfS7Fb5H0iqQtaVpcg9w6JW1N629PsTMlrZf0Uro9o8o5vb1gm2yR9Kqkz9dqe0m6X9IBSdsKYiW3kaSb03tul6QrqpzX/5D0gqTfS3pY0t+neE7S/ynYdv9W5bxK/u1qvL3+vSCnTklbUrwq22uQ2jD276+IOKkm8geIXwbOA/4OeB6YXaNcpgPz0vxU4EXyw0/cAny5xtupEzh7QOxbwPI0vxy4vcZ/xz+Rv8CkJtsL+EdgHrBtqG2U/q7PA5OA5vQenFDFvP4ZODXN316QV66wXQ22V9G/Xa2314DHvw3892pur0Fqw5i/v07GPf26GeohIvZFxHNp/iCwk/xVyfVqCbAqza8Crq5dKvwT8HJEdNUqgYj4DfCXAeFS22gJsCYijkTEbqCD/HuxKnlFxOMRcTTdfZr8NTBVVWJ7lVLT7dVPkoBrgdVjse5BcipVG8b8/XUyFv1iQz3UvNBKygFzgWdS6LPpq/j91e5GSQJ4XNLmNOwFwLkRsQ/yb0rgnBrk1a+N4/8Ra729+pXaRvX0vvsk8FjB/WZJ/1vS/5L07hrkU+xvVy/b693A/oh4qSBW1e01oDaM+fvrZCz6ZQ31UE2STgd+Bnw+Il4FVgDnAy3APvJfL6vtkoiYR37E05sk/WMNcigqXbj3AeAnKVQP22sodfG+k/QV4Cjw4xTaB7w1IuYCXwQelDStiimV+tvVxfYCruP4nYuqbq8itaFk0yKxEW2vk7Ho19VQD5Imkv+j/jgiHgKIiP0R0RcRbwA/YIy+1g4mIvam2wPAwymH/ZKmp7ynAweqnVdyJfBcROxPOdZ8exUotY1q/r6TtBT4F+CjkTqCU3dAb5rfTL4v+G3VymmQv109bK9Tgf8M/Ht/rJrbq1htoArvr5Ox6NfNUA+pv/A+YGdE3FkQn17Q7IPAtoHLjnFep0ma2j9P/iDgNvLbaWlqthR4pJp5FThu76vW22uAUttoHdAmaZKkZmAWsKlaSUlaBPw34AMR8VpBvEH537FA0nkprz9UMa9Sf7uabq/kvcALEdHdH6jW9ipVG6jG+2usj1LXYgIWkz8a/jLwlRrmcSn5r2C/B7akaTHwI2Briq8Dplc5r/PInwnwPLC9fxsBZwEbgJfS7Zk12GZvAnqBNxfEarK9yH/w7ANeJ7+n9anBthHwlfSe2wVcWeW8Osj3+fa/z/4ttf1Q+hs/DzwHvL/KeZX829Vye6X4SuC/DGhble01SG0Y8/eXh2EwM8uQk7F7x8zMSnDRNzPLEBd9M7MMcdE3M8sQF30zswxx0TczyxAXfTOzDPl/x4ppZmP02p0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# superimpose the distributions of the number of hits for both datasets\n", "nhits_antares = []\n", "i=0\n", "while i < len(dataset_antares): \n", " for row in dataset_antares:\n", " nhits = dataset_antares['nhit', i ]\n", " i=i+1\n", " nhits_antares.append(nhits)\n", "\n", "nhits_icecube = []\n", "i=0\n", "while i < len(dataset_icecube): \n", " for row in dataset_icecube:\n", " nhits = dataset_icecube['nch', i ]\n", " i=i+1\n", " nhits_icecube.append(nhits)\n", "\n", "plt.hist( nhits_antares , bins=50 , range=[0,200] , label = \"Antares\" , fill=False , color = \"red\" , histtype = 'step' )\n", "plt.hist( nhits_icecube , bins=50 , range=[0,200] , label = \"IC-40\" , fill=False , color = \"cyan\" , histtype = 'step' )\n", "plt.legend()\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.9.5" } }, "nbformat": 4, "nbformat_minor": 5 }