radius_measurement_shuffled_over_order_of_events.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "75ea49c9",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import glob\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib import rcParams\n",
    "rcParams['figure.figsize']=(11,6)\n",
    "rcParams['axes.labelsize'] = 20\n",
    "rcParams['axes.titlesize'] = 20\n",
    "rcParams['xtick.labelsize'] = 20\n",
    "rcParams['ytick.labelsize'] = 20\n",
    "rcParams['legend.fontsize']=20\n",
    "\n",
    "rcParams[\"font.family\"] = \"serif\"\n",
    "rcParams[\"font.serif\"] = \"STIX\"\n",
    "rcParams[\"mathtext.fontset\"] = \"stix\"\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "be542fe4",
   "metadata": {},
   "outputs": [],
   "source": [
    "#injected radius=11.86984\n",
    "\n",
    "radius_accuracy={'10m':0.00084,'50m':0.0042,'100m':0.0084}\n",
    "\n",
    "network={'CE40+2A#':'radius_data_files/output_eos895_CE40_2Asharp',\n",
    "         'CE40+CE20+A#':'radius_data_files/output_eos895_CE40_CE20_Asharp',\n",
    "         '3A#' : 'radius_data_files/output_eos895_Asharp'\n",
    "        }"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "80c1ebfb",
   "metadata": {},
   "source": [
    "## Plotting Histogram for number of years to reach radius accuracy of 10 m"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "4d08c9b2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAncAAAGVCAYAAACcrtd2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAsqklEQVR4nO3de5xdVXn/8e+DciuXQALWJFSoRBAwiiSAqKioofAzSVuVitQqlBKUm4DWKiQFTLTeCCgIJhCItlK0VuoEEEEEL4iQEEHxThQUEhUCDMhFQZ7fH2sdcjg5M5mVOWv2ZK3P+/Var82cyz7PnPkyeWavtfcxdxcAAADKsFHTBQAAAKB3aO4AAAAKQnMHAABQEJo7AACAgtDcAQAAFITmDgAAoCDPbroAANgQmNmukg6WtLWkVZK+7O6rm60KANbGkTtgFDKzbc3sVDO73MxWmNlPzOwRM/O28Roz29fMzmm77ddmdn0c3zKz5Wb2aLxvp0Fe741mdmeX283MjjWzPjP7hJl91sz+sYff54vN7Aoz+52Z3Rv/+yVDeN44M/uOmdkQX+frZnZc/G8zswfMbP8hPnczM/uMpBsk3SdprrtfsK7Gzsy2N7PPxJ/JH8zsBjN79RBe72Qzu34otQFANxy5A0YZM3uVpC9LepakoyUtcffHzGxjSQdKWiTpLyXJ3W+SdFN8zoslHeXuX+vY39YKjclArzdB0gWSHu5y93xJb5Y02d0fNLOxkn5gZtu7+9mD7HNbSVMk3eDujw3wmBdLulLSpZJukrS/pP8n6TVm9lJ3//lA+5c0S9IrJB0g6RuDPK7lhZJ+Gv97O0nbSPrZup5kZltIulrSbpL2dvcVQ3gtmdmWkv5b0hcVfl4vkPRxSVea2a7ufvcAz5ss6T8k3TiU1wGAbmjugFHEzHaQdLnCUfWXtDcT7v6EpCvMbC9Jyzue+mDc/rFzn+7+kJmdJ+nPXV7PJJ0t6YeSduq47wWSTpR0ors/GPd1v5mdKeljZvYFd181wLfyEknXSPprSXcO8JjTJL3K3X/Z9pofkfRvko6PYy2xyT0qfvleraO5M7OtJE3UmuZuF4X3697BnhctkPRySW8YamMX/bOkU2PzLUlLzWwTSRdLmipprebOzDZVaOx+kvA6ALAWpmWBUSI2WudL2krSmQM1E+6+UtK8zpsH27e7n+/uv+ly17GSrlD3BmxW3H694/ZrJW0s6ZjBXnMwseH6ZntjF306bnce5OmHSPqSpB9IOtjM9ljHy+0q6RFJ98Svd5H0M1/HZy+a2b6S/lHSUne/ch2v0enitsauZSNJT2ngo3JzJZ2pNY06AKwXmjtg9HiepOnxvz+3jsd+QdIdQ9mpmZ01wO27Sdp3kNc6IG5/2nF768jSq4by+t24+8Pu/qkudz0Qt7/o9rzYAB+v0AR+Mt78ngEeu9jMXNJSSVtIeip+fZGkfeM6xMXxsW80s5Vm9ndtuzgybn9pZl+NawLvMrOLzGzcur6/jlo2kfRGSW909991qfW18XnXDbbftsdPMLOvxO/hF2Z2YPu+zOzBeBS0dduOZnaBmX3ezH4Qv4et2u5vra38rJl9yMyuM7O5ZvbseP82ZnZUXMf5ajObbWarzGxevP9v4/5PMLNPmtnNQ/k+AOTBtCwwerw4bp+QdNdgD3T3oUwpyswmKqyZO6nj9k0kfUzS4e7uA5yXsIOkh9z9GdO57v6EmT0iaceh1JDoRQpHty4a4P79JK1097vM7BJJH5X0NjM7tcsU8fEK07b/pnAksHUkcrHCUb/5WjONPVHSeIUGu+U1cfsrSafHuv5WoRne18xe4u5PrusbMrOXxtcaL+m3Xe7fVtLJCj+nIXH3lWb2doUjrlvomVPT34+3nxL3/3yFKfID3P3Xcer/Tkm7mtn+7v6UwjT3uZImufsKC2cG/1QhhxcqTLO/WWFd5JGSfqTQNG9lZtsp/LGxa/y5PEvSd4f6vQDoPY7cAaPHbnH74FCahgFcaGa3xvEzhcZkhy6PO03Sues443Osup9kIUn9Cs1K62zSbdqHpC3j47buvK91NGgA75T0AXe/bYD7T5T0KUly98cVprE3Vpf1efHo4H2xzh+7+33x6wmSbo9fPxwfe46k7TuOJk5UmCI91d3/6O5PuPuX4mvuLumtg3wfiu/NKyS9K365q6QbzezvOx42X9Ip8fsZMnfvV3gvxuuZjeEhki6JTZskfUTSZe7+6/i8uxVOsHm51hydNUkrJLXy0DqZZUp8zjcVpsIl6TZ3/6i7z3T3dys0+Zu29hX/GPhAyvcCoLdo7oDR49G43T6epbk+/sXd94xjV0nPUThZ4mkWLgGyRedZtV3cpzVNWqcttKYReL/CdGr7WBLvu63Lfa/stkMze2t8vY8PcP+OCkfgvtV28/kKRzrfFc9Q7WaS4hR2nNbdWVLnWj/Fxq/dnyXd09YktbTW3+01wOu17/MGd5/l7gdIeq2kJyXNafueDlNoNH+wrn0N4FOS/iDpfbbm8OtbJf1X3L9JOkjSy+M09eI4Ff2Ews/mubHOBe4+KZ4RPUlrmuXN217ribi9paOGHyocCb3IwqVfxrn7UM5gBpAJ07LA6LGs7b93VvgHc1jiP9btzcTWkmZLequZbdb20GeFu8Nt8SjSPZKmmNlG7Q1OPPI2RmvWxV2ktU+6eKlC4/EPChf8bffDjq9lZi+SdISkNw1yosNxCs3q9zumkR9TuLTJP8fXlJk9T2vev60lTTazsxWOUG0t6Roz+7OkF7eOaHVxp8JlUzq1vp+xAzyvK3e/zsy+LGlGrPGvFBqxt3T8LDaK928WnuZrnQHdts/VFs6Efp+k15nZCkmPxJNuJGl7hRN0vjjYpWvi671S4QjgNQpTtJ8c7PFtNfzJwqV4FilcuucQM3u7u18xlOcD6D2aO2D0WKowHbaLQqNz0uAPHxp3/4r09JGvFylcK2+g6djWNelM0lUKl+3YRc88qeIFcfv1uP9fS3pGg9Q29brU3e8crL7Y5HxA0ps7T0Roe8zWCick7Oruj3bc97pYy0lmdl6c0l4paU+FKelvK5w48geFk0A+Hr9WfNxAbpR0pJlt6e5/aLt9k7hNuTRKy2/ikKTXKZxA88gAj31MYc3bTuvY53xJJyg0eN9UWFPY8oDCEcjduz3RzLZ19wfiEcSFknZrnVU9wDrMbvvYWOG9PSSOj0u63MymuXtn0w9gBDAtC4wSca3S8QoL94+PC/G7MrMDzGx6+02D7dvMNlL4R//bCicldI4rFRb7t76WpM/H7es7djctbhcP+g0NQTzh46OSjnH3h9pu3zhO07Ycp/BxX4927kPhZIJbFZqgN0mSuz8Zm8otJd3t7rfHr7dWWH93ZxxPr220tc+AXaxwRHN6x+2ttZF9bc/dOjY56zJFa97XK9T9Z/H9OPaT1Lk+by3x7NsLFX4ub9OaKfHWtRG/L+lwCxdIfpqZvUzhRA5J+neFafi7432tZQFDOQDwCkmv9+CLCn8QPKxwQWoADeDIHTCKuPvVZvYPCtNi3zWzDyqsn7pb4cSB3SW9Q2ER/Tvanrp13G6lDhY+VeIsSf2xgfpel8fcK+mP7v70fe7+UzM7V9J7zew/3b0/ntn5XknnuPs6P+FhMGY2XmEK8KsKzWzrrk0U1qd9ND5uS4UmpLPJatXpZvY/Ckfq5prZZe7+p3j3i/TMaeAXqstFgs3sBEmfNLPj3f3cePP3JF0S93lNnALdRKHR/JS73xqf+1yFNX23S3pZvO2DCkcN57v77fG2f1E4inZmrPtedbmQspk9FO9f6+c0iE/Fur7eZRp3tsJR2BstfIzaryTtI2kPhfV4krStwpT3O83sNkkzFY4cTjWz1yt8gkgrW+1TyC3vMbOr45T6gwpH8jqv8wdghNDcAaOMu/+vmV2j0NAcoPCP9rMkPaQwdftZd79aksxsb4V1bXvGp3/ezH6qMB23qcI/2nvE579L6d6tsPbsC2b2c4UzPs+OY73FM2q/odBs7dblIb9R+KiuveJrjVNoMs9092dcBNjCtela13l7gaRrzezYeJLCHgpNV8sLJX1lgNf7rdqmaWPT+A6F6c6vmNmP412fVziRo+UxhfWJ7dcd/L3C+rUbzeyLcb8/k3TwMM6EHpC7/8LMfqIuR1Pd/WsWroN3hsKauAcVjhq+oe0kkuMU1tjNlfQ1hWnejRQukfIGhY+7OzY+draZ7eTun2l7mWmSrjOz6xSm8c9U+Og1AA2wgdcuA8D6MbMxkiZLWpZ6iQ+ki9Pbl0nad5ATUgBUguYOADZAZmatRs7MLpR0vbv/V8NlARgFaO4AYANkZksVLtVyl8I6uFfmmPIFsOHhbFkA2DAtU1iL+JjC9QFp7ABI4sgdAABAUZLOlt1uu+18p512ylQKAAAAhuqWW265z92377w9qbnbaaedtGzZsnU/EKPOknhZ0xkzmq0DZSJfyIl8Ad2Z2V1db0+Zlp06darT3G2YWteHZRYeOZAv5ES+gO7M7BZ3n9p5OxcxrsT0rtf2B3qDfCEn8gWk4cgdAADABmigI3dcCgUAAKAgNHcAAAAFobmrhNmaRclAr5Ev5ES+gDQ0dwAAAAXhbNlKcAkB5ES+kBP5AtJw5A4AAKAgNHcAAAAFobmrxIwZfHQP8iFfyIl8AWm4iHEl+Pge5ES+kBP5Arrj48cq19fXdAUoGflCTuQLSMOROwAAgA0QR+56ZNo0qb9/+PuZPFlatGj4+wEAAGhHc5eov1+aM2f4+1m4cPj7WJ/XmzVrZF8XdSBfyIl8AWmYlk20zz69a+6WLBn+foaKBcnIiXwhJ/IFdMe0bOWOOqrpClAy8oWcyBeQhuauEiM9DYy6kC/kRL6ANFzEGAAAoCA0d5VYuTIMIAfyhZzIF5CGadlKTJwYtixIRg7kCzmRLyANzV0lxo9vugKUjHwhJ/IFpKG5q8SGPKXBhaNHvw05Xxj9yBeQhuYOo96GeuFoAACawAkVAAAABaG5q8SUKWEAOZAv5ES+gDRMy1Zi+fKmK0DJyBdyIl9AGpq7SlT+kcDIjHwhJ/IFpKG5qwRTGsiJfCEn8gWkYc0dAABAQWjuKnH66WEAOZAv5ES+gDQ0d5U444wwgBzIF3IiX0Aa1txV4rTTmq4AJSNfyIl8AWlo7irBlAZyIl/IiXwBaZiWBQAAKAjNXSVuuSUMIAfyhZzIF5CGadlKTJ0atu7N1oEykS/kRL6ANDR3ldhrr6YrQMnIF3IiX0AamrtKMKWBnMgXciJfQBrW3AEAABSE5g4AAKAgNHeVmDAhDCAH8oWcyBeQhjV3lVi1qukKUDLyhZzIF5CG5q4S99zTdAUoGflCTuQLSENzVwmmNJAT+UJO5AtIw5o7AACAgtDcVWLWrDCAHMgXciJfQBqau0pccEEYQA7kCzmRLyANa+4qsWBB0xWgZOQLOZEvIA3NXSWY0kBO5As5kS8gDdOyAAAABaG5q8SSJWEAOZAv5ES+gDRMy1Zi5sywdW+2DpSJfCEn8gWkobmrxPTpTVeAkpEv5ES+gDQ0d5VgSgM5kS/kRL6ANKy5AwAAKAjNHQAAQEFo7iphFgaQA/lCTuQLSENzBwAAUBBOqKgElxBATuQLOZEvIA1H7gAAAApCcwcAAFAQmrtKzJgRBpAD+UJO5AtIw5q7Slx+edMVoGTkCzmRLyANzV0l+vqargAlI1/IiXwBaWjuKsGUBnIiX8iJfAFpWHMHAABQEJq7SixcGAaQA/lCTuQLSMO0bCWOPjpsZ81qtg6UiXwhJ/IFpKG5q8RRRzVdAUpGvpAT+QLS0NxVgikN5ES+kBP5AtKw5g4AAKAgNHeVWLkyDCAH8oWcyBeQhmnZSkycGLbuzdaBMpEv5ES+gDQ0d5UYP77pClAy8oWcyBeQhuauEkxpICfyhZzIF5CGNXcAAAAFobkDAAAoCM1dJaZMCQPIgXwhJ/IFpGHNXSWWL2+6ApSMfCEn8gWkobmrxLJlTVeAkpEv5ES+gDQ0d5VgSgM5kS/kRL6ANKy5AwAAKAjNXSVOPz0MIAfyhZzIF5CG5q4SZ5wRBpAD+UJO5AtIw5q7Spx2WtMVoGTkCzmRLyANzV0lmNJATuQLOZEvIA3TsgAAAAWhuavELbeEAeRAvpAT+QLSMC1bialTw9a92TpQJvKFnMgXkIbmrhJ77dV0BSgZ+UJO5AtIQ3NXCaY0kBP5Qk7kC0jDmjsAAICC0NwBAAAUhOauEhMmhAHkQL6QE/kC0rDmrhKrVjVdAUpGvpAT+QLS0NxV4p57mq4AJSNfyIl8AWlo7irBlAZyIl/IiXwBaVhzBwAAUBCau0rMmhUGkAP5Qk7kC0hDc1eJCy4IA8iBfCEn8gWkYc1dJRYsaLoClIx8ISfyBaShuasEUxrIiXwhJ/IFpGFaFgAAoCA0d5VYsiQMIAfyhZzIF5CGadlKzJwZtu7N1oEykS/kRL6ANDR3lZg+vekKUDLyhZzIF5CmiuZu2jSpv783+xozpjf7GWlMaSAn8oWcyBeQpormrr9fmjOn6SoAAADy44QKAACAgtDcVcIsDCAH8oWcyBeQhuYOAACgIFWsuQOXEEBe5As5kS8gDUfuAAAACkJzBwAAUBCau0rMmBEGkAP5Qk7kC0jDmrtKXH550xWgZOQLOZEvIA3NXSX6+pquACUjX8iJfAFpaO4qwZQGciJfyIl8AWlYcwcAAFAQmrtKLFwYBpAD+UJO5AtIw7RsJY4+OmxnzWq2DpSJfCEn8gWkobmrxFFHNV0BSka+kBP5AtLQ3FWCKQ3kRL6QE/kC0rDmDgAAoCA0d5VYuTIMIAfyhZzIF5CGadlKTJwYtu7N1oEykS/kRL6ANDR3lRg/vukKUDLyhZzIF5CG5q4STGkgJ/KFnMgXkIY1dwAAAAWhuQMAACgIzV0lpkwJA8iBfCEn8gWkYc1dJZYvb7oClIx8ISfyBaShuavEsmVNV4CSkS/kRL6ANDR3lWBKAzmRL+REvoA0NHfAepg2Terv782+Jk+WFi3qzb4AAKC5q8Tppz9zi+Hp75fmzOnNvkr4UHTyhZzIF5CGs2UrccYZYQA5kC/kRL6ANBy5q8RppzVdAUpGvpAT+QLS0NxVgukM5ES+kBP5AtIwLQsAAFAQmrtK3HJLGEAO5As5kS8gDdOylZg6NWzdm60DZSJfyIl8AWlo7iqx115NV4CSkS/kRL6ANDR3lWBKAzmRL+REvoA0rLkDAAAoCM0dAABAQWjuKjFhQhhADuQLOZEvIA1r7iqxalXTFaBk5As5kS8gDc1dJe65p+kKUDLyhZzIF5CG5q4STGkgJ/KFnMgXkIY1dwAAAAWhuavErFlhADmQL+REvoA0NHeVuOCCMIAcyBdyIl9AGtbcFWDaNKm/f/DH7Lhj2O6zz+CPmzxZWrRoZGoaqjFjerOf0erxx9f9cxmKXv3s1seCBc28LtbWy//3msxUO/IFpKG5K0B/vzRnTm/2tXBhb/bTy5pKd8IJvdlPr35264Mps9FjNP4+GC7yBaRhWhYAAKAgNHeVuPnmMIAcliwJA8iBfAFpmJatxLx5YdvX12wdKNPMmWHr3mwdKBP5AtLQ3FVi772brgAlmz696QpQMvIFpKG5qwQnNyAnpsyQE/kC0rDmDgAAoCA0dwAAAAVhWrYSrQXJ6zqholcX1C39wsN4JrOw7VzwXuIFdTHyBsoXgO5o7vAMvbqgLiCVeUFdABjtaO4qwSVQkBNHVJAT+QLSsOYOAACgIDR3AAAABaG5q8TcuWEAOcyYEQaQA/kC0rDmrhJLlzZdAUp2+eVNV4CSkS8gDc1dJWbPbroClIwTdpAT+QLS0NxVohfXrgMGwpQZciJfQBrW3AEAABSE5q4SV10VBpDDwoVcZBj5kC8gDdOylTjvvLA96KBm60CZjj46bGfNarYOlIl8AWlo7ipx4IFNV4CSHXVU0xWgZOQLSENzV4njjmu6ApSMKTPkRL6ANKy5AwAAKAjNXSVWrw4DyGHlyjCAHMgXkIZp2UoccUTYcjFQ5DBxYti6N1sHykS+gDQ0d5UYO7bpClCy8eObrgAlI19AGpq7Sixe3HQFKBlTZsiJfAFpaO4a8vjjvftIsDFjerOf0pX+nvfy+5s8WVq0qDf76pXSvz8MzbRpUn9/b/ZFDlAqmruGnHBC0xXUp/T3vJff32i89ETp3x+Gpr9fmjOnN/siBygVZ8tW4qSTwgBymDIlDCAH8gWk4chdJVasaLoClGz58qYrQMnIF5CG5q4S8+c3XQFKtmxZ0xWgZOQLSENzV4lJk5quACVjygw5kS8gDWvuAAAACkJzV4lLLgkDyOH008MAciBfQBqau0pcemkYQA5nnBEGkAP5AtKw5q4Shx7adAUo2WmnNV0BSka+gDQ0d5U47LCmK0DJmDJDTuQLSMO0LAAAQEFo7ipxxx1hADnccksYQA7kC0jDtGwlTj45bPv6mq0DZZo6NWzdm60DZSJfQBqau0rsvHPTFaBke+3VdAUoGfkC0tDcVeKss5quACVjygw5kS8gDWvuAAAACkJzBwAAUBCau0ocfngYQA4TJoQB5EC+gDSsuavE/fc3XQFKtmpV0xU0Y9o0qb+/N/uaPFlatKg3+xqNhvNetfK1zz7SmDG9q+nxx8M+h6v0nx02PDR3lbj44qYrQMnuuafpCprR3y/NmdObfS1c2Jv9jFbDea9Wrw7bceN6V48knXBCb/ZT+s8OGx6au0r0+pci0I4pM+TE7y8gDWvuAAAACkJzV4lzzw0DyGHWrDCAHPj9BaShuavE1VeHAeRwwQVhADnw+wtIw5q7ShxzTNMVoGQLFjRdAUrG7y8gDc1dJQ46qOkKUDKmZJETv7+ANEzLAgAAFITmrhI33xwGkMOSJWEAOfD7C0jDtGwl5s0L276+ZutAmWbODFv3ZusYql59MkEvPy0BA+P3F5CG5q4Se+/ddAUo2fTpTVeQplefTICRwe8vIA3NXSV69RFJQDdMySInfn8BaVhzBwAAUBCaOwAAgILQ3FVi5sw1i96BXjMLA8iB319AGpo7AACAgnBCRSW4hABy2lAugYINE7+/gDQcuQMAACgIR+4ArIWL/I483vMNV69+dpI0ebK0aFFv9oV60dxVYu7csOV6URiK1Iv8kq/h48LKAxvt+erlz27hwt7tC/WiuavE0qVNV4CSkS/kRL6ANDR3lZg9u+kKUDLyhZzIF5CG5q4SvVoPAnRDvpAT+QLScLYsAABAQWjuKnHVVWEAOZAv5ES+gDRMy1bivPPC9qCDmq0DZSJfyIl8AWlo7ipx4IFNV4CSkS/kRL6ANDR3lTjuuKYrQMnIF3KqKV9cEBm9QHMHAMAowQWR0QucUFGJ1avDAHIgX8iJfAFpOHJXiSOOCNu+vmbrQJnIF3IiX0AamrtKjB3bdAUoGflCTuQLSENzV4nFi5uuACUjX8iJfAFpWHMHAABQEJo7AACAgtDcVeKkk8IAciBfyIl8AWlYc1eJFSuargAlI1/IiXwBaWjuKjF/ftMVoGTkq0y9/LSEMWPW/7nkq1nTpkn9/U1XkU+Jn+RBc1eJSZOargAlI19l6uWnJQwH+WpWf780Z07TVeRT4id5sOYOAACgIDR3lbjkkjCAHMgXciJfQBqau0pcemkYQA7kCzmRLyANa+4qceihTVeAkpEv5ES+gDQ0d5U47LCmK0DJyBdyIl9AGqZlAQAACkJzV4k77ggDyIF8ISfyBaRhWrYSJ58ctn19zdaBMpEv5ES+kFMvL9Y9Wi6ITHNXiZ13broClIx8ISfyhZx6ebHu0XJBZJq7Spx1VtMVoGTkCzmRLyANa+4AAAAKQnMHAABQEJq7Shx+eBhADuQLOZEvIA1r7ipx//1NV4CSkS/kRL6ANDR3lbj44qYrQMnIF3IiX0AamrtKjBvXdAUoGflCTuQLSMOaOwAAgILQ3FXi3HPDAHIgX8iJfAFpaO4qcfXVYQA5kC/kRL6ANKy5q8QxxzRdAUpGvpAT+QLS0NxV4qCDmq4AJSNfyIl8AWmYlgUAACgIzV0lbr45DCAH8oWcyBeQhmnZSsybF7Z9fc3WgTKRL+REvoA0NHeV2HvvpitAycgXciJfQBqau0rMmdN0BSgZ+UJO5AtIM6qbu2nTpP7+4e9nzJjh7wMAAGBDMKqbu/5+/mIDAABIwdmylZg5MwwgB/KFnMgXkIbmDgAAoCCjeloWvcMlBJAT+UJO5AtIw5E7AACAgtDcAQAAFITmrhJz54YB5EC+kBP5AtKw5q4SS5c2XQFKRr6QE/kC0tDcVWL27KYrQMnIF3IiX+vn8celffYZ/n74IIAND81dJXrxPzgwEPKFnMjX+jnhhKYrQFNYcwcAAFAQmrtKXHVVGEAO5As5kS8gDdOylTjvvLA96KBm60CZyBdyIl9AGpq7Shx4YNMVoGTkCzmRLyANzV0ljjuu6QpQMvKFnMgXkIY1dwAAAAWhuavE6tVhADmQL+REvoA0TMtW4ogjwravr9k6UCbyhZzIF5CG5q4SY8c2XQFKRr6QE/kC0tDcVWLx4qYrQMnIF3IiX0Aa1twBAAAUhOYOAACgIDR3lTjppDCAHMgXciJfQBrW3FVixYqmK0DJyBdyIl9AGpq7Ssyf33QFKBn5Qk7kC0hDc1eJSZOargAlI1/IiXwBaVhzBwAAUBCau0pcckkYQA7kCzmRLyANzV0lLr00DCAH8oWcyBeQhjV3lTj00KYrQMnIF3IiX0AamrtKHHZY0xWgZOQLOZEvIA3TsgAAAAWhuavEHXeEAeRAvpAT+QLSMC1biZNPDtu+vmbrQJnIF3IiX0AamrtK7Lxz0xWgZOQLOZEvIA3NXSXOOqvpClAy8oWcyBeQhjV3AAAABaG5AwAAKAjNXSUOPzwMIAfyhZzIF5CGNXeVuP/+pitAycgXciJfQBqau0pcfHHTFaBk5As5kS8gDc1dJcaNa7oClIx8ISfyBaRhzR0AAEBBaO4qce65YQA5kC/kRL6ANDR3lbj66jCAHMgXciJfQBrW3FXimGOargAlI1/IiXwBaWjuKnHQQU1XgJKRL+REvoA0TMsCAAAUhOauEjffHAaQA/lCTuQLSMO0bCXmzQvbvr5m60CZyBdyIl9AGpq7Suy9d9MVoGTkCzmRLyANzV0l5sxpugKUjHwhJ/IFpGHNHQAAQEFo7gAAAApCc1eJmTPDAHIgX8iJfAFpaO4AAAAKYu4+9Aeb3SvprnzljArbSbqv6SIKxvubH+9xXry/efH+5sX7m99Ivsc7uvv2nTcmNXc1MLNl7j616TpKxfubH+9xXry/efH+5sX7m99oeI+ZlgUAACgIzR0AAEBBaO7WtrDpAgrH+5sf73FevL958f7mxfubX+PvMWvuAAAACsKROwAAgILw2bJojJk9W9KL3P3WpmsphZltJOlvJO0j6VFJP3f3rzRbFQCUycw2d/fHmq6jE81dBzPbRNKNks5x98UNl1MUM/umpFd13PxPkm4d+WrKY2aTJf2npMckvd3df9FwSUUws/MkvWuAu1e7+3YjWU+JzOy1ko6QdLekbST9laQPufuNTdZVEjN7iaR3SnpI0l8rvNenjsbGZENgZltLOk7SeySN63L/JElzJT0s6U+SNpE0291/PxL10dyt7XRJezVdRGnM7BWSNpV0ZtvN/ZK+2ExFZTGzV0m6QtLVkt7i7k82XFIRzGyMpOmSLpD0C0nti5TfofCHIIbBzF6v8HtgD3dfFW97taSvm9me/JEyfPH9/F9Ju7v77+OsyRWSvmNm+7v7o81WuGGJf4y8VtLJkjbvcv/zJd0kaZ67nxVve5+kb5rZAe7+2+w1ckLFGmb2SoUu/O8kHcGRu94xsyWS3u/uP2q6ltKY2XhJtyn8RT7F3fsbLqkYZvbPku5w92913L6xpN9KOtTdr2mkuEKY2SWSdum86KuZPSDpDHc/u5HCChEbuVsl/dDd39p2+8sl3SDpve5+5gBPxyDM7LuS9nN367j9vyQdLOk57v7neNuzJN0r6fPufnzu2jihIoqHWD8g6ZSmaymNme0p6UBJp5jZ4WY2seGSSjNH0vaS5tPY9dzFnY1d9HqFo3jXjXA9Jdpc0u5mNqF1g5ntoDA9u6qpogqyq6Q9JP204/abFP4gPHnEKyrHnzpvMLNxkt4s6Rutxk6S4n9fL+kIM9smd2E0d2t8QuEfSdYf9N5MSU9JOkzSxZLuNrOLzWyzZsva8JnZppLeHr98rpktN7NHzOz7ZvZOM7PBno/B+cBTG4dK+jLT3z1xvsJ6pOvNbM/4e+Fjkj6nMJWI4Xlh3D7efmNsNn4raYKZbTviVZVrP4UlSD/uct9PJG0h6aW5i6C5k2Rmb5L0K3df3nQtJXL3DyoEemeFo6P9kg6XtLi5qooxReG9fVhhvd0USTtIWq7wj+b7myutTLH5+HtJ/9N0LSVw96slHSJpR0nfk/RdSV9w93fQPPfEQ3G7Q5f7VsftpBGqpQat9/nBLvc9ELc75i6i+uYuTgX8k8JfisjE3Z9y91+6+0cUpgl+JektZrZfw6Vt6FpT3Be5+3c8eEDh7M7fK0yF/0Vz5RXpYIXpGKZke8TdL5P0SUmmcFTj38ws+z+AlbhNIa8z4rqvdq2v/ziyJRVtbNw+3OW+1rKZ8bmLqLq5i9cE+6Skk9rnxpGXu/9Oa6YS922ylgK0cnt3+43u/ieFI3lbSnrBSBdVOKZke8iC2ZL+UtLzJF2mMLX1TdbnDl+89MaZCkeL5pvZFma2jZkdKWnP+LC7B3o+kt0Xt1t2uW+LuF3d5b6eqrq5k3SspGslrTKzzeJ0y6bxvo3jbVwuJo+bFP5a3KTpQjZwd8bt2C73rRrkPqwHM9tS0gxxCZ9eOlXSCZKOjH/4vVnSxxWakdMbrKskpypck22KpJsVrof5iEIPsMzd72+wttLcE7fbdLlvbMdjsqm9uXuTwrqkx9pG64yihfHr2c2UVrynFI46cQ2r4fmJwpqa53e5r9U4rxi5coo3Q+EfxesbrqMkb5H0vdaRUHd/SmGt6CpJL2+ysFLE5RqfdvdXuvse7j5DoXl+tqQFDZdXmm8p9A4v7HLfrgonttyQu4jaj0odI2nrjtvGS/qypHkKF3nkcHUer1VYcPq1huvYoLn7Y2Z2qaRDzWxTd29fO7ObpO9L+k0z1RWJKdneu0/SVu03uPtTZnazulwgFsNnZjtJep+kpeLEtp5y94fNrE/Sa81so/jHSmsZ2OskXebuD+auo+ojd+7+Y3f/XvtQ+MdQklbE22juhsHMXmZmvzOzL7TWz8SPZfmwpEO4MnpPfFBhivvDrRviRw29TNK7B7mcBxLEa1MdLM6S7bVPSdq//eQqM/tLhaN2n26sqkLFT135ksIFdWfwh8qwbCk9fVHzdqcqzJy0X6z4RIUTWOaMRGG1H7lDfr+U9G2FBdLLzexqSb+WNLP1UUMYHne/x8z2lfQRM/uawvWVxkra391/0Gx1Rfl7hbPdrm+4jqK4+2VmNl3SB83s1wpH9J8n6R3u/tVGiyuMmR0g6SJJ35F0rLs/tI6noIv4cW5vUljDKEkXmtlX3f1SSXL3FfGPlXnxIv6mcHR6P3cfkWUyfPwYAAAFM7PDJe2i0Dhf6e63N1oQsqO5AwAAKEjVa+4AAABKQ3MHAABQEJo7AACAgtDcAQAAFITmDgAAoCA0dwAAAAWhuQPQE2Z2tpm5mT1iZneY2e1m9kS87edm9jMzuz9+fWt8ztj42HMaLl+SZGZ3mtnD8Xs53cx2H6HXPcXM5pnZD+P785qReF0AZeITKgD00sUKV75/TArNksIHlB/o7nfGz1c8Og4pXLl9o7gdLa509xNH8gXd/cOSZGY3Srp8JF8bQHlo7gD0ikv6QKux6/qA8CHa55vZq+LXqyU9f4TqG6oB6x8BjzT42gAKwbQsgF65wN1/N8THfjBrJQBQMZo7AD3h7j9OeOxPJMnMNjWzN5jZh1v3mdmLzGyumf3IzPY0s381s1+Y2e/N7KPxMS83s6vM7AEzu9nMJrXv38zGmdm5Znaxmd1mZl8ys+euz/dlZq8zs1VxLdxFZrZzvH1zM/toXEf4CjN7rpmdZGbLzez5ZnZOrPmd8fFHxZqOM7MLzaxvfeoBgHVhWhZAI2JDdpykd0u6S9Ip8a5HJE2StLuk2ZI+J2k/SSdKOtXMtpN0q6QjJO0g6duS/kPSIXG/20j6nqS3uPtyM9tK0i8lfcfMXuLuSVOf7n6tmb0v1nGfu6+Itz9mZndJ+ri732Bmb5D0NkkvlfSvkn4u6UeStjSzyZLOkvQcd3801vSVlDoAYKg4cgegEe5+Rzxx4b6O238l6dr45Wfcvc/d71NoriTpL9z9HHdf5e5LJS2XtH/bLk6RdJu7L4/7e1jSlZJ2lnToepb735J+JelIM9ui7fa3SvpsfJ0rJF0Vb/+6u5/l7ge4+yck7SJpC0n7tNU0bz1rAYBB0dwBaFq3I2lPdmwl6Y64/WPHY++QtH3b1wdL2sXMFrdGvP82SWPXp0B3f1Lh6OBYSUdKTx95fMjdV7Y99Im4vaVjF9+QtErS18zsQ2a2hbt/Y31qAYB1YVoWwAbB3Z8y63rFlKf0zD9Ud5a0MMPlTD4n6d8lnWxm50v6R0mLh/JEd3/AzKYqHAE8RdLbzOwf3P2mHtcIABy5A1CcexXW663FzLZd3526+x8lfUzhun2HKBwhHNJJEWa2WTzC9zpJ75S0paRvm9lu61sPAAyE5g5ATq3fMSN5keKbJU0zs79pv9HMXiDpQ8Pc94UKzeNZkpbFhm8oDjWzXd39SXdfIOkVCjMn04ZZDwCsheYOQBbx0yjGxC/HDPLQzSVtHx/fMi5uN2vb39bxP8d3PH/beP+m8eszFNa+XW5m55vZMWa2QNJlkuYnfyNt4gWaF0h6jsKncXTaqrPuNu9t++97FdYOLh1OPQDQDc0dgJ4zszcpND+thuxcMzuxrQGTme1hZrMVGqW/kHS2me1oZm+TNCs+7P1mNs3MXijpzHjb/vE6eJub2QmSXh1v/7CZ7ejutytcOuVbCmez/rtC0/W37t46KWM4lki6XeEs3fbv+VitORv3E2bWeWbuv5jZ/5nZ6ZLOl3SUu9/Yg3oA4BnM3ZuuAQBGhfhZuNe7++GDPOb9kh5x93MyvP5rJF0n6QB3v77X+wdQB86WBYBnmtB5g5mZu3u8gPLhkqZmeu1tMu0XQEVo7gBgjRskPcfMPibpfkn/pzCle62Z/UhhCnm+u/+hly9qZu+RtKnCp1tcK+mBXu4fQF2YlgWAQZjZREnXSNpO0qclzXX3p5qtCgAGRnMHAABQEM6WBQAAKAjNHQAAQEFo7gAAAApCcwcAAFAQmjsAAICC/H+kzu7F4emyDgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 792x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAncAAAGVCAYAAACcrtd2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAuJklEQVR4nO3de5xdVXn/8e8XLyggwXCpCVbRgIo1VWGIVavFn4bGCqmXVpFiG9RECzECaotIJEjUFpVYQTCDSLQVsdqqE7RIRVHwQkiiqHhBELRAqkBgsEgU5fn9sfaQk5PJZOacNdl77fm8X6/12jP77LPPc56zIc+syz6OCAEAAKAddqo7AAAAAORDcQcAANAiFHcAAAAtQnEHAADQIhR3AAAALUJxBwAA0CIPrjsAAMCOYXtPSS+W9FhJv5D05Yj4Ya1BAciOnjtggmw/0vbbbF9s+wbbP7R9j+3oaIfafobtszr2/dz25VX7mu31tn9dPbbfGK/3Uts3jbLfto+zPWT7vbY/avtvMr/Xh9h+te2LbF9r+zrbt3W912W2Z9o+zfbGat9w9R4vt/1V22s6nreg6/wn2/5elcNrbb/V9s6jxDK3eq/vt31e9XqT8gdqlduv2v6DcR5/iu1Pd/z+GdtLx/G83Wy/x/aPqvf/HdsLbY/5/2bbD7d9R9fn0NneO8r7WSzpJ5IOkPTOiPgghR3QTvTcARNg+7mS/lPSgyS9TtLqiLjX9kMkHSbpfEl/IEkRcZWkq6rn/LGkhRHxxa7z7S7p62O83kxJ50n61SgPnynpryTNjoi7bE+X9F3be0fE+8c45yMlHSzp6xFx7xjHPVrSkKSnSjpN0nERcYdtS3pK9V4Pqd7rrZJOtb2HpCWSzomIt3ad78GSBrte5qOSHinpIkn7SDpa0rskPU7Soo7nvkTSJyX9aUSsqWL4vKSLbL8iIn4/xvt4gaQfVDGO159Ieq6k4yS9fRzHP0nSjzp+f4KkT4z1BNsPUnoPv5D0MaXetKOVcrS3Uh625ZWSfi7pQ5Lu7Nj/MEmnK12jI69jSSskvVHSSyPiM+N4PwBKFhE0Gm0cTdKjJd0t6f8kzdrGMTMl/a+kQzv2fVVSdO7res7fS/rDUfZb0r9LulzSTV2PHVCd841d+0+Q9FtJM8Z4H4dWz91vjGMeLOnq6riXb+OYh0r6b0nLOvadVj1n2Tae87SR80n6M0nv73p8X0kbq3PsVe17kKSfSfpM17FPr4572XY+t5C0YIKf9SerPN4uaZdxHL9W0tEd8f5G0lO385y/k3R8176nVq97t6SdxnjuuaPFJeklkm7ufK6k11Y5OGOy/tug0WjNagzLAuNQ9X6cK+kRkt4XETeMdlyk3qHl3bvHOndEnBsR/zPKQ8cp9ezcNMpjI71aX+raf5mkh0g6dqzXHIc3ShqQdIWkT412QET8VqmY3GL3WCeNiO9ExL9Xvz5N0ildj98iaaRn6fHVdp6kx2jr93qNUiF4/FivOVG2H6NUyJ8raU9JC7ZzvLVlz91jlArfn2znpWZKOqtzR0Rco5TzR0jaa4znHhsRvx5l/5GSPhUR91ex7Sbp3ZLuq7YApgCKO2B8HiPp8Ornj23n2E9Kun48J7W9Yhv7D5T0jDFe63nV9kdd+0fmUD13PK8/hsXV9qMRMVbBdq2kz43nhLb/0vahI79HxL9ExP+Ncuidku6X9NPq95H3+oPOg6oC5keSnlUNi+eyWNI5SoVXSDqxGkLdgu0FtqOKdVdJV1e/j8R9z8hcSduPt32T7QeGWiPi3TH6cPKdkoaVeg1HNdpnUhVyR2jLYny+UpH4PUnvs31jNVfvEtuzt3X+ai7ku6r5e7+v3utDqscOtH1VdY49qn272D7D9oed5pL+l+0ndp3zBbY/Uc3RHKp+3qt67KG259v+lO2Tqp+vq+Zs2vaTba+y/SbbS53mcD5tW/EDUx3FHTA+f1xt71MaItymiLgtIm7e3glt76s0Z657/0MlnaE0ZLetwurRku7uLg4i4j5J9yjN3+qJ7WmS9qt+vW6sYyP59jjOuYukheMM4SlKQ7Ajxc2jq+1doxx7p9L/x/Yd57nHVBVIhyv1fl0v6WJJsyT95SiHX6Q0N+4lkn5Z/by3pJMlfa36+aDq2D2qGJ+0ndcfmc/44ZHetwk4Qqkn81sd+w6ttrdLOkmpN/QFSsP6a2w/XqOIiPsi4mRJX1bK75rq2lKkRRj3K80hvasq+r4m6WsR8Vql+YqPlnRlNW9Ttv9I0n9J+lJELJP0Ukl/Kun91Us+QekPkr+qYn6WpEuVekB3Uvoj5ysR8b6IOF2p1xbANlDcAeNzYLW9KyJ+1+M5PlythvyO7R9LulGbC5dOp0o6OyLuGONc0zX6Igsp9frMkCTbD7O9R2eTtFt13O7dj1WLHjp7XMaKYSyv73iv31daNPCi7T3J9v5K/9D/fcfu6dV2tPc7XG1n2H7wKO9nj+rxXUZ57GGjnO9vJX2yGnKWNhcfb+4+MCI2VQXoXpJ+EhG3V7/vLenH1e8bq2PXKy202aqY7/Jnkn4tabsrbUexxZBsZaTo/ceI+GVHMX680uKLt23nnCNTDJaM7Kh65O7smEqwUJIi4uJq+1ulBR17KS06klKR9r+qhqqr/4Z+qrSwRxHxfUkrq2Pvi4iTImJxRDyr+gPmiZIOrYpfSfonjV7sAxCrZYHxGpnftLftXSPinh7O8dqIuHzkl6rw+FrnAbafI2nX6FpVO4rbtblI67arNhdlJykVi6MZrffjedpyOHA/dQ2HjtOHqh4aSQ/0Rp471hOq3r1/k/SKiLit46GReEZ7v7tW2zuUeoK+so3Tf7BqnU6T1BnjTkpF5dyOY74i6fuSnmn7WRHxjVHOvb+2HIafJemb3QeNFHrbYnsfpWLyxTHGKuZtPHcPSS/U5p66ESM9u909yZcoDTkfpLFdrvReFtg+NSJ+obSid1XHMS+UNN125749lK6v3SSpKij/sIr1kZJepvSHTedw933Vdt0ocfyr0mfzWNtviIhrtxM3MKVR3AHjs7bj51mSvtvvCashrQd6aJxui3KKpFd29So9KD2c9kXEJkm3SDrY9k6dPTVVz9s0bZ7M/xFtvRDh6ZI+IOnlkjZ0PfY9SfcqrQjeTalw6VtE/Nb2+7T5H/AtVD0yH1BarHJ118O3VNs9RnnqSK/erUq9g88Z5ZgrlBYTfKFr/8+7fn+hUsFxyeYOoi1e481Kw4kjMX9XaS7mLpLC9ourhx4h6TDbJyktfLhwlJi2UH1uH5b0hoi4aXvHj+IlSu//W137R841XR1Fe0TcZ/t2bX5vo4qIsH26Uu7eUF2vL5T0zo7DZkn6bkQsGOtcth8n6fVK1+aFSkXifmM9p8Nipc/4NEnX2D5V0rt7GLoGpgSKO2B8rlaaf/YEScdo61WiPYmIz0mS7ccqzbU6TNseCh3pzbFSz8tAFU/noooDqu2XqvP/XF1FjDff+PfqbRUSti9SuoXG39n+4DYm/k9IRPygOvfuknbvmpf4j5Iui4jRVuZeolRYPUmpJ2kkRiu9/291LMy4cpT3IknXRcRWj3V5k6R5ke5P2Pn8PST9j6QX2z4gIkYK579Q+n/olUq9o5dJ2lnp8/gLpWH3bS6K6Dj/TpLeK+nMiLhie8dvw5GSPj1KsfNNSW9QmmvXPX/yoUq9kttziaT1Siuwr1DK96aOx2+TdGD3HxrSA5/1PUrTGi6X9PqI+HT12DheenOsEbHc9heV7u+4XCnX47kHITDlMOcOGIequHmD0kTyN9h++raOtf0824d37hrr3NU/7quU/uF85ijtC0rzlUZ+l6SPV9sXdJ1uZEhx1ZhvaPuWKS0SOEhbzn/rjv1xTt988MCucZ77gbmGtt8s6eaI2OKmv9Xqyn2VioJbtfV7na00j+2CcbzmmGw/S9LuktZ0PxYRdyn1qlnSiR37b1YqbGZKurwqlB8i6XeSroyImzpXA9ue7q5vnqgK1DMkDXUO2VePvWKkt9Zp9eru24h9H0nPV7onYrfPKs1LfHHXc2Yo9fB+drRzdqoW9bxL6WbTH9PW19YapSJ7iwUzVVznKA3/Lla6rUxnD/iuGn8Hw6lVLFcrTR24UptXrwPoQs8dME4Rcantl0s6W9I3bL9DaY7YzUr/qD9Z6ca0M6rtiJF/lB/RfU6nb5VYIWk4Iu7W1sNqsn2bpN9ExAOPRcSPbJ8t6c22/zUihqu5TG+WdFZE/LjP93qL7blKw7pn2R5QKkKuUypwH690/7lXVm0873VnpWG5V0l6S7XvREmvkPQ52533vNtXaRL986uhweMl/ZvtgYhYWxVFpysVCx/t571W3q7U87Wt1cmfUlqE8BrbKyJipBfsQKUe1Rur35+ktLhii+Fn2wdJukrptjF/Ve2zpH+pznFXVWBKqYh8nKSHR8Qnq33fkPRE2/tHxC+7YnuZ0vD6VV37FenbU94i6QO2L+jolTxRKXcru5+zDZ9T6gG+W1vPiXuvUi/vh2z/qdI1/ASlPzReHRH3V9emJL3N9kqlhSN7SZrp9A0iP1fqiZPSQo9ur7b9nojYGBG/s/2/GkevKDBVUdwBExAR/2H7v5X+cXyeUo/Eg5T+0bta6b5wl0qS7UOU5rU9rXr6x23/SOn2HTsr9YT8UfX8bfaOjeGNSnOqPmn7OqVi6P3avMKzLxHxXdvPULqJ7xFK/8Dvo9SLeJNSwXNoNX9rhqRXK82jkqQ32n6+0j/AVir6nqhU9H0yIn5v+3WS3lcdPzBKCAtHiq2I+JTtX0tabvsGpd6ymyW9KiJ+0+t7tL23pH+W9OeS9rF9Y0dBNXLMQUrFi5SK+P+y/Q8R8R9Kn9+1HcORT9Lmew12uktp7mDnza/fpdQbLG3dKyltubDjRqW5fZtGOW5bQ7KSpIg4z/adks6wfYtSj+z9kp473txVBdXnJd3QXQBHxC+q4n+FUsH/QqXh4L+JzbfJeYdSEXuk0tSBJUq3FFqpNBy+rDpGko6yvVHSBR0La/ZRunXLx5VWI9+h6g8EAFvztv9QBdBGTvexmy1pbdfcqVaqepN+Uq30RA+q4eS1SnMSu3sOATQMxR0AYFS2XQ2LH63US/va7T4JQO0YlgUAbMX2e5TmGK5V+oaW0YbOATQQq2UBAKMZ+ZqxR0l6WYzjK/UANAPDsgAAAC0yoWHZvfbaK/bbb79JCgUAAADjtW7dutsjYu/u/RMq7vbbbz+tXbt2+wdiUq1enbZHHFFvHCUjh3mQRwCoj+2fjbp/IsOyAwMDQXFXv5Fv7WFEvXfkMA/yCAD1sb0uIrZa7MRq2QIdzpfu9I0c5kEeAaB56LkDAAAo0LZ67rgVCgAAQItQ3AEAALQIxV2B7M0T2dEbcpgHeQSA5qG4AwAAaBFWyxaI2070jxzmQR4BoHnouQMAAGgRijsAAIAWobgr0BFH8HVP/SKHeZBHAGgebmJcIL7yqX/kMA/yCAD14evHWmRoqO4IykcO8yCPANA89NwBAAAUiJ477FBz50rDw3nONXu2dP75ec4FAEDbUdwVaHAwbRctqjeOsQwPS0uX5jnXyPvNqYQcloA8AkDzMCxboBImsc+Zk7e4W706z7lGlJDDEpBHAKgPw7ItsnBh3RGUjxzmQR4BoHko7go0GcOUUw05zIM8AkDzcBNjAACAFqG4K9Ctt6aG3pHDPMgjADQPw7IF2nfftGUSe+/IYR7kEQCah+KuQDNm1B1B+chhHuQRAJqH4q5ADIP1jxzmQR4BoHmYcwcAANAiFHcAAAAtQnFXoIMPTg29I4d5kEcAaB7m3BVo/fq6IygfOcyDPAJA81DcFWgyv9537lxpeLj/80yb1v85JhNfkZwHeQSA5qG4K9BkDoMND0tLl07e+ZuCocQ8yCMANA9z7gAAAFqE4q5Ay5alht6RwzzIIwA0j2MC3xs0MDAQa5lkUzs7bSfjK5/mzGnesOzgoLR6dd5zTmYOpxLyCAD1sb0uIga69zPnrkCnnlp3BOUjh3mQRwBoHoq7AjEM1j9ymAd5BIDmYc4dAABAi1DcFWjdutTQO3KYB3kEgOZhWLZAA9XUSSax944c5kEeAaB5KO4KdNBBdUdQPnKYB3kEgOahuCsQw2D9I4d5kEcAaB7m3AEAALQIxR0AAECLUNwVaObM1NA7cpgHeQSA5mHOXYE2bKg7gvKRwzzIIwA0D8VdgW65pe4IykcO8yCPANA8FHcFYhisf+QwD/IIAM3DnDsAAIAWobgr0KJFqaF35DAP8ggAzeOYwPcGDQwMxNq1aycxHIyHnbaT8ZVPc+ZIS5fmP28/Bgel1avznnMycziVkEcAqI/tdREx0L2fOXcFWrmy7gjKRw7zII8A0DwUdwViGKx/5DAP8ggAzcOcOwAAgBahuCvQ6tX556BNNeQwD/IIAM3DsGyB5s9PWyax944c5kEeAaB5KO4KdPjhdUdQPnKYB3kEgOahuCsQw2D9I4d5kEcAaB7m3AEAALQIxR0AAECLUNwVyN78zQDoDTnMgzwCQPNQ3AEAALQICyoKxG0n+kcO8yCPANA89NwBAAC0CMUdAABAi1DcFeiII1JD78hhHuQRAJqHOXcFuvjiuiMoHznMgzwCQPNQ3BVoaKjuCMpHDvMgjwDQPBR3BWIYrH/kMA/yCADNw5w7AACAFqG4K9DgYGroHTnMgzwCQPM4JnAX0oGBgVi7du0khoPxGPm6p8m4geycOdLSpfnP24/BQWn16rznnMwcTiXkEQDqY3tdRAx072fOXYEWLqw7gvKRwzzIIwA0D8VdgRgG6x85zIM8AkDzMOcOAACgRSjuCnTrramhd+QwD/IIAM3DsGyB9t03bZnE3jtymAd5BIDmobgr0IwZdUdQPnKYB3kEgOahuCsQw2D9I4d5kEcAaB7m3AEAALQIxR0AAECLUNwV6OCDU0PvyGEe5BEAmoc5dwVav77uCMpHDvMgjwDQPBR3BeLrfftHDvMgjwDQPBR3BWIYrH/kMA/yCADNw5w7AACAFqG4K9CyZamhd+QwD/IIAM3jmMD3Bg0MDMRaJtnUzk7byfjKpzlzpKVL85+3H4OD0urVec85mTmcSsgjANTH9rqIGOjez5y7Ap16at0RlI8c5kEeAaB5KO4KxDBY/8hhHuQRAJqHOXcAAAAtQnFXoHXrUkPvyGEe5BEAmodh2QINVFMnmcTeO3KYB3kEgOahuCvQQQfVHUH5yGEe5BEAmofirkAMg/WPHOZBHgGgeZhzBwAA0CIUdwAAAC1CcVegmTNTQ+/IYR7kEQCahzl3Bdqwoe4IykcO8yCPANA8FHcFuuWWuiMoHznMgzwCQPNQ3BWIYbD+kcM8yCMANA9z7gAAAFqE4q5Aixalht6RwzzIIwA0j2MC3xs0MDAQa9euncRwMB522k7GVz7NmSMtXZr/vP0YHJRWr857zsnM4VRCHgGgPrbXRcRA937m3BVo5cq6IygfOcyDPAJA81DcFYhhsP6RwzzIIwA0D3PuAAAAWoTirkCrV+efgzbVkMM8yCMANA/DsgWaPz9tmcTeO3KYB3kEgOahuCvQ4YfXHUH5yGEe5BEAmofirkAMg/WPHOZBHgGgeZhzBwAA0CIUdwAAAC1CcVcge/M3A6A35DAP8ggAzUNxBwAA0CIsqCgQt53oHznMgzwCQPPQcwcAANAiFHcAAAAtQnFXoCOOSA29I4d5kEcAaB7m3BXo4ovrjqB85DAP8ggAzUNxV6ChobojKB85zIM8AkDzUNwVaKoNg23aJM2Z0/95Zs+Wzj8//TzVcjhZyCMANA/FHRpvyZI85xkczHMeAACajAUVBRocpFDpFznMgzwCQPM4JnAX0oGBgVi7du0khoPxGPm6p8m4geycOdLSpfnP2wSDg9Lq1ennyczhVEIeAaA+ttdFxED3foZlC7RwYd0RlI8c5kEeAaB5KO4KxDBY/8hhHuQRAJqHOXcAAAAtQnFXoFtvTQ29I4d5kEcAaB6GZQu0775pyyT23pHDPMgjADQPxV2BZsyoO4LykcM8yCMANA/FXYEYBusfOcyDPAJA8zDnDgAAoEUo7gAAAFqE4q5ABx+cGnpHDvMgjwDQPMy5K9D69XVHUL5+czh3rjQ8nCeW2bOl88/Pc64djWsRAJqH4q5AfL1v//rN4fBwvu/gLflbHrgWAaB5KO4KxDBY/8hhHuQRAJqHOXcAAAAtQnFXoGXLUkPvyGEe5BEAmofirkCnnZYaekcO8yCPANA8zLkr0Kmn1h1B+chhHuQRAJqH4q5ADIP1jxzmQR4BoHkYlgUAAGgReu4KtG5d2o7chiLnDXWnTctznibatEmaMyf9fM89abvrrr2dq815mojuaxEAUD+KuwINDKRtRNrmvKFumy1Zsvnn+fPTdmionljaovtaBADUj+KuQAcdVHcE5Zs1q+4I2oFrEQCah+KuQCNDYejdihV1R9AOXIsA0DwsqAAAAGgRijsAAIAWobgr0MyZqaF3Cxakhv5wLQJA8zDnrkAbNtQdQfk2bqw7gnbgWgSA5qG4K9Att9QdQfkuuKDuCNqBaxEAmofirkAMg/Vvzz3rjqAduBYBoHmYcwcAANAiFHcFWrQoNfTu7LNTQ3+4FgGgeSjuCnTeeamhd5demhr6w7UIAM3DnLsCrVxZdwTlO/bYuiNoB65FAGgeirsCMQzWv3nz6o6gHbgWAaB5GJYFAABoEYq7Aq1enRp6t2ZNaugP1yIANA/DsgWaPz9tI+qNo2TLl6ft0FC9cZSOaxEAmofirkCHH153BOU75JC6I2gHrkUAaB6KuwIxDNa/pUvrjqAduBYBoHmYcwcAANAiFHcAAAAtQnFXIDs19G7+/M2LAdA7rkUAaB6KOwAAgBZhQUWBuO1E/7gFSh5ciwDQPPTcAQAAtAjFHQAAQItQ3BXoiCNSQ+9OPz019IdrEQCahzl3Bbr44rojKN/VV9cdQTtwLQJA81DcFYjFAP075ZS6I2gHrkUAaB6KuwIxDNa/OXPqjqAduBYBoHmYcwcAANAiFHcFGhxMDb275JLU0B+uRQBoHoZlC/S616XtokX1xlGyc85J23nz6o2jdFyLANA8FHcFWriw7gjKd9hhdUfQDlyLANA8FHcFYhisf4sX1x1BO3AtAkDzMOcOAACgRSjuCnTrramhd3fckRr6w7UIAM3DsGyB9t03bSPqjaNkxxyTttyEtz9ciwDQPBR3BZoxo+4Iyjd9et0RtAPXIgA0D8VdgRgG69+qVXVH0A5ciwDQPMy5AwAAaBGKOwAAgBahuCvQwQenht6dcEJq6A/XIgA0D3PuCrR+fd0RlO+GG+qOoB24FgGgeSjuCrR2bd0RlO/MM+uOoB24FgGgeSjuCsQwWP/237/uCNqBaxEAmoc5dwAAAC1CcVegZctSQ+8uvDA19IdrEQCah+KuQKedlhp6d9FFqaE/XIsA0DzMuavJ3LnS8HBvz505M23nzEnbadPyxDSVHHlk3RG0w6mn1h0BAKAbxV1NhoelpUvrjmLqOuqouiNoB4ZkAaB5GJYFAABoEYq7Al1/fWroHTnMY9261AAAzcGwbIFOPDFth4bqjaNk5DCPgYG0jag3DgDAZhR3BZo1q+4IykcO8zjooLojAAB0o7gr0IoVdUdQPnKYB0OyANA8zLkDAABoEYo7AACAFqG4K9CCBamhd+Qwj5kzN99UGwDQDMy5K9DGjXVHUD5ymMeGDXVHAADoRnFXoAsuqDuC8pHDPG65pe4IAADdKO4KtOeedUdQPnKYB0OyANA8zLkDAABoEYq7Ap19dmroHTnMY9Gi1AAAzUFxV6BLL00NvSOHeZx3XmoAgOZgzl2Bjj227gjKRw7zWLmy7ggAAN0o7go0b17dEZSPHObBkCwANA/DsgAAAC1CcVegNWtSQ+/IYR6rV6cGAGgOhmULtHx52g4N1RtHychhHvPnp21EvXEAADajuCvQIYfUHUH5yGEehx9edwQAgG4UdwVaurTuCMpHDvNgSBYAmoc5dwAAAC1CcQcAANAiFHcFmj9/80R29IYc5mGnBgBoDoo7AACAFmFBxQTNnSsND/d/nmnTen8ut+/oX5NyuGmTNGdO3VFsafZs6fzzt38ct0ABgOahuJug4WFWWiKvJUvqjmBrg4N1RwAA6BXDsgAAAC1CcVeg009PDb0jh3kccURqAIDmYFi2QFdfXXcE5SOHeVx8cd0RAAC6UdwV6JRT6o6gfORwbONd5HHAAWm7rWPHuzADAJAPxV2BmrayskTkcGy5FnmwMAMAdjzm3AEAALQIxV2BLrkkNfSOHOZBHgGgeRiWLdA556TtvHn1xlEycpgHeQSA5qG4K9Bhh9UdQfnIYR7kEQCah+KuQIsX1x1B+chhHuQRAJqHOXcAAAAtQnFXoDvuSA29I4d5kEcAaB6GZQt0zDFpOzRUbxwlI4d5kEcAaB6KuwJNn153BOUjh3mQRwBoHoq7Aq1aVXcE5SOHeZBHAGge5twBAAC0CMUdAABAi1DcFeiEE1JD78hhHuQRAJqHOXcFuuGGuiMoHznMgzwCQPNQ3BXozDPrjqB85DAP8ggAzUNxV6D99687gvKRwzzIIwA0D3PuAAAAWoTirkAXXpgaekcO8yCPANA8FHcFuuii1NA7cpgHeQSA5mHOXYGOPLLuCMpHDvMgjwDQPBR3BTrqqLojKB85zIM8AkDzMCwLAADQIhR3Bbr++tTQO3KYB3kEgOZhWLZAJ56YtkND9cZRMnKYB3kEgOahuCvQrFl1R1A+cpgHeQSA5qG4K9CKFXVHUD5ymEepeZw7Vxoe7v88s2dL55/f/3kAICeKOwBTzvCwtHRp/+cZHOz/HACQGwsqAAAAWoTirkALFqSG3pHDPMgjADQPw7IF2rix7gjKRw7zII8A0DwUdwW64IK6IygfOcyDPAJA81DcFWjPPeuOoHzkMA/yCADNw5w7AACAFqG4K9DZZ6eG3pHDPMgjADQPxV2BLr00NfSOHOZBHgGgeZhzV6Bjj607gvKRwzzIIwA0D8VdgebNqzuC8pHDPMgjADQPw7IAAAAtQnFXoDVrUkPvyGEe5BEAmodh2QItX562Q0P1xlEycpjH9vK4aZM0Z06e15o9Wzr//DznAoA2o7gr0CGH1B1B+chhHtvL45Il+V5rcDDfuQCgzSjuCrR0ad0RlI8c5kEeAaB5mHMHAADQIhR3AAAALUJxV6D581ND78hhHuQRAJqH4g4AAKBFWFBRIG7f0T9ymAd5BIDmoecOAACgRSjuAAAAWoTirkCnn54aekcO8yCPANA8zLkr0NVX1x1B+chhHuQRAJqH4q5Ap5xSdwTlI4d5kEcAaB6KuwLl+iL2qYwc5kEeAaB5mHMHAADQIhR3BbrkktTQO3KYB3kEgOZhWLZA55yTtvPm1RtHychhHuQRAJqH4q5Ahx1WdwTlI4d5kEcAaJ4pUdzNnSsND+c517Rpec7Tj8WL646gfOQwD/IIAM0zJYq74WFp6dK6owAAAJh8LKgo0B13pIbekcM8yCMANM+U6Llrm2OOSduhoXrjKBk5zIM8AkDzUNwVaPr0uiMoHznMgzwCQPNQ3BVo1aq6IygfOcxjR+Zx06Z834jRhIVRADBZKO4AFGHJkrojAIAysKACAACgRSjuCnTCCamhd+QwD/IIAM3DsGyBbrih7gjKRw7zII8A0DwUdwU688y6IygfOcyDPAJA81DcFWj//euOoHzkMA/yCADNw5w7AACAFqG4K9CFF6aG3pHDPMgjADQPxV2BLrooNfSOHOZBHgGgeZhzV6Ajj6w7gvKRwzzIIwA0D8VdgY46qu4IykcO8yCPANA8DMsCAAC0CMVdga6/PjX0jhzmQR4BoHkYli3QiSem7dBQvXGUjBzmMdXzuGmTNGdOnnPNni2df36ecwGY2ijuCjRrVt0RlI8c5jHV87hkSb5zDQ7mOxeAqY3irkArVtQdQfnIYR7kEQCahzl3AAAALULPHQCgGHPnSsPDec7FPEe0FcVdgRYsSNtVq+qMomzkMA/yiB1teFhaujTPuZjniLaiuCvQxo11R1A+cpgHeQSA5qG4K9AFF9QdQfnIYR7kEQCah+KuQHvuWXcE5SOHeZBHAGgeijsAaIBcN0TOuUiAxQvjQ57QNBR3BTr77LRdvLjeOEpGDvMgj/nkuiFyzkUCLF4YH/KEpuE+dwW69NLU0DtymAd5BIDmoeeuQMceW3cE5SOHeZBHAGieRhd3ueYxTJvW/zmaZN68uiMoHznMgzxivHLNKWzb/8+BydDo4i7nPAYAQH1yzSkEsH3MuSvQmjWpoXfkMA/yCADN0+ieO4xu+fK0HRqqN46SkcM8yCMANA/FXYEOOaTuCMpHDvMgjwDQPBR3BWIeYv/IYR7kESVr4iKPXDFJ+W6I3PabNLfx/VHcAQCmpCYu8sgZU64bIrf9Js1tfH8sqAAAAGgRirsCzZ+fGnpHDvMgjwDQPBR3AAAALeKIGP/B9m2SfjZ54WAUe0m6ve4gkAWfZTvwObYDn2M7TPXP8bERsXf3zgkVd9jxbK+NiIG640D/+Czbgc+xHfgc24HPcXQMywIAALQIxR0AAECLUNw1X0PumoMM+Czbgc+xHfgc24HPcRTMuQMAAGgReu4AAABahK8fAwBMWbYfLOkpEfGdumPBxNjeSdKfS5oj6deSrouIz9UbVTPQc9dQtl9q+1u277Z9k+3zbO9Td1yYGNvPsf0V28O2N9j+iO296o4LvbP9UNvrbC+oOxZMnO2v2o6RJuk+SU+pOy5MjO3ZktZLerukCyPiPRR2m9Fz10C2F0laIunTkr4s6cWSXivpT2wfEhGbagwP42T7+ZLeIukMSXdKmi/prZKmK32mKNMySQfVHQQmzvazJe0s6X0du4cl/Xs9EaEXtp8r6fOSLpX0ioj4Xc0hNQ7FXcPY3lXSyyUNjBRxtk+V9BVJz5b0Ikn/UV+EGI9quOBISS+JiHur3d+y/SxJh9YWGPpi+08lHVh3HOjZSZJeExHX1h0IemN7hlLHxy8kvZrCbnQMyzbPMyW9rbN3LiLuk7Sy+nVWLVFhoizp+I7CbsROkq6sIR70yfbuSj2vJ9cdCybO9tMkHSbpZNsLbO9bc0jozVJJe0s6MyKG6w6mqSjuGiYivhQRV43y0J3V9ic7Mh70JiJ+HxH3dO6zfaCkh0haWE9U6NN7lf5h6S7YUYb5ku6XdJSkCyTdbPsC2w+rNyyMl+2dJf1t9eujbK+3fY/tb9t+vW3XGV+TUNyV4ymSbpX0hboDwcTY3tn2UZK+IekHkpgzWRjbL5N0Y0SsrzsW9CYi3iFpV6XRj7cqzbVbIGlVfVFhgg5W+gx/pTTf7mBJj1ZaWHGu0rA7RHFXhOqvlaOVJo7+pu54MGHHSXqhpI2SXi3pmmqIDwWwPVPSq5QWxqBgEXF/RPw0Iv5J0hMl3SjpFbafWXNoGJ+RofSPRMSVkdwp6e8l/VJpyH2X+sJrDoq7MnxA0sqIYK5WgSLizIh4laQnKH1Vzh8q9Rig4aqFMf8i6YSI+H3d8SCfiPiFNg/xPaPOWDBuI/8N3ty5MyJ+q9STt5ukA3Z0UE1Ecddw1b20Nkk6u+ZQ0KeqOHiz0v+gnlhzOBif4yRdJmmD7YdV87N2rh57SLWPuw6U6ypJv5H00LoDwbjcVG2nj/LYhjEem3Io7hrM9guU5todH3wJcCtExK+UFsf8T92xYFxepjSX596O9qPqscHq91PqCQ0Z3K/0xxYL1crwQ0l3S3r8KI+NFOg37Lhwmou/OBvK9vMk/YWkN3UWdtXy/SdFxGW1BYee2X6MpGlK92lC8x0rqXt+5AxJ/ylpudKNVG/ufhKK8f8k3SXpizXHgXGIiHttXyTpSNs7d81BP1DSt8UfzpIk0yHUPNXdtz8m6SNKf1mOmKZU8B0aEbfVERvGx/Yekj6lNA/kQxHxK9uPkHSRpKGIWDnW89FctvdTmoh/TESsqjcajIftP5H0OUmXSzoxIm6xvb+kT0h6Y0R8o874MH5VB8e3Jf1rRLyp2vdUSV+TdHhEXFFnfE1Bz13D2B5Qut3JrpJOG+WQj1PYFeG3Ssv1T5S0yPanlQr1kyPimlojA6aen0q6Qukm8ettXyrp55LmR8SGMZ+JRqkK82dI+ifbX1S6vdR0Sc+JiO/WG11z0HMHAADQIiyoAAAAaBGKOwAAgBahuAMAAGgRijsAAIAWobgDAABoEYo7AACAFqG4AwAAaBGKOwBZ2H6/7bB9j+3rbX/f9n3Vvuts/9j2xur371TPmV4de1bN4UuSbN9k+1fVe1lm+8k76HVPtr3c9veq/By6I14XQDvxDRUAcrpA0nERca+UiiVJj5V0WETcZHsnSa+rmiRZ6Y9M1xDrtnwhIo7fkS8YEe+SJNvflHTxjnxtAO1DcQcgl5D01pHCbtQDIu6XdG71/cmKiDskPX4HxTde24x/B7inxtcG0BIMywLI5byI+MU4j33HpEYCAFMYxR2ALCLiBxM49oeSZHtn2y+y/a6Rx2w/xfbptq+1/TTbb7H9E9u/tP3P1THPsn2J7Tttr7G9f+f5be9p+2zbF9i+xvanbT+ql/dl+/m2N1Rz4T5ie1a1/+G2/7maR/hs24+yfYLt9bYfb/usKubXV8cvrGJabPvDtod6iQcAtodhWQC1qAqyxZLeKOlnkk6uHrpH0v6SnizpFEkfk/RMScdLepvtvSR9R9Ixkh4t6QpJ75b019V595D0LUmviIj1th8h6aeSrrT91IiY0NBnRFxm+x+qOG6PiBuq/ffa/pmk90TE122/SNLRkp4u6S2SrpN0raTdbM+WtELSPhHx6yqmz00kDgAYL3ruANQiIq6vFi7c3rX/RkmXVb9+KCKGIuJ2peJKknaJiLMiYkNEXC1pvaTndJziZEnXRMT66ny/kvQFSbMkHdljuJ+QdKOk19jetWP/KyV9tHqdz0u6pNr/pYhYERHPi4j3SnqCpF0lzemIaXmPsQDAmCjuANRttJ6033VtJen6avubrmOvl7R3x+8vlPQE26tGWvX4NZKm9xJgRPxOqXdwuqTXSA/0PN4dEbd2HHpftV3XdYovS9og6Yu232l714j4ci+xAMD2MCwLoAgRcb896h1T7teWf6jOkjQ4Cbcz+Zikt0s60fa5kv5G0qrxPDEi7rQ9oNQDeLKko22/PCKuyhwjANBzB6B1blOar7cV24/s9aQR8RtJZyjdt++vlXoIx7UowvbDqh6+50t6vaTdJF1h+8Be4wGAbaG4AzCZRv4fsyNvUrxG0lzbf9650/YBkt7Z57k/rFQ8rpC0tir4xuNI20+MiN9FxEpJz1YaOZnbZzwAsBWKOwCTovo2imnVr9PGOPThkvaujh+xZ7V9WMf5dq9+nNH1/EdWj+9c/X6a0ty3i22fa/tY2yslfUbSmRN+Ix2qGzSvlLSP0rdxdHtEd9wd3tzx821Kcwev7iceABgNxR2A7Gy/TKn4GSnIzrZ9fEcBJtt/ZPsUpUJpF0nvt/1Y20dLWlQddpLtubafJOl91b7nVPfBe7jtJZL+rNr/LtuPjYjvK9065WtKq1nfrlR0/WVEjCzK6MdqSd9XWqXb+Z6P0+bVuO+13b0y97W2P2t7maRzJS2MiG9miAcAtuCIqDsGAGiE6rtwL4+IBWMcc5KkeyLirEl4/UMlfUXS8yLi8tznBzA1sFoWALY0s3uHbUdEVDdQXiBpYJJee49JOi+AKYTiDgA2+7qkfWyfIWmjpM8qDeleZvtapSHkMyPi/3K+qO03SdpZ6dstLpN0Z87zA5haGJYFgDHY3lfSf0vaS9IHJZ0eEffXGxUAbBvFHQAAQIuwWhYAAKBFKO4AAABahOIOAACgRSjuAAAAWoTiDgAAoEX+P78pgyVbw8SIAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 792x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAncAAAGVCAYAAACcrtd2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAo4klEQVR4nO3dfZheVXnv8d+NImDEYAJWgi/UREXqVAvJWPVQX4NphbQ9aE1zrE1qiTSmkQS1ihkNJljfmlSJWAbTpLZNc07bo51QxVR8b8GQRFQUpUkFC+SoJDgiElG4zx9rjzwOk2SSWbPXrPv5fq5rXTuzn/3suZ/8xniz19p7zN0FAACAGI4qXQAAAADyobkDAAAIhOYOAAAgEJo7AACAQGjuAAAAAqG5AwAACOThpQsAgNzMzCSd1Yz7Jd3m7n9XtioAaAdX7gDIzM4ys8+Y2aCZ7TGzvzazE0fxvt80szWj/B4PN7OfmNkzmq+faWb3mdkh/yPTzM4ws41m9i4zu8zMPmFmv3WAY58o6TOSrpC01d3fdbiNnSUDZrbycN4HABMBV+6ALmdmL5b0RknvkXSXpLmS3iJpiqTfOcTbL5T0XDN7h7v/4BDHnqr0b86u5uunStrl7j87RH3PkPRZSc919xubfU+VdIOZvdTdv9Bx7NMlfVrSrZLOcPd7D1HTgZwv6VxJO4/w/QBQDM0d0MXM7ChJ8yT9bkcjdJ2ZPVfSCw7x3tObYx4haZFSc3gwp0n6trvvb75+qqRvjaLM+ZLuG2rsJMndbzazb0k6T9IXmnomSfpYU88rjrSxM7On6BCfHQAmMqZlge5mki4coRE6StIXD/HeCyVdIOmnkl5vZo84xPGnSfpmx9ejbe6OkzTVzH7t50WbHSPpKZL2dBz3uuac6939v0dx3ocws6MlvU/pswFAlWjugC7m7ve7+z2d+5qpzaOVpiZH1KzHe7akv5G0WdI0pSuAIx17i5m5pPdKepmZefP1qyX9WfP1gubYdzXH/3LHKTZI+rGkT5jZi5qrjZdK+rykyzuO++MHv6Vda2Z3m9nXzWzFaNb1NS6W9GF3/95oDjazmWZ2U/MZtpjZr3a8tqxZUzi3Y98LzOxvmvV8Xzazi5ubP4Zen2Rma83sA832WjP7vY7Xf7n5PN8ys0ea2SYz+56ZzW3WCfaZ2XvN7HXNa/2j/NwAAqG5AyApXQ0zs/mS/kPSNyTtP8jhi5SaoAckvb/Z94bORqXDGZJOknStpDc0fz5J0g8kndP8eXNz7FMkPV5pvZ8kyd2/Kmm2pGMlfUrS9ZL+S9LL3P3upvZpzXvvb77PcyWdLOmfJa2S9MFRfP5nS3qsu2851LEdtW2X9Nrmy3ubWod8S9IGdx9ozv9ySSslXeDucyW9W6lJvbTjPVcqTTW/3t2XSdoo6R/M7LTm9WcrNbFPbc61TdK3JU2S9DJJCyW9yd0/qNSoThGA7uPuDAaDIUnLJf2tpN2SXNJ3JD16hOMeIemmzteUrqK5pLMPcv7bJL2k+fPk5vhpw445StKUA7x/saSfNO+7SdKvdrw2q9n/sRHet0Op6Tv1ILU9StLHJT2qY59LWjnKv7vPKU1PP7Fj3yZJvc2fj5H035LO6Xh9klIDvV/SpGbfRzo/g6QXNnW8smPf3zf7zhxWw0VKVzhP79j3otI/VwwGo/3BlTsAkiR3X+Puf6B0Vahf0hMkLRjh0FcoPWLkhx37fn71bqRzm9kjJZ2iB++UfbJSU/P/htXwgLvvG+H9f6R0le8UpUecnCbpc2bW0xxyf7O9bYRv/3GlpvGZI9XWeI+kS9z9Rwc55mBWK92gdmFT7/FKVxKvb15/htIVyT9qHumyUelq4g1K6xAfL0nu/mp3/51mivUsSS9v3n9cx/f6aXPsjmE1/IvS38N1ZrbczI52908f4ecBUDHulgXwC9z9fjN7g6TXSHpa52vNtOvrJT3azG7oeOlhkn4mabaZPdPdv9IcP19pXdzQf0h+pVlvd7TSNOs+M/uOu/+qDsDMXqU0XXmiu98l6QIz+7qkD0haozRle0tz+EjTkHsO8prM7FylKeKvmNmxw15+eLPvfnf/6YFqVJou3iZpkZmtUnqczP9xd29en95s3+HuNxzoJM3f76uUGtGPSPpHpSuWh+Tuu8xsVvOev5C0wMxe7u43j+b9AOLgyh2Ah/C0lu0upanETmdJ2uvup7n7szpGj6R3Ncdc1HH8gKRnSXqHpBuVmpZnKd2IMfTaiA8j7vBKSV9vGruh+i6TdJ3S2jo1td6sdEVwuKG7eHcf4PznKT3X795hQ5Le2vz5yoMV2DRxq5WmWv9E6eaSzgcnf7/Znj7S+83sMc0fL1f6+3uL/+L6vUMys2Pd/ZuSepvPM0PS9o5zA+gSNHcAHqL5LQ+TJf3TsJdWSFp3gLddpjRl+PtmNjTN+CN3v0XSiZJucPdbmq8fK2ln8/XPp1LN7CgzG36F7U5Jx4/w/b6kZhq2aa42Suptbq7o9HSlq3dDU6QysykdN3+skvScEYYkrW/+vOoAn7nTVZK+prR28T5373xMy9eU1gu+zcwmd77JzF4n6QVN3Rco/T0NXSWc1GxHM8uyzMwe7e73uvu7lB5AfbxSQw6gi9DcAV3MzE4ws38zszc268SG1ot9SNKfuvuujmN/XdLzJH1ypHN5enzIZ5UakdXDXn6GUoMz5DSlmyKG+0dJ3+18pp2kv5I0bdgjQSZJeqlSQznkMkn/Ken9zeNSZGZPkPS76niWXzMNe6ekP2/q3u3u1w0fzTlva74+0FW/zs/vSmsPpypNqXa+dqfSVOnTJH25eWTJRWZ2laTnKz18eajpO8fMXmpmc5Qe4CxJLzSzlzR/Hspp+BTyMZL+tOPr25UaysO6Agigfqy5A7rbfZLuVrratMjM/knSA5IuHlo3J0lmdp7S+rZHSrrUzNb5sAcFm9mFSjcRSNIfmtnDJC1x90FJv6LmmXTNM+eeopGbu5sl3SFpcGiHu3/JzJ4vqc/MXiZpr9INCJcq3Tk6dNyPzOx5Sk3bNWb2NaWG6ZXu/pmO73GX0o0cR/Sg40PYImmf0lW84fqUpmf/SGna9L+VnuG3pmkMbzKzdyutsfsrpQZxoaQnSnqx0nP+LpY01OR92Mz+3t0/0fE9VjWPTfm20t/5bzdXSgF0EXtwvS8AYCzMbJ6kWe5+0SEPBoBxQnMHAGNgZububmbHKT1A+Xe4WgagJKZlAeAINTeOfNXMblV6HMwnaewAlMYNFQBw5PYprRN8ktJv6VhRthwAYFoWAAAglMOalj3xxBP91FNPHadSAAAAMFo7duy4091PGr7/sJq7U089Vdu3b89XFcZsy5a0PffcsnUgHzKNiVwB5Nas933o/sOZlp05c6bT3E0sQ8/YZ3Y9DjKNiVwB5GZmO9x95vD93C1buXPOKV0BciPTmMgVQFu4cgcAAFChA12541EoAAAAgdDcAQAABEJzVzmzBxdqIwYyjYlcAbSF5g4AACAQ7patHI9ViIdMYyJXAG3hyh0AAEAgNHcAAACB0NxV7txz+XVG0ZBpTOQKoC08xLhy/EqjeMg0JnIFkBu/fiyogYHSFSA3Mo2JXAG0hSt3AAAAFeLKHUZl9mxpcHDs5+npkdavH/t5AADA4aG5q1x/f9ouWpTnfIODUl/f2M8zVBcOX+5MMTGQK4C2MC1budyLtHt78zV3W7aM/TzdiIX3MZErgNyYlg3q/PNLV4DcyDQmcgXQFpq7yjH9GQ+ZxkSuANrCQ4wBAAACobmr3B13pIE4yDQmcgXQFqZlK3fKKWnLIu04yDQmcgXQFpq7yp18cukKkBuZxkSuANpCc1c5pnniIdOYyBVAW1hzBwAAEAjNHQAAQCA0d5U788w0EAeZxkSuANrCmrvK7dxZugLkRqYxkSuAttDcVY5f9RsPmcZErgDaQnNXOaZ54iHTmMgVQFtYcwcAABAIzV3lVq5MA3GQaUzkCqAt5ofxu3Bmzpzp21k4MqGYpW2uX2nU2yv19Y39PP390pYtYz9PN8qdKSYGcgWQm5ntcPeZw/ez5q5yb3976QqQG5nGRK4A2kJzVzmmeeIh05jIFUBbWHMHAAAQCM1d5XbsSANxkGlM5AqgLUzLVm5ms4ySRdpxkGlM5AqgLTR3lTvjjNIVIDcyjYlcAbSF5q5yTPPEQ6YxkSuAttDcHabZs6XBwbGfp6dHWr9+7OcBAADoRHN3mAYH8z3kFwAAIDfulq3ctGlpIA4yjYlcAbSFK3eV27OndAXIjUxjIlcAbaG5q9ztt5euALmRaUzkCqAtNHeVY5onHjKNiVwBtIU1dwAAAIHQ3FVu0aI0EAeZxkSuANpCc1e5K69MA3GQaUzkCqAtrLmr3BVXlK4AuZFpTOQKoC00d5VjmiceMo2JXAG0hWlZAACAQGjuKrdlSxqIg0xjIlcAbWFatnJz56ate9k6kA+ZxkSuANpCc1e5c84pXQFyI9OYyBVAW2juKsc0TzxkGhO5AmgLa+4AAAACobkDAAAIhOaucmZpIA4yjYlcAbSF5g4AACAQbqioHI9ViIdMYyJXAG3hyh0AAEAgNHcAAACB0NxV7txz00AcZBoTuQJoC2vuKnfVVaUrQG5kGhO5AmgLzV3lBgZKV4DcyDQmcgXQFpq7Qvbvl3p785yrpyf2dM/s2dLg4NjP09MjrV8/9vOMt8hZdjNyBdAWmrtCli7Nd67+/nznmogGB6W+vrGfJ/rfEwAAEjdUVO/qq6Vbby1dBXLq76cRjYhcAbSFK3eVu/zy0hUgt9e+Nm0XLSpbB/IiVwBtobmr3NlnS9/8ZukqkNP555euAOOBXAG0heauckuWMNUTDXnGRK4A2sKaOwAAgEBo7iq3d296rAriuOOONBALuQJoC9OylVu4MG1zPTNv8uQ858n5HL9cNeWU69l70kOfv3fKKWnrnuf8mBjIFUBbaO4qN2VK2uZ4DlxOOZ/jNxHlevae9NC1WCefnOe8mFjIFUBbaO4qt3Fj6QqQG1N3MZErgLaw5g4AACAQmjsAAIBAaO4qt2xZGojjzDPTQCzkCqAtrLmr3O7dpStAbjt3lq4A44FcAbSF5q5ya9aUrgC5bd9eugKMB3IF0Baau8rNmFG6AuTG1F1M5AqgLTR3QGG5Hvg8/GHIAIDuRHNXuU2b0nb+/LJ14MgNf+DzkWbKL6af2Fau/MUtAIwX7pat3ObNaSAOMo3pkkvSAIDxxpW7ys2bV7oC5EamMb397aUrANAtaO4qx3RsPGQaE9OxANrCtCwAAEAgNHeV27UrDcRBpjHt2JEGAIw3pmUrt3x52g4MlK0D+ZBpTDNnpq172ToAxEdzV7np00tXgNzINKYzzihdAYBuQXNXubVrS1dQj1wPC5akyZPznGckZBoTU7IA2kJzh64x/GHBAABExA0VAAAAgdDcVW7BgjQQB5nGNG1aGgAw3piWrdy+faUrQG5kGtOePaUrANAtaO4qt2FD6QqQG5nGdPvtpSsA0C1o7io3dWrpCpAbmcbElCyAtrDmDgAAIBCau8qtW5cG4iDTmBYtSgMAxhvNXeW2bk0DcZBpTFdemQYAjDfW3FVu8eLSFSA3Mo3piitKVwCgW9DcVW7OnNIVIDcyjYkpWQBtYVoWAAAgEJq7ym3blgbiINOYtmxJAwDGG9OylVu9Om0HBsrWgXzINKa5c9PWvWwdAOKjuavcrFmlK0BuZBrTOeeUrgBAt6C5q1xfX+kKkBuZxsSULIC2sOYOAAAgEJo7AACAQGjuKjd37oMLtREDmcZklgYAjDeaOwAAgEC4oaJyPC4jHjKNiUegAGgLV+4AAAACobkDAAAIhOaucqtWpYE4yDSmc89NAwDGG2vuKnf99aUrQG5kGtNVV5WuAEC3oLmr3IoVpStAbmQaEzfKAGgLzV3lentLV4DcyDQmpmQBtIU1dwAAAIHQ3FXu6qvTQBxkGlN/fxoAMN6Ylq3c5Zen7Zw5ZetAPmQa02tfm7aLFpWtA0B8NHeVO/vs0hUgNzKN6fzzS1cAoFvQ3FVuyZLSFSA3Mo2JKVkAbWHNHQAAQCBcuavc3r1pO3Vq2TqQz5Fmun9/vseo9PRI69fnOVcus2dLg4N5zlXi891xR9pOm9bu9wXQfWjuKrdwYdrygNQ4jjTTpUvz1TARpxAHB6W+vjznKvH5Tjklbd3b/94AugvNXeWmTCldAXIj05hOPrl0BQC6Bc1d5TZuLF0BciPTmIamZQFgvHFDBQAAQCA0dwAAAIHQ3FVu2bI0EAeZxnTmmWkAwHhjzV3ldu8uXQFyI9OYdu4sXQGAbkFzV7k1a0pXgNzINKbt20tXAKBbdEVzl/Php5Mn5zlPLjNmlK4AuZFpTEzJAmhLVzR3OR9+CgAAMJFxQ0XlNm1KA3GQaUwrV6YBAOON5q5ymzengTjINKZLLkkDAMZbV0zLRjZvXukKkBuZxvT2t5euAEC3oLmr3Pz5pStAbmQaE1OyANrCtCwAAEAgNHeV27UrDcRBpjHt2JEGAIw3pmUrt3x52g4MlK0D+ZBpTDNnpq172ToAxEdzV7np00tXgNzINKYzzihdAYBuQXNXubVrS1eA3Mg0JqZkAbSFNXcAAACB0NwBAAAEQnNXuQUL0kAcZBrTtGlpAMB4Y81d5fbtK10BciPTmPbsKV0BgG5Bc1e5DRtKV4DcyDSm228vXQGAbkFzV7mpU0tXgNzINCamZAG0hTV3AAAAgdDcVW7dujQQB5nGtGhRGgAw3mjuKrd1axqIg0xjuvLKNABgvLHmrnKLF5euALmRaUxXXFG6AgDdguaucnPmlK4AuZFpTEzJAmgL07IAAACB0NxVbtu2NBAHmca0ZUsaADDemJat3OrVaTswULYO5EOmMc2dm7buZesAEB/NXeVmzSpdAXIj05jOOad0BQC6Bc1d5fr6SleA3Mi0HrNnS4ODh/ee3t6R9/f0SOvXj70mAKC5A4AjNDiYrxnv789zHgDghgoAAIBAaO4qN3fugwu1EQOZxkSuANpCcwcAABAIa+4qx+My4iHTmMgVQFu4cgcAABAIzR0AAEAgTMtWbtWqtOXZaHFMhEz37z/w89hKmTy5dAVjMxFyBdAdaO4qd/31pStAbhMh06VLS1cQz0TIFUB3oLmr3IoVpStAbmQaE7kCaAvNXeUm2tQZxo5MYyJXAG3hhgoAAIBAaO4qd/XVaSAOMo2JXAG0hWnZyl1+edrOmVO2DuRDpjGRK4C20NxV7uyzS1eA3Mg0JnIF0Baau8otWVK6AuRGpjGRK4C2sOYOAAAgEJq7yu3dmwbiINOYyBVAW5iWrdzChWk7MFC2DuRDpjGRK4C20NxVbsqU0hUgNzKNiVwBtIXmrnIbN5auALmRaUzkCqAtrLkDAAAIhOYOAAAgEJq7yi1blgbiINOYyBVAW1hzV7ndu0tXgNzINKZD5bp/v9TbO/bv09MjrV8/9vMAqBfNXeXWrCldAXIj05gOlevSpXm+T39/nvMAqBfNXeVmzChdAXIj05jIFUBbWHMHAAAQCM1d5TZtSgNxkGlM5AqgLTR3ldu8OQ3EQaYxkSuAtrDmrnLz5pWuALmRaUzkCqAtNHeVmz+/dAXIjUxjIlcAbWFaFgAAIBCu3FVu16605TELcZDp+Mv1wODJk0d/LLkCaAvNXeWWL0/bgYGydSAfMh1/uR4YfDjIFUBbaO4qN3166QqQG5nGRK4A2kJzV7m1a0tXgNzINCZyBdAWbqgAAAAIhOYOAAAgEJq7yi1YkAbiINOYyBVAW1hzV7l9+0pXgNzINCZyBdAWmrvKbdhQugLkRqYxkSuAttDcVW7q1NIVIDcyjamtXHM9oDm3nh5p/frSVQDdgeYOAAIp8YDm0ejvL10B0D24oaJy69algTjINCZyBdAWmrvKbd2aBuIg05jIFUBbmJat3OLFpStAbmQaE7kCaAvNXeXmzCldAXIj05jIFUBbmJYFAAAIhOauctu2pYE4yDQmcgXQFqZlK7d6ddoODJStA/mQaUzkCqAtNHeVmzWrdAXIjUxjIlcAbaG5q1xfX+kKkBuZxkSuANrCmjsAAIBAaO4AAAACobmr3Ny5aSAOMo2JXAG0heYOAAAgEG6oqByPVYiHTGMiVwBt4codAABAIDR3AAAAgdDcVW7VqjQQB5nGRK4A2sKau8pdf33pCpAbmcZErgDaQnNXuRUrSleA3Mg0JnIF0Baau8r19pauALmRaUzkCqAtrLkDAAAIhOaucldfnQbiINOYyBVAW5iWrdzll6ftnDll60A+ZBoTuQJoC81d5c4+u3QFyI1MYyJXAG2huavckiWlK0BuZBoTueYxe7Y0OJjnXD090vr1ec4FTCQ0dwCAagwOSn19ec7V35/nPMBEww0Vldu7Nw3EQaYxkSuAtnDlrnILF6btwEDZOpAPmcZErgDaQnNXuSlTSleA3Mg0JnIF0Baau8pt3Fi6AuRGpjGRK4C2sOYOAAAgEJo7AACAQJiWrdyyZWm7dm3ZOpAPmcbU7bnu3y/19o79PJMnj/0cQHQ0d5Xbvbt0BciNTGPq9lyXLi1dAdA9aO4qt2ZN6QqQG5nGRK4A2kJzV7kZM0pXgNzINCZyBdAWbqgAAAAIhOaucps2pYE4yDQmcgXQFpq7ym3enAbiINOYyBVAW1hzV7l580pXgNzINCZyBdAWmrvKzZ9fugLkRqYxkSuAtjAtCwAAEAjNXeV27UoDcZBpTOQKoC1My1Zu+fK0HRgoWwfyIdOYyBVAW2juKjd9eukKkBuZxkSuANpCc1e5bv0l5JGRaUzkCqAtrLkDAAAIhOYOAAAgEJq7yi1YkAbiINOYyBVAW1hzV7l9+0pXgNzINCZyBdAWmrvKbdhQugLkRqYxkSuAttDcVW7q1NIVIDcyjYlcAbSFNXcAAACB0NxVbt26NBAHmcZErgDaQnNXua1b00AcZBoTuQJoC2vuKrd4cekKkBuZxkSuANpCc1e5OXNKV4DcyDQmcgXQFqZlAQAAAqG5q9y2bWkgDjKNiVwBtIVp2cqtXp22AwNl60A+ZBoTuQJoC81d5WbNKl0BciPTmMgVQFto7irX11e6AuRGpjGRK4C2sOYOAAAgEJo7AACAQGjuKjd3bhqIg0xjIlcAbaG5AwAACIQbKirHYxXiIdOYyBVAW2juAAAYg9mzpcHBPOfq6ZHWr89zLnQvmjsAAMZgcDDfo276+/OcB92NNXeVW7UqDcRBpjGRK4C2cOWuctdfX7oC5EamMZErgLbQ3FVuxYrSFSA3Mo2JXAG0heaucr29pStAbmQaE7kCaAtr7gAAAAKhuavc1VengTjINCZyBdCWCT0tm+vZQZMnj/0cE9Xll6ftnDll60A+ZBoTuaJWPMevPhO6ucv57KCozj67dAXIjUxjIlfUiuf41WdCN3c4tCVLSleA3Mg0JnIF0BbW3AEAAARCc1e5vXvTQBxkGhO5AmgL07KVW7gwbQcGytaBfMg0JnIF0Baau8pNmVK6AuRGpjGRK4C20NxVbuPG0hUgNzKNiVwBtIU1dwAAAIHQ3AEAAARCc1e5ZcvSQBxkGhO5AmgLa+4qt3t36QqQG5nGRK4A2kJzV7k1a0pXgNzINCZyBdAWmrvKzZhRugLkRqYxkSuAtrDmDgAAIBCau8pt2pQG4iDTmMgVQFto7iq3eXMaiINMYyJXAG1hzV3l5s0rXQFyI9OYyBVAW2juKjd/fukKkBuZxkSuE8/+/VJv79jPM3ny2M8B5ERzBwDoSkuXlq4AGB+suavcrl1pIA4yjYlcAbSFK3eVW748bQcGytaBfMg0JnIF0Baau8pNn166AuRGpjGRK4C20NxVbu3a0hUgNzKNiVwBtIU1dwAAAIHQ3AEAAATCtGzlFixI240bS1aBnMg0JnLFaOR69p4k9fRI69fnOVcuuT5fzs82e7Y0OJjnXBPl75zmrnL79pWuALmRaUzkitHI+ey9/v5858ol1+fL+dkGB6W+vjznmih/5zR3lduwoXQFyI1MYyJXAG2huavc1KmlK0BuZBoTuQJoCzdUAAAABEJzV7l169JAHGQaE7kCaAvNXeW2bk0DcZBpTOQKoC2suavc4sWlK0BuZBoTuQJoC81d5ebMKV0BciPTmMgVQFto7gAACCjXA4MnTx77OdAumrvKbduWtrmeaI7yyDQmckXbcj4QGXWhuavc6tVpOzBQtg7kQ6YxkSuAttDcVW7WrNIVIDcyjYlcAbSF5q5yuX4fHiYOMo2JXAG0hefcAQAABEJzBwAAEAjNXeXmzk0DcZBpTOQKoC00dwAAAIGYu4/+YLPvS7p1/MqZEE6UdGfpItAKsu4eZN1dyLt7dHvWT3L3k4bvPKzmrhuY2XZ3n1m6Dow/su4eZN1dyLt7kPXImJYFAAAIhOYOAAAgEJq7h+ovXQBaQ9bdg6y7C3l3D7IeAWvuAAAAAuHKHQAAQCBd+7tlzewoSS+V1Cvpx5Judvd/KVsVgMNlZse5+72l68D4I+vuQdZj05XNnZn1SPpbSfdKerW7/+ew12dIWiXpbkn3SXqEpBXu/r22a8WRM7MXSVoo6TZJJ0h6gqRL3f3ajmPIulJm9mhJSyRdJGnqCK+PKlt+Bia+UWT9PyW9SdLpkvZJ+jdJbyXr+hwq62HHPkLStZIuc/eNw17r7qzdvauGpN9QCvufJT18hNefLGmvpGUd+94k6SZJjytdP2PUOb9E6R/5kzv2PV/SPZKeQtZ1D0kvkrRa6aq7j/D6qLLlZ2Dij1FkvUjSjZJWSnqnpG9Icklfk3QsWdczDpX1CMe/s8l6wbD9XZ918QJa/bDSyZK+J2mXpMkHOObvmh+Kh3Xse1jTKFxW+jMwRp31JknbR9h/l6QLyTrGkPQfB/g//FFly89APWOkrCVNkvSpYU3c0ZK+2Pyf/nlkXd840P+uhx3zPyR99ADNXddn3W03VPRJOknSGncfHP6imU2V9HJJn3b3+4f2N3/+rKSFZnZCO6VijI6TdLqZTRvaYWaPV5qe3UPWYdw3fMdos+VnoDoPyVrSc5SmX/cP7XD3n0q6ovlyusS/7RUaKeufa6Zu3yLp4hFeI2t10d2yZnaMpFc3Xz7OzHaa2T1m9mUzu8DMTOkfimOULusPd5PSfyX+WjsVY4w+pLTG4rNm9iwzO1bSeyR9RGlKnqzjGm22/AxUzt0/5e5fGuGlu5rt0Hpqso7lfUoXa0a64YKs1UXNnaQzlUK9W9LW5uvHS9qp1Ai8uflakn4wwvuH/rF40rhWiSzcfaukVyjldZ3SZf7/7e5/6O4/E1lHNtps+RmI6xmS7pD08eZrsg7CzM6T9G1333mAQ8ha3dXcndJs/9rdv+jJXZL+RGkd3sWSpjTH3D3C+4emcU8e3zKRi7t/VNL7JZnSf6n9mZkN/Y+arOMabbb8DATUzNK8StIr3f0nzW6yDqBZZvMHSrMwB0LW6q7mbmju/bbOne5+n9KVvEcp3Ump5s/DTWq2e8elOmRlyQpJvyTpiUoLb58j6XNmdoqkO5tDyTqe0WbLz0BMH5B0hbt/sWMfWVeueTbt+5XugL3/IIeStbqrubul2U4Z4bU9zXZXsz1hhGOG3nd7vpIwjt4qaamk17j7d5UW2L5X6XL8Sj2Y4wkjvJes6zbabPkZCMbMFkjaL2ndsJfIun6vk3SN0g1xxzbrqI9pXju62fdwkbWk7mrubpL0Q6Xn3wz3iGZ7o9ICzdNGOOZpSv9o/Pu4VIfcXinpumZ9ndz9AaV1lXskPVfS50XWUY02W34GAjGzlyittbvQm2dfdCDr+p2ntD7+3o7xzea1/ubrFSJrSV30Gyrc/V4z2yxpnpkd07EWQ5KeLunLkr4jaUDSi8zsqKYhGLoc/GJJH3X3H7RcOo7MnZKO79zh7g+Y2TZJx7n73WZG1gEdTrb8DMRgZi+U9FuSLups7JolGKe5+zVkXb3Fkh49bN/Jkv6v0oOP/1XSbfzbnnTTlTtJeoeknyg91VqSZGbPlPTrkl7f/KPwVqUreX/a8b4LlR6A2NdapRirD0g6y8yeM7TDzH5J6ardB5tdZF2/R0mSmR09bP9os+VnoB4jZm1mvyFpg9LdkW81sxXNeK/SeuqvNoeSdT0ekrW7f8Pdr+scShdlJGl3s29oTX3XZ901V+4kyd1vN7NnS3qXmX1S6Tk4UySd5e5fbY7Z3TQEq83sWUp3Wh4v6TnuvrtQ6ThM7v5RMztH0jvM7DtK//A/UdIfuvsnmmPIulJm9nylaZozm10fNrNPuPtmafTZ8jMw8R0sazObqfS4k0mSLhnh7X/v7t+XyLoGh/rf9WiRtWQPXZoAAACAWnXbtCwAAEBoNHcAAACB0NwBAAAEQnMHAAAQCM0dAABAIDR3AAAAgdDcAQAABEJzByALM/tLM3Mzu8fMdpnZjWb202bfzWb2LTPb13x9Q/OeKc2xlxUuX5JkZreY2d3NZ1lpZqe39H0vNrPVZva15u/nBW18XwAxddVvqAAw7jZIep273yulZknSkySd7e63NL/f8bXNkNKT449qthPFx939wja/obu/U5LM7FpJV7X5vQHEQ3MHIBeX9Jahxm7EA9Iv8f5Q8/tA5e57JT25pfpG64D1t+Cegt8bQBBMywLI5Up3/+4oj33HuFYCAF2M5g5AFu7+jcM49iZJMrNjzOxlZvbOodfM7BlmtsrMvm5mzzKzN5rZf5rZ98zs3c0xzzWzq83sLjPbZmYzOs9vZlPNbJ2ZbTCzr5jZP5nZ447kc5nZi81sT7MW7q/NbHqz/zgze3ezjvB5ZvY4M1tmZjvN7MlmdllT8wXN8ec3NS0xsw+b2cCR1AMAh8K0LIAimoZsiaTXS7pV0sXNS/dImiHpdEkrJH1E0nMkXSjprWZ2oqQbJC2U9HhJX5D055Je0Zz3BEnXSXqlu+80s+Ml/ZekL5rZM939sKY+3f0aM3tTU8ed7r672X+vmd0q6b3u/u9m9jJJr5L0a5LeKOlmSV+X9Cgz65G0VtJj3f3HTU3/cjh1AMBoceUOQBHuvqu5ceHOYfu/Lema5su/cvcBd79TqbmSpEe6+2Xuvsfdr5e0U9JZHae4WNJX3H1nc767JX1c0nRJ846w3H+Q9G1JrzGzSR37f1/S3zTf518lXd3s/5S7r3X3F7r7+yQ9VdIkSb0dNa0+wloA4KBo7gCUNtKVtJ8N20rSrmb7k2HH7pJ0UsfXvynpqWa2cWg0r39F0pQjKdDdf6Z0dXCKpNdIP7/y+EN3v6Pj0J822x3DTvFpSXskfdLMLjWzSe7+6SOpBQAOhWlZAFVw9wfMRnxiygP6xf9QnS6pfxweZ/IRSW+TtNzMPiTpf0naOJo3uvtdZjZT6QrgxZJeZWa/5+5fylwjAHDlDkA431dar/cQZvaYIz2pu/9E0nuUntv3CqUrhKO6KcLMjm2u8L1Y0gWSHiXpC2b29COtBwAOhOYOwHga+jemzYcUb5M028xe2rnTzJ4i6dIxnvvDSs3jWknbm4ZvNOaZ2dPc/WfufoWk5ynNnMweYz0A8BA0dwDGRfPbKCY3X04+yKHHSTqpOX7I1GZ7bMf5Ht388eRh739M8/oxzdeXKK19u8rMPmRmi83sCkkflbTmsD9Ih+YBzVdIeqzSb+MY7vjhdXd4Q8efv6+0dvD6sdQDACOhuQOQnZmdp9T8DDVk68zswo4GTGb2K2a2QqlReqSkvzSzJ5nZqyQtag57s5nNNrPTJP1Fs++s5jl4x5nZUknPb/a/08ye5O43Kj065fNKd7O+Tanp+m13H7opYyy2SLpR6S7dzs/8Oj14N+77zGz4nbl/bGYfM7OVkj4k6Xx3vzZDPQDwC8zdS9cAABNC87twP+vuCw5yzJsl3ePul43D93+BpM9IeqG7fzb3+QF0B+6WBYBfNG34DjMzd/fmAcoLJM0cp+99wjidF0AXobkDgAf9u6THmtl7JO2T9DGlKd1rzOzrSlPIa9z9Rzm/qZldJOkYpd9ucY2ku3KeH0B3YVoWAA7CzE6R9G+STpT0QUmr3P2BslUBwIHR3AEAAATC3bIAAACB0NwBAAAEQnMHAAAQCM0dAABAIDR3AAAAgfx/ZEvxMNKnAmAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 792x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "frac_error=radius_accuracy['10m']\n",
    "for keys,value in network.items():\n",
    "    detector=keys\n",
    "    filepath=value\n",
    "    listfiles=glob.glob('{}/nevents_*.txt'.format(filepath))\n",
    "    if(len(listfiles)==0):\n",
    "        break\n",
    "    number_of_years=[]\n",
    "    for file in listfiles:\n",
    "        events,time,err=np.loadtxt(file,unpack=True)\n",
    "        for i in range(len(events)):\n",
    "            if (err[i]<frac_error):\n",
    "                number_of_years.append(time[i])\n",
    "                break\n",
    "    mean=np.mean(number_of_years)\n",
    "    stdev=np.std(number_of_years)\n",
    "    plt.hist(number_of_years,bins=30,density=True,alpha=0.3,color='blue')\n",
    "    plt.hist(number_of_years,histtype='step',density=True,bins=30,alpha=0.7,color='blue')\n",
    "    plt.axvline(mean,linestyle=':',linewidth=2,color='blue')\n",
    "    plt.axvline(mean+stdev,linestyle=':',linewidth=2,color='blue')\n",
    "    plt.axvline(mean-stdev,linestyle=':',linewidth=2,color='blue')\n",
    "    plt.xlabel(\"Time [yrs]\")\n",
    "    plt.yticks([])\n",
    "    plt.title(str(keys)+\":\"+str(mean)[0:4]+\" years\")\n",
    "    plt.show()\n",
    "    \n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d2d28b9f",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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
}