{
"cells": [
{
"cell_type": "markdown",
"id": "b138f6d9",
"metadata": {},
"source": [
"# rf503_wspaceread\n",
"'ORGANIZATION AND SIMULTANEOUS FITS' RooFit tutorial macro #503\n",
"\n",
"Reading and using a workspace\n",
"\n",
"The input file for self macro is generated by rf502_wspaceread.py\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:32 PM."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "1e7f3925",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:13.263234Z",
"iopub.status.busy": "2026-05-19T20:32:13.263118Z",
"iopub.status.idle": "2026-05-19T20:32:14.216531Z",
"shell.execute_reply": "2026-05-19T20:32:14.216039Z"
}
},
"outputs": [],
"source": [
"import ROOT"
]
},
{
"cell_type": "markdown",
"id": "50b72896",
"metadata": {},
"source": [
"Read workspace from file\n",
"-----------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "d74fb114",
"metadata": {},
"source": [
"Open input file with workspace (generated by rf503_wspacewrite)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "e82f853d",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:14.218753Z",
"iopub.status.busy": "2026-05-19T20:32:14.218595Z",
"iopub.status.idle": "2026-05-19T20:32:14.414417Z",
"shell.execute_reply": "2026-05-19T20:32:14.413724Z"
}
},
"outputs": [],
"source": [
"f = ROOT.TFile(\"rf502_workspace_py.root\")"
]
},
{
"cell_type": "markdown",
"id": "35ce05ad",
"metadata": {},
"source": [
"Retrieve workspace from file"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "4eefa2a9",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:14.416096Z",
"iopub.status.busy": "2026-05-19T20:32:14.415969Z",
"iopub.status.idle": "2026-05-19T20:32:14.610661Z",
"shell.execute_reply": "2026-05-19T20:32:14.610099Z"
}
},
"outputs": [],
"source": [
"w = f.Get(\"w\")"
]
},
{
"cell_type": "markdown",
"id": "8c21f62a",
"metadata": {},
"source": [
"Retrieve pdf, data from workspace\n",
"-----------------------------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "6c1e9928",
"metadata": {},
"source": [
"Retrieve x, and data from workspace"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "c14b491c",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:14.618636Z",
"iopub.status.busy": "2026-05-19T20:32:14.618486Z",
"iopub.status.idle": "2026-05-19T20:32:14.791072Z",
"shell.execute_reply": "2026-05-19T20:32:14.790518Z"
}
},
"outputs": [],
"source": [
"x = w[\"x\"]\n",
"model = w[\"model\"]\n",
"data = w[\"modelData\"]"
]
},
{
"cell_type": "markdown",
"id": "3d6b3cd9",
"metadata": {},
"source": [
"Print structure of composite p.d.f."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "c6a7cfcb",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:14.792884Z",
"iopub.status.busy": "2026-05-19T20:32:14.792756Z",
"iopub.status.idle": "2026-05-19T20:32:14.902728Z",
"shell.execute_reply": "2026-05-19T20:32:14.902038Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0x563ffe70b010 RooAddPdf::model = 1/1 [Auto,Clean] \n",
" 0x563ffe8ae1e0/V- RooChebychev::bkg = 1 [Auto,Dirty] \n",
" 0x563ffe10a1b0/V- RooRealVar::x = 5\n",
" 0x563ff5706f40/V- RooRealVar::a0 = 0.5\n",
" 0x563ff56ff780/V- RooRealVar::a1 = 0\n",
" 0x563ffe9543b0/V- RooRealVar::bkgfrac = 0.5\n",
" 0x563ffe947260/V- RooAddPdf::sig = 1/1 [Auto,Clean] \n",
" 0x563ffe8ce0e0/V- RooGaussian::sig1 = 1 [Auto,Dirty] \n",
" 0x563ffe10a1b0/V- RooRealVar::x = 5\n",
" 0x563ffe0a9870/V- RooRealVar::mean = 5\n",
" 0x563ff5d716a0/V- RooRealVar::sigma1 = 0.5\n",
" 0x563ff5d70490/V- RooRealVar::sig1frac = 0.8\n",
" 0x563ffe8db850/V- RooGaussian::sig2 = 1 [Auto,Dirty] \n",
" 0x563ffe10a1b0/V- RooRealVar::x = 5\n",
" 0x563ffe0a9870/V- RooRealVar::mean = 5\n",
" 0x563ffdf996b0/V- RooRealVar::sigma2 = 1\n"
]
}
],
"source": [
"model.Print(\"t\")"
]
},
{
"cell_type": "markdown",
"id": "ce340b80",
"metadata": {},
"source": [
"Fit model to data, plot model\n",
"---------------------------------------------------------"
]
},
{
"cell_type": "markdown",
"id": "7bcd17a8",
"metadata": {},
"source": [
"Fit model to data"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "f23eb51e",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:14.904383Z",
"iopub.status.busy": "2026-05-19T20:32:14.904252Z",
"iopub.status.idle": "2026-05-19T20:32:15.105981Z",
"shell.execute_reply": "2026-05-19T20:32:15.105421Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Fitting -- RooAbsPdf::fitTo(model) fixing normalization set for coefficient determination to observables in data\n",
"[#1] INFO:Fitting -- using generic CPU library compiled with no vectorizations\n",
"[#1] INFO:Fitting -- Creation of NLL object took 7.09194 ms\n",
"[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model_modelData) Summation contains a RooNLLVar, using its error level\n",
"[#1] INFO:Minimization -- [fitFCN] No discrete parameters, performing continuous minimization only\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"model.fitTo(data, PrintLevel=-1)"
]
},
{
"cell_type": "markdown",
"id": "7124a7ab",
"metadata": {},
"source": [
"Plot data and PDF overlaid"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "33e363bf",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:15.107799Z",
"iopub.status.busy": "2026-05-19T20:32:15.107660Z",
"iopub.status.idle": "2026-05-19T20:32:15.287613Z",
"shell.execute_reply": "2026-05-19T20:32:15.287071Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"xframe = x.frame(Title=\"Model and data read from workspace\")\n",
"data.plotOn(xframe)\n",
"model.plotOn(xframe)"
]
},
{
"cell_type": "markdown",
"id": "1fa74600",
"metadata": {},
"source": [
"Overlay the background component of model with a dashed line"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "c1ef8d4e",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:15.289550Z",
"iopub.status.busy": "2026-05-19T20:32:15.289421Z",
"iopub.status.idle": "2026-05-19T20:32:15.404117Z",
"shell.execute_reply": "2026-05-19T20:32:15.403584Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooAbsPdf::plotOn(model) directly selected PDF components: (bkg)\n",
"[#1] INFO:Plotting -- RooAbsPdf::plotOn(model) indirectly selected PDF components: ()\n"
]
}
],
"source": [
"model.plotOn(xframe, Components=\"bkg\", LineStyle=\"--\")"
]
},
{
"cell_type": "markdown",
"id": "f535d1f3",
"metadata": {},
"source": [
"Overlay the background+sig2 components of model with a dotted line"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "8fd3caf7",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:15.405630Z",
"iopub.status.busy": "2026-05-19T20:32:15.405489Z",
"iopub.status.idle": "2026-05-19T20:32:15.510589Z",
"shell.execute_reply": "2026-05-19T20:32:15.510075Z"
}
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Plotting -- RooAbsPdf::plotOn(model) directly selected PDF components: (bkg,sig2)\n",
"[#1] INFO:Plotting -- RooAbsPdf::plotOn(model) indirectly selected PDF components: (sig)\n"
]
},
{
"data": {
"text/plain": [
""
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"model.plotOn(xframe, Components=\"bkg,sig2\", LineStyle=\":\")"
]
},
{
"cell_type": "markdown",
"id": "56adf8c3",
"metadata": {},
"source": [
"Draw the frame on the canvas"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "67d77935",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:15.512665Z",
"iopub.status.busy": "2026-05-19T20:32:15.512518Z",
"iopub.status.idle": "2026-05-19T20:32:15.726125Z",
"shell.execute_reply": "2026-05-19T20:32:15.725584Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Info in : png file rf503_wspaceread.png has been created\n"
]
}
],
"source": [
"c = ROOT.TCanvas(\"rf503_wspaceread\", \"rf503_wspaceread\", 600, 600)\n",
"ROOT.gPad.SetLeftMargin(0.15)\n",
"xframe.GetYaxis().SetTitleOffset(1.4)\n",
"xframe.Draw()\n",
"\n",
"c.SaveAs(\"rf503_wspaceread.png\")"
]
},
{
"cell_type": "markdown",
"id": "0c3f8cfe",
"metadata": {},
"source": [
"Draw all canvases "
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "d09a08b6",
"metadata": {
"collapsed": false,
"execution": {
"iopub.execute_input": "2026-05-19T20:32:15.727865Z",
"iopub.status.busy": "2026-05-19T20:32:15.727734Z",
"iopub.status.idle": "2026-05-19T20:32:15.908977Z",
"shell.execute_reply": "2026-05-19T20:32:15.908358Z"
}
},
"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
}