{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "6c251746",
   "metadata": {},
   "source": [
    "# rf305_condcorrprod\n",
    "Multidimensional models: multi-dimensional pdfs with conditional pdfs in product\n",
    "\n",
    " `pdf = gauss(x,f(y),sx | y ) * gauss(y,ms,sx)`    with `f(y) = a0 + a1*y`\n",
    "\n",
    "\n",
    "\n",
    "\n",
    "**Author:** Wouter Verkerke  \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:30 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f80a7927",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:43.856960Z",
     "iopub.status.busy": "2026-05-19T20:30:43.856850Z",
     "iopub.status.idle": "2026-05-19T20:30:43.870269Z",
     "shell.execute_reply": "2026-05-19T20:30:43.869800Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "#include \"RooRealVar.h\"\n",
    "#include \"RooDataSet.h\"\n",
    "#include \"RooGaussian.h\"\n",
    "#include \"RooPolyVar.h\"\n",
    "#include \"RooProdPdf.h\"\n",
    "#include \"RooPlot.h\"\n",
    "#include \"TCanvas.h\"\n",
    "#include \"TAxis.h\"\n",
    "#include \"TH1.h\"\n",
    "using namespace RooFit;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c367d677",
   "metadata": {},
   "source": [
    "Create conditional pdf gx(x|y)\n",
    "-----------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "58b709ca",
   "metadata": {},
   "source": [
    "Create observables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "cbc6c5e3",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:43.871771Z",
     "iopub.status.busy": "2026-05-19T20:30:43.871654Z",
     "iopub.status.idle": "2026-05-19T20:30:44.083457Z",
     "shell.execute_reply": "2026-05-19T20:30:44.082772Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar x(\"x\", \"x\", -5, 5);\n",
    "RooRealVar y(\"y\", \"y\", -5, 5);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7217bd61",
   "metadata": {},
   "source": [
    "Create function f(y) = a0 + a1*y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "49f7db54",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:44.084885Z",
     "iopub.status.busy": "2026-05-19T20:30:44.084767Z",
     "iopub.status.idle": "2026-05-19T20:30:44.292614Z",
     "shell.execute_reply": "2026-05-19T20:30:44.291938Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar a0(\"a0\", \"a0\", -0.5, -5, 5);\n",
    "RooRealVar a1(\"a1\", \"a1\", -0.5, -1, 1);\n",
    "RooPolyVar fy(\"fy\", \"fy\", y, RooArgSet(a0, a1));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a9cba1ea",
   "metadata": {},
   "source": [
    "Create gaussx(x,f(y),sx)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "12f47dcb",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:44.294110Z",
     "iopub.status.busy": "2026-05-19T20:30:44.293991Z",
     "iopub.status.idle": "2026-05-19T20:30:44.501939Z",
     "shell.execute_reply": "2026-05-19T20:30:44.501603Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[#0] WARNING:InputArguments -- The parameter 'sigma' with range [-inf, inf] of the RooGaussian 'gaussx' exceeds the safe range of (0, inf). Advise to limit its range.\n"
     ]
    }
   ],
   "source": [
    "RooRealVar sigmax(\"sigma\", \"width of gaussian\", 0.5);\n",
    "RooGaussian gaussx(\"gaussx\", \"Gaussian in x with shifting mean in y\", x, fy, sigmax);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1824189a",
   "metadata": {},
   "source": [
    "Create pdf gy(y)\n",
    "-----------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "56b547c6",
   "metadata": {},
   "source": [
    "Create gaussy(y,0,5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "8fd60fc6",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:44.517997Z",
     "iopub.status.busy": "2026-05-19T20:30:44.517855Z",
     "iopub.status.idle": "2026-05-19T20:30:44.733234Z",
     "shell.execute_reply": "2026-05-19T20:30:44.732508Z"
    }
   },
   "outputs": [],
   "source": [
    "RooGaussian gaussy(\"gaussy\", \"Gaussian in y\", y, 0.0, 3.0);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4622cf94",
   "metadata": {},
   "source": [
    "Create product gx(x|y)*gy(y)\n",
    "-------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a2d6cd48",
   "metadata": {},
   "source": [
    "Create gaussx(x,sx|y) * gaussy(y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "15c961a3",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:44.734946Z",
     "iopub.status.busy": "2026-05-19T20:30:44.734826Z",
     "iopub.status.idle": "2026-05-19T20:30:44.942619Z",
     "shell.execute_reply": "2026-05-19T20:30:44.942003Z"
    }
   },
   "outputs": [],
   "source": [
    "RooProdPdf model(\"model\", \"gaussx(x|y)*gaussy(y)\", gaussy, Conditional(gaussx, x));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ed57b4a9",
   "metadata": {},
   "source": [
    "Sample, fit and plot product pdf\n",
    "---------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "91b552f7",
   "metadata": {},
   "source": [
    "Generate 1000 events in x and y from model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "fa3fa5a1",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:44.944685Z",
     "iopub.status.busy": "2026-05-19T20:30:44.944547Z",
     "iopub.status.idle": "2026-05-19T20:30:45.165924Z",
     "shell.execute_reply": "2026-05-19T20:30:45.165308Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_54:2:2: warning: 'data' shadows a declaration with the same name in the 'std' namespace; use '::data' to reference this declaration\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, y}, 10000)};\n",
      " ^\n"
     ]
    }
   ],
   "source": [
    "std::unique_ptr<RooDataSet> data{model.generate({x, y}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "188e0fc2",
   "metadata": {},
   "source": [
    "Plot x distribution of data and projection of model on x = Int(dy) model(x,y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "0a66aa64",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:45.167494Z",
     "iopub.status.busy": "2026-05-19T20:30:45.167375Z",
     "iopub.status.idle": "2026-05-19T20:30:45.376288Z",
     "shell.execute_reply": "2026-05-19T20:30:45.375632Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_55:3:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(xframe);\n",
      "^\n",
      "input_line_54:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, y}, 10000)};\n",
      "                             ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:344:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(initializer_list<_Tp> __il) noexcept\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:312:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(_Container& __cont) noexcept(noexcept(__cont.data()))\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:323:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(const _Container& __cont) noexcept(noexcept(__cont.data()))\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:334:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(_Tp (&__array)[_Nm]) noexcept\n",
      "    ^\n"
     ]
    }
   ],
   "source": [
    "RooPlot *xframe = x.frame();\n",
    "data->plotOn(xframe);\n",
    "model.plotOn(xframe);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a7c1845a",
   "metadata": {},
   "source": [
    "Plot x distribution of data and projection of model on y = Int(dx) model(x,y)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "306ea71d",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:45.377954Z",
     "iopub.status.busy": "2026-05-19T20:30:45.377831Z",
     "iopub.status.idle": "2026-05-19T20:30:45.586523Z",
     "shell.execute_reply": "2026-05-19T20:30:45.585939Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_56:3:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(yframe);\n",
      "^\n",
      "input_line_54:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, y}, 10000)};\n",
      "                             ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:344:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(initializer_list<_Tp> __il) noexcept\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:312:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(_Container& __cont) noexcept(noexcept(__cont.data()))\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:323:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(const _Container& __cont) noexcept(noexcept(__cont.data()))\n",
      "    ^\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/bits/range_access.h:334:5: note: candidate found by name lookup is 'std::data'\n",
      "    data(_Tp (&__array)[_Nm]) noexcept\n",
      "    ^\n"
     ]
    }
   ],
   "source": [
    "RooPlot *yframe = y.frame();\n",
    "data->plotOn(yframe);\n",
    "model.plotOn(yframe);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f20f3e10",
   "metadata": {},
   "source": [
    "Make two-dimensional plot in x vs y"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "65d108f5",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:45.588186Z",
     "iopub.status.busy": "2026-05-19T20:30:45.588065Z",
     "iopub.status.idle": "2026-05-19T20:30:45.795597Z",
     "shell.execute_reply": "2026-05-19T20:30:45.795125Z"
    }
   },
   "outputs": [],
   "source": [
    "TH1 *hh_model = model.createHistogram(\"hh_model\", x, Binning(50), YVar(y, Binning(50)));\n",
    "hh_model->SetLineColor(kBlue);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2ef09515",
   "metadata": {},
   "source": [
    "Make canvas and draw RooPlots"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "56325b51",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:45.819182Z",
     "iopub.status.busy": "2026-05-19T20:30:45.819048Z",
     "iopub.status.idle": "2026-05-19T20:30:46.024742Z",
     "shell.execute_reply": "2026-05-19T20:30:46.024159Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_59:2:3: error: use of undeclared identifier 'xframe'\n",
      " (xframe->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "  ^\n",
      "Error in <HandleInterpreterException>: Error evaluating expression (xframe->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "Execution of your code was aborted.\n"
     ]
    }
   ],
   "source": [
    "TCanvas *c = new TCanvas(\"rf305_condcorrprod\", \"rf05_condcorrprod\", 1200, 400);\n",
    "c->Divide(3);\n",
    "c->cd(1);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "xframe->GetYaxis()->SetTitleOffset(1.6);\n",
    "xframe->Draw();\n",
    "c->cd(2);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "yframe->GetYaxis()->SetTitleOffset(1.6);\n",
    "yframe->Draw();\n",
    "c->cd(3);\n",
    "gPad->SetLeftMargin(0.20);\n",
    "hh_model->GetZaxis()->SetTitleOffset(2.5);\n",
    "hh_model->Draw(\"surf\");"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "703310aa",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "dc96bb18",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:30:46.026401Z",
     "iopub.status.busy": "2026-05-19T20:30:46.026277Z",
     "iopub.status.idle": "2026-05-19T20:30:46.238875Z",
     "shell.execute_reply": "2026-05-19T20:30:46.238523Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779222646196\" style=\"width: 1200px; height: 400px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779222646196() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(22412,'WkwIUhwAjFcAeAHtnG1zHDeS57+KomNfYhhIPBVQeEVJ1nruZIthySNpJyYmWmRT6hPVzW22bHk2/N0vfglUdZOU9jy74ztPnB1uqhIPlQkg8c9EAqj/WPx1/9P1arP8sFqMixePlpsfljcvV2+eb5bXN++2+4VZXH6/Wf/7x9UfHy9GaxaXD9f7m/b07M3/Wp3vSV9Q7Nn1fr3ddOJ/rjcXi9GbxeX8pvE/PsfrSwy8D86HYBaXT9eb1aPt1Xa3GKWTz/c/Xa0O5Mv1xf5dI5+sr656YYSFnApbq9VXl/tvlru3681itCekfLd+++5O0sPtfr/9cLvYi+317YRXl2uEcGZx+frweNoeefGrm/1yD5dSKHOLOm0UlZ/slh9Wd+Um7U7D53K3GzQXnZKnV1L9uGd45cPt7mK1e77+W++9o8RvtherNq6vpP37uv/7yrW+fd3/fbXfnr65OVt/Wl297zX228+Ri1GkJNp+q4YftM8OdW4nLMY/tITvb9ViOL8/1DkmZz5/ulWjveVPhzq3E2Y+U2P221e9PdqY++RizCc+ibPBivfDUPzqD2jo4QWv37eu6i+ARNCJXIx/sCfWupSzDM4Gb8UOEVW52v747eNHvd+PiVffX2sGw/r66PnlnPr1/HT65ubWi07f3Nx61+mbm0O10zc3h5rff/qgE4I+/unw+OnD8lNr0Pc/zY8v3q32y8XotWHv1v3p9OZ6db7/brlfb1srvv344c1q155frM/ffzo8/tQen27f9sSn27eHtL+13LPlxdlyvWECmcXlo9325ubdct1fOJNn2w5Gx5qNejX6oNTfbC/Wl+vVxWK8XF7drMzi8l9364tPt8mfDuTpm5tH2+3uqPxXF+v98g1Tfr/7yAuerD+tLm61e3r12W79Yb1f/7C6uYd6T9c3gOoEuJ1c7naL8c9/MYvt9Z6Hn83i8qtPq/Obxbj5eHVlFpffNoTeXXob/3q+3Vycb3e76932AgB9sd4j12J3+ZnMbz9+OFterfb7CUXpz29Xn/b3Ux//8fnZ09PXi3HxL9OjWVw+3n58c7V6+PHychrP71b75XpDZ/aueHWz/tvq+5sp//VtUnO/Wy2vFqODuWY3Op0k/hvM4vLlenOx/fHF9vpVG/KZfn1Md0TjNa3A1yvwu2vJjzMUPHq36PP90XK/vzcMp/t9s3U08NXD1f7H1WrTwfwWpV37ZLf98GJ7vRjlBN16dbHcg5VKvJ4IjN1pI+Rns3j/zfaH1bPr5b9/nDXm/XcrOuZ24uXX67fvntKGbrdUe5f783dT775//m7741c/rDb75/vl/uPNrKTvTz/ut6jJXPKb1ebjw+Wu0SjR6TlqONe4/G61vHi2ufppqnH5cr1/t/24P9bYSYu/Xt50HZxSjkv9+Y45/4c5DSDqF52Gl6s3igrrzdsveQ4KF1fLm5s+ZSjXXJXjhGt0dWFHF6PpvyqjNdbY6jSVJz9azbU1aJ6Lsca5jq2plyR96GV5zmMuxomYHGrhOQfjJFWBYTD9V0VGKc5IdkYGV8WNQzL6fxU/inWm/6qEUdxg+q9KHCV6039V0ihDNv1XZRidDab/quTRuWL6r0oZXXCm/6qzWlgK9W11cpt0oxRvpIiRlKvzo+RkJEP66kLjm6KRaKuLo2RPjpEYq0OqLmSw1Q2jxGgkOCPeVpdH8VEbmFN1ZZQhGmcTklRvR/HOKOvkq5dR/IAERpKt3o3irZGYDS/2vvHJudWlq1zjE331cRQZjNYPsfo0ig1IoCL7YRRbjFbwrvo8igvar86W6suh12OswY4yWKPNLVKDjAyxyuxjDW50Toz2DyRdlZuQ3tYQRu1iRs2FGuKofYoYUmpII68pwZRUw9DY0FFhqCGPylFUzBrKKGEw2ghxNVoVQlWsRhklJKM6W6PTfuiEH0uensM4hOk5jlGm5zS6YXoeRj89ZhS7v6aM1gQUONZkR2tSbs9MnSLt2Y00SqdVTcwg8aXlMInQDq0eIfLQiAQPaUySziTfuSh3naRSE+yVcEMd4K+EDXVAAAjJtg5IoERMdZjmMOM7IIHm2FwHJIDIsQ4qQIwmSR2mmeyHOkzsxdehjE1CY2u2YxOQZxmjNoNnN2btBp6Zwg1dbM1opXYDRJu+8LM1M03SRDB1m1S2ZiaudgME01Z729haOnD5DCENkSRANPCSAnwVrzmSEkRoRHAQDcVQKltLg7Gibx601CCkZ30OA89Fn503toq1o4iYEBUUxIpqaHAmuCrWNVAoDK9UsX7MgxHVGrKDtkNiMh4yjs5l4wcTQxWbRonWDEzQUsUCGUnLSuJVTE9nvDfOQdIlJtGrUkWxNTKHXBUBxsR4MQ7Kaa8CRqlUAVt55udiFcCV0aTjea/EJnHTZwFdrZhSjFhbRYYRgKM9kquAruhdKkZA9DJKLjCWKs6OEsCsYCSB7+BY67cM6UatZgeTSxXnR+eVj60CtOaIWTBMcwFbxRuFMkfhpFPdSW4yg65HJsBZGYswJYtUx5AAie1XnQWbvOm/6hgToT36q87GhrfAenUMSvYKAbY6O4zMzUHMwJvz6HxpYslQnS0jtQDrWB1DghZQuTqRcRDTKlcnTnG5166OQYEflauTQM1WuTqJ1GyVq2M8QHQZqF2dDE3aiWmeG6dci8o7scXcAUdN6Op0SJC28WVIWlsl++p0SLSxzofqXJhbq5wZlN5ahHag+KFX3TC1t1XOc4Op7Kc5DLbZUP00jVEm7BYzuVvxRrfJrPnWVc98nvKVPvIJlG6zmvJRqmdi9+JK5tmDUBJsa+5GlBrswTGyNdjmE1Adys01oQ5+ClTAnGBRarD0RkMv2hd0Miv2a3sCkxkaaxx9DaiOdXP7gthGk29dDThLiIgFVdrBqzcviFfKBpVfVA6K0hpRI0dlpSaotzWIQn1rlwD1VCEdQ0MFW4PDzrQmuT4kzWesAXeoNwmTFpi0vUnaZPSjN6nlgyROm9zy8T56k+gCl6cmteJlapOW9hYSBlrYi5Jt1Grwk9HTql6tXhvvGrwa3rmiWt65Ht2hg1YjKtl1hLfEIvNQN/ow9EgUy2HwW35o7WltrBGVbHxU5FgSIrdRq7EME6Xv0ra3BtZYtOlQ1tVkGQK1ijXZ5qVTJ1narINWk6XFOmg12eat4zGghW2GqYgJ+GozrNNZm9yHrSZsCvlt2GpSCFMNbCzRwqaB2qSEadEZ1l4nXl9HExEQy6JK2Au3/piaJcfLh9S9dlCPmm129MZhVdDG1h1OFYFRqwmVnNpHPRCrtw8tSUBWm2Gd7oa/D39SzGrt065WL11nWKebeWsKklyZmgdvvPTWOqXa/Owl/WERlXzrFDyShHfOHNPn1h3tWVWDx2ECwiZ2HToQ9m6vA1rXRYxSB5Suk1SOc3s18zD8NG/o7o2KUIfu4MBJC7fWtcw8YaCiQwYDe79Gqdk2JNBBqbmDoIJEtvOYk9PFadiTuzZ2eWuetFGRJqOLXfeUx9zZ+JgNkRuLroQUpiCmtMMKJRUJ+ysPQEhOh6PW3IwdnQDZupqxo03ttXCeKeUxYxFOriog7CG6dWi9UzMmdMJ1XgtIom0NjDMY2Ukqd0+4i9S1sfdCzV0bJ05trkIhEmvGzghnu2NXEwpl7G1jWmSWjL2LmBbZt05qelPzASCR2Ku9QKZWtqN2N9TZT+uDVviA2cpIlyhT1XC8RslBEZv36HuDInabfDUHXSl1PQkTYNPOcBAH6cJBnPaaWRwaFmZxWuFJHF4UJ2mU6D3UTHTNLBf7uGhD4qxG+t6IRB2dc1s19g6tOTaJeofXHJFIcwmFRJXoQCJRx8esa8huD3JqBq0NW81tHdkGrWZdSHbgzG0h2aStOTXz3lSh5nRkZOGvy8nDW1v/TK3WFeVctS0pp7yhh0ImgYZu/SaRhqbYU6cMTbEbYuehKfbU0KE5Q3MvDKxsphV2HppiTx3BQrNlolG5KXYzfzm3YZuK5jZqzWjkrGPWrETO6gN1cVhxqm61wWXN2TFT+yi3qT/xaFN/5qFTv7Mo3fZotY7UUw90oO78u3OgZip3kO5EA8VpNDtGz+PZMXp6zWyAVKE7RuuLSofoThwQmn4rRxCNBhdAutti62rpMN24lI7TvSVlwuk+1KXjdOvM0n2GiWoeQx/p0oG6Z/bomuoUQrHcoOPbrxbAEj1uv1qcjmCzocWFMXqW1eKpqmoeWBoPtbg0Rqsr8MCifxhDYaEugZfmMWTDwjOkWlwZQzLRGgm5Fm/HAHwaiVKLlzEEo8KEWrwbgzcKT0Mt3o/BGZbfiUhCGIM1RNySq8XH0RcTC7GwWnwafTaphQiLH0Y/GJbIAxzy6KNhATzAoRA2ScnIkGoJdvTepMHIUGoJMnrHap61cQlu9LrkkxxrCX50xbAWzkMtIYwua5CgENaIhK2GSLSwFiJpyQyDkQKHgQjBUIwUOLBYNRlNllpCGZ0zufl1ReNoRpd+Qy2AYzE5EpOqBWjMJifjxNcSiR6YnI2TWEsMoyRTiF7lWqIuCAmEOakFXIymuE4x5TU84Xgn/h0hl0bpsrmkRiX8KlOGTqEyprDOJ1RDzIBYhJNSS/JQQuBCSRVGtVbJ2HJ9z02NbIGLkoZG0shcS8o0SyythCQgawhNKDkQ+jSESp0MtQyiJOEwJQn5GEJCjdQeEnGdJEpsiF4QAi8gYzaisQZIogstXqjk0Ej4ksti3ghRDnobZIREjFhLtkoSKVBSGokY5BK71qBHI30jESPUQkCOXJoPSUDWiKP5kKmRiAE5NBIxfC0AJIURA7Io6RHD11KaVMSexNXCIioTfexk6yviTSK1gJJEpxEDso2gRwzINoIeMYjBqUKJRwxI1SiNriqpKiUB/S61sKQPRohiWWJpBO2gUfgMTTAKwJhooupGAuEIzfejEIFOxtmB8kQwiUhPdGx0nmjCMkZCISJO+UHpiDhKY1WAHeNsJF/DqRKRB1rsWLopVlKU1G0Mct1IML3FSsSKVxJZNDcoiSie3DgS0kESJRNkQhAijgTwiCdOZFYSMTS3jETz6BVCjM5qXY0OQgqMgLNCnzqHkGx/FLpMXV0Nt0qhy/B1XadpE94uMUjQCbE1iNdpBFV3t21zFGQhZITk4JcK09RfMRL2LMAK2zZGMvw9Hkun4c8ijK0l3ZWAJsLYYooZeTyx2E4jj2dx2OKixCet1x0BjZO2SO4A8Arh2Uxv+Awss3EDgosF3uEPjTwAPPwb3osF4uEPDX9AHv5YAPgHj7FocUz6IwRMSaORJ0QMjYDnA/xDwgxJps/hHwaMVKOVf8aEtd0q+iMUDFyj6Y9oMX+6/ZaQJwrGUcOiCf7RYUgbDf/uE+seF+MRwxipj0GHP37xYCSzcIB/TCM2UmnlN4wJechXfnlEd5Sm/bGMbHER99VQebJEEBsNvyQjMWvyA/ySbhg2Gn7Jj/QNO3QBfimMbC0qDb8UR3gpDb+UNNCtNO1Lw6hRe/wI5ZdHNv+gNRyfyljY/YOG32DZVWo0/YtZsAgM3sEQw8ACRxPgqD4zXaSISEIYMSMtAZ6YByZNLrpNKBYDQZSVBKdcARbYkqBs2RaALQnKln075qjVzTaxGArH1iIJjCymAtXVBNhiLHTHkRLoFuYCZWTnkzi/xWCgjZQQ2GIyUL9DAjF5di0pQf9jNlAwTUAwDEfocrBbYDEdqJiWQDCMB96ZJtBBUwxOE5BUg3C9PwRJNR7S+0ORXo0IfcqOBaKrGel9qliOIVHNowSSYkqAJU1AUoxJwxm2GknAIehjS4JuAikSZCPArG4DMbU1gV2XaSdIE3SHhn2H6R1spWBTdCMiG/CU3SDVZmj2VnTrgaZlA6DqhhAdkQc2SkU3hHQDfTAALDtCUtC3wQCw7AlJoVmDAWDZFhL2caB1l4gwDgJNdItUKY3Aug1Bk3DxoFm5IE/SHRu2hti1pDwAzN4QW8JKIx+bQ6rMmU1GYXeITRbaA+CyQeRUt7PJyEcsBcRnOxf+RPe6qgO4ukUEwhdrAFx2iXTPGRr+RPdA/CJmoL+Ip+hEEAPgskvkQHho5c8OIzrtDIArxPdAeGj6h6CKThNvAFzBwLC1VbxpG2ludDprggFwRQ0MCh8MgCsYGJ1EOP3QcXQe/WdJAJ1G9jCVhj8GRqdUMgCuYGAC/BIbsyK+nV1gTAFcCXZ0OsEGA+AKW/IgODT82ZPX+ZYNACysI3S6ZQMACyuJoKsDE+HPWkJnXzEAsLCaUAdFl5girCc0zmBNpD9YUSjtDIAsrCmU9ibAn1WF0tEAyBJZo+JaJAMgSyRKqg6YAZAlEr3FEykmIE9kzxbHUXRrVSLyQXsDQAtb9koHA0BLRD7ykwGgJSIf9MBmtkhEPuhiPPIk5GO1Yg14LQn5oJ0BriUhH3QwHnkS8kFHA1hLQj7owXj6I7X+Q4fYpJSEfNY4XGPkScgHLQbkloR80N44lQf5oNk2Z6sV+aCTAbaFkIzSrIagkY98PHFoxrcYdAYQl6GNL2ddBHmG2POZg9Cp03gk0KwXqY/XBY3+QbO5y0Gihx8vOUv2L3pgb/Hkarvce7cwiys9PRajWfywGP9cXDDFsfRLpjiwJZviiinemsIWr3emeG+KD6b4aIpPpngwKpviiynBmhLElOBMCd6UEEwJLB45IwK2ZVNCMSVaU6KYEp0p0ZsSgynsgcdkSgQTsymxmJKsKUlMSc6U5E1JHDSJpqRkSgJLsympmDJYUwYxZXCmDN6UIZgysBRNpgyDKQM4XEzJ1pQspmRnSvam5GBKZpGaTAFLFb+LKcWawh50caYUbwoHXEo0pbCEHUwp4L7uqxOHAYmtRkIsOxAaGbXsxQDBFrfBAr6WAzsW2LVgrQVgLX61BVoteGp1VQySWuDT4iVbgNPiJ1gg0uIJW8DRgogW98CChRYAtKCexZm14J0F5Cy+qwXerB5cAsgsHoEFwiy4ZXE/LYhlUTmLt2lRNgsqWdTMtsM71ECrLIEZC+bYSA2NyuFUWmy8Jfxi8SktFt7iTVpMuyXaYrHplkiLZUViseiW4IrFlFtcSYsNtywQLMbbEkKx+JAW083hYxws/lCDAwQWD9LiUltMtcWBtNhoi+to8T0tDqPFC7Rq2tVpwCxbnCWLQbb4jRZTbAs1cBstltjiMFrGnKgFTgN/OOTAmLdTTYw5gQjsPX9Y8TPm6igSZRD1ENt5CMZc/UNhzNvBAMZcz0sRGsCq8kcPUsCDMWfdj3nkDzwYc3UI24ksxlzdQWHMWaxju/hDDcacRTj2iEAGPBhzdf5YUGNY+EMNxlwPbKmfx/IYs8AfajDmGgwUxlzPx+mJF2HM9YAby1iAmj/UYMz1hJ2uQtvJMMZc3Tc9Iad+my4K1WFrR0gYc3XXdHWmfpquu2TIf/n555/Nr3Wak9PtXzzN2a4s/CdXTPqJTY5x7zYUnE46t5QH8wvu3JWAPLoPcfcqxOGSiHA14/4tkW+Wu/er3dGtk5Zw9MqeMN+keLH6tD/dvOUYNgdQIVumPbH0geZfrd9uOKHc6aP3k/1ky2n3pCeMl5/W98+On+73p6Rzfvti/cP6Zr3d3CzGKHAk5+iFT5dvVtNNGPgp3TgEOCj97PLyZqVXVMDZnjiL7VXu9fn7p6vNWy7Y2BPLIWcdg6mqtoXj53er7a+m8+ZzkYk9B5Nf/9O0UIfzv9DCf/unaeE8QH/nGD5c7o4uGD1c7ial0CPYXNliml6dPW9z4vFu+WO7ltHoZ9f7wxWQRvRbII3oF0GeXe8ft/P27ZIZh+CZRDqNnl3vOyTQiGfX+yd6KasXfbLuFwPunaCngCZerPfcLJvoF9vtlR6gJ6FdYXm03ey3H3c3/fLC6b6LcwcxT/d7JrGC1H+CBe7vBAPmCs3nUsSJ3mCD4kpEo77aXHy12237nS5mtpJaHFZPPm7OOyyQCXmEYpB9CMnlSkwvTPshe2HmPeTReD9dvV1tLo7v3TAELfUIYXnRIXHiPV3G4xUTLMwF6b+uimZx+TV3IVY3dzC8pz6/Xp5zF0B5zzfkjtowX4/racg4l7stzVx0Sp6K3mGt5e62+yjxcOno6/UNCnksD0m8r4uTLM2eyk2MW+9MRafUXvCONJT6Zr1Zf/j44d9Wu+3hqgcZty4tKsS3Wy9nu9XlavevTw+lW/pRx7WE42Yi6XHqoZ0t9fHq8uvFGC2DNqe8XIzD7ZRXC/YBjoq87glny2PlO1ve0i2Yz0kHzpp0/7Lm2fLiVtPpu7Plxf3Ln2fLi8/c/zxbXqDsrw7d01Ne30rBNvbLSjBcn7/vV5XOltftluWrDhpzwusF25OLy+fnu9Vq82R5ruiDeMDaUfdDMg+O1Jak4/GYah3NH4pAHmYPKZMCNbVqZXYfQNFFOgl6kYZExYySuywKMSxOm2gvp/tOlPwaQrQd+936+vHqfP1heXUzXypSSO7ejptdhqPWaYE7zdO04/ahJJp41MCZnk1Vc1JwQNql4LnWDIFapzWnzKVpz/R+mtPAZ/vHzWa1+47mUZKppq+9WYx/5iLRgwf8EfdAWPQ/CA/yRKcH4YEkKGcpEHqRA6VJx+nO3i4+52Xb87J9EOyD8CBoSeX9K/35C+CyWl6sdlhrvTOl3TZTT9b7J5PSxK40ehOLUZxzVJe0186XV1qZ0f8f2/WGxMkReLS8PiZfrD/M7uSQc5EcFDH++GH5dsWLZoB/tNxcXK1evlvfvF/tvltu3vZrzC394fZTT2uj11JVkqNrnH9ab6/Wmym1319sRR+td+dXd9G+Z3H/FKGPDOArXO6vPl2/OnZ7psTXx4mvP1dySrxVkoLfLD89Xr/Va/Yo4bPd/t320fLDarfs6HN/uXa2vPiH3b9jxL64Yjtb6oXTz30SgJUUPXSEWJAT+PTco/67Yxnvr8LuwDWd8Rmw/v2yPuD1yy7ri5yUWEopg3dp8BjVL97d90M4GQKXc50lFnLvGn/KJy5b/sue4IYWmL4D4FM4iYlMSckVLvhON/65eH8/9/Zt/88I+sXL//cFvfMdgPuCHn8U4L4oU+7nBT3c9udzAVzo91Ik5OJyzOIDZvGLHxBo1/9TDqn45CQFQsZHr9QPCJxYgpXBZvGF0On8vi99T2AI3nqXSkwlB+bh4TMAJ9Yy3gz5kEJ2Ucdp/sLAZ7Onbw7YE+98+y/GlCNzePoGAV+10LG31jkfRW9wt08YnHjptYZsfYJh+9KAPSlNX6wNYWjxhVvfLPisNLe+YvClEl9mPX/n4HPcf8XvHtBZv3/3YP6yzN/93YM7y+tfXv+LX0r4qxxHEO9/SaHl/9LPJfzjbtd3H+TW5fq/3IvI/n9h4n+38NM3eg5L3L/vczxB7ElyfFPDF7YXf6sG/r6cv1H7Lie2xCGENGT+ORjjux8L0pCcdf9XjbsP3dKyr6YL8l9s2lPUL6+klIvlBNqvbdrvi/pPatjxx3437L8pw04wpsUr9CNIdz+R1PJ/N+z/D9fuvxv2/65hz7acuALc25Q4EfNbteyfEfQ3atrdicRm24MrnJmZF9q/Bduehm6gQ/AcG1n8ctt+f2n+6y7b74v6T2rb217F7x8r7J+D/eWL7v6huF9h0e7/D7Zd838rtp1tlOfnu/U138jVXRS+7nfVvu632azO94cNRP18ZNvG0u/r/fy/AXyuZcg=').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779222646196', 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_1779222646196();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gROOT->GetListOfCanvases()->Draw()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "ROOT C++",
   "language": "c++",
   "name": "root"
  },
  "language_info": {
   "codemirror_mode": "text/x-c++src",
   "file_extension": ".C",
   "mimetype": " text/x-c++src",
   "name": "c++"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}