{ "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": [ "[root://eospublic.ce...: 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": [ " [root://eospublic.ce... Downloaded: 0%| | 0/9000000000.0 [00:00]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dataset_name = [\"root://eospublic.cern.ch//eos/root-eos/benchmark/Run2012B_SingleMu.root\"]\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": "iVBORw0KGgoAAAANSUhEUgAABG4AAALlCAYAAACVaqBLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAB7CAAAewgFu0HU+AABeK0lEQVR4nO3de7xtVV03/s9Xj4DgNdE0gY6CCWmmiYR5QbppHjNQMy9PiaJmlqZPKketR7yUJFqWWl6yvJT6iA/kBUstFZUswZ+VJYiox87xCt64g8D4/THndi+Wa629z9mXNc/Z7/frtV9rzDXH/K6xjlNgf86YY1RrLQAAAAAMz/XmPQAAAAAAJhPcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIHaNO8BwFqoqr2T/ER/eEGSa+Y4HAAAAPZ8109yy7796dbalatRVHDDnuonkpw170EAAACwId0jydmrUcijUgAAAAADZcYNe6oLFhqf+MQncpvb3GaeYwEAAGAP99WvfjVHHHHEwuEFs/ruDMENe6rvr2lzm9vcJgcccMA8xwIAAMDGsmrrrHpUCgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAA7Vp3gMA2F1t3nr6kn22nbRlHUYCAADsqcy4AQAAABgoM24A1tBSs3LMyAEAAGYR3ABMsZxHoQAAANaSR6UAAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAANlO3BgwxrCdt/LGcO2k7asw0gAAIAhMuMGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoDbNewAAzLZ56+kzz287acs6jQQAAFhvZtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAbKduDAHmmpLbQBAAB2B2bcAAAAAAyU4AYAAABgoAQ3AAAAAAO1YYObqrpVVT2oql5QVf9QVRdWVet/3rAL9R5QVadW1Y6qurJ/PbWqHrATNfatqmdW1Seq6ltVdUlVnVNVL62qg3aizp2q6tVVdX5VXV5VF1TVR6rqN6tq2esaVdUjqup9VfXVqrqiqrZV1Zur6sidqHGLqnp+Vf1HVX23qi7q28+vqlsstw4AAABsRBt5ceKvr0aRqqokr07yxLFTt01ybJJjq+q1SZ7UWmsz6hyc5PQkdxw7dWj/8/iqelRr7b1LjOf4JK9KsvfI2/skuU//c1xVPai19s0ZNfZJckqSB42d+tH+51FVdWJr7YVLjOUeSd6Z5DZjp+7S/zy+qn6ltXb2rDoAAACwUW3YGTdjtid5/y5e+6IshjafSvLIJEf0r5/q339ikqkhR1XdKMl7shjavC7JzyX5mSTPTXJJkpsmOaWq7jKjzv2TvDZdaPP1JE9N8tNJfinJqX23I5OcWlWz/rd/fRZDmw8lOab/Tscn+Xy6++YFVfX4GWO5bZJ3pwttrk7ykiT37X9e0r/3I0ne0/cFAAAAxmzkGTcvSHJWkrNaa1+vqs1JvrgzBarqkCTP6g/PTnLf1trl/fFZVfWuJGckOTzJCVX1N621z08o9Yx0s2qS5FmttZNHzn28qj6U5CNJ9k3y8iQ/O2Esm5K8Ml2oclGSe4191j9W1auSPDldePK/krxpQp2jkjyqP3x3kmNba9eMfadPJjkoyUuq6h2tte9M+E5/mOSH+/ajWmunjJz7aFWdneTtfZ8XJnnchBoAAACwoW3YGTettee11t7TWlvJI1NPz2L49ZSR0GbhMy5L8pT+cFOSp40XqKobJPnd/vCcJC+bMNaPp5sFkyRHV9XdJ4zl2CSH9O0XTwmInpnk2yPtSRaCqGuSPHkktFkYy4VJTugPb55uFs51VNUPpwuGkuR9Y6HNQp1TkryvP/yN/hoAAABgxEaecbMi/do2v9Ifntta+9dJ/Vpr/1pVn033GNQxVfXUsbVu7pfkZn37ja21a6d85BuS/Gbffki6WS+jjhnrO2ksl1XV2/s6d66qO7TWPjfynW6U7hGtJPlAa23HlLGcmm5Wz036sYyHTQ9Ocv2+/TdTaiyM8/593wene0QM2Embt56+ZJ9tJ21Zh5EAAACrbcPOuFkFt0u3AHHSPQ41y8L5A5JsHjt3nwn9Jjk7yaV9+94Tzi/U+Wxr7WvLGMukOkdkcVHjqWNprV2VZCGoOqKfNTRpLDPrLDEWAAAA2PDMuNl1h420z12i7+j5w3LdtXSWVae1dnVVfT7dbkyj1yzMlDlgF8eSKcfLqfOL6e6hOyT5zIQ6350VIrXWvlpVCzN3xscyU1UdsESXW+9MPQAAABgiwc2uO3CkPe2RogXbp1w3enzplEV+x+vcJcktq2rv1tqV/fsHJKlVHMuu1PnM2PFyaizUudOEsSznOgAAANijeVRq1914pH3JEn0vHWnfaEqdpWrMqrPaY1mtOjvzncZrAAAAwIZnxs2u22ekfdUSfa8cad9wSp2lasyqs9pjWa06O/OdxmssZakZOrdOt907AAAA7LYEN7vuipH2Xkv03XukffnYuYU6S9WYVWe1x7IadfZdRo3ROuM1Zpqx41WSpNv0CwAAAHZvHpXadRePtJd6zGe/kfb440MLdZbzqNC0Oqs9ltWqszPfaTmPVQEAAMCGIrjZdaMzPpba4Wj0sZ7xRXUX6uxXVTdbZp0LRhYmXouxrFadpWqM1rHYMAAAAIwR3Oy60V2UDl2i7+j5c3alTlVtSnLwpBqttUuyGHys+VjGzl+d5PwpdW5aVVO35a6q26TbCnzSWAAAAGDDE9zsui8m+UrfPmqJvvftX7+cZNvYuY+NtGfVOTyLjxWdOeH8Qp07zgpLxj5jvM5ZWVxQeOpYqmqvJEcuXNNaG1+EeLnfadZYAAAAYMMT3Oyi1lpL8s7+8NCqOnJSv/79hdkp7+yvG/XhJN/t24+p6avqHjfSPm3C+b+f0nd0LPsmeXh/+JnW2nmj51trFyf55/7w56tq2qNOD8niTJlJY3lXkmv79mOn1Bgd57X9NQAAAMAIu0qtzMuTPCHdn+Mrquq+rbXv745UVTdM8or+8Oq+/3W01q6qqj9P8gdJDkvyjCQnj/apqnsmOb4/PKO1Nmmb69OSfD7d41TPrqpTWmufH+tzcpKbj7QneWmSX+q/06uq6iGttWtGxrJ/kj/uD7+T5K8mfKevVdXfJfn1JPevqoe11t4x9p1+Ncn9+8M3t9a+NmU8MNHmrafPewgAAABrbsMGN1V17ySHjLy1/0j7kKo6brR/a+0N4zVaa+dV1UuTbE33KNOZVfXHWQxQTkhyt777ya21z00ZzslJfi3JjyV5SVUdkuRt6bbIPjrJc9L9b3V5kqdNKtBa+15VPTXJu9PNhjmzql6U5BPpwponJHlo3/1jSd48pc4Hq+ptSR6R5MFJPlBVL0/3WNhPJHlukoP67ltba9+e8p2em+QBSW6Z5K1VdXiS9/TnHpTk9/r2BUl+f0oNAAAA2NDqB5/c2Riq6g1JHrPc/q21iY8wVdX1krwuyeNmXP76JE9srV07rUMf1rw3yR2mdLkoyaNba++Zcn6hzhOSvDLJXlO6fCLJltbahTNq3DDJO5I8cEqXa5O8sLV24hJj+el0j3BNW3Pna0mOaa3926w6u6J/zGt7kmzfvj0HHLCcDa7YnZhxs3O2nbRl3kMAAIA92o4dO3Lggd/fgPnA1tqOWf2Xyxo3K9Rau7a1dnySLenWvPlKugV+v9IfP7C19vhZoU1f5/x0s3NOSHJ2useQLkvy2SR/muQuS4U2fZ3XJbl7ujDpC0muSPLNdLNsfivJvWaFNn2Ny1trW5I8OskHknyj/07bk7wlyb2XCm36Ov+WbpbOi5L8V5JL+p9P9+/deS1CGwAAANhTbNgZN+zZzLjZ85lxs7rMyAEAgJUx4wYAAABggxHcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQG2a9wAAmL/NW09fss+2k7asw0gAAIBRZtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGKhN8x4AwLjNW0+f9xAAAAAGwYwbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoDbNewAA7B42bz195vltJ21Zp5EAAMDGYcYNAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAzUpnkPANh4Nm89fd5DAAAA2C2YcQMAAAAwUIIbAAAAgIHyqBQAq2I5j8BtO2nLOowEAAD2HGbcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbhZJVW1V1UdX1X/WFVfraorq+qSqvpsVf11VR25zDoPqKpTq2pHX2NHf/yAnRjLvlX1zKr6RFV9qx/HOVX10qo6aCfq3KmqXl1V51fV5VV1QVV9pKp+s6o27USdR1TV+/o/lyuqaltVvXm5fyYAAACwUS37l2+mq6oDk5ye5CfGTu2V5Mf6n8dW1Z8m+b3WWptQo5K8OskTx07dNsmxSY6tqtcmedKk60fqHNyP5Y5jpw7tfx5fVY9qrb13ie90fJJXJdl75O19ktyn/zmuqh7UWvvmjBr7JDklyYPGTv1o//OoqjqxtfbCWWMBAACAjcqMmxXqZ56Mhjb/meS4JPdM8otJXpDk0v7c05M8Y0qpF2UxtPlUkkcmOaJ//VT//hOTTA05qupGSd6TxdDmdUl+LsnPJHlukkuS3DTJKVV1lxl17p/ktelCm68neWqSn07yS0lO7bsdmeTUqpp1D70+i6HNh5Ic03+n45N8Pt3994KqevyMGgAAALBh1YzJGyxDVT00yTv6w48nuU9r7ZqxPnfvz90gybeT3Kq1dvXI+UOSnJNuBtTZSe7bWrt85Py+Sc5IcniSq5Mc2lr7/ISxnJjkef3hs1prJ4+dv2eSj/Sf86HW2s9OqLGpH8shSS5K8lPjn1VVr0ry5P7wMa21N02oc1SSD/eH705y7OifS1Xtn+STSQ7q/0xu31r7znidXVVVByTZniTbt2/PAQccsFqlWQWbt54+7yEwUNtO2jLvIQAAwC7ZsWNHDjzwwIXDA1trO1ajrhk3K3evkfaLx0ObJGmtfTLdTJgkuXm6R5ZGPT2Lj609ZTS06a+/LMlT+sNNSZ42/hlVdYMkv9sfnpPkZRPG8fF0s2CS5Og+UBp3bLrQZuH7/EBAlOSZ6cKWhfYkz+pfr0ny5PE/l9bahUlO6A9vnm4WDgAAADBCcLNye420vzCj32gA8v11Y/q1bX6lPzy3tfavky7u3/9sf3hMf92o+yW5Wd9+Y2vt2injeMNI+yETzh8zpe/oWC5L8vb+8M5VdYfR8/0jWz/XH35gRsp4arpZPdPGAgAAABua4Gblzhtp335Gv4P715bkcyPv3y7dAsRJ9zjULAvnD0iyeezcfSb0m+TsLK65c+8J5xfqfLa19rVljGVSnSOyGE5NHUtr7aokC0HVEf2sIQAAAKAnuFm5t2Zx1sgJVXX98Q5VdbckCws3vK21dtHI6cNG2ucu8Vmj5w8bO7esOv3aOguzf65To58ps7AYzJqPZez8piR3mNURAAAANhrbga9Qa+2Cqjouyd+lW+/mrKp6ebqZODfq3/u9dI9U/XuS/z1W4sCR9lILF22fct3o8aXLWOR3e5K7JLllVe3dWruyf/+AJAuPYK3GWHalzmeW6J/k+4sPz3Lr5dQBAACAIRPcrILW2mlVdXi6UOZxSd441uXr6XZ7em1r7dKxczceaV+yxEeNXnujKXWWqjGpzkJws9pjWWmdWbYv3QUAAAB2bx6VWgX92iyPSvLLWZyxMuqHkzwy3QLC4/YZaV+1xEddOdK+4ZQ6S9WYVWe1x7LSOgAAALChCW5WqKr2S/JPSZ6b5BZJXpJujZe9k9w0yS8m+ViSeyR5d1X97liJK0bae2W2vUfal4+dW6izVI1ZdVZ7LCutM8uBS/zcYydqAQAAwCB5VGrlnp/kvn37+Nba6GNSVyX5QFV9KMn7kxyd5E+q6kOttf/s+1w80n+pR4X2G2mPP4K0UGc5jxtNq7PaY1lpnalmbDGeJPnB3dIBAABg92PGzQpUlw48tj88byy0+b5+J6c/6A+vN3JNct3Fe5dacHd00d/xNV4W6uxXVTdbZp0LRhYmXouxrLQOAAAAbGiCm5X54SQ/1Lc/tUTfT460Dx1pf2bK+5OMnj9n7Nyy6lTVpiQHT6rRWrski+HJmo9l7PzVSc5foi8AAABsKIKblbl6pL3UY2c3mHLdF5N8pW8ftUSNhUeyvpxk29i5j420Z9U5PIuPJ5054fxCnTtW1awttUc/Y7zOWVlclHjqWKpqryRHLlzTWlvOwsoAAACwYQhuVuZbSS7q2/fsZ7NMMxpgfHGh0VprSd7ZHx5aVUdmgv79hdkp7+yvG/XhJN/t24+p6Yu8HDfSPm3C+b+f0nd0LPsmeXh/+JnW2nmj51trFyf55/7w56tq2uNSD0lykxljAQAAgA1NcLMCrbVrk5zeH/5Iup2lfkBV3TzJH4+89Z6xLi/P4iycV1TVdbbF7o9f0R9e3fcfH8tVSf68PzwsyTMmjOOeSY7vD89orZ01YbinJfl83352VR08oc/JSW4+0p7kpf3rpiSvqqrrj41l/yz+mXwnyV9NqQMAAAAbluBm5V6Q5LK+fWJVvauqHlpVd6uqe1bV05P8e5If7/v8c2vt/aMF+hkrC0HH4UnOrKpfq6rDq+rX0j2KdHh//uTW2uemjOXkJAuzX15SVa+pqqOr6siqena6na02pdt2+2mTCrTWvpfkqUmuTTcb5syq+p2qOqKq7l9V70jy5L77x5K8eUqdDyZ5W3/44HS7az24/06PTfKvSQ7qz29trX17yncCAACADat+8IkbdlZV/XyStybZf4muH0zysEkhRVVdL8nrkjxuxvWvT/LEfqbPtLEckuS9Se4wpctFSR7dWhuf9TNe5wlJXplkryldPpFkS2vtwhk1bpjkHUkeOKXLtUle2Fo7cdZYdkX/eNb2JNm+fXsOOGCpza1YT5u3nr50JzakbSdtmfcQAABgl+zYsSMHHvj9jZMPbK3tmNV/ucy4WQWttX9Kt/7MCenWmrkgyffSzWz5YpK3Jzkmyc9Pm1nSWru2tXZ8ki3p1rz5SroFfr/SHz+wtfb4WaFNX+f8JHfrx3J2useQLkvy2SR/muQuS4U2fZ3XJbl7ujDpC0muSPLNdLNsfivJvWaFNn2Ny1trW5I8OskHknyj/07bk7wlyb3XIrQBAACAPYUZN+yRzLgZNjNumMaMGwAAdldm3AAAAABsMIIbAAAAgIHaNO8BAHsWj0EBAACsHjNuAAAAAAZKcAMAAAAwUB6VAmAwlnrUzq5TAABsNGbcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBAbZr3AABguTZvPX3JPttO2rIOIwEAgPVhxg0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICB2jTvAQDAatq89fSZ57edtGWdRgIAACtnxg0AAADAQAluAAAAAAZKcAMAAAAwUNa4AXbKUuuHAAAAsHrMuAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABmouwU1V7VtV+844/5Sq+mhVnVNV762qB63n+AAAAACGYN2Dm6r65SQXJ/lKVd14wvm/TvLyJD+T5I5J7p/knVX1rPUcJwAAAMC8zWPGzf2TVJK/b61dPHqiqu6d5Lj+8LIkn0pyRd//RVV1p3UcJwAAAMBcbZrDZx6ZpCX50IRzT+xfv5Lknq21HVV1YJKPJTkgyW8meeq6jBKAPdLmracv2WfbSVvWYSQAALC0ecy4uVX/+rkJ5x6QLtR5RWttR5K01rYneUW6WTdHrcsIAQAAAAZgHsHNLfvXS0bfrKofT7J/f/iusWvO7l83r92wAAAAAIZlHsHNNf3rD429f5/+9YLW2rlj577dv+6zZqMCAAAAGJh5BDdf7l/vOvb+lnSPSX10wjU37V8vXKMxAQAAAAzOPIKbj6Zbr+Z3qmr/JKmqe6Rb3yZJ3jfhmsP616+t/fAAAAAAhmEewc1fJLk2ye2SfKGqzk5yRrodrr6d5P9OuOZn083G+fd1GiMAAADA3K17cNNa+/+SPDNdEHOjJD+Vbu2a7yV5Qmvt4tH+VXXTdI9RJckH1nGoAAAAAHO1aR4f2lr706r6pyQPS3LrJF9N8tbW2mcndL9fkrP69j+tzwgBAAAA5m8uwU2StNY+neTTy+j3ziTvXPsRAQAAAAzLuj8qVVUH9T/X34lrrrdw3VqObaWqav+qelZVnVlVX6uqK6vqK1X1b1V1clXdcxk1HlBVp1bVjv76Hf3xA5a6dqTGvlX1zKr6RFV9q6ouqapzquqlO/NnWFV3qqpXV9X5VXV5VV1QVR+pqt+sqmWHflX1iKp6X1V9taquqKptVfXmqjpyuTUAAABgI5rHjJtt6RYnvkuSzyzzmtsl+Vx/3dxmCc1SVb+a5C+T3GLs1G36nyOS3CHJMVOurySvTvLEsVO3TXJskmOr6rVJntRaazPGcXCS05PccezUof3P46vqUa219y7xfY5P8qoke4+8vU+S+/Q/x1XVg1pr35xRY58kpyR50NipH+1/HlVVJ7bWXjhrLAAAALBRzWNXqaTbDnw9r1tTVfUbSd6WLrT5RpLnJ/mFJHdPt7DyU9MtrPy9GWVelMXQ5lNJHpku7Hlkf5z+/NSQo6pulOQ9WQxtXpfk55L8TJLnJrkkyU2TnFJVd5lR5/5JXpsutPl6P/6fTvJLSU7tux2Z5NSqmnUPvT6Loc2H0oVWRyQ5Psnn091/L6iqx8+oAQAAABvWIGevTLAQ2Fw711FMUFWHpQs5rpfko0l+ubX23QldX1FVe02pcUiSZ/WHZye5b2vt8v74rKp6V7ot0w9PckJV/U1r7fMTSj0j3ayaJHlWa+3kkXMfr6oPJflIkn2TvDzdNuvjY9mU5JX997koyb3GPusfq+pVSZ6c5L5J/leSN02oc1SSR/WH705ybGvtmrHv9MkkByV5SVW9o7X2nQnfCQAAADasec242Vm36V8vntlrPl6RbmbKhUkeMiW0SZK01q6acurpWQzRnjIS2ixcd1mSp/SHm5I8bbxAVd0gye/2h+ckedmEz/94ulkwSXJ0Vd19wliOTXJI337xlIDomUm+PdKeZCGIuibJk0dCm4WxXJjkhP7w5ulm4QAAAAAj5hncTF2nZUFV3aCqDk33mE+STNoufG76sf1cf/jKPozY2RqV5Ff6w3Nba/86qV///sL3P6a/btT9ktysb7+xtTZtdtIbRtoPmXD+mCl9R8dyWZK394d3rqo7jJ7vH9la+HP5QGttx5SxnJpuVs+0sQAAAMCGtubBTVVdM/qz8HaS/xo/N6HvFUn+O916MS3JO9Z6vDvpV0fapyw0qurmVXWHqhpfqHiS26VbgDjpHoeaZeH8AUk2j527z4R+k5yd5NK+fe8J5xfqfLa19rVljGVSnSOyuKjx1LH0M5AWgqoj+llDAAAAQG89ZtzU2M+095f6OSXduixDsrCd9XeTnFNVj66q/0jyrSTnJbmwqr5QVc/rZ6FMcthI+9wlPm/0/GFj55ZVp7V2dbqFgX+gRj/GA9ZrLGPnN6XbdQsAAADorcfixM8fO35eutkzr063A9M0Ld2Mm68m+Zcpa63M24/3r9vSrXXz2xP63C7JiUkeVlX3b619Zez8gSPtaY8ULdg+5brR40uXscjv9nTbsd+yqvZurV3Zv39AFsO11RjLrtRZ1hbxVXXAEl1uvZw6AAAAMGRrHty01q4T3FTV8/rmq1pry/olfcB+qH89NMlPJvlOkq1ZXLvlJ5K8IN022ndOtw33fcbWn7nxSPuSJT7v0pH2+AyehTpL1ZhUZyG4We2xrLTOLNuX7gIAAAC7t3lsB/7Y/nWpmRi7g/36173T7Z70S2OLC59dVQ9K8p504c3PpFuEd3Stnn1G2tN2nVpw5Uj7hmPnFuosVWNWndUey0rrAAAAwIa27sFNa+2N6/2Za+iKLIY3p0zaEaq1dm1VPTNdcJMkj8x1g5srRtp7LfF5e4+0Lx87t1BnqRqz6qz2WFZaZ5bxx7PG3TrJWTtRD+D7Nm89feb5bSdtWaeRAACw0c1jxs2e5OIsBjf/MK1Ta+2/q+rL6XaPuseEGguWelRov5H2+CNIC3WW87jRtDqrPZaV1plqxhbjSZIf3C0dAAAAdj9zDW767bLvmeT26dZFuf5S17TWXrDW49oJ27O4CO5yFuG9bZJbjb0/et1SC+6OzjIZX+NlR5KfTrJfVd1siQWKF+pcMLIw8WqPZbTO2btYBwAAADa0uQQ3VXWrJH+a5GG7MIYhBTf/ncUZNEuFTgvnrx57f3SB5kOXqDF6/pwJdR460u8HHttKkqralOTgSTVaa5dU1fZ0YcpKxzKp36w6Vyc5f4m+AAAAsKFcb70/sKpunuRjSR6R5Abptp/emZ8h+chI++CpvTq371+/PPb+F5MsbBF+1BI17jtSY9vYuY+NtGfVOTyLjyedOeH8Qp07VtWsLbVHP2O8zllZXJR46liqaq8kRy5c01pbzsLKAAAAsGGse3CTbrvsQ9KFMO9P8oAkt0xy/dba9Zb6mcN4Z3lXku/17YdM61RVRyW5RX/40dFzrbWW5J394aFVdWQm6N9fmJ3yzv66UR9O8t2+/ZiavsjLcSPt0yac//spfUfHsm+Sh/eHn2mtnTd6vrV2cZJ/7g9/vqqmPXb1kCQ3mTEWAAAA2NDmEYT8SpKW5D2ttQe01t7fWvvmhCBi8Fpr30zyV/3hL1TVI8b7VNWNk7x85K3XTCj18iw+QvWKqrrOttj98Sv6w6vH6i2M5aokf94fHpbkGRPGcs8kx/eHZ7TWJu26dFqSz/ftZ1fVpJlEJye5+Uh7kpf2r5uSvKqqrvMoWVXtn+SP+8PvZPHPEQAAAOjNI7g5qH991Rw+ey08L8n/9O03V9Urquroqrp7VR2X5BNJ7tqf/8tJYUk/Y2Uh6Dg8yZlV9WtVdXhV/Vq6R5EO78+f3Fr73JSxnJxkYfbLS6rqNf1YjqyqZ6eb4bQp3bbbT5tUoLX2vSRPTXJtutkwZ1bV71TVEVV1/6p6R5In990/luTNU+p8MMnb+sMHJ/lAVT24/06PTbcGz8K9sLW19u0p3wkAAAA2rHksTnxJkr2TfH0On73qWmsXVNUD0j02dUiS3+l/xv11kt+dUeq56XacelySu2Ux9Bj1+iS/P2MsF1fVliTvTXKHJE/sf0ZdlOTRrbV/n1HnvVX1pCSvTPLDWZztM+oTSY5trV0zrU6673KTJA9McnT/M+raJC9srU2ahQQAAAAb3jxm3Hy6f/3ROXz2mmitnZNuVs0zk/xbkm+lW5x3R5L/m+RnW2vH97NZptW4trV2fJIt6da8+Upf4yv98QNba49vrV27xFjOTxf8nJBuG+7vJLksyWfT7eR1l9bae5bxnV6X5O5JXpfkC0muSPLNdLNsfivJvVprFy5R4/LW2pYkj07ygSTf6L/T9iRvSXLv1tqJS40FAAAANqpa76Vlqurh6WaTnNpae9i6fjgbRr8g8vYk2b59ew44YNr6yIzavPX0eQ8BdgvbTtoy7yEAADAwO3bsyIEHHrhweGBrbcdq1F33GTettbcneWuSY6tq63p/PgAAAMDuYt3XuKmq+6bbQehHk/xhVT0k3WMz56Z7pGem1tpH1naEAAAAAMMwj8WJP5xuO/AFd+9/lqNlPmMGAAAAWHfzCkFqTp8LAAAAsNuYR3AzviU0AAAAABOse3DTWjtjvT8TAAAAYHe07rtKAQAAALA8ghsAAACAgZrrDk1VdZMkD0tyzyS3TrJvkse11r400udHktwsyRWttS/MY5wAAAAA8zC34KaqfjvJHya58cJb6bb73m+s61FJ/i7JFVV1QGvtW+s3SgAAAID5mcujUlV1YpI/T3KTJFcl+eSM7v83yVeT7J3koWs+OAAAAICBWPcZN1V1tyR/0B/+bZKntNa+W1XXTurfWru2qk5J8tQkv5DkdeszUgCYbPPW05fss+2kLeswEgAA9nTzmHHzlHSPRX28tfYbrbXvLuOaj/evP7F2wwIAAAAYlnkEN0elW8vmlTtxzbb+9barPhoAAACAgZpHcHOb/vWzO3HNlf3r3qs8FgAAAIDBmkdwc1X/eoOduGYh7PnO6g4FAAAAYLjmEdzs6F/vtBPX/GL/ev4qjwUAAABgsOYR3Hww3eLEj11O56q6fZLj062L84E1HBcAAADAoMwjuHllkquT3KuqTpzVsaoOT/L+JDdKt87Na9Z8dAAAAAADse7BTWvtvCQvTDfr5g+q6t+q6lkjXR5QVSdU1T8n+bckt0s322Zra+2r6z1eAAAAgHnZNI8Pba29sKpukOQ5Se6R5PB04UySnDzStfr3X9Ba+/P1HSUAAADAfM3jUakkSWvt/yQ5MsmpSS5PF9KM/nwvyT8kuU9r7fnzGicAAADAvMxlxs2C1trZSR5WVZuS/HiSWyW5fpJvJvnv1trl8xwfAAAAwDzNNbhZ0Fq7Osl/znscAAAAAEOy7o9KVdXN1/szAQAAAHZH81jj5qtVdVpV/WpV7T2HzwcAAADYLcwjuNkryYOTvC3JN6rqr6vq56uq5jAWAAAAgMGaR3Dz6nSLD1eSGyd5TJL3JflyVb2squ4+hzEBAAAADM66BzettScnuU2SByV5S5LL0oU4t07ytCSfqKpzq+r3q+rg9R4fAAAAwFDMY8ZNWmvXtNbe21r7X0l+OMmjk7w3yTXpQpwfS/L8JOdV1cer6rer6pbzGCsAAADAvMwluBnVWrustfbW1tqD0s3E+e0k/9KfriQ/neTP0z1K9d45DRMAAABg3c09uBnVWvtma+0vW2v3TnL7JM9N8t/pApxNSe4/z/EBAAAArKdN8x7ANK21L1XVqUn2TXLbJDeb74gAYPk2bz195vltJ21Zp5EAALA7G1xwU1W3SfLIJI9Kcrex01es/4gAAAAA5mMQwU1V3STJw9ItUnzfdI9wVX/62iQfSvJ3Sf7fXAYIAAAAMAdzC26qaq90W4I/OskDk+y1cKp//VS6sOatrbWvrv8IAQAAAOZr3YObqvrZdGHNQ5LcZOHt/nVbkrck+dvW2rnrPTYAAACAIZnHjJt/StKyGNZ8M8nbk/xda+1fpl4FAAAAsMHM61GpK5K8K92jUP/YWrt6TuMAAAAAGKx5BDfHJTm1tXbJHD4bAAAAYLex7sFNa+1N6/2ZAAAAALuj661l8ar6k/7nVlPOX7+qDqqqg5aoc/uq+v+q6pNrM1IAAACA4VnrGTdPS7cQ8V8l+caE84cm+XSSa5cYyw2T3LWvBQAAALAhrOmMm51QS3cBAAAA2FiGEtwAAAAAMEZwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYqE3r9DlPrqpvTHj/VguNqvo/M66/1YxzAAAAAHuk9QpufmvGuda/Pm89BgIAQ7B56+lL9tl20pZ1GAkAAEO2HsFNrcNnAAAAAOxx1jq4OXqN6wMAAADssdY0uGmtnbGW9QEAAAD2ZHaVAgAAABgowQ0AAADAQAluAAAAAAZqvbYDBwZgOdsPAwAAMBxm3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAO1ad4DAAAm27z19Jnnt520ZZ1GAgDAvJhxAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4WSNV9ZKqaiM/91vGNQ+oqlOrakdVXdm/nlpVD9iJz923qp5ZVZ+oqm9V1SVVdU5VvbSqDtqJOneqqldX1flVdXlVXVBVH6mq36yqTTtR5xFV9b6q+mpVXVFV26rqzVV15HJrAAAAwEa17F/AWb6q+skkT9+J/pXk1UmeOHbqtkmOTXJsVb02yZNaa21GnYOTnJ7kjmOnDu1/Hl9Vj2qtvXeJ8Ryf5FVJ9h55e58k9+l/jquqB7XWvjmjxj5JTknyoLFTP9r/PKqqTmytvXDWWAAAAGAjM+NmlVXV9ZK8Ll0o9o1lXvaiLIY2n0ryyCRH9K+f6t9/YpKpIUdV3SjJe7IY2rwuyc8l+Zkkz01ySZKbJjmlqu4yo879k7w2XWjz9SRPTfLTSX4pyal9tyOTnNp/12len8XQ5kNJjum/0/FJPp/u3ntBVT1+Rg0AAADY0AQ3q++pSe6R5Nx04cVMVXVIkmf1h2cnuVdr7W2ttbNaa29Lcu/+/SQ5oZ9VM8kz0s2qSZJntdae2Fr7YGvt4621P0ryi0muTrJvkpdPGcumJK9Md19c1I/lFa21T7TW/rG19tAkf9F3v2+S/zWlzlFJHtUfvjvJL7TW3tl/p79OF/z8T3/+JVV1synfCQAAADY0wc0qqqoDszgr5reSXLWMy56exUfWntJau3z0ZGvtsiRP6Q83JXnahM+9QZLf7Q/PSfKy8T6ttY9nMUg6uqruPmEsxyY5pG+/uLX2+Ql9npnk2yPtSRaCqGuSPLm1ds3YWC5MckJ/ePN0s3AAAACAMda4WV1/keRGSd7YWvvwUgsS92vb/Ep/eG5r7V8n9Wut/WtVfTbdY1DHVNVTx9a6uV+Sm/XtN7bWrp3ykW9I8pt9+yFJPjl2/pixvpPGcllVvb2vc+equkNr7XMj3+lG6R7RSpIPtNZ2TBnLqelm9dykH8sPhE0AzLZ56+lL9tl20pZ1GAkAAGvFjJtVUlUPT7emy7cyfSbKuNulW4A4Sc5You/C+QOSbB47d58J/SY5O8mlffveE84v1Plsa+1ryxjLpDpHZHFR46ljaa1dlWQhqDqinzUEAAAAjBDcrIJ+jZY/6w9PaK1dsMxLDxtpn7tE39Hzh42dW1ad1trV6RYG/oEa/UyZA9ZrLGPnNyW5wxJ9AQAAYMPxqNTqeEmSWyf5lyxjQeIRB460pz1StGD7lOtGjy9trX1nGXXukuSWVbV3a+3K/v0DktQqjmVX6nxmif7fV1UHLNHl1sutBQAAAEMluFmhqrp3ksen27HpSWNrzyzlxiPtS5boe+lI+0ZT6ixVY1KdheBmtcey0jpL2b50FwAAANi9eVRqBapqrySvTTdT5U9ba5/eyRL7jLSX2oHqypH2DafUWc4uVtPqrPZYVloHAAAANjwzblbmOenWdPmfJM/fheuvGGnvtUTfvUfal4+dW6izVI1ZdVZ7LCuts5TxR7TG3TrJWTtZEwAAAAZFcLOLqurQJM/uD5/SWrt0Vv8pLh5pL/Wo0H4j7fFHkBbqLOdxo2l1VnssK60z04xtxpMk3U7rAAAAsHsT3Oy6p6ebUfKFJPtW1SMm9LnzSPtnq2phwdx390HPaPiw1GK7ozNMxtd32ZHkp5PsV1U3W2KB4oU6F4wsTLxQY7XGMlrn7F2sAwAAABue4GbXLTzmc/skb11G/z8Yad8u3cK8o7soHbrE9aPnzxk795kkDx3p96+TClTVpiQHT6rRWrukqranC1NWOpZJ/WbVuTrJ+Uv0BQAAgA3H4sTz9cUkX+nbRy3R977965eTbBs797GR9qw6h2fx8aQzJ5xfqHPHkdlBk4x+xnids7K4KPHUsfQLOx+5cE1rbTkLKwMAAMCGIrjZRa2141prNesn112w+OiRc9v6Gi3JO/vzh1bVkZmgf39hdso7J2w5/uEk3+3bj6npC7wcN9I+bcL5v5/Sd3Qs+yZ5eH/4mdbaeaPnW2sXJ/nn/vDnq2raY1cPSXKTGWMBAACADU9wM38vT/eoUJK8oqqusy12f/yK/vDqvv919LNV/rw/PCzJM8b7VNU9kxzfH57RWpu049JpST7ft59dVQdP6HNykpuPtCd5af+6Kcmrqur6Y2PZP8kf94ffSfJXU+oAAADAhia4mbN+xspC0HF4kjOr6teq6vCq+rV0jyId3p8/ubX2uSmlTk6yMPvlJVX1mqo6uqqOrKpnJ3l/uiDl8iRPmzKW7yV5apJr082GObOqfqeqjqiq+1fVO5I8ue/+sSRvnlLng0ne1h8+OMkHqurB/Xd6bLo1eA7qz29trX17yncCAACADc3ixMPw3CS3SvK4JHfLYugx6vVJfn9agdbaxVW1Jcl7k9whyRP7n1EXJXl0a+3fZ9R5b1U9Kckrk/xwFmf7jPpEkmNba9dMq5Puu9wkyQOTHN3/jLo2yQtba6+ZUQOAFdq89fSZ57edtGWdRgIAwK4w42YAWmvXttaOT7Il3Zo3X0m3wO9X+uMHttYe31q7dok656cLfk5Itw33d5JcluSzSf40yV1aa+9Zxnhel+TuSV6XbrvzK5J8M90sm99Kcq/W2oVL1Li8tbYlyaOTfCDJN/rvtD3JW5Lcu7V24lJjAQAAgI2sfnCdW9j99Ysib0+S7du354ADpq2RvLEs9TfvwMZjxg0AwOrYsWNHDjzwwIXDA1trO1ajrhk3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAM1KZ5DwAAmJ/NW0+feX7bSVvWaSQAAExixg0AAADAQJlxA3uIpf7WHAAAgN2PGTcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABmrTvAcAAAzX5q2nL9ln20lb1mEkAAAbkxk3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZq07wHAADs3jZvPX3m+W0nbVmnkQAA7HnMuAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZq07wHAADs2TZvPX3JPttO2rIOIwEA2P2YcQMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgdo07wEAAGzeevrM89tO2rJOIwEAGBYzbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZq07wHAACwlM1bT1+yz7aTtqzDSAAA1pcZNwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgbIdOOwmlrMVLgAAAHsWM24AAAAABkpws0JV9VNV9Zyq+oeq2l5VV1bVJVV1XlW9oarus5P1HlBVp1bVjr7Wjv74ATtRY9+qemZVfaKqvtWP55yqemlVHbQTde5UVa+uqvOr6vKquqCqPlJVv1lVy56tVVWPqKr3VdVXq+qKqtpWVW+uqiOXWwMAAAA2Io9KrUBVnZHkvhNO7ZXkDv3PY6rqzUke31q7akatSvLqJE8cO3XbJMcmObaqXpvkSa21NqPOwUlOT3LHsVOH9j+Pr6pHtdbeu8R3Oz7Jq5LsPfL2Pknu0/8cV1UPaq19c0aNfZKckuRBY6d+tP95VFWd2Fp74ayxAAAAwEYluFmZ2/avX0kXUHw0yf8kuX6Seyb5vb7Pr6f7s37UjFovymJo86kkL0ny+SQHJ3lWkrv15y9I8vuTClTVjZK8J4uhzeuSvC3J5UmOTvLsJDdNckpV3bO19p9T6tw/yWvTzcj6epI/TPJvSX4oyROSPCTJkUlOraqjW2vXTvlOr89iaPOhJH+W7s/qJ5I8p/9uL6iqr7bW/mraHwwALMdSa4FtO2nLOo0EAGD11IzJGyyhqt6T5E1J/l9r7ZoJ5/dPcmaSH+vfum9r7aMT+h2S5Jx04c7Zfb/LR87vm+SMJIcnuTrJoa21z0+oc2KS5/WHz2qtnTx2/p5JPtJ/zodaaz87ocamfiyHJLkoyU+Nf1ZVvSrJk/vDx7TW3jShzlFJPtwfvjvJsaN/Rv2fzSeTHJTk20lu31r7znidXVVVByTZniTbt2/PAQccsFql58bixAArI7gBANbSjh07cuCBBy4cHtha27Eada1xswKttQe11t4+KbTpz1+YbtbNgodNKfX0LM5+espoaNPXuSzJU/rDTUmeNl6gqm6Q5Hf7w3OSvGzCeD6ebhZMkhxdVXefMJZj04U2SfLiSQFRkmemC1sW2pM8q3+9JsmTx/+M+j+bE/rDmyc5fkodAAAA2LAEN2vvwyPtg8dP9mvb/Ep/eG5r7V8nFenf/2x/eEx/3aj7JblZ337jjMeX3jDSfsiE88dM6Ts6lsuSvL0/vHNV3WH0fP/I1s/1hx+YkTKemm5Wz7SxAAAAwIYmuFl7e420J4Upt8viWjlnLFFr4fwBSTaPnbvPhH6TnJ3k0r597wnnF+p8trX2tWWMZVKdI7K4qPHUsfSLNS8EVUf0s4YAAACAnuBm7R010j53wvnDljifKecPGzu3rDqttavTLXr8AzX6mTILi8Gs+VjGzm9KtwsXAAAA0LOr1Bqqqusl2Try1tsndDtwpL3UwkXbp1w3enzpMhb53Z7kLkluWVV7t9au7N8/IMnCI1irMZZdqfOZJfon+f7iw7Pcejl1AAAAYMgEN2vr6ekeG0qS01prZ0/oc+OR9iVL1Lt0pH2jKXWWqjGpzkJws9pjWWmdWbYv3QUAAAB2b4KbNdJvh31Sf/iNJL81pes+I+2rlih75Uj7hlPqLFVjVp3VHstK6wDAqtm89fQl+9gyHAAYGsHNGqiqOyU5Ld2f75VJHt5a+/qU7leMtPea0mfB3iPty8fOLdRZqsasOqs9lpXWmWX88axxt05y1k7UAwAAgMER3KyyqrpdkvcnuXmSa5I8srU2a5eni0faSz0qtN9Ie/wRpIU6y3ncaFqd1R7LSutMNWOL8STJD+6WDgAAALsfu0qtoqr6kST/lORHkrQkj2utnbbEZaMBxFIL7o7OMhlf42Whzn5VdbNl1rlgZGHitRjLSusAAADAhia4WSVVtX+SDyS5ff/WU1prb1rGpaO7KB26RN/R8+fsSp2q2pTk4Ek1WmuXZDE8WfOxjJ2/Osn5S/QFAACADUVwswqq6qZJ3pfkx/u3trbWXrXMy7+Y5Ct9+6gl+t63f/1ykm1j5z420p5V5/AsPp505oTzC3XuWFWzttQe/YzxOmdlcVHiqWOpqr2SHLlwTWttOQsrAwAAwIYhuFmhqto3yelJfqp/6w9ba3+83Otbay3JO/vDQ6vqyEn9+vcXZqe8s79u1IeTfLdvP6amL/Jy3Eh70mNcfz+l7+hY9k3y8P7wM62180bPt9YuTvLP/eHPV9W0x6UekuQmM8YCAAAAG5rgZgX6GSOnJblX/9aftdZ+fxdKvTzdo0JJ8oqqus622P3xK/rDq/v+19HPVvnz/vCwJM+YMN57Jjm+PzyjtTZp16XTkny+bz+7qg6e0OfkdIsvL7QneWn/uinJq6rq+mNj2T/JQsD1nSR/NaUOAAAAbFh2lVqZtyb5xb79wSSvr6o7z+h/1fjslCRprZ1XVS9NsjXdo0xnVtUfpwtQDk5yQpK79d1Pbq19bkr9k5P8WpIfS/KSqjokydvSbbN9dJLnpPvf/PIkT5tUoLX2vap6apJ3p5sNc2ZVvSjJJ9KFNU9I8tC++8eSvHlKnQ9W1duSPCLJg5N8oKpenu6xsJ9I8twkB/Xdt7bWvj3lOwEAAMCGVT/4xA3LVVU7+4f3pdba5im1rpfkdUkeN+P61yd5Ymvt2hljOiTJe5PcYUqXi5I8urX2nlkDraonJHllkr2mdPlEki2ttQtn1LhhknckeeCULtcmeWFr7cRZY9kV/eNZ25Nk+/btOeCApTa3Gr7NW0+f9xAA9njbTtoy7yEAALupHTt25MADv79x8oGttR2z+i+XGTcD0Ycxx1fV/0vyxCT3SLJ/kgvTLfb7mtbaPyyjzvlVdbckv53kV5Mcki582Z4u0Pmz1tqXllHndVX18SRPTfJz6bY4vzTdDlJ/l+SvWmtXzyiR1trlSbZU1aPSrZfzk0luluTrST6a5JWttY8vNRYAWC9LheSCHQBgvQluVqC1Nm0B4JXUfG+6gGUlNS5N8pL+ZyV1/itdiLQirbW3JHnLSusAAADARmNxYgAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgou0oBACzTUtuFJ7YMBwBWlxk3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFB2lQIAWEVL7Txl1ykAYGeYcQMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjbgQMArCPbhQMAO8OMGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQNkOHABgQJbaLjyxZTgAbCRm3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCiLEwMA7GaWWsDY4sUAsOcw4wYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEBZnBgAYA+z1OLFiQWMAWB3YcYNAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBsjgxAMAGtNQCxhYvBoBhENzAQCxnBxAAAAA2Fo9KAQAAAAyUGTcAAPyA5cwE9TgVAKw9M24AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyUxYkBANglSy1gbPFiAFg5M24AAAAABkpwAwAAADBQHpUCAGBNLPUoVeJxKgBYihk3AAAAAAMluAEAAAAYKI9KAQAwN3amAoDZzLgBAAAAGCjBDQAAAMBACW4AAAAABsoaNwAADJYtxQHY6My4AQAAABgoM24AANit2ZkKgD2ZGTcAAAAAAyW4AQAAABgoj0oBALBH8ygVALszM24AAAAABkpwAwAAADBQHpUCAGBDW+pRqsTjVADMjxk3AAAAAAMluAEAAAAYKI9KAQDAEuxMBcC8CG4AAGCFrJMDwFrxqBQAAADAQJlxAwAA68DjVgDsCjNuAAAAAAbKjBsAABgA6+QAMIngBgAAdhMetwLYeAQ3AACwhzBrB2DPY40bAAAAgIEy4wYAADYQj1sB7F4ENwAAwPd53ApgWDwqBQAAADBQZtwAAAA7ZTmzcmYxYwdg+QQ3AADAuvI4FsDyCW4AAIDBsYgyQEdwAwAA7HbM2gE2CsENa66qDkry1CRbkhyU5Mok5yd5e5K/aK1dNsfhAQCwh1rpWjyJ8AeYP8ENa6qqtiT5uyQ3HXl73yT36H8eX1UPbK19YR7jAwCAWTyyBcyb4IY1U1U/mW5Wzb5JLkny4iQfSnLDJI9I8oQkd0xyelXdo7V2ybzGCgAAu8IjW8BaE9ywll6eLrS5OskvttY+PnLug1X1uSQvSXJokv+d5AXrPkIAAFhjHtkCVkJww5qoqnskuV9/+Pqx0GbBy5I8NslhSZ5WVS9urX1vnYYIAAC7jdUIf5YiHIJhEtywVo4Zaf/NpA6ttWur6k3pHqG6ebqg5wNrPjIAAOAHrEc4tBwCJLguwQ1r5T7966VJPjmj3xkj7XtHcAMAABva7hIgWbia9SK4Ya0c1r+e31q7eka/cydcAwAAMFcrDZB2lwBqOSzCPV+CG1ZdVe2TZP/+cMesvq21b1fVpUn2S3LgTnzGAUt0ue1C46tf/epyy87V1RddOO8hAAAAe5gDnvzGPeZz/vU5P7fmn7ESY797Xn+16gpuWAs3HmkvZ4vvheDmRjvxGduX2/GII47YibIAAAAM0YF/Oe8R7JRbJvnSahS63moUgTH7jLSvWkb/K/vXG67BWAAAAGC3ZcYNa+GKkfZey+i/d/96+U58xlKPVe2V5NAk30hyQZJrdqL2ert1krP69j2SfG2OY4FJ3KMMnXuUoXOPMnTuUYZud7lHr59upk2SfHq1igpuWAsXj7SX8/jTfv3rch6rSpK01maundP7wnLrzVNVjR5+bZnfDdaNe5Shc48ydO5Rhs49ytDtZvfoqjweNcqjUqy61toVSRZW2p25iHBV3TyLwc2y160BAACAjUBww1o5p389pKpmzew6dMI1AAAAQAQ3rJ2P9a/7Jbn7jH5HjbTPXLvhAAAAwO5HcMNa+fuR9mMndaiq6yX5jf7wO0k+tLZDAgAAgN2L4IY10Vr7RJKP9ofHV9U9J3T7vSSH9e0/a619b10GBwAAALsJu0qxln433eNPN0zy/qr6o3Szam6Y5BFJntj3Oy/Jy+YyQgAAABgwwQ1rprX2qar6tSR/m+QmSf5oQrfzkmxprV084RwAAABsaNVam/cY2MNV1Y+mm32zJd324FclOT/JKUle2Vq7bI7DAwAAgMES3AAAAAAMlMWJAQAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW5gjqrqoKp6aVWdU1WXVtW3quoTVfWMqtp33uNjz1NVP1VVz6mqf6iq7VV1ZVVdUlXnVdUbquo+O1nvAVV1alXt6Gvt6I8fsFbfgY2rql5SVW3k537LuMY9ypqpqv2r6llVdWZVfa2/x75SVf9WVSdX1T2XUcM9ypqoqr2q6viq+seq+urIv/M/W1V/XVVHLrOOe5SdUlW3qqoHVdUL+v/mvHDk391v2IV6K74Hq2rfqnpm/7vWt/r/L5zT/y520M6Oab1Va23eY4ANqaq2JPm7JDed0uWzSR7YWvvC+o2KPVlVnZHkvsvo+uYkj2+tXTWjViV5dZInzqjz2iRPav5Fwyqoqp9McnaSTSNvH91a+/CU/u5R1lRV/WqSv0xyixnd3tlaO2bK9e5R1kxVHZjk9CQ/sUTXP03ye5PuMfcou6qqZt0Pb2ytHbfMOqtyD1bVwen+/3DHKV2+m+RRrbX3Lmdc82DGDcxB/wvI29OFNpckeW6Sn0nyc0le13e7Y5LTq+pGcxkke6Lb9q9fSfJnSR6W5Igk90zyv5N8uT//60nesEStF2XxX6KfSvLIvtYj++P051+4CuNmg6uq66X7Z+OmJN9Y5mXuUdZMVf1GkrelC22+keT5SX4hyd2TbEny1CQfSPK9GWXco6yJqtqU64Y2/5nkuHT/vv/FJC9Icml/7ulJnjGllHuU1bA9yft38doV34P971LvyWJo87p0v3P9TLrfwS5J9zvZKVV1l10c55oz4wbmoKo+lOR+Sa5Oct/W2sfHzj8zyUv6w+e11l6wviNkT1RV70nypiT/r7V2zYTz+yc5M8mP9W/dt7X20Qn9DklyTrpfos/u+10+cn7fJGckOTzdPX5oa+3zq/x12ECq6mnp/lb43CSnJXl2f2rijBv3KGupqg5L9wvD3kk+muSXW2vfndJ3r0mzF92jrKWqemiSd/SHH09yn/F/71fV3ftzN0jy7SS3aq1dPXLePcouq6rnJzkryVmtta9X1eYkX+xPL2vGzWrdg1V1YpLn9YfPaq2dPHb+nkk+0n/Oh1prP7v8b7p+zLiBdVZV90gX2iTJ68dDm97L0v2DKkmeVlU3WI+xsWdrrT2otfb2SaFNf/7CJL838tbDppR6ehYfV3nK6L9E+zqXJXlKf7gpydN2edBseP10/4W/SfutJFMf4RvhHmUtvSJdaHNhkodMC22SZMYjp+5R1tK9RtovnvTv/dbaJ9PNQkiSmyc5dKyLe5Rd1lp7XmvtPa21r6+gzIrvwf53qN/tD89J9zvW+Fg/nuT1/eHRfag5OIIbWH/HjLT/ZlKH1tq16WZGJN2/TO+3tkOC7/vwSPvg8ZP9s8a/0h+e21r710lF+vc/2x8e018Hu+Ivktwo3d/QfXipzu5R1lJVHZpuin2SvLIPvHe2hnuUtbbXSHvWWomjsxP2Xmi4R5m3VbwH75fkZn37jf3vWJO8YaT9kJ0d73oQ3MD6W9i159Ikn5zR74yR9r3XbjhwHaP/sTfpX263y+JaOWdMOD9q4fwBSTavbFhsRFX18CQPSvKtJM9c5mXuUdbSr460T1loVNXNq+oOVTVroeIF7lHW2nkj7dvP6LfwFzQtyedG3nePMm+rdQ/eZ0K/Sc7O4rpPg/y9S3AD6++w/vX80WeJJzh3wjWw1o4aaZ874fxhS5zPlPPuYXZKVd0s3SLaSXJCa+2CZV7qHmUtLWyf/N0k51TVo6vqP9KFi+clubCqvlBVz5uxuYB7lLX21iQX9e0Tqur64x2q6m7pFtJOkre11i4aOe0eZd5W6x5cVp3+d7KFGWiDvI8FN7COqmqfJPv3hztm9W2tfTuLye+BazkuSL6/c8/WkbfePqHb6L048x5Ot4vApOtgOV6S5NZJ/iWLz54vh3uUtfTj/eu2dGvd/G2S8V1IbpfkxCQfr6ofmVDDPcqa6oPu45Jcnm69m7Oq6jeq6siq+vmqel662Qd7Jfn3dDtLjnKPMm+rdQ8uHF/aWvvOMuvcsqr2ntlzDgQ3sL5uPNK+ZBn9F4IbW4KzHp6ebovFJDmttXb2hD47cw9fOtJ2D7NsVXXvJI9Pt0vEk9rObYHpHmUt/VD/emiS307ynSRPSnKrJPskuUeSf+j73Dnd9rLj/73tHmXNtdZOS7fbzuuT3DXJG9PtIvWBdMHiZekCm3u31r42drl7lHlbrXtwoc7O/N41qc7cCW5gfe0z0l7OzihX9q83XIOxwPdV1VFJTuoPv5Fu955JduYevnKk7R5mWapqrySvTVJJ/rS19umdLOEeZS3t17/uneSaJL/UWntNa+2C1tqVfeD9oCyGNz+TH1zo0j3Kmut303lUkl9O98/TcT+c5JGZvAGGe5R5W617cKHOzvzeNanO3AluYH1dMdLea2qvRQvT9C6f2QtWoKrulOS0dFspXpnk4TO2b9yZe3h0mql7mOV6Trrny/8nyfN34Xr3KGtp9P46ZdJOJ/2uJaOLaT9yRg33KKuuqvZL8k9JnpvkFukePT0s3f100yS/mORj6WaIvbuqfneshHuUeVute3Chzs783jWpztwJbmB9XTzSXs4UvIW/2VvO9D7YaVV1uyTvT7ft/DVJHtlam7Xq/s7cw/uNtN3DLKnfavnZ/eFTWmuXzuo/hXuUtTR6f/3DtE6ttf9O8uX+8B4zarhHWQvPT3Lfvn18a+2E1tq5rbWrWmsXtdY+kOToJB9KNxvnT6pqdK0m9yjztlr34EKdnfm9a1Kduds07wHARtJau6KqLky3QPEBs/pW1c2z+A+Q7bP6wq7oF838pyQ/km4r0Mf1z8TPMrpA3Mx7ONddIM49zHI8Pd3fin0hyb5V9YgJfe480v7Zqrp13353H/S4R1lL29Mtmp0sb8HM26Zb/2aUe5Q1U1WV5LH94XmttTdO6tdau7qq/iDdzJvr9dc8vT/tHmXeVuse3JHkp5PsV1U3W2KB4oU6F7TWrpzRby4EN7D+zklynySHVNWmGVuCHzp2Dayaqto/3QKFt+/fekpr7U3LuPQzI+1Dp/b6wfPuYZZjYZry7dNtZ7uUPxhp3y7dwoLuUdbSf2dxBs0PbLE8ZuH8+L/n3aOspR/O4iLan1qi7ydH2qP3mnuUeVute/AzSR460u8HHm9NkqralOTgKTUGwaNSsP4+1r/ul+TuM/odNdI+c+2Gw0ZTVTdN8r4sbmu7tbX2qmVe/sUkX+nbR83qmMVp2l9Ot3UurAf3KGvpIyPtg6f26iwE418ee989yloaDQqX+kv6G0y5zj3KvK3WPfixkfasOodn8UmHQf7eJbiB9ff3I+3HTurQbx36G/3hd9I9gwwrVlX7Jjk9yU/1b/1ha+2Pl3t9vy3zO/vDQ6vqyCmfc2QW/wbknTu5nTMbVGvtuNZazfrJdRcsPnrk3La+hnuUtfSuJN/r2+O7RX1fv1PfLfrDj46ec4+yxr6V5KK+fc9+JsE0o7/IfnGh4R5l3lbxHvxwku/27cf0jxJOctxIe6llA+ZCcAPrrLX2iSz+R9zxVXXPCd1+L93q/0nyZ621703oAzul32b5tCT36t/6s9ba7+9CqZdn8W/mXlFV19kysT9+RX94dd8f1tPL4x5lDbTWvpnkr/rDX5i0DlNV3TjXvadeM6HUy+MeZQ30u5qd3h/+SLqdpX5Av5bi6F/cvGesy8vjHmW+Xp4V3oOttauS/Hl/eFiSZ4z36X8XO74/PKO1dtZKB74WSjAK66+q7pZuGt4N061a/kfpZtXcMMkjkjyx73peksNbaxdPqgM7o6r+Xxb/hviDSZ6WblHiaa5qrZ03pdaLk2ztDz+V7j/+Pp/u0YETktytP/fi1tpzVjZyWFRVJyZ5Xn94dGvtw1P6uUdZE1V1yyRnJzko3S8Lr05yarpZDj+R7v5a+Bvgv2ytPXlKHfcoa6Lfoe+TSfbt33p3kjemW/h9nyRHpvtvgIP68//cWvv5CXXco+ySqrp3kkNG3to/ycl9+8wsBuBJktbaG6bUWfE92IfpZyf5sf6t1yZ5W7otv49O8px0u05dnuRnWmv/voyvuO4ENzAnVfXLSf42yU2mdDkvyZbW2vnrNyr2ZFW1s//A/1JrbfOUWtdL8rokj5tx/euTPLH/2z9YFTsR3LhHWTNVdVi6x6YOmdHtr5M8adqsWfcoa6mqfj7dIu/7L9H1g0ke1lr79oQa7lF2SVW9Icljltu/fxR6Up1VuQer6pAk701yhyldLkry6Nba+MyzwRDcwBxV1Y8m+d0kW9JtdXdVkvOTnJLkla21y+Y4PPYwqxncjNR8YLoZYvdI9x+HFyY5K8lrWmv/sCvjhFmWG9yM9HePsiaqar8kv5XkYel+GbhRkm+k+9vk17TWlrU+nXuUtVJVt0j3CMgvJblTkpulmyX2tXT32FuSvGuptWnco+ys1QpuRuqt+B7s/5n920l+NV3ovle67cPfm275gC8td7zzILgBAAAAGCiLEwMAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYKMENAAAAwEAJbgAAAAAGSnADAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAANQVZurqlXVG+Y9FgBgOAQ3AAAAAAO1ad4DAAAgSfLlJIcl+e68BwIADEe11uY9BgAAAAAm8KgUAMAuqqr79evSnFhVP1NVH6qqi6vqgqr6i6q6Yd/vAVV1ZlVdWlVfr6o/rqrrj9WauMZNVT20qs6oqm9U1RVVtb2q/rGqjpkwnmX3XcZ3e0M/nvGfzTtba08x8r/3+M+J8x4bAHsuwQ0AsGFV1W9N+UV81s/jJ5T66ST/nO4xp9ck+Z8kv5XkdVX1q0lOTbK9P/edJM9KsnU540vyjiR3SHJakj9J8k9JDkxyzK72XWtVdb2qOqaq/rKq/qMPq66qqouq6gtV9c6q2lpVP7aGY3jNyP9mR+/ktT83cu0r12qMALAc1rgBADayu+7CNWdNeO8BSY5prb0zSarqBknOTvKoJPdPclRr7az+3POSnJ/k6VX1x621q2d81uOTXJXkJ1trF4yeqKpbrKDvzvhK/x0WfHlW56p6YJKXJTl0wukbJLlxktsleXCSF1fVGUme01r7lxWMcZI3JXli3/71JB/aiWv/10j7zSPts5L8xMjxp3dtaACwfIIbAGAju2v/ekGSn13mNZ+Z8N6HF0KbJGmtfa+q3pHkLknevRDa9Ocurqr3JHlckgOSbFvi877X/1xHa+2bK+y7XN9rrf3XcjpW1QlJXpyk+rfOTPLuJJ9K8s0k+yT54ST3SrIlyR2THJXk/6QLv1ZNa+3Mqvp8koOTPKyqfru1dvkyvsMNkzy0P/xsa+3fRmpemuS/Rvqu5pABYCLBDQCwIVXV9ZLcuT/8z+WGE1N8asJ7X+1f/33GudtmdnDz9iQnJfmvqnpbkg8n+Vhr7Tsr7Lvqquo3+s9PkguTPLq19v4p3U+tqmck+eV0Qc9aeVOS56eb5fMrSd62jGuO6fsn151tAwBzYY0bAGCjumOSffv2SkKbJLlowntXL+PcDZao+5IkT0jytST/O93slQv6NWJut4K+q6qqbpvk1f3hpUnuOyO0SZK0zruS3D3J69doaG9OsrCF6q8v85qFfi3J3676iABgJwluAICN6q4j7UGuVdKHG3/VWjs8yS2THJtuoeMHJzl9dGeqnem7Bv53khv27d9vrZ2z3Atba1e01k5Zql9VHVFVr6uq86rqkn6HrnOr6lVVdYcptb+Y5GP94S9W1a2W+IwfTvIL/eEZrbUvLfd7AMBaEdwAABvVXUfaK51xs+Zaa99srf19a+3XknwwyWFJDllp35WqbqGX3+gPL8kqz56pqk1V9RdJ/i3dAsx3SLJfutlSd0zy5CT/XVVPmFLiTf3rpiSPXOLjHpnFpQTeNKsjAKwXwQ0AsFHdtX9tSf57juOYqqruX1Wbxt67QZIf6g8v35W+q+xOSfbv2x9trV28yvVfn25r9ST5h3Q7Ph2R5B7pHg3773SPnL22qn55wvVvz+J3X+pxqYXzl6XbWh0A5s7ixADARvWT/evXk2xexg5BX13h7ky74v8muayqPpbkS+kCil9I8uNJ/m9r7X92se9qustI+/9bzcJV9dAszuZ5Qmvtr8a6nF1Vf5vk9HS7gv15Vf3D6BbrrbWLquqdSR6R5O5VddikR7mq6seT/FR/+PdrEEABwC4R3AAAG05V3TrdttRJcussb42bhydZci2WVfbsdNtkH5FuB6ZLk5yf5DeT/PUK+q6m/UfaF8zqWFV3yuJW4eO+2G+3PerZ/etpE0KbJN0aOVX1O+m2ad+c5H5J/mms25vSBTdJN6vmORNKjc7G8ZgUAINRrbWlewEA7EGq6peSvHcnL7t9v9jthlBVb0jymCRfaq1tntHvuUle1B8+vrU2dY2bqro6ybRFko9urX14pO9tk+zoDx/ZWpu5lXdVXZAuRPo/rbUXjp27fpLtSW6T5H+SbG4j/xHcr9PzpSQHptuq/cDW2jWzPq+/bqHG81trJy7VHwB2hTVuAICN6K4j7V9ordUyfjZMaLOTRh8p2m8V6x4+0n5rVbVZP1mc+XPr8UJ9CPOW/vCgJEeNdblfutAmSd6ynNAGANaL4AYA2IjuOtIe5Fbgu5HRdX9uOatja23TaBiW5Pkzus/cunuGfae8/8aR9vgixR6TAmCwrHEDAGxEd+1fv9Fa+/o8B7IH+I+R9k9N7bXzRh+penSS/1zmdd+e9GZr7dNV9R/pFqV+WFX9Tmvt8qq6YZKH9t3+o7W23M8BgHUhuAEANpSq2jfJIf2hX9JX7r/Tzbq5RZL7VNV+ExYZ3hWjM3laa+2/VqHmG5P8SZKbJHlwup24fqU/Tsy2AWCAPCoFAGw0P5nF/wba6eBmqbVWJqy9skfrF/ldCDxunOS4VSr9qZH2L65SzbckWdgq/NfHXkfXwQGAwTDjBgDYaO460t6VGTeT1mV5XpLvJnn5LtTbE/xJkicluWGSP6qq97XWzl9Jwdba+VX1mSQ/nuQRVfW81tr/rLDm16vq/UkemOT+VXXnLIZC72+tfW0l9QFgLQhuAICN5q4j7Z0ObiZt+1xVz0vynY26JXRrbUdV/XaSv0732NFHq+qRo9t7T3HzJc6/KN0smH2SnFpVv9Rau2BSx6raO8nxSf66tXbFjJpvTBfcbErytiz+97DHpAAYJMENALDR3LV/bUmu18+6WMoXWmuXrd2Qdn+ttb+pqtsmeUG6Lbk/VFUfSfKudAHZN5NUut2ifjLJsUmOGClx+YSab62q+yd5TJK7J/lMVb0myRlJLki3/fjBSe6T5CFJfihLBzDvSvKdJDdLcqf+vYuSvHOnvjAArBPBDQCwYVTV9ZIsBDWV5OxlXnpQEsHNElprL+p3bnpZkjskuW//M8uZSU5orf3blPPHJ/l6kt9Lsn+S5/Y/k1yabq2aWWO8oqpOSfKEkbdPaa39QHAEAENgcWIAYCP5sST77uQ132itbV+LweyJWmvvTnJYuhkwr03y6XSzY65OcnGSLyV5b5ITk9yptXbv1tqZM+pd01o7Id1aNy9Lt2jxt9MFNBen29Xq79LNyrnNMgOYN44de0wKgMEy4wYA2DBaa+emm2nDGmqtXZPktP5ntWqel+QZq1TrzLgPANhNCG4AAJjlBmPrAH22tfa9uY1mjqpqvyS3m/c4ANhYBDcAAMzyI+ked1pwuyTb5jOUubtHkg/NexAAbCzWuAEAAAAYqGqtzXsMAAC7tapqSb7UWts877EAAHsWM24AAAAABkpwAwAAADBQghsAAACAgbLGDQAAAMBAmXEDAAAAMFCCGwAAAICBEtwAAAAADJTgBgAAAGCgBDcAAAAAAyW4AQAAABgowQ0AAADAQAluAAAAAAZKcAMAAAAwUIIbAAAAgIES3AAAAAAMlOAGAAAAYKAENwAAAAADJbgBAAAAGCjBDQAAAMBACW4AAAAABkpwAwAAADBQghsAAACAgRLcAAAAAAyU4AYAAABgoAQ3AAAAAAMluAEAAAAYqP8fdtBDuSCYhk8AAAAASUVORK5CYII=\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'$E_\\mathrm{T}^\\mathrm{miss}$ [GeV]')\n", "plt.ylabel('Events')\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 }