{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "b58d8771",
   "metadata": {},
   "source": [
    "# rf108_plotbinning\n",
    "Basic functionality: plotting unbinned data with alternate and variable binnings\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:29 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "a82aeace",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:11.343844Z",
     "iopub.status.busy": "2026-05-19T20:29:11.343731Z",
     "iopub.status.idle": "2026-05-19T20:29:11.356653Z",
     "shell.execute_reply": "2026-05-19T20:29:11.356061Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "#include \"RooRealVar.h\"\n",
    "#include \"RooDataSet.h\"\n",
    "#include \"RooGaussModel.h\"\n",
    "#include \"RooDecay.h\"\n",
    "#include \"RooBMixDecay.h\"\n",
    "#include \"RooCategory.h\"\n",
    "#include \"RooBinning.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": "e1c9ac6c",
   "metadata": {},
   "source": [
    "Setup model\n",
    "---------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "da284b56",
   "metadata": {},
   "source": [
    "Build a B decay pdf with mixing"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "2bcaa774",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:11.357903Z",
     "iopub.status.busy": "2026-05-19T20:29:11.357788Z",
     "iopub.status.idle": "2026-05-19T20:29:11.711137Z",
     "shell.execute_reply": "2026-05-19T20:29:11.710323Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar dt(\"dt\", \"dt\", -20, 20);\n",
    "RooRealVar dm(\"dm\", \"dm\", 0.472);\n",
    "RooRealVar tau(\"tau\", \"tau\", 1.547);\n",
    "RooRealVar w(\"w\", \"mistag rate\", 0.1);\n",
    "RooRealVar dw(\"dw\", \"delta mistag rate\", 0.);\n",
    "\n",
    "RooCategory mixState(\"mixState\", \"B0/B0bar mixing state\");\n",
    "mixState.defineType(\"mixed\", -1);\n",
    "mixState.defineType(\"unmixed\", 1);\n",
    "RooCategory tagFlav(\"tagFlav\", \"Flavour of the tagged B0\");\n",
    "tagFlav.defineType(\"B0\", 1);\n",
    "tagFlav.defineType(\"B0bar\", -1);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dae89a6d",
   "metadata": {},
   "source": [
    "Build a gaussian resolution model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "28319b83",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:11.712877Z",
     "iopub.status.busy": "2026-05-19T20:29:11.712761Z",
     "iopub.status.idle": "2026-05-19T20:29:11.920295Z",
     "shell.execute_reply": "2026-05-19T20:29:11.919623Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar bias1(\"bias1\", \"bias1\", 0);\n",
    "RooRealVar sigma1(\"sigma1\", \"sigma1\", 0.1);\n",
    "RooGaussModel gm1(\"gm1\", \"gauss model 1\", dt, bias1, sigma1);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "85d0ae3e",
   "metadata": {},
   "source": [
    "Construct Bdecay (x) gauss"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "bc94858b",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:11.922363Z",
     "iopub.status.busy": "2026-05-19T20:29:11.922246Z",
     "iopub.status.idle": "2026-05-19T20:29:12.129713Z",
     "shell.execute_reply": "2026-05-19T20:29:12.129123Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBMixDecay bmix(\"bmix\", \"decay\", dt, mixState, tagFlav, tau, dm, w, dw, gm1, RooBMixDecay::DoubleSided);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3afd690e",
   "metadata": {},
   "source": [
    "Sample data from model\n",
    "--------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fbaf058a",
   "metadata": {},
   "source": [
    "Sample 2000 events in (dt,mixState,tagFlav) from bmix"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "e8a010e3",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:12.131893Z",
     "iopub.status.busy": "2026-05-19T20:29:12.131775Z",
     "iopub.status.idle": "2026-05-19T20:29:12.339580Z",
     "shell.execute_reply": "2026-05-19T20:29:12.338963Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_64: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{bmix.generate(RooArgSet(dt, mixState, tagFlav), 2000)};\n",
      " ^\n"
     ]
    }
   ],
   "source": [
    "std::unique_ptr<RooDataSet> data{bmix.generate(RooArgSet(dt, mixState, tagFlav), 2000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4e12b20f",
   "metadata": {},
   "source": [
    "Show dt distribution with custom binning\n",
    "-------------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f3e4105c",
   "metadata": {},
   "source": [
    "Make plot of dt distribution of data in range (-15,15) with fine binning for dt>0 and coarse binning for dt<0"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0fafdb65",
   "metadata": {},
   "source": [
    "Create binning object with range (-15,15)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "1a319dea",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:12.341449Z",
     "iopub.status.busy": "2026-05-19T20:29:12.341333Z",
     "iopub.status.idle": "2026-05-19T20:29:12.549038Z",
     "shell.execute_reply": "2026-05-19T20:29:12.548297Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBinning tbins(-15, 15);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1dc26fdf",
   "metadata": {},
   "source": [
    "Add 60 bins with uniform spacing in range (-15,0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "872b09df",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:12.550979Z",
     "iopub.status.busy": "2026-05-19T20:29:12.550863Z",
     "iopub.status.idle": "2026-05-19T20:29:12.758470Z",
     "shell.execute_reply": "2026-05-19T20:29:12.757736Z"
    }
   },
   "outputs": [],
   "source": [
    "tbins.addUniform(60, -15, 0);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "87bb5dc0",
   "metadata": {},
   "source": [
    "Add 15 bins with uniform spacing in range (0,15)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "06480f75",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:12.760353Z",
     "iopub.status.busy": "2026-05-19T20:29:12.760239Z",
     "iopub.status.idle": "2026-05-19T20:29:12.968051Z",
     "shell.execute_reply": "2026-05-19T20:29:12.967331Z"
    }
   },
   "outputs": [],
   "source": [
    "tbins.addUniform(15, 0, 15);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8d8e8a85",
   "metadata": {},
   "source": [
    "Make plot with specified binning"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "2615767d",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:12.969956Z",
     "iopub.status.busy": "2026-05-19T20:29:12.969841Z",
     "iopub.status.idle": "2026-05-19T20:29:13.178346Z",
     "shell.execute_reply": "2026-05-19T20:29:13.177595Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_68:3:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(dtframe, Binning(tbins));\n",
      "^\n",
      "input_line_64:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate(RooArgSet(dt, mixState, tagFlav), 2000)};\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 *dtframe = dt.frame(Range(-15, 15), Title(\"dt distribution with custom binning\"));\n",
    "data->plotOn(dtframe, Binning(tbins));\n",
    "bmix.plotOn(dtframe);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "50744ea8",
   "metadata": {},
   "source": [
    "NB: Note that bin density for each bin is adjusted to that of default frame binning as shown\n",
    "in Y axis label (100 bins --> Events/0.4*Xaxis-dim) so that all bins represent a consistent density distribution"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e007b368",
   "metadata": {},
   "source": [
    "Show mixstate asymmetry with custom binning\n",
    "------------------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6c6c89af",
   "metadata": {},
   "source": [
    "Make plot of dt distribution of data asymmetry in 'mixState' with variable binning"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "81cfe7aa",
   "metadata": {},
   "source": [
    "Create binning object with range (-10,10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "f72fdfe9",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:13.179942Z",
     "iopub.status.busy": "2026-05-19T20:29:13.179828Z",
     "iopub.status.idle": "2026-05-19T20:29:13.387653Z",
     "shell.execute_reply": "2026-05-19T20:29:13.386845Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBinning abins(-10, 10);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "efbdaf10",
   "metadata": {},
   "source": [
    "Add boundaries at 0, (-1,1), (-2,2), (-3,3), (-4,4) and (-6,6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "5cda331c",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:13.389292Z",
     "iopub.status.busy": "2026-05-19T20:29:13.389170Z",
     "iopub.status.idle": "2026-05-19T20:29:13.596777Z",
     "shell.execute_reply": "2026-05-19T20:29:13.596046Z"
    }
   },
   "outputs": [],
   "source": [
    "abins.addBoundary(0);\n",
    "abins.addBoundaryPair(1);\n",
    "abins.addBoundaryPair(2);\n",
    "abins.addBoundaryPair(3);\n",
    "abins.addBoundaryPair(4);\n",
    "abins.addBoundaryPair(6);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0f7f9119",
   "metadata": {},
   "source": [
    "Create plot frame in dt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "3e246d19",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:13.606318Z",
     "iopub.status.busy": "2026-05-19T20:29:13.606195Z",
     "iopub.status.idle": "2026-05-19T20:29:13.811506Z",
     "shell.execute_reply": "2026-05-19T20:29:13.810808Z"
    }
   },
   "outputs": [],
   "source": [
    "RooPlot *aframe = dt.frame(Range(-10, 10), Title(\"mixState asymmetry distribution with custom binning\"));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "26622de7",
   "metadata": {},
   "source": [
    "Plot mixState asymmetry of data with specified custom binning"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "f91be7e2",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:13.813393Z",
     "iopub.status.busy": "2026-05-19T20:29:13.813276Z",
     "iopub.status.idle": "2026-05-19T20:29:14.021771Z",
     "shell.execute_reply": "2026-05-19T20:29:14.021141Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_72:2:2: error: reference to 'data' is ambiguous\n",
      " data->plotOn(aframe, Asymmetry(mixState), Binning(abins));\n",
      " ^\n",
      "input_line_64:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate(RooArgSet(dt, mixState, tagFlav), 2000)};\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": [
    "data->plotOn(aframe, Asymmetry(mixState), Binning(abins));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a068ab5c",
   "metadata": {},
   "source": [
    "Plot corresponding property of pdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "84316005",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:14.023239Z",
     "iopub.status.busy": "2026-05-19T20:29:14.023124Z",
     "iopub.status.idle": "2026-05-19T20:29:14.230480Z",
     "shell.execute_reply": "2026-05-19T20:29:14.230044Z"
    }
   },
   "outputs": [],
   "source": [
    "bmix.plotOn(aframe, Asymmetry(mixState));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0904589b",
   "metadata": {},
   "source": [
    "Adjust vertical range of plot to sensible values for an asymmetry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "3aa8b066",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:14.246410Z",
     "iopub.status.busy": "2026-05-19T20:29:14.246268Z",
     "iopub.status.idle": "2026-05-19T20:29:14.453935Z",
     "shell.execute_reply": "2026-05-19T20:29:14.453232Z"
    }
   },
   "outputs": [],
   "source": [
    "aframe->SetMinimum(-1.1);\n",
    "aframe->SetMaximum(1.1);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8253c857",
   "metadata": {},
   "source": [
    "NB: For asymmetry distributions no density corrects are needed (and are thus not applied)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ac010da5",
   "metadata": {},
   "source": [
    "Draw plots on canvas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "e144741c",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:14.455784Z",
     "iopub.status.busy": "2026-05-19T20:29:14.455666Z",
     "iopub.status.idle": "2026-05-19T20:29:14.663674Z",
     "shell.execute_reply": "2026-05-19T20:29:14.663025Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_82:2:3: error: use of undeclared identifier 'dtframe'\n",
      " (dtframe->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "  ^\n",
      "Error in <HandleInterpreterException>: Error evaluating expression (dtframe->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "Execution of your code was aborted.\n"
     ]
    }
   ],
   "source": [
    "TCanvas *c = new TCanvas(\"rf108_plotbinning\", \"rf108_plotbinning\", 800, 400);\n",
    "c->Divide(2);\n",
    "c->cd(1);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "dtframe->GetYaxis()->SetTitleOffset(1.6);\n",
    "dtframe->Draw();\n",
    "c->cd(2);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "aframe->GetYaxis()->SetTitleOffset(1.6);\n",
    "aframe->Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9ca8635c",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "1676b606",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:29:14.665435Z",
     "iopub.status.busy": "2026-05-19T20:29:14.665317Z",
     "iopub.status.idle": "2026-05-19T20:29:14.892963Z",
     "shell.execute_reply": "2026-05-19T20:29:14.892347Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779222554834\" style=\"width: 800px; height: 400px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779222554834() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(20387,'WkwIsRsAo08AeAHtnG1vHLmRx7+KMchLxmDxqcnmK9leZ3PnXQu2N7azCIKxNLLnLM8oo/GuN8F+98OvyO4ZSXaQPSTABhfAI3fxoatIFv9VLJL9t8Wf9z9drTbLD6vFuHjxcLn5YXn9cvXm+WZ5df1uu1+YxcV3m/VfPq5+/2gxWrO4eLDeX7enp2/+Z3W2J31BsadX+/V204n/Xm/OF6M3i4v5TePfPsfrSwy8D86HYBYXT9ab1cPt5Xa3GKWTz/c/Xa4O5Mv1+f5dIx+vLy97YYSFnApbq9VXF/tvlru3681itPdJebZ+++5W0oPtfr/9cLPYi+3VzYRXF2uEcGZx8frweNIeefGr6/1yD5dSKHODOmkUlR/vlh9Wt+Um7VbD53I3GzQXnZKnV1L9uGd45YPt7ny1e77+a++9o8RvtuerNq6vpP3/uv//yrW+fd3/f7Xfnry5Pl1/Wl2+7zX228+Ri3EoiabfqOAH7bJDlZsJi/G3LeG7G7UYze8OdY7Jic0fblRoL/nDocrNhJnN1JT99lVvjTblLslIWisuJpcl2MH6KNksjl7w+n3rqP4CSOScyMX4W97gUs4yOBu8FTtEFOVy++O3jx72Xj8mXn13pRkM6uuj55dz6tfz08mb6xsvOnlzfeNdJ2+uD9VO3lwfan736YNOB7r4p8Pjpw/LT61B3/00P754t9ovF6PXhr1b96eT66vV2f7Zcr/etlZ8+/HDm9WuPb9Yn73/dHj8qT0+2b7tiU+2bw9pf225p8vz0+V6w/Qxi4uHu+319bvlur9wJk+3HYqO9RrtavRBpb/Znq8v1qvzxXixvLxemcXF73br8083yZ8O5Mmb64fb7e6o/Ffn6/3yDRN+v/vICx6vP63Ob7R7evXpbv1hvV//sLq+g3lP1tdA6gS3nVzudovx+z+ZxfZqz8PPZnHx1afV2fVi3Hy8vDSLi28bPu8uxOY/X11u92/Wm8168xb0fLHeI9bnM7/9+OF0ebna7ycIpTu/XX3a30199Pvnp09OXi/GxW+mR7O4eLT9+OZy9eDjxcU0nM9W++V6Q1/2nnh1vf7r6rvrKf/1TVJzn62Wl4vRwVyzGy3QL9eb8+2PL7ZXTLBj+vUx3YHsUODrFbDd1ePHCQIevlv0ef5wud/f6f6T/b5ZOFr26sFq/+NqtekQfoPSPn282354sb1ajHIfnXp1vtyDkEq8nghM3Ekj5GezeP/N9ofV06vlXz7OmvL+2YoeuZl48fX67bsnNKFbK9Xa5f7s3dSt75+/2/741Q+rzf75frn/eD0r5/uTj/st6jGX/Ga1+fhguWs0ynNyhvrNNS6erZbnTzeXP001Ll6u9++2H/fHmjpp79fL6657U8pxqe9vGfF/mquApf+iq/By9UbRoGv8lxySh5fL6+s+VRjh5qAcJ1yhpAs7uhhN/1UZrbHGVqepPPnRaq6tQfNcjDXOdWxNvSTpQy/Lcx5zMU7E5FALzzkYJ6kKDIPpvyoySnFGsjMyuCpuHJLRf1X8KNaZ/qsSRnGD6b8qcZToTf9VSaMM2fRflWF0Npj+q5JH54rpvypldMGZ/qvOamEp1LfVyU3SjVK8kSJGUq7Oj5KTkQzpqwuNb4pGoq0ujpI9OUZirA6pupDBVjeMEqOR4Ix4W10exUdtYE7VlVGGaJxNSFK9HcU7o6yTr15G8QMSGEm2ejeKt0ZiNrzY+8Yn51aXrnKNT/TVx1FkMFo/xOrTKDYggYrsh1FsMVrBu+rzKC5ovzpbqi+HXo+xBjvKYI02t0gNMjLEKrOPNbjROTHaP5B0VW5CeltDGLWLGTUXaoij9iliSKkhjbymBFNSDUNjQ0eFoYY8KkdRMWsoo4TBaCPE1WhVCFWxGmWUkIzqbI1O+6ETfix5eg7jEKbnOEaZntPohul5GP30mFHs/poyWhNQ4FiTHa1JuT0zdYq0ZzfSKJ1WNTGDxJeWwyRCO7R6hMhDIxI8pDFJOpN856LcdZJKTbBXwg11gL8SNtQBASAk2zoggRIx1WGaw4zvgASaY3MdkAAixzqoADGaJHWYZrIf6jCxF1+Hgp2ikK3Zjk1AnmWM2gye3Zi1G3hmCjd0sTWjldoNEG369lcxTdJEMHVnJkxc7QbqMG21t42tpQOXzxDSEEkCRAMvKcBX8ZojKUGERgQH0VAMpbK1NBgr+uZBSw1CetbnMPBc9Nl5Y6tYO4qICVFBQayohgZngqtiXQOFwvBKFevHPBhRrSE7aDskJuMh4+hcNn4wMVSxaZRozcAELVUskJG0rCRexfR0xnvjHCRdYhK9KlUUWyNzyFURYEyMF+OgnPYqYJRKFbCVZ34uVgFcGU06nvdKbBI3fRbQ1YopxYi1VWQYATjaI7kK6IrepWIERC+j5AJjqeLsKAHMCkYS+A6OtX7LkG7UanYwuVRxfnRe+dgqQGuOmAXDNBewVbxRKHMUTjrVneQmM+h6ZAKclbEIU7JIdQwJkNh+1VmwyZv+q44xEdqjv+psbHgLrFfHoGSvEGCrs8PI3BzEDLw5j86XJpYM1dkyUguwjtUxJGgBlasTGQcxrXJ14hSXe+3qGBT4Ubk6CdRslauTSM1WuTrGA0SXgdrVydCknZjmuXHKtai8E1vMHXDUhK5OhwRpG1+GpLVVsq9Oh0Qb63yozoW5tcqZQemtRWgHih961Q1Te1vlPDeYyn6aw2CbDdVP0xhlwm4xk7sVb3SbzJpvXfXM5ylf6SOfQOk2qykfpXomdi+uZJ49CCXBtuZuRKnBHhwjW4NtPgHVodxcE+rgp0AFzAkWpQZLbzT0on1BJ7Niv7YnMJmhscbR14DqWDe3L4htNPnW1YCzhIhYUKUdvHrzgnilbFD5ReWgKK0RNXJUVmqCeluDKNS3dglQTxXSMTRUsDU47ExrkutD0nzGGnCHepMwaYFJ25ukTUY/epNaPkjitMktH++jN4kucHlqUitepjZpaW8hYaCFvSjZRq0GPxk9rerV6rXxrsGr4Z0rquWd69EdOmg1opJdR3hLLDIPdaMPQ49EsRwGv+WH1p7WxhpRycZHRY4lIXIbtRrLMFH6Lm17a2CNRZsOZV1NliFQq1iTbV46dZKlzTpoNVlarINWk23eOh4DWthmmIqYgK82wzqdtcl92GrCppDfhq0mhTDVwMYSLWwaqE1KmBadYe114vV1NBEBsSyqhL1w64+pWXK8fEjdawf1qNlmR28cVgVtbN3hVBEYtZpQyal91AOxevvQkgRktRnW6W74+/AnxazWPu1q9dJ1hnW6mbemIMmVqXnwxktvrVOqzc9e0h8WUcm3TsEjSXjnzDF9bt3RnlU1eBwmIGxi16EDYe/2OqB1XcQodUDpOknlOLdXMw/DT/OG7t6oCHXoDg6ctHBrXcvMEwYqOmQwsPdrlJptQwIdlJo7CCpIZDuPOTldnIY9uWtjl7fmSRsVaTK62HVPecydjY/ZELmx6EpIYQpiSjusUFKRsL/yAITkdDhqzc3Y0QmQrasZO9rUXgvnmVIeMxbh5KoCwh6iW4fWOzVjQidc57WAJNrWwDiDkZ2kcveEu0hdG3sv1Ny1ceLU5ioUIrFm7Ixwtjt2NaFQxt42pkVmydi7iGmRfeukpjc1HwASib3aC2RqZTtqd0Od/bQ+aIUPmK2MdIkyVQ3Ha5QcFLF5j743KGK3yVdz0JVS15MwATbtDAdxkC4cxGmvmcWhYWEWpxWexOFFcZJGid5DzUTXzHKxj4s2JM5qpO+NSNTRObdVY+/QmmOTqHd4zRGJNJdQSFSJDiQSdXzMuobs9iCnZtDasNXc1pFt0GrWhWQHztwWkk3amlMz700Vak5HRhb+upw8vLX1z9RqXVHOVduScsobeihkEmjo1m8SaWiKPXXK0BS7IXYemmJPDR2aMzT3wsDKZlph56Ep9tQRLDRbJhqVm2I385dzG7apaG6j1oxGzjpmzUrkrD5QF4cVp+pWG1zWnB0ztY9ym/oTjzb1Zx469TuL0m2PVutIPfVAB+rOvzsHaqZyB+lONFCcRrNj9DyeHaOn18wGSBW6Y7S+qHSI7sQBoem3cgTRaHABpLsttq6WDtONS+k43VtSJpzuQ106TrfOLN1nmKjmMfSRLh2oe2aPrqlOIRTLDTq+/WoBLNHj9qvF6Qg2G1pcGKNnWS2eqqrmgaXxUItLY7S6Ag8s+ocxFBbqEnhpHkM2LDxDqsWVMSQTrZGQa/F2DMCnkSi1eBlDMCpMqMW7MXij8DTU4v0YnGH5nYgkhDFYQ8QtuVp8HH0xsRALq8Wn0WeTWoiw+GH0g2GJPMAhjz4aFsADHAphk5SMDKmWYEfvTRqMDKWWIKN3rOZZG5fgRq9LPsmxluBHVwxr4TzUEsLosgYJCmGNSNhqiEQLayGSlswwGClwGIgQDMVIgQOLVZPRZKkllNE5k5tfVzSOZnTpN9QCOBaTIzGpWoDGbHIyTnwtkeiBydk4ibXEMEoyhehVriXqgpBAmJNawMVoiusUU17DE4534t8RcmmULptLalTCrzJl6BQqYwrrfEI1xAyIRTgptSQPJQQulFRhVGuVjC3X99zUyBa4KGloJI3MtaRMs8TSSkgCsobQhJIDoU9DqNTJUMsgShIOU5KQjyEk1EjtIRHXSaLEhugFIfACMmYjGmuAJLrQ4oVKDo2EL7ks5o0Q5aC3QUZIxIi1ZKskkQIlpZGIQS6xaw16NNI3EjFCLQTkyKX5kARkjTiaD5kaiRiQQyMRw9cCQFIYMSCLkh4xfC2lSUXsSVwtLKIy0cdOtr4i3iRSCyhJdBoxINsIesSAbCPoEYMYnCqUeMSAVI3S6KqSqlIS0O9SC0v6YIQoliWWRtAOGoXP0ASjAIyJJqpuJBCO0Hw/ChHoZJwdKE8Ek4j0RMdG54kmLGMkFCLilB+UjoijNFYF2DHORvI1nCoReaDFjqWbYiVFSd3GINeNBNNbrESseCWRRXODkojiyY0jIR0kUTJBJgQh4kgAj3jiRGYlEUNzy0g0j14hxOis1tXoIKTACDgr9KlzCMn2R6HL1NXVcKsUugxf13WaNuHtEoMEnRBbg3idRlB1d9s2R0EWQkZIDn6pME39FSNhzwKssG1jJMPf47F0Gv4swtha0l0JaCKMLaaYkccTi+008ngWhy0uSnzSet0R0Dhpi+QOAK8Qns30hs/AMhs3ILhY4B3+0MgDwMO/4b1YIB7+0PAH5OGPBYB/8BiLFsekP0LAlDQaeULE0Ah4PsA/JMyQZPoc/mHASDVa+WdMWNutoj9CwcA1mv6IFvOn228JeaJgHDUsmuAfHYa00fDvPrHucTEeMYyR+hh0+OMXD0YyCwf4xzRiI5VWfsOYkId85ZdHdEdp2h/LyBYXcV8NlSdLBLHR8EsyErMmP8Av6YZho+GX/EjfsEMX4JfCyNai0vBLcYSX0vBLSQPdStO+NIwatcePUH55ZPMPWsPxqYyF3T9o+A2WXaVG07+YBYvA4B0MMQwscDQBjuoz00WKiCSEETPSEuCJeWDS5KLbhGIxEERZSXDKFWCBLQnKlm0B2JKgbNm3Y45a3WwTi6FwbC2SwMhiKlBdTYAtxkJ3HCmBbmEuUEZ2PonzWwwG2kgJgS0mA/U7JBCTZ9eSEvQ/ZgMF0wQEw3CELge7BRbTgYppCQTDeOCdaQIdNMXgNAFJNQjX+0OQVOMhvT8U6dWI0KfsWCC6mpHep4rlGBLVPEogKaYEWNIEJMWYNJxhq5EEHII+tiToJpAiQTYCzOo2EFNbE9h1mXaCNEF3aNh3mN7BVgo2RTcisgFP2Q1SbYZmb0W3HmhaNgCqbgjREXlgo1R0Q0g30AcDwLIjJAV9GwwAy56QFJo1GACWbSFhHwdad4kI4yDQRLdIldIIrNsQNAkXD5qVC/Ik3bFha4hdS8oDwOwNsSWsNPKxOaTKnNlkFHaH2GShPQAuG0ROdTubjHzEUkB8tnPhT3SvqzqAq1tEIHyxBsBll0j3nKHhT3QPxC9iBvqLeIpOBDEALrtEDoSHVv7sMKLTzgC4QnwPhIemfwiq6DTxBsAVDAxbW8WbtpHmRqezJhgAV9TAoPDBALiCgdFJhNMPHUfn0X+WBNBpZA9TafhjYHRKJQPgCgYmwC+xMSvi29kFxhTAlWBHpxNsMACusCUPgkPDnz15nW/ZAMDCOkKnWzYAsLCSCLo6MBH+rCV09hUDAAurCXVQdIkpwnpC4wzWRPqDFYXSzgDIwppCaW8C/FlVKB0NgCyRNSquRTIAskSipOqAGQBZItFbPJFiAvJE9mxxHEW3ViUiH7Q3ALSwZa90MAC0ROQjPxkAWiLyQQ9sZotE5IMuxiNPQj5WK9aA15KQD9oZ4FoS8kEH45EnIR90NIC1JOSDHoynP1LrP3SITUpJyGeNwzVGnoR80GJAbknIB+2NU3mQD5ptc7ZakQ86GWBbCMkozWoIGvnIxxOHZnyLQWcAcRna+HLWRZBniD2fOQidOo1HAs16kfp4XdDoHzSbuxwkevDxgrNkv9GDeovHl9vl3ruFWVzq6bEYzeKHxfh9ccEUx9IvmeLAlmyKK6Z4awpbvN6Z4r0pPpjioyk+meLBqGyKL6YEa0oQU4IzJXhTQjAlsHjkjAjYlk0JxZRoTYliSnSmRG9KDKawBx6TKRFMzKbEYkqypiQxJTlTkjclcdAkmpKSKQkszaakYspgTRnElMGZMnhThmDKwFI0mTIMpgzgcDElW1OymJKdKdmbkoMpmUVqMgUsVfwuphRrCnvQxZlSvCkccCnRlMISdjClgPu6r04cBiS2Ggmx7EBoZNSyFwMEW9wGC/haDuxYYNeCtRaAtfjVFmi14KnVVTFIaoFPi5dsAU6Ln2CBSIsnbAFHCyJa3AMLFloA0IJ6FmfWgncWkLP4rhZ4s3pwCSCzeAQWCLPglsX9tCCWReUs3qZF2SyoZFEz2w7vUAOtsgRmLJhjIzU0KodTabHxlvCLxae0WHiLN2kx7ZZoi8WmWyItlhWJxaJbgisWU25xJS023LJAsBhvSwjF4kNaTDeHjnGw+EMNDhBYPEiLS20x1RYH0mKjLa6jxfe0OIwWL9CqaVenAbNscZYsBtniN1pMsS3UwG20WGKLw2gZc6IWOA384ZADY95ONTHmBCKw9/xhxc+Yq6NIlEHUQ2znIRhz9Q+FMW8HAxhzPS9FaACryh89SAEPxpx1P+aRP/BgzNUhbCeyGHN1B4UxZ7GO7eIPNRhzFuHYIwIZ8GDM1fljQY1h4Q81GHM9sKV+HstjzAJ/qMGYazBQGHM9H6cnXoQx1wNuLGMBav5QgzHXE3a6Cm0nwxhzdd/0hJz6bbooVIetHSFhzNVd09WZ+mm67pIh/+nnn382/6rTnJxq/+JpznZR4e9cLOknNjm+vdtQcDri3FLuzS+4dUMC8ugWxO0LEIerIcKFjLt3Q75Z7t6vdkd3TVrC0St7wnx/4sXq0/5k85bz1xxNhmyZ9r6lDzT/cv12syAq0uij95P9eMsp96QnjJef1nfPjJ/s9yekc3D7fP3D+nq93Vwvxqinp8k5euGT5ZvVdP8Ffko3DgEOSj+9uLhe6cUUcLYnzmJ7lXt99v7JavOWazX2vuWQs47BVFXbwrnz29X2l9NB87nIxJ6Dya//bVqow/l/aOEf/21aOA/QLxzDB8vd0bWiB8vdpBR0meoi0/Ty9HmbE492yx/bdYxGP73aH65+NKLf/mhEvwDy9Gr/qJ23b1fLOATPJNJp9PRq3yGBRjy92j/Wq1i96ON1vxdw5wQ9BTTxfL3nPtlEv9huL/UAPQnt6srD7Wa//bi77rcWTvZdnFuIebLfM4kVpP4OFrhfCAbMFZrfLh3RSChuRNxX6qvN+Ve73bbf5GJmK6nFYfX44+aswwKZkEcoBtmHkFyuwvTCtB+yF2beQx6N95PV29Xm/Pi+DdK11COE5UWHxIn3dAWPV0ywMBek/7oqmsXF19yFWF3fwvCe+vxqecZdAOU934s7asN8Ka6nIeNc7qY0c9EpeSp6i7WWu93uo8TDZaOv19co5LE8JPG+Lk6yNHsqNzFuvTMVnVJ7wVvSUOqb9Wb94eOHP65228NVDzJuXFVUiG+3Xk53q4vV7ndPDqVb+lHHtYTjZiLpceqhnS310eri68UYLYM2p7xcjMPNlFeLdt1nLvK6J5wuj5XvdHlDt2A+Jx04a9LdK5qny/MbTafvTpfnd698ni7PP3Pr83R5jrK/OnRPT3l9IwXb2O8qwXB99r7fVDpdXrW7la86aMwJrxdsTy4unp/tVqvN4+WZog/iAWtH3Q/JPDhSW5KOx2OqdTR/KAJ5mD2kTArU1KqV2X0ARRfpftCrYyQqZhQuMkIpxLA4beTL6b4TeV9DiLZjv1tfPVqdrT8sL6/nS0UKyd3bcbPLcNQ6LXCreZp23D70SBOPGjjTs6lqTgoOSLsKPNeaIVDrtOaUuTTtmd5Pcxr4bH+/2ax2z2geJZlq+trrxfg9F4nu3eOPuHvCov9euJcnOt0L9yRBOUuB0IscKE06Tnf2ZvE5L9uel+29YO+Fe0FLKu9/0Z8/AS6r5flqh7XWO1PabTP1eL1/PClN7EqjN7EYxTlHdUl77Wx5qZUZ/f/arjckTo7Aw+XVMfli/WF2J4eci+SgiPH7D8u3K140A/zD5eb8cvXy3fr6/Wr3bLl52y8vt/QH2089rY1eS1VJjq5v/mG9vVxvptR+cbEVfbjenV3eRvuexb1ThD4ygK9wub/6dPXq2O2ZEl8fJ77+XMkp8UZJCn6z/PRo/VYv16OET3f7d9uHyw+r3bKjz93l2uny/J92/44R++KK7XR5/qX1GispeugIsSAn8Om5R/13yzLeXYXdgms64zNg/Z8r+oDXP3ZFf7hfYimlZOdsJly24A7/52/s++zu22yttT7GTNFbt/dTvu8032ZPcEMLTO/yKdyPidqSkivDUS4X7u/m3rzk/xlBv3jn/66gt+7/3xX0+GMAd0WZcj8v6M3PBHCR39nsvfcuh+ySmph+0f/2dwTUa7cuyeCiKy7nHIfPfDjgviVYGWwWXwidHj4c8IXvCAzBWw/rVHJgHh6u/9+3lgEvpQxDCtlFHaf5ywKfzZ6+NWDvh8IIWmud81F0EPu3B/iWxe289lEBez80rbA2BL6NgGlpHzW4X+7k3PhWwWelufH1gi+V+DLr+fsGn+P+L/zeAYD3n+8dzN+T+cXfO7i1vP7H63/pCwl/luMA4p0PKLTsf/QrCf+8u/XdA7lxtf5Pd+Kx/y8M/H/s+/RdnsMC95d9gifYeGzhf7Xm/a6cv1LrLvdtcgHzmWMIGRX9Fdn2eMO247j/w5b9rvX+11r2u6L+m9p13LH/2PVfk10nFNOiFZ/79lHL/rXYdeIsz8926ys+nadhFj7/c9k+/7PZrM72hwijfleqxbn0Azw//y/h0w0t').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779222554834', 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_1779222554834();\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%jsroot on\n",
    "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
}