doc/v632/rf304__uncorrprod_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussian.h"

#include "RooProdPdf.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "RooPlot.h"

using namespace RooFit;


void rf304_uncorrprod()

{


   // C r e a t e   c o m p o n e n t   p d f s   i n   x   a n d   y

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


   // Create two pdfs gaussx(x,meanx,sigmax) gaussy(y,meany,sigmay) and its variables

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

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


   RooRealVar meanx("mean1", "mean of gaussian x", 2);

   RooRealVar meany("mean2", "mean of gaussian y", -2);

   RooRealVar sigmax("sigmax", "width of gaussian x", 1);

   RooRealVar sigmay("sigmay", "width of gaussian y", 5);


   RooGaussian gaussx("gaussx", "gaussian PDF", x, meanx, sigmax);

   RooGaussian gaussy("gaussy", "gaussian PDF", y, meany, sigmay);


   // C o n s t r u c t   u n c o r r e l a t e d   p r o d u c t   p d f

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


   // Multiply gaussx and gaussy into a two-dimensional pdf gaussxy

   RooProdPdf gaussxy("gaussxy", "gaussx*gaussy", RooArgList(gaussx, gaussy));


   // S a m p l e   p d f ,   p l o t   p r o j e c t i o n   o n   x   a n d   y

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


   // Generate 10000 events in x and y from gaussxy

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


   // Plot x distribution of data and projection of gaussxy on x = Int(dy) gaussxy(x,y)

   RooPlot *xframe = x.frame(Title("X projection of gauss(x)*gauss(y)"));

   data->plotOn(xframe);

   gaussxy.plotOn(xframe);


   // Plot x distribution of data and projection of gaussxy on y = Int(dx) gaussxy(x,y)

   RooPlot *yframe = y.frame(Title("Y projection of gauss(x)*gauss(y)"));

   data->plotOn(yframe);

   gaussxy.plotOn(yframe);


   // Make canvas and draw RooPlots

   TCanvas *c = new TCanvas("rf304_uncorrprod", "rf304_uncorrprod", 800, 400);

   c->Divide(2);

   c->cd(1);

   gPad->SetLeftMargin(0.15);

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

   xframe->Draw();

   c->cd(2);

   gPad->SetLeftMargin(0.15);

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

   yframe->Draw();

}

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

RooDataSet.h

RooGaussian.h

RooPlot.h

RooProdPdf.h

RooRealVar.h

TAxis.h

TCanvas.h

TRangeDynCast
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Definition TCollection.h:358

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

gPad
#define gPad
Definition TVirtualPad.h:305

ROOT::Detail::TRangeCast
Definition TCollection.h:311

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

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

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

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

TCanvas
The Canvas class.
Definition TCanvas.h:23

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

rf304_uncorrprod
Definition rf304_uncorrprod.py:1

xmlio::Title
const char * Title
Definition TXMLSetup.cxx:68