{
"cells": [
{
"cell_type": "markdown",
"id": "f3b17c86",
"metadata": {},
"source": [
"# rf310_sliceplot\n",
"Multidimensional models: projecting pdf and data slices in discrete observables\n",
"\n",
"\n",
"\n",
"\n",
"**Author:** Clemens Lange, Wouter Verkerke (C++ version) \n",
"This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, May 19, 2026 at 08:31 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "4ac7f6c9",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:03.923545Z",
"iopub.status.busy": "2026-05-19T20:31:03.923407Z",
"iopub.status.idle": "2026-05-19T20:31:04.876579Z",
"shell.execute_reply": "2026-05-19T20:31:04.875933Z"
}
},
"outputs": [],
"source": [
"import ROOT"
]
},
{
"cell_type": "markdown",
"id": "da195efe",
"metadata": {},
"source": [
"Create B decay pdf with mixing\n",
"----------------------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "4c50d1b9",
"metadata": {},
"source": [
"Decay time observables"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "44eb99c5",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:04.878783Z",
"iopub.status.busy": "2026-05-19T20:31:04.878648Z",
"iopub.status.idle": "2026-05-19T20:31:05.042545Z",
"shell.execute_reply": "2026-05-19T20:31:05.041830Z"
}
},
"outputs": [],
"source": [
"dt = ROOT.RooRealVar(\"dt\", \"dt\", -20, 20)"
]
},
{
"cell_type": "markdown",
"id": "a05e47ea",
"metadata": {},
"source": [
"Discrete observables mixState (B0tag==B0reco?) and tagFlav\n",
"(B0tag==B0(bar)?)\n",
"Define state labels of discrete observables"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "5c2a8574",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.044137Z",
"iopub.status.busy": "2026-05-19T20:31:05.044016Z",
"iopub.status.idle": "2026-05-19T20:31:05.160758Z",
"shell.execute_reply": "2026-05-19T20:31:05.160065Z"
}
},
"outputs": [],
"source": [
"mixState = ROOT.RooCategory(\"mixState\", \"B0/B0bar mixing state\", {\"mixed\": -1, \"unmixed\": 1})\n",
"tagFlav = ROOT.RooCategory(\"tagFlav\", \"Flavour of the tagged B0\", {\"B0\": 1, \"B0bar\": -1})"
]
},
{
"cell_type": "markdown",
"id": "29b2b81f",
"metadata": {},
"source": [
"Model parameters"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "d4ef5185",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.162987Z",
"iopub.status.busy": "2026-05-19T20:31:05.162864Z",
"iopub.status.idle": "2026-05-19T20:31:05.308859Z",
"shell.execute_reply": "2026-05-19T20:31:05.289379Z"
}
},
"outputs": [],
"source": [
"dm = ROOT.RooRealVar(\"dm\", \"delta m(B)\", 0.472, 0.0, 1.0)\n",
"tau = ROOT.RooRealVar(\"tau\", \"B0 decay time\", 1.547, 1.0, 2.0)\n",
"w = ROOT.RooRealVar(\"w\", \"Flavor Mistag rate\", 0.03, 0.0, 1.0)\n",
"dw = ROOT.RooRealVar(\"dw\", \"Flavor Mistag rate difference between B0 and B0bar\", 0.01)"
]
},
{
"cell_type": "markdown",
"id": "ad2f0573",
"metadata": {},
"source": [
"Build a gaussian resolution model"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "2ecdbb06",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.320397Z",
"iopub.status.busy": "2026-05-19T20:31:05.320255Z",
"iopub.status.idle": "2026-05-19T20:31:05.458624Z",
"shell.execute_reply": "2026-05-19T20:31:05.458076Z"
}
},
"outputs": [],
"source": [
"bias1 = ROOT.RooRealVar(\"bias1\", \"bias1\", 0)\n",
"sigma1 = ROOT.RooRealVar(\"sigma1\", \"sigma1\", 0.01)\n",
"gm1 = ROOT.RooGaussModel(\"gm1\", \"gauss model 1\", dt, bias1, sigma1)"
]
},
{
"cell_type": "markdown",
"id": "a4cb2a3f",
"metadata": {},
"source": [
"Construct a decay pdf, with single gaussian resolution model"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "88a254ba",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.460409Z",
"iopub.status.busy": "2026-05-19T20:31:05.460284Z",
"iopub.status.idle": "2026-05-19T20:31:05.595538Z",
"shell.execute_reply": "2026-05-19T20:31:05.595014Z"
}
},
"outputs": [],
"source": [
"bmix_gm1 = ROOT.RooBMixDecay(\"bmix\", \"decay\", dt, mixState, tagFlav, tau, dm, w, dw, gm1, type=\"DoubleSided\")"
]
},
{
"cell_type": "markdown",
"id": "608416f0",
"metadata": {},
"source": [
"Generate BMixing data with above set of event errors"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "6c125878",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.597168Z",
"iopub.status.busy": "2026-05-19T20:31:05.597046Z",
"iopub.status.idle": "2026-05-19T20:31:05.770214Z",
"shell.execute_reply": "2026-05-19T20:31:05.769713Z"
}
},
"outputs": [],
"source": [
"data = bmix_gm1.generate({dt, tagFlav, mixState}, 20000)"
]
},
{
"cell_type": "markdown",
"id": "018bcc6c",
"metadata": {},
"source": [
"Plot full decay distribution\n",
"----------------------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "e5bf2cdd",
"metadata": {},
"source": [
"Create frame, data and pdf projection (integrated over tagFlav and\n",
"mixState)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "288746e3",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.772286Z",
"iopub.status.busy": "2026-05-19T20:31:05.772162Z",
"iopub.status.idle": "2026-05-19T20:31:05.958246Z",
"shell.execute_reply": "2026-05-19T20:31:05.957588Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt integrates over variables (tagFlav,mixState)\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"frame = dt.frame(Title=\"Inclusive decay distribution\")\n",
"data.plotOn(frame)\n",
"bmix_gm1.plotOn(frame)"
]
},
{
"cell_type": "markdown",
"id": "ee1bbd59",
"metadata": {},
"source": [
"Plot decay distribution for mixed and unmixed slice of mixState\n",
"-------------------------------------------------------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "0b287a29",
"metadata": {},
"source": [
"Create frame, data (mixed only)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "7ae4e80a",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:05.970372Z",
"iopub.status.busy": "2026-05-19T20:31:05.970238Z",
"iopub.status.idle": "2026-05-19T20:31:06.088546Z",
"shell.execute_reply": "2026-05-19T20:31:06.088019Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 3787 events out of 20000 total events\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"frame2 = dt.frame(Title=\"Decay distribution of mixed events\")\n",
"data.plotOn(frame2, Cut=\"mixState==mixState::mixed\")"
]
},
{
"cell_type": "markdown",
"id": "bbbf5d1d",
"metadata": {},
"source": [
"Position slice in mixState at \"mixed\" and plot slice of pdf in mixstate\n",
"over data (integrated over tagFlav)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "f2be5300",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:06.090233Z",
"iopub.status.busy": "2026-05-19T20:31:06.090108Z",
"iopub.status.idle": "2026-05-19T20:31:06.221123Z",
"shell.execute_reply": "2026-05-19T20:31:06.220247Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState)\n",
"[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt integrates over variables (tagFlav)\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"bmix_gm1.plotOn(frame2, Slice=(mixState, \"mixed\"))"
]
},
{
"cell_type": "markdown",
"id": "bf2a01ae",
"metadata": {},
"source": [
"Create frame, data (unmixed only)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "560afd04",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:06.222757Z",
"iopub.status.busy": "2026-05-19T20:31:06.222623Z",
"iopub.status.idle": "2026-05-19T20:31:06.337982Z",
"shell.execute_reply": "2026-05-19T20:31:06.337431Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 16213 events out of 20000 total events\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"frame3 = dt.frame(Title=\"Decay distribution of unmixed events\")\n",
"data.plotOn(frame3, Cut=\"mixState==mixState::unmixed\")"
]
},
{
"cell_type": "markdown",
"id": "a2692622",
"metadata": {},
"source": [
"Position slice in mixState at \"unmixed\" and plot slice of pdf in\n",
"mixstate over data (integrated over tagFlav)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "842a79d4",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:06.339734Z",
"iopub.status.busy": "2026-05-19T20:31:06.339590Z",
"iopub.status.idle": "2026-05-19T20:31:06.581519Z",
"shell.execute_reply": "2026-05-19T20:31:06.580991Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState)\n",
"[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt integrates over variables (tagFlav)\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"Info in : png file rf310_sliceplot.png has been created\n"
]
}
],
"source": [
"bmix_gm1.plotOn(frame3, Slice=(mixState, \"unmixed\"))\n",
"\n",
"c = ROOT.TCanvas(\"rf310_sliceplot\", \"rf310_sliceplot\", 1200, 400)\n",
"c.Divide(3)\n",
"c.cd(1)\n",
"ROOT.gPad.SetLeftMargin(0.15)\n",
"frame.GetYaxis().SetTitleOffset(1.4)\n",
"ROOT.gPad.SetLogy()\n",
"frame.Draw()\n",
"c.cd(2)\n",
"ROOT.gPad.SetLeftMargin(0.15)\n",
"frame2.GetYaxis().SetTitleOffset(1.4)\n",
"ROOT.gPad.SetLogy()\n",
"frame2.Draw()\n",
"c.cd(3)\n",
"ROOT.gPad.SetLeftMargin(0.15)\n",
"frame3.GetYaxis().SetTitleOffset(1.4)\n",
"ROOT.gPad.SetLogy()\n",
"frame3.Draw()\n",
"\n",
"c.SaveAs(\"rf310_sliceplot.png\")"
]
},
{
"cell_type": "markdown",
"id": "d0f530f8",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "2009161d",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:31:06.583274Z",
"iopub.status.busy": "2026-05-19T20:31:06.583151Z",
"iopub.status.idle": "2026-05-19T20:31:06.765059Z",
"shell.execute_reply": "2026-05-19T20:31:06.764500Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"
\n",
"\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"from ROOT import gROOT \n",
"gROOT.GetListOfCanvases().Draw()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.12"
}
},
"nbformat": 4,
"nbformat_minor": 5
}