{ "cells": [ { "cell_type": "markdown", "id": "3f68b11e-4a92-4998-8d50-c257d5cc01ae", "metadata": {}, "source": [ "# Benchmark example of Rn220 fitting" ] }, { "cell_type": "code", "execution_count": 1, "id": "771c69a1-540e-4690-8dd3-075782460667", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Using Normal as an approximation of Binomial\n", "Using aptext package from https://github.com/XENONnT/applefiles\n" ] } ], "source": [ "import os\n", "import numpy as np\n", "import jax.numpy as jnp\n", "import matplotlib.pyplot as plt\n", "\n", "import appletree as apt\n", "from appletree.utils import get_file_path" ] }, { "cell_type": "code", "execution_count": 2, "id": "6a222673-1eb2-414e-9868-708b2ca64d4b", "metadata": {}, "outputs": [], "source": [ "# set GPU memory limit\n", "apt.set_gpu_memory_usage(0.2)" ] }, { "cell_type": "code", "execution_count": 3, "id": "9d6eb0bb-24c1-4418-83e5-19be539955e9", "metadata": {}, "outputs": [], "source": [ "# Load configuration file\n", "# You should change the configuration file for XENONnT\n", "config = get_file_path(\"rn220.json\")\n", "\n", "# Initialize context\n", "tree = apt.Context(config)" ] }, { "cell_type": "code", "execution_count": 4, "id": "a17146d1-e0a6-4d3b-99e4-5c165b17ab2f", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "========================================\n", "LIKELIHOOD rn220_llh\n", "\n", "----------------------------------------\n", "BINNING\n", "\n", " bins_type: equiprob\n", " bins_on: ['cs1', 'cs2']\n", "\n", "----------------------------------------\n", "DATA\n", "\n", " file_name: data_Rn220.csv\n", " data_rate: 2000.0\n", "\n", "----------------------------------------\n", "MODEL\n", "\n", " COMPONENT 0: rn220_er\n", " type: simulation\n", " rate_par: rn220_er_rate\n", " pars: {'s2_threshold', 'rf0', 'fano', 'rn220_er_rate', 'py3', 'py2', 's1_eff_3f_sigma', 'p_dpe', 'drift_velocity', 'g1', 's1_cut_acc_sigma', 'py4', 'py1', 'elife_sigma', 'rf1', 'gas_gain', 'py0', 's2_cut_acc_sigma', 'w', 'nex_ni_ratio', 'field', 'g2'}\n", "\n", " COMPONENT 1: rn220_ac\n", " type: fixed\n", " file_name: AC_Rn220.pkl\n", " rate_par: rn220_ac_rate\n", " pars: {'rn220_ac_rate'}\n", "\n", "----------------------------------------\n", "\n", "========================================\n" ] } ], "source": [ "tree.print_context_summary(short=True)" ] }, { "cell_type": "code", "execution_count": 5, "id": "3c429d1f-b69f-4db7-90d3-e0608322982d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "With no backend\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 100/100 [02:34<00:00, 1.55s/it]\n" ] } ], "source": [ "result = tree.fitting(nwalkers=200, iteration=100)" ] }, { "cell_type": "code", "execution_count": 6, "id": "4a29843d-9570-4996-8b4e-23fc80584763", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlAAAAGwCAYAAABmTltaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB38klEQVR4nOzdebxkZWHn/89Zaq9bVXdfuu/thWZrNhUUETU6ISDDJBKdjCb+IkYSRoWJgKNiougkcUhwYoxLJMaMJDNmYpwZHcWIIWgwKoKyKFs3Szd0993X2qvO9vz+qNu3+9ILXc3tvtzm+/Z1X3ad89Sp55xazpfnec5zLGOMQURERESOmL3aFRARERFZaxSgRERERNqkACUiIiLSJgUoERERkTYpQImIiIi0SQFKREREpE0KUCIiIiJtcle7AieiKIoYGxujo6MDy7JWuzoiIiJyBIwxlMtlhoaGsO3DtzEpQB0DY2NjDA8Pr3Y1RERE5Cjs3r2b9evXH7aMAtQx0NHRAbTegFwut8q1ERERkSNRKpUYHh5eOo8fjgLUMbC32y6XyylAiYiIrDFHMvxGg8hFRERE2qQAJSIiItImBSgRERGRNilAiYiIiLRJAUpERESkTQpQIiIiIm1SgBIRERFpkwKUiIiISJsUoERERETapAAlIiIi0iYFKBEREZE2KUCJiIiItEkBSkRERKRN7mpXQEREXhw8z6NarWKMwbZtHMfBcRxs2yaKIowxS3+rzbZtXNfFdV1s28ayrOe9TWMMhKu/byvFcl/cbTAKUCeg8XqTkh8uPTbGEAXBsjJNL8THwrL3+wIYCCOz788YzH7bKNZ8Gs0yjl/HqZaxm3VMs4lZKBEBJpUiwiFWm8NgU7VilK0Yc1GTmUaDyIQExsdE0PDiEEZYUUjC84kFPpgIjAEiLAMFHFrVs7EICZwanuMThjFCxyG0HALLBtvCsiKMifCsED80WFGEhQELQsAYl1iQhsjBxya0LCJa2w6wMFjEMLiWB5aPHdqkaYJlE7MsYrZNZFsEJokfpnDxiawqtvEWf/DBslIYWj+yHgGBCXGMhYsBA5EB2zjEHYOLwWATYYiskNCKwAZjWVhW66hHWK03BbCAmBXiEGFjiIUhaVMl4dSwDNgmwIoMfhAj8F2MsTHYhMbBxcLC4OEBPk4U0iAgjGJYro9F61i5tF7bbr0F2FHrs2FZFqHl4FmxpdoYIuJWQMwKwHaxrFZztgFs28LEDZHjgAWBE8NYMbDYb4/2lTetowAYLMvGiQy2BZYb4gQhjgkJAc+ita9BDCeMYQN+6+jSeqttksRZPHzYUevAWWbxAO49kvbeBfv2xiyeHC1YPFrW4jEwhHZEYHtEhETGxg5DUg0f11jY2Iv7sPgClo1jObi42BhsE+IQtrZpORhsjGUv7nVEYIdElgHL4BBhDBgrxCx+dsFtVTMKl76je2teTqQJndYnyWAICfAJW5/7xf06lL1rQtsldBwA/GYCY+x9B8vYWOw7Lq330gbLXrYd61nbBLBMRMxrtj70xuBEAU5YJZMo45omERaBE8O3YwSOC9gY28G1HWLGwjXO4vtgCExr/wCM1TraGHDYr6773t3Wd5E4FjaW3Tqm+2cWe/+aWhAu28byY+YYC2NbhK7T+g5by566dIzNsqXLj7QVQWRbPHvbrc/ivkeu72OH0eIaa/FdtQhofW7234IhwDPB0qO4FWIR7fsc7lePEIsQZ7FoawtxK87iL9PyGluLHzAr5ECtGlhAZCIi02Dvr5S1d18tiMzSI8DCscLFmttwhCHUWqyr9awj6+Bg42JZNlHd50N/9Jkj2t6xogC1xnzzq19kz9PbAKjZLmUroOY1cA04dQ8TGeJ+DHfxm2mAJi5E9rIvih3GCJ3WSdaw+MNk9v3w2hjCxVOcbQeLn32DZ5zWyWbxc50zNTqM3zrlLH73HCvAdgPcKKLDhHQA63BbX+LWN4OECZbqE1kWkd2qh734Y4oFNaqtk8vijsQCSHoOFhFuYEgYh+W/hTa25WBbrY91FO37MQptg2ci9p6yMRbGWFgWrQBhRVhuhAltTGQTGAefeCuYYZZO8NgRdtwHDI6xWyc0a/GHxZjF+ljYOITYrfBqTOs8svdnJqR1orJaJS3s1j4fcCpa/Hlb/OWMsDC2wViAZRFGMUK/sJgqWuUd18dOB1iEGKt1jD3bpenEgDiQACxytN5/i+U/oNGy3/h9P5gOECNYCicAxliY0F363CzVPLIwgYXlt7bhEGJH4dLJx1r2GovLLGtxM/ve71bgYPEEb5FZrE9ERGRFGFo/YEufI0IaVIkIaJ0m925zbzBvpSkTAUvfkP3sfem92WXx4NiRgxUlsBd/LkMH6gWLyNp3Ym/90LfqFOITmQbGgtB2CHGIrBDsOsYKF09Ore9LPHKIW62TW7j4GTBRDCuyW59RywfHYNkOxo4tvZKxLdJ+HYdg8T2ySZkUeZNc+gwGlk2098S5WNO9/1Gx9/VjTZ+E3wRj4cTL2HaAMYufVCtsfe73/i5YEY5nln1eDpXRjA1BJobtxDC2wcfGNxlqwQi+lcQ1IekoJBYGOFHQOvV7AVHk4RERWWHrO4lLmtZ+W4BtWt+V0LUJnMWkv/juR8a0vlxY2FYZY0dYUStm2Y4DVqvuITaO1QpWtuUvffNa713r99IYC5cIzwkJTUCiES0evdZ7FwvqRFaTyIphmTg2No6JFn9NraXt+ZZLZFvEjIVluQRWDIPV+s4bgxWYxc+hTTXdQaUjS+Ta+z57ROSjCq7tYVvB4mcE3ChGwk+26mw7VKw0ez+yWAF2FBDDxw4NicAnYRb/a8JYOMYQUsTGxzEh9mJQax0dGxPYhFGcIFz8dlmtz1hgtc4TERZu5IKVIyDeCsNAsLhnFkHrP8xdC8cy1KIM8Sgibhpgh0v/gWgAK3SxvRi27S5uO8JYrc93PR4ndFq/AZYJCUKPqtXEJ4CoSbYncfAP33GkALXG7Ejmub9jC5FfJworeFWHzkqDyEniO0m6nDTPxEMaltc64RqDa4fETAXPsbEsh6TnkvRCUn4Dx/hgQcyKiLk2nnEhcgmMTdxu/TwFvotrBUT4dGWfwVv8QvVRx40sSlEcHAffdogsB7sZEAY2nuUQOGmq5EjZIVnHxYkgbeI0TGvbNhGWCWli40QhxguIewksy8ExFjELItMkHtax6g0qdoLAilGO2YRemWQUkqFJK0G0ugN8HCzLwrINCauJ51iUslm8RJYG+daBtH3ibpW0WyNuGZJWgtB0EYU+jaYhHvlkTBzXhSAyeMYstlwYYjSx7NbJr+kHRKGNCVsnedvxMfEmYdDAtgxROkaUjGHHLLAiiBnAwbJ9LCvCJUGcCDuoAhExExAP0mTcXhwWfyAsm3iuA8txWs9NJIjScUomBXYBaGJMfbEFD6qBSzXeA24SY7lkw4i+xRaMdNwma1l0xR3cWJx6aPAs8PwGWCGdKYfkfs3ylu1gOTa27SxmueVnTAPY1vIoEoaGRsVvnXf3Fn/Wf3jajo3jxvZ1iyyWiyVdEmkXvxHQrIWtgGtZRCai4tWICJf9V6xtWThxh3jCwYnb2IEBz5BMJkgmk0vlXMvFsRwc28WxbWLJBE4qdkTfuRebyESEUYg5SDpyLAfbWpnurBcCE0WExSJRuXzAOqezE6ejYxVqJWuFAtQaY5opco0OkqaHWAgd/gJOPknMSZC1LYpU6fEdfCZI+jXCeABRHZw4lm0RNSNcvwMTG6fYaROFcayYT2QHLEQ2SdvDTYQYt9UV0mUXaTo14lWD8W3ybpoBt5NcWOXJcAPjTgwnZrUaVgMXYwzhQoNSMo9jItKBx0g9Igx9TOQRWYaGWyZmGjimSWjH8GIdFLt66fV8Oq0suSDA8gNqToV6o0i26GElOsie3Mfp4SzxsIkxDr6Tpm5iOPkhCvE5GtTwghheOSQMICBiPt1J3EpxWrNMhEcQjbZauwILmzSBNUyYThFiMGGNKO6SK7hkrThO08aK0iRsyMbBNJpE5TomkaMMBLU5wtQcHb29JJKt/9J37Rh5t0B+eIR4Prf0vjlOAoxF3RvDmIhkvBfHbp3gQwPjCzUwhiC/AexDfy33ntJc26Jc96mHETF7+cksbdv0xvZt49lBJt0RJ55qrQ/8kGYtIPAO1mR/dGzHJntyAss+MU6yLza2ZWM7L46xLZZt43Z2QmfnaldF1iAFqDVmeGGBjtknqbgVGjGX2vwcbsYiSNjMx3yaVoIwqmEn4lRzncSdJF0VsFJJEnMhmE7CTJl6ppeT413EYzVct0q+kcIJe0nly5QIoG5RDSzCwCJeylIzLl3ddbrx8OwsU/YQp/V18JJMAtMImQsDFmwDUY3N3gN0DqQYH50iavjEqLJuYwE3nsP2yhA6hE4n6e4ePCsFwGhkk16YpVGp4qaTjIVN6tMlesMk2fX90AhJNmcoWRnCjgJOcx7jOFhZQz3cRSXqwo7WUa17ZE/qJD2QYbaaZTj06ciGFFkgJE4hMUg2m2XDQCeOs+/jPz8+RtfQBiKvSLoBUaKfKEjidrfqZ4IA7+mnmV4/Qm2hSGFqO+5Aiv6OdeRSA4d9z6IooFZ/BrBIdb0K23YJI8PTczUiY7Ati+HNaZwjCBzVhSbFus+o53Fqf5Zc8vm1orgxBzfvPHdBERFZRgFqjenKjzGbHycXxjgpNIwVfPo7+6hhYUUxKrUmuYzFoFmH7dfJmgRhUMZ1cxQ39uBHVYJswGCym4nQJh03uKSZnexhMBcx3ujAIkaiWmHBG+PkaoNMpkKxp5tEWKVkzsGELsNd2VavWQMajkXSjnOOmSG2eztjsY1MzZcZGHkZoR/Re/YriNXHwYREqQFMLEPd2049yDG/UGIs8uieeho320e912dXrUzvWJ3h/Hrczf3YiRR23KdhdWM6+omNPU4zvQ4710HDrzLvd7Ew6+GFDbJ9EYlMgXjFYVNvnsD1eaI5Tn/ylRQSLt3ZBGEwx2ixhL04nqA2N00y30lz1xPE3W6mcxtgrE69s0pzbA4A68mn2LFuHbmf/xjXhv7hMzDeOHPVDuaq1UO+X1FYI/DHiSU2t7o96k2giWXBcFcaIkN5rkF5ut7qItsvQ7UGpi/f3nQM4oUYL09lV/iTJSIi7bDMC+F60RNMqVQin89TLBbJ5XLP/YQ23P/AY8TiWTr6u9j+nW/zSHEXb7j0LfTlknzriW2smy0yOJKnnHY4KXcSwd13M9eZY8wvMFsuM9JRYdyuYbuD5FMzBD6UmpvoHZ9lwp4llY3TORswbZXoHH2KnlNOZff6M+h1Y4w+dQ/Zao5sfm/uNsS9eVzTpDuawJ+tE1gdTG9cRz7TQzDnkV6fJ5bKQHYAbId6pUyjMoo/9ziVah5SaWZ8jwFvN9WSQ7JgMVANGJuNkd/Sw6aeHqyFndQrIc1EAtMsk+reSDKTpRiO8+jcevq9iJGhGGF5Dj+ZocNJYFlQCevMhWVG4n3YNkSxOSwbXDuLa2XAsggCn2alTKbQRdlNMBN5+KM+QRYSTpy0k8QbHWU0X2Akk6ajZ4h4No/nT+LYGfYE8QPeo71XtgRhiSiqEo8NHlhm75VijkUmnwDLWjZMyLIOHqD64zFSL5LuFRGR462d87cC1DFwLANUbXaMytw4u594gunpOUrU2HjmRh6bLbG+HJLKd5JIVulNn0Yw7TP1zPcxQ6+m1hxnfXKean8HielxpoM95FND+PYW+ryAubkq5e4UZwyexNM7HyeX6YQn/pXeV7yOuco0ZuZR4mEOp2ML+b7WwEpjgUn1gl8lmKvB8Dk8PTFOszRLN3nCDsh0djNW9gjCkHhjkngyRT16nEb2PFLxkOTko+yMhWyenif9uosZ8Dp4+F8fIXlSJy8d7iCsVyiVXZrrT6E+PU5vdzdWLMETE0/yxEKck02GjRsKFKfGCd0Y+VwnHV0pSkGRWlBlXXY9jcY4frMK4SDgUA9rjFZHSbspitMT5HtbXXBJJ0WX30F2IIObbnWNze4ZZdZyOGXd8m66avUpSu4wccumO35gQ26jOdHaZuLw3XsiIvLC0c75W114a8xsVOKx4g46eywsk2BuYY4oniY3so746DypkfV4/hgTXppZt0Rndydkn2B9/Rk6UjbFicewY5sg/RLy+XX0uTXMdJ1nEjYXbuhifOyn5OohmdL9dG45lWkrIOH00MyfRabjdJyYQzK3r9WlPjNFY/cYQb6H8OnHqOzazYaB9djpBQr9Q5SaFboThoGsQ5Q7j0r9SXzOpTI2x8Yo4t5Np/LK+TKJ7gsxtTw/+853ibosTs1vYrbRQfGZacJ1W5h+6FGGhwZ5dLpC0i0yE/mc43bSM5Illg5J+A6R65LtiTPXmAELRgoj1Ou7icUzZLLrAAiigKniHs4bPIfqwjx29wCpjtaXJKr5RF60FJ4mpmdp+gGnbFq37D2o10cJYwN4oWEwdeBXqFZ7BtfNEo93H6uPgYiIrDIFqDWmGnbQ4/WyMDXOaL6PvpLFSWygqzRNPdmLPWGoJ+PEbZtO7xGirCHldJDrzzE38nI2zpTYZQ+yLjlHf6oHy7LYMzdHuhN+/kST6dIgfadspbj9Lp5ubiLc06A3ERJGZYp75ugZcqjMgD87jtf0MQs1okKOZCLEj7v4QRXTYcj3D2GAyVrEls4YYTEiWtiG8Z+mYjpZ/8zPmShk6S2tIzv9FMmTfZ78xqN09uUZOXkryajBk0/sJL1pC7ueeoyThgcYe/ghhnoyPBbN0D9TwMrXqDVtJufmiKXSOI7D9NO7SLop0m6K+ebjOE4O49YJmAFj2F3Zw3B2PbVntlOemSbX20eNVmsRkcHtTuGVoBqENB2XDZs2LDv+xhjCsMZo2MXp2eXzkPj+PJ43SzI5hOOkj9MnQkREVoO68I6BY9mFN71tBxUrwZ4H7uOhoM4FHVme6U2QnvMJYxlisTjFyhj5KE5U/AmJrq2k/BmqsRSzVki93k0inQLm6Mutx1Q8Hq8liJVnSbgRw84I61KP8vhklkLaUDjzPDo6PCqzHkEjT2e6SOBHJAtd2MVZatks+f7WGJ/HnnycPgzdW04F4KnZKpu60rDQxMrFKVcepmZlyc5VKM4/QWzdhczP7mEgijMzW6Gedenty9I7MMTCZJOp1BS7nq6yuTvHZLFB36ZhGpToWghJO3k6B5J49Rq+59E00FUoLB2nevNp4rE+HHtfkNlZeZrh9DCu7dCsVendsKk1wd6zRMawvdrg9GzqgHW12jM8E3ZySrYDZ3GAUhQ1qdd348YKJOI9K/hui4jI8aQuvBOYe9II5e3bqZw8QuGJp5je2OCszBB2Zx9WkKBjUwe7HqwxMpzAerCAfd6lRI/9DRNbXsr5bood9SyZfCelyceJsZ6fLuzmpR0JYp1dRM2Ak88awXnsZ5SHswyevpV87wDTY/eRzw5gx2bpOu10rGqJqF4nvvUc5sdHKfQthpSfzzJ84auol4tMlBp0JW0s4xOEDfz6NFhNfGNjxh6mgo83s4f+Zx5mLt5N2WuQpEKqdwvjj4ySzifYvnucf7tpA8W0zciWAbIJw/zMFI3YAH0ngWU1qTdmyPR3EvM84tkmAPX6M2Szgzi2A7SWjVXGGenrJOW25syN51wazd0HPcZP1ZpsSsWZm5tcnIF8nwnPpisd0qzPLS5pDRnPZLas/JstIiIvWApQa0y6tkD9iYfYvO4UHrPrJKtVRpvdVBZK5DpyzP0sIF3dSfD4KM3oDAqNccpTe/C7PLY10mzoT/DkYzPEkh65xCyXn7eFxx7dTiGeIzWQI1f5MbviwxSyebrzWaLZMawdPyPWAaQ7KT/5fWJdvTg9BSpTj2BsQ63m8/R8Bbe8m+nRnxHv7MIkob8rTlStE+uDTLqPp3c9RerxSUy1wsArXo3npOhonERq82kkp+9lyzn/nvmJGqmRGLsmKrh9cZ60BqFh4c9bTEcVGqUkXcNZ5sLWDQiKgQULTQpdBZoBNL1JHKcPN0ot3eXBYOjJnU4hWXjO47u74bGh4GB7u8jlzsB1l89EPFttMJhJHuLZIiLyYrFmroe+6aabePnLX05HRwd9fX1cfvnlbN++fVmZRqPB1VdfTXd3N9lslje/+c1MTk4uK7Nr1y4uu+wy0uk0fX19vP/97yd41o12/+Vf/oWXvexlJBIJtmzZwq233nqsd++IzXiGicw69iRCYt1b2BbEcHPdJPMFUgMBuDUcs0CsvJ6eU07GmDphPEuuq4+enhTFZsTGUwpkHMO60CFWDbGnPWws3Fgd6gvUTJ741DjUFyA2Qfz0C3FOOouwO0b+jNfSsf4lpNMbISjQ3X8mYSPL5DNFNqw7nd74EGMP1+iZNtS2N6k9XKX8VI3HvvcDZnf65HN56qedQz22mR0/eIzt9gDj86PMNQd5+IkKP5tp8ODOaR4Id3Pq0BZG1m/kgjNO4/TNJ7FxMM3moZewqXcDw7lhCs0kp64/m6GOIYZzwwxl+hhK97Kp83SGc8NLfyO5kSMKTwt+gA243i7i8Z4DwtO05zOY0O0/RERkDbVA3XXXXVx99dW8/OUvJwgCfu/3fo+LL76YRx99lEwmA8B1113Ht771Lb761a+Sz+e55ppreNOb3sQPf/hDAMIw5LLLLmNgYIAf/ehHjI+P8/a3v51YLMZ//a//FYCdO3dy2WWX8a53vYsvf/nL3Hnnnfz2b/82g4ODXHLJJau2/3vVSlXO7exhMj7PaP1JtiTnyO++n419XRD20JfJEZY6iLvryJ56BlM/u42hLa9mLt5F5/zT+LYhPmGzrvI0uazD7FMPk0nn8WcfoPj4NirxzYynK5DIUtvxIFWrSuT1YMemwJqnPDe2VJfK3AyzkzuIvADGdtF042xLJ5jOGPx8SLPhESSakLIYL3Rxdu9JjM4UKVJhcsfTbLUdBs7qZ/ejPyU3cCYz9RJ9Ax5J18WvZ9iUyUBYY3a2BsDc1DT5rg5mZxsAFGdm6MCiu7ubyBjGSk/hJ06CWvOAe6IficAY1lnjBw1PAOUgojetACUiImt4EPn09DR9fX3cddddvPa1r6VYLNLb28vf/d3f8e///b8HYNu2bZx++uncfffdvPKVr+Tb3/42/+7f/TvGxsbo7+8H4JZbbuGDH/wg09PTxONxPvjBD/Ktb32Lhx9+eOm13vrWt7KwsMDtt99+RHU7loPId9/zryQ3Rzy4bZzSzJOcecbp+BMnsX7TEKPzVTYXkuy8++v0nvJyqhs3kfnxrTS6TyOZ7CFM5KgncpjKLD3b7sDtHGJidILIyZFIlInSHWRy63h4psbZQwYvn6K8UCRW7sOvjJPv7mrdONa2adbqgCGRSrGtHjAYNcmvX8fupkvVtzit2yJWc6DXsH3XNs476VzM+OPUx+rM1nbx5GSKk/79LzG9+0GmZxz8DSeRzdgMZAaomgTF2hhbCsuvgCuObiO/7jQAgsCnXlwg293busO6N053spNCInfUNzqt1XYeMjwB7Kg12Zxe/TuAi4jIsfGiGEReLBYB6OrqAuC+++7D930uuuiipTKnnXYaIyMjSwHq7rvv5qyzzloKTwCXXHIJ7373u3nkkUd46Utfyt13371sG3vLXHvttYesS7PZpNlsLj0ulUorsYsHlSh0EuyoU9th6A7nqD86SV82gbV7luaox5RTJdcTUkrHOSmYYc/UkwTrLiNMBPR1+Iw/9TBRbBQ73km1eyvFhRy5vgLzzQa9SRuvadNVHsVKbyZtuqjWA9x1W4jsHjLrt2CMwfd96mOj5Pv6W0OoR8cZyKYYzfaxtZDmgd0zWIU4o9U9lCeKbEifwdRdP8Kv+rC+mx19I5QGT8KpzWCHJ5Mf6eCUDevpjbdad56sNeiJfHrNLmwrvniMp4k7cToXpxwozk0y1NeHxQRR5OFF0yQii3p94aiPbTzei+vqFikiIvLc1mSAiqKIa6+9lgsvvJAzzzwTgImJCeLxOIX9LmUH6O/vZ2JiYqnM/uFp7/q96w5XplQqUa/XSaUOvLT9pptu4r/8l/+yIvv2XLr71/OT6XvYRJlGb46HEyNscjfTKIdkNsWZLs7QEatzak+G2ekas84Wzj7rVBYe/R4TlZMpvORU9uz6HulSgvmFMmW7k/6hUwmrk2TqZbaFddadfBLVQidRGBBfdzbEy6QYwrIsLMsim80SFQp09/Yy7wd0V6s0mwEdqTjNSoCVCJmrTtMZBqyzc6S6MjRO6WDBTjHu7+DsDZfw6NgMM9UcXjzktI3r6V4MT80oImpOQVQjmzljab/LczUGNpyKs3gbk0bcJZNuTXBZqTxOV9erjulx96II9+gatkRE5AS0JgPU1VdfzcMPP8wPfvCD1a4KAB/60Ie4/vrrlx6XSiWGh4ePyWvNz1SZ37WdSsWnmh5gXb1EX8cY882Q9V6GSnGMuDvL0w/upDL5IIVcB/M/vp3Z2SImYZGavJ9m5DH31DYKZ2ymmeojVtlJBzUmnniIejqH2+ORbp7N7PjDxBMnkciU8GtFojBJuVGiXJzFGIvHR+8nEQT0hQFPj5YYOqmPXQs+leYY65x+at0xnMEB5vf8hEb/S+lycvQ1Evx8IosJi0TM0JN2SQWj1AIwJuTRhT2cnOliIvSp1Z5e2u96fYJmsxVe65UyVsyiVvMxJiKZPPBecytttOGzIXXgfe9EROTFac0FqGuuuYbbbruN73//+6xfv35p+cDAAJ7nsbCwsKwVanJykoGBgaUy995777Lt7b1Kb/8yz75yb3Jyklwud9DWJ4BEIkEicXzGxtjVOZJOkt0Dw5zWGyPdPJlkLkes6ZHuSBMFhlS+k1LmZeQqzzC49ZU0ZqboyvWyEAvIFAokavOE0RSu10mtDF7Wxw4iGvEk6XVbsTrrVK1O6oOnM1a06cieTBgLiSUCxht7sBJZvFwnAwkX+5lRHu9eR4YZJgcH2eE0iAWGXX6W9GAv841duNEQ/vgkox0ezeIOEplZml6BU9b34frNZftXKLwCzyqzIbuFlNs63mHgkUp5rSv/gMbCKF1D6559aI4pA9hHObZKREROPGtmGgNjDNdccw1f+9rX+O53v8umTZuWrT/33HOJxWLceeedS8u2b9/Orl27uOCCCwC44IILeOihh5iamloqc8cdd5DL5di6detSmf23sbfM3m2stvpgguTmM0j0w8DGAnY+jteTZOO5g2QGs6R7LPyOONn+DIn4Ak7/Opphle5XvxSrK0205WzCdJJ0fz9jiQEGzjiZ9EkjRLkc82GEXZ/AyWxgV8NjV7FOfS5iZscctbk9PPPkNoqTm2iOd9A9U8D7WYNGsBG2N8gnNpOczsN4jORYkkLfadjlCln7PBJhL930cZbrcVLfOrqTr+D0obPJWD30rz+XdHoj6fRGSu4wg4k49aC+FJ4AZiZ20DNwEtD6HLTijIiIyOpZMy1QV199NX/3d3/H//t//4+Ojo6lMUv5fJ5UKkU+n+fKK6/k+uuvp6uri1wux3/6T/+JCy64gFe+8pUAXHzxxWzdupXf/M3f5Oabb2ZiYoIPf/jDXH311UstSO9617v47Gc/ywc+8AHe+c538t3vfpd/+Id/4Fvf+taq7fv+Sn6Dk3o2MzXRGqhuDNS8gJMycfx6nTBq0vBTbDRzjDtpmmOTGN/jm1/+Cl2DG3nsqbtJ1rcRm7KYay6wvsNmz855Zn7+U8Y6tzDoGR7dNke+q4eh6k42veKlzE4+Te/g6eypTtOcTzDSM0LMGKKBNG7/AMnJURIbNwLwVL3MlpNPpzvVIJO+ENt2qc6M4Z4+zPjEv3Lqlv/Aj56c4cyBLspT41j2vgxfDSMGE3Fmn7XPJopw3NYtV0rTk+T7+hEREVlNayZAff7znwfgda973bLlX/rSl3jHO94BwJ/92Z9h2zZvfvObaTabXHLJJfzFX/zFUlnHcbjtttt497vfzQUXXEAmk+GKK67gD/7gD5bKbNq0iW9961tcd911/Pmf/znr16/ni1/84gtiDigAK0rz+Pw0fZku5hYqWA6ARVTzsZMOofGwrALV3Q9hihZWUKfq2hTOfjkXnnMyD2/bQU/ow3qX7bMRW09xKM5E1B+ZZ9Mmj5GgSezUETp5knA+STpWwcmfgh0P6CHH/HyZ3u483pNPkjj5ZAC8xbqVZ+eoWGW6Uyky6c3Ytku4sACOi2+XiCVa3aSRMTgW8KwuMXtx9iZrv1mcTGTYf1KnMAhx3OM7F9NU06c3vma+KiIichysmbPCkUxXlUwm+dznPsfnPve5Q5bZsGED//iP/3jY7bzuda/jgQceaLuOx0Of08PTzk7yQLFcp7czSU9HHK/UJEiGhIBLjNkf/oDul70OtyfN/KRLYsJmyp2huqtILChjmwTpjaeQTvdTiSKadgHXyxLmN7ApnSEWexm1+iip1Hqa5Rpz9TnsoI/ebBNv507ii12oxhjCMKRWKrIwuUA61yDm9GAim8Bv4k1OEnXYxKw8kZlmtlgjG7NYmBwn3z+wtF+VICTrHNijPDs5RmffvnKrMQypEkb0aQZyERHZz5oJUNISVGosjO0gysdIj5aYn5+g0ZgmFhnqUZVqYw+99hiVZJKpTC/O00/xs7ksuW6Lsj/D3GA38zvvJ3JSFPwiDzxVYfvDjzNTjnhlrMnJwzmIHCg7WAmbyPMxQYBxQvxmk/T8PHZvHwAmCKjOz1GZnqRZq2DFMgTVIl41xGMOf8cO3E2bqO18CjtcIOwYZGJyGteyiBc6sG1nab+mvYCNqThBFLQm69y7v36RRKI1YDzwPJyYgoyIiKw+Bag1xu+N6JoaoLMzTY9dJ5HuptA7TCyWoRLMMz33CMkwomPaovukUyjd+xCdqWFesTXP/O4nySRs8r1J8oVzWOjow5vcxUu7suzZMc2GpzLUmg0SiT7qxQpWNmS2NkmJBXqy3eze8xC9G4YIm03C0VEAatOT0NPDhi2n8Njj99I3cAa5nj68PaNktp6FnUrCfCeVhM/G4dPYOVPFAjKFzAH7ZlkWE5UJBjKLXX1hhGXva3KqzM1q/JOIiLwgKECtMY3pOaLJCVKJOI1iGu+ZR4kH63CSLuW5cbzGBPNPzdGVtJn51jcZm5gm547zo3/azab1gyxM7Ka2a47qSIkzT08x7XpUNp5O4sl7YShN12t/GYBw9HE6hk6iON3ESUyRz29k16NPUk8CmX1XyDWrKdZtPol6fRfTXp4N/UnCYhHLdXCyGfzJKZzuTqi0rnzcG4cq8w0CP1raTqnpsVAJKVcaZOse4OHVfbLdyaUyxkTLBp0fD80oImZr+gIREVlOAWqNGRraxLZde4hiAQ07R1dnE5PMQFeScG4ns16N4b4c6152Nl6jysy6TWSjftKTj7B+uEC3V8afbdKT96lO7mH+mR14s3M052fxBgcoPrUdgNrcLkzDMDlTIpYLKc40aY7uxutIkDThUn2sUpH5p35GGFVpTMepUKce+Tj9QwSTs/h7xgl7DEGq1apkABNEhMai0Jde2s58zSHOBDGmiWdbH0srMU8qfdbxO7gHMdbw2agJNEVE5FkUoNaYMJai0NtFKtUg1d9PtjaD1RjHqqVYKE5RmNlF8oL/QMKb4slZi42nbOXxn47Rt3Ej1VwX1jM/p6Orj9mz1tER62DOCcgvTBHbsBFz3tl46S5MOSDIdOLluphpljhlw0kU/QBT78TLJPD6u5bq07A8TK5KLH4Kc5WnOH3kJJrlCaKO1jipMOMRdBiGCqe1LgQwhtp8k02ndS/br3pjN4n8MPHAIp0eWVy68Tgd1cM72psTi4jIiUsBao2ZGl8gV6vhlyokrAWqu3/O+tdfhpuLU33oSfqGNjOYtdlV6qU3s0DUaNKMIlLlItXIEOtxyfhbqFsOiWpER+9GYmO7GCh0s677TCzLwq/V8dNNjKmQtWtkqDExMcnmvhRBs0mG2lJ9ArdOf9e5lGd9crk0w7lh8D1I9wLgpWoESUPMjjFerBMDCh3xZaGk2pwn4aRwnOMzm7uIiMjzpQC1xoR2kwV3ho2TNRrzO8kOnoFZGGVhtALj48RPjlHcvQuT3oJb3cPM6BPUTY7CYA9BoZ+o9BTzzQX6My+lXp3DclwafpFCahOPzz/emgG84uNFUzRqDTYPnU86nYbQJ2bnyA904u53JVwzHsO2Y0Rhc9+A7wNabFqPG35EdabBGfu1PhkT8Ux5klO7T6Xm10gcIkQ1azXiqfRB1x1LmvNcREQORgFqjdnzzBMYk2WiPkF642uZXTdCozvJxPTPmOndzNBQgkZqPbmgQWr9ebjbZ3CaLgu1NOWogTVXBK+TzLTPk+NFenKdjFVDZucqdNc76YhnibwGnuvT4UJHbl9oCYNgWXjaX6VSIZPNQLMCsX1X2EVRgGW1BoL7jQA34Sxvfao+STK1AceymKnPMNxx8Jsw10oLFPqP/U2D9zfZ9OnXBJoiInIQOjusMb39GSjNkD/pHFJuGtIW1dHHsatl1nW7pPw6SatMaXaa0akSQb1AJpiGRBK/MU6s5BGlIyJ/Dqc7w6YNA0z82MLtbbJ+XR9hqYm1IU8jLNFO+0slDMnGXKhOQee++xT6/hyZxDkAeGWfbO9+V/A1J0kkerH8fVfWHW680cHW1cKIyaZ/xPVsx7wf8LL8gdMtiIiIKECtMWZqHHoHKO9eILJcKpUF1p0yQmmyROaUIfKME7fzREETKzVIcTZNxnMZPGmEqR0LOF1ZNndvJTW8gV27dmNZFlFosBc/CZEXEcs5mKrBWpzQ0m80cOOJ1o33DmE2ihh2LQjMsi68KAqw7TiNik8ste/jFoZNwrBGItGP8ZtHfTwmmz6b0sdm7NQmNCZLREQOTgFqjdl05mt4fGw7BA7xDgvbOETNJhMzE5xR6eSJ2jzG90m4NvWuLAtBk1oYEs5UqNGkEkK8GVKt1Fkfd6jOPAWxDHF3//mWDJ43TT7fajnaMzbB4OAQpjR90DoFXohHk0y8AEFxaXlUrWKlWiGkXGyQzsUJI4PnzeL5c2QzJx+7AyUiInIMKUCtMaWESzHeTTLeIEqlsLs2MleZZGqwn5PXn8KGuXkSXhISBaqhRSybxC9WKdiG/nScZj3GcCpOvCPNfHWB4ujPaGCzfuAlS/cb9P15XDeDbbfCT6URsCkTZ7588O61arFJI16hO3kS1Ipsm9tG2k0T7Rmj0RUyuvMp5nyPwZTL2PxOIqsfN9bJXGkXoTFMeBFu4NCV7Dro9vdOfyAiIvJCoQC1xvgzz3DGgkN5106SqVNobP8ZVjZOf7NMvjRHfWoa10mR3fQKqjsmMIkYFtBsTlDyIzallt8KJVyoEBlDvG+QqOTh5OI0gilsO40xBmM8jPGIoiZR1Pr/vYwxmMin7lWIOxbGeJSa83Tlh+hJdeOlm8w3ffI9A7il3QwkY9idF7Kxe9+4oj0Nj3M7XBKHmWG8Vlwgne9c8WMpIiJytBSg1hgrnyPwbCq2TefIBjKxLkYjn+GpR+js7aUxvkBy3TkYIIqWt9qYsEzcXZwRPIqwsKiWaqQ6sq2xUH5ELO9A0BrGVKvtxLYSmKiE580RhAt43r5xQc1aDSseUlmYJ5Fr4nlzzNZ3M9S1Ds+bY35uBmvAxo6VcGP9xBN5nEZjWZ28yBw2PAF4jTqZggKUiIi8cChArTGddo5ndvyIZKGXrnXrKU7XGJ+b5zW5BOy6G5PpJspsorx7DqtoiKUMQTWkMVUjFY8wkUtg6tS8CRLZDAsLNXo7Rghm6xBEy17Lsiy8ikNH3+k4JMhk4iST++ZwapYmyPf0UpybxI43SSQGiLldJJODlGbrZFL9uAWbZHIQq1JlqtykL7d8YLbm+BYRkbXo+N6ZVZ63hJsg7jrk+tYvLYvN7yQRlSG3DtJdRFWfIANOR4pEXwoyEeMsMLxpE07Gwu3OYDKGWI9NlLQga+N2p4gNtLrWwrAOONh2gjAIcVyXanGBdKGwrC4mMti2Qz0yJG2LsZlHGeg6jcp8g1jCIZGy2X8qBC+ISLjOcThKIiIix5ZaoNYYUy1j5nZQ746z+6lHqM7VyRS6Ye98RYtTDxBBgMGr+NiJBWynE6KIqN4g3tcHlRKN+iSelyKVbj1nYaqGZUGxPE481kcysZFn9uxm3fo8pYUGlr18uoHynAdWjXk7ZDjuYOZnifVspTpdJ9uZpDZVw3UHiCKDbVlEGgguIiInCAWoNebpiSazjU7OfvmFZDNxfh5bYHOsCiwfWxSVPGaJ052I8bAX0NudIpxfwEnksLNZqJQYe+Y+0sle8l2L3XIG8n1pAidBKpkmnc4Qm4nRN5Rl3i6R711+KxUTxsn3pahUSuQSw/hOctkcUL63QCp2OqV6QEfSpVhfPuGlOYJAtXes1sHUwoiko05AERE5/hSg1pik8ShYAX55humioT5ToxxNYnV3Q6OIH87guFNUiiXmoxiZIM7OqTF+obOP4uwYSTuiPjlOaXqKualxgkYBrxIyPzHGzm176JrvoFp/injMg3CBpptlbm6OYrGEiSeX1aVULBLG4iQzSeYac2xKdWMig2WBCUOwLSzLotTwGSqkDghQM35Az3PcKqU8N0O2u+eg66Y9n3WJ+PM7oCIiIkdBAWqNGTp5A6NPP0rnwBBPVBuc6hQwc7N0bjwN5nZQb3jEEn3UQ0NUNwQxFyce56StZ1J6qkwq3k+QSjOw5RQWGvfjTbl0n7aV0EnSv8lm3YYBRsceYWjwF7Asi6dnqnR1ZbC8Bp1dy+dpsr0mhe5u7KniUitRZaFJtpAkmJnG6W6Vj4zBsQ9sKaoEEb3pg99bb6/D3X8vMuAeZLsiIiLHmgLUGjMxtov5+Xm8Rx9lPvBp1gwZu0o0N4eZmyMMKoRMYWNRKHTgxkI6s5nF+8i1wkazXiOwPbKmgznTwMlmKU/VSKRdwtDDwj3sPemeba4+x3ldZ0JxlDCIcGI2YaOB3Zt87ieLiIisQQpQa0zSjTGyZQtzGzZzfibJwlSNgjMO3d0E1e0E3Vsoz8VxwgTpfJzR8cfoTHTx8+mfk6pNEvMsygmHtBPgLNgkO1vzK3lNj1QhQbn8ENnsSQD4YYR7BGOM/Mgn6SosiYjIi4cC1BozMTrGFC5mZpZK0aVW9MjYZfAmMZN7iLw+omqdZDbOQtNnZ7VGCpfhzCnU4iUcu4fIjTPVrDFf3EEtkaY4Os/sxCydVhdBcRux/KVY5XnmGxUKWZ9npqCyMEuHs/wqvMrCLEmrRi4WO/BWK5bF3ikMjIFKMyCT0MdNRERODDqjrTEbe7qxaw3O2LAO3wvx0gGZ0IPufhoL4OQHIG9TNRBreDSSAUmTYEPSphJLkHBj1OIO07UFegodRLbDqes6ebpSp5mexmtWyecqANjsYX1qPY1ald5ClmRi+RxOpRhkkjZ9XadDdQYyvVBaXPmsQDVbaTLcufwqvucSBj6Oo4+oiIi88OjstMZUEr0UansIZus0qz62YxPUA6BOULIh5uN2p8CL8OoBVtiATA/ezEPEKxXqlWmyJ21mfM9DuMkNpOMFoDWlQMJ26Uj3MpQbaS3DMJIbYb42RqF/8IBxUQu1GIX8YOtBaYIwtQl7ca6oKAqxrH1XyBkD9n4DvktBSNY9/Dyu5dlZcr29z/OIiYiIrDwFqDWmPLOTPiYoT5WplCLSWYtmKgf1CUKriBWUqI7tZnqqwVR9nj0Lu3nFYMTuso9dzNCYXyCTrTEXdpIqFnEKLvMTY4zO7mQo7uOFDvMTYwAUK5PM12KUpqf2VcCYpbmeSjNTmMVuOrs4RWnKIZ2LMT9hUZ98nETiVBrFMYoLrVA1H+y7jcszXsBIzGH+MIPVLcvGtjVzuYiIvPAoQK0xfdkMhZN+CWwHa6pGoW+xW6y4h3r/RmLZPmw/RimWYCB8iNKeUzjn7A1UrICE38nC7mcwm7Zy2oYCkz/ZRkdvH50DQ9iTo/T29hAEaTo7hyh7ZYbym+lMteZg6hwYoh5GjDd9bAsCP6DRaWHyi1MbRFVKUYGwKwVAOVugo3Oo9e9wcYbz/L4uvFgY0ZVNHb8DJyIisoIUoNaYVMaCg7XKBE1wWy08zapPLJUhqllYzVno/jeEO36Mu2EL7H6aqhfSn3agWMY5qwdjDFE4CpxKPN4HwFxjjg25Da3tRYYnqg2Sjs3mdOs1KvUKiZ5OYonFVqWazUIYp5BKEFaq1PIJOlKL6xIBABtTy28k/Hw0woi45oASEZFVogC1xiSSB44bmvZ8jB/RCCxcL2Ci6mPSMFmcIkz1M9X0KfmGhh8wF1lU/QAaHnNBRNxxmR57jFqyj7kwASFUmz6zfkSy6VMKQppBxJmZ5dMUBJ5Hdu8tYJ4lnJ/D7swfk/1f2mc/oP85ZjEXERE5VnQGWmMatZBmWAOgMtegYSKm/IDeUgNjmhinSTI0FL0FslEDJ9lFxrUJLYjmZunr6mbKsQnHp8ikU6StCYw1QD7lkXFak22mXJu0bZF1bfoTMebLRzAOaf+xTFGEZbeC3pHc7+5oBJEhbh9+ELqIiMixogC1xiRzGZK51lgiY2A+63BWKoMVPEKtax12PcJOzTI1vp3OQg5/ehqrmaVS2kUwO0ums4v67Cy18YeJZT3yya2Ua1PEoymSVppYrBObiI5YjIzTCk6HupkvQLHmM16skirVqTkNFmxDNF+n6VikalWqzQAvjNjSlz0ux0dEROR4UIBaa3JDS//c2+hjWRaN5gR1PyIVpaFZodJssrl7mIWZBWZni7iTYyTsTrxkhpG0od5j4S2cSzZ7CrumnmBd4SU4ToN4vJvd5T2s71gPtKYj4DBXyi3UPU5NVylFg4wMdmLZFl41RdBlSKcz7J6rMZhP4jpqLRIRkROHAtQa4/sLNJtTeN4MlXKThSDJXCMiHtUAmzCq4VWfpmZBd6aHn1UnyaYT0H0GgdtJNLyVfD4knLGgsTggvFbm5OGN+P4YlmVjMNhWK/DUFhbI5AuHrI8xYIVNTKwVnvZp/TsyRuFJREROOApQa5DrdpDJbMEJG4QZh650AhZ8sPpJFTawa/xu0vF1TM4lGciPkG3mCDo3UquVCaMGvl/CBF3EU635mYIwJBaL4fut7e/fZed7zUMOFl9izLJGqiCs47g9K73bIiIiLxgKUGvMfDlktlgDdlKb95nPOjRiNsnxaarzBivfYHTbPHYixmR1nFi5zNzCPI2eXpIkWBjfTc/QKUT1x3ETrSvr9p9h3BizNDlmO/YfKx7486RipzzfXRUREXnBUoBaY+Jxh0w6BoBfD+lIOmQSLlbCJhZ3sawF5onxinPOoZJeT291js6kQyXmEM4/Tu/A+eS70zBfx+9eByy/Um6iOsFgZvCI62NZHDBGqjVZub30bxERkRONAtQaE0/E6XBaV7SNViy29KRxATvj4iSSRPNP05NO4k5PM9OYoNcuUTJQizfxK70MJCeoFw2m6WHZFo1mHdeyMSYCDH7kE3fih63D/owBrxkSy7Su2IsaDazkkT9fRERkLVKAWmv8GUxjFIDIC3EaDsZvEDFL5Ht4xUeIZ3tJbFyHXfRIJbI0SxXi3QUakxXsoRhBUMXq6aRZGWOsWKQ/n6bR2EM83g3B/LKXO9wUBgBBWGVuZjcd+Ri1mkUwNoHb33fMdh/AiyJczUIuIiKrSAFqjUmnR0inRwDIVhfo6irA5KMwdApm1zM4ZoFY97l05NdR272bhmuTSg+SGzmD2tw2UqkRGo1xYk4HQWkSmk1yuTypKANND7s+D05m3wvW56B6YIuSVZ+HahxTfIpkephMZiMAXgziHRuO6TGY9gJ6Y/roiojI6tFZaI2pV8o0KmUAKnN15mM1nNFthEOvoLTrSULHUK9UCLw6Hd3ddMYPvGFvFHnY5XGwY63xS7YLtsNCc4GOZOfye+3ZzgH33jPGgNN6TuTEsDrWHdN9fjYvMiQ1NYKIiKwiBag1xrKSWIuBxnIcLLsDyGGaCYKwRqzzDDLpTrKd3TA5AUkboohmEBJz9nV7RVjYqTw0G0TxDkgVKPklRnIjS2XCwMfOdEGqsKwOzUqFeOcApHJEiY7l1+wdZtJNERGRE4UC1BqTzMRIZmIYY8h6HoW+NJQdqrO7yW4ZJo5LPRdnvuaTJcD4AW5fH9ue2kVvs8j8xBhNb4pororTHaM8N0vJA9tvUqxOMl/b95GoFYsMnLTlgDo0axU6enoPXkFddiciIi8CClBrTKPq06z5TPo+yVJAcbqGmTPUJst462JMFiOqwRSTMyGNuTlicYNjRewaH2v12NXr+E0fv+JDqkylaihnDM16nY74OmynsPRa2a4ClXkf8JfVoTTTBKs1CWdtLiB18sGvujtWNxIWERFZbQpQa8zeFqi5WpNON0a+N83EXY/ivuzf0pz7EeXE6SRT82TsPEG9QW/Gwc13spAq0p1ycfOdNJtVwtAmGszRZXWx4eTBZZNpPhcTxsn3pvGCiHwQJ5Ha72O033bqfkg67hxkCyIiImubAtQaNzpfo7pQwY1lCSyHhp0mMFViJkZo4lSNgxXF6Uok6Boewu3spFYtkk50sGtxG+2Ep/0t1DwKqWd9hPZrdZqv+fRmE0e5ZyIiIi9cClBrTLiwQFgqMVNsUJr1GOpxSK7vw7UmiQ1soF5Okiy4JOw0QT3Oulwcsg71uIXT0dHaSG0W8qfjl6eIu+1Pern3Vi+VZoPB7L7nP7vLLggj4q6ulhMRkROPzm5rTAOfsWYRJ+6zZYNLwp8gveUMwqhJzD1wygKAZq1KLBHHchfzsl+HeIZmxSfZcfQZ2vPmSCT2TZoZVSrYudxRb09ERGStUAvUGuNkuthY6KU+XcHNJqmWPbKFNLXKBI7dDRw4e3jo+cTd2AHbiiKD8zzmUzLGx7b3bTecnyc2PHzU2zsSQWTQJOQiIrLa1AK1xiT9IvOTT5Ke20Os/DSmOIqdX0dkfCyrFWYiE+17wuL4pmXjnCyL2foshXjn86rLAbd5MWbZ6xyLi/CmfZ/e+IFhUERE5HhSgFpr0l0sdAxT6BrB7tsCuSEix8EKfaomSTpmU2oY8ukYUa2GvXfc0/6MoepX6Yhnn9dUA4ZnPfc4TKLZCA1pzUIuIiKrTGeiNWr/3NNojBP3HGYo0J12qTXBtizSjSpOoXBE21jxSomIiJzANAZqjWkGIdWmT7wZ4TR8al6A3fQJA0PVC3H8gGZgKNV9cjELy7JaLUXPCjd7u9+OdgqD/bchIiLyYqMWqDUmjAxeaPDDCC+ICELww4ggYnGZT19HnMgY3IONto4isGwMBs/ziMWOfjzRAe1Nug+eiIi8SKgFao1Jx10603EKUUghmyBIuXSkY4QJlyAdJ2OFWJnMsnSzrKWoMgmZHowxNOoNeob6DnyRw/C9Jm4sfvCxU+rCExGRFwkFqDWm1JjDeEWaTY/qPHjOPK67kYBnAAhNA8tOY4JDbMCvQSqJVw7IJJPYdnuNkLvnFghTaSjPMBWlSdWbS+s8PyS+3+M9TQ/qK/sRa+5/haGIiMgqUYBaY2ajDJvynRS9OrHKLKm+rbhRGjvWmgnchB5M+US1EH9iAu+ZZwimp7GGR5Ztp1kL6Mq1//aXvCYvGeinVpvDTXWzMdW6VYuJIvxEjHhqv1u3JIKl9SIiIicSBag1pln3KVZCyrMNgnqNSoeNPb6LWrqD8nQFE0Cz3MDv2Ug62Q2pXuopcK0sTNUwCz5eKiAKzfMaQG5MhL3fGKtwfh6ns2vpccMPScQ0xE5ERE5MClBrTF8uSSHmYlmQqiWxsklyVoXZdI7hRAyLDhYmAuIpQ24gR7ovDVGcQl8agNA2BB2dNJtNjN/emKV6GJFYDF1135CKOUvromqVeHf30uP5mkdnuv377ImIiKwFaiJYY7pircy7bLy2ZRFEBsv3iCVaoSUYGyM+vP6A54dBlQCXuNN+uJn0fPoWbw5crIfk04e+gs8LIpL7BSwREZETyZoKUN///vf55V/+ZYaGhrAsi69//evL1htjuPHGGxkcHCSVSnHRRRfxxBNPLCszNzfH2972NnK5HIVCgSuvvJJKpbKszM9//nNe85rXkEwmGR4e5uabbz7Wu/b8GIPF4szglgUGoloNJ58/oGgYNSj5DbqSXW134RnTmqATwAsjEq4CkoiIvDitqQBVrVY555xz+NznPnfQ9TfffDOf/vSnueWWW7jnnnvIZDJccsklNBqNpTJve9vbeOSRR7jjjju47bbb+P73v89VV121tL5UKnHxxRezYcMG7rvvPj7xiU/wsY99jC984QvHfP+ej70NUlGxjJXvABNhuc/qoY1CjG3jR+FRtUAtpzmfRETkxWtNjYG69NJLufTSSw+6zhjDpz71KT784Q/zxje+EYC//du/pb+/n69//eu89a1v5bHHHuP222/nJz/5Ceeddx4An/nMZ/i3//bf8t/+239jaGiIL3/5y3iex3//7/+deDzOGWecwYMPPsgnP/nJZUFrf81mk2Zz3+X7pVJphff8yEWlCtbQeqAILF6ZtzdelScg3XPI5049vYNYMnnI9UUvoDvdWq/4JCIiL2ZrqgXqcHbu3MnExAQXXXTR0rJ8Ps/555/P3XffDcDdd99NoVBYCk8AF110EbZtc8899yyVee1rX0s8vq+F5pJLLmH79u3Mz88f9LVvuukm8vn80t/w8PCx2MXl9g6CKo8fNhQFvoe7d1/CJritaQXCMFw2B1QUhsRTaToHhg75l+/pJ9fTRxR5WNbRz2AuIiKy1p0wAWpiYgKA/v7+Zcv7+/uX1k1MTNDXt3zmbdd16erqWlbmYNvY/zWe7UMf+hDFYnHpb/fu3c9/h55DVKlgp2IQBZhYatm6/Wce92o1EulM68Fi6LKwqFarZDKZpXKlmWlyvb2Hfr39Rq03m9PE4l3LCzxrFnJNSi4iIieyNdWF90KVSCRIJI7vhJFhqUyUrOF2b2Ku6NOZdcFbXGeipQjlNxqk84UDnh8EwbL74JkoxLYPPSh80vPpT7TKG+Pj2ProiIjIi9cJcxYcGBgAYHJyksHBwaXlk5OTvOQlL1kqMzU1tex5QRAwNze39PyBgQEmJyeXldn7eG+Z1VSZLVGfr1J6+mkaDOKUFpgMoN8ts1Bs4JYrTJg5OvwGs7OzlItFwlirC88uFqk3HeYDi1Q69RyvtFwjNAwm9jVY7t/CZHz/wAHrIiIiJ7ATpgtv06ZNDAwMcOeddy4tK5VK3HPPPVxwwQUAXHDBBSwsLHDfffctlfnud79LFEWcf/75S2W+//3v4/v+Upk77riDU089lc7OzuO0N4eW7c7Ru6GTXE+CzIZBBk4dpnOwh0KhQKGzQC6fhbrHxs0jdHd3U8jn6e7upru7m85CoVWuUKCn59DjptoVzMzgHKb7T0RE5ESzpgJUpVLhwQcf5MEHHwRaA8cffPBBdu3ahWVZXHvttfzRH/0R3/jGN3jooYd4+9vfztDQEJdffjkAp59+Om94wxv4nd/5He69915++MMfcs011/DWt76VoaEhAH7jN36DeDzOlVdeySOPPMJXvvIV/vzP/5zrr79+lfb6IOZ2QHYQY8yygeCm2cRKxAn9gFjyIF2KhxiY5DcauM+jC9J4HnZcs46LiMiLx5rqd/npT3/K61//+qXHe0PNFVdcwa233soHPvABqtUqV111FQsLC7z61a/m9ttvJ7nfpflf/vKXueaaa/jFX/xFbNvmzW9+M5/+9KeX1ufzef7pn/6Jq6++mnPPPZeenh5uvPHGQ05hcLyFzzxC0IgIJ8aI7CpepUI0X8dPzuPtmMPemMdU6wd5YgCHGLdUmZ+lMDDURi2eNYnB87innoiIyFpkGaPrpVZaqVQin89TLBbJ5XIru/Eowg8MlcefJhos0N3dzc7pEgPZIpVHJshs7eVnUwHn2nHiIyMsTIxTGBiE4h7I9FLzJpgJLEZyI0ubnJ8Yo/M5AtTOWpNN6VYrVa32DNP1HjZ0t67i83btIj6yb3teEDFdabKu0N44KxERkdXUzvl7TbVACTwzX+epsRLRTA2/aZGbjwjCBn5gUS97xOeahOYgPbOhB26CRq1Byu06cP1hVIKQjNPapjGG58rcCzWPzsPcJ09ERGStU4BaYzZ0ZyhEkOjwqXZk6O7upunN4Ngp5kfHyXQlyAQ2LLTKL81Cvhh6Sl6JddnNbb3mtBewMbV4k+KgSDPqIJs49Een7of0dhzfaR1ERESOpzU1iFxa/Jk54n37rqKLwjqOs2+cl3WYG60EUUDM3tc6VCsVSWWfu5tx742Hfb9IxUuSTx2+handGxWLiIisJQpQa5BpNLHT6X2PTYRl7ZsE02CWzTre+sfBA02zWiGZzbbz6kQGXEcfHRERefHSWfAEEQUB1uKUBiYKwV65t/bZ2UuNSyIi8mKnMVAniMaePSSGhoAm0egEQZQFYwimp/AWJgAHKrtgdhqvsWvpecHUJJ4fHnbbnh/ixVotXH59jLDRxKu2rrAzQYhTyC8rr+s6RUTkRKcAtcb4/gLN5gS1mku9XqJWK9NojOFNl0metAG/MU99chvmvNcQuBFB0CSMB5iuDUCEZ0UE3RHQuqIuMBFBT3TI1wuNIWpGBMlWs1PKPRW7lCTekznkc5SfRETkRKcAtcbEYgUSiQHS6Y3U67Ok092ARWRNk80NE4ZdxK1pMrktAASNZ0ht/gUAoiggmZwind4IQGlmmt7BdbiHmUV8d8Njc8Elvl+XoGVVD1n+SKY5EBERWes0BuoEYSyIojq2ndw3UfjsU0T7TZgZhjUsZ9/klmHgHzY8AQSRWRae4PBddJOlJoN5TaApIiInNgWoE0gYNoisGA4OVGcJrDh2Mrvf+gqWfWzDTd0PScWd5y4oIiKyhilArTHG87BiB+l5NQZjAgwWDhHUZvGcDhLpfWOVosjHttu76a+uuBMRETmQxkCtMcGuJ3EdYG4HFEuE9Qkaz9xPOBYQOjl8O0/CL0HPyXiTE+Q79k2S2ZqVvJWIwiDAtp+7pajd4UzKWyIi8mKgALXGmFgGcr2YzhSYOcJymdiml0BvHGd9gsjpxM3sAcvCYJbmhgJoBA0ysV4AyjPT5Pv6j6oOh2qVmq006cy018IlIiKyFqkLb40KggDXbeXfqFTBXrxrdBAFONbBW5aKzSKFRGHpsfUck23O+QGF2IHbOlSrVKkRPOctXkRERE4EClBrjLV41Vyz2SSR2HfD3r03DQ5NuBSgnn1PvMhEOEfQbbdX0Q/pfNZ4qzAyOLY66kRE5MVNAWqNiQ0OAuB5HjHXPaA/zQ99bOu531ZzlNNdFus+ObUyiYjIi5wC1BpljCGcmFwKVHuFtSqxdDs3Bz60g411Kjd8OhIHDp0rNXw6khpSJyIiLw4KUGtZ4C916e0VVaq4udxBi+/fpffs7r122AfpwputePRkEwcpLSIicuJRk8EaEy4sEOyZJHQa+NXqYjOR2ReImk2sZPKYvb7u0iIiIqIAteY4hQLu+jiWUydWLhPfsIGgBl5pXxlrhWa/bCcsaVi5iIi8mKgLb42KKvumLjhSRztwfH8Hy2YNPyTu6qMkIiIvHjrrrTFmsVkoLBZxu7r2LT9IOFqJwHTg6x+4bLzYYDB/7LoNRUREXmjUhbfG1GtzBNEszWiCWq3VAhWLddKk2ipwiH63PfM1JkoNXLtVrrTQYMGtHvJ1vChiPgixasGy5ZVmcNDyK9VtKCIishYoQK0xVtRBRy5HMN0knd4IgIkiLGqHfM50uUk27jKYS7Gxp3Vz4YUgSaEnc8jn7G54bIrHiGnSTBERkQMoQK0xfjMknXNhv9uw+M0GsWQSaC4bpLT3yrxSw6c3H5GN7wtMz9W9Fxpz0PA0XqzT9KOlx5ExxBz1BIuIyIuLAtQaE0WGaGqSWG/v0jKvXieRzeIHxQPKz1c9CqkYxeYUXfGjm2DTGMPuuTp+FDGYT5LO62MjIiIvbjoTrjXGYIIAk0wuDSj3m03ShU4832BM6y8MAizbZqba5KTeLJVSBI619BzMvgHpBxNFhigyPDNXIzKG4c60rrQTERFZpAC1xlQaITNln3rYoGbVW8uKDcrxBvWGT7HcxJmr06iWqXkRUSxg91ydsUodP2YRr+//nPohX2e86eE2Ika60rp5sIiIyLMoQK0x2ZRLH3HK2RTd3WkAFvwUhe40tVqcRDHJSHeaYlhiKp7m5IE8AJFtM5gukEgsf86hhDWHkbRuzSIiInIw6pNZY8JiCadQOHSBxW65WjMgndx3nzwT1XGcQwemZ9OsBCIiIofWVoDyfZ93vvOd7Ny581jVR55DVKng5PPLlh3sirqJss+6Qmrf86JGWwFKREREDq2tABWLxfg//+f/HKu6yBFw8rmDDv6OIg/bjoNl0QxCXOdZBYzBstTgKCIishLaPqNefvnlfP3rXz8GVZEjYWezBEGA67rUvZAd0xVG5+vsmJ5n92zIZLnJzpkq63LxZc9r57Yu1SAkbStsiYiIHErbg8hPPvlk/uAP/oAf/vCHnHvuuWQyy2ez/t3f/d0Vq5wcKJGL8cT4PLbtkDMem3uzLIRlkgWDXU2QGOxhZCDH/ETlkNuIwhDrMAFpxg9Yn4wfcr2IiMiLXdsB6q//+q8pFArcd9993HfffcvWWZalAHWM+cBA1qVQKGDvF4KisIFVczHZ1jinvbOQH0zgNXHjhw5IkQFHo8hFREQOqe0ApQHkqyuXjDFbNcvC015WFGEdMPgJwijE3m/8U+B5uHFNUSAiInK0ntdAl72zXssLhDEHbXmab86TT+y7cq8VoA7dAqXGJxERkcM7qgD1t3/7t5x11lmkUilSqRRnn302/+N//I+VrpschYMNFq/5NZJuculxFIY4ruZQFREROVptn0U/+clP8pGPfIRrrrmGCy+8EIAf/OAHvOtd72JmZobrrrtuxSspR+gQTUfW4v9ERERkZbQdoD7zmc/w+c9/nre//e1Ly37lV36FM844g4997GMKUKvsYEGpnSkMRERE5Lm13YU3Pj7Oq171qgOWv+pVr2J8fHxFKiVHbyXCkoa1iYiIHF7bAWrLli38wz/8wwHLv/KVr3DyySevSKXk0OqVMqXpSRYmxpf+EtnsAeX2D1LtdN/pogAREZHn1nYX3n/5L/+Ft7zlLXz/+99fGgP1wx/+kDvvvPOgwUpWVirbQa63n0J397Lltdoc8Py78IpBSCF24FQIIiIisk/bLVBvfvObueeee+jp6eHrX/86X//61+np6eHee+/lV3/1V49FHeWIGDAGgyGKwqO+791CEFI4yFxSIiIiss9RXct+7rnn8j//5/9c6brICqmWxoisIrVaEwC/MY7VMXTEz7c0EZSIiMhhHVGAKpVK5HK5pX8fzt5ysgqs1nQF9coU3QPnLM31FAssUqmRVa6ciIjIieOIAlRnZyfj4+P09fVRKBQO2kJhjMGyLMIwXPFKypFrdeFFS+FJg8JFRERW3hEFqO9+97t0dXUB8L3vfe+YVkhWVtkv0xHvOOLyylsiIiLP7YgC1C/8wi8AEAQBd911F+985ztZv379Ma2YrIxis8i67LrVroaIiMgJpa1LtVzX5ROf+ARBEByr+shRMk0P6yA3CDbGYLdxRZ7Gj4uIiDy3tq91/zf/5t9w1113HYu6yPMQVarY2Q6iyMOydKNgERGRY6ntM+2ll17KDTfcwEMPPcS5555LJpNZtv5XfuVXVqxycqAgCHDdA9+2qNbAXp8mnJkl7u6bsqCdWcgjY3TLYRERkSPQdoB6z3veA8AnP/nJA9bpKrxjr9lsEj9IVx20roI0JgLL2W/pkY8Kn/UDumNqvRIREXkubZ8toyg6FvWQI9RsNikUCsdk29UwojceOybbFhEROZEc3f0+FjUajZWqhxwhYwy2ffi3bf9uu3a68EREROTItB2gwjDkD//wD1m3bh3ZbJYdO3YA8JGPfIS//uu/XvEKynMzJoL9gtL+3XbtdOGJiIjIkWk7QH384x/n1ltv5eabb142FufMM8/ki1/84opWTo5MFDVwnORqV0NERORFo+0A9bd/+7d84Qtf4G1vexuOs2+w8jnnnMO2bdtWtHJyZMKwjm0nlh6rC09EROTYansQ+ejoKFu2bDlgeRRF+L6/IpV6Ifjc5z7HJz7xCSYmJjjnnHP4zGc+wyte8YrVrhYLQUi53ly2rNGoYgUxEvUmo55FLQpYqDdpBA2qgUP4rPIlr7X+2ULdx0VEROSItB2gtm7dyr/+67+yYcOGZcv/9//+37z0pS9dsYqtpq985Stcf/313HLLLZx//vl86lOf4pJLLmH79u309fWtat0KrkN3KrFsWc1EuMkEThzCpEt34FJIJRitzHBSoY+YvfzKuoV4a72IiIgcnbYD1I033sgVV1zB6OgoURTxf//v/2X79u387d/+LbfddtuxqONx98lPfpLf+Z3f4bd+67cAuOWWW/jWt77Ff//v/50bbrjhgPLNZpNmc1+LTqlUOmZ1SyYPPdbJ86bBymG7ra7VMAoPCE8iIiLy/LU9BuqNb3wj3/zmN/nnf/5nMpkMN954I4899hjf/OY3+aVf+qVjUcfjyvM87rvvPi666KKlZbZtc9FFF3H33Xcf9Dk33XQT+Xx+6W94ePiY1e/ZM7/vz5iAoOGRSGcPWSaKQixb46JERESej6Oadvo1r3kNd9xxx0rX5QVhZmaGMAzp7+9ftry/v/+Qg+Q/9KEPcf311y89LpVKxzREHY7fbBI/TCtV4Hm4cXXfiYiIPB9tt0Bt3ryZ2dnZA5YvLCywefPmFanUWpNIJMjlcsv+VosxEdZhJtpsBaiD3QpGREREjlTbAerpp58+6P3ums0mo6OjK1Kp1dTT04PjOExOTi5bPjk5ycDAwCrV6ggc4RV0gddUgBIREXmejrgL7xvf+MbSv7/zne+Qz+eXHodhyJ133snGjRtXtHKrIR6Pc+6553LnnXdy+eWXA60pGu68806uueaa1a3cIZjIgGUD5jnnfYqCENvRDYNFRESejyM+k+4NE5ZlccUVVyxbF4vF2LhxI3/6p3+6opVbLddffz1XXHEF5513Hq94xSv41Kc+RbVaXboq74XG1OrYmX6gfETlLUuDyEVERJ6PIw5QURQBsGnTJn7yk5/Q09NzzCq12t7ylrcwPT3NjTfeyMTEBC95yUu4/fbbDxhY/kIRVWrY6zsw/rGbPkFERET2absvZ+fOnQcsW1hYoFAorER9XjCuueaaF2yX3bMZ38M4BoIY9mEGkIuIiMjKaPts+yd/8id85StfWXr8a7/2a3R1dbFu3Tp+9rOfrWjl5Mh53jROvBPHak2i6Uc+ju08x7NERETkaLQdoG655ZalOY7uuOMO/vmf/5nbb7+dSy+9lPe///0rXkE5MsYEhAZiVmvm8fnGPF3JrlWulYiIyImp7S68iYmJpQB122238R/+w3/g4osvZuPGjZx//vkrXkE5AoszGARRgG21MnEjaNCXXt379omIiJyo2m6B6uzsZPfu3QDcfvvtS7c8McYcdH4oOX5CEy514YmIiMix03YL1Jve9CZ+4zd+g5NPPpnZ2VkuvfRSAB544AG2bNmy4hWUI7A4K0EQBcRsd3GRpioQERE5VtoOUH/2Z3/Gxo0b2b17NzfffDPZbOvGtePj47znPe9Z8QrKkWiFpTAKcRdboAxHNjO5iIiItK/tABWLxfjP//k/H7D8uuuuW5EKydFohaXABCTUhSciInLMHdU9PZ566ik+9alP8dhjjwGwdetWrr322hftzYRfKKIwxLZjq10NERGRE17bg8i/853vsHXrVu69917OPvtszj77bO655x62bt3KHXfccSzqKM/BGB/bTuA168RTqdWujoiIyAmv7RaoG264geuuu44//uM/PmD5Bz/4QX7pl35pxSonzy0M6/hehVy8B78+QaI7s9pVEhEROeG13QL12GOPceWVVx6w/J3vfCePPvroilRKjkwYNmg0RonHO7HtGGEY4LhH1SsrIiIibWg7QPX29vLggw8esPzBBx+kr08TNx4vYdik3thNJrOFZ89YEJlI0xiIiIgcQ203V/zO7/wOV111FTt27OBVr3oVAD/84Q/5kz/5E66//voVr6AsFy4s4M9PU2+Mkklvwpt5BqczD+yb+6nYLJJL5FazmiIiIie0tgPURz7yETo6OvjTP/1TPvShDwEwNDTExz72MX73d393xSsoy1m5DL47QWf2dUvLgtozy8pUvArrO9Yf55qJiIi8eLQdoCzL4rrrruO6666jXC4D0NHRseIVk4MLgjLZzKmHLWMwWJa68ERERI6Vox5xPDU1xfbt2wE47bTT6O3tXbFKyaG5sQKRidh/ovHIRIRRSGhCwihsrRcREZFjpu0AVS6Xec973sP/+l//iyhqnagdx+Etb3kLn/vc58jn8yteSdlnoTZKsT62fKHlEDdxZptzTNQmyMazq1M5ERGRF4m2A9Rv//Zv88ADD/Ctb32LCy64AIC7776b9773vfzH//gf+fu///sVr6Ts05UZpiszfNB188lx1mXXHX4D6toTERF53toOULfddhvf+c53ePWrX7207JJLLuGv/uqveMMb3rCilZP2HNENhI1uMiwiIvJ8tT0PVHd390G76fL5PJ2dnStSKREREZEXsrYD1Ic//GGuv/56JiYmlpZNTEzw/ve/n4985CMrWjlpz3NNnhlFobrwREREVkDbXXif//znefLJJxkZGWFkZASAXbt2kUgkmJ6e5i//8i+Xyt5///0rV1N53kLPx43HV7saIiIia17bAeryyy8/BtWQ48H3mgpQIiIiK6DtAPXRj370WNRDnidjzHMOIg88j3S+cHwqJCIicgJrewyUvDAFR9C6FAUBjnvUc6eKiIjIIgWoE0SzViOeSj9nOd3iRURE5PlTgDpB+M0Gbuw5xjcpPImIiKwIBagTyHO2LmkSTRERkRWhAHWCeK45oERERGTltD2i+Prrrz/ocsuySCaTbNmyhTe+8Y10dXU978qJiIiIvBC1HaAeeOAB7r//fsIw5NRTTwXg8ccfx3EcTjvtNP7iL/6C973vffzgBz9g69atK15hERERkdXWdhfeG9/4Ri666CLGxsa47777uO+++9izZw+/9Eu/xK//+q8zOjrKa1/7Wq677rpjUV85hCO6kbCIiIisCMuY9kYWr1u3jjvuuOOA1qVHHnmEiy++mNHRUe6//34uvvhiZmZmVrSya0WpVCKfz1MsFsnlcsflNecnxiinA0ZyI4csszAxTmFg8LjUR0REZK1p5/zddgtUsVhkamrqgOXT09OUSiUACoUCnue1u2l5HjSIXERE5Pg5qi68d77znXzta19jz5497Nmzh6997WtceeWVS/fJu/feeznllFNWuq5yGOrCExEROX7aHkT+l3/5l1x33XW89a1vJQiC1kZclyuuuII/+7M/A+C0007ji1/84srWVEREROQFou0xUHtVKhV27NgBwObNm8lmsytasbVsNcZALUyMU0r7GgMlIiJylNo5fx/1nWWz2ezSXE8KTyIiIvJi0vYYqCiK+IM/+APy+TwbNmxgw4YNFAoF/vAP/5Aoio5FHUVEREReUNpugfr93/99/vqv/5o//uM/5sILLwTgBz/4AR/72MdoNBp8/OMfX/FKynPTIHIREZHjp+0A9Td/8zd88Ytf5Fd+5VeWlp199tmsW7eO97znPQpQqyAMAmzbAYLVroqIiMiLQttdeHNzc5x22mkHLD/ttNOYm5tbkUpJe7x6jUQms9rVEBERedFoO0Cdc845fPaznz1g+Wc/+1nOOeecFamUtKdZqxFPpVa7GiIiIi8abXfh3XzzzVx22WX88z//MxdccAEAd999N7t37+Yf//EfV7yC8tyMiRa78J6DpdnKRUREVkLbLVC/8Au/wOOPP86v/uqvsrCwwMLCAm9605vYvn07r3nNa45FHeU5HPFtXI5uyi8RERF5lqOaB2poaOiAweJ79uzhqquu4gtf+MKKVExERETkhartFqhDmZ2d5a//+q9XanMiIiIiL1grFqBk9WgOKBERkeNLAepFwkSRBpGLiIisEAWoE8CRDCIPfA83FjsOtRERETnxHfEg8je96U2HXb+wsPB86yJH6Ui68ALPw40njkNtRERETnxHHKDy+fxzrn/729/+vCskx0bgeaRyh38PRURE5MgccYD60pe+dCzrIc/DkXThhUGA4x7VrBUiIiLyLBoD9WJhDJYGkYuIiKwIBSgRERGRNqlPZ42pV8o0KuVlyxxHb6OIiMjxpDPvGpPKdpDKdqx2NURERF7U1IV3gjC6UbCIiMhxs2YC1Mc//nFe9apXkU6nKRQKBy2za9cuLrvsMtLpNH19fbz//e8nCIJlZf7lX/6Fl73sZSQSCbZs2cKtt956wHY+97nPsXHjRpLJJOeffz733nvvMdijlRWZCNtaM2+niIjImrZmzrie5/Frv/ZrvPvd7z7o+jAMueyyy/A8jx/96Ef8zd/8Dbfeeis33njjUpmdO3dy2WWX8frXv54HH3yQa6+9lt/+7d/mO9/5zlKZr3zlK1x//fV89KMf5f777+ecc87hkksuYWpq6pjv4/MRmhDXVo+siIjI8WCZNdb3c+utt3LttdceMPP5t7/9bf7dv/t3jI2N0d/fD8Att9zCBz/4Qaanp4nH43zwgx/kW9/6Fg8//PDS89761reysLDA7bffDsD555/Py1/+cj772c8CEEURw8PD/Kf/9J+44YYbjqiOpVKJfD5PsVgkl8utwF4/t5pfoxbU6En1HHT9wuQEhf6B41IXERGRtaid8/eaaYF6LnfffTdnnXXWUngCuOSSSyiVSjzyyCNLZS666KJlz7vkkku4++67gVYr13333besjG3bXHTRRUtlDqbZbFIqlZb9HW+BCXAs57i/roiIyIvRCROgJiYmloUnYOnxxMTEYcuUSiXq9TozMzOEYXjQMnu3cTA33XQT+Xx+6W94eHgldqktYRTi2IcJUGuroVFEROQFbVUD1A033IBlWYf927Zt22pW8Yh86EMfolgsLv3t3r37uNchiAJcS2OgREREjodVPeO+733v4x3veMdhy2zevPmItjUwMHDA1XKTk5NL6/b+/95l+5fJ5XKkUikcx8FxnIOW2buNg0kkEiQSiSOq57GiQeQiIiLHz6qecXt7e+nt7V2RbV1wwQV8/OMfZ2pqir6+PgDuuOMOcrkcW7duXSrzj//4j8ued8cdd3DBBRcAEI/HOffcc7nzzju5/PLLgdYg8jvvvJNrrrlmRep5rASRxkCJiIgcL2tmDNSuXbt48MEH2bVrF2EY8uCDD/Lggw9SqVQAuPjii9m6dSu/+Zu/yc9+9jO+853v8OEPf5irr756qXXoXe96Fzt27OADH/gA27Zt4y/+4i/4h3/4B6677rql17n++uv5q7/6K/7mb/6Gxx57jHe/+91Uq1V+67d+a1X2+0hpHigREZHjyKwRV1xxhQEO+Pve9763VObpp582l156qUmlUqanp8e8733vM77vL9vO9773PfOSl7zExONxs3nzZvOlL33pgNf6zGc+Y0ZGRkw8HjeveMUrzI9//OO26losFg1gisXi0ezqUXmm+Mxh18+Pjx2nmoiIiKxN7Zy/19w8UGvBaswDtau0i5HcyCHXL0yMUxgYPC51ERERWYtelPNAyXOwrNWugYiIyAlDAUpERESkTbrufY2pV8o0FwfO769UnWKhFjvk84yJjmW1REREXlQUoNaYVLaDVLbjgOWlkk8hpzFOIiIix4O68ERERETapAAlIiIi0iYFKBEREZE2KUCJiIiItEkBSkRERKRNClAiIiIibVKAEhEREWmTApSIiIhImxSgRERERNqkACUiIiLSJgUoERERkTYpQImIiIi0SQFKREREpE0KUCIiIiJtUoASERERaZMClIiIiEibFKBERERE2qQAJSIiItImBSgRERGRNilAiYiIiLRJAUpERESkTQpQIiIiIm1SgBIRERFpkwKUiIiISJsUoERERETapAAlIiIi0iYFKBEREZE2KUCJiIiItEkBSkRERKRNClAiIiIibVKAEhEREWmTApSIiIhImxSgRERERNqkACUiIiLSJgUoERERkTYpQImIiIi0SQFKREREpE0KUCIiIiJtUoASERERaZMClIiIiEibFKBERERE2qQAJSIiItImBSgRERGRNilAiYiIiLRJAeoEEJkIy7JWuxoiIiIvGgpQJ4AgCnAtd7WrISIi8qKhAHUCCKIAx3ZWuxoiIiIvGgpQJ4DQhLi2WqBERESOFwWoE0AQBTiWWqBERESOFwWoE4BaoERERI4vBagTQBAFClAiIiLHkQLUCUBdeCIiIseXAtQJIDShApSIiMhxpAB1ArCwNJGmiIjIcbQmAtTTTz/NlVdeyaZNm0ilUpx00kl89KMfxfO8ZeV+/vOf85rXvIZkMsnw8DA333zzAdv66le/ymmnnUYymeSss87iH//xH5etN8Zw4403Mjg4SCqV4qKLLuKJJ544pvv3fBnMaldBRETkRWVNBKht27YRRRF/+Zd/ySOPPMKf/dmfccstt/B7v/d7S2VKpRIXX3wxGzZs4L777uMTn/gEH/vYx/jCF76wVOZHP/oRv/7rv86VV17JAw88wOWXX87ll1/Oww8/vFTm5ptv5tOf/jS33HIL99xzD5lMhksuuYRGo3Fc91lEREReuCxjzJpsvvjEJz7B5z//eXbs2AHA5z//eX7/93+fiYkJ4vE4ADfccANf//rX2bZtGwBvectbqFar3HbbbUvbeeUrX8lLXvISbrnlFowxDA0N8b73vY///J//MwDFYpH+/n5uvfVW3vrWtx5R3UqlEvl8nmKxSC6XW8ndPqhdpV2M5EaO+euIiIicyNo5f6+JFqiDKRaLdHV1LT2+++67ee1rX7sUngAuueQStm/fzvz8/FKZiy66aNl2LrnkEu6++24Adu7cycTExLIy+Xye888/f6nMwTSbTUql0rI/EREROXGtyQD15JNP8pnPfIb/+B//49KyiYkJ+vv7l5Xb+3hiYuKwZfZfv//zDlbmYG666Sby+fzS3/Dw8FHumYiIiKwFqxqgbrjhBizLOuzf3u63vUZHR3nDG97Ar/3ar/E7v/M7q1Tz5T70oQ9RLBaX/nbv3r3aVRIREZFjaFWnr37f+97HO97xjsOW2bx589K/x8bGeP3rX8+rXvWqZYPDAQYGBpicnFy2bO/jgYGBw5bZf/3eZYODg8vKvOQlLzlkHROJBIlE4rD7ISIiIieOVQ1Qvb299Pb2HlHZ0dFRXv/613PuuefypS99Cdte3nh2wQUX8Pu///v4vk8sFgPgjjvu4NRTT6Wzs3OpzJ133sm111679Lw77riDCy64AIBNmzYxMDDAnXfeuRSYSqUS99xzD+9+97uf596KiIjIiWJNjIEaHR3lda97HSMjI/y3//bfmJ6eZmJiYtm4pN/4jd8gHo9z5ZVX8sgjj/CVr3yFP//zP+f6669fKvPe976X22+/nT/90z9l27ZtfOxjH+OnP/0p11xzDQCWZXHttdfyR3/0R3zjG9/goYce4u1vfztDQ0Ncfvnlx3u3RURE5AVqTdyB9o477uDJJ5/kySefZP369cvW7Z2FIZ/P80//9E9cffXVnHvuufT09HDjjTdy1VVXLZV91atexd/93d/x4Q9/mN/7vd/j5JNP5utf/zpnnnnmUpkPfOADVKtVrrrqKhYWFnj1q1/N7bffTjKZPD47KyIiIi94a3YeqBcyzQMlIiKy9rwo5oESERERWS0KUCIiIiJtUoASERERaZMClIiIiEibFKBERERE2qQAJSIiItImBSgRERGRNilAiYiIiLRJAeoEYNBcqCIiIseTApSIiIhImxSgTgAW1mpXQURE5EVFAeoEoC48ERGR40sBSkRERKRNClAnAHXhiYiIHF8KUCcAdeGJiIgcXwpQIiIiIm1SgDoBqAtPRETk+FKAOgGoC09EROT4UoASERERaZMC1AlAXXgiIiLHlwLUCUBdeCIiIseXApSIiIhImxSgTgDqwhMRETm+FKBOAOrCExEROb4UoERERETapAB1AlAXnoiIyPGlAHUCUBeeiIjI8aUAJSIiItImBagTgLrwREREji8FqBOAuvBERESOLwUoERERkTYpQJ0A1IUnIiJyfClAiYiIiLRJAUpERESkTe5qV0DaU2wWKTaLy5ZFJlql2oiIiLw4KUCtMflEnnwiv9rVEBEReVFTF56IiIhImxSgRERERNqkACUiIiLSJgUoERERkTYpQImIiIi0SQFKREREpE0KUCIiIiJtUoASERERaZMClIiIiEibFKBERERE2qQAJSIiItImBSgRERGRNilAiYiIiLRJAUpERESkTQpQIiIiIm1yV7sCJyJjDAClUmmVayIiIiJHau95e+95/HAUoI6BcrkMwPDw8CrXRERERNpVLpfJ5/OHLWOZI4lZ0pYoihgbG6OjowPLslZ026VSieHhYXbv3k0ul1vRbctyOtbHj4718aNjffzoWB8/K3WsjTGUy2WGhoaw7cOPclIL1DFg2zbr168/pq+Ry+X0hTxOdKyPHx3r40fH+vjRsT5+VuJYP1fL014aRC4iIiLSJgUoERERkTYpQK0xiUSCj370oyQSidWuyglPx/r40bE+fnSsjx8d6+NnNY61BpGLiIiItEktUCIiIiJtUoASERERaZMClIiIiEibFKBERERE2qQAtYZ87nOfY+PGjSSTSc4//3zuvffe1a7SmnfTTTfx8pe/nI6ODvr6+rj88svZvn37sjKNRoOrr76a7u5ustksb37zm5mcnFylGp84/viP/xjLsrj22muXlulYr5zR0VH+v//v/6O7u5tUKsVZZ53FT3/606X1xhhuvPFGBgcHSaVSXHTRRTzxxBOrWOO1KQxDPvKRj7Bp0yZSqRQnnXQSf/iHf7jsXmo61kfn+9//Pr/8y7/M0NAQlmXx9a9/fdn6Izmuc3NzvO1tbyOXy1EoFLjyyiupVCorUj8FqDXiK1/5Ctdffz0f/ehHuf/++znnnHO45JJLmJqaWu2qrWl33XUXV199NT/+8Y+544478H2fiy++mGq1ulTmuuuu45vf/CZf/epXueuuuxgbG+NNb3rTKtZ67fvJT37CX/7lX3L22WcvW65jvTLm5+e58MILicVifPvb3+bRRx/lT//0T+ns7Fwqc/PNN/PpT3+aW265hXvuuYdMJsMll1xCo9FYxZqvPX/yJ3/C5z//eT772c/y2GOP8Sd/8ifcfPPNfOYzn1kqo2N9dKrVKueccw6f+9znDrr+SI7r2972Nh555BHuuOMObrvtNr7//e9z1VVXrUwFjawJr3jFK8zVV1+99DgMQzM0NGRuuummVazViWdqasoA5q677jLGGLOwsGBisZj56le/ulTmscceM4C5++67V6uaa1q5XDYnn3yyueOOO8wv/MIvmPe+973GGB3rlfTBD37QvPrVrz7k+iiKzMDAgPnEJz6xtGxhYcEkEgnzv/7X/zoeVTxhXHbZZead73znsmVvetObzNve9jZjjI71SgHM1772taXHR3JcH330UQOYn/zkJ0tlvv3tbxvLsszo6OjzrpNaoNYAz/O47777uOiii5aW2bbNRRddxN13372KNTvxFItFALq6ugC477778H1/2bE/7bTTGBkZ0bE/SldffTWXXXbZsmMKOtYr6Rvf+AbnnXcev/Zrv0ZfXx8vfelL+au/+qul9Tt37mRiYmLZsc7n85x//vk61m161atexZ133snjjz8OwM9+9jN+8IMfcOmllwI61sfKkRzXu+++m0KhwHnnnbdU5qKLLsK2be65557nXQfdTHgNmJmZIQxD+vv7ly3v7+9n27Ztq1SrE08URVx77bVceOGFnHnmmQBMTEwQj8cpFArLyvb39zMxMbEKtVzb/v7v/57777+fn/zkJwes07FeOTt27ODzn/88119/Pb/3e7/HT37yE373d3+XeDzOFVdcsXQ8D/abomPdnhtuuIFSqcRpp52G4ziEYcjHP/5x3va2twHoWB8jR3JcJyYm6OvrW7bedV26urpW5NgrQIksuvrqq3n44Yf5wQ9+sNpVOSHt3r2b9773vdxxxx0kk8nVrs4JLYoizjvvPP7rf/2vALz0pS/l4Ycf5pZbbuGKK65Y5dqdWP7hH/6BL3/5y/zd3/0dZ5xxBg8++CDXXnstQ0NDOtYnOHXhrQE9PT04jnPA1UiTk5MMDAysUq1OLNdccw233XYb3/ve91i/fv3S8oGBATzPY2FhYVl5Hfv23XfffUxNTfGyl70M13VxXZe77rqLT3/607iuS39/v471ChkcHGTr1q3Llp1++uns2rULYOl46jfl+Xv/+9/PDTfcwFvf+lbOOussfvM3f5PrrruOm266CdCxPlaO5LgODAwccKFVEATMzc2tyLFXgFoD4vE45557LnfeeefSsiiKuPPOO7ngggtWsWZrnzGGa665hq997Wt897vfZdOmTcvWn3vuucRisWXHfvv27ezatUvHvk2/+Iu/yEMPPcSDDz649Hfeeefxtre9benfOtYr48ILLzxgOo7HH3+cDRs2ALBp0yYGBgaWHetSqcQ999yjY92mWq2GbS8/lTqOQxRFgI71sXIkx/WCCy5gYWGB++67b6nMd7/7XaIo4vzzz3/+lXjew9DluPj7v/97k0gkzK233moeffRRc9VVV5lCoWAmJiZWu2pr2rvf/W6Tz+fNv/zLv5jx8fGlv1qttlTmXe96lxkZGTHf/e53zU9/+lNzwQUXmAsuuGAVa33i2P8qPGN0rFfKvffea1zXNR//+MfNE088Yb785S+bdDpt/uf//J9LZf74j//YFAoF8//+3/8zP//5z80b3/hGs2nTJlOv11ex5mvPFVdcYdatW2duu+02s3PnTvN//+//NT09PeYDH/jAUhkd66NTLpfNAw88YB544AEDmE9+8pPmgQceMM8884wx5siO6xve8Abz0pe+1Nxzzz3mBz/4gTn55JPNr//6r69I/RSg1pDPfOYzZmRkxMTjcfOKV7zC/PjHP17tKq15wEH/vvSlLy2Vqdfr5j3veY/p7Ow06XTa/Oqv/qoZHx9fvUqfQJ4doHSsV843v/lNc+aZZ5pEImFOO+0084UvfGHZ+iiKzEc+8hHT399vEomE+cVf/EWzffv2Vart2lUqlcx73/teMzIyYpLJpNm8ebP5/d//fdNsNpfK6Fgfne9973sH/X2+4oorjDFHdlxnZ2fNr//6r5tsNmtyuZz5rd/6LVMul1ekfpYx+02XKiIiIiLPSWOgRERERNqkACUiIiLSJgUoERERkTYpQImIiIi0SQFKREREpE0KUCIiIiJtUoASERERaZMClIiIiEibFKBERI6BjRs38qlPfWq1qyEix4gClIisee94xzu4/PLLAXjd617Htddee9xe+9Zbb6VQKByw/Cc/+QlXXXXVcauHiBxf7mpXQETkhcjzPOLx+FE/v7e3dwVrIyIvNGqBEpETxjve8Q7uuusu/vzP/xzLsrAsi6effhqAhx9+mEsvvZRsNkt/fz+/+Zu/yczMzNJzX/e613HNNddw7bXX0tPTwyWXXALAJz/5Sc466ywymQzDw8O85z3voVKpAPAv//Iv/NZv/RbF4v/frv2ENP3HcRx/Tut7WJszQWjG4ItU1EpSq0M7FNTolHibRkwICSZeOkhnD5otKFIjOolrdCgIEnYwaoaQh9KCiBQC19gOY7HGIvGQOH8H+X1pGb+f+2X9fr/xesBg+3z/vD/fz2G8eH++n616/f39wMYtvFQqRXt7Ow6Hg5qaGgKBANls1jre399Pc3Mz0WgU0zRxuVx0dnby5cuXX7toIvKPKECJSMUYHh7m+PHjXLx4kUwmQyaTwePxUCgUOHXqFC0tLczNzTE5OUk2myUQCJRcH4lEMAyDmZkZ7ty5A0BVVRUjIyO8e/eOSCTC1NQUly9fBsDn83Hz5k1qamqsen19fRvmVSwWaW9vJ5/PMz09zZMnT0gkEnR0dJSct7i4yKNHj4jFYsRiMaanp7l69eovWi0R+RnawhORiuFyuTAMA7vdzq5du6zxW7du0dLSwpUrV6yxsbExPB4P79+/Z9++fQDs3buXa9euldzz2/epTNNkYGCAUCjE7du3MQwDl8uFzWYrqfe9eDzO27dv+fDhAx6PB4C7d+9y8OBBZmdnOXbsGLAetMbHx3E6nQAEg0Hi8TiDg4M/tzAisuXUgRKRivfmzRuePXuGw+GwPvv37wfWuz5/OnLkyIZrnz59yunTp9m9ezdOp5NgMMinT59YXl7edP2FhQU8Ho8VngC8Xi+1tbUsLCxYY6ZpWuEJwO128/Hjx7KeVUR+D3WgRKTiLS0t0dbWRjgc3nDM7XZb33fs2FFyLJlMcvbsWXp6ehgcHKSuro7nz5/T3d3N169fsdvtWzrP7du3l/y22WwUi8UtrSEiW0MBSkQqimEYrK6uloy1trby8OFDTNNk27bN/+29evWKYrHI9evXqapab9g/ePDgb+t978CBA6TTadLptNWFmp+fp1Ao4PV6Nz0fEfnv0BaeiFQU0zR58eIFyWSSXC5HsVikt7eXfD7PuXPnmJ2dZXFxkcePH3PhwoW/DD979uxhZWWF0dFREokE0WjUern823pLS0vE43FyudwPt/b8fj9NTU2cP3+e169f8/LlS7q6ujh58iRHjx7d8jUQkV9PAUpEKkpfXx/V1dV4vV7q6+tJpVI0NDQwMzPD6uoqZ86coampiUuXLlFbW2t1ln7k8OHD3Lhxg3A4zKFDh7h37x5DQ0Ml5/h8PkKhEB0dHdTX1294CR3Wt+ImJibYuXMnJ06cwO/309jYyP3797f8+UXk97Ctra2t/duTEBEREfk/UQdKREREpEwKUCIiIiJlUoASERERKZMClIiIiEiZFKBEREREyqQAJSIiIlImBSgRERGRMilAiYiIiJRJAUpERESkTApQIiIiImVSgBIREREp0x+Y3BFJTX165gAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "logp = tree.sampler.get_log_prob()\n", "\n", "for _logp in logp.T:\n", " plt.plot(_logp, lw=0.1)\n", "\n", "plt.xlabel(\"Iteration\")\n", "plt.ylabel(\"Log posterior\")\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 7, "id": "92fed717-cfc6-4bd2-878e-c1d99c7d0002", "metadata": {}, "outputs": [], "source": [ "logp = tree.sampler.get_log_prob(flat=True)\n", "chain = tree.sampler.get_chain(flat=True)\n", "mpe_parameters = chain[np.argmax(logp)]\n", "\n", "par_manager = tree.par_manager\n", "par_manager.set_parameter(par_manager.parameter_fit, mpe_parameters)\n", "parameters = par_manager.get_all_parameter()" ] }, { "cell_type": "code", "execution_count": 8, "id": "a6bd6937-27ea-49d7-9245-87b225d844b3", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjQAAAG4CAYAAABSEX09AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABVIUlEQVR4nO3deVxU5f4H8M857KiAiCyyuy+ZYi6plUtmWrfNbpve0uzaBpWZpmYudSvten83W7i5FC6VmlqWtlim5QKioGgqu7KJLG6IqMgw5/z+QAgIjWHOM3OG+bxfr3n9fgzM93k43vLTs0qqqqogIiIismGytTtAREREZC4GGiIiIrJ5DDRERERk8xhoiIiIyOYx0BAREZHNY6AhIiIim8dAQ0RERDaPgYaIiIhsHgMNERER2TwGGiIiIrJ5DDRERERk85pNoLl06RJCQ0MxdepUa3eFiIiILKzZBJq3334bN998s7W7QURERFbQLAJNRkYGUlNTMXr0aGt3hYiIiKzA6oFm586duOeee9CuXTtIkoRvvvnmTz8THR2NsLAwuLq6YsCAAdi3b1+d70+dOhXz58+3UI+JiIhIbxyt3YGLFy+iV69emDhxIsaMGfOn73/55ZeYMmUKFi9ejAEDBmDRokW48847kZaWBl9fX3z77bfo3LkzOnfujLi4uL9s78qVK7hy5UrN14qi4OzZs2jTpg0kSdL0dyMiouZDVVVcuHAB7dq1gyyLGw8oLy9HRUWFJrWcnZ3h6uqqSS3dU3UEgLpx48Y67/Xv31+NjIys+dpoNKrt2rVT58+fr6qqqs6YMUMNCgpSQ0ND1TZt2qgeHh7qG2+8cc025s6dqwLgiy+++OKLrya98vLyhPwdqKqqevnyZdUdDpr11d/fX718+bKw/uqJpKqqCp2QJAkbN27E/fffDwCoqKiAu7s7NmzYUPMeAIwfPx4lJSX49ttv63x+xYoVOHLkCP7zn/9cs436IzTnz59HSEgI2t8fhZYBHTX9fbR0LiMReb98jgmz3oV/aHth7RzduxObP/0QQaNfhEubQE1rXziehOK4tZjz7/cR1kHcs96z81cse/8/GBX1BryDwjSrm3UgDnvWLcFb//cB2nfspFldANj123Z8/N5C/PP1d+Ef2kHT2rUd2bsT33zyAV5sE4QgZxdh7QDAgUsXsPZ8MT6e/jQ6BwcIbava1oTfsWDlRrz81v8hKFz7/43t3/0bVn/8Hpb8azo6hwdrXr++rbsT8M7ilVjyxivoFCamvV/iEvHOks+x+N//QqcOYdrX3xmL+e8vxruLPkIHjf+5qW/nb9vxwX/exW3PzIZXuzDN65/OTkPc8n+jpKQEnp6emtcHgNLSUnh6emIcAuFs5qqQCij4Avk4f/48PDw8NOqhfll9yul6Tp8+DaPRCD8/vzrv+/n5ITU1tUk1XVxc4OLy53+RtwzoCK8OvZpU0xKunD8FAAjrdgPCuvYU1s7ZopMAAPeAjnD31/YvBEPpaQBA9xt7oVvP3prWrq3oZD4AwL9jd/h37K5Z3QunCgEAN/TqjR439tasLgAUnDwBAAjr1hPh3cT/+XZydUdHF3dh7QDA6UoDACCiczgiuogL4bXlFZ8BAHTq0Qudetyoef1TBVXPr0+PzujTo4vm9es7UVAMAIjo3hl9uncW00Zh1b9bIm7sjj49e2hf/2TVPzc39uqNnr16a16/toL8qn+O2rbvjrbtuwlrxxLLE9wgw1kyL9A46Ga4wjJ0HWhMNWHCBGt3gYiIyGwOkgQHM4OTA6SqiSc7oetA4+PjAwcHBxQVFdV5v6ioCP7+/lbqFRERkViyBDiYORAkA3YVaKy+bft6nJ2dcdNNN2Hbtm017ymKgm3btmHgwIFW7BkRERHpidVHaMrKypCZmVnzdVZWFg4ePAhvb2+EhIRgypQpGD9+PPr27Yv+/ftj0aJFuHjxIp588kkr9pqIiEgczaac7IjVA01iYiKGDRtW8/WUKVMAVO1kWrFiBR555BGcOnUKc+bMQWFhIXr37o0tW7b8aaEwERFRc+GgwZSTgzZdsRlWDzRDhw7FX+0cj4qKQlRUlIV6RERERLbG6oGGiIiI6uKUk+nsNtBER0cjOjoaRqPR2l0hIiKqg1NOptP1LieRIiMjkZycjISEBGt3hYiIiMxktyM0REREesUpJ9Mx0BAREemMBPOnUOwrztjxlBMRERE1HxyhISIi0hlOOZmOgYaIiEhnuMvJdAw0REREOlMVaMwdobEvXENDRERENo8jNERERDrDKSfTMdAQERHpDBcFm85up5yio6PRvXt39OvXz9pdISIiIjPZbaDh1QdERKRXsvTHtFNTX7J9DdBwyomIiEhvOOVkOrsdoSEiIqLmgyM0REREOsNdTqZjoCEiItIZBhrTccqJiIiIbB5HaIiIiHSGi4JNx0BDRESkMw7QYMpJ1aQrNoOBhoiISGdkDUZoZDM/b2u4hoaIiIhsHkdoiIiIdEaTXU72NUDDQENERKQ3miwK5pSTfeDllERERM2H3QYaXk5JRER6Ze7FlFpMWdkaTjkRERHpDKecTGe3IzRERETUfHCEhoiISGdkSTL7HBl7O4eGgYaIiEhnJAcJkmxeIJHsLNBwyomIiIhsHkdoiIiIdEZ2kCCbOULDKSciIiKyLgcZkmzmJIpkX7dTMtAQERHpjCRLkMw8SEaCfY3QcA0NERER2TyO0BAREemM7CBBNnOERrazERoGGiIiIp2RZPPX0Eiqfa2h4ZQTERER2Ty7HaGJjo5GdHQ0jEajtbtCRERUB6ecTGe3IzS8bZuIiPRKcpA0eTXW/Pnz0a9fP7Rq1Qq+vr64//77kZaWJvA31J7dBhoiIiKqsmPHDkRGRiI+Ph5bt26FwWDAyJEjcfHiRWt3rdHsdsqJiIhIr6pGWMxcFAyl0T+7ZcuWOl+vWLECvr6+2L9/P2677Taz+mEpDDREREQ6o+UamtLS0jrvu7i4wMXF5bqfPX/+PADA29vbrD5YEqeciIiImrHg4GB4enrWvObPn3/dn1cUBZMnT8bgwYNxww03WKiX5uMIDRERkc5IkgTJzMspJaXq83l5efDw8Kh5/69GZyIjI3HkyBHs3r3brPYtjYGGiIhIZ2QHGbKZa2hkterzHh4edQLN9URFReG7777Dzp07ERQUZFb7lsZAQ0REpDOmbrtusIba+M+rqooXXngBGzduxG+//Ybw8HCz2rYGBhoiIiI7FxkZidWrV+Pbb79Fq1atUFhYCADw9PSEm5ublXvXOAw0REREOmPpEZqPP/4YADB06NA67y9fvhwTJkwwqx+WwkBDRESkM1quoWkMtRlcZMlt20RERGTzOEJDRESkNxpMOcGEKafmgIGGiIhIZ2RJgmzmOTSyZF+BhlNOREREZPM4QkNERM3KufzjQuqWnMwWUrchkoNs/uWUin2NWTDQEBFRs2AwVECSZGz/YJa1u2I2TS6nVOxrysluA010dDSio6NhNBqt3RUiItKAk5MzVFXBvOlTEB4arHn9lPRMLFgUrXld0obdBprIyEhERkaitLQUnp6e1u4OERFpZPSIYYi4Uftbonft2WexQKPJwXocoSEiIiJr4hoa0zHQEBGRcBnpacLbyM3NEd6GpcgO0GANjUadsREMNEREJIzBYIAsy3jxmaes3RVq5hhoiIhIGCcnJyiKgjdffAphQQFC24o7cBiL134rtA1LkWQJkpkH65n7eVvDQENERMKNuu1m9OnRRXg7i9d+i5T0TCG10zPFnG/TEFnW4HJKI9fQEBER2RyDoRKyLGNC5MvW7gpZAQMNEREhNUPM6EN23gkhdRvi5OQIRVHgf+vjcPLy07z+ldN5KN7zpeZ1G6LJtm1zL7e0MQw0RER2rHpUY/wL06zdFc04e/nBxUf7g/VUwxXNa16LJtu2zfy8rWGgISKyY9WjGrJ/BCTnlprXVy4WQz2ThtTj4rdUZ2SfACQJuZv/I7wt0h8GGiIiguTcCpKrgFPTy0shyzKeePUt7Wtfg/+Au+Hs0UbzuuVnC1G8/2fN6zZEkmVIspkjNGZ+3tYw0BAR6ZzI0Y2M7DwAEoy5O4W1oQKYNmsOQkJChbUBAHvidmP1qpUo3Pu90HYsQXbQYJcTp5yIiEgPKiqrDqV7Yvo7wtvyv+0JOAtYSHsxLxlnkr7H8BEj0bNXb83r17d65XLMnTMHYaHah6e09HQsePddzeuSNhhoiIh0ytmx6lC6iS/PREBQiJA2Du/fh02rl8PZyw+uAhbSVpQUAbDs1Qd3jhyJiIjemtffvTvWcoFGg0XB4AgNERHpyc1DbkfnG3oJqV1pqMCmNTJyNy0UUh+Axa8+SE0TE57SMzKE1G2IJGuwy4lraIiI6HpSjlnmEsSsEwXC23B0cgZUBX0efh4t27bTvH5x+iGkbl2PGa/PQaiAaaDaYmOr1tBMfMr2743iomDTMdAQETVSpaGiak3LtDct2m7OMXEjAwUncgEArdq2g2dguOb1y06dBACEhIaiU2exVx/k5ORAURTc8/QUtGmn/fRZYXYmtqyI1rwuaYOBhoiandxj6ULqni4qgKIoGD1pCtoEBAlpo7bMg/uw77t1ePuV54S2I0kydkS/Lqy+LMt4fpLlRk3atAuGf1hHzetWlJdrXvNaqg7WczCzhlGj3tgGBhoiajYMBiMkWcaCqc8Lbaf7zUMQ3OUGoW1Ui9+0FnMnPIBw/7ZC6scdzcDSzb9i0D0Pw9PHV/P6+Zlp+H3XVvgPuEvI2TC1leVn4lzqXqyYZ/t3OfGkYNMx0BBRs+Hk5ABVUdB24Dg4eWq/BflSfgpKDv+IopxjmtduyJmCqnuQRvW/ERGdxKw/qaisxLLvdyBu8zoh9QEAkozCvT+Iq1/P3CfuRZi/j+Z103ILsGDtj5rXJW3YbaCJjo5GdHQ0jEb7GpIjsgctw/vCza+D5nVVoxHnj/6Ez9+conlta3F2dISqKJgzdy5CQ8M0rx8fvwfLli7FvOkvIzxY7DRdbMJ+LF25Gnf264mIjtpvc999OMNigUaWZchmLuo19/O2xm4DTWRkJCIjI1FaWgpPTwHHfROR1VSczRNSt7LsNFRFwUPPT0NbAYtO60s/lIhf1q9Cau5JYW1kFZ4CAIwceSciIiKEtLFs6VKEBQehayftQ2ZtWVdv9k7NFbM7LP1EoZC6DeGUk+nsNtAQUfNjMBgBSUL+lv8Kbadtu2C0C9d+0Wl9hbnHIcsSnlywTHhbaWmpQupmZmZAlmVMiHpFSP36ZFnCxIUxFmmL9IWBhoiESs3Jt1hb+afPAKqKG8Y8gxY+As5USTuA7F2b8b/XX9S89vU8GvUqfAPFjAgdTYzHrxvX4KmJE4XUr9ZhyL1w89J+XUttJXmZOHFgJ8a+OB3+Ap5X3vFMrF+ySPO6DeEIjekYaIhICIOqQpYlTHjzI4u33cKnHTwCwjSve/H0SaiKgnmTHkZYgPY7guqLO5yGpRu3wjcwGEHhnYS0UZyfB0VR8GKfrghq5a55/YSC09iQcQLHdmzSvPa1+AUGI6i99s+rvPyy5jWvRZI0OFhPYqAhIjKbkyRBUVTMvHswQtpYZp1aXGYePttzFHuXzhXazqibIxDRpb3QNoCqKbRl8jZ8MPMF4W3dFuyHHj5eQmqvT8/F9FlzECz4pOD42N344rOVeG96lNB2SJ8YaIhIqBE92qNXsPZbqK/ls7jDwo7Z3xsfjxWfil/PUq16G3rU9FkIDBITBpIS47Fupfg1J8EWOCn4+LFMqIoqtA1L4ZST6RhoiKjZGXHHSNzYW8yOnRWfLrPYuqDsgmIAQGBwKMI7dhbSRv6JqnupjpdcEFI/+3wZZFlG1NOWOyn4jVdfQpiALeKpmccw//3FmtdtCAON6RhoiIgaqcJQYfF1QbIsY2bU02LbkIBXdxwQ2kaLTkPh4O4ltI2Ks3koz9uPUcNvQ5+ePTSvvzM+wWKBRnaQIZsZSMz9vK1hoCGiZic9PU1I3bzcXFh6RkNRFKE7hKp3Bom6mqAsPxNnU+JxMeM3zWsT1cZAQ0RCpReesVhbx4vPWeQixK7DH4B7a7FbkAHgdHYacg/sFr9DyAJXE9z36OPw8RW7liozNRk7fv4BqRnHhdTPOJYtpG5DJFkyf5eTLGnUG9vAQENkh3IrxN8afKKiHLIk4blVlrvDp9qt9z4i5KLF3PRkHI79FanbN2pe+3rueHAsWrcVEwZy0lOwd/sWBNx8N1wEjNBcyM/EuZR4fLv2M81rN0SWZYx/YZpF2hKJa2hMx0BDpCMnszKF1i/MPQ4ZwMJTuULbqa3/3Q/Do4340QwAyM9MRWr8b9i16Uuh7Tw1ZSYCgrS/K6i+g/vi8P26L7D1q9ViG5JkFMR/L7SJll2Gw8G9tdA2Ks7m4nJOAibPeF3IFvHM9HR8/N5CzeuSNhhoiHTAUGGAJMsWO4H25ZtvQJBHC6Ft7MsvxrqUHOz7XuAtztcw56HbEdpW+78849NzsOyXBAQEhSC0g5hdR7UVnMiFqiiYNVvMNnQAiI+PR8wnyzC5fw8Ee2r/v4m9J4qxLjUHZWnbNa99LUNHjESPG3trXjdhT6zFAg1HaEzHQEOkA07OTlAVBQ9HvgrfIHGXHqYlJWDrulUYEhaAG3y9hbVT7cujWcJ/p9qqf7+RvTujd7j2Vx9UGI2Qt+/HW1Oe07z29YSEhqJzFzFnuOTkVG3bHhoWgB4CQiAArEvJRrfhD6BF67ZC6lc7nZuB7IRfkSloUfjxzAwhdRvCk4JNx0BDpCMRtw5H+249hbaxdd0qofXr8w0MRqCgY/vrKz4h5pbtas4ODlAUBaMfGidsTUttx1OPYN9vv+CZf4o/wyXzbKmQulnnLkCSZKRYaN2RLMuY+vwki7RF+sJAQ0RCGIxGSLKMj14Tf2x/fWn5p4TUzSw8A1mW8eP6L4TUvxbffqPh3FL7BbsAUFaQidK0eEz9ZZ+Q+tUecvNBWwcnoW2kGy7jlysleG7qLLQL1n6KLjszHZ9++B/N6zZEcnCA7OBgdg17wkBDpCP5x8UOaVePYBwT9F/jtRVdLIeqKMJ2HDUkNz0ZR2K346n/bRDaTlTPjghsof1FjvUlFp/BV1kFKE74UXhbYaMnwtU7QPO657MOoyBuE3q5tES4o5vm9ev75UoJBg+7A1179tK89oG9cRYMNFxDYyoGGiIdqF4UbInRDFkCXvl5r/B2AECSZeE7jhoyb9Z0hIVov24nbm8ClsasxC0BvujuLf7CTYNRxVfHLXPNgqt3ANz9tB/VKD9boHlNooYw0BCZSMSCw8KCk1AVBQ89Pw1t24lbQJt+KBG/rF8Fn4hRcGoldlHw5aIslKTtwZ1/HwtvC6w3AYDs9BTs2bYFYSHB6NpZ+3U72bli1+jU5+RQdTCayF1OsbG7sXLFSqR+8baQ+taQlSlmUXCO4BHU2jhCYzoGGqJGMhoNwhcc9ho8DOECFwVXGgyAJON00hZhbdQmyzJ+2iD4DJUG2pzwTJTQNrJKy4TWr5Z/8RIA4I6RI9FL0GWbALBy+XLc/uRL8PbX/kLH7CMHkLB5jeZ1G1IJFbIsY/ZLz1ikPZEkWYNdTmZ+3tYw0BA1koODExRFwZsvTkRYoLZrDeKSDmPx2k3CD9Y7d6oQUBXMmTsXoaFhQtuKj9+DZUuXYmJAEAJcXIW2Ve1wWSk2nS4WdojblTPZqMjdj5nxv2te+3rS0sSMOAB/bNtu7R+EtiEdNK9/tvAEAOBkZYXmtes7Z6yEoijwiHgQji213yJuOH8SF34XfA3FVRyhMR0DDZGJRt06AH26a3uomsFQiSXydxY7WG/kyDsRESHuv/irLVu6FAO8WqOzu9hD/KpVqgo2nTkl/BC3N577B8LaiZ9G23ngCJZ/u1X4tm1JlrFhvrjrAiQA/ys7Kax+fY4t28LJU/tziFQLhDJqOgYaIh1wcnKEqijo9+jzaOWr/bB/tcLUJCT/vB5paanC2qiWk5MtvI36HCUZUFUMvf9ReAnYWZWXmYr9v/2MsHZ+6Bou/rDA7JNFUBQFEaP+jlbeYg6lK8pKR9qebQi/6ym4tdF+l1PJ8cM4GfstWnQbDYcWYraeV7tyKhNXcvfi7K7FQtuxBEmWzB+h4eWURM3DmTxtb9wtKRK/26SVbyBaB4YLq3++IAeyLOOpiROFtVFfzuXLFmsrr/wyJFnGb9+sFdaGLEsYP/v/hNWvT5JlJG0Ruw0dANzaBKCFgF1Ol89U7XJyaNEGjq3EjmoZL57BFVWFZ8RDwqacSn//RvO6DeEaGtMx0FCzYzAYIEkyNv/fdGt3pdEMhkpIkoztH8yydlc0JQN4J1vsuqCG9H80Ep4CRrpOpibh6M/r0GHIvXDzEn/hZkleJk4c2Cns9wGAE0cTkLL9GyR/9paQ+lUklCZ+LrB+XW5BveDso/1/GJQXpFgs0JDpGGio2XFycoKqKngsajp8A7X7SyA1KRE/rVuJ1OM5mtWsll98CqqqYM4jIxEq8I6l+NRsLNsaX3VOi6BtwNXi4vdiacxKjH38Cfj6az+V0ZDUo0ew5Yfv4eEbiNZB2v+Fdq4gB5Ik49gOyywMrRbW51a0bd9NWP2UbV/jjVdfQliw9qEpLuEAlqxag5f6dkNQK7FrqQ4UnsaalGwYzotZr2MoLRRStyGS7ABJNvOkYDM/b2vsNtBER0cjOjoaRqPR2l0hQfrcOgwdut+oWb3KSgO2bpDxxPR3NKtZ38iIrohoHyisPgAs2xqPsJAQIee01JZ57BhkWcbqzyx7d5Qkyfjlg9eEtvHi9NcRKOBo/fqSEuKxduWnOHtC2+nT2s4XV+1CGjX8NvTp2UNIG0tWrcFtwf7o4eMlpH41g1HBmtQcnNn5P6HtWITsUPUyt4YdsdtAExkZicjISJSWlsLTU/yJn2T7HB2rtm179XsUTi21XXB6uSgVZck/Iy2/WNO69WUWnLp6Tkuk0HZq6zTsfrhbYHoGAM7kpOHEgV3C27l1+B3ofmNv4e0YDAZIq5YLD2jNhZND1aLwgWOj4Omn/X8YnDmRhYT1SzWvS9qw20BD1FROrXzh7KXt0HxFaQFkScJTH4pbzFrba6/PQYjoKafY3fhs5Upk/PqN0HYa8mbUkwgL9Ne8blzSESxetxnHM9I1r92Q4sKTUFUFcyf+HeEBYnY5xR1Ox9JN24TUtpbwm26FX4fumtc9cTTRcoFGlqte5tawIww01GydOK7tYtSTuVmAJOHU9g80rVvbvKfHIryduIscY39PxdKvtyAkNBSdu3QR1g4A5ObkQFUVzHvtVSH3KjUkbl8Clsaswqhb+6NPd+2n1AyVlVi64XtMjxJ3WnRDRg3ohYjO4na/NbdA0xxIDg5m35bN27aJBBP9X7c5x49BkmW8P1PM8fcTJs+Af5C2oxtHD+zF5tUrMHpQH0R01f601mqGykp88o2MZyeJPaitNlH3KjVE9F1LTo6OUBQFQSOehGtr7UeA6ivNOYLifZuFt0PUHDDQkMWoV+9Ceu2Fpy3SntZbXau37PoHhSCkg7Z/QReeqNo5lZJ9QtO69Z0oPgtFUTDz7sEIaSN27VhcZh6+iD+KCc+Kv0G8vtTjuULqZudX7XJxbe0Pt7YhQtqorfyc5XbVkM5wUbDJGGjIYqSrdyF1vGcS3H20P5a82tnM33Fi10bNt7oaKw1I+WUDFkx9XrOatcmyhAnzFgmpXV9IG0909hd7amvumfNQFAXTXpuN4BDxO4IAID4uFms/X4knZs4X14gkI3O9wPpWkpohZidVdp7YkN5sybIGgYZraIiE8r1hIDxDxK7fOLFro+ZbXcvOVh1DP3nG6wjWeEHt/r3xWL3iU4Tccg9cPcXtCCrJTcfpI3F4btUPwtqoLzgkFB07i/3zrpaXmwNFUTBrzlyEClj0HB8fj5hlS/HKTMuEtMR98fh8+SdC26iorIQsyxj/gri7nMh0PCnYdAw01OyolQahZ5EEhYSiQydt/4LOOnYMkGTk7rbMegm38MFwcBM75VRxLg+VRUfx4rP/FNpOQ+4YORK9e4u5fDNm2VIMGzESN/TqLaR+fZ8v/wSpOeIudsw/dQ6KoqDvw8+jla/2W50L0w4iZet6zesS1cdAQxZXVpgttP7lkqpTd1+Y/joCNVy8m7g3Fl+vXoWpz4vb4TJ5+usIEridel9cLNZ9vhKXs2KFtVHfmy8/g/AgcVOMtcXuP4TFq7+2SFuWYDAYIMsSnnznY+FtBUfcAh9BpxEz0DSBpMEaGolraIhQflr73SLlZ/MBScbvy9/UvHZDbh2m/eFnGz5fgWmz5iBE4+mGhL3xWBWzTMjoT20ncnIAqJg7exbChV99EI+ln8Rg9JBB6HNDV6Ft1bZ49ddIT0sTUjsnR/trL67HyckJiqJi0pSZaBcsZhHy74n7sPGL5UJqkxm4KNhkDDRUR2WloWrqY9NCYW3c+o8X4SXgFM9qJ5KTkPTjWhzP1HZ7eP7VnUgi1oQcO54JWZYxNdIy55uEhYaiS5fOQtvIuvqXf8qxbKHt1JaenQtZlvH0U2JvE89MFxOY6svLrXqGA4eNQJcbeglrh4GGmgMGGqrD0dEJUBX0efh5tGyr7TRBcfohpG5djw59b4V/R+1P8axmNBpw6CcZMwQcfibLMl58RtwZLrePGYvWPmJOhAWA7LRkHNi1DU/+0zJb52VZxhOvzLVIW7XdMXEyWvtrf9Hi8UP7cODHDZj8nOXXBZF94aJg0zHQUIOCeg+GT7j28+mpW9fjTJ64i/YA4MLpYiiKgkejXoVvoHYn1KYmJeDndatw3zOvwKedtiffph3Yi9jNX2Lb16s1rXst3W8fgxatxd6vdDonA1kJv+LpsFAEuLkKbava4ZLz+LqgEF36D0FgZzEXLSZ+vw5vTHsBYcFiLxEFgLjEJCxZ9aXwdkiHOOVkMgYashjFWLX7aPP/TbdIe76BwQgK1+4AvOL8qnVFPu2CERDWUbO6AHD6ZB5URcEb0ycjLEj7kYVqO/fuQ8wXG5C8zXILZwf6eKNLq1YWa+/rAvGH0Y0adiv69BQ3ylgbAw1R4zDQkMXIDk5QVQXBg++Bi6e4Q91K8zJw6ugefDBT+xNqJVnGp3Mna1632qhhQ9DnRjEjC9U+/XydkFGm+jIPJWLH158LbYOo2eLBeiZjoCGLUYxVC47zYi1z1krHv02CW5sAzeqdO34YJ3ZtxLxnxyK8nbb3+MQeSsbSr7YgNfOYpnXryz5RdWqriFGm+k6fFHuvElFzZo3LKXfu3ImFCxdi//79KCgowMaNG3H//feb1QdLYqChBp3Pz9K85qWzpwBVwcCxUfAUuMspPyUJh7esg1ubALTw125r8uUzBQCA0YP6an6BpMFQiU82/ozxUVM1rdsQ0aNMelCcKyYYnivkMf7UfF28eBG9evXCxIkTMWbMGGt3x2QMNFRH5dVTdndEvy6sDU+/QHgHtRdW/9zJXECScXilmPNuUrK1H3k4ceo0FEWBR7cRcGzRWvP61a6czsHF7H2IfHUWgoLFnkOTlBCPdatihLZRX6WiQJJlrHtHfDAkEkqWzZ8yMvHzo0ePxujRo81r04oYaKgOR8eqdS5ho5+Cm7d20zUAcO7YQRTt/QFb3pupaV1LkmUJE+a8J6a4JKE05Rcxteu5dfgd6Nazt/B21q2KQfbFS8LbqXbqSgVURcFjUdPhG6j94urUpET8tG6l5nWJ/kTDXU6lpaV13nZxcYGLi4t5tXWIgYYa1KbbzWgVrP3Ba4Xx3+EOR2+0lsX9Ty/HWI69xlLc/uRL8NbwLJLsIweQsHkN5r32KsJCtF1QG7cvAUtjxGwJr82SC3UNBgNkAG+kWOYQutr63DoMHbrfKKQ2Aw1ZgiQ7QDIz0FR/Pji47r9T5s6di3nz5plVW48YaKhBl4q1P+L98tmqNShdHN0RKIs9l2SvsRSt/YPQNkS7tS5nr66fGDViOCJ69dSsbrWlMauEL9atXqh7PEPbU5QbUlR4EgoA1w7DILl5CW8PAIwleag4kWiRtohsRV5eHjw8PGq+bo6jMwADDdVTffVByudvWbsrTWaECkmWsWH+NCH1U9MzNK+Zefy4xRbryrKM116wzEnBACC7eUFuIfYQv2rq5RKLtEMknKTBGhqp6vMeHh51Ak1zxUBjY05mZQqtX1JcCKgKfAaOhZOntluTL+cno+TwFk1rNsQBElRFgXuHW+Gg4chAxdkcGAqPYMKz2p9vU+2BZ1+Bj4anG9eXeTARv371Oea+Lv5yyp2xsYhZsQqXjmwU2g5Rc6TllJO9YKCxEUplJSRZxtLZL1mkPSdPP7h4a7uo0nC+6gTXYsWgad36TisGQJJw6dguIfXnzXgF4aHaho7YvYlYuuJz9Bw8DKFdtZ/Oqu3Xrz7HqJF3IKK3uMsOq8WsWInps+YgSHB4qpYYH4+VMcss0hZRc1NWVobMzD/+ozkrKwsHDx6Et7c3QkLE3PauJQYaGyE7OkJVFATf8SRcWms7clJb6fHfUbz/BxRsEbOTRwKwrqJISO36nps6C+003Jp8KDEeGz6LQVhIMLp20nadS1ZO1dqWAsEjcKfzLX/Y3bA7RuLGXr0t1h4DDTULVjgpODExEcOGDav5esqUKQCA8ePHY8WKFeb1xQIYaGxM6y790TJQ+91HtRUnfo/XXp+DEI3/q3pffDxiPl2GNyc/jfBgbW/yri12/+9YvPprtAsOQXhH7Z5VbtYxyLKMCc9P1qxmbZIsY9kcMbWJyMZY4RyaoUOHQlVV89q0IgYaatCIkSPRq3eE5nVjPl2GsKB26NohTPPa1dKzciHLMma/9IyQ+lqP/AB/jP7MmzVd8y3htcXtTcDSGG47JqLmh4HGxlwuzhVa/8q5qnUu6Wnanx1y/FgmZFnGE1PnaV67IfNmTtM0HMTtS8TS5asweNgd6NpT+/UnGz6LubolXMz5KdUYaIj0zxp3Odk6BhoboVRWApKMjHXzhbclyzKenfSUsPr+A+6Cs4e427YvFhzHmSOxGHX7MM3Pi1m6fJWm9awlVUBgrS8rR/uzjIjshoYnBdsLBhobITs6AqqCmx+LgoevuIsdC1KTcPindZg7dy7CNF5Ds2v3bsQsX4HCvT9oWvdaUjO0XWCbnVu1oDYrU/swcDLPMn/5GwwGyLKMJyeJmY4jIo0w0JiMgcbGhPW5Bb4dugtt4/BP63DnyJGIiBCwhiYmBnMfvwdhfuJGaHYdycCKn/dgwnMval5b5NocQMyhfbXlFxRCURQ88dJ0+AeJ3Up99MBefL+G01tEZBkMNNQgEVMS2VenIO7s2wMRHcWeaRCzJRZ+/UbDqZV2welS4XGcTY7Do1Gvwlfjw++OJsbj141rMOGZKE3rXku/20agUw+xa3UAMNAQNZEky5DM3OVk7udtDQONjTl7Ikto/XMnq3YITZw4UVgbqXmFwmoDQMbJYkCSUZTwo5D6fW4ZjvbdtT/8bttXX+CNV19CWLD2t0RXi0s4gCWr1girT0QakTSYcpI45UQ6pFRWQpJk/Pz+TIu0p/UOIQDYGRePFZ+vwcT/rNC07rU8FjUdvoHahYPUpET8tG4lTmRpPy1UfPXAu1HDb0Ofnj00r18bAw0RNUcMNFddOp0HBxc3a3fjmipKT0FVFTz8/DS0DRI3XZOelICt61cJ2SEEADGfrcZ4nwD4O4m77fXo5TJ8V3IafW4dhg7dtZtWqaw04OcNMj6YKe4uJyIiAIAk1VwuaVYNO8JAc1Xm+v+zdhcapfetwxHeTexdP1vXr9J8hxDwxy6h/i090cnVXfP6tX1Xclrzmo6OTlAVBT43j4Wjh5+mtcsLUixycScR2QhJ1iDQcA2NXXrupq7o6K3f69X3nzyN1UePI/+42F0wBTlXj/cXsEPI0k4c1zaUFeefAAC0CLsJrr4dNK0NAGCgISJqMrsNNNHR0YiOjobRaAQA3Brij37t2lq5V9dmMCpYk5KN6FmWCRqz58xFqMbn0MTHx+OTZUuRe6Vc07r15VeUQ5JlvD9TzI6hirMnNK9ZWWqZCzuJyDaokgzVzBEWcz9va+w20ERGRiIyMhKlpaXw9PS0dnf+kpODDFVRLHax48iRI9Fb43NoDBUViAHwbkG2pnWvZe7cuQgNC9Os3u5duxCzfDkKfhZzEzkRUQ1OOZnMbgONrRo9dBD69OgitI3Fq79GmoBzaE4WFEABLHYB48g779T8cMCYmBj0efh5tGyrbagsTj+E1K3rNa1JRGRPGGhsTEpmttD61TdV//MpcefQ2PoFjEG9B8MnvJvmdRloiKiGJJm/S8mGdjmVlJTAy8vLrBoMNDbCYFQselP1TaP/jlbevprWLMpKQ0rcNuHH+1fvpiIislmyXPUyt4YOvfvuuwgLC8MjjzwCAHj44Yfx1Vdfwd/fHz/88AN69erVpLoMNDbCyUGGoijw6nkHHN1bC2un/NRxXMw+gP0/bhBSX5Zlix3vT0Rkq5rzouDFixfjiy++AABs3boVW7duxY8//oh169Zh2rRp+Pnnn5tUt1GB5oMPPjC58JNPPolWrVqZ/DlqmMGoAJKEksNbLdLevFkzERaq7QF+cfH7sDRmOZ59ZRbaBYs7HPBQ4l589XkM0lJTNa2bk52taT0iIntUWFiI4OCqdZTfffcdHn74YYwcORJhYWEYMGBAk+s2KtBMnjwZQUFBcHBo3L0QeXl5+Nvf/sZAoyEnBxlQVbwx7QWEBQcKaycuMQlLVn2JUXfcjojeTRv2u56lMcsxaNgIdO2pfe1qBoMBG1evEHofFRGRUM14l1Pr1q2Rl5eH4OBgbNmyBW+99RYAQFXVmqNUmqLRU06JiYnw9W3cmgoGGXFGDbsVfXp2F9rGklVfClnnkp2TW/V/M9M1r13bqcICKIqCN1+ahPCgAM3qxh44jMVrNuJ8vvYXhJadOql5TSKyYc040IwZMwZjx45Fp06dcObMGYwePRoAkJSUhI4dOza5bqMCzdy5c9GyZctGF33ttdfg7e3d5E5Zw/FzF+DupN8lRXnnLwIAUjOPC20n83hO1TqXp58TUl+WZcyZ/IyQ2vWNHnKzplvcKyoNWLJWxo7o1zWrSURkb9577z2EhYUhLy8P//73v2vyRUFBAZ5//vkm1210oDHFzJmWuRFaS6//tt/aXfhLsixj/IszLNLW3AkPINxf25OT445mYOnmX/G6gFOIa4uPj8eny5ZqXtfZ0QmqqsAz4iE4ttT22VwpTkdZ2i+a1iQiG9aMR2j27NmDyZMnw9GxbgR54YUXEBcX1+S6jR6SKC4uvu6UU2VlJQ4cOID+/fs3uTPWNLaNP0JcXK3djWs6eqkMm0tOY+SDY+HdVtuLEWvLTk9B/PYtGNX/RkR00j50LN38q5BTiOsTEWiquQX1grNPuOZ1GWiIqJoqSRrsctLnOTTDhg1DQUHBnzLF+fPnMWzYsCavo2l0oAkICKjTgZ49e+KHH36oWal85swZDBw40KwFPdZ0U0sP3NhCv5dTVqoqvis9i5+/Wm3trhARETWZqqqQGghbZ86cQYsWLZpct9GBRlXVOl9nZ2fDYDBc92dIO46SBFVR0O/R59HKN0hYO4WpSUj+eT1Sc7VfpJpVeErzmkREzVIznHIaM2YMAECSJEyYMAEuLi413zMajfj9998xaNCgJtfXdBVsQ4mLtNXKNxCtA7Wf7qh2viAHsizhyQXLhLUh4p6o2nJycoTWJyISrhlefVB9EbSqqmjVqhXc3Nxqvufs7Iybb74ZkyZNanJ9/W7roToqVRWSJGP7B7Ms0t6kKTM1P/wuaW8cvlv3BSYJvCeKiIj0afny5QCAsLAwTJ061azppYY0OtBIkoQLFy7A1dW1Zv6rrKwMpaWlAFDzf23ViYorcJMvWrsb13S60gBVVfBc5/Zo18Ltrz/QRIfOlGBDbj4GDhuBLjdof/jdprWfIWTkRLh6+2teu1pp9hEUxm8SVp+ISLhmOOVUzdSd041l0hqazp071/k6otZOlWst8rEV7xXYxjTFYD8fdPUUu3h5Q26+kMPvTuZVHazn3bU/WgZ1/oufNg8DDRHZsuZ8l1NRURGmTp2Kbdu2obi4+E/rb4Xvcvr111+b1ICt6AMPeMLJ2t24piJcQTLKkHVB7ChSbtlFyLKMN18Wc7AeERE1gqTBbds6DTQTJkxAbm4uZs+ejYCAAM0GQxodaIYMGaJJg3oVCFe0g7ipHHMpUCEDmH3wqEXae/C5aWjbLljTmumHErF9wypcKs7VtG595WcLhdYnIqKm2717N3bt2oXevXtrWrfRgUZRFCxcuBCbNm1CRUUFbr/9dsydO7fOKmUSR4YEBcC9rj7wkcSt5c40lmNnRQl6DR6GsG49Na1dWWnA9q9kpK99R9O6RETNTjNeQxMcHCzkmJdG/8349ttvY968eRgxYgTc3Nzw/vvvo7i4GDExMZp3iq6tp2MLhDgKPNH4CrATwMmsTM1LnysuBFQFr70+ByECrz7YFx+PmE/FbTsnIhKuGQeaRYsWYcaMGViyZAnCwsI0q9voQLNq1Sr873//wzPPVF0s+Msvv+Duu+/GJ598AtnceT4dOA8DnKDf36MMlRZpx6iqkGQZi2e/KKyNESNHoldvsVcfMNAQEenTI488gkuXLqFDhw5wd3eHk1Pd9atnz55tUt1GB5rc3FzcddddNV+PGDECkiTh5MmTCAoSd3KtpezEOWt3oVEKlAqIzDbnYISqKHhu6iy0C9Z2FOVQYjw2fMYRPSKiv9TMR2hEaHSgqayshKtr3akOJyenP11/QOJIAGIuFVikrcHD7kDXntqfQ8NAQ0T015rz5ZTjx48XUtekc2jq371QXl6OZ599ts5pf19//bW2PbSQB9y8Eejg8tc/aCVplZextbwET7w4HX5B2p7gW1vygX34fu1KYfWJiIiOHTuG5cuX49ixY3j//ffh6+uLH3/8ESEhIejRo0eTajY60DSUqP7xj380qVE9utG5Jbo5a3sMs6YuA1tRAr+gEAS37ySsmaITYrdUExFRIzTjKacdO3Zg9OjRGDx4MHbu3Im3334bvr6+OHToED799FNs2LChSXUbHWiq72BorgqMFXA16PMPH6jqnyzLWPhqpLW7QkREojXDyymrzZgxA2+99RamTJmCVq1a1bw/fPhwfPTRR02uy8spr1paZhuHsc2d+HeEB7QVVj/ucDqWbtqGrEztb8Q+mWcb10sQEZE4hw8fxurVq//0vq+vL06fPt3kuo0KNGPGjMGKFSvg4dG4O4TGjRuH9957D76+vk3umKVNnzEDnTt3sXY3rik+fg+WLV2KUQN6IaJzuLB2Kior8cl3v2L2S88Ia4OIiP5CM55y8vLyQkFBAcLD6/5dlpSUhMDAwCbXbVSg+fbbb3Hq1KlGFVRVFZs3b8a//vUvmwo0wcHB6NJFv4EmJyfbIu04OzpCURTMe3YswttpeyN27KFkLP1qi6Y1iYiao+Z8OeWjjz6K6dOnY/369ZAkCYqiIDY2FlOnTsUTTzzR5LqNCjT1b9pujqIiuTalttGD+iKiawfN6zLQEBE1QjMeoXnnnXcQGRmJ4OBgGI1GdO/eHUajEWPHjsXrr7/e5LqNCjRNuWnbnGEja3jQvY2ut22nGi7h5/ISpOacFNpOVkHjRuKIiIiawtnZGcuWLcPs2bNx5MgRlJWVISIiAp06mbeDt1GBprnftA1UbdvuruNt2waokCvO48l3PrZ2V4iISLCqg/XM26Wk14P1du/ejVtuuQUhISEICdHuXDXucrrqpM63bZ9TKqEoKqKmz0JgkLiLHZMS47FuJU/zJSKyJlWteplbQ4+GDx+OwMBAPPbYY/jHP/6B7t27a1KXgeaqJRcsc6WAuW4Zdge639hbaBsMNEREJMrJkyexdu1arFmzBgsWLMCNN96IcePG4bHHHjPrbkgGmqteHtEPnf3aWLsb17Q3Kx8r4g5buxtERGQBiqpCMXOIxdzPi+Lj44OoqChERUUhKysLq1evxsqVKzFz5kzcdttt2L59e5PqMtBcNbRLGAZ20PdC5hVxh5GVmS60jfwTPPyOiMja1Ksvc2voXXh4OGbMmIFevXph9uzZ2LFjR5NrMdBclXnqHFq4OFm7G9d07NQ5yLKMmVFPW7srREREZouNjcUXX3yBDRs2oLy8HPfddx/mz5/f5HomBZoffvgBX3/9Nby9vTFx4kR07dq15nvnzp3Dgw8+2OShImt7Zd0v1u5CowxzbI3WEBe8ctXLSDReEFafiIj+mqJWvcytoUczZ87E2rVrcfLkSdxxxx14//33cd9998Hd3d2suo0ONKtXr8YTTzyBUaNGIS0tDR9++CE++eQTjBs3DgBQUVFh1lARNU4nuQUCZIHn5VQCiWCgISKyJlVVoZq5Bsbcz4uyc+dOTJs2DQ8//DB8fHw0q9voQLNw4UL897//xYsvvggAWLduHSZOnIjy8nI89dRTmnXIWsa28UeIi6u1u3FNRy+VYXNJ0y/tIiIi0oPY2FghdRsdaDIyMnDPPffUfP3www+jbdu2uPfee2EwGPDAAw8I6aCl3NTSAze2aNzlm9bCQENEZB+a85QTAHz22WdYvHgxsrKysGfPHoSGhmLRokUIDw/Hfffd16SajT5JzsPDA0VFRXXeGzZsGL777jtMmzYNH374YZM6QERERH+mmvnSq48//hhTpkzBXXfdhZKSEhiNRgBVt3AvWrSoyXUbPULTv39//Pjjj7j55pvrvD9kyBBs3rwZf/vb35rcCT04UXEFbvJFa3fjmgoNFdbuAhERWUhzHqH58MMPsWzZMtx///1YsGBBzft9+/bF1KlTm1y30YHm5ZdfRlxcXIPfGzp0KDZv3oxVq1Y1uSPW9l4Bz18hIiISLSsrCxEREX9638XFBRcvNn1godGBZsiQIde9pHLYsGEYNmxYkztibf/s3B7hHvq9nPLg6XNYn33C2t0gIiILaM67nMLDw3Hw4EGEhta9l3DLli3o1q1bk+s2OtBUVlbCaDTCxeWPLcNFRUVYvHgxLl68iHvvvRe33HJLkztibTf7eeOmtt7W7sZ1MdAQEdkH5erL3Bp6NGXKFERGRqK8vByqqmLfvn1Ys2YN5s+fj08++aTJdRsdaCZNmgRnZ2csWbIEAHDhwgX069cP5eXlCAgIwHvvvYdvv/0Wd911V5M7Q0RERM3bP//5T7i5ueH111/HpUuXMHbsWLRr1w7vv/8+Hn300SbXbXSgiY2NxUcffVTz9apVq2A0GpGRkQFPT09Mnz4dCxcuZKAhIiIyk6pWvcytoVfjxo3DuHHjcOnSJZSVlcHX19fsmo0ONPn5+ejUqVPN19u2bcODDz4IT09PAMD48eOxfPlysztkLTlll+DuqN+rrfIvXrZ2F4iIyEKa8y6n2tzd3c2+8qBao/8Gd3V1xeXLf/ylGh8fj4ULF9b5fllZmSadsoZ/JSVbuwtERETURI0ONL1798Znn32G+fPnY9euXSgqKsLw4cNrvn/s2DG0a9dOSCctwXfgI3DxCbZ2N67p0olknEn6wdrdICIiC2jOu5xEaXSgmTNnDkaPHo1169ahoKAAEyZMQEBAQM33N27ciMGDBwvppCU4e/rCVceBxlBS9Nc/REREzUJz3uUkiknn0Ozfvx8///wz/P398dBDD9X5fu/evdG/f3/NO2gpJ7bw6gYiIiJbZdIq2G7dul3z0Junn35akw5Zy4xpU9G1c2drd+Oa4uLjsfTTGGt3g4iILECFBrucNOmJGFFRUXjzzTfh7a3d+W8mb+uZP38+/Pz8MHHixDrvx8TE4NSpU5g+fbpmnbOkEcOH4ZbBg6zdjetioCEisg+KqkIxM9GY+3mtnThxAkFBQQCA1atX49VXX4W3tzd69uyJH374AcHB5i37aPRt29WWLFmCrl27/un9Hj16YPHixWZ1hoiIiMy/aVuPN2537doVoaGhGDt2LMrLy5GXlwcAyM7OhsFgMLu+yYGmsLCwzmLgam3btkVBQYHZHSIiIiLriI6ORlhYGFxdXTFgwADs27dPs9olJSVYv349brrpJiiKgrvuugudO3fGlStX8NNPP6GoyLzNLyYHmuDgYMTGxv7p/djYWJvetk1ERKQX1QfrmfsyxZdffokpU6Zg7ty5OHDgAHr16oU777wTxcXFmvxOBoMB/fv3xyuvvAI3NzckJSVh+fLlcHBwQExMDMLDw9GlS5cm1zd5Dc2kSZMwefJkGAyGmnNotm3bhldffRWvvPJKkzvSVCUlJRgxYgQqKytRWVmJl156CZMmTbJ4P4iIiDSjwdUHps45/fe//8WkSZPw5JNPAgAWL16M77//HjExMZgxY4aZnQG8vLzQu3dvDB48GBUVFbh8+TIGDx4MR0dHfPnllwgMDERCQkKT65scaKZNm4YzZ87g+eefR0VFBYCqU4KnT5+OmTNnNrkjTdWqVSvs3LkT7u7uuHjxIm644QaMGTMGbdq0MalOekYGWrRoIaiX5svKybF2F4iIyAaVlpbW+drFxQUuLi513quoqMD+/fvr/D0uyzJGjBiBPXv2aNKP/Px87NmzB3FxcaisrMRNN92Efv36oaKiAgcOHEBQUBBuueWWJtc3OdBIkoR3330Xs2fPRkpKCtzc3NCpU6c/PRxLcXBwqLkH4sqVK00+XfH5Fydr3DMiIqKmUaBCMXNZb/Xn6+8emjt3LubNm1fnvdOnT8NoNMLPz6/O+35+fkhNTTWrH9V8fHxwzz334J577sHixYuxc+dOpKSk4IknnsDUqVPx+OOPo3///tixY0eT6jf5NsaWLVuiX79+Tf14jZ07d2LhwoXYv38/CgoKsHHjRtx///11fiY6OhoLFy5EYWEhevXqhQ8//LDOIX4lJSUYMmQIMjIysHDhQvj4+Jjcj392bY/2HvodoUk6fQ7rj5+wdjeIiMgCtLxtOy8vDx4eHjXvW2sAoj5PT088/PDDeOqpp7B9+3a4u7s3OcwAZgQarVy8eBG9evXCxIkTMWbMmD99v3qR0uLFizFgwAAsWrQId955J9LS0mquG/fy8sKhQ4dQVFSEMWPG4O9///ufUuZfGejfBje1NW2aytIYaIiIyFQeHh51Ak1DfHx84ODg8KedRkVFRfD399e8T7///jsCAwMBAKGhoXBycoK/vz8eeeSRJtc0eZeT1kaPHo233noLDzzwQIPfr71IqXv37li8eDHc3d0RE/PnQ+b8/PzQq1cv7Nq165rtXblyBaWlpXVeREREemLpXU7Ozs646aabsG3btj/6oCjYtm0bBg4cqPnvFxwcDFmuiiBHjhwx+1A9QAeB5nqqFymNGDGi5r36i5SKiopw4cIFAMD58+exc+fO6277mj9/Pjw9PWteWjxEIiIiLVVPOZn7MsWUKVOwbNkyrFy5EikpKXjuuedw8eLFml1Pemf1KafracwipZycHDz99NM1i4FfeOEF9OzZ85o1Z86ciSlTptR8XVpaylBDRER275FHHsGpU6cwZ84cFBYWonfv3tiyZYvJSzisRdeBpjH69++PgwcPNvrnG9quRkREpCda7nIyRVRUFKKiosxq11p0HWgsvUiJiIhID7Tc5WQvdL2GxtKLlIiIiPSg+rZtc1/2xOojNGVlZcjMzKz5OisrCwcPHoS3tzdCQkIwZcoUjB8/Hn379kX//v2xaNEim1qkREREROJZPdAkJiZi2LBhNV9XL9gdP348VqxYYfOLlIiIiExlVKpe5tawJ1YPNEOHDv3LqwpseZESERGRqbSYMrK3KSddr6EhIiIiagyrj9AQERFRXYqqwsgRGpPYbaCJjo5GdHQ0jEajtbtCRERUR9XVBeYGGo06YyPsdsopMjISycnJSEhIsHZXiIiIyEx2O0JDRESkV9zlZDoGGiIiIp3hLifT2e2UExERETUfHKEhIiLSGaMGu5zM/bytYaAhIiLSGQXm71KysyU0DDRERER6Y1RUGM1MNOZ+3tZwDQ0RERHZPI7QEBER6YyqwS6nv7onsblhoCEiItIZo1r1MreGPbHbKafo6Gh0794d/fr1s3ZXiIiIyEx2G2h49QEREelV9cF65r7sCaeciIiIdIa7nExntyM0RERE1HxwhIaIiEhneJeT6RhoiIiIdIa7nEzHKSciIiKyeRyhISIi0hlOOZmOgYaIiEhnFEWFYuYuJXM/b2sYaIiIiHRG0WANjZ3lGa6hISIiItvHERoiIiKd4Roa09ltoImOjkZ0dDSMRqO1u0JERFSHUVVhNDOQmPt5W2O3U068y4mIiKj5sNsRGiIiIr3iLifTMdAQERHpjBEanBSsSU9sh91OOREREVHzwREaIiIineEuJ9Mx0BAREekMdzmZjlNOREREZPM4QkNERKQziqLCyF1OJmGgISIi0hmjBoHG3M/bGgYaIiIinWGgMR3X0BAREZHN4wgNERGRzhgV80dYjIpGnbERDDREREQ6wykn03HKiYiIiGye3Y7QREdHIzo6Gkajvd12QUREescRGtPZ7QhNZGQkkpOTkZCQYO2uEBER1VF9Do05L3s7h8ZuAw0RERE1H3Y75URERKRXRlWDKSc7u8uJgYaIiEhnuIbGdJxyIiIiIpvHERoiIiKd4QiN6RhoiIiIdKZSUeFgZiCpZKAhIiIia+IIjem4hoaIiIhsHkdoiIiIdEbRYITG3g7WY6AhIiLSGaOqmn2OjL2dQ8MpJyIiIrJ5HKEhIiLSGS4KNh0DDRERkc4w0JiOU05ERERk8zhCQ0REpDMcoTGd3Qaa6OhoREdHw2g0WrsrREREdRhVBUZFMbuGPbHbKafIyEgkJycjISHB2l0hIiIiM9ntCA0REZFe8WA90zHQEBER6YxRUSFzDY1JGGiIiIh0plIBJLNv29aoMzbCbtfQEBERUfPBERoiIiKd4ZST6ThCQ0REpDPV59CY+xLl7bffxqBBg+Du7g4vLy9h7ZiCgYaIiIhMUlFRgYceegjPPfectbtSg1NOREREOqP3Kac33ngDALBixQphbZiKgYaIiEhntDyHprS0tM77Li4ucHFxMau2HnHKiYiIqBkLDg6Gp6dnzWv+/PnW7pIQHKEhIiLSGaOimn0OTfUIT15eHjw8PGrev9bozIwZM/Duu+9et2ZKSgq6du1qVr9EYaAhIiLSGVVVoZoZaFS16vMeHh51As21vPLKK5gwYcJ1f6Z9+/Zm9UkkBhoiIiJC27Zt0bZtW2t3o8kYaIiIiHRGUVSzL5cUeTllbm4uzp49i9zcXBiNRhw8eBAA0LFjR7Rs2VJYu9fDQENERKQzqqrWTBmZU0OUOXPmYOXKlTVfR0REAAB+/fVXDB06VFi718NdTkRERDqjKqomL1FWrFhRE7pqv6wVZgAGGiIiImoGOOVERESkM3pfQ6NHDDREREQ6oypVL3Nr2BO7nXKKjo5G9+7d0a9fP2t3hYiIiMxkt4EmMjISycnJSEhIsHZXiIiI6mhowW1TXvaEU05EREQ6wzU0prPbERoiIiJqPjhCQ0REpDNanCMj8hwaPWKgISIi0hstDsazs0DDKSciIiKyeRyhISIi0hlFVSGZuUtJ4S4nIiIisiZV1WANDQMNERERWRMXBZuOa2iIiIjI5nGEhoiISGcUBZDMPlhPo87YCAYaIiIindHi6gJ7W0PDKSciIiKyeRyhISIi0hlVqXqZW8OeMNAQERHpjKKoGqyh4ZQTERERkU3hCA0REZHO8Bwa0zHQEBER6QwDjek45UREREQ2jyM0REREOsPLKU3HQENERKQznHIyHQMNERGRzvC2bdNxDQ0RERHZPI7QEBER6YyqqGYfjMcpJyIiIrIqXk5pOk45ERERkc2z2xGa6OhoREdHw2g0WrsrREREdXCXk+nsdoQmMjISycnJSEhIsHZXiIiI6lCurqEx92VP7DbQEBERUfNht1NOREREeqUqRqiKeUsizP28rWGgISIi0hkGGtNxyomIiIhsHkdoiIiIdEZVFA1GaBSNemMbGGiIiIh0RjUaoZp5rIi5n7c1DDREREQ6o6oarKFR7SvQcA0NERER2TyO0BAREekMdzmZjoGGiIhIZxhoTMcpJyIiIrJ5HKEhIiLSGY7QmI6BhoiISGd4Do3pOOVERERENo8jNERERDqjKEbAzBEahVNOREREZE1cQ2M6TjkRERGRzeMIDRERkc5whMZ0DDRERER6YzRClc0MJLyckoiIiKxJVc1fFMzLKYmIiIhsDEdoiIiIdEZVFPNHaOzsYD0GGiIiIp1RNTiHxt4WBXPKiYiIiGweR2iIiIh0pmrKybwpI045ERERkVVxysl0nHIiIiIim8cRGiIiIp3hCI3pGGiIiIh0RlGMkBhoTMIpJyIiIrJ5HKEhIiLSGdWoAJKZIzRG7nIiIiIiK+JdTqZjoCEiItIZVTGaP0LDNTREREREtoUjNERERDrDERrTMdAQERHpDAON6ew+0KiqCgBIK7lg5Z5c3/HSiwCAk0o5KlRxK9dPqRUAgAOpmSi7XK5p7ZTsPADAoYMHcfHiRU1r15aRng4AOHA0HWWXLmtWN+VYDgDgyplsKAZtn43h/EkAQNLvySi7eEnT2rWlZh4DAGQmH8LlS+L+DAAg91jVn8NhwX/etR27+md/LPkwygX8fvlZVc8v6XAyyi6J+3OqlpqZBQBIOyLuzys7MwMAcDorBYZy7X+nkvyq3yH5dAkuGSo1r1/bsZJSAEDxsWQhv8uprDQAf/y9IZTRALNbMRq06InNkFSL/MnoT3R0NKKjo3HlyhUcP37c2t0hIiIbcezYMbRv315I7fLycoSHh6OwsFCTev7+/sjKyoKrq6sm9fTMbgNNtZKSErRu3Rq5ubnw9PS0dndsSmlpKYKDg5GXlwcPDw9rd8em8NmZh8+v6fjsmu78+fMICQnBuXPn4OXlJayd8vJyVFRUaFLL2dnZLsIMwCknyHLVRi9PT0/+w91EHh4efHZNxGdnHj6/puOza7rqvzdEcXV1tZsQoiVu2yYiIiKbx0BDRERENs/uA42Liwvmzp0LFxcXa3fF5vDZNR2fnXn4/JqOz67p+Oz0ze4XBRMREZHts/sRGiIiIrJ9DDRERERk8xhoiIiIyOYx0BAREZHNs+tAEx0djbCwMLi6umLAgAHYt2+ftbukO/Pnz0e/fv3QqlUr+Pr64v7770daWlqdnykvL0dkZCTatGmDli1b4sEHH0RRUZGVeqxfCxYsgCRJmDx5cs17fHbXl5+fj3/84x9o06YN3Nzc0LNnTyQmJtZ8X1VVzJkzBwEBAXBzc8OIESOQkZFhxR7rg9FoxOzZsxEeHg43Nzd06NAB//rXv+rcQcRn94edO3finnvuQbt27SBJEr755ps632/Mszp79izGjRsHDw8PeHl54amnnkJZWZkFfwuCaqfWrl2rOjs7qzExMerRo0fVSZMmqV5eXmpRUZG1u6Yrd955p7p8+XL1yJEj6sGDB9W77rpLDQkJUcvKymp+5tlnn1WDg4PVbdu2qYmJierNN9+sDho0yIq91p99+/apYWFh6o033qi+9NJLNe/z2V3b2bNn1dDQUHXChAnq3r171ePHj6s//fSTmpmZWfMzCxYsUD09PdVvvvlGPXTokHrvvfeq4eHh6uXLl63Yc+t7++231TZt2qjfffedmpWVpa5fv15t2bKl+v7779f8DJ/dH3744Qd11qxZ6tdff60CUDdu3Fjn+415VqNGjVJ79eqlxsfHq7t27VI7duyoPvbYYxb+Teyb3Qaa/v37q5GRkTVfG41GtV27dur8+fOt2Cv9Ky4uVgGoO3bsUFVVVUtKSlQnJyd1/fr1NT+TkpKiAlD37NljrW7qyoULF9ROnTqpW7duVYcMGVITaPjsrm/69OnqLbfccs3vK4qi+vv7qwsXLqx5r6SkRHVxcVHXrFljiS7q1t13361OnDixzntjxoxRx40bp6oqn9311A80jXlWycnJKgA1ISGh5md+/PFHVZIkNT8/32J9t3d2OeVUUVGB/fv3Y8SIETXvybKMESNGYM+ePVbsmf6dP38eAODt7Q0A2L9/PwwGQ51n2bVrV4SEhPBZXhUZGYm77767zjMC+Oz+yqZNm9C3b1889NBD8PX1RUREBJYtW1bz/aysLBQWFtZ5fp6enhgwYIDdP79BgwZh27ZtSE9PBwAcOnQIu3fvxujRowHw2ZmiMc9qz5498PLyQt++fWt+ZsSIEZBlGXv37rV4n+2VXV5Oefr0aRiNRvj5+dV538/PD6mpqVbqlf4pioLJkydj8ODBuOGGGwAAhYWFcHZ2/tPNs35+figsLLRCL/Vl7dq1OHDgABISEv70PT676zt+/Dg+/vhjTJkyBa+99hoSEhLw4osvwtnZGePHj695Rg39c2zvz2/GjBkoLS1F165d4eDgAKPRiLfffhvjxo0DAD47EzTmWRUWFsLX17fO9x0dHeHt7c3naUF2GWioaSIjI3HkyBHs3r3b2l2xCXl5eXjppZewdetW3pzbBIqioG/fvnjnnXcAABEREThy5AgWL16M8ePHW7l3+rZu3Tp88cUXWL16NXr06IGDBw9i8uTJaNeuHZ8dNVt2OeXk4+MDBweHP+0mKSoqgr+/v5V6pW9RUVH47rvv8OuvvyIoKKjmfX9/f1RUVKCkpKTOz/NZVk0pFRcXo0+fPnB0dISjoyN27NiBDz74AI6OjvDz8+Ozu46AgAB07969znvdunVDbm4uANQ8I/5z/GfTpk3DjBkz8Oijj6Jnz554/PHH8fLLL2P+/PkA+OxM0Zhn5e/vj+Li4jrfr6ysxNmzZ/k8LcguA42zszNuuukmbNu2reY9RVGwbds2DBw40Io90x9VVREVFYWNGzdi+/btCA8Pr/P9m266CU5OTnWeZVpaGnJzc+3+Wd5+++04fPgwDh48WPPq27cvxo0bV/P/89ld2+DBg/90REB6ejpCQ0MBAOHh4fD396/z/EpLS7F37167f36XLl2CLNf917uDgwMURQHAZ2eKxjyrgQMHoqSkBPv376/5me3bt0NRFAwYMMDifbZb1l6VbC1r165VXVxc1BUrVqjJycnq008/rXp5eamFhYXW7pquPPfcc6qnp6f622+/qQUFBTWvS5cu1fzMs88+q4aEhKjbt29XExMT1YEDB6oDBw60Yq/1q/YuJ1Xls7ueffv2qY6Ojurbb7+tZmRkqF988YXq7u6ufv755zU/s2DBAtXLy0v99ttv1d9//12977777HbrcW3jx49XAwMDa7Ztf/3116qPj4/66quv1vwMn90fLly4oCYlJalJSUkqAPW///2vmpSUpObk5Kiq2rhnNWrUKDUiIkLdu3evunv3brVTp07ctm1hdhtoVFVVP/zwQzUkJER1dnZW+/fvr8bHx1u7S7oDoMHX8uXLa37m8uXL6vPPP6+2bt1adXd3Vx944AG1oKDAep3WsfqBhs/u+jZv3qzecMMNqouLi9q1a1d16dKldb6vKIo6e/Zs1c/PT3VxcVFvv/12NS0tzUq91Y/S0lL1pZdeUkNCQlRXV1e1ffv26qxZs9QrV67U/Ayf3R9+/fXXBv89N378eFVVG/eszpw5oz722GNqy5YtVQ8PD/XJJ59UL1y4YIXfxn5Jqlrr6EgiIiIiG2SXa2iIiIioeWGgISIiIpvHQENEREQ2j4GGiIiIbB4DDREREdk8BhoiIiKyeQw0REREZPMYaIiIiMjmMdAQ0TVlZ2dDkiRIkoTevXtrWnvFihU1tSdPnqxpbSKyPww0RM3Ejh07MHz4cHh7e8Pd3R2dOnXC+PHjUVFRAQAoLy/HhAkT0LNnTzg6OuL+++9vdO1ffvmlzuV88+bNqwkjjo6OCAsLw8svv4yysjIAdYNQ/Vd8fDwA4JFHHkFBQQEvQyQiTThauwNEZL7k5GSMGjUKL7zwAj744AO4ubkhIyMDX331FYxGIwDAaDTCzc0NL774Ir766iuT6rdp0wZt2rSp816PHj3wyy+/oLKyErGxsZg4cSIuXbqEJUuW1PzML7/8gh49evypFgC4ubnBzc0Nzs7OTfmViYjqYKAhshGKouA///kPli5diry8PPj5+eGZZ57BrFmz8PPPP8Pf3x///ve/a36+Q4cOGDVqVM3XLVq0wMcffwwAiI2NRUlJiVn9cXR0hL+/P4Cq0ZZt27Zh06ZNdQJNmzZtan6GiEgkTjkR2YiZM2diwYIFmD17NpKTk7F69Wr4+fkBAPz9/VFQUICdO3darX9ubm4101tERJbGERoiG3DhwgW8//77+OijjzB+/HgAVSMwt9xyCwDgoYcewk8//YQhQ4bA398fN998M26//XY88cQT8PDwEN6//fv3Y/Xq1Rg+fHid9wcNGgRZrvvfTdXrbIiItMQRGiIbkJKSgitXruD2229v8PsODg5Yvnw5Tpw4gX//+98IDAzEO++8gx49eqCgoEBInw4fPoyWLVvCzc0N/fv3x8CBA/HRRx/V+Zkvv/wSBw8erPMiIhKBgYbIBri5uTXq5wIDA/H444/jo48+wtGjR1FeXo7FixcL6VOXLl1w8OBBpKSk4PLly9i0aVPNFFi14OBgdOzYsc6LiEgEBhoiG9CpUye4ubnV2Tr9V1q3bo2AgABcvHhRSJ+cnZ3RsWNHhIWFcacSEVkd19AQ2QBXV1dMnz4dr776KpydnTF48GCcOnUKR48exVNPPYUlS5bg4MGDeOCBB9ChQweUl5dj1apVOHr0KD788MOaOsnJyaioqMDZs2dx4cKFmikgrQ/Nq3bmzBkUFhbWec/Lywuurq5C2iMi+8VAQ2QjZs+eDUdHR8yZMwcnT55EQEAAnn32WQBA//79sXv3bjz77LM4efIkWrZsiR49euCbb77BkCFDamrcddddyMnJqfk6IiICAKCqqpA+jxgx4k/vrVmzBo8++qiQ9ojIfkmqqH+TEZHNy87ORnh4OJKSkoSN4gwdOhS9e/fGokWLhNQnIvvANTRE9JcGDRqEQYMGaVrziy++QMuWLbFr1y5N6xKRfeIIDRFdU2VlJbKzswEALi4uCA4O1qz2hQsXUFRUBKBqXY2Pj49mtYnI/jDQEBERkc3jlBMRERHZPAYaIiIisnkMNERERGTzGGiIiIjI5jHQEBERkc1joCEiIiKbx0BDRERENo+BhoiIiGze/wPYVWlTEbqcsQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "key = apt.randgen.get_key()\n", "key, model_hist = (\n", " tree.likelihoods[\"rn220_llh\"].components[\"rn220_er\"].simulate_hist(key, int(1e6), parameters)\n", ")\n", "model_hist += tree.likelihoods[\"rn220_llh\"].components[\"rn220_ac\"].simulate_hist(parameters)\n", "data_hist = tree.likelihoods[\"rn220_llh\"].data_hist\n", "\n", "bins = tree.likelihoods[\"rn220_llh\"].components[\"rn220_er\"].bins\n", "\n", "apt.plot_irreg_histogram_2d(*bins, (model_hist - data_hist) / jnp.sqrt(data_hist))\n", "plt.yscale(\"log\")\n", "plt.xlabel(\"cS1 [PE]\")\n", "plt.ylabel(\"cS2 [PE]\")\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 9, "id": "28f95026-9be1-40c3-89b2-1de67d1996cc", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "p = 1.0\n" ] } ], "source": [ "from scipy.stats import chi2\n", "\n", "t = jnp.sum(((model_hist - data_hist) / jnp.sqrt(data_hist)) ** 2)\n", "p = 1 - chi2.cdf(t, 30 * 30)\n", "print(f\"p = {p}\")" ] }, { "cell_type": "code", "execution_count": 10, "id": "a127627e-d0f4-43a4-98af-7206519d07a4", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkcAAAG4CAYAAABPb0OmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABKoklEQVR4nO3de3RddZ3///e5n9zTpG3SQks7yMUitkgvFPDHrcplBi84js4XtYALR40KdFRABcY1Ijd1EMkIzvoKyxkY0RkBccQRi1JFLqVaRymX6hSo0KSXNPecnNv+/eGXuNP3q7rbJjlp8nys1QX5ZF8+Z++Tk3f25/35vGNBEAQGAAAAMzOLV7oDAAAAkwnBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQAjBEQAAQMiUCY4GBwftsMMOs49//OOV7goAADiITZng6Nprr7UTTjih0t0AAAAHuSkRHG3evNmeffZZO/vssyvdFQAAcJCreHC0bt06O/fcc23u3LkWi8Xsvvvuc9u0t7fbggULLJvN2ooVK+zJJ58c9f2Pf/zjdt11101QjwEAwFSWrHQHBgYGbPHixXbRRRfZeeed575/zz332Jo1a+y2226zFStW2M0332xnnnmmPffcczZ79my7//777cgjj7QjjzzSfv7zn//Z8w0PD9vw8PDI1+Vy2bq6uqy5udlisdiYvjYAADA+giCwvr4+mzt3rsXjY/ysJ5hEzCy49957R7UtX748aGtrG/m6VCoFc+fODa677rogCILgiiuuCA499NDgsMMOC5qbm4P6+vrgs5/97F7Pcc011wRmxj/+8Y9//OMf/6bAv61bt455PBL7f0HJpBCLxezee++1t73tbWZmls/nrbq62v7jP/5jpM3MbPXq1dbd3W3333//qP3vvPNO+81vfmNf+MIX9nqOPZ8c9fT02Pz5823r1q1WX18/pq8HAACMj97eXps3b551d3dbQ0PDmB674sNqf8rOnTutVCpZS0vLqPaWlhZ79tln9+uYmUzGMpmMa6+vryc4AgDgIDMeKTGTOjjaVxdccEGluwAAAA5yFZ+t9qfMnDnTEomEdXZ2jmrv7Oy01tbWCvUKAABMZZM6OEqn03b88cfb2rVrR9rK5bKtXbvWVq5cWcGeAQCAqariw2r9/f3229/+duTrLVu22MaNG62pqcnmz59va9assdWrV9vSpUtt+fLldvPNN9vAwIBdeOGFFew1AACYqioeHD311FN22mmnjXy9Zs0aM/vDjLQ777zT3vWud9mOHTvs6quvto6ODluyZIn94Ac/cEnaAAAAY2FSTeWvhN7eXmtoaLCenh5mqwEAcJAYz9/fkzrnaDy1t7fbokWLbNmyZZXuCgAAmER4csSTIwAADjo8OQIAAJggBEcAAAAhBEcAAAAhBEcAAAAhBEcAAAAhBEcAAAAhBEcAAAAh0zY4YhFIAACgsAgki0ACADAh3hR/p2t7qPzt/ToWi0ACAABMEIIjAACAEIIjAACAkGSlOwAAAA5uUXOJ9je/aKLx5AgAACCE4AgAACCEYTUAwJgZy6naGBvqnkQV9d6p7aKedzK+P3hyBAAAEEJwBAAAEDJtgyPKhwAAAIXyIZQPAQDsgdypPziQvCG1b/KIwyMdr7j5d3/2HJQPAQAAmCAERwAAACEERwAAACGscwQAGFcHY/7OgfTvQNYV2hcqf0fl6kTdLvI5mmtd29lHXbHf++Zasq4tK/r3xrd/YfSxCjm3zVjhyREAAEAIwREAAEAIw2oAMI4OxiGlqCaiErsarokyzdts7MtmqOMV3rzUtaXE8dQQ09DCGZH6UrVlt2xX16H/3Se4tlI65tqa1fHEkNdQvX81ffPSri3XVOPaUoP+HPkG35Yc8G19bSe6tv55o78u5xJm/+X3HQs8OQIAAAghOAIAAAghOAIAAAiZtuVD2tvbrb293Uqlkj3//POUDwGmIfKBooua+xN52vgJi13TQz//TKTzRvXgc9e7tpX/54uurXHDDtfWv2ima6vdtDPSdomhkmsrVSX22s+wwZl+u+qd/nhmZj0L/LYNL/hthxv8c5DcDJ+HlN3tw4F8nd9O5Q0Vq0RbtT9eIuePV6gv++2G/XbBHk3lXM5e/NSnKR8yltra2mzTpk22fv36SncFAABMItM2OAIAAFAIjgAAAEKmbc7Rq3p7e62hoYGcI2AKmSq5RDLPR6xFo3J11DUY63ITUbdTawFFpXJ1VA5Nw/N+UZ2uY/zaO4FY3S/T43Ne1NpAisrdKfrTyrV81HYF0Rbz3fvDtg3iG+I3eqws8nfiakNxuPpCpO2s6O9Jqjbv2sol8UwmJnKTkv61FYdHvxfKgzl78eJ/JOcIAABgvBEcAQAAhBAcAQAAhFBbDQD2Q9S8JpU3pNbfedOJn3NtBVF7S9XZOvWsG/x2Ih9Inle8DrV2T1bkOu06stq1lf6/FteWKPicEpWrkxdpI9ku35bI+baXT/XJOrWv+PPuPtyfd3DYPycoi9+OKm9oeKY/RzkVLZU3aPT5PDGRfxOInCEzs1RVMdJ5VJ5PJuPPnc/7F91QN+Ta0gm/lpLqd1y0peJ+31lV/sKWRWJTbz4z6uviwLC96LYaGzw5AgAACCE4AgAACGEqP1P5gQk1mabZRy1NoabPR9W/wNdVqH3BD1Wo7aJS085LKd+WaxL7Zn3bjM1+GnXvPP+3dCAqYmT9qJ/1HOF/zcSKvn9xMWu8lIk4RJUWw1tZMdU9I9rEsJraLl3jp6bHxJT4YsFfmLpaPxZYmxl2bcGeNTLMrCzazMyqU74/pbKYUi+GwZJx//pmZvzwVlIMg+VKKb9v2u/bWxRvLqEm6a/DEVXbXdu2PeqWDPcX7JaT72cqPwAAwHgjOAIAAAghOAIAAAiZtlP529vbrb293UolP54KYN9FzSUa6+nuyV39rm27mE4+e12na4tKnWPXytn7fbwdx/tp5/k6v53K3+k/1LclfQqTlUU+UFQ7lvi/m4s1e6lhsYfBw0R+Ua1PJlLT0xNp/3msppyrvJy4yP3JpPxUd5Xns6PX55S1NvS6toK4qElR26Mq5ftcn/LnzZX8r+BZWf9eK+0l56gq4c+TikX7nZaJ+2vTIN5I8b3VLtlDXdznVKmp/Eo65vuSifnX1pQcfW2GytGWMtgf0/bJUVtbm23atMnWr19f6a4AAIBJZNoGRwAAAArBEQAAQMi0zTkCJjOVv6MkRYmI4ubf7fd51fFU7o8qV/ETkUu08v980bU1btjh2lSf1b6m1gzyW8lyFSoPKRCfgKpshDLY6tuqO3xb/7xoxwtEzsywr+IhS1MMt/o8k0S1z8cod6ddW6xO5QP5v5uTWb9dJuPPUZ32a++oUhBqPR6Vq5MQ6/HMyPjcmKbMoGvrzvv3S1qs23N4/U7XpnJyZqV9PlBC5OSUAv/a1Jo/s9N9rq1W1EbZWx5RQSw0lRW5OkpWLCpVF/fXNRf490xjwq9pVBA/TNmYfy+o15IL/LpJsxL+Wg+WR/dlIDl+OcM8OQIAAAghOAIAAAghOAIAAAihthq11TAGoq7xcyB1xaLmIRXevNS1VW3xC+YMLZyx39sNN/n8AlXfS6nuEGuz9Pq2rmP8WkBqraLO030uUc532fKNkbpnCb8cjQ3P97kTNujzPWIlfw2CBv/a4glR16rJ51js7PIZVemsz4WpqfKdVjk99Vmfz5JJRFsrRtXiUjW/Wqt8Hk1R5ODUJPw1rU/5/ilqDZxhkbfSnPJ9Uarj0XJjhsv+HC2pHtemcoFUPk9CrAOUEmv+1IjcHTOzksjlUtT+edHHdMQ1kqpFvtKguDZZ8VrqRC5XTrw/siqXa4/X299Xtjcs6qS2GgAAwHgjOAIAAAhhWI1htYPCRAxbKaqsxVhPlS82+6GTQr1/RJ364VP+gCcsjnTefjEFXhlu8H8vqantpYxo8zN+rcrP2rfev/BtYvay5ZvER5MYhhBP5C1Iin3F1HYlJqbUV9X6Yat0SkyfF0NPaugkX/JDGgsbuyJtV530Qxp5UdYiK4bLqpP+ddQn/cVX085bMn6ISg2rzUj5ad4NCT+kpEpGKGpafJ16wwh7Tv0202Uu1NR0Re2rho7UsJrSGPf3o1dMnU+YLuGRFcNgCfPvNzX8Vifeq6WI0UBCjOZVx/xrHgx8/+pi/gOlIF5fQYQm+T3a+vrK9vpF2xlWAwAAGG8ERwAAACEERwAAACGUD8FBIWrekMrpiepA8ppyb10R6Ry1m3yZglyLz+9IDPmx+l1/d2Kkc6hp7KJKgczzUSUxCjV+7D9W9kkH5azPG8iJ45Wr/Hb5hD9HPOOvQTLt29S0eDXBubHG570MDPv8jtb6XtdWFNPiD63xU7hzJVEGIRstV0eVq2gW+TvDIgmsWkyLb0j4chqqTEOTyLdR2zWK46nyEKqMRJ3IOVLUlHOVL9OopsWLXJvuss+1axKvo6CWGhC5RBmRP6aygVJiqYGEeB3q6UST2Ffl3+xtf9WfuLg2DXH/ubOz7D8UVC5RKuJzlbqYf805kYdULfKQUuJa1+5xPLW8xFjhyREAAEAIwREAAEDItA2O2tvbbdGiRbZs2bJKdwUAAEwirHPEOkcHbKzXFooq6hpEB5I3pHKE1LpEqtRFesCPhw81+79HZm3wOR+dK/zxxHI01rvQt4lV961UJdYHEnk+Qcq3pZt9fkd+yOekJFLRxv9nzfA5ODOqRH5M0echNGVFX8RaQI1pv50qYTEz7a99XFzAsti3PunPoXJ/VBkKlQ+k1sbR2/nr0pz0pUdUeYgasa5OVuS4qLaUemMJKp9KrQWk1uiJi2QxlUuUFbksSvScHHHeiOcoiV+hdXH/81EIol2/qpjft2jRynrsjepjddzn25VEH9W51fEyIm9oOPD3XV3XuLgDBZGbtOd6SH19ZVt49DbWOQIAABhvBEcAAAAhBEcAAAAhrHOEAzbWNc6i7qtyf2yzON6Jn3NtKr8o2+mTenatnO3aSik/Zp5v8OfNizplarsXzvX5RYXGaDkGQVrU8qqKVrOqusbnx8TF2iKppO9LMeNzUhbO8LXB0nG/b1muQuS9ptbneykZUUOsKu5fW1NSrefjPwLnprojnVfn5fi+qJwjpTHu84vSMi/H3yNVe0vl6jSJa1UQWacFkTdULWrOKSlxf1WeiTqvWlOnJF6HWmenLLZTeTCDgb9vtTGff9Mb+PysjIn6YRbt5y0e8X2fEn0u73WdI3/M/kDkvYnrNSyug8r9UesSZcV9UtdL9U+/38TaTvK+jz5e1Dyu/cGTIwAAgBCCIwAAgBCm8jOV38z0UFZUqmTHg89d79pOPesG15bqFY9T6/1U1qh6F/rH48kh/xYfmhlxmq5fYV/OBR6cJ4Y/Gvzj7fKAmAJf77cLRHmORNKfOCuGt+qr/PDgUMGf96im7a4tKaZr16X88fKihEVjyk9trxU1SmaI4S019DRY9vdSlbCoi/tzZOP+upTEUFGzLJ0hynOIshb1Ylp8Xgy71IjhN1USQw2D1YnyCIPivRF1yCsb838Pq6EPNQyz5zTqvVFlMlLivFGng6up5ElxnZWoU8kzYvr8YNnfc0UdT702tZ0asopamsNMDyWqYUh1TwbFudW+SklslhAfqTnxM6fkxFIZqqRLaY9r099XttNft5Wp/AAAAOON4AgAACCE4AgAACCEqfxTXNRp8VGn40elSnuk1NR7oW+ezzURaSBS/6G+LRCD4cVanz8RJMVAesnvG2sU08Rn+NwVVUbBGnxeTk1aHC/rc2v6C/66zK3p9ccTJSzqRd6QsiCzw7XtOc5vpvN8lFlJ37+ukn8fqO1SItdETZ9PiDwpVepCTXfX26lyFX479ZelKomh8m2UhDyib8smVO6Pf6+p6ekqv0jlrZTFtcqYP57aV5HHE/lFippefiB5SCofqF9MQ1c5ViofaLAcLa+pUI6WH2Titak8MzNdbkWcxlJi+QeV56PeRwOBz8dKRyxnonLr+so+kTMurnVZXOvePfYdLB5YWZU/hSdHAAAAIQRHAAAAIQRHAAAAIeQcjYMDKZ0x1tQaRAeSS6RexzGX/5Nrm/Urn/cy3OTfbgMtYln7Jn/eZL9vyzf6tmKdyA2ZKZa1F2PwGbFmUDrlx7SrUn67xqzPJcqXfA5EY8bn6jSnfb5STdL3eWbKXwS1PlBc5McoKn9nVrLPtdUl/GsriCQwVf6iWpTxOCTZLfqi1lzx+Qp1YruoBQSqRR6Iyg1RUuKj8kBKU0Rdf0eJi1y2lCjnoNoKYm0bnYfkqT5HzTlS12pIXKuoa/ToNZd8/wbKYm0wcZ37VJ6O6ov42VLlMPrEz756P6v3fZ9Y32tvJWhKIm+or5yR2+5J5dapPhYC9Vr8eQfEecuif6rEjvo8UcfrKo0utZTLF81sk9tuLPDkCAAAIITgCAAAIIRhtXFQqSE0pbj5d9E2PGGxa+o5stq1veFDfghNPcTdvtQ/GlZPwnOzxKPdBlHlWjy6jqX8AZNZv291lR/aqcv6Yau0qFiuhsHmV/vq8wMlfxXmZ/126rH8TDGUlRPTZ9Uw2NzUbteWFo+t4+IC1ovp+Gq4MSdKhajlAuri/rxZWUFeDQupISW/b1xM11bHU9Twka4W7+9RbdzfX1XVXE31VtPOo5bJUNRQlhouU8Ng6hoMi/Oq7frKfru6uL8fw+L69alhRDlt3DXJ4Z8BMdSjtlN///eUVVkKUQJE7KuW59hzevnepMXQWLe8zlW+L3spw5EL/OdsVpS6UXYUfbkN9ZqHy/tfzkkNyam2wZL4fSE0JEcP8athu7HCkyMAAIAQgiMAAIAQgiMAAICQaZtz1N7ebu3t7VYqRV9+vFJT9FUpDpVLVHjzUteWiphLlJvhx7RFConlG0SbH7q2Yp24rgk1X1ZMqa/3+UCZlO9MVdrnfMyr73ZtafFCZqT9tPOM2i7pt6tL+FydhoTfTuUXKXOTPm+oIHJr6mL+vOocs0T/VEmB6ni0fKBSwt8PtV1GlMlQ+TEZ8dqiTvWOWopDlZxQJSz0vtFymFQZClmKQ2yXlNPsxTRx0WfVP7VdThyvJN4H6qqoq6y261QHFP3LianaKvdETVkfEFPbVWmJAZF/UyPyb7rL/vNPTS/PRcy1UdPQu0v+HHpf9T6I/ms5ak6P0lfyuVJDJf+a1fGqE/5nc1fev+Yqsd1A0efvFUXukMrvqkmOvp/5wWj5Vftj2j45amtrs02bNtn69esr3RUAADCJTNvgCAAAQCE4AgAACJm2OUf740Dyiw6kZIeiyoJ0zfNj7vkGsRS9r1ZhQ7N9m1jKxnKHiUYxNpyo9mPNTfX+xCqXqDnr83ca076ERSru8xOaxItT+QnzMrtcm1pHKBsXpUJEflFd3PdPUTkQam2h5rgoeSKOlxDpMSoNpDHhf9Sj5q6oPJpqkSOkqLwctbZQQuQSZcTHk1ynR63nI65LJubzKXrKPj9L2V3291ddg66Sf7/UxP1r6y7546XFdcmJ66fW5FE5ZXmR21UQuR0qHyglctm6RB6NKvGwt1IXvi/R8m3yEbcbFLlJat2eqH1R1Lo6PUW/LlGvaFNlfdR5Ve5jV77GtZmZ5csiN1GUHvpdX7Nrq0/77YpiDahk3Pd7d86/vjpxvJ1Dvt+q1NL/7vT9q8n6z8razOhzFAf8OccKT44AAABCCI4AAABCCI4AAABCyDmaICpfSa1fFNXL57S4NrUuUaHOtw22+gSF4gy/cywtapdlxFpAdT4Hp6nat83M+nwgtQbRrHS/a9uzpo6ZWXXcj0lXi1wdZVay17UlRO7PgqTPTVJUjkZKHC+l1sERx6sTeSoq3yYr1suJSuUDqeOp9XeiHi+qwbK/lyr/KTLRlUGxbpLKp1LrAxXE8XZErBdWEAlBg6Jencol6hN1u1RuXLeox6XqbsXlekg+F0vl1qj8IpXno9YB6yrWujalp6Tqivm+DInzqtemqDyfvqK/zj1535esWLcnL+6lygVS6/aonBxV4zFX1GsuJUQ+ULHkz51O+mN2DPhfDoHo447dfrt4Qqy/JdpKRX/vtg3OcG1JkaPa0+tz3LqGR/elPBQtZ3B/8OQIAAAghOAIAAAghGG1CXLqWTe4tpQoAdLzvpWurZTyjzrFTFEbbhLTfuv8o85YjZiyLh5rVmWiledoykQrxTE36/dV035rxWP5uSlfYqNeTJ9Xx1PT57OiXEBW7JsRQ2NqynVJTnff/yGvlCx/4c+hymSo4aOo50iK6d9xMc1enUMOKamhDjFEpYbkCmpDISXO26nmtkekpokPiKEnNRzVJ4a31DIPaohKHU/1Jeq0c3WOnWKcPWoZD1WqRpXYWN+7MFL/6pP+51xNd9+Z90Nyg0U/rJYU+/YW/DUYEkNUasirKunf44WSHwYbyPu+DAz7NjXcpYbFYuIzJzfkj2dmlhJpDmporDgshvmK4n1UEmtgDPuf/2JK9DsnPjvEqg5xcY7ki2IZBvGSy3v8PosNjV8Iw5MjAACAEIIjAACAEIIjAACAEHKO9kKV+1DT8aOWBSm9dYVr67rkRNem0gnU1PtSvRjMzfq2+hk+H6g6LUociKXf51b76e6NKZ8/US+m2c8X5TnUNHs1fX5Wwp+3Ruyrps+nzV8DtV212lfkEmVE3pDKrVEyorxEWeTgRC6xIc6rpvdHzVcaFteqLPrSJ6asqxIlaiK1mp5eEDkRUfNoFJUPpKZ/q5yewcDnpKhcnazIW+sW+Sc9osRGSuS3qb4kxD3alm9wbaocTklc02E5xdy3DZV8X+JyyQl/jt68nwKv8neK4n68WPRTutU5SqKkRVe/v86NNf5zqHvA54Cp/qk8HZkbVxLXfthf05iY1h7kxXs8IU4yJHJ3hvVnznCVmOJfjPb5lOwT5xHpiuKj0kqZaPsmRakq9UhGbZdv9G11T4/euZSP20viFGOBJ0cAAAAhBEcAAAAhBEcAAAAh5BztRdT8opzIJRqc6cdjReqAiaVErH+BGOBt9PkOGbEuUX2VXzekscqPw8+u8uU5VBmPQ8W6RE1icFjlVKh1iVTekCqFUCfKgmTVwLeQFTlMdSL3R/1doNYqaoj7nIphscaPWm+oEIicHnEOlZuk1hGqjvn8giGxnVpbSOU19YhLmhH3sk/kA6VlGQq/ncoHUrk1LxcbXZtaa6er5Ne8Ue+hqKUu1DlUPlBP0ee4qON15X0eUlPa/8yoPJ9eUcJC5cvtGPbXYHfO59ao9XIKoqyFyrdRJSx6xDkKYq0c1VYuR8vzUSUohnv9vUzX+c+I7V31rq0k1sGJDYhrINJ0kgMqOUb8/A6IfK9Z/n0VH/LHS3eLHCuxhl12p28zM4uJHChFrYtX+0q0tcDUOnvpAfEZkxavRewbVcr/mrKmp0f/LBWLlA8BAACYEARHAAAAIQRHAAAAIeQc7YXKL+p/9wmuLV8jck18SoANzfZtxQaxJk+TH0Otqfa5OlUpn2dxZOMO19aY8usczU73ubYZIpdoVtKvN6TyO9IiT6UuLl6H2FfVM0uJXImsWPcnJdpUXo7Kt0mJ9YvUmkGlIFquk1q/KOpaRSo3aVict6/sr2nUtYUGRa6JskvUylL1vfIiv6i77PNy1L5qTZ60eB+odYS6RJvKG3p52K+ho+r9qbWAhkRRJ5UPlBN5Q3lxnX/dNce1ZUXdrnxJ1FET9bjE29m6xbo/Kn+nVPSfV4lktDyfTL3/HBoeEAWwxDo9MVGfS+X05Gt9X0T5MQu2+fPG1VJFVf5i1Wz1G6qc0KLYNynyizJdft+aV/w1GBS/A7Ji37wvf2etP1cLBpn1L/DJROp3ksrfyYikw+EGkRcl8otqX/C5rMld/iT9i2a6tsSQ+Mzv9T8PXcf4n/WJxJMjAACAEIIjAACAEIbV9iJ5xOGuTT1yHJjr9y35J9JWnO2nnlY3+keTcxr8UFa1eAR/SHW3a2vNiH3FtPiWVI9rq4v7vrQm/XZZMYSmhsbUdO2MGC5riquSDv4c1WK6u5oCr6a2K2qqfEEMZWXF8JuSE0Njg2L8I+qQVzymyl/415sTw0IF88frK/vH7zkxhKZLXfghjG5RJkOVZVDDap0FP+VaDW/pIS/f52LEIcPuvL8G6YS/b7uG/GtT0+LV1PZc3ve5qdYPbXf0+GuQFMNgA/3+wyTYJdrq/H1TVW6Snf5eFrJiSDgvdt7q943NEGV4dovyNeI3Tc0rvq2cFO97PzJjjZv9/cjXiZ+tVt9WtdPv2/wb/3mghnrUMJaifle0PuE/i0tV/r1brPJ9jnpeM7Om//tz16Z+nw0t9EPPDZ1iavzjv4p2YnGOrDieGn5TfaneKda62bMvET/v9wdPjgAAAEIIjgAAAEIIjgAAAELIOdoHIh3DcoeIsepaP7bcXOdzeuozfjz2iHo/Hb8+6bebnfb5RU0JP90zG/f9a1TbiWn2jaLcR52aVys0iCn1ipruXi3elv0iD0lJiVydvrJ/beq8iirF0VcWU6RFgkefyOkZECUnyuJvFDVVvizOsb3oc1dUfpG6v+p4Kg+pq+jXpvh9LtpU+X4xLb4q4fvSmfOvY7Aorl/BH69UFstpiCnw3f2i/MWwf69lq/3Pb1GVyegX09jFtPjBbWJtD1FOQ6rx1zQlpsCXh/21ihfFUhL+R9qqt/jjDfvbYflG31b7grj2Yiq6yi9KiPSWxq3+2vcuFPdclKWY1e5zbdTyK5muaJ8lKtcmU7/UtVVt8aWSVA5N6odP+bYTFrs2lZNT3Pw72UWVSxTx1cl+P/jc9a7tTSd+zu8sro3cVyyJo/r3k+d8uS5FHW+88OQIAAAgZNoGR+3t7bZo0SJbtmxZpbsCAAAmkWkbHLW1tdmmTZts/fr1le4KAACYRGJBoBajnz56e3utoaHBenp6rL7+jwPtp551g9t265v9uH7mL3zuj1qraG61b2tM+Tyk11R3uraEyHtR6xctSIvyIWL9opKIiVVpj6a4f2uoXJ2UOJ5ag0iWyTDfpsp4lMW6SYWIb92c2E79VdAjFmNRuURqDaftJZ9XIktiiBIbvSVfmkLlDan3QS7w78ndRb/s/q68718h8Fdh25Cvo9CS9eVmduX9OVSOUF/er8mTivvrkhP7dg9EW9tF5gPl/L1UaxUFQyLtMu+vS1Dl+5zs9vuq8hIijcvKaVHSpk/8bPlLL/OGRLpX5PIQySGxHpfI6VHrzvQd6q/9nB/6zzCl+/hZru2xu//etamcF5mX0+zf41FLWtRu2rnXfoapvBolap/V8c4+6grXtk85R2Lbh8rRcnomE3Ud9rxee/v9PRam7ZMjAAAAheAIAAAghOAIAAAghHWO9kLV1Sm3+DH82qzP/TmqYbtrm532yQPz07tcW41YW0jXOPP9axELhwyLelfVIpdISYm6Yqr+WG1crUvkr0tG1PxS1DlU3pB6bap+W07k1gyIXJ0uUS8spfKGxHYJcd6Oos8l+n2hye8rcpj6RB7SkFgjqbfgt1M1yfqLPvenV+QDxcXreLZntmtTugd9jpDK8+nr9duV8+K9Ie6v5fy9jBf8duqvPlWCLdXnt6wSKTPFGr9dUaREpf2Pqom3kBVrxGsTS4ip/KKqXX5DlTekPHb3Za5N5Veqml+qTpaZf/9FpfKLVJ7JQwew9k5RrCOUvf8Jv6/I3YmaXySpdYAi5v0c0HmnkEpfB54cAQAAhBAcAQAAhDCsthevvNFPVa6q8c/M3zBrq2tbVO3Xyo869b5ODKtlRd2SOvGsviYupui7FrOMuO1qqryaUp+Nq6E23xc1vV+dY0CU4kiIEYcuMR+6IEpsqKn3hcC/3u0lX+OgRwyXKXlxvJeGm12bGpLrFHUZhkr+tRXFUGBObLcr5/scV8N0YghNDQXmCv619e7276JADGXFEmJoZ0iMZYnSGWp4S1wCi6mhJ79ahSxXURBVPNQU+HSffx3FnO9z/Rb/M61KMqgp5v0L/Jhc7Qv+hajtGjf4zw1FTZVX5RdUCYuf/ODySOeIOu1cTSVXfRnz7X7+Gdd2ICbT690XUctuHIxT/scLT44AAABCCI4AAABCCI4AAABCyDnaCzHT2+qrfD7QUaLcx6Lsy65N5RKpkh2KKuNhIrdG5giJ6filiGU3VI7QYNknaaiSIsoOMaU5JfbdUfJT1nNienpf4KcRD5R9bo3KJVJ5Q/1i+vz2gs8R6i/6/qnp870Fny+SFEkznYM+J0WV08iX/L3M5f15h/r9NRBvDRO314Ky+Hupz58jLnN/RJ5Zxp8k3eU7k/bVdWzYr3pgmS7fVtMZbWr7oFhKom6rzxsabvKvt/kxvzyHyiVSbSrvJWruSu3jrklOT1fnGOvclfHIhdmTui6qRIac8h+xLwfyOsZ6evlEXNPxOuZUx5MjAACAEIIjAACAEIIjAACAEHKO/p+3NrzPkrFQnscVJ7ptXt/sF085TKxV1Bj365U0xX1+UaMou6HyhhSV56NyiVSbyiVS+gKfX1RQu4pz9IkcnLzI+RgU+UV5sX5RR7HBtaXFOkK7Sj7nY1iskfQ//Ye4tsaUv28dOZ9zpPTmfb7Sy32+z4Wif20qBWywzx8vkfK5NaVuf/3iw/5vnsSgWP+p2d/f7MvRPhLqXxRrAYlyGpkev93uo6K9x1V+0SHf9zl+/YtmujZV6iIryoIkd/W7tp/8wOeVvPHtX3BtP7334/6AEUVeG0eUyYi6ds9Y5+CMdd7KWOc/TcR5D0TUNY0mykTlOx2seHIEAAAQQnAEAAAQEguCiPO6p6je3l5raGiwnp4eq6//4xDKgn+5yW371qW/cG2rmx91ba0JP1yhpqyrKvW1cVFlPYg25V+V8Rg231YQt3xAvAsKqiq6ELXCvZpmnwv8sFBnwQ9HqbIguwu+rMWgqFzfLabUq2n2XTnfVhCl3Lv6xdIAYkq9mhafSPn7kd8lKpuLS5/oEUNyojpHutfvLEY5LTng27K++oWcFp/p8UN8umr7GBPVztVUbzWlPurxlMlWDqIS5wAmm739/h4LPDkCAAAIITgCAAAIITgCAAAIYSr/XsTEdOhDMt2uLSumk6uIsz4mSjoIu8uDri0ljqhyiZS+spj+LRJa1NT7vrLPhcmJshuqFEe3yDlSZTxU2/Z8nW8b9m2DooxHWeRJdfT7fVUpjmTCX6vdu31eUynnX29MlHiJ9as8JNdk1R3iHSOaqsRU9IJIrRGrI9is5/37ZaDVX4P6Lb6cRqrX57ypKfCqtMLK//NF11b7TV8To/Dmpa5tsNXnsjXv8vlFE1HSQYma0zMRJTsmIr+IvCZMJzw5AgAACCE4AgAACCE4AgAACCHnaC+Cep9nkRElQOriPomkQaxVlIr5Sz0UDLu2bMzngfSV/Xn7xLpEZdGmcoR65XpDPr9D5xz5hJauos/L2Vn0eT5qzaAdIr9ox7BPpOkviPMO+XyluCiNkiuI17bb72s5379YwecwidQkS/X4vzNiMr/It/Uf6tvqXvBtCVG7Jcj5/kUtsdH8K7HQUcR1f0ysLXT2UVe4ttqIZTLUvqkf+n39CmLRjXUu0USgL0BlRAqObrnlln0+8IUXXmh1df4XHwAAwGQWKTi69NJL7dBDD7VEQizHK2zdutX+6q/+iuAIAAAcdCIPqz311FM2e/bsSNsSFAEAgINVpODommuusdraiLWKzOxTn/qUNTU17XenKuGtq260ZPKPOTbxd/s6W0opYmW6/rKvOzUY+AyKgsqZEWv3RF2XaEDkCL1SmOHa1FpFu0UuUSnwuTW9JX+ttg76c6Tj/vX+b2+za0vF/Zo8u4f8OXJ5n0tUGPavo9wtFv4RUxESA6JR5A1ld4g6eb1+u2FR6qemw7+2GpGHNNzg+9L82Ha/odB9/CzXptYWUlSdMrWO0JtO/Jxre+jnn4l0jsm0Xs5kWh9oMl0XYLqLHBztiyuvvHK/OgMAAFBpkafyb9/+p/9qLRaL9uSTTx5whwAAACopcs7RnDlzbNu2bSN5R8cee6x9//vft3nz5pmZ2a5du2zlypVWKkUrazHZ3P+jT1p9/R/HQRbccaPbZrjsh3EGxTBTIfDXQJX7GAzUEJo/3gtFP0SVEsfbUfLjOGqafVfRD5EOi2E6Nc1+d15MgRe6836Ib7jkz9HRJcaexDBiqeCvS1zMqY//3g+/qan3YlUGS/f4toRfbcHEpZJlN0pVfgKDKtnR8rCfep8VU+AfjDjEoqbFqynwUYdsDmS4p1JDSpNpiGoiyowAGFuRnxwFe/wif+GFF6xQKPzJbQAAAA42Y7pCdizm/+IHAAA4mFA+BAAAICRyzlEsFrO+vj7LZrMWBIHFYjHr7++33t4/zGF+9b9TRSDyVErmn4zlAp9D8mLJ71wQ25n5tu6yz5npKDS6tsGyn55eErHu1mE/VX6XyBvKi0SaF/t8rlM85odO80VRomTIlygZGvBtQUnE54P+usSHxXaiLypHSLVV+TQf2/Avl7m2N779C66tdtNO1za00F+rYpV/v/ziq/4cKkdITamPWv5CmYj8ogM571jn25C/A+BARA6OgiCwI488ctTXxx133KivGVYDAAAHu8jB0Y9//OPx7AcAAMCkEDk4OuWUU8azHwAAAJNCLIg4/75cLttNN91k3/3udy2fz9sZZ5xh11xzjVVVRSuzMVn19vZaQ0ODnWpvtWTsj+sY/faWE9y2q09d59pOq33GtTXGh1zbrrLP81F5SM8Oz3VtL+Rmurahkl9zqajyn/p9LkxZrCM0kPf5QIOiPEep7HN/BnrEe2DI9yXRJ9YqKvi+pPr84Rpe8HlcQ83+eFW7/HaltFg3KRVtCLjp//7ctal8FpU3FFVRrGkUdS0ghdIUAKaDV39/9/T0jFqncCxEnq127bXX2qc+9Smrra21Qw45xL785S9bW1vbmHYGAACg0iIHR9/4xjfsn//5n+2///u/7b777rMHHnjA7rrrLiuXxbQuAACAg1Tk4Oill16yc845Z+TrVatWWSwWs1deeWVcOgYAAFAJkROyi8WiZbOj62WlUilXQmSqCBI+Fau/5OuFvSzqng0k/BpEfSWfl/Ncbo5r6xXb9RX9eVWNs5d6Gl1bVdrfn519vt5aPu/fCsUd0fLJ4kWfvzPjN75tqMXvmxyIdArL1/g4vu73vr5ctjPn2gr1Pncq1euvy0M//4xre5PIOZK1y0Te0FibiPWGAAD7uM7RBRdcYJnMHxN3c7mcffCDH7Samj/+sv3Od74ztj0EAACYQJGDo9WrV7u297znPWPaGQAAgEqLHBzdcccd49mP/dbd3W2rVq2yYrFoxWLRLrnkErv44ov3+Tj393xj1FTAo/7xn9w2Lw81uLZD0r4tV/bDOFuGZ7m2zmE/9XDHcK0/b58/R3/OD90VS376/I6ORtdmohRHsldMvY/7ocXsDjH1ftCfItfk22p/74/X8LzfOdfihxFVyY5is79WyV39/sSPR5sqH9WDz10faTumzwPAwSlycDRZ1dXV2bp166y6utoGBgbsda97nZ133nnW3OxrigEAAPw5kWarnXfeeftUWPb888+37du373en9kUikbDq6j8kJw8PD1sQBBZxXUsAAAAnUnB0//33244dO6y3t/fP/uvp6bEHHnjA+vvF8Iawbt06O/fcc23u3LkWi8Xsvvvuc9u0t7fbggULLJvN2ooVK+zJJ58c9f3u7m5bvHixHXroofaJT3zCZs70K0oDAABEEWlYLQgCO/LII8elAwMDA7Z48WK76KKL7LzzznPfv+eee2zNmjV222232YoVK+zmm2+2M88805577jmbPXu2mZk1Njbar371K+vs7LTzzjvP/vqv/9paWsS88X2Q8BVALF/2l6tHTL3vLPhcop15nx/zfI/PQ1LT7JXBPp+XYwO+f6ndPpeoVOMX7kyK1xvEfX5RIu+3S/jZ87ot75/oqfwiNR0/cp7PiZ9zbYWFSyPtO9bILwKAg1Ok4OjHP/7xPh/4kEMOibTd2WefbWefffZev/+lL33JLr74YrvwwgvNzOy2226z//qv/7Kvf/3rdsUVo9ebaWlpscWLF9tPf/pT++u//mt5vOHhYRseHh75el+GCwEAwNQXKTg65ZRTxrsfUj6ftw0bNtiVV1450haPx23VqlX22GOPmZlZZ2enVVdXW11dnfX09Ni6devsQx/60F6Ped1119lnP/vZce87AAA4OEUuH1IJO3futFKp5IbIWlparKOjw8zMXnzxRXvjG99oixcvtje+8Y320Y9+1I499ti9HvPKK6+0np6ekX9bt24d19cAAAAOLgf9VP7ly5fbxo0bI2+fyWRGrfK9N9ndvu2Vfp9L1JL1w3IqN+nZ7tmurWfQ5ysN9fu+BX1+3aR0l88lSokRwsDvajOeU7lEPg9Jleyo3ulLdiSGfFvVFnEBBZVLpNYHikqVAFGirkE02fOGWEsJAMbepH5yNHPmTEskEtbZ2TmqvbOz01pbWyvUKwAAMJVN6uAonU7b8ccfb2vXrh1pK5fLtnbtWlu5cmUFewYAAKaqig+r9ff3229/+9uRr7ds2WIbN260pqYmmz9/vq1Zs8ZWr15tS5cuteXLl9vNN99sAwMDI7PXAAAAxtI+BUff//737Tvf+Y41NTXZRRddZEcfffTI93bv3m3veMc77OGHH96nDjz11FN22mmnjXy9Zs0aM/tDods777zT3vWud9mOHTvs6quvto6ODluyZIn94Ac/OOB1jP4ckTZku8QaRDtq61xb55Bf02gw72uh9Yu1ihKv+LZ4wfclLtYbynb5tkQh2mrhtS/4hY56jqz226kaZ5t97bKiOMdY58IcSL7NgfRlMuX5kF8EAGMv8rDa3XffbW95y1uso6PDHnvsMTvuuOPsrrvuGvl+Pp+3Rx55ZJ87cOqpp46U/Aj/u/POO0e2+chHPmIvvviiDQ8P2xNPPGErVqzY5/MAAABEEfnJ0U033WRf+tKX7GMf+5iZmX3rW9+yiy66yHK5nL3//e8ftw6Ol/b2dmtvb7dSyc+0AgAA01fk4Gjz5s127rnnjnz9N3/zNzZr1ix7y1veYoVCwd7+9rePSwfHS1tbm7W1tVlvb681NDS478/aMODafnuCv1xqCE2VACkU/L7xTj9tPy2m49f+3g+NNTw/6Nr6F/ilARo37HBtavr88gu+5NqevHONazuQIaWzj7rCtUUtCwIAwESJHBzV19dbZ2enLVy4cKTttNNOs+9973v2V3/1V/b73/9+XDoIAAAwkSLnHC1fvtwefPBB137KKafYAw88YDfffPNY9gsAAKAiIgdHl112mWWzohK8/SGp+oEHHrD3ve99Y9YxAACASog8rHbKKaf8yQK0p5122qgp+VNRMOAvV8+QDxiDwJfnyPf5qfzpvN8u6VOdrPmx7a6t2OxznWq/+bjfzh/O3nTi5/w5dvWLLceWmvKvTMR0/APB9HkAmNoiB0fFYtFKpdKoumSdnZ1222232cDAgL3lLW+xk08+eVw6CQAAMFEiB0cXX3yxpdNpu/32283MrK+vz5YtW2a5XM7mzJlj//RP/2T333+/nXPOOePWWQAAgPEWOefo0UcftXe84x0jX3/jG9+wUqlkmzdvtl/96le2Zs0au+mmm8alkwAAABMl8pOjl19+2Y444oiRr9euXWvveMc7RtYIWr16td1xxx1j38MKKdSnXFt8KOHa+vr82kLlgt8u0ecvdb1IwZm9rtO1qfyipMgRUvlFySMOd20P/vwzrk2tX6QcSL4NuToAgINB5CdH2WzWhob+WH/r8ccfH1XGI5vNWn//+Cf1jpX29nZbtGiRLVu2rNJdAQAAk0jk4GjJkiX2r//6r2Zm9tOf/tQ6Ozvt9NNPH/n+7373O5s7d+7Y93CctLW12aZNm2z9+vWV7goAAJhEIg+rXX311Xb22Wfbt771Ldu2bZtdcMEFNmfOnJHv33vvvXbSSSeNSycrIdVbcG3Jfj8dP1/th9+s7Kfox/zhLFHwZUF2rZzt2hq+8Zhre1AMUR1IeY6oQ16TqSI9AADjYZ/WOdqwYYP98Ic/tNbWVnvnO0f/klyyZIktX758zDsIAAAwkSIHR2Zmr33ta+21r32t/N4HPvCBMekQAABAJUXOOXrVddddZ1//+tdd+9e//nW74YYbxqRTAAAAlbJPT47MzG6//Xa7++67Xfsxxxxj7373u+3yyy8fk45NRrGyb4v3+2n76W6fc6T5nCOVXxQ1p0flF411jhD5RQCAqW6fnxx1dHSMSsR+1axZs2zbtm1j0ikAAIBK2efgaN68efboo4+69kcfffSgmsoPAACg7POw2sUXX2yXXnqpFQqFkXWO1q5da5/85Cft7//+78e8gwAAABNpn4OjT3ziE7Zr1y778Ic/bPl83sz+sDr25ZdfbldeeeWYd7BSVHmOWLHGtWV3+PyiQp0/Xt1Lvq3h+cH96ptZ9FyiqbJ+0WTvHwBg6tjn4CgWi9kNN9xgV111lT3zzDNWVVVlRxxxhGUymfHo37hpb2+39vZ2K5VKle4KAACYRPY5OHpVbW3tQV2XrK2tzdra2qy3t3ekeC4AAMA+J2QDAABMZfv95GiqK27+nWurf9HXPRuaKdY06vNN5Ul+pSd7/s5k7x8AYOrgyREAAEAIwREAAEDIJB/sqRw1jHPqWb523NDMtGuLiQlwaljtoZ9/xrWpKetR+wcAAA4cT44AAABCCI4AAABCCI4AAABCyDnaB4OtKdeWGvDb1W/Ju7aqLbtd25u+vP8lMSinAQDA+ODJEQAAQMi0DY7a29tt0aJFB3UJFAAAMPambXDU1tZmmzZtsvXr11e6KwAAYBIh52gvVE5PwwmLXVuh3uchlaoSrk2VI1E5QmcfdYVre/C56yPtCwAADty0fXIEAACgEBwBAACEEBwBAACEkHO0F1Frq6V6C65NrWlUFOeIml8EAAAmDk+OAAAAQgiOAAAAQhhW2wdquCzqMBhDaAAAHBx4cgQAABBCcAQAABBCcAQAABASC4IgqHQnKqm3t9caGhqsp6fH6uvrK90dAAAQwXj+/ubJEQAAQMi0DY7a29tt0aJFtmzZskp3BQAATCIMqzGsBgDAQYdhNQAAgAlCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABBCcAQAABAybYOj9vZ2W7RokS1btqzSXQEAAJNILAiCoNKdqKTe3l5raGiwnp4eq6+vr3R3AABABOP5+3vaPjkCAABQCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCCI4AAABCpm1w1N7ebosWLbJly5ZVuisAAGASiQVBEFS6E5XU29trDQ0N1tPTY/X19ZXuDgAAiGA8f39P2ydHAAAACsERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABACMERAABAyLQNjtrb223RokW2bNmySncFAABMIrEgCIJKd6KSent7raGhwXp6eqy+vr7S3QEAABGM5+/vafvkCAAAQCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACDnog6OtW7faqaeeaosWLbLXv/719u1vf7vSXQIAAAexZKU7cKCSyaTdfPPNtmTJEuvo6LDjjz/ezjnnHKupqal01wAAwEHooA+O5syZY3PmzDEzs9bWVps5c6Z1dXURHAEAgP1S8WG1devW2bnnnmtz5861WCxm9913n9umvb3dFixYYNls1lasWGFPPvmkPNaGDRusVCrZvHnzxrnXAABgqqp4cDQwMGCLFy+29vZ2+f177rnH1qxZY9dcc4394he/sMWLF9uZZ55p27dvH7VdV1eXve9977Ovfe1rE9FtAAAwRcWCIAgq3YlXxWIxu/fee+1tb3vbSNuKFSts2bJlduutt5qZWblctnnz5tlHP/pRu+KKK8zMbHh42N70pjfZxRdfbO9973v/5DmGh4dteHh45Ouenh6bP3++bd261err68f+RQEAgDHX29tr8+bNs+7ubmtoaBjTY0/qnKN8Pm8bNmywK6+8cqQtHo/bqlWr7LHHHjMzsyAI7IILLrDTTz/9zwZGZmbXXXedffazn3XtDMUBAHDw2bVr1/QKjnbu3GmlUslaWlpGtbe0tNizzz5rZmaPPvqo3XPPPfb6179+JF/pX//1X+3YY4+Vx7zyyittzZo1I193d3fbYYcdZi+99NKYX1zsm1f/CuApXuVxLyYP7sXkwb2YXF4d+WlqahrzY0/q4CiKk08+2crlcuTtM5mMZTIZ197Q0MCbfZKor6/nXkwS3IvJg3sxeXAvJpd4fOzTpyuekP2nzJw50xKJhHV2do5q7+zstNbW1gr1CgAATGWTOjhKp9N2/PHH29q1a0fayuWyrV271lauXFnBngEAgKmq4sNq/f399tvf/nbk6y1bttjGjRutqanJ5s+fb2vWrLHVq1fb0qVLbfny5XbzzTfbwMCAXXjhhWNy/kwmY9dcc40casPE4l5MHtyLyYN7MXlwLyaX8bwfFZ/K/5Of/MROO+0017569Wq78847zczs1ltvtZtuusk6OjpsyZIldsstt9iKFSsmuKcAAGA6qHhwBAAAMJlM6pwjAACAiUZwBAAAEEJwBAAAEDKtg6P29nZbsGCBZbNZW7FihT355JOV7tKUd91119myZcusrq7OZs+ebW9729vsueeeG7VNLpeztrY2a25uttraWnvHO97h1rrC2Lv++ustFovZpZdeOtLGvZg4L7/8sr3nPe+x5uZmq6qqsmOPPdaeeuqpke8HQWBXX321zZkzx6qqqmzVqlW2efPmCvZ46iqVSnbVVVfZwoULraqqyg4//HD7x3/8Rwun6HI/xse6devs3HPPtblz51osFhupfPGqKNe9q6vLzj//fKuvr7fGxkZ7//vfb/39/fvUj2kbHN1zzz22Zs0au+aaa+wXv/iFLV682M4880zbvn17pbs2pT3yyCPW1tZmjz/+uD300ENWKBTszW9+sw0MDIxsc9lll9kDDzxg3/72t+2RRx6xV155xc4777wK9nrqW79+vd1+++32+te/flQ792Ji7N6920466SRLpVL24IMP2qZNm+yLX/yizZgxY2SbG2+80W655Ra77bbb7IknnrCamho788wzLZfLVbDnU9MNN9xgX/3qV+3WW2+1Z555xm644Qa78cYb7Stf+crINtyP8TEwMGCLFy+29vZ2+f0o1/3888+3p59+2h566CH73ve+Z+vWrbMPfOAD+9aRYJpavnx50NbWNvJ1qVQK5s6dG1x33XUV7NX0s3379sDMgkceeSQIgiDo7u4OUqlU8O1vf3tkm2eeeSYws+Cxxx6rVDentL6+vuCII44IHnrooeCUU04JLrnkkiAIuBcT6fLLLw9OPvnkvX6/XC4Hra2twU033TTS1t3dHWQymeDf//3fJ6KL08pf/uVfBhdddNGotvPOOy84//zzgyDgfkwUMwvuvffeka+jXPdNmzYFZhasX79+ZJsHH3wwiMViwcsvvxz53NPyyVE+n7cNGzbYqlWrRtri8bitWrXKHnvssQr2bPrp6ekxMxspHLhhwwYrFAqj7s3RRx9t8+fP596Mk7a2NvvLv/zLUdfcjHsxkb773e/a0qVL7Z3vfKfNnj3bjjvuOPuXf/mXke9v2bLFOjo6Rt2LhoYGW7FiBfdiHJx44om2du1ae/75583M7Fe/+pX97Gc/s7PPPtvMuB+VEuW6P/bYY9bY2GhLly4d2WbVqlUWj8ftiSeeiHyuiq+QXQk7d+60UqlkLS0to9pbWlrs2WefrVCvpp9yuWyXXnqpnXTSSfa6173OzMw6OjosnU5bY2PjqG1bWlqso6OjAr2c2r75zW/aL37xC1u/fr37Hvdi4vzv//6vffWrX7U1a9bYpz71KVu/fr197GMfs3Q6batXrx653uozi3sx9q644grr7e21o48+2hKJhJVKJbv22mvt/PPPNzPjflRIlOve0dFhs2fPHvX9ZDJpTU1N+3RvpmVwhMmhra3NfvOb39jPfvazSndlWtq6datdcskl9tBDD1k2m610d6a1crlsS5cutc9//vNmZnbcccfZb37zG7vtttts9erVFe7d9POtb33L7rrrLrv77rvtmGOOsY0bN9qll15qc+fO5X5ME9NyWG3mzJmWSCTcrJvOzk5rbW2tUK+ml4985CP2ve99z3784x/boYceOtLe2tpq+Xzeuru7R23PvRl7GzZssO3bt9sb3vAGSyaTlkwm7ZFHHrFbbrnFksmktbS0cC8myJw5c2zRokWj2l772tfaSy+9ZGY2cr35zJoYn/jEJ+yKK66wd7/73Xbsscfae9/7XrvsssvsuuuuMzPuR6VEue6tra1uYlWxWLSurq59ujfTMjhKp9N2/PHH29q1a0fayuWyrV271lauXFnBnk19QRDYRz7yEbv33nvt4YcftoULF476/vHHH2+pVGrUvXnuuefspZde4t6MsTPOOMN+/etf28aNG0f+LV261M4///yR/+deTIyTTjrJLWnx/PPP22GHHWZmZgsXLrTW1tZR96K3t9eeeOIJ7sU4GBwctHh89K/HRCJh5XLZzLgflRLluq9cudK6u7ttw4YNI9s8/PDDVi6X960m6wGnkx+kvvnNbwaZTCa48847g02bNgUf+MAHgsbGxqCjo6PSXZvSPvShDwUNDQ3BT37yk2Dbtm0j/wYHB0e2+eAHPxjMnz8/ePjhh4OnnnoqWLlyZbBy5coK9nr6CM9WCwLuxUR58skng2QyGVx77bXB5s2bg7vuuiuorq4O/u3f/m1km+uvvz5obGwM7r///uB//ud/gre+9a3BwoULg6GhoQr2fGpavXp1cMghhwTf+973gi1btgTf+c53gpkzZwaf/OQnR7bhfoyPvr6+4Je//GXwy1/+MjCz4Etf+lLwy1/+MnjxxReDIIh23c8666zguOOOC5544ongZz/7WXDEEUcEf/u3f7tP/Zi2wVEQBMFXvvKVYP78+UE6nQ6WL18ePP7445Xu0pRnZvLfHXfcMbLN0NBQ8OEPfziYMWNGUF1dHbz97W8Ptm3bVrlOTyN7Bkfci4nzwAMPBK973euCTCYTHH300cHXvva1Ud8vl8vBVVddFbS0tASZTCY444wzgueee65CvZ3aent7g0suuSSYP39+kM1mg7/4i78IPv3pTwfDw8Mj23A/xsePf/xj+Tti9erVQRBEu+67du0K/vZv/zaora0N6uvrgwsvvDDo6+vbp37EgiC05CcAAMA0Ny1zjgAAAPaG4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AgAACCE4AjAlPfCCy9YLBazWCxmS5YsGdNj33nnnSPHvvTSS8f02AAqg+AIwKTyyCOP2Omnn25NTU1WXV1tRxxxhK1evdry+byZmeVyObvgggvs2GOPtWQyaW9729siH/tHP/rRqKKV//AP/zAS2CSTSVuwYIFddtll1t/fb2ajg6o9/z3++ONmZvaud73Ltm3bRsFRYApJVroDAPCqTZs22VlnnWUf/ehH7ZZbbrGqqirbvHmz/ed//qeVSiUzMyuVSlZVVWUf+9jH7D//8z/36fjNzc3W3Nw8qu2YY46xH/3oR1YsFu3RRx+1iy66yAYHB+32228f2eZHP/qRHXPMMe5YZmZVVVVWVVVl6XR6f14ygEmI4AjAhCqXy/aFL3zBvva1r9nWrVutpaXF/u7v/s4+/elP2w9/+ENrbW21G2+8cWT7ww8/3M4666yRr2tqauyrX/2qmZk9+uij1t3dfUD9SSaT1traamZ/eAq0du1a++53vzsqOGpubh7ZBsDUx7AagAl15ZVX2vXXX29XXXWVbdq0ye6++25raWkxM7PW1lbbtm2brVu3rmL9q6qqGhnCAzA98eQIwITp6+uzL3/5y3brrbfa6tWrzewPT4ZOPvlkMzN75zvfaf/93/9tp5xyirW2ttoJJ5xgZ5xxhr3vfe+z+vr6ce/fhg0b7O6777bTTz99VPuJJ55o8fjovyVfzUsCMPXw5AjAhHnmmWdseHjYzjjjDPn9RCJhd9xxh/3+97+3G2+80Q455BD7/Oc/b8ccc4xt27ZtXPr061//2mpra62qqsqWL19uK1eutFtvvXXUNvfcc49t3Lhx1D8AUxfBEYAJU1VVFWm7Qw45xN773vfarbfeak8//bTlcjm77bbbxqVPRx11lG3cuNGeeeYZGxoasu9+97sjw3yvmjdvnr3mNa8Z9Q/A1EVwBGDCHHHEEVZVVTVqOv2fM2PGDJszZ44NDAyMS5/S6bS95jWvsQULFjDjDICZkXMEYAJls1m7/PLL7ZOf/KSl02k76aSTbMeOHfb000/b+9//frv99ttt48aN9va3v90OP/xwy+Vy9o1vfMOefvpp+8pXvjJynE2bNlk+n7euri7r6+sbGeYa6wUeX7Vr1y7r6OgY1dbY2GjZbHZczgegsgiOAEyoq666ypLJpF199dX2yiuv2Jw5c+yDH/ygmZktX77cfvazn9kHP/hBe+WVV6y2ttaOOeYYu+++++yUU04ZOcY555xjL7744sjXxx13nJmZBUEwLn1etWqVa/v3f/93e/e73z0u5wNQWbFgvD5NAGCSeOGFF2zhwoX2y1/+ctyeLp166qm2ZMkSu/nmm8fl+AAmDjlHAKaNE0880U488cQxPeZdd91ltbW19tOf/nRMjwugcnhyBGDKKxaL9sILL5iZWSaTsXnz5o3Zsfv6+qyzs9PM/pCHNHPmzDE7NoDKIDgCAAAIYVgNAAAghOAIAAAghOAIAAAghOAIAAAghOAIAAAghOAIAAAghOAIAAAghOAIAAAg5P8H+aPJ+VOYpIcAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from matplotlib.colors import LogNorm\n", "\n", "key = apt.randgen.get_key()\n", "key, (cs1, cs2, eff) = (\n", " tree.likelihoods[\"rn220_llh\"].components[\"rn220_er\"].simulate(key, int(1e6), parameters)\n", ")\n", "\n", "plt.hist2d(\n", " cs1, cs2, weights=eff, bins=[np.linspace(0, 100, 100), np.logspace(2, 4, 100)], norm=LogNorm()\n", ")\n", "plt.yscale(\"log\")\n", "plt.xlabel(\"cS1 [PE]\")\n", "plt.ylabel(\"cS2 [PE]\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.16" } }, "nbformat": 4, "nbformat_minor": 5 }