doc/master/rf804__mcstudy__constr_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussian.h"

#include "RooPolynomial.h"

#include "RooAddPdf.h"

#include "RooProdPdf.h"

#include "RooMCStudy.h"

#include "RooPlot.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "TH1.h"

using namespace RooFit;


void rf804_mcstudy_constr()

{

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

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


   // Observable

   RooRealVar x("x", "x", -10, 10);


   // Signal component

   RooRealVar m("m", "m", 0, -10, 10);

   RooRealVar s("s", "s", 2, 0.1, 10);

   RooGaussian g("g", "g", x, m, s);


   // Background component

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


   // Composite model

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

   RooAddPdf sum("sum", "sum", RooArgSet(g, p), f);


   // Construct constraint on parameter f

   RooGaussian fconstraint("fconstraint", "fconstraint", f, 0.7, 0.1);


   // Multiply constraint with pdf

   RooProdPdf sumc("sumc", "sum with constraint", RooArgSet(sum, fconstraint));


   // S e t u p   t o y   s t u d y   w i t h   m o d e l

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


   // Perform toy study with internal constraint on f

   RooMCStudy mcs(sumc, x, Constrain(f), Silence(), Binned(), FitOptions(PrintLevel(-1)));


   // Run 500 toys of 2000 events.

   // Before each toy is generated, a value for the f is sampled from the constraint pdf and

   // that value is used for the generation of that toy.

   mcs.generateAndFit(500, 2000);


   // Make plot of distribution of generated value of f parameter

   TH1 *h_f_gen = mcs.fitParDataSet().createHistogram("f_gen", AutoBinning(40));


   // Make plot of distribution of fitted value of f parameter

   RooPlot *frame1 = mcs.plotParam(f, Bins(40));

   frame1->SetTitle("Distribution of fitted f values");


   // Make plot of pull distribution on f

   RooPlot *frame2 = mcs.plotPull(f, Bins(40), FitGauss());

   frame1->SetTitle("Distribution of f pull values");


   TCanvas *c = new TCanvas("rf804_mcstudy_constr", "rf804_mcstudy_constr", 1200, 400);

   c->Divide(3);

   c->cd(1);

   gPad->SetLeftMargin(0.15);

   h_f_gen->GetYaxis()->SetTitleOffset(1.4);

   h_f_gen->Draw();

   c->cd(2);

   gPad->SetLeftMargin(0.15);

   frame1->GetYaxis()->SetTitleOffset(1.4);

   frame1->Draw();

   c->cd(3);

   gPad->SetLeftMargin(0.15);

   frame2->GetYaxis()->SetTitleOffset(1.4);

   frame2->Draw();

}

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

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

g
#define g(i)
Definition RSha256.hxx:105

RooAddPdf.h

RooDataSet.h

RooGaussian.h

RooMCStudy.h

RooPlot.h

RooPolynomial.h

RooProdPdf.h

RooRealVar.h

TAxis.h

TCanvas.h

p
winID h TVirtualViewer3D TVirtualGLPainter p
Definition TGWin32VirtualGLProxy.cxx:51

TH1.h

gPad
#define gPad
Definition TVirtualPad.h:308

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:24

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

RooMCStudy
Helper class to facilitate Monte Carlo studies such as 'goodness-of-fit' studies, that involve fittin...
Definition RooMCStudy.h:32

RooPlot
Plot frame and a container for graphics objects within that frame.
Definition RooPlot.h:43

RooPlot::SetTitle
void SetTitle(const char *name) override
Set the title of the RooPlot to 'title'.
Definition RooPlot.cxx:1203

RooPlot::GetYaxis
TAxis * GetYaxis() const
Definition RooPlot.cxx:1224

RooPlot::Draw
void Draw(Option_t *options=nullptr) override
Draw this plot and all of the elements it contains.
Definition RooPlot.cxx:597

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

RooProdPdf
Efficient implementation of a product of PDFs of the form.
Definition RooProdPdf.h:39

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

TAttAxis::SetTitleOffset
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
Definition TAttAxis.cxx:298

TCanvas
The Canvas class.
Definition TCanvas.h:23

TH1
TH1 is the base class of all histogram classes in ROOT.
Definition TH1.h:59

TH1::GetYaxis
TAxis * GetYaxis()
Definition TH1.h:326

TH1::Draw
void Draw(Option_t *option="") override
Draw this histogram with options.
Definition TH1.cxx:3068

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 CodegenImpl.h:64

rf804_mcstudy_constr
Definition rf804_mcstudy_constr.py:1

m
TMarker m
Definition textangle.C:8

sum
static uint64_t sum(uint64_t i)
Definition Factory.cxx:2345