{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "dcb015fd",
   "metadata": {},
   "source": [
    "# gr007_multigraph\n",
    "Allowing to overlay different graphs can be useful for comparing different datasets\n",
    "or for plotting multiple related graphs on the same canvas.\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "**Author:** Rene Brun, Jamie Gooding  \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:37 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "bd5b87ca",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:51.490303Z",
     "iopub.status.busy": "2026-05-19T20:37:51.490156Z",
     "iopub.status.idle": "2026-05-19T20:37:52.633872Z",
     "shell.execute_reply": "2026-05-19T20:37:52.633428Z"
    }
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import ROOT\n",
    "\n",
    "ROOT.gStyle.SetOptFit()\n",
    "c1 = ROOT.TCanvas(\"c1\", \"multigraph\", 700, 500)\n",
    "c1.SetGrid()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3ce00136",
   "metadata": {},
   "source": [
    "Initialize a TMultiGraph to hold multiple graphs\n",
    "This ensures the entire dataset from all added graphs is visible without manual range adjustments."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8a4fc13a",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:52.636218Z",
     "iopub.status.busy": "2026-05-19T20:37:52.636036Z",
     "iopub.status.idle": "2026-05-19T20:37:52.756128Z",
     "shell.execute_reply": "2026-05-19T20:37:52.755471Z"
    }
   },
   "outputs": [],
   "source": [
    "mg = ROOT.TMultiGraph()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "286b057b",
   "metadata": {},
   "source": [
    "Create first graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "0a996cc2",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:52.758183Z",
     "iopub.status.busy": "2026-05-19T20:37:52.758062Z",
     "iopub.status.idle": "2026-05-19T20:37:53.016399Z",
     "shell.execute_reply": "2026-05-19T20:37:53.016025Z"
    }
   },
   "outputs": [],
   "source": [
    "n1 = 10\n",
    "px1 = np.array([-0.1, 0.05, 0.25, 0.35, 0.5, 0.61, 0.7, 0.85, 0.89, 0.95])\n",
    "py1 = np.array([-1, 2.9, 5.6, 7.4, 9, 9.6, 8.7, 6.3, 4.5, 1])\n",
    "ex1 = np.array([0.05, 0.1, 0.07, 0.07, 0.04, 0.05, 0.06, 0.07, 0.08, 0.05])\n",
    "ey1 = np.array([0.8, 0.7, 0.6, 0.5, 0.4, 0.4, 0.5, 0.6, 0.7, 0.8])\n",
    "gr1 = ROOT.TGraphErrors(n1, px1, py1, ex1, ey1)\n",
    "gr1.SetMarkerColor(ROOT.kBlue)\n",
    "gr1.SetMarkerStyle(21)\n",
    "\n",
    "gr1.Fit(\"gaus\", \"q\")\n",
    "func1 = gr1.GetListOfFunctions().FindObject(\"gaus\")\n",
    "func1.SetLineColor(ROOT.kBlue)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e577a751",
   "metadata": {},
   "source": [
    "Add the first graph to the multigraph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a4021d91",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:53.022696Z",
     "iopub.status.busy": "2026-05-19T20:37:53.022549Z",
     "iopub.status.idle": "2026-05-19T20:37:53.130025Z",
     "shell.execute_reply": "2026-05-19T20:37:53.129433Z"
    }
   },
   "outputs": [],
   "source": [
    "mg.Add(gr1)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ae9b5754",
   "metadata": {},
   "source": [
    "Create second graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "ae2667d0",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:53.151376Z",
     "iopub.status.busy": "2026-05-19T20:37:53.151230Z",
     "iopub.status.idle": "2026-05-19T20:37:53.262449Z",
     "shell.execute_reply": "2026-05-19T20:37:53.262089Z"
    }
   },
   "outputs": [],
   "source": [
    "n2 = 10\n",
    "x2 = np.array([-0.28, 0.005, 0.19, 0.29, 0.45, 0.56, 0.65, 0.80, 0.90, 1.01])\n",
    "y2 = np.array([2.1, 3.86, 7, 9, 10, 10.55, 9.64, 7.26, 5.42, 2])\n",
    "ex2 = np.array([0.04, 0.12, 0.08, 0.06, 0.05, 0.04, 0.07, 0.06, 0.08, 0.04])\n",
    "ey2 = np.array([0.6, 0.8, 0.7, 0.4, 0.3, 0.3, 0.4, 0.5, 0.6, 0.7])\n",
    "gr2 = ROOT.TGraphErrors(n2, x2, y2, ex2, ey2)\n",
    "gr2.SetMarkerColor(ROOT.kRed)\n",
    "gr2.SetMarkerStyle(20)\n",
    "\n",
    "gr2.Fit(\"pol5\", \"q\")\n",
    "func2 = gr2.GetListOfFunctions().FindObject(\"pol5\")\n",
    "func2.SetLineColor(ROOT.kRed)\n",
    "func2.SetLineStyle(2)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "277ef65e",
   "metadata": {},
   "source": [
    "Add the second graph to the multigraph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7fed05f8",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:53.277358Z",
     "iopub.status.busy": "2026-05-19T20:37:53.277213Z",
     "iopub.status.idle": "2026-05-19T20:37:53.384572Z",
     "shell.execute_reply": "2026-05-19T20:37:53.384223Z"
    }
   },
   "outputs": [],
   "source": [
    "mg.Add(gr2)\n",
    "\n",
    "mg.Draw(\"ap\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "28fb4173",
   "metadata": {},
   "source": [
    "Force drawing of canvas to generate the fit TPaveStats"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "a9fd84f4",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:53.407353Z",
     "iopub.status.busy": "2026-05-19T20:37:53.407209Z",
     "iopub.status.idle": "2026-05-19T20:37:53.689893Z",
     "shell.execute_reply": "2026-05-19T20:37:53.668686Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779223073649\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779223073649() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(34939,'WkwIkysAe4gAeAHtnW2THDeOoP+KonY+zAuVS/CdmbcXIcvWeu4kS2HLY2m1vouSulqqdaurp7pkyzPh/37xgMyq6pbk1dpj7/jGDncrwWSSIAgCIACy/7r4v7tvL1bny5erxbh4eHt5/vXy8ovV08/OlxeXLza7hVmcfn6+/vOr1R8/XIzWLE4/WO8u29P9p/+xerajfEG1+xe79ea8A/97fX6yGL1ZnO5bGv/6tr7e1YH3wfkQzOL07vp8dXtzttkuRungZ7tvz1YH8Iv1ye5FA++sz856ZZAFnCtbq5+vTnf3ltvn6/PFaAdKPl0/f3Gt6IPNbrd5ebXaw83F1YJHp2uQcGZx+vjweKs90vCjy91yRy+1UucKdKtBfHxnu3y5uo43ZdcGvq93dUD7qnPx3CSfH1OGJj/YbE9W28/Wf+nUOyq8tzlZtXl9JIvxph1iTiVHK8H6aL0OkxdhkCwlR3E1Vh+ZoUduMcrg7XF9RkyxH5w7qk7t3ebW08sH69ers68WoxM/SBAfXHC2Jm1ut3mPt4vRpzqEHHMOUqp4JfKVtn0sQ67ivE2xDeDQ8jvfLcabLg/ikxQfU82psASutAtnfX5o6hhcjKkmszj905UPQoZ8fzp8crUAqmqNmS67zaOv3jEHSpn3eA/fWesyQ8++phSdUujQxeOv3jo9vQPevm2u59eKnvWpRJEUi/U+OZ3ds803n3x4u3HS42Pg0ecX+gIGfXz0/MW+9OP9062nl1cauvX08kpbt55eHj679fTy8OXnr1+ytG/awZcCQp9/2wrcELxi+Pnrl8vXMKzIXEFhGULRCg9frHbLxeiZ1wcv1v3p1uXF6tnu0+VuvWmD++TVy6erbXt+uH721evD47ft8e7meS+8u3l+KPtLe/tgefJguT5HQpjF6e3t5vLyxXLdG9yDDzZd2h4vXbipwYdVe29zsj5dr04W4+ny7HJlFqf/ul2fvF6Mu+2rGfp2D916enl7s9ke1f7oZL1bPkWi9Q/urF+vTq6Mem74wXb9cr1bf726fEOo311fojNmfdLB5Xa7GJ98aRabix0P35nF6UevV88uF+P5q7Mzszj9pCmgZ4I+eLjegcfi5auz3fr5dnnxgtJPXr18sDxb7XazNoBsn6xe794s/fCPnz24e+vxYlz8Zn40i9MPN6+enq0+eHV6Ok/bp6vdcn0OzfqYH12u/7L6/HJ+//gqqG8/XS3PFqOjc33dYAmDKzHD21+sz0823zzcXLAGj+HHx3CXzIcKH6/QQ50ZvpnlyO0Xiy4sbi93uzfIfWu3ayqb8T36YLX7ZrU67zrpCqQkvbPdvHy4uYD34aBHJ8sdIl+BxzOA5L7VAPnOLL66t/l6df9i+edXe8746tMVdLlaePrx+vmLuwyhq1/l0eXu2YuZuF999mLzzUdfr853n+2Wu1eXe0b96tar3QZ22Ne8tzp/9cFy22CY5dYz2G3/xemnq+XJ/fOzAz9/sd692LzaHXPmzK0fLy87r80lx7WeXLNK/ma2D5rsnbbPF6unuvbX58+ZuXdZWLfPlpeXfWlQr1lcxwUXsOrCji5G038mGa2xxk5OS3nyo9W3dgr6zsU4xf03dkq9JuW51+W5jKUaJ2JKmCrPJRgnaRI6DKb/TCKjVGekOCPZTeLGnIz+P4kfxTrTfyYJo7hs+s8kcZToTf+ZJI2Si+k/k+TR2WD6zyRldK6a/jNJHV1wpv9MzmplqXxvJydXQTdK9UaqGEllcn6UkowUQD+50PpN0Ui0k4ujFM8bIzFODqw6ksFOLo8So5HgjHg7uTKKjzrAkiZXR8nROJvAZPJ2FO+Mdp385GUUn8HASLKTd6N4ayQWQ8Pet35Kad9CKtf6iX7ycRTJRr8PcfJpFBvAQFH2eRRbjX7g3eTLKC4oXZ2tk68Hqsc4BTtKtkaHW2UKMjLFirOPU3Cjc2KUPoCQqjQkvZ1CGJXEzJoLU4ij0hQ0pE4hjTRTg6lpCrl1A6FCnkIZtUdRNKdQRwnZ6CDETdEqEspiU5RRQjLKs1N0SocO+LGW+TmMOczPcYwyP6fR5fk5j35+LDB2b6aO1gQYOE7Jjtak0p5ZOlXasxsZlC6rKbGCMDH1CxYR3KFABCi5AYk+pHWSdCX53ov2rotUpkT3Crg8ZfpXwIYpgwCAFDtlMFAgpinPa5j5zWCgb2yZMhgAlDhlRSBGk2TK80r2ecpz9+KnXMeGobFTsWNDkGcZow6DZzcWJQPPLOEmXexU4EolA0BbvvRnp8IySTPA0m1Y2amwcJUMACxbpbaxU+2CyxcAaRJJAkATXlIRX9XrG0kJIDQgOIAmxWAqO9Umxqq2nLVWFsqLPofMc9Vn542dxNpRREyIKhTEinJocCa4SaxrQqEyvTKJ9WPJRpRreB10HBKT8YBxdK4Yn00Mk9g0SrQms0DrJBaRkbSuJJpieTrjvXEOEJKYBFVlEpWtkTXkJhHEmBgvxgE5pSrCKNVJkK088+PiJAhXZhPC067EhnHjZ0G6WjG1GrF2EskjAo7xSJkE6QrfpWoEiV5HKZWOZRJnRwnIrGAkId+RY41uBdCN+pnNptRJnB+d137sJIjWElELhmUuyFbxRkWZo3LSpe6kNJyRrkcqwFkZq7Akq0yOKUEktp/JWWSTN/1ncsyJMB79mZyNTd4i1ifHpBSvIsBOzuaRtZnFZFouI9sjRUvy5Gwd+QphHSfHlMAFfDw5kTGLaR9PTpzK5f715JgU+uPjyUngy/bx5CTyZft4cswHEl0yX09OcsN27rTsB6e9VsV37hZ1hzhqSE9OpwRsW79MSRurFD85nRIdrPNhci7sR6s9Myl9tCDtkOIHqro8j7d9XPYD5mM/r2Fkmw2Tn5cxzITeYiV3Ld7gtpj1vXWTZz3P7xU+sgkUbqua+lEmz8Lu1RUsewtCQWRbMzeiTMEeDCM7BdtsAj4HcvsvgQ52ClBAnaBRpmChRpNejC/oYlbZr+MJLGZgtHH0U4B1rNuPL4htMO+tmwLGEiiiQRV29NWHF8QrZIPiL4oHVRmNqJLjY4VmUW+nICrq27gEUc8nlKNo+MBOwaFn2pBcn5JmM04Bc6gPCZUWWLR9SDpk+KMPqb1HkjgdcnuP9dGHBAlcmYfUqtd5TFrbW0A60MpeFGyzNgU/Kz391KvWa/M9Ba+Kd/+hat79d5BDJ22KsGTnEVqJVfZT3eDD1INRrIfJb+9DG08b4xRhydaPohxrAuU2a1OseYa0LR17G+AUqw4dyLopWaZAteKUbLPS+SZZxqyTNiXLiHXSpmSbtY7FABe2FaYoJsRXW2EdLjrkPm1TQqfwvk3blFSEKQe2LuHCxoE6pIRq0RXWmhOvzTFEEESzKBP2yo0e87DkePuQutWO1OPLtjr64NAqcGMjh1NGYNamBEvO4+M7JFYfH1ySEFlthXW4K/4+/UllVhufklqtdF1hHW7qrTFIcnUeHn1jpbfRKdTWZ6/pD5uo5BtRsEgS1jlrTJ8bOdqzsgaPeRaEDe0pd0HYyT5luK6jGGXKMF0H+Tjux6svD9PP8HI3bxSFKXcDh560chtde1lmGajSoSADO12jTMU2SaCTMpUuBFVIFLufc950dJrsKZ0bO75TmblRJU2BFzvvaR97YmNjNoncuuhMSGUqokq7WKGmSsLe5EEQ8qaLozbcgh6dBbJ1U0GPNrbXymUPaR97WYSRqwxI9wBdOzTqTAUVOst1mkVIwm1NGBdkZAf5uFvCHaXOjZ0KU+ncOPfU1ioQKLFn7B1hbHfZ1ZCCGfvYWBaFLWMnEcui+EakxjdTOQhIMPaqL8Cp1e1Suyvq4uf9Qat8kNnakW5R5k/D8R6lBJXYtKPtBpXYbfFNJehOqfNJmAU24wwHdMAuHNBpzezRYWBhj06rPKNDQ3HGRoFOoaaip8J2sc+LDiTu2UjbjWDUpXNpu8ZO0KnEhlEn+FQiGOlbXCFRMTqAYNTlY9E9ZNcHJTWF1qZtKm0f2SZtKrqR7IKztI1kw3Yqqan3xgpTSUdKlv51O3lotdFnHrXuKPefti3l/C53V8iMUO7ab0YpN8aeiZIbYzeJXXJj7HmguRlDeypkdjbzDrvkxtgzIdhotpdwVGmM3dRfKW3a5qqlzVpTGqXonDUtUYraQB0ddpzKW21y2XN2mak0Km3pz320pb/vQ5d+76J23aOfdUk9U6AL6t5/Nw5UTZUupDvQhOI8m11G7+ezy+i5mb0CUobuMlobql1Ed+AgoaFbPRLRcHBFSHddbN1Uu5huvdQup/tI6iyn+1TXLqcbMWu3GWaoWQx9pmsX1P1l964pT4EU2w0I336mirCEj9vPVJ3OYNOh1YUxerbV4vlU2TywNc5TdWmMVnfggU1/HkNloy6BRssYimHjGdJUXR1DMtEaCWWq3o4B8WkkylS9jCEYRSZM1bsxeKPiKU/V+zE4w/Y74UkIY7AGj1tyU/Vx9NXEii9sqj6NvpjUXITV59FnwxY500MZfTRsgDM9VNwmKRnJaarBjt6blI3kOtUgo3fs5tkb1+BGr1s+KXGqwY+uGvbCJU81hNEVdRJU3BoRt1WOeAuniictmZyNVHrIeAhyNVLpgc2qKXCyTDXU0TlTml1X1Y9mdOuXp4pwrKZEfFJTRTQWU5Jx4qca8R6YUoyTONUYRkmm4r0qU426IcQR5mSqyMVoqusQS17dE442se9wuTRIt801NShhV5maOwTLmMo+H1cNPgN8EU7qVJMHEhwXCioyyrUKxvbW97epgc1xUVNuIIMsU02FYYlllIA4ZA2uCQUzrk+Dq9RJnmoWBXGHKYjLx+ASaqBSSMR1EC+xwXuBC7wiGYsR9TUA4l1o/kIFcwPpl7ds5o3g5YDaSEZA0IhTLVZBPAUKSgNBg7f4rtXp0UDfQNAIU8Uhx1uGD4hD1ohj+ICpgaABmBsIGn6qCEgqgwZgVdCDhp9qbVjhexI3VTZRBe9jBxut8DeJTBUpiXcaNADbDHrQAGwz6EEDH5wylHjQAFSOUu+qgspSEuDvOlW29MEIXiyLLw2nHTAMX4BxRiEwZhivupGAO0Lf+1HwQCfjbKY+Hkw80jMcG1xmGLeMkVDxiFM/KxxBR2G0CmLHOBt5r+5UieADLHasXRUrKApqGIO3bsSZ3nwlYsUrCC76NigIKp63ccSlAyYKJsAEIngcceDhT5zBoiBo6Ns64s2DKrgYndVv1TsIKHSEOKvQ1DmQJPxRIZmauupulQrJsHVdhxkT1i4+SKQTaKsTr8MgquZuC3NUcMFlBObIL0Wmsb/KSLpnA1YJ2xgp9O+xWDpM/2zCCC1pVAIYD2PzKRbw8fhiOww+ns1h84vin7ReIwLqJ22e3IzgFdyzBWr4glgmcIMEF4t4p39g8EHA03+T92IR8fQPTP8IefpHA9B/8CiL5seEHiGgShoMPiGiaAR5nuk/JNSQFGhO/yGjpBqs/RdUWItWQY9QUXANhh7Rov40/JbAJwrKUd2iif6jQ5E2mP67TawxLuYjhjHyPQqd/rGLs5HCxoH+YxrRkQprf3lM4MN77a+M8I7CjD/WkRAXfl91lSeLB7HB9JdkxGfN+0B/SQOGDaa/5EdoQ4Qu0F8KI6FFhekvxZG+FKa/lNTRrTDjS3lUrz12hPZXRoJ/wOqOT3WsRP+A6S9bokoNhr6oBQvCyDs6RDGwwdECelSbGRKpRKQgjKiRVkCfqAcWTakaJhSLgsDLSoHTXhEsdEuBdktYgG4p0G6J27FGrQbbxKIoHKFFCphZVAWsqwV0i7LQiCM14C3UBcxI5BM/v0VhwI3UELpFZcB+hwJ88kQtqQH9URswmBaAGIojdDyIFlhUByymNUAM5YF1pgUQaPbBaQGYqhOu00PAVP0hnR4q6VWJQFMiFqCuaqTTVGU5ikQ5jxpgiipBLGkBmKJMmpwh1EgBBkGfWwo0CKSSoBhBzGoYiKWtBURd5kiQFmiEhrjD3AahFHSKBiKKQZ4SDVJuBia2oqEHhlYMAlUDQhCiZAKlogEhDaBng4AlIiQVfssGAUtMSCrDygYBS1hIiOMAa5QINw4IzXDzVCkMwhqGYEiYeMDsXMAnacSG0BBRS+ojgIkNERJWGPwIDikzF4KMQnSIIAvjQeASIHLK28UU8MOXgsQnnEv/ePc6qyNwNUSEhK/WIHCJEmnMGZj+8e4h8auYDL3wp+hCEIPAJUrkkPDA2j8RRnjaGQSu4N9DwgNDH5wquky8QeAKCobQVvWmBdLc6HTVBIPAFVUwMHwwCFxBwegiwugHjqPz8D9bAuA0EsNUmP5RMLqkkkHgCgom0F8iMCviW+4Cc4rAlWBHpwssGwSuEJJHggPTPzF5XW/FIICFfYQut2IQwMJOIujuwET6Zy+hq68aBLCwm1ADRbeYIuwn1M9gTYQe7CgUdgaBLOwpFPYm0D+7CoWjQSBLZI+KaZEMAlkiXlI1wAwCWSLeWyyRagL4RGK2GI6ioVWJ4AfsDQJaCNkrHAwCWiL48T4ZBLRE8APOBLNFIvgBV+PBJ4EfuxVrkNeSwA/YGcS1JPADDsaDTwI/4GgQ1pLADzgbDz1Sox88RJBSEvhZ4zCNwSeBH7AYJLck8AP2xik+4AdM2JxQK/gBJ4PYFlwyCrMbAgY/3mOJAzO/1cAzCHHJbX7JdRHwybG/Zw0Cpw5jkQCzX+R7rC5g+A+Y4C6JRB+8OiWX7DeamLe4c7ZZ7rxbmMWZZo/FaBZfL8Yn1QVTHVu/ZKpDthRTXTXVW1MJ8Xpnqvem+mCqj6b6ZKpHRhVTfTU1WFODmBqcqcGbGoKpgc0jOSLItmJqqKZGa2oUU6MzNXpTYzCVGHhMpkZkYjE1VlOTNTWJqcmZmrwhXbimaGpKpiZkaTE1VVOzNTWLqdmZmr2pOZia2YomU3M2NSOHq6nFmlrE1OJMLd7UEkwtbFKTqchSld/V1GpNJQZdnanVm0qCS42mVraw2dSK3Ne4On4YJLFVT4glAqGeUUssBhFsMRsswteSsGMRuxZZaxGwFrvaIlot8tTqrhhJahGfFivZIjgtdoJFRFosYYtwtEhEi3lgkYUWAWiRehZj1iLvSOPGfOYXX2jiEoLMYhFYRJhFblnMT4vEsrCcxdq0MJtFKlnYzLbkHb6AqyyOGYvMsZEv1CuHUWnR8Rb3i8WmtGh4izVpUe0Wb4tFp1s8LZYdiUWjW5wrFlVuMSUtOtyyQbAob4sLxWJDWlS3xVK36GyLDW5JILBYkBaT2qKqLQakRUdbTEeL7WkxGC1WoFXVrkYDatliLFkUssVutKhiW/kCs9GiiS0Go2XO8VpgNPCLJAfmvGU1Mec4ItD3/GLHz5yroYiXQdRCbPkQzLnah8Kct8QA5lzzpXANoFX5pYkU9MGcs+9HPfKLPphzNQhbRhZzruagMOds1tFd/OIL5pxNOPoIRwZ9MOdq/LGhRrHwiy+Yc03YUjuP7TFqgV98wZyrM1CYc82P04wXYc41wY1tLIKaX3zBnGuGne5CW2YYc67mm2bIqd2mm0I12FoKCXOu5pruztRO032X5PLld999Z36qbE4yyd+ZzdlOXnxPHmfP2CRZe3tOxTnDuZXc2Ddw7cgH4NGxjusnOg5nXYTDD28edrm33H612h4dnmkFR032gv2BkIer17tb589JvyY1GbC9tIOFBvr+bP38fIFXpMFH7fP6zoac9qQZxsvX6zdzxG/tdrcoJ337ZP31+nK9Ob9cjFHokTdHDd5dPl3NB3roT+HWQ6AHhe+fnl6u9KQNcrYX7tH2ivf62Vd3V+fPOSdkB0uSs87B/KmOhbTz65/tzuZ0832VuXsSkx//Ykao0/kDRvhvv5gR7ifovziHHyy3R+ekPlhuZ6aAZMqLLNOzB5+1NfHhdvlNO3zR4PsXu8NBjwb0sx4N6Mc97l/sPmz59u2sHEnwLCJdRvcvdl0kMIj7F7s7erasV72zpiYS6HoGPRW08GS944DcDD/cbM40gZ6CdlDl9uZ8t3m1vexnF27tOjrXJOat3Y5FrELqe2RBOyhDhfcTBqwVht/OJTEWIE5EDAp9dH7y0Xa76UfTWNkKanXW751X58+6WOAl4JEUA+xTyFsOvvTKjB+wV2bdAx7N993V89X5yfHpGrBrpUcSloYOhXPf85lCmpjFwr4i5OmsaBanH3MWYnV5TYb30s8uls84C6B97w/6HY1hf8qvl4Hjvt5VbPZV5+K56rWutd71cR8VHo4Wfby+hCGP8aGI9jo6yTLsud7ccaPOXHUu7RWvYUOte+vz9ctXL/9ttd0cjnrw4srZSxXx7dTLg+3qdLX917uH2q38iHCt4HiYYHpcehhnK/1wdfrxYoyWSduXfLEY89WSRwviAEdVHveCB8tj5nuwvMJbdL4vOvSsRW+eOX2wPLkydGj3YHny5hnWB8uTtxxjfbA8gdkfHcjTSx5fKUE39rNKdLh+9lU/qfRgedEOiz7qQmNf8HhBeHJx+tmz7Wp1fmf5TKUP6CHWjsgPyDo4YluKjudj/upo/VAF8LB6KJkZqLFVq7N9iRRdpCHoQRoKVWZUTm0CqYhhc9rAL+bzTrz7GEB0HLvt+uLD1bP1y+XZ5XwMqVkIXcC5vclwNDqV2deGp2XH44NJtPBogHt4r6qakYIB0s4277/ai0D9pg2n7msznrl9htOEz+aP5+er7acMj5osNW32cjE+4SDRjRv8EndD2PTfCDfKDKcb4YYkIGepEHqVA6RFx+XOXq2+f1dsf1fsjWBvhBtBa2rfP9GvLxEuq+XJaou21jNTSrY9dGe9uzMzTexMoyexmMX9G+Ulpdqz5Zl+zOz/r836nMLZELi9vDgGH65f7s3JXEqVElRi/PHl8vmKhvYC/vby/ORs9cWL9eVXq+2ny/Pn/TR2K/9g87qXtdlrpYrJ0XHNP603Z+vzubQfX2xVb6+3z86uS/v+ilOmIH2kAB9hcn/0+uLRsdkzFz4+Lnz8tppz4ZWaVLy3fP3h+rneFgAT3t/uXmxuL1+utssufX7C7ZrKlPkKg+smjuq4d23XkBsQ6EhgAc6yp789It81xfjmJqwfp+/HkfUM/f4scjs53w8i63n5/SnkYxFCtw3uKuNvTrrlBQRptzZ8L/HucSb3X+czuW87uAgFj44o6vrra1G/e3ND+H2Hhq+Zp9qCmobvvDlCvVx7FPaoznjsC/52s4zUvbLVTkdFnXMcVL2+2/7kYrM+5zoNtSQeXXfJJny16pLlPR7Zm1xZoVtxO7ho7OD5xU+iPBs7FKBSjR1q1IPWj9+jWTFuqCYOyeQhmGrqkEwZskmDN2GIRrQlzGzu+PhRM3hHj3i/jXHUjbmfuOfLVzrF87zN8NG0NTK/e3libB0tz+tOkjem7WiKDheMXJ+zR/P9Ai4h1h7tLxMA+uSCU8soi09O1i+bOf/JBbcNKEc8/PaiO1Y+uThVyagT/8mHdxZjxpZ8sb7886vldrUYubnCSbUhJFeCWh/dPF6MN/uu8d7yNfbyoeDBctvXxvjEDr7aGLOVUDL/wzfOp1KssxavbVBWclJKsq4EG0sIKNAHy+097lB4oseI5xLuSNASN0hxLtuaJEX8bphWS05nXw0oHHOvb+wL58ZkQ/HelpRSgmFtStZH7x23VAQL59pcpAabUrE+B6t4VuedVOeijTnZZOyAA7eEmkosvirHi8sxx+R9sdWTYz5IqF58ylli8SSjDpKDC6nULNG6Sh1nXY4SnbgUXaGO89F558Sn5ELQkuzESfI5SAiV3r34kmu0KQfrKwPxydaK+zoU9frbIYivKVRvY4hYrXYI2XtvXfAhRqIAdoi+Ju9trsnZ7Bl9Eh88QKpC6MQOqUZXvfPWl+AzJbnEXHy1PkbvKl+VEkOKNUhNJQRarjWUUqwUcoVw8Q74cnOWnGwm2jC4AA1STL5WHNWDLylbHOg2ZhzoYFdrskWqC3j8h2yjUxonrteQoRSXY0whF4/H3A02u5QTl8e4SHLSAHFx4dqQAqlNQ6jWhlRCCVKIug1Zso+OqIPXJmpMObuYigvZB+MHZ3N11toSXLLZ+CFk61IsIVefAzVy8OKjslQyYbAuZogusaScTRi8pFprscVmB5wkluq8lBSK54vq8NA7G+AGMXFwXpyDCSXhmI5DDPQqKboQXTRxKLGWkMVFidkiK4V1JtUXTzcmDaHkUpKrORCuTEOuTKmLNQUJJg+21grJS/TEMPLgayjZRWsz97WYPKRck1e+i57o9VCjEx9sdizZasrg0He1uFQdAaEyhOicSz4XV32kIOVYJdjik4JwbcWdTtwEUW8hRi7BOclE4QbnrC/VSq7VEWIbPKwVvYc1xJk6hORdkCrJIwXrEEN0ElJ2zlehzVhrSZ6h+0TkbUguZfFSo8sEMiiI0VXLEvHE24ZYk3cSRaqNBNjgPR9iyFwPRLLowJqpKVkuhiK8OPCvDY717In3dUbJxZcCe9bBhsy0Sqk212zKUBILOxSPQAhQJ+ZcKgNI2Xvo512ogeFVom5l4HahkBlYcrkyB67CaRZRFJnHFAGcz8kHJ8xjtj4Kk5m9R6narBS0qUgOFVZgQZQCR3gb4ZUkyTIFrvhY4KZQfPBibYw+oZwLMjL5ijwj3XmIUiSFnK0PPijP2lCJtFRbyTYLQ6lRQrKlukoAJgwRittqY7E5OhMGB+M7l2MphDjDYK1Eccm5EpJnZUllLVZJhdiqH0LIGXbKLpdCDSnJlUpc0RLZdUP1iaBU8dVnUh2HVHO24mFI7vRyA5PqJXuXCMK5AZEScpUcQyVJdrAxSarJp5xAQ4aSsreSi89ZI0ZDqraEEDOX9hCRG6JLyQYwj4Xw4+Az3CHRBV8Jrw3O+1hjlpRyIkw3IHJSTYhI8mHsUAut+FQcgTs7lJxSScXGkG2pYm5aVftI0mZcfcJm8mz9l+OtGWV/PN+tnm+5KQYrg43lq7PlmzuRXv6O0NH7GUS32Vg+WHK9126Ff/cJi6umGIPzKcBddogxoqtdUzNoEEdRCi7ahHJDi986Ozu08xk37RxuxtEXNP7XxW8u9JqiBf/8j4e4Ts6fm/X5zjzs49G697lp7H8ytPX2Eh/N7c355W55jmv5cvVsw/1wloDdD2zu3mpJKG1uSn5EU5+tn79cHrXlvoMa8+SMiycz6l/+fvX64rfcjfb73/729c0n4PDl7/79n59oC1/+7u2lv4MIR/ZgvzeqmYRYr8vtg+WWbQBXJJ3+aYU37Yil8J3P5O/XJd3evLzYXK7bhXuUvRH5fLD8eoWP653bJCKGR7Y0tu+7LWneHlnSytNvXK03b3SjLS5aQviYTwTJBAN4sJXsbYte0hvx2Pxik8bqQq2hhjzvkbk4zlHfFQk2YHhhaku7WGwQDPvHM1SAHrn+rt3zNkPVR7q0zmNs4Om7vq3+4/khuLI82XxztO+/vdniRVuerLmsiDHvrzh8uqW7w2aF2/zabuWdMZMeLJ3jJSB95DlsW5nmhCJ8eghu9q3z3c358xXrSK8lO9xE2AKe6/P/4l1c17bVd5e71et3iSFQf8uG/p+evVj/n7+67278+z/fOD85vfEvN3TPAphp6n1JMVPm2K+6DyPzsvumu1P6vXkWCt/lorLmnm6eVvjo/nb9fH2+97xyPWGJNtVSXMlvLKQfQppZXtz4lxs36pD8jX+6eHlDN2QJx+c/MmkQmTf4719uNK3UacPWtP6jE0fVyEwcF2PaE0fKL5o4+6v/VH9d/aWuJQ15t3h4k7Porm48vcVlz1v085Xgz4PldoXs/OviN9vV6WKU8t1bOtbuCCpunm+XL98wyT6WO++Sg1GFdpeEezfi7J/aFxwJqBYW/ZFq9Ypn8b1cVJ88W50RxiIhbXH6jiScnoHzNkfckcB/rdkfiCxiU10X/aKzdlBmT9f4MJsB9r3+vEetptpld5odC3HuLtHEPHEbMmaMUujl6sP15cXZ8ttDeJWg0J5RoaMmLB3uvby3Obm7fNrh70koer/Z+vbverYakf5rOVZHswUrz1e5ani3u15lEEIOP/9UvSMz6v2m6i9/11MFsXXF762knkX3Pelwb5sqprxP03/LFF1L5jrK9LLqkP/ofLddYzorb16+evnN0eO8rimeL/fV5/5i736fWfDgoHeNQ/EEzMkR9NBzsPpG77NXL+lDhct+Z4GM+GFxlr2OVQ3Xrrnt09cy/vXaW+3tg/U5GvSj7fY+F/LyHvj+16vt6dnmG3LEcAdst4iyq179wzEBse47JNZh1O8KS1wZ3EfvEebqES71+2uIQH9pBANnhrpmWhlOG9tcMR+9R6DLDnxAa3h3aItG+eGZMt6RSnw9k1jV+99/4FG3vz1k1WJYelPx9SDWDwk8OqW1TgDXLQ6OX0GdSxAu8aQ+NCODbQHD95gRN4jxQ0kmGxL2jbArwD2MX9Xhd+ek988ffrzYnMVju2eGj607JfaRdceiATxKRTxymvyNw48hEjHchx8LUAs/4uY4cje9X/gRP81R+DFzcX4u2TsXctLY5n++zo/Cj1eFxiEUmOZAYA4u4bhOsbgoxfjBhiguOF+zc3pEYLAcBLAh8IcGOCI/eOt9cYUoFxdX+EG8cylWjbJw4nyIEpMvhFscJwKOopnvxGiuRIDzeyu9Z5DzJmh4wpWZwxKpmpvihxCt+JhTtvj/b4ol6iXFSrWZs6M3y+A5LiFSaq0EQG6mwdlYBfLn4kM2NwM+bWdtLLHGFIO56XSegsOnX33J5qYMwUfJPtsixXp1WRcXbLQES2pVMZpKDjaHWBNnWnCFh+pLlVqjJA7quEFiDNXlEFItmcBeKjVacTGrgxy/P07zXFwqqbKIgydCQUTNclDGD6ngi3e5ViGK6ocSCdYWgpkcw/FDzUG8T+LFExsMygSEpkLKHJwKg+Vr8eKKJ94VBlsI5lZihZHziIPNqSRbeB8jsT0bQ/HZxpqDbwWueJdLtbVk0U+sd8GmAsUriBHI88Vam2yuURHLLlVLkCJFBltLSJYjVwR8tRdbS7TRJx9jK8DFHmqojpiDCYM4/nBGTsEHTrCFgehWxhtpXUPMM24hnOqCo0bIhJ5TgqotTkkoMXL+SUKhUcgZGJktrhAEIp6DZz9kl3MhcJljjc7nkEOJGrj0HNCOKZXACcA4sKYYjBTLEco0wDrOxuILB4rSQCS+hFIZahCiVZwVi7HYkFIlBGYJ5UVnYXEusRl8LMSWYw7Csco8pJKyEDyuidNbZYD/xZXoovdCXNLTpAaRss/km6RQeVeIvwUxRYOdpDC4mjmyVAcXLOFYBsLhyjpEbneAzMGyQuqQYyXeGVO0yRZTCelIIDBXIjFCsYMUX6vzVgicCzk0PjnrUopBD+cNoUriL3/EwIVHaKSaLMxNBgInwIaUQpUgRSCrFtTERBEu0kN1hKdicY6DSXoQb0ghhpJiyS7p0bkhJvK+PK5H7ST4nIlpRe9bmwTlapBAUgN9EBSuKSec4gy1BIkB5yUHvDWCml0NMTuEKnftDC57QnO5SOY4VhlqyDYU8dYTgDNliJUgXEnR5QzsXHWFdRNtDEw0YfpCPJiBZqKUjtCoZwn2uKazHPvKPhMDTwRora0kQSTCmMTcQkWT+MT5af4qiQ+OKGPmEFccsiXO56ol3gtLB+YvWSQPp81g6UAZDtrgqSGkEggCij8+QkEI0VaXiWppDecqwc/CSqiORn0OlehrikSCIykUTgL5GcmLN6mFQkkYgU9ZFjWXEIMjkF2EsG5J8HeBXTxXMA01p+xjzTiPlbe8j6F6h1jV43pDTfz1FYr1BL8fCLOW6CVnPaKYBjIRSiXuy2KUMuScfCTeaKPemDPk4lypJYYAwVwcHEMV510NORtXNVRSqw+WOIrxngB9ZXmTyJKNL4NPwYXivCXDJviheI6T+uqy9w4lgT3BVjqr1vxJ45qJeDlGSvujH7NpNbuYZviNwOY7tXMzKohde0LENgcYieQU1purOTr+jo/Pg5XoLFkvgZOPJQ3FuxJykRA473rT+zKgnGrJVrLlhqI8ZJciMp3Fp9T5aQOlF6RJz3HNHxMivSDVbm7oxwRIL0gTnxtyPyLSeuGPGvI/piFSI2eMwo9pCCt/biheC/z+9smF/fIPTy7ky9+//sOTC/fl7x/eW+5ejONnf/7t69/94cmF//L3F5tvfvvaeKAwQwEozlB8VwSY0HzbMfx/FAG2Q86VtAKbSOeQOQBchpJ98JlcsUpS1vx3wyRS0aIpuSSg/92w7wn/6p/K2od/s2Zh7gO+7MF+geFf3Uoewp7vDv86fJ1/n+Ff9owa/tUrJn4x4d/gMWGTlOyqqJS88icRf0j898JqiJMgpyqlHsbLweVfVojzp6CN7GmDgm4hTt39/0oatycNpopTthH8IL+Sxu9Jo0abkuYGLqFfSRP2pGnmq9JGnWO/0ibuaaOGfFtRaYjuF02afajnaibFtcQKBXW/8lNkV3j/a3YFNuzbr8d4vyDwP1B2xdvCFT9/yP4d17W832z9Y2RX2MGHo4NtMhTdpf78c/VreoU6xI5uDZjzljTup3/A/BAB5OlnT1b62dIrGg/uI5Azj/6aXzHfHrTenPfkvPfMryDnQTgH0lIoiOC33ArNstB8Cy3ToH9QK+L9Eiz4ak6yoC1OT/JzJcniyyPbgZSXhXbwt8h7+WGZnZok8ZYcd83YO9wG8Gs+ZzvS/Pd0C9tR2pmmdO3zOfthkn2GBDcy/3cIyV9tjv9Ujx1diNGO0+9vwfj5ddqvZsd/Ol3/MGYHZ0uvZDjuL275x7A8rt8vsTfANNVTdThXgK0vOLyoyvLj9fMXZ+vnL3a3N+fnq2e7Qyb8nfXr1UnLRz1dnl2uvvt/sRbVZQ==').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779223073649', 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_1779223073649();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "c1.Update()\n",
    "\n",
    "stats1 = gr1.GetListOfFunctions().FindObject(\"stats\")\n",
    "stats2 = gr2.GetListOfFunctions().FindObject(\"stats\")\n",
    "\n",
    "if stats1 and stats2:\n",
    "    stats1.SetTextColor(ROOT.kBlue)\n",
    "    stats2.SetTextColor(ROOT.kRed)\n",
    "    stats1.SetX1NDC(0.12)\n",
    "    stats1.SetX2NDC(0.32)\n",
    "    stats1.SetY1NDC(0.82)\n",
    "    stats2.SetX1NDC(0.72)\n",
    "    stats2.SetX2NDC(0.92)\n",
    "    stats2.SetY1NDC(0.75)\n",
    "    c1.Modified()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f56ef6f8",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "3ca1f10b",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:37:53.706516Z",
     "iopub.status.busy": "2026-05-19T20:37:53.706363Z",
     "iopub.status.idle": "2026-05-19T20:37:53.827076Z",
     "shell.execute_reply": "2026-05-19T20:37:53.826340Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779223073823\" style=\"width: 700px; height: 500px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779223073823() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(34939,'WkwIkysAe4gAeAHtnW2THDeOoP+KonY+zAuVS/CdmbcXIcvWeu4kS2HLY2m1vouSulqqdaurp7pkyzPh/37xgMyq6pbk1dpj7/jGDncrwWSSIAgCIACy/7r4v7tvL1bny5erxbh4eHt5/vXy8ovV08/OlxeXLza7hVmcfn6+/vOr1R8/XIzWLE4/WO8u29P9p/+xerajfEG1+xe79ea8A/97fX6yGL1ZnO5bGv/6tr7e1YH3wfkQzOL07vp8dXtzttkuRungZ7tvz1YH8Iv1ye5FA++sz856ZZAFnCtbq5+vTnf3ltvn6/PFaAdKPl0/f3Gt6IPNbrd5ebXaw83F1YJHp2uQcGZx+vjweKs90vCjy91yRy+1UucKdKtBfHxnu3y5uo43ZdcGvq93dUD7qnPx3CSfH1OGJj/YbE9W28/Wf+nUOyq8tzlZtXl9JIvxph1iTiVHK8H6aL0OkxdhkCwlR3E1Vh+ZoUduMcrg7XF9RkyxH5w7qk7t3ebW08sH69ers68WoxM/SBAfXHC2Jm1ut3mPt4vRpzqEHHMOUqp4JfKVtn0sQ67ivE2xDeDQ8jvfLcabLg/ikxQfU82psASutAtnfX5o6hhcjKkmszj905UPQoZ8fzp8crUAqmqNmS67zaOv3jEHSpn3eA/fWesyQ8++phSdUujQxeOv3jo9vQPevm2u59eKnvWpRJEUi/U+OZ3ds803n3x4u3HS42Pg0ecX+gIGfXz0/MW+9OP9062nl1cauvX08kpbt55eHj679fTy8OXnr1+ytG/awZcCQp9/2wrcELxi+Pnrl8vXMKzIXEFhGULRCg9frHbLxeiZ1wcv1v3p1uXF6tnu0+VuvWmD++TVy6erbXt+uH721evD47ft8e7meS+8u3l+KPtLe/tgefJguT5HQpjF6e3t5vLyxXLdG9yDDzZd2h4vXbipwYdVe29zsj5dr04W4+ny7HJlFqf/ul2fvF6Mu+2rGfp2D916enl7s9ke1f7oZL1bPkWi9Q/urF+vTq6Mem74wXb9cr1bf726fEOo311fojNmfdLB5Xa7GJ98aRabix0P35nF6UevV88uF+P5q7Mzszj9pCmgZ4I+eLjegcfi5auz3fr5dnnxgtJPXr18sDxb7XazNoBsn6xe794s/fCPnz24e+vxYlz8Zn40i9MPN6+enq0+eHV6Ok/bp6vdcn0OzfqYH12u/7L6/HJ+//gqqG8/XS3PFqOjc33dYAmDKzHD21+sz0823zzcXLAGj+HHx3CXzIcKH6/QQ50ZvpnlyO0Xiy4sbi93uzfIfWu3ayqb8T36YLX7ZrU67zrpCqQkvbPdvHy4uYD34aBHJ8sdIl+BxzOA5L7VAPnOLL66t/l6df9i+edXe8746tMVdLlaePrx+vmLuwyhq1/l0eXu2YuZuF999mLzzUdfr853n+2Wu1eXe0b96tar3QZ22Ne8tzp/9cFy22CY5dYz2G3/xemnq+XJ/fOzAz9/sd692LzaHXPmzK0fLy87r80lx7WeXLNK/ma2D5rsnbbPF6unuvbX58+ZuXdZWLfPlpeXfWlQr1lcxwUXsOrCji5G038mGa2xxk5OS3nyo9W3dgr6zsU4xf03dkq9JuW51+W5jKUaJ2JKmCrPJRgnaRI6DKb/TCKjVGekOCPZTeLGnIz+P4kfxTrTfyYJo7hs+s8kcZToTf+ZJI2Si+k/k+TR2WD6zyRldK6a/jNJHV1wpv9MzmplqXxvJydXQTdK9UaqGEllcn6UkowUQD+50PpN0Ui0k4ujFM8bIzFODqw6ksFOLo8So5HgjHg7uTKKjzrAkiZXR8nROJvAZPJ2FO+Mdp385GUUn8HASLKTd6N4ayQWQ8Pet35Kad9CKtf6iX7ycRTJRr8PcfJpFBvAQFH2eRRbjX7g3eTLKC4oXZ2tk68Hqsc4BTtKtkaHW2UKMjLFirOPU3Cjc2KUPoCQqjQkvZ1CGJXEzJoLU4ij0hQ0pE4hjTRTg6lpCrl1A6FCnkIZtUdRNKdQRwnZ6CDETdEqEspiU5RRQjLKs1N0SocO+LGW+TmMOczPcYwyP6fR5fk5j35+LDB2b6aO1gQYOE7Jjtak0p5ZOlXasxsZlC6rKbGCMDH1CxYR3KFABCi5AYk+pHWSdCX53ov2rotUpkT3Crg8ZfpXwIYpgwCAFDtlMFAgpinPa5j5zWCgb2yZMhgAlDhlRSBGk2TK80r2ecpz9+KnXMeGobFTsWNDkGcZow6DZzcWJQPPLOEmXexU4EolA0BbvvRnp8IySTPA0m1Y2amwcJUMACxbpbaxU+2CyxcAaRJJAkATXlIRX9XrG0kJIDQgOIAmxWAqO9Umxqq2nLVWFsqLPofMc9Vn542dxNpRREyIKhTEinJocCa4SaxrQqEyvTKJ9WPJRpRreB10HBKT8YBxdK4Yn00Mk9g0SrQms0DrJBaRkbSuJJpieTrjvXEOEJKYBFVlEpWtkTXkJhHEmBgvxgE5pSrCKNVJkK088+PiJAhXZhPC067EhnHjZ0G6WjG1GrF2EskjAo7xSJkE6QrfpWoEiV5HKZWOZRJnRwnIrGAkId+RY41uBdCN+pnNptRJnB+d137sJIjWElELhmUuyFbxRkWZo3LSpe6kNJyRrkcqwFkZq7Akq0yOKUEktp/JWWSTN/1ncsyJMB79mZyNTd4i1ifHpBSvIsBOzuaRtZnFZFouI9sjRUvy5Gwd+QphHSfHlMAFfDw5kTGLaR9PTpzK5f715JgU+uPjyUngy/bx5CTyZft4cswHEl0yX09OcsN27rTsB6e9VsV37hZ1hzhqSE9OpwRsW79MSRurFD85nRIdrPNhci7sR6s9Myl9tCDtkOIHqro8j7d9XPYD5mM/r2Fkmw2Tn5cxzITeYiV3Ld7gtpj1vXWTZz3P7xU+sgkUbqua+lEmz8Lu1RUsewtCQWRbMzeiTMEeDCM7BdtsAj4HcvsvgQ52ClBAnaBRpmChRpNejC/oYlbZr+MJLGZgtHH0U4B1rNuPL4htMO+tmwLGEiiiQRV29NWHF8QrZIPiL4oHVRmNqJLjY4VmUW+nICrq27gEUc8nlKNo+MBOwaFn2pBcn5JmM04Bc6gPCZUWWLR9SDpk+KMPqb1HkjgdcnuP9dGHBAlcmYfUqtd5TFrbW0A60MpeFGyzNgU/Kz391KvWa/M9Ba+Kd/+hat79d5BDJ22KsGTnEVqJVfZT3eDD1INRrIfJb+9DG08b4xRhydaPohxrAuU2a1OseYa0LR17G+AUqw4dyLopWaZAteKUbLPS+SZZxqyTNiXLiHXSpmSbtY7FABe2FaYoJsRXW2EdLjrkPm1TQqfwvk3blFSEKQe2LuHCxoE6pIRq0RXWmhOvzTFEEESzKBP2yo0e87DkePuQutWO1OPLtjr64NAqcGMjh1NGYNamBEvO4+M7JFYfH1ySEFlthXW4K/4+/UllVhufklqtdF1hHW7qrTFIcnUeHn1jpbfRKdTWZ6/pD5uo5BtRsEgS1jlrTJ8bOdqzsgaPeRaEDe0pd0HYyT5luK6jGGXKMF0H+Tjux6svD9PP8HI3bxSFKXcDh560chtde1lmGajSoSADO12jTMU2SaCTMpUuBFVIFLufc950dJrsKZ0bO75TmblRJU2BFzvvaR97YmNjNoncuuhMSGUqokq7WKGmSsLe5EEQ8qaLozbcgh6dBbJ1U0GPNrbXymUPaR97WYSRqwxI9wBdOzTqTAUVOst1mkVIwm1NGBdkZAf5uFvCHaXOjZ0KU+ncOPfU1ioQKLFn7B1hbHfZ1ZCCGfvYWBaFLWMnEcui+EakxjdTOQhIMPaqL8Cp1e1Suyvq4uf9Qat8kNnakW5R5k/D8R6lBJXYtKPtBpXYbfFNJehOqfNJmAU24wwHdMAuHNBpzezRYWBhj06rPKNDQ3HGRoFOoaaip8J2sc+LDiTu2UjbjWDUpXNpu8ZO0KnEhlEn+FQiGOlbXCFRMTqAYNTlY9E9ZNcHJTWF1qZtKm0f2SZtKrqR7IKztI1kw3Yqqan3xgpTSUdKlv51O3lotdFnHrXuKPefti3l/C53V8iMUO7ab0YpN8aeiZIbYzeJXXJj7HmguRlDeypkdjbzDrvkxtgzIdhotpdwVGmM3dRfKW3a5qqlzVpTGqXonDUtUYraQB0ddpzKW21y2XN2mak0Km3pz320pb/vQ5d+76J23aOfdUk9U6AL6t5/Nw5UTZUupDvQhOI8m11G7+ezy+i5mb0CUobuMlobql1Ed+AgoaFbPRLRcHBFSHddbN1Uu5huvdQup/tI6iyn+1TXLqcbMWu3GWaoWQx9pmsX1P1l964pT4EU2w0I336mirCEj9vPVJ3OYNOh1YUxerbV4vlU2TywNc5TdWmMVnfggU1/HkNloy6BRssYimHjGdJUXR1DMtEaCWWq3o4B8WkkylS9jCEYRSZM1bsxeKPiKU/V+zE4w/Y74UkIY7AGj1tyU/Vx9NXEii9sqj6NvpjUXITV59FnwxY500MZfTRsgDM9VNwmKRnJaarBjt6blI3kOtUgo3fs5tkb1+BGr1s+KXGqwY+uGvbCJU81hNEVdRJU3BoRt1WOeAuniictmZyNVHrIeAhyNVLpgc2qKXCyTDXU0TlTml1X1Y9mdOuXp4pwrKZEfFJTRTQWU5Jx4qca8R6YUoyTONUYRkmm4r0qU426IcQR5mSqyMVoqusQS17dE442se9wuTRIt801NShhV5maOwTLmMo+H1cNPgN8EU7qVJMHEhwXCioyyrUKxvbW97epgc1xUVNuIIMsU02FYYlllIA4ZA2uCQUzrk+Dq9RJnmoWBXGHKYjLx+ASaqBSSMR1EC+xwXuBC7wiGYsR9TUA4l1o/kIFcwPpl7ds5o3g5YDaSEZA0IhTLVZBPAUKSgNBg7f4rtXp0UDfQNAIU8Uhx1uGD4hD1ohj+ICpgaABmBsIGn6qCEgqgwZgVdCDhp9qbVjhexI3VTZRBe9jBxut8DeJTBUpiXcaNADbDHrQAGwz6EEDH5wylHjQAFSOUu+qgspSEuDvOlW29MEIXiyLLw2nHTAMX4BxRiEwZhivupGAO0Lf+1HwQCfjbKY+Hkw80jMcG1xmGLeMkVDxiFM/KxxBR2G0CmLHOBt5r+5UieADLHasXRUrKApqGIO3bsSZ3nwlYsUrCC76NigIKp63ccSlAyYKJsAEIngcceDhT5zBoiBo6Ns64s2DKrgYndVv1TsIKHSEOKvQ1DmQJPxRIZmauupulQrJsHVdhxkT1i4+SKQTaKsTr8MgquZuC3NUcMFlBObIL0Wmsb/KSLpnA1YJ2xgp9O+xWDpM/2zCCC1pVAIYD2PzKRbw8fhiOww+ns1h84vin7ReIwLqJ22e3IzgFdyzBWr4glgmcIMEF4t4p39g8EHA03+T92IR8fQPTP8IefpHA9B/8CiL5seEHiGgShoMPiGiaAR5nuk/JNSQFGhO/yGjpBqs/RdUWItWQY9QUXANhh7Rov40/JbAJwrKUd2iif6jQ5E2mP67TawxLuYjhjHyPQqd/rGLs5HCxoH+YxrRkQprf3lM4MN77a+M8I7CjD/WkRAXfl91lSeLB7HB9JdkxGfN+0B/SQOGDaa/5EdoQ4Qu0F8KI6FFhekvxZG+FKa/lNTRrTDjS3lUrz12hPZXRoJ/wOqOT3WsRP+A6S9bokoNhr6oBQvCyDs6RDGwwdECelSbGRKpRKQgjKiRVkCfqAcWTakaJhSLgsDLSoHTXhEsdEuBdktYgG4p0G6J27FGrQbbxKIoHKFFCphZVAWsqwV0i7LQiCM14C3UBcxI5BM/v0VhwI3UELpFZcB+hwJ88kQtqQH9URswmBaAGIojdDyIFlhUByymNUAM5YF1pgUQaPbBaQGYqhOu00PAVP0hnR4q6VWJQFMiFqCuaqTTVGU5ikQ5jxpgiipBLGkBmKJMmpwh1EgBBkGfWwo0CKSSoBhBzGoYiKWtBURd5kiQFmiEhrjD3AahFHSKBiKKQZ4SDVJuBia2oqEHhlYMAlUDQhCiZAKlogEhDaBng4AlIiQVfssGAUtMSCrDygYBS1hIiOMAa5QINw4IzXDzVCkMwhqGYEiYeMDsXMAnacSG0BBRS+ojgIkNERJWGPwIDikzF4KMQnSIIAvjQeASIHLK28UU8MOXgsQnnEv/ePc6qyNwNUSEhK/WIHCJEmnMGZj+8e4h8auYDL3wp+hCEIPAJUrkkPDA2j8RRnjaGQSu4N9DwgNDH5wquky8QeAKCobQVvWmBdLc6HTVBIPAFVUwMHwwCFxBwegiwugHjqPz8D9bAuA0EsNUmP5RMLqkkkHgCgom0F8iMCviW+4Cc4rAlWBHpwssGwSuEJJHggPTPzF5XW/FIICFfYQut2IQwMJOIujuwET6Zy+hq68aBLCwm1ADRbeYIuwn1M9gTYQe7CgUdgaBLOwpFPYm0D+7CoWjQSBLZI+KaZEMAlkiXlI1wAwCWSLeWyyRagL4RGK2GI6ioVWJ4AfsDQJaCNkrHAwCWiL48T4ZBLRE8APOBLNFIvgBV+PBJ4EfuxVrkNeSwA/YGcS1JPADDsaDTwI/4GgQ1pLADzgbDz1Sox88RJBSEvhZ4zCNwSeBH7AYJLck8AP2xik+4AdM2JxQK/gBJ4PYFlwyCrMbAgY/3mOJAzO/1cAzCHHJbX7JdRHwybG/Zw0Cpw5jkQCzX+R7rC5g+A+Y4C6JRB+8OiWX7DeamLe4c7ZZ7rxbmMWZZo/FaBZfL8Yn1QVTHVu/ZKpDthRTXTXVW1MJ8Xpnqvem+mCqj6b6ZKpHRhVTfTU1WFODmBqcqcGbGoKpgc0jOSLItmJqqKZGa2oUU6MzNXpTYzCVGHhMpkZkYjE1VlOTNTWJqcmZmrwhXbimaGpKpiZkaTE1VVOzNTWLqdmZmr2pOZia2YomU3M2NSOHq6nFmlrE1OJMLd7UEkwtbFKTqchSld/V1GpNJQZdnanVm0qCS42mVraw2dSK3Ne4On4YJLFVT4glAqGeUUssBhFsMRsswteSsGMRuxZZaxGwFrvaIlot8tTqrhhJahGfFivZIjgtdoJFRFosYYtwtEhEi3lgkYUWAWiRehZj1iLvSOPGfOYXX2jiEoLMYhFYRJhFblnMT4vEsrCcxdq0MJtFKlnYzLbkHb6AqyyOGYvMsZEv1CuHUWnR8Rb3i8WmtGh4izVpUe0Wb4tFp1s8LZYdiUWjW5wrFlVuMSUtOtyyQbAob4sLxWJDWlS3xVK36GyLDW5JILBYkBaT2qKqLQakRUdbTEeL7WkxGC1WoFXVrkYDatliLFkUssVutKhiW/kCs9GiiS0Go2XO8VpgNPCLJAfmvGU1Mec4ItD3/GLHz5yroYiXQdRCbPkQzLnah8Kct8QA5lzzpXANoFX5pYkU9MGcs+9HPfKLPphzNQhbRhZzruagMOds1tFd/OIL5pxNOPoIRwZ9MOdq/LGhRrHwiy+Yc03YUjuP7TFqgV98wZyrM1CYc82P04wXYc41wY1tLIKaX3zBnGuGne5CW2YYc67mm2bIqd2mm0I12FoKCXOu5pruztRO032X5PLld999Z36qbE4yyd+ZzdlOXnxPHmfP2CRZe3tOxTnDuZXc2Ddw7cgH4NGxjusnOg5nXYTDD28edrm33H612h4dnmkFR032gv2BkIer17tb589JvyY1GbC9tIOFBvr+bP38fIFXpMFH7fP6zoac9qQZxsvX6zdzxG/tdrcoJ337ZP31+nK9Ob9cjFHokTdHDd5dPl3NB3roT+HWQ6AHhe+fnl6u9KQNcrYX7tH2ivf62Vd3V+fPOSdkB0uSs87B/KmOhbTz65/tzuZ0832VuXsSkx//Ykao0/kDRvhvv5gR7ifovziHHyy3R+ekPlhuZ6aAZMqLLNOzB5+1NfHhdvlNO3zR4PsXu8NBjwb0sx4N6Mc97l/sPmz59u2sHEnwLCJdRvcvdl0kMIj7F7s7erasV72zpiYS6HoGPRW08GS944DcDD/cbM40gZ6CdlDl9uZ8t3m1vexnF27tOjrXJOat3Y5FrELqe2RBOyhDhfcTBqwVht/OJTEWIE5EDAp9dH7y0Xa76UfTWNkKanXW751X58+6WOAl4JEUA+xTyFsOvvTKjB+wV2bdAx7N993V89X5yfHpGrBrpUcSloYOhXPf85lCmpjFwr4i5OmsaBanH3MWYnV5TYb30s8uls84C6B97w/6HY1hf8qvl4Hjvt5VbPZV5+K56rWutd71cR8VHo4Wfby+hCGP8aGI9jo6yTLsud7ccaPOXHUu7RWvYUOte+vz9ctXL/9ttd0cjnrw4srZSxXx7dTLg+3qdLX917uH2q38iHCt4HiYYHpcehhnK/1wdfrxYoyWSduXfLEY89WSRwviAEdVHveCB8tj5nuwvMJbdL4vOvSsRW+eOX2wPLkydGj3YHny5hnWB8uTtxxjfbA8gdkfHcjTSx5fKUE39rNKdLh+9lU/qfRgedEOiz7qQmNf8HhBeHJx+tmz7Wp1fmf5TKUP6CHWjsgPyDo4YluKjudj/upo/VAF8LB6KJkZqLFVq7N9iRRdpCHoQRoKVWZUTm0CqYhhc9rAL+bzTrz7GEB0HLvt+uLD1bP1y+XZ5XwMqVkIXcC5vclwNDqV2deGp2XH44NJtPBogHt4r6qakYIB0s4277/ai0D9pg2n7msznrl9htOEz+aP5+er7acMj5osNW32cjE+4SDRjRv8EndD2PTfCDfKDKcb4YYkIGepEHqVA6RFx+XOXq2+f1dsf1fsjWBvhBtBa2rfP9GvLxEuq+XJaou21jNTSrY9dGe9uzMzTexMoyexmMX9G+Ulpdqz5Zl+zOz/r836nMLZELi9vDgGH65f7s3JXEqVElRi/PHl8vmKhvYC/vby/ORs9cWL9eVXq+2ny/Pn/TR2K/9g87qXtdlrpYrJ0XHNP603Z+vzubQfX2xVb6+3z86uS/v+ilOmIH2kAB9hcn/0+uLRsdkzFz4+Lnz8tppz4ZWaVLy3fP3h+rneFgAT3t/uXmxuL1+utssufX7C7ZrKlPkKg+smjuq4d23XkBsQ6EhgAc6yp789It81xfjmJqwfp+/HkfUM/f4scjs53w8i63n5/SnkYxFCtw3uKuNvTrrlBQRptzZ8L/HucSb3X+czuW87uAgFj44o6vrra1G/e3ND+H2Hhq+Zp9qCmobvvDlCvVx7FPaoznjsC/52s4zUvbLVTkdFnXMcVL2+2/7kYrM+5zoNtSQeXXfJJny16pLlPR7Zm1xZoVtxO7ho7OD5xU+iPBs7FKBSjR1q1IPWj9+jWTFuqCYOyeQhmGrqkEwZskmDN2GIRrQlzGzu+PhRM3hHj3i/jXHUjbmfuOfLVzrF87zN8NG0NTK/e3libB0tz+tOkjem7WiKDheMXJ+zR/P9Ai4h1h7tLxMA+uSCU8soi09O1i+bOf/JBbcNKEc8/PaiO1Y+uThVyagT/8mHdxZjxpZ8sb7886vldrUYubnCSbUhJFeCWh/dPF6MN/uu8d7yNfbyoeDBctvXxvjEDr7aGLOVUDL/wzfOp1KssxavbVBWclJKsq4EG0sIKNAHy+097lB4oseI5xLuSNASN0hxLtuaJEX8bphWS05nXw0oHHOvb+wL58ZkQ/HelpRSgmFtStZH7x23VAQL59pcpAabUrE+B6t4VuedVOeijTnZZOyAA7eEmkosvirHi8sxx+R9sdWTYz5IqF58ylli8SSjDpKDC6nULNG6Sh1nXY4SnbgUXaGO89F558Sn5ELQkuzESfI5SAiV3r34kmu0KQfrKwPxydaK+zoU9frbIYivKVRvY4hYrXYI2XtvXfAhRqIAdoi+Ju9trsnZ7Bl9Eh88QKpC6MQOqUZXvfPWl+AzJbnEXHy1PkbvKl+VEkOKNUhNJQRarjWUUqwUcoVw8Q74cnOWnGwm2jC4AA1STL5WHNWDLylbHOg2ZhzoYFdrskWqC3j8h2yjUxonrteQoRSXY0whF4/H3A02u5QTl8e4SHLSAHFx4dqQAqlNQ6jWhlRCCVKIug1Zso+OqIPXJmpMObuYigvZB+MHZ3N11toSXLLZ+CFk61IsIVefAzVy8OKjslQyYbAuZogusaScTRi8pFprscVmB5wkluq8lBSK54vq8NA7G+AGMXFwXpyDCSXhmI5DDPQqKboQXTRxKLGWkMVFidkiK4V1JtUXTzcmDaHkUpKrORCuTEOuTKmLNQUJJg+21grJS/TEMPLgayjZRWsz97WYPKRck1e+i57o9VCjEx9sdizZasrg0He1uFQdAaEyhOicSz4XV32kIOVYJdjik4JwbcWdTtwEUW8hRi7BOclE4QbnrC/VSq7VEWIbPKwVvYc1xJk6hORdkCrJIwXrEEN0ElJ2zlehzVhrSZ6h+0TkbUguZfFSo8sEMiiI0VXLEvHE24ZYk3cSRaqNBNjgPR9iyFwPRLLowJqpKVkuhiK8OPCvDY717In3dUbJxZcCe9bBhsy0Sqk212zKUBILOxSPQAhQJ+ZcKgNI2Xvo512ogeFVom5l4HahkBlYcrkyB67CaRZRFJnHFAGcz8kHJ8xjtj4Kk5m9R6narBS0qUgOFVZgQZQCR3gb4ZUkyTIFrvhY4KZQfPBibYw+oZwLMjL5ijwj3XmIUiSFnK0PPijP2lCJtFRbyTYLQ6lRQrKlukoAJgwRittqY7E5OhMGB+M7l2MphDjDYK1Eccm5EpJnZUllLVZJhdiqH0LIGXbKLpdCDSnJlUpc0RLZdUP1iaBU8dVnUh2HVHO24mFI7vRyA5PqJXuXCMK5AZEScpUcQyVJdrAxSarJp5xAQ4aSsreSi89ZI0ZDqraEEDOX9hCRG6JLyQYwj4Xw4+Az3CHRBV8Jrw3O+1hjlpRyIkw3IHJSTYhI8mHsUAut+FQcgTs7lJxSScXGkG2pYm5aVftI0mZcfcJm8mz9l+OtGWV/PN+tnm+5KQYrg43lq7PlmzuRXv6O0NH7GUS32Vg+WHK9126Ff/cJi6umGIPzKcBddogxoqtdUzNoEEdRCi7ahHJDi986Ozu08xk37RxuxtEXNP7XxW8u9JqiBf/8j4e4Ts6fm/X5zjzs49G697lp7H8ytPX2Eh/N7c355W55jmv5cvVsw/1wloDdD2zu3mpJKG1uSn5EU5+tn79cHrXlvoMa8+SMiycz6l/+fvX64rfcjfb73/729c0n4PDl7/79n59oC1/+7u2lv4MIR/ZgvzeqmYRYr8vtg+WWbQBXJJ3+aYU37Yil8J3P5O/XJd3evLzYXK7bhXuUvRH5fLD8eoWP653bJCKGR7Y0tu+7LWneHlnSytNvXK03b3SjLS5aQviYTwTJBAN4sJXsbYte0hvx2Pxik8bqQq2hhjzvkbk4zlHfFQk2YHhhaku7WGwQDPvHM1SAHrn+rt3zNkPVR7q0zmNs4Om7vq3+4/khuLI82XxztO+/vdniRVuerLmsiDHvrzh8uqW7w2aF2/zabuWdMZMeLJ3jJSB95DlsW5nmhCJ8eghu9q3z3c358xXrSK8lO9xE2AKe6/P/4l1c17bVd5e71et3iSFQf8uG/p+evVj/n7+67278+z/fOD85vfEvN3TPAphp6n1JMVPm2K+6DyPzsvumu1P6vXkWCt/lorLmnm6eVvjo/nb9fH2+97xyPWGJNtVSXMlvLKQfQppZXtz4lxs36pD8jX+6eHlDN2QJx+c/MmkQmTf4719uNK3UacPWtP6jE0fVyEwcF2PaE0fKL5o4+6v/VH9d/aWuJQ15t3h4k7Porm48vcVlz1v085Xgz4PldoXs/OviN9vV6WKU8t1bOtbuCCpunm+XL98wyT6WO++Sg1GFdpeEezfi7J/aFxwJqBYW/ZFq9Ypn8b1cVJ88W50RxiIhbXH6jiScnoHzNkfckcB/rdkfiCxiU10X/aKzdlBmT9f4MJsB9r3+vEetptpld5odC3HuLtHEPHEbMmaMUujl6sP15cXZ8ttDeJWg0J5RoaMmLB3uvby3Obm7fNrh70koer/Z+vbverYakf5rOVZHswUrz1e5ani3u15lEEIOP/9UvSMz6v2m6i9/11MFsXXF762knkX3Pelwb5sqprxP03/LFF1L5jrK9LLqkP/ofLddYzorb16+evnN0eO8rimeL/fV5/5i736fWfDgoHeNQ/EEzMkR9NBzsPpG77NXL+lDhct+Z4GM+GFxlr2OVQ3Xrrnt09cy/vXaW+3tg/U5GvSj7fY+F/LyHvj+16vt6dnmG3LEcAdst4iyq179wzEBse47JNZh1O8KS1wZ3EfvEebqES71+2uIQH9pBANnhrpmWhlOG9tcMR+9R6DLDnxAa3h3aItG+eGZMt6RSnw9k1jV+99/4FG3vz1k1WJYelPx9SDWDwk8OqW1TgDXLQ6OX0GdSxAu8aQ+NCODbQHD95gRN4jxQ0kmGxL2jbArwD2MX9Xhd+ek988ffrzYnMVju2eGj607JfaRdceiATxKRTxymvyNw48hEjHchx8LUAs/4uY4cje9X/gRP81R+DFzcX4u2TsXctLY5n++zo/Cj1eFxiEUmOZAYA4u4bhOsbgoxfjBhiguOF+zc3pEYLAcBLAh8IcGOCI/eOt9cYUoFxdX+EG8cylWjbJw4nyIEpMvhFscJwKOopnvxGiuRIDzeyu9Z5DzJmh4wpWZwxKpmpvihxCt+JhTtvj/b4ol6iXFSrWZs6M3y+A5LiFSaq0EQG6mwdlYBfLn4kM2NwM+bWdtLLHGFIO56XSegsOnX33J5qYMwUfJPtsixXp1WRcXbLQES2pVMZpKDjaHWBNnWnCFh+pLlVqjJA7quEFiDNXlEFItmcBeKjVacTGrgxy/P07zXFwqqbKIgydCQUTNclDGD6ngi3e5ViGK6ocSCdYWgpkcw/FDzUG8T+LFExsMygSEpkLKHJwKg+Vr8eKKJ94VBlsI5lZihZHziIPNqSRbeB8jsT0bQ/HZxpqDbwWueJdLtbVk0U+sd8GmAsUriBHI88Vam2yuURHLLlVLkCJFBltLSJYjVwR8tRdbS7TRJx9jK8DFHmqojpiDCYM4/nBGTsEHTrCFgehWxhtpXUPMM24hnOqCo0bIhJ5TgqotTkkoMXL+SUKhUcgZGJktrhAEIp6DZz9kl3MhcJljjc7nkEOJGrj0HNCOKZXACcA4sKYYjBTLEco0wDrOxuILB4rSQCS+hFIZahCiVZwVi7HYkFIlBGYJ5UVnYXEusRl8LMSWYw7Csco8pJKyEDyuidNbZYD/xZXoovdCXNLTpAaRss/km6RQeVeIvwUxRYOdpDC4mjmyVAcXLOFYBsLhyjpEbneAzMGyQuqQYyXeGVO0yRZTCelIIDBXIjFCsYMUX6vzVgicCzk0PjnrUopBD+cNoUriL3/EwIVHaKSaLMxNBgInwIaUQpUgRSCrFtTERBEu0kN1hKdicY6DSXoQb0ghhpJiyS7p0bkhJvK+PK5H7ST4nIlpRe9bmwTlapBAUgN9EBSuKSec4gy1BIkB5yUHvDWCml0NMTuEKnftDC57QnO5SOY4VhlqyDYU8dYTgDNliJUgXEnR5QzsXHWFdRNtDEw0YfpCPJiBZqKUjtCoZwn2uKazHPvKPhMDTwRora0kQSTCmMTcQkWT+MT5af4qiQ+OKGPmEFccsiXO56ol3gtLB+YvWSQPp81g6UAZDtrgqSGkEggCij8+QkEI0VaXiWppDecqwc/CSqiORn0OlehrikSCIykUTgL5GcmLN6mFQkkYgU9ZFjWXEIMjkF2EsG5J8HeBXTxXMA01p+xjzTiPlbe8j6F6h1jV43pDTfz1FYr1BL8fCLOW6CVnPaKYBjIRSiXuy2KUMuScfCTeaKPemDPk4lypJYYAwVwcHEMV510NORtXNVRSqw+WOIrxngB9ZXmTyJKNL4NPwYXivCXDJviheI6T+uqy9w4lgT3BVjqr1vxJ45qJeDlGSvujH7NpNbuYZviNwOY7tXMzKohde0LENgcYieQU1purOTr+jo/Pg5XoLFkvgZOPJQ3FuxJykRA473rT+zKgnGrJVrLlhqI8ZJciMp3Fp9T5aQOlF6RJz3HNHxMivSDVbm7oxwRIL0gTnxtyPyLSeuGPGvI/piFSI2eMwo9pCCt/biheC/z+9smF/fIPTy7ky9+//sOTC/fl7x/eW+5ejONnf/7t69/94cmF//L3F5tvfvvaeKAwQwEozlB8VwSY0HzbMfx/FAG2Q86VtAKbSOeQOQBchpJ98JlcsUpS1vx3wyRS0aIpuSSg/92w7wn/6p/K2od/s2Zh7gO+7MF+geFf3Uoewp7vDv86fJ1/n+Ff9owa/tUrJn4x4d/gMWGTlOyqqJS88icRf0j898JqiJMgpyqlHsbLweVfVojzp6CN7GmDgm4hTt39/0oatycNpopTthH8IL+Sxu9Jo0abkuYGLqFfSRP2pGnmq9JGnWO/0ibuaaOGfFtRaYjuF02afajnaibFtcQKBXW/8lNkV3j/a3YFNuzbr8d4vyDwP1B2xdvCFT9/yP4d17W832z9Y2RX2MGHo4NtMhTdpf78c/VreoU6xI5uDZjzljTup3/A/BAB5OlnT1b62dIrGg/uI5Azj/6aXzHfHrTenPfkvPfMryDnQTgH0lIoiOC33ArNstB8Cy3ToH9QK+L9Eiz4ak6yoC1OT/JzJcniyyPbgZSXhXbwt8h7+WGZnZok8ZYcd83YO9wG8Gs+ZzvS/Pd0C9tR2pmmdO3zOfthkn2GBDcy/3cIyV9tjv9Ujx1diNGO0+9vwfj5ddqvZsd/Ol3/MGYHZ0uvZDjuL275x7A8rt8vsTfANNVTdThXgK0vOLyoyvLj9fMXZ+vnL3a3N+fnq2e7Qyb8nfXr1UnLRz1dnl2uvvt/sRbVZQ==').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779223073823', 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_1779223073823();\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
}
