{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# CEvNS Examples\n", "\n", "Examples for using the `CEvNS` code. " ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "from __future__ import print_function\n", "\n", "%matplotlib inline\n", "import numpy as np\n", "import matplotlib as mpl\n", "#matplotlib.use('Agg')\n", "import matplotlib.pyplot as pl\n", "\n", "#Change default font size so you don't need a magnifying glass\n", "mpl.rc('font', **{'size' : 16})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Import the CEvNS module" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import CEvNS\n", "#help(CEvNS.xsec_CEvNS)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Differential Cross Sections\n", "\n", "Calculate and plot the differential cross-section for a fixed neutrino energy (3 MeV) and a few New Physics models, involving new vector and scalar mediators" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Max. Ge recoil energy (E_nu = 3 MeV): 0.26468447404989304 keV\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa4AAAEqCAYAAABN1LRUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsnXd4lMX2xz9n0zshJAQIvXekCSIIgkqvImLDBl6Vq17kZ7votetVryKKimIXRRTpiAoICFKUKr33QCBAet/5/TEbDCGQTbKbN5vM53neZ3fnnXfmu2+ye3ZmzpwjSikMBoPBYPAUbFYLMBgMBoOhKBjDZTAYDAaPwhgug8FgMHgUxnAZDAaDwaMwhstgMBgMHoUxXAaDwWDwKIzhMhgMBoNHYQyXwWAwGDwKY7gMhnKCiGwTke5W6ygriMgrIvKI1TrKEiKyTkSaW62jpBjDVYEQkVtE5E8RSRaRWBH5UUSudpw7KCJpjnO5x7tOtntQRE6KSFCesntFZFme11eLyO8ikiAiZ0RklYh0cPmbLCEicrOIrBWRFBGJczx/QETEam2FoZRqrpRaVhp9Of7mvUqjr+L0LyKRwB3AlDxl4SKi8v2PJ4vIw27Q95XjM5YoIrtF5N5C6h8UkUwRqZKvfJNDc51Crv9JRJ4voHyQiJwQEW9H0RvARfU8DWO4KggiMg6YCLwMVAVqAe8Bg/JUG6CUCs5zjC1CF95AgV8AIhIKzAfeASoDNYDngIwivxE3IiKPAm8DrwPR6Pv0D6AL4FuM9rwLr1XxKKX7ciewUCmVlqesDXAm3/94sFLqbTf0/wpQRykVCgwEXhSRdoVccwAYmftCRFoCAU729xlwewE/sG4Hpimlsh2v5wI9RKSak+2WTZRS5ijnBxAGJAPDL1PnINCrgPIngO/zlb0NTMp37RPAGaCSo+xeYJnjeXvgXBH0NgWWAeeAbcDAfH2NB7YACcC3gP9l2noC2AckAduBIZe5RynAsEK0VQdmAqfQXzQPFXAfH3foy0Ab9IPA/znKUoCP0UbxR4euxUC4M3ov9/7z/g0BBTTIc91nwIv52nFKUwH34EvADqQ5/q8ec1J3/vvSFtjouOY7x3t5sbD7XFD/BWhcCtyWr+xfwG8WfP4aA7HATYV8/iYAf+QpewP4t+NvWedy9wVt4BKAbnmuDwfSgdb5+voFGFXa98Gl99RqAeYohT8y9AayAe/L1Dn/pZevvDaQCoQ6Xns5PoSd8l8L/JDniyev4QoF4oHPgT6X+kJ01PUB9gJPoUc51zq+2Brn6Wud4wNcGdgB/OMy7Q131LUBI9Bf0tWKeY9swHrgGYe2esB+4IZ892ITUBMIyFO2Bm0YagBxwAbgCsDP8SX7H2f0Xu79U3TD5ZQmZ/5fnNR9/r447t8h9CjdBxgKZAIvFuE+X/T/muf8KaBDvrIvgCnF/AzNR/+QKuiYf4lr3kN/dpTj3gYXdj+BXegfbl7AEfTnTwF1CrsvwEfA1Dxt3gdsKqCvScCb7v7ecedhpgorBhHAafX3dMGlmC0i5/Ico5VSh9AfusGOOtcCqUqpNQVc/wzwT8f6wnmUUonA1egP4EfAKRGZKyJVC2ijExAMvKqUylRKLUV/aYzMU2eSUuq4UuoMMA89BVQgSqnvHHXtSqlvgT1AxwKqViHfPXKsyZ1zrP11AzoAkUqp5x3a9jvez8352pqklDqiLpymekcpdVIpdQz4DVirlNqolMoAZqENhrN6nX7/heCUJmcogu7c+9IJPeqapJTKUkr9gDbI4Px9vhyV0D948tIGPZ2W93/8TSffX3+lVKVLHP0vcc0DQAjQFf2jzpmp8S/Ra3PXATuBY3nOFXZfPgeGi0ju9OIdjrL8JKHvj8di5uArBvFAFRHxLsR4DVZKLS6g/Gu04fgCuMXx+iKUUltFZD562mhHvnM70OsOiEgT4Cv0mtvIfM1UB44opex5yg6hRwW5nMjzPNVxTYGIyB3AOPQvVtBGsUoBVS+6R0qpqxxtHEX/2q0NVBeRc3mu80J/6eflSAHtn8zzPK2A18FF0Ov0+y8EpzQ5g5O6896X6sAx5RgC5Dvv7H2+HGfRRiNXnx96JNNZKfVnEdopEUqpHGCliNwG3I8e7VyOL4EVQF305y0vl70vSqmVInIKGCQi69CGbmgBfYSgR4oeixlxVQxWo+e6BxdW8RJ8B3QXkRhgCJcwXA7+A4zmQkNzAUqpneipqxYFnD4O1BSRvP+btbjwl6dTiEht9C/SsUCEUqoSsBUoyENwNfoX8aACzuVyBDiQ79d2iFKqb756xUpyV0S9hZEKBOZ5HV0cTZfh/Hssgu689yUWqJHPmaCm49GZ+1zYPd4CNMrzuoXjmr/yVxSRdiIy1fH8SRG5poA6PxbgjZh7/FiIFtCDhPqFVXLMcBwA+qJHaXlx5r58gR5p3Q78rJQ6ycU0BTY7obnMYgxXBUAplYCexpssIoNFJFBEfESkj4i85sT1p9DOEp+iPzg7LlN3L3qR/aHcMhFpIiKPOgwfIlITPdIqaLpxLXp95DGHxu7AAGC6c+/2AoLQX1anHP3eRcHGEqXUObSn43sicqOIBIuITUTaONoBPZWVKCKPi0iAiHiJSAsXuvU7rdcJNgG3ODT2Bi76Mi4hJ9FrLFA83auBHGCsiHiLyCD+nlp05j7n7b8gFnLhe74C2OaYBs3Pcf7+obUPaJK/glKqj7rYGzH36JO3rohEid5WEezQfgP6/33pZfTm5R7gWqVUSr5yZ+7LF+i1stEUME3oGHm2QztoeCwV2nCJyBIR2ezYK/Gb40sq7/neovc9bRGRNSLS2iqtJUUp9SZ6KmcC+gvmCPoX8uw81ebl+yU5K8+5r9EfiMuNtnJ5nr+/7EHPqV8JrBWRFLTB2go8WoDOTLT7cB/gNHqB+w7HKK1IKKW2A/9Df0meBFoCqy5T/zX0PXoM7axwEr0P6HHgd8e0zwD0WskBh76paI/EElNUvYXwMFrrOeBWLvw7u4JXgAmOaau+FFG34+88FP0lfQ64Db2WmeHkfT7fv4iML6CLL4C+edZ72gAt8/1/J4lImENzlKPeQGBlEe5DgW8PPS14FD1l+QbwiFJqjlMXK7WvoOlMZ+6LUuog8Dv68ze3gOYHop2mjhfh/ZQ55MIp5oqFiIQ5RiOIyBC0F1Ubx+twtHfb1UqpHSJyFfChUqq4v4ANBsNlEJG1wAdKqU9d1N7LQJxSaqITdTcBk4EYpdR/XNF/WcRxj+9RSm21WktJKDMjLhGJEZF3RGS1iKTKZXaLi0hNEfledBSGRBH5QURqFbXPXKPlIDTf6fpAfO60mFLqd6CWiLQtaj8Gg+FiROQaEYl2TBWOAloBi1zVvlLqKWeMloN0oFl5NloASqkrPd1oQdnyKmwA3ITep/AbcH1BlUQkED1XnAGMQg/LXwR+FZFWBcwLXxYRmYaeC7ejpzxy2QNUFpEuSqlVIjIQ7Y1TB+0ebjAYSkZjYAbaA3EfcKNSKtYKIUqpTlb0aygeZWaqUERsuS7QouN6fQTUdczZ5q33MPAmekPqXkdZXbShecyxloOIbEB7oxXEFUqpC1yWHX0OUUr1y1PWDb1eE4Ke9+4FPKGUmlfCt2swGAyGYlJmDFdeCjFcS9AhbrrkK18OoJQqlveUwy03A73bP76A837o/TMdcg2mwWAwGEqfMrPGVQSaoz3S8rMNaOZsI6IjRecNNDkM7Ul2Jk+dvOefBpYao2UwGAzWUpbWuJylMtrFND9n0EElnSUc+FZE/NHrW3FA/3w7+Z8Xka7o+7Qa7bpbICIyBhgDEBQU1K5Jk4u2ghgMhgrIrl27AGjcuLHFSso+69evP62UiiysnicaLih413yRogs44nxdduOoUmp0Edr7EPgQoH379urPP0stqozBYCjDdO/eHYBly5ZZqsMTEJFDztTzRMN1Fj3qyk84BY/EDAaDwTImTJhgtYRyhycarm3oda78NEPnAbIMERkADGjQoIGVMgwGQxmiVy/LEkWXWzzROWMu0ElEzscpc2xU7kLBIU5KDaXUPKXUmLCwYkYA+vMT2DYbMpJdK8xgMFjGpk2b2LRpk9UyyhVlasQlIjc6nuamuO4jOkz/KaXUckdZbhTqOSIyAb3e9QI69t6U0tTrUpSC1e9B/B7w9of6PaHpAGjcGwKK4nNiMBjKEo888ghg1rhcSZkyXOj0GXl5z/G4HOgOoJRKEZFrgbfQuWsEWIIOYum5QxUReGANHF4NO+fDjnmwawHYvKFOV2g+GJoNMkbMYDBUeMrkBmRPJc8a1+g9e/aUrDGl4PgGbcC2z4Uz+8DLFxr1htY3Q4PrwNvXJboNBoP7MF6FziMi65VS7QutZwyX63G5O7xSELsJNn8LW7+HlFMQUBlaDIN2oyC6pev6MhgMLsUYLudx1nCVtalCQ0GIQPUr9HH9C7DvV9gyHTZ+CX98BDEdocM90Gww+PhbrdZgMBjcihlxuYFS24CcegY2f6O9EeP36lHYFbdC+3ugcl33928wGArl999/B+Cqq66yWEnZx0wVWoBL17iKglJwYDn88THsXAAoaDoQujwENdoVernBYDCUBYzhshBLQz4lHoe1U+DPTyEjQXskdnkYGvTSU44Gg6FUMSMu5zGGy0LKRKzC9ETY8DmseR8Sj0F0K+jxb2h0gzFgBkMpYpwznMdZw+WJkTMMzuAfClf9Ex7aBIPeg4xE+GYETO0Fe5fo6UWDwWDwQIzhKu94+2qHjbF/woBJkHwSvhoKn/aBw2usVmcwGAxFxhguFyIiA0Tkw4SEBKulXIyXj97z9c/10PcNOHsQPrkBZozSzw0Gg8FDMIbLhZQ4yG5p4O0HHUdrA9b9SdjzM7zbEX75j14XMxgMhjKOcc5wA2XCOcNZEo/Dkuf1frDgaOjzqt7IbBw4DAaXkBsZvk2bNhYrKfsYr0IL8SjDlcux9TD/XxC7WcdB7PcGhNexWpXBYKhAGK9CQ9Go0Q7uXQq9X9UR6id3gpVvQU621coMBo9m8eLFLF682GoZ5Qoz4nIhlkXOcDUJx+DHx3R6lRrtYcgUqGKyOhsMxcHs43IeM+KyAI9wznCGsBpw8zS48ROdTuWDq3U0DrvdamUGg8FgDJfhMrQYppNb1u2qR2BfDtLOHAaDwWAhJq1JGWLw5FXsi0sG0WmdbTZBAJE8j7nn8jw/Xy4g6Oc2xzXkeZ73vOQpK7Atx3ObgKh/0TOsKaMOTiFzYifeC3+MvwLaX9SWTf6+xssmeNkEm+hHLxFseR9tXFCWW/+Ca84/p4CyC6+z2QQfm+DjZcPbS/D1suHtZcPHS5ddWO4os9nw8Ra8bbqeGE9Kg8EjMIarDHFD82jiaqWfj8aklMKuQKFQCpSjTCkdscmulKMsT508ZbnnyX2ep57dcSK33b/P5+njfF3F0sA+7PBtwbizL/FY/ARmBt3Mt4G3kI2Xow91vh27XT9m2xV2uyJHKXIueK7P5+Q/77jGqmVX74sMn37t523D38cLfx+vPM9t+HnnfSz4XO7rQF9vgny9CfLzItjPm0A/bwJ9vLDZjLE0GIqKcc5wAx7pDu8smamw8P9g01c68vywqRAS7dIulMOQ5TiMYH7Dl9fY5R52pcjKUWTnKDJz7GTn2MnKUWTZ7WRl28m2K7JyyxznM3OUo16ecrsiM9vuqKPbysy2k56VQ3p2DhlZdtKzc0jP0mUZjnMZWXYyc4q2BigCgT5eBPp5E+ynjVqgb+5zb4L9vAj19yE0wIcwx1Ep8O/nYQE+hPj74GWMX5lm165dADRu3NhiJWUfs4/LQsq14cpl09cwfxwEhMPIr3V25gpOjsPoFWTkUjOySc7IJjUzh+SMbFIysknJzCElI5vUzGySM/RzXZ5NSoaul5iWRUb2pQ2iCIT4eRMW6EN4oC8RQb5UCfajSoiffgx2vHY8Dw/0NaM8Q5nFWcNlpgoNxaPNLRDdEr4ZCZ/0gcGTtTNHBcbLJgT4ehHg6+XSdtOzckhIy/r7SM3i3AWvM0lIy+JsahZxSRnsiE0iPiWDrJyLf5TaBCKC/YgO9adamD/VKwVQLcyfapUCqO54rBrih7eX8dtyFfPmzQNgwIABFispP5gRlwspN/u4ikLyKZhxu9603O0xHf/QZr70rEYpRUJaFqeTMziVlEl8SgankzI4nZzJqaQMYhPTiT2XRmxCOskZF24ytwlEh/pTOyKIOlUCqRMRdP557cpBLjfM5R2zj8t5zIjLApRS84B57du3H221llIjOBLumAsLxsGK13Sk+UGTdToVg2WICJUCfakU6EuDqMvXTUzPIvZcOscT0jiRoA3akbNpHIpP4adtJzmTknlB/ehQf+pWCaJxdAiNo0NoVDWERlWDCfH3ceM7Mhj+xhguQ8nx9oWB70Dlujpgb+ppuOkL8AuxWpnBCUL9fQiN9qFxdMF/r4S0LA7Hp3IwPoVD8SkcOJ3KvlPJzPjzCKmZOefr1agUQOPoEJpXD6VVTCVax4QRFepfWm/DUIEwhsvgGkSg66M6wvzcf8Jn/eHW7yC4kJ/7hjJPWIAPLWPCaBlzYUQYu11x7Fwau04ksetkErtPJrHrRBLLd58iR++3IDrUn9Y1w2gVU4m2tcK5olYl/H3MVKOhZDhluETkcDHaVkA/pdTWYlxr8FSuuBWCqugElZ/cAKPmQViM1aoMbsBmE2pWDqRm5UB6Nat6vjwtM4dtxxPYfDSBLUfPseVoAj9tOwmAr5eNNjUrcWW9ynSqF0HbWuFmzcxQZJxyzhARO7AQOOVkuzbgNqCDUmpD8eV5JhXCHb4wjqyDr4Zpd/k750OlWlYrMlhIQmoWfx46w9oDZ1izP56txxKwK23IOtatTPfGkVzbJIq6VYLKXQSTI0eOAFCzZk2LlZR9XLqPy2G4Oiml1jnZuTeQCbQ3hqsCc2w9fDkE/MJg1Fy9BmYwAEnpWfx56Cy/7z3Nsl2n2BOXDEDtiECubRJFv5bVaFsr3Ow5q2C42nD9G/hEKRVbBAH/BqYqpU46e015wRiuPBzfBF8MAt8gPW0YUd9qRYYyyJEzqSzbFcevu06xcu9pMrPtVA/zp1+ravRvVZ1WMWEeOxL79ttvARgxYoTFSso+JnKGhRjDlY8Tf8HnA7XxuvsnnTbFYLgESelZLN5xkvmbY1mx5xRZOYo6EYGM6FCLG9vFEBniZ7XEImH2cTmPq0dcnwCfKaVWuEJceaVCbkB2luMb4bMBEFoN7loEQRFWKzJ4AAmpWfy0/QQz1x9l7YEzeNuE65pV5ZYra3F1gyoeMQozhst5XG24UgB/4DDwBfClUmpviVWWU8yI6xIcXAVfDYXIJnra0D/UakUGD2JvXDLf/nGY79cf5WxqFo2rhjCmWz0GtK6Or3fZjdZiDJfzuDoDclXgXuAgMAHYJSIrRWS0iHh4ul9DqVGni96YfHKrjnGYlW61IoMH0SAqmH/3a8aap3ryv+GtAXj0u810e+1XPl55gPSsnEJaMJQXnDJcSqlkpdSnSqkeQB3gaaAyMAWIFZFvRKSPiJTdnz2GskGjG2DIFDi0EuY8gGXJtwwei5+3F8PaxbDoka58elcH6lQJ5IX52+nxxjKmrztMdhHTyxg8jxI5Z4hIR+AO4CYgAogDpimlxrtGnmdipgqdYOVbsPhZuOZx6PGU1WoMHs7ve0/z3592sfnIOepVCeKpvk0v2BRtJadPnwagSpUqFisp+5SqV6GI+ACvAP8CUEpV6K3wxnA5gVIwdyxs/AqGfAitjauwoWQopfhl+0le+2kXe+OS6dkkimcHNqdm5UCrpRmcxNVrXJfqpIGIPA/sQhutJOCTkrRpqCCIQL+3dBbluWPh0GqrFRk8HBHh+ubR/PhwV57q24TV++Pp9eZyPli+73zsRCv47LPP+OyzzyzrvzxS5BGXiFQGRqCnCDuiYxIuBj4HZimlKvyKuxlxFYHUM/DxdZB2Fu5bYeIaGlxGbEIa/5mzjZ+3n6R97XDevKkNtSJKf/RlvAqdx6UjLhHxEZEhIjILOA5MBkKAJ4FaSqneSqlvjNEyFJnAyjByOmRnwow7IDvDakWGckK1sACm3N6ON29qza4TSfR+ewUz1x+1WpbBBTg7VXgC+B64GvgI6KiUaqGUek0pddxt6gwVgyoNYcj7OrbhoiesVmMoR4gIQ9vGsOhf3WgVE8aj323mmTlbycw2noeejLOG6zfgRqC6UuqfSikzD2ZwLU0HQJdH4M9PYOM0q9UYyhk1KgXw1T1XMrprXb5YfYhbPlpDXJKZIPJUnN3HNVgpNUspleVuQYYKzLVPQ91usGAcxO20Wo2hnOHtZePf/ZoxaeQVbDueyJDJv7PvVLLVsgzFoFhehSJyhYj8ICKnRSRbRNo6yl8Wkd6uleg5iMgAEfkwISHBaimeiZc3DJ0KvsEw8x4TWcPgFga2rs6M+zqTnpXD8A9Ws/nIObf2t3DhQhYuXOjWPioaxfEqvBrtRbjf8TgWR94tEXkRaKGUGuxypR6E8SosIbt/hq+Hw5X3Q59XrVZjcCdZaXD2IJw5AAlHtHdp2jlIP6fP5eLlC/5h+giqAuF1/j58AorV9YHTKdz+8VrOpGTy5T0daVe7sgvekKEkuG0DsoisBOKBwYAXeRJGishQYKJSqkKnuzWGywX8+Dis/QBu+Q4aXW+1GoMryEzRDjjH1sOxDTpXW8Lhi+v5hWoD5RMAOKK/Z6dDeoI+yPOdJTaIbAo1roCYDlC/J1RyPtPwycR0bv5wDaeTMpg2+kpaxVQq0VssiPfeew+ABx54wOVtlzfcabhSgaFKqUUi4gVk8bfh6gb8pJQq3k+gcoIxXC4gKx2m9oTkk/DAGv0r2+BZKAUnt8HexbBvCRxeAzmZ+lx4HajeFqKaQnhdnR27Um0ICNdTxpfCbofUeDh3SI/UTu3UKXOObYC0M7pOZFMdE7PVTVC1eaEyj59L46Ypq0lKz+bb+zrRJNq1WQvMPi7ncdZwXeY/5JKkA5faxVcNMAs8hpLj4w9DP4Ip3eDHx+BGE5DFYzi9F7Z+D399B/GO7EdRzeHK+6Bud6jRVu/fKw42GwRH6iMmz/ebUnB6N+z5Bfb8DKvfhVUTIboltLsL2txyySnF6pUC+GZ0J4Z/sJq7P/2D2WO7EBXiXzx9hlKhOCOuuUAloIejKAtop5TaKCI/A6eVUre4VqZnYUZcLmT5a/DrSzBiGjTtf75YKUVSVhJxKXGcTj9NUmbS+SMxM5GkzCQycjLItmeTbc8mR+WQbc/Gruz42Hzw8fLB1+aLr5cvPjYfgn2DCfMNI8zv7yPcL5yowCj8vc2XWKFkZ8C2WbDuIzj2JyBQ52poMUyPfkKrl66elHhtPDd+BSe2QGAEdLwPOj8AfiEFXrL1WALDP1hN4+gQpo/phL+Pa0KumhGX87hzqrA1sAqdm+t7dIqTd4DWQDugg1JqV1EFlyeM4XIN6dnpHDq3j0M/3M2hrAQONh9AbEY8calxxKXGkZadVuB1NrER7BOMv7c/3uKNt80bL5sX3jZvBCHLnkVWThaZ9szzj6lZqSgK/iyE+4UTHRR9/qgVUou6YXWpF1aPqkFVsVXkbD7JcfDHVL3/LuUUVGkEbe/QBqu0jVVBKAWHfoff34HdP0JgFej+BLS7E7x8Lqq+aOsJ/vHVega1qc7EEW1ckmHZGC7ncWt0eIf7++tAN7SDhh29SXmcUmpjkRssZxjDVTSUUhxPOc7O+J3sOLODnWd2svvsbmJTYi+oFyU+1IhsQVRgFFGBUVQNrEpUYBRVAqoQ5hdGqG8oIb4hBHoHFvkLJ8eeQ3JWMgkZCSRkJHAu4xxnM85yIuXE+SM2JZYTKSdIzvp770+AdwB1QutQr1I9mlZuSrOIZjSt3JRg32CX3JsyS3IcrHob/vhYO040ukFPBdbroQMol0WOrodfntG54KJbwsB3oXqbi6pN/nUvr/+0i5eGtODWK2uXuFtjuJynVNKaiIg/OqHkOaVUqqMsWClVoXf1GcN1eVKyUth8ajOb4jaxMW4j2+O3k5iZCOjRUt3QujSq3Ii6YXWpG1qX2qG1qb3xWwJXvgW3z4b6PQrpwX0opTiTfoYDCQfYn7CfAwkHOJBwgD3n9hCXGne+Xp3QOjSNaErryNa0q9qOhpUa4mUrB9l+MpLgt//Bmg8gJwNa3QzdxkNEfauVOYdSsGMeLBwPKafh6n9B9ycvcAix2xWjPl3HugNnmDv2ahpHFzy1aHA97pwqnKSUeugS54LRXoVditRoOcMYrgtJykxi3Yl1rI1dy6a4Tew6uwu7siMIjcIb0TKyJU0rN6VJ5SY0DG9IgHcBi+hZ6fBeJ7B5wf2/g7df6b+RQjiddprt8dvZHr+dHfE72Ba/jZOpJwEI9gmmTVQb2lVtR/uq7WlRpQXetuL4RlmE3Q6bpsGS5yElDlrepKfcPMVg5SftLCx6CjZ/rVPr3PgJBEedP30qKYM+b/9GlWBf5o69Gl/vCjwdXIq403AlAa8opV7OVx4MLEJHizf7uCqw4cq2Z7P19FZWx65m9fHVbDm1hRyVQ4B3AK2qtKJNVBvaRrWlVWSrok2p7fkFpt2oQ0N184wk28eTj7P+5Ho2xG1gw8kN7E/YD0CITwidqnfiqupX0aV6F6oFV7NY6WWI2wFzxmqni5iOelN4jXZWq3ING6fpEGMB4TpLQZ6pw8XbT3LvF38y7rpGPNSzYbG7eOONNwAYP94z/metxJ2GqzcwB/iHUupTR1kg8BNQD+imlNpXdMnlh4pouJIzk1l5fCVLDy9l5bGVJGUmIQjNI5rTuXpnOlfvTJvINvgUsCBeJKbfCnuXwINrIbzk6w+lzZn0M/xx4g9+P/47q46tOj8iqxtWl2tirqFnrZ60imxVNhw+sjNh5Vuw4nXwD4UbXtF7o8rqGlZxid0C02/Ro7Cbp0G97udPjf16Az9tO8HCh7rSsGrxpgzNGpfzuNs54w50epNhwBL0SKshcI1Sak+RGyxnVBTDdTrtNMuOLGPJ4SWsjV1Llj2Lyv6V6RbTjS41utApuhOV/F0cieDcEZjcUTs/3pHYAAAgAElEQVQBjPzatW2XMkop9ifsZ9WxVaw8tpI/Tv5Btj2byIBIetTsQc9aPekQ3aHkxr44nN4L39+lXclb3Ah9/lu+N4EnHoevhul9Zzd/Aw17AXA6OYNeby6nSXQI34zuVCwvQ2O4nMftzhki8hjwLLAVqAV0V0qZkN6Ub8OVkJHAL4d+YcH+Baw/uR6FokZwDXrW6sm1ta6lTWQb9zsh/PYmLHnOckcNV5OYmchvR39jyeElrDy2krTsNEJ8QuhRqwf96vXjyugrS8fBY/N0mD9OryMOehea9HN/n2WBtLPw+UC9kfnW73SmAuDL1Qd5es42Pry9Hdc3jy5ys8ZwOY9LDZfIJectJgI3Az2BbbmFSqkKnaWtvBmutOw0lh9dzoL9C1h5bCXZ9mzqhNahT90+9KzVk0bhjVyy38VpstLh3Q4QEAZjVuhoCuWM9Ox0Vh9fzeLDi1l6eCnJWclEBkTSu25v+tXrR7PKzVx/z7PSYcGjsOkrqN1FRy4Jq+HaPso6KfHwWT8d8PfunyC6Bdk5dnq//Rs5dsVPj3QrsqOGMVzO42rDZYdL7M7UUTDznlNKKQ9yl3I95cFwKaXYfGozM/fM5OeDP5OanUpUQBR96vahb72+NK3ctHSNVX62fAc/3AuDP4A2I63TUQqkZ6ez4ugKFuxfwIpjK87/cBhYfyCDGgwiKjCq8EYKI+kkfHsrHP0Duj2mPQbLg/t+cUiMhY96gM0HRi+F4Eh+3RXHXZ/+wXMDmzPqqjpFaq5Pnz4A/Pjjj24QW75wteF6lksbrotQSj3nbN2ygIjcBXwCDFFKzc5TXh/4HIgCUoDRzmR/9mTDdSb9DPP2zWPmnpkcSDhAoHcgN9S5gf71+tOuaruysxfJbtdfLimn4Z9/Fju1haeRO1U7f/981p9cj5d40TWmKzc2vJEuNboUz8U+djN8M1JPlQ35AJoNcr1wT+P4RvikD1S/AkbNQ9m8GDFlDUfOprL8/3oY93g3USobkMsDIlIb+AY9cvxvPsP1C/CtUmqqiFwHvAs0UYXcNE8zXHZlZ/Xx1czcM5Nfj/xKtj2b1pGtGdZwGDfUuYFAn0vFVLaYA7/B5/2h53+g6zir1ZQ6hxIP8cOeH5izdw7x6fFEBUYxuMFghjYcSo1gJ6f49i/XHnX+lbSzS7XW7hXtSWyZAT+M1huUuz/B8t2nGPXJOl4d2pKbO1boHT9uw9Ujrs5KqdUuUXbpPmKAx4H26LiHAUBdpdTBAurWBN4CrkMbnMXAI0qpApL7XLZPG/Czo9//oXOJzXaci0Qny6yslMpylO0Cbi1s1OUphispM4nZe2fzzc5vOJJ0hEp+lRhQfwBDGwylQXgDq+U5x9cjdCy6hzcXP+K4h5Nlz2LFkRV8v+d7Vh1bBcDVNa7m1qa30rl650u71u+Yrz0HIxrAbT9AaBneS2YVP9wHf82AOxeianVi0ORVnEvNYumj1+Dt5dyo64UXXgDg6aefdqfScoGzhsvZ8e5KEYkVkSki0ltE3OGf2wC4CTiLjntYII49Y0uBJsAo4Ha0K/6vIhJUxD7HAauUUusLOFcLOJ5rtBwccpR7NPvO7ePFNS/S87uevPbHa0T4R/Dfrv9lyfAlPNbhMc8xWgC9ntVhiH6fZLUSy/Cx+dCzdk/e7/U+P9/4M/e1vo/t8dv5x+J/MGj2IL7Z+Q0pWSkXXrTpa5hxux5h3bnAGK1L0fd1CKsJc8ciOZk82KMBh8+ksmjbCaebWLJkCUuWLHGjyIqHsxPiNYAhwCBgNpAhIouAWcBCpVSiC7SsUEpVBRCRe4FLpb0djd7o3FgptddRfwuwB7gPeNNRtoFLG5krgFDgRqBrETR67M5Lu7Kz4ugKpu2YxprYNfjafOlTtw+3NL2FZhHNrJZXfKKaQvMhsPZD6Dy2fO81coLooGgebPMgY1qO4adDPzFt+zReXvsykzZMYnCDwYxsMpJaB1bD7Aeg3jU6XYxfOQ8IXBL8Q6H/m3qP16q36dX1/4gJD2DamsP0b1UGot9XUJwyXEqpE8D7wPsiEgL0RxuxDwB/EVmONmJzlVLHiyOkCC70A4E1uUbLce0BEVnl0PSmo6zt5RoRkfuB2sAeh3dcNPChiFRTSr0PHAaqi4hPnlFXbUe5x5CZk8mC/Qv4dNunHEg4QNXAqjx0xUMMazSMyv7lZGqt+xM6F9Tvk+C6561WUybw8fKhf73+9K/Xny2ntjBtxzSm75rOtB3T6JWSwt11O9Ji5PQK49RSIhr0guZDYcUbeLUYxsiOtXj9p13sjUumQZQx+lZQZNcYpVSSUuobpdTNQCR6JHYAnZfriIisFZEnXKwzL83Rm57zsw1weuiglHpfKVVNKVVHKVUHWAOMcRgtlFKngHXAnQAO5wwBCppWRETGiMifIvLnqVOnivB23ENSZhKfbP2E3jN788zvz+Dn5cd/u/6XH4f9yOhWo8uP0QKIbAwtb9RJDJOtv/dljVaRrfhvt//y8xVPcW9CImsCgxhJLHcvfZCVx1ZS0R20nKL3Kzp/15Lnual9TbxtwjfrPOo3bLmiRD6dSqkspdSPSql/KKVqAFcDy4A7XCHuElRGr4Pl5wwQ7uK+/gHcJSK70fnHbr2UR6FS6kOlVHulVPvIyEgXy3CeuNQ43lz/Jtd/fz1vrX+L+pXqM+W6KczoP4O+9friY7MgfFBpcM3jOi9UBV7ruixH1xM5eywP+dXml2E/Mb79eA4lHuL+xfczfN5wFuxfQLY922qVZZeQaOj0AGyfTWTyTm5oEc3364+SnpVT6KURERFERESUgsiKQ5l0h3escX1EAV6FIpIJ/E8p9WS+8peAx8vC5mcrvAqPJx9n6l9Tmb13Njkqhxtq38CdLe707PWrovL93bD7Z/jXVghwcYxET+bcYfiop54WvHcJBOsfVlk5WSw4sIBPt37K/oT91AiuwV3N72JIwyH4evlaLLoMkp4AE1tBzY781vE9bv94He/f2pY+LY1ji6twiVehiNQVkSUisl9E3nQkjsw9t84VQovBWfSoKz/hFDwSKzVEZICIfJiQkFBqfR5NOsqzvz9Lvx/6MXvvbIY2HMqCIQt47ZrXKpbRAujyMGQm6TTyBk16Aky7CbIzdPy94L9nA3y8fBjcYDCzBs1iUo9JRARE8OLaF+n7Q1+m75xOZk6mhcLLIP5hcPUjsOdnOvvup0qwL/O2FGtJ31BCCpsqfA/4ARiONhZLHM4ZAFbNOW1Dr3PlpxmwvZS1XIBSap5SakxYWJjb+zqSdIRnVj3DgFkDmLdvHsMbD2fh0IVM6DSBmJAYt/dfJqnWWkeNX/uBjrtX0bHn6FFo/B4Y8YVeCywAm9joUasHX/X5iinXTaF6cHVeWvsSfX7ow9c7viYjJ6OUhZdhOo4B/0p4r3mXfi2rsWRHHEnpWZe95Mknn+TJJ5+8bB1D0SjMcFVVSk1WSq1XSt0JzEcbrzCKEALKxcwFOolIvdwCEakDdHGcK9ccTjzMhJUTGDBrAAv2L2BEkxEsHLqQp658iuigokeuLndc/Qgkn4Qt061WYj3LXoW9i/VepHrdC60uIlxV/So+7/05H13/ETHBMbyy7hX6zuzLtB3TjAED8A2C9nfBzvncWD+bjGw7S3fGXfaS1atXs3q1W+M3VDgKM1wX5EdXSr0CzEDn4CpeVrXLICI3isiNQG561T6OsmvyVPsIOAjMEZFBIjIQndjyCDDF1ZrKCrHJsTyz6hkGzh7IooOLGNlkJIuGLeKJjk9QNaiq1fLKDnWv0SOvVZN0PMOKyp5fYMVr0OY2aH93kS4VETpV68RnvT9j6vVTqRlak1fXvUrfmX2ZsWsGWfbLjzDKPR3vA7HR4uh3VAn2ZcmOyxsug+spzJFht4hcp5T6JbdAKfWGI1r8G27Q812+1+85HpcD3R39p4jIteiQT1+iXdSXoEM+JbtBk9OIyABgQIMGros8cSb9DFP/msr0nXoEMbLJSO5peQ9VAir2RttLIgJXPQQz79GjjUaX2sdejjl3RMfYq9pCj7aKiYhwZbUrubLalayLXcekjZN4Yc0LfL7tcx5s8yC96/YuG5maS5vQatC4L7LlG3o1HMLCHXFk5djxcTIElKHkXNarUET8AJRSF80RiEgNpdQxN2rzWFzhVZicmcwX27/g822fk56TzsD6A7m/9f1UDza79QslOxMmtoDoVnDb91arKV3sOfD5AB3x/b4VEFHfZU0rpVh+dDmTNk5iz9k9NA5vzENtH6Jrja7Wprixgj2LYdowNnaayJBlUUwf04lO9Qp2eTf5uJzHWa/Cy464CjJYec4Zo+UGMnIymL5zOlP/msq5jHNcV/s6xrYZS71K9Qq/2KDx9tXTY8tegfh9Lv3yLvOsngyHVsGg91z+vkWE7jW70y2mGwsPLGTyxsk8uORB2ka15aG2D9GuarvCGykv1O8BoTG0jJuPr9e9/Loz7pKGKyamgjpLuZFi7+NyRGivCfjnP6eUWlpCXR5NcUZcdmVn9t7ZvLfpPU6mnqRztc483PZhmlcpyIHSUChJJ+Ct5toLrPcrVqspHU5s1TnKGl4PI77S06ZuJMuexaw9s/hg8wecSjtFt5hujGs3jvqVKsgPhV+egdWTuTvya07nBDF37NVWK/J4XB0dPm/D9URkNdpB4jd0SpHFwC95HiskJdnHJQgzds0gKjCKqddP5cPrPzRGqySEREOzwbBxGmRYuvRZOuRkwax/6LxaA952u9ECHZX+psY3sWDoAh5p+wgbT25k6NyhPL/6eU6nnXZ7/5bTfAjYs7k5ZAtbjyWQWIhbvMF1FHnEJSJLgcbAq8BO4KJdikqp5S5R56EUd43rbPpZKvlVqnjrBe7i8Fr45Hro/1aRPes8jpUTYfF/dLT3pv0tkXA2/SwfbP6AGbtm4Ovlyz0t7+H2ZrcT4F1OA/kqBZPacM6/Jm0O3M/Ho9rTs+nFHr6PPPIIABMnTixthR6H20ZcQAfgIaXUO0qpX5RSy/MfxWjTAIT7hxuj5UpqdoSo5rDhS6uVuJezB/WerSb9LTNaoP9/n7zySWYNmkXn6p15Z+M79J/Vnzl752B3OvmDByECzQYRdnI1lb3TWbM/vsBqmzZtYtOmTaUsrnxTHMN1lAJGWQZDmUME2t4OxzfAyW1Wq3EPSsH8cWDzgj6vWa0GgDphdZjYYyKf9f6MqIAoJqyawIj5I1gTu8Zqaa6n4Q2IPZtbq+xn3YEzVqupMBTHcL0MPF6MbMMGQ+nT8iaw+cDGr6xW4h62z4F9S+DapyGshtVqLqBd1XZM6zeN17q9RmJGIqN/Hs1DSx/iSOIRq6W5jpodwS+Unj5b2B6bSEZ24dHiDSWnOPm4vkRvCD4oIvNE5It8x+eul+kZWBFk11AIQRHQpC9snq73d5UnstK1Z1vVFtBxtNVqCsQmNvrU7cPcIXN5uO3DrIldw6A5g3h7w9ukZqVaLa/kePlAve40Tl5LVo6dXSeSrFZUISiOV+GdwJNAJaAt0LWAo0JSmkF2DUXgijsg7QzsWmi1Etey9gM4dwiuf1FPFZZh/Lz8uLflvcwfMp/edXoz9a+pDJg1gPn753t+IsuG1xGQdpKGcozNRy/+0dqoUSMaNWpkgbDyS3G8Cg8BfwL3KKXOuUWVh2NFPi7DZbDnwFstoGozuG2m1WpcQ/IpmHQF1OkCt3xrtZoisyluE6+ue5Vt8dtoE9mGJ658guYRHrr948x+mHQFL8lozjW7ndeHt7ZakcfiTq/CCOA9Y7QMHoPNC1rfDPuW6i/88sCyVyA7TY+2PJA2UW34ut/XPH/V8xxOOszI+SN59vdniU8r2DOvTBNeF4Kr0iNwP1sKGHEZXE9xDNdKoKmrhRgMbqXlcFB22D7baiUl5+xB2PA5tLsLqjS0Wk2xsYmNIQ2HMH/IfO5odgdz9s5hwKwBfLHtC8+KQC8CtTrRPHs7e+KSSM+60EFjzJgxjBkzxiJx5ZPiGK6HgdEicquIRIiILf/hapEGQ4mp2gwim8LWcjBVuOJ1EC/o+qjVSlxCiG8I4zuMZ+agmbSKasXrf77OiPkj2HByg9XSnKfWVYRlxFJVxXPgdMoFp3bv3s3u3bstElY+KY6R2QG0BL4A4oCsfEc5c91yHuNVWMZpOQwOr9ZpPzyVM/th0zc6EkhoNavVuJR6YfV4v+f7TOwxkaTMJEYtGsWElRM4k+4B+6NqdQKgg20Xe+IqQIgxiyksH1dBPI912Y/LNEqpecC89u3bl03f5IpO86Gw9EXY9gN0edhqNcVj+evaBfvqR6xW4hZEhJ61etK5Wmc+3PIhn2/7nF+P/MrDbR9mWMNheJVV78mqzVFevjS3HWTvSeMS726KbLiUUs+6QYfB4H4i6kP1tnq60BMNV/w+2DIdOj2ggwiXYwJ9Anmk3SMMqD+Al9a+xAtrXmDWnllM6DyhbHofevkgUU1pF3eEj82Iy+2Y9ShDxaLljTrJ4um9VispOmveA5s3XPVPq5WUGvUr1efj6z/m1a6vEpsSy8j5I3lpzUskZiZaLe1iolvRWF084mrTpg1t2rSxSFT5pDgbkB8XkXcucW6SiPxfyWUZDG6i2WD9uGOOtTqKSuoZnaKl5U3lfrSVHxGhX71+zBsyj5FNRjJj9wwGzBrA3H1zy9bm5ehWhNgTSI0/SlbO30GFJ06caCLDu5jijLjuArZc4twmx3mDoWwSVgOqXwE7PSyKxh8f631bnR+0WollhPiG8OSVTzK933RigmP498p/c+/P93Iw4aDV0jTRLQFozAEOxacUUtlQEopjuGoBey5xbj9Qu/hyDIZSoHE/OPYnJJ20WolzZKXDug+hQS/t1l/BaRrRlC/7fskznZ9hR/wOhs0dxpTNU8jKsXjvV3QLAJrJIQ7F/x2H8bbbbuO2226zSlW5pDiGKxW4VBjqGCCj+HI8G+MO7yE07qMfd/9orQ5n+es7SImDzmOtVlJmsImN4Y2GM2fwHHrU6sG7m95l+LzhbIzbaJ0ovxBywmrR2HaEw2f+NlxHjx7l6NGj1ukqhxTHcP0G/J+I+OUtdLx+1HG+QmKC7HoIVZtDpVqeM1247kOdELNed6uVlDkiAyN545o3mNxzMqnZqdzx4x28sPoFy5w3bFUaUNcWx5EzaZb0X1EojuF6FmgI7BaRl0TkARF5CdjtKH/GhfoMBtcjoqcL9y+DjDLuunx8I5zYAu3v0roNBdItphuzB83m9ma38/2e7xk0exA/Hfyp1J03pHJ96soJDps1LrdSnHxcm4EewCHgceBdx+MBoLvjvMFQtmncB3IyYP+vViu5POs/A+8AHWvRcFkCfQJ5rMNjfN3vayIDIhm/fDz/XPpPYpNjS09ERH2CSCXpzInS67MCUqx9XEqpdUqpbkAIel0rRCnVXSllcnkYPIPaV4F/WNmeLsxIhr++hxZDIaCS1Wo8huYRzfm639eMbz+edSfWMWjOIL7c/iU59lLITly5HgB+iQfPF3Xu3JnOnTu7v+8KRHFCPp1HKZUGmMlcg+fh5QMNroO9i8FuB1sZ3Iu/dSZkJkPbUVYr8Ti8bd6Maj6KXrV78dKal3jtj9dYsH8Bz131HI0rN3Zfx5XrAxCZeZTUzGwCfb155ZVX3NdfBcWpT6uILBWRJs426ogSv1REPDfngqH806Cn9tY7udVqJQWz/jMd0b5mR6uVeCw1gmswuedkXr/mdWJTYrl5/s28u/FdMnPcFAu8Ui3s4kVt2wlOJKS7pw+D01OF3dHTgs4ixbjGYChd6vXQj/uWWqujIOJ2wPEN0G6UccooISJC7zq9mTNoDn3q9mHKlincNO8mtpy6VByFEuDtS2ZQDerKCWIdhmvYsGEMGzbM9X1VYIoyPzJbRPY7c6A3KJehWCwGQwGEVtNu5vuWWK3kYrbM0Dm3WpgvPFdRyb8SL3d9mck9J5OclcxtC2/jtT9eIzUrtfCLi4CqXI/acpLj5/QqSnx8PPHxHpjZuQzj7BrX58Vs/3Qxr/NIRGQAMKBBgwZWSzE4S/0eep9UZgr4BlmtRqOUdsqo1x2Co6xWU+7IdZ2fuGEiX27/kl8P/8qzVz3LldWudEn7PpVrU/3QBn5LqrCxGNyOU4ZLKWXiDzqBycflgTToCavfhYOroNH1VqvRHFkHCYfh2n9braTcEuwbzIROE7ihzg08+/uz3PvzvQxrOIxH2z9KiG/JVji8w2tSRRI5l2jycrmLMuhKZTCUIrU6g7d/2Vrn+us7ralJP6uVlHs6RHfg+4Hfc2fzO5m1dxaD5wxm+ZHlJWs0tDoAOeeOuUChoSCM4TJUbHwCtNfewZVWK9HYc2DbLGjUG/yMb1NpEOAdwKPtH2Va32mE+oYydulYHlvxGGfSzxSvwTAdylWStOHq2bMnPXv2dJVcA8ZwGQxQp6t2iU87a7USOLwGUk9D88FWK6lwtKjSghn9Z/BA6wf45dAvDJ49mB8P/Fj0sFGhMQD4puiIHU8//TRPP/20q+VWaIzhMhhqdwEUHFpttRLYOR+8/PTmaEOp4+Plw/1t7mdG/xnUCK7BYyse4+FfH+Z0WhH8zEKrAeCfHucmlQZjuAyGGu20sbB6ulAp2DFfezr6BVurpYLTMLwhX/X9ikfbPcqqY6sYPGcwC/cvdG705RtEpi2AoKwzZOfY6dOnD3369HG/6AqEMVwGg4+/Xuc6ZLHhOrFFexM26W+tDgMAXjYv7mxxJ98N/I7aIbV5/LfHGbdsnFOjrwy/CKpIAonp2aSlpZGWZiLjuZISGy4RqScin4rIhyJSL095SxExE7sGz6B2FzjxF6RbmAR0x3wQ29+JLg1lgnph9fi8z+f8q92/WH50OUPmDGHRwUWXvSbTvwpVSOBcqptCS1VwXDHi+gLwBWKBj0RkkIjsBNYALV3QvsHgfup0AWW3dp1r94/aPT+oinUaDAXibfPm7hZ3892A74gJjuH/lv8f45aNu6TnoT0oUu/lSssqZaUVA1cYrgil1K1Kqf8A9wAzgC+BKKXUTS5o32BwPzEdwOYDR9ZY03/SCT3ia9DLmv4NTlG/Un2+7PslD7d9mGVHljF49mB+OvjTxRWDIqkiCSQYw+UWXGG4zk/eKqUOAgeUUi8ppUwKUIPn4BMA0S3hqEUp5XI3QBvDVebxtnlzb8t7+bb/t1QLrsb45eMZv3w8Z9P/3k7hHVKVyiSRmJxG//796d/frFu6khLl43LQWERmA9sch90FbRoMpU9MB9j4FeRkg5crPhpFYO9iCK6qjafBI8j1PPx066e8v/l9/jjxBxM6TeC62tfhGxqFTRSpifGMHz/eaqnlDleMuLoDi4AoYDxQX0TOishPIvKcC9r3GERkgIh8mJBg4QK/ofjEdICsFDi1o3T7tefoEVf9a00KEw/Dx+bDmFZj+Lb/t1QNrMq4ZeN4bPljpAf6A5CZXMzoG4bLUmLDpZT6Qyn1gVJqtFKqLToH13XAD0C1krbvSSil5imlxoSFhVktxVAcYtrrx6N/lG6/xzfqqB1mmtBjaRTeiGn9pvFgmwf55dAvDNs9lSWBAWQlx9O9e3e6d+9utcRyhbMZkJ3+GaiUylRK/amUmqKUGlN8aQZDKRNeBwKrlP46194lgPyd2NLgkfjYfPhH638wvf90Iv3CeaRqJIvSZ5Btz7ZaWrnD2RFXooiYyXdD+UZETxeW9ojr4G96bSsoonT7NbiFxpUb83XX/3H/2QR2yx62xW8jMSPRalnlCmcNVxAQkPtCRGwisl5E6uatJCL+IhLqSoEGQ6kS0x5O7y69gLvZGdpQ1rm6dPozlAo+QZE8cC6BW5O64CVe7Dq7ixfXvOjybMsVleKucQlwBRCer7w1YFYjDZ5L7jrXsQ2l09+x9ZCd7gj0ayg3+Ot17lppQrOIZkQHRTNj1wyGzxvOprhNFovzfNwRq9C4RRk8l+hW+vHEltLp7+Aq/Vj7qtLpz1A6eHmTZgvCPzuRm0fczKP3PMrHN3xMtj2bUYtGMXH9RDJzTDio4lLKm1UMhjJOYGUIqwWxpWS4Dq2EqOa6X0O5ItUrhIDsRB544IHzZTMHzuT1P1/n460f89ux33j56pdpXLmxhSo9k6KMuIqYTc1g8FCqtSqdEVdOFhxZp+MkGsodWV7B+NlTSU3VB0CwbzDPXfUc71z7DvFp8dy84Gam/jXVeB4WkaIYriUislpEPgAeRBsyH/fIMhgsJLoVxO+DjGT39nN8I2SlmvWtckq2dyB+9jT69u1L3759LzjXvWZ3Zg2aRY+aPXh7w9vcuehODiUeskip5+Gs4RqDjgKfA9wCTESvZa0Skd0iMlNE/gOYHZQGz6daK0DBya3u7eewIxK9Wd8ql+R4BxJIGpfKPRnuH87/rvkfr3Z9lf0J+xk+bzjTd053LlllBcepNS6l1NS8r0WkIdAG7VnYBugMDMmt7kqBBkOpk+ugEbsFanVyXz/H1kOl2hAc5b4+DJZh9w0mkCPkqBC8LxHDQUToV68f7au255nfn+GltS+x9PBSnu/yPNFB0aWs2HMollehUmqPUuo7pdRTSqm+SqnqQDTQF3jKpQoNhtImtDoERsCJze7t59gGqNHOvX0YrMM3kCBJJ8de+G/5qkFV+aDXB0y4cgKbTm1i6JyhzNs3z4y+LoHL3OGVUnFKqUVKqf+6qk2DwRJEdCQLd3oWJp2EhCPGcJVnfEMIIh27E4YL9OhrRJMRfD/ge+pXqs9TK5/i0eWPXjJZZUXGqalCEVlahDaVUqpnMfUYDGWD6Faw5n3IzgRvX9e3f2y9fszd8Gwod9j8ggkinb5DR1C7SrDT19UKrcVnvT/js22fMXnTZNafXM/zVz3PNTWvcaNaz8LZEZcN7YyRezRBpzOpgw4FVcfxujEeugFZRO4SESUig/OVPyUiu0TEnrUFdtcAACAASURBVP+coRwT3RLsWRC/1z3tH1sP4vX3epqh3CH+wfhIDtcNGMKdd95ZpGu9bF7c0/IepvefTpWAKoxdOpbnVj9nQkY5cMpwKaW6K6V6KKV6AG8DWUAnpVQ9pVRnpVQ9tINGluO8RyEitYHRQEF525eg1+5WlKoog7VEOjaFnt7lnvaPrYeqzcE30D3tGyzH5qtHWfGxRzh9+nSx2mgU3ohv+n3DXc3vYubumQyfN5wtp0ppc3wZpjhrXC8ATyul1uUtVEqtBZ4FXiyuGBGJEZF3HPvFUh0joDqXqFtTRL4XkQQRSRSRH0SkVjH6tAEfA/8EMvKfV0qtVUrtK2q7Bg8noiEgcMoNhstuh+PGMaO8I/4hADw97kFuvPHGYrfj6+XLuPbj+PiGj8myZ3HHj3fw3qb3yLJnuUqqx1Ecw9UQOHWJc3FAg+LLoQFwE3AW+O1SlUQkEFiKnrIcBdzu0PWriAQVsc9xwCql1PpiKTaUT3wDIbw2nNrp+rbP7If0BLO+Vc6x+TnWtVSOS9rrEN2BmQNn0rduX97///bOPL6q6lr835UJIjIKMsiQAEKkIKlQpAjik6cRGURBqg8I1Iqi1eprtVYLiq9U5VmqqK2iYjHUn6A4IAhSgVJF1IZaCDPUIjI8EARBZnKzfn+cc+NNcpPceUjW9/M5n+Tss88+655971lnr732WmufZezisbV20XIoims7cGslx24FvghZGvhAVZur6tXA61XUGw+0B4ap6tuqOh8YCrTzlU1EPhORA5VsbUTke8AIwhglGjWYZjnRGXF53exbdo9820bCkOoqLi2JjOICqJ9Rn0f6PcLv+v+OHUd2cP2C63lty2u1zm0+lCC7DwOviMh6YB6wD2iOowBygFGhCqOqJQFWHQp8oqqlM+equl1EPgKuAX7vll1UVSMichuOstvmJnluATwvIi1V9dkQPoJRk2jWGT5fDp5iSI1gPOp9GyAlDZpacNWaTKprKiSCistLXlYeuc1ymfTRJH7zyW/4YNcHTO4zmaaZTSN+rUQk6BGXqs4B8oDDwP3AH9y/3wB5qjo3ohL653uAv3g8G4AugTaiqs+qaktVzVLVLBznjFtMaRmAo1g8p+HQF5Ftd98GaNopOm72RsKQXqeu80/A7+PB0bxec5674jl+1etXfPJ/nzD8neH89cu/RuVaiUZAiktEfui7r6pLVfUSHFf4FkCmqvZV1WVRkNEfTXDmwcpzkIrJLcNCRCaKyC4cr8kXRWSXiFSIxSIit4jIahFZvX9/ZVOARlLRLMf5G+l5rn0bHI9Co0aTnuEorkFX9OO2226LyjVSJIVRF4xi7uC5nHvWufzsrz/joVUP1Xi3eQnENioiHhzHi3eAt4BlqhpVlxYRuRl4AchW1S/KHTsNTFPV+8uV/xa4T1XjmmesZ8+eunr16grlJSUl7Nq1i2PHjsVBKiNotAQO74K6jaBug8Rt06VevXq0bt2alJRo5Ic1guXMvs2kP3sx73WewlU33hn963nO8Ic1f+Cl9S/Run5rHun7CLnn5kb9upFERP6hqtV6LQX6gD8PJ4juNcDbwCkReQ9HiS1S1SMhSxoah3BGXeVpjP+RWEJw4MABRITOnTvbwyVZ2KtQ52xonBWZ9k4dha9PQ5MOEVVcJSUl7N69mwMHDnDuuRa0NxFIS68DwP59+9i5cydt2rSJ6vXSU9O5u8fd9Gvdj1+v/DVj3xvL+G7jubX7raSn1KwMVIEuQN7rzgddBTQDJuBEgX8O+EpElojIBBFpFUVZfdmAM89Vni7AxhjJUAERGSIizx8+fNjv8W+++YbmzZub0kom0utC8cnItVd84rt2I0hKSgrNmzensu+eEXskPROAx597hTFjxsTsuj2a92DekHkMaT+EGUUzyF+Uz/bD22N2/VgQinPGt6r6qqregKPErsVxkZ8E7BSRT0XkVxGWszzvAL1FpL23wF2ofIl7LC6o6gJVvaVhw4Z+j3s8HtLTa9abT40nrS4Un6LSpErBcuaEE+opCm/A6enpFBdbJt2EIdV1vomDq/rZGWczpe8Ufn/Z79l5dCcjF4xk7ua5NcZtPqxXf1U9o6qLVXWCqp4H9AVWAPmhtikiI0RkBOANKzDQLfONMPkCznqx+SJyjYgMBeYDO4EZoV47FkgleXmMBCWtjjMvFanU6mdOQHqmE4E+wth3K8FIc0yF0fIqDIQr2l3BW0PfokeLHkz5dAq3L7udAydCCz+VSETUZqWqH6vqfaoasEu6H153twnu/h/d/Yd9rnMMuBzYCswGXsEZ9V2uqlHOt27UKlLdh09xhWhgwaPqmB1dE5JRw/F+d+I8yml2VjOeHfAsD1z8AIV7C7l2/rUs+zJWDuDRIWzFJSLfisjfROQJERktIl0kjFc/VZVKtsvK1ftSVYeragNVra+qw8p7H8aa6ua4jCTEu9bKEwHF5TntvH2nRXZ+y0hQUtPwkALEb8TlRUS4MedGXhvyGi3rteTuv96d1NHmIzHiOh94HGcB8o+ApcC3IvKRiDwVgfaThurmuBKdlStX0qdPHxo2bEiTJk245JJLKCwsZNasWYgIP//5z8vUf/vttxGRoFM2JBXeeYpIjLi8bZjiqjWcIY0R/b/HL37xi3iLAkD7hu155epX+EnXn/DG1jf40cIfseHrDfEWK2jCVlyux+FCVX1YVYeoaivgOqAO8NOwJTRiwpEjRxg8eDB33nknBw8eZPfu3Tz00EPUqeOYOzp06MDcuXPLTP4XFBTQqVOneIkcGyTFUV4RUVyud6J37sOo8ZyRDPpc0IIhQ4bEW5RSvG7zM/NmcqL4BKPfHc3MdTPxRCE0VbSI2ByXiHQXkd+IyCacEdgCwKKIJglbt24F4MYbbyQ1NZXMzEyuvPJKLrzQSXTYokULunXrxpIlSwA4ePAgq1atYujQoXGTOWak1YmMqbD4lOtRGNf18UYMKZZ0du89wJYtUcrrFgbeaPOXt72cJz97kvHvj2fvsb3xFisgwv4FicijwHAcU+EbwFBV3RZuu8mIiAwBhnTsGFhml4cXbGDjnuiu3e7SqgEPDak+vFCnTp1ITU1l7Nix3HDDDfTu3ZvGjctGz8rPz6egoIBBgwYxZ84crrnmmtIRWY0mtQ6cPu5MsofjuVd80lGC5v1XayiWDJ54o5BGhbeyYsWKeItTgYZ1GvK7/r9j/ufzefTTR7nunet46IcPkZeVF2/RqiQSI677gK9xPP8KqTxXV40nmee4GjRowMqVKxERxo8fT7NmzRg6dCj79u0rrXPttdeyYsUKDh8+TEFBAfn5Ia96SC7S6jg5lcJ1iS8+ZfNbtYxiSUcSwDmjKkSEYR2H8fqQ18lukM09f7uHiSsncuxM4oami4TN4mKcNVc9gBuB74nITuAfwGeqOjUC16iRBDISiiUXXHABs2bNAmDz5s2MHj2au+++m7w85+0rMzOTQYMGMWXKFA4cOMAll1zC4sWL4yhxjPA6aHhOQ2qIC4dLPFByxua3ahmelPS4u8MHStsGbZk1cBYz1s7ghXUv8I99/+CxSx+je7PEm/GJhHNGoao+p6rj3fxX9YEbcDIUdwi3fSM+5OTkMG7cONavL5s9Jj8/n2nTpsU0hE3cSfNRXKFiHoW1Eo9kJPyIy5f0lHTu+P4d/CnvT5RoCWMXj+W5tc8lnONGJJ0zCkRkG/AXoL6qzlDVWyLVvhFdNm/ezLRp09i1axcAO3fu5NVXX6V3795l6vXv35/333+fO++MfrTrhCE1AorL69xhI65ahSclAyE5Rly+XNT8IuYNnUdeVh5/WPMHblpyE7uP7o63WKVEMnLGW8B/AL8FnhCRn0Sw7aQgmRcg169fn08//ZSLL76YevXq0bt3b7p27cq0adPK1BMRBgwYQJMm/oLz11BS0hxvQE8YmXy8I65USx5ZmyiWDCZc2pKJEyfGW5SgqZ9Rn6mXTuXRfo+y9dBWRrwzgnf//W68xQICzMcVdKMiDYBFwNMxyoicUFSWj2vTpk1ccMEFcZDICJuvNjnehee0r76uP775Ek4ehhbdIitXOew7llhsfDyPlBP7yXnws3iLEha7vt3FAysf4J9f/ZNB7Qfx64t/Tf2M+hG/TqD5uCJpKswSkX4iMgq4HSd24KxItW8YcSU1I/w5Lhtt1TqKU9LZtucIa9asibcoYdG6fmteynuJn+b+lPe2v8eId0bw2b74KeNImgrX4gTEHYPjlLENR4EZRvITruLynP4u6KpRayiWdB55bw933313vEUJm7SUNCZ0n8DLA18mRVL48ZIf88w/n+FMSRgm9FBlCfYEEflvnKzHZZaCq2ryLV4yjEBJy3DXcnkgJTW4c1UdxZXZuPq6Ro2ihFRq2nLz7s26M2/oPB799FFmFM3g4z0f81i/x2jTILoZnn0JZcQ1DfivSAtSE0hm5wyjGsLxLPSeY6bCWofW0Cgp9dLrMaXvFB7v/zjbj2xnxIIRvP2vt2OWqDJUU+F/icg+ETktIl+KyFMi0jWikiUhyRw5w6iGSCguc4U3ahhXZV3FG0PeoMs5XZj00STu+ds9HDkd3TB2ELriygJW4gTTXQpcDxSKyISqTjKMpCUcxWWu8EYNpuXZLXnxyhe566K72HpoKxID42ioIZ8eVtUp3h0RSQF+BTwjIl+q6qKISGcYiYI3onsoa7nMVFiruW/AOTT7ySPxFiOqpKakcnO3m8nvkk9GDL7noSiuM8DHvgWqWgI8IiLn4QTdNcVl1CxEICXdiTcYLJ4zzrk1dL7DqJoftM2kVZ8+8RYjJsRCaUFopsIvcQLq+mN+FccMI7lJTQ99xGWjrVpL4ZcnWLVqVbzFqFGEorjeBCaKyOV+jnWAJIooaZRh5cqV9OnTh4YNG9KkSRMuueQSCgsLAcjKyiIjI4MDBw6UOSc3NxcR4YsvvoiDxLHhmWeeoWfPntRpfSHjfnpvtfUr3CvPGUhND/he5eXl8eCDD1Yonz9/Pi1atCiThdpIfKYu+5oHHngg3mLUKEJRXA/jmAr/IiKLReRnInK9iEwCpuJEha+VJLM7/JEjRxg8eDB33nknBw8eZPfu3Tz00ENlEkVmZ2fz6quvlu6vW7eOEydOBHyNyZMnM3ny5EiKHRNatWrFxIkTuWnUSNDA3stK75W7hmvd5s8Dvlfjxo1j9uzZFVyLZ8+ezahRo0hLswzKRu0maMWlqsdVNQ/4GXAO8CQwF0ehbQPuiKiESUQyu8Nv3boVgBtvvJHU1FQyMzO58sorufDCC0vrjBkzhoKCgtL9l19+OWrJJF988UXy8vK47bbbaNy4MZ06dWLjxo1Mnz6dtm3b0rRpU958882oXLs81113HcOGDeOcc85xFFEAyqv0XpV4AOXlOW9WuFd79uxh+PDhNGvWjOzsbJ566ikAhg0bxsGDB/nwww9L6x46dIiFCxfWnuSdhlEFIYd8UtU/qmovoAVOMsnOqtpDVXdFTDojZnTq1InU1FTGjh3L4sWLOXToUIU6vXv35siRI2zatAmPx8PcuXMZPXp0VOQpKipi9erVjBgxggMHDtCtWzcGDhwIwOeff86kSZOYMmVKNa34Z/DgwTRq1MjvNnjw4MpP9EbMCGCeq/RebShy7tWb88vcq5KSEoYMGUL37t3ZvXs3y5Yt48knn2TJkiVkZmYycuTIMi8Jr732Gjk5OXTvnnhJ/Qwj1oRtc1DVr4CvIiBL7WPxr2Dvuuheo0U3GPhYtdUaNGjAypUrmTp1KuPHj2fv3r1cffXVvPDCCzRv3ry0nnck0b9/f3JycjjvvPOiIvbatWu5//77GTBgAABdunTh1KlT3HXXXQB07dq1zFxPXl4ep04566V2795No0aN+NOf/sSKFSu4446yRoCFCxeGJpS473mewDIZjxkzhoLZf6Z/bntyOncuc68KCwvZv39/6VxW+/btGT9+PHPmzCEvL4+xY8cyaNAgnn76aTIzMykoKGDs2LGhyW0YNYyAFJeIBDNvpao6IER5jDhywQUXMGvWLMBJLDl69GjuvvvuMvNaY8aM4dJLL2X79u0Bma0GDx7MypUrATh58iQATz75JAB9+/atVIkUFRUxY8aM0v2NGzeWGQ1t3LiRnJyc0v0lS5YAsG/fPq644gpeeuklli9fzsiRIwP56IHhVVwBusSPGTOGS/v1ZfvWLuTnl80YvWPHDvbs2UOjRo1KyzweD/369QOce9OsWTPmz59Pr169KCwsjJlp1Igsk69qxrkTnoy3GDWKQEdcKVAmjWdnHBPhF8A+oDlONI3/A7ZgBEYAI6F4kZOTw7hx48ooD4B27dqRnZ3NokWLmDlzZrXt+Comr2NGdQ4aO3bs4PTp03Tq1Km0bM2aNdx773cefUVFReTm5pY57/jx41x77bU8/vjjdOvWjTVr1tCiRYsK7Q8cOLDM/JEv/fr1Y/Hixf4Fk8BNheDeq3ZtWLT8I2a+8kaZY23atCE7O5tt27ZVen5+fj4FBQVs2bKFK6+8sszI10geurasQ6ty31UjPAKa41LVy1T1P1T1P4DpOIuQe6tqe1X9oaq2B37olk+PnrhGtNi8eTPTpk1j1y5ninLnzp28+uqr9O7du0LdmTNnsnz5curVqxcVWdauXUu3bt1ISXG+nkeOHGHHjh1lHEXWrl1bZr6npKSE0aNHk5+fT15eHuCMePyxePFijh496nfzp7SKi4s5efIknpISPJ4STh77NmCX9JnTH2P56zOpd/bZZcp79epFgwYNmDp1KidOnMDj8bB+/frS5QfgKK6lS5fywgsvmJkwifng8+MsXbo03mLUKEJxzvgNMElV/+5bqKqfApOB0GbMjbhSv359Pv30Uy6++GLq1atH79696dq1K9OmTatQt0OHDvTsWW2S0pBZu3ZtmdFUUVERHTt25KyzzgIcJbV+/foydX75y1/Svn17JkyIfLjMKVOmkJmZyWNTp/LnN98ls0XHgB1DOrRtSc8eFd+2U1NTWbBgAWvWrCE7O5umTZty880347uUIisriz59+nDs2DGGDh0asc9jxJanPjgYsiOR4R8JNgy9iJwErlXVCq+mInI18IaqZkZIvqSkZ8+eunr16grlllY9OsyYMYMlS5Ywb9680lFa1Ni/xfEuPKdjYPW/2uQ4cjRpH125XOw7llgUPnE9d01/l4ysXqxYsSLe4iQ8IvIPVa32rTiUX/l24NZKjt2KM+9lGDFh79693HHHHezevZvLL7+cyy67jBtuuCF6F0xJg5IgIleUFDtxCg3DiBihuMM/DLwiIuuBeXznnDECyAFGRU685EJEhgBDOnYM8G3cCJsWLVpw5kwMU4enpMGZAKOFqLqKyyJdGEYkCSVyxhwgDzgM3A/8wf37DZCnqnMjKmESkcyRM4wA8Y64AjGxe0dmqaa4DCOShPSLUtWlwFI3D1dT4ICb2sQwajapaYAb9snrHl8ZXrd5MxXWah4bci7N75hRfUUjYMJ6FXSVlUXNMGoPXrNfSfF3IaAqwzviMlNhraZj0wxade4cbzFqFFF2wTKMGoav4qoOb4QNMxXWav6y5SgLFiyItxg1ClNchhEMwSguj3fEZabC2szzq77xux7SCB1TXIYRDF7FFUjYp5JiJ75hdSZFwzCCwhSXYQRDsKZCm98yjIhjisswgiEl1RlFBWQqNMVlGNHAFJdhBEug0TMsaoZhRAVTXEbC8Mwzz9CzZ0/q1KnDuHHjqqyblZVFRkYGBw4cKFOem5uLiPDFF19Ue728vLzSRI6+zJ8/nxYtWlQeAT4YxWUehbWe6dc1Z/bs2fEWo0ZhistIGFq1asXEiRO56aabAqqfnZ1dJsnlunXrOHEiwHBMwLhx45g9ezblA03Pnj2bUaNGkZZWidJJSfvOY7AyLNyT4XJew3TatGkTbzFqFKa4IoiIDBGR531TUyQTL774Inl5edx22200btyYTp06sXHjRqZPn07btm1p2rRpVLPwXnfddQwbNoxzzjknoPpjxoyhoKCgdP/ll1+ukJV5z549DB8+nGbNmpGdnc1TTz1VemzYsGEcPHiwTFLJQ4cOsXDhwqqzO6ekgnqqFs573BRXrWf++m+ZO7fWRsKLCqa4IkiyxyosKipi9erVjBgxggMHDtCtWzcGDhwIwOeff86kSZMCzis0ePBgGjVq5HcbPHhwROTt3bs3R44cYdOmTXg8HubOncvo0aNLj5eUlDBkyBC6d+/O7t27WbZsGU8++SRLliwBIDMzk5EjR5ZRfq+99ho5OTllklRWICUNSqpRXN7j5gpf65ldeJhnn3023mLUKOx1MI5M/ftUNh/cHNVr5DTJ4b5e9wVUd+3atdx///0MGDAAgC5dunDq1CnuuusuALp27Vo677No0SL+93//F4/Hg8fj4cMPPyQ19buH9MKFCyP8SfzjHXX179+fnJwczjvvvNJjhYWF7N+/v3Qeq3379owfP545c+aUZkkeO3YsgwYN4umnnyYzM5OCgoLqsw2LO+JSBRH/dbyKS+wnZhiRxn5VRilFRUXMmPFdMNCNGzeWGR1t3LiRnJwcvv32Wx588EE++ugj6tSpw6FDh8oorVgyZswYLr30UrZv317BvLdjxw727NlDo0aNSss8Hg/9+vUr3e/bty/NmjVj/vz59OrVi8LCwurNod5RlHoqV0ylcQptxGUYkcYUVxwJdCQUC3bs2MHp06fp1KlTadmaNWu49957S/eLiorIzc1l7969HD16lJ07d9KxY0caN25cob2BAweWmTvypV+/fixeXCGBdki0a9eO7OxsFi1axMyZM8sca9OmDdnZ2Wzbtq3KNvLz8ykoKGDLli1ceeWVNG/evOqLepVRiafyOSw1U6FhRAub4zIAx0zYrVs3UlKcr8SRI0fYsWMHF154YZk63bt35/zzz+dnP/sZI0aMoEePHqxbt65Ce4sXL+bo0aN+t8qUVnFxMSdPniw1P548ebJyl3QfZs6cyfLly6lXr16Z8l69etGgQQOmTp3KiRMn8Hg8rF+/nsLCwjL18vPzWbp0KS+88EL1ZkL4Lp1JVfNcFhneMKKGKS4DcJRSbm5u6X5RUREdO3bkrLPOAhxHh/Xr15fWuf3221mzZg033XQT06dPj4gMU6ZMITMzk8cee4w///nPZGZmBuQM0qFDB3r27FmhPDU1lQULFrBmzRqys7Np2rQpN998M+W9PrOysujTpw/Hjh1j6NCh1QvqayqsDHPOMFyeH9mSefPmxVuMGoW9DhoATJo0qcx+37592bz5O8eRlJQUjh8/DsDWrVtLTYrffPMNnSOUa2jy5MlMnjw5oLqVLTBOS0srsy6rVatWZdZ6VcaKFSsCui5Q1lRYGSUeJzSU2LthbadJvVSaNm0abzFqFKa4jKB58MEH2b59O5mZmXTu3Jmnn3463iLFFgkg0G5JcfUZko1awdx/HqHxrFnVRoMxAscUlxE0c+bMibcI8SVQU6HNbxnA62uOkPGNKa5IYnYMwwgWr/mvKlOhFtv8lmFECVNchhEsIo4ZsLo5LlNchhEVTHEZRihUF6/QAuwaRtQwxWUYoVBVvEJVG3EZRhQxxRVjyqfQMJKUqkyFWgJozOMU2ncrMZk9qhWLFi2Ktxg1ClNcLiLyYxFRERnmU1ZXRN4WkU0iskZElohI+1CvUbduXb7++mt7wNQEUlIdBwx/xGHxsary9ddfU7du3Zhd0wiMzIyU0oX8RmQwIzwgIu2A8cAnfg4/q6pL3Hp3AC8Cl4dyndatW7Nr1y72798fsqxGgnD8IBSfgK/9RIf3nIZvv4J6JZAeu76uW7curVu3jtn1jMCY9fdvaPTHP3L77bfHW5QaQ0IpLhFpDdwH9AS6A5lAtqp+4aduG+AJ4ApAgKXA3ar6ZZDXTAFmAncC03yPqepJYIlP0SfAPcG070t6ejrZ2dmhnm4kEu8/CJ88B5O+qnjsy09g3kgY/QZ07BF72YyEYuGGo2Qce80UVwRJNFNhR2AkcAjwH1ocEJGzgOVADjAWGAOcD/xVROpVdl4l/Bz4SFX/EUDdO4H5QbZv1ETqNgTPKThzsuKx00edvxn1YyuTYdQSEmrEBXygqs0BRORm4MpK6o0H2gOdVfVfbv0iYBtwK/B7t+wzoG0lbXwfaACMAPpVUqcUEbkf6AQMCPTDGDUYr1I6fRTSy80rnXIVV52zYyuTYdQSEkpxqWpJgFWHAp94lZZ77nYR+Qi4BldxqepFVTUiIrcB7YBt4mSybQE8LyItVfVZn3r3AMOB/1TV40F8JKOm4lVKp76FeuUCqJaOuExxGUY0SDRTYaB8D1jvp3wD0CXQRlT1WVVtqapZqpqFM4d1Szml9XPgRuAKVf0mPLGNGoNXKXmVlC+njzl/65ip0AAwL+JIk1AjriBogjMPVp6DQMV0vCHiOotMA/6NM38GUKyqFZI/icgtwC1AQ+CoiGwpV6UhcDiAsqbAgfClDxp/ssSqnUDPqa5eVccrOxZIv1TeJw9f6LfYOXZO5ccCJ179kuh9Asn0W9nxN9znR3jtxK9fAu2rcPukXUC1VDUhN+BmnFeVLD/HTgOP+in/LY5iiafczwdaXknZ6kSSOxbtBHpOdfWqOh5Ov8SrT+LZL4neJ/Hsl9r4Wwm0r2LVJ8lqKjyEM+oqT2P8j8RiyYIgyiurGw8iJUso7QR6TnX1qjpu/RKd+tYnsW0nXv0STF9FHXG1ZMLhehW+gJ91XCKyHMhQ1b7lylfgfKb+sZIz0ojIavVjijTih/VJYmL9knjEqk+SdcT1DtDbN/ySiGQBl7jHkpnn4y2AUQHrk8TE+iXxiEmfJNyIS0RGuP8OACYAtwP7gf2q+je3Tj1gLXACmIgzF/YboD5woar6cfUyDMMwagKJqLgqE+hvqnqZT722lA35tAwn5NMX0ZbR+QWlnQAACs9JREFUMAzDiB8Jp7gMwzAMoyqSdY7LcBGRZSKy1k278qGI5MZbptpMpFPhGJFDRB4QkS0iUuKbvsiIHSLSQURWishWEfmniITkyGGKK/m5TlW7q2ouTqirWXGWx3BS4Vzg9skCnFQ4RvxZBlwNfBBvQWoxzwGzVLUT8EvgFfGzMrs6THFFEBFpLSJPi8jHInLcTUyZVUndNiIyT0QOi8gREXnTnbcLClX1XbXeIETRayyx7hNVPalu/jaXT3ACQhs+xOm38qmqfh6u7LWJSPaTiDQDegMvA6jq++6hoHP/mOKKLPFIy4KIvCIiu3A8K0eHIHdNJi594oOlwvFPvPvFCIxI9lNbYI+qnvE5dQeVZ/ConFiE56gtG5Di839VIavuAjxAR5+ybKAY+LlP2Wc4cb/8bW38tHsz8G6870MibfHsE+B+4GPgrHjfh0Tb4twvK4Bh8b4HybBFsp9wRlZbyp33Ps50R1By2YgrgmiYaVkAb1oWb9lFqtq0km2nn3ZnAleISESiu9YE4tUnPqlwBqqlwqlAAvxWjACIcD99CbQSkXSf89q55UFhiis+RCQti4g0FpGWPkXDga9wouQbwRGRPgFLhRNhItYvRlSptp9UdT/wd2AcgIh41+AGkn2+DMma1iTZiVRalsbAXBGpC5TgKK3B6o7BjaCISJ8EkwrHCIiIpTASkYk40XiaAV1F5Bmgp6ruDVtKI9B+mgC8LCL3AseBUaE8r0xxxQ9/nRWUW6iq/hv4QWTEMYhMn+wK9hyjWsLuFwBVnQJMCV8coxKq7SdV3Qb0CfdCZiqMD4mclqW2Yn2SmFi/JAcx7SdTXPFhA45NuDxdgI0xlsVwsD5JTKxfkoOY9pMprvhQk9OyJCvWJ4mJ9UtyENN+siC7EcbSsiQe1ieJifVLcpCI/WSKK8JYWpbEw/okMbF+SQ4SsZ9McRmGYRhJhc1xGYZhGEmFKS7DMAwjqTDFZRiGYSQVprgMwzCMpMIUl2EYhpFUmOIyDMMwkgpTXIZhGEZSYYrLMAzDSCpMcRmGYRhJhSkuI+6IyDgRUZ/ttIh8LiKPuEky4ynbZN+QN+X3jeoRkfHl+td3OyYiqWG0PcunrRU+5ZPdsrByDorIfBE5KCJ1Kjle3/0Ms9z9iT7y7Arn2kblWCJJI5G4HtiFE5jzWuB+9/874yjTi8B7cbx+TSAXOAb8p59jJ1TVE2b7e3G+L0fCbMcfLwNDgcHAG36OjwDOcusB/AlYCkwCukdBHgNTXEZisUZV/+X+/76InA/8RETuUtWSeAjkZjRO6jdnEamjqqfiKEIusEFVP4lS+6ei2PZC4GsgH/+KKx/4ElgBoKq7gd0isj9K8hiYqdBIbD4DMoGmvoUi0l1E3hGRQyJyQkQ+EpF+fuq8JSJfu3W2iMj95epcJSIfu8cPi8jbItK5XJ2QTYMByuk1aZ0vIu+KyFER2SEiD4pIShjtdRWRJSJyFHjNPXajiGwWkZMisk5EhorIinImthHu+RVGC27dj4O8BwJcCBQFc140cPv7qIg847231d1TVT0NzAEGikj572FboD8wWy1aeUwxxWUkMlnAYZw3XgBE5CJgFU6a8PHAcPf4UhHp4dbpBXwMdAD+GxgE/B5o7dPOVcC7wFHgR8BtQFdgpYicF67ggchZjreA5cAw4G3gYWBsGO3NB/6GY+Z6QkSuAF4BNrvn/g54EuhU7ry3gT3AreU+T2ech/SMgG7Ad5wPnA1sEJG0clvIc1vBIiL5OAkNp6rqHapaEsQ9fRlIx/me+DIaJ31HQbTlN8qhqrbZFtcNGIeTeK4zjvm6MXATUAzcUa7uMmATkOFTluqWve3ufwDsBM6q4pqrgW1Amk9ZNnAG+L1P2WTnZ+J/v4r2q5XTtz3gx+XOXwf8JYz27irX3ipgPW4qI7fsIrfuinJ1J+O8MNTzKfs9cAjIDLJvR7rX8Lf9MwLfnVnAF37KvfchDfil2683h9JHbvkG4NNyZZuAVVXItStev6mavtmIy0gkNuM8YA4CM4EZqvqM96CIZOK89b8OlHjf3HHeepcCl4rIWTjpwl9R1eP+LiJOttaLgLmqWuwtV9XtwEfuNUImEDn9nPZuuf31QNsw2nvLR55UoCfwhrpPVQBV/QzY7ufc53EcDm50z6+LM/orUNUTVX/6CnwfR4FcAvyg3ObNrIuI9BKR4yKyxjVnrhKRjkFeyx9P4IxeR6jqiz7XC/aeFgC9RKSTV14gBxttxQVTXEYicS3OA+1qnIfH7a6Jx0sTnDfiSTgKzne7A2ek1hjne12VQ0VjnAfU//k5tte9TjhUK2f5+SscZe3LKcC7FCCU9nw/W1McU9dXfmTdV75AVffgmBonuEXXuzIEayYExzHj36q6SlVXl9s+96nXA1iqqrmqmgN8gfN5w+VGnNHS0nLlwd7TPwMlOM4YuH9PAXMjIKMRJOZVaCQS69X1KhSR5TgT+o+LyBuqegz4Bufh8Qcqf9M95Napap7qEM4ooIWfYy3wmVMLkWrl1OC8JENpz9dZ4ADOA/lcP6c2x/GKK88fgWXuXM+twIequjEImb3kAoF4/PXAUTBe/gW0CuF65RkA/AVYLCJXq+pRtzyoe6qqu0VkKTBaRP4HZ77rHVU9FAEZjSAxxWUkJKp6SkTuxXnzvx14XFWPiciHOOtjPqvs4S8iK3EfMP5MW247/wCuF5HJ6q4jEpF2QB/g6TBlD0jOWLWnqh4RWQ0Mdz+vArhKKRs/iktVl4vIJpy5rUuAUcHKLSLNcV4EAvEo7AFMc8/LAvrhOEyEywbgMhzHl/dEZKCqfhviPX0Zx8HlUZxRrJkJ44QpLiNhUdV3RKQQuEdEnnGV0M9xnC+WiMhMHJNYU5w5q1RV/RVwD45H3cciMg3HbNgeyFVV72LmSTjzSgtF5I84nm8P4zglTIuA+IHIGcv2HsIZebwlIs+7507GMY1W9tB+DpiOM2Irs4bJVS7bgYdVdXIl53/f/VssIr39HC9S1ePiRKX4HjBRRCbh9NVI/W5NX1io6iYRuQz4K47yukpVvyX4e/oWziLn/8Yxu9rC9Dhhc1xGojMRx8Q1AUodCn6AY857CudhPB3ohvMQQlULcUYJO3FGT4uAe/GZ91LV93Dc5BvhrHN6DsdLrK87xxMWgcgZy/ZU9X2cUdMFOA/g+4Bf4Ciuw5Wc9rr7d5ZWXMBcz/27t4rL5rp//wdneYLvtgrIcI9fCJwALlDVzm79x6r7TMGgqltwnDHaAX8RkQbB3lP3xel1nPnR/+fr2GPEFvFxMjIMoxYhIq1x5pJ+q6q/8XN8PI5DRqfyox8RuQX4LdCuMu/NIOSYAAxX1Svc/bNxRjQ/VNW11Zw7C8cU2BFnmUK44aPCwl1wnYrjFTtAVVtXc4oRAjbiMoxagIhkisizIjJcRPqLyI+B94HjOPEYfet2EZEhOKbTtysx2fUHnghXabn0wMeBw3WgWISzGDsQ2uE4nyyLgCzh8mscWfKrq2iEjo24DKMWICIZOK7bvYFzcILefgg8oKrry9VdgeOksgr4r0iYTqOFO9fmDcX0rWsSjBsi0pLvPFpPq2rcQ13VRExxGYZhGEmFmQoNwzCMpMIUl2EYhpFUmOIyDMMwkgpTXIZhGEZSYYrLMAzDSCpMcRmGYRhJhSkuwzAMI6kwxWUYhmEkFf8fuJg3TjPtnw8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Specify a single neutrino energy\n", "E_nu = 3 #MeV\n", "\n", "#Specify a nucleus - Germanium\n", "A_Ge = 73\n", "Z_Ge = 32\n", "\n", "#Calculate differential xsec\n", "E_R = np.logspace(-3,0,500)\n", "dsigdE = CEvNS.xsec_CEvNS(E_R, E_nu, A_Ge, Z_Ge)\n", "\n", "#Specify New Physics parameters - Z'\n", "m_V = 1.0 #Mediator mass, 1 MeV\n", "gsq_V = (2e-6)**2 #(g^V_q)*(g^L_nu)\n", "\n", "#Specify New Physics parameters - S\n", "m_S = 1e0 #Mediator mass, 1 MeV\n", "gsq_S = (1.5e-6)**2 #(g^S_q)*(g^S_nu)\n", "\n", "#Calculate differential xsec including New Physics\n", "dsigdE_vector = CEvNS.xsec_CEvNS(E_R, E_nu, A_Ge, Z_Ge, \\\n", " gsq_V, m_V)\n", "\n", "dsigdE_scalar = CEvNS.xsec_scalar(E_R, E_nu, A_Ge, Z_Ge, \\\n", " gsq_S, m_S)\n", "\n", "#Maximum recoil energy\n", "Emax = CEvNS.ERmax(E_nu, A_Ge)\n", "print(\"Max. Ge recoil energy (E_nu = 3 MeV): \", Emax, \" keV\")\n", "\n", "#Do some plotting \n", "pl.figure()\n", "pl.loglog(E_R, dsigdE, label=\"SM\")\n", "pl.loglog(E_R, dsigdE_vector, label=r\"SM + $m_{Z'} = 1\\,\\, \\mathrm{MeV}$\")\n", "pl.loglog(E_R, dsigdE_scalar, label=r\"$m_{S} = 1\\,\\, \\mathrm{MeV}$\")\n", "\n", "pl.axvline(Emax, color='k', linestyle='--')\n", "\n", "pl.ylim(1e-42, 1e-39)\n", "\n", "pl.legend(fontsize=12.0, loc=\"best\")\n", "\n", "pl.title(r\"CEvNS on a Germanium target $(E_\\nu = 3 \\,\\,\\mathrm{MeV})$\", fontsize=12)\n", "pl.xlabel(r\"Recoil energy, $E_R$ [keV]\")\n", "pl.ylabel(r\"$\\mathrm{d}\\sigma_{\\nu N}/\\mathrm{d}E_R$ [cm$^2$/keV]\")\n", "pl.show()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Recoil Rates\n", "\n", "Calculate the recoil rate in a Germanium target, assuming a neutrino flux from the Chooz reactor.\n", "\n", "First, let's load the neutrino flux (need to choose `source=\"CHOOZ\"` or `source=\"SNS\"`):" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAAEoCAYAAAAXPoWWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3XmcHGW59//PdzKZLJOFTFayQ8IWdogaRWVRUNAACsLxIJsQRI/+5PGgqIACgjznHAQVRTaV5SAiKJv6yE7Y1QCyhC0sCQkQCCQkkIUsc/3+qJrQ6XTPdPd0T/dMf9+vV79muuquqqu7Zvrqe6m7FBGYmZlVW0O1AzAzMwMnJDMzqxFOSGZmVhOckMzMrCY4IZmZWU1wQjIzs5rghGTWAUkfk/Rshfb9OUnzJb0raWdJcyV9sgLH2UrSo5LekfT/lXv/nSXpKEn3VTsOqy4nJKtp6Qf065KaM5YdK+nuMu0/JE1ur0xE3BsRW5XjeDmcA3w9IgZExKMVOgbAd4C7I2JgRPy8gsfJS9KnJN2TJsVFkmZK2r8asVhtckKy7qAR+GY1DiypscKHmADMrvAxOjyOpF6VPLikg4FrgSuAscBI4AfA9Eoe17oXJyTrDv4HOFHSJrlWStpa0m2SFkt6VtIhGevulnRsxvP1TUOS7kkXP5Y2mR0qaQ9JCySdJGkh8Nu2ZRn7mCvpREmPS1oq6RpJfTPWz5D0fBrPTZJG54i5j6R3gV7p8V/IUeYySWdmPF8fh6RJ6f53SZ+PlvSmpD1y7OdOYE/gF+nr3DLd968k/VXScmBPSYMlXZHWXuZJOkVSQ8b7dr+k8yS9LelFSR9Jl8+X9IakI/OcHwHnAj+KiEsjYmlEtEbEzIiYkVX2HElLJL0kad+M5aPT93Jx+t7OyFjXR9JPJb2aPn4qqU+67ub0Nbc9WiUdlStOqz4nJOsOZgF3Aydmr0ib8m4DfgeMAL4IXCBp2452GhEfT3/dMW0yuyZ9PgpoIalVHJdn80OATwObATsAR6Xx7AWcna7fFJgH/D7Hsd+LiAEZx5/UUbxZ278AnARcJak/8Fvgsoi4O0fZvYB7eb9p8Ll01b8DZwEDgfuA84HBwObA7sARwNEZu/oQ8DgwlOT9/j3wAWAy8CWShDeAjW0FjAOu6+BlfQh4FhgG/Dfw6zSZAVwNLABGAwcDP5b0iXTdycA0YCdgR+CDwCnpa5+evuYB6XYLgTs6iMOqxAnJuosfAN+QNDxr+WeBuRHx24hYGxGPAH8k+fApVSvwwzRprMxT5ucR8WpELAZuJvkwBDgM+E1EPBIR7wHfAz4saWIn4skpIi4B5gB/J0l+Jxe5ixsj4v6IaAXWAIcC34uIdyJiLvAT4PCM8i+l7/M64BqSJHNG+j7dCqwmSU7ZhqY/X+sgnnkRcUm6/8vT1zRS0jjgo8BJEbEqIv4FXJoR22FpHG9ExCLg9Ky4kbQlSXPhoRExv6M3xqrDCcm6hYh4Evgz8N2sVROAD6XNSG9LepvkA2pUJw63KCJWdVBmYcbvK4C2msFoklpRW9zvAm8BYzoRT3suAbYDzk8TYDEyP5iHAU1kxJ7+nhn36xm/rwSIiOxluWpIb6U/N+0gnvXvaUSsSH8dQPKeLo6Id/LEtsF7nv6+vplU0mDgRuDUiLi3gxisipyQrDv5ITCDDT8k5wMzI2KTjMeAiPhqun450D+jfCGJqjNT4L9KkiSB9U2KQ4FXSthXu7GnzWM/BX4NnCappcj9Z77ON0lqSRMylo2ntLizPUtyng4qcftXgRZJAzOWZca2wXuernsVIO0D+x1wV0RcVOLxrYs4IVm3ERHPkzQVZV5H82dgS0mHS+qdPj4gaZt0/b+Az0vqr2R49zFZu32dpM+kXH4HHC1pp7Rj/cfA39MmsGL9C9hPUoukUcAJWet/BjwcEccCfwEuLDXotJnsD8BZkgZKmgB8C/jfUveZse9I93WqpKMlDZLUIOmjki4uYPv5wAPA2ZL6StqB5DxelRa5GjhF0nBJw0iad9viPgtopkqjNK04TkjW3ZxB8gEDQNqMsw/wbyTfihcC/wX0SYucR9K38TpJv8RVbOg04PK0ue8QOiki7gBOJenHeg2YlMZWiiuBx4C5wK0kyRgASQeQDKo4Pl30LWAXSYeVeCyAb5DUyl4kGeTwO+A3ndjfehFxHUkf1ZdJztPrwJkkTWmF+CIwMd32epI+vtvSdWeSDHx5HHgCeCRd1rbdNGBJxki7zrxHVkHyDfrMzKwWuIZkZmY1wQnJzMxqghOSmZnVBCckMzOrCU5IZmZWEyo9k3GPMmzYsJg4cWK1wzAz61YefvjhNyMie9qvjTghFUDSdGD65MmTmTVrVrXDMTPrViTN67iUm+wKEhE3R8RxgwcPrnYoZmY9lhOSmZnVBCckMzOrCU5IZmZWE5yQzMysJjghFUDSdEkXL126tNqhmJn1WB72XYCIuBm4eerUqTNK2X7VmnW8sez9m3lGnvu/5Zt4vb352DNna89XbsP9Rp7luY+XWSbybJuvTHvHyFem2PemPVKe5Wij9W2/CyFBg5KfAqT3lzVkrEuei4aG5PdeEg0NSZnGhgYaGtKfSvZhZu1zQuoCT7+2jM9d8EC1w7AqamwQvRpE39696N/Ui+Y+jXxsi2Ec9ZGJTBja3PEOzOqAE1IXGN/Sn598YccNluX99t7Bt/pittmwzPuFtMHyjo+3QU0i77a5y2dv06mYiqhk5K39pcuD92uXscG6IAJaI6mtJb/H+vXrWpM6XGsEra1Ba7q+NYJ1bc9bg3Xp87Xr2n5vZe26YNWadaxYvY7Fy1dz5YPzuOyBuey9zUi+tudkdhq3SeEv0KwH6rEJSdL3gSOBLYDPR8QNGesmkdw9dATJHTJnRETFpmAYOqAPB+06tlK7t27q9WWruPLBeVz193kc+Mv7+fS2ozjxU1sxecSAaodmVhU9eVDDHcB+wD051l0IXBYRWwLfAa6SG/mti40c1JcTP7UV9560Fyd8cgvunbOIT/30Hk67aTYrVq+tdnhmXa5mEpKksZLOl/SgpBWSQtLEPGXHSbpO0lJJyyT9SdL4zDIR8feIeCHHtsOBaSQ1JCLitnTVrmV9QWYFGtCnkRM+uSUzv7MnX/zgOC57YC77/exeZs1dXO3QzLpUzSQkYDJwCLAEuDdfIUn9gTuBrUma5A4naZa7S1IhvcPjgVcjYk3GsnnpcrOqGTagD2ceuD1Xz5jG2tbgCxc9yGk3zWbpijUdb2zWA9RSQronIkZGxH7Ate2UmwFsDhwYETdExI3A/sAE4CslHtvNdVYzPjxpKH874eMc9qHxXP7gXPY45y6ufGgea9e1Vjs0s4qqmYQUEYX+t+0PPBQRz2ds+xJwP3BAAdu/DIyW1Dtj2YR0uVlNGNCnkTMP3J6/fONjbDVqIKfe8CSfPf8+Hnl5SbVDM6uYmklIRdgWeDLH8tnAlI42johFwD+AowAk7U1SQ3o4V3lJx0maJWnWokWLSo3ZrCRTRg/i6hnT+NVhu7B05RoO+tUDnHz9Eyxd6WY863m6Y0JqIelnyrYYGNL2RNIpkhYAHwYulbRA0qh09fHA0ZKeA/4HOCwi95UrEXFxREyNiKnDh3d4w0OzspPEvttvym3f2p0v77YZV//jZfY5byb3P/9mtUMzK6vumJAg9yw5G/QDRcSZETE2IvpExLD094XpujkR8ZGI2DIidoqIf7R3MM9lZ7VgQJ9GTv3sFG78j48yoE8jX/r13zn7r0/z3tp11Q7NrCy6Y0JaQlJLyjaE3DWnTvMdY62WbD92MH/+xsc47EPjueieFznoVw/wytsrqx2WWad1x4Q0m6QfKdsU4KlKHNA1JKs1/Zp6ceaB23Px4bsy780VHPCL+/inr1uybq47JqSbgGmSNm9bkF5Au1u6ruxcQ7Jatc+2o7j+P3ZjYN/e/PslD3HNPz1Y1LqvmkpIkg6WdDDvz5qwb7ps94xilwBzgRslHSBpf+BGYD5wUYXicg3JatbkEQO44Wu7MW3zoZz0xyc497bnyDNGx6ymqZb+cCXlC2ZmROyRUW48cB7QNmT7DuCEiJhbyfimTp0as2ZVbA5Ws05Zu66V71//BH+YtYAvfnA8Zx64Hb0afM23VZ+khyNiakflamq274go6L8nIl4GDqpwOOtJmg5Mnzx5clcd0qxojb0a+K+DdmDYgD5ccPcLLFm+mp9/cWeaGmuqIcQsL/+lFsB9SNZdSOI7n96aH3x2Cn+bvZATr32M1tbaaQUxa09N1ZDMrDy+/NHNWLV2Hf/9t2dpaW7ih9On+DbqVvMKSkiSShm6E8BnIiLXND/dipvsrDv66u6TeOvd1fz6vpcYNqCJr++1RbVDMmtXoTWkscBfgUInc2sAvgQ0lRJUrYmIm4Gbp06dOqPasZgVShIn77cNS5av5pxbn2PskP4cuPOYaodlllcxTXZndDTFThtJjST3KTKzKmpoEP918A688vZKTvrj42w2rJkdx21S7bDMcip0UMOpJNf5FCQi1qbbvFJKULXG1yFZd9a7VwMXHLYLwwb04bgrZ/HGslXVDsksp5q6DqnW+Tok686eenUZB/3qAbbedCC/P24afRp7VTskqxOFXodU8rBvSSNK3dbMut6U0YM495AdefTltzn7r89UOxyzjXTmOiRXFcy6mX2335Sjd5vIZQ/M5dbZC6sdjtkG2h3UkM4Tl0/fMsdiZl3gu/tuzT/nLubb1z3OdmMGM3qTftUOyQzooA9J0jpgJlk3v0tNi4i6+EvOuA5pxpw5c6odjlmnvfTmcj7783vX3yK9sZcnbbHKKVcf0hzgmIjYM/sB1M39kz11kPU0mw1r5qzPbc8/5y7h1/e9VO1wzICOE9LlwLA86y4scyxm1oUO3HkMn9xmBD+7Yw6v+o6zVgPaTUgRcXZE/DPPurMqE5KZdZUfTt+W1gjOuLkiN1s2K4objs3q2LiW/nxjry342+yF3PXMG9UOx+pcUQlJkifCMuthZnxscyYNb+YHNz3JytXrqh2O1bFia0iPViSKGuepg6wna2ps4MwDt2f+4pX8auYL1Q7H6lixCakub6jiUXbW03140lCm7ziaC2e+wPzFK6odjtWpYhOSJ74z66G+v9/W9JL40Z89wMGqw4MazAyATQf34+t7TebWp17nnucKvfWZWfnUZUKS9GlJsyQ9LukhSTtWOyazWnDsxzZj4tD+nHbzbFavba12OFZn6i4hSRoCXAUcHhE7AN9Kn5vVvT6NvfjB9Cm8uGg5lz8wt9rhWJ0pNiEVfJO+cpI0VtL5kh6UtEJSSJqYp+w4SddJWippmaQ/SRqfUWQS8FZEPA0QEQ8A4yXtUvEXYtYN7LX1SHbfcjg/v3MOb737XrXDsTpSVEKKiF0rFUgHJgOHAEuAe/MVktQfuBPYGjiS5DbqWwB3SWpOi80BWiTtlm6zPzAQmFip4M26m1M+sw0rVq/jp7d7MmHrOgUlJEm/kfTxSgfTjnsiYmRE7Adc2065GcDmwIERcUNE3AjsD0wAvgIQEUuBzwNnSXoY+ATwFLCmki/ArDvZYuRADvvQeK76+zyee/2daodjdaLQGtKhJLWMlySdLmlyJYPKFhGF9q7uDzwUEc9nbPsScD9wQMayeyJij7TG9x1gNPB0GUM26/ZO+OSWNPdp5My/+F/DukahCWkkcCwwFzgFeFbSfZJmSKqlq0W3BZ7MsXw2MKXtiaRNM9adCtyZmcQySTouHZE3a9EiD4W1+tHS3MQ3P7EF9zy3iLue9Tx3VnkFJaSIeDcifpveB2kiyYd4C3AR8JqkqyXtK6nao/ZaSPqZsi0GhmQ8P0PSM5KeJ2nOOybfDiPi4oiYGhFThw8fXt5ozWrcER+eyPiW/pxzy7O0dzNPs3IoOoFExPyI+HFETAGmAb8h6Yf5M/CKpHPKHGOxcv3XbDDlUUTMiIitI2JyRBweEW+3t0PPZWf1qqmxgW/sNZnZry7jtqder3Y41sN1qkYTEf+IiK8DY4DzgBHA/ylHYCVaQlJLyjaE3DWngnguO6tnn9t5DBOH9uent89xLckqqlMJSdJkSWcAz5IkondIakzVMpukHynbFJKRdCVxDcnqWWOvBr6x1xY89doybnUtySqo6IQkqUXSVyU9SJKIvg88B3wJGBURM8ocYzFuAqZJ2rxtQXoB7W7pupK4hmT17oCdRrPZsGZ+5lqSVVCh1yH1lvQ5SdcDrwK/JLmY9HvA+Ij4dERcHRGrKhWopIMlHQy0XZy7b7ps94xil5CMBLxR0gHpRa83kswwcVEnju0aktW1pJY0madeW8Yts11LsspQId92JL0FbEIyWu33wOURMavCsWXHkC/QmRGxR0a58ST9WXuTDGa4AzghIuZ2NoapU6fGrFld+rLNasbada188tyZDBvQh+u++pFqh2PdiKSHI2JqR+UaC9zfvcDlwJ8joiozGkREQTcHjIiXgYPKeWxJ04Hpkyd36fXAZjWlsVcDX5o2gTP/8jRPv7aMbTYdVO2QrIcp9DqkAyPi+molo2pzH5JZ4uBdx9KnsYH/fWhetUOxHqikUXaSdk5n0X5T0tq2mbIl/VjSp8sbopnVik36NzF9x9Fc/+grvLOqLr+fWgWVMsruo8CDJDNq/y5rH63A8eUJrXZ4UIPZ+w6fNoEVq9dxw6OvVDsU62FKqSH9X+AWkut9vpW17hGgx91XyE12Zu/bcdwmbD9mMFc+NM9DwK2sSklIuwC/iuQvMfuv8U3AE76Z9XCHT5vAc6+/yz/nljwBitlGSklIq4D+edZtCvS4di032ZltaPqOoxnUt9GDG6ysSklI9wEnSOqVsaytpnQMyR1bexQ32ZltqF9TLz638xj+Nnshb69YXe1wrIcoJSGdStJs91j6ewBHSrqLZPbv08sXnpnVqkM+MI7Va1s9uMHKppTbTzwGfBx4HTiZZDaEr6erd4+IZ8sXnpnVqm1HD2a7MYO4ZtYCD26wsijpOqSIeCQiPkEyn91YYFBE7BkRj5Y1uhrhPiSz3A6dOo6nX1vGk68sq3Yo1gN09n5IqyLi1YhYUa6AapH7kMxy23+nMfRpbOCaWS9XOxTrAQqay07SXsXsNCJ63MAGM9vY4H692W/7Tbnx0Vc5eb8p9Gvq1fFGZnkUOrnq7bw/ki7fJKeRrgvAf5VmdeKQqeO4/tFX+H9Pvsbndxlb7XCsGys0IUFyN9g/po/llQnHzLqbaZu3MGFof657eIETknVKoQlpT+AIkts6fAG4nuSeSG6aM6tzkvjsDpty4cwXeXvFajbp31TtkKybKvT2EzMj4hhgFMnkqSOAWyS9LOlsSdtUMshq8yg7s/btM2UU61qDO595o9qhWDdW1Ci7dFTd7yJiX2A88DNgP+BJSb+oRIC1wKPszNq3/ZjBjBrUl1t9e3PrhM4M+34LmJs+AhhShnjMrBtqaBCfnDKCmc8tYtWaddUOx7qpUu6HtJukC4HXSG5r/i7wGeDwMsdmZt3IPlNGsXLNOu6b82a1Q7FuqqCEJGmypNMlvQDcA2wFnAiMiojDIuKWiGitZKDlJGk/SQ9LelTSE5KOqHZMZt3dtM2HMrBPI7c95WY7K02ho+yeA5YBfwKOBdrmnB8haUR24Yh4sTzhlZ+kBpI73X4kIp6SNAF4TtKfIuLdKodn1m01NTaw59YjuP3p11nXGvRqyHfJolluxTTZDQKOIrlIdk4Hj7KSNFbS+ZIelLRCUkiamKfsOEnXSVoqaZmkP0kan1kk/dmS/tyEpD/Mc+ibddI+247kreWreeRl37jPildoDenoikbRscnAIcDDwL3APrkKSepPcj+m94AjSQZbnAncJWmHiFgeEeskfQG4QdJyksEYn48IJySzTtp9y+E09Wrg1tkL+cDElo43MMtQUEKKiMsrHUgH7omIkQCSjiVPQgJmAJsDW0XE82n5x0lqbV8BzpXUSHLbjM9HxD2SPgDcmCYs98aadcLAvr358KSh3PnMG5z8mSnVDse6mU7N9t1VihgwsT/wUFsySrd9CbgfOCBdtBMwOiLuSdf/E3gF2Ll8EZvVrw9u1sILi5azZLkbHaw43SIhFWFb4Mkcy2cDbV/X5gOjJU2BZAQhSZOgbyxoVga7TkguSXx0vvuRrDg9LSG1ALn+CxaTXrgbEa+TNO1dI+kxknn5/iMict7QRdJxkmZJmrVo0aIKhW3Wc+wwdjC9GsQj896udijWzRQz23d3keteyhuMP42Iq4GrC9pZxMXAxQBTp071fZrNOtC/qZEpmw7i4XmuIVlxeloNaQnvD+fONITcNaeCeHJVs+LsOmEIjy14m7Xrus318lYDSk5IuS6IrQGzSfqRsk0BnuriWMzq1s7jN2HF6nU8s/Cdaodi3UhnakizyhZF+dwETJO0eduC9ALa3dJ1JfFs32bFWT+wwRfIWhHa7UOStH87q/uWOZZ2STo4/XXX9Oe+khYBiyJiZrrsEuDrJNcVnULSn/QjkpF1F3Xi2NOB6ZMnTy51F2Z1Zcwm/RgxsA8Pz1vC4R+eWO1wrJvoaFDD9cBMsgYFpAaWP5x2XZv1/IL050xgD4CIWC5pL+A84EqSuO8ATujMPHURcTNw89SpU2eUug+zeiKJXScM4ZGXPdLOCtdRQpoDHJNeXLoBSfMrE1JuEVHQTI3p8O2DKhyOmXVgl/FD+H9PLmTRO+8xfGCfaodj3UBHfUiXA8PyrLuwzLHULI+yMyveLmk/kidatUK1m5Ai4ux0ap1c686qTEi1x4MazIq33ZhBNPVq4BFfj2QF6mnXIVWEa0hmxevT2IvtxvgCWStcUQlJ0phKBVLLXEMyK80OYzdh9qvLaG31JCfWsWJrSI9WJAoz65G2GDmAlWvWsXDZqmqHYt1AsQnJ9yQ2s4JNGj4AgOffKPmqC6sjxSakuqx3uw/JrDRtCemFRU5I1jEPaiiA+5DMSjNsQBOD+/V2QrKCOCGZWcVIYtLwZjfZWUGckMysoiYNH8ALi5ZXOwzrBopNSF06XVCtcB+SWekmjxjAonfeY+nKNdUOxWpcUQkpInbtuFTP4z4ks9J5YIMVqugmO0kNkhqzln1K0n9K2rl8oZlZTzBpRJqQ3I9kHehotu9crgbeA44AkHQ8798KYo2kz0TE7WWKz8y6uXFD+tHUq4HnXUOyDpQyqGEa8NeM598GLgUGA38CTi5DXGbWQzT2amDisP688IYHNlj7SklII4BXACRNBjYDfhER7wC/BbYvX3hm1hNMGj6AF11Dsg6UkpCWAUPT3/cA3oyIx9Pn6+jiW5t3BY+yM+ucySMGMG/xClavba12KFbDSklIDwDflfRZ4AQ2bL6bDCwoR2C1xKPszDpn0vABrGsN5r3lZjvLr5SE9B2gBbiJpDZ0Wsa6Q4EHOx+WmfUknmTVClH0KLuImANsKWloRLyVtfqbwMKyRGZmPcbmw5sBX4tk7Stl2DcAOZIREfFE58LpGpJGs2FTYzPJ4IwREbG4OlGZ9VzNfRoZPbivpxCydpWckLqziHgV2KntuaTvAh9xMjKrnEkjBrjJztrVbSZXlTRW0vmSHpS0QlJImpin7DhJ10laKmmZpD9JGt/O7r8M/LoScZtZwkO/rSPdJiGRjOA7BFgC3JuvkKT+wJ3A1sCRwOHAFsBdkppzlP84MBD4SwViNrPUiEF9WL56HStXr6t2KFajulOT3T0RMRJA0rHAPnnKzQA2B7aKiOfT8o8Dc4CvAOdmlT8GuDwi1lYkajMDoKV/EwCLV6xmTFO/Kkdjtajb1JAiotAr6vYHHmpLRum2LwH3AwdkFpQ0CDgIN9eZVVxLc5KQlixfXeVIrFZ1m4RUhG2BJ3Msnw1MyVr2ReDhdCi7mVVQW0J6ywnJ8mg3IUnaTNIdkl6UdK6kvhnr/lH58ErSQtLPlG0xMCRr2TEkE8PmJek4SbMkzVq0aFGZQjSrP0NcQ7IOdFRDuoBkBu8vkHzQ3yFpYLqudyUD66TIsUwbFYr4YERc2e6OIi4GTgceaWpqKlN4ZvVnfR+SE5Ll0VFCGhkRv4yIhyPiKODPJElpMLk/9GvBEpLkmW0IuWtOHfJcdmadN7hfbxoES1Y4IVluHY2y65P5JCLOlrQGuINkqHQtmk3Sj5RtCvBUKTuUNB2YPnny5M7EZVbXGhrEkP5NriFZXh3VkJ6TtHfmgog4B/gdMKliUXXOTcA0SZu3LUgvoN0tXVc015DMymNIsxOS5ddRDenfci2MiHMlXVOBeNol6eD0113Tn/tKWgQsioiZ6bJLgK8DN0o6haRp8UfAfOCiEo/rGpJZGbS4hmTtaLeGFBHvRcR7bc8ljclY90olA8vj2vRxfPr8gvT56RlxLQf2Ap4DrgSuAl4C9oqIkuYtcQ3JrDyGNPd2H5LlVexMDY+S3MK8KiJio5Fyecq9THLBa1m4hmRWHi3NTTw87+1qh2E1qtgLYwtKCD2Na0hm5dHS3MSSFauJqNVBulZNxSYk/xWZWcmG9G9iXWuwbKWnjrSN9cSpg8pO0nRJFy9durTaoZh1a23TBy12P5Ll4IRUADfZmZVH2/RBHmlnuTghmVmXGer57KwdxSak+RWJwszqwpD+brKz/IpKSBGxa8eleh73IZmVR4ub7KwdbrIrgPuQzMqjf1Mvmhob3GRnOZV8C3NJ44BxQN/sdRFxZ2eCMrOeSZKnD7K8ik5I6aSlVwEfbFuU/oz09wB6lSU6M+tx2i6ONctWSg3pUmA8cALwDNDj/7I8dZBZ+bR4xm/Lo5SE9AHgqIj4Y7mDqVURcTNw89SpU2dUOxaz7m5IcxMLlqyodhhWg0oZ1LCAOqgVmVlltPTv7RqS5VRKQvoxcJKk5nIHY2Y935DmJpatWsuada3VDsVqTNFNdhFxpaStgbmSHgKWbFwkjixLdGbW47TN1vD2ijUMH9inytFYLSlllN1RwPeAdcAubNx85xnBzSyvtvnslqxY7YRkGyhlUMPpwPXAMRFRF3fa8ig7s/JpSacPeuvd1TCyysFYTSmlD2kocEG9JCPwTA1m5ZRZQzLLVEpCug/YptyBmFl98HzZCLmrAAAVlUlEQVR2lk8pTXbfBP4gaQnwNzYe1EBEePiMmeXUNuO357OzbKUkpKfTn1fkWR8l7rfLSGoC/hv4DMmgjHkRsV91ozKrD02NDQzs0+hbUNhGSkkcZ9D9R9L9GGgCtoqIVkmbVjsgs3oyxNMHWQ6lXId0WgXiaJekscBJwFRgR6AfsFlEzM1RdhxwHrA3yWSvtwMnRMTL6fr+wHHA2LamxYh4rQtehpmlnJAsl+5yP6TJwCEk/VX35iuUJps7ga2BI4HDgS2AuzJmlpic7ue7kv4p6QFJB1QyeDPbUEv/3h5lZxspOiFJOknS+XnW/VzStzsf1kbuiYiRaT/Pte2UmwFsDhwYETdExI3A/sAE4Ctpmd4ks5W/EBEfAL4MXCppUgXiNrMcWpr7sGT5mmqHYTWmlBrS0cDjedb9K11fVkWM2tsfeCgins/Y9iXgfqCtFjSPpA/synT9MyRx71y2gM2sXS3NnmDVNlZKQhoPzMmz7kWS2ki1bAs8mWP5bGAKQES8CdwCfBogHdCwHfBErh1KOk7SLEmzFi1aVJGgzepNS3MfVq5Zx/L31lY7FKshpSSkFcCYPOvGAu+VHk6ntZDjuihgMTAk4/lXgW9KeoLkWqoTI+LZXDuMiIsjYmpETB0+fHjZAzarR2OH9ANgvu+LZBlKSUj3At+WtMGsiOnz/6SdQQddJNeQdG1QIGJuRHwiIraPiB0j4qr2dihpuqSLly5dWtZAzerVxKHJGKO5bzoh2ftKuQ7pNOAB4DlJ/wu8QlJj+hLJPHdHlSu4EiwhqSVlG0LumlNBfMdYs/IaP7Q/APPeWl7lSKyWlHId0mOS9gTOIbk2qAFoJZnj7qCIeKy8IRZlNkk/UrYpwFOl7tSzfZuV1+B+vWlpbmLuW64h2ftKug4pIv4RER8HBpL0Gw2MiD0iYlZZoyveTcA0SZu3LZA0EdgtXVcSz/ZtVn4ThvZ3Dck20KkLYyNiZUS8GhEryxVQPpIOlnQwsGu6aN902e4ZxS4B5gI3SjpA0v7AjcB84KJOHNt9SGZlNnFoM/NcQ7IMBSUkSXemty0viKSGdJstSg9tI9emj+PT5xekz09vKxARy4G9gOdIrjO6CngJ2Csi3i31wK4hmZXfhKH9eXXpSlatWVftUKxGFNqHtAdJ81yhVMI27YoIdVwK0jnrDirXccF9SGaVMHFoMxGwYMkKJo8o20eFdWPFDGq4QVIx1xh19xnB1/MoO7Pym5COtJv7phOSJQpNSJeXuP83S9zOzHq49dcieWCDpQpKSBFR9vnpuhM32ZmV3yb9ezOob6MHNth63eX2E1XlQQ1m5SeJicOaXUOy9ZyQzKxqJnjot2VwQiqAr0Myq4yJQ/uzYMkKVq8t9A4z1pM5IRXATXZmlTFhaDOtAa+8XfFr660bcEIys6qZ2Db02/1IhhOSmVXRhHTo97w3nZCsjAlJ0hWS5ki6O50N3MysXcMGNNHc1MuzfhtQ3hrS9cCewFnAeZKOKeO+q8qDGswqo23ot2f9NihjQoqI6yNiQUTcBnwcOFrSoeXafzV5UINZ5XjWb2tTzia7iZI+Jukw4Gsks2xfVq79m1nPNGFof+YvWcHadR76Xe9KuYV5Po8BL5Pc0nw+MIckMZmZ5TVhaH/WrAteW7qKcS39qx2OVVHZElJEuD3LzIo2anA/ABYuc0Kqd0UnJEn9gH8DdiK579G/gN9HhBuBzaxomw7uC8DCpauqHIlVW1EJSdJ2wF+AcVmrfihpv4iYXbbIaohn+zarnJGDnJAsUfCgBkm9gT+k2xwMDCK5I+zn0yJ/kFTOPqma4VF2ZpUzqG8j/Xr3YuEyJ6R6V0wCORTYApgaEY9lLL9B0ovALODfgSvKGJ+Z9XCS2HRwXyckK6yGJKkB2B14AHhcUkPmA3gSuA/YQ5IqF66Z9UQjB/V1k50VXENaC0TG7/nsARxRxH6rRtJc4D2gbZrhX0TEpdWLyKx+jRrcl3+8tLjaYViVFZo4zgAOSsv/Psd6AV9If/9DGeLqKodGxL+qHYRZvRs1uC+vL1tFa2vQ0OBGlnpVUEKKiNMkrQOOB34cEWsy16eDGY4BLo2I08sfJkgaC5wETAV2BPoBm0XE3BxlxwHnAXuTJMvbgRMi4uVKxGZmnTNqUF/WtgZvLV/N8IF9qh2OVUkxUwf9ERgO/HeOdf8XGAlcV46g8pgMHAIsAe7NV0hSf+BOYGvgSOBwksEYd0lqzip+haQn0pnKx1QmbDPryKj0WqTXPbChrhXc1xMRT0k6B/huenuJW0n6lfYmqbGcXeHrkO6JiJEAko4F9slTbgawObBVRDyfln+cZCqjrwDnpuV2j4h5ae3u+yTJ9MMVjN/M8hiVXov02tJVbDfGl1fUq6ImV42I7wMnAEOAE4Fvk9SavhURp5Q/vA2OXejMi/sDD7Ulo3Tbl4D7gQMyls1Lf64lad77UHqtlZl1sbYakod+17eiZ/uOiJ9HxARgBDAiIsZFxM/KH1rJtiUZhp5tNjAFQFKzpE0y1h0GPJndN5aWPU7SLEmzFi1aVJGAzerdsAF96NUgFi5d2XFh67FKHp4dEW+WM5AyaiHpZ8q2mKRmB0l/1x8l9SIZ9DCf90cJbiAiLgYuBpg6dWrkKmNmndOrQYwY2IeFS9+rdihWRTV/vVCJciWO9WNJI+JFYOdCd+a57Mwqr23ot9Wvct7CvFYsIaklZRtC7ppThzyXnVnljRrUl9fcZFfXemJCmk3Sj5RtCvBUKTuUNF3SxUuXLu1UYGaW38hBfXl9mZvs6llPTEg3AdMkbd62QNJEYLd0XdFcQzKrvE0H9+Xd99byzqqNxhZZnehWCUnSwZIOBnZNF+2bLts9o9glwFzgRkkHSNofuJFk4MJFJR7XNSSzCvPFsdatEhJwbfo4Pn1+Qfp8/XRFEbEc2At4DrgSuAp4CdgrIt4t5aCuIZlV3vs36nOzXb3qVqPsIqKgWRfTOesOKtdxPcrOrPLabmXugQ31q7vVkKrCNSSzymurIbnJrn45IZlZTejbuxdD+vf29EF1zAmpAB7UYNY1fOfY+uaEVAA32Zl1jVGD+7qGVMeckMysZmw62DWkeuaEVAA32Zl1jZGD+vLmu6tZvbbQu81YT+KEVAA32Zl1jbYb9b3xjmtJ9cgJycxqxpDmJgDeXuHpg+qRE5KZ1YwGFXTtu/VQTkgFcB+SmVnlOSEVwH1IZmaV54RkZmY1wQnJzMxqghOSmZnVBCckMzOrCU5IBfAoOzOzynNCKoBH2ZmZVZ4iotoxdBuSlgJzcqwaDGRXn3ItGwa8WYHQOpIrlq7aTyHbdFSmvfX51tX6OYHynJdKnZNCylXqvHT3c1Lqfnry/8qEiBjeYamI8KPAB3BxocvzLJtVS3F3xX4K2aajMu2t767npFznpVLnpJrnpbufk0qel57+v+Imu+LcXMTyfGWroVyxlLKfQrbpqEx767vrOYHyxFOpc1JIuZ54Xvy/UngsZecmuy4kaVZETK12HPY+n5Pa43NSm7rivLiG1LUurnYAthGfk9rjc1KbKn5eXEMyM7Oa4BqSmZnVBCckMzOrCU5INUzS9yU9K6lV0oHVjqceSZok6T5Jz0l6VJI726vM/xe1R1JfSTdIelrSvyTdImnzYvfjhFTb7gD2A+6pdiB17ELgsojYEvgOcJXk25pWmf8vatOvImKbiNiJZNj4pcXuwAmpQJLGSjpf0oOSVkgKSRPzlB0n6TpJSyUtk/QnSeOLPWZE/D0iXuhs7PWknOdJ0nBgGnA5QETclq7atcIvo0cp9/+O/y/Ko5znJSJWRcQtGZs8BLiGVEGTgUOAJcC9+QpJ6g/cCWwNHAkcDmwB3CWpuQvirHflPE/jgVcjYk3GpvPS5VY4/+/Upkqel28ANxYbUGOxG9SxeyJiJICkY4F98pSbQfLNYKuIeD4t/zjJHHhfAc5Nlz1C/g+2nSNifhljrydlPU85uLmueJU+J1aaipwXSd8DtgQ+UWxAriEVKCJaCyy6P/BQ24lLt30JuB84IGPZLhExLM/DyahEZT5PLwOjJfXO2G5CutwKVO7/HSuPSpwXSScCBwH7RsSKYmNyQiq/bYEncyyfDUzp4lgsvw7PU0QsAv4BHAUgaW+SGtLDXRNi3fH/Tm0q6LxI+hbwRWDviHi7lAM5IZVfC0mbbLbFwJBidiTpFEkLgA8Dl0paIGlUGWK0ws/T8cDRkp4D/gc4LDy9SaUUdE78f9HlOjwvksYCPwE2Ielb+pekWcUeyH1IlZHrA6vovoeIOBM4s/PhWB4dnqeImAN8pGvCMQo7J/6/6HrtnpeIWEAZ+lddQyq/JSTfKLINIfe3DKsOn6fa43NSm7rsvDghld9skjbXbFOAp7o4FsvP56n2+JzUpi47L05I5XcTMC1z2oz0YrPd0nVWG3yeao/PSW3qsvPi208UQdLB6a+fIOns/hqwCFgUETPTMs3AY8BK4BSSttcfAQOBHSLi3a6Ou974PNUen5PaVGvnxQmpCJLyvVkzI2KPjHLjgfOAtmHCdwAnRMTcSsdoPk+1yOekNtXaeXFCMjOzmuA+JDMzqwlOSGZmVhOckMzMrCY4IZmZWU1wQjIzs5rghGRmZjXBCcnMzGqCE5KZmdUEJyQzM6sJTkhWFElHSQpJb0sakrWuMV13WoWOfWB6V8pitjmtnelRrAtImpH+XeR6LJfUq4R9Xpaxj7szlh+VsXzLHNvtkbH+k0Ue80ZJiyX1ybN+YPp6Lkufn5JxrAXFvcL65IRkpRoMnNTFxzwQKCohAZeS3FnUqmcnYDnJech+fCQi1pW434XpPr6WY907wOE5lh+RrivF5ST3APpsnvUHA/3TcgC/TeP7a4nHqztOSFaqW4Fv1PKtoyX1iYgFEfFQtWOppnzf6LvQTsDsiHgox+OxTuz3vXQfue7J8yfgS5LW38VUUj/gIOCPJR7vz8BbJEktlyOAl4G7ASLilfRvb1GJx6s7TkhWqrZbSJ/cXiFJO0q6SdISSSsl3S/pY1llLpM0N8e2d7c1x6TNIEcCYzKaQeZmlD0tXbadpFskvQv8IVeTXUbZLST9RdK7kuZJ+oGkhqyyn5b0YBr7Ukk3SNqqkDeowNdeTCzF7G+D9yFd90VJz0haJekJSftnvccHp9vumOdcPFjI687aTsAOwOPFbttJVwITgI9mLPsc0Is8Camj9zciVgO/B/aVNCxr2/HA7sCV4RmrS+aEZKV6DfgFcJykCbkKSNoFeIDk9sczSL6dvgXcLmnXIo/3I5Kmj0W839zzuRzlbgRmAvuTTJffnuuBO0maAm8ATidJem3xfxr4C/AucCjwVWA74D5JY9rbcQmvvaNYit3fBu+DpL2Bq4Bn0m3PAX4KZPaz3AC8Cnwl67VsRfJhe1F7rzmPLYABwGwlfYyZj6L7joowD7iHDZvtjiB5nze6f08R7+/lQG+Sv4dMXyK5LcMVZYq/PkWEH34U/ACOIrlB12SSf963gd+k6xrTdaelz+8AngaaMrbvlS67IWPZZcDcHMe6G7g7q9yCPHGdlh77m7mW5yl7dNbyJ4BbM57PAuYAjRnLNgPWAOd28D4V+toLjaXY/WW/Dw8AT5LeciZdtkta9u6s7ZcCzRnLzgWWAP1K+Hs5JD1Grsejnfg7zPc3k/n3+eU07r7ApsBakvv57JGW+WSx72+6fDbw96xlTwMPtBNrzr9bPzZ8uIZkJYuIxcBPgCOym7HS9vrdgWuB1rZvxSTfIm8HPl6hsK4vouxfsp4/CYyH9XfJ3AW4JiLWthWIiJeA+0leW04lvvb2Yillf+vfh7QmMhX4Y6SfkOlreQR4KWu7i0k65r+YbtuXpKZ2RUSszPea27EzyYf/bsAHsh4Ht7NdOVwL9AGmA4eRDIK4I7tQCe/vFcAHlY7ik/RBYGtcO+o0JyTrrPOAxcAZWctbSL5hnkpSo8h8fB0Ykt1HUiavFVF2cdbz90i+TUMymkp59reQ5PXlU8prby+WUvaXGfcwkmamN3LE+nrmk4h4laS57/h00RfS45fSXAfJgIYXI+KBiJiV9XihxH0WJCLeIWmGPJykue6qiGjNUbTY9/d/gVbeH9xwBMn5uqYCL6OuNFY7AOveIuJdSWeT1JT+J2PV2yT/tL8kzzfHjA+HVUBTjiJDSdrxiwqpyPL5LEn3lWsU4Sjaj6uY116IUvaX+T68SfLhOiLHpiNJRoZlugC4I+07+Qpwb+QeyVaInYAORzlKWgbsFBEvSvo88LmIyDVsu1hXkNQ+G0hrfTkU9f5GxCuSbicZxXcGSX/STRGxpAzx1jUnJCuHC0iuD2obeUdELJd0L7Aj8EgHH8DzgJGShkXEmwCSJgFbkfR9tHkP6Ffu4HNJ438Y+IKk0yK9ViYdwPER4PwOti30tRcaS8n7i4h1kmYBB6WvJQDShLMZWQkpIu6U9DRJ39FuJM1dRZM0kiR5FzLCbj4wDngR+CcbfrnpjNtIRhm+HRGzcxUo8f29nGSQyNkkNVA315WBE5J1WkS8l35TvDhr1bdIRjrdIunXJM1Iw0j6ZnpFxHfTcteSjKK7StK5aZnvkXyzz/QU0CLpqyQDDlZFxBOVeE2pU0m+Xf9Z0gUko8VOJ+n0/0kH2xb62gvV2f39kOTaseslXZxuexpJ82OuD+ALgZ+RnIONhklLmkjS/3R6RJyW55g7pz/XSpqWY/3jEbEi/X0BaZ8ZSX/MK+2+mgKlXyTy1YwyFfv+Xg8sA/4PSVPo38oRb71zH5KVy29JRqStl3aaf4CkeevnJB+IPwO2J/nnbyv3PEkH9xiSNv/vkHxAPJd1jEtJrgP5MfAP4OYKvI7M+P8GfAbYhORb9oUko6k+mva1tLdtQa+9iFg6tb+IuI2kprMNyYfpScB/kiSkpTk2uTb9eVlEvJdjfXP6c2E7h90p/XkG8GDW4wE2bKadD4yTNJhkSHp2n2RFFfv+pgM8riXpZ/xd5sAXK50yBt2YWR2RNBZ4HjgrIn6UtW4GyUCGLdMvDNnbHgecBUzIqOV0JpYfklyztBlweURk17azy19GMnx7Msmw/lKnH6qY9KLgXsCvgU9ExNgqh1TzXEMyqwOS+kn6laSDJO0u6WiS/pUVJDXPtnJTJE0naZq8IVcySu0OnFeOZJSaTzKi75cdJaMME0gGa2w0lLtGnEwSX76phiyLa0hmdUBSE8mw5GkkoxeXA/cC34+IJzPK3U0yaOMB4N87apqslrQPq236nnci4tnqRZObpE1JmqEBVkdEV0+f1O04IZmZWU1wk52ZmdUEJyQzM6sJTkhmZlYTnJDMzKwmOCGZmVlNcEIyM7Oa4IRkZmY1wQnJzMxqwv8PWJjr6dCKk6UAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Initialise neutrino_flux interpolation function, using default settings\n", "CEvNS.loadNeutrinoFlux(source=\"CHOOZ\")\n", "\n", "#Plot neutrino flux\n", "E_nu = np.logspace(-1, 2,100)\n", "\n", "pl.figure()\n", "pl.loglog(E_nu, CEvNS.neutrino_flux_tot(E_nu))\n", "\n", "pl.title(r\"Neutrino flux from Chooz\", fontsize=12)\n", "pl.xlabel(r\"Neutrino energy, $E_\\nu$ [MeV]\")\n", "pl.ylabel(r\"$\\Phi_\\nu$ [cm$^{-2}$ s$^{-1}$ MeV$^{-1}$]\")\n", "pl.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we calculate (and plot) a whole bunch of recoil rates.\n", "\n", "Choose which rates to plot by specifying `SM`, `magnetic`, `scalar`, `vector` in `interactions`:" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "interactions = [\"SM\", \"scalar\"]" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAisAAAFTCAYAAAD4N0wZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3Xl8TFf/wPHPSSIRCUISsogIEUFiaUNpKEXVEpEEJfb2x1Nb0fbp011RbbVqqV310YUKtUtJPaXV2osmIUEstYTYgywESc7vjzsZWSbJDIkszvv1ui+Ze8495zt3xsyZe88ipJQoiqIoiqKUVmYlHYCiKIqiKEpBVGNFURRFUZRSTTVWFEVRFEUp1VRjRVEURVGUUk01VhRFURRFKdVUY0VRFEVRlFJNNVYURVEURSnVVGNFURRFUZRSTTVWFEUpMkKIWCFE+5KOo7QQQnwmhBhvQv4zQohOxRmTkj8hxF9CiMYlHYeSl2qsKI+VEKK/EOKAECJFCHFRCBEhhGijSzsjhLijS8va5hpR5hkhxGUhhE22fcOEENuzPW4jhNgthLglhEgUQuwSQrQolif5iIQQ/YQQ+4QQqUKIK7q/RwkhREnHVhgpZWMp5fbirqekv9SNqV8I4QgMBhZl25fv+780yPV/8JIQ4jshhO1jqrc0NNK+BCaXdBBKXqqxojw2Qog3gFnAp0BNoDYwH+iZLVsPKaVttm2MkcVbAOPyqbcK8DMwB6gOuAKTgLsP9USKkRDiTeArYBrghHaeRgD+gKWJZVkUeYDlwGM8L0OBzVLKO7p6jXn/lwY9pJS2QDOgOfBuCcdToCJ+PTcCzwshnIuwTKUoSCnVprZi34CqQArQp4A8Z4BO+aSdBrrq/rYDrgG22Y57B0gE7HT7hgHbdX/7ATdNjLchsB24CcQCgbni/DdwCLgFrAQq5lPOO8ApIBk4AgQXco5SgV6FxOYCrAGu6s7L2Fyxva2L7S5aI+4M8JZuXyrwX7QvywhdXFuBasbEXNhzz/4aAhLwzJb2HTAlV16j48p23FIgE7ije0/9x8i4c5+Xp4BIXf5VuucypbDznF/9BuL8DRho7Pu/sPNLAe/JgtKBvrq6s7a76P5vFPZ/EPgC2GTC+6+g18ANWKs77jowt5DXs7D/gzleT2M/M3THXdDFGAd0zHXcr8CQh/2sU1vxbCUegNqejA3oAqQb+lDJlifHB2WutF3AsGyPd/DgS/EM0En3QThFty97Y6WK7sPxe6ArBr4Ac9VVATgJvId2NaOD7oOtQbb6/tJ9aFcHjgIj8imrjy6fme5LIxVwfoRzZAYcBCboYqsL/AO8mC22KN0Xg3W2fXvRGgKuwBXgb7RfzVZoX6ofGRNzYc8d0xsrRsdV2HvFiLj150V37s6iXY2rAIQA97K9f4w5zwbfq9niuQq0MPa1Lej8Uvh7ssD0bOVX0ZX5amHnFagFHAa+MuH9Z/A1AMyBaGAmYANUBNrk93oa8XxzvJ4mfGa8CsQDLrp9dYB6uY6bDcwoys8/tT36pm4DKY+LPXBNSpleSL71Qoib2bbhuv0JaF9oWSoCuftwTABe0/UV0JNSJgFt0L48FwNXhRAbhRA184mhFdovsKlSyntSyt/QbiOFZsszW0qZIKVMBMLRLpnnIaVcpcuXKaVcCZwAWuZTrwO5zpGun81NXT+C54AWgKOUcrIutn90z6lfrtjipe72g84cKeVlKeUFtA/tfVLKSCnlXWAdWgPB2JiNeu5GMjquwhgZd9Z5aYV2dWW2lPK+lHItWiMhizHnuTB2aF+wYPz7PyvO3Oe3sPdkoe9ZIYQZsBytEb+I/K0XQiSjfalfAT7KllbgeSngNWiJ1oh5S0qZKqVMk1LuLCAGY/8P5n6fZ2foMyMTrSHcSAhRQUp5Rkp5KtdxyWivnVKKqHvayuNyHXAQQlgU8oEdJKXcamD/BXQfPEIIP7RfdDuyZ5BSxgghfka7FH00V9pRtD4ECCG8gWVo/Qeyf/hlcQHipZSZ2fadJecH36Vsf9/WHZOHEGIw8AbaLzjQPoAdDOXFwDmSUj6rK+c82q9Vd8BFCHEz23Hm5DwX8QbKvpzt7zsGHus7URoRs1HP3UhGx1UYI+LOfl5cgAtSaj+lDaQbc54LcwOorPvb2Pc/GD6/hb0njXnPfqKLZ2wh9QdJKbcKIdqhNW4c0G7FQCHnpYDXwAo4a2RjzdjnY+h9np2hz4ylaFd7JgKNhRBbgDeklAnZjqvMg+erlBLqyoryuOwB0oCghzw+AaglhKiM9oEzWkqZZiDfR8Bwcn6o5SClPIZ2S8KngLrcdL9Es9RG+/AzmhDCHe1X5xjAXkppB8SQ94pQlj1o998L6nAZD5yWUtpl2ypLKbtlyyPzO7gYYi7IbaBStsdODxuXATmeo5FxZz/mIuCaa4SVW7a/CzvPxpzjQ4CX7u+ieP8X9J4sMF0I0Q+tYd5bSnnfmAqllH+g/T/5MtvufM9LIa9BPFC7gM6wuc+nMf8HC3sNDH5mSCmXSynboDW8JPB5ruMaot2yUkoR1VhRHgsp5S202zTzhBBBQohKQogKQoiuQogvjCgiAe2X0TpgmZRyRT71nETrlKj/9SiE8BZCvCmEqKV77Ib2wb03n7r2of36+o8uxvZAD8BgnQWwQfswvKqr92XybyAhpbyJNkppvhCitxDCVghhJoRopisLtFsVSUKIt4UQ1kIIcyGETxEOwzYp5kJEAf11MXYB2hVNiIB2BaZutsemxr0HyADGCCEshBA9yXnLqLDznLt+Qzaje85F8P4v7D2Zb7oQojnaSLggKeVVI+rKbhbwgu49CAWfl4Jeg7/QGohThRA2QoiKQgj/bPXkPp9F8X8wz2eGEKKBEKKDEMIKrfF4B+19gC5mK+BptE62SinyxDZWhBDbhBDRQogoIcSObP8ZlWIipZyBdon4A7QPtHi0X2Hrs2ULFznnWVmn238BaAD8IaX8pJCqJvPgyx20e9DPAPuEEKlojZQY4M184rwHBKJ1xr2GNrx0sO6KjNGklEeA6WhfjJcBX7ROfwUd8wXaOfoPWn+By2jzdLwN7JZSZqB9aDdDG+1wDfgGbbTJI3uYmAswDi3Wm8AAcr7Oj+oz4ANdf55/mxq37jUOAf5PF99AtD4Rd3XphZ3nHPXnU80PQDchhLWuTGPe/wXFm+97spD0nkA1YGe2/1cRhdWpK/eq7nl8WNh5Keg1yHacJ3AOOI/WATdL7tezKP4PGvrMsAKm6sq8BNRA68SbJRCtT0/220JKKSBy3rJ9cgghqup+7SCECEYbdaAaLIryhBJC7AMWSim/LcIyPwWuSClnFVWZSvHRvQf+T0oZU9KxKDmVmSsrQohaQog5Qog9QojbQggphKiTT143IcRqoc1WmiSEWCuEqJ09T1ZDRadKMYauKEopJIRoJ4Rw0t0GGgI0AX4pyjqklO+phkrZIaV8RjVUSqeyNBrIE3gJbYz/DqCzoUxCiEpo8zPcBYag3UOdAvwuhGgipUzNlvdHtHvKmUA3A8UpilJ+NQB+Qhuxcgqt8+nFkg1JURRDysxtICGEWdYwNiHEMLRe5x5SyjO58o0DZqBNHnRSt88Dbbz/f3T3jXOXPQxtpsXuxfssFEVRFEUxVZm5DZRrvH1BAoG9WQ0V3bGn0Tp65Tck9L9oPd7tHy1KRVEURVGKWlm6DWSsxsAGA/tj0aaCRghRDW29jaxLvr3QRl4k5j5ICPEv4F8ANjY2T3t7exdHzIpSpty+f7ukQ1DKoUoVKhWeSSn1Dh48eE1K6Vh4TuOVx8ZKdbSZI3NLRBu+h+7flUKIrOmXrwAB0sA9MSnl18DXAH5+fvLAgQPFErSilCUHEw6WdAhKOfS0y9MlHYJSBIQQZ4u6zPLYWAHDMxvqZ6rUrWdRVJNoKYqiKIpSjMpMnxUT3EC7upJbNQxfcSmUEKKHEOLrW7duFZ5ZURRFUZQiVR4bK7Fo/VZyawQceZgCpZThUsp/Va1aJJOEKoqiKIpigvLYWNkItBJC6NeZ0E0e569LM5m6sqIoiqIoJafMzLMCIITorfuzIzACGIW2xsZV3QqhCCFs0FbMvIO2BocEPkZb9ruJlDLlYetXHWwVRVPmOthKMEsxwyzD7BHWpFaKW0WLiiUdgmIEIQSVKlXC1dUVIfIuyC6EOCil9CvKOstaB9tVuR7P1/37B9AeQEqZKoToAMxEWxZcANuA8Y/SUFEUpQy7A5WtKuNY0xEzs/J4Qbl8qGxVuaRDUIyQmZnJxYsXuX79Og4ODo+lzjLVWJFS5m3CGc53Dm3ulCIhhOgB9PD09CyqIhVFeYws7ltg72SvGiqKUgTMzMxwdHQkPj7+sTVW1P9cI6gOtopSxmWCuYV5SUehKOVGhQoVyMjIeGz1qcaKoihPBEP31hVFeThCCB5nn1fVWFEURVEUpVRTjRUjqKHLiqIoilJyVGPFCFl9VqxtVE91RVHKlxHDRjD5o8mltjxjDB06lA8++MCovHXq1GHr1q3FHJFS1AptrAghzj3EdlYI4fM4nsDjdOJKCsv3nXus9+kURSn/9uzaQ6f2nahVoxa1nWvzQvsXOHhAm8vGx8uH37f9XsIRFg0fLx/sK9tz/dr1HPv9W/pTpWIVzpw5UzKBZbN8+XL8/PywtbXF2dmZrl27snPnTkBr6FhbW2Nra6vfxowZY1S5derUoWbNmqSmpur3ffPNN7Rv317/eOfOnTz77LNUrVqV6tWr4+/vz/79+4v0+ZVVxgxdrgVsRpt8zRhmwEDA8mGDKq0qWZrz3rrDbD16mam9fKlRWU1gpCjKo0lKSuKlkJeYMXsGIb1DuHfvHrt37sbKyqqkQytQeno6Fhamz37hXsedVT+tYsSoEQDExsRy586dog7vocyYMYOpU6eycOFCXnzxRSwtLfnll1/YsGEDbdq0ASA8PJxOnTo9VPnp6el89dVXvPfee3nSkpKSCAgIYMGCBbz00kvcu3ePHTt2lPr3weNi7G2gyVLKl43ZgP8j2wrH5UFWn5VqFulMCGjErpPX6DJrB7/EXCrp0BRFKeNOnjgJQJ++fTA3N8fa2pqOL3TEx9eH4S8PJz4+nr69+uJs78ys6bMAmDFtBk0aNsHFwYUWzVoQviFcX56Plw+zZ86mtV9ratWoxdCBQ0lLS9OnR0dF07ZVW1wcXLS0uw/SCio3q+yZX86ktV9rnKo7kZ6eXmB5hvTr34+wH8P0j5cvW07ogNAceY4ePUr79u2xs7OjcePGbNyYc6WUyMhInnrqKSpXrkzfvn1zPL+EhAR69eqFo6MjHh4ezJ49u8B4sty6dYsJEyYwb948QkJCsLGxoUKFCvTo0YNp06YZVYaHhwcREREA3Lx5EwcHB1JSHsxF+tZbb/Hll19y8+bNPMceP34cgNDQUP37oHPnzjRp0sSouss7Y5rFHwLxxhYopUwXQnwIXHjoqEoZKWU4EO7n5zf8lTYetK3vwOs/RTFi2UF6P12Lj3o0onLFCiUdpqIoRnr7329zOPpwsdbh29SXz7/8vNB8nvU9MTM349X/e5VeL/WiRcsWVKtWDYDF3y5mz649zFkwh+c7Pq8/xqOuB1u2baGmU03WrVnH8JeHExUbhZOzEwBr16xl7ca1WFW0ovPznflx6Y/83/D/4969e/Tv05+Rr43k1ZGvsil8E68MfoXxb443qlyA1T+tZtW6Vdg72JOZmVlgeYa0aNmCFctXEHcsDs/6nqxdtZYtv23h44kfA3D//n169OjBK6+8wv/+9z927txJz549OXDgAA0aNODevXsEBQUxfvx4xowZw4YNGwgNDeXtt98mMzOTHj160LNnT8LCwjh//jydOnWiQYMGvPjiiwW+Dnv27CEtLY3g4OBCX7P8uLi4cOGC9tVnZ2dHw4YN2bt3r/5KjJ+fH+3bt+fLL79kypQpOY718vLC3NycIUOG0K9fP1q1aqV/HyhGXFmRUn4ipbxoSqG6Yy4/fFilW/2alVk70p8xz3uy9u/zdJm1g33/XC/8QEVRlFyqVKnC/377H0IIxo4aS91adenbqy9XLl/J95jgXsE4uzhjZmZGrz69qOdZT9/HBWDEqBE4uzhTvXp1unbvqm+Y7d+3n/vp9xn92mgqVKhAUEgQT/k9ZXS5WWXXcquFtbV1oeXlp1//foQtC+O3bb/h1cALF1cXfdrevXtJSUnhnXfewdLSkg4dOhAQEEBYWJg+/f79+4wfP54KFSrQu3dvWrRooT2//fu5evUqEyZMwNLSkrp16zJ8+HBWrFhRaExZU8cXdmsrKCgIOzs7/bZ48WJ9WvbGCkBaWlqePo6TJ09mzpw5XL2as2dFlSpV2LlzJ0IIhg8fjqOjI4GBgVy+XG6/Sk1i0g1HIcQuYCHwk5TybvGEVDZYWpjx7xcb8Ly3I2/8FE2/xXv513N1eeMFL6zUTJmKUqoZc8XjcWrg3YCF3ywE4HjccYa/PJy3//023y791mD+5cuWM2/2PM6dPQdASkpKjk6rNWvW1P9tbW3NxQTt9+bFixdxdnHOMUGem5ub0eUCuNZy1f9dWHn56de/H107deXM2TP0G9AvR1pCQgJubm45lkZwd3fXNwISEhLyLKDn7u4OwNmzZ0lISMDOzk6flpGRQdu2bQuNyd7enmvXrhXaF2f9+vX59llxdXXVx3ngwAH++eefPHX7+PgQEBDA1KlTadiwYY60hg0b8t133wFw7NgxBg4cyPjx4/UNtSeZqUOX7wPfAwlCiBlCCO9iiKlMedq9OpvHtqVfi9os+uMfes7dxdGLSSUdlqIoZZRXAy/6D+rP0SNHgbwz7547e46xo8by5cwvOZNwhvjL8TRq3MioUYpOTk5cTLiYI+/5+PMmlZs9noLKK0ht99q413Hn119+JTAoMEeai4sL8fHxZGZmPnjO587h6qo1kpydnblw4UKOOs+d0xpXbm5ueHh4cPPmTf2WnJzM5s2bC42pdevWVKxYkfXr1xeaNz8uLi6cP3+e5ORkBg0axLx586hYMe9AjEmTJrF48eIcV2Fy8/b2ZujQocTExDx0POWJSY0VKWV7oCFag2UwECuE2C6E6CuEKLedNgqbFM7GyoLPQnz57xA/rqXcpefcXSz64xQZmWqIs6IoBTsed5w5s+Zw4bz2xXU+/jyrV66mRUvt1oZjTUfOnD6jz3/79m2EEDg4agvILft+GUdijxhVV8tWLbEwt2DBvAWkp6ezcf1G/W2ehym3oPIKM3fhXMJ/CcfGxibH/meeeQYbGxu++OIL7t+/z/bt2wkPD6dfP+0KTOvWrbGwsGD27Nmkp6ezdu1a/vrrLy2eli2pUqUKn3/+OXfu3CEjI4OYmBijhv9WrVqVyZMnM3r0aNavX8/t27e5f/8+ERER/Oc//zHqObm4uPDPP/8QHBzMwIED9THn5unpSd++fXN0/j127BjTp0/n/HmtsRcfH09YWBitWrUyqu7yzuRJ4aSUcVLKNwBXYChgDiwHzgshpgoh6hZtiCXP2IUMOzasyZbxz9G+gSOfRRwjdPFe4hNvP6YoFUUpi2xtbTmw/wAdnuuAU3UnOrbrSKPGjfjk808AePOtN5k2dRpuNd2YPXM23g29eW3ca3Rs15F6tesRGxtLq9bGfaFZWlqybOUyli9dTm2n2qxdvVZ/ZeNhyi2ovMLUrVeXp57O27/F0tKSjRs3EhERgYODA6NGjeKHH37A29tbn7527Vq+++47qlWrxsqVKwkJCQHA3Nyc8PBwoqKi8PDwwMHBgWHDhmHs7ONvvPEGM2bMYMqUKTg6OuLm5sbcuXMJCgrS5+nRo0eOeVayd8h1dXUlLi6Odu3a8f777xdY14QJE3LMuVK5cmX27dunb6y1atUKHx8fpk+fblTs5Z141AnOhBBPATOA53S7MoF1wGtSynI1ttfPz08eOHCg0HxSSlYfPM+kcO1XycTAxvR6ylUtpKaUGwcTjPv1XFpY3LCgrme5+x1V7lS2UrOElyVxcXE0aNAgz34hxEEppV9R1vVQ0+0LIayFEK8IIf4C9gOOwDjABRgJPAv8WGRRljFCCPr4uRExri2NnKvw71XRjFz2N4mp90o6NEVRFEUpc0xqrAghfIUQc4EEtFFBZ4FOUsrGUso5UspLUsrFwAjAv+jDLVvcqlci7F+teLerN9uOXabzzD/57ZgahqYoiqIopjD1yko0EATMAtyllH2klIYWrTgJ7HnU4MoDczPBq+3qsWF0GxxsLXnluwO8t+4wqXfTSzo0RVEURSkTTG2s9EFrpEwqaKI4KeVRKeXz+aWXNYWNBjJGI5cqbBjjz6vP1SXsr3N0n72Dv8/dKMIoFUVRFKV8MnXo8hopZUZxBVNaGTsaqDBWFua8260hYcNbcT9D0nvBbqb/L477GZmFH6woiqIoTyjTl8wEhBBNgQZAntlupJQ/PGpQ5V2ruvZEjG/LpI1HmPPbSbbHXWVm36Z41lA94RVFURQlN1On27cDNgFZg++zxuJmH/+sGitGqFKxAtNfasoLjWrw7trDdJ+9k3e6ejOkdR3MzNQQZ0VRFEXJYmqflU8Be7Q5VQQQDHRAG6b8D9CySKN7AnTxcWbL68/xbD17JoUfYfCSv7h4605Jh6UoiqIopYapjZUX0Rose3WPz0spt0spBwNb0eZaUUxUo3JFlgxtwSfBPhw8e4MXZ/7JxuiEkg5LURRFUUoFUxsrzsA/uk62aUD2ThZrge5FFdiTRgjBgGfc2TyuLfVq2DI2LJLXwiK5dft+SYemKIqiFKHt27fTsWNHnn/+edatW1fS4ZQJpnawvQRkrb19FmgNbNc99iyimJ5oHg42rHq1NQu2n+KrbSfYfzqRaX2a0La+Y0mHpiiKojyitLQ0pk+fTkREBJaWliUdTplh6pWVnWgNFIClwEdCiEVCiHnANGBLUQb3pLIwN+O1jvVZN8ofGytzBv33LyZujCXt/hM3alxRlFLGx8uH37cZmgtUMcbu3buxtramR48eBAcHc+lSuVpCr9iY2liZBPyi+3saMA/t1k8osBF4rehCU3xrVWXT2LYMfbYO3+0+Q/fZOzh8/uEnplMURTFk0YJFtHu2HQ5VHBgxbESx1BEUEMSUSVPy7N8UvglPd0/S00vfrN5z587Fz88PKysrhg4dmic9MTGR4OBgbGxscHd3Z/ny5YWWefnyZU6ePEl4eDjDhw9n4sSJRR94OWTqpHCnpJQ7dH/fl1K+KaWsJaWsLqXsL6W8XjxhlqyimMH2YVWsYM7EwMYs/b+WpN7NIHj+LuZsO0G6mkhOUZQi4uzszFvvvMWgIYOKrY4BgwawcvlKpJQ59q9YvoKX+r2EhcVDTftVrFxcXPjggw945ZVXDKaPHj0aS0tLLl++zI8//sjIkSOJjY3l0qVLtGnTJs926dIl7Ozs8Pf3x9LSko4dO3LkyJHH/KzKpodadflJU1Qz2D6KtvUd2TL+Obr5OjP91+P0WbSHM9dSSyweRVGKzvdLvicoIIjXX3ud2k61ae7TnGNHjzF/7nwaeTaijmsdNq7fWGz1BwYFEhAYQLXq1Uw67njccXwb+LL6p9UAXEy4yMB+A/Go5YFvA18WzFugzxsQGMCNGzfYvXO3ft+NGzf4ZfMvhA4IBWDUqFGMGjUq3/oyMzOZPHkybm5uuLi4EB4ejqWlJTduFM/SJSEhIQQFBWFvb58nLTU1lTVr1vDxxx9ja2tLmzZtCAwMZOnSpTg5ObFz5848m5OTEy1btuTo0aNIKYmKiqJevXrFEnt5U2hTVgjxmwnlSSllx0eIRylA1UoVmB3anI4Na/Dh+hi6frWDDwIa0r9lbYRQE8kpSlkVExND5N+RjH9jPF/O+pIhA4bQK7AXo8eNJvpoNIsXLWba1GkEBgUWWlaf4D7s3b3XYFqrZ1uxat2qIok5KjKK/n36M332dLp260pmZiZ9e/WlW0A3lvywhAvnL9CzW0/qe9Wn0wudsLa2JrhXMGE/huHf1h+AdavX4dXAC98mvgDMnz+/wDonTpzI9u3b2bt3LzY2NnTt2pWaNWtSrZpxjayAgAB27txpMK1Nmzb8/PPPRj//48ePY25ujpeXl35f06ZN+eOPPwo8zt7enuDgYNq1a4eZmRlLliwxus4nmTHX3czIOUNtA8AJOANcBmoCdYCLQFzRhqcY0rOZKy09qvPWqkO8vy6GbUev8HmvJjhWtirp0BSlTPg84hTHLqUUax3eTra83dW4X80xh2J48603ad+hPQANvBtw9+5dRo3RrjI0atQoT5+Obi90Y/Ovm/OUVVSNkYLs3rWbpd8v5ev/fs1z7Z8D4OCBg1y7eo133n8HAI+6Hgx5ZQhrflpDpxc6AdB/YH/6BPdh2sxpWFtbE/ZjGKEDQ42q8+rVq8ycOZPIyEhcXV0B6Ny5M/v378+Rr3379mzfvt1gGaY0RgqTkpJC7qvtVatWJTk5udBjR48ezejRo4sslidBobeBpJTtpZTP61ZR/gq4D7SSUtaVUraWUtZFGyF0X5euPAbOVa354ZWWfNSjEbtOXqPLrD/59cjlkg5LUZSHEBsTS5duXfSP447F5Xh87NgxvBo8+AWfmJiIvUPeWxOPy5JvlvBMq2f0DRWA+HPxXLx4Ebeabvpt+hfTuXLlij5Pa//W2DvYsyl8E6f/Oc3fB//mpb4vGVXntm3b8PT0xNPzwSwZiYmJ+Pr65njs4OBQBM+wcLa2tiQlJeXYl5SUROXKao234mBqj6aPgQ+llH9l3yml3CeEmAhMATYUUWxKIczMBC/7e+Dv6cD4FVEM/+EAoS1r82FAQypZlr7OaopSWhh7xeNxOHf2HPfu3cOz/oMv4UPRhxj7+lj949jDsfpbJUsWL2Hh/IWkJKcwfsx4Zs2dlaO8kMAQ9uzaY7Cu1v6tWbtx7SPHPGvOLGZOn8k7b73D1GlTAXCt5Yp7HXeiYqMKPDZ0QChhP4Zx4vgJOnTqQI2aNYyq89q1a7i4uOgfZ2RkEBERwaRJkwBYtGgRs2fPJjk5mREjRrBw4cI8ZXTt2pUdO3YYLL9t27ZEREQYFQuAl5cX6enpnDhxgvpbnlPeAAAgAElEQVT16wMQHR1N48aNjS5DMZ6pHWzrA1fzSbuCmhiuRHjVrMy60c/yaru6rNh/ju6zdxIVf7Okw1IUxQgxh2No5NMIMzPt4zgpKYn4c/H4+Pro8xw+fBifJj5ER0Xz559/MvjlwUybOY2YmBhyj1Jcu3EtF69fNLjl11BJT08nLS2NzIxMMjIySEtLK3Aosa2tLWs3rmX3zt189MFHAPi18KNKlSrM/HImd+7cISMjgyOxRzh44GCOY0MHhLL9t+18/+339B/YP0fa0KFDDQ4RBmjYsCG7d+/m5MmTJCUlMXbsWE6dOoWvry+RkZH8/vvvDBs2jDlz5nDo0KE85wUgIiKClJQUg5uhhkrWecnIyMhzXmxsbAgJCWHChAmkpqaya9cuNmzYwKBBxTei6klmamPlNPBqPmmvovVjUUqAlYU573ZtyPJhrbh7P4NeC3YzWw1xVpRS7/ChwzRp0kT/OPZwLHXr1aVSpUqANgLmaOxRmjRpwsb1Gxk4eCCHow/j28QXC3MLrKweva/aF599QQ27Gsz4cgYrw1ZSw64GX3z2RYHH2NnZsX7Ten7d8isfT/wYc3NzVq5ZyaHoQ/h6++Lh6sGYkWNIupXzVol7HXeeafUMt1Nv0y2gW460+Ph4/P39DdbXsWNHQkNDad68OX5+fjRr1oxKlSrh7e3N2rVrefnll4mKiqJZs2ZYWBTNeZkyZQrW1tZMnTqVZcuWYW1tzZQpD+aKmT9/Pnfu3KFGjRqEhoayYMECdWWlmIjcY94LzCxEP7QVlo8Cq3nQwbY34A0MkFKuLIY4SwU/Pz954MCBkg6jULfu3GfChhg2RCXwVG07ZvZthru9TUmHpZQjBxMOFp6pFLG4YUFdz7olHcYje+etd+jUuRMzp81k2sxpfDblM5aGLS3psIrEvXv3aNOyDYcOHaJChQqF5l+4cCGbNm0iPDyc119/nS5dujB16lTmzJnDxIkTWb169WOI+skWFxdHgwYN8uwXQhyUUvoVZV0mNVZ0QXRCm8nWD6iA1rF2P/CRlHJbUQZXXIQQFYEVaCOb7qI1ukZKKf8p6Liy0ljJsiHqAh+sjyEzU/JRYGP6PF1LDXFWioRqrJSMw4cOM2LYCOLj43nW/1lmzZmFk7NTSYdVZCpb5d85de/evTg7O+Pm5sa2bdsYMGAAGzdupFWrVkRHRzNkyBDOnTtH27ZtWbhwIc7Ozo8x8ifT42ysmNwLU0q5FdgqhDADHIBrUsqyeK9hgZRyC4AQYgzwDdChZEMqWj2bueJXpzpvrIziP6sP8dvRK3wW4ks1G7V4lqKURb5NfPl6ydd8+99v+XLmlyUdzmMVGRlJ9+7duX//Pl5eXnz33Xe0atUK0OY3Wbp0KV9//TVz5swp4UiV4lBonxUhRGtD+6WUmVLKK4+roSKEqCWEmCOE2COEuC2EkEKIOvnkdRNCrBZC3BJCJAkh1gohameLPS2roaKzFyj7P7sMcLWzZvnwVrzb1Zttxy7z4qw/+fN4fn2kFUUp7U4cP5FjGPOTYuTIkVy/fp2kpCQOHDhAt245+7vExcXh7e1dQtEpxc2YKys7hRBX0BYqXAdsk1LeL96wDPIEXgIOAjuAzoYyCSEqAb+h3d4Zgjah3RTgdyFEEymloTnqX6McD7k2NxO82q4ebeprQ5wHL/mLoc/W4Z2u3lSsYF7S4SmKYoKgkKCSDqFU6t27d0mHoBQjY0YDuQKTAXdgPXBNCLFSCNFPCFGlWKPL6U8pZU0pZTegoCkah6NdJQmSUq6XUm4AAtHizzOSSQjxLuAFvFsMMZcqjV2qEv5aG17211Zx7jFnJ7EJahVnRVEUpXQzZgbbS1LKBVLKLoAjMALtasVC4IoQYosQYoQQwqXAgh6RCbebAoG9UsqT2Y49DewCembPKIT4N9AL6CqlvF1UsZZmFSuY81GPxvzwSktu3blP0LxdLPrjFJmZpnW0VhRFUZTHxaR5VqSUyVLKMCllP7SGSzDa3CsfAvFCiH1CiHeKIU5TNAZiDOyPBRplPRBCvAGEAi9IKfOdQU0I8S8hxAEhxIGrV8tPX4/nvLRVnDt61+SziGP0/2YvF27eKemwFEVRFCUPUyeF05NS3pdSRkgpR0gpXYE2wHZgcFEF95CqA4bWC08EqoHWWReYDtih9WWJEkIYHJMspfxaSuknpfRzdHQsrphLRDUbSxYMfIovejfh8PlbdJn1JxuiLpR0WIqiKIqSgzGjgZYIIZ4rLJ+Uco+U8m0pZaPC8j4Ghu5p6CcYkVKel1IKKWU9KWUz3ZbvmHAhRA8hxNeGpm8u64QQvOTnxuZxbalfw5ZxK6IYtyKSW3dKog+1oiiKouRlzJWVvmhXH04LISYJIUr7+j830K6u5FYNw1dcCiWlDJdS/iv3cuDlibu9DT+92po3XvDi50MX6TrrT/acul7SYSmKoiiKUY2VmsAwtHV/PgDihBA7hRDDhRCl8ds7Fq3fSm6NgCOPOZYyxcLcjLEd67Nm5LNYWpjR/5u9fBZxlLvpGSUdmqIoivIEM2Y0UIqU8lsp5fNAHbTOtNWBRcBFIUSYEKKrbkbb0mAj0EoIoZ/kTTd5nL8uzWTl+TaQIc3c7Ng0ti39WtRm0R//EDxvNycuJ5d0WIqiKMoTytTRQPFSyk91/VJaAUuAjsDPwAUhRLHO/yyE6C2E6A08rdvVVbevXbZsi9GuAm0QQvQUQgSiTfgWj9bAMtmTcBsoNxsrCz4L8WXxYD8uJaURMGcn3+06jalrSSmKoijKo3qU0UB/SSnHoE0aNxOoAbxeVIHlY5VuG6F7PF/3eFK2uFLR1vg5DixFWyX6NNBBSplSzPGVOy80qskv49vybD17JoYfYci3+7mSlFbSYSmKopQb27dvp2PHjjz//POsW7eupMMplUxeyDCLrqPtYGAg2uywyRQ8s+wjk1IatWSwlPIc2mRvRUII0QPo4elZ2vsWF48alSuyZGgLlu07xyebjvDirD/5LKQJXXzKz2qviqIoJSEtLY3p06cTERGBpaVaZDY/Jl1ZEUJUF0KMFELsAeKA99CuYAwEnKSUw4shxhL3JN4Gyk0IwaBW7mwa25Za1SoxYtlB3l17iNv30ks6NEVRlDJr9+7dWFtb06NHD4KDg7l06VJJh1QqGTPPSgUhRLAQYh2QAMwDKqOtpVNbStlFN6utujfwBKjnaMuakc8ysn09VuyPJ2D2Tg6dz3cCYEVRyoBFCxbR7tl2OFRxYMSwEQXm9fHywb6yPdev5ZzawL+lP1UqVuHsmbOF1hcUEMSUSVPy7N+wYQNOTk6kp5e+H0Fz587Fz88PKysrhg4dmic9MTGR4OBgbGxscHd3Z/ny5UaVe/nyZU6ePEl4eDjDhw9n4sSJRRt4OWHMlZVLwGq0GWoXAy2llD5Syi+klAnFGp1SKllamPF2F2+WD2vFnfsZhMzfzYLtp8hQ6wspSpnk7OzMW++8xaAhg4zK717HnVU/PbjrHxsTy507xi/XMWDQAFYuX5mnw/7SpUsZMGAAFhYP3UOh2Li4uPDBBx/wyiuvGEwfPXo0lpaWXL58mR9//JGRI0cSGxsLwKVLl2jTpk2e7dKlS9jZ2eHv74+lpSUdO3bkyBE1w4YhxjRWdgC9ARcp5WtSSoPT0pdnT9rQZWO1rmdPxLi2dG5ck89/OcaAb/aSoNYXUhSTfb/ke4ICgnj9tdep7VSb5j7NOXb0GPPnzqeRZyPquNZh4/qHmnnBKIFBgQQEBlCtejWj8vfr34+wH8P0j5cvW07ogNA8+S4mXGRgv4F41PLAt4EvC+YtACAgMIAbN26we+dufd4bN27w888/M3iwtmLLqFGjGDVqVL4xZGZmMnnyZNzc3HBxcSE8PBxLS0tu3HiouT8LFRISQlBQEPb29nnSUlNTWbNmDR9//DG2tra0adOGwMBAli5dCoCTkxM7d+7Mszk5OdGyZUuOHj2KlJKoqCjq1atXLPGXdcbMsxIkpVwnpXxi519XfVbyZ1fJknn9n2Ja7yYcOn+Lrl/tYNOhiyUdlqKUKTExMUT+HUnP4J6cvnCaxj6N6RWojRGIPhrNf979D9OmTjOqrD7BfXCr6WZw6xPcp0jibdGyBclJycQdiyMjI4O1q9bSN7RvjjyZmZn07dUXH18f4v6JY2PERhbMWcDWX7dibW1NcK/gHA2edavX4e3tTdOmTQGYP38+8+fPzzeGiRMnsnXrVvbu3cuRI0f49NNPqVmzJtWqGdfgCggIwM7OzuAWEBBg0vk4fvw45ubmeHl56fc1bdpUf2WlIPb29gQHB9OuXTvefvttPvzwQ5PqflKYfK1NCNEcbWK459AWAmwppfxbCPEp8KeU8pcijlEp5YQQ9PFzo0Wd6oxbGcXo5X+zPa4WHwU2xtaq9F3OVRSrXz/C7ErhXySPIrNGY+6+MKnwjEDMoRjefOtN2ndoD0AD7wbcvXuXUWO0KwuNGjXK04+j2wvd2Pzr5jxlrVpXrIMy9fr170fYsjD8n/PHq4EXLq4uOdIPHjjItavXeOf9dwDwqOvBkFeGsOanNXR6oRP9B/anT3Afps2chrW1NWE/hjFkyBCj6r569SozZ84kMjISV1dXADp37sz+/fvz5G3fvj3bt2/Ps//nn3828RnnLyUlhdw/ZqtWrUpysnGTaY4ePZrRo0cXWTzlkUnfJEKINsBW4B9gOTAmW3Im2vwnqrHyhKrjYMPqEa2Zve0E834/yV9nEpnVtxnNaxv3S0dRnlSxMbF8Ne8r/eO4Y3F06dZF//jYsWN4NXjwqz0xMRF7h7y3Ix6nfv370bVTV86cPUO/Af3ypMefi+fixYu41XTT78vIyKC1f2sAWvu3xt7Bnk3hm3ja72n+Pvg3G9ZvMKrubdu24enpSfbpJBITE/H19c2RLzExEQcHh4d5eiaxtbUlKSkpx76kpCQqV65c7HU/KUz92TsV2AIEAebkbKz8jTbvSrnzpM+zYooK5ma82bkBbes78vrKKHov3MPrneozsr0n5mZGTZOjKMXO2Csej8O5s+e4d+8envUffL4cij7E2NfH6h/HHo7Ft4n2Rbxk8RIWzl9ISnIK48eMZ9bcWTnKCwkMYc+uPQbrau3fmrUb1xZJ3LXda+Nex51ff/mVeQvn5Ul3reWKex13omKj8i0jdEAoYT+GceL4CTp06kDNmjWNqvvatWu4uDy4kpORkUFERASTJj14XRctWsTs2bNJTk5mxIgRLFy4MEcZXbt2ZceOHQbLb9u2LREREUbFAuDl5UV6ejonTpygfv36AERHR9O4saFl6pSHYeoMtk8BC6TWhTv30I9rgGORRFXKqD4rpmvpUZ3N49rS3deZL/93nNCv93L+xu2SDktRSp2YwzE08mmEmZn2cZyUlET8uXh8fH30eQ4fPoxPEx+io6L5888/GfzyYKbNnEZMTAy5O/6v3biWi9cvGtzya6ikp6eTlpZGZkYmGRkZpKWlGTV8eO7CuYT/Eo6NjU2eNL8WflSpUoWZX87kzp07ZGRkcCT2CAcPHNTnCR0QyvbftvP9t9/Tf2D/HMcPHTrU4BBhgIYNG7J7925OnjxJUlISY8eO5dSpU/orK5GRkfz+++8MGzaMOXPmcOjQoTznKSIigpSUFIOboYZK1jnKyMjIc45sbGwICQlhwoQJpKamsmvXLjZs2MCgQcaNrlIKZ2pjJQ2olE+aM6CGyyh6Va0r8FW/Zszs25QjF5Po+tUONkar0e6Kkt3hQ4dp0qSJ/nHs4Vjq1qtLpUraR21mZiZHY4/SpEkTNq7fyMDBAzkcfRjfJr5YmFtgZWX1yDF88dkX1LCrwYwvZ7AybCU17GrwxWdfFHpc3Xp1eerppwymmZubs3LNSg5FH8LX2xcPVw/GjBxD0q0Ht0vc67jzTKtnuJ16m24B3XIcHx8fj7+/v8GyO3bsSGhoKM2bN8fPz49mzZpRqVIlvL29AVi7di0vv/wyUVFRNGvWDAuLRz9PU6ZMwdramqlTp7Js2TKsra2ZMuXBXDHz58/nzp071KhRg9DQUBYsWKCurBQhU28D7QTGCyGy31jMusLyf8BvRRKVUm4IIQhuXgs/9+qMWxHJ2LBIth+7wqSejalcsUJJh6coJe7t997O8bi1f2sOHnpw9cHMzIzLNy4D2hBZgPPx50lJScGxhiMVK1Z85Bje+/A93vvwPaPyxhyPMbjfwsKCpLSc/TacXZz5dum3BZZnqJPwvXv3SEhIyPfKCuQcLbRw4UI6dOign64+JUVbBu7cuXMkJydTo0aNRz5PEydOLHDCturVq7N+/fpHqkPJn6mNlQ+BXUA02kRxEhgihJiBthJyi6INTykv3KpX4qdXWzP395PM3naC/We1zrdPu1cv6dAUpcwYMGgAI4aNID4+nskfTWbWnFmFH1QGWVpacvTo0XzT9+7di7OzM25ubmzbto0JEyawceODeWiGDh3KkCFDOHfuHO+//36e/ipK2WPSbSApZTTakOXLwPuA4EEn23ZSyriiDa90UJPCFQ0LczPGd/Ji1QhtNMBLi/Yya+tx0jMySzgyRSkbfJv48vWSr3mp30usWL0CJ+cnczHRyMhInnrqKezs7Hj33Xf57rvvaNWqlT69adOm+tlwN2zYgLOzcwlGqxQFU/usIKX8W0rZEW19oFpAFSnl81LKSCGEbZFHWAqoDrZF62n36mwe25aeTV2YtfUELy3aQ3yi6nyrKMY4cfxEjmHMT6KRI0dy/fp1kpKSOHDgAN26dcuTJy4uTt+HRSn7RO61GQrMLMRsKeXYfNJsgS1SSsM9osoBPz8/eeDAE7faQLHaEHWBD9bHICV8HNSYoGauCKGGOJd2BxMOFp6pFLG4YUFdz7olHYZSiMpWal6SsiQuLo4GDRrk2S+EOCil9CvKuky9svKyECJPLyxdQ+UXwC3vIYqSv57NXIkY15ZGzlV4fWU041ZEcevOE7uyg6IoimKAqY2VPsBHQoiXs3YIISoBEYAH8HwRxqY8IWpVq0TYv1rx785ebDp8kW5f7eCv04klHZaiKIpSSpjawfYXYDiwUAgRIISwRmuo1APaSylPFUOMyhPA3EwwpkN91ox8FgtzQb+v9zD9f3HcV51vFUVRnngP08H2B7QhzD8BfwANgA5SyhNFHFupoUYDPT7N3OzYNLYtvZ+uxZzfTtJ74R7OXEst6bAURVGUElRoY0UIYZZ7A74EvgHqAC8Ax7OllTtqNNDjZWtlwRe9mzJ/wFOcuZZK99k7WHPwPKZ0BlcURVHKD2MmhUsn7zpAWQSQfZUqaWSZilKobr7ONHOz4/WVUby5KpodJ67ycZCPmvlWURTlCWNMw2Iy+TdWFKVYudhZs3x4K+b/fpJZ207w97mbzA5tTjM3u5IOTVEURXlMjGmsbJFSGl5vXFEeA3MzwWsd69O6nj3jVkTRe8Fu/v1iA/7Vti5mZmpOFkVRlPLOmD4mO4UQF4UQi4QQXYQQ6hq8UiL86mgz377Y2ImpEccYvOQvriSllXRYiqIoSjEzprHiinYryB1YD1wTQqwUQvQTQlQp1ugUJZeqlSowt39zpob4cuBsIl2/2sHvcVdKOixFURSlGBXaWJFSXpJSLpBSdgEcgRFofVgWAleEEFuEECOEEC7FHKuiACCEoF/L2vz8WhscK1vx8rf7+fjnI9xNzyjp0BRFUcq0PXv20Lp1a9q1a0doaCj375eOGcVNnRQuWUoZJqXsh9ZwCQZOo827Ei+E2CeEeKcY4lSUPDxrVGb9aH+GtHbnvztPEzJ/N/9cTSnpsBRFUcosd3d3fvvtN/744w/q1q3Lhg0bSjok4CEmhcsipbwvpYyQUo6QUroCbYDtwJCiCq60UJPClV4VK5gzqacPiwf7ceHmHQLm7GTVgXg1J4uiKMpDcHFxwdraGgALCwvMzErH9GlFFoWUco+U8m0pZcOiKrO0UJPClX4vNKrJL+Oeo0mtqry1+hDjVkSRnFY6Ll8qSmm3aMEi2j3bDocqDowYNiJHWmJiIv1f6o9TdSca12/MTyt+MindkNU/reb5ts/jVN2Jum51eb7t8yxetLjQHxkvvvgiEyZMyLN/w4YNODk5kZ6ebsSzffzat29PxYoVsbW1xdbW1uBKxQ9j7ty5+Pn5YWVlxdChQ/OkJyYmEhwcjI2NDe7u7ixfvtzosk+fPk1ERAQBAQFFEuujMmkCNyHEbwUkZwK3gIPAf6WUlx8lMEUxlVPVivw4rBULtp9k5tYTRMbfYHa/5jSvXa2kQ1OUUs3Z2Zm33nmLbb9u486dOznS3hz3JpaWlpw8d5LD0YfpE9wH3ya+NGzU0Kj03ObMmsOsGbOYPms6HV/oiK2tLYeiDzF75mzuvXoPKyurfOMcOnQo7733HpMmTUKIB9MWLF26lAEDBmBhUXrnJJ07dy7Dhg3LN11KSVRUFM2bN8+xPzo6Gh8fH8zNzfMc4+LiwgcffMCWLVvyvG4Ao0ePxtLSksuXLxMVFUX37t1p2rQp9vb29O7dO0/+1atX4+TkRFJSEkOGDGHp0qVYWlo+xLMteqZeWRFoawG1RxsdVFH3b3ugIdrKyx8CMUKIRkUWpaIYKWtBxJ9ebUVmJvRZuIf520+SmaluCyml1/dLvicoIIjXX3ud2k61ae7TnGNHjzF/7nwaeTaijmsdNq7fWGz1BwYFEhAYQLXqORv2qampbFy/kfc/eh9bW1ta+7ema/eurFi+wqj03G7dusUnkz9hxlczCAoJonLlygghaNqsKf/9/r/6hkpCQgK9evXC0dERDw8PZs+eDUBQUBCJiYns2LFDX+aNGzf4+eefGTx4MACjRo1i1KhR+T7XzMxMJk+ejJubGy4uLoSHh2NpacmNGzce/gQWgTNnztC5c2d++eUX/b4dO3bQsWNHjh49avCYkJAQgoKCsLe3z5OWmprKmjVr+Pjjj7G1taVNmzYEBgaydOlSnJyc2LlzZ54t6+pUaGgoEydOLLIrQEXB1MbKDCANeFpKWU9K+ayUsh7QQrd/ElAfuAp8UqSRKooJnnavzuZx2pwsX/wSx6Al+9ScLEqpFRMTQ+TfkfQM7snpC6dp7NOYXoG9AIg+Gs1/3v0P06ZOM6qsPsF9cKvpZnDrE9zHpLhOnjiJubk59evX1+/zbeLL0SNHjUrP7a99f3H37l269+ieb52ZmZn06NGDpk2bcuHCBbZt28asWbPYsmUL1tbWvPTSS/zwww/6/D/99BPe3t40bdoUgPnz5zN//vx8y584cSJbt25l7969HDlyhE8//ZSaNWtSrZpxV2ADAgKws7MzuBV0y+Tdd9/FwcEBf39/tm/fnifdw8ODNWvWMGDAAH7//Xf++usvQkJCWLZsGT4+PkbFlt3x48cxNzfHy8tLv69p06bExsYWeFxYWBj79u1j8uTJtG/fnpUrV5pcd3Ew9ZrZFGCilDIy+04p5UEhxCRgipTSVwgxDW2xw3IlOjoaDw8PKlasmGezsrIyuL8o0gxd/lMKV9Vam5PluQMOfLQxli5f7WB6n6Y8712jpENTStisv2dx4kbxLhRfv1p9xj813qi8MYdiePOtN2nfoT0ADbwbcPfuXUaN0a4QNGrUKE9/jG4vdGPzr5vzlLVq3apHiju71JRUqlTNOZ1WlSpVSElOMSo9t+vXrmPvYJ/jdk2n9p2IOxrH3bt32bJlC1ZWVly9elXfN6Vu3boMHz6cFStW8OKLLzJkyBC6d+/OnDlzsLa25ocffmDIEOPGdVy9epWZM2cSGRmJq6srAJ07d2b//v158rZv395go+Lnn382qq7sPv/8cxo1aoSlpSUrVqygR48eREVFUa9evRz5nnvuOZYvX07v3r0xNzfn66+/pkuXLibXB5CSkkLufpZVq1YlOTm5wOMGDRrEoEGDHqrO4mRqY8ULuJZP2lXAU/f3KcDmYYMqrezs7HjuuedIS0vLsSUmJpKWlsbdu3fzpN29e/eR67WwsMjTkLG2tsbW1hYbGxuTttzH2NraUrFixRz3f8sTIQR9W9TmafdqjFkeycvf7ecVfw/e7toAKwvVCFRKh9iYWL6a95X+cdyxOLp0e/AldezYMbwaPPiFnJiYiL1D3kv/Rc3G1obkpJxfbsnJydhWtjUqPbfq9tW5fu066enp+gbL1u1bAfCu501mZiZnz54lISEBO7sH639lZGTQtm1bANq0aYOjoyMbNmygZcuW7N+/n7Vr1xr1fLZt24anpyeenp76fYmJifj6+ubIl5iYiIODg1FlGuOZZ57R/z1kyBDCwsLYvHkzr732Wp68tWvXxsLCgszMTOrUqfPQddra2pKUlJRjX1JSEpUrV37oMkuSqY2VM8AwIMJA2r906QAOwPWHjqqUcnd35/vvvzfpmMzMTO7du5dvYyZ3w8aYtDt37nD79m1SU1NJTU3lypUr+r+ztszMTKNjrFChAlWqVHnkzcbGptQ2erLmZJkacYwlu06z/0wic/s3x92+3LWpFSMYe8XjcTh39hz37t3Ds/6DL9BD0YcY+/pY/ePYw7H4NtG+UJcsXsLC+QtJSU5h/JjxzJo7K0d5IYEh7NlleDm31v6tWbvRuC92AM/6nqSnp3Py5En9F/zhQ4f1nWcLS8+t5TMtsbKyYlP4JnoG9zSYx83NDQ8PD06cyP/K1+DBg/nhhx+Ii4ujc+fO1KxZ06jnc+3aNVxcHsxfmpGRQUREBJMmTdLvW7RoEbNnzyY5OZkRI0awcOHCHGV07do1R5+Z7Nq2bUtEhKGvx5yEEAZHPp06dYoXXniBzz//nMqVK9OtWze2bt1K48aNjXp+2Xl5eZGens6JEyf0t+mio6MfqqzSwNTGymRgmRDiELAGuALUAHoBPkB/Xb5OwL6iCrIsMzMz018NeVyklNy9ezdPAyZrS0lJyfF3UkhXm1sAACAASURBVFJSnu3SpUscP36cW7dukZSURFpa4f09zMzMsLOzo1q1anm26tWrG9yflZbV0a44VaxgzsTAxrSuZ89bq6LpPnsnn4b4EthUTb6slJyYwzE08mmkn88iKSmJ+HPx+Pg+6Kdw+PBhugd2Jzoqmj///JPBLw/Gw8ODmdNncuvWrRyX+01pjGRJT08nPT2dzIxMMjIySEtLw8LCAhsbG3oE9eCTSZ8wd+FcDkcfZvPPm/l1+68AhabnZmdnxzvvv8Mb495ASkmnzp2oVKkSMYdjuJ16G4CWLVtSpUoVPv/8c8aOHYulpSVHjx7lzp07tGjRAtAaK1OmTOHQoUPMnDkzRx1ZQ3i/++67PPU3bNiQDz/8kJMnT1KjRg3effddTp06pb+yEhkZye+//86wYcOoW7cun3/+eZ7za0xjJLubN2+yb98+2rVrh4WFBStXruTPP/9k1qycjcyEhAQ6duzI+++/r38OaWlpdO7cmR07dlC3bt08ZWe9bhkZGTlet6zXLiQkhAkTJvDNN98QFRXFhg0b2L17t0nxlxYmNVaklGFCiGtoHWnfAyoA94EDQGcp5VZd1jcANfd5CRFC6BtIhnqJP4x79+6RnJxssGGT1aC5desWN2/e5MaNG/rt7NmzJCYmcuPGDTIy8n9LmJub6xs69vb2ODg44OjoiIODQ75/V61a9aEmLHqxsRONXaowNiySsWGR7Dl1jQkB/8/encdFWbV/HP8cNtlXEQVBBBQVd820ckvbNM20LO1pfVosbbWyXfulZWamldpiT7ttppVptpiamVYuaSKCICDusu8wzFy/P27EDWQbZBjP+/W6XwMz9z1zaQtfzn3OdWJwc9G3hbRz798d/9K1a9eK72P/jSUiMgJ3d3fAGJ2Ni42ja9euvPfue/znlv/w1edfMfKakTg5Op11qW9NzXppFjNnzKz4/ovPvuCJp5/gqWefYs68OUy8ZyKRoZH4+/sz5/U5p4ycVPf66R6a/BCtglsxd85c7vnvPbh7uBPeNpznpz/PRRddhKOjI8uXL2fy5Mm0bduWkpISoqOjmT59esV7hIeHc9FFF7F9+3ZGjhx5yvunpaVx4403VvrZQ4YMYdy4cfTo0cNYrv3YY7i7u9OhQwcAli5dyu23387ixYsZPXo0Tk71//s1mUw888wz7N69G0dHRzp06MA333xzxkqbgIAAXn31VcaMGVPx3E033YSnpyctWlQ+z2769OmnjAp98sknTJ06lWnTpgHGZOM77riDFi1aEBAQwMKFC5vsyIqqa6dPpZQDxu2edBGp+T0HG6GUegqj2247YLSIfFPdNb1795bNmzc3eG32SETIz88/JcicfBwPNFlZWWRkZJCenk56ejrHjh2rclTH0dGxIsBUFWpatmxJq1ataNmyJb6+vqeM3pjMFub8nMDCtUlEB3nx5vgetAtqmvdzz7UtB7c0dgm14pTlRETUmb+ZNjVPPPYEQy8fymuvvMYrr73CS9Nf4uPPPm7ssqzGq1n9/vsrLS2lW7du7NixA2dn52rPf+utt1ixYgXLly8H4OGHH+bKK69k5syZvPHGG0ybNo0lS5bUqyZ7Fh8fX+nyZqXUFhHpbc3PqnNYqfTNlPIUkSaxOYtS6kKMycLvAXN1WLFNIkJhYWFFcDk5xFT1dUZGRqX3g11cXGjZsuUZR5FvW1ZmNMckiof6B3PrgOiK32q1yumw0jj+3fEvE+6cQFpaGhddfBFz35hLy1YtG7ssq6lvWKnOpk2baNWqFaGhoaxevZqbbrqJ7777jr59+wLGnI5bb72Vffv20b9/f9566y1atWrVoDU1ZecyrNS2g+3rIvJAFa95Aj8CF1ujsErevzUwBegNdAPcgLYiklLJuaHAa8BlGI3sfgEeEpF9x88RkT/Lz22IcjUrUUpVrFxq06ZNja4xm81kZ2dz9OhRDh8+XOmRkpLCpk2bOHbsGCKCo4cfAVc/ystmeHb+p5j//JTWQc0JDQ0942jdujWhoaEV+2do2rnSpWsX3vnfO7z/3vvMfs3uukM0uG3btjF8+HBMJhPt27fngw8+qAgqYPQh+fjjj3nnnXd44403GrFS7XS1nWB7u1LqsIi8ePKTSikPYBUQarXKzhQFjMVo578euLyyk5RS7sCvQAnGbR7B6A+zRinVVUQKGrBGzQY4OjoSEBBAQEAAHTuefauqsrKyilBz8NAhluzKZ13MINw79qFVyiqOJG1n8+bNHDt27IxrAwICKgJMeHg4kZGRFUd4eLgOM1qD2JOw55RlzFrN3Xvvvdx7771nPSc+Pr5iDotmO2obVq4HvlVKHRKR96EiHKzCaLU/wMr1new3EQkq/8w7qSKsAHcBEUC0iCSWn78D2APcg9GFV9MAo4dNcHAwwcHB9ASuHg5/7s3ggc+3kdLuep55+Dlu7tuGkpIS9u/fT1paGmlpaad8nZqaypo1a8jPP/UOaEhICJGRkURERJwSZCIjI6028Vk7/4waPaqxS7Brle2ZozW+2q4GWqWUugt4Vyl1DFiN0XMlEhgoIkkNUOPxz67pJN6RwKbjQaX82mSl1AbgGnRY0apxYUQAPzw4gEe/2s5z38ayITGdWWO6ndFM6mQiQnp6OklJSSQlJbF3796Kr3/88UcOHTp0yvkBAQG0b9+e6OhooqOjK76OioqyyuoOTdM0e1LrLSpF5COlVEvgS2AnEAYMEpGG7V1dczHAt5U8H4sxMlQrSqm7MRreERYWVr/KtCbD38OFRbf05n8bkpn5w26Gvb6eN8b3oGcVOzgrpQgMDCQwMPCUe+DHFRYWkpycTFJSEomJicTHx5OQkMCPP/54Sj8IBwcHwsPDK8JLp06diImJISYm5pSOnpqmaeeTasNK+RLl080GWgM3AkOAhOPn2cAyZn+gsu0zM4GKnzRKqWeACUAg0Fkp9SbQW0QOn3yRiLwDvAPGaqCGKlqzPQ4Oijv7R9A73J9Ji7cy9q2NPHpFNHf3j8DBoXYTs93d3StCx+lyc3PZs2cP8fHxFSEmPj6e3377jcLCworzQkJCKt6jc+fOxMTE0KlTpybbPlvTNK2majKyUoYxSbUyCvjnpO+lhu/Z0Cqr95SfLiIyHWPibbWUUiOAEVXdAtDsW/dQX1Y80J8nl+5g5g+72ZiUwZyx3QjwtM7tGm9vb3r16kWvXr1Oed5isbBv3z5iY2PZuXMnsbGxxMbGsnDhwlN6z4SFhVWEly5dutCzZ086dOigN8DUNM1u1CRY/B9VhxVblIUxunI6PyofcamWiCwHlvfu3fuuOldlLoPibCjKhqIs4+uSPDAVQmkhmArKHwuhtODEo7kULGVgNoHFbHxtMZU/mkEsoBzBwQGUg/G1cgAHR+NrR2dwdgdnN+PR5aSvnd2hmSe4+Z10+BuPLh6gl3VX8HFzZv74nnzy5z5e+H4Xw15fz9wbetAvsuEmyh6/JRQeHs7w4cMrnjebzSQnJ1eEl+NB5pdffqG0tBQwRnJ69OhREYJ69+5NdHT0eR1gRES3KtA0KznX/z1V2xROKdVPRCrfFauRlK8GepdK+qwopX4FXETkktOeX4vx5x1Y1889pSmcqRgKjkL+Mcg/Uv51+VFw1AgkRVlQlGM8lp59W+4Kx0OEi4dxODUDB6czD0dnI5CgjMAilvLwYj7pa4sRdkyFYCoyjtIC49FczW7QDs7g0Rw8g8CrpXF4tgSvIPBqBd4h4BsGbuffPIpdB3OZtHgrKRkF3H9pOx4Y0g7HWt4WaghlZWXEx8ezZcuWimPbtm0Vt5I8PT3p06cP/fr146KLLqJv3774+1eW66vX5JrC5TgRFhqGk7MtDPxqVWnopnCa9ZSWlpKWlkZkZOQZrzVKB1ullBljw8LvgGXAahExWbOI2qomrDyEMaemvYjsLX8uHGPp8hMi8modPm8EMKJba/e7/nm8oxFISnIqP9nVFzwCwT3A+EHu6ls+YlH+6Op74vlmXuUjHR7Go5ObMUJyLljMRogpySsf7ck8KWBlQWEmFKRD/mHIO2I8FqRzxiCbqw/4tgG/NuWP4RAQCc2jwTvYbkdnCkrKePbbnSzdeoB+EQHMG9edFl7nbrPKmjKbzezevZstW7bw559/snHjRnbs2FGxT1N0dDT9+vWrODp16lSj0ZemFlYoBB8nHwKDAuu0n5R2buiw0jRYLBYOHTpEs2bNaN68+RmvN1ZYaQlci7HsdxBGs7VVGMFlpYjkWrOgamo5vgB+CMbk2PuAY8AxEVlXfo4HsB0oAp7B+On6AuAFdK3PdgC9w31k86xRxmiDZ6Dx6NECPMsPj0BjJMRemU1QcAzyDkHOfshKhezUE4/Z+6DspH18XLygeTsIjIbm7Y3HoBgj1NhJiPlqcxrPfrsTL1dn5t3YnYsiz/wP19YUFBTw999/s3HjRjZu3Mgff/xBRkYGAD4+PvTv359BgwYxaNAgunfvXml4aXJhRcAh3wEHs0PTuql9nnF1sr3Ar51JKYW7uzshISGV3gpq9L2BlFJewNUYweVKwBVYhxFcvhORg9YsrpLPr6rYdSIy6KTzwji13f5qjHb7KfX5fL03UDUsFuOWWEYipMfDsYQTj3kn/avh6gMtu0Krbicem7crv63V9MQfzuPeT7eQkl7Aw0PbM3FwVK1XCzUmESExMZGNGzfy+++/s27dOhISEgBj8u+AAQPOCC9NLqxoTUKv4F7Vn6TZvEYPK6cV4wwMxQguI4CWwGZgmYjMPNu1Tc1Jq4Hu2rPHVtrJNDHFuZCeAId3wKEdcGg7HN11YiTGyQ1CekJoHwi9EFr3AY+m0+W1oKSMp5b9y7f/HGRA+0Bes+JqocZw6NAh1q1bx9q1a1m7di3x8fHAifASfUE0fQf2JTwqXE9a1axGhxX7YFNh5Yw3UqofMAoYISKdrPKmNkaPrFiZuexEgDn4D+z/ywgxljLj9YAoI7iE9oG2A8CvrU3fPhIRPvsrjWnLY/F3d+GN8T24ILxuE1htzcnh5ddff+V4aA8KDqLvwL70HdCXPv374OPn08iVak2ZDiv2odHDSnnjNwcRKTvpuSuAzsCvIrLNmsXZGh1WzgFTERzcBml/QtpfxmOhMZ8CnzCIGAgRg4zw4tmiMSutUuzBHCZ+upW0rCIevyKau+rQRM7WLf9zOZvWbWLTuk389ftf5Ofmo5SiU/dO9B3Ql4FXDKRj14561EWrFR1W7IMthJUvgBIRuaX8+wnAgvKXTcBwEfnFmgXaAn0bqBGJQPoeSF4He9dC8voTK7GCOhvBpd3l0OYiYzm3jcgtNvHE1ztY+e9hhnRowatju+Hr7tLYZVnNyXNWysrK2PXPLjb9tolNazexc9tOLBYLLUNaMviqwQy5eghde3XVq3C0aumwYh9sIaykAlNE5PPy75MwJq9OxmhJ31JEBluzQFuiR1ZsgMUMh/4xgsvedbBvk9Ezppk3RA2B9ldBu8vAvfFvv4gIH21MZfqKXbTwcuXN8T3oUcXeQk3N2SbY5mTl8NtPv/Hryl/Z9NsmTKUmmgc159JhlzL06qF0u6Dbed2cTquaDiv2wRbCShFwuYisV0pFAQlAdxHZoZS6HFgsIra/drOOdFixQSX5RnBJWAUJPxoN+ZSDMdcl+iroOBL82zZqidvTspm4eCtHcot58qqO3H5x05+UWtPVQPl5+fy++ndWf7+aP9b8QUlxCQEtArj0qksZMnwIPfr20MFFq6DDin2whbByBLhHRL4pb8z2ooi0KH9tCPCtiHhas0BbosOKjbNYjPkuCT8Y4eXwv8bzwT0g5lroNMpoXtcIcgpNPLpkOz/vOsKVMS15+bqu+LjZzm2r2qrL0uXCgkI2rN7AL9//wu+rf68ILlddexXDrxtOu07tGqBSrSnRYcU+2EJYWQa0wtgAcCawWURuK3/tHuBhEelgzQJtgZ6z0kRlpUDsN7DrGyPEAIT0MkJLzLXgG3pOyxER3vs9mZk/7CbY140FN/Wkc0jTXD1T3z4rRYVFbPh1A6uWreL31b9TZiqjfaf2DBszjCuvvZLmQXY7QKudhQ4r9sEWwko7YAUQBewFhh5vtFa+J0+qiNxuzQJtiR5ZacIyk43QErvMWB4NEN4fuo2DTiONrQ/OkS2pWUxavJWM/FKeHdGJ/1wY1uRuC1mzKVx2ZjY/ffcTK5asIHZbLA4ODvQd2Jdh1w1j0BWDcHXTXU3PFzqs2IdGDysnFRIgIhmnPdcFOCwix6xVnK3RYcVOZCTBzq/hn8WQlWxsHNlxhBFc2g44J510MwtKeeTLf1gbf4yru7Zi5piueDZrOpvsNVQH25TEFFZ+vZKVX6/k8IHDeHh6MHTEUK675To6du3YIJ+p2Q4dVuyDzYSV85UOK3ZGxOjjsv0z2LnMWBLtHQLdboSetzb4/BaLRVi4LolXf4onvLkHb/2nF+2DmsZGbg3dbt9isbB101ZWLFnBL8t/oaiwiJgeMVx/y/UMHTFUj7bYKR1W7EOjhxWl1F7gWhHZXslrnTH2B4qwYn02RYcVO2YqgvgfjOCS+IsRZKKGQu87jD4ujg036rFpbwaTFm+joKSMmWO6cE33kAb7LGs5l3sD5efms+LrFSz5cAnJe5Lx8fVh5I0jGXPLGFq3aX3O6tAang4r9sEWwooF6Csif1XyWm/gTxGxu3WIeoLteSZnP2z9CLZ8CPmHjdGWnrdCz5vBO7hBPvJIbjGTFm/l75QsbrsonKeGdcTFyXabqDXGRoYiwpaNW1jy4RLWrFqDWITBVw3m5gk307ln53Nej2Z9OqzYB1sJKxeKyN+VvDYBmCEiTWf3uVrSIyvnGbPJWAK9+X1IWg3KEToMg773QVg/q+9TZDJbmPnDbt77PZmeYb7Mv6knrXzcrPoZ1tLYuy4fPXSULz/4kq8//pq8nDy6XdCNmyfcTP/L+uu+LU2YDiv2oVHCilLqYeDh8m9DgGNA6WmnuQH+wOcicpM1C7QlOqycxzL3wpYPjBGXoixo1d0ILTHXgpN12+iv2HGIx5dsx9XZkTfG9eCiKNtbxtvYYeW4woJCvvv8Oxa/u5iDaQcJaxvGLffdwvDrhuPs0nT72JyvdFixD40VVq7B2E0Z4FZgJUZgOVkJsAtYJCKF1izQluiwolFaCDs+h00LjR2jPVtCnzuh1x3gYb1BxcSj+Uz4ZAt7j+Xz6BXRTBgQaVObIdpKWDmurKyMNSvX8NHCj4jbEUdQcBC3TryVa268hmauzRq7PK2GdFixD7ZwG+h94P9EJNmaRTQVOqxoFSwWSPoVNs03Hp1cjaXPFz8A/taZY15QUsaUr3fw/Y5DDO0YxKtju9lM11tbCyvHiQgb125k0dxF7Ni8g4AWAdw84WbG3DwGN3fbvKWmnaDDin1o9LByvtNhRavU0TjYtAC2fw6WMogZDZc8DC3rP+lTRPjwjxSmr4gjxM+NhTf1olOwtxWKrh9bDSvHiQhb/tjCormL2PzHZnz9fbnp7pu4/rbr8fSy2x1BmjwdVuyDTYQVpdStwDggDDi92YGISKSVarM5OqxoZ5V7yBhp2fw+lOZDuyug/yMQ1rfeb70lNZP7Pt1KdqGJGdd24bpejbtk19bDysn++fsf/jfvf/yx5g+8fb258b83cuMdN+Lt2/ihTzuVDiv2odHDilLqWeB5YGf5UXL6OfbYbl8vXdZqpTAT/l5kzGspyoSwi2DAoxB5ab1WEKXnl3D/4m1s3JvBuD5hTB3RCVfnxln50pTCynG7tu/ivXnvse7HdXj5eHHbpNu44fYbdIM5G6LDin2whbCSAiwTkYerO9ce6ZEVrVZKC4zVQ3+8AbkHILQvDH4S2g6sc2gpM1t49ecEFq5NomtrHxbc1JPWfu5WLrx6TTGsHJcQm8D8mfPZ8OsGWrRswd2P3s3V11+Nk1PT2e7AXumwYh9sIazkAdeIyK/WLKKp0GFFq5OyEtj2Mfz2KuQdhDYXw+CnIPySOr/lT7GHmfzldhwdFfNu7MHA9oFWLNggIpSYS8gtzSWvNI+80ryKr2OPxVJiLsFkMWEymzBZTJRaSiu+N4sZQRARBMEiFkQEC8aji6MLzRybVRwuDi64OLrg7OBMM8dmODs608zBeHR3csevmR9+rn44O1hvgvGWjVt4Y8Yb7Ny2k7bt2jLxyYkMvHxgk9tU0p7osGIfbCGsfA+sFpHXrFlEU6HDilYvpmLY+iGsn2N0xg3vb4SWNhfV6e1S0guY8MkW4o/k8dgV0dw7MLLKH7QiQlFZEZnFmWQVZ5FVknXi62Lj68ziTLJLsk8JJSaLqdo6HJUjzg7OxuFoPDoqRxyUAwqFUgqFMr4v/9pkMVFSVkKJuYQSSwml5lLMYq72s7xdvPFv5o+fqx+BboG0cG9BkFsQQe5BBLgG4FjLTShFhDU/rOHNl95k39599Ozbk0emPUKHLh1q9T6adeiwYh9sIaxEAUuB2Rj9VjJPP0dELFarzsbosKJZhanIaDC3fg4UHIXIITB0GrTqWuu3Kiwt49ElW/ghbjf9op24vo83OaUZHCk8wrGiYxwtPMrRwqNkFGVQbC6u9D1cHFzwc/XD39Uf32a++DTzwcvFq+LwdvHG28X7lOf2Zu7F1dHVCCZW2qXabDFTaiml1Fx64tFcSkFZwYlAVZJZEazSi9MpMZ+YNueoHAn1DKWtT1vaehtHsEdwjeorKyvjm8Xf8NYrb5GTlcOIG0Zw3+P30TzI9hry2TMdVuyDLYSV40GkqotEROz2xq8OK5pVlRYaE3HXvwrF2dBlLFz6NPiFn3JaUVkRh/IPcSD/AAfzD3Kw4KDxmH+QA/kHyCjOOOOt3ZzcCHIPooV7CwLdA2nu2hx/N3/8mhmhxM/VryKguDu51/rWhy3MWRERckpzOFJ4hKOFRzlUeIiU3BSSc5MpKisCjCDWxrsNEd4RRPtFE+0XjYezR5XvmZ+bz3vz3uOz9z7DxcWF2++/nfF3jdeN5c4RHVbsgy2ElWlUHVQAEJHn61mTzdJhRWsI5oIMjqyfyb5/P2efo2J/aA8O+IdysCiDgwUHySw+dQDTycGJYI9gWnm2IsQzhGCPYFp6tORYtitv/nwEc6k388b249KOQQ1Wsy2ElapYxMKRwiMk5yaTnJNMcm4yqXmpmCwmFIowrzA6+nekg18Hov2icXM6s1lcWnIa86bPY+2qtbRq3YoHnn6AoSOG6vksDUyHFfvQ6GHlfKfDilZXZZYyDuUfYl/ePuPI3UdaXhr78vaxP2//KXNDXEQILrMQ7NWa4JALCfYOI9gzmBDPEFp5tCLQPRAHVfmOzGmZhdzz8RbiDucy+bL23DcoqkHa9NtyWKmMyWIiKTuJuKw4dmfuJikniTIpw1E5Eu0XTbfm3ege2J0g91MD3uYNm5kzbQ4JuxK44OILmPLiFMKjwhvnD3Ee0GHFPuiw0sh0WNHOxiIWDhUcIiUnhZTcFPblGsEkLS+NA3kHKJOyinPdnNwI8wojzDuMUK/QU75uUZCNw5rpELccvFoZ81m6jAWHygPK6YpKzTy5dAff/HOQK2KCeHVsdzybWffubFMLK6crNZeSmJ3Izsyd/HPsHw4WHASgpXtLugd2p1vzbrTzbYeTgxNms5llny5j/sz5FBUWccu9t3DHA3fo/iwNQIcV+2ATYUUp1QN4FhgA+AJ9RGSrUupF4DcRWWXNAm2BbgqnnazQVEhqbirJOcnGHIkc43ZDam7qKZNYPZ09jSDiHVYRRo4/BrgGVH9LYd8mWPUEHNwGIb3gypkQ2qdGNYoI729IYcbKONo29+Dtm3sRGWi9NvNNPayc7mjhUbanb2d7+nZ2Z+6mTMpwc3Kja0BXegf1pmvzruRn5jPvhXms/HolIWEhPDb9MS4ZUvfl59qZdFixD40eVpRSlwC/AHvLHycBvcvDynSgs4iMOtt7NGV6ZOX8ISIcKzpWEUSOHym5KRwqOFRxnoNyIMQz5MQKFJ+2hPuEE+4djr+rf/3nOFgssOML+GWasdy583Vw2fPgU7N2+xuTMpi4eCumMguv3dCdoZ2sM4/F3sLKyYrLionNjGX7se1sO7aNPFMeLg4udA/sTu+g3pgTzcx5eg7Je5IZPGwwU2ZMoXkLvWrIGnRYsQ+2EFZ+BzKAUYAjUMqJsDIamCsiYdYs0JbosGJ/RIQjhUdIyk4iMTuRxOxEkrKT2JuzlwJTQcV57k7uRiA56Qj3DifMO4xmjudgpUhJPmyYB3+8Dihjd+eLHwKX6rvXHsguYsLHW9h5MIdHL4/mvkFV92OpKXsOKyczW8zEZ8fz95G/2XJ0C7mlubg4uNDFvwuF2wtZ/vJyXBxceOyFx7jy2iv1BNx60mHFPthCWCkERovIKqWUI2DiRFgZAPwoIna7D7sOK02XiJBelH5KIEnMTmRv9l7yTHkV5/m7+tPOtx0RvhFE+EQQ7hNOW++2tHBvYRs/iLLT4JepsPNr8AmFK16EjiOqbd9fbDIz5esdfPvPQUZ2C2bWdV3rta/Q+RJWTmYRC/FZRnDZfHQzuaW5OCtnTHEmkpcn07N1T55+6Wndm6UedFixD7YQVjKBO0VkaSVh5QZgnoi0tGaBtkSHlaYhoyjjjJGSxOxEcktzK87xbeZLlG8Ukb6Rpzz6ufo1YuW1kPoHrHgUjsYaTeWumgXNo856iYjw1rq9zPpxN52DfXjnll608qnb7xbnY1g52fHg8ufhP/nryF8UlhVSll1G4ZZCxl44lhtG3WAb4baJ0WHFPthCWPkOY1Lt4PKnTEAvEdmmlPoJSBeR8dYs0JbosGJbisuKScpOIiErgYSsBOKz4knMSiSrJKviHG8X74owcnIwqdEEV1tnLoO/34U1L0JZMVx0P/SfDC5VNz0D+GXXER78fBvuzZx4++Ze9AyrfUA738PKyUrNpWxP387PNoWbJgAAIABJREFUiT+TkJcAjuB02Imb+t1E/4j+ODnYbZ9Mq9NhxT7YQljpBmwAUoAlGKuC3gC6Ab2AC0Qk3poF2hIdVhqHiHC44HBFIDkeTlJzU7GU7+7g5uRGO992tPNrd0owCXQLbPqhpDp5R+Dn52DH5+DdGq58ETqOPOutoYQjedz54WYO5xbz0rVdGNOrZhN2j9NhpXJZRVnM/3Y+cRKHc3Nn3HFnaNuhDGo9CH9X/8Yuz+bpsGIfGj2slBfRE3gFY+myI2AB1gOPiMg2axZna3RYaXiFpkISsxMrAsnxI6/0xLySEM8Qov2iae/fnvZ+7Yn2i6a1V+sqG6WdN06+NdT+Khj+KviEVHl6VkEp9326lY17M7h7QARTruyAYw0byOmwcnbxsfFMnTsVU7QJr65eOCgHegf1Znj4cNp4t2ns8myWDiv2wSbCyknFuAL+QLaIFFqzKFulw4p1HSs8Rlym0VF0d+ZuErIS2Je7Dynf0cHdyZ32fu0rjmj/aKJ8o/B0sV6/ELtjLoNNC4xbQw5OMHQq9P5vlQ3lTGYLL3y/i482pjIoOpDXx/XA29W52o/RYaV6xUXFvPnSm3z97ddEXBeBRz8PSqWUGP8YhocPp6N/R/sf9aslHVbsQ6OHFaXUVuAjYLGIHLVmIeeSUioS+BBoARQAd4lItSlEh5W6ERH25+2vCCZxmXHEZcSdsgFfa8/WdPDvYAST8hGTEM8QPVpSV5nJ8P1DsHcthF4II9+AwOgqT//0z1SmfhtLWIA77916AW2bn33eiw4rNffHmj94/uHnKSwrZNRLo0j1TiW3NJe23m25NvJaugR00aGlnA4r9sEWwsoPwJDyb3/BCC7fiEjle8/bKKXUz8AXIrJIKXUZ8CbQQar5y9BhpXomi4nknGTiMk4Ek/jMePJN+QA4KkcifCOMjeT8O9DRvyPR/tF4uXg1cuV2SAS2fwarngRTIfR/FC55GJxcKj19094M7v1kC2aLMP+mnvRvF1jlW+uwUjvpR9J56r6n2LppK6NuGcWFd1/Ij/t/JL04nXa+7RgdOZqO/h0bu8xGp8OKfWj0sFJeRAtgPPAfoCeQhzHZ9hMRWWPN4k76zNbAFKA3xmReN6CtiKRUcm4o8BpwGaAwQtVDIrKv/PVAjA68/iJiKn8uHriputEVHVZOVVRWREJWArszdleMmuzJ2kOppRQAV0dX2vu3PyWYRPlFnZsmatoJ+ceMtv07l0BgRxi1AEJ6VnpqWmYhd320mYQjeTx3dSduu7htpefpsFJ7ZWVlLHh5AR8t+IgOXTow4+0ZJDkk8V3yd2SXZNPJvxOjI0cT5Xv2Jej2TIcV+2ATYeWUi5XqCNyMEV5Cgf0iYvXZY0qpQcAXwBaMSb2XU0lYUUq5A9uBEuAZQIDpgDvQVUQKlFK9MG5jRZ903U/AWyKy9Gx1nM9hJackxxgpyYirCCYpuSkVq3G8XbwrQkmHACOYhHuH4+hQ98ZjmpUl/AjLH4L8I9D/ERjweKWjLPklZTz0+T/8EneEW/q14bmrO+HkeOrtOB1W6m7tj2uZ9tA0lFK8tPAlel7SkzX71/B98vfkmfLo1rwbo6NG08br/JuIq8OKfbC5sAIVE23HADOBYBGx+k8npZSDiPFTUSl1J/AulYeVB4E5QLSIJJY/1xbYAzwuInOqCCs/Awt1WDFkFWexK2MXsRmx7MrYRVxGXMWutAAt3FvQyb8THQI6VIyYtPJope+7NwVF2cZtoe2LIagLXLsQWnY54zSzRZi1ajdv/7aXAe0DeXP8qRNvdVipn/0p+5n838kkJyTz4LMPMv6u8ZSYS/gl7Rd+SPmBwrJCBoYMZEzUmPPqFqkOK/bBpsKKUupSjFGV0YAn8BfwsYgssF55lX7u2cLKasBVRC4+7fl1ACIysIrbQAnA+PPxNlBmcSa7Mnadcpy8UV+YVxgdAzrS0d84OgR00P0i7MHulbD8QSjKhIFTjLksjmeuAvri7308vWwnbZt78L/bLiDU39iLSIeV+ivIL2Dqg1NZu2otI8aO4MmZT+LSzIUCUwHf7f2On9N+xtXRldGRoxncevB5MUqpw4p9aPSwopTqjDFXZTwQAqQCn2CElD3WLOwsNZwtrBwGvhWRe057fgFwvYgEln+/GvhcRN4tn2C7AGhv7xNsjweT2HRjxGRX5i4OFxyueL2Ndxs6+XeiU4BxdAjogLeLdyNWrDWowkxY+Zgxl6VVd7j2LWhx5iTPP5LSufeTrTg5KN65pRe92vjrsGIlFouFd+e8y7uvvUuXnl14ZdErFXsLHcg/wKfxn7IrcxetPVtzS4dbaO/XvpErblg6rNgHWwgrFiAH+AojoKy3ZjE1rOFsYaUUmCMiT5z2/HTgCRFxKv++HcbS5eZAIXC3iPxVxefdDdwNEBYW1is1NdW6f6AGklGUceqISRXBJKZ5jBFM/DucV8PN2kl2fQvfPwwleTD4KbjoATjtt/i9x/K544O/OZhdzKzruhIadLiKN9PqYvWK1Ux9cCpePl7Mfm82Md1jAGPZ/5ajW/gs4TMyijMYEDKAse3G4ulsn72GdFixD7YQVq4HvhOREmsWURs1CCuvisiTpz0/A5hyPKzUla2OrBwPJsfnmOzK2MWRwiMVr4d7h9MxoCMxATqYaFXIPwYrHoa45dDmYmOUxTfslFOyCkqZ8MkW/kzO5Lo+Htx4oaeep2RFCbEJTL5jMhnHMnh29rNcNfqqitdKzCV8u/dbVqWuwsPJg3HR4+jXsp/d/f3rsGIfGj2s2IJqwsoRjL4vZ70NVIfPHAGMiIqKumvPnnNyt6tK2cXZxGbEsjN9Z0U4OT2YHL+N0ymgEx39O+qOr1rNiMD2z2Hlo6Ac4eo50OW6U04pLbPwzDf/8uXm/VzczpWJl/nQzMm+fmA2pqyMLKbcM4WtG7dyy323MOnJSTic1H04LS+ND+I+ICkniV4tenFbx9vs6hcPHVbsg02EFaVUD4wNDAdg7MDcR0S2KqVeBH4TkVXWLLCSzz9bWPkVcBGRS057fi3Gn3VgfT77XI+sFJoKK0ZMdqbvZGf6Tvbn7wdAoYxbOTqYaNaWmQxL74b9f0GXsTB8Nrj6VLwsIjz3/Xo+2ZBHVJAzU672xc/D/id/nitlpjJmPzebJR8t4bKRl/H83OdxaXZiiblFLKxKXcXSxKW4O7tzR6c76B7YvRErth4dVuxDQ4SVWt0WUUpdgtFkbS+wGJh00ssWYALQoGGlGt8Bs5VSESKyF0ApFQ5cDDxxlusanclsIiErwQglGUYw2Zuzt6KPSSuPVnRu3pnro6+nc0BnOgV00sFEaxj+beH2H2D9q7DuZdi3CUa/A236AaCUYlQvD1r5OjLvxxye+DKDJ0f4Ed68+j2FtOo5OTsx5cUphLQJYd4L88jOyOaV917B08v4791BOTAsfBhdArrw9s63mfvPXAaGDGRc+3G4Ork2cvWa1jBqO2fldyADGIXRnK0U6F0+sjIamCsiYWd7jzoXqtTx8eghGKHoPuAYcExE1pWf44HRFK6IE03hXgC8MJrC5dfxs616G8hsMZOSm8LO9J38m/4vsemxxGfFY7KYAPB39ScmIIbOzTvTuXlnYgJiCHALqPfnalqtpf0NS++E7H3Qf7KxzNnRuWI10N6jJl5ankVRqfD41b50DdXdia1p5dcref6R54lsH8nrn7xesVLoOJPFxLKkZfyQ8gNB7kHc2/XeJt1MTo+s2IdGvw2klCoERovIKqWUI2DiRFgZAPwoIm7WLPCkz66q0HUiMuik88I4td3+aox2+yn1raEut4FEhAP5B9iZsZPYdON2zq6MXRSWGRtVezh70CmgE50DOhPTPIYuzbvoBmuabSnJgx+egH8+gZBeMOY9thRnVryckW9m+rdZHMwq4/7LfbikfYP8L+C8tXHtRh6/63H8Avx4/ZPXCY8KP+OcuMw43t75NvmmfMa1H8elrS9tkv8P0WHFPthCWMkE7hSRpZWElRuAeSLS0poF2pKahJX0onQjlJTfyolNjyWrJAsAZwdnOvh3ICYghi6BXegc0Jlwn3C9s7DWNMR+A8sfABGS+j9MdsSJKWAFJRZe/j6L2AMmbuvvxYgeZ9+1WaudXdt38eDND2KxWJj30Tw69+x8xjm5pbks2rmIHRk76N2iN/+N+S9uTk0rOOqwYh9sIax8hzGpdnD5Uyagl4hsK99fJ11ExluzQFtQ1W2gvNI8dmXsqpj8ujNjZ0UvEwflQKRvJJ0Dym/lNI+hvW97nCvpEqppTUb2PvjqdjiwmaOdrmF/33uR8v2FSsuE13/KZmNiCSN7uHPzJV44NMHf7m1VWnIak26aROaxTOZ+NJde/c78wW4RCz+m/shXiV/R0r0lD3R7gJYeTef3Rx1W7IMthJVuwAYgBWOn5WeBNzB2Qu4FXCAi8dYs0JbEdI+Rpz57qiKcpOSmVLwW6hVacSunc/POdPTviLuze+MVq2kNpayUw9/dR8sdX1EYEMXeoc9R4tMaMPYUev+3PH7YUUj/aFcmDvXB2VEHFmtJP5LOvTfcy8G0g8x5fw4XDriw0vPiMuOYv2M+ZjEzofMEugV2O8eV1o0OK/ah0cNKeRE9gVcwli47YqwCWg88IiLbrFmcrXFr6yZR06IIdAs0QknAiQmwvq6+jV2epp0zWw5uwSd1I+FrX0ZZTKT2f4SsqCGAMU9r2ZYCPv0jn66hLjw+3Bc3F32r01oy0zO578b72Ld3H7PencUlQy6p9LxjRcd4Y/sbpOWlMTpqNFeHX23z81h0WLEPNhFWTirGFfAHskWk0JpF2Zrjt4FCIkLu2rJjC0EeQY1dkqY1quOrgZzzjxKxejqeR3ZyrMNw0i66v+K20K+7Clm4Opc2zZ14eqSf7sViRdmZ2UwaP4nE3YnMfGsmg64cVOl5JeYS3t/1PpsOb6J3i97cGXOnTS9v1mHFPthUWDkf2Wq7fU07107ZyNBSRvDm92n1z2cUBEaTdNk0TJ5GoN+aUsLsldn4uDvw7Cg/gn3rteOFdpK8nDzuv+l+4nbEMWPBDIZePbTS80SEVamr+HLPl4R4hvBAtwdo4d7iHFdbMzqs2IeGCCvVjs0qpX5VSnWo6RsqpRzKr2lXv9I0TWsSHJw42OcuEi9/Adfs/XRcOgGvA1sB6BnejOdH+1FcauHpLzPYc7i0kYu1H14+Xsz/fD6de3bm6YlPs/bHtZWep5TiqvCrmNxzMlnFWTz/1/MkZiee22I1rZ5qciN5EEZTtZpSdbhG07QmLif8YuKuXUCZmy/tVj5O0D+fgwjtWrowY2wAri4OTF2axdaURtsH1e54eHow7+N5dOzSkSfueYLff/m9ynM7B3TmuQufw8PJg1lbZvHPsX/OYaWaVj81nfX2jVJqb00OYA9G51i7oZQaoZR6Jycnp7FL0TSbVuIbyu5RC8hq25/Wf71DxC/P41BaSLCvEy9e70+wnyMvLc/i1112Pc3tnPL08uSNT9+gXcd2PHbXY2xcu7HKc4Pcg3imzzMEewTz+vbXWX9g/TmsVNPqrto5K0qp9+v43lNFZF8dr7VJes6KphlOmbNSGRGCdnxJyF/vUuwTStLlz1PiG0ZRqYVZK7LZkVbKfy7y5Nreen8ra8nJyuHeG+4lNSmV1z54jT79+1R5blFZEW9uf5PYzFjGRI2xmZVCes6KfdATbBuZDiuaZqg2rJTzOrCVtqtfwMFsInnwE+SEX4LJLLzxcw4bEooZ3duD8f08beIHpT3Izszmnuvv4UDqARZ8voCuvbtWeW6ZpYxFsYvYdHgTl4ddzrj24xr9n4MOK/ahUSbYapqm1VVeSE/iRr9FsW8oUT89R8utH+PsAA9e7sNlnd1YurmARevysOhfmqzC19+XBZ8vIDAokIdufYjE3VVPpHVycOLuzndzWehl/LTvJz6M+7Bil3dNszU6rGia1qBMnkHEj5hHRtRQQja/T/iaF3GymLhnsDcje7qzakchb/6cg9miA4s1BAQGMP/z+TRzbcak8ZM4sO9Alec6KAfGR49nePhw1h5Yy3ux72G2mM9htZpWMzqs1ICeYKtp9SNOLqQMfpIDF/yXgMTVtP/+YZyLsrjlYi/G9fVk3e5iXv0hG1OZDizWEBwazJuL36S0pJSJ4yaSfjS9ynOVUlwXdR3XRl7LhkMbeHvn2zqwaDZHh5UaEJHlInK3j49PY5eiaU2XUhzucRNJl03DLTOZjsvuwz1zL9f18eT2AV78mVTCS99nUWzStyKsITI6knkfzSPjaAYP3PQAeTl5VZ6rlOKaiGu4od0N/HXkL97Z+Y4OLJpN0WFF07RzKrvtAOJHzgUsRH/3AN77/uTq7h7cN8Sbf9NKeeGbLApKdGCxhi69uvDKe6+wd89eJt8xmdKSszfluyr8Ksa2G8ufR/5kUewiPYdFsxk6rGiads4VNW/P7lHzKfEOIerHp2m+6zuGxLjz8JW+JB4xMXVpJjmF+gelNfQd0Jdpc6exddNWpj08DYvl7H+vw8KHMTpyNBsPb+R/u/6nA4tmE+ocVpRSEUqp95VS7yilIk56votS6lnrlKdpmr0yeQQSP3IeOaF9aPP7XEI2LeSiKBemXO3Hgcwynv06g4x8fSvCGq4cdSWTnprET9/+xJsvvVnt+SMjRjIqYhS/H/ydD+I+0IFFa3T1GVn5CHABDgHvKqWuUUrtBjYBXaxRnKZp9s3i7EbS5S9wNGYULXd8RcTP0+jVWnjmGn8y8y08sySTwzlljV2mXbj1vlu57pbr+GjBR3z5wZfVnn9NxDWMaDuC3w78xie7P0H35NIaU33CSoCI3CQiU4H/Al8CHwMtRGSsVarTNM3+OTiSdvEDpPWbiG/KBqKXP0I3/3ymjfansNQILPsyTI1dZZOnlOLRFx6l/2X9mf3s7LO25T9+/ujI0QwLH8av+3/ls4TPdGDRGk19wkrR8S9EJAVIFpEZIlJQ76psjF66rGkN72iXMSRd/n+4ZqXQ4ZuJdHY5wAtj/BGB577OJCVdB5b6cnJyYsb8GURGR/LkvU+SmpR61vOVUlwfdT2XhRmN41akrDhHlWraqeoTVqKVUt8opWYopcYDdntTUy9d1rRzIyf8YhJGvIaD2UT0t/cTU/wv06/zx9lR8fyyLNIy9C2h+nL3cOfV91/FycmJR25/5KxLmsEILOPaj6Nvy74sSVzCugPrzlGlmnZCfcLKIGAV0AJ4FIhUSmUppX5USj1vjeI0TTv/FAZGE3ftfEweLWi3cgqdj/zMtNH+OCiYtiyTg1k6sNRXcGgws96dxf7U/Tw98WnM5rNPZHZQDtwZcyedAzrzwa4P2HZ02zmqVNMMdQ4rIvK3iLwlIneJSE/AC7gMWAq0slaBmqadf0yeQey+Zh55wd0JX/cKPdKWMu1afywCU5dmcjhbB5b66tm3J1NmTOGPNX/UaIWQk4MTk7pOItw7nAX/LiAhK+EcVKlphmrDiqrhNpwiUioim0XkbRG5u/6laZp2PrO4eJJ45YtkRgwkdNNb9En+hKmjfCk1C1OXZXI0Vy9rrq/R/xnN9bddz8cLP2bFkurno7g6ufJwj4cJcA1g7j9z2Z+//xxUqWk1G1nJVUrppciapp1z4uhM8qXPkB59FcFbP+KSPe8y9Ro/CkuEaUszycjTgaW+Jk+bTO+LejPj8Rns3Laz2vO9XbyZ3HMyLg4uzN46m/Siqvcd0jRrqUlY8QDcjn+jlHJQSm1RSrU9+SSllKtSytvaBWqadp5zcCR1wGSOdB5N0M6vGRT/Os9d401ukYVpyzLJKtCBpT6cnJ2Y+fZMAoMCeezOx8hMz6z2mkC3QCb3nEyJuYTZW2eTV3r2SbqaVl91mbOigB6A32nPdwOq/7dc0zSttpQD+/tN5GDPmwncvZLLY2fz3AhPMgssTFuaSU6hDiz14evvy6x3Z5GbncszE5+pdsItQKhXKA91f4j04nRe3/46JoteWq41HGvvDVSj+S2apmm1phSHet/O/gvvxn/vGq7a+SJPD3PnaJ6ZacuyyCu22+4J50R052gen/E4f/3+F+/Oebdm1/hFc1fMXezJ3sNHcR/ppnFag9EbGdaAbgqnabbjSLcbSb3kQXz3bWRE3Es8NcyDg1llvPRdFiUm/cOyPq658RpG3DCCRXMXseHXDTW65sKWFzKy7UjWH1zPT/t+auAKtfNVTcPKef1/AN0UTtNsS3qna0i95GF8923i2rgXmTzUnYTDJl5ZmUWZ+bz+31W9TZkxhXYd2/Hc/c9xaP+hGl0zKnIUvVr04vOEz/k3/d8GrlA7H9U0rKxWSm1USr0FTMQIL84NV5amadrZpXcaQWr/R/BJ+5PxyS9y7yBXtqWWMv+XHCz6dkSdubq58vI7L1NmLuOJCU9QWlJa7TUOyoG7O99Na8/WLPh3AQcLDp6DSrXzSU3Cyt0YOyybgfHAXIy5KRuUUglKqa+VUlOBoQ1XpqZp2pnSO15NyoBH8U77m3v2v8jNFzrzW3wxH6zP0/Mn6iEsIoypc6YSuy2W1/7vtRpd08yxGQ92fxBnB2fmbZtHgcnutonTGlG1YUVEFonIJBG5RES8gWjgBmAWkAj0A6YCLzRopZqmaZXI6DCM1AGP4r1/C4+lv8Soro6s+KeQpZv1D8v6uHTYpfznnv/w1Qdf8dO3NZuL0tytOfd3u5/04nTe3/W+Doya1dR6gq2I7BGRr0TkKREZJiLBQEtgGPCU1SvUNE2rRkaHq0gd+BheB7YxrehVBrd3ZPHGfNbEFVV/sValSU9Oomuvrrz4xIscPnC4Rte0823HmKgxbD66mTX71zRwhdr5wiqrgUTkqIisEpGXrfF+mqZptZURfSX7+j+Mb9pfzHZcSPfWjixcncO21JLGLq3JcnJ24v9e/z8sZgvPPfhcjfqvAFzZ5kq6BnRlccJiUvNSG7hK7XxQk72Bfq3FsfpcFK1pmlaZ9I5Xs//CewhIXstCvw8J9XNk9sps9h7VDcvqqnV4ax594VG2btzKp+98WqNrHJQDd3a+E09nTxbuWEhxWXEDV6nZu5qMrDhgTKg9fnQABgHhGG34w8u/j0Y3hdM0rZEd6XYDh3rcRKs9K/kgbAmezWDGd1kcydE7NdfViLEjuHTYpSx4eQHxO+NrdI23izcTOk/gSOERPtz9oZ6/otVLTSbYDhKRwSIyGJgHmIC+IhIhIv1EJAJjkq2p/HWbp5R6SikVr5SyKKVGNXY9mqZZ18Hed3C00zWEx33FR9E/UWYWpn+bRW6R7nJbF0opnnr5KXz9fXlm0jMUF9VspKSDfwdGRYxi46GN/H7o9wauUrNntZ2z8gLwrIj8dfKTIvInMA2YbqW6GtpqjAnBvzV2IZqmNQClSLv4fjKihtIp9n3e77KeY3lmZi7XXW7rytffl2mvTSN5TzJvvPhGja8bETGCjn4d+TjuYw7m6/4rWt3UNqy0A45V8dpRIKq2BSilWiul3ihvOleolBKlVHgV54YqpZYopXKUUrlKqaVKqbDafqaI/CkiSbW9TtO0JkQ5kDJoCtltLuaCnQtY0G0rCYdNzPspWzeNq6O+A/sy7r/j+OJ/X/DX+r+qv4ATDeOaOTZjwb8L9IaHWp3UNqwkA/dU8do9QEodaogCxgJZwPqqTlJKuQO/YsyZuRW4GSM8rVFKedThczVNs3cOjuwd8iy5wT0YGjeHFzvH8mdSCZ/+kd/YlTVZE5+cSFjbMKY/Np3CgsIaXePn6sedMXeyP38/y5KWNXCFmj2qbVh5HhihlNqplJqmlLq3/HEnMBzjVlBt/SYiQSIyDPjqLOfdBUQAo0TkGxH5FhgJtOGkAKWU2qqUSq/iCK1DfZqmNWHi5ELS5S9QEBjNDXtncX9kIt9sKeCX2Jr9oNVO5ermynNzjH2D3nzpzRpf1y2wGwOCB/BDyg8kZic2YIWaPapVWBGRz4ErgBzgSWB++WM2cIWIfFHbAkSkpjPeRgKbRKTi33IRSQY2ANec9FxPEWlexZFW2/o0TWv6LC7uJF75EiU+ITx49GVGBh/knTW5/Jume7DURfc+3bnhjhv48v0v2bppa42vGxc9Dn9XfxbFLqLErP/utZqrSZ+Vfid/LyK/iMjFGMuWWwJu5a34G7rHSgyws5LnY4FODfWhSqm7lVKblVKbjx2rarqOpmm2zuzqzZ6rZmJxceOV4pl0887mlZXZHMjSS5rrYuITEwkJC+GFyS/UeHWQm5Mb/+30Xw4XHubrxK8buELNntRkZOV3pdQhpdTbSqkrlVLOYIyIlHeuPVdrAf0x5rWcLhPwq80bKaWeUUrtx1hyvUgptV8p1bKyc0XkHRHpLSK9AwMDa120pmm2w+TZgj1XzcTJVMSHzV7GTxXy4nd6SXNduLm78cwrz5CWksbCWQtrfF2ngE4MCR3CT/t+Ynfm7gasULMnNQkrIcD/YcwN+QZIV0p9oZS6USnl3aDVnamyKfy1bkQnItNFpLWINCu/PdRaRGq28YWmaU1asX8ESZc/j0feAb4OmEtefjGzVmRhKtMrhGrrgksuYPR/RrP43cX8u+XfGl83tt1YWri14L1d7+nutlqN1KQp3GERWSgiVwKBwASM0PAWcFQp9aNSaoJSKriBa83CGF05nR+Vj7hYjVJqhFLqnZycnIb8GE3TzpG8kJ6kDJpCUMa/LAt+j90HS3hnba7usloHDzzzAC1ateCFR1/AVFqzZcnNHJtxZ8ydpBelsyRxSQNXqNmD2k6wzRORz0TkRozgci3GcuZngTSl1J9KqScaoE4w5qbEVPJ8J2BXA30mACKyXETu9vHxaciP0TTtHMqKGsL+PnfT/thvvN96Cb/uKmLFdr1CqLY8vTyZMmMKexP28um7Nds7CKC9X3sGtx7M6rTVerNDrVp13nVZRExzy+SgAAAdB0lEQVQi8oOITBCREOASYC1wi7WKO813QF+lVMTxJ8qbx11c/lqD0SMrmmafjnS7gaOdrmFQ+jc8G/QrH67P4x+9S3OtDbh8AAOvGMii1xZxaP+hGl83JmoMns6efLL7EyznbPqj1hTVOaycTkQ2isgUEan1yhyl1HVKqeuAXuVPXVX+3MCTTnsXo+nct0qpa5RSI4FvgTTg7XqWf1Z6ZEXT7JRSpF00ieywftyR+z/G+OxkzqpsDuoVQrX22AuPATD7udk1vsbD2YOx7cayJ3sPfxz6o6FK0+xAncOKUipPKbVOKfWaUuo/SqlOSqm67rr8Vfkxofz7BeXfP3/8BBEpAC4FEoCPgU8xbkFdKiK6HaWmaXXj4EjykGco8o/gpbK5dFD7eOn7LApK9G/6tdEypCV3PXIX635cx28/1XzbtYuDLybSJ5IvEr4gtzS3ASvUmrL6jKy0A17BaAh3A/ALkKeU2qCUer02byQiqopj0Gnn7RORMSLiLSJeIjJKRFLq8WeoEX0bSNPsm8XZjcQrZmBxceNDt1cx52Tw2qpszBY94bY2xt85noj2Ebzy7Cs17r3ioBy4ucPN5Jvymb9tfgNXqDVV9ZmzclhEvheR50VkhIgEA6OBZsBEq1VoA/RtIE2zfybPQBKvnIGbKZdl/q8Rl5rHJxvyGrusJsXZxZknXnqCQ/sP8fFbH9f4unDvcAa3Hszn8Z+TmKVb8WtnqvecFaVUN6XUC0qpOIyRluVAt3pXpmmado4VNW9P8qVP0yo/kcXN32H5tnzWxBU1dllNSs++PRkyfAgfzv+Qo4eO1vi6ayOvxcPJgzlb5jRgdVpTVZ85Ky8ppRIwJr7mAyNFpFv5SEtlbfE1TdNsXk74xezvO4Ge+RuZ5fc1b63OIf5QaWOX1aTc//T9mM1mFry8oMbXeLl4cXfXu1l/YD0bD25swOq0pqg+IytTgAyMibB/A3a7cY6es6Jp55ejXa7jWMcRXF+0jNvcf2PWimzS88yNXVaT0bpNa8bdOY7vv/qeuB1xNb5uXMdxhHiGMHvzbMwW/fetnVCfsHIh8CHQHpgNHFFKJZa34p9ilepshJ6zomnnGaXYd/H95LTuzVPmd+lm2snL32dRYtITbmvqjvvvwC/Aj9eef63GnYGbOTbjoZ4PkZCVwHdJDdo+S2ti6jPB9m8ReUtE7hKRnoAX/H97dx4fVX3ucfzzzcIqKIIVKEJEBbciaLUgWlFAEatyKyotmoqIWEu91WqvWrW4dPH2iiJSWwRFvdi6VYiirOJSAasiWMSFfVcWwyY7ee4f58TODQk5k5nJTCbP+/U6r5Bzfr8zz8yTYX5zfsuhH/A6cFSyAnTOubTIyWNJj7vYdUgrRtV9CG1Yzshpm31J/ogOanwQg28ezJzZc5jx2ozI9c4rOI8OzTow4sMRbN/jKwq7QDIG2D4laSEwBWhkZn8xs2sTD80559KrpM5BLOr1e3Ly83m+0QMsWLiOF9//Ot1h1Rh9ftyHtu3bMvy+4ezeFW3cjyRuPvVm1u9Yz9MLos8octktGSvYvgScDfwWeFDSwCScM6P4mBXnaq/djZqz+Nx7OWTfRp49+CFenFXMu4v9TsFR5OXlcdNvbmL18tU898Rzket1+lYnuh3RjSc/fpLNu/z/XZeExoqZvWRmq8xsKvB9YICkyxMPLXP4mBXnarevDz+eZd1upd2uTxnZaAwPT9nEyo2+JH8Unc/qTNdzujJ6+GiKNxZHrjek4xC27dnGE/OfSGF0rqZIRjdQgaQzJfUHridYAn9soud1zrlMUnzU2aw55Sp67nmLQXmv8t8TfUn+qH5x1y/Y8fUORj0wKnKd9oe2p9eRvRj3yTg27NiQwuhcTZCMbqB5BPfyuZJgYO1CgkaLc85llbUnX0Hxkd/nRnuG9lvnMGLKZkp8wG2ljjzmSH545Q/5+//+nVXLVkWu97OOP2NPyR4e++ixFEbnaoJkdAMdbGbfMbNe4cyge8zMr9s557KPcljW7b/Y2fRIHq33CBuXLeWFf/qA2ygG3jCQvPw8xgwfE7lOm8Zt6HN0H577/DnWbFuTwuhcpkvGlZWs5wNsnXOlSvLrs/jce8nNz2dcw2G89u6XfLDUB9xWptnhzehb2JeJL0xkxZIVketdd9J1CPHovEdTGJ3LdJEaK5LOkPQ7SfdI6lRBmaaSCpMbXmbwAbbOuVi7GzVnSc+hHL7vSx5rOJKHJxezZpMPuK1M4fWF1Klbh8cejN6t07xhcy5vfzlFi4tYunlpCqNzmazSxoqkfsAbwK3AHcB7kq4LjzWXdKOkt4AvAO/+cc7VCttadGBl1xv43r653JTzV+5/ZRM7dvuA2wNpelhTLr3qUiaPn8yyRcsi17vmO9dQN7cuI+eOTF1wLqNFubJyK/ARcCrQGhgA3CbpJoKZPw8QLLk/FrgoNWE651zm2XD8haw7/mKu5hVO2zzDV7iNoPCnhdStVzeuqytN6zfliuOuYPKyySzZtCSF0blMFaWxcgxwp5l9EK6n8jTwS8L7AQE/AFqEg2snpjBW55zLOCtP/xlbW3Tk/jqj2b74YybM8QG3B9KkaRP6Xd2PKROmsPizxZHrXXH8FdTNrctTC55KYXQuU0VprNRn/zsqTwl//tLMXjX/KuGcq61y8ljc8zeUNGzK2PoPMnXmShasjra0fG3Vf3B/GjRswOgHR0euc2i9Q7n4qIspWlzk667UQlFnA5VtjJR+dViWvFAyl88Gcs4dyL56B7Oo1300ztnBqHrDGfHaBjZv35fusDLWIYceQr+B/Zj68lQWfbIocr3CEwrZW7KXZz55JoXRuUwUtbHypqQ5kp6WdCvB2BQDasVoMp8N5JyrzM5D27Ks26/oYJ8zZM/TPDR5M/tK/KJzRfpf25+GjRrGve7KOa3P4dnPnvU7MtcyURor1wKPA9uBi4HfAS8AAiZLelXSvZIulvTt1IXqnHOZbVPbbnz5nb4U5k6mYPUbvPiej1+pSONDGtO3sC/TJ06Pa1Xbq064ii27t/DSopdSGJ3LNJU2VsxstJkNMbMzzKwx0B7oB/wBeB/oCPya4O7L0Vf6cc65LLTqe9ey9fAT+WPdx/jwn58yb8WudIeUsfpd3Y/cvFzGjRoXuU7Hb3XkpMNO4ukFT7O3xNe2qS3iXsHWzBaa2fNmdruZ9TazlkBzoDdwe9IjdM65miQnj6U97kL1GvBYvYd4bNJavtrm41fKc1jzw+h9SW+Kni2K647MA04YwOptq5m+YnoKo3OZJCnL7ZvZOjObZGb3J+N8zjlXk+1p2Iyl3e/iCL7gNyV/YdhrxT5+pQJXXnclu3ft5tnHn41cp9sR3fj2Qd/m2c+i13E1W5QVbF+PY/NmrnPOAdtansSa066hV867nL5uAn+dtS3dIWWkgqMLOOu8s3hu7HNs2xbtNcrNyaVvu76898V7LNnsi8TVBlGurOQQDKYt3Y4FugEFBGuwFIS/tw+PO+ecA77scDnFBWdye/5fWTXnA7/hYQUKf1rIlk1bGDMm+sygPkf3IS8nj+c/ez6FkblMEWWAbTczO9vMzgaGA3uAzmbW1sy6mFlboEu4f3hqw3XOuRpEYlm3W9jduAV/rvcwf5u6jOKvffxKWR2+24FO3+vEsGHD2Ls32qDZZvWb0aN1D4oWF7FzrzcCs128Y1buJVh6/5+xO83sXWAocF+S4nLOuaxQUucglpx7Nwfn7OC3JY/wyOSvKPFFv/fTf3B/VqxYwYQJEyLXubTdpWzZvYUpy6dUXtjVaPE2Vo5h/6X3S60Djk4snMzkK9g65xKx89C2rDrjBrrkLOCMtS9Q5PcP2s+ZPc6koKCAhx9+OHKdU5ufSkHjAl74/IUURuYyQbyNlaXA4AqODSZLl9/3FWydc4na2K4XG4/qzk35L/L57PdZ9OWedIeUUXJzcxkyZAhvvfUWc+fOjVRHEn2O7sOH6z5k5ZaVKY7QpVO8jZW7gQslzZc0VNJPw5/zgQsIuoKcc86VJbHizBvZ2agFw/MfYcxrK9mxu1bcsSSyq6++mgYNGjBixIjIdS5oewFCvLLklRRG5tItrsaKmf0NOA/YDNwGjAx/bgLOMzOf9O6ccxUoqdOA5T3u5DBt4ZYdjzL6De9ajtWkSRMKCwsZN24cGzZEu7Ny84bNOa35aby85GXMxwJlraqsYDvNzLoSTFtuDtQPl+L3NVacc64S2w9rx5rOg+mZ+wEFCyfw9mc70h1SRhkyZAi7du1i9OjRketceNSFrNy6knnr56UwMpdOVV7B1sxKwpVr/Tqmc87FYd2JP6S4dRduz3+Gt2fM48vNfo+bUieccALdu3dn5MiRkacx92jTg3q59ShaXJTi6Fy6JGW5feecc3GQWN7tV5TUP4RhepjRU77w5fhj3HDDDaxatYrx48dHKt8wvyHntD6HScsmsXvf7hRH59LBGyvOOZcG++odzPIev+YIraPf+jFM8OnM37jgggto3bo1o0aNil6n7QVs3b2VWWtmpTAyly7eWHHOuTTZ1uIkvuj0Yy7Le5Mt705jyTqfzgzBNOaBAwcydepUliyJdu+fLi260KhOI18gLkvVusaKpHqSxkv6RNJcSZMltU13XM652mntKYVsaXos9+WP4ZlJS9i117uDIJjGnJOTE/l+Qfm5+Zx9xNnMWDHDu4KyUK1rrIQeNbPjzKwj8DIQfdi5c84lU04eK3rcToPcfdzy9QieecenMwO0atWK3r1788QTT0QeaHtewXls3bOV2Wtnpzg6V93S3liR1ErSCEmzJG2XZJIKKih7hKQXJG2WtEXS3yW1jufxzGynmU2O2TUb8Csrzrm02XVwK1Z3HcLpuQtoOf955q3Yle6QMsKgQYNYu3YtEydOjFS+S4suNMpvxORlkysv7GqUtDdWCO4ndBlQDLxdUSFJDYDXgWOBnwBXEtyraIakhgk8/s+B6HfOcs65FNjY/nw2tjmTW/KfY9qUuWzd6atC9O7dm5YtW0YeaJufm8/Zrc9mxsoZ7Nnn43+ySSY0Vt4ys8PNrDfw/AHKDSK4AtLHzMab2QTgIqANMfcrkjRH0oYKtiNiTyjpNqAdwSq8zjmXPhIrz7qJvfUO4Z59I3hyekX3jK098vLyGDBgAJMmTWLlymj3/jm3zbnBrKC1Pisom6S9sRLHonIXAbPNbFFM3aXAO8DFMftONrNmFWzf/LVLuhm4BDjfzLYn59k451zV7at3MCvPuZWjtYZzlo9h5sKd6Q4p7QYOHIiZ8fjjj0cq36VlFxrkNWDGyhkpjsxVJ2XSvRQkXQM8BhxpZsvKHPsCmGBmg8vs/xNwqZkdFsfj3AT0B3qYWXElZa8Frg1/PRGYH/VxUuxggns0pft88dSLWvZA5apyrLz9zYBoNx9JvWzNZWVlsi2Xyc5jIufMlFzGuz9bc5nI+WpiLtubWaNK4oiPmWXMBlwDGFBQzrHdwB/K2X8fsDeOx2gVPsZiYG64vR+xbqRy1fRajcqE88VTL2rZA5WryrHy9nsuU5/LyspkWy6TncdsyGUV9mdlLhM5X03MZSrymBexTZMpyrsMpLhOYLYq3joZ6uUMOV889aKWPVC5qhxL9muVbNmay8rKZFsuUxFbTc9lTcwjZM57Mt66WZvLmtQN9CUw3pLQDZRAfO+b2XdT/Tgu9TyX2cNzmT08l9khFXlM+wDbOHwMnFDO/uOBBdUUQ/QbVbhM57nMHp7L7OG5zA5Jz2NNaqwUAZ1jl8YPF4/rGh5LOTPzN1KW8FxmD89l9vBcZodU5DEjuoEk9Q3/2R24DrgeWA+sN7M3wzINgXnADuAOgvEr9wKNgA5mtq2643bOOedc6mVKY6WiIN40s24x5VoDDwI9CQbJTgd+UXZ8i3POOeeyR0Y0VpxzzjnnKlKTxqzUOJKmS5onaa6ktyV1THdMLn6S6kkaL+mTMJeTY8dOuZpF0u2SPpNUIqlPuuNx0Ug6StI/JH0u6UNJPmuohqrKe9AbK6n1QzM7ycw6AsOAsWmOx1Xdo2Z2XJjLl4HR6Q7IVdl0oDfwVroDcXH5MzDWzNoBvwLGScqGNbNqo7jfg7WmsSKplaQRkmZJ2i7JwtlE5ZU9QtILkjZL2iLp7+F4mbiYWewSxI2rGLoro7pzaWY7zSz2nvOzCW6q6RKUpvflu2a2ONHY3YElM7eSDgM6A08CmNnU8NApKX4ajuS/T6vyHqw1jRXgaOAyoBh4u6JCkhoArwPHAj8BrgSOAWaEM5LiImmcpFUEM5euqELcbn9pyWWMnwMTEqjv/i3duXSpk8zctgbWmNmemKrLw/0u9dL+Pq1py+0n4i0zOxy+WSn33ArKDSL41tzewjs8S/oIWAgMJujOQdIcKn6jdLLwDs9m1j/mMe8HLkjKs6nd0pLLsOxtQDuCafYucWnLpUu5pOa2HN4FVH1SnctK1ZorK2ZWErHoRcDs0hc6rLsUeAe4OGbfyWbWrIKtvP8QxwA9JTVN5Hm49OVS0s3AJcD5ZrY9Oc+mdsuA96VLkSTndgXQUlJ+TL024X6XYsl+n1ZFrWmsxOEEYH45+z8mWNo/EklNJLWI2XUJsA74KrHwXBySkksASTcBPwJ6mtmmJMTm4pO0XLqMU2luzWw98E/gKgBJpWttfVA9IbqIUvY+9cbK/g4l6Jcr6yugSRznaQIUSfqXpHkEl8B+YL6wTXVKSi4ltQIeAA4h6HudK+n95IToIkrW+xJJd4TjyLoAoyWtktQ8CTG6qoma2+uAAZI+B/4I9Pf/TzNOpFxW5T1Ym8asxKO8N0Bc/aNmtgQ4NTnhuAQkI5er4q3jUiLhXAKY2X3AfYmH45Ko0tya2ULg9OoJxyUgSi7jfg/6lZX9FRO0DstqQvktRpe5PJfZw3OZvTy32SNlufTGyv4+Juh3K+t4YEE1x+IS47nMHp7L7OW5zR4py6U3VvZXBHSOXU49XPyma3jM1Ryey+zhucxentvskbJc1qobGUrqG/6zO8FgreuB9cB6M3szLNMQmAfsAO4g6H+7F2gEdDCzbdUdt9uf5zJ7eC6zl+c2e6Q7l7WtsVLRk33TzLrFlGsNPAiUTo+bDvzCzJalOkYXjecye3gus5fnNnukO5e1qrHinHPOuZrHx6w455xzLqN5Y8U555xzGc0bK84555zLaN5Ycc4551xG88aKc8455zKaN1acc845l9G8seKcc865jOaNFeecc85lNG+sOOeccy6jeWPFuRSQdJUki9l2S1os6XeS6qU5tqGxS2eX/d1VTtKgMvmN3b6WlJvAucfGnOuNmP1Dw315CcY+QdJXkupWcLxR+BzGhr/fERPPqkQe27mqSuiP3jlXqUuBVQQ38voP4Lbw3z9PY0yjgUlpfPxs0BH4GuhRzrEdZrYvwfN/QfD3siXB85TnSeAi4AfAi+Uc7ws0CMsBPAFMA+4ETkpBPM5VyhsrzqXWXDNbFP57qqRjgIGS/tPMStIRkJmtImhA1ViS6prZrjSG0BH42Mxmp+j8u1J47leAjUAh5TdWCoEVwBsAZrYaWC1pfYrica5S3g3kXPWaA9QHmsXulHSSpCJJxZJ2SHpH0pnllHlJ0sawzGeSbitTppekWeHxzZLGS2pfpkyVu30ixlnaXXGMpImStklaLukuSTkJnO9ESZMlbQOeC4/9SNKnknZK+pekiyS9Uab7pG9Yf7+rAmHZWXG+BgI6AB/FUy8Vwnxvk/RI6Wtb2WtqZruBvwHnSyr7d9gaOAt42vwuty6DeGPFuepVAGwm+GYLgKSTgZnAocAg4JLw+DRJp4RlTgNmAUcBNwIXAMOAVjHn6QVMBLYBlwM/BU4E/iHp24kGHiXOMl4CXgf6AOOBu4GfJHC+CcCbBF0YD0rqCYwDPg3r/g/wENCuTL3xwBpgcJnn057gg/kvkV6AfzsGOAj4WFJema3KY1XiJakQKALuN7MhZlYSx2v6JJBP8HcS6wpAwFOpjt+5uJiZb775luQNuAowoD1Bd2sT4GpgLzCkTNnpwCdAnZh9ueG+8eHvbwErgQYHeMz3gYVAXsy+I4E9wLCYfUODt375vx/g/JXGGXs+YECZ+v8CpiRwvv8sc76ZwHxAMftODsu+UabsUIJGYsOYfcOAYqB+nLm9LHyM8rYPk/C3MxZYVs7+0tchD/hVmNdrqpKjcP/HwLtl9n0CzDxAXKvS9Z7yrXZvfmXFudT6lOBD5StgDPAXM3uk9KCk+gTf7p8HSkq/oRN8u50GfF9SA6ArMM7Mtpf3IJIaEnxQP2tme0v3m9lS4J3wMaosSpzlVJtY5vf5QOsEzvdSTDy5wHeBF83sm+4KM5sDLC2n7iiCQaM/CuvXI7jK85SZ7Tjws99PJ4JGQ1fg1DJb35gYT5O0XdLcsKtqpqSj43ys8jxIcJWqr5mNjnm8eF/Tp4DTJLUrjRc4Fr+q4jKQN1acS63/IPgQ603wgXF9ePm+1KEE33zvJGjUxG5DCK7INCF4rx5oUGwTgg+lteUc+yJ8nERUGmfZ8SgEDbRYu4DSadtVOV/sc2tG0I2xrpxYvyy7w8zWEHQjXRfuujSMId4uIAgG1y4xs5lm9n6ZbXFMuVOAaWbW0cyOBZYRPN9E/Yjgqsi0MvvjfU3/FyghGFBL+HMX8GwSYnQuqXw2kHOpNd/C2UCSXicYlPlHSS+a2dfAJoIPjJFU/I22OCxzoHEnxQTf9puXc6w5MWNkqqjSOC2+2U1VOV/sgM8NBB/C3yqn6uEEs1nK+hMwPRy7MRh428wWxBFzqY5AlJk6pxA0KkotAlpW4fHK6g5MAV6T1NvMtoX743pNzWy1pGnAFZLuIRi/UmRmxUmI0bmk8saKc9XEzHZJuoXgG/71wB/N7GtJbxOsXzGnog98Sf8g/FApr9siPM8HwKWShlq4zoekNsDpwIgEY48UZ3Wdz8z2SXofuCR8vgYQNkSOpJzGipm9LukTgrEqXYH+8cYt6XCCxl+UmUCnAA+E9QqAMwkGvSbqY6AbweDlSZLON7OtVXxNnyQYpPx7gqtV3gXkMpI3VpyrRmZWJOk94GZJj4QNj5sIBtBOljSGoLujGcEYlFwzuxW4mWAmzCxJDxB0CbUFOppZ6QJzdxKME3lF0p8IZqzcTTCw9IEkhB8lzuo8328IrjC8JGlUWHcoQbdXRR/UfwaGE1yZ+X9rjIQNiqXA3WY2tIL6ncKfeyV1Luf4R2a2XcHqsCcAd0i6kyBXl9m/19xJiJl9IqkbMIOgwdLLzLYS/2v6EsHCczcSdKn5YoEuI/mYFeeq3x0E3RfXwTeDQk8l6Kp5mOADeDjwHYIPHszsPYKrASsJrpK8CtxCzDgWM5tEMKX5EIJ1SP5MMLvjjHDMRkKixFmd5zOzqQRXR44j+ND9L+CXBI2VzRVUez78Odb2X1SuYfjziwM8bMfw5z0EU8ljt5lAnfB4B2AHcJyZtQ/L/6Gy5xQPM/uMYEBtG2CKpMbxvqZhY/l5gvFOz8QOznYukyhmIL1zztVokloRjA35rZndW87xQQSDatuVvcoh6Vrgt0CbimZdxRHHdcAlZtYz/P0ggisXXcxsXiV1xxJ08xxNMKU80aX7ExIugpdLMJutu5m1qqSKc0nnV1acczWSpPqSHpV0iaSzJA0ApgLbCe5/FFv2eEkXEnSLja+gO+Ys4MFEGyqhU4gZhBsOgn2VYIG8KNoQDCCenoRYEvVrglgKKyvoXKr4lRXnXI0kqQ7BNNvOQFOCGwu+DdxuZvPLlH2DYKDxTODHyegWS5Vw7EzpMvhbw+6etJHUgn/PRNttZmm/zYCrfbyx4pxzzrmM5t1AzjnnnMto3lhxzjnnXEbzxopzzjnnMpo3VpxzzjmX0byx4pxzzrmM5o0V55xzzmU0b6w455xzLqP9H1Rd4TbHN0HBAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Specify a nucleus - Germanium\n", "A_Ge = 73\n", "Z_Ge = 32\n", "\n", "weight = 1.0\n", "\n", "#Calculate differential rate\n", "E_R = np.logspace(-3,2,500)\n", "\n", "#Vectorize the differentialRate functions\n", "diffRate_CEvNS = np.vectorize(CEvNS.differentialRate_CEvNS) #Vector exchange\n", "diffRate_mag = np.vectorize(CEvNS.differentialRate_magnetic) #Magnetic dipole\n", "diffRate_scalar = np.vectorize(CEvNS.differentialRate_scalar) #Scalar exchange\n", "\n", "\n", "#Do some plotting\n", "pl.figure(figsize=(8,5))\n", "\n", "pl.fill_between(np.linspace(0.05,1), 1e-2, 1e3, alpha=0.2, facecolor='green')\n", "\n", "\n", "if (\"SM\" in interactions):\n", " pl.loglog(E_R, diffRate_CEvNS(E_R, A_Ge, Z_Ge),'k',label=r\"Standard Model CE$\\nu$NS\")\n", "\n", "if (\"magnetic\" in interactions):\n", " pl.loglog(E_R, diffRate_mag(E_R, A_Ge, Z_Ge, 2.2e-12), label=r\"$\\mu_\\nu = 2.2\\times10^{-12}\\,\\mu_B$\")\n", " pl.loglog(E_R, diffRate_mag(E_R, A_Ge, Z_Ge, 2.9e-11), label=r\"$\\mu_\\nu = 2.9\\times10^{-11}\\,\\mu_B$\")\n", " pl.loglog(E_R, diffRate_mag(E_R, A_Ge, Z_Ge, 2.0e-10), label=r\"$\\mu_\\nu = 3.0\\times10^{-10}\\, \\mu_B$\")\n", "\n", "if (\"scalar\" in interactions):\n", " pl.loglog(E_R, diffRate_scalar(E_R, A_Ge, Z_Ge, (1.0e-6)**2, 1e-3),label=r\"$m_\\phi = 1$ keV; $g_\\phi = 10^{-6}$\")\n", " pl.loglog(E_R, diffRate_scalar(E_R, A_Ge, Z_Ge, (1.0e-6)**2, 1e0),label=r\"$m_\\phi = 1$ MeV; $g_\\phi = 10^{-6}$\")\n", " pl.loglog(E_R, diffRate_scalar(E_R, A_Ge, Z_Ge, (5.0e-2)**2, 1e5), label=r\"$m_\\phi = 100$ GeV; $g_\\phi = 5\\times10^{-2}$\")\n", " \n", "if (\"Zprime\" in interactions or \"vector\" in interactions):\n", " pl.loglog(E_R, diffRate_CEvNS(E_R, A_Ge, Z_Ge, (4e-6)**2, 1e-2), 'm-',label=r\"+ $m_{Z'} = 10$ keV; $g_{Z'} = 4 \\times 10^{-6}$\")\n", " pl.loglog(E_R, diffRate_CEvNS(E_R, A_Ge, Z_Ge, (1e-5)**2, 1), 'r-',label=r\"+ $m_{Z'} = 1$ MeV; $g_{Z'} = 10^{-5}$\")\n", " pl.loglog(E_R, diffRate_CEvNS(E_R, A_Ge, Z_Ge, (1e-1)**2, 1e5), 'b-',label=r\"+ $m_{Z'} = 100$ GeV; $g_{Z'} = 10^{-1}$\")\n", "\n", "\n", "pl.ylim(1e-2, 1e3)\n", "pl.xlim(1e-3, 1e1)\n", "\n", "leg = pl.legend(fontsize=12.0, loc=\"best\", fancybox=True)\n", "leg.get_frame().set_alpha(0.75)\n", "\n", "pl.title(r\"CE$\\nu$NS on a Germanium target (Chooz Reactor $\\nu$s)\", fontsize=12)\n", "pl.xlabel(r\"Recoil energy, $E_R$ [keV]\")\n", "pl.ylabel(r\"$\\mathrm{d}R_{\\nu N}/\\mathrm{d}E_R$ [events/kg/keV/day]\")\n", "pl.tight_layout()\n", "#pl.savefig(\"CEvNS_rate_scalar.pdf\", bbox_inches=\"tight\")\n", "\n", "pl.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }