{
"cells": [
{
"cell_type": "markdown",
"id": "fa480ac3",
"metadata": {},
"source": [
"# Background_synt256\n",
"Example to illustrate the background estimator (class TSpectrum2).\n",
"\n",
"\n",
"\n",
"\n",
"**Author:** Miroslav Morhac, Olivier Couet \n",
"This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Thursday, February 29, 2024 at 12:27 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "158858df",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-02-29T12:27:37.894497Z",
"iopub.status.busy": "2024-02-29T12:27:37.894155Z",
"iopub.status.idle": "2024-02-29T12:27:38.990402Z",
"shell.execute_reply": "2024-02-29T12:27:38.989242Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Info in : created default TCanvas with name c1\n"
]
}
],
"source": [
"Int_t i, j;\n",
"const Int_t nbinsx = 64;\n",
"const Int_t nbinsy = 64;\n",
"Double_t xmin = 0;\n",
"Double_t xmax = (Double_t)nbinsx;\n",
"Double_t ymin = 0;\n",
"Double_t ymax = (Double_t)nbinsy;\n",
"Double_t** source = new Double_t*[nbinsx];\n",
"for (i=0;iGetTutorialDir();\n",
"TString file = dir+\"/spectrum/TSpectrum2.root\";\n",
"TFile *f = new TFile(file.Data());\n",
"auto back = (TH2F*) f->Get(\"back3\");\n",
"gStyle->SetOptStat(0);\n",
"auto *s = new TSpectrum2();\n",
"for (i = 0; i < nbinsx; i++){\n",
" for (j = 0; j < nbinsy; j++){\n",
" source[i][j] = back->GetBinContent(i + 1,j + 1);\n",
" }\n",
"}\n",
"s->Background(source,nbinsx,nbinsy,8,8,\n",
" TSpectrum2::kBackIncreasingWindow,TSpectrum2::kBackSuccessiveFiltering);//kBackOneStepFiltering\n",
"for (i = 0; i < nbinsx; i++){\n",
" for (j = 0; j < nbinsy; j++)\n",
" back->SetBinContent(i + 1,j + 1, source[i][j]);\n",
"}\n",
"back->Draw(\"SURF2\");"
]
},
{
"cell_type": "markdown",
"id": "68d5cd21",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "82b58a87",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2024-02-29T12:27:38.994627Z",
"iopub.status.busy": "2024-02-29T12:27:38.994264Z",
"iopub.status.idle": "2024-02-29T12:27:39.797761Z",
"shell.execute_reply": "2024-02-29T12:27:39.787242Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO29XYxk13HnebKrq8CWeqq6KcmEjGk1DGhEwZI4WlBYjDjqzizJpKiVBrBFGxK8+2DMeOfBC8zClJ/UFDITIvUkixgDMw+LtbDAzgMFj7QGbA/J5pqZ2T2Q/KAGRNnWmhIBQSRgWaLp7iqQ7EKXsnMfoisYdb7uuefj3rg344dCIfN+nrxf8b8RceIMVquVEgRBEARBsHGi7QYIgiAIgsAXEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgODnZdgMEQRAaYjAYwIfVatVuSwShQwzkhhEEoZegLEDgcSdyQRBqIR4FQRD6gEsWWKGKQeSCIPgRoSAIQveoJQtciFwQhBBEKAiCwJ0sssCFyAVB8CNCQRAEdmjKoAH7LXJBEFyIUBAEoWWKOgxqIXJBEExEKAiC0Ch8ZIELkQuCQBGhIAhCQfjLAhciFwQBEKEgCEJOmk8vKIrIBUEQoSAIQjzddRjUQuSCsM6IUBAEIZQ1kQUuRC4I64kIBUEQ7Ky5LHCBckGOhrAmiFAQBEEpkQU1Wa1W4loQ1gQRCoKwpvQs6zCEvG4AiUQIa4IIBUFYC9bQYYCyAD+U+MkiF4TeI2E2QeghIguyLBnRBrUGh1pYN0QoCELnEVkQt65nSmLb1BqcAmF9EKEgCB1DZEGurRVF5ILQG0QoCAJ31jDrEChk0T2blWCEIJiIUBAEXqyhwwAolzrg2lHiMrX2qNbmPAr9Q4SCILSJyILmKxe1VStJ5ILQUUQoCEJziCzgUNCwVhskGCEI7d+0gtBj1jO9gJUsMIlrlQQjhLWF420sCB1lPR0GzGWBCZ928mmJIHiQy1QQIhFZ0NHfm9hycS0I60ZXb3VBaBiRBb35vQx/i8gFgTPsbhhB4IDIgh7/3ly/TvIchTWhz48DQQhnDbMO10QWmJQw8BKMEHrMej0gBAFYQ4fB2soCk04cAZELAh86cMMIQiIiC9puDi/KHRMJRgi9RB4iQt8QWdB2c7jTwFGSYITQJ+SxInQbkQVtN6d7dPS4iVwQ2qKTN4ywzqxh1iHQUfPGkCaPpAQjhB4gjx6BNWvoMADEbVCOVg6pBCOE7iLPIIERIgtEFjRAbw5yb36IwBy5zoTWEFkgD/pWaPewi2tB6BzynBIaQmSByAImMDkRGZshckEoCosbRugl65l1KLKAP309NSIXhEL084YRmmc9HQYiC7oFnCZuJ0uCEQJzeN0wQlcQWbAOv7fruM7XOpw+kQtCRvp/wwjpiCxYh9/bdcLPF9szK0UXBJ7wuk8EJqxhegFb4yFYSTlf5ircTroEIwRW8Lo9hFZYQ4eByIJukfd8redJF7kgRLOON8yaI7Kg7eYIoRQ6X57NcrtCJBghcIDXXSFkR2RB280RQmnsxIVsn9vFI8EIoUV43QxCIiIL2m6OEEqLJ04uFUDkghCI3DDdZg2zDgF51ncLVnquVhs4NJhSIhjB6gcKDJFLpEusocMAYGVmhEqYn6+4VnH7LVIBWmgMXpe+QBFZwO3RLFjp3PnqSjsbRuSC4EJuGC6ILJDHdyfowflKbDm3Hy55jkJpeF3x64PIAm5PW8FKL89Xrt/Sp2OiIXJBoPT2QufGemYd9tLM9Jg1OV99/XVSdEEoRD9vmNZZT4fBmpiZ3rC256v3HQckGCHkhdf13VFEFqzD7+06cr4QOQIRiFxYZ+SGqY3IgnX4vV1HzpeHcseE29GWYISQBV6XNU/WML1AzEy3kPNViwaOErcTIcEIIQVeVzMH1tBhIGamW8j5SkSOWxZELqwP637DiCxouzlCKHK+ctHkkeR21iQYIUTA6yIujciCtpsjhCInrhytHFJu5zF7MILVrxPy0uezK7Kg7eYIociJaxI5yCUQ10KP6eENg/qgfz/NhTz4uoXIgnZp97BzO+mS5yhUcqLtBgi1QSW0hpKoi3jOl5y4NQRPuunybIW8F+FqtVqtVoPBgMmvE7LAS9vG4Qox9EPbyttnt5DzxRk4KXJqrEieo+CihzeMdrl362IVM9Mt5HxxxnV2uJ2sfrenW09gwUr/Qw+cXWHilO4Wcr44E3F2mDwTJBghMKf/QgHgcLGKmekWcr44k+Xs9NtCp5PxsHB4AgvRrItQAJq8WMXMdAs5X5wpfXa4nWIm1jS7kBK50FF4xcayEBhgyxs5k1h1t5Dz1QkKnR3PZrldD9zakxHJXegQfbgKUworxV2sYma6hZyvTtDYaQrZvlwqVqRnxHqSJ/Qwn8+trqTJZDI4YjKZmKtYZymlRqMRzBqNRpV7Xx2nVstDXGHilO4Wcr46AfPTxKENFCa+eglGrCd55KFVFY5Go8ViMR6PlVLz+XyxWMxmMzD88/l8d3d3OBxOJpPJZAKLoVyArcGK0+nU3HJlY+J+lPbAkleKTiDegk7A6jTVagOHBq8JcqjZknRiwN7jV21Tg8FAM//D4XA+nytDWICkgK+wTU1S4NegnxT1FKAfzN8icEOeKZxhJQtM4lrF5LcwaQaStz3y+OVJ6jkGHWC+908mk+l06skYohqCqgEqGnBhVBghhGQqVV7c3O7GNYe54VlzOnd2utJOP/34FVZELnAjNUcBYgfD4VCbjnZ9Pp9PJhNq5uEz9RCg80AptVgstK0Nh8PFYhHRtpQgqITNWoR59HrNkbPD5LHA7WhnPCySuMCNk0W3jqc5ItUgZLN+Ep9fEoloBvN9dN0MD2fk7Jhw++FMXAs0zzFLe+QJzIdSBZfABzAcDkEbzmYzddyLkMLKS5ZdaPsSbZsLeR/ljJydWjB5JmTviZBI3ktFnsAcKCUUIHyAEYfRaDQej+MiCByQizUOMTyckbOTSL8tdDoSjOgNLZRwphkJAM1asGYkmDkQrSAXqx8xPJyRs1MObgeQyQMqu5CSJ3BblBIKZpRB67bgEgrmiovFIlfMIgtysQJieDgjZ6cVmDwT+u3qkCdw85QSCtD1EQ08FFyCGkpKKQhDgD6Yz+fT6RR9BtqK8MEs3dg663axiuHhjJwdJvTbQqcjwYiOkic91Sx+oIxyTFotBFjFOktbsbLaknmhNHx79C8vt3M949cKOTuJrOFx4/aTpUxTtyh+9ZhVEwLn+lf00NYt0d2LVQxPJ+j32Vl+/4Mb9/1NM/tq8khyO2vc2pORHv+01unhkW33cuEvF0QWdIJ1O019FQrt7pQ/4lroBC30eug33CJnEr3uBHKaeg+388jkAZU3q4Pb47c3iFAoAl6vDe9X7E0nkNO0zjCxYdyutIztEbmQHREKBWnsYhV7wxmRBYEsv//BtpvQBNx6RjBpBiJFFxhSdqyHZuB8HWCYWeWzCmb0WuwNHzxnR06TQGFyPVDhwqFJ2YsuKMldSKYPHoXVcdpujoUUbSvvo5yRs9MteL5UMGkVtytWii7woQ9CoSuEXKxieDgjZ6dbeM4XK7i1iok1zR6jEbkQjQiFpqEXqxgezsjZ6RYR54ubwWDSHm5ZFNmDESIX6iJCoTnogwyjJBI8Y4LIgnY5eH3n4PWd8OWznC+2FpFbe5ggwYi2EKFQipAHmVysbSGyoFuUPl/cTjq39jB5QEkwoi36IBQGx2mrDdqH8AeZXKylEVnQUZo/X9zuQSbt6berQ57AlfRBKDTf66GE4ZGLNRciC7pOi7mHbC0it/YwoUQwItcG+0QeoTCfz/3HdzAYaENFwyqAOYr0aDSCWRGDQmWnScMjcqEuIgs6inniEA4njkMbKNzaw+QBVSIYIY9fkzxCgY4KbWIaexhIejgczmaz4XA4nU6pVhgMBovFYjwej8fjxWLR8DnjYHhE21YisqBbcLit4uB2GzJpT49dHfK2ZpIkFNAr4FlmMpksFgttIgiL+Xw+Go3m8zloBdymUmo2m00mk8lkMpvNcGIJOD+/5GJFPKdJYAjn26oubC0it/YwQfIcS5AkFEajEbz3uxaYz+fT6RSMvQZdC9wJoAbgMzoh4IMZm4ijc8+v9bxYO3ea1pw1OV/cfgu39jB5RkmeYwlSQw/w3j8cDq1zIb6ghR5AENCJ8BmmLxYLbWvD4dD0SYTQm+dXvy/W3pymNaHE+Xrj2/86vWGNwe02ZNIebndriTxHJoe6eQoOCkXNf14Cz1bP7M2qFwWaBjJmUqdo5ny9dusMfHh7ri2WhNuFyq09A2aDS+VqTz+ewHGU6h45n88Xi4U16JDOykuJPfKhW9pWvAXdovXz9fL80w3sJRfc7kEm7el3FkW3nsC5KCUUIOiglJrP5+BUwA9COjwv1tbNjFALOV+JsLWI3NrDBAlGRFMw9LBYLGi3Sfg6m80wJIFpCjRrwZqR4MqBWHPadYVJEKFbyPkqB7cDyK09EozoOqU8CmY4YDgcrlYrTRzQzzDLLLqwWCw4lF1iSzPaVt4+u4Wcr1bg9n7JpD39dnWsg3ehuRLO1NhDJSUMSUynU/QZFO0e2WPyXqxiZrqFnC8msLWI3NrDBAlGBFIw9OBhMplAcUb4OhwOqYNhNpvt7u7i4a7MiOzliYkjzhUmTuluIeeLP9xOBLf29DUYUWJrHGj595g1FUJm+enfSYrDJRdMMyN0gh6fL9rZ4T2jvyi9uyaPJLezJu3JiPmO2tdHa99+j+rjSUpBZEFHWasT12Oh0O5OPXBrDxM8h8UlC2ptpKP0YZhpAbHGqiFs1rMLt2esc5LBC8/9Nn6+frjTYkuKwu08cmsPkwgyDUZorAzabWqTtJOjIOTFH6tenz48XUGSDDR6rA9MuKl2Ju0plC4QiFWmcDgsTBCh0D08ZsaDyIW2EFkgUNq1iCYc2kBppj2aMvDslMlpapc+CAUmPqtC5DUzIhdKI7JACITb9cDNImZsT3h6gQmrY9IWfRAK2onstG5oxsyIXMiFyIJuwc0WAkxaxdbVUbc9KbLAv1kOh6UV+iAUuku7ZkbkQl1EFnSLuCBd83BrVbfaU0gWeJqxhopBej00B8/Mdkjf7WtBsRR4ni/BRcT54nbNt9uew5+e06bwPD6td0ZYw9tfhEIpumVmRC5063wJWc4XfUfM2rpIWm+PphVab486LgvU0dt8k7LA37YW994kEnrIQz+c0usTjOjH+VofSp8vbid9bdsTGEdg4vxfn2BEH4RC87Ku32amf3Kh3+erxzR/vrg98fvdnuj0Am4WmkMbipIn9DCfz63WejQaoddIGwESVrHOoiuGjPVQ2hO1nk7p7gYj1vN8dZeXT5ymf8p24hqDg7Od0qf2DAzS0wu43dFMTlN28ggFHAeSMhgMFovFeDyezWbD4XA6naIggKEjh8OhOYuuCKNRN3zoxcxQUC603RAncr46iuei4nDiOLSB0sX2NJl1yOQZxe005SLJb0OHilbHj9FkMplOp7PZDF0Co9FosVhY3dp0FmwTV9S+Bv2kOs4o0yktWOEWiZDz1S08N9qf/eW/p1//zSf+j8Yak33hBijUHkxj3Hz3K3HtiY4j5IXD+eLQhrwkeRRGoxG895uzzEGiNUtP1wJ3AqwCn6m8wInpyNtnNO1GIjwnTmBIn260Pjn/M4J+Avxc1GEQTucusE6QmsyINn6xWNDpYPWtU1waYj6fg2thOBzSFYfDobbxQCSFLTtWh1B25MR1izU5X9x+S5PtCfEWcHuN5taeTtNQr4fJZAJpB1m2Fiile/m04kBeubAmZqY3yPniZoGytycuiEBdHRyOD7f2dJriQgHzGMbjca4Igv+sM/EQ9p44uSBmplvI+TLhZoGyiHX6LXGDHI4JhVt7ukhZoYChhMlkEp6NKHQIv1wQM9Mt5HzVgtsxCRQuVofB8vsfbKs9jcGtPR2ioFAAlWB2WKAZCTCFZi1YMxK0rAWBG5pcEDPTCUQW5IKJBbK6OlrsjMDW9cKkPR2i1FgPkJTg6dZIsx2pUDCXXywW4o1giDWzfcW76MI6w7YnAhRZMj93CD5WR+uJwKQzAp/jA3BrD3/yCCtaCIFOMbMXIU2BVlnA4ksoHQaDAX41t1yJqMUS1C04wa3uwrrRoQIh/2n2qDblf9v9WtE9ljsmzR9tTZdrew9vDw09bNz3N1naZsLtaizUHm4/M51SQsH1WonLwCrwmaoEZdRxqqy2xKTQR5/IZWZELjRDh2SBSZ+EQuldRD/rKtvTjFAIb0+n6d+vy5OjYFZNqDxMsIpZUwG+rlYr6ywrpoiuXEWglItVx/WMEPxIbgF/spyIjK9A3C4Mbu3pn2nPS8ujR3p0gOQlFKJ5M9Pdl11WiCzoIuGXfTOeUW63IZP2SJ6jnz4MMy144PP2iXmOch+G4DlxQodwWaC2AqbW9pToG5nSnhbh0AaGiFDoFXxkgRWJRLhgfuKELFBx0PpppQ04eH2Hznp7441RDA6IBhPhwoRS3SOFJuHW4c0P9qJc22wStj0VhVwMDGjvRFZnmdttyKQ91NXRbks4IEKhe3jMTIdYH7kgsqD3eGQB1QcANwvE7SKU9jCkD6EHJvdbIfrtlO5fMKLf5ys7Zt/IomTxJ+dKL5DrwQ835z+39jRJH4RCn7pHrqeZ6a5cWM/zlZEfr3aqF4olS0JoA1mH3CwQk/awzXNk0p4m6YNQ6C5iZij85YKcL85EnB3rE19TBg2cWW4WiEMbKNKe1hGh0BxiZkLgIxfkfHEmy9lZ2YYmafHkcruumAgXRNrTFpLMWApJYUuh+VRHOV+cyXV28IoKyTpsESbxU7Z5l9KehhGPQh7k7bME5bwLcr46QeLZcXkLPC+CTN4RJRjhR9rTMH3wKGgdk5rZo/ZBzEw50LsQvQU5X90irt9vrT6KLri9IzZzfb5260zgkkwOCyLtaYY8QmE+n1sP0Gg0gjvWHLgBVgFg7OnAFU1Kew7FzHAgPBIh56sTeE5l4GlKlwUeuF0qpS3Qa7fO4J9nMbZCilt7ekYeoUBHhUZgIOnxeDwejxeLBT2RMJD0cDiczWbD4XA6nVKtMBgMXCs2gJgZtlgTF+R8dQLPafLznj/9c1zRrwxKNFuJBXIg7VkrkoQCegWss8DYTyaTyWQCZxHVAAiL+Xw+Go3m8zloBVxRKTWbzWDF2WymbMNY50LMTOfQ5IKcL4ak31aoBl75jX9TwmEQDrd3VibNQKQ960BqpgzYfjDzdFOTyWQ6ndIpEEQAkz8YDEBDwCxwMMxms9FoBH4IuuJgMBgOh+FaYRCQqeRZRmCI68RpWkFohZTb6pHn/kgp9a2H/ndt+mq1Ql+CUurlX/9McjPfoh+3f8Zf8fL80/Tre0Z/kWWza0s/LjBKaq8HMPbgP6DTrXYdloFZNPkANQSohOFwSNcaDofaxgMxn1/y9tkJwk8cLmbOEgqR5bbSXvs+e/k/0q/ffPA/ZGhoMZiYAQ5toDA5LAi39nSaUr0eTNMekpYYiBmtpOBiIgs6Qbqb2pq7IGTBenbue+a/qJqywJVe8NnL/1FTCfzh9jxhctlzi9Fwa0+nKVhHQXMMZMR/o8plwZlybh7xLty6tbO1tZeyhfCzA1rh+w//L9aNaFP6ejqYvLNSi8iqPUzg1p4uUkooRMcLhD7RfPRH5EI40WfnxsFbPejWRxaYiIX2w+SwINza0yGaq8yIWQs0I4HOgq9WhVHOOSFkhE9SiMgFkxK5BQM5vEopfgeBiUVkK6SYtKdDlMpRwCRHnKJlKdJZVCiYqQyLxSJjfoOQi/TcgtKsc+5CrrNjzS049//8Gf7lbXbXYXKl1Q3PXz/coX8vPPfbhdrDBG7t4U8poQCmHTtAwgf8CpWUQB/M5/PpdIoaApZBZaBtR2gL/rLAxTrIhUKygGYdAnmb7WH7qav4+cbBTmP7TYRbAh23e5PJYUG4tYctBUMPs9lsd3cXzwSUSYDPk8kEaifAV61Mgrmif0dysrPDJ4iQC+3ndJpcZ4d5egHVB9tPXd3//IUWG1MXVkdSsbnyuTn/ubWHLXmEgrVqwmg0Wq1WZtUEuop1rn9FE+0Ei26oS/9kgYvV0eBSXfxdvZQFP1re/S82/kmbeN8z/0Wp8620JztMLBA3i8ihDRRu7eFG8WRGv6X3zJW8hEKsjyyw0pU8R89pClxdm+Jf/eX5p9uqx/ej5d2t7LcZxEL7YXJYEG7tYUJzvR6EVlhzWeCCm1xIP01HyiBmgGahAbidESYWUYRUJyiVzCg0yR88f2dIre6mHLZCW3mO6afJzDq89ff//Nbf//Pmsw6FujCJjXK7SLi1h8lpYoIIhZ4AWkFkQQSl5UIJWbAyKNFyKzhW042DnQ71R+ADt54RTJqBMGkPt9PULn0QCtoztO3mFMc0PEI6ueRCu7Lg8KfnUhofAugDlAi0H2NRGttRYzBR80yagUh7GNKHHIV+93qQJIMmqZu7kOXsaFesnFYX6Z0kmQTCNVpv1csnTmtT/mUr7ThO64dFg1t7mqQPHoU+IUkGHHB5F7KcnUqHQaYfUQTra30z7/qHb77t8M231VrFc74qOfX1a9qHcjA/6W3BzfnPrT1NIkKhTUQWMEeTCyVkgZxoKyAL6F/lKuF3011//F0w/64nfgPiwAo3C8SkPdzuEW7taYA+hB66Qj+CCI8890dKqW8++B/abkhmXGenVjDCfLB24sy+8NxvK7Xznsb3C2NVxxF9N1ERQN8R8TNdYOfJF7bIQ3LnyReUUnu/X8oxb21Pi3BoA4XJYUG4tacc4lEoRWPeAuwb2SQgF7pL3bNjDUb0zFuQfSigSm4cnHH9uVZBK57xbjr19Wt0I1t7J+EPp4A+aBJuVw4314K0p2H6IBS0h3VbbdA+NO8taEUxdIJcZ4fKha7LAgHYefIFqgk8D5Dm5YJiY4G4WURu9xq39mSnD0Kh+Yc1B1kguGjg7DCUBSde3YG/thvCCGuqAToMTMO/8+QL21/7Hv1atn0BiIX2w+SwINzak4uyQmEymeDrlzl2w2g0cs2az+e4YutjTDOXBT9ensW/rkcEImB+dvqEq3dD3c4IQAOPVOoq8FOpCayqAv5iWlYfbhczE4soQqoZCgqFyWQynU6Hw+FsNhuPx4vFggqC0Wi0WCzG4zHMoqcZRqCGFYfD4XQ6bVIriOHpBHJ2msfsiRDeZdFzWxUFDLnLnNPptUx+i84GbhaRW3uEEhQUCqAS5vP5aDSaTCYgCGDWfD4HlTCZTCaTCZxjVAO7u7uwzGg0ms/noBUKNVJkQS3a8li0ZWaEQDQnf1u3ldV+U63AIZqQSHYL/Z9mj2pTzPpLIe1hAhPh0jOay1EAd8J8Psf/1E8AkgK/jsdj/AyL0bnRiCxI5EfLu/EvpW+bHzlNXeTU16+1Kws0KWANFpTYb1tkPLYvnzhN/+I2wsRCyyOiBAWFArgQ0MCDn4DKBQ3wN8AsLUjhWsWP2JtOIKepN4QXMIjAn4SoBREasOWb+5v4v13EQlthclj6QdkcheFwuLu7i70W8UrCGARi5jN6MPuvW7tHir1hiMiCbmGer1qY4iBQLqAsMKsaUJp8ud/c36R/iodKUPxuHCYWmlsWRacpKxQWi8VwOISMRXVcDQyHw+gtm/3XpS87TzjIAhwWGT8IGuhz/rO//PcwJeP5osrg1NevUZOvDURJJ67evR+4/XJawVQGVt45/UGhBtSFiUXkZqHFKKRTsITzdDqFdEX4OhqNdnd3IXtxOByaTgWh6wwY16gWlWAycBegTTlf75z+4B/Hv6qOJAIoA1oLeefJF7QqyJpWCCG7PkhxD3iOZJNQC82qPUxgcli6SCmPQt1UA5xuLmZuSkA4FE6QIELrvPHtf33w+g78uZaJ6zkCxj4wXqAZb0/IQKWNOZmuEjRvQWIQgdsFz6093FwLTNrTIZquzAj23uzIAEEK/CpCIRDaDQH+Cu3Ic2txeyqZrKE7ISXoE5F1iJpg58kXbu9t1F29NHllgQduFohJe7hZaP6PLG6UEgpg13d3dzVXAf2PUQn4gF9pd4n5fA71GAq1s39k6bUodQu6TrgswLwED/4Sh9rc5bmDux990f/SXzQJsTFZYMLWInJrDxOYHBb+FMxRWK1Wg8EAekUCs9kMHQOz2Qw6RJizJpMJFGeEr1qJBRM52YkE5haUK5xQlIj4N3PM8xUIZg94cA2yrA24zA0mHRAQbhaRW3uYpAtwy+pgS9k7f7VaaR4FZDQa4VzrWA8qOOignWDRDX6iUw7NwX/f86d//vKvf6ZMMwWdxBRRmpwPauDmv73fXCx8fIR6u97OucE4WfDGuZtvf+VUznYEwM0Cme358cqipP/g+elXPz42p+eCm4Xm0AbOFH9F8Jv5lLlCCJx7IjTP9lNX9z9/oe1WBOE5cYG41IDWR7EBDrcPs2wn3W3QvFZgaxG5tYcJTA4LN/owzLRgIj0RukW5ghMgF9CdsPPkC6gSWkwjCASVQcOpBoHc3tuAv5CFud2A3NrDxBPM7bAwgW/QUUhBLnfOxLl5qIfg1NevWQMHiJZe4PH8V9b2j1AMy2jb++0AACAASURBVHMHlcvAsJObb3vTNevtZaofRjgVKqXA3Y+++E9fuzdwa/LOaoWbq4NJM5ggQqF7UDPzoaf/77abI1SQN/pjLYr8j+9PbuVxipY7hA+3dn6hjjQBN4r28ORmEbnB5JjIaaL0IfTgGuuhH3CogtxdUqr6xNHA+dIKGWkJBxxCBoDZTZFJ+OCNcze1KRhECA8lpCP3rx8mD3M5TaofQqFPYz3UMjNmHaEbB2fMvglCOerKAmstI0+BI2r1teGR6KzXbp3Bv7/+H74Y3v68qoKhJqDQIaYOtw+bVwYumFhEbnArAsGkGa3QB6HQXdLfPkEZ4F+JRgpI6ZRDRBs0OWXLzcBTGSAc1ICHrr/elIbJ8eEmXJpEhEJzdDqI0MUqyFrc4fDNt9WKiEecr4jix9ANgYYPQBlo+qCWXCg6YBLbgAKixRE4SwQNbhZI2mOF/+M6O5LMWIp1KGDAvOhhXVlQ+nxBTwRz7EScq01v0Zdg2n5uagDpkA6opPkEukee+yOlzsLnX9m4Dh9+vDzbVnv8cGgDhclhaQDxKOSh096CEG4c7Jh/qo1swUTgpb/o+UK/gsvS11UAubYTAkMPgUkXvQW1aOWh8ePlWfhj0h4P3FwLTNpTjj4IheZ7PfReFjRMw6NInPr6tVrnK3yQZYggmHEE68JZzHzcRthGDSrJKA5Cij20DjcLxKQ93Cx07x/7xYXCaDQC+23WY/bMms/naPhxVEkXpXs99EMWRITPO8pgMDj19Wv+35vraJh22swwsOYcWFcvFGvwpBT8zoU/oaX+f7za+Z0Lf1KiDXUxOyYkSoS8smC1v6H9Zdz4sR1xtYjc2sMEJoclO2VzFOCojcdjpdR0OqURndFotFgsrLNg6MjhcDiZTCaTyXQ6VWQQ6tL0OLegspxfF7Ger3AdUHlMKjeVNy6QRSvU9RNYhwUqjTnexK2dXzQ8CEUg5XRAjTYwe/hwaw+TdAEObShBQY8C+AlWqxXY+9lspsiwkKASYBYcXJQCMMD0fD4fjUbz+Xw4HIJWKEE/vAW1YOhaiOjY2dj50rokxFE0LbFW4YRWMB0DVt9AOZVQy53QmLcgDm7vrEzaw83V0TMKCoXFYjEcDvErjCsN6gHkAnUSDIdDHJBaHTkhAFiMzo2mT7LA01/RTDCkQwFpMNEN/u6XnhOXAvz28HQBXKzd/gja3/XDHe2vrbYh3HonelwszGWBCTeLyO3Jya09/aCsow9kAVj60WiEuQhWq79YLHAWzVpAbVF31OkeBxFUVNfErb2TO0++oNq3IxU0duK0rENrr0WNhiVCt5INOWgCF/YjubNcslcGLrg9xJg4/xFu7ek0pTwKGGIYDAbT6XQ6ne7u7qILATQBpZYIGHjBxfokCzJyRy4cp0W/QsNuHq0QMsVMLfS4EMp1TexilwTqOTixs2y3MQh/J8Hdj75496Mv0q8RG+HmWpD29I+yvR4Wi8VsNoPOCJBqgL4EGpWoy8pLnqb3HWrnQCV4tIKnWELd0f8CZUGKcMFah9Y/uljdzUY3yYR/fUPAml7As/Rh54IIVlkQoRVSHno/Wt5da3oI3B7C3NrTRcoKhfF4rIUbwKkwHA5Np4LQMGD5QjLIoPix9U95jXpj3oJcVhy8CIFZC3V3ylwTINwUQAir/Y3B9lIpBf+ZYPUWmFIgzpFgEvfq/KPl3eZfi+0pB7f2dIhSQqFuPgF6Gmi2I51Vd4MCUmnPrK74iHf6OFkQMYpESE2C1umEw6A3UH0w2F4Otpcb52612B4XoAly+RI0uDnbpT29oaxHwUxapOmNdK7WRUKEQvNA7kKKz7+ZpJDw+kVN0hVZ0EDI4PbeBqYpwAczayGkv+Lm/mZKukNRdwjVJVvX97au7/mXN/VBLi+CFW7OdmlP1ykoFMbj8WKxwARGrQeEIt0j4QN+hRUxHXI6naYkNAi1aLfiDbr9o5vRjFboiixQVcoAy/vTvyw7VQ6VAF8rjxvOqpVwoMUdQsIQ4Y4HTRDA163re4fnT6sjw09fVWtJgUK6gdurs7SnoxS0CpPJBMw8lkuCmkv4eXd3F8/TbDZDnwGsCGWXlFFiwUROdl52nnxh/9EPb3/te6p8F6Ptp66euUspb29P7LhoLXi89/v/sgFxwFkKmDDJLTD7QcCH1f6GOpoSfmAxCyE7G+duLV/Zcs3VvAXw9XD7tHXh1Wp196Mv/tPX7h0MBmd//+/ytjMC6mzn8BrNtj2Cn7Kvj+gVUEbsAOovucIKnhVNtJMtukEdr0y8ub95uH1oXcxlYkEl0MGTwnfn2pGrSgFKhO2nrg7UtnWVwCIH6XRRELx84rRS6j23X2+7OU5u722s9jfw3T06ewD9CrUUg393cAwH28vV/kZlBMEE3AnwXxHHACgGbWIIhTprcLOI3NrDRLiwpYnRI2mppfBZlXPXmVNfv+bql4idESjR7vGM9RUi3vvNoZWybFajExEEZeupSOe+fOI0/rXVQg9grXO5OsCaBoYkPDulszCIoI6bf5zo+ewh8L2FRljO/u5LIavEwe09ikl7JM/RTx+GmV5bXP0VKe+c/gA/x9nCRK3gMeQluiCGw1wc8KxVEA3+hFy/pdabd0j5B00WaDogRBZYPQfW91RsvCsDo5xW4GYR2bZHoHAcq00oCphGDEaYQ/ZFDOIHPSYqowNo9ZuJI1jhqQy6rgNaxB+P2Lq+58onqIupFTzqAUXD2d99abB9rJ1A6/WguFlEbu2RYARFhEJv8ecNaHIhha29k6AtGssksNKkAnjj3M23v3Iqbt0OaYK4vgZtYc12DAwQ5OXs776k1FuHrnVN4MEzcOt7/vTPX/71zzTZGICJheaWd9kuIhS6x2AwuOuPv1t3rTfO3VRKaeZtc3/z1s4vqAuBfo7oo2j1GRStecDTQ0DpkDIA6ho2WD5cLqRUR4CsQ+0DNiNLzwi/trh1dofugrMI6CjcrDK39rRCH4QCk+BWIaxDKUbnDVhfhakgaLeOgh8OmgD0locOyYLdTz3bdhMquL23gYIAQTvtr5qQIhpO7Cw951HEQTNwe5Xn1p4m4WsVwulT98jSIyybOgCmoA2OC0ZYK0B7Ah8h7gQOssBka++kVlOwQ8pAs3Cukv7nfvOvXvmv/yp645XmmZZuDKfJAAec09X+BngXbp3dUfuN7bwe209d3f/8hbZbUQRuzn9u7WmSPgiF7lJaFlDQloMU0Mzw4fahp9xCHJVqYPVu/enLUxyYbO2dPNw+3HjlrrYb4qTFt96QSIRVXVnVw4md5cZPXr911lmPy0UWYYGH0V+UqV16rBUAblaZW3saQIRCc2SRBYFBB5AFdz/6olLOtzcQB4H7rcWtnV9ULrO5v/nGzrHOnJtve7NEY1Kw2DOipZbnDnhqBQ6+8bqJC8rrm8HOC/6gQAnM2AcfTn39mlL1xnmvBbdXZ2lPW0gdhVKkj7AcnYigBQI8D1bqQsjoTsDhKCFT0hwZwSpQDt98W2UGQGm6UroA+9+v9jdoqSX4e/x//s9tN/AO6SaWw1kYbDetTuqSsTAahW2RA2lPw4hHIQ/NBBF2nnxhq84pw1fe23sb6rgOyBtlcGFNjYSEytY1AcDcAABs32jBYUCzE6xNjfAuUA7Pnz6xs1Q1yytDMqMrETIEbXVza/54xPLXfrrx/747ZEee0tHh76yVZdRT4PbeLO1pmD4IhebVXKAsQJkfcQPX1QQumnSPh2c5tKsSMoqD7IeXrSaoBGwwWNNAARFO4CmrLGqkGf7KDZodK2ppBQ+V40rc88i1n33z/qIW6D1/+ud1V+HmbJf2NENDoYfJZGKa89FoNBgMBoOBdVCowRE4/LSL1XHytfoO4UEElwOwkGMwkOW5g+W5g+rlkmnGSxFIh0og0zhC222xMNhe4pDKOFGzc411SaDnkR63iKOXcsArf+/y134K//FA4THUjmToHgu8Dt042HGN2uqaxc3ZzrY9PaMJjwIMNq1NHI1Gi8ViPB4rpabTKRViMMb0cDicTCaTyQTWrZQLucgVRPCIBnQw0H4BLdY0dOFKdWQlCFTw6yYHrcBTCpgMtpebP7kzHCUtgaxpBZgV8WqeCOTnLnN0SUgJi7iuqK3re+gxu/lbb576k7imHaOcBfKM8O6Bm0Xk1p6e0YRHYXd3V5syn89BJYAUgHOMUgCWn8/no9FoPp8Ph0NTZ+QiS8phpcMAQ/U7T74ApY5rNRKWLzc8ElKZb6gtFrGLw+3DLDWd+LsKzPddViohzjS2UhHZBM77YHu5ce5W9KDVlPCzM9hemvEI6nQBLeUZpy0RJq/OiLRnHSguFCCsAJ4DZD6fq+NOguFwCBMBujwsRudGky4LEH+UAQQB/oVvdufJF4oWPDZJGWQ5cC3wQFA/RIRW4BBHqIzg8NQELlx1/j31/5V7fEWEVk4sHZJwXQkp6iFuhAuPhNr/dHRbLIiz3Q+39vSDskIBPAez2cycbi68WCxwFs1agM8RQiGjLHABsgDMHooDl1EPN/a4ZF59ENhNMXqznmVMreChckTgttB+Jnw9sbPkLA62ru9p77swBebeODhj/sEsavlyORJKHCL0LqRvCk9idDsLeVze86d/biYecrOITIQLwq09naasUNjd3R2Px2auImgCirmMh4EXXKxclUPQBIoUDAhZRfuqjV5vHczeXLEW0ZogrldCrV24Dh0TTeAnl8ZKwW8ao5PmKIfnT1stX7g5NCMv1ohMXk9M3pqMEWSXC9bUwnZTpE1cro53Tn/QRnPYuV46TcFkxsFgAAmJ1rnD4TB6y37DX/SyKDeS8vLcQXb/QcrqcSMp160DbWoFJuUO7XYiYdjDvNC+iD9e2ZPRDs+f3vzJ65pjAFMU66Kta7WFWcx8YtGFjES35OYHIw9yONbu00z65kEb3jn9wT+Of7XttijFz/XSRUoJBdAHk8kEQgb0/2g0Gg6HplOhc1C73kCmoYbfJLf4vgu7Tu8ZcfODr5/6m+ZS53gGDqxgdwNV1R2f5hO81Ysh6/tuoeMWYaRpkYPbexvQFYL+b6Ylh+dPK1VcKFjZ/tr3Wtmvi7Z8CS6YCKkuUrZ7pNbfAb5aTxWmIGBGAgYjzKyFJkEFkN2RULe2gWn7S6uBOKcCQNuWa0TK7KSMRNwMVh3QSteD1f5G84Mo1jXSplbA1beu7y3vOxUtGkLGxsxyXs68+hL5dr9Savupq3U30rBFNJ0HmkRgYqHXefjHRErlKECnRwR6McxmM9oTkuYnLhYLGoygs9oVCiZFnQdgXwtlHUbQbglF1FKFfLn0oW9N0Q80UWaO/bnf/Ku6jdHyDdMzDKwEGjNWHTtrNYCeMkyDxcQXOKQb525FH9vKlniuVe3gn3n1JfijX7VV7r34Dfx8zyPX7nkkKC8hMDzv77p5+ObbTn39WoRXwFyFm1Xm1h7+NDooFBp7+IDpCxingK/j8XixWGC0YjqdpiQ0ZKSQRNDUQOtZctnrKdUVOqavpZZ3wZUiF2hvaO++8J2aoFawbselCUr0NSixteZJFCuH50/TPNnlfaeipViEckIRAPmhmibQvt5413utn2/+VvXwqtozqtIi0noP2h8uY9UKOFGb6xcW3PIKubWHLa2NHjmbzRaLBfRTmE6ns9kMZcRkMhkOh7u7u4PBAEo0+vtGuno95CWXSmDiKvBDi0Dk6qlY68eaWsEVqQl898XpIU95q3chpGs+LGO6E8wOip6NmBY9u43vUEIGJd23YR7JRLlgtscUtVt7J//hN95Pp7z6kQ+rIxHw2ic+RNWAItEH6k5AAp0KGuaDsa6rwFyeToHPldvk1hOBW3vY0txYD5q2HY1Gq9VqNptBPEKLLMznc5xbWUGh9FgPKXRCFihSusB0J5gWOkUxZDkI6S5x17ouL0JgAR+tQz/oA/gc0s/QdCfQbERtimezrumaf6X697Ak8dRr1+3yvsgsHNoebQrV2agbNDVgftWmAFowwhqeCOEdk7+FDykW0fQimFohEG4Pam7t4UZrHgVgNBp5kg/8cxnSFVkApDsJGqiGRE2CZiEOHkgK4deNLHh+pnWWv76hn1rm30pIiQLonW/9i255k0TLhTuamBxM9PTElXTElngylE2TDw3QWoKAbvjp7/2yOb3Su+kJFmDfRf8Wam05ZS1ur/Lc2sOEloVCD+iELChd67CZKkm3P/AG/ZDYLcK0oC6rM9heYk6clhyHU6gHm9Y3NEHHQJaAwr/9zFPpGzGBwQ9LAEcpY85momsET2UJ14IGlREQfaCcefWl1z7xIfhs9S7AdL9HAZ9CLtv8zukPEvMT83Z65Ob8F9eClSZGj+w0tCoAZylgpckqh7CvE2ElibQj6UmfRHHg+nrib99eo5UOPN3wXI4H7dgG1jKy5h9oZZH8q2jW6MfLs0qpX9m4XrlrjXJSAHEpAHM6Trl1Nmokw9jKSMeiOXtqgzRmed+pje/fjGtPBIfnT6tX1e29DRAQ7/ru9yAYcXg8f3HrvlOvqQ9tHnVSxecSvYMCDXnRDk0RTx5uFlq6UFJEKATROYnQFrf3NgK1AgUPb0SHi1xaQR23wbnqK1h7MYT3Ufz45/7b89/4n2h5ACsgF/gQ7SSAFVPkgko4d8v7Tt3e26AHeuv6Xi6toOUzAiALzGqk4Dlw+RUQvHGYPKCyvJkwsdBSdIHSB6HAxGfVDOG3YoTBzgK2MFoxhMgFzbWAHDywd9e3Q5/s1KKkJ/eB+T998Ud1V7x1dmeg7FLgR8u7z/3mX73yX/9VYttKUKLMg7blaMUQpxWsYQjQCim1HSM4PH/6hqpQCeG4ujKmUM5byc0qc2tPK/RBKGgnsk+6IePdeGJn2fB4S7WCEZS6RaBBNJz427ff/sAbWuKCqRs8b+cwq+FajV3pg1BOFlTuNEIupHgXDs+fVntKKbVxFIDYur63PHdHOmRUDHhrnNjfvG1kLSzvO/WOv/zr/WS5cM8j11Rs+gXcTeAmbH60Nm6v8tza0yR9EAp9IuPdGLGpulWlw5vRgFx4a48feANEQ62QBA6zVHd3teAvCJBWlIGVlHhE9magYjBxebmQzf3NpUOJWi3xy79TWyXcuve1rRffEbKkmQt8a+cX5iwte5dSYjQWOuoe7dLJwUKvczBCej20SdGeCJWYxtsa6dzc30w3bw3UXaCPaXgWg4MB8uqxL+XND77ejCZQzAohd528ZZECoX0X/d0iQlRC4E6jg4a37n0N/6ujDg7hvTm0wg95ufvRF0MW29zfNHtYrFYrVuNLrZtKUOJRaJjm3Xd1d219nJnu1pQYcJPeBaXUwcdOwqsSWJr8PpPyNJ+mcPjm21I7CxYgJbUwIhih5ZxinqM1hcIzCpp2T72lWowbofLWSBmnDXjrOZC7Untp3jn9wWuTD7Ay0uvjWhChUIoWNYEf08npeXK5iO6QlkK4XNBiKEcvSb8o0qwuYL6O7x8V8998mz6CgH+soHZpKwwBt8zh+dMn1BK7wroaE9cHgaqEEzvLiO6LeOK2SJsTMXtkUAqNBe/yH0gwoi36EHpoZqwHP+0GEWpxe6/2OEkuWnGnR9e2OvjYyYOPndzaO4nRB1VgXEpa3xBqFWhFnfPuzkPo0BLu0YDYklKmKSISRL33Ws8IGFwqfO9xlR+RwLjAG+duFn0ELc8dZLxxKjtlyHCUrdMHj0LRXg8ebzlnNeDnTgh/b6Oygz7Fle9d17sQHYCgpAQjQp62/sMSXnHBOnR1A+qKTyqiC9d4ieHQ31jaxwB+OPTGQc8IOIvL+05tfX/Pk+SYC/91W3pEeHQt+CWCNSsTbvZ3Tn/wj+NfhSlg+7WvEXB7lefWnlyU9ShMJhN80TdHbRiNRq5Z8/kcV8Thp1ukdAnkFgk08JVvQnW9C1mOZLR3wezfsdrfuPnB1yNy3wKHo6ycUhd4q/7R8m7r3z3/6/+XuP1coCCAAY3wL04cuIjzMcC5O7FzTLqFS9i4ItxwK1XeUNrIUqgSrM+iBh5HcKMtzx1QveLSLtRt4zmecfWkj+1IKkA3QkGhMJlMptPpcDiczWbj8RgGlca5o9FosViMx2Nz1nw+h9GlZ7PZcDicTqdNagVsCZMrrxlyucTjghFNDmCNLM8drPY3tvZOmm0+eGCv7ohTsJFf/I8/9ywTcZBh2IgbB2f+2Wf+LuP4CEW5+VtvgjjAoQ5dwxOgVqACImXXtY4MNdWa2fbYthM7S/yDKbf3NuqOE5ESgCjRhzkQ/MlUu0B7XAOTar80cZgJD3210EwoKBRAJczn89FohMNMg8mfz+egEiaTiTZLKbW7uwvLjEaj+XwOWqFQI01ZgBcczyvPOipglg54WbQCbCRxWL+4XUdXsa3s2JaF8MOLUuCffebvijYpF5osUEQBwAccyoh6EXB5OsIyfEhXDH7QelEzVvn6a5Kx+Clbh+Xm/mblz1ztb2iCYOPcLfgzFy7d0XGtXvAao2zowRpTwP/UTwCSAr+Ox2P8jNoivT11ZUFg399yROiAcgmGnodXXtd69LOy0rWw8yuvHNsRkQjWqo4qeTBrD4PtpeYkYOUq2P+0b65VFiibaICvqABcG6QbQelQSzEcfOwkvNnHHUbNsL1VNrHKTNaVCyGBg9Jug8Dt490U+EgJjKrU0gp19RPPF7yuUzCZ0XrCQDpYrf5iscBZVGHgKqbs8IN5JfiBv7cg+6ZCDDYs43/ewS3qr39cu5VVpAwxZc1zvHGwo969P/jptr6j4wNYI4VUwtb1vcPt/J3KyqFpApp2QKd7VvRvnMYg6F4qR1VWR1WWFQbRj766BowINGaFYOIqqOt+a+xwZTw+fc0rbIXmukfCSzy4B0ATUGqJgIEXXIyzLMgbNfDvKHDJom7P5oMR0XmOKmuJOv+vjsuGa4atvZNm4qH2ih+ekFjpSzAXQJVw5tWXDj7mOxf+FAFPLoJ1CiLFNJuh0JOHW55jp2mieyRkNSqlZrMZThwOh9Eb9Bt+npdFi0+ciIGOrHcsdHzaOHdr4/s3b53dwXSEwG0mBiNU4yUdc9FKcapEQCRRU21dzG/+Q9wJ1sUOPnbyrv/+C5gLImBzXy3v23zHX/41TAGfgQdNOsSVaQrJcESnl+uWsU7E67lJB8OxW7XA0LKVXoe2vCkMXxQ7R1mPAvRynE6n4/F4tVqh22A4HJpOBc5E9JdjVeQ/Sxvofe6JAQ+2lyV6A6o2vAumX8ETjDh4YM+sVoTL0+NPHQn4GT4U9THQzgVm7zuzJ16lpddCBuZ01/LWiZiasLm/ad11dI9KzF2IW51ihsw9l6Urvq4JiPSkyBM7S+0MLs8d4Gbh2isXPnBtmVt6Js93yE5QUChgL8fValXZvxGzFswkBjNroTHqphBy0AQUarOz1PnBRw/as+ZflKMfOqt370fu8ShxIcTSeKocuq6Qw/On41RCSLEEredhLVv7D7/xftcsMOouJXHzt/Sy0NqKyujgoAUgXvvEh6yZjHQKdRuEdFCs24lRkbETs1g73Jpnmco0Q9oNITDRkvbyGGwvyymGBpRBSoK5BCOiKRh6AJVgzVucTCa7u7s0P3GxWNBgRItC4Z5Hrr31xeGrZKUGXBTKLlRK3T5uz0J2lD0FMkt5x3BQImztnUzPR/cHg0AunL74o8rt3Dg4Q7+eueuGNosmG/qDCL69vOu9WpqCpwvDy7/zXvWmUkq92/AxYANQZLz2iQ+94y//mi722ic+RD9rIQazTFOE7a9Es3AlDN6d8o7bhxvGxjf3N2/bhpja3N883D7E6XXDFuHawoMrhMfBW1ALCUbUpZRQQEuv+RJGR8As2k8SlxyPx9PpFOsoQD2GQu20Yi0b1wlx0DxmwLVumAY/R4iGOLmAT9u6uQuYUqe5FvzB74MH9u76tr4AaAVUDJV5JDcOdipNoqYbMqI5D0LUxsu/8973/F/H3v5vfN/ugaDiwA9Kh5A0BQ2rpGjXwoEmMK9ebQperlr4DBa7cwd5K0TlaKy6vbdBd4/3zub+pmoq2SL7i4H0jAikbDLjYrHQchHm8zmIg9lstru7iy6g2WyGPgMQEFB2SRklFkwq/Uj3PHLtZ9+8v3KZM6++9OKVz9EXr5+frVirFTCIgObq8PxpzTxrVge+lvMxAIlaKiLpkrahaKqjlnWP2Xb7n1bbfxE5/LFZxo5y40DfoGllb6q3e7ZPuxdq63piCtbt3GnS8bpJWvcHV1Vmsw10Sb9KMK07TqnrS6AjNfAh3fK5RmDJiPWgRXcpqqTWabr70Rf/6Wv3ml8hQkFnWVmf4R8TKZWjMBqNVjZoLsJqtZrNZrPZjOY5AvP5HOdWllrSdlHZtnsvfkObAipBaXGHkrjKlmnQDANakGfr+p4pBZqs1aPF2nMp/ZQkj5RUx7g9KqW29k7uf1rtfzpDxYXKH/7T3/tlpdSrH/lw+DZN4513hAXcpt/HADt97RMfotkJ5dAGe6RZh0X3Swm/I1DmZs9wytLZh5W0glOppSng1yzpC4KVlkeP9GceZMlLsNr+ey9+48UrnwPF8OKVz6njqddFy8cCnvpF1N4fbp+Gr4H1ee4s3EjvfKpUPE+TiDeelP6ErhpNb5y7qY5SC83WmCHhcKCPwK2dX9Ray0q4TwXkwmlVnccARF/VrnQE083w2ic+ZB3A0AwTLO87paUmpNOiPYNhFc334FpvxtDxOM43EJGc6GoYHcCiAQeM9aDBB61knHZraI6EjIhrwUofhpn2cO/Fb5xR6sa73nvvxW/ceNd7f/bN++955NoZYxnl7bXVAH5nAE2J3/zJ69YpHkpEHCrHVk7vkRUtFyojEW+cu/n2V1KT4PwlgKIJ0Qq1nAoIWPQzr74UHnfY3N9c3re58f2bkHVIYwfHsg5VRZKBNUwQnprgohVxAL0SYMBlBWNE7W/ePvrsKUBSueXbe5Gh0zTEswAAIABJREFUt8rNxi1c6PCaXUNxR6v9jdvB2zGDDtqHOEQlWOmnUKAOA0CrJK+9IVmrwpVrHq3gO9heVlp6iukt0Dri1wJT6qxz6z626KZyhU6zyAXrm66L6FTHXKz2N5a/9tNcW6P5Af/wG++H1IoQ8DicefWl19SHsJ9CrSxC7T5ChVFXJTBxgIMsUEqF2zO6buWviA7hWa/Vzf3NFo+bp66U2SrzEZSlO3cK4lqg9FMoaCoBdcA9j7y1jJmpAOCL1413vZcmqcUlrOG6+Pnw/GlPHAGN/eZPXi8XQaB3oOdurDVghLYkfT0abC9P7CyXr2xFq4cUuXDzgzV0GAU6pMWtq2Ht+BBBin79h994P3gINo+Xk6C/MSRdQ3Mn1CJwRT6ywOOip0Gu1f6Gon0QlMKvdCAVusEs/RVxa5sOrZAxfOBvJ5an1PZI14KD49FY0Q6VEqP3SZ4jpQ9CwdrrATIPNDUQUmlOW6ZWhmBcOiFGEKzKIFouoHFNF+bp/RGyPK2ar4Xs6haxMkaW8g8MkSWDQaVln7hEQGAuJ+ZFhqsELMBcSTllAAkElYtFRA2oCbSOW21+9e+ospEamuYo1w1B1Wne7b0NeFzQBNL0R5D1rm9mdF9RCaofQkE7kRFVt8xKtIGvbnwGBVaOOEIu911igmGWNtCWRDcmgizBCNdghrXQFBt0oTxzl/0i/Onv/fK7//Pf31kywRWBysBv9bX4ztbeyeV9uunyj4xQiEqtYL7uU2g0rUTzImDicdFAiWCt4N6DOjTr7Frog1BoElbKQBk6oPTdmGKk4SXDPIDL+05FxyMaLiCdMhzlwccU5AekxLCUTbGZFRcor37kw+/67vcgzyB6pyGYPhVaSRBtW1EjZ2qCzZopBQ3U+nQNLpWyelFC9ojduKL3Yl2X1VBq6xyMEKFwB7O7l0lRleD3KnPT43mNdLTtrPRzuBIUXL0oG2D/02r7LyoGMzQHl7oz/SjXIfD433jXe0/sLEErxLb3rU1FrNWw50CbQvMArH0XtbKGmGpgUsKdEHdwaMShme6L/uBLYHgx7tWi1vKNPSTXTSUoEQoqqjxtA9w6uzNQrP110cGIw/OnN3/yujaiz9b349+z41pippXVIjrVcWvvpFJ3XAtKqVtnd+JSHUOujcPzp0+opVLq1Y98GD5E44o7wGGsHMqoNNgfwVNOWAsinNhZLhu8v7Ib9cANpugJ/4oZfxFbTeBnfVwLZYeZbh1IaTQ/U8xR9Yr2jfTUX9OGoGRyM/iJKyd3eP60OcQcLT1ZtCWakzz6eRc9ejWtwVDUTYXR9yxbo+GDoiMEhuP6aebFgOMo1vUNRKwCNDCUIh9K/0D6bOTzYFyf4Sj74FHIcpIadiTwudazEP5Ob0ZY4BFDD0dKFF/zybs6I2jB7BTvQi3XApUIOGCEUurggb3tv1D7n1Z3fXsnpBp0AxX+kZQiPIGdDuLwqIRCe3TRlg4IvFw1p4L/+LSbFtDRB2Pv/Qp8PQrz+XxwhDYEpUbdsR5crgXEWrCWSUiCM+liH2VEinch5bWj8uXv1r2vmRMD/QquYo5beydR0LDKlm34VZimFJj1+8zivq6FwQdgugGsU+jEQM+B5irg7C3QWoiC3tU3AT83abAZugqi6atrgalHAUaPHA6HMJLkdDpVxojVgWjFlypVAkBH3us3HssUnWAYnbhgTk/pVZgyvpT1Xc2qElSmGo4HHzu5tacgfSGEjXO3ir76hds/f4JCLaeCVpJIawOe0MqxlXPBQQRkLL1gHVQ2caj3WvRACvjpq2uBaS7GaDRaLBbYNu2rn5AEE1dZRj8RuQsoOHBwnZ+/r82hq+u+sKZ05Kv70MGSU1bFENeL8vYH3nDNCrRe9HEMpQIO/t1H7vrj7yqlrANGBIqGu/77L8DBgB8Qf3UmulOXcEG2XnwHfCjdKdG/gPVQa4MmqKNDDRU87yzzypYiL/raBVCosAEHcUCpm41Ir9jAFc3QYcZuTVCVNcum6nL9/2xuvB5K/5IcmYYeFovFcDjEr+BLqBxvuhKMZfzw6udhSqCDoRZWP0SLwQt040e4tRNDABFrKSOPAfPtU1Ido6nrW4YMx8o8RxQHhQaXUkodvvm20r7x5bmDSpWwub9pxhHoFjS3sycWgDGF6ARDK3ziCGYYJQLPbwnfeIpKoLnYq/2NQoNdCU3CMfQAgoCOMQ2f5/N5rYGnzXCRVeVZiz1Xgk4C5chpoHPRnVB6rKlyG4+jVoEmTR/QXpS0bx72KszZ0CqswYjKUShLjy9lllvYfNubrlnNY5VKUKuAikh/YmZ2twFbh4HVitcaULFyR65Z/aicKJSDo1AAPJqgMmEEFqjl/NFSGazQokxo+/2reMRELprRB3QvcekCcYkLak8pdErfdwqHLkyXC3VHAag1BKWGa8yICA63D99q9jljboP64K3SBQGs9jc8gqBcdWRussDEWrPSvF+ar1zedfCR1b9AQPPwFQomGHqoPOuJqafUx+AaXAoUgDXhUfMl5JIIDB0GtUjxLuAzVHvqe+SCJzuBNsnamPCocKVTgRLoYPCPL9UAVEJ5BguuhRlEyPWibMJBHHgOl/buHmj401/6a1VOrLUvz0+gQ8hmHKZOw/VsvOeRaz/75v2iEtLhKBRcvoTwuANcGXF+Bc9XilmmybqAqRJCXoWZa4LE8Y1S8qQgHnF4/vSJnSV4Fw4+dlKpN0787dvpYpUq4Y7NPndr+cqW1pi2Sju3gtWtsrm/Sa24Vgg55fhYnQp0SIhwOKiBEKxDtQEhJpNenNl9CVTNZDTernaavzdFNEQ8JMW1EA1HoQDUzUgwqSsXPLKAzkJPw70Xv2H1GdCvNELhT1BgLg6sRPv/03tk3d7b2LiTu3AILYHAhNYYrIZEnf+aWbK+/FGjmL2ofmIkYnnu4NbOL079TfyQ0+pIImiaAMnYJU8DXAtxv52/OMCXZm14ZY9WaLJ5FDiVWfSBVfr7B3mq9cNzPRvXeVSnRFgLBfwMvR7q1lGAqyGXd8GFphKobqA5jNpiXZQFLlrJLlSYwUBMDo4ZQaWA6zMFI+hNPrhTUh239k6mDK9wJ72AyKAGRhjyYBVwXZEFcbSeaoDjPYYnlnpIGYbeVFGlH4+iEurCtHvkeDxeLBaoFabTKe0tGYKmGUExQN/I9OZVhiRcX3tM4hgN4c9cLX0BOiKa3SIC0WLn9I25sRhE9OjVuWjlV5vw6aboJ/xajVYDrcsIFTXIRWCNRWs8Zev6HtSLbPIlqq+FFLPD1KMAzoPd3V2cUquIgsuzFO1diANVAn+5kKUPZ6JrAV9KKvsjFHr99TwWT+wst/ZOaqWQIOswLsSukbFbhAf/cWvMqWD+zM2qgR9bxGPz0MiFLKPWrBdirVSkwfbycPvOC4Cn6loEkM/omiuuhUCYehSUUpPJZLVazWaz2WxW63RWxp8SvQshI1IiaHQZFoSmNaCgeTQBM7q1tWo0WbMOoURP5autucDtvY3D86fNcadygYMy0NEZVEvF//0FHAEIT2hDknpsVVF3Ah4cWpMqegTOQmgvxNmdB65xFsLRxlE0/yK2ie2B4ki5equmqKKMd/E9j1yrXEZcC36YehSQiHzGkM6TKbkLSEgSA6Yxlq62FNgSVUeyJLbZ5WC4/YG3Oils7m9qj0woyAOf2w2cU27vbag6r/ueEjq5oErFKhpggRM7y+X+hqpfqyD74NRFgdbW2peWUWvtvmtN1Hdt0O8wiBYHdc1t3Uxhc5nbextafczoGsxdKf8geY5+uAuF7Ji5CyprMMJTuKkVreDJtbQuoC0Z0l/DjykXNHEQ9xiymrHo/ntWI5pu4VJGr65FrdIL7dY6bCAHwq8v6QiKucZDymgF24pNaEcju9JNUQyuEePC8QcgKKISrPANPRTClbuQJdXR5WMoPW413Rp8dm3fFCtWEaCtjsGIxHgEfg0ZAhhAhyr2VKRzzQdZ3KPNjB1olialiwFuX5toDdWX8MkfPJAtO6zyQHFAK3Go+eSp579WFCDu7byS9MBByMZDFra6FrCrp5ZmGC0349KWa0UitIdh+IoUCUZQiguFyWSCQzFpcYTRaDQ4Du0AOZ/PrdPLEScXXrzyOfgLWdhjwmutHpFM4CklaS6Jf3TXtA2Vu1ved4r2RACoXKA51fDB+ujJO/yPRmOmDnSGtq8QHZBFMWztnYw+hu1qgvBBkqB5cekFIYARTfQc5EopiN519Oq0b0L62GxFf7j1rSbikUuDEbna1l3Khh4mkwn0bJxMJvP5fDqdUs8/DBFpDv6klJrP57u7u7AibETVr6NQFy13oa4PqsRAlIAnZqGFBqz1pP0Fp+PqTNcKSdAxGvy4whBg56xOBVa4+mvQukZxwZHE8aUCu2awcg9Y3UWgBrTpZp3gWvjNf11xoHWFYNXTIaXagSL9EeAdIL2XU6Lw0irWqOMPQNdjLTwMAUgwQilVNnFjMBgMh0Pas3EwGIzHYzD59LO5oiJnaDQaLRaL0k3Vth+duHDvxW+EDDGljhtac8BJ61WO0+kH5ZAIlVvzt83Mf3QVmgQpYLoQQB/AdPrZw/KVrXQXgsvg0ZS3WqMY0NDD/qMf3v7a9+hcFAeaUIC1Nl65C/dL9xWteKJ7UZoVjW5+8PXECo91qfzV2vGhygDsilnukNobzfzkykIIgY8gCCTugKBQ0NR/Ssm1kJa867vfs/o4rZjPSXiRwx4QtbTCnUYGvz32LyOyeOjB7LZAdYOnU8N4PMbPICZqlVKoi3leoxMXwl0LppfM79inxtu6inYjmQvT7pr0gxZo0HZqLhA+KmalMtCoqxJqRc3p3PC35/AEBejYGb16+F7ighG39zYOtw/Z5hYg4JrGRnrKA5tJBoG1hHM1svkgQuscnj+d/foJOXr4lnLm1ZcqUxboo898ssHcex65BrohpP+kWu9gRFmhsFqtTIcBiAOw+piIMBqNUAfABzMkUVQouEiRC2buAk7E6S4fvnW6x23gmWLdAtw/KAKsG6F3pmcAC5ACVsVQVyUEUi5qTgdAClzFFAfmuqYPgC7TmNlGyZJdu+TCGlnQaKujXf9kQeIPAWuNNjsxfaHywFpLpLjecMw3onsvfgNG54EnFf2sgrXCnab2y1sQQqO9HsDWUvfAdDodj8dQsHl3d7cVKVAJ5i4kdovwDE2pmWHtfZ2abVeOIX7FVVwhCde66rifwOxXabaTLvDaJz7k2hFSVzdYvQXZbaq2QdPMmz0RUjom8MmxWJ470PpV3vxgnnJ4lZgn0eM/AKLHEai1et5CRs2QmF2Y8hs1rQAf0tvjmfvqRz5Mv5rvMB7nKCVLT/X1cS00VEcBExJnsxlMAScB+hugc8Tu7i43sUajTSlFF6zxCJiI2Qz+XgwhSQae8EFlQQV1PELhTxumX19T1RIhhCad4ZVVolXNGj6a+bcmarkkQsPSQXMnaIeidNZCYOin0q7X0g3+hTmLABOXDU4c/D0RLHWAH1JSHaOTLsPztzRfwplXX1LqfpxSt+hC/5ISNIoLBei/oJTS8hZHo5GWoDAej0FMsMKau6BqyoXEDhGmwbaaeY8LwZUXGbJHcxZo9rq1mEx3QsNhco84wBxDs0tCtBU/sbN0rZlFGVidGSmjRbCNR2QEBFyHlEHEq3mKYkjsFpGXRLlgxRXkRbQYRKBo6LdKUKWFAvZyrBVTwIwE2ltSRZVzLke4XEhRCS7bT+WCNcPRmtxQ2fehMlphbrkWLSbQLc8d0D6Km/ub2hBEtKpjuiHfOHer9AhPVqh68DfAKgtqVXikXhOr6c1obzLWP2auErau7906u5NlBEXYVMSKcRbammBIf0iuxvz8ffcP1DKkJ5ffO6ttAXqr0cVqJS4AfXUtlBUKHpVg9pzUFjNnsRIKilwT5Yaj1FJ2tT4/NFThUgbaBml+ogrIPLAOKanJEU+1Vz6p9RHJBHEFDyL24l8gug2b+5uBYsXs1Wm6Xqz2FbVCrTERKkmva5SyemNUBhFy7SLFQqt8gg81ULSrw/XVhScZ3J8GHiER7jSsjypBFU1mREs/OQ5MhwRGmqOwWCywSyTMxZ4RULWpXFMj0HIXcg1HaWJNKrR+ULYbw5OI4Nqmq48lXR6jDzi3seK+IR7yWsbV1dOhkEqoe4hSjmR4xqWZuKDI4zhEBGhjKUW0NpquZB0CWYobMt+jpyUqkxICteHv1+16eCqHswG6QtApdcvP9JXiOQqLxWKxWNAp8/l8Pp9jrUbMS6BJDDAXkhuUUnWDFw2QJXehLvhmb3YOdnVScPUd8iQqKm9QQ1s9S/JwOC6VEDfOZLnRKc1X+cQdJZZpAq1QKxqy2t+4/YE3Nl65Cx0MGX0GKfCXAgAHw6wR52AoMQJkdFhEw/Ps8i/sz9Bq+LHGn4JCYTQa+e0lOgyUuy4Tz6CDh0S5UFnS0WW5PXPNz1qKA25EM/9njJKOrd8/d16OHQZ+RcZTPrGzPFYNmtjXogEFaJjlFT6TKEkZ1i88GIGAMkOtULe3YTSVCRDcYCgLXETHI/KmOiZWgKZoudWe9IXAGrVmhiOkNNYt/9wb2h9m2i8COiQREJq7EN4twkylCcSvf7Uej65Z2geX/0A5REYDWN/+3zIhO0uqEjRw4h31UEA08MnG8BPhWqgEe1RCpkIuQ+IaD5oJHVIGVhJTHZW3cnZEM7L4GDx1GD11mfzSwcxXgClrpRjWbpjp0mi5C+lDV4dgFk+0LhM+3VQMWiwDPRNtxfBoNNpV9bncsJMazauExD0GDRD1gTfo1+W5A9MSlCvQlGW0xrxgpJ9JvD+dvD8kOkEEm1G3MT9/3x1TrUkB7anlKUFLF3PpCXzK3fPINe25F53z2Dna9yj0CdOFENGLEj9gDCJwiKlArCLaWqjEXKB5dJ+B2wHgH3nS9dWDf8wI+NBWmcX0/SYOR6kcKiGitFHeUowZ6YcaqCQxEmFOTPE1pvgVXJlbniceQAMW2rsQfr734jeUzfN6zyNr4VoQoZANT6Ahe0nHSjwlySK2Y2YqZNcN1ofOxrlbS8c7iksWmIsFNiDuHb3WWhlzJyvLTocTHozAwTDhq6foQq03y4zRiji6KwgyxgEzZgz8/H33/9IPr6VoBfigNaaykJQrocrTgdzcgrY1M0mcKob1QYRCNkIUAAQjAhc2JUKia8GagkAxcyHN0ENGxRBoTjyaIFAuWOlKVkEIWuWoiC3UlQv069beycQcglpZcukpCx1VBtaXYPyqZSPF2emMcgFI9C5YGxNSerJW2VkMKIT0ESNbuB+Wh/Syfuc59rOMFE+oyyGxF2WIXDC7OKZ7ArRtwgeMFHoe9HWf7KAArP8jmp1XE1x/8v1nf//vPAt4THVgS7Kkl0e7GeKCESf+9u1xuzMxfziWXtZmhVxXnZYFrq7/gcHBaCOdqBW2ru95ahjENcZ6Hl2/2tVd3LWA1kjzKWfdvhakOPPqS/B217/6jOJRaA566TRTdEEZLgGK9W7x9yyybsdK9DsfSgEOmiCOwJKOlYeILhAhGqIrS0b0osyLSye5Cj1pR7LTyiBkonWuab08KUeVJA4xZb6iRGzEbIyJ9XnlUgnU12LVAYHvUZW90PtHP3s9wFiUgNnBcjQauWbN53NckQ5hVYgso1cHQqMGZn9Iuhh+vuEY0poug+4ElTZerQuQC/S/laK1IOviN89xhyjukEYfk7gRtPNCj1KlTupEfwRX7yTrvemB3r/WeKK5tWgzFn08zZ4IzZtS6xMspI+YMka/c5E44F9X6KFQgCGth8PhbDaDUtDUDI9GI6gVbc7CIaxms9lwOJxOpw1oBRUlF6zpC/Dff+HeMCqeWj0Nrv5C9AHUVvGlxmpFF6L5xL1EudCuaPCLKjiYPMWBZpA0vzf901b06HjrLsyJrilxptqvwKyzXAKlYblw5viAOIqUnw/xj7YibnjSQ6EAKgEGn5xMJuDYB5M/n89BJcCoE3SWUgoqRsOK8/kctEIDDcYARC25gJqglqSlt42ZmUiXdBVfMpfMyGD7WL8AHPq5W5rAKmVab3zrDUjBLxfouIXWMQwbQzPJWha9smXmV2baa1CdEdiYBtDERPY3ilq/xezgoDk2EpWTh4z92FnRQ6GgqgpCUz+BNooEjkqliLYo08Y7aGkvEXKh1nTl7k+suS7NHGDtKVZOblOz2l3bVolnFCUoR+j6mkJvghEhlBMN5s2ibPob4+Lmqy3F7MhHw+ohXkAr1py7FEIiO7hAeLMDFY/ngUP3FeiPidAKfn9tv2MQPRQKq9XKDBmAdLBafRizyhxUwrNKRqyZjCgXAjeiuRbCPQ3We0nL8XEt01booUUypl9YVQIIAtfXjHRdLlinH54/DcoAP+TFDBZQPe2xN67ontbl2PThmQ4JRW5PT1MxuO5xE3qEi5/KdBBzuravWrt2vdhYcXVVaMa50mOt0P9eD2BuQTpo41iqo5SF5lsVQq2iC3FU3nXWZaIfMR3CLwhyjY5DM/ZbyV1QUb0o6w4YYbUoKR3wSh+08Gs7sDORaxee11+z611gk5TXuNIcvSw+hsqWeN5GtK/WJbXiMfBVk0FxEkQIpyGhAMNGK6Ugb6CxnUKSwWw2w4nD4bCZvadDB5dSJeWCBj5EGlblW9f3DrdbCDBHOwmiuy96gg7+vYRvNpxEuYDQ9UOyCxM74GncOrszUJafcHj+9OZPLKWmrdM9nm0tsYC+kdPceNOiB3oa/IZTuxld8Qv/vqxpE6aGcG2hLuYRU4YUsDY4xGXCShBQR0JfnQpNCAXqQl8sFtPpFNIJy+0R+i8opbQdDYdDtv4DDW1wKZU8dHV4VUcaWK1cWDMJIQ/9dnPUS4xDmPG9NqV5Ke/ZKaNXA9GnNUUxeH4yRh80TaB9dVlH+kLvt/0R5irQzvlNpkeduNayGu/KHflbGK0taISFqi5zSXx2YVA1Y86gf2uwR3h+whP4fReeUv0VBC7K5ihAWQKl1Gw2Wx0xHo+L9jzEXo7WZAVzYfhgZiSYWQtNYgqCuF6UNF+h1sXt8oV6FvDQYjd3Wt2B4WjFqkD9iejttNszIuXy8B/Aw/OnQfiCetDspTqeKKBNceG3xObNUnlDebZpNqYyR0Ed/1GeSMcZo1OAZ8vh3LDVeDDzPNTx0+HymlR2/K6Lq4e5B3gCK6V+ePXzGVvSCcoKBXitX61W1NxOJpOiPQ9BJViTEM2ODIvFggYj+AgFF43VaHKlJijyCLA+rbTReJsXB40pg/Rkw3ItTPn5cVphed+piLVMEi8Y109+9SMf1j6jevBv0GOSw3sf4AeP2g6x/R4do2kUj1ff/9m1rgur7ac/2SWeKC6FlDJ6Qoi2MBPA6cuVOaIvfG6sSh4fCgoFsMq0w6E2i1rl5557LstOcZuT41Crj24G+IBfoQQTdqSEegxZWpUXiErEXayB4TTP64v/cda6wyDONIaY/BL9Fcv1a0BSSjq27l2Iu5xclwHVCkhIVoHVHxBuTWlCg2aeAzfiCi5o6RHWFc2muox33hwFE8+Ptc5K6YtotfHhWwv0LqyPXCiYowA+g5AQw3A4vHLlSsZkvcVioeUizOdzUACz2Wx3dxfP7mw2Q58B6AnwgiijxAIT0nMXskT4ooOaHm6d3VH7dz5nHF/KgzlMQKVWCBygoa7tTx8R0UNKH43oASNysXV9Ly7bURtByjWYtStsj19rXepauP1MnVHZTDttrmt1A3jyE+lc6tW3/lLzJ5gb8VPiseDC2hsiZEkVFYq10kqyeSuUTWa0uhOUUqPRaDweg4V+7rnnrly58thjj2XZ42g08p8wWMAVVqB1mbgFHQBr7oKqebHW1Qrm46mZZ4HfdmYxrhEWHT979l5Lf9DFGMoF9Cu0pRiixz6um9155vi4QRqmcTVnud7dPeEATU/U7Qlp/gSzeWeOumbgAtq+6FetqVovTXOPmO5ntqE02q7NlmhLFmId5EIpoVBpa9HT8PjjjyulvvzlLxdqiRW/COApEfzUvVipa65WnwikqKMy0F6GWNYsjn1tL9FDOnka00wdhUTvgjLkQmMuh/TOES6ngqrZrdHEs3BlByKPv8HTl6HSnWBmFJ0xRjmy7s6jZpp8SQiHSQcE+gTuH8V7PfgXAHfC5cuXizZjTYjIXdD6RIQH9hp7aaiLWfw4bjuFMiIrt1Y6ZSGwGR7MxAVMaMiVz1gO/OHwQStUoBzG0uMPoMu4igRUWlbN2NNkYWt6hHXFM7ay63Su1uXBlWKpfbVmTcLfz755f61UgFpEbzB7/whBlRMKgfWPH3roIfgPtu1LX/pSofb0Gzh6NHehyUQbmvYcLSCiq+5YzWrdJMEsGZHhZJELmOuX2KUwYkXUCppo8Cc/pryyU9L70eCvhqtO86678hatFjSwi6NnlilWlM1dd8ZRiJDed1Z54Uql9PxY89SEn6wUI905Az8wwCoAbTctMwU9CtCJQEtmHI1Gg8EABAT0dLh48SIc2ccee+zxxx8XrVAXvDq16XUHjADqjkXJxLVQ9y2cQ3GFELlAv1p7nKJWSJELEWv5+0R4VCOdlXLxNNnz1ho+0F7QXT8q5B4x16XZA1qWg2tr2Eh/EqLf8VDpw0ikc1KAYsoCdfSY7as+QAZFfxuMpDAcDsHBAP0gsDeBGVAfDAYXL17sSvFEDmiDT7qWUbFZNhH9I0KeMmZO+M/fV6/DdIQyqLW8hxt//C/O/Lsf5doaUPlzrHWIXcS5Z7KEPN713e/Rr9oruLk8DZ+n7DfaI/VLP7xGGxOuXTwLV6YuWj8rh9vAFR/Rdueabu7OtXxdPFbffG5gLpQ5l1Y/xLmtSwrzFcv6CA15AveA4jkKs9kMyjaDShiPx9SdYGYnXLx4ET8/99xzueor9JXAazRLJCLjrdtYpwnOBRkpeRuZt/xAHK6AvblfUbO1AAAZU0lEQVRAFqdUonfBn4VgvVy1bMEbtipkriiGFX+OoTJ8GK5tmuLMzGEs7Qu0Piv8E9sdMcETROi9tyCE4sNMQ3dEBCMRDz300MWLFx988EFcEoIO0P3hueeeGwwGDz30EKQviFyIBgVvrsQF1z2sTadPIi2JoTJ9OuSJr1k1M8mAvzgw8Rchpn+Vm4r2zCcePdfrr8dktptFj64sGuDXlokWEMqw7pUuFqtzhS5J76zADIlcsiDOfnuyC8P1RHb8siB8OxHh3S5Sto6CC9AEWu2Exx9/nLoTMAYxHA4feuihy5cvU1UhBKJd9NG9KM3prqhEYM6aNj36JThiLeaE9GB0DYpokqtgUV0qPQouUsIQcV0oa+3O1TmiUg8pRzUka96idddasoJyOCFcnolE6trvlLqK2QmMIwguinsUXFjdCZid8OCDD+LnxWJx8eJF6B/RCXAoLApkcVJomiesYk4vRK4K0Nq7QrmOSe0OHlGLXC2s1EDgVwj3LsS1wd8M61zq6NasoMeO5vWHx52FkG6Qlasr20u/5yBYMyVNAs2/ObeuwQ68hVvPIfDgiiMoknvo34I4FTTa8SiY5ZUef/xxdDBAReeLFy8+9thjICbgKy4M/SMaa21dsAg0hSZ1AvgZh7uEYSnCS1+ng5d4ir4OHzPCivV9jrMUoHjaSWdFv9PDhyzehUL1DTX/hzUFDxe2vl7ji7K2pJl8F4HpUDHPi+cqvWEUTDxTVbbI+hNcB8G/dxfZkxA9NQ0r122dcG8BPO5K+BJ6759gkbFJx3qAs/7YY49duXLlypUruAyGHi5evHj16tVnnnlGKfXJT36ypSZboONEKOPSGQwG4/HYav41Uw1dRUqfF+2GSZcLSN2OElbvayXRWe4pgI35+bc+8kuf/W7ddRMbXCsKALrBIyByVa2gHgWYtfmT183Uek+ufmXXAG35aPAnm9rOLxTMBcyJIf0ONN1jlUclsNY5RqxdDHCiNphC0XaGkxJHWJNOCtlpLfRAeeyxx6D7A3aF+PKXvwzGErtFgEp49tlnr169qpR64oknHn744cFgkGuQiHRgAAvX8BbKWxmarmUOrVkCM3ehyRpNlNIPymjyxjsSN1I3wdAfkigRjMBZpmtdOSoIVQbjzQJHEW1G4n61Jgjq9udsMo2AYsYBozMK2611GNIfIfypJSohDhZC4cEHHwQdgIEGbQGUC0888YRS6plnngEPxKVLl5544olnn3222fY6gdiBOTg1jjUFF/poNEIdYA6KEVjUsgRZ5ELIY8XV0astTE1QLvbRgFwISVxIaYmnZwS+uAf2cFG2i8GzcLRlrfyl1h4EZpiAegJuHB/oWdnEjbmLZlwIIRMrt9C8RJBuijxhIRQokH8wHA6/9KUvfelLX6K9KMGdcOnSJYw4QKYC+Bg4A1Z/Op2Cy2GxWOzu7jIcwxrJIhc8czl4EdpNjWzSu1CuJVZcvf5cUQZtRZehtaYKZsFlv/3hD6uD5Iy3MGKhy94cNNk1BIMrO6EtPLJAHb23hGynAfUwmUywkZpvmFWieiHaSWb08OUvf/nixYs0QQGDC+BOoGmM4Eu4cOECfB0MBs8++yzD/hFwYeFVAtfc7u4uQ3U8IHUXVFruQniGVAS1OsKxTY2MTjBUyd0XzZY0k/YRGKHXsgi1Op7KpioqCY8XuPoxatmXqsr2cxDEfKiVW1Au8TACyDGHfPP5fD6dTmnbeCaq54WdUFAkEvHcc8899NBDV65cefDBB8GdADmMyMMPP3zhwgVwMEC3iMcffxy+Xrp0iU/PiNFopInQ8XgMFw03XHUXCt2xnmzw9M70XaFQf4S4ZsS1JASzH6DWRcI0w+p43wFta+YWKgnUCp0w8GYtk4jB4ouiKYO6tp+JSlBKgUrAMDG86YHthwVgorkipLdj6Hk0Gk2n0y4KBXahB8qDDz54+fJlUAAPP/ywOt7NATwNly5dgq8QgHjssccwd4FPnqMHMyPBzFponZRIRKBvM+XR3JUSC37yDuwEUwLLOGZsiYfKU3zj+JiH5oraRJolELJ9uqK565BWqZIREIpZpMScq5iFEirTC1SXSw6YD2T60GaVqF4C1kJBEe/CpUuXNHfCE088obkTVqsVxB0ef/zxCxcuQKiCA2ZYS7tWmAsFRVwLKXKhxFOsu8rASpYEQ/yfks1Q2RK/MrvhGIfQnEiLDVgDCtbtANbMgBDC19KWzC4OzDQCV1kzq3Qod2f5WauswxUZfACh73j8E9UT4S4UEIwpANSdcPnyZUhypMtfvXoVcxdaRxtxezKZLBYLVJowFy84cHO11FIn9EFQVy7EvfdYE875+4S1SHYceb0LKaGEkJZULuN3+Gu1laxDJFhxpRZGiIZay2fHlXtYd5lyhHgLAumuUwGB9lP3QIcS1ePojFDQoO4EyEWgGQnQnRJDD5cuXYKL++LFi5pbQin1rW9961vf+lbR1kKfSUiBGQwGcFVR3TAcDnd3dyHDEYNhfNAii4lyIXCYOJjFXB9o75r0s6sITzgZQwCJxzBQLlh3amYveloS4iGgSkLTCtq5CHQYqBxaIU4Ntztkoou18hbUArITlFKz2QymQPkccDlMJhM4MtbivJ2GYzJjCHilWt0JICAgDAGVHL/4xS8qpa5evfqpT33q6aefhowH4JFHHvnoRz/62c9+NlfbrGYeHQbKHe7iGXRQjqyi9G4RVmhCFkN9cObVl34eVoVXC6gn5mbmSjCM6CxAW+JvBmoFtOLUWtOJdSsSepIWtb1YT02T5TrgAtYUgFbrkFVugbJlHeL0Eh0UWfVoCAQL72oFdjuUqJ5CV4UCAr0oTXcC9Jy8dOnS1atXqTK4ePHipz71KbxG/+AP/kAp9e1vf7uZ1vpFAEOJUElduUD9CpUFYlvP4q6V9RayqXblAlroaOES0gyPK8V1EKzSwZqiaAoCV8aDttaLVz6n1J0RpWtdVKaZx+na1kwd4FEGLRLeTbGcRe+iSqjr68WMBNpbUnXzOd/V0APy+OOPa9fcJz/5yQsXLoA74Stf+YpS6lOf+hRGHMC1gPzhH/7hRz/60aYa20/QMxmd6qjRYslY6rKu9BPUJb0wX3okwhooiWhGZUvoxrUfbroZQrDGGlSYd8cksLxxIOEJBBxqHSoyiOI6xxHC8aiEHiSqh9B5j4IG5CXAfxAHTz/9NEQc8DMu3LA7oX9ojoQ470Ld52b2t/y6s6LR3obLvdMHciPH8IxWrOdIcxuE5CuYS1q7TdYlwlRbQwmF9pVCoLegUEyhl6Cl1zo+QNABAg1YU8FMVJ9Op+BUYJuoHkLHAkWVXL58+cqVKxCJuHTp0le+8hX8gfBVKXXhwgUIWAwGgy984Qtf/epXW2xwd/E/a3LlLoQ8kesOz5PID69+/n0XnsqyKZVmpK1ywVo36Zd+eC17MzzxC7+Y8+QouGorlUg2BKxhBc9akH+gZSG0iCu9oHKtnj35y6ENC4zQEkyLxQKnm0kMOJdhonogfb5cnnnmGZqOAFy6dOnChQsPP/zwAw888J3vfAfmXr58GftesirpyJbAB026XMDnsmeZcHOVhbxCQeUbQBnQggIwt1IoRDTGDC7QWZoLwd8HpPRoSbWMOpUOTNQA4nEYcLP9WLFY2Qwkmk9zFjXMmtHljD+s0N2gA8Dr2soONVTPPPPM1atXsQoTdScMBoMLFy489thjV65cgY6X5giWQi3okc8iF+hXq3TI8t4ZQiGhkNIzQh0JAmvqwK2zO+FCQSV4aJoZA8lDuGPADxOJEJ51iMszeZ6bgyMo0nhQCeCf12ZpgyPAYl3RCj2Gy4VVDqtrznQnoCMBvvIcXKorWB9YKU8xUyiYE5sZw1cVEApIuV58dQ9LpXZxbdClFZocYdnUChysfiVWWcDH9tdiMBhoroLBYAAmH6QANf84SxkvFSApungEekbnez1Uslqtnn766S9+8YtPP/00XnDf+c53vvCFL8BnEARPPPEE9KukE4U4XHUXsnSLcD30OddlCqT1n2AOZFDZHu10uJYv/dNc3WvL7TEFsyeCq6hRd+sY+qvFUCeBJin6MThCz+hbrwcrDz/8MK2w9MADDyilIOhw+fLlhx56CO5GdCR0VMCCVDcbzyccGB2JMPusq+NmgNWgeYl4Yv9N7lqbaLbEo9iyNyyCtgZBcFE3jtBprD/NM9IBPKA8gyN0N7rfD/rvUTB54IEHwJ0AUQboSwlX9ic/+UnqV6C4BhLlgzU1F8OBUIecPq0wHDibzaDCdAM/EFVCrqILlTRsa9PxjIHUFn6VkN0Yu0ZFMhezzuUgDiodBuGb6q5TARmQwRFoBwFARAB/1lEofPWrX7V2iYS719XlAcdt4gYOXGadhdlAWIcc1QCOlQ59fEErFG0qPi5V8oAR5vTKF1xu1teFa/yIuuvmbRJu/MyrL/3sm/drCzRjm+kgiv49tl7UyBpHoPTA/IdgDo6glOpoLYF1Zh2FAvLQQw/BaNTgVHD5EhCeyheKftDAHsItHGh9XDbpXego4UMc4fLalAgvRd0uDNG22eoPiBhPuTFx4JEF6uiSDtlOuF+hi6oCXmBgALzVaoUPz+FwaDoVBOastVBQSl25cuXSpUtPPPEE5Cioo8EjrGjBMz63LjgMTJ0eFw4s00YfzcgFJrHzFCrHXXR9zbsvF/4QQKUPoKNBBGrvy1n0biU0YFhztVpVBjTxmWM+grpefqA3rLtQUEejRTz77LPPPvus627UrtfBYLBYLKgzjScdCgdmyV0IT33vShjCxOVdsP7MM2TU5sZ0UuBLf2WWQ1sqwZQF6vjgCP7Vu2XRC+EZHMH0XC4WC/qSI0KBIWvR6yEEf39IrUOwUmo2m3Xi8u1EOHBAOotH94xAaNUd5pn5KZidEVofbFORAIFnaFDXCKKtENIZYdXZegat4B8cAR6bUFABl8ElezM4Qs8QoRAEXvqdUwmdCAdaH80qrSOlCjZIjZV0LEGhUZ08BL7ouwZUbDe4sFZ9FNtlsVhoD5/5fA4P0tlstru7i+eCPk6xIhN87e7gCD1DQg81QLvVCZVgpVvhwJTchVoWiH8kwmwhnZKxG4LZNTFls+0mGfjTC1TNdMKSLX0LyAHUJkJGFIW+rGO/J216W4xGIzONY7Va0YfParWazWaz2cx8nM7nc5wrKoEJ4lEIAqVxt94/JpPJ7u4uLVfSuXAgzV0odPA76k6wNts6GnK4h9+VK9B6gCCEut4CngEFazUUuG3N1GNlDI4APZw5yIVK/I8abg+iNUeEQjVgQbvoBOt6OFDLXQiPRHjKLXTC5pnEqRkzEBP38zn0PjDRlAGGqxj6CSpxjWWMuAq+YTUUdVQTpZnKacJaIaEHJ5PJBFJvPBm8/JnNZlCQEfo0a+HA4XAIwUKev1F7jidGImrV9G0sGEF3VG6nrrIEgQtzoLKbIizWxZIDylsNBRdwzZLBEYTSiFBwAsWPobJhJ248iO1pE/sXDkwvulBZ10+1NOBCA0EQq2IIOSANY8oCdbyPIh9nQC781VAwEQH8f3QWk2ooQo+R0IMTMLFttyIP/QsHUq0QfZqKRiKgMwUdr9kz5FIcie33JCry7IwQGFbgmXwQDVh9KHEIH3Z3d7vS8UroByIUhO6h6YN0uaBI8N60vh4zXwkdtSGuednhkKhhTS8QrGCmEXyF0ROs48QKQiEk9CB0jAEZWQqIC0bUemmOMPPR6YeeTEzldQPU2kvEWtF40gtUmT6KHc1UsGKmMXpSGQShBCIUhC7hcSlnkQuthOpd5Qq45Q0EUpl1qIm8QhZ9HV64O1ENRegBIhSELhH49I/uGVG/RTnphHtAwyoLlAyOkA9IYKRTtFxFEQpCaUQoCP1Be3mtJRcaMLcRu0DXAvUxtNhPIdBb0CfPf+uMx+PFYkFzFKArFp2LPSMYVkPpE1lySAc8CmjWQoSC0B9WWesu5GvXsc0GhhVcYyU0THgQQSgE9JmcTqdwCqD7A9UNzKuh9Iy64wZDv1ZaSVN1pHQmRYSCUARafx6hC2D5+tJu0mjvAn1rjxjjoPK9X1ug9cCHChgcIZzGNIR1cATlvcC4DY5AsVZDofVOVquV1uYuVkPpKFAH0xzsCmpgWA/+aDSCqndwHYKYa6KtWZHukUJBXOnZUMwK+4UPyvd6R60QsSPr0MnWKZ7V/VNawbSv/oOzYlmfwFr52HOByeAIQiK04BW9/LTStwhohd3d3Q6XwzJ9iYKQDjyjrbPAcTcej3GK9rUE2JiUy/59F55634WnlFLvu/AUnQJ/eRpaDNe939GHgOb+NWe5LjBteXi9a6LFQgeZHYFT6AWjXTz+5xhesaXaWhIJPQhF8KhmM0pXOrA6OD641Cp5wAg6RfvAhMCsQ9XZxEPP4AiVF5gMjiBUAvGp3SPQT0AdBtpgvCsjKmSlE+4rDREKQkEwTkzvDetDWQv75WVlvE+nyAXFTxYoW9ahqtNHsYv4B0fQgAtMBkcQNGDIzYExiAbEp+DeoR1P6PUG02FduBo9O4JrFcbs7Z5WKOOoENYduLrgxsBbi86iC3viFA1Q90bQFm4l6BB4I4f/ru4+CszYgecCA/cv9SSvGol8Ca3z/PPPP//88/gV41PwH51MMNd8ItHLTCmFm6LPN0C7uhCchTvN+wOLIsmMQhHG4zEM0g1fQaGj6OaT94sv34PYASMa8C7UzToU+FxgQjrw0Ih4C8e0gI9//OM48fnnn8dQwnQ6pbcSZLYqpebzOVxC8yM0l+disYA0Rq1V4D01HVSwGB22o1tOLAk9CEWYTCYuB+9wOCwaaAgnV+5C9lb5+yiGq4RaS7b+w3PB5wITaoH9DM3+n2jC/etq42IsFouPf/zjH//4x+H1/fnnn1dEDWho9htiCru7uxCbGI/HNH8WK1zRgIW/eXR1V19KvjTnvBDWG6UUxPxMXzG3zHP/rZG9qXVvSVbHqnVCLifNb0y9vtZghNAw1hJGeFIqQ5PaiphbAMoAv668fRZW5Nqw7hG0wsq4nNTx4EUvryV54gj58UeCzbnazcwEl9lOtNNWsS62Pxrzce+/wLSLrd38GGF1dL6GwyGeMq0nIXxdHZ0sq4xwqQpNF9K55qnHa8PsYetakUochg+xXMgdIhRBu5nNDsc0o5ihDNfeGKyzQjZilQX+PQq1sLqjPBcY/YomqqnGChas94Wm9tRR3glUNaCrQDqhZ0Vq78HHAKmI/oxF+AxeBFiSOiq6lYqYjjybhCKYvkQqBbS5bFUCnUKlg2fFcGVQuVMhBKtQ8F9gNNVRVELrqIAuANqtZPVZzmYzPLO0vpZ2ipXbtalJB60nZO0f1iPk2SQUxKxrZs5ttEFpaLY/RBOwNf/WqDBdAJ+S3TWl/gusc5dfLwnMETGvT2q88WKGKfSitSYi0DiUmbAS+1P6jHSPFArSs+r0q6OOlOF9KVcsx0dA+AzGUYKeXX7rjKcnCxZHwn4EdGG4krVNWT/Dwl28zhtAhIIg1CCi3ALPR49nuFvoNY5jGU8mE+gd3r1yckLHwUqaVntPexjSz3RJzGzA0ZuoOO5YH8X2kDoKghAKljToQckBVoNxCGsLrUFkAmUMlK3StiYdoJ4BDNAA4mAwGICzYTwe0zEeYeyG7D+k34hQEIQgOup+98NkMA5hzbFeXZos0ObSUZpgzAVQAOPxeD6fr1ar2dGIz+IJS6eHzz5ByI6pErquG8ApMhwO4VUMntSYhKGOR0wmk4lW7FYQcmEmGQDgM8CrbjAYYDgMFxBHVzOIR0EQqumfjYQO4lD4FuvL4lNYxkoQGmM0GkH6IRp+KI2M6bQADB9DVxSV0BjdfisShBbpulNBA4K7+IwWj4LQJHCN0Sma/0BoEen1IAiRrI/hlFc3oTTQrQavNOm8ygoJPQjrBeRFmxl8awUcBNP8u4bBXSwWEowQGmB0RNsNEY4hQkFYIzBtCuq3T6fT9dQK8CDGvmfquETQcsWxmkKjTRQEgQ29CrIKgh8t+m4G49cH0Ex0CnYnM+fSWYIgrBsiFIQ1AvP14CuYw2grSK1pRxOv/CFhWutGEIS1RZIZhXXBNHtY7i3CFmIUAzKwIGG7c1pBRkMQBKESyVEQ1otcxg8z/kaj0WQyGQ6HWucuQRCEfiBCQVh34vr+aR0BzJ4CgiAI/UCEgrAuuHwJcXEH5Yhi1G+XIAgCa0QoCOtFRlsuIXxBENYBEQrCekGFQvYKAeJREAShf4hQENYIGI4Wzfl0Oo0rOJgxiiEIgsAc6R4prBHgPKClhFJ8AHH9KgVBELqFFFwS1pH0UkJa7SYZX1EQhL4iQkEQYgBlgFUdNd0gCILQG0QoCEIkoBXwq9xKgiD0EhEKgpCEDIggCEK/EaEgCIIgCIIT6R4pCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOBGhIAiCIAiCExEKgiAIgiA4EaEgCIIgCIITEQqCIAiCIDgRoSAIgiAIghMRCoIgCIIgOPn/AWo7cEMl4sEgAAAAAElFTkSuQmCC\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
}