Skip to content

Instantly share code, notes, and snippets.

@DavidMStraub
Created October 13, 2016 20:31
Show Gist options
  • Save DavidMStraub/d122e350743d602928647f83c3799429 to your computer and use it in GitHub Desktop.
Save DavidMStraub/d122e350743d602928647f83c3799429 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": true
},
"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",
"%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": true
},
"outputs": [],
"source": [
"flavio.config['implementation']['B0->K*0ll subleading effects at low q2'] = 'B0->K*0ll deltaC9 shift'"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"observables = [\n",
" ('<dBR/dq2>(B0->K*mumu)', 1.1, 2.5),\n",
" ('<dBR/dq2>(B0->K*mumu)', 2.5, 4.0),\n",
" ('<dBR/dq2>(B0->K*mumu)', 4.0, 6.0),\n",
" ('<dBR/dq2>(B0->K*mumu)', 6.0, 8.0),\n",
" ('<FL>(B0->K*mumu)', 1, 2),\n",
" ('<FL>(B0->K*mumu)', 1.1, 2.5),\n",
" ('<FL>(B0->K*mumu)', 2, 4.3),\n",
" ('<FL>(B0->K*mumu)', 2.5, 4),\n",
" ('<FL>(B0->K*mumu)', 4, 6),\n",
" ('<FL>(B0->K*mumu)', 6, 8),\n",
" ('<FL>(B0->K*mumu)', 6, 8.68),\n",
" ('<FL>(B0->K*mumu)', 4.3, 6),\n",
" ('<S3>(B0->K*mumu)', 1.1, 2.5),\n",
" ('<S3>(B0->K*mumu)', 2.5, 4),\n",
" ('<S3>(B0->K*mumu)', 4, 6),\n",
" ('<S3>(B0->K*mumu)', 6, 8),\n",
" ('<S4>(B0->K*mumu)', 1.1, 2.5),\n",
" ('<S4>(B0->K*mumu)', 2.5, 4),\n",
" ('<S4>(B0->K*mumu)', 4, 6),\n",
" ('<S4>(B0->K*mumu)', 6, 8),\n",
" ('<S5>(B0->K*mumu)', 1.1, 2.5),\n",
" ('<S5>(B0->K*mumu)', 2.5, 4),\n",
" ('<S5>(B0->K*mumu)', 4, 6),\n",
" ('<S5>(B0->K*mumu)', 6, 8),\n",
"]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"bins = {}\n",
"bins[0] = [(1, 2), (1.1, 2.5)]\n",
"bins[1] = [(2, 4.3), (2.5, 4.0)]\n",
"bins[2] = [(4.3, 6), (4.0, 6.0)]\n",
"bins[3] = [(6, 8.68), (6.0, 8.0)]"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def bin_name(bin_list):\n",
" q2min = min([b[0] for b in bin_list])\n",
" q2max = max([b[1] for b in bin_list])\n",
" return \"{0:0.2f}--{1:0.1f} GeV$^2$\".format(q2min, q2max)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def observables_bin(i):\n",
" return [obs for obs in observables if obs[1:] in bins[i]]"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"fit_parameters = ['B0->K*0 deltaC9 c_0 Re', 'B0->K*0 deltaC9 c_- Re']\n",
"nuisance_parameters = [p for p in flavio.default_parameters.all_parameters if p not in fit_parameters]"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"fits={}\n",
"for i in bins:\n",
" fits[i] = flavio.statistics.fits.FastFit('low q2 subleading fit bin ' + str(i),\n",
" par_obj = flavio.default_parameters,\n",
" fit_parameters = fit_parameters,\n",
" nuisance_parameters = nuisance_parameters,\n",
" observables = observables_bin(i),\n",
" input_scale = 4.8,\n",
" )"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 11min 8s, sys: 312 ms, total: 11min 8s\n",
"Wall time: 11min 8s\n"
]
}
],
"source": [
"%%time\n",
"for fit in fits.values():\n",
" fit.make_measurement(N=200)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"CPU times: user 18min 44s, sys: 9min 17s, total: 28min 2s\n",
"Wall time: 41min 17s\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAGZCAYAAAAggY4CAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzsnXl8W8W597/nSLLkVbZsZ98sZw/Bke2EENbEsVkLJXYU\nKFyWEi8plMIt1M593957+97SrLSlt21QFCi0tJDEDmULENsUKBRIsJQAJYTEVjayepES77Z03j8c\nq1G8yYlseZnv5+OPrZk5M49s+fzOPPPMPJKiKAgEAoFAEAzkYBsgEAgEguGLECGBQCAQBA0hQgKB\nQCAIGkKEBAKBQBA0hAgJBAKBIGgIERIIBAJB0BAiJBAIBIKgIURIIBAIBEFDiJBAIBAIgoYQIYFA\nIBAEDXWwDQgUkiSlAcnnXs4F8hVFcQTRJIFAIBD0wJAQIUmS9ECyoijrzr3OBIqByUE1TCAQCATd\nIg2FA0zPzYJ2KIqiOvdaD9QARkVRDgbTNoFAIBB0zZCYCSmKUipJUsp5RXPbijsKkCRJscANwEGg\nsV8MFAgEgqGBDpgEvKMoSlUgOhwSIgSgKMru817+BMjpoukNwJ/73iKBQCAYstwN/CUQHQ0ZEWpH\nkqRsYIuiKM920eQgwIsvvsiMGTP6za7+5LHHHuNXv/pVsM3oM8T7G7wcOXKEJ598kmeeeSbYpvQZ\nQ/Hvt23bNj755BMeeOAB7rnnHjh3Hw0EQ0qEzq0NVSmKsq2bZo0AM2bMIDk5uZtmgxe9Xj9k3xuI\n9zeYiYqKIjIycsi+Pxh6f7/nnnuOjz76iDfeeIN9+/a1FwdsKWPI7BOSJCkZoF2AJEnKliQpKrhW\nCQQCweDl+eef54UXXuD1118nPDy8T8YYEjMhSZISgM8ARZIkAAmoURTFGlTDBAKBYJDywgsv8Nxz\nz/Hmm28SERHRZ+MMCRE6tyl1yMzqBAKBIJi8+OKLWK1Wtm/fTmRkZJ+ONSRESODLXXfdFWwT+hTx\n/gY3t956a7BN6FMG+9/v/fff5/e//z1vvfUWUVF9v6IxJDar9oZza0dlZWVlQ2rxUDD4OXz4MJWV\nlcE2o085cuQIAOPHjw+yJYLOiIuLY/To0TQ0NHQqQDabjZSUFIAURVFsgRhTzIQEggHA4cOHmTFj\nBvX19cE2RTCMCQsLY+/evUyYMKHfxhQiJBAMACorK6mvrx/S+9cEA5u9e/dyzz33UFlZKURIIBiu\nDOX9awJBZ4iIMoFAIBAEDSFCAoFAIAgawh0nEAgEQxCrtW2vvs1mY/HixWRmZgbZos4RIiToUxRF\nQamtxVNVhae6BsXjQY4IRwoPRwqPQA4PQ9Jqg22mQOA3LpcLvV4fbDO6xW63YzAYvMIjyzJOp7Nf\n9v30FiFCgktCURTcR47Q9PEntHz5JZ7KStzf7MNTW4untg5PbS20urvvRKVC0mqRdFrk0FBUcQY0\nV16F2mhEbUxAbTQijxjBuSOZBIOQ0tJSLBYLW7Zs8au91WpFkiQURcHlcvH444/3qr4n1q1bB8Cu\nXbswGo2sXr26R/vT09O9rxMTEykuLmbSpEl+jVdQUOC1V5Ik5s6dy5IlSygoKOhxbGh7vxaLBZvN\nRnp6OmvWrGHOnDneervdztKlS6mpqSEnJ4fU1FRKSkq8ImQ0GqmoqPC55mKIjY0lOzub++67bwmQ\nepHdbDz/hRAhQa9pPXSIpo/+QWPhFpq/KcdTUwOShGrkCGR9FHJkBJpRI5AjItq+Itu+pIhwJEnC\n09SE0tiE0sl3T1097tOVNBQW4q6qhnObqSWtFtXIeDSmZDTTpxMy/wpCkpLELGqAU1paSnFxMU6n\nE4fD4dc1VqvVR1iKiop8btY91ffEhW3NZjNms7lbgXQ6ndhsbXszo6Oj/RYfm81GTk4Oa9euZdGi\nRd5yh8OB2WzGbrf7ZXd2djZmsxmDwUBubm4HMTGZTKxdu5bo6GjvOO2i6XQ6qampuWQB6iuECAn8\nouWbb2h47XUaNm+m9dhxkCTUE8ajSzWhmZJIyBQjcoBP2VVaWnBXVtF68hTuk6dxnzxF61f/pPGd\nd1DWNIFGg8Y4Cd2t30F75XxCTCYknS6gNggujbS0NNLS0igqKqKsrMyva9asWUNJSYn3dWZmJtnZ\n2d6bdU/13eFyuSgpKeHMmTNe19TKlStJSUnh4MGD3YpLb8SnncWLF1NUVMTChQt9yhMSEli2bBl2\nu93vvvR6PZmZmVgsFpYsWdKhvry8nCeeeML7uv395eTkUFhY2Cu7+xMhQoJucVfX4PpBHg0ffowU\nqkM7exbht91MyIypyKGhfTq2pNGgHj0K9ehRPuWKx0PrkW9p/uYALd8coHbDM5x96pegVqNJmIju\nllvRpaWhMc0RLrxBhsvlwuFwdLjZO51Odu/eTUJCQrf1/jztOxwOH9eU0WgEoKKiotci0x35+fnE\nxsZ2EKB2MjMz/Z69tZObm0tGRoaPiLYTExPTof26detYsWJFlzYMBIQICTpF8Xio/8tLuP7f/wOP\nQuSyTEKvXYCkDv5HRpJlNBPHo5k4HtIXtonSt8dp+eYAzd/sp9Zq5eyvn0aOiSb0jjsIvfkmQubN\nQ1Kpgm26oAcqKio6LY+OjqaiooKuzrpsr+9JhPR6PVVVVT5l5eXlSJLkFaOuKC4uJjY2lqqqKsrL\ny3sUkKKiIp91pM5YuXJlh7KCggLmzZtHeXk50dHRZGdne+vS0tLQ6/Vs3LjRZx2sqKgIs9ns009h\nYSHJycksXLgQm82GLMsD0iUX/DuKYMDRvHs3NQ89ROvBw+gWXEHkktuQo/r2OPdLQZJlNOPHohk/\nlrC061A8Hlr2l9No30PDtm3UPfcHpMgIQm+7jdCbb0K7YAFSSEiwzb4kPA0NtB440KdjqCdP7vPZ\n7oVUV1d3Wm4wGKiuru4yKq29/mKwWCykp6d3OwtKTEwkMTHRexO3Wq3k5eV1m6a8oqKCxMTEbse+\n0K2WkZFBQUGBd13HbDaTmJjos55kNpuxWCw+IuRwOHxmRna7HbPZTExMjDd4w+3uIUAoSAgREng5\n3/WmHjeGmPxHCUns/ulwICLJMiHTphAybQqKeQmtBw/TaNtD0zvvUP/nvyCFhaJNmk3Ejx8nZP78\nQemyaz1wgNM33tynY8S/vZ2Q2bP7dIxgY7PZePfdd71BB11x4Qxi8eLF5Obmsnbt2oCFPdtsNkpL\nS9mxY4e3LD09HYvF4iNCubm5WK1Wr/vR4XB0mMWZTCY8Hk9A7OprhAgJ2lxvL72M62c/+5fr7bqr\nhoT7SpJlNMZJaIyTiMi8jdajx2iy7aZxZxmVWWZUI+IJf/BBwpZmoRo5Mtjm+o168mTi397e52P0\nNwaDodPy6upqDAZDj/VFRUWdPlS0h0ZfOPNYuXIlNput14nbEhISALp1ARqNRsrLyzuUt0cMtgtf\nSkoKq1atoqysDKPRyLZt27xuR5fLxdy5c32uN5lMJCQkYLFY2LBhAyUlJT4uu8GGEKFhjtLYSOUd\nd9D8+ZfoFlxBxJLvoBqAG9oCgSRJXrdd+G030/LNARo+/Jgz65/izJq1aGfPJPzRx9AtWjQg1r66\nQw4NHZKzlPYn+gsX3p1OJ4mJiV6XWWf1RqOxV2seeXl5WCyWHgXI5XKRkJCAzWbzju9yuXrsPysr\ni6Kiog7l7RGDeXl5HWY+0NFF1xm5ubmsXr2aDRs2+GXLQEacHTeM8dTXU3nrrTTv3Uf0D3PR33/3\nkBWgC5EkiZBpU9A/eC/x6/6HyDszcdc4qX7gQU5cnoRr1WpazyVgE/Qfer0eo9HYYX1HkiSSkpK6\nre+NAFmtVgoKCryiYrfb2b17d5ft586d67Nm1B7M0N2Yq1evprq6mm3btvllU2pqaqeBGZ3tr8rJ\nyaGmpoaCggIWL17sV/8DFSFCwxRPQwOVN99MS8VBYh7JQzt7VrBNChpyWBhh119D7P/9CYb/+wRa\nUxJ1zz3HyQVXU3XHd2n5am+wTRz0XBiR1o7dbicvL8+nLD8/H4vF4n1ttVpZs2aN3/U9UVhYiNPp\npLy8nNLSUgoLC7FYLN5Z2IU26fX6DlFuq1evZu3atT2OVVJSQn5+Pps2bfIpd7lcHQTHZDKRlZXl\n09blcnW6l0iv13v3IA3EiLfeMLB9DoI+w5m9nJaDh4n58cOEGBOCbc6AQTNhPJrvjScy83YaPvqE\nuh2lnErPQJt0GVFP/oIQkynYJg4q7HY7mzdvprCwEIfDwYoVK0hJSWH58uVA25pKaWmpzzXLly9n\n/fr1bNq0iZqaGqqrq1m1apXf9d3hcrkwm80d1o0kSfJGunVm0+OPP+496qe8vJyMjAzve+gOk8nE\n/v37KSgoYNmyZRgMBvR6PXFxcezYsYP169f7tN+8eTMrV65k/fr16PV6JEnqcpy8vLwuQ9oDwZ49\ne4iKivKuf23fvp1Tp07hcDioqqrqNLz0xRdfHDdixIiGjIyMzp86OkHqKu5+qCJJUjJQVlZWNmyT\nhzX+/UOq7ryLyDuzCFt0bbDNGdAorW4ad35G3VvFuE+eImTGNKJ+/nO08+cHdBybzUZKSgrD+XMp\nCC7nfwYPHTpEZWUlixcvJiEhga+//prm5mb0ej35+fkkJib+8cknn/zo/Ot/8pOfzDl9+nTEXXfd\ntbcHEfI5O06444YZnjNnqPnBDwiZPpXQ668OtjkDHkmtInTBFcT+7D/QZ9+P5+xZKjOXcuq662h8\n770uN08KBIOZO+64g6SkJO/rsWPH+qSd12q1HTYdrV27dvf8+fOP9nYsIULDjJrlD6LU1xN13/eQ\nZPHn9xdJltHNTcbw03yiH8pGaXVTdfe/cXrRIpr37Am2eQJBnxIZGYlGo+HDDz/k2LFjzJkzpzJQ\nfYu70DCiobiExo8+IdK8BFVs5/stBN0jSRLapNkYVv470T/MxVNbx+mbb6Xqju+KaDrBkOfqq69m\n0qRJbN++PWALySIwYZjgrq7B+ciPCLlsJrqrArueMRyRJAnt7FmEzJxOwz8+pe617Zy8+hoiHnyQ\nyEd+iBwdHWwTBYKAsXPnTqKioggNDUUOsAdFzISGCTXfvx+ltZWoe+8alMfUDFQklYqwaxYQ9/Of\nEn7zDdQ9/zwn5l7BWctGlKamYJsnEASEOXPmUFFRwdtvv823337LokWLDh86dEj76KOPXmxiOy9i\nJjQMqH/tdZp22Yh68F5U0QM7LfFgRdJqibj1RkKvXUDd629x5udPUvu73xH95JPobr1FCL9g0DFv\n3jzvzyEhIdx8880cOnSIsrIyxo4d2zBx4sSmX//615+df01ubu7B3o4jZkJDHM/ZszgffwJtchK6\neSnBNmfIo4qKIuruZcT+VwHqsWOozltB5c034T5+PNimCQQDEiFCQ5zGHcUodXVELr1DPI33I+rR\no4h5OAf9iuW0HjzCyWuuo27zFhHSLRBcwLB1x7V6WoNtQr9Q/9yzaIyTRDRckNCZLidkaiJnt7yC\n899/TP1zz2L4w/OoxowOtmmCYYDVagXaNqIuXryYzMzMIFvUkWErQv/50f/hvogHuHHSTcToOqbF\nHQp46upo+udeIm7v27wzgu6Rw8PRP3APupQ5nPnTy5y85lr0T/4PYcuWidmpoM+w2+0YDAav8Miy\njNPpDFj+o0AxbEVoYvRECr/ZwpZ9L3NZ/Ezun7mcxOj+z5/SlzSWvgstLeiSB/cBh0MF7eWXEfvf\n/8HZra/g/PET1D/3LDF/eAH12DHBNq3PaT93bdeuXRiNxh5TY0PbU7wkSd7MoOdnEvWnvjdkZGR0\nSKlwIS6Xiy1btgBt58c5nU7WrFnTZbbXCykoKPDaK0kSc+fOZcmSJRQUFPj9+7BYLNhsNtLT01mz\nZo3P4aV2u52lS5dSU1NDTk4OqamplJSUeEXIaDT6lQK9vxm2InTtuKuYOGsCthO72XX8Mx5770dM\niBrPXdPuYf7oK1HJgz+hW/2zVtTjx6GKiw22KYJzyOFh6O+/2zsrOnXttUSvXwcTJwbbtD7jwpus\n2WzGbDZ7b+idYbVafYSlqKjIp5+e6ntDYWFhhwNLOyM/P5+8vDzvTTwvL4+lS5f2KF42m42cnBzW\nrl3rkyHV4XBgNpux2+1+2Z2dnY3ZbMZgMJCbm9tBTEwmE2vXriU6Oto7Tvvp306nk5qamgEnQDDM\nAxNC1aFcNe5KHkl9CPP0TCQk1uxaxX3v3MO2/UXUtdQF28SLRmlooPmLf6JNTuq5saDf0c6eRex/\nr0SbNJuahx/BVVAQbJP6BJfLRUlJCWfOnPGWrVy5ksLCQg4ePNjldWvWrCErK8v7OjMzk40bN/pd\n3xv7OsvX0xkOh4OSkhLv68TERL/Ea/Hixaxbt85HgKAtO+uyZct6Za9eryczM9MnlcX5lJeX+4zT\n7nrLycmhsLCwV2P1F8N2JnQ+siQzI246M+Kmc7z2BDuP7eKPXz3PX75+kVsSbuU7ibcTFxoXbDN7\nReMHH6A0NQtX3ABGDgsj6vv/Rsi0KTQ9/+dgm9NnOBwOHzdQe96eiooKn0Rx7bQLw4V1TqeT3bt3\nk5CQ0G19b572t27dSk5ODgV+PAS88847Pq/Ly8s75Bm6kPz8fGJjY1m4cGGn9ZmZmb2eveXm5pKR\nkdEhuyxATEzH9e1169axYsWKLm3oDVVVVVitVqxW6zZFUWyX3CFChDowOmIUt0/9DosmXs+nx3ex\n3fEmr5b/ldnxl7H8slwmRk0Ktol+Ub/Rgmr0KNSjRwbbFEE3SJJE6NVXEiXL8Av/E7MNFvR6fYeE\ndu1ZSdvF6EK6ypETHR1NRUVFl2Hu7fX+ipDdbic19eI2/LfnHDp/ZtQZRUVFPQrVypUrO5QVFBQw\nb948ysvLiY6OJjs721uXlpaGXq9n48aNPutgRUVFmM1mn34KCwtJTk5m4cKF2Gw2ZFkecC45IUJd\nEKmNZPGkRVwz7irKTtj59NhOfvjuQ8yMnc6KpB8OaDFSFIWmL74i7Lqrgm2KwE/UI3o3025sdnOw\nsm/dxZPiwtGFBH5t1GKxkJ6e3uksCOiQursdg8FAdXV1l4EA7fX+UlZWxvLly3G5XH5fA23rURs3\nbsRisTCxh7W8iooKEhMTu22zZMkSn9cZGRkUFBR43Wpms5nExEQfN5vZbMZisfiIkMPh8JkZ2e12\nzGYzMTEx3uANt7tDBoagI0SoB7RqLQvGzeeKMXPZc+oL/n7kQ3747kPMipvJQ0k/YlzkuGCb2AFP\nZSVKXR3qieODbYqgjzhYWcf9lo/7dIznc69k+pjAhvPabDbeffddbLaAeHIumqKiIr8yo3ZGdnY2\n2dnZZGRkUFZWxhNPPBEwu2w2G6WlpT7BDunp6VgsFh8Rys3NxWq1et2PDoejw8zSZDLh8XgCZltf\nIUTIT1SyiuRRc0gaMZvdp/bwwZGPeKg0jzkjL+fhpMeID4sPtoleWvcfAEA9SrjihiqT4sJ5PvfK\nPh8j0KxcuRKbzUZkZGSXbQyGzjdWV1dXYzAYeqwvKirqdP9Ve2i0yWQi+rxTzi/2FIv8/HzS09NZ\nunRpl7M6o9FIeXl5h/LS0lKKi4u9YpySksKqVasoKyvDaDSybds2r10ul4u5c+f6XG8ymUhISMBi\nsbBhwwZKSkp8XHaDCSFCvUQlq0gZlUzSiMv57LiNvx/5iJziB7nVeBtLpy4lShv8A0JbDxwAWUYV\nP3CEURBYdCGqgM9S+pq8vDwsFku3AgT/Cly4cOHd6XSSmJjoveF3Vm80Gntc87BarVRUVHjXc2pq\naoA2gWzfu3MhLpeL7OxsNm3a5B2z3c6SkpIuZ1VZWVkUFRV1KE9LSyMtLY28vLwOMx/o6KLrjNzc\nXFavXs2GDRt67VIcSAgRukjUspr5Y+dhGpnEx8c+5a2Db/LWwTfJmrKU2xK/S5gmLGi2tRw4gCo+\nDkk9+Pc6CYYGVquVgoICr4DY7XYkSepUMPR6PUajkerqah+RkSTJm3K6q3p/Ft0vnDHY7XasViur\nVq3q8pr2QITzx3Q6nV5bumL16tVYrVa2bdvml7CkpqZ2GpjhcDhISPDNI5eTk0N+fj4FBQXceeed\nPfY9UBnW+4QCgVat5foJ1/JI6g9IHjmHzfte5vs77uO18ldpcbcExaaWj/+BevSooIwtEFxIYWEh\nTqeT8vJySktLKSwsxGKxeG/edrudvLw8n2vy8/N99sJYrVbWrFnjd31v6Mwdd6FNJpOJnJwcH7fb\n5s2bSUlJ6bD/50JKSkrIz89n06ZNPuUul6uD4JhMJrKysnzaulwu7HZ7h371ej2LFy+mqKhowEW8\n9YZhOxNSCOxpxuGacG4wpjN/7BW8f/gDnv3CytZvNvP9Wcu5bvz1yFL/6X3riZMibYNgQOByuTCb\nzR3WaCRJ4plnngH+Ncs4n+XLl7N+/Xo2bdpETU0N1dXVPjOVnur9xWq1snXrVgCWLVtGbm4uixYt\n6tSmgoICn6N3XC6XX5tVTSYT+/fvp6CggGXLlmEwGNDr9cTFxbFjxw7Wr1/v037z5s2sXLmS9evX\no9frkSSpS3dfXl5elyHtgwVpuB0tL0lSMlC2OP933JN1GbMmqvrkEMnK+krePfQee6v2MSp8JCuS\nHsY0Ijng41yIp66O41OnE3X/3YQuuKLPxxMEBtuBcq549AnKyspITu77z4lAcCE2m42UlJRuP4Pt\nbYCUQG1WHVLuOEmS/D6nXC1LbHi9gfVb6/nmaODTOsSFxWGekcX3L78PjUrDf/3jpzz63sM4XP4d\nEXKxtJ57KhLuOIFAMBgYEiIkSVKmJEnZwFZJkvwKGbo7TcvDt4cC8PQrDfyxuIHahsDH1I+PGscD\ns+9l2YwsXE1nePRvP2Tj589Q21wb8LEA3MeOAYj8QQKBYFAwJNaEFEUpApAk6Rl/r5GQmD5BzbTx\nKj75qoVXPmriy4Ot3HGVjvkz1AF10UmSxPTYaUyJmcynx3ax49A7vHuklOWX5bBoQlpA14uUhoa2\nMbXagPUpEAgEfcWQmAmdR6+VQ5YkFswK4T/vCWfmRDUvljby620NHD7lDngqZpWsYsG4+TycnEdi\ndAK/sf+ah9/No9zZcTPbxaI0NLb9oBkSzxcCgWCIM9RE6KKJDJO5PyOUH94eiqvOw5rN9fzipXqK\ny5pw1gbWTRepjWTJtO9y3+x7aHI38+/v/QjLng0BcdEpDQ2g0SDJ4k8rEAgGPuJx+QKmT1Dz07vD\n2XvYzc59Lbz5aTOv/qOZqeNUzJuuYU6iGl1IYFx1k/QTyZ3zIJ8e20Xx4R28d/Rv5MzO4/rxCy/a\nHag0NCCFaAJin0AgEPQ1Q02E/PafPfXT3xAR5Xs21g1L0rlxSToqlcRlCWouS1DT0KRgL29l59ct\n/KmkkZffgySjmnnTNEwZqyJEc2mC1O6iuyx+JjscpfzK9hRFB7bweEoBk/STet1fmwiFXJJNAoFA\n8NJLL/HSSy/5lPXF8UBDTYT8VoQf/88jTE+a1mO7UK3EgpkaFszUUH3Ww659Lez8upXPvmlAlmFc\nnEzCKBWTRqkwjlIRGyVd1CwmShtF1vQ7SHbO4a3yd/jR3x7mtsTvctf07/XqCCClsVGIkEAguGTu\nuusu7rrrLp+y8/YJBYwhIUKSJKUBybTNhFZKklSsKMq7gR7HEClzQ6qWjJQQjld7KD/mxnHCzd7D\nrbz/edsRPZGhEgmjVCSMlpk0UsW4OBVhOv9FyRidQJ4pm4+Pfcp2xxu8e6SYFUkPc9WYq/0SN+GO\nEwgE7VitVqBNPBYvXkxmpt9bKfuNISFCiqKUAqXAuv4YT5IkxsSqGBOr4prZbWW1DR4cJzw4TrQJ\n09u7mmk6d3RcTITE+HgVE0fKTBypYuKI7oVJJau4etwCZsfN4u2KHazdtZopMZP5j3k/JbaHNONK\nQwOSRoiQQDDcsdvtGAwGr/DIsozT6eyQEjzYDAkRGghEhMrMTpCZndD2K/V4FE7UeDhW5eHoaQ+H\nT7kptjXT2NzWfkS0xIQRKsbGyoyIkRkZIxOvl1Gr/iVOep2eZTOX8nXVPt4sf5sVJbmsSHqo+8AF\ntRplECSyEgxfMjIyOqQu6Ayr1epzTtv5WUT9qfeHgoICJk+ejKIoPjfsi7Wpp7Har5UkyZs2oqCg\ngNWrV/d4vdVqxWKxYLPZSE9PZ82aNT4Hl9rtdpYuXUpNTQ05OTmkpqZSUlLifU9Go7FX6c/7CyFC\nfYQs/2u2lDq1rcyjKJx2Khw66ebgSTeHT7n558FWGs4JkyxBbJTEyBiZEdEyo2Jk4qNl4iKnkH35\neEoPFfMr21NsP/ga/zHvP4nRdTwVQQoLQ2ls6sd3KhD4T2FhoV+HflqtVp+bfFFRkc/Nuqd6f8jI\nyGDjxo1MmjQJu91Oampqt+mvL3ZMm81GTk4Oa9eu9Tlx2+FwYDabsdvtftmdnZ2N2WzGYDCQm5vb\nQUxMJhNr164lOjraO056ejrQlnaipqZmwAkQCBHqV2RJYmRMm8jMm97mMlMUhbMNCidrPN6vU04P\nXzha+dsehfP3y0aGpREVN4n9nnf5/lvZzI26k9T4q4iP0hETHoIhQosUHoHS2BikdygQdI3L5cLh\n8O/sxDVr1niTzgFkZmaSnZ3tvVn3VN8TVquVlJQUb2oGk8lEWVnZJdnUFe3pFhYuXOhTnpCQwLJl\nyzpN09AVer2ezMxMLBZLp/mJysvLfdKNt7vecnJyKCws9Huc/kSIUJCRJImoMImoMJkpY33rWtwK\nVS4PzlqFmlqFmloPNWenUeUcR432XT6Rnue9w3/n7NGFKK1tEXQqZSb6m39KyD4D4Vo1Uz1OYjQK\nerVClEohSg1RauXcV9vPWrGvVdAPbN26lZycHAoKCrpt1y5WF6bMdjqd7N69m4SEhG7r/Xnaz8/P\n73BT7u66nmzq6tr8/HxiY2M7CFA7mZmZvZq9QVtG1YyMjA6ZZQFiYmI6tF+3bh0rVqzo0oZgI0Ro\nAKNRSYwyqBjVweumA5byz8q9vHngLaLj/sw0lZlI9yyw25D37aR8xj3UNbn5oiWS2vpWGlvcNLZ0\nvlakVkkKJ+zrAAAgAElEQVSEalRo1Sq0GhmdRsWkVifhKoUIFUSoFcJVEK5SCFNBmHzu+7nX4SqF\nULnNnSgQdEa7u8sfusqPEx0dTUVFRZfHabXX9yRCLpcLp9OJ0+n0Ro9VVFT0mFn1YsYsKiryusS6\nYuXKlR3KCgoKmDdvHuXl5URHR/tkg01LS0Ov17Nx40afNamioiLMZrNPP4WFhSQnJ7Nw4UJsNhuy\nLA84l5wQoUHMrLgZTIqawBvlb7Gn6gVmxc0kLSqBGWV/ZU/BI3h0oT7tPYpCY4ubhua2r/N/bmhx\n09TipqnVQ2OLm3KPnqZmN42tHprOCZinh7P0NCqZELV07rtMiEpGo5YZ1+oiVIZQWUGnglAZdLJC\nqKrtu06GUFXb9/afz28TIkEfpHwa1DS1NnK09mifjjEuYhxatS4gfZWVlbF8+XK/NjtWV1d3Wm4w\nGKiurkav13db3xPtglJdXe29uZeWlmI2m9myZctF2dTdWImJid3ac6FbLSMjg4KCAu+6jtlsJjEx\n0Wc9yWw2Y7FYfETI4XD4zIzsdjtms5mYmBhvIEV3a17BQojQICc8JBzz9Ey+rPyK7eVvc2jiAe65\nPAq5ob6DCMmSRFiImrCQi/uzt7o9NLV6aD731dTqPu9nD82tbppbFVrc/2rT4vZwTIpuK2v2eOta\n3G3t3B7/DrloFzSNSmoTOLWMRiUzvl3gVAph5wlYmKpt9tb+PVzlO5tTDXJRO1p7lMfe+1GfjvGr\n658mMXryJfdTVFTUZWbQYFBdXY0kST4zs7S0NNLT0zl48GAHl1t/YrPZKC0t9YkeTE9Px2Kx+IhQ\nbm4uVqvV6wp0OBzedOntmEwmPIMgUlaI0BBAkiRmx89ikn4C2z/bwjO5k0hyP0uS+yF0qvCeO/AT\ntUpGrZIJD2CWCI9HodndJk4t58Sp2fvzv4Ss7WeFllaPt31zq4dvpei2103/Er7mc3XdoVHJ6DQy\noSEqdBoVEz1niFRBpFohUqUQqYZotUKMRsFw7itSNXBmZOMixvGr65/u8zEuFYfDQXR0tPe1PyfT\nGwyd58Kqrq7GYDD0WF9UVNTpFob20GiTyQTQ4aYNbSLQmQj1NGZXGI1Gyss7npJfWlpKcXExNltb\nctKUlBRWrVpFWVkZRqORbdu2eX9XLpeLuXPn+lxvMplISEjAYrGwYcMGSkpKfFx2gwkhQkOIyJBI\nvhd9PeW/Wsvm+zVUtP4X10XnMEY3NdimdYksS+jkNiEIJIqieGdo7W7G9tlaU0vbLK6huc312NDi\n5lhzdNvPDe5z62duLpykyZJEhE5NhFZNuE7NdHcNI0M8jNIqjNYqjAxRiAtRUPeDUGnVuoDMUvqa\nkpISKioqvFFlNTU1QNs6SPs+mQtpF4cLF96dTieJiYlekeis3mg0+rXmoShKr/bMdGdTZ2LWTlZW\nFkVFRR3K09LSSEtLIy8vr8PMBzq66DojNzeX1atXs2HDhj45062/ECI0xFBGxjJ/p5NRd87n5chD\nvHHqKUxRN5OivxVZCuyNfiAjSRJajQqtRgWhvT9BQlEUmlo91De1UtvUSl2Tm7rGFmqbWqltbCvb\n1RyBy9VCQ/O//OySBFE6DfowDTM4w1idh4k6hQmhHibohl8k4oVP53a7HavV2m0QgF6vx2g0Ul1d\n7XPDlySJpKQkgC7r/RWVlJSUDiIkSRLJycm9tqm7MVevXo3VamXbtm1+CUtqamqnQRAOh4OEhASf\nspycHPLz8ykoKODOO+/sse+ByjD7lxj6eAx6FK0Gw4la7p19NwsnXsfuM2/xjvMXnG2tCrZ5gwZJ\nktBpVBgitEyIDWfGmChSjbFcP2Mkt5rGcuf8iSy/PpEf3zSdn9wyndyFidw5fwI3zh7NrHF6InUa\n/qlEsrkylP8q1/LAl6GkfRbGLV/que/gCP73kIZXT6nYfUamtjXY77b/6MwdZ7fbycvL8ynLz8/H\nYrF4X1utVtasWeN3fU+sXr2a4uJin+uzsrK8s6yLsakrSkpKyM/PZ9OmTT7lLperg+CYTCaysrJ8\n2rpcrk73Eun1eu8epIEW8dYbxExoqCFJuEfFoTpWiSzJXDv+ahL0kyja91e2nfwZ18TchzEssKfg\nDndC1Crio1TER3UeSdbQ3EpVbTOVZ5uoqm2israZ4rPh1Jxs9uYeUTkj+8/gIGG1Wtm6dSsAy5Yt\nIzc3l0WLFlFRUdHhFIXly5ezfv16Nm3aRE1NDdXV1T6zp57qeyItLQ2Hw+HdsyRJEps3b/bWX4xN\nXWEymdi/fz8FBQUsW7YMg8GAXq8nLi6OHTt2sH79ep/2mzdvZuXKlaxfvx69Xo8kSV0GduTl5XUZ\nPj5YkAKdwnqgI0lSMlD2YslzfqVyGIzE5D6JEq7D+csfe8saWxt548B2/lm5l+SRc5ijeRC1LFI+\nBJNWt4fqumZOOBvZs9vOlp/dS1lZWZcuIYGgL2lP09DdZ/C8VA4piqLYAjGucMcNQdxj4lAdr/Qp\n06l1ZE67g+9MvoXPT3/Ja1U/o6q5b/eZCLpHrZIZEaXj8gnRXDml+9PRBYKhyrAVoSOVA2/TVqBw\nj4pDdfx0h3JJkkgeNYecOQ+ikmT+evIX/PPs3/wKmxUIBIK+YNiK0Ob3GvnSMTRXhN1j4pBrzkJD\n56dpx4fFsXzOA6SMMvFRzct8UPcUje7afrZSIBAIhrEIGUepsbzZwMdftQTblIDjHh0PgOpEZZdt\n1LKamxJv4M4ZSzl85iivnP5vjjfu7y8TBQKBABjGInT7VSFcOVPDi6WN/LG4gcbmoeOSco8bAYD6\n0Ike206LnUqeaTkGXQxvnHqKL86UCPecQCDoN4atCKkkie8t0nFvuo7d5a2sermOgyeGxjqRZ2Qs\nnsgw1AcO+9U+ShvFvbPvZv7YeXzs3MrHjf9Li0ckxhMIBH3PsBWhdq6YrmHlXeGE6ySeKqrnnc+a\n8Ph5qOaARZJonTwB9X7/RAhAlmQyEhaTOe0O9lV9w5vVP8fVcqoPjRQIBAIhQgDE62V+nBnGYlMI\nr3/czP++2oCzduCfPtsdLVN7J0LtXBY/k+VJD9DqaeWvp37OoYbP+8A6gUAgaEOI0DlUKonbF2h5\n5I5QTtZ4+GVR/aAWotbJE1AfPAbNvQ+8GBEeT/acB5ikn8g7p3/HZ87XUZTB+7sQCAQDFyFCFzB1\nnJrHl4bh9sDvXmugvnFwuuZap05AcntQO769qOt1ah3LZmSxcOJ12M68wXu162ny1AXYSoFAMNwR\nItQJhkiZH94eiqvOw4Y3GmhuGXxC1Dp5PADqb3rvkmtHkiSuHX81d8+6kyNnv+W1yv8nTlkQCAYR\nVqsVq9XKihUrOk0pMRAQItQFowwqVnwnjKOVbja91YDbPbiESIkIo3VMvN8Rct0xOSaRnDkPolVp\nefXUKvbXfRoACwUCQV9it9sxGAxkZ2ezYcMGli5dypkzZ4JtVgfEKdrdkDBKRc7NoWx4vYEXihu5\nN12HehDlhW6dMgHNJcyEzidGF833L7+PNw68xd9OP8fp5kPMj84cVjmKBJdGQUEBkydPRlEUDAYD\nmZmZ3ba3Wq1IkoSiKLhcLh5//PFe1fdE+/XtJ2IXFBSg1+u7bL906dJOZxMpKSns2rWrx/EKCgq8\n9kqS5E3qV1BQwOrVq/2y12KxYLPZSE9PZ82aNT4pHOx2O0uXLqWmpoacnBxSU1MpKSnx/p6NRmOv\nEvn1F0KEemDGBDXfv0HHH95p5JmmBpbfFIouZHAIUcvMBML/8jZ4PCBf+qRXo9Lw3anfYUzkaN6p\nKKZB5eDqiB+ilcMCYK1gKJORkcHGjRuZNGkSdrud1NRU3O6u9+VZrVYfYSkqKvK5WfdU3xPr1q1j\n6dKl3vxBLpeL/Px8nnnmmS6viYuLo7S0lJiYGG/Z5s2bWbZsWbdj2Ww2cnJyWLt2LYsWLfKWOxwO\nzGYzdrvdL7uzs7Mxm80YDAZyc3M7iInJZGLt2rVER0d7x0lPTwfaMsDW1NQMOAEC4Y7zizmTNfzg\ntlAcx908/Uo9Z+sHR6RYy5xpyK5aVAePBaxPSZK4Ysxc7rnsLr49e4w3qn6Os+VkwPoXDD2sVisp\nKSneG77JZKKsrKzba9asWUNWVpb3dWZmJhs3bvS7vieKi4u99kBbgrie8vKkp6ezcOFC5syZw5w5\nc0hISCA2NrbHG/vixYtZt26djwABJCQk9ChgF6LX68nMzPRJrnc+5eXlPuO0Z4HNycmhsLCwV2P1\nF0KE/GTaeDWPZoZRU6vwVGE9la6BL0Qts6egyBIhu78JeN/G6ASWz7kfgNdO/YKjjXsDPoZgaJCf\nn+99Im+nuxu3y+XC4XD4iAS0Pc3v3r27x3p/qK6uZt26dT5lktS9h+PC9NyrVq3q0QWYn59PbGws\nCxcu7LQ+MzOT6OhoPyz+F7m5uZSUlHS6vnP+LK2ddevWsWLFii5tCDbCHdcLxsereDwrjN++Ws/6\nwnoeWxLGyJiBq+NKeCitUyai2b2PhiWLer6gl8SGxrI86X4Kv36Ft049zYKYZcyKHJgf9KFGS1Mr\nzqOuPh0jepwejfbSbhEulwun04nT6cRqtQJtWUu7y0ja1YwkOjqaioqKLs82bK/3x+W0Zs0a0tPT\nKS4uxmKxeL/8xW63M2/evB7bFRUVdRDgC1m5cmWHsoKCAubNm0d5eTnR0dFkZ2d769LS0tDr9Wzc\nuNFHBIuKijCbzT79FBYWkpyczMKFC7HZbMiyPOBcckKEekmcXubHWWH8elsDz7xRzxNLwwnTDdw1\nouY5U9F+8kWf9a9T6/jerGXscJTw0bGXadJ+wxzNchGw0Mc4j7rY9u9v9ekYS355E/GJsZfUR7ug\nVFdXe2+kpaWlmM1mtmzZ0uk11dXVnZYbDAaqq6u7DB5or/eHtLQ0iouLycjIYPLkyWzdurXDzKo7\nVq1a1aX951NRUUFiYmK3bS6cYWVkZFBQUOB1q5nNZhITE33cbGazGYvF4iNCDofD636DNqE0m83E\nxMR4gze6W4cLFkKELoLIMJncW0NZt6WOZ99u4Ae3haKSB6YQtcyZRvjmHUg1Z1Bionq+4CKQJZkb\njRmMCIvnzfK3qYpayzURj6BThffJeIK2WcqSX97U52NcKtXV1UiSRGpqqrcsLS2N9PR0Dh482Ksb\nfyBxOBzY7XZqamrIz89n6dKlrF69mieeeMKva3ty3V0sNpuN0tJSduzY4S1LT0/HYrH4iFBubi5W\nq5Xdu3czZ84cHA4HRqPRpy+TyYTHM/CXDYQIXSQjomWW3xTKb19t4JWPmsi6RhdskzqlJWkqACF7\nvqHp+tQeWl8ayaNMGEINbNlbxBtN/8PimEeJ1ozq0zGHKxqt+pJnKf1B+43xwhsktN1wOxMhg8HQ\naV/V1dUYDIYe64uKijoVifbQ6CVLlpCfn++dyWzYsIHFixdjNpt9Iua6orCwsNP30xlGo5Hy8vIO\n5aWlpRQXF2Oz2YC2MO9Vq1ZRVlaG0Whk27ZtXrejy+Vi7ty5PtebTCYSEhKwWCxs2LCBkpISH5fd\nYEKI0CUwbbyarGu1bHm/iTGxKhbM1ATbpA64x8Tjjo9Bs7vvRQhgkn4i2XMe4KWvtvDqqV+QZshj\nXOjMPh9XMDBJSEhAUZRe7U9pv8GfOXPGx73kdDpJTEz0ikRn9Uajscdx7HZ7BxdZZmYmP/nJTygp\nKWH58uXdXr9582by8vL8ei9ZWVmd7i1KS0sjLS2NvLy8DjMf6Oii64zc3FxWr17Nhg0bcLn6dn2w\nLxm4q+qDhGtna7hqloaX/9bIPw8OwHThkkSzaRoh9v6LXovRxfDg5fczPmocb53+DV+dfb/fxhYM\nPFJSUjoEG0iSRHJycqft9Xo9RqOxw/qOJEkkJSV1W++v0HUW3GA0Gv2a4VRUVHQ5G7uQ1atXU11d\nzbZt2/xqn5qa2mlghsPh6FCWk5NDTU0NBQUFLF682K/+ByJChC4RSZJYdp2WmRPVWLc3sO/IwBOi\n5vmXo/niANKZ/juAVKvWctdMM3NHp/BhzV/4pKZInMQ9TFm9ejXFxcXe11arlaysLO+Mxm63d5hZ\n5Ofn+0SrWa1W1qxZ43d9d5hMJux2OwcPHvQpt9ls3nWXzmxqx+l09iqsuqSkhPz8fDZt2uRT7nK5\nOgiOyWQiKyvLp63L5cJut3foV6/Xs3jxYoqKigZcxFtvEO64AKBSSTx4ow7Lmw0882YDD98WSuKY\ngfOrbVqQhOT2EPLpFzSlz++3cWVJ5qbEG4jRxfCOYwet2iPM1/0AtRzSbzYIgk9aWhoOh4OCggKg\n7cFt8+bN3vqKigpKS0t9rlm+fDnr169n06ZN3mN1zg/r7qm+J7Zu3covfvELJEnyRtWdL2Kd2dRO\nTEyM3zMhaBOW/fv3U1BQwLJlyzAYDOj1euLi4tixYwfr16/3ab9582ZWrlzJ+vXr0ev1SJLUpYsw\nLy+vx022Ax2pq5j7oYokSclA2YslzzE9aVpA+25uUfj96w0cOe3mke+GMXHkwAlTjr3jx7QkTeHM\nf/vnyw40X1fto2jfXxkVPpJFUf+OThURFDsGKgf2fsFj99xMWVlZl24qgaAvsdlspKSkdPsZbG8D\npCiKYgvEuMIdF0BCNBJ5t4Yy2iDz21fr+bZy4MTkN1+VhPajPRCkh47psdO4f/Y91DQ6eaNKpA4X\nCARtCBEKMLoQiR/cFkZslMxvX23AVTcw1kGarpqD6lQ16gNHgmbD2MixPJh0H5Ik89rpX3CyqWPo\nqkAgGF4IEeoDwrQSK74TiiTBH95pxO0JvsuzOXk6ii6EkH/sCaodbZFz9xEfFsebp39JRX1AZvQC\ngWCQIkSoj9CHyzxwg47yY27e+KQ52OaANoSmubPQfuTfAY99SagmlH+77HtMN0ylpNLC52dKujwP\nTCAQDG2ECPUhU8aquW2Blh1lzXzhCH7odvOCJEJsXyPVNwbbFNSymiXTvstV467kE+dW/lGzGY8I\n4RYIhh1ChPqYxSYNSUY1f3inLWoumDRdbUJqaSXk48+Dakc7kiSxeNIibkm8ia9q3+Oj+l/TqrQE\n2yyBQNCPDJzNLEMUSZK4N13H06/Us+H1Bp5YGkZMZHC03z1hFC2J49C9u5OmtJ6Poe8vUkcnExkS\nQeG+V6hrWcPCqMfQysPz8NO9e0VeJkFwCNZnT4hQP6ALaQvdXr+1nt+91sC/Z4URpg3OqdtNaVcQ\n9tLb0NIKmp7//Iqi0FLjpv5oM+46N55m5dyXB3dT23dPs4KnSUFxK6jDZFThKtThMqowue11hMpb\nrtGrkNUd3/u02Knce9ndvPTVFt5qXUV69OOEq3uX7GswExVtQBcaxj333BNsUwTDmLCwMOLi4vp1\nTLFZtR85Ue1mfWE94+NVPHRbKGpV/wuReq+DuDsLqH7m/9B85eU+dZ4WhcbjzdQfbqb+aBMt+2XO\nnqiltdF3PUsVokKlkZE1KlRqGZWm7TWyRGtjKy2NrbQ2tNLa1HEdTJIkQmN0hMWGoUpoRTdKg3ak\nBt0oDRq9isqGKl7850tISNxg+DHRmpF9+vsYSJw6/i1nnP7lwxkoaI9/S8Ka/8PR3MeonZHUY/vK\nk98C8G1TGJ85qvmP2y/j0BefkJSUxMiRw+dvPVCJi4tjwoQJXdb3xWbVITMTkiQpAcgCKoAEwKoo\nyoA6WnaUQUXuLaH89q8N/Lm0kXvTdX2Wl6QrWqdPwj0mHl3pTpqvvJzGE81U/v0stbZW6irrUc6F\nk4cZQokYEcHEK8YTOSqCyBERhERokFWy3zYrHqVNkBpbaGlspaWhhUZXE3VV9dRV1XP2s2ZO1rjg\n3HOQWqsiLDaMWydnUqbZRUncb7hm0nJG6hL66tcxoBgxeiwjRo8Nthm9Y/plTNn8HKOPHeXIkp5n\ncWERbSdlzBs3idoPD/Lk82+gP/4PHnroISIjI/vaWsEAZMiIELBVUZRUAEmS9MBWICO4JnVkylg1\n/5au4w/vNGKIauY787X9a4AkUXf9fJwfVXHkZ1XUHHahCVUzYlo841LGEjUygogR4agvMa0zgCRL\nhIRpCAnrOsWFx+2hvrqB+nPCVHu6DteXZ5hcPZvJzKZet4/qscfRTBxLyLgoNGMjkUOG0sd2kCNJ\nuOZfi+Hd7RxZ8QSo/DuqSpYkYk/v5MjON0l+/GkhQMOYIfHfLEmSCe/zNCiK4pIkKVWSpEmKohwM\nnmWdkzpVQ81ZD3/9RzMj9DJXzOifPET1R5o4/cFZahxX0RqvYGioY/YdMxkxPQ6VOjjn3MkqmYj4\ncCLifQMRmutbqDpaxQcH3ufMcSdxHzZDMyBB1KhIlMkxaCfHEjIuEkkWQZ7BxDX/WkZu+zPh+76k\nbmbPLjmA97a/wgdvbOZHqzex42sXxV8cJ3326D62VDAQGRIiBKQCFzrTqwEjcLDfrfGDxckhnKjx\n8PJ7jUwapWJkTN/cSD2tClX/OEvNjibOHD9LSEQI4+eOY+JffkrouAW4L7umT8a9VELCNIyeOoqs\nKZm8XPMn3qx9mdvdZqbWzKD6YA1Vnx3j7N8Po9aqiTXG0DTFgC7RgCqqn2eWAuqmXUZLtAH9x+/7\nJUI7Pyjho5I3+e/f/JGwiEiOnlV48tUvuXxCNCP1of1gsWAgMVREqLMwKmcX5QCc+PYk4ZHBDQNe\nMEHhi88bePpPsPzGUDSdRI1dCi1nWzlidVF3up6YSdGMujaemIkxyCqJk5MnIL33Nq1pt0A/r0v1\nlgXKdTSebeIv9c+xIPw6FqYuJjJlLGdP1+I85OLg4YOc3fUlAGGxoUjGaHRT41DHihtaf6HMmE34\n+zs4lv4d6ObjVPraVv7x7lv8+Oe/weWswuWswhTjYd83x/n+U6/y639L6fd1UoH/HDkS+LMnh4oI\nOYELE3xEnyvvlOd+9UfCI8N8yq5Jv4prbrgq8NZ1QYhawnydlk1vNfLqx01kXq0N2D+gu8nDkY0u\nWhpbSDLPInKEb+oET8qVqD8sRTpUjjJpckDG7CskSSIt6gbCVOGUnn2bBk8dN0XdTtTISKJGRjJh\n3jhaGlpxHXVRfchJ1eenqNt1nIiR4Ugz4gidGoukGyof9YHJ2dkpxPy9FE3lSVriO49yUxQ4VP4N\n9z6c7w1QANCGyCycOZI37Mf44OtTXDdDRMkNBF5//XXeeOMNn7KzZ88GfJwhEaJ9bk1oo6Ioc88r\nqwaSL1wTCmaIdlfYD7Sw6a1GvjM/hBvnXro7ydPiwbHaydkTtcy910TkyE5y97hb0Xz/NjzXZOBe\n/uglj9lffFL3ES/X/BFTaCr3GB5AJXUUF4/bw+lvqvh2z3GqDlQjyRIjpsXRYhqF1hiDJIsn7UCj\nqj3L7Ltv5PBD+VRn3NZlu2NH2tJUjxnfMeJx264jHKqq46+PXUdkaP+skwp6h8gn1AWKotg5z/Um\nSVI0UD4QgxI6wzRZw83zQnj9k2b2VFzasTWKR+Hbp+twHT2DadnszgUIQKXGc/Vi5A+LwT1w8h71\nxPzwq7jfkMOeBhvPVj1Ds9LxcFhZJTNyRjzJd17OtY9eyeSFCZw9VUvVnz+n6ulPcJVW0FpVHwTr\nhy7uiEjqjdOI/OLi70uLLxtFS6vCys3BP2RX0H8MCRE6x1JJkh6XJCkTyAeWBtug3nDTvBDmJKp5\nfkfjRSfDUxSFU5ZmTn59msuXzCJmQvcnDniuzUCqqUL6cnClU5gTlkJ23MN807SXjZX/S6On6wNZ\ntRFaJl05gQV587ji+ynET4mlcdcxTv52J40vfk7j/ipxgneAqL08mYjPyy46cWJUqIbrZsTzmaOa\n8pOBd/sIBiZDRoQURdmtKMp6RVGKFEVZOVhmQe3I586Yi9fLWN5s4GxD70+UPv66kyOffcvMW6Yx\nYlrPR28oU2ehjBqL/MGOizE5qMzQzWJF3KMcaT7E70//kjpPXbftJUlCPzaKmTdP47rHFjDrtuk0\n1zVT9ZcvcP52J7WfHsXTyQkPAv+pnZlESHUlmtMnLrqP1AQD+jAN//PKFwG0TDCQGTIiNBTQaiRy\nbwmlqQU2bW+k1e3/E+Xp985w7NUaJi9MYJxpjH8XSRKeazOQP/4bNDddpNXBI1E7hYfjf0ylu5L/\nPbUOl7vLOBQfVBoVY5NGM395KvPuTyZqVCRn3inn1C8/xvnWfuGqu0jqp80CIPzrLy+6D5Usc+20\neL4+fpavj50JlGmCAYwQoQFGbJRM9s06HCfcbP2gyS9XUU1ZHYderGTC3LEkXDWxV+O5r81AqqtF\nKvv4Yk0OKuNDJvJI/BPUK/X85tQ6qlor/b5WkiSix+u5PHMW1zwynwnzxtH85SlO/nYnTX/+nMaK\nGuGq6wWt0QaaRo9l947Xef/tv150P7PHRWMID+HnfxWzoeGAEKEByOQxau68XsuHX7bw9y+6D1Q4\nu68Bh+U0o2aMYNoNU3of4j0+AY9xKqr337kEi4PLKM1ofhT/EwB+c3odJ1uO97oPXZSOKQuNXPuj\nK5l123SazjZT9ac91D1rp/EbsW7kLzsiY/j930u4PHXBRfchyxLXTh/BgZO1fHnEv9mtYPAiRGiA\nsmBWCNcnadj6QRP7v+16reKw1Yl+XBSX3T7jovcYea69Aemzj6B28Lo/YtVxPDLiCXRSKL85vZ6j\nzRe3qU6lPueqy04l+a7LkWSJqpe+oHaTjaaDNQG2emhh+/h9fvn15/wpKgqD/tLScMwaG0V8pJZV\nr/0zQNYJBipChAYwS67WMnmsCuv2RipdHQMVmmtaaahpYMLcccjqi/9TehbeCB438t/euhRzg45e\nFc0j8Y8TozLwu9NPcbj54EX3JUkScZNjmXufiZR75gBQ+cIemv78OS0nagNk8dBh96d/5/nf/IIn\nV65iFBLabw9fUn+SJHHl5DjKT9Vy8LT4fQ9lhAgNYFSyxIM3hhIaApY3G2hs9nUJ1VW0BRPox0Vd\n2h6Hl94AACAASURBVEDRsXjmX4/q7VcuOrx2oBCuiuAH8Y8RrxnJ707/EkdT+SX1J0kSsQkxXPH9\nFC7PnEV9dQOnLJ/h2baXVmdDgKwe3OzZ+RHP/up/+K+nnydiTioAusOOS+535tgowrUqtn56aYIm\nGNgIERrgRIRK5N4aStUZD38sbsRznkjUVTSijdKiiwzAKQs3fBfp6EGkr/Zccl/BJkwO4wdxjzJG\nM44Nlb/mQNM3l9ynJEmMmjmCBXnzmHHzVKocNZz67U5cJRV4modvaPcXn33Mpqf+m/96+gViR4zG\nHRFFiyGO0MMVl9y3WiWTPMnAa7ZvOdtwaZu4BQMXIUKDgDGxKu7PCOXzilbe2vmvEwKa9krox1zi\nLOgcyuwUlDHjkd95JSD9BRudHEpe3I+YEDIJS+Vv2Ne4NyD9yiqZ8Sljufqh+SRcNZH6T49S/dtd\n1H95clgGLzirT/OfT79A3Mh/pWFomJAQkJkQQMqkGNwehdds3wakP8HAQ4jQIOFyo5pb54ewfWcz\ntgMtKG4F17Ez6McGRoSQZdwZ30X+6F04MzQikrSylpzYhzGGTMZa+Vv2NgZukVsdomLydQlctWIe\n+jGR1BTtpeGFPbScHF7rF9dk3Eb8KN99aY0BFKEInYaZY6P444cVuD3DT+SHA0KEBhE3pIaQMkXN\nn4obcXzZiKfFEzgRAjxptwAgv/tmwPoMNiGyluVxDzFFNx1r5e/4siGw7sbQ6FDmmGeT/L0kmmqb\nOWX5DOdb+/EMY/dR4wQj2uNHkQK0AXqeMRZXfQsf7jsVkP4EAwshQoMISZK4J03HyBiZ0u21bVlG\nRwcwLXJUNJ4FC1G989dBH6BwPhpJw4OxK5ilm82zVRvY0xD4s/LiEg0syJ3L1LREGnefoPK3O6mz\nHRuWLrrGCQlIHg+6S4yQa2dMTCjjDKH8vuTS1/YEAw8hQoOMkHNH+0S7WqgL10CAE+F5brwD6dgR\npC/KAtpvsFFLau6PzSYpNJnnqzZiq98Z8DFklcykKydw9Q+uIC7RgPP1b6h71k7L6e7PtRtqNJ5L\n06A7dOnBCe3MM8ZyqLKe/SfEwaZDDSFCg5CYSJnpIRpOq2Usu78N6NO2MnMOnvEJqLYXBqzPgYJK\nUvNvhgdJCZvHH6ufpawPhAhAG6ll9ndnMvc+E62NrVRuLKPus8D+nQYy7ohIWmJi0R4LXBbOaaOj\niNSpebUs8Jk9BcFFiNAgJWy0h8TQED781sWWrwPoK5ckPLdkIX36AZzq/fE3Ax2VpOJ7MfeTGnYF\nf6p+Flv9rj4bK2ZCNPOzUxmTNArnm/tp3fxP3PWDb61o/1d7+M3/Z+8sw6O6ujZ873NGMpZJMnEX\n3L1oBUqVGrQFWqpQe0tdqL11b7+6lwoVtLRFqlAvpbgHSYhDiLuOnO9HQt/SIslkZpIJ574uLmBy\nZu8VO8/Ze6/1rEfvbNV7GsMi0BUXeCwGWRL0jrHy9eb9OJytd5hX6bioIuSnaK0yAQ4XU3tG8GV6\nEav3V3hsbNcpZ4HBiPzNYo+N2ZGQhHSIEG2qXe+1uWStTK+zutP/oj6U5ZRT9tZ6v7L/Sd+5lRcf\nup1JV/ynVe9rDI1AW+TZRILesVZqG52szyz16Lgq7YsqQn6KJlCmsaaRc1JCGRlj5e3NeWRXHrm5\nW6sIMOAafy7S90uhwUNjdjAOCtFA41A+Kp3tVSECiOgRxohrh2IMMVA8ZwuVP2agdPAn+vSd23jh\nwdu477l3iElIbtV77aFhaEs8K0KR1gBsZh3v/Zzu0XFV2pfjVoQU/Ht/XmuVsdfaUVwK1/aLIdKk\n54V12VR7qHrfedaFUFuN9PO3HhmvIyIJiWnBVzHQMISPSmezuda7yRgBgQEMmTaALicnUf17LjUf\nbMZR1jGtfzJ27+CFB2/l3mffIjYxpdXvt4dGoCsu9GiWpRCC3rFWdu+vot7uPy3pVY7OcStCKzN/\nxO70v/35g2gDNQA01trRayTuGBpPrd3FKxtyPVPUFxGNMnQM0vKFnSpd+59IQuLSkKsYYBjCnNJ3\nvS5EQhIkj0lk6JUDaahppOTt9dSlFnl1ztaSuSeV5x+4mXueeZO4pK5ujdEYGo5cV4tc49ni3T4x\nVhqdLlbt6VhfMxX3OW5FKKcqjw+2fkRlg3+2L9BaZQAaq5tsfMKMOm4ZEseOkhrm7/LMgbBzwsVI\nORmdLl37n8hCZlrIVQwwDGZO6btsqdvk9TmDYq2MuGYotuQQShftoGLFXhRX+2/PZe/dzXP338Ss\np14nPrmb2+PYQ8MBPL4lF2LWEx1k4KPfPJf+rdK+HLcidNug26lx1PLu5vfJq/I/XyptYJMINdT8\nz0uuT6iZab0iWb63mFV5bbfeUfoOwpWQgrxsQZvH6ug0CdHVzXVEb7PVB0KkDdDQb1Jvup2aQs3q\nPOo/3orzb9/P9sAcaGXWU6+T0KVHm8ZpDI0AQOfh5ASA7lEWMgqraVC35DoFx60IxVhiefWUNwgO\nCObDrR+ztdC/WglrAg9dCR3kjCQbY2KDeGfLPjLa2mpACFwTLkas+x3y89o2lh8gC5nLQqbTzzCQ\nD0reZlvdZq/PKYQgcUQ8g6f1p7qohrJ3NtDYjgWZtrDINgsQgD3EhiJJaD2Ypn2QbpEW7E6FDVlq\nllxn4LgVIYDggGBeOvk1+ob15os9S/kuYyUupf23RFqCpBGYQo2U5Ry64hFCMKNfNHGBATy3Npvi\n2rY9WbtOOh0Cg5CXzGvTOP6CLGQuD5lOP8MAPih5x6Omp0cjJDGY4TOGoDfpKP1gM/Xpfn6DlTU4\nLFY0XjDDDbXosQRo2KCmancKjmsRAtDKWh4Z8SQz+l7Lmv1r+XTHfOrsHTNj6Z9Yhuko2lOC6x9n\nCTpZ4s5hCWglwbNrs6lty7aFPqDpbOiH5VBe0saI/YMmZ4UZ9AjoxXvFb5BWv9sn8wYEBjDk8oEE\nJwRROncbNRv9u1jYZTQi19Z6fFwhBHE2I7/s9PwqS8X3HPciBE0/1OemnMejo54gv/oA7255n4Ka\nju/YGzzIiL3OTnnOvwtVg/QaZp2QQGm9nZfW5+BoQ8ac66xJIMvIS+a3JVy/QiM0XGW7jiR9F94p\nea3NHVpbPK9OZsDFfYgZFEX5st1U/pTpt3Y/ToMJqdY7vnlxIUbyy+vVc6FOgCpCf6N/WH9ePuU1\ndJKO97fOIb3MNzcedzEm6tEH6inYdfh01RhLALcNiSe1pIZPUtvwVG0OxHXmJKRvFkO1f2YTuoNW\naJlh+w+x2jjeLn6F3MZsn8wrSRI9z+xG17HJVP2ajfLlbo9nztnt3k+AcBpNyHXeEyGnS2Hn/uPn\n57GzoorQP4g0RfLK2DdIDIxn7o4FrMvvuOnJQgisQ/UU7io64tNy71AzV/SJ4rvMUn7Jdd8uxnne\nFHA6kZYvcnsMf0Qv6bku9CbCNBG8WfwS+fb9PplXCEHSqAT6nt+LA9sLcCxMRXF4Rojyc7O484rz\nqPPSKuUgTpMZ2UtzhAcGoNNIbM3xHwsklcOjitBhMGgMPHPiC0xIPoev937LdxkrOmzCQvBQMw1V\njZTnHdk77tSEEE6JD+a9rfvZ627GXJAN1/hzm9K1vXzz6mgESAauD7uFIDmY14teoNDuu7OIqL4R\nDLi4D8XppTTM24arsW3bTwfysnn8junccO8TGIwmD0V5eFxe3I6TJEFMsIEV2/z73ExFFaEjIguZ\na/pdx7V9r2PN/nUs3LmYRmf71nAcDnOKHr1ZR8FRqu6FEFzZJ4r4wABeWJdNRYN71j7OC6ZBfS3S\nt1+4G67fYpJM3BB6K0bJyOvFL1DiKPbZ3GFdQxl0ST/Kcyup+2Qrrnr3vn8H9uXw2O1Xc9MDz9Kj\n7yAPR/lvnEaj11ZC0LQll1tah0tt++3XqCJ0DCaknMsDwx8kozyTD7d9TFVDx2qqJSSBdVgABTuP\nvCUHTRlztw2Jx+FSeHmDm4kKYRG4TjkLecncTmtsejQsciA3ht2OBg2vF71AudN3W0EhicF/1RLV\nfLQFZytT7wvz83j8tqu58b6n6dFvsJeiPBSn0eRdEbIZqbc7yS4+vlbmnQ1VhFrA0MhhPHfi/1Hd\nWMPsLR9woLpjpYYGDzHRUNVAadbRazJsBi23DYlnT2ktc7bnu5V15Zx0OVSWI32/xN1w/RqrHMSN\nYbfjxMnrRS9Q5fTdwXhQrJUhlw2grqKe6jlbcNW1zPuwMH8fj916FTfc+wS9Bgz1cpT/w2UwItV5\nPkX7IDHBBgSwNdfztUgqvkMVoRaSHJTCK6e8hklr4v2tc9hd0nH63Zu7BWCJMJPxe9Yxr+1hM3F1\n32hWZpfyTYYbdT9RsbhOPgN58UfH5WoIIERjY2bYHdS56nir+BXqXb6rKwuMtDD0soE0VDVS++k2\nXMdwTa+uquCxW6/kulmP0XvgCT6KshkhIbzoh6fTyASbdGQUdqzdCZXWoYpQK7AZQnn5lNdJCU5m\n/s5FrN3vva6crUEIQfiFJsqyyv/loHA4xiaEcE5KKJ+kHmD9gdY/yTsnT29aDXXSpnctIUwTzg2h\nt1DiKOLdktexK75zZDeHmxh0ST+qi2pomL/jqFlzJnMg9//fbPoMGu6z+HxJmEXPur3HRxF1Z6XF\nIiSEGOvNQPyFAE0AT49+nvNSLuCbjO9ZmfVThygmDBpgxBJhZu+vWS26fkrPCIZGBfLaxtzWe8xF\nxuAaNwF58cfgxe2Wjk6MLo5rQmeS3ZDJnJJ3cSq+K5y0RgcycHJfynMqcHyWesQ6IiEEkbEJPovr\nUBQUIbw6g82ip6iqwatzqHiX1qyELvJaFH6GJCSm953B1X1msCrvD5akLcPpat/KbSEJwieZKM0s\na9FqSBKCGwfGEmdxz2POefFVUFuD9NXxVTf0T1L0XbnKdh076reyoOxjnz6QhCQG0//C3hSllSCW\np3WIh6FDUBTwsgiFWfRU1TuocTNjUKX9aY0IefenyQ85v8sF3DnkbrYV7WBe6kIaHO37RBY00Ig5\n3NTi1ZBOlrij2WPuuXU51LemGDIsEtdp5yF/8Sl4uHGZv9Hb0I9Lgq9kTe0fLK3w7RZlWLdQek/o\nwb5N+VT/nuPTuY+JgtfvGqEWPQBZxcf3z6A/0xoR6mCPWR2DE2NP4pGRj5NXtY852z+hurH9fhmE\nJIg4uBpqYcZQkF7DXcMSKKhp5O0tea16mnZeeAU0NiAtO3485Y7EUNNwJlon82P19/xQ5duW6NH9\nI0k+MZHKHzOp3daRMje9f8uwmZtEKLNITdP2VzyemCCECBRCDBBCTBRCDPD0+B2R/mH9efbE56lq\nrOb9rXMoqWs/i/mgQSbMYSYyWrgaAogLDOA/A2P4c38lS9NbUYRpC8N1xgVNbR6qjuzYcLxwkmUc\np1vOZmnF56yu+c2nc4f0sxLWy0b5kl00tGA71nd4dymk00hYjVqyitSVkL/iMRESQliFEN8D5cBG\n4DNgoxCiRAgxw1PzdFSSrMm8dPKryELm/S1z2NdO3VoPng2VZJRRnttyYRgWZeWCrmEs2FXApoKW\np7w6J10OTifyl3PdCbfTcWbguYwyncSCsk/YXOsb38Hi0mKuvPlS6OogKM5Kxfzt2Es6QMKID86E\nAELNeta05uFJpUPhyZXQPcB1iqJIiqJIwHggpflPl+NBiMKN4bx0ymvYDCHM2fYpaaXp7RJH8GAT\n5nATu1ekt2p77cLu4QyMsPDqxlxyK1tYAxQUgmvCxUjLFkCpeiMQQnBh0FQGGAbzUel77K7f6dX5\nSspKuPrWy7nl2jsY2G8QAy7sg86oo2b+9jb7zLUV4XCgSN6vArGZdZS2c1t0FffxZGLCekVRMg/+\nR1GUH4AkRVHKFUW5pwXv7xQE6gJ54eRXSA5KYl7qQrYVbvd5DEISxFwZSMW+SvZtarnBoyQEMwfF\nEm7U8ezabMrqW1b74px4Geh0yPNnuxtyp0ISEtNCrqabvjvvlbzJvsZcr8xTVl7G1bdczk3Tb2Hc\nmFMB0Bq09L+oD/UVDbB8T7tmzMm1NThNFq/PYw7QUuOmH6JK+9MaETrW3sKQv/9HCJFI0yoIIUQg\ncNxUlOllPU+Pfp7+4X35fM8S1rdDOwhLdwPR/SPZ88NeGlrxlGjQyNw1LAGnovDc2uyWZcyZLTgv\nvhppxVJEdsfuweQrDjbFC9dE8Hbxq5Q6PPvjX1ZexlW3TOPGq29i/EmnH/Ixc5iJXmd3J39bAbXt\n2J1VrqnCaTJ7fR6TXkO93aU2uPNTWixCiqK8e4xLFgohSoUQaUKINGCRoijvCiGsNJ0ThbQlUH9D\nlmQeHvEE5ySfy1d7v+X3vD98HkPYVXqEEOxZ0bptQZtBy13DEthf3cjrG3NxteBp2nXmJIiIRn7v\npaazABX0UgDXhs5EFjJvF79Krcsz5zTlFWVcfetlXHfFfzjt5DMOe01U3whiB0dT+U06jfntY2sj\n11T7RITMehmAMnVLzi85pggJIcY2r2SOiqIom4Akms6GrlcUZWjz6xXN50TH3V6NJCRm9L2WKd2n\n8kPWT/zgY3cFrUUmZkoQ+dsKKMlsXcZektXAzYPj2FBQxaepB1owmRbHVTcjbVmHWL/KzYg7H4Gy\nletDb6bSWc57JW/g8IC9z50P38Y1067nzLFnH/W67qd1wRRmpHpRqtvtH9qCXF2F0+z97ThTgAZA\nPRfyU1qyEhJAefMK500hxMQjXdgsOIubz4NUaDqovqTnNK7uM4Pf8/7g64zvfCpEtlFmghOCSP1q\nD85WblcMirBwZZ8ovs4o4fvMY28nKcPG4Oo/FM37L4Pdd15qHZ0IbRQzQm8kqyGDT0vntLlB4itP\nvsFZp0445nWyRqb/hX2w19pxfrnL5+dDGh+thEz6JhEqqVbte/yRY4pQs6A8B2TSlHotms97VFrB\n+V0uYOaAm1mfv4Ev9yz1mc2PEIKYa800VDa02Enh75yWZOOsZBsfbs8/duq2EDivvgUO7EP6+jP3\nAu6kpOi7Mi3kajbVrWN55ZdtGstoMLb82mADfc7rQeHuYmo3+KY1+UGatuN8sBLSNa+EqtWVkD/S\nku24icAaRVFOUxTl3eaVTtY/rrlLCHGnt4LsLJyWeDp3D72H7cWpLNr1OQ6Xb7ZIAiJ1RJ5rJXt1\nLpUHWn8+cGmvyL9St/dVHT11W0ns0mTns+A9qOxIRZPtz0DjEM6zXsgPVd/ye/XPPps3vHtY0/nQ\nd3ux+7ABnFTrm5WQJAmMOlldCfkpLdmOG68oyudHu0BRlOcURXleCPF0S86PjmdGx4zhv8MfYm95\nBnN3LPBZy/DIM4IwhRpJXb4bl7N120EHU7dtBi3Pr8uh6hg9bJyXXAuAPPcdt+PtrJxiGc9J5nF8\nVj6PbXWbfTZv9/FdCLAGUP/5riM6bnsUp9NnKdrQtCWnroT8E49WkjXXA13ryTFbgxBiUnvN3RoG\nRwzh0ZFPsK96P5/umO8T41NJI4i7zkpVQTV7f8k89hv+gUEjc+fQBGrtTv5vbQ6NRxMyazDOydOR\nvvsSkbG7DVF3Ts63XkRfwwDmlL5LVkOGT+aUtTJ9z+9JZX4VNWu9vy2nqShDKAr2oGCvzwWg18rU\nqrVCfklLRKi1qdWLj5a84A2EEJOEENcAi/xlJdYntA9PjHqKgppCnwmRKSmA6AuCyVyV0+psOYAI\nk467hiWQWVHHG5vyjpq67Tr7IpTYROS3ngdfPHn7EZKQuCxkOrHaeN4peY0iR+G/rqmprWF3+i6P\nzmuNDiR2cDTVP2Xi9HIPHl1hU31SY0SUV+c5iCyBQ/0580talB3Xmht7s2tCsvshtZ7mc6p38TOn\n7+4hPXhi9FMU1hbxyY551Du83y478gwrIUnBbP9yJ41upLR2CTZy0+A41uZXHj11W6PBee0dSLu3\nIf30dRsi7pzohI5rbDdikky8W/zaITVEtXW1XHP7lWzftc3j83Y9JRlJI8H33l2B6QuaRSjcNyIk\nCYHD6Ve//irNtESE3gaOVaj6T2xuxOIJ/M4aqFtwd54c/QzFtSU+ESIhCWJvtOByKWxf5l7a7pDI\nQK5oTt3+NuPIqdtK38E4TzwNec7rUN0+BZMdGZNs5hrbTCqdlX91Zj0oQOeefj6TJni+j6TWoKXb\nqSkc2F5AQ2aZx8c/iK4wH4fZgsto8tocf0eWBA6XKkL+iOZYFyiK8oMQ4hkhxHRFUd5r4bg+XQm5\nw2233YbVaj3ktalTpzJ16lSfx9I1uCtPjXmGe3+fxSfb5zGtz1QCNAFem08XrCFhho30Vw6Qs24f\nCcNiWz3G6Uk2imob+WhHPjaDlqFRh18sO6+8Cek/k5HnvoPz2jvaGnqnI1wbwVW2a3mr+BUWFcxl\n+aNfcfb4c5lywaVemzO6XyT7NuVT93UauuuHIGTPm4zqCg/4bBUETSuh/WUdwDm8EzFv3jzmzZt3\nyGsVFZ5v2XJMEWrmYmC9ECJFUZT7jnZhs01Pm1ckzWc8Kfx7i000v7ZCUZQf3R3/xRdfZNCgQW2I\n0LOkBHXhqdHPct/vs/h4+1ym9ZmKQWPw2nxB/Y2Ejw9kz8p0guOtBEa2Povpkl6RFNfZeXVjLv8d\nmUTX4MPUr9jCcE6ZjvzR67jGn4uS1NUD0Xcuugf0YoJhIvfOuouzTz6XSyZO8+p8Qgh6ntmNP99d\nj/xnHpZR8R6fQ1eY77PzIGhaCdnt6pmQJzncQ/nGjRsZPHiwR+dp0SOQoigZwKnA9c3OCdMPd07U\nXMS6EnirrYE11yTdoyjKvf/4c/C1wwmQX6/HU4JSeGr0M5TVl/PxtrnU2eu8Ol/sJBvmUBNbP0/F\n4YbtvyQE/xkYS7LVwHNrszlQc/jDbtc5kyE6Hvmd51VfucNQX1/HnIc+ZNTo0ZSfXERavfczCi0R\nZuKGxlDzSzaOCs9vAesK8n2+EmqJx6FKx6M1BqYbadpmy6LpjKisWZAWNNv5fAfsBVa2ZYXSRvzu\nTOifJAel8PSYZylvqOCj7Z9Sa/feFoOkFSTcHER9ZT27v0tzawydLHHHsHgsWpln1mQfvoZIo8Fx\n3Z1IqVuQfvqmjVF3Pj5eNIexo8fx/BUv00XfnfdL3zpsxpynSTkpCVkvw3cedj5XFHRFB2gMj/Ts\nuEdBkqCV5W8qHYRWbQY39wYaD5wGfE5TAsJFwGSaBOA0RVHu9XiUx0AIMU4IcRdNK6F7hRBjfR2D\nJ0myJvPMmOepbKji4+3eXREFROqIn2Zj3+Z88ja5Vz9i0Wm4+4TEphqidYevIVL6DcE5ZjzyB69A\nefu1P++IzJh2HVdOmY4sNFxluxazZObd4teo85Dr9pHQBmjofmoXCnYWeTRJQVNWgtTY6POVUHv2\nTlJxH7dOJBVFWakoykWKooQ0O2SHNNv6tItxqaIoPzS7NshH2arzKxKtiTxz4nNUNFR6vY4odLSF\n2MHR7Px6DyVu3owiTDruHJZAZnkdr27MxXmYTCXnjNsA0Lz1nLot9zfE31pgGyXTvzLmvElkn3AC\noyw4fs7y2E3ckNW0sqpL8F1+ksOpoPFCgoWK91G/ax2YhMBEnhj9FMV1JcxLXYjd6T1n6vBr9YQk\nBrHls+3UuOkv1jXYyK1D4tlYUMX72/b/+6YWFILj+ruQVv+E9NtKD0TdOTmYMbe7YSdLKrxrBCuE\nIOWkRMpzKmjM8ozXnyEzDafBSGNEtEfGawkNDifxtpYbu6p0HDwuQkKIRCHEAk+Pe7zSJagrj458\nnP3V+SzY+ZnXTE8ljSDutkD0Zh0b52+jsdY9H66BERau7R/DjzllfLb73+cayqhxOEeNQ37nOSg7\nbprttpruAb2YGDSZX6p/4I/q37w6V2gXG4HRFpw/Z3lkPEPGHuqSujYd1PiIBofrr5YOKv6Fx35K\nmsXnTZqSEy701Lgq0NPWiwdHPExWRTaLd33htTYQGqNM8p02HA0ONi/ajqslrb0Pw0lxwUztGcHn\naUV8n/VvoXFefxdIMpo3nla35Y7CGPMpjDKdxGflc8ls8F7bdCEEyaMSKMupoDGv7XUghsw06pK6\neCCyltOoipDf0mYREkIE/k18rgM2tTkqlX/RP2wA953wAHvK0lmStqzNjdGOhD5MS8otYVTuq2LH\nV7vdPic4JyWUM5NsfLgtnzX7/3FjCwzCccMspLW/If38rQei7vjYHXYWL1/U6q/nxKDJxOsS+aDk\nbSqdlV6KDsK6h2K0GRG/57ZpHFFfT8C+HOqSunkospbRYHeqIuSnuC1CzeLzFFBGk/g8BwQDT3so\nNpV/MDRyGHcOuZvtRaksT//Ga9lA5i4BJEwPJX/rATJ/z3ZrDCEE03pHMiLGymub8kgtrj7k48qI\nk3GedDryuy9AiffTkdsTu8POrQ/MJC8/95AkhJagERqusl2HCxdzSt/xWqKCEILEEXEU7i7GXux+\nVp4hZy/C5aI22bdFyepKyH9xS4SaG9iVAbNoFp/mItIK/LxgtKMzOmYMtwy6jU0Fm/k2Y4XXhMh2\ngpnoc4NJ/zmTAzvcEwlJCG4YEEOPECPPr8shu+LQVHPnNbeDTo/m9ac67bacw+Hg9v/eTJekrtxy\nze1ujWGVg7jSdg0ZDel8VdG2rqxHI7pvJHqzDu0q91dDhow0FEmmPj7Jg5EdG/VMyH9plQgJIWYI\nIUqAZ/m3+Kj4iLHx47ih/42szV/Hj9k/e22eqHODiOwTwfalOynf5963WCNJ3D40nkiTjqfXZFP4\n94QHixXHjfcgbViNtHKZh6LuODgcDm5/6BYS45O4tY2+eV303TnXOokfqr9jS91GD0V4KJJGIv6E\nWPZvO+B2qwdDRhr1sQkoeu95H/4Th9OF06VgClBFyB9pkQg19+spAd6hyS1BFZ925syks7iq93R+\nz/uDVXmrvTKHEILomQYCI81smrfN7dRtg0Zm1gmJBGgknlidSWn9/1LNlaGjcY49G3n2S3BgPCll\newAAIABJREFUn6dCb3ecTid3PnwrcdFx3H79Xa3ehjscJ5tPZYBhMJ+WfkiB/ShtNNpA7KAYJFmi\nZkO+W+83pu2kLtm350F19qYtSrMqQn5JS1dCk4FFqOLTobig60Qmd5/Cyqwf2VywxStzSFqJxLuC\n0Zt1rP9kC7Vl7rk3WPUa7hueiMOl8OTqLCr/1gXTec3tEGhF8+LD4PT/7phOp5O7H72dqIho7vzP\nLI8IEDQ9FEwNvgKrbOX9krdocHm+gFkboCGiZxjObQWt3uqVaqoxZuymus9Aj8d1NMprmx5qooK8\nZ/ir4j1aamB6MfAMMNhfOpceL1zSYxqDIgeyNO0r9pS65/92LDRmmeR7bMhaiQ2fbqbeza2aMKOO\n+0ckUdXo4Ok1WdQ2P8FiNOG4/WHEnh1Iiz/yYOTtw56M3URFRHP3zHs9JkAHCZACmG67gVJnCQvK\nPvbKmWBU3whqS+uw72tdNp45dQvC5aKqr49FqLk5Y7QqQn5JawxMM5vtcFKEEBNVMeoYCCH47wmP\n0N3WjUW7Pie3Ms8r8+iCNKTMCkVxKmz4ZLPbxazRZj33DU+ksLaRZ9dkU99ci6T07I/rwiuQ572H\n8EJHUV/Ss2svj66A/kmkNpqpwZexoW4tv9f87PHxQxKC0Vv01G4taNX7LJvW0hgaQWNU6/tTtYXy\n2kaMOhmjmpjgl7Q6O05RlE2KonwODG1OVFDFqJ2RhczjI58m2hzFvNSFFNUWe2UefaiWLrPCsNc5\n2PDpFuz17m2dJVgNzDohkazKel5Yl4292fDUOXk6SrdeaJ57ACrVHd+jMcg4jBPNY/mifCFZDZ5t\n1S0kQVSfcBq3F6G01JpaUbCu+Y2KYaPAS+J7JMpr7QQZdT6dU8VzuF0n1GwaOhsYr4pR+6OTdTw9\n+nksOguf7JhHZYN3ChsDInV0uSuMuvJ6Ns3f6lYfImjymbtzWDy7Smt55aDhqUaD487HoaEezUuP\ngEv15j8a51kvbCpkLX2baqdn26dH9Y3EXmenPr1ljucBORnoC/OpGDbao3G0hLKaRnrHWo99oUqH\npM2OCYqiLP6bGN2JH7T27qyYdWaeGfM8AsEnO+Z5rQWEMU5PlzvCqSqoZvPCbTgd7glRn1Aztw6J\nY1NBFW9tzmtqShYWgeP2h5E2/IH0+ScejrxzoREarrRdi0Nx8HHp+x510bBEmLFEmNFuadmWnHXN\nbzgNRqr7ebbrZksor7UTHayeB/krHvOOaxaj52nqK6TupbQTNoONp0Y/S3VjDfN2es9525wcQJdb\nwynPrWDr4lRcbnYUGxQRyI0D41i1r4IPtuWjKArKoBE4L7oS+dO3ENu9UxPTWQiSg5kWcjW7Gnbw\nS7VnO6lE9Y2gaE8Jrvpj/wxZ166icuAJKFrfbos5XS6q6uxqZpwf43GbW0VRnlUUJcTT46q0nFhL\nLI+NeoID1QV8tvsLr/nMWbobSJkZTnF6CduX7MTl5vbZiBgr1/SPYWV2KXN3NqUGO6fOQOk1AM3z\nD3ZIt21FUXjl3RcpLvXO+Vtr6BnQm5PN41lW8Tm5je7ZLB2OyD4RuJwu6lKP/jlqykox7tlB5bBR\nHpu7pVTU2lGAmBC1jYO/ovYT6qR0C+7OfSc8QFppOl/v/dZr9j7WvkaSrwujILWI7V/udHtFdEp8\nMFf0jmL53mLm7ypAkWQcdzwKKGievR8cHad+SFEUHvu/h8gv2E9IUMd43jrHej5R2mg+Kp1No4fq\nhwIseoLirOjSjn4uFPzLdyiyhoqhvj8PKqxs+lwTbCafz63iGVQR6sQMjhjCzIE3s+HAJlbt846r\nAkDwEDPJ14dRsLOIrYt3uC1EZyTbuKx3JEvTm4Uo2IZj1pOI3duQ33/Zw1G7h6IoPP7iI9TU1fDE\nfc8g+bBnztHQCC2Xh8ygzFnKlx5shBeSGERZdvmRH2IUBduK5VQMPxFnoO+TA/LL6zDrNYQF6n0+\nt4pn6Bi/QSpeY3zCaUzuPoUfsn5iR1Gq1+YJHmImZWY4ReklTckKdveSFc5KDj1EiFw9+uG85nbk\nrxYh/fi1h6NuHYqi8NTLj1NZVclT9z/XYQToIBHaKM63Xsyqml/YXucZB42QxGDsdXYchYe3bDKm\npWLIyaBk/ASPzNda9pfXERVs8FpNlor36Vi/RSpe4ZIe0+gb1psv9iz1WjErQFB/E11uiaA0q7xN\n6dtnJYdyWa8mIVqwqwDn6RfgHDcB+c1nEOm7PBx1y1AUhWdee5LS8hKefqDjCdBBRplOpE9AP+aV\nfUSls+35QdaYQIQsaDhC62/b98toDI2gqv/QNs/VWhRFIb+8nlN6Rvh8bhXP0TF/k1Q8ihCCh4c/\nQYwlmvmpiyitK/PaXNbeRrreHkHF/io2zt2Co8G9s5yzUkKZ1iuSJenFLNxdhOO6O1Hik9E8fQ9U\neC/+w6EoCs++9hSFRQU889//Q5Zln87fGoQQTAm+HIC5pR+2+SxQ1soExQSiz/y3CIn6eoJ/XUHJ\nuLOgHb4mZTWN1Nud9IhRSxT9GVWEjhO0spYnRz9LgCaAuakLvFZDBE1Zc13viqC6qIb1n2zGXude\nmvjZKaFc2iuSL9OLWJhZgX3Wk02FrM/eD3bvpJ4fjkZ7IxZzIM899GKHFqCDWORALg2+kp0NO/it\n5qc2jxecEExZzr/PhYL/+BG5rpbSU9tnKy6/vB6AXtFqoao/o4rQcUSgLpAnRz9Dnb2OBTs/w+Hy\nXsaZOTmAbndHUFdWz7qPNtFQ457X3ISDQpRWxKJSCfuspxC7tiG/4btGeHqdnv9cNdMvBOggvQx9\nGWM6hSXln5Fv39+msUISg7DXOf51LhTy/TKq+g2mMTK6TeO7y/7yOqxGLUEm1bLHn1FF6Dgj2hzN\ngyMeIa9qH8vSvvJa6jaAMUFPt3sjaKyxs/6jTW67b09ICeWSnhF8kVbEQk00jpseQP7xa6QF73s4\n4s7FuUGTCNWE8VHpbByK+ytHa2wgkizR8LctOf3+XCw7NlMy/hxPhOoW+8vrVOfsToAqQschvWy9\nuH3wnWwt2s4vub95dS5DjI5u94XjaHSybs4makvd2wY8p0vYX0I0L6I/9mnXo5n3LtLK5R6OuPOg\nEzouD5lBgf0Ay9vQFlzWyFhjLOiz/5foYPt+KQ6TmfIRJ3ki1FbjdCkcKK9nbO/IdplfxXOoInSc\nMib2RKb1vIxfcn5juxdTt6HJ9LTb/eEIAWs+2EBFK/vUHOScLmF/JSu83/007KdfgPz6U4i13hVS\nfyZGF8fZ1vP4uXolmQ173R7HFGairrmhoVRTTeg3X1By+nk+beP9d/JKa7E7XQxJ6hjFwiruo4rQ\nccxF3SbTN6wPS9KWsa/Ku6219aFauj4YjinEwLqPNlG42z27m7NTQrmmXzQrskt57YSpOE44Cc1z\nDyBSN3s44s7DKebxxGrjmV/2MQ7FvXPAgMAA6iubEgHCvl6MsNspPG+KJ8NsFekFVZj0GrpHqZlx\n/o4qQscxQggeGv4YkaYI5qcu8lr7h4NoLTLJ94cQ2sXG5kXbyFnnXs3S2IQQbhoUx+r8Sp4/5Roa\nevRH8/idiKz0NsX39kdv8MOvK9o0RkdEEhJTgi+n0HGAH6u+d2sMg1WPvc6BUlVD2NIFlI47C0dI\nqIcjbTnpBdWc3DMcSVKLVP0dVYSOc3SyjidGPYMkyczfuchrrtsHkXQScXeYiTjVyq5v09i9Mt2t\n5IgRMVbuGJrAluJanjr9ZmqjE9A8fCsUuJcJ9u7Hb7Fu01pGDz/Rrfd3dGJ1cZxsPpXvKpdTaG9d\nx1SAAGvTtpvluxVoKisomHipp0NsMeW1jRRVNTCyW1i7xaDiOVQRUiE4IJhHRzxOcW0JS9KWezVj\nDpo6d8ZNsRE31Ub26ly2Lt7hVk+igREW7jkhkb1VjTw24W6qzUFoH74FylvWiO0g7819lz83rub1\nZ95Gr+u8HmRnBJ6DVQ5iYfmnrf4eBwQ2iZDxp58pHz3O5y28/87egmqEgGHJtnaLQcVzqCKkAkBy\nUAp3DL6LHcWp/Ja7yidzRpxqJeXGCIrSStjwyRYaa1u/CusVauL+EUnk17t48Pz7KXNJaB69DWoP\n73X2Tz6YN5vf1/zKG093bgEC0Et6Lgq+lLSGXayrbZ2hrd6iAxSctU4KJk3zToAtJK2girgQIxaD\ntl3jUPEMqgip/MXImFFc0mMaP+X8ws5i33i0BQ8y0e3uSGpKaln7wQZqy1qfwp0SZOChUUlUugT/\nnfggxZW1aB6+FWqqj/q+OQve55c/fuKtZ99F305ZXr6mZ0BvBhuG8WXFola1BJcE6J11lCf0pi65\nmxcjPDq1jQ4yCqu56IT4dotBxbOoIqRyCJO7T6F3aE++2LOUA9WtPztwB3NKAN0faKr3WPP+Bsrz\nWm+8GWsJ4OFRybg0Gu6f9BD7K2rRPHQzVB8+2eLjRXP44beVvPXc7ONGgA5yQdDFKIrClxWLWvwe\n6YflBNSXUpHS34uRHZud+ytRFNT6oE6EKkIqh9Bkdvo4oQYb83cuotZe65N5AyKaU7htRtZ9tInc\nDftafW4RYdLx0KhkjAY99517LzudOjT33QCH6X4aGx3L28+/R0DA8Vdxb5EDOS/oQtbV/sme+has\neIsKkD94FW2gkQadxfsBHgFFUdiSXU5SuBmbuXNvnR5PqCKk8i/0mgAeH/U0dpedz3Z5rz34P9Fa\nZJIfCCF2UDQ7v97DjqW7Wt2XKCRAy0Mjk0kMMfHIuJn8FpSM9t7r/5U1d8qocRiOQwE6yDDjSJJ1\nXVhU/unRa4ecTjQvPQwBBpzRCQhN+90yckpq2V9ex8zT2m87UMXzqCKkcljCjGHcP+y/ZFVksyLz\nB5/NK2kEEdfrSZoRxoHUQtZ+sLHV50Rmncw9JyQwOjaIV4ZOZmHXk9DMuhaR7b5jQGdDEhIXBV9C\nsaOIn6qOXBslffkpYsdmHLc9hNMJQtt+t4w/04sJs+gZ0aX96pNUPI8qQipHpG9YP2b0vZY/969l\na+E2n85tG2GhxwNROBqd/Dl7PUVpJa16v0aSuK5/DBd3D2dh91N4bcSlKA/MROze7qWI/Y9obWxT\n7VDVckoc/96yFGk7kT99G9fEy1D6DsZpdyG07eMkXlzVQFpBNdeN66p2Ue1kqCKkclQmJJ9D//C+\nLEv/mv3V+T6d2xinp/uj4QTHWdk0fyvpv2S26pxICMEF3cKZOSiWVTF9efz0W6h7fBZi81ovRu1f\nnBE4AaNk4vPy+Yd+oL4OzQsPoiR2xTn1GgBcdidC1z4i9OfeYsx6Daf1jWqX+VW8hypCKkdFCMED\nJzxCuDGMhTs/o8besvobT6ExysTfHUj0BcFk/JrFpvlbW9Ukz+l0MiomiPtHJJEdGs/9595L8UtP\nIf740YtR+w96KYCJ1slsr9/Ktrr/+e/Js1+EkiIcdzwK2qZ6HKfd2S7bcdX1drblVnD5mCR07Xgm\npeId1O+oyjHRy3oeG/kUdpfDp4kKBxGSIHpCMF1vi6RiXyWr313fIifuL79ezCPPPwhAD5uJR0en\n4AoJ497z7if9w/eRFn7gs8Z4HZn+hkH01Pfm8/IF2BU7YvVPyCuW4pxxG8T8rx7H2ehCaoftuHWZ\npchCcMGQOJ/PreJ9VBFSaRFhxjAeOOFBcipz+d6HiQp/x9rHSI+Ho9Cbdaz9cCPZa/OOuD235Nsv\n+OSzj7jzxll/vRZl1vPomBRiwkN4ZMLdrFqzBc0z97bYXaGzIoTggqCLKXeW8WvBUjSvP4VrxMm4\nxp/71zWKojSvhHwrQo0OJxsyy5g4LE51SOikqCKk0mL6hPZlRt9rWdMOiQoH0YdqSXnQRtg4C7u/\nS2PLZzv+tT237PslfLzwQ95/+SMCzYda/Vt0Gu4bnsjI+BBeOXkGn4hIpLtnQF6WDz+LjkeENopR\nxjF8X/MdVVYDjv/cC39LAHA5mla/QufbW8bm7HIaHE4mD0/w6bwqvkMVIZVWcXbSBPqH92NZ+tcU\n1PjGUeGfSBpB/JRQUm6MoDSrjD/eXkdJZpNp6VcrlvHBvNlNAmSxHvb9Wrkpc25ar0iW9jqVx4Zc\nQu19NyH98p0vP40Ox9nL8sHl5Ntbx0PgoV+76qKm1aJs9Z27RL3dye97iugbG0SU2sa706KKkEqr\nEELwwLCHsBlCWLhzMfWO+naLJXiQiV6PRWOyGdjwyRZmvzCb2Z+8zQcvf3xEATqIEIKzU0K5f0Qi\nOVFduOuC/5Lx8RzkN56BxgYffQYdB+mrz7B+/hljK7rymyGVckfZIR8vyShF1snoYn3XRG7VniIa\nnS4ev6h9rYJUvIsqQiqtRq8J4NGRT1Jjr+XLPcu83vrhaOhCNCTdH8yB7hnM/fYj/jP4duT6lp8d\n9A4189SJXbCF2XjgvHv5fl8V8t3XQH6uF6PuQCgK8idvoXnneZznTOGkgbeglwL4rmr5IZeVZJQR\nkhiMkH1zyyiraWRtRilXjkkm3IerLxXfo4qQiltEmaK4e+g97C7dw6q81rUF8DRCEgy/fACvLv4/\njFoTf85eT/aa3BaLY4hBy4Mjkzg1ycbs4VN4rfvpuO66Bunnbzp39py9EfmFh5AXfYjjipk4p99C\ngGxgvOVM/qxZRZGjEABHg4Py3AoaU4J9FtoPOwow6mQuHZXoszlV2gdVhFTcZmjkMC7uNpkfs38m\nszyrXWMJjQwlqm843R8LJ/QUM7u/T2f9x5tbbPmjkSSu7BPNjQNj+TN+APed/wBF772J5sGbIC/b\ny9G3A1UVaB68GWn1z9jvfgLXxGl/JSKMMp+ERQ7km4qlAJRml6O4FAJ8JELZxTXsyq/k9rN6YtBp\nfDKnSvvRKURICDFOCHFX85+FQoik9o7peGFqz0tJCkrks91fUNlw7NodbyPpJOKnhtLtrijqyutZ\n/fa6Vjlyj44N4rExKTSEhHPX5Cf4XROK9pZpyJ++Aw3td/7lUQ7sQzvrWkRuJo7HXkMZNe6QD+uE\njjMsE9hYt4799jxKMkoxBAUgh3g/OcClKKzYfoDoIAOnq+4IxwV+L0JCCCswSFGU5xRFeQ5YABzZ\nkVHFo8hC5pERT6CRNCza9TlOV+vbdHuDwB4Gej4RSWSfcHZ+vYeNc7dQX9kyEYkPDOCJMSkMig7i\n5UEX8tolD9OwdAHam6chNv7p5ci9i9i9He3dM8DlxP7cbJSe/Q573QmmkdjkUL6uWELJ3lJESrBP\nPNvWZ5RyoKKehyb1RZJUj7jjAb8XIWAI8PTf/r8SSBZCJLZLNMchVr2V/57wEPnVB/g2s+Pov2yQ\niJppoOutkVQX1fDHW+vYv+VAi1ZFRq3MjQNjuX5ADKsNUcya/goZCb3QPnIr8rP3Q0mRDz4DzyL+\n+BHNAzeiRMdhf3Y2RB3ZgUAWGs60nkt6UTq1pXUEpIR4Pb7Cynp+SC3g4hPi6R/vu/MnlfbF70VI\nUZQfgMF/e2lo08tKVvtEdHzSPaQH1/W7gfX5G9jixULWyvLWb/lZ+xrp+UQUYd1sbF+6k00LtlFf\ncexVkRCCk+KCefLEFHR6Pff1n8wXNz4P2zehvXEy0qIPoc43Tf/ahMOBNP89NM/ej2voaByPvgqB\nQcd82yDDULrl9UYRCvqkY1/fFurtTj5fn0eISceN49V+QccToj3Ta72BEOJ7YIGiKO8d4eODgA0n\nnngiVuuhtSRTp05l6tSpPoiyc6IoCg+uvo8dxanM6H8VEaZwj47/589reePJd3j/67fQaNw7sC7b\nWEPenDIcDU6SRsWTMDwOuQVWNHani/m7Cvgmo4QeQXpm7vmOiG8XgUaDa+zZOM+cBLEdrKq/oR5p\n5TLkLz6F4gJck6/GOXk6SC179myoaeTX11eRHb+X5EknE6FPaXNI+3MzAYiO+9+xrdPlYt7qHA5U\n1PHBdSNICjO3eR6VtjNv3jzmzZt3yGsVFRX8+uuvAIMVRdnoiXk6rAgJIa4BUoB/BiiaX1uhKMqP\nh3mPoijK7KOMOwjYsGHDBgYNGuThqFUaHPXM/PEG7C4H1wy4igCNZ2o81v66nhcffJVX5j9PWGRY\nm8Zy1LrIX1ZG4cpK9GYdXcemENknvEVnHqnF1by5eR81didXxBs4ZeO3yN8vQVSW4+o/FOdZF6IM\nHQVyO2Z11VQjfbMYeel8qKrANWY8rkmXoyS0TkS2L9lJUVoxf05ZiSXYxCmWu9oc2j9FSFEUlm3a\nz/a8Cl67cgiDEr2/7afiPhs3bmTw4MFwPIhQaxFCjAOsiqJ8fozrVBHyMvur93PLTzNJCkrk4h6T\n2nygvf73DfzfA6/w8rznCY9qmwD9nfoCO4Vz6ijcXYw1JpDup3UhKPboTgsAtXYnc7bn82teOf3C\nzFzdI5SoTb8hff0Z0p4dKGGROM+4oMkA1OrDs43yUuRlC5C+/gwaG3GdOgHnBdMgMqbVQ5XllLNu\nziYSLg8lv2cOX+xZwsTI+wnVxR/7zUfhnyL02+5CftlVxMOT+nJGv+g2ja3ifVQROgLNwhLcfD50\ncEW0QFGUfx0gqCLkG/7MX82Tax5nfOI4RsYOd3ucDas28dx9L/LS3OeIjInwYIT/o3JXHfkfVVFV\nUE1k73C6jk3BEHTsFdzGgio+2LafygYHk7qHc1ZyKNqM3U1i9NsKcLlQhozE1WcQSu8BKAldQPaw\nC7WiQH4u8lefIX2/BCQJ1xkTcZ47BWzuCbbL6eLPd9cjaSW6PhKKIhRe2/AWkaYIxphua1O4fxeh\nrbnlLN24j2vHduHqk9q+1afifbwhQn5fCdZcE7QeUJqfuAVQpijKu+0a2HHO8KgRTOp6IZ+nLSbG\nEk2CtfVP0JtWb+bZe1/gZS8KEDSlc1seD6DkDwP7F5Wz6o01JAyPJWlUAhr9kX9FBkVY6GXrwqLd\nhczfWcCqfRVc0y+eLjc/gPPKm5BWLkNa9xvyh68hHHYUowmlZz9cvQY2iVKXnn81jGsxZSWItFSk\ntJ2I9FRE+i5EZTmKJRDnpMtwnX0RHMM371jkrMujuqiGnv+NQUgCgWB07EiWpX9FH91+grVtX7Fk\nFdWwfNN++scHcdWJyW0eT8V/6RQrodagroR8h9Pl5LZfbqK4tphrB07HorO0+L35uQe49dK7eOnT\n54iKi/RilIfirHdx4JtyCr6tRNbKJI2KJ25IzDGTFzLK65i9dR9ZFfWMTwxhco8IjAff09iASEtF\n7NiMtGMzYtc2RH0tik6P0rUXSlAI6PWg1YNOBzo9ik4H2qZ/U1+HlL4TkbYTUdJkpaNYg1G69ETp\n2hNX114ovQeCwdjmz7++soFVb64hul8kETfo//d1cTl5ef1rdAlOYaj+P26Pvz83k9LqRr7JcBEd\nbODjG0ai8ZEfnUrbUbfjPIAqQr6lrL6MG3+8HltACJf3vRRJtOyGoygKVRVVBAb5zrX57zSWOiid\n52Df5ny0Ri3JoxOIHRiNdJT20k6XwneZJSzcXYhJ22QDNDTqMPE7HYiMNETqZqRd26C6EuyN0NgI\n9kZEY0Pzv5v/1mhQkrs3CU6XXihde0JY5CH9fjzF1sU7KM0uo9fT0WiMhwrvb7mr+CXnNy6JfhqD\n7N73JT09jcVr8wiKjGPezFGYA9RGdf6EKkIeQBUh35NasoN7f5vF8JgTOC1p3LHf0IFoKLJTMtfO\n/m0HCLDoSR6TSHT/SKSjPL0X1Tby4fZ8NhZUMSjCwpSeEcRZOr4TdElmKRs+2ULi9DBCR/571Vpn\nr+OFda8wwHIWg6xnt3r8qjo7s7/6g7oGJ5//9yIi1R5Bfoc3REhdB6t4nV623lzdZzqr9/1JavGu\n9g6nVejDtETfYqT3Y7EExVpJ/Wo3v7/2Jznr8nDaD29RFGbUcefQeG4dHEdeVT2zfk7n1Y257Kvq\nuH2Kqgtr2PbFToLjrdhGHL5Ox6A10D+8H7vqfsSltM6eac+BKt79eS/1jU4evbCfKkAqf+H3iQkq\n/sG5Keez+sAqlqQtI8wYSpgxtL1DahWGKB2xd+iw5ekoWWRn13dpZPyWRUz/KKIHRGGyHXoeI4Tg\nhGgrgyMt/JxTzpdpRdy1L42RMVYmdQsnyqw/wky+xeV0kf1nLnt/zcIYYiD+9qCjptQPjhzIhgMb\nya3bToLx2M3m7A4XK3ccYENWGV0jzNwyYQhBJp0nPwUVP0fdjlPxGXWOOmb+eAMuxcU1/a9Cr+kY\nN2J3qC+wU7CygrI/anHUOwhOCCJmQBQRPcMOm8Rgd7r4KbeML9OKKK93MDo2iIndwog0td/XoGJ/\nJanLd1NVWE3EaVaizwtG1h97c+TtTbOx6q2cZL7jqNcdKK/jiw15VNTZue2MHkwcGsfevXsB6NKl\ni0c+BxXfom7Hqfg1Bo2Bx0Y+QXVjNV/sWYqiKGSlZ2NvtLd3aK0mIEJLwqWh9HsxlqRrwkA0OQz8\n8uIfpH69m8r8qkOu18oSpyXaeGlsNy7vE8W2omru+CmNtzbnsa/Kty0iHI1Odq9IZ837GwDo+UAM\ncRfbWiRAAAMjBrCnNI1aZ8VhP64oCqvTi3n/10w0ssTHN4xk0rB4n7hwq/gf6nacik+JNsdw99B7\neOzPR5i3YiFfPrqUl+Y+T3S8f/aOkXQStuEWbMOhvtBC8e9VFP1aTN6G/VgizET3iySseyjG4KYz\nEJ0scUaSjbHxwazMLmVZejG/5JbTy2bitMQQhkQGInuxhUHx3hJSv95DY3UjMZNCiBhvRdK0br6+\n4b35PnMle2r+ZEDg6Yd8rLLOztKN+8gqruHSUYlcN7YruqNkFKqoqCKk4nOGRg5jJKN54p4neeit\n+/xWgP5JQLiW2IkhxJwXTMX2Wiq/c7Lnx73sXpGOKcxEeFcbYd1DsUYHopMlzkoOZXxCCGvzK1mR\nVcpLG3IJCdAwNiGEsfHBBHswfbmxtpHd36eTv62AkMQgutwdRkCEe+MbNAZ6hfZkb/XrigFgAAAc\nKUlEQVTP9FdOQwhBo8PF1pwyftlVhEYWvHrFEIYm2zwWv0rnRRUhFZ+zdetWPr7/EyY9eh5rnevp\nX9eXEEPnMa4UsiCov4mg/hBTZ6YytY7yLTXkbcon848ctEYtYV1thHULxZYczKjYIEbFBpFVUceK\nrFKWpRfxxZ5ChkVZOS0xhO4hRre2sux1dkoyyyjJKKVwVxGKAolXhWEbZW7z1tjAiP5sK9pOZtVu\ncvaFsCGzlHqHkz4xVp6/dBBWo5p8oNIyVBFS8Snbt29nypQpLFiwgKTuScz88QYW7FzM9P5XoJM7\n341LNkgEDzYRPNiE4lKo3ttAxZYaKtZVsn/LAYQssISbMYeZMIebuCDcxKThSawpq2FFdhmP/JFJ\njFnPkEgLgyIC6RJsQDqCgLicLir2V1Kyt5SSjDIq9leCAiabkZDRJiLPtKK1euZX3uCKQ6tY+Srj\nW2rzTmfi0FgmD08kOlhNvVZpHaoIqfiMHTt2cPHFFzN//nz69u0LwKMjH+e2n29hadpXTOp+fqc+\nvBaSwNI1AEvXALiwKcOuYnstdTmNVGfUcGBnIS67CwB9gIYpYSYaLDr2FdSyL6ucDKeCXhJEm3RE\nGnWEG3XoJYGiKFQX1VKW3dQnSWvQEJIUQtB4HYG9jOhtnvk1d7kUduY4+W27nW2ZDoJje2EIW8ec\ni/5LZBv96lSOX1QRUvEJLpeLmTNnMnfuXPr16/fX6/GBCdw2+A6eXfc0sZYYhscMa8cofUtAhJaA\niP/dvBVXKI0lDur2NTb9yWtEznIQVecgwuHECThcCvaiWsoVKAdkSaCVJfSBegwnW4gcZCQwWY9o\nY3KDoiiUVinkFjnJK3KRW+Qku9BFVa1CtE3i0nEBdE8czGsbV7OzbAORlrFt+2KoHLeoIqTiEyRJ\nYuXKlciHaWUwOmYMaWV7WLL3S6LMkW45bncGhCTQh2nRh2kJGmA66rUVNS525jjZmeNgV46T6noF\nCkB8YyfI4iDMKhEeJBFmlQgLEtgsEkKAS2nq/qAoTf92uZS/XiuvbhadYhd5RU5qmw0ezAZBXJjE\n8B5aBqRoSIiQmlesVuIssXyVuYRT4lURUnEPVYRUfMbhBOggl/e6ko2F6/ls1+dcO2A6Fn3LHbeP\nR6wmieE9JYb31KIoCpW1CoXlLooqFIrKXRSWu8g84GTtLjuNjpaPawsUxIXJjB2oIzZUJi5MwmoS\nR9wm7RXak5VZP1Jjr8GkPbpwqqgcDlWEVDoEsiTz6MgnmfnjDSzctZgr+k5DI6k/ni1BCIHVJLCa\nJLr+o4mqoihU1CiUVTc5o0iiyXxbiKZ/SxJIQiBE04rHqG/dNl6v0B58l7mCtflr1NWQiluoVWQq\nHYbggGAeHvEo+dUH+HrvtxxvllLeQAhBkFkiKVImKVImIUImPlwmLkwmJlQmKkQmIrhp6661AgQQ\nqA/8a0tORcUdVBFS6VB0De7GzIE3s6lgC+sPbGjvcFRaQK/Qnuwtz6TGXtPeoaj4IaoIqXiUzMxM\nNmxom3iMiz+Vc5LP49uMFWRX5HgoMhVv0Su0B07Fydr8Ne0dioofooqQisfIysri3HPPpb6+7Yac\nV/W5mvjAOBbuWkxF/eGNMlU6Bn9tyWUtbe9QVPwQVYRUPEJ2djbnnHMOb775JqNGjWrzeBpJw+Mj\nn0IraVmw8zPsTv9z2j6e6BXak71lGdTaa9s7FBU/QxUhlTaTk5PDOeecw+uvv87o0aM9Nm6g3soj\nIx6jqK6Y5enfqIkKHZiuIV1wKk5SS3a0dygqfoYqQiptIjc3lwkTJvDKK69w4oknenz85KAUbhl4\nG1uLtrE2f73Hx1fxDCEBwQTqLGwt3tLeoaj4GaoIqbjNvn37mDBhAi+99BInn3yy1+Y5Ke5kzku5\ngO8yVpBVnu21eVTcRwhBUlAif+b/0d6hqPgZqgipuM3s2bN54YUXGDvW+0WKV/a+ikRrAot2fa4m\nKnRQEgLjKagpVM+FVFqFKkIqbvPQQw8xbtw4n8wlSzKPjXwSrawmKnRU4gLjUFDYXbarvUNR8SNU\nEVLxGwL1Vh4d+TjFdSV8uWepmqjQwbAZQjBoDOwsSW3vUFT8CFWEVPyKJGsydw2ZRWrJLn7O+bW9\nw1H5G0II4gJjWXNgdXuHouJHqCKk4ncMjx7BFb2u5Nfc39lWuL29w1H5G3GBseRV7cPpcrZ3KCp+\ngipCKn7JxK4X0j+8L0vSlpNXta+9w1FpJs4Si91lJ7cqt71DUfETVBFSOSrFxcW8+OKL7R3GvxBC\n8OAJjxFtiWJ+6iIqGirbOyQVINT4/+3deXhU5b0H8O+bSWaykRkSEgJCSQJY3EJZC9ULXkErXFuB\ngIjKtShoeSgVK2DoY63thSYsdblel4jlylZkNRK4hCKLFLVoQLZCBQzBQFQgmYxZZkgyee8fc8YO\nMWkmMDPvnJPv53l4yJycOfnOM8tvznve8ztJAICyan4xIP+wCFGLysvLMWrUKKSnp6uO0qwoUxQW\n3LoQkRGRWHN8LercdaojtXuxkbGINkXjfA2LEPmHRYiaVVFRgdGjRyM7OxtjxoxRHadFNosN//Wj\nP8DuqsQ7n73LGXOKCSGQGNMRZdVlqqOQTrAI0XfY7XaMGjUKc+bMQVZWluo4rUqzpuHpQfPwj4qT\n2HV2j+o47V5STBJOVp5QHYN0gkWIrlBZWYlRo0bhqaeewvjx41XH8dug1MGYctMj2HfuQxy+cFR1\nnHYtKSYRFU676hikEyxC9C2Hw4HRo0dj1qxZuO+++1THabMxvcbhBymZKDi1FV84ODtLlcSYRNTU\n16C6rlp1FNIBFiH61vHjxzFz5kzcf//9qqNcFSEEfjPk9+jW4Tq8fWI9LtVeUh2pXUqKSQQAfFnD\n40LUOhYh+tbQoUMxadIk1TGuSVREFHJuW4R4cxxW//1tfhtXICnaU4Q4OYH8wSJEhhNv7oDc2xbD\nLd1Y/fe3cbnhsupI7Yol0gKLyYJyV7nqKKQDLEJkSMmxKZh/ay7srkqs+8dGtpEJMYvJDFeDU3UM\n0gEWITKsdGs6fjPktyhxnMXmU1t4DlEImU1mOFmEyA8sQmRomcl98asBs3Hk4jHsPLtbdZx2g0WI\n/MUi1M7U1tZiypQpqKmpUR0lZIZ1G45Hb56KD859hP1ln6iO0y6wCJG/IlUHoNBxOp249957MW7c\nOMTFxamOE1L39hqLclc58k+/g3hzPG7qdIPqSIZmMZnxdS1nx1HruCfUTjidTowZMwZjxozB9OnT\nVcdR4mc3PYKbk2/CO5+9ixLHWdVxDM1sMuMyG8qSH1iE2gGXy4WxY8finnvuwYwZM1THUSZCROB3\nQxbgewnd8fbx9fiy+kvVkQzLbDKjrpFFiFrHImRw3gJ09913Y+bMmarjKBdlikLObYuRFJOIlcfW\n4ELNBdWRDMlsMqPOXa86BukAi5CBXb58GVlZWbjrrrswa9Ys1XHCRmxULJYMewFWSwJWHPszLrK9\nT8DVNzYgMsKkOgbpgCGKkBBihBAiSwgxTQixTgjRT3WmcNDQ0IDJkyfjySefVB0l7MSbO2DJsBcQ\nFxWLFUdXodzJs/sDyVnvRGxkrOoYpAOGKEIAdgA4LaVcqv28XnGesBAXF6fbZqShkGCxYsnwFxEd\nGYPlR1fz8gMB5GxwokvcdapjkA4YpQhlSCkP+9zmpwn5xWaxYcnwF2A2RWHFsVWodFWqjmQIznon\nEswJqmOQDhjiPCEpZYnPzZ4AJrR2nyeffBJWq/WKZZMmTdJ9F2lqu8ToRCwZ9gJ+9f4sLD+6Gj/L\nnAyrhR+g16K2wYkO5g6qY9A1WLNmDdasWXPFMofDEfC/I4zST0sIkQ7gcQD9AGRLKT9tYb3+AA4c\nOHAA/fv3D2VECnMXay/gqfdnIUJE4OFbHkICC9FVy/loMR7sMxlje4+7Yvnp06cBAL169VIRi67R\nwYMHMWDAAAAYIKU8GIhthu2ekBBiGjx7NU2rpNCW7ZBS7vIulFKeAZAthJgDYKcQIk1K+U3IAisk\npYQQQnUM3UuOTcGiYc9j9vuzsOLYavzslsmIN8erjqU77kY36tx1SOCeEPkhbI8JSSmXSimzpZTz\nmvzzLtsFePaAhBC5PnfdAMAGYKCS4CHW0NCABx98EB9++KHqKIaQGpeKRcOeR527Dm8dXQmHK/DD\nD0bn7RkXzyJEfgjbItQGGQBG+NzuCc/EhCI1cULH7Xbj4YcfRnp6OoYOHao6jmF0je+KxcNegLvR\njWVHluNi7UXVkXSlqq4KAGC12BQnIT3QfRGSUu4EkCeEmKoN4U0DcIfRh+LcbjemTJmC7t27Y/78\n+RyOC7Cu8V3xwu0vIzoyBsuOrEDpN+dUR9KNc1VliBARSE9IUx2FdED3RQgApJRvav+WSiknNpmu\nbThutxuPPPIIUlNTkZOTwwIUJEkxSXjx9peREpuMFcdW42TFKdWRdOHcN+eQGtcZlsho1VFIBwxR\nhNqTxsZGTJ06FcnJyVi4cCELUJDFm+Pxx+EvoactA28fX4/DXx9RHSnslVadx4CUwapjkE6wCOnM\n9OnTYbPZsHjxYhagELGYLMi9bQl+0DkT+acK8OG5v6mOFLaq66phd9nRJ7GP6iikE2E7RZua9+ij\nj2LQoEEsQCFmijDhuSELsPrESqw7uRbV9TW4M+0OPg9NnKs6DwDok8iLBpJ/WIR0ZvBgDnOoIoTA\nQzf+J2wWG944mofa+hr8pNd/wMRu0d8q/eYcEswdkBybrDoK6QSLEFEb3dPzp0iwWPH8gSWwuyox\noc84ntSqKa06h+4J3VTHIB3hMSGiqzCs23Dk3LYQ5c4KLD20DGVVZaojKdfQ2ICyqi8xtMu/qY5C\nOsIiRHSVbki6Ef9zx6uIN8dj2ZEVOPjVIdWRlDpVcRpu6UZmp76qo5COsAiFISklampqVMcgPyTF\ndMJ///ur6Ns5EwWnt6Lg1FY0NDaojhVyjbIRu794H2nWHkizpqmOQzrCIhRmpJSYPXs2XnzxRdVR\nyE9mkxnPDZmPmf2ewOELR/G/R1a0u55zB776FBdrL+GJfr9SHYV0hkUojEgpMXfuXFy6dAnZ2dmq\n41Ab3dnjLiwe/kfU1Nci79CfUFx5RnWkkHA1uLDn7F70TclETxsv0UBtwyIUJqSUmDdvHr766iss\nW7YMJhOn/epRL1tvvDLidXSN74JVx9Zg7xf70CgbVccKqr+WfoD6xnrM6jdbdRTSIRahMCClxDPP\nPIPS0lK89dZbLEA6l2BOwJLhL+K+70/Eni/2YtmR5bhUW646VlDYXXbsL/sEE66/D0kxSarjkA6x\nCCkmpcSzzz6L4uJiLF++nAXIIEzChAdvmIyFwxbDWe9C3qE38bfzH8MoVzL22nFmF2KjYjG217jW\nVyZqBouQYq+++ipOnjyJlStXIjKS5w4bTZ/EG5A3cilGpY3G9jM7sPzoKthdlapjBcRZxxc4Uf4P\nTL35MXbMpqvGIqTYxIkTsWrVKhYgA7NERmNa5uOYf+sfUHnZgdcOvoGPy4p0vVdU567Dts+3o2t8\nFwzvfrvqOKRjLEKKderUCVFRUapjUAhkJvdF3sg3kZlyM7YVb8dbR1ei3Km/Y0XuRjfWndgI++VK\nzB04DxGCHyN09fjqIQqh2KhY/HbIfCy4NQdVddV47eBSvFeyC7X1taqj+UVKiXdPFaDEcRa/GfIc\nMmw9VUcinWMRIlLgluRMvDHyTWT1noCPy4rwUtEr2H32fbgaXKqjtUhKie1n3sPRi3/HUwPnoG8y\n2/PQteOBCCJFLJHReOjGyfhJz59i06mNKCh+F/vLPsHQ636IIV0HwxJpUR3xW1JK7Dv3IfaXfYyf\nZ07HbdexSSkFBotQiBQXFyMjI0N1DApDVosVU25+BGN6jcGGkxvwf2e2YH/ZJ/jRdUMwuOtAmE1m\npfnKnRUoLN6O0/Zi3P/9SRidcY/SPGQsLEIh8NJLL2HXrl3Iz8/nlTipRR2jEzEt8zGM7T0O60+u\nxfaSQnxUth8/um4IMlNuRgdzh5DmqXPX4a+lH+Cj8/vRwRyPXw9+Bj/sMiSkGcj4WISC7OWXX0Zh\nYSHeeecdFiDyS6eYTpjedwayek/AK4dfwq6ze7CzZDfSbWm4Jfkm3JDUJ6hDdVJKHL90An858x5q\n6msx4fqJyLp+PCym8BkeJONgEQqiV155BVu2bEF+fj6io3kyH7VNSmwKfjd0AarrqvBB2QfYUvwu\n3j21BVs/L8T1ib1xS/JN6NWxJyIjAvc2vlBzEYXFf8EZRwm+n9gbTw14GqlxXQK2faKmWISC5LXX\nXkN+fj42b96MmJgY1XFIx+LNHfDjtLvx47S7cbH2Ivaefx+FZ7Zi7YkNsJgs6JN0PbondIfNYoUt\n2garJcGvwtTQ2ADHZQcqXQ7YXZX4rOIznLYXIzG6I54d8hwGpg4KwaOj9o5FKAjy8vKwceNGFiAK\nuOTYZGT1Ho+s3uPxxTdnse/8X7GzdAcOXzh6xXodzB1gi7bCZrHBFm1FXFQcquuqUXnZgUpXJSov\nO1BdV/3t+gICXeJT8US/WRjW7XZEmXgCNYUGi1CAHThwAGvXrkVBQQFiY2NVxyED+15CDzyQ0AMP\n3PAQ6hvrUe68hK9rL+BC7dfavwv43HEKZx1nUV1f4ylMFiuut/VBSmxndI5NQUpsZ6TEdkZSTFJA\nh/WI/MVXXYANGDAAhYWFMJvVTqul9iUqIgqpcV14/IZ0hx0TgoAFiIjIPyxCRESkDIsQEREpwyJE\nRETKsAhdpT179qCurk51DCIiXWMRugrr1q1DdnY2nE6n6ihERLrGItRGGzZswPPPP49t27bBarWq\njkNEpGs8T6gNNm3ahEWLFqGwsBAdO3ZUHYeISPe4J+Sn/Px85ObmYtu2bUhMTFQdh4jIELgn5IfN\nmzdj/vz5KCwsRFJSkuo4RESGwT2hVrhcLuTl5WHbtm3o1KmT6jhERIbCPaFWREdHY+vWrapjEBEZ\nEveEiIhIGRYhIiJShkWIiIiUYREiIiJlWIQ0u3fvRlFRkeoYRETtCosQPM1In3jiCaSmpqqOEhBr\n1qxRHSGo+Pj0raCgQHWEoDL68xdohitCQojXhRAJ/q6/d+9e/PKXv8TWrVvRrVu3YEYLGaO/Cfj4\n9G3Lli2qIwSV0Z+/QDNUERJCjAAwAYBffXX27duHGTNmoKCgAN27dw9uOCIi+g7DFCEhhLeldYU/\n6x86dAjTp09HQUEBevToEcRkRETUEsMUIQATpJQ7AQh/Vl6wYAE2b96MtLS04KYiIqIWGaJtjzYM\nt87P1aMB4Be/+AXsdjvsdnvwginicDhw8OBB1TGCho9Pv0pLS1FVVWXYxwcY+/k7ceKE98foQG1T\nSCkDta2AEkJMA9ATQNOAQlu2Q0q5SxuGG6jtBUEIcRrASCllSQvbfQDA6qAFJyIyvgellH8OxIbC\ntgj5SwiRBcB7hTkBIA/AXADvSSkPNbN+EoAfAygB4ApRTCIiI4gGkAZgu5SyPBAb1H0RakoI0Qgg\no6U9ISIiCh+GmZgghLAKIebAM1T3tBAiTW0iIiJqjeH2hIiISD8MsydERET6wyJERBQA2iQpaiMO\nx/kQQrwOYK6U8hvVWYi8hBBZUsqNqnNcCyFEOoDxAIoBpANYKqV0qE0VGFrxSYRnZq7NiJ8f2rmY\n/bWbgwA8LaU8E5Btswh5+JzwOsAoM+u0x2SD5w1yJ4AcKeWnalMFTjDfGOHASB9uQogiKeVA7Wcr\ngPVSyrsUxwooIYQbQEc9P0/N0Z6vx6SUi7XbWQAWSil7BWL7HI5D2/vO6cgOAKellEu1n9crzhMw\n2nPWX0q5WHtzrIXnMRqGlHKj9tzp+puiEKIffB6Dtgc00IAzWP1qGaZDAwHk+tx+D0BGoJ4/FiGP\nNvWd05EMKeVhn9tG6lEU1DdGmNH763IgvvsFrwJAhoIs1EbaZ+MAn0WDPIsDM2LU7otQG/vO6UqT\nF0lPeC5zYQjBfmNQQNmaWVbZwnIKQ026z8wF8Figtm2IBqZNtbHvHIw2hutLOyD8OIB+8LQ3KlEa\nKICC+cYIFn9fmyEPFlyV+O41vmzaciPR9bCpP7TX7zop5Z8Cts32PDGhrX3nwsHVfohp3STmAUgL\n56J7NY9Pu4+UUr4ZmpShpfcD3toxoTeklIN8llXAc0yvRFmwANNahul6Asm/oo0aWaWUmwK63fZc\nhJoyUt857x6QlDLb5/bn8HQYN8w37WC9McKJET7chBCnpJS9tZ9t8HyZGNTK3XTBZ5ZmLoBFMODe\nrBCiPzxfhLxXK5gGYG0gXpOGHI5rK+8URPyz79xCAxSiDAAjfG73hGdiQpGaOIGnvTHgLUCBfGOE\nA58PNwlgnhBCzx9uE4QQswGcgWeigtGOT+4EsFh1lmDQvsAWAZBCCMAzKmHXZm5e+/a5J2RcQoip\n3h8BjATwhyaz5XTLZ8/O+wL2vjGS1KUiorZiESIiImXa/RRtIiJSh0WIiIiUYREiIiJlODuOyACa\n6VK90UjNXMm4ODGByACEEKe9XY21Uw52ertWE4Uz7gkRhTHtXKHH4ZmK7m0CukNKuUn7nYTn/K8r\nulQLIaxCiDQDnO9GBsciRBRkQoh18JxE6/eF6YQQGfA01v0cwFTfE3CFEOOEEAsBzIan7dREfLcP\nmwOeE11Lri09UXBxYgJREGnHaioAvNGG+/SH5wz116SUE5t2gNA6RPQDUKz9zobmr4XVtGkoUdhh\nESIKrnQATwPo6NPBokVaX7X3ALzeSqfi9dp6QMtdqo12kUYyIBYhoiDRukcXaVcSXQRgoR93Ww+g\nXEr561bWq8A/ryRbhOavzXPQ36xEqrAIEQVPus9QWg48e0OzW1pZCDESnqazuS2t4+OAt3GrlPLT\nJtuxgRf4I53gFG2iINBmrn3SZEJBLoA5aOHaQEKI9QDGtfT7Vv7eD+CZoFAET5fqPBYh0gMWIaIg\nEEJMbXqRPe38HTuAhVLKec3cpwKeobjeIYpJpByH44gCTLti77qmy32ODc0VQiQ0c1cbeByH2hkW\nIaLA+1fDaTnwXPuouUkKlfBjRpt2qXYiQ2ARIgqglvaCvHz2hh4TQqQ1+XURPFfEbc0VY+hCiH5C\niBwhRJb2v7VtqYnU4TEhogASQsyWUi5pZR3vsaE8KeV0n+UjAPwFQM+WJhUIIXKllNlNtnVGSpmo\n3U7XtnvXNT8YohBgESIKECHEXHj2UopbWxVANjxdD64oOEKIHHguxT5S22vyLk8HkAXgjSYz7rIA\n5PpOZhBCNAKwtXWGHZEKLEJEAaLNbmvrUNgGKeXEJtu5A8DP8c+CVg5Pi55NzfzNEfB0V2hahPpL\nKQ+1MQtRyLGBKVGAeIfEArCdXQB2+bnuTiFEpbdjtlaUAM+xJRYhCnucmECkc1LKQQDu1PagvEOB\nrQ0JEoUFDscR6ZwQwuo9fqS17PlcSpmkOBaRX7gnRKR/dp+TX7Ph6dpNpAvcEyLSOe0SERUAkuBp\nXPpmK3chChssQkREpAyH44iISBkWISIiUoZFiIiIlGERIiIiZViEiIhIGRYhIiJShkWIiIiUYREi\nIiJlWISIiEgZFiEiIlKGRYiIiJRhESIiImX+H+qo95w0jLXhAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f8b9fcea390>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%%time\n",
"plt.figure(figsize=(4.2,4.2))\n",
"for i in bins:\n",
" flavio.plots.band_plot(fits[i].log_likelihood, -4, 2, -4, 2, n_sigma=1, interpolation_factor=1, steps=40,\n",
" col=i, label=bin_name(bins[i]))\n",
"flavio.plots.flavio_branding()\n",
"plt.legend()\n",
"x=np.linspace(-4,2,100)\n",
"plt.plot(x, x, c='k', ls='--', lw=0.6, label='NP-like')\n",
"plt.axhline(0, c='k', lw=0.2)\n",
"plt.axvline(0, c='k', lw=0.2)\n",
"plt.xlabel(r'$\\Delta C_9^0$')\n",
"plt.ylabel(r'$\\Delta C_9^-$')\n",
"plt.tight_layout()\n",
"plt.savefig('DeltaC9-binbybin.pdf')"
]
},
{
"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