{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# ShapeFactor" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import logging\n", "import json\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "import pyhf\n", "from pyhf.contrib.viz import brazil\n", "\n", "logging.basicConfig(level=logging.INFO)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def prep_data(sourcedata):\n", " spec = {\n", " 'channels': [\n", " {\n", " 'name': 'signal',\n", " 'samples': [\n", " {\n", " 'name': 'signal',\n", " 'data': sourcedata['signal']['bindata']['sig'],\n", " 'modifiers': [\n", " {'name': 'mu', 'type': 'normfactor', 'data': None}\n", " ],\n", " },\n", " {\n", " 'name': 'bkg1',\n", " 'data': sourcedata['signal']['bindata']['bkg1'],\n", " 'modifiers': [\n", " {\n", " 'name': 'coupled_shapefactor',\n", " 'type': 'shapefactor',\n", " 'data': None,\n", " }\n", " ],\n", " },\n", " ],\n", " },\n", " {\n", " 'name': 'control',\n", " 'samples': [\n", " {\n", " 'name': 'background',\n", " 'data': sourcedata['control']['bindata']['bkg1'],\n", " 'modifiers': [\n", " {\n", " 'name': 'coupled_shapefactor',\n", " 'type': 'shapefactor',\n", " 'data': None,\n", " }\n", " ],\n", " }\n", " ],\n", " },\n", " ]\n", " }\n", " pdf = pyhf.Model(spec, poi_name=\"mu\")\n", " data = []\n", " for channel in pdf.config.channels:\n", " data += sourcedata[channel]['bindata']['data']\n", " data = data + pdf.config.auxdata\n", " return data, pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "INFO:pyhf.pdf:Validating spec against schema: model.json\n", "INFO:pyhf.pdf:adding modifier mu (1 new nuisance parameters)\n", "INFO:pyhf.pdf:adding modifier coupled_shapefactor (2 new nuisance parameters)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "data: [200.0, 300.0, 220.0, 230.0]\n", "expected data: [100. 100. 120. 90.]\n" ] } ], "source": [ "source = {\n", " \"channels\": {\n", " \"signal\": {\n", " \"binning\": [2, -0.5, 1.5],\n", " \"bindata\": {\n", " \"data\": [220.0, 230.0],\n", " \"bkg1\": [100.0, 70.0],\n", " \"sig\": [20.0, 20.0],\n", " },\n", " },\n", " \"control\": {\n", " \"binning\": [2, -0.5, 1.5],\n", " \"bindata\": {\"data\": [200.0, 300.0], \"bkg1\": [100.0, 100.0]},\n", " },\n", " }\n", "}\n", "\n", "data, pdf = prep_data(source['channels'])\n", "print(f'data: {data}')\n", "\n", "init_pars = pdf.config.suggested_init()\n", "print(f'expected data: {pdf.expected_data(init_pars)}')\n", "\n", "par_bounds = pdf.config.suggested_bounds()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "initialization parameters: [1.0, 1.0, 1.0]\n", "parameters post unconstrained fit: [1.00004623 1.99998941 3.00000438]\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/srv/conda/envs/notebook/lib/python3.7/site-packages/pyhf/tensor/numpy_backend.py:334: RuntimeWarning: divide by zero encountered in log\n", " return n * np.log(lam) - lam - gammaln(n + 1.0)\n" ] } ], "source": [ "print(f'initialization parameters: {pdf.config.suggested_init()}')\n", "\n", "unconpars = pyhf.infer.mle.fit(data, pdf)\n", "print(f'parameters post unconstrained fit: {unconpars}')" ] }, { "cell_type": "code", "execution_count": 5, "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": [ "(\n", " obs_limit,\n", " exp_limits,\n", " (poi_tests, tests),\n", ") = pyhf.infer.intervals.upper_limits.upper_limit(\n", " data, pdf, np.linspace(0, 5, 61), level=0.05, return_results=True\n", ")" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "expected upper limits: [array(0.74138115), array(0.994935), array(1.38451391), array(1.92899382), array(2.59407668)]\n", "observed upper limit : 2.1945969322493744\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmUAAAGxCAYAAAAjyMGDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAB+H0lEQVR4nO3dd3gVxdfA8e/cdEIIvfdepVcRQRBpggooCIqCIK+FZsMG2BsiqKAgooDYQZAmoKD03ov03lvokDbvHxN+IiZw783u3U1yPs+TJ7l3d2cOCZrD7Ow5SmuNEEIIIYRwlsfpAIQQQgghhCRlQgghhBCuIEmZEEIIIYQLSFImhBBCCOECkpQJIYQQQriAJGVCCCGEEC4QsKRMKTVGKXVMKbUxheNKKfWxUmqHUmq9UqpaoGITQgghhHBaIFfKvgaa3eB4c6BU0kcP4LMAxCSEEEII4QoBS8q01vOBUzc4pQ0wThtLgaxKqXyBiU4IIYQQwlnBTgdwjQLA/mteH0h67/D1JyqlemBW04iMjKxetmxZ24Lavn07Z8+eJTw8hAoVKgGKhIQEgoKCbJtTCCGEEOnTqlWrTmitcyV3zE1Jmde01qOAUQA1atTQK1eutG2utm3bMmnSJO6+J44fv9vJpUstyJNnMi+99BL9+79i27xCCCGESH+UUntTOuampy8PAoWueV0w6T1HjR8/nqBgxfY4gDPExX3Hc89d4vbb3wEeYvfuL6hfvx4rVqxwOFIhhBBCpGVuSsp+BR5OegqzDnBGa/2fW5eBlilTJoqXz87aZdB7AJw4Aa++CnXrXgS+4ciRHpw7t4zcuT8FtrBz5062bdvmdNhCCCGESGMCdvtSKfUd0BDIqZQ6AAwEQgC01p8DM4AWwA7gIvBooGK7kaNHj7Ln7xiIhU/fhxwhMODVf47XrQvr1iUC44BxvPFGHiZNOsvRo0eJiIhyKGohhBBCpDUBS8q01h1vclwDTwYoHK9lypSJiPBQ4iIukdgeMje+8fnvvnuUBx6AiIiyQHcefngj9eo1oWfPngGJVwghhBBpk5tuX7pSVFQUZ878SeYXgIIw8C9YvwvOnUv+/Lx5oXlzgENcufIaR45M4syZ4cAytNacPHkycMELIYQQIs2QpMwreSiWFbgM589AzSowaNDNrwoLg9mzNc8/vxGow6xZtSlcuCDLli2zNVohhBBCpD2SlHmhceNH2Pg08BEQBrH1oXQT769XynwuWXIFXbtepmrVt4ENbNy4kZiYGOsDFkIIIUSaI0mZF/LnL0imaOAKcByoDe9vh0txvo1TsiR88gmEhv6K1rfQsWN9WrRoZEPEQgghhEhrJCnzwvjx4+n3RTbzYp/5tOs4NLgfRo/2b0ylYPz4M7z33jrgIeLj/+bbb78lMTHRkpiFEEIIkbZIUualGuXzQGb+l5QRBKvWw6rN/o9ZpQrcdpsGvuGnnyrSqVMn/vhjZuqDFUIIIUSakybbLAXaDz/8wEMPb4VI/knKFOhOsK2kNXN06JBAzpzQpEkvIIRly6KpWLEikZGR1kwghBBCCFeTlTIv5MmThxIloqEkpvzt1TuMQTB3N7z/C6xenbo5lII77wSldnHhwl20bNmAxx7rnLpBhRBCCJFmSFLmhYYNG7Jlc2/y3A9U5d/ftUR4qQf0fca6+SIjYfLkWAYO/B0YwdmzMRw4cMC6CYQQQgjhOpKUeS0PRaOBE/xzCxPAAwn3QO3e1s5Wvz6ULXseeJKBA8tRsWI5Tp06Ze0kQgghhHANScq8cObMGXLleoHlPYGpwLTrTigIIzbDkfNgx8OTvXod4f33L5A9+ydAIgkJCdZPIoQQQghHSVLmhaioKGrXrkiZ+kBB4Bhw6d/nXLgItRrB4MHWz1+sGPTooYFBrF9flwoVyrBu3TrrJxJCCCGEYyQp84LH42HatLE8+xFmsz/A/utOCoEDsXDW5jJj8fHLyZlzH/nybbV3IiGEEEIElCRlXstDkWggN+a7tu+/Z+j7YHUhe6OoVg0WLIgjd+4OaP0S7777NsePH7d3UiGEEELYTpIyLz300JO0qgjMAvIBe5M/b+YOeG8c7NhhXyyml6Zm06Z3GDToFX74YaR9kwkhhBAiIKR4rJcaNGjA2bM/My3yMonZMIVkk3PJlMjY1QVG2pwrVawI69drSpb8CKjK3r0VKVSoEB6P5NpCCCFEWiO/vb3UvXt3pkypRNE6mFuYKSVlEZDYGWp2C0xcpUuDx3OKc+dacdttlejV6+nATCyEEEIIS0lS5oPExFzk00A8sAhIqe9lIXhtEZy/DFoHJrbMmaFfv3M89th24EpgJhVCCCGEZSQp89JPP/1ESMgMFvUBYoA1SR8pOHAASpaHH38MSHgoBX36QJUqc4AmjBnzCbt27QrM5EIIIYRINUnKvFSuXDl69KjKnd2BCKAwpixGSiUwMsPJCEgMDViI/3Pq1EKef74vH3zwauAnF0IIIYRfJCnzUsWKFfnssy507YXZT1YYuAykVI3CA/H3w/yQgIX4P9mzw9KlCXz00e/ASukAIIQQQqQBkpT5IDExJ5kvYW5fFk56M5l6Zdf6YiW88RGcPGlvbNcrWRLCw49x6dLt3HFHFUaNGhXYAIQQQgjhE0nKvBQfH094+CN0aY6pVZYNiAbO3/i6hOMw4BkYP97+GJOTmHiRrFk3kj37SmcCEEIIIYRXpE6Zl4KDgxkw4P+AT3h7H1xSQG9untbmBrpDtXvsjjB5kZEweTIo9QWQgy1bHqZ06dIEBQU5E5AQQgghkiUrZT545ZUBvPIKlKya9Ia33738MPAvcGprl+kAAPv2vUvNmpUZNGiQM4EIIYQQIkWSlPkgISErBw4osp4E4oBzwJfApptf++dcyF8I9qbQnikQCheG99+P48knpVemEEII4TaSlPmgc+eHqFVLseBVzFOXmYAjwJ/AnptcnAPicsDFi/bGeDNPPAF5845E6/5MnDgRHajqtkIIIYS4IUnKfPDQQw/x7LN5eOx9zEb/IKAQJkHbeJOLo+F0OzgYZneU3pky5T3atWvHlClTnA5FCCGEEEhS5pMWLVrQr18lWrfGFJAFKJL0eQspF5K9xsszwA15UJs28Ouv0KbNeqdDEUIIIQSSlPkkPj6ew4czc3EvcCjpzav1yi5c894NLP8O2raDo0ftidFbSsHdd4NSAzl06CWGDRsmtzKFEEIIB0lS5oOZM2eSP/8kBvYC/kh6swBQIunrv70YpB6UexHy5LElRL+MGvUOr7zyAvv373c6FCGEECLDkqTMB1WrVmX48HsYMgSy35v0ZijwEFAU2OrFIJlhYxDM2G5XlL4bMABWrbpC4cKTnQ5FCCGEyLAkKfNBwYIFeeKJe2jRAkqXve5gccyGfy/bKT3WF3r0sDhAP3k8ULo0QG9++KEHX375pdMhCSGEEBmOJGU+OngwiE2bIPhv4FLSm7uAuUlfe3MLEzh8BnafgkQvHg4IFK1h3LjRTJjwCYluCkwIIYTIAKTNko/KletJ48awcDLwCOa2Zd6kg5GYW5i3ejFQYziZ759q+26gFHz/vSY0dBcez2agotMhCSGEEBmGrJT5aOTI9+ndGwZMwGzyB1NENhdmf9k+btqkHAAFa47AyN9hzx57YvVHVBSEhZ3j0qUW9O7dnWPHjjkdkhBCCJEhSFLmo44du9OwoeL2ukDINQeK8E8yts3LweLgqXvhpZcsDdES27bt58svx/D77zOcDkUIIYTIECQp89GxY6fZuDGKA2uAa5+gLIzph5kZ757CBAiBhDbQyCUb/q9VuTLs3JnIgw9OB6R+mRBCCGE3Scp8NHDgQBo1usDoYcCCaw4UB1oCpYCdQKyXA5aGj7eYTfZuY2qp/czq1V358MMPnQ5HCCGESNckKfPRY489xldflWP8OCjS85oDmYGaQCUgHpOYeWnjDqjbDNautTJS63z11dd8/PE7nDlzxulQhBBCiHRLkjIfVa9enVatylOkCJQscN3Bs0kfYXhdGgOAEFi9ArZssSxMSw0ZAitWnCE6erXToQghhBDpliRlPjp//jzLlmk2bYIL84Fz1xzcDkzGbPrfBiR4OWgExD0FUdWtjdUqISGQO3c8Wt/H++8/x7Zt3j7JIIQQQghvSVLmo7Vr11Knzk/MmQNLvwSOXHPwanPyLJjCsr60kgyCdxeCm9tPHjsWwwcfDGHMmBFOhyKEEEKkO5KU+ahixYpMnfoU7dvDyPlAyWsO5sTULLsCBOHbLUxg0Y9QshQcPmxVtNbKkwdWr07knXf2OB2KEEIIke5IUuajrFmz0qrVnRQoAJWLAddW5FeY1bL9QDFMaQxfnqosAyXamgKublWoECg1haNH32DmzJlOhyOEEEKkG5KU+WH58hNs2gRr5gBrrjtYBIjBtF86DfhSED87/F0adl+wJExb9ekzkM6dO3D+vDftC4QQQghxM9L70g/t279Ko0Zw7Bh4NkFi1WsOVgYqYFbNfsfcwszj/dga6PUJ9KgAHTtaGLTFPvpIc/x4JJkzxzkdihBCCJEuSFLmh+++G0+uXI3Jnx/qjIWNx685mOmarwtgbmHe7tv4f/4AJ6OgQwd3NSy/Vt68kDfvYaALa9YMokqVqii3BiuEEEKkAXL70g/16t1BqVKZiYyE4tmTOWEL8BtQFjgE+FpztQ3UG+DehOxa8+ZNpVq16vzwww9OhyKEEEKkaZKU+WH79u1Mn56ZXbtg/0/AyetOOAYsxewrA+97YV4VBeM2wZFzkOBtrTOH3H47DB3qoU0bH+7RCiGEEOI/JCnzw1dffcU99xzlzBnYOB04dd0JV+uVXQSy43NpDIBLp6DSLfDjj6kK1XYeD/TunUhExMPExh6Wjf9CCCGEnyQp88Pjjz/OihW3U6kS/LwW04T8WgUx39l9mFuYe4DLPk6SGc5khtDI1EYbGPHxB2jUqDyPPfaY06EIIYQQaZIkZX4oUqQIVaqUITgYSiS3pywEs8l/LyYpS8S0YPKFB+Lawa4cqQw2QIKD4f77Y2jfXp4dEUIIIfwhSZkfYmJi+Omnkxw8CDMnAH8lc1LRpM8FgEj8uoUJMGQhTPzFv2sDrXdvaNv2e2AhWvtSNVcIIYQQkpT54cCBA9x//88sWgTr10DooWROugN4DNNuqQxmpSzW97mOLIB298HKlamJOJAS+Pnnttx55x3ExUkNMyGEEMJbkpT5oVSpUqxfP5gWLWDcOKj2XDInXS1noYFbMAnZFj8mqwz5/g+qVvM32sDzeI5x+fLfnD592ulQhBBCiDRDkjI/hIWFUalSLTJnNq+LZU3hxKnA95inMbMC6/yYLBQO54GJ/iR0DrnvPpg//wi5c693OhQhhBAizZCkzE9Tp25i7lzYvh3WfAwcTuYkD7ALs1pWOelrXwvJJnn2TXguuRU5l/J44Ny5R3j99Ze5cuWK0+EIIYQQridJmZ9eeeVThg0zycfpPUByTcSLAHHAEUxSBv6tlgH798D8FZCY6N/1Tli8+CCvvfYOc+fOdToUIYQQwvWkfoGffv11OlFRpcmePZYf/oKGY5M56WoR2b1APUyStha4jX/2nHmrKeQua5LAtOKuu2DrVk3JkmkokxRCCCEckoZ+xbtLkSJFyJ7dtBYqlVItsSxANkwRWTCrZaeAA35M6IHp22DVbrjsayFaB5UsCdCdTZsWc+nSJafDEUIIIVxLkjI/rVmzho8/Nl9P+ByCJ6dwYk3+qVlWHrM2uda/OfVpqFsBxia3Kudi+/cfpnr123jrrbecDkUIIYRwLUnK/DRnzhx6997P+fMQGwuRKdVKrQfUSfo6HJOYbcTsNfNVVqA2VKjqx7UOKlQIPvsskd69SzodihBCCOFakpT5qWfPnpw8+SCRkfDyy9D81RucfJl/nrqsDFwBtvoxqYK4hrDIjyK0Tnv0UciV6wW0Pi5FZYUQQohkSFLmpyxZspA9e2FU0ob90sn1wARTDmM48EfS62KYvWZr/Z976GwY8Zn/1zslIeEY99xTmb59+zodihBCCOE6kpT56fTp0wwZsplNmyAmBr7rS/LlLhRQCPMEJpjveGVgJ3DWv7mPLIZeveCAPw8MOCgoCMqXP0zp0uedDkUIIYRwHSmJ4aeLFy/yzDO/kjkzPPYYZM9Kyt/NIsBmIAazL6wysADYANzqx+Q1oEITKFjQj2sd9s47ANOAE0BOZ4MRQgghXCSgK2VKqWZKqa1KqR1Kqf7JHC+slJqnlFqjlFqvlGoRyPh8kS9fPmJiJtO9u6kdNus3oEIKJxdJ+rwn6XNOoCDmFmZKDwjcSASsvwKL9t38VHc6yZw5nfn888+dDkQIIYRwjYAlZUqpIMzuquaYZxA7KqXKX3faK8CPWuuqQAdgRKDi85XH4yE6usT/9pRFh0PuyBROzo158nLvNe9VAY4Dh/wMIBEe7Hx15Snt+fLLWYwYMZj4+HinQxFCCCFcIZArZbWAHVrrXVrrWEyr7jbXnaMx2+ABovE/ZQmIr776g3HjzNdffAExbwMJyZzoAe7FlMe4qgIQhN9tl/DA/tNwPI1uz/rsM1i+PIxguYEuhBBCAIFNygoA+695fSDpvWsNAjorpQ4AM4CnkxtIKdVDKbVSKbXy+PHjdsTqlbFjf/lfIdfChaFodSClchVlgFzXvI4AymL2lfm5WKTbAQ38u9Zp2bJBePhm4uIG8/fffzsdjhBCCOE4tz192RH4WmtdEGgBjFdK/SdGrfUorXUNrXWNXLly/WeQQPntt9/4/ffcgOnz2G0QJtlKTjwmAbv2ickqwCVgm/8xfLkGVq4D7c/eNBfo1u1V7rjjdi5cSK6juxBCCJFxBDIpO4gpDnFVwaT3rtUN+BFAa70EsxPLtY/ohYeHo1Te/70unQNIqfe2B5gKrLnmveJAZlJVs+zseqhZBebN838MJ/XuHc+IEUXJlCmT06EIIYQQjgpkUrYCKKWUKqaUCsVs5P/1unP2AY0BlFLlMEmZc/cnb2Lx4sX06xfDlSvm9YvtgekpnOwBCvPvzf5BwC3ADsDfvWHFIfs9ULmKn9c7rHp1uOee5SiV0jdOCCGEyBgClpRpreOBp4BZwBbMU5ablFKvK6VaJ532DNBdKbUO+A54RGv33pjbtGkTo0Yd4tQp87r9faCK3OCCopjyXNcmYFUwq2sb/AwiBE5VgflH/bzeJb799lGaN29KQkJyT0oIIYQQ6V9A95RprWdorUtrrUtord9Kem+A1vrXpK83a61v1VpX1lpX0VrPDmR8vurWrRvnz/ciXz7z+vXXoOiNNt5fTdiuXS3LDeQnVbcwAQaMgZEjUzeGs05w6dLfnD592ulAhBBCCEe4baN/muLxeDBZ1T9KRZPyvrL8QAhw+Lr3KwNHk3nfBxt/h7fehbS60NSxI8ybd5icOY85HYoQQgjhCEnKUuHs2bP06fMHf/1lXs+YAXO6YRKs5AQBvYEm171fCdOiaXkqgmkK9d8z/SXTIqVAqXjOnOnBuKvF34QQQogMRJKyVAgODmbMmIVs2WJely8PzR4l5bIYYJ62vF4mzGrZevzf8J8JJm2DY+fTbnkMgE8/XcSjjz7Cjh07nA5FCCGECChJylIhU6ZMnD27kJ49zeuiRaHvi5im4yk5D/wM7Lzu/bqYbgAr/I/nygmoUg1mzvR/DKf16wcrVkRTsmQOp0MRQgghAkqSslT7956yYlm48WpXOPA3sP2693MCpTFJWZyfoWSBUwoS0/BKWUQEVKsWAwzk4sWLTocjhBBCBIwkZak0dOj3vPHGP6+f7AzquxtcEIwpm7s3mWN1gYuY25j+CIIrD0JccT+vd5ExY4ZTvHgRTl2tNyKEEEKkc5KUpdLq1etZvjzkf6+ffBIK3nWTi4oAR4DL171fFMgLLCHlJzi98PESWOdvo3OXqFUrkVatIkhMTMU3QgghhEhDlItrs3qlRo0aeuXKlQ5HUQ5zT9Jo9yNM3HKD03cDY4EHMbcsr7UemJTCMW9NhSw74NBBiIz0cwzXmMN/H1cVQggh0ial1CqtdY3kjslKmSX+2VeWmAi54zCNxlNSEMhH8qthFYAozGqZv2rCbX0gPDwVY7jEzp1P8t577zgdhhBCCGE7ScpSaf78+bRrt5OTJ83rPXvgs86YRlIpCQEeB8omcywIqI1ZTfO3mGxemB8GF+P9vN5Ffv55G2+88Rr79u1zOhQhhBDCVpKUpdLZs2fZvPkiZ86Y10WKwAvvYvaH3UxC0sf1qmMSt6X+x3XuPHR5AVakosSGG/TuDdu3Z6Zw4axOhyKEEELYSpKyVGrVqhWbN/eleNITj0FB8OzTQPabXHgYeA+zIna9CKAqpkn5Wf9jmzIKpk3z/3o3CA+HfPlOAm9z4sQJp8MRQgghbCNJmSXy/OuVvgBZbtbCMQcQD+xJ4XgdzJ4zf1svhULik9Cgi5/Xu8xrrw2mQoWynD2biixVCCGEcDFJylLpypUr3HvvGL7//p/3PvgAzo3ixmUtQjGb/ZOrVwZmpa0csBKI9TO4zPDpirTddumqli0TePrpvISEhNz8ZCGEECINkqQslUJDQ9m//ywxMf+898gj0ORV4GbJUFHgICknXXUxtczW+h/flO+gTHmI9Texc4kaNeCVVzYREbHK6VCEEEIIW0hSlkpKKVau/Ot//S/BNCa/vSHmScobKYJZTTuQwvFCQAHMhn8/a6jqzBCcl/89iJDWLVjQnQEDXnU6DCGEEMJykpRZIhdmd76RmAgJe4Cb7SsrDDQAolM4roB6wClgq5+hlYITLSBLNj+vd5k//vibr7/+jJhrlyaFEEKIdECSMguMGDGCZs3U/14rBR88hdkPdiPhwB2YTf8pKYtJ2lJRTPb4RRg5Dw6ktCKXhrzwAvz9dxhZs4Y6HYoQQghhKUnKLODxeFAqnKttGpWCyVOB+l5cHAvswjyJmZwgzJOY+zD7z/wRB/3awGuv+Xm9i0REQKZMh0hMfJ+9e1N6SkIIIYRIeyQps0DPnj2ZOfNuPNd8NxvfDoUKenHxLmAcN064qgJhwGI/AwyBhFbQ8lE/r3ehRx55i0aNGhCb1p9gEEIIIZJIUmaZwv96tX8/RG0B4ry8bM8NzgkHagKbgCN+hlcBJt9sj1sa0rVrPG+8UZbg4GCnQxFCCCEsIUmZBQ4fPkytWuP55Zd/3lu0CDZ/gdmkfyOZMLVnb3Yn7lZMcva7/3F+Nxd6P8P/brOmZQ0bQqdOf+Lx7Hc6FCGEEMISkpRZIDo6mmzZshMW9s97zZrBi98DOb0YoAiwn+T7YF4VgdmjtoMbr6rdQOwB+GwEbNvm3/Vuo3UsX3/diVGjRjkdihBCCJFqkpRZIFOmTMyaNY4WLf55L2tWqF+Nm9cqA1NENg44dJPzagNRwBxuXpg2ORWg0CAoU8aPa11IKZg0aRE//zwWnR7aFgghhMjQJCmzTOH/vLN/OfC3F5cWBx4D8t/kvBCgIeahAG/GvV4w7LoMf+xOH62XAMaPh1mzcqGUuvnJQgghhItJUmaRZ54ZxJ13/nvT+bdfgPLmiclwoCDerapVwdwS/YMb3+5MSTx0vg/eecePa10oOhqUmsLFi39x+vRpp8MRQggh/CZJmUWKFy9OpUr/Ls3//XdQqq+XAxwBZpJyvbKrgjAFZ08A63yNEgiGYwpUuB/XutTly1C+fDNeeuklp0MRQggh/Cb1BCzy5JNPAtMxmZWRLx+UzQvbvOk7GQMsA8ph9pjdSDlMT8x5QCXMbU0f6HsgoYZv17hZeDg8++xlqlQp7XQoQgghhN9kpcxS/95XdvgwnPoNOOnFpUUwvS53e3GuAu4EzgHLfQwxyahVsG69f9e60VNPQf363zodhhBCCOE3ScossnbtWvLlm8Dv19QRO3sWFo4FDnsxQASQF+/LXRQFSgILgEu+RGrs/w2qVoH01Kno/PmVDBp0P+vW+XNfVwghhHCWJGUWyZMnDy1bViPHNc3FS5aEGRuBil4OUgw4gOmH6Y0mwGVgkQ+BXlUBKnWH3Ln9uNal4uPhk08m8ttvM5wORQghhPCZ7CmzSL58+Rg9ehBmF74RFARV/lspI2XFgPXAaUyV/5vJi9lTthSoBWTxYa5o2JgVjl2BIhE+XOdiWbPC1q2J5MyZjjJNIYQQGYaslFmq0H9aGC2ZA6Hzvby8BPAM3iVkV90BJAJ/+XBNksR46PUezJvn+7VulTMnwGscO7ZfCsoKIYRIUyQps1DLlk/TvPm/31uwALP65U1+4MFs4vdFNqAGsBpTJsMXCqaNhDFf+Xidyy1btp+iRUswffp0p0MRQgghvCZJmYXuvrsNbdpk/td7gwdD28/xPtnaDAzD7BXzVgNMWYw/fLgGwAOJXaDlcz5e53LVqkHPnkFUqFDE6VCEEEIIr0lSZqGePXvyxBP/biwZFASlc6RwQXLCMXvK9vlwTWagLrDFx+sAssKo1T5e43IhITBkyGWKFZvidChCCCGE1yQps1hsbAESrml/dPEiLBwBbPVygEKYqv3e1Cu7Vj3MRv9fuXlXgOvMmwkVKsMlP0pruNmePR8wYMALJF6/0U8IIYRwIUnKLDR16lTCw39lw4Z/3gsPhw0LgVNeDhKCScx8TcrCgFaYfWULfbw2HM554NgxH69zucWLz/L++0PYcO0PRAghhHApScosVKFCBQYOvPNftco8HtixG3N70VtFMb0wL/oYQGlMiYz5gC8JVjG40BHyFPBxPpd74AHYtSuSypWLOh2KEEIIcVOSlFmoePHiDBzYjUKF/v1+VBjky5z8NckqBVTH59uQADTDrJr9iimV4aVTl2Ds8vS1WhYUBPnznwGGcSm93ZsVQgiR7khSZrHY2LzExPz7vVmz4PJXeF+pvwBwN74Vg70qEpOYHQBW+nBdPDzVDAYN8mNOl+vf/21q164he8uEEEK4miRlFitX7mGefPLf7125AiFX8O12ZCJw3M8gbsEUov0diPHymmCIbwi3tvZzTherVesKrVvnJDbW26xYCCGECDxJyiz2yisD6NTp39/W1q3hze+BrD4MtBQYDpz3IwiF2fSvgel4V7gWoDosSrj5aWnNfffBm2+uJzz8itOhCCGEECmSpMxijz7ajRYt/rtjvqKv7Riv1j3d42cg2TAtmLYDG72/bPx8+HCon3O6WgxLljzDn3/+6XQgQgghRLIkKbNYXFwce/bk+k8PzAmDAV9qmebFbNj3tTTGtWoD+YGZeH3r9PwGeO4Z2LkzFfO6UGIiPPbYVwwc+LLToQghhBDJkqTMYqNHj6ZYsdUcOfLv96MzQ1SkDwMFYVbLUpOUeYDWmJZNs7y8pjJUfgdKlEjFvC7k8cDPPycyffrtTocihBBCJEuSMovdcccdfPllMyKvS8Deegsa/J+PgxXDFJ09k4qA8gK3AusAb1a/wmHtJVh35OanpjXlykHmzCPQ+jRae7vRTgghhAgMScosVqZMGbp2bUF09H+P+byvrDzwIJAplUE1AHIAU/GuLEcs3P8gjB+fynld6OjRM9SrV4kffvjB6VCEEEKIf5GkzGJaa/bsCePQoX+/f/IkfPM44Evz72hMlf6QVAYVgrmNGQP84d35O3fBwXS4WpYrF2TPfpTgYHkSUwghhLtIUmaDihX78MEH/34ve3YoXwFT3NUXx4DFeF/WIiVFgFrAMmDLTc5VkPAI5GmSyjldyOOB6dPjaddul9OhCCGEEP8iSZnFlFJ89dUndOly/fswdRIEl/NxwD3AbOC0BcE1xXQLmAycvMm5CkathhMnLJjXhRITh/Lrr9+SkJAOC7MJIYRIkyQps0H79t2oUiX8P++HBUOJrPi26lUs6XNqnsK8Khhoj/mp/8hN95ct/RaKFYcLFyyY22VmzjxLmzadmDLFlzolQgghhH0kKbPB8ePHWbIkx3/e//Zb2P4MvlXpzwlkxpqkDExXgfuAo8AMbpwgloBKbflPzbX0oHlzmDw5kjZtGjkdihBCCAFIUmaLESNGcOutB7m+1WK5clCrDb6tlCnMatkeH6+7kVLA7cBabvzgQQHYUhaC/7vol+Z5PNCmzQWCgj5zOhQhhBACkKTMFh07duS335qg1L/fr1oVnnsdyOLjgEUxBWBTU6/sercDxTGrZYdSPi3mIgwYDdu3Wzi3i/z887vcf39bqVsmhBDCcZKU2aB06dI0bVqHkGRKWVTIhUmwfHEL0B/fGprfjAdoi3ka9EfgUgrnXYYPe8Hnn1s4t4vExJxjz54VnDp1yulQhBBCZHCSlNkgLi6OBQsuszuZfWDPPArqGx8HDMFs0rdaJGbj/1ngFyC5vWOZQHeBzn1smN8FunaFZcviyZEjs9OhCCGEyOAkKbNBbGwsDRoM5vvv/3usy8NQyJ+95ZuBr0k+cUqNQphSGdsw9dBSOOerjRbP6xIeDyh1mIsXR7J//36nwxFCCJGBSVJmg8jISH7/ffR/apUBtG8Pt9/rx6AJmM3+h1MXW7JqAxUw1f5TeMpzzLfQ9TEb5nYBraF27Wfp3j2d/gGFEEKkCXbcFBNA48btgeR/yRcLAy7iW0/Lq/XKdmIKwFpJYdowHQF+Brrzn/1rF47DnJ1w9ixk8fVBBZdTCgYOjCNPnupOhyKEECIDk5Uym2zatJ/p0yP+8/65c/B6S2CVjwNmBvJikjI7hAEPAPHAeOD6grG1oHD/9JeQXdWuHdx22ySsvz8shBBCeEeSMpsMHz6chx76b8n8qCh4awimVpivSgD7Abt6aecGHsSU3viGfz8lGgSLD8D6wxAXZ9P8Djt7diuDBj3Axo3pdAOdEEIIV5OkzCbPP/88ixbVS/bYS30huogfg5YCSpNy+QorFAHux1T8/x64NgE7A/Uqms4E6VF8PAwZMonZs2c7HYoQQogMSPaU2aRo0aJAOWDBf45duACFz8MGDxDky6BJH3YrDdwDTAImYspmBAFZIK4kFCgcgBgckD077NqVSM6cpZ0ORQghRAYkK2U2OXPmDN9+eyrZWmWTJ8OGt4CTfg5+LhWBeesWoBnwNzAV0+JJQWwLOPrftp7pRs6cAG9xIT12YRdCCOFqkpTZ5NSpU3Tq9DPz5v332B13QPf38b3dEpgHBD7E2pZLKanDPz0y5/zz9vAFMH9+AOZ3yMSJS8mfPzf79u1zOhQhhBAZiCRlNilUqBCbN39Fx47/PZYvHzzYHvCn0ffVchi7UhGcLxoCNTGFZReat5aMhNZt4IpdDxw4rFYtaNcuCx6P/OchhBAicAL6W0cp1UwptVUptUMp1T+Fc+5XSm1WSm1SSqXZLeXBwcGUK1efiP9WxQAgPAZTDNZXeTDtkewqjXE9BTQHKgK/Y1bqboO734KwsADFEGCFCsGXXx6hYEFZKRNCCBE4AUvKlFJBwHDMr/jyQEelVPnrzikFvAjcqrWuAPQJVHx2mD3772RbLQF8+BYETfNjUIUpjbGLwJXU8mA2/pcEpgHH4LcLEJsQoPkdsnPnS0ycONHpMIQQQmQQgVwpqwXs0Frv0lrHYgoutLnunO7AcK31aQCt9bEAxme5UaO+5rXXkn+8cuBAqNHXz4GLYzoCHPU3Mj8EY0plFAJ+hhOz4b5usHVrAGMIsIED/+Lxxx/jSnq9TyuEEMJVApmUFcCUPr3qAP9tGFQaKK2UWqSUWqqUapbcQEqpHkqplUqplcePH7cp3NQbMWIEy5dXSvZYxYpQp5qfA5fEpLPR/kbmp1CgM+an9DvM+AYWp9TEPB14+23YsKEBYen1Pq0QQghXcdtO5mBMidSGQEfgC6VU1utP0lqP0lrX0FrXyJUrV2Aj9EHu3LmJiiqa7LG4OIjdiH8NxjMDVfGtd6ZVQjHtmKqBToCZc03R1fSocGHIl286sMPpUIQQQmQAgUzKDmJufl1VMOm9ax0AftVax2mtdwPb8K8hkSvs37+f9947xv79/z3m8cCXrwDr/Rz8HLCCf1fcD5Qg4G6gAfz0DdxzD1y86EAcAXDxYgLt2zdl+PDhTocihBAinQtkUrYCKKWUKqaUCgU6AL9ed85kzCoZSqmcmBtlgSr+YLlDhw7Rv/9iNmz477GgIFi4FGjk5+BHgOnAXv/jSxUF3AEhJWD6dGjSBE76WwzXxTJlgsuX9xIff8rpUIQQQqRzAUvKtNbxwFPALGAL8KPWepNS6nWlVOuk02YBJ5VSm4F5wHNa6zT7q75atWqcPfsVLVokf7xmFSia28/Bi2BWrBxOWeNKQ41GsGoV1K8P6bHe6tSpifTufd7pMIQQQqRzAd1TprWeobUurbUuobV+K+m9AVrrX5O+1lrrflrr8lrrSlrrFApKpA0hISFERaXcR3HnTohYBlz2Y/BQoDCBq1eWktqQvSvMmQOHD0PduiS7Mpj2fcaqVX+SmBioOiRCCCEyGrdt9E93vvhiAWPHJn9s2zbY8h3g7wOkJTBlMQLRC/MG5uyE8wrmzjWv69eHSZOcjclqs2ado0aNRkyb5k9xOSGEEOLmJCmz2fjx01IsINuoEYxayL8ff/BF8aTPh/y83iJ6D7RsALt3w5IlUKYMtG0LTz+dfloxNW4Mn30WxR131HM6FCGEEOmU0lo7HUOq1KhRQ69cudLpMFIUGxtLaGhJ/l2i7R/rj0Llz/0cPBG4hGm75KREiN4M2z+DXDkhNhb694ePPoJq1eCHH6BkSYdjtMxw4AmngxBCCJFGKaVWaa1rJHdMVspsFhoayo2Wwv5eAB5/C7B6cD4hA/DAmYqw+IR5GRoKQ4bAlClm9exqYpYe/P77G/Tv/7zTYQghhEiHJCmz2Zo1a3jqqWOcOJH88T/mQNCaVExwDPgm6bPD3hwN317TQr51a1izxnQv6NABevaES5eci88Ky5cf4fvvv+LMmTNOhyKEECKdkaTMZgcOHOCbb/ZzKIV9X8OGQZtPUzFBGKbgvAuKzq+cCu9+8O/3ihSBv/6C55+HkSOhTp203S+zb1/Yvj0v0dGB7nElhBAivZOkzGatWrUiJuYDbrkl+ePh4XBLnlRMEA3kwPF6ZQC0hnsG//ftkBB47z1TZPbgQaheHUaNgrRYXSIiAkJCNpKYOJXz56V2mRBCCOtIUmYzpRQ32lMWHw/LvgQ2p2KSEsAenGm5dK3M8PV6SEzh2ZEWLWDtWqhVCx5/3JTOWLcuoBFaIiEBqlfvyLPPPut0KEIIIdIRScoC4PnnJzFmTPLHgoNhw0JStyesBBBPSg94BtT+TVCqQsqV/QsWhD/+gLFjYft2s2rWrx+cc7jWmi+CgqBDhws0aVLA6VCEEEKkI5KUBcD8+ZvYuDHl47t2QESTVExQFFPd3w3VTTLDqVg4doMkUyl4+GGzt6xbN1M6o1w5+PlnSCsVWl54Adq1W+Z0GEIIIdIRScoCYOnSVQwZEp7i8SAPlMuVignCgK6YFTOnZYNzD0O+lLtL/U/27Gbz/5IlkDMntG9vbnHudLp1lJeuXJnO6NGDiImJcToUIYQQ6YAkZQGT8r6yFSvg2CjgdCqnuILz+8qABA0jl5HiE6fXq1MHVq40K2YLF5oSGq+95v5bmlu2QPfur/HTTz85HYoQQoh0QJKyAJg7dy6tW59JMclQChJigIupmOQ48B6wJRVjWEXDu49A9x7eXxIcDH36wN9/w913w6BBUKyYeWrTrQ85VqkCK1cG8dhjdzodihBCiHRAkrIAuHDhAnv2JHA8hcbjNWrA+NlAavaN5wDCge2pGMMqCuLqQt22vl9aoAD8+CMsXWqe0uzf3yRn778PFy5YH2pqVa+egFIfOR2GEEKIdECSsgC4++67Wb/+DYoXT/mc6vlBpWYSD1ASU0TWDfW/KsGKCP8vr10bZsww+82qVzcb64sVg8GD4WJqVhRtMH785zRoUJeEhASnQxFCCJGGSVIWMDfe+f7DWAj9mtQ9QVka06D8YCrGsNC0NfDae6YWm7/q1IHffoNFi8ztwueeM8nZkCGm8bkbRETEkinTcU6dOuV0KEIIIdIwScoC5IknxvLiiykfz5wZcuUhdRv1S2CW29xwCxNI3AuD+sNifxuuX6NePZg9GxYsgEqV4Jln4K23Uj+uFdq1g99+O0OuXG7oDi+EECKtkqQsQOLiQomLC07xeKdO0GcIEJqKSSKANkClVIxhpdKQ9wWof5t1Q9avD7//bpqdjxgBly9bN3bqnODUqU/YtcsN/a6EEEKkRZKUBcgXX4xm8ODyNzynRn5SXwC2CpCammdWCoIjETDThpW7Pn3gxAn49lvrx/ZHYiLUqvUKTzzxf06HIoQQIo2SpCygbryvbMSLwHepnCIR00dzTyrHsUoiPPU0fPCBtcM2bAi33AJDh7qjC4DHAx99FM/77zdyOhQhhBBplCRlAbJu3TqqVl3AkiUpn1OnFuQsk8qJFPAbsDSV41jFA3v2wU6L+3IqBb17w4YN8Oef1o7tr7vvhltu+d7pMIQQQqRRkpQFSPbs2cmXLw/qBnUv+vaFZo+mciIFlAJ2YZqUu8EDkO8+64d98EHTnmnoUOvH9tfRo+vo3fsedu/e7XQoQggh0hhJygKkUKFCzJgxkjp1bnxetTxAaks9lEoaY18qx7GKgi/XwPET1g4bHg49e8LUqe7plxkfD6NHT2XhwoVOhyKEECKNkaQsoG58bzIhAV5vDfyVymmKAUG4pjQGwP75kD8/7LM4Ufy//zMtmj75xNpx/VWgABw6lMhDD5V1OhQhhBBpjCRlAfTcc29Tt27KZTGCgqBPb/AUS+VEYUBR4Ggqx7FSUSjawqxuWSl/frj/fhgzBs6etXZsf0VHA7zHBTf2hRJCCOFakpQFUNmyZalb98b1Kga+CuXqWjBZe+AhC8axSlbYXQ1iLU7KwJTHOHcOvvrK+rH9NXToJIoXLyKJmRBCCK9JUhZA3bp1Y8iQpjc9r0IEcCWVk4WTymaa1ktIhNe+gdWrrR23Rg1T8f/jj80tYDeoU0fz8MMFiItLTYsGIYQQGYlfSZlS6o6kz/mUUkHWhpTelSbxBg3DN2yAH7sC2yyYai4wyYJxrJIIY16EN9+0fug+fWDXLpg+3fqx/VGnDnzwwVayZnVNywEhhBAu5+9KWTOlVEHgc+AjC+NJ12JiYsiZ890bbkovWxZ6DwQKWDBhLLCJ1D/NaZUgSHwAOrxi/dD33guFCrmrPAZcYfXqF1mwYIHTgQghhEgD/E3KsgIvAM8DshTgpejoaDp2bEmFCimfExIC774CwTktmLAUkAC4qWRWAfhmq/XDBgfDU0/BvHmwfr314/tDa3j44XG8+OLzTocihBAiDfA3KXsdmKK13opp7CO8oJTik0/G0KTJjTd7JVyBIjGk/jtbBAjBVaUxAKb/Dq3uBau3Wz32GEREwLBh1o7rL6Xgu+8SmTatudOhCCGESAP8Ssq01ge01r8rpfporftbHVT6FsH58ze+N/nTT7BzCHAylVMFAyUwSZkL+kNelXgJFi4Bq4veZ88OXbrAhAlw/Li1Y/urUiXImnUk7rmHLIQQwq1S+/RlP0uiyEAGDx5MliwHuHQp5XPuugue+AjIYsGElTA1a93ScgmgFIQ/A8VKWD90r15w5QqMHGn92P7as+cQjRtXZulStzQkFUII4UapTcpcVnTB/W677TbefrsO8TdIkvLlg0fvxxSBTa0KQAvMbUy38MDRizBpM1y8aO3Q5cqZpHb4cIh1yeJUzpxw/Pgujh074nQoQgghXCy1SZmLboqlDbVr16Z//w5ERd34vPAzEGxFWQwwe9PcVN0fIAG6NYdXX7V+6D594MgR+PFH68f2R+bMsG5dLK1by38uQgghUnbTpEwpdU4pdTaZj3NA/gDEmO5cvFiYEzdpzj16JCT+jHl6MrXmY4qX3OCWacAFwYVSUKSi9UM3bQplypjyGNoleZBSoPW7bNy40elQhBBCuJQ3K2XVgBZa6yzXfgDNAem67IfSpf+P55678Tl9+kD7j7Gm50JxzJrmLgvGslJD2J3P+mE9HujdG1atgkWLrB/fX2+9tZxq1apy8OBBp0MRQgjhQt78yv8IOJPM+2eQwrF+eeutd+jcOeXG5ABFi0KTmliza68gEIE1XQIs9vUK+HGi9eM+/DBkywZDhlg/tr86d4YvvihHrlw37n8qhBAiY/ImKcujtd5w/ZtJ7xW1PKIMoEuXR2ncuNRNz7uwEfjbggk9QElgB66rKhezBB5oZ9pLWSkyEh5/HCZPNu2X3KBoUejSZQOhoVb8UIUQQqQ33iRlWW9wLMKiODKU2NhYtmzJx+Wb9EKY+CWoxRZNWgq4ABy2aDyr3AJln4OKNuwte+opCAoyjcrdQmsYP74nX3/9tdOhCCGEcBlvkrKVSqnu17+plHoMWGV9SOnfrFmzKF9+LmvX3vi877+D6laV5i0FdAbyWDSeVSLg70hYY0O1iAIF4IEH4Msv4UxyN+AdoBRMmLCUCRO+dDoUIYQQLuNNUtYHeFQp9adS6sOkj7+AbkBvW6NLp2rXrs348d0ocZPiqfnzQ63CFk0agbmFeeOtbM5IgB7PwrffWj90375w/rxJzNzi2281s2dXcToMIYQQLnPTpExrfVRrXQ94DdiT9PGa1rqu1lqqYfohd+7cdO78MDfb7335MhyZCey0aOIzwO9AjEXjWcUDa/6CPxdYP3T16nDbbeYW5o0K9gZS9uyg1BhiYw+RkGBFzRMhhBDpgdcFF7TW87TWnyR9zLUzqIxgz54IVq++8TmhoTDnG0wabIV4YCHWPDxgJQWJj8Atj9ozfN++sHev2fTvFjt3XqR48bJMmjTJ6VCEEEK4hBVVsIQfnnrqNbp2vfG33+OBgwcgUzOLJs0B5Aa2WDSelULg85Wmb6XVWreG4sXhIxcVcClaFBo3jiV//uxOhyKEEMIlJClzyIABAxg5svRNz4vKDFXyWjhxWWAf5klMl9k0F3LlgePHrR03KMgUk128GJYts3ZsfwUFwdixV7j11k1OhyKEEMIlJClzSK1atahdu+pNz9u5E05/BxyzaOJymOr+Wy0az0p5IFcNe1bLHn0UsmRx12oZwNmzH/D99xOcDkMIIYQLSFLmkEuXLjF7dhD79t34PI8HDq7Gus35eTG3MC9aNJ6VcsOBhhCWzfqho6Kge3f4+Wdu+j0PpFGjDtCxY2e2bXNhuwUhhBABJUmZQ06ePMldd33D9Ok3Pq9oUViyBbj5nU7vKOD/gPoWjWex2AT4YBpssWHf29NPm+Ktn35q/dj+6t4dVqwoTunSN+/wIIQQIn2TpMwhBQoUYN68z3nggRufpxSUzQVRoRZOfrWfpktKRPxLIgzpAS++ZP3QRYpA27YwapSpXeYG0dFQo8YuYKrToQghhHCYJGUOUUrRsOGDZPfi4bsZ00F9BcRZGMBYYLKF41nFAwn3wH397Bm+Xz9T3d9tXY4GDnyCvn37Oh2GEEIIB0lS5qCNG/cyYUL0Tc9TKmmlzMp9YFmB7bhztaw4fG/Tvq86dczHsGGQ6KLm7DExBzl1ajNaa6dDEUII4RBJyhz0ww8/0KXLGWJjb3xey5bwwXfAzfM375UDrgC7LRzTQjMXQdsH4dIl68fu2xd27IBp06wf219Dh8LYsWYFVQghRMYkSZmDnnzySXbt6kBIyM3PrZHf4smLAaG4r7r/VRdg+lTYuNH6oe+7DwoXhiFDrB/bXyYXm83+/TOIiYlxOBohhBBOkKTMQXnz5qVw4Wp4szgyfhgEj7Vw8hCgFCYpc9FtvP8pBmHPQfnK1g8dHGyexPzrL27a6iqQjhyBkiXvZtiwYU6HIoQQwgGSlDlIa82YMfuZP//m5xYuDPmKYm0CVQu40+IxraLgrIbx6+0pJvvYY6Z2mZtWy/LmhU8+0TzySAOnQxFCCOEAScocpJSif/8JTPCioPtjj8FTb2DtT6wIUAUItnBMiz33qKnGb7WsWU2NsO+/d1cx2R49NEWKSIV/IYTIiCQpc9jatasZPty7H0ODIlhbFgPgLLAK03rJhc7ng5xl7Bm7d2/z2W13C7dtG8dTTz3CFTuWCIUQQriWJGUOy5+/CMHBxbw695P+4Pna4gC2Y+qWHrV4XKvcCgfK2jN04cLQoYMpJuumvfV79sTx9dffsnbtWqdDEUIIEUCSlDls48aNvPqqdxXmW7aAElZvNyqDqfBvQ1sjq0zZAl//xE1Lh/jjmWfM937UKOvH9tedd8L+/SHUrl3S6VCEEEIEkCRlDtu2bRvvvLOL3V7UC3vwQej8uMUBZAYK4d7SGEDiLnj0fpg82fqxq1aFxo3NLUw7kj5/KAXZsl0EPuHChQtOhyOEECJAJClzWKtWrbh4cSiVKnl3fu08QIzFQZTD3L48ZfG4VikGUQ9D87vtGf655+DQIfjuO3vG99eTT75Lgwb1pcq/EEJkEJKUOSw0NJTQ0PJen/9aN1C/WBxEWcwtzAMWj2sVD5wrDj/btJrXtClUqgSDB4Ob8p8GDa5w3335iI93Yy8sIYQQVpOkzAU+/XQJw4d7d+5LL0LZey0OIBvwHHCLxeNa7LWP4Y03rB9XKXj2WdM9YNYs68f31wMPwMsvryckxEWZohBCCNtIUuYC06Yt4rffvPtRtGoFbVraEEQmG8a02N6N8NMUexqJd+gA+fPDBx9YP3ZqaH2QP/54ifXr1zsdihBCCJtJUuYCM2bMYOrUil6fXzwBOGxxEJeBcYCL2g79x11QoT94bPhbGxpq6pbNneuu1kuXLsH993/Ehx8OdjoUIYQQNpOkzAU8Hg9Q2uvzP34O+N3iIMKA08Ami8e1UghM3AJ7TkBCgvXDP/64ab304YfWj+2vTJlg9uxERo5s7HQoQgghbCZJmQvs3r2bRx7Zire1Qkd+DuW7WByEwjyFuRu4ZPHYFoo7ChVLw6+/Wj92dLRpvfTDD7B3r/Xj+6t6dQgPH4zWbmxSKoQQwioBTcqUUs2UUluVUjuUUv1vcF5bpZRWStUIZHxO8Xg8/PbbPg4e9O78evXgzpo2BFIO05zcxTXLyA6JJSFPXnuG79PHbPx3W+ulDRs2UqNGKTZv3ux0KEIIIWwSsKRMKRUEDAeaA+WBjkqp/9SCUEpFAb2BZYGKzWmFCxfmyJHFtPRyA7/WELETs6plpYJAVmCDxeNaKQgutYS9me0ZvlAh89TjF1+4q/VSvnwAhzlx4oTToQghhLBJIFfKagE7tNa7tNaxwPdAm2TOewN4D7P1PENQSmGKhUV6eT78OBTr01YF1AdKWTyuDT6cCzNm2DP2s8+a1ksjR9ozvj9y5oRVqy7RoME5p0MRQghhk0AmZQWA/de8PpD03v8opaoBhbTW0280kFKqh1JqpVJq5fHjx62P1AG//jqNBg2CvG7189sMKGN1yyWAGkBdG8a12KoJcO99cOaM9WNXqQJNmphbmFeuWD9+asTHv8HixYudDkMIIYQNXLPRXynlAYYAz9zsXK31KK11Da11jVy5ctkfXACYVjqZ8fbuVKlS0NCuftVXgG02jW2VW+H2t8zmfDs8+ywcPuy+1ktvvLGM229vwL59+5wORQghhMUCmZQdxLS+vqpg0ntXRQEVgT+VUnuAOsCvGWWzf5s2bZg//y3y5/fu/Ph4ODsXe5KnFcC3mBIZbpUV5l6APTH2DN+0Kdxyiykma0exWn/16AETJ5alUKFCNz9ZCCFEmhLIpGwFUEopVUwpFQp0AP5X2EBrfUZrnVNrXVRrXRRYCrTWWq8MYIwOq+b1mUFB8NdPwA4bwqiQ9HmjDWNbKCEW7ukCEyZYP7ZS8PzzsHkzTJ1q/fj+KlAAWrfehFILnQ5FCCGExQKWlGmt44GngFnAFuBHrfUmpdTrSqnWgYrDzV54YTwtW3r3I1EKtm6Bog/aEEg2zJqmm5/CBAiGDatho00lPB54AIoVg7ffdlejcoARI3rw9ttvOx2GEEIICwV0T5nWeobWurTWuoTW+q2k9wZorf9TClRr3TBjrZJBvnwFKFYsp9fnZ84MDYrYFEwl4Bhw1KbxraAgsSvkaG7P8MHBZrVs+XKYN8+eOfy1YsXfLFgwLWkvohBCiPRApfX/qdeoUUOvXJmecrcngM+8OvP8eWjeGRaGYSq/Wek88CFwJ1DP4rEtVigLLOkABfJZP/bly2a1rEIF+N3q1lapcOUKhIXdzTU7AIQQQqQBSqlVWutk98u75ulLcVU1r/s6RkbCwU2ADWUhyAz0Ik2Ux9i/FAoXgo027IELD4d+/eCPP2DFCuvH91dYGMBUzp5dyNmzZ50ORwghhAUkKXMRrTWVK79P377ena8U7NwGee+0KaBsmIKyblcUct9pCqzaoWdPyJoV3nnHnvH9dfo0FCvWmMGDBzsdihBCCAtIUuYiSinatGlP7dpBPlwDtxW2MaipwB82jm+FTHCkDmy3qdBrVBQ8/TT88ot5GtMtsmWD/v1jadOmktOhCCGEsIDsKXOl6sBqr848ehTq3QW7SvNPKQsr/QjsBfoB3ueKjrg9BJ4qDe3aWT/2iRNQpIgZe+xY68dPnU7AN04HIYQQwguypyyNuXy5MpcueXduzpyQIzP23WasCFwA9tg0voX+Ggf9nrOn2GvOnKZw64QJsGeP9eOnxrFj3/Hqq09y/vx5p0MRQgiRCpKUucyOHTvInHksEyd6d35QECxbANm8rzvrm1JAGO6vWQbQHJp9AB6b/lY/84wZ221buHbvTuTttz9j7ty5TocihBAiFSQpc5miRYvSv38XKvmwTUgpuLUQYEc7oBCgHKbcb5wN41spK3y7DU57ucroq4IF4eGH4csvzW1jt6hdG/buVbRuXdbpUIQQQqSCJGUuExwczJtvfkblysFeX7NtG/z5FLDVpqCqJn24PSkDLpyEGvVh+nR7xn/+eVMjbOhQe8b3V8GCicBrXLx40elQhBBC+EmSMhdKSAjm779LeL03qmhRuP1OTG0xOxQBmgGZbBrfSpGw/xSctWl7VenS0L49jBgBMTH2zOGvUaO+o1ixwsS4LTAhhBBekaTMhcaOHUu5clvZ4WWz8dBQmPwdRBa3MahEYDdgU9kJywRB3CMQb+OdvP794exZk5i5Sa1amnbtshIfH+90KEIIIfwgSZkLNWnShK+/7uxTMdRgD1TPin23GA8CYzF7y9KAIUtgg00PJ1StCs2bm1uYbrpbWKUKDB++i5w5jzgdihBCCD9IUuZChQsXpkuXJ8ie3ftrliyB+U8Cu2wKqiCQFbChlZEd1v4EVavBvn32jP/ii3D8uNn07y6azZv78uOPPzodiBBCCB9JUuZSR47kYdEi73881apBz+eAXDYFpDA1y3Zi6pa5XWUo/Sjks6FJOcBtt0H9+vDBBxAba88c/nr99d/p0+dJYt0WmBBCiBuSpMylBg16n1atwNuGC2FhMOI9yGdny6VKgAY22TiHVaJhSwFYbWPpipdegv37Yfx4++bwx+DBsH59dUJDQ50ORQghhA8kKXOpp59+munT7/Q6KQNTyb56AnDWpqDyALmB7TaNb4P/eweGD7dn7GbNoGZNePNNiHNRuZCCBSFnzlnAKhISEpwORwghhJckKXOpChUqUK9ec5+q0x85AtNeANbbFhZ0BDrYOL7F1vwJX33j/YqjL5SCQYNM2yW39cNMSIDWrVvwwgsvOB2KEEIIL0lS5mILFgSxYIH35xcoAD9NhpBatoUE2TCNydNKH/vWUOZZk0DZoXlzqFUL3nrLXXvLgoKgRIljFCwoK2VCCJFWSFLmYr17j+aNN3y7pl0baFDGnnj+ZxPwGWmiwj/h8MMm+PuIqcRvNTevln30EfTps9npMIQQQnhJkjIXGz9+At98U8Kna2JjIec2wKZSEABEAMdIMzXLEs5B9QowapQ94zdrZvpPum21DEDr2fz++1AOHTrkdChCCCFuQpIyF6tQoQK5c9f26ZqgIPhtBGBT4VQAimJuY662cQ4rZYYrZaBQaXuGv7patncvfP21PXP469AhaN68H5988onToQghhLgJScpc7MKFC4waFcuaNd5fExQE69dC8U62hWX+1lQF9gAnbZzHQgl3wSKb9pUB3HWXO1fLChSA2bM1gwbd4XQoQgghbkKSMhfzeDw88cQkfv3Vt+sKF4ZWNq0K/U8VTEFZHxJGp322CD74yDyZaDWl4LXXTAeBr76yfvzUaNQIwsLeRNvxCKoQQgjLSFLmYhEREezZs4EBA3y7Tms4NhVYaUtYRhagMVDSxjksdmEzPN8Pn55o9UXTplCnjlkts+OhgtTYsGE+1aqV5O+//3Y6FCGEECmQpMzlChYsj1K+bfZXCg5uguBjNgV1VX3M/rK0ojxk6Q0169kz/NXVsv373bdaljcveDyHOXHC7r8UQggh/CVJmctt27aN554L5qiP7YLmzYUWfe2J6V9OkHY2/HvgbDYYucqeYrIAd94Jdeu6b7UsVy5YufIS9evvdToUIYQQKZCkzOWOHz/OJ5/sxNe7TkFB0CIQtxbXAVOxr7WTDQa9BS1b2TP21dWyAwdgzBh75vCXUhAf/yqTJ/8k+8uEEMKFJClzudq1a3Pu3FRuv933axePBCZbHdF1qmKq+6+1eR4LnUuEk4lw+bI94zdpArfeCm+/7a7VMoAJE/Zy7733M3/+fKdDEUIIcR1JylwuODiYED/7JhXOCzlzWhzQ9bIDxTC3MBNtnssqteFkcwgJtWf4q3XLDhyAL7+0Zw5/PfggTJuWhQYNqjgdihBCiOtIUpYGTJ48n0cfzezzdW+8Ad372xDQ9aoBMcDuAMxlkZ2nYcTvsNummBs3hvr1zWqZXSty/ggJgZYtz6LUYLmFKYQQLiNJWRqwd+9eli5VnD/v+7UtSgF230IrC0SRZgrJAhAPfdvCs8/ZM/zV1bKDB2H0aHvmSI2JEz+gTp3qXHHb/VUhhMjAJClLA3r16sWWLS+Q2ffFMsa+AR67k4IQoA/g311WZwRDQmu4vYd9U9xxBzRoYJ7EvHDBvnn8ER19hYiIo5w4ccLpUIQQQiSRpCwNUEoB1f26tmVLqNwC+/d7BSV9dtGtupsqDZ9uhXibvjdKwTvvwJEjMHSoPXP4q0kTmDfvKAUKXHQ6FCGEEEkkKUsj+vefTl8/6o7dcw/060dgftK/Al9insZMI7bvg1tbwKJF9oxfrx60aQPvvQduW5RSKoEzZ55nzpw5TocihBACScrSjPPnNZcvR/l1baNCoA5bHFByCgDHgQMBmMsqobB6FWy2sfvQ22+b25dvvWXfHP565pnJ3HtvG86cOeN0KEIIkeFJUpZGfPrpp3z22QN+XfvRW5gVrFhLQ/qvipj9ZWmlwj9AKMT/H8TdYt8U5cvDo4/CiBGwZ4998/jjlVdg/vyKREdHOx2KEEJkeJKUpSlN/GoP1LUrdH4T+3/aYZjEbCNpa29ZELy1AFausa/90qBB4PHgc3N5uxUtCtWqrQBmS4kMIYRwmCRlaUjnzj/Rrp3v15UvD30fBoItD+m/agJxwJoAzGWhQ6ugZjX47Td7xi9YEHr3hm++gXXr7JkjNV57rQtdunRxOgwhhMjQJClLQ6pUqUP16nn9ujZXAkRvwv5N+PmB9vj7sKhzikPm1lDVxrIeL7wA0dHw4ov2zeG/IwQH7yUhIcHpQIQQIsOSpCwNefbZZ3nppU5+XTtrFpz5CbMR324VAJtaGNkmGM5Xg6822zdFtmzw0kswcyb8+ad98/hjwAAYM2YXQUFSTFYIIZwiSVkak5jYyK/SCm3bwqczgNyWh5S8jcDEAM1loXe+h4e6QKJNtcueesrcynzhBfv2r/lDKYAD7Nr1AnPnznU6HCGEyJAkKUtj7rrrQ9q2VT5flzUrdGkM4YHYVwZwAdgA7A/QfBY5dxSmzLCvJ2ZEBLz+OixfDpMm2TNHanTrNoIePbrKbUwhhHCAJGVpTPfuPXn88dJ+XXvqCBSYD5yyNqZkVQHCgSUBmMtKFYBekC2/fVM8/DBUqGBuZcbF2TePPz77LJG//qpEUFDQzU8WQghhKUnK0pj777+fBx980K9rlYKDC4FAFJINw2z23wKcDsB8VvHAuUR4fwHst2mVLyjIFJTdtg3GjLFnDn+VLQsFCkwDFhAba3dhOyGEENeSpCwNOnGihl9lFQoVgqNHIUcN62NKVi1AAcsCNJ+FPuwP9RvYt5J1991w663w2mvua1YO0K3bPXTo4F+xYiGEEP6RpCwNuv/+D+jSxb8fXZbM0L48gelPGQ3cDhQNwFwWi78FyrU3BV/toJTph3n4MAwbZs8cqVGx4imqVEkk0a4nHoQQQvyHSutVvGvUqKFXrlzpdBgBtWjRIoKDX6B2bd+7aGsNDe6CheeBu6yPLT0JD4ZtT0EhGzsQ3XMPzJsHO3ZArlz2zeOf3MB2IIvTgQghRLqhlFqltU72npWslKVBt956K7Vr+1HaH7NCU60cZA9UaQwwT2IuBtLYA32X46HDQPj0U/vmeOcdc/vylVfsm8N/x1i48HEmufExUSGESIckKUujli/Pw9Sp/l07bBj8Xy9r47mhfcBs4O8AzmmRxb/DVxPsqylWrhw8/TR88QWsdmEj9wEDfuCNN16RvphCCBEAkpSlUe+88yN9+vhftuDBigSmNAZAGSAbaa88BsDdEPHY1eKq9hg4EHLmNMmZ23KfceM0CxcWRdn5DRBCCAFIUpZmDR48mOXL7/H7+vEfgWckEIiqBx6gDnCANFdMlnBYdAC+WWlfiYysWc1tzMWL4dtv7ZnDXwULQmTkTBITf+P06bRU20QIIdIeScrSqBIlSpAjR0u/r3/gAWj3DKZkRSBUwdQuWxqg+aykodu90LmzfStZjz4KNWrA88/D+fP2zOEvraF583Z07NhBbmMKIYSNJClLwyZOjOWTT/y7tkoVGNofPIFqHB4G1ADigbRWZUFB7K1QsYN9tzE9Hvj4Yzh0yBSWdROloEOHC3TqFMinQ4QQIuORpCwNmzJlESNH+p9VRQdB+cOYpyMDoTHQkbT5t64cTDgDJy7aN0XduqYF04cfmhIZbvLoo/DQQ9NQ6qjToQghRLqVFn89iiTDhw9nw4aufl+/ezds/BzTCikQrv5tOw1cCdCcFjpzCVr0gM8+s2+Od9+F0FDo18++OfwXw3ff3cs777zjdCBCCJEuSVKWhkVFRaFUE7+vr1AB5i6A8FoWBnUzp4CPgVUBnNMqHli5DOYutm+KfPng1Vdh6lSYOdO+efw1b95Spk37hvj4eKdDEUKIdEeSsjTu44+307u3/9c3qg9tylkXz01lB4oAiwjMk58W0x3gYnN75+jdG0qVgj59wG09wYcOhfnzLxAc7LLAhBAiHZCkLI3bu/c427Zl8fupQK0hZBGwwtKwbqwRZh9bIOe0SjDM2A7fL4a9e+2ZIizMJD/btpnN/26SKRMEBe3lwoWXWbTI9zZfQgghUiZJWRo3ePBgZs78P7+fClQKDq2HsEDu3y4ClMCslqXBvWXEw0MtzEqWXVq0gJYt4bXXTNNyt3n66WG0aNGMM2fOOB2KEEKkG5KUpXGm0rr/+8rA7F3qNsiScLzXCJOQpbVisgDBEN8cajxi7zQffQRXrsCLL9o7jz9efVUzbVohoqOjnA5FCCHSDUnK0oEXX5xFixb+/yhDQ6HTLQS2flhBoB9QMoBzWqksfLIdztm40leqlHkKc+xYWOqyorvFisFtt20BPiMuLs7pcIQQIl2QpCwdyJu3IEWK5EtVtfl10yB4BJBgWVg3F5n0OY3eATt6Guq1hs8/t2+Ol182T2T27AluzH2+/vpZKlQow7lz55wORQgh0jxJytKB3r1789lnT6aq2nyxYlCxJoHf4/UnMAKwsSirbUJg027Yss++KaKi4NNPYd06GDzYvnn8VabMZSpXjuXy5ctOhyKEEGmeJGXpRmMupKIyf7Nm8P04IJNlAXmnHCYRXBLgea2gQHeElSXt64kJcN995uO112DrVvvm8UfduvDTTwfJlWuZ06EIIUSaJ0lZOvHYYyOpXTt1P84yOaFSGIGtH5YHqIBpVB6odk9W8sDi/fDSNzB/vn3TfPopRETAY49Bogt7hx469Dj9+z8jRWWFECIVJClLJ5o3b8Gjj5ZP1YrNxo2w4UVgrVVReakhplF5Wi17peH956Hfs/atmOXLZ3piLlwII0faM0dqLF58iKFDh7F69WqnQxFCiDRLaTvvuwRAjRo19MqVK50OwyW+BB7z+2qtYegn8MYxOB1iXVRemQTsAXoBwQGe2wrH4d6aMKmLfVNoDXfeCcuXw+bNULCgfXP5Sms4dMhDgQLLgBpOhyOEEK6llFqltU72f5QBXSlTSjVTSm1VSu1QSvVP5ng/pdRmpdR6pdQfSqkigYwvrYuNbc38+f5nU0pB317wxB0WBuWtpsCTpM2EDCAX/LIHZmyDkyftmUIpGDUK4uPN05hu+veUUlCgQCLwEEuX/iW3MYUQwg8BS8qUUkHAcKA5UB7oqJQqf91pa4AaWutbgJ+B9wMVX3owbNjX3H57HLt2pW6c+sHg+cuamLyWGQjD1EpLi1X+k7TvAHc1gwSbSosULw5vvgnTp8P339szR2qsWfM3des25LPPPnM6FCGESHMCuVJWC9ihtd6ltY4FvgfaXHuC1nqe1vpqcYSlmBKjwkudOnVi6tSXKFQodeNsWgHBy4HzloTlvQRgJDAnwPNa6GIpyHErqSpPcjO9e0PNmtCrF5w4Yd88/qhaFcaMgW7dijsdihBCpDmBTMoK8O+mOgeS3ktJN2BmcgeUUj2UUiuVUiuPHz9uYYhpW/78+WnVaiAhIdlSNU7PnvDnOszqVSAFAYWA1UBMgOe2SnmYmwM22vjXMigIvvwSYmKgb1/75vHXo49Cpkz/R3z8Sa5cScPLnkIIEWCufPpSKdUZs1v4g+SOa61Haa1raK1r5MqVK7DBudzZs5cZNqwsGzf6P0ZkJNQtAXcUI7AV/gFuAxRgY3kJu8Unwv1vwlNP2TdHpUqmJ+Y338Bvv9k3j78uX95P/fpledGNjTuFEMKlApmUHcSsg1xVMOm9f1FKNQFeBlprreWf2T5KSEjg+edXMn166saJi4Njw4F5loTlvWhMOr4GOBrguS20dQv8MgvO2NhC6uWXoVw5ePxxcFuXo/BwaNz4BHXqOB2JEEKkHYFMylYApZRSxZRSoUAH4NdrT1BKVcXsKmqttT4WwNjSjWzZsrFr1y5eeKFoqsYJCYE76kCeVO5P88vtQDiw3IG5rVIPzj8CF238LywsDEaPhv374aWX7JvHX2+9BfffPw444nQoQgiRJgQsKdNaxwNPAbOALcCPWutNSqnXlVKtk077ALOT6Sel1Fql1K8pDCduoECBgsCDqR5n2DAY1C/18fgsE/Ao0NKBua0SBGfj4cmpMDPZnZHWqFfP3CYdPhz++MO+efx3ku+/b07v3r2dDkQIIVxPisemU6+++n+cO/c5Q4embpwLVyDfE3CuIM7sQLyUNG+YA3Nb4U9Qf8GWLVCmjD1TXLgANWqYW6Xr10POnPbM468BA+D334vz++8byJQp0M1VhRDCXVxTPFYEztmzoZw5kz3V48yZCefGADtSH5PPLgGfAoGumWalehD5KETms2+KyEj47jtTtLZrV3cVlQV49VWYP/8ImTIdcjoUIYRwNUnK0qlhw4bx1VevpHqcVq1g1DgIKW1BUL6KAEpjKtal1conoXC+MHSZDAcP2pcwVakC770HU6eC2+q2hoRAcPBFzp/vxIgRn5LWV+eFEMIukpSlax05ezZ1P+LgYOj+EHSsbFFIvmoChGIq1qXh3+Vzl0DxUqaEhV1694bmzeGZZ0hVSRS7fP/9cp56qherVq1yOhQhhHAlScrSsa+//o3cueGoBaUlSu8H5qZ+HJ9FAo2AXcBmB+a3Sm6Irw55K9o3hVLw1VcQHQ0dOsClS/bN5Y+uXWHZMkWNGoFuFSGEEGmDJGXpWN26denTp7klYx3fCVkOEPhismDqluUBdjswt1U8kNgE+i6By/H23cbMkwfGjoVNm+DZZ+2Zw18eD9SsmQh0YM2a2ZxwW48oIYRwmDx9me5dwGQ0F1I1yuXLMGcvtHaqCfZlTO2ytC4ByiyEznfAK6nf8peiZ56BIUNgyhRo3frm5wfSmTNQpEgQ99zzIF9/Pc7pcIQQIqDk6csMTOtMrFxZnwMHUjdOeDi0Kg2lsgDxloTmYwBJn08Cpx2Y3ypBsPUU7LT5z/D226Y5eNeu5gEDN4mOhgkTEvjgAymPIYQQ15KkLJ07duwYtWvPYeTI1I91+DAcfRPTMNwJccCXwHTS9KZ/7oFZ+eHERfumCAszZTIuXYKHH4YEJ24730DLlpAr10i0/oldu3Y5HY4QQriCJGXpXJ48eZgyZRL9+qW+omi+fPBgB8jtRHkMgBBMw/IdwN8OxWAFBYfPQ/uh8Prr9k1Tpgx8/DHMnQsffGDfPKnx4oudqVmzGkeOSCsmIYSQPWUZRi/gE0tGGr8OHp5syVC+SwA+B2IxTbtCHIrDCr9D9u2wYwtky2bPFFrDAw/AL7/AwoVQu7Y98/hrxw6YPj0vvXptR6nMTocjhBC2kz1lgsmTC/H559aM1aYY5FsBnLVmPJ8EAS2AM8A8B+a3UkO43B1O2vjvIqVg1CgoUADatgW3LUiVLAm9ex9Bqcc5d+6c0+EIIYSjJCnLIH78cQ0jRoRaUorhxAk4MQfYlvqx/FIMqJ70dVpe6A2Gi8HQ9gcYMdI84WqHrFlh8mQ4fRruvReuXLFnntTYtetbypYtxLhx8jSmECLjkqQsgxg+fDhr1ryIUqkfq3hx2L0L7umc+rH81gpoCljw53Ha+tXwZE+wMx+pUsXUL1u6FHr2dF9/zMKFoWnTc9xySzr4gQohhJ8kKcsgsmXLRlDQwyQmWjNegQLwfhMIvoAzq1VXf3fvBWY7ML+VCgJd4VQ5e6dp1w4GDICvv4Zhw+ydy1fBwfDVV4lUqfIycIzY2FinQxJCiICTpCwDWb06hrJlM1vWF/HENtDDgO3WjOeXPcBiYL2DMVihMLw8D75dAsuW2TfNwIHmFuYzz8BsVyaz+/nggxo0aFCfixdtrBkihBAuJElZBlK0aFHy5y/LhdQV9/+fGjWgx+OQpbA14/nlNqAQpnZZjINxWCBRQ5cu0LY92LVQ5PGY26QVKpinMrc7mVCnoHTp/ZQpc4TgYPnfkxAiY5GSGBnSfcAvlo324WJ4do5lw/nuNPAZkA/oQtr+p8ZJKJ4Z1g2CzKH2TbN7N9SsCblymX1m0dH2zeW/p4iN/ZCQkBCUFZshhRDCBaQkhviXS5eeZ+ZM68a7Jz9Ezcb0p3RCNkyZjL3ABodisEoO2BUGD/8Cy5fbtyG/WDH4+WdTJ6xTJ/dV/Ac4c+ZTbrutJB999JHToQghREBIUpYBvf/+LFq2NKslVog5CbFrgFT210yVysADQCUHY7DQL7NModexY+2bo2FDU/F/+nR4+WX75vFXVBSUKbOfEiWOOh2KEEIEhCRlGdATTzzBvHmjKFbMmvGqV4fDB+HWRtaM5xcFlMP8jT6HqfiflhUGWkBkVXun+b//g8cfh/feg2+/tXcuX13d/9amzcfAYq64scCaEEJYSJKyDChXrlzcfnt3TLEva2TLBh82BU5ZNqR/LmH2lzm5x80KCqgF3WbAugNw8KB9U338MTRoAF27mj6Z7nOZ2bObUbp0cbZtc6pisRBC2E+Ssgxs8ODi9O5t3Xh7FmHaazp5GzMCcytzBc51HLDQuVi4tRk0vMO+JzJDQ2HSJNPyqHVre0ty+KtEiXNUrBhDdHSc06EIIYRtJCnLwI4cCWH//jyWbfJu0QL6vQShea0Zz2+NgTzAZJxfubPAhTqQeDtcsfFB6Rw5YM4cyJMHmjeHDS57YKJECZg+/SJ58nQjMfECly5dcjokIYSwnJTEyMASExPxeFYCtS0d942/YMA8nG2BdAIYDUQB3YBwB2OxSMOi8H4lKF8GIiPtmWP3bqhfHxITYcECs3rmNt27F2LHjuLMmjWb0FAb64YIIYQNpCSGSJbH4wFqcexYQw4ftm7cB/JD1HjAyYfmcmKexsxDuvlb/udmqHcb9LLwlvP1ihUzK2ZxcdCkCRxw8lZ0Cho12k+jRoeRfEwIkd6kk19Xwl+XLl2ifPk1lpZEyBoN+cMg0untP8WAdkAoEIczPTqtlAniW8LpOvY2FC9fHmbNgtOn4c474fhx++byx4MPwoAB24C27Nu3gwtWtagQQgiHSVKWwUVERDBs2HBeeKGeZWPmzg2bN8KoXpYNmToXgC+A5U4HYoEK8MtBeHoG/PKLfclZ9eowbRrs3Qt33QVnztgzT2pcuTKNRo0q07nzg06HIoQQlpCkTNCpUyfKlHnf0jE9HuhYEeoex9mG5WCeyMwG/OaCWCwyfALcd595atIut91mxt+4EVq1Arf1Bw8Lg3ffvciAAScBqWEmhEj7JCkTABw5UoKePQuwZ491Y8bFwbn5ELnJujH94gHaYvaX/YSze92sUhroAHvz2TtNs2YwYQIsXgxt29pXlsNf7dtD1aqLgHuZMuVnzp8/73RIQgjhN0nKBADx8fF8++1pliyxbszQUPh9Dvz2CwQ53U86FOgIhAHfAmn9d7cCyppG8EPnwA8/2DdV+/bwxRfw22/Qpg24cQvXnj0zad/+ft5++3WnQxFCCL8FOx2AcIeCBQty8OARoqLaAPMsGzdPHrNA9Vx1ePd7wMkSC9GYxGwG4MIG3P7QQL+XIXwLNG4MOXPaM0/Xrmb/Wo8eZp7p001tM7coWhRmz9bUqbMWuEy6qIEihMhwZKVM/E9UVBTwGStXhhEfb+3YJ6eA5yfA6X1J+TF1y6KBRNL+E5mAbgqXO8FPFjWYT0m3bjBxIqxda/ab7d9v73y+atgQwsPncPny3Tz55OMcOXLE6ZCEEMInkpSJf1m//gq1asUyfLi14775Jvw4BaKzWTuuXxQQD/yIWTVLdDacVAsBnQ+emAHtB0LPnljWpeF699xjymUcPAi33gpbttgzT2qsW/c7Y8eOZtmy35wORQghfCJJmfiXW265hS++GEW3bg0tHTd3bmh7F4xoCbihvEIQ5onMFcA00n5iluTnv+DXhXDZxocRb78d/vrLbPqvX999vTJr14ZduxJp02YgsJGYmBinQxJCCK9IUib+o1u3x8iceTwJCdFcsfiXe55j4PkE2GntuD5TQFPgNmA1MIX0kZg1gsP3QtcZcPYCnD1rzzRVqsCiRZAtG9xxh1k9c5PcuQH2sWFDHYoVK8TEiROdDkkIIW5KkjKRrNjY3DRsmJf+/a0d99Zb4fHHoXota8f1iwLuABoC6zC3MtODYPhxE5RuAg0b2VfGokQJWLgQSpc2dcy++86eeVKjSJEL3HvvBerWTScF6oQQ6ZokZSJZoaGhNGjQlho1rG1WHh4OIz6FWV2hXHbgtKXD+05hkrKmQBVHI7Hc0RJwsjSctfihjWvlzQt//mmS7QcfhKFD7W0B5assWWDMGE3+/C+idT8+/3wEF91WBVcIIZIo7ab/g/qhRo0aeuXKlU6HkY6dBioChywfudMj8MMkSHgCU3XfLTYBZUg3BWNK54ChlSFPJFSrZs8cly+bpOyXX6BzZ/j8c4iMtGcuf61YYfabDR8+lP/7Pxu7ugshxA0opVZprWskd0xWysRNZGPSpK70sqGP5cvPw6DXILdNtbX8cgRT9f97TBPzdGDbCWjTCdo+YN9TmeHh8PPP8PrrpgNAnTqwbZs9c/mrZk1YuhR69hwHHOby5ctOhySEEP8iSZm4qQ0bglmyJBdWd7ApXx5e6Qu/dYLIs8Ala8f3S17gbmAHpvK/y9oK+UVBXAc40Aw+WQGJiXDJhu+1xwOvvmoq/x8+DDVq2Nub0x+1aoFSqzl9ugoVK5Zg9OjRTockhBD/I0mZuKmXX36ZxYu3kzlzaVvGL5sVMn0Lnqm2DO+76sA9wB5gHHDOyWAsEgnxOaHvLLjlQahdB+yqFNG0KaxeDeXKmX6Zzz2H5cWIU8vjOUadOoeoXPlP0s2SqBAizZOkTNxUcHAwISHRXLz4BZ995rF8I3dEBHw5CkZ8BMFu+RtZBWiPaV7ussr1qbUpGPZFw34ba5kVLgzz58MTT8DgwaY10+HD9s3nq+ho+OYbqFlzAlCPYcNeYdq0aU6HJYTI4NzyK1CkAT/+uIsnnkhk0SLrx777bni8KXzZGtiGO2qGlQd6J30GOEa6aMtEKTjTGOp8CR/PgQEDIM6GxaKwMBg+HMaPN5vsq1UziZrbxMWtZPz4d/juu7edDkUIkcFJUia89vDDD7N8+RLq129v2xzFz2H2cq21bQrfZE76fAz4HJhMurnbdTEOen8E73wA23bYN0/nzqbqf1SUKTT76qtYXpQ4NUJCYOHCREaOXAI8xKFD2/j777+dDksIkQFJSQzhh3hWrbqTbdv+pGNH60efOBH25oXn/4AEt/z1TATmA39iHgZ4ANOmKT2Igerl4Id2cGyreXJSKeunOXsWnnrKrJyVKwejR0O9etbPk1r33x/Jn3+GsmfPATJlyuR0OEKIdEZKYgiLBfP221l59dVM2FFVoG1b6Hcr/HA3BP8MnLJ+Dp95MEVmO2JKt43CPKGZHmSFVYehwosmSRozxp5psmSBceNg5ky4cMH0zXz6aTjnsgcphg27wLhxZ8mU6TXgAjt3Ot0TTAiRUUhSJvwybtw3/PnnBsLDm9s2R5EEiD4KBdyyWgamqGwPIArY63AsFruSVA7kk0uw+rDZmG9HXbNmzWDjRrNqNnw4VKxoymi4Rb580KxZAvA+c+YUo1SpUsycOdPpsIQQGYAkZcIvkZGRFCxYHK0nMnBgcX75xfo5atSAfXtgw3twZ3HguPVz+CUH8BjQKOn1Hsyes7TOA1SHdSeh5udQ+Ta45157poqKgo8/Nr0zM2WC5s3hoYfgxAl75vNXrVrHGTRIc8cdw4BdHD16lLS+5UMI4V6SlIlUuXJFMXt2TubMKWDL+JkyQbYIeLcieEYCK2yZxnehmP96NDAT8xDAXNLNQwCJCo5XhRV5YfZO08/SjlvV9erB2rVm8//335uCwuPHmwK3bhAdbZ5ODQubRVxceZo0uYXOnW3YSCmEEMhGf2GBCxcuEBGh8XhaovV8WzaJJyTAhx9CeG14dj7EueSXNgAXgFnAeswqWiugmKMRWa7OcTgyE5Yuhjx57Jlj/Xp47DFTPqNyZXjnHXOr046/T/5ITISvv4a8efPRosVXxMXdwZkzZ8iZ0019woQQbnejjf6SlN1Mnz7mn/LipmLjLrNxwwqKFksku41PJsZchvVbIDESs2LlFpeAk0AC5gnNMGfDsdRlUOehWCkokAXi4yDMhu+9Bo4dg927zcpcdDQULw7RWayfK7UOH45kx45LVK9Rg0wREU6HI4SwQpUqMHSorVPcKCkLtnVmkaEogoBITFZy0bZ5woCgS5ApEixux5k6EUABTFBXE7J40sd/ZeGgw2HXaThwBhL2Q768ULKEtdMoIE9uyJULDh+CvXthzRrImQOKFYPISGvnS40s0RcoWBAiInYCRTlzJoFMmTIREhLidGhCiDRKVsqEpbTWKJUIvMWkSa9Rr14iefNaP8+ZM2az+JSt0H0InDyDqbzvkltdAJwFPgHyYx4KKOpoNNaJB1ZC3rLwThdoXQT27zO3HK12/jwMGwbvv29KZzz8MLz2GhQpYv1cqREXB4ULh1GzZh1+/fVPp8MRQriY3L4UARcTE0ORIgVp397D6NH2FqJq3hJWbYWTnc0GdddIBFZiis6exyRlDUk/yVmSnCvg5EzYucOsZtnh5El491345BPz0EHHjtC7N1Stas98/ti40XyuWPE+YmL6MXToHHr16kX27NmdDUwI4SqSlAlHbN68mTx5MpEjxyBOnhxLRIR5mtJq8fFmH9IxD3SfBCunAdVxz23DOExytgiTnD2NeSAgvbgEbIeqTeHl2+DgXChVwpS5sNr+/fDee2bD/dUCtL17wz33QLBbft7ADz+YxHHDhnFUqPAQly9fJiwsDOWWpxaEEI6RpEw4rlWr6hw6tI4VKxIICrJvnm8mwEOdIbIHXMhv3zx+iQN2YQrQAiwACpK+ntRMhODP4JaaMPNHyB1pkier94LFxMBXX5mVs927oXBhePJJ8/SmWxamdu++unJYiV69srN06TmWLFlOkJ3/AQghXE/aLAnH9e37Pn36DCYoyFRctevfAp07wapVsOND6FkdgtcAbuktHcI/CdkVYDkwFvgK2IjZq5XWeSC+J6yuCIU+gjYjIGcuLC8unDUr9O0L27fD5MlQogS88AIULAiPP24emHb635v/3MrdQO3af9Gy5UaCgp4A1vL555+zYcMGB6MTQriRrJSJAEtkzpyePPvsl/zySyLFi9s3k9ZQuSokZIMjzeHUJcw+L7f8UyQOWA0sBs5gnt5sB1j8RKOjYoAFUKk99G0K5eJh7Sro0gWsriKxfr3pEjBhgimnUbasuYXYsSOUKmXtXKlx9iwUKODhqaea8847P6N1GAcOHKBQoUJOhyaECAC5fSlcZdq0aXzyyRCmTr2V0NDh7NlzmoIF7dkTlJBgntQMywwfzoLXH4aEVkBJ6+fyWyLmtuYa4C4gC7ATk9BUAMIdi8xyYb9D4mqYuByal4OtWyB3UgkMq5w8CT//DN99B/Pnm+S8enWTnD3wgFlNc9qpUyauHDmiWL68DrVrz2HKlB9p3bp90hPMsvdMiPRKkjLhWomJZ7jlljIULXqaadNibZ1r61Zzy+uup2HsPljzN3AZyIO7SmkA/IpZRQvBJGaVgCK45+EFf2ngHJAFckRA8FiI1LB1IwR74MoVCLOw6O6BA/DjjyZBW7nSdAe47TaTnLVs6Y7SGkeOwNix0KNHJrJla8m33+ZlyJCFTJs2g7x21JMRQjhKkjLhWlprpk2bRliYpmnTncTGvs+4cUfo3BnCbV4hat4efp8GUa/A6XhMghaGOxI0DRzEJGYbgVhMUvZo0vErpI+OAUeAS5CjPNxTGqY+CU/2hAGvmsOJieCx6Hbz9u2mv+Z338GWLea90qXhrrugaVNo2BAyZ7ZmrtSYMgVGj/YwZUprPJ72DBu2m127jjF06FBZQRMiHZCkTKQZv/zyI/fd9wAzZ+ajWbPDXLhgVk7suLV57BisWwcN74A/dkPX++BEPMQ9YP1cqRIL7E76ukzS6w8wK3ylk97LjTuSydSIBeZDRDG4oxnUzQHvd4CRn0GHDv9s3E9tXqI1/P03zJ4Ns2bBn3/CpUsQEgK33moStKZNTQ00qxLC1HjuOdi0KYoZM7oDjXn++dnkzJmf559/3unQhBB+kKRMpBlaa5YtW0bt2rVRai1vvdWHoUMXsWtXAlFR9s49ahQEBUO+2+CHjfBtf4ivjLl96CaXgGXANuBQ0ntZgJaYBC0ek6Cl9coLZ4D5kL0+3HUbFDsPo/vDlF+gTh24eNEkTaldUb1yBRYuNEna7Nn/tLrNkgVq1oTatf/5sKsZuy/uu0+RN28eRozoATSmU6fPaNKkKY8+apZRZU+aEO4mSZlIs+bNm8dff/3FoEH3AT/Rt+9wlDrNkCH2znviBNx3H9S9GzLVhD82w6K3IfEO3PV05DlgO+bBgFsxLZ02A79g+nAWSvooCNhQuDegDgMLociDUKc0sAZ+fh/WboSKZUxdsJgYuOUWUlUL78gR+P13WLwYli0zT3XGJ5UrKVLknwStZk2oWBGyZbPgz+Yjrc2K4aVL0LSpom3bAvTp05IrVypTsODLvPnmmzz++BMkJiayZ88eihYtiscNy35CCEnKRPrRq1cv4CQff1wa+IsOHRZw223xPPmkvfNu3Qrdu0PnvnA6F/w6D5aOgMQ2QD7MHi+NO56UPAysBfZj9mwlJr3fB8gK7AVOY2555gRCAx2gRQ4DW0A1grK5IGgubJ4Cf2yBGkVg9jTYvBleftkkMFcTGV9dugSrV5sE7erH3r3/HM+bFypUgPLl//moUAFyONC14dQpePVVaNculEaNqrN1awnKlv2GsWMH8fDDz3DkyHm++OILOnXqRHE769EIIVIkSZlIlxISEmjZsgWNG5fnuefKEh+/iJIlv6N//3h69jSbxHfvNqsbVu9JW7oUBg6C/oPhQAJ89w3MHAz5XobDIZiE4RimSXqItXP7JBZzi/MAUA9To+3qk51XZcPsT7s/6fgxzO3PrDgbu6/OYGJPqkkWOQcSt0PXr6BsTpj2Ppw8ACuWmuOzZpkkrWlT36c6csQkaps3w6ZN5vPmzaaB+lW5c0PJklC0qCkkW7ToPx+FC0NoAJLhEydM4d4774SiRRV//JGTJk2OM39+W267rQl//nmZHj0+YdKkKVSsWJHdu3ezbNkyWrRoQZYsWewPUIgMSJIykSGcPn2aZ599lnbt7qJ58yj27p1L0aKDGTmyBD16nOXYseO8+y507WpuO8XHm8TNil+OmzbB1KlmU/aFeOj1HIwfAT0nwa5zsHoinFgOid0wic9+4AJQNvVz+ywBs1J2LOnjOCZ565R0/BtgR9LXkZjkrADQIum9qw8dZE46HoF7HzJI4J+9dauBcxB1J5TOAXuGQngQPDMKCkfDh09Cvtww8Qdz+hdfmD1krVub13v3mluVKeUqWpvenNcmart2wZ495v2EhH/OVQoKFDDJWb58ZrUtb95/f503r0nsQixOjC9cMH/nQ0LMqt8HH8CwYZEUKFCU0aM9dO++gd27n6do0aqMG7eZ118fz6JFS8iTJy9LlixhwYIFPP3000RERHDy5EliY2PJmzev7GMTwkuuScqUUs2AYZj/TY7WWr973fEwYBymnfRJ4AGt9Z4bjSlJmUjJ6dOn+eWXX7j99tspUaIEy5fPpWHDFvz6az+aNAnnr78W0bDhbObOLUKjRudYt+4UgwfDwIFmhePIEbOfqF49UyohNtYkcd5sLI+LM7/ESyYVqR07FubMgaGjYG8MPPckrF4Ij46HfWdh0acQswcyPQkxl0EvBy4CDZMG3Jf0uXDS54uY/4rsKItxCDiBKV579SMSaJt0fAQmmbvKg3nA4OpTq79hkqHwpI8IzG3Sq7GfwqzAhSZ9dmqr00VMMpo16fUCE094fZOk7X8DcpaCNi+ZHp7DHoSK1eGd4eb1/z0ITe6Aqw9BDhoEtWpBi6Tkdc4cKF7ctICKj4c1a+D0aTh0yCRqe/bAvn3m79mRI+ZYcrJmNf08s2c3SeHVr69+ZM0KUVEpf/hS9+3iRbO6XLas2Zc3Zw6MGQPjxgUTEpKHt95K5JVXDhMb24WQkPwMHLiaN96YTWzsHIKD8zBkyE/8/PMsFi1aglKKyZMns27dOgYOHAjAmjVrOHXqFI0bNwbgxIkTAOTMmdP7IIVI41yRlCmlgjDPi92JuZmyAuiotd58zTlPALdorXsqpToA92qtb1igQJIy4YvExES01gQFBbFz504mTJhA165dKViwIHPmzODxx3syY8ZgypYNZ8KEKXTuPIYtW7pStmww48atpUuX5ezYUZkSJa4wceIR3n33LDNnJpIzp7kdNmECDB9ufhkuWmTKLTz3nFmZWLPGPNn38MOm1c7SpWbFpFkzGDfOJID33w9Vq8NDj8D2HfD0mxCWG17vDudioM0g0FEw8QWTJJZ9Gi6GwI7PITEYglrBZQ8kTMNs7L8NU3B2ISaxqoxJglYBUZiSGgBbkl5frXa/K+n11Ur7B/ln1ewEJklMwDzpeR6TYNXCJGCjMCtxlzH77MDcxr0vKZa3MC2mrgoBqmG6GSjgS0zCGZJ0fjBmRbFi0nXzk45f+1EI85BDLOb/Mp7rPnImxR6LSSg9SXOppK+jkmKPS/rzqKTYE5P+bOFJ40YAeZOOTwZVBLLVh6xRsPdlKNgQqneCrGHw1f1QvyO06QmZQ6FnbXjwaejeD8I0tK4NffvD4z0h/jLcdzd07my6D+zaZVpGlS4N0dEmaVu1CjJlMk+Lnjhhkjhv/vcdEmJaWoWEmH9cREWZv49K/bPyFxxs/sERHW2ax3s85nXWrGbOoCCz0nf1CeiICDP38eOmQG/Tpma8P/80t3UHDgzB44lixIh45s27xG+/1SEoKJqBAzexdOkx/vqrO0FBUfTrN43Vq/exePFHBAVloW/fEezZc4SpU7/G48nEoEEfcerUWUaM+AyPx8NHH33ElStXePnll1FKMXr0aJRS9OjRA6UUP/30E+Hh4bRp0walFNOnTydz5sw0amR67s6fP5+oqCiqVq0KwIoVK8iSJQtlypimtJs2bSIqKorChc2/IHbu3EnmzJnJk/TY7YEDB4iMjCRb0tMdx48fJyIigsxJBe5iYmIIDw8nPDwcrTUXL14kNDSUkJAQtNbExsYSHBxMUFAQWmsSExNRSslDGBnMjZIytNYB+QDqArOuef0i8OJ158wC6iZ9HYz537+60bjVq1fXQtjh+PHjesGCBfrSpUtaa63XrVun3377bX3u3DmttdbTpk3TzZo102fOxGitL+gxY4bpYsUK6XPnVmmtV+u3335cA/ry5Yla65/1oEFtNaATEj7XWn+qX3qpsQ4O9mitX9daD9DPPltbZ8oUorXup7XupZ9+upLOmjVUa91N79hxn77//iI6T55wrXVHPXHirbpRo1y6UKEIrfU9+tVXS+sqVbLokiUzaa1b6Ac65NUlSkbo0mUz6e0nb9Xlb4nUufOF6FIVIvQfu6roPAVDdPbcwbrkLeF63NrSOmuuIB2dM0iXrBauP1xcVGeK9ujMOTy6ZK0w/fIfBXRoJqUjs3t08bph+ukZeXRQCDoiq9LF64fprpNzakCHRytd9PZQ3f67bBrQIVHowg1CdOMPo8zrzEoXuTNUV3w0XAPaE4KOLuHRBe8I0YAOCkcXbhaio0t7NKBRZoywHEoDWgWjC7YI0Xgwx6/9UOZ4gZYh/z0GGg9ahaAL3B2c/PEgczx/6xsfz3d3CuMHoVUoOm+roBTnJwido/4Njit0jupBGnXdn0uhg0LN6+K1Q7W65s+vgtDKgw6NMK9L3BKq1TXXezzm+uAQ8zpn3v+Of/U8QIeEJhObfLjiQynl6HGPx+Pq+Kz6Hs+bN8/23y3ASp1CThPIlbJ2QDOt9WNJrx8Camutn7rmnI1J5xxIer0z6ZwT143VA+iR9LIMsNXm8HNiEkThLvJzcR/5mbiT/FzcR34m7hSIn0sRrXWyHX/TZCc9rfUozE2SgFBKrdQpLTUKx8jPxX3kZ+JO8nNxH/mZuJPTP5dA3sg+iNn5cVXBpPeSPUcpFQxEYzb8CyGEEEKka4FMylYApZRSxZRSoUAHTMWka/0KdEn6uh0wVwfq/qoQQgghhIMCdvtSax2vlHoKs5k/CBijtd6klHods+ntV8xzV+OVUjswD853CFR8NxGwW6XCJ/JzcR/5mbiT/FzcR34m7uTozyXNF48VQgghhEgPpDiKEEIIIYQLSFImhBBCCOECkpTdhFKqmVJqq1Jqh1Kqv9PxCFBKjVFKHUuqaydcQClVSCk1Tym1WSm1SSnV2+mYMjqlVLhSarlSal3Sz+Q1p2MShlIqSCm1Rik1zelYhKGU2qOU2qCUWquUcqxNkOwpuwFvWkOJwFNKNcA0wxmnta7odDwClFL5gHxa69VKqShMI6d75L8V5yjTITxSa31eKRWCabbVW2u91OHQMjylVD+gBpBFa93K6XiEScqAGtcXqw80WSm7sVrADq31Lq11LPA90MbhmDI8rfV8zNO5wiW01oe11quTvj6H6aZZwNmoMrakji7nk16GJH3Iv8IdppQqCLQERjsdi3AfScpurACw/5rXB5BfNELckFKqKFAVWOZwKBle0m2ytZg27HO01vIzcd5Q4HlMu3vhHhqYrZRaldTK0RGSlAkhLKOUygxMBPporc86HU9Gp7VO0FpXwXRQqaWUktv9DlJKtQKOaa1XOR2L+I/6WutqQHPgyaRtMgEnSdmNedMaSggBJO1bmghM0FpPcjoe8Q+tdQwwD2jmcCgZ3a1A66T9S98DdyilvnE2JAGgtT6Y9PkY8Atm+1LASVJ2Y960hhIiw0vaVP4lsEVrPcTpeAQopXIppbImfR2BeWDpb0eDyuC01i9qrQtqrYtifp/M1Vp3djisDE8pFZn0gBJKqUigKeDI0/2SlN2A1joeuNoaagvwo9Z6k7NRCaXUd8ASoIxS6oBSqpvTMQluBR7C/Mt/bdJHC6eDyuDyAfOUUusx/8Cco7WWEgxC/FceYKFSah2wHJiutf7NiUCkJIYQQgghhAvISpkQQgghhAtIUiaEEEII4QKSlAkhhBBCuIAkZUIIIYQQLiBJmRBCCCGEC0hSJoQQQgjhApKUCSGEEEK4gCRlQogMSykVoZT6SykVlPQ6Ianw7Ual1E9KqUxJ7xdUSk1RSm1XSu1USg1L6vJxdZzzSZ9DlVLzlVLBzvyJhBBpmSRlQoiMrCswSWudkPT6kta6ita6IhAL9ExqITUJmKy1LgWUBjIDb10/mNY6FvgDeCAg0Qsh0hVJyoQQ6ZJS6k+lVNmkr3MopZLrZdcJmJLCEAuAksAdwGWt9VcASQlcX6Dr1ZW060xOGlcIIXwiSZkQIr0qCWxL+voWYMO1B5NuPxbXWu+5/sKk24/Nk66pAKy69rjW+iywL2mO620EaqYydiFEBiRJmRAi3VFKFQEOaq0Tk966BVh/3Wk5gZjr3otQSq0FVmKSri99nTtpJS1WKRXl67VCiIxNNqMKIdKjyvw7CasO/HDdOZeA8Ovf01pXufYNpdRmoN1172UBCgM7Upg/DLjsW8hCiIxOVsqEEOlRFZISLqVUKaAN192+1FqfBoKUUtcnZtf7A8iklHo4abwg4EPga631xetPVkrlAE5oreNS+4cQQmQskpQJIdKjyoBHKbUOGABsBrokc95soP6NBtJaa+BeoL1Sajtmn9pl4KUULmkETPczbiFEBqbM/2+EECL9SEqeqmmtz93kvGpAX631QxbOPQnor7XedtOThRDiGrJSJoRIV5I22OubJWSYk1YD864Wj7Vg7lBMPTNJyIQQPpOVMiGEEEIIF5CVMiGEEEIIF5CkTAghhBDCBSQpE0IIIYRwAUnKhBBCCCFcQJIyIYQQQggXkKRMCCGEEMIFJCkTQgghhHCB/wdV85mwD2pwQAAAAABJRU5ErkJggg==\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}')" ] } ], "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.7.8" } }, "nbformat": 4, "nbformat_minor": 4 }