{ "cells": [ { "cell_type": "markdown", "id": "8a5192fc", "metadata": {}, "source": [ "# rf709_BarlowBeeston\n", "Implementing the Barlow-Beeston method for taking into account the statistical\n", "uncertainty of a Monte-Carlo fit template.\n", "\n", "\n", "Based on a demo by Wouter Verkerke\n", "\n", "\n", "**Author:** Harshal Shende, Stephan Hageboeck (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": "92701bc5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:56.908076Z", "iopub.status.busy": "2026-05-19T20:34:56.907954Z", "iopub.status.idle": "2026-05-19T20:34:57.864755Z", "shell.execute_reply": "2026-05-19T20:34:57.864142Z" } }, "outputs": [], "source": [ "import ROOT" ] }, { "cell_type": "markdown", "id": "5a515edb", "metadata": {}, "source": [ "First, construct a likelihood model with a Gaussian signal on top of a uniform background" ] }, { "cell_type": "code", "execution_count": 2, "id": "9a663a52", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:57.872645Z", "iopub.status.busy": "2026-05-19T20:34:57.872481Z", "iopub.status.idle": "2026-05-19T20:34:58.074272Z", "shell.execute_reply": "2026-05-19T20:34:58.073638Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#0] WARNING:InputArguments -- The parameter 'sigG' with range [-10, 10] of the RooGaussian 'g' exceeds the safe range of (0, inf). Advise to limit its range.\n" ] } ], "source": [ "x = ROOT.RooRealVar(\"x\", \"x\", -20, 20)\n", "x.setBins(25)\n", "\n", "meanG = ROOT.RooRealVar(\"meanG\", \"meanG\", 1, -10, 10)\n", "sigG = ROOT.RooRealVar(\"sigG\", \"sigG\", 1.5, -10, 10)\n", "g = ROOT.RooGaussian(\"g\", \"Gauss\", x, meanG, sigG)\n", "u = ROOT.RooUniform(\"u\", \"Uniform\", x)" ] }, { "cell_type": "markdown", "id": "8ceec8e8", "metadata": {}, "source": [ "Generate the data to be fitted" ] }, { "cell_type": "code", "execution_count": 3, "id": "3950ab6c", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.076221Z", "iopub.status.busy": "2026-05-19T20:34:58.076083Z", "iopub.status.idle": "2026-05-19T20:34:58.265363Z", "shell.execute_reply": "2026-05-19T20:34:58.264780Z" } }, "outputs": [], "source": [ "sigData = g.generate(x, 50)\n", "bkgData = u.generate(x, 1000)\n", "\n", "sumData = ROOT.RooDataSet(\"sumData\", \"Gauss + Uniform\", x)\n", "sumData.append(sigData)\n", "sumData.append(bkgData)" ] }, { "cell_type": "markdown", "id": "87540901", "metadata": {}, "source": [ "Make histogram templates for signal and background.\n", "Let's take a signal distribution with low statistics and a more accurate\n", "background distribution.\n", "Normally, these come from Monte Carlo simulations, but we will just generate them." ] }, { "cell_type": "code", "execution_count": 4, "id": "1e7a2b4f", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.267423Z", "iopub.status.busy": "2026-05-19T20:34:58.267292Z", "iopub.status.idle": "2026-05-19T20:34:58.391809Z", "shell.execute_reply": "2026-05-19T20:34:58.391134Z" } }, "outputs": [], "source": [ "dh_sig = g.generateBinned(x, 50)\n", "dh_bkg = u.generateBinned(x, 10000)" ] }, { "cell_type": "markdown", "id": "ec647f80", "metadata": {}, "source": [ " Case 0 - 'Rigid templates'" ] }, { "cell_type": "markdown", "id": "dcb0d126", "metadata": {}, "source": [ "Construct histogram shapes for signal and background" ] }, { "cell_type": "code", "execution_count": 5, "id": "975a09ba", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.393679Z", "iopub.status.busy": "2026-05-19T20:34:58.393530Z", "iopub.status.idle": "2026-05-19T20:34:58.508152Z", "shell.execute_reply": "2026-05-19T20:34:58.507628Z" } }, "outputs": [], "source": [ "p_h_sig = ROOT.RooHistFunc(\"p_h_sig\", \"p_h_sig\", x, dh_sig)\n", "p_h_bkg = ROOT.RooHistFunc(\"p_h_bkg\", \"p_h_bkg\", x, dh_bkg)" ] }, { "cell_type": "markdown", "id": "a6d73129", "metadata": {}, "source": [ "Construct scale factors for adding the two distributions" ] }, { "cell_type": "code", "execution_count": 6, "id": "ca323e36", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.510284Z", "iopub.status.busy": "2026-05-19T20:34:58.510155Z", "iopub.status.idle": "2026-05-19T20:34:58.613742Z", "shell.execute_reply": "2026-05-19T20:34:58.613096Z" } }, "outputs": [], "source": [ "Asig0 = ROOT.RooRealVar(\"Asig\", \"Asig\", 1, 0.01, 5000)\n", "Abkg0 = ROOT.RooRealVar(\"Abkg\", \"Abkg\", 1, 0.01, 5000)" ] }, { "cell_type": "markdown", "id": "08866bb7", "metadata": {}, "source": [ "Construct the sum model" ] }, { "cell_type": "code", "execution_count": 7, "id": "bd7a6a29", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.615620Z", "iopub.status.busy": "2026-05-19T20:34:58.615483Z", "iopub.status.idle": "2026-05-19T20:34:58.808764Z", "shell.execute_reply": "2026-05-19T20:34:58.808056Z" } }, "outputs": [], "source": [ "model0 = ROOT.RooRealSumPdf(\"model0\", \"model0\", [p_h_sig, p_h_bkg], [Asig0, Abkg0], True)" ] }, { "cell_type": "markdown", "id": "d9603dca", "metadata": {}, "source": [ " Case 1 - 'Barlow Beeston'" ] }, { "cell_type": "markdown", "id": "e34bf8b5", "metadata": {}, "source": [ "Construct parameterized histogram shapes for signal and background" ] }, { "cell_type": "code", "execution_count": 8, "id": "bb15ff1b", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.810592Z", "iopub.status.busy": "2026-05-19T20:34:58.810460Z", "iopub.status.idle": "2026-05-19T20:34:58.923765Z", "shell.execute_reply": "2026-05-19T20:34:58.923110Z" } }, "outputs": [], "source": [ "p_ph_sig1 = ROOT.RooParamHistFunc(\"p_ph_sig\", \"p_ph_sig\", dh_sig, x)\n", "p_ph_bkg1 = ROOT.RooParamHistFunc(\"p_ph_bkg\", \"p_ph_bkg\", dh_bkg, x)\n", "\n", "Asig1 = ROOT.RooRealVar(\"Asig\", \"Asig\", 1, 0.01, 5000)\n", "Abkg1 = ROOT.RooRealVar(\"Abkg\", \"Abkg\", 1, 0.01, 5000)" ] }, { "cell_type": "markdown", "id": "5c74a898", "metadata": {}, "source": [ "Construct the sum of these" ] }, { "cell_type": "code", "execution_count": 9, "id": "4ac1c1fc", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:58.926025Z", "iopub.status.busy": "2026-05-19T20:34:58.925891Z", "iopub.status.idle": "2026-05-19T20:34:59.029665Z", "shell.execute_reply": "2026-05-19T20:34:59.028976Z" } }, "outputs": [], "source": [ "model_tmp = ROOT.RooRealSumPdf(\"sp_ph\", \"sp_ph\", [p_ph_sig1, p_ph_bkg1], [Asig1, Abkg1], True)" ] }, { "cell_type": "markdown", "id": "5cf3e35c", "metadata": {}, "source": [ "Construct the subsidiary poisson measurements constraining the histogram parameters\n", "These ensure that the bin contents of the histograms are only allowed to vary within\n", "the statistical uncertainty of the Monte Carlo." ] }, { "cell_type": "code", "execution_count": 10, "id": "b844d4ce", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.031938Z", "iopub.status.busy": "2026-05-19T20:34:59.031804Z", "iopub.status.idle": "2026-05-19T20:34:59.144268Z", "shell.execute_reply": "2026-05-19T20:34:59.143535Z" } }, "outputs": [], "source": [ "hc_sig = ROOT.RooHistConstraint(\"hc_sig\", \"hc_sig\", p_ph_sig1)\n", "hc_bkg = ROOT.RooHistConstraint(\"hc_bkg\", \"hc_bkg\", p_ph_bkg1)" ] }, { "cell_type": "markdown", "id": "61eabac5", "metadata": {}, "source": [ "Construct the joint model with template PDFs and constraints" ] }, { "cell_type": "code", "execution_count": 11, "id": "f95f1530", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.146144Z", "iopub.status.busy": "2026-05-19T20:34:59.146009Z", "iopub.status.idle": "2026-05-19T20:34:59.292381Z", "shell.execute_reply": "2026-05-19T20:34:59.291728Z" } }, "outputs": [], "source": [ "model1 = ROOT.RooProdPdf(\"model1\", \"model1\", {hc_sig, hc_bkg}, Conditional=(model_tmp, x))" ] }, { "cell_type": "markdown", "id": "3f2d97d8", "metadata": {}, "source": [ " Case 2 - 'Barlow Beeston' light (one parameter per bin for all samples)" ] }, { "cell_type": "markdown", "id": "6c708a66", "metadata": {}, "source": [ "Construct the histogram shapes, using the same parameters for signal and background\n", "This requires passing the first histogram to the second, so that their common parameters\n", "can be re-used.\n", "The first ParamHistFunc will create one parameter per bin, such as `p_ph_sig2_gamma_bin_0`.\n", "This allows bin 0 to fluctuate up and down.\n", "Then, the SAME parameters are connected to the background histogram, so the bins fluctuate\n", "synchronously. This reduces the number of parameters." ] }, { "cell_type": "code", "execution_count": 12, "id": "28b186b4", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.294696Z", "iopub.status.busy": "2026-05-19T20:34:59.294550Z", "iopub.status.idle": "2026-05-19T20:34:59.398229Z", "shell.execute_reply": "2026-05-19T20:34:59.397807Z" } }, "outputs": [], "source": [ "p_ph_sig2 = ROOT.RooParamHistFunc(\"p_ph_sig2\", \"p_ph_sig2\", dh_sig, x)\n", "p_ph_bkg2 = ROOT.RooParamHistFunc(\"p_ph_bkg2\", \"p_ph_bkg2\", dh_bkg, x, p_ph_sig2, True)\n", "\n", "Asig2 = ROOT.RooRealVar(\"Asig\", \"Asig\", 1, 0.01, 5000)\n", "Abkg2 = ROOT.RooRealVar(\"Abkg\", \"Abkg\", 1, 0.01, 5000)" ] }, { "cell_type": "markdown", "id": "626a0950", "metadata": {}, "source": [ "As before, construct the sum of signal2 and background2" ] }, { "cell_type": "code", "execution_count": 13, "id": "fd85402a", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.416573Z", "iopub.status.busy": "2026-05-19T20:34:59.416430Z", "iopub.status.idle": "2026-05-19T20:34:59.520045Z", "shell.execute_reply": "2026-05-19T20:34:59.519486Z" } }, "outputs": [], "source": [ "model2_tmp = ROOT.RooRealSumPdf(\"sp_ph\", \"sp_ph\", [p_ph_sig2, p_ph_bkg2], [Asig2, Abkg2], True)" ] }, { "cell_type": "markdown", "id": "377cbaf8", "metadata": {}, "source": [ "Construct the subsidiary poisson measurements constraining the statistical fluctuations" ] }, { "cell_type": "code", "execution_count": 14, "id": "c7cae29c", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.522190Z", "iopub.status.busy": "2026-05-19T20:34:59.522066Z", "iopub.status.idle": "2026-05-19T20:34:59.635764Z", "shell.execute_reply": "2026-05-19T20:34:59.635118Z" } }, "outputs": [], "source": [ "hc_sigbkg = ROOT.RooHistConstraint(\"hc_sigbkg\", \"hc_sigbkg\", {p_ph_sig2, p_ph_bkg2})" ] }, { "cell_type": "markdown", "id": "c2932e91", "metadata": {}, "source": [ "Construct the joint model" ] }, { "cell_type": "code", "execution_count": 15, "id": "1cfd3362", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.637523Z", "iopub.status.busy": "2026-05-19T20:34:59.637395Z", "iopub.status.idle": "2026-05-19T20:34:59.741283Z", "shell.execute_reply": "2026-05-19T20:34:59.740763Z" } }, "outputs": [], "source": [ "model2 = ROOT.RooProdPdf(\"model2\", \"model2\", hc_sigbkg, Conditional=(model2_tmp, x))" ] }, { "cell_type": "markdown", "id": "ea844132", "metadata": {}, "source": [ "************ Fit all models to data and plot *********************" ] }, { "cell_type": "code", "execution_count": 16, "id": "ad207c0b", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:34:59.743002Z", "iopub.status.busy": "2026-05-19T20:34:59.742874Z", "iopub.status.idle": "2026-05-19T20:35:00.335202Z", "shell.execute_reply": "2026-05-19T20:35:00.334653Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Minimization -- p.d.f. provides expected number of events, including extended term in likelihood.\n", "[#1] INFO:Fitting -- RooAbsPdf::fitTo(model0_over_model0_Int[x]) 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 807.675 μs\n", "[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model0_over_model0_Int[x]_sumData) Summation contains a RooNLLVar, using its error level\n", "[#1] INFO:Minimization -- [fitFCN] No discrete parameters, performing continuous minimization only\n", "Minuit2Minimizer: Minimize with max-calls 1000 convergence for edm < 1 strategy 1\n", "Minuit2Minimizer : Valid minimum - status = 0\n", "FVAL = -2388.31039421315609\n", "Edm = 3.25299602361765145e-06\n", "Nfcn = 60\n", "Abkg\t = 0.0614547\t +/- 0.00211669\t(limited)\n", "Asig\t = 0.833778\t +/- 0.1898\t(limited)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Minimization -- p.d.f. provides expected number of events, including extended term in likelihood.\n", "[#1] INFO:Minimization -- Including the following constraint terms in minimization: (hc_sig,hc_bkg)\n", "[#1] INFO:Minimization -- The global observables are not defined , normalize constraints with respect to the parameters (Abkg,Asig,p_ph_bkg_gamma_bin_0,p_ph_bkg_gamma_bin_1,p_ph_bkg_gamma_bin_10,p_ph_bkg_gamma_bin_11,p_ph_bkg_gamma_bin_12,p_ph_bkg_gamma_bin_13,p_ph_bkg_gamma_bin_14,p_ph_bkg_gamma_bin_15,p_ph_bkg_gamma_bin_16,p_ph_bkg_gamma_bin_17,p_ph_bkg_gamma_bin_18,p_ph_bkg_gamma_bin_19,p_ph_bkg_gamma_bin_2,p_ph_bkg_gamma_bin_20,p_ph_bkg_gamma_bin_21,p_ph_bkg_gamma_bin_22,p_ph_bkg_gamma_bin_23,p_ph_bkg_gamma_bin_24,p_ph_bkg_gamma_bin_3,p_ph_bkg_gamma_bin_4,p_ph_bkg_gamma_bin_5,p_ph_bkg_gamma_bin_6,p_ph_bkg_gamma_bin_7,p_ph_bkg_gamma_bin_8,p_ph_bkg_gamma_bin_9,p_ph_sig_gamma_bin_0,p_ph_sig_gamma_bin_1,p_ph_sig_gamma_bin_10,p_ph_sig_gamma_bin_11,p_ph_sig_gamma_bin_12,p_ph_sig_gamma_bin_13,p_ph_sig_gamma_bin_14,p_ph_sig_gamma_bin_15,p_ph_sig_gamma_bin_16,p_ph_sig_gamma_bin_17,p_ph_sig_gamma_bin_18,p_ph_sig_gamma_bin_19,p_ph_sig_gamma_bin_2,p_ph_sig_gamma_bin_20,p_ph_sig_gamma_bin_21,p_ph_sig_gamma_bin_22,p_ph_sig_gamma_bin_23,p_ph_sig_gamma_bin_24,p_ph_sig_gamma_bin_3,p_ph_sig_gamma_bin_4,p_ph_sig_gamma_bin_5,p_ph_sig_gamma_bin_6,p_ph_sig_gamma_bin_7,p_ph_sig_gamma_bin_8,p_ph_sig_gamma_bin_9)\n", "[#1] INFO:Fitting -- RooAbsPdf::fitTo(model1) fixing normalization set for coefficient determination to observables in data\n", "[#1] INFO:Fitting -- Creation of NLL object took 1.6126 ms\n", "[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model1_sumData) Summation contains a RooNLLVar, using its error level\n", "[#1] INFO:Minimization -- [fitFCN] No discrete parameters, performing continuous minimization only\n", "Minuit2Minimizer: Minimize with max-calls 16000 convergence for edm < 1 strategy 1\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10425.6) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=0.024352\tAsig=0.0100451\tp_ph_bkg_gamma_bin_0=0.351482\tp_ph_bkg_gamma_bin_1=0.352595\tp_ph_bkg_gamma_bin_10=0.342683\tp_ph_bkg_gamma_bin_11=0.368891\tp_ph_bkg_gamma_bin_12=0.38197\tp_ph_bkg_gamma_bin_13=0.365529\tp_ph_bkg_gamma_bin_14=0.358098\tp_ph_bkg_gamma_bin_15=0.358584\tp_ph_bkg_gamma_bin_16=0.367344\tp_ph_bkg_gamma_bin_17=0.356263\tp_ph_bkg_gamma_bin_18=0.362211\tp_ph_bkg_gamma_bin_19=0.340063\tp_ph_bkg_gamma_bin_2=0.35054\tp_ph_bkg_gamma_bin_20=0.352376\tp_ph_bkg_gamma_bin_21=0.351589\tp_ph_bkg_gamma_bin_22=0.351615\tp_ph_bkg_gamma_bin_23=0.344488\tp_ph_bkg_gamma_bin_24=0.365687\tp_ph_bkg_gamma_bin_3=0.355197\tp_ph_bkg_gamma_bin_4=0.347776\tp_ph_bkg_gamma_bin_5=0.348629\tp_ph_bkg_gamma_bin_6=0.357616\tp_ph_bkg_gamma_bin_7=0.351061\tp_ph_bkg_gamma_bin_8=0.342678\tp_ph_bkg_gamma_bin_9=0.357292\tp_ph_sig_gamma_bin_11=0.3339\tp_ph_sig_gamma_bin_12=0.341347\tp_ph_sig_gamma_bin_13=0.333434\tp_ph_sig_gamma_bin_14=0.328932\tp_ph_sig_gamma_bin_15=0.328666\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10425.6) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=0.327409\tAsig=0.26768\tp_ph_bkg_gamma_bin_0=0.63877\tp_ph_bkg_gamma_bin_1=0.639509\tp_ph_bkg_gamma_bin_10=0.632897\tp_ph_bkg_gamma_bin_11=0.650254\tp_ph_bkg_gamma_bin_12=0.658769\tp_ph_bkg_gamma_bin_13=0.648049\tp_ph_bkg_gamma_bin_14=0.643155\tp_ph_bkg_gamma_bin_15=0.643476\tp_ph_bkg_gamma_bin_16=0.64924\tp_ph_bkg_gamma_bin_17=0.641941\tp_ph_bkg_gamma_bin_18=0.645868\tp_ph_bkg_gamma_bin_19=0.63114\tp_ph_bkg_gamma_bin_2=0.638144\tp_ph_bkg_gamma_bin_20=0.639364\tp_ph_bkg_gamma_bin_21=0.638841\tp_ph_bkg_gamma_bin_22=0.638858\tp_ph_bkg_gamma_bin_23=0.634106\tp_ph_bkg_gamma_bin_24=0.648153\tp_ph_bkg_gamma_bin_3=0.641235\tp_ph_bkg_gamma_bin_4=0.636302\tp_ph_bkg_gamma_bin_5=0.636871\tp_ph_bkg_gamma_bin_6=0.642836\tp_ph_bkg_gamma_bin_7=0.63849\tp_ph_bkg_gamma_bin_8=0.632894\tp_ph_bkg_gamma_bin_9=0.642622\tp_ph_sig_gamma_bin_11=0.626988\tp_ph_sig_gamma_bin_12=0.632002\tp_ph_sig_gamma_bin_13=0.626672\tp_ph_sig_gamma_bin_14=0.623623\tp_ph_sig_gamma_bin_15=0.623442\n", "\n", "Minuit2Minimizer : Valid minimum - status = 0\n", "FVAL = -2282.22793579973768\n", "Edm = 6.50103677974530594e-05\n", "Nfcn = 730\n", "Abkg\t = 0.0613644\t +/- 0.00219573\t(limited)\n", "Asig\t = 0.850329\t +/- 0.233134\t(limited)\n", "p_ph_bkg_gamma_bin_0\t = 0.998158\t +/- 0.0370021\t(limited)\n", "p_ph_bkg_gamma_bin_1\t = 1.0004\t +/- 0.0370335\t(limited)\n", "p_ph_bkg_gamma_bin_10\t = 0.980402\t +/- 0.0356998\t(limited)\n", "p_ph_bkg_gamma_bin_11\t = 1.00613\t +/- 0.0401139\t(limited)\n", "p_ph_bkg_gamma_bin_12\t = 1.00446\t +/- 0.0453421\t(limited)\n", "p_ph_bkg_gamma_bin_13\t = 0.991203\t +/- 0.0395193\t(limited)\n", "p_ph_bkg_gamma_bin_14\t = 0.996054\t +/- 0.0385246\t(limited)\n", "p_ph_bkg_gamma_bin_15\t = 1.00908\t +/- 0.0369959\t(limited)\n", "p_ph_bkg_gamma_bin_16\t = 1.03008\t +/- 0.0382325\t(limited)\n", "p_ph_bkg_gamma_bin_17\t = 1.00779\t +/- 0.0364465\t(limited)\n", "p_ph_bkg_gamma_bin_18\t = 1.01976\t +/- 0.0396035\t(limited)\n", "p_ph_bkg_gamma_bin_19\t = 0.975108\t +/- 0.0359113\t(limited)\n", "p_ph_bkg_gamma_bin_2\t = 0.996261\t +/- 0.0365836\t(limited)\n", "p_ph_bkg_gamma_bin_20\t = 0.999962\t +/- 0.0369391\t(limited)\n", "p_ph_bkg_gamma_bin_21\t = 0.998374\t +/- 0.0370495\t(limited)\n", "p_ph_bkg_gamma_bin_22\t = 0.998427\t +/- 0.0380223\t(limited)\n", "p_ph_bkg_gamma_bin_23\t = 0.984048\t +/- 0.0361184\t(limited)\n", "p_ph_bkg_gamma_bin_24\t = 1.02675\t +/- 0.0373117\t(limited)\n", "p_ph_bkg_gamma_bin_3\t = 1.00565\t +/- 0.0386186\t(limited)\n", "p_ph_bkg_gamma_bin_4\t = 0.990685\t +/- 0.0377676\t(limited)\n", "p_ph_bkg_gamma_bin_5\t = 0.992405\t +/- 0.0376505\t(limited)\n", "p_ph_bkg_gamma_bin_6\t = 1.01052\t +/- 0.038689\t(limited)\n", "p_ph_bkg_gamma_bin_7\t = 0.99731\t +/- 0.0359322\t(limited)\n", "p_ph_bkg_gamma_bin_8\t = 0.980392\t +/- 0.0362349\t(limited)\n", "p_ph_bkg_gamma_bin_9\t = 1.00987\t +/- 0.0372551\t(limited)\n", "p_ph_sig_gamma_bin_11\t = 1.09457\t +/- 0.29733\t(limited)\n", "p_ph_sig_gamma_bin_12\t = 1.07037\t +/- 0.213545\t(limited)\n", "p_ph_sig_gamma_bin_13\t = 0.890724\t +/- 0.195929\t(limited)\n", "p_ph_sig_gamma_bin_14\t = 0.947374\t +/- 0.309135\t(limited)\n", "p_ph_sig_gamma_bin_15\t = 1.14132\t +/- 0.81204\t(limited)\n", "[#1] INFO:Minimization -- p.d.f. provides expected number of events, including extended term in likelihood.\n", "[#1] INFO:Minimization -- Including the following constraint terms in minimization: (hc_sigbkg)\n", "[#1] INFO:Minimization -- The global observables are not defined , normalize constraints with respect to the parameters (Abkg,Asig,p_ph_sig2_gamma_bin_0,p_ph_sig2_gamma_bin_1,p_ph_sig2_gamma_bin_10,p_ph_sig2_gamma_bin_11,p_ph_sig2_gamma_bin_12,p_ph_sig2_gamma_bin_13,p_ph_sig2_gamma_bin_14,p_ph_sig2_gamma_bin_15,p_ph_sig2_gamma_bin_16,p_ph_sig2_gamma_bin_17,p_ph_sig2_gamma_bin_18,p_ph_sig2_gamma_bin_19,p_ph_sig2_gamma_bin_2,p_ph_sig2_gamma_bin_20,p_ph_sig2_gamma_bin_21,p_ph_sig2_gamma_bin_22,p_ph_sig2_gamma_bin_23,p_ph_sig2_gamma_bin_24,p_ph_sig2_gamma_bin_3,p_ph_sig2_gamma_bin_4,p_ph_sig2_gamma_bin_5,p_ph_sig2_gamma_bin_6,p_ph_sig2_gamma_bin_7,p_ph_sig2_gamma_bin_8,p_ph_sig2_gamma_bin_9)\n", "[#1] INFO:Fitting -- RooAbsPdf::fitTo(model2) fixing normalization set for coefficient determination to observables in data\n", "[#1] INFO:Fitting -- Creation of NLL object took 741.894 μs\n", "[#1] INFO:Fitting -- RooAddition::defaultErrorLevel(nll_model2_sumData) Summation contains a RooNLLVar, using its error level\n", "[#1] INFO:Minimization -- [fitFCN] No discrete parameters, performing continuous minimization only\n", "Minuit2Minimizer: Minimize with max-calls 13500 convergence for edm < 1 strategy 1\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=7.28945\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=7.28945\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=7.28945\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=7.28945\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=7.28945\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=7.28945\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=7.28945\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=7.28945\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=7.28945\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=7.28945\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=7.28945\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=7.28945\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=7.28945\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=7.28945\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=7.28945\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=7.28945\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=7.28945\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=7.28945\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=7.28945\tp_ph_sig2_gamma_bin_9=1\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=1\tAsig=1\tp_ph_sig2_gamma_bin_0=1\tp_ph_sig2_gamma_bin_1=1\tp_ph_sig2_gamma_bin_10=1\tp_ph_sig2_gamma_bin_11=1\tp_ph_sig2_gamma_bin_12=1\tp_ph_sig2_gamma_bin_13=1\tp_ph_sig2_gamma_bin_14=1\tp_ph_sig2_gamma_bin_15=1\tp_ph_sig2_gamma_bin_16=1\tp_ph_sig2_gamma_bin_17=1\tp_ph_sig2_gamma_bin_18=1\tp_ph_sig2_gamma_bin_19=1\tp_ph_sig2_gamma_bin_2=1\tp_ph_sig2_gamma_bin_20=1\tp_ph_sig2_gamma_bin_21=1\tp_ph_sig2_gamma_bin_22=1\tp_ph_sig2_gamma_bin_23=1\tp_ph_sig2_gamma_bin_24=1\tp_ph_sig2_gamma_bin_3=1\tp_ph_sig2_gamma_bin_4=1\tp_ph_sig2_gamma_bin_5=1\tp_ph_sig2_gamma_bin_6=1\tp_ph_sig2_gamma_bin_7=1\tp_ph_sig2_gamma_bin_8=1\tp_ph_sig2_gamma_bin_9=7.28945\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=0.024352\tAsig=0.0100451\tp_ph_sig2_gamma_bin_0=0.351536\tp_ph_sig2_gamma_bin_1=0.352649\tp_ph_sig2_gamma_bin_10=0.342736\tp_ph_sig2_gamma_bin_11=0.36962\tp_ph_sig2_gamma_bin_12=0.383863\tp_ph_sig2_gamma_bin_13=0.366897\tp_ph_sig2_gamma_bin_14=0.358549\tp_ph_sig2_gamma_bin_15=0.358659\tp_ph_sig2_gamma_bin_16=0.3674\tp_ph_sig2_gamma_bin_17=0.356318\tp_ph_sig2_gamma_bin_18=0.362267\tp_ph_sig2_gamma_bin_19=0.340116\tp_ph_sig2_gamma_bin_2=0.350595\tp_ph_sig2_gamma_bin_20=0.352431\tp_ph_sig2_gamma_bin_21=0.351643\tp_ph_sig2_gamma_bin_22=0.35167\tp_ph_sig2_gamma_bin_23=0.344542\tp_ph_sig2_gamma_bin_24=0.365744\tp_ph_sig2_gamma_bin_3=0.355252\tp_ph_sig2_gamma_bin_4=0.34783\tp_ph_sig2_gamma_bin_5=0.348683\tp_ph_sig2_gamma_bin_6=0.357671\tp_ph_sig2_gamma_bin_7=0.351115\tp_ph_sig2_gamma_bin_8=0.342731\tp_ph_sig2_gamma_bin_9=0.357348\n", "\n", "RooAbsMinimizerFcn: Minimized function has error status.\n", "Returning maximum FCN so far (10415.9) to force MIGRAD to back out of this region. Error log follows.\n", "Parameter values: \tAbkg=0.327407\tAsig=0.267678\tp_ph_sig2_gamma_bin_0=0.638805\tp_ph_sig2_gamma_bin_1=0.639545\tp_ph_sig2_gamma_bin_10=0.632932\tp_ph_sig2_gamma_bin_11=0.650729\tp_ph_sig2_gamma_bin_12=0.659993\tp_ph_sig2_gamma_bin_13=0.648946\tp_ph_sig2_gamma_bin_14=0.643452\tp_ph_sig2_gamma_bin_15=0.643524\tp_ph_sig2_gamma_bin_16=0.649276\tp_ph_sig2_gamma_bin_17=0.641977\tp_ph_sig2_gamma_bin_18=0.645904\tp_ph_sig2_gamma_bin_19=0.631174\tp_ph_sig2_gamma_bin_2=0.638179\tp_ph_sig2_gamma_bin_20=0.6394\tp_ph_sig2_gamma_bin_21=0.638876\tp_ph_sig2_gamma_bin_22=0.638894\tp_ph_sig2_gamma_bin_23=0.634141\tp_ph_sig2_gamma_bin_24=0.648189\tp_ph_sig2_gamma_bin_3=0.641271\tp_ph_sig2_gamma_bin_4=0.636338\tp_ph_sig2_gamma_bin_5=0.636906\tp_ph_sig2_gamma_bin_6=0.642872\tp_ph_sig2_gamma_bin_7=0.638525\tp_ph_sig2_gamma_bin_8=0.632929\tp_ph_sig2_gamma_bin_9=0.642658\n", "\n", "Minuit2Minimizer : Valid minimum - status = 0\n", "FVAL = -2291.45108203063864\n", "Edm = 0.00021711856797120699\n", "Nfcn = 1108\n", "Abkg\t = 0.0614871\t +/- 0.00221998\t(limited)\n", "Asig\t = 0.834474\t +/- 0.202991\t(limited)\n", "p_ph_sig2_gamma_bin_0\t = 0.997805\t +/- 0.0474115\t(limited)\n", "p_ph_sig2_gamma_bin_1\t = 1.00004\t +/- 0.0474649\t(limited)\n", "p_ph_sig2_gamma_bin_10\t = 0.980078\t +/- 0.0456101\t(limited)\n", "p_ph_sig2_gamma_bin_11\t = 1.01069\t +/- 0.0488754\t(limited)\n", "p_ph_sig2_gamma_bin_12\t = 1.01199\t +/- 0.0483997\t(limited)\n", "p_ph_sig2_gamma_bin_13\t = 0.983813\t +/- 0.0467055\t(limited)\n", "p_ph_sig2_gamma_bin_14\t = 0.994769\t +/- 0.0483032\t(limited)\n", "p_ph_sig2_gamma_bin_15\t = 1.0095\t +/- 0.0473249\t(limited)\n", "p_ph_sig2_gamma_bin_16\t = 1.02967\t +/- 0.049202\t(limited)\n", "p_ph_sig2_gamma_bin_17\t = 1.00742\t +/- 0.0467328\t(limited)\n", "p_ph_sig2_gamma_bin_18\t = 1.01937\t +/- 0.0509499\t(limited)\n", "p_ph_sig2_gamma_bin_19\t = 0.974793\t +/- 0.0458586\t(limited)\n", "p_ph_sig2_gamma_bin_2\t = 0.99591\t +/- 0.0468523\t(limited)\n", "p_ph_sig2_gamma_bin_20\t = 0.999605\t +/- 0.0473385\t(limited)\n", "p_ph_sig2_gamma_bin_21\t = 0.99802\t +/- 0.0474749\t(limited)\n", "p_ph_sig2_gamma_bin_22\t = 0.998073\t +/- 0.0487512\t(limited)\n", "p_ph_sig2_gamma_bin_23\t = 0.983718\t +/- 0.0461774\t(limited)\n", "p_ph_sig2_gamma_bin_24\t = 1.02635\t +/- 0.0479702\t(limited)\n", "p_ph_sig2_gamma_bin_3\t = 1.00528\t +/- 0.0495745\t(limited)\n", "p_ph_sig2_gamma_bin_4\t = 0.990344\t +/- 0.0483736\t(limited)\n", "p_ph_sig2_gamma_bin_5\t = 0.992062\t +/- 0.0482299\t(limited)\n", "p_ph_sig2_gamma_bin_6\t = 1.01014\t +/- 0.0496947\t(limited)\n", "p_ph_sig2_gamma_bin_7\t = 0.996957\t +/- 0.0460029\t(limited)\n", "p_ph_sig2_gamma_bin_8\t = 0.980068\t +/- 0.0463103\t(limited)\n", "p_ph_sig2_gamma_bin_9\t = 1.00949\t +/- 0.0478068\t(limited)\n" ] } ], "source": [ "result0 = model0.fitTo(sumData, PrintLevel=0, Save=True)\n", "result1 = model1.fitTo(sumData, PrintLevel=0, Save=True)\n", "result2 = model2.fitTo(sumData, PrintLevel=0, Save=True)\n", "\n", "\n", "can = ROOT.TCanvas(\"can\", \"\", 1500, 600)\n", "can.Divide(3, 1)\n", "\n", "pt = ROOT.TPaveText(-19.5, 1, -2, 25)\n", "pt.SetFillStyle(0)\n", "pt.SetBorderSize(0)\n", "\n", "\n", "can.cd(1)\n", "frame = x.frame(Title=\"No template uncertainties\")" ] }, { "cell_type": "markdown", "id": "6f3d3b93", "metadata": {}, "source": [ "Plot data to enable automatic determination of model0 normalisation:" ] }, { "cell_type": "code", "execution_count": 17, "id": "28ac45a5", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:00.336981Z", "iopub.status.busy": "2026-05-19T20:35:00.336854Z", "iopub.status.idle": "2026-05-19T20:35:00.504419Z", "shell.execute_reply": "2026-05-19T20:35:00.503878Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model0) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model0) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model0) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model0) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model0) WARNING: argument CurveNameSuffix is duplicated\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sumData.plotOn(frame)\n", "model0.plotOn(frame, LineColor=\"b\", VisualizeError=result0)" ] }, { "cell_type": "markdown", "id": "38db0072", "metadata": {}, "source": [ "Plot data again to show it on top of model0 error bands:" ] }, { "cell_type": "code", "execution_count": 18, "id": "b7014d7c", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:00.506240Z", "iopub.status.busy": "2026-05-19T20:35:00.506105Z", "iopub.status.idle": "2026-05-19T20:35:00.610582Z", "shell.execute_reply": "2026-05-19T20:35:00.610040Z" } }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 18, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sumData.plotOn(frame)" ] }, { "cell_type": "markdown", "id": "1a507073", "metadata": {}, "source": [ "Plot model components" ] }, { "cell_type": "code", "execution_count": 19, "id": "d63b299c", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:00.612290Z", "iopub.status.busy": "2026-05-19T20:35:00.612166Z", "iopub.status.idle": "2026-05-19T20:35:00.809051Z", "shell.execute_reply": "2026-05-19T20:35:00.808516Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model0) directly selected PDF components: (p_h_sig)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model0) indirectly selected PDF components: ()\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model0) directly selected PDF components: (p_h_bkg)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model0) indirectly selected PDF components: ()\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 50 will supersede previous event count of 1050 for normalization of PDF projections\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model1) WARNING: argument CurveNameSuffix is duplicated\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 19, "metadata": {}, "output_type": "execute_result" } ], "source": [ "model0.plotOn(frame, LineColor=\"b\")\n", "p_ph_sig_set = {p_h_sig}\n", "p_ph_bkg_set = {p_h_bkg}\n", "model0.plotOn(frame, Components=p_ph_sig_set, LineColor=\"kAzure\")\n", "model0.plotOn(frame, Components=p_ph_bkg_set, LineColor=\"r\")\n", "model0.paramOn(frame)\n", "\n", "sigData.plotOn(frame, MarkerColor=\"b\")\n", "frame.Draw()\n", "\n", "pt_text1 = [\n", " \"No template uncertainties\",\n", " \"are taken into account.\",\n", " \"This leads to low errors\",\n", " \"for the parameters A, since\",\n", " \"the only source of errors\",\n", " \"are the data statistics.\",\n", "]\n", "for text in pt_text1:\n", " pt.AddText(text)\n", "\n", "pt.DrawClone()\n", "\n", "\n", "can.cd(2)\n", "frame = x.frame(Title=\"Barlow Beeston for Sig & Bkg separately\")\n", "sumData.plotOn(frame)\n", "model1.plotOn(frame, LineColor=\"b\", VisualizeError=result1)" ] }, { "cell_type": "markdown", "id": "74b3d971", "metadata": {}, "source": [ "Plot data again to show it on top of error bands:" ] }, { "cell_type": "code", "execution_count": 20, "id": "546d7806", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:00.810761Z", "iopub.status.busy": "2026-05-19T20:35:00.810629Z", "iopub.status.idle": "2026-05-19T20:35:00.970722Z", "shell.execute_reply": "2026-05-19T20:35:00.970198Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model1) directly selected PDF components: (p_ph_sig)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model1) indirectly selected PDF components: (sp_ph)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model1) directly selected PDF components: (p_ph_bkg)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model1) indirectly selected PDF components: (sp_ph)\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "name": "stdout", "output_type": "stream", "text": [ "[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 50 will supersede previous event count of 1050 for normalization of PDF projections\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n", "[#0] WARNING:InputArguments -- RooAbsReal::plotOn(model2) WARNING: argument CurveNameSuffix is duplicated\n" ] } ], "source": [ "sumData.plotOn(frame)\n", "model1.plotOn(frame, LineColor=\"b\")\n", "p_ph_sig1_set = {p_ph_sig1}\n", "p_ph_bkg1_set = {p_ph_bkg1}\n", "model1.plotOn(frame, Components=p_ph_sig1_set, LineColor=\"kAzure\")\n", "model1.plotOn(frame, Components=p_ph_bkg1_set, LineColor=\"r\")\n", "model1.paramOn(frame, Parameters={Asig1, Abkg1})\n", "\n", "sigData.plotOn(frame, MarkerColor=\"b\")\n", "frame.Draw()\n", "\n", "pt.Clear()\n", "pt_text2 = [\n", " \"With gamma parameters, the\",\n", " \"signal & background templates\",\n", " \"can adapt to the data.\",\n", " \"Note how the blue signal\",\n", " \"template changes its shape.\",\n", " \"This leads to higher errors\",\n", " \"of the scale parameters A.\",\n", "]\n", "\n", "for text in pt_text2:\n", " pt.AddText(text)\n", "\n", "pt.DrawClone()\n", "\n", "can.cd(3)\n", "frame = x.frame(Title=\"Barlow Beeston light for (Sig+Bkg)\")\n", "sumData.plotOn(frame)\n", "model2.plotOn(frame, LineColor=\"b\", VisualizeError=result2)" ] }, { "cell_type": "markdown", "id": "71e63ef7", "metadata": {}, "source": [ "Plot data again to show it on top of model0 error bands:" ] }, { "cell_type": "code", "execution_count": 21, "id": "593ed241", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:00.972203Z", "iopub.status.busy": "2026-05-19T20:35:00.972081Z", "iopub.status.idle": "2026-05-19T20:35:01.233358Z", "shell.execute_reply": "2026-05-19T20:35:01.232839Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Asig [normal ] = 0.8337778710558245 +/- 0.18981418763032815\n", "Asig [BB ] = 0.8503293171822359 +/- 0.2357828123142156\n", "Asig [BBlight] = 0.8344736023878351 +/- 0.20295518584642697\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model2) directly selected PDF components: (p_ph_sig2)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model2) indirectly selected PDF components: (sp_ph)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model2) directly selected PDF components: (p_ph_bkg2)\n", "[#1] INFO:Plotting -- RooAbsPdf::plotOn(model2) indirectly selected PDF components: (sp_ph)\n", "[#1] INFO:Plotting -- RooPlot::updateFitRangeNorm: New event count of 50 will supersede previous event count of 1050 for normalization of PDF projections\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "Info in : png file rf709_BarlowBeeston.png has been created\n" ] } ], "source": [ "sumData.plotOn(frame)\n", "model2.plotOn(frame, LineColor=\"b\")\n", "p_ph_sig2_set = {p_ph_sig2}\n", "p_ph_bkg2_set = {p_ph_bkg2}\n", "model2.plotOn(frame, Components=p_ph_sig2_set, LineColor=\"kAzure\")\n", "model2.plotOn(frame, Components=p_ph_bkg2_set, LineColor=\"r\")\n", "model2.paramOn(frame, Parameters={Asig2, Abkg2})\n", "\n", "sigData.plotOn(frame, MarkerColor=\"b\")\n", "frame.Draw()\n", "\n", "pt.Clear()\n", "pt_text3 = [\n", " \"When signal and background\",\n", " \"template share one gamma para-\",\n", " \"meter per bin, they adapt less.\",\n", " \"The errors of the A parameters\",\n", " \"also shrink slightly.\",\n", "]\n", "for text in pt_text3:\n", " pt.AddText(text)\n", "pt.DrawClone()\n", "\n", "\n", "print(\"Asig [normal ] = {} +/- {}\".format(Asig0.getVal(), Asig0.getError()))\n", "print(\"Asig [BB ] = {} +/- {}\".format(Asig1.getVal(), Asig1.getError()))\n", "print(\"Asig [BBlight] = {} +/- {}\".format(Asig2.getVal(), Asig2.getError()))\n", "\n", "can.SaveAs(\"rf709_BarlowBeeston.png\")" ] }, { "cell_type": "markdown", "id": "7165fbeb", "metadata": {}, "source": [ "Draw all canvases " ] }, { "cell_type": "code", "execution_count": 22, "id": "bc999203", "metadata": { "collapsed": false, "execution": { "iopub.execute_input": "2026-05-19T20:35:01.235079Z", "iopub.status.busy": "2026-05-19T20:35:01.234954Z", "iopub.status.idle": "2026-05-19T20:35:01.421896Z", "shell.execute_reply": "2026-05-19T20:35:01.421454Z" } }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "
\n", "
\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%jsroot on\n", "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 }