doc/v632/rf504__simwstool_8C.html

#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;


void rf504_simwstool()

{

   // C r e a t e   m a s t e r   p d f

   // ---------------------------------


   // 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);


   // Construct poly(x,p0)

   RooRealVar p0("p0", "p0", 0.01, 0., 1.);

   RooPolynomial poly("p", "p", x, p0);


   // Construct model = f*gauss(x) + (1-f)*poly(x)

   RooRealVar f("f", "f", 0.5, 0., 1.);

   RooAddPdf model("model", "model", RooArgSet(gauss, poly), f);


   // C r e a t e   c a t e g o r y   o b s e r v a b l e s   f o r   s p l i t t i n g

   // ----------------------------------------------------------------------------------


   // 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");


   // S e t u p   S i m W S T o o l

   // -----------------------------


   // Import ingredients in a workspace

   RooWorkspace w("w", "w");

   w.import(RooArgSet(model, c, d));


   // Make Sim builder tool

   RooSimWSTool sct(w);


   // B u i l d   a   s i m u l t a n e o u s   m o d e l   w i t h   o n e   s p l i t

   // ---------------------------------------------------------------------------------


   // Construct a simultaneous pdf with the following form

   //

   // model_run1(x) = f*gauss_run1(x,m_run1,s) + (1-f)*poly

   // model_run2(x) = f*gauss_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"));


   // Print tree structure of model

   model_sim->Print("t");


   // 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");


   // B u i l d   a   s i m u l t a n e o u s   m o d e l   w i t h   p r o d u c t   s p l i t

   // -----------------------------------------------------------------------------------------


   // 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"));


   // Print tree structure of this model

   model_sim2->Print("t");

}

d
#define d(i)
Definition RSha256.hxx:102

f
#define f(i)
Definition RSha256.hxx:104

c
#define c(i)
Definition RSha256.hxx:101

RooAddPdf.h

RooCategory.h

RooDataSet.h

RooGaussian.h

RooPlot.h

RooPolynomial.h

RooRealVar.h

RooSimWSTool.h

RooSimultaneous.h

RooWorkspace.h

TAxis.h

TCanvas.h

TFile.h

w
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void w
Definition TGWin32VirtualXProxy.cxx:138

TH1.h

RooAbsArg::Print
void Print(Option_t *options=nullptr) const override
Print the object to the defaultPrintStream().
Definition RooAbsArg.h:320

RooAddPdf
Efficient implementation of a sum of PDFs of the form.
Definition RooAddPdf.h:33

RooArgSet
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:55

RooCategory
Object to represent discrete states.
Definition RooCategory.h:28

RooGaussian
Plain Gaussian p.d.f.
Definition RooGaussian.h:24

RooPolynomial
RooPolynomial implements a polynomial p.d.f of the form.
Definition RooPolynomial.h:25

RooRealVar
Variable that can be changed from the outside.
Definition RooRealVar.h:37

RooSimWSTool
The RooSimWSTool is a tool operating on RooWorkspace objects that can clone PDFs into a series of var...
Definition RooSimWSTool.h:38

RooSimultaneous
Facilitates simultaneous fitting of multiple PDFs to subsets of a given dataset.
Definition RooSimultaneous.h:41

RooWorkspace
Persistable container for RooFit projects.
Definition RooWorkspace.h:43

RooFit::SplitParam
RooCmdArg SplitParam(const char *varname, const char *catname)
Definition RooGlobalFunc.cxx:998

x
Double_t x[n]
Definition legend1.C:17

RooFit
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition JSONIO.h:26

rf504_simwstool
Definition rf504_simwstool.py:1

m
TMarker m
Definition textangle.C:8