doc/v630/rf312__multirangefit_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussian.h"

#include "RooProdPdf.h"

#include "RooAddPdf.h"

#include "RooPolynomial.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "RooPlot.h"

#include "RooFitResult.h"

using namespace RooFit;


void rf312_multirangefit()

{


   // C r e a t e   2 D   p d f   a n d   d a t a

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


   // Define observables x,y

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

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


   // Construct the signal pdf gauss(x)*gauss(y)

   RooRealVar mx("mx", "mx", 1, -10, 10);

   RooRealVar my("my", "my", 1, -10, 10);


   RooGaussian gx("gx", "gx", x, mx, 1.0);

   RooGaussian gy("gy", "gy", y, my, 1.0);


   RooProdPdf sig("sig", "sig", gx, gy);


   // Construct the background pdf (flat in x,y)

   RooPolynomial px("px", "px", x);

   RooPolynomial py("py", "py", y);

   RooProdPdf bkg("bkg", "bkg", px, py);


   // Construct the composite model sig+bkg

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

   RooAddPdf model("model", "model", RooArgList(sig, bkg), f);


   // Sample 10000 events in (x,y) from the model

   std::unique_ptr<RooDataSet> modelData{model.generate({x, y}, 10000)};


   // D e f i n e   s i g n a l   a n d   s i d e b a n d   r e g i o n s

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


   // Construct the SideBand1,SideBand2,Signal regions

   //

   //                    |

   //      +-------------+-----------+

   //      |             |           |

   //      |    Side     |   Sig     |

   //      |    Band1    |   nal     |

   //      |             |           |

   //    --+-------------+-----------+--

   //      |                         |

   //      |           Side          |

   //      |           Band2         |

   //      |                         |

   //      +-------------+-----------+

   //                    |


   x.setRange("SB1", -10, +10);

   y.setRange("SB1", -10, 0);


   x.setRange("SB2", -10, 0);

   y.setRange("SB2", 0, +10);


   x.setRange("SIG", 0, +10);

   y.setRange("SIG", 0, +10);


   x.setRange("FULL", -10, +10);

   y.setRange("FULL", -10, +10);


   // P e r f o r m   f i t s   i n   i n d i v i d u a l   s i d e b a n d   r e g i o n s

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


   // Perform fit in SideBand1 region (RooAddPdf coefficients will be interpreted in full range)

   std::unique_ptr<RooFitResult> r_sb1{model.fitTo(*modelData, Range("SB1"), Save(), PrintLevel(-1))};


   // Perform fit in SideBand2 region (RooAddPdf coefficients will be interpreted in full range)

   std::unique_ptr<RooFitResult> r_sb2{model.fitTo(*modelData, Range("SB2"), Save(), PrintLevel(-1))};


   // P e r f o r m   f i t s   i n   j o i n t    s i d e b a n d   r e g i o n s

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


   // Now perform fit to joint 'L-shaped' sideband region 'SB1|SB2'

   // (RooAddPdf coefficients will be interpreted in full range)

   std::unique_ptr<RooFitResult> r_sb12{model.fitTo(*modelData, Range("SB1,SB2"), Save(), PrintLevel(-1))};


   // Print results for comparison

   r_sb1->Print();

   r_sb2->Print();

   r_sb12->Print();

}

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

RooAddPdf.h

RooDataSet.h

RooFitResult.h

RooGaussian.h

RooPlot.h

RooPolynomial.h

RooProdPdf.h

RooRealVar.h

TAxis.h

TCanvas.h

mx
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t Int_t Int_t Window_t TString Int_t GCValues_t GetPrimarySelectionOwner GetDisplay GetScreen GetColormap GetNativeEvent const char const char dpyName wid window const char font_name cursor keysym reg const char only_if_exist regb h Point_t winding char text const char depth char const char mx
Definition TGWin32VirtualXProxy.cxx:228

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

RooArgList
RooArgList is a container object that can hold multiple RooAbsArg objects.
Definition RooArgList.h:22

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

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

RooRealVar
RooRealVar represents a variable that can be changed from the outside.
Definition RooRealVar.h:37

RooFit::Save
RooCmdArg Save(bool flag=true)
Definition RooGlobalFunc.cxx:564

RooFit::PrintLevel
RooCmdArg PrintLevel(Int_t code)
Definition RooGlobalFunc.cxx:572

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

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

rf312_multirangefit
Definition rf312_multirangefit.py:1

Range
Ta Range(0, 0, 1, 1)