Logo ROOT  
Reference Guide
rf511_wsfactory_basic.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_roofit
3/// \notebook -nodraw
4/// Organization and simultaneous fits: basic use of the 'object factory' associated with
5/// a workspace to rapidly build pdfs functions and their parameter components
6///
7/// \macro_output
8/// \macro_code
9///
10/// \date July 2009
11/// \author Wouter Verkerke
12
13#include "RooRealVar.h"
14#include "RooDataSet.h"
15#include "RooGaussian.h"
16#include "RooConstVar.h"
17#include "RooChebychev.h"
18#include "RooAddPdf.h"
19#include "RooWorkspace.h"
20#include "RooPlot.h"
21#include "TCanvas.h"
22#include "TAxis.h"
23using namespace RooFit;
24
26{
27 RooWorkspace *w = new RooWorkspace("w");
28
29 // C r e a t i n g a n d a d d i n g b a s i c p . d . f . s
30 // ----------------------------------------------------------------
31
32 // Remake example pdf of tutorial rs502_wspacewrite.C:
33 //
34 // Basic pdf construction: ClassName::ObjectName(constructor arguments)
35 // Variable construction : VarName[x,xlo,xhi], VarName[xlo,xhi], VarName[x]
36 // P.d.f. addition : SUM::ObjectName(coef1*pdf1,...coefM*pdfM,pdfN)
37 //
38
39 if (!compact) {
40
41 // Use object factory to build pdf of tutorial rs502_wspacewrite
42 w->factory("Gaussian::sig1(x[-10,10],mean[5,0,10],0.5)");
43 w->factory("Gaussian::sig2(x,mean,1)");
44 w->factory("Chebychev::bkg(x,{a0[0.5,0.,1],a1[0.2,0.,1.]})");
45 w->factory("SUM::sig(sig1frac[0.8,0.,1.]*sig1,sig2)");
46 w->factory("SUM::model(bkgfrac[0.5,0.,1.]*bkg,sig)");
47
48 } else {
49
50 // Use object factory to build pdf of tutorial rs502_wspacewrite but
51 // - Contracted to a single line recursive expression,
52 // - Omitting explicit names for components that are not referred to explicitly later
53
54 w->factory("SUM::model(bkgfrac[0.5,0.,1.]*Chebychev::bkg(x[-10,10],{a0[0.5,0.,1],a1[0.2,0.,1.]}),"
55 "SUM(sig1frac[0.8,0.,1.]*Gaussian(x,mean[5,0,10],0.5), Gaussian(x,mean,1)))");
56 }
57
58 // A d v a n c e d p . d . f . c o n s t r u c t o r a r g u m e n t s
59 // ----------------------------------------------------------------
60 //
61 // P.d.f. constructor arguments may by any type of RooAbsArg, but also
62 //
63 // Double_t --> converted to RooConst(...)
64 // {a,b,c} --> converted to RooArgSet() or RooArgList() depending on required ctor arg
65 // dataset name --> converted to RooAbsData reference for any dataset residing in the workspace
66 // enum --> Any enum label that belongs to an enum defined in the (base) class
67
68 // Make a dummy dataset pdf 'model' and import it in the workspace
69 RooDataSet *data = w->pdf("model")->generate(*w->var("x"), 1000);
70 w->import(*data, Rename("data"));
71
72 // Construct a KEYS pdf passing a dataset name and an enum type defining the
73 // mirroring strategy
74 w->factory("KeysPdf::k(x,data,NoMirror,0.2)");
75
76 // Print workspace contents
77 w->Print();
78}
const Bool_t kFALSE
Definition: RtypesCore.h:92
bool Bool_t
Definition: RtypesCore.h:63
RooDataSet * generate(const RooArgSet &whatVars, Int_t nEvents, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none())
See RooAbsPdf::generate(const RooArgSet&,const RooCmdArg&,const RooCmdArg&,const RooCmdArg&,...
Definition: RooAbsPdf.h:58
RooDataSet is a container class to hold unbinned data.
Definition: RooDataSet.h:33
The RooWorkspace is a persistable container for RooFit projects.
Definition: RooWorkspace.h:43
void Print(Option_t *opts=0) const
Print contents of the workspace.
RooRealVar * var(const char *name) const
Retrieve real-valued variable (RooRealVar) with given name. A null pointer is returned if not found.
Bool_t import(const RooAbsArg &arg, const RooCmdArg &arg1=RooCmdArg(), const RooCmdArg &arg2=RooCmdArg(), const RooCmdArg &arg3=RooCmdArg(), const RooCmdArg &arg4=RooCmdArg(), const RooCmdArg &arg5=RooCmdArg(), const RooCmdArg &arg6=RooCmdArg(), const RooCmdArg &arg7=RooCmdArg(), const RooCmdArg &arg8=RooCmdArg(), const RooCmdArg &arg9=RooCmdArg())
Import a RooAbsArg object, e.g.
RooFactoryWSTool & factory()
Return instance to factory tool.
RooAbsPdf * pdf(const char *name) const
Retrieve p.d.f (RooAbsPdf) with given name. A null pointer is returned if not found.
RooCmdArg Rename(const char *suffix)
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...