Organisation and simultaneous fits: using RooSimWSTool to construct a simultaneous pdf that is built of variations of an input pdf
Author: Wouter Verkerke
This notebook tutorial was automatically generated with ROOTBOOK-izer from the macro found in the ROOT repository on Tuesday, March 19, 2024 at 07:16 PM.
%%cpp -d
#include "RooRealVar.h"
#include "RooCategory.h"
#include "RooDataSet.h"
#include "RooGaussian.h"
#include "RooPolynomial.h"
#include "RooSimultaneous.h"
#include "RooAddPdf.h"
#include "RooWorkspace.h"
#include "RooSimWSTool.h"
#include "RooPlot.h"
#include "TCanvas.h"
#include "TAxis.h"
#include "TFile.h"
#include "TH1.h"
using namespace RooFit;
Construct gauss(x,m,s)
RooRealVar x("x", "x", -10, 10);
RooRealVar m("m", "m", 0, -10, 10);
RooRealVar s("s", "s", 1, -10, 10);
RooGaussian gauss("g", "g", x, m, s);
[#0] WARNING:InputArguments -- The parameter 's' with range [-10, 10] of the RooGaussian 'g' exceeds the safe range of (0, inf). Advise to limit its range.
Construct poly(x,p0)
RooRealVar p0("p0", "p0", 0.01, 0., 1.);
RooPolynomial poly("p", "p", x, p0);
Construct model = fgauss(x) + (1-f)poly(x)
RooRealVar f("f", "f", 0.5, 0., 1.);
RooAddPdf model("model", "model", RooArgSet(gauss, poly), f);
Define two categories that can be used for splitting
RooCategory c("c", "c");
c.defineType("run1");
c.defineType("run2");
RooCategory d("d", "d");
d.defineType("foo");
d.defineType("bar");
Import ingredients in a workspace
RooWorkspace w("w", "w");
w.import(RooArgSet(model, c, d));
[#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooGaussian::g [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::x [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::m [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::s [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::f [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooPolynomial::p [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::p0 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooCategory::c [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooCategory::d
Make Sim builder tool
RooSimWSTool sct(w);
Construct a simultaneous pdf with the following form
model_run1(x) = fgauss_run1(x,m_run1,s) + (1-f)poly model_run2(x) = fgauss_run2(x,m_run2,s) + (1-f)poly simpdf(x,c) = model_run1(x) if c=="run1" = model_run2(x) if c=="run2"
Returned pdf is owned by the workspace
RooSimultaneous *model_sim = sct.build("model_sim", "model", SplitParam("m", "c"));
Splitrule for p.d.f model with state list parameter m is split with constraint in categories (c) [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: list of prototype pdfs (model) [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: list of splitting categories (c) [#1] INFO:ObjectHandling -- RooSimPdfBuilder::executeBuild: processing prototype pdf model [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: configured customizers for all prototype pdfs [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode run1 [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode run2 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooSimultaneous::model_sim [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_run1 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooGaussian::g_run1 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::m_run1 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_run2 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooGaussian::g_run2 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::m_run2
Print tree structure of model
model_sim->Print("t");
0x7f34c194fbb0 RooSimultaneous::model_sim = 1 [Auto,Dirty] 0x7f34c18ffa00/V- RooCategory::c = run2(idx = 1) 0x7f34c1885d00/V- RooAddPdf::model_run1 = 1/1 [Auto,Clean] 0x7f34c189b040/V- RooGaussian::g_run1 = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1984360/V- RooRealVar::m_run1 = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1797620/V- RooPolynomial::p = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c18e3a50/V- RooRealVar::p0 = 0.01 0x7f34c198dd60/V- RooAddPdf::model_run2 = 1/1 [Auto,Clean] 0x7f34c1993320/V- RooGaussian::g_run2 = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c19958d0/V- RooRealVar::m_run2 = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1797620/V- RooPolynomial::p = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c18e3a50/V- RooRealVar::p0 = 0.01
Adjust model_sim parameters in workspace
w.var("m_run1")->setVal(-3);
w.var("m_run2")->setVal(+3);
Print contents of workspace
w.Print("v");
RooWorkspace(w) w contents variables --------- (c,d,f,m,m_run1,m_run2,p0,s,x) p.d.f.s ------- RooGaussian::g[ x=x mean=m sigma=s ] = 1 RooGaussian::g_run1[ x=x mean=m_run1 sigma=s ] = 0.011109 RooGaussian::g_run2[ x=x mean=m_run2 sigma=s ] = 0.011109 RooAddPdf::model[ f * g + [%] * p ] = 1/1 RooAddPdf::model_run1[ f * g_run1 + [%] * p ] = 0.505554/1 RooAddPdf::model_run2[ f * g_run2 + [%] * p ] = 0.505554/1 RooSimultaneous::model_sim[ indexCat=c run1=model_run1 run2=model_run2 ] = 0.505554 RooPolynomial::p[ x=x coefList=(p0) ] = 1
Build another simultaneous pdf using a composite split in states c X d
RooSimultaneous *model_sim2 = sct.build("model_sim2", "model", SplitParam("p0", "c,d"));
Splitrule for p.d.f model with state list parameter p0 is split with constraint in categories (c,d) [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: list of prototype pdfs (model) [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: list of splitting categories (c,d) [#1] INFO:ObjectHandling -- RooSimPdfBuilder::executeBuild: processing prototype pdf model [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooMultiCategory::c,d [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: configured customizers for all prototype pdfs [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode {run1;bar} [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode {run1;foo} [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode {run2;bar} [#1] INFO:ObjectHandling -- RooSimWSTool::executeBuild: Customizing prototype pdf model for mode {run2;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooSimultaneous::model_sim2 [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooSuperCategory::model_sim2_index [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooMultiCategory::model_sim2_index_internalMultiCat [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_{run1;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooPolynomial::p_{run1;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::p0_{run1;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_{run1;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooPolynomial::p_{run1;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::p0_{run1;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_{run2;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooPolynomial::p_{run2;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::p0_{run2;bar} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooAddPdf::model_{run2;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooPolynomial::p_{run2;foo} [#1] INFO:ObjectHandling -- RooWorkspace::import(w) importing RooRealVar::p0_{run2;foo}
Print tree structure of this model
model_sim2->Print("t");
0x7f34c1afdef0 RooSimultaneous::model_sim2 = 1 [Auto,Dirty] 0x7f34c1b025b0/V- RooSuperCategory::model_sim2_index = {run2;foo}(idx = 1) [Auto,Clean] 0x7f34c1b029b0/VS RooMultiCategory::model_sim2_index_internalMultiCat = {run2;foo}(idx = 1) [Auto,Clean] 0x7f34c18ffa00/VS RooCategory::c = run2(idx = 1) 0x7f34c18c7060/VS RooCategory::d = foo(idx = 0) 0x7f34c1b02de0/V- RooAddPdf::model_{run1;bar} = 1/1 [Auto,Clean] 0x7f34c17e1710/V- RooGaussian::g = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1892820/V- RooRealVar::m = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1b0e030/V- RooPolynomial::p_{run1;bar} = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1b105c0/V- RooRealVar::p0_{run1;bar} = 0.01 0x7f34c1b10c40/V- RooAddPdf::model_{run1;foo} = 1/1 [Auto,Clean] 0x7f34c17e1710/V- RooGaussian::g = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1892820/V- RooRealVar::m = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1b15810/V- RooPolynomial::p_{run1;foo} = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1b19590/V- RooRealVar::p0_{run1;foo} = 0.01 0x7f34c1b19c10/V- RooAddPdf::model_{run2;bar} = 1/1 [Auto,Clean] 0x7f34c17e1710/V- RooGaussian::g = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1892820/V- RooRealVar::m = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1b1e7e0/V- RooPolynomial::p_{run2;bar} = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1b20d50/V- RooRealVar::p0_{run2;bar} = 0.01 0x7f34c1b213d0/V- RooAddPdf::model_{run2;foo} = 1/1 [Auto,Clean] 0x7f34c17e1710/V- RooGaussian::g = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1892820/V- RooRealVar::m = 0 0x7f34c18d6070/V- RooRealVar::s = 1 0x7f34c194cd30/V- RooRealVar::f = 0.5 0x7f34c1b26060/V- RooPolynomial::p_{run2;foo} = 1 [Auto,Dirty] 0x7f34c1783690/V- RooRealVar::x = 0 0x7f34c1b2b290/V- RooRealVar::p0_{run2;foo} = 0.01