{ "cells": [ { "cell_type": "markdown", "id": "08677d85", "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:** 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:34 PM." ] }, { "cell_type": "code", "execution_count": 1, "id": "ef8a28c6", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:53.094515Z", "iopub.status.busy": "2026-05-19T20:34:53.094391Z", "iopub.status.idle": "2026-05-19T20:34:54.049962Z", "shell.execute_reply": "2026-05-19T20:34:54.049319Z" } }, "outputs": [], "source": [ "import ROOT" ] }, { "cell_type": "markdown", "id": "943f712b", "metadata": {}, "source": [ "B-decay with mixing\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "a0b778d8", "metadata": {}, "source": [ "Construct pdf\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "9677d8de", "metadata": {}, "source": [ "Observable" ] }, { "cell_type": "code", "execution_count": 2, "id": "1ad6e629", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.051649Z", "iopub.status.busy": "2026-05-19T20:34:54.051504Z", "iopub.status.idle": "2026-05-19T20:34:54.212455Z", "shell.execute_reply": "2026-05-19T20:34:54.211877Z" } }, "outputs": [], "source": [ "dt = ROOT.RooRealVar(\"dt\", \"dt\", -10, 10)\n", "dt.setBins(40)" ] }, { "cell_type": "markdown", "id": "684900a1", "metadata": {}, "source": [ "Parameters" ] }, { "cell_type": "code", "execution_count": 3, "id": "debfe7ee", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.214040Z", "iopub.status.busy": "2026-05-19T20:34:54.213921Z", "iopub.status.idle": "2026-05-19T20:34:54.334716Z", "shell.execute_reply": "2026-05-19T20:34:54.334099Z" } }, "outputs": [], "source": [ "dm = ROOT.RooRealVar(\"dm\", \"delta m(B0)\", 0.472)\n", "tau = ROOT.RooRealVar(\"tau\", \"tau (B0)\", 1.547)\n", "w = ROOT.RooRealVar(\"w\", \"flavour mistag rate\", 0.1)\n", "dw = ROOT.RooRealVar(\"dw\", \"delta mistag rate for B0/B0bar\", 0.1)\n", "\n", "mixState = ROOT.RooCategory(\"mixState\", \"B0/B0bar mixing state\", {\"mixed\": -1, \"unmixed\": 1})\n", "\n", "tagFlav = ROOT.RooCategory(\"tagFlav\", \"Flavour of the tagged B0\", {\"B0\": 1, \"B0bar\": -1})" ] }, { "cell_type": "markdown", "id": "6f2f5137", "metadata": {}, "source": [ "Use delta function resolution model" ] }, { "cell_type": "code", "execution_count": 4, "id": "d3a941f8", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.336289Z", "iopub.status.busy": "2026-05-19T20:34:54.336170Z", "iopub.status.idle": "2026-05-19T20:34:54.451911Z", "shell.execute_reply": "2026-05-19T20:34:54.451313Z" } }, "outputs": [], "source": [ "tm = ROOT.RooTruthModel(\"tm\", \"truth model\", dt)" ] }, { "cell_type": "markdown", "id": "2476f96b", "metadata": {}, "source": [ "Construct Bdecay with mixing" ] }, { "cell_type": "code", "execution_count": 5, "id": "974db23e", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.453463Z", "iopub.status.busy": "2026-05-19T20:34:54.453344Z", "iopub.status.idle": "2026-05-19T20:34:54.589897Z", "shell.execute_reply": "2026-05-19T20:34:54.589284Z" } }, "outputs": [], "source": [ "bmix = ROOT.RooBMixDecay(\"bmix\", \"decay\", dt, mixState, tagFlav, tau, dm, w, dw, tm, type=\"DoubleSided\")" ] }, { "cell_type": "markdown", "id": "5738bdb6", "metadata": {}, "source": [ "Plot pdf in various slices\n", "---------------------------------------------------" ] }, { "cell_type": "markdown", "id": "f3f8111d", "metadata": {}, "source": [ "Generate some data" ] }, { "cell_type": "code", "execution_count": 6, "id": "e165c9ae", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.591492Z", "iopub.status.busy": "2026-05-19T20:34:54.591371Z", "iopub.status.idle": "2026-05-19T20:34:54.744793Z", "shell.execute_reply": "2026-05-19T20:34:54.744186Z" } }, "outputs": [], "source": [ "data = bmix.generate({dt, mixState, tagFlav}, 10000)" ] }, { "cell_type": "markdown", "id": "b31e6b5a", "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": "f35540a7", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.746364Z", "iopub.status.busy": "2026-05-19T20:34:54.746241Z", "iopub.status.idle": "2026-05-19T20:34:54.979920Z", "shell.execute_reply": "2026-05-19T20:34:54.979323Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 4542 events out of 10000 total events\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (tagFlav)\n", "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt integrates over variables (mixState)\n", "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 5458 events out of 10000 total events\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (tagFlav)\n", "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt integrates over variables (mixState)\n" ] } ], "source": [ "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=\"c\")\n", "bmix.plotOn(frame1, Slice=(tagFlav, \"B0bar\"), LineColor=\"c\")" ] }, { "cell_type": "markdown", "id": "257a9e4d", "metadata": {}, "source": [ "Plot mixed slice for B0 and B0bar tagged data separately" ] }, { "cell_type": "code", "execution_count": 8, "id": "98a71f19", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:54.981265Z", "iopub.status.busy": "2026-05-19T20:34:54.981143Z", "iopub.status.idle": "2026-05-19T20:34:55.193656Z", "shell.execute_reply": "2026-05-19T20:34:55.193060Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 933 events out of 10000 total events\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState,tagFlav)\n", "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 1461 events out of 10000 total events\n", "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState,tagFlav)\n" ] } ], "source": [ "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=\"c\")\n", "bmix.plotOn(frame2, Slice={tagFlav: \"B0bar\", mixState: \"mixed\"}, LineColor=\"c\")" ] }, { "cell_type": "markdown", "id": "ea09248b", "metadata": {}, "source": [ "Plot unmixed slice for B0 and B0bar tagged data separately" ] }, { "cell_type": "code", "execution_count": 9, "id": "75f68838", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.194976Z", "iopub.status.busy": "2026-05-19T20:34:55.194858Z", "iopub.status.idle": "2026-05-19T20:34:55.318826Z", "shell.execute_reply": "2026-05-19T20:34:55.318198Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 3609 events out of 10000 total events\n", "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState,tagFlav)\n", "[#1] INFO:Plotting -- RooTreeData::plotOn: plotting 3997 events out of 10000 total events\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsReal::plotOn(bmix) plot on dt represents a slice in (mixState,tagFlav)\n" ] } ], "source": [ "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=\"c\")\n", "bmix.plotOn(frame3, Slice={tagFlav: \"B0bar\", mixState: \"unmixed\"}, LineColor=\"c\")" ] }, { "cell_type": "markdown", "id": "ae3a658c", "metadata": {}, "source": [ "B-decay with CP violation\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "a177005c", "metadata": {}, "source": [ "Construct pdf\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "d899ea7c", "metadata": {}, "source": [ "Additional parameters needed for B decay with CPV" ] }, { "cell_type": "code", "execution_count": 10, "id": "6bfd64a9", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.320204Z", "iopub.status.busy": "2026-05-19T20:34:55.320085Z", "iopub.status.idle": "2026-05-19T20:34:55.430541Z", "shell.execute_reply": "2026-05-19T20:34:55.429941Z" } }, "outputs": [], "source": [ "CPeigen = ROOT.RooRealVar(\"CPeigen\", \"CP eigen value\", -1)\n", "absLambda = ROOT.RooRealVar(\"absLambda\", \"|lambda|\", 1, 0, 2)\n", "argLambda = ROOT.RooRealVar(\"absLambda\", \"|lambda|\", 0.7, -1, 1)\n", "effR = ROOT.RooRealVar(\"effR\", \"B0/B0bar reco efficiency ratio\", 1)" ] }, { "cell_type": "markdown", "id": "2ab1a675", "metadata": {}, "source": [ "Construct Bdecay with CP violation" ] }, { "cell_type": "code", "execution_count": 11, "id": "f2df63bb", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.432141Z", "iopub.status.busy": "2026-05-19T20:34:55.432017Z", "iopub.status.idle": "2026-05-19T20:34:55.555121Z", "shell.execute_reply": "2026-05-19T20:34:55.554493Z" } }, "outputs": [], "source": [ "bcp = ROOT.RooBCPEffDecay(\n", " \"bcp\", \"bcp\", dt, tagFlav, tau, dm, w, CPeigen, absLambda, argLambda, effR, dw, tm, type=\"DoubleSided\"\n", ")" ] }, { "cell_type": "markdown", "id": "d5eebcb3", "metadata": {}, "source": [ "Plot scenario 1 - sin(2b)=0.7, |l|=1\n", "---------------------------------------------------------------------------" ] }, { "cell_type": "markdown", "id": "3391f297", "metadata": {}, "source": [ "Generate some data" ] }, { "cell_type": "code", "execution_count": 12, "id": "bc9445e5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.556720Z", "iopub.status.busy": "2026-05-19T20:34:55.556576Z", "iopub.status.idle": "2026-05-19T20:34:55.671660Z", "shell.execute_reply": "2026-05-19T20:34:55.671088Z" } }, "outputs": [], "source": [ "data2 = bcp.generate({dt, tagFlav}, 10000)" ] }, { "cell_type": "markdown", "id": "f33cf62e", "metadata": {}, "source": [ "Plot B0 and B0bar tagged data separately" ] }, { "cell_type": "code", "execution_count": 13, "id": "8e54e5ec", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.682529Z", "iopub.status.busy": "2026-05-19T20:34:55.682400Z", "iopub.status.idle": "2026-05-19T20:34:55.790849Z", "shell.execute_reply": "2026-05-19T20:34:55.790167Z" } }, "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" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "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=\"c\")\n", "bcp.plotOn(frame4, Slice=(tagFlav, \"B0bar\"), LineColor=\"c\")" ] }, { "cell_type": "markdown", "id": "776b7b38", "metadata": {}, "source": [ "# Plot scenario 2 - sin(2b)=0.7, |l|=0.7\n", "-------------------------------------------------------------------------------" ] }, { "cell_type": "code", "execution_count": 14, "id": "e6a1c96b", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.792359Z", "iopub.status.busy": "2026-05-19T20:34:55.792236Z", "iopub.status.idle": "2026-05-19T20:34:55.900298Z", "shell.execute_reply": "2026-05-19T20:34:55.899632Z" } }, "outputs": [], "source": [ "absLambda.setVal(0.7)" ] }, { "cell_type": "markdown", "id": "8d923611", "metadata": {}, "source": [ "Generate some data" ] }, { "cell_type": "code", "execution_count": 15, "id": "d14926f5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:55.901785Z", "iopub.status.busy": "2026-05-19T20:34:55.901658Z", "iopub.status.idle": "2026-05-19T20:34:56.009226Z", "shell.execute_reply": "2026-05-19T20:34:56.008514Z" } }, "outputs": [], "source": [ "data3 = bcp.generate({dt, tagFlav}, 10000)" ] }, { "cell_type": "markdown", "id": "88d377b5", "metadata": {}, "source": [ "Plot B0 and B0bar tagged data separately (sin2b = 0.7 plus direct CPV\n", "|l|=0.5)" ] }, { "cell_type": "code", "execution_count": 16, "id": "19eb3cf3", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.010629Z", "iopub.status.busy": "2026-05-19T20:34:56.010497Z", "iopub.status.idle": "2026-05-19T20:34:56.118894Z", "shell.execute_reply": "2026-05-19T20:34:56.118215Z" } }, "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" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "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=\"c\")\n", "bcp.plotOn(frame5, Slice=(tagFlav, \"B0bar\"), LineColor=\"c\")" ] }, { "cell_type": "markdown", "id": "e2a48127", "metadata": {}, "source": [ "Generic B-decay with user coefficients\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "3f05c49c", "metadata": {}, "source": [ "Construct pdf\n", "-------------------------" ] }, { "cell_type": "markdown", "id": "1016f5f6", "metadata": {}, "source": [ "Model parameters" ] }, { "cell_type": "code", "execution_count": 17, "id": "ea86a76b", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.120382Z", "iopub.status.busy": "2026-05-19T20:34:56.120264Z", "iopub.status.idle": "2026-05-19T20:34:56.223917Z", "shell.execute_reply": "2026-05-19T20:34:56.223287Z" } }, "outputs": [], "source": [ "DGbG = ROOT.RooRealVar(\"DGbG\", \"DGamma/GammaAvg\", 0.5, -1, 1)\n", "Adir = ROOT.RooRealVar(\"Adir\", \"-[1-abs(l)**2]/[1+abs(l)**2]\", 0)\n", "Amix = ROOT.RooRealVar(\"Amix\", \"2Im(l)/[1+abs(l)**2]\", 0.7)\n", "Adel = ROOT.RooRealVar(\"Adel\", \"2Re(l)/[1+abs(l)**2]\", 0.7)" ] }, { "cell_type": "markdown", "id": "58ac7e81", "metadata": {}, "source": [ "Derived input parameters for pdf" ] }, { "cell_type": "code", "execution_count": 18, "id": "edde823a", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.228079Z", "iopub.status.busy": "2026-05-19T20:34:56.227957Z", "iopub.status.idle": "2026-05-19T20:34:56.400262Z", "shell.execute_reply": "2026-05-19T20:34:56.399636Z" } }, "outputs": [], "source": [ "DG = ROOT.RooFormulaVar(\"DG\", \"Delta Gamma\", \"@1/@0\", [tau, DGbG])" ] }, { "cell_type": "markdown", "id": "56a4c322", "metadata": {}, "source": [ "Construct coefficient functions for sin,cos, modulations of decay\n", "distribution" ] }, { "cell_type": "code", "execution_count": 19, "id": "441095c6", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.401764Z", "iopub.status.busy": "2026-05-19T20:34:56.401639Z", "iopub.status.idle": "2026-05-19T20:34:56.519993Z", "shell.execute_reply": "2026-05-19T20:34:56.519374Z" } }, "outputs": [], "source": [ "fsin = ROOT.RooFormulaVar(\"fsin\", \"fsin\", \"@0*@1*(1-2*@2)\", [Amix, tagFlav, w])\n", "fcos = ROOT.RooFormulaVar(\"fcos\", \"fcos\", \"@0*@1*(1-2*@2)\", [Adir, tagFlav, w])\n", "fsinh = ROOT.RooFormulaVar(\"fsinh\", \"fsinh\", \"@0\", [Adel])" ] }, { "cell_type": "markdown", "id": "bc4a499b", "metadata": {}, "source": [ "Construct generic B decay pdf using above user coefficients" ] }, { "cell_type": "code", "execution_count": 20, "id": "17b33a01", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.521457Z", "iopub.status.busy": "2026-05-19T20:34:56.521339Z", "iopub.status.idle": "2026-05-19T20:34:56.653284Z", "shell.execute_reply": "2026-05-19T20:34:56.652641Z" } }, "outputs": [], "source": [ "bcpg = ROOT.RooBDecay(\n", " \"bcpg\", \"bcpg\", dt, tau, DG, ROOT.RooFit.RooConst(1), fsinh, fcos, fsin, dm, tm, type=\"DoubleSided\"\n", ")" ] }, { "cell_type": "markdown", "id": "32ca67ad", "metadata": {}, "source": [ "Plot - Im(l)=0.7, e(l)=0.7 |l|=1, G/G=0.5\n", "-------------------------------------------------------------------------------------" ] }, { "cell_type": "markdown", "id": "9b9a9784", "metadata": {}, "source": [ "Generate some data" ] }, { "cell_type": "code", "execution_count": 21, "id": "c69441c5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.654746Z", "iopub.status.busy": "2026-05-19T20:34:56.654620Z", "iopub.status.idle": "2026-05-19T20:34:56.839992Z", "shell.execute_reply": "2026-05-19T20:34:56.839257Z" } }, "outputs": [], "source": [ "data4 = bcpg.generate({dt, tagFlav}, 10000)" ] }, { "cell_type": "markdown", "id": "d579f7ef", "metadata": {}, "source": [ "Plot B0 and B0bar tagged data separately" ] }, { "cell_type": "code", "execution_count": 22, "id": "f19da528", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.842112Z", "iopub.status.busy": "2026-05-19T20:34:56.841981Z", "iopub.status.idle": "2026-05-19T20:34:57.240729Z", "shell.execute_reply": "2026-05-19T20:34:57.240015Z" } }, "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": [ "Info in : png file rf708_bphysics.png has been created\n" ] } ], "source": [ "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=\"c\")\n", "bcpg.plotOn(frame6, Slice=(tagFlav, \"B0bar\"), LineColor=\"c\")\n", "\n", "c = ROOT.TCanvas(\"rf708_bphysics\", \"rf708_bphysics\", 1200, 800)\n", "c.Divide(3, 2)\n", "c.cd(1)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame1.GetYaxis().SetTitleOffset(1.6)\n", "frame1.Draw()\n", "c.cd(2)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame2.GetYaxis().SetTitleOffset(1.6)\n", "frame2.Draw()\n", "c.cd(3)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame3.GetYaxis().SetTitleOffset(1.6)\n", "frame3.Draw()\n", "c.cd(4)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame4.GetYaxis().SetTitleOffset(1.6)\n", "frame4.Draw()\n", "c.cd(5)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame5.GetYaxis().SetTitleOffset(1.6)\n", "frame5.Draw()\n", "c.cd(6)\n", "ROOT.gPad.SetLeftMargin(0.15)\n", "frame6.GetYaxis().SetTitleOffset(1.6)\n", "frame6.Draw()\n", "\n", "c.SaveAs(\"rf708_bphysics.png\")" ] }, { "cell_type": "markdown", "id": "920c1204", "metadata": {}, "source": [ "Draw all canvases " ] }, { "cell_type": "code", "execution_count": 23, "id": "18c02f09", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:57.242872Z", "iopub.status.busy": "2026-05-19T20:34:57.242728Z", "iopub.status.idle": "2026-05-19T20:34:57.436313Z", "shell.execute_reply": "2026-05-19T20:34:57.435747Z" } }, "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 }