{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "7fe30cfa",
   "metadata": {},
   "source": [
    "# rf701_efficiencyfit\n",
    "Special pdf's: unbinned maximum likelihood fit of an efficiency eff(x) function\n",
    "\n",
    "to a dataset D(x,cut), where cut is a category encoding a selection, of which the\n",
    "efficiency as function of x should be described by eff(x)\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:34 PM.</small></i>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "26dde609",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:16.059496Z",
     "iopub.status.busy": "2026-05-19T20:34:16.059386Z",
     "iopub.status.idle": "2026-05-19T20:34:16.074698Z",
     "shell.execute_reply": "2026-05-19T20:34:16.074171Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "#include \"RooRealVar.h\"\n",
    "#include \"RooDataSet.h\"\n",
    "#include \"RooGaussian.h\"\n",
    "#include \"RooConstVar.h\"\n",
    "#include \"RooFormulaVar.h\"\n",
    "#include \"RooProdPdf.h\"\n",
    "#include \"RooEfficiency.h\"\n",
    "#include \"RooPolynomial.h\"\n",
    "#include \"RooCategory.h\"\n",
    "#include \"TCanvas.h\"\n",
    "#include \"TAxis.h\"\n",
    "#include \"RooPlot.h\"\n",
    "using namespace RooFit;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c26982a3",
   "metadata": {},
   "source": [
    "Construct efficiency function e(x)\n",
    "-------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b8685809",
   "metadata": {},
   "source": [
    "Declare variables x,mean,sigma with associated name, title, initial value and allowed range"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "b19a3e6f",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:16.075953Z",
     "iopub.status.busy": "2026-05-19T20:34:16.075843Z",
     "iopub.status.idle": "2026-05-19T20:34:16.391697Z",
     "shell.execute_reply": "2026-05-19T20:34:16.391067Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar x(\"x\", \"x\", -10, 10);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "73acfd12",
   "metadata": {},
   "source": [
    "Efficiency function eff(x;a,b)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "dc26f4ec",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:16.393304Z",
     "iopub.status.busy": "2026-05-19T20:34:16.393186Z",
     "iopub.status.idle": "2026-05-19T20:34:16.601963Z",
     "shell.execute_reply": "2026-05-19T20:34:16.601309Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar a(\"a\", \"a\", 0.4, 0, 1);\n",
    "RooRealVar b(\"b\", \"b\", 5);\n",
    "RooRealVar c(\"c\", \"c\", -1, -10, 10);\n",
    "RooFormulaVar effFunc(\"effFunc\", \"(1-a)+a*cos((x-c)/b)\", RooArgList(a, b, c, x));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2fe6c783",
   "metadata": {},
   "source": [
    "Construct conditional  efficiency pdf E(cut|x)\n",
    "------------------------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b530fd30",
   "metadata": {},
   "source": [
    "Acceptance state cut (1 or 0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "c206116a",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:16.603580Z",
     "iopub.status.busy": "2026-05-19T20:34:16.603463Z",
     "iopub.status.idle": "2026-05-19T20:34:16.809253Z",
     "shell.execute_reply": "2026-05-19T20:34:16.808648Z"
    }
   },
   "outputs": [],
   "source": [
    "RooCategory cut(\"cut\", \"cutr\", { {\"accept\", 1}, {\"reject\", 0} });"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fea81594",
   "metadata": {},
   "source": [
    "Construct efficiency pdf eff(cut|x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "55169641",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:16.811166Z",
     "iopub.status.busy": "2026-05-19T20:34:16.811049Z",
     "iopub.status.idle": "2026-05-19T20:34:17.017020Z",
     "shell.execute_reply": "2026-05-19T20:34:17.016244Z"
    }
   },
   "outputs": [],
   "source": [
    "RooEfficiency effPdf(\"effPdf\", \"effPdf\", effFunc, cut, \"accept\");"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0fd1314b",
   "metadata": {},
   "source": [
    "Generate data (x, cut) from a toy model\n",
    "-----------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "49d6db43",
   "metadata": {},
   "source": [
    "Construct global shape pdf shape(x) and product model(x,cut) = eff(cut|x)*shape(x)\n",
    "(These are _only_ needed to generate some toy MC here to be used later)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "0399fc76",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:17.018559Z",
     "iopub.status.busy": "2026-05-19T20:34:17.018441Z",
     "iopub.status.idle": "2026-05-19T20:34:17.227325Z",
     "shell.execute_reply": "2026-05-19T20:34:17.226697Z"
    }
   },
   "outputs": [],
   "source": [
    "RooPolynomial shapePdf(\"shapePdf\", \"shapePdf\", x, RooConst(-0.095));\n",
    "RooProdPdf model(\"model\", \"model\", shapePdf, Conditional(effPdf, cut));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "415c99ed",
   "metadata": {},
   "source": [
    "Generate some toy data from model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "6257e57e",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:17.228944Z",
     "iopub.status.busy": "2026-05-19T20:34:17.228826Z",
     "iopub.status.idle": "2026-05-19T20:34:17.437950Z",
     "shell.execute_reply": "2026-05-19T20:34:17.437361Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_60: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, cut}, 10000)};\n",
      " ^\n"
     ]
    }
   ],
   "source": [
    "std::unique_ptr<RooDataSet> data{model.generate({x, cut}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "00d8d71d",
   "metadata": {},
   "source": [
    "Fit conditional efficiency pdf to data\n",
    "--------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "be946cda",
   "metadata": {},
   "source": [
    "Fit conditional efficiency pdf to data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "e45d54f2",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:17.439237Z",
     "iopub.status.busy": "2026-05-19T20:34:17.439124Z",
     "iopub.status.idle": "2026-05-19T20:34:17.648722Z",
     "shell.execute_reply": "2026-05-19T20:34:17.648166Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_61:2:16: error: reference to 'data' is ambiguous\n",
      " effPdf.fitTo(*data, ConditionalObservables(x), PrintLevel(-1));\n",
      "               ^\n",
      "input_line_60:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, cut}, 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": [
    "effPdf.fitTo(*data, ConditionalObservables(x), PrintLevel(-1));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d9388536",
   "metadata": {},
   "source": [
    "Plot fitted, data efficiency\n",
    "--------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5d94d2a4",
   "metadata": {},
   "source": [
    "Plot distribution of all events and accepted fraction of events on frame"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "3c220ac8",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:17.650082Z",
     "iopub.status.busy": "2026-05-19T20:34:17.649966Z",
     "iopub.status.idle": "2026-05-19T20:34:17.860215Z",
     "shell.execute_reply": "2026-05-19T20:34:17.859659Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_62:3:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame1);\n",
      "^\n",
      "input_line_60:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, cut}, 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",
      "input_line_62:4:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame1, Cut(\"cut==cut::accept\"), MarkerColor(kRed), LineColor(kRed));\n",
      "^\n",
      "input_line_60:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, cut}, 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 *frame1 = x.frame(Bins(20), Title(\"Data (all, accepted)\"));\n",
    "data->plotOn(frame1);\n",
    "data->plotOn(frame1, Cut(\"cut==cut::accept\"), MarkerColor(kRed), LineColor(kRed));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1a5fb02a",
   "metadata": {},
   "source": [
    "Plot accept/reject efficiency on data overlay fitted efficiency curve"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "b3e3c4dc",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:17.861807Z",
     "iopub.status.busy": "2026-05-19T20:34:17.861693Z",
     "iopub.status.idle": "2026-05-19T20:34:18.071208Z",
     "shell.execute_reply": "2026-05-19T20:34:18.070689Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_63:3:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame2, Efficiency(cut)); // needs ROOT version >= 5.21\n",
      "^\n",
      "input_line_60:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{model.generate({x, cut}, 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 *frame2 = x.frame(Bins(20), Title(\"Fitted efficiency\"));\n",
    "data->plotOn(frame2, Efficiency(cut)); // needs ROOT version >= 5.21\n",
    "effFunc.plotOn(frame2, LineColor(kRed));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5a305959",
   "metadata": {},
   "source": [
    "Draw all frames on a canvas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "6abe53e4",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:18.072782Z",
     "iopub.status.busy": "2026-05-19T20:34:18.072657Z",
     "iopub.status.idle": "2026-05-19T20:34:18.281850Z",
     "shell.execute_reply": "2026-05-19T20:34:18.281129Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_65:2:3: error: use of undeclared identifier 'frame1'\n",
      " (frame1->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "  ^\n",
      "Error in <HandleInterpreterException>: Error evaluating expression (frame1->GetYaxis()->SetTitleOffset(1.6000000000000001))\n",
      "Execution of your code was aborted.\n"
     ]
    }
   ],
   "source": [
    "TCanvas *ca = new TCanvas(\"rf701_efficiency\", \"rf701_efficiency\", 800, 400);\n",
    "ca->Divide(2);\n",
    "ca->cd(1);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame1->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame1->Draw();\n",
    "ca->cd(2);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame2->GetYaxis()->SetTitleOffset(1.4);\n",
    "frame2->Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "298e9273",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "e5d4c361",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:18.283135Z",
     "iopub.status.busy": "2026-05-19T20:34:18.283021Z",
     "iopub.status.idle": "2026-05-19T20:34:18.512658Z",
     "shell.execute_reply": "2026-05-19T20:34:18.511899Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779222858453\" style=\"width: 800px; height: 400px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779222858453() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(20381,'WkwItBsAnU8AeAHtnG1vHLmRx7+KMchLxmDxqcnmK9leZ3PnXQtrb2xnEQRjaWTPWZ5RRuNdO8F+98OvyO4ZSXaQPSTABhfAI3fxoatIFv9VLJL9t8Wf95+uVpvl+9ViXDx/uNz8uLx+sXr9bLO8un673S/M4uL7zfovH1a/f7QYrVlcPFjvr9vT09f/szrbk76g2NOr/Xq76cR/rzfni9GbxcX8pvFvn+P1JQbeB+dDMIuLJ+vN6uH2crtbjNLJZ/tPl6sD+WJ9vn/byMfry8teGGEhp8LWavXVxf6b5e7NerMY7X1Svlu/eXsr6cF2v9++v1ns+fbqZsLLizVCOLO4eHV4PGmPvPjl9X65h0splLlBnTSKyo93y/er23KTdqvhc7mbDZqLTsnTK6l+3DO88sF2d77aPVv/tffeUeI32/NVG9eX0v5/1f9/6Vrfvur/v9xvT15fn64/ri7f9Rr77efIxTiURNNvVPCDdtmhys2ExfjblvD9jVqM5veHOsfkxOYPNyq0l/zhUOVmwsxmasp++7K3Rptyl2QkrRUXk8sS7GB9lGwWRy949a51VH8BJHJO5GL8LW9wKWcZnA3eih0iinK5/enbRw97rx8TL7+/0gwG9dXR84s59ev56eT19Y0Xnby+vvGuk9fXh2onr68PNb//+F6nA1386fD48f3yY2vQ95/mx+dvV/vlYvTasLfr/nRyfbU623+33K+3rRXffnj/erVrz8/XZ+8+Hh4/tccn2zc98cn2zSHtry33dHl+ulxvmD5mcfFwt72+frtc9xfO5Om2Q9GxXqNdjT6o9Dfb8/XFenW+GC+Wl9crs7j43W59/vEm+elAnry+frjd7o7Kf3W+3i9fM+H3uw+84PH64+r8RrunV5/u1u/X+/WPq+s7mPdkfQ2kTnDbyeVutxh/+JNZbK/2PPxsFhdffVydXS/GzYfLS7O4+Lbh8+5isPLn1cXF+my92px9Ajyfr/dI9dm8bz+8P11ervb7CUDpzG9XH/d3Ux/9/tnpk5NXi3Hxm+nRLC4ebT+8vlw9+HBxMQ3md6v9cr2hJ3s/vLxe/3X1/fWU/+omqbnfrZaXi9HBXLMbLdAv1pvz7U/Pt1dMr2P61THdYexQ4OsVoN2V46cJAB6+XfRZ/nC539/p/JP9vtk3WvbywWr/02q16QB+g9Iufbzbvn++vVqMch+Nenm+3IOPSryaCAzcSSPkZ7N49832x9XTq+VfPsx68u67FT1yM/Hi6/Wbt09oQrdVqrPL/dnbqVvfPXu7/emrH1eb/bP9cv/helbNdycf9luUYy75zWrz4cFy12hU5+QM5ZtrXHy3Wp4/3Vx+mmpcvFjv324/7I/1dNLdr5fXXfOmlONSP9wy4f80RwE7/0VH4cXqtWLBevPmS94CmvHwcnl93ScK5Zp7cpxwhZIu7OhiNP1XZbTGGludpvLkR6u5tgbNczHWONexNfWSpA+9LM95zMU4EZNDLTznYJykKjAMpv+qyCjFGcnOyOCquHFIRv9V8aNYZ/qvShjFDab/qsRRojf9VyWNMmTTf1WG0dlg+q9KHp0rpv+qlNEFZ/qvOquFpVDfVic3STdK8UaKGEm5Oj9KTkYypK8uNL4pGom2ujhK9uQYibE6pOpCBlvdMEqMRoIz4m11eRQftYE5VVdGGaJxNiFJ9XYU74yyTr56GcUPSGAk2erdKN4aidnwYu8bn5xbXbrKNT7RVx9HkcFo/RCrT6PYgAQqsh9GscVoBe+qz6O4oP3qbKm+HHo9xhrsKIM12twiNcjIEKvMPtbgRufEaP9A0lW5CeltDWHULmbUXKghjtqniCGlhjTymhJMSTUMjQ0dFYYa8qgcRcWsoYwSBqONEFejVSFUxWqUUUIyqrM1Ou2HTvix5Ok5jEOYnuMYZXpOoxum52H002NGsftrymhNQIFjTXa0JuX2zNQp0p7dSKN0WtXEDBJfWg6TCO3Q6hEiD41I8JDGJOlM8p2LctdJKjXBXgk31AH+SthQBwSAkGzrgARKxFSHaQ4zvgMSaI7NdUACiBzroALEaJLUYZrJfqjDxF58HQp2ikK2Zjs2AXmWMWozeHZj1m7gmSnc0MXWjFZqN0C06dtfxTRJE8HUnZkwcbUbqMO01d42tpYOXD5DSEMkCRANvKQAX8VrjqQEERoRHERDMZTK1tJgrOibBy01COlZn8PAc9Fn542tYu0oIiZEBQWxohoanAmuinUNFArDK1WsH/NgRLWG7KDtkJiMh4yjc9n4wcRQxaZRojUDE7RUsUBG0rKSeBXT0xnvjXOQdIlJ9KpUUWyNzCFXRYAxMV6Mg3Laq4BRKlXAVp75uVgFcGU06XjeK7FJ3PRZQFcrphQj1laRYQTgaI/kKqArepeKERC9jJILjKWKs6MEMCsYSeA7ONb6LUO6UavZweRSxfnReeVjqwCtOWIWDNNcwFbxRqHMUTjpVHeSm8yg65EJcFbGIkzJItUxJEBi+1VnwSZv+q86xkRoj/6qs7HhLbBeHYOSvUKArc4OI3NzEDPw5jw6X5pYMlRny0gtwDpWx5CgBVSuTmQcxLTK1YlTXO61q2NQ4Efl6iRQs1WuTiI1W+XqGA8QXQZqVydDk3ZimufGKdei8k5sMXfAURO6Oh0SpG18GZLWVsm+Oh0SbazzoToX5tYqZwaltxahHSh+6FU3TO1tlfPcYCr7aQ6DbTZUP01jlAm7xUzuVrzRbTJrvnXVM5+nfKWPfAKl26ymfJTqmdi9uJJ59iCUBNuauxGlBntwjGwNtvkEVIdyc02og58CFTAnWJQaLL3R0Iv2BZ3Miv3ansBkhsYaR18DqmPd3L4gttHkW1cDzhIiYkGVdvDqzQvilbJB5ReVg6K0RtTIUVmpCeptDaJQ39olQD1VSMfQUMHW4LAzrUmuD0nzGWvAHepNwqQFJm1vkjYZ/ehNavkgidMmt3y8j94kusDlqUmteJnapKW9hYSBFvaiZBu1Gvxk9LSqV6vXxrsGr4Z3rqiWd65Hd+ig1YhKdh3hLbHIPNSNPgw9EsVyGPyWH1p7WhtrRCUbHxU5loTIbdRqLMNE6bu07a2BNRZtOpR1NVmGQK1iTbZ56dRJljbroNVkabEOWk22eet4DGhhm2EqYgK+2gzrdNYm92GrCZtCfhu2mhTCVAMbS7SwaaA2KWFadIa114nX19FEBMSyqBL2wq0/pmbJ8fIhda8d1KNmmx29cVgVtLF1h1NFYNRqQiWn9lEPxOrtQ0sSkNVmWKe74e/DnxSzWvu0q9VL1xnW6WbemoIkV6bmwRsvvbVOqTY/e0l/WEQl3zoFjyThnTPH9Ll1R3tW1eBxmICwiV2HDoS92+uA1nURo9QBpesklePcXs08DD/NG7p7oyLUoTs4cNLCrXUtM08YqOiQwcDer1Fqtg0JdFBq7iCoIJHtPObkdHEa9uSujV3emidtVKTJ6GLXPeUxdzY+ZkPkxqIrIYUpiCntsEJJRcL+ygMQktPhqDU3Y0cnQLauZuxoU3stnGdKecxYhJOrCgh7iG4dWu/UjAmdcJ3XApJoWwPjDEZ2ksrdE+4idW3svVBz18aJU5urUIjEmrEzwtnu2NWEQhl725gWmSVj7yKmRfatk5re1HwASCT2ai+QqZXtqN0NdfbT+qAVPmC2MtIlylQ1HK9RclDE5j363qCI3SZfzUFXSl1PwgTYtDMcxEG6cBCnvWYWh4aFWZxWeBKHF8VJGiV6DzUTXTPLxT4u2pA4q5G+NyJRR+fcVo29Q2uOTaLe4TVHJNJcQiFRJTqQSNTxMesastuDnJpBa8NWc1tHtkGrWReSHThzW0g2aWtOzbw3Vag5HRlZ+Oty8vDW1j9Tq3VFOVdtS8opb+ihkEmgoVu/SaShKfbUKUNT7IbYeWiKPTV0aM7Q3AsDK5tphZ2HpthTR7DQbJloVG6K3cxfzm3YpqK5jVozGjnrmDUrkbP6QF0cVpyqW21wWXN2zNQ+ym3qTzza1J956NTvLEq3PVqtI/XUAx2oO//uHKiZyh2kO9FAcRrNjtHzeHaMnl4zGyBV6I7R+qLSIboTB4Sm38oRRKPBBZDutti6WjpMNy6l43RvSZlwug916TjdOrN0n2GimsfQR7p0oO6ZPbqmOoVQLDfo+ParBbBEj9uvFqcj2GxocWGMnmW1eKqqmgeWxkMtLo3R6go8sOgfxlBYqEvgpXkM2bDwDKkWV8aQTLRGQq7F2zEAn0ai1OJlDMGoMKEW78bgjcLTUIv3Y3CG5XcikhDGYA0Rt+Rq8XH0xcRCLKwWn0afTWohwuKH0Q+GJfIAhzz6aFgAD3AohE1SMjKkWoIdvTdpMDKUWoKM3rGaZ21cghu9Lvkkx1qCH10xrIXzUEsIo8saJCiENSJhqyESLayFSFoyw2CkwGEgQjAUIwUOLFZNRpOlllBG50xufl3ROJrRpd9QC+BYTI7EpGoBGrPJyTjxtUSiByZn4yTWEsMoyRSiV7mWqAtCAmFOagEXoymuU0x5DU843ol/R8ilUbpsLqlRCb/KlKFTqIwprPMJ1RAzIBbhpNSSPJQQuFBShVGtVTK2XN9zUyNb4KKkoZE0MteSMs0SSyshCcgaQhNKDoQ+DaFSJ0MtgyhJOExJQj6GkFAjtYdEXCeJEhuiF4TAC8iYjWisAZLoQosXKjk0Er7kspg3QpSD3gYZIREj1pKtkkQKlJRGIga5xK416NFI30jECLUQkCOX5kMSkDXiaD5kaiRiQA6NRAxfCwBJYcSALEp6xPC1lCYVsSdxtbCIykQfO9n6iniTSC2gJNFpxIBsI+gRA7KNoEcMYnCqUOIRA1I1SqOrSqpKSUC/Sy0s6YMRoliWWBpBO2gUPkMTjAIwJpqoupFAOELz/ShEoJNxdqA8EUwi0hMdG50nmrCMkVCIiFN+UDoijtJYFWDHOBvJ13CqROSBFjuWboqVFCV1G4NcNxJMb7ESseKVRBbNDUoiiic3joR0kETJBJkQhIgjATziiROZlUQMzS0j0Tx6hRCjs1pXo4OQAiPgrNCnziEk2x+FLlNXV8OtUugyfF3XadqEt0sMEnRCbA3idRpB1d1t2xwFWQgZITn4pcI09VeMhD0LsMK2jZEMf4/H0mn4swhja0l3JaCJMLaYYkYeTyy208jjWRy2uCjxSet1R0DjpC2SOwC8Qng20xs+A8ts3IDgYoF3+EMjDwAP/4b3YoF4+EPDH5CHPxYA/sFjLFock/4IAVPSaOQJEUMj4PkA/5AwQ5Lpc/iHASPVaOWfMWFtt4r+CAUD12j6I1rMn26/JeSJgnHUsGiCf3QY0kbDv/vEusfFeMQwRupj0OGPXzwYySwc4B/TiI1UWvkNY0Ie8pVfHtEdpWl/LCNbXMR9NVSeLBHERsMvyUjMmvwAv6Qbho2GX/IjfcMOXYBfCiNbi0rDL8URXkrDLyUNdCtN+9IwatQeP0L55ZHNP2gNx6cyFnb/oOE3WHaVGk3/YhYsAoN3MMQwsMDRBDiqz0wXKSKSEEbMSEuAJ+aBSZOLbhOKxUAQZSXBKVeABbYkKFu2BWBLgrJl3445anWzTSyGwrG1SAIji6lAdTUBthgL3XGkBLqFuUAZ2fkkzm8xGGgjJQS2mAzU75BATJ5dS0rQ/5gNFEwTEAzDEboc7BZYTAcqpiUQDOOBd6YJdNAUg9MEJNUgXO8PQVKNh/T+UKRXI0KfsmOB6GpGep8qlmNIVPMogaSYEmBJE5AUY9Jwhq1GEnAI+tiSoJtAigTZCDCr20BMbU1g12XaCdIE3aFh32F6B1sp2BTdiMgGPGU3SLUZmr0V3XqgadkAqLohREfkgY1S0Q0h3UAfDADLjpAU9G0wACx7QlJo1mAAWLaFhH0caN0lIoyDQBPdIlVKI7BuQ9AkXDxoVi7Ik3THhq0hdi0pDwCzN8SWsNLIx+aQKnNmk1HYHWKThfYAuGwQOdXtbDLyEUsB8dnOhT/Rva7qAK5uEYHwxRoAl10i3XOGhj/RPRC/iBnoL+IpOhHEALjsEjkQHlr5s8OITjsD4ArxPRAemv4hqKLTxBsAVzAwbG0Vb9pGmhudzppgAFxRA4PCBwPgCgZGJxFOP3QcnUf/WRJAp5E9TKXhj4HRKZUMgCsYmAC/xMasiG9nFxhTAFeCHZ1OsMEAuMKWPAgODX/25HW+ZQMAC+sInW7ZAMDCSiLo6sBE+LOW0NlXDAAsrCbUQdElpgjrCY0zWBPpD1YUSjsDIAtrCqW9CfBnVaF0NACyRNaouBbJAMgSiZKqA2YAZIlEb/FEignIE9mzxXEU3VqViHzQ3gDQwpa90sEA0BKRj/xkAGiJyAc9sJktEpEPuhiPPAn5WK1YA15LQj5oZ4BrScgHHYxHnoR80NEA1pKQD3ownv5Irf/QITYpJSGfNQ7XGHkS8kGLAbklIR+0N07lQT5ots3ZakU+6GSAbSEkozSrIWjkIx9PHJrxLQadAcRlaOPLWRdBniH2fOYgdOo0Hgk060Xq43VBo3/QbO5ykOjBhwvOkv1Gj+ktHl9ul3vvFmZxqafHYjSLHxfjD8UFUxxLv2SKA1uyKa6Y4q0pbPF6Z4r3pvhgio+m+GSKB6OyKb6YEqwpQUwJzpTgTQnBlMDikTMiYFs2JRRTojUliinRmRK9KTGYwh54TKZEMDGbEospyZqSxJTkTEnelMRBk2hKSqYksDSbkoopgzVlEFMGZ8rgTRmCKQNL0WTKMJgygMPFlGxNyWJKdqZkb0oOpmQWqckUsFTxu5hSrCnsQRdnSvGmcMClRFMKS9jBlALu6746cRiQ2GokxLIDoZFRy14MEGxxGyzgazmwY4FdC9ZaANbiV1ug1YKnVlfFIKkFPi1esgU4LX6CBSItnrAFHC2IaHEPLFhoAUAL6lmcWQveWUDO4rta4M3qwSWAzOIRWCDMglsW99OCWBaVs3ibFmWzoJJFzWw7vEMNtMoSmLFgjo3U0KgcTqXFxlvCLxaf0mLhLd6kxbRboi0Wm26JtFhWJBaLbgmuWEy5xZW02HDLAsFivC0hFIsPaTHdHDnGweIPNThAYPEgLS61xVRbHEiLjba4jhbf0+IwWrxAq6ZdnQbMssVZshhki99oMcW2UAO30WKJLQ6jZcyJWuA08IdDDox5O9XEmBOIwN7zhxU/Y66OIlEGUQ+xnYdgzNU/FMa8HQxgzPW8FKEBrCp/9CAFPBhz1v2YR/7AgzFXh7CdyGLM1R0UxpzFOraLP9RgzFmEY48IZMCDMVfnjwU1hoU/1GDM9cCW+nksjzEL/KEGY67BQGHM9XycnngRxlwPuLGMBaj5Qw3GXE/Y6Sq0nQxjzNV90xNy6rfpolAdtnaEhDFXd01XZ+qn6bpLhvynn3/+2fyrTnNypv2LpznbNYW/c62kn9jk8PZuQ8HphHNLuTe/4Nb9CMijOxC3rz8cLoYI1zHu3gz5Zrl7t9od3TRpCUev7Anz7Ynnq4/7k80bjl9zABWyZdr7lj7Q/Mv1m82CqEijj95P9uMtZ9yTnjBeflzfPTF+st+fkM7B7fP1j+vr9XZzvRijnp4m5+iFT5avV9PtF/gp3TgEOCj99OLieqXXUsDZnjiL7VXu9dm7J6vNGy7V2PuWQ846BlNVbQvHzm9X219OB83nIhN7Dia/+rdpoQ7n/6GFf/y3aeE8QL9wDB8sd0eXih4sd5NS6BFsrmkxTS9Pn7U58Wi3/Kldxmj006v94eJHI/rdj0b06x9Pr/aP2nn7drGMQ/BMIp1GT6/2HRJoxNOr/WO9iNWLPl73ewF3TtBTQBPP13tuk0308+32Ug/Qk9AurjzcbvbbD7vrfmvhZN/FuYWYJ/s9k1hB6u9ggfuFYMBcofntyhGNhOJGxH2lvtqcf7Xbbfs9Lma2klocVo8/bM46LJAJeYRikH0IyeUiTC9M+yF7YeY95NF4P1m9WW3Oj2/bIF1LPUJYXnRInHhPF/B4xQQLc0H6r6uiWVx8zV2I1fUtDO+pz66WZ9wFUN7zrbijNsxX4noaMs7lbkozF52Sp6K3WGu52+0+SjxcNfp6fY1CHstDEu/r4iRLs6dyE+PWO1PRKbUXvCUNpb5Zb9bvP7z/42q3PVz1IOPGRUWF+Hbr5XS3uljtfvfkULqlH3VcSzhuJpIepx7a2VIfrS6+XozRMmhzyovFONxMeblo133mIq96wunyWPlOlzd0C+Zz0oGzJt29oHm6PL/RdPrudHl+98Ln6fL8M3c+T5fnKPvLQ/f0lFc3UrCN/a4SDNdn7/pNpdPlVbtZ+bKDxpzwasH25OLi2dlutdo8Xp4p+iAesHbU/ZDMgyO1Jel4PKZaR/OHIpCH2UPKpEBNrVqZ3XtQdJHuB71IQ6JiRuEaI5RCDIvTRr6Y7juR9zWEaDv2u/XVo9XZ+v3y8nq+VKSQ3L0dN7sMR63TAreap2nH7UOPNPGogTM9m6rmpOCAtIvAc60ZArVOa06ZS9Oe6f00p4HP9vebzWr3Hc2jJFNNX3u9GH/gItG9e/wRd09Y9N8L9/JEp3vhniQoZykQepEDpUnH6c7eLD7nZdvzsr0X7L1wL2hJ5f0v+vMnwGW1PF/tsNZ6Z0q7baYer/ePJ6WJXWn0JhajOOeoLmmvnS0vtTKj/1/b9YbEyRF4uLw6Jp+v38/u5JBzkRwUMX7/fvlmxYtmgH+43Jxfrl68XV+/W+2+W27e9KvLLf3B9mNPa6PXUlWSo8ubf1hvL9ebKbVfXGxFH653Z5e30b5ncesUoY8M4Etc7q8+Xr08dnumxFfHia8+V3JKvFGSgt8sPz5av9Gr9Sjh093+7fbh8v1qt+zoc3e5dro8/6fdv2PEvrhiO12ef2m9xkqKHjpCLMgJfHruUf/dsox3V2G34JrO+AxY/+eCPuD1j13QH+6XWEop2TmbCZctuMH/+fv6Prv7NltrrY8xU/TW3f2U7zvNt9kT3NAC07t8Cvdjorak5MpwlMt1+7u5N6/4f0bQL974vyvordv/dwU9/hTAXVGm3M8LevMjAVzjdzZ7773LIbukJqZf87/9FQH12q1LMrjoiss5x+Eznw24bwlWBpvFF0Knh88GfOErAkPw1sM6lRyYh4fL//etZcBLKcOQQnZRx2n+rsBns6cvDdj7oTCC1lrnfBQdxP7lAb5kcTuvfVLA3g9NK6wNgS8jYFraJw3ulzs5N75U8Flpbny74Eslvsx6/rrB57j/C792AOD952sH89dkfvHXDm4tr//x+l/4PsKf5Th+ePvrCS33H/1Gwj/vZn33P25crP/TnWjs/wvz/h/rPn2T57C8/WWf3wk2Htv3X61xvyvnr9S2y32bXMB45hhCRkV/RZY93rDsuO3/sF2/a7v/tXb9rqj/plYdZ+w/Vv1XZNWJw7RQxWe+e9Ryfy1WnRjLs7Pd+oqP5mmIhU//XLZP/2w2q7P9IbqoX5RqMS79+M7P/wsVtwnN').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779222858453', 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_1779222858453();\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
}