Skip to content

Instantly share code, notes, and snippets.

@DavidMStraub
Created October 13, 2016 18:13
Show Gist options
  • Save DavidMStraub/dacbba6d098d595b118f26a8da59ca78 to your computer and use it in GitHub Desktop.
Save DavidMStraub/dacbba6d098d595b118f26a8da59ca78 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": [
"(-1.6933130732911199e-10+6.3442916109792175e-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='Bs',\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_bsbs': (Cr + 1j*Ci) * CSM }"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"observables = ['DeltaM_s', 'S_psiphi']"
]
},
{
"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 + ' future', [obs] )\n",
"fit['global'] = my_fit('CVLL Bs mixing', observables)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": false,
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 2min 50s, sys: 168 ms, total: 2min 50s\n",
"Wall time: 3min 31s\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+naQAAIABJREFUeJzs3Xl0W/eV4PnvDztAEARBUhQp7qSszdZC0kmcymZtdpJa\nEsuWo67UdPW0ZcnTPV2V03JJypw5s1ZZsjTtSU9125aq69RM15THluVTlVRSdpl0FidxYlOkvMiS\nLIrQvlAkSHADQCy/+eOBFEVxJyBwuZ9zeAg+PLz3o0jx4rfdq7TWCCGEEGLuMWW6AUIIIYQYmwRp\nIYQQYo6SIC2EEELMURKkhRBCiDlKgrQQQggxR0mQFkIIIeYoCdJCCCHEHCVBWgghhJijJEgLIYQQ\nc5QEaSGEEGKOsmS6AamilKoEHgfagErgqNY6OM65m4Da5JcPAnu11v7pXkcIIYRIpwUTpIFjWut6\nAKVUDnAM2Dr6pORztVrrQ8mvtwFvAzXTuY4QQgiRbgtiuFsptQEYrhSS7PnWK6Uqxji9Hjgw4usG\noEopVTHN6wghhBBptSCCNEbgDYw6FgCqRp+otW4E6kYcetA4rC9M5zpCCCFEui2U4W7vGMe6xzmO\n1vrkiC//DHh6JtdRSuUBjwAXgPAU2yqEEGJhcwAVwFta687ZXGihBOluwDfqmDd5fFxKqZ3Aa1rr\n/zLD6zwC/L/Ta6oQQohF4g+Bv5vNBRZKkG7idm94iA9jhfaYkiu8O7XWb8ziOhcA/vZv/5ZVq1ZN\np72LTiye4E/+nybuL/Hy9fVFKb32c//T99n/v/xFSq8ppi9TP4cb3WFebDzH9//gflYUee75/eeS\ny5cv8+d//ue89NJLmW7Konb69Gm++93vQjJGzMaCCNJa6xal1PCQdPLx+eQ889DCsu4R26xqk697\nI/n1TuDVya4zhjDAqlWrqK2tHecUMeTJTjc/OXmNFavXYjGnbjlEtsfDmgfWp+x6YmYy9XNYmdC8\nclZBbgW1tRX3/P5zicfjITs7W/4ezR2zngZdEEE66Qml1B7Aj7EA7IkRz+0H3gcOJ/dBNwFaKQWg\ngC6t9dEpXEfMwjfWFfP/vXeRT64EWV+em+nmiAXCbFKU52dx5npPppsiRMotmCCdXAw2tCDs+Kjn\nto947GeCVe0TXUfMzvKl2ZT4nPz6XIcEaZFSFflZfHSpK9PNECLlFsoWLDEPKKX4g7pSmi8EGIjE\nMt0csYBUFri53h2mPyy/V2JhkSAt7qlHHigiFte83zZ6O/rMffNbj6fsWmLmMvlzqCjIAuDMdcng\n+7u/+7uZboJIoQUz3C3mhyU5DjZU5PL6+5dYU+KhINsx62v+rgTpOSGTP4dirxOP08prv7lEXWVe\nWu5x6dIlOjo60nLtVLl8+TJVVVU0NzdnuikLXn5+PmVlZWm/jwRpcc899+R6vvvir/k/3/yM//Fb\na3BYzZlukpjnTCbF19cV8fr7l7nWFaI415nS61+6dIlVq1YxMDCQ0uuK+cvlcnH69Om0B2oJ0uKe\ny3HZ+MEf1fPfHnmPoz89z7/ZshyTsdJeiBnbtLqQn5y8xg/ePMPBHRtSeu2Ojg4GBgYkJ4IAbu+D\n7ujokCAtFqaqJW7+18fX8WevtPCj5qv8QV1Jppsk5jm71cw31henrTcNkhNB3HuycExkzFdWLuHp\nh2t4o+kKzRdSt5BMLF6bVheSZbfwgzfPZLopQqSEBGmRUX/8lSrqK328/E4rVwIy3ydmx2418411\nRfzys1tc6wplujkiTRobG7lw4UJKr3n06NHJT8oACdIio0wmxfM7NpCfbecHb52lT/a5ilnaKL3p\nBc3v9+P3+6moqEjpdTdv3jwnA7UEaZFxLruFH/xRPQODcX7w1lkGY4lMN0nMY3armW+uL+bds+20\n3uzNdHNEih08eJCnnnrqruONjY1s3759jFdMTWVlJcFgkGBwbu21lyAt5oRlPhf/4Q9r8d/q48hP\nW0lonekmiXls85pCCnMc/Pnff4KW36UZ2b59O++8886UzvX7/ezevRuTycQzzzwz7nktLS2YTCYe\neeQR/uqv/mrabWpsbKSmpuauY/v27ePYsWP4/f5pX3Okbdu2sXfv3lldI9UkSIs5Y21ZLv/bE+to\nagvw6m8uZbo5Yh6zmE3seKic09d6ePfsrUw3Z94JBoM8+eST/Nmf/dmUzq+srGTXrl3s2rWLtrZx\nKwTT0NCAUoq33nprzN7wZI4dO8bTT99ZTXjTpk0cOHCALVu2TPt6o1VWVs460KeaBGkxp3xtVSHf\n+/pK3vzoOo2nbmS6OWIeW1vq5YGSHA7/+FOZQpmmrq4utm3bRjAY5I033pjSa5qamnj88cfHDdLH\njxv1ijZv3jzjdjU0NODxpLdmeFVVFSdPnpz8xHtEgrSYc7Z/oZztny/jv/7qAh9KZSMxQ0opvvNQ\nOR29EY79VkZmpqqlpWV4UdaBAwem3JtWSlFVVcX58+fHvGZtbS1vv/32jHu8fr+f6urqGb12Ompr\na2loaEj7faZKgrSYk/7k0ZWsK8vlPzec41Jnf6abI+apEp+Lh1cV8lc/a6WrfzDTzZkX1Ijsf9u2\nbUMpNeX548rKSpRS9PTcWdu7q6uLyspKGhoaZtyT7u7upqqqakavDQaD7N69m/r6ekwmEyaTiZqa\nGpYvX35XW30+35hvNDJFgrSYk8wmxX/4w1oKcxy88E9nCcgfWDFD364vQQHP/+OpTDdlzmtsbGT9\n+vV3HJtKb7qlpWU4+FZWVtLU1DT83PHjx9m4cSPNzc0ope66/lRNNNc92es2bdrE/v37aWpq4sSJ\nE+Tm5tLa2sq5c+fuGj73er0EAnMnuZIEaTFnuewW/vJfPogGXm48RyIhq3TF9GU7rfxBXQk/O93O\nZRmVmZAaI4f+tm3byMvL4/Dhw+O+rqmpaXiIvKqqajig+v1+6urqAOMNwGzmo2dq9+7dHDp0iPLy\ncgA2bNiA1vquHvRcJbm7xZxW4HHwvz+xjn/7Nx/w5sfX+ca64kw3ScxDG1cX8pMPr/Ef3zrLoX+R\n/tzb4cE4FzrS+4agIj8Lhy11FeSOHz/Otm3bxnzuwIED7Ny5k6effnrMhVsjg/vIeWm/38/GjRsB\nePvtt9m6deuM2zeTHu5Q4pOHH374rvaOtwBtNsPq6SBBWsx5dZU+dnyxgtd+c5H7S3Ioy8vKdJPE\nPGOzmPjGumJe++0lrnUNUJzrSuv9LnT088cvv5fWe/zNrodYWZyalc7BYHDMXvSQbdu2sW/fPp57\n7jmee+65u573+XzDj6urq/nggw945513hgM0GCuzDx06dNd9jxw5Ql1dHSdOnODZZ58dtw1VVVV0\nd3dP59uiubn5roIoR44cmTDpSSAQIC8vPTXJZ0KCtJgXdm9azi/OtPPyO+f5nx+7H6tZZmrE9Dy8\nagn/2HKV//jWWQ58J7WlLEeryM/ib3Y9lPZ7pMqRI0d4/PHHJ9wjvG3bNg4ePMj+/fvv6IUeP378\njkBYVVXFgQMHeP7554ePDe2PXrdu3R3XfO2111BKsXHjxkkDcGVl5YQ96c7OzruO1dbWcuDAgeGv\ng8Egx48f59ixY+Ne58SJE7PKXJZqEqTFvGCzmHjuyfX8q5ff440PLvPkF8oz3SQxz9itZr6+rpjj\nH1zmZjBEYU7qS1kOcdjMKevlplswGGTv3r3s27dv0nOVUuzdu5cXX3wRgH379nH06FE2b97MwYMH\nqaiooLa2lqeffnp4jvrQoUO8+uqreL1eDh8+zJ49e4avt337dnbu3EleXt6U8mZXVVXR09Nzx5uE\nlpYWXn31VV5//XX8fj/PPPMMdXV1PPXUU8NJVg4fPjycqOTYsWMT7rVua2u7YwQg05SkzJs5pVQt\ncOLEiRNSY/Ye+dtf+vlPb3/Gvt9bPW/+CIq5IzQY59//XQtfX1fEnm+unvLrmpubh4dk5f96arS0\ntNDQ0MCzzz5LS0sLjY2NdwTwsbzxxht0dnayc+fOtLTJ7/fz/PPPD78JGc9kvw9DzwN1Wuvm2bRJ\nxgzFvLLjixUsX5rN0Z+2EhqUilliepw2M4+uLeLvm65wqyec6eYsalVVVeTm5nL8+HEaGxvvSvc5\nlscee4zm5lnFvAm9/PLLkrtbiNkwmxQHvrOBvkiMNz64kunmiHlo85pCLGYTf98kvz+ZlJOTw1NP\nPcW2bdvYs2fPlNN97t27Ny0lJf1+PzU1NSkvgTlbEqTFvFOc6+RffaWaxk9vcjMovSExPS67hQer\nfPyo5YpUyJqHKioqqK6u5sKFCym9bmNj44yKfqSbBGkxLz35UDk5TivH3peczGL6HqrJp70nwqkr\nc6t2sJiajRs3przHOxcDNEiQFvOUw2rmv39kBR+0BTh3ozfTzRHzzKpiD16XlTc/up7ppggxIQnS\nYt565IEiyvNdvPLeRRm2FNNiMim+UJPPWx9dIxaXMpZi7pIgLeYtk0nx7DdXc769jw/a5k5CfDE/\nfHF5Pr3hGO+33Z0EQ4i5QoK0mNfqq/JYX+bltfcvEZUekZiGsjwXxV4nx34j6xrE3CVBWsx7e39/\nDZ29EX52uj3TTRHziFKKh5bn03whQDgaz3RzhBiTBGkx71UWuPni8nz+seUqgzHpTYupe7DKRySW\n4P3zMuQt5iYJ0mJB+JNHV9ITivLOpzcz3RQxjxR5nRR5HfzijIzCiLlJgrRYEErzsvjyCqPKkQxd\niumoq/Dxs9M3ZZX3HNDY2JjyJCXpyE52L0mQFgvGn359JQODcRo+uZHppoh5pLbSR184xoeXpler\nWKSW3+/H7/enPEnJ5s2b53WgliAtFowir5Nv1ZXwkw+vS/ENMWWVBVnkZtn4xRmZKsmkgwcPjpn1\nq7GxcVb1nSsrKwkGgwSD8zO7nNSTFgvKv/xKFT9svsJbH9/gW3UlmW6OmAdMSlFbnkvDJzf400dX\nopTKdJMyyu/38/rrr5Obm4vWmmAwyLZt22hubmbbtm1puWdjYyM1NTV3HXv77bfp7u7G7/fP6vrb\ntm1j7969vPTSS7O6TiZIkBYLyhKPg20PlvEPzVfYcv9SsuzyKy4mV1fpo/HTm5y93ruo65S3tLRw\n5MiRO+opB4NBKisr01oi8tixYzz//PN3HNu0aRObNm3i+PHjnDhxYlbXr6ysnHWgzxQZ7hYLzh99\nqZJYPCFz02LKVhRlk2W38NNFvjtg586ddwXLnJwctmzZktYSjg0NDVMuVTlTVVVVnDx5Mq33SAcJ\n0mLBycu28636Ut76+AahQVnpLSZnMZuoq8jlzY+uLeo88M3NzWN+/7t27UrbPf1+P9XV1Wm7/pDa\n2loaGhrSfp9UkyAtFqTv/k4F4Whc9k2LKauv8nEzGOZ8e1+mm5IxW7Zsoba2lkOHDt0xPLxx48a0\n3bO7u5uqqqoZvTYYDLJ7927q6+sxmUyYTCZqampYvnw5PT09d5zr8/k4f/58Kpp8Ty2YIK2UqlRK\nPauU2qaU2qOUypnCa+5aBaGUekkplVBKxZVSHyil1qenxSKdCnOcfHlFAW9+dJ2I7JsWU7B6WQ4u\nm3lRD3m/9dZb7N69m9dee43q6mry8vImHCI+fvz4uM81NjZO6Z5tbW3TbufQ6zZt2sT+/ftpamri\nxIkT5Obm0trayrlz5+4aPvd6vQQC868Qz0JaVXNMa10PkAzQx4CtY52YDM4+4GWllFdrPfItVyuQ\nA6hRx8U88+8eWcHjP3iXn51u55G1RZlujpjjrGYT68tzefPDa+x8uGbyF0wgEgtzpe9Kilo2thJ3\nCXaLI+XX3bNnD3v27KGnp4e9e/fy+OOP09raOua5E62ET/eWp927d3Po0CHKy8sB2LBhA1prenp6\n0j6/fS8tiCCtlNoADE+kaK2DSql6pVSF1vrC6PO11seTrxtrPb7SWvemrbHininOdfHF5QX85MNr\nPLy6EJtlwQwciTSpq/Dx63MdXAkMUOJzzfg6V/qu8L2f/UkKW3a3F772A6q9s3szMcTv95OXl3dH\ncPN4PBw4cIDXXnttzNc0NjayefPmca851bn9mfRwhxKfPPzww3ccV0qNG6BnM6yeSQsiSAP1wOif\ncgCoAi5M8Lqx3gbmKaUeA4LAFuBlrfX8XLsv+HePruDJ/+uX/OJMO5vvX5rp5og57v7SHCwmxbtn\n2tnxxYoZX6fEXcILX/tB6ho2zj1S5eWXX+b73//+XcfPnz8/nEjE7/fT0NDAiRMneOmllwgGg3g8\nHvbt28fnPvc5HnvsMYLBIEeOHMHr9RIMBjl8+DA7d+4kJ2f82ceqqiq6u6eX7a25uZna2to7jh05\ncmTCpCeBQIC8vLxp3WcuWChB2jvGse5xjk/mpaHet1IqgDFsXj/zpolMKsvL4gvV+fzjyWt8ddUS\nrGbpTYvxOaxm1pTk8E8fXptVkLZbHCnr5d4LDQ0N1NTU3JHxKxgMcvDgweGUmn6/n+3bt/P6668D\nt3vKu3btoqWlBYAnnniC119/naamJurr6+ns7Jw0iUhlZeWEPenOzrsrlNXW1nLgwIE72nr8+HGO\nHTs27nVOnDgxq8xlmbJQgnQ3xhzzSN7k8WkZNTzeBtQqpTwyPz1//cmjK9jxn37FL860s2mN9KbF\nxOorffz1z9vo7IuQ57Znujn3xHe+8x3q6+s5dOgQubm5dHV1oZTi6NGjw8PHGzdu5OjRo2zZsgW4\nPR/d1tbG5s2bOX78OFu3bsXj8dDV1YXH4+H8+fNTyuBWVVV111xyS0sLr776Kq+//jp+v59nnnmG\nuro6nnrqKSorK9m1axeHDx8eTlRy7NixCeei29ra0rpKPV0WSpBuAp4edcyHEWQncsekSXJuu1Fr\n7YPhue1JJ1a+973v3TWcs2PHDnbs2DHZS8U9UFHg5gvV+fyo5RpfXrFE5qbFhDaU54KCd8+08636\n0kw3557Ys2cPAOvXT7yZZag32tLSQm1t7XBO7PGC43PPPXdXcpSxPPnkk7z66qvs3Llz+NiGDRvY\nsGHDHT3mkcbK8z2edO7FfuWVV3jllVfuOJbKRXMLIkhrrVuUUsND28nH50cMW28AuseYWx79Fq8N\n+IsR13kcaJisF/3CCy/cNT8i5pY//fpKvvOXMjctJpfttLKyyMOPW64umiA9Vbt37+all15CKcWW\nLVvwer089thjgJEf+/Dhwxw/fpyGhgb8fj/f//73p5Sp7LHHHuOZZ55JW7tffvll9u7dm5Zrj9Uh\na25upq6uLiXXXxBBOukJpdQewI8xh/zEiOf2A+8DhwGUUpuAWoye9H6l1Nta63eSPeeW5HWCGAvP\nRl5HzFPl+Vk8VGPMTX9lpfSmxcTqKn288t5FekNRsp3WTDdnzli/fj0vvfQSb7zxxnBwHmnPnj00\nNjZy8ODBaW+D2rt3L0ePHr2jN50Kfr+fmpqatKY1TacFE6S11ieBoV33x0c9t33U141AI3BojOsM\nPScWmD99dCVPSm9aTEFdpY+//dUFfvnZLb6+rjjTzZlzvN7x1+RONPw9kYqKCqqrq7lw4UJKA2pj\nY+O0hsbnGulOiEWjbERvOhpPZLo5Yg7zZdmoXuLmR83pTUgyX020AGusHvZ0rpvqHu98DtAgQVos\nMv/ukRV09Q/y3rmOTDdFzHH1lT4+vtzNQCSW6aaIRUyCtFhUKgrcbCjP5Z8+vE5iEVc7EpOrr/IR\njWvea5U3dCJzJEiLRWf35uVc6w7x0aVpb6MXi8gSj4PyfBf/0HQ5000Ri5gEabHorCvLpabQzY9P\nXst0U8Qc92BlHh9e6iYSk0pqIjMkSItF6emNy/nsRi+tN6SWihhffZWPSCzBJ5dl1EVkhgRpsSh9\n6b4CirwOfvKh9KbF+Iq8Tkp8Tn72aXummyIWKQnSYlEymRR//JVqmi90cb07lOnmiDmsvjKPszck\ndb/IDAnSYtH6+rpiPE4rb310PdNNEXNYfaWPSHRxz0kHg0H8/qlX7G1sbGT37t3s27fvnt53IZIg\nLRYtm8XEkw+V88vPbtETima6OWKOKvE5yc2yZboZGdPS0sKmTZvYunXrlF+zadMmqquraWycefLG\nmdx3IZIgLRa1bQ+WopTinU9vZropYo5SSrGqOGfyExeoDRs2DNeUno7ZFh2a6X0XGgnSYlHLcdn4\n3Q3LaDx1k8GYpAoVY1tVnJ3pJohFSoK0WPS+81A5vaEov5ZUoWIcpb6sTDdBLFILpgqWEDNVlpfF\nhopc3vroOl9ZWYBJjS4zLhY7tcC7M36/n4MHD1JfX09rayvV1dX4fD4aGhp48cUXx33doUOHqK6u\nRmuN3+9nz549dzyvteaNN95Aa00gEKC7u5tnn332jvs2NzcD0NTUxObNm9m0aVN6vsl5SoK0EBjJ\nTXb/9ft8fLmbdWW5mW6OmOcSoRCx1ta03sNSU4PJ6UzJtZ544gmef/55Nm7cSDAYpK6ujtbWVurq\n6sZ9zfbt29m9e/dwRSy/38/WrVv553/+5+FzWlpaqK2tHa5sdfToUXbv3s1LL70EwK5du9i6dSt7\n9uxh27Zt1NTU0NDQMG9rP6eDBGkhgHVlXsrzXbzzabsEaTFrsdZWbj36jbTeo+DNn2B74IGUXMvv\n9+Pz+QDIycmhra1twrrOfr+fxsZGXnvtteFjlZWVBAIB3nnnneHAXVVVdcc1du7cya5du9i3bx8V\nFRUcOXJk+L5D5zc3N0uQHkGCtBAYK3i/84UKnv/xp3T2Rchz2zPdJDGPWWpqKHjzJ2m/R6ps376d\npqYm1q9fT3NzM9XV1RMGyoaGBqqqqu46XlVVxdtvvz1hvWmv1zsciCsqKvD7/bz22mvk5uYSCAQI\nBAKp+JYWDAnSQiRtXVvED946wy/OtPPt+tJMN0fMYyanM2W93Hth8+bNKKU4evQowWCQEydO3JP7\nPv/885w4cYJXX30VYPizuE2CtBBJWXYLj64t5hdn2vn92hLMJllAJhaHtra2OxZ0Taa+vn7MbGLN\nzc3s2LFjwtcOzXkHg0H27dtHInF762N3t1HI5OTJk6xfv37K7VnIFviaRSGm51v1pQT6B6XWtFhU\n2traJky/qbWms7Nz+OsNGzawefNm3njjjeFjzc3N5Obm8u1vf3v4WCAQoKfndt7zo0eP8vTTT1Ne\nXk4gEEApdcfzXV1ddHV10dbWNuZ9FyPpSQsxwspiD5UFWfz09E02VMgCMrE4PP7441RXV6OS2w+r\nqqrYvHkzL774Ii0tLezbt49gMMj+/fvZv38/Ho+HV199lf379xMIBNBa09bWxgcffHDHdY8ePUpD\nQ8Pw80qp4S1dlZWVHDx4kL1797JlyxYAjh07xq5du8jNzR33vouNBGkhRnnyC+Uc/NGndPZGyMuW\nBWRiYWtpaeHYsWP4/X7Ky8vp6ekhEAhw4MAB9u/fz3PPPXfHtqqRnnvuuXGvO5X9zqP3VQO89dZb\nw4/Hu+9iIsPdQoyy9YEibBYTvzx3K9NNESLtmpqaqK+vp7y8HACPx0NFRcVw71ZklvSkhRjFZbfw\nYJWPd8/c4vc2LJMMZGJB27lzJ0ePHmXfvn3k5+eTk5PD+fPnyc/Pn7CnLO4NCdJCjOG7X6pi91+/\nz+mrPawpWbwVkMTisHPnzkw3QYxDhruFGMO6Mi9FXgc/P9Oe6aYIIRYxCdJCjEEpxeOfK+eEP0Bf\nOJrp5gghFikJ0kKM4xvri9HArz6TEpZCiMyQIC3EOHKzbNRW5PLzM+1orTPdHCHEIiRBWogJ/Isv\nVnC1K8TFjoFMN0UIsQjJ6m4hJvBgVR4ep5Vfn7tFRUFWppsjMuz06dOZboKYA+7l74EEaSEmYDGb\neHRtEW99fJ3vfKEckxTdWJRyfXk4XS6++93vZropYo5wuVzk5+en/T4SpIWYxCNri3jtt5f49FqQ\n+0u8mW6OyIDiZaX8+Ge/pStwZ7GHHzVf5UJHP//HH9ZlqGV3unz5MgClpVJqNd3y8/MpKytL+30k\nSAsxidXLcijMcfDeuQ4J0otY8bJSipfdGfx6nSW82NhK2fI15M+BPO9DBShqamoy3BKRKrJwTIhJ\nKKX4vQ3LaPIHiETjmW6OmENWFhtBseVCIMMtEQuVBGkhpuCRtUWEowlaLnZluiliDvG6bBR5HZzw\nS5AW6SFBWogpKM3LorIgi/fPL+4C9OJuK4s9/Pa8JLwR6SFBWogp+sb6ZXx8JShD3uIOK4o8XO8O\nE+iLZLopYgGSIC3EFD28upDBWIKPLndnuiliDllemA3AJ1eCGW6JWIgkSAsxRSU+F2V5Lj5ok/lH\ncVue24bXZeUTefMm0kCCtBDT8Oi6Yj681MVgLJHppog5QilFTWE278u8tEgDCdJCTMPG1YWEowk+\nuSK9JnFbdaGbtlv9xOLy5k2k1oIJ0kqpSqXUs0qpbUqpPUqpnCm8ZlsqriMWj4oCNyU+F79plVXe\n4raaJW4GYwlab/ZluiligVlIGceOaa3rAZKB9RiwdawTk8HZB7yslPJqrXtmch2xOH1jfTH/9y/a\nGIwlsFkWzPtcMQsVBW7MJsUnV7qHE5wIkQoL4i+MUmoDMFzwV2sdBOqVUhVjna+1Pq61PjryNTO5\njlicHl5dSCSW4JSs5hVJNouJEp+L01fld0Kk1oII0kA9MHrJbQComuR1o0sazfQ6YhGpLHBT5HXQ\nJFmmxAjVS9y0XJCMdCK1FkqQHqvqQfc4x+/FdcQCt/WBYloudslCITGsakkW17tD9IWjmW6KWEAW\nSpDuxphjHsmbPJ6J64gF7uHVS+iPxDh7vTfTTRFzRNUSNxo4c61n0nOFmKqFsnCsCXh61DEf0DbJ\n6/Sor2d0ne9973vk5Ny5CHzHjh3s2LFjktuL+WpFkYc8t40T/gBrSmQDgIAirxOH1cynV4PUV+Vl\nujniHnnllVd45ZVX7jgWDKZubcKCCNJa6xal1PCQdPLxea31heTXG4BurbV/1EvvmJOe7DrjeeGF\nF6itrZ3dNyHmFaUUW+4v4s2PrvHdL1VgUqOXN4jFxqQUVQVZ/Ka1g//my7KMZbEYq0PW3NxMXV1d\nSq6/UIb2fRSqAAAgAElEQVS7AZ5I7mveBuwFnhjx3H5geE+0UmqTUupZjJ70fqXUxileR4hhX121\nhO6BKBdu9We6KWKOqFzipq1d9kqL1FkQPWkArfVJ4GTyy+Ojnts+6utGoBE4NJ3rCDHSA6VeXDYz\nH13upmqJO9PNEXNAZUEWPz55jc7eCHnZ9kw3RywAC6knLcQ9ZTGbWFOSw0eXZF2hMFQWGG/WzlyX\nxWMiNSRICzELj6wtpq29j56QbLsRRkUst8MiSU1EykiQFmIWvlCTj0ZqCQuDUorK/Cya2iS3u0gN\nCdJCzEJ+tp3y/Cw+uiSZpoShosDNhQ5ZTChSQ4K0ELP08OpCProcJJEYve1eLEaVBVl0D0S51RPO\ndFPEAiBBWohZemh5Pv2RGH7pPQmgoiALkMVjIjUkSAsxS2uW5WC3mDhzTealBfiyjMVjn0mQFikg\nQVqIWbKYTawo8vDpVfmjLIzFY+V5WTRLlTSRAhKkhUiBr65awmc3eqUqlgCgLM/Fxc6BTDdDLAAS\npIVIgbrKPAZjCc5LSkgBlOdn0dEbkf3zYtYkSAuRAsuXZpNlN8uQtwCMIA1w7oaUMhWzI0FaiBQw\nmxQrijyyeEwAsDTHgc1iksVjYtYkSAuRIl9bVUjrzT4i0XimmyIyzGRSlPpcfCY9aTFLEqSFSJHa\nSh+xhKb1psxLCyjxuTh1RYqviNmRIC1EilQVuMl2WDh9TYY4BZTmubjWFZIV/2JWJEgLkSImk2Jl\nsYfTMi8tgFKfi1hCc1m2YolZkCAtRAp9dVUh/lv9hGVeetEr8bkAaL0p89Ji5iRIC5FCdZU+4gkt\nW28EboeF3CybrFEQsyJBWogUqsjPIsdplXlpARhD3icvSnpQMXMSpIVIIaWMeekzEqQFsCzXybWu\nUKabIeYxCdJCpNiXVhRwsaOfwZis6l3slnqddPRG5HdBzJgEaSFS7IFSL7GE5qLUl170luY40MCV\ngKzwFjMjQVqIFKspzMZmMXFOVvUuekVeBwCXOuUNm5gZCdJCpJjFbKJqiZvzsqp30fM4rThtZi7J\nqIqYIQnSQqTB56vzOXezF611ppsiMkgpxdIcB5ckoYmYIQnSQqTB2jIvwYEoHb2RTDdFZNjSHCdn\npRqWmCEJ0kKkwf0lOQCSyEJQ5HVwo1u2YYmZkSAtRBrkuGwUZNu5IHORi94Sj4PecIz+cCzTTRHz\nkARpIdKkoiAL/y3pSS92S3KMFd5XumReWkyfBGkh0uTBqjwudgyQkMVji1qhxw7IXmkxMxKkhUiT\nFcUewtE4N4PhTDdFZFCW3YLLZuaqBGkxAxKkhUiTlUUeAC7cknnpxUwpxRKPQ3rSYkYkSAuRJjku\nG/myeEwASzx2PpPypWIGJEgLkUYV+VlckMVji15hjoP2Hpn2ENMnQVqINKqt8HE5EJLMY4tcfraD\nrv5BolINS0yTZaYvVEqtB6oAX/JQG9Cmtb6QgnYJsSBUF7rpj8ToHoiSm2XLdHNEhhRk29Ea2nvC\nLPO5Mt0cMY9MOUgrpTzALuBJoBLwAwGgO3mKF/AppSqBJuAY8JrWWvLhiUWrujAbMLbfSJBevPKz\njW1Y17pDEqTFtEwpSCulngW2Ay8DT2it/ZOcXwlsBt5RSr2ktf6rWbdUiHmo2OvEbjFxOTDAA6Xe\nTDdHZEie24YCrndJelAxPRPOSSulcpRSrwEntNYPaq3/arIADaC19mutj2qt6wG/UurFZE9ciEXF\nZFIs8zll+80iZzGbyM2ycV1yeItpmmzh2HZgp9b6nZneQGvdCOxLXkuIRWfNMq8EaUF+tp1rEqTF\nNE0YpJO94eBsb6K1DsqQt1isqgvdXOsKkUjICu/FLD/bzvmbsldaTI9swRIizSoK3ETjmo4+qS29\nmOW5bQT6BzPdDDHPpD1IJ+e1dyqlNqb7XkLMReX5WQAyH7nI+dx2uvsHicuIipiGtAfp5FD3UaBa\nKXUuXfdRSlUqpZ5VSm1TSu1RSuXM5Fyl1EtKqYRSKq6U+iC5H1yIGSv0OLBZTNzoloxTi5kvy0ZC\nQ2evjKiIqZtxMpOxKKWeGm/uWWt9VCn1RCrvN8qx5GpykkH3GLB1Bue2AjmAkj3eIhVMJkWhxyE9\n6UXO5zb2yd/sCQ/XmBZiMqnuSW+Z5PnmFN8PAKXUBmB4DCm52K1eKVUxg3OV1rpXArRIpRXFHlnZ\nu8jluY2EJlK6VExHqoP0EyOGiZ8bYx66NcX3G1KPkf1spABG2tLpnpunlHpMKbVJKXUgmZhFiFmp\nKnBztUtyeC9mLpsZu8UkQVpMS6qD9F6gBjgCVAOvjwjae4C6FN9vyFipnLrHOT7ZuS9prd9I7u9+\nFWMoXIhZqVzipi8coycUzXRTRIYopfC5bbT3yIiKmLqUBmmt9aER2ca2a619wIPAaxhzvk+n8n4j\ndHO70McQL7fzik/53FEFQtqAWsmWJmZraIX3DelFLWq5WTY6ZOGYmIZxF44ppTwTzctO9vwQrXUz\nxlz0oWQO8HRo4u43AD6MIDvlc5Pz1Y3JNxdorYNKqUnHJ7/3ve+Rk3PnYvIdO3awY8eOKTZfLHTL\ncp0ojPnIFUXyni8ttIZYDAYHITqIikaNr+MxiMVQsTjE48bX8TgkEsZrkp/ViMeYkv0XpUCBRhmP\nzSYwW8BsTn5Y0GYzWCxgtaJtNrDawGYDq9V4zQhel41LnZJ9biF55ZVXeOWVV+44FgzOOgfYsIlW\ndz8NHJ7F82M5Ms3zp0Rr3aKUGh7GTj4+P9QrTgbf7mQvf9xzlVJdwF+MeO5xoGGyNyMvvPACtbW1\nqf2mxIJit5rxuW3Sk56uSATCYVRoACIR1GDE+BwJYxroh0gYFTGCMtHYxNcym4xgmgywWpmMYKyA\nkY9RgDaWl2p9+wMN8QQqPjrYT/A+3mYFmx3tdJDI8rC8vZv+9gixi8tQTifK7TY+jwrmYv4Yq0PW\n3NxMXV1qZncnCtLPK6WenOD5KqYQpJNbnLZjBMIZ5wCfgieS895+jMVhI7d77Qfe53Z7xzw32XNu\nST4XxPge07ltTCwihTkOWTQ0UiQCoQFUf78RhAcGUKEQpr4eVDgM4dDdAdBqMYKe3Q52O3i8JOzJ\n3qvVers3a7EaAXI4KFvu6tWmTCJhBOtoFAYjEI2iBgeTXw+iIhGIhDH1Blna086SSzcJvzPiDYXZ\njCkrC+V2Y/LlYcp2o7I9mLw5mLKy0tNmMW9MFKSPAl3AB2M8p5ji/HJyi9PRZNaxl7XWy6ffzCnd\n5yRwMvnl8VHPbZ/GuY1AYzraKBa3FUU5NLV1ZroZ947WRuDt60X19UFfL+aeblR/H4RCEIvfPtds\nApcL7XCC200iPx8cTnA4jGMOhxGUzebMfT/jMSV74VYruIxa0SPfWox8PFDSxdv6Av/+m3Vk6xiJ\n/n50bx+6r49EXx/xSxeJ9vUZgR9QVpsRrJcswZSbi8nnw+T1oubiv4NIi3GDtNZ6V7IXvBnj9+yO\nYd/k0PAdMpzMRIg5rTTPxY9PXkVrvbCGN2Mx6OtF9fSgenswd3Wi+nqhv/92T9ikICsLneVGLylE\nO11GUHa6wOUE++JI7uF1WgEIDEJOfg6mnLsTI2qt0X39JHqC6GCQRHeQ+JUrRD/9dHi+3JyXh7l4\nGaYlBZjz81GOxfHvtxhNmHEs2Qs+DpDcN5wDtGmtTyZ7nKNtASaqdpWWZCZCzAeleS4GYwm6BqL4\nsmyZbs70DQXjYNAIxoEOVG8PDIzYUuR0oN3ZRiDOcqPdbsjKAqfr9mKsRSzXZQTpW72DVOaPPZSt\nlEJluzFlu2HZsuHjOhYjEewh0dFBouMW0c/Ooj80BgRNHg/mZcuIhcOYCwrS/42Ie2bKaUGHgrJS\nqkIp9SLwttb6jVGnPZFcbNUMNCTPGTkPna5kJkLMectyjaHQjt7I3A/SAwOoYDcq2I25ox3V2wv9\nfbfHbl0utMdDoqQU3NlodzZkZxtDvmJcHqfxJ7drBtWwlMWCOc+HOc8HK+4DINHXR6Kjk3jHLeLX\nrhH57DMABlpbsVRUYC4qwpSfj5I3SPPWlIN0cq/w08BujD3Fb49x2l7gdYwh8i3A0JB5M0ZikGom\n7mkLsWAtTeZr7uiNcN/S7Ay3ZoS+PlR3ANWdDMjBbhhMJl2x29CeHHThUnR2NtqTA55sY2GWmDaH\n1cg6FhhITVIbk9uNye3GUlFuXL+slMStDkw2G9FPPmGwuRlls2FetgzL8vswFy2VgD3PTBqkk1Wg\nvo8ReF8Ftmit/WOdq7U+lHx4NPmBUqoW2ISRzGQT8Mzsmy3E/OOyW8iyWzJbBSkcRnV3oQKdRkDu\n7oJIslfncqJzvCSqayDHi87xgtOZubYuUB6Hha4UBenRTE4nprJS7NXVaK1JdAaIX7tG/NIlwn4/\nyuHAUrMcS2UF5vz8tLRBpNZEyUyewugZAxwcvUJaKfXYGMPdd7lHyUyEmBfys2109N2jIB2LGUPW\ngU7Mt26iugK354/tNnSuj0RlNXi96NzcRbN4K9M8TsuMhrunSymFOT8Pc34erH2AeCBA/OIlYudb\niX7ysdEDv+8+LFVVmLLn0MiOuMNEPekjyY8TgEoG7SEK+DNg0iA9xjWFWLTy3fb0pYUMh1Adt1Cd\nHUZQDgaN1dVmE9qbS6K4BHJz0d5cYzGXyIgcp5WrXfc+f7fZ58Ps82Fdv45Eezuxi5eIfmwMiZsL\nC7Hefz/m0tKFtfNgAZgwmYnWet94TyqlJitLeZfkanEhFq3lSz38/MzN1Fysv98Iyh23sNy8bmx5\nAnC70T4fibIKo4fsyZGV1XOIx2HhUiBzRTaUUpgLCzEXFqLraolfvkKstZVwYyMqy41t7QNYqqtR\nsghwTpgoSL88yWv3jj6glNqY5qxiQsxrS70OOnojM9sr3ddrBOTrV1GdHRAyspdpjwddWIjOy0fn\n5RlJQMSc5XFYCM6RamjKbMZSUY6lopx4Z4DY2bNEfvMbBptOYF2zGuvKlShZl5BREyUzGXNx2CTP\nH8SoeiWEGMMSj4NoXNMfieF2TNJT6etF3WrHcuMaqqMDwmGj2IPXS2JZKQwFZdsc384l7pDtsNAX\niU9+4j1mzvNh/uJDWNetJXauleipT4l+8gnW1WuwrF6FSYJ1Rky0cOxZjNzV57XWh5NbqY5h1IRu\nAHaOUXiiK5n3+vVRJR+FEEB+th2A7oHo3UE6EkHdake138Ry/YqxyEthzCeXlkJ+AdqXJ3uRRxtZ\nBEMPbeROFshQ6nbazjnCbTcTiSUYjCWwWeZOu4aYsrKwrV+HddUqomfPEj19muinp7CuXYt1zRqU\nZco7d0UKTPSv3YaRXWwot3UD0KS13prcVrU/+TFMa70VhrOT1U5l9bcQi0l+trGCuqt/kBKvAxXo\nRN28gfn6FVSXUf5ce7LRRcvQBQXo/PyFvSc5HjcKbSQrWzFUmCIaRUWNx6ZQv1HhKh5LVqFKVp8a\n+tATVKEaMjJYK5NRBctqlJfEYiXhcILNhrZYjZEJhwPtcBg5wx2OlP4M3Hbjz25POEq+256y66aa\nstuwrX0A64oVRM+cZvDkSaKnT2N78HNYqiplgdk9MlGQ1kNBNtmLrsPY54zWulkpNW4dLq11o1Iq\nJ9mrPjKVutNCLAa5OkJJ1zWsv27HZosYqTbtNnTBEhKV1eiCJQtjb7LWEAnDQAgVDhkFNULGY1Nf\nT7IylFEx6i5mC1gtaKstWcnKSsKdTKBiNhv1m0fWdTaZble4UurualeJxO1KVVoP15JWMePNgCk8\ngKmvF2JRowZ1ZNB4QzCSxQIOJ9rlJJGTi84y0p1qtxtcWdPqqbsdxp/dYCg2p4P0EGW3YVu3Dkt1\nNdEPPyTyi58T/eRj7F/+MmafL9PNW/CmOm6xGWPYe2SwnfDta3Il92Gl1DalVJcsKBOLkY7HSdxs\nJ3b1CvELF0j09bEu0EaooJTE51eilxSC15u+MorpFI8bK8z7+4Y/m3q6UQNG2cmhSk6AEUwdDrTT\nRSI7B52frGpltxs91qHSk1bbPat0NfQHLDHWk7GosTAvHDZ6+OGwUUaztxvTjeuogX7j+wfjZ5dM\nkxr3FUy67zzLbnx/PXNk8dhUmdxu7L/zO1juu8Vg0wlCP/whtgfWYl2/TqpypdFEQbp6xONdGMPd\nE1JKVYyei9ZaH0/m+34OeE561WKh05EI8avXiPnPE79+Ax2NorKyMBcXYV26lE8CZURLcqlfUZLp\npk7NUGGNvj7o6cHc3ZnM5d1/e6jZbEFnuSAri9iyMqNko9OFdjmN1eb2ud9jvIPFCtlWyM5mxCz3\n7YCutTE6MPQmpa8Pc+AWlnNnjZrSYHz/ubnElxQaIyTZHlCKbPvtnvR8ZC4owPHIVmJnzjD4ySfE\nLl3E8fDDmLzeTDct5SKxMD2DPfQO9tIb7aNvsJe+kZ+jfQzGo0QTgwzGB+mL9hJNRLl25lrK2jBR\nkG5QSjUBORjJSx4HY74ZI5HJsTFec1Ap9RJGgK/DWHhWD3iT19iMrP4WC1Cit5f4lSvE2tqI37oF\nWmPy+bCuWoV5WfEdf8Cy3bfSlhZy1vr7UT090NONOdCJ6glCX9/t551OdHY28eJSI5e32w1u9+Lb\n9pXsPeNyoZNVp4ZDbn+/kW61qwtzx00sH39o9LptdnRBAVlLinAOhuZdT3okZTJhXb0ac3ExkV/9\nmtCPfoTtoYew1tRkuml3iSfiDMQGGIgOMBDrp2/QCK790f5kwO1LBtxebg20DwffvsE+oom7f0YK\nhdPiJMuahdPixGa2YjFZsZqs2M123FY3g87U/Wwn2oLVAtQrpSpHbbcKAOMlOXki+dGGkQq0BSPL\nWFvyekIsCEZe5E7ily8T8/tJBINgNmMuLMRWX4d52bJxt6yU+Jz4b/Xf4xaPojX09hppQ7u7MHXe\nMoLz0ByxzSisES8uQWd7jMIaUuVqarKy0FlZsKyEGA8YAToQQN1qx3zjKvaTTTx88TrOhqtEYuuw\nlCzDtGTJvCx8YfJ6cTyylcGWFiLvvkvswgUcX/4yapYjJwmdIBwLE4oNjAiwIUKxAUIjHg/EBgjF\nQoSiIQZiAwQjXYTjEcKxEKFYmHAszGBi/BSsdrMdl9WFy+LCZXWRbc2mwJWP2+omy+bGbXXjtmaR\nlfxwWV04LU5MauKf1cfxj/lL/vOs/g2GTDonPXo/9CTBdsIsZULMZ1prI53ipUvEzp9Hh0Ioux1z\ncTHWtQ9gLiycUpYmr9NKb/geD3VGwqhAAAKdmG+1o7q7by+OyspC53iJrSwFj0cKa6Sa2QwFBeiC\nAmKr10A0yrnOX7HCayJ+8QKx06fBZsVSVYWlsnLeFb5QFgv2Bx/EvHQpg7/9gNAPf4j10a2E7Sb6\no/3Jjz4GogP0R/voGz5mHO8MdxjBNhloQ7EQ4VgYPcGyJ6vJisNsx2Fx4rDYcZiNz9k2DwUWBw6L\nA6fFaXw2337ssDhxWZzDgdlsmvtz6VOpgrURY+vVpHPJEqDFQqO1JtHRQezCRWLnW43A7HJhKSvD\nXFpi1Oqd5qKvHJeV3nQns+jrRXV2GsU12m/cHrJ2OtHeXGJrHjByeHu90ju+16xWQvmFXCnz4dpa\nYxS+uHyZWGsrsTNnUe4sLMvvw1JdlZEEIgmdIBQzAmZ/srfaH+0ffjwQvX18IDZAIBIwerOxEAln\nH5Uf3iDe8gKf3J/NgOvuEGM323GOCJZOi4tcey7F7mJcFifOER92swOn1YnT7MBuceC0OHCYHfMi\nuKbKVFZ3K6BbKXUeY/HY27L/WSx08Y4O4hcvET3fiu7vRzmdWEpLMZeXYcrLm9UeUa/TSl8kRiKh\nMZlSsKpba+gJojo6MN+8jgp0GtubAJ2TQ7xoGdqXZyRCcblmfz8xay6rmb6IMZIxXPhi7VoS7e1E\nP/uM6EcfEf3wJOaSEqxr7jcqWU1Ca000HiUUD3Gt/zr9sX5CsRD9yd5pf7SfUDxEf/R2gO2KdBGK\nh5JDy2HjcXz8AjAmTMkAavROh3qpS5wFxuN8J+5SC0UnLrCxUxFd8zlsBUuMQGtxzpve61wyleHu\nRqXUIWADxjyzGmsVtxDzXaKri9iFC0TPtaL7+1AOB+bSUixlpZgKClKWvMHrsqI19A/GyXbMMHtT\nsBt16xbmm9dQnZ3GXLLJhM7NJVa1HPLyJDvZHOa0mWnvuR0MtdaE4mH6PWb615fRf5+HyPnzJM6+\nR/TDtxjIttJZkUegwEF/fIDuSBehuDHnGoqHCceNINtzzRjwzLp4d5Uzq8mK0+wYMQTswGl2kmPL\nxpEMtkPHnCOGiUces5lsU/p/oMsG4Zfvo5uuoR+uQOUUpOzfbrGZynD3Y8BvtdZ3FdQYcc6zGMlP\nDqeycUKkmw6HjcB89gyJQBfKZsNcVoql7EFjMU8a9i97XUbg7A3Hph6ko1EjZejNm5ivXTbyeJst\nRlBesQrty4fcXCPphrjn4jpOOB4iFB8gHB8gFB9I9kqNz8bx28935XYQ1RH+9c9iyaHl0NhzsEs1\n+Y4EFTfi5P9ck3A7GVheiCrOw2fPxem6HXQdFgfd9m7sFjvVVVXDgXboeYvp3v1uKJsN/ZWHUL/+\nAH76K/TWr6J8uffs/gvJVH5qW7TWz0x0gtb6EIBS6gDwF7IXWsxlOpEgfvUq0TNniF8z9jOai4ux\nr1mDubg47atsvc7bQXpCvb1GytCrl40h7EQC3NnEyqvQhYXgy7tniT8WslgiSjgeIpwIDwfV8IgA\nG0mE7jg2kAgO916Nj9CYW3WGGD1YJw6zA4fZidPsxEYW0bCX+qolxjysOTlHO2I+duhrq8n4fdGB\nbtSHp+DmLVB56AfXoXI8d9zLH74IQKW3PH3/YFOkLGb0Fx9E/fw9aPwl+tGHUdnuTDdr3knpWyut\n9b5kKlDpUYs5Jx7oInb+PLFzn6EjEUy+XGzr12OpKJ/1lpHp8DiN/3Z9g6MWjyUS0NmB6cZ1zFcv\nG8lCzBZ0fj6x9XVGdrKsu4cxF6tYIkY4ESISDxnBMjEUSENEEuHhx0YANj4PJHqIxMPGNp14iEg8\nQkyP/2bJoizDwdVhtg8/9lhzko+HnnMMDxmPPm5Wd7+ROvZ+NyfaQzz2zakHU+XzwsO/g77Rjmr6\nCPVPP4XKMvT6NSj73KyEpixm9Jc+Bz/9FarhXfTXN6Ic8yyxTYZNJUhPNznrcaXUY7K4TMwFOhYj\n5vcTPXWKRFcXyuHAUlmJpbIiYxmSspIZp0KDRoEI1X4T80U/6uYNI7uX00V8WQm6cCk6v2BBDmFr\nrRlMDI4aHh64Y1j49kc//fFg8ng4GWTDEwZXEyZsw0HV2KJjN9nJtmRTYC/AnjzmMDuSj+8OxHaT\nPW1DxHaLIhIdMyHppNTSJehvbITWC6iPTqOuXEM/uB5VtizFrUwNZbfBV76AbnwX089+TeKRr0lx\njmmY0upupZRnqkPYWmu/UmrbLNslxKwkenqIfnaO2Jkz6FjUGM5e+wDmpUsznjTCZoIlA11knWzH\n+mk/RKNoj4fYytXopUWQM/fTKw4NEYcSoRFDxEMB1giykWg/uiuAKRDA1NWBvasXe1cfjmAIS2gQ\nRyiOMxzHGUqQG45TFE5gjSYwJTTmBJgSoLTGlABTQhO3mEhYLSRsZrTVCjYLympFOeyEigrRvlxU\nng9TXj7Klwc+H4miJWhf7pzLjW6zKCKxKVTvGocymeC+KnTZMlTTh8bc76Wr6CV5c7KUpHI54aF6\n9E9/BafPwer7Mt2keWMqP82XgaPAk9O47uT7BYRIMa018WvXiH78MfHr11F2O5bl1VhqajDNgWHi\n+NB+a38bn7/2GdgLiH15PXpZiZHN6x5I6MTwUHAkHjZ6rokwkXjodqBNzs0a860DDCR6CI/owQ71\nYlVCk9sVJa9zkPzOQfI7BsnrHKS00ziW3RPFNCoODWbZiGU70Q47uNyYnA5MuU4sTidWhwurzYFF\nWVBmsxGITApMZuNzNIaORtHRQfRgFD1oVLFKhMJkXe8m8dkVEj096P6BO8pXKocd85ICBu6rJlG6\njETZMhKlJcRrKiFDQ68Oq4lITKO1nlWvUjns8KXPoS9dRb3fgjrbin5wfQpbmjoq3wcra+DkKXRR\nISo3J9NNmhemugXroFLqX2ut/8sUr1s1y3YJMWV6cJBY63mipz4h0deHyZeL/Qufx1xWlvHqPIne\nPmJ+P7HWc0bwcDiw1izn1KpcEitL2LBy6kOU0USUULz/9jBwbOCOedjwcPA1HofiQcKJyPAcbCQe\nJpIYfw+sQmEz2UcM+Tqwmx1kmVxUBiwUXtb4rkbwXo3ga7uB/WY3KnZ7Xl15PZjz8rAU5GFelYfJ\n68WU48Gc48HkycbkyZ5SRrbZ0vE4id4+I596R4D4zXbi7bfIuniF+ActJILJQUGTCUtxEQN1a4mt\nuo/4+vtJlJfek1633aLQGqJxjc0y+/upsmVobw68/iPUL3+LzvOhli1NQUtTbPV9qOvtqF++T+Ib\nGzP+/3M+mOq4yHagSSlVrbX+/kQnJmtPz62xJbEg6XCY6NmzRD/+BB2PYSktxfbQQ1NK/JDWdmlN\n/Pp1Yx78xk2wWjFXVqJLixjIzaIrHiL24QVuxs7Q0nXl7vnY2AD9ie7hOdhwzOjljjcHq1DYzXbs\nJmP+dXiO1eIkx5SbPJZ8zuQYEYDtw4HYYbZjM9kxKRMq2IP51BnMH5/G814TUf9FdDI5inJnYSla\niqV6JZYvFWJesgRzQR7mPN89CcBTocxmzN4czN4crKV3VxpLhMPEb7QTvXiJ6IVLuJo+JPbDN42i\nKN4cwl/6PNFNXyH24Ia0rQewWY0/keGoxpaiWyiPG/2lz6FaPkG9+xv0Q/Wo8rlVaU2ZzejPb0A3\nvAtnz8uw9xRM6ddDa92mlNqMURnrCeAAcGz0PLVSqgKjOta4e6qFmC0dChE9fYbop6dAayzV1VhW\nrffBr88AACAASURBVExbCsWhRU4jMzWNTJk4kDweCvVgvnAVu/8a0d4gwSzFlSIr13wm+rrDxLtG\nrOZeCpeA/3rBCLIOsyO5TWco37CTbKtn+LgzuU1n6PHQVh6HeeoJJsb55lDXbmBp/i2Wj07hbDpJ\n/PpNAEzZ2aiqCrK+uRVLWSnWkmWYPPdmWD6dTA4HpooyrBVl8FXjWCIcIdp6nsEzn6E+OIntH/8Z\nlZVFZPNXGNz8VeIb1oIldb0+h8VYFxGOJvA4U3ddZbGg69fBrU7Ue01omHuBOscD1eXw8Rl09b3d\nWTEfTfk9nNa6WSlVhRGEjwJHlFJD1a4CGEPcmzGKbLyTjsaKxS0xMED01KfEzpwGpbAsr8G6YgXK\n4Rj7fJ0gHI8M5xU2gutQ9RyjdzowfNxIndgVCRCORQjFw0ai/7iR0Smhx1+JmxsyUX4jTlFnDBtW\nAkW59D+wHOXzUWZxsiKZQtE1FGwtDv7uVwMMDlr4txuXYTPZ7+lqVyMof4Sl+UMcvzlBItBl/Hsu\nK8Z6Xw3WR7dgra7EXDD9vOTzlclhx37/auz3r8a97Q+IXblK+INmTL96H9s//BOm7GxCW79G9He3\nEl81+97f0BD3bBaPjUcphf58LQrmbKBm1X1w4Qqc+gxqH8h0a+a0aQ20aK27gS3JXvUuYBNGacpu\noAnYqrVuTHkrxaKU0AkGogP09gcInWxh8NNTRIgRrMgjUOKlV31KX9v79Ef76Qx3GsE4nixRF48Q\niUemVEnHOXJ/q8VBnsOX7M3ah4+PTJM4/HVXCOep85i6uiDHgf5cNVSVT2nParbpJu2DUezmsd9g\npJTWmE+dwfrPP8PZ+AsSnQEjKJcUY6tdh/W+5djuq8Ykeb0BI8hZS0uwlpbg/vbvEbt4ifAHLfDO\nu9iP/whrVQVde/4N8dq1M76HxWwE6Wg89UEaxgjUJhOqtDgt95oJZbcZi8hOnUXfV4VyZ35h51w1\no9kQrXUDRrENIcaktWYwHjGGhJNl6Yar6UT7k8eMIuxGkO0gFL3d4w3FQkQGQxTeCFN+aQCl4XKJ\nk+tFduJRE45LdqMGrMWFy+Iky+Iix+W5KxXi0HzrWMctM0z0rwPdqPdPQfstyPGgf+dzsGx6W7vs\ns9yCMxWmVj/Wt39K1o8biHd0YsrxYK9dj23lfUZQngMr3uc6pRTWinKsFeW4t/0+kY9O0f+Tt3D/\nd89iX/cAHf/Dn5Iom34v1Zr81Yum8XdgOFDH4qjfNqN9XtT/z96dR0d2X4ed//7ee7UBqMK+oxf0\n3uyFvbFJ0SRFcZVoW7ZELWacOMrYlOMczTlRIsWSxxPPJBPbipxRHDuJbHqJk1Foi6ak2AllmU1a\nFCk2yd5JNhu9AWg00Ni3KqD2937zx6tCA2jsKKAKwP2cU6dQVQ+vfo0G6r7fdm9xAV2I7WiGy9fg\nejvcvS/frSlYswZppVQzbj7u9tVrjigEWmsSdiITWN1gGssE2GxN2NsB1g3Cw5PK1U2UrZtjiDib\nncktVecnYBVRESjHbzUQsAKUj6Sp/aCbolgafc9mjAN7CRSXErD8FFmBGbM4rTQdHkOdfx91qwdK\nStzg3FS/pCFhn8cguRJDnd29eP/mFYr/+m+xb3WjiorwHj2E/54jeHbtyPse8bVMGQb+QwfwHdxH\n/PQ5xr7zVwR/7lkST/80iV/6++hpKTrn4sn0pJMr1JPOcgP1YdT3/w715mn0Yw8WzBSGskzYuhmu\ntKL375GV3rOYryf9FaXUUdxe8x9IwC5ctrYzPdHoRADNFk93Vwe799mSdHE7RjQVzczBZva/puMT\nXy8kwAYsPwFPEQErQKkvRF1x7UTeYf+UurB+JteI9Vt+DDVzsNBj4xinzqO7elA1++DQ/rzvp9TR\nGOrdD1DtnRDwo+87ClualvVh57MUyfTSMk7dIZnE88rrlP3F90hevoryevDcfYDgJz+O967dBZnc\nYi1ThkHg+FH8hw8SfeWHqL/+Af7/9beM/7N/QuonH1/QFq5skE6vcJAGUB4P+v5jqBM/cldU79mx\n4u+5YNs2w5XrcPMWbN2U79YUpFn/erXWbcA/BlBKHUYC9qpI2knCyTDhxCjhZJjR5CiRRJhwMsJ4\napyx1BiD8f7bK4ozQ8RxOz7neS3DcrfpmO4+WJ+VSYto+Sj1lbp7Y61J+YdnCKzzBdjl0u03UW+d\nQXu9qPuOoDblP82hbr+JOv0uGAp99CBsy83ea6+llt+TTiTx/vXfUPzH38IZHkHv3E7oH/49fEcO\nYUh+5BWnPB6KP/o4gZ+4j8hf/g/0//PvSLTdIP6FX5o3UGfnpFdiNGUmqqoCvW836sIHbpayopXZ\nCbFYKlgCtVWoK61oCdIzWugWrHPcDtiPIgF7Xtk52bHUOOOpsUn3Y4QTbvAdTYzQG+0hkoxM3BIz\nFFy3lEWxp5giT5F7swJUBipospoosoom9WozW3MyAdUNxO5+2EIutK7TaYy3zkL7TXcV6uEDKE9+\ne386kUS9dQbV3QubGt2KQ97cFTFY1nB3KoX3ey+5wXk0jPf4UYqfegKrvgCTV2wARjBI6T/6+3g2\nb4Jv/SVlnd30/NavwxxTC54VXjg2oz074HKr23M9tH/13nc+9XXodz9Ap9My6jODRf9EMqu3XwEJ\n2Fn/+uT/RXl/GWnHxnbSpLVN0k5ia3vG4z2Gh6A3SNAbJOQNUl9cx87ynZR4Sgh5gxR7SijxlhD0\nllDiKcFnru4WndWkR8Ko106iYzHUPYdQBXA1rXv6UCfPgNboB+5FNdXnPDuP11SLn490HDwv/5Dg\n7/0R9uAQ3vvucYNzbU2OWyeWoujRD6N8XsL/7c+p/Wf/kt7f+b9n3VttGasfpJXHg967A/XBFfSe\nnYVTjaq22q0A1zcADXKhOd2yLltmCdjNuHunN0zA3lG+k+baLZiGhWWYWMqDx7QyvdwiijwBAlYR\nxZn5W69ZmGXlVpvu6ET9+BSqpAQefQgVyn+tWf3eJdTFy1Bbjb7vKCqwMlukTEPhLGJK2mi7Qfm/\n+Fekb3Zi3X2Asi98HquhfkXaJpYu8MCHUF4vo3/y36j9Z/8nvb/7b2Yc+s52sh1nFXvSADub3QIX\n19pg/57Vfe9ZqFAJuigAPf0SpGeQs7GFjdzDfmzLIxzYKRvyF0Pf6kG9/o67d/PY3QWxslO/34K6\neBl9eD9q944VzW1rGaBxP6QNY575yx/+mJLf+BpUVlD+L/4p3h2SGr+Q+Y8fRWtN+I//K9Zrb5J+\n+CfuOCb7X76ao90AyutFb93kLoIskCAN7py56huYI6vBxrUiEwCzBOxyrfViKmmJdUr3D6JeO4lq\nqIV7DhXEtiB96Srq/Rb04QOo3dtX/P3MzJyk7cwxdek4+P7ov+H/k/+O98ghQp/7eVkQtkb4jx8l\n/vZpjK/9HsP3HoVpIzJKKZTKQ08aoK4GrrW5uxYKZAEZJcXovoF8t6Igrfino9b6Fa31P5YALSAz\nB/3qj1GVFXDvkcII0FfbUBcuog/etSoBGm7PSaZn+5COJ6j5wlfw/+nzlHzipyn95X8kAXoNUUoR\n/LmncSIRfH/632c8xlSQjxhNTSWg3OHlQhEsgXjCLT8qpsj/J2SOKKWalVJfVko9rZT6UqYa16KP\nXcx5xOJord0edHEA7j9WGEPcI6OoMxfQe3ehVrEij5n5y7Nn+ZT2//5zJN77gNJf+SWKP/b4ul04\nuJ5ZNdUUP/4I/ue/A6k7g4+h1Kz//ytJeb1QFoLBoVV/71ll04KOjee3HQVoRYa7M4HtM8D1VSy2\n8YLW+tik938BeGIJxy7mPGIxunshMgYf+YnCKWv4/mUoKoIDqzs/Zxq3h7vvaFN3L77vvkTJx5/C\nf0jWOqxl3n17GP/+32Lc7MLZtnXKa4aRp540oEMlqEIKiNmtV/bMO2I2shXpSWutR7XWzwHblVJX\nV+I9JsskW5n4dddajwLHMqUzF3zsYs4jFs94rwVVUY6qqsh3UwA3kxid3ei9O1d92D3bk3b0nZ/S\ntf/uP6OKAgQ+8tCqtknkXnYFvtl6447XlJr5/39VeDxQSEPL2ZV0i9nysEEs+5NJKfVLs72WCdRt\ny32PBTiGWy5zsmz5zMUcu5jziEXQQyPo/kHYVUA/yuvt7hV8HvZmq1nWjqvuXmJvvk3xk4/JHPQ6\nYJQUY4RCGK3td7yW1wkMr3fGIfi8ub0nLb/tKEC56D48Ps/rZ3PwHvMpm+G5kVmen+vYxZxHLEY4\n4t7X1+a3HZNFxqGiLO/ZzSYzunvdVd379+a7KSJXFPkb156NdsjzZcJU2eC8QimH17JcfDp9Win1\nKdxgfAJ4edo89LUcvMd8RoDpY6hlmecXc+xizjOhq7ObkmBwwY3diHR3L6qvH3W9FeUpjGQuun8A\nNTwKMwxFrrRbXTHiQ/103kgRLrq9gE6hKUmnGbh0GelHr33O2BhDg0PESoOk26f+nkUHe+jrGqOt\nPJyz9+u82bWwAy9lZiHz8Ls/Ez08gu7rR/f0osYLaK58ibo6u3N2rlwE6V8F/hJ4DLdX/cuZBVdn\ngb8AtgN/lIP3mctp4PPTnqsAWhd57PAizjPhP3z99ymZVrT8sY89yuMfe2zuVm8k2fSIaQcKY80Y\neD2QLqAhP0CXl6G8HpxCWnkrlizd3QuA07w5zy2ZJhZH11QWTl86kXTvfYVxAb8YL3//BCe+/8qU\n58ZyuChv2UFaa/31zJfPZW4opY4Aj+Kuin4U+JXlvs88bTinlJoYks58fT2b5SyzIGxEa902z7Ht\nc51nNl/73X/DAVmFOyddGkJdu4GqrS6chWOpNCoyjm6sR63yh8OwEcVf4aFpSz3lxVP/DCs3NaE6\nOimvqy2IfeRi6SJvnCQWDDJy/OgdWWuKKk1qGstp3lae8/dt3rZl1te0baOCJehdO1BzHLeadNpG\n19fC3tXbBpkrn//CL/H5L0xdmvXe+ff46AM/lZPzr9Tq7rNa669rrZ8AvrIS7zGDT2f2NT+N27v/\n9KTXvgo8vcBj53pNLFVlOSoUhA8u57slt21pAsNEnXsv3y2ZYuSf/xNSV68z9r3/me+miGVI3+om\n+uqPiP38p+asiLXqevvdrU611fluyW3DI6hSSUkxk9VYMfOHq/AeaK3PA+czD1+c9tpnFnHsrK+J\npVNK4Rw5gPrhm9DbjyqADwjl96EP70OdOo/e0byqPfzbW2/uHHC0jx8h9oVfhN/7Izybm/AfO7Jq\n7RK5obUm/K1vY1ZVkvgHn5n5mFVuU5ZquwnBIKq0cNbR6KERt851vhtSgFbk8k4pVaqUelYp9Uhm\nr7EQbsnHinJ47xK6ULZabNsCZaWod86h06uXSCGbxMSc5S8w+fc+hf/4UUb/9FskLl5atXaJ3Ii9\n/iapq9cZ/ZdfmnWe1XHAXOVMctq2obsPvbVpVd93LjqRdDONVeZ+2H89WBfJTMTa4Rw/hB6NwFtn\nCiJQK6XQ9x6B8Sjq1TfcD4xVkE0Hac1WAUspev/tb+DduZ2R3/3PjP6Xb+EUUoYoMavYW6eIfOvb\nBB68H/vYoVmPc3QeRsFbOyCdhs2Nq/zGc+jpc++rK/PbjgK1XpKZiDVCVVWgH/4QuruvcAJ1eSn6\nkQdgbBz18o/cTGQrLFtYY7aeNAB+H/1/9A2iX/2nJM6/y8Bv/Btib59G5ytLlZhX9PU3Cf/p/4f/\n/nvp+e1fn/NYx9EYq9iT1qk06r1LsKUJFcx/7fYJN2+hKssLpyJXgcnFnPTjzL3FajWSmeTVya6T\n9JX3YRoWHsPCVBZe00ORVUTAE6DIKqLIU4Tf9GMa+S8qkW+qoQ794Q/BaydRJ0+j7zua92IbqrIc\n/fhDboWuv30N/ZGfWNE5u+xw93y1pDEMUj/zMUYfuJfqf/k1wn/8X4m/dYrQz38Gs0p6HoUk+upr\nRP78RQIPP0jPb/7avN1kR89zkZZrV65DOo2++66CmfvViSS6tx99z90F06ZCs16SmeTVye63uRj4\ngLSTxtZzz2t6DA8BK0DA8hOwityvPQECVoAiK4A/85rfChAw/QSsAH7Lj9/y4zN9E7e1HuxVYx36\nI/fDD99EvXYSffwwatpe81VvU7AE/cSHUX/3Y7dHfWQ/NG9ekQpUaXsBPelJdGUFff/xa1ivv4Xx\nm99g4Dd+k8C9xwh8+AE8W1Y/ram4zYnGGPvuXxN77Q0SP/9pRr/wi25i7jlordGwaj1pPTaOunQV\nvXtHYfVYO7pAa9hUQMPvBWa9JDPJq9968GscOeKuwHW0g+3YxO0446kxxlPjjKXGGEuOEUvHiKVj\nRNNRxlNjRFNRxlJjDMeHGIwNEkvHiKfjxNIxUs7cSTYsZeGz3IDtNb14Da97b3rxZe7d53z4TG/m\nuMzX1u1g77f8+E0/PsuH3/TjMTyrVhZR1deiH/8wvP42vPwjOLQPleekDyrgd3vUp99FvXMOrra5\nFxDlud0e4miNUov/kE4/eB/Df/mn+L79PYxv/xWxN07iad5C4MMP4D922C1DKFaFdhzib77N2Hf/\nJ04iQezLXyD5yZ+aN0ADZK7RVqUnrW0b9cY77gK2u3au/BsukHYc9JXrqE0NIHnqZ7UukpkUEkMZ\nGKaBx/QQ9C59uDTtpImmo0RTUWKZ+7gdJ56OE0/HiNsJ4ukYsXSchJ0gYWfvE4wmRhhLjpOwEySd\npHtvJ0mkE6R1et72+0235z7Rizf9E717v+UnkHl94nHmmNujAG6v31hAHl5VVYH+yccw3jmHPn3B\nLWd55CAqj3+0yuOBDx1F72xGvXMO9YMfovdsh327c1ZiM2XPsWhsPsVFJP7R3yPxDz6L9ebblP/Z\nXxD+L98i8sJ3Cdx/H4EP/wRWTf63uK1nyeutRP78RdI3buI/fpS+X/8iehE/c3tiTcLKXxCr0+/C\naMQdJSqQErGA24uOxnAOFs7weyFakX3SWuuzuD3pryulvrwS77HeWYZFyBsi5A3l9LzZXn4806t3\ne/YxYunoxONYOkYsFZ24SIimo4wkhuiPDWQuEuLEMxcFc/FPBPNsIJ8U1DOP3V6/D1+zl5LiIBXn\n38PqaME+tBdr8yZ8hjfT6/diGatbCENVVaA/+hG42op69xK030QfPgCbGpadCSyZdvBay/xoskzS\nD91P/0P3Y9zswvu9l+B7LxF9+VU8O7fjP34U34F9mBWytSVX0j29jP+vHxB/+zTW5ibG/uD/ZfTu\nfYs/T2ZWzGOubHjSrR2othsrMhq0HFprdMs1VEMtlOX2M269WTfJTMTCmIZJsVFMsWf587+2YxNL\nRxlPRSeG8bPBPjoR5McnXoumogwnhhhNjE4ZBUjayYkevqfGYce1cSq/8wPCQYv2rUWES92rf7eX\n78Nn+CYCt8/0uc/N+Dh7rHfiucmv+U3fvMP7yjBg9w70pkbU6Quok6fh3SL0vj2wpXHJC94SKb38\nID2Js6mR+P/+LPFnfwHPD39M2V98j8jzf0nkW9/GamzAe+AufAf24dm2Ne+L9NYanUySuPA+sTff\nJnnxEkYoRPTXvkjqJx+HJf4ss2sSVjJI67YO1DtnoXlLwaT/nNB6AyJjOA8el170PFY8SEsyk/XL\nNExKvEFKljGsn3W7hx8nYceJd3aQPn2G9PAwMX8p4d2NRIuMTC8+MXFszI65Q/zJUfrjiYkh/+wt\nPc9CPgNjIsC7PX83sE++95t+N7hv81FSU0zF1V6K/+46ZlER7N2JuX0bfl/xgoJ+ViIXPemZ+H2k\nPvoI/R99BCJjeN4+Q8X/OkHsx28R/ZsTqEAA767tePfsxrt3F2Z93aqtQVhLtOOQunqd2FunSJw5\nj47H8WzbSvTX/zmpJx526zEvQyoTpFfkdwDQ19pQpy9A8xb08UMFFQh1PIF+r8VdfyKjPPNadpDO\nZBV7df4jhZjdHT38PY3o3fdht7eTOHUKfT6CtWM7nn37MAILX52aclKZHnuc2KTA7t67i/TidpxY\nZgg/O9w/mhxhJDlKPJa9KLgd+CmFgMemqStGzfdfxjYUtxr83Grw43isKYHdb7rz8+5wv29ifv+G\nncYsMzg/1DNpZCBzMWC432Mpa3kBNFhC6rEP0/vYh8FxMD+4gvXOGawfvknkxe9B2sYIhfDu2Yl3\nzy48O7dj1lRv2KCttSbdeYvEmXPE3j6NMziEWVVJ7OefJvXkozg5TACSDdLWCvSkdVsHqn8QvWcn\n6tC+ggrQALx3CRQ4xw4WXtsKUC560l8D7snBeYSYQimF1dyMuXkz6ctXSJ47i93WjrVtG9buXRgl\n8ydk8BgePF4PkJs9z452iNuJ23P3Y8OkWi5jX75CqjPJWG2I0aYKwmUeYnYiMwUQYyQb9NPuKMCI\niuGUJ/jzttkTk5jKdFfoG+4Kfp/hxZOZo3efz8znTzrGa3jxTXztjhB4jcwxe7fi3beLxP/28xCP\nY124iHn6PMbrbxE/dRa0RhUVYW1qxNPUiLUpc6uvQ1mrux5gNei0Tbqzk+S1VlLX2khda8UJh1GB\nAPEnP0Lqo49i371vQau1F2uiJ53DIK0dB33pKup6O/qRB1D79+Ts3Lmie/rR7TfRHzqG8smK7oXI\nxV/esFLqS8BfzlfSUYilUKaJ5669WNu3kWppIfX++6SuXsXavBlrz27MitUrjGEog6LMnnYAShqh\nbj/6Q0nS16+RunwZfXkcVezH2nkQq3nrjBcTv/69S/SMxvjnT2wmbru9+YQTJ27HSNgx4hNfx0k6\nCZJOgkTmPukkiKYjjDjDJOxk5rkkSTtJSs9fH9tjeNxAXuTF+xEfvkfvJhQ/RFNbhIbWERqu9VP2\n7hmKXvkhANowSFWXYtdWQm0VRl013ro6fHX1FJWU4cUs+N63E4uRvtVD+lY36a5u0p1dpNo7IJkE\ny8LTvIXoz34M++jdpA8fWPZw9nwm5qRzNNytozHUj0+hWjvQ+/cUZoCOxdHvnEXVVqO3FVh97QKW\niy1YTwAopR5VSh3RWn9n+c0S4k7K58N799149u0jff06yQvvkv7B32LW1eLZvQezoT6PbfPiuesu\nPHfdhd3fT+ryZVLvvUfq/HmMmmo8u3Zjbmqa6JEOjiXwWabbAzZ8hHK0M8bRDgnHDezuFrw4CScx\nEezjTpyknXnOiZO03eCf8MS5vt/LB3uDJJwaknYCFY1R2TFCXecYtX0Jant7qT17g4qhFEpDEhgJ\nGAxVeBkt9xGp8DNeUUSsoohkeQl2aTGquBgrECCgvATwuDfl3vvxUIR34rF782ItYOvedNpxcEZH\nsYdGcIZHsIeHcYZGSPf2ke7qxhkezvxHKcyaaqyGemK//AvYB/dh796x4kF5ulQOF47pm12ot8+D\nZaLvP4oqL1v2OXNNaw3vnANl4Dx4b8Ff1BWSnI1haa1fyVS/+hLwh1rrcK7OLcRkyrLw7N6NtXMn\n9s2bJM+fJ/7aaxhlZXh278LcvDmvw7NmdTVmdTU6lcLu7CR1+TKJN94Ay8JsasSzYyexRIrqoD/n\n720og4BZRMAsglwE/ntuB/6EHWfMSTAUHYHODoyODsxb3Xj6+mnouInn+hi+U0N4YlN782lLES3x\nMFZsEQmaRIoNwh6DpDd7U6Q8BimPgWOAqQ28GPgcE6828GoTr6MoimsCCQd/zMYXS+ONpfHEU1hj\ncazRcZQzaerA50WVlWLVVBP96Sdwmrdgb9+Ks2XTrFWpVlMi7bbV71n6Vj6dSLo7Dm52QVODu0Cs\nsztXTcyt91vQ/YNuoiBJXLIoOf0ky6zk/h2l1NNKqWFZUCZWkjIMrC1bsLZswe7pIXnhAom330Gd\nPY+5dQueHdsxyvLXq1AeD1ZzM1ZzM05kjHR7O+lrV4m332D/2W6CWxpROz3omloo4DnfKYEfwF8P\nFXvh4J3HpoH0WATV24MaHkYND8PoCIGRYYqGh6kdGcZ/qwM7lcRJJtHJJDqVhGQKUmmU7aCVQivQ\nRuY+8zjpN4n7TWIBk4GAwViZwbhfES0uYqislJFyD8OZWzRgTswlm+oKXqMdb/wtfNcmz/G7c/sT\nc/yT5vezz2W38XmnPV7unv1Eyg3SPs/ie5TacaCtA3X+Imjc+d0tTQW7CEtfaUW3XEMfuxtVU5Xv\n5qw5uVjdvXX6XLTW+kWl1Fal1G8BvyW9arHSzLo6AnV1OJEI6Wvu3HD66lWMykqs7duwNm/Oa7Yl\nI1iC98B+vAf24wwP09byfR5IRLFOvgGmia6uxm7cjK6vB1/ue9irqiSILgky25K4pRYDtYCSzC3L\n0c7EnHx2+D6ZGcpPOIkZHt8e/o/rCJF0hMHEYGZ6IDkxTeAwd3W2yYv6fNngPmkP/0RAN32Tjrl9\n35dwML0xks44CduPd4Fb93TfgLu1KhyBLZvQh/cXdM9Ut95AX7iIPrAXtWdHvpuzJuVkdbdS6pu4\nObqPAtuAY0AZoHBzesvqb7EqjGAQ7+HDeO6+G7uri9TFiyTfOUXq3HnMLZuxtm3HrFy9hWYztrG8\nnJaKLew6WEtqezGquxuz8wbW2VOgNbq8HLtpC7quHkKSjWkuboKbAH4zAOQuo1baSU0E88Skef3b\n8/3Tn3f368ecMCPJEXe7XuYiIZG5iJiuch985VTm35HZr+83M3v0J7bwudv2SiM2ddeHCA7HoaKM\n+KFdWFUW/lgb/qR/YstfzI7jM/I/nA+gb3Siz7yL2rEV7r4r381Zs3JSBStza8VNBXoON8tYq9b6\nXA7OL8SiKcPA2rQJa9MmnPFxt3fd0kL62nWMsjKsrVswN2/GKF79yluO1sSSNgGvCcUl6B07Se/Y\nCYkEqrcHs7MD64P34b3zUFyMXduArq5GV1WBv4AqGK1jluHBMjwUW7mpu5zt8Wfn9l++fIs3rvfw\nr5/efns736R9+tmv07e6CLVfxzsQpj0ArZt9dJcNkO5rgb4732e8exyA8u7yiSDvm3bzm947Rsqp\nOwAAIABJREFUHxvTj7mdwS+b7MdUC8+upm90ok+dR23dhHPvkYIdil8LchGk/63W+is5OI8QK8Io\nLnZXhR88iH3rFumWSyTfex/OX8CsqcHcugWraRNqlRYUxZI2GijyTvvQ8/nQm7eQ3rwFbBvV34fq\n68Xo6Ua1XnWPKQli1zegq6rRlVUge03XhCk9fg/4bPCmfNxbc+cgo06nSd+4Qer9i+hwDWpzGd6n\nDmBuuV02NZukxw3ot3PuX71+jZgdp6Kx/HainkzAz2bmG02OkpiUoCdpJ0nM0NOfzmN43MBtZCvo\nZbLyTbsYqOoYpepKH/aWOuLbN+MbuOgm+MkU3snm9F/Nintr2bxBWin1CHB6tnllCdBirVBKYTU2\nYjU2opNJ0jc7SV9pIXnqNMnTZzDr690edkPDiq4OH0u4ecqL5lrZa5rounp0XT3OwUMQj6EGBlAD\nA5i3OuHqZQB0MIRTW4eurJKe9hoST9n4rKn///bAIOnWVtLXr0M6jdnYgOdD92HW1Nzx/dkkPaFp\nKXlDYXd6ZMfW7YtqTzZJT7Z4TtyOTWThm5yRb3qmvrgdZzQxQiQ1RiIdp/pKP/bNEd5u9HIt0AHv\nvzPreyrU7cA96X4iPW+2QM/kHv604ydumfK7C6m8t9Ys5JNIASNKqevACeBl2Qst1jrl9eLZvg3P\n9m04sRj2jRukrlwh8eM3UR4PZmMj5pbNmLW1OS9IMRpzg3SJfxEXAv4AumkTumkTDochGkUNZoJ2\nzy24dsU9rqQEu6YeXVGBLi+HYGhFMmaJ5YmlHAIeEycaxe7oIHX5CjocRhUF8Ozbh7WteUEZ9XLl\njiQ9i6QTSRJvvold0ovvF+/Hs2sXtrZJpBN3VtjL3OLp+IyvDSeGGIwNuYv70lNz8c/HY3jcoD49\nyE8K/gErMHHzW4FJpXj9BKwiijwBvIa3YHr5835KZPY/fx04jDvnrGZa0S3EWmUEAhh79uDZs8dd\nHd7eTvrKFdLt7W7Arq/DbGzCbKhH5SDpxUjU3UccWkyQnq6oCF20Gb1ps7sOOR5DDQ6iBvox+npu\nD4+bFrqsFKeyBl1eji4tg5ISCdz5ojUMDxNqv8qR/h5iL7aAaWA2NOA5fhyjrrZggsNC2T09JN5+\nB9JpAk8+iVnvJhUylUmRp4giT1Fu3icT9G/n2I+6aXYn8u9PDvyT5vfTUUYSwwzHh4inb6f0nSvo\nuxctRQSsAEWeokwQnxrQ3ce+zPoFE1NZE/dtI205+TfDwoa7Pwm8rbX+1TmO+TKgtda/k7OWCZEH\nRjCI98ABvAcO4IyMkO7sJH39GumTJ8Ew3DnspibMpsZFFfqYLBukg8sJ0tP5A+jGJnRjkxu0UykY\nGUGNjmAO9GF2dcCVS+6xpoUuDeFUVKFLy9zAHQotueyimEcijurrw7x1E9XXD8kEZZ1hVHUNvgce\nwGxsyMnF32rTtk3q3fdItbRg1tbie/hhjKLcBOSZ5DzoO/ZEEI+lY4ynxommo4ylxhhPjWduY4wl\nx4ilYwzFB4kkI5kLgdjE1ICj79yuN3I9d8UfF/Ip8bjW+lfmOkBr/XUApdRvA78p+6LFemCUleEt\nK8O7fz/O+Dh2Z6ebjvTMGTh92t2D3dSI2dSEsYitUqOxFB5T3TEnmVMeD1RXo6ur3ZXj4OapHh1B\njY6iRoYxentQ16+6vTul3JXmoRB2RSUEQ+hQCIql171o0ShqaBA1OOCOaoTdj0NdVkZ6+050bS0n\nrX72NpVhNW/Na1OXyhkZIXHyLZxwGN+992Ht3bPmRgDcUrsllHiXN61gaxvbsUk7adJOipST5szZ\nMzzJkzlpZ64zjn0lkxZUetRiXTGKizF278azezc6kcDuukWq9Tqpix+QvPAuRjCI2dCA2VCPUV09\n5zz2SDRFiW+ZZSiXwuuF6hp0dQ0a3B53Og3hMCo8igqPYgwPYV25DMnMUKBpoktKcErL0cXFUFSM\nLi6B4mIIBDZ2ANcaouOoSATCYcyBPtTIMMRi7uslQeyaOvSefeiamilJasbTfZT41t7IhU6n3Yp0\nFy9ilJQQ+OmPY27wmtCmMjFNE695ezSkKpC7zGoLCdKLzfzwolLqk7K4TKxXyufD2taMta0ZnU5j\n9/Zi37xJur2d1OXLKMvCqKnBrK/HrK/HCE69Uh+JppY3H51LlgUVFe5CM7idZyseg0jEDd6RMYzR\nYVR/H8Sit7/XNCEQQAcCOMEydFEAAkXoQBEUBdygVMDpThcskXB7x7Go+zOJhN2fR2QMbHcRIJaF\nLisnvaUZKirRFZVzbo+LJm2KfWvrZ5Nuv0Hy3XfRsRje/Qfw3H1wXZYwLTQLWt2tlAotdAhba92m\nlHp6me0SYk1QljWxrct3333YQ8PYt25hd9wgee4cnDnj9rLr6zHq6zBrahiOpXI7H70S/AF3nnty\nrxvAtt2ANT4G4+OoaBQjMooxNIDqiroBbTKPB3x+tM+HUxKc+BqvF7xedOYerxc83tWbF9fanbdP\nJCCZQCUSt7+OxzHCo6ho1O0VZwMxuG0OBnEqqtFbt6GDIXcF/SLXJ8SSNiVrJEjbAwMkz57DGRzE\nbGrC99GPLmp6RyzPQn5L/gB4DvjsIs5bubTmCLG2mRXl7vDf/n1uFayeHuzubtLt7egrV8AwCJ4f\npay2FkYqobRsbQ0ZmyYEg+iguz93SgAHN4jHYxCNoeIxN/DF46h4DGN8DDU44A6lJ2dJnmEYbu/b\ntNCW5b5f5t7x+tzXs8cpA4zMz07jBl7tZO7dm5GIu21Kp1HptBtw05mbnpZdXCk3CPv9EAhgNzRm\nRgWK0EVFECjKSfKYZNohaevCGU2ZhTM2RurCu6Q7OjAqygl87GOYdXX5btaGs9AtWF9TSv2i1vqP\nF3jebctslxBrnvJ4JlKT+o4fd+sdd3cTPvUDDg/cwPzhMHg9bvawzG3N5+o23VSnFJdMKbBxRzDX\n2g3UmQpYKplwe7aplBtAbdsNqpnAasRjGImY+33OncEYpSbdjMy92x6nuAQsz+2g7/G49z6f25P3\n+W/35lfhgimbzKY0kL+CL3NxhodJtVwm3dGB8vnwPfRhrG3Na25h2Hqx0Eu5zwCnlVLbtda/NteB\nSqlSkFStQkxnlJaiQiF+XNNB4z2NpBp9qL5ezJ5bGO+/6wYfnzeTPSyTq3u9JiNRyu2VZnqmM1XM\nmvzc3DWp1paxhA1AaaCwetJ2dw+plhbsnh5UcTG+48exdu6Ueec8W9BPX2vdqpR6DDihlPo08NvA\nC9PnqZVSW4EXgFn3VAuxkY3F06RsTWlJAF1Via6qxrlrvzscO+QmI7kjaFdUoquq3Fzda214XNwh\n25MOFUBPWjsO9o0OUpcv4wwPY1SU43/4YcwtW1DG+kuxuRYt+BJJa31WKbUNNwg/B/yhUipb+WoI\nd4j7MdyCG6+uRGOFWOsGIu7CqvKiaR/QloWuqUXX1N4O2sNDqP4+zN5ujA8ugu2Ax3KDdkUFuqIK\nysvAyv+HvVi4bE86lMeetBOJkG5rI93Wjo5GMevrZc65QC3qt0RrPQI8nulV/zLwKG6ZyhHgNPCE\n1vqVnLdSiHViYMwN0qVF82SYsix3ZXV1jRu0bRs1PIwa7HeD9rWrkHIziOlQCMrdXN26omL9DpGv\nE+OJNApWfXW3TiRJd3SQbmvDGRxEebxYO3fg2bULo3xj73UuZEv6LdFan8AttiGEWIT+cByAsuk9\n6fmYpjvkXVWFs3uvu1gqEkYNDaGGBjH7ezE62t2JXMtCl5ehy8qhrBxdXgErmK5RLE4knqbYZ2Ia\nK38hpdNpd3fBjQ7sri7QGrO+Hv8jj2A2NeW8eIzIvTmD9GL2R88nl+cSYq3qGY1T4rfweZb54agU\nhErRoVL01uZM9rCU29seGsIc6MXo6oSrmUIbPp/b0y4rh9IydFnZovf2ityIxNMruv1KJ1Ok22+Q\n7uzE6e5Gp9MY5eX4jh/H3Lp1yTnnRX7M95vyWaXU9eXOMWdqUm8D/mg55xFiresZiVFVsvy9tjOy\nPLeHyHfvcZ+Lx9zAPTzk9rZbr0HSLfCB14MuLYVQGbq01B02LwmujyxhBSwcT+d00ZjWGmd4GKe3\nl8TpM9jDwySqqzEqKvAcPoK1ZTNGMDj/iURBmvOvUWv9nFLqWaXU48AfLLY8ZWa191eAa1IhSwi4\n3jdGVXCFgvRM/AF0fQBd34DDfve5aBQ1mqmQNdiP6ulGXb/mvqaYKLShQ6Vu4pJQqZurW1b75sRo\nLE19mX/+A2ehtUaPjGD39WH39uH09aNTSXerlGXhPXyYovvvxyguzmGrRb4sJJnJc0qpZtzV3Bp4\nGXc+unWGLVgh4BhwBHgCd4bsH2utc1dcU4g1bDCS4O7NZfltRCaDlhu4M9Ipt1pTeBQ1Ooo5PIDR\n1no7zaehoLgYHQyhS4IQCqFLSqTnvQTheJpddQuvvKSTSZyhIZzBIezBQZz+fnQyCaaJWVWF5+6D\nmLW1GJWV+Nvcj1oJ0OvHQvdJtwFPZBKVfAb4t8C2zJasyTkHsqu8XwA+rbXOXVFNIdY4rTWDY4nV\n7UkvlOVxi0JUVE4rtBFHRcKZwhKjmMNDGB03IB6//b1+vxuwi4vRRcXufXGJu1gtB2k015twLEXF\nLKv7nWgUPRrGCY/iDI/gDA7iZEpdKo8Xo7Icz/79blCep9qaWB8WuwVrFHeP9HPZ5zKBGwnIQsxt\naCxJytaFGaRn4/e7uaynF9pIpdzgPTaGGougxsYwwsMY3bduz3kDWJlKWUXFmQpZAbcn7w9AwO8W\n8thAPXFHa2LjMartcdI3b6LDEZxwGCccRkci6FTmZ2eaGGVlmJs24amswqyqRIVCkppzA1r2X8dM\nwVkptXWx89dCrHfdo26d4cqVWji2mjxTe94Adva1ZNKtszw25pZ3jEYxI6MwPIRxKzo1iAN4rMzF\nQMC9z+bSzqQN1ZNza3sKOHGL47jFQxJuJS0S2a9jMD4O0Sip0QgPXblBU1E9ie4SlMeLKg1h1tRg\n7NiJKivFCIVQJSUSkAWQgyA9i18FfmWFzi3EmtQ55NZirgktfdHQmpAtQ1lWPhHAp+TeTqfcYfRY\nLHPvloRU8TjGeARjZDhTNjJ157lV5vyWB+31uOUtswUzLAs9UTXLXUSFYWQqZqmJe52toAVMzNZl\nKmKpbMEO23aDbvbecSYqaKlkyv03ZAqCqFTKLRSSmKGyl8fKjB4UY9fVM1jZxLs3ivnMUx+ieHcj\nyr/OfxfEsi0pSGcWiH0Nd5EYwPSVMNuQIC3EFDcHo5QGPAS8G3we0fJAicddgMbURS325OOyPdN4\nApWI366YlQ2KySQqlcKIjU+UozRse+Jr0mlywswEesOtoKU9FlgenEAxOuhxLxI8HndawOefuMfn\nu6M+dm/nCL2hGJWb6yVAiwVZak86u9/52zO8poBnl3jeJcusQP8U0Ao0A8/NNk8+17FKqW8Cn8f9\n7DgLPKu1Pr/y/wKx3t0cjFJbKh/MC2YY7py1P3BH2cvJbOaQ7QVrnbl33MIlzqS+fbbE5cT7Kjcg\nm+btXniOjIy7owPrYspDrIqlBulTWuuvz/aiUqpiieddjhe01scy71+Ku8L8iSUcew0oBZRkSBO5\ndLUnTFOFpOdcVdnh7gIxEk1S7DOXn3FObBhL/e0dmetFrfVXlnjeJVFKHWbSBXamV3wsk0xlsccq\nrXVEArTIJa01vaNx6UlvcCPRFGXzFVcRYpKlBulWpdSh2V5USv3nJZ53qY7hlsucLFs+c7HHViql\nPqmUelQp9duZoXEhliUcSxFN2tRJkN7QhseTlBVLkBYLt9QqWK8opb6klHoON3nJdJ9hdReOzZTC\naWSW5+c79pvZ7WNKqSHcofBjM3zPhC9+8YuUlpZOee6ZZ57hmWeembvVYsO4Oeiu7Jae9MY2Ek2y\nqy6U72aIHHr++ed5/vnnpzw3Opq7tCFLXd39TdxAfBqonPZyGTMHwqW8z7PAdu5cK6Iyz72cKf4x\nAkyfBy9j5mH5OY+dtr+7FTgyXwWvb3zjGxw5cmTuf4zY0Nr6x1BA7XrffiXmNDSepDoki8bWk5k6\nZGfPnuXo0aM5Of+S90lrrWddHJYJ4sumtX5u/qMA92Lh89Oeq8ANsgs+NjNf/Ur236a1Hs3kKxdi\nWVr7xqgO+WTB0AbmOJqR8eT63ycvcmqpc9Ivz/P6ry7xvEuitT7HpN67UqoMuD5p2Ppwdm55nmNb\ngd+c9NqngBOyiEws1/udIzSUy8rujWwkmsTRUFcm9ZzFwi25Jz3PEPCzwGqXpvy0UupLQBvuHPKn\nJ732VeCdSW2a8dhMz/lc5rVR3MVkk88jxJJ0j8Q5vi0fOxNFoRgaczOSSU9aLMZSg/R14KuZ3LKn\nZnj9l1nlIJ1JOJJNOvLitNc+s4hjXwFeWaFmig0onrIZjCSolx7UhjY47gZpWTwoFmOpQfps5n4E\nNyBPVsadC72E2LBuDkbRyDDnRjc0lsBnGQT9G6fql1i+JQfpbMaumSilZkoXKsSG1DE4DkBDmfSg\nNrKhsSQVJV6pbiUWZakLx+ZbGPZbSzyvEOtOa98YQb9Fib+AyyyKFTc4lpCc3WLRlhSkM/O2c71+\nbmnNEWL9uXBjmE2Ss3vD648k2FkXzHczxBozb5BeSorPPKQFFaJg3RyKSmENwYAsHhRLsJA56WNK\nqbtxs3wt1JxpNIXYKOJJm77ROE2HJEhvZOOJNNGkTX25BGmxOAsJ0ke5vZp7IbIpO4XY8Nr6x9Ag\nPekNbnAsASA9abFoCwnSrcDXFnFOBXx5ac0RYn253jcGQKP0oDa0/rAbpBskSItFWkiQPrGIHNoA\nKKWk2oQQwLXeCNUhH37J2b2hDUQSeC2DcilTKRZpIau7F9OLXs73CLHuXLgxzGYZ6t7w+iMJqoI+\n2SMtFm3eIK21blvsSZfyPUKsN1prbgyM01xdku+miDzrHY1LmVKxJEtNZiKEmEfnUJRo0mZrdXG+\nmyLyrC8cZ3dDKN/NEGuQBGkhVkhLt1skbmuVBOmNzHE0/ZEETVKqVCyBBGkhVkjLrTCVJV6CAUkH\nupENjiewHU1jhazsFosnQVqIFXKubUh60YK+UXf7leyVF0shQVqIFaC1pl0WjQnc+WhDSSITsTQS\npIVYAR2D7qKx5hrpSW90vaNxKoM+LFM+bsXiyW+NECvgvZsjKGB7jfSkN7ru0Tj1pdKLFksjQVqI\nFfBuxzCNFUUEvAtJ6ifWs97RGHsbZfuVWBoJ0kKsgDNtQ+yslV70Rpe2HfrCCTZXyrSHWBoJ0kLk\nWCSWoms4xo66YL6bIvJsIOJuv9osq/zFEkmQFiLH3u8cAWCH9KQ3vJ7ROACbK2X7lVgaCdJC5Nh7\nN0cI+i3J1SzoHonhtQyqg/K7IJZGgrQQOfbmlX521gWl4pGgeyROXakfw5DfBbE0EqSFyKF4yuZa\n7xh7pZiCwC2ysq+pLN/NEGuYBGkhcuhi5whpR7NHgvSG52hN13BU9sqLZZEgLUQOnWkbosRvSZ5m\nwUAkQTzlyCp/sSwSpIXIoTcu97OnPoQh89Eb3s3BKAA7aiVIi6WTIC1EjsSTNtf7xmSoWwDufHSJ\n36KyxJvvpog1TIK0EDly4eYwtqNl0ZgA4OZQlE0VRbLKXyyLBGkhcuRs2zBBv0VjuRRTEG4ltIOb\nZWW3WB4J0kLkyBuX+7irsVR6ToJY0qZ3NM6uehlVEcsjQVqIHBiPp2nrl/3RwnVzyF00tktWdotl\nkiAtRA6c7xjG0XBXY2m+myIKwI2BcSxD0Vwte6TF8kiQFiIHTrcOUlHipSbky3dTRAHoGBinsaII\ny5SPWLE88hskRA78+Eo/extCMh8tALghi8ZEjkiQFmKZRqNJbg5GuatBhroFpG2HrqGozEeLnJAg\nLcQynW4bQgN7G2XRmICu4RhpR8vKbpETEqSFWKY3r/TTVBGgskTmowW0D4yjlKzsFrkhQVqIZXAc\nzRuX+zm4SeYfhautf4zG8gABr5Xvpoh1QIK0EMtwpSfCaCwli4TEhPb+ce7eXJ7vZoh1Yt0EaaVU\ns1Lqy0qpp5VSX1JKzbuKRyn1dC7OIzauN6/24/eY7JRKRwJI2Q4dg1EpsiJyZj2Nx7ygtT4GkAms\nLwBPzHRgJjhXAH+glCrTWoeXch4hXr3Yw/6mUtkPKwC38pVbZEWu7UVurItPFqXUYUBnH2utR4Fj\nSqmtMx2vtX5Ra/3c5O9ZynnExjYaTXK9b0yGusWE9v5xDAU7ZNGYyJF1EaSBY8DQtOeGgG3zfN/0\nzBNLPY/YgE5eG0BrOCCLxkRGW/84jeVF+D1mvpsi1on1EqRn+pQcmeX51TiP2AD+59kutteUUFHs\nzXdTRIFo7RvjSHNFvpsh1pGCnpNWSj0LbGfasDRuD1gDL2utX8UNpNP/Msoyzy/Gks7zxS9+kdLS\nqXNQzzzzDM8888wi316sFfGkzbs3R/jZo035boooEPGUTedwlL//QHO+myJW0fPPP8/zzz8/5bnR\n0dGcnb+gg3Rm3nghTgOfn/ZcBdA631vk4jzf+MY3OHLkyHxtFOvIW9cHSKYdjjbLVhvhau8fR0sl\ntA1npg7Z2bNnOXr0aE7Ovy6Gu7XW55g0JK2UKgOua63bM48PK6VmurydMic933mEyPqrM500VRRR\nVxrId1NEgWjtH8NrGTRXF+e7KWIdKeie9CJ9Win1JaANdwHYpye99lXgHeB3AJRSjwJHcHvSX1VK\nZYfN5zuPEKTSDudvDPP4/rp8N0UUkNa+MbZWFct2PJFT6yZIa63PA+czD1+c9tpnpj1+BXgF+Ppi\nziMEwJn2IaJJm2PbZIGQuK21b4yPHmzIdzPEOiOXfEIs0qsXe6gO+dhUUZTvpogCMRJNMjiWZK/M\nR4sckyAtxCKkbYdXLvZw3/ZKlJq+zV5sVNd6xwA4sEmCtMgtCdJCLMLptiHGEzb3bKvMd1NEAbne\nG6G82EutLCQUOSZBWohFePViDzUhH5srZahb3Hatd4wdtSX5boZYhyRIC7FAadvh7z7o5VhzhQx1\niwlp26F9YJwP7azOd1PEOiRBWogFutAxTCSelqFuMcXNoSjJtCM53MWKkCAtxAK9erGXyhKvJKsQ\nU1zpjuAxFbvrpYa0yD0J0kIsgO1oTlzs4agMdYtpWrrDbKspwWvJx6nIPfmtEmIBTrcNMhpN8aEd\nVfluiiggjtZc7g7zwK6afDdFrFMSpIVYgG+/dYPaUr8MdYspOoeijCdsKU8pVowEaSHmEU/anGkb\n4v6dVTLULaZouRXGYyr2N0kSE7EyJEgLMY/XL/cRTzky1C3ucOlWmO21QXweM99NEeuUBGkh5vHi\nOx1srymhttSf76aIAnJ7Plr2R4uVI0FaiDmMjCd5r3OUD+2UXrSY6uagzEeLlSdBWog5vPx+N1pr\n7t0uCUzEVJcy89H7pPKVWEESpIWYhdaaF97u4NCWckIBT76bIwrMpa5RdtbJfLRYWRKkhZhFy60w\nHYNRHt4re2DFVLajaemO8NCe2nw3RaxzEqSFmMVfne2kvNjLgSbJySymau8fI56yObZN5qPFypIg\nLcQMYsk037/QzUO7qzEM2RstpvqgK4zfY7JH8nWLFSZBWogZvHKxl0TK5sE9sr1G3OmDW6PsqQ9i\nmfIRKlaW/IYJMYNvv3WDu5pKqQ7K3mgxVTLtcLUnwsN3yXy0WHkSpIWYpmNwnCs9ER7cLb1ocadr\nvRFStuao7I8Wq0CCtBDT/I/TnRT7LI5ulQ9hcaf3O0cJBTxsrwnmuyliA5AgLcQkybTDX53t5IFd\nVVIfWMzoYuco+xpDsqBQrAr5FBJiktcu9RKJp3l4r8w3ijtFYiluDIzz2IH6fDdFbBASpIWY5Ftv\ntrO7PkhDeSDfTREF6IOuUTRwzzZJEytWhwRpITJuDIzTcivMR6QXLWbxfucoDeUBakKy6l+sDgnS\nQmR87/RNSvyWZJESM9Ja837XKA/tljSxYvVIkBYCiKds/vpcFw/sqsYjCSrEDLpH4gyNJTkuFdHE\nKpJPIyGA75+/xXg8zSOSoELM4t2bI3hMxeGt5fluithAJEiLDc9xNH/2eitHmyuoLZW5RjGzCx3D\n7GkIEfBa+W6K2EAkSIsN78dX++kZjfPRg7KtRsxsPJHmcneEJw805LspYoORIC02vD/+u2tsrylh\nZ51kkBIze7djBNvRUnBFrDoJ0mJDu9Q1Skt3hI/dLb1oMbsz7UM0VxdTWyr758XqkiAtNrRvvnKV\n6qBP8nSLWSXTDu/dHOHJgzLULVafBGmxYXWPxDjVOsgTB+okD7OY1Qddo8RTDh/eI/ujxeqTIC02\nrD8/2U7Aa/KQfPiKOZxuG6Ku1M/W6uJ8N0VsQBKkxYY0Gk3yvdOdPLavDr/HzHdzRIFK2w5n24d4\n8mADSsloi1h9EqTFhvSdUzdxtOax/XX5boooYC3dYcYTNo/skyQ3Ij8kSIsNJ5Gyef7kDR7YXU0o\n4Ml3c0QBO9U6RHXIxy7ZnifyRIK02HC+f+EWkViKj8lqXTEHx9GcaRviiQP1MtQt8kaCtNhQbEfz\npz9q5dg2SQEq5na5J0xE8rmLPFs3QVop1ayU+rJS6mml1JeUUqUL+J6nZ3jum0opRyllK6VOKaUO\nrUyLRT681tJL72icj90tvWgxt1OtQ1SUeLmrcd6PEiFWzHrKFP+C1voYQCZAvwA8MdOBmeBcAfyB\nUqpMax2e9PI1oBRQ054Xa5zjaL554ip3NYTYXlOS7+aIAuZod6j7yYMy1C3ya130pJVShwGdfay1\nHgWOKaW2znS81vpFrfVzk79n8um01hEJ0OvPay19dAxG+dljTfluiihwrX1jjERTfOQuWf0v8mtd\nBGngGDA07bkhYNs83zfTJXKlUuqTSqlHlVK/rZRqzkkLRV45juY/vXyFuxpD7K4P5btmWXwOAAAX\nqklEQVQ5osCdaRsiFPBwYFNZvpsiNrj1Mtw901/SyCzPz+ebWut2AKXUEO6w+bGlN00Ugtda+rg5\nFOX/+Phd+W6KKHBaa063DfHIXbWYki5W5FlBB2ml1LPAdu4cllaZ517WWr+KG5CnV0goyzy/KNkA\nndEKHFFKheYa/v7iF79IaenUxSXPPPMMzzzzzGLfXqwAx9H8x5evsK+xlF3SixbzuDkUpS+c4GFZ\n1S0W4Pnnn+f555+f8tzo6GjOzl/QQTozb7wQp4HPT3uuAjfIzvkWkx9k5rZf0VpXZN5/VCk107z1\nFN/4xjc4cuTIApsqVtsPW3rpHIryD39mX76bItaAM21DFHlNqYwmFmSmDtnZs2c5evRoTs6/Luak\ntdbnmDS0rZQqA65PGrY+PMvc8vSxrFbgNyed51PACVlEtna5c9FX2ddUKlmjxIKcbhvm7s3leKx1\n8fEo1riC7kkv0qeVUl8C2nDnkD896bWvAu8AvwOglHoUOILbk/6qUuplrfWrmZ7zucx5RnEXnk0+\nj1hjXr/c5/aiH5RetJhfXzhO51CUf/LYznw3RQhgHQVprfV54Hzm4YvTXvvMtMevAK8AX5/hPNnX\nxBqnteY/nbjK3oaQ9KLFgpy/MYxlKO7dUZXvpggBrJPhbiFmcvLqADcGxvn4kcZ8N0WsEafbhtjX\nVEqxb930X8QaJ0FarEtauyu6d9SWsLdBVnSL+Y1Ek1zpjvDTRyTZjSgcEqTFunSqdYjrfWP8zJEm\nSesoFuRM2xBKwYO7q/PdFCEmSJAW69J/evkyzdXFHNgkxRHEwpxuHWJvYymlRd58N0WICRKkxbpz\n/sYwLd0RPn6kUXrRYkEisRQt3WF+6pCsXxCFRYK0WHf+/d+0sLmyiENbyvPdFLFGnLsxjNbw0N6a\nfDdFiCkkSIt15XTrIC23wnziWBOG9KLFAp1pH2ZnXZDKEl++myLEFBKkxbqhtebf/00LzdXFHJZe\ntFigeMrmYucITx5syHdThLiDBGmxbpy8NsC13jE+eWyTzEWLBXvv5ggpW/OwDHWLAiRBWqwLWmt+\n7weX2VkXlBXdYlHOtA2zqaKIxoqifDdFiDtIWh2xLrx+uZ+2/nG+8lN7pRctFixtO1zoGOaZD23N\nd1NWVTKZ5MSJE3R1dQHwuc99jr6+Pl566aWJxx6PJ2fv92d/9mckk0nuvfdeDh48mLPzbgTSkxZr\nnuNofv9vL7OnIcTeRulFi4Vr6Q4TTdp8eIMNdV+6dGkiQFdXr3zyFqWUXDwvkfSkxZp34mIPHYNR\nfu3jd+W7KWKNefv6INUhHzs3WAGWRCIBQFVVFT/7sz+74u/3C7/wCyv+HuuVBGmxpiVSNv/hb1o4\nsrWc3fWSo1ssXDLtcKp1iJ+7b8uG6uV997vfZWBgAICBgQGee+45fu7nfm7GY5PJJG+//Tatra0k\nk0m8Xi+NjY3ce++9BINB3njjDS5dukQoFOKzn/3sxPe99NJLdHV10djYyFNPPTXrcHdXVxcXLlyg\nv78fgFAoxPHjx2lslKQyWTLcLda0F97uYDia4rP3bs53U8Qac6FjmFjS5smD9fluyqrau3cvVVVu\nKU6fz8ehQ4fw+/0zHvud73yHlpYWkskkPp+PZDJJW1vbxNz13r17AQiHw0QikYnvyw6lb9++fdZ2\nXLp0aSKYJ5NJkskkAwMDvPTSS7S0tOTk37oeSE9arFkj40n+5LXrfGRvDXVlgXw3R6wxJ68OsLWq\nmK3VJfluyqras2cP4XCYgYEBgsEg99xzz4zHDQ4OEolEUErxiU98gsrKSrq6unjppZcmgnJlZSVe\nr5dkMklXVxd79uyZCNAA27Ztm/HcyWSSN954Y+KYRx99FGCiZ/7666+zffv2BS1ee+ONNwiHwwA8\n/vjjM37PhQsXCIVCNDc3A0xcZCQSCR577DGCwTunO6Z/T75IT1qsWX/y2nU08DNHpbSgWJzxRJoL\nHSPyuzOHyspKPve5z03MJ1+6dIl33nnnjuOyvenW1lbgdi+6qalp1iCbHd4GeOihhya+Pn78+IzH\nzKarq4tEIsFTTz3F3r17efvtt+845rvf/e6Udre0tNDU1MRTTz3FoUOHuHDhwrzfk0/SkxZrUsfg\nOC+eusnT9zQRCuRuq4jYGE61DmJrzeMHNtZQ92KdO3duShALhe5c97F9+3YuXLgwEZw7OzuB28F7\nJpOHxicHcq/XO9EzD4fDNDTczgL30ksv8eCDDxIMBmlpaSEcDuPz+SZWp1dVVXH+/Pk73usTn/jE\nlIDb2Ng4ZXjf672z6tn078kn6UmLNelrf3WRsiIPj++XD1mxeCevDnBXQylVQcnVPZu2traJAP34\n44/z7LPP8tRTT91xXGVl5cRwcUtLC4ODgwBs3bp11nNPHl5OpVITX2fnpuHOC4Jt27Zx6dIlwO21\n7927l0QiMXGuYDA4JfjP9d4ej4cTJ05w4sSJOS8mCoEEabHmvHqxhzPtw3zm3s14LfkVFovTMThO\nS3eET8liwzll53l9Pt9EwM0Gyemyc8/Z4ebZ5qKzJu/N/tGPfjTx9eTh6un7t7Pz3dkgHgwG8fl8\nE4E5EonMOLc8m8ce+//bu/fYKq47D+DfHxA7hYKNMUlatqlfDabdKviB6LbBeWDoAyUrQSCtGu02\nbQxZadttoxBgtdp/qqbOQ1rtXwtx9qWoUh0DaZuXEky2SrSkNbaBpFmnwr6QsJCkGDC4LnF4/PrH\nnLkeLnMfM3funbnj70dC4s6dOffMvdfzu+fMOb/TjrVr16K3tzfnY8LA7m4qKaPjk/jJL99Ga20V\nltcvCLs6VIJ+NXgc1XPLseovbwi7KpFmjwCfnJxEV1dXshva5vz/kiVLcOjQoeS2bEG6rKwMK1as\nwOuvv45EIpG8n21ra2tzvZ89d+5c9PX1JVu/1dXVyZHgo6OjOSVm6evrw7x589DY2Oja1R01bIZQ\nyVBVbOs+gFkzBd9ZUTut5rZSMI6f+RP6E6dx/231mDWTlz+3vyF7mz0fury8HOXl5bjuuuuwdu3a\n5H4jIyPJ/8+dOzcZ1EUkpxHRjY2NWLNmDaqrq5OvsXDhQqxZswaLFy92PWbJkiVIJBLJ8hctWoTy\n8nK8+OKLGBoawvLlyzE+Po7u7u60r7t06VIkEgk8++yz6O3tRVtbW9ZjwiSqGnYdSpaINAMYGBgY\nQHNzc9jVib1dfe/h8ReG8ODXF+PmG7kUJXm3fe9h/P6DcfziR7fG8lbJ8PAwAKChoSHkmkxvg4OD\naGlpAYAWVR3Mp6z4fUsplt4bncC/vvx73PH56xmgyZcPxs7jNyOn8N1b62MZoCme+E2lyLt46TK2\n/PwA5s8pwze/xME+5M9zB46j4hPX4M4mppyk0sEgTZH3n68l8O7oBDbdXo/ya2aGXR0qQR+e/Qj7\nDo/ivjZ+h6i0MEhTpP166EP8x69H8Nctf4H666fXSkUUnF37j2HuJ65hhjEqOQzSFFn/d/ws/nnn\nm2itq8JdzeyiJH8OfzCO346cwvdXL8a1ZWxFU2lhkKZIOnHmPH709AA+UzUbG29vwAxOtyIfLqvi\nZ/uO4rPVs/GNmz+d/QCiiGGQpsgZP38Bf//f+1E2awZ++LXFHIlLvr1xeBRHTk5g651fwIwZ/KFH\npYdXP4qUCxcv44dPD2Bs4mM8+PVGLp5Bvk1euISevmNora1CU01V2NUh8oVBmiJDVbHl5wfwzvvn\n8IPVN+HT87lGNPn34qETGD9/Advu+kLYVSHyjUGaImPHq8PYd3gUHbfXY8miirCrQyXs5PhHeOHg\nCXzryzVYVDU77OoQ+cYgTZHw0qET+K/XEtiw/Eb8VUN12NWhErd7//9jTvks3NeWeaEHoqhjkKbQ\nDR0/i0d++TvcclM1vnEz14em/Lw/dh5vDI/i/tsaMLucC/1RaWOQplCdGp/Egz8bxGcWzMbfrqjj\nylaUt18NHkfl7DLcybn1FAMM0hSajy9exj883Q9VxQ9Wc6oV5c9uRX/vVqb/pHjgVZFCoar4p56D\nOHpyAt9ffRPmz4n+4usUfWxFU9wwSFModva9h9feOYn72urQwJzcFAC2oimOGKSp6PYnTuFfXnoH\nq794A25ZvDDs6lBMsBVNccQgTUWlqvjxs2+h8VPz8M0vfTbs6lBMnP/4Et4YHsXf3FLLVjTFCoM0\nFdWbx8bwh3OTuKtlEWYGlEv5+V/sDKQcyk+Yn8O7oxNQBZbVLwitDlHx3HPPhV0FChCDNBXVSwdP\nYMEny7D4U/MCK/MFBulICPNzODo6gbJZM1BTPSe0OkTF888/H3YVKECxmekvIrUA7gaQAFALoEtV\nz6bZdyWAZvNwGYAtqnrEaznkzeSFS3jlrfdxx+ev59KTFKgjJ/+IGxfMxqyZbHdQvMQmSAPoUdVW\nABCRCgA9AFan7mSea1bVx83jdQD2AGjwUg5597+HT+JPH1/CV27iYDEK1pGTE7i18bqwq0EUuFj8\n7BSRJgBqPzYt31YRqXHZvRVAp+NxL4A6EanxWA551POb91C7cA5Xt6JATUxexIdnP+KiLBRLcWlJ\ntwI4nbLtNIA6AEedG1V1r4i0ODYtszbrURFZlWs5xrUAMDQ05Lvi08X4+QvY99t+fPWLN+Dtty4F\nW/a5c3j7rYOBlknehfU5JP7wR0y8fwQXT87B4OAHRX/9KDl27BjGx8cxODgYdlWmNUdMuDbfsuIS\npCtdto2l2Q5VdV5JHgaw0U85AGoA4N57782pkgS8WaBy133ttgKVTF6E+Tnc1RXaS0dOS0tL9p2o\nGGoA7MungEgHaRHpAFAPRxe0/ZTZtkdVX4UVSKtS9qk027OV/4yq/rvZ5LWclwF8G1Yr+6NMr0VE\nRNPGtbAC9Mv5FiSqqfGv9Jh7yU+q6jLHttOwBogdTXPMSgAVqro7n3KIiIgKJRYDx1T1ABxd0iJS\nCWDEDqwi0mSmVtnPN5vjdpvHHSIyL1s5RERExRSLljQAiMhSAO0AjsAaSLbDEaSfAdCnqk+YYD2C\nqS50AXBGVRdkK4eIiKiYYhOkiYiI4iYW3d1ERETFJCK1IrJZRNaJyEMm+VW2Y9Z5fh22pInILz9p\ndEVknaruKkb94shjCmSmOS4QEelPzU6pqq7ZKU1wrgKwA0Clqp7L+XUYpL3xkSO8EtaHswrAT83g\nNMpTULnaKT/FulDRFI/vec77Uu78zgQSkUsA5nv57kd6nnREecntvQdAk6oeEmtBiR5M5Qin/ASV\nq518ckujKyKtIlLjdqGyW88isr14tYwXL++518+HPMk5y2UKzysL8Z60Bz5ye9ep6iHH4zOFq930\nEVSu9sLVcNrIdKHKhEug+eflPff7+VB2XrNT+sYg7Y2nL33Kr9V6AOsLU61pJ+fPQVX3AnDN1V6w\n2k0fRbtQUZKX95yfT+H4ynLpB7u7vfH8pTf3TjcBaAIwH5m7Qig3QeVqJxeFTsdLefHynvPzKZx+\nXH0dqYI1RiYTz4PAGKRR2IuSGaC0VUQ2A9hr7gdxwIyLEHK1kwtVzXWpiqJdqCjJy3vu9/OhLFT1\ngMlICcA9yyWAMZcBqp5v9TBIozAXJbsFrapbzaadAB6F1VX7qs+qxlohg4MZ4X3Kmaud8lPMCxVZ\nvLzn2falvK0XkYcwlZ3SeTtzG4A+AE8AV8wwUQDbRMRucGTFKVgeichhVf2c+X8lrNbdMvM4+Qdi\nPpROx3PtALoB1LIlnb9cPwfzuBnWtIe95nEHgG5+DvnLNR2veWxfqDoBPIapnhHywON7zjTHJY5B\n2iOPfyD324eZYx5JGe1NPgWVq52IKMoYpImIiCKKU7CIiIgiikGaiIgoohikiYiIIopBmohCZZb8\nqw2gnCazdOBa82+dyU9gv0aniFx2DOhMPX6Hef4hEZmXb32IgsAgTVQiRKQjpsFjC4JJcNKlqo+r\n6m7HfPg6wEoqZHIWPGle7wpmIZYzAAZU9QlPSwlan0vWtYSJ/GAyE6I8mItzF6zpYJWwVtiyk6oI\nrFSwnfkuUWrmdo9kCh5mHvImWAHPzm2+R1V3m+c0ovOSW7Is75f1vGBNxbsiLayq7nJpoQ8AWCki\nS1PSxbYD2A9gpY/6PwPrO7DBx7FEGXEKFlEATBDdDmvls3cd22thBYZH7PnzPsquBbBRVbeleb4O\nVqAYAdDhDOQishbAcgCbEcE1nE3imQ1u5+bhvCpUdVxEhmH9QNoBoNdtfXHzOSmA1aq6IaU8AbDV\nuUawh/O4A9aPjce9HkuUCbu7iYJht/CuSHlpsp5tAfCYScDixw7z7yomm1o/gH9T1XtSg7Dp9m1C\nllZ4iDbB5dw8nte42dwMK6B3AjgjIvtduqFVVZ8C0G7fOjCt8d58TsL0UGyI6e0IChGDNFHh2d3f\n7V4PNEHGdWlNkw61F8D2LIuG9CDPIFRAV3V1+z0vVT2nqn9n0sVWwgryPWmO7cZU/veKgH7A7ID1\no4MoMAzSRIXXAquLdcDHsRtg3ed20wNr0ZB/zFLG6QxlhCZDC9bzeZmR3cn7yXbAxpVriTs9BuAB\nEalBcKtC7QVwT0BlEQHgwDGigjILq3TCGjz2Py7P24tOjMEKKHtUdZdjl/Ww1sB2K3clgI4cqjEQ\n0UUV1sN6b5LyPK+HYQVKu6xmXP0jQADrNoSIjAHYYoJ58jm/TJl5TyUjcmKQJgqOANgkIiOwRnXf\nA+AUgHa3UdUishHASlW1W19dIjIsIpoyhcitpWePdk7XnZsU0QANuI/qzue8elIWtamDCfYmeD4K\noElE7CVLH7EPND+WNgKoM8sPPumzC/yIy8hxIt84upsoACKyDtZI5Hpn8BCRbgD9bqN+ReQ0gLud\nAVxEOgE0qepXzePLqnrVbSlz7Cl7uc5SY4Jie+qo7hic1yuw7qVz3XIKBO9JExXWVgCPmik6SaZb\ntwJWy61DRO53TA8aNPtUIH2Sj0p7Py8cWbe6zf1Y53OvmBHRqXXd7DgmXbaurPukWA/3Eeu+zitC\nxgBUhV0Jig92dxMVkLlPCVhBydnlbd+7fMZnt+oYpqZ9pSUim52teFXdKiIPw2rtHU3Z/ZU0c7l3\nwrqnnmlQVC77OKVLYOLrvIjiii1pouKoS3ncj6n7pq7cknGkHJ/2WGcxLtsSAOqdG0z38840ZbQj\ne+s2l32cr5VuSlg+5xUFlQhutDgRgzRREYwBaHVuMGlCz8Bl7rS9KIR9bJoEGY/CSshRk+5FRaQz\nTct4EI6pSaZbvSLDALNVyD7POpd9bOm6ugGf5xWhbvwqMEhTgBikiYKxAFbLuNLluScBVDoyXNnT\nizYA2OYMwqaV6ZxP3Y+UAA8AqroX1lzfntSsWmKt+PQQHKOXUyRwZWu1I8tAp3Zkn2edyz622nQ/\nCPyel1k8A0jfjb/MZYT9TutQvcdkIXOTyz5Oac+NyA/ekybKk2n52lOHukRkh/OCbu4DV8AaQDYC\n0+JU1b0mKD8lIn2wWtyJlGDSA2se9VVTuFR1m4jsMa+psILvKVNGpjzh+zE1NWkdrB8R6c6tGVYr\n+9WU7RV2d3wu+zi2rUSWbvE8zsvuxk/ORy9yN34trB9VRIFhkCbKkxnAlHEQkyNhRur2g8i8elKP\n+ecanExg9LqyVQJWy74GwJksA9fSBdUNsFZ+ynUf21UJTNz4PC+7G/8pIJRu/HbkML+byAt2dxNF\nmGmJaqZ7tD7KPACTeCWHpStXwcpznWS660c87mPLuCxlnsLuxt8Ea648UWDYkiaKvgdgzbd+IMAy\nBwD8NN2TphW6EVYrecQMmmqA1fVeq6qfy2WflDLzXm0qi7C78bdHdKUxKmEM0kQRZ+ZaD4vIHTm0\nfHMtM+OaySYIZezGz2WfFDl1dechlG5882Pl7nS3NIjywe5uohJgBkzVpZmOVSrmF3Lkc4jd+Oth\nrRlOFDjm7iai2BCR/bAWLXFtRTu66DthdYcPIH03ftp9Cn4iRAaDNBERUUSxu5uIiCiiGKSJiIgi\nikGaiIgoohikiYiIIopBmoiIKKIYpImIiCKKQZqIiCiiGKSJiIgiikGaiIgoohikiYiIIopBmoiI\nKKIYpImIiCKKQZqIiCii/gxHAOuFVyz2rgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fbe4e2dbe48>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(5,5))\n",
"xmax=0.2\n",
"x0=-0.1\n",
"ymax=0.1\n",
"flavio.plots.band_plot(fit['DeltaM_s'].log_likelihood, x0-xmax, x0+xmax, -xmax, xmax, n_sigma=1, col=1,\n",
" contourf_args={'alpha':0.5}, label=r'$\\Delta M_s$ ($1\\sigma$)')\n",
"flavio.plots.band_plot(fit['S_psiphi'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=1, col=2,\n",
" contourf_args={'alpha':0.5}, label=r'$S_{\\psi \\phi}$ ($1\\sigma$)')\n",
"flavio.plots.band_plot(fit['global'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=1, col=0,\n",
" label='global')\n",
"flavio.plots.band_plot(fit['global'].log_likelihood, x0-xmax, x0+xmax, -ymax, ymax, n_sigma=(2,3), \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-Bs.pdf')"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import scipy.optimize"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"collapsed": false
},
"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.7332560555767309"
]
},
"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