Skip to content

Instantly share code, notes, and snippets.

@DavidMStraub
Created October 13, 2016 18:12
Show Gist options
  • Save DavidMStraub/9f7d185d68bb057e8765cd6bf99a6fb7 to your computer and use it in GitHub Desktop.
Save DavidMStraub/9f7d185d68bb057e8765cd6bf99a6fb7 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import flavio\n",
"import flavio.statistics.fits\n",
"import flavio.plots\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"from numpy import cos, sin\n",
"%matplotlib inline\n",
"plt.rc('text', usetex=True)\n",
"plt.rcParams['savefig.dpi'] = 100\n",
"plt.rc('font', **{'family': 'serif', 'serif': ['Computer Modern']})"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(-5.8997126213669525e-12-5.6117423595563851e-12j)"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"CSM = flavio.physics.mesonmixing.wilsoncoefficient.cvll_d(\n",
" par=flavio.default_parameters.get_central_all(),\n",
" meson='B0',\n",
" scale = 4.2)[0]\n",
"CSM"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def wc_fct(Cr, Ci):\n",
" return { 'CVLL_bdbd': (Cr + 1j*Ci) * CSM }"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"observables = ['DeltaM_d', 'S_psiK']"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def my_fit(name, observables):\n",
" return flavio.statistics.fits.FastFit(name,\n",
" par_obj = flavio.default_parameters,\n",
" fit_parameters = [],\n",
" nuisance_parameters = flavio.default_parameters.all_parameters,\n",
" observables = observables,\n",
" fit_wc_names = wc_fct.__code__.co_varnames,\n",
" fit_wc_function = wc_fct,\n",
" input_scale = 4.2,\n",
" )"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"fit={}\n",
"for obs in observables:\n",
" fit[obs] = my_fit('CVLL ' + obs, [obs] )\n",
"fit['global'] = my_fit('CVLL Bd mixing', observables)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false,
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 2min 49s, sys: 160 ms, total: 2min 49s\n",
"Wall time: 3min 30s\n"
]
}
],
"source": [
"%%time\n",
"for f in fit.values():\n",
" f.make_measurement(N=200)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAekAAAHqCAYAAAAgWrY5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xl0lOd99//3dc+MZkajZbSLRbvYvCIh4zWNzZa2SZcY\ng+MmT9MmxuAup49/tQPk9Pdr+pwnAYz7+KRpa0B2mqdp6rKIpkmaxEZy0jRxbAMS2GBsQBqJfdEy\no220zMz1++MeCQktCGmkkUbf1zk6mrlnub+A0Geu674WpbVGCCGEENOPEe0ChBBCCDE8CWkhhBBi\nmpKQFkIIIaYpCWkhhBBimpKQFkIIIaYpCWkhhBBimpKQFkIIIaYpCWkhhBBimpKQFkIIIaYpCWkh\nhBBimrJGu4BIUUoVAE8AdUABUK619o3w3JWAG0gFVgPbtNY1U1WrEEIIMRYxE9LAfq11GYBSKhnY\nD6wZ4bmHgBKt9XGlFOHnFk9JlUIIIcQYxUR3t1KqBOjfKSTcgi5TSuWP8JJCrfXxAfdbJq86IYQQ\nYnxiIqSBMqD5pmPNQOFwT9Za1w+4WwSsm5yyhBBCiPGLle5u9zDHvCMcB/qvYW8ESoAUoH6E56UB\nnwo/3jXBOoUQQsQ+B5APvKG1bprIG8VKSHsxB4EN5A4fH5bW2gNsUUq9AFQppfK11q3DPPVTwPci\nVqkQQojZ4vPAv07kDWIlpI8Az9x0LBVzpPcgfS1orfWW8KEDwA7MLvO3hnnveoCd39pD4YKFkap3\nRtr2119l6998I9plRN10/3t44/0r1NS38HdfXEac1TKp53ruued4+eWXJ/Uc09358+f5+te/zq5d\nu6JdyrQgPxNw6tQpvvCFL8AIPbS3IyZCWmtdo5Tq79oO367tu/YcHljmDbeeC4GVA15ehDlw7MgI\nb98FULhgIXfevXQSqp85EpOSZv3fAUz/v4fUHD/H9x7H55jHp+6ZO6nnSk5OprS0dFLPMd0lJSWR\nmJg46/8e+sjPxCATvkQaEyEdtk4p9TzgwWwVDxwMthV4D3hJa12llNqtlHoaUMAqYMUIXd1CzDhz\n3E4Wz03iX9+un/SQFkJMrpgJaa31MeBY+G7FTY+tv+n+qwPulk9yaUJMuUcXZ7LrrbMcP9fCvbkp\n0S5HCDFOsTIFSwgxQFlhKsVZCfw//1JN3bX2aJcjhBgnCWkxZp/+/SeiXcK0MBP+HmwWg//ntxaT\n4orjT/7pPS61+CflPE899dSkvO9M85nPfCbaJUwb8jMRWUprfetnzWJKqVLgaMVPfz6tBwsJMRxv\nRw9f/8FJFIpvb3yAtAR7tEuasc6dO0djY+OQ4+fPnwcgJydnqksSUZSenk5ubu6wj1VXV7Ns2TKA\nZVrr6omcJ2auSQshhnK74vjKp5fwv//jJM9++z1e2/AAiU5btMuacc6dO8eSJUvo7OyMdilimoiP\nj+fUqVMjBnWkSEgLEeMykhy88OklfOMHJ3nmtXd55Y+X43bFRbusGaWxsZHOzk7+5V/+hSVLlkS7\nHBFlffOgGxsbJaSFEBM3PzWezZ+5g53/eYovlb/Dri8tJzPJEe2yZpwlS5bIHGAxpWTgmBCzRF66\ni7/6vTvp6g3ypT3vcL6pI9olCSFuQUJaiFkk2+3kr37vTmwWxdPl73LmSlu0SxLitlRVVVFfXx/R\n9ywvn77LZUhICzHLpCXY+erv3onbZWPja+9y8sKI+9AIMa14PB48Hg/5+fkRfd9Vq1ZN26CWkBZi\nFkpy2tjymTuYm+LkL757lHON0vUtpr8dO3bw9NNPDzleVVXF+vXrh3nF2BQUFODz+fD5fBMpb1JI\nSAsxS8XbrTz3m4tJctj40+8cpqm9O9oliWli/fr1vPXWcJsCDuXxeNi0aROGYfDss8+O+LyamhoM\nw+BTn/oUr7766ojPG0lVVRXFxcVDjm3ZsoX9+/fj8Xhu+z0HWrt2LZs3b57Qe0wGCWkhZrEEh5W/\n/O3F9ARD/Nl3DuPvCUS7JBFlPp+PJ598kq985Stjen5BQQEbN25k48aN1NUN2R24X2VlJUop3njj\njWFbw7eyf/9+nnlm8I7EK1euZPv27axevfq23+9mBQUFEw76ySAhLcQsl55o5y9/azGXvX7+53eP\nEgiGol2SiKKWlhbWrl2Lz+fj4MGDY3rNkSNHeOKJJ0YM6YoKc8+jVatWjbuuyspKkpKSxv36sSgs\nLOTYsWO3fuIUkpAWQpCX7uLPVi/kxAUff3Pwg2iXI6Kkpqamf1DW9u3bx9yaVkpRWFhIbW3tsO9Z\nWlrKoUOHxt3i9Xg8FBUVjeu1t6O0tJTKyspJP8/tkJAWQgBwd46bpx7M49CJK7x7duga1SL2KaX6\nb69duxal1JivHxcUFKCUorW1ddDxlpYWCgoKqKysHHdL2uv1UlhYOK7X+nw+Nm3aRFlZGYZhYBgG\nxcXFLFiwYEitqampw37QiCYJaSFEv5V3ZnHnvGS+dvADWv290S5HTKGqqiqWLh28idBYWtM1NTX9\n4VtQUMCRI0f6H6uoqGDFihVUV1ejlBry/gOfd/OgsIFGu9Y9mrq6OlauXMnWrVs5cuQIR48eJSUl\nhbNnz3LmzJkh3edut5vm5uZxnWuySEgLIfoZSvHlRwvp7g3ytYr3o12OmEIDW9F91q5dS1paGi+9\n9NKIrzty5Eh/F3lhYWF/oHo8nr6doKiqqhq1FV1aWtr/3EjatGkTO3fuJC8vD4CSkhK01kNa0NOZ\nrN0thBgkLcHO5x/Op/xntfz81FUeXZIV7ZJmlK6eIPVTMO88P92FI84SkfeqqKhg7dq1wz62fft2\nNmzYwDPPPDPswK2B4T7wurTH42HFihUAHDp0iDVr1ox4/srKylGvV4+nhdu38Mljjz02pN6RBqBN\npFt9skhICyGGeHhBOkc9zXz9+ycoyUshOV52zRqr+sYO/mj3ryf9PN/Z+CCL5058tLPP5xu2Fd1n\n7dq1bNmyhW3btrFt27Yhj6empvbfLioq4vDhw7z11lv9AQ1mCO/cuXPQ6zweDxUVFRQUFLB7924O\nHDgwYg2FhYV4vbe3Ml51dfWQzVD27Nkz6qInzc3NpKWl3dZ5JpuEtBBiCKUUX3ykgC37jrP9ByfZ\n9rmSaJc0Y+Snu/jOxgen5DyRsGfPHp544olR5wivXbuWHTt2sHXr1kGt0IqKikFBWFhYyPbt23nx\nxRf7j/XNj7733nv7j/l8PtatW9d//XrLli2jLvVZUFAwaku6qalpyLHS0lK2b98+6JwVFRXs379/\nxPc5evTohFYumwwS0kKIYbldcay9L4fv/aqeUxd9LJmXHO2SZgRHnCUiLdyp4PP52Lx5M1u2bLnl\nc5VSbN68mVdeeQUwg7W8vJxVq1axY8cO8vPzKS0t5ZlnnukP3J07d7J3717cbjcvvfQSzz//PAD7\n9u3jvvvu669hLF3MhYWFtLa2DvqQUFNTw969ezlw4AAej4dnn32WZcuW8fTTT/cvsvLSSy/1L1Sy\nf//+Ueda19XVDeoBmA6U1jraNUxrSqlS4GjFT3/OnXcPPzJRiFgVDGn+uuID4qwG3332IQxj5G7R\nWFZdXc2yZcs4evSo7CcdATt37iQlJYWnn36a8vJyfD4fq1atGnH0N8DBgwdpampiw4YNk1KTx+Ph\nxRdf7P8QMppb/Tz0PQ4s01pXT6QuGd0thBiRxVD8j0fyqb3Wzn8euxjtckSMeOaZZ6itreXgwYP4\nfD5qa2tvOTDs8ccfp7p6Qnk3qt27d0/Ltbulu1sIMapFc5J4sDidb77xMZ9YlInbJYPIxMQkJycP\nOwjtVjZv3kx5eXnEW9Mej4fi4uKIb4EZCdKSFkLc0ucezCUU0nz9P05EuxQxi+Xn51NUVER9fX1E\n37eqqmpcm35MBQlpIcQtuePjWH9/Lv/98XWOeqbXikxidlmxYkXEW7zTNaBBQloIMUafXJJJcVYC\n//v7H9ATkJ2yhJgKEtJCiDExlOKPf6OQa63dfPeX41tLWQhxeySkhRBjNj81nt+8Zw7/9F91NEzB\n0pdCzHYS0kKI2/L7y+aTmhDHXx84Tigk6ywIMZkkpIUQtyXOavBHnyjko8tt/LBG5k4LMZkkpIUQ\nt+3O+ck8sjCDb/70I5rauqNdjphFqqqqIj4FC6C8vDzi7xkJEtJCiHF56sFcLIbiawdl32kxNfq2\nn5yMRUdWrVo1LYNaQloIMS4JDhtfeDifw3XN/OKja9EuR8wCO3bsGHZOc1VV1YR3ryooKMDn8+Hz\n+Sb0PpEmy4IKIcbt/qI03j7TyDf+4wTL8n8Dl0N+pcxkHo+HAwcOkJKSgtYan8/H2rVrqa6uZu3a\ntVGtraqqiuLi4iHHDh06hNfrHXWrzbFau3YtmzdvZteuXRN+r0iR/1FCiHHr23d6677j7HrrDH/5\n20uiXZIYp5qaGvbs2TNoFyifz0dBQcGkbmwxVvv37x+0TzXAypUrWblyJRUVFRw9enTC5+jb0nI6\nke5uIcSEpCXaefy+HA68e44PL06vrkIxdhs2bBgSgsnJyaxevXpabDxRWVk56l7QkVJYWMixY8cm\n/TxjJSEthJiw1Xdlk5MWz98c/ICgzJ2ekaqrq9F66L/dxo0bo1DNYB6Ph6Kioik5V2lpKZWVlVNy\nrrGQkBZCTJjFUHzxEwU0NHbw70fOR7scMQ6rV6+mtLSUnTt3DuryXbFiRRSrMnm9XgoLC8f1Wp/P\nx6ZNmygrK8MwDAzDoLi4mAULFtDa2jrk+ampqdTW1k605IiRkBZCRERxViKfXJzJPx46TVO7zJ2e\nad544w02bdrEvn37KCoqIi0tbdRu34qKiiHH+vZlfvLJJwfNZS4rK+PZZ58dNhTHoq5ufGvF19XV\nsXLlSrZu3cqRI0c4evQoKSkpnD17ljNnzgzbfe52u2lunj47vUlICyEiZv39ORiG4huy7/SM9Pzz\nz3P48GG8Xi/r16/niSeeGPG5SqkhxwoKCgB48cUX+69jl5eXc+DAAV555ZUpuaY80KZNm9i5cyd5\neXkAlJSUoLUe94eFaJCQFkJETILDxrrlufzqdCPHGlqiXY4YA4/HMyS0kpKS2L59Oy0tw/8bVlVV\nsWrVqiHHfT4fLS0t5OXl4fV6efXVV9mwYcOEB56Np3Xbt/DJY489Nui4UmrUDwsT6VqfDDIFSwgR\nUb+xOIOfn7rK9h+c5F//9GEMY2iLK5Z1B7q40H5h0s8zP2E+dqtjwu+ze/duvvrVrw45Xltb279A\niMfjobKykqNHj7Jr1y58Ph9JSUls2bKF5cuX8/jjjwPmCOxVq1ZRU1PDunXrIjZ1q7CwEK/Xe1uv\nqa6uprS0dNCxPXv23HLRk+bmZtLS0m67xskiIS2EiChDKT73QB7bfvghb564zG/eMzfaJU2pC+0X\neO7nfzHp53n50W9S5C6+9RNvobKykuLi4kErefl8Pnbs2NG/TKbH42H9+vUcOHAAoH8U+MaNG6mp\nqel/Xd/CIikpKaxevZrdu3fzwgsvDDlnRUUFLS0taK2pra1l06ZNo7a2CwoKRm1JNzU1DTlWWlrK\n9u3bB/2ZKioq2L9//4jvA3D06NEJr14WSRLSQoiIWzw3idL8FP7ujY95dEkWDpsl2iVNmfkJ83n5\n0W9OyXki4XOf+xxlZWXs3LmTlJQUWlpaUEpRXl7e3y28YsUKysvLWb16NXDjenRdXd2gbu8jR47w\n2muvkZ+fz1e+8hXKysqGhLTP5+Pw4cP9AfrSSy+NqTu8sLCQ1tbWQV3VNTU17N27lwMHDuDxeHj2\n2WdZtmwZTz/9NAUFBWzcuJGXXnqpf5GS/fv33/K6eF1d3bQY0d5HQnqMTp73kZbjJzPZgTHMgAkh\nxGBP3p/LV/e/z753GvjDT0yfa3yTzW51RKSFO1Wef/55AJYuXTrq8/pamDU1NZSWlvavcz0w9Orq\n6rj33nsBs/W7bNkyysvL2bBhQ/9zkpOTqaurY+vWraxfv77//Lfy5JNPsnfv3kHvVVJSQklJyaAW\n80DDrfM9mqmcjz1WMnBsjPa9d47Ne4+z6duH2fmfpzjW0EJomIn/QghTttvJijuy+PZ/1dHS0RPt\ncsQE9U3P2rNnD9XV1Rw9erT/WrTH42HTpk0opQZN23K73WzZsmXIVK59+/axceNGdu/ezVtvvTWm\n8z/++OOTvjzp7t272bx586Se43ZJSI/Rt754H3/3h2VseKyYzp4AL//0Y7667zhvfXiV7t5gtMsT\nYlr6vdJ5KAXf/q/psziEGJ+lS5eya9cuVq9ezeOPPz6oS7igoIBdu3bR1NQ0qEW+b9++Icf6rnPn\n5+ezbt262xoQtnnz5knbTrJvjvd0WAJ1IOnuHqMkp43SojSWF6Xx+Yfz+eC8lz1vneWff+mh4vB5\nfrdkHivvzMJqkc89QvRJdNr49NK5HDx8nvX355KT5op2SWKC3G73uF/r8Xioq6vj1Vdf7b/2Pdbu\nbjCDvaioiPr6+oiHaVVV1W13j08FCelxUEpxT24Kf/9H93GppZNvvfExr7/TwM9OXeUPHsznntzx\n/xALEWvW3JVN1cmr7PzRh/zdF++LdjligiYyqKqgoIBt27ZF7fyjmY4BDdLdPWFzU+LZ9rkS/nnT\nQyTHx/G3P/mI//OTj7js9Ue7NCGmBbvNwmfL5vNeXTMfX545Kz0JMR1ISEfIguxE/umZB9j25FIu\ntvj56r7jfOcXdXhlwIwQPLIwgzluB//nx6eiXYoQM4qEdAQppXjsjiwq/uIT/OnqRbxX18wL/3aM\nA++dp7M7EO3yhIgai6H47LL5HD/n5YPzt7dylBCzmYT0JLDbLHz+4Xy+/9xv8NSDebzxwWVe+Ldj\n/PjYJRkJLmat+4rSmJ8az8s/kda0EGMlIT2JEp02nl21kIq/+AT3FaRy4PB5Xnj9GG9+cJmeQCja\n5QkxpQylWHvffD682MpRz9BlHIUQQ0lIT4GMJAfbPlfCvj9/hLtz3PzrrxvY/G/H+NmHVwkEJazF\n7FGSl0JBhouXf/JR//rPQoiRSUhPoXmp8bz0+VL+7c8eYUF2Iv/3vz1s2Xuc//74GsGQ/MISsU8p\nxeNlOZy92s67tdKaFuJWJKSjIC/dxTf/sIx/+ZOHyEt38erP6/jqvuP8+myjLDUqYt7dOckUZyXw\nzZ9Ka1qIW5GQjqKirET+/o/u4zsbHyQr2cGuqrP87++fpO5ae7RLE2LS9LWmPdc7+NXp69EuR4hp\nLWZCWilVoJR6QSm1Vin1vFIqeZTnrgw/9wWl1D6lVMFU1nqzxXOT2P3l+3nlj++jJxDib/79BK/9\nvBZfp8yxFrHpjnlJLMhO5B8rT0treprz+Xx4PJ4xP7+qqopNmzaxZcuWKT1vrIqZkAb2a613aq0r\ngHJg2J29w+FdGn7uTmAvcGgK6xxRSX4qr//Zw7zw6SUcrW9h897jvPnBZUJyvVrEGKUUv1c6j7pr\nHXJtehqrqalh5cqVrFmzZsyvWblyJUVFRVRVVU3peWNVTIS0UqoE6E8yrbUPKFNK5Q/z9DJg4Oaj\nlUDhCM+dclaLwdrluRz8n5/ggeI0vvd2A1//wUkutnRGuzQhIuqu+ckUZLj4x0Ono12KGEFJScm4\ndp0qLS2NynljUUyENGbwNt90rBkYstO81roKWDbg0H3mYV0/adWNQ3J8HNueLGHXl5bT3hXg/zvw\nAf9RfUGmbImYoZTid0vncfpKG8caWqJdjhDTUqyE9HDbTnlHOI7WeuAO5F8BnpmMoiJhaV4Ke//8\nEf7goXy+f+QCXzt4gnNNHdEuS4iIWJqXwvzUeP7hzY+jXYoQ01KsbFXpBVJvOuYOHx+RUmoDsE9r\n/dpkFRYJdpuFP1m9kBV3ZvNX+4/xv/79BOvvz2XVXdkYSkW7PCHGzVCK3ymZyytVZzl10ceSeSOO\n9xSTxOPxsGPHDsrKyjh79ixFRUWkpqZSWVnJK6+8MuLrdu7cSVFREVprPB7PkH2htdYcPHgQrTXN\nzc14vV5eeOGFQeetrq4G4MiRI6xatYqVK1dOzh9yBouVkD7C0NZwKlA30guUUiuBJq31wbGc4Lnn\nniM5efAvkKeeeoqnnnrqNksdv8Vzk/jenzzMP1ae5ntvN/D+eS8bHi0iOT5uymoQItKWF6Zx8MgF\n/v7Nj/mHP14e7XImLOT3Ezh7dtLPYy0uxnA6J/w+69at48UXX2TFihX4fD6WLVvG2bNnWbZs2Yiv\nWb9+PZs2berf29nj8bBmzRrefPPN/ufU1NRQWlpKfn4+AOXl5WzatIldu3YBsHHjRtasWcPzzz/P\n2rVrKS4uprKysv/5M8Xrr7/O66+/PuiYz+eL2PvHREhrrWuUUv1d2+HbtX3XmcMDy7xaa0/4fmn4\ndQfD9zcAe7XWI252+/LLL094MEQk2G0WnvutJTy4IIO/PvA+f3XgAzY8WsQ9ucP27Asx7RmG4rfv\nnct3flFH/fV28jMSol3ShATOnuX6b/72pJ8n46c/Ju7uuyf8Ph6Ph9RUsyMyOTmZuro66uvrRwxL\nj8dDVVUV+/bt6z9WUFBAc3Mzb731Vn9wFxYWDnqPDRs2sHHjRrZs2UJ+fj579uzpP2/f86urq2dc\nSA/XWKuurh71Q87tiImQDlunlHoe8GAOJFs34LGtwHvAS+E50UcArcyuYgW0aK1n1FDCB4rTef3P\nHmbLvx3jb3/yEb9bOo/PLpuPYUj3t5h5Hl6Yzr8fOc8/HDrNzj+I/ofhibAWF5Px0x9PyXkiYf36\n9Rw5coSlS5dSXV1NUVHRqEFZWVlJYeGQMbkUFhZy6NCh/pAejtvt7g/i/Px8PB4P+/btIyUlhebm\nZpqbbx7/K2ImpMODwfoGhFXc9Nj6Abc9xMiAudQEO7u+tJzv/tLDrrfOcPZKG8+uWkCS0xbt0oS4\nLTaLwafumcOB987T2NZNeqI92iWNm+F0RqSFO1VWrVqFUory8nJ8Ph9Hjx6dkvO++OKLHD16lL17\n9wL0fxeDxURYzWaGofjibxTy918s42KLn//3wPucvjxir70Q09ajizOxWQz2vdMQ7VJmlbq6Oh5/\n/HE2bNjA888/T1JS0qjPLysro65u6HCf6upqli8ffUxB3zVvn8/Hli1bBgWz12uO8z127NhIL5+V\nJKRjxLKCNL73pw+Tlexg+49OUXniiiy3KGaUeLuVR5dkUnH4HP6eQLTLmTXq6upGXX5Ta01T041V\n4UpKSli1ahUHD94Yc1tdXU1KSgqf/exn+481NzfT2nqjwVBeXs4zzzxDXl4ezc3NKKUGPd7S0kJL\nS0v/B4CbzztbSUjHkPREO69teIAnlufw3V/V80+/8NAri5+IGWTVXVl09gT5yfHL0S5l1njiiSco\nKirCYrFgsVhYsGABzz77LGCO0N6yZQs+n4+tW7f2h+revXs5fPgwr776KuXl5ezfv5/Dhw8Pet/y\n8nIqKyupqKhg586d+Hy+/ildBQUF7Nixg82bN3Pw4EEOHjzI/v37qayspLm5ecTzzkYxc01amKwW\ng+d+awkLspPYHl5O9M/XLMQt07TEDJCR6GBZfirf/WUdvy8DISddTU0N+/fvx+PxkJeXR2trK83N\nzWzfvp2tW7eybdu2QdOqBtq2bduI7zuW+c43z6sGeOONN/pvj3Te2UZa0jHqMyXz2PXl5TS2dfO1\ngyfwXJftL8XM8Km7s7ns7eKd2sZolxLzjhw5QllZGXl5eQAkJSWRn5/P6tWro1yZ6CMhHcPumu/m\nu88+RKorjm/84EOOemR6g5j+FmQnkp/u4rWf10a7lJi3YcMGtNZs2bKFl156ifLycrZs2YLH4xm1\npSymjnR3x7iMJAevPfMAL/xrNd968zR/8FAea+6eE+2yhBiRUoo1d2ez52e1eK63UzDDFzeZ7jZs\n2BDtEsQopCU9CzhsFr75P8r4g4fz+d7bDXzv7XrZo1pMa8uL0kh22jjw7rlolyJEVElIzxKGofjz\nNYt4/tNLOHTiCv9QeYaegIz8FtOTzWLw2B2Z/KjmIm3+3miXI0TUSEjPMk8sz+XFp0p4/7yXv/3x\nKTq7ZT6qmJ4euyOLQEjz4+OXol2KEFEjIT0LfWJRJt/6YhnnmjrZ/qNT0lIR05I7Po5l+Sn826/r\nZWEeMWtJSM9S9+amsPvLy2np6OEbP/yQlo6eaJckxBCP3ZHFZW8Xx8+NujW8EDFLRnfPYguykyh/\n+n6e/fZ7fP0HJ9n8mSVkJDqiXZYQ/RbPTSIrycE//3cdS/Mis/XfRJw6dSraJYhpYCp/DiSkZ7m8\ndBevbniATd9+l20/+JCtv3MHGUkS1GJ6MJTi0SWZHDxyHl9nD8lRWjkvPT2d+Ph4vvCFL0Tl/GL6\niY+PJz09fdLPIyEtmJvi5NWnH2DDq++y7YcS1GJ6eWRRBgcOn+cnxy/xuQfzo1JDbm4up06dorFx\n6Cpo58+fByAnJ2eqyxJRlJ6eTm5u7qSfR0JaAJCZ7KD86ft5+tV32P6jU3z1d+8gLWHm7ukrYkeS\n00Zpfgr73j3Hkw/koVR01vPOzc0d9pdy39aOxcXFU12SmAVk4Jjol5nsYM+X7wet2fGjU3hlMJmY\nJh5dksmlFj8nLviiXYoQU0pCWgyS7Xay+8v30xMIseM/ZXqWmB7umJdMekIc3/3vumiXIsSUkpAW\nQ8xLjWfXl5bT5u/lb3/yEf6eYLRLErOcoRS/sTiTd2qb6OiSBXjE7CEhLYaVl+7iW1+8jyu+Lv7u\nzY8JBGUJURFdn1iUQW8wxKETl6NdihBTRkJajGjx3CT+z+dLOX25jVd/XktIVn0SUZSaYOfu+W4O\nvCebbojZQ0JajKokP5W/eeIe3jnbxH7ZkUhE2UML0zl7tZ0LzZ3RLkWIKSEhLW5p5Z3Z/MVvLuLH\nxy9z6MSVaJcjZrHSvBQcNoM33pdNN8TsICEtxuRzD+bzuQfz+N7b9RxraIl2OWKWstsslBWk8oPq\ni7LphpgVJKTFmP35mkWU5KXwj5VnON8k3Y0iOh5akMFVXxcfXpQ50yL2SUiLMbMYir/9fClZyQ6+\n+cbHMocUx6qHAAAgAElEQVRaRMWSuUkkO20c+kAuvYjYJyEtboszzso3/7CMrt4g/1B5RqZmiSln\nGIr7i9L4yfuXCIaky1vENglpcdvmuJ28+FQJp6+08fqvG6JdjpiFHliQjq+zl+r65miXIsSkkpAW\n41KSn8pf/vYSKk9e5RcfXYt2OWKWKcxwkZlk580PZGETEdskpMW4PX5fDp9cnMk//9LDuaaOaJcj\nZhGlFA8Up1N14gq9AbnkImKXhLSYkP/1xD3McTv5h0NnZI1vMaXuL0qjsyfIu7VD93gWIlZISIsJ\nsdss/O3nS/F29vJPv6iTuatiysxPjWduipOK985HuxQhJo2EtJiwnDQXf/X7d/FubRM/OyXXp8XU\nub8ojZqGZrp6pRdHxCYJaRERq+7KZuUdWfzr2/U0NMr1aTE17i9Ko6s3xLtnpctbxCYJaRExf73W\nvD79StUZuqVlI6bAHLeTuSlOflh9IdqlCDEpJKRFxMRZDXY8VUJjWzd7ZccsMUWW5adQ0+CVhXVE\nTJKQFhFVkJHAX3xqMVUnr3L8nGzEISZfSX4qHd0B3j/vjXYpQkSchLSIuLXLc7gnx81rP6+jVdb3\nFpOsIMOFO97GL2TQoohBEtIi4pRSbHtyKSGtZVqWmHSGUpTkp1B18or8rImYIyEtJkVaop2tv3cn\n1fUtvFfbFO1yRIwryUvhels3nusys0DEFglpMWlW3JHNfYWpfPdX9dLtLSbVkrnJ2K0Gv/xYurxF\nbJGQFpPqa2vvIaThe2/XR7sUEcPirAZ3zk+WDTdEzJGQFpMqLcHO859ewjtnm6ipl9HeYvIszUuh\n9mo7LR090S5FiIiRkBaT7jfvmcPdOcn88688ssiJmDT35rrRwK/PXI92KUJEjIS0mHRKKb72+D20\n+Xv5QfXFaJcjYpQ7Po6CDBdvvC9d3iJ2SEiLKTEvNZ4vfqKQn7x/mctef7TLETHqnhw3Jy7I6mMi\ndkhIiynzhUcKSE2I47u/qpf5rGJS3JvrpqM7yMmLvmiXIkRESEiLKeOwWdjyO3dy8oKPw3XN0S5H\nxKCCzAQSHVbePi3XpUVskJAWU+rhhRkszXWz950GegLSJSkiy1CKu+a7+S9ZIlTECAlpMeW++nt3\n0dLZyxsyp1VMgjvmJdHQ2IGvU6ZiiZlPQlpMudx0F+uW5/Kjmot45RepiLA75iWjgWqZly9igIS0\niIovPVqExTA4ePhCtEsRMSY90U5mkp0jdbJmvJj5JKRFVCQ5bWxcUcwvPrrGuSbZFEFE1h3zknlb\nFjURMUBCWkTN4/flkJFk59+PSGtaRNaSuUlc9nbR1NYd7VKEmBAJaRE1VovBs6sWUl3fQu219miX\nI2LI4jlJAFQ3yFQ/MbNJSIuoWnP3HOamOKk4fD7apYgY4nbFMcftoNojg8fEzCYhLaLKYij+fM0i\nTl7wceqSrBIlImfxnCTeOSvXpcXMJiEtou7RJZnkp7tkpLeIqEV916Xb5bq0mLkkpEXUKaX4k9UL\nOX2ljY8utUa7HBEjFmUnAnD8nHR5i5krZkJaKVWglHpBKbVWKfW8Uip5DK9ZOxW1iVt7ZGEGOanx\n/LBGtrIUkZGaYCc90c7xBglpMXPFTEgD+7XWO7XWFUA5sH+kJ4aDfAOwXymVNGUVihEZhmLjygWc\nuOCTkd4iYhZmJ/LOWVnURMxcMRHSSqkSoH/vQ621DyhTSuUP93ytdYXWunzga0T0PXZHFnPcDn5Q\nLa1pERmL5iRyrqmDjq5AtEsRYlxiIqSBMuDmCZHNQOEtXqcmpxwxHhZD8eVHiznW0MKlFn+0yxEx\nYGF2ElrDiYveaJcixLjESki7hznmHeG4mMbW3D2HZKeNN0/IDlli4rLdDlx2CyfOS0iLmSlWQtoL\npN50zB0+LmaQOKvBkw/k8avTjbR39Ua7HDHDGUpRnJXIO2cao12KEONijXYBEXIEeOamY6lA3S1e\nN+Zr0s899xzJyYMHjD/11FM89dRTY30LMUafvS+Hf/pFLT/78Bq/Uzov2uWIGa44K5GfHL9EKKQx\nDLnCJSLr9ddf5/XXXx90zOeL3MJMMRHSWusapVR/13b4dq3Wuj58vwTwaq09N710zP9jX375ZUpL\nSyNRrriFFFccDy3I4NDJK/zWvXOwWmKlw0dEQ3FWAp09QeobOyjMTIh2OSLGDNdYq66uZtmyZRF5\n/1j67bcuPD96LbAZWDfgsa1A/5xopdRKpdQLmC3prUqpFVNbqriVjSsX4Ovs5Wi9zHEVE5Of7gKQ\nZWfFjBQTLWkArfUx4Fj4bsVNj62/6X4VUAXsnJrqxO0qzExg8ZxE3jp5hfuL0qJdjpjB4u1WspId\nfHSplU8vlcsnYmYZd0grpZZiTnHqG7BVB9T1dTELMVFfeKSQv9p/nAvNncxPjY92OWIGK8hwUVMv\n21aKmWfMIR1emWsj8CRQAHgw5yL3jaB2A6lKqQLMgVz7gX1aa1mMWYzLJxdnkuy08daHV/nDRwqi\nXY6YwfLTXVTXtxAMaSwyeEzMIGO6Jh2+flsFtADrtNZpWusyrfUarfX68Nea8LE0YBPmoKy3lFJP\nT175IpbZrAaP35fDr0430t0bjHY5YgbLz3DREwjR0NgR7VKEuC2jhrRSKlkptQ84qrW+T2v96jAj\npIfQWnu01uVa6zLAo5R6RdbIFuPxmZJ5dPUGZQCZmJDcNHPw2Jkr0rEnZpZbtaTXAxu01m+N9wTh\nQVpbwu8lxG2ZlxrPwuxEfnn6erRLETOYy24lPSGO01faol2KELdl1JAOt4YnPG9Ba+3TWr860fcR\ns9MTy3P58IKP5vbuaJciZrDcdBfHZNtKMcPE0jxpEaNW3JmF1aJ4W5Z2FBOQmxbPucYOtJbN78TM\nMekhHb6uvUEWDBHjleCwUZqfyq/PSkiL8ctJc9HWFaCpvSfapQgxZpMe0uGu7nKgSCl1ZrLPJ2LT\n2uW5XGj2c6G5M9qliBkqJzzXvvaaXJcWM0dEQ3q06VbhoL7lyHAhhvNgcTouu4V3pDUtxikj0U6c\n1aD2anu0SxFizCLdkl59i8erI3w+MUvYrAYr78zmndomuaYoxsUwFPNSnNRelZa0mDkiHdLrlFJB\npdRhpdS2Ya5Dn43w+cQs8ql75nC9tZu6a9ISEuMzPzWekxdkow0xc0Q6pDcDxcAeoAg4MCC0nwci\ns3eXmJWW5qWS5LRxxCNrMIvxmZfi5JLXTygkvTFiZohoSGutdw5YbWy91joVuA/YB6wBnonk+cTs\nYjEUn1ySydH6FunyFuMyN8VJTyDE1dauaJcixJiMGNK3WsZzrMt8aq2rw+G9BnPlMSHG7ZOLM7nq\n6+KyV37Jits3N8UJQP11uWQiZobRWtK3avWOp1W8ZxyvEaJfWWEadqvBUdl2UIxDWoI5wttzXTba\nEDPDaFtVvqiUenKUxwuBl251AqVUMua63bUTWQNcCACHzcLdOW5q6lv4nZJ50S5HzDCGUsxxO6Ql\nLWaM0UK6HHNrysPDPKYYY0s6vPZ3eXjVsd1a6wW3X6YQN6y5ew7bfniS9q4ACY4xb4kuBADZyU5O\nX5bdsMTMMOJvOK31xnAreBWggUqtdf9PtlJqyEr1SqmnR9pIQ2tdrpRaF4GaxSx3f3EaWsOHF30s\nL0qLdjlihslOdvCxhLSYIW61C5ZPa12htT4I3KeUelwptTT8WNUwL5HFTMSky0p2Mtft5ITMdxXj\nkO124u3spaM7EO1ShLilMU/B0lpXhcPaq5R6RSn1+DBPk8VMxJT4xOIMPrjglalY4rZlJzsAON8k\n68CL6W/MIa2USgovSFKJOfd5OLKYiZgS9xel09zew1WfTMUStyerL6SbZYS3mP5uGdJKqaVKqX1A\nPWbwrtZal4Vb1YPIYiZiqtybm4Kh4NQlubYobo/LbsVlt3Cp2R/tUoS4pdEWM3k6vLXkfuBNrXWq\n1vpZrbUn/Phw3d1DyGImYjK4HFYKMhIkpMW4pCc6uNgi3d1i+htt/sqe8NdRQN20DaUCvgIMaU3f\ngixmIiLm4YUZ/MfRC2itUUpFuxwxg2Qk2mXLSjEjjLqYidZ6xJavUupWI7mHCM+ZFiIilhWm8s+/\n9HDJ62deSny0yxEzSEaiXVatEzPCaNekd9/itZtvPjDMaG4hJs09OW4MBacvy/7A4vakJdppbu+R\n3bDEtDfaYiae0V44wuM7GHnktxAR5YyzkpMWz9mr7Tx2R1a0y5m9tIZgEEIhCAYgGIJQ0FwCyQhf\nhlDh9oDR932Y4wpAgXXyV5FLdcURCGlaOntIS7BP+vmEGK8R/zcopV7AXJ+7Vmv9Unj1sf2Y06gq\ngQ0DVyALawlPtTqgta6fpJqF6FdWkMbbZ65Hu4yZqbcXenqgtwfV02Pe7ulB9d64bfF3QKDXDN5g\nEBUKmoHcF8TBIES6NaoAmw2sNrTNZoa2zUbIEW/eH/jYgC/tjAeH48YHgVGkJsQBcL21S0JaTGuj\nfWStA+q01hXh+5XAEa31GqVUKbA1/NUvPIIbpdRKpVTpcNO0hIiku3LcVBw+T3tXLwkOW7TLib5Q\nCLq6UJ2d0NWJ8vuhsxPV5cdobxsUwAyXrQqIi4O4OLTN/I7dARYLGAYhi9UMQYsBfbcNwwzS8G1t\nsZqPg9nK1jp8LvO2GubYoK9AwPwAEehF9faGbwcw2luhN3wsED4+3J/BbkfHOwklJqPjXeCMR8fH\nm7fj4yEurj+Yr/q6WDw3eTL+JYSIiNFCWveFbLgVvQxYGX6gWik14sIkWusqpVRyuFW9Z5gWtxAR\ncXeOG4Daa+3cm5sS5WommdbQ5Ud1dIDf3x/ElrZW87i/C7puWtzFagWnw2xlJiSgbXHocAgTF75t\ns0FfINsm/4PO7bS7b/ncQMAM7J4e88OJ32/+XXT6MfydqOYm8HeaLf8+Vitp8U7KLl6h870AvXFL\nMFLcGMnJKItlHH8iISbPWC/+rMLs9h4YtqP+/wmP5H5JKbVWKdUi21SKyTAvxYnLbqG+sSM2QjoY\nhM4OVEcHqr0dOtqxtHpRHe3QeVPY2Kxm+DqdkOQmlOUEhwMdHw8Op3l8CkI3qqxW88vhhKTk/l9K\nQ345dXdBRyfK32n2LPg7cdkb6alvoDt0zXyOUhhJSRgZmRgpKWZwp6RgxMvMARE9o4V00YDbGzG7\nu0ellMq/+Vq01rpCKZWvlNoGbJNWtYgkpRS5aS7ONc6ghSlCIehoR7W2DgjiFrOF3DlgFSxDgcuF\ndiWgM7PM7lqXK9xt6wRrjAdwJNkdYHegSQXMEL90Lp7jWQn8zqcKCflaCXlbCHl9aK+X3oYGdG8P\nACouDsPtxsjMuhHcbre0usWUGC2kK5VSR4BkzCtVT4B5vRlzIZP9w7xmh1JqF2bAL8MceFYGuMPv\nsQoZ/S0ibGleCpUnrkS7jKGCQWhvR7X5UK2tWLzN0BYO5r7BVjYr2pUALhehlFSId/Xfx+kEWaRl\n0ridVs43d6JsNizpaVjSB297GuroIOT1Emrxon1egg319J4IT/ezWLCkpWGZP59gewdGagz04ohp\nabQpWDVAmVKq4KbpVs2MvLznuvBXHea2lDWYq4zVhd9PiIhbOCeJ13/dQEd3AJd98qfvDBEMmuHb\n1oZq9WFpaUK1tUFH+41+V7sdnZQI6RmECorQCYmQlGi28ERUuONtfHRl5FXHDJcLw+WCefP6j+lA\ngJDXR6jxOsHrjfSePEnXuXNgGPjr6rDk5WGZMwcjNVVWwRMRccvfaDfPh75F2I66SpkQk2FhdiJg\nbj24eG7S5J6suxvl84a/fBhN18xA7gtjpwOdmITOykInLkAnJkJikjkoS0wrbqcNr//29pRWVmt/\nq9u2GLTWOI4fJ9jYBHFx9B47Rs/hw6i4OCzZWVjy8s3Qdrkm6U8hYt0tQzq8itiRsVxLloAW0ZCb\n5sJQcNnrj2xId3SYYextwdJ0HdXqu3HN2GpBJyZBWjqhwmJ0UpIZxrE+UCuGJDmttHcFCIY0FmN8\nrV6lFEZiIkZiIo6iInQoRKipieDlKwSvXCHwy1+C1hhJSVhy87Dm5mBkZkorW4zZWPoGFeBVStVi\nDh47JPOfxXRisxpkJTu41DLOrQe1hvY2VEsLqqUZS3MjyueF3nAryx6HTnYTmp9jjiBOdkNCglwv\nnuGSHFY04PP3kuqKTE+HMgwsGRlYMjLgnrvR3T0Er14heOUqgbpaek98gHI6sRYVY83LxcjIkMAW\noxpLd3eVUmonUIJ5nVkNN4pbiGia63Zy0TvGkPZ3opqbzRby9asob8uNQHa50G43oeKF4Hajk5PN\n6T0i5iQ5zV6Pls7IhfTNlD0Oa24u1txctNZmK7vhHIHas2Zgx8ebgZ2fhyU9fVJqEDPbWLq7Hwfe\n1VoP2VBjwHNewFz85KVIFifEWN2Vk8IPqi8MfaCnB9UyIJBbWm4s+OFwoFNSCBUvRKekgDtFrh3P\nIkkO89dfc0cPRRmTf81YKYUlPR1Lejq20hJCjY0EGs4ROHOa3g/ex0hOxrZ4CdbCApRDBhQK01i6\nu1drrZ8d7Qla650ASqntwDdkLrSYannp8Xjbu+i53oij1Wt2W1+/Cu3h0bs2K9qdQig3D1JS0O4U\nc4rTbBIImF37ofDa2zpkTgW76b4KhTfLCIVuLNXZt/ynYaD7bivDnMttWAZ/VwOWDZ3GXbmJ4ZD2\n3ebgsUhQSvV3i+tlpYTC3eHdh9+j+/B7WOfPw7r4Dixz50h3+AwR0iG6Al34A51c6bgcsfeN6HwV\nrfWW8FKg0qIWky7k9xO6fp1QYxP5759mxekP6P3xGRIS7OhktznCetFis5XsipFryIEAdHdDb4+5\njnWP+d386ulf69rwd0BPLyrQCz3hta5DoVu/fyQpFV561A5xcYTiXRBnR9vt4TXB7QOWKA0fm4Id\nsPo4bQaGglZ/75SdczhKKSxzsrHMySauu5tAfT2BOg9db76BcrmwLV6MbcEC1Gz7UDlFekO9+Hs7\n6Qx04g/46ewNfw+Ej/X68Ydvd4af5+1uxh/ooiv85Q/66Q50o8PTPLy1vojVN5b/Eam3+Z4VSqnH\nZXCZiCQdDBJqbibU2Ejw0kWCTc3ojg4AlMtFaloSZzIKaLjrLpLuzDU3hJhJAr1mN3xXF6qry1zG\n0t9lbozR0X7jWO8IgdK3E1RcHFhthJzxkJB8Y6eovgA0jP4Wr765BTzc974PNlqHW9ah8A5YfS3v\nwV+qr/UdDIZ32eo2P0h0d5sfHNp8qO7wBh/BYVqwFqu5rGmCi2BSSnjFNVd4cZf4Me1wNVZKKRLs\nVlq7pr4lPRJlt2NbtAjbokUEG5sIeOroPX6cnmPHsBbkE3fPvRjJs3tDEK21GayBcHiGQ7UvWPtv\n9x83v7d0N/eHalfQT1egm65AFwE98r+/QmG32HFanTisDhxWB06rA4fFidueMuh+/22rg3P2c/w3\nv4rIn3dMo7uVUklj7cLWWnuUUmsnWJeY5XRXF8Hr1wlevUbw0kVCLS1mEFgsWFJTsebmYqSlYqSl\nYcTH49CayzW/5KLFxd3TMaADvebqYx0d5vfOToxWL6rLb4Zz4KZfFBYrOOxoh4OQKxFS09B2h3kd\n3eG4sRmGLfx9mvQSDLegf9+xIe34QCC8I1e3OXagpwe6u1AdHRg+L9aL58z1yvt6AAwDnE60y0Uo\nyY1OSDCnwSUmjfvShctuoTUK3d1j0TcfW99zrznQ7PQZOs/WYpk3j7jS0hk30CykQ+Hw7KAz0ElH\nb2f/7b4W6sBWrD/Qibe7JRyqZrh2B7vxB/yE9Oi9Qg6LA7vVjsPi6A9Xh8VBmjOt/7bT6sBuuRGs\nNz/XYXUQZ4nDULf/wVAnRW771rGE9G6gHHjyNt437dZPEeKGkN9P6No1glevErxwgZDP7C5SLheW\n9HSs+fkY6WnmmsnDtKYMpchIiONaW89Ul35Db6+5Jnd7uznHuqPdDOL2drOLuk9cHDreRSgpGZ2W\nYe6B7HD0hzAOx+yYb923OUZ4A4uBm2P0/woOhcDvNzcdCf+9WlpbMK5dQXk6b7TG4+LQqakEM+eg\n09Ig2T2m3pQEuzXq3d23ouxx2O64A+uiRQQbGug99RH+H/4Qy/z52MvKMFImf0nSnmDPgDDtML/3\nB+uA0A0/1tzVhD/gpyvgxz9Md/CQPyPqRkvV4sA+ICiT7ck4LHbswzzmGBjEEwzW6WqsU7B2KKW+\nrLV+bYzvWzjBukSMC3V2Erp61QzlixcJtZodNUZiIkZmJrYlSzAyM25rpaZUVxwtHVMQ0l1+aGsz\nN8hoa8PwNg0fxH0tvux55jaR8S5zfrWMIB87wzC7ul0udEYmAP3tXq3N8G5tA58Xy7XLWE9+YAa3\nxYJ2pxDMzIa0NHRq2rAffOLjLLR3T8+W9M2UxYK1sBBLQQHBhnP0fvABnd//Ptb8fOLKyjASEwc9\nv69b+EYLdUCwDmm9mt+bu/vCNXytNXx7tC5hq2HFaXHgsDpxWp04bU4cFgeZ8ZnmfWtfV7BzwNeA\n51sd2C12GSA3grGO0lgPHFFKFWmtvzraE8N7T8vfthgk1N5O8Oo1M5gvXiAUHnVtJCdjZGViu+su\nM5QnMDgmLz2ejy63Rapks6u1vQ3V1opqDYdxW9uN68KGgU5IIJSSJkEcDUqBK8HckGTOHAKLl5gt\nb58P1dyI5eoVrHVn4NQJAHRyMqHMbHRqGjotHZxOnDaD623dtzjR1AqGgviDfddYB3wP9g1gCt8v\n6sRafxHne/9F1399i2tZDhpy42m39vQPagrq4IjnMZQxODQtThxWO267G6fL2d8N7LQ6cVgGBqtj\n0DVamzELen2iaEwhrbWuU0qtwtwZax2wHdh/83VqpVQ+5u5YI86pFrNDqLWV4LVrhK5eI3DxQv8g\nLyMlBcvcudgyM7BkZqLs9oidMzPRzttnm2//hd3d0NqKag2vx+1tNsO47zqxxYpOTCDkTkXPy4XE\nRHNN7nhXRAcyiQgwDHOKXUoKgaIF5rH2NlRTE6q5Ccvli3DmY/N4fDxF1y2cIgkdCKDGObK8r8Xq\n7fb1h2v/4KXg4C5f/4DbLd3NdAW7za++667BbnpDo3e/2424/q5fR5Kd+LIk5ly2M/9cO4UtPXQV\nzqV7YS5x8QnhgL0Rrje6lJ1YDau0XmeAMf9Uaq2rlVKFmCFcDuxRSvXtdtWM2cW9CnOTjbcmo1gx\nfYXa2gheuULw3DmC16+j/X5QCiMlBWtODkZmBpb0DJR98lqYmYl2vP7e0ddi7u42VxhracHSeM1c\nj9sfXqnMYq7HHUpJRefkh8M4ybxmKr/MZq6ERHRCIjov37zW3d1lhnZTE+6G4+RfPsXF79bTnZZA\nZ24WbekuOumio7eTjnDXcMfAgUx9wRo0BzI1nm9GE8JVN/KlGbvFjsMSvn7a991qJ9XuDh8LH7cO\neNzSd83VPmCgkx2LGuZaewno3l44XYc+XQfHQuh758OiIgniGe62Pjpqrb3A6nCreiOwEnNrSi9w\nBFijta6KeJVi2tF+vxnKV64QOH/ebCkrhZGairUgH0tGJkZ6GmoKu30zEuPQesBazIEA+LyovkBu\naTa7sMG8Zux2E8gvlPW4ZyCtNd0hs/XZHTK7gruCnXQFzdvdQT/+0OBjnSEf3eFw7Qr66Qp2Ebw/\niMMfxNPcQ9qpHpLeCxCwKBrT4/BmJxJIS+rvDnZaHMRb40kNT71xWOzYLQ68hhe7EUdBUf6AEL4R\nyFM1kEnZbHDnIigugFOn4chxVP15Qo8sRyXILlwz1bj6d7TWlZibbYhZon9K1JWrBC+cJ+T1AuY1\nZeu8eRjZWVgyMqY0lAfVpzWpwS7meS/Te8SPtccX3kJSm93V7mQCOfnmetx9i5uIKRfSIXpC3f3h\n2heWXSEzWLuCfrpCXeFwNcPUDNeuG63XkNmCHWmkMIDdGNAatZjXWxOtiaTb0/vvOywOzlyCI1eD\n/Onv5hNvdeD0a+IvNeGov4zR1g2heHRhHuTNHzHoPL0NABSk5U3K39ntUvY4WHoXzJ+Lfq8G9aND\n6LKlqOL8aJcmxmHEkFZKFWCux10/deWI6UJ3dxO8dm1oKCck3Bh9nZU5oYFeE6qvq4tgU5O5YcGV\nK4Samklu97PkykWCV3MJLchDL1iMTk2FhES5djwBWmt6Qt3hlqv/Rgs2HKjd/SHbdeN+0E9XqJWu\nUHd/wHYHu+gOjT5Iy6Zs5hzXAd3CdouDBFuCGayGPfyYI/zYzWHsIM6wj7nl2tPczi9aWliYXGh2\nCycAGTnoe+5FNzajztajTp2BD05BWiq6KB9y5pit1mlOpafC6k/C8ZPwzlFUwwVCD5WhnLIu+Exy\nq5b0FqXUMsxW824J7Nilu7sHt5RbWgBQrgQsWZnYFi/CyMyMyub1Ohgk1OIN79N7iVBzM7oj3G3t\ncGBJSyXu3ntwJ6Xws2unKC4tYm7x7S6UNw2FQtDWivJ6za0zw9+Vz9e/NKjq7jGXCO0xv1RvD87G\nK4QCvQR7ewkFetGBXnRvLzoUImgxCNkMghaDgE0RsCoCFoXfaafXCt3WEL0WTY+h6bZqug1Nh0PT\nmmChLclKa6KV1iQbbYlWQhbz0oCBgT3c9Wu2YO1meFqdJBsp/d3Cdosdh3EjWPu/h4+NeL11Etmt\nCg30BDR2241LHUopyEiDjDR0IAiXrpiB/V41HDEgdx56yUJUcuLIbz4NKJsVyu6FednoI8dRP3wT\n/cAyVO68aJcmxmjEkNZae4BNAEqpEiSwY06wuYXgxQsE6uoINZujopXLhSUzE9uihRgZGRgJ0ekW\nDnm9BC9dItDQYNYW0mAxMNwpWPMLzIVN0tIGfWiwAda4Wtqm0TKPI2r1YTQ0oK5dgcZGVGMjqvE6\nzoazBH0+Qm3t6I5OlB7apdvrjCMQZyFoVfTaDHqt0GNT9FoV3VbosUKvQ9GbqAhYDXpt5veQYcES\n1Em2I+cAACAASURBVNgCGnsgiL0X4gKKuADEdfSQ0KtJCYA1EMIa0FgCISyBENauXiz+Yeafu+Ix\nkpOwpqbStnA+oexM9JwsQlkZhLKz0Blp0/4av91q1tfVq7GP0DhWVgvkzjODudMPDRdQH59F1Z+H\n7Cx0ShIqcXpfPlFzsmDNo3D0ffjFO6j8HEIPLkNNx9X5xCBjnYJVw43AXokE9oyku7vNwV4XLxE4\n14D2+1E2G5Y52dgWFJst5SiFsg4EzO71S5cINjSYvwwtFixZmcSV3YeRno6RMvxqYwMl2KO3OIXZ\nLdxDV8gcrNTV007o0nlUQz22cxewn79Eoucczqs+7O03un17bQZtyXH4kq20JFnwZltpTUqkPcFN\nu8tKR4KV9gQLnS4bxDtxWuzEE4dLxZnfiSM+fNuJDSc2XMrWf9uBDeeA+3asWG5zIJPu7SXU1k6o\ntc38ajO/B70+gk3NuN6rIdjcbP67han4eGx5ObQuX0pw8UKCSxaiszKmVXDHWc2/h55ACLh1YKl4\nJyxZgF5UBOcuot4/hXr/Q8iZg56TNa27kpU9Dh4qg/rz6Or3MVrbCT320LSuWYxj4Fh49HYVSGBP\nd1prQs3NZou0vp5QUxNojeF2Y83LwzInGyMz85bBN1lCHR0EL10m0FBP6Np1CAZRrnis+QVY5s7F\nyMy47U/6zjgL/t7b2+2pL1xvXHMdMKAp6Ker/76//zmdIV/4umwXuqOD5AteUi63knXFT9a1brKu\ndDOvsQdbwGwJd8UZXM9y0JDlxLc4m7bsJLqykgmmJhEX7yJe2XFhJx4b85V9UPi6MAPYQfTmtSqb\nDUtqCpbU0ZegDPn9BJtaCDU10XvhEoH6c8T/x08J/d+9gLminDUvh9b7SwguCQd3WvQuTdjCP169\nwdtba1kZBuTnoHPmon/+Nup0rTlA646FsLDI7GaeplR+DiQlon91GPWfVegVD6Fu8e8qomdCP0kj\nBHYB5txpCewo0F1dBC9fJlBfT/DyZXR3N8oWh5GdRdx9ZViys6NyXRnCHxoaG/u7sbWvFRQY6enE\nlZaYwZyUNKb3MldlCi9bGP7uD/qxJtZyRYf41fXTQwc6hbrwB310h7rDA5m6xzRS2Kqs4eunTtxd\nBjnnu1h8rp05DT6yahtJuNaGCr+8NyWBYFYaenE6OjsTsrOJz55LmjudPCP2uxYNpxNjvhPmz8V+\n7939x4NeH4GGc/TWn6O34TzOvd9Ht4cXuHEnY8vPw3d/CcE7FhG4905z/fIpYA1fV7/dkO6jLBZU\nQS56fja6qxd14mM4U4e+5w7+f/beO0yys7zTvt9z6lQOXdW5Z3qme3JUmBmhNBLKAQwGhIRF2jUL\ntvnM2sbGBryOu+sliF1sf/68Bnm9xmAECCEsTBAKCAkUZ0ajGU0OPdPTOVXOVef9/jjVYXKH6q7q\n6ve+rrrq1KlTp57urq7fed4n0dFesQvgyyFCdXDHTfDiLnjyOeQN1yBWLq+0WYoLULbLPeVhVwZp\nmlZC1Xj8dnQUsDp72VavQm9tRWtoqNiXhczmKPT1kew6SqbvDIVMmrwhyDT6SW8Jkaj3kBEFMsXj\nZIYPkOm3BDVdzJAppInmxj3W3ETziEt2ZaqHFPB4jzaRjOSYSFayhDZgr7tAopO1Pb7fqTtxx7K4\nj3VjP9pF4NXXKXQfpThi/X6x2zHal2PbfDXGve3YlrWiNzejOcvXQa2W0OsC6HVbJ4TbWuUJl0S7\nm8Kpblz/8h2rCY5hYF+3hujdt1K44RrM5W3zZpcxR5EeRxh2xPq1yDUdiNffRLz6Ohw6jty2xYoH\nVyHC5UTecj1iz3544RVkJApXbFLNT6qMeVmTuYhgB6WUM5mkpbgEMp+ncPw4uf37kckkwrCjtzZb\nseWWlrKXRo0vCSfySZL5JIlCgng+QTKfJJ6Pk8gnSOSTJPIJRjOjpDMJ/F1DBE+NUijkiHt1xoIG\n4ZCduLc0p3gQ60ap1WGp8cN4drAlsHb8hrc0+cY+pTTHMaX8ZvLxP/4sRiyp8al7W2f1ZSMGhzCe\n/QX+HzxJ/uQpa5/LhbliOY6rr8S2YjnGinb05saq9ZIWA0II9PoQen0I5/arAOszVhwYJPvmIXL7\nD+D6m6/A//p7jNWdRD9wH/nbboYyXwTZtPKI9DjC44adb0GORRB79iN+/hK0NlsZ1fPYbW+2CF2H\na64Cvw/2H0JEYpg3Xas+21XEvAdOpgq2Yu7IdJr8kSPkDxxAFgrYVqzAdsP1aKHQZf+xxsfNTbQ6\nzCdJFkqiW9oeF+Gx7FhpWs7kIPXCRZr1OzQ7bsONx+bGrblYNlxgeVcCl+kju20z5qp2nN66Uh9h\nBw7dPinCugO7ZpStI5PbZjJazM5MoItFbL98leA/f4vcwcNg09E2b8T/6x/EWN2J3tigvIsFQAiB\nrbUFW2sLnjtvxcxkyL15kPQLL+H+r19CfOnvybzjLnLvejtm54qyvKde+tgVZ5bGcFnGl5PlmT7E\nq68jfvwM8rodiJbG8r5RmRDrV4PXg3x5N9ovX7O6lKnPfFVQvdkNM6QUC38vcBLoBB6WUkbnemw1\nUJRFUqNDpA7sI3vkEDkzT2pFE/GORlL2JKnEy6Qi4/NcLREOZ8dKcdrMRLP/i4msQJTaHrpwG27c\nugu34aLO3oLb5sJtc+OyufCU7t02F55xUba5sWnWx0gODCN274N4HLYsQ169xcqGXUAMXVCYplck\nIlGMJ36C59vfxxwdw+xYif8/fgDHVVegLbDdivPRnE6cO7bh3LGNwtAw6RdehB89jePb38dYu5ro\nB95L/tadMAcPdbzHe9Esjyd9LqK9DdkQQry0G/HcL61Y9ca1VSmAYlkLvOUq5Mt7wOexupYpKs68\niHRpXOUDwIkFHLbxqJRyx5T3fxS4qwzHzojxiTjZYnZy+k0hNZHcNHVOq5X0lD5rbmskFy4dkyZT\nyGILx2nsjtIwmiNnaPS1ORlocVDIaXDUek+7Zp8Q2fExcl7DQ6Oz3uo7rDvPGknn0l1WC8SSGDt1\n55w8WRmNIXbvRwwNQ0MIef0tiFBdReaVWiJ96WP0Q0dp+OrXyby2BwD7W7bjvuUmjI7yeGeK8mNr\nasR336/ifefbyO7dT+r5X+L+iy8gHvo7Mu+8m9y73465YuaJT+NzWOZJo4FS7PfWG+DQMatca3jU\nStSqwgxw0b4MMjnY+ybS4UBsXFtpk5Y88xWTjgIPCyE+JoT4ipRyXv/SpWYrE/9mUsqoEGKHEKLj\n3IS1mRw7le8efZQX+DnZYpZoNkyumCNbzJEzc+SKWXJFS5hzxdwlM4XB8lwdtsmuSw7dUXrswGd4\naXI1EogWaOkexRMJoAU2Uti+Am1lO06HG1fpdZecirMAyFzOSpLpOgNeN/Km6xDLWio6TNzQIX+h\nb1zTxHjyZ/i/9m0Kp06Tqw/hfee9uG68Hq3KG1EoJhGGgfOabTiv2UZhcIj08y/CE0/i+NbjOK/Z\nxvCnfhtzBt20tHn2pMcRQsCmdcg6P+LFXYif/hz51uuqcvCFWNtpdbTbvQ/pdCDKFFpQzI45i7QQ\n4qNSyn+80HNSyodL86fnmx1Y4zKnMj4+89Qcjp3gWPgoqWgCu27HrtvxO/zYNTv2UrKTo7Tfrjuw\na9bjs4TYNhmDNTTjkstdUkrEd/8dhBdx2xZY3VGVy2OcOA1d3citG60lvCpINrnYcre+/xDuv/wi\n2oZ11P32x7Bv3VwV9ipmj625Cd/978L7rreTfvEVkj/6Kf4P/RbJP/0U+TveOq1zjHvS5Y5JXwzR\n1oK86xbE8y8hnn4BeddbFzwkNB3Elg2QycJLu5EBn6qjriDl8KTvBC4o0iX2lOE9LkfdBfZFLrJ/\nJsdO8NErPsLWq7Ze6pCyIYRA3rYT8bMXkYeOI+oC0FCFvag72uHAUUQ8UT0DLC52MZOz2lr6P/wg\nekP9AhqkmG+EYeB+605c111D7OvfQv7J/yB76CiZj38EbNVXmy78XuTtNyF++nPrf/zOmyo2Pe6S\nbNuKiESt8qy33VGVy/NLgXL81u8XQrwXS4yfBp46Jw59vAzvcTkiwLkqVlfaP5djJ+jt6cfrW9hm\n+nLjGrRXX0c++gRi41rELGJu841sDCJefR1p06vianukL0JyJMHpU2fXUet9A7gLBaKDQ2i5C/Sh\nVtQE8t47ydT5SX792xgvvEz4v3wSeYkhGImMSWZsiP4zObqEe1bv2XOmd3a2drYjXnwNvv1vVolW\nFfbRls0NVmeyx3+I3H5lpc1ZNPT29JftXOUQ6U8D3wXuwPKqf7OUjLUH+Dawmkt72uVgF/Ab5+wL\nYWVvz+XYCf72ob/De0786I57b+fOe++YmaUzQDgdmDvfgrb/EPLAEYjGYNP66vpnXt5mDRx44yDy\n5uuqy7YpTAyqqMKogaJ8CCFw3XAdttYW4t/6Hv5P/imJP/5dims6K23aeQivB/mWqxEv7ULs2oe8\n5sqqC8EIj8eKpe8/hGxtQbRVZ2OWSvLUj5/m6R+fXWWcKHXTKwdzFmkp5UOlzYdLN4QQ24DbsTKm\nbwc+Ptf3uYwNrwshJparS9snxhPBSsliESll1+WOvRhf+Ju/WrDl7vNYuwp56oz1z3y6Bzaug4ZQ\n1TRHkA0hxE+es8ZHbt5Q0WWxhkEfnkiMlR1nd6myDY7gsdmor6vD1tZaIesUC0ZbK8VNG4j+wz+R\n/4svkvovv0/+AhfUsXQRZ8hGa3sLnavmlsTVuWrlrF4n21oQz70Eo2G4dtucbJgXVq1EGnZk3yBy\n21arYYtigt/4xEf5jU989Kx9+/fu556dv1KW88/LZZuUco+U8iEp5V3AZ+bjPS7A/UKITwkh7sPy\n7qcmrH0WuG+ax16QodQwsWyMonmZ+p55QnS0I+++FYom8sXXkE88iXzyOeSe/cjuXmQme9lzzJtt\nAT/yyk2Iw8cR3/8xvLgLOTxaEVtMU14wLF1csRzhdDD2uf9J/HtPYMYTC2+cYkHRg0GCn/odnDu2\n4f6r/4XoGzjvmIkFlgqusIiWJuRbrraSMAeGK2fIpdhxBeg62u59lbZkybEQLs9XF+A9kFLuBfaW\nHj52znMPTPfYi/Hwvn/k0eSjCARuw43P7sNn+PDZfXjtXutx6d5r+PDbfQSdQZy28g0KEKE65Dvv\nsjzW4VHk8CiibxBOnLIOCAVhWQu0Ni/4MHqxYQ1yxTJr6ft4F6K7Bzxu5OoOa9DAAmWw5osSQzv/\nG1e2NhP9/texf/Mx+PbjpH/2PK5bb8Zz561oC5xroFg4hGHg/8AD5A4fpfG/f5mhv//CWc+bJZXW\nK7zKLFatgJOnELv3Ie+9tfqWvQ0DNq9D7t6HjMQQddMbhKOYO7XUzGRe+d1tn2TZhjZi2RiRbIRY\nLko0G2UoNchQaoh4Pk48F8eUZ9dyuG0ugs4QQWeQkDNYug8RdFjbQWdwomPXdBEeN3jcVnY1WF50\n/yDiVA/y0FHYfwjpcVvxo7YWa2l8Af7px2ftsnGtdQFxrAtx4AjsOwiNDch1ndDWMq9x63xRTkw2\nOhdZFyD7/3yE3Pvvw/7IY/Ct71tifctNeO66TYl1jSIcDrzvfgex//sN9L37KU4JW42XR2tVUOIo\nt1+JePJn0NUNqzsqbc75dLTDkRNoe/Yhb9tZaWuWDDXRzGQh6Ax0sq310vEiKSXJfJJoLkokE2Yk\nPcJIepjh9DA98W5ORroIZ/eQzE8mFQgEvpLXHXKGCDlD1LtChJz1NLjqqXfW4zIu7YUKpwM6VyA7\nVyCLRRgagd4BON0Dx7rAMKCl0RLslsYFKfcQjfXQWI/MF6CnD3H8FOKXr4HdjlzbCZ0r5qWRw6VE\nehxZFyD78Y+Qe3BSrFM/ewH3rUqsaxXntTtIPfs8dZ/7G0Yf+epEyeB4fXSlPWkAEQzAynYrCXPF\nMst7rSKEpsHm9chX9iBHw4j6yldzLAVqpZlJVSCEwGv34rV7Wea9eNejTCFTEu+SiKeGGUkP05vo\n4c2R/YxmxiiYhYnjXTYX9c56Qs4gXrvP6ptteHAbnlIP7SnbhgdHS5M1Hm/HlcixCPQNIE73IF/Z\nYwXfGuutuuv6INQH5/XLQBg26FxhXUTEEtaS3pGTcOAINDUi29uguRHhL0/Xr3xRTowfvBwTYv3+\n92L/5mOIbz1O6tnncV67HcfmTdg3rkNzqySZWkBoGr73vYfwF/8a46nnyN99G2BdWEN1eNIA8qrN\niJ4+OHQMrthUaXPOp73Nyj3Zsw955/QaxijmRq00M1lUOG1OlvvaWe5rv+DzpjSJZiMMpoYYTg0x\nlBpkMDVEb6Kb/mQ/qXySZN6aYnWhFqSa0PAYnsmbzY17pQf/cjsNY3nqRobxncnjLAqcNieO+iac\nLdY8ZBpC8ybawu+Fq7Ygt2yc9K737AdpgsuFbGu2PP6mxllnrmfzEodtZl+4MuAn+/Fft5bBv/N9\n9H9/iswvXgZNw+hciX3zBuybNmB0rKy6WKFi+hgdK8EwECOTSY3jIyqne2E33wiX05pJfbqnKkVa\nCIFcvxr56uvIZEplei8AtdLMpKbQhFaKY4fYENpw0eNMaZIupInnYsRzCRKluHgiFyeeT0zsH0kP\nMZQaoiufJCWSJIMpzDoTZ7pIIJqn7kwB/5t5fAUbDpsTGfAiG0OIpgZsjU143XX4DK91s3vxGT68\nhmfWQzmETbfiWx3tyEIRRkZhYNjKvh1PgqsLIJe1QFODdeEwzTh2NF3A55ydXTLgJ/uxD5P92IcR\nfQPYXt1D8MnnSD31HMknfoxwu7FvXIdj80bsmzagV0HzFsX0KZzpgXye4pWT050KpeXuahFpwFrl\nOnICmc1VTZnlWbQ1W+GCnn5Yv7rS1tQ8tdLMZEky1WNumUF4V0pJupAmlotNJMBFs1ESYwNk+3vJ\n9/eS7++neKKXTDHLMWeBsYBG1G8Q89soGBqa0AgYfgIOPwG7dauzBya2A6Vtv913STEXNh1amqzb\nVZuR6QwMDiN6BxDHT1nL4roOTQ3I9mXQ2oRwXTxjPpOXNPrm7u3Kthby73obQ+96GxSK6AcPY3tl\nN7Znf0Hs698CKdFbm3Fs2mh52mvXVOcXqmKC3LHjYLdTXL9mYl+hyjxpwKrSABgLQ2v1NQ8RhmGF\nqLq6kUqk552aaGaimBlCWGVkbsNNi6dl8omVwNWTD814nOLgEMWBAdJnTpKKjZIeypDw2gjX2Rh2\nCwZceQayg3TFThPNxYjnE2ctwQsEfrtviogHCDoC1rYjQNAeoM5Rh8fmRghhCXDJywZrDCb9Q4gz\nfYhXXwek5WW3t1nCHqo7a/hIKmfispd5SdqmU7xiM8UrNltedjSGvmsvxsu7kb94mdQzz4FNx75m\nNfZNG7BvWo9t+TK1NF5l5I6ewL6qw0qkLDG+3H25ZMOFRHjc4HDAWKQqRRqA5a3I1/Yi05lLXjQr\n5s58ZXfvwfKkHxJC/OF8vIdi/tF8PjSfD2PNapzcSCCRoDgwSPFMN8WBAeRYBmGzoTVsRmtsRG9s\nQAbriJoJxrJhwtkIY5kxwqXt/lQfpxLdvDEaPU/MbcJGwO7HW1oZ8No8pW03Xr8Xz5VefGYQfziH\nbyCK6+BhbPsOWd5rYz2ypQkaQ6SzRVzG/H7hyoCfwu03U7j9ZtJSop3qxvbqHuxPPU/i338C33sC\n4fVgdHZgdKyYuKms8cqRP9VN/ugx0h984Oz91ehJAwTrEIMjsHl9pS25MONL3r0DsKaj0tbUNDXT\nzEQx/2heL9oaL8aa1UgpMcNhin39FHt7KBw+Qn7/ftA03MEgvsZGOpta0Zu3IhyO885VMAuEs2HG\nsmFGM2OMZS0xj+cSxPIxRjOj9CR7SRSsBLminNLpTQfRKKlLwLK4TuNhQWg3OISdG0cETurZl21F\na2rGXt+I1+7HY/PisXlmXJN+WYTA7FxJrnMlw+97N+Ry6PsPYdu1F+PVPaSe/bnVfAbQ6kMYKydF\n27ayHc1VfWMKawlpmqSefo7E4z/AtnwZuXe/7azncwVLpJ3zfGE3YwoFsFdXCdZUhN2O9LghFq+0\nKTXPvIt0qWZaUWMIIdBDIfRQCLZsRkqJjEQoDg9T7Omh0N2NPHwYAC0QQGtsQG9sQmtuQnO5sGk2\nGl2NNLoaL/teUkrSxTTxXJxYPk48lyCejxMrLa/H83GG01Higz0MHB9C5k+R3/0msligoAuiAYNo\nnUEkYMP0efAaPjw2Dx7Di9fmnRBwb+mxlRjnw627Zz7H226nuP1KituvZPg3/wNIidbbj37oKPrB\nI+i73iD5w58gszkQAr25CWNluyXaHSswli9Xse0yUYzFiP3ffyV34BDZD9xP9Lf+w1lL3WBVAwA4\njCoLTcTiyPWrq3oejHC7ENH4BepLFOWkHHXSt9VKVzHF7BFCIIJBtGAQY906AMxkEnNwiOLwkCXc\nx08Almjrzc1ozU3oTU2Xba4ihMBtc+O2uWnm4jG606Mpfu2NXbzjbWvZ0OQmPzZIZuAUhb5T5McG\nyY2lSekFIgEbI0Ebg740ffYoiUKCRD5BQRbOOp8u9LNE22f48Nl8Zz82fPhsfuz6RX4GITCXt2Eu\nbyN/5y1kAIpFtFNn0A8dQT90DM/ufWR2vw6FIgiBra0VW0m4jY4V2Ja1VV1ji2one/AwsX/6OkhI\n/vVfUbhux4WPK5joWnUtd8t0BpHLQaDKW2+6XdZkPsW8Ug5P+gvANWU4j6LG0DwetFWd2FZ1wrXX\nItNpigMDFAcHKXSfQR49CkKghULozU2WcDc0IGyz+1iGU9YMaZ/TZtU4N7RiNLTCluut5cNwGDE8\nhD7YhxgIQ78El4tiSxtmfQPZkJ+4vUA8HyVeiBLPx4gVosTzUaKFYfpT/RzNHyVRiJ+9/A7YNQd+\nw4fP8JfE24/f8E/ZZ+136k40Xcdc3YG5uoP8r9xtCXc+j3byNPqho9gOHcW9Zz+Zl18D0wRdx7as\nFWNlydte2Y6tpXlBOsctNmShSOKJH5L6ydPYN21g5KE/R9afOz5+kmxh5nX18060tIRc7f2xXU5k\n/1Clrah5yiHSYSHEp4DvXm7co2JpI1wubJ2d2Do7cVx3XSl7fJDi6dMUTpwkf/AQaBp6fT1aU6O1\nPN5QP20vciyZAyDgusDH2maDxkZkYyOFTZuhkEeMjMDICPpgH3rXCQzA63LTEgpRbGxCBldBfd1E\nC8lxpJSkiskJMY/lLSGPFSLE8lEi+WH6U33E8nGyZubs3wHWqoDH5sFt81jL7jaPtVJQ58FzUxOe\nWztx296Dp2AjcHoY19HTlnC//ibpX7w0MbpJqwugNzZYt4Z69MYGbKXHwuuZ+VL9IkWaJoXuM2Re\n2UXmtdcxEwnSn/go0fffd97f7lys5jdVttQ9OGx9Xqu9UYjLCZkMUsol81mrBOUowboLQAhxuxBi\nm5Tye3M3S7EUmMweXzMZ0x4cpFBaGs8fODgjT3solsVh0/DYp9H4xGYgW1qhpZXClq2QzSLGRmF0\nFH14ANv+N6BYBF1HBgIUG1sgFEKGQginqxTH9tLCxdu/AuTMLLF8lFg+QrwQI1mIkywkSBYSJAox\nEsUwfalekoUkqUKKrHn+yFGtU8O91o3nPVsIFLfT3pulZShHaDhD3VCC4JnTOPbtQ0ukJ14jnE70\nxvqSgI8LeQg9WIcWrFu0CWtmJkOht49CTx+Fnl4KZ3op9PYjs1k0v5/0224n9457MKc5oMIq2ase\ngZHhKOLQMeQVG6tf+HJ5cNir385FTtkSx6SUzwghAiWv+qtSShWsUEybs2LaG6wua2Y0SnFoiGJ3\n99medkMDWnMzenMTWn39RD3yUDxH0G3M7kvD4UC2tkFrGwW2WgIdjSLCY+jDg9i6T8GRg9axLjcy\nGKTY1IwM1kPd+d72OHbNQYOjiQZH07TMKJh5ksXEhJBbtzipYoJEafvMmjEOdyRJ5E2SBQcF2QQ0\n4UwXaRjJ0TiSo21U0jxWpHnwDMHdx/CE02jmZIqPdNghGMBWV4cRDKLX1aEFA2h1dejBgCXkXm/F\nar2llBRHRifFuKeXQk8fxeER6wBNw9bSTGrrRor33k5x0zprutUMJ6xl8rL8dfWzRJom4uXdEPDB\nhjWXf0GlSaURbpdKHJtnyprdXcrk/pIQ4j4hRFgllCnmghYIoAUCGGutIWpmOGx52mfOTJR8WXXa\nDegtzcT7o4TcZUqw0vUJz7mwuvSFmUkjxsZgbAx9ZPDC3nZ9PTIYBOfsPFWbZhDQggSM6bUclVKS\nNTMkCnHLMy/ESeSt+xOFOHsLMeLFMVLpGLaRMezDEQKRHHWRfOk2RGigj+DhAr5oDn2KkJu6RjHg\nRfg82FxuDKcb3eVCOJ1oTifC5UA4nRM3zVXadjmtyWymROZykM8jc3lkLocc387nkLk8TNmWuTwy\nnabQN0Cht9cawQoIjwejfRmp22+muHYVxTWdmJ0roAwx+fR8NL+ZLYePQzSOvPuWxdEIJ5W2lrwV\n80o5srs7zo1FSykfE0J0CCE+B3xOedWKcqBN8bSllJhjY5gDgxS6T5N/8wCBl09yo8eBVj+KbGhE\nNjSAP2BN/ioHTheybRm0LbO8bdOESOTC3rbbTbGxGRmqR4ZClh3z8MUrhMCpu3Dqrml566Y0SRWT\nJPKxiXj6kUKMWD5KIhclP9aLNjyCMRLBORqnLpzFk8ziyqRwZEw8cYknI3FlTByZAka6cJaHPmMM\nA2E3EIaBsNsRDjupDWsx776F4hpLkGVjffn+hueQzpkEPQvRLuLSyGgcsf8wcvM6a2TlYiCVxmxr\nruoysVqgLNndQoh/wOrRvR1YBewA6gCB1dNbZX8ryooQAr2+Hr2+HmPzJqRp8svj/8Zb67EypQ/s\nB1OCYUPWN1hi2dAAdcHyiaWmne9tp9OIsVHE2Bja8ADizOmJDG3p92M2NiPrgsi6IHi98yY+Dd7M\nCAAAIABJREFUFzVZaHhtPrw234Xj6Z2Tm6Y0SRUSxEsiHilE6c5HiecjxApRovkR4vkY6XQMkUzi\nTBdxZoo4sybOjImu2dCcbmwuDzaXB8Ptxe72Y3f5cbj9uNwB3PZSzbrNjdvmQRczW66eK8mcSXto\nYd/zXOTwKOIXr1qJYpvWVdSW6SLzeWQiAT7Vu3u+KcsUrNLtJFYr0NexuoydlFK+XobzKxSXJVeU\nnMaNbcsq8huaoFBARMKIkWH0wX60o4fhYBFsOjIYRIYaLNEOBsFWxhpklwu5bDly2XJMrrBKv2JR\nRDhsedx9PXDUavKCYSADdRQbGmFcuD0zmJQyz2hCw2v48Rp+Wl3LL3lszsyRKMRKCXHxKQly1n2i\nmCBeGCNZ6CeVOU4qkcLEPO88Tt052VzG5i21iLWS9LyGtyToXryl7Pi5dpBLZky8s5yaVg7kkROI\n19+EhhBy51umPe2t4pzpty6Cl7dW2pKapxwi/UUp5WfKcB6FYtYMRDNIoNFfakFqs5WWvBsxN2yy\nvNloBG1ctE+dhCOHQYCsq0PWN1rx5FC9NdygXNhsELLOK8GSpXweImHE2Bj66BC206fgcGmZ3O5A\n1gUoNjRDXZ0l3IsgE9uu2QnZGwjZG6Z1vClNMsW0JeLFyQS5ROmWLMSJFUbpT/VZCXP5JHmZP+88\nTt1VEnGrY5zP8OG1efGW7n2GD4/NGrNq1xznJRUmsyY+58ILo8wXEC/vRvT2Izetgy0bFkccepxT\n3YjmRquhiWJeuaxICyFuA3ZdLK6sBFpRDfSGrf7YTb6LJLJoGgRDmMEQ5tr1Vq1xPIYYGcE22IfW\newaOHwNA+v3Q0GAtkzc0gKPMyTGGAY1NyMYmCpTmhWcziHAYIhH0kSFsJ45CtlSO5XQi6+ooNjRB\noM5KTCu3TQuMJjTcJW/48o1hLbLF7GRyXKmcbSJhLh8nWhzhdOIUiXyCZCF5nqduCMNq/VoScI/N\ni1YPAzSxZ6SJgGGNVvXb/Tgu1kGuDMhYAvH8y1aN8c5rEYvMG5WxBHI0jLz5WhWPXgCm40kLICKE\nOAE8DTylaqEV1UbfWBpdEwQ90/xyFQL8AaQ/QH5VKa6WTCJGhhEjw9gG+xEnTwIg/T6onyLas8zc\nviQO52TdNhutfek0IhKGcNjyuI8ehpzVsAWnE+kPUKxvKP0cfvB4Z1yCtJhw6A4ceiP1jsvL+rkJ\ncvF8jHjB2k7kY0QLI/Qke3HWR3k9uYs9B88WdKfuwG/3Twh3wO7HZ/fhtXnw2S3P3FtaijelecmZ\n6eNIKeF0D2LXG9bf765bEH7vrH8fFeP0GetCc9niurhYrFxWpEv1zw9hTRreA4gLZXQrFJXk1EiC\nloATTZvDtb3Hg/R4kCs7yAGkU4iRkZJo9yG6uiaPa5iSjOaZpy9alwvpcpVqtzdb+5JJS7ijUfSx\nEWynTkK61MRECPB4kX4/xVC9Jd4+n2XfEms4cdkEOaAnnOYzLx3mf39wK6uabUSyEcI5a6xqJBux\n7nMRBlMDDKQHieXiJPLJs0asAiT7U7h0F62jzXgNT0nAvRPb/oKN4HCaQHcYd6qAc0UHthuuRRiV\nzyqfKTKeQB7vQm5Ys3ji54uc6Sx3vwd4RUr56Usc84eAlFJ+qZzGKRTT5UBPlLZgmT1clxvZvgLZ\nvoIc2yGTsTzt0RErrn36tHWc02mJdsnbxuebP1EsXUiwbDkT40DyeSs5LR5HxKJo4TFsR49ArrRc\nrutInw+zLoT0+cHnQ3p9VpLaYoqDlplYxvoN1nuc+O0u/HYfK2i/5GtMaZLIJ0sT2OJEczEOew+T\nyCdxN7uJ5WIMZ4YZCHdj7xvF1T+GK5LG1AThoEHPchdx8SrGK4/jN3yWkNt9+Azv2bcJb92Kq8/n\n8vt0kaYJr+4FpxO2bKi0OUuG6VzK3Sml/PilDpBSPgQghPg88D9UXbRioemLpLllw/S6es0apxO5\nvB25vB2TqyGfR4yOTGaQ971hZbw67BNetqxvgEDd/HqyhjG5HA+TkdhMGkrCLWJxtPAIoveMJeoA\nmob0eMHrpRisL4m3F7w+K+GtxomlLZEOeaaf3a8JrRS39k3sa0jUA7B6eTuFnjMUBk5gDg1ZLW2v\naobOlSQbfcRIE8vFiOZixHIxIrkosVyMwfQgQ+lhTsS6iOcSpIvp897XoTssr9xuTV7zG77StrXP\nb/jwle7dNtf8tOo8cgIZjljNVpbA56NaKHfHsc+U2oIqj1qxYMTSeaKpfPk96cthWP2/ZUsr5har\n3EqEx6xJW0P9aAcPQNG0arVD9ZMx7XLWal8Kp8tqwNLYdAHxTiDiMUQibnneJ49BZsowEJd70vv2\nei3v3e0Gl7tm4t6xTAFDF7in0+v9IshcjsKZHgrdp0npujU/vKkJ+/XXY1u+HFGqFHAB08t7h7yZ\nP0/Mo9kokVyESC7GUGqA7sQZYqX56udOZNOFPinedj8Bw7ofj62PC7zf7selO6cl6DIcRR48ityy\nAdFw8aliivIzHZGe6V/kMSHEe1RymWKhODmUAFh4kT4Xm80SxMYmzE1boFi06qNHrUlb2tEjcPDA\nZK12faPVjSwYsrzhhWJCvBvPFu983vK8E3GIx9AjY+i93ZBKTUzesl7vRHo8mF4/0m2Jt3S7F52I\nR9N5/E7btL1OKSUyk0EmEpijoxR7eykODpLrH0AL1mG/7nr09uVzHl5iaAb1znrqnfXTsilZSBHJ\nRojmokRyUSLZyISwD6QH6EsNcDhyjGgudt7MdEMzpoi2b8oSeykxzvDgkw6CLx7C6XFh27J+Tj+b\nYuZMK7tbCOGf7hK2lLJLCHHfHO1SKKbNkf4Yhi5oq6uymk1dt5a8Gxow1284u1Z7oA/t5HE4bC09\nS7+/VE8dQgZDVlx7oTGMiQ5qwGTM2zSt5LRUCpFKWvfJBFoihhgasDzwqSLucFgJby4XpsuLdDis\n2nOHY3Lbbge7o6Ix8UgqT53r7Isjmc1iJpLIZMIS41QKcyyMTCWRqZS1MgKg62ihIPZt23Dm82hu\nN8bqhe++JYSYENPll5nIJqUkVUhNJMRFslHC2fDE9kB6gJ5kH4l8kkQ+QaaYxciZbDkQx5Etsm+r\nn9xzP7ZE3O7FMy7o9lKinN1rxddLN6/hxW24p5X5rrg40xHprwAPA++bwXkvfwmoUJSJo/0xlgXd\n2PQq/zK4UK12Io4YHUWMWSMytVOlDHK7YS2RB4NWg5Vyd0abqd0ej5W0VqpqPssDnyri6dTEvZaI\no42NILJZK4mtWDz/3CWxlg4HpttjCblus7zxi9ykroOmg23KfqFZdphF632KpW3TRBRL+0xz8t40\ncZ44xWazQPrpCDKVQiaSZ9toGAiPG83jQV+xEuHxILweNI8H4fdPNB/RTpyY119/uRBC4DE8eKYh\n6ADZWIToM0+S7ogT3XklO12SWNaKo8dyMaLZKCPpIfqT/cRzcRL5BKY8u5TN6lo3RcANL95xITem\nirq1HO+0TW/5fSkx3RKsLwgh/pOU8v9M87yr5miXQjFt9p+J0Nm4COtNhQCfH+nzIzs6J7qRWb2/\nR9GHBtCOH4P8Iaszmu8cb7sCvb8vyFkibnGWiI9TyFsNWrI5S7izGchNbmvpJMSjiMK40I7fCuee\naW7oOmganp5R2pv9QAi9bRma12MJsceD5vEiHJXPqK4UZjxO8bkX8OoeGt97H9o0VnZMaZLMJ4lm\no0RLy+/WtvU4Voqn942Lei5xXjmbTbNZCXClunT/FBGf+thv9+O1e5eElz7dxLEHgF1CiNVSyj++\n1IFCiACoRjSKhSFXMOkNp7ll4zxndi8UhoFsbkE2t2Bu3DzZGW1sFDE6im14ADHubRs2ZLAk2MGg\ndV/OlqblxmZYNw9nfTWPb5/fyXsK54r2uLdcLCKKBev3pGlIbdy71ixvW9Osbd1mbWvaxIXN04n9\nvO+aZbh2rpynH3hxYkYiZH72HMJux3nvvWhu97RepwltwjNe7rt0r3colbPl4kSyUSLZsJUcl40Q\nLdWpD6T66U8OcDR8lHguQbaYPev1AjEh2H7HuDdeujlK3eMc1mOP4Vm0gj4tkZZSnhRC3AE8LYS4\nH/g88Oi5cWohRAfwKHDRmmqFopwcHYhRNCUdi9GTng5TOqPJjlVWk5VczsoiL43ItPqQlzqReTzW\nEvm4eAcCiyaR65KML2tfgNkMyiwUTeKZAo2+Kr6oqQCF093kXtuF8Hlw3X0Pwjl/7Wc1oeF3BPA7\nAqxgxWWPzxQyVlJcScjD2TCRTIRINkwkazWdORk5SSwXO0/QNaFZfd3tpYQ4+/j22TF1b2kJ3mPM\nfXhLuZi2FVLKPUKIVVgi/DDwVSHE+OSrMawl7juwBm48Ox/GKhTnsr87gqELVtZP72q/JrDbJ73t\nDZusfYmEJdxjo5Zw9/VaNduaQAbqrISwCo7IrDbGUnkk0OxXIg1QHB4mt/cNzJER9PZ2nDffjLBX\n13K/0+akxdZCi6flssdagm6JdzgTniLuUeK5GEPpQQaSAyRyCRL5xHllbAAum2synl4S74l4ut3q\nAT+xbfdiaPOTMzKjSwUpZQS4s+RV/yZwO9aYygiwC7hLSvlM2a1UKC7Ci8eG6Wz0Vn/S2Hzj9Vr1\nzO0rrGXjYhERjUyOyBwcRIwnONls1uSvumBp0lbdkmsdOpa0suqblrgnbUYi5N7YR7GvDy0UxHX3\nPehti78ntyXorbR4Lv+zjGe9jyfDjSfHRXNRYtkYsVyUodQAvYm+iQS5cz318ff0lYQ70ZUs288y\nK39eSvk01rANhaKinBhMcN0aVUxwHrpuZYeHrN9NEaxl8kgYURrYofX1Tkz+wpgq3MFJ4a5RxpJW\neKBpiXrSZjJJfv+bFLq60LxenLfeir5y5ZLMrJ6a9d46DVEHyBYyREuiPp4kF8tGiWSjxHJR9sm9\nZbPvkiI9k/roy1HOcykUAIPRNGPJHKubK1BTvBix25FNzcimZsyJEZlZS7gjYWuZvLcHjpWE225Y\nS+SBwOSkLa+vJmLco8k8brs+p25jixGZzZI/cJD8sWMIux3HjTdiW7Nmcc2yrgIcNidNNidN7gsn\nrO5hD//Ev5TlvS7nSb9PCHFirjHm0kzqVcA/zuU8CsVU9pwKA7C+RYn0rHE4JuPb60sjMjOZiwu3\nJqz+3uMTtsbF2+1ZVMvlI4kcDd7qirnOF1JKzKEhCidOUujpQWga9m3bMDZsQCxkpzvFrLikSEsp\nHxZCfEwIcSfwlZmOpyxle38GOK4mZCnKza6To7SH3Phc6oumrDidkz3JxxPT8nlENGrVMcei6GOj\naEODkCsN67DpVh13SbSlPwA+7/zM3i4DQ/EsjTUu0mYqReFkF8WuLsxEAs3nw75tO8aa1fOata0o\nL9NpZvKwEKITK5tbAk9hxaNPXqAEyw/sALYBd2FVR/yWlLKr7JYrljRSSl4+PsL2TtXsf0EwDGs4\nSMM5k7bSKUQsVpq0FUUbG0Xr6Z5sn2nYrLGY4+Mxq2RM5nA8x1vXT3fkxeJAmibm2Bjm4BDFwUGK\nQ0MIXUdvb8fx1reiN9VIL4ElxnTrpLuAu0qNSh4AvgisKpVkTS1THM/yfhS4X0oZLbO9CgUA/ZE0\no4kcG9v8lTZlaeNyI11uZLNVFlMEq+1mMmHNt47HEYkYWngMrb8P8qXuYZqwarq9PmvZ3Dsp4vM9\nbMQ0JSOJHMvqFrc3KaVERiKWIA8OYQ4PI/N5hGGgNTVZ8eaVK6uulEoxM2ZaghXFqpF+eHxfSbhR\ngqxYSHadHEMIWN+qRLrq0LTJdqelXRNVqJk0Ih63pm2VJm1pPWcgNWWGsmGbmK41MWHLNXk/165q\nY6k8BVOyrNoGslwGaZrIaJTiyCjmkCXMMpsFXUdvbMC46ir0lha0UEglgtUQc26pciFxFkJ0zDR+\nrVDMhKfe7GdVoxePozq6AimmidOFdLrg3BnXhYI1bCQeR6RSkEqix6OIoUFEOg2FKc0mbHppNKYH\nXC5L0F2uaU/YGohZNa7toeoTaVksIlNpZCZt3adTmPE45lgYMxq12qEKgVZfj7F5M3pzM1pjI6IG\nMu4VF2a+vuE+DXx8ns6tWOIUTcnB3ih3bL585yHFIsFmK9VnBy/cyzuTseLfyWRp0lYSPRaF0VHL\nE89fYAiH3bCy1+0l4XY4kQ4HyZ4kbbFRmjJRimNpS+B0HWEr9fe22crqiUoprRBAoXi2+Jburcdp\naxJX9pz+1IaB8HjQm1sw1q9HC4XQgkGVlb2EmJVIlxLEvoCVJAZQd84hq1AirZgnDvVFSWaLbGk/\n92OnqFmcTivrPBi6sIgX8pDJInLZ0qSt7MSITJHNoqUSEAmj5bJohwa4IZyk8PMiF52vpWkI3TYx\nnEOUxmKKiR7iEkyJLBbJ9PVC0SR98JAlyFPGYcrS/YUQTifC7UK43OhtbQiXC+H2oLlcCI/beqzi\nyUue2XrS4/XO37nAcwL42CzPq1BclleOj+Cy66xuqt2OWIoZYjPAayCZ/ExMzWid2pn5Oe0wulnE\n/d6tyHzeEtVCATlxb1qTtgpFZLFg7SsUrfnU4/sQlmBrAj2fB6Fh6+goTd/SEJr13PhYTKFZk7ks\nIS4JsFqiVkyD2Yr0a1LKhy72pBBC1cUo5o2n3xxg87IAurZ4mmcoqoe+SJo7t7ROa0bydLAfP27d\nr1lTlvMpFFOZbeAlcqknpZSfmeV5FYpLMhTL0DWcZFtHsNKmKBYhuYLJSCxLZ6On0qYoFNNitiJ9\nUghx1cWeFEL871meV6G4JC8cHkITcNVKJdKKmdMfSSOhduePK2qO2U7BekYI8SkhxMNYzUvO5QFU\n4phiHvjxG31saPOr0ivFrOiLWPXYHQ3Kk1YsDmblSQsh/gH4YyAM1J9zW8352d4KxZyJp/Mc6oux\nrUOlPChmR184TcBtqH7vikXDrN0RKeVFvylLIq5QlJXnDg1impIdql+3Ypb0jKVYHnRX2gyFYtrM\nNib91GWe//Qsz6tQXJTHd51hQ5ufoEfVjipmR89Yiq0r1EKfYvEw67Y6pYYmF0PVSSvKymg8y6G+\nGNetqa3JRYqFI5svMhzLqvp6xaJitsvdJ4DPCmvI+2sXeP43ATU/WlE2nj4wgCYE16xSS92K2dEb\ntjK7VzeXpz5aoVgIZivSe0r3ESxBnkodZzf7USjmzL/tPsPW9jqV1a2YNT1jKQSoGmnFomK2y917\npJSalDJ0gZsGPFZOIxVLm1PDCU4OJdm5Ti11K2ZP92iKljonLru60FMsHmYr0pdLDPvcLM+rUJzH\nj/b24XHoqoGJYk6cHkmyol550YrFxaxEWkr5zGWef3125igUZ1M0JT94vZdrVzdg6GqQvWJ2SCk5\nM5ZiuyrfUywyLvutN5sWn6otqKJc7OoaJZzMcaNa6lbMgZF4lnSuyNoWlTSmWFxMJzizQwhxJdYI\nyumy4/KHKBSX5xu/6KK1zqnKZhRz4tRIEoB1LZeqHFUoqo/piPR2JrO5p4OgAtndQohO4L3ASaAT\neFhKGb3Ma+6TUqoktyplOJZhz6kwD16/klK5n0IxK04NJ6lzG9T7HJU2RaGYEdMR6ZPAF2ZwTgH8\n4ezMmROPSil3AAghAsCjwF0XOlAIcR8QAr4ihKiTUsYWzkzFdHliTw+6JrhxrVrqVsyNUyNJNflK\nsSiZjkg/LaV8eCYnFUJsm6U9s0IIcTVTvHcpZVQIsUMI0SGlPHXu8ePes+oxXr0UiibfffUMN6xt\nwK1qoxVzQErJqZEkv3bdykqbolDMmOmky87Ei57La+bCDmDsnH1jwKrLvE6toVYpLxwZJpzMcdum\n5kqboljkjCZyJDIFNrSpeLRi8XFZkZZSds30pLN5zRy5UMf8yEX2KxYB//TzE6xr8bFSzf1VzJET\ng3EANi4LVNgShWLmVPU6ohDiY1jzqc9NRBtPTntKSvksliCfWwBZV9pfFj75yU8SCJz9T/7ggw/y\n4IMPlustFCUO9EQ4NhDnd+5aV2lTFDXAscEETX4H9V6VNKYoP4888giPPPLIWfui0UvmLM+Iqhbp\nGcTCdwG/cc6+EFbS2yXfYrq2fPnLX2bbtgUNtS9ZvvrscRr9Dq5WHcYUZeD4YJxtHaqJiWJ+uJCz\ntmfPHrZv316W89dEC6dSh7OJpW0hRB1wYjxpTAhxdalE61xUTLrKGIymee3kKHdtaUHT1J9HMTey\n+SLdoym2tqvIl2JxUtWe9Ay5XwjxKaALK5Hs/inPfRZ4ldL4TCHE7cA2LE/6s0KI8WVzRYX59svd\nOA2dm9Y3VdoURQ3QNZykaEol0opFS82ItJRyL7C39PCxc5574JzHzwDPAA8tjHWK6RBJ5njs1W7u\nvqIVl12vtDmKGuDIQAyXXVczpBWLlppY7lbUBt96+TQAd21pqbAlilrhaH+ctc0+dBU6USxSlEgr\nqoJ4Os+3Xz7NbZua8bmMSpujqAGKpuT4YJwb1jVW2hSFYtYokVZUBY++2k2haHLPla2VNkVRI5wZ\nTZHJm2oOuWJRo0RaUXHi6Tzf+EUXt2xsps5tr7Q5ihrhSH8MQxdsVJ3GFIsYJdKKivPNF09RMCW/\nclVbpU1R1BCH+2OsbvbhMFQSomLxokRaUVHCyRyPvHSKOzY3U+dRXrSiPJhScqQ/xk4Vj1YscpRI\nKyrKv7xwEhC8XXnRijJyZjRFMltkW6fqNKZY3CiRVlSMoWiG777azT1XtOB1qoxuRfk43GfFozer\noRqKRY4SaUXF+PwPDuCy27jnCpXRrSgvB/tUPFpRGyiRVlSEAz0RXjw2wnt2LMdlr5nGd4oqoGha\n8eibN6jWsorFjxJpxYIjpeTzPzjI8pCbt6ovUkWZ6RpOkM4VuWZVfaVNUSjmjBJpxYLzzIFBjg3E\nef/1K9WkK0XZOdQbw2noqj5aURMokVYsKJlckb/+8SGuWhlk83KV1KMoPwd6o2xo9WHT1debYvGj\nPsWKBeWfXzhJJJXn/devrLQpihoknStwdCDOHVtVMqKiNlAirVgwukeSfOMXXbz9qjaaA85Km6Oo\nQQ70xCiakhvWNlTaFIWiLCiRViwIUkr+6+P7qfPYVeMSxbzxxpkwbXUu2oLuSpuiUJQFJdKKBeFn\nhwZ5syfKB2/oULWrinlBSsm+7gi3bGqutCkKRdlQIq2Yd+LpPF/8wUGuWlHH1R1qbKBifjg9kiKS\nyqulbkVNoURaMe/89++/SSZv8uGbOittiqKGef10GLddV/OjFTWFEmnFvLLr5Cg/PzzE+65bQb3X\nUWlzFDXM3tNhtrbXqdIrRU2hPs2KeSOVLfCX39vPhlYft2xUncUU80c4mePUSJJ7VVKiosZQIq2Y\nN/7+6WNE03k+8tbVaEJ1FlPMH3tPhxECrl+j4tGK2kKJtGJeePn4CN99tZv739KuaqIV887rp8Os\na/ERcNsrbYpCUVaUSCvKTiSZ488f28eW5QHu2NJSaXMUNU4mX+Rgb5R7rlBL3YraQ4m0oqxIKfnj\n7+ylaEo+eota5lbMP/vPRMgXpRpNqahJlEgrysoP9vSy51SYX795FUGPWnpUzD97ToVZHnKzLKS6\njClqDyXSirJxajjB//zRIW5e38iOzlClzVEsAQpFkze6w9ypwiqKGkWJtKIsZHJFPvXNPYS8dj5w\nY0elzVEsEQ71xUhmi6oVqKJmUSKtKAt/9t03GIpl+cSd63Cq3tyKBeKVE6M0+52sa/FV2hSFYl5Q\nIq2YMz/a28vzR4b58M4Olqu4oGKBKBRNdneN8bar2hAqQVFRoyiRVsyJrqEEn//BQXaua+Sm9Sq7\nVrFwvNkTJZUrqjI/RU2jRFoxa5KZAr//r7tp9Dn48M6OSpujWGK8fHyUtqCL1c1qqVtRuyiRVsyK\noin5g3/dTTSV53fuXqdmRCsWlGy+yJ5TY7zj6mWVNkWhmFeUSCtmxVefPcYb3RF+6/Y1tARclTZH\nscTYezpMtmBy59bWSpuiUMwrSqQVM+ap/f187YUuHrhuBVeuULN7FQvPS8dHWd3kVYmKippHibRi\nRhzui/Lfvv8mN6xt4N4rlBejWHgSmQL7zkR45/bllTZFoZh3lEgrps1AJM3vfn03y0Mufv3mVars\nRVERXjs5iikld2xWWd2K2keJtGJaJDJ5PvG1XRi6xu/dvR67TX10FJXhxWMjbFkWoN7nqLQpCsW8\no75pFZclVzD5z1/bxVgiy+/fu17N7FVUjOF4hqMDcd7zlhWVNkWhWBCUSCsuiWlK/uiRPRwfjPN7\n96xnWVAl6igqx0vHRnHYNG7ZqBrnKJYGSqQVF0VKyf/70yO8cnyU37xtDetb/ZU2SbGEkVLyi6PD\nbO8M4bLbKm2OQrEgKJFWXJSvvdDFIy+d5gM3dnDNqvpKm6NY4hzpjzMYzfCBGzsrbYpCsWAokVZc\nkO++2s0/PHOMd+9Yrmb1KqqCnx8eotnvZFuHqs1XLB2USCvO44ndPXzph4e4e2sLv7pNtV1UVJ5k\ntsBrJ0e57y3tqvRPsaRQIq04i5/s6+NzTxzgtk3NPHj9SvWFqKgKXj4+QtGUvO0qddGoWFookVZM\n8NP9/fzl9/azc30jH9rZoQRaURVIKXnu0BBXrQzSoGqjFUsMJdIKwBLoP39sHzeubeAjN69CUwKt\nqBK6hpN0j6b4oEoYUyxBVB2Dgid29/C5Jw5w47oG/tNbV6NpSqAV1cNzh4YIee1cu6ah0qYoFAuO\nEuklzrdfOs2Xf3KY2zY186GdHcqDVlQV6VyBl4+P8OGbOtHVxaNiCaJEeokipeSfnz/JV549ztuu\nbOWBa1eoGLSi6nj5+Ci5osk7rlYTrxRLEyXSSxDTlPztT4/wrZdO855rlvPOq5cpgVZUHVJKnj04\nyJUrgjQFnJU2R6GoCEqklxi5gskffnMPr54Y5cM7O7hdjftTVCnHBhN0j6b45L0bKm13nCOXAAAb\n8klEQVSKQlExlEgvIaKpHP/5a7voGk7w23euVa0+FVXN028O0BJwcu1qlTCmWLookV4i9I6l+MTX\nXiOeKfCZX9nEmhZfpU1SKC7KWCLLaydH+b17NqhqA8WSRon0EuDNngi//409uAydP33XZloCrkqb\npFBckmcPDmG36bxddRhTLHGUSNc4P9rby+eeOEBHo5ffuWsdfpdRaZMUikuSK5j87NAg79y2DI9T\nfUUpljbqP6BGKRRN/u6nR/nWy6e5eUMjH97ZiaGrBnOK6ufFY8MkMwUeuG5lpU1RKCqOEukaJJrK\n8clv7OZwX4wP3tjBHZubVYmVYlEgpeTJ/QNc3RFkechdaXMUioqjRLrGODYQ5/f/dTfpXJE/evtG\nNi4LVNokhWLavNkTpS+c5k/etaXSpigUVYES6Rri31/v5Qs/OEBrnYtPv30jjX7VAEKxuHhyXz8r\n6t1cvTJYaVMUiqpAiXQNkMkX+bNH3+D5I8PcvKGRD93Yid2m4s+KxUXPWIr9PVH+4r6tKjyjUJSo\nGZEWQnQC7wVOAp3Aw1LK6EWOvR3YVnp4DfBpKWXXghhaZs6MJvmDf91DfyTNR29ZxU3rmyptkkIx\nK378Rh8hj507VBc8hWKCmhFp4FEp5Q4AIUQAeBS469yDSs9tk1I+VHp8H/AUsGYBbZ0zUkp+sq+f\nL/zgAAGXnT979xZW1HsqbZZCMStG41leOj7KJ+5ch01VISgUE9SESAshrgbk+GMpZVQIsUMI0SGl\nPHXO4TuAzwMPlR4/Day6yLFVSSKT508ffYOXjo9y47oGPnRjBy57TfwpFUuUn+zrx2lo/Op2Ne1K\noZhKrXyz7wDGztk3BqwCTk3dKaV8Rgixfcqua6zdi0Og954O8yff2UsyV+S3bl/D9WtUX2PF4iaW\nzvPc4SE+dGMnbketfCUpFOWhVv4j6i6wL3KR/Ugp9055+EfAb1zuDT75yU8SCJxdzvTggw/y4IMP\nzsDM2ZPJFfnqz47zyEunWNvs47O3raHB51iQ91Yo5pMn9/UjgAeuW1FpUxSKGfPII4/wyCOPnLUv\nGr1gOtSsqGqRFkJ8DFjNlKXs8adK+56SUj6LJcihc46pK+2/3Pm/I6X8P5ez5ctf/jLbtm273GHz\nwu6uUf7ie/uJJHPc/5YV3HtFqxo6oKgJktkCzxwc5L1vWUHAba+0OQrFjLmQs7Znzx62b99+kVfM\njKoWaSnlw9M8dBfne8MhrEzvC1LK8B6VUn5vlubNO4lMnv/2+Jv8/PAQ61p8fOreDbTUqeEYitrh\n6QMDFIom77+ho9KmKBRVSVWL9HSRUr4uhJhY2i5tnxiPM5cSyyLjZVZCiG2l132v9PhjwLellLGF\ntv1CmKbkx/v6+NufHCFbMPnwzg5u3dSMpmpHFTVEJl/kp/sH+NXty6lXoRuF4oLUhEiXuF8I8Smg\nCyuR7P4pz30WeBX4UqmeehcgSw0TBBCegdc+rxzsjfJX33+TE0MJrl1dz69dt4KQV32BKWqPZw8O\nks4V+eCNnZU2RaGoWmpGpEvJYOMJYY+d89wDU7a7gKorxByKZvjSDw/y/JFh2kNuPvuOTWxo81fa\nLIViXsjmi/z4jX52rmtUIRyF4hLUjEgvVkbjWb72wkke33UGu03nwzs7uGVjM7pKDFPUMM8eHCSZ\nLfC796yvtCkKRVWjRLpCRJI5vvHLLr7zSjc2TfCOq5dx19YW1ZREUfNk8kV+uLePm9Y30hZU4ygV\nikuhFGGBiafzPPLSKb754ikA7t7ayr1XtuJRTRwUS4Sn3xwglSvye/dsqLQpCkXVo5RhgUhmC3zn\n5dN8/RddFEzJ7ZuaedtVbfhdRqVNUygWjFS2wI/e6OPdO9pVLFqhmAZKpOeZXMHk8dfO8I/PHSed\nK3LrpmZ+5ao26jyqcYNi6fHk/n5yBZP/ePOqSpuiUCwKlEjPE0VT8tP9/fx/Tx1lNJHlpvWNvGub\nqgdVLF0SmTw/2TfA/deuVC1tFYppokR6Hnjl+Ahf+uEhzoyl2NYR5A/u3UBbUC3tKZY2P3qjHykl\nH9r5/7d3/8FR13cex5+fKlnauiGSQC2hDklAk6oYfs+1FVsJdYRW5bdXO553V9TOdeauc6h4P/rH\n3RxFbed+dKYVsb3ajmcBhXO0mVaCp0BtgRBJrSYISUBIy48QQiJiFsjn/vh+d/my2U12N7v57m5e\nj5nMZL/7+X738/3uj/f381vjokUSpSCdRu+fOsual95h3+HTXHdNkH++6wYmXxP0O1sivus6G2Lr\nH47xtc9N4mo19YgkTEE6Dc5+dIGfbG/hF789zNWfGMW35k9hZtlYjKbxFAHgpYZ2rvyY4V7N0S2S\nFAXpITh9NsSWPUd4/reH6L3Qx10zSrlj6gQKrsy6Cc1EfHOs6xyvNx3nmzXXEdRoBpGkKEinoO3E\nBzy17QBvHjiJMYbPX1fCndNKNce2SAybdh+h6JMFLJ+j9aJFkqUgnaDafe281HYl7xzt4mjnOYo+\nMYq7pk/kS58dz1WjVToQieXgsR7q2zr5zqIbCYy6wu/siOQcBekEbd5zhM/eNJYp1wT5SnUps8rH\ncuUVqtYWicday/O/O8y1xZ/g9qkT/M6OSE5SkE7QP955AzfefKPf2RDJGfVtnRw8/gE/+IuZWjBG\nJEUqCibI6EqJJOz8xT427Hqfm68tYlZ5sd/ZEclZCj0iknbb3jlOR08vq796g99ZEclpCtIiklbd\n587z0t6jfLFyPGXjr/I7OyI5TUFaRNJq4673MQYeUSlaZMgUpEUkbd471sOO/Sf51vzrNf2nSBoo\nSItIWlzss/xsRxtl4z7JnTMm+p0dkbygIC0iaVH3h2Mc7fyQ7yy6SUOuRNJEQVpEhqzzg1421x/l\nS5/9FFWlY/zOjkjeUJAWkSF77s3DBEZ9jMfuVGcxkXRSkBaRIWk41El9WyerFlZplSuRNFOQFpGU\nnQtd4Oc7D3HTZ8ZQc8M1fmdHJO8oSItIyjbtOsLZ3gv8y9KbMUadxUTSTQtsiEhK3jvWw7Z3j/Pt\nOyr5dNHH/c6OJCgUClFXV0d7ezsA999/PydOnKC2tjbyeNSo9DVbPPvss4RCIebMmcPUqVPTdtyR\nQiVpEUna+Yt9/PcbrVSMv4qls6/1OzuShKampkiAHjduXMZfzxijWpYhUElaRJL2y31/5Hj3Rzz7\n0J9pTHSO6e3tBaCkpIS7774746933333Zfw18pmCtIgk5Y+nz/FyQztf/3wZkz8V9Ds7koQtW7bQ\n0dEBQEdHB+vXr+eee+6JmTYUCrFr1y5aW1sJhUIUFBRQWlrKnDlzCAaD7Ny5k6amJgoLC1mxYkVk\nv9raWtrb2yktLWXBggVxq7vb29tpbGzk5MmTABQWFjJ79mxKS0szeAVyj6q7RSRh1lp+uqOV4mCA\nv7q13O/sSJKqqqooKSkBIBAIUF1dzejRo2Om3bx5M83NzYRCIQKBAKFQiLa2tkjbdVVVFQDd3d30\n9PRE9gtXpVdUVMTNR1NTUySYh0IhQqEQHR0d1NbW0tzcnJZzzRcqSYtIwnbsP8n+P/XwX/fNJDDq\nCr+zI0mqrKyku7ubjo4OgsEgs2bNipnu1KlT9PT0YIxh0aJFFBcX097eTm1tbSQoFxcXU1BQQCgU\nor29ncrKykiABigvj30TFwqF2LlzZyTNvHnzACIl8x07dlBRUZFQ57WdO3fS3d0NwPz582Pu09jY\nSGFhIWVlZQCRm4ze3l5qamoIBvvXBkXv4yeVpEUkId3nzvOL373P56aUMLui2O/sSAYVFxdz//33\nR9qTm5qa2L17d7904dJ0a2srcKkUPXHixLhBNly9DTB37tzI/7Nnz46ZJp729nZ6e3tZsGABVVVV\n7Nq1q1+aLVu2XJbv5uZmJk6cyIIFC6iurqaxsXHQffymkrSIJOQXvzsMWP7p7hv9zooMg7feeuuy\nIFZYWNgvTUVFBY2NjZHgfPToUeBS8I7FWzXuDeQFBQWRknl3dzcTJkyIPFdbW8stt9xCMBikubmZ\n7u5uAoFApHd6SUkJ+/bt6/daixYtuizglpaWXla9X1DQfznV6H38ppK0iAzq3fYz/Oa9Dv729krG\nXhXwOzuSYW1tbZEAPX/+fFauXMmCBQv6pSsuLo5UFzc3N3Pq1CkAJk2aFPfY3url8+fPR/4Pt01D\n/xuC8vJympqaAKfUXlVVRW9vb+RYwWDwsuA/0GuPGjWKuro66urqBryZyBYK0iIyoAsX+/jZzjau\nuybIV6ap5+1IEG7nDQQCkYAbDpLRwm3P4ermeG3RYd6x2du3b4/8762ujh6/HW7vDgfxYDBIIBCI\nBOaenp6Ybcvx1NTUsHjxYurq6hLexy+q7haRAb369jGOn/mIJ/58Oh/TmOgRIdwDvLe3l/Xr10eq\nocO8/1dVVdHY2BjZNliQLigo4JZbbmHHjh20trZG2rPD5s6dG7M9OxgMsnv37kjpt6SkJNITvKOj\nI6GJWXbv3k1hYSGVlZUxq7qzkUrSIhLX6bMhXmo4ytLZ1zLlGo2JziexZgELbwuPhw4EAgQCAcaP\nH8/ixYsj6VpaWiL/B4PBSFA3xiTUI7qyspKFCxdSUlISeY1x48axcOFCrr/++pj7VFVV0draGjl+\naWkpgUCA2tpampqamDNnDj09PWzYsCHu61ZXV9Pa2sqWLVuoq6tj7ty5g+7jN2Ot9TsPWc0YMx3Y\n++KvXueGm6r9zo7IsPph3QHebT/D5r+bq2Uo4zh48CAAkydP9jknki0aGhqYMWMGwAxrbcNQjqWS\ntIjE9PaRLna1nOLbd1QqQIv4REFaRPr5sPcCP9neStWEQu64ecLgO4hIRihIi0g/z715mA97L7Jm\nRbVWMBLxkYK0iFym4VAnO987yd8v0DrRIn5TkBaRiJ5z5/np9jZuvrZIY6JFsoCCtIgA0Ndneeb1\nFi70Wf5tuaq5RbKBgrSIYK3l5785xO+PdPGvy6ZSEtTUnyLZQEFaRKht/COvvXucR796A5+bMvjM\nTSIyPBSkRUa4Nw90sHHXEf761grumjHR7+yIiIeCtMgI9m77GZ55vYUvXDeOb3ypwu/siEgULbAh\nMkK9c/QM//nr/VRNKGTtPeooJpKNFKRFRqD6tk5+VHeAqtJCfnDfLK68QpVqItlIQVpkhNmx/wQ/\nfqOVWWVj+f69Mxh1pQK0SLZSkBYZQV59+0889+Zhbq0cz5oV1Vyh9aFFspqCtMgI8dLeo2yuP8rX\nPz+Jv5l/ndqgRXKAgrTICPCr3/+JzfVHefC2yfzlrerFLZIr1Bglkud27D/B8789zH1fKFOAFskx\nCtIieWxvWyc/fqOVL1aN55s1U/zOjogkSUFaEvbK/77gdxayQq5ch6b2M/yw7gCzysZmbMGM559/\nPu3HzEUvv/yy31nIGvpMpJeCtCTslzkSnDItF65D5we9/Mev91M5oZDv3zsjY7249YPseOWVV/zO\nQtbQZyK98qbjmDGmDFgKtAJlwHpr7Zk4aecBRcBYYD7wXWvtW8OVV5FMe6P5JNbCv39d46BFclne\nBGlgk7V2JoAxZgywCfhynLRbgWnW2ka3CnATMHlYcimSYX19lu3NJ7h96qcJfnyU39kRkSHIi1ts\nY8w0wIYfuyXomcaYSXF2KbfWNnoen85c7kSG1++PdNF5NsSimZ/xOysiMkT5UpKeCXRGbesEyoFD\n0Ymttd5tFcCyAY49GqD1wHtDymA+6Onu5p239/mdDd9l+3XY9OZhCs+F+PBYCQ3HMzthyZkzZ2ho\naMjoa2S7I0eO0NPTM+KvQ5g+E9DU1BT+d/RQj2WstYOnynLGmIeBGmvt7Z5t9cAaa+3mOPuUAQ8C\n04DV8dqkjTFfA55Lf65FRCTP3Wut/Z+hHCCrS9LGmJU4Jd3oOwnjbttqrX0N6MLpBOZV5G6PyVrb\nBqx2A/w2Y8wka213jKS/Bu7FKZF/lMp5iIjIiDIamIQTP4YkX0rS04CnrbWzPNs6gelRVduRErS1\ndrXncQtOSfy14cu1iIjIwPKi45hbVV0UfmyMKQJawgHaGDPNDcbgtFPP8+xegdNxrH54cisiIpKY\nvChJAxhjqoEaoA2nI9k6T5DeCOy21n7PffyN8G7uPmuienuLiIj4Lm+CtIiISL7Ji+puEREZHsaY\nMmPMw8aYJcaYVe7kUYPts2Q48paPVJKOIZkpRqP2ewp4JE4vcZG8kMr3wxizxFr74nDkLxOSnHY4\npd+PXGGMqY+e3dFaG3N2Rzc4jwXWAUX59NuYwlTU092Hs4BH3RFGg7+OgnR/yXwIPfvMAzYCM6J7\nlOcyzYnuGK4vZC4YiT/SSZ5z0r8fuSKZkTRR+10Ers7V9z+WRN9n97kHrLVPuo+XAI9baxOaijqr\nx0n7IdYUo8aYme446kNx9glX90TPepYPNCe6I6Hr4D43PeoLuZU8uQ7Jfj/CpWe3liknJXPOqfx+\n5JikZnf0yOzUd8Msyfd5JrAWeNJ9XAeUJ/qZUJt0fwN9CONZZq3dxgj4IDIC50RP8jqEv5BhkS9k\n5nI4rFL5fkBufzeSOedUr0+uKIqxrSvO9nyW8PvsxoYZnk2znM2J3bQpSPeX1IfQU82dj5L6wUly\nTvRcMmxfyBwwEn+kkznnfL8+Sc/umKeSep+ttd7J/h8BHkj0hUZMdXcmphgNV3PnUztLlKR/cKLm\nRL+agavAcsWwfSH9kskpePNAMuec79ennv6f57E4fTUGkm+dn1J6n93v2UZr7Y8TfaERE6SttesT\nTJrMh7AGuNqdHMXglKyWGmPqon6os0qWzInuu0wHplS+kH7J0PfjspdIOlPZI5lzTvX65ARr7Vvu\njI5A7Nkdga4YHSVzubkjlqTfZ7fW9VS8RZ/istbqL+oPOOD5vwjY43k8DSiLs18fMMnv/KfxOkzz\nnru7rTPWOeL0eF4b9bgPuM3v8xjO6+B5fh6w2O+8Z+h6JP39cD8LhX7nfTjOeaC0+fAHVAOrgCXA\nd73fA5ymv1Wex/OAh4GLbtqc/z1I8TMxHZjnebwy0e+DhmDFkMwUo+62MTh3VWuBp3G61x8a5mxn\nhDHmgLV2ivt/EU7pcpb7OHLX7N4lrvU8VwNswPmgZmVJOhmJXgf38XSc4Sbb3McrgQ35cB0g6Sl4\nw8PR1gJPcKl2Iqckec5x00r+SPQz4VnEKRxsDXDaWluc0OsoSMtANCe6Y7i+kCIiXgrSIiIiWUpD\nsERERLKUgrSIiEiWUpAWERHJUgrSIjKs3KUOy9JwnGnukomL3b8l7vj88GusNcb0eTo0Ru+/zn1+\nlTGmcKj5EckEBWmRLGWMWZmnweNR0jO5yXpr7ZPW2s320gQR5eBMqmOtXY0zJPLR6B3dYZOngb3W\n2u8lMzzOfV8GXUNZJB1GzIxjIung/jivxxmOVYSzwlV4liGDMxXqWjvEJTrdsdUtAwUPdwzygzgB\nLzy3+FZr7Wb3OZulY5IHXM41kfPCGQp32bSs1toXY5TQ9wLzjDHV9vJZAGuAPTiTbSRrI85nYHkK\n+4okRUOwRFLgBtGncFb+OuzZXoYTGNZ4J7xJ8thlOOvPPhbn+XKcQNECrPQGcmPMYmAOzixPWbd+\nszvxy/JY55bEeY2x1vYYYw7i3CCtA+psjPW93ffJAl+21i6POp4BVlvP2shJnMdtODcbTw6aWGQI\nVN0tkppwCe+yOYndWcceBZ5wJ0BJxTr3rx93NrN64EfW2hXRQdit9p3GIKVwHz1IjHNL8rx63M3T\ncQL6WuC0MWZPjGpoa619BqgJNx24pfG6oZyEW0OxPE+bIySLKEiLpF+4+rsm2R3dIGNjVQe705HW\nAU/ZgRft2MQQg1AG9avqTvW8rLXd1tpvutO1FuEE+U1x9t3ApQURxqTpBmYdzk2HSMYoSIuk3wyc\nKta9Key7HKedO5ZNOKvo/MMgx+gc4Bi+GaAEm/R5uT27I+3J4YDN5Wt5ez0BPGSMmUT6VqTaBqxI\n07FEYlLHMZE0chcWWYvTeez/YjwfXnCiCyegbLXWvuhJsgxnDepYx52Hs3rOYPZm6YIOy3CuTcQQ\nz+sRnEAZPtZ0+t8EGHCaIYwxXcCjbjCPPJcq95hDHkomMhAFaZHUGeBBY0wLTq/uFcApoCZWr2pj\nzAM4y9WFS1/rjTEHjTE2aghRrJJeuLdzvOrciCwN0BC7V/dQzmtT1KIu5bjB3g2ejwPTjDHhNXzX\nhHd0b5YeAMqNMauAp1OsAm+L0XNcJG3Uu1skBcaYJTg9kSu8wcMYswGoj9Xr1xjTCSz1BnBjzFpg\nmrX2dvdxn7W2XzOUu++p8HKZucYNijXRvbrz4LxexWlL3zxoYpEUqE1aJL1WA4+7Q3Qi3GrdMTgl\nt5XGmG94hgc1uGnGEH+Sj6JwumR4Zt3a4LbHep971e0RHZ3Xhz37xJuta9A0UZYRu8d6SueVRbqA\nsX5nQvKXqrtF0shtpwQnKHmrvMNtlxtTrFbt4tKwr7iMMQ97S/HW2tXGmEdwSnuHopK/Gmcs9ws4\nbeoDdYpKJI1XvAlMUjovkZFCJWmRzCiPelzPpXbTmGJNxhG1f9x9vYeJsa0VqPBucKufX4hzjBoG\nL90mksb7WvGGhA3lvLJBEenrLS7Sj4K0SPp1ATO9G9xpQk8TY+x0eFGI8L5xJsh4HGdCjknxXtQY\nszZOybgBz9Akt1p9zAAdzOYz+DjrRNKExavqhhTPK4uq8ceiIC0ZpCAtkppinJJxUYznngaKPDNc\nhYcXLQce8wZht5TpHU9dT1SAB7DWbsMZ67spelYt46z4tApP7+UorVxeWl05SEenGgYfZ51ImrCy\neDcEqZ6Xu3gGxK/GnxWjh/0Lzq52hTsLWSyJpPGKe24i6aA2aZEkuSXf8NCh9caYdd4fdLcdeAxO\nB7IW3BKntXabG5SfMcbsxilxt0YFk00446j7DeGy1j5mjNnqvqbFCb6n3GMMNE/4Hi4NTVqCcxMR\n79ym45SyX4vaPiZcHZ9IGs+2eQxSLT6E8wpX40fGow9zNX4Zzk2VSMYoSIskye3ANGAnJs+EGdHb\n9zHw6kmb3L+YwckNjMmubNWKU7KfBJwepONavKC6HGflp0TThPWbwCSWFM8rXI3/DPhSjV9DAuO7\nRYZC1d0iWcQtidqB2mhTOOZbuBOvJLB05Xycea4j3Or6liTThA24LOUQ+V2N/yDOWHmRjFFJWiT7\nPIQz3vqhNB5zL/DdeE+6pdAHcErJLW6nqck4Ve9l1topiaSJOuaQV5sahN/V+E9l6UpjkkcUpEWy\njDvW+qAx5rYESr6JHnPANZPdIDRgNX4iaaIkVNU9BL5U47s3K0vjNWmIpJOqu0WykNthqjzOcKxc\ncXUmez77WI2/DGfNcJGM09zdIpKzjDF7cBYtiVmK9lTRr8WpDt9L/Gr8uGkyfiIicShIi4iIZClV\nd4uIiGQpBWkREZEspSAtIiKSpRSkRUREspSCtIiISJZSkBYREclSCtIiIiJZSkFaREQkSylIi4iI\nZCkFaRERkSylIC0iIpKlFKRFRESylIK0iIhIlvp/3AKJvfxjPJcAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fbbd5c0d208>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(5,5))\n",
"xmax=0.3\n",
"x0=-0.1\n",
"ymax=0.1\n",
"flavio.plots.band_plot(fit['DeltaM_d'].log_likelihood, x0-xmax, x0+xmax, -xmax, xmax, n_sigma=1, steps=30, col=1,\n",
" contourf_args={'alpha':0.5}, label=r'$\\Delta M_d$ ($1\\sigma$)')\n",
"flavio.plots.band_plot(fit['S_psiK'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=1, steps=30, col=2,\n",
" contourf_args={'alpha':0.5}, label=r'$S_{\\psi K_S}$ ($1\\sigma$)', interpolation_factor=1)\n",
"flavio.plots.band_plot(fit['global'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=1, steps=30, col=0,\n",
" label='global', interpolation_factor=1)\n",
"flavio.plots.band_plot(fit['global'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=(2,3), steps=30,\n",
" interpolation_factor=1,\n",
" contourf_args={'alpha':0.1}, contour_args={'alpha':0.4}, col=0)\n",
"flavio.plots.flavio_branding(x=0.77, y=0.09)\n",
"plt.legend()\n",
"plt.xlim([x0-xmax,x0+xmax])\n",
"plt.ylim([-xmax,xmax])\n",
"plt.axhline(0, c='k', lw=0.2)\n",
"plt.axvline(0, c='k', lw=0.2)\n",
"plt.xlabel(r'$\\text{Re} ( C_{VLL}/C_{VLL}^\\text{SM} ) $')\n",
"plt.ylabel(r'$\\text{Im} ( C_{VLL}/C_{VLL}^\\text{SM} ) $')\n",
"plt.tight_layout()\n",
"plt.savefig('CVLL-Bd.pdf')"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import scipy.optimize"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"c2min=scipy.optimize.minimize(lambda x: -2*fit['global'].log_likelihood(x), [0,0])"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def get_sigma_2d_chi2(chi2):\n",
" chi2_1dof = scipy.stats.chi2(1)\n",
" chi2_2dof = scipy.stats.chi2(2)\n",
" return np.sqrt(chi2_1dof.ppf(chi2_2dof.cdf(chi2)))"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"1.6300810234462912"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"get_sigma_2d_chi2(-2*fit['global'].log_likelihood(([0,0])) - c2min.fun)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.5.2"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment