diff --git a/SiPM_peaks.ipynb b/SiPM_peaks.ipynb index 7209677bb8435ffea826ae74a938a34941330823..b3dbcb01cdb33c89ce79fe1cc262ff1947a58182 100644 --- a/SiPM_peaks.ipynb +++ b/SiPM_peaks.ipynb @@ -214,90 +214,6 @@ "plt.show()" ] }, - { - "cell_type": "code", - "execution_count": 8, - "id": "b6692d1c", - "metadata": {}, - "outputs": [], - "source": [ - "from scipy.optimize import curve_fit\n", - "\n", - "# Define the Gaussian function\n", - "def Gauss(x, a, mu, sig, c):\n", - " y = a*np.exp(-(x-mu)**2/(2*sig**2)) +c\n", - " return y\n", - "\n", - "parameters=()\n", - "def fit(i,j):\n", - " xdata = np.asarray(ds[\"ADC Channel - Charge distribution\"][i:j])\n", - " ydata = np.asarray(ds[\"Counts - Charge distribution\"][i:j])\n", - " param, covariance = curve_fit(Gauss, xdata, ydata, p0=[1000, i*2, 200, 0])\n", - " fit_A=param[0]\n", - " fit_B=param[1]\n", - " fit_C=param[2]\n", - " fit_D=param[3]\n", - " print(\"mu=\", fit_B, \"\\nsigma=\", fit_C, \"\\n\")\n", - " fit_y = Gauss(xdata, fit_A, fit_B, fit_C, fit_D)\n", - " plt.plot(xdata, ydata, c=\"darkblue\", label='data')\n", - " plt.plot(xdata, fit_y, '-r', label='fit')\n", - " plt.xlabel(\"Canal ADC\")\n", - " plt.ylabel(\"Cuentas\")\n", - " plt.grid()" - ] - }, - { - "cell_type": "code", - "execution_count": 9, - "id": "f15890ec", - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "mu= 19.157181333696382 \n", - "sigma= 49.24544222362323 \n", - "\n", - "mu= 479.0099909259773 \n", - "sigma= 51.310704584285276 \n", - "\n", - "mu= 939.6075158534509 \n", - "sigma= 52.7541516809037 \n", - "\n", - "mu= 1396.856500246315 \n", - "sigma= 58.13163757023065 \n", - "\n", - "mu= 1851.4849042741075 \n", - "sigma= 66.47690197710251 \n", - "\n", - "mu= 2292.507518508694 \n", - "sigma= 90.54113833966873 \n", - "\n" - ] - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa4AAAEYCAYAAAAEZhLyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABeAklEQVR4nO2dd3hUVfrHP2dmkkx6J0AghCJFEAVBBUWxCzbsrq4N67rqYt1Vf5ZV13VX3dW166qo2FZcdG3YQEAUkab03hNI723K+f1x74T0zCQzmZnk/TzPfe7Mveec+85Mcr/3nPOe91VaawRBEAQhXLAE2wBBEARB8AURLkEQBCGsEOESBEEQwgoRLkEQBCGsEOESBEEQwgpbsA3oCaSlpens7OxgmyEIghA2LF++vEBrnd7SORGuLiA7O5tly5YF2wxBEISwQSm1s7VzMlQoCIIghBUiXIIgCEJYIcIlCIIghBUiXIIgCEJYIcIlCIIghBUiXIIgCEJYIcIlCIIghBUiXEK3o7i4hvfe2xBsMwRBCBAiXEK348ILP+E3v/mUXbvKgm2KX1i2bB8ffbQ52GYIQsggkTOEbsf69YUAuN3dI0nq+PGzAND6jiBbIgihgfS4hG5HdbUTgLo6V5At6Tz791cG2wRBCDlEuIRuhcvlpqioBoAlS3KDbE3n+eSTrcE2QRBCDhEuoVuxdOm++tdXXPFFEC3pPGvXFnDttV8F2wxBCDlEuIRuRU5ORaP3t946n8rKuiBZ0zl+/TU/2CYIQkgiwiV0K7Ru7JDx1FPLeeKJ8Ewp43C4g22CIIQkIlxCt8fpDE8BCPehTkEIFCJcQrfC5eoeLvAACjdf8jIb+DspiHehIHiQdVxCt8LpdBOJky95BQdWpnI1SgXbKt/Zu7ecKWzkFIyFx1ezFLdbY7GE4YcRBD8jPS6hW+F0upnCBiazjZPZzEmEZ8SJfv1eYirrKSeKlfTlDNY3czwRhJ6KCJfQrXA6NVPZQDU2arEylQ2ocOxyAceyne/JZkPmGI5gNyce82awTRKEkECES+hWOJ1uDiOH7xnIT2RxOHuCbVKHsOFiGPmsoi9RE4/AjpPInduCbZYghAQiXEK3YvvWIg5mP2vJYBV9GU0uSoeXV6HbrRlCAZG4WEsG+5OzABhOXpAtE4TQQIRL6Fa8//jXxFHHIRcdx3p6EUcdrj05wTbLJ6qrHYxkPwDryKCy30DcKEaIcAkCIMIldDM8vZIJV5/EdlIAmP/6wmCa5DPV1c76z3Hvm1eQ1j+VHSSLcAmCibjDC90Gt1vTnxIAVPaAeuHKpiiIVvlOSUktWZSwnzhGH5XNoEGJfHtVKiPtpcE2TRBCAulxCd2G2lon/SjFhcLWry87SQZgYJgJ12mnfUh/SthNEjabwmq1EDF4IL0d4fU5BCFQiHAJ3YbaWhf9KGUf8djskdQSQQ4JZFMcbNN8YuvWElO4ErHZjH/RiuTe9HKVQm1tkK0ThOAjwiV0G/btqySTUvaQiFKK6Ggbu0msHz4MJ7Lqe1zGv2hNeh/jxN69QbRKEEIDES6h23DzzfPoRynWrP4ALFx4MfuIJ4PwijiRQDUJ1LKrgXDV9c409lu2B9M0QQgJRLiEbkNlpYN+lFIYnQrA6NHp7CeeDMqDbJlvZJk9xIY9LndmPwBqNotwCYIIl9BtmHZiLxKpYcJ54wCIjLSSODSLdCrB5Qqydd7Tv4FwWa1GuCprP6PH5dgdXmvSBCEQiHAJ3YaoQmOdU/zwgfXHdK8MrGgoKAiWWT5z5rgYAHJIqD8W0yuJSiJw7RXhEgQRLqHbYC0qBEBlZNQfq04w1nKxf38wTOoQiTXGeq1r7j6R+PhIABIS7ewnHr1PFiELggiX0G2IKDF7Venp9cdqktKMF2EkXHlrtlNOFP/36In1ke0TEiLZRzyWvPD5HIIQKIIuXEqpVKXUNUqpOUqpLUqpaqVUqVLqe6XU1UqpFm1USk1USn2ulCoy6/yqlJqhlLK2ca0zlFLfme1XKKV+Ukpd0Y59VyillprlS836Z3T2cwv+J7KsxHiRllZ/rNYULmeYxCt0uzW9qGA/cY2OJyREsp84bEXS4xKEoAsXcAHwCnAk8BPwFPAhMAr4N/Af1SShklLqbGAhcCwwB3gWiAT+CbzX0kWUUjcBn5jtzjKv2ReYqZR6opU6TwAzgT5m+VnAIcAnZntCCGEvNyNLNBCuxVsMp4y5bywJhkk+U1XloBcV5DURrgEDEii0JWIryA+SZYIQOoRCrMJNwFnAZ1ofyD+hlLoHWAqcB5yLIWYopRIwRMQFTNZaLzOP3wfMA85XSl2stX6vQVvZwBNAETBOa73DPP4Q8DNwu1LqQ631jw3qTARuB7YC47XWxebxx4HlwBNKqU89bQnBJ7qyhCpLFDHR0fXHdhZpqrFRszM8Fu56hGsbqY2OR0XZqE1KI7agDJxOsIXCv64gBIeg97i01vO01p80FC3z+D7gRfPt5AanzgfSgfc8omWWrwH+z3z7uyaXmQ5EAc82FBpTjB41397QpI7n/V88omXW2QE8Z7Z3VfufUOgqYqpKKY1IaHRMWRQFxBJfFx6LkCsrW+5xAVQnpmFBQ770uoSeTdCFqx0c5t7Z4NgJ5n5uC+UXAlXARKVUlJd1vmhSpjN1hCASX1NKeVRj4bJYFIXEEOcIj0XIVRV1pFPZbI4LwJFiDIHuWrqpq80ShJAiZIVLKWUDLjffNhSPYea+2X+v1toJbMcYAh3kZZ1coBLop5SKMa8dC2QCFeb5pmw290O9+jBClxBXU0ZlTGKjYw8/fAxFxNA7MjyC09bm5mFFt9jj2l5tDIE+cddnXW2WIIQUIStcwGMYjhSfa62/bHDcc2dqLTmR53hSB+okNtn7co1GKKWuU0otU0oty5ehnYBTV+civra0kSs8wEknDaBExRJbGx5Dhd+9vxyAS2Yc1+zc1mLj3zUhTHqPghAoQlK4lFK3YDhGbAAuC7I5HUJr/bLWepzWelx6k5up4H+2bCkmjUrsmRnNzhXoaNwFhUGwync+ec3wD4of3K/Zuf1OI6JGiruyS20ShFAj5ITLdDN/GlgHHK+1bpo9r2nvqCme4yUdqFPaZO/LNYQgkreriDjqiGpBuIqIIYUqtNvdQs3QopcZyd7SuwUBdkXhRpHkCo/eoyAEipASLqXUDOAZYA2GaO1rodhGc99sfsmcFxuI4cyxzcs6fYBYYI/WugpAa10J7AXizPNNOcjcyyx5iKCLS4x9cnKzc0VEY8eJqzz0eyoe4bL27tXs3GmnD6GYaJKlxyX0cEJGuJRSf8RYQLwKQ7RaCxEwz9yf1sK5Y4EY4AetdcPZ+LbqTGlSpjN1hGBRUmLsE5t3koswhthc+aEfaDeFKgAiMpoPL7/00smUWONIcsocl9CzCQnhMhcPP4axsPdErXVbd5jZQAFwsVJqXIM27MAj5tsXmtR5HagFbjIXI3vqJAP3mG9fbFLH8/5es5ynTjbwe7O919v7bEIXUWqO7iYmNTvlEa5wmOdKoZpioomOi2p2LirKRk1sInF1IlxCzyboy+/NWIEPYUTCWATc0iTCE8AOrfVMAK11mVLqWgwB+04p9R5GRIyzMNzeZwPvN6ystd6ulLoT+BewTCn1PlCHsZi5H/Bkw6gZZp0flFL/AG4DflVKzcYIK3URkALcLFEzQgdLmSFclpSkZuc8wqULm06Xhh4pVFFENCn2lv81K6ISSKoM/Z6jIASSoAsXxpwUgBWY0UqZBRgxAwHQWn+klDoOuBcjJJQd2IIhMv/SWuumDWitn1FK7QDuwFgfZsFwAPk/rfUbLV1Ua327Umo1Rg/rOsANrAAe11p/6tOnFALKAeFqeY4LwqPHlUoVRcTQ195yrOiK6ESySiQLstCzCbpwaa0fBB7sQL3FwFQf63yCEWjXlzozaSCaQmhiKS8z9slJzc6dfdWR8DroonDpccVgb6XHVR2dQKJTvAqFnk1IzHEJQmexVRg9Lmtq8x5Xv9FZxoswGSosJIYWhssBqI5NJEbXQXV1F1smCKGDCJfQLVClpTiwYEtoHiqJ6BiqsaGKQn+o0NPjao2a2CTjRWHofxZBCBQiXEK34OdvN1KKncio5kNsNpvFEIPi0O5xuZ0ukqluU7hq45KMFyJcQg9GhEvoFiRRQwnR2GzN/6StVkUx0bgKQlu4XEXFWNFtCpcz0RwKLRDPQqHnIsIldAuSqKaEaOLiIpudKy+voxQ7vyzYEgTLvMfj9eiMT2q1jCMxBQAtwiX0YES4hG5BirWW+P7NwyQBVFc7KSGaGEdoh0py5RvCdfx5h7ZeJskQLtd+yTgg9FxEuISwp7raQbyrCp2Q0OL5ujoXpdhJpKaLLfMNzwJpR0JSq2XcScZQoStPelxCz0WESwh78vKqSKIG1UKAXTggXEmEtgv5wo9+AaDO44DRAjtzaygjii1LNrdaRhC6OyJcQthjCFc1thbWcAE4HG5KiDZ6XM2DqoQMn721DICamJZ7jgA//phDITHkrtvdVWYJQsghwiWEPXk5ZUYurl4pLZ739LgicUFN6A4XeiLDO+NbSwMHL7xwEoXEkhXr6CqzBCHkEOESwp6a/cbcUER6y8I1fnwfSsx4hfXpT0KQ3pG1lGLnuBMHtlpm7NgMSlQMUVVlXWiZIIQWIlxC2KNKS4CWA+wCXHzxcOIz04w3nvQnIciQZDeVUfGMHJnWahmlFBVR8dgrQvdzNKW0tJbHH1+K2x26w7RCeCHCJYQ92hQjawspTTwkeFzlQ7jHFV1dRlVUfLvlqqITiK4Jnx7XrbfO5667FvL22+uCbYrQTRDhEsIeSzs9LoAauxnDMIR7XLG15VTHtD6/5aEuNpHYugpwu7vAqs5TXGzMK15++RdBtkToLohwCWGP1UxpYktreY4LwBFj9mRCuMeV4KigNq594aqNTcSKDmkRbkjDSPclJaHrHCOEDyJcQthjqTBS2Uektd7jqnfcCNGbvdaaRHdlm+GePFRGm+7yYZBfDODLLw8kvpw06b0gWiJ0F0S4hLAnwnRUiEhPbbVMcnYGAHUhGnGiqqKWFKpwt7KIuiHzVhmJJH/9bmOgzfILVVXO+tdr1oTm9y+EFyJcQtgTUVmOG4UlqfVhtvjeyTixUJtX3IWWeU/pngKsaFRK68OdHgrN6PG7Vm5vp2Ro4nKFx9ycELqIcAlhj62qnHKiwNL6n3N0TASl2HEXl3SdYT5QsWs/ANY2eo0ePGlPoqtD37Pw559zicTJbSzgdAyvwosu+iTIVgnhTvOse4IQZkRUlFFmiaYttwa73UopdlJCVLiq9xjCZctIb7esR7jslaE5X9eQI454m2f5H7/nRwCO4wY+/BDcbo3FotqpLQgtIz0uIeyxVZZRaYtts0x0dIQRPaMsNG/2tblGmpKoPu0LV7EZBcQeBtEzUqjkOn7i17FTySWBO1kAGBH9BaGj+EW4lFLDlVK3KqWuV0q1788rCH4ksqqC6qi2hcvT41Ih6lXo2G84LURntpxTrCEurJRgx14Vmp+lIWexjgjcrDz6QvRvfsMpbCKBaiorRbiEjuOTcCml7ldK5SqlUhocOwlYCTwBPA+sUEq1P1AvCH5Aa40qK8ER23pEdQC73UYpdizlodlL8SSRjMvK8Kp8ETEhP1SYl1fJGaxjJ0lsihuI88yzicTFqWwS4RI6ha89rinABq11wwUkfwU08ADwAjAQ+IN/zBOEtnE43MS7qrBntP2sFB8fSQl2XIUhuvapyBCuhAG9vSpeSAxRId7jOurIt5nIThYyiNS0GGwTj6KSCI5mRyMXeUHwFV+FKxtY73mjlMoEDgee11o/orW+CZgHTPOXgYLQFtXVThKpwdVGKhCA0aPTKSWaqJrKLrLMN1RxMeVEYU+I8ap8ETHsX78nwFZ1DueOnfShnL39R3HLLWOJT4llKVlMZIf0uIRO4atwJQMNH1mPxuhtfdrg2HIgq5N2CYJXVFXUkkgN7vi2g9NaLApHTDxRdVUhGePPVlZMicU70QKjx5XgqAigRZ1j//5KJrATgH7nn4zNZiE2NoIlZHEYOVSXhK7tQujjq3DlA5kN3h8POICfGhyL7EC7gtAhagpKsKLRiUntlq20xWBBQ1nozXNZS0uoiGx7ns7DrbceThEx9YknQ5GrrprLOPZQg42izIMA4+HhF/oSgRu1MTyifgihia8Cswo4Syk1Sik1BLgI+F5rXd2gTDaQ6x/zBKFt6vKNSBg6KandsjvLzD/3EPQsjKoopaYdBxMPM2YYwpVMdUj2HgGKimo4mP1sIB0iIuqPP/LRDQBEblrfWlVBaBdfhevvQCLwC7DRfP2k56RSyooxfLjMXwYKQls4TG+8tsI9eSjFDoAuDr2wT3GOcmpivBMui0VRiNF7DMXPAkZYpxHksZ6MRgFNrCOGUYcV67q1wTNOCHt8Ei6t9SLgDOAjYA5wvta6YZKdicBe85wgBJyty415lMh2vAoBkgcYa6QcBaF3s09wVlIT7Z1wORyu+ugZeRtC00HDVldDNsWspxdW64HbTGJaHBtIp3rZquAZJ4Q9Ps9Faa3naq3P01qfr7We0+TcIq31GK31bP+ZKAit8+8nFwEQ3bv1dPceTjh3NAB1eYUBtclntCbJXUltXJJXxdPTY+oD7ZbvCM1R+X6VOVjQrCOjUWinlJRo1tCb/mW7WLQoNEVXCH3EiUIIaxIxplcj09uPqq7MeTBHfoit5SovJwK3V0kkAeLiInnomTMBqN6bF0jLOsyAqr0ArKcXl1wyotG5zaSRRQkfvb8mGKYJ3YAOB9lVSvXD8DCMaum81nphR9sWBG8ZkgoUQt+D+7Vb1pJi5LpyFpUE1igfyV27iz5AvtPudR1PaKi6ffkBsqpzZFftwYmFMRccTXx8ZKNzW0jDiqZ3TWjaLoQ+PguXUuoU4J/A8HaKWjtkkSD4wJA0BYVgS20/AaM1NQkAV1FozXF98NISbgE++aGYa7ysY8swhkYtxSHWezRJK81hB8n0G9R87nELxrE+laE5zCmEPr7GKjwKY7FxEvAsoICFwCvABvP9J8BDfrVSEFrBXldJjSUSIiPbLRudEEMVEegQE65kt7EY93f3nOB1HU+2Z0tJaH0WD9kUs4NkHnro6GbntmCIbmrx3q42S+gm+DrHdTdQA4zXWnviEc7XWt8AjAIeAU4CxDlD6BJi6ioot3oXcaI+tUlJaK3jiqgw7Dn1ksO9rmOPjaKYaBZ89GugzOowBQVVDKCYnSQTGdl84OWBf51JKXb6Vu0PgnVCd8BX4ZoA/E9rndO0DW1wP0Yswz/7yT5BaJPYukoqI9pOaeIhOtpGCXZUaUlgjfIRW5nRa1Kp3idVsNutFBJDaghGz3j+Hz/Sh3J20PLw7YUXDWcLqSQV7u5iy4Tugq/ClQjsavC+Dmh611gMHNsZowTBW6IdVVTZvOtxxcTYKCEaFWKpTSLLzR5gcvvzdB7sdlvIhn1KrzHm3Xa2IlyRkVa2kEZCngiX0DF8Fa48aPTXmAcMblImAswUrYIQYOIcvvW4SrFjKw+tocLIylIqVRREteig2yJRUUaPKxSFK6XMcLpoTbgiIixsI4W4gr246iRKvOA7vgrXJhoL1RLgZKXUUAClVG/gPGCzf8wThLaJdVZSFemdcMXFRVJCNNaK0Opx2StLKbV69xk8WK0WikJ0qDCjxsjmfP+rv2nxfGSklZ0kE4GbRf/5uStNE7oJvgrXXOC4BhmQn8boXa1USv2M4VmYDjzlNwsFoQ3iXNVURcZ5VbZXrxhKiSaiMrSEK7q6jLII7z5DQ5IH9yHdUhMAizqHff9enFgYd3bLziY2m6W+N5ZQLC7xgu/4KlwvYcxfOQC01ouBC4DtGF6FucDvtNZv+tNIQWgRrYl3VVET5V1vxWaz4IxLMJJJah1g47wnuqqMKrt3cQobUh2dSLy7ClyuAFjVcWLyc9hLAjEJLS+oVkrVC1fUPnGJF3zH1yC7ZVrrn7TW5Q2OzdFaj9JaR2utR2itX/a/mYLQAjU1RGon1XbveysqORGb2wnV1e0X7gLcbo2tvIQS5X0SSQ/lEXEhGSE+viiHXSqFiIjWYxC89s11AFh272q1jCC0hq8LkI9VSrWZ3Vgp1V8pJV6FQuAx82rV2tvOftwQqxn2iZKSABjkO1u2FJNCFZuLVPuFm/DNSuP58bv//uJvszpFYsk+9tradu1P6JNMATHsWSzxCgXf8XWocD5wZTtlLjfLCUJAqTHj9OkE74fZXPFmINsQEa7aGiepVDHmhPYiqDWn0FyJUrgpdKKs5+WUklBewGZH279JenoMO0nGuW1H1xgmdCt8FS5vHgsVEDoTCILXfPzxFoYPfw2nMzSz6jbl3Rd/AuDrn713b68XuRDJglxXVEoEbuIG9PG5ricnl70qdJxNdv2wHhvuVl3hPXiEawChNcwphAeBSGsyAChvt5QQclx77Zds3FhEYWFozP+0h73WiPGXPqS313VWbq8D4JcFmwJik6+4zAzOpLSflqUpUy41vPYiK0NDhAFqNm4FWl/D1RCPcJUUh8ffmxA6tBsdXil1f5NDk5VqseNlBbKAi4HvO2+a0NV4Ev65XOHRYU631gJw56OneF3nh/XGuqfNP2/n0IBY5Rse4VKpvgvXbY+cAm+DrTR0ei2v//lzjgFm/PPcdsvuJJlYHGzfkkPS+KZxDAShdbxJa/Jgg9camGxurbEX+FOHLRKChifFel1daLlXt4bHmy5xQIbXdX5zw1HwIsS7Q+MpXxcawmVN9z5OoYfI9BTcqAMho0KAvg5j8fGp17Tvn7WTJAAq1m4GES7BB7wZKjze3E7AmL+a2eBYw+1YYCQwQGu9PBDGCoElJ8cYenv66RUsX74vyNa0j6XMuGHHmEkVveGev50KgC1Uomd0RriiIykmmsiKEj8b1XGyKSaXeGxx7bv33/fKxQB8/sJ3AbZK6G602+PSWi/wvFZKvQF81PCY0D2YM+dAlK6nnlrOU08tR+s7gmhR+1jKSqnBht2Lm6QHW3wsdVhDJuyTzjM8I2OyfHfOsFgURcQQVRE6Pa6RsZWUR/bBm0+TMmYYAPuWrgusUUK3w9cFyFdprf8XKGOE4HHuuR8H2wSfiagoo1T5GM9ZKUpVNBEhIlyOHCMnVXy278IFUGKJwV4VOsLVp66Agjjvhm5TD8qknCjxLBR8xps5rhZRSsVgRIpvcXm81lqWxIcRUVFWamtd9KIcC5p9+B6CqKuJrCyl3BqL9zNcBmWWGCKqKwJik6+sXbCe0USTnOR75AyAEkscGaHiDu9208dRxK9x3g3dxidEsYYkDomvDLBhQnfDZ3d4pdRlSqk1GC7vuzDiFDbdtvnTSCHwpKZGM4w8tvIYa3mCTEpCfj1XVU4BxbQcD68tKq3RRFUHf8WG261Jp5J84mjFU7ddSq2xRNeEhnCt+WYNkTj5aV+E13Wq0zPJdBYF0CqhO+JTj0spdSXwGuACFgG7Aaf/zRK6mrKyWu5jEXEY65yuZwkzZszj2WdPCrJlrRPjqKQA33sqFRGxpNYEv8dVV+cinUpiOrD42EN5VDwxISJcqz7+mVHAymLvf5OiuF4cVLQxcEYJ3RJfhwrvAIqBY7TW6wNgjxAEXC43FRV1TGEDszmEZKqZxlpO/9/WkBauZKop69Xf53pVEbH0r80JgEW+YQhXBTXx/TrcRm1MArFlleB0gq3DI/9+4dPn5/NbYIcXi489FCdkkOSqhIoKiPM9tYvQM/F1qHAI8IGIVveiosJBNsUMoIT5DGYBgxjJfmIcoZek0IPWmiSqSczy3hXeQ3VUHNEh8Nlqa40eV12C74uPPTjik4wXIRB70eNk4U3UDA+lSWbUk507A2GS0E3xVbiKgNpAGCIEj7KyWo5mOwCuCcewhCwsaEZWhe5UpaPORRLV1MX67kRSY48lzhl8h4DaagepVFGX2HHhKjUTUNaY3onBJJsiComhkiiv69RkZBovRLgEH/BVuD7FCPnUsZlkISSZN28Xo9iPyxrBFY+dx1KMzDUHV20PsmWtU1dcRgRunHGJvteNTsDuroPa4D6DOfIKseHGkeT74mMPw44aAkDptuBnEh5ACTtJxmbz/rbi7m/8rbm2he7fmhB6+CpcdwNRwItKKRmQ7iZceeVcRrKP6v6DmHBsNqVEs5cEBjtDN3qGw4zx54zzvceVU2XMBdXlB9ebrWyrkY7ElZLW4Tbs/Y2htpce+cYvNnWGARSzg+T6mJfeYOufSR1WajeFbu9eCD18Fa4PgCrgGiBXKbVCKTWvhe1b/5sqBJKR7Mcx7EBOqA30Yjj5QbSobVwFhui4Enzvcf2yy/CcXPdjcG+WN138HgDu1I4L1z5tPD/uWb7FLzZ1FJfTRTZF7CQZl8v7ZRTpGXHsJpF5MyUut+A9vrohTW7wOhY4rJVy4RFeXAAgEifZFOMYOxqA3/3uUDa8kM6lrAStIQRHhp2mcLkTknyumzWqP6wBV2Fwe1zpGPNsOq3jwnXO9UfDQzBhYItxALqM2pw8YnGwk2Ti4iK9rpeVFc9OkkkpD/4cnRA++BryyeLlFtz/IsEnThmmsKCJGm7Ml8TFRbKRXiRRA/mh2etyF5UY+w70uKbfPgmAij15/jTJZ9Ix1pJFdyBOoYe0vkmUWGNJdwd3LZdji9F73W9P44cffuN1vYkTM8mJSGWQNXTCVgmhTyASSQphRnKpOZc1YAAA9903ge2Ynm47dgTHqHZwm1HVdVKSz3X7jTQ82XLW7fWnST7j6XF1ZgEyQHFkIonVwY3355mjuuiuUzn4YN96kPEjD6KXqxTq6gJhmtANEeESSPUM02RnAxAfH0nmhJHGsRB1U14ydwMAEem+u5Lbe6cD4CgM7s0+nUrKiCK9n/frnlqiOCqJxJoS/xjVQarWGtkFYg4+yOe65Sl9sKBh925/myV0UzoSq9CilLpZKbVEKVWqlHI2ODdGKfW8Umqof80UAkV+fhWplXm4lAUyM+uPf7C0BjgwBBRqrF60CYCoXh1YA5VoDC8GO0L80CQnpZEJpKV1LMCuh9KoJJJrS/xjVAdxbd1OCXbSDurrc93qXrKWS/ANn4RLKRUJfA08BQzGCLTbcOZ+OzAduNRP9gkB5umnV5BNEbt0UqOQQUWuKEqwU7Fmc+uVg0gSNZQRRXS8j2lNAOLicKGwVQY30G5qXQkV8R13zPCQ644jsbbEJ28+f6K1Zvv8X9hBMn36xPpcv7a3IVyylkvwFl97XHdiZDv+M5AB/LvhSa11CbAQONUfxgmBJzs7gWyKcWY2jvk3depAI+ZciD4FJ1FNCdG43R1wYLVYqLDGYA9yOpA0RzGlMR2PmuHh13wrCdTy0/ytfrDKd77+eie9a/LYQQrp6b73Hh29M3GjcGwR4RK8w1fhuhRYrLV+SGvtpmW39+1ghl7wAqXU+UqpZ5RSi5RSZUoprZSa1U6diUqpz5VSRUqpaqXUr0qpGUqpVr0ZlVJnKKW+M4c3K5RSPymlrmjnOlcopZaa5UvN+md4+9nCAbvdRjbFZBxxcKPj9903gZ0kY90bmvMOyVRTTHSHnvAByiPjiakOrnClO0spj+t8j2s/8QBEFAXHA1S73WSbi499iZrhwRptJ5d49i2VEKiCd/j6VzYQWNJOmSLAl8fI/wNuwlgT1q6bl1LqbIxe3bHAHOBZIBL4J/BeK3VuAj4BRgGzgFeAvsBMpdQTrdR5ApgJ9DHLzwIOAT4x2+sWOCpr6EsZ2vQo9JCaGs1OkonM3RMky9omK8FNXXQ8Awb47g4PUGlPIKY6iC7YFRXE6xoqEjoe7snDPlO4LHnBWQuV4KggjroDnqg+snt3GTtJZvt3v/rZssCwcWMRWstS1WDiq3DVAEntlMkCSnxo81ZgKJAA/K6tgkqpBAwRcQGTtdZXa63vxBC9H4HzlVIXN6mTDTyBIajjtNa/11rfCowGtgK3K6UmNKkzEbjdPD9aa32r1vr3wOFmO0+Y7YY91rxcLGhU/8apNVJT7ewlEXttJVQGPyBtU+IcFR2KmuGhLj6J2Nog9rhyjdiCVYm+R7dvSqHVCHv18M0fd7qtjmDbYyQ7P+P3J3aovtOp2UkyWTr0E0rOmrWO4cNf44031gbblB6Nr8K1CjjFdNJohlIqEWN+a6m3DWqt52utN2vvHmHOB9KB97TWyxq0UYPRc4Pm4jcdI77is1rrHQ3qFAOPmm9vaFLH8/4vZjlPnR3Ac2Z7V3lhb8jz3P1fAmDJbOwNlpRkJwczDmBu8AO4NsTt1sTWlOFK7LgbuSsxmQRnZdAcGjzfaXVS54cK/7v4WgAyKMfhcHW6PV+x7THmQTOOHNGh+nfdNZ6dJNOfUpx1oZeXNiengrffXseMGfO47LLPAbjqqrkdm18V/IKvwvUy0B942+z91KOUSsIYWksGXvSHcS1wgrmf28K5hRhxFCcqpRrmVWirzhdNynSmTljSB6PXYevXWLgsFkWuOQRVsWlHV5vVJsVF1aToSuyZGR1uw5WUTCpVVFY6/GiZ97j3Goks61I73+PqN9aIeNKb8qB8nsi9Ro/LnZXdofp9+sSxkyQicbFtcejNc5188gf89ref8/TTKxod//rrHcExSPA55NO7GOJ0HpCP2btRSi0DcoGzgee11p/718x6hpn7TS3Y5sRwDLEBg7yskwtUAv2UUjEASqlYIBOoMM83xeMf3i3WqvXGcAmPyMpsdi4HYyhuw7w1XWpTexTvLSIaJ6pXx3srruRU4qmlvLDCj5Z5z19+/18AalM7Lr71RESQRyx9KaOiouuFy7V1O0VEE5HW8R7wbU+dB0DuknX+MstvrF9vRGkZRAFTWY8d4zv++uvQ9LjtCfjsAqS1no4x/LYOY9hOAWOBLcDVWuub/WphYzyTGq3NqnuOJ3WgTmKTvS/XaIZS6jql1DKl1LL8EI33B0aPy43C2qd3s3OeHldybXAjTDTFk3vKlpHe8UZSDaeIqt1dn7pFa01UcR412Jj9rX/+NmrT+9KfUioquj5s0p7vV7ODFCIjOx6iNGmM8Xzp2BjcKPdNueOO79AaLmcZG3mcz3iNn3maXpTz5JPL2m9ACAgdCvmktZ6ptR4DxAH9gHit9SFa69f9al0Yo7V+WWs9Tms9Lj29EzfYAHNQXC3l9sRGi489PDfrfKqxYcsLrbxcG38wOs8ZI/q3U7J1LOmGcNXmdH2g3bo6F30oJ5d4jpvs9cqRNokYPID+lJCfX+2X9nxhIEXsIJnIyI5HkEs+fCRuFJZtwVmL1hpPPrmMg8jnFWazgEFcwiUMppBXmB1s03o0nYpVqLWu1lrnaK27yu2sae+oKZ7jJR2oU9pk78s1wpbkuhLKY1t2Yx5yUDI5JGALkpt1a7jzjF5K2vCOC5fNHGYs2tz1gXarqpxkUkouCTzwwIT2K3hB4sghZFHCjz/m+KU9b1k4fweDKGQzaaSldSCKiYktNpocWwoRO0JLuAAe43NqsXEJl/AuY3iAUziLdUxkO7W1oedM0hMItyC7G819s/klpZQNY52ZE9jmZZ0+GHnF9mitqwBMEd4LxJnnm+KJItpsziwcSXeWUhHX8loiu91GLglEFIRWjyuirAQAa6+O92Qj+xhOEU/d/7U/TPKJqioHWZSQH51GVJSvKfFaJnrYIBKpIX9r1/5Wleu3EoWLTaQTG+t9Hq6WKEjuR2J+aC14H0gh57KGz4eeRZ45dP4cR5NHLHczn5tvnhdkC3smvsYq3OblFqjHJs9fyWktnDsWiAF+0FrXellnSpMynakTdmit6eUuozqxNeGykkMCEfmh1eOKqigxXqR2fPGuxyMxlSo/WOQbleW19KeEISce6r9G+xu9T8e2rnUYiM/dAcAdL3qfg6s1qvsOoF/NvpBZ3FtSUsOD/TfgQnH2p39D6zvIyoqnikhe4iimsIH134THounuhq89LguGM0bTLRnINrfIDrTrLbOBAuBipdQ4z0GllB14xHz7QpM6rwO1wE0NFw0rpZKBe8y3Td33Pe/vNct56mQDvzfbC/v5vNoaJ70ppya55Z5LTEwEuSRg2R9aPS67R7hSOh7nz5ZhDBWm0vWLq2t37SUKF86+HR/qbIYpXPa8rh36tO824gtGjBze6baq+w0khWoqd4XG31tq8r84cfc8liSNxn7QQAD++McjAEi6+VqsaC5yiINGMPBpnEJrnd3aOaXUEOBfGENvXgfZVUpNA6aZbz2ubROUUjPN1wVa6zvM65cppa7FELDvlFLvYUSyOAvD7X028H4Tm7crpe40bVumlHofqMNYzNwPeFJr/WOTOj8opf4B3Ab8qpSajSHIF2GEs7q54WLmcKVqz35ScFGX0vJaoszMOHJIIJEaqKiAuLgutrBl7JWllKpoEiMiOtxGn8G9qMYWlB6X0wwmq7MGtFPSB7IMJ4/Yoq7tHUfv2U45UUT061wyTADnwMEAlK9YS1wnk2t2lqoqB0eyi0zKuKNkJEebx2+8cQw33jgGp9PNkmfu4qg9i9mzp5x+/eKDam9Pw289I631FuBcjDVQD/hQ9TDgCnPzCN6gBsfOb3Kdj4DjMBYcnwfcDDgwRObiliJwaK2fwRC3tcDlwHXAPuBKjyi2UOd2jOgY+8zyl5v1z9RaP+vD5wtZitfvAMDRyiJYpRTJBxtPmqEUPSO6qpRiS8eC63qIiLBSaU9gaErXR87QO4zhPOugbP812qcPbmXBnp/TpREddn27nE2kYY/u+ENEPUOMhdRFP63ufFudpKiohrNYhwMLX9C8N2mzWdg99gTGsYc/X/ufIFjYs/HrkJ4ZeulrwOsBb631g1pr1caW3UKdxVrrqVrrZK11tOmK/0+tdavxbrTWn2itj9Nax2utY7XW47XWb7Rj20yzXKxZ7zit9afefrZQ56s3jHjJ+y0JrZZxppsLZHO61lutLWKqSymxdr73VxmVQFxd1+fkspix/SLN4Se/YLNRldyLjNoCdu7smuDBZWW1DCWfzaQRFdXxNVwe9MCBuFDM/tvn1NQE11uvutrJWaxlAYOYfusxLZY57okbAbDO/ZytW0u60DohEHNRTg4M+QkhTKbViBox5ao2XLL7mEM2+0Jj3gEgtqaMElvnhasqOoEER9cLV2TOHgqJIbZ35yPDN8TRP5vBFLJjR+CDB9fVuUhL/CfZFLOJdL8IV3RSLDtJZjj5QQvF5cG9aRMHk8f/GInVqlosE3/koWwlhTNZx5Ah/26xjBAY/CpcSqk04BwgtHxahRaxmuuzkkdkt1omor8Rw9C1J3R6XHE1ZZRFdH5OoVDFEldbzs8/d90waG2tk70LVrKFVOLjO+c+3hQ9dChDKSAvL/Dzdq++upohFGJFG0OF9s679cfERLCG3oxiH9XVwe1xuT43wpR+xnD69m35ISk6JpK5DGMyW4lA1nN1JT79tSml7m+jnf4YsQoTgbs7aZfQBUQW5VFFJDEJrQ8VxvbPoBYrzh176Nyskp/QmoTaUsriWrfZWzYUWDiHCg4+4m1amer0O8XFxvDaQgYxLs4P80INsA4fRgqVLJ27losu6ryXX1tUVTk4BEPw9yZlo1TLvRJfsNksrKY3U9nA9pIKCKLDw7rnPySWJLaRyi23jG213PKkQ/h9yY8cya4utE7w9THpwXbOlwGPaK3/3jFzhK4ktqKIfFsCA9q46aSlx7KPeBJ3hUhCyYoK7K5a3L06H5y2MCKJVEdVlz4t15WWM4ASNpGG1erfkXr7aCOtyPczF/DhGYeQnGznhBP8E1KqKW635hD24cRCv5OP8EubsbERrKYPNty8cc9/eeR/v/dLuz7jdnM8W/mEEcyZM63N3+mqN27CdfYrXNkvtNY6dnd8Fa7jWznuBoqBDWaUdiHE0VpTs303u4mjLafstLRo9hFPfG5ozHG5cnKxAlEtRLP3laLIJKiCXnRdhHjHeiPgyib8H78y6pCDARhKPuef/z+AgPUktYZDyGUTaTz376l+afOgg5I59JLJ8M7bbPvkB4wlk12Pc+UvpFLFfIZwVVJUm2UnnXUIq+0DGLpnJcuX7+Pww2V6vyvwNa3Jgla2RVrrNSJa4UNJSS29Ka9P+94aaWnR5JKANUQWIeevNqJ5JRzU+cW7pfYkADK6ULj0xsAJFwMHoq1WhhL4bASffLKVQ9jHr/QhIaHtm7svnHfPmdRhrR+GDAaf3fkqAJXjj2GyF0GQP6sZyFHs5LhxrwbaNMGkXeFSSkUqpZYqpb5VSrU6KG+W+1YptaStckJoUFFRZ0Yob3uuyNPjiiwOjdQs+381Fu+mjuq8K/ld/zgTOJCTrCuwbDJCZ26m85mPmxEZiRo4kGENhGvFisAMYa36fiuDKGI1/l0oPHRkBoXpAxjFPq677iu/tu018+ezhVR6jRvWflngW4YQgZtJbA+wYYIHb3pcvwUOx4gw0aqPqta6DngcOAK41D/mCYGioqCMZKq97HHFE11eDI7guigD7PzJ6LEMObrzzgeDJhg3pgzKcbm6ZiFyxNrVbCWFSvzXS2nEiBGM5IBY7doVGNf4QzB64KsDsPKlfMBQRpPLK6/8isPR6tLMgJC7p5Tj2MZ8BuN0ereQe9rfr6AOK8cTWrnEujPeCNe5wDZvshprrediZAi+oLOGCYGleofh3t6ecNntNiL7m/NJ+4M/AR1RlI8LRcpQP8T5yzAcPHpT3mXZbB3LVrCSzs/PtcrYsQwnjxiMhJLnnPNxQILWTk0x8phFTDjS721HTjiCAZTQi/IuX8+18+OFJFHDfIbgdHr3MOOKimYJWRxP6KVk6a54I1xjgO98aHMhRhgnIYRx7DKCsea2I1wA/cYbGWF0CIR9iq0ootASD9bOL3glOpq66Dh6U86UKR92vr32KCtjkCufVfTl3XfPCMw1xozBimY0B9bdFRXV+P0yo6u2kRuZyutzr/J72zEnTALgKHZ2uXDF/PQ9APMZ7HUvfNq0IcxnMGPZCyUlAbRO8OCNcKUBvjxq7wf8GxJA8DuuvcaN7co/nthuWc8i5MotwV+rUrxhJ7lu/60oq01K6zLnjF2fGjfFuoNHc/HFAVpnNdZYczSmgXClpT1HRUWdXy9zaM029g88xK+OGR6iJ46nDitHsYuKiq4TrpUr95OwfDEbSGcfCZxzzkHtVwKyshKInnIyVjSu774LrJEC4J1wVQO+xNeJA/z/iCf4F7P3NPG8w9stGjnAGNqq2hbctVwul9srT0hfcKamd5lzxr8ufR6A7SlDAneRfv0gLY3xTYLXnHfe//x2ifLNOxlAMfkDD/Fbmw2JS09gFX1N4fKv4LbGZ59tZfzYmaSsW8Z8BpOb+zumTfNOuACKhx1GDTZWPvl++4WFTuONcO0GxrVb6gDjQJaRhzqWvDzcKGKzvZhvyTCjxwd5qNDhCIBwpfUio4uE6zi2sZF09jgDGINEKZg0iclN5lu++moHixb558HjnslGfIHaMeP90l5TlFIsIYsj2MWWdV0zr7plSwlj2UsCtcxnCAkJvoXjOv28g/mBAdi+XxAgC4WGeCNc32Hkx2pXvJRShwMTgfmdtEsIMM49ORSoWBJT27+JRsTGUEAMKi+4a7nqap1kUMF+nwYA2iaifyZ9ukK4XC4msZ0FDGLJkgA/AEyezECKGUBRo8PHHvseu3d3zsuwrKyWETkrKSeKqfdf2Km22uLqWbcQi4PnL3+RJ574OWDX8RARYal3rljAIGJifFvRM3p0OvMYwmHkQmFhIEwUGuCNcD0LaOADpdSI1goppYYDHwAu4Hn/mCcEiojC/RRHJWOxtB9jLirKai5CDq5XoWNfHnac7CXRb20mjRxEIjX0ig7s2nm9YgVJ1LCAQQG9DgDHGwFuTmjBPfuiiz5l+vS5HfI0/OabnQwc+AqnsIncYeOx2APk0g/YTjkJBxZOYyN33rmA2lonq1blBWzZwqpV+UxmK2vJ4Njz2x8+b0pCQhSlYyYabxZIryvQtCtcWuuNwEPAAGClUmqWUmq6UuoUc7tKKTULWAlkA3826wghTGJVMSX2ZK/KRkVZ2Uc8tvzgCpd7hzECvZsk/zU6wAh4lV5TGNAEjDXv/AcnFubi3aLWTjFyJLtJZBprm5368cccXn99DaWltT43e/LJH5BUtIchFNLr0rP8YWmrRKYls5hsTsO4ldx332LGjHmTxx5bGpDrvf7KSiaxne8YxLZtHctnljfgYKpUJMyXAadA41XIJ631Q8D/AQq4BHgF+MLc/m0eswD3aq0fCYypgj9Jri2mNMZ74colHkeQU5u4dxrCtcufwmWmvO+niwPnCKA1lv9+yDwGU9QVMfYtFuanH8VpbCSJllOc5OR470m5e3cZNTWGd98F/ApA0mWBXaqplGIuwxhDDv0o4emnVwCwfLn/H55eeGEVh7OHOOqYzxBmzz6zQ+04LREs0tmUfvylny0UmuJ1rEKt9aPAUOBhjDmsDeY2H6NHNlRr/ddAGCn4GbeblLoSymK9CzsUEWFhH/HElRcY0VWDRNVGI6SOX3tcpnANoJgXXljlv3YbsmgRUbu28765vHHWLP8EpW2LC2bfTyQuch5P5Oyzm3sx5uZW1r/WWrNsWcvzl1u3lpCV9TJ//OMiAC7iF/ZkHgzZ2QGxuyEfMBqAS1hJXZ0RQaO1pI6d4cYbv2k0vzVwYFKH2vnvfzczn8Ek7t4MeXl+tFBoiq9BdndqrR/QWp+ktR5pbidprR/UWndN6AGh0+St30UEbhZs9m5ep6LCQS4JROEK6gLLb15fRA02/vDwqf5rtE8f3BYrWZTwpz8t8l+7DXn8cWrik3mHMTz++HFceunBgblOA6InTYAjjyT6+acZPqS5F+ZDD/3I7Nkb2b+/kszMFxk/fhbvvLO+WblJk94F4F//WsGh7GUMOURf9duA2w+wjTQWM4DLWY4xzQ4WiyIvr5Lvv/fv0ozJbGU1vbnloY7/bd188xjmYzwkVHwapDiLPQT/JgQSwoLqrcYan+Io74YKDz449YALehBd4pPK8thDInkFflwmaLVSm96bLEoA/D9c+OWX8OmnLJt4ITVEcO21o/3bfmsoBffcA9u383D6Cl588WSGDEmqP71w4R4uuOATxo+fVd/7uvTSz8jPbzy02LBn9ge+p5IInJde1iUfAeDfHMlI9nMyRozK//xnI4MG/ZtJk97zSyir6moHUTiYxHbmMYSrr+742rR//vN4lpNJGVFseH52p20TWkeEqwdStdWYK7rq3lO8Kp+REcvR5kLlvF+DF4/tsNRadpPEddf5+eafNaBeuB566Ef/tbttG1x5JQwbxqfDzyU62ubz+qBOceaZcOaZRDx4P9ePKGX9+un8/e/HNiqye3fjpQDnnPMxr722mtraxr3xIeRzKSt5nfFE9Q5AZPsWWLHiMt5hDHtI5D6+xdPr8oSBKi/v3ENGTY2TmJinOY5txOBgLsNITOz472O1WnBhZREDiV/+Q0BiRAoGIlw9kL/fZkRRcGd4H9m7PM6I4vXE7Z8ExCZv6FWVzx5LMgcf7N8bp33YIA6KMm7gjz/uhzVDVVXw6qtw5JFQVwcffsiOfbX07h3rlxT3XqOUYUd2Npx6Kra/P8adV7cdtWPx4r1cffWXnHHGHAYNeqX++ON8Ri02HuFEoqN9zT/bMcaMyWDXvpt5Z8j5TGI7l7G80fnExGeYPPm9DrdfWFgNwBQ2UI2N/5U+RWxs5x4sfvnlCuYxhGHkU7NVZk8CRdf8BQohRV+MRai+CFdFXIrxYl+QhgqrqkisyGdHlP+jNaisLHo7i7Hgxo0Fp9ONzdaBZ7q1a2H6dFixApxOGD8e3nyT4oxs3n//C6ZO7XwOMZ9JTzfWFV13Hdx7L9x3H7XDhvPNRjf5xFKLDScWHFhxodCYwvpNw0Y0u0jir5Gn8shz5xMV1XW3jYyMWO5Y/xKuY9Yxc81nrKrsy2r61p9fsGAPeXmVREdHEB/vvejs2lXG9OlzATiNjXzHYKYkdH5h++jR6XyNESrqloPu5KmKtzsthkJzpMfVA+lDGUVE44rwfgHp0VMOpooIDknr2vxI9WwxFtPujgpAavSBA7G4nPXDhTNmzOtYO716QVQU3HEHzJsHP/3EF9sjSUl5FoDFi4O0nKBXL/joI1i5Eu67j8ghA+mtKpjMVs5gPRfwK5eznGtZynX81MK2lA8ZzXe9JnLNNV00R9cAi82K9cMPsCQmsDjiFaaxGs+wIUBGxgskJPyLZcv24XK5+etff2LPnrajofzud1/z7be7yKaI4eRjO91/np5XPvEbtpPM2axl9uxNfmtXOID0uHogfSgjlwRcLu/H4KeePpgdlgQGRbe8LijgbN4MwLIy/8UprGe4Ean94kM1j/0CH3+8lWefPcn3dtLTYeHCRoemTv1v/eu+fbtgDVdbHHaYsQFjgQ8+2MiLL/7CvHnehRY9/iA/5EDrKJmZsHgxcdOmMeeXN1lJX/7LKFaSyV4SKSCW08a/xP8+O5+/3PMNX8xZx5TTh/DPp1aQl3cjAC+9tIrTTx9McnIUtRU12HBxBusAOPkfN/jP1H7xfMxIbmAJaVd+zBVXjPJb24KBCFcPpC9l5BLvdaI8D1WJaahghX0yhSsgKe+HGdEsDnLnAans2VPOwIEvs2rVFSQm+i+s0aefnuu3tvzBBRcM44ILhnHZZZ8za9a6RucWLbqYSZOM+aMffriEuXO3c9BB3nmhBozsbNTPP8Nrr5F8/z95OK8Fl/PTHzSS1PxsbHcDRM4A4HqAm41iDUdCC4gh7SDvI8G3R2JiFC8ykhl8zylsYvTomfz665V+a18Q4eqR9KGcRQwkxscQR46UXqRsD9LQx+bN7FfxXHfr0f5vOy0NUlIYUJ0DGOE4d+woIynpGbS+o8PNNg2rNGhQUieMDBx//eukZsJ19NGZXHTRMD76aAsTJvRlwoS+rdTuYiIi4PrreXLNEN549kdGsY8MykmjCjsOrGhsuLDhxmoOJ/720hHMerv5GjUPv//jJMORxU8MH57CIgZSRDTTWMMVqwOT/qUnI8LV09CavqqcHJ3AjDN8C/hak5JO/63LcLu1V8F5/Ylev4GNOi0giQtRCoYNo98O/81B5eZWcPfdAVrQ7Gc8btuZmXHcfPMYfv21AKUU773XsdBHXUFGRgzl2PmR7HbLPvo2QMvCe8klI+j72Ol+tS07O5GH/nIcc+4dxYX8yu/ompxi3jJr1joOPjiVsWMzgm1KhxHnjJ5GURGR2knS8GyfvcOcqb1IoZry/M6lxvAZtxvnipWsoq9PnmM+MXw4A2uaC9eYMW+ycqXvw6NjxrzJG28cCHKblBS4SOqdJT09BoCHHz6aP/7xSN5+27838kBw662HN/tbWLDgIjIzffMMfO659jOAd4R77jmKi7/4K/HUci6rcTrd3HDD16xZkx+Q67XFV1/toFev59i1q4zExH9x2WWfc/jhb7F/fyXTpn3Egw8u7nKbOosIVw/DscsIlZMy0nfXbEum8dRasGabX21ql61biaitZhV9A7eG6LDDiCwu4ORRjfMwrVqVx9ixb7Fhg285lvbvP+DEUlMzg337fucXMwOB3W5D6zu46qrwGdKKjY2ktPRmyspuYeHCi3n99dM49tj+7NnjvZNFr14xJCXZA2fjqSeQE53BlSzjb39byksv/cI113R9KKg//nEh+fnVDBjwMmVlB3p/Z5/9ER9/vIU//9mPi+67CBGuHkblFsODzNLPi8zHTUgcbqQAyft1u19tapdVq4wdff05FdGYcUae1IMrW/5sxxzzHtu2lbBihdH70lpTUtI89NTSpbksWdK45xYVZevStU89BaUU8fGRTJrUjyuvPOC5N3HigWHBiy8eXv/6rrvGM2vWVJYvv4zHHz+O7767KNAGsufEczmRLbz5f4Z3qd1uDew1W2DVqpYD/v7004E1mZ7F2OGC/Df1MGq3G3EKPb0nX4ge2A8A5+69frWpXVauxIGFdQRwTP6ww8BiYXjlDqAf6enR5Ocf+GcuLKzmkENmUlXlZM6cs1mxYj8PP7yEnJwb6NMnDofDxciRM9m8ubhRs8OGpQTOZqFFFi++pNH7qCgrb7yxlgcemFif2bir5neWHX42oz99hdtZyPWcj8WieOSRHxk7NoOpU7sgqSigcHMc25jEduKpZQupfMxI9pNQX2bv3gpSU6O7xB5/ID2uHoZzpyFcEQP6+VzXPtBcx7O/5RQYAeOHH1hNH2rxLZ26T8TEwMEHc352KYMGJfLDD5eQkRHTqEhVlRG/75xzPubhh5cA8OKLvwBQWFjTTLQAvv02sHmrhPaZOXMKWt9RL1pdSURmH2YyjstZTiYl7NlTwX33Leb00/9LWVktDkfgFvS/8sqv7J39FSt5ivm8xIN8zU0s5iX+y27LY/yNT4nCiPu4b19lO62FFiJcPY1du8knlrheiT5XjRnQFxcKW14XCldNDXrJEuYzGGg89ON3JkwgbcMKtm68iiFDkvnss/bXXT300I84HC7efbe5u/VLL51MZmYAFkwLYcP06aNIfPR+AP7C3EYPN4mJz3DRRZ8G5Lp1dS5+vu4hMi6YQjLVXM5FxPMIMTzKwdxBxBW/5S4WsIjn6UMpxcU1PPfcSnbt6mLHqw4iwtXDsObsZhdJHVpYG5cUTS4JbPzmlwBY1gpLlqBqa/mOwbz44skBnUzn5JOhrAyWGunhvXUEiYz8J7fd9l2z4+E09CIEBqvVwm/unsJTTOIKljOJxo5Nc+Zs5oMPNjJnzma/Xtf1+JO8zId8xVBGcTtvMY6/PzsFUKwnA157DT1nDodEFDKfl5hx8TvcdNO3nH32R361I1CIcPUwIvft7bBwWa0WdpLEAEq6LmXDV1/hVhYWMZCjjuoT2GudeKKxpusrw/OrswuGA26vEDY8wolsIo23eYdUGg/LXXjhJ5x77sfMn+9d6K12ee01ov/vj/yH0ZzFlZRjp2/fOG688TC++eYCFi/+DQBq2jTU3C/IpJRveJkEqpvlYwtVRLh6GLacPewiuX7tjq9UpvYhi+L6nEgBRWuYPZvN/Q7DERPPyJEBzgOVkgITJsCHH4LW2O021q27iscfP47nnz+J/PwbeeONKY281lrioYeOJjf3dzJMKNRTSRSXcAlpVDKXf5NIcy++1asLUOqJZlFMvOV//9vCh3e9hfPa69k7aiK/5Te4MLwYU1PtKKU48cQBTJx4wKM48vhjefaUPzOUfN7hHXL3lvH666tD3stQvAp7EqWlJFDLLpI6vJC317gR9P/yZwpKaoiLC3C6hqVLYfNm5h32OwYmJHYs1YivXH453HADLF8O48YxYkQqI0akNjg9klNOyaZPnxcaVfvii/NYuTKPtWsLmD59FL17BzmgrhBSfP75uezfX8X5V1UwhzdYxtNczsWNIn/84Q9GVoLLLvucykoHV1wxErvdu1v0p59u5aqz32YlT7GbeMauORUHNkaMSGH9+qJWHzSVUvzpyz+x6TY3p//zXv7K50yfbuGss7bw8cfndPpzBwrpcfUgHFuNNUoTLjii42307U8kLj55uQvCGT39NMTH89+ow7tOCC66COx2eOaZVov07h3Lq6+eSk6Osdj1kEPSOO20gdx995HMmnW69LSEZkyZMogrrhjJ54zgzSv/hR0nP/Acn/NvfstysimiYaqWG274mvj4f/Hdd94NH5515ofM4l0yKOd8LqcEY0TlxhsPw2az8NBDbcf4HPqPe3gt8mjuYgHn80ujBfShiPS4ehBla7aQCliyB3S4jcEnHgqvwxsPf8Gld50cuF7XsmXw3ntw112se8vBySd3kRgkJcFNN8GTT8Ltt8PolvNPTZ9uRJkoLPy910/FQs9GKVUftHnK1lgOX/Qu17OEtzCi8FcSQS4J5BNLDRHUOq2UHP8ae4/qR2a/eFAKl1YU/2Y6P9gGc9ZZQ7juuq/46adc7uVbprCR6zmXFRxY6mKsMbzNK/tmT/o9I77dzev8h5v0GP9/AX5E/uN6ENUbDY+mqIM6nok35TAjBUgWJRQXB2i4sKAALr4Y+vSh9MbbyPnbmwwd2oUpNe6+G954w+h9/fADJLd+7ZQU8RwUfCexbwp/4SQe5QQOJZcj2cVQ8ulDOalUYsdJMnXYKce9tozyPZEUFVUzvupa8mevBlbXt3USm/gzX/MWY3mZo+qPP/nkZM491/t0LW+9P43RaetYwVM8uuGfUHYNJCS0XzEIiHD1IBybt1GHlfghnUgImJUFwACKKS2tpb8/cws6nTB3LvzhD7B3L8ybx6PPGeuj2nOI8CspKfCf/8App8BRR8FTTxmvrV0frkfonrz44slMntwfm83C1KmDyMx8sfXC5ebWApmU8A7vsI5e3MB5gOKuu8bz0ku/cs01h6B8iJGWmhrNI69dyEXTC/im7GX2Tb2QxTOe5pxzh3Z5Noj2UF3m1tyDGTdunF62bFmwzWDb2JOoXbma2J2bycrq+JNUTWwir1cdzKiFHzBpku8ROJgzB+bPh+pqqKkx9vv2wdq1UFICgwbBrFkwYQLHH/8+3323G7f7dp/+Cf3CwoWGs8bOnUava9Qo6N0b4uMhOhosFsN9vul2//2Q6PsCb6HnUllZxymnzGbs2AyefXalV3UicbKAFxjJfsZzCxvpBUBFxS1ER0d0WGyOOeZdjlz8Dk/yKXdyOsNeeZRrrml5yDyQKKWWa63HtXROelw9iOg921mv0pjSr3PzRe6Bgxm8tpDc3IqONfDjj/DWW8bN37OlpsKFF8Kpp8IZZ0CkMQRZUFDNtGlDul60AI49FjZsgE8+ga+/Nl6vXg0VFVBVZbjru93GvuF2++0iXIJPxMZG1sdYHDEihd///ttmZaZOHYjVauGTT7YC8A/+x1Hs4gJ+y6EXTmLjfzYCEBMT0an/lz59YvkHx3Iku3iMz3nju5MhCMLVFiJcPQW3m5SiveTGHdPpbn/EqOEMXTuXd7eWdqyBv//d2LwgP78quNl37Xa44AJjE4Qu4MYbx3D00ZlYLIrRo98AYMiQJD777Dy01nz77S4sb7/FCTN/5HGOYzaHot8/k4qKOs49d2inH/KefvoEli3bx9U7LmAU+zjr3f/D8dDZWLMHhMyQobjD9xRyc4ly1VKcltXppiIOHk4WJaxeuhO3O3BDzVprCgqqSUsTBwihZ3Hoob045JB0NmyYDsDNN48FDM/Ek6J2ccK7j6CPP567mVJf57PPzuPqqzufU61v3zjmzJlGBXbO5Qoi3Q5WDD6WaVPe73Tb/kKEq6ew2YiFVtW3467w9QwdigXNrx8tZerUDzvfXivMm7cLl0t3TZQOQQhBhg1LISfnBm6+2XRPX7sWzjoLBg5EffABLqxcc43/E4AedpgxX7aRXlzJRRzJbs766p846px+v1ZHEOHqIeQsMCZ8K/pkd76xoUONHfl8+eWOzrfXArW1Tk466QOgiz0KBSHE6NMnzhj+q6qCKVOMOeG5cyE1lbq6W3nppVMCct2lSy/l7ruPZA6H8Agncg1LeSFqMnt2N48gX1ZW26WR5UW4egg/vrWQWqyc8bvjOt9YA+EaNCgwTggNPasuvHBYQK4hCGFFTAw89hh88QUMMEZOIiKsAZt3Gj++T30aofs4lSc4lltYzN6pF0NdHVde+QXPPWf8n44bN4sBA14OiB0tIcLVQ8iq2M0m0jnuhOzONxYXB337cvoQB7W1gUmEt2LFgXTjQfEoFIRQ5JJL4NBDu+xyo0enM3PmaQwcmMSdnMFDnMSRa75gUdRwvn7jB266yfB+bCmJaiAR4eohDCzbyfaEbP81OGIEA6v2kptbSWlprf/aNVm9Oh+A44/35wpnQRB85YorRrFt27U8/fSJPMCpXMSlHM4e1vIE17IEm/pbfVmlnuDpp5cH3CYRrp5AcTFp1QVsj+u8R2E9hx1Gn4JtKLeTk0/+wH/tAhs2FLJ6dQEAc+ZM82vbgiB0jFtuGYvWd/AfDmM0t7GSTF7mQ9aZAuZJ1TJjxvwWM4L7ExGunsBqI67ZzsRs/7V52GFY62oZTj4//7yPBQt2U1PjH4+j//73QDbYjiS8FAQhcDz66CS2ksaJXMc0rqCSSF7mQ/bzZz7n39zNt/z7kqfY9cNaY0F+ABDh6gn88gsAe1M6Hly3GWMM99zD2AvA5Mnvc/31X3W62aKiau699/tOtyMIQmC4++4jyc+/kcFDUviYUYxlBuO5heeZSBYlPMpcvuVlso4eRWViOtrt9rsNEjmjB6CXLGGfJZGqpHT/NTpsGNjtHF6zl7c5HIA331xHv37x/OUvkzrc7I8/5vjLQkEQAkRaWgybN1/DjTd+zQsv/MIy+rMMYz46hUrGsJeDKCC2vI5HndoTwc1vSI+rB+BcsIjv3QM49jg/OjrYbDB+PKfE7m10+NFHf+Knn3Lr3zscrnaT4b3zznoefvhHtNZ89dVOwPBmuvvuI/1nryAIfufxx4/jpZdOxu2+nVdfPZVbbz2cImL5lqG8yEQuWPIKkZH+z6og0eG7gKBGh9+7F/r1YwZncenSFxk/vo//2r7vPvRf/0qC60EqsDc6dfXVh/Dyy6eQnv4cRUU19O4dy5o1V1JV5aCkpJYNG4pISbFz4okDUOqJZk17Eu4JghB+LFy4m8LCGs45x/t8YE2R6PA9mUWLAFhMNn8dlebfto89FvXII5TPHU/29TvZufPAyvlXX13Nq68eSHa3b18lDzywmNmzNzVKCz537nnNmj30UD8OaQqC0OUce2xgl7HIUGE3R3/6KfnEsoJMoqMj/Nv4hAlG+pEvv2TVqsvbLf7cc6saiRbAaac1j3X4yiuBCWEjCEL3QISrO+N0oj/7nM8Yzu13BmC+KC4OTj4Z5swhKTEKre+gtPTmDjd3553j+emnS/07nCkIQrdDhKs78+23WEqK+R8jOeQQPw8Tejj3XNixA1asACAhIYrS0pvZs+d6+vc/kLBy+vRRjarZ7VZ++unS+vcLF17M3/9+HEccIaIlCELbyBxXd+bFF6mISeKzqhFc3ysmMNeYNg1uuglefBFeeQUwxCshIYpdu64nJ6eCHTtKmTgxk9tuG8dnn21j9+5yrrxyJIcf3psff7wEt1szcWJmYOwTBKHbIcLVXdm4Ef3JJzzrOpY6bGRkxAbmOikpcNll8MYb8PDD0Lt3o9N9+8bRt28cACNHpjFyZOOe31FHScoSQRB8Q4YKuyNaw623QkwM/+BYhgxJCqyn3h13GNe8/fbAXUMQBMFEhKu7oTU88AB88QVbr7iNfOJ4+OFjApsa5KCD4E9/gnfegaefDtx1BEEQEOHqPrhcsHgxnHOOMWR35ZUc9KwxPDhmTK/AX/++++Dss2HGDCNn0LJlhk2CIAh+Rua4QpmLLjLSdbtcxuZ2H3jd8H1pKezcCTU1hov6Y49RfsMfYOazAAwdmhx4W202+OADeOghePJJePddiI2FrCxITTXWe0VEGJvFAmPHGj1DQRAEHxHhCmV27oS6OrBajc1iOfA6KurAsawsOPNMGD8eTj0VEhN5/V+Ge/qJJ2Z1XQbhiAijtzdjBsydC0uXGiGnioqMz1FRAQ6HMZzZr1/X2CQIQrdDYhW2g1KqH/AQcBqQCuQCHwF/1lp7la+6q2MVut0aq/VJ8/XtXSdcgiAIfkJiFXYQpdRg4AegF/AxsAE4AvgDcJpS6mitdWEQTWyG0+lm4sR3AGOIUERLEITuhghX2zyPIVq3aK2f8RxUSv0DuBX4C3BDkGxrkYiIf9S//vnn3wbREkEQhMAgXoWtYPa2TgF2AM81Of0AUAlcppQK0Mpe33niiZ/rXz/22CQSEiTtvSAI3Q/pcbXO8eb+K611o9zTWutypdRiDGE7Cvi2q40DWL++kB07Srnkks8oKamtP37mmYO5664jgmGSIAhCwBHhap1h5n5TK+c3YwjXUAIkXIcf/hYVFXU4HG6cTjcOh5vKSgcWi6K0tLbFOp99di5Tpw4KhDmCIAghgQhX6ySa+9JWznuOJ7V0Uil1HXAdQFZWVocMGDEiBZdLY7NZiIiwYLNZsNkUubmV5OVVsXVrCVarhYsvHsZVV41i1ChJwCgIQvdHhCtAaK1fBl4Gwx2+I23MmnW6X20SBEHoDohzRut4elSJrZz3HC8JvCmCIAiCBxGu1tlo7oe2cv4gc9/aHJggCIIQAES4Wme+uT9FKdXoe1JKxQNHA1XAkq42TBAEoScjwtUKWuutwFdANvD7Jqf/DMQCb2mtK7vYNEEQhB6NOGe0zY0YIZ/+pZQ6EVgPHImxxmsTcG8QbRMEQeiRSI+rDcxe1zhgJoZg3Q4MBp4Gjgq1OIWCIAg9AelxtYPWejdwVbDtEARBEAykxyUIgiCEFZKPqwtQSuUDO4NthxekAQXBNkLwCfnNwgv5vbxngNa6xXBAIlxCPUqpZa0lbhNCE/nNwgv5vfyDDBUKgiAIYYUIlyAIghBWiHAJDXk52AYIPiO/WXghv5cfkDkuQRAEIayQHpcgCIIQVohwCYIgCGGFCJcgCIIQVohw9XCUUv2UUq8ppXKUUrVKqR1KqaeUUsnBtq07oJQ6Xyn1jFJqkVKqTCmllVKz2qkzUSn1uVKqSClVrZT6VSk1QyllbaPOGUqp75RSpUqpCqXUT0qpK9q5zhVKqaVm+VKz/hkd/azdAaVUqlLqGqXUHKXUFvP7L1VKfa+UurppiqMG9eQ360q01rL10A0jYPB+QAMfAY8B88z3G4DUYNsY7huwyvw+yzGyC2hgVhvlzwacQAXwKvC4+Vto4INW6txkni8AngP+Cew2jz3RSp0nzPO7zfLPAYXmsZuC/b0F8fe6wfwOcoC3gb8Cr2FkOtfAbEynNvnNgvg7BdsA2YL448OX5h/9zU2O/8M8/mKwbQz3DSMFzkGAAia3JVxAApAH1ALjGhy3Y6TX0cDFTepkAzXmDSy7wfFkYItZZ0KTOhPN41uA5CZtFZrtZXfmc4frBpwAnAlYmhzvDewyv7fz5DcL7iZDhT0UpdRg4BRgB8aTW0MeACqBy5RSsV1sWrdCaz1fa71Zm3eZdjgfSAfe01ova9BGDfB/5tvfNakzHYgCntVa72hQpxh41Hx7Q5M6nvd/Mct56uzA+FuIoodmRNBaz9Naf6K1djc5vg940Xw7ucEp+c2CgAhXz+V4c/9VC/+k5cBiIAY4qqsN68GcYO7ntnBuIVAFTFRKRXlZ54smZTpTRwCHuXc2OCa/WRAQ4eq5DDP3m1o5v9ncD+0CWwSDVn8TrbUT2I6RQ2+Ql3VyMXrO/ZRSMQBmDzoTqDDPN0V+9xZQStmAy823DcVDfrMgIMLVc0k096WtnPccTwq8KYJJR34Tb+skNtnL7+4bjwGjgM+11l82OC6/WRAQ4RIEQWgDpdQtwO0YnoKXBdkcARGunkzTp7qmeI6XBN4UwaQjv4m3dUqb7OV39wKl1E3A08A64HitdVGTIvKbBQERrp7LRnPf2rj4Qea+tTkwwf+0+puYcywDMRwDtnlZpw8QC+zRWlcBaK0rgb1AnHm+KfK7myilZgDPAGswRGtfC8XkNwsCIlw9l/nm/pSm0QCUUvHA0RgeUUu62rAezDxzf1oL547F8PL8QWtd62WdKU3KdKZOj0Ip9UeMRb6rMEQrr5Wi8psFg2AvJJMteBuyALmrv+/JtL8AOR/fFrMORBaz+vt3us/8fpYBKe2Uld8sCJvk4+rBmIuQfwB6AR9jhCQ6EmON1yZgota6MHgWhj9KqWnANPNtb+BUjGGjReaxAq31HU3Kz8a4Eb0HFAFnYbhQzwYu1E3+aZVSNwP/wriBvQ/UYSyM7Qc82bD9BnWeBG4D9pjtRgIXAakYDzLPduqDhylmrMCZgAtjmLAlT74dWuuZDepMQ36zriXYyilbcDegP/A6kIvxz7MTeIoGT3Wyder7fRDjSbm1bUcLdY4GPgeKgWpgNXArYG3jOmcCCzBiIlYCPwNXtGPblWa5SrPeAuCMYH9nIf57aeA7+c2Cu0mPSxAEQQgrxDlDEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuAShB6OUmqmU0kqp7GDbIgjeIsIlCB1EKTVcKfWMUmqNUqpUKVWnlMpRSn2mlLq6Sbr2bodSqr9SymUK36PtlNVNtlqlVL5SaoVS6t9KqSlKKWs7bcQqpWYopeYppfLM77tEKbVUKfUXpdSgtuoL3QeJnCEIHUApdT/wAMbD348YAVkrgAyMYLqDgOVa63HBstEblFIzgSuAgVrrHT7W/TNwP0YYpP1Af22kq2+prOdG82dzb8XI2DsSI1xSJMZ3eKnWull6DqXUURjx+TIxYvV9C+RgpAAZgxGE1g0cpbVe4cvnEMIPW7ANEIRwQyl1D8YNeDdwgdb6pxbKnIGRNbdbYvaOpgNlwCzgRozAsv9tq57W+sEW2srACGh7AfCNUmqcbpBGRCk1HCOTQRzwJ4wgtM4mbQwE/oYRrV3o5shQoSD4gDkX9CDgAKa2JFoAWutPaZI7SSl1pVLqQ6XUNqVUtVKqTCm1WCn121au9Z05rGZTSt2jlNpsDrHtVkr9TSkV2UKdaUqpWUqpTUqpSnNbrpS6pWnetU4yBSOS+fvAC+axazvSkNZ6P3Ax8B1G0Od7mhR5BkOQ/qa1/ltLvTqt9Xat9YUYvV+hmyPCJQi+cRUQAXyotV7TVkHdOHkgGDf4AcBCjAj875nv31JKPdxGU+8AN2OkQnkBI/r4XcBLLZR9DBgL/IRxw38To6fyNPBGW/b6yHXmfqb5PSzHSEo6oCONaa3dwCPm298opRTU96ROwkgZ8ncv2mn6nQvdEBkqFATfOMbcf9uBuqO01lsbHjB7TV8Af1JKvai13ttCvcHASK11kVnnXuAX4HKl1N26cUr501u4hgUjdc3lSqlnW+sleotSKhOYCmzSWv9gHp4JHA5cg5GIsSN8j5HmvhdGgsTtHPi+l2utSzrYrtDNkB6XIPhGH3O/x9eKTQXFPFYHPIfxEHliK1X/6BEts04l8DbG/28j549WruHG6HGBkciys0zHcK6Y2eDYOxj53Ka35x3YGmZvyZO4NN3cd/j7FrovIlyC0EUopbKUUs8ppTYopao8ruHAh2aRzFaqLmvh2G5zn9zkGqlKqceUUr8qpSoaXGN5O9fw9jNYgKsxPPje9Bw3hfUToC9wemcu4WmyE20I3RwZKhQE38gFRuCjAJhrjJZiCM0i4CuMtPAujGGxK4AW1321MkTmcVCo790opZIwsuMONK/1JkYaeSeG6/kfWruGD5yKMS/3ZQvDmjOB8zDmv/7na8NKKTuQYr7NN/e55r5Tgit0L0S4BME3vgdOwBjWe9WHercBqcBVWuuZDU8opX6DIVyd5RoM0fpzU7dzpdQEDOHqLB6njFMbrM1qymlKqf5a692tnG+NYzDuSfsbrCn73tyPU0olaq1LfWxT6IbIUKEg+MbrGK7w5ymlDm6rYJPIGUPM/YctFD3OT7YF9BpKqd7AGRhrt15tZVuM0Quc7mPbFuBe8+07nuNa6+3AN4AduNOLdrp1tBLBQIRLEHzA7Ak8iBHp4TOlVIuRMZRSp2F4C3rYYe4nNyl3KkZPyR+0do0xwN1+aH86Ro/oba31NS1twJUY81NXe7tuTCnVC2NpwGRgF9A0fNQtGGJ5t1LqdqVUs5Eic/7wPWBCBz+bEEbIUKEg+IjW+lHz5vkA8LNS6gcah3w6FjiIxk4Vz2OsAftAKTUbI1zRKIxFyv8BLvKDaW9i9EqeUkodD2w27TgDI6JFh69hrqvyCOy/Wyuntd6ilFqAIUJTgM+atPOg+dLCgZBPx2A8CCzFCPlU0KTN9abAzwaeAP6glGoY8ulQjLBRGiN6htDNEeEShA6gtX5IKfUBRqij4zFEyY7hzr0K4wY6q0H5X00xeQTD686GsRbrXKAEPwiX1jpHKTUJYxHyMRiOFBtMG7/p5DVOwpg/W+lFLMBXMITrOpoIF4bYg+E6Xw7sxBDcD4GvTNf9Zmitl5ihn67FCC11OoajSxWwBXgSeNkcWhS6ORJkVxAEQQgrZI5LEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCtEuARBEISwQoRLEARBCCv+HxQ2gk/QavdsAAAAAElFTkSuQmCC\n", - "text/plain": [ - "<Figure size 432x288 with 1 Axes>" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - } - ], - "source": [ - "fit(0,250)\n", - "fit(250,500)\n", - "fit(500,750)\n", - "fit(750,1000)\n", - "fit(1000,1200)\n", - "fit(1200, 1400)" - ] - }, { "cell_type": "markdown", "id": "12ee35ee", @@ -310,110 +226,11 @@ }, { "cell_type": "code", - "execution_count": 10, - "id": "8f6d8519", + "execution_count": null, + "id": "2e7fe648", "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAApkAAAHyCAYAAACgdNVWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB1NklEQVR4nO3dd5wcdf3H8dfncqkXUmihk9CRXkVqABWkl9AVBAVRQEpEARXBHypKE0FEEIlIlx6qaDiqKXRRhAQSekIgJCHJpd19f398Z7J7e7t7W2Zmd27fz8fjHt/d2ZnZ781tcu/7fuf7/ZpzDhERERGRKDXVugIiIiIi0vMoZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiESuudYVaAQrrriiGz58eKzvMX/+fFpaWmJ9j0ak6xofXdt46LrGR9c2Hrqu8ajmur7wwgufOOdWqrYOCpkJGD58OM8//3ys79Ha2srIkSNjfY9GpOsaH13beOi6xkfXNh66rvGo5rqa2TtR1EHd5SIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5ApIiIiIpFTyBQRERGRyClkioiIiEjkFDJFREREJHIKmSIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5Ap6eIcTJgAhx0GLS3Q1OTLww+HiRP96yIiIlJzCpmSHkuWwNFHwx57wD33wIIFPlQuWAB33+23H320309ERERqSiFT0sE5OPZYeOABHyo7Ojq/3tEB8+fD/ff7/WrZoqnWVhEREYVMSYmJE2HsWB8wi2lr8/tNmpRMvXKptVVERARQyJS0uOwyHyBL0dbm909amlpbRUREYqaQKenw0ENdQts1fBfDMYdBnfft6PD7Jy0tra0iIiIJUMiUdMjTinkFZwIwg2El7R+7Aq2t77N6131r1doqIiKSEIVMSYf+/btsaqcXAM0sLWn/2OVpbb2dI1iT92llt8771qq1VUREJCEKmZIO++7rR2lnWUozAL1o77xvU5PfP2l5WjGfZScAXmXzkvYXERHpKRQyJR1Gj+7SOhm2ZDqs8779+vn9k5an9XQeAwEw8gzyqUVrq4iISEIUMiUdtt8e9t+/UzALWzLDsAn41w84ALbbLuka5m1tHcPxQJ6QWavWVhERkYQoZEo6mMFNN8GBBy6b4LxTyGxqggED/Os33eT3T1qe1tZQEznTGdWqtVVERCQhCpmSHr17w623wrhxcOihzGIFANr7LwejRkFrK9x2m9+vFvK0toY6hcxatraKiIgkRCFT0sXMh7k771y26Y8nvQB33FH70JbT2vqZLZ95CVcfra0iIiIJUciUVFq8OPP4yitrV48uslpb9x46ftlm69OnPlpbRUREEtJc6wqIVGLOnFrXoIigtXXirKxNV18FJ9auSiIiIklTS6akUu4Uk//5T23qUSr1jIuISKNRyJRUcjkzAm26aW3qUaom/UsTEZEGo199kko5qzfWPYVMERFpNPrVJ6mU25IpIiIi9UUhU1IpbS2Z7e3d7yMiItKTKGRKKqklU0REpL4pZEoqpa0lc8iQWtdAREQkWQqZkkppC5lqeRURkUajkCmpVO+h7d13Oz9PWygWERGplkKmpFK9h7aTT+78vN7rKyIiEjWFTEmlem/JfOSRzs8VMkVEpNEoZEoqpS20ffObta6BiIhIshQyJZXSEjLDUeVLltS0GiIiIolTyJRUqvfu8lCfPrWugYiISG0oZEoqpaUlUyFTREQalUKmpFJaWjJ79651DURERGpDIVNSKS0tmf3717oGIiIitaGQKak0YUKta1CagQNrXQMREZHaUMiUVDr9dF8+9lht69EdhUwREWlUCpmSamaZx/XShZ5dj5/+tHb1EBERqSWFTEm1pqxPcL20ara3Zx5vvnnt6iEiIlJLCpmSatktmfUy4fnSpZnHZjB8eM2qIiIiUjMKmZJq9RIss2WHzKYmOOII6Nu3dvURERGpBYVMSbWFCzOPm+rk0zxnTuZxU5OfK3Px4vTM7SkiIhKFOvm1LFKZ7JCZ3XVeS9/6VuaxmV/1x7nO92qKiIj0dAqZkmrZIbNevPFG5nFTU2ZpyXrs2hcREYmLQqakWj22ZGZ3i5tllpZcvLg29REREakFhUxJtVVXzTyux5CZ3ZKpkCkiIo1EIVNSaYcdfHnAAbWtRz6FQqa6y0VEpJEoZEoqjR/fdVu9jC5Xd7mIiIhCpvQg6i4XERGpHwqZkjqF1iivl7XLc1sy1V0uIiKNSCFTUid3vsmzz86/vVbUXS4iIqKQKSkULtu4006+POIIX9ZjyAR1l4uISGNSyJTUCUPmwQf7slcvX9ZLyNx2287P1V0uIiKNSCFTUicMmc3NvgxDZr3ck7nLLp2fq7tcREQakUKmpE5uyAynLqqXlszcsKvuchERaUQKmZI6hVoy6y1k7rGHL9VdLiIijUghU1Kn3kNmOPDnscd8qe5yERFpRAqZkjr1HjLDlsxwcnh1l4uISCNSyJTUqfeQGYbJ8F5RdZeLiEgjUsiU1Kn30eW/+IUvw5ZMdZeLiEgjUsiU1Kn30eW51F0uIiKNSCFTUqdQS+Y119SmPt0ZPNi3as6cWeuaiIiIJEchU1KnUMh89dXa1Kc7/frBsGHw/vu1romIiEhyFDIldQqFzHrWvz8sWlTrWoiIiCSn5iHTzFYws2+b2b1mNsXM2sxsjpk9Y2bfMrO8dTSzHc3sYTObFRzzqpmdYWYFI4eZ7WdmrcH555nZBDM7rpv6HWdmE4P95wTH71ft9y2VS2PI7NsXFi6sdS1ERESSU/OQCRwGXA98EZgA/Ba4G9gU+BNwp1k4TtczswOBp4BdgXuBq4E+wBXA7fnexMxOBcYG5705eM/VgDFmdmmBYy4FxgCrBvvfDGwGjA3OJzVQaOBPPevXTy2ZIiLSWJprXQHgTeAA4CHn3LJJaMzsPGAicChwCD54YmaD8IGvHRjpnHs+2P5TYBwwysyOdM7dnnWu4cClwCxgW+fctGD7z4FJwGgzu9s596+sY3YERgNvAds55z4Ltl8CvABcamYPhueS5KglU0REpP7VvA3IOTfOOTc2O2AG26cD1wZPR2a9NApYCbg9DJjB/guBnwRPv5vzNicAfYGrs0NhEBx/GTw9OeeY8PkvwoAZHDMN+H1wvuO7/w4lamkMmWrJFBGRRlPzkNmNcI2UpVnb9gjKR/Ps/xSwANjRzPqWeMwjOftUc4wkII0hs08feOaZzLrmIiIiPV3dhkwzawaODZ5mB70Ng/LN3GOcc0uBqfjbANYp8ZiPgPnAGmY2IHjvFmB1YF7weq7JQblBSd+MRCoMmWG4TMM9mY8/7ieLv+GGWtdEREQkGfVwT2YhF+MH6TzsnHssa/vgoJxT4Lhw+5Ayj2kJ9ltQ4Xt0YmYnAScBDBs2jNbW1kK7RmLevHmxv0e9eP75YcDGvPbaBGbPbgNglVW+SJ8+HbS2Tor0vSq7riMBco7z25544h3WW29q9RXrARrpM5skXdf46NrGQ9c1HvVwXesyZJrZ9/GDbv4HfKPG1amIc+464DqAbbfd1o0cOTLW92ttbSXu96gXr73my732+iIrreQfb7cdvPcekV+DSq7rkCHwjW/kr8tmm63NyJFrR1K3tGukz2ySdF3jo2sbD13XeNTDda27jsZgaqArgf8CuzvnZuXsErYiDia/cPvsCo6Zk1OW8x6SkNmzfTk466fT3JzpRq+ljg5fv0L3iTbX5Z91IiIi0aurkGlmZwBXAa/hA+b0PLu9EZRd7ocM7uMcgR8o9HaJx6yK7yp/3zm3AMA5Nx/4ABgYvJ5r/aDsco+nxG/2bBgwwA+mCTU3+3sea+388335+ef5X+8846uIiEjPVTch08x+hJ9M/WV8wPy4wK7jgnLvPK/tCgwAnnPOZU8YU+yYr+XsU80xkoDZs2Ho0M7bPv8cXn8dlizJe0hibrvNl3MK3c0rIiLSIOoiZAYTqV+Mn+R8T+fcJ0V2vwv4BDjSzLbNOkc/4KLg6R9yjrkRWAScGkzMHh4zFDgveHptzjHh8x8H+4XHDAdOCc53Y3ffm0Rv7lxYbrnO2x4N5h+4Pe96T8kr1F2ehpHwIiIiUaj5HWLB2uE/x6/g8zTwfevapzjNOTcGwDk318xOxIfNVjO7Hb+SzwH4qYruAu7IPtg5N9XMzgZ+BzxvZncAi/ETu68BXJa92k9wzHNmdjlwFvCqmd2FX7ryCGB54DSt9lMb7e2F722s9TyU4Ue3UJjs6Mi/XUREpKepecjE30MJ0As4o8A+T+LXEAfAOXefme0G/Bi/7GQ/YAo+EP7Oua5Rwzl3lZlNA36An3+zCT+46CfOub/ke1Pn3Ggz+ze+5fIkoAN4EbjEOfdgWd+lRKajo3CIq/XE7GGILFSPWnfni4iIJKXmIdM5dwFwQQXHPQvsU+YxY4GxZR4zhqyAK7WXxpD54x/DL35RH4OTREREkqA7xCR1ioXMWt/zGIbI3JB54YWdXxcREenpFDIldeo5ZOYueRkK66WQKSIijUIhU1Injd3lZr7OCpkiItIoFDIldeo5ZIYhMl/9evVSyBQRkcahkCmpU8/d5cVGlytkiohII1HIlNRJQ0umQqaIiDQ6hUxJnWIhs9a6C5majF1ERBpFnf6qFiksX8gcNMiX9dKS2XXRKg38ERGRxqKQKamTL2TeESwk2rt38vXJpoE/IiIinkKmpE6+kDlggC9rvXa5QqaIiIinkCmpky9kht3TtQ6Z4T2X+brLFTJFRKSRKGRK6hQLmfUysEYhU0REGp1CpqROvpAZPq91S2ZIIVNERBqdQqakTj13l4d0T6aIiDQ6hUxJnRdfhHnzOm+rt5CplkwREWl0CpmSKkuX+rK1tfP2ersnUy2ZIiLS6BQyJVXa2vJvr7d7MhUyRUSk0SlkSqoUCplp6C5vaqqfllYREZG4KWRKqoQhc6+9Om8PQ91f/pJsfXJ95Su+POywrq+pJVNERBqJQqakShgyv/nNztvDkHnXXYlWp4uNNoIhQ2DDDbu+ppApIiKNRCFTUiUMmf37d96e7x7IWsg3vVJIIVNERBpJnfxqFilNGDL79eu8Pd89kLXgXOG6KGSKiEgjUciUVCnUklkvA2qcU0umiIgIKGRKyhQKmfUS3jo61JIpIiIC0FzrCoiUI80tme3t8PHHydZHRESkVhQyJVXqPWRed13h1555Jrl6iIiI1Jq6yyVV6r27vBT1EohFRETipJApqVLvLZmlWLSo1jUQERGJn0KmpEo9h8xSW1MVMkVEpBEoZEqqtLX50dt9+3beXg8hc8mS4q///ve+rGnIdA4mTPDrXra0+FFKLS1w+OEwcWL9LP4uIiKpp5ApqbJoEfTp03WaoHq4J3Px4uKvh8G4ZiFzyRI4+mjYYw+45x5YsMCHygUL4O67/fajj+4+LYuIiJRAIVNSpaPDzzeZa5ttkq9LrjBkHnxw/tdrGjKdg2OPhQce8KEyt+m3owPmz4f77/f7qUVTRESqpJApqdLenn8eyqFD4aCDYPPNE6/SMmED4F575X+9piFz4kQYO9YHzGLa2vx+kyYlUy8REemxFDIlVQq1ZILfXst7M8OWzN69878ehszuutVjcdllmVFTgZmsSBv9uu7b1ub3FxERqYJCpqRKoZZM8NvrIWT26ZP/9Zq2ZD70UJeLszIz2ZWnuu7b0eH3FxERqYJCpqRKsZbMpqbaDgDqLmSG22sSMnNaMUPPs11Z+4uIiJRKIVNSJQ0tmbnTK4Vq2pKZO7FolrxDfIrsLyIiUgqFTEmVNNyTWZfd5fvuWzCdv8IWnTc0Nfn9RUREqqCQKamShpbM7rrLazLwZ/Togq2TbeRs79fP7y8iIlIFhUxJlTTfkxmG45oE4e23h/33zxs028m6oP37wwEHwHYF7tUUEREpkUKmpEqxlsw+fWrUShio65BpBjfdBAceCC0tOMtcxJ9zvq/cgAH+9Ztu6rqkkoiISJkUMiVVirVktrT4RWtqpdSQWbPFdHr3hltvhXHjaD/ksGWbH+erMGoUtLbCbbcVnuhTRESkDM21roBIOYq1ZA4Y4EOmc7VpiCs1ZL71VjL1ycsMtt+eu0bdDndnbb/jjppVSUREeia1ZEqqdNeS2dFRo9HbZN63u5D5058mU59inn++1jUQEZGeTiFTUqVYS2a/YIXEhQuTq0+2Ulsy64FuuRQRkbjV0a89ke61txduyWxuzuxTC92FzHoKdvVUFxER6ZkUMiVVOjoKtwiG4XPp0uTqk627FX/qqSVz0KDOz9V9LiIiUaujX3si3SvWkhlur9eWzHoKmWus0fn5U0/Vph4iItJz1dGvPZHuLV2a6RbPpZBZmrlz4fjjO2/TAj8iIhK1Ovm1J1KaxYsLd0eH4bPW3eWFppmsl5D573/XugYiItII6uTXnkhpFi0q3FJYDy2ZvXsXHlRTL4NtCl0/ERGRKClkSqosXlzfIbNYgKuXlszsltZ6mLNTRER6pjr5tSdSmlK6y2sVMou1skL9hMzseuy5Z+3qISIiPVud/NoTKU0pLZm1vCczDSEzO4T371+7eoiISM9WJ7/2REqzcKG6y6uVHcILDVISERGpVp382hPpXlsbTJsGa6+d//Vad5d3FzLrZeBPdsjUICAREYmLQqakxqJFvlxxxfyv17q7vK2tePdzvbRknnVW5rFaMkVEJC518mtPpHtheCw0GfvAgb6cMiWZ+uRKS8gcPz7zWCFTRETiUie/9kS6113I3GADX86dm0x9ci1cCP36FX69XkJm6A9/6Nxd7lzt6iIiIj1Pnf3aEyksDJmF1i4PW+XS0F1eq/tGs225ZeeWzHqok4iI9BwKmZIaYQgq1JIZbl+yJJn65ConZC5cGH99utPc3Lkls1bhXEREeiaFTEmN7rrLa92SuXBh8ZCZPbq8VnXs6Mg8bm7u3JKpkCkiIlFSyJTU6K67PAyfaegur1Vr62WXZR4rZIqISJwUMiU1umvJrIcpjIoN/MluyVy8OP765PPss5nHvXr57vIdd/TPaxV8RUSkZ1LIlNTo7p7MsKXw5ZcTqU4X3XWXZ6tVyMweQR5ex298w5dqyRQRkSgpZEpqdNeSGbrvvtirkld3Uxhlq6eQWet7WUVEpGdSyJTU6O6ezFpyztev1MnNa9U1nTvwJ7tUyBQRkSgpZEpqlNqSWQtheCs1AH/+eXx1KSa7JTOsq0KmiIjEoS5CppmNMrOrzOxpM5trZs7Mbi6w7/Dg9UJftxd5n+PMbKKZzTOzOWbWamb7Fdm/l5mdaWavmlmbmc0ys4fNbMcovm8pT7h2ed++ta1HPuUG4GnTYqtKUflaMsN7XT/4IPn6iIhIz1UvbUI/AbYA5gHvAxuVcMwrwH15tr+Wb2czuxQYHZz/eqAPcCQw1sxOc85dnbO/AbcDo4A3gKuB5YEjgKfM7FDn3P0l1FMiEt7HmD2BeK4DDoB3302mPtnCoFZqS+b8+fHVpZjskBnWNZwYftw4GDky8SqJiEgPVS8h80x8+JsC7AY8UcIxLzvnLijl5EHL42jgLWA759xnwfZLgBeAS83sQefctKzDjsQHzOeAPZ1zC4NjrgWeAa43s3HOuRp1fDaeCy7wZbGWzObm2nT7ltuSWasVf7JDZjga/2tf8+WqqyZfHxER6bnqorvcOfeEc26yc9l3jEXq5KD8RRgwg/edBvwe6Ascn3PMd4PyJ2HADI6ZBNwBrIQPoZKQ8eN9Wawls1ev2oTMUlsyn3nGl/UQMsN5O8Npl773veTrIyIiPVddhMwKrWZm3zGz84Jy8yL77hGUj+Z57ZGcfTCzfsCOwALg6VKOkeQ0FfnUNjdnAl+SSh35vu22vgzvL01a9p9xLS2+rMd7XEVEJP3qpbu8El8JvpYxs1bgOOfcu1nbWoDVgXnOuY/ynGdyUG6QtW1doBfwtnMuX7tYvmMkIcXmmKxVd3l3E8WHwlbYWrVkhiFzl10y0y2VOreniIhIOdIYMhcA/4cf9PN2sG1z4AJgd+CfZralcy4cWjE4KOcUOF+4fUjWtkqO6cTMTgJOAhg2bBitra2Fdo3EvHnzYn+PWhswYGcWLGjm/fefYc6c/Ely5swNmT9/KK2t4yN5z1Kv68yZfYEvMWXKG7S25vtbJqN3712ZPPl9WlvfLrpfHD77bEtgCHPnzqa19WUgDJ4jARL9DDXCZ7YWdF3jo2sbD13XeNTFdXXO1dUX/redA24u87hmYHxw7OlZ21cLtr1f4LjeweuLsrbtGGx7psAx6wevv1FK3bbZZhsXtyeeeCL296i1E090bpVViu/z7W87t9pq0b1nqde1tdU5cO6qq7rft6XFubPOqq5eldplF1/PkSM7b/dRM9m6NMJnthZ0XeOjaxsPXdd4VHNdgeddBJkuzfdkduJ8t/afgqe7Zr0UtjoOJr9w++wqj5GYOVf8fkyoXXf5WWf58pFHiu8H/r7NWtw3CjBzpi+7u44iIiLV6mm/aoJfobSEG5zvNv8AGGhm+SZpWT8o38za9hbQDqxjZvluKch3jMSsoyMzIrqQWo0uL2dehFqGzP/9z5dbbFGb9xcRkcbR00LmDkGZe7PbuKDcO88xX8vZB+enLHoOGADsUsoxEr+OjvptyQznmNxmm+73rWXIDP3sZ7V9fxER6flSFzLNbGsz61JvM9sTP6k7QO6SlNcG5Y/NbGjWMcOBU4BFwI05x/whKC8KpjQKj9kOv+rPTODuCr8NqUCp3eW1CHBbb+3LUsJbrYJwtnBkuYiISFzqYnS5mR0EHBQ8XSUov2RmY4LHnzjnfhA8vhxY38yew68SBH50eThn5U+dc89ln98595yZXQ6cBbxqZnfhl5U8Ar9U5Gmu82o/4JeUPAQ/4fpLZjYWWCE4phdwonNubqXfs5SvlO7yWq7406dPactK1kNLZnfXUUREpFp1ETKBLYHjcratE3wBvAOEIfOvwMHAdvhu697ADOBO4GrnXL7J03HOjTazf+NbLk8COoAXgUuccw/m2d+Z2VH4bvMTgNOAhcBTwEW5QVbiV0pLZq3uyVy6tPR1yz/4AP78Z7jhhnjrVEzudTzkELjnntrURUREeqa6CJnOr0F+QYn73gBU9OvZOTcGGFPG/kuBK4IvqbFS78lsb/eBNMnWuvb20tctr5X58zOPc6/juutqUnYREYlW6u7JlMZVand5uG+Sli6t/5CZPQI+N2T26pX8NRMRkZ5NIVNSo9Tucki+yzwNITP7muRex6YmhUwREYmWQqakRqnd5ZD8wJr29tLvyayVJUsyj3NbhBUyRUQkagqZkhrldJerJbOrYtdEIVNERKKmkCmpoe7y6nQXMgE++yyZuoiISM+nkCmpUU53edIhMw2jy4tdkzCc77hjMnUREZGeL5KQaWYbmdmZZvYdMxscxTlFcpXSXd6njy8XLYq/PtnKmSfz5JMz9UxSGDJvzl0Pi0x4D9c2FxERqVZZIdPMzjezj8xs+axtXwZeAi4FrgFeNLMVoq2mSGnd5QMH+jJ7TsgklNNd3tJSm2Udw5CZr57PPptsXUREpOcrtyXza8D/nHOzsrb9CnDAz/DrfY8ATo+meiIZpXSXt7T4MumQuWRJ6cGxVoNsioXMV19Nti4iItLzlRsyhwOvh0/MbHVgG+Aa59xFzrlTgXFk1iEXiUwp3eW1CpmLFkHfvqXtW48h849/TLYuIiLS85UbMocC2a2YO+FbMbPX/n4BWKvKeol0UUp3eS1DZqn3WdYqZIbzZOYLmXvs4ctNN02uPiIi0rOVGzJnAqtnPd8dWAJMyNrWp4LzinSrnrvL096S2a8fbLghfOELydZJRER6rnLD4MvAAWa2qZmtBxwBPOOca8vaZzjwUTTVE8kop7t83rz465Nt0qTSVxmqdcgsdO9o797JT/0kIiI9V7kh8zfAYOAV4I3g8WXhi2bWC9+F/nxUFZQacA4mTIDDDvOpranJl4cfDhMn+tdrVK16bMlcssSHsyeeKG3/l17y38t//hNvvXIVa8kMtytkiohIVMoKmc65p4H9gPuAe4FRzrlHsnbZEfggeE3SaMkSOPpof5PePffAggU+ES1YAHff7bcffXTnhbATUspclP37+3LhwvjrE1q8uLz9HwzuYL434X8l3YXMl1+GBx5IrDoiItLDlb1GiXPuUeDRAq89DWxVbaWkRpyDY4/1SWPBgq6vd3T4JsL77/f73Xpr9/3XEVq0CJZbrvg+YVdwki1yYTf5OeeUtn+vXv6YJIMwZG4hqPeViUREpGfQAB3JmDgRxo7NHzCztbX5/SZNSqZegRkzuh/BHQaoJBtaw5A5bFhp+4cDhJIOmYce6kuFTBERSULFv27MbA38SPO8Y2qdc09Vem6pkcsu8wGyFG1tfv877oi3ToFPPoGpU/1XMbVYu7y7buhctVpfPff9c511lubLFBGR6JQdMs3sq8AVwEbd7FriSs5SNx56qNOw5zkMYiSt3MrRbEzOotYdHX7/hMyZU9p+Zr47uhYtmaWuXR6Om+puEFNcCoXM3r2Tn/pJRER6rnLXLt8BP/H6EOBqwICngOuB/wXPxwI/j7SWkoycVsxH2ZuX2YoLuKCk/eNUaoADH/qmTImvLvneD0pvyVxvPV/26xdPfbpTqAX1ySd9+fe/J1cXERHpucptSzkXWAhs55wL1yd/wjl3MrApcBHwZeCu6KooiQmHZgcM3+TmKDC4J2f/OJXbtXznnfHUI5+wbqUG4Xvu8eUKK8RTn+4UmoFq/HhfPvNMcnUREZGeq9yQ+SXgAefch7nncN75+LXNL4yofpKkffft1IdbNGQ2Nfn9ExJ2f998c2JvWbJyu8tXXrnzcUlpaYGhQ2GLLfK/HtY/6XqJiEjPVG7IHAy8m/V8MdCSs8+zwK7VVEpqZPToTq2TYcjMq18/v39CwpBZqy7mYsrtLg/DXJIDf2bM8PdbjhxZeJ+w/gqZIiIShXJD5sfA0Jzn6+bs0xtIrh9VorP99rD//l26wbu0ZPbvDwccANttl1jVyh3BnaRyu8trEeZuucWXxSaAV0umiIhEqdyQ+SadQ+V44CtmtgGAma0CHApMjqZ6kigzuOkmOPBAaGnBgonWO8KPSVMTDBjgX7/ppkQnYg9bMgutu51t4MB465Kr3O7y8I6EJFsyBwzofp+7gjup11473rqIiEhjKDdkPgrsZmbLB8+vxLdavmRmk/AjzFcCfhtZDSVZvXv7lXzGjcN2+CIAjmDt8lGjoLUVbruttLQXoTCQlfK2J52UWcM8CeV2l4fTLCXZYhgOMgpHtufzpS/5UuuXi4hIFMoNmX/E32+5BMA59yxwGDAVP7r8I+C7zrmboqykJMwMtt8eO/sHALy6zkF89t48P/F6gl3k2cKWzFKCXO/etZmMvZxplnr1SraOYci85prC+4T3u44ZE3t1RESkAZR1h5tzbi4wIWfbvUCRO70kre6+25dvvw277Qavvlq7upTTXZ50gAvr1t2Sl9mam5NtyQzn2C82cCpc7vKVV+Kvj4iI9HzlTsa+q5mt1c0+a5qZRpf3AOFgEYB//7t29YDyBv6EAa7QfJBRW7TIl+WEzKSDcCmrDCV4i62IiDSAcrvLnwC+2c0+xwb7SYq9+GKta9BZOS2ZYRDNWiEzVosX+zJsCSxFrVoya7WUpYiINJ5yf+WU0tZhUGyCRUmD6dNrXYPOyhn4E4bMpFoKw5BZzy2ZYchUa6WIiCQljnaNtYHPYzivJCjBFSNLUs7An6QnOw8DeT3fk1lKdzn4EebDhsVfHxER6fm6/ZVtZufnbBpp+ZtDegFrAUcCWv045cJQVy8q6S5PKsR9+9u+LGdWp3ptyVx3Xb86kIiISLVKGV1+QdZjB4wMvgr5ADin4hpJXQi7gOtFuQN/so9JSjkDjZJuyVy40JfdtWQ2N2ueTBERiUYpIXP3oDRgHDAG+Eue/dqBT4E3nHMJDbmQuKS5JTPp7vKdd4ZnnoFNNin9mKRbMg87zJfdtWQ2N9ffz15ERNKp25DpnHsyfGxmfwHuy94mPVO+lsz29vImHI9SOEF4PXaX9+0LO+5Y3jFJt2SGumvJTHoiexER6bnKGvjjnDveOfdAXJWR+pEvZP7pT8nXI/Tcc74sZXBN0t3llYTvpFsyQ+ouFxGRpJS14k82MxsADMUP+OnCOfdupeeW2svXZTpnTvL1CB1yCNxzDwwd2v2+SXeXVxIya9WSqe5yERFJStkh08y+AfwI2LjIbq6Sc0v9yNeSmdQKOvkMHAjDh5e2b9Ld5R0d5Y0sh2RbMmfNyjxWS6aIiCSlrCBoZt8E/owf5PM08B6gX0k9UL7WrFqGzPb20lerSUN3eZItmdOmZR5315KpezJFRCQq5bY2/gD4DNjZOfd6DPWROhG2ZB5+ONx5p39cy5DZ0VF6kEtDd3m935PpnFYHEhGR6pS74s96wN8UMHu+sCWzX7/MNrVk5tfRUf6a4M3NMGFCPPXJlR0WS7knE2pzv6iIiPQs5YbMWcCiOCoi9SVsycxuoUtLS2bSQamSlsxZs5IbYFNOyAzvLVWXuYiIVKvckPkgfllJdaT1cIsX++D0xS9mtqWlJTMN3eVf+1py3dHZ79Pdz7BWqyWJiEjPU27IPBfoC1xrZgNjqI/UienTfXg66aRa18Qrp0s6Dd3lffokt3RnJSFT0xiJiEi1yh348zdgAfBt4GgzmwzMzrOfc87tWWXdpIZuuMGX5QSUOPW07vI+fXyQqySgVqO7n6G6y0VEJCrlhsyRWY9bgC0L7FfDOCJx6ajhivSVDPyZOhV23TW+OoUqCZl33eXLa6+F730v+jplyw7b6i4XEZGklLusZFOJXzVa4VrilJaWzPnzffnNb8ZWnWWWLoXJk8sPZeHclW+9FXmVuignZIYtmeE1FBERqVSCHXWSdmkZ+DNvXrx1yXbddb4MWybrUSUtmfvsE199RESkMShkSrfC+zJr2V1eTkvmwASHpFUaaE8+2ZerrRZdXQrJDpnd/Qzb2nw5eXJ89RERkcZQdsg0syYzO83MxpvZHDNbmvXaVmZ2jZltEG01JUlhEDn8cF/Ww2CQcloy997bl2GQi1O592KGRo/2ZRKBuJyWzI02ircuIiLSOMoKmWbWB3gc+C2wLvA5kD3b31TgBOCYiOonNRBOX7PFFr4MQ2Ytp7UppyXTzLcQJhGKKx0ZHq6kFLYcxikMmeuuC5tuWnzf3XbzdTvooNirJSIiPVy5vyLPBnYHLgSGAX/KftE5Nxt4CtgrispJbYTzN/bp48t6mDuxnJZM8EEpiQBXaUtmktMshe/x5z+XVt8vfEGjy0VEpHrlhsxjgGedcz93znWQf6qiqcBaVddMaiYMk2EL5uOP+3KVVWpTHyivJRNg+eX90o1xqzRkhoE5iftcw8BYzjyjCpkiIlKtckPmCGB8N/vMApavrDpSD8KQGbZkbrmlL9MyTybAsGF+1aK4VTriPgx8SbRkLlzoywEDSttfIVNERKJQbshcCAzpZp+1yL8KkKRE2F0etmQmvRZ4PkuWZOpTigEDMuEqTpXeQpBkyFywwJflhEwtKykiItUqd8Wfl4Gvmlkf51yXlZfNbDD+fsznIqib1Ehud3lTkx9Mk9QyjbmWLoWXXiqv1bB372TWBl9pJV/+4hflHReGzCRah8OQ2b9/afv37p1MQBcRkZ6t3JbM64A1gVvMbFD2C2Y2BBgDDAWujaJyUhu5A3/Ah6JatWT+61++fPnl0o8J1waPWxgWDzmkvOPCrv96bclUd7mIiFSr3GUlb8MHyUOBmcB3AczseeAj4EDgGufcw9FWU5L0/vu+zO6eXroUfvWr2tSnkmmCevdOJmSGLZHl1tEsudbhcJS9ustFRCRJZf/6ds6dgJ8L87/ASvh5MrcGpgDfcs6dFmkNJXFHHunLeln1pZ5DZhgSKxll3qtXst3l4dyc3VFLpoiIRKHcezIBcM6NAcaYWX989/gc59z8KCsmtfPZZ75Mcg3wYioJcEndk1lNyFy6FF57Ldr65LNggQ+YpYb13r0VMkVEpHpVrV3unGtzzn2ogNmzDBvmy759a1uPUCUtmUndk1lNyAR44IHo6pLPBx/AJZeUN5BHLZkiIhKFqkKm9Eynntq5BNhzT79UY1qkobs8CW+8Uf4xNb0n0zmYMAEOOwxaWvxfGC0tcPjhMHFi5ROTiohI4srqLjezt0vc1Tnn1q2gPlIHwpbD7Clv1lqrssAShUoGx4Rdvs75ATZxqfeQmT1DQKlq1l2+ZAkce6xv3l24MHPD6oIFcPfd8PDDsP/+cNNN5U2aKiIiNVFuS2YTfqBP7tdQYHjw1aeC80odyRec+vWDRYtqU58w8Bx8cOnHhBkk7ha5nhgyP/oI3nsv4Z+3c5mAuWBB1xFRHR0wfz7cf7/fTy2aIiJ1r9wpjIY750bk+RoKbAA8CrwFbBxHZSUZ4e/37ODUt2/tJugOg1x29313FDK9SkLmo4/6MlyzPhETJ8LYsZmh8IW0tfn9Jk1Kpl4iIlKxyFocnXNTgEOA1YGfRXVeSV4YnLIH3NRDS2Y5QS4MV40eMqvRXNHcExW67LLMhJ7daWvz+4uISF2LtFvbObcQeBw4KsrzSrLa2zOThYc++MBPCfS//9WmPlBe6FFLphe2Sl96afnHJhoyH3qoSxf5fAZwMT9iKTkXt6PD7y8iInUtjnsnlwKrlLqzmY0ys6vM7Gkzm2tmzsxu7uaYHc3sYTObZWZtZvaqmZ1hZgV/1ZvZfmbWamZzzGyemU0ws+O6eZ/jzGxisP+c4Pj9Sv3e0qq9vWtouvNOX954Y/L1qaQlMwyZcc+VGUXIfPHFaOqST1i/DTcs/9hEg3OeVsyfcz7ncjE38/WS9hcRkfoSacg0sxWBg4H3yjjsJ8CpwJbAByW8x4HAU8CuwL3A1fjBRlcAtxc45lRgLLApcDNwPbAafkL5vG08wfYxwKrB/jcDmwFjg/P1WB9/3HV0cS1b7MI8Uc68nWlqydxjj2jqkk8l9dtpJ1/OmBF9fQrKnsogcAdHAPA5y5W0v4iI1JeyQqaZnV/g6+dmdiPwJplQVqoz8YOGBhGshV7k/QcF524HRjrnvuWcOxsfUP8FjDKzI3OOGQ5cCswCtnXOneKcOxPYHD9IabSZfSnnmB2B0cHrmzvnznTOnQJsE5zn0uC8PdINN3TdtuWWvhw6NNGqAJmws0rJ7ePpuCfznHN8OWdOdPXJVcna6r/8pS+PSvKml333XVbJDowX2Jp3GA7Ae6zZed+mJr+/iIjUtXJbMi8o8PUT4LjgfBc5535T6gmdc0845yY7V9KcJKPw66Xf7px7PuscC4M6QNegegLQF7jaOTct65jPgODXKSfnHBM+/0WwX3jMNOD3wfmOL6G+Pca11/pyrbWSf+9Zs3y5wgqlH5N0S2YlqxKtv360dcmnkhBckwFeo0cva53sRQfb8sKyly7hh5337dfP7y8iInWt3F+Nuxf42g3fMriicy7OkeVhx+KjeV57ClgA7Ghm2R2rxY55JGefao7p0QYN8uV99yX/3mHoKWf+7STvyay0qzyJqR4racncc09fbr119PUpaPvt/UTr3XWD9+8PBxwA222XTL1ERKRi5c6T+WSBr6edc6855+JeJyQcvvBmnrotBabiVzFap8RjPgLmA2uY2QAAM2vBT8M0L3g91+Sg3KCSbyCtwtB2e967XuO1aJFvvCpn5Z6kWjIXLarvkFlJS2ZTk5/4Pu6A3omZX8nnwAMLV2rAAP/6TTfFu4yTiIhEotuQaWZ9ghHW/zSzgm1JwX7/NLPxxfar0uCgLHQXW7h9SAXHDM4py3mPHq+Wq/gtWlTeoB9IJmSOH++na6w0jA0cGG198sk3sX4pBg2K917RvHr3xt1ya96X3KGjoLUVbrtNS0qKiKREKTPhfR0/4GV/51zBX9nOucVmdgnwMHAMfmR2wzKzk4CTAIYNG0Zra2us7zdv3rwI32MkQKfzffppH2DHLtuTMHXq+jQ1rURr63MlH/Pf/w4BtmTixJdYvLjytFTsut5332qEDdqVXJOVVwYYyV57Tae1NZ4JSG+5ZW1gBC+//CJLl84t+bi5c9fjs8+G0dr6bCz1gvzXdvFiw99909kP176QfedP90FTior2/wLJpmsbD13XeNTFdXXOFf0CHgQmd7df1v5vAA+Vun/OsSMBB9xc4PVJwevbFHj9teD1jbO2zQy2rVDgmHnB6wOC5y3B888L7L9i8PqMUr+vbbbZxsXtiSeeiOxc667r3DHHdN72ySfO+Q7eyN6mZN/8pnNrrlneMU8/7ev6979X997FruuNN1Z/TVZbzblvfavy47sT1m/8+PKOO/dc53r1cq6jI556OZf/2ra0ZOqc/fWTn8RXj54myv8LpDNd23jousajmusKPO8qyHG5X6Xck7kV0FpGbn0KP6VQHN4Iyi73Q5pZMzACPxn82yUesyo+VL7vnFsA4Jybj5+vc2Dweq5wTHCXezx7iqVLu672Eo7s/trXkq/PrFmZgUelSqK7PIo5Q3v3jv++USh/9Ptyy/n7OZNer37+/Pzbk66HiIhUr5RfPSsC5UzLPAMoY7KZsowLyr3zvLYrMAB4zjmXPQlLsWO+lrNPNcf0GPlCJsBmm5V/b2QU3nkHRowo75gkQmYl0xbl6tMnmQBVbiBeLpj//PPPo69LJRQyRUTSp5Rfk21AOUMUBgJx/Uq4C/gEONLMtg03mlk/4KLg6R9yjrkRWAScmj2BupkNBc4Lnl6bc0z4/MfBfuExw4FTgvPVYIHFZBQKmb16ZUYrJ6m9PTO5eqnC/eMcIZ2z1HZF1lrLh+i4lXv96i1k5q5AJSIi9a+UgT/vAdt2u1fGtsC7pe5sZgcBBwVPwzVdvmRmY4LHnzjnfgDgnJtrZifiw2armd2OX4HnAPxURXcBd2Sf3zk31czOBn4HPG9mdwCL8RO7rwFc5pz7V84xz5nZ5cBZwKtmdhd+6cojgOWB01zWxO49TT2GzHJbDZNY8SeKc6+8Mrxb8r+W8mSH4CFDyju23kJmLT53IiJSnVJ+dbfiQ1+3QdPMtsEPQX6ijDpsiV8t6Dhgr2DbOlnbRmXv7Jy7Dz/89CngUOA0YAk+EB4Z3LBKzjFX4YPof4Bj8aO+pwPfDANsnmNG41f1mR7sf2xw/P7OuavL+P5Sp95CZkdH5SHzH/+Ivj6hMGReeGHl52hpgXnzoqlPruyfVbn3tIb7zy19QHqs8i11KiIi9a2UX91X40dT/83MNi60k5ltBPwNv674NaVWwDl3gXPOinwNz3PMs865fZxzQ51z/Z1zmznnrnDOFYxAzrmxzrndnHPLOedanHPbOef+0k3dxgT7tQTH7eace7DU7y2t6jFklntPYRgyb4zxpoYwZH43dyHTMgwcWHiwS7Wyf1Zhy2SpBgezxSY+V2Ygd2WpKG5NEBGRZHXbXe6ce8PMfo5fo/yloOt4HPB+sMvqwJ74VsW+wPnOuTfynUvSoVDIbG5OX0tmnMKQWc3c4GFLpnPRL2IT/qwuvrj8c4chM+mWzLXX9mu6H3igvzZxBXAREYlfKfdk4pz7uZktBX4GHA0clbOL4busf+yc+1W0VZSkFWvJrMUAjJ4cMgcO9N9fuHRmlCpd7QcyS4i3tUVXn1I0N4eT1GvlSBGRtCspZAI4535pZrcAJwA7AeEckh8BzwA3OucSGCcrcXLOt4DlC5mvvQaffJJ8naoZ+BOnMHBXGzLBt2ZGHTIrWbc8FNZl0aLi+0XppZfgrbdgR7+wFGuvDf/5T+b1Qn/8iIhIfSrrv+wgRP4sprpIHQiDSb5f5rUImFBZS2YSy1uHLZnVBJ+WFl/Omwcrrlh9nbJVEzLD+VCTnJ9y6619GV7Pxx/3q0gefbR/vnixQqaISJpEMJ209CRh61y+X+bhaj9JD8KoZOBPFKvxdGfJEh9+q5mUPWzJjOPew7S1ZIbC+q66Khx1FFxxhX/+0UfJ10VERCqnkCmdFAuZu+7qy6SDRyUtmQDbbZdpHYvDkiXVt5hmt2RGrZqQ2bu3vyeyFivtTJ/e+flbb/nyrLOSr4uIiFROIVM6KRYya9W6Vck9mQArrRTv4JEoQma9tmSa+S7zWrRkPpgzSdg3v+nLcieUFxGR2lLIlE5KCZlJt25V2pIZ95RLUYbMOFoyx43zZaX30vbuXZvZBMJ7MEPbbOMDb5wT64uISPQUMqWTf//bl5Mnd32tliGzkta45uZ4Q9LSpdF1l8fRkhneyzh+fGXHm8Ff/+pnHIhb9ntstVXX1xctgg8/jL8eIiISHYVM6eSpp3z5Rp7p9NPYkhlnyFyypPrRznG2ZIbBvNLW3LlzYeZMuPnm6OpUSPZgsmLBPYnAKyIi0VDIlE622MKXF1zQ9TWFzM6iHPgTdUumc/DCC/5xtbMBfPpp9fXpTqlB+DvfibceIiISHYVM6SQMZWELW7ZahcxKB/7EvUJRPY8uHzMm87ja+1KrmaKpVNk/p2Kh+Prr46+LiIhEQyFTOqnXgT/1eE/m3LnVh8zevf3XggXR1Cl0wgmZxz+rcvmEJOYczQ7C6hIXEekZtH6GdFKvIbOS1rQXX4T334e334Z11om2TtOmwUMPRXOuPn38ajZx2Wyz6o5PoiUzO2SOGBH/+4mISPzUkimdFAuZtVhqECoPma+84stHHom2PgBvvhndueKej7La1tYkQubbb2ceH3xw19erbY0VEZHkKWRKJ8XW465FS+bSpf6rT5/Kz3HqqdHVJxTl0pp9+tR3yIxzrtHQeef5crvt8r8+c2b8dRARkWgpZEonYUtmvmBSixV/Pv/cl4MHJ/eepYjyvsHp0+Md0FJtS2S+QWBRC+9JLRSIR43yZdS3PYiISHwUMqWTersnc+5cXw4alNx7liLKlszQa6/Fc29mtUtr/uEP0dSjmKef9mWhFuvdd4eRI2GNNeKvi4iIREMhUzpRyCxNHCOgN9sMzjgj2nM++2zlx4arP1W6YlAlit0WEfdsASIiEi2FTOnk9tt9WS8hc84cX/bk7vLf/CbzeOLE6M5breWWS+69wm7wAw4ovM+iRfDcc9FP9yQiIvFQyJROwpBTL6PL67UlM8ru8qOOyjwOr3E9SHK+yvnzYa214HvfK7xP2KV+5pnJ1ElERKqjeTIlr3whs6nJd2cqZEYbwLIH1tRTyBw+3Ldexz0Z+xtvwIwZ/nEp94++9Va89RERkWioJVPyKhQs+vWrTXd5JSHz8MOjrUu2KFsyw6UlAZ54IrrzQvVh+MQTq58CqTuzZ5e3f5KzG4iISOUUMqUsSYfMsCWzknsyi3W9VivKlsy4Q1w1Wlp8CLz77vjeo62ttP3WWsuXSS8GICIilVHIlLIkHTJnzPBd9NmtfaXafntfVnJsd+JcBjJK1U5fNGCAL8N5KuMwdmxp+229tS8VMkVE0kH3ZEpZBgyAefOSe7+pU2HddSsLS/37w957w6xZ0der1Na3WoiylTWOgJ7r8stL2y+8X1Xd5SIi6aCWTCnLqqvChx8m937t7dUtKdmnTzytjuG9ovUoey7JagNn2JIZp1Lvtw1XLlLIFBFJB4VMKctGG8FLLyU3vU17e3XLIvbuDS+/HH19zzor2vNlu+666o7PDmFrr13duZJoydx55/L2V3e5iEg6KGRKWdZYw3cVL1mSzPt1dFQ3hU44YOVvf4umPqGhQ315wQXRnhfgO9+p7vgwZF50UfXLMCbRklnqZ2mrrXy5yirx1UVERKKjkCnLlDItT7jqT1L3JHZ0VNeSGXrjjerPkW333WGTTeBnP4vmfE89Fc15IBMyV1yx+nMl0ZJZ6lKRo0f7crvt4quLiIhERyFTlglblNZfv/A+SS8tWW13eSjq7vJqW1hzldtlXEwYMqOY2L2eWjKbmmDECN2TKSKSFhpdLsuEv+xPPLHwPv37+zKpkPn447DCCtWfJ8qQ6Rx8+mk04TeUO3reucqnHwpDWDUDpkJJtGSGn7t//rP7ffv2hZkz462PiIhEQy2Zskw4CrtYOEmyuzycKunTT6s/V5Qr9Fx/vV9He8qU6M6Z6+GHKz82/DnWe0vmwoVNjB4NEyb4LvA99uj+mKVLM2uYi4hIfVPIlGXCFqViK9CEryU18CcqUbZkPvSQL+OcL3S//So/Nsru8nxr2Efl7rvXWDZH5qRJpR0zYgQsWBDtHw0iIhIPhUxZppSWzPA+xPb2+OsThbDLOcpQUu0qOnG77z5fljqgppjslsyof+bt7eVfyMcf9+Wdd0ZbFxERiZ5CpixTSktmeB9iEi1JUYSasDUvygnZ6z1ktrb6ctttqz/X6qv7kfQAn31W/fmyNTVV3rysuTJFROqfQqYss2CBL8PBPfkk2ZIZRZBN81KEJ51U2XHPPefLaufIDB17rC/nzo3mfKHssD5yZGnHjBjhyyQGJImISHUUMmWZTz7x5UorFd4nyZAZxXuEQW3gwOrPFYqiG7oU119f/jFxLKEZ3pd5zDHR3tua3ZJ5/PGlHXPTTb6M4n5TERGJl0KmLBN2h4ar2eSTtpbMiy7yZZSjpJNsFS13gNWsWb786U+jq0MYMsePz3TFRyG7JbPU6ZYGD/ZlGlumRUQajUKmLFPK1DdpC5nh/aVRtvCFg0+ScMYZ8NZbpe//3//6cpNNoqtD9qTzH3wQ3XmXLs2kzGL3AWdL8+0PIiKNRiFTlgm7gYtNW5O27nIz30oWRzdyEq65Bg45pLR9f/1r2HNP//jzz6OrQ3bI/MY34PTTM/fvVuOGG9ZZ9lghU0Sk51HIlGXKCZlJjC6P6j369El3KCk10J17bubxV74S3fvnLp/5u9/BVVdFd34oPtgsm0KmiEh6KGTKMmHILGUKo7S0ZEK6WzKh9LCdPShn7bWje/98f3RE/fMvZbUfiGdKKhERiYdCpiwTDjIppSXzhRfir0+ULZlpCCXDh+cfdFXr1W3yrdFe7VyhM2Z0fp7bWlpIGDK1tKSISP1TyJRlSukuD8PFeefFX596D5mbbRbt+aZO9aPDjzmm8/bu5qdsb4drr422Lrnnj9q++1Z2XBgy77knmT90RESkcgqZskwpITPJNcvDcHPhhdWdp08feO+96uuT66WXoj8ndG0lXH/94vuPGQPf/W7m+Xe+E2198v3Mq23JfOedyo7LbvG8667q6iAiIvFSyJRlnn/el8VCZpIDLsIQEc6NWKkpU2DcOHjsserrlK3ULt5y5Qa4FVfMPP7tbzMr+oQmTOj8POqWxzhCZvbn6Morqz+HiIjUH4VMWeaWW3xZbOBPkr/Y//hHX778cjTne/31aM7T3Bzv7QLbbNP5+dy5/uutt+DMM2GnnTq/nrsyUNQhM9+tBlGGzP33r/4cIiJSfxQypYtiLXTLLZdcPaKeSL3UuRiLWbLE31YQ5QpCuU47LROwwQ9yGTwY1luvtOOTaMmcORO+9KXyJorPlv0zLffe2yFDup5DRETqj0KmLLPRRr4s1l2+yy6+3Gef+OsTLjVYT2EinLMyzpDZ1JRZc70SUa+tvs46XbddeqlfZjKKAUflhsxBg3zZ1lb9e4uISHwUMmWZ9daDrbfufr/ddoM5c+KvT7g04hZbRHO+KFr4wmBT6uTh1Zg0qfBr2S2duaJuydxtt8KDr4otQVqqcrvewxbpW26BTz6p/v1FRCQeCpmyTHt7/jkRc62+enT3SRaz666+zB45XY3TT6/+HEm0ZIa23bbwayefXPi1OKYc2m+//NsruS/yoosyj1dYAdZdt7zjs1czqrS7XkRE4qeQKcs88gj85z/d73f//TB/PvzrX/HWJwxLxbrvS/HQQ9XXJZRkyOyOme9q3njjztuj7i4H38K9885dt1cSMs8/P/P4nHPKb8n83e8yj8NbKkREpP4oZAqQGdxRyn1u8+f78rXX4qsPZEJmtVMFRXn/aBgyk+guB7jkkuKvH3VU11HzLS3x1OXOO+Ggg/yAn9BVV5U/xVT28peV/AGRPYArisFcIiISD4VMASpbXafaFsbuRBUyo/T5574MB5/E7Qc/KP76nXd23Zbd0helVVeFe+8tf0WiYg48sLo61XrJTRERKUwhU4DKuljjDn9hneopZL76qi+TnMqpHMcfD8svH+97rL125ce+/Xbn5yNGVHaesOs+qkFhIiISPYVMASobLFLthNzdieqeTMiEkmqXlzzrLF8m1ZJZruyu6LgMH9512+9+V9qSo1ttFU0dzj038ziJ71lERMqnkClAeSHzT3/yZVz3/oXCOpUy4r074So5I0dWfo6nn848TrIlM3fZyGKSCFybbtp12+mnw4kn+vd///3Cx2Z3re+554yK65D9h0cp4VZERJKnkClAeSFzxx19Gfcv9/b26LrKwwEiud215QinVIJkWzK33x4OPzz+luNyZF+L0F/+Ag8+CGuu6Wcq6M6mm1Y+2Wp2yMwerS4iIvVDIVOA8kJm1Ms9FvLLX0Y352PUg5SimIS8HHfcAU891f1+a60Vf13Ah8h8c3U+9pgvr7jCh+Lx4/29tTPyNFo2N1fe7Jo9qvzXv674NCIiEqOYxwdLWpQT5sKBJR99FE9dQlF2/cY9Ej4JxW5PuOwyf6/k/vsnU5cBAzLLkGb7/e99+fjjvnzwQX97xQ03dL7dAKBXr8p/wD3h5yki0tOpJVOA8kPm4MHw4Yfx1Qdggw2iO1d2KLnjjurOdfDB1R1fqS237LptjTV8ecgh/ivJeSNLmZC+udkHTMisex9aa60Fkb13vqmcIuecv0H2sMN84m9q8uXhh8PEiRqBJCKSQyFTgPK7pfv2rWy1l3KssAJ8+cvRnCs7fB15ZGZS9Uq89FL19alEvnsyV1nFl7WY5qlfv+73KVavTTapfILNIUM6Pw+DbGyWLIGjj4Y99oB77vEfIOd8effdfvvRR2sUkohIFoVMAcoPmX36xB8ylyyJbtnA3O7Vs88u7/jsSb+Tvh+zmHvvhSuv9INtklbKQKRCk6VXu9Rn7lygsQ6Kcg6OPRYeeMCHytxvqqPDL4N1//1+P7VoiogACpkSCEPmsGGl7Z9ES+bixdF1/+aGzFJGP2f7858zj+slQ6y8su8u//73a/P+pfxsLrgg//ZqQ+Fyy8Hmm2eeRzHNVUETJ8LYsd03f7e1+f0mTYqxMiIi6aGQKUBmdZ0rryxt/yRC5pIl0YXM3POUGxTnVD7bTiw++gjefLO2daimlbmUrvbuvPxy5nGsLZmXXeYDZJaPWYk28nwTbW1+fxERUcgUr9x1wpMKmVF1l+dOnl7OmtevvNL9GuJJW2UVP/iqlvbZp/Jjq5kUP2SWmbN1/vzqz1fQQw91+sA4YBgfsxtPdt23o6P6ewFERHoIhUwB4JNPfFnqSjZ9+sQ7T+bcuTBlSnTh4bDDOj+fPj3/3I35XHxx5vE++/hxH7Xy3//6r3rQt68fBPXVr/rn665b2h8eq68eXctjOCn80qV+qqTp06M5bydBK+aPuJiv81dW5wMAJrE9t3EkW/JS3v1FRBqdZpsTAN5915frrVfa/nG3ZF5+uS/vvz+a8/Xq5bvMw8G/ixf71sBSus2zB0Xdc09tB/5svHHt3jufLbf097eecgqcemqyUygB/N//+T8Chg/3y1rusAP8618Rv0n//rBgAb/hR11eOprbAFhIX/qxKLO/iIiktyXTzKaZmSvwlbc9w8x2NLOHzWyWmbWZ2atmdoaZFewkNrP9zKzVzOaY2Twzm2Bmx8X3ndVGGL5KDVBxh8yZM6M/Z75AObeEWXSyQ2ZU3fc9SVMT/OEPsMkmpbVQRj3J/nLLweef++cffxzduZfZd99uRxbNIbh3oanJ7y8iIqlvyZwD/DbP9nm5G8zsQOBuYCFwBzAL2B+4AtgJOCzPMacCVwGfAjcDi4FRwBgz28w5V2d36lUu7PouNUT17ZvpYo+zPqVM+F2N/fbrfrnGcFAU1Nf64Wly2WV+NPj99/sWx6g98IAvY/n5jB4NDz8MRW7deImt2JvH/Iim0aNjqISISPqkPWTOds5d0N1OZjYIuB5oB0Y6554Ptv8UGAeMMrMjnXO3Zx0zHLgUH0a3dc5NC7b/HJgEjDazu51zUXfO1UQlITPOlsywtSvK7td8LWi5Sx3mE9X66Y3iL3+Bt9+GUaNgs838tpNP9n8wRDW5frawFRNiCpnbb+/X67y98C6fsoLvJj/gANhuuxgqISKSPqntLi/TKGAl4PYwYAI45xYCPwmefjfnmBOAvsDVYcAMjvkM+GXw9OS4Kpy0ckNmEpOxh+8TlULdtGut5cPmzJl+SsRQWxvsvvtIDRYu07HH+vkxN900sy2p2xRjCZlmLLr+pqK7LOwzCA48EG66Sc3dIiKBtLdk9jWzrwNr4TuzXgWecs7ltj3tEZSP5jnHU8ACYEcz6+ucW1TCMY/k7JN6lbRkxjm6PAyESdwD+d57vofzk09g6lT/3m1tflnLbCutFH9depqOjmQz1/z5cMUVcPrp0U7QvmBJ8Sb1haefA79ZO7o3FBHpAdLekrkK8FfgF/h7M8cBk81st5z9NgzKLtNXO+eWAlPxgXudEo/5CB9q1zCzmO8aTMbixT4M1Ms8meFApCjX5D7kEF9ecknX1yZN8gEz9N57XWeiqeXURWmVdKPehx/CWWf5WyiTtHCYAqaISK40t2TeCDwN/Af4HB8QTwVOAh4xsy85514J9g2nrS60bku4fUjWtlKOaQn267LenJmdFNSFYcOG0draWvy7qdK8efOqeo8pU9ahuXkNnnyym1EwgZkz12P+/FVobX2m4vcs5v33vwCszKJFC2ltHR/JOU880Tj88N6stNJiYGTB/VpbW5k+vR+wQ6ftS5e2EvOPsaFU+5nNGNlly4QJ/2XgwGiGmn/+eTNz5/YGvthp+5lnvsEVV/i/RV9//W1aW9+N5P2qFd11lVy6tvHQdY1HXVxX51yP+sIP1nHAvVnb3gy2rVfgmGeD17+UtW1xsK25wDEfBK+v2l2dttlmGxe3J554oqrjzzzTueWWK33/s892rl+/qt6yqBNPdA6c+81v4jm/7xTP//XlLzv35ptdt0u0qv3Mhl57revP6txznZs8OZLTF/yczJ6defzTn0bzXlGI6rpKV7q28dB1jUc11xV43kWQydLeXZ7PtUG5a9a2sDWy0EJ84fbZFRxTZ6taV2bx4vLufwwH/kQ552G23r1h+eXh7LPjOX8x//gHTJuW/PtKZTbZpOu2X/0K1l8/3vcdODDzeOFCP6JeREQyemLIDKfxbsna9kZQbpC7s5k1AyOApcDbJR6zanD+951zXbrK06jckNm3rw+Y2XNIRumTT7oOvElSuFSiSCG9esFdd/nHt9zil9V85JHix4iINJKeGDLDG+myA+O4oNw7z/67AgOA51xmZHl3x3wtZ5/UKzdkNgd38776ajz1mTHDL/soUs8OPdSXH37oy1deKbyviEijSWXINLONzawlz/bhwNXB05uzXroL+AQ40sy2zdq/H3BR8PQPOae7EVgEnBqcNzxmKHBe8PRaeojFi8ub+Px///PlYV3WSYrGjBkwbFg8565Ec5qHyDUwM9hnn+jO9+c/w0svwb33Fn4/ERHx0vqr8wj8ijtPAe/gR5evC+wL9AMexg8AAsA5N9fMTsSHzVYzux2/ks8B+KmK7sIvNUnWMVPN7Gzgd8DzZnYHmWUl1wAucz1ktR8ovyVz0CBfLojpZoGPP4aVV47n3JV4s8tEVpIWUXZhv/UWHH88bLll/tejnJtTRCTt0vpf4hPAg/hgeTRwFrAb8AxwHLCfc67TVOHOufuCfZ4CDgVOA5YExx4ZjKYi55ir8EH0P8Cx+CmJpgPfdD1o3XLwLZPl/II8L2jLPf746OviHMydmwmycfjkEzjjDL9Ay3PPFd/3j3+EESPiq4ukR74W7SOOyDxWS6aISEYqWzKdc08CT1Zw3LNAWZ1nzrmxwNhy3ytN2trgP/8p75iwlTGO5QIXLfIDipZbLvpzh1ZYwa8ME5o8ufBo5JNOiq8ekoz99oPDD/dLXpZiyhT/b+LjnKk2w0UCst1+O9wR9IP07VtdPUVEepK0tmRKhD7/3JfZa013p1cv32qT75duVPWJM2Tmyr0f9Sc/gb/+FY4/fmr+A6SujBsHaxdZdOehh+C440o/3/rrw0EHdf0Do9DnPVxN6vvf91NviYiIQqbg5/gD331cjj594lm/vB5C5nrrwde/Dsce+05ylZCK7b47nHNO9/uVMq/r/PmFXysUMv/2t8zjzz7r/j1ERBqBQqYsW6O73K7v3r17TkvmaqvB9dfDbsGq9+WMtJf6EH5eVlqp8D5NTTB7dvHz7Lln4deGDy98XhER6Uz/NcqylsxyQ+a8eXBtDJM41SJkAnz72/4+zdVX12TsaXTkkXD55Z1bFfOZNKn46xMm5N/+2GNw2mmV1U1EpBEpZMqylsx+/So7tr092vrMnevLOEeXF7LVVvD++7Diism/t1SnVy8488zuB9989avdB81Cx6nFUkSkdPovUypuyVxjDV/+7nfR1mfWLF9qAIVUopQ/erbfPv/2QvNffv/7FVdHRKRhKWRKxS2ZYRh8991o6xOet5Zrl0t6ldOy7hy8844fxHb99YWXhbzyymjqJiLSSBQypeKWzDCULl0abX0+/dRPjzRkSLTnlcZQ6lyVM2f67u/hw/0AtmrnQz355OqOFxHpaRQypeLR5eEv8zCkRmXWLB8we/WK9rzSGLbfHq65Bi65BC67rPB+pS5bWupMA1dfXdp+IiKNQiFTKu4u32EHXw4YEG197rzTj/AWqYQZfPe78IMfwFlnVXeu3/8e/vvf0vbt1Qt++EP/+KOPqntfEZGeQCFTKp4y6LrrfBkOAIrKxx/rfkyJzne/W/mx3/qWn5i/VOFgtdVWg6FD45lHVkQkLRQypeIpg8J7JsOW0CiE93cWmxBbpBzXXANTK1wdtE+f8vbPngh+9uzMIDYRkUakkCnMmQMtLeXfA9nc7L+ivCdz0SJfVjJnp0gha68Np5zSdXu4wlMhZuW9T+5gtagHxYmIpIlCpjB/PgwcWNmx/ftH25IZBlaFTImSGZx7btftt9/edduvflX5++S2fEY9KE5EJE0UMoUFCyofvBN1yKxm9SGRYlZfHZ59tvO2VVbxXdqbbZbZFrbob755+e+xzz6+1TQU5b8NEZG0UcgU5s/33eWV6N/fh9SoVDpnp0gpdtwRXnut87ahQ+Gpp2DXXf3zXr3gs89g/Pjyz9/U1Hmg0aOPVl5XEZG0U8iUqloyBw2CefOiq4u6yyVuG27YdduQIbDFFv5xr17+eaV/6Ky2Wubxxx/DtGmVnUdEJO0UMqXqkBmOTo+CQqbErVcv2HlnuOuuztvD5SirXQTgyCMzjy+5BEaMqO58IiJp1VzrCkjtzZ8Pw4ZVduyQIfDQQ/7esyi6uBUyJW5m8PTTXbc3BX9yl7rCTyHVHi8i0lMoZEpVLZnh5NO77gqTJlVflxtv9GWzPpmSsAsvBOfguOOiP7dz5U+HJCKSduouF15/vfJRsLvs4svnn4+mLn/+sy/D+TJFkrL88n798Sha0Z94ovNKWK+/Xv05RUTSRiGzwU2f7ss5cyo7Pup1y/fd15df+Uq05xVJ0siRsPHGmeebbBLtADkRkTRQyGxw//ufLytd33nUKF/ut1/1dfnkE39/J2TujxNJq8cf7/z8lls6P7/mmsy/PxGRnki/yhvcddf5crnlKju+b1/YaKNoBv1ccUX15xCpF3vs0fn5okVw330werS/R/OUU2DbbWtSNRGRRChkNrittvJlOBF1JVpa/Ah1Ecl4+OHOz8ePh4MPhssvhyVL/Db9uxGRnkwhs8F1dPiymmlXogqZues+i6RZ376+xTJ0222Zx4sXx/jGzsGECXDYYf4fZ1OTLw8/HCZO7FwpEZEYKWQ2uCgmoI4qZGp+QemJDjqo67bYBgEtWcLGF13k++rvucfPT+acL+++228/+uhMU6qISIwUMhtcGDKrmZdy8GC/1nO1wqB72GHVn0ukXtx7L5x6audtq66aeXzxxRG9kXNw7LGs+OyzPlSG3RShjg7/1+D998Oxx6pFU0Rip5DZ4JYu9WU1o7lHjIB33skE1kqdc44v//KX6s4jUm+KTcl17rnw7rsRvMnEiTB2LL26m2S2rQ3Gjo1m9QQRkSIUMhtce7tvQaxmNZJhw3xYffbZaOqkJSWlp+mup+C99yJ4k8su67SqwmJ6M4uhLKRv133b2vz+IiIx0uJ9DW7p0uruxwQYOtSXu+1WeQ9cds+elt+TnmbrrYu/vvPOEfReP/QQV3acSgdNHMMtDOPjZS85cv5RdXRkJqUVEYmJQmaDa2+vfp3wlpbq6/Gf/1R/DpF6tcoqPtdFtcjAL3/p/7jrtIhCWxtncCUAN/CtTvv/m03ZjNc6n6TStWRFREqk7vIGF0VLZvYvzkonVP/gg+rqIFLvzGCddQq/fv31cOutpbVo/vjH8L3v+cfz5/tzm8t0B/yHTTvtfxejup4kihUURESKUMhscFG0ZGaHzLPOquwc++9fXR1E0uCtt2CXXfK/dtJJcMwx/t/TnDm+oTFf6My+f/Oss2DgwO7ft0t3eVMT7LtveZUXESmTQmaDi6IlM4p7KMNR7iI9XSn/XoYMgSOO8KHzlFM6v/Z//5d5XGrPwcPswyPszUes4jf06+fXtxQRiZFCZoNra6s+ZA4e3Pn5L35R/jk23ri6OoikxVprlbbf2LG+/MMfOm9/5pny3/MFtmUfHmE3nvTd5AccANttV/6JRETKoJDZ4F57DTbcsLpz5K57/pOflH+OKCZzF0mDa66Bm26Cp58u/9hJk+D11yt/78lsAAce6CugaRxEJGYKmQ1u5kxYe+3qzlHp76rXX89M1ffZZ/CNb8DUqdXVRaTeLbec/6yH0xbdemv3x7z8Mrzyil8ZslR/3OaPebcv/ettWsNVRBKhkNngPv0UVlih+vPkdul95zu+fOEF2GsvWLy46zE77AA/+IGfvmjRIvjCF2D48OrrIpImRx0F110HN99ceJ+ttoItt4Rf/7rwPuefn3l8yilw0vPfYfbsrvvNnVtpTUVEyqOQ2cAWL4Z586IJmSef3Pn5ddfBG2/4wQt//3smdIJvtXQu88tu02C2lY8/RqQhnXiiH+TjXPcTt+ez5ZZw4YWZoDpsmC9z75cOX/v884qrKiJSMoXMBjZzpi+jCJkA553X+fmJJ/opWwDGjPHLTu61Fyy/PJx+etfjjzkmmnqIpNn48eXt/9xz8NJL/vFRR8HZZ/+Pc84pvP/SpTBoEDz8cOV1FBEphUJmAwtbH6Oak3mzzTo/zx3YsN9+vlUT4Kqruh6/zTbR1EMkzXr3hokT4Utf6n7f8eM779fUBPvsM73TLZerrpr/2CeeqK6eIiLdUchsYA8+6Mt890tW4ogjir+e7/4wEelqu+3gySe7bj/pJN+lfvrpsNpq8MUvdn+uqVNhwYKu2zW4XETippDZwHbbzZdf/nI05zPzS9yJSPV694bDDvOPzz3XD9K79lr//Le/LX0p1r59fW/F5Zd33n7JJb4888zMPZwiIlFSyGxgm2/uBwYUW0+5XAMGlN9iucMOmidTJJ++fX35hS/4e5mraX3Mt3Srcz6watCdiMShylWrJc0WLozufsxsgwb5X4a5ay4Xcu21fhk9Eenswgv9WuX5AmK51lyz67Z8A/BERKKilswG1tYWT8g0g46O4vvssIMvv/Md2GKL6Osg0hOssw60tuafiqhcffvCb37TeVv2ALw33/TlF78IBx1U/fuJiChkNrC4QmYo/KUFcOednV8LR5IfdVR87y8inZ19duHXNtzQ/58wcSLcf39MFXAOJkzwN5u2tPjh8C0tcPjh/o1L7f4QkVRQd3kDiztkrr++v9fyk0/8Sj5nnOFHzI4ZA+uuCzvu2HXdcxGpnRNPjPHkS5bAscfCAw/4e3XC7o4FC/x6mQ8/7O8LuOkmLXsp0kOoJbOBffgh9OsX73sMGQLrrQfNzXDFFfDii37AUUsLHH20plERSdp++/mvfG65JfP4lVcifFPnMgFzwYKu99N0dPipKe6/3++nFk2RHkEhs0FNmQIvv+xX4RGRxjF2rP/qzpZbwk47+Xk2qzZxon/TfBN2Zmtr8/tNmhTBm4pIrSlkNqj33qt1DUSklv785+73ee45P/jotdd87nvyyQrnwr3sMh8gszhgFkO77tvW5vcXkdRTyGxQ8+b58vrra1sPEamN44+Hp57yk7Hvu2/xfTfbDLbfHkaOhIEDob299PeZMgUefKAdOjq4ku+zCh8xkxX5I99hBWbxJut3PqCjAx56qOzvR0Tqjwb+NKhZs3y5xx61rYeI1M4uu/gv8L3UBxxQ2nHNzf5+7oUL4Xvfgx/+0C9Pu36QF19+2Y/lGTUKDj4Ypi66m3dYizO4EoBTuZo5+HmZprAeGzC58xvktHqKSDopZDaoTz/15fLL17YeIlIfsid833//7u/bXLjQl9dc478AjjgChg7NLH950UWZgeJr8+6yY+/kiGWP/8Jx7MLTLMe8zMnjnPZCRBKjkNmgZs2CXr2imeRZRHqGmTP9PZdDhsAvfwn77OMnZy81891xR9dtS5YUP+ZOjqCdXgxiLrNYnvuaDu2+/15EUkEhs0H98Y/+vipNISQioRVX9F8Av/51Zvull8IPfhDf+97NqMyTfv1g9Oj43kxEEqOBPw1owQI/QbqISClGj/ZTGT3+ePzv9adNrsBtu138byQisVPIbEBz59a6BiKSNsOHw5e/nP+155+P7n1OnHQSw0eoi0WkJ1DIbEAKmSJSqYkTfbnXXn6985/8BLbaqut+b72VeTxjBmywAfxgtOOte17h62s9xZP99lr2+saD3u907Lvvwve/7yeEf/rpGL4J0DrqIgnQPZkN6PXXffmFL9S2HiKSPttt55ec3HRTn8tCf/sbvPQS7LqrH3m+zjqZ11ZeGd54A8CALfjrweDcrnz9WJ/xBg9eg5EjO7/PVVf5ctdd/RK4DzwAa6wR0ZggraMukgiFzAY0ZYovn3uutvUQkXTafPOu20aN8l+5jjkm/znM4K9/zTzfYQcYPz7/vqutlnn8wAN+ft+WltLr20nuOuq5ctdRv/VWjZAUqZC6yxvQP/7hy0GDalsPEenZFi3yjYGlGDOmtP0OOMCvOvSPf/iJ5M1g8uTuj1tG66iLJEYhswE9+qgv9ce5iMSpT5/OXerFbLhh97kv21e+As884x9vsIH//8wsM/vRQw/BDTf4hsu33oJ77vGvH3IoXLbgZJ5hJwA+YhX+y8b536QW66jrXlHpQdRd3mC6mxhZRKRW+vWr/hyXX+6/QqeemlmdCODeD77IvXwRAIexGh8B8D1+z2/4Ia+yOavxIW+wIXt2/JNeSa6jrntFpYdRS2aDueceX4brFYuI1Aszf8/47rv7hruBA+E3v4Ejj/Svn3RS+efMDphd3o9Mq+A1nMJA5rMj/2I477AXf+dCfsbvF3wzkflBu9wrGgbMUO69omrRlBRQyGww4X/Wv/1tTashIpLXuuvCuHF+icrPP4ezz4Yf/tC/ds45/rVTT81/7L/+FW1d/o/zOdVdzVe/6gNweztMn555/aWX/BK9kajne0XVhS8VUshsIH/4Q+bxZpvVrh4iIuXYaiufY0aM8K2cV13lJ4B/6aXMPptt5keox/kH9Je/PJJVV/WBc9VVYeutYYUV4LXXfIvreed13n/BApg9u8STX3aZD5ClSPJe0SVL4Oij/ZD+e+7x35RzmS78Pfbwr9fiXqw0ht801rkKCpndMLM1zOzPZvahmS0ys2lm9lszG1rrupXr/PMzj3U7j4ik2Tbb+Mna//pXuPBCeOIJv/300/3v6Vdf9bMPrbEGnHkmnHBCZqAQwDh2L3juPfhnt++f3aK52Wbwox/Br36VGYA0YoTPDkOH+gB6+eVw1lmFGyo/evAFOjo6B4xPWIE/8a2uO3d0+JFNcavnLvx6Dr+FpLHO1XLO6avAF7AuMANwwH3AxcC44Pn/gBVKOc8222zj4vbEE090u4//NDt3/PGxV6fHKOW6SmV0beOh61rcP//R4T488GTn+vd3t3Kk+ztfdv9lI9eOuQs4393DQW4pTcv+v3zh+Q533HGZ/z+j+BoyxLk+fTLP//hHX36f37oOcO2Yu4pTOh2zKa92PklTk5s717kTTnBu1qzM9/ePfzi3447OLV4cwcUaP965lpbSvqmWFucmTKjobcr+zHZ0OHfkkc4NGFC8Tv37+/06OiqqV6RqUOdq/i8AnndR5KgoTtJTv4DHgkB5Ws72y4Pt15ZynnoImZMnZz7Dd90Ve3V6DP3Cjo+ubTx0XUuweLH/Rd7S4lxTU+df8k1Nzg0Y4J7c4wL3/tRMUvvsM+e+/e233DPPODdypHP77dd97orja2P+407iWndwr/s7bV999a77/vnPzk2a5Nzyy/u8+Oqrzn34Yf7LMX16nut02GFdr0+hr6Ym5w4/vKIfR9mf2YTCb6RqUOd6CJnqLi/AzNYFvgpMA36f8/LPgPnAN8ys0nUnEnPDDbD++v7xPvvAoYfWtj4iIjXVu7fvSx83zv+HmH1v3KhR0NrKrv/8GasPz9xXNGQIHHPMu+y0k++aHzvWJ4L2dnjkEbjvPj/D0OTJvgd5wgRYb73oq/46X+A6vsO97Qd02v7BB133PeEEvwzorFn+ftXNN/erJzU1wcUXw1pr+a79Pn1glVX8gKobbvDbfvc7+M/Yt/lVxw+ZzwBO4AaO5hbO5ZcYjpv4Bh+zEtMZRhv9eLJjZzoefJj2dn9Nwjp1dPhlSB9+OMKLUK/3rxaTxjpHQPNkFhbesPN351ynG1Gcc5+b2bP4ELoDlHADT4I6Ovx/Eq+84m+Yz3bffTWpkohIfTGD7beHO++s6jRNTbD33l23b799ZiWiGTPg4499MAVYvNiH1jPPhC228P9XZ/sF5/FjfgnAYGYzhyGcw6+4mHOrqmvIOTg3z6l+n9WccvrpAM8DcB6/6rLvceRZymkBFaaKkYAPuhtuCE8+Cb16wc9/7m/3HDoU3ngD+vaFTTeFvcauRHvHibzHmjzNLoziLhbThwfZjz0Yx1q8yzBm8B5rslfHY3zwwCe8eDUMGOC/99ZW2HVXf1/vuHGw/PI+iB9xhA/Is2bB44/7n09zM7z7rv95PfmkP65fP+jfP1P7JUv8LbI77QQrrug/Ws75Fa/mzfPbeOghlnQ08RGrsybv8SS78Qw7813+wAr4KQocYFDwnlvn0reIikJmYRsG5ZsFXp+MD5kbUOOQ+b3vwaOPbkufPv75G2/k3++99zTgR0QkacOG+a/cWT1OOcWHmBdfhJYBjgE/Oo0hj/+NQQs/5nhupDdLWJFPl+3/K87jb32OZvyIo/n+o/swf4HRvz8MHuzH5my6Kdx2m1/VbautYOZMv3zwI4/41lWAvfaCqVPhzUK/2Wpo+vTMgKr2dvjxj7vu88YbcHdO5+LT7Lrs8ZNBYO1kIXBa500339x1t7PP7vz8nHOK13ellfw1zrXaavDZZ7kNl/PznuOnXMTOPM0z+Mmrm1nCUnozbP50ZgSBsl8/P5vBhx/CJpvACy8Ur1c9UcgsbHBQzinwerh9SL4Xzewk4CSAYcOG0draGmXdOpkyZSM++mhFoJ3lllsCdF424/DD3+Nb35rKlCkdTJkSWzV6pHnz5sX6s2tkurbx0HWNT5zX9vPPwb5/EBvNf52W555j2KIZNPkxAAB0mOH69GHXnaaz8rl9mTrtyU7HDx/uW832399/ZTvlFF8uXmz07u2KtoZNnjyQ9defx4IFvdjq1z+j46n/cQ+Hchh/Ywiz6cciOjAeYy9u5Wh25Sk25TU241U+3XUX/nfeeZjB448PY4UVFtHUBB991I8pUwby+efNvPjiUA4++APWX38eP/3ppgCsvfZcZs4cwIIFzQwatIQBA5ay3HJLGTx4CSNGzOeNN5ZjwIB2XnxxCGuvvYDJk5djbx7hUb4GwIrM5FDu5u98lamsA8BavMO7rM2avMt7rFXdDyfLWmvNZ/bsPsyd25sBA+YDne+YW3nlhfTq1cHQoYZzfWhvN5ZffjHzZixlPgP998s03mE4AFvyEovoy8rM4GOG0YfFDGABa9gHzHCr0K9fOwsX9mLaNIdzxsKFc2htfamkutbF/wVR3NjZE7+A6/Ct198u8PovgtfP7e5c9TDwRyqj6xofXdt46LrGJ5Fr29HhB30cdlhmYFJLix9UM3Fi/O+frZzBKgMGJDe6PKEBSZGqQZ018Ke+hS2Vgwu8Hm6fHX9VRESkIWTfKzpvnu83njfPL4G03XbJ1mX77X2zaPYNiPn07w8HHJBc/UaP7r5OoX79/P61lsY6R0Ahs7DwzsYNCrwejNcueM+miIhIepnBTTfBgQdmRuBna2ryo2kOPNDvl9SolHoNv8Wksc4RUMgs7Img/KqZdbpOZrYcsBN+LN34pCsmIiKSiBKme+K225IdVVqv4beYNNY5AgqZBTjn3gL+DgwHTsl5+UL83b5/dc7lHzImIiLSE9RTF36oHsNvd9JY5yppdHlx3wOeA35nZnsCrwNfxM+h+SaQZ4IFERERiV1Ec50mKo11roJaMosIWjO3Bcbgw+Vo/HrmVwI7OOc+LXy0iIiISONSS2Y3nHPvAcfXuh4iIiIiaaKWTBERERGJnEKmiIiIiEROIVNEREREIqeQKSIiIiKRU8gUERERkcgpZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiEROIVNEREREImfOuVrXocczs5nAOzG/zYrAJzG/RyPSdY2Prm08dF3jo2sbD13XeFRzXdd2zq1UbQUUMnsIM3veObdtrevR0+i6xkfXNh66rvHRtY2Hrms86uG6qrtcRERERCKnkCkiIiIikVPI7Dmuq3UFeihd1/jo2sZD1zU+urbx0HWNR82vq+7JFBEREZHIqSVTRERERCKnkCkiIiIikVPITDEzW8PM/mxmH5rZIjObZma/NbOhta5bvQiuiSvwNb3AMTua2cNmNsvM2szsVTM7w8x6FXmf/cys1czmmNk8M5tgZsfF953Fz8xGmdlVZva0mc0NrtnN3RyTyLUzs+PMbGKw/5zg+P0q/V6TVs61NbPhRT7DzsxuL/I+ZV0nM+tlZmcGP7e24Of4sJntGMX3HSczW8HMvm1m95rZlKD+c8zsGTP7lpnl/X2nz2z3yr22+syWzsx+bWb/NLP3sur/kpn9zMxWKHBMej6zzjl9pfALWBeYATjgPuBiYFzw/H/ACrWuYz18AdOA2cAFeb5+kGf/A4GlwDzgBuCS4Ho64G8F3uPU4PVPgN8DVwDvBdsurfU1qOLavRx8D58DrwePby6yfyLXDrg0eP29YP/fA58G206t9XWL+toCw4PXXy7wOR4VxXUCDPhb1v8hlwQ/x3nBz/XAWl+3bq7pyUHdPwRuAX4F/Dn49++AuwjGIegzG++11We2rGu7GBgfXM+LgauAScH39AGwZpo/szW/wPqq+IP5WPDDPi1n++XB9mtrXcd6+MKHzGkl7jsI+BhYBGybtb0f8FxwXY/MOWY4sDD4xzc8a/tQYEpwzJdqfR0qvHa7A+sH/5GPpHgQSuTaATsG26cAQ3PO9WlwvuHVfN91eG2HB6+PKeP8ZV8n4KjgmGeBflnbtwt+rh8Dy9X62hX5nvcA9geacravArwbfG+H6jObyLXVZ7b077tfge2/CL63a9L8mVV3eQqZ2brAV/EB6vc5L/8MmA98w8xaEq5a2o0CVgJud849H250zi0EfhI8/W7OMScAfYGrnXPTso75DPhl8PTkuCocJ+fcE865yS7436UbSV278Pkvgv3CY6bh/y30BY4vob41Vea1rUQl1yn8+fwk+LmFx0wC7sD/fEfFVN+qOefGOefGOuc6crZPB64Nno7Mekmf2RJVcG0r0XCfWVj2ecvnzqBcP2tb6j6zCpnptHtQ/j3PP/rP8X/VDQB2SLpidaqvmX3dzM4zs9PNbPcC967sEZSP5nntKWABsKOZ9S3xmEdy9unJkrp2jXy9VzOz7wSf4++Y2eZF9i3rOplZP3zrxQLg6VKOSZklQbk0a5s+s9HId21D+sxWbv+gfDVrW+o+s82l7CR1Z8OgfLPA65PxLZ0bAP9MpEb1bRXgrznbpprZ8c65J7O2FbyuzrmlZjYV2ARYB38fXXfHfGRm84E1zGyAc25BNd9EnYv92gUt86sD85xzH+Wpw+Sg3KCK76OefSX4WsbMWoHjnHPvZm2r5DqtC/QC3nbO5QsLqb22ZtYMHBs8zf6lqc9slYpc25A+syUysx8AA4HBwLbAzviAeXHWbqn7zKolM50GB+WcAq+H24fEX5W6dyOwJz5otgCbAX/E31vyiJltkbVvJde11GMGF3i9p0ji2jXq534B8H/ANvj7qIYCuwFP4Lso/5lza0ycP4shBV6vZxcDmwIPO+cey9quz2z1Cl1bfWbL9wP87W5n4APmo8BXnXMzs/ZJ3WdWIVN6NOfchcH9RDOccwucc685507GD5Dqjx/pKFK3nHMfO+fOd8696JybHXw9he+tmACsB3y7trWsT2b2fWA0fvTtN2pcnR6l2LXVZ7Z8zrlVnHOGbxA5BN8a+ZKZbV3bmlVHITOdumsdC7fPjr8qqRXerL5r1rZKrmupxxT6q7CnSOLa6XOfJegi/FPwNKnP8ewCr9cdMzsVuBL4L7C7c25Wzi76zFaohGublz6z3QsaRO7FB/IVgJuyXk7dZ1YhM53eCMpC90SEo9EK3bMpEHZBZHfZFLyuwb1HI/A3t79d4jGrBud/v4ffjwkJXDvn3Hz8vHEDg9dzNeLnvsvnuMLr9BbQDqwT/LxKOaZumdkZ+PkGX8OHoHwLL+gzW4ESr20x+syWwDn3Dj7Eb2JmKwabU/eZVchMpyeC8qt5VlpYDtgJf0/M+KQrliLhyPvsf4zjgnLvPPvvih+x/5xzblGJx3wtZ5+eLKlrp+vdWb7PMZR5nYIpUJ7D/5x2KeWYemVmP8JPHv0yPgR9XGBXfWbLVMa1LUaf2dKtFpTtQZm+z2zuxJn6SscXmoy9lGu0MdCSZ/tw/Ag5B5yXtX0Q/q/scia6HUEPnYw95/scSfeTscd+7egBE1tXcG23JmcS7GD7nsH364Adq71OlDax9aBaX69uruVPg+/heWD5bvbVZza+a6vPbGnXdANgcJ7tTWQmY382zZ9ZCw6UlAkmZH8OWBm4Hz9dwRfxc2i+if8H/Gntalh7ZnYB/sb0p4B38Ev4rQvsi/9H+TBwsHNucdYxB+GXSFsI3A7MAg7ATwNxF3C4y/lHY2anAb/D/+O7A79M2ChgDeAy59wP4voe4xRci4OCp6sAe+FbH8I56T7J/t6SunZmdhlwFvB+cN4+wBH4+5dOc85dXdU3noByrm0w5cv6+H/v7wevb05mnrqfOucuyvMeZV0nMzP8BNCj8IM5xgb7HoH/93Koc+7+yr/reAXrMI/Bt/pcRf77oKc558ZkHXMQ+sx2q9xrq89saYJbD34FPANMxX+mhuFH4q8DTAf2dM79N+uYg0jTZ7bWSV5flX8Ba+Kn6Pko+NC8A/yWrL88Gvkr+Id6G/4/n9n4SYNnAo/j53azAsfthA+gnwFtwL+BM4FeRd5rf+BJfJCdj1979rhaX4Mqr98F+L9mC31Nq9W1A74Z7Dc/OO5JYL9aX7M4ri3wLeBB/Apf8/CtGO/if1nsEuV1ws+dfGbwc2sLfo4Pk9PqVI9fJVxTB7TqMxv/tdVntuTruilwNf72g0/w91POCb7/CyjQYpymz6xaMkVEREQkchr4IyIiIiKRU8gUERERkcgpZIqIiIhI5BQyRURERCRyCpkiIiIiEjmFTBERERGJnEKmiIiIiEROIVNEJIXMbIyZOTMbXuu6iIjko5ApIg3HzDYys6vM7DUzm2Nmi83sQzN7yMy+ZWZ9a13HOJnZmmbWHoTUX3azr8v5WmRmM83sRTP7k5l9zcx6dXOOFjM7w8zGmdnHwfWebWYTzewXZrZOtN+hiNQDrfgjIg3FzM4Hfob/I/tfwPP4pe+GASPxawa/4JzbtlZ1LIWZjQGOA0Y456aVeeyFwPn45QBnAGs655YW2Df8JXFhUPYChgCb4Je364O/hsc4597Mc/wO+LWPV8evg/xP4EOgBdgK2BHoAHZwzr1YzvchIvWtudYVEBFJipmdhw9L7wGHOecm5NlnP2B00nVLStDqeAIwF7gZ+B5wAHBPseOccxfkOdcw4CrgMOAfZratc+7jrNc3Ah4DBgLnAJflhlkzGwH8GhhU+XclIvVI3eUi0hCCexcvAJYA++QLmADOuQeBvXOO/aaZ3W1mb5tZm5nNNbNnzezrBd6rNehabjaz88xsctDN/J6Z/drM+uQ55iAzu9nM3jSz+cHXC2b2fTOL8v/qrwFrAHcAfwi2nVjJiZxzM4AjgVZgTeC8nF2uwofHXzvnfp2vtdQ5N9U5dzi+VVlEehCFTBFpFMcDvYG7nXOvFdvRObcoZ9MfgLWBp4DfArcHz/9qZv9X5FS3AqcBTwfnaAN+CPwxz74XA1sDE/Dh7CZ8C+CVwF+K1bdMJwXlmOA6vAB81czWruRkzrkO4KLg6VFmZrCshfLLwELgNyWcJ/eai0jKqbtcRBrFzkH5zwqO3dQ591b2hqA18hHgHDO71jn3QZ7j1gU2cc7NCo75MfAKcKyZneucm56177553qMJuDHY/+pCra+lMrPVgX2AN51zzwWbxwDbAN8GflrhqZ8BlgIrA8OBqWSu9wvOudkVnldEUkwtmSLSKFYNyvfLPTA3/AXbFgO/x/+xvmeBQ38UBszgmPnALfj/ezsNLCrwHh34lkyAvcqtdx4n4AfujMnadiuwGDihu1HihQStkJ8GT1cKyoqvt4j0DAqZIiLdMLO1zOz3ZvY/M1sQTucD3B3ssnqBQ5/Ps+29oBya8x4rmNnFZvaqmc3Leo8XunmPUr+HJuBb+JHcN4XbgxA8FlgN2LeatwhPWcU5RKQHUXe5iDSKj4CNKTOsBXM4TsSHwqeBvwNzgHZ81/BxQN55NQt0E4eDX5a1GprZEGASMCJ4r5uAWcG+Q4DTC71HGfbC30f6WJ6u/THAofj7NR8o98Rm1g9YPng6Myg/CsqqwrGIpJdCpog0imeAPfBd2zeUcdxZwArA8c65MdkvmNlR+JBZrW/jA+aFuVMFmdmX8CGzWuGAn72y5r7MtbeZremce6/A64XsjP99MiNrzs5ngnJbMxvsnJtT5jlFJOXUXS4ijeJG/PRFh5rZF4rtmLPiz3pBeXeeXXeLqG6xvoeZrQLsh58b84YCX8/iW1dPKPPcTcCPg6e3htudc1OBfwD9gLNLOE+PXmVJpBEpZIpIQwha2C7Ar1DzkJnlXdHHzPbGjxoPTQvKkTn77YVvgYxCoffYCjg3gvOfgG9pvMU59+18X8A38fdTfqvUeTnNbGX8dE4jgXeB3CUqv48Ptuea2Wgz69J7FtzvejvwpQq/NxGpU+ouF5GG4Zz7ZRB0fgZMMrPn6Lys5K7A+nQesHMNfo7Nv5nZXfglETfFT9h+J3BEBFW7Cd/a91sz2x2YHNRjP/xKPBW/RzBvZRiG/1RoP+fcFDN7Eh8YvwY8lHOeC4KHTWSWldwZH9on4peV/CTnnK8HYfwu4FLgdDPLXlZyC/zSlA6/6o+I9CAKmSLSUJxzPzezv+GXU9wdHyD74afgeRkfdm7O2v/VIPhdhB993Yyf6/IQYDYRhEzn3Idmtgt+Qvad8YN0/hfU8R9VvseX8fd7vlTC2uDX40PmSeSETHwwBz/d0efAO/hwfDfw92C6pS6cc+OD5SVPxC9fuS9+ENUCYApwGXBd0L0uIj2IOafZJkREREQkWronU0REREQip5ApIiIiIpFTyBQRERGRyClkioiIiEjkFDJFREREJHIKmSIiIiISOYVMEREREYmcQqaIiIiIRE4hU0REREQip5ApIiIiIpH7f3WO0mR4X75JAAAAAElFTkSuQmCC\n", - "text/plain": [ - "<Figure size 720x576 with 1 Axes>" - ] - }, - "metadata": { - "needs_background": "light" - }, - "output_type": "display_data" - } - ], - "source": [ - "from scipy.signal import find_peaks\n", - "\n", - "x=ds[\"Counts - Charge distribution\"]\n", - "\n", - "peaks, _ = find_peaks(x, distance=200)\n", - "\n", - "plt.figure(figsize=(10,8))\n", - "plt.plot(peaks, x[peaks], \"or\", markersize=12) \n", - "plt.plot(x, \"-b\")\n", - "plt.xlabel(\"Canal ADC\")\n", - "plt.ylabel(\"Cuentas\")\n", - "plt.grid()" - ] - }, - { - "cell_type": "code", - "execution_count": 11, - "id": "94385c3f", - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([ 131, 361, 599, 821, 1056, 1278, 1520, 1750, 1950, 2172, 2387,\n", - " 2605, 2862])" - ] - }, - "execution_count": 11, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "#Posición de los picos\n", - "peaks" - ] - }, - { - "cell_type": "code", - "execution_count": 12, - "id": "b1915eef", - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "distancia entre el pico 1 y 2 = 238\n", - "distancia entre el pico 2 y 3 = 222\n", - "distancia entre el pico 3 y 4 = 235\n", - "distancia entre el pico 4 y 5 = 222\n", - "distancia entre el pico 5 y 6 = 242\n", - "distancia entre el pico 6 y 7 = 230\n", - "distancia entre el pico 7 y 8 = 200\n", - "distancia entre el pico 8 y 9 = 222\n", - "distancia entre el pico 9 y 10 = 215\n", - "distancia entre el pico 10 y 11 = 218\n", - "distancia entre el pico 11 y 12 = 257\n" - ] - } - ], - "source": [ - "dist=[]\n", - "for i in range(1, len(peaks)-1):\n", - " d=peaks[i+1]-peaks[i]\n", - " dist.append(d)\n", - " print(\"distancia entre el pico\", i, \"y\", i+1, \"=\", d )" - ] - }, - { - "cell_type": "code", - "execution_count": 13, - "id": "61992f52", - "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "promedio de la distancia entre picos: 227.36363636363637\n" - ] - } - ], - "source": [ - "print(\"promedio de la distancia entre picos: \", sum(dist)/(len(peaks)-2))" - ] + "outputs": [], + "source": [] } ], "metadata": {