{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "\n", "
\n", "**This is a fixed-text formatted version of a Jupyter notebook.**\n", "\n", " You can contribute with your own notebooks in this\n", " [GitHub repository](https://github.com/gammapy/gammapy-extra/tree/master/notebooks).\n", "\n", "**Source files:**\n", "[spectrum_simulation_cta.ipynb](../_static/notebooks/spectrum_simulation_cta.ipynb) |\n", "[spectrum_simulation_cta.py](../_static/notebooks/spectrum_simulation_cta.py)\n", "
\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Spectrum simulation for CTA\n", "\n", "A quick example how to simulate and fit a spectrum for the [Cherenkov Telescope Array (CTA)](https://www.cta-observatory.org)\n", "\n", "We will use the following classes:\n", "\n", "* [gammapy.spectrum.SpectrumObservation](http://docs.gammapy.org/dev/api/gammapy.spectrum.SpectrumObservation.html)\n", "* [gammapy.spectrum.SpectrumSimulation](http://docs.gammapy.org/dev/api/gammapy.spectrum.SpectrumSimulation.html)\n", "* [gammapy.spectrum.SpectrumFit](http://docs.gammapy.org/dev/api/gammapy.spectrum.SpectrumFit.html)\n", "* [gammapy.scripts.CTAIrf](http://docs.gammapy.org/dev/api/gammapy.scripts.CTAIrf.html)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Setup" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import astropy.units as u\n", "from gammapy.irf import EnergyDispersion, EffectiveAreaTable\n", "from gammapy.spectrum import SpectrumSimulation, SpectrumFit\n", "from gammapy.spectrum.models import PowerLaw\n", "from gammapy.scripts import CTAIrf" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Simulation" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Define obs parameters\n", "livetime = 10 * u.min\n", "offset = 0.3 * u.deg\n", "lo_threshold = 0.1 * u.TeV\n", "hi_threshold = 60 * u.TeV" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# Define spectral model\n", "index = 2.3 * u.Unit('')\n", "amplitude = 2.5 * 1e-12 * u.Unit('cm-2 s-1 TeV-1')\n", "reference = 1 * u.TeV\n", "model = PowerLaw(index=index, amplitude=amplitude, reference=reference)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Load IRFs\n", "filename = '$GAMMAPY_EXTRA/datasets/cta/perf_prod2/South_5h/irf_file.fits.gz'\n", "cta_irf = CTAIrf.read(filename)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "NDDataArray summary info\n", "energy : size = 500, min = 0.005 TeV, max = 495.450 TeV\n", "offset : size = 45, min = 0.050 deg, max = 4.450 deg\n", "Data : size = 22500, min = 0.000 m2, max = 4033200.000 m2\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEOCAYAAAAg4qj5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VdW99/HPj5Awz5NIQFCjglJBU8SxioqoVWyrrUMr\ntVZaH+1t721vq70+V2/tYG/vU3t9Wq1YLepjxaFV0YsiWucBCQUZVSKChHkMYQok+T1/7BU8hOSc\nk+Fk5yTf9+t1Xjl77bX3+m0h/thrr72WuTsiIiJxaBd3ACIi0nYpCYmISGyUhEREJDZKQiIiEhsl\nIRERiY2SkIiIxEZJSEREYqMkJCIisVESEhGR2CgJiYhIbNrHHUBL17dvXx86dGjcYYiIZJW5c+du\ncvd+qeopCaUwdOhQioqK4g5DRCSrmNnKdOqpO05ERGKjJCQiIrFREhIRkdgoCYmISGyUhEREJDYZ\nT0JmlmNm88zsubA9zMxmm9kyM3vMzPJCeYewXRz2D004x82h/EMzOy+hfEIoKzazmxLK692GiIg0\nv+YYov19YCnQPWz/GrjT3aeZ2R+Ba4F7ws+t7n6kmV0e6n3NzEYAlwPHAocCL5nZUeFcfwDOBUqA\nOWY23d2X1LeNTP8HEJHWrarKKSuvoHTXPkp372Pb7r2U7o6+7yyvYF+l85uZH9Z67I/GH0VuTjs6\n5+XQvVMuPTrl0rNzHj065dKvWwe6dmjdb9Jk9OrMLB+4EPgF8C9mZsA44MpQ5UHgNqIEMTF8B3gS\n+H2oPxGY5u7lwCdmVgyMCfWK3X15aGsaMNHMlta3DXf3Jr1wEclqO8orWLttN2tK97BlZzlbd+7j\nZ88tyUhb//XiR/U+5p6rTmBYvy4M7dOFjrk5GYiq+WQ6xf4O+DHQLWz3Aba5e0XYLgEGhe+DgFUA\n7l5hZqWh/iDg3YRzJh6zqkb5SQ1sY1Ni0GY2GZgMMGTIkHpftIi0XPsqq1hXuoe1pXtYs203q7ft\nrvMupSn909kF5OUYuTnRU5CKKmdfZRW/e2lZvc91/SP/qLX8vqsLGTmoBwO6dyD6N3zLl7EkZGZf\nBDa4+1wzO7O6uJaqnmJfXeW1Pc9KVj9V+58VuE8BpgAUFhbqLkmkiZVs3cXStWWU7t5Hh/btOKxP\nZwb17MSabXtYvW0Xu/dVUlkFuTlGXk47cnPakdu+3YHbOe3Ia2+4w559VezeV8me8Nm9r5KNZeVs\n2rGXP772cZPE/P2zC+jZOZdenfPo2fmzLrMenXLp3rE97XMa9oj9B+ccVec+d2fbrn2sLd3Duu27\nWb1tDys37eRPb35Sa/3rHjpwdpfbLhrB2CP6cFT/brRr1zKTUibvhE4FLjazC4CORM+Efgf0NLP2\n4U4lH1gT6pcAg4ESM2sP9AC2JJRXSzymtvJNDWhDRDJkY1k5Mxau5dbpi+MOpU63XjSCgT06Mahn\nJw7t2ZHeXfJaxJ2EmdGrSx69uuQx4tDu+8tv+eIIILqrW7l5F0vWbmfx6lLufX35Acff9uyBXYi/\n/spIzjq6P/27d8x88Gmy5ngcEu6EfuTuXzSzJ4C/JgwaWODud5vZDcBId/9uGDTwZXf/qpkdC/yF\n6DnQocDLQAHRXc1HwNnAamAOcKW7L65vG8liLywsdM0dJ1J/u/dWMvzfX4g7DP7zK5+jb7c8+nXt\nSL9uHejTNW9/l1hr5O6s2LyL2cs3c9PfFtZZ7wfnFHDO8AEce2j3jCRcM5vr7oUp68WQhA4HpgG9\ngXnA19293Mw6Ag8Do4nuTi5PGHTwb8C3gArgB+7+fCi/gOjuKgd4wN1/Ecrr3UZdlIRE6m9jWTmf\n/8VLB5Tdfslx5PfqRI9OuRx7aHc6tM/uB+rZwt1ZtWU3ry/byC1PL6q1zvfPLuDiUYdyRL+uTdZu\ni0pC2UxJSKR+yvbs40t3v03xhh0c3rcLd35tFMcP7hl3WBLs3lvJO8s38a2ptf9/7ZYLh/Ol0YPo\n07VDo9pREmoiSkIi9fOTJxfwWNEqCvp35dHJY+nbyP+ZSeZUVFbx9sebufqB92rdP23yWMYe3qdB\n5043CbXejlERaXYvLVnPY0XRmxN3X3WCElAL1z6nHWcc1Y8Vd1zIB7dP4J6rTjhg/1vFm+o4sglj\nyHgLItIm7Civ4KdPRQ/C/+2C4RQM6JbiCGlJOubmcP7Igay440LWlu7mr3NLmDhqUOoDG0lJSESa\nxD2vFrOhrJzjB/fkW6cNizscaYSBPTpx47iCZmlL3XEi0mglW3dx3xvRC5S3XjSCnBb6YqS0PEpC\nItJov5n5IXsrqrj4+EM5YUivuMORLKIkJCKNsmx9Gc/MjyYl+fGEo2OORrKNkpCINMrvXykG4Otj\nh5Dfq3PM0Ui2URISkQb7eOMOnn1/Dbk5xvVnHhl3OJKFlIREpMGmvLacKodLT8xnUM9OcYcjWUhJ\nSEQapGzPvv0vpn779MNjjkaylZKQiDTI9PejwQhjhvVu0okvpW1REhKRBpn2XnQXdMWYwSlqitRN\nSUhE6m3R6lIWri6le8f2nH/cwLjDkSymJCQi9fbYnOgu6Msn5NMxV+sCScNlLAmZWUcze8/M3jez\nxWb2H6F8qpl9Ymbzw2dUKDczu8vMis1sgZmdkHCuSWa2LHwmJZSfaGYLwzF3WVge0Mx6m9msUH+W\nmfVK1YaIpGfPvkqenr8agK99Xl1x0jiZvBMqB8a5+/HAKGCCmY0N+/7V3UeFz/xQdj7Rst0FwGTg\nHogSCnArcBLREt+3VieVUGdywnETQvlNwMvuXkC0HPhNydoQkfS9vHQDZXsqGDmoB8MHdo87HMly\nGUtCHtkRNnPDJ9kKehOBh8Jx7wI9zWwgcB4wy923uPtWYBZRQhsIdHf3dzxame8h4JKEcz0Yvj9Y\no7y2NkQkTU/Ni+6CLhmd+Wn+pfXL6DMhM8sxs/nABqJEMjvs+kXoDrvTzKpXvRoErEo4vCSUJSsv\nqaUcYIC7rwUIP/unaENE0rBt115eWroegIuO17/fpPEymoTcvdLdRwH5wBgzOw64GTgG+DzQG/hJ\nqF7b3O/egPJk0jrGzCabWZGZFW3cuDHFKUXajv9ZuBaA0wv60r9bx5ijkdagWUbHufs24FVggruv\nDd1h5cCfiZ7zQHRXkviUMx9Yk6I8v5ZygPXV3Wzh54YUbdSMd4q7F7p7Yb9+/ep5tSKt1zPzol+X\n5lhxU9qGTI6O62dmPcP3TsA5wAcJycGIntUsCodMB64OI9jGAqWhK20mMN7MeoUBCeOBmWFfmZmN\nDee6Gngm4VzVo+gm1SivrQ0RSaFk6y7eW7GFjrntOO/YAXGHI61EJpf3Hgg8aGY5RMnucXd/zsz+\nbmb9iLrG5gPfDfVnABcAxcAu4BoAd99iZrcDc0K9n7n7lvD9emAq0Al4PnwA7gAeN7NrgU+By5K1\nISKpVU/Tc87wAXTrmBtzNNJaZCwJufsCYHQt5ePqqO/ADXXsewB4oJbyIuC4Wso3A2fXpw0RSa66\nK+4SdcVJE9KMCSKS0tK12/lwfRk9O+dyxlF6TipNR0lIRFKqniHhwpEDyWuv/21I09HfJhFJam9F\nFX+dGyWhL5+grjhpWkpCIpLUrCXr2bSjnKMHdOOEIb1SHyBSD0pCIpLUI7NXAnDV2CGEOYJFmoyS\nkIjU6eONO3j74810ys3RXHGSEUpCIlKnv8z+FICJow6lu94NkgxQEhKRWu3ZV8mTc6M5gq866bCY\no5HWSklIRGr13IK1lO7ex/H5PRiZ3yPucKSVUhISkVrtH5CguyDJICUhETnI4jWlzPt0GwBf1LpB\nkkFKQiJykOoBCd88ZSid8zI5z7G0dUpCInKAHeUVPB2W8L7ypCExRyOtnZKQiBzg6Xmr2bm3kjHD\nenPUgG5xhyOtnJKQiOzn7jwSuuKu0l2QNAMlIRHZr2jlVpau3U7vLnlMOO6QuMORNiCTy3t3NLP3\nzOx9M1tsZv8RyoeZ2WwzW2Zmj5lZXijvELaLw/6hCee6OZR/aGbnJZRPCGXFZnZTQnm92xARuO/1\n5QBcOWYIHdrnxByNtAWZvBMqB8a5+/HAKGCCmY0Ffg3c6e4FwFbg2lD/WmCrux8J3BnqYWYjgMuB\nY4EJwN1mlhOWDf8DcD4wArgi1KW+bYgIrNi0k1lL1wNw9Sl6N0iaR8aSkEd2hM3c8HFgHPBkKH8Q\nuCR8nxi2CfvPtmjK3onANHcvd/dPgGJgTPgUu/tyd98LTAMmhmPq24ZIm/fAW5/gDpedmE//bh3j\nDkfaiIw+Ewp3LPOBDcAs4GNgm7tXhColQPXUvIOAVQBhfynQJ7G8xjF1lfdpQBsibdq2XXt5oiia\nJ+7bpx8eczTSlmQ0Cbl7pbuPAvKJ7lyG11Yt/KztjsSbsDxZGwcws8lmVmRmRRs3bqzlEJHW5ZHZ\nn7J7XyVnHNWPow/RsGxpPklfhTazf0njHDvd/d5kFdx9m5m9CowFeppZ+3Ankg+sCdVKgMFAiZm1\nB3oAWxLKqyUeU1v5pga0UTPeKcAUgMLCwoOSlEhrUl5RydS3VwBw3enD4g1G2pxUd0L/CnQFuiX5\n/LC2A82sn5n1DN87AecAS4FXgEtDtUnAM+H79LBN2P93d/dQfnkY2TYMKADeA+YABWEkXB7R4IXp\n4Zj6tiHSZk2fv4aNZeUcc0g3Tjuyb9zhSBuTalKoh939Z8kqmFmXOnYNBB4Mo9jaAY+7+3NmtgSY\nZmY/B+YB94f69wMPm1kx0d3J5QDuvtjMHgeWABXADe5eGdq+EZgJ5AAPuPvicK6f1KcNkbbK3bn/\nzU8AuPa0YVq+W5qd6UYgucLCQi8qKoo7DJGMeGPZRr5x/3v069aBN39ylt4NkiZjZnPdvTBVvZQD\nE8zsGDM728y61iif0JgARSR+970R3QV985ShSkASi6RJyMz+ieh5yveARWY2MWH3LzMZmIhk1ofr\nynj9o410ys3RPHESm1TPhK4DTnT3HWGKmyfNbKi7/ze1D3cWkSzxpzeiKXouK8ynZ+e8mKORtipV\nEsqpnvXA3VeY2ZlEiegwlIREstaGsj08MTd6OfVbp2pYtsQn1TOhdWY2qnojJKQvAn2BkZkMTEQy\n56G3VwJw3rEDGNq3rgGuIpmXKgldDaxLLHD3Cne/GjgjY1GJSMbsq6zi0feiNYM0RY/ELWl3nLuX\nVH83s15Esw1UH7M7g3GJSIa8uWwTm3fu5cj+XSk8rFfc4Ugbl+qZEABmdjvwTaIJSKtfLKqeEVtE\nssjT81cDcMmoQ/VyqsQurSQEfBU4IiyZICJZamd5BS8ujtYMmjhqUIraIpmX7izai4CemQxERDLv\nxSXr2L2vksLDejG4d+e4wxFJ+07oV8A8M1tEtGIqAO5+cUaiEpGMeHpeNKH8xNG6C5KWId0k9CDR\nUtgLgarMhSMimbKxrJzXPorWx7pw5MCYoxGJpJuENrn7XRmNREQy6rkF0V3QOcP707uLZkiQliHd\nJDTXzH5FtB5PYnfcPzISlYg0uafnRaPiNCBBWpJ0k9Do8HNsQpmGaItkiU827eT9klK65OVwzvAB\ncYcjsl9ao+Pc/axaPkkTkJkNNrNXzGypmS02s++H8tvMbLWZzQ+fCxKOudnMis3sQzM7L6F8Qigr\nNrObEsqHmdlsM1tmZo+FFVYJq7A+FurPDpOvJm1DpDWrvguacNxAOuVpyQZpOVIt5fDFVCdIUqcC\n+KG7Dye6g7rBzEaEfXe6+6jwmRHOM4JopdNjgQnA3WaWE1Zm/QNwPjACuCLhPL8O5yoAtgLXhvJr\nga3ufiRwZ6hXZxuprlEkm7n7Zy+ojj405mhEDpSqO+43Zraa5DNm/xJ4rmahu68F1obvZWa2FEjW\nGT0RmObu5cAnYQnuMWFfsbsvBzCzacDEcL5xwJWhzoPAbcA94Vy3hfIngd9b9Gp4XW28kyQukaw2\nf9U2Vm7eRb9uHTjliL5xhyNygFRJaD3w2xR1lqVqJHSHjQZmA6cCN5rZ1UAR0d3SVqIE9W7CYSV8\nlrRW1Sg/CegDbHP3ilrqD6o+xt0rzKw01E/WhkirVN0Vd/Hxh5LTTtP0SMuSagLTMxvbQFgW/K/A\nD9x9u5ndA9xONLDhduD/AN+i9rstp/YuQ09SnyT7kh2TGPNkYDLAkCFacVKyV2WVM2NRNBH+xFHq\nipOWJ91pexrEzHKJEtAj7v43AHdf7+6V7l4F3MdnXW4lRLN0V8sH1iQp3wT0NLP2NcoPOFfY3wPY\nkuRcB3D3Ke5e6O6F/fr1a8ili7QIc1duZWNZOfm9OjFyUI+4wxE5SMaSUHgGcz+w1N1/m1Ce+Kr2\nl4jmpYPoHaTLw8i2YUAB8B4wBygII+HyiAYWTHd3B14BLg3HTwKeSTjXpPD9UuDvoX5dbYi0Ss8v\nWgvABSMHasZsaZHSfU+oIU4FvgEsNLP5oeynRKPbRhF1g60AvgPg7ovN7HFgCdHIuhvcvRLAzG4E\nZgI5wAPuvjic7yfANDP7OTCPKOkRfj4cBh5sIUpcSdsQaW2qqpwXQlfc+ccdEnM0IrWz6AYhjYpm\nxxENke5YXebuD2UorhajsLDQi4qK4g5DpN7+8elWvnz32xzaoyNv3TROd0LSrMxsrrsXpqqX7qJ2\ntwJnEiWhGUTv7LwJtPokJJKtnl8YdcVNOE5dcdJypftM6FLgbGCdu18DHA90yFhUItIo7s6MhaEr\nbqS64qTlSjcJ7Q6j2SrMrDuwATg8c2GJSGMsWr2d1dt2079bB04c0ivucETqlO7AhCIz60k0pHou\nsAONKhNpsWYsqu6KO4R2ekFVWrC0kpC7/6/w9Y9m9gLQ3d0XZC4sEWkod9//POj847R4nbRsaXXH\nWeTrZvbv7r4C2GZmY1IdJyLNb+naMlZs3kWfLnmMGdY77nBEkkr3mdDdwMnAFWG7jGhmaxFpYV4I\nXXHnHXeI5oqTFi/dZ0InufsJZjYPwN23Vq/dIyItywy9oCpZJN07oX1h3R0HMLN+QFXGohKRBlm2\nvoziDTvo2TmXsYf3iTsckZTSTUJ3AU8B/c3sF0Qvqv4yY1GJSINUvxs0fsQAcnMyOj+xSJNId3Tc\nI2Y2l+iFVQMucfelGY1MROqtesJSjYqTbJEyCZlZO2CBux8HfJD5kESkIZZv3MEH68ro1rE9px6p\nFVQlO6S8Xw8zJbxvZlrdTaQFez4MSDh3xADy2qsrTrJDuqPjBgKLzew9YGd1obtfnJGoRKTe9q8d\npK44ySLpJqH/yGgUItIon27exaLV2+mSl8NpBeqKk+yR7sCE1xK3zexU4ErgtdqPEJHmVH0XdPbw\nAXTMzYk5GpH0pd1xbGajzOw/zWwF8HMg6eg4MxtsZq+Y2VIzW2xm3w/lvc1slpktCz97hXIzs7vM\nrNjMFpjZCQnnmhTqLzOzSQnlJ5rZwnDMXWFJ8Qa1IZLNql9QvUDLNkiWSZqEzOwoM/t3M1sK/B5Y\nRbQa61nu/vsU564Afujuw4GxwA1mNgK4CXjZ3QuAl8M2RAvlFYTPZOCeEENv4FbgJGAMcGt1Ugl1\nJiccNyGU16sNkWy2ettu3l+1jU65OXzhqP5xhyNSL6nuhD4gejfoInc/zd3/L1CZzondfa27/yN8\nLyO6cxoETAQeDNUeBC4J3ycCD3nkXaCnmQ0EzgNmufsWd98KzAImhH3d3f0dj9Yof6jGuerThkjW\neiHcBY07pj+d8tQVJ9klVRL6CrAOeMXM7jOz6pdV68XMhgKjgdnAAHdfC1GiAqr/6TaI6E6rWkko\nS1ZeUks5DWhDJGvtX7ZBXXGShZImIXd/yt2/BhwDvAr8MzDAzO4xs/HpNGBmXYG/Aj9w9+3JqtYW\nQgPKk4aTzjFmNtnMisysaOPGjSlOKRKfdaV7KFq5lQ7t23HW0eqKk+yT1sAEd9/p7o+4+xeBfGA+\nnz1nqZOZ5RIloEfc/W+heH11F1j4uSGUlwCDEw7PB9akKM+vpbwhbdS83inuXujuhf369Ut1mSKx\nmbk46oo78+h+dOmQ7hsXIi1HvV+rDs9m7nX3ccnqhZFq9wNL3f23CbumA9Uj3CYBzySUXx1GsI0F\nSkNX2kxgvJn1CgMSxgMzw74yMxsb2rq6xrnq04ZIVpqhFVQly2Xyn06nAt8AFprZ/FD2U+AO4HEz\nuxb4FLgs7JsBXAAUA7uAayBKemZ2OzAn1PuZu28J368HpgKdgOfDh/q2IZKNNpaV896KLeTltGPc\ncHXFSXbKWBJy9zepexDD2bXUd+CGOs71APBALeVFwHG1lG+ubxsi2Wbm4nW4w+kFfeneMTfucEQa\npD4vqx5mZueE753MrFvmwhKRVPYv2zBSXXGSvdJKQmZ2HfAkcG8oygeezlRQIpLc5h3lvLt8C+3b\nGecOHxB3OCINlu6d0A1Ez3i2A7j7Mj5790ZEmtnzi9ZRWeWcVtCXHp3VFSfZK90kVO7ue6s3zKw9\nqd/JEZEMmf5+9GbBxccfGnMkIo2TbhJ6zcx+CnQys3OBJ4BnMxeWiNRlbelu3vskGiB67gh1xUl2\nSzcJ3QRsBBYC3yEa6nxLpoISkbr9z4Lqd4MOoZtGxUmWS3eIdvXEn/dlMhgRSe3Z0BV3kbripBVI\n907oYuAjM3vYzC4Mz4REpJmt2LST90tK6ZKXw7hjNDZIsl+6c8ddAxxJ9CzoSuBjM/tTJgMTkYM9\ntyC6Cxp/7CFaQVVahbTvaNx9n5k9TzQqrhNRF923MxWYiBxs+v6uOL2gKq1Dui+rTjCzqURzrl0K\n/AnQb4FIM/pwXRkfrd9Bz865nHakZneX1iHdO6FvAtOA77h7eebCEZG6TH9/NRCNistrX+8J8EVa\npLSSkLtfnulARKRuVVXOM/NDV9znNCpOWo+kScjM3nT308ysjANnSDCiSam7ZzQ6EQFgzootlGzd\nzaE9OjL28D5xhyPSZJImIXc/LfzUjNkiMXpqXtQVN3H0INq1q2uFFJHsk+7AhIfTKRORprdnX+X+\nWRK+PHpQzNGINK10n24em7gRXlY9MdkBZvaAmW0ws0UJZbeZ2Wozmx8+FyTsu9nMis3sQzM7L6F8\nQigrNrObEsqHmdlsM1tmZo+ZWV4o7xC2i8P+oanaEGnJZi1ZT1l5BZ/L70HBAHVKSOuSNAmF/2mX\nAZ8zs+3hUwasB55Jce6pwIRayu9091HhMyO0MwK4nCjZTQDuNrMcM8sB/gCcD4wArgh1AX4dzlUA\nbAWuDeXXAlvd/UjgzlCvzjZSXINI7P72jxJAd0HSOiVNQu7+q/A86Dfu3j18url7H3e/OcWxrwNb\n0oxjIjDN3cvd/ROi95HGhE+xuy8PS0lMAyaamQHjiBbaA3gQuCThXA+G708CZ4f6dbUh0mJtKNvD\n68s20b6daa44aZXS7Y57z8x6VG+YWU8zuyTZAUncaGYLQnddr1A2CFiVUKcklNVV3gfY5u4VNcoP\nOFfYXxrq13Wug5jZZDMrMrOijRs3NuwqRZrA9PlrqKxyzjy6P326dog7HJEml24SutXdS6s33H0b\ncGsD2rsHOAIYBawF/k8or224jzegvCHnOrjQfYq7F7p7Yb9+ejNd4uHuPDk36or7ygnqipPWKd0k\nVFu9es+k7e7r3b3S3auA+/isO6wEGJxQNR9Yk6R8E9AzYTbv6vIDzhX29yDqFqzrXCIt0vslpXyw\nrgyAccM1Y7a0TukmoSIz+62ZHWFmh5vZncDc+jZmZonzzX0JqB45Nx24PIxsGwYUAO8Bc4CCMBIu\nj2hgwXR3d+AVonnsACbx2UCJ6WGbsP/voX5dbYi0SI/O/hSAyWccTof2GkMjrVO6dzPfA/438FjY\nfpEUK6ua2aPAmUBfMysh6r4708xGEXWDrSBapRV3X2xmjwNLgArgBnevDOe5EZgJ5AAPuPvi0MRP\ngGlm9nNgHnB/KL8feNjMionugC5P1YZIS1O2Zx/PhmUbvvb5wSlqi2Qvi24S0qxs1tXdd2Qwnhan\nsLDQi4qK4g5D2phHZq/k355axEnDevPYd06OOxyRejOzue5emKpeujMmnGJmS4juIjCz483s7kbG\nKCJ1ePS9qCvuypOGxByJSGal+0zoTuA8YDOAu78PnJGpoETasoUlpSxavZ2enXM579hD4g5HJKPS\nXpTE3VfVKNLzFJEMeHROdBf05dH5WsJbWr10ByasMrNTAA+j1P4JWJq5sETaptLd+3g6zJh9xRgN\nSJDWL907oe8CNxDNMFBC9LLpDZkKSqSteqJoFbv2VnLqkX00Wam0CakWtfu1u/8EOMvdr2qmmETa\npMoq56F3VgLwzVOGxRyNSPNIdSd0gZnlAkknKxWRxvv7Bxv4dMsuBvfuxLhjNEOCtA2pngm9QDRF\nThcz205Y1hst7y3S5Ka+/QkAk04eSo5WT5U2ItWd0C3u3gP4n4RlHPb/bI4ARdqCD9eV8VbxZgAu\nK9SABGk7UiWhd8LP7ZkORKQtu/f1jwG4+uTD6NEpN+ZoRJpPqu64PDObBJxiZl+uudPd/5aZsETa\njpKtu5g+fw057YzrTj887nBEmlWqJPRd4CqgJ3BRjX0OKAmJNNKf3viEiirnS6MHMbh357jDEWlW\nSZOQu78JvGlmRe5+f7K6IlJ/W3buZVqYIeE7X9BdkLQ9SZ8JmdmPAdz9fjO7rMa+X2YyMJG2YOpb\nn7BnXxVnH9OfYw7RWB9pe1INTLg84XvNd4UmNHEsIm3Ktl17uevvxQBcf+YRMUcjEo9UScjq+F7b\n9oE7zR4wsw1mtiihrLeZzTKzZeFnr1BuZnaXmRWb2QIzOyHhmEmh/rIwSKK6/EQzWxiOucvMrKFt\niMThvjeWA3B6QV8Kh/aOORqReKRKQl7H99q2a5rKwXdLNwEvu3sB8HLYBjifaLntAmAycA9ECYVo\nRdaTgDHArdVJJdSZnHDchIa0IRKHzTvK+fNbKwD4l3OPijcYkRilSkLHm9l2MysDPhe+V2+PTHag\nu79OtLx2oonAg+H7g8AlCeUPeeRdoKeZDSRaw2iWu29x963ALGBC2Nfd3d/xaGnYh2qcqz5tiDS7\ne19fzq4BSCIbAAAPnUlEQVS9lYw7pj+jh/RKfYBIK5VqdFxTL2YywN3XhnOvNbPqCbIGAYnrFZWE\nsmTlJbWUN6SNtY29KJH6WL99Dw+9swKAfz5Hd0HStqW9qF2G1fZ8yRtQ3pA2Dq5oNtnMisysaOPG\njSlOK1I//zXzQ/bsq2LCsYcwMr9H3OGIxKq5k9D66i6w8HNDKC8BEifMygfWpCjPr6W8IW0cxN2n\nuHuhuxf269evXhcoksyi1aU8MTe6ib/p/GNijkYkfs2dhKYD1SPcJgHPJJRfHUawjQVKQ5faTGC8\nmfUKAxLGAzPDvjIzGxtGxV1d41z1aUOkWbg7P/+fJQB8+7RhDO3bJeaIROKX7vLe9WZmjwJnAn3N\nrIRolNsdwONmdi3wKVD9AuwM4AKgGNgFXAPg7lvM7HZgTqj3M3evHuxwPdEIvE7A8+FDfdsQaS4v\nLlnPu8u30LNzLt8bVxB3OCItgkWDy6QuhYWFXlRUFHcYkuV2llcw/s7XWb1tN7ddNIJvnqqVU6V1\nM7O57l6Yql5LGZgg0qr97qWPWL1tN8cN6s7Xxx4WdzgiLYaSkEiGLVpdyn1vRKum/upLn6N9jn7t\nRKrpt0Ekgyoqq/jpUwsB+NapwzQkW6QGJSGRDPr9K8UsKCnl0B4d+ZfxejFVpCYlIZEMmffpVn73\n0jIA/uurx9O1Q8YGo4pkLSUhkQzYWV7BPz82H4DrTh/GKUf0jTkikZZJSUikibk7tzy9iBWbd3HM\nId340XlHxx2SSIulJCTSxB5+dyVPzVsNwF1XjKZD+6aeB1ik9VASEmlCc1du4d+fWQxECeioAd1i\njkikZVMSEmkiq7bs4iv3vANEw7EvPv7QmCMSafmUhESawNade5n05/cAOPnwPtx8gWbIFkmHxoyK\nNFLprn2Mvn0WAMcc0o17rz6RXM2KIJIW/aaINMK2XXu56v53929PvWYM3TvmxhiRSHbRnZBIA23Y\nvodJf57D0rXbOaxPZx69biyH9OgYd1giWUVJSKQBPlxXxnm/e33/9mOTT1YCEmkAJSGRepq1ZD3X\nPfTZGlNzbzmHPl07xBiRSPaK5ZmQma0ws4VmNt/MikJZbzObZWbLws9eodzM7C4zKzazBWZ2QsJ5\nJoX6y8xsUkL5ieH8xeFYS9aGSDq279nHTX9dsD8BXThyIB/cPkEJSKQR4hyYcJa7j0pYee8m4GV3\nLwBeDtsA5wMF4TMZuAeihEK0ZPhJwBjg1oSkck+oW33chBRtiNTJ3Xlm/mo+d9uLTJuzCoBbLhzO\n/71iNB1zNRuCSGO0pO64icCZ4fuDwKvAT0L5Qx6tQ/6umfU0s4Gh7ix33wJgZrOACWb2KtDd3d8J\n5Q8BlwDPJ2lDpFbvr9rGxD+8dUDZzB+cwdGHaCYEkaYQVxJy4EUzc+Bed58CDHD3tQDuvtbM+oe6\ng4BVCceWhLJk5SW1lJOkDZH93J23ijfz9ftnH1D+66+M5LITB9OuncUUmUjrE1cSOtXd14QkMMvM\nPkhSt7bfeG9AedrMbDJRdx5Dhgypz6GSxaqqnBeXrOPuVz9mQUnp/vLvnHE43/nCEfTukhdjdCKt\nUyxJyN3XhJ8bzOwpomc6681sYLhDGQhsCNVLgMEJh+cDa0L5mTXKXw3l+bXUJ0kbNeObAkwBKCws\nrFcCk+yzt6KKp+ev5sdPLjig/Efjj+IbY4fSo7NePhXJlGYfmGBmXcysW/V3YDywCJgOVI9wmwQ8\nE75PB64Oo+TGAqWhS20mMN7MeoUBCeOBmWFfmZmNDaPirq5xrtrakDZo194KHnjzE4665fmDEtDS\nn03gxnEFSkAiGRbHndAA4Kkwaro98Bd3f8HM5gCPm9m1wKfAZaH+DOACoBjYBVwD4O5bzOx2YE6o\n97PqQQrA9cBUoBPRgITnQ/kddbQhbci2XXt58O2V3PnSRweU//arx3PR8Ydq3jeRZmTRoDOpS2Fh\noRcVFaWuKC3ehu17GPPLlw8qn/KNEzln+AANOBBpQmY2N+EVnDq1pCHaIhlRtmcf9762nN+/UnxA\n+aPXjWXs4b0Jd+UiEgMlIWm19lZU8ZfZK7nt2SX7y847dgDfG1fAcYN6xBiZiFRTEpJWx92ZsXAd\nN/zlHweU//X6UzjxMM3UJNKSKAlJqzJ7+Wa+NuXdA8qmfONEzh0xQN1uIi2QkpC0CsvWl3Huna8f\nUPbLL43kq4X5tNdoN5EWS0lIstr67Xu4c9ZH+ycWBfjnc47i26cPo0sH/fUWaen0WypZacvOvfzx\ntY+Z8vry/WXfGHsY/3R2Af26aWkFkWyhJCRZZdOOcqa+teKg4dYv//ALHNGva0xRiUhDKQlJVli5\neSdf+M2rB5U/e+NpjMzXcGuRbKUkJC1WRWUVbyzbxGNzVvHC4nUH7HvyuydTOLR3TJGJSFNREpIW\nxd35cH0ZT89bwx9f+/iAfZedmM/kMw6nYIAWlBNpLZSEJHbuzoKSUp5ftO6gxFPt3ZvP5pAeHZs5\nMhHJNCUhicXmHeW8s3wzN/5lXp11nvjuyRQe1ksvmYq0YkpC0iw2bN/D+yWlXPdQ8hnJH71uLJ8f\n2ksvmIq0EUpC0qT2VVaxcvNOijfs4KP1O/jtrI+S1v/xhKM59Yi+jBzUQ0spiLRBbTIJmdkE4L+B\nHOBP7n5HzCFlFXdnY1k5n27ZxcrNu/jhE++nfeyNZx3JKUf24YQhveiYm5PBKEUkG7S5JGRmOcAf\ngHOBEmCOmU139yXJj2w7yisq2bC9nLWle1i3fQ/rSnfzyxkfNOhcv/vaKEbm92BYny660xGRg7S5\nJASMAYrdfTmAmU0DJgKtJgm5O+UVVezeW8nufZXs2lvJnvBz194Ktu+pYOvOvWzZuZf/fnlZo9v7\n3rgjObJ/V47oF3065ekOR0TS0xaT0CBgVcJ2CXBSUzfyzsebueK+d1NXzBI/veAYhvTuzODw6d4x\nN+6QRKQVaItJqLY+IT+ggtlkYDLAkCFDGtRIeUVlg45rbrdeNIJenfPo1SWPXp1zOaR7R/p07UCO\nus5EpBm0xSRUAgxO2M4H1iRWcPcpwBSAwsLCAxJUusYe3oe3bhqHAdWvuRiGWciCduC2mR1UF6PW\n/dXHVbNazqUkIiLZoC0moTlAgZkNA1YDlwNXNnUjHXNzGNSzU1OfVkSkVWlzScjdK8zsRmAm0RDt\nB9x9ccxhiYi0SW0uCQG4+wxgRtxxiIi0dZobRUREYqMkJCIisVESEhGR2CgJiYhIbJSEREQkNube\noHcx2wwz2wisBHoApXVUq2tfbeU1y5JtJ37vC2xKO/C6JbuO+tTVNScvb8g1N9X11hVTQ+o11TXX\nta+tXHNL/nudbH9jrrnA3XukjMzd9UnjA0yp777aymuWJduu8b0o09eha473mpvqeutzzanqNdU1\n17WvrVxzS/573RzXnOyj7rj0PduAfbWV1yxLtp2szYaqzzl1zenty+ZrTlWvqa451X+PptCSr7kl\n/71Otj/j16zuuCxhZkXuXhh3HM2prV1zW7te0DWLBiZkkylxBxCDtnbNbe16Qdfc5ulOSEREYqM7\nIRERiY2SkIiIxEZJSEREYqMklOXM7BIzu8/MnjGz8XHH0xzM7HAzu9/Mnow7lkwysy5m9mD4870q\n7niaQ1v5s03UFn+HEykJxcjMHjCzDWa2qEb5BDP70MyKzeymZOdw96fd/Trgm8DXMhhuk2iia17u\n7tdmNtLMqOf1fxl4Mvz5XtzswTaR+lxzNv/ZJqrnNWfV73BTUxKK11RgQmKBmeUAfwDOB0YAV5jZ\nCDMbaWbP1fj0Tzj0lnBcSzeVprvmbDSVNK8fyAdWhWqVzRhjU5tK+tfcWkyl/tecLb/DTapNrqza\nUrj762Y2tEbxGKDY3ZcDmNk0YKK7/wr4Ys1zmJkBdwDPu/s/Mhtx4zXFNWez+lw/UEKUiOaTxf9g\nrOc1L2ne6DKjPtdsZkvJot/hppa1f7FbsUF89q9fiP5HNChJ/e8B5wCXmtl3MxlYBtXrms2sj5n9\nERhtZjdnOrhmUNf1/w34ipndQ2amfYlTrdfcCv9sE9X159wafocbTHdCLY/VUlbnG8XufhdwV+bC\naRb1vebNQGv6Za31+t19J3BNcwfTTOq65tb2Z5uormtuDb/DDaY7oZanBBicsJ0PrIkplubSFq85\nUVu8fl1z27jmlJSEWp45QIGZDTOzPOByYHrMMWVaW7zmRG3x+nXNbeOaU1ISipGZPQq8AxxtZiVm\ndq27VwA3AjOBpcDj7r44zjibUlu85kRt8fp1zW3jmhtKE5iKiEhsdCckIiKxURISEZHYKAmJiEhs\nlIRERCQ2SkIiIhIbJSEREYmNkpBII5hZpZnNT/gkXYaiuZjZCjNbaGaFZvZUiK3YzEoTYj2ljmO/\nbWYP1ygbEJYmyDWzx8xsi5ld0jxXI62Z3hMSaQQz2+HuXZv4nO3Di42NOccKoNDdNyWUnQn8yN2T\nzkxuZr2AZUC+u+8JZTcCI939O2H7/xGtdfR0Y+IU0Z2QSAaEO5H/MLN/hDuSY0J5l7Dg2Rwzm2dm\nE0P5N83sCTN7FnjRzNqZ2d1mtjisozTDzC41s7PN7KmEds41s781Is7Pm9lrZjbXzJ43swHuvhV4\nG7gwoerlwKMNbUekLkpCIo3TqUZ3XOLKmJvc/QTgHuBHoezfgL+7++eBs4DfmFmXsO9kYJK7jyNa\nVXUoMBL4dtgH8HdguJn1C9vXAH9uSOBm1gH4b+Ar7n4i8P+A28PuR4kSD2Y2OMTyekPaEUlGSzmI\nNM5udx9Vx77qO5S5REkFYDxwsZlVJ6WOwJDwfZa7bwnfTwOecPcqYJ2ZvQLRvP/hec3XzezPRMnp\n6gbGPhw4FngpWhuRHKKZniGaWPMuM+tKtOT04yEWkSalJCSSOeXhZyWf/a4Z0Z3Hh4kVzewkYGdi\nUZLz/plokbs9RImqoc+PDFjg7qfX3OHuO83sJaLVTi8Hrm9gGyJJqTtOpHnNBL4XlmXHzEbXUe9N\nolVV25nZAODM6h3uvoZoHZpbgKmNiGUJ0WqmY0IseWZ2bML+R4F/BXq6+5xGtCNSJyUhkcap+Uzo\njhT1bwdygQVmtojPnsHU9FeirrFFwL3AbKA0Yf8jwCp3X9LQwN29HLgU+K2ZvQ/MA05KqPICUVfh\ntIa2IZKKhmiLtFBm1tXdd5hZH+A94FR3Xxf2/R6Y5+7313HsCmoM0W7i2DREW5qE7oREWq7nzGw+\n8AZwe0ICmgt8jmg0W102Ai+bWWFTB2VmjwGnEj2TEmkU3QmJiEhsdCckIiKxURISEZHYKAmJiEhs\nlIRERCQ2SkIiIhIbJSEREYnN/wfIgavTd8YMBAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "aeff = cta_irf.aeff.to_effective_area_table(offset=offset)\n", "aeff.plot()\n", "print(cta_irf.aeff.data)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/deil/Library/Python/3.6/lib/python/site-packages/astropy-2.0.dev18336-py3.6-macosx-10.12-x86_64.egg/astropy/units/quantity.py:1005: RuntimeWarning: invalid value encountered in true_divide\n", " return super(Quantity, self).__truediv__(other)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "NDDataArray summary info\n", "e_true : size = 500, min = 0.005 TeV, max = 495.450 TeV\n", "e_reco : size = 500, min = 0.005 TeV, max = 495.450 TeV\n", "Data : size = 250000, min = 0.000, max = 0.165\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEOCAYAAACetPCkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xlw3Od95/n387v67sZ9gyAAHiJFHbTo+4wd2Y6P+Iyv\ncTKbZKOaqmRmtrZmd2sm2fFWKrPJbmV2Kll7K+PyNfY4PuI4jjWOY8vOyLIlS6JIUQfvmwBxA40+\nf92/69k/ugFBMkWDEkFc35eLRXYD6H5gUvzweb7P832U1hohhBBitYz1HoAQQojNRYJDCCHEDZHg\nEEIIcUMkOIQQQtwQCQ4hhBA3RIJDCCHEDZHgEEIIcUMkOIQQQtwQCQ4hhBA3RIJDCCHEDbHWewA3\nk1LqvcB7M5n07+3es2u9hyOEEJvG0SPH5rTWnav5XLUVe1Xdc+igfvixB9d7GEIIsWkkrJYjWutD\nq/lcWaoSQghxQyQ4hBBC3JAtFRxKqfcqpT67uFhY76EIIcSWtaWCQ2t9v9b6vpaW3HoPRQghtqwt\nFRxCCCHWngSHEEKIGyLBIYQQG9RGPS6xpYJDiuNCiK1Ca01EtN7DuKYtFRxSHBdCbAWhDgl1iBd6\n6z2Ua9pSwSGEEJuZ1roRGlGIH/nUwvp6D+maJDiEEGID0FoTEhHqiHrkUYs8yn55vYd1TRIcQgix\nziIdERIRRAFe6OEGNcp+mYK3uN5Du6Yt1R1XCCE2k0hHaDRBFOBHAW5Yox7WKHgFrlZKnCk6HGhb\n71H+IplxCCHEOliqZwRRQD30qIV1il6BolfgbDHPuZLNExOt6z3Ma5IZhxBC3EJL22y11vhRQKAD\nykGZWuByoTSLGwT8w4Ve/MigUN6Y5zi2VHAsXeQ0Mjq83kMRQohfEOlGHSPQIYEOKHll3LDK6cIM\nJV/xtyf78L2QwqLLpfML9A9uzKMFW2qpSs5xCCE2qkZYhHjNbbYVv0rRL3C5PM+FUpwHLnVTyLv4\nfsjk1SIju9pp70yt97CvaUvNOIQQYqNZWpoKooAgCigHFephjVpY40Ipz1jF4chEKzPzPjNTJXw/\nIpuL0d6VZl93ab2Hf00SHEIIsQZWBkaoo2Zg1BkrT1MNAsYriscnu5jMG0xdLeB5IdlcnJHhNAe6\nC/QmJ9iVia33t3FNEhxCCLEGltqGuGGdWlijUM8zXi3wTN7h0Su9uK5Pfr7KuTPz9PZlWJivsmdf\nB51pl9FMhe5EnB3p3vX+Nq5JgkMIIW6iSEdEaPzIx48CSn6JWuByvrTAlYrD4+NtzM2UsWyTq2MF\nRne30dGV5tDBFl7VO0NfMmA020nWzpGyk+v97VyTBIcQQtwkkY4ImjONSuBSC12ulKcp+R5/e7aH\nqmcyN1vl8sU8QRCRTDl0dKW5c9BlKFPhrjaHrN1CW7ydmOEQM531/pauSYJDCCFeJq01gQ6IdNQM\njBpXK1OUfJ+/PddKxbOYmvG4eH6K7p4MqbTD6N5OXjFYYjA1y51tCTJ2O22xNmzDJm7GUEphKnO9\nv7VrkuAQQoiXIdIRkY7wowAv8in5JQpegUtll2NzLUzmLZShmBgvMDzaTmdHjJ6sTX8mzyvaAjri\naToT3TiGTcJKYCgDUxko1Hp/ay9qwweHUur9wLuBLuAzWusfrvOQhBACeO5An988zFePalwuzXC1\nGvCjyz1cnfI5f2YCzw9Jp2PEEza3dZUYTFUZSCkG0h1k7CwpK4llWFjKRCm1oUMD1ik4lFJfAN4D\nzGitD6x4/p3AXwAm8Dmt9Z9prb8DfEcp1Qr8OSDBIYRYVytPgBf9EvWwzuXSNGPVkAcudjMxVef0\nyXECPyKZtLn9zh5euctnODPHKzoSZOwu0nbmeYFhqM1zHnu9ZhxfAj4NfHnpCaWUCXwGuBcYBw4r\npb6rtT7R/JQ/an5cCCHWhUYTNpemqs1DfBOVKYq+xzfPtJJ3bSYnSpw5OUdnV5J8vsbrXt/HUGuV\nN/aU6Uok6U32Yhs2jmHjmA4GalOFBqxTcGitH1JK7XzB068CzmmtLwAopb4OvE8pdRL4M+D7Wuuj\nt3SgQgixgtaaKArxdUAtqFHyi0y5FZ6czzC1aBKGAVcuLrJzuIXuviy9rRFvGpyiP2mzI91LykqR\nsVOYRuOvXkuZaDZmI8Pr2Ug1jn5gbMXjceDVwL8EfhXIKaV2aa3/6lpfrJS6D7gPYHDH4BoPVQix\n3fiRTxAFlPwKRb/AmcIMj0y38NjJVjwv5OmjZ+kbyDE82kZ3b4q3D0+zIw2j2T7Sdoa0lcRsLkut\ntNHrGdeykYLjWv/vaa31XwJ/+cu+WGv9WeCzAPccOrj5IlwIseEs35mhQypBlXpYY9qdZaJa4aGJ\nVk5PxJiZmmNxsUZ3b4bBoRbuHo0Yys5xe6tNzkmRtFLYhrVht9a+FBspOMaBlVOFAWDiRl5A2qoL\nIW6WIAqWaxn1sM60O03Rq/PtCxnK9TSXrvqcOz1BS2sc0zA4cHcvfS0eh7ry9CZjDKR6sQyLtJ1a\nLn4rtflmF9eykYLjMLBbKTUMXAU+BnziRl5Aa30/cP89hw7+3hqMTwixDWh0Y/lIKdCNJSo3qDLt\nVnh0NsXTZ0Nct8aJh8+Q3dFNMuWwd38n79szwVA6xlB6JykrRdJKYG6y3VKrtV7bcb8GvAXoUEqN\nA5/SWn9eKfUHwA9obMf9gtb6+A2+rsw4hBAviUYvH+TzI59Fr0DRW+TYQplHJ9o5ejwBRJx4+io7\nhtt47XvuYqDL5O07pxlIBQymdpCy0yTNOIZh/kItYytRWm+9csA9hw7qhx97cL2HIYTYJJZ7TEUB\nRb9MLawxXpllrOLx0FgnFycjLl9YwK36tLTGGdnTyauHCgylq9zZ1kLSStEaa1neZmuw+ZalElbL\nEa31odV87kZaqhJCiFtqqcdUqCMqQZVaUONKZYp8vcY3T/dSrinmZsqcOz1PZ1eSmhtw+x2d7Gyr\ncnebS1ciRUe8E8uwiBsxDMPclKFxo7ZUcMhSlRDihijQujHjqPhlaqHLP03YjBWynD5d4NTxKbp6\nswC89lAreztKvKW3SM7O0hYfIG7GsQ0LNFs+LFbaUsEhxXEhxGr5kU+gQ0p+mbJf4kJxlifmEjxy\nJoFX9zl7dp6do+109WTo7Y7xuv5pdqRM2pwO0k6GmOE0dkqhrn2YYAvbUsEhhBDXo9GgISSiHnpU\nQ5dZd4Yr5SJ/d66Ls5dqHHvsHIQBbb2tjOzp5J17ZhnNFNmV7SNhJcnZGQzDxNwGS1IvZksFhyxV\nCSGuR2tNPfLwQ5/Z2ixuUOWbFwzmqx2cuehy+sQsfUNtFIs1fvWtvexpX+TOVou2WJqskyNmOljN\ndiHbNTRgiwWHLFUJIV6MF3p40XOB8cXTJuOLWY4/Pc3Y08fJDvXgeSFveG0H+zuL3NtXI+N00B5r\nwzIsbMPalO1B1sKWCg4hhHihpfsyil6Jolfg5OI0D17t5OdPutTrZcZOjNOzb5DRPZ0MdWvePTxN\nfypNR7yLhBXHMZ3GspSExrItFRyyVCWEWCnUIV7kUwvrLHp5xsqLfO9iLyfOVDh15BJ4HpkdPezd\n382v3z7HzjTsyg6QsJK0Ornl4vd2Xpa6li0VHLJUJYSAxvbaeuRRDz1m3CnGKot8+XgXp087nP75\nYfB9MruGeOVrd/CxAxMMp4sMpIaIW3GydhpDGc+1HhG/YEsFhxBCLJ0An6/NUwtdvn7B5+JCF0eP\nzDD+5Fla9+ykXg/4zQ90cGf7JK/s6CNlp8nZGUzjud5SEhovToJDCLElhDok1BFlv0LJL3JycYKf\nTLTw48c86rUq409doOvACLff2cvunhq/NlihK9FGLtZC3IxhG7YsSa3SlgoOqXEIsf1EOiLUIeWg\nStkrcaY4wU+nUjzweJJKuczlR0+S2bOT17zvHm4b1LxvZI7eRIauxChxK758kE+s3pYKDqlxCLF9\naK2JiKiFdbzQY8adZsIt8p3zbZy+HHLqxBSR65LePcSBu3p41x1l9uY8RjP9JKwkKSvZ6C0loXHD\ntlRwCCG2j5CIIAooeEVqQZV/vOpyuZDjyDMVTj41QUdvjlrd5hPv7+L2tjyv624jbWVpjeWueYWr\nWD0JDiHEphLpqHEzX+BS8stcLk/y2KzBPz6VoLjocuqnJ0gP99O/I8fwUIo3900ylM7S6rSRtBJY\nhoWBzDJeDgkOIcSmEeoQN6hRCapcKV/l0VnN3x9OsJh3Of+zx0nvHuI17z/Evh2aD42W6U7E6Iof\nIGY6xMzYeg9/y5DgEEJseFprvOa5jLnaLJPuIl85lePZU1WO//gIuFXYMcze2zr54MECt+UUo5lB\nklaCpJ3ElBnGTbWlgkN2VQmx9Sy1DFmoLeCGVb59qcpYIcuRp4qce/IyrXuHqdUCPvL+Pva3F3lj\nd5q0nSHrZHAMG0uZja644qbZUsEhu6qE2DqeC4w8bljlO5fLjBeTPPpkxOmnTtHW20qsLcevv7OL\nOzoKvLFHk7b7aY+1YipzuYstyGG+m21LBYcQYvNbWfwuB43i9+E5gx89k6RcqnH60bO0Dvcxsrud\nkUGH9+ycoz+VXW5KaBu2bLFdYxIcQogNJdAhXuiR9xaZdRf41vk0T50JOfqDo1Apw8AgI7va+Mir\nytyWU+zO7iRhxUnbqW1x3/dGIMEhhNgQ6mEdPwqYrc0yV1vkG+cSHDvjcOTHz0C5BH0D3PPOu/jQ\nq1xuaymzN9foYpux09v6Nr71IMEhhFhXSy1DCn6JelDjRxNzXCqmeOiJKmcfO0e8p5NaJss///gQ\nB7vyvLG7g7SdpSWWw1QGphzku+UkOIQQ6ybSjZYhlaDKeHmCx2YDvvloksWFCpceeork3l0cfOUA\nd+30+NCwS1usm454B47p4Bj2eg9/29pSwSHbcYXYPPzIx4t8FmoLzNUX+JsLDkfPxjj2kxNQKsLg\nELff2cOH71pgT85kILWDhBUnYcXl5Pc621LBIdtxhdj4lmoZc7U5Zmt5vnUhzpNnYxz+p1NQLGL0\n9XHPr+7nA/dU2Jsrsjc3QNJOkrFSGMqQHVMbwJYKDiHExhXqkDAKWfSKuEGV74/Pc7mY5CeHa5w+\ncolERyuuE+NDHxjmzu4Cb+1NkbaztMfbsF5wLkOsL/mdEEKsuUhH1EOPSlBlsjrFkbk63zmaYX62\nwvkfP4G9axd79nVw22icf7Y3T28yR1e8m7gZwzGd9R6+eIHrBodSygHeBbwR6ANc4FngH7TWp9Z+\neEKIzWypKWE5qDBWnuCxWc3fPR4nv+Bx8cGfEt+zi3s+/Eb2jdh8dHeBnmSMrsQe4obcyLeRvWhw\nKKX+CPgQ8BBwBHgAiAN7gP+kGr+j/0Zr/eytGKgQYnNZLn7XF5itLfC1s0mOnPR5+ruPNA7yDY9y\n+x3d/NbrSuzJ2YxkhkhY8cYFS1LH2NCuN+N4Rmv9Jy/ysf9bKdULDK7BmIQQm9jSuYylWsbfX5nj\n7EKaH/80z8Tjz+KMjuB5Iff95gAHOxd5bVcfKStNzslgKUtmGZvA9YJDKaVsrbV/rQ9qrSeBybUZ\nlhBis1k+yOeVmsXvac4XUjzwcMTFx58guaMPe+dOfvujvdzduchrOhOk7S5yTgZTyRWum8n1guN3\ngf+slPoe8DXgx1rr6NYMSwixmaw8yDdRneDIXI2/f7qVQt7l4uGztO3Zwb47erh9h88Hhmt0JhoH\n+WKmgy0H+TadFw0OrfV7lVItNOoc/yvwJaXUt4Gvaa0fuVUDFEJsXEuzDD8KyNfz5OsL/Lcriicu\ndfDIT85BsUBioJd9d/TwsYNz7M7GGM4MyEG+Te66u6q01ovA54HPK6W6gN8A/kopldFa35Lj2Uqp\nEeAPgZzW+sO34j2FEL9cpCPcsHGNaz2s808TV7lcTvHDxyOOP3GceC5NLUjxjncOcVfPIvd0ZMnY\nWTJOGltZ0mNqE1vVOQ6lVA54N/A+oB34u5fzpkqpLwDvAWa01gdWPP9O4C8AE/ic1vrPtNYXgN9V\nSn3r5bynEOLmCHVILaxTDVymq5McW6hwpZLkgaNpjh2ZIJm0wTB5+zt2ckdPgV8bqJC1O+iIt2MZ\nlixNbQHX246bpBEUHwdeDXwP+HPgRzeh1vEl4NPAl1e8nwl8BrgXGAcOK6W+q7U+8TLfSwhxk3iR\nj9/sL5X3Fvn7K3D0Uhf5hSpP/fAp4gO97BhqYWhkhA/sucLOdJbuRC9JK4FjOlIA3yKuN+O4AvwY\n+CLwYa21d7PeVGv9kFJq5wuefhVwrjnDQCn1dRrBtargUErdB9wHMLhDdgkLcTMt3cpX8svUwxqP\nzo5xsZTix0dNnnrsFMzNQDbHzuFWPvQmuL3lKne07SRlpRr3ZShTttluIdcLjp1a6zI0TpArpXZp\nrc+t4Vj6gbEVj8eBVyul2oH/ABxUSv1brfWfXuuLtdafBT4LcM+hg3IzvRA3QahDQh2x6BWoBS7/\nNDHFlUqSHz2R4Zkj48TTcYhC9t97kPe+VnGgtcQdbYOkrRRZ2Wa7ZV1vV9VSaLwb+H8ABxhWSt0N\nfEpr/YGbPJZr/XNEa63ngX+xqheQtupC3DRLW2yrgctkdZIjczW+/VQb+fkKJ354FGfnECO72tg1\nupO7exd550CMrL2T9ni7NCXc4lbzO/vHNGoc/x1Aa31MKbVrDcYyzvNPog8AEzfyAtJWXYibI4gC\n/Mhnsb7Iopfn2xdNDp9r4ef3H4XFBejoYnRXG7/xxojbW6YZTOXoSfYRN2PEDEeWpba41QSHr7Ve\nfMEfhLVYCjoM7FZKDQNXgY8Bn1iD9xFCvAg/8pfPZMzXF/j7K3D4XIKffu8ZmJ+D9g72v+e1fPgN\nEftaXA60DpKyUsStOHEjJoGxTawmOE4qpT4CGM2/1P818OjLeVOl1NeAtwAdSqlxGktfn1dK/QHw\nAxrbcb+gtT5+g68rS1VCvESBDikHVWpBjZ9NX+F8KcX3H4k4+fjTOG0teLE4//w3R7m7c5E39bST\nsvtpdXJYhiV1jG1GaX39yYNSKgX8e+Dtzad+APyx1rq6xmN7ye45dFA//NiD6z0MITaFxl0Zdaph\njanqJE/Ol/j6kx3MTZc4/YPHMUd3sW9/J/t2p/jNfc27MhLdxM04tjQl3DISVssRrfWh1Xzu9c5x\n/J9a63+nta4A/1vzx4YmMw4hbkygQ/zIJ+8VKHh57r8S8tjFLh7+3jFYmIeuXvbf3sXH3+Cxr2WB\n3dkdJO0EKSspJ7+3sRedcSiljmqtX3GLx3NTyIxDiOsLdEgYBeTrBRa9Bb5z2eeRs0l++reHIT8P\nXT0ceMMefuvNdfa3GOzKDpK0EnJXxhZ2U2YcgKmUauXa22TRWi+8lMEJIdaP1pqQiKJXohbW+Mnk\nJS6U0tz/s5DTjx5FdXSi43E+9Jt3cU9vnl/pbSFtZ2mNtWAqQ0JDANcPjtto3Px3zfMVwMiajOhl\nkKUqIV7cci0jcJlyJ3lirsI3n+xgfqbM6R8+gT06wr4D3dw2OszH987Tt+Leb6lliJWut1T1pNb6\n4C0ez00hS1VCPJ8b1qgFNabcSY7NF/ivRzqYmSpx7h8egV17ufNgLwd2x/jkbSW6Ey10xbtwTAdH\nGhJuGzdrqUoIscktHeRbqOdZ9PL87SV4/Fwnj3z/mUbxe2iEu17Ry8deV2N/zmU4vYOklZC7MsR1\nXS84/uKWjUIIcVMtNSUs+CVqQY0Hp8a4VErxw8c1x584SbKzlWoqzSc/upO7uxZ5c28raTtDm9OC\naZiyY0pc1/WC4/VKqSNa62de+IHm2Y6PAnWt9VfXbHQ3SGocYrtb2i21UF/EDavcf6XA5UKShx63\nOX/0NLHONlCKe9/euFzpHf0eWad3ub+U3JUhVuN6wfEZ4H9XSt0BPAvMAnFgN5AFvgBsmNAA6VUl\ntq+l3VJu4FL2K1wpT3J4TvP9p1IU8i7nf3ac5MgORne3MzqS5iN7ptiRztKdbNyVETPkrgyxetfr\njnsM+IhSKg0cAnoBFziptT59i8YnhFgFP/LxIp+F+gJztTzfvRzjyQsWj/34FJQK0NHJrr0dfOA1\nPvtz0+xvHSRppcjaaQzZZitu0C8tjjfbqz+49kMRQtyolUtTtdDle2N5LheTPHy0xsknz5PuylFO\nJPjQB0c50F3g7X0OGaeT9lgblmFJ63PxkmypPzVS4xDbxYsFxk+fUJw+eopkRwvYNu+4d5A7ugu8\nvb9Gxu6kPd6OrSQwxMuzpf70SI1DbHVLtYxqUKXsV7hUmuCxWYPvHYmzmK9w/ifHiI+OsGtvB7uG\nk3xkzwyDzXu/E1ZcahnipvilwaGUOqC1fvZWDEYI8eIiHVGPPLzQY7Y2w2S1wNfP5Dh+rs5TD52C\nchn6d7D/QDe//so6+3Iz7GsdIGklydgpucZV3DSrmXH8lVLKAb4E/LXWenFthySEWGkpMOqhx4w7\nxXilwJdPdHLipMHZf3oIXBf6d3DHOw7yodcF7G8psq9lQIrfYs2spjj+BqXUbuB3gCeUUo8DX9Ra\nP7DmoxNim/MinyAKmKvNUQkqfPkMXFro5ujhq8wcPo45MkJYqvD+T9zBXT2L3LtU/JZ7v8UaWtWf\nKq31WaXUHwFPAH8JHFSNjmf/Tmv97bUc4I2Q4rjYSiIdEUQBlaDKQn2Rw3N1jp7NUMgvMPPsBdru\n3sfdrxzg7sESH9hRJud00B5vxzFsCQyxplZzA+CdwG8D7wYeAD6vtT6qlOoDfq61Hlr7Yd4YaXIo\nNrNIR7hhjUpQZax8lZ/PRHzrUYfFBZfLDx0jtW8Pdx/q544hn98Y8WiPt9AR78AxbDn5LV6ym93k\n8NPA52jMLtylJ7XWE81ZiBDiJvFCj3rkMVubZaqa58uncjx9vMzJB58Et4qzezd3H+rnk/fMsjub\nYDA9QMKMkzDjUscQt8xqahxvus7HvnJzhyPE9hTqkDAKma/nccMqf3OxxIWFFn76yCyTh48THx3B\n80J+6yN9vKJzlld39pCy0+ScLJbslhK32Gq24z5D4+KmlQo06h1/orWeX4uBCbFdBDpsnssoc744\nwcPTDt99WFMqzjH56DOk9u3hzlf0cvtQxIdHXDoT3XTEO4jJBUtinaxmqer7QAj8dfPxx5o/F2ls\n0X3vzR+WEFub1ppQh1SCKiW/xJnCBA9NJfmHhw0qZZeJwydo2b+Lt/zO27h9oMYHd7q0xXJ0xAfk\ngiWx7lYTHK/XWr9+xeNnlFIPa61fr5T65FoNTIitauVBvhl3hvFqnr8+1cGJU0XOPnoG6h6p3SPc\ndU8/H71zmtFMnIHUAAkrLhcsiQ1hNcGRVkq9Wmv9GIBS6lVAuvmxYM1GJsQW5DfPZczW5nCDKl89\nFzJW6OTYsTkuPTNG12078L2IT7yvjQNt07yqo4uUnaHFyWLIBUtig1hNcPwu8MVme3WAEvC7zcuc\n/nTNRvYSyDkOsVF5kY8f+cy6s7hhlf9yWjFeSPD0sRnGnnqGluF+4h2tvPtt7RzoKPCGbpO0PURb\nrE0O8okN57p/GpVSBjCitb5DKZWjce5jZcuRb67p6G6QNDkUG83SbqlSUKHkFTm5OMlDky385Ik6\ntVqZsSfP075nB3tv72JPf8R7hwv0JbJ0xruIW3Ecw5YdU2LDuW5waK0jpdQfAN/UWhdu0ZiE2PSW\nutjWghr1sM6UO8nlcpFvnOrl7NkCp45eAa9OZrif/Xd089478uzOwu5sPwkrSdpOSY8psWGtZv77\ngFLq3wDfACpLT2qtF9ZsVEJsUiuvcK2Fdaark1wqF/nrkwNcOOtz+sgxCAIyg90cPNTP+w7MsyuT\nZzTbR0K62IpNYjXB8TvNn39/xXMaGLn5wxFi84p0tHyF64w7QyWo8pUzNlPFfs6enufc8Sn69vRR\nc30++K4O9rXN87quDGk7Q1ustXEjn5zLEJvAak6OS6VZiF8iiAICHVL0ShT9RZ7NT/PQ1U5+/qSL\n71c4/9QVenb1sPu2Tka6fX5tR4HeZIbOeBcJK4FjOhgYEhpiU1jNyfEk8D8DO7TW9zVbrO/VWv+3\nNR+dEBtcqEMiHVEJXGpBjfHKJBfLVf7m5A6uXJjn9NNXIQho6e9g995O3n5bnpGMt7w0lbZTmIaF\noQwUEhpic1jNUtUXgSPA65qPx4G/ASQ4xLYV6pBQR1T8Cm7ocrk8xYWSx7ee6eLSeZ9Lzx6BICA9\n0MWdB/t454ESuzKz7GvpJWmlyDkZTGViGqYEhth0VhMco1rrjyqlPg6gtXaVzKfFNhbqkHroLe+W\nWvRqfP7ZFuaLcOHsHFdOXKV7tIdq1ef97+rhtvYib+p2yDjttMfasQ0Lx3BAIaEhNqXVBIenlErQ\nbHSolBoF6ms6qhWaBw3/P8ADHtRaf/VWvbcQLxTokDAKqAYuJb/IxdIij0y3cO5SlcCPuHJpkZE7\nd7BnXwe7O8q8Y7BEZzzdOJdhxoiZMdkxJTa91QTHp4B/BAaVUl8FXg/8Dy/nTZVSXwDeA8xorQ+s\neP6dwF8AJvA5rfWfAR8EvqW1vl8p9Q1AgkPccpGOCHVI2a80mxJO8shMih89kaFWq3H2ySv07urh\nbb+2h73dVd4zNE9nPENnfICYGSNmOFL4FlvGanZVPaCUOgq8BlDAv9Zaz73M9/0SjQuivrz0hFLK\nBD4D3EujjnJYKfVdYAB4pvlp4ct8XyFuWKhD3OZBvpnaNJdKBb5xuoez54qcOHwZwoD2oW723dHD\nB/ZNMZy2GUr3kbSSpKxko/AtoSG2kNU2wIkD+ebn71dKobV+6KW+qdb6IaXUzhc8/SrgnNb6AoBS\n6uvA+2iEyABwDKQtqLi1vNDDi3zmanNUgwpfOWswWezh6afnufjMGJ3D3biuzwd+rZ19bTMc6mgl\nbWVojbVgK0t6TIktaTXbcf8v4KPAcSBqPq2BlxwcL6IfGFvxeBx4NfCXwKeVUu8G7r/OOO8D7gMY\n3DF4k4cmtpOlOoYfBcy407ihy5fPmEwVkxw/vsD54+foHGwn1dXGu97WwW3tRd7co0nbfXTE2jAN\nSy5YEltTIoDYAAAgAElEQVTaav459H4a5zbWuiB+rf/KtNa6Avz2L/tirfVngc8C3HPo4AtvLBTi\nl1rZX6oW1il6i5xYnOInE+389AkXz6tw8cmLdO3p57bbu9nd6/OenQV6Exk6E90kpPgttonVBMcF\nwGbtd1KNAyunCgPAxI28gLRVFy9VpCO8yMeLPGbcaSqBy8VSkW+cHOTcmXnOHbsCvk9msJu9+7p4\n9/48o9mQ0Ww/STPROMgnTQnFNrGa4KgCx5RSP2ZFeGit/9VNHsthYLdSahi4SuOK2k/cyAtIW3Xx\nUviRT6BDCl6RolfgeH6Wh652MF1uhsaJaQb29FKt+rzvnZ3c1rbIG7oTpO0M7c37MmRpSmwnqwmO\n7zZ/3DRKqa8BbwE6lFLjwKe01p9vtnD/AY3tuF/QWh+/wdeVGYdYtSAK8HXQDIxFnl6Y52cTHTz6\nZJIwqnDu1EUGd7bR2ZfjPW/Ncltrkdd0alJ2T+OCJalliG1Kaf3LywHNA4A7tNan135IL989hw7q\nhx97cL2HITao5XYhQRU3cBkrT3GxXONbJ3q5dG6e88cnIAjo3d3bbEoY8O6di/Qm03QluklaCWJm\nDFOaEootJGG1HNFaH1rN565mV9V7gT8HHGBYKXU38Mda619/ecMU4tZbbhcSeUxVJyl4Ll88kWO2\n0MLFc7NcOjVF30gXNddn7/4u3rVvnpFMxO5so8dU1k437/6WpoRi+1rNUtX/QeOMxYMAWutjzTrE\nhiNLVeLFhDrECz1qkcd0dZKSX+OLJzLMlpNcPDvH+ZPTDI520DnYznvemmVPS5nhzBy7sj0SGEK8\nwGqCI9BaF14wJd+Q212lOC5eSGuNrwP8yKfgFSh4ixybL/CT8U6OPVMi0przZ2bZuaeTXbd1sauz\nyr0DJXqTaVJWmvZ4O7aysA1bmhIK0bSa4HhWKfUJwGzexfGvgEfWdlhCvHyBDomikHKzlnGlMsWl\nUp2/O93P2MUFzpyeR0cRPf05dt3Wxdt3zTCSiRjN9pKwEiTNBDHDkZYhQrzAaoLjXwJ/SGMr7l/T\n2PX0J2s5qJdKlqoEvHhgfPtUH2MXFzh96jREIV29WQ7c3cc7ds8xnJlkT66XhJUkZ2eW7/2WcxlC\n/KJV7arabGRX1fYU6YhIR9TC+oq7Mup85VSOhYrD+KUFzp2dZ3BHC9Wqx7vf2spoS4XXddtk7Qwd\n8Q4sw2p0skXJLENsKzd1V5UQm0Gko+VzGUW/SMEr8Gy+xMMT7Tx7qkoUljl5Ypadwy0MjbSxs9Pj\nVwcX6Eum6U50k7KSxJvtQmSWIcT1bangkKWq7UejCXX0vKWp8co0l8o17j/Xy5XLJc6dmSMINT29\naUb3dvK20XmGMx57cz2NpSkni6kMTGXKLEOIVdhSwSG7qrYPjSaMQgIdUvLL1AKXS+VpLpZC7j/T\nw9XLec6cOU8YRHR2pdh/Rw+/sqfASOYq+1q6SVopWp2cdLIV4iVYzQHAAeD/Bd5Ao636z2hc5jS+\nxmMT4ppCHaK1xg0bXWwnq5MUvDpfPd3CYtVhYizfqGUMZimXPN7xK52MtBZ4XRe0xrrojHdiGzaO\n6WDI6W8hbthqZhxfpLGb6jeajz/ZfO7etRqUENei0ctdbP3QZ9FbZKG+yLH5Gj8b7+DEqSLKUJw6\nMcvO4VZ2DLcy2BHwq4PT9CXT9KW6SFkpEmYc0zAxlbne35IQm9JqgqNTa/3FFY+/pJT6n9ZqQC+H\n1Di2Jq01Ec8FRt5bZL6W59hCnZ+OdXLqVIRSRc6dnWdouJVfefsudrTVedvgFL2JFD3JUdJWirgV\nX65lCCFeutUEx5xS6pPA15qPPw7Mr92QXjqpcWw9WmsCHTRqGV4ZN6xyoTTNxRL8w5kuJscLXLyw\nQBRpOjpT7Bhu480jCwynXfbmuklaSVqcHJYysZpnM4QQL89qguN3gE8D/4lGq5FHms8JsaaCKCDS\nEZXQpR7UGa80zmV8/XQnxZrF3EyRsSsFRne3Uy7Vedvrcwzl8ryqQ9Ea66A93oGtrGYtQ0loCHGT\n/NLg0FpfAaQTrrhllgKjGtaoBTWuVicpeB7fPNtGsdbK7EyVC2cn6BvI0t6e4NV3xhnKRby6s0Jr\nLEVXvAvbsImbMZRSsjQlxE22ml1V/4XGLqrF5uNW4D9qrWXWIW6qpZPfjStcfQrN4vepQo0nplq5\nPBmhVJ1L5xcY3NlCT2+W3qzL63rn6U3G6E12krbTxK04ljLlXIYQa2Q1S1V3LoUGgNY6r5Q6uIZj\nEtvM0vZaL/LxI59Fr0A9dDm1OMvlSpL/frabuZkyly7miUJNOuOQzca5rbvIYMplJJMi6yRJ22kc\nw8ZSJkoaEwqxZlYTHIZSqlVrnQdQSrWt8utuOdlVtfmEOiSIAkIdUfYruKHLxdI0Zd/j787143om\nlXKVqckSQztbKCzWeM0r2+nLuby6s0RnPEFvshvHsElZSQxlYCkL6X4uxNpZTQD8R+ARpdS3aBTH\nPwL8hzUd1Usku6o2j+cFRlChFrhcLs9Q8j2+c66bqm+ykPe4eG6Kzq408bjF3l1p+rImr+qaojOe\nZDC9E9uwSdspzOauKaWU3JkhxBpbTXH8y0qpJ4C30vh33Ae11ifWfGRiS9JaN+/8DptdbD0WvTzT\nbpFTBc2RqS7GpkNsS3Ph7Dy9/Vla25LsGU1xT2+e/qTLaLadrJ0jaSWJGTa2sqQ5oRC30GqXnNqA\nitb6i0qpTqXUsNb64loOTGw9L6xlFLwCtbDG+eI841WLh8c6mJn3uXplEd8PSSQsWtuS7O936U+5\n3NUGrbE22mMdOKZD0oxjSHNCIW651eyq+hRwCNhLo9WIDfxX4PVrOzSxVSwFhh/5+FGj7Xk9rHO5\nPM+kG/HI1Q6mCjZjlxbw6iHxhMXe27vZ31OmJ7nAnW3Q6sRpi3cQMxySVuJ5tQxZmhLi1lrNjOMD\nwEHgKIDWekIplVnTUYktYalVyPJBvsClFrpcrcxSDX1+MJbB9U3mKzGuXp6nrT3FwnyFg3fk6E6X\nOdiRpyMWZ0e6B1vZZOx04wS4YUlzQiHW0WqCw9Naa6WUBlBKpdZ4TGILWGoVEuqIelinHnmUvSIL\n9RJjFZdnFnJMl+KYhubKxQW6+7KkUhYj/Sb7Ogv0JGoMZ7K0ODniZpy4GVs+m2FgyK4pIdbRaoLj\nm0qp/wy0KKV+j0a7kc+t7bBeGtmOu740uhkYjeJ3PfTwQo+iv8h8rcyVisvxfJaxxR5Mpbl4uUJ7\nR4rdu3Nk4wF72ufpitfZnU3SEushZaVJmHFs037uQJ9c6SrEulvNrqo/V0rdCxRp1Dn+vdb6gTUf\n2Usg23HXT6QjAh0SNWcYXuRT9ArUwxrnS4tMVmM8MdnNzKJi7FIerx4QhBGxmEUmGSNmheRsn6wN\nKTuO2ewxtdQyRCmjUcuQzBBi3a1qV1UzKB4AUEqZSql/prX+6pqOTGwKS9trIx1RjzwCHVL2y9RD\nl0ulBdww4MdjnXiBQaluMTOZp70zxex0ib0HumlJRuzpKNKdqHFPRwpLWbQ4rTiGTdxwMJWJ2QwN\nmWkIsTG8aHAopbLA7wP9wHdpBMfvA/8LcAyQ4NjmnrdbSge4QQ0vrFP0i0xWy1wuGxyf72CmHMMx\nIy5fyNPRlcayTbq64gy1VulM1tmbrbEjnSFlpYlbCeJmDNuwG0VwZUhNQ4gN5nozjq8AeeDnwP9I\nIzAc4H1a62O3YGxig/plgXGloji5kGOqGMexIsYuF2jvTNM/1ErSiRhqrdKW8NibrTGYSpOye0hY\nSeJmHMewsQ0LpQzMZmDIdlshNpbrBceI1voOAKXU54A5YIfWunRLRiY2nEhHRGjC6Pm7pUpeES/y\nuFQuMOtaPDuXY77iMDXtUq8FaA22Y9KV8WhN+uzKlWmPh+zMtBAzYmTsLKZhEjdjzx3ok6UpITas\n6wWHv/QLrXWolLooobE9aa3RPL+WEUbNWkZUY6ySpx6G/GyylUhD1TeZmffJtSaYnS6zdzRFJlZn\nuKVCe6zOnW0pbMOkxWnDNmwSVrx5oK/Za0r6TQmxoV0vOO5SShWbv1ZAovlYAVprnV3z0Yl1tTIw\nlmYYjS62JephfTkwHp7OoTVUPJNLUwaptI1fr9DWbTJ4m81QS4E2x+Ou9gSWytAaa8UybJJmAlMZ\nz9Uymr2mJDSE2NheNDi01nJt2jamtSakcbGSH/kEUbOWEdUp+yWm3AqTbsSZxSz5qo1jRpwdh/Y2\nC8uIGB2K0Zt1aY35jGaqDKQypKwECTNFzIjhmA52MzCWl6YkMITYFDbkvRpi/SwFhtYRfhQsd7H1\nm7WMWuRxtVJmtmZyoZhhruIwU3Jwqz6mCUk7pDXpkY0FDKUrtDghA6kMCStG2spiGAYxM4ZlWI0r\nXZeWpSQzhNg0NnxwKKVGgD8EclrrD6/3eLaipaJ3FIWERNRDjzAKqIYuXlhntlagFoacL8Jc3eHM\nXA+FqsHcTBnfqxKFFXr6s+TSipgV0JWs0xrzGEpbpO0E7bFWDGUsX+m6NNNQSknPKSE2oTUNDqXU\nF4D3ADNa6wMrnn8n8BeACXxOa/1nL/YaWusLwO82L5ISN1Gko0YNoxkYfthod14NqoQ6ZNErUvTq\nnC5CqBU/H2vHMjVhpFicr9LW3jjINzySI2mH9GZdWmIeB9sDDBS9yRy2ESNmxbGVhd2cZZiGhQI5\nnyHEJrXWM44vAZ8Gvrz0hFLKBD4D3AuMA4eVUt+lESJ/+oKv/x2t9cwaj3HbCXWIhuUZhh/6BDqk\nFrrUwzrVoELJ95iruVwsJ5ivOZybS1PzIJOAsXGX9q40jhWxZ1eazpRLNu6zI1VtzjIcMnaKmBkn\nZsaImzFMZTYuXDLM54WG1DWE2HzWNDi01g8ppXa+4OlXAeeaMwmUUl+ncajwT2nMTl4SpdR9wH0A\ngzsGX+rLbGnLS1L6BUXv0CWIAqpBhWrgcbVaph4qLpTS5GsOY4sJqq6mVKyRn9MoBTErojXpkYkF\n9Kddco7PUMoiZdu0xVoxUKSsVHNpysZASWgIsUWsR42jHxhb8XgcePWLfbJSqp3GHecHlVL/thkw\nv0Br/VngswD3HDqob95wN6+loNA6el7XWi/08HUjKIIooOBVKPk+V8rga8WlYivFms1kwcarBVQr\nVQBMy2CgL057sk7KqdCfcsk6PoMph5SVpDWWw8AgaTcCwzFsVPN8BkrJSXAhtoj1CI5r/a3xon/R\na63ngX+xqheWtupAo7151AyLpUN7ERovbDQhrAZV/KjOYr2CryPGKzXydYcz+TSmoXn6okkiYRGP\nK/JzHl09KRbmXUYHLJJOnd50jazjc3uLgcahO5HDNhwSZgLTMIkZDqYysA0boLF7CiQ0hNgi1iM4\nxoGVa0kDwMTNeGFpq/6L9YsgCho/dEgtcPEjj2pQZdGrk/fq5OsW8/UEF/Jp6oHBlckA01KkEorp\nGZeB/gQJx6dnp6IzXSVrB/QnXfqSSVKWRdJKYpuxRi3DaLQMsQzrF2YZIKEhxFaxHsFxGNitlBoG\nrgIfAz6xDuPYMpbOXqA1vg6a9YsAP/KphTX8yMOPfIpelXy9znxds+g55L0MlxZS+JHB1ck6TkzR\n0+UwM+8TsyL2DxvErRrdqRppO6Av6dObSJGwWkmYSWzDIWbGcEwbA2N515QEhhBb21pvx/0a8Bag\nQyk1DnxKa/15pdQfAD+gsZPqC1rr4zfp/bbVUtXSnd5LM4yIRsfaSGtqQY1A+1SDCrWwThBFzNVq\nlAPNRDUJQKFmY5sRrm/S1p4k0lDzobdT05d1aU94pK2A7oRH0rKbS1I2aTuNgbl8LqNxB7gUv4XY\nLpTWW6+OfM+hg/rhxx5c72GsqaU+Uks7pFZe1xrqEDeoNtuD1Fio1wh0xETVpOjbjBcTxKyIMzMp\nsomQ1oTHdCnOjtYq+arD7vYSCSukO14nYRkMpXMAZOwslrKJmY0LlpZbhhiNf3+YNHpNSWgIsfkk\nrJYjWutDq/ncDX9y/EZshxnHcmA0d0stXdfqhR6hjnCD6nIdI+/VWfTqFH2ToucwXm4ExnQpTjbu\ns6erwmw5RjYeEOk6/WmX0VyZrrhPwrIYSLYDkLYzWIaNYz5X9DZQWM3AMJqBISfAhdgejPUewM2k\ntb5fa31fS0tuvYeyJvSK/y3NNMIoJIgC/CigHtZxQ5d65LHo1akEPjOuQy0wmXZjmEpzOZ+kJeHT\nma5TrNnsbKuQdXz2tJXoitfpSQR0xBMMpFqay1JZbMPBMZ3mJUv28vKUQj23LCWhIcS2saVmHFvZ\n0pJihAa9tN02ItQhgQ7xIw8v8poNCUOKvkfRV8zXHRwzwjEjir5Na8IjjBTFmk0u4dPi+LQ6Hgkr\npM2JkbAsck4Sy7BJWWmUUiTMOKYyl+/LMF5wB7gsSwmxvWyp4NiqS1WapdCIlg/zLZ3P8Jvbbf3I\nxwvr1MOAcuATRBGLXhzHjJgsx0k5Id2pGjPVOINZl6JvMZSqkrBCWh2TuGnTlWhcsZK2M1jKWl6a\nsgwLEwPTaJzHkOK3ENvblgqOrXiOY2URXKOfq2k024U0wqKGG9ZwA5+S7zNTC6kFDgt1BwPYma0y\n7caImxG5mE9HvE5/0qXFMYmZDt2JDBpIWekVB/jM5QuWXljLkNAQYnvbUsGxlWg0aJ7bbtsMjzBq\nzjRCn1CHjX5TOqQaeIRaUw18TKVwQ5OUFVILDQqezWDKxTIiWhyPrKOJGSY5xyHnJDGVRazZiNAx\nnEbxu9kqpLEs9fxOthIaQmxvEhwbzC8Ehm4c7FvZZ6rRybaG1lHznEZAvl5HKZh0HZxmQMzVY/Qm\nahR9m5QdkLIibMMgY8eIGSYZO46hTBJW41xHzIphKgNTGShlSA1DCHFNWyo4NkONI9JR4+dmkXup\nVhE2d0l5kU+kI+phjQhNLajiRz4L9RqVIGDSbdQZrlYTlOsWE4UUSSfkwphDOhMjl1ZMz/n0dVks\nFGG02yNuRxR9i46YR8oOMJVPZGoyNgSRTxgFmIZFGIVopVGGwmjWVUxlEKFRaJRu3NanVeNjEihC\nbE9bKjg2Yo1jqUahmyXupRnEUk+poHk9q68Dgshvnsdo3I0R6Ih8vY4bBFQCzUwtRjUwKfsWQWQw\nVYzjWBHjMxHZlgSd2YCpBYPhXoOY5dPZG9ES94lbIUkzZFc2wos0vYk0lmFhmw6msrAMu7FMZdpo\nrRs1Da0xXnAPuASFEAK2WHBsJNeaWWh43m4o3exYq4moh/Xm/Rg+lcCjHPjUgpD5uiLSBgueQzUw\nmavESNgh+aqDbWoqnkkmY+EHEXMlCydmELdd0k5A3A7JOT5JK6A/CaDoTaZQSmEbDrbx3K4px7BB\ns3y4Tze33QIYKwND6hxCbHsSHDfJ0q6nRv+olecsnt8SxAvrzcCoowE3qFILQ+pRSNHzcMOoERD1\nGAqbom9R800WXZt0LKAeGFQ8i3pgEEYK09CYKiLhQH/OJWGFxMyQrB0QN0M64iFx06ItFsdSiriV\nwDGWWoY0DvYByzuo0BpDNReq1HPxoBpbqZ77tRBi29pSwbGeNY7l+7tXBEaow+WfGx1qA/zII9IR\nlcAl0pqC5wGw6PlUA5P5ehI3MMnXbGJmRLFu4/omFa9R25jMW7RmND2ZGlOlOAMtLnMVh8Gci21E\npJuBMZyBehjRm0xjKxONbmy3bQaGYTRahyga92UYysBUJqq51VYptTzTkNAQQqy0pYJjPWscblhr\nzCyad18ES0XuqE6kQ9ywTi0I8ZqnuqsBjYDwGv/iL/oJql7jt2O2HMMxI+ZDg3QsoODatCQ9ZkoO\nXS0+3ek6cxWH3R1lbCOiI1EnbobEm8HRFbeJgL5kGgNwzBiWsprbbE1iVgytNbZqNidURqNBodYo\nZbzoTEMIIWCLBcd6Wmoy6EfecnAEkU85qBFpzaJXR2uYr2vqoUUtMlioO9QDk6pvUqpbxK1Gi/NQ\nK/KuTVvSw/VNsnEfrRXtaR9TaQq1xrJV0gqImxExMyTVXKLK2DaGUnTGk0Rak7JSyzOKlUtTjTs0\nGomw1N32hTUNKYwLIa5FguMmOVe8Qj0MqIYBZV8zV3eIdKNXlAHMVlJUvcZSUL7auFI1iiJqbkAy\n5TTaiASNdiLtWU0mFhBEBo4ZoZQmYYdYhiZuN9qdJ60Ay9DLZzNSlo1txGh1koQ6JGEmm6e+l3ZM\nOaD18qE+jcZsLmEpeMH+KVmeEkK8OAmOm2Ss4lL2LfJenIlygjBSVDyL+bKJYRiUijXicUWtFqCj\ngHo9ACDwQ7TWKKXQWuM4JobSlOoW3Zk6g9kq+brDjnSVkm/Rn3JxjIiEaWAbBh3xFEEUkXUSGMrE\nUo1b+MzmbXxxM0akIwxU87nG/RlLUWDw/IN+6rm0WH5OCCFW2lLBsZ7F8c8/1sGj33qUN33i9cxO\nlZicKDGwI4fWMDtTZmCwBccxKSxUiSVs9o7EGJ+BTMamXo/Y2emRdx16MjXyrs2BniKtMY96aLA7\nVyJjB7TGPNJWowFIznHI2g5eFNESSxPp6P9v7/5j5CjvO46/Pzu7e7fn3zYxBduAKa4brApcG4OU\ntlIlh9AkNKRQ4qT5Iy2FpoFE7R9RUil/IKUVtI1S6iRNYkjiRm1NIxQlpBFFhapBiSjYISZ1sJw6\nRAJDqYFixz/vbm+//WNm7/bWu+fbePf2PPt5SSftPDPzzHNf7c335nlmnmEoGaaAGEqGKRdKSKKo\nJH2dK1NTn5cKpWnJqnFdt9WPY2b5kavE0c/B8ae+9wKlNauIWrB/909Zvf5iTp0c5yc/fpVrf20t\nx4+NUh4rcunlyzhyZJTRasLFF9RYPnKSYhIsKo+zaKjK8uExrrrgCCfGi1w0Mg5ApVhkqDBMqVCg\nnCQsKFaIgJFsqhCpwHAynM4ppQLlQokaQblQImDaS5cmYgJIpzRJSCYfTUxIehIXJw2z/MlV4uin\njZsu5o3/O0l5uMjbbt3E8wdf4xfXrWDDhuWM1wpcclH67MXykdO8aaGolE6zcmSUpBAsK49RLATl\nQo3hJH0WZEkJVlZGKBXSV7QWVSSoZQ/tJQwlwwRBuZAOdqeD30Uiew6jkgxRjQmUPVdSlyghIiav\nQuoD4mZms+XE0SXr15ZZ8MsFjo1WuXDhKG+9fIIXjh9n6dAYR8dKLCmPM14r8AuV04zXCiwoVVkx\nJMZqNZaUyiwuD1GLYDgpMVobZySpIEFBCSWVKCdDAFkXVDpzbY3a5BhFqVBiolalUhphrJZeqRSV\ntLyV1lcBZnYunDi65B2XH+aq5Zdw4OhLXLpwBUfHxrluJfz02EnWLloKlCgXhoAR0vO2qCQjJEqo\nkT4oSASV4sjkXU8lpU9zS9nkgsTkdOe1qKW30WZXFOVCCQolqjGRfsbjC2bWG04cXbJqZCHVWpX1\nS1aRqEglGeFE9ThXr7iUggqcqp6inAyxYmg5x8dPkNSn9QBKhWECJl/NOhETFChMzkxbVMJE1FhQ\nHGG8No5UIEEkFAgFxYbupmm31DppmFkPOHF0ybKhpZyunuT10dOsWbCKclLO7lqqkajIotJiioUi\n1VqVkWKFE9WTjBQrk89TlLOrhUQJpShObptIk+uqMZF2ScUEkSWF5jGK+tiFmVmvOHF0yaUL10x+\nPnzqNRaXFkJpIaWs2+hn48fSq4hCAhFcWHkTkN7lVH9+QvWpP8TkuEX9qqHYkBCcHMysn3KVOObD\ni5zGJsZYWbngjPLFpUXTlmtRm5xYENLxiGLDZ3BXk5nNT7m6FzMivhURdyxduqRvbajPBXU2zV1M\njUnCCcPM5rNcJQ4zM+s9Jw4zM+uIE4eZmXXEicPMzDrixGFmZh1x4jAzs444cZiZWUfmfeKQdJOk\n+yV9U9L1/W6Pmdmg62nikPRlSYcl7Wsqv0HSAUkHJX18pjoi4hsRcTvwAeA9PWyumZnNQq+nHNkJ\nfBb4ar1AUgJ8DngrcAjYLelhIAHuadr/DyLicPb5E9l+ZmbWRz1NHBHxhKTLmoq3AAcj4nkASQ8C\n74qIe4B3NtehdP6Ne4FHIuKZXrbXzMzOrh9jHKuAFxuWD2Vl7XwY2ArcIumD7TaSdIekPZL2vPrq\n691pqZmZnaEfs+O2msEvWpSlKyK2A9vPVmlE7AB2AGzavLFtfWZmdm76ccVxCFjTsLwaeLkbFUu6\nUdKOI0eOdqM6MzNroR+JYzewTtJaSWVgG/BwNyqeD9Oqm5nlXa9vx90FPAmsl3RI0m0RUQXuAh4F\n9gNfi4gfdel4vuIwM+sx1d82lyebNm+M7z31H/1uhpnZeaNSXPr9iNg8m23n/ZPjZmY2v+Qqcbir\nysys93KVODw4bmbWe7lKHGZm1nu5ShzuqjIz671cJQ53VZmZ9V6uEoeZmfVeP+aq6hlJNwI3Aicr\nxaX7m1YvAY7+nMsXAK91sanNxzrX7Wda32qdY9G+zLH4+ZbP11jMtnwQYrFu1keIiNz9ADvOVtbJ\nMrCn1+07l+1nWu9YOBaORfv1sy0f5Fi0+slrV9W3ZlHW6XI3dVr32bafab1jMfM6x6J92SDEYrbl\ngxyLM+RyypFuk7QnZvkoft45FlMciymOxZRBiEVerzi6bUe/GzCPOBZTHIspjsWU3MfCVxxmZtYR\nX3GYmVlHnDjMzKwjThxmZtYRJ45zJOkmSfdL+qak6/vdnn6RdLmkL0l6qN9t6QdJCyT9ffZd+L1+\nt6efBv270Civ54eBThySvizpsKR9TeU3SDog6aCkj89UR0R8IyJuBz4AvKeHze2ZLsXh+Yi4rbct\nnVsdxuV3gIey78Jvz3lje6yTWOTxu9Cow1ic9+eHVgY6cQA7gRsaCyQlwOeA3wKuBN4r6UpJvyLp\nX5p+Vjbs+olsv/PRTroXhzzZySzjAqwGXsw2m5jDNs6Vncw+Fnm3k85jcT6fH86Qq7mqOhURT0i6\nrKuPL8oAAAUfSURBVKl4C3AwIp4HkPQg8K6IuAd4Z3MdkgTcCzwSEc/0tsW90Y045FEncQEOkSaP\nveTwH7IOY/Hc3LZubnUSC0n7Oc/PD63k7gveBauY+s8R0hPCqhm2/zCwFbhF0gd72bA51lEcJK2Q\n9AVgo6Q/63Xj+qhdXL4O3Czp8/R2+on5pGUsBui70Kjd9yKX54eBvuJoQy3K2j4lGRHbge29a07f\ndBqH14Hc/GHMoGVcIuIE8Ptz3Zg+axeLQfkuNGoXi1yeH3zFcaZDwJqG5dXAy31qSz85Dq05LlMc\niykDFQsnjjPtBtZJWiupDGwDHu5zm/rBcWjNcZniWEwZqFgMdOKQtAt4Elgv6ZCk2yKiCtwFPArs\nB74WET/qZzt7zXFozXGZ4lhMcSw8yaGZmXVooK84zMysc04cZmbWEScOMzPriBOHmZl1xInDzMw6\n4sRhZmYdceKw81Y2J9Le7OcVSS81LJf73b65IGl19q6Htzf87sez6b33SvpKm/0K2TMIlzWVf0HS\nRyRtlvTFufgd7Pzj5zgsFyTdDRyPiE81lYv0e17rS8NmIKmYPTh2LnX8DfBYRHy7oey7wF0RsXcW\n+74SEX9Zbw/p1BlXR8Qrkp4Abo2IV86ljZY/vuKw3JF0haR92QytzwBrJB1pWL9N0gPZ5wslfV3S\nHklPS7quRX1FSZ/O1v9Q0h9m5VslPZ7tf0DSVxv2uUbSdyR9X9Ijki7Myr8r6S+yk/JdktZJeiqr\n+5P1dkraJekdDfX9s6S3N7VLwE3Av50lHiVJ92XHeFZSfTLGXaRTY9RtBfY1JIpvA7fOVLcNJicO\ny6srgS9FxEbgpRm22w78VURsJj1JPtBimzuAwxGxBbgGuFPSJdm6XwXuzI73ZknXSRoC/ha4OSI2\nAf8AfLKhvsUR8RsRcR/wGeBTWd3/27DNA2Sz7Upalh330aZ2XZG1a2ymQAB/DBzKjnEt8CeSVkXE\n08AiSeuz7baRJpO6PcCvn6VuG0CeVt3y6icRsXsW220lnXOovrxMUiUiTjVscz1pUqj/d74EWJd9\n/s+I+B8ASXuBy4DTwAbgsazehLQLqO7Bhs/XAvUriX8C/jz7/O/AZyStAN5LOvdR85sFLwJencXv\neD1whaT3Z8uLSZPOS1lbtkm6N2vHnzbsdxi4eBb124Bx4rC8OtHwucb09yUMN3wWsOUs/7UL+FBE\nPD6tUNoKjDYUTZD+TQn4YUS0+2/9RJvySRERkv4ReB/p+6rf12KzU0z/XdoR8EcR8Z0W63YBDwHP\nAk9GxBsN64azY5hN464qy71sYPyNbDyhALy7YfVjpF1NAEi6ukUVjwIfygaPkbReUmWGQz5H+ia8\nLdn2ZUkb2mz7dEN7tjWt+wrwUeB0RBxose8BYO0M7WjX/jdLGgbIZnAdB+5mejcVwC8B+2ZRvw0Y\nJw4bFB8D/hV4nOndRncCb8kGvZ8Dbm+x7xeB/wb2StoHfJ4ZrtYjYhS4Bfi0pGeBH5B2SbXyEeBj\nkp4GVgJHG+p5GfgxaQJpdZyfAS9KOlvy+DvgBeAHWfs/S9p9VreLtOut+ZW3v0k6QG42jW/HNesj\nSQuAk1nX1PuBd0fEzQ3r/gu4KiKOtdn/d4ENEXF3D9r1GPCW+Xgrs/WXxzjM+usa4L6sC+0Npu6k\nehtwP/DX7ZJG5iHSwfpuuwT4qJOGteIrDjMz64jHOMzMrCNOHGZm1hEnDjMz64gTh5mZdcSJw8zM\nOuLEYWZmHfl/sNeSKNVqN1cAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "edisp = cta_irf.edisp.to_energy_dispersion(offset=offset)\n", "edisp.plot_matrix()\n", "print(edisp.data)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# Simulate data\n", "aeff.lo_threshold = lo_threshold\n", "aeff.hi_threshold = hi_threshold\n", "sim = SpectrumSimulation(aeff=aeff, edisp=edisp, source_model=model, livetime=livetime)\n", "sim.simulate_obs(seed=42, obs_id=0)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/deil/Library/Python/3.6/lib/python/site-packages/gammapy-0.7.dev4616-py3.6-macosx-10.12-x86_64.egg/gammapy/data/obs_stats.py:209: RuntimeWarning: divide by zero encountered in true_divide\n", " self.background))\n", "/Users/deil/Library/Python/3.6/lib/python/site-packages/gammapy-0.7.dev4616-py3.6-macosx-10.12-x86_64.egg/gammapy/stats/poisson.py:385: RuntimeWarning: divide by zero encountered in log\n", " m = n_off * log(n_off * temp)\n", "/Users/deil/Library/Python/3.6/lib/python/site-packages/gammapy-0.7.dev4616-py3.6-macosx-10.12-x86_64.egg/gammapy/stats/poisson.py:385: RuntimeWarning: invalid value encountered in multiply\n", " m = n_off * log(n_off * temp)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "*** Observation summary report ***\n", "Observation Id: 0\n", "Livetime: 0.167 h\n", "On events: 165\n", "Off events: 0\n", "Alpha: 1.000\n", "Bkg events in On region: 0.00\n", "Excess: 165.00\n", "Excess / Background: inf\n", "Gamma rate: 0.03 1 / min\n", "Bkg rate: 0.00 1 / min\n", "Sigma: nan\n", "energy range: 0.01 TeV - 501.19 TeV\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlsAAAJgCAYAAAC9eKYfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcXFWZ//HPN3vIShYgZA8EZN/CjhIWkVUQFREV5Keg\nM6AwgorKqKiMOBMRGJcRBERcEEEUIoKAhH1LIEAgLAGCCUnIvq/deX5/3NtJpem9q/reqvq+X696\ndd1zt6c6dW6ePveccxURmJmZmVlpdMo6ADMzM7NK5mTLzMzMrIScbJmZmZmVkJMtMzMzsxJysmVm\nZmZWQk62zMzMzErIyZaZmVUUST+QtFDSvHT5I5JmSVopaZ8inuf9kl4t1vGscjnZshaRdIakyenF\naq6kv0s6rMTnDEk7lvIcZlZ+JM2UtCa9HtW9fpquGw5cBOwaEdulu0wAzo+I3hHxXDvOu8U1KSIe\niYid2/NZmjnfZ9Nznlaqc1jHcLJlzZL0FeAq4L+AbYERwM+Bk7OMy8yq2klp8lT3Oj8tHwksioj5\nBduOBF7q+BDb7SxgcfqzUZK6dEw41lZOtqxJkvoB3wPOi4g/R8SqiNgQEXdFxFcldZd0laQ56esq\nSd3TfT8r6dF6x9v0l6GkX0v6maS/SVoh6SlJO6TrHk53eT79q/UTkgZJmihpqaTFkh6R5O+wmQEg\n6WjgPmD79LrxB0krgc4k15I30u22l3S7pAWS3pL05YJjdJb0TUlvpNelKZKGN3JNGi9pdrrfJZJu\nqxfP1ZKuSd/3k3R9emfgnfRWZ+cmPstI4HDgXOBDkrYtWDde0mxJX09vld6Ylp8oaWp6jXxc0p4F\n+1xS8JlelvSRdvyqrZX8H5U152CgB3BHI+u/BRwE7A3sBRwAXNqK438SuAzYGpgBXA4QER9I1++V\n/tX6R5JbA7OBwSQtbN8E/LwpMwMgIu4HjgPmpNeNT0ZE73T1XhGxQ/oH2l3A88BQ4CjgQkkfSrf7\nCsl16XigL/D/gNWNXJMK/QE4XlJfSJI24DTg9+n6m4AaYEdgH+AY4PNNfJwzgckRcTswHfhUvfXb\nAQNIWu3OlbQvcAPwBWAg8Evgzro/foE3gPcD/Uiuub+VNKSJ81sROdmy5gwEFkZETSPrPwV8LyLm\nR8QCkkr8mVYc/88R8XR6/N+RJG2N2QAMAUamrWuPhB/uaVat/pK24NS9zmnhfvsDgyPiexGxPiLe\nBK4DTk/Xfx64NCJejcTzEbGouYNGxNvAs8ApadGRJEnak2mr1HHAhendgfnATwrO2ZAz2Zyo/Z73\n3krcCHwnItZFxBrgHOCXEfFURNRGxE3AOpI/homIP0XEnIjYmCaKr5P8cWwdwMmWNWcRMKiJPgHb\nA28XLL+dlrXUvIL3q4HejW0I/A9J69c/JL0p6ZJWnMfMKsspEdG/4HVdC/cbSXKbcVOiRtJKXneb\nbjhJK1Bb/J6kVQzgDDYnSyOBrsDcgnP+EtimoYNIOhQYDdxScNw9JBX+MbogItbW+1wX1ftcw0mv\nx5LOLLjFuBTYHRjUxs9preROddacJ4C1JH+t3dbA+jls2fl0RFoGsArYqm5DSdvRDhGxguRW4kWS\ndgMelPRMRDzQnuOaWVWZBbwVEWObWL8DMK0Nx/4T8GNJw4CPkHTDqDvmOmBQE3cJCp0FCJgqqbD8\nTGBq+r5+q/4s4PKIuLz+wdL+X9eR3DJ9IiJqJU1Nz2EdwC1b1qSIWAZ8G/iZpFMkbSWpq6TjJP03\nST+FSyUNljQo3fa36e7PA7tJ2ltSD+C7rTz9u8CYuoW08+eOSq4+y4Ha9GVm1lJPA8vTzuU90w7x\nu0vaP13/K+D7ksYqsaekgem6La5J9aVdKSaRdFh/KyKmp+VzgX+QJGJ9JXWStIOkw+sfI71WnkbS\nMX7vgteXgE81cZfhOuCLkg5M4+4l6QRJfYBeJMnZgvQcZ5O0bFkHcbJlzYqIK0k6jV5KUllnAecD\nfwF+AEwGXgBeJOmz8IN0v9dIRjLeT9I/4NH6x27Gd4Gb0mbv04Cx6bFWkrS4/TwiJrXjo5lZ+bpL\nW86z1dggni1ERC1wEkkC8xawkCTB6pduciVwK0lytBy4HuiZrvsuW16TGvJ74Gg230KscybQDXgZ\nWEJyp6ChDuqnAGuA30TEvLpXGkdn4NhGPtdkkn5bP02PPwP4bLruZeDHJNfNd4E9gMcaid9KQO5f\nbGZmZlY6btkyMzMzKyEnW2ZmZmYl5GTLzMzMrIScbJmZmZmVUK7m2Ro0aFCMGjUq6zDMrB2mTJmy\nMCIGZx1HnvlaZ1YZWnq9y1WyNWrUKCZPnpx1GGbWDpLebn6r6uZrnVllaOn1zrcRzczMzErIyZaZ\nmZlZCTnZMjMzMyuhXPXZMit3GzZsYPbs2axduzbrUEquR48eDBs2jK5du2YdiplZrjnZMiui2bNn\n06dPH0aNGkXyvOzKFBEsWrSI2bNnM3r06KzDMTPLNd9GNCuitWvXMnDgwIpOtAAkMXDgwKpowTMz\nay8nW2ZFVumJVp1q+ZxmZu3l24hmJTKhRMnIxRElOa6ZmZWGW7bMzMzMSsgtW2YlVqyWqFK1lJmZ\nWWmVRcvW+PHjGT9+fNZhmFmBV199lXnz5mUdhplZ7pVFsmVmLTdz5kx22WUXzjnnHHbbbTeOOeYY\n1qxZ0+C2U6dO5aCDDmLPPffkIx/5CEuWLAGSP3C+/vWvc8ABB7DTTjvxyCOPdORHMDOrKE62zCrQ\n66+/znnnncdLL71E//79uf322xvc7swzz+RHP/oRL7zwAnvssQeXXXbZpnU1NTU8/fTTXHXVVVuU\nW+lMmDCBCRMmZB2GmRUoRr0si2Rr0qRJTJo0KeswzMrG6NGj2XvvvQHYb7/9mDlz5nu2WbZsGUuX\nLuXwww8H4KyzzuLhhx/etP7UU09tcv+dd96Z7bbbrvjBV7GJEycyceLErMMwswLFqJfuIG9WYll0\nbO/evfum9507d270NmJLjtG5c2dqamqKFpuZWbUpi5YtN62bFV+/fv3YeuutN/XHuvnmmze1crXE\nvHnzWL58eanCMzOrGGXRslXXfHfxxRdnHIlZy5XD5KM33XQTX/ziF1m9ejVjxozhxhtvbPG+y5Yt\nY/Xq1SWMzsysMpRFsmVmLTdq1CimTZu2abmpP1L23ntvnnzyyfeUF/aRHDRoUIN9tszMrGWcbJmZ\n5YQHApnlTzHqpZMtsypw3nnn8dhjj21RdsEFF3D22WdnFJGZWfVwsmVWZBGBcvZonZ/97GdFP2ZE\nEGXQL62c1A0Ecv9Us/woRr0si9GInmfLykWPHj1YtGhRxSchEcGgQYMYPXp01qFUFM+zZZY/nmfL\nLGeGDRvG7NmzWbBgQdahlFyPHj0YNmxY1mGYmeVeWSRbblq3ctG1a9eqae1xvTQza5myuI3opnWz\n/HG9NDNrmbJItszMzMzKVVncRjQzqwYeCGSWP8Wol27ZMjMzMyuhkiZbkvpLuk3SK5KmSzq4lOcz\ns+ogqbOk5yRNTJdHS3pK0uuS/iipW1rePV2eka4fVXCMb6Tlr0r6UEH5sWnZDEmXFJS3+hytNWHC\nhE0DD8wsH4pRL0t9G/Fq4J6I+Fh6YdqqLQdx07pZ/mRcLy8ApgN90+UfAT+JiFsk/R/wOeAX6c8l\nEbGjpNPT7T4haVfgdGA3YHvgfkk7pcf6GfBBYDbwjKQ7I+Ll1p6jLR+qbsCBR3haJYsIVqyrYcGK\ndSxcsY4FK5Of373r5TYf89ITdmHkwF6MHrQVwwdsRfcunYsWbzHqZcmSLUl9gQ8AnwWIiPXA+lKd\nz8yqg6RhwAnA5cBXlEzXfyRwRrrJTcB3SRKhk9P3ALcBP023Pxm4JSLWAW9JmgEckG43IyLeTM91\nC3CypOmtPUdU+sy2Zk1YX7ORmYtW8cb8lcyYv5If3/daSc/3g79Nb7D8sg/vxh7D+rHH0H507Zxd\nz6lStmyNARYAN0raC5gCXBARq1p7oHKYz2dCweNZLvY11qpAhvXyKuBrQJ90eSCwNCJq0uXZwND0\n/VBgFkBE1Ehalm4/FHiy4JiF+8yqV35gG8+xsDBoSecC5wKMGDGi1R/aLG82bgwWrFzHrMWreXPh\nKr522wvtOt6FR49lcJ/uDOqdvPr17EqPrp3o3qUz3bp0IiJYX7ORdTUbWbG2hsWr1rNo1ToWrFjX\naLL1nTtf2mL52yfuygd2GsyO2/RuV6ytVcpkqwuwL/CliHhK0tXAJcB/Fm7UkguQm9bN8ieLeinp\nRGB+REyRNL6uuIFNo5l1jZU39KdvU9s3d/7NBRHXAtcCjBs3zn+RWYdZu6GWe6bN48I/TgXgCx8Y\nww6De/O129+bHJ196Ci6d+m8Kcnp3qUT62o2smpdDT99cEabY7ji1D0YtvVWbNu3O9v260Gf7l2K\n+gzZz79/zKb3q9fXMH3ucl6cvew9tya/N3Hz8lc+uBMn7DmEHQaXPvEqZbI1G5gdEU+ly7eRJFtb\n8AXIzFrhUODDko4HepD02boK6C+pS9ryNAyYk24/GxgOzJbUBegHLC4or1O4T0PlC9twDrPMvLFg\nJUf9+KEG1/3y4Tcb3e/Gx2a2+Zz/cfRO7LhNb3bYphejBvaiR9fi9Ztqja26dWG/kQPYb+QAPnto\n8kSPucvW8OjrC/lqQevblfe9xpXp7c0rTt2DD++9PVt1K01aVLJkKyLmSZolaeeIeBU4Cmh77zcz\nq3oR8Q3gGwBpy9bFEfEpSX8CPgbcApwF/DXd5c50+Yl0/T8jIiTdCfxe0pUkHeTHAk+TtFKNlTQa\neIekE/0Z6T4PtuYcbfl8HgxkxTDxhTmc//vn3lN+6Qm7MKh3900tXBccNZbhA7aia2dtuj136V+m\nNXrcS0/YhV7duzC4d3eGD9iKYVv3pFf38piuc0i/nnx83HA+Pm4462s28tgbCzn7xmc2rb/kzy9y\nyZ9f5MyDR/LZQ0YxpqC1qxj1stS/pS8Bv0tHIr4JnF3i85lZdfo6cIukHwDPAden5dcDN6cd4BeT\nJE9ExEuSbiX5A7AGOC8iagEknQ/cC3QGboiIl9pyDrMsPDZj4aZE66P7DuO8I3Zg5MBedO60+Zbd\nKfsMbWx3Pn3QyJLHmLVuXTpxxM7bMPOKE1i7oZa7X5zLV259HoDfPPE2v3nibQDu/vL72XX7vk0d\nqsVKmmxFxFRgXCnPYWbVKSImAZPS92+yeTRh4TZrgY83sv/lJCMa65ffDdzdQHmrz9Fa5TAYyPJr\nxdoNXHBL0mr1xcN34OvH7lzUflGVqEfXzpy67zBO3XcY0+cu57irH9m07vhrHuHI923D4HlPMrjL\nunxO/VBMblo3yx/Xy+LzYCBrj1898hYLV65j3xH9+eqHnGi11i5D+jLzihOYu2wN1z38Fjc89hb/\nfGU+MIZeC1/mjKVr2L5/zzYd24/rMTMzK3PL127g+kffAuAbx++yxW1Da50h/Xry7ZN2ZcqlRwMQ\ntTWsGrAzNbVtH8NXFsmWH2Fhlj+ul2b5ceszs1i5roaDxwxk/1EDsg6nIgzs3Z2ZV5zA8Od/xaA3\n/s6IgW16CA5QJsnWxIkTNzWvm1k+uF6a5UNN7cZNUzZ87rDR2QZTgbqsX0HvRQ1PmtpSZZFsmZmZ\nWcMeeGU+7yxdw8iBW3Hk+7bJOhxrQFl0kDczqwYedGBt8afJswH4zEEj6eS+WkVXjHrpli0zM7My\ntWjlOu6f/i6dO4mT9258/izLlpMtM7Oc8KADa607n0+eGnX4ToMZ3Kd7xtFUpmLUy7K4jeimdbP8\ncb0sPs+zZa11+7PJLcSP7jss40gqVzHqpVu2zMzMytCM+SuY9s5y+vbowlG7uGN8npVFsuWmdbP8\ncb00y9Z9L88H4JjdtqNH184ZR2NNKYtky/P5mOWP66VZth58JUm2PN1D/pVFsmVmZmabLVu9gSn/\nWkKXTuKwsYOyDseaURYd5M3MqoEHHVhLPfT6Amo3BgePGUjfHl2zDqeieZ4tMzOzKjTpVd9CLCdO\ntszMcsKDDqwlIoLHZiwE4P07+RZiqXmeLTPLjOtl8XmeLWuJNxas5N3l6xjUuxs7b9sn63AqnufZ\nMjMzqzKPvp60ah264yAkPwuxHJRFsuWmdbP8cb00y8ZjbywC4NAdfAuxXJRFsuX5fMzyx/XSrOPV\n1G7kybpky1M+lI2ySLbMzMwMXnhnGSvW1TBq4FYM7d8z63Cshcqig7yZWTXwoANrzuMzNvfXso7h\nebbMzMyqyKNpsnWYk62y4mTLzCwnPOjAmrJmfS3Pvr0UCQ7eYWDW4VQNz7NlZplxvSw+z7NlTXlm\n5mLW125kj6H96L9Vt6zDqRqeZ8vMzKxKPJg+oscPni4/ZZFsuWndLH9cL8061oOv+HmI5aqkyZak\nmZJelDRV0uS2Hsfz+Zjlj+ulWcd5a+EqZi5aTb+eXdlneP+sw7FW6og+W0dExMIOOI+ZmVlF+mfa\nqnX4ToPp0rksbkpZgbLoIG9mVg086MAa41uI2SmHebYC+IekKZLObWgDSedKmixp8oIFC0ocTj5N\nkJjgh4mamVkDVq6r4am3FtFJScuWlZ9SJ1uHRsS+wHHAeZI+UH+DiLg2IsZFxLjBg/0lMrPq5UEH\n1pBHX1/IhtpgnxFbs3UvT/nQ0XI/z1ZEzEl/zpd0B3AA8HBrj+OmdbP8cb0sPs+zZQ3xLcRs5Xqe\nLUm9JPWpew8cA0wr1fnMzMwqTURsml/riJ2dbJWrUrZsbQvcoaQvUhfg9xFxT1sOVNd857/2zPLD\n9dKs9F6as5z5K9YxpF8PdhnSJ+twrI1K1rIVEW9GxF7pa7eIuLytx/J8Pmb543ppVnp1Uz6M33kb\n5IFUZcuTdZiZmeXUP91fqyJ4ni0zs5zwoAMrtGjlOp6fvZRuXTpx6I4Dsw6nahWjXjrZMrMOJekr\nLdhsVUT8suTBmOXYpFcXEAEHjRnIVt3833U5821EM+toXwV6A32aeF2UWXQZ8jxbVuif6SjEI3f2\nHJRZyv08W8XipnWz/GlHvbw5Ir7X1AbpdDFVx/NsWZ31NRt5+LXkqSpHvm/bjKOpbrmeZ8vMrCER\n8bVibGNWyR5/YyEr1taw07a9GTFwq6zDsXYqi2TLTetm+dOeeinpfZKOktS7XvmxzezXQ9LTkp6X\n9JKky9Ly0ZKekvS6pD9K6paWd0+XZ6TrRxUc6xtp+auSPlQYQ1o2Q9IlBeWtPodZW/39xXkAHL/H\nkIwjsWIoi2TL8/mY5U9b66WkLwN/Bb4ETJN0csHq/2pm93XAkRGxF7A3cKykg4AfAT+JiLHAEuBz\n6fafA5ZExI7AT9LtkLQrcDqwG3As8HNJnSV1Bn5G8jzXXYFPptvS2nOYtdWG2o3c+7KTrUpSFsmW\nmVWUc4D9IuIUYDzwn5IuSNc1OWtjJFami13TVwBHArel5TcBp6TvT06XSdcfpWRmyJOBWyJiXUS8\nBcwgeXbrAcCMdFLm9cAtwMnpPq09h1mbPPXmYpau3sCO2/Rmp209a3wlKIsO8mZWUTrXJUwRMVPS\neOA2SSNpJtkCSFufpgA7krRCvQEsjYiadJPZwND0/VBgVnquGknLgIFp+ZMFhy3cZ1a98gPTfVp7\njoX14j4XOBdgxIgRDX42DwYygL+9OBeA43ffLuNIDIpTL92yZWYdbZ6kvesW0sTrRGAQsEdzO0dE\nbUTsDQwjaYnapaHN0p8NJW9RxPKmzrFlQcS1ETEuIsYNHuyh/NawmtqN/OOl5Bbicb6FWDGcbJlZ\nRzsTmFdYEBE1EXEm8IGWHiQilgKTgIOA/pLqWuqHAXPS97OB4QDp+n7A4sLyevs0Vr6wDedoNQ8G\nsqdnLmbRqvWMHtSL923nW4h5UIx6WRbJ1qRJk9y8bpYzba2XETE7IuYBSNpa0p6S9pW0L7CmqX0l\nDZbUP33fEzgamA48CHws3ewskg74AHemy6Tr/xkRkZafno4kHA2MBZ4GngHGpiMPu5F0or8z3ae1\n52g1DwayulGIx+2+nR88nRPFqJfus2VmmZD0feCzJH2u6pKTus7ujRkC3JT22+oE3BoREyW9DNwi\n6QfAc8D16fbXAzdLmkHS2nQ6QES8JOlW4GWgBjgvImrTuM4H7gU6AzdExEvpsb7emnOYtVbtxuCe\nlzwKsRKVRbJV13znWZXN8qMI9fI0YId01F+LRMQLwD4NlL9J0n+rfvla4OONHOty4PIGyu8G7i7G\nOcxaY8rbS1iwYh0jBmzFbtv3zTocK6KyuI3opnWz/ClCvZwG9C9SOGZl7+50FOJxe/gWYqUpi5Yt\nM6tIPwSekzSNZLJSACLiw9mFZJaNjRuDv0+rm/LBtxArjZMtM8vKTSSzrb8IbMw4llzwQKDq9dys\nJby7fB1D+/dkz2H9sg7HChSjXjrZMrOsLIyIa7IOwiwP7vYoxIrmZMvMsjJF0g9Jpk4ovI34bHYh\nZcuDgapTRPD3Tf21fAsxb4pRL8si2XLTuln+FKFe1o0qPKigrLmpHypa3YADJ1vV5bV3VzJn2VoG\n9+nOPsM9ZiRvilEvyyLZMrPKExFHZB2DWR48NiN5jOZhOw6iUyffQqxEZTH1gx9hYZY/ba2Xkk4s\nxjZmleLxN5Jk65AdBmYciZVKWbRsuWndLH/aUS//R9I7NPwA5zr/BXhyPat4NbUbeerN5FGah+44\nKONorFTKItkys4ryLnBlM9u83hGBmGXthXeWsWJdDaMH9WL7/j2zDsdKxMmWmXWoiBifdQx55cFA\n1efxGb6FmHfFqJdl0WfLzMysEj02YxEAh+zgW4iVrOTJlqTOkp6T5P4XZmZN8GCg6rJ2Qy1T/rUE\ngIPdspVbxaiXHXEb8QJgOtDmR5i7ad0sf1wvi8+DgarL5JlLWF+zkV2H9GVAr25Zh2ONKEa9LGnL\nlqRhwAnAr0p5nkoxwY9osCojaXdJp0k6s+6VdUxmHeWxdMqHQ3d0q1alK3XL1lXA14A+jW0g6Vzg\nXIARI0Y0uI0fYWGWP+2tl5K+A4wHdgXuBo4DHgV+U5QAzXJuU+d4T/lQ8UrWspVOSjg/IqY0tV1E\nXBsR4yJi3ODBgxvcZuLEiZua8cwsH4pQLz8GHAXMi4izgb2A7sWIzSzvlq3ZwIvvLKNLJ3HAqAFZ\nh2MlVsrbiIcCH5Y0E7gFOFLSb0t4PjMrL2siYiNQI6kvMB8Yk3FMZh3i6bcWszFg7+H96dXdszBV\nupL9C0fEN4BvAEgaD1wcEZ8u1fnMrOxMltQfuA6YAqwEns42pGx50EH1eGZmMmv8QWPcXyvvilEv\nW51sSdoaGB4RL7T77GZWtSLi39O3/yfpHqCvrytWLZ56K0m2DhjtW4jVoEW3ESVNktRX0gDgeeBG\nSc09bmOTiJgUEX6wrJltosSnJX07ImYCSyUdkHVcWfI8W9Vh1boapr2zjE6CfUdunXU41oxi1MuW\n9tnqFxHLgVOBGyNiP+Dodp25FSZNmuTmdbOcKUK9/DlwMPDJdHkF8LN2hlXWPBioOjz3r6XUbgx2\nH9qP3u6vlXvFqJctTba6SBoCnAb4SmBmxXBgRJwHrAWIiCWAZ3a0ivf0W8kjejwKsXq0NNm6DLgX\nmBERz0gaA7xeurC25KZ1s/wpQr3cIKkzEACSBgMbixGbWZ7V9dfa3/21qkZLk625EbFnXYfWiHgT\naHGfrfZy07pZ/hShXl4D3AFsI+lykglN/6sYsZnl1bqaWqbOWgrA/m7ZqhotvVn8v8C+LSgzM2uR\niPidpCkkE5sKOCUipmcclllJvTh7GetqNjJ2m95+HmIVaTLZknQwcAgwWNJXClb1BTqXMjAzq1yS\nOgEvRMTuwCtZx5MXHghU+Z6e6Skfyk0x6mVztxG7Ab1JkrI+Ba/lJI/aMDNrtXTm+OclNfxAVLMK\n9bTn16pKTbZsRcRDwEOSfh0Rb3dQTGZWHYYAL0l6GlhVVxgRH84upGy19+Helm+1G4MpM5cATrbK\nSTHqZUv7bHWXdC0wqnCfiDiyzWduBTetm+VPEerlZUUIo6LUDThwslWZps9dzop1NQwf0JMh/Xpm\nHY61UDHqZUuTrT8B/wf8Cqht89nMzFJpy/kmkg4FzgAeangPs/K26RbiKD8Psdq0NNmqiYhflDSS\nJrhp3Sx/ilEvJe1NkmCdBrwF3F6M2MzyaHN/LT+ip9q0dJ6tuyT9u6QhkgbUvUoaWQHPs2WWP22t\nl5J2kvRtSdOBnwKzAEXEERHx02LHaZYHEcEzm0YiumWr2rS0Zeus9OdXC8oCGFPccMysCrwCPAKc\nFBEzACT9R7YhmZXWGwtWsWjVegb17s6ogVtlHY51sBYlWxExutSBmFnV+ChwOvCgpHuAW0gmNa16\nHgxUuepatQ4cPQDJX/dyUox62aJkS9KZDZVHxG/aHYGZVZWIuAO4Q1Iv4BTgP4BtJf0CuCMi/pFp\ngGYl4Pm1qltLbyPuX/C+B8njNZ4FnGyZWZtExCrgd8Dv0j6gHwcuAao22fJgoMpVl2z5eYjlp8Pm\n2YqILxUuS+oH3Nzms7aSm9bN8qeY9TIiFgO/TF9Vy/NsVabZS1bzztI19O3RhZ2365N1ONZKxaiX\nLR2NWN9qYGybz2pmZlYl6vprjRs1gM6d3F+rGrW0z9ZdJKMPIXkA9S7AraUKqj43rZvlj+ulWcvU\n3UI80P21qlZL+2xNKHhfA7wdEbNLEE+D3LRulj/FqJeSRgJjI+J+ST2BLhGxojgRmuXDU+4cX/Va\ndBsxfazGK0AfYGtgfSmDMrPKJ+kc4DY299MaBvylmX2GS3pQ0nRJL0m6IC0fIOk+Sa+nP7dOyyXp\nGkkzJL0gad+CY52Vbv+6pLMKyveT9GK6zzVKx+m35RxmC1as480Fq+jZtTO7D+2XdTiWkRYlW5JO\nA54mGS10GvCUpI+VMjAzq3jnAYcCywEi4nVgm2b2qQEuiohdgIOA8yTtSjKK8YGIGAs8kC4DHEfS\nv3QscC7wC0gSJ+A7wIHAAcB36pKndJtzC/Y7Ni1v1TnaYtKkSR4QVGHq+mvtN3JrunZuazdpy1Ix\n6mVL/+VywvmvAAAgAElEQVS/BewfEWdFxJkkF6f/bNeZzazarYuITa3kkrqwuW9ogyJibkQ8m75f\nAUwHhgInAzelm91EMn8XaflvIvEk0F/SEOBDwH0RsTgilgD3Acem6/pGxBMRESTT2xQeqzXnMPP8\nWga0PNnqFBHzC5YXtWJfM7OGPCTpm0BPSR8E/gTc1dKdJY0C9gGeAraNiLmQJGRsbiEbSvLsxTqz\n07Kmymc3UE4bztFqEyZM2DTwwCqD+2uVv2LUy5Z2kL9H0r3AH9LlTwB3t+vMreBmdbP8KUK9vAT4\nHPAi8AWSa8qvWrKjpN7A7cCFEbG8icefNLQi2lDeZDgt2UfSuSS3GRkxYkSDB/JgoMqybPUGXpm3\nnG6dO7H38P5Zh2NtVIx62WSyJWlHkr/mvirpVOAwkgvLEyQzP5uZtVXd7bfrWrOTpK4kidbvIuLP\nafG7koZExNz0Fl5dS/xsYHjB7sOAOWn5+Hrlk9LyYQ1s35ZzbCEirgWuBRg3blxzCZxVgMlvLyYC\n9h7enx5dO2cdjmWouVuBVwErACLizxHxlYj4D5K/QK9qakdJPSQ9Len5dNTQZW0N0k3rZvlThHr5\nYeA1STdLOiHts9WkdGTg9cD0iLiyYNWdQN2IwrOAvxaUn5mOGDwIWJbeArwXOEbS1mnH+GOAe9N1\nKyQdlJ7rzHrHas05rMq5v5bVaS7ZGhURL9QvjIjJwKhm9l0HHBkRewF7k3Q+PagtQU6cOHFTM56Z\n5UN762VEnA3sSNJX6wzgDUnN3UY8FPgMcKSkqenreOAK4IOSXgc+mC5D8ofhm8AM4Drg39NzLwa+\nDzyTvr6XlgH8G8ntzBnAG8Df0/JWncPsSSdblmruL8keTazr2dSO6Uieleli1/TlpnMz2yQiNkj6\nO8m1oSfJrcXPN7H9ozTcRwrgqAa2D5IpJho61g3ADQ2UTwZ2b6B8UWvPYdVr1boapr2zjM6dxL4j\nt25+B6tozSVbz0g6p36fCkmfA6Y0d3BJndPtdgR+FhFPtTlSK6oJBR2KLw7nwNbxJB0LnA4cQdJf\n6lck8/hVLQ8GqhzP/msJtRuDvYb1o3f3lo5FszwqRr1s7htwIXCHpE+xObkaB3QDPtLcwSOiFthb\nUv/0OLtHxLTCbVoyQiev6hKWYiYrE6Q2Ha89yVNbz2nWTp8FbgG+EBHrMo7FrKjcX8sKNZlsRcS7\nwCGSjmBzs/rfIuKfrTlJRCyVNIlkJuZp9dZ5hI5ZFYqI07OOIW/8cO/KsXl+rYEZR2LtVYx62aK2\nzYh4EHiwNQeWNBjYkCZaPYGjgR+1PkQ3rZvlUVvrpaRHI+IwSSvYsh+nSLpA9S1CeGXJ82xVhrUb\napk6aykA+49yf61yV/J5ttppCHBT2m+rE3BrRHhIoVmVi4jD0p99so7FrBRemL2M9TUbed92fei/\nVbesw7EcKFmylU4ZsU8xjuWmdbP8aW+9lHRzRHymuTKzcvPUm4sA99eyzcri+YaeZ8ssf4pQL3cr\nXEgnNd2vXUGZ5cDjbyTJ1sFj3F/LEmWRbJlZ5ZD0jbS/1p6SlqevFcC7bJ6V3awsrd1Qy5R/LUGC\ng5xsWcqTf5hZh4qIHwI/lPTDiPhG1vHkiQcDlb8pby9hfc1Gdtu+L1v3cn+tSlCMeumWLTPLytOS\n+tUtSOov6ZQsAzJrr8ffWAjAITu4Vcs2c7JlZln5TkQsq1uIiKXAdzKMJ3NFeLi3ZeyxGUl/rUN2\nHJRxJFYsxaiXZXEb0U3rZvlThHrZ0B97ZXFNKhXPs1Xelq/dwAuzl9Klk9h/lEciVopi1Eu3bJlZ\nViZLulLSDpLGSPoJLXjmqllePfPWYjYG7DW8v5+HaFsoi2TLTetm+VOEevklYD3wR+BWYA1wXhFC\nM8tE3S3EQ91fy+opi2TL82yZ5U9762VErIqIS4DxETEuIr4ZEauKF6FZx6rrHH/wDu6vZVsqi2TL\nzCqPpEMkvQy8nC7vJennGYdl1iaLVq7jlXkr6N6lE/uM6J91OJYzvqlsZln5CfAh4E6AiHhe0gey\nDSlbHgxUvp5IH9Gz/6gB9OjaOeNorJg8z5aZlbWImFWvqDaTQMzaadMjetxfyxrgli0zy8osSYcA\nIakb8GVgesYxZaq9D/e2bEQED7+2AIBDPb9WxSlGvSyLZMtN62b5U4R6+UXgamAoMBv4B1U+GtHz\nbJWnNxeuYvaSNQzo1Y09h/ZrfgcrK8Wol2WRbJlZ5ZD0o4j4OnBERHwq63jM2uuhV5NWrfePHUSn\nTso4Gsujsuiz5Xm2zPKnHfXyeEldAT+E2irCQ+ktxMN3GpxxJJZXZZFseZ4ts/xpR728B1gI7Clp\nuaQVhT+LG6VZaa3dUMuT6UjE9491smUNK4tky8wqyqUR0Q/4W0T0jYg+hT+zDs6sNZ56azHrajay\n+9C+DO7TPetwLKfcZ8vMOtoTwL6AW7Hq8WCg8lPXX8u3ECtXMeqlky0z62jdJJ0FHCLp1PorI+LP\nGcRk1iYPvTYfgMN32ibjSCzPnGyZWUf7IvApoD9wUr11AVRtsuV5tsrLzIWreGPBKvr06OJH9FQw\nz7NlZplpa72MiEeBRyVNjojrixpUmfM8W+Xl/unvAnDEztvQtbO7QFeqYtRLfzvMrENJ+hpARFwv\n6eP11v1XNlGZtd59LyfJ1tG7bptxJJZ3ZZFseZ4ts/xpR708veB9/bm2jm17RGYdZ8mq9TwzczFd\nOonxO7tzvDWtLJItz7Nllj/tqJdq5H1Dy2a59OCr89kYcNCYgfTt0TXrcCznyiLZMrOKEo28b2jZ\nLJc23ULcxaMQrXll0UHezCrKXulM8QJ6FswaL6BHdmFlz4OBysPaDbWbHtHj/lqVL9fzbEkaDvwG\n2A7YCFwbEVeX6nxmVh4ionPWMZi1xxNvLGL1+lp2GdKXYVtvlXU4VgZK2bJVA1wUEc9K6gNMkXRf\nRLxcwnOamZUtz7NVHu56YQ4Ax+62XcaRWEfI9TxbETEXmJu+XyFpOjAUaHWyVQlN6xO0Zb/fi6Px\nrin1t21qv4a2rStrybbNnb+pONu6vVWGSqiXeeN5tvJvXU0t972U9Nc6ca8hGUdjHaEY9bJD+mxJ\nGgXsAzzVwLpzgXMBRowY0RHhFEVhgtFUwtPSY7TlvK1JhJrbti2xtOR8LT1/S4/X0fvnUSV+JrNy\n8fBrC1mxroZdhvRlh8G9sw7HykTJRyNK6g3cDlwYEe958GxEXBsR4yJi3ODBDc9V4nm2zPLH9dKq\n0cT0FuKJe7pVy1qupMmWpK4kidbv2vNwWc+zZZY/rpdWbdZuqOX+dMqHk/bcPuNorJyULNmSJOB6\nYHpEXFmq85iZmXWEB1+Zz6r1tew5rB8jBnoUorVcKftsHQp8BnhR0tS07JsRcXcJz2lmVrY86CDf\n/jrVtxCrUa7n2YqIR/GjN8zMrAIsXrWeB15JbiGevPfQjKOxcuPH9ZhZ2ZB0g6T5kqYVlA2QdJ+k\n19OfW6flknSNpBmSXpC0b8E+Z6Xbvy7prILy/SS9mO5zTdodok3naAsPOsivvzz3Dhtqg/E7D2bb\nvlX9oIOqU4x6WRbJ1qRJk9y8bpYzGdXLXwPH1iu7BHggIsYCD6TLAMcBY9PXucAvIEmcgO8ABwIH\nAN+pS57Sbc4t2O/YtpyjrTzoIL9umzIbgI/vNzzjSKyjFaNelkWyZWYGEBEPA4vrFZ8M3JS+vwk4\npaD8N5F4EugvaQjwIeC+iFgcEUuA+4Bj03V9I+KJiAiSx42d0sZzWAV5ac4yXp67nP5bdeXoXf3g\naWu9ski23LRulj85qpfbpk+sqHtyRd3/hkOBWQXbzU7Lmiqf3UB5W87xHpLOlTRZ0uQFCxa06gNa\ntv40OflanLzX9nTv4kd7WuuVRbLlpnWz/CmDetnQAJ1oQ3lbzvHewhZM4Gz5s3ZDLX+Z+g4AHx/n\nW4jWNmWRbJmZNeHdult36c/5aflsoPB/x2HAnGbKhzVQ3pZzWIW4c+oclq7ewJ7D+rHb9n2zDsfK\nlJMtMyt3dwJ1IwrPAv5aUH5mOmLwIGBZegvwXuAYSVunHeOPAe5N162QdFA6CvHMesdqzTnaxIOB\n8iUi+PXjMwE46+BRqEjPkLXyUox62SEPojYzKwZJfwDGA4MkzSYZVXgFcKukzwH/Aj6ebn43cDww\nA1gNnA0QEYslfR94Jt3uexFR1+n+30hGPPYE/p6+aO05rDI8+68lvDx3OQN6deMET2Rq7eBky8zK\nRkR8spFVRzWwbQDnNXKcG4AbGiifDOzeQPmi1p6jLeoGHFx88cXFOqS1w42PzQTgkwcMp0dXd4yv\nVsWol2WRbLlZ3Sx/XC+Lr27AgZOt7M1cuIq7X5xLl07i0weNzDocy1Ax6qX7bJmZmdVz7SNvsjHg\nI/sMZUi/nlmHY2WuLJKtHM3nY2Yp10urVPOXr+W2dG6tLxy+Q8bRWCUoi2SrDObzMas6rpdWqa59\n+E3W127k2N22Y8dtemcdjlWAski2zMzMOsKcpWv4zZNvA3D+kTtmHI1VirLoIG9mVg086CB7V93/\nGutrNnLinkPYfWi/rMOxHChGvXTLlpmZGTBj/gpuTftqXXTMzhlHY5XEyZaZWU540EG2/vueVwH4\n1IEjGD2oV8bRWF4Uo16WxW1EN62b5Y/rZfF5nq3sPPTaAv7x8rv07NqZLx81NutwLEc8z5aZmVk7\nrd1Qy1k3PA3AhUePZdu+PTKOyCpNWSRbblo3yx/XS6sUP5/0BgBjt+nN/ztsdMbRWCUqi2TL8/mY\n5Y/rpVWCZ/+1hJ89OAOAH5yyO107l8V/i1Zm/K0yM7OqtGLtBi645TlqNwbnfmAMB44ZmHVIVqHK\nooO8mVk18KCDjrNxY/C1215g1uI1AFzsqR6sEZ5ny8zMrA2ufuB1/j5tHgD/vOhwunXxf4dWOv52\nmZnlhAcddIw/Pzubqx94HYAbz96fMYP9/ENrnOfZMrPMuF4Wn+fZKr2/Tn2Hr9z6PACXnrALR+y8\nTcYRWd4Vo16WRbJlZmbWHhHBzx6cwYR/vAYk82l9/v1jMo7KqkXJki1JNwAnAvMjYvf2HKuu+c5/\n7Znlh+ullYt1NbV84/YX+fNz7wDwreN34fPv93xa1nFK2Wfr18CxxTiQ5/Mxyx/XSysHU95ezM6X\n3rMp0bruzHGc84ExSMo4MqsmJWvZioiHJY0q1fE7yoSCCnlxxHvKmtunPdu09xiF69t6vrr96j57\nqbTk99zQNo2tK8bv18zK18p1NfzPPa9w0xNvbyqb+KXD2H1ovwyjsmqVeZ8tSecC5wKMGDEi42ia\n5v/AE00lPY1tV/9319L9mtumocSqvf9OrYm1qbjMWsuDDorjwVfnc+kd03hnaTKH1r+P34EvHzWW\nHl07ZxyZlaOKmGcrIq6NiHERMW7w4MFZh9NmLf3P1f8Jm5mVxuJV67nwluc4+8ZnNiVaf/vyYXzt\n2Pc50bJMZd6yZWZmCQ86aJuI4M7n53DBLVO3KJ9x+XF08bMOrZ2KUS/LItly07pZ/rheFp/n2Wq9\nd5au4dAr/rlF2aSLxzNqUK+MIrJKk+t5tiT9ARgPDJI0G/hORFxfqvOZmVl1+evUd7Zozbri1D34\nxP7DPdLQcqeUoxE/WaxjuWndLH9cLy0rq9fX8N07X+LWybMBOGbXbfn+Kbuzbd8eGUdm1rCyuI3o\npnWz/HG9tCy89u4KjvnJw5uWL//I7pxxwAi3ZlmulUWyZWZm1S0i+OMzs7jkzy9uKrvnwvfzvu36\nZhiVWcs42TIzywkPOmjYirUb+NYd07jz+TkAfGy/YXzv5N3Yqpv/C7PSK0a99DfVzMxya9o7yzjx\nfx/dtHzlaXtx6r7DMozIrPWcbJmZ5YQHHWwWEdz0+Ey+e9fLm8oeuOhwdhjcO8OorBp5ni0zy4zr\nZfF50EFi3rK1HPTDBzYtf+agkXzrhF08C7xlItfzbJmZmbVGRPDXqXO48I+b5876+af25fg9hmQY\nlVn7lUWy5aZ1s/xxvbRimrV4NZfd9RL3T58PwBE7D+aKj+7pubOsIpRFsuWmdbP8cb20Yli+dgPX\nP/IWVz/w+qYyzwRvlaYski0zM6ssy9Zs4KbHZ3Llfa9tUf70N49iG7dmWYVxsmVmlhPVMOjgjQUr\nuenxmfzmibe3KL/1CwdzwOgBGUVl1jjPs2VmZrm3Zn0t909/ly/94bn3rPvDOQdx8A4DM4jKrOM4\n2TIzy4lKGnSwen0Nj89YxOd/M7nB9fde+AF23q5PB0dl1nqeZ8vMMuN6WXzlPOhg7YZapr2zjMlv\nL+GKv7/S4DbfPnFXPrLPULbu1a2DozNrO8+zZWZmHW5D7UbeWriKl+cs32JOrIZc9MGd+PDe2zNy\nYK8Ois4sf8oi2aqkpnWzSuF6uZmkY4Grgc7AryLiioxDKooNtRuZtXg1Mxet4o35q7j87ukt2u+a\nT+7DYTsOYoBbsMyAMkm2yrlp3axSuV4mJHUGfgZ8EJgNPCPpzoh4uek9s7NxY7ByfQ0r1tawbPUG\n3l2xlneXreXd5ev4yf2vNX+AAl/90M7sM7w/ew7vT+/uZfFfilmHc80wM2ufA4AZEfEmgKRbgJOB\nVidbS4Ydwvwu2zHqkr8VOcTi+P4pu7PLdn3Yebs+9OnRNetwzMqGky0zs/YZCswqWJ4NHFh/I0nn\nAucCjBgxosEDHfyhj3LPS/NKEGLLfe/k3diubw9GDerFiAFb+eHPVvU8z5aZWfYaeqZMvKcg4lrg\nWoBx48a9Zz3ABUeP5cyDR4JACCk5eKdOQkDy9BrRSSBtLuuUPtZGBft10ub96/br2a0zvbt3oXf3\nLnTu5EfhmHUUJ1tmZu0zGxhesDwMmNOWA+0ypG9RAjKzfCmLZMvz+Zjlj+vlJs8AYyWNBt4BTgfO\nyDYkM8uTski2zMzyKiJqJJ0P3Esy9cMNEfFSxmGZWY6URbLl+XzM8sf1crOIuBu4O+s4zCyfOmUd\nQEtMnDhx05w+ZpYPrpdmZi1TFsmWmZmZWbkqabIl6VhJr0qaIemSUp7LzMzMLI9KlmwVPMLiOGBX\n4JOSdi3V+czMzMzyqJQd5Iv2CAuAWQ89xAR5Er6sNPa7b215S9cX6zxtPV979+2I72rm9eHww7M9\nv5lZmVBEgxMZt//A0seAYyPi8+nyZ4ADI+L8ettteoQFsDPwKtAPWFbvkIVlzb0vLBsELGxl+A2d\nvyXr8x53S2JsqKw17/Mad+FyseJuLIbWxNhQWbnGXbfcPyIGtyGmqiFpAfB2utjYv3dDy3n5LjQX\nZ/3lco+78L3jbt36So97ZIuudxFRkhfwceBXBcufAf63hfte21RZc+/rlU1uQ+zvOX9L1uc97pbE\n2Ja4632GXMbdxO+5zXG39btSqXG39N/Sr9b9Dltwrejw70K1xd3QZ3Dcjrs1cZSyg3x7HmFxVzNl\nzb1vaP/WaG7/xtbnPe7Gtmlv3C09d2tiask2rYm7cLlYcbdk/2qKuyXHtfdq7nfY3LUiL9e7So67\n8L3jbt36aoi7WaW8jdgFeA04iuQRFs8AZ0QHz6wsaXJEjOvIcxaD4+5YjtvKXbl+Fxx3x3Lc2ShZ\nB/nIzyMsrs3gnMXguDuW47ZyV67fBcfdsRx3BkrWsmVmZmZmnkHezMzMrKScbJmZmZmVkJMtMzMz\nsxJysmVmZmZWQlWbbEk6RdJ1kv4q6Zis42kpSWMkXS/ptqxjaY6kXpJuSn/Pn8o6npYqp99xoXL9\nTlvplet3o5zqoq93HavsvtNtmZE16xdwAzAfmFav/FiSx/3MAC5p4bG2Bq4vw7hvy/vvnuSpASel\n7/9Ybt+ZrH7HRYi7w77TfuXzO9DEsXy9K9Fn8PUus7jL4nqXeQBt/Af5ALBv4T8IyVxebwBjgG7A\n88CuwB7AxHqvbQr2+zGwbxnGndXFpzWf4RvA3uk2vy+X70zWv+MixN1h32m/8vUd8PUu08/g6102\ncZfF9a5kk5qWUkQ8LGlUveIDgBkR8SaApFuAkyPih8CJ9Y8hScAVwN8j4tnSRpwoRtxZa81nIHlk\n0zBgKhnfsm5l3C93bHSNa03ckqbTwd9pKz1f77Lj613HquTrXSX12RoKzCpYnp2WNeZLwNHAxyR9\nsZSBNaNVcUsaKOn/gH0kfaPUwbVQY5/hz8BHJf2CfD4zr8G4c/o7LtTY7zsv32krPV/vsuPrXceq\niOtdWbZsNUINlDU6PX5EXANcU7pwWqy1cS8C8vbFavAzRMQq4OyODqYVGos7j7/jQo3FnZfvtJWe\nr3fZ8fWuY1XE9a6SWrZmA8MLlocBczKKpTXKNe5C5foZHLeVq3L9DpRr3IXK9TM47gxVUrL1DDBW\n0mhJ3YDTgTszjqklyjXuQuX6GRy3laty/Q6Ua9yFyvUzOO4sZd1Dvy0v4A/AXGADSdb7ubT8eOA1\nkpEL38o6zkqJuxI+g+P2q1xf5fodKNe4K+EzOO78vZR+EDMzMzMrgUq6jWhmZmaWO062zMzMzErI\nyZaZmZlZCTnZMjMzMyshJ1tmZmZmJeRky8zMzKyEnGxVGUm1kqYWvC7JOiYASTMlvShpnKQ70thm\nSFpWEOshjez7eUk31yvbVtJ8SV0l/VHSYkmndMynMbM88PXO8sLzbFUZSSsjoneRj9klImraeYyZ\nwLiIWFhQNh64OCJObGbfrYHXgWERsTYtOx/YIyK+kC7/FrgtIv7SnjjNrHz4eufrXV64ZcuATX9p\nXSbp2fQvrvel5b0k3SDpGUnPSTo5Lf+spD9Jugv4h6ROkn4u6SVJEyXdLeljko6SdEfBeT4o6c/t\niHN/SQ9JmiLp75K2jYglwOPACQWbnk4yG7GZ2RZ8vbOO5mSr+vSs16z+iYJ1CyNiX+AXwMVp2beA\nf0bE/sARwP9I6pWuOxg4KyKOBE4FRgF7AJ9P1wH8E9hF0uB0+WzgxrYELqk7cDXw0YjYD/gt8P10\n9R9ILjhIGp7G8nBbzmNmFcPXO8uFLlkHYB1uTUTs3ci6ur/AppBcTACOAT4sqe5i1AMYkb6/LyIW\np+8PA/4UERuBeZIeBIiISPsXfFrSjSQXpTPbGPsuwG7A/ZIAOpM8PwuSB5NeI6k38Ang1jQWM6te\nvt5ZLjjZskLr0p+1bP5uiOQvq1cLN5R0ILCqsKiJ494I3AWsJblAtbW/g4AXIuL99VdExCpJ9wMn\nk/zF929tPIeZVQdf76zD+DaiNede4EtK/7SStE8j2z0KfDTty7AtML5uRUTMAeYAlwK/bkcsLwND\nJR2QxtJN0m4F6/8AfBXoHxHPtOM8ZladfL2zknCyVX3q92G4opntvw90BV6QNI3NfQbqu52kiXsa\n8EvgKWBZwfrfAbMi4uW2Bh4R64CPAVdKeh54DjiwYJN7SJr8b2nrOcysovh6Z7ngqR+saCT1joiV\nkgYCTwOHRsS8dN1Pgeci4vpG9p1JvaHQRY7NQ6HNrGh8vbPWcMuWFdNESVOBR4DvF1x4pgB7koym\nacwC4AFJ44odlKQ/AoeS9KEwMyuGZq93kr4p6VcN7JvJ9U7SpyT9o9jntOa5ZcvMzDKTtvJsS9JR\nvc6vI+L8bCJqGUmTgIOADUCQTDT6J+An6S1As03csmVmZlk7KSJ6F7yKnmhJKsXo+/Mjog8wBLiI\nZGTg3XUd7DuaEv5/PYf8j2JmZrmUztz+qKQJkpZIekvScQXr+0m6XtJcSe9I+oGkzgX7PibpJ5IW\nA9+V1FnSjyUtTI91vqSQ1EXSx9NbgIXnv0hSs/2eImJVREwCPkwyt9YJ6f7fTftPIamHpN9KWiRp\nqZJZ6rdN102S9ENJTyt5PuJfJQ0oiOMgSY+n+z2v5NE+FOx7uaTHgNXAmPSzvylpRfo5P1X4+yzY\n95A0jmXpz0PqHff76e9whaR/SBrU8n89K+Rky8zM8uxA4FVgEPDfwPUFLUc3ATXAjsA+JJOSfr7e\nvm8C2wCXA+cAxwF7A/sChQ9rvhMYLWmXgrJPA1s89LkpEfEvYDLwnrmxgLOAfsBwYCDwRWBNwfoz\ngf8HbJ9+pmsAJA0F/gb8ABhAMtv97do8Sz3AZ4BzgT4k/cGuAY5LW90OAabWDyZN5v6WbjsQuBL4\nW9rhv84ZJLPgbwN0Y/NM+9ZKTrbMzCxrf0lbbepe5xSsezsirouIWpLkagiwbdoqdBxwYdqyNB/4\nCeljbFJzIuJ/I6ImItYApwFXR8Ts9BmDm6aCSPtZ/ZEkwSKd02oUMLGVn2UOSVJU3waSpGbHiKiN\niCkRsbxg/c0RMS0iVgH/CZyWttJ9Grg7Iu6OiI0RcR9JQnd8wb6/joiX0glUa4CNwO6SekbE3Ih4\nqYF4TgBej4ib09/PH4BXgJMKtrkxIl5Lf3e3kiSp1gZOtszMLGunRET/gtd1Bevm1b2JiNXp297A\nSJI5sebWJWkkc15tU7DvrHrn2b5eWf31NwFnpC1nnyF5DE5rO7sPBRY3UH4zyaSpt0iaI+m/JXVt\nJJa3ST7bIJLP+fHCZJTkcUFDGto3TdY+QdJyNlfS35Q+aLue7dPzFHo7jb/OvIL3q0l+79YGTrbM\nzKwczSJ55M6ggiStb0QUzrJef7j9XGBYwfLwwpUR8SSwnuQ24Bm04hYibHoo9H4k00FsISI2RMRl\nEbErya29E9nyuYmFsYwgaQlbSPI5b66XjPaKiMIJWrf4nBFxb0R8kCQhewUoTF7rzCFJ5AqNAN5p\nwUe1VnKyZWZmZSci5gL/AH4sqa+SR+fsIOnwJna7FbhA0lBJ/YGvN7DNb4CfAjUR8WgD699D0lbp\nef9KMsHp3Q1sc4SkPdJbg8tJkqnC6S4+LWlXSVsB3yOZlLSWZH7CkyR9KO3g30PSeEnD6p8jPc+2\nkj4sqRdJMrqy3nnq3A3sJOmMdIDAJ4Bdaf1tU2sBJ1tmZpa1uyStLHjd0cL9ziTpuP0ysAS4jS1v\nr20NNH0AACAASURBVNV3HUmC9gLJ42/uJunjVJiM3AzsTstatX4qaQXwLnAVyWN8jo2IjQ1su10a\n33JgOvAQW070fDPJsxTnAT2ALwNExCySB05/k6Tz+yySZyI29v93J5JpKOaQ3M48HPj3+htFxCKS\n1rWLgEXA14ATSzWrfbXzpKZmZlaV0mkk/i8iRhaU9QTmA/tGxOsdFMck4LcR0dBs81YB3LJlZmZV\nQVJPScent82GAt8B6rei/RvwTEclWlYdSjGjrpmZWR4JuIxkioc1JPNMfXvTyuTRQWLL+bfM2s23\nEc3MzMxKyLcRzczMzErItxHNzDqIpJOAk/r06X3O2J12zDocM2unZ6dMXRgRg5vbzrcRzcw62H7j\n9onHnpqUdRhm1k49u/SfEhHjmtvOtxHNzMzMSsjJlpmZmVkJOdkyMzMzKyEnW2ZmHUTSSZKuXbp0\nWdahmFkHcrJlZtZBIuKuiDi3f/9+WYdiZh3IyZaZmZlZCTnZMjMzMyshJ1tmZmZmJeRky8zMzCpe\nlpO4O9kyM+sgHo1olp2NbMzs3E62zMw6iEcjmmXj/7N331Fyl3e+599P5VzV3dW51ZJaCZGDDDiH\nAQaMZcAJg+0J9pidc67n7v3j7t6ZM7M75+yZ3bn3nrtp7nh3lrmDPR4HbBwwDLYxDhibZCRAKOfY\nQZ2rqivX7/fsH1UtBAgTpO7q6v68fHSk+qn7V19oWf3h+zy/7+NYh5pba9r7K2yJiIjIsuVYB8d1\nKLuVptWgsCUiIiLLkmtdHOtSdisUa6Wm1aGwJSIiIsuOa11q1qHi1INWttK8vZIKWyIiIrKs1Dta\n9aBVcIpkq7Mczk02rR6FLREREVk2rLX1oOVWKTllspUMR7JTHMoGm1aTwpaIiIgsC68IWrUSmcos\nR3OTHMr52TbS1rS6FLZERESk5b26ozUftA5m/fx2uJ3MnIaaiogsexpqKrJw5oNWsVZkpjzN3sxp\nnp8K8dN9bRw/WWD7b081rTaFLRGRRaKhpiILo2YdqrZW36NVzXJ8boojuRAvjSWpVh1Gh7Ok05Gm\n1edr2juLiIiInKf56fAlp0yumuX43CQHs162j7SRm3M4dniaRDLI1ZcnmlajOlsiIiLSkqy1Z4JW\ntpLlxNwEh7Pw3Eh9j9bRQ1PE4wGuuzLOtb3TTatTnS0RERFpSY51KDY2wx/Onmb3rI9fH0szPTHH\nwQOTRKMBLr04yWWd01yS8jetToUtERERaTk161B1q8xVcwznJzmaM+wYS1Ep1xgdydHbG+cdV8TY\n0jPBZW0heiPdTatVYUtERERaSs06VJ0Kc7UCx3Jj7Jwx/PJoNzNTeQ7smyAWC3LJxUku75riklSQ\nvkgPUb82yIuIiIi8ofmOVr5WJFOZZeeM4RdH0kyO5zi4b5JozM/113Vyff8UV3f46Y/0EvNHCXoD\nTatZYUtERERagmMdHLdWP1i6muHE3DTPDneSy5Q4fmSGdGeEd72jjev7Jri8PUxvpIeoL0LIG8Rr\nvE2rW2FLREREljzXujhufUN8rppjrDDDz4dj5IqGkVMZ+gYSvPOKANd1j3NxKklnuJuoL0LAG8Br\nPBhM02pX2BIREZElzbUuNbdG3imSrWQ4nhvn+0eS7Dtu2bv7RH2P1qYol3VMclEyvqSCFmjOlojI\neTHG3G6M+UdjzA+NMTc1ux6R5WY+aBWc0pmg9d3DSXYfcTiwd4Jg0Mf113Xyrv5Jru6I0B3pXVJB\nCxS2RERewxhznzFm3Biz61XXbzbG7DfGHDLG/DmAtfZBa+0XgT8C7mxCuSLLlmtdatZpDC3NMFqY\n4kcnYhwZ9TA2nGVgMMmNH+jkhtWn2ZKO0RPpI9YIWj7jXRJBCxS2RETO5avAzWdfMMZ4gS8DtwAX\nA3cZYy4+60P+qvH7InIBnOlo1YrMVmY5MTfOD4+F2XXCz+EDE8TiQSLRAJvTOYbiITpCaYLeAH6P\nH+8SizdLqxoRkSXAWvsE8OqzPa4FDllrj1hrK8D9wG2m7j8BP7bWPr/YtYosVxaLg0vVrVKsFdg+\n5eX4TISjh6YIh/2EIn4uXgPv7PKwNjFA3B8j4gvj8/gwZml0tOZpg7yIyJvTD5w86/Up4Drgz4Ab\ngKQxZr219h/O9cnGmHuAewBWDa5a4FJFWptrXSpulcniJIdzY3xrfzsv7S6x7/kdpHrb+MiH0tw6\nOMdQopeu0BoC3gAWu2SWDV9NYUtE5M0519/i1lr7d8DfvdEnW2vvBe4FuGbLVfYC1ybS8iz2zNJh\nrppnrDDC81MFvr+rk2OHpzj+wmGIJ7n6ulWsTU6SCATxewJ4jAdr7bn/H7pEKGyJiLw5p4CzW1ID\nwMhbuYExZiuwdWjd2gtZl8iyYDC41qXsVshWM5wqzPHUSCfTk3mOHzhNuK+HW24Z5MY1o1zZnqIz\n1EU8EMfn8S3prhZoz5aIyJv1HLDBGLPWGBMAPg089FZuYK192Fp7TyqVXJACRVpVzToUnRIzlQzH\ncsd5dvw039jTy54DRfbtGadzoJ0P3zLITWtOc2V7iu5wD1F/FF9jKvxSDlqgsCUi8hrGmG8BTwOb\njDGnjDFfsNbWgC8BjwJ7ge9Ya3e/xftuNcbcOzubufBFi7QoxzrU3Bpz1TwTxdPsnMnzyOFejhzJ\nsPOFYeLxIB94bzc3rx3jqo4OeiJ99TlaHj8e0xoxxlirrQMiIovpmi1X2SeffbzZZYg0Xc06VJ0K\n2eoco4Ux9mXyPHGyk2NjLnt3nibVFuI916Z4T/8klzc6WmFfiIDH39SzDueFfant1totb/Rx2rMl\nIiIii86xDqVaiUwly4n8GDumXX52qJfJ8Tn27honGg1wyWWdXNk1waZklHSo88wcLU+LLcwpbImI\nLBJtkBd5mWtdSk6Z2coMx+fKPH60n7HhWXbvPE08EeS97+vnQ2tO847GZPiIL9xSS4dna72KRURa\nlDbIi9TNH8OTqcxyYi7Lz451c+rELLt3jJBIBrnlpn7u2DjMe7rTrIoNEvdHCXoDLRm0QJ0tERER\nWUTzA0vnqnmOzU3yi+E2jg+XObBvgq7eBB98Xxc3DJ7m4lRn4wieIN4ldM7h26GwJSKySLSMKCuZ\ntZaarVF0SuQqOU7kT/OdA2mOnKyw4/kRYrEg7313D+8dGOfiVJJ0qIugx4/XeFq2ozWvtasXEWkh\nWkaUlc5aS8WpMluZ4cGjYba9MMu2n+7C7/fy0Zu7uPuiSd7XM8Dq+BpSgQQhX2hJPHV4vhS2RERE\nZMFV3Sq56hzjxTG+ut/HL54tcuTpPXja2vjIh1fx4dXTbEj20RFsJ+QN4jGell46PJuWEUVERGTB\nnDmCp5JjtDDG1/YHefqFIod2nCK1cQ0fvqmHW9eeZmOyj1QgSaCFN8K/HoUtEZFFoj1bshLVrEOx\nVmK2MsNDJzw8t6vMvu3H6FzTxW03tXPL4AzrEwO0BduI+iKYxv+Wk+UVHUVEljDt2ZKVxLEOxVqR\n2XKGE3Mn+ep+l1/t8HBw/yQDG3tZtzHN+/qmWBProL2xdGiMwZjlFbRAnS0RERFZANZaqm6NbHWW\nh08afrWtzN6nD5IY7GZwqJ27r5nius6NpIJJ/Ma3LEPWPIUtERERuWDm52jlqnOMFEb4+oEAT73o\nsvelUXo3r2L1UDs3X5rjyvYYcX+0vhF+GQctUNgSERGRC8ixDkWnxEx5ikdOwra9DnteGqV3sI2L\nLunmtovHeWdXG13hboKNpw6Xu+X/TygiskQYY7YaY+6dnc00uxSRBVF1q8zVCowXT/PNw1Ueedry\n4s93E09FWD3Uzq0XTXJ9Z5LOUCcxX2RFBC1Q2BIRWTTaIC/LVX1YaYW5ap7xwhjfO1bmqT0+9u2Z\noPfiQS69oodbL5/jnV0ROkKdRH2RZTGs9M3SMqKIiIicFweXkltmpjzNT4aLPL4zzHO/PkK6N8mG\nizr5+GWTvKu7na5QNxF/BN8KClqgsCUiIiLnoebWKNSKjJfGefhEjp/sCPPcLw8Q60rSP5jkfRty\nXJOO0hZorx+/swIX1RS2RERE5G2pNTbDT5Yn+deTOX78YpjnfrGPcEeKofUd3HZtlZv6A/REekn4\nY/g9/maX3BQKWyIiIvK2FGtFJkoTPHA0x0NPedjzi+cIr17FZVf2cPe1Gd7TnaYzvPL2aL3ayuvl\niYg0iZ5GlOVifjr86eIYD53I8PieMHt3T9CxeYiNm9N88NIKl7eHiPsTBD0BPCs8bqzsf3oRkUWk\npxFlOXCsQ9mpMFGa5FtHSnz/KR/P/GQXoZCPvoEEt2xxuHkA+qN9pAIJAp7Ash9a+ka0jCgiIiJv\nyvx0+JnyLCfmxvnFC2F2/vYI7YNdbLiok09ek+W9PQm6wj3EfBF8HsUMUGdLRERE3oSXg9YMx+dG\n+eGxGC8+d5K2njaGNnRw65V53tUdIx3qIuoLK2idRf8mRERE5A3VrEO2kmMkP853D8d4/JkcwXCA\nofXtfOraOT7Qk6Ir3EXMH13Rm+HPRZ0tEREReV3WWspOmanSNIdzJ/na/giPP5Pj6G92csXVfdx1\nff5M0IoqaJ2TOlsiIiLyuhzrkKvmGS+N86/HQzz5XI6jv9gGq4f4zJnxDl1EV9BZh2+VwpaIiIi8\nxvwerenyNMdyo9x/MM4zL8xx5Jm9BDZu5K6PD3B9F3SE0oS9IQWt30FhS0RERF6j3tGa43Rxgu8f\nifKb32Y5+usdhNev40/u7OQTQ2UGousIe0MrfrTDG1EMFRFZJBpqKq2i4lSYKs9wPHeK+w8F+M32\nIkef2Ytn1Wru2NrPh/oLdIY6CXg1Q+vNUNgSEVkkGmoqraDslJkuz3JyboQHjgT51XMl9j++E09P\nL5/+xGo+vn6GDYl+2oIpfNoM/6YobImIiAhQ72jNVrKcyo/yvSOhetD65Q5Id/LxO9bysfXTbEoO\n0BFqJ7BCD5V+O7RnS0RERKi6VTLVHKOFMX5wLMDPni5w+OfboLefO+/awGc2zbAhuYb2QIqAN9Ds\ncluKOlsiIiIrXM2tka3kGC2M8t2jPn786wKHf/wUdHZz12cu4nMXZdiYHFTQepvU2RIREVmhrLVU\n3Ppm+KO5Eb6xL8lvnp3m5G924Nm0mT+4cxUfXzfLuvgq2gJJ/Fo6fFvU2RIREVmharZGpppjrDDO\ng0ejPL1thpO/fA66e/j0J1bzifVzbEgO0h5q05OH50GdLRERkRXIWkumsXT47cMBfvZkjmO/epHg\nxZv57+7u5pNrS/RHh4j7YzpU+jzp356IiMgKU3PrHa1HTh7i/u1tPPPIDqJ9Xbz7cx/kD7eMc21n\njHQoTcQbVjfrAlDYEhERWUHmg9ZYYYQHXuzgmV/so2Oojyu29PPZy4a5Or2KdKhDQesCUtgSkZZl\njAkAHwbeC/QBRWAX8CNr7b5m1iayFM3P0RoujPC9o36efGwf0bY4A6uTvHdohotSbST8cQIev4LW\nBaSwJSItyRjzV8DHgSeA7cBjQAjYCPyfpv6d4t9ba3c1r0qRpWH+qcOZSobh/CgPHA7y6xcqRNvi\nrN/Uwe3XVbmpL0B3uJeYP6qnDi8whS0RaVU7rbV/8zq/95+NMb3AqoUuwhgzBPwlkLTWfmKh30/k\n7ajaGtnqHJPFcR4d9vDsHpc9v9rDpR+8hI+8o8aNfX66I73EFbQWhEY/iEirMsaY1/2uYK0dtdb+\n9m3e+D5jzLgxZterrt9sjNlvjDlkjPnzxvscsdZ+4e28j8hiqLhVZisZRguj/OAE/HS7hxcf3QGR\nKDdfa7ihz0N3uFdPHS4ghS0RaVVfAE41gtGNxpgL+ffZV4Gbz75gjPECXwZuAS4G7jLGXHwB31Pk\ngqu5NXLVOaZKE/zklMvPnzfseGQbhMPceufl3LHaw+rYKlKBRH2fFtqntRAUtkSkJVlrtwKbgCeB\n/xE4aYz5r8aYd12Aez8BTL/q8rXAoUYnqwLcD9z2Zu9pjLnHGLPNGLNtYmLqfEsUeUNVt8pMJcNo\nfpgfHK/xyHNedvzoeYgn+P1PXcOdFw3THanv0VJHa2EpbIlIy7LWzlpr/8laeyNwFbAP+AdjzNEF\neLt+4ORZr08B/caYDmPMPwBXGWP+4nfUeq+1dou1dktnZ8cClCfysqpbrT91mB/h/iNeHvq1w0s/\n+A2EQtxy1zX84aXDXNmxmqQ/rj1ai0BRVkRanjEmCdxKvdPUAfxgId7mHNestXYK+NM3dQNjtgJb\nh9atvaCFiZyt4lbJVLKMFkb5/lEvP3+2wv6fPw/pLm69+2o+s3mES9vrs7TU0Voc6myJSEsyxkSM\nMXcZYx4CDlCftfVfgFXW2i8twFue4pVPNw4AI2/lBtbah62196RSyQtamMi8qlslW8lxujDKg8cN\njz1bZd+jz0GyjVvvvprPXjzK5e2DdIbSBD2BZpe7YijSikirOgH8HPgK8InGPqqF9BywwRizFhgG\nPg3cvcDvKfKmVd0qmUqOkcIIPzwBjz7r1INWR5qP3n0ld20e5dK2VbSH2gl6dKj0YlLYEpFWtcZa\nOwf1SfLGmPXW2kMX4sbGmG8BHwDSxphTwF9ba//JGPMl4FHAC9xnrd39Fu+rZURZEDXrMFcrMFEa\n519PWn78lMveR7dB3wB/+Ll13LlhhqH4OtoCSQJedbQWm7HWNrsGEZG3zRhzK/B/AAFr7VpjzJXU\nw9EdTS7tdV2z5Sr75LOPN7sMWSZePutwlAeOWn7yVJX9j/4WOrv55Ocu43ObZ9iQXEV7sA2/8amj\ndQGFfant1totb/Rx2rMlIq3ufwGuA2YBrLUvAuubWpHIIqlZh2x1jvHiGA+fcHjstw77f7oNunr5\n5Ocu4+6LptmQWEV7IKXzDptIYUtEWl3VWjv7qmtLsmVvjNlqjLl3djbT7FJkGai5tcZThyN895jD\nQ0+67HnkGWjv4La7LuNTGyfZmOynLZjUeIcmU9gSkVa31xjzKcBjjFlrjPm/gGeaXdS56GlEuZBy\ntTyTxXEeOVnjx0857H34SUi18/E/uJI/uXSKy9uHSIfSBL1BdbSaTGFLRFrdl4BrABf4PlAC/l1T\nKxJZQNZaik6J0fww3zlW5ftPuOz916ehp49Pfv5q/vDiaTYkB0kFkwTU0VoSFLZEpCUZY/43AGtt\n3lr7H6y1VzV+/Lm1ttDs+s5Fy4hyIZTcMjPlWX5w3OWR39TY/8hT0N3DJ//g8jOb4dsCClpLicKW\niLSqm9/4Q5YWLSPK+So6JWbLGX41eowHn6iy/1+fhJ4+7v7jS/njSzJsSq2mPdim8Q5LjOZsiUir\n8hpj2jj3MTpYa199kLRIy3KtS9kpM12e4ZejJ7n/hU4O/fxxWD3Ep+7axJ0bplkbGyTlT+A3+ta+\n1OgrIiKt6iJgO69zZiEwtLjliCyMl4PWLI+PneS7O7t48qe7YdVqPnvnGu7cMMP65GB9vIM6WkuS\nwpaItKo91tqrml3EW6EJ8vJWOdahWCsxVZ7mkZMT/GBbkm2PPIOnq4vbbx/ijnWTDCXqe7Q03mHp\n0p4tEZFFoj1b8la8HLSmeOTkBA9uj7Ltp7sI9PVwyWXdfHrTOJuS/bQFUwR01uGSps6WiLSq/7vZ\nBYgslJp1KNaKTJWm+NGpKR7eEee5X+0n2tfFuo0d3H5dlYuS/bQH23SodAtQZ0tEWtW7jTGXnes3\njDFRY8znjTGfWeyiRM7X2UHrx42g9czP9xNpizO0vp2PXlvjpv4A6VAHQW8Qj9G38qVOnS0RaVVf\nBv6nRuDaBUwAIWADkADuA77RvPJE3rrXBK2X4jzziwNE2xOs31TvaN3UF6A70qug1UIUtkSkJTUO\nnP6UMSYGbAF6gSKw11q7v6nFvQ5tkJffZT5oTZYmeeTkDA+/EOW3P99DJJ1i/aYO7riuyk39QbrD\nvST8MQWtFmKsXZLntYqILFvXbLnKPvns480uQ5aQmnUo1ApMlib54fEMj2wPsv1nuwl1p9m4uZM7\nrq/y+2cFLZ9HvZKlIOxLbbfWbnmjj9NXS0REpIlqbo18rcBEaZzvHi3yk+1+djyxj2hfF5s2d/LR\nd5S5sS9AV7iHuD+qoNWC9BUTERFpkrOD1neOFvnRM7D7iV2EervZcFGa268tcUNfSB2tFqcFXxFp\nacaYS5tdg8jbUXNr5Gp5ThfH+NqBCt/7RZXdP/otno40Gzal+dh1ZW7oC9OjoNXy9JUTkVb3D8aY\nAPBV4JvW2tkm1yPyppScMtlKhp8O53lyV4gjO46Tvmozl1/Tz2cuG+bqjn46Qh1EfRFthm9x+uqJ\nSEuz1r4H+AywCthmjPmmMebGJpd1TsaYrcaYe2dnM80uRZqs4laZKk/x4PFpHng6wPOPvURqoJMr\ntvTzR1cM847O1XSHu4j5ogpay4C+giLS8qy1B4G/Av4D8H7g74wx+4wxH2tuZa+k43oEoOpWma1k\nePjEFD941s9Lj71EdKCXS67o5c5LRri8faA+Gd4b1GT4ZULLiCLS0owxlwN/DNwKPAZstdY+b4zp\nA54Gvt/M+kTOVnEqTJVnOJob5puPe9n7+AtE1qziyi39fPaaCd6R7qMj2EHYG1JHaxnRV1JEWt3f\nAy8AV1hr/4219nkAa+0I9W6XyJIwH7SO5Ib5570p9j62ndDgAJdf3cdnrp7gus4e0qE0EV9YQWuZ\nUWdLRFqatfZ9v+P3/mUxaxF5PWWnzFR5msPZEb6+r42nnpkkunE9l1/dy2e2THNdVzfpUJqwTx2t\n5UhhS0RamjFmJ/DqozAywDbgb6y1U4tflcjLyk6ZydI0B7PD/POeNM8+e5rRZ3by/j/+PT575Rhb\n0n3qaC1zClsi0up+DDjANxuvP934OUt9HMTWJtQkgmtdym6FydIUBzLDfG13F89vH2d0+z4iF23g\n01ec5up0Dx2hDnW0ljmFLRFpde+21r77rNc7jTFPWmvfbYz5bNOqkhVtPmhNlaY4lB3h2/vT7Nw5\nxakXD5O4aIjP3p7m2s4QnaFOIr4wXuNtdsmygBSjRaTVxYwx182/MMZcC8QaL2vNKencNGdrZXhl\n0Brl2wfaeWFnhqMvHCW1fhWfvT3N1tUl0sE0YQWtFUGdLRFpdV8AvmKMmQ9YOeALxpgo8LfNK+u1\nrLUPAw9fs+WqLza7FlkYjnUoOWUmihPsy4zy9T29vPTiOCe37Sc0OMAdH+7i9weyDMb6ifjC+BS0\nVgSFLRFpWcYYDzBkrb3MGJMEzKuO6/lOk0qTFajiVslV5xgvjvG1A4Zf/jbE/mefo31tL+/6xPV8\n7IpZ3tMDXaH1xHXW4YqiZUQRaVnWWhf4UuPXGZ2LKM1ydtD6+kH41bYy+7cdJb2uj0su7+HjV87w\n3p4EXaEuBa0VSF9tEWl1jxlj/j3wbSA/f9FaO928kmSlsNZScStkq3OcLozy1f1+frO9yP7tx+ha\n18NFl3Rz26VTvLsrQTrURcwfVdBagfQVF5FW9/nGz//mrGsWGGpCLbKCuNal6lbJVHOMFcb4lwN+\nnt5RZv/OUfo29rFxcydbL57kXV1x0uEu4v4oPqNvuyuRvuoi0tKstWubXYOsPK51qbhVspUcI4VR\nvrY/yNMvFNm34xQ9qzvYcFEnH948xTu74nTOLx0anw6WXqG0Z0tEWpoxJmKM+StjzL2N1xuMMR9p\ndl2yvJWdMrPlDCfmhvnH3WGeer7Avh2n6F2TZsOmTm66aIZ3doboCKWJ+aP4FbRWNIUtEWl1XwEq\nwLsar08Bf9O8cmS5q7hVZitZTuZHuG9PjGefz7L/+eN0DrSzZl07v7c5x7u7/XSFe0j64/g9fgWt\nFU5hS0Ra3Tpr7X8GqgDW2iKg72yyICpulXw1z/G5Ef5xV4Jnts9y+MXjdK3tqi8dXjbHB3t99EZ6\nSQbiBDwBHcMj2rMlIi2vYowJ0ziM2hizDig3tyRZjuaD1lhxjH/clWL781Mc33GMzg39rN/Uya2X\nZXl/T1AdLXkNhS0RaXV/DfwEWGWM+QbwbuCPmlqRLDuOdag4FTKVDCfnMhw5HufUqSxDVw+x6eI0\nH90wwpZ0uj7ewRfReAd5Bf1pEJGWZq19zBjzPHA99eXD/95aO7lY7984Fuj/ob5v7HFr7TcW671l\ncbjWpexUmC5Pcyh7mu8e7OTFp/fSuybN2g0d3Lp+lC3pNJ2hLqL+qI7gkdfQQrKILAchYAbIAhcb\nY953PjczxtxnjBk3xux61fWbjTH7jTGHjDF/3rj8MeC71tovAh89n/eVpcexDoVakfHiOLtnhvn6\n3i62vzhL75o0my/r4Y7NY1yb7iAd6iLqiyhoyTmpsyUiLc0Y85+AO4HdgNu4bIEnzuO2XwX+Hvja\nWe/jBb4M3Ej9icfnjDEPAQPAzsaHOefxnrLEONahWCsxWZpk7+wo397fw0svTXH0+SO872PXcOum\ncbak20iH0lo6lN9JfzJEpNXdDmyy1l6wTfHW2ieMMWtedfla4JC19giAMeZ+4DbqwWsAeJHfsVpg\njLkHuAdg1eCqC1WqLJCadSjVSowXT7Nndoxv7+tnz+4pjr54jMTqHu64dIp3difp0hE88iZoGVFE\nWt0RwL8I79MPnDzr9anGte8DHzfG/L/Aw6/3ydbae621W6y1Wzo7Oxa2UjkvZ4JWaZy9mdM8sL+P\nPbunOPL8EWL9nXzsw91c3xUnHepU0JI3RX9CRKTVFYAXjTE/56yRD9baf3uB3+dcz/Bba20e+OML\n/F7SJDXrUDyzR2uUb+4ZZN/u0xx/8SixVd184iM9fHQoQ1d4Q/2sQ48Po7Fu8gYUtkSk1T3U+LHQ\nTgFnr/8NACNv5QbGmK3A1qF1Os5xKaq5NQpOkfHiaXbNnOb++aC14xipNb18/NZObl2TYUOir97R\nMgpa8uYobIlIS7PW/nNjqOmgtXb/Ar7Vc8AGY8xaYBj4NHD3W7mBtfZh4OFrtlz1xQWoT85DZ/BM\n5wAAIABJREFU1a1SqBU5XRxj18wkD+xbxaH9ExzfM0z72l7uuCXNLaszrIv30hZs01mH8pZoz5aI\ntLRGt+hF6oNNMcZc2XhK8Hzu+S3gaWCTMeaUMeYL1toa8CXgUWAv8B1r7e63Wqsx5t7Z2cz5lCcX\nWM2tnQlaL01P8cC+AQ4dmOLw7hE613TxqY908NE1OTYl+0mH0oS8QQUteUuMtbbZNYiIvG3GmO3A\nh6gPFL2qcW2ntfay5lb2+q7ZcpV98tnHm12G0DiCp1bgdGGUF6dn+e6eXo4cmOTozpOk+tPcfksX\nt66ZZWOyl45QB0FvEK/xaPlQAAj7UtuttVve6OO0jCgira5mrc28qtOg/4qUNzQftMYKo7w4NcsD\nu3o4emiCY7uHaV/V2Vg6nGVDQkFLzo/Cloi0ul3GmLsBrzFmA/BvgaeaXNM5aYP80lFxKszVCowV\nRnhhOscP9vZy8vg0xw5O0j3UzW03prhxIMP6RC/twXYFLTkv2rMlIq3uz4BLqI99+CaQAf5dUyt6\nHdbah62196RSyWaXsqI51qkvHRbHeGkmyw/393Hi6DSH9o3TsyrF7Tem+P1VOTYme0mHOgj5Qgpa\ncl7U2RKRlmatLQB/2fgh8jvVrEPZKTPa6Gg9uK+P40cmObh3nO6+JB+9IcWNq7KsU0dLLiB1tkRE\nFomeRmyu+cnws+VMY+mwh+NHpji4b5Ke/mRj6TDL+kQPHcF2QgpacoEobImILBItIzbPfEcrU8ly\nIj/Gg/t6OXF0moN7x+nqifGR32vjhoEMG5I9dAQ76kHL41XQkgtCy4giIrKs1dwaJadMppLh2NwY\n/7IvyfDxGQ4dmKanP8mtH2rjQ/0zrE/0NJYOA3iNgpZcOOpsiUhLM8YMGGN+YIyZMMacNsZ8zxgz\n0Oy6ZGmoWYeyWyFbzXIqP863D8bZta/A3j0TdHZFuOWD7dywaobNqV46Q52EvaH6MTwaWioXkMKW\niLS6r1A/G7EX6AceblxbcrRna/FY7Jmlw9lyhpNzp/nGgSgv7i1z6MAkPb0xtt7Qzk2D02xK9jTm\naDU6WgpacoEpbIlIq+u01n7FWltr/Pgq0Nnsos5Fe7YWh8XiuPXN8DPlGY7mRrhvT4LndhTYuWOM\n7p4YH3pvmvf3T7E+0UlHsD6wVB0tWSgKWyLS6iaNMZ81xngbPz4LTDW7KGkOi62fdegUmS5PczA7\nwlf2tLNjV5YD+ydZszbF0IY0t6wZ5+JUD12hLkK+kA6WlgWlsCUire7zwKeAMWAU+ETjmqwwZzpa\nTplMOcPxuXHu39/BvgM5jhyeZvWaFKuHOnj/hgwXpbrrh0oraMki0NOIItLSrLUngI82uw5prvmg\nVXRKzJZnOTZ3mn/Z28auvTkOHZxiYFWS1UPtfGBDhvd0G3W0ZFGpsyUiLc0Y88/GmNRZr9uMMfc1\ns6bXow3yC8PaxtJhrch0eYaD2VHu293Bzj05Dh+aZvXaNoY2pvm9jbO8t9tDb6RbQUsWlcKWiLS6\ny621s/MvrLUzwFVNrOd1aYP8hWetpWZrFGslZiozHMmN8Y19aQ4czHL0yDSDa1IMrm3n/UPTvLPL\nQ0+kk1QgiU9PHcoi0jKiiLQ6jzGmrRGyMMa0o7/bVgTXutTcGnmnSKY8y6Hsae7f38nBw3lOnsiw\nbkMH/YNtvHv1FNemDd3hNAl/goA3gMeo1yCLR38hiUir+9+Bp4wx3wUs9c3y/2tzS5KF5lqXmq3v\n0cpWMhybm+C7hzo4eqLEqRMZBlenGBhM8e7Vk1zfaeiLdJEMJAn5Qni1qCOLTGFLRFqatfZrxpht\nwIcAA3zMWrunyWXJAprvaBWcEpnKLEdz4zxwMMWhYxWOH52hf1WCnv4k162a4poOS0+4k0QgUT/v\nEI+WD2XRKd6LyHLQDuSttf8VmDDGrG12QbIw5udo5Z0iM+VpjufGefBIkqOjltOjOVatSTEwmOT6\nwRm2pB36Ih0kAgmC3iAeo6AlzaHOloi0NGPMXwNbgE3Uj+nxA18H3t3Mus7FGLMV2Dq0Tlnw7XCs\nQ9WtkalkmS1Psz8zyfcO9HPiVIFTx2fp7o2TTIXZkJ7jolSBrnAHMX/8zB4tBS1pFnW2RKTV3UF9\nzlYewFo7AsSbWtHr0NOIb59rXRzXoeJWyFUzHMnN8NDhPoZHipw6PktXT4yOzihXDBZ5Z3eWjYlO\n0qFO4v4oAY8fr/FgUNiS5lBnS0RaXcVaa40xFsAYE212QXJhOdah5tbqk+ErGQ5kJnjwUDcnRioM\nn8rQ25+gPR3l8lVF3tE5y1AiTSrYRsQXxuvx1btaClrSRApbItLqvmOM+f+AlDHmi9SP6vlvTa5J\nLpCaW6PiVpmr5cmUZzmcm+SBA71MzTicODZDOOzjovURruqZ4coOw2B0LTF/lLA3hK8RtESaTWFL\nRFqatfa/GGNuBLLU9239z9bax5pclpwnay2OdepBq5onU6kHrUeOppnNWQ4fnCIWC3DNlSmu6pnh\n8jbDYLSbqD9af+rQ41XQkiVDYUtEWl4jXD0GYIzxGmM+Y639RpPLkrdpPmhV3Sr5WoFcNcvJ/Aw/\nO9nGyKyfo4cmiMb8XH15giu6ZrmszTAY6ybiixD2BtXRkiVHfxpFpCUZYxLGmL8wxvy9MeYmU/cl\n4Aj1wabSgl4ZtIpkKxlO5Kd49GSc49Nhho/PEAh6ueqyJJd3ZbiszTIY6ybqixLxhc8ELe3RkqVE\nnS0RaVX/AswATwN/AvwPQAC4zVr7YjMLk7dn/pzDqltjrpZnrpJlOD/N4yNRTudCTJzO0d2XoL/D\nsrkzw+aUw6pYNyFvqL50aLx4UNCSpUdhS0Ra1ZC19jIAY8x/AyaBQWttrrllydvhWvfMHq18rVAP\nWoUZHh8LM5wJMzJWpr0jSiRs2Nw5zcWpKmtj3YR9EaK+CH6PH68Ol5YlSmFLRFpVdf4X1lrHGHN0\nqQctDTU9t1pjtEPZKTNXnSNTmeVQdoYfHethdMpw/MgU8UQQj8fQHvfTFSozEGnDa7yEPEG8xgMK\nWbKEac+WiLSqK4wx2caPHHD5/K+NMdlmF3cuGmr6StbWj96putXGgdJZZiszHMjO8OjxbsZnvYwO\nZ+jsjhFPhujtDnBRZ5ZL2+IEvSFi/hh+b72jpTMPZSlTZ0tEWpK11tvsGuTts9bi4FK1NcpOhXx1\njlw1y9FchsdPdTExF+T06Axt7RGCQR89HYZN6QyXtBWJ+ftI+BOEfCH8xqflQ1nyFLZERGRRnf3E\nYamxdJitZjiSy/CLk2nGsiHGhjOk2iMEAl6621w2dWTZlCyyJtZGzB8n6A3iM16deSgtQWFLREQW\n1WuDVpZjc1meGOlgMh9kcnyO9nQUf8BLT6LCuvY5NiVKrIm3kQq2E/aG8HvU0ZLWobAlIiKLptZ4\n4rD8iqCV4YmRNsZzQcZOF4knQ/j8XrriFYba5tiQKLImniIVaKt3tM4KWhrzIK1AYUtERBactRYX\nF+espw5z1Swn5rL8ZrSNyXyQ8akqybYwfp+hO15mbSrPhkSRtfEUqUA7Id/8PC2Pgpa0FIUtERFZ\nUPOb4WtujYpTIVfNkavmODGX5cnTSabyAcYmHeKJED6vJR0tszqZZ128wOpYkmQgRdAbJOQJ4Dea\nEC+tR2FLREQWzNlBq+yUKTol5qo5TuazPDOeYDof4PQMJJMB/B6HjmiF1ckC6xNF1sSSpALtBL2h\nV5x5qKAlrUZhS0REFoRrXVzr1jfDuxWKtSJz1SzH57I8OxFnpuhnPOejLWHxehzawhVWJYoMxQsM\nROMkAykC3kB96VBBS1qYwpaIiFxQ890sx61RdisUakVKtSJT5Vn2ZUr87FgPuaLh9EiWeBKqFS9t\nCQ8Rv0MqUCEVCBL3RfEYD0FPoHHmYWOPlrKWtCCFLRGR82CMGQL+Ekhaaz/R7Hqa7Uw3qzGstFAr\nUKwVmChlOJCpMloMkyuaM+MdPB5DMmYYSBVZk8gzFPfSF2kj6A0R8UUbIx4aG+K1KV5alI7rEZEV\nyxhznzFm3Biz61XXbzbG7DfGHDLG/Pnvuoe19oi19gsLW+nSNz+otObWKDllCtUCuUqWXDXLcH6a\nnTMOJ/MRDkzGmZ0q0N4RJRD00Zm0rG4rsD4xx/q4oTeSJOANEvFFzszS8hgPHrSEKK1LnS0RWcm+\nCvw98LX5C8YYL/Bl4EbgFPCcMeYhwAv87as+//PW2vHFKXXpcqyDtZaqrVF16uccFmp5CrU8w/k5\n9mcDTJUCHJ+JAtDeGSPgc2kLV+hNFBmMFlgd89EbSRH0hoj6ogS8jScPPV48eLR8KC1NYUtEVixr\n7RPGmDWvunwtcMhaewTAGHM/cJu19m+BjyxuhUvb/N4s1579tGGRfDVPrlpguDDH3tkYM6UAJ2fD\n+DyWk6eK9PVFaItU6IuXWBUrsDrqoyeSIugJ1oOWx18/isdTP2C6vlVLaUtal5YRRUReqR84edbr\nU41r52SM6TDG/ANwlTHmL37Hx91jjNlmjNk2MTF14aptEsc61GztTMjK1wrkqjkylVkmyzkO5/K8\nMJVkPB/k+EwEr7GcGi7i9XlIR8v0J4qsihUYjProiSQJeAJE/I2g1Vg+VNCS5UKdLRGRVzrXd3b7\neh9srZ0C/vSNbmqtvRe4F+CaLVe97v1aQc06uK6DYx3KbpVSrUTByZOp5Jgulzg+Z9g/08Zc2cdY\nxk8wACdOFQgEvKwb8DGQmGMwVmAgGqI7lMTvfbmj5ff48Hp8ClqyrChsiYi80ilg1VmvB4CRC3Fj\nY8xWYOvQurUX4naLzrUuLpaqW6XqVqk4VYq1AvnaHJOlPDOVMsfmwhycjlGueZnIeIhHLKOnS0Rj\nQdb3uXRGi6yN5+mNhOkKJfF7AkR8kVd0tDyNgKWgJcuFlhFFRF7pOWCDMWatMSYAfBp46ELc2Fr7\nsLX2nlQqeSFut6jmRzo4bn0TfKlWJl+rn284WswxXipzIBNl70SCfMXH+KwhGjaMjZcJhXxs7Hfo\njRdZHcuzKhqhO1w/gifqjxH0BPB7fPiM9+UDpo2Cliwf6myJyIpljPkW8AEgbYw5Bfy1tfafjDFf\nAh6l/gTifdba3U0ss+le8bShW6VQK1Jw8uSqc5wuFJgqW47PxTkyHaXqepjNWZIxw/hUhWgsyMae\nEt3REn2REhsSEdqD8x2tMP75pcNG0NLSoSxHClsismJZa+96nes/An50od+v1ZYRrbW4uFhrqVmH\nqlul7FQoOgVy1TkmikUyVZeRQoSxuRB+r0ux6mVVusp4LkRvp2VDR47OUJnucIW18XrICnnDhHwh\nPXUoK4aWEUVEFkkrLSNaa7HYVwStilOl4pTJN4LWbMVhuBBmdC6EwVKqeemOl8gU/bRHKlzalWEg\nWqA/WmNdoo2gN0jCnyTiCxNs7NFS0JKVQJ0tEZFF0kqdLYvFxdb3aTUmw9fHPJSYrZQpOjUmSiFO\n50N4PZZCxUdXrMxUPkB7pMKG9jk6ghU6Q176Iwn83gBRX4ygN4jf2+hoGQ8eU/9vfgUtWc7U2RIR\nWSSt0tma72q51m10tWpU3CpFp0DJKZGtVhguBBnOhwGYK/toj1aYKfpJhKpsbJ+jK1SmK+ylP1o/\nfifqixHyBl9eOlTQkhVEYUtERM6Y36d1pqPV2BRfdkpU3ArZapnJkpfpcoCA16XqGHoSJebKPtrC\nVTa15+gJl0iH/PRHkvg9fiK+KEFvEF/jicP5sw5BQUtWBi0jiogI8MqgNb9Pq+pUKTklKm6ZuWqR\nmXKZqXIIr7Hkax66omWmigFS4QprYgXagxVSgSA9kQQ+j5+oL0bAG8BnvGc2w3u0R0tWGHW2REQW\niTFmqzHm3tnZTLNL+Z0sgLVgwcEFa3Fch5LjYAyEvC4110NnpEy24icWrDEUz9MRqtIeDNETqT91\nOL9HK+Dxnxla6qmnLAUtWVEUtkREFslS3rP16uXDqq1RditUnQolp0SuWqLsOGQqUKh5aQ9W6kHL\nXw9a7cEaqUCQrnAcn/ESnp8K35if5fF4gXrIUtCSlUZhS0RkhXtl0HJxXIeaW38CsepWqbhlatal\nUKsyXgrSFqyQqfiJ+GqsjuVJBVySgSDd4The4yPiixHwBPB7/Hg99eVDA2eWD0VWGu3ZEhFZoSwv\nLxWePeKh4lSoNDbFF2p5irUauWqF0WKAkNchV613tPoiRSI+D3F/kHQohqfR0fJ7fPg8XrzGg8+8\n3NESWakUtkREFkkz52zNj3NwrINj5zfB16g4VcpOmYpbpuQUma2UyFTKzFZgOB+m5HgZzsQ4fCpE\nLBGiXKoRDPnY3Fug7HhIh8rE/TBRzNEdTlJxymfe0+P1gFs7M+LBYzwYa7DGKnzJiqKwJSKySKy1\nDwMPX7Plqi8u1nu61j0zM8taS8WtUmt0sMpOmbJbouyUyFcrTJdLzNUcxktBxgohqo6HyXyQY8NV\nItEA5VKNgS5IhopE/DXi/hqDMYPjunSFY/g8PvyeAD6Pj4A3UC/ANPZo6WBpWcEUtkRElqEzg0kb\nQevMxne3Wt/87pQpO0VKtRJztSoz5TLZKkyXQ5wuBqk6HmYKAU6edoglQriOy0AXhH0OyVCV9mCF\n9QkXg6EzHCPoDRHw1KfDh7xBALzGgxcPHk99z5bRk4iyQilsiYgsI2d3ss5seLfOWQdJvxyystUK\nhVqNmYpDvupjpuJnqhTEWsiW/IxMe0il/NQcw2B7Ga/HEgvWaAtUWR93sEBnKIrf46s/cWg8BDwB\nrLX1pUNjMI2OlgKWrGQKWyIii2Qh92y51n1FF8u1bv08w/lulluh4pQp1ApkqhWylQpFx2WyFKTo\nBMhWfcwUAvi9LnNlH8NTXuJRDxbL2o4Cfq8l7HNIBKpsTNYwGLrD0cYcrQABT5CQLwSA17y8OV6D\nHkQ0+kFEZNEs5Jwt17q4roPTODC64lYpuxVKTpGSU2yErDnGS0XGiyWmyoaT+QizFT9TxQBT+SAe\nY5kr+xid8dGWMHg9ljXtBcL+eshqC1ZIB8t4jKEnHCXkDRH0BAl4goR9IQzU52p5vPg8vnrQMubl\n2VpaQpQVSp0tEZFloOiUqDbmYtXcKjVbo+yUKDplcpUKmWqFQg1GCmFKTpDZUn0De77ipVDxEQvW\nKFR8TMz56UpV8QD9ySJ+j0vY5xD31Yj5ayQDHnrCMfyNOVo+j5+gL4ilEbSMB2+jo4WClgigsCUi\nsixkKhnKTpmqW6HoVMhUKsxVKxQdw2QpyEwlRrbkJ1/xkSv7yJc9RIMu5ZoHxzUEvC7GWK4ayDBT\nDBD0ubQHK0R8NaI+h75IENd66AhFCXvDZ4KWz+PDb3x4jeflEQ/n6GiJrGQKWyIiy8BwfoJ8rcp0\nucZkOUDZ8TFXDTGRDzJb9FOo1IeLFvIVKuUysXiQTNUCDl0pi9dj6U8UyVV8JEJV+iNFIj6HqM+l\nJxzBsZaOYBS/x/+KoBX0BPCcHbQaG+NfTV0tWckUtkREloGT+SKjxRAnsgkctx5sTs2GyWVLVEpF\n/AEv5XKNWtUhkQpTLtcIBLx0pSxhv8OqRIF8zUd7qEJ/tEjU5xL2+kiHogC0+SP4jA+/x0/AG8Rn\nvC+PePD4oPEE4is6WqDlQxEUtkREloV/ejbNM08cYfOlPYyO5Fi3oYOTJ04xsCqFay3lco1NQ0Ey\nRT/GQNDn0hMv4TGW/liRsuOhLVihJ1wi5jN4jZeecISK65IKxvF7AngwhH0Rgt4AfuPD4/FirYvf\n48e17plQZTTAVOQVFLZERBbJQo5+eOYnu1h7+SDBsJ9EMsjkRJ54PEi5XCMU8pHuijE6VSUcNvQn\nCwR9Lj6PZV1ijkzVT8JfY03cxeAj7veTDEQAQzIQxtdYOvR76p2tM8uG1uJYFz/gWOdM6DozZ0tE\nADDW2mbXICKyolyz5Sr75LOPX9B7fuyHx4iEDcPDeVJtYY4cnKSnL0E0FiQZq39Md6xEIlQl6nPw\nGMtQvETVdYn5/ER8PkJeH16PIe6L4jFeAt4gYV8Yv/HhWAefx4cHc2asg89T/+91xzp4jZf57yfq\nbMlKEfaltltrt7zRx6mzJSKyDMxMF7BtYSLRAB0JWP2uFD2xEo6tEvXXAEj4q6QCVcqOh4vbosxV\ng42QFTgz/T3gCWKMafzajzEG17qEvEGM8RDw+KlZB5/xnnlvb+PXClki56awJSKyDFy50ce6tlku\nbXM4kDFE/TWCHpfhQphLUhUslqrrMhhNnOlIWR+0B1NA/SnCkDcMQKDxpKGlvhxYn5sFfo+/3uE6\nK2iJyBtT2BIRWQZuGJyiPRiiM9RG3D9HZ6iLXDXH5pQl4U9TdkuAIeaLU3HLGAzpUARP41id+T1W\nxpgz+7EC3iDY+mHW/rM6WtZadbFE3gKFLRGRZaA/EmO0kGdNLEbYG8EYaAu2UawVCXpDBL0hvMaD\nBVxcEv54PVg1niB0sfiNj4A3QNWtYjwevHiwxuKbD2KN91LQEnlrFLZERJaBtmCKjlA7XuMl4A0Q\n90WZqWSI+xNEfPXlwXytQNgXIuQN4vV463utGpvavcac2Y8138UCXvFUoVfLhyJvi8KWiMgysDq2\nirJTJlPJkfDXHz/sDneSreaI+aMUa0W6w504jRBlAS+eM+2q+W7V/H4s7csSuXA0CEVEZJEYY7Ya\nY+6dnc0s2HukAokzR+oARHyRxnt7KDtlvKbe0fIZ75mANf/zfDdLRC4shS0RkUVirX3YWntPKpVc\nkPsHvUEC3sArrs13qELeIMHG8TpnO3v/lbpZIgtDYUtERERkASlsiYiIiCwghS0RERGRBaSwJSIi\nIrKAFLZEREREFpDCloiIiMgCUtgSERERWUAKWyIiIiILSGFLREREZAEpbImIiIgsIIUtEZHzZIy5\n3Rjzj8aYHxpjbmp2PSKytChsiciKZoy5zxgzbozZ9arrNxtj9htjDhlj/vx33cNa+6C19ovAHwF3\nLmC5ItKCfM0uQESkyb4K/D3wtfkLxhgv8GXgRuAU8Jwx5iHAC/ztqz7/89ba8cav/6rxeSIiZyhs\niciKZq19whiz5lWXrwUOWWuPABhj7gdus9b+LfCRV9/DGGOA/wj82Fr7/LnexxhzD3APwKrBVRes\nfhFZ+rSMKCLyWv3AybNen2pcez1/BtwAfMIY86fn+gBr7b3W2i3W2i2dnR0XrlIRWfLU2RIReS1z\njmv29T7YWvt3wN+94U2N2QpsHVq39jxKE5FWo86WiMhrnQLOXusbAEbO96bW2oft/9/e/cbKUZVx\nHP/+oNY2RFqCgWgLaU1rQxsjlbaQEE1Mmor4DyziRXkBqaDSSvQFQRNfNEEjKsFaVEopUo16kTRE\nQIKNrRGCqbQFbrF/UkVeSIu1atoq5Y9SHl/sWZm77Oy/e2d3Z+/vk0w6c+bMmeeZ3eQ+mTPdibh2\n+vRpYx3KzErExZaZ2RvtAOZKmi1pMjAEPNDjmMyspFxsmdmEJmkY2AbMk3RA0oqIeBVYBWwG9gH3\nRsSecTjXRyStP3r02FiHMrMSUUTuYwhmZlaA8xYtjN89/tteh2FmYzR10vQnImJRs36+s2VmZmZW\nIBdbZmZd4mlEs4nJxZaZWZf4fyOaTUwutszMzMwK5GLLzKxLPI1oNjG52DIz6xJPI5pNTC62zMzM\nzArkYsvMzMysQH4RtZlZl1RfRA28OHXS9H2ZXdOAYw3Ws21vBf7RwemzY7Tbp7a91e16uXQaf6P4\nmu2v194oB38G7cfXbP8gfgYAc1s6a0R48eLFi5cuLsD6vO166zVtO8fjnO30aRRvo+2cXDqKv5Uc\nWo3fn4E/g6JyyFs8jWhm1n0PNtiut17bfzzO2U6fRvE22s7Lq1PNxmg1/npt/gxa48+g/fH8bkQz\nszKRtDNaeBdbvyp7/FD+HMoeP5QvB9/ZMjMrl/W9DmCMyh4/lD+HsscPJcvBd7bMzMzMCuQ7W2Zm\nZmYFcrFlZmZmViAXW2ZmZmYFcrFlZmZmViAXW2ZmA0DSJZLulHS/pGW9jqcTkt4h6S5Jm3odS6sk\nnSLpR+naf7rX8XSijNc9qwzffRdbZmY9JumHkg5L2l3TfpGk/ZKekfTlRmNExC8i4hrgKuCTBYZb\n1zjl8GxErCg20ubazOXjwKZ07T/a9WBztJNDv1z3rDbj7+l3vxUutszMem8jcFG2QdLJwPeBDwLz\ngSskzZf0Lkm/rFnOyBz61XRct21k/HLotY20mAswE3gudTvRxRib2UjrOfSjjbQff6+++035RdRm\nZj0WEY9KmlXTvAR4JiKeBZB0D/CxiPgG8OHaMSQJuBl4OCKeLDbiNxqPHPpFO7kAB6gUXCP00Q2M\nNnPY293ommsnfkn76OF3vxV988UwM7NRZvD6HROo/FGf0aD/F4ClwGWSPldkYG1oKwdJp0taByyU\n9JWig2tTXi73Acsl3c74vLuvSHVz6PPrnpX3GfTjd38U39kyM+tPqtOW+8qPiFgLrC0unI60m8M/\ngb78Y0lOLhFxHLi628F0KC+Hfr7uWXnx9+N3fxTf2TIz608HgLMy2zOB53sUS6cGIYeqQcil7DmU\nNn4XW2Zm/WkHMFfSbEmTgSHggR7H1K5ByKFqEHIpew6ljd/FlplZj0kaBrYB8yQdkLQiIl4FVgGb\ngX3AvRGxp5dxNjIIOVQNQi5lz6Hs8ddSRO70uZmZmZmNke9smZmZmRXIxZaZmZlZgVxsmZmZmRXI\nxZaZmZlZgVxsmZmZmRXIxZaZmZlZgVxsmZlZX0vv7htJyyFJBzPbk3sdXzdIminpfkkXZ3J/QdL+\ntH53znEnpd+pmlXTvk7S9ZIWSbqjGzlMZP6dLTMzKw1Jq4EXIuKWmnZR+Zv2Wk8Ca0DSpPSDnGMZ\n4zvAloh4KNP2GLAqIkZaOPZQRHyzGg+VV9+cGxGHJD0KXB4Rh8YSo+XznS0zMyslSXNY8eE0AAAD\nQ0lEQVQk7Za0DngSOEvS0cz+IUkb0vqZku6TtFPSdkkX1BlvkqRb0/6nJX0mtS+VtDUdv1/SjzPH\nLJb0iKQnJD0s6czU/pikr6dCZpWkuZIeT2PfVI1T0rCkD2XG+7mki2viEnAJ8Osm1+NNktakc+yS\nVH1B9jCVV9tULQV2Z4qrh4DLG41tY+Niy8zMymw+cFdELAQONui3FvhWRCyiUlhsqNPnWuBwRCwB\nFgMrJZ2d9r0HWJnOd46kCyS9GfgusDwizgN+AtyUGe/UiHhfRKwBbgNuSWP/LdNnA3A1gKTT0nk3\n18Q1J8X1n0YXAvg8cCCd43zgi5JmRMR24C2S5qV+Q1QKsKqdwHubjG1jMKnXAZiZmY3BnyNiRwv9\nllJ5z151+zRJUyPipUyfZVQKqepdoGnA3LT++4j4K4CkEWAW8DKwANiSxj2ZyvRc1T2Z9fOB6h2r\nnwFfS+u/AW6TdDpwBZX3/Z2oif1twN9byHEZMEfSlWn7VCqF2sEUy5Ckm1McX8ocdxh4ewvjW4dc\nbJmZWZkdz6y/BiizPSWzLmBJk7tDAq6LiK2jGqWlwCuZphNU/n4KeDoi8u4KHc9p/7+ICEk/BT4F\nXJX+rfUSo3PJI+CzEfFInX3DwCZgF7AtIo5k9k1J57CCeBrRzMwGQno4/kh6Puok4NLM7i1UpgEB\nkHRunSE2A9elB8iRNE/S1Aan3AvMkLQk9Z8saUFO3+2ZeIZq9t0N3AC8HBH76xy7H5jdII68+M+R\nNAUgIvYA/wVWM3oKEeCdwO4WxrcOudgyM7NBciPwK2Aro6f0VgIXpgff9wLX1Dn2DuBPwIik3cDt\nNJgBiohXgMuAWyXtAp6iMl1Yz/XAjZK2A2cAxzLjPA/8kUrRVe88/wKek9Ss4PoB8BfgqRT/96hM\nbVYNU5kWfbDmuPdTeUjeCuKffjAzMyuYpFOAF9O04ZXApRGxPLPvD8C7I+LfOcd/AlgQEasLiGsL\ncGE//mzGoPAzW2ZmZsVbDKxJ05tHeP1/IH4AuBP4dl6hlWyi8sD+eDsbuMGFVrF8Z8vMzMysQH5m\ny8zMzKxALrbMzMzMCuRiy8zMzKxALrbMzMzMCuRiy8zMzKxA/wOyqtq7q0MI+gAAAABJRU5ErkJg\ngg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "sim.obs.peek()\n", "print(sim.obs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Spectral analysis\n", "\n", "Now that we have some simulated CTA counts spectrum, let's analyse it." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:gammapy.spectrum.fit:Running SpectrumFit\n", "Source model PowerLaw\n", "\n", "Parameters: \n", "\n", "\t name value error unit min max frozen\n", "\t--------- --------- ----- --------------- ---- ---- ------\n", "\t index 2.300e+00 nan 0 None False\n", "\tamplitude 2.500e-12 nan 1 / (cm2 s TeV) 0 None False\n", "\treference 1.000e+00 nan TeV None None True\n", "Stat cash\n", "Forward Folded True\n", "Fit range None\n", "Backend sherpa\n", "Error Backend sherpa\n", "/Users/deil/Library/Python/3.6/lib/python/site-packages/gammapy-0.7.dev4616-py3.6-macosx-10.12-x86_64.egg/gammapy/stats/fit_statistics.py:55: RuntimeWarning: divide by zero encountered in log\n", " stat = 2 * (mu_on - n_on * np.log(mu_on))\n" ] } ], "source": [ "# Fit data\n", "fit = SpectrumFit(obs_list=sim.obs, model=model, stat='cash')\n", "fit.run()\n", "result = fit.result[0]" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Fit result info \n", "--------------- \n", "Model: PowerLaw\n", "\n", "Parameters: \n", "\n", "\t name value error unit min max frozen\n", "\t--------- --------- --------- --------------- ---- ---- ------\n", "\t index 2.229e+00 4.879e-02 0 None False\n", "\tamplitude 2.774e-12 2.984e-13 1 / (cm2 s TeV) 0 None False\n", "\treference 1.000e+00 0.000e+00 TeV None None True\n", "\n", "Covariance: \n", "\n", "\tname/name index amplitude\n", "\t--------- ------- ---------\n", "\t index 0.00238 -1e-14\n", "\tamplitude -1e-14 8.9e-26 \n", "\n", "Statistic: 460.387 (cash)\n", "Fit Range: [ 5.01187239e-03 5.01187225e+02] TeV\n", "\n" ] } ], "source": [ "print(result)" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaMAAAEOCAYAAAAkF3jEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8ZPV97//Xd2ZUZkZdmlHbImlVtoEBU0w1xcACu0sP\nYAyLjUniGzu5SRwnTpzrJL6OfX9x8rOT+Maxsc0ui5eOvUs1zYDpLL1J2pVWddT7aNRmPvePGS1C\nqOtIMyN9no8HD+vU+QzH0ptzzrcYEUEppZSKJlu0C1BKKaU0jJRSSkWdhpFSSqmo0zBSSikVdRpG\nSimlok7DSCmlVNRpGCmllIo6DSOllFJRp2GklFIq6jSMlFJKRZ0j2gXEi5ycHCkqKop2GUopFVcO\nHjzYISKe2fbTMJqjoqIiXnvttWiXoZRSccUYUzeX/fQxnVJKqajTMFJKKRV1GkZKKaWiTsNIKaVU\n1GkYKaWUijoNI6WUUlGnYaSUUirqNIyWQX9/PyIS7TKUUipmaRgtg9bWVt577z26urqiXYpSSsUk\nDaNZGGN2GGN+2tvbu6jzDA8PU1tbywcffEB/f79F1Sml1MqgYTQLETkgIn+Ynp5uyfkGBwepqqqi\nurqaQCBgyTmVUireaRhFSV9fH++//z5HjhxhZGQk2uUopVRU6UCpUdbZ2UlXVxder5e8vDwcDr0k\nSqnVR//yxQARobW1lY6ODnJzc8nNzcVm05tWpdTqoX/xlpq/k6SO9+a0azAYpLm5mXfffZf29nZt\nDq6UWjU0jJbay//F2kdvZMPLf4urp2pOh4yOjlJfX897771Hd3f3EheolFLRp4/pltrpf0Zn3yDp\n7+1h03N/TE/uaTRX7CKQXjbrocPDw9TU1OByuSgoKMCqFn1KKRVrjD4KmpsTTzxRFjrT66FDh+jv\naMZbez+5NffgGB2gO+90fOU3EUjfMOfzpKSkUFhYSEpKyoLqUEqp5WaMOSgiJ862n94ZLZNQgpuW\n8htoK76c3Jr78NbcS2bLLXTnn0Vz+Y0MpZXMeo6BgQEqKytJT0+nsLAQp9O5DJUrpdTS0zBaZqGE\nFHwVu2grvgJvzb3k1t5Hpu9ZuvI/i6/iRoZSi2c9R29vL729vWRmZlJQUEBycvIyVK6UUktHwyhK\ngomp+DZ+kbaSK8mtuQdvzf1k+p6lu+BsfOU3MpS6ftZzdHd309PTQ1ZWFgUFBSQmJi5D5UopZT0N\noygLJqbRvPFmWkuuIvfwPXhr7yez+Xd0FZ6Lr/wGhlPWzXi8iBztOJuTk0N+fj4JCQnLU7xSSllE\nwyhGBBPTad70ZVpLribv8F14jvyarKanI6F0I8Mpa2Y8XkRob2+ns7MTj8ejozkopeKK/rWKMcGk\ndJo2/yGtG64m9/BdeGt/Q1bTU3SuOR9f+RcYcRfOeHwoFKK1tZX29na8Xi+5ubkaSkqpmKdNu+do\nMU27u7q6aGpqWtCAqI6hLvIO34nnyH6MjNG55sJwKLny53S83W4/Gkp2u33en6+UUosx16bdGkZz\ntNAw6g2MMjA8RkF6Mt3d3bS2tjI4ODjv8ziGOsk7tA9P3QGMhOhYeyEtZV9gxJU3p+Ptdju5ubl4\nvV4NJaXUstEwsthCw+jHTx/i3x6v4pJj8rnlzBKOWZNOX18fra2t9PX1zft8CYF28g7tI6f+IRCh\nc+02fGXXM+rKndPxGkpKqeWkYTQFY0wJ8HdAuohcZYy5DLgE8AI/FpHfTnfsQsOoqSfAbc/Xsu+V\nBgaGxzilOItbzizh3I1ehoYCtLa20t3dPe9BURMCbeRV/4qc+ocB6Fx3UTiUnN45Ha+hpJRaDjET\nRsYYO/Aa0CQi2xd4jl8A24E2Edk6ads24EeAHbhVRL4/h/PdKyJXTVjOBH4gIjdPd8xi3hkB9A2N\nctcrDfzy+Vqae4co8bi5+YxirjxhDSY0RmtrK52dnYRCoXmdN2GwlfxDvyK7/hEwho51F9NS+nlG\nnZ45Ha+hpJRaSrEURn8BnAikTQ4jY4wXCIhI/4R1pSJyaNJ+ZwEDwJ6JYRQJuirgfKAReBW4jnAw\nfW9SKV8SkbbIcZPD6F+BO0Tk9em+x2LDaNxoMMTD7/j42XM1vNvUR5Y7kS98Zj03nrqejGQ7bW1t\ntLW1EQwG53XexMEW8qrvIKfhUcTY6Fi/nZbS6xhNzpnT8drQQSm1FGIijIwxa4DdwHeBv5gijK4G\nvgJcLCJDxphbgMtF5OIpzlUEPDgpjE4F/kFELowsfxNARCYH0eRz3Rt5TGeA7wOPi8gTMx1jVRiN\nExFequni1udqePLDNhIdNq48oZCbzyihJMdFe3s7ra2tjI6Ozuu8iYM+8qvvILvhUcQ4aF+/g5bS\n6xhLzprT8eOh5PV6tUm4UmrRYmWg1B8C3wBSp9ooIvcYY4qBO40x9wBfInyXM1eFQMOE5UbglOl2\nNsZkEw7G4yPB5Qc+B6RH7sh+MsUxO4AdpaWl8yhrdsYYTt2QzakbsjncPsCtz9Vy/+uN7HulgXM3\nernlzBJO2bqVrq4uWlpaGB4entN5R1z51H3q6/hKryO/ei/eIw/gqXuQtqKdtJZey1hS5ozHB4NB\nfD4fra2t2k9JKbVsluzOyBiznfAdz/8wxpwNfH26d0bGmDuBi4ENItI+zT5FfPLO6GrgQhH5cmT5\nBuBkEfmald8FrL8zmkrnwDB7X6pnz4tH6PSPsLUwjVvOLOGirXkM9PXS0tJCIBCY1zmTBhrJr76d\nrMYnCdkTaS+6lNYN1zCWlDGn4202Gzk5OeTl5ekwQ0qpeYv6YzpjzPeAG4AxIBlIA+4XkS9M2u9M\n4L+Ag0C/iHx1mvMVYdFjuoVYjjAaNzQa5IE3mrj1uRoOt/spSE/mptOLuPbkdcjwID6fD7/fP69z\nJg3Uk1+1l6ymJwnZk2grvpzWDX9AMHFuE/YZY46Gkg7IqpSaq6iH0aRizmaKOyNjzPHAPsLNq2uB\nvUCNiHxrinMU8ckwchBuwHAe0ES4AcPnReQ9q7/DcobRuFBI+F1VGz97tpYXazpJSXJwzUlr+eLp\nRaQ7wo/T+vv7Zz/RBMn9deRX7SGz+XeE7Mm0FV9B64arCSamzel4YwxZWVnk5eXp1BVKqVnFSxid\nDvSJyDuR5QTgJhH52aT99gFnAzlAK/BtEfl5ZNvFhN9N2YFfiMh3l+I7RCOMJnq3qZdbn6vhwbd9\nCHDR1jxuObOEDZkOfD7fvDvQJvfXkl91O1nNvyPocIVDqeRqgolTvt6bUmZmJnl5ebhcrvl9GaXU\nqhFTYbQSRDuMxjX3BNj9whF+9XI9/cNjnFyUxZfPLObU9Sm0trTQ29s7r/Ml99VQULWHTN+zBB1u\nWkuupK3kKoIJc5/aPC0tjby8PFJT5x5kSqnVQcPIYrESRuP6h0a569UGfvn8EZp6ApTkuPnSGcVc\ntCmLno42enp65nU+Z+9h8qt2k9nye8YSUmgruYrW4isIzSOUUlJSyMvLIz19bu+hlFIrn4aRxWIt\njMaNBUM88m4Ltz5Xw1uNvWS6ErjhM+u56rhcRvs7FxBK1RRU7iGj9XnGElJpLbmatpIrCDnm/ijO\n6XSSl5dHZmYm4a5cSqnVSsPIYrEaRuNEhFePdPOz52p44oNWEuw2Lj+ukC+cXIBrtHfeoeTqqSK/\n6jYyWl9iLCGN1g1/QFvx5YQczjmfIykpidzcXLKzs7HZbPP9SkqpFUDDyGKxHkYT1bQP8PPf13Lv\nwUaGx0KcU+HhxlMKWZsYmPc7JVfPhxRU7ia97WVGE9Np3XAN7UWXziuUHA7H0VEddKghpVYXDSOL\nxVMYjevyj7D3pTr2vHiEjoERNuenceMphRyfHWKgb56h1P0BBZW3kd7+KqOJGbSUXkv7+p2IY+7N\nu202Gx6PB6/Xq32VlFolNIwsFo9hNG5oNMiv32ji1t/XcqhtgLy0ZD5/UgFn5htCw/PrPOvuepeC\nyt2kdRxkNCmTlg3X0l60E7EnzfkcxpijzcKdzrnfYSml4o+GkcXiOYzGhULCM1Xt/Oy5Gl443Ik7\n0c7lx+VxwfoE3AzN61zuzncoqNpNWsfrjCZl0VJ6He3rt88rlECbhSu10mkYWWwlhNFE7zb18vPf\n13LgrWZCIpy/0cPFG5JY4xyb13lSOt+ioHI3qZ1vMpKUTUvZ5+lYdwlin99jOJfLRW5urrbAU2qF\n0TCy2EoLo3G+3gC3jXeiHRrj+DVpbC9zckw22OcRCikdb4RDqettRpJzaCm7no61F807lBITE8nN\nzSUnJ0db4Cm1AlgSRsaYY+fwWaMi8sF8iotHKzWMxg0Mj3H3qw384vlaGrsDrM1MZnuZizMLE0hy\nzDGUREjteIOCyttI6X6XkWQvvrLr6Vy3DbHNb8Rvu91+tLGDjhauVPyyKoz6gTeAmf4arRWRonlX\nGGdWehiNGwuGeOy9Vn76XA1vNfSQnuzgwg0uLihJIjN5js2yRUjtOBgJpfcZduZG7pS2gW1+cyON\nD8yam5urjR2UikNWhdGzInLWLB806z4rwWoJo3Eiwmt13fzs2Roe/6AVh83w2fVOLil1si59jncq\nIqS1v0Z+5W2k9HzAsDMPX/kX6FxzwbxDCcKNHXJzc0lLm9sI40qp6NN3RhZbbWE0UW2Hn1/8vpZ7\nDjYwNBrihPwktpe5ONabOLfGBiKktb1MQeVu3L2VDLvy8ZXfQGfh+WCbfydYp9NJbm4uWVlZ2thB\nqRhnaRiZ8G/8VqAACADviUjnoquMI6s5jMZ1+0e44+U6bnsh3Im2KMPBjjI3p69LJsE2t1BKb32R\n/KrduHurGXIX4iu7ga7C8xYUSgkJCXg8Hjwej06NrlSMsuoxXRHwDWAb4cnv2gnP2loG9AA/AfbK\nKri90jD6yPBYkN+80czPnjtMdZufrGQbF5W5uKDERUriHFrAiZDe+gIFlbtx9R1iyL0GX/mNdBWe\nA2b+oWSz2cjOzsbr9eqEf0rFGKvC6G7CU4I/IyKhSdvygeuBDhG5bXHlxj4No08SCXei/e9nDvFi\nTTfJdsO5xU4uKXORlzKHOxUJkdHyPPlVu3H11RBIWYev/Aa6C85eUCgBpKenk5ubq51olYoR+s7I\nYhpGM3u/uY+fPF3JQ++2IQInFyZxaYWb8uw59DOSEBm+5yio2o2z/wiBlPX4KnbRnX8WmIX1NdL3\nSkrFBqvfGV0BPC4i/caYvwFOAP5ZRN5cfKnxQcNoblp6h/jJUx9wzxst+EdCVGQnsLPczUmFSbN3\nopUQmc3PkF+1B+dAHYHUYprLb6Qn/8wFh5K+V1IquqwOo7dF5FhjzGnAvwD/BvyViHxm8aUuH2NM\nCfB3QLqIXGWMuQy4BPACPxaR3053rIbR/AwMjfLL333I3lebaPUHyXPb2V7u4pwiJ8mOWYJFgmQ2\n/46Cyj0k+xsYTCvBV76LnrwzYIF3OdpfSanosDqM3hCR440x/0y4Jd0d4+tmOCYZeBZIAhzAvSLy\n7bl/hY+d6xfAdqBNRLZO2rYN+BFgB24Vke/P4Xz3ishVE5YzgR+IyM3THaNhtDAjo2Pc/UIle15p\noqpzlJQEwwUbXFxc6iLTOct7IQmS1fQ0+VV7SPY3MphWSnPFLnpzT1twKAGkpqaSm5ur06MrtQys\nDqOHCbem2wacCPiBV0XkUzMcYwC3iAwYYxKA3wN/JiIvTdjHCwREpH/CulIROTTpXGcBA8CeiWFk\njLEDVcD5QCPwKnAd4WD63qSSviQibZHjJofRvwJ3iMjr030fDaPFGRkZ4bGD1ex9zccrTcPYbXDm\nOic7yl2sn60TbShIVtOT4VAabMafXoavfBe9uacuKpSSkpLwer06Dp5SS8jqMEoBLgbeFpEPjTEF\nwKdE5JE5FuMiHEZfEZGXJ6y/GvgKcLGIDBljbgEuF5GLpzhHEfDgpDA6FfgHEbkwsvxNABGZHEST\nz3Vv5DGdAb5P+H3YEzMdo2FkjcHBQV55v4Y732jn6SMBhoPCp3IT2Vnu5lO5s3SiDQXJbnqc/Krb\nSRr04U+voLliF33eUxYVSna7nZycHJ30T6klEBOt6SJ3LgeBUsLvZP56in2+AZwG3AN8FThfRAam\n2K+IT4bRVcA2EflyZPkG4BQR+eo09WQD3yV8J3Ur4Tu8XYTvqN4UkZ9MccwOYEdpaekt1dXVc//y\nakY9PT18WFPP/ve7efjQID1DIdalO9hR7uLMtU4S7DOF0hjZjb8lv2ovSYEWBjI24avYRZ/npEWF\nEkBGRgZer1ebhitlkZgIownFZAAPAF8TkXen2H4n4TuvDSLSPs05ivhkGF0NXDgpjE4Wka9Z/R30\nzsh6IkJ7ezv1jc08XdvPgcpB6vvGyEy2cVGpiws2uEidoROtCY2S3fAYedV7SQq0MZC5meaKm+jP\n+fSiQ0mbhitljZgKIwBjzLcBv4j8YNL6Mwl3rD0I9M9wV1OERY/pFkLDaOkEg0F8Ph+tra282TLM\n/io/b7WOkGQ3nFMUfq80UydaExolu/4R8qvvIHGonYHMrZFQOn7RoeRwOI42DdepLJSav6iHkTHG\nQ3iuox5jjBP4LfB/ROTBCfscD+wj3Ly6FtgL1IjIt6Y4XxGfDCMH4QYM5wFNhB+3fV5E3rP6+2gY\nLb3h4WGampro7u7mSM8oB6r8/L5+iGCkE+3Ocjcbc6Z/p2OCI+Q0PEJe9R0kDnXQn3UszRU3MZBz\n3KJrM8aQmZmJ1+vF7XYv+nxKrRZWDQdUCPwfoBB4BPg3ERmLbLtPRK6c4dhjgd2EW7bZgLtF5J8m\n7XM60Cci70SWE4CbRORnk/bbB5wN5ACtwLdF5OeRbRcDP4x8zi9E5LuzfemF0DBaPgMDAzQ0NDA4\nOEhXIMgjhwb57eFBBkaF8qwEdla4OXmGTrQmOEJO/UPkVf+KxOFO+rOPo7liFwPZ0zb+nBe3243X\n69Up0pWaA6vC6DHgAPAScDPhkbt3ikj3bP2MVhoNo+XX2dlJU1MTo6OjDI2FeOpIgAerBmn1B/G6\n7ewoc3FOsRPnNJ1oTXAYT92D5B3aR8JwF305x9NcfhP+7GMsqU9Hd1BqdlaF0ccCxxizC/g6sBO4\nT0ROsKLYeKBhFB2hUIiWlhZaW1sJhUIERXi1KfxeqbJzFPeETrRZ03SiNWNDeOoOkHfoThJGuunL\n+TTNFTfhz9piSY3jozt4vV5cLpcl51RqpbAqjN4HjheR4QnrtgH/CbhEpMCKYuOBhlF0jYyM0NTU\nRFdX19F1lZ0j7K/080rTMDYDZ6xLZke5m6KMqRsa2MYCeOr2k3voLhJGeuj1nExzxS4GMzdZVmdK\nSgper5eMjAx9hKcU1oXRXxEeaeF3k9afCPyLiJyz2ELjhYZRbPD7/TQ0NOD3+4+uaxkY46HqQZ6q\nDTAUFI71JrKzws1x03SitY0F8Bz5DXmH7sQx2kev9xSaK25iMKPCsjoTExPJycnRR3hq1Yt6a7qV\nRsMotkx8nzRuYCTEb2sGebh6kO6hEOvSIp1o103didY2Noi39tfkHr4bx2gfPbmfobn8JgIZ5ZbV\nabPZjrbC00d4ajWyejigdYRHRygiPOgpACJyxSJqjCsaRrFn8vukcaMh4fn6IfZX+anrHSMj0on2\nwhIXqUmfbOxgG/XjrX2A3Jp7cIz205N7Os0Vuwikl1parz7CU6uR1WH0JrAHeAc4+lsvIk8upsh4\nomEUu0ZGRmhsbKS7u/tj60WEt9vC75XebB0h0Q7nFrnYXu4if4pOtLbRAXJr78d7+B4cY366887A\nV7GLQNoGS+vVR3hqNbE6jF4RkZMtqSxOaRjFvv7+fhoaGggEAp/YVt87yoGqQZ6tDxAMwUkFSeys\ncLMxO+ETdyn20QG8NfeSW3Mf9jE/3fln0Vy+i6G0Ykvr1VZ4ajWwOoxuANYDjwFHW9aJyNuLKTKe\naBjFBxGho6OD5uZmxsbGPrG9eyjcifaxw4MMjAhlWQnsKHfxmcJk7LZJoTTSHw6l2vuwjQXoLvgs\nvvIbGUotsrzulJQUPB6PdqRVK47VYfQd4MtADR89phMROWtRVcYRDaP4EgwGaW5upr29nan+Pz40\nFuJ3R4Y4UOWnxR/E67JzSbmL86boRGsf6SO35h68NfdjCw7RXXA2zeW7GE5dZ3nd4x1pc3JydCw8\ntSJYHUaVwLET+xutNhpG8SkQCNDQ0EB/f/+U24MivNY8zP5KPx92juJKMFxQEu5Em+36eCda+3Av\neTV34an9NbbgCF2F5+Irv4HhlLWW161j4amVwuowuofwxHgdVhQXjzSM4ltPTw8NDQ2MjIxMu09V\n5wj7qwZ5uXEIY+D0tcnsrHBTPKkTrWO4h9zDd+E58ptwKK35XDiU3IVLUrvb7cbj8eh0FiouWR1G\nTwHHAi/z8XdG2rRbxY1QKERraystLS0fawo+Was/3In2yZpwJ9pjvOGZaI/P+3gnWsdwN7mH7sR7\nZD9GRulccz6+shsYcS/NwCQ6nYWKR1aH0XlTrdem3SoeTdcUfDL/SIjHawZ5qHqQrqEQayKdaM9a\n5yRxQidax1AXeYf24ak7gJEgnWsvxFd2PSOu/CWp3xhzdEbalJSUJfkMpayyFJ1e20RkKLLsBHJE\npGHRlcYJDaOVp7+/n/r6eoaGhmbcbzQkPN8wxIFKP0d6x0hPinSi3eAibUInWsdQ54RQCtGxdhst\nZdcz4spbsu/gdDrxer1kZWVhs00/K65S0WJ1GL0GnCYiI5HlJOC51dT3SMNoZRqf+ry5uZlgMDjr\nvu+0jbC/ys8bLeFOtGevd7Kj3E1B6kedVxMC7eQd+hU59Q+DCJ3rLsJXdj2jTu+SfQ+Hw0F2djZe\nr5fExOknIFRquVk+AoOIHDdp3VsiYs1sZXFAw2hlGxsbo7Gxkc7Ozjnt39A3xoEqP8/UhTvRnliQ\nxI5yN5tzPupEmxBoI686EkrG0LHuYlpKP8+o07OUX4X09HS8Xi9paWlL+jlKzYXVYfQk8K8i8nBk\neTvwlzpqt1pp/H4/9fX1DA4Ozmn/ngmdaPtHhNLMBHZWfLwTbeJgC3mH9pFd/0g4lNZvp6X0OkaT\nc5byq5CcnIzH4yE7Oxu7feq5npRaalaHURmwD8gGBOgArheR6sUWGi80jFaXjo4OmpqaphzFYSrD\nY8Lv6gIcqPLjGwjicdm4pMzN54qdOBPC73ISB1vIq76DnIZHEWOnPRJKY8nZS/lVsNvtZGdn4/F4\nSE5OXtLPUmqyJZlCwhiTASAiPYuoLS5pGK0+Y2NjR0dxmKvQeCfaKj8fdIzichg+V+LkkjI3OZFO\ntImDPvKr9pLd+BhiHLQX7aSl9FrGkrKW6qsclZaWhsfjISMjY8k/SymwbnK9a4G7ZJqdjDFFQIGI\nvLDAOuOGhtHqNTg4SENDAwMDA/M67lDXKL+p9PPShE60O8rdlGSG+wgl+ZvIq9pLduPjiC2BtuJL\nad1wDWNJmUvxNT4mMTERr9dLdna2jhyulpRVYfSXwA3AK8BBoB1IBkqBs4E+4K9FpNKCmmOahpGa\nakK/uWiLdKJ9ojbA0Jiw1ZPIzgoXx+clYTOGpIFG8qtvJ6vxSUL2RNqLL6Ol5BqCSelL9E0+YrPZ\njo4c7nQ6l/zz1Opj2WM6Y4wDOB84HcgHAsAHwMMiUmtBrXFBw0jB7AOwzsQ/GuKJmgAPVfvpDIRY\nk2pnR7mbs9aHO9EmDdSTX3U7WU1PEbIn01Z8Oa0briaYuPShBDr5n1oaOu24xTSM1ESBQID6+vp5\nP7oDGAsJLzSEZ6Kt7Ql3ot1W6mJbpBNtcn8d+VV7yGz+HSGHMxxKJVcTTFyepto6criykoaRxTSM\n1FS6urpobGyc96M7CHeifbd9hANVgxz0DZNog7OLnGwvd1OY6iC5v5b8yj1k+Z4h6HDTWnIlbSVX\nEUxYniGAxkcO93g8OuyQWjANI4tpGKnpLObR3bjGCZ1oR0NwYn54JtrNOQk4+2spqNpDpu9Zxhxu\n2kquorXkSkLLFEoALpfr6MjhOuyQmg+rGjCcJCKvWlpZnNIwUrNZzKO7cb1DQR49PMijhwbpGxE2\nZDrYUe7m1DXJpPYfJr9qN5ktzzOWkEJrydW0FV9BKGH55jvSYYfUfFkVRrcCnwbeAx4FHhORuXe6\nWEE0jNRcLebR3bjhoPBMXYADlX6aB4LkOCOdaEucZA8epqByNxmtLzCWkEbrhqtpK76ckMNl4beY\nXUZGBh6PR4cdUjOyegSGrcBFwIWEm3Y/RTicXhKR6SeGWUE0jNR8WPHoDsKdaA/6hjlQ5ee99o93\nol03coj8yt1ktL3EWEIaLRv+gPbiywk5lreJtg47pGayZO+MjDFu4DzC4XSKiJywsBLji4aRWggr\nHt2NO9Q1yoEqPy80hqe8OG1NeCbaYzhEQdVu0tteYTQxndYN19BedOmyh5Ldbj/aZ0mHHVLjtAGD\nxTSM1GJ0dnbS2Ng457HuZtI+GOShaj9P1AQIjAlbPAnh/krJhyms2k16+2uMJmbSUnot7et3II7l\nD4bU1FS8Xi/p6enaZ2mV0zCagjGmBPg7IF1ErjLGXAZcAniBH4vIb6c7VsNILVYwGKSpqWleY93N\nxD8a4slIJ9qOQIiCSCfa7Wk1rK/eTVrHQUaTMmkpvS4cSvYkSz53PhITE4/2WdJhh1anqIeRMWYt\nsAfIA0LAT0XkRws81y+A7YRnm906ads24EeAHbhVRL4/h/PdKyJXTVjOBH4gIjdPd4yGkbLK4OAg\n9fX1+P1+S843FhJebBziQJWfw91jpCXZ2LbByTVZNZQd2UNaxxuMJGXTUnYdHeu2I/blbwVnjDn6\nCM/lWt6GFiq6rGpN90PgVyLyygIKyAfyReR1Y0wq4bHtLhOR9yfs4wUCItI/YV2piByadK6zgAFg\nz8QwMsbYgSrCwxU1Aq8C1xEOpu9NKulLItIWOW5yGP0rcIeIvD7d99EwUlZrb2+nqalp1hlm50pE\neL9jlP2Vfl7zDZNgg8+ud3KT9zCfathDatfbjCTn0FL6eTrWXRyVUAJwu914vV4yMzP1Ed4qYOVA\nqdcCWcDkexd5AAAgAElEQVSdwD4ReXeBBf0G+E8ReXzCuquBrwAXi8iQMeYW4HIRuXiK44uAByeF\n0anAP4jIhZHlbwKIyOQgmnyueyOP6QzwfeBxEXlipmM0jNRSmO8Ms3PV1B/pRHskwEgIPp2XyC15\n1ZzedgepXe8ykuzBV3Y9nesuQmzRGfLH4XDg8XjweDw67NAKZnXT7g2EQ+lawAC/Au4UkZo5FlME\nPAtsFZG+Sdu+AZwG3AN8FThfRD7R9GiaMLoK2CYiX44s30C4hd9Xp6kjG/gu4TupWwE/sIvwHdWb\nIvKTKY7ZAewoLS29pbp61cwlqJbZwMAA9fX1BAIBS8/bOxzisUODPHJ4kL7hECXpdv6koJoLevaR\n2v0ew04vLWVfoHPthVELJWPM0T5LqampUalBLZ2lbNr9acJ/yI8VkVk7FRhjUoBngO+KyP3T7HMn\ncDGwYbpOtdOE0dXAhZPC6GQR+dq8vtQc6J2RWmoiQltbGz6fz7JHd+OGg8KzkZlom/qDZCcb/rSw\nissH7iSt9wOGnbn4IqGELXoNDZxO59E+Szrs0Mow1zCa09U2xtiNMRcZY3YDDwE1wDVzOC4BuI/w\n+5jpguhMYCvwAPDtudQzQSOwdsLyGqB5nudQKiYYY8jNzWXLli1kZlo7wV6S3XB+iYsfXpjDN0/P\nIC81gW8fLuOUzr/nxzn/i4Ajg6K3/5WtT+8iu/4RCFkbhnM13i/r7bffpqGhgeHh4ajUoZbfbO+M\nziHcIGAn8Abh90b3T2xwMMOxBtgNdInI/5xmn+OBfYSbV9cCe4EaEfnWFPsW8ck7IwfhBgznAU2E\nH7d9XkTem62++dI7I7Xc+vr6qK+vX7I/yDXdo+yv8vN8wxAgfMXzHjcH7yHLX82QqwBf+Q10FX4O\nbNEdVSEtLe1onyUVf6xqwPAc4fdD9853TDpjzBnAc8A7hJt2A/ytiDw8YZ/TgT4ReSeynADcJCI/\nm3SufYRnls0BWoFvi8jPI9suBn5IuAXdL0Tku/Opc640jFQ0iAgtLS20tLQQCi3NyFsdEzrRDo6F\n+GLGW/wPcx+ewGGG3GvwlX2BrsLzoh5KSUlJR/ss6bBD8cPyd0bGmM8A5SKyJ9IQwC0i9YusM25o\nGKloGh4epr6+nr6+vtl3XqDB0RBP1EY60Q4Gudb1Bn+ecB+5w7UMudfQXL6L7sKzwUQ3CHSq9Phi\ndWu6bxGednyDiJQbYwqBu0TkjMWXGh80jFQs6O7upqGhYVEjgs8mGOlEu7/KT033CJcmHuQbSfdT\nMFpHIGU9vvIb6C44G0z0GxjoVOmxz+owehM4HnhdRI6PrHtbRI5ddKVxQsNIxYpQKERzczNtbW2L\nGhF8NuOdaA9U+TnYHOAS+yv8jfMBCscawqFUsYvu/LNiIpR0qvTYNdcwmmsbzmEREWOMRE6u43ko\nFSU2m401a9aQnZ1t2YjgUzHGsMWTyBZPIs39qRyoOofzjpzM+bzEXw8+QMnBf2IwtRhf+S568s+I\naiiNjo7S3NyMz+fTqdLj1FzvjP4aWAdsA/43cDPhRg0/XNryYofeGalY1dHRQVNTkyUjgs+mbzjE\nY4cHeezQAGeOPs/Xkx5gnTTjT9uAr3wXvXmnQ4w8LtOp0mPDUjRguAi4gPAIDI+JyCOLKzG+aBip\nWLZUwwpNZyQoPFsf4KHKAU7wP8ufJ97PelroTy2ldeMuenNPi5lQcjgc5OTk4PF4dKr0KLCqafdv\nReQCSyuLUxpGKh4s1bBC0wmJ8GbLCA9W9lHa+TR/mvAA600rvalltG36In3eU2ImlADS09Pxer06\nVfoysiqM3hhvsLDaaRipeCEitLa24vP5lqxv0lRqe0Z56MNePM1P8FX7A6yztdORUkH35pvo854c\nU6GkU6UvH6vCqAb4+nTbpxviZyXSMFLxZmRkhPr6enp7e5f1czsHgzxa3Yf7yKP8kXmANaaDFlcF\nvVu/yID3pJgKJZ0qfelZFUadwG8IvyeaTETkSwsvMb5oGKl41dPTQ0NDAyMjI8v6uYHREE/X9mGv\neoQbQ+FQanBuom/rFxnK/XRMhRLoVOlLxaowel1ETrC0sjilYaTi2XL1TZpKMCS82jDA2AcPcvXI\nfRSYLg4nbaZn6xeh4NPLWstc6FTp1tJ3RhbTMFIrQSAQoK6uzrIpz+dDRKhq9xN4Zz+X+O8jz3Tz\nYcIW2jbdhHt97IWSzWYjMzNTp0pfJKvCaOtCZ3ZdaTSM1Epi9ZTn89XSM0j/2w9wbs99eE0P7zi2\nUl+2i9wNJ8TkIzKdKn3hrAqjB0Vk+ywfNOs+K4GGkVppxsbGaGhooKurK2o1+AcD9L11P59pv5cc\n08tBcwxVJTdSvPEEEmyx90c/ISHhaJ8lHXZobqwKox7C04VPuwuwRURK5l9ifNEwUitVf38/9fX1\nDA0NRa2G4EiAgbfu57iWe8milxc5lrfWfIHNW07AnRh7oyeMT5Xu9Xp12KFZWBVGn53DZ42IyIvz\nKS4eaRiplWx83iSfz7fsDRwmMqODjL77ABub7iFd+ngmdBwv5V3PMVs/Ra47NhsTOJ1OvF6vDjs0\nDcuHA1rtNIzUajA0NER9fT39/bNO5rykbGMBHB/cR3Hd3aTIAE8Fj+eJ7Os4ZvOxlGfH5pA+drv9\n6CO8pKSkaJcTMzSMLKZhpFaTrq4uGhoalmXw1ZnYxgZJqbqPgtp7cIcGeDz4aR5wX8OWTVs5qTAJ\ne4w2JtBhhz6iYWQxDSO12gSDQRobG+no6Ih2KdhG/WQdvg/v4Xtwhvw8FjyRvQlXU16xmXOLnCQ7\nYvPxmA47ZP3kel4RaZu0rkJEKhdRY1zRMFKr1XIPvjoT2+gAnsP34am5l6Sgn0eCJ/FTrqJ4w0Yu\nLnWR6YzNP/g2m43s7OxVOeyQ1WFUCfy9iNwdWf5L4GYR2bzoSuOEhpFazaI1+Op07KMDeGvuJefw\nvSQGB3k4eDL/EbyS3LVl7Kxwsz49dptdr7Zhh6wOo3zgp8AQkAt8APyliCzNFJMxSMNIKRgeHqah\noWHZB1+djn2kj9yae/DU3I8tOMQjoc/wb6NXkOItZme5m0/lJsbsH/zVMuzQUkyu9yfAN4EQcJ2I\nPL+4EuOLhpFSH+nu7qahoYHR0dFolwKAfaSX3MP34K29HxMc5lFO4wfDlzOWto6d5W7OWJcck51o\nIdxnaXzk8JU47JDVd0aPAz7gT4E1wC+AZ0Vk2uklVhoNI6U+LhgM0tTURHt7e7RLOcox3EPu4bvx\nHPk1JjjC47bT+V7gMnqTCrm41MX5G1ykxmAn2nErcdghq8PoMhH59YRlB/BNEfnO4sqMHxpGSk3N\n7/dTV1cXEw0cxjmGu8k9fBee2t9gQqM8nXAW/zhwKa22PM4tdrK9zEVeSuw+GltJww5p026LaRgp\nNT0Roa2tjebm5pho4DDOMdRF3uE78RzZD6Exfp98Nt/u3UGd5HJSYRI7y91szInNTrSwMoYdsvrO\nqB8Y3zERSAAGRCR9UVXGEQ0jpWYXrdllZ+MY6iTv0D48dQdAQrziPpf/1buDqtEcKrIT2FHu5uQY\n7kQL4HK58Hg8cTfs0JLeGRljLgNOFpG/XUhx8UjDSKm5i7UGDuMSAu3kHdpHTv1DICHeSP8c/9i3\nk7cHM8l129le5uKcYifOGO1EC+BwOI72WUpMjN27unFL/pjOGPOSiHxmQQfHIQ0jpeYnFhs4jAuH\n0q/IqX8YRHgv+0K+79/B891ppCQYzt/g4uJSF1kx2ol2XDwMO2T1Y7orJizagBOBz4rIqQsvMb5o\nGCm1MLHYwGFcQqCN/Oo7yK5/BIzhkHcb///QDh5rScFm4Ix1yewod1OUEduNCGJ52CGrw+iXExbH\ngCPAzyYPEbSSaRgptXCxNoLDZImDLeRV30FOw6OIsVGXfzH/HbqUX9c7GQoKx3oT2Vnh5rgY7kQL\n4ZHDs7Oz8Xg8MTPskLams5iGkVKLNzw8TH19PX19fdEuZUqJg77wnVLDo4hx0Lx2O7fbLuXu2iS6\nh0KsS3Owo9zFmeucJNhjN5QA0tLSjg47FE1WTa73H3zUiu4TRORPF1Ze/NEwUso6XV1dNDY2xlwD\nh3GJ/ibyq/eS3fg4YhJoWb+D3yRdxp21idT1jpGRbOOiUhcXxngnWoCkpKSjww5F4xGeVWG0a6aD\nRWT3AmqLSxpGSlkrlqaomE6Sv4n8qtvJanyCkD2R9vU7eSr9cu6pTeCNlhGS7IZzipzsKI/tTrQQ\nHjl8fNghp9O5bJ9rVRg5RCS6s2vFCA0jpZbGwMAAdXV1DA0NRbuUaSUNNIRDqempcCgVX8ZBzxXc\nV5vAs/UBgiE4OdKJtiI7IabfK0F45HCPx0NGRsaS12pVGL0uIidEfv4PEfmahTXGFQ0jpZaOiNDS\n0oLP5yOW32Mn9ddTUL2HzKanCdmTaSu+nKrCK9lf7+Cxw4MMjAhlWQnsLHdzSmES9hgdnHXccowc\nblUYvSEix0d+PhpMq5GGkVJLb2hoiPr6evr7+6NdyoyS+4+QX3U7mc2/I+Rw0lZ8BXXrr+SJJgcP\nVg3S4g/iddm5pNzFeUVOnAmx/V5pKUcOX4o7Iw0jDSOllkVnZyeNjY2MjcX2W4LkvloKqnaT6XuW\noMNNa8mV+Iqv5OX2BPZX+vmwcxRXguGCknAn2mxXbPUBmorVI4dbFUaDwCHAABsiPxNZFhE5dtGV\nxgkNI6WW19jYGA0NDXR1dUW7lFk5+w6TX7mbzJbfM5aQQlvJVbQWX8GHfYnsrxrk5cYhjIHT1yaz\ns8JNcYx3ogXrRg63KozWz3SwiNQtoLa4pGGkVHT09fVRX1/P8PBwtEuZlbO3moLK3WS0vsBYQiqt\nJVfRVnwFvpEkHqwa5KnaAENB4RhvIjvL3RyfF9udaAEyMzMpKSlZ8PHa6dViGkZKRU8oFMLn89Ha\n2hrTDRzGuXqqyK/aTUbri4wlpNG64Q9oK76MvlAyj9cM8nD1IF1DIdZEOtGetc5JYox2otUwijEa\nRkpFXyAQoK6uDr/fH+1S5sTVU0l+5W4y2l5iLCGNltJraC+6jGFbMi80DLG/0s+R3jEykmxsi3Si\nTUuKrcYOGkYxRsNIqdgxPpFfMBiMdilz4ur+gILK3aS3v8JoYgatpdfQvn4nQXsy77SNsL/Kzxst\nIyTa4ZwiJ9vL3BSkxkYn2pgJI2PMsSLytjHmGBF5Z8EVxTkNI6Viy8jICA0NDfT09ES7lDlzd71H\nftVu0ttfYzQxk5bSa2lfvwNxJFPfO8qBqsGjnWhPLAh3ot2UE91OtLEURj8Efgz8iYj8zwVXFOc0\njJSKTT09PdTX18fsOHdTcXe+Q0HVbtI6Xmc0KeujULIn0T0U5NFDgzwa6URbmpnAzgoXnylMjkon\n2pgII2PMt4Fs4HrgDqBDRP5pwVXFMQ0jpWJXLE/kN5OUzrcoqNxNauebjCRl01J2HR3rtiP2RIbH\nhKfrAjxY5cc3EMTjsnFJmZvPFS9vJ9qYCKPIiXYC24BHRWT/giuKcxpGSsW+eBjnbiopHW9SUHkb\nqV1vM5KcQ0vZ9XSsvQixJxIS4bXmYfZX+fmgI9yJ9vwSF5csUyfaWAqj74rI3xljviMif7/giuKc\nhpFS8SFexrn7BBFSO9+g4MPbSOl+l5FkL76y6+lctw2xhTudVneNcKBqkBcbPupEu6PcTUnm0nWi\njZkwUmEaRkrFl6GhIerq6hgYGIh2KfMjQmrHQQoqbyOl+32GnbmRO6VtYAu3sGvzj/FQ9SBP1AYY\nGhO2ehLZWeHi+LwkbBY3dtAwijEaRkrFp46ODhobG+OmGfhRIqS1v0Z+5W2k9HzAsDMPX/kX6Fxz\nwdFQ8o+EeLw2wMPVfjoDIdak2tlR7uas9dZ1otUwijEaRkrFr9HRURoaGuju7o52KfMnQlrbyxRU\n7sbdW8mwKx9f+Q10Fp4PtvA7o7GQhDvRVvmp7RkjfUIn2vRFdqLVMIoxGkZKxb/e3l7q6+sZGRmJ\ndinzJ0J664sUVN6Gq+8QQ+5CfGU30FV43tFQEhHebQ+/VzroGybRBmcXOdle7qZwgZ1oYyKMjDF2\n4MvAGsKt6Z6fsO1bIvK/F1xhnNEwUmplCAaDNDc309bWFu1SFkaE9NYXKKjcHQmlNfjKb6Sr8Bww\nH7Wua+gb48EqP8/UBRgNwYn5SeyscLN5np1oYyWMbgVcwCvADcAzIvIXkW2ran4jDSOlVha/309d\nXR2BQCDapSyMhMhoeZ78qt24+moIpKzDV34D3QVnfyyUeoaCPHp4kMcODdI3ImzIdLCj3M2pa5Jx\nzKETbayE0dvjcxYZYxzA/wVygOuAl8ZngV0NNIyUWnnithn4RBIiw/d7Cqp24+yvJZCyHl/FLrrz\nzwLz0fui4aDwzJEAB6r8NA8EyZnQidY1QyfaWAmjD0Vk46R1/wu4EPCKSNmCK4wzGkZKrVzxMt35\njCREpu9Z8it34xyoI5BaTHP5jfTkn/mxUAqJcNA3zP5KP+93jOJyGD5X4uTiMjeeKTrRxkoY7QX2\nisijk9Z/GfgvEYn96QotomGk1MoXt83AJ5Igmc3PkF+1B+dAPYNpJfjKb6Qn74yPhRLAoa5R9lf5\nebExPGLFaWvCM9FumNCJNibCSH1Ew0ip1SGum4FPJEGymp4mv2oPyf5GBtNKaa7YRW/uaTCpAUOb\nP8jDh/w8URMgMCZs8SSwo9zNp/OTyM7Kin4YGWO+ISL/X+Tnq0Xkngnb/llE/nbBFS4jY0wJ8HdA\nuohcFVl3GXAJ4AV+LCK/nekcGkZKrS7xOBr4lEJBspqeJL/6dpL9TfjTy/CV76I399RPhJJ/NMST\nNQEeqvbTEQhRkGrn2uM8/OEFx5GcsLBx8KwKo6Mt5ia3nluu1nTGmF8A24E2Edk6Yf024EeAHbhV\nRL4/h3PdOx5GE9ZlAj8QkZtnOlbDSKnVJ15HA59SKEhW0xMUVO0hadCHP72C5opd9HlP+UQojYWE\nFxvDM9HW9IzxzYs28kef3bCgj51rGM3WC8pM8/NUy0vlNuA/gT1HPzjc/+nHwPlAI/CqMWY/4WD6\n3qTjvyQiM3Uo+FbkXEop9TF2u51169aRlZUVl6OBf4zNTtfaC+kqPI/sxsfJr7qdslf+Fn/GRpor\nbqLPc9LRUHLYDGeuc3LG2mTqh5I494S1S17ebGEk0/w81fKSEJFnjTFFk1afDBwSkRoAY8ydwKUi\n8j3Cd1GzMuFeX98HHhGR162rWCm10qSkpLB582Z8Ph8tLS3x2wwcwOagc91FdK45n5yGR8mr3kvZ\ny3/DQOZmmst30e858WgoGWM4rsBNhitx6cuaZfunjDF9xph+4NjIz+PLxyx5ddMrBBomLDdG1k3J\nGJNtjPkJcLwx5puR1V8DPgdcZYz542mO+0NjzGvGmNdWxG26UmrBjDEUFBSwadMm3G53tMtZPJuD\njvXbee/c26k75s9JDHRQ/vJfU/H8n5Ha/josc+DGRWu6yJ3Rg+PvjIwxVwMXisiXI8s3ACeLyNeW\nqgZ9Z6SUmqitrY2mpiZCoVC0S7GECY6Q0/AIedV3kDjUQX/WMTRXfJGEsnOWpTXd8s1da61GYOJD\nzDVAc5RqUUqtQl6vly1btpCenh7tUiwh9kTaiy7l3XP3Ur/1ayQNNlPx4l+Q/9gt0HhwyT8/XsPo\nVaDMGFNsjEkErgVW7ZToSqnoSExMpLS0lOLiYhyOhY2KHWvEnkh78eW8e+5eGrb8CQl9R6C3fsk/\nN+bDyBizD3gRqDDGNBpjbhaRMeCrwGPAB8DdIvJeNOtUSq1eWVlZbNmyhezs7GiXYhmxJ9FWciUN\nl++HTZcu+efFfJSLyHXTrH8YeHiZy1FKqSk5HA6KioqONgOPyzmTpiAOJ9iW/r4l5u+MlFIqnqSl\npbFlyxa8Xm+0S4krGkZKKWUxm83G2rVr2bhxI06nM9rlxAUNI6WUWiJut5tNmzZRUFAwr9lVVyMN\nI6WUWkLGGPLz89m8eTMpKSnRLidmaRjNwhizwxjz097e3miXopSKY8nJyVRUVLB27Vpsy9AgIN7o\nv5FZiMgBEfnDldKxTSkVXSuts6xVNIyUUmqZrcTOsoulYaSUUlEy3lk2Kysr2qVEnYaRUkpFkcPh\noLi4mNLSUhITl36qhlilYaSUUjEgPT2dzZs34/F4ol1KVGgYKaVUjBifWbaiooLk5ORol7OsNIyU\nUirGjM8sm5eXt2o6y2oYKaVUDDLGUFhYyMaNG3G5XNEuZ8lpGCmlVAxzuVxs3LiRwsLCFd1ZduV+\nM6WUWiGMMeTl5bFp06YVO6SQhpFSSsWJ8SGF1q1bh91uj3Y5ltIwmoWOTaeUijUej4fNmzevqCGF\nNIxmoWPTKaVi0UobUkjDSCml4thKGVJIw0gppeLcxCGFEhISol3OgmgYKaXUCpGens6WLVvIycmJ\ndinzpmGklFIriN1uZ/369ZSXl5OUlBTtcuZMw0gppVag1NRUNm/eTG5ubrRLmRMNI6WUWqFsNhtr\n1qxh48aNMT/wqoaRUkqtcG63m82bN5Ofnx+zA69qGCml1CpgjKGgoIBNmzbF5MCrGkZKKbWKOJ3O\nowOvxtJdkoaRUkqtMuMDr27evDlmBl7VMFJKqVVqfODVtWvXRn16Cg2jWehAqUqplc7r9bJ582ZS\nU1OjVoOG0Sx0oFSl1GqQlJREeXk569evj8r0FBpGSimljsrJyWHLli3LPj1F/I87rpRSylIJCQmU\nlpbS1dXFwMDAsnymhpFSSqkpZWVlLdvUFPqYTimlVNRpGCmllIo6DSOllFJRp2GklFIq6jSMlFJK\nRZ2GkVJKqajTMFJKKRV1GkZKKaWiTsNIKaVU1BkRiXYNccEY0wtUT1qdDkw1nPfk9TlAxxKVNpvp\nalzq88x1/9n2m2n7XP/9T7cuWtclWtdkPscs9Losdr3+rix8v1j9XVkvIp5Z9xIR/WcO/wA/ncu6\nqdYDr8VS3ctxnrnuP9t+M22f67//GdZF5bpE65osx3VZ7Hr9XbH+msz3ukTrd0Uf083dgTmum2l9\nNFhVy3zPM9f9Z9tvpu3z+fev12R+xyz0uli1Phr0d2Vun7Mk9DHdMjDGvCYiJ0a7DvVxel1ij16T\n2LQc10XvjJbHT6NdgJqSXpfYo9ckNi35ddE7I6WUUlGnd0ZKKaWiTsNIKaVU1GkYKaWUijoNoygz\nxpQYY35ujLk32rWsZsYYtzFmtzHmZ8aY66NdjwrT34/YY4y5LPJ78htjzAVWnVfDaBGMMb8wxrQZ\nY96dtH6bMabSGHPIGPM3M51DRGpE5OalrXR1muf1uQK4V0RuAXYue7GryHyui/5+LI95XpNfR35P\nbgKusaoGDaPFuQ3YNnGFMcYO/Bi4CNgMXGeM2WyMOcYY8+Ckf7zLX/KqchtzvD7AGqAhsltwGWtc\njW5j7tdFLY/bmP81+VZkuyUcVp1oNRKRZ40xRZNWnwwcEpEaAGPMncClIvI9YPvyVri6zef6AI2E\nA+lN9D/SltQ8r8v7y1vd6jSfa2KM+QD4PvCIiLxuVQ36S2e9Qj76L2wI/5ErnG5nY0y2MeYnwPHG\nmG8udXFq2utzP3ClMea/iK0halaLKa+L/n5E1XS/K18DPgdcZYz5Y6s+TO+MrGemWDdtz2IR6QQs\nu6BqVlNeHxHxA19c7mLUUdNdF/39iJ7prsm/A/9u9YfpnZH1GoG1E5bXAM1RqkV9kl6f2KTXJfYs\n6zXRMLLeq0CZMabYGJMIXAvsj3JN6iN6fWKTXpfYs6zXRMNoEYwx+4AXgQpjTKMx5mYRGQO+CjwG\nfADcLSLvRbPO1UqvT2zS6xJ7YuGa6ECpSimlok7vjJRSSkWdhpFSSqmo0zBSSikVdRpGSimlok7D\nSCmlVNRpGCmllIo6DSOlLGCMCRpj3pzwz4xThywXY8wRY8w7xpgTjTEPRGo7ZIzpnVDradMc+2Vj\nzO2T1uVGphpIMMbcZYzpMsZctjzfRq1k2s9IKQsYYwZEJMXiczoiHQ8Xc44jwIki0jFh3dnA10Vk\nxlHkjTGZQDWwRkSGIuu+ChwjIn8UWd5LeB6oXy+mTqX0zkipJRS5M/lHY8zrkTuUjZH17siEZq8a\nY94wxlwaWX+TMeYeY8wB4LfGGJsx5v8aY96LzIH1sDHmKmPMecaYByZ8zvnGmPsXUedJxphnjDEH\njTGPGGNyRaQbeAG4ZMKu1wL7Fvo5Sk1Hw0gpazgnPaabOANmh4icAPwX8PXIur8DnhKRk4BzgH8x\nxrgj204FdonIuYRnoC0CjgG+HNkG8BSwyRjjiSx/EfjlQgo3xiQBPwKuFJFPA3uB70Q27yMcQBhj\n1kZqeXYhn6PUTHQKCaWsERCR46bZNn7HcpBwuABcAOw0xoyHUzKwLvLz4yLSFfn5DOAeEQkBLcaY\npyE8jn/kfc4XjDG/JBxSNy6w9k3AFuAJYwyAnfCIzRAeGPPfjTEphKeYvjtSi1KW0jBSaukNR/43\nyEe/c4bwnUjlxB2NMacA/omrZjjvLwlPBDhEOLAW+n7JAG+LyJmTN4iI3xjzBOFZV68FvrLAz1Bq\nRvqYTqnoeAz4moncihhjjp9mv98TnoHWZozJBc4e3yAizYTnl/kWcNsianmf8KyqJ0dqSTTGbJmw\nfR/wV0CGiLy6iM9RaloaRkpZY/I7o+/Psv93gATgbWPMu3z0jmay+wg/MnsX+G/gZaB3wvY7gAYR\neX+hhYvIMHAV8G/GmLeAN4BTJuzyKOFHiHcu9DOUmo027VYqxhljUkRkwBiTDbwCnC4iLZFt/wm8\nISI/n+bYI0xq2m1xbdq0W1lC74yUin0PGmPeBJ4DvjMhiA4CxxJu/TadduBJY8yJVhdljLkLOJ3w\nOyulFkXvjJRSSkWd3hkppZSKOg0jpZRSUadhpJRSKuo0jJRSSkWdhpFSSqmo0zBSSikVdf8PMCk1\nwfkFsT0AAAAASURBVDWvBh4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "energy_range = [0.1, 100] * u.TeV\n", "model.plot(energy_range=energy_range, energy_power=2)\n", "result.model.plot(energy_range=energy_range, energy_power=2)\n", "result.model.plot_error(energy_range=energy_range, energy_power=2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Exercises\n", "\n", "* Change the observation time to something longer or shorter. Does the observation and spectrum results change as you expected?\n", "* Change the spectral model, e.g. add a cutoff at 5 TeV, or put a steep-spectrum source with spectral index of 4.0" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": true }, "outputs": [], "source": [ "# Start the exercises here!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## What next?\n", "\n", "In this tutorial we simulated and analysed the spectrum of source using CTA prod 2 IRFs.\n", "\n", "If you'd like to go further, please see the other tutorial notebooks." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "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.6.1" }, "nbsphinx": { "orphan": true } }, "nbformat": 4, "nbformat_minor": 1 }