{
"cells": [
{
"cell_type": "markdown",
"id": "bb99e077",
"metadata": {},
"source": [
"# tprofile2polyRealisticModuleError\n",
"Simulate faulty detector panel w.r.t. particle charge\n",
"\n",
"\n",
"\n",
"\n",
"**Author:** Filip Ilic \n",
"This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, March 19, 2024 at 07:11 PM."
]
},
{
"cell_type": "markdown",
"id": "be089e92",
"metadata": {},
"source": [
" Arguments are defined. "
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "97aead64",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:55.002530Z",
"iopub.status.busy": "2024-03-19T19:11:55.002114Z",
"iopub.status.idle": "2024-03-19T19:11:56.306242Z",
"shell.execute_reply": "2024-03-19T19:11:56.272315Z"
}
},
"outputs": [],
"source": [
"Int_t numEvents = 1000000;"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "e1f263c1",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:56.327220Z",
"iopub.status.busy": "2024-03-19T19:11:56.326820Z",
"iopub.status.idle": "2024-03-19T19:11:56.642263Z",
"shell.execute_reply": "2024-03-19T19:11:56.641110Z"
}
},
"outputs": [],
"source": [
"TCanvas *c1 = new TCanvas(\"c1\", \"4 Malfunctioning Panels\", 800, 400);\n",
"c1->Divide(3, 1);"
]
},
{
"cell_type": "markdown",
"id": "04690b30",
"metadata": {},
"source": [
"-------------------- Construct detector bins ------------------------"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "311929eb",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:56.649883Z",
"iopub.status.busy": "2024-03-19T19:11:56.649502Z",
"iopub.status.idle": "2024-03-19T19:11:57.099175Z",
"shell.execute_reply": "2024-03-19T19:11:57.086405Z"
}
},
"outputs": [],
"source": [
"auto th2p = new TH2Poly();\n",
"auto avg = new TProfile2Poly();\n",
"auto err = new TProfile2Poly();\n",
"\n",
"std::ifstream infile;\n",
"TString dir = gROOT->GetTutorialDir();\n",
"dir.Append(\"/hist/data/tprofile2poly_tutorial.data\");\n",
"infile.open(dir.Data());\n",
"\n",
"if (!infile) // Verify that the file was open successfully\n",
"{\n",
" std::cerr << dir.Data() << std::endl; // Report error\n",
" std::cerr << \"Error code: \" << std::strerror(errno) << std::endl; // Get some info as to why\n",
" return;\n",
"}\n",
"\n",
"std::vector> allCoords;\n",
"Double_t a, b;\n",
"while (infile >> a >> b) \n",
" allCoords.emplace_back(a, b);\n",
"\n",
"if (allCoords.size() % 3 != 0) {\n",
" std::cout << \"[ERROR] Bad file\" << std::endl;\n",
" return;\n",
"}\n",
"\n",
"Double_t x[3], y[3];\n",
"for (Int_t i = 0; i < allCoords.size(); i += 3) {\n",
" x[0] = allCoords[i + 0].first;\n",
" y[0] = allCoords[i + 0].second;\n",
" x[1] = allCoords[i + 1].first;\n",
" y[1] = allCoords[i + 1].second;\n",
" x[2] = allCoords[i + 2].first;\n",
" y[2] = allCoords[i + 2].second;\n",
" th2p->AddBin(3, x, y);\n",
" avg->AddBin(3, x, y);\n",
" err->AddBin(3, x, y);\n",
"}"
]
},
{
"cell_type": "markdown",
"id": "2f879a08",
"metadata": {},
"source": [
"-------------------- Generate particles ------------------------"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "8d48c652",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:57.130746Z",
"iopub.status.busy": "2024-03-19T19:11:57.130370Z",
"iopub.status.idle": "2024-03-19T19:11:58.794142Z",
"shell.execute_reply": "2024-03-19T19:11:58.792461Z"
}
},
"outputs": [],
"source": [
"TRandom ran;\n",
"for (int j = 0; j < numEvents; ++j) {\n",
" Double_t r1 = ran.Gaus(0, 10);\n",
" Double_t r2 = ran.Gaus(0, 8);\n",
" Double_t rok = ran.Gaus(20, 2);\n",
" Double_t rbad1 = ran.Gaus(1, 2);\n",
" Double_t rbad2 = ran.Gaus(2, 0);\n",
"\n",
" Double_t val = rok;\n",
" // -------------------- Malfunctioning panels -------------------\n",
" if (th2p->IsInsideBin(4, r1, r2)) val = rok - rbad1;\n",
" if (th2p->IsInsideBin(20, r1, r2)) val = rok - rbad2;\n",
" if (th2p->IsInsideBin(13, r1, r2)) val = rok + rbad1;\n",
" if (th2p->IsInsideBin(37, r1, r2)) val = rok + rbad2;\n",
"\n",
" // -------------------- Fill histograms ------------------------\n",
" th2p->Fill(r1, r2, val);\n",
" avg->Fill(r1, r2, val);\n",
" err->Fill(r1, r2, val);\n",
"}"
]
},
{
"cell_type": "markdown",
"id": "212eefcb",
"metadata": {},
"source": [
"-------------------- Display end state ------------------------"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "4f03685a",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:58.800969Z",
"iopub.status.busy": "2024-03-19T19:11:58.800395Z",
"iopub.status.idle": "2024-03-19T19:11:59.023949Z",
"shell.execute_reply": "2024-03-19T19:11:59.022769Z"
}
},
"outputs": [],
"source": [
"c1->cd(1);\n",
"th2p->SetStats(0);\n",
"th2p->SetTitle(\"total hits\");\n",
"th2p->Draw(\"COLZ\");\n",
"\n",
"c1->cd(2);\n",
"avg->SetStats(0);\n",
"avg->SetTitle(\"average charge\");\n",
"avg->Draw(\"COLZ\");\n",
"\n",
"c1->cd(3);\n",
"err->SetStats(0);\n",
"err->SetContentToError();\n",
"err->SetTitle(\"error\");\n",
"err->Draw(\"COLZ\");"
]
},
{
"cell_type": "markdown",
"id": "5baadd9a",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "63ecd29e",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-03-19T19:11:59.028918Z",
"iopub.status.busy": "2024-03-19T19:11:59.028577Z",
"iopub.status.idle": "2024-03-19T19:11:59.479093Z",
"shell.execute_reply": "2024-03-19T19:11:59.477720Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAF0CAIAAAAmXahBAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO29e6xcx33nWS2KtLP2iI+YJDSIbMCgxBlHoQJQGKztmH2oLAP6MRiPiVlaBr0LONoFjAwYWDuDgRcY3b4OEAdwQmOFHQRYMAE2IWJxgOvx7tgJdzjS7UvDziIQd0M6tvdKXAwsBushJfOh+Hk1nN4/frqVuvU6v6r6VZ065/w+uBBazT6PPn2q69PfX506k8ViIRiGYRiGYZg07ut6BxiGYRiGYYYASxXDMAzDMAwBLFUMwzAMwzAEsFQxDMMwDMMQwFLFMAzDMAxDAEsVwzAMwzAMASxVDMMwDMMwBLBUMQzDMAzDEMBSxTAMwzAMQwBLFcMwDMMwDAEsVQzDMAzDMASwVDEMMwQmk0nXu8AwzNhhqaJnPp9H/JPWJcxms9lspi3bNI1/JQzDMMyQmM/nk8lkMpnIHgF6h6Zp5vO5fCCEaJpmMplANwEvg2W72e+xwlJFz9GjRyP+qZWmaaBRQVOJXg/DkDCbzdTvevW0lI/ltzz8L3QGk8lkPp+/8MILWh8gX7+0tCSfNF+jvViuXO6SfDH8DlGfgRdDNyNfLH/AyNdrv2cYpivm8/nRo0cXi8VisVheXoZTfXl5Gf5pPp8vLy+DV00mk2eeeWaxWEynU/nzezabLRaLTt/B+FgwpEynUyHEdDpdLBZLS0twkJeWluT/wj+trq6qr1wsFtpnsbS0NJ1O4TWw+Orq6tLS0urqqvzg5ErgBQxTDDgP4bF6NqrPwDmsvhiawOrq6mKxmE6nzz///GKxeOaZZ6AzEELAM3BWwwNtDRK5clgnLC4bGjyYTqewLfUZc4VyQXXfuE0xNQDn+erqKrQvOHvlqStbwWJrDyKbD5zPTEk4qSIGfknI3xBwlOEXhvw1LIQ4evSo/A3h+lm8trY2n89hcbnypmmgqcBK4AUcXDGFaZoGTjw141lbWxNCzGYzsKLl5WU1RoIHMoWaz+dPPPHECy+8AEsBTzzxhBDimWeegf9dW1sz1wDIX+3qb3HYGbkI7BU0Rrmg3GHYSSEEPFC30jSNbHQM0zlwDstihYo1xGU6hKUqF+q39tLSktYYoBuYz+dqj6Khfum7tMnsMximAFBEg7qDfBJOVGla8pXz+Vyt6AFN04CHyfPcsy3rGlqBTcjBiOa/ag/ktoQQMhJmmA6Bk9McYmvlhRdekP9luuL+rndgyHj6gMlkApVvT4/S2oUsFgtobGtrawsunDMFAcvRvujl2QinLpzh8Ng8meVJu7S0pP20kP/rWcN0OpU9zWQysZ7/chPmMMSmaSAtFkIsLy9LY5M5lkvFGKYkcBLCKEA4512vfP7553/1V391Op2ura09//zzpXaQ0WGpyoX6ra19QcvahEgLb9WLPpA/ZRiGBDi9IWqdTqcwlhaSJzmUcK5ceWQGPyBMa2trzzzzzNra2gsvvPD888/Djw35Gs8a4J9gB1ypktwEeJvqVfB7RvZV8OTS0pLcHP9KYSrBLERo9W7giSee0E5armB0Q+ExXGNAbA4/V+t38E/qSEM5UF1sNhJ1JXJo7WJzvKEcCCxH5mprYJheA6PUAdlk8iHH8ML1H7k3xzDMGLDH5gzDMIWRRcPPfe5zzz//PAxazwfM3SCEgAtKsm6LYZiRwFLFMEwtwBjb3DolSS/BMwzDqLBUMQzDMAzDEMBTKjAMwzAMwxDAUsUwDMMwDEMASxXDMAzDMAwBLFUMwzAMwzAEsFQxDMMwDMMQwFLFMAzDMAxDAEsVwzAMwzAMASxVDMMwDMMwBLBUMQzDMAzDEMBSxYwOeWv3+XzOtyhhGIZhqODb1DDjommatbU1OO2bpoG7vzEMwzBMOpxUMUNmNpuBNqnp1HQ67XCXGIZhmKHCUsUMGZCqpmmkVEm49scwDMPQwlLFDB9rmY9rfwzDMAwtLFXMkJltYvUna4LFMAzDMHHwQHWGYRiGYRgC7idcF0QCQhkUDA8mkwnhVhgGWCwWGxs7/a95y1ter+FnAzcNpjA/+9kD/hdA04DzsGkaqIZrV3UUwNU0BLcOJg+tTUOkdRxkSZXnSvXJhPOwLew/cdl88sbKYXhw8Mj5O3sPeF7AiM2TqhdSxU0Dz6EL56zPXz1+Cv714W23zH9dOXY67271jclkgpGqpaUlIcRsNoPzEP6rKU5W/PObcOtQOXjkvPX59Usn4YG1WxHccWwF0zREWsdBllSp11Ktra3BT5+RD1hxneW0a+M2g0H+IBabv8tl55H7LOWmYeKSpzhOXHzW+jzLluQtb3ndfBJORSEETDKyurpaeK+4aZi45Cka7jgKk2Wg+urqqhYjTwxybLc2+KythKZplpeX5WM4Oefz+Ww2KzyrAjcNAJInphg/+9kD2p/YbAvSY9RfGuX30GwaYpStQyZPTE+hHFOloca5I0lx5W+C8jp18Mj5m48cFkLcPmspHY4cbUoqcCnoQjrpP0bYNGQ01YlOnbj47Mv39nS19WqB3xXyVIQW0e3kbVoRcCStY/dT1+BB+W/vDvusoZJFqtTZq3OsvzZoy3xCiF2vXrMOq8LQYfvsHGuNQ2U+n6+trcl5Fors1BbG1jRoy3xCiJfv7bEOq0IC+8NqBciASihNQ2z+6uhkZ8RomoZQvqs1bj5yeN9LxH1KK9CLsVqlQylVsjGMp1W06lSKHiUCLXZUamWOQNQ0q2ma6XQKv8hLnqUjbBrkOkVIt7FZPWhnYyex0AibhsulaoCDq3Qylv8GDHk0lY8xB1cu5E9zKAuO59u8ADW7lAkHV0xJatYpDQ6uomGpCqNHOqUxcruy/iBmoyKkXzqlwsEVk5UIl+qkAmjCwVUELFUoUlyqwwqglRGWBZl89NelTDi4YgjpUS7VCgdXePjefyiiT6bo62N3vYpqkHyWM90SrSDRC8J1fK3wJFVMt0T/cC3wi5c7jnywVGFBWg5QuEpIPl8cw+AJuhyvfLLlmhqUYWpDDbdgipys9Hc0S80Ukqo+Tt22+6lrVPltkJCVhPA9MnH0sWkIIQ5dODekwp/JiYvPso11Th9bx/4TlwcvK2N4j9EUkqrFJmU2l45VNap1o1C0eJm9qkN61zSEkTatHDudMndUVWhFSfaqbuld69BUY0hDV82KIXuVFR6oruMyjPVLJwdcZePR6wyGYadTVsCreIQW4yeHYVRyDaAHHsBuwmOqtmAalfZMdFgVsWDuYMx0xN5FVpMfPOD/63oHB4VpVOozJcMq5Fj1aMx32rvIiptGSUyjGnCK0/c329o0ElsHJ1Vv0jufMNn30uX0sY27n7rGeRWjMcKAyoQjK8akXz6Rj/0nLnNeBXBSJUSgUfV9ZFWrM/HodUYFY1SgGn0fWYWZ5aF3kRWTD4xRJf5GLXANIAaMMPHodYClKsCo4iadkkuR21jWXwbsVYzInFHlm2Yza57EFwYyInNGlbVckLXjYK8atVQhIxnzNRg9KnZuxQ2fxyzFkdWYQU6aoL0GE1aVLCZG2A9y99irRgsykonrAsyv3ExhVdzuId/4mNVqvFLVC10gD7cifgD14kAxtNQ/iCp0rHqr7UXEZhxZjZAx60IQoz1QY5QqkgAmTnd6Oh6LvWo8RM/qKStufR9ZFQp71UiIDmD6fulPdLlwnF41xhnVU9Ka6Hv5VcWAJ9zqHVU1DcG3E2ZqoqrWEeEWVFZRyXD1cUpSKCOdUZ1kRrWexk5MVdTWNCIwk63+hlX11z1HRW2tg7/zR8VsNpvNZvP53HzSs9S4yn97nl6HB+uXTnY1U22HzZLkrunyGDJDQpWJ6LCqv3M4peRz8l2zkA0VmdCsXzoZ9wXe9wpgHDLb62PEBeY0m82OHj0qn2yaZjabwX9dC45IqkwbSPeqAoZU4W8j9qqBkcMGCoRVuedVj4C9aniYQtDJd3JXFUCq+Rd651XSnKbTqfYkeJVrwbFIlfQANawKWkPKsKocI7HgXI8YHRW6iPobSy7LXjUYpAd0IgQ5hnDJ6KjAEHIzpmKvGhJSBdSwKm4NoWQKt6DjiNir0EWkjcleo1qv2vF3/1r7E0I0TTOfz7VQSisFWhmFVHkMIGtY5TqHKgyfImCvGgAuAwh1Het6PGFVneZBtVd1vjsmFI8EZP0OH/DV1nV61cb/9wvanxBiPp/DgCo1lFpbW4MnR13+M/t+7ZnKbwNOhfm7JzRPPnjk/KuP/7L6DHtVr/HfIDmOPg6rSk/LchxJpltabxsc6lXpyVMl1wAGYRZG6vQqE9CmpmlAquC/S0tL8IxHqgZ+Q2V/r79+6SR85MhbEbtuNrzr1Wt39lYxDpHknsoqrV8Ee55ev3XmIOEWg7jv1Z1tL/nrEvvRQ1p7/avHT6WbwcPbblU48okKv0EeunCuw/kpuGmk4O/4ZceB/Obv782Go3dbXdA8RN0eEETTEEL8tVbpg//1X/f35vpjdqonBOUoQXmVWVYP/dUS9HrCqJlnqGJE2Rwl64j1yo2N86o+EhSlDGMsRyvk8VJf8qoIBjv5p9+oooer95S45Nk6RN2Ei4ApVDW9oUjwgJ4KhAyTQvcfOZNCTw9LJZRvHf7OPnq4OhW9qACaQ9StDNWrhjn5Z1w3nzi4aiQ/WaywV0VTfnpDfDcfV70yi2L9nQuUBPaqaAq3jrhuHv/NP87ZqjwM0qsGWP5DdvDWsKrVq6Kvy0ivGKYTVPszYyptlLoGe1UvwHTwhSWgjzfGCZ3wk72qfpAdvDWsav0yJ7xpYLGwKm7qTmRMJRmeVw1NqtK7dmReZY1/zaaV6YyJnqQqK+xVlRPRtVMZjxlWZfKMkpNUBcFeVTPpX9S0P5IzzaoQPUlVNMiB/AX2pBiDkqrQTt31+gonWUhpsYmZc5C6sVdVS1CnHmcAcqmSEyukjFWPHlAVt1SFXqXd3Qwzt+HwCO3UUyYgHEkFMPQH/5C8ajhSldKdE4467OqyPlpSWj57VYXU0J33aGSVf1dTlLGGD0Ki3d0MHozNq1K6cx6urpEyUcJgvGogUkXekXvCqp5Od1tbrZApRmJH3scxT1Z6ZHVl0O5uNpvNlpaW1BcMXrDIO/KIG2zUSa/vhdwtA5GqW2cObv/eDyMW9Myt0FoEdP1MqTZ/QmKdScE/Sp2pmTgr6uldlnMTOkS9Eu5/7Nvan3Dc3WxU3Fg5vLEbMxWkjmduhb5//0cTOkR9qAxEqoQQN1YOx3mVhwKDq3K3wPIl/A4nWGdcJKZN+MU7GVYVR+KAqsSNFuY/XflF7U847m42Nm6fPRDnVR783+ok38nVVgCj6emk8ya9l6p3LH9H/d8Ir4qYCLS1AphDlcr8AEJO+Mn0gkRd6EtYVWxe9Z7GVC60u5upzOfzYcuWVtWK8KqIiUAHXEqzxlSV3L2tMP2eUR2MSnoVfK60eRVhWFVmVKN1h4vpkRpTab7LWMk3ZzR0/LL7r3loVM37Rkgn2ZiHuQI8M5vN5O1jazCqTK0D5EYqDvySpM2rqH4Dd3K1YPkBVWpM1Xf1HNqM6nFe5R/n7vEqjyfJRtV6ioQ2P/wkVSkNMsLDVKOCMI+9qpXyM6oHkcPJZFhF7hb4SapS3k7cbmtGVYlXVU6x1hHnVf7v9kSvaq2HhFYAC0xSFffr3TSqXntVj8t/aodtdt5xeZVLkiqcuSoHLgmLGKW+7aGN5N1h4lH77MJhVY+GVYkQuenF22FaUTtss/OOy6tCL1oa2GxVgxkORUKPpUpDKwIKIbZ/74eEpUDNqzATK1RyGUiZ2p8ZUwFZw6qf/nCn/y/fpmvGIwopXpUjrBoJ1sJf1rCKm4YV06K0IqAQYmP3TsJSoNYL9CiDKbOrrsJfvq23No3E1tFXqfJ31ernhPSqiOHqhCD1i9DSzOg4ZYi6y6gALgJ2Tkr/3YtCVe6x6ilD1D1DqXpxbMeD+h2I9KqI4eqE5LgGMDR2cs2kgBmlPqShVJK+SpUVV+dNlVdZi4D+ttRhWNVhwsy1vw5x9dOdjFj3lMw6DKtGMi6e0XB1267nqfIqay8wsApgDnqqWb2UKk/yYRYBAYxXYcIq1at6OrW6lXwxFcBhVa8JFbK+BDDI/cwUU0WvkyHHLAICGK/ChFWqV/XUFaykTPiJian6eKzuJ1wX3JtTbM4pl2mi3ujuefv3fvjGu96evgP7XrocFLruevWaEG++HhMRv+V/fOOrb/0SPH5n87XEEVHpA6rwo9T9lvmO5e+8tvSexJ3pKWWahmjroQ9dOAd9/NXjp2roy9WwCqNr6m5X8hZIkJ/LCCnWNKK7543dO3fcvpu+A7tevRY3bxM+05Lqtn7pZKKOlLSZPpqTB7Kkqmma5eVleAyz9KrTn5TEFVYJRF6FvIdgsYsBX5l/+OLn/uji5/5o/dJJ+MPPpxBHSkwFcO1Po1jTCJKM0F68tvmu1P2BP/x8ColbxFPbxFS1UaxpYLptV1glEHkVUguKjQY5eOT8rlev7Xr12o2Vw/JPIPYz7jq+xJjKT++UiyypMqffzXEDhNCY6sbKYe0jQeZV65dO4k8Uz4vNaOqV+Ydd63ln8zXXP8ml4DVyta7tBhXso6v7mMKfyjjDqjJNA0n5UGTl2Gmr5ViHW125+AnXeh479ifW503b89tV6NuPnkkhyKjGGVZV1TRUbp89oH2hIfOqoI7D3Ir6T9ozb3z/Idd6tj943fVPcmegy6Cd1TN6JoVBjk+XUJb/AGuTMKfEzTqfm7/zhrwqsRQIRcDdT12zGsn6pZOgQa/M/yhlKyZWu+r8fjL44WXj9Coga9PAZyFlioAeS1g5dvrKxU9cufh/km8RHmh25QquakuPxulVgEukSFoHvtvef+KyRxQgr0osBUIR0LWh22cPgDy98f2UjViw2lU+oUEWOpE74P9caoN+oDrc+Vw+hgcLg4g1B8VUniIg4CoF4udWcF0MCAt6YicSXpl/GP7kRs0dDp11PUjOrDdO5tqfh3xNowzpFcCVY6dBdFyxEwmHLpyDP22jEUQPUefCXxDWpiG6aB2eIiDgKgXi51ZwXQwIW/TETiQcPHIe/oQQamVQBeM60UPU4/SoR4FWlqv/YLChvI1U54R6FR7pVS6tieMjP30yYqkrFz/xr//l/7Z+6SThxbquUeqhhT+VTq4E1Ea/yv/NNCrWRaamEdpzFx4gBVrT7YzkL9/b8/K9PWXeb7RRjdnAMjWNlM441KvwSK8Cl+p2egU4RFa7ImfYhT+AUqrUu3JCaZxw5RE9MWYRq1dRTQT6K3c+nbJ4HJ5RKSYpMykAcfNKFPYqdTys2By3IZTbx+Yma9NI7ImDPIPQxj4n/svENUQQdKxSZlKIZmxelbVpxHXbmKWsXkU1Eei+f7kat2Dc1YVA0LFKmUkhYnPRi3RCL+epwtNaBBRpeZV1HLpnKHoZwKvyDbSyFv5ExbW/+Xw+nU7lY/X5SsLU8uTuvE3rUktyHZL7fsZc+BsArUVAkZZXaR3H/hOXazCG3Dcz7tG4KMAU/ZmCa6l+SFXuYMP0KszcClrDsOpUJ2GVUEqByNeXjKmAHJ/p29//Te3PfI2/PfQOkvvPxIVVHrQVdq5T2tbhf6udyJRtjIQCmmJ6FWaj6tdyJTolkaVAzItdvUZrYBb9lkseK+gmZrPZ0aNH1Sdb6xs9kKrE3hcTVgn33ZddXqI+70qnXvnoRwJ2NA3rF/GVi5/weFVcIV/GVBVOKP/Db7xP+9NeAL85mqZZW1sbgFp12/u6pER93q9TnVQAJf6r7eJGgHFMVQnpvS8mrBLuuy+3dhx+nYquAIZi3Qf/1XZxgRPtbA60WH+Ny3loZZVD0vuk6rWl99y7voNkVa1nA6YUqA5Ir6HY5we8Kq7Gb45SdxX+RMW1PxW4yAiqgVVdSxFH+sAm2hHrMAmnXLNfLLodtw6AV1GNXmejqgfaSlPr709MKVB+D7emUwXGrbf2COBVVIex8sKf9dc4dBbmFP8w9b9nbT2QKiHE7bMH7l3fEa1WQVmX9CrrcHWkTqn/FFcBRF4AiJm4AYZYaa0ocYh6ekxFMlvVDzZ2+f/kK7UBsDKy6mTSf0IIL2fDr8qqYkid0v4pX1gVVN3TXtzJEHUJyWeKbxpDJb0jD8pUpFdZh6tjdEr7p2JhlQtrKTBxiHpiTEUiZ61NA1oHyJN1ftoqyn+TTaLXABIQrVbIIiDgyqugYdSfTplcufiJlMgK8MRUoYx2/k+T9KaR2AeTeAPsQ+djp6y07hLsdmJkRRVTjXb+TyuJrYPQqzDpkSuvgjvxlSl7pVwAaAK7nRhZ9TTugjiqaRpQKHxNo5BUxU3dtvOLV9T/lad1SmolQrxKDatgkvQ4nco9XB35PS4jK1WtkD84/BNTBdX+VKPSPuIREtc0tE+cyqviwiqYmb1CnQpCRlYRekReQAQy3cSwR8S1DlVfCHtivFepYRXMn061D7TgJzQXWyeywsdUVBNTaZ9jgZuIzBWEUutoLW7UXv5zeZUIV6vQAe9qXtXHdMpK0CxWVhILf2xUVNB6VcpuJOpUt8PVVdK9MGUNbFSEEHpVqAqodcBEneq8AijBvxFXVNYvo0qhXqmSPa7Hq0SgWgUVAYXiVRH3nHnz3n8FLwBEonqV5+xUR6lTFf5Uo5IpIKtVKLLnJvSq6BHrERutOdNyHVsPJIU/l1GxWoUiO+8cXoUcQg5eFbHRajMtoexbXEwVjcuoqlWreqVKxe9VIrkgiCGxhBexeNzNaoIIHWJljamQtT/NqO499NOgTTNWcuRVxUIvOR48X1gVLTolkz/OqDKRqQ5YANnH5QirEid8j1g8Whb7lVEBlUqVmV60epXAqVVQWKW+5pWPfqSrmTxdkNyz+f/6nZZvcJKJqawZlYTDKjymJVB5VZB/WK/7q4diu5QeU7UaFTsWHrPzpvKqoLBKvqZfJocHMxw+fWKqVqOq07EqlSorGK8SQuAjK/8ZT94eMjlZ6Lf5Y8f+RFvW41UkhT/TqDimooXcq/xrqNCiqAh6awWMikmHPK/ye1W3d0cW6AsAQ11HO3T+raQf5z5mVECfpEqgIw1/ZIUZsQ6f6I+++X7t+SAx6tHY9tZE1xVTtdb+/BmVhMOqRIqNWyeZPcFctp7h6gDttKhW2KhoQU4BFd3fYyzk9tkDu5+6lhLqm1upZ7g6IN+d6VXmMySjxHpkVKJOqfL3r+q/+n8TeNTKXwRUjUr1Khh1Hho4VThWHTh04dyfP/D78FjTRxilnl74sxoVx1TR+D2GxKv8MiGNyvp8bYRqn/VduO5mkxhTBRkV+1Y6tF5l7XrAqLQnh1EBdL0L669x2sKf36gq9K0apaoVvFcJxEAr64yxZkZFQuejstTan4ZWBEwv/CEzKgmHVenQ5lXa4qZR1XwdXwrmDZithyIazqjIwfTftHVArevRjKrCW6MCiaPUdz91bcftu/J/95+4rK2QtvBXoTO1Ut2M6sieNcirhE2trEVAq1ElFgFTaL0AMLpXsy5obW+ebwdP7c81ZzrHVBr4poH8rNO9yrohV0aVsiGN2iqAGq4cK3E9SKMap3il32/AQ1z3b7U3a0aVuCEV8gog4dwNplelbCLCqGoTr6pnVPcT6lXCUCutCIjJqGQtL9qrOgyrPDEV8KNvvh+aR2LhTzMqTEwFjDCsytE0CCexBAPwGFXitoJuq5w1GPOX89RbBKYU/jijCgLZOvD9N4lMaEVAl1ElhlWdD3iXSNHxvCPoOBILf33PqIC6yn+hfWqEVwlHQdBvVNbnMXqUb6x6ynwK0BnIAVUq0quoUI2KY6poQvtv9fUktSpyobGukCSsyrrDhKPHQo2KDYwE8iIgeaXPaiS1DVdXgR1O7DhSjKoqA6tLqiKI8yqxqVZqERA5jirfwPParhb0f1m4an+hQ6k0gsT69hs7/X+hWx8SiV6FHztV4EK53ESET4kxVW5DQjaN2Ww2m83w9zWrkIhQJN2r5Br8X5LyXwsPVye8rbIWU6kDqlTUQxrxiZTMqFqbRmLH0Xup0giNTGVk1WpU0WGVMDyMqgIY9LUua3+tS8E7TS/8aUbFMVW3lJkUdNhE3MRGUtKokMxmM/jv0aNH5YM+elUEJHlVtaPR/WS6GQ6stnKjKkBFUhU9ngY5KaiHlGv9co+RKnCzGpMfffP9EYcxMaOSjHBklR+qyxHy3aePMKzqcLg65s3GFQGpMrxoG/vlX/uS9ieEaJoGvGo6nQohZrPZ0tISyX4WI8UPEr0KaVSEYVV6BTC6QofXx9C3SZXh1WNjFUlVCtHd8O2zB/BGJV8ZVAG0FvW6Gq4uuw3rgCoNj1eZtT+/UXFM1RXRXnX1+KncGVW+9UffcifHUjUMTv/Lf/ek9ieEaJpmPp9LtRoh0V4Ft0zORyV3Vg41HtjtaE+qR4xSqEWq0sOJuMFVJPNRdT77VCut1/2ZBB0ZqoxKwmGVhPBSPgBpAylD44MwLwDsJKzKNJqK3KgInWw+n8OAqqZpqNZZEtpL+US2kU/RVUKzI+tkuHrQ/gd9KLSFv0qcrBapIiHUq/76U/9t6CZcqhHhVSkqFnfpX2i3hywCthoVx1SdU+wmNrlXTr71lHnnW5etIaPyAAFV0zQ9lSoqQr0qOqbq1+zq5t66RqljlvW/phIfIqGKyT8JYwm8V0UYlYpZASycV+ENKSKmUjEPo1r7c03yySDxNw3CAlkBr+rvkPYce165UQkh5grwzGw2q0qwPK0jX4HMLwRxRlV+SLvrAsDQ4xax55gi4FCNSlQy+efdzzx29zOP3fftt9337belbwvpVW973zfe+3qMBkVUDF+Zf7jzOwCq3QZmQBXgf7OYST7jYqq7n3ksYqk+4m8aV4+fIrQfpFelbxG/BpfNpFQAy7idP6zKZ1RB06X2HU/ruLFy+MbK4Y3dO0mGN2mq4RGCfS8lyRw+rHLZT80TVkn86pbJqGhnWIymovIfCBCoVaJgtXoVPJtWPQMAACAASURBVHn1+Kk4rwLSwyrk6zu5AFDFegyjp01vZTxGhaRHXtWvsCrTNKFsVMWAryZQq0TBwnjVjZXD65dOxnlV5/MvBDkHvP3EfXYdQ/l4eEYlqpIqIcTtswfU/jtFsDxepf5vnFd58pvaBq1D7S+925AHDWp/SKMKjakgs4zawYFDGFn5vUr+b39n8iyM/17LVEa1cuw0G5UVrddIESy/V8n/jfaqXoMfUKXiOoaCzqjWL52sx6hEbVIFWHORCMHCj6967+ufTomsUoiTsGJ5wNve9w3zyUwZFetUK2W8imTl9U9YlRhTWd9gPqMiWc+AsX7DRwhWUB0wVK0g+OlwwirkgKr0mKpYRlWVTgHdS9XP/eHln/tD/ZP2OxDerkyvMtcsvweDvArCKutIqdrCKq3bwA+o0ozqtdMfFSFGFRRTmUa15+l188QYFYcunDO7/Kxe5Q9d8kEoDWV+b3RrVNYTY1TsP3HZ9AN/r4G3K6tXaZYg+/LckVU9t1XG061RWc+NknQvVYDVq1rPJ0x8ZXrVvauPaq+J8yoPplcRjlXHz6eQeN2fCnyJa58ISUZlLfnteXr9J4/+MH3lA8DqVSS6Y/UqkhArKKzy+EFoWNVJyRK+YXIYlbXkN3KdUrF6VWuvgYmvrF6lqUCcVwWFVR45wIRVhHcAjMa8EpDEqFwlP1qd0u6MqT7pWaoWqRI2rxIhnu4RLHPKBr9XIdXKf2UcPq/CvDJ6rHr0V7AaUx26cO7//ru/I7YeN79RIWMqq06xUWnki6y01cL/lpzUqitSan/qAdn5xSvv+vvf0lZFZVTaMxxQmUREVioewbJ2zx6vqnyUFb5MZq39xc1QtbGhH1UqozKfJA+otDtjAs0mHq/qWKo0kXKVAkMjUFOwpFdJLTC9SgUfWXU+V4KLlJhKMyrzBekZlSug0p4ZbQXQPOyZIivr50s+jftg0IzKfEG6USEDqqEe4VbMvjMustIwBUtdrbQBjxYgvarzywDzoRnVQ0++AI/hSGY1qvQ1a2h3xpRAcFWvVFlJjKw0PPVBzau0/gnjVSXDKknK/TQwA6pMo4KYStJqVK0xlXVMulwtPqZ65b63+/+Q6+kL1o8+0ygr9X+jN1F4uDqyaZDMpABGBTGVhMSozCcjdnVsTcParUb3Gv76oCoHWh+Pz6vKD1fHmEfKEPXcGVVQyc/zZs224LrduHlnzLW1NXi+e6nyTxttQutVEjOdSvcq4Q2rVFuy3lm5Zqxf5feuPnrrzMGU1bpKfvB4bIW/0KaRdfS6ua2UTWS6m16HyINgzagOXTiXOOKeS34aQa0jsRTowRSCRK/CW0sld1YOAoxKxlRA4mV6WUt+//u//++0P+G4M+Z0OgXN0gZaqXQ5o7q/skNVCmyltQ7oVyuSWzKT89ixP4n7LpYxlVxci6la8cRUrSW/sRmVcDSN1s+OvBTo2mJ6j47cK4+O5JhbIeJ9mUalxVQp8Jh0K2braO1ESUqBGPx1QIxaIcMqz87nmF09MaYyMyogxX6KlfxUtDtjyv/WPqaqFarI6vbZAx55Uv/J2gekXBW4JazCZVoe4m6ljMQ6K5VJXFjlD6g8jHZYlR+qyAq/SNawKh3k7sWpJ/5KxoiwKjqgYuuyQhVZwSIeeZL/ZO31/V5VZmRV0AWALslrHaVuLqjFVEBEWGUt+ZWZNEG7M6YcSgXPeG6RWbtUiWylQI0Ur/rRN9/vH65OMnMVXACI/Bq1vsw/oMo6OJ0qpsKMSR9hTJUOiVf5TyqSwVUFKFNtJI+pOKDKRKZSoEaKVxVj/dJJpIVEqJ5/KFUKiQFVV5XTHkiVSC4F+mMqCcaroiMr8Cr/sCqqWUMTp6dCfqcjwyrkVX5sVNGklAIjppIK9apDF85hFvGfdVQVwMSYyjqUygQZVnHJLzcppcDWmErS6lUutdr91DVMBdAvB1QVwPSB89KorDGVEGL/icv4sKqTkh8JnUlVRE2nTGTVitWrqEZWkXylRqzEHEolvDGVVVLNmApZ8ms1qlFVACM+vpRSYEr2WRKqHQhdj9WoPDEVcv1UY9I7/1xKEtGn5hu9HkTuyIrENlJiKsKMqsOSHwlZpGq2CfmaXV7lbySYmMp8padPsnpV64RVmCAKMoaUQccRMZXVqBJBBlRjI1/TEDnnXtfWnyms8pPeNFJAZlQYOKCykrVpxEVW+FkA/GEVYPUqZFjVyo2Vw/CXvioNz4Aqq1G5YiqJP6zqb0AluZ98ja1TYyUCXvWTT+lnz+2zB6yiffvsgXtXA9Z/7+qj2w79FTy+evyU68vuva9/Whui9KNvvl+Iz/pX/it3Po2fLFR2HkFfuK4XuwZUuYwKM5rq1pmD0pPUmMo/DZXGeAp/uZuGcOiL6zT2nN6t649YthX8EO+4phGxz7Ahzagwo6lWjp22zllFNQ3VwCjQNPafuGw6h6fXCF3/wSPnQQjWL5102di+ly7ffCTGe/D7I9+jKiKtUnJj5XBoTBWXUVk/BQm5Ufk3lwn6pEqOjfdM5JBezcGXAvExVegi5hCrTLOryx/o72y+VmBGIiSuoxR0ld94jErgmgbJVOZZJ7KK28NM3qDGV+RNw2pUGFxvlqehcoFpGulxRVApMGKySswi5hCrTJcByuxq/dLJTCGWMIyqNaYCTHnqe8lPhV6qoGHM53P1jjmTrfz01x9P3xCmFBhdO5fG0PpNrXqVf2TVKx/9CIl1eXbpk/c+FLQq1xwK+Iv+YLg6xFShJb8goxrAsCpM0/jWBz9Jsq3WUmB61BQqMf7XkyiRZyWh79e1KvxFf6pCZS35DUDLrE1DbG0dN7/8+M0vE3QcraXAlBFX4FWtw7E1r/IbD5UPudaTElORMICSn0rGgerqHXMWW3nrH7xIsgnrVYFCaRURMZUkzqusUOmUhKQHSh9KpR7e6GmoRoinafzSn/0x1VbyRVZxg6tcpxltyESbV0UPpVLfLAdUeLT7rKlNY9/HXtz3MZqOozWySrmnSoRXuZyGPGEKXZt1QJWr8IeMqeDIy+MzmIBKkkWq5KyjOVZuEjd6HUOQV4FafXff59XnyXVKEtQVmQOqPEYVOjfVrTMHI8akj6rwJyncNIQ7skrs1OO8SntxplHn1tVGxFSmUYXOTcVj0vGUbxr55l7He5VUK0138hXszNUGxVSEl/u5Sn6Jq829wlYylv9cow5z1HGs60yJqeJWokZW+XRKRe08Qmt/Itv3eyaj6nsFsLVpZPo46lmtqmK5L+LTap2hCxJe7qdSz2dRFa1NI1PXaF1t+q1/g1YCXiXNJp9OSdT1x23LNCpkTOWn1wGVhP7qv65wXRWYDlwPiPxx/97XP/3dfZ9H3u+FhIjUAXbPulTwFOrKxZKi+KQJr9z39pKb6y/wWcddNOdfbWjuVXJCBLljoe/XalShMZV2JWZh7+GmgQQ68hweA9cDIkev33wku0up3Fg5DG88IqZKz6jMi/JK6lTuptGPGdXxyDCDJKbSQPYHKTcKjCPod7/HqBJBGtU4C381EGcYmHVWe/ua0AsDrx4/lSOj6nuSNHhkp04SU6kg5xAvf0MbiMRCJ2uwGlViTDWMgEpSOqna8/T6z4m8nrjn6fVbZw5CfAJTxfzr/f9L4jq1PMZPVx2MKzBQB1T5jSo0ppLgAyo2Khdl+l0Zn0DTePneHsJ1ttKhe2HiNI9RRd/pj3UqnTKdLsQn4EBUaiUnr2ol4mbDVLgm61JHqVNlVBoD0ymgdFJ168xB+CuzORgl+l/f+O/lX/SqIPpq7RW6/clefuvgmsgPNMWofvKpwym1XTlWA8ZtyP/NOt9gEIWnC4em8fC2W/IvcYWVN43WHcg0zVWBd93VLPPFyDdduAtCxcF4VYdGBaQMz4+IqeRHWWD02BAm/0RSTK20q29SBKvVqyr/aiMv/OHTuxQSdUoI0TTN8vIyPD569CiIFKhVyauNkBTrI7WmkSJYrUXAypsGQFj4K/N+B69TGiX7yGJe1blRtUIbU5X5BDvRKaDjgerSq7IOcIbOw7xrhOpVyBKhZ9B6Jd9urjJHq1EF1f4idCo0piK85gAma4bHq6ur6vP1JFUa5IPKrbiahupVyBKhZ9B6JU1DeO/Y4zGqoNpfmWgq9yZqxnonlhwQlgJdg9brMSpXEdBvVEExVQHL6UqkVApJ1WQygQeLxcL6ArCr3GplvRsXIAWr1a6s46uq+ppTew4YUEWbUZlv/5+9sHzjn5/d/4V/71okyKiideqfPvHF1tc0TSOTqhpiqtamAacW1WdnPVH9TUMKVqtdue48GLiPeTG9inBwuvlmT1x8VgiCgWuu9Q8bf+uAHpRKrVz9MfIKvlZqNipA9SoYUEWYUVkP7+6nru1IX7V7/Z1QSKqgSbzzK19951e+Ktz3yCNRK1gJnAc7dmyZENbfeQCY+Ore1UevHv+ruPkPu8LfK+NjKs2o/tkLy0KIz//9s55F8Ebl0ak9T6//p39wUwjx+sc/4HrN//zCZ7RnTM0CqapBpwBoGicuPgtnpuumwoRqZS3VYZoGJr4Cr+pX0/AbFT6m0t4s4aUA5so1og84/MCAWTchuJW/OmoAWoecYtvacRKqlWsrVF6lrZN2hZnwGBU+pnLNOGqdsT0Uv04dPHL+zt4DrS8jpJvyH0atBEVwZaqVq95hxSNY0qvq7DbUju1t7/sGSWds1SnRZlQY/NEUnAZgVOmsra1pY6o893AtD0atBNF9l8XWPjioaXgES3pVnU1DbG0d3/vuL5GsUHuGyqhaj2HKmSAbwmQyWSwWs9lM3tW4nl8dKq1qJYjuu2xuhfaqwGqNSg2rbqwcTs+oPBO4JxpVqyTBhwVGVZIux1SBWonMwZVwqBWy85CYggVelbhv+YCeA2NUmJgKY1S3zhw0Pyx/TOWPpuTjdKOS5qTVEaoyKok8OWmDK9fNWzS1Cm0apmDVbFQAtI5Wo8LEVNaASmw1qogZepHRVArytjDqHfekWomtgxHrQWoTYXBlrsqlVuSRVW2AV7UaFSamymRUrdGUfFzeqETnA9UB8uDKejZsbOxM9CqJFKyVuOVLkX6LN0A1KqlTApFRuYwKE01JqDKqPkIeXFlfqTlQStNIn5qhGOQZlXrQojOqTNGUdbihLPmpItUjCgRX5tzfI/Eq8oxKHQUfZ1TIaErSiVGJSqQK8KuVwAVX/mkatMgqqN7RU1q9yh9TuXRKxBpVkE6JcRuVxK9WgmLElRZZjaFp3P3MYymjqVwBlbAZFeajyVrpsw43BKOSYa0cStUvu/LfZwYfXLleY66fthRYJzt23I0eTeXRKWEzKnnPHOQKTczPoiujElVJFVCgJmiqVVzn4enkhoHHqFxYK4AAstKnwkalklgTxNTjTLXipmHFFVD5F4n7aPLNqSFHqYvNSh84FiB6ZVcpNUHMEGarWsV5VbUDqkgINaqgtWm4jn+HRiUqlCpJdE0QmVuq1cCUekf9eL7NPTGVNCqrTlljqje+/9D2B6/DYxlThUZTEjYqF9E1QWT3rFYDh900PGGVJ6byGxW+8Jc1msKgDShUL/rrkU5pZK0JatXAMZQCTTwxlXpwzFmvXEZlfhwR0ZSkW6MSJaVKRlBxS+UIrtTIKrTeMeDf4n6dEujCX0Q0JfEY1QPPfd0zq0IfibOWTIPZATWy4qYhaQ2okEaVL5qq/xKBUOIEKMdgdnXNcrWhpcAexVT+CqCGX6fwxEVTEo9RmcPjWlEnHIFn1BDX9cOj9OSfD/2bfxu3BmRwFYGmVgP+Xa5hjanSjcr/WWD0F59R/YcF8Q0+yyObxsf+3f8UtwZkcBWBplbjaRrWmEoeSddxaDWqktHUAJqGUFrHvo+9GLcGTHAlwgXLVKuRRFbWmApjVP7CX0o0JUFmVMimoU04Ip9sneaw3OSfcUmVRutg9mikWmE6j979FkdeBghG5Rk+lTIlFTJNHFvVb7FYkMhK62D2aKRaDbJptA5XB9KNyk/uSl9PWSwWWWf1JFwtxqt6FFPhkUfVE1ClTKCAtFXyqp91whGxeXMzz1w89Y6p8uCpCSZeCAoDrcZw6ZMWU7XqlIg1qqDK7NiMihxPTTCx24a60hiahhZTteqUwBmVtTDHLlUSf00wZbXSq8SwrgrUKoBaTIXRKRFrVEGHMdGofu9XZ+aT6oQj8kk5X7Rn1uheSpUkR3DVGln17rc44AmrMEYVQehANzYqQnIEV62RVU+bhieswhhVHKxTHUIeXGEiq4HFVPmMKtRK0zOq/+H5mfbM7/3qzJNIVVH+y0o+tVo5dneQP8rVmGrbob/C6BQ+poq7aICNKgc51WqYQ6zUmAqMqvVtBhX+2KXqQZoQ1S2Z5QoHOcRKjanAqFoHpAcZVdwRy3etnznhCIxSl2OtXAv2Xqp+fOXdbz76KM3NtFU2NnauHLsrEFda9ZE3dYrUqIR3qiorrFOZUDt78r780IVzV49vKQUOqWkgdUqED6WiuskBUy1QDdRKgf2NqbQKIFKnRHhGFWqiuadO0DIq+F/Mvcbvy7I7+fnxlXfDX+4NbWzs3NjYOYwOA7oKiKmQAZUINKo3vv9Q0C6xUZHz8r098Jd7Q4cunDt04dwwmsbdzzwmNmOqfEbFjIT9Jy5DatVfl9KAmCqfUYWGhZ1PRuWhUFJFculfAYVysbGx80vTLbdk7jX/4rV/ggmoRNrlfq1EG9UDz32ddk86hKSI1mHvfujCOSH2DGZupKvHTyE/kWqNakhJGFVhrhPerAYOxaturBxOmYCKkGijKnM69SCpKhZKtUJeXizPr+34f//Fa/8E+eI6jYqRFAulWhlAR373M489vO1WVqMawFEaHrk72l57IbBjx92N3TvxRpUygUIrNWdUQKVjqiIUitZ4Hn/hX1mfv3qccCMd8LtPLOWo+oXCRhVNRHdO25f/7Le32/+h501DoCc4rcFimQrZ9apLO4jnxyrP7bMHkFI1cqMSxaTq+j/+h/DAP6N64TjKZU5MVtioVL78a78JD/wzqhfuyJ3mxDAFufnlx+FB9Izq5LjNiWGjEqKYVOV2qdeW3rOxYf8nNqcI8sVUbFQaHboUm1MEHFOVJLdLeeZTYHkKhY0K6Kz8RxtKHbpwTojs8jSY0bh+0o3KNasCGxUG8m67gDwN5ionP2xUA+PgkfO7ut6HYUBiVK5ZFXpkVKK8VGUq8Km6w6NB/fiHVXFG1RX5OmxpPMObkLAkbFTDQ/0xkKl1jOH3BmdUKoWkquRgKSlYbFehpBjV9geve/41k1H9h3u7c6y2JCW7araraBI/pvIh9wCaRmG4dcSRaFT+mwXlMKrcTaOQVL229B7zyXcsfyfrRmnjq8HX/uKMyu9SQIpRvf7xD1ifnzwZvcq6cJ1XWX8S0P5AH/xv8Wijyv2l4Vr/RHwy63aLYe1xc09SUCC+YjA3XkwxKtf6J5PoVWLpckoFq2mJPLLF8RUhGJGSII3KJU+jpZhs8Q90FeSsCh5yiNTgf9EF4eovc8gWtw6JdVaFoJgq6A7WSKMivCs2FTXOU5VVtnj0lRVMTBXkUoBpVCxPieSTLf6BbgUZU5F4D8tTClllK651DDvERRpVhPeYRlWhPLkoJFXmfUUiOleXbAHWu875PQAfXw3sy04bq+43qgiXEkLcOnNQCCHEwYhlR4V57sWdbJ6lrKe3fyv4H+jD7jb8RkX+STEa5ukXcb55+mOrb7X23xxfiTajinagAXyfFJKqv3nyCDz4O1+6BA9INKsV07SsijDm+MpqVHiRsl5I+Nu/cobqPonWeHLj4A9IVl4D3/rgm8NffunP/hgeUGmWH+RWxhxfuYwK/3FE6GwQ1vU/vO0W1fo756WvfxwePPKB5+ABiWb5MU3LpQijbR0uo8K7lCs+pEqkrPO/Z71KUVJIqqRLebDeJZfctFo1a1SjrzSj8rgU8uY2v/0rZ+L2BFnbHZJOAdKlPOTum11b0TYxqh/omlH5j3bu7wrk+oekU4B0KQ/Ws5HWtDCaNZ7WoamJX4NyX1VQw81zNGocU6VSINByadbgU3ppVKpLIeXJBKlT0QPjhqdTiRQItFybGEBE70calXZIK5Enk+HpVCK5Ay2PZg27dUg70VyqEnkyKalTQO1SZdKhZg2Mz373qTcffTdSpACrThFewsk6haRDzRoSK8dOy7eZyaIIV8s6haRDzRoMt88ekG8zk0VFy5NJeZ0C7utkq4TsWP/53Jv4W/lgDKLrfUjYqBjGBRtVNDcfGZrxMCpdGZXoY1IlYZ3qFtYphukK1qloWKeGTYc6BfRSqlinuoV1imG6gnUqGtapYdO5TgFZpGo2m83n8/l8Tr5m1qluYZ1KJF/TYAbPsHUqa9NgnRo2legUQC9V0Crm83nTNLQtJLdRsU55YJ1KJ1/TYIbNsHVKZG4abFQDpiqdAuilCn5w0K6Tdapb2KhIyNE0mMGT1ahms5kQomkasJnZbNY0DTxZkkxNg3Vq2FRoVKL8jOoSzHSgAq1T1nvUIGGj8lC5TiHvy1YzckZ1CWY6UGYMpOgUpmmAPM1ms8lkslgswGzArpqmid40IXJGdQlmOlDBOjV0UnQqd69BL1XWNolUKBUePtUtletUH7E2DVYoxqRMvQ8CKiHEdDqVT0q1KrAD6p5YN4pUKBXWqf5yY+Vw6zxVdaZTKlnKf9A2opsl61S3sE5lIr1pMIMnk06ZgajYVJnZbNaJSKmQNA3WqWFTXqfU4rj2vKdEnqX8l1Id5+FT3cJGlRUeU8V4yBdQmYHotz74STAqeU7KfqITu0psGmxUw6Yro5LFcfX55eXl0lIVB+tUt7BOMUxXdHJ9n/whLjavvOvjlB+sU8OmgE65clyzOA6Yz6hUIVWsU93COsUwXdHhdAmaPJW/6C8R1qlhUyydsua4anFcPo+5jKNjqeLhU93COsUwXTH42afywTo1bGoYja4VxyVN06ytrXmGVXUpVX0JqH73iSWS9VQF6xSTztXjp7reBXoOXTiXexOsUymwUXXF/hOXc2+iBp0CzOK4rIz753LrRqridOq1pfcgXxmnU6Y8bWzsFEJsbJx5cu23Vo6djlhntezYYTl34f2S8PrHP2A++cBzX6daP6NSQG60TYB5HLpw7uFttwbWNKwHk9a0rEfsxMVnCTcxSOJ06vbZA+R7onJjZcteSfPY/dS13JsujPZOAVrTKrAJJFpGpf6vf9BhaanC6BRenkwwOtWaPGlu8eTab0XvT7Wo3+DyK940LULNEjbTYs0Korw8aWhuMdTERX2bcEDMw0IeaJmmxZolwehUYXnS0Dr+jd2U35z1cPDIefl4/dJJYTss5A5UYBOElJ5R/edn35ZPpsiTialTQWU7lz0M0qg05He3+bVeXrNEHtOSJXDtXhz+GUcKIC88UQdLdi5PKi57GKpRaci3rx2xTjRLjMy05Izqux55WT7ZrTypuLr2oRqVhhQssCtJJ5qVYytxFJIqdZqHHGx/8PrvPhizoF8RxmBUKtb4SiO3ZokMgZYcWiiUCaOhQN75PJy5m4aIVTS/IozEqFTM+EqjgGaJkQVauVsH3p9U/J33SIxKxYyvNMpETZUEWlVMqYDknV/5qvnkj6+8WwSGXkgJGJtRaXjiKw3r8Kx3LH/H+uK4eNLUrMmTAYub8iTn4+ndZeQuPP13eiilMWyjwhwEV3ylETo8K05/zRY6Eb8ZsZ6hYr3zyb6XLguHAbjA9NCaUQ1sWJUqT62v8R9bl866DnKc/ppLTSYRqwmjRqmyypMVMCokQYHKsI0q6KctJr7CQytbGtY53FRqiKYSyRF+BK1z2EYVSmt8Fb02lUFeZUlO623jJGBUGIKijhFmVB5a46sgaGUrK11KFV6erGCMKq4yNWyjSgEfX4VCIlvWOdzU/1UvhYUKIH7lhcl9YX+BiQPGBjK+SlmzxjhlCy9PVjBGVcnonMGAjK8iqFC2CklVoj+Z+I0qZZSPy6hOXHx2YJeOp0AbX3lwyRZw56e7Qleo3osDRqx3q1aF5SZxc9aYakhNg+TjoI2vkBvSiGgaFZLoTyZ+o0pxqcHHVJjaX9BKyAVLxfNR5m4aNZb/WrEaFflwaQZJvvjKJHpOUfOusaL6sIoKKm/jwl8E+eIrJhSrUZGEUoM3qhzki6+6pX9S5cqoduy4m+5VYyj85btWKJ9d8fzsKVw9firdq9ioEuF6a7e4MqobK4cTvcpvVAMbq56DgdlVz6QqaGR6KGMwKlqee/QL1uffsfxfpa+cRaoeRmJU7D1DBT8yPZSRZFQktT/Jnb12yxzGULZCUnX9H/9DePDQv/m30SvxGNU7lr+TePkYG5ULlzllYmwuZZ38sypGYlRMhdz5g4fhwa5ff9n/Sg8eozp45HxKOjISo4rDZU6Dp5BUpbgUwBlVbgrLk8rYREqlWpdimM5JcSkgX0bFiBGbk4d+lP/YqKh44/sPPfdo1zuxyZhdqi+MJ6YqU/sjGd/GIMlqVOOJqfafuCzyy1P64LZKuK/rHWinEqNaOXb6je8/1PcvxO0PXt/+4PXcm/C/YOPgD+Av624wKnGXniGNauXY6UMXzvW9aVw9foov0BsSlRjVvpcu990VbqwcrnCOzWopJFWvfPQjcQvWYFQrx06DUR2+8vl8O1OGExefPXHx2QJqpSFFisSlrLdh7inVduQYo4Km0XedksAbqfYTQdL3/VeJvm6uBqO6ffbAYIqPB4+c3/Xqtb6rVZn9rzqpCjUq603oUpA69cb3H6JdcyfIKQ/gQQG14lCqpyCNSguoem1XcucPXTiXKbUaku7UTKjKBPW1GKMCnYIr5mDUUa/DKnnp38Ej5zOlVn3XNZVKperHV96dNaMSiJgKJluSOiVjql73HBIpWORe9cBzX2eX6jWtRjWwgMokX2Q14INWA/teutx5RnX77AHatmetQAAAIABJREFUCQiqAt5aDgfqtXeq1DhQPVSnIiZT8BuVnLtyGAEVYM75Cc+sHDsNXkXyZrc/eP31j7/5+IHnvp6+QiaFCC3wG5VsGgOTA/PtqF5F+GblJzKwA1gDoToVOplCq1FBsVI1qgFcHGcK4ptedemkIDUhKWp9t6vqpCrIqOLmpvIYlUunBjCayoW8cRuhWgEw8onVqhNypCxwnvhtAGpn5JvuCrUUSKtBOdY5ZoKMinzmblOnTPafuDykIhdM8QXviFaDcqyzJHVJFdKoUub5jDAqk971HP5b06g3xN3+4HXafE4OKme7KoDntMRczO+KqYYaUAFITUycDcH60XBwRQLSqDLN8+kyqkHGVNq/wiFNnA3B6pr9Da4KSdVkMoEHnllAMUaVadp09V51plIMOKaSyFKgyBBZATmCqzs/7f1UMbJpJM4Cmm75VEbVu58crchSYL54iXzNA2gaQmkdnllAMUaVGE25jEpenIgcRDWwsEpsvvFMkRVAvubcTaOQVC0Wi3d+5aueF7QalV+noKu+/y/2/cffcb7GalR+nRoG+Dsoa5GVCDwmsMiep9fhf2+dOWi+hmuCGovFIrEr9RuMXLn/ZVajGnZABeDfmhpZBS1obsv1WXBNUGOxWOx+6prnBa1G5dEp6Kd3vXrNr1xWo2rVqcHHVNorZWQlogRILuKSzh7VBMuV/1756EdcXuUxqtZoSnbP9//FPs/Loo3KFVMN7+c4oHqVSEutwK48aiWi7GpIk1QBcXWl1tMPv07TqNRzgPt4idrqgwRI+7D8Vxem1ASH96V0++wBl1d5jMrvSbJv3vWqz9hErFGNEPUuikECpFkULOVXK/zKPRvKR/djqlxGhYmmJNKoduy4u7HRHu6pfYYYbkYlQmIqbRFNraIPEQdX5CCjKSQ5jKovPzmi0yZVraKlk4OrFFxG1RpNSVSjso4KMo1KnYzUY1SemKovFcAIX5SlQPjflIFWlQRXs9lMCNE0TdM02pPwXysdS5VpVPhoSuLPqMTWmCpIp8YwmspFXGQlFcr1T+TB1XhIjKYwfTMHVEgiIivPC/IFV0PFNCpkNGVFTVlcIHWKiY6srOQLrlqR8jSZTBaLBTzZNM18Pp9tYl2wS6nSjCpCp0z2PL2uDauKNqoBEBFTaYtrR8yjVshJRD1qJTi4ckBV6TNjFTWmGptRJb5HLY1LT5VaZxzl4ArQjCpOp4IKf0FGNarRVK7F1Q+lVa1alcivVphNRCADqul0Kp+cz+fyX10LdiZVqlEFVfo0kEOpNDkQCKPCxFSVlzngHjvw+ON/9c8j1mCWAgXFtAtcE0RCW+nTkEalfb4k3XblTUMkT5EgbBpkXWfQceCaoB/VqPCVPg28UWk3H6TKqCovAq5fOokfdmZFKwUKRzUw6CDkqwn+zZNHzCfVUEp9fj6fz+dzj1T9ba6VD5meyYHqYFTp0ZTVqP7j7/wDeABGFaFTgF+qrh4/9aNvvv+9r39aVDw4FL55Lz/2We35CMEyD6PYeiS3P3jdU/vz4wquAO00AOWCk6r1DPmbJ48UOMOjkU1D6yMJB6G7AKPKoVMqlTeNdK8CzLeprjZlK6FnArx+Mpn8nS9d8i8om4bsNqALaZrGM16kJLJ1yIHqYFQplT7htgSZMIFRxelUa0x1Y+XwwSPn4WV1epWMhcwjGSFY5oelrjYlYWo9etpq4fWYpiGE+Jsnj6yurs5mMxlNAbKxgHJZl+1Aqn585d0p0ZTElVGBVBUwKiEESJWosvNQv21Nr5LgBcvqVWLzkKZIlcRjV/KsGLZUZY2mJGWMCqi8aeTzKkGqbkETZ+ClqmmatbU1OAmhnwCp8vwQL4YmVfteuhwdTan4pSrFqESbVIFRqS+rzas04/EcVbxgWT81taKXWLbzH0Mt38JLlVb1k8EVuJSnjRSSKvnYv7mgio+n8Hf0g/9Ha7LixyNV8LVVuVSZX+Ier5JgBMulVsI7Sj0If3AFDEaq5GP/fhIaz8PbbpkfItX6TYGov2lQeZVwv9nc69cISqrkb27ZbYjNUbr+MkdukK0jqEv2qMCdvQc2du/UdEqQGhWsTX1ZVVJlLc9hDi9GsFxCTDUQCnkk8VIV3XGUm/zT/4LQATT+oVTRARUQalSishEk1m9wzECo5x79Ajzw2JU5eh2gMirRNph9YLQ2DfIAKV9AZbWT+puGoPMq65slPMKtg9ldWEeNaMiRuTXEVIC/dYT2x/6+f9er19ZXtnT8hFf5SaPSnq9ncJXrYGK8Smqi5whjLrFMoXUwezE6nlIhbjBy6xwKGlRX+alGZVJJ5+H6Bj904dzV49gB5tKuhE2wXF5Fi38w++ApMxI5t1HJrVTfNE5l9SpakMViFfPXualZUPirpPbnIS7bCB0PFGpUnpjKZVRADV7lOqSwb/gZp9SDYB7w3F4lEIPZC9CZVEVf25XbqFwxld+oKgFxa9jgC/esgmW9KjATowquRMELu2iNimpVmUDeNZlwWwWOCe2G5Ph0OXCkNrsqdpeSAhlVPfiPaqhXSayCZV4VmIkOg6uiA9UFxXXyGKmCDjguoMJLlVr+017WCfjJikiiOxAs8Kqf+8NCX3Y/+dSbjWQwY6rk7hVzKcIhREHr7EXToD0yOQ61Z1vycdCYqmqRrSPdpZAx1fqlk3H2g4+pXK/sKlkJupkMidTCZwFeVcaS1WNbYEzVfXGLRZN44za8UQn0dJQqiUYlups/Jmi7EUfG5LlHv/Dco18AP/vJpw5L3WHiKOMc6i2BCdeJOf160TRoj4wsLBKuc4QkCgfeqAR1iIKPqTq5VXDoRknM787eA3f2HpBJUuelT3JKS1UKQUZFSETh79CFc4W7kKvHT/35A7+PeWW+HSvvVW/8+L/w/xXen8rJ0buHhjEnLj6bONF/KPh3na9sV96ruGkAQUZFS2jhb/+Jy4XVCi80WXessFe1No3E1pFFqvx3xomj1ahunTmYaFTWmCplKFXJwTGwLbxXkYRVJhxZ+cnRNDCYeQlJNx+9kmJeFfrzhtartMnWObLy0FXTWL90UjOqUMGyVvSih1IV86pQh8s6SmlIeRW9VMmpsQibB8aorM8nqkPcl+Ar8w/Lx7m9Su0z3va+bwi0VwmiIqB1JZm8qu+61to0MnW6WVcboSxAbq9Sm0bEEejjZ9FfWptGXKfbGlNlGjHt2lvP0Ct1V3N7VWIkRiJA5koylQLL61oWqZLTj1pfkDisysQ0qp1fvBI6bZIZU8nvqaCYKussNdpqtTXL//3zB36/Va3gxSRetfupa/I+EpK+C1AOWptGDqwnJNUE34nDlTKVAj1NA7k4PEg0lavHT1nfYN8FKAc5mkaEUUXYhulJshdPueIvUykwcbWEExYcPHLeenwGEFkRSNVcQWxOyzufz48ePSpfM9kKZiY6iSemspb8dn7xyrbrbw17DwZ+o3KNUhdbYyoJrVd5OkUIq4ACXqUua/WqkatVRNP41gc/SbV1a71Ju0NL9JojlnKdtIRe5WkacXc1JsHqVSNXK0zTEFtbx80vP37zy4+TbN0s+Ymt+VB0ghVnVC6TIPQqKksj8Sr/Ueq7VxFIVaOgPq/eOmexFcw1jYDfqLRndn7xys4vXpH/C2EVRhpcF/2FDqXyfFFSzSsYtJ4CXqW6VNbIqo9+FtE0funP/phk060BVfoJSXhVXbpXkV8akji4CmIq+b9ZI6s++hmmaYitrWPfx17c97EXkev3xFStARVePlzlvNCMyvP6dBPC6FSQx6R7lfp+rZEVVSmwEz/LMlDd2loicBmVK6CCB4QxFS3wNX3l4ieiOwDMLevVsErgSoEiyqtci3Bk5aG1aaSfe60BVeLmIgp/GE5cfBbaRVzrwDSNuGNL+G3AkZWH1qaB7yBdRtUaUGkvRm5O4t/D6Bmq9p+4fPDIeXCjCMfCHLdodYuwFusiQ4qsSk/+KcHMAmqVKuuYdDWg0qTq1pmD/rkutZgKM5TKNe2nrP195KdPWhe8evzUlYuf+OS9D1m3KBB91eXHPuu5NSGsLbRkKXcgaFLQ7Q9eN/1JYt6XVKTNEappGZxUrSv86a8/XvMMh9amkeIr1h7as8LQCRHijErdysPbbrle9vK9Pa4tYjb64hO/8fgL/wq5G0jwW9eW8mRv1vsQ0H7ok8nkrX/QEupU3jSErXWk3NzXKkmeFd5YOewPnDRJai38uaRK7rB/GLtn8BZJXS90znTVeKjuaR36GbViahmmaYi01lHvPFVIo9JKftaYCh/DRBtVCuqvc8zP1suPfda/Ni2sAjx5lfxCxx+o1lfy6PXyZBqT7ll/AWTrwGz9xSd+w7+q0LdANWhdhUevl4dkTLqH6MHpKWGMzK46SXTUo4ffAf8rs5YCi1GpVJlG5S/5efBfBqimPnGX+8llrUPUqwI/2wIVVKXAUdlYXEku+rJT5ObiMhtR9lYtWQk9UH6oSoGjsrHWztVMQYJKfuayrn9ScyOSy/06p4b7OgNUpcCu3k6NUmU1KvNlplGljKaq57vpj7f9qedfISLyh1Vi62WAKq4hVkFhFbzGU/tT4cgqNykBVah45dOj9Jtz79hxV7SFVSK8pUdcL4kcd8+RFS1Wo9KewQdUQTfFa8U1QSiVh/nnj8ikF6FhFf79uiKroN3ris6kCj9bVWvJr5U9T69bXUHGVOoXWVxMhXzZlYufCF25C49XtfZ86V4VCntVJgrMi5bSzXeiCB6vijs4+KUi3i97FQlIoyLZljQktZvPF1PBVpD320khIqyKKwIi6WkpsLqkSo2pgkp+6Rf9CYRRuQZUBdX+/FkUFa6wCoj2qqCYSiW6FMj6ZSWl5KctguzFsxpJyRJhhLXgB1dFTA8RXQpk/bKSUvLTFkFeAxhhVHg5QK7cY11Zp2hHelVcLBddCuzQvQpJlZzAzf8yzajMFwQFVCrmyCozpoq7wV+O7zWPdal9T0pYJdpmW8iUV3FkpWFtGq0nVe4x6eaGckxzQIu6h+RhlYrrfaW8Xy4FWjFbh7WzVH0i95h0M6aKzqjqHIOVqCM58qp+lQILSZWcwM3zGmlUroDKY1SYmMpqCSRfW4WHqKtvxD+4yh9WAZpXqV2OecSiYyoVwomsdty93/+Xsp9lwDQNFaqAylyt58no9SMXTB9QBcCwKsA/uColrPIsnjiLKeFEVgNoGgLXOqRRUQVUGtawKqhT98yGUAzkFuMOl7aU9b5+iQZJWApsbRqJraPL8p91WBVtQOUCYqrEoVQi5Ks59JsRXyJsHbTeSpBXCSFunz3wn3/xR/AXt0Wrlo08ssIQOg0V+baKLU5C66D1UDxeBf+7cuz0y/f2wF/cJqxaVsPBrBPVqMx/JS97gR6RDKVCLhjkDVTjriJkpdWr1i+dvLP3APzF7VVfJgitZUwVxFRxRoWJqcw1pxtV5UAHgAmrhHe2BelVqmDd/cxj8EDaVahgIUuBozUtsx/Nemtkz5PDmBBBJfEWNBL/oZN2FSpYyFLgaE3L2mFrz5CU/KwddrpRVegBOZBv03wg7SpUsJClwG6PcBVSdf9f7Iso+eGRawYtOHzl83FGpY1Sh5Xga3+El/5ppIdVYusQK60fteZV0qskEYJV5p42s9lsNpvB4/l83jSN/N9ekKnkJ7Z2zPJxulFlupVNBFnDKrH5Tq2plUqEYPE9bVrZ9eq1TCU/sbVjhk3c2XugzOV+LjBBVHpYBUcvPayyrsR8JkKwKr8qsHupAqMyn0fqVGtMpa4cPvUavpgw1T3Pa0zLsXpVUFgF+L1KcyDTqyR4wco9en0+nwshZrMZ3FZsNpvNN6HaRFbyjUnX1gy9ePnWsXLsNKF+qcOqAKtXpYRVVq/SHMizZrxg8eh1D2BU2pNUY9K1Hhp68bhu23o3G7yQUalbmbk9XV6lvgvPbuAFq+bR6x1L1esf/0DWQVTqyt+x/B0hxC/806+oL4gu/IXGVITkDgBcXgXqg/cqCcauskZW8/lci6ZArUhWno98AZVwGJV2Y77exVTFNqptCI4e3qskGLviyMrkxsrhfNNQWY1Ko84L96iIDquE7VOAY4X0KgnGruqMrLpPqjSCSn7+mEozqm1Xf6K9oNhQquhvwKAZrTxFwKCwSrTNthDhVYA/vkq8qFAI8frTv6z9aS+A2l/iVjokq1Fl2lYNeIqAtHYS4VWAP75KvKhwDPTCqEL7+zg/KDBNaByhXgX446usdgvjRrSf3+pgEit1SRXhVX5mRiWB5hdqVOT3Uc6B6VUpXeOfP/D7cnHNeKK9SmIVLGspEM8DZ/5S+xOboRQ0g6ZpZAvpl10RTkPlNyr43yEZFWB6VeJ7lItbD6AkQtqsgmUtBTKCdBoqv1GZiUscPUq50sMq6zGMXrlVsKylwHSgy5jNZkePHpVPwnjcpmk83UdFUhVqVJ6YymNUwP4Tl0PDG5WI2l++Ueoarrwq7v2CV1lFJ92rJJpgpUdWKrINQDOAViGHWPUCQsXBZFTmy0LXj99hqhmqMLjyqug363qb6V4l0QSLvUqDcNIEZEaFnGMdsN53OQVkCkU1XD1x8dacL+WYaIKV4lXWEoe8nmk6ncpXYq5wqmJGdaqr/ACrUZm1PxHrGYVxVQDxM54ndsk56oAeQK0IzwdIp2RYZU10O8TfNGjnSUcaVVcj1qkwx6q7SD+2/mMoITmYoFbDCxE9eFoH7TzpGKOSBHmVRocxVdxY9RxDlKi8SgJqFX0+WEsc0HFol4rL3+SeHqSKGdUjOmNXTNWaUUngA4jwqtCYyvV9mu8OgLRh1S/84f/q+VfaVGmE+JsGodkgjSp9EzX3+rRhlX+pHF41NjytI994ZL/3FJubqvMB14ngD2Od7xR+h2vDcKHE4f9NXlH5Lx28UQF4sQ0aUPWRnz6Jf3EmSGauEm1GBaheRRJWMeSEGlWmsEq7uhAoLGFUM1ddPX4KeRjVRUg2zdASOvtUaBEwmpRYi7ACGCE9mK1X7lVyGK4cNyKEWF5eBq/y1AFrkaqgztgaU3mMylr7kwTlN518LZqZVmg/JF8f9GbxLyb3KpYzCfnE30EZVb4bMYk8A6qim0bQ/uNfTO5VbGYSkp44ej7PVq8ivPVyv8B/LrReRWtmcwWxOd/hYrFQB5NYqUWqErFOdoUhrgjYyfRUJv5hVYlh1dve942g0h7nVXUSbVRx9cFKan/+YVWJYZV1tk8PnFdVS7RR9UKPrHFRkHlEhFUp906uMK+KoCKpQvbEZkylGRWm8BdHxLchLJJ+6V/cACzVq4LCqlCjAqi8ip1MI/FaPEmcJ2G23jtRUL0qKKwKNSqAagRb745zblK64fR7zmCKgHF76Fmq2mmoRJRRUelpPUJWkVSRYBqVv/YnwYdVlcRUSCLyKjgUt88e2PbQRuiyhfOq7a9v9//l3oEekdKvI/vySmIqJBF5lTwOK8dOR9wm2VxPPvBNo9e3xYyG6i5++MFVNYRbKddLtloLvCBitFnhwVWtTSOx4+iZVPljqoiMSv0IXV4Fo9SH9BvRY5DqP7229J5tD22EqlWiV3FMRYV6xkYYFX7AU1zTKDlDFR7PezHvmhx0g2RR3KswNE2zvLwMj3t3W0wg8cK6TBkVDKjqPD4pNmGVef/pICoftB5EXVIV2qEmGpVK9AwLKQQV9awvxsxWJcOq1iAB3r7ZiiLUSsKSREJoH5xoVCpZp63KF25hZquSYVXrbsi3b74yVK3MdWZ6PYb5fK7ObSi23nVgkFBlVMD+E5cxDhGxIdpR8yrR1uJaMPQu0VZSvKoqD6ti8k8kWkyFMSp/7U9+EnueXhdeJZdfZ/XU/oJ6I7MIaOqjalTyULy29B75ArxaRU8KOkIDI2kaGum9r4yR4DTzeJXHNggJekdBe2IWAc1tae/ROgYLqVZdDVr/wewXtT/zNXD1eFVGhWwdQcOok3fqbw0GvuhcXlVJT597DJZqVFKM4qacGEZeVcXknyrIbjU9o7J+Zom3r9Eg/8ZMny/U1d+4Miqx1asEWq14UlAk+KaBPJ3SB6dbC3ODv0GKq2kEWSNGreK8KvHL5Odn39b+zNfIKaTruYMTvnVEQFj4I5zbXSB8ovBYdfXdafvmyaiKeVVt+lVX+c+DGlMRVv3EZkwl0bzqva9/OjGmKnbXPz+esMpqVP6vCYxahQ6uGmFMRQ7ttOmaRpy4+Kw5rsj6ylaqGlDlCatcb9B/wWCrWtU5yYKcLboeqcKD6VlpC3/a70b1Nn939h6g3VYimSTMalS0N5yuTZgw1ChV/s4Vb1Se2p9W+NOg/dlRD9YrAT0ZlbUIqNKqVjx5FS3+3pfEqLTCn4bpVcPAeiUgUhldB8SvVkFelfWYy3pfH0ep4yGxHK3wt+P2Xeu/1kxc72aGVZ6MKrEIKEK8qkLrqlGqTGRMFT3Jp4r5MWz/3g+1Z979e29+5dXWhZgVQPydlYUQlx/7rNpDeIxKw+VVok2tkF7FypVIplv7mTeWkWuurWmY4O+sLIR48YnfUJtGq1EhwzmPWtWZV/UXT/9Ka1QuUu7rUpjEPcSPTCcRzfqPp0o/pAogn+TTGlNJpFcBobU/5FdkvtsqI/EbFf5njUetOK/KDZVR+WMqiWZa/ZqeKgj/W8PPGspe1SHklTjPgFG1CBi3xUwCQV4BzHHbac8aeuRVlUqV2vVCTBVqVK7an7/wp6F51TBQwyqMM7UWAVVcauX3KjYtPJ5r04DcRgVYb42cY0Ml0cIq2t1zRVatXsWmhcfseskn+Wy9BMf0qgg6H4ZlonYWyNtOwwOSIqBo+3DroVKp0qDKqKI/g+iZFHKMUpfhVvSXPgyuyjd0zKpWnFflIFPVD09tYgRE7xV+5irtZUj1saoV51WZKD9tukY+N4pInqjCqpLC18e8ql6pgn532/W35ri1HyamUnln87V3Nl9L324mgoZVSW6sHH7jXW/HvDIorJKYamX1KhasUMxr04AUo4pOj64eP1WzBAQNq1KJvtEnBlOtXF5V87GtE9nv5rj+Tv0G00apW1m/dLLC0evwfR4tKOuXTiLTuPSwSri9qlrBqleqgDijstb+XJ+BOUrdRT1elTgSS5UwpFdJgrxKGGrFeRUVOYxKA1/jG0zfHyFh0amYplacVxFCaFTpVlShV0WgdaChVc4cXlUnVc+ornW6KRmV+jHgYyq51MEj5+FzRXpVtV+I2x+8/ujJ35XnKDgrxqsSa4WqWrFXqcQ1jUxVP7wfqMkWLIU856uaoUplx467O3bclW9HhLTiiEUkqlqxV2nEtY5MRoWc0Fjd+v4Tl+GbE6kUuXUhrgJ4Y+XwjZXDshMUm/uJ8SqqjLBHXlXdjOouSKp+IrzwpyK9ijayKnMB4PYHr29/8LrrFA/yqtCwSiLVimq+9f98d5v/j2QrWUlvGolGRTJsXHpVHz0AdEp0Oj5MqhWVHw+gaYjk1kE4+iflK0t6Ve8iK9Ap4T6SQV6V+PapPs3WppHYOmov/wFBRmXW/uKs1nqdYGhklQ9NxTzDqkCnxNaTUgurRGAdMNqrxKZa8X1sSOjk7jFWDwuNrIoZjKeiJ3VKGO8FHpQJqySgVoO/I1AZ0vvgCAnQYir1MT6yQu55XOykLuWf2Ust1Jj7Jv81qA5YiVdlpQdSlZhRxRX+/CC9qsANavw9k6pTmNOx1asILxhE3puZ8ZDe+5LPbhDkVVnxvyNNp2gNL/HtY27MzPihNSrzFyBmlLpJUCmwKzSdwhzJVq/qhQxR0QOpSglFVOKMyrUUnG1WryrWo3jqhlKnhPuEtj7/xrve7ler9CIgQ0Xi+KREo/JMNX7owjlrK6hhQJWqU4I0MKtzdolxQigucZm66/cnRFbW3Ss5Tsi6e6pOicCO487eA361oioC1k8PpEqE9N9a7a/1NHVd+oc8v8GrOi8Fqqg6JXA/Ecyb/4ReEhgBOxkJ+TTFdekfcovgVTVEVhJNp4RXg+LKeSRFwKoOWn9J6bzjlsVbEXhVV3ph1g01nRK4jsN8vyQTn3rohZBRStVsNoMHcJ9z+b81QFX4c1GDV5k6lRK6eryKw6pQKmwaxaY1r0ERrDrFwVINVNg0VHIP/ezcEqw6ldJxeLxqJGEVmVQ1TbO8vAyPZ7MZ+Q3PI/rv6DQ16FY2khSvir4A0LNg0IhI652qMXlVnFeNysZyN42SNbU4FavBqzQwu9RJWFXhscpH7qYR13mnx1RBY087MQzPIPegjsPayWLyqoh33RcVI5Oq+Xw+nU7VZ6CRUK0fg1r7yzE+3cP+E5cPHjl/7Jn/RqpVgVHqJvAjQ57xNx9JLdK7vCplxPqojEoUaRqhXlX47nvmVsqnRBBNqdtNuX0hklBDGpVRiSJNI7QnxkxMFTdK3QoMsbqz94DcblBKRHjnmcSASsPlVWMYsX5/3GLqed80jfZP2jNCCHP2toipR15bek/ElYChRhUUU5ktFk6a9UsnHzv2J/LJ3ILlOVNvPnJ430s+ATp45Dy8i1tnDlrf8hvvert15Nn+E5fhWMV9LkOlk6YRRLRRBS1o9QNzOqsCduXZxMPbbiEvtbt6/BR+V13j9EdOaNMQXbQOSVDhLzSmMgMeeesYtU/JrSAp65cdx42Vw9a3fGfvAavzyQXXL53E70AnMRUUo5umQZ6uQKRUedao1sXly6haArL/Rhb+8DeoCUUrHkvBSrEr1dK2bOi3nOfl7bMHdj91rdWrWnF5lQTvVYOPqbpqGivHThPOb5Q1xVGnXSCxK5fB+NcJR8zvVdF6JBfE29jgPSy0aQii1oHvvDsvMEk7gS4sLr7SsL4pfy4VpDseXF5FvqEcwAk5m80mk4k8D+fz+dGjR5eWluilyg+0ENrqeCuy9kdS+Ata0PPTRAZX8L8Yu7LKk3W1/hfgvyAwL4Y6oKZWMqxikGRtGhivSi/8hS64iW9wAAAPqklEQVTomXlB2O4c7Fk/ybSiHCPVSdamgem8I+5IYxI0LsLzYu2exxi7IpxWFAmm44A6oKZWQd2TyC+7t7/498wnZUClVqibplldXfWfopRSJbeU9QoOfCgSYVRZFWH/icvKufWP4ESJy65CBR8ZVrkqgBIzsgoqAg4+pnJRpmm0kmJUuYfDSx18eNubm4vWvqBSnWgrAkZkThELjtbzKmkaKqFGlfuH5Rb32nvAf98YP60dR6juAK4KoMSMrPBFwALx4e7P/D/aM7e/+PeappnP57PZLPTMzJJUdUjJ4IRoxv1/pK0N+YuqNVnVSC8CAq2lQBejNaqS0BYBo0lXhBMXn33YuAFXPjXBD66Kw+9VozWqkvg7b/yXefoo9cRO6k192RwJLtfmd8HbZ998fWjHQUVX240GjCoiN+3H5J8aZvcMtb/Ewl/cTArk0M6fJtuS8F4M6J9bQUO7JJCnraoKV6REElMNbG4n9VhhhpGF2s/ADlffcZlTSuEveiYFcgivSRRoy/TPraChdW2Yaas6HOWmjVL3DKLS6KVUiZr67w4bEvK2ACrpkywArqkWXJ9LPZ8XU5I6rYJqWFjQsi4h45hqnFTbcXSy6c4vETCZKwilTt06RW1fpcokNKbSaliemGrvi39pLh53EiDP5k5O+qCwSmz1qq6+IBavb/P/dbJXnWOGVUFRk5bZDDWmMskRVnUFNw0r5vd24Ziq/it7QkdrBYVVYmvv1sloqtamkdg6eixVMvzYdvUnhaf6rAeMfqkVQEEXVomtd1/2FAE5piqP6lWFragv2mHF5VWZwqpeH6ueonbVJFf81QO+Aphp0nP81uUODO/eNT2WqhyMRMhcXhV3RYlZClQtio1qGERYRS+SLTPVyzpHl2pRbFS9w2UtoVF9t6OvyInrOEyxc/luvygkVZNNaFcLHXZ6TBUdyWLaBpwcoa2o9ZeE3GfMDSxN/HkVsgIoAa8a2DdFGTI1DbGpC4kxVdaZFKRVkF+x2OorOeYazbS5MZOpdciL+eUzETFV7kIerJ/8SxXTcaSvHAnsQ7Xzf8ZRSKoWm5CvWQ1CiuVMVokuNhCK6hQkrAMKw6vgc+GYqpV8TUOMaTgUQPU2rWFV4sq1IiDHVBjytY5OCn8FBlS5srSg91jScjSvMn23d4y0/KeOUi88k0IO92pdpzasykPocHUVrQ7IRtVHVJ+IcLLBuIK/CJj4NgdzlEZLPTMppNDacQTdlgMeRIijthu9NioxJKka3nCoIP2Kuys4bVglhHjjXW/v71fMUOkkpqo/G5N76Cpxml5FFVYx9VB+fDp/Q5rc2XtgMEXAgUhVtFGlxFTDaBumV6WEVcI9hRXTCdEd+dhKh6EMJm2Sk+7M5/PWOXgGRpBRycpa5TFV0CyghJWTlLCKdk+6ZSBSFdf9x1FtOOk/KT0VQPK8iqmHwfT9nZD1SsAaaJpmeXkZHsubcuS4pXGd4MdFpFP/DFUqw5vpoBgDkSqR5lXVlg67kveUJLak4DIYUrwqKKYapMBpXlX4ssHczOfz6XSq/i/kVfJ/O9mrkqR4VYUxVQQ5wqoI+iWdHoYjVSK8O+/qU8ynStY1Y050V1gVekjrNKrZbAa3HBejrHGI8O48cSaFvlQMW4dVAa68KuiodmtUd/7gYe3P+rKxJVVAkFfR3A45M7T3AYwg9CgNxqjEwKRK4Dp17QY1QcgsFNk24iapQkJ7IqYXAes0KtlPgEhBjUPe0Wk8YDr1ktWufJNUaesnJ0IZO8+odv36y9qf9gJoDrPZTBYBR0XuOmCcZBTrOLodzzQkoxLDkyqB7toLz6QQTWj45GkerV8cqleFDlev06jEZkyldRXj7DmQXXuvh6hn2ueUsKpzo8IAJT9oKRDrqnXAMeD/eoTsp/Ih6n5Cr3OMG1YVOlx9YEYlhLi/zGbklLiZJjnUuHXmYOWqVC03Hzm876XgL4tKjMpT1GiaRnYYZXeqhcJN4+rxU5m0Q6pDH1UMw8Pbbr18b0/oUvUblfx1UeHPjMKt4/bZA1lnWCjpYTtu393YvRP/+jt7D+x6tfTsEsMzKjGAGdVd+Lv5vsRUQGg2a74+bvwgMqyqxKiEo8bRbAKPZc9RiV2Vbxr+br7XMVUEyGFVgMyr1IPjOZ71G1XllG8dHnodU5UEGVYN0qhEsaSqE8jzqtABVR46+VmAJCis6taoMOMxZTqlFjhGOKZKJV9exVRC50OVe0qOsKov9lCyV+rwmORuGgMcU6VidvnqKPVexFRxuMIt5HhMObjKH3HVk1H5AYWSA9XHOaBKwwxR1GFDrFwurGGVFY6pekrrl2RfYqoIe9A6jujZqvwdR18sM46BS5WwdfxBn+jeF/8S/+LCcyUELRJaAdQuBjQPY1+MinFh9vrRMyn0ekBV6D6bg9bNI8lG1Ws0r9px+266B/RFxWgxj9uwjUqMQapEd91/3+ei9UyywEY1DLrq++t0DrxTZr3dMlMDJSdb1+hWO7JOrzB4oxIjkSqhSED0EPXoAVVZf6D4z1Fr8wj9sjCHq7NRDQlpAJ0MUc80SRWQW26sw9XZqAaD/KpMGaIerREkHYerAhixV9EVQLmtMRiVGI9UifpUIPGavjL39DbDqtoOI5POwDygjBRaw6qBHUmmw7wqB8gx+GpfQ9XRjMSoxKikSihCMOAh6hrpzePmI4dDJwJlekfKoKheD6gC4vb84W23kHMrMP1F9aro9GicA6pCJwIdBoWkarJJmc2R0/fRUYmk38GGcdH3pjEwQofql7yxzwjpe+uoQSbiZhDo9sY1vWawk396SImpuv3BEXeim0sNLNPuNVU1jf5GTcwgqap1jDNtUhl5uIBkXOU/4NaZg1zGCoIlbCRcPX4quoY1TiGLnoSC6R03Vg7HJU89tTH5azx6WFXfJSxuUsMxShVQg1cVi1hTmgcb1dgI8qrBDCSKkEI2qhFSQ0UvlE6m1x+AUcF/jx49GrTgkG9T0wp4VWs1MO7kyHRKmXcSuLFyGPNLSFsQczcG1qnRAqqU+9bLBYi4Ic/KsdOtEz2wTo0Z8KrWb904/SombZiOw+xu1i+dbP1l3juduvnlx80n5R1jp9Np0NrGm1RJ8JFVRIqbNfjNOqsCGxWT237yTVKVtRbJRsWIzPaTr+Mgv7OhRu+MSgix72Mvan9CiKZp5vM53Cs2aG0sVULUUQrEE100hAVbVez22QNsVAyAHGU1jAFVre9i5dhpNipGgvGqSgZURVcAkb2GEGL90sk+GpULuF0seFXQgqMu/6kgS4GtlLzRdw76pVOth/pmmf0YOq4i2mAGVGHol05x0ygDshQ4eHqkU5gO+ubmmCowqqCx6pxUbcGMrHKfK4WnA1E3Z/pTv4yKKUnKhYG9w/SnfhkVU5hh3zZY66TMPrFHRoVnrhC0ICdVOq7IqqrfIuR5GOsUgyFi3HffYZ1iMNQfWe24fXdj907adQ5SpxLhGdXt9GuUVRB39h7QCuRsVB3Su6ZhRlZD0iztvbBRdUvvWocWUNXsWEGYvYZgo3JQKKmqZErcIAbsVRLWqc7pY9MQIxhNxTpVA31sHUMq/LlgnfLAY6qyUOCcU+vcoc2Y7+vEdEV5Gxu8/zGDIbeQadcABm2Oew0k5cp/vd4EpLjWKcJal2pFPVlD30XoVFWhbyGCvn/W5en1EYN62bc++MnQBUMnqQp9CxFFyYh3EUSvP+hO6PURk9//mToOSehbiJiqijsOPDxQHQUY/TC+r8aQTjPFgBxoIvLqSDE41mKokN+03HGMCpaqMKq6xCPoGsAheSFTIVnvbJMVKVKDUUOmNurpOIKuARyYF5aBUqrgls5ic9Ys9cHAkKdaDY3ED/+8qIHxNA0pKL2wK86lOmc8TUPUpFatcMcRDdmYqqZplpeX4THM797aNqwFTvNJ5MuQm0jZqMnNLz9unnwvff3j1le2PuPaRNzabqwcRm5iAB9EzfSladBu4urxU+YQJeugpS//2m+2PuMCuQntSdg3jFEN4IOomYimIfJ/meT+ROCb2ew4kH1EdMeBWduNlcMpHUe/PohMkCVV6i1y1tbW4A7PA/7NoZIeXCXezs+6M0wljLlp1BZccTRVFWNuGqK7iof1PoDccVARmVT5Z3BfXV2F59XmMdlK3HYrBzR/nFtngNCmAW1h8E2j87vcdL4DTESvIbjjGPrWh0dkUoW5b3PTNGrj6eM0bin0onDOkBPaNBaLxWQyGVvrCJ1PgRkAEb2G4I6D6RtZrv6TPzVC70Q4JDT3n0xSRymqK5RrY/oFNw1Ana98In5T/m+cbKkR1ESgBksxtcFNA3B9z6d3HNxrlKHET+ShZrZMt2BO3cpDIG4aTCZaT/vKm4bg1sHkIXfHUXu7YhiGYRiG6QV87z+GYRiGYRgCWKoYhmEYhmEIYKliGIZhGIYhgKWKYRiGYRiGgBJSpd7UCXkjguhNwBS9mebkzbf/sHLMPC4p68+386LIR1DgRCrPMFpHr5uGyL//8ICbRhDDaBoi8/73ummIQbaORWam06ncynQ6LbOJ6XS6urpKu5XV1dWlpSXadZorz3SIsu78oshHUOBEKs8wWkevm8Yi8/5z04hjGE1jkfPs6nvTWAy0dWSXqoXyToQQ0+k0x4ckNwEPcpwKS0tL0+k0R6tbGPtPTtadBwp8BAVOpPIMoHX0umks8rcObhpxDKBpLHKeXQNoGoshtg76GdXV+XC1ZHJ1dRXyvcQIzrUJ9facJPzpn/6puiHY5/onzTMptvPkH4EVqhOpPINpHYNpGqLU/nPT8DOYpiEG1Dq444iAXqpaD415dyeqTagVWZJP6EMf+pB8LPcZ4kRa4JgUOKty7LwK+Ufg31bv7mgxmNYxvKYhMrcObhp+BtM0RKnWMZimIYbVOrLc+8+FOiQt0ybgs3HdCD195SLPRy7HG+Y7n3KvX91Qvo8AKHAilafXraPXTUOUah3cNOLoddMQ2c6uITUNMaDW0b9AkmEYhmEYpkJ4niqGYRiGYRgCWKoYhmEYhmEIYKliGIZhGIYhgKWKYRiGYRiGAJYqhmEYhmEYAliqGIZhGIZhCGCpYhiGYRiGIYClimEYhmEYhgCWKoZhGIZhGAJYqhiGYRiGYQhgqWIYhmEYhiGApYphGIZhGIYAliqGYRiGYRgCWKoYhmEYhmEIYKliGIZhGIYhgKWKYRiGYRiGAJYqhmEYhmEYAliqGIZhGIZhCGCpYhiGYRiGIYClimEYhmEYhgCWKoZhGIZhGAJYqhiGYRiGYQhgqWIYhmEYhiGApYphGIZhGIYAliqGYRiGYRgCWKoYhmEYhmEIYKliGIZhGIYh4P8HNlLAtPvPdP4AAAAASUVORK5CYII=\n",
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"gROOT->GetListOfCanvases()->Draw()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "ROOT C++",
"language": "c++",
"name": "root"
},
"language_info": {
"codemirror_mode": "text/x-c++src",
"file_extension": ".C",
"mimetype": " text/x-c++src",
"name": "c++"
}
},
"nbformat": 4,
"nbformat_minor": 5
}