{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "79d9f6e4",
   "metadata": {},
   "source": [
    "# rf708_bphysics\n",
    "Special pdf's: special decay pdf for B physics with mixing and/or CP violation\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": "6ad86782",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:51.102489Z",
     "iopub.status.busy": "2026-05-19T20:34:51.102364Z",
     "iopub.status.idle": "2026-05-19T20:34:51.118223Z",
     "shell.execute_reply": "2026-05-19T20:34:51.117719Z"
    }
   },
   "outputs": [],
   "source": [
    "%%cpp -d\n",
    "#include \"RooRealVar.h\"\n",
    "#include \"RooDataSet.h\"\n",
    "#include \"RooGaussian.h\"\n",
    "#include \"RooConstVar.h\"\n",
    "#include \"RooCategory.h\"\n",
    "#include \"RooBMixDecay.h\"\n",
    "#include \"RooBCPEffDecay.h\"\n",
    "#include \"RooBDecay.h\"\n",
    "#include \"RooFormulaVar.h\"\n",
    "#include \"RooTruthModel.h\"\n",
    "#include \"TCanvas.h\"\n",
    "#include \"TAxis.h\"\n",
    "#include \"RooPlot.h\"\n",
    "using namespace RooFit;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "71caae09",
   "metadata": {},
   "source": [
    "-------------------------------------\n",
    "B-Decay with mixing\n",
    "====================================="
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8c560740",
   "metadata": {},
   "source": [
    "Construct pdf\n",
    "-------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "acd27fd0",
   "metadata": {},
   "source": [
    "Observable"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "9154aca6",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:51.119818Z",
     "iopub.status.busy": "2026-05-19T20:34:51.119702Z",
     "iopub.status.idle": "2026-05-19T20:34:51.442260Z",
     "shell.execute_reply": "2026-05-19T20:34:51.441493Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar dt(\"dt\", \"dt\", -10, 10);\n",
    "dt.setBins(40);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bb92a534",
   "metadata": {},
   "source": [
    "Parameters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "71709110",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:51.444147Z",
     "iopub.status.busy": "2026-05-19T20:34:51.444013Z",
     "iopub.status.idle": "2026-05-19T20:34:51.650235Z",
     "shell.execute_reply": "2026-05-19T20:34:51.649459Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar dm(\"dm\", \"delta m(B0)\", 0.472);\n",
    "RooRealVar tau(\"tau\", \"tau (B0)\", 1.547);\n",
    "RooRealVar w(\"w\", \"flavour mistag rate\", 0.1);\n",
    "RooRealVar dw(\"dw\", \"delta mistag rate for B0/B0bar\", 0.1);\n",
    "\n",
    "RooCategory mixState(\"mixState\", \"B0/B0bar mixing state\");\n",
    "mixState.defineType(\"mixed\", -1);\n",
    "mixState.defineType(\"unmixed\", 1);\n",
    "\n",
    "RooCategory tagFlav(\"tagFlav\", \"Flavour of the tagged B0\");\n",
    "tagFlav.defineType(\"B0\", 1);\n",
    "tagFlav.defineType(\"B0bar\", -1);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e32504ae",
   "metadata": {},
   "source": [
    "Use delta function resolution model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "29866b54",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:51.652002Z",
     "iopub.status.busy": "2026-05-19T20:34:51.651882Z",
     "iopub.status.idle": "2026-05-19T20:34:51.860893Z",
     "shell.execute_reply": "2026-05-19T20:34:51.860407Z"
    }
   },
   "outputs": [],
   "source": [
    "RooTruthModel truthModel(\"tm\", \"truth model\", dt);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0947079b",
   "metadata": {},
   "source": [
    "Construct Bdecay with mixing"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "dcd4eab7",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:51.863039Z",
     "iopub.status.busy": "2026-05-19T20:34:51.862920Z",
     "iopub.status.idle": "2026-05-19T20:34:52.071804Z",
     "shell.execute_reply": "2026-05-19T20:34:52.071258Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBMixDecay bmix(\"bmix\", \"decay\", dt, mixState, tagFlav, tau, dm, w, dw, truthModel, RooBMixDecay::DoubleSided);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c68bc146",
   "metadata": {},
   "source": [
    "Plot pdf in various slices\n",
    "---------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "435c4dc8",
   "metadata": {},
   "source": [
    "Generate some data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "c300db68",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:52.073506Z",
     "iopub.status.busy": "2026-05-19T20:34:52.073381Z",
     "iopub.status.idle": "2026-05-19T20:34:52.282614Z",
     "shell.execute_reply": "2026-05-19T20:34:52.282147Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_65: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({dt, mixState, tagFlav}, 10000)};\n",
      " ^\n"
     ]
    }
   ],
   "source": [
    "std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5bab2ae9",
   "metadata": {},
   "source": [
    "Plot B0 and B0bar tagged data separately\n",
    "For all plots below B0 and B0 tagged data will look somewhat differently\n",
    "if the flavor tagging mistag rate for B0 and B0 is different (i.e. dw!=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "0b62c672",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:52.284082Z",
     "iopub.status.busy": "2026-05-19T20:34:52.283964Z",
     "iopub.status.idle": "2026-05-19T20:34:52.494471Z",
     "shell.execute_reply": "2026-05-19T20:34:52.494035Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_66:4:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame1, Cut(\"tagFlav==tagFlav::B0\"));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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_66:7:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame1, Cut(\"tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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 = dt.frame(Title(\"B decay distribution with mixing (B0/B0bar)\"));\n",
    "\n",
    "data->plotOn(frame1, Cut(\"tagFlav==tagFlav::B0\"));\n",
    "bmix.plotOn(frame1, Slice(tagFlav, \"B0\"));\n",
    "\n",
    "data->plotOn(frame1, Cut(\"tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bmix.plotOn(frame1, Slice(tagFlav, \"B0bar\"), LineColor(kCyan));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e19fed54",
   "metadata": {},
   "source": [
    "Plot mixed slice for B0 and B0bar tagged data separately"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "d01e1d22",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:52.495958Z",
     "iopub.status.busy": "2026-05-19T20:34:52.495840Z",
     "iopub.status.idle": "2026-05-19T20:34:52.706284Z",
     "shell.execute_reply": "2026-05-19T20:34:52.705856Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_67:4:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame2, Cut(\"mixState==mixState::mixed&&tagFlav==tagFlav::B0\"));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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_67:7:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame2, Cut(\"mixState==mixState::mixed&&tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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 = dt.frame(Title(\"B decay distribution of mixed events (B0/B0bar)\"));\n",
    "\n",
    "data->plotOn(frame2, Cut(\"mixState==mixState::mixed&&tagFlav==tagFlav::B0\"));\n",
    "bmix.plotOn(frame2, Slice({{&tagFlav, \"B0\"}, {&mixState, \"mixed\"}}));\n",
    "\n",
    "data->plotOn(frame2, Cut(\"mixState==mixState::mixed&&tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bmix.plotOn(frame2, Slice({{&tagFlav, \"B0bar\"}, {&mixState, \"mixed\"}}), LineColor(kCyan));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8d3e667b",
   "metadata": {},
   "source": [
    "Plot unmixed slice for B0 and B0bar tagged data separately"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "3868f5db",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:52.707769Z",
     "iopub.status.busy": "2026-05-19T20:34:52.707640Z",
     "iopub.status.idle": "2026-05-19T20:34:52.918774Z",
     "shell.execute_reply": "2026-05-19T20:34:52.918302Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_68:4:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame3, Cut(\"mixState==mixState::unmixed&&tagFlav==tagFlav::B0\"));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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",
      "In module 'std' imported from input_line_1:1:\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/initializer_list:45:11: error: 'initializer_list<value_type>' (aka 'initializer_list<pair<RooCategory *const, std::basic_string<char> > >') layout not recognized. Must be a non-polymorphic class type with no bases and two fields: a 'const E *' and either another 'const E *' or a 'std::size_t'\n",
      "    class initializer_list\n",
      "          ^\n",
      "input_line_68:7:1: error: reference to 'data' is ambiguous\n",
      "data->plotOn(frame3, Cut(\"mixState==mixState::unmixed&&tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
      "^\n",
      "input_line_65:2:30: note: candidate found by name lookup is 'data'\n",
      " std::unique_ptr<RooDataSet> data{bmix.generate({dt, mixState, tagFlav}, 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",
      "In module 'std' imported from input_line_1:1:\n",
      "/usr/lib/gcc/x86_64-redhat-linux/14/../../../../include/c++/14/initializer_list:45:11: error: 'initializer_list<value_type>' (aka 'initializer_list<pair<RooCategory *const, std::basic_string<char> > >') layout not recognized. Must be a non-polymorphic class type with no bases and two fields: a 'const E *' and either another 'const E *' or a 'std::size_t'\n",
      "    class initializer_list\n",
      "          ^\n"
     ]
    }
   ],
   "source": [
    "RooPlot *frame3 = dt.frame(Title(\"B decay distribution of unmixed events (B0/B0bar)\"));\n",
    "\n",
    "data->plotOn(frame3, Cut(\"mixState==mixState::unmixed&&tagFlav==tagFlav::B0\"));\n",
    "bmix.plotOn(frame3, Slice({{&tagFlav, \"B0\"}, {&mixState, \"unmixed\"}}));\n",
    "\n",
    "data->plotOn(frame3, Cut(\"mixState==mixState::unmixed&&tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bmix.plotOn(frame3, Slice({{&tagFlav, \"B0bar\"}, {&mixState, \"unmixed\"}}), LineColor(kCyan));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8f824d02",
   "metadata": {},
   "source": [
    "-------------------------------------------------\n",
    "B-Decay with CP violation\n",
    "================================================="
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5b0edb1f",
   "metadata": {},
   "source": [
    "Construct pdf\n",
    "-------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d2f180db",
   "metadata": {},
   "source": [
    "Additional parameters needed for B decay with CPV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "3bd14523",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:52.920288Z",
     "iopub.status.busy": "2026-05-19T20:34:52.920170Z",
     "iopub.status.idle": "2026-05-19T20:34:53.137028Z",
     "shell.execute_reply": "2026-05-19T20:34:53.136458Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar CPeigen(\"CPeigen\", \"CP eigen value\", -1);\n",
    "RooRealVar absLambda(\"absLambda\", \"|lambda|\", 1, 0, 2);\n",
    "RooRealVar argLambda(\"absLambda\", \"|lambda|\", 0.7, -1, 1);\n",
    "RooRealVar effR(\"effR\", \"B0/B0bar reco efficiency ratio\", 1);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b4fc485f",
   "metadata": {},
   "source": [
    "Construct Bdecay with CP violation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "b8ed3476",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:53.138860Z",
     "iopub.status.busy": "2026-05-19T20:34:53.138741Z",
     "iopub.status.idle": "2026-05-19T20:34:53.347882Z",
     "shell.execute_reply": "2026-05-19T20:34:53.347297Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBCPEffDecay bcp(\"bcp\", \"bcp\", dt, tagFlav, tau, dm, w, CPeigen, absLambda, argLambda, effR, dw, truthModel,\n",
    "                   RooBCPEffDecay::DoubleSided);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "688a1ed5",
   "metadata": {},
   "source": [
    "Plot scenario 1 - sin(2b) = 0.7, |l|=1\n",
    "---------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9654838e",
   "metadata": {},
   "source": [
    "Generate some data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "18b1d054",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:53.349658Z",
     "iopub.status.busy": "2026-05-19T20:34:53.349518Z",
     "iopub.status.idle": "2026-05-19T20:34:53.558413Z",
     "shell.execute_reply": "2026-05-19T20:34:53.557896Z"
    }
   },
   "outputs": [],
   "source": [
    "std::unique_ptr<RooDataSet> data2{bcp.generate({dt, tagFlav}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c5fbc63e",
   "metadata": {},
   "source": [
    "Plot B0 and B0bar tagged data separately"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "e5f9ba1a",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:53.560214Z",
     "iopub.status.busy": "2026-05-19T20:34:53.560097Z",
     "iopub.status.idle": "2026-05-19T20:34:53.769328Z",
     "shell.execute_reply": "2026-05-19T20:34:53.768756Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 4495 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcp) plot on dt represents a slice in (tagFlav)\n",
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 5505 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcp) plot on dt represents a slice in (tagFlav)\n"
     ]
    }
   ],
   "source": [
    "RooPlot *frame4 = dt.frame(Title(\"B decay distribution with CPV(|l|=1,Im(l)=0.7) (B0/B0bar)\"));\n",
    "\n",
    "data2->plotOn(frame4, Cut(\"tagFlav==tagFlav::B0\"));\n",
    "bcp.plotOn(frame4, Slice(tagFlav, \"B0\"));\n",
    "\n",
    "data2->plotOn(frame4, Cut(\"tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bcp.plotOn(frame4, Slice(tagFlav, \"B0bar\"), LineColor(kCyan));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d16bd42b",
   "metadata": {},
   "source": [
    "Plot scenario 2 - sin(2b) = 0.7, |l|=0.7\n",
    "-------------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "71c1f026",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:53.771171Z",
     "iopub.status.busy": "2026-05-19T20:34:53.771053Z",
     "iopub.status.idle": "2026-05-19T20:34:53.979860Z",
     "shell.execute_reply": "2026-05-19T20:34:53.979171Z"
    }
   },
   "outputs": [],
   "source": [
    "absLambda = 0.7;"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c1039240",
   "metadata": {},
   "source": [
    "Generate some data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "dd2e6997",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:53.981916Z",
     "iopub.status.busy": "2026-05-19T20:34:53.981799Z",
     "iopub.status.idle": "2026-05-19T20:34:54.190595Z",
     "shell.execute_reply": "2026-05-19T20:34:54.189991Z"
    }
   },
   "outputs": [],
   "source": [
    "std::unique_ptr<RooDataSet> data3{bcp.generate({dt, tagFlav}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ab2ff074",
   "metadata": {},
   "source": [
    "Plot B0 and B0bar tagged data separately (sin2b = 0.7 plus direct CPV |l|=0.5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "2f182a04",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:54.192647Z",
     "iopub.status.busy": "2026-05-19T20:34:54.192516Z",
     "iopub.status.idle": "2026-05-19T20:34:54.401672Z",
     "shell.execute_reply": "2026-05-19T20:34:54.401105Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 3617 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcp) plot on dt represents a slice in (tagFlav)\n",
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 6383 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcp) plot on dt represents a slice in (tagFlav)\n"
     ]
    }
   ],
   "source": [
    "RooPlot *frame5 = dt.frame(Title(\"B decay distribution with CPV(|l|=0.7,Im(l)=0.7) (B0/B0bar)\"));\n",
    "\n",
    "data3->plotOn(frame5, Cut(\"tagFlav==tagFlav::B0\"));\n",
    "bcp.plotOn(frame5, Slice(tagFlav, \"B0\"));\n",
    "\n",
    "data3->plotOn(frame5, Cut(\"tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bcp.plotOn(frame5, Slice(tagFlav, \"B0bar\"), LineColor(kCyan));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "99961e70",
   "metadata": {},
   "source": [
    "---------------------------------------------------------------------------\n",
    "Generic B decay with  user coefficients\n",
    "==========================================================================="
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e67af61b",
   "metadata": {},
   "source": [
    "Construct pdf\n",
    "-------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7d3acc73",
   "metadata": {},
   "source": [
    "Model parameters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "95339e66",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:54.403426Z",
     "iopub.status.busy": "2026-05-19T20:34:54.403313Z",
     "iopub.status.idle": "2026-05-19T20:34:54.612029Z",
     "shell.execute_reply": "2026-05-19T20:34:54.611438Z"
    }
   },
   "outputs": [],
   "source": [
    "RooRealVar DGbG(\"DGbG\", \"DGamma/GammaAvg\", 0.5, -1, 1);\n",
    "RooRealVar Adir(\"Adir\", \"-[1-abs(l)**2]/[1+abs(l)**2]\", 0);\n",
    "RooRealVar Amix(\"Amix\", \"2Im(l)/[1+abs(l)**2]\", 0.7);\n",
    "RooRealVar Adel(\"Adel\", \"2Re(l)/[1+abs(l)**2]\", 0.7);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "725081e5",
   "metadata": {},
   "source": [
    "Derived input parameters for pdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "58cadac0",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:54.614034Z",
     "iopub.status.busy": "2026-05-19T20:34:54.613921Z",
     "iopub.status.idle": "2026-05-19T20:34:54.822667Z",
     "shell.execute_reply": "2026-05-19T20:34:54.822053Z"
    }
   },
   "outputs": [],
   "source": [
    "RooFormulaVar DG(\"DG\", \"Delta Gamma\", \"@1/@0\", RooArgList(tau, DGbG));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "20444024",
   "metadata": {},
   "source": [
    "Construct coefficient functions for sin,cos,sinh modulations of decay distribution"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "83954752",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:54.824757Z",
     "iopub.status.busy": "2026-05-19T20:34:54.824626Z",
     "iopub.status.idle": "2026-05-19T20:34:55.035328Z",
     "shell.execute_reply": "2026-05-19T20:34:55.034702Z"
    }
   },
   "outputs": [],
   "source": [
    "RooFormulaVar fsin(\"fsin\", \"fsin\", \"@0*@1*(1-2*@2)\", RooArgList(Amix, tagFlav, w));\n",
    "RooFormulaVar fcos(\"fcos\", \"fcos\", \"@0*@1*(1-2*@2)\", RooArgList(Adir, tagFlav, w));\n",
    "RooFormulaVar fsinh(\"fsinh\", \"fsinh\", \"@0\", RooArgList(Adel));"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8f62661d",
   "metadata": {},
   "source": [
    "Construct generic B decay pdf using above user coefficients"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "ee5b9907",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:55.037552Z",
     "iopub.status.busy": "2026-05-19T20:34:55.037432Z",
     "iopub.status.idle": "2026-05-19T20:34:55.246229Z",
     "shell.execute_reply": "2026-05-19T20:34:55.245618Z"
    }
   },
   "outputs": [],
   "source": [
    "RooBDecay bcpg(\"bcpg\", \"bcpg\", dt, tau, DG, RooConst(1), fsinh, fcos, fsin, dm, truthModel, RooBDecay::DoubleSided);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f3bbc66a",
   "metadata": {},
   "source": [
    "Plot - Im(l)=0.7, Re(l)=0.7 |l|=1, dG/G=0.5\n",
    "-------------------------------------------------------------------------------------"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "84ac0742",
   "metadata": {},
   "source": [
    "Generate some data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "1abb9fc0",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:55.248325Z",
     "iopub.status.busy": "2026-05-19T20:34:55.248209Z",
     "iopub.status.idle": "2026-05-19T20:34:55.456995Z",
     "shell.execute_reply": "2026-05-19T20:34:55.456363Z"
    }
   },
   "outputs": [],
   "source": [
    "std::unique_ptr<RooDataSet> data4{bcpg.generate({dt, tagFlav}, 10000)};"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8abc058a",
   "metadata": {},
   "source": [
    "Plot B0 and B0bar tagged data separately"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "01b02d71",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:55.459065Z",
     "iopub.status.busy": "2026-05-19T20:34:55.458949Z",
     "iopub.status.idle": "2026-05-19T20:34:55.668208Z",
     "shell.execute_reply": "2026-05-19T20:34:55.667919Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 4991 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcpg) plot on dt represents a slice in (tagFlav)\n",
      "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 5009 events out of 10000 total events\n",
      "[#1] INFO:Plotting -- RooAbsReal::plotOn(bcpg) plot on dt represents a slice in (tagFlav)\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "input_line_130: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": [
    "RooPlot *frame6 = dt.frame(Title(\"B decay distribution with CPV(Im(l)=0.7,Re(l)=0.7,|l|=1,dG/G=0.5) (B0/B0bar)\"));\n",
    "\n",
    "data4->plotOn(frame6, Cut(\"tagFlav==tagFlav::B0\"));\n",
    "bcpg.plotOn(frame6, Slice(tagFlav, \"B0\"));\n",
    "\n",
    "data4->plotOn(frame6, Cut(\"tagFlav==tagFlav::B0bar\"), MarkerColor(kCyan));\n",
    "bcpg.plotOn(frame6, Slice(tagFlav, \"B0bar\"), LineColor(kCyan));\n",
    "\n",
    "TCanvas *c = new TCanvas(\"rf708_bphysics\", \"rf708_bphysics\", 1200, 800);\n",
    "c->Divide(3, 2);\n",
    "c->cd(1);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame1->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame1->Draw();\n",
    "c->cd(2);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame2->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame2->Draw();\n",
    "c->cd(3);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame3->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame3->Draw();\n",
    "c->cd(4);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame4->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame4->Draw();\n",
    "c->cd(5);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame5->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame5->Draw();\n",
    "c->cd(6);\n",
    "gPad->SetLeftMargin(0.15);\n",
    "frame6->GetYaxis()->SetTitleOffset(1.6);\n",
    "frame6->Draw();"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4f07f925",
   "metadata": {},
   "source": [
    "Draw all canvases "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "d73ebdc5",
   "metadata": {
    "collapsed": false,
    "execution": {
     "iopub.execute_input": "2026-05-19T20:34:55.670121Z",
     "iopub.status.busy": "2026-05-19T20:34:55.670010Z",
     "iopub.status.idle": "2026-05-19T20:34:55.902883Z",
     "shell.execute_reply": "2026-05-19T20:34:55.902367Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "\n",
       "\n",
       "<div id=\"root_plot_1779222895843\" style=\"width: 1200px; height: 800px; position: relative\">\n",
       "</div>\n",
       "\n",
       "</div>\n",
       "<script>\n",
       "   function process_root_plot_1779222895843() {\n",
       "      function execCode(Core) {\n",
       "         Core.settings.HandleKeys = false;\n",
       "         \n",
       "Core.unzipJSON(28384,'WkwIMx0A4G4AeAHtnW2TVDeS778KUbEvNR1KPR3p6FUDZj33YtNh8AA7MTFRdFdDXZqq3urCxrPh777xS+lUVT/g8Mw1dyEuE1Ptk3o4mZJS/0ylpMN/zf6+/eVysZq/W8zG2bMH89VP86vni1dPV/PLqzfr7czMzn9cLf/z/eLPD2ejNbPz+8vtVXt68ur/LE63pM8o9uRyu1yvOvG/l6uz2ejN7Hz3pvG/7uL1MQbeB+dDMLPzx8vV4sH6Yr2ZjdLJp9tfLhZ78vnybPumkY+WFxe9MMJCToWt1eqL8+13883r5Wo22iNSfli+fnMj6f56u12/u17s2fryesKL8yVCODM7f7l/PG6PvPjF1Xa+hUsplLlGHTeKyo8283eLm3KTdqPhu3LXG7QrOiVPr6T6Yc/wyvvrzdli83T5j957B4nfrc8WbVxfSPvvy/7fF6717cv+3xfb9fGrq5Plh8XF215ju76LnI0iJdH2azWGQftsX+d6wmz8U0v48VothvPHfZ1DcsfnL9dqtLf8ZV/nesKOz9SY7fpFb4825jY5G/ORT+JssOL9MBS/+BMaun/By7etq/oLIBF0Imfjn+yRteJK9D4NeYgpDAOqcrH++fuHD3q/HxIvfrzUDIb15cHz813qt7un41dX1150/Orq2ruOX13tqx2/utrX/PHDO50Q9PEv+8cP7+YfWoN+/GX3+OzNYjufjV4b9mbZn46vLhen2x/m2+W6teL79+9eLTbt+dny9O2H/eMv7fHx+nVPfLx+vU/7R8s9mZ+dzJcrJpCZnT/YrK+u3syX/YU78mTdwehQs1GvRu+V+rv12fJ8uTibjefzi6uFmZ3/+2Z59uE6+cuePH519WC93hyU/+ZsuZ2/YspvN+95waPlh8XZtXZPrz7ZLN8tt8ufFle3UO/x8gpQnQC3k/PNZjb+9W9mtr7c8vCrmZ1/82FxejUbV+8vLszs/PuG0Jvzwea/v7p888vV8lTB89lyi0x35Hz//t3J/GKx3U7wSUd+v/iwvZ368M9PTx4fv5yNs3+bHs3s/OH6/auLxf335+fTQP6w2M6XK3qx98GLq+U/Fj9eTfkvr5Oa+8NifjEbHcw1u9Hijzz/M7Pz58vV2frnZ+vLF22sd/TLQ7pDGa9pBb5dANxdPX7eYcCDN7M+0R/Mt9tb/X+83TYjRwNf3F9sf14sVh3Fr1Har48263fP1pezUY5Qqhdn8y0gqcTLiaANx42QX83s7XfrnxZPLuf/+X6nKm9/WNAx1xPPv12+fvOYNnSDpWo7356+mXr37dM365+/+Wmx2j7dzrfvr3ba+fb4/XaNfuxKfrdYvb8/3zQa7Tk+Rf92Nc5/WMzPnqwufplqnD9fbt+s328PVXVS32/nV135ppTDUn+9Ycf/MG8BKP2ot/B88UrhYLl6/TGXQXHiYn511ecK5ZqPcphwia7O7OhiNP1XZbTGGludpvLkR6u5tgbNczHWuKtja+olSR96WZ7zmItxIiaHWnjOwThJVWAYTP9VkVGKM5KdkcFVceOQjP6/ih/FOtN/VcIobjD9VyWOEr3pvypplCGb/qsyjM4G039V8uhcMf1XpYwuONN/1VktLIX6tjq5TrpRijdSxEjK1flRcjKSIX11ofFN0Ui01cVRsifHSIzVIVUXMtjqhlFiNBKcEW+ry6P4qA3MqboyyhCNswlJqrejeGeUdfLVyyh+QAIjyVbvRvHWSMyGF3vf+OTc6tJVrvGJvvo4igxG64dYfRrFBiRQkf0wii1GK3hXfR7FBe1XZ0v1Zd/rMdZgRxms0eYWqUFGhlhl9rEGNzonRvsHkq7KTUhvawijdjGj5kINcdQ+RQwpNaSR15RgSqphaGzoqDDUkEflKCpmDWWUMBhthLgarQqhKlajjBKSUZ2t0Wk/dMKPJU/PYRzC9BzHKNNzGt0wPQ+jnx4zit1fU0ZrAgoca7KjNSm3Z6ZOkfbsRhql06omZpD40nKYRGiHVo8QeWhEgoc0Jklnku9clLtOUqkJ9kq4oQ7wV8KGOiAAhGRbByRQIqY6THOY8R2QQHNsrgMSQORYBxUgRpOkDtNM9kMdJvbi61DGJqGxNduxCcizjFGbwbMbs3YDz0zhhi62ZrRSuwGiTV/42ZqZJmkimLpNKlszE1e7AYJpq71tbC0duHyGkIZIEiAaeEkBvorXHEkJIjQiOIiGYiiVraXBWNE3D1pqENKzPoeB56LPzhtbxdpRREyICgpiRTU0OBNcFesaKBSGV6pYP+bBiGoN2UHbITEZDxlH57Lxg4mhik2jRGsGJmipYoGMpGUl8SqmpzPeG+cg6RKT6FWpotgamUOuigBjYrwYB+W0VwGjVKqArTzzc7EK4Mpo0vG8V2KTuOmzgK5WTClGrK0iwwjA0R7JVUBX9C4VIyB6GSUXGEsVZ0cJYFYwksB3cKz1W4Z0o1azg8mlivOj88rHVgFac8QsGKa5gK3ijUKZo3DSqe4kN5lB1wMT4KyMRZiSRapjSIDE9qvOgk3e9F91jInQHv1VZ2PDW2C9OgYle4UAW50dRubmIGbgzXl0vjSxZKjOlpFagHWsjiFBC6hcncg4iGmVqxOnuNxrV8egwI/K1UmgZqtcnURqtsrVMR4gugzUrk6GJu3ENO8ap1yLyjuxxdwBR03o6nRIkLbxZUhaWyX76nRItLHOh+pc2LVWOTMovbUI7UDxfa+6YWpvq5x3Daayn+Yw2GZD9dM0RpmwW8zkbsUb3Saz5ltXPfN5ylf6wCdQus1qykepnondiyuZdx6EkmBbczei1GD3jpGtwTafgOpQblcTau+nQAXMCRalBktvNPSifUEns2K/ticwmaGxxtHXgOpYt2tfENto8q2rAWcJEbGgSjt49eYF8UrZoPKLykFRWiNq5Kis1AT1tgZRqG/tEqCeKqRjaKhga3DYmdYk14ek+Yw14A71JmHSApO2N0mbjH70JrV8kMRpk1s+3kdvEl3g8tSkVrxMbdLS3kLCQAt7UbKNWg1+Mnpa1avVa+Ndg1fDu6uolndXj+7QQasRlew6wltikd1QN3o/9EgUy37wW35o7WltrBGVbHxU5FgSIrdRq7EME6Xv0ra3BtZYtOlQ1tVkGQK1ijXZ5qVTJ1narINWk6XFOmg12eat4zGghW2GqYgJ+GozrNNZm9yHrSZsCvlt2GpSCFMNbCzRwqaB2qSEadEZ1l4nXl9HExEQy6JK2Au3/piaJYfLh9S9dlCPmm129MZhVdDG1h1OFYFRqwmVnNpHPRCrtw8tSUBWm2Gd7oa/D39SzGrt065WL11nWKebeWsKklyZmgdvvPTWOqXa/Owl/X4RlXzrFDyShHfOHNPn1h3tWVWDx2ECwiZ2HToQ9m6vA1rXRYxSB5Suk1SOu/Zq5n74ad7Q3RsVoQ7dwYGTFm6ta5l5wkBFhwwG9n6NUrNtSKCDUnMHQQWJbHdjTk4Xp2FP7trY5a150kZFmowudt1THrvOxsdsiNxYdCWkMAUxpR1WKKlI2F+5B0JyOhy15mbs6ATI1tWMHW1qr4XzjlIeOyzCyVUFhD1Etw6td2rGhE64zmsBSbStgXEGIztJ5e4Jd5G6NvZeqLlr48SpzVUoRGLN2BnhbHfsakKhjL1tTIvMkrF3EdMi+9ZJTW9q3gMkEnu1F8jUynbU7oY6+2l90ArvMVsZ6RJlqhoO1yg5KGLzHn1vUMRuk6/moCulridhAmzaGfbiIF3Yi9NesxOHhoWdOK3wJA4vipM0SvQeaia6ZpaLfVy0IXGnRvreiEQdnXNbNfYOrTk2iXqH1xyRSHMJhUSVaE8iUcfHrGvIbg9yagatDVvNbR3ZBq1mXUh24MxtIdmkrTk1895UoeZ0YGThr8vJ/Vtb/0yt1hXlrmpbUk55Qw+FTAIN3fpNIg1NsadOGZpiN8TOQ1PsqaFDc4Z2vTCwsplW2Hloij11BAvNlolG5abYzfzl3IZtKprbqDWjkbOOWbMSOasP1MVhxam61QaXNWfHTO2j3Kb+xKNN/R0PnfqdRem2R6t1pJ56oAN159+dAzVTuYN0JxooTqPZMXo3nh2jp9fsDJAqdMdofVHpEN2JPULTb+UAotHgAkh3W2xdLR2mG5fScbq3pEw43Ye6dJxunVm6zzBRzWPoI106UPfMHl1TnUIolht0fPvVAliix+1Xi9MRbDa0uDBGz7JaPFVVzQNL46EWl8ZodQUeWPQPYygs1CXw0jyGbFh4hlSLK2NIJlojIdfi7RiATyNRavEyhmBUmFCLd2PwRuFpqMX7MTjD8jsRSQhjsIaIW3K1+Dj6YmIhFlaLT6PPJrUQYfHD6AfDEnmAQx59NCyABzgUwiYpGRlSLcGO3ps0GBlKLUFG71jNszYuwY1el3ySYy3Bj64Y1sJ5qCWE0WUNEhTCGpGw1RCJFtZCJC2ZYTBS4DAQIRiKkQIHFqsmo8lSSyijcyY3v65oHM3o0m+oBXAsJkdiUrUAjdnkZJz4WiLRA5OzcRJriWGUZArRq1xL1AUhgTAntYCL0RTXKaa8hicc78S/I+TSKF02l9SohF9lytApVMYU1vmEaogZEItwUmpJHkoIXCipwqjWKhlbru+5qZEtcFHS0EgamWtJmWaJpZWQBGQNoQklB0KfhlCpk6GWQZQkHKYkIR9DSKiR2kMirpNEiQ3RC0LgBWTMRjTWAEl0ocULlRwaCV9yWcwbIcpBb4OMkIgRa8lWSSIFSkojEYNcYtca9GikbyRihFoIyJFL8yEJyBpxNB8yNRIxIIdGIoavBYCkMGJAFiU9YvhaSpOK2JO4WlhEZaKPnWx9RbxJpBZQkug0YkC2EfSIAdlG0CMGMThVKPGIAakapdFVJVWlJKDfpRaW9MEIUSxLLI2gHTQKn6EJRgEYE01U3UggHKH5fhQi0Mk4O1CeCCYR6YmOjc4TTVjGSChExCk/KB0RR2msCrBjnI3kazhVIvJAix1LN8VKipK6jUGuGwmmt1iJWPFKIovmBiURxZMbR0I6SKJkgkwIQsSRAB7xxInMSiKG5paRaB69QojRWa2r0UFIgRFwVuhT5xCS7Y9Cl6mrq+FWKXQZvq7rNG3C2yUGCTohtgbxOo2g6u62bY6CLISMkBz8UmGa+itGwp4FWGHbxkiGv8dj6TT8WYSxtaS7EtBEGFtMMSOPJxbbaeTxLA5bXJT4pPW6I6Bx0hbJHQBeITyb6Q2fgWU2bkBwscA7/KGRB4CHf8N7sUA8/KHhD8jDHwsA/+AxFi2OSX+EgClpNPKEiKER8HyAf0iYIcn0OfzDgJFqtPLPmLC2W0V/hIKBazT9ES3mT7ffEvJEwThqWDTBPzoMaaPh331i3eNiPGIYI/Ux6PDHLx6MZBYO8I9pxEYqrfyGMSEP+covj+iO0rQ/lpEtLuK+GipPlghio+GXZCRmTX6AX9INw0bDL/mRvmGHLsAvhZGtRaXhl+IIL6Xhl5IGupWmfWkYNWqPH6H88sjmH7SG41MZC7t/0PAbLLtKjaZ/MQsWgcE7GGIYWOBoAhzVZ6aLFBFJCCNmpCXAE/PApMlFtwnFYiCIspLglCvAAlsSlC3bArAlQdmyb8cctbrZJhZD4dhaJIGRxVSgupoAW4yF7jhSAt3CXKCM7HwS57cYDLSREgJbTAbqt08gJs+uJSXof8wGCqYJCIbhCF0OdgsspgMV0xIIhvHAO9MEOmiKwWkCkmoQrveHIKnGQ3p/KNKrEaFP2bFAdDUjvU8VyzEkqnmUQFJMCbCkCUiKMWk4w1YjCTgEfWxJ0E0gRYJsBJjVbSCmtiaw6zLtBGmC7tCw7zC9g60UbIpuRGQDnrIbpNoMzd6Kbj3QtGwAVN0QoiPywEap6IaQbqAPBoBlR0gK+jYYAJY9ISk0azAALNtCwj4OtO4SEcZBoIlukSqlEVi3IWgSLh40KxfkSbpjw9YQu5aUB4DZG2JLWGnkY3NIlTmzySjsDrHJQnsAXDaInOp2Nhn5iKWA+Gznwp/oXld1AFe3iED4Yg2Ayy6R7jlDw5/oHohfxAz0F/EUnQhiAFx2iRwID6382WFEp50BcIX4HggPTf8QVNFp4g2AKxgYtraKN20jzY1OZ00wAK6ogUHhgwFwBQOjkwinHzqOzqP/LAmg08geptLwx8DolEoGwBUMTIBfYmNWxLezC4wpgCvBjk4n2GAAXGFLHgSHhj978jrfsgGAhXWETrdsAGBhJRF0dWAi/FlL6OwrBgAWVhPqoOgSU4T1hMYZrIn0BysKpZ0BkIU1hdLeBPizqlA6GgBZImtUXItkAGSJREnVATMAskSit3gixQTkiezZ4jiKbq1KRD5obwBoYcte6WAAaInIR34yALRE5IMe2MwWicgHXYxHnoR8rFasAa8lIR+0M8C1JOSDDsYjT0I+6GgAa0nIBz0YT3+k1n/oEJuUkpDPGodrjDwJ+aDFgNySkA/aG6fyIB802+ZstSIfdDLAthCSUZrVEDTykY8nDs34FoPOAOIytPHlrIsgzxB7PnMQOnUajwSa9SL18bqg0T9oNnc5SHT//Tlnyf5NT+rNHl2s51vvZmZ2oafHYjSzn2bjX4sLpjiWfskUB7ZkU1wxxVtT2OL1zhTvTfHBFB9N8ckUD0ZlU3wxJVhTgpgSnCnBmxKCKYHFI2dEwLZsSiimRGtKFFOiMyV6U2IwhT3wmEyJYGI2JRZTkjUliSnJmZK8KYmDJtGUlExJYGk2JRVTBmvKIKYMzpTBmzIEUwaWosmUYTBlAIeLKdmaksWU7EzJ3pQcTMksUpMpYKnidzGlWFPYgy7OlOJN4YBLiaYUlrCDKQXc13114jAgsdVIiGUHQiOjlr0YINjiNljA13JgxwK7Fqy1AKzFr7ZAqwVPra6KQVILfFq8ZAtwWvwEC0RaPGELOFoQ0eIeWLDQAoAW1LM4sxa8s4CcxXe1wJvVg0sAmcUjsECYBbcs7qcFsSwqZ/E2LcpmQSWLmtl2eIcaaJUlMGPBHBupoVE5nEqLjbeEXyw+pcXCW7xJi2m3RFssNt0SabGsSCwW3RJcsZhyiytpseGWBYLFeFtCKBYf0mK6LZ66xWZbfHDLAQKLB2lxqS2m2uJAWmy0xXW0+J4Wh9HiBVo17eo0YJYtzpLFIFv8RosptoUauI0WS2xxGC1jTtQCp4E/HHJgzNupJsacQAT2nj+s+BlzdRSJMoh6iO08BGOu/qEw5u1gAGOu56UIDWBV+aMHKeDBmLPuxzzyBx6MuTqE7UQWY67uoDDmLNaxXfyhBmPOIhx7RCADHoy5On8sqDEs/KEGY64HttTPY3mMWeAPNRhzDQYKY67n4/TEizDmesCNZSxAzR9qMOZ6wk5Xoe1kGGOu7puekFO/TReF6rC1IySMubprujpTP03XXTLkv/3666/mU53mjL91mrPdVfiNuyX9xCbntzcrCk7HnFvKvd0LblySgDy4CHHzDsT+dohwJ+P29ZDv5pu3i83BdZOWcPDKnrC7QvFs8WF7vHrNGWwOoEK2THtk6QPNv1i+XnFCudMH7yf70Zpj7klPGM8/LG8fGj/ebo9J5/z22fKn5dVyvbqajVHgSM7BCx/PXy2mKzDwU7pxCHBQ+sn5+dVC76aAsz1xJ7ZXuZenbx8vVq+5WWOPLIecdQymqtoWzp7frLa9mM6b74pM7DmY/PKLaaEO57/Qwv/4Ylq4G6B/cgzvzzcHN4vuzzeTUugRbO5qMU0vTp62OfFwM/+53cdo9JPL7f7uRyP69Y9G9BsgTy63D9t5+3a7jEPwTCKdRk8utx0SaMSTy+0jvY3Viz5a9osBt07QU0ATz5ZbrpRN9LP1+kIP0JPQ7q48WK+26/ebq3554XjbxbmBmMfbLZNYQeo3sMD9k2DAXKH5XIo40qtrUFyJaNQ3q7NvNpt1v8zFzFZSi8Pq0fvVaYcFMiEPUAyyDyG53IXphWk/ZC/MvIc8GO/Hi9eL1dnhhRuGoKUeICwv2idOvKdbeLxigoVdQfqvq6KZnX/LXYjF1Q0M76lPL+en3AVQ3rurcQdt2N2L62nIuCt3XZpd0Sl5KnqDtZa72e6DxP1to2+XVyjkoTwk8b4uTrI0eyo3MW69MxWdUnvBG9JQ6rvlavnu/bv/WGzW+6seZFy7ragQ3269nGwW54vNvz/el27pBx3XEg6biaSHqft2ttSHi/NvZ2O0DNou5flsHK6nvJixD3BQ5GVPOJkfKt/J/JpuwXyXtOesSbdvaZ7Mz641nb47mZ/dvvV5Mj+74+LnyfwMZX+x756e8vJaCraxX1aC4fL0bb+qdDK/bNcrX3TQ2CW8nLE9OTt/erpZLFaP5qeKPogHrB10PyTz4EBtSTocj6nWwfyhCOR+9pAyKVBTq1Zm8w4UnaWjoBdpSFTMKLnLohDD4rSJ9ny670TJbyFE27HdLC8fLk6X7+YXV7tLRQrJ3dtxO5fhoHVa4EbzNO2wfSiJJh40cEfvTFVzUnBA2m3gXa0dBGqd1pyyK017pvfTnAY+6z+vVovNDzSPkkw1fe3VbPwrF4nu3eOPuHvCov9euJcnOt0L9yRBOUuB0IvsKU06THf2evFdXrY9L9t7wd4L94KWVN6f6M/fAJfF/GyxwVrrnSntth31aLl9NClN7EqjN7EYxV2O6pL22un8Qisz+v9rvVyRODkCD+aXh+Sz5budOznkXCQHRYw/v5u/XvCiHcA/mK/OLhbP3yyv3i42P8xXr/v95ZZ+f/2hp7XRa6kqycH9zb8s1xfL1ZTa7y+2og+Wm9OLm2jfs7h4itAHBvAFLvc3Hy5fHLo9U+LLw8SXd5WcEq+VpOB38w8Pl6/1fj1K+GSzfbN+MH+32Mw7+txerp3Mz/6w+3eM2Efv353Mzz62XmMlRQ8dIBbkBD4996D/bljG26uwG3BNZ9wB1l9v6QNev++WvshRiaWUMniXBo9R/eilfT+EoyGklAZniYXcvL/vh3zksrXWyuAioU8u+E8fAPCDPYqJXB9i89Wn3Nn4pztyr1/zv0PQj976vy3o9Q8A3CHowdcA7hBlyr1b0P01f74TwE1+L0VCLi7HLD5gFj/65QDu/buUckjFJycpaL/tX6lfDjiyzmZuRbscsuPDANOXA+7+kIBLJSfJ3oYi3jPT9vf/j6xlvBnyIYXsIvZ9/zmAo3gtm7k/fWrAHnmnd7O9jzHlyHunTw/wMQvG1lrrnI+CjO2bAvbIS681ZOsT7NoHBuxRUG2x1oYwtOjCtU8V3CnqtY8X3CHt7mMGdzHefdzgLt6f8GMHdNXXjx3sPifzT3/s4MbS+vfXv/PzCH+Xw8jh9U8ntLzf+4mEP+5Gffc7rl2o/9utKOz/F2b9q1WfPsizX9b+c9/eCWKPkkspJV/YUvxcjfptOT9Tmy5HtsQhhDRk/rM3wDe/DKRhuP/HBt2Hbl/ZS9NF+O805ymiIimlXCxnzj61Ob8t6BdpzPF8vxrzz8aYE3hpsYlbXztqeV+N+f/gGv2rMf+/NebZliNXgHibEidfPldrfoegn6k5d0cSmz0PrnA2Zreg/hzseRq6WQ7Bczzkdy/Pby/BP+3y/LagX6Q9b/sRX79E2L/1+vsX1/1jcH/w4tz/hj3XvK/2/Ks9/4I/jHtHJPtfDrkPKXwZIfc7BJ2C6m9nn2HI3R3Z5AJR8hxDyLixny7kfldE/vcF3cP/SND9Lnk/y7A7IY2vK/XPZqWuFyb0FMHtlbrmfbXsXy37F2zZb4ezP0/DflvOf3md/mnt+h8cdv9jjfrtaPZB2P0uE7k36b8VeP8UJv22qNcX6ndJ+1kadM6UfDXon41Bj7+xVNe8rwb9q0H/gg36HRHtz9Oi3yHoZ2rS/+jQ+x9r029HtH+3Tf+t4PunsOm3Rf1CbTobWl9t+mdj07my8LFFuuZ9LjadGwpPTzfLS/7dOb2gwD+cc9H+4ZzVanG63d/N0X+Sqd0Q0X+65tf/BvwKSOA=').then(json => {\n",
       "   const obj = Core.parse(json);\n",
       "   Core.draw('root_plot_1779222895843', 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_1779222895843();\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
}