{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2edfdbb9",
   "metadata": {},
   "source": [
    "# df003_profiles\n",
    "Use TProfiles with RDataFrame.\n",
    "\n",
    "This tutorial illustrates how to use TProfiles in combination with the\n",
    "RDataFrame. See the documentation of TProfile and TProfile2D to better\n",
    "understand the analogy of this code with the example one.\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "**Author:** Danilo Piparo (CERN)  \n",
    "<i><small>This notebook tutorial was automatically generated with <a href= \"https://github.com/root-project/root/blob/master/documentation/doxygen/converttonotebook.py\">ROOTBOOK-izer</a> from the macro found in the ROOT repository  on Tuesday, May 19, 2026 at 08:09 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "67fe67c1",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:17.713156Z",
     "iopub.status.busy": "2026-05-19T20:09:17.713029Z",
     "iopub.status.idle": "2026-05-19T20:09:18.677923Z",
     "shell.execute_reply": "2026-05-19T20:09:18.669419Z"
    }
   },
   "outputs": [],
   "source": [
    "import ROOT"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ae487f1f",
   "metadata": {},
   "source": [
    "A simple helper function to fill a test tree: this makes the example\n",
    "stand-alone."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d618c0c6",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:18.679740Z",
     "iopub.status.busy": "2026-05-19T20:09:18.679563Z",
     "iopub.status.idle": "2026-05-19T20:09:18.803986Z",
     "shell.execute_reply": "2026-05-19T20:09:18.803589Z"
    }
   },
   "outputs": [],
   "source": [
    "def fill_tree(treeName, fileName):\n",
    "    d = ROOT.RDataFrame(25000)\n",
    "    d.Define(\"px\", \"gRandom->Gaus()\")\\\n",
    "     .Define(\"py\", \"gRandom->Gaus()\")\\\n",
    "     .Define(\"pz\", \"sqrt(px * px + py * py)\")\\\n",
    "     .Snapshot(treeName, fileName)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b36823df",
   "metadata": {},
   "source": [
    "We prepare an input tree to run on"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "c8136ea8",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:18.812425Z",
     "iopub.status.busy": "2026-05-19T20:09:18.812279Z",
     "iopub.status.idle": "2026-05-19T20:09:20.593897Z",
     "shell.execute_reply": "2026-05-19T20:09:20.593310Z"
    }
   },
   "outputs": [],
   "source": [
    "fileName = \"df003_profiles_py.root\"\n",
    "treeName = \"myTree\"\n",
    "fill_tree(treeName, fileName)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "09bfa943",
   "metadata": {},
   "source": [
    "We read the tree from the file and create a RDataFrame."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "3a3aeca5",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:20.595310Z",
     "iopub.status.busy": "2026-05-19T20:09:20.595190Z",
     "iopub.status.idle": "2026-05-19T20:09:20.709624Z",
     "shell.execute_reply": "2026-05-19T20:09:20.709018Z"
    }
   },
   "outputs": [],
   "source": [
    "d = ROOT.RDataFrame(treeName, fileName)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8b3ee2e8",
   "metadata": {},
   "source": [
    "Create the profiles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b631b386",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:20.711093Z",
     "iopub.status.busy": "2026-05-19T20:09:20.710975Z",
     "iopub.status.idle": "2026-05-19T20:09:20.960941Z",
     "shell.execute_reply": "2026-05-19T20:09:20.960548Z"
    }
   },
   "outputs": [],
   "source": [
    "hprof1d = d.Profile1D((\"hprof1d\", \"Profile of pz versus px\", 64, -4, 4), \"px\", \"py\")\n",
    "hprof2d = d.Profile2D((\"hprof2d\", \"Profile of pz versus px and py\", 40, -4, 4, 40, -4, 4, 0, 20), \"px\", \"py\", \"pz\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f9ac02ec",
   "metadata": {},
   "source": [
    "And Draw"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "c58158f9",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:20.967419Z",
     "iopub.status.busy": "2026-05-19T20:09:20.967292Z",
     "iopub.status.idle": "2026-05-19T20:09:22.151032Z",
     "shell.execute_reply": "2026-05-19T20:09:22.150589Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Info in <TCanvas::Print>: png file df003_c1.png has been created\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Info in <TCanvas::Print>: png file df003_c2.png has been created\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saved figures to df003_*.png\n"
     ]
    }
   ],
   "source": [
    "c1 = ROOT.TCanvas(\"c1\", \"Profile histogram example\", 200, 10, 700, 500)\n",
    "hprof1d.Draw()\n",
    "c1.SaveAs(\"df003_c1.png\")\n",
    "\n",
    "c2 = ROOT.TCanvas(\"c2\", \"Profile2D histogram example\", 200, 10, 700, 500)\n",
    "hprof2d.Draw()\n",
    "c2.SaveAs(\"df003_c2.png\")\n",
    "\n",
    "print(\"Saved figures to df003_*.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c61faea",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5aca09e0",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:09:22.158294Z",
     "iopub.status.busy": "2026-05-19T20:09:22.158166Z",
     "iopub.status.idle": "2026-05-19T20:09:22.343947Z",
     "shell.execute_reply": "2026-05-19T20:09:22.343442Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779221362328\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779221362328() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(24795,'WkwIRSQA22AAeAHtnW2PHLeV77/KoJGXnAIPn1mFvNCDtc69kiVYcixtEFy0ND1SX81Mz/a0bDmBv/vF75DV3TOSnIdNdmPsDTxCHxZZPCQP/+eJrPx58X92P12vrpaXq8W4ePFgefXD8ub71evnV8vrm3eb3cIszr+7Wv/Hh9XvHi5Gaxbn99e7m/br6ev/u3qzo3xBtafXu/XmqhP/e311thi9WZzv3zT++XN9fakD74PzIZjF+eP11erB5mKzXYzSyee7ny5WB/L79dnuXSMfrS8uemWYhZwrW6vNV+e7J8vt2/XVYrQDJd+u3767U3R/s9ttLm9Xe7G5vl3w8nwNE84szl8dft5rP3nxy5vdckcvtVLnFnWvUTR+tF1eru7yTdmdge/r3R7QvupcPL+S5sczwyvvb7Znq+3z9Z/67B0VPtmcrdq6vpTFeBoH2/6XQ7Sx8NZXlMvgg3U1SE4hxswKv3SL8XPV3WKUQUJOMdgaQpaaqb3b3Ht982z9cXXxfjH6UFrZ5wsWY6pDrLXWUpMruSoft97gJA/ZpVSSi1ZsgtHd/nW/8JTBlDpkCTFGn0qNrMZ3t16ODH13eNsxCWt09vtbDULmJb8/NLldsBhPW8E8BbvNy/efn22dkL/8GPmyEnzKRZK3yUaBg8P7X73/7DL01/P0c2s6P16Mp3awNrosNqcSopSoK3ax+fGbhw+awLw6Jl5+d60PdKWOfn+/L/16/+ve65tbL7r3+ubWu+69vjk0u/f65tDyu4+X7OBTAOK7n9pvGWyVUKv3IXqfEM3vPl4uPy7GXoufdii1lFhddDHUkmHzxbvVbrkYPcv77N26/7p3c716s/t2uVtv2jC/+XD5erVtv1+s37znbbRev3n/U/v5ePO2Fz7evD2U/ak9fbY8e7ZcXwEJZnH+YLu5uXm3XPcX7slnmw6vx3uVJW30YZs+2Zytz9ers8V4vry4WZnF+b9t12cfb5M/Hch7r28ebDbbo/pfna13y9eA2G77gRc8Wn9cnd0a9/zqZ9v15Xq3/mF18wmOP17foCZmFdLJ5Xa7GP/wR7PYXO/48bNZnH/1cfXmZjFefbi4MIvzb5rOeSOogBfrHXwsnm03AOjJu/XNbvN2u7w8WX1cXl5frKj0zYfLZ8uL1W436wPm8ZvVx92npQ9/9/zZ43uvFuPiN/NPszh/uPnw+mJ1/8P5+byO3652y/UVk9in4OXN+k+r727m569uk/r029XyYjE6OtfHjZYwuBJVnr5fX51tfnyxuWZ3msWBfnVMd2w+VPh6hSbq0vHjjC8P3i06iDxY7nafzP693a4pbcb38v5q9+NqddW10i1KZ/jRdnP5YnMNIiBSL8+WO0BfiVczwca51wj52SzeP9n8sHp6vfyPD3tBef/tinm5XXj+9frtu8cMoStgFdrl7s27eXLfP3+3+fGrH1ZXu+e75e7DzV4239/7sNsgHfuaT1ZXH+4vt41Gdu69Qfr2Lc6/XS3Pnl5d/DS3OP9+vXu3+bA7FtRZeL9e3nTRm0uOa/3hjl3yD7N+QJ0vWj/fr14rGKyv3rJyX7KxHlwsb276TqFes7mOC64R1YUdXYym/00yWmONnZyW8suPVp/aKegzF+MU923slHpNynOvy+8ylmqciClhqvwuwThJk9BhMP1vEhmlOiPFGcluEjfmZPS/Sfwo1pn+N0kYxWXT/yaJo0Rv+t8kaZRcTP+bJI/OBtP/Jimjc9X0v0nq6IIz/W9yVitLpb2dnNwm3SjVG6liJJXJ+VFKMlIg/eRC6zdFI9FOLo5SPE+MxDg5uOpMBju5PEqMRoIz4u3kyig+6gBLmlwdJUfjbIKTydtRvDPadfKTl1F8hgMjyU7ejeKtkVgML/a+9VNKa8tUudZP9JOPo0g22j7EyadRbIADZdnnUWw12sC7yZdRXNB5dbZOvh5mPcYp2FGyNTrcKlOQkSVWnn2cghudE6PzA8lUlcakt1MIo04xq+bCFOKocwobUqeQRl5Tg6lpCrl1w0SFPIUyao+ibE6hjhKy0UGIm6JVJlTEpiijhGRUZqfodB464cda5t9hzGH+Hcco8+80ujz/zqOffxYEu7+mjtYEBDhOyY7WpNJ+s3WqtN9uZFC6rabEDhJf2xM2EdKhzSNEyY1I9CGtk6Q7yfdetHfdpDIlulfC5SnTvxI2TBkGIKTYKcOBEjFNed7DrG+GA31iy5ThAKLEKSsDMZokU553ss9TnrsXP+U6Ng6NnYodG4P8ljHqMPjtxqLTwG+2cEMXOxWkUqcBom1f+rNTYZukmWDrNq7sVNi4Og0QbFudbWOn2oHLFwhpiCQBooGXVOCren0iKUGERgQH0VAMobJTbTBW9c1Za2WhvOjvkPld9bfzxk5i7SgiJkQFBbGiEhqcCW4S6xooVJZXJrF+LNmISg2Pg45DYjIeMo7OFeOziWESm0aJ1mQ2aJ3EAhlJ60riVWxPZ7w3zkEyJSYxqzKJYmtkD7lJBBgT48U4KKezChilOgnYym/+XJwEcGU1mXjeK7Fx3ORZQFcrplYj1k4ieQTgGI+USUBX5C5VIyB6HaVUOpZJnB0lgFnBSALfwbE2bwXSjdrMZlPqJM6Pzms/dhKgtUTUgmGbC9gq3iiUOSon3epOSuMZdD1SAc7KWIUtWWVyLAmQ2P4mZ8Emb/rf5FgTYTz6NzkbG94C65NjUYpXCLCTs3lkb2YxmTeX0fna2JI8OVtHWgHWcXIsCVJA48mJjFlMazw5cYrLvfXkWBT6o/HkJNCyNZ6cRFq2xpNjPUB0ybSenOTG7dxp2Q9Oe63K79wt6g44akxPTpcEblu/LEkbqxQ/OV0SHazzYXIu7EerPbMofbQw7UDxw6y6PI+3NS77AdPYz3sYbLNh8vM2RpjQW+zkrsUb3TazPrdu8uzn+bnSRzaB0m1XUz/K5NnYvbqSZW9BKAm2NXMjyhTswTCyU7DNJqA5lNu3hDrYKVABdYJGmYJlNhp6Mb6gm1mxX8cT2MzQaOPop4DoWLcfXxDbaJ5bNwWMJVhEgyrt6KsPL4hXygblX5QPqjIaUSVHY6VmqLdTEIX6Ni4B6mlCOYqGBnYKDj3ThuT6kjSbcQqYQ31IqLTApu1D0iEjH31I7TlI4nTI7TnWRx8SU+DKPKRWvc5j0treQtKBVvaiZFu1KfhZ6WlTr1qvrfcUvCrefUPVvPt2TIcu2hQRyS4jvCVW2S91ow9LD0exHha/PQ9tPG2MU0QkWz/KcqwJltuqTbHmmdJ36djbAKdYdehQ1k3JsgSqFadkm5VOm2QZsy7alCwj1kWbkm3WOhYDUth2mLKYgK+2wzpddMh92aaETuF5W7YpKYSpBLYukcImgTqkhGrRHdZeJ15fxxBhEM2iQtgrt/mYhyXH7kPqVjuoR8u2O/rg0CpIY5sOp4LAqk0JkZzHRzsQq48PKUlAVtthne6Kvy9/Usxq49OpVitdd1inm3prApJcnYdH31jpbXRKtf3Za/qDE5V8mxQskoR1zh7T32062m8VDX7mGQgb21PuQNinfcpIXWcxypQRuk7SOO7Hqw8Py8/wcjdvlIUpdwOHnrRyG117WGYMVHQoYGCf1yhTsQ0JdFGm0kFQQaLY/ZrzpLPTsKd0aez8TmWWRkWagix22dM+9pONjdkQuXXRhZDKVESVdlihpiJhf+UBCHnS4agNt6BHZ0C2biro0Sb2WrnsKe1jj0UYuSqAdA/RtUObnamgQmdc57WAJNLWwLiAkZ2kcbeEO0tdGvssTKVL49xT26tQsITP2DvC2O7Y1ZhCGPvY2BYFl7FPEdui+DZJTW6mcgBIOPaqL+Cp1e2o3RV18bN/0CofMFs7UhdlbhqOfZQSFLF5j743KGK3zTeVoJ5Sl5MwAzbjDAd24C4c2Gmv2bPDwMKenVZ5ZocXxZkbJfoMNRU9FdzFvi46kLgXI31vhKOOzqV5jX1CpxIbR33CpxLhSJ8SConK0YGEo46PRX3Irg9KagqtLdtUmh/ZFm0q6kh24CzNkWzcTiU19d5EYSrpSMnSv7qTh7e2+ZlHrR7lvmlzKednuYdCZoZy134zS7kJ9jwpuQl2Q+ySm2DPA83NGNrPQsazmT3skptgzxOBo9keIlGlCXZTf6W0ZZurlrZqTWmUomvWtEQpagN1dvA4Vbba4uJzdszUOSpt6899tK2/70O3fu+idt2jzTpSzzPQgbr3340DVVOlg3QnGijOq9kxer+eHaPn1+wVkAp0x2h9Ue0Q3YkDQjNv9QiikeAKSHddbN1UO0y3XmrH6T6SOuN0X+racbpNZu02w0w1i6GvdO1A3R/26JrKFEzhbjDx7W+qgCVy3P6m6nQFmw6tLozR41aLp6mKecA1zlN1aYxWPfCA05/HUHHUJfDSMoZicDxDmqqrY0gmWiOhTNXbMQCfRqJM1csYglFmwlS9G4M3Ck95qt6PwRnc70QkIYzBGiJuyU3Vx9FXEyuxsKn6NPpiUgsRVp9Hnw0ucqaHMvpocIAzPVTCJikZyWmqwY7em5SN5DrVIKN3ePP4xjW40avLJyVONfjRVYMvXPJUQxhd0SBBJawRCVvlSLRwqkTSksnZSKWHTIQgVyOVHnBWTUGSZaqhjs6Z0uy6qnE0o65fnirgWE2JxKSmCjQWU5Jx4qcaiR6YUoyTONUYRkmmEr0qU43qEBIIczJVcDGa6jrFltfwhOOd2HeEXBqlbnNNjUrYVabmTiEypuLnE6ohZkAswkmdavJQQuBCSWVGpVbJ2J76/jQ1sgUuasqNZJBlqqkwLDLBjSQgawhNKJkJfRpCpU7yVLMoSThMSUI+hpBQI3WGRFwniRIboheEwCvIWIxorAGS6EKLFyqZG0m/PMWZN0KUg9kGGSFhI061WCWJFCgpjYQNnhK71qBHI30jYSNMlYAcTxk+JAFZI47hQ6ZGwgZkbiRs+KkCkFSGDciqpIcNP9XauCL2JG6qOFGF6GMn21wRbxKZKihJdBo2INsKetiAbCvoYYMYnAqUeNiAVInS6KqSKlISkO86VVz6YIQoliWWRtAOGoEv0ASjAIyZJqpuJBCO0Od+FCLQyTibqU8Ek4j0TMdGl5kmLGMkVCLi1M9KR9hRGq0C7BhnI881nCoRfqDFjrWrYiVFSU1j8NSNBNNbrESseCXhRZ8GJWHF8zSOhHTgRMkEmWCEiCMBPOKJM1mUhA19WkeiecwKIUZnta1GByGFjoCzypw6B5OkPypTpqauhlulMmXYuq7TjAlrlxgk6ATbGsTrNIyqudvSHBVeCBnBOfilzDTxV4ykexywStrGSKF/j8XSafrHCSO1pFkJaCKMLaZY4McTi+00/HicwxYXJT5pvWYENE7aIrkZ4BXCs4XZ8AVYJnEDgosF3ukfGn4AePpveC8WiKd/aPoH5OkfDUD/waMsWhyT+QgBVdJo+AkRRSPgeab/kFBDUphz+g8ZJdVo7b+gwlq2ivkIFQXXaOYjWtSfpt8S/ERBOWpYNNF/dCjSRtN/t4k1x8V6xDBG2qPQ6R+7OBspOA70H9OIjlRa+8tjgh+ea39lRHaUZvyxjqS4iPtqqDxZIoiNpr8kIzFrngf6S5owbDT9JT8yN2ToAv2lMJJaVJr+UhzpS2n6S0kD3UozvpRHjdpjR2h/ZST5B63h+FTHSvYPmv6yJavUaOYXtWBhGLyjQxQDDo4W0KPazEyRIiIFYUSNtAL6RD2waUrVNCHnZ0aySVrgtFeAhW6BYO2WtADdKp7ThLwde9Rqsk0sisKRWqSAlUVVILpaQLcoC804UgPZQl0gjGQ+ifNbFAbSSA2hW1QG4ncoICZP1pIazD9qAwHTAhhDcYTOB9kCi+pAxLQGjKE8sM60gAmaY3BaAKcahOvzIXCq8ZA+H4r0qkSYUzIWsK5qpM+pYjmKRCWPGnCKKgGWtABOUSYNZ0g1UoBB0NeWAk0CKRIUI8CspoHY2lpA1mXOBGmBZmjIO8zvIJWCTtFERDHgKdkglWZociuaemBoxQComhBiIkomUSqaENIEejYALBkhqchbNgAsOSGOzykNi5pwZ4dkA8C2xBAMzXSLVCkNw5qGYEiYeNB4LvCTNGNDaoisJfUBYHJDpISVhj+SQyrMhSSjkB0iycJ4AFwSRE5lu5gCf8RSQHzSufRPdK+LOoCrKSIQvloD4JIl0pwzNP0T3QPxq5jMfBFP0Y0gBsAlS+RAeGjtnwwjMu0MgCvE90B4aOaHoIpuE28AXEHBkNqq3rREmhud7ppgAFxRBYPABwPgCgpGNxFGP3QcnUf+cQmg00gOU2n6R8HolkoGwBUUTKC/RGJWxLezC6wpgCvBjk43WDYArpCSB8Gh6Z+cvO63YgBgwY/Q7VYMACx4EkG9AxPpH19Cd181ALDgTaiBoi6mCP6ExhmsicwHHoXSzgDIgk+htDeB/vEqlI4GQJaIj4ppkQyALJEoqRpgBkCWSPQWS6SaAD+RnC2Go2hqVSL8QXsDQAspe6WDAaAlwh/PkwGgJcIfdCaZLRLhD7oaDz8J/vBWrAGvJcEftDPAtST4gw7Gw0+CP+hoAGtJ8AedjWc+Ups/ZIgkpST4s8ZhGsNPgj9oMSC3JPiD9sYpP/AHTdqcVCv8QScDbAshGaXxhqDhj+dY4tCsbzXIDCAuua0vZ10EfnLsz9mD0KnTWCTQ+Iu0x+qCRv6gSe5ykOj+h3POkv1Gz+ktHl1sljvvFmZxoafHYjSLHxbjH6oLpjpcv2SqA1uKqa6a6q2ppHi9M9V7U30w1UdTfTLVZ1N9MdVXU4M1NYipwZkavKkhmBpwHjkjArYVU0M1NVpToxjO4NboTY3BVHLgMZkawcRiaqymJmtqElOTMzV5UxMHTaKpKZmawNJiaqqmZmtqFlOzMzV7U3MwNeOKJlNzNjWDw9XUYk0tYmpxphZvagmmFpzUZCpYqvhdTa3WVHLQ1ZlavakccKnR1IoLm02t4L7m1YnDgMRWIyGWDIRGRi25GCDYYjZYwNdyYMcCuxastQCsxa62QKsFT616xSCpBT4tVrIFOC12ggUiLZawBRwtiGgxDyxYaAFAC+pZjFkL3llAzmK7WuDN6sElgMxiEVggzIJbFvPTglgWkbNYmxZhs6CSRcxsO7xDC6TKEpixYI6NtNCoHEalRcdbwi8Wm9Ki4S3WpEW1W6ItFp1uibRYPBKLRrcEVyyq3GJKWnS4xUGwKG9LCMViQ1pUt8VSt+hsiw1uOUBgsSAtJrVFVVsMSIuOtpiOFtvTYjBarECrql2NBtSyxViyKGSL3WhRxbbSArPRooktBqNlzYlaYDTwD4ccWPN2qok1JxCBvucfPH7WXA1Fogwobv6hBWuu9qGw5u1gAGuu56UIDaBV+UcPUtCCNcfvRz3yD32w5moQthNZrLmag8Ka46yju/iHFqw5Tjj6iEAGfbDmavzhUKNY+IcWrLke2FI7D/cYtcA/tGDNNRgorLmej9MTL8Ka6wE33FiAmn9owZrrCTv1QtvJMNZczTc9Iad2mzqFarC1IySsuZpr6p2pnaZ+l+Tyx59//tn8s05zxl86zdnuXvzCOc5+YpPT29srKs4HnlvJyf4Fdy59QB5d7Lh7p+Nw20W4Y/LpdZcny+371fbo+kwrOHplL9hfCXmx+ri7d/WW09gcTYZsD+1gmQN9frF+e7UgKtLoo/fz+NGGQ+5JTxgvP64/PTJ+b7e7RznHt8/WP6xv1purm8UYhR55cvTCx8vXq/lKD/0p3XoI9KD00/Pzm5XetQFne+Geba98r9+8f7y6estNITtYDjnrGsxNdSycQr/bbHcxHzffV5m752Dyq1/NCHU5/44R/vuvZoT7Bfob1/D+cnt0U+r+cjsLBVOmssg2vXj2vO2Jh9vlj+02RqOfXu8ONz8a0S9/NKLf/3h6vXvYztu323IcgmcT6TZ6er3rkMAgnl7vHuntsl710brfC/jkBD0VtPBsveOK3Ey/2Gwu9AA9Be3myoPN1W7zYXvT7y7c23V27iDmvd2OTawg9QtY4P5GMGCvMPx+Y6lT3IgY9GLeV1dnX223m345jZ2tpFanq0cfrt50WOAh5BGKQfYl5Ck3YXplxg/ZK7PvIY/W+/Hq7erq7Pi6Ddy10iOE5UWHwrnv+VYhr5hhYV+R+euiaBbnX3MXYnVzB8N76fPr5RvuAmjf+6t+R2PY3/PrZfC4r3ebm33VuXiueqdrrXd33EeFh7tGX69vEMhjfijifZ2dZBn2XG/uuM3OXHUu7RXvcEOtJ+ur9eWHy39fbTeHqx48uHX7UiG+3Xp5tl2dr7b/9vhQu5UfTVwrOB4mnB6XHsbZSh+uzr9ejNGyaPuS7xdjvl3yckEe4KjKq17wbHksfM+Wt2SLzvdFh5616NNbp8+WZ7eGztw9W559eov12fLsMxdZny3PEPaXh+npJa9ulaAb+10lOly/ed9vKj1bXrfroi87aOwLXi1ITy7On7/ZrlZXj5ZvFH1gD1g7mn5I9sGR2FJ0vB5zq6P9QxXIw+6hZBagJlatzvYSFF2kIehFGgoVMyrXFKEUYnBOG/n9fN+JZ19DiI5jt11fP1y9WV8uL272l4oUkru14/Ymw9HotMKd4WnZ8fgQEi08GuCe3quqZqRggLTbzftWewjUNm04dV+b8czvZzgNfDa/u7pabb9leNRkq+lrbxbjH7hIdHLCP+JOBKf/JJyUmU4n4UQSlLNUCL3KgdKi43Jnb1ffPyu2Pyv2JNiTcBK0pvb9T/rnj4DLanm22qKt9c6UTtueerTePZqFJnah0ZtYrOL+icqSztqb5YU2ZvX/12Z9ReFsCDxYXh+TL9aXe3Myl1KlBEWM310u36540R7gHyyvzi5W379b37xfbb9dXr3t97Fb+f3Nx17WVq+VKidHtzd/v95crK/m0n59sVV9sN6+ubiL9v0R105h+kgBvsTk/urj9ctjs2cufHVc+OpzNefCWzWp+GT58eH6rX4vACF8ut292zxYXq62y44+/0R3TTFl/ojBXRNHddyX3DVwgwk6AizIGXv606Ppu6MYP3XC9EI9twH7DfpPrilzg16fOzbspxeUj8EEBhrdlcd/1yT2q7q/NI3d6313vd2cy5lalADc0T3fzfnJ9Z9Oflhtbz7cnFx/pMrR3Dej4j85+3+7C/zNm9UFSiChGL7gwXb39Uv3NvvIP6rrNMNzx6NftcuLN/F6jbeekNiXh2v4L/dX7V+2Ctz+Pn+03t503+Xxcv7FN0QQdZ2Yy9XD9c31xfLoqjpAuodilQmCAIer4082Z4+Xrzv9C074X7dIP/1LL9Jeue8NhR7O+IW4xGGRQIu2RrynLZGW/Zcv0ReiCH/dEv3pX3qJmFAFtl/5Et0JfBxFRazaMV9d7bZrTEgXW8GLmw+XPy7ILOmnXZR0d2g+vxGHakssPol3pIcWrepHrZv1ywlVYnUhgilPlh9xAvkwSYuNdK/wUPDNZns5uxzIdY9stA9OnD//cAkXt5NbiayXJrdAdXJbfvDVuVyCSzUSWtc788YOEqTmUEtw2edaTBi8cyFn60pJxK7DEJx3RXIo4nLw1EgxJf3QTY6asxhiIhOQfKh6gM4OJPWKz7XaSGIhDzbkWJI4G0PNxstQa4k+22ir8ybUwcboQnYuWZ+riXnwPpZsU8o5cvo0DCln8cGFxAdsjHg/OEklSHVC6kpyHJK3UqvjoyicbxscXQYqRdKLg4TgbCXnkpPxwQ4pkyXy3uYoJogbkvN8MiWnlIIJ3vEFoJxLriXobVc3WF9TSTmH6MSkVAYbLN+KQU6qybYOudqQvORSrJjs45Cqr1ZnNTpTYh7E2RqT9UHzcTEPNVofoq1BMyPWy5Csy7n47DSbUvKQQqqejjU1EQYbUwrJWdcTL0NONoZAEozbjVUGK9G7nG0ib+ViGirfFarJkq0U5/0QfSQLFUj5uJKGElN0vJZ8j7gwWFI/LmRNSCUZvC0pI0iSTbVuyDY5nyU6EpClxsH5ZEvwKZKtZqjBOxtr4pNDxaTkhhy9OF/FkVqJGa5yIJuUSAjFlAYrPgYEJmZnfKlDiZ4JKzFmUtpxsMHVUiUUvYUsfsheXK5OssvBOMc7UrLFJ5LuzqYhB1dt8q4W8l4xD9kHK8G5lEilyhAjW6Om5MnhitTBpWhddLWQkClh0NdHm8VZZ5IfMtPpixfYMyEP0ecoMaYQSfU6P9hakbSQUtHDr0ONzuufzSmYNBSXJXhbXMxRM3ODz9ZmEnGRg5CDeGc9j3MoyZky1FI9e7Gk5IMzfvAl1mKLjxGRMiQurLWpSLHW5uz5NoA1YXASo9TsYqkiXj9Rc+sjZgQV+abZpwmMX/rmzZ1w6rPlDytCDV/80hmJmyObF2Xyn7R48TfcwGCtddalqJ9swvmwgy+lqKCQUwtq4OJ8zLUTaKOfXWu+SC2x+hySy9V7/Y7ZS2nffRpy60DYO3wHBPemP0qVCWeuQ0lW7UTXH9XeqtrsQ6H/V/tHvrVyvgQOJs1+TtOqzMvvroiD8+v5u+XZ5scjF+3BZkvAY3m25rsyKIX9Ur7e0jc2UfviEJ9e09XA4/98qgsNxcOe58Kguhvk2ee5Dlmobt8/3ly9XWHqCq7k4aNxLTO1vvobv6F0R54eL3cr9ZP+gutx7HUcO2BfHPNg0SG3Bj7Pw3HAaz9uHvagYY8W/tVSzHw+5oNSTYm3BSbA8nS7fru+OphRKfPttyqollg/Se7+PXPRjZiT356oFaNxqL9aEv7bJ8QW8XzbLDufY/yHTMiT1fLq5OS3J3znTYrPKjhflJJ/MQH5Z83HTye/PdEv34kPBEP/R0vI893ZycPVDydNSmoN/v/PSJsRxOREBmuJov9qZWT/lT79Mt+X/pnto30u+un17RT25yLrWD77SM6co3m23K44IfLnxW+2q/PFKPHn2zxoV+0zfV0XtBOL+tk+DSvdX1/x6q+226d8X5Dn0E9/WG3PLzY/4qNylGS7Jaz0lzxBGTA4xGOPRw5FNU/w1A7BZ8lSOWCVY40UFTzBoD6GfmZmkFBLytVlKd5lc+qGaq13KUotNiShkctRKlZ0TDbzRbUheO9KqVWK4yTZqR8kxxQ4HudDSN6koXonHgdPXImB10jGE8UUluKqN6cypFC8+kg116SvqaHGGHOJzuVkTp0fks/e5xI83q45lTL4HKrn+FaKrphTV4caqw3FFcu5u9M61FSLTdmVWLH1XR2sDdV76zLHvOyQQvTR5eiLK3pDyQ21+OSLd47jmIIrhdWdcrScez11MlgfcvIi1YfMVLmBb3HmXH3xFkbCkErxkmKttko0p2KHzE2n4F3FkzCnuGQ4dV7wiwOfKBr4Iq1NrqYaYB/f2aXA5xiLS9UbHE481phdrZ7bVmGoVrLF9fMcezx1CUfDVc93ZoMvJg+24LCXWqRIMKchDTH4EFzOPoZg8lADxkktuQZO1sY6RFds8nyD1OWIAxRcrhykzDZXZ059HaCstTVVx5FEP6SIp5dsKRyFdXhuyUkKNWB0m1NnhyQ+2Zg9LHM7dbA4k6liYevVhaGGmIPLySafxPhBaomppBCt43yns0MIEqxkJyyQOXVxID2FWR9zslXlIklKOTLhkphpJ9VK8tVayznB0zh4/DNf8N5LziYMKaVaMktCrAPeaomhSogpepvMqYQhZHHFI5fZ8RL1nwsrVBIXJIdoc2DixeWMQJbBOWdzzr7YlHjtEDKKJ9ocLSd9Uc6Voblig0Tmiu1QrJVI1Zh4sx8sH3kNWTIxKjF2qDG5bIvzLtcgjkbeO2tDdKXEUPTVzvnknA1EdwKfE+PW+FCCC4QXLN8jzg2k1lf7qNhfhJmOK854Ew3HIg0nPfWbYVxQrnwYSs9p6gUW3GVnXE3GZ07Nc9nZmiTJZLGmSDRFz+zqGWA9OKkHNjn1ygfAxHGQ1XEcVm9uOD2IqQeBOQSdqimpmOKcyVJM1IsEzgQmjf2RuZrZ7+LoGV+98Zsct6N1mWFcxFSTTDTetME5nRQ9uHM42vBKcwRA9Kv2CV5UJZHDnxbjqUf+iwvFRhuqpuLbM2KDrDZbitO6Wb0wcL+F9/7wx3944uv1RXcU29fL1c/8UgoRtx4rGb3/OU/sH+jVM1+aQpRBt53VUEdIxFGbU18LITkfpfgaAnzj1H+uevsmt0d1+RRq4ZA0b+/eOR+QrlJjsI6w0LFPD1AUb2NMEovTVrPjrrG8YgNIziH2Y5++RwJCTgTu7vj0jOzv9Onnpeo+/U7zin+DDebghQU8iiEcfNsv+/R6CuRfwaefv5T8mQzqr8V7i1JSCDYD1ta6uzagGoUcMFpf8385oBENvjF80b4xfHW1erM7HGPSb1c3a1G/8vvz/wN3o1Gw').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779221362328', obj, '');\n",
       "});\n",
       "\n",
       "      }\n",
       "      const servers = ['/static/', 'https://root.cern/js/7.11.0/', 'https://jsroot.gsi.de/7.11.0/'],\n",
       "            path = 'build/jsroot';\n",
       "      if (typeof JSROOT !== 'undefined')\n",
       "         execCode(JSROOT);\n",
       "      else if (typeof requirejs !== 'undefined') {\n",
       "         servers.forEach((s,i) => { servers[i] = s + path; });\n",
       "         requirejs.config({ paths: { 'jsroot' : servers } })(['jsroot'],  execCode);\n",
       "      } else {\n",
       "         const config = document.getElementById('jupyter-config-data');\n",
       "         if (config)\n",
       "            servers[0] = (JSON.parse(config.innerHTML || '{}')?.baseUrl || '/') + 'static/';\n",
       "         else\n",
       "            servers.shift();\n",
       "         function loadJsroot() {\n",
       "            return !servers.length ? 0 : import(servers.shift() + path + '.js').catch(loadJsroot).then(() => execCode(JSROOT));\n",
       "         }\n",
       "         loadJsroot();\n",
       "      }\n",
       "   }\n",
       "   process_root_plot_1779221362328();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779221362338\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779221362338() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(58503,'WkwItWYAh+QAeAHkvWuTJMdxpvtXxsr2YzIs3CPCw73K9AEXYqlzAAJGgCKwMNmxBqYH6IOZ7tmeBjiQjP997fHI6ssA5JG0q7OSLckeVmRlVkbGxS+vv+75z4f/5+7n15fXF68uD8fDFx9cXP908eZPl998fn3x+s33N3eH7fDij9dX//3Hy7//8HCs2+HF+1d3b9anT7/5fy+/veP4gdM+fX13dXO9N/7vq+vnh2PbDi/uf+n4z792r792g9a6tt63w4uPr64vP7h5eXN7OMre/Pzu55eXD80/XT2/+341P7p6+XI/mc7SPJ9ca15++eLuk4vb766uD8daOPKHq+++f+fQ+zd3dzevnp72xc3rpwe+fHFFJ3Q7vPjq4eN76yM//OWbu4s77hLBOU9a760WF390e/Hq8t1+c+ydB78/7+kD3Z96Pnz+SS5/PDL85Ps3t88vbz+/+qd99B4d/OTm+eWa1y/lcPzNKHX9Z/ZRh/OrX/2V41/q4fhrp//64S/vbt775s1nV28vX/5wOLbujNTdza8fOBwtyoiI8DD1GdmPJ7+gzejbwy88PXA4/qbPIvxEaJjZYMz/+OQnWCl/fPiBx006wO//w5ML+uRH/uHhkqcH8p6ccX7Qu5svf/j1Mc3H/v/+mlVUpTebLtaq1SFPf/+rH351Dvaf59tfm9Hz14fjb2qpKuLWe41oZj6Zl5c3f/79hx+sZfHV48aXf3ydX+R8PPr8p/ujv7v/9N43b5780HvfvHnyW+998+bhsve+efNw5R/fvmKf/gYx8MefH31+++ri7eG4Hz5//OL7y7uLw7Exf599f7V/eu/N68tv7/5wcXd1s57j9z+++ubydn3+4urbH94+fPx5ffz45rv94Mc33z0c+6f17WcXzz+7uLpmZ2+HFx/c3rx58/3F1f6D983PbnYp+XjLMWer/bDbPrl5fvXi6vL54fji4uWby+3w4r/eXj1/+7T580PzvW/efHBzc/vo/N8+v7q7+AZZdHf7Iz/w0dXby+dPnvv805/dXr26urv66fLNL8Txx1dvkPZnTbA3L25vD8ev/3E73Ly+48NftsOL3769/PbN4Xj948uX2+HF75fq+FaR5F9c3dGPw2e3N8hB/fDZ91dv7m6+u7149ezy7cWr1y8vOe33P7767OLl5d3dWbAzkr+/fHv3y6Mf/v3nn3383leH4+G/nD9uhxcf3vz4zcvL93988eI8k3+4vLu4umYY90H48s3VP13+8c35+6+eNvPbP1xevDwclZvn16stvaiPycr+09X185s/f3Hzmg34uP3V4/YuZB9O+N0lKmVfH38+i5APvj/scuKDi7u7X4z/e3d3S/vyfF++f3n358vL6129PGnlGH90e/Pqi5vXh6MUFtWXzy/ukN7Z+OrcQP2+txryl+3wwyc3P11++vriv/94v1R++MMl4/L04IvfXX33/cc8wq5Jc9le3H37/Xlwf/j8+5s///any+u7z+8u7n58c786f3jvx7sb1sf9mZ9cXv/4/sXtarN63vuW9Xd/xYs/XF48//T65c/nK1786eru+5sf7x4v1fPy/d3Fm33xnY88PuvrdwyM/2VmDHLmr5oxf7r8JsXB1fV3zNxfM5Y+eHnx5s2+VzhvGU+PD7xmqR7qUcfY9r+THOtWt3rSPMqndqz5bT31/E7HOI37a+rJ9jM5Pvdz+exHj01FNu+n4LP3TcVOwg37tv+dRI4SuonrJlNPosdpW/7vJO0oVbf97yT9KDq3/e8k4yijbfvfSewo07f97yTzqLVv+99J/Kga2/53kjhq123/O2nNkyW4vp5Unjb1KNE2CdnE/KTtKG6bOM120r7ua2OTUU86juKNbzYZ46T0au9kryedRxljk66btHpSP0ob+YBuJ42jzLFpNXpyavUoTbe8tbVTk6O0SQ82sXpqepRWNxm+8cOtrfu4r2sZKl33Ge3UxlFkbnl9H6dmR6mdHmSX2zxKjS0vaHpqfhTtOa5a49TiYdTHOPV6lFm3fNyQU5cjU5x9buPU9agqW44PTYbKVydbPfV+zCFm1rSf+jjmmNINiVO3Iz8TfQs79bluw0D1eep+zDtKdvPU4yh9bvkQoqdRsxO5xE5DjtJtyzV7GprjsDfaMfz8uR9nP38exyHnz3bUef48j+380VnY+8/EsW6dBTxOVo91M1+f2Toh67MeeajcVidjB0mL9Q2biNWRlw8aPlfDuIesm1jupLbfJe+em1ROxu2zofM0uX82aj9NOkBDvJ4mPcjGsNM872Hmd9KD/Kb6adIDGj5OMzswxmZymued3OZpnm8v7TTjuHq41ZPX4+ogn+U48jH4rEfPYeAzW3hJl3pyVmUOA421fblfPTnbxM4Ntu7qVT05GzeHgQbbNkd7q6fYBVdzGrIkknQaS3hJIL6i5TdiRqOvRlcaS4qxqOoplhiL/OWZZ03huOfnPvkc+VnbVk9S61FEtj5SKEiVXKFdt64nqbqEQjC9cpLajj43yVXD1z2fQ4ZtjeY4qvrW5jb6SaodZdRtskHjJBWRYXmuGD/F9tSttU2VJkOyGaMqJ0nZOthDehJBjMnWZFNamqOKMLI4CbKVz/zpOAnCldlk4PldGavHaz0L0rXKFrFJrSeReUTA8TziJ0G6su4sNkGix1E8uLGcROtROjKrb2LId+TYGjenqce8rM7N4yTajtryPvUkiFYfqIWNbS7IVmlbijLlZMutruKrz0jXRypAqxxD2JIhJ2VKEInr76QV2dS2/e+kzInwPPl30jqWvEWsn5RJ8ZYioJ60ziN7c8o2+WU/aovVLZknrXHkKoT1OClTwirg4pOKHKds6+KTiqZc3q8+KZPC/bj4pNK5cl18UhlcuS4+KfOBRJfJ1SeVuXp7vqnfP1zeNbK/59ui7hBHq9MnzSmht+u+TMl6VvF20pySfFht/aTa758278yk7E9LpxUp/jCqOs/Puy72+wfm4nbew8i22k/tvI1ZTOgtdvKuxVd7beb8vuqpsZ/P32f7kU2Q7bWrOX/IqbGx99Oz6fcWRDaRbcvcGHLq9cEwqqdel03A5bT0/kpaD3YKrY46QaOcemU0lvTi+Xpu5pT9+TydzUwbbTzaqbN0qt4/X5e62nxf9dQxlugiGjTbyr32x+vSslV79l+yH5zK00gqOS7O1lnU11OXFPXruQRRzyUcR9FwQT11Rc+sR9J9SpbNeOqYQ/sjodI6m3Z/pHxk1sf+SOt7JInmI6/vsT72R2II1M+PtE6P8zPl2a3S5AZ5cpNsrlk79XZWenlpS6235vvUWyre+wtT895fx3DkpJ0GS3JfI/zKCLmf6tV+mHp6NOJh8tf3fT3PesbTYEmu+2SXRxhdXrN2GjHPrfytfPb1gKcR+ei0qp6sMgWpFU9Wl5XONVZ55py0k1WeOCftZHVZ61gMrMK1w7KLhvhaO2xvez7yPm0nQ6fw/Zq2k6UIyxW4bskqXCswH8lQLbnD1s9Jy5/jEekgmiUX4X7yGo/zY8lj98F2qx2px5Vrd+wPh1ZhNa7h0FwIzNrJWJLn5+M6JNb+fKwSQ2StHba3d8W/T7+lzFrPl0OdVnrusL291NtaIKZxfjzujZW+ni5ba3/uZ7YHJ8raGhQsEsM6Z4/l5zUc63MuDT7OsyBc3T7NXRDuw36arLq9i0NOk0W3N7l43D9vfvkw/Tze3M2b7MJp7gYOd8qT19OtL/0sA1M6ODJwH9chJ69LEuSknHwXgikkvN7POd/s3Vmyx/fVuPf35OfVmJLGWYv72st73A82NuaSyOsW+yLkZE5Ele5ihTNTEu4/+SAI+WYXR+txHT16FshVT44eXcs+T/b7Vt7jXhZh5OYC5PY0du2wRufkqNCzXOdnEZKstiWMHRm5N7l4t4T3Lu2rcR+Fk++r8XyntVdp0SV8xv1GGNu77FqdYjHuz8a2cFzGfYjYFt7WIK11c/IHAUmPW+oL+rTO3aX2rqi9nf2DdfKDzM4bpYtyvrQ/9lG8p8Tmd/J3e0rstflO3tNT2tdJPwtsnrM/dIfe9YfurJ+57w4P1u+7s04+d4cfGufeZGMfoaWiT467uM9LPsi4X0b5u4Me7dLZl9e4D+jJx+rRPuAnH/QovwUKGdmjhyY92uWjpw+56wO3pdDWtJ18+ZFr0k6ejuQuOH05kqu3J7el3tdSOLk9UrLcP93Jh19d43N+6vQo7y9dLuX5u7lDIecOzV37nbs018I+D8pcC3tJbJ9rYZ8fdC5j6H4UJp7N2cP2uRb2eSBwNNeXrChfC3upP/c1bedTfc3aUhruOWdLS7inDbR3B48z19aaXHzOXWbmGPna+ud7rK1/f4/c+vstYtc9edkuqc8jsAvq/f67cZBqynchvTeWUDzP5i6j7+dzl9Hnn7lXQLmgdxmdPxS7iN4bDxKacYtHIpoVHAjpXRdXPcUuptddYpfT+5PEWU7vUx27nF6DGbvNcG4ti2Gf6dgF9f7ljq7lmqJTuBsM/Po7BcKSdbz+TqE5g0uHhvbjaLjV0rg0l3nHNZ6nUDuOmh54x+mfxx446tL5UT9233A8u51C49htG3WT7qdo9dgRn5sMOUWTY+9bdqafoumxty3F0zxFa8euG+63gST0Y68biJvpKdo4tthGgIWdotmx+WYLIow2j21uuMiTO/ixjQ0HeHKHADYx22TaKXo9trbZ3GTGKbocm+LN4xtH12NLl098nKK3o8aGL+zzFL0f1RMkCGCNAWw1B2jhKUDSbJtzk+AOE4RgxibBHXBWN2clyyl6HFU3X3ZdJI62pes3T4FwjM0HmNQpEI2+uW0q7RQD9GBz31TGKUY/im0BeuWnGOkQAoSpnAK5OLbQvcWWT3hC+U3sOyCX1Uq3OWy1DLtqi7m3WDJb4OcD1YAZgEWoxCms0RKAi2xmZ3LVZnOsb9v+ra3mAi7C5mrykH4Kcx5LKk9JE0B2A5rI5gT63IBKVeYppmQTOCybQD4bkNBq5giJ6N4EJd5AL4DAA8nomyTWQBN0YeGF2ZyryX35Fmd+E1AORhvJSJNujFN4zSZIQTZlNekG34JdJ+ixmm016UY/BYAc3/L4NAFkN1Een6atJt2gOVeTbrRTICA5mW7QjGw2utFOEatXYE+ip8CJctDHvbnGCrxJ5BRISdBpukFzzWCjGzTXDDa6AQaXC0oa3aCZKyrR1WzmkpLO+o5T4NL3TUCxKlgaoB1tFrzTBoxCYJzboOqbdOCI/L4dBQTaNq2T80EwQaTP7bHafm4Dy2zSA0Sc82e2B93JNloFsbNpHXyfcKoM+kNb6jF2VZxNyWaGMfhWj4DpCyuRKi2b9CW/7dmkK41vxxFIh55k02gaHQFxBMADTzw3PZt0I7+NI2geowLEqDWvTXSQpnAjxFkwpqp0kvBHMGRp6ibcKsGQYevq3uaZsHbBIJFOdDtBvL1NR9PcXWGOoC9ARvQc+ZWdWcs/ZSS3xwELwjabOPdvWCx7m/vjhBFayqgEbRDGhSk6/WlgsXub/jScw4WLgk/WlhGBxEkXkjsRvAI864xGc8QygRskuFTEO/enTX8Q8Nx/yXupiHjuT5v7I+S5PxqA+/eGslg4JuPRO6pktelPHygaQZ5P7t8NNSTOmHP/PlFSq533d1TYilYxHj1QcKvNeIyK+svwm9GfISjHhEWN+w9Fka42999t4oxxMR+jHwfXo9C5P3bx3MRxHLj/sCM6Mtt5v3k0+sP3eT8/snayzfOPOBLiAvdNqNwqCOJqcz+TI5g133fuZxkwXG3uZ+3I2BCh69zP+pHQYra5n40j98o29zNLoDvbPJ/NY6L22BF5Pz8S/KOdcLzFMYj+0eZ+sxJVWm3GF7VQ6TDyjhuiGHBw8gB3TJuZIUqJyIF+RI2sA9wT9cCm8cgwoVQUBCgrBzTvimDhthzI2xIW4LYcyNsSt2OP1gy2SUVRKKFFDjCzqAqWbh7gtiiLjDhyBmsLdcFiJPIJzl9RGKxGzhBui8pg+T0cAJMnaskZjD9qgwWWB+gYiqPv/SBaUFEdLLE8g46hPLDO8gADdMbg8gA9TRBuHw+hp4mH7OORkj6VCGNKxIKupxrZxzRlOYokVx5n0FNUCWIpD9BTlMmSM4QaOYBBsM8tBzIIlJLAN0HMZhiIrZ0HiLqcI0F5ICM0xB3Ov0EoBZ2SgQjfkKdEg3I10ya2kqEHHs03BGoGhBgInwRKJQNCGUCfGwKWiJAE621uCFhiQhI81twQsISFhDgO7YwSAePQoXN7IVXZpsMZhuCRMPFo47nQH8uIDaEhopacjwAmNkRIONv0j+BQLmYnyChEhwiy8DwIXAJEmmvbN6d/YClIfMK53B90b1/qCNwMESHho24IXKJEGXOmzf1B95D4IdtkvMBTciPIhsAlSqRIeNp5fyKMrGndELgCvoeEp834AKrkNmkbAldQMIS2om0rkKZHzV3TNwSupIJhwfcNgSsomNxEGP20RxIGJXAJaNuRGGa2uT8KJreUbQhcQcF07mcEZkXa4i4wpwhc6fWoucHmhsAVQvJIcNrcn5h87jffEMCCH5HbzTcEsOBJ9PQOtsH98SVy98WGABa8iTRQ0sUUwZ9InKFug/HAo8i2bghkwafIdts698eryPbYEMgy8FExLWxDIMsAJU0DbEMgywC9xRKJrdOfQcwWw1EytCqD/tFuGwJaCNlnu28IaBn0j+9tQ0DLoH+0J8FskUH/aMfW6I/RP7yVuiGvxegfbd0Q12L0j3bfGv0x+kd7bAhrMfpHe26N8bA1fqwhgpRi9K9uimlMf4z+0ZYNyS1G/2i3TbM/9I82YXNCrfSPtm2IbQGSyTbeEG36x/dY4rSZ39hYMwhxmWt+4boI/Zlj/549SNv2NhYJbfxFrsfqos36o01wFyLR+z++gEv2X5Kpd/jo5c3FXYOL9zLZY2Nsh58Ox69D+xaK62dbKLLFt9DYotUtCPE23aK1LVrfoo0tmm3RkFG+RYstet2iyxZdt+hti9636DiPcESQbb5Fjy1G3WLIFkO3GG2L0bcgBj5si4FM9C1GbGF1C5MtTLewtoVBNBlbmG1hyFLfwmKLWbeYssXULWbbYvYtJq6obTHnFhM5HFt43cJlC9ctvG3hfQvHSbUtkKUpv2OLqFsQgw7dItoWEFxibBG4sHOLQO5nXB0cBklcEwmpRCASGa3EYhDBFbOhInwrhJ2K2K3I2oqArdjVFdFakac1vWIkaUV8VqzkiuCs2AkVEVmxhCvCsSIRK+ZBRRZWBGBF6lWM2Yq8qwi5iu1aEW81iUuNK7AIKiKsIrcq5mdFYlWWXMXarCy2ilSqLLO6yDtcwaqqADMVmVMHVyQqh1FZ0fEV+KViU1Y0fMWarKj2CtpS0ekVpKXikVQ0egVcqajyiilZ0eEVB6GivCsQSsWGrKjuiqVe0dkVG7xCIKhYkBWTuqKqKwZkRUdXTMeK7VkxGCtWYE3VnkYDarliLFUUcsVurKjiGlyB2VjRxBWDsTLnoBYYDfwDyYE5X6wm5hwgAn3PP3j8zHkaiqAMkhbi4kMw52kfCnO+iAHMefKlgAbQqvyTRAruwZzj96Me+Yd7MOdpEC5GFnOe5qAw5zjr6C7+4QrmHCccfQSQwT2Y8zT+cKhRLPzDFcx5ErbSzsM9Ri3wD1cw5wkGCnOe/LhkvAhzngQ33FgENf9wBXOeDLv0QhczjDlP8y0Zcmm3pVOYBtuikDDnaa6ld5Z2WvpdMv0f//KXv2z/XmzO8bfYnCuJ4m/wOHfGJvzt22tOPFOe15Fn9z/wTvYGzUcZGu8mZzykrQjJIr/MW/nk4vaHy9tHeTDrwKOf3A/c53Z8cfn27r3r7+BjQ02mub6spTIG+f3Lq++uD6Aiq/3o9/n6oxto7pYM44u3V78kjb93d/cex6FvP7/66erN1c31m8NxCHfkm0c/+PHFN5fn3Bzul+11h84dsv3pixdvLjNpBjm7H7zvdst+X337w8eX19+R8lNLheScc3C+NJ8FHvq7l929PNPN70853x5i8lf/aZ4wp/Pf8IT/7T/NE95P0L9yDt+/uH2U8vT+xe15UTBkuRbZpi8/+3ztiQ9vL/688jFW+9PXdw+5H6uxp3+sxp4B8unruw8X336lvUGCZxPlNvr09d0uEniIT1/ffZRpYvupH13teQG/YNBzQh58fnVHrtu5/cXNzcsk0HNg5a58cHN9d/Pj7Zs9d+G9u70770jM9+7u2MQppP6GLNB/pTBgr/D4e1LS3iIjomSG3W+vn//29vZmzzJjZ2czT+dWH/14/e0uFviS5iMpRnOfQr4lF2Y/meenuZ/Mvqf5aL4/vvzu8vr544QbereOPpKw/NDDwfO9z+mB/MRZLNyfyPjtS3E7vPgduRCXb96R4fvRz19ffEsuQN77Pmfv0TPcJ+ztx+jj/XlPe3N/6vnw+dR3bp3nvfvcjw4+ZBv97uoNC/JxfzjE7+3dscpjn88733iNzvnU89H9xHd6w1mfXF1fvfrx1X+7vL15SPXgiydplCniV9bLZ7eXLy5v/+vHD2ev448Gbh14/Jj09PHRh+dcRz+8fPG7w3FUJu3+yJ8Ox/n0yJcH4gCPTvlqP/DZxePF99nFk7XFze8PPdw5D/0yffSzi+dPHp2x++zi+S/TUT+7eP4rGamfXTxnsX/5MDz7ka+eHEE37rlK3PDq2x/2TKXPLl6vvM8vd6Fxf+CrA+HJw4vPv729vLz+6OLblD50D7H2aPhpsg8eLVsOPZ6P81WP9g+n0HzYPRw5L6C1rNY5t6+QogcrPRNpOJgyI8gQpZUiBud0Nf90znfiu9/RkHyOu9ur1x9efnv16uLlm/ukohTJu7Wj9ybDo6fLE955vDz2+PlYJHnw0QPet+9V1TJSMEBWmvL9VfciMK9ZjxP3Z/M859/ncZbwufn76+vL2z/weJzJVsuffXM4fk0i0bNn/CP6THD6n/Vnfm7bs/5MjJZWTuj7KQ+tPPT4uNanp99/53X/zuuzXp/1Zz3PzHv/O/3zjwiXy4vnl7do68yZymG7b310dffRedGMfdFkJhazeP9NrqUctW8vXubFzP7/dXN1zcGzIfDBxevHzS+uXt2bk9M9xHtKjL9/dfHdJT90L+A/uLh+/vLyT99fvfnh8vYPF9ff7YnV6/j7N2/3Y2v21tHsyaP8zX+4unl5dX0+uqcvrlM/uLr99uW70n7/isRTOv1IAX6Jyf3bt6+/fGz2nA9+9fjgV7925vngkzM58ZOLtx9efZeJ/yzCT2/vvr/54OLV5e3FLn3+Hd21lCnnagTvmjip4/6au4bcYIAeCSyaZ9mzf/to+N5RjL90wjIznmzATIXnA7nveWD//8fSgjus9q4d/neN0n027t8aqd2x/f717c0LfZ5GIzLsIZn32c2LZ6//6dlPl7dvfnzz7PXbZxfXz5+9/pkzH43yMh/+J8f5X+/s/v7by5eIe5mW0/Lrns7uqv61HM19CN7mxWdRvMue/9TuLZ7DN1d45h218uVDJv2X99nzX64TyPV+8dHV7ZvdT/n44vyJwh961lGvLj+8evP65cWjxHSE5r3YzVWBw/+QKP7JzfOPL77Z23/D4f6XTdLP/6En6V6R3xsFO3TxNzCI/5iT9Fcwg3/ZJP3Tf+hJQj6njPs3TRJXr420thRVKPLY/+/76B2Y4xEGUtNq+e313e0VBqOOdeCLNz+++vOBONKyPmnqO23Ka3jpMkyGUBXE02rm1Ld57qxemkY3iZFlWT65eIvLdzj+ZkdCdh/w4cDvb25fnR0MBm3HMd4NXxlxrRW+qtuBGgcrhlWLSu/uXby3aNO3Wrp6TPNmPn0O22qx3roNHaOJA5PX4mPGaDbbHJO6KVLq1N61hkoYoYKibtHqIHRhoO9lVOnioVHnIIhRpoyos/U5ZBAcKNFqizFtxCDxsrzTNS1tRJtqVb01sU3LmM2NeI/ODLGW6d569BZTiH+3Ujl5+nAZQP+tqMymptyFCEornafqFqM2ANJWrFurjYCAEZtrxa3OkNEkCIX0UufQGVzQXTwLeHy+5vtp0PBh1JcOfX2ASZPBQxnFZ3fXkDE62brv/Fd6ac7jznCzIf+4HV7L4Ui44fCTHI5fi5ToQ+aoTbUJNHORMlzFaxsSvcb+k9JKjyZzeDQlKNNLHdqwwUWVkFMrzdqQLmOYts699JAsjsNPyr1aMXMb6lpnM/LKRcqMHibdhkyiS1qaMTU2hnSiQFI8pkSdajrIO96kljn6sDpqa04a+ur0mK1W6XOo02mpJXq32WX4JIyUZ/VRq07vtRlhGS1alWXg6r0ny34WGWpjulWLpjxGOyT34fBT4zFq0cESVdWhRqq8UAtJVdxMWifqWWpvw1kufRDXidKdsRKPEVPz2aNYaBtdR69G6LSMMdRYMN6I70ot1oy1p9X72J+AwN8ICY9JivV6/FZyMTef7p0NuA7PQlkiGd5D3CbP0oEV2nb4qfMsWiQq86rTM7TO+HvrHm4+MyBa+lA1YpOdEFzdojThozVht7dNuU3rNbRKq0yIzKJiPfeums1Ne5nSI6w2q6wzsTKsmTUbPSBa1G0WtrWbqmrPFLYordc2Y3Yf2uy8vqWWGnOMOueIGatAQdFeu3dts82ctEGdob4dfhqH49faSm0yQmcbnQCfSkmBVZEkBKVzFltvJq7dqycbsgzRIX2q9zFj64GsYG1otQz6jRI6qvYxtEGHRTJEtRBtHTnAeHLu1KlzwA3V0vtovVsbbWjbmhVzBOOcY1j1rWvxpoG8UTfzrdVC5DV0TmMet95KZ8vL6FNmb5wR3toMM5dRneE27ZAqmCKCmrNoiPWYMuDmcKC3qjKr9EqgHsprkWAzVA2EKHtPWnMCsayA3NFGaTBG1dby8cHe0Vll5CRqLZEyjhU6kgBUrNr0EK0xSbSfpUqzLq5MRS6F2dqstXr1KW1rwiLsrYcMJVLaRhnexapVM51ta72E8bgzhsHQaUxMcxVtMfrcei8RtVmv033GZOo8omqN4ZMharU4hewGY4QMH7V4jUEkfvYx4fsXHa7aVKePyHoXtViV6Tqn1oEmYN35nL1Jn3AhZkEz+WQhr7mREmLu2kabRIbZLDUGKr8Nq7Z58SlMlKIhRVMQRO8txMRsrsXaJ7K1IUXYxvNwbI4amMyDl4gurYd39eTdFuQON7UJbUx7kSCgK8PQVKzUSlkxNmDGyVsOhzdVgvtjG1HUZxM1lsfsm0bh932EzImqtV5mbmafvZlabo82m6p7wPVQKbOzNaKHNyd7hB1WvWqvTojZePap9Mxihm2YNigsb3Mgk7emxaui212iw7X2Mmqd1fqU8CyHU6K76kBbQ6Nos7CBY3raEp0fbaOx0CViQngYRSQaMst8ki/R2Q2ztvBqDr2m1D6quzjyKbcHFR+KzxSOc0wzYx78cOzAtD/52g9DpyHohvaRkrCWMZBMs8/ZYDRE6WpaZYSod2E3D3dXaCEyG7UuSjedI3Q42nYbOY51ehjGjm+tlQhkfU29KdsYSFNRHgc9tE0to7miRyTVdTfkuCBVPAkPNguWSquurj0Gl0STkI5W9YhtsnJ1donhqNUtZmkxxELZuuTcjGKOnGmdY7pZKzLqZMl4hSkxatE+YCu0UGoZTSvuzuZutar1zaNU9o83hJmQeF8wv+pss/J/m6U4MHXU+WhjG1pmeMVwsgapD4HY60TfiMFaS7ns4ioDCxSKVpkGlSOg0bRU9FIo2NdV+C3rS6g09aYe0WQp0ZRFFiaKzXHWt4/tLGllsP5rY7l1NAk0T8au9loFC3ZTK+IVi6ZVFjaJaI5IVjdpw9rs26xlzsGwzDYcHluUGlW7YM2jaL0WmaN3mcgGs82stCrNxXRWiDTeCmKim9HdMGa1topwQCZNxt+mmcxoBu0nZokx+4wW5j1ZR6PYFHZC6EhyVC9YslH7bF3G5rM0rBl4JGFwu2apjL5jVs8+kQxRu8OFmj0JocUtRhvdap862JHDhvH8SNzYOkanW8McdqMEVhSXGCNYOH30rXlpiJtaVcWogiPFWrQxMLQYrl7LUO+4HV6hhDUraadV9+kdEzVKzEoNRQYjCTXFag2rHRUgKVkDmgSeTWA3QGRgD7iJ9zQ1mhdv1lNOJkdpMG9Bb21CITUMOzHxWrvaFB4Org+ugNTusc1ZoH7Njj1ssPKkKMYrs6pQeucsEYF9yMqbtrkX9Gmj331CxvNSu4zeWJbJFWq91Cp4G70ntbjOItZGzluDQ6SjjOo1MN1mspykTE3driPZu5K+DgSz4Wx5qVqwI2YKsbRJtdTae2jD3s012easbTCguYxn8THc2FZ6lvqj94oKn2m/TLqhuzFiuaFd6xhVmXEx1HeXpiyGCR9JpQwd6CKsXwskoLfo0WzO6ZaqxbUPE3VtDgkRg3MqAzC941pImSx0nDzuhvSWejiOUXGGKvK7laHWZ3eryXBDFLGkrWGerTU3zau0jmUKPbyMKd4cLQrn0cuQsJAuiiOUBt+oKj2sVhb8NjpWvVe85c6QDsefCdHuLSldrRbM4BG9WqAAvBfrklPUxDEZrHigzqw5BgISEcIdwoZRyOXCNM5m+NJti4pzLUqx0p7cRSmD5WhWfSZTzMtoOlo1FERSzEqrmK2KL5s1gYqoRZ9NHDtfWiuCs9Y7bh3eG3Pv3Zhohd2mHSdsIoIE+qloLYquqC1S+o8io1WWYDSqyCI6dIaPsBYQ5LyWKW2mm8tAst1kjtlFfEwc+KHYYLh/6lAmWyvV5hxzmkr3htkhPuowVSyojjTurl5bm1rNdm0gbWDF1sG2zOVSO7Ja0WCey0UOxxFYv3z6WqJITDyEaW2SeVn67N3nCMN5XKocW6TrmMlUriUUs7CaVyPzVsuc02av1XXA/UbXN/bG6Pjk22xFVJqohLphP2oZSFo1mdMz+7N0HYAJKYAyRbAortAY0xYhM8pERAve5WITloG5aM19JGUZPQCA0a3CCUYS6JQw7W49+aaBKu/eUSfQHHkY0ZEuZ/Iym5SOpGX6BnxfQJ/Jwm/abP9Rbz4UDMCTvjvS38Q0Y4vATQR9GNbR7kM3l1JHpCFWa9SOhlNMcEFONhjTFTO0s9LPukV8zordiRfL9E/BYcJ084Gex98AuTAkF9aCouAG7pWxZw1YCdtnhCUFtvTaMPjTvGVvFK0Oz9UcFwGN3wvjrQgM8ZbLRQ9Ha5PlAv7hZUpISJPmin+FDwqMNPGLJjz/dBbbqNZjNPKyC4aRNQnMI8VebsJ2aV1aIyPai4XhCY7hCKRoRZCLHEoB3qQMHJ5uNjDKBb2JxqhuWAJQhYuZIkis9mRIY18NZ3wDPuocxWPMitgTQ+o1nBjMTWG1s1nYCiBt2Fs50wNlrOmxtgoPvCg/kdAa5qZgVqiMGqFMJMTdomyOOZo7VrpKAeFURBs2W9QybDY0JeTkNtAaAFuxSM24eGw3GWuwpMwBJIJlQpqGlzEg8QJsgPsjvursPuuOQWB5q+FFVOvpSjjO2rDhY+qAZm2li6nibfY0l7DRvHamfkYCD9IOR5s57eBFTEAdPuDjVqj7iv0E5ojtgpSsLEnFtW9RazeURCsDe8GgkvfkIksiGA3piJlsoGAzEr6ZQXq6F5kVVRTC8oZnjneVywuXfAvEHuhO4P3AvW5aRmAMIqKTBM7ecPVILxiAzAqaTPus+NNb9NINiT6rKER7vOQJ7NBHFXTktCJ41S1m1Ea3UtJiolrknrfSOv5BsGah+5fe1XtF6rdkiCuIMk8OCJjM8EKRbSBlNMlab3Pwe60n9bsiFbtMjQrbmW3Q6TH2k2VOWFHkCAhtz8Sk2ksMVbBkxKphgc45sbeB7dh8ivprrQ3V2vD1m6HkcYbGoOBBQY6xhRsaDWugYzlHpLWY6KUoSGFVjbZgv4bJNZjqVuH6P/YlEqIM/FvzBtzbMQbGCMMnn9oz57bUaBUgHXfdcVEkrcKGQ9QHhrCDXovGdCeRQ4upjOGjqXdXnD8BeGSfqSvm5mg+JBy92fvmo2iY19kiZOzmgeiUUYHr8KgEAKWCuOL2wOHvTAEYHhAIG9qBfiTGjOmZrTCLAmU6oHn3FO6zGYBJJMBmVjAucGVmRXSaYUz2ijnQyQ3BWMLLqnPMRjLLbPgFMlhLBjU+pOArg7iNLmSR4K7g24uNSZKDeKmY7risLffWKBHDnZAB/ZJaiwx8tTp0plMUBXx6DLDIzKUgStD6nNHdZQ0Yih6PYIwhDDHBiwiAQaIbfaSo1oipw0AT2dFYoq3G7G1uXoQIiscYE+8K/HP03obZxK4gSb2sQvASVVBdTxcPrURwqaePLUQiQaKxmKi4ytNzt2BeI2pCgcY7gk9taCU+08GipxTwCkySSDtXDKsLWxLbj5zGMqeq515BZ04vKK2EiQ2XDs/SAr06B9uGJIXS8YA7Nmhm67RC6EaRdKSuhJUBfjF56AFsXnE9p1sAqTipeyW0VzAWPH3DgCKaMIHb0kptXtTMWh3Y6KEcGOg6/OxOxk4LFI+2SRIYgnqwG6yLocrTxiBqJFi5mH4rDwOJ5E0rGk7INy7YbKJjAoeloOzmdeCM90y5acBQlRykPkl+wQHp+JLWuRN+NC7vcHWv5DfNXiqIHbB5B93rgao2kDxW3gTFAdHrrRPWmAEaBvBreP/W0IoGGmbSsBqQz1HwOo14ifnQc2jgYeHg8mm4sUXxXRLhB43ijrlixualNWqvOyBclsUoI0ZNW7/1TAMtxL+IpXSsUFa7D1wJY4hTtQy8ccd6VLIanUtyg8QEVE2Tk5HA1VHSvkapxII6eT0pXquXCa4NzgNaiMjRGgTheuYTTUwskNW09nAmepC0lcY1mZ6EAXDmp6nlg8yiNicRMpwv6tKUWVUngZo+yF6aBW0K4tLGMKrv4olrA5PBskhoeQqQJVtqjLSXwHrEo1bysjwKA5PrK51TEQBsbXWiWzJ4NhEGjpZID1j6LLY0IDKKBLii2g39ZGL4ZzrwDGQkAISmlYJLPPALsBLSi/aGiYbMHiReFlfgJgNg3GygqfGAgrShjIHUotWWTZa+EznuhfQnHMeGR5WRQOVlF6CAnrG2tZpYSX00LOngqzwcpeIhuLoMHVmjAXhTpSXYv6IqHrzKROYUi+ymYWQxqsjsVGji3iSijzEaypqXoGBGA4s2HkXbHNFGrQSz0qodprWO4V3N2Fz43EShAz8xoY3aJroDqz3Xm4kNaS1skOqI2yiKnrQgZ9G09GqNIJEGIq7NIg3QcqR9RH2P0tPhYx83UnGxOrq0iskAhNcLVje8RWITmTxaYiq4ciN4ngGhcAGublrBggx0E8AzbBr5kCAs1kfXOgS3k21AuBHn3UfmKaJ8EEyGj49b5cV69RA8XmoSzQT5BrDhJOc38EUIUM5uJKdi/ej0iNljuGYvWp8VUxlfmyzk0iexSZnYnWkeRVcCoZ5ZilZYWzGsRg9piakT/a1Y6Zm1VkuG/4mHRJtInnf/i3ktqPUaaq6Gzc4oMdnGyk+gMCN/w4cqMjLKbMCPkyw7SemOu5OvabHoWKxhvE2lDiJ/VILQYkFkn7AmS0O9oECcKI97Z7hUBpBXx+zNYtAI75jOniSBM0OeWFQ9xmhtd5q0jYnaNMFgZ/sS+LA6rPEohKPYVnhMmAcKTOHuGDYML16c45YOmSydWbWFY17vcOZg9YsEllhaqMq+DMOzBSkHAdRo/EDH5Zy9YFcCC3Rf+X1aDHQS23l54Fpg1fYRKmO3xvtwIb5WMR6nFO1CAKKCOeASNEQjsGIDYwZ4xGRL25F8zo5crgO/K0aq4FkyEOETn536Oo7JOlnz5qPFBjjeB5HrvlwxVvCYQF6YUys6PgyxrQGB4hdLB3wfq9Kd0SBw7QWruZoIATwS4YkDgQ04CpAoIbHXCawQDZs0sYMqqr2ztCWQAdgauRyJhpELyjUdNT4zHXUW7dpwER0hyAkziA+YNmgPIfh26Js2SKpMY4J44MRZyUTsWghwsW9MO0ES7KAwJ5wmpOx3L4QYGysb0BvIYkjMaS2MMQGPAMdT9aGZZtyLTPeBJE5QXgfKFA1cR50R/CbKaKphBg/NwEszcAVH1Bs9NwYiGknLmRZs4OOYvaMnMkBsvuViapEZwB12ELYUUAq4tSCrjK4RsSTjHTvJWhUF5aJcAwY7ITQoNmv5xBCXMQnWECKqpWJ6106dC5/bLCYNILF2iDgJo1U3ljix5qyTUVhp+OrEZCYmLxA9bmnSQX6xfpJBQmauNpgNSr2vTrxrGLitAQekFSijgWIQQMTYtOkYZzATZlrWNpCwMj2jvRmIMwZNpuJKTMMIZhoC82tgLTRzadhCKym4SjFCcpNAG+QJFq712ZyAfeYOGzEHVjuRI9LvEd9IcJ/4vbvQcxnAM6MPkKNE0MScSKciz/F3jSi0mwy0GwH2ZhP1r8thUUNUV6AKKuFi0vWurP4xLWN35MyHNA9KEYwoSG9Tx+OitoEVIr8NPMqyYk7lcQexz8DESjuZgFqzhlW3hC+0mdoczZyGM3gu1kKfHcqEQz2aIJqjN8r/RQFewSpmsUpCPaT/N1xYoq5DkqnhoKzIBJ4WDkJGCD1zudHlNomFwm5IR0qKAYWEWwrEX1kzUQjBJo0Fe4+AAIgpU5mAG4XoISEN1FcfkEgGkS2Zyj6FOELvKwIZtBp/Hm/dJrYiKCoHsAgyfk2QkeoJowQOURA7bDisxN6HY3f0HDMsgGG4hbAkEoaB2QO0xqKjekMvzQiyVK0jyzDgJlf2Lp4Sa2iUBiED1GqMSpVF9jJYNBQ9T6g1HFMOSFJZIbiNM8J4EuwKw+qv4HBoNZMMCOP1wv8QJdgRyDoiWcbNsgYAEqILuEGm0RND1snWHihS6b0MpfYDQHOWYyDAxqSN1hDSBtmQSDVw4qrrUAhQx5iTKB2KslWspTm8Uv9gYlACNw3tigkAHi6D6hGgRlSXI0KJX9kUOKhtQLG94/JW6GkdHocaFIo6QHlBjWEoAcnN6aCJGJTIzcEuStfrqQEktaB0m4NCZb0LJ4zvpr2ifllMEQJfBesQO6SVSO4NeC6+LOyhtpQmvlhDg0Wt3lpiINOJn1gTcQcNSGB/lsqOq6ZJ20Coy2gQltjXWVoAW4YYTCjxE6kAN9hcYDewHCRdMRvNcd0BEnndYnLGzEFhIiAChIvXOWFY8BtWIQPCU6EcCh4QcY9RvRGT8FEG9NYgNHn2GEMAuqgJAUOhwh2osHqoK7FAqRmBwymeaDf6KFrMIUDDWWIB4H5Ew6FepS6wZ4y4iPU0dyR5i+gWCmJYJdwd1pRgWCj0A8D1Hn1GVAHyHITQKxinaYZskqbaBiqJYoNecpsQ6I8Bw4nQUSWK0UZvWTwL6N+JmwKfsBK0EbytSb9MXw2o36qO6Q0m3LtrR4HTa4NLkB4kqrkHcRH0aMPc9tIAXSHTDCVkjiWM1gkjCNZSoTWBiwQbxSQ5cDFGn0TrEoKCJmKCd+GIvYRoCF1SPRFqAQXmClxcABvxBKapdDDg8EHXokBNIdwXSuwIzzSjGLNOsOgBYyec8G4C+X0aVkMtAZ5BoJa4OSU2GTFCvhjleEXwV6CzwTCgkkzBQQSI7VUT0oScghGAQQmcoLMETOAK8kHREsUQ64aCAIDKeKYqmm2k68NaiAp6UWHvUYnDoN41aHTEGAamG7IaDZmB+FSJkERhR8h0auMmR5mQVVCQAxdHa4OFp4N6OzqSFJWCJrJyDYGiGN5G1Zo1XEfRwcMDSaLuA5Q0EYYGTZyNALmvQ4qJjjmsvagOuA0dMhKCg5Cah2K5Gfby4nHO1qyO5bXdLzHOHrMCwRGgzwBY2qhDRiB2qbtYnAC9tJhqFNxsBTJRwjFJg0gaE7bEANjnbTxQzHDMqs8YQWWskmYhZMvZsY16mQMrmnWSKAgEBcLWDfSewkVOULY16AkTIoALEQx4jqz9dBwQPRUxQ8iPCoZ4H4PQVQdlY9rTQgH579RnaVj9uoLakEURbzB8lSI6WfYGvU/0YXrPijUD6mrlKWB2UqSxhAEetAZ2g/0+4Ih1r43yTygsB+UyADBcUPjVbtAMQE7ylo6NQJjTwzLwggaFDIwbnRuvjhE+8P5kpok3h1UBkgTVH431DaY2qlAtSr3o9Mpibt1hUuVgWQV0n/mSH2J/VP8ZQxVQibEBSRzJFBtJCpVOlQVXIp1JKxaAoEmMAxVLBcIC0jWJgwM2YRU1wMJaG3KDGtMF4HAgdCarBYaHOgq4DUEtQQJNIGn2aT4E6wMaU7WYuY7zAMkBwT2xgsCkJxwyOBFAVwNMFpZDMsdSk/UCOVwHYdUdRCPCLajdCQaP3k8BCy82KwphawcBhimYkggawAUfREiTOdELpMraiRgDMFSin4BlbY6sckQsKZqNsHDCK6kNayUAgT7Ll2AVPGFoXJNYnRO8Y5VS+NGxtZIICaNK8ARYlFIgFlNSCLmWimpifQqKzYRLgM0kjW8HurTS4JVq10mcGXLKcAgMsF7NmChcZknkSIyl7WUCXNU0SJMYULCiIW1Y4gEpOODnQFuApM7bxp78l1SBai0TB8j+wDICvamEBU0WyTwmJGjMIhJL0v+3XqFINyLoOHCCcxsETQdx/VEIbxD8TbIMShXqKXtBUjEPQb5M/EADyyHu6h2yUwXXp6RaJ3rA1LZAlGKeQhbrAyia1y1jwDIcCAf86gnT2H32pLEwQFWJJpnBZKOUmhE3rAPPFdpKRkSJIsFFZG+nsd6hNBjINiYJ+pCHAyuiGFzFNxm1MSRoQyVVABcaab2qfwG2Jd0sa5thglBKCzoq8Qe2wiAAPYgd5EzXQli8CQH8iaytJNdQoyomXAOYtxDurEszViDEAby73sFiW0YF4A0mmdedF8ylmYIGUCcwn9W0iH5bBW8i/4VzWqlE5ggQp2JYi4Ofh7PkXdV8MWUdOA6ETSmI9rCIREsl6gzthfhbfpUiNWDEQx1LWsfs2gJ+NwG41tm3WFnDIdUlawrsgvAQqAn7Qyx8clOAFYKLM2ANE62l/BzhR6jsuOuWlJ+RalQAxDsYCSIG4ho8T1UYl4Hx7OHE5vFOvAMSQ99hcLNOWy2saXgPPmpb8VponBPGdtJqYTnirQn0K8RF9REOqQp3k62OptAMlnZq7FnRWomseUflYmVUNVR8AxWS9AOIfBM5kyQQAJgJ+UICm4m6mwU/r5NN5eyOXAFug5CWaVDZsfhUg/chBikDmDmICoprQHCCw6oGxtDA9hfrbhx4KyfsdT6ShTLQaljLTgAlJQhAUtKzZ9ZtheUNGYZkFAh/kLGhGSHPwCJTPzTtRI+N9ZZ2FOTAxvlEvrBUIcJGrZBRbUVgw6JX14nNgqVFitOEVxza8LQzkoM92ClrNyqbGUuNFAAfyGVAFvxI6ylCFrzsPluKKhYA6CsAIrznCX4/ob602chJSfc9KWKgLQRBADTCoLkmxMNLAmdBglg0cpJ2NHlF0nTi6TG/GLaBMAM4BTpPO7TqrKPl6wMKKCQ9m0YYGAlqyWpK2UR5bXjUairqLVcm0UNt4AHh6w2mJSTQZXWIUUoQBik0BEXqzkWTM+Z30SozTSQNv2qdYJ7rblrC9qksrdp8AIDjP6Mdq/bGGsBTgxy3kkt4j4GW5ovUxUxTSLXAgh/WosJ5AW5gt6sNzNzaMzyskPwG+FtgEs3eYBWTFTd4FUVpk/QBpFprislpgH59GnG/3HyNqBtZhgrfAl1NBkN4uFMXEY6ajAHTzdDYTF/FPGxBzIqXqcxiAxoN+Y14K51MiAmqMjGQedFjScgHgThbzSXg8C77gNwM9zTK6JZWWSefgmATPhksDWwoIGBULkIUJ5khbLUwNFhJkhRO9muHadkxmRdLbTJbQUoYH7/G4iMmAedwAHLuu3FGkMqEL0l8LaUaLLFMoMLhALSsE3eIop5eQDoIZJP1BuAqpOTgnDkBscaK7XhXTVtvDeChFaxV4BF4BWxgjM5G+ATnj43T0dvgdI3KskNKA+IG0QI9X5A1xhyxq5VsBMUH8AccMheAYpfBk4XEziZguciQyE4AfkCsloyNeao9KK1QhaDvpUibcHEE3przKmIyOACbkAtmSMFgowEDI4FzY5HYR9IWaXsrWEGiZZXsK+F5bH5YhUYG7RzRM11PnHTNRSLMFBMsJeh9oGOt7TmBQydRJTJgdAXw0XgVrQT/KTvgcKlJOQv4v14GmZsJAI6Jdu8LZAQDqjoT2yZltU4WA6g73hlQFEBcS5dkEJybbNrMj1IvJHZlEkrW1oZExKUkE4yWeVyERAhXIUzT6m9eHaAJLC+H0o2Yq2knr5dFjmOE8gV/WeQEzXxfgAxAZzISiXvDd8FEQvULlh/OsNdhtrROMJLLSYHaz/CTuuJiFsorVGrJmNaYY2amND9MHmjIpFNEBtMlBi/FMiFvjJf5FczhCpUw4yiMJKl7sPMwihenAF8ZUwuYFA+fJUY6pEMoBmuO3jtBXTQce5nItY2kebDGJsxPKJmpSlYanQ0AKO6UBWbhNkzCKnA9SSNrmUtLEgUyG/UNx6lpAm+8T7kRyXNCnYZlkVsdFxBqB1ZXn4v4rPVwJP5JJm8S5TNN1IO8gHSbACCsZugYZtcOltfByiYNJgvCskwjwIudUr8ML0GNJJixcYAt0INGYgHBnyRUw7mryOeB4UPOR890PDhm3CdKta4BfEXiB8BPYrkgX1DAEvsgeYg4HmjbAspyiZArGvAhGD9dcDeBYyJeMHiVZGlIH4ef+Ph1BtFAkBqajAK3DQQXmji5dguWk4RoJ3pZSU6szDcxkglBmFAZmoxMFSGEz+ahz9gQgZLH80xd1w2eRFpgCR4nGXJOwW5On0lQXQtdqpTWTuJlRlcUchCPqLBZyMTJIsM4VT26Bwsk5YuSok3S4eEnPiL9IQIimDXmQyYxBjfhSPBCMs1FQfoFaoRQ83rZ5oDk3RqJDcI85Rg27pApx/mxYdpUKLi4x8IoAx8YJQ9+Wh9L0rX6JHSlzhlYjLMGZ/DRS6+KrCPgRaXR14rNMRupHHxUslQgQMCCM8kz0HNA5Ief8mNkQpxMklZ6nxQ++fT13dXN9V5YhoKHNH9ZXPXjqzcUb7y+eEURmnMz63x//U6px88ufrqkDFpWYv218jYUlX1UpYdCFf+TNXqohUR0g/8optCg3hKFkSSppUgqg3WQRTEoj3Q+2ZQ9fF81qUCmgWVDWntWkvhSfv/hB1Q8A+bnhww0euQV56+w4fgPFjOcIkox7Vdl7BVzHdKfUxPuq/uv9nuRMMQrE84FmlYFEHr699dU6OTT599fPL/586PiUR/c3FKK7eL51Y9vVh2c+4n85pZ7U8FlTdab81xQi+zXi/BSfIQv9wq8lH95t/zcfQXeh/q4+6L5+Ob6u0sK8wiD/snF7XdX1wzYqpl7dU21j3eWyL9iNX18cXf5lqf5tZVEx/enfFwj6XHdqL/6zCXZ308e/DwOj0vx3T83X+7lDPc6hv/iNcx4fnz16upuFRxZE0wJk09vr767un4o+WIz+A+u1xjxi7LT/5ax2AuuPPu7Z1lxJSvk/YtXwv/2AcmwNHFuyE/jf8mAfHJ5cf3s7bO/e0ZIiwysXDl/dZn8B1sh/14D8vOzv3v2mxwR8/Z/+oh8fvf82YeXP+2rJKJTzvL/6E1zHhGWyTOIqPKfeUT+cTvcvL47HL++/vHly+2v/ZPlibJm9l5+7tPXT6tr/1rRTwyf+8Jz5/Kxn13cXlK8/p8P/+X28gXm9V+e9iFv9f6PL16surhLoa/2YVXBe//qmp/+7e3tp/QcsUT7058ub1+8vPkzBbU4RF3RF4/K5maNrJ8Px990LH0y90BGiM8f9u+onzVqr9AfILTwuoX9u7c/4z0GvEaMmz5JNTm8eO/2llJ7/+KKTQTgYChZpsOnE/SAavMJwFmpSQHY4KvWx5OKTbCHyNhS0n8AG7gG6JeoXsPjO9OSWrEJDO3wojSLUmXYQSzr5PCin2JjmMDVDHheaa4/LtnUioEpdCgseG7rZq2D+/esgdM2MtXIUqBmCtE2ugN2Chyu+aILeqgz8EtxaUk4Wo/aKAhEAMeSOsl51KOYhDR9WtZ2ajC5qPVKBZxBOZysQ+XNCDkRYqByUKJeeMwDbgzP8aRm08oZJZHUIIBxDeBZ6BwkSQNCZZUvkvxCpmToguQUmDmGh03OS15ElYmMKDRerdMKbInhlD5iPfGkRG584LhBmsqHaqDb8OqoOLPPOFmp5J5liQrcnLUQCMuT/kRshNRjnuVJzSaeXzpGeV7dNwoIQYdTgjVudIGSH8MmoFUmCRFxIWWIwk2WHNQieNgN0mD2aOBzacvwjwPte4E5TpWNjKbENohYE59Nzmmlv1qAvaw1ClI0WWuRuEeW5iBB9PxULdOnDJY68CBjmTUQKFZAIoWu2gtPajYBrBvsXO0GNY4HFcp+TSpJEAPICRFib5DYKAa1jRLUSKlGSbnMcYSOSzB7QNEkCFqMCgmkiHsA2DtR005NB5jlUDaIK4Oak4tE7QkyAUGUSUKQJP1S80NUZ6cYF2QTIyBF1Bm0PN+Bg6+tXYWQJKgEaZG1VsLmsTLqCjBociW7C6ghaRPkLE+IbIw4eAYBj2lUWeLhKkEDEo0gP+YAUFVDtCGyGnFzUJ8sUkXZkppFap4UbQKsJB+kkfvPCzIZRaoLgdaAvzb2BqkTEcSG88U9RApgjnbK22hnMQB3ZkGFLKRGDiKJuHAzFP5S9oM6WpQ0WUiJki3fWoAU86pYfEVKwyQfmJd1gs0TFANcP5frgPWY1PDk5ZcsVtVJ/ZvUJIKYZkowh5IzgLkgYkTNKfvCAJErUUmgcyWFICEbZAArOUkPuYrJYCS/tFKmdkFDCxjKslYQuEnjJ1018+W0kNEITClECpGkUAaICMPHIiGF7cOiJONkrmoTT2o2Gcx8qvcZdWqIEVA5ambRFQBrVt1AA2TK5iCBmaQEVjk56aSybrB+zEjMYTaTYVMBCLMcA7QS5q2rQQTSxqMp5V+CanKNGAqQIMkXs1O8sEPLK00pEpEqAsozT8aipdAIHBINsv6J+CULPOeRej5wOzR5P2SxDYWO3aFBMZwUm5mQ+cAloddRPwWkeE7S6wmhD2hqRG4JXcFLQ2YaBfHYCKCcQVoCCd+MBDwAJdnIUMeMf50jUykh4aUspeRPp2ZC0ivnL8s1tYIQb0RVO8JyKaoOJZ7KagD2bAUmgNocoKzZWbKgQDlgKEPMyTyD7gPyPfS2PeRJbQ4ECnsBgUipP2pbrkDWQD6T/kmVgwbPjK9gyfFOL1g2lL6DxEUCLWcYiUwNxMVJXu/UKYzm0I5WMjtx3oDNTc1AuB1EOWCZGalFBEidAinJcXUAQ1hijSpyZBhnTbCBEUE8jE3JFRWGFKG0ztvfmrNIeRIKQcAGooBOlsHJOnhZSStJsxlAycVBaRPSzXiDaWLhbBmHQ59FEBDHgHkkQKGIU7QEzKqE0PvS1ZCo4CoQqc53v1Huz8aUTMjK6Z4lDSbILYzqGZJ8bF61Qu2WzGiILE/Fl8hl5ptce5JI4SaZYsogjah0wcu6CzqeChISGQuZJUFU+Dw7MxC6lpNYREmVnB9g5Yp1gxXGDHZo7qTAYre1KCQeEDEj4pXzQ90pEhJmM4iAaOlJ9iKxdZiCpDaTRkKmM5kjY7CXKAFFXSlURqfCERVaoISRXEeWB7x4aKDEGUkjbBTUhHMPa4Zce6aESno7/ZBERbiv8AMTik6xDnBX8z16DheV2gPEPCGsZBJu0gqojZDJTMIGgidDvSZSb0YnjgyZObd4JyUyXJG6u6iqxNwaQRFT36xQ9iIIZpPQkoYfVEuHcUwGBpbRk2JNxD+z8FLWvMlyqKDNrZNEDes6JQvKDwknZP4RV6FKLYWkWJO8sRt8WsmHtVWNp6+ySmYkY0OaggwIBX6SCcPEZW6JkzUws3hKp9rXhE2JAkxaPDY9hWXIPtBMDrOohPkocpZkDNIhKUeZQgtaEfYfuw/7y2ARxZyVagOV+BsHqEyHAQIRoiV7pVLfAj5CZmOsJP45SIuheBmOTxIKU+3y5k4nxwGb2YYTtiGOE6QmGEFyYu6VREH4Y4T98iX3RIL7oBNkz6C2DR6Sk2O2OCA1c/WhyJCahYRcYWOK8CUbktTwDrrHBgumlvfR9b7qalLDj4K6FDojZXmRCp6WasK5QRAtlgR1NEmSJGJNbhZGMWke0EIJRToKuUYfSZ6XzHzTlBdZ75GCAthf8B5IqAkqJeWzw4SlsszEcUnObKMkB0EqVg/hH9YiFWbRTHAWoHFaFrLgvZu9FZhb1JgQtFEGSYkM9x5EnkaWDxuZVUgRDc8CUVnHgXRiiqauNGbSsoL6d1SKm3s5nyy8NwKWmmk4pS2YKF4JXJBXRoiRySSDsXbKJhI4bryUGIUSKGCjxGGmBqJ+yFhGnSd9NCNQZGVlyUCWdVbyq7CEccXwqzBZ2JGkyzUv0HexFiAYWIaxYjhscCeSnAojg8GNeGfGkXxMkDzKaJJuCGlPsdTxzCbiGH5/JZufDGucDMsSkNCWMQ/J8S+4mtD6LSvKIBreKdZkBa4UxO5GMgklmitswyy6WrPsSPKfMkUGcnHW8sQupKDhqrCWJRMxc8k3J+sd+mxS9GBhY39SZpNdSm00yo6SIazNSAuZZNZAD8tHpZaOTUl2kEFKyxosFLM1CkRmQSmC7FRZhNFMXSVsUJYyFXsyRS9TP4O8CErhIsNJnFdqJFHADbMF81OSMAgnOMmB1HSOSbovHhoJUs0zgVqDZCgSzxbXwuB0TRYFhbngn7ZJmixyjrc6F2hUpLT2wO7pnWywBr8wyAInkkyaS5t16kpIp9xMJeG8ipH7pjWr/EHQ64MqZCokP5MPRLFXKlWl54goaW4V7y0rTCvpkdBnWR4QIaj6Q+Ztur/YMNYr5VJIYchyM/As8NRXvR04W8SCyXW2uXgwT6s1Kb9LAaOkPGp6hdQeVAF9hwEzC2mLULgoK5H8Acg1TtyWUmJofAr6QdmiutXM5QFvTogWUw6NWmZU2F0ls6VlNhQbyGZyUJ2tHlBNKERG7SgUA3Vys4Qw7PvkF1JTMq38CW2PQtyRhZmQ5ewxzEuHJ4xCIdpMsluDPVBhW2UGeRrUULJh/Q5YXezKQQHXQHwINPDJmqUb3lZZW8i4GMJcAj/NcmsM2JVJzkriGCWNjI4Cz1D0gbR05wwSFG0fs0gTKWnYZL9C36JQG2405UIAU/At0vvAOVhZH+QdwQvHFlwcolzRkJQo2o45CRgCQ5yYNyhClmuDbBf4USkrnpZsyuyJCi9ikB0LxQBrnVA64BmOAjY9rCLqOKL9cuVinEEThxYBCYzEMOqcUaEizW3WnyOHIatgm1GwF7+KBBNLRhRcPQrodpxm2GlUeMg84WT+QmAXamzUFAZkAZC+CR8KPm2+NpcMMOygZWuQbjwoEglNKgTCYltYCFnOZkwlUiCWq7pKgqqxftRnltmA5EelT8JxrB/qtKAVgM9IfydnjzouJHtwQrLV4YEJpXlUyKBvBV4/HjgEoKwe7JCdGA4qAMHGpSIWOUMOx4HNwelULIIWmJoSK2EyiCThYLyRXVgdBkcyHOCHALThn2LPQEiAbqDABPk+ZSgmOZk9/XfsApCAoNRbvjsY4IJCx+QoJ7kaBG9Scoia33BRzqDW2bVISJYapBSbpCIItCEmEhIa1dhAcvAKOiXsAEcobMa2pmgDb0wn+X+Q/komsadZRtY6CDSM1fT+YKRV6j2Qzkw1VNDdRJAov0OmzDQ4O3OSP0aSPtUxcKIsJDLzn/QWReq1zm6kdwT88wUFVAtTaPBYdy3fD00RaWoXZGETXkWeNaE8ssBdDVaS1YZyBIOTVDOkgiVlDMpiWiEBvwuckYGB4Qdtm+xPiANIIlxM2NBpc0WUgENngEa4FpUtz9KE9OUQ2yHza3eyXigFMgtirsH6zeQos5L1oTqsvWR1U7Otkm6MuU1tGgVISX2a8imHl3IRBJvXC48h+SdvpkP/Ti65UWOLUtQx6gQ0yupJlEKDMJsQq1OgMsl5HQuUhNfWyYZOyuR5uTz8fyK3UB4ivQPWFTiskvuKaeIdQJPccPK+GJP1igSqSLgwLaQG4liSJuKgQOApg7pvwBhKEqSmoAGtB95VqIyZ621Z3YlqU9TyhoJN1SmSzBhWp4Aa6fWw8HhhgZOg03VqZg1QCQ4MNFMBwXnTKc/74oxWXgSCbCJ1PEiToTTCFHxR6rNUoiAOvT6z6UH7qDzJAQpMUUwA6RxM5kh2LikKCFHS6dDJFIliq0SnBD0JZpX0opUFwwKecBvzXdsk3zQKVzF9vFq+FmIizghTUwDbB2YeRdl66mk0aJAyhBlI3h/l5FtT8D/iHGnoAc9kISI0M+JfMH1g07HDsA0rLFyE4gQFq1TW4+UnlZfAU7sQsieFAzukTcxRTDZQItzVR8tkLZhGZiZOI0vQ81UdC1caECQyTZhVCXWK3HsqquN2p7pAy1I4ZJu4RdXII3Eh7U8rDFogJBh2lrUVgO/JsgWszppTpHujljDzs9YO1ciwskgYCbQUBdjYpYbAZ+LZCJTdRVA3dBDgmeBxoLizKAYhCuo9p2TwTiEJqIBU+oZUHYUUQEQxNdV65isM8kESg5Usi04QDNM4lGJxY1WEpvguJcnSwABBwOmhYHsaLTAScxa5gpxQakQMMWh9q5YJXjgVQKZT6IXaC3VmnWowWLJ7WMMMV1J7ydwj84Yxo/R0JlmRYkaRU/K6WV4EeihyQVJQCm9Ku4OTW3LsxHGC+L2aaCWAChVv8nsKxafnKdWhmYb0PaxEQXwiWvxYKiWMgwqC5I2S7YTizDrcT606+j1ND41GFRSK/TNje0wqcx2A93H9eT88kDT1vqhXSx0lOppvBxlUIKQYbxJBydogSwTL0qm926nVmqocjBEsEF8YMDKzcHjhCQ5UI3+Kuou8ywGZESiRLHpGPvEkhyizJQaBFEQPqaKoHfYzJYSz+u4qUo1tQzzBp6zqRjD/GmVr5szMQGsYY/ivlEjMOg2DypckHSb9mVeRUFAYAqgnAxduqkD/BYYcZI+Q5YGmSGEGGgaG2QfCzKOAb2S9M7iPmRlCLpVh4hJ5il6MNDBqdDdwbQSkQ1dnY3CPCTrHGxQwGK0m4qe4RvkynNqyKiUqWiiHJ8kUzvJtjRQdims7RtEqZCvTiLUlQINjRxEKMsAwipi9JP3ywhRsZCjD1IGeSFFWArGdboBC1HFIXie4B9U0KZH4oMTOnwhokKI2wHrRGBCsAf2UAG6lChXVj6kmToSRyh7UY6f4TkpuSs8pJd80ExQIAydISqIi7yhAs2Z96GFUw1AKqFrWWcpEynyZBpAo6PvM1CsIpFmdHg3VcLgUNKpmJiXQJoZPJtLi/wt5Zk5gk+KXI3ENpzCsZ70NhgMYhBQ/8jvAZJHeJEDQ76TUDiomV1Ay4oFNpJLQz+ttSBwDWCNjUWAOY7egQkgvnxSMxCADqzcMaKQXUpTNrbwoh6lOWYRF1ng9i6cpT6l7JHcmya0fpZQIAUOqP3SlVAp5ScmizvR9xF8Sn9PFono9sU08f97hhUYg7YyyaIKfnuGZmckPVEckabYWkkgw5JF6I9UOoRnSvsiUYC1Up2w4STTOegLFpn6ooiXODIXz0klEJysaV7xnpEAGl1sWTVViHAtsox4lJjTOPeu6UTGEddzmqqwiVN/hrRWExRA+Bl6P7jeKwJFqT2oaCiKNAwMAHclkJloTWCAkfpA0kK+6yiRq6m1Te43Kg5mdSQ4zD6hOhjCLFg/DJ2V29xWGNKodCCfrnDoluohhhNgkXQ9HwYiq9N1Cq6QdUL6p9trADSsRVsBZo8jJ8GKUh4reNd9d4HDwlVKRk52daUi8lYOkMxKqJN202gyDigIWWQeW6vnEw0Cwl5UXQnZlFlBN2ygdGSRlvnpnkP5MqfpJlcb1MpWawWFkOtVqdFJpBgOOWjkrdoWFxjt/yEZIYCRLgFOrR/9Hc/fSY7dxhGH4rwiz7hyQ7OZNy1yMBEGQAF4EMJCFbI3sCSSNIMmXGPB/D56vSY3myF44SIBAc9MZDnnY3eyurvrqrR3+AokBSXJZTX67ISavSBTaimQADSsRuxxhboekCozyYlTAUSbpYeycP8H1DTWuILA8nkPOz9R8kUoVQbyxD2SpoEMnf8PRS/Jgh5gag6Gg5eBdZdfNDWpm3kXiRzsvttMs5SApZuoGKTmR/BRYtqXJ1jQ2UMX7DLSto3DL5JEJnYnbcJ6Fdo6ET3VwFHchGqi6LrkhNjUoQvsSAme1u0DtEkGQP+pdyoE9WDwt0fVVqZVk/HLQb2gIWaBwvgEJIPtwD4x1GABQOU5OLj1u2KryjlzsJKIGX4jUIht+21XygA0TDxymkMkwE4WT5HlJwQPewULg3W89SW8e5pYdtLCj54V6iYcPGbUuaBVQQ2PKgXU3Tw/RNUayOckqQnEB27DoJ9b5PmG4pggAwMBAEiOUlzplZFF8saIwSy8Tdg6V87u4Qww5Onn5nxcpIeKMeMIYrZfKoJbtCnU4cWhbmb1R+epTfOACKHLD1LsSqAWhsry2cVmSI9uqxP1NScBevkakPGEIrQzxTyzPGQA0peTBwC6BS4AwXhfCD0nTgJQdYaj0nXCQ9OjQdJBj5JNxW4S3jqjNX5nqXYyPcXInKwb7shtB+AwMK2yx0BLgNqvCBJw8qqwNEsZA2KyToa8ymZrkmvD9icI6NWlSuMoohUKjCdIJivvYsgQyZT3fbUBGAO11NfbysAgAtllkVhUeE1/jJEaSn5Nxaxnh9qkY1tZWShIsgBmcSgsLjjA8pgUVCddUsTEBc8imbnlAUW/WuEHRx7xJEZOQLVg8sq3sQfjT4jVWpwg1l98qOrdwvAcQvajBzlHUv1e4nU0iEM47SheVIa9bVVBuTEqkwg20Uovnb75IWsRXALvr+X/Zy43sBDl5UKliYotSFD08NyK+x1yUQb10aqRiE5KiwijdlUDjqNAa2VJX3hbRW1ndq/TzJAlBe8+zzdmOH4Q1YrPU4wuKVIxrE+rQa7yUQrjyZuGd1BmiQbEbhh0eWJ3SpkUo0cRmCJANop2ZIb1quAwLWrt3ozXIfcx2hhBgzy46KfnKoOXah4ZDeWvSv+JjHu0r5TWmPGGLvSO3kD2z2vA0ugl2r1oRNiR88IuE2nlTLStBvaTNIy+moGaV0t+SCoiEk/xOhXZmsCdVFPrWfVwrjQN2e+Yn/KsFU3O2iVFcxaQnJmqagGbaGYAkZHtyMutlAYqvtp/VGnc9fIBWUaR4BTcm9yjPtaH1S5DsJtGOFyOtxyQtZCnwVuXCTYt0SfVU4h9QVBMelTCLQ5XRyJMF62DoVBUI8faoQMDA5PoSoTBEA9AS/WdbMD0m2q9U+rT5DCt/o4VT9idemnGU66fqET+9wmD8pbC5otxeGEa110Q3OuqZ/czWZfcYgwqDDBSsStaYOVaBJPnNEwM6Y9AGLTmd1u5hvti9rGPfd+QcIJ12ArjANmDwI/yVfApdcgGyWns4ejFhYdURzk4hQDTsjo1vjGTWIyhHcUGvGppaYWYf7OPM9wm7ywmsciVbUwnAaOEBsqtW1YTzl9HVqVucvEAV9ANWtE3oMcuQMUj02ZptD4+BGiozdASbgjOE2xNLjmCgXmhUuE6i7+lqF5IgWkC5xgiLH4aYo9dJoNe8mu0bu46bLpAzcST6goDzxVuoohHIduayB1shtuiYjE1ptkg41GXoQsOWeiOW6dkzQdLA6eJxHKdhdf/UhQFEpv4iZpu/Y1+qR7SFGYXuJ/9+4pgG4z3CobB00kLJMrILcoCAXnIwE9CwFw8/ixFrC4iduKutcgymMCjwkbBjtsbQaegQwWE7B/7v3Az0wEdAomgTweq6I2i4RIkiF1kdywVMBsaC0yEGV8aOzYoY8RR+BsFo0qlVqBEQnWf19DLR0i9CbNrICzT1Gotr9py2/BPz05wIaSJDOVjrizl0XeEIdoYQyzkuC0IQe40Lu5NzyrMV1ayAroIlFsil1wBuj4FNiGPg7wKarF5FPzSHMlSKd8YwAvxS45KowqA2IM0t5touchtk53PTHLE6HCkEqpSdskyvJJb0u2tcQfViyoMSXVZO1KqOLruXXJ1SmvconmspxDKWqWAAQbYx4MBuXcZkJrIAkrJP4uXfFfpCaFkZJYsqIRa0zZpR4+4wbdIVqcYmjgLzu8u1XvX0rv4K3oIAO48TVxsrpS1H8RXwYmF2ZT9VJdo9BSFCAtCGkEgFsIFG970o3AVK8ObZULOrLYpwitYRUi6XTZZzyyLR4pOy/I1mDet3vGsiyfgrgBccKB7pkM0US7MTz6Jlwleaj2JD5JLtRGHnsY8MRmFL4L9MEZGiB6+7qkTVkQ8Py5Ni1qtdNd2Emly0NnaQhCFgDnEojlgslYQBfD8uAE4//1daITrJiigR7M/hhlID22JkbbJ41WZDS1Yb5IWyb40kZ24qECVOBk4ooB5xViDkKIEM8oWVtnF4wPGT9lM9LpdUF6MRsMczJ8Furbu9ksRvy9u60My6MDUmcV9E13xflFVq3EAxBtHOhlKSciLIJJPbUogxjDzFc+X+8xZG4sAUHqPqgaYURSWCA56MOmtDNZwhMFP0DjKMxnGuHqmKoL5QMtvjqX06Xexx1MdQASNiNarZo8qbONWkKMvojlHBTQd8SHwqBOt5kqOmZpDyaHTNPkxlU5GOWtTiUm3vXYr/oPdvxyhpGO14DA8ezbXh7KVIGr594CSLbbXE+Y9HEdk6SqZQ/dgDHFQWKd4Q1rz6LhzT/DepjKI4S+ygRc0DaTM1i0wV21PBIVHgRsNczfABCMPaCL1uI7H5EDAKhhkhKjpoilOAsnNwEiVnS6FkVAiKTIWIM1PUN4WnUeVb7FmqiHlVQBrQLwakbANTc6MSQMulxhBXj6UFVuqRJ4hIBEjBauk03TYdhP7dDlVqBAEWthUAd5u6XBtOeydqDThAbVvLuPrQoaBh6ygUaI3mRkwCv/B9StWL0tptI6nTnHTm3WXalDYR61N5y2bIkmV6xEUNOYKyFjqbgR+JxBW0SfSvVZ5phWaPMvHqAQr5ifuLpzIqsomK5C1FlYEPxGIlUgV+Zu1ZJtwpGOpOFgjhXkFDhRqmXYVJMenw0QO+2ZvFKl0Byzle+Oxgj0YAJCI771zrc/qGk0tCTIautkenne2pe18VJ0yIb1Ro0GTEN5hZYDSVjf3deQH7N64gzdbA4wxputY5roGQazHrrRxJ+RjIaNYK2tdRIVKdbH2VABPWYYT4M6b0Jj63V5h8dZ6z0ye1bVw+CseZAkfLmeyMDVZFKgY/o+2sTAC6Pb4SCMTQGiqJA+EpAXqKYZG+8RpwicGnb90kuII20TwsW9rZZqtvV9ilWJETisUi+p7KplSVJrU1FY/hTwg7AvtMEHK0A7UjWVRrtYzZNhhxzTybeW1RaYLDXJ7H6hQc9rST2knUv3J/QJyIh7FWuT0ZcDYQPZvEcilVxuKZMIh8K9IM9RTI/UcAfY4WxO/0n+J4MWS65xt3XfoWkbDuEe8J7T57RYAsc/EAOhffeNQHCj8l9WV0b80EnRKgqRSk9m1K662EESynFJihkeL8i8VmAbRvFjY0Z1rWCQOEWMW5s+Kqo4hEKbSdJ/IK3KTREHhhOYQJjqUdLZqSWm3syBLipGChtjlV8JjoFXHFfldh3e5ssiu3lY/blJmrUpjKCQl2G/ZyF2j/PObkBYvcAClMcWXsWmpcl+Bpukih7TbIkXvGpQdQnkl3tO2GdiUP410NtkXtUGIoFsyuCuwK94VyiO/XrwJHL5KTSoL2XnWfZCigevBX0oJUaW4p5JG9vDYir8sOwC5nFLhd4wLpt2/Tudk+wM0J1eP27APZ3pnSliUytj61354eMotSg0mnSAddsZuYTZzg3MYCATGbekXtWToDT5XQLBXORH3hI/syzHnGCr6f/uFYgEbgdneTxnSnVZMWwXsQMvCLaJJ4X9Xj2CSOzJcWzgx0PgmSIxjqNjbNHiAb1sH9S2HJCmxrvk3zrnZSNJgy8uaJ84MaCZQc8XKnjc6sz+IEXKfg7uYDtcTcVA5leh4WY9PCLALLJbAe7Tqm9qnLu2I2aXUemGHZ5N6YfxrXxraulcRWrgY1m5zCET5VriNRjd3ASP0jsiEysJJaxETS82nCdcKWHtRZ1QaEZZCvbQh+NrgoRCeo12FKngrvHlQ+YG6eY74GzlJjfoxskKBKNVzB7GyY10D/V2Q4TkZKT0lk5nNJIbBPtAFLqkdxIhJUKf4hOATqq37CijZku3to0dLrEUrQNmZ2uEY2JU1X8IjDap5IQQwUM48coYV/yVI8K7gmw2mEzLwEum3P16hyerch1gsFU1i5BUAnVqDKiF1XHT/GInCVJ2NU/oiKS32r3vcgNLzKvEYtol1DOMlGCqG4crLhxOQDbmNIGiLMfpz5KQe1tYkcgMcufWt7jWyqssBSAU0oIQmS4mjg4wuXdH87VaABWV0ampd4cezcJXClTr3EuWQbSMECRsqOT+qQzGj5RtRcnEAmwFU4xWiVVWjDgqLjBW2VCJxQjD/hYfCAqqWuNUUx+kRJUuIFMTQ1R7iWU2D0mtlEgUIG41GIBKUb3ZXBDTmPcCuZmoIoedjKT0jdOI8jOpllairRE2vgCtpUyXdmuE2lwdv8KbSJDc5xItlxF4/6BNrUeG2ARFN+fV4+hTbhs1riCf82gPqbN1fQJnULeeyJktd9m0Cbfnv3+uC6/Io8/HP/8rA96UaQJH8cqsC+DIPjiNE/c/qj9Ph+hN+dR54/9dfSW48S0f1+On/76C/PMwz9Oo9Svsf8jSPMpVOpRd+d13IPOa/39yiDespxvi75q1pqmctSWj6dxXmmfsVHScL9fU6l5XMuc9nKWrZ87/v9MpextNzC2G/0UXqrU7vcVHZHjmV3MN102ctS5lLLUtYc5FpDP8V285Qi4+a7oA5zX1OZSytrmYu6w74wrH3YKwBlO7tXVr8cypab7TemudzgQzefPx0XzRHiGd4kx36Jm6VMQ6Hi2PPRCp9/rpnza0L3r3EyKB5lobnaXBjqvlDYOaeTyAYq1Ua/cEzInylVxc5cVEvtZS013XMMt6usJ8NKyxoICqOUJeeech2VA0qESkVqgepffjAlU6iksfJHesR1PnT9VaqM87dS0+FuO63AfYkgX8SJxe0KlTFEvxfBWKnUCg/5kDuZj85NN8Ycepxg4R7SZex9rc02HkrqTjmlunX7WnbBAxdy0cGtKPGeO0ZGztic+ui5UvEbfrpClqqenfC+nWCRTVfkjuxlH30o5scG3emiC+WbakfpHFWfXMnYN177qDyHUf/uVjwrfYhkDE2Lzp6HQodhOdHyCZwp96UWqoiYsqEc2oUgfiuNx7xUBrY3bCT3Hn98Of/r49vDIcgDnOlrG7QR04N33mVUM5C0rCCYJ+X89JaAxbWEKLZMb51tBGZcf3LNfuNu1VOoq+V75ZLKtKvJrFJMAUbpvimd5Hl1M25bmHBwTZvU3hkkhWlaAYzzCufd9vdgxHv8ZX84uioCqIOkWjNf9Y/7w0FPBQ5+MTRZVWeILQoHXcaup001RB3qnOcUfF7x/K4ztb8jd6MmVZUNfWdjyDkF/Yg5yPRj55eC2qStGVXUWP5O52QQmfdMmH0MnZfq33WnOYi97j5lTw2Fbszg9OllFzgb1M86XwjgGMDkGGOhFdoziDI2drdxPIyPr+l/utO4NQoz6WXck1BOuU/tq2anokrhtRpHFJfkRo5YW2FTLUU6XPUkL+nPfkM/d0GX7LN4JnJXbJmtaFj6zebh6DNARlSmLdkDhXBu8wBrB+mqetIT8GEp+rmB21fMfiwFaiZh/HdXJE/2q2MEi3/j4/YbzCQqbTEPh77RGKUdY6eP2Y/v8jQGLG882j7M+vzOmSu34rkYhrKp2Gpo0rhmnrCpAGXNn5nbPIj1mALOh+G8WH9Y/N7qLG8yw0xtIQ26KFWod/ww++gTm/lTochMMTo3C7PluNeub48DQ320mmaMMxOEWWXJ4jUU0tEh8zbcgE6sq0GrWzLwdc5QVsP6w7C/7qCzxXqr9hlIZ/RbIrCxQnog5CtnQe5f+/Lv9N5+b4/+tbdRP/HH7aW73EtWQbfDaZfhM2QC5QjOYiSDPb/ijsyRQx6UsfRUwStfaL8l5+69ccwu2jc34nZki+YHb9jLzujhs54fdsS1D+7hvNaxvlxm/MWe8s4yxGL2aOQ+dVv1ajdMrl1EnsreCBrHDJcVxOrWjQkzhaXfp/vph/fl+8qbcTZ171l/4l589pX/nAFOS+ZqK9/P7iRarT/B/aduVGoUV0+DX29mzz/u78G9PPz00Vu+3gueh2kfs3Cf/vM1tv01Czfn0l7OnvPmsCugrV2Bx+Yjgq2XrpC1XnrMqM1LV1BaG5o/vH17//Yv988P9NkXr4IQLTcvvnj17Iebp1MwjO++ffX9jzdPheAqBz+0EETOgSv7EdcV1uUC1CAprAGSIqd9++p7oN9//PQJQfLvt19+/vrZm3ff3MPb/iJk9K9f/vP2q/d/+v2BPH1grb48WKt/vnv9/EC1nqe7hp7i4uJM/tJ1/otYXM2FifubSaw0OY01GbgHFVcAfVfTfqLo3U4q7s8dHVhu59tSx9rs5+QH3zZICOEP6UdA4y8+UHEFG5NPCgAxubsPVNzI/FDXZxissHo/UHEPlq48GsXdH1Nx3dh/SMU9e+rgxb6/e//y9tdADCfvRf99ROF9oMP+MhUXn/n/gor7t7f3L+5e3j65f/HkzY9Pvrt9++7bd0/e/PDk2evnT97869c0RZ7Gh5Y4//vZPcih//2PObki+YoEpuLQMF2zFANXfPH5V2/v3rx/dzyxf7z7+puXd19/8/53969f3371/ubpi2cv392Wmxef3f1w+7wjePPST/8GmBDXxg==').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779221362338', obj, '');\n",
       "});\n",
       "\n",
       "      }\n",
       "      const servers = ['/static/', 'https://root.cern/js/7.11.0/', 'https://jsroot.gsi.de/7.11.0/'],\n",
       "            path = 'build/jsroot';\n",
       "      if (typeof JSROOT !== 'undefined')\n",
       "         execCode(JSROOT);\n",
       "      else if (typeof requirejs !== 'undefined') {\n",
       "         servers.forEach((s,i) => { servers[i] = s + path; });\n",
       "         requirejs.config({ paths: { 'jsroot' : servers } })(['jsroot'],  execCode);\n",
       "      } else {\n",
       "         const config = document.getElementById('jupyter-config-data');\n",
       "         if (config)\n",
       "            servers[0] = (JSON.parse(config.innerHTML || '{}')?.baseUrl || '/') + 'static/';\n",
       "         else\n",
       "            servers.shift();\n",
       "         function loadJsroot() {\n",
       "            return !servers.length ? 0 : import(servers.shift() + path + '.js').catch(loadJsroot).then(() => execCode(JSROOT));\n",
       "         }\n",
       "         loadJsroot();\n",
       "      }\n",
       "   }\n",
       "   process_root_plot_1779221362338();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from ROOT import gROOT \n",
    "gROOT.GetListOfCanvases().Draw()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.12.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
