{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Multi-bin Poisson" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import logging\n", "import json\n", "import math\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "import pyhf\n", "from pyhf import Model, optimizer\n", "from pyhf.simplemodels import uncorrelated_background\n", "from pyhf.contrib.viz import brazil\n", "\n", "from scipy.interpolate import griddata\n", "import scrapbook as sb" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def plot_histo(ax, binning, data):\n", " bin_width = (binning[2] - binning[1]) / binning[0]\n", " bin_leftedges = np.linspace(binning[1], binning[2], binning[0] + 1)[:-1]\n", " bin_centers = [le + bin_width / 2.0 for le in bin_leftedges]\n", " ax.bar(bin_centers, data, 1, alpha=0.5)\n", "\n", "\n", "def plot_data(ax, binning, data):\n", " errors = [math.sqrt(d) for d in data]\n", " bin_width = (binning[2] - binning[1]) / binning[0]\n", " bin_leftedges = np.linspace(binning[1], binning[2], binning[0] + 1)[:-1]\n", " bin_centers = [le + bin_width / 2.0 for le in bin_leftedges]\n", " ax.bar(\n", " bin_centers,\n", " data,\n", " 0,\n", " yerr=errors,\n", " linewidth=0,\n", " error_kw=dict(ecolor='k', elinewidth=1),\n", " )\n", " ax.scatter(bin_centers, data, c='k')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "tags": [ "parameters" ] }, "outputs": [], "source": [ "validation_datadir = '../../validation/data'" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/srv/conda/envs/notebook/lib/python3.7/site-packages/pyhf/infer/calculators.py:352: RuntimeWarning: invalid value encountered in double_scalars\n", " teststat = (qmu - qmu_A) / (2 * self.sqrtqmuA_v)\n" ] } ], "source": [ "source = json.load(open(validation_datadir + \"/1bin_example1.json\", encoding=\"utf-8\"))\n", "model = uncorrelated_background(\n", " source['bindata']['sig'], source['bindata']['bkg'], source['bindata']['bkgerr']\n", ")\n", "data = source['bindata']['data'] + model.config.auxdata\n", "\n", "init_pars = model.config.suggested_init()\n", "par_bounds = model.config.suggested_bounds()\n", "\n", "(\n", " obs_limit,\n", " exp_limits,\n", " (poi_tests, tests),\n", ") = pyhf.infer.intervals.upper_limits.upper_limit(\n", " data, model, np.linspace(0, 5, 61), level=0.05, return_results=True\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "expected upper limits: [array(1.07644221), array(1.44922838), array(2.01932904), array(2.83213651), array(3.84750318)]\n", "observed upper limit : 2.381026330918668\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmUAAAGxCAYAAAAjyMGDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAACID0lEQVR4nOzdd3hURRfA4d/d9ISE3kOH0HvovTdBVEQREBGxgIAifmJXsPcugooIWFF6772F3nsn9BAgpG7u98cEBQSySW7Z7J73efaJu3t35lAkJ3NnztF0XUcIIYQQQtjLYXcAQgghhBBCkjIhhBBCCLcgSZkQQgghhBuQpEwIIYQQwg1IUiaEEEII4QYkKRNCCCGEcAOWJWWapv2oadoZTdO23+Z9TdO0LzRN269p2lZN02pZFZsQQgghhN2sXCn7CWh/h/c7AOXSHo8D31oQkxBCCCGEW7AsKdN1fRlw4Q6X3A38rCtrgFyaphW2JjohhBBCCHv52h3AdYoCx657fjztteibL9Q07XHUahohISG1K1SoYFpQBw8eJCYmBn9/H6pWLQcEc/bsOQIDAwkNDTVtXiGEEEJ4ng0bNpzTdT3/rd5zp6TMZbqujwZGA0RGRupRUVGmzTVgwAC+/fZbKtR1ErV8N6mp/hQsmEqXLpX58svngQb07z+Se++9lw4dOpgWhxBCCCGyP03TjtzuPXc6fXkCKHbd8/C012z1zTffULJiLrYdgnkH4Pz5JE6cSGHEiC3Ag1y8WIL588dy8OCrwPdcvryVnj17sn79ertDF0IIIUQ24k5J2TTg4bRTmPWBWF3X/3Pr0g73dquPfgI6PwMlS8GOHZA7t3ovVy44fNjJk09uAPpz8GB1Fi36ndjYV4EJ7Nq1lL59+7Jv3z77fgFCCCGEcHuW3b7UNO1XoDmQT9O048DrgB+AruujgFlAR2A/cBXoa1Vsd6LrOr/9oFa9ki6Cozo48v33Oh8f9bV6dTh50omuzwXmsn8/TJ/uw8iRGtCDFSt8Wb06igEDBhASEmLVL0MIIYQQbk7Tdd3uGLLE7D1lADVqVGLnkV0kdwQioGQuWNwTAlOgUKH0P5+aCo60NcnXXvPnyy91zpz5Fj+/u1m0aCuaptGiRQszfwlCCCGEcAOapm3QdT3ylu9JUuYKJ49N8+OHTf/+XoVNhvAU2LQR/P0zNlpMzLXbnw5atMhBbGwONm6cB1Rm9+7dlClTBj8/PyN/AUIIIYRwA3dKytxpT5kb86FGoTxwEVWoA7hUA1IaQGomfgev7UeDVKZPv8Qvv5wEqpCSUplmzSJ57LEH/rk2uyfNQgghhHCNJGUuGDt2LM80Og8/AlMAHSgFe8PhwUmwew8kJ2du7Bw54N8yazsZMyaOp56aDNTh3LkRlCpVnJkzZxrwqxBCCCGEO5OkzAWFChWibNkwKAGcA878+97UKKheG15+Oevz+PpCly5Qvz5AFBcvvk716scJD38V+IZdu1bz4YcfEhsbm/XJhBBCCOFWJClzQYcOHdi96zmKdwc04PqW6qGQ1Ax86xk/b9myMHUqVK++CRjI7NlNeOWV4SQnLwHgwoULpKamGj+xEEIIISwnSZnLSlCjBFAS2IG6hXlNXXh/B0zZBYcOmRfB0KFODh9OJV++rkAEffs2oXFjE7JBIYQQQlhOkjIXpKamUrz4/5jeByiAaqt+/qZrdHhgKFSpCgcOmBdL4X9atO+jZ8+dPPbYJuB+dH0uQ4c+y5o1a8ybXAghhBCmyZa9L63mcDho27YFO478zppSQCMg7L/XJVWCQA0Cb1Fc1gzduwM4gUmcOjWJ8eMdlCt3gvr1x+J0BpKamiqlNYQQQohsQlbKXPT99xP57U8fKMotEzJQr1+qB/f8ASdOQ3y8dfEVLgxHj6by6KN/AsX4++9ulClTkgNmLtsJIYQQwjCSlLnMhyIhRQmNQ92+HA8cu/WV649AxVrw1FMWhgcEBUFAAEAMhQpNoWnTU5Qs+SawjZUrV3Ls2G0CFkIIIYTtJClz0WeffUZQ0FEuf4zq2HmYG09hXs8PLlcHvaZl4f1HkyYwYUIqPj7j0fVq9OvXloce6mhfQEIIIYS4I0nKXFSjRg0efrgMLQYCgUA51CnM21WkqAc/X4SJW+HcOauivDVNg7lzr/Lll9uBely9OoUhQ4Zw9OhRewMTQgghxD8kKXNR8+bN+fHHXvTui1opqwxcAdLJax55D4qXhA0bTA/xjkqUgBo1ANaxdu09fPfdlxw58re9QQkhhBDiH3L6MgNSU4uTNx5V0T8C9bu3HVW77DZSwsG3CgTktyJC17RoAceP6+TL9ywwhffeK8+JE/58/vnnOBySpwshhBB2kO/ALkpJSSE0dCADHwBtJhAANAAKpvPBEEjoAH3mQFxi5ntkGi3fP2U7lnL+/GjOnPkDh2MtAMnuEqQQQgjhRSQpc5Gvry8vvjiAoUOh9ANpL7YC6rj2+Y3HIaI+DB1qVoSZ9+GH8NtvZ4CGnDzZmtKlizNr1iy7wxJCCCG8iiRlGfDKK+/z7LM+NKh/3YtJwFYgLp0P+8DJHBCT07z4skLT1NeEhIVUr36aChX+AGJITEy0NS4hhBDCW0hSlgG67sOJEwXJcZx/2yxNAf4G1rkwQGv4MxDWHAddT/9yO5QuDTNm6JQuPQ4oy5NPNuPuu7ugu2vAQgghhIeQpCwD3njjDYoVO8noZ1CrYwDV076uQ3U8SkeSEzq/A7XqwIULpoRpoAvUqLGWOnXWoGnzAIi3sk2BEEII4UUkKcuATp068eWXdfljGlAv7cUyqBIZ8cA+18Y5lwz7z8H5GFPCNNSQIfDKK2eB9mzc2JhixYqwdOlSu8MSQgghPI4kZRlQt25dBg5sz32doOi1U5e+qJplAGtdHCgcrvSBb/cbH6OZQkJW0qZNLNWrTwUSSE29XeVcIYQQQmSUJGUZoOs6p0/nYutWyL0TSEh7o2ra10P8u9csPQ74dCV06gPTphkeqinKl4dff9XJletTdL0qHTvW56OPPrI7LCGEEMIjSFKWAcnJyRQu/Bzvvw/bfwBOpb1REugNaEBUBgZMhTkLYfpCgwO1QHz8fvLlW09Y2BTgos3RCCGEENmfJGUZ4O/vz5gx7zBwIHw5ByiR9oYPam9ZRWAT4GrtVT9I7QNLykJsQvqXu5PgYJgwAfr3XwlUZOrU//HUU08RF5debRAhhBBC3IokZRnUr99zNGzoQ9tI1MrYNQlAYtrX7RkY0B/2X4B7R8HIke5bKuN2VH2zU+zc+SFr1/6Cv/8Zu0MSQgghsiVJyjLo4sU4Nm/Ox9HN4H/9xn5/VE9Mf2B9xsddNB3e+xCOHzckTMu9+CKsWXMJP7/qJCR8yauvvsqlS5fsDksIIYTINiQpy6CPPvqIyMjTzJkNziXAtQOIDqAK6tblSeBEBgduDMlPwmkf42K1mr8/wGWWLh3Mu+++zZo1c+wOSQghhMg2JCnLoB49evD338148UXo+ws3/g5WAXTUHrOMbPhHjZMcAg9Mgol/QHau0dquHezbp9O27SBgJlFRUdLkXAghhEiHJGUZVLlyZbp0aULevFA7/KY3iwC5gRBgG6qgbAYd3A29HoCvv85yqLYqVQrgDKdP30WzZg14/vln7Q5JCCGEcGuSlGWQ0+kkKgqOHIGtfwGbr3tTA+qgSmSk3PSeqwoCvSCkSRYDdRMFC8K4cSkMGzYH2EBiYqL00RRCCCFuQZKyDEpJSaFOnbcYNw7WLQAO33RBQ+BeoBhqw39mit6XhecWwJoD2Xfj//W6dYPw8ANAAx59tB7du98viZkQQghxE0nKMiggIICZM8fw8MOwZjVE9LvFRalAWeACqsp/JsQnQ4vWcO992a9Mxu3oejK1a2+hdu0NaFq03eEIIYQQbsXX7gCyo44dHwaeAFKpXhD23txaaRWwGAhErZaVycQkGiQ0hLyVrtUCy/40DYYOBbW8WIuoqDeJitJ54okn0DzlFymEEEJkkqyUZcLBg8eZPTs3J07A3jHAkZsuqJD2tTCwB4jN5EQRMCcFJmyFq1czG627Os3o0U/x7rsvcuXKFbuDEUIIIWwnSVkm/PDDD3TpcoHAQDi5Hbh80wX5gEKo05c6sCFr8z32LpQsDceOZW0cdzNqlM6KFRcJDX2E1NSLHPeEDXRCCCFEJklSlgn9+/dn3br25MoFm3aj6pPdrAqqYXlJYCPqNGYmJeaH5GLg45/5MdyRwwHFigH8zWefRVC5ckUOHDhgd1hCCCGELSQpy4SSJUtSs2YtfHygaBjkC77FRZXTvuYCrgC7szBhXrh4F3y0OQtjuLl77z3LM88kUrr02vQvFkIIITyQJGWZkJiYyJQpF9i1C2bMgKTvgKSbLsoNPAp0QiVmmeiHebNPF0KT9rBxY9bHcjclS8KbbyajaT05depRevbswYULF+wOSwghhLCMJGWZoOs699zzLZMmgZ8f5MwB3GojfnHAD1VQ9ggZ74d5Mw1WR8H6zVkcx81t3DiW2bMncfToNrtDEUIIISwjSVkmBAYGsmnTFAYPVn0e3/kZtRp2Mx1YeO1DwIosThwMzqdgcVgWx3FzHTvC4cMp1KjRH9jLyZMn7Q5JCCGEMJ0kZZlUo0YHcuZUtbVqFb7NRRpwDLXRPxLYBZzN4sS+8PsOeOl7WLgw/cuzq7AwgH3Mnl2b0qVLsnjxYrtDEkIIIUwlSVkmRUVt5ccfcwIw7iPw+fk2F1YBzgOlUaV6VxkweSq8/yq88rrnVPu/nfr1rzBggJP69ffYHYoQQghhKknKMumvv/7iiSdicTrVJvUCpVC3K29WEbVidhCoCWwh88Vkr3FA6gPg2/vWU3qS3Lnhk09SCQp6isTE53jzzTe46nmVdIUQQghJyjLrueeeIzr6HhwOeOopePB5VPJ1sxBUm6XtQANUFrXGgAByw4pT8P4y2LTJgPGygcWLP2HEiDdZunSu3aEIIYQQhpOkLJPy5ctHvnzl/+lLGVkEcN7m4hpAUSAYdTsziluf1syEV16Exk3g9GljxnNn7dvDnj3QocPbwEmSkm6uQyKEEEJkX5KUZVJcXBxffHGQDWktlN59GJh+m4urAPejTmA2BpIxpG4ZQGo9yPMA5MxrzHjurmxZgA1s3lyTiIhSrF0rxWaFEEJ4BknKMsnhcDBkyO/MTbuT1vMBCCqXzofOovpiRqBuYRqx0JMLjheH4QsgPt6A8bKJXLnOUK7cWYoUOWN3KEIIIYQhJCnLpKCgIM6eXc2LL6rnw4dD/U53+MB+4GvgMGq1LB5VKsMgn/8JRYrxz8qdpytZEubPT6ZYsYfQ9XlERUXZHZIQQgiRJZKUZUG+fDXQtH9391fPAyTe5uISgD9qw3/xtMcqbr8PLaPyQ2JR0AIMGi/buML48R2pU6cOy5YtszsYIYQQItMkKcuCuXOX8s47IQCcOQOfdwVudxLSD1UeYyeQglotuwQY1UkoCOLvha8PGDReNtK9u5OvvtJo3Hi73aEIIYQQmSZJWRYsXryYDz6IR9ehQAF47mXUCtjtVEWtpO0HygEFgJVAqnEx/bgWOvcCb7qbFxgIAwfqOBwDOXt2OCNGjMDpNGoJUgghhLCGJGVZ8OabbxIT0+2fshgfjoDcpe/wgVKoshjbUTXNGqM2/+81MKhUmD0T5i0ycMxs5M8/3+fdd0ewe/cuu0MRQgghMkSSsiwICAhA00r+8zwlBco7USUvbsUH6AHclfa8MqqR+QqMK80fpJqW7y9v0HjZzIABsGOHk8qV3wdS0D29D5UQQgiPIUlZFly8eJHnn9/A8uXq+bx5sOZl4OQdPlQMVa8MVJLWEDgOHDEwsAAYuxm+nQVbthg4bjZRujTABKZObUynTh2lLZMQQohsQZKyLPD39+err5axLW2zfoMGMPg9VC2yO9kMzEv775qoW5orDA4uFQb1gYGDDB43G7l0aS0XL64lOVmSMiGEEO5PkrIsCA4OJi5uEwMGqOe5c8PQJ1D9Lu/kDKp47FXUqcz6qM3/0QYG5wDnPVD0UQPHzGZ694bly2PImbMfKSlXZcVMCCGEW5OkLIscjlI3PA9IgLCj6XyoKurE5c6053VQtzSN3pxfFP44AtP3wPnzBo+dTfj4AEyjb98IOnbsQEpKit0hCSGEELckSVkWTZz4N4MHB/3z/Icf4NKPqIr9t1MIyIs6hQkQhDqJuQ9V8d9gD/aFRo0hIcH4sbOLjh1P0KnTGXx9pVSGEEII9yRJWRbt3r2bxYvh2iG/3r2h79eo6v23o6FWyw6jCsgC1AVCgYUYdxIzzdWyEFof/PyMHTc76dEDnn9+N3APBw7s5PLly3aHJIQQQtxAkrIsGjlyJNu2df6nVlnx4nB3S9TJyjupgmq9dG2bkz/QDDgG7DE4yDIQVRJme2G1/5slJMymZcvaPPxwT7tDEUIIIW4gSZkhSt749BiQXu3SfEBf1K3Ma2oCeVCrZQZW+b/mkc+hbQdIvF1/Ti8QGAiff57AyJHnuPM9ZiGEEMJakpRl0dmzZ+nRYxHz5//72sQx4LPAxQHi+He1zAdoharyv9XIKJXzsbBuKxw/bvzY2UnXrlClymrgbiZP/kNOZQohhHALkpRlUXBwMOvWneDcuX9f+/BDaPWuCx++CnwCrL/utYpAYWAxqnG5kcpBbD84Kn/qAOzZM59u3R7gk08+sjsUIYQQQpKyrAoJCeHAgeX06PHvayVKQMNyLnw4GAgHtvHv5n4H0BqIBcxoKu4D/afCd9+Dt/fsLl8e5s6FF17Yi+GnK4QQQogMkqTMECW5/rilrsOZpbjWaLwKcA44fd1rZVDNy5cBJpSxOLAWnuwPs2YZP3Z207o1+PlN5MqVJ/jss8+kV6YQQgjbSFJmgB9/HEfbtr7/PNc0mD2ef+uQ3Ukl1J/Ctpteb426vbnaqCivUx58+kJ4pAljZ1O//jqGYcOGEhVlxvKkEEIIkT5Jygyg6zqpqcE3FGeNWgdF+7jw4RCgNLCDG++gFUUlbKuBK8bFCoAGzhLQbxqcPvtvjTVv9thjsHGjTp06C+0ORQghhJeyNCnTNK29pml7NE3br2na8Fu8X1zTtMWapm3SNG2rpmkdrYwvs/r168eCBX0IDPz3tTx5oE5RFwdoA/RBFZW9XksgGVhuRJT/tWkrlCoNv/5qzvjZiaZBtWoAL7JixVC+/vpru0MSQgjhZSxLyjRN8wG+Bjqg1oB6aJpW6abLXgH+0HW9JvAg8I1V8WVdxA3Pzp6FM3+hapalpyCQ+xav50PVLlsPxGQ1vluPn1QF8pU2YexsbMyYT/nyy7eJj5c6ZkIIIaxj5UpZXWC/rusHdV1PAn4D7r7pGh0IS/vvnMBJC+PLNF3Xad78Gz66rrJCYCBsnImqOeaKY8BfwM0nIpuj/pQWZz3O/3CAsz2878qBBC8yZgysWHGGoKC5docihBDCi1iZlBXlxnWj42mvXe8NoJemaceBWcCgWw2kadrjmqZFaZoWdfasq1mPeTRNo3DhkoSF/ftaaCgcOQXUcnGQq6jN/je3QgoD6qGKyZ6++UPGWLQHWnaHOXPMGT+78feHfPmcpKY+wPDhD7F4sRkZsRBCCHEjd9vo3wP4Sdf1cKAjMF7TtP/EqOv6aF3XI3Vdj8yfP7/lQd7Kr79O4fHHQ254rUAolMzl4gBlgCBuXcm/ERAIuNolIKMcsHQZLF9r0vjZVFxcEjNm/M7cuRPsDkUIIYQXsDIpOwEUu+55eNpr1+sH/AGg6/pqVCqSz5LoDHFjxdg1ayBhHHDZhY/6onba7QFu7k0ZDDQF9mF8s/K0uVMfg7037/DzcqGhsGpVKu++Oxe1sCuEEEKYx8qkbD1QTtO0Upqm+aM28k+76ZqjqO6PaJpWEZWU2X9/0gWTJ0+maNE9nLxuF5yugx6L6yUtqqFOW+6+xXv1gPzA7LRrjOYHk3bCl9Nhqwl9N7OrsDDQtBMcPdqK7t3v4eLFi3aHJIQQwkNZlpTpup4CPA3MBXahTlnu0DRthKZpXdIuew7or2naFuBX4BE9m5RYL1q0KG3alLuhdVGDBjBxPqqXpSuKoWqW3epPxQd1Q/cippXIwAnP9oEhz5g0fjZ24MBeliyZyd69O+wORQghhIfSsknOc1uRkZG6+1Rh/xlVcOxfFxMgz/sGdlb8C9gJDADyGjXodY7DE61h1AMmjJ3NXbkCOXL0ACby36JyQgghRPo0Tdug6/ote+q420b/bC6C1NQbX/l9HPj/RMayshRuX5esLWr/2ewMjumqcPh+D2yKBinTdaMcOQB+ZeLEzrz00kt2hyOEEMLDSFJmoLZtX6RbtxtfCwmB/AXI2D6wn4FJt3kvFGgB7EfdBDaBMxVadYaHepozfna3Zs1MVq+eRFJSkt2hCCGE8CCSlBmoY8e7ads2+IbXevWCZz4F/DMwUATqXOr527xfB9UFYA5gRl6gQUx+8CkhfTFv5bPPYM6c/fj7z7A7FCGEEB5EkjIDPfPMMzz5ZLX/vB5ZhIzdaqya9nX7bd6/tun/ErAsIxFmQH1YUBBOGd0M3QP4+EBAgM7Vqw/x4IOtWbVqld0hCSGE8ACSlBksKakMyTfdqvzmRdRZUlflBEqgCsneLpkrAVQHVmFa0ZDYROjxAXz6qTnjZ3fx8Yls2rSEHTuW2h2KEEIIDyBJmYHWrVtHcPCvLFp04+v16kCBChkcrBrq9uWpO1zTBvBDNaQy6Tbj0pnw2Tcg26f+K29e2LLFSf/+E1DLlkIIIUTmSVJmoDJlyvDii10pUeLG14cOhbsezeBglVH9DQrd4ZocqFK7hwCzyme1A58nIVX+ptxSYCDATtaubUe/fo/idN7cUV4IIYRwjXyrNVDevHkZOfJVKtxiVaxeUf7bPulOAlHFZNMrhxWJKk47N4PjZyCOQ1fgzUVS6f9ONm5cw+LFf3Hq1J2WNoUQQojbk6TMYCkppbj5+7Kuw8huqMQpI+KA6cDhO1zjADqh+msuyeD4GfDBK9CkKcTcrn6al3vqKdi69RJFi8qJTCGEEJkjSZnB+vQZQMOGPje8pmkw7FkIrpLBwfxRJzA3p3NdOFALWMOd96BlQWo9KPYw5M5tzvieIEcOSE0dyBtvPMySJUvsDkcIIUQ2I0mZwR555BFee63sf14fMgSatMngYH5ARVRbpfSKz7YGgoHJqI4ARssHO/LCT5uldtmdxMU5+f33iUybNs7uUIQQQmQzkpQZrE2bNjzySLNbvlc1GMho3a9qqAKxe9O5LhjoApwGTKzQMPgjqNcAEs3Yv+YBQkNh1apUPv54BbfvlSWEEEL8lyRlBktNTeXo0fycO3fj65cvw8fdgIz2Ti+JOmXpyib78kANYAVwLIPzuOgycOKq7C27k9y5QdP2c+JEZ4YOfYaUFDOWLoUQQngaScoMFhsbS4kSb/PTTze+HhoKX3wLVMrggA5UW6VcuFaLrD0QhrqNaUZtsXIQfR8ck8oP6Vq8eCVjxnzDjh1m1SsRQgjhSSQpM1ju3Ln54Ye36NTpv+89/QRUrpyJQZsBHUi/PAaoUhpdgQvAgkzM5QIdeHwSjB5jzvieolcv2LcvmerVpQ2TEEKI9ElSZoJHH/0fFSv6/Of1xEQoFUvG95WByoROu3htKaAesA44mIm5XLB5Fjz5JOzcac74nqJQIYAh/PXXW6xfv97ucIQQQrgxScpMcP78JVauLPif1w8fhhkvkf6m/VtZC3yL63vHWwN5gSlAfCbmS09dCB0EBUqaMLaHSUhIZtiw13n//TfsDkUIIYQbk6TMBGPGjKFx45Ncvnzj6xER8O1EIKN9MEFt4gfY5uL1fsA9qJ35czIxX3p84VJueGE+xMaaML4HCQyEhQtTmTjxLOa0XRBCCOEJJCkzwf3338+8eV3x97/xdU2DJ3pA7jyZGDQ3UALYguvNx8OBJmmf2ZWJOV3w43goWgwOHTJnfE9RujQEBKwnPv4JJk6caHc4Qggh3JAkZSYoU6YMbdq0JiDgv++dPg1F9wEJmRi4OnAeOJGBzzRFNTWfTub2sqWnBATWgKBgE8b2QN9+O47evXuzVRqJCiGEuIkkZSZZvTqFLVv++/rOnbD9B+B4JgatBPiiWi+5yhd1GzMRmIHrq2yuCoPzrWDSEYPH9VCDB8PSpQ6qVbuc/sVCCCG8iiRlJnnggQ/46KP/vt6wIYxdAJTJxKCBwCOoTfwZURBoCexG3co0wUt/Qa++EG/GoQIP4usLTZo4gW5s3Tqf06ddPVIrhBDC00lSZpI//5zE22//9/5lYCDc1xgcmf2dD0etfmVUg7TPLsSU3piXT8Pvv8HGjcaP7YmuXDlFy5YdGTJkkN2hCCGEcBOSlJmkXr0GFC/+38bkAPt3Qt5VQGomB19Lxk9UOoDmqNOYZmxnKgUpQ0AvZsLYHihHDpgwIYXPPvOzOxQhhBBuQpIyk5w8eZIJEwK5dOm/723bBucXkvl+1ReA9WS8/lgZ1Kb/lWQ+IbyTIBgwE9auB93ovWseqH17KFToF3R9FHv27LE7HCGEEDaTpMwkmzZtonfvDdyq7eH998N3K1HFXTOjOuAEMtpSUQMao05w7s7k3OnYtgLq14VZs8wZ3xO9887T1KpVg0NSV0QIIbxaZnYnCRc0bdqUnTtHULbsa/95LygImmZmo/81hYH8qE37kRn8bCUgD7ACqIhr/TQzogwEdoaq9Q0e14M9+qiTkJAwSpQIszsUIYQQNpKVMpOEhoZSsWJz/G6zZWj7UgjIbKV9DbVadgx1KzMjHEAj4CRgxsKMDyTUhjdWmDC2hypcGJ55JgaH4xEuXYpFl3u/QgjhlSQpM9GcOceYNOnW7+3dC75HgORMDl4V1XopMycpqwM5gOWZnNsFY+dCtUjV71O45vjxGVStWpKvv/7a7lCEEELYQJIyE3311S+MHHnr3+IXXoCXfkf1qMyMnEAPoEAmPuuLKpFxiIx1B8iIANh3Ao5lpkiulypSBDp1iqV+fdlVIIQQ3kiSMhN9//33rF5d/ZbvaRo0CDdgkhgyd4qzNqoYrVm3GXNBwuOwR9ovuczhgG++0YmMfAOIJjXVjCOyQggh3JUkZSYqVKgQwcEVbvv+onGgTcvCBMnAN2QusQoE6qAalZ/NQgx34oDh8+H7cZCUZNIcHuk0H37YiLvv7iKJmRBCeBFJykx09uxZ3nvvNDt33vp9PQVy+mRhAj+gAqo0Rmb2ltVH3cpclYUY0nF+J/R/BP74w7w5PFFo6CFy5DhAYmKi3aEIIYSwiCRlJkpISODFFxexbt2t33/rLejxehYnqQ4kAHsz8dkQoBaqtEZsFuO4nVLgeBiqZrRfp5d74gn45ZfdBAXNszsUIYQQFpGkzETh4eFcurSYRx65/TVZ3ldWCnWSMrONxhsAOrA6i3HcjgappWHwHJBFH9dpmnpER/fmvvvac/y4nJgQQghPJ0mZiTRNIzS05h2vGf8SMD0Lk/igymMcBDKT9ORO+/wG4GoW4kjHstVQqBisMvFWqSeKjb3MqlUL2bp1g92hCCGEMJkkZSabPn0Zr756+yOIkdUhpFAWJ2kEPAMEZOHzycBtbrMaIi8kFABN+m9nSIUKcPBgCh07ym1MIYTwdJKUmWzVqlX88EMytztE98470KpXFifJgdofllkFgQhgLWDWKckASLgfppq1d82DBQUBfMOMGS+xaNEiu8MRQghhEknKTDZy5EhOnuyF4w6/0/WLkvnK/tecBn4EzmXy802AeGBjFuNIxyfLYPALcCGj7aG8XEoKDB/+Ph9+OMLuUIQQQphEkjKT+fr6AuVu+35KCnzcDViaxYmCUb0wN2fy88WAEqjyGJkpr+Gi5HPw1ccwLSv12byQry/MmpXK5MlXyXoGL4QQwh1JUmaypKQknn56GZMn3/p9X1948glwlMjiRKFAWdQpzMzWG20MXAK2ZTGWOykI+iDI19DEOTxU8eIQGLiexMT/MW+e7DETQghPI0mZyfz8/Jg5cyt79tz+mrdGQI2mBkxWA7iMOomZGWWBQqgOAWYWks8Fz86Fw8dA102cx0O99dZndOzYgUOHDtkdihBCCANJUmYyTdM4dGg/w4drd7yuZhiqCGxWlAeCgE2Z/LyGWi07D+zOYizp2L8dypWFSZPMnccTPf88TJ+eg1Kl/O0ORQghhIEkKbNEEGrT1q2dOAE/9CDzBWCv8QWaAqWzMEYlIA+wHFVU1iyFwFEfSlUxcQ4PFRYGHTpcAnpy9OghdFluFEIIjyBJmQUWLVpEhw5Xib1NOYiiReHN91DV+bOqAVA7C593oOqWRZP526AuzpPUEr68w21dcWe7dy+lUqXyfP3113aHIoQQwgCSlFkgKSmJU6ccnD17+2teewHKVTBowgRUk/LMqo46OLDcmHDu5OcV0Ol+OHLE/Lk8TfnyMGxYMl27FrA7FCGEEAaQpMwC7du3Z9Omjylb9vbXJCdDhXjURv2s2gT8CZzJ5Od9UStuh1FlNszkhLmzYc0ak+fxQJoGb7wB4eHPAGdJTpZSGUIIkZ1JUmaZqnd899gxmP4CsMuAqaqh/mQ3Z2GM2qitcCsMiOdOcoFzCMSXN3keD6br0fTtW5O+ffvaHYoQQogskKTMIi+//Cu9et3+BGapUvD9L4ARG99DUG2TtgDOTI4RANQD9qC6BZjJH15cCMvXqGK6ImM0DcqVO0FExHlSb9fPSwghhNuTpMwiAQHBBAXlvO37mgb9ekDpogZNWAOIA/ZnYYy6gB+w0oiA7uzUHmjaAMaONX8uT/TSS/Daa4txOMys/CuEEMJMkpRZ5LXXXmPMmHZ3vObiRSiyHzCiaXc51IrZ0SyMEQxEoir8xxgQ050UAUdnqN3G5Hk8WiJRUV156KEHZH+ZEEJkQ5KUWaraHd89exZWfEnWVreu8QEGAFlNchqgisqavVqmQWpteHWVyfN4uMOHD7NixSyOyHFWIYTIdiQps0hCQgLVqo3ms89uf03ZsjB3JVDToElD0r5mZZtRGOpW6CaMORmajllroUod2CZ34TKlWzfYvfsKZcvKb6AQQmQ3kpRZJDAwkOrVa1KkyO2v0TRo2xBK5DZw4sXAj1kcoxEqsbOibEUg7DsOR7Jy29XLBQdDamo/vvnmHc6fP293OEIIIVwkSZmFxo//m+7dQ+94zdGjELoc4/ZwhQDHURX6Mysvqv3SeiDeiKDuIBiSHofdRiamXmjv3hieeeYVfvzxB7tDEUII4SJJyiylkZpamTu1Krx6FXZNA04ZNGUV1P6yzVkcpzGQhErMzOaAEUvhx4mQkNUm7V6qQgWIitIZNkxqjAghRHYhSZmFZs6cSVhYFLt33/6a8uVh8yGgokGTBgPlga1AVr4/F0ad6FyDSs5Mdvkg9OsFP/1k/lyeqlo10LTXOXt2LgcOHLA7HCGEEOmQpMxCERERPPZYYwIDb3+NpkGVcCgWZuDENVG3HfdmcZzGwFXUpn+zhYPWG2p1smAuD5aamkLr1l3o2fNB9Dst0QohhLCdJGUWKleuHJ999galSt35uq1bgd8xbl9ZGaA5arUrK0oAxVHlMSy4K6aXgaHzpcp/Vjgc8OmnSYweXQhNu31HCSGEEPaTpMxiul6FixfvfI2vL8QfBy4ZNKkDlZQZsXm+CSqurQaM5YKVUVCkJKxebc18nqhlS6hWbQbwO1euXLE7HCGEELchSZnFHn10GLVq+dzxmkqVYPVW1MqUkfahellmRVnUitsKslb/zFW54HIoJElLxyz74ou+VKpUngsXLtgdihBCiFuQpMxi999/P889F3HHE5gAZfNAkTtXz8i4pcB8ICtbizTUatkFYIcRQaUjABIehEUWHC7wdI0bx9OxI/j7+9kdihBCiFuQpMxiHTt2ZODATqS3vWfJErjyKcb2nKwFnAOOZXGcCkB+YDnWrJYBHyyFV9+GS0bd0vVCtWrBqFEnyZHjO7tDEUIIcQuSlNkgJqY0p0/f+Zr8+aFYMYwtP1EF8Ac2ZHEcB+ok5hmyfqLTRQkn4a1X4O+/rZnPkx058hLt2zfi4MGDdocihBDiOpKUWUzXdYoXH8Y779z5usqV4a9pQEEDJ/cHqqJuO2a1Mn8VIBdqtcyKSgtFgQFQpoUFc3m8ZHbvXs++fVbcfxZCCOEqScospmkaX3/9BQ89lP5vffl8UPAONc0ypRaqoGxW93r7oFbLTgCHshqUiwrA4Dlwwchbul6oRAnYty+Zdu0W2B2KEEKI60hSZoOHH+5HvXoR6V43aRKcexWINXDyIsAzqJWnrKoBhALLDBjLRZujoGgxmD/fujk9kZ8fwJdMmvQq69db0TtLCCFEeixNyjRNa69p2h5N0/Zrmjb8Ntd01zRtp6ZpOzRN+8XK+KySkJDA+vVFuXz5ztdVqQJN7zN4cg31p+4k67cwfYGGwGGyfnjAVQWBypA3q4VwBVev6gwd+g6ffvqe3aEIIYTAwqRM0zQf4GugA1AJ6KFpWqWbrikHvAg00nW9MmpNx+OsXr2aunUXplsQtUIF+OoLIKfBATiBrwAj7l7VRt0OtWq1zA8SOsKEExbN58GCg2HRolTGjZP2S0II4Q6sXCmrC+zXdf2grutJwG/A3Tdd0x/4Wtf1GABd189YGJ9lateuzV9//Y9atdK/tmI+yJNgcAA+qHZJ28j66U5/oD6qMG10FsfKgC8Xw8OPw/nz1s3picqWBT+/yVy9OorNmzfbHY4QQng1K5Oyotx4k+s4/93ZFAFEaJq2UtO0NZqmtb/VQJqmPa5pWpSmaVFnz541KVzzhIWFce+9j5MvX/rXfvstXHgP41ouXVMLlZAZcQCvDhCAOolpkZRLMPEnVc9NZN2jjw6iXbvWxMXF2R2KEEJ4LXfb6O8LlEN1auwBjNE0LdfNF+m6PlrX9Uhd1yPz589vbYQGOXRIY86cgHSva9sWuv8PMLoIe3EgH1mvWQYQhFoH3QlYlSMXhNRnIKiaRfN5uNdeS+GXXwoSEmL0cV8hhBCusjIpOwEUu+55eNpr1zsOTNN1PVnX9UOo0qTlLIrPUl999TX33JOE03nn68qWhVeGohIfI2mo1bLjqCKwWVUflVKvMGAsV4XA0Lmwaw/ptq0Sd1apErRqtRP4gKtXr9odjhBCeCUrk7L1QDlN00ppmuYPPAhMu+maKahVMjRNy4e6nemRZccHDhzI2rVd0223BFA8EMLM2NheA3gYtWKWVSGoTf9bMbY1VDr2bIDKlWDKFOvm9GRTprxKyZLhHDhwwO5QhBDC61iWlOm6ngI8DcwFdgF/6Lq+Q9O0EZqmdUm7bC5wXtO0ncBi4Hld1z1yK3fp0qWpVq0ZDhf+BEaPhktjgHRKaGRYMFAa4/4WNEStwK0yaDxXFAP/VlAp0sI5PVjt2k5atXISEuJrdyhCCOF1ND2b3/eJjIzUo6Ki7A4jU6ZPf4uwsFdp1uzO1x09Cu9Pg29Oo24RGikRlf6WRDUaz6ppwBZgCBBmwHgueqwmjOmS/nXCVc8Cn9gdhBBCeBxN0zboun7LpQR32+jvVYYN+4nPP0//uuLF4fFuGJ+QgTpAsBtYa9B4jYFUrF0tA35YBC3aw8mT1s7rqS5c+JQePVqyadMmu0MRQgivIUmZjWbMmM3PP7t2ejT4MgRvNSEIB1AT1b8yq/0wAfIA1YAo4IoB47lI12HlGtghPbYNoeuwevUytmxJp8KxEEIIw0hSZqNy5cqRI0d1l66dMQOu/o3x9cpAJWUasNGg8ZqgugZY+f08DyQPggvSfskQefPC7t1OHnlEkjIhhLCKJGU2On36NB9/nIQrB90efhg+mIU5+7TCUIVHNqOSqazKB1QB1gFW1iL1hf8tgHkLITXVwnk9VGAgwARWrBjJihVW1joRQgjvJEmZjS5dusSwYcvS7YEJauWiRyMTg6kDlEVt/DdCEyAZWGPQeC46ugHatYZJk6yd11M5ndC//5u8+uoLdocihBAeT86926hMmTKcPTuffPnauHR99B7IvRBimqP6VxqpHMaW6S2Aaju/DlUqw+jit7dTGvy7Qe0WFs3n4Xx8YOpUJ4ULB9sdihBCeDxZKbORw+EgX75GuPrHcOAAxG3A3OKsJzGuHlpT1MqbUSc7XeGApCrwyjKp8m+UiAgIDV2A0/k1u3fvtjscIYTwWJKU2Wzx4jUMG+baRrF774W/N2BMBf5buQSMRvVeMEIhoDzqFmaCQWO66Lf5UL4KHPTIfhD2GDz4GRo3bsiFC0Yc0xVCCHEzScpstnnzZr7/Po6LF9O/1t8fWpaBAKNvXV5zbcP/RiDFoDGboRKydQaN56oQOH4ZThvR11MAMGhQCp9+mpfcuS2sCiyEEF5EkjKbPf3008TEvEiuXK5dv3YlBE0A4k0KqA6qvphRd6mKoBK91Rh3iMAVYRD/KGz3s3BOD1ehAvTuvR9N+4CUFKOydiGEENdIUmYzPz8/NK2ay9f7+kKIjvF9MK8pC+RCFX81SlNUEml1NywNXpoP334PyckWz+3Blix5jYiIkhw6dMjuUIQQwqNIUuYGRoxYyqefunZt48YwaxnqdKMZHEAkEI1xNcaKoRqfrwKSDBrTRed2wID+MGWKtfN6stKlnZQufZGUlKt2hyKEEB5FkjI3sH79QTZvdn2jWLWCUDAYMOt0YR1gKBBi4JjNUEneBgPHdEUZ8HkUKqfT9F24rnhxWLAgjnLlxtkdihBCeBRJytzA9OkzGDeuhsvXL1gAsW8D500KKCDtoWNMhX+AEkBJYCWqqKxVNHAWh2fnQYLFJ0A9XXz8Rwwd+gDbt2+3OxQhhPAIkpS5jaouX1m6NNRsiHkrZaBWtb5BtV4ySlPUIYJNBo7ponmLoVA4bN5s/dye6vJlnYkT/2L+/Bl2hyKEEB5BkjI3EB0dzd13r2fePNeuL10apv4JWn4TgwpG/e1Yj3HJXynU/rIVGFdyw1X5IaUoaGaVE/FCBQqopuXPPivF4IQQwgiSlLmB3Llzc+DAVZdqlV2TPwSqhmLc7cWbaagN/6eA4waO2QxVpHazQWO6Kgji7oXZsRbP6+Fy5wYYw549o9m0yYYlUCGE8CCSlLmBwMBAtm/fRvfuri/jLFkCW4cBR00LC6oB/hhX4R+gDFAUWI71q2XAW/Nh6HC4KgcHDZOaCl27DmDgwCfQpbeVEEJkmiRlbiME1cHbNbVrwyNDUDXFzBIAVAd2YFx5DA1oDsRi/WoZEHcMPvtQHZYQxnA4YMIEJ3/9VRBN0+wORwghsi1JytzE8uXLqVjxKK72ew4Nhe8+ghwFzY2LekBXVIJmlLJAOLAM61fLSoA+CArUsnheD1e7NhQuPANd/4XTp0/bHY4QQmRLkpS5iYIFC1KuXAmSMlBc1Qeonox5LZdANT+vCvgaOOa11bJL2HISk9wweDZER9swt4cbPvxRIiNrcjEjGySFEEIAkpS5jYiICKZNG0c11zsusWkTrHwd2GdaWEoSag/YYQPHLIM6iWnT3rL1K6B4CVi82Pq5PVm3bokMGJCD0NBQu0MRQohsR5Iyt1IFpzPI5atr1oRPx6BuB5rJAaxBtUkyyvWrZRsNHNdVxcC/AYSXtmFuD1anDrz44j58fMbIpn8hhMggScrcyFdfjSJPniSXK8/7+MAzj0HJIubGhS9QC7Uid9HAcUsDxVGrZVY3DPeDqy3hWxf38ImMiYp6loYNa8n+MiGEyABJytxItWrV6N+/BvEZ2CN2+TKUPARcMC0spXbaVyN7V15bLbuMPatlwBfzoHM3OHvWnvk9VXBwArGxe4iOPmF3KEIIkW1IUuZGmjZtykcfPZ9WkNM1ly7Bks+BPaaFpeQCIlDJk5F7wEqhVstWYP1qGeBMhDmzISrK+rk9WaVKsH17PDVqzLc7FCGEyDYkKXMzul6Hc+dcv75oUVi9ARz1zYvpH3WAwoCRhVc1oAVqtczIVThX5YeUIRBf0oa5PZzDAU7nq3zyyfMcOXLE7nCEEMLtZSop0zStZdrXwpom3QSN1KvXqzRqlLHf0vq1oH5xkwK6XlmgFxBm8LilgBLYtlpGAAydC0tXqur0wjgnTybz+uuf8PPPY+0ORQgh3F5mV8raa5oWDowCPjUwHq/Xs2dPnn++Ehk5uBYTAz6LMLfl0vUuAjEGj9kCuALYdBvxyGZo3hh++82e+T1VsWKweXMqr7ySgQJ8QgjhpTKblOUCXgD+B7h4VlC4omPHjjz22H1kpFtNYCBsmAKcNCuq66SgUnGj63uVTHusQNVFs1oJ8L0barWwYW4PV6YMaNoHnD49i4MHD9odjhBCuK3MJmUjgKm6ru8B5IaPwc6dK09GvncFBcHp05C7uVkRXccX1Q9zO2ofmJFaoHps2rFa5oCUmvDSMjK0Silc43Q6ad68K3379rE7FCGEcFuZSsp0XT+u6/oCTdOe0XV9uNFBebsmTV7jmWcy9pkcIdDKqkKodVGp+HqDxy2B2l+2EntWy4DJy6FiTdi1y575PZWPD3z+eTLfflvM7lCEEMJtZfX05VBDohA3+OijT3nxxYxVhI2Lg0PfAZtNCelGeVHlMaIwfmP+tdUyoxM+VwXD4bNwXMprGa5tW6hU6TdgPvEZKcYnhBBeIqtJWQZ2PglXderUiQYNmmboM8HB4BcPmlUrTPVRuwmNTl6Koyr9rwQSDR7bFSGQ2B+2Btswt1fQeeedbkRG1pLETAghbpLVpEx235jA6XSyenVB9mWg0bimweqV0KibeXHdoBTwLGpzvtFaoGqhrTVhbFdo8MZi+PoHuGpkTTYBQGTkJVq0cEpvTCGEuEm6SZmmaZc1Tbt0i8dlwOyui14pJSWF5s2/YcyYjH+2a3mMrbh/OxoQmvbfRs9XDHV7dCVg02LKlcPw9GMwfrw983uytm3hq6/2ERw80+5QhBDCrbiyUlYL6Kjretj1D6ADUMHc8LxTQEAAc+ZM57nnMl6X96/hwFTjY7qtP4A/TRi3Jer25SoTxnZFOPAIVGlv0/xeYO/ex+nW7S5iY2PtDkUIIdyCK0nZp8Ct/tWMRQrHmqZFi3YULFgtw5+7pzMUrmJCQLeTF9V30+iG6IWAKsAaVFFZO5SEQXPgsl3ze7gLFy6ybNl8du7caXcoQgjhFlxJygrqur7t5hfTXitpeEQCgJiYGMaNy8mxYxn73PPPwyP9zInpluqg/hatM2Hs5qhbo8tNGNtFm7ZBsdIwfbp9MXiq+vXh8OEkGjT4zz8vQgjhlVxJynLd4b0gg+IQNzlz5gyPPLKEBQsy/tmOpYCzhod0a2FAJWATxp+WzAfUQJXesOsOVx64WhQCctk0v4cLDgZdH8rvv3/BmTNn7A5HCCFs5UpSFqVpWv+bX9Q07TFgg/EhCYBy5cqxY8cU+mSiAPoXw8ExEevOxtZHJWSbTRi7WdrXpSaM7QofSO4Kvxp9e1b84/jxOPr0eYbPPpPdEEII7+brwjXPAJM1TevJv0lYJOAP3GNSXF7P4XBQqVJnIAcZ3VQ1eDBcrghzdKypJBeO+ptQ3oSxcwG1UcVkG6H2sNngp3WgLYaXB6lejsI4xYrBsmU6tWqFpn+xEEJ4sHRXynRdP63rekPgTeBw2uNNXdcb6Lp+ytzwvNvOnbt5+eU8Ga6V1bgxDO5F1qvQZUR1INCksZugfnxYYtL4rkiAcd/B1Gk2xuDB6tYFX983iYtby6lT8s+KEMI7ufxtW9f1xbquf5n2WGRmUELZv38/H3xwnD17Mv7ZCD8I3GF8THe0E5hlwrihQD1gG3DahPFdEQapA8GvoU3ze4HU1CSaNm1Bz549pLCsEMIrWbmWIjKoXbt2XLo0npo1M/7ZyZMg4U+s3SB/HnUK04zEqSEQANj540AovLoYNu4G6RBkPIcDhg+P57XXCqFp0sFNCOF9JClzYwEBAQQFNcnUZ/v0gQ+nAjmNjemOaqNuM5rRHikYlZjtAY6bML6LYs9CvRrwwQf2xeDJ7r8fmjX7A1iO0+m0OxwhhLCUJGVubtq0TQwYkPHu2PnzQ7824Gvln3AwUA3YCsSZMH79tDnsXC0Lg5QWUKWtjTF4vFRGj76Hhg3rkZhoR1d6IYSwhyRlbm737t3MnZu5xthnj0GhZVjbP7IBquCrGcVkA4DGwEHgkAnju6o+jNwCzlQbY/BwRYqcp0iRs8TLfWIhhBeRpMzNDRs2jAMHXiI444tlxMTA6eWAlYfZ8qMSs4ImjV8HtfF/EdbVYbuFLUegfkeYJqcxTXHXXTB58lFy5VpsdyhCCGEZScrcnMPhAOpm6rN16sC2Q0ApQ0NKXztUlX8z+KEKyh4D9po0h4txbNoO2zJxMla47syZx3jyyT5cunTJ7lCEEMJ0kpRlAyNHLmHw4Ix/zuGA8oWhRiHDQ0pfAmrDvxl7tWsCeYCFgF23EH3A2Q/2lLVpfi9x4MAFxo+fyOrVq+0ORQghTCdJWTZw8WIC585lrtr5sWNw+Xtgv7ExpeswMBtVu8xoPkBL4AywxYTxXeWA8Vvh2ylkqpacSF+DBnD0qJN27Q7bHYoQQphOkrJs4OOPP+aXX7pm6rP584NfApBkaEjpi0A1FF+BOXu/KgNFgMVAsgnjuyoJBvWCV1+zMQYPlzcvwHMsXTpRqv0LITyaJGXZRub2lQUGwq5tULyBweGkx4HqVXkaOGDC+BrQBriEOXXRXOUPzgegQl8bY/AC58/H0bHjw7z55ht2hyKEEKaRpCwb0HWd9u0n8sormR+jSwSqVIWVqqJOSq40afxSQDnUalwmSoYYpjh8tAH2n5dK/2bJmxdmzEjlo4/y2x2KEEKYRpKybEDTNEqUqEyBAj6Z+nxyMkweACw1Nq50+aIKvmqYd4uxNepQwXKTxndRfBLUawLPPGNvHJ6sRQsICXmPlJQ1nDlzxu5whBDCcL52ByBc89133wP7gGUZ/qyfHzzcEz4/asOCUgPUbUyzFARqoIrV1gNymTjXnTjgQiHwLW7T/F4jha5dW3PmTAVWr16Lj0/mflARQgh3JCtl2YiuNychIXOffedt6HqPsfG45NrfsFjgsklztECtxtldZ7QZTAmCy9IZyFSPPRbH0KG5JSETQngcScqyiZSUFEqV+o4338z8GG3DgfOGheS6ROBrzLt9mhO1SrYFa7sX3MLJS/DgCPj+e3vj8GRdu8KDDy4ApqPrNrZ1EEIIg0lSlk34+vrSp89jNGjgn+kxvnwatCnGxeSyAKAKsBm4YtIcjYFAYIFJ42fArEnwzfcg+YK5ZszoRWRkDWJjY+0ORQghDCFJWTby5ptv0aVLk0x//p23oMljBgaUEQ0xr1E5QBDQBFUk96BJc7hCA+4B3762tub0CnnzXiIo6Ji0YBJCeAxJyrKZixcbcPp05j7bti0M7mZsPC7LB1RAJWVm7bmqi7qVOR/72i8BBMH6U/DFCthpRkcDAahq/8uXx1CsmHSFF0J4BknKspGkpCSKFv2IDz7I/BiVfCFwg3ExZUgjVEJ2xKTx/VCb/qOBHSbNkQHPPwHtOkCS1d0UvIimwdWrw3j22T4cPXrU7nCEECJLJCnLRvz9/fnqqy956KGgTI+xcB4kzABijIvLZcWAoagWTGaphiqTsQjri+XeJKUhRDwK/pnfBihccOpUAj/+OIG5c2faHYoQQmSJJGXZTN++j1G7dvNMf753b5i4HMhtWEgZc62vulm3MB2ogrIxQJRJc7iqmMoN5x2QTf9mKl0a9u9PpX//Q3aHIoQQWSJJWTaTmprK+vXl2LMnc5/PmRMeaAAFQoyNK0PmAGMwb99XWVQLpqWAG7Q9euglaNYcnE67I/Fc+fMDfMyOHT+ydetWu8MRQohMkaQsm0lJSaF589F89VXmx4g+CblnAScMCytjigHnMG/flwa0RSVkGW+AYLjzqXAyCS6bVTxXAJCSksrddz/O4MED7A5FCCEyRdosZTP+/v7MmDGDSpXuBTJXCiAsDM5tB/IARY2MzkUVgfyohKky5vxoUBioCawFIoG8JszhqspwqAocuAq1c9kYh4fz9YXffnNSrJidf9hCCJF5lq6UaZrWXtO0PZqm7dc0bfgdrrtP0zRd07RIK+PLLlq0aEXBgi0y/fmwMDgdDSUzX/IsaxxAU+AsYGbJiJaoHzvmmziHKzR1p/aRifD+B7K/zEyRkVCw4DR0fSxHjph1zFcIIcxhWVKmaZoPqtlOB6AS0EPTtEq3uC4UGIJa4xC3kJSUxIQJoazNwu+Qjw88WBmwa59TZVTtMjP/lENRBWV3Y29B2TTbl8DLr0jtMiu88soT1KpVg9OZLeonhBA2sHKlrC6wX9f1g7quJwG/AXff4rqRwPtAJltvez6Hw8HAgVOYODHzY+g6LH8HmGFYWBnjALoBD5k8T31UQdm52FtQFiASfAdCUBGb4/ACffok8/LLIeTPn9PuUIQQwmVWJmVFgWPXPT/OTTuaNE2rBRTTdf2OBYc0TXtc07QoTdOizp49a3ykbs7X15dNmzbz2Wf5Mj2GpkHb5lC4nHFxZVghVHskHfN6EvkBbYDTwCaT5nCVAxJzwRMz4MABm2PxcBERMHToCRyOV0hJsblgnRBCuMhtTl9qmuYAPgGeS+9aXddH67oeqet6ZH51Ft7rlC5dBocj8/vKAF57DZ5+2qCAMusiMArIZIkPl1RGnfhchFusvy6YDeUiYNEiuyPxfDt2fEyFCsVZs2aN3aEIIUS6rEzKTqC+NV4Tzo1FGUKBKsASTdMOo248TZPN/rd29epVXn31CvPmZW2c7hVRq0h2CQWSgSWYt1qmAe2BOGC5SXNkRGkIbAml/rOjUhgtPBxKljyPr680LRdCuD8rk7L1QDlN00ppmuYPPAj800lY1/VYXdfz6bpeUtf1ksAaoIuu63bXZXdLgYGBjB69jqwuAIz5ELTvsW8FyQe1Gf8UsNfEeYoC1VF/q+xoMXU9P4hvDK+usjkOL5AzJyxYkERk5KeYl/ULIYQxLEvKdF1PAZ5GbbneBfyh6/oOTdNGaJrWxao4PIXD4eDo0RO89lrhLI3Tpw888ib2Vqyrhmr7tBRzv2+2Qv2Nt7tERpqJy6FaPdhk9143L+B0zuHttzsxfvx4u0MRQojbsnRPma7rs3Rdj9B1vYyu62+nvfaaruvTbnFtc1klu7OAgAAga/vKKlWCdweBj58xMWXKtdWyk4CZG+DDgEao2mjuUMLKH3YfhYPuEIsXmD9/DsuXT7Y7DCGEuC232egvMu7SpUt0776X33/P2jiBTog4hNpzZZfqwD1ASZPnaYhKzuZgf4mMYEh+HFYE2hyHF/DxgZkzdb77bgf2/kUXQojbk6QsGwsNDeXgwQRisrhH6uhR2PUTsM+IqDLJB5WYmX0b1R9oDUQDW0yeyxUO+HwtjPgadu+2OxjPFhICmraXkycf48cff7Q7HCGE+A9JyrIxTdOIitrGk0+WyNI4VavCqijwr2VQYFmxBfgdc/eWVUFt/F8IJJo4j4v0eHjzf/DRx3ZH4h0++eQ3hgwZKNX+hRBuR5Iyj9Aiy/0UG9SGThHGRJMlSahjIGa2RXKgSmRcwT1KZARBah/Iea/dgXiHkSMhKsqHggWv2B2KEELcQJKybC4mJobq1RcyenTWxklNBedcYIMhYWVeTdSeryWYu1pWDHW7dBVwzsR5XFUQPlsPS/bByZN2B+PZgoKgfPk44AGWL1+E02lXA1ghhLiRJGXZXK5cuahYsQYFCmRtHIcDLuwHvwvGxJVpvkBjVEOu/SbP1QbVhmkWblHCKjUVOrSB7g+Q5ZVPkb4NGzbQtGkrRo0aZXcoQggBSFKW7Wmaxm+/TeOee8pmeawli+GhYQYElVW1UHXLFmDuCckcQEvUrdKdJs7jKg0S6kOxu1VvUmGu2rVh/Hjo1887W7UJIdyPJGUeIjm5OVevZm0MHx/oVxOwu3+zL9ARtZJldnISCRRElTROMnkuV1SEP6/C+hOyWmaFXr0gMPBJEhP3cuGC3cvEQghvJ0mZBzh37hx58kzg+++zPtaW6eDzBaofpZ3KAWUxPynzAToBl4BlJs/lIqcOdz8H7TuoW5rCXLoeQ+vWkTzwQHd0yYSFEDaSpMwD5MuXjyFDniTSgNbtNWpAw/bYn5QBOIF5mH/4oDjutekfiI6HQ5chTuqcmk7T4IknLjNoUA40uW8shLCRlt1/MoyMjNSjoqQbk1IZIzZHnb8KRT+BRLsPpenAT6hEaTAQYOJcV4AvUfXLemP+Cl16dPBxwLr+UCtr7U2FyzRgOikp7fD1tbMZrBDCk2matkHX9Vsuo8hKmQc5fDiSU6eyPk7eYGgZDMRmfaws0VDV9+OAtSbP5Yab/p069B4Pb7wptzGtoTN1ag+qVKnA2bNn7Q5GCOGFJCnzEDExMZQpMwEjTvfHxMD84cCarI+VZcWA8sBKIIsHGdJ1/aZ/N6j0D7BzJbz1Nmxxh5ZQXqBEicsUL36e5OR4u0MRQnghSco8RO7cufn55x95+OEQA8aCv/+C0p0NCMwIrVAnI82uvn/9pn93qPQPUBO0gZBa0O5AvEONGjBv3kWKFJHaZUII60lS5kF69uxD6dKdDBmrc2d4sqkhQ2VdAaADqtq/2YoDNVCb/t3hDpYGKbmg92RYsRpS7C5X4iXi4t6lT5/WLFmyxO5QhBBeRJIyD6LrOgsWlGHxYmPGi7gMDoPGyrK6qOTMCq1Rlf5n4xaV/gF2bYcmjeCbb+yOxDs4nbB+/WK2bnWX/wGEEN5AkjIPomkazz47lffeM+bo4K7N4L8FcJftNZeAPzF/BcvdNv0DFAI6Q+kWdgfiHcLCYOPGVAYPnonbbDAUQng8Sco8zKRJfzNlijHfuQcPhilrgCBDhss6H2AfsMiCuSJRidAcIMGC+VxRCwYsgHNXIFHyBNMFBgJsYOPGh/jwww/tDkcI4QUkKfMw5cuXJyjoXkPGCg6GthFQOhfm9qB0VQjQENgFHDd5Lh/gLuAyqgenmzh2ASpGwvDhdkfiPX7++W+++up9YmPtrhEjhPB0kpR5oIkTUxkxwpixTp6E+K9xn9t4DYBgVKJk9n6vcKAeEAUcNXkuV/nCuQKQWsTuQLzH++/Dhg1x5MzpLn8JhBCeSpIyD7R69R6mTAk2pKF14cJQtSz4mFlNPyMCgGbAYeCABfO1BHIC07C/Ufs1rWCCDicv2x2IdwgIgHz5EtD1e/n++y+5etXsgnlCCG8lSZkH+vjjj9m48UWMaOPncMDcWdC1S9bHMkxtoClgxWpRAOo25jlghQXzuehCPHR6Re37E9bYuHE/jz8+mJ9++snuUIQQHkqSMg8UEBAAdDV0zL5VgdOGDpl5vqgVrGCL5isHVAGWAWcsmtMFm7fC9EXStNwqtWvD6tXw1FPyGy6EMIckZR5qwoTN1KkTgNOgpuK/vwuO8bjPLTyAaOBnzG+/BNAetWo2Hfc49ADQBE4/ACfd5XSoF6hXDzTtJU6fnsq2bdvsDkcI4WEkKfNQOXLkoHDholy4YMx4gwZBz1dxr78xDuAQYEV9zxxAO+AYsMGC+VzhA/Gp8NBv8NXXGLKHUKRP11O49977eeCB+0iVTvFCCANpejb/lzwyMlKPioqyOww3tRy1+coY0Zeh+GeQ4k7fh2aiTkc+iWombiYdGI8qx/E0EGbyfK7aDEyBFSugUSObY/ESGzaAj08NatRYh2r/IIQQrtE0bYOu65G3es+d1j2E4Rpx+XJew1ZQglOh/F4gxpjxDNECdVtxDuaXyNBQm/5TUcmgu/w8Ux0cT4Az3O5AvEft2lCjxmZgMLt27bI7HCGEh5CkzIPNmjWHfPli2LHDmPEuX4Zdv2BNKQpXBaMSs0Ooav9my5M23x5UEVt3oEFqYej5N6zaAFeu2B2Q95g7dxSVK1dm6tSpdocihPAAkpR5sFq1ajFoUBdy5DBmvPBw2LsXItoZM55hIlErWKUtmq8+qgXTLNynLyhw/CQ0bQQvv2x3JN6jZUt45x0HbdvmtDsUIYQHkKTMgxUqVIiPPvqVkiVDDBuzTBl4rgHucwIRVEukSFSpDCvi8gG6AHG4VQsmcoCzExRqb3cg3sPPD4YPdxIU9BCJiYeJiXGne/tCiOxGkjIPp+sBbN9en/PnjRszaS34jgEMKrdhmEPAN6h+lWYrglox25A2r7uoBiM2wpZTIK0araPr0XTsWJ177ukqJzKFEJkmSZmH279/P1WrLuS334wbs3g4VKwEJBo3piHCgAvAQovma4HaYzYFcKNaYQkp0Ow+aNUakpLsjsY7aBo89tglBg70xeGQf1aFEJkj/3p4uHLlyjF+/Lfcd59xf9RdusDiaRDsbtto8qJWrzYDJyyYzx+4B7gEzLNgvgyILQY+VcDX1+5IvEePHnD//YuAb4mTNgtCiEyQpMwL9Or1JIUKGVevDCBvMNwfjuoJ6U6aAiFYUyIDoBjQCNgI7LVgPldVhHUl4c+dUlTWasuWDaJkyaKsW7fO7lCEENmMJGVeICUlhb//LsXq1caN6XTCnOHAfOPGNEQg0ApVed+q0h3NgQLANKxp+ZQB/b6FGpFw6pTdkXiPypWdNG+eQJEimt2hCCGyGUnKvIDD4WDAgOl8951xY/r4wLix0GGQcWMapgbQCyhj0Xy+qNuYV1FlMtxIHLDvNJxyo0bqni5vXvjzz0TCwweg61dJTk62OyQhRDYhSZkXcDgcLF++iu+/r2bouO3awRtdDR3SGA6gLKoCv1V1xAoDzYDtaQ93UQDi+8HvkpRZTtej6N+/Cn369CG7t7MTQlhDkjIvUa5cOXx9uxo+boFkKDgfa8pQZNQB4FNUr0orNAaKolowudPvhwPeWw69n4ElS+wOxntoGpQte4hy5az6CyiEyO4kKfMiX3wRz4cfGjum0wmXtwAnjR3XEEVRfTGnY01NNR+gK5CcNqc7LY4kwy+/wDgDS6OI9A0fDm++uRxN+wWn090K+wkh3I0kZV5k7doTLF5sXHV/UBX+z5yCio0MHdYYgUBH4DRg4CGHO8qPOmiwF1Waw10EQGpf2FITElPsDsb77N3bl2rVyrJ27Vq7QxFCuDFJyrzI2LFjmTXrNcPHDQmBoQ1wqwKq/6gIlAeWAFZ1wKkHlABmAxctmtMVwbDpFPQdD99+a3cw3iV37mRy5jyOpskxWCHE7UlS5kX8/f2Bh0lJ8TF87D2/gc93uF/rJVCrZT7Afovmc6BuYwJMxb36hAK/joVnn4MTVhTYFQDkzw8rV6ZQt+4rwGXZ+C+EuCVJyrzMrFkbKVzYh6NHjR23bRtocR9ul4AAkBMYDNSxcM7cQDtUX0x3qyHaDHgCzhqfm4s70DSA7Xz1VX3uvfce2WMmhPgPScq8TKVKlWjduq7hPRHbtIHfv4CQYGPHNcy1rXTHsK7Aay2gHKrArjvdtfKFxFzQ7Q/4azpcuWJ3QN5F03YC20mSxqRCiJtIUuZlSpYsya+/LqJs2fyGj507ENpouOdJTFA9KscCCyyaT0PdxgwC/gTc7HvwgQNwf1d4/327I/EuAwfC338fIChonNzGFELcQJIyr+THyZNdOXTI2FETE2H5p6BZddIxo8JQDcs3AkcsmjMEuBc4j9r4707ygN4D/JrZHYj30TSIiXma9u3rMHfuXLvDEUK4CUnKvFBycjJVq/7B668bO25gICyYB/e/aOy4hmoO5ELVEbOqNERpVKP0TcBWi+Z0VTkYsQrm7Ia97tRQ3Qv4+DiJidlETMwmu0MRQrgJScq8kJ+fHz/88BOvvFLV8LFr1IC32qrDjm7JH+gEnANWWDhvM6AYMAO1auZGnDp0fQBatIJ4q9pSCcLCYPXqVB588AvgsNzKFEJIUuatunbtSkTEAHMGvwBhY4HD5gyfZeWAmqhq/1bxAe5D/R83CetW6VyU2ACCOqOaqwvL+PgARLNoUVPq1avN+fNulrELISwlSZkX27KlCu+9Z/x34fBwKF8M/Nz5B/+7gQYWz5krbd5oYKHFc6enMBwoCP2mwZEjIIs21vL3P4au7yUh4aLdoQghbCRJmRdbuHAdI0bonDK4XENQEKxeDgMfNHZcU+wAVlk4X0WgLqrtkxvu4fp1EZSrIBX/rda4MaxdG0fRos+j605SU92x4J8QwmySlHmx/v37Ex39F4UKmTP+/xpA4B7cs6DsNXtQdcSOWzhnG6AgMAVVpsOd5IPkehBs/HZDkQ6HA3R9Mk8/XYNBgwbJHjMhvJAkZV4sNDSUnDm7ACVNGX/VAkj4FevaG2VGByAUmIx1dcT8gPuBZOBv3CtpdQAtYNAK2HoKzp2zOyDvommQI8d2QkK22x2KEMIGkpR5uejoU7Rv78OMGcaPfc898Nc0yFPN+LENE4Qq8Hoe64rKAuRD9eQ8DCy3cF4XXUmCJg9CnXpwyd1W8zzce+/BBx8sR9MmkpLiZidChBCmkqTMy+XPn5+LF0OJizN+bIcD7u0MwxsD7nwnpjRQD9Wj0sp2SDWAqsAS4KCF87roUlmgBgS6a+ssD6V6ZOocOdKX6tXLMmfOHLtDEkJYRJIyL+fr68uaNZt44IFWps1R7CT4/Yi6XeeuWgMPACbtr7slDbgLtWo2Cbho4dyuKAaHq8Gg2dIf0w65c6dQqNBxcuU6YXcoQgiLSFImAND1vpw+bc7YBfJCyXxY1wg8M/xQJyPB2s33Aahk0An8jlsmrqMXQ4kIGD3a7ki8S1gYLFjgpH79/wHbuHz5st0hCSFMJkmZAOChh6bSsqXDlPpULVvC9nVQpqTxYxvuOPAFYOU+63yo/pjRqIr/7narNxRiCkNKAbsD8T7qVuYFfvutCeXKlWav9MISwqNJUiYA6NnzYYYNa4pZ5ZH8feGF2rhvlf9rCqPKVczA2hWz8qhWTFtQe9vciQP0LvDyHth5FlP2H4o7i4yMpU2beIoWdbeMXQhhJEnKBAB33XUXfft+kNb2xRzLvgPH71hXeiIzfIB7UG2QpmLtqlUzIAKYCxyxcF4XXUyAxgOhfEU4edLuaLxL2bIwfnwcISF3k5h4jDNnztgdkhDCBJKUiX9cvVqZiRPDuWrS3q/XXoUPxqOagruzfEBb4ACw3sJ5HajbmLmBP4BYC+d2UUxeuFoUfOREpk320Lt3NVq2bEZiYqLdwQghDCZJmfjHxo0b6dXrOJMnmzN+uXLw3H1QP9yc8Q1VBygLWH3qMBC18T8ZlZi5W5mqQhDTHnpOg8txIHmB9QYMuMizzyYREOBufzmEEFklSZn4R6NGjVi+fAEPPWTuju6ym4CZpk6RdRrQA2hpw9wFUAVtTwCzbJjfBQv3Q9l60K+f3ZF4n+bNoV+/g0BXdu/eylWzlraFEJaTpEz8Q9M0GjduhaY9a+o8RXNAqRy4V3uhW7m2v+4YKom0cn9ZJaAxsBGIsnBeV/nAmSJwuYTdgXiv2NgFNGlSh0GDBtodihDCIJKUif+YMCEP7dr5mHYS8913YeVUCAsyZ3zDHUPtLVtt8bwtgTKo1bJjFs/tikYwzR++WQ8xMXYH431y5oTPP0/i5ZdP4373uYUQmSFJmbiFYJzOcNO+0WoaFA6FwRHAbnPmMFQDoAKqN+ZRC+d1APcBOYFfATdNfJ7+HoqVhOnT7Y7E+zz0EJQuPRvowcSJPxMfH293SEKILJCkTPxHz549WbBgLXnzBpo6z/bx4Dsb9/8hX0Pt8coJ/Im1m/+DgYdQt3onAm74PVfPBwllQbeyRZW4webNk+jduw/fffeN3aEIIbLA0qRM07T2mqbt0TRtv6Zpw2/x/lBN03ZqmrZV07SFmqbJjhUbaJoGFOTixYc4dMi8eT7/HCbNA83XvDkMEwh0R7WKsrq4a37gQeACqhWTuyWx/uDsAo8uht1npYaZHWrUgCVLYNCglbhlry4hhEssS8o0TfMBvgY6oLYx99A0rdJNl20CInVdr4Zq0fyBVfGJG6WmplKnzkIGDDBvjuLF4e660LcGkGDePIYpDPQDmtswd0nUat1hYBru14oJOB8PdXtBtRoQHW13NN6naVPw8ZnMhQtdGTz4aTmVKUQ2ZOVKWV1gv67rB3VdTwJ+A+6+/gJd1xfrun7tX5I1QHaoaOWRHA4HH330Oe++28H0uUKWguN73G8F6FaKoP6vuYS1+8sAqgEtgK3AEovndtHlCHDWBqcUl7XNihWz+P77UWze7I7HdoUQd2LljaOi3HiG7DhQ7w7X9wNm3+oNTdMeBx4HKF68uFHxiZvcfffdQGlu88dgmM4d4HAcZKt94pOBU8CTqL1mVmmK2vC/FMgF1LRwblfkg4v5oPV4+LMD5AmEokXtDsq7dOkCBw86KVToY6A+uu6XtiVBCOHu3HKjv6ZpvYBI4MNbva/r+mhd1yN1XY/Mnz+/tcF5mXPnCvPcc6XYv9+8Odq0gamjoH5J8+YwXCfAidr4b+UKnwZ0RuXK04GDFs6dAXvOQb2W0Kkz6G54q9XTFSoEMI05c5rSoUM7rlyxujWFECIzrEzKTgDFrnsenvbaDTRNaw28DHTRdV2auNgsOTmZ7747xfLl5s6jafBkQdDm4pb7pf4jH9AFtd5rdWFZH9Shg7yojf+nLZzbVRrEt4LEtnDFnRvQe7iYmLWcP7+OlJTLdocihHCBlUnZeqCcpmmlNE3zR50nm3b9BZqm1QS+QyVkZyyMTdxG4cKFOX78JH37NjZ9rvN7IWQf1vebzKwqqNuJm1CV960UCPQE/IBfAHf8nlsCdgdBp19g6ky4dMnugLxPjx6wenUsuXI9QErKec6ePWt3SEKIO7AsKdN1PQV4GpgL7AL+0HV9h6ZpIzRN65J22YdADuBPTdM2a5o27TbDCQvlypULeBGz74AMGgR7dkOhwubOY6jmQBugsg1z50LVMLuKSszc9ATr8u1w7z3wyqt2R+KdfH0BljN0aAXq1YskNjbW7pCEELdh6Z4yXddn6boeoet6GV3X30577TVd16el/XdrXdcL6rpeI+3R5c4jCquMGnWUEiV8TG2n4+cHRfLCh625xY1tN+UAGqFWrpIAqxciigD3ow4d/JoWg7sJg9QecKAmpLh7v1MP1rv3OZ54Io6cOd20NYQQwj03+gv306BBA3r3bmNaP8zrHZ4O2o/AOfPnMtQUYBxg9UJEBHAvcAT3LC4LUBpmHYEHf4c3R4CU0LJenTrwwgvngcbs2zeLlStX2h2SEOImmp7Nj0ZFRkbqUVFSj8caTlQGYO6Rv5gY+PkPGHERLrjpLblbOgP8gLqt+CgQYPH8G1G7NCugVs98LJ7fFQeB8fDHn3D/fXYH473at/dj16687Nt3BH9/f7vDEcKraJq2Qdf1yFu9JytlIgN82LKlB3PnmjtL7tww5An44W4gO52/LYBKhs6g+lFYfauuFtAe1eR9ig3zu6I0MBB+TIB46QZkm59/TmbatFj8/ZfZHYoQ4jqSlIkMeeaZ5QwZ4mfJbcyIVAj8Bthp/lyGKQt0BPYBi22Yvz7QEtiG9aU6XJUP5uyHZh9A6zZwLrvdpvYABQpA9erxQCdGjXqcr7/+2u6QhBBIUiYy6IcfxrJy5Y84LPibU7Ys3NsZSpQxfy5D1UGdyLSr2n5ToDGwAZiHeyZmwPrdsGILHJfiN7bR9STmzh3D3LmjSLXiJy0hxB3JnjKRKbrejri4eeTIYf5cUSeh4Q+QnB2/Z+jABVShV6vnnQ2sA5qhema6oxSoXQzm9Qa/FAgNtTsg7+N0QmIiBAcPITb2dXx9/QkJCbE7LCE8luwpE4br2dOPrl01S1ro1CgA9XcBq8yfy3BLgFFY37xcQ+0vq4Hqk+muB+18YUM0VOkNVarBqVN2B+R9fHwgOBh0/XO6dy9PmzatZNVMCJtY2ZBceJDWre/l8uU4dH0JZvc69vGBvAlQNDX7lC/7RyRqf9dEoA+qrphVHKhWUMnAfNT/7fUsnD8DovNAznBIsvrEqviHpsGgQWe5eDEYh+M0kJ2qOAvhGeT2pciCS6gSGeY3X0xJgeg4qDYKLmanMhmg6pb9iCrs2hd1StNKKajG6XuAVkATi+fPgNK54c+7gItQq5bd0XizYixe/BpJScVo166d3cEI4VHk9qUwSRizZ/dg9GjzZ/L1hWI54c1qwGrz5zNUTuBhVN2wX7C+uKsvqoF5FWBh2sNNfxY7GAMNu0HLVtIr0066fowRI57khRcG4HQ67Q5HCK8hSZnIkh9/PM633wZj1b/b0cvAfyXZp2n5NXlRiVln7Nk04IOq+l8LWA7MwW0Ts8RWkHA3LJf9ZbbRNJg2zcmMGYfx8RmN0+mU5EwIC8jtS5ElFy5cIEeOnfj7N8WK7/KJiXDoOHSaqVZVsq2dQDHA6tOGOiohW4sq2dEZt/3RzEeDAbkh7DSMHInpexfF7Q0eXJ19+woxbdp0/Pz87A5HiGxNbl8K0+TJkwd//8YkJ/fkqAUnDAMCoEIZ+O0+8NuLe/Z5TE8cquL+eMDqHpDXTmU2BTYBf6O6Z7khpw5f/gJjfoO4OLuj8W5VqmyhWrWd+Pll55+EhHB/kpQJQ9x990k6d3ZYdhvT/xwk/4Kqw5XdhAAPAueBCYDVBxc0VNX/1sB24A/UCU131A7O3A9PzYf4RNlnZpfHH4f33z8G1Gbfvkn88ssvdockhEeS25fCEHPmzCEp6Tc6dx5n2W2mOXNgmQ+8mx3rl4E6Dfk7EA70xPoG5qCS2lmonpQPAm7cm7roMsh5GtavU3W1hD2eeMLBX38FsW/fMXLnzm13OEJkO3e6fSlJmTBQMuqI317LZtR16DoOpi0DSlg2rXG2A38Bd6MKvdphEzANlRz2ANw14TkIBWMgagKEh9kdjPdKToZ9+6BSpceAr7h0KZGwMPkDEcJVsqdMWMSP775ry8svWzejpkGOheDzGxBv3byGqQI8xb8JmR0/I9UEugEnge9Rt1XdUWk4XRvqjoFfF8OECXYH5J38/KBSJYDv+fHHilSsGMHBgwftDksIjyBJmTDUjh0aGzcWwMouLR99AL/8AYXzWzenoa4Vk40GRqN6ZVqtMqrjQAIqMTtsQwwuir4CvZ6HAYMhNtbuaLxbZOQhOnaMpUSJQ3aHIoRHkNuXwlDx8fEEBFzG4aiJWnqxzvoT0ORtSMwDZMdT+ydQG/99UDXNrK78DyohnAjEoG6pVrchBlc4gRjo3xa+7ACJV0HuoNnJh0uXhjN8+HlGjnyLvHnz2h2QEG5Lbl8KywQFBeFwFODixVF88ok1DcuvKaxD6g+oJuDZUVFUGyYNGAsctyGGPMBjQHFgMrAY9ywy6wPkgzEboWJfqFgZjh2zOyhv5mTVqrcZN240e/YstzsYIbItScqEKX777QT/+x9s3mzdnOHhMO4n+N8L1s1puALAo0AgMA51S9NqQUAv1D63pahaZm5cD+5QMFwMhyNuWm/NW7RvD4cPp9KwYT9gCsuXLyc52V1rrQjhnuT2pTBFamoqu3Zto3LlgcBKS+fWdXjwT/hjMVDI0qmNcxlYAbTBnrZMoFbIVqB6ZRYHHkDVWHNT/j7wXiO4sgKGD1cb0oU9Dh2CiAgHw4c/z8iR79kdjhBuRW5fCss5HA4qV64O/ML69aGWVmTXNCi5CRw/oPZGZUehQAdUQhYHbLMhBg1ogjqZeQJ1AOCsDXG4KMkJQz+DN0bA1u12R+PdSpWC335LZejQqcA2EhMT7Q5JiGxBkjJhquPHHTRuHM/IkdbO+/xz8MUXUL+qtfOaYgWqltkc7GmJVAV4BEhCnQ61I0F0VS1IfRL6rYXd52CvdSXzxE3uuw9y595Namod7r67BgMHDrQ7JCHcniRlwlTh4eGMHz+RF1/sa+m8+fLBwCdhbi+o5odbl3hIV2ugHrAG+AV76rEVA55A3Q7+C5iJ++4zywdbTkPN16FiJRg71u6AvJuuJ1Kv3m6qVFmKWnIVQtyOJGXCdN27dydnzm9wOqtw3uLCpGEBkH8J+E/DfZOI9PigbmV2Bg4BY7CnllkYasWsAbAe+BG4aEMcLkrIA6nNYYoGMfFYehJY/MvHB958E556agdQmTlznqVPnz5ckkamQvyHJGXCIoE8+GAROnZ0kGJxcjRxPCyaC7WKWTuv4WqjCrwGo05n2sEHaIfa9H8eGIWVXbUyxhdoAtOOQLVvoH4LGDXK7qC8XSx79nzG5s1/ExBwyu5ghHA7kpQJy/To8QSPP94LHx9r5y1YEBpFwvzeEL4POGrt/IYqAfRDJWYpqP1ddqwAVQQeB3KhbqkuxJ79bi46fh7WnYJZB8FpYbcJ8V9DhkBU1BUCAiJJSvqS4cNf4OxZNz5BIoSFpCSGsMEDJCb+QUCAtbMmJEDV6nAuJ1zsZO3cpogCZqDqid2FPaUzkoHZwEagJOqkZg4b4nBF2j91jYpDvxDIHQBdu9oakddbtgxat9b466/v6Ny5v93hCGGJO5XEkKRMWG79+iXcd19bJk1Kpm5da+c+fRqcAdD+V9h2xtq5DZeKKu66FAhHJUS5bIplE2rzfwBq71sFm+Jwkc9EKOKAg9vB1646cAJQnRiKFQsG3mHmzNJUqlSFUqVK2R2WEKaROmXCrYSHl6dy5SYUKFDS8rkLFoQiuWDafRA6mex9K9MBtADuB86g9nfttymWmkB/1CrZb8AUVHNzN+V8EI51gpbjYetx+OMPOQhgl2LFAK6SnPwMTzxxH0OGPGx3SELYRn5GFJYrXLgws2cvBI6g6w24fDna8mbSfslQMB6KBsBua6c2XmWgMDANyGljHAVRidlSVG21g6im5mVsjOl2fIAwWH4Uaj8FKTOgeEmob/HKrfiXnx+sWZNMSsoKoC/nz7/AqlX7uOuuu9A0ze7whLCErJQJG5XgjTfupm5dBzEWV94vWhR2bIMtX0LvasA53LPxtqvyoMpV5E97vhA4aUMcvkAr1GEEP2A86rZmkg2xuCilFtAbBm6Gzadg61ZIlcMAtggPh5IlAX7iq69qcM89XTl0yK7lXyGsJ0mZsFXr1g/RsWN3cuYMtnxuf3/VL3FkLfD/AVhieQjmuAJsRrVFWoXae2a1cOBJoD6qptm3wBEb4nCFAygDG6Mh8mOoXRdefNnuoMRLLyWyYEEqpUt3A5Yza9Ys4qzs1yaEDWSjv3ATc7lw4S5CQlIsP5Wp6/Dll+BfFZ5dBQnZtcjs9a4CU4E9qNuHXVH9NO1wGLXH7CLQELUPzl2bhevAFihTG0Y9BNVD1W21XLnsDsy7nTgBJUtqPP/8QN5550u7wxEiS+T0pXB7iYmJREaWo3r1Y0yYYF8cq49B675wtTxQwL44DKGjymbMRZ3KHIB9a+OJwPy0eHKhOhSUtymWDAhfCM6DcOgAlv+wIG60fDmULx9EgQLPsmtXVw4cOE2nTp1kv5nIduT0pXB7AQEBPPnkCzz22GBb4yjhAyHbocAxW8MwhgbUQRV5bY/6v90J2NHdJgBVS60PapXsV2Ai9rSLyoDjleBsXXhjOVxJkgbndmrSBAoUiAfe4fPPm9Cz5/3Exmb3ujZC3EhWyoQbeoM9e94kIgLs+CH4zBkICIMHJsHczagyD57yw/gaYBFqM34d7PmxzAmsRe3hcwKN0x7uekszTd4LcOFL+Pln6NXL7mi8W3IybN8ONWuWAN5k5Mgj3HdfNypVqmR3aEKkS1bKRLaybl1HqlRx8NVX9sxfoADkDITfOkPOiaiK9Z4iAiiG+jX9AETbEIMPam/Z06h2TUuBr1H739zY+Rygt4TPL8La43DoEFy+bHdU3snPD2rWBDjC8eOP8OGHbzJnzid2hyVElslKmXA7qampfP311/Tp4yAs7FlULx874oAvvgD/4vDBUTgSa0sYxtOB7cAc1IGAtkADG+M5hCqbcQ4oh9pvlsfGeFygAfn/hNw67Npuz4qu+Ne5c5AjBwQGNmH27A78+ec+Pv74Y3Lnzm13aEL8h2z0F9lWSso8nnuuM88+m5RWv8gesQnQ/HHYvBeVxFjcVN0UV4EFQDVU30on9v26rr+lmQLUBppi34lRV5wAx1V45H54tSmsmAndukFgoN2Bebcvv4RRo4LYvHkifn73EBcXR0hIiN1hCfEPuX0psq19+4rx88+BLFmSz9Y4cgZC64JQNyeEBdkainGCgS6ohAxgHvAz9hSdvf6WZk1gA/A5KmmMtyEeVxSF1HLw42YoNwx694Yx4+0OSgwaBFu2xOPndy9OZ3Xq1CnP888PszssIVwiK2XC7Z09e5b8+XWgCxcvrrW1ZpTTCccvQ/exsG4+9m2WN8M6YDEqCaoKtATsuvtzHrVqtg11crMRqhCtv03xuOIwBJSCAfWg1iVIiVOJmo8nrKpmUwkJ8MEHUKNGUbp0eYvExPtZvXo9zZo1k1IawjZy+1J4hEOHdlGvXk3efTeRfv3sjWXESBj5FugDwWlnv0mjJaD6Vq5B7T3rDNSwMZ5TqNOie4EQ1C3N2rh9116fSZD3KmzZDIVCVYFiyQHs98MPeXjssQusWjWfBg1a2x2O8FKSlAmPEB8fzzPPDOH55zXKlh1tayy6Drt2QWJuGDgLVs8DyuLeKzkZcQm1UtUIyAtcBgKxr2zFUVQ/zyOo4rONUMmiu5bR0IE4CMwFD1WEZS/Bqy/Bww/bHZh3i4+HadOge/ccaFpfPv88jFOnnLz99ts4HJ6y5C3cnSRlwgP9yPvvP8lddyVTubK9kezeDRUrQo4OcKWevbGY5lfUXrMmqGTIjuRTBw6gbrGeQO2Jq5v2sL51quuuADOh/j3w/uNQLSecPAlSUst+AwfC0aP5mT59HNCe/fsPUKZMGbm1KUwlSZnwOGfPnqVSpfL07x/CO+8ctzscli2DiCrwcRR8Ng1SdLJ/m6brHUatVB1DJUD1UMmQHYcedNSK2UpgH2q1rCaqrEc2qIAQvhVOTIa9+6BsGbujEU6n2vcXE1OG8PBjDB06hJEjP7A7LOHBJCkTHunMmTPkzBlKQMC37NjxEr6+iZR3g36KTVpC1BZIGIhnlM64RkfdRlyBSoaaog4D2OkMsArYioqvEurWZhE7g0pHHLAPSreAp+tAzALIn0udGhT2SUiAiROhQYNgKlXqwd69rXjttam8/fbblCkj2bMwjiRlwuO1bFmfw4c3s3dvIr42bwI/dw727YMzueHZuXBoBVABt9+cniGnUDXEQlCb8HehkiG7KpdcQh1O2IBqfl4CdSCgIu677yyN41coXhB+GgfNSsKWLVCtmhwMsNv06fDYYz5s3vw8hQs/zZ49V/Dx8aFs2bJ2hyayOUnKhMeLjo7m8OEDNGiwEl1/jRMnkggPtzsqmDoDunaGPL3hgqf+sL0GVU8sBSgPRAJlsKdUSAIqMYsCYlC3V6ujErT8NsTjqhTAF8r4wMHX4dU34M1X7Q5KXLu1CT706FGAefMuc+rUKfz8QnA6nfhIvRORCZKUZcUzz8DmzeaNLwwXHX2EffsPU6sW5LC5kLcOXLwIYWFwOg4OnYDkRMCTymiAqsh/Oe2RijqpWdDWiFSCdhnVuQBUvbNQ1J44d12F0oGroAVC/jDIqcGFU1CurHQKsFtiIsRdhTy5fYH8bNoUS2hoblk58zQ1asBnn5k6xZ2SMk+6oSIEAHnyFKJYeCohIX7AUa7GJxMUaM/tIA3InUv9d5FQuOyAmCRw+MFVe1p6msMHVaoiJyoJuvZ7raMKwYZg/aGAwLSHE7WP6zKqv6YG5Eh7uFsJEw0IUb9tZ+LgzFXQLsGxK1DIAc54SEmBvHndN6/0VAEB6gEp6ESTMycEBScCOrqej527oilapAi57KxuLbI9WSkTHu3q1bOULVuGtm0T+Okn98iCrlyB4BD4ZSP06whJDVCnBz3RSWACKlHLBdRC9drMZUMsqahTmxtQe+CcqNOalYEqqJU9d8x0dP6JK8cU8D0NKzdDpQJw+DAUKyZdA+x24AC0auXgk0/acO+9wzl1KoJp02bQrVs38uTJY3d4ws3I7UvhtXRd5++//yY8PIx69RZx6dLnrF0bT+vW9m+kPn0ahgyBGnfBKh+YuQlSjwIReNapzRRgNyoZOpT22tOoQwGp2LP3LD4tpu3AQVTikw+VnFXBvgML6XECF4G8UK0AHH8b6teBmZPV20lJ4O9uq39eQtchNVUlyOPHh/Hww5fYtu1TqlTpx/79pzlz5gz16tWTfWhCkjIhrvnoo9d5/vkR7NjhS6VKKXaHc4M334M3XoTCwyE6EPsSFjPFoMpp1EGt/kxF3d6sgipnkcOGmOKAnagE7Ujaa4VQJzcj0v7bHVfQUlErfoFQsR60Lw5jesF778LAAXYH5910XRWVrlABNM2fYcMK8cUXJzh/fh2hobU4cOAAuXPnllU0LyVJmRBpEhMTmT9/PnfdVRkYyeef/0zhwk66d7c7MrVXaPVqaNQY5u6HgU/D4SOgP4h7JgVGWA1sBM6ifo0lUSclq9gUzyX+TdCu1SQOBcqlPUqjDgy4ozhgBeSsBXe3hlq+8MPL8P0YqFvX7uC8W0yMOi/WogVAWbp0cbJnTxJ79uwGcrB9+3aKFy9OWFiYvYEKS0hSJsQt6LpOnTo1qVzZl3HjYoH9REdD4cJ2R6Z8+insPwpF74VftsGOv4BwVGLgac6gEqEdqB6iHVArQYtQv97iWH8s6QqwH1WH7QCq/pkPKnG8lqTlwX0T5hOgzYNGQ6BzXfDbD3N+hZ/Gus/fcW+1bh2cOQN33eUL1KF8+T1EREQwffoiIIioqCgqVapEcLA79w8TmSVJmRC3kZqaypUrVwgLC+XIkV8oU6Y3332n0a9fqt2h3SA+HsKLQ/VOcKkhbDiBunVVBnXC0FPoqH1TvqhE7bu05/5AKVQiVBF1mtNKKagWU3vTHufTXg9DJWnXHrlx3yRtB/ishk7vQeuycGox7F4Pv/+O7QWXvZmuw+LF6mRno0YBJCTUJVeuVQwceB8ff/wTuh7IX3/9RZMmTShY0O46M8IIkpQJ4YLz588zatQoevduS/Hic4iK+po//zzNCy+AO2z9SE5WrWBCQ2HqQujaGsr1g/3FQE9E1eXytPpniai+m/vSHrFAb1Qyehp1urME1idD51GrZ4fTHtdqoWWnJG01BByF+96CpsVh8dfgmwITxqu3k5PBz827IXii5GRYuFCdqq1c2Y99+8oTEbGdMWP689hjL3HhQhg//vgj3bt3p3jx4naHKzJBkjIhMuHTTz/hzTdf4+TJ3gQHz2Lp0qPEx0O7dvaf3HQ61f6zqlXhqgNe/gzGvga5noWLOVH7izRUoVRPoaP2nuVBraQtSXuA2vdVIu1RA2tbK12L6zD/TdJyAEWvexTBniburlgCfk5o+hg0CIc/h0L1yvD7RPX23r1QosS1Wl3CKikpqvVW8eKQPz/Mn5+btm1jWLLkCZo1e4CNGx18+un3jBw5kpIlS6LrOprd/0CJO5KkTIhMio+PJyhIfRft3Lkpe/fuYPfuImjadjZvhtKlVbV+ux0+DNOmwcCBsP0sDHsRFo6H4FdVfkYsaoO6J93qTEUVgz1y3SMReAG192s1qmBskbSHVatW1ydpx4ET/Hu7EyAvKkELT4urAO5XxBZU4/kcUL4V1CsKfz4Mbe+Gn3+AsAD47jto0gQqVbI7UO9z+jTkzq3Kn0ybBk895cP69W0pUqQhY8eeY8SIv1m9eh2FChXi2LFjXL16lXLlyuFweNpx7uxJkjIhDJCQkMDhw4epUKECur6fUqXqEhnpw6RJl4FEVqyAKlXAHQp6b92qNhP36QvrTsATfeDgDggaBhfiUfXCAgFP2vCtozLQa2U1JqMODzjTngeienPek/Y8Nu1aK8pGxaNutZ5IexxPi/WaPKjitQXSvhZEJZHu8j3UiarrlhO0cCjhB4dfhg4D4OlnICIHDHtKdaVr3lyt5MbHQw47Spx4ufnzYdw4+Pnn/DgcNXnxxQt89NFGrl5djJ9fNf76ayG7du3i5ZdfRtM0WVmzgSRlQhhM13VWrVoFQKNGkcTFLSVXrg4MG1aCd989R2rqZf7+W32DyucGhUiXL1enve67D/ZfgNaNQA+C8KdhUzTEL0c17Pa0Nn4pqAMD0aikKAhonfbex6jEKC8qGcrPv3vBzKajym+cRO2NO5P29foVNb+0uPKlxXjtkQf3WFm7jEpog4Hz4PM7RDwIjVtDWAx83Ae+mwD9H4JTp9SKTteuIHvVrbV3L+zYAfek/TDy9NPBzJuns3dvf6Ayffr8xcGDMSxfvhbQmDNnDj4+PrRp08bOsD2aJGVCmCw5OZmVK1dSpEgRIiLKsH37X1St+gBjx9bmkUfOEB19jG++gUcegTJl7I4Wjh6Fy5ehcmVIdkKRolC1GdR4FLafgSUvQ3JloH7aB06hEgJP2fitA1tRtxnPpH2NAeoCHVHJ3HeodlB50h65USuLoSbGlZQWy82J2qWbrsvJjUlarusegdh/uCAW2AJUhZACUPAoHPwOHvsWmjSAM5vh25EwayaULw8HD8K2bdCmDUgVCPNdf4jj++/VD2wvvRQIlKZx41P4+gazZMnLQFm6d/+UAgVK8NVX3wAwZcoUChUqRP369W87vrgzScqEsFhycjKbNm2idOnS5MuXj4ULJ9OuXTeWLx9MgwaXWLp0Bf/73z7GjdOpUAHOn4cLF9QeNTu6sDidEBen9sclJkKfPtCyHVRuBVGH4JnmUKYbpDSCY+cgdSaqj2Vx1K2tBNSKid3JQFYkoZKxYNRG/enAhbTHtbaprYHGqKRjCuq0Zc7rvhbGnK4ESajk7ObHOdQ+uuv5c2OSlistvtDrHlYn16moxDIH6pDGUWANBN8DpYuCtha2jYORM6FyKdg4E34bDevWqr1T69appO3hh1UykZCg9lPJFinjxcRAbCyULKmeDxsGefNqvPhiEaAEJUtuoXHjEkyY8DRQnFq1nqd9+/a8884nAHzwwQfUrl2bVq1aARAdHU2+fPnwk6O8/5CkTAg3kJCQgK+vL76+vixevJiRI0fw++9vkT//CUaNmshTT03j2LH6hIcfZ/bs40ydCh99pPblnDunfrotVMj6k5/x8TBzptovV6EC7NkHTZtAvxehWENYtxl+ehxK9IPYcnAxGliGWmUriErYYlC34bLjv8vX9qpdQCU3uVCrWNNQydmVtGsA7gOqopKO6agEKEfaIwTV/DwXKslKRt1OzUpioaP2q11Me8Re99/XHjcnbaBW065P0q7Fd/MjGGuK9iagkszCqN+PPcBW8O8OxXJB8lw4Pg9emAWFw2DBdzD/d9hyFIqEwu8TYONG+OorNdzGjXDpkto+AOqHDmk5aYwTJ1SPz2LF1PMhQ6B2bXj44SB0vRAhIUcYMiSCd99tT2pqIQICXuaFFx7krbeG43TmpVWrHjz55JM8+OCDpKSkMHr0aJo1a0blypVxOp1ER0eTP39+Ajz4mK8kZUK4uYMHD7J8+XJ69+6Nw+Hgq68+5Z133uXYsdH4+Bzm5Zd/5f3315OY2BAfn2hGjTrK4sUp/P67+vz69WqlrV079dyKlQRdVwni6dPw22/QubNa6Vu8HB68H179CvKUh0Xz4YfnoPkISCwCRzbDqb8g9V7UPq4zqE3ktVGJQBwq0cmL9VX8M8OJ2l91CRVzCGoz/3LUr+PaIwV4FLW6uAV1EEFDJWZBqASoa9oYx1GHMYJRCdS1R0HU74mO66uS8WnxXf+4dNPzK6jVrFsJSPs1Bd30CLzF8wD+PeXrj3GHKFJQfy+u1eE7gNon2Fg99VsCjkNQ93XIHwJbv4Fze+H1vyF/MIx5CaIPwfxVkDcIPn5PJW0ffaQ+P2mS+tqtm/q6e7cq/VGqlHp+7e+6SJ/TCUlJEBSkfpAcOxaqV4d69dTveefO8MQToTz0UHGio3NSpMgqvvmmOU891YwTJ3wID3+N774bxOOP9+L48UQ6dHiC99//kI4dO3Hq1Cnee+89Hn30UapVq0ZMTAwLFy6kcePGFCpUiISEBM6fP0/+/Pnx93eHjZe35jZJmaZp7YHPUf+rfq/r+ns3vR8A/Iz65/k88ICu64fvNKYkZcIbbNy4kS1bttC3b18APvzwQ+bOncGCBV8BJ+nTZwRLl+7i8OHHgXP06jWXdevOsXdvSeAcI0ac58wZ/Z+VhIkT1T+cacMRFaWquteooZ6fO6eSOiPKfZw5AytWQMuW6mTqypXw3nvw2ZcQkAd+HAuvPwOfzgJyw/w/YdYn0HU0JIXA7jlwZDLkfR5iHZC4HdiGSmD8UaUwjgINUf+ynEOtEJVBJS5xqFWpXFn/tWSajlqx8k17nEMlFldQSdPVtK9dUYnHSmD+LcYZilqtW4IqWXEtEfJPe/RM+7oT9Xvif93DD1XDTUubPz7tNb+0mJyoxCwuLZ64mx4JaZ+59khw4dftx42Jmt9NMfnf9Nr18dzuv689fLj9KuPFtBivnS7elvZrSNsG5ZgNfslQoi/kDoS9H4GfL9z3vno+YQgE54CXR6vyH//rrvZdfjdePR/cHyLKwsiRaryXX4ayZf/9/2nsWFXTrWVL9XzBAihaFCpWVM/371cFqa8VpfbW27Gpqerfh6AgyJlTJW2//65KrVSooPa+DhkCQ4dqNGmSiy1bgmja9BS//FKJTp1Ks3p1Mg0bzmb27B60b1+LZctO0azZxyxYMJJWrRqzdOkeHnzwVaZOHUvduk1Zu3Y7r7/+Jl988QURERFs2bKFcePG8fzzz1PYwt5jbpGUaZrmg2pQ0gb1c+B6oIeu6zuvu2YAUE3X9Sc1TXsQuEfX9QfuNK4kZULAyZMnuXDhAlWqqE7ekydPJjo6mgEDBgAwbNhznD59nPHjPwJiad36ERIS4lmx4lUgliZN3sPXV2fx4nuBOBo0+JuwMAdz59YA4mjefDPFizv4+eecQALdu1+kXDknb7+t/v0YOFB9w3n6aRXP66+r5w8+qJ5/9pl6fm0lb+xY9Y9ugwbq+S+/qHpXNWqommtjxsBDD6mDCIsWwZdfqlWNMmXgm1Hw/vvw6wwIyQufvgM/fwXfLweCYeInsORXeHYyOP1g4SjYORs6/wCJGmz9Cc6sg/JvQZIG0VMgfh8EP67uKiYsRh1suA/1TX8N6rZgm7TnG1Hf8BuiEpxtqJWcGmnP96CSsAppfzgH0l6/1rP0CCqhCE97foJ/T1qC2tzvh7qlGI9aEbrWeqoMapUrOu1zCfy7CpcC9En7zHJgE+oXdI0P8GraZ2amxX29IOB/aeNMRnVQuD4JygU8nPZrm5f2e3T96lEw6vZsImrl89rqW2pa/NdWBZNQq3PJaa87yTxH2sOHf5M7n7Q5fa5779r7AWlxOK97z5H22rWTpBrqVrWDf8uSnEj7bMm09zepa/0rQ6A/XJ0PgQWhaEv1fOcoyBsBVbtAkB/MeRdK1YIWD0GgH3wzCGq1gC6PqOf/6w6t74GHnwR/X3isK3TvDb36go8OfXtA7z7wUA9wJkP//tC7N3TqpPaC/u9/8MADKgm8fBnefVedtK5XT+0PGzVKrVBVrar2jP36q/p/sXx59XzaNPXZEiXU84ULoXFjlUheuKAKVdevr4rXXmuuXquW2u938SLs2qW2N4SGqueHDqmxg4PV/CdOqP93AwJUfKdPq7n8/ODKFTVHkSLqB8O4OPWZggXVLeerV9VrefOqpDUhQW2pyJVLfd25U/VyLVIEoqNh8mRo317Nt307fP45PPec+vdm6VJ4/nmNsWNDqVw5J3//nUqfPtGsXVuDSpXyp/0FfygLfyFdc6ekDF3XLXkADYC51z1/EXjxpmvmAg3S/vvaz5PancatXbu2LoTImNTUVD0hIeGf51u3btU3bdr0z/NJkybp06ZN++f522+/rY8aNeqf5/3799ffeecdXddTdF2P09u0aa6/8MJAXdcP6v9v7+5j7KjqMI5/n7svfdm21NICpaUvhi0oSNdCgaSISkRAGqkRAhFQgoKSglATY9VgA8FEJRiQEBJiq9EFKpVaqhgDoQU0okBreWlBKJEoxKSpkVAUqMDjHzMty+12l0rdGfY+n2Syd+acOfO7e/bld2fOzLGf9KGHzvQll5xp+/e27/OECWO9cOF826tsr/DYsaO8aNFJtpfZvsmdne1evPhjtq+x/T0Dvvzyj9i+wq+++k0Dvuqq421/1du2XWzAV189z/ZCb916ngFfd92xts/3Qw/NN+Abbzza9tlevfp4A162bK7tM3zDDT0G3Nt7lO0FXrLkEANeseJI26f6ooumGfDyn8/21n+d4PkLJhX7r3i/N2451nPnjTXgH9za7Xv/0uNZHxhlwN+5ZaZXbjrM02eNMOAlvdO8dH23D5jRYcCLe6f42gdmet8p7QZ8We9kX3nvQd5nvzYjvLB3fy+++0B3TWhYwhf0TvKXf32AR46TJXzezRP9xV/u584uWQ18zi0TfN6qiW4fidXAZ936Hn/m9gludBTrZ9w23p9ePt4I04ZPun6cFywfb8DqwMd8bbTnXFrE3mjHh18w0icuLd5boxNPP6XTB364Y2f9qSd2+EPXdxXlI/A+s9rcNa1RlLfhCT1tnvvd0UV5B24fXbRDkcZ5XHfDs5eMLOq3v7l9Z/n7Gu6+qHNne83lXTMa3v+jxfcO7Vo+Yj95zMGNXbbvWNrH4I59+i/b2WY/7WbZO0ujMUh5P33+lv7r52em79LRMXB5Z+fg5W1t8po1a/b639tmwMPeTU4zlGfKTgdOtv2Fcv1c4BjbF/ep83hZ57ly/Zmyztamti4ELixXD6H4bPr/NJEiQYx6Sb/UT/qkntIv9ZM+qaeh6Jfptif1V/BuGEa7C9s3ATcN1fEkPezdnWqMyqRf6id9Uk/pl/pJn9RT1f0ylMMKnwcO6rM+tdzWbx1J7RRDXv9BRERExDA3lEnZQ0C3pJmSOoGzKJ7009dqiqGqAKcDazxU11cjIiIiKjRkly9tvybpYorB/G3AMtsbJV1JMehtNbAU+KmkzRT3v5w1VPENYsgulcYeSb/UT/qkntIv9ZM+qadK++Vd//DYiIiIiOGgxR5VFxEREVFPScoiIiIiaiBJ2SAknSzpz5I2S1pcdTwBkpZJ2lI+1y5qQNJBktZK2iRpo6RLq46p1UkaKelBSY+UfXJF1TFFQVKbpD9J+lXVsURB0rOSHpO0QVJl0wRlTNkA3s7UUDH0JB1PMYnLT2wfXnU8AZImA5Ntr5c0FlgHLMjvSnUkCeiy/ZKkDorZMi+1/YeKQ2t5kr4CHAWMsz2/6niiSMqAo5ofVj/UcqZsYEcDm13MHbMdWA6cVnFMLc/2/RR350ZN2P677fXl623AE8CUaqNqbeWMLi+Vqzum9M6n8IpJmgqcCvyw6liifpKUDWwK8Lc+68+RfzQRA5I0A/gg8MeKQ2l55WWyDcAW4G7b6ZPqXUsx9fsbFccRb2XgLknryqkcK5GkLCL2GkljgNuBy2y/WHU8rc7267Z7KGZQOVpSLvdXSNJ8YIvtdVXHErs4zvYc4BRgYTlMZsglKRvY25kaKiKActzS7cDNtldWHU+8yfYLwFrg5IpDaXXzgE+W45eWAydI6q02pACw/Xz5dQvwC4rhS0MuSdnA3s7UUBEtrxxUvhR4wvb3q44nQNIkSePL16Moblh6stKgWpztr9ueansGxf+TNbbPqTisliepq7xBCUldwMeBSu7uT1I2ANuvATumhnoCuM32xmqjCkm3Ag8Ah0h6TtLnq44pmAecS/HJf0O5fKLqoFrcZGCtpEcpPmDebTuPYIjY1f7A7yQ9AjwI3Gn7N1UEkkdiRERERNRAzpRFRERE1ECSsoiIiIgaSFIWERERUQNJyiIiIiJqIElZRERERA0kKYuIiIiogSRlERERETWQpCwiWpakUZLuk9RWrr9ePvj2cUkrJI0ut0+VdIekpyU9I+m6cpaPHe28VH7tlHS/pPZq3lFEvJslKYuIVnY+sNL26+X6y7Z7bB8ObAe+VE4htRJYZbsbmAWMAb7d3Jjt7cA9wJlDEn1EDCtJyiJiWJJ0r6RDy9f7SupvLruzgTt208RvgYOBE4BXbP8IoEzgFgHn7ziT1mRV2W5ExB5JUhYRw9XBwFPl6yOAx/oWlpcf32v72eYdy8uPp5T7HAas61tu+0Xgr+Uxmj0OzH2HsUdEC0pSFhHDjqTpwPO23yg3HQE82lRtIvBC07ZRkjYAD1MkXUv39NjlmbTtksbu6b4R0doyGDUihqPZvDUJOxL4WVOdl4GRzdts9/TdIGkTcHrTtnHANGDzbo4/Anhlz0KOiFaXM2URMRz1UCZckrqB02i6fGn7n0CbpObErNk9wGhJny3bawOuAX5s+9/NlSXtC2y1/Z93+iYiorUkKYuI4Wg20JD0CPAtYBPwuX7q3QUcN1BDtg18CjhD0tMU49ReAb6xm10+Ctz5P8YdES1Mxd+biIjho0ye5tjeNki9OcAi2+fuxWOvBBbbfmrQyhERfeRMWUQMK+UAew+WkFFUWg+s3fHw2L1w7E6K55klIYuIPZYzZRERERE1kDNlERERETWQpCwiIiKiBpKURURERNRAkrKIiIiIGkhSFhEREVEDScoiIiIiaiBJWUREREQN/Bckuvt9+h17HAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(10, 7))\n", "artists = brazil.plot_results(poi_tests, tests, test_size=0.05, ax=ax)\n", "print(f'expected upper limits: {exp_limits}')\n", "print(f'observed upper limit : {obs_limit}')" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[1.0, 1.0, 1.0]\n" ] } ], "source": [ "source = {\n", " \"binning\": [2, -0.5, 1.5],\n", " \"bindata\": {\n", " \"data\": [120.0, 145.0],\n", " \"bkg\": [100.0, 150.0],\n", " \"bkgerr\": [15.0, 20.0],\n", " \"sig\": [30.0, 45.0],\n", " },\n", "}\n", "\n", "\n", "my_observed_counts = source['bindata']['data']\n", "\n", "model = uncorrelated_background(\n", " source['bindata']['sig'], source['bindata']['bkg'], source['bindata']['bkgerr']\n", ")\n", "data = my_observed_counts + model.config.auxdata\n", "\n", "\n", "binning = source['binning']\n", "\n", "nompars = model.config.suggested_init()\n", "\n", "\n", "bonly_pars = [x for x in nompars]\n", "bonly_pars[model.config.poi_index] = 0.0\n", "nom_bonly = model.expected_data(bonly_pars, include_auxdata=False)\n", "\n", "nom_sb = model.expected_data(nompars, include_auxdata=False)\n", "\n", "init_pars = model.config.suggested_init()\n", "par_bounds = model.config.suggested_bounds()\n", "\n", "print(init_pars)\n", "\n", "bestfit_pars = pyhf.infer.mle.fit(data, model, init_pars, par_bounds)\n", "bestfit_cts = model.expected_data(bestfit_pars, include_auxdata=False)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs4AAAD8CAYAAABnwSKjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAU/ElEQVR4nO3dbbBdZ3ke4PsJFlA7jGyCcYw/EDRuGVNNDFUdEjodKrXEuJ3YmVJq6oCTOiPS4IyZ8MeJZorcVNOk08SYhrijxExMogIuCcUCksZIdJj8MMSmwPFHAUEssMbY5sMCxhMqu29/nCV8EJLOez72WXsfXdfMnr3Wu9ba+3llP7NvLa29V7XWAgAAnNwPjV0AAADMAsEZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgw6LBuaqeXVWfrKrPVNV9VXXjMP6iqvpEVR2oqvdV1TOH8WcN6weG7ZsmPAcAAJi4njPO302ytbX240kuSXJZVb0iyW8luam19mNJvpnk2mH/a5N8cxi/adgPAABm2qLBuc37zrC6YXi0JFuTvH8Yvy3JlcPyFcN6hu3bqqpWq2AAABjDaT07VdUzktyT5MeSvDPJF5M83lp7ctjloSTnDcvnJflKkrTWnqyqw0l+JMnXjnnN7Um2J8kZZ5zx91/ykpesbCawjtxzzz1fa62dPXYdC+lZOD79CrNlJT3bFZxba08luaSqzkzygSQr7sDW2u4ku5Nky5Yt7e67717pS8K6UVUHx67hWHoWjk+/wmxZSc8u6Vc1WmuPJ/lYkp9McmZVHQ3e5yc5NCwfSnLBUNhpSTYm+fpyCwQAgGnQ86saZw9nmlNVfyvJP03yQOYD9GuH3a5J8sFh+Y5hPcP2/a21too1AwDAmuu5VOPcJLcN1zn/UJLbW2sfqqr7k7y3qv5Dkv+d5NZh/1uT/FFVHUjyjSRXTaBuAABYU4sG59baZ5O87DjjX0py6XHG/ybJv1yV6gAAYEq4cyAAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHRYNzlV1QVV9rKrur6r7qur6YXxnVR2qqk8Pj8sXHPNrVXWgqj5XVT89yQkAAMBaOK1jnyeTvLW19qmqek6Se6rqzmHbTa21/7xw56q6OMlVSV6a5AVJPlpVf6e19tRqFg4AAGtp0TPOrbWHW2ufGpa/neSBJOed5JArkry3tfbd1tpfJzmQ5NLVKBYAAMaypGucq2pTkpcl+cQwdF1Vfbaq3lVVZw1j5yX5yoLDHsrJgzYAAEy97uBcVT+c5E+SvKW19q0ktyT520kuSfJwkt9eyhtX1faquruq7n7ssceWcigwAj0Ls0O/wmR0Beeq2pD50LyntfanSdJae6S19lRr7f8l+f08fTnGoSQXLDj8/GHs+7TWdrfWtrTWtpx99tkrmQOwBvQszA79CpPR86saleTWJA+01n5nwfi5C3b72ST3Dst3JLmqqp5VVS9KclGST65eyQAAsPZ6flXjlUnekGSuqj49jP16ktdX1SVJWpIHk7wpSVpr91XV7Unuz/wvcrzZL2oAADDrFg3OrbW/TFLH2fSRkxyzK8muFdQFAABTxZ0DAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQIdFg3NVXVBVH6uq+6vqvqq6fhh/blXdWVVfGJ7PGsarqt5RVQeq6rNV9fJJTwIAACat54zzk0ne2lq7OMkrkry5qi5OckOSfa21i5LsG9aT5DVJLhoe25PcsupVAwDAGls0OLfWHm6tfWpY/naSB5Kcl+SKJLcNu92W5Mph+Yok727z7kpyZlWdu9qFAwDAWlrSNc5VtSnJy5J8Isk5rbWHh01fTXLOsHxekq8sOOyhYezY19peVXdX1d2PPfbYUusG1piehdmhX2EyuoNzVf1wkj9J8pbW2rcWbmuttSRtKW/cWtvdWtvSWtty9tlnL+VQYAR6FmaHfoXJ6ArOVbUh86F5T2vtT4fhR45egjE8PzqMH0pywYLDzx/GAABgZp222A5VVUluTfJAa+13Fmy6I8k1SX5zeP7ggvHrquq9SX4iyeEFl3QAsEpuuvPzY5cAcEpZNDgneWWSNySZq6pPD2O/nvnAfHtVXZvkYJLXDds+kuTyJAeSPJHkF1azYAAAGMOiwbm19pdJ6gSbtx1n/5bkzSusCwAApoo7BwIAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBmXVn586dY5cAdPrzd/+XsUsAluBU/4wVnFl3brzxxrFLADr9xR//7tglAEtwqn/GLhqcq+pdVfVoVd27YGxnVR2qqk8Pj8sXbPu1qjpQVZ+rqp+eVOEAALCWes44/2GSy44zflNr7ZLh8ZEkqaqLk1yV5KXDMb9XVc9YrWIBAGAsiwbn1trHk3yj8/WuSPLe1tp3W2t/neRAkktXUB8AAEyF01Zw7HVV9cYkdyd5a2vtm0nOS3LXgn0eGsZ+QFVtT7I9SS688MIVlMFquunOz49dwqpYL/OYJnoWZod+hclYbnC+JclvJGnD828n+TdLeYHW2u4ku5Nky5YtbZl1wPfMzc1l//59SZKbb357tm7dls2bN49c1fqhZ1lN+nWy9Ov0mfUTOgt7dtM//+VTtmeXFZxba48cXa6q30/yoWH1UJILFux6/jAGEzU3N5e9e/fmyJEjSZLHHz+cvXv3Jskp2dgwzfQrzBY9+7Rl/RxdVZ27YPVnkxz9xY07klxVVc+qqhcluSjJJ1dWIixu//5932voo44cOfK9vx0D00O/wmzRs09b9IxzVb0nyauSPK+qHkrytiSvqqpLMn+pxoNJ3pQkrbX7qur2JPcneTLJm1trT02kcljg8OHDSxoHxqNfYbbo2actGpxba68/zvCtJ9l/V5JdKykKlmrjxo15/PEfbOCNGzeOUA1wMvoVZouefZo7B7IubN26LRs2bPi+sQ0bNmTr1m0jVQSciH6F2aJnn7aSn6ODqXH0ywn79+/L4SRnnrnxlP3GL0w7/QqzRc8+TXBm3di8eXM2b96cX/3wLbn++reMXQ5wEvoVZouenedSDQAA6CA4AwBAB8EZAAA6CM6sO6/+uevGLgHopF9htpzqPSs4s+5c9sZfGbsEoJN+hdlyqves4LxKdu7cOXYJwBLoWZgd+pVp4efoVsmNN96osTllPPKtv8lNd35+7DJW5MYbb8zGV/7rscsAOviMZVo44wwAAB0EZwAA6CA4AwBAB8F5hfbs2ZNNmzYlSTZt2pQ9e/aMWxBwUnNzc7n55rcnSW6++e2Zm5sbtyDghHzGMm18OXAF9uzZk+3bt+eJJ55Ikhw8eDDbt29Pklx99dVjlgYcx9zcXPbu3ZsjR44kSR5//HD27t2bJNm8efOYpQHH8BnLNHLGeQV27NjxvYY+6oknnsiOHTtGqgg4mf37930vNB915MiR7N+/b6SKgBPxGcs0EpxX4Mtf/vKSxoFxHT58eEnjwHh8xjKNBOcVuPDCC5c0Doxr48aNSxoHxuMzlmkkOK/Arl27cvrpp3/f2Omnn55du3aNVBFwMlu3bsuGDRu+b2zDhg3ZunXbSBUBJ+Izlmnky4ErcPTLCTt27MjBgwfzwhe+MLt27fKlBZhSR78AuH//vhxOcuaZG7N16zZfDGRdm9k7fT7/H+Rf/ca7s3//vhz88C154T/7t9m6dVseff7m2ZwP64LgvEJXX311rr766lRVHnzwwbHLARaxefPmbN68Ob/64Vty/fVvGbsc4CT0K9PGpRoAANBBcAYAgA6CMwAAdBCcV8nb3va2sUsAluDVP3fd2CUAnfQr00JwXiU7d+4cuwRgCS5746+MXQLQSb8yLQRnAADosGhwrqp3VdWjVXXvgrHnVtWdVfWF4fmsYbyq6h1VdaCqPltVL59k8QAAsFZ6zjj/YZLLjhm7Icm+1tpFSfYN60nymiQXDY/tSW5ZnTIBAGBciwbn1trHk3zjmOErktw2LN+W5MoF4+9u8+5KcmZVnbtKtQIAwGiWe+fAc1prDw/LX01yzrB8XpKvLNjvoWHs4RyjqrZn/qx0znr+C9w+E6bcsT0LTC/9CpOx4i8HttZakraM43a31ra01racsfGslZYBTJiehdmhX2EylhucHzl6Ccbw/OgwfijJBQv2O38YAwCAmbbc4HxHkmuG5WuSfHDB+BuHX9d4RZLDCy7pAACAmbXoNc5V9Z4kr0ryvKp6KMnbkvxmktur6tokB5O8btj9I0kuT3IgyRNJfmECNQMAwJpbNDi31l5/gk3bjrNvS/LmlRYFAADTxp0DAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOhw2koOrqoHk3w7yVNJnmytbamq5yZ5X5JNSR5M8rrW2jdXViYAAIxrNc44/+PW2iWttS3D+g1J9rXWLkqyb1gHAICZNolLNa5IctuwfFuSKyfwHgAAsKZWGpxbkr+oqnuqavswdk5r7eFh+atJzlnhewAAwOhWdI1zkn/YWjtUVc9PcmdV/Z+FG1trrara8Q4cgvb2JDnr+S9YYRnApOlZmB36FSZjRWecW2uHhudHk3wgyaVJHqmqc5NkeH70BMfubq1taa1tOWPjWSspA1gDehZmh36FyVh2cK6qM6rqOUeXk7w6yb1J7khyzbDbNUk+uNIiAQBgbCu5VOOcJB+oqqOv899aa39eVX+V5PaqujbJwSSvW3mZAAAwrmUH59bal5L8+HHGv55k20qKAgCAaePOgQAA0EFwBgCADoIzAAB0EJwBAKCD4AwAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdBCcAQCgg+AMAAAdBGcAAOggOAMAQAfBGQAAOgjOAADQQXAGAIAOgjMAAHQQnAEAoIPgDAAAHQRnAADoIDgDAEAHwRkAADoIzgAA0EFwBgCADoIzAAB0EJwBAKDDxIJzVV1WVZ+rqgNVdcOk3gcAANbCRIJzVT0jyTuTvCbJxUleX1UXT+K9AABgLUzqjPOlSQ601r7UWvu/Sd6b5IoJvRcAAExctdZW/0WrXpvkstbaLw7rb0jyE6216xbssz3J9mH17yW5d9ULWVvPS/K1sYtYBethHuthDn+3tfacsYtYSM9OJXOYDvp18tbD/yfrYQ7J+pjHsnv2tNWupFdrbXeS3UlSVXe31raMVctqWA9zSNbHPNbLHMau4Vh6dvqYw3TQr5NnDtNjPcxjJT07qUs1DiW5YMH6+cMYAADMpEkF579KclFVvaiqnpnkqiR3TOi9AABg4iZyqUZr7cmqui7J/0zyjCTvaq3dd5JDdk+ijjW2HuaQrI95mMPkTXt9PcxhOpjD5E17fT3MYXqsh3ksew4T+XIgAACsN+4cCAAAHQRnAADoMEpwrqrnVtWdVfWF4fmsE+z3VFV9enhMxZcLF7uVeFU9q6reN2z/RFVtGqHMk+qYw89X1WML/ux/cYw6T6aq3lVVj1bVcX+btOa9Y5jjZ6vq5Wtd42I65vCqqjq84L/Dv1vrGhfUomdHNOs9ux76NZmdntWv45r1fk3WR89OrF9ba2v+SPKfktwwLN+Q5LdOsN93xqjvJHU/I8kXk7w4yTOTfCbJxcfs88tJ/uuwfFWS941d9zLm8PNJfnfsWheZxz9K8vIk955g++VJ/ixJJXlFkk+MXfMy5vCqJB8au86hFj073XOY6p5dD/3aOY+p6Fn9OvVzmOp+HWqc+Z6dVL+OdanGFUluG5ZvS3LlSHUsVc+txBfO7f1JtlVVrWGNi1kXt0NvrX08yTdOsssVSd7d5t2V5MyqOndtquvTMYdpomfHM/M9ux76NZmpntWv45n5fk3WR89Oql/HCs7ntNYeHpa/muScE+z37Kq6u6ruqqor16a0kzovyVcWrD80jB13n9bak0kOJ/mRNamuT88ckuRfDP/88v6quuA426dd7zyn3U9W1Weq6s+q6qUj1qFnx3Mq9Ox66ddkOnpWv47nVOjXZP307JL7dWK33K6qjyb50eNs2rFwpbXWqupEv4n3wtbaoap6cZL9VTXXWvviatfKD9ib5D2tte9W1Zsy/7f7rSPXdCr6VOZ74DtVdXmS/5Hkokm9mZ6daXp2OqxZz+rXmaZfp8Oy+nViwbm19k9OtK2qHqmqc1trDw+n9h89wWscGp6/VFX/K8nLMn/t0Fh6biV+dJ+Hquq0JBuTfH1tyuuy6Bxaawvr/YPMXy83a2b+tu+ttW8tWP5IVf1eVT2vtfa1Cb2fntWzY5n5fk3Wtmf1q34d2cz37HL7daxLNe5Ics2wfE2SDx67Q1WdVVXPGpafl+SVSe5fswqPr+dW4gvn9tok+9twFfqUWHQOx1yn9DNJHljD+lbLHUneOHzz9xVJDi/4p8uZUFU/evTavaq6NPP9OtYHhJ4dz6nQszPfr8lU9ax+Hc+p0K/JOujZZffrSr6xuNxH5q9H2pfkC0k+muS5w/iWJH8wLP9UkrnMfyN1Lsm1Y9R6nNovT/L5zP+tfMcw9u+T/Myw/Owk/z3JgSSfTPLisWtexhz+Y5L7hj/7jyV5ydg1H2cO70nycJIjmb+26tokv5Tkl4btleSdwxznkmwZu+ZlzOG6Bf8d7kryUyPWqmenew5T3bProV875zEVPatfp34OU92vQ40z37OT6le33AYAgA7uHAgAAB0EZwAA6CA4AwBAB8EZAAA6CM4AANBBcAYAgA6CMwAAdPj/YYY4lmNo7H8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "f, axarr = plt.subplots(1, 3, sharey=True)\n", "f.set_size_inches(12, 4)\n", "\n", "plot_histo(axarr[0], binning, nom_bonly)\n", "plot_data(axarr[0], binning, my_observed_counts)\n", "axarr[0].set_xlim(binning[1:])\n", "\n", "plot_histo(axarr[1], binning, nom_sb)\n", "plot_data(axarr[1], binning, my_observed_counts)\n", "axarr[1].set_xlim(binning[1:])\n", "\n", "plot_histo(axarr[2], binning, bestfit_cts)\n", "plot_data(axarr[2], binning, my_observed_counts)\n", "axarr[2].set_xlim(binning[1:])\n", "\n", "plt.ylim(0, 300);" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "## DUMMY 2D thing\n", "\n", "\n", "def signal(m1, m2):\n", " massscale = 150.0\n", " minmass = 100.0\n", " countscale = 2000\n", "\n", " effective_mass = np.sqrt(m1**2 + m2**2)\n", " return [countscale * np.exp(-(effective_mass - minmass) / massscale), 0]\n", "\n", "\n", "def CLs(m1, m2):\n", " signal_counts = signal(m1, m2)\n", " pdf = uncorrelated_background(\n", " signal_counts, source['bindata']['bkg'], source['bindata']['bkgerr']\n", " )\n", " try:\n", " cls_obs, cls_exp_set = pyhf.infer.hypotest(\n", " 1.0, data, pdf, init_pars, par_bounds, return_expected_set=True\n", " )\n", " return cls_obs, cls_exp_set, True\n", " except AssertionError:\n", " print(f'fit failed for mass points ({m1}, {m2})')\n", " return None, None, False" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "nx, ny = 15, 15\n", "grid = grid_x, grid_y = np.mgrid[\n", " 100 : 1000 : complex(0, nx), 100 : 1000 : complex(0, ny)\n", "]\n", "X = grid.T.reshape(nx * ny, 2)\n", "results = [CLs(m1, m2) for m1, m2 in X]" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "X = np.array([x for x, (_, _, success) in zip(X, results) if success])\n", "yobs = np.array([obs for obs, exp, success in results if success]).flatten()\n", "yexp = [\n", " np.array([exp[i] for obs, exp, success in results if success]).flatten()\n", " for i in range(5)\n", "]" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD8CAYAAACfF6SlAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAClWElEQVR4nOydd5wb5fHGv7MrXXHv3cYFGzCYYnrvPTEdDAkBAqEEEgghAVIggRBCwo+SUBJaSELoJXEoppjeTLUBY9x77z7bdydpd35/7EpaSbvSSqezdaDn85F92jK7avO878w7z4iqUkUVVVRRxTcLxpa+gSqqqKKKKjY/qs6/iiqqqOIbiKrzr6KKKqr4BqLq/KuooooqvoGoOv8qqqiiim8gqs6/iiqqqOIbiILOX0QeEJHlIvKFZ1s3EXlZRGa4/3d1t4uI/FlEZorIZyIy2nPOWe7xM0TkrNZ5OVVUUUUVVYRBmJH/g8BRWduuAiao6nBggvsc4GhguPs4H7gbHLIArgX2BPYArk0SRhVVVFFFFZsfBZ2/qr4JrM7afBzwD/fvfwDHe7b/Ux28D3QRkb7AkcDLqrpaVdcAL5NLKFVUUUUVVWwmREo8r7eqLnH/Xgr0dv/uDyzwHLfQ3Ra0PQcicj7OrIGaduauvYe0A8CpQxYM7BJvuWUQgiuhFUn9HZHC91eoploK7M88Nthavn3FXdd5hVLMjZXwKs2Az7ac71f+47VoW0H2iq2bNyX/p1XKfZUL+a4tSAnvv/8Z6vm3FFsff9a8UlV7Fnk7GTji4Ha6arUV6thPPou9qKptdhBbqvNPQVVVRMqmEaGq9wD3AAzaoZP+9Mk9UAUbsDEDz+toNGXaAWzHZfldBbPIL5mJ5esAVSGB4XudLuYmX1tW6tjscxSjSAdkYhER28cxKxGsohy2oJjuu5ZpCSwMOkhzEXfmvhaf90uRvK+yk5F5HeezdO4w987I863wv69a/N8vASJFerLgb1jxBNA+DwFEEKQ49i0bBDAD3HwUE6OI+3JsmTmvRVWxsNEi37WIx5bZd8a8ok72warVFu+O9x2X5qCu35weLb3elkSpzn+ZiPRV1SVuWGe5u30RMNBz3AB32yLgoKztr4e5UFJ6yC4QoWqw63K2mdjUSCL1Q1cFEccpFQsLgwi2ayNtzwr8+cNaq13Oti7mJgzUdbHqObeQSyzuvgzJdbyF4HVYkrVNgQZNv8cdJZNscyHODznrvpI282G9XQukSSD9MjLfLyhluZqQQIioZn0viiUR710kLae3lTIa2qRKHeQ4RsNn2+aEpv7XFAGo66Y3EacDNUXZSn8vXFuq6e1FwkYxdMu+P20VpS71HAckV+ycBfzXs/177qqfvYB1bnjoReAIEenqJnqPcLflRfJHFDSyLgQLg2aNYKlgqzMTaLKjbLDraMh6FIaQcGlDFWx17GuRb+Faqx3rrHZssGqJ2ya2QtLxl/ZhiPMacV+jQkINGu0aGqx6Gqz6omxZGNhuwCg5e7J83v8GrUs9gqBJO5oc8eebjeVivV2beqTvIO1aS/tWOJ9bHAPbvS8baFaTDVraWEg9/5fq+JPnNwEJTQY5nRG3KVt+RXYCxSbp9BUbxXJf6QaNFWXLwsbGdpy+Jm2FC7Vkw/bYqqI4FPy2i8gjOKP2HiKyEGfVzh+Ax0XkXGAecKp7+PPAMcBMYBNwDoCqrhaR64EP3eOuU9XsJLLf1bFKGo+lYWPQrIV/PPkIIB1Savn9JKEYNGmNr6cIChcFQ0jkcVzZBNDRbMxry0aKyqwEzwgkRQAtRUPWbKClsBFiPp9lg+28jx2NRFH2yuV6FIgBMdeZdTa2vONPwnYJwA8bNEYHCT8DcGyV5vBb09Y3CQWdv6qeHrDrUJ9jFbg4wM4DwAPF3Jylwgarjg5moRBD6yKIGLLzDOVCdrioeDLIDy8Z5CeCEmwXFRoqHtkhodZCqSRQbqyzLTob5RlwtDaKJYBKhIUWPZNpq2hxwndzYIMVJizDZicJP1JoDUJoTTJoq0TwTSKBdbYzqm0LJJB0nG2dBL4JaBPOPywKkcTmIIdsQmhtMmgtIoDykkGSCKokUDqqs4AqyomvlfMvhCByaE1SaO3ZweaYFVRJII0tTQJtbRZQJYDKxTfK+Qdhc5NCa84OWmNWUCWBXFQCCVQJoIqWoKKdv42RcsxbIum7uUihtcig3ETQGvmBKgmUjrYyC2hLeQBLNfW+ft1R0c7fizBJ381FEH73Us5re8mg0omg3CQA5SWCbwoJVDoBQHUWUGmoaOdvqxC3jQD5glz4OeVkQZYi1Blxao1E0ZWvXluWOupCjhRA+r6yr12IDJK2LDUQlIhhYfjcV5hZgVOoZJBI2hIbI4/iRj4iULcYLuHWRkTExsxjq9BswNKkLcEU25GPKPD+B80G0oV1YAbIRwQhmwQcWQ6nDsFw5T5KLxJNF2U12CYdpDhZjXLAbxbgFFA56+DBKRgrRoqhNZAkANuVc0jel7RQvsKotiYpGhXt/BVotk1iGLQzi3faqpDQtI5Oox2l0Y4S8UgfhB2xq0JzqpDKkWaw1KCGhK/TzkcGqtBkR5NyaYCSsE1qjXheRwu5s4LkfdmpelcloSZRTRA1CpdqJYkgSQIxNUl4ip8SahJRi1qj8FQ4ezYQsw3ipN+zhBqY2NQS7rP0koClQjx1X4Ltlo/VaHGOdr1di6BEJfl6xK0wpWhbDtKSHOIqdmxUE1XoGOI9KzeSswBVJZ4lmGCjRFS3eMXwJo1hepy1jWIgREsooBQkw1YV4VHRzt+BUyEasw1qzeIUPS2P40//r1gqRFwnGzauH09VCWfaimuEWik81c+8TlJfKPP+mu0o9UYstANyiECzZA6cv+JEiGh4Ww4JKBGfkXkCk6jmn01k3JdVT0qrKMOWIxWR1CMKiwatde4rY6ukBOeKsYUrXpcpDOHoECUwiBapGpvh+N3/VZOzgC0TClpnO3XoUZ8PPwEYqltUC8fwEeWwUWzVomcmJkZZX4uNpD63rzvaCGUKsSI/kKSWTK76Szi5gQ1WXcbDUjPQVrGyIuJ7Xw6KlXbLp5RjFfnx5hO8K9ZWMAQrhNxG5hl5bJWg7uN/Rmm2IFfpNf3ceT+3hDPJ9w5vSRWcfPdlFU28lS3mJiJHicg0t7PhVT77txKRCW7Xw9dFZIBn33gRWSsiz2adc6iIfCIik0TkbRHZOmv/SSKiIrJboftrI87fwfp4Xd5HW4cibLKLT4j5/ZgVoVGjLb8pQDWpQFoWa0Xb2lzOSrXcziRtr8GOfGNGlKWisl15cRARE7gTp7vhSOB0ERmZddjNOM2vdgSuA2707PsTcKaP6buB76jqzsDDwK881+wIXApMDHOPbcL5p0fx+eElgoZEHQnb8BmVl9aww0ZybDnJZGGjnTlLKAQ/9+d9nj3rKGQrH4pRLs1V8k9jo13LWqudr0x1MVAVYhotUm3UPdfnebzoyKX/bC0pz51UES3qvnxsBWFzEYAFgUqXW9LJBo3tFS06dl+KBPRmxB7ATFWdraox4FGcTodejARedf9+zbtfVScADT52Fejk/t0ZWOzZdz1wE444bEFU/FAk/f0t/isb1wiGxDO8ho3hhJCsTHudogVW53j06b3w6zMQduWPPwHkvs789sT3R+DX/KTwElIn9GFmOY3sngXhlo264bUsW97XGH65qCOnHUmtW3Fgk9lnIOwS0WRsXzO+F5nkt96uDbU8NJm5ySGAPOdsjlyAjUMA2Z9ljZQ3Rl4KLDQntWuhbMT5rRazHNTCwtQttsy1h4h85Hl+j9uMCvy7F+6Zdf5k4ETgduAEoKOIdFfVVXmueR7wvIg0AuuBvQBEZDQwUFWfE5Gfhbn5inb+qsImq4Z6s9QfidBsR90lgYqtwar5QWGjNCmI20ksSQPZyeRgBNcFZLuIcD9Kr70OZlPGHaUt5rcVtIRUMUh4hJg1y/FnI3u1UCb88iu5tsKRgEMA6VYiufcVngSEuGsr+c75zQfD1Qj4k2+Yz7K1SSCBQwDJO3H6FtigW7YwzOnRkX7Hs9+9YuoBHFtW2eL/FkYxM7+Vqlowtp4HVwB3iMjZwJs4Ta8KLRH7CXCMqk50nfwtbuvbW4Czi7l4ZTt/BBuDjVbhL0J7M0iG1dGnL3WGWJgUSkOxdQGtbSdnCWmRP6bg2UB4O4UriJNkUoDYQpFAsmFN4fsrTAItczwNdqTVCCAow1IJhWH5fpLFFoRVYAgoqKthCqq6GGfkj4h0AE5S1bVBBkWkJ7CTqiZj+o8B44GOwA7A6+6srg8wTkTGqOpHvsaocOdfDFpGEMUjmxQqkQxaSgSlVheXo5q4HBXE5ZaOaM1q4S2xLLQSCCAf2nhF8IfAcBEZguP0xwJneA8QkR7AalW1gasp3O9kDdBZREao6nTgcGCqqq4DUv2EReR14Ip8jh++Rs4/DPwIolyE4DdDaAkhlMOJt9RGOWQm8oeFQtxDlQRaFZWuD9RWCUBVEyJyCU67WhN4QFWniMh1wEeqOg44CLhRRBQn7JNqhCUibwHbAh3cDornquqLIvID4CkRsXHI4Pul3mNFO3+77EvvchE0YygHKZRrdlCOWUG5iKBKAmmETQqXgtYMBfmhkmcBm1MYzlLJ25e6GKjq8zitbb3brvH8/STwZMC5+wdsfwZ4psB1DwpzfxXt/AE2JQp/4O0i5W+71hqk0BpksLmJoKWzga8bCbT2LKBKAGm01VlApaLinX8YFCKIcpJDNilsaTJo6aygJZLZLZkNVEkgHKphoExUCaB8+Fo4/0IIIodykEJrkUFLZwWlkkAp57ZkNvB1JIHqLKB1USWA8qCinb+lTlVtS2tSvHUuXlt+pFCIEIJsJVFMUjmfrVKIIGlPpGXOfINVl7LVMVIaEXhJwHtfQQhLAkHvWSkkkG2rQeu+drMA7wLIYn9GQTLRwBYvFMsmgKBq5mJhY5RUfd4WUdHOX4F1sVpHBz6g8rtdJB58vkLcNlLa9ABRsYiawV+UfLOEmGUQtyOpH1TUsKgxCssA+80OEuqI1SVtRbADbRUKD1kqxGwzZctAqfXYKoYIkn0GkrbWxtthil0yCZhYKVE4Q5U6yS9bHUQCqhDHyKiojqqVYysMCSR7Fnghqhjy9ZoFKLlyCpkKsOGwzrboKE7xX+o7pkqkhRr8LUUyEVyLmVX7XUUYVLTzT2pWWmogavvq5m9KFBIv80onQ1xNsPITgB82xmtyCozitjMiqjWL0233mx0kMFAb6kLY8pJBh0gTzXam4qiN0weh1odM8oWFnP4HmdLVipJQg4ZEHSLFzyIsTyG/DWzSGtrT7PtZepFNAnGnDQwZ7z8mov79FBqs+kACSDv+9ImO3k9agrpcJNCas4B8BJB2/JKx1YYSVPMhpnaGo7eBOEpUt+wswESwsCte4bMS0SaE3QAsu9RbzZVhjmvxnBckDx23TTbGa0KtSsoPwcIsWh56kxX1UaN0qprzVlD6CMf5S1YkSSD4vPBwnHe8CC0WR0yuPsfxJ5HI8xVusOpzpvDBWlHiK2xXrmV/xYrFhUE+pVB/AbXMzzIsgj6t4vVZyw/nW1F1/KWgwkf+SQTpp5QOvxlDvhBSOJuZBFBKQtkr6Rw2eRw08LKRvBr9SSQdeZ00k7/PQKatUpPLMaLE7Gjo5LCgqIrP6/TX5MlGS5PCbXcWUB6nKASP7rek86+6/JahTTh/R9JZaMxy2PUtdNbZKDchtJQMvOGh4KSxK00mudud60vohHFMo0Twzzvkm0eUSgJhl4oGO3gNRW6p67kdxtr5Ot/CdspJApsnGZwZ8iwV6ZBYrq0t6YBbg3gSarRYtrytoOKdf3Ka7ucAsskgG+UghyQhCOqTkE2ufAhrqyZlq86MF7WKKWhJaXJE7l0Vld3FLJkjKEQCqVBRli2RcD/yUlcYFSaBtNS0975AiEhx3Z+SXeGiPk3WwxJJOVYFtfYswCAZ+vESQBhZvFxYOI7CSwDJ5vC1xpaNHKe1Z6vzgGJR8c4fCIz3FkJjIoohNlHDUaS31WkibqtRNDEozoqapC0FEraZWn1SzAxBEZqsCFHDWa2iUDABmo3MWUEzhtipoZAjXe0jUVxw+agQJ4KJldKBtxAsjdCsUpRDd4hA6WCGd275SCChJoqdSvvaCAk1aXa7lRVTJ5AggqoQxXLeN3Blv0ObqPhZQC4B5G/5WQjNOM4i/b1wJKO3dD2A5Yq0t5nkZQWhop2/bYtvs5SibKhBs5Vrw2/WUIgQlOBewtkho0JkoBjE8iSxiwkRbbTSycSweYJgIhAsIr6i4sWHd6SkkJA/CTgJcSsgWVxssVjKVrKeYQvmA1pjFtAalcEJ95GNLU0ATiOeLZ16bnuoaOcP0BgrENqpKV/cvxRCCEKxZJB7fnpkXxwRFJ8wLragrJTwTvlIID9KrRguR1K4kmcBrY0tTQDlgo0RquXp1wFtfrbUGIvmfbTYfiLq+ygWmxLRjEdx59akHsVgo1WTeoSFtw9yGBS77LOUZaJhexB7UWrSriXVnQ1a1+KloaX0EC6EzdU3eJ1tpaqCqwAROUpEponITBG5ymf/IBF5TUQ+FZHPROQYd3t3d/sGEbkj65waEblHRKaLyFcicpK7/UIR+VxEJonI2z7N4nNQ8SP/liKIAFo6Y2jpyqNSZwalriBKEkAx2kNhE8VQ/Gxgc8wEqrOANDanQNzXZRbQEoiICdyJ03BlIfChiIxT1S89h/0KeFxV73ad9fPAYJwG7L/G6c61Q5bpXwLLVXWEiBhAN3f7w6r6V/faY3DaOh6V7x4r2vkX20qwGPiRQksIoZxkUEyIqNjwUEvCQmGXjBbj2CudBL6OuYAqAWwW7AHMVNXZACLyKHAc4HX+CnRy/+4MLAZQ1Y3A2yKytY/d7+M0ecHtALbS/Xu955j2hFgJW9HOH6A5kf8WayPl+yKXc5bQEjKoVCL4JpLA13UWUCWAsqCHiHhbJd6jqve4f/cHFnj2LQT2zDr/N8BLIvIjHId9WL6LiUgX98/rReQgYBZwiaouc/dfDFwO1ACHFLr5inf+hbA5yKEcs4RSyaClRNAaYaFKJIHWTgq3hASqBNB2CMBSKSYntVJVd2vB5U4HHlTV/xORvYF/icgO7ojeDxGcRvDvqurlInI5cDNwJoCq3gncKSJn4ISUzsp38Tbv/AshiBxaSgrZhNASMmhNImjN2UCpq4RagwQ218qgUkNBlRgGqhJAq2IRMNDzfIC7zYtzcePyqvqeiNThNGJfHmBzFbAJeNp9/oRrIxuPAncXusGKXu1j20JzLFK02JkXqmDZQiJhYNmSstWciOQ8wtqKW5m2gKJXGKVsJQzWN9eyKV7cKiLvyqGN8SgJW2hKmDQlTOyA9yvsiiFbIWEZrIvX0hByhVHQCiFbIWYbxCwj1Z+hmBU/3mMdWybNVoSEGr7fi7Arg5KqlxbCGqtdUd+xbME4RwZbXPlwyWur0GqgpFhastraz1Q5VwPlE4fzuzcL530r9ieZvQrIViWhNnG1scukxV9h+BAYLiJDRKQGGAuMyzpmPnAogIhsB9QBK4IMqtO04H84jd9xz/3SPX+459BjgRmFbrDiR/6q0ByLOBIDAfnfmqj/6EUV4onkiEPAUiwgGvHXr8k3S1CFmJW2lfzCRk0rpzq30KwgactR4xRQxbINakyr6BmBqlPR650VxGyTOjNBjek/e8w3G0jYgp2UdVbBVlgXryViaNEhobidqeCZsMHEpsawMxrOFDMTSLqduG1iYlNn+Mtk5AsFJZ2rA8EGVtvtMVC6FjkL6GA00pxSiXU+ywRQSyLw+xo0C3Ccffq+1LMn29TmDANl9wVIkpJQ3OgxSQDtRTKKCG0UQ5WoVPRYtCioakJELgFexBFGfUBVp4jIdcBHqjoO+Clwr4j8BOctPdt18IjIXJxkcI2IHA8c4a4UuhInPHQbDlGc417yEhE5DIgDaygQ8oE24Pzx0a7JRiwe9mU4thKWQTQSXhPGnxTStmoK2Momg6hppR1/ypZDCLWSdhph8gTp0WHmm+PIR8QKShZ4cwO24jr+zJNsNVC1iswLOE14sq9vYWCjmB7XVlyIJ23QwiChJlHxX1seFArK1fN3PksbYa3VrqgwUKNd42iLpm7LsRVXk5qA+0rdX1YuIFdS2/3uI77Ceps/DJT5YaqrtVvMmjwBEj4icTbObMDYwh3CbIwS5cpzoarP4yzf9G67xvP3l8C+AecODtg+DzjAZ/ulxd5fi6hWRH4iIlNE5AsReURE6txpzkS3sOExd8qDiNS6z2e6+wcXcaWW3GaOLdVyjTDE0ZUsctYatzKbr3iRz5RfkVkmiWRiYxFFZZsSNTT79gZwYHu2hykey6cCmrD9r1F8AZgQD/FZekNB6bvK1fNPwukhEK5ILNPxp7fa5A//pO4tJwzkf1/5TJU7DJSNfEObYgM2+T6tqkTD5kXJXlBE+gM/BnZT1R1wpjZjgZuAW1V1a5zpRzIhcS6wxt1+q3vcFkMsHvF9lGTLKi53EIQgx+uHJAkkNFdn3w+lVBaHQTAJlE7YxZCAXURfqga7jg1FOMrNJe1brsrgcqGYPEAVbRctHQJHgHoRiQDtgCU460ufdPf/Azje/fs49znu/kMlZP83VSf5m0gYRT1sm6JH5eUghTCJ5Oz7SiWi48Umjo1A0vDbmk9iIt9bZeTpuwu5JBB0tKrQbId7bWFJoLgZg9MWNPcO1Xe2UmgW4DfCVy1eJRQgpqbvfTl3XRitJQuRz0mUIg8dBKMqy7xZUbLzV9VFOGtM5+M4/XXAx8BaVU0GDhfiFDuAp+jB3b8O6J5tV0TOF5GPROQja/1GR5veBitmYsWKWy5m20nt8cxHIuEQSTGIxc0UmSQfYRFEBoVshVlBlEzGZb9GQwo7n1wSEF9blgqNVriVQmkSEDdXkGnLdsNUxeoHbQqhTxSWBGIacWSvk80LQjjYIAJIqs5mf5bNGi1aJyiOiZVaxeS9r/BftnJrAzXYEY8UtGY8SpWIjuP0A/A+TNji8f5vGkqe24lIV5zR/BBgLc6a07xaEmHgVsjdA1A7tL9asWRfW+eLUSwBUJO56sJrK4gAIr4JXCFhme6qI3V/6Omvf9CKIz8kCUBQaiLhhLC8BJC5ekiw1XRWTLgJOHXvNXV8gVVD2fUDlif1mO0WwxaPpZPDzSnn5ecuwjeaMdhg1RIRZ4VPPhRKICtCk0YxUUQ1Za+Q6/GvCxASmKmGktlusbi6AKGJCIa7+qVOwt2XH8q5GiiZCPYuPW1JbwAbpz9A8htqAc2qdK4A32+pEXpQ0tbRksDeYcAcVV0BICJP42Suu4hIxB3dewsbkkUPC90wUWecooVgqGDH/Z29hFytY8X8X6JZE+x085FC2ulnIjs0FIYMFMmYCYQtPPMnAv8G5JBeNRRm6WiSCAoVkYUnAWcUGrZgLD8JOA1cNlhm6OWhwcc53cEANtjOd6wYmYhsAvAPGjkolgBsdwa2Qc0WF4WVmwDK6Z+zf4Hf0GKwLYaWxPznA3uJSDs3dp8sOHgNONk95izgv+7f40ivPT0ZeDW5prUUaMIo+MiHZBipmHBSdk4hH0rJF5SSNC6muKwYSepiVgkVFw7Kj7DhoNbIBxQjG13MiiDILQwLfV6FJYJbG1VJ6M2HlsT8J+Ikbj8BPndt3YNThHC5iMzEienf755yP9Dd3X45kKNvnY1e9R24Y//jOGnwKIgbEDfoU9MJI2E6zwvdYxGkkE0GYQghDAnA5ieCUMd+jUggDMKSQLG9A4pdEVQqAbSEBKoEUIUfWvRJquq1wLVZm2fjyJlmH9sEnFKMfQNh++69WbppA8yYQtQweGfsBdgoCxvWs6BhLfPXr+PFuTN4a9FcBGgfrWFDPAbR4LCQlwDyhY+yCSAoVOQlAP98QRpeAgibJyg2NBScH/A5NmRFcVhdobBaQmELxsKEg4qtFA5zXDFiccVqBLVEH6jUMFA5C8I2hybQlgoB2SpFNT9qy6joxbxLN27g4IeTEwfBMA1+8cbLDOrUmUGdujCoS2eOHdqHBQ3reGvRXPq278i7Z1zIqsZNLGhYy7z1a5nfsI4XZk/ny9WuVlIWKWTPBFpKBpVKBGGE58LmBsqZF9gSJBD2uFIko4shACheJbSlCqHlygNsjuYw1RxA66KinT8KRsJJMdkRpdmyeGzq5zmHGSKgQlOTxR/ee5NBnTozsHNndu7Vj2OHbsuMNSv5cvVyRvfqx98OP54Za1cybfVKZqxfwbS1K5myahlNlvMlLoUM2sKMoJTZwOYmgbDaQYVWBoUd3RdDApU0C6gUAoDWnwVUCaD1UNnO34MkCWTDjmhKZG11UyN//fSDjP2miEMOtrAxHmPC/Fls07UHp4wYRYcaxymd8r+H+XDVAnbvNYAjBg1n+tqVTF+7khlrV7IpkenY/MigEAlAmggKkQBsPiLYEiRQabOA5LFbehZQCgFA6RLRVQKoorKdvw0SSzt9rcldHJSPFAAsVSxVUJi2fBVXvfoSABJV+nfoxPCuPZiycjkkDLbp1Iszt92Fukh6lLygYS3HP/cvVjVtYljnbtSZEWauXU2zncghgmJIANJLR/PVthRDBF65amjZbCBpS6T4vEA2CXhtlSMUlHzPwpJAPseetFWOXED2Z9mSMFCh70XqvBCzgCDxta8TAbRg4eA3FhXt/MWGyEawakEjmUSQhB8hgEsKrgC5WDi6iKLYUUWjoHFh4ZoGFq5pgKhj46Gpk3j4q8kM6tiZEV17MKJbD4Z27sbqJucH/IORezB2m52wbJvpa1YyacUSJq1czONzJ2Vc25sbCCICVYg1p48zDMWM5KpgehFEBKqO/n6y/sAQxTTtkmYDddE4sYSZEnIzxClEy1YaDUMC9WacmG24cgqO5EGNK4FdLAmoQtw2sNwFaoJSY9gYogVDQdmzAFVIaNKWIKpExCo5DJRQp9ArKekcwSLiymEUHwaqwxS3zE4hikVU/CXIU+cEEIBXIjpdnJUpD93WCUBVsdzixnLAkUcvT8JXRI4CbsepZ7tPVf+Qtf9W4GD3aTugl6p2cfeNB/YC3lbVb/nY/jPwfVXt4D6/ELgYp3xiA3B+VrP4HFS08wcQFcwmsKOAAXY082P2IwRwScHj+JO2jJgjF2HXeuzE0zbsqDJ3/Vrmrl/LS/NmulsNULjzk4m8OX8u23TryU49+3D04BHs0WcAj838DICf7bY/UcNg0srFTF65hMWb1vvOBhzJiswaSdsGO2ZimFp0aCiz7tJpeqJZstVhZwNN8cxlnbZCczxCbTSzUjoMCTRaWbZwpKbrzbStsCSQLqNKOjOh2RZqjQRGEbMAgFqJpRx/0pbTd8DCFC0qDFQvzWnHDyQrftE0AUDYWYDzGtWzAjuOCUoJ8tB+FdX+8tBtaSVQNgGU0/GXEyJiAncCh+PI3HwoIuO8DllVf+I5/kfALh4Tf8IhhAt8bO8GdM3a/LCq/tXdPwa4hQKKCxXv/JMQC9QAI57r7LMJAUCaXYefvR0BSzPCRckQEZBBBGTZXdCwjgUN63h+zvTUtm51dal59Tade7F/v8GcP9J5W5c3buCJWZ9x86Q3MkhAfSWNXe12LS5RnD4387mjo0NOo5mwJJBhC8VWwfQRdyumejjVA8E2iGY1mslPApnk5kXCNjKa1hROCCsJNX1lmBNqYLpONuwsIK4Rf1uYRMh8jwsRQJA8dByTqOYf/UMmAeS2f3FskUeDv5wrgTYHAahWpuN3sQcwU1VnA4jIozhyOEGj8dPxLJtX1Qluk/YMuKTyJ+AM4ATP8es9h7UnhNRvRTv/Pp06AM4XdWTvXqxs3Mjyho1k15n4EkIkz2vP2hWKCALsrW5qAneAe96Ep6gxTLbt2pOde/Zjpx59WR9zfoy1ZoSnjjyTFxdM57HZk1jRtCH4/jwongiSECcMFKDGWdySUSGWiORNEIcngWAZCiicFA5jq9gG80n4OcywuQBfez5x+1L6Bjv3Fk5LJ5kIbi/NIc/IRFsigI5bvvNXDxH5yPP8HlebDDxCli4WAnv6GRGRrXA00l4Ncc1LgHGquiRbFFlELsYpoK3BUVfOi4p2/pticeoAFP5w3BFs06cni9et59OFS5i0cAkT5y5g2vKVGeckiUESkDEbz4AgHv/izRtkJ5BTZGATXA8dN1LXiUUtPlu1lM9WLfUcYNC9tgNrmhu5fKcD+PGo/ZiwaAaPzvqUt5bOCd3DNLuauBAZWJaJZRVOFBebIG4pCVhqsClhBK4OCqsbpJDRESwbxZJAkDpPsc3kU/by+N6gWYAfYahK0fLQ8VSHs68vAWxQm/oyK4Fmt0QtgJWqulsZLjsWeFJV88b2RKQfTqHsQX77VfVO4E4ROQP4FQVaOW5x6syH9U3pL+Av//Myv3/+dSbPX8Iu/fvyyyMP4gd77YYRByMOFx+wF4duM4yete0wEmAmxCGA7N+zOiEkIy6pGYPEJOeRhJEQjIQgyQRyli0swEsYrgxFxgNYtHE9333xcQ586m/cN+UDdu05gL8fNJYduvZxTLl5ADterLyEFJSYDist0Zyn+bsXYWQkGhMRrICOXUkU+pGlJSMyk5ZexEM0c8mUifDvsqUKUclPpn7yEJZvxy4n6VsI2dIQfjMPd6EaG+zi5aE1p9FPYenqJMolCdGachCKs5qvQpEUskzCK3KZjbHAIyFs7gJsDcx0e/y2c+VysvEo6T4qgajokT8AFhgJmDp3GVPnLuPffApA704dqImYGDHo1q0dF+67BzURxxHMX72WSQuW8PhHn/PBooXYtSTDnRgJEDtZOBacQ8hOJAuCjY1GNP3rSYBYgkbIDBElkcwZeHSI5q1ezx/ef4tbJr3F/gO3YvKqJQBct+tR9Kht78wGlsxxvtg+BJC7ekiwLGe1UGpLwK87SQDBMwEhbplE3FU0hQZV+WcCQsw2ibpJ1CBbYWoEHALQnJlA0pWFqQ/wzgLUVc00PI7DxmCTXRuqOth7jGJgYWN6bFkYxImEWoOfOQNw5x4eW153XUw9gLry0DVqYXpIrdjeAJU+A2gGoqpEK68XwIfAcBEZguP0x+LE6TMgItviJG/fK2RQVZ8D+njO3eB2RkREhqvqDHfXscAMHxMZqGjnb1gQDQiNrljlxsyjsHr1Jvb8zZ2M7NeLnbfqx85b9WXvIYN4c8pszISwfY+eXHnUAXw8bzEfz1vE5IVL2KBx/L6P+QjBaMqdKGmNOjOMDBtJp59lx5NAjjXBhBnzSCYMVm7axDFbbcvRg7dhQcNaHpk+mSdmfs7KeEOGiSBCsH1G2UFhofy1A0LCzrxG6eEgIW5HSG4tpUbAa2tTorbFVcLJhLAivmPzMDH+nKWjGPi9Q2ErgzPzANmLMbNsFkEANgZN7kq1LV0M1poEEAfiFTYDUNWEiFwCvIgTgH5AVaeIyHXAR6o6zj10LPBotsKxiLwFbAt0EJGFwLmq+mKeS14iIofhvB1rKBDyAZBKLo5o13ugDh97eeD+gjPKGsFWZbch/fn5sQeybb+emIZBwrL5cslyrnzqBeauWpu2FxCFCLqO3yojv7qDvMlnj40aw+SIwVtzxsid2KfvVvxj6sdcO/EVIDVxKdjHwK+uIEyiOEwlcZgVQmE0hAolhQtVCjvH5HfyYaQiCuUCwsT4wxwTtioYiksEF1sV3BJJiHLVArRmDmDYwKUftzQG33GbPrrLXWeGOvatw25u8fW2JCp65I868fxsJJ10vu+RHQFizorpT6YvYuz0h2lfW8NOg/syeuv+jB7cn5WrNmLE4KJD9+KonUbw0bxFfDR3IR/PW8TS9RtS1wo7Q/ALF2mN+lYh+80OYlGLZ2dP49nZ0xjSuSvNCefCe/QewM37Hcuj7mxgRePG1DlhqozDSEsUDgmFSwyXIykcNhRUrllAEEqZBfihWG2g1pKFaKkqaKXPAKooDhXt/EXBjGWOmq0a8SWEJPIRQ2MixrvT5vHutHmpbQYwd8VqFq9ez7dGbcPY3XcEYMbSlYy5618YcejSro61m9I/mmxCSM4MjKxiMcgtQkvODHyXl3rOn7NuTeoOEzFY1LCen+96ID/ZZT9emj+DR6ZN5p0lcwPlqdsKCbQkFFRoVVDYXEBLCSDMccVoA1UJoIrNgYp2/tkjfzuaSwZehCGG7O+cHYEXP57Oix9PxxBhRL8ejB7Sn44d6jBcn3Hf+SfSp0tHPpm7iI8XLOKjeYv5aulyLFsDbfrNCsCfDPISQVT5ZPliTn/uMYZ27srYbXfi5BHbs1fvQez9xF3EbIuoYRC37RQRZJNAdiiokkigGLmIfCTQkllAoSWhxUhFl0shtEoAVbQ2Kjrm36H7QB115GUAWNH82fygeH0SVo1fEjfAVhYlHr/7SHYbNpDRQ/szsHtnAJ6b9BVXPP0CADsP7MuXS5YTS1i+Nv1yBtn5guxcgW+ewD2n1jQZ1rk7X65ejgATTj2XqWuW8/C0Sby7ZF66HjbLuQfpDFVSTqCl+YCW5gLach7gm54DqMb8i0Nlj/xtMGKKXSOY8Twj/mjhEb//jMGfELK/j+Pe+5L/fOhUZffq1J7RQ/uzsmETRsxZcvrID8YSSyT4bNGyVM7gk/mL2aDOTYUJEWWrl+abDTTHbb60nOY09ZEoE+bN4uQRO3Ds4G1Z0LCWp2dN4dHpk1myqaHgLAAqbybQ2rMACCaBYkXiSj2mNfIApcwAoDQSKKcgXCXBtiV0K9S2jooe+XfsMkB3POyyvMfYPiN6L4JmDH4jdL/ZQdCxSSdeF42w14hBjN66P7sO7s/I/r2JmAZXPjWecZOn0rtTB0b1783H8xazKp77w/TOCkqaDQBElVrT5MjBIzh5xA7s138rzpvwFK8unEXnmjpstdmgmT/wfLLTm2smUJ0FVMYMAEqfBZSDAMo1+i/HyL/98L468s/nhDr2o2NurI78WwsRMehKhPWxWE7lUtKpG3lyAHaNYMZsjJg6wnARsKOCVWPkzBSCZgcZeQR11vQLoCZYNRBLJHh9ymxenzIbgPqaKKOG9WXakhUYMThk2FCuOeFQAGYuX8XH8xbx0bxFvDxjJk1NCSKuBEuiPng2kJ0kziGBuNAcsxn31VeM++or+nRuz8rGTaAG5+2yOz/YfndenD+Dp2Z+wTsrZmOr5u09EI87a8MNQzFM/wKtsDMBVYgYDpkYhp0hNFdsPkAVLHXWwZuibIxHESltFqAKGxI1RA0LE9u3EC1sMlg1WekriCimR6Ct2ESwKk4Fsqu+aXhslTIDsDWpf5R5XznnlJgHKH0GkP4ON9gmHY1q0/bNjYoe+W+7zQ761bQvaG6Os3LVBs67+EEaG2PsvvcwhgzqwcrVG1ixqoFVqzewYtUGmmOOI0qN9m0lutHOqHBXE2KdzQy5y1CzA1XMZjKr5QWauwhqekM0mTbMOpPt+/di1yH92XVIf3bZqh8damvZ97d307CxmQNHDqVDbQ1vfTWH1WYzagbPBnxrCExFrFzpCa21wYTtu/fitO1GMWbodnSprWfpxgYemzWZ2z57K32PHgKwbUAz3w/TtDHz1SrgTwJ+Xy1T7AypaS8K1whkGjREqTGsvASQhJcEkpIJGbZQ2pnxQOcYTALqkZt2IECdkWkr1CxBmog5Ljpjew2JDNIMNwNQaiWeI3xXIxZGgNgfbK4ZQG4Zm0LJBJC0NbQ68i8KFT3yX7tmI3ff8iLde3aka88OxNY1QY3B/rtvzZhjds44tqkpzlHH34Jda/Cdb+3G0EE9WLmigVUrNrBq1QaWr1jPtOlLwYJog4VVn18TJjuPIBaOLpDnGFWoWa/EO7gjcp8loFaTxWezlvDZrCXcH/0IQ4TB3brSsLEZAb49eluO2nkb4pbFx7MX8cqsWbw6fTaL1q3PWDXklxcA75LRbO1mA623mbJqOde8PYHfvfs6hw4dyklb78BW7buiCQOJ2By71Xa8s2Qua2ONGFHLdfxZsyzLwFYlGg2fE0g7/ixbamDamiFHkURjLFqAALJWULkzgYhoUbmATDft2gJitkmt6e+AgmYBaafv+WzQHFth8gCbtMaVYsj4lhHHpNZTjxxmBiAolho5ZBZTk1oSW3AG4HYZ8PH+pcwAcr+tLYMiGYq3X2dU9Mi/U/v+utf2Ti8DBewO6aF4u3Y19OjRke7dO9CjR0fq62t4ZvwkAH500aHst/cIundrT8TV+1m4aDVnnnMvAKefticLVq7jvQ9n0UiuQ/PLI9ji/yVTwKoDRDJyBtl5gpwVP+6hIjBqYB8O3n4YB+8wjGG9u/PRvEV855+P+56bkxeI5jof574UNRVqsz7faDIIAAM7dOatky8kZlm8uGAaf5r8Ogs2rvVxDIqIklTQDZMTiKbi/X5jPKgrkDMIUykMYGBTG8l0GPlnAUqdGfe5L2efIflzAdkEECxQrdQb/jOJIAKIEKTZr74OOx8BmFmzBa+tQqN/aM0ZgL/zT8/GpKgcgNf5DynDyL/d8H669S3nhTr28zHXt+mRf5ty/tQEj9btmkzdHQWICJ27tKd7tw7U1JhM/coRUfvrHd9jmxF9adjQxKtvf8WLr03hi68WA/4hILvGHd/5/JgUJ5SUTRjZyWMvGdiRYFsD+nehY30tXyxeRtd29Txwzkn8Y+In/PezqcTNrARwVFETX23WlPNP7suWonBH8dt17cnJW49i7IidiBgG9301kbu/fDercX2m808iHwmI2JiGX4zZ/Ym721uaFPZz/kn4k4BSawSNfDXlMMMmg0tx/uBPAMU6/yT8SKClzh9aiwAKO38InwSuOv/SUdGSzl4IIDEr9ciGEbMzHgKo7YSOZs1ennL8qHLRj//FFb94jPc/mMWRB23PXTd9h++fsg8AZlwzHo5tBctHN1ndXikiKWnpZKjIjGnqAWTut/HVJxaFBUvX8uXcZRgx6N6+HYmEzY1jjuT5i87ihJHbEbUkFVoy4k5bSl9bgGF5fmFxyRSac+Wmp65ZwfUfvsohz9zLc3O/4oytd6HG8Jn2+jgSR1La/yukGuCpslBoil1o2Z2d5yvsLxntELmfDLMXySWhQfDKRGebUhwXl09o0k8e2k9qOjBM4kG2NLRjywiU+Q6r7JlcClosyiEH3ZpS0FU4qPCRfz/da+QFzi/VKMxTmjUzUAGr1sxxXGKlE0719TXse/C2fDV9CfPmr2KnUQM57+wDGP/6FF57ZxobNjqjGAXi7fxiP5KxKXvmkB3+sdxZhFVLXur1nnfwdkP54RF7M7JvL2avWM1J9z/MprjDMlbEqTLOyd5Y5Iz2c1YJ+cwGurarZW1iEwLcuMcxPD57Mp+sWAQqeZeHQu5MIDn6Tz/3/p17fr5ZgCE2tdHg6xezLNRp/p6+VvK+/Jxs4RlAo68zTfbPLZTozdyvmD4ziVqJB4ziM5E5A0jaypJHwaCDWdzqnPLPAPwTvn4jjDAzgHImfL9JI//KplcFiVvOEksfAV67NtOz5swIakzMJgs1BTUEsTXD8ds1Bo2NMV553mnAbgLtTJOOHer42cVHcukPDuXdiTN58ZUveP+zubDJCbPYpmNLJbMZdnYxWnbS2Luc1IxBok5Q0xnxq4dXnERv+rzXps7m9a9mc8jIYYwa3JemTXEMYKcBfZi0bCmaANtQbHewmlzWSqGloh4JCee5wZp1cagxGNipMwf1G8ZpW+/Mc/OmctOnr7Ngw1rn3gNIILtYTNUgYWkqxGCr0zQ9aHlocyISSAC2GjTGhPoa//1hisOSBOA0f484yztRbFuoNf3DKoWLwurpFHWWVTrdjjMT8IWWembuFywMlzY0NYNIhEzCZvcG8NoCXFqRzSYHEZwETr66zG2+1w4hA1HO4autEqrx0dcBFT3y71zfV/cZfDZaE/7DyCaEbGTPDlLnZeUMhg/vwxFH7MDBh45EEE7+7p0kEjZdu7RjVWPuDyFnxJ+dA8gzIyiYH/DZvlWPLvzv8rOYvWI1d7z2Pi/OmJH6EYQtHCs0E6ivNzl/+z24YIc9MA2Df077mFsnv0Wj5TjRQjMBCM4LFCoSa0kuIB8J5EsGt6QorKUFYeWUhKgkOYjNWQRWjiKvuq376+A/XhDq2GknXZv3eiJyFHA7zrjyPlX9g88xpwK/weGwyap6hrt9PLAX8Laqfstz/L+B3XB0+z8ALlDVuIh8B7gSh0UbgItUdXK++28Tzj8IYUihpWRgmgYDBnZj3tyViMC//n0Rzc0J3np3Ou++P5Mv5y3Lia2WiwjyrhiqEY4cNZwfHrY3Q3t146ulK7jztfd4aeYs3+OLrh72HN+rUzuuGL0/u/bqz9HP3UfcDqcZ5EUpJNBaBAClk0BLqoI3Z0XwN5EAKsn5i4gJTAcOx2ne/iFwuqp+6TlmOPA4cIiqrhGRXqq63N13KNAOx7l7nf8xwAvu04eBN1X1bhHZB5jq2jka+I2q+jaMT6KyE74FGspKLJHzyIbRHM95ZNqwfBPJycSxNiZYMN3R0TEMg4f//S5r12zkjFP34u7bv8czD1zEMQdsj9FsYzQ7Ds4vYZx8+O4PSBRnJ5GNhKeOIKa8+PF0TvjjP/n5oy9QZ0a4+eRj6BapS5/nOT7Zszi7b3Hq9bq9ilPwJIiXr9/Ez995gWPGPUgsBnXU8vBhZ3D4gOEAoXoOByWG802xmxORvAnhfMngxkQ0VSHsh3z9g5NS0X4oJhGcDb8kb7HHJCuCC8EvCexrzyquN3Ap+IYmgPcAZqrqbFWN4fTVPS7rmB8Ad6rqGoCk43f/noAzgs+Aqj6vLnBG/gPc7e8m7QDvJ7fnQ+W/o80BI7SAEb0fAWTPELIJwDs7yM4bJGcGRsxGsRn/zKc8/9xkOnasY/fdh7L3PsNZt9Jp/DJkqx5ccO5BvDtxJu9NnMWyho25OQCPhIQ3R5AvP+CVmMgtJFNe+PArXvp4GiOH9GH9uiYM4J7zTmT+6rU8PHEy09eu8hyfv3AsR1AuLs4sIG7QHLchCv3ad6RHXXvuOehk3l82jxs+fpUvVi/NIIBiBOQKSUXkE4xriVCcNw+Quy9YGiKMRHQ+eeiWSkO3hix0WGx+GQjPtTeHDLQGD1R80ENEPvI8v0dV73H/7g8s8OxbCGSPxEcAiMg7OKGh36jq+DAXFpEocCZwqc/uc0nPDgJR+c4/CEWQQiFCKJYMjJjNxlWbeH38F7w+/gvsGgMT6NGxHQP7d+PyHx0JP4IZM5fx7sSZPPnMR6yLx/ISgRdeIsjVHMoMCyV/C1ZE+XzWEgwcSYmV6zZy4i7bc/oeOzFxzgIenjiZV2bNIpFwHG82CTivNTdBnCIASJHArHWrOeq/D3D6iJ24bOf9+N8x5/DM7C+4euILNFuuxEYe7aBEwvANA8XikYJaQUGhoHzVwYUIAPzDQIUIAEpTB600AtjyCeA2i5UtDDNFgOHAQTgj9TdFZJSqrg1x7l04IZ+3vBtF5GAc579fIQOVHfZRhVjM/xGE5rj/Iwv5Qkb5QkX5QkST3p/N9878K+ecfQ/3/O1VmjbFOO2kPUgkLIxmm71GDWa/XYZSWxMJDA1BZlgoqH4AMmsHkiEeI+FISvzq4Rc55MZ7ueWFtxjQpRO3j/0WJ2y/nW8IyRsOcl5jOiSUEQ5KHhM3sGLCQ9MmcfDT93DXZ+/RpaY+5fgNz7KZoJBQUI1ALB4pORTUGmGgfCEgyB8GakshoGLDP1uqBqANhX8WAQM9zwe427xYCIxT1biqzsHJEQwvZFhErgV6Apdnbd8RuA84TlVX+Z2bcXwlJ3x7dRqgd930IH0G9WDx3BWsWLyWlUvWsmLxGmZ9sTD4xJr8P9gM5EkI50so+yWS/ZLHNZ1qaWpyPO4f/zSWXXcdQnNznE8mzePdiTN559M5rFy9IXW8NxnsTRRnbPepHcjel50cPmDbIXwwawEbiHPy6B3Ye9ggHvp4Eh8vWJx7POnZQODqoGRCOEvvp0+7jjx+9On8dcr7PD5rMpbn+1VsM5mWrAja3KuBtuRKoHImgdtyArgsCd9h/XXAjReFOnbWab/Ol/CN4DjzQ3Gc/ofAGao6xXPMUThJ4LNEpAfwKbBz0nGLyEHAFVkJ3/OA7wOHqmqjZ/sg4FXge6r6bpj7r2jnv+3wkfrVDCc5rqqIO6KMNcUZ/8i7rFq6nl32H0Hn7h1YvmgNi+euZOn8lSyavYKPX5/qbzQMMZRACNlk4EcERrsIo0YNZJ99hrPX3lvTr19XPvhoNj//1RMA9OjfmeUrnRxPEAnk7CuwUshvqej39tuFiw7di071dXy1dAX//mgy//tiKo3xRChF0RQJeFcQuSQwtFM3btr3aHbvPYCv1iznxk9e480lszPuv5wkUCUAB1UCqCznD6mVObfhxPMfUNUbROQ64CNVHSeOQ/s/4CicsswbVPVR99y3gG2BDsAq4FxVfVFEEsA80sngp1X1OhG5DzjJ3QeQKPReVLTz79muj/7w2CvpMagnPfp3o++QnvQZ2J0uvTrRoXM7OnT2n842N8WY+tEcVi5dx877jqC2Lsqa5etYMG0x0z5bwFeTF/D5ezPTJ4QlhHgCGhvBsh1J6Lo6tKNPaX2IWYFtwNA+XagxDWZ+tYSOA7rw+H9/wrz5K3lj4gzeeHc6s+auAELMBmzFSDghomRvgHh7A0QCSaA+GuHYnbdl7D47sV3fXrw1Yy7nPfqMW1gHYqfPS7RPawQlSaDQLOCorUZw1a4HMbhTV15bOIvzXnkKCxvc/gCltJQ0DAvbdvsUi2Ka6f4AxZBAtqRzUhraW41cDAF47QlKrWERNfxfR6GloLYKcTVSxWIRrFSvgWIJQCGjXlhQuhqb8kpFpGxlEUCyz4AX3l4DxRGApnooJO1EsUPdVzY6SAJXcJdtB7Xc+dcOHaD9rrs41LFzz/xFm67wrWjn38norntFjnSeRCKI+0s36p34ZG27Grr36UL3Pl3o0bcL3ft2oe+QnvTs15X2nerp3rcLPfp1wTQzHe/GhkY+fWsaC2YsY99jdmL1svXM/Hw+c79awoJZy1gwYxmNG5szScGyoNlntFJfB3XpOKbfzCBnVhARjE2Oc0l+39t1rOPQ43dh30NHMmrUQEzTYNHiNfz+zy/w+VQnVOhLAqqI66eSexXHaTd3S7/ufDOBXbbqhx2FT+ctpmtdPTeefASPf/A5b0yb44RtBGKdNXWBsLOAqBh8b9vR9GnfiRs+fA1QTBEsMy1gVnpfYed60UjaVlgCSK8g9nobpc5MZBAABJNAkgACbRmJEghAiYiVY8vETtkKSwCZjt/7zYDu5sa8NlK2XALIdPyZtoonACWSrDWW9DYBagKF7YJRL4mUNMS2g5ZUnX8RaDvOXwSJFpHscV+XiNC5Z0d6DejOwG36MmzUILYaOYA+Q3rRd3BPTFfy2RtWmvHZfJ7959ssmr2Cg44bzYKZy1gwbTFffTzHIYVs1NdDXe7sIYgIJGGB7a9vkuhSR5du7dln3+Hsf+C23Hz7eJavaGDfA7dh1x0H8fq705n0xYIMIlD8hZPjHZ1GM2HDQQrsNKQvt449lj6dOzJnxWr+8sp7vPDFdKw63E5ouaEg31lAynmnr71jjz7cedAYrv9wAi8vnp7a3hICMEQzjitEANnqkdm26n3OzzcLqDNjPrYcex3zhIj8CMDAChBxU2olLT8RhgCsHGedtmWgdC2iCMzWtFaRny3v/RYiAEk6f5/XGPXMBsJAUGo8tqrOvzi0HeevrtJJbcjVAoVelwi1nTvSd0hPBm3Tl6226cvwnQez1TZ96d6nM1F3tG5bNobphBoa1mzkP/e9xrj7XmPDOs+0uL4+9xfrG/rxeFyfIY4Cdn0UTCNHbuLU7+zNWd/Zl7q6KJOnLOTP905gxpzlgV3IHPE4QSPp0E8hErBqAANMQzhs5NZcdPBejOjTg3dmzOOiR/9Lo0c6OW9COKqOYJEj5ZPCDt1786f9jma7br14fOZkfvPhyympCPAngUIEINi+TWaCSUCpjSYIctjto0Ejff/ttYGyzY677FAEAZh5JJ1rJFOiuRABtDPc/qA+tpJBpbA1AO2MpkBb4pHBTiIfARjYge0kTWyiUojwvbYyw0VV518cWuT8RaQLztKiHXD8zfeBacBjwGBgLnCqW3IsODoXxwCbgLNV9ZN89lPOP3WPhe81RQ4hnH82kuEkwxB6D+rOyD2GcdLFRzBkZH9WLV3L2pUNDNthIFcc939MmZiWUcA0HXtBuYOcsI/bRjLrHhTQqInWZs4YkkRQVxfl4CN34Nyz9qdzp3b884n3uO/Rd3D00nxsRcCqdc4Nkxi2DdAoqd+5IcIpu49imz49uPaFVyFEi8kUARgKNZrjMyIi/Hjnfblkp72ZuW4ll7z1H2asW5na7z8LUCIBbSQNsQMJwp8AlJpIsJM1JDgZ7EcA+Zx/mN4AXgIwsAI1+L0j/yTyEUBUEkRzQkjufXnk58IQQISEq4JaeOSfRBAB5Bv5R7CJbOmR/5AB2vfaH4U6dt45V7Vp59/Sdf63A+NVdVtgJ2AqcBUwQVWHAxPc5wBH46xhHQ6cD9xd3KXCfSm0uRltbs4vgB4Au7HRedjKkrkrmfD4RC4+6Hdc9727WbVkLcN2GMjqZevYetQg6trVcsF1J/OTW89k0Ii+joGgGoTsWgPLx8GpIqpOeCugjqCpKc4L//2Us878G0/95yPmzVqB2WRjmgZmltcQAENy6gUgt1YAnDX/KR/lvtW2Ko9N/IzfjnsVw4JtunbnjuO/Re+69u65mXIRjh23LiDZsyDjY1MSanPLp29x5kuP0bW2PiURkXprfOsCBNv253NHKdQf/vUAgmUH6+YDgfUAfrUAcfXTzc/cELYWQH01+F1pZp+vc746gISaPj0VcgNeYWoAEiRteW+utEFjahjnOV0VVAWzSJvJzskVHLyoaJQ88heRzsAkYKh6jIjINOAgVV0iIn2B11V1GxH5m/v3I9nHBV2jk3Rzwz4tmJ34jcZDEkNyJpDE6ENGMvayoxm193A2NTSyeM4KBo3oS01dlPdf+pwn757AlA89yxrzzARUBCK5vQaS8CaJc1YKeUJCVo3BaafswVGHj+K2v7/KR5/Nc94uDxn4rRDymwnYUdweCDgdwlIGnNs8avQIbjjhSBrjca4b/yrPT5memgn4zQIUReudRvLuBmdG4N5Ot9p61saaUNNiZNfezG1YndFBLHMW4PT9ldwJUwrhl4QqppHsNJa+7+xRd9jloKZYRNxwRam9AbydwQyvTDhCOyOW9ysbNAMwsN1wkec1+mjpQ+EZgLi2TE9YJmjUn0Rw+MchNCP1DBJujLD4JaDOjMEEttuqOvIvBi0Z+Q8BVgB/F5FPReQ+EWkP9PY49KVAb/dvP62L/tlGReR8EflIRD6K00xLHD+AxmKph3uB0OcmZwJJfPLql/x8zP/xk6Nv4r0XJrPVNo7jXzZ/FTvuPZybn7mMUy8+PG0gz0xAmmJIPAEx5xFUZQwEzgQAzJjNgpkrqI1GuP3aU7n5yhPZfkifvNXD4F81bMSd0X+kCSIbnYfZDKZ73PhPpnPCXQ8xf9Vabj3xWG4/6Vi617ihMh/ROEEwGk1oEudhCyScDmIAq5sbsVWp0Rr+fvCp/O+Yc9iuS6/U68ycBQi2bWBZQiLh/xkWqgxOQ7Bsk1jCpDkecbqe+ZgMWxFsqUmzHaHZjtBkRf0iekC4GYA6wZ/UQzHYWKCSN2gGYGPQpDVssmswUCfWHmCj0AxAMWjWGveONLX8NB+Cq4AFC5M4BnEMEpgkRwTFVwALCUyaK1ysoBLRkncsAowG7lbVXYCNpEM8ALgzgqK8t6reo6q7qepuUVquBphhOxZzQkJFIpsEvvpoDjdf/CBn7nQ19//2aVSVdh3r2Li+kR69OtKjWzuG7TCAI0/fm2hNJD8JNMczfpClkMDEibP4/tn3cO+9rzFy237c839nct5p+wL4kkA+RdHkNjPu/BzTEhDO//OXrOHMux/j/15+m0NGDOWM3XYKlIpIh4IMjISR2Uoynv7qNVsJfvzmONpHanjm6LM4ffjOGW9TZhjIWXWSTx4iCLlhIMdWUzwaKA0RXhJCUo9SVUHTISDxPBwUknrIv1/YEEIKIkwIaINdX1RENb8MhN8KolIlIIoP8/pCSbU4Lfho42hJ2KcP8L6qDnaf74/j/LemTGGf/t0G6k9P+RUzPpnDvCkLSMQL68YX/TrCrh7KgjckZBjCrodsz7HnHMDuh+2AKiydt5L+Q3uxetk6/vPAGzz/r3fYuN4lkJCJYchdLlooHNSuXQ3HH78rk6Yu5PMpi+jWrT09+3Vm6oylgCfsU0A6whsKSu2LZB4zpH835q9eR9yyGD6gB0vWNbC+qTlvKChfXUD3unbcuv+3OKD/EP4390uuePdZYnbmZx52RVC5q4JLqQhujWrg1q4Ebo0+AJurAnjnrRa2POwzeID2/eWloY6dd/7P23TYp6Wrfd4CzlPVaSLyG6C9u2uVqv5BRK4Cuqnqz0XkWOASnNU+ewJ/VtU98tnfbvhInerKO9x83t289I836DmwO7sdsRMzPpnN3C/KRwjlIAGA3oO6c/T39ueo7+xL5x4dadrUTF27WjZtaOKJu17h0T+/lD64FUkgiR/88BDGnrwn730wiwcffzcUCQTlA8C/PkDqhOd+dDZ10Qi/evZl3pw113dVUJjqYAEuHLUnI7v15kdvjEN8HHsxS0KDSKBKAMFoqwRQdf7FoaXOf2ecpZ41wGzgHJxQ0uPAIBydiVNVdbW71PMOHB2LTcA5qvqRn90kOhvd9IQRZzJ89FCmvDuNlYtWc/iZB/CzB34IQKw5ztwv5jPjkzn8+4anWblodcmvJeN1lYEIojUR9vv2aI79/oFsv8cwrITFnC8XcfuVjzHri4X0H9qThbOWp08Ou0y0SBKor3dmAqectiedO9VnkEAh2Ygw9QHJ5yP79+L3px7J8N49ePyTz/nDy2+yMRYreRaQxFYdu3DAwME8ND13VXA2CZRTG6gtE0A5dIDaIgFUnX9xqOwiL+mme0WOyNned2hvho8ewvDRQxmx61C23mUw525/OWtXrOeky47lkNP3Zcanc5jxyRxmfDKbOZ8vIB6k/58HpZIAZBLBkO37c+zZB3LIKXtQ376ORbOW0X9Ybz54dQpP3PkKX3hrBsCfCArMBoohgYkfzuL3f3oO260BKJYEgmYBUdPk4qP25vv77sqy9Rs4/R+PsXT9hoxjA2cBAQRw1a4HceGoPRk/bxpXTnye9bGswqgtQABQXC2As71KANC6BFB1/sWhTTr/fDh47L4cedaBbD16KJ26dQBgU0MjJ/Y8F9uy2e2InaipizJv6iKWzF6GbYWrKCwHEbTrWMehp+7FmHMPYsDwPti2YhjCsoWreeO/H/PvW8cTa/I4jxJmA2FIoK4uwsqNTWw9tBfnff9A/v7Yu0ydvsSXBIoJBSWf7ziwDyfsuQO/+d8rji5MznEhwkAAUUcD5tztd+fKXQ9k6cYGfvTGOCavzZRFb808QJUAyksAFe/8txqofa8M6fwv/lnV+bcWSnH+XvQZ3JPho4fSrW8X/nvniwDcPOEadjxgJOCEjRbNWMKk16Zw9+X/AKDngO6sWbY2by6hHEQwap/hfPu8g9nn6J0xIwaJeIKn73mNt5+bRL/BPVm9Yj1TPpidJqcQs4FiSABg332H89MrjqFz53a8/9HsHBIoZhYAuSQA0KN7e/5y+hhufPkNPlm4uOQw0C49+/GXA8fQu10Hvv/Kk7y9vLBMdJUAvlnx/0pz/q5e/+04lS73qeofsvafDfyJdJOXO1T1PnffeGAv4O0sPf8hOP2AuwMfA2eqakxELgfOAxI4S/C/r6rz8t1/xTv/Pc3DU4JrpSD1+tzq2foOdQzarj9bbTeAQdsNYKvt+rN62Vpuu/BeAB786nZ6b9WDxTOXMm/qIuZPXcjkN77k01e/yLEFjnBcqWQgNTV06tKOfb89mv2O342d9tuWSNTESliYEZMN6zbxzvOTeeeFyUx6ezrx5BLQbCKojToy04kECGi7+tTC9UIkUFcXZcxJuzH25D3o3Lkdr779Fb+5aRxGXLEjThGZXZcZHiomIbxd317cdua36NelE/e+8yF/efM94rZdXBgIIGrTKVrLj3bah//76C2aNIH4aPqEDQNFIwlUxf0oFcPNkZdCAHVmHBtBVVytG2cNfKkEYKtTuSySKZ9QKgHUStyRwA6oFA5LAKpgIc4SYJ8Cr2IIoANNznsGgVo/2QgigOR68tEV5PxFxMRp5nI4Tk3ThziNW770HHM2sJuqXuJz/qFAO+CCLOf/OI6G/6Mi8ldgsqre7bZvnKiqm0TkIpyVlKflu//Kd/5yqPPEMIomAVUF2+fHn8fWQafuzeDtBzJouwEM2rY//Yf34cUHX+fWC+/BAP42+WaWzFrG1InT+ezNqUz/eLbvLCEfIagqJKzMunQROvXvwd7HjuaA40azy4HbEYlGUqGht5+fxI0X/h1VqK2P0rQpliaBWMxx/F50aI+2T6+xDksC2MrjD70HwJChPZkzewXxjiZWnRkYCgqcBdggNrSrjfLzEw7kxL1H8eWiZVwxbjyzVq4OJADn/CwS0NQ/zsuL1vC7/Y/gpsmvsrSxAS/CEoDXnogSMR3nWAoBZNhCqXVlq4sjgKR2ZuZvstZIi7oFEUChEX7y/uqNmK9yZiECaG80YWWVBUXdXgMZdkIQgCTF3TyvM4IVStcnmwAUR0kEYNfKcv574zRkP9J9fjWAqt7oOeZsApy/u/8gPJ283EUzK4A+qprIvobnvF1wZhH75rv/NtMQ0zvaDg0/x5/cbmbrxzh4/fH3gPdSzyNRk7r2dWDb1HWsZ87n8xm641bsPcb5zGNNMe687EFeuP9VzIhJtCZC06Zm32KyFCHYmitIosr6hSt48YEJvPLoe7TrWMceh49i/+N2ZbdDt2e/Y3bm3x9fz5QPZrHH4aP4+I2pvPPcZN5/+XM2bPJxVhs2IhFHcE5rIqlCMbs2mioUS5KAEbOJxZp5/KF3U1IOe+41jBtuOo3nn53E/fe+zlIrjhFzQkFmXDOazKcdX5oEjAQpqenG5ji/ffQV3pwyh2vGHsZ5u+/K1S+8jJEgs5F80lqNZjaPj2TaBxjetQeHDhrGvv234odvPs1HK9JtPa2YmUEAyWKwXBJI21MFyxYipgY2ik8WgvmTgMcWELcNakybTYlo0c3hMwuWlJhtUmc6ryeoMXzhpvCOzSa7hnZGc85PKV9DeEFJqJlzThwTQzMF5wo3g/dWGqdPTLi2/IXtgiwlvzNlKvByjRoBFeQ+6CEi3hWL96jqPe7ffooGe/rYOElEDsCZJfxEVRf4HJNEd2Ctqia/nL4qCTgN3F8odPNtZ+QPSIDD9kPgqL8EWwCaJcbWsVsHdthvW3Y8YCRvPjWRrz6Yyc4Hbc/vn7+aGZ/M4cUHX+eF+1/F//3NI1ATjWTMSoz6eura1zoE8K1d2PPIHalrV5uSmrYSFpPfmc6dVz3K4jkrMm3V1kI7j3BYgZyACqjhTOvr62v47ln7ctIpe9DUGOfu+1/jf699kZMPCJoFqDiqotno2rEdzWKxVprp16UTGzXOmk2NBcJA6jtMGdalG/ceeQIDOnTmV++/xBNzJmXsL14iWqmJps8ptUVk0ladmXaMYWYAEii/4DSH8X5lWjIDCBr9g/8MICKJgNCM+o7+IXgGYLjaRbmvU0NLOidH/9nOf9etFrR85D9ooPa/4rJQx8659Ip8I/+TgaNU9Tz3+ZnAnt5Rvoh0BzaoarOIXACcpqqHePYfRObIvwdOYe3W7vOBwAuquoPnnO/i1FMdqKp5EyVtv0Z5C6Fh9QbeG/cRf7vin0ydOAOA5QtX8fjN/yMSNbns7h9w86vXMnCbfsUZzvod2Y2NbFq5ljcfe4c/nH8/p21zBb89825ee+oDEnGL5qY4XXt1YvWy9T62NENR1CsdESQbIS5ZNTbGuPevr/GDc+5jxoyl/Oyyo/npRYdnSEUAOaqhXsVQP2GP1Q2b2LCxGSMGt512LE99/3S26d0jQx7CuaekNIQglvjamrV2Ncf/51+8t3Qef9z3aL4zbHTGfj+F0CBdID80JyIByqAEykEEwU8R1NkeooWoD7xqoF4UkoEoNEpuCCEDEQYNVn3hg0pEafIPmx2LgIGe5wNIJ3YBUNVVHgd9H7BrAZurgC5uc/gcmyJyGPBLYEwhxw9tyfkXGfJpSZK4WCSvtXjmUh685jEu3vMX/Oncu9lq5ACue+ZnGDlzWZ+wDzjb4sHJQLuxkaY163l//GRuvvhBfnLUTSxfsIoh2/Xn/N+eROduHbjmwQsYuftQ5wTvbMVT55BPOygbC+av4meXP8w997/Ohx84K2xSDj9AMA7AbM4zo3R3/e6ZVzENg0fPHstR2w33JQCJ+Tv+JNbHmzj3laf46+cTmbBgFuqj95MJPzln97oBo+FSCMAvIRqGAALmiUWhEAEY+d5QH+STzS7aVp5XYxJu2TW0CQL4EBguIkNEpAYYC4zzHuDK3yQxBkcSPxCuVtprwMnuprOA/7q2dgH+huP4l/tbyETbCftspoRvOW116dWZ3lv1YNqHs4hETYaMGsSMT+a4e91zRDxE4EkaFlhBlFwyGq2NcOaV3+akiw9n9dL1mFGDrj078cJDb/PAjf9zOo75rQ5KXjErFGRHDIhkvh41JLV6yKoxOGHMaIgaPPXsJ4FhIIB4Pdg13tfqXscVcbSj0KNDO2773rfZZVA/7n5rIre//q4znc8KA9l1NppM/ma/3TUemWERLthhD/4142M2JtJEmi0PbZr+K0yKXQoqKHU1AQniopeBNmc0X0x+G+qMRCAx5esH3MEMHvwVtwTU6S3s1dsXSTZfye+w/cI/RpZsNTgJ32K6eIET/km2jYAyJXzLFPYBEJFjgNtwlno+oKo3iMh1wEeqOk5EbsRx+glgNXCRqn7lnvsWsC3QAWfEf66qvigiQ3GWenYDPgW+64aNXgFGAUmttPmqOibf/Ve+8zcOA5GSR/KqmuwWQVIMfkvYOuXyb3HujWfwnzvG8+A1j9GU6gXs/ZnnIh8JeKuIt9t9KJf/5SwGDOvN4jkr6D2wG+vXbOSe3z7D6//52DnISwIBBABg1UVRt3Ulbg4A0vUB1/z2RA7cfxvefH8GN/1lPGs9ye1sErBqwIqAmoJoOg/grQmImia/OvFgRvbtxRn3PUajWFnHuARQY6OmOj8lZ41gmgjcZZ+je/bj8aO/w9TVy/n+64+zsindrDybAETUN+9ffC2AUl9GAoDsb4WUWANQXgIwUAzXQVtq0jlsYZfvcYrhjvSTvcVKXf+fjP2XxfkPHKgDLv1JqGNn/+ynbbrIq+LDPlLCKD3jfBHHhmluUVvP3TuBZ//2Mif++BjunXwzux2xk7snP/nmk6D2Sk1P/XA2Pzzgeu77zVN06dHB4ShbOeHcA9P36ZWVDsgFAJhNcYyE7TyyJKQBrrv2ae742wT23nUo99/yPUa5/QOcYzLDQGYMInGINqnT1tfTOQyc53HL4tonXuGsvz5Bc8Kig0QZ1LVzThjIiBkYloERM5yMsjd+70rsfrJiMedNeIphnbvx9JHfY2inbqlDsqWhVcshC+3YKlYSOjgEVEs69Zt8hO8Gln1f+VA4P5Bpy8YgoRESGkGRUBLQ+e2ZqUWfUEj+OZ+ldL+gKsKj4p3/1wWbGhq548d/57IDr6FpYzO/f+5qvvPLE0Odm2pNGYAkAcRjCZ6682XO3fMaXn7kXbr07EifrXpyzOl70bVXJ065+DAiauWSgIswfQSSBPDMYx9w6aX/wraVP98wlh7dOoROBiefJ/sAJJ9vijnLSX95zME8ee4Z7D9sq8A8QGo5XjyXAF5fNJvTxj9CXSTKk0d8j117plfD5baIpEwEEJwDKJ4A/B19KQRQnIPPRDX5+/VGxYd9WiLvUKmI1kQYe+XxvPfsR8z8dC49B3THtm1WLV5T8NywuQCAoTsM4ILfncqO+45g5ZK19OjbhfnTl/Lnqx5lygezi8oFpLennaddY9C+fS277DmUN96ehl1rUBM1afQk7sIWhnnDQP26duIvZ41heK/u/N+rb3P/ex+Hk4jOkoce2KEzfzvkBH713st8uiZz+fTmloMopxpoUAjo6yIBUWr4pyzyDtWwTxWtiXgswb+uf5KZn84FYOxVx/PQ7Du57j8/Z9/jdycSDZ7EFupE5g0Fzf5iIVcefwvXn/3XlKpp74HduPnpy7j0T6fToV0k9CwgvT1zFrBxYzNvv+osUth95EAevvs8Ro/oHxgG8usaBplhoMVr1vPdOx/lpS9n8vPDDuDmMUdRj5l3OSiQ0yVswYZ1HDvuQT5ZsQhNGOzYPb24ohwzgCD4zQDydQQLQrEzgODwT/4ZQKHZQaEZQMvCPz7XKzH8Uw44oUkJ9WjrqDr/CsBTtz7LY3/8L1vvPJhrn/gpD8+7mzN/fXLg8YXCQJBJAu8+N4kL9vstD1z3NFbCwrJsjjhtTy641g07hcgFeMNAzr40ARgxGzNms37FRpqa49x2/Wl879S9SdZLlUIAjfEEVzz0HLeNf5u9hg6iW7Q+65gQBBBPrynZr+9g/nv02Vy1y8GpSHhLCSAo/BOEYsM/paC1CKAQwhBANfxTWag6/wrA4lnLePCax/ju0Iv55bf+wOTXp9C+S/rHdNh396d959wfV5h+xEkSiDcneOIvL3HeXtfy8iPvArDHISM55rv7MmBYL/r371zSLCC1PWYze/ZyfnjBg7z6xlTO+85+/N9vTqFXh3YZNQF+eQDwJwAjDve+/iHH3vwgy9ZvwEzAiF7dwxMApGYB7y6dxz+nfsIF2+/FrfuOocZwHH9rEUAlxv9bgnLF/0NfbwuO/r8pqMb8KxzDRw/hzok30twY453/fMD4v7/O5Nen5MhGhFUWTeYEho0ayAW/O4VR+4xgU0MT0doIj9z+Ik/ePYG4d0xQpGx0cjnokUeN4sc/PpK773uNZ16eDOS2jiwmDwBw5kGjufyw/bhu/Ks88ekXJTWIuXCHPblqt4N4d+lcLnzjaRriDoG2NAdQro5g34T4f2vJPx8wZFaLY/B1AwbqwIsvD3XszF9cXo35V9F6mPHJHH64x1WM//tr7H7UzvzxpV/xj+l/ZtB2mXpOYWYBkJ4JzPp8AT8/7hZ+d87f2LBuI9GaCN/72bH87bVfMGr0oPQswBMGgsKzgGQY6MXxn/P9c+7lv89+itFss9WAbi0KAwH85/0pvD9nPr/71uFctN8eJc0A/vrFRC5983/s1nMgB/Qdktr9dZsBfF3CP9UZQOuh4kf+e8qhJVXkJpFTmVuptkIUjUVro+x73G4ccPLe3HjmX4g3xzno1L1BhHeenkisKZayJXV1oe7NqK8nWhvh+AsO5YwrjqW2LoqI8KeL7mfCI++64nD1jq1iZwFRQeI23TrX88CjF/P5Fwu47rbnWG85XruoGYA6S/ujKvx27BGM2WMk97/zEX989a2SZgAD2ndi4Yb1ybcLdeWcS5kBZCpzKxHDxjQyf1fF9wR2Vr9HDYuIx1ZxMwCl1kikNPgj4lTkJr8Wpc4ABJsoFjViBaqu5J8BKDUSx3JX55sFbOUf/WuqKmK/IXNaPvLvP1C3ujDcyH/6NdWRf+vDtgPUMfNDbTtXkqFUW37yDrbtXKMctpKVw3kQb47z+uPvcd2pt6RW7xx5zsH84qEf89jie/jRnecxaNv+oIo2NobOCTSvbeCJP7/IObv+klceeRfbtvnBdadwxHf2pX2twNp1zj0XOwtoSiAJm7WrNnL/XRPYffeh/O2P32Vwzy6BeQDwnwGI7fxvW8o1D7/II29O4sw9d2ZE5xJyAIrr+IW9+gzi+THn0LumkyPrXOQMIP2ROQVZqgZxy8SyMz1ZsTOAZKFXzDaJW+nrh58BKKbYJNRI2YqrQdxO2yp1BqAYxIjQqKUkqh1nHXcLxRQhgUmjRgO//vlG/6lyuLa/+Gazo22M/F1UkqRzBoqcAZTTFqrseMB2HHXOIex/8l5EaiLcf/W/eeLmcRm2wuYExDTZeqdBXPiHM9hhnxE0bWzmg5c+4/YrH2WTt6I2zCwgXaAKwKidB/HrG08hEjW57pZnef/jOaHyAAr4+ZmhfbozY/UqrPqQfYKjrjy0RxZiVPc+PHzUaSzZ2MCp4x9ineXEo8POAKKRROaLJPnSlVqf/EDxMwAApT5DHrrw6N8QK8ApKvUeeeiWxP9BaSfB8tD+o3/X5fvcV60k8moF5c4AMjud7Tu4OvIvBm1j5F9FINS2mfz6FG466y+cMehC/vOX55n24Ux3Z/pHGWZ5KCKoKjMnz+dnx9zE365+lEhNhP2P34173vw1g4d2LzoX4MXnk+ZzyTn3sWTxWg7dZ1uAcPUAAfJHs5euQuLw7R235Y6TvkXULFALEJfMWQDw+aql/GDC02zVsQsPHHoKtabDHmFnAKr+RK34T+SKnwGk7SURdvQfNIbwKmu2NP6f3d3Li+LW/wtWwHtZReug6vy/Rli/qoG/Xv4PPnvDaRP6w1vP5sL/+15G0VgoEsCZOT1z10tctM81rFi4mm69O3P7c1ew7zE7hVoSKvEEfh572dJ1XHbRg/zpthcwmm06d6onkghXEJZzjwAGdIrUcvjI4dx58reoKUQAPnh/6QJ+/Ob/GN2rP3/a+9ii6wCCELPK0wsAcucWYQggaFKf3SaydAKQgpLO4QlAA5VLk9gcyd+k/lSYR1tH23H+RnG3Km4CNWBnWa9fdNJ3M9lShRN/fAy3vP5beg3qkbnPjwBcb+ENBc6ftpiL9r2GSW99RU1tlB/+7mSiNfkrgyWWcFtV4uuBmmwb3ZSgvr6Gu248g19cegy1mhRvyyUAURBLAzXwHnl/Mr95+hX233oIfzvtOOoikTwEYDgCulm2Xpw3nT989Doxy8KU9HtaiABs268/QHpDMYVgzVZQKFJ8v7L5CMBPg18VV5kz95xSV/k0amnNaLzvWfLvSBF6/kG2KgkicpSITBORmSJyVZ7jThIRFZHdPNuuds+bJiJHerb/RESmiMgXIvKIiNS52//tHvuFiDwgIgVHGG3D+Zcqw+xHAPlIIa+pgPOKJKVWseV3nmHw1yv+yXWn3sLAbfpx14d/YK9jMztd+RKAmydR1RQJbFi9gV+d8meeuONFuvXqzO8fuZj+Q3rSrsZz3awwkBFLIM1JEkgSgZu4javTN3hdEy+N/4wjD96em359Eh0jUffcXAIwmz0E4P7Yxf3HiMMTH3zOr596ib2GDOKesccTNYxAApAmDwF47P31i4n89K3nSaiN4fl88hGArZIiAOftCueJ/Eb/thrELP/vQPGdwARbjdR9JT+CmF28XEV5l3+mm1V61zjUSzzUzzJz9C/Y7jym0ghAREzgTuBoYCRwuoiM9DmuI3ApMNGzbSRO85ftgaOAu0TEFJH+wI9xmr7vgJO9Guue9m8c/f9RQD1wXqF7rHznbxhICU4R0hLMeB4tkXXOsWWalWErSQA+tt5+5gN+uMfVLJu7nKsf+jGde3TMONc3DJRcJeVZLZXYuIkHfvsMfzj/PrYeNZA7X76SO178OX37dQqcBRixBEbMQpoSSFM87bQ9NQEPP/weN/3hWXbZYSB3/P50enVMVwSDRxo6DpFmiGxSJAHiUwvwn4+/5KrHXuCDOQuxmu2s/en3VhBIGNAkYOE8PKP5AXVdee6Y77NbzwFh3n1sNUhYzsP5PDK5vZgaAEtNGouUjfDDRqsGRbDUSD1s9TZAyEWp4Z/iq38lpeGv7t8bi7CRTQBOZ4CKyxfsAcxU1dmqGsNpwHKcz3HXAzcB3mz2ccCjqtqsqnOAma49cLpZ17utHNsBiwFU9Xl1AXyA0+IxLyre+ZejHaO4M4dvqq0ls5dx2QHXctXRN7BuZQPgNKD3Iuyy0Ncefosrvv0nNjU00WdQd+566UpGH7Bt/jAQjssJKgp76aXP+cUvnqBfn85cdq6zuit7KSi4ISAgEnPGjn7FYM9NnsZfX3ofgG26dqdj1iqnzPi/uMVfHj0gYH2smahh8tcDT2Jghy5AmPi/Yy8e4LiLKwLz/z6UVgAmnoeDUgvA8qE08bfM+2qZ9s8Wcf49ROQjz+N8z77+gFdKdqG7LQURGQ0MVNXnsuz6nquqi4Cbgfk4HbvWqepLWTajwJnA+EI3X/HOv4ryIN4cZ+r7TqP5I846kAem3MLuR+6ccUzYKuEZE6dz0QHXM+3TudS1q+V3D13EiecfXLI2kBGz+fijOVx26UPc+peXMNxRu3clUJhaAO/z2ojJvd87kQfPPIlu0Tp3X7g6gHUxpy+wgfDAwafQMeoQyOasAt6SFcBBaGn4p5xoteSvghEL9wBWqupunsc9YS8jIgZwC/DTIs7pijMrGAL0A9qLyHezDrsLeFNV3ypkr+r8v4GY8u50Vi5eww3PXsX3fzcWw0x/DcKuBlqzYDk/+9bNjP/XW4ghnPXzb9FrQNe8CqFJBBHArFnLaVixAdM0uOHnx3Hg3iN8l4KGIYD4Jotr/vsyw3t05x9nnkyPGlcVNCQBzG1Yw4WvPcOgDl2564ATiLhJYD8C8EM5lEA3hwqoH1or/LOllT83MxYBAz3PB7jbkugI7AC8LiJzgb2AcW7SN+jcw4A5qrpCVePA08A+yYNE5FqgJxCqUKHq/L+BWDRjCZfu+yuev28CY688nj+9/Gu69+uacUwYAoit38Dtl/+bP1/xbwzT4IaHL2arbfrQtXOt7ywgDAEAtDdMunRpx29/PoaTvjW6JAIAeHvKXH748H/ZqmsX/vm9U+hd1949NhwBTFy2gF+8N572kVraR4NXtQQt/yymF0Cp6/+92Byj/0LYnNr/Fa778yEwXESGiEgNTmJ2XHKnqq5T1R6qOlhVBwPvA2NU9SP3uLEiUisiQ4DhOHH8+cBeItJOnPjuocBUABE5DzgSOF1VQy2bqjr/byhiTXFuu+hebjrrDrbeZQjDRw/NOSZsHuCFf7zFVSfeSoeOdfzl+Z/xt9d+yXa7DimZADZubObnP32Ed96bwaU/OJQLzzrQtxYgHwEk/35/6nwueOg/9OvUsSQxuCdnfsHJzz/E2k3p19LS9f/FisD5obXDP60p/lYIX4fRv6omgEuAF3Ec9OOqOkVErhORMQXOnQI8DnyJE7u/WFUtVZ0IPAl8AnyO47+Toaa/Ar2B90RkkohcU+geK17e4Zsu6bw50LlHx1QieOTeI/hq4gxsO/29CCsN0Xvrflz/6I8YtE1fbFv5y1WP8eIj72W2i/RpFRkkDW0Ywg9/fDgnjNmVJ//3MX++71VfOYh8gnDJv0cM6snslatpEiu1zzk2nBhcx2gNfzloDA9N/4QJi5wK6mwJCD/5ByiPDHQxEtCVIP+8OaWfk7IP5ZB0ru8zULf+Tjh5hy9uqco7VNHGkXT8/Yb15uYJ1/CH8b+ka+/Oqf1hE8HLZi7m0iNu5O1xn2CaBpf96XQuufE0TNvj5HwSwUEzANtW7rjtJe6651XGj/8M8JeDCDMDmD5/BbGERRezljtO+TaDOnR2jw03A4jbNl1r23H7fsexXddevq+/2PBPOeL/fmgL4Z8qtjwq3vm3dGaSLFbSEtU8fe25tspyb+W05XmUgsWzlnHrhfey7Z7Dufujm9jxwJEpW2ETwY2r1/P78+7l/uueRlU5+jt7M2qPIejGTWkF1CIIAOCpRycyc9ZyjCaLU44ZTd+6erA0FAF4/zZi0LtTB3Yf1J9/n3UKQzt0QeJgNBcmgKYmm/MmPMW6WBP3H3QKveo7ZIR/VEFtiMfLE/9PWIKVVT3sRwCFwj+qYNmCrY6tloR/VJ2avWTdXkvCP2utehJqkHAL0fwQNvzToLVVXaASUPHOH9tGLatlMsy27X5zW2ZLLSvDVkvkoctlCzzS1Z5HqbZeevA1LtnjKhpWb+Cm8b/kuB8ekWErbB7gyb+8xK9P+wuNDc1cedf32XG3raiLNaKNbjigSAIwmy36dGzP97+7H3fe/F0GmnVENiRSy0IBXwLIXgI6a+EqznrgSaKGyUNnn8rwTt2IbILoenEqkPMQwIr1jZz7ypN0qqnj3oNOpt6MYsUMbEtQW1B1/o81m74OzY8Agkb/cStCLGHSHI8UrF4NIoCELcRtE0sNErZB3Pa/r0LYYNVh2ULCUzCWUAO7gK2g0b+gRLBp1gjNGmGT1mRITWfYKEgAjq04xan0VtEWnH8SpXxrgySdS9DgD7z+FralQX0ASiCApK15Xy7kR3tdzfvPfkzHrh1ybIUlgA/Hf8qPD76O9asa+P1/fsrfJ91E397t0YTrkUMSgMRsxFJWrGjgp5c+RIeOdfzhT2PpYtb4FoJBfgKYsXglZ933BAD/OO9kBnTthFgQ2VR4BjB1+Sp+/MY4ute2p2/7jp7aosxiqnjIJaGQL/zjSBfEEmlbYZd/CoplJyt6Pfdlm2wICAsFjf4FdStoMx+WGjQkih39KybJhjJpWzEiBcnEz1YkVdtbnpG/qBsuDPFo6/jaOv+yJ7Lz2Cv6WuW0VQr5BF889WfjhiZ+e9LN/Pt3TwGw7e7DSsoDLJq1jMuPvJGFM5bSpWdHbn/11/Tv65GYCEkA4Py8Z85YxrW/fJL+A7px/e9Ppl0qpBOOAJKyELNXrOas+55g/qp1RAwjVTEcJgcwYc4cDnnmXmatXpsWGMpAsCMqPv6fDtkkESb+76+Z72rqUFz4x8DKK11dzNr/bEVRLxLqT5hBo/9qoKdlaDvOv4rNjmT+wIyY/PKRy7j9resZuG2/9P58eQCP59mwZiOXH3kjc75YQMeu7bn5f1fQq2f79LGFCCCRuapm0qfz+OON/2Pk9v0ZtcOAnEpgyE8AScxduYbv3vMY81atdV8Q7rGFCaC5yRlx/mr0YXxn60zBvCSKjf8Xw/vFxP+DEEQAfghWICnsgr0EkO/oyl13+PXE19b5l0Mvp9WuVc572wy2rITF9afdSk19Dbe9eR2j9ts2Y38+eegkNq1v5KoxN7Ni4Wq69OjEH5/8MbXegV4+AvCZ3Lw24UvO+s5f+fjTeUB6FZD3b79+AIaPnHPHulp+e/xhDO/Xo6g6AEOEwR268ptdj+DQ/sO9L945V4pbARSctAz/GScJIJ/SZSFr2aN/JTgpmxzJh03+5hNgy9fFy2/0XyWLlqHFzt+VGv1URJ51nw8RkYmuFvVjbnUbbrXaY+72iSIyuLg7LeFWg87Z0raCHHYpiqN5ehYUayvf8TM+mcOl+/2aNUvXceP4X3LQaftk7A8iAO/qo3WrGrhyzM1sbGiktr6G2rpoQT0goznuuIuYnaEBrMDSpevQiLDPToMZ0LeLrxqolwCSwnBmExlSzoYIh44cxvXHH5Yh5Zy+F38CsGLCJW+M44s1S/nzPsezY7e+KScpnqbrYQvALDuPk/X5aPKFf7yyyen/nb6+Xlthwj92Ujo5dW/Om2dkhZbChX8Em0w5Z1VHy79QY5hcCFZS0rnIM6soz8j/UtwSYxc3Abeq6tbAGuBcd/u5wBp3+63ucYXhyhWXMpJPSR0nz61UW9ByW14SaoEMtp+t5PNl81Zw2QHX8NXEGex3/B45p+YQQHIlk7ukFdtmyYLVXHP6HbTrUMc1953LNqO3yu0LkHxtHgIwbEWaLUg4aw3FVjCgQyTKVZcfy2+vGEONWzQWRADJvw0bIptwpJxtWLexiRuff4Od+vfljN128pWBDkJjPMG5Lz/J8sYN3HfgqfSp64gYufH2cAJwQswyiVu51y2++tdxjDaS6qsTMSxMoxQ36TRZT2BiqSAuiZRmC2wMLAxsnGWjFgaNGi04ifUf/RskXFtlgSZzP4UfbR0tcv4iMgA4FrjPfS7AITglyAD/AI53/z7OfY67/1AJ4e1aor/v3pNjwzQr11YLtPxTtkRazZbX3oa1G7n66N/zx3PuBKBLz045wnA58I7Ym5qY+uFs/nTx39l+z63505OX8psHz6fW60iCZgAKRsLGaLbcdntKY2OMm256luFb9+bH5x4S2BMYMuP/ohBtAjPm/P38x1/x1oy5XH7wvvTp1CF8+CdusLKpkXNefhIUtu/aBzsefqWPHwHYAYnPYtQ/k7YUwcbAxqApYJQfLvnr2EkQoUlrfbuBQfjRv3NfJhamO0tpSfjSsVVFcWjpyP824Oeko7LdgbWurgVkalinNKrd/evc4zMgIucn9bHjhFtRUsXmRTyWINYUJ1oT4aaXfsVvn76CuvZpCYhCBWF2YyNvj/uE+37zFNGaCDvsOYxf3vN9IupJ7BZRB/DBxFn8+7H3GHPUThx+4MiiqoC9f1//9AREhKsOOcDZHjAD8COA2etXs88Td6ekH1pL/8cPpRR/VSK+YaqfWxwlO38R+RawXFU/LuP9oKr3JPWxo4TTlKliyyAeSzDurhfZ7cid+b9Xr6Vbny4Z+wsRwFN3vsxzf38DEWH3Q7bn53/5HoYVXgrCiwfvfYPJny/gih8eTq8eHXMIAArXACxas54rnnieP774pu81kqP/IMRsC00YHNDXEclrDf3/coi/+WFzCb9VZR8qBy0Z+e8LjHG1qB/FCffcDnRxW4xBpoZ1SqPa3d8ZWNWC61dRAXju3glcc/wfGTCiH7e/fT2DtstoVlSwHuCuqx/jg5c/x7Zs9v/WLhz9nX2cJHAyEZyHALJ1gG647j/c8dcJrFy0DihyBZC7743PZ7N0/QaMONSYZlHhH4CjthrBPw49jWMGOSuiWkoAYVEO7Z9KwJYe/Ys6340wj7aOkp2/ql6tqgNcLeqxwKuq+h3gNeBk97CzgP+6f49zn+Puf1UrWVK0itD4cPwkfnrIb4jURPjRX87N2R9EAHZjI7Zlc+MP7mP2lIU0N8aYNnFG+oAiCWDVqg2MHzcJgB7t0w6iUAIYMmsAInHhnjNP4NdHHJyxLzAB7CGAl+bPYOrq5Vy5y8HUGI7jD9sAxg/V0X/bhYgcJSLT3BWOV/nsv1BEPnclmN/2NngXkavd86aJyJGe7V1E5EkR+UpEporI3u72U0RkiojYbkOYgmiNdf5XApeLyEycmP797vb7ge7u9suBnDejiraLmZ/O5dL9fs1NZ90B5C4bzUcATRubufaMO1m/egPXPnQxw7bpw5hznLh7sQQAMGRwDx667wcce9iowASwEQ9WAbVVmbF4JafuNoo9+2fOZJzrZY3+va8nJtzw4WsM6tiFs4bvntqeTQAtHf0Xn/zNRLHKn0FoDW3/LT36LwdExATuBI4GRgKne527i4dVdZSq7gz8EaetI+5xY4HtgaOAu1x74ERXxqvqtsBOpFdafgGcCPjHLH1QFuevqq+r6rfcv2er6h6qurWqnqKqze72Jvf51u7+2eW4dhWVg+XzV7Jy0WoMQ/j1Yz9h7M+PC3We3djI6mXruOb0O6hvX8v1j1zCRdefzCkXH+Yc4EMA+bBw5gq+/GoJl51/KMMG9wy9AigJIw53vfIeC1ev4zffPqy48A/w9vz5vLpgFpfstA9dzXQlc6kEUAnJ33L3/N2cHb+2EPYAZrr+MIYTGs/4Qajqes/T9qTLFY4DHlXVZlWdA8wE9hCRzsABuANqVY2p6lr376mqOq2YG/zaVvhWseVgmAaxphjfv+F0Lrv7B6mloIXi/3OnLuZ35/yNjl3bs2rpWr5/9RiO/d5+zs4sAigU/7/xd/+lYUMT1/1sDO3qa0IlgCEd4mnelOC6ZyYwtGc3Ltpnj4x9BcM/wA3vv87qpk0M6NAZ9Tj5MATgh0oI/wShkhq7b2b0SK5MdB/ne/alVje68K58TEFELhaRWTgj/x8XOHcIsAL4u1tYe5+ItKdEVLzzL1WCOXW+K5+cenwTbLVEHjrZY6AFthJxiz987w7+/funOea8Q7numZ9RVx9FLQu7yb/Dk93odG2a9OZX/OWKf9O9TxeWLVjFJb8/lUNP2t2pE2hqch5r18OmRqTJIQQ/Ali7dhO/u/6/9OvbhSt/6IRMMzSAmm3MRpuadTaRjTaRRjsn/v/OjHn879OpfGvHbam1zYx9zvXc0X9ckIQgcUEaDbBg1trVHPL4/Xy2aqnzvhZBALF4hFjcJJ4wsO38K4zyhX+SWv4J9xG3ciuIw4/+lYZ4DY1WhEYrQjxPNbIX/sSgCMpGu5ZGO0osQGq60Oi/wap3ewwI8YDaiKJRXJHXyuTKRPdxTwHruZdTvVNVh+GEy39V4PAIMBq4W1V3ATbSgvB5xTt/oEW6+Tmql+W2VaoMc2va8ttWjC3v+5Psg1DCe/bgrx7hlvP/yq6H78hv/3OlszGPrSQBvPTwuzz8f8/Re2B3Fs9Zzok/OBjTILM3Q3MM1m/ISwCff7aABx54g/XrG0l2UkwSgCScSl9xf+zRBkUSmhP///241zjptoeIW1ZWaCgz3CMKgjg9AZodArBVqbWjfHvIds5bmYcAciGoGiQsA8ut+C0m/KPqaAW5rxYQGq0IzVY4J5k5+ncqeg1R120LcTWI2Wlb4Uf/ipG6o2SFrkmzRktSbbfcauEK1PhMrW504V356IdHSRfEBp27EFjo9vIFp1jWX1UwBFq2vmxzwrbBLJLd8+n5l8uWq11TVFVtOW0VkIfenAJ32dcGeOG+CaxavIZVi1en9zU1IfX+CTu7sRGjvp5//eF/9BnUg0NO2ZM/X/4QViyRK3CjCk3NUBc8en30kfexkvIRkbR6jNeS44SgpsGmuavpJH3dAfWG9c3YUag1TYYP7MFnK5ZhJJwewEZcsOo1w5ok/44ZaL3NqdvsyPX7Hsba5kbeWjwXTRiI2+vXipmpPsCJhBHQA1iwbAPDsBBxCCC7929jLJrT9zdmmZiGnfNKY7ZJrVoZb+WmRE1g31/nLJdCJHOrhSPPEFTtm8QGqy6n369kPbNxZB/MLKGGtVa7PP1+FVuNsmoblhEfAsNFZAiO4x4LnOE9QESGq2pyeduxQPLvccDDInIL0A8YDnygqpaILBCRbdz4/qE4Td5LQtsY+ZeAb8wq0nK+znyzhRb0LPjg+U+YNWkuAF17d3FILkQ/gNsu+xefvTOdi248jX2O3okr/nIWow/aLvOgRKLg6h8zZrP10F7c+ptT6N61ve9rEZzRflAB2C/GHMTfzz6JIR27ONsKNfNw1cYe++oz5q1fy693P5SIuLmPgBlAOT9KkUxS8sJPPdQv/JMc/Yv49QZwYHtshRn9B/tpKboVY/FCcJsProrBJcCLOCtyHlfVKSJynYiMcQ+7xF2eOQlnBeRZ7rlTgMdxHPt44GLVVPn7j4B/i8hnwM7A7wFE5AQRWQjsDTwnIi8WusevrfOvogS0sjz06VefwP1f3kqnHp2A/Ms/wakgvv6svzJ7yiKuvvcHbL/nMH7zjwvZ8/BR6YPd338hAohYysht+/HHX51Iu7r8yU+/BPA9r31AwrK5bey3aBd2wpwQYrbF9e+9yoiuPThnZP7l12Gdf9jkb5C9Ziv8hH+jVeO0pwywlf0xBxFAYWJQl7ByERT7bw3XL+qEBsM8CkFVn1fVEao6TFVvcLddo6rj3L8vVdXtVXVnVT3YdfrJc29wz9tGVV/wbJ/k5hd2VNXjVXWNu/0Zt+6qVlV7q+qR2feTjbbj/ItUqcwb7ijFyeU5p+jQSp7XUk5bRaPc71kW3vnPh7TrWM/Z152W2laIADas28QvTr6NaZ/OpWe/rqxYspar7zmXbsmuYtFoqBVAM2Ys47rfPsOwIb343U+Pw8yKVShg1YhvARjAshUNXPnYeEb06sGvjk0Xf0kcHy+kaFKoLi68Mn8WL8+byWU770vfdk4XM7/Rf5CTNfKMvL3wEoCtft8Lv/BNGkHJX/UZr6u6s6UiXbBm/e9FpEhtTk3KOVfuBKCi0Tacfwna9IC/YxQp3fn7nVeC85Vy2/I7rwTV0da2NX/qQsbd9SLH/OAwhu44KLW9EAFsamjiV6f+mS8mzqTPoO5EayKceNGhzntYRO7mgw9mc/Pt49lz9BCuvujIDK+hBniFRf1G/+9OmcvfXpvISaN34OQdtneOawKSDWIUUm4t62387bsT+GzFUtpF82jwq/hG3vwceaHkr7phFE3dk1vj4I6ui+v65fTrTYqzJt+2WiNR1E/JGf1LysWr52EhJfwshZhGXHKqMkCxqHzn3wJt+qQ0cUqT3rVVqm6+eOyQlGMucURcVlteTf8Ks5X9/v/r+qfYsGYDF91ydsaxhQigaWMz155+B5PenAYCBx6/G2ZttKj1/wAvPTuZv//rbfp270Q7M4qRXP7pcbr55B/ueuV9Hn1vMpMWLHHsI5hxwWjEIYG4gIDh1eSPCws3rOf05x5j1sq1ed8z9R2xh4d39K8YJNQgbpuYokR8+gyEgRP7dwgg/TDZENC8vXCIx+0zAKlHodU6wcs+HQJososhsyqgDTj/cqxWSWnUV21tdlvZ9jas3ciD1z7OkFGD6DOkV8ZxhQiguTHO/13yIImYxaevf4mV8A8TFCKAh/7+Flf84jHi65udzl4+AnDgX/1LTLn+v68ye+VqZ4moW/0r7jof7zvmJ//Qra6eq3Y9iFrTDEz8lrfy17krv4/Sb/Sff92/eB6lIZMYWm4v01YVxaDinX8VXz+8cN8EztnuMpbOWZ6zrxABrF62jvEPvc3BJ+1J70HdidZGStL/0U0JOnSo5fbrTmXwwO7uvnDVv86x8LvjD+e2E47xnJMp/ZABd9+23Xpy4ag9uWCHPX1fZznQUt0fP5RT9C0IFSH5YDvfgzCPto6q869is8O2lYbVGzAMYehOW+XsL0QAT/z5RdS2ueP5n3HKRYdlHuSj/xNEANEEbDuiL2eeslde+eegv2csWslh223NuXvvWlj7x8W7i+fzv1lTuXjHvRnUsUvG6L8QWqr744fiR/9VfF1Qdf5VbDGcff1YbnvzOnoOyGnolpcAVi5Zy4v/fpd2neo44QcH0a5DXWYjeBfJ0X8Q1q9v5D/Pfsoh+23LgH5dgfzib5C59v9f73zK+C+mc/nB+7HboP6Ba/+zxd+uf/81YpbFb/fMJK5CoZ9isDlH/0EoZdnnN0XuuRJQdf5VbDE8+7eXERHOu/EM3/35CODxP4/Htmw6dG6XFn9Loojwz1OPvk88bvHdk/fMK/7m/dvr5K954mUWrFnLrSceS4/27dzj8od/lm/ayK0fv8PBA4Zx5KARvq8xCK0x+vdDayt+loqvgdpnxaDq/KvYYlg+fyWP3zyOg8fuyw77but7TBABrFi0hpcefhfbtjn5wkOorYtmjv5Dhn/WrNnEsy9M5siDtqdvr85FJX+NOGxsjnHpo89i2zZbdesSOvzzjymf8NfPJzJp5eLA0M/mGv23puJnECp19C8arsArTJFXpaPq/KvYonj8T+NYsWAVF916FkaASIwfAdiNjTx2+4vYttKpWwcOPtGtni0Q/skmAIDHH36Phx59j41r0hoyxSR/Zy1cxRG3/Z1PZy/O2eeLuGCp8of332TZpg0FDs5FWxr9b6lmL1UURsU7/3Jo9KgrmFbJtlpqr63aam6Mce9VD1Hfvo7u/bsF22luduWmNSUzvXzBKl559D3isQQfvuajb+UT/smAbWNsSrBqRQP/vP9NNmxwSKaU5K+1yaIGk29vv21q0WGh0T9A//ad+Ofhp7J1B2fZqxUzM4qpyjH6t2zB9hRnlWP03/Kvv2KptxAtE8WM/tMlbFUUg4p3/iS15UuUYVbLSksBV4gtwJFvznpUhK2krHNr2fIpYX398fc4f+crWLFgVX6DsTgkEpCwIJ7AatjAY7eNxzCEk887yHOcf/jHG/+XeAJjUxyjKU5kbRPGxhh77zGMU49OK+QWk/wVCw7Zeig3n3g0hw8els4PFOj8tanRYlT3Pvxur8Nx9KAVtSX1sC1/51jM6D9umcQSEWKJCHaRH2XQ6H+THc0RYgs/+ndqcpvsGprtKI12DVaJxW0KbsFYdZ1/sah8559EKQ4oj3Ry2a5fqm6+n70y2yrWafv2Bii3rYAZQCJuUde+lt2P2jnAos+P27ZZMnMpE56YyLHnHMgd43+OGQn3lRZbM0qMjGaLA3Yfxg/OOYBuXdqHTv4mtQoEeO2LmcxfuZYfHLg7ZlOI8A+wprmRmz5+gz37DOKEoTu4LzOz+Cneggbw4rGnQDwRQbXlo38QmuxICT8lzXp1zl/NRdpaa7VLOf7k/ZQFCkazHerR1vG1df55nVU5nX+ha/khj5Mvp62iUe73LAgB93zmNafw26evYPD2A/zP8ylTtRsbefTWFzBMg2E7DOCQE92m6UHJ32T4J8uWAI/f9yaRiMnp33byB6GSv54+IpatPPj6R+w4qC97Dky/hkKj/8e+/JxPVyzmF7scSqea2uwXDfi//cXH/qWkEEncDnYTFqWO/nORyOrGFT70U5mjfhE5SkSmichMEcnpuCUil4vIlyLymYhMEJGtPPvOEpEZ7uMsn3PHicgXnuenuPLQtojkl4910XacfxVtCyWQxeN/+i8b1mzkhv9dTZ/BPcNeiCVzVjDhsfewbeW8Xx9Plx4dcw8L0fx98fxVvPbql4w5dmc6eQTYCiV/vfjvR1+yesMmTtt7R2cUGbD6J/MVwK/ee5GutfVcsN1evnZLif0HqXOoG65p+eg/bascKPYrs9HKJsrKgYiYwJ3A0cBI4HQRGZl12KfAbqq6I05Xrj+653YDrgX2xGkEf62IdPXYPhHIXinwBXAi8GbYe/zaOv8t1cEqFFpZN79VbBV7nXzHBwj1rVvZwFVH30Bt+1r++NKv6Tkwt/gr9zoGdmMjD1z/HxobGunQuZ6f3OzWDfis/AmCAmoaPPbYROrraxhzzM5A8Og/BZuMoXQsYfHCp9Po3aUjppV7eBLZo/8py1Zy4ZtPcccX7waeU6zmT5AzNQJ084PgLw8NIJiSO4srdd2/aeTayjf6tyrbfe0BzFTV2aoaw2nTeJz3AFV9TVWTS8zex2nXCHAk8LKqrnb1+l8GjgIQkQ44jV9+l2VrqtvdKzQq+t3LQCnKnkHnlNlW0UQTdHwpYmrltJX3MiVIOpeA2Z/N5+pjfk/Hbh244r6LPHuy9IST/5vO57JuZQP3XPskhmGwx2HbM3K3Ic7+rPCPAMStDFtJV6g1EWbPXs6zz37K8hUNGXHdoNG/mTSv6f/+9OybnHHfY6lbzB79B+HlubNotOJExcSQpEyxIiWqcUL226VEjMwWjmFG/5YKVkAj+UItHL1ISjo7d+K5R8DALqE3gHgCWVtkvU8PEfnI8zjfs68/sMDzfKG7LQjnAsmmLfnOvR74PyCot2VoVH4PX1f7vlQZZjWMzBhzCySdy27L68wq2FapENPMTPyG/CxnfDybq46+gbXL12ftUUDce1OIRlO27MZGXvr3uxx26l4M3WEgi+asSJ8Wi0GNu2qlOY7URtF4wnltpoEAdk0EDMGI2dx6y/h031+cEb8VlYy/k31+RR1Nf6sGMBzXFsfGSEC0nUkjucN/iQlaoxgJwY54nZbQw+zI40efwb1ffcAjMz7FiNgZztqv128sHqEmmpld9vb6TROAEDFLcZJOQradkRs622jV0N4MP8OCXDetFN8bIAkDULRsq31EnbafIbFSVUPF1/NeU+S7wG7AgQWO2xkYpqo/EZHBLb1uxY/8S9XfT53vasqnHpVkK6mXX+G2yvaeFfFZTvtwFsvmrUBEOP2q4+nUPRnHT4/y/Gz95YqHqamLcv6vjmPY9gMCV/+IglieH7lIRuFXPQaHH7p9hiPwU3I0Y4qo0w8gtbpH4JCRw3jnygsZ0N5pWZl35U9yRhA3WNG4iaWbNnDlTgfRvaZ9WaN6LUNxN5I/9JNcfSSptT8b84R38oV+1lrtHFnuylvpvwgY6Hk+wN2WARE5DPglMEZVmwucuzewm4jMBd4GRojI66XeYMU7/yq+2Ri0XX/O+MWJ/GH8L+jYtX3eY+3GRhbMWMrjt4/nkFP25Pbnr+D0S91WpgXW/nthxGz23W8Ev/zZt9hr92FA8YVfUxcvp31tDWN2ypStyKv74+LX771M59p6xg7fOUPsLYmW6P23RPKhHGqfrVHxW6H4EBguIkNEpAYYC4zzHiAiuwB/w3H8Xn3zF4EjRKSrm+g9AnhRVe9W1X6qOhjYD5iuqgeVeoNV519FRWPelwv5zUk3M2i7Afz++V/QrlN9al+Q7s9jt49n4cylNG+KMfZHR7Dt6MHOjjzaP9myD2++8RXLlq3jtJP3CIz9J+FX+LVkbQMfzF7AcTuNxCda4jk+lwRmrlvFxysWctSgbQB8CaCSUC69n68TVDUBXILjyKcCj6vqFBG5TkTGuIf9CegAPCEik0Qk2dh9NU5s/0P3cZ27LRAicoKILMSZHTwnIi8Wuseq86+i4vHxy5/xu9NuZdhOW3HDs1dT3yF49Gg3NhJvTvCXKx6mXcc6mjY187M/n0ldO9dB+RCA3+hfGxM89dSH7LzjILYZ0QcofvQ/7pOpDO7RlR0HOOfnTfx6Qj8AL8ydzg7d+jCgfWff19lS2YdS0Zpqn6WKvZVV68dWjJgd6lEIqvq8qo5Q1WGqeoO77RpVTTr5w1S1t6ru7D7GeM59QFW3dh9/97E9V1V38Dx/RlUHqGqta/PIQvdXdf5VtAm8/9wn3HDG7QzZYaBvA5hsfPbOdF565F1q29XSZ1APzr/2hILnGFmzgeefm8yGDU2MPTmz61bY0f9Ln8+gKZ7guJ228x395wv9jJs9lbNffYxljQ1AuNH/5gj9lAPfoNBPRaPq/KtoM3jnPx/yveE/Zso76eXMQYqfAPf/5ik2rmtk1ZI1LMte/ZOET/FXMvTTvK6J//3vU3r26EhyMU0xo//GDTGueeplHvngs/R+n8SvX+hneeMGXp8/l7hndVk2AWyp0b8fqqGftofK+fZUUUUIrF/VgGEaDNou35Jp99jVG7n3mifo2b8b69dshGaP088qAMtu+pLE3x94k0sv+SeWlTnNDzv6f27SV8xe5AjWZZJE4dBPn3YduWLHg+hV3yHwNWYTQKUlfttk6Ocbgqrzr6LN4ZI/n8Mtr/+WaE1wmUpy9D/h8YlMevMrvn/NCRxy6p7p6l/IWwGcHP0nnX7nTvV0NJ3rFRv733VIf767587pbZ7Rf77QT/tolIt33JujBm6T2lbJyd9iRv/V0M+WR8U7/5bIJkNaijn1+LraSkpMu7ZaIsOcYauFks7lspWyZ1m88/REOnXrwB5H7+JsD1j1k8QdP3uYmtooJ1xwKEeM3Yujv7OPU/lkWdDYCA0bwbYDl3127lzPo/+4kBOP3y1rX7p8VizFbFKiG2zEcgnBNXf49lvzsyP3p1NdphaNd/SfCv0kyxhiwqyVa5ixZiVHepw/FA7/+I3+YwkDO6BKNwxUwfY84rbRAq0/RbBptk2abZNEi20ppk8xXYnmkJgV6tHWUfHOHyhZU95XVriSbZUq6WzbmUIuSf38cthyt5XktMtoSz2v6ZMJn7Nq8WoO++7+eW0lR/+LZi/n0VtfYOsdBzHr8/mcf80J9BvUzdsxBdY3gDvKz172uW5dI5MnzefEMaOpc3VukiN+xXXymiwcg+gGxWxM39ezH06lJhLhqB1GFEj8qqf8VUDhhTnT2bP3QLqa+WscCkExiNsGCStNAOFDPxFPCZUrwWxFaLbCzUKyQz9JKQfFQDFIYBBTM/VVKSb0Y2BjYhclM1GFg7bh/KE06eKgc8psq2wyzKWM2MspNV1GlFNSO7tngW3ZTHj4bfY8djSdugXHw7148o6XmD9tCZ26dSQeT/CzO8/OrP5VoLHJ91wjZvP44xPp2qU9Rx62fe5rkXT9a1KnPtKUHv1/uXA5M5et4rgdt0vbDKH1D8L4OTMwDYPDBw7P2FNo9O918l57lhY/yk5r7nhtCgk1sbNUPQuHftya3ozTxJ0LFOvB1VHUrjr+ktB2nH+R2JKOb7NiM+jsl3SdVu4N8Mq/3iASjXDASc4yzKDQT3L0H48l+PMV/6Zn/65M/XgO244ewn7H7pJ5cCIRmPidPHk+X321hLEn75m58idIV8/KTAo/+9FURg/uz8CunYMTvyY59r5ctZwpK5fRq12HnEbv+eL/+ZQ7i5Vhdmz5n5PdzasQJI8Mg5dIwuQEqj6/ZfjaOv8qSkClSk37YM7n8/nRXlcz/oHXCh6bJIAp78/kv/e+yh6H7sDff/8f3vjvxwXPTYZ+jJjNP//xFn37dGHnHQeFu0lJNm+BZz+ZypK1DQzqmC7aykn8BvjFY5/+B3+Z7Eg9ZxOAF97Rv+Lf/tG5WHpHmNCPEszXYT/lMKt+wn5liunvWyxEFaM5HurR1tF2nH+R6pJ5BcRKcUx5zila+CzPaymnraJR7vcsCGX6LL/6YCYJzxLMQolfgPt/+zSzvljAiRccRvc+nWnXsS69aiji/O/b7B2YOHEW55xzDx9/Ojcl+SAWOZ5RgaymVCxbu4HDb7qP92bOD7w3SQQQQDKW5IOg0b8d0H1LKE6GGYL0/J0b9dPzzxf60Tx0YVJcODb5Vn1TJvnlRttw/qUqSwY5mXI6/xKcr7jSxq1tqyQNfr972MK2kudlo+egHlzwpzPpN6x3wdNT4Z/mBDeedx81dRGu/eeFPDHtZnbcd4Tj+COZq2SyE78ACxc4EitduzjryiNN6jZ00XS+Vhy5Z8hc9ikxMERoX1vjH/qxATd57CWBqJi8MOYcLt5x74KvMwlVwfKJ+0fzdZgJsoW4rRwV780J0Fh0cZdgYaRy7ck+AzVilfCzFBKuC6sSQPEo2fmLyEARec3tQTlFRC51t3cTkZfd3pMvJ9uPiYM/u/0sPxOR0eHu0EBKHN0m5YQxjNSjVMnjctoCR6q6tWyV7GDxOO0KtJX9/keiEU7+ybfY9fAdU8flG/2nVv/MWsZdVz3G8B0dmYgd9h6evEjee0jquXz78FE88uCF9OvSwZETToDhlgyIgpoOIWcXfRkiPPvTs/jpIftmbE/bNzAsAyxIrVwUiKuNrcq+fZ37DRv6sbOnH0DMKq3gy8agyYpkNF8v9qNMh34cAkg+aiXhm6MIVwsgJDBTJFBFeLTkHUsAP1XVkcBewMVuj8qrgAmqOhyY4D4Hp5flcPdxPnB3mIuUoxtVUpP+m2KrpfYq1Va2vaVzV7B8/kpG7Z/dGrUwXnn0PV59YiIiws77edbR5xF7S+LDD2YjAueedQDgJH69TjEItipTFi7j6B23IWrm++m5ljyyD+8ums+uvfpTa+Y6781X+CXJfjxliAKm37HyRBSr6d9iUbLzV9UlqvqJ+3cDjmxpf5w+lf9wD/sHcLz793HAP9XB+0AXEelb6vWrqAJg8ptfstOB22VsCzP6B3jyzpcQEbYeNTCwWjg79GPEbFasaOCJJz7g8EO2Z9ut+6SP9Wn1mP33sx9OpUv7evYfPjhw1Y+f1s87i+dRa0YY3bOf8xpD6vr4FXxtDrR2tW+rJX1VkVgi1KMQROQoEZnmRjuu8tl/gIh8IiIJETk5a99ZbvRkhoic5W5rJyLPichXbrTlD57ja0XkMfdaE8N0+irLXMm90C7ARKC3qi5xdy0FkgHZYntaVlFFQXz+5lS69u7CwG36hT4nSQBzv1xM48YmIlGT4TsOLHBWJh595H1Wr9nID886EAhu9J4d+nl3+jxWNmzkuFEjM7YnkSP34JLCh0sXkrDtVOinJQir9RMW5Wjy8nWDiJjAnTgRj5HA6W5kxIv5wNnAw1nndgOuBfbEaQR/bTJ8Dtysqtvi+Nt9ReRod/u5wBpV3Rq4Fbip0D222PmL003+KeAyVc1ouqqqWamrUPbOF7chcpzCqzeq+GZj8htfsnbFenoN6pGxvdDKH7uxEVVlyvszWbNiPUvmr0pr/WSFfvwSv42NMR7819tsM6IPA/p1JSwsW3lh8jQO2m5IjtyDF9mj/43xOHdOep9PVyz2t+sJ/ZRL7XNzyjwHoQ1rAO0BzFTV2aoaAx7FiX6k4GryfwY5y5yOBF5W1dWqugZ4GThKVTep6mvuuTHgE5wWj5AZcXkSOFQKxFlbNCcUkSiO4/+3qj7tbl4mIn1VdYkb1km2JwvV01JV7wHuAegk3ao5/CryYsnsZZzW/wLfoj5tbkZqgx0swKdvfMVuh+6Q95hsGDEbu8bghXGf8vb7M1i9eiPUGu4+xa5x1/fHM1f9JP9+9K3JvDVtLhtjMQzbe4xgR32+8nGBqHLrx+84z93jNWEgkRKq1bcg1sfr6BTNraTeYNXRwfSvsG6j8It07BlwbJhzM6IkItIF+DZwe/Y5qpoQkXVAd2Bl0EVastpHgPuBqap6i2fXOOAs9++zgP96tn/PXfWzF7DOEx6qooqSUWo1t93YyBfvzQDgtIsPKzohbdvqOH6gT69OoUM/81eu5Z0Z87Ds/PftF/sf2LEz/dp39Dk6GGFlnsuN1tb4b5W4v63OzC/MA3okoxTu4/zy31AuRCQCPAL8WVVnl2qnJfPDfYEzgUPc/pOTROQY4A/A4SIyAzjMfQ7wPDAbmAncC/ywBdeuoooURu23LX//8lYGjMhdP1Ao/DPz8wU0N8U47tyD2Grb3PODQj/JZZ9mzOaiHxzMvTefSXu3VaSf1n82enZsz88O3Z9+XTo55xSSeY4LtabJKyd/n7O32T39+jwhnpau+tnSzd3bIFaq6m6exz2efaEiHQEodO49wAxVvc3vHJccOgOr8l2kJat93lZVUdUdPT0on1fVVap6qKoOd3tUrnaPV1W92O1nOUpVPwp5nVJvMcNGS2SOv8m2WmqvnLay7SWxZvk6+g/vy44HFLnkUxUrFmfax3MB2GGPYTlx/zB45bUv6dy5Hd89aU9n5OhB0KqfGjX4/oG7MWanbXNX/Shu4Vjm6L/Zsvh4+SL26VdYXiI77l/uIqiW2vMWeRWyFSbuX4Hx4Q+B4SIyRERqgLE40Y8weBE4QkS6uoneI9xtiMjvcBz7ZVnneCMuJwOvaoEfXOVXRpQodQxp/feUxHFSW75EGeZy2QKPRHGl2UrKOnvtVYitnPff/V4snL6EVUvWBDr/nNG/x+NoczMfvTQJgJ33HU4xMGI2Csz8aikvvvwFp4zZjX517TE3WaCaN/SzZE0DH8xewJid0stUFUXdYieJCxJzqoJTni0uvLt4PiO796KLUVznqngiQjzREt38TDRZkWye84Vf6KchXkNCJeNha2nJXQWskhRBWxeqmgAuwXHaU4HHVXWKiFwnImMARGR3EVkInAL8TUSmuOeuBq7HIZAPgetUdbWIDAB+ibN66BM32nKee8n7ge4iMhO4nHR9VSAq3/lDaVLHkFc6uZR7KOoaeU0FDHfKbKvY98y3z0AF2Eqe53ORlK3P35zKjgdsl3tM8tA84Z/P33Z6Ao/aa2vf/flW/ahbq/TvR96jpibCvntvjZGASEPhz/J/n0xlSI9ujOrvrIYWGzAyyp9AwWzKLPYyRNizbzoq4B/68UmAqxBPpENDLVnyqTgE4P0ow4V+lFojAZm1wliaR4guj620068s5w/gRkJGuNGOG9xt16jqOPfvD1V1gKq2V9Xuqrq959wHVHVr9/F3d9tCN9qynSfacp+7r0lVT3GP3yNMLqBtOH8o2jHmdTDldP6FruWHPK+lnLaKRivLMKdQzs/StTX5zS/p0b8b/TxFVwHGcjZN/3g28ViCTt060Kt/1+JCP4ZT7rpg4WqWLlvHjqMGOpIPIQq+Xvl0Bk3xBGN22g4jDhohx4fJ/7d35kFyVHUc//y690oiuSBQIQdJIAmQcAmSpADRiBCQIpbiAZamFK8qKIHSEvnDsyyrrMKzSi3xwLMABQpTEVg1KJaUiSQCOci1Icnu5thdQg4Xsrsz3T//eG92Z2e658pk53qfqs7OvO7+vX7zJr/u+b33+z5kOAQEsOlgD28khvKGfjxPI5pqjBfblYkgyk0YW5l6/vmIEoJLUew32UOtPHUZHb+q6f9CthqnMul/juqknA5+jJW2Xnx2C3/5zXNojliEDg4iLdlPp8lEwM6Nr9LS2kzv/iMF1SdDAdoyeoD1C/c9Sk/vsbQKFT8BQbNxTv6QEqSmgSahf2CIp17enmY0vj4vaaaBJsOQO9qfoOPoYbPQe7NxmVHTPuMmL6kKkkPvP5Ncxxbby7ncdK51BqKmglbfs35tUbfOXyTXshF1hEj5HK3nxT+VFyvAkuu6TsG6AQc6DvHAHQXJRUWyZ0sX1966NO9x3mCCsDUtNBIq+OYa9h8YuXGoZ+f254mEfPmxv5oXLRgxt4gFXTJZ19kNUfkAaaSigVEfdTGOH1KSztFqoLkWjYkiUKEpTiBXtKj5/oFz/ydF7YR9SlH2jDunzLaKFizLITVdtK1yrjNQRnLWXYo8dNw5GeVzFs2kbXx8YpfG/Fzv3HGA06ZM4IN3XhddTYyWi58Ihz3tWWdO5LOfeiezZk0laBvdx3GhnxRTxo/DS/lXHfmjZIeDJre28YnFlzNvUnxmcdxTtOeFoz6yQuL+cat1+ZK9dm5U3D990FfxIuL7I8J4xZGSdG6Ix7yyUxvOv0Qp4GFJ4VqxVaqev+eNdoIp2eNy2LJlJd1IymhLMtsUYWvRVQt58KUHuPz6iyMsjBB1A+jaYxIhb/ro1aN35Ij7y1CAhOAPBKDQ0uLzoQ9cycJFZ+MV6I+8BNxzw1U89blVtCZ9vEHI1PPHRpdS0z5bm5r4yrLlXDNzjgn9ZGAGfaP1/Jv8UhylcCLiJhFESEYXQkKbSKhvl3Q0y7j7oiX9IEyob9cYLtMNQEMYHCpsq3Gq3vmfjM49jOjAD2/1aiu1JoC1Vaq9LFsnq8FfJlvD9nLY2r6+g+Ov97P05ssLMTZq6+4wKiTTpk+mqdnPOaCXuYSfKHihcqjrCMlkwOxZp4/aHzflM8VLnQeZPH4cy86djSBIKEggsRrRPW/00/tmPxefkWdwu6xhkfLaCtSnyVOaPcXPeEYo1lZ/2Ibvnv6Lpuqdv8NRKEEyYEP7Syy56a14edYqzJz6+dqBIwwODOH5HjPmnRl5Tj4Z3yAIOXDwKLOnFyb0lsrqfX7nPo6fGODGxQti6s1uy6a+Q1w0LZ/zry9qWOStKnHO31FXrFuzkcnTJnL+kvwJW5k3gIN7+wA4Z0HhTjVT7qGz6/CoJ/9CpB4SQcDarbu57oLzaPb9UVIPcWzu6eG8yaczoXkkPl+Ixn+l9P2LwTn5scE5f0dd8UL7yyQTSZa+p7BVQtPZs6UbVWXqmZNG74iI+2eGflJ0dr3OW8a3RoYxcg36PrNpB6e1tfL2c89JOybbSCruv6mvB1Xl3EmnZx2TTiHyzpUSeRtZ1rGKCJXwxImCtlqn+h8DHI4ieOPYm3xpxbfY9d/ixA7DEyfYt+MgIkL7I/82hUNDEJEXkIuHfvYPHnzoObTVw0/o8Bz/XHgJWNfRxX2PP8MLHd0F1fN81z4W/eoHDARJYGS+f7k4MdTMuJb8SW5vJpsZ35TIKGthfFPtD4jWO+7J31F3bPrnK5zoL14bvmvnIQBmnhsd84fouH966CeMSTLLN+ibDENWv7yN/phZJJlx/6EwsI5/7IhS+HTULs75O+oOzxM+/MWVXPO+QtfOMHTuNMtL3Pud282Mn3z1RIR+fN/jK/ffwoqr43WGRtlI899toc+qZZex5JyZw+VRoZ8UK+dcyHeuvTFy39gt6u6oVZzzd9QdYagsv/1qbv7Mu4s67+DePoJkyNwLZmTP+ClQ4jkIQi69aBaXXDyivRO3sHvWuWHIJ695G7dfcUnOOlJx/xmnTeT9CxYzsSX3amWZ1MKgb7GcsgXd65iqd/56EjLAMCIFPLydrK2UZHID2Srb51/Ovsxja92f/8tF15zPhEm5pI8FHRxCk6adQTKkt9usf5Frxk9s6CcRIINJOjsPM3vm1LztyRwADkKlfesurp0/l/F+M5JglNpZZuhnc58JUy0646zhstEzfkzyU0oFO9dHn3/QVxGUROCRCIwEc+ldqdZeyP8SLSTCbKnpuBk/UeWC4sfITxR9ZWGIDg4WtNU6Ve/8sdLApWrwZ2nVnKyt9HPLbatUSedTbKtU5dBh/f2RCkr/zDKvI4+tdWs20tTcxNtuiHuKto5UxNhNJNEwZM8r+1FVZqecfxHqjaIgCaX71T7OsdM945Z2jEKBp7fuZFxzM8sXzENCc1PwBqNDP5t7ewBMsldWpq8O57Clq3kW05UjMg8mC1dQ+69HMvQIwviQVC6ZB7FWUtc3GPoMhH5JNxOPEL/q1PwNIrJCRHaISIeIZOnri0iriDxq968XkTlp++635TtE5IZ8Nu2iMett+aN2AZmcVL/zT1GKA4o7p8y2inZmcceX8pRdTls5qylBz7+cdedoZxTb1+/iaN/xmGzftOzeNMKBgeG4/5yF2Us6RpLxvRCgc+9hJk4cx6SJ42JPy1zcBSBsgo2d++k53s+NixcMJ/h6CSLVC44ODtB57CgXTTsra5/nm+saaeJIunCpXTP64zILqKhGL+sYj0baCtQrwYXbm1IO6adKISI+8CPgRsziK7eJSOZqQ3cAR1T1POB7wLftuRdiVv5aBKwAfiwifh6b3wa+Z20dsbZzUjvOv0jK6XzKTrVKJ5dTz78ADf6yEFNPGCrPP/kfPL+Igc9Ekt2buzh+uJ+hwYiZNFFx/wivs3d3Lzu3HWBy2jTRQpK98I1rfGrLDloyrzvj9FTc/1/d+xhIRoShvNK0cqJIPaVHUayefy6SYXHuqMr8fSZXAh2q+qqqDgGPACszjlkJ/Nq+fgx4lxjNkpXAI6o6qKp7MOueXxln056z3NrA2nxvvguUanaSItIH7KtA1WcAr1Wg3krRSO1tpLZCY7V3oaqedjIGROQZzGdWCG1A+pziB9Uu4i4itwIrVPWT9v1HgSWqeldaXVvsMd32/W5gCfA1YJ2q/s6W/wJ42p6WZTPt+PNs+SzgaVVdnOviq3rYX1WnVaJeEdmgqldUou5K0EjtbaS2QmO1V0Q2nKwNVV1RjmupBeo27ONwOBwVZD8wK+39TFsWeYyINAGTgMM5zo0rPwxMtjbi6srCOX+Hw+EoPy8A8+0snBbMAO7qjGNWA6vs61uBZ9XE4VcDH7azgeYC84H/xNm05/zd2sDa/FO+C6zqsE8FebDSFzDGNFJ7G6mt0FjtrZq2qmpSRO4C2jHL8fxSVbeKyDeADaq6GvgF8FsR6QBexzhz7HF/AF4BksCdqhoARNm0Vd4HPCIi3wRetLZzUtUDvg6Hw+E4Nbiwj8PhcDQgzvk7HA5HA9Jwzl9EZonI30XkFRHZKiJ32/KpIvJXEdll/06x5SIiP7Rp05tEpPhVQqoAmyH4ooisse8j08FzpZzXAiIyWUQeE5HtIrJNRJbVc9+KyL32e7xFRB4WkbZ66lsR+aWI9No58amyovtTRFbZ43eJyKqouhqNhnP+mAGUz6vqhcBS4E6bIv0lYK2qzgfW2vdgUqnn2+3TwE/G/pLLwt3AtrT3cengkSnnNcQPgGdU9XzgEkyb67JvRWQG8DngCpvQ42MGDeupb3+FkThIp6j+FJGpwFcxCVFXAl9N3TAampRqY6NumClR7wZ2ANNt2XRgh339U+C2tOOHj6uVDTPvdy0mBXwNJjP+NaDJ7l8GtNvX7cAy+7rJHieVbkOB7ZwE7Mm83nrtW2AG0AVMtX21Brih3voWmANsKbU/gduAn6aVjzquUbdGfPIfxv7svQxYD5ylqgftrkNASi0r9R8sRbctqyW+D3yREYHg04GjqpoShUlv03B77f5j9vhaYC7QBzxkQ1w/F5EJ1Gnfqup+4AGgEziI6auN1GffplNsf9Z0P58qGtb5i8hbgMeBe1T1ePo+NY8HdTEHVkRuBnpVdWOlr2UMaALeCvxEVS8D3mAkJADUXd9OwYiAzQXOBiaQHSKpa+qpP8eahnT+ItKMcfy/V9UnbHGPiEy3+6cDvba8kDTtauYq4BYR2YtRAVyOiYvHpYPHpZzXAt1At6qut+8fw9wM6rVvrwP2qGqfqiaAJzD9XY99m06x/Vnr/XxKaDjnLyKCyX7bpqrfTduVnmqdnh69GviYnUmwFDiW9pOz6lHV+1V1pqrOwQwGPquqHyE+HTwu5bzqUdVDQJeILLRF78JkSdZl32LCPUtFZLz9XqfaW3d9m0Gx/dkOXC8iU+yvpettWWNT6UGHsd6AqzE/EzcBL9ntJkzscy2wC/gbMNUeL5gFFHYDmzEzKyrejhLb/g5gjX09D6MX0gH8EWi15W32fYfdP6/S111kGy8FNtj+fRKYUs99C3wd2A5sAX4LtNZT3wIPY8YzEphfdneU0p/AJ2y7O4CPV7pd1bA5eQeHw+FoQBou7ONwOBwO5/wdDoejIXHO3+FwOBoQ5/wdDoejAXHO3+FwOBoQ5/wdDoejAXHO3+FwOBqQ/wPD3wayeAd9rQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "int_obs = griddata(X, yobs, (grid_x, grid_y), method='linear')\n", "\n", "int_exp = [griddata(X, yexp[i], (grid_x, grid_y), method='linear') for i in range(5)]\n", "\n", "plt.contourf(grid_x, grid_y, int_obs, levels=np.linspace(0, 1))\n", "plt.colorbar()\n", "\n", "plt.contour(grid_x, grid_y, int_obs, levels=[0.05], colors='w')\n", "for level in int_exp:\n", " plt.contour(grid_x, grid_y, level, levels=[0.05], colors='w', linestyles='dashed')\n", "\n", "plt.scatter(X[:, 0], X[:, 1], c=yobs, vmin=0, vmax=1);" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "application/scrapbook.scrap.json+json": { "data": 225, "encoder": "json", "name": "number_2d_successpoints", "version": 1 } }, "metadata": { "scrapbook": { "data": true, "display": false, "name": "number_2d_successpoints" } }, "output_type": "display_data" } ], "source": [ "sb.glue(\"number_2d_successpoints\", len(X))" ] } ], "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.8.6" } }, "nbformat": 4, "nbformat_minor": 4 }