{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "5efbd85a-a55f-41a9-a34e-452151c0a9a7", "metadata": {}, "outputs": [], "source": [ "import pandas as pd\n", "import servicex as sx\n", "import awkward as ak\n", "from func_adl_servicex import ServiceXSourceUpROOT\n", "import matplotlib.pyplot as plt\n", "\n", "plt.rcParams['figure.dpi'] = 200 # make figures bigger and more readable" ] }, { "cell_type": "code", "execution_count": 2, "id": "1caae859-3d03-44c1-ba5c-88da8b576dd6", "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "", "version_major": 2, "version_minor": 0 }, "text/plain": [ "cms:DYJetsToLL_M-50_...: 0%| | 0/9000000000.0 [00:00]" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "", "version_major": 2, "version_minor": 0 }, "text/plain": [ " cms:DYJetsToLL_M-50_... Downloaded: 0%| | 0/9000000000.0 [00:00]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dataset_name = \"cms:DYJetsToLL_M-50_TuneCP5_13TeV-amcatnloFXFX-pythia8/RunIIAutumn18NanoAODv7-Nano02Apr2020_102X_upgrade2018_realistic_v21_ext2-v1/NANOAODSIM\"\n", "sx_dataset = sx.ServiceXDataset(dataset_name, \"uproot\")\n", "ds = ServiceXSourceUpROOT(sx_dataset, \"Events\")\n", "filtered_missing_ET_5 = ds.Where(lambda event: Zip({'p4': Zip({'pt': event.Muon_pt,\n", " 'eta': event.Muon_eta,\n", " 'phi': event.Muon_phi,\n", " 'mass': event.Muon_mass}).ToFourMomenta(),\n", " 'charge': event.Muon_charge})\n", " .Choose(2)\n", " .Where(lambda pair: pair[0].charge * pair[1].charge < 0)\n", " .Select(lambda pair: (pair[0].p4 + pair[1].p4).mass)\n", " .Where(lambda mass: 60 < mass and mass < 120)\n", " .Count() > 0\n", " ).Select(lambda event: event.MET_pt).AsAwkwardArray().value()" ] }, { "cell_type": "code", "execution_count": 3, "id": "60ef3cd2-a47a-4dfc-8633-e6ec1fb045e0", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "filtered_missing_ET_5" ] }, { "cell_type": "code", "execution_count": 4, "id": "4cd0af13-76b9-4859-bf65-3a8eb19aa18e", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABG4AAALeCAYAAAD711clAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAB7CAAAewgFu0HU+AABv40lEQVR4nOzdebxdZXno8d8jkSEoQwEBSTAKVCiDA4OgjE5YgwooilgVBRFtnSpKHO4Vh9YI2qJAKyAV5Vq5YkGEYJFbAYWLQLzWapkESZswCSpDIIAhz/1jvZuzstnTOexzzkr27/v57M9+373e9Zxnh6VwnrxDZCaSJEmSJElqnqdMdwKSJEmSJEnqzMKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUDOmOwFpMkTEWsCOpXs38Ng0piNJkiRJWv2tAWxS2r/MzEeGEdTCjVZXOwLXTncSkiRJkqSRtCuwcBiBXColSZIkSZLUUM640erq7lbjmmuuYfPNN5/OXCRJkiRJq7k77riD3XbbrdW9u9fY8bBwo9XV43vabL755syaNWs6c5EkSZIkjZah7bPqUilJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUDOmOwFJ0uphzrwFTzrGovlzh5CJJEmStPqwcCNJGsgwCjOSJEmSxselUpIkSZIkSQ3ljBtJUmP0m9XjUipJkiSNGgs3kiSXQUmSJEkN5VIpSZIkSZKkhrJwI0mSJEmS1FAulZIkrTIGWdLlPjiSJElanTjjRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGco8bSRoBHvctSZIkrZqccSNJkiRJktRQzriRJK1W+s0u8tQpSZIkrUqccSNJkiRJktRQI1m4iYjLIiLH+dq3R7xXRcS5EbEkIh4p7+dGxKvGkdPMiPhIRFwTEb+PiKURcX1EfDEithxHnO0j4qsRcXNELIuIuyPixxHx7ogYeIZVRBwaERdHxB0R8XBELIqIsyJi93HE2CgiPh0Rv4iI+yLi/tL+dERsNGgcSZIkSZJGlUulBrMC+HX7hxERwFeBo9oubQEcBBwUEacBR2dmdgseEVsBC4Dntl3atryOjIjDMvOiXklGxBHAKcBatY/XBvYqr8Mj4oDM/F2PGGsD5wAHtF16VnkdFhHHZeZn++SyK3A+sHnbpZ3K68iIeF1mLuwVR5IkSZKkUTaSM26AdwA79nm9qTb+3zLztg5xPsdY0ebnwJuB3cr7z8vnRwFdixwR8TTgQsaKNqcDLwNeDHwCWAqsD5wTETv1iLM/cBpV0eYu4P3Ai4A/B84tw3YHzo2IXv/cz2CsaHMpcGD5TkcAt1A9M5+JiCN75LIFcAFV0WY5cDywd3kdXz57JnBhGStJkiRJkjoYyRk3mXlrvzER8dZa95sdrm8NfLR0FwJ7Z+ay0r82Ir4PXA7sAhwbEV/PzFs6/KhjqGbVAHw0M0+oXbsqIi4FfgzMBE4EXtohlxnAyVRFlfuBl7T9rH+NiFOA91IVT/6iy3faBzisdC8ADsrMx9q+08+ALYHjI+K7mXlvh+/0N8CmpX1YZp5Tu/aTiFgIfKeM+Szwzg4xJEmSJEkaedFjBc/IKjNS/ptqydNSYNPMfKhtTKsQArBHZv60Q5zdgatK9+TMfF/b9acCvwU2AK4HdsjMFR3ifBV4d+nukpk/a7t+CFUhBOBjmTm/Q4yZwBJgQ+BXmbljhzELgFcDjwFzMnNJhzGHAt8u3WMy80tt1zcFbgPWAC7OzI77/ETEvwL7l5+1RWbe1WncREXELGAxwOLFi5k1a9Yww0uN0u8UJa3MU6UkSZI0GZYsWcLs2bNb3dmdfqeeiFFdKtXPy6iKNgDf7VC0CeB1pXtDp6INQPn8xtI9sNxXty9V0QbgG52KNsWZtfbBHa4f2GVsPZeHGCvu7BAR29SvlyVbLyvdS3o8YOdSzerplstrqYo2AF/vEqOe5xrlHkmSJEmS1GYkl0oN4G219hOWFAHPZqywc3mfWJdT7V8zC5gD1Jdp7dU2rpuFwIPAusCeHa634tyYmXf2yaU1c2dPVt5weTfGNjXumktmPhoRPwVeCewWEU/NzD92yKVnnLZre1Lt7SNJk26QGUrOypEkSVJTWLhpU2aeHFS6/w1c1mHYdrX2DX1C1q9vx8qFm4HiZObyiLiF6jSm+j2tfFvrgMabC136g8R5JdXzsw1wXYc49/UqImXmHRFxP7Beh1z6KkuhetlsvDElSZIkSWoaCzdP9HqqmS0AZ3U5xnt2rd1vzdriLvfV+w922eS3Pc5OwCYRsVZmPlI+nwW0lmANI5eJxLmurT9IjFac7TvkMojF/YdIkiRJkrRqc4+bJ+q3TArg6bX20j7xHqy1n9YlTr8YveIMO5dhxRnPd2qPIUmSJEmScMbNSsrym31L96eZeVOXoWvX2o/2CftIrb1Olzj9YvSKM+xchhVnPN+pPcYg+s3S2Qy4dgJxJUmSJElqDAs3K/sLxmYhfaPHuIdr7TX7xFyr1l7WJU6/GL3iDDuXYcSZOUCMepz2GH31O1btiQd4SZIkSZK06nGp1MreWt4fAf53j3EP1Nr9lvmsW2u3Lx9qxRlkqVC3OMPOZVhxxvOdBllWJUmSJEnSyLFwU0TELsCfle6FmfmHHsPrsz36nW5UX9LTvqFuK866EbHBgHHurm1MPBm5DCtOvxj1OG40LEmSJElSBxZuxtQ3Je61TApWPkVp2z5j69evn0iciJgBbNUpRmYuZazwMem5tF1fDtzcJc76EdH1SO6I2JzqKPBOuUiSJEmSJNzjBoCIeCpwaOneDfygzy23ArcDzwT26TN27/J+G7Co7doVtfY+wE+7xNiFsWVFV3a4fgXwZuC5EbFZZt7ZJU491/Y411JtKLxmGTe/U4CIWBPYvXVPZrZvQnwFY0vO9qH7krNeuUjStJozb0HP64vmz52iTCRJkjTqLNxU/hzYpLT/OTOX9xqcmRkR5wPvAbaNiN0z8wlFl4jYnbHZKednZrYNuQy4D1gfeHtEHN9hDMDhtfZ5Ha5/j6pw0xr7hKJLRMwE3li617WfmJWZD0TEv1H9Wbw8ImZ12QD4YMZmynTK5fvAP1LN5noH3Qs3h5f3FeUeSV30KyJIkiRJWn25VKpSXyb1zQHvOZFqqRDASRGx0pHWpX9S6S4v41dSZqt8pXS3A45pHxMRewBHlO7lmdnpiOvzgFtK+2MRsVWHMScAG9banXyxvM8ATomINdpy2Rj4QuneC3ytPUCZ7fOt0t0/It7QPiYiDgH2L92zeswQkiRJkiRppI184SYiNgQOKN1fZeb/G+S+MmOlVejYBbgyIt4UEbtExJuolv/sUq6fkJm/7hLqBKA1++X4iDg1IvaLiN0j4mPAD6kKKcuAD3bJ5Y/A+6lmr6xXcvmriNgtIvaPiO8C7y3DrwDO6hLnR8DZpfta4JKIeG35Tu+gWsq1Zbk+r8cGzp+gWnIG8O2ImB8Re5bXfOCfy7W7gU92iSFJkiRJ0shzqRS8CVirtAedbdPyCeAZwDuBFzBW9Kg7gx7FibJEaS5wEbANcFR51d0PvCUz/71HnIsi4mjgZGBTxmb71F0DHJSZj3WLQ/Vd1gNeDexXXnUrgM9m5qk9clkcEa+hWsK1GXBsedXdCRzYZTmWJEmSJEnCGTcwtpHuY4wt8RlIZq7IzCOAucD5VBsWP1rezwdenZlHZuaKPnFupir8HAsspFqG9BBwI/D3wE6ZeeEA+ZwO7AycDvwGeBj4HdUsm/cAL8nMe/rEWJaZc4G3AJcAvy3faTHVTJk9M/O4AXK5GtgR+BzwK2Bpef2yfLZDGSNJkiRJkrqIznvhSqu2iJhFOSZ98eLFzJo1a5ozkibOzYmbx1OlJEmS1G7JkiXMnj271Z09rBUmzriRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaE8VUqSpHEaZN8h98GRJEnSMDjjRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGco8bSZpGg+yVIkmSJGl0OeNGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ3kcuCRJk6DfUe+L5s+dokwkSZK0KnPGjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ82Y7gQkaXU2Z96C6U5BkiRJ0irMGTeSJEmSJEkNZeFGkiRJkiSpoVwqJUnSNBhkGd2i+XOnIBNJkiQ1mTNuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhZkx3ApIkqbM58xb0vL5o/twpykSSJEnTxRk3kiRJkiRJDWXhRpIkSZIkqaFcKiVJE9RvGYskSZIkPVnOuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaqgZ052AJEmamDnzFvQds2j+3CnIRJIkSZPFGTeSJEmSJEkNZeEGiIiNI+KjEXFlRNwZEY9ExO0RcXVEnBARewwQ41URcW5ELCn3Lyn9V40jj5kR8ZGIuCYifh8RSyPi+oj4YkRsOY4420fEVyPi5ohYFhF3R8SPI+LdETHwLKuIODQiLo6IOyLi4YhYFBFnRcTu44ixUUR8OiJ+ERH3RcT9pf3piNho0DiSJEmSJI2ikV8qFRGHAP8ItBcRNi+v3YBtgAO73B/AV4Gj2i5tARwEHBQRpwFHZ2b2yGMrYAHw3LZL25bXkRFxWGZe1Of7HAGcAqxV+3htYK/yOjwiDsjM3/WIsTZwDnBA26VnlddhEXFcZn62Ty67AudT/TnW7VReR0bE6zJzYa84kiRJkiSNqpGecRMRbwPOpira/Bb4NPAKYGdgLvB+4BLgjz3CfI6xos3PgTdTFXveXPqU612LHBHxNOBCxoo2pwMvA14MfAJYCqwPnBMRO/WIsz9wGlXR5q6S/4uAPwfOLcN2B86NiF7/7M9grGhzKVXRajfgCOAWqufmMxFxZI9ctgAuoCraLAeOB/Yur+PLZ88ELixjJUmSJElSm+gxCWS1FhHbURVW1gJ+ArwmM+/rMnbNzHy0w+dbA9dTzVxaCOydmctq12cClwO7UBUqts3MWzrEOQ74VOl+NDNPaLu+B/Dj8nMuzcyXdogxo+SyNXA/8ML2nxURpwDvLd23Z+Y3O8TZB7isdC8ADsrMx2rXNwZ+BmwJ/AF4Tmbe2yHOmcDbS/eNmXlO2/VDgO+U7tcz853tMZ6MiJgFLAZYvHgxs2bNGmZ4CRhsY1hpurk5sSRJ0tRYsmQJs2fPbnVnZ+aSYcQd5Rk3J1EVbe4BDu5WtAHoVLQpPsTYcrP31Ys25b6HgPeV7gzgg+0BIuKpwAdK93rgSx1+/lVUs2AA9ouInTvkchBV0Qbg850KRMBHqIotrXYnHy3vjwHvrRdtSi73AMeW7oZUs3BWEhGbAn9Ruhe3F21KnHOAi0v3beUeSZIkSZJUM5KFm4jYlmopEsDJpRgx3hgBvK50b8jMn3YaVz6/sXQPLPfV7QtsUNrfyMwVXX7kmbX2wR2uH9hlbD2Xhxib5bJDRGxTv16WbLX+XC7pUR08l2pWT7dcXgusUdpf7xKjnuca5R5JkiRJklQzkoUb4JBa+/HZIBGxYURsM+BpR8+m2oAYquVQvbSuzwLmtF3bq8O4ThYCD5b2nh2ut+LcmJl3DpBLpzi7MbapcddcygykVqFqtzJrqFMuPeP0yUWSJEmSpJE3qoWb1nHW9wHXR8RbIuIXwO+Bm4B7IuI3EfGpMgulk+1q7Rv6/Lz69e3arg0UJzOXU20M/IQYJcfWJi6Tnkvb9RlUp251inNfryJSZt7B2Myd9lwkSZIkSRp5o3oc+J+V90VUe938ZYcxzwaOA94QEftn5u1t12fX2v02HFrc5b56/8FOm/x2iLMTsElErJWZj5TPZwGtJVjDyGUica7rEGeQjZgWA9t3yKWnsvlwL5uNJ57UiZsPS5IkSZpuo1q4+ZPyvi3wPOBeYB5je7fsCHyG6hjtHaiO4d6rbf+Zp9faS/v8vAdr7fYZPK04/WJ0itMq3Aw7l2HFGc936jazqZvF/YdIkiRJkrRqG9WlUuuW97WoTk/688w8NTPvzsxHMnMhcADwgzLuxTxxE961a+1up061PFJrr9MlTr8YveIMO5dhxRnPd2qPIUmSJEnSyBvVws3DtfY5nU6EKrNr6kdmv7lHjDX7/Ly1au1lbddacfrF6BVn2LkMK854vlN7jH5m93ntOs54kiRJkiQ1zqgulXqAsVk3P+g2KDP/MyJuozo9qr0Q8ECt3W+Zz7q1dvvyoVacQZYKdYsz7FyGEWfmADHqcQZZVvW4HkeVA/DEU9claTT126tp0fy5U5SJJEmSJmJUZ9zU90cZdBPeZ7R9Xr+v30a59Y132/dmacVZNyI2GDDO3bWNiScjl2HF6RejHsc9ayRJkiRJajOqhZv/rLXX6DO2dX152+f1U5S27ROjfv36icSJiBnAVp1iZOZSxgofk55L2/XlwM1d4qwfEV1Pd4qIzYH1uuQiSZIkSdLIG9XCzY9r7a26jqo8p7zf1vb5rUDriPB9+sTYuxZjUdu1K2rtXnF2YWxZ0ZUdrrfiPLdXsaTtZ7THuZaxDYW75hIRawK7t+7JzPZNiAf9Tr1ykSRJkiRp5I1q4eb7wB9Lu/20qMdFxD7ARqX7k/q1zEzg/NLdNiJ2p4PyeWt2yvnlvrrLgPtK++3RfXOWw2vt8zpc/16XsfVcZgJvLN3rMvOm+vXMfAD4t9J9eUR0W+p0MGMzZTrl8n2gdXT6O7rEqOe5otwjSZIkSZJqRrJwk5m/A75Wuq+IiEPbx0TE04ETax+d2iHUiYwtoTopIlY60rr0Tyrd5W3xWrk8CnyldLcDjumQyx7AEaV7eWZe2yGX84BbSvtjEdFpJtEJwIa1didfLO8zgFMiYqWlZBGxMfCF0r2XsT/Hx2XmncC3Snf/iHhD+5iIOATYv3TPKvdIkiRJkqSakSzcFJ8C/ru0z4qIkyJiv4jYOSIOB64Bnl+u/2OnYkmZsdIqdOwCXBkRb4qIXSLiTVTLf3Yp10/IzF93yeUEoDX75fiIOLXksntEfAz4IVUhZRnwwU4BMvOPwPupZq+sV3L5q4jYLSL2j4jvAu8tw68AzuoS50fA2aX7WuCSiHht+U7vAH4KbFmuz8vMP3T5Tp8A7i7tb0fE/IjYs7zmA/9crt0NfLJLDEmSJEmSRlo8ceXO6IiI7aiW6GzdY9g/AUeXwkinGE8BTgfe2SPGGcBRmbmi24CI2Bq4CNimy5D7gbdk5oU9fg4R8S7gZGDNLkOuAeZm5j09YqwDfBd4dZchK4DPZuZxfXJ5EdUSrm577twJHJiZV/eKMxFlmddigMWLFzNr1iAHXEkr63eMsrQ68DhwSZKk4ViyZAmzZz9+APPszOx3ivVARnnGDZl5PdWsmo8AVwO/p9qcdwnwv4GXZuYR3Yo2JcaKzDwCmEu1583tJcbtpf/qzDyyV9GmxLkZeAFwLLCQahnSQ8CNwN8DO/Ur2pQ4pwM7UxWTfgM8DPyOapbNe4CX9CralBjLMnMu8BbgEuC35Tstppops2e/ok2JczWwI/A54FfA0vL6Zflsh8ko2kiSJEmStLoY6Rk3Wn0540bD4IwbjQJn3EiSJA2HM24kSZIkSZJGjIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDzZjuBCRJ0vTpd3qap05JkiRNL2fcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaiiPA5c0svodgyxJkiRJ080ZN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqqBnTnYAkSWquOfMW9B2zaP7cKchEkiRpNDnjRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaasZ0JyBJklZtc+Yt6Hl90fy5U5SJJEnS6scZN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lHvcSFot9dtzQ5IkSZJWBc64kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaFmTHcCkiRp9TZn3oK+YxbNnzsFmUiSJK16nHEjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNNbKFm4jIAV+XDRDrVRFxbkQsiYhHyvu5EfGqceQzMyI+EhHXRMTvI2JpRFwfEV+MiC3HEWf7iPhqRNwcEcsi4u6I+HFEvDsiBj5FLCIOjYiLI+KOiHg4IhZFxFkRsfs4YmwUEZ+OiF9ExH0RcX9pfzoiNho0jiRJkiRJo8rjwJ+EiAjgq8BRbZe2AA4CDoqI04CjMzN7xNkKWAA8t+3StuV1ZEQclpkX9cnnCOAUYK3ax2sDe5XX4RFxQGb+rkeMtYFzgAPaLj2rvA6LiOMy87N9ctkVOB/YvO3STuV1ZES8LjMX9oojSZIkSdIoG9kZNzX/COzY4/WOHvd+jrGizc+BNwO7lfefl8+PAroWOSLiacCFjBVtTgdeBrwY+ASwFFgfOCciduoRZ3/gNKqizV3A+4EXAX8OnFuG7Q6cGxG9/rmfwVjR5lLgwPKdjgBuoXpmPhMRR/bIZQvgAqqizXLgeGDv8jq+fPZM4MIyVpIkSZIkdeCMG/htZv5qvDdFxNbAR0t3IbB3Zi4r/Wsj4vvA5cAuwLER8fXMvKVDqGOoZtUAfDQzT6hduyoiLgV+DMwETgRe2iGXGcDJVEWV+4GXtP2sf42IU4D3UhVP/gL4Zoc4+wCHle4FwEGZ+Vjbd/oZsCVwfER8NzPv7fCd/gbYtLQPy8xzatd+EhELge+UMZ8F3tkhhiRJkiRJI88ZNxP3IcYKX++rFW0AyMyHgPeV7gzgg+0BIuKpwAdK93rgS+1jMvMqqlkwAPtFxM4dcjkI2Lq0P9+lQPQR4A+1dietQtRjwHtrRZtWLvcAx5buhlSzcFYSEZtSFYYALm4r2rTinANcXLpvK/dIkiRJkqQ2Fm4moOxt87rSvSEzf9ppXPn8xtI9sNxXty+wQWl/IzNXdPmRZ9baB3e4fmCXsfVcHqKa5QKwQ0RsU79elmy9rHQvycwlXXI5l2pWT7dcXgusUdpf7xKjnuca5R5JkiRJktTGws3EPJtqA2KolkP10ro+C5jTdm2vDuM6WQg8WNp7drjeinNjZt45QC6d4uzG2KbGXXPJzEeBVqFqtzJrqFMuPeP0yUWSJEmSJGHhBuCQiLixHJ39QET8OiK+ERH79bhnu1r7hj7x69e3a7s2UJzMXE61MfATYpSZMrOmKpe26zOAbdquteLc16uIlJl3MDZzpz0XSZIkSZKEmxMD/Flbf+vyeltEfA84PDPvaxszu9butqSoZXGX++r9B7ts8tseZydgk4hYKzMfKZ/PAlpLsIaRy0TiXNchTr8YrTjbd8ilr4iY1WfIZuONKUmSJElS04xy4eYh4PvAv1HNIFkKbALsAxwNbES1d8z5EfGKzPxj7d6n19pL+/ycB2vtp7Vda8XpF6NTnFbhZti5DCvOeL5Te4xBLO4/RJK0qpgzb0HP64vmz52iTCRJkppllAs3W3SZ5XJJRJwE/AB4AVUh5z3AV2pj1q61H+3zcx6ptddpu9aK0y9GrzjDzmVYccbzndpjSJIkSZIkRrhw02tpUmbeFRFvoDqie02qY73rhZuHa+01+/yotWrtZW3XWnH6xegVZ9i5DCPOzAFi1OO0xxhEv+VVmwHXTiCuJEmSJEmNMbKFm34y8zcRcQkwF9g6Ip6ZmbeXyw/UhvZb5rNurd2+fKgVZ5ClQt3iDDuXYcSZOUCMepxBllWtpMdx5QA88eR1rW76LauQJEmSpNWBp0r1Vt90d4tau1406LdJbn1mSPu+LK0460bEBgPGubu2MfFk5DKsOP1i1OO4X40kSZIkSR1YuOmt27SNekFn2z4x6tevn0iciJgBbNUpRmYuZazwMem5tF1fDtzcJc76EdH1ZKeI2BxYr0sukiRJkiQJCzf91I8Kv73WvrXW36dPjL3L+23AorZrV9TaveLswtiyois7XG/FeW6vYknbz2iPcy1jGwp3zSUi1gR2b92Tme2bEA/6nXrlIkmSJEmSsHDTVUQ8B3hF6f4mM29rXcvMBM4v3W0jYvf2+0uM3RmbnXJ+ua/uMuC+0n57dN+Y5fBa+7wO17/XZWw9l5nAG0v3usy8qX49Mx+gOhod4OUR0W2p08GMzZTplMv3gRWl/Y4uMep5rij3SJIkSZKkNiNZuImI15TlR92ubwp8F3hq+eiUDsNOpFoqBHBSRKx0pHXpn1S6y8v4lZTZKq3TqrYDjumQyx7AEaV7eWZ2OinpPOCW0v5YRGzVYcwJwIa1didfLO8zgFMiYo22XDYGvlC69wJfaw+QmXcC3yrd/cvpXCuJiEOA/Uv3rHKPJEmSJElqM6qnSp0EPDUi/gW4imoJ0zJgY2Bf4GhgozL2CjoUbjLzpoj4IjCPainTlRHxBaoCylbAscALyvATMvPXXXI5AXgT8KfA8RGxNXB2yWc/4ONU/5yWAR/sFCAz/xgR7wcuoJoNc2VEfA64hqpY8y7g9bXvc1aXOD+KiLOBQ4HXApdExIlUy8J2BD4BbFmGz8vMP3T5Tp8AXgVsAnw7InYBLizXDgA+XNp3A5/sEkOSJEmSpJEXT1y9s/qLiEXAswYY+i/AkZl5b5c4TwFOB97ZI8YZwFGZuaLbgFKsuQjYpsuQ+4G3ZOaFXa634rwLOBlYs8uQa4C5mXlPjxjrUM02enWXISuAz2bmcX1yeRHVEq5ue+7cCRyYmVf3ijNRZanXYoDFixcza9Ygh1xpVeJx4NJoWTR/7nSnIEmS1NOSJUuYPfvxQ5hnZ+aSXuMHNaozbt5OtTnuHsBzqGbarAe0Tmj6v8A3MvOqXkFKMeaIMnPnKGDXEuseqs1+T83MH/RLJjNvjogXAH8JHAJsTVV8WUxV0PlyZv7XAHFOj4irgPcDLwOeCTxIdWrTt4CvZebyHiHIzGXA3Ig4jGofmucBGwB3AT8BTu7351LiXB0ROwIfAA4E5pRLt1LtD3RiZv6uXxxJkiRJkkbZSM640erPGTerP2fcSKPFGTeSJKnpJmvGzUhuTixJkiRJkrQqsHAjSZIkSZLUUBZuJEmSJEmSGmpUNyeWJEmrkEH2tXIfHEmStDpyxo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaasZ0JyBJkjQMc+Yt6Hl90fy5U5SJJEnS8DjjRpIkSZIkqaEs3EiSJEmSJDWUS6UkNU6/5Q6SJEmSNCqccSNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhpqxnQnIEmSNBXmzFvQd8yi+XOnIBNJkqTBOeNGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqqBnTnYAkSVJTzJm3oOf1RfPnTlEmkiRJFWfcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaqhGbE0fEWsCewMbArZl5zTSnJEmSJEmSNO0mvXATEc8C/rJ0/zYz7227vjvwXWDz2mf/D3h9Zv73ZOcnSZIkSZLUVFOxVOog4BjgpR2KNk8HvkdVtInaa2dgQURM+YygiDg+IrL22neAe14VEedGxJKIeKS8nxsRrxrHz50ZER+JiGsi4vcRsTQiro+IL0bEluOIs31EfDUibo6IZRFxd0T8OCLePZ4/z4g4NCIujog7IuLhiFgUEWeVQtugMTaKiE9HxC8i4r6IuL+0Px0RGw0aR5IkSZKkUTUVhZtXAElVoGl3FPCM0v4K8DrgH0r/z4C3T3ZydRHxPOBD4xgfEXEq8AOqAtUWwJrl/SDgBxFxakREnzhbAf8POB7YFdgQWBfYFvgw8B8R8eoB8jkC+BnwbmArYG2q5Wd7AV8FftKvYBIRa0fEBcC3gVcCmwFrAc8C/gK4MiL+xwC57Ar8EvifwE7AesDTS/t/lu+0S784kiRJkiSNsqmY0fKc8v6zDtfeSFXUOS8zP1g+uyAiNgEOAd4AnDHpGQIR8RTgdKo/k98yVlDq5XNUxSeAn1MVXm6hKpp8FHhBuX438MkuP/dpwIXAc8tHpwNnA8uA/YCPAesD50TEHpn5H13i7A+cRlWMuwv4G+Bq4E+AdwEHA7sD50bEfpm5ost3OgM4oLQvBb4M3A7sCHy8fLfPRMQdmfm1LrlsAVwAbAosB/6ufEdK7L8GnglcGBE7Z+ZtXXLRamrOvAXTnYIkSZIkrRKmYsZNqwByV/3DiFgPeGHpfr3tnrPL+/MmMa9276ea7XIDAxSLImJrquIMwELgJZl5dmZem5lnU222vLBcP7bMqunkGKqZNQAfzcyjMvNHmXlVZv4t1ayX5cBM4MQuucwATqb653l/yeWkzLwmM/81M1/P2EymvalmznSKsw9wWOleALwiM88v3+mfqAo/rX2Hjo+IDbp8p7+hKtoAHJaZx2bmT8rr2NrP2BT4bJcYkiRJkiSNvKko3Dy9vK/R9vlLymePAZe1XVtc3v9k8tIaExGzGSsgvAd4dIDbPsTYjKX3Zeay+sXMfAh4X+nOAD7Y4ec+FfhA6V4PfKl9TGZexVghab+I2LlDLgcBW5f25zPzlg5jPgL8odbupFWIegx4b2Y+1pbLPcCxpbshcER7gIjYlLHC0MWZeU77mPLZxaX7tnKPJEmSJElqMxWFm/vK+zPbPt+3vP8iMx/scu/Dk5LRE/0D8DTgG5l5Wb/BZc+a15XuDZn5007jyuc3lu6BHfa62RfYoLS/0WP50pm19sEdrh/YZWw9l4eA75TuDhGxTf16WbL1stK9JDOXdMnlXKpZPd1yeS1jRbr2mVR1rTzXKPdIkiRJkqQ2U1G4+VV5P6j1QUSswdj+Npd2uGeL8n5Xh2tDFRFvpNp35fd0n4nS7tmM5Xh5n7Gt67OAOW3X9uowrpOFQKu4tWeH6604N2bmnQPk0inOblSbEPfMJTMfBVqFqt3KrKFOufSM0ycXSZIkSZLE1GxOfB6wD/DWiLgL+AnwVqpTipKxWSB1rdOG/rvDtaEpe7R8uXSPzcy7B7x1u1r7hj5j69e3A24db5zMXB4Rt1CdyFS/pzVTZtYEc6FLf5A4r6R6frYBrusQ575eRaTMvCMi7qc6bao9l74iYlafIZuNN6YkSZIkSU0zFYWbU6mOp96OaiPeY2rXLsjMhR3uOYiqqHPZJOd2PNUv+P+X8Z1eNbvW7rakqGVxrT277Vqr/2Bm3jtAnJ2ATSJircx8pHw+C2gtwRpGLhOJc11bf5AYrTjbd8hlEIv7D5EkSZIkadU26YWbzHwkIl5GderRa4CnUm3++7+Bv2ofHxF7A39GVbi5ZLLyiog9gSOpTmw6OjNzHLc/vdZe2mdsff+ep3WJ0y9Gpzitws2wcxlWnPF8p/YYkiQ10px5C/qOWTR/7hRkIkmSRsVUzLihLJl5Q0SsRXVS1O/KXimdLAb2K+1Os3GetIhYEziNaqbK32fmL8cZYu1au98JVI/U2ut0iTPIKVbd4gw7l2HFGc93ao8xiH6zdDYDrp1AXEmSJEmSGmNKCjctZXnPHX3G3MrK+8BMho9TLd36b+DTE7i/ftrVmn3GrlVrL2u71orTL0avOMPOZRhxZg4Qox6nPUZfPU69AuCJB3hJkiRJkrTqmfRTpSLinyLijIjYfBz3bNK6bxLy2Rb4WOm+r8dR5L08UGv3W+azbq3dvnyoFWeQpULd4gw7l2HFGc93GmRZlSRJkiRJI2cqZtwcTrVfzZfoM9umZr3afUcMOZ8PUc0G+Q0wMyIO7TBmh1r7pRHROqHoglLoqc/26He6UX1JT/uGukuAFwHrRsQGfTYobsW5u7YxcSvGsHKpx+m1TK1fnE0HyKUex42GJUmSJEnqYEqXSjVEa3nOc4BvDzD+f9Taz6baULd+itK2fe6vX7++7dp1wOtr437aKUBEzAC26hQjM5dGxGKqIsiTzaXTuF5xlgM3d4izM7B+RGzW7UjwMgNrvS65SJIkSZIkpmCp1AS1Nrh9pOeo6XMrcHtp79Nn7N7l/TZgUdu1K2rtXnF2YWxZ0ZUdrrfiPLc2O6iT+s9oj3MtYxsKd82lbOy8e+ueDptMD/qdeuUiSZIkSZJobuHmJeX9rmEHzszDMzN6vVh5w+L9atcWlRgJnF+ubxsRu9NB+bw1O+X8DkeOXwbcV9pvj+476h5ea5/X4fr3uoyt5zITeGPpXpeZN9WvZ+YDwL+V7ssjottSp4MZmynTKZfvAytK+x1dYtTzXFHukSRJkiRJbYa+VCoi/meXS++NiN/2uX0tqiVBr6Xa36bJMzFOBN5F9Wd4UkTsnZmPn44UEesAJ5Xu8jJ+JZn5aER8hWo51nbAMcAJ9TERsQdj+/xcnpmdjrg+D7iF6s/uYxFxTmbe0jbmBGDDWruTLwJ/Xr7TKRFxcGY+VstlY+ALpXsv8LUO3+nOiPgW8FZg/4h4Q2Z+t+07HQLsX7pndVtOJUmSJEnSqJuMPW6Ooyq61AXwnnHECKpjpbsVGKZdZt4UEV8E5lEtZboyIr7AWAHlWOAFZfgJmfnrLqFOAN4E/ClwfERsDZxNdUT2flRHl88o/Q92yeWPEfF+4AKq2TBXRsTngGuoijXvYmwvnSuAs7rE+VFEnA0cSlU8uyQiTqRaFrYj8AlgyzJ8Xmb+oct3+gTwKmAT4NsRsQtwYbl2APDh0r4b+GSXGJIkSZIkjbzJ2py4vuQnO3zWzcNUJ0/9X+CLmfmLYSc2ZJ8AngG8k6pIc3aHMWfQoziRmQ9ExFzgImAb4KjyqrsfeEtm/nuPOBdFxNHAyVSnOp3UYdg1wEH1WTQdvJOq+PNqqsLRfm3XVwCfzcxTe+SyOCJeQ7WEazOqItaxbcPuBA7MzCVIkiRJkqSOhl64ycyV9s2JiBVUxZsdMvO6znetmjJzBXBERPwLVbFlV2Bj4B6qzX5PzcwfDBDn5oh4AfCXwCHA1lRHli+mKuh8OTP/a4A4p0fEVcD7gZcBz6Q6Bet64FvA1zJzeZ8Yy4C5EXEY1T40zwM2oNpv6CfAyZl51QC5XB0ROwIfAA4E5pRLt1LtD3RiZv6uXxxJkiRJkkZZPHG/3CH/gIhFVIWbV2Rm+9HR0qQomysvBli8eDGzZnXba1nTYc68BdOdgiRNmkXz5053CpIkaRosWbKE2bNnt7qzh7XCZLKWSj0uM+dM9s+QJEmSJElaHU164UaSJGmU9JtV6IwcSZI0HlNauCknJr0N2INq09p1gFfVl1BFxA5UJxc9mJmXT2V+kiRJkiRJTTIlhZuIeArwBarjrJ/C2AlTSbUJb91sqqOjl0fEszPztqnIUZIkSZIkqWme0n/IUJwK/DWwBnA78N1uA8spTL8pY98wJdlJkiRJkiQ10KQXbiJiX+CI0v1bYE5mvrHPbedQzcrZb/IykyRJkiRJarapWCp1dHm/KDM/OeA915T37SchH0mSJEmSpFXCVCyV2oNqL5szxnFP66zzzYafjiRJkiRJ0qphKgo3zyjvt47jnuXl/alDzkWSJEmSJGmVMRWFm2XlfeY47tmyvP9hyLlIkiRJkiStMqaicNOaafOCcdxzQHm/bsi5SJIkSZIkrTKmonDzQ6oToo6KiL4/LyJ2Bt5KtS/Ov05ybpIkSZIkSY01FYWbk6mWS+0InB4RXfetiYjXUxVr1gTuB06bgvwkSZIkSZIaadKPA8/M2yLi/cDpwOHAKyPigtqQIyJiJvBy4DlUs3MSOCoz75vs/CRJkiRJkppq0gs3AJl5RkQk8BVgC+DdVMUZgA+W9yjvjwBHZ+Y5U5GbJEmSJElSU01J4QYgM/8pIn5IVah5LbB125DbgO8DJ2TmoqnKS5IkaSrNmbeg5/VF8+dOUSaSJGlVMGWFG4DMXAIcAxwTEesBzwDWAH6XmfdMZS6SJEmSJElNN6WFm7rMvJ9qA2JJkiRJkiR1MG2FG0mrp35LACRJkiRJgxtq4SYi/ucQwyXwEPAH4AZgYWY+OsT4kiRJkiRJjTbsGTfHMXZa1LDdHxGnA8dl5kOT9DMkSZIkSZIa4ymTEDMm6bU+8GHghxHx1EnIW5IkSZIkqVGGPeNmvyHHWxPYENgeeDPVEeJ7AEcDJw35Z0mSJEmSJDXKUAs3mXn5MOPVRcTfAP8GvAR4ExZuJEmSJEnSam4ylkpNirIx8Qmlu9105iJJkiRJkjQVVpnCTbGovK8/nUlIkiRJkiRNhWHvcTPZlgAfmu4kJEmSJEmSpsIqVbjJzN8DX57uPCRJkiRJkqbCqrZUSpIkSZIkaWRYuJEkSZIkSWqoVWqplCRJ0upuzrwFfccsmj93CjKRJElN4IwbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkPNmO4EJEmSND5z5i3oeX3R/LlTlIkkSZpszriRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNdSM6U5AkiRJwzVn3oK+YxbNnzsFmUiSpCdrJAs3EbEe8GpgV2AXYAtgE2Ad4F7gOuAi4IzM/N0A8V4FHAXsVuLcDVwDnJaZ/zpgTjOBvwQOAbYG1gQWAwuAr2Tmfw8YZ3vgfcDLy/daClwPfKt8n+UDxjkUeAewE7AhcCfwE+CUzPzpgDE2At4PHAjMAQK4Ffhe+U59/2zVPIP8MiBJkiRJGo7IzOnOYcpFxMuBSwYYeg/wF5l5cZc4AXyVqmjTzWnA0dnjDzoitqIq0Dy3y5D7gMMy86JeyUbEEcApwFpdhvwUOKBXwSQi1gbOAQ7oMmQFcFxmfrZPLrsC5wObdxlyO/C6zFzYK85ERcQsqsIXixcvZtasWZPxY0aShRtJWj0440aSpOFasmQJs2fPbnVnZ+aSYcQd5T1uFgPfBD4AHAzsAbwEeBNV4eIxYGPg+xGxU5cYn2OsaPNz4M1Us27eXPqU612LHBHxNOBCxoo2pwMvA14MfIJqxsz6wDk98iAi9qcqEq0F3EU10+VFwJ8D55ZhuwPnRkSvf+5nMFa0uZRqtsxuwBHALVTPzGci4sgeuWwBXEBVtFkOHA/sXV7Hl8+eCVxYxkqSJEmSpA5GdcbNGpn5WJ8xBwLnle65mfn6tutbUy1BmgEsBPbOzGW16zOBy6mWYi0Hts3MWzr8nOOAT5XuRzPzhLbrewA/Lj/n0sx8aYcYM0ouWwP3Ay9s/1kRcQrw3tJ9e2Z+s0OcfYDLSvcC4KD6n1NEbAz8DNgS+APwnMy8t0OcM4G3l+4bM/OctuuHAN8p3a9n5jvbYzxZzriZPM64kaTVgzNuJEkaLmfcDFG/ok0Z8z3ghtLdu8OQDzG2R9D76kWbcv9DVHvNUMZ9sD1ARDyVasYPVIWXL3XI4yqqWTAA+0XEzh1yOYiqaAPw+U4FIuAjVMWWVruTj5b3x4D3tv85ZeY9wLGluyHVLJyVRMSmwF+U7sXtRZsS5xygtfzsbeUeSZIkSZLUZiQLN+PwYHlfu/5h2dvmdaV7Q7fNesvnN5bugeW+un2BDUr7G5m5okseZ9baB3e4fmCXsfVcHmJslssOEbFN/XpZsvWy0r2kR2XwXKpZPd1yeS2wRml/vUuMep5rlHskSZIkSVIbCzddRMR2wPNL94a2y8+mOrEJquVQvbSuz6I6Walurw7jOlnIWBFpzw7XW3FuzMw7B8ilU5zdGNvUuGsumfko1SbHALuVWUOdcukZp08ukiRJkiSJET0OvJuyL80WwGuolg21Zo58uW3odrV2e1GnXf36dlTHYY8rTmYuj4hbqI7mrt/TminT2sBlvLnQpT9InFdSPT/bUB2f3h7nvl5FpMy8IyLuB9brkEtfZQ+bXjYbb0xJkiRJkppm5As3EXE4vZf0fBH4Vttns2vtfpsNLe5yX73/YKdNfjvE2QnYJCLWysxHyuezgNYSrGHkMpE417X1B4nRirN9h1wGsbj/EEmSJEmSVm0jX7jp4d+BozPz6g7Xnl5rL+0T58Fa+2ld4vSL0SlOq3Az7FyGFWc836k9hiRJkiRJwsINwPeo9pABWAfYCngj1UlN34qID2bmhW331DcrfrRP/Edq7XW6xOkXo1ecYecyrDjj+U7tMQbRb5bOZsC1E4grSZIkSVJjjHzhpixRurf20bXA2RHxVuAbwPkRcURmnlkb83CtvWafH7FWrb2s7VorTr8YveIMO5dhxJk5QIx6nPYYffU49QqAJx7gJUmSJEnSqmfkCzfdZOZZEXEA1eybkyPi/Mz8Q7n8QG1ov2U+69ba7cuHWnEGWSrULc6wcxlGnJkDxKjHGWRZlSRJGqI58xb0vL5o/twpykSSJPXiceC9nV/e1wX+vPZ5fbZHv9ON6kt62jfUbcVZNyI2GDDO3bWNiScjl2HF6RejHseNhiVJkiRJ6sDCTW9319rPqrXrpyht2ydG/fr1bdcGihMRM6j23nlCjMxcyljhY9Jzabu+HLi5S5z1I6LrkdwRsTnVUeCdcpEkSZIkSVi46WeLWru+nOdW4PbS3qdPjL3L+23AorZrV9TaveLswtiyois7XG/FeW6vYknbz2iPcy1jGwp3zSUi1gR2b92Tme2bEA/6nXrlIkmSJEmSsHDTzyG19i9bjcxMxpZRbRsRu9NB+bw1O+X8cl/dZcB9pf326L6j7uG19nkdrn+vy9h6LjOp9usBuC4zb6pfz8wHgH8r3ZdHRLelTgczNlOmUy7fB1aU9ju6xKjnuaLcI0mSJEmS2oxk4SYiDo+ItfuM+RDw6tJdxMozSQBOpFoqBHBSRKx0pHXpn1S6y8v4lZTZKl8p3e2AYzrksQdwROlenpmdjrg+D7iltD8WEVt1GHMCsGGt3ckXy/sM4JSIWKMtl42BL5TuvcDX2gNk5p3At0p3/4h4Q/uYiDgE2L90zyr3SJIkSZKkNiNZuAGOA26LiNMi4m0R8ZKIeF5E7BkR74mIK4C/K2MfBd6VmcvrAcqMlVahYxfgyoh4U0TsEhFvolr+s0u5fkJm/rpLLicArdkvx0fEqRGxX0TsHhEfA35IVUhZBnywU4DM/CPwfqrZK+uVXP4qInaLiP0j4rvAe8vwK4CzusT5EXB26b4WuCQiXlu+0zuAnwJbluvzaqdstfsEY/sDfTsi5pc/2z0jYj7wz+Xa3cAnu8SQJEmSJGnkxRNX76z+ImIRK2823M0S4J2ZeUmXOE8BTgfe2SPGGcBRmbmi24CI2Bq4CNimy5D7gbdk5oW9ko2IdwEnA2t2GXINMDcz7+kRYx3gu4zNNmq3AvhsZh7XJ5cXUS3h6rbnzp3AgZl5da84E1WWei0GWLx4MbNmDXLIlQbR7/hYSdLqwePAJUkanyVLljB79uOHMM/OzCW9xg9qxjCCrIJeBrwc2I9qidKmwEbAw8BdwL8DFwLfycyHugUpxZgjIuJfgKOAXYGNgXuoNvs9NTN/0C+ZzLw5Il4A/CXVvjpbUxVfFlMVdL6cmf81QJzTI+Iqqtk3LwOeCTxIdWrTt4Cvtc8c6hBjGTA3Ig6j2ofmecAGVH8uPwFOzsyrBsjl6ojYEfgAcCAwp1y6lWp/oBMz83f94kiSJEmSNMpGcsaNVn/OuJk8zriRpNHgjBtJksZnsmbcjOoeN5IkSZIkSY1n4UaSJEmSJKmhLNxIkiRJkiQ11KhuTixJkqQeBtnTzH1wJEmafM64kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaFmTHcCkiRJWjXNmbeg5/VF8+dOUSaSJK2+nHEjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkO5x42kx/Xbq0CSJEmSNLWccSNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhPA5ckiRJk2LOvAV9xyyaP3cKMpEkadXljBtJkiRJkqSGsnAjSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaasZ0JyBJkqTRNWfegp7XF82fO0WZSJLUTM64kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqqBnTnYAkSZLUzZx5C/qOWTR/7hRkIknS9HDGjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQI1u4iYgXRsTHI+IHEbE4Ih6JiKURcVNEnBkRe40z3qsi4tyIWFJiLSn9V40jxsyI+EhEXBMRvy/5XB8RX4yILccRZ/uI+GpE3BwRyyLi7oj4cUS8OyJmjCPOoRFxcUTcEREPR8SiiDgrInYfR4yNIuLTEfGLiLgvIu4v7U9HxEaDxpEkSZIkaRRFZk53DlMuIi4H9h5g6FnAkZn5aI9YAXwVOKpHnNOAo7PHH3ZEbAUsAJ7bZch9wGGZeVGvhCPiCOAUYK0uQ34KHJCZv+sRY23gHOCALkNWAMdl5mf75LIrcD6weZchtwOvy8yFveJMRETMAhYDLF68mFmzZg37R6yW5sxbMN0pSJI0bovmz53uFCRJYsmSJcyePbvVnZ2ZS4YRd+DZF6uZLcr77VQFip8A/w2sAewBfLiMeSvVn9FhPWJ9jrGizc+B44FbgK2AjwIvKNfvBj7ZKUBEPA24kLGizenA2cAyYD/gY8D6wDkRsUdm/keXOPtTFYmeAtwF/A1wNfAnwLuAg4HdgXMjYr/MXNHlO53BWNHmUuDLVH9WOwIfL9/tMxFxR2Z+rUsuWwAXAJsCy4G/K9+REvuvgWcCF0bEzpl5W5dcJEmSeur3Fw8WdiRJq7JRnXFzIfBN4F8y87EO1zcGrgT+tHy0d2b+pMO4rYHrqYo7C8u4ZbXrM4HLgV2oihfbZuYtHeIcB3yqdD+amSe0Xd8D+HH5OZdm5ks7xJhRctkauB94YfvPiohTgPeW7tsz85sd4uwDXFa6FwAH1f+Myp/Nz4AtgT8Az8nMezvEORN4e+m+MTPPabt+CPCd0v16Zr6zPcaT4YybiXHGjSRpdWThRpI0FSZrxs1I7nGTmQdk5nc6FW3K9XuoZt20vKFLqA8xNmvpffWiTYnzEPC+0p0BfLA9QEQ8FfhA6V4PfKlDPldRzYIB2C8idu6Qy0FURRuAz3cqEAEfoSq2tNqdfLS8Pwa8t/3PqPzZHFu6GwJHtAeIiE2Bvyjdi9uLNiXOOcDFpfu2co8kSZIkSaoZycLNgC6rtbdqv1j2tnld6d6QmT/tFKR8fmPpHljuq9sX2KC0v9Fj+dKZtfbBHa4f2GVsPZeHGJvlskNEbFO/XpZsvax0L+lRHTyXalZPt1xeS7XsDODrXWLU81yj3CNJkiRJkmos3HS3Zq3dqZjybMb2yrm8T6zW9VnAnLZre3UY18lC4MHS3rPD9VacGzPzzgFy6RRnN8Y2Ne6aS9msuVWo2q3MGuqUS884fXKRJEmSJGnkjermxIPYp9a+ocP17fpcp8v17YBbxxsnM5dHxC3ATm33tGbKtDZxGW8udOkPEueVVM/QNsB1HeLc16uIlJl3RMT9wHodctEkcA8bSZIkSVq1WLjpICKeAsyrffSdDsNm19r9Nhxa3OW+ev/BTpv8doizE7BJRKyVmY+Uz2cBrSVYw8hlInGua+sPEqMVZ/sOufRUNh/uZbPxxJMkSZIkqYks3HT2IaplQwDnZebCDmOeXmsv7RPvwVr7aV3i9IvRKU6rcDPsXIYVZzzfqT1GP4v7D5EkSZIkadXmHjdtynHY80v3t8B7ugxdu9Z+tE/YR2rtdbrE6RejV5xh5zKsOOP5Tu0xJEmSJEkaec64qYmI7YHzqP5cHgHemJl3dRn+cK29ZpcxLWvV2svarrXi9IvRK86wcxlGnJkDxKjHaY/RT7+lVZsB144zpiRJkiRJjWLhpoiIZwM/BDYEHgPenJm9TkR6oNbut8xn3Vq7fflQK84gS4W6xRl2LsOIM3OAGPU4gyyrelyPo8oBeOKp65IkaVT125x/0fy5U5SJJEnj51IpICKeCfwf4JlAAu/MzPP63FYvHPTbKLc+O6R9b5ZWnHUjYoMB49xd25h4MnIZVpx+Mepx3LNGkiRJkqQ2I1+4iYiNgUuA55SP3peZ3xzg1vopStv2GVu/fv1E4kTEDGCrTjEycyljhY9Jz6Xt+nLg5i5x1o+Irqc7RcTmVEeBd8pFkiRJkqSRN9KFm4hYH7gY+LPy0bzMPGXA228Fbi/tffqM3bu83wYsart2Ra3dK84ujC0rurLD9Vac5/YqlrT9jPY41zK2oXDXXCJiTWD31j2Z2b4J8aDfqVcukiRJkiSNvJEt3ETETGAB8MLy0d9k5hcGvT8zEzi/dLeNiN07jSuft2annF/uq7sMuK+03x7dN2c5vNbutIzre13G1nOZCbyxdK/LzJvq1zPzAeDfSvflEdFtqdPBjM2U6ZTL94EVpf2OLjHqea4o90iSJEmSpJqRLNyUGSPnAS8pH305Mz85gVAnUi0VAjgpIlY60rr0Tyrd5WX8Sspsla+U7nbAMR3y3QM4onQvz8xOpyWdB9xS2h+LiK06jDmBavPlVruTL5b3GcApEbFGWy4bA60C173A19oDZOadwLdKd/+IeEP7mIg4BNi/dM8q90iSJEmSpJpRPVXq28ArS/tHwBkRsUOP8Y+2z04ByMybIuKLwDyqpUxXRsQXqAooWwHHAi8ow0/IzF93iX8C8CbgT4HjI2Jr4GyqI7L3Az5O9c9qGfDBTgEy848R8X7gAqrZMFdGxOeAa6iKNe8CXl+GXwGc1SXOjyLibOBQ4LXAJRFxItWysB2BTwBbluHzMvMPXb7TJ4BXAZsA346IXYALy7UDgA+X9t3ARIpmkiRJkiSt9uKJK3dWfxEx3i/9X5k5p0uspwCnA+/scf8ZwFGZuaLbgFKsuQjYpsuQ+4G3ZOaFXa634rwLOBlYs8uQa4C5mXlPjxjrAN8FXt1lyArgs5l5XJ9cXkS1hKvbnjt3Agdm5tW94kxEWea1GGDx4sXMmjXIAVerv37HoUqSNIo8DlySNAxLlixh9uzHD2CenZlLeo0f1EgulRqmzFyRmUcAc6n2vLmdaoPf20v/1Zl5ZK+iTYlzM9XsnGOBhVTLkB4CbgT+HtipX9GmxDkd2JmqmPQb4GHgd1SzbN4DvKRX0abEWJaZc4G3UJ249dvynRYD/wzs2a9oU+JcTTVL53PAr4Cl5fXL8tkOk1G0kSRJkiRpdTGSM260+nPGTWfOuJEkaWKclSNJ6scZN5IkSZIkSSPGwo0kSZIkSVJDWbiRJEmSJElqKAs3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUDOmOwFJkiSp6ebMW9Dz+qL5c6coE0nSqHHGjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQM6Y7AUmSJGlVN2fegr5jFs2fOwWZSJJWN864kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhpqxnQnIEmSJI2COfMW9Ly+aP7cKcpEkrQqccaNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmh3JxYWk302/BQkiRJkrTqccaNJEmSJElSQ1m4kSRJkiRJaiiXSkmSJEkNMMiy50Xz505BJpKkJnHGjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQM6Y7AUmSJEmDmTNvQc/ri+bPnaJMJElTxRk3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDWUhRtJkiRJkqSGsnAjSZIkSZLUUJ4qJUmSJK0m+p06BZ48JUmrGmfcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaiiPA5ckSZJGSL8jwz0uXJKaxRk3kiRJkiRJDWXhRpIkSZIkqaEs3EiSJEmSJDXUyBZuIuIZEXFARHwmIn4QEfdERJbXmROI96qIODcilkTEI+X93Ih41ThizIyIj0TENRHx+4hYGhHXR8QXI2LLccTZPiK+GhE3R8SyiLg7In4cEe+OiIH3NYqIQyPi4oi4IyIejohFEXFWROw+jhgbRcSnI+IXEXFfRNxf2p+OiI0GjSNJkiRJ0iga5c2J7xpGkIgI4KvAUW2XtgAOAg6KiNOAozMze8TZClgAPLft0rbldWREHJaZF/XJ5wjgFGCt2sdrA3uV1+ERcUBm/q5HjLWBc4AD2i49q7wOi4jjMvOzfXLZFTgf2Lzt0k7ldWREvC4zF/aKI0mSJEnSqBrlwk3dYuB64JUTuPdzjBVtfg4cD9wCbAV8FHhBuX438MlOASLiacCFjBVtTgfOBpYB+wEfA9YHzomIPTLzP7rE2R84jWom1V3A3wBXA38CvAs4GNgdODci9svMFV2+0xmMFW0uBb4M3A7sCHy8fLfPRMQdmfm1LrlsAVwAbAosB/6ufEdK7L8GnglcGBE7Z+ZtXXKRJEnSFOp36hR48pQkTaVRLtx8BrgWuDYz74qIOcCt4wkQEVtTFWcAFgJ7Z+ay0r82Ir4PXA7sAhwbEV/PzFs6hDqGalYNwEcz84Tatasi4lLgx8BM4ETgpR1ymQGcTFW0uR94SdvP+teIOAV4L7A38BfANzvE2Qc4rHQvAA7KzMfavtPPgC2B4yPiu5l5b4fv9DdURRuAwzLznNq1n0TEQuA7ZcxngXd2iCFJkiRJ0kgb2T1uMvNTmXlhZj6ZJVMfYqz49b5a0ab1Mx4C3le6M4APtgeIiKcCHyjd64Evdcj1KqpZMAD7RcTOHXI5CNi6tD/fpUD0EeAPtXYnrULUY8B7a0WbVi73AMeW7obAEe0BImJTqsIQwMVtRZtWnHOAi0v3beUeSZIkSZJUM7KFmyer7G3zutK9ITN/2mlc+fzG0j2w3Fe3L7BBaX+jx/KlM2vtgztcP7DL2HouD1HNcgHYISK2qV8vS7ZeVrqXZOaSLrmcSzWrp1surwXWKO2vd4lRz3ONco8kSZIkSaqxcDNxz6bagBiq5VC9tK7PAua0Xdurw7hOFgIPlvaeHa634tyYmXcOkEunOLsxtqlx11wy81GgVajarcwa6pRLzzh9cpEkSZIkaeSN8h43T9Z2tfYNfcbWr2/HynvpDBQnM5dHxC1UpzHV72nNlJk1wVzo0h8kziupnqFtgOs6xLmvVxEpM++IiPuB9TrkIkmSpIbqt4GxmxdL0vBYuJm42bV2tyVFLYu73FfvP9hlk9/2ODsBm0TEWpn5SPl8FtBagjWMXCYS57q2/iAxWnG275BLTxExq8+QzcYTT5IkSZKkJrJwM3FPr7WX9hn7YK39tC5x+sXoFKdVuBl2LsOKM57v1B6jn8X9h0iSJEmStGpzj5uJW7vWfrTP2Edq7XW6xOkXo1ecYecyrDjj+U7tMSRJkiRJGnnOuJm4h2vtNfuMXavWXtZ2rRWnX4xecYadyzDizBwgRj1Oe4x++i2t2gy4dpwxJUmSJElqFAs3E/dArd1vmc+6tXb78qFWnEGWCnWLM+xchhFn5gAx6nEGWVb1uB5HlQPwxFPXJUmSJEla9Vi4mbh64aDfRrn12SHte7MsAV4ErBsRG/TZoLgV5+7axsTDzqUeZ+GTiLPpALnU47hnjSRJ0mrCU6ckaXjc42bi6qcobdtnbP369ROJExEzgK06xcjMpYwVPiY9l7bry4Gbu8RZPyK6nu4UEZtTHQXeKRdJkiRJkkaehZuJuxW4vbT36TN27/J+G7Co7doVtXavOLswtqzoyg7XW3Ge26tY0vYz2uNcy9iGwl1ziYg1gd1b92Rm+ybEg36nXrlIkiRJkjTyLNxMUGYmcH7pbhsRu3caVz5vzU45v9xXdxlwX2m/PbpvznJ4rX1eh+vf6zK2nstM4I2le11m3lS/npkPAP9Wui+PiG5LnQ5mbKZMp1y+D6wo7Xd0iVHPc0W5R5IkSZIk1Vi4eXJOpFoqBHBSRKx0pHXpn1S6y8v4lZTZKl8p3e2AY9rHRMQewBGle3lmdjot6TzgltL+WERs1WHMCcCGtXYnXyzvM4BTImKNtlw2Br5QuvcCX2sPkJl3At8q3f0j4g3tYyLiEGD/0j2r3CNJkiRJkmriiRNARkNE7AlsXftoY8aKGVfSVpDIzDO7xPk8MK90f05V1LiFaj+aY4EXlGufz8yPd4nxdKqNgP+0fHQacDbVEdn7AR+nOqFpGfDizPz3LnFeDVxAVZC7C/gccA1VseZdwOvL0CuAfTPzsS5xvg0cWrqXUhWcbgd2BD7B2F47R2fmqV1izAZ+BmxCVbT6EnBhuXwA8GGq4tDdwAv7nRI1XmW20GKAxYsXM2vWIPskN1u/Tf4kSZJWJ25gLGlVs2TJEmbPfvwcn9nD+j13lAs3ZwJvH3R8ZnZcwhQRTwFOB97Z4/YzgKMyc0W3ARGxNXARsE2XIfcDb8nMC7tcb8V5F3AysGaXIdcAczPznh4x1gG+C7y6y5AVwGcz87g+ubyIaglXtz137gQOzMyre8WZCAs3kiRJqzYLN5JWNZNVuHGp1JOUmSsy8whgLtWeN7dTbfB7e+m/OjOP7FW0KXFuppqdcyzV7Jt7gYeAG4G/B3bqV7QpcU4HdqYqJv0GeBj4HdUsm/cAL+lVtCkxlmXmXOAtwCXAb8t3Wgz8M7Bnv6JNiXM11SydzwG/ApaW1y/LZztMRtFGkiRJkqTVxcjOuNHqzRk3kiRJqzZn3Eha1TjjRpIkSZIkacTMmO4EJEmSJKldv9nGzsiRNCqccSNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaig3J5YkSZK0yum3eTG4gbGk1YMzbiRJkiRJkhrKwo0kSZIkSVJDWbiRJEmSJElqKPe4kSRJkrRa6rcPjnvgSFoVOONGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIZyjxtJkiRJI6nfHjjgPjiSpp8zbiRJkiRJkhrKwo0kSZIkSVJDuVRKkiRJkrrwSHFJ080ZN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lHvcSJIkSdIEeaS4pMnmjBtJkiRJkqSGsnAjSZIkSZLUUC6VkiRJkqRJ5JHikp4MZ9xIkiRJkiQ1lIUbSZIkSZKkhnKplCRJkiRNI0+mktSLM24kSZIkSZIaysKNJEmSJElSQ7lUSpIkSZIazpOppNHljBtJkiRJkqSGsnAjSZIkSZLUUC6VkiRJkqRVnEuppNWXM24kSZIkSZIayhk3kiRJkrSa6zcjB5yVIzWVM24kSZIkSZIaysKNJEmSJElSQ7lUSpIkSZLkBsdSQznjRpIkSZIkqaGccSNJkiRJ6ssNjqXp4YwbSZIkSZKkhnLGjSRJkiRpKNwnRxo+CzdSQwwy9VSSJEmSNFos3EiSJEmSpoT75EjjZ+FGkiRJktQYLreSVubmxJIkSZIkSQ3ljBtJkiRJ0irD5VYaNRZuJEmSJEmrFZdbaXXiUilJkiRJkqSGcsaNJEmSJGmkuNxKqxILN5IkSZIktXG5lZrCwo0kSZIkSeNkYUdTxcKNJEmSJElD5nIsDYuFG0mSJEmSpoGzdjQICzeSJEmSJDWQs3YEFm4kSZIkSVplOWtn9WfhRpMuIrYE3g/MBbYEHgFuBr4D/ENmPjSN6UmSJEnSastZO6s+CzeaVBExF/gWsH7t45nAruV1ZES8OjN/Mx35SZIkSdKoc9ZOs1m40aSJiOdRzaqZCSwFPg9cCqwDHAq8C3gusCAids3MpdOVqyRJkiSps0Fm7fRj8WfiLNxoMp1IVbRZDrwyM6+qXftRRPwaOB7YFvhr4DNTnqEkSZIkadJZ/Jm4p0x3Alo9RcSuwL6le0Zb0ablS8D1pf3BiHjqVOQmSZIkSdKqwsKNJsuBtfbXOw3IzBXAN0t3Q8YKPZIkSZIkCQs3mjx7lfcHgZ/1GHd5rb3n5KUjSZIkSdKqxz1uNFm2K+83Z+byHuNu6HBPXxExq8+QLVqNO+64Y9Cw02r5/fdMdwqSJEmS1FhLliyZ7hR6avvdc41hxY3MHFYsCYCIWBtYVroLMvOAPuOXAusCP83MPQb8GT64kiRJkqSm2jUzFw4jkEulNBmeXmsPcsT3g+X9aZOQiyRJkiRJqyyXSmkyrF1rPzrA+EfK+zrj+Bmz+1xfk+qY8d8CdwOPjSP2VNsMuLa0dwXunMZcpE58RtV0PqNqOp9RNZ3PqJpuVXlG1wA2Ke1fDiuohRtNhodr7TUHGL9WeV/Wc1RNZg6yuPE3g8abThFR79454HeTpozPqJrOZ1RN5zOqpvMZVdOtYs/ofw07oEulNBkeqLUHWf60bnkfZFmVJEmSJEkjw8KNhi4zHwZaRyT1PP0pIjZkrHCzeDLzkiRJkiRpVWPhRpPl+vK+dUT0WpK3bYd7JEmSJEkSFm40ea4o7+sCO/cYt0+tfeXkpSNJkiRJ0qrHwo0my/dq7Xd0GhARTwHeVrr3ApdObkqSJEmSJK1aLNxoUmTmNcBPSveIiNijw7APA9uV9pcz849TkpwkSZIkSasIjwPXZPoA1fKndYAfRsTfUs2qWQc4FDiqjLsJ+NK0ZChJkiRJUoNFZk53DlqNRcRrgP8FrNdlyE3A3My8eeqykiRJkiRp1WDhRpMuIp5FNftmLtXx4I8CNwPnACdn5kPTmJ4kSZIkSY1l4UaSJEmSJKmh3JxYkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaaRhGxZUR8MSKuj4gHI+L3EXFNRBwTETOnOz+tfiLihRHx8Yj4QUQsjohHImJpRNwUEWdGxF7jjPeqiDg3IpaUWEtK/1WT9R00uiLi+IjI2mvfAe7xGdWkiYiNI+KjEXFlRNxZnrHbI+LqiDghIvYYIIbPqCZFRKwZEUdExL9GxB21f+ffGBH/FBG7DxjHZ1TjEhHPiIgDIuIz5b8576n9u/vMCcR70s9gRMyMiI+U37V+X/63cH35XWzL8eY01SIzpzsHaSRFxFzgW8D6XYbcCLw6M38zdVlpdRYRlwN7DzD0LODIzHy0R6wAvgoc1SPOacDR6b9oNAQR8TxgITCj9vF+mXlZl/E+o5pUEXEI8I/ARj2GnZ+ZB3a532dUkyYiZgMLgB37DP174MOdnjGfUU1URPR6Hr6RmYcPGGcoz2BEbEX1v4fndhlyH3BYZl40SF7TwRk30jQov4B8h6posxT4BPBi4GXA6WXYc4EFEfG0aUlSq6MtyvvtwJeBNwC7AXsAfw3cVq6/FTizT6zPMfYv0Z8Dby6x3lz6lOufHULeGnER8RSq/2+cAfx2wNt8RjVpIuJtwNlURZvfAp8GXgHsDMwF3g9cAvyxRxifUU2KiJjBykWb/wAOp/r3/SuBzwAPlmsfAo7pEspnVMOwGPjhBO990s9g+V3qQsaKNqdT/c71YqrfwZZS/U52TkTsNME8J50zbqRpEBGXAvsCy4G9M/OqtusfAY4v3U9l5memNkOtjiLiQuCbwL9k5mMdrm8MXAn8aflo78z8SYdxWwPXU/0SvbCMW1a7PhO4HNiF6hnfNjNvGfLX0QiJiA9S/a3wDcB5wMfKpY4zbnxGNZkiYjuqXxjWAn4CvCYz7+syds1Osxd9RjWZIuL1wHdL9ypgr/Z/70fEzuXaU4E/AM/IzOW16z6jmrCI+DRwLXBtZt4VEXOAW8vlgWbcDOsZjIjjgE+V7kcz84S263sAPy4/59LMfOng33TqOONGmmIRsStV0QbgjPaiTfElqv+jAvhgRDx1KnLT6i0zD8jM73Qq2pTr9wAfrn30hi6hPsTYcpX31f8lWuI8BLyvdGcAH5xw0hp5Zbp/62/S3gN0XcJX4zOqyXQSVdHmHuDgbkUbgB5LTn1GNZleUmt/vtO/9zPzZ1SzEAA2BLZtG+IzqgnLzE9l5oWZedeTCPOkn8HyO9QHSvd6qt+x2nO9CjijdPcrRc3GsXAjTb0Da+2vdxqQmSuoZkZA9S/TfSc3Jelxl9XaW7VfLGuNX1e6N2TmTzsFKZ/fWLoHlvukifgH4GlUf0N3Wb/BPqOaTBGxLdUUe4CTS8F7vDF8RjXZ1qy1e+2VWJ+dsFar4TOq6TbEZ3BfYIPS/kb5HauTM2vtg8eb71SwcCNNvdapPQ8CP+sx7vJae8/JS0daSf0/9jr9y+3ZjO2Vc3mH63Wt67OAOU8uLY2iiHgjcADwe+AjA97mM6rJdEitfU6rEREbRsQ2EdFro+IWn1FNtptq7ef0GNf6C5oEfl373GdU021Yz+BeHcZ1spCxfZ8a+XuXhRtp6m1X3m+uryXu4IYO90iTbZ9a+4YO17frc50u132GNS4RsQHVJtoAx2bm3QPe6jOqydQ6Pvk+4PqIeEtE/IKquHgTcE9E/CYiPtXjcAGfUU22bwP3l/axEbFG+4CIeAHVRtoAZ2fm/bXLPqOabsN6BgeKU34na81Aa+RzbOFGmkIRsTawceku6TU2M//AWOV39mTmJcHjJ/fMq330nQ7D6s9iz2eY6hSBTvdJgzge2Az4v4ytPR+Ez6gm05+V90VUe938L6D9FJJnA8cBV0XEMzvE8BnVpCqF7sOBZVT73VwbEW+LiN0j4uUR8Smq2QdrAv9OdbJknc+optuwnsFW/8HMvHfAOJtExFo9R04DCzfS1Hp6rb10gPGtwo1HgmsqfIjqiEWA8zJzYYcx43mGH6y1fYY1sIjYEziS6pSIo3N8R2D6jGoy/Ul53xb4S+Be4GjgGcDawK7AD8qYHaiOl23/722fUU26zDyP6rSdM4DnA9+gOkXqEqrC4kNUBZs9M/POttt9RjXdhvUMtuKM5/euTnGmnYUbaWqtXWsPcjLKI+V9nUnIRXpcROwDzC/d31Kd3tPJeJ7hR2ptn2ENJCLWBE4DAvj7zPzlOEP4jGoyrVve1wIeA/48M0/NzLsz85FS8D6AseLNi3niRpc+o5p05TSdw4DXUP3/abtNgTfT+QAMn1FNt2E9g6044/m9q1OcaWfhRppaD9faa3YdNaY1TW9Zz1HSkxAR2wPnUR2l+Ajwxh7HN47nGa5PM/UZ1qA+TrW+/L+BT0/gfp9RTab683VOp5NOyqkl9c2039wjhs+ohi4i1gX+D/AJYCOqpafbUT1P6wOvBK6gmiF2QUR8oC2Ez6im27CewVac8fze1SnOtLNwI02tB2rtQabgtf5mb5DpfdK4RcSzgR9SHTv/GPDmzOy16/54nuF1a22fYfVVjlr+WOm+LzMf7DW+C59RTab68/WDboMy8z+B20p31x4xfEY1GT4N7F3aR2TmsZl5Q2Y+mpn3Z+YlwH7ApVSzcf4uIup7NfmMaroN6xlsxRnP712d4ky7GdOdgDRKMvPhiLiHaoPiWb3GRsSGjP0fyOJeY6WJKJtm/h/gmVRHgb6zrInvpb5BXM9nmJU3iPMZ1iA+RPW3Yr8BZkbEoR3G7FBrvzQiNivtC0qhx2dUk2kx1abZMNiGmVtQ7X9T5zOqSRMRAbyjdG/KzG90GpeZyyPif1DNvHlKuedD5bLPqKbbsJ7BJcCLgHUjYoM+GxS34tydmY/0GDctLNxIU+96YC9g64iY0eNI8G3b7pGGJiI2ptqg8Dnlo/dl5jcHuPW6WnvbrqOeeN1nWINoTVN+DtVxtv38j1r72VQbC/qMajL9J2MzaJ5wxHKb1vX2f8/7jGoybcrYJto/7zP2Z7V2/VnzGdV0G9YzeB3w+tq4JyxvBYiIGcBWXWI0gkulpKl3RXlfF9i5x7h9au0rJy8djZqIWB+4mLFjbedl5ikD3n4rcHtp79NrIGPTtG+jOjpXmgo+o5pMP661t+o6qtIqjN/W9rnPqCZTvVDY7y/pn9rlPp9RTbdhPYNX1Nq94uzC2EqHRv7eZeFGmnrfq7Xf0WlAOTr0baV7L9UaZOlJi4iZwALgheWjv8nMLwx6fzmW+fzS3TYidu/yc3Zn7G9Azh/ncc4aUZl5eGZGrxcrb1i8X+3aohLDZ1ST6fvAH0u7/bSox5WT+jYq3Z/Ur/mMapL9Hri/tPcoMwm6qf8ie2ur4TOq6TbEZ/Ay4L7SfntZStjJ4bV2v20DpoWFG2mKZeY1jP1H3BERsUeHYR+m2v0f4MuZ+ccOY6RxKccsnwe8pHz05cz85ARCncjY38ydFBErHZlY+ieV7vIyXppKJ+IzqkmQmb8Dvla6r+i0D1NEPJ2Vn6lTO4Q6EZ9RTYJyqtmC0n0m1clST1D2Uqz/xc2FbUNOxGdU0+tEnuQzmJmPAl8p3e2AY9rHlN/FjijdyzPz2ieb+GQIC6PS1IuIF1BNw1uHatfyv6WaVbMOcChwVBl6E7BLZj7QKY40HhHxL4z9DfGPgA9SbUrczaOZeVOXWJ8H5pXuz6n+4+8WqqUDxwIvKNc+n5kff3KZS2Mi4jjgU6W7X2Ze1mWcz6gmRURsAiwEtqT6ZeGrwLlUsxx2pHq+Wn8D/I+Z+d4ucXxGNSnKCX0/A2aWjy4AvkG18fvawO5U/w2wZbn+b5n58g5xfEY1IRGxJ7B17aONgRNK+0rGCuAAZOaZXeI86WewFNMXAn9aPjoNOJvqyO/9gI9TnTq1DHhxZv77AF9xylm4kaZJRLwG+F/Ael2G3ATMzcybpy4rrc4iYrz/h/9fmTmnS6ynAKcD7+xx/xnAUeVv/6ShGEfhxmdUkyYitqNaNrV1j2H/BBzdbdasz6gmU0S8nGqT9437DP0R8IbM/EOHGD6jmpCIOBN4+6Djy1LoTnGG8gxGxNbARcA2XYbcD7wlM9tnnjWGhRtpGkXEs4APAHOpjrp7FLgZOAc4OTMfmsb0tJoZZuGmFvPVVDPEdqX6j8N7gGuBUzPzBxPJU+pl0MJNbbzPqCZFRKwLvAd4A9UvA08Dfkv1t8mnZuZA+9P5jGqyRMRGVEtA/hzYHtiAapbYnVTP2D8D3++3N43PqMZrWIWbWrwn/QyW/8/+S+AQqqL7mlTHh19EtX3Afw2a73SwcCNJkiRJktRQbk4sSZIkSZLUUBZuJEmSJEmSGsrCjSRJkiRJUkNZuJEkSZIkSWooCzeSJEmSJEkNZeFGkiRJkiSpoSzcSJIkSZIkNZSFG0mSJEmSpIaycCNJkiRJktRQFm4kSZIkSZIaysKNJEmSJElSQ1m4kSRJkiRJaigLN5IkSZIkSQ1l4UaSJEmSJKmhLNxIkiRJkiQ1lIUbSZIkSZKkhrJwI0mSJEmS1FAWbiRJkiRJkhrKwo0kSZIkSVJDWbiRJElazUTEmRGRHV5zpju36RAR+3b58zhuunOTJKkfCzeSJEl6XEQ8JSIOjIh/jIhfRMRdEfFoRNwfEb+JiPMjYl5E/Okk/fxTa4WV/cZ578tq9548GflJkjTVZkx3ApIkSZo0twP71/q39RocEa8GvgRs2+HyU4GnA88GXgt8PiIuBz6emf93OOkC8E3gqNJ+K3DpOO79i1r7rFr7WmDHWv+XE0tNkqSpZ+FGkiRp9fXHzPzVIAMj4ljg80CUj64ELgB+DvwOWBvYFHgJMBd4LrAP8D+BVw0r4cy8MiJuAbYC3hARf5mZywbIfx3g9aV7Y2ZeXYv5IPCr2thhpStJ0qSzcCNJkjTiIuJtwPzSvQd4S2b+sMvwcyPiGOA1VIWeyfBN4NNUM3xeB5w9wD0HlvGw8mwbSZJWae5xI0mSNMIiYgvgq6X7ILB3j6INAFn5PrAzcMYkpHUWkKX91gHvaY1L4H8NPSNJkqaJhRtJkqQJiohNa5vhvjsi1ouIT0TElRFxd9nU97/LRr+bTXe+Xfw1sE5pfzIzrx/0xsx8ODPP6TUmInaLiNMj4qaIWBoRD0bEDRFxSkRs0yXurcAVpfvKiHhGn5+xKfCK0r08M/9r0O8gSVLTWbiRJEmauOfX2gn8Gvgc8GJgY6oNfWcDRwNXR8QzpzrBXqLa7OVtpbuUIc6eiYgZEfEPwNXAkcA2wLrATKr9cd4L/GdEvKtLiG+W9xnAm/v8uDcztgXAN3sNlCRpVWPhRpIkaeKeX2t/EXgGcA7Vviw7A28C/qNc3xL4+6lMbgDbUxWYAH6SmQ8MMfYZwHtK+wdUJz7tBuwKvAv4T6rC1mkR8ZoO938HaG1K3G+5VOv6Q8B3n0TOkiQ1jpsTS5IkTdzzau11gDe2LR36fxHxA6qZOJsCB0fEBpl57xTm2MtOtfb/G1bQiHg9YzN53pWZX2sbsjAi/hewAHgp8JWI+EFmLm8NyMz7I+J84FBg54jYrtMyroj4M+CFpfu9IRefJEmads64kSRJmrjn19of7bTfSykknFy6M4AXTEFeg9q41r6718CI2D4idujyWrdt+MfK+3kdijZAtT8O8FelOwfYt8Ow+rKnbrNu6p+7TEqStNqxcCNJkjQBEbE28Kelex3wlR7Df1Vrb9x11NR7eq29tM/YXwC/7PLatTWonFK1c+l+p1fAMoPmntLdo8OQHwJ3lPZbyp48jyv9t5TuHcD/6fMdJEla5Vi4kSRJmpgdgTVK+6uZ+ViPsX+otR+qnUQ10GvyvgL1ZUXts2Ymapda+9sDfLdWIesJp26VP9N/Lt0tgX3ahuxLtfkzwD/3+WcgSdIqyT1uJEmSJqa+v833+4ytFyXuAj7dYcyngPuAE59cWuPyu1p7k14DM3Ol/26MiOOocm7X8+juHmZ2+fwbwIdL+63AZbVrLpOSJK32LNxIkiRNzPPL+z2Z+V99xu5e3h8Brs/Mhe0DIuJTwL2ZedzQMuzvF7X2C7uOGp81au23MHaqVj9/6PRhZv4yIn5BVSh7Q0T8VWYui4h1gNeXYb/IzEF/jiRJqxQLN5IkSRPTmnFzV69BZR+WA0r3x5n54KRmNT7/STXrZiNgr4hYdwj51WfxZGb+quvIwX0D+DtgPeC1wP+mOnJ9vXLd2TaSpNWWe9xIkiSNUynGtI7S7rbEp+UVwNal/Y1JS2oCMjMZK3o8HTh8CGF/Xmu/cgjxoNrnpnVU+Fvb3uv74EiStNqxcCNJkjR+z2ZstsfsiHhap0ERsSZwQuneCHx7CnIbr78DlpX230bE1r0G95OZN1OdsgVwaERs+WTilZh3UZ0wBbB/ROzAWFHoh5l555P9GZIkNZWFG0mSpPGrb0w8g7HNcx8XEWsBZ1LNzHkMeE9mrpiS7MYhM5cAf1m66wE/iYh9B7h1wx7XPlfe1wbOjYiuGx9HxFoR8d5yvHovrdlKM4CzGVvy7zIpSdJqzT1uJEmSxu/55f33wL3AcRGxOfA9qk12nw+8D9i+jPtkZl46pRmOQ2Z+PSK2AD5DdQLWpRHxY6rTsv6Dat+aoDox6nnAQcButRDL2uJ9OyL2B94O7AxcFxGnApcDd1MdPb4VsBdwMPAn9C/AfJ/qz3oDxv5c7wfOH/cXliRpFWLhRpIkafxaM27+A/hb4ALg3eVVtxz4eGaeQMNl5ufK6U1fArYB9i6vXq4Ejs3MqztcO4Jq4+YPAxsDnyivTh6kmpXUK7+HI+Ic4F21j8/JzGXd7pEkaXXgUilJkqTxe355/4/MvIRq5sh5VIWKR4BbgVOBHVeFok1LZl4AbEc1C+Y04JdUM2SWAw8A/wVcBBwHbJ+Ze2bmlV1iPZaZxwJ/RlUM+jnVbKTHSqz/BL5FNStn8wELMO2bO7tMSpK02nPGjSRJ0jhExPrAs0r3lwCZeS1VsWOVl5mPURWhzhtSvJuAY4YU60qqJVuSJI0MCzeSJEnj8/xa+5fTlcSAnlpOYGq5MTP/OG3ZTJOIWJfqJDBJklY5Fm4kSZLGp7W/TVIt92myZ7JycenZwKLpSWVa7Qo0dnNoSZJ6cY8bSZKk8Xl+ef9NZi6dzkT0/9u5YxuEgiCGgl6JbsjohE5phJqWhALIvjnNZJdt/GQdAJxvdvfqGwAA/sbMvJM8krx293n1PQDA2SxuAAB+NDO3JPfvs/1/GwDgABY3AAAAAKUsbgAAAABKCTcAAAAApYQbAAAAgFLCDQAAAEAp4QYAAACglHADAAAAUEq4AQAAACgl3AAAAACUEm4AAAAASgk3AAAAAKWEGwAAAIBSwg0AAABAKeEGAAAAoJRwAwAAAFBKuAEAAAAoJdwAAAAAlBJuAAAAAEoJNwAAAAClhBsAAACAUsINAAAAQKkPljELl8FzvyMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.hist([filtered_missing_ET_5], bins=100, range=(0, 100))\n", "plt.xlabel(r'$p_\\mathrm{T}$ [GeV]')\n", "plt.ylabel('Jets')\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.12" } }, "nbformat": 4, "nbformat_minor": 5 }