doc/v636/rf301__composition_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussian.h"

#include "RooPolyVar.h"

#include "RooPlot.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "TH1.h"

using namespace RooFit;


void rf301_composition()

{

   // S e t u p   c o m p o s e d   m o d e l   g a u s s ( x , m ( y ) , s )

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


   // Create observables

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

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


   // Create function f(y) = a0 + a1*y

   RooRealVar a0("a0", "a0", -0.5, -5, 5);

   RooRealVar a1("a1", "a1", -0.5, -1, 1);

   RooPolyVar fy("fy", "fy", y, RooArgSet(a0, a1));


   // Create gauss(x,f(y),s)

   RooRealVar sigma("sigma", "width of gaussian", 0.5);

   RooGaussian model("model", "Gaussian with shifting mean", x, fy, sigma);


   // S a m p l e   d a t a ,   p l o t   d a t a   a n d   p d f   o n   x   a n d   y

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


   // Generate 10000 events in x and y from model

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


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

   RooPlot *xframe = x.frame();

   data->plotOn(xframe);

   model.plotOn(xframe);


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

   RooPlot *yframe = y.frame();

   data->plotOn(yframe);

   model.plotOn(yframe);


   // Make two-dimensional plot in x vs y

   TH1 *hh_model = model.createHistogram("hh_model", x, Binning(50), YVar(y, Binning(50)));

   hh_model->SetLineColor(kBlue);


   // Make canvas and draw RooPlots

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

   c->Divide(3);

   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->cd(3);

   gPad->SetLeftMargin(0.20);

   hh_model->GetZaxis()->SetTitleOffset(2.5);

   hh_model->Draw("surf");

}

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

RooDataSet.h

RooGaussian.h

RooPlot.h

RooPolyVar.h

RooRealVar.h

kBlue
@ kBlue
Definition Rtypes.h:66

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

TH1.h

gPad
#define gPad
Definition TVirtualPad.h:308

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

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

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

RooPolyVar
A RooAbsReal implementing a polynomial in terms of a list of RooAbsReal coefficients.
Definition RooPolyVar.h:25

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

TCanvas
The Canvas class.
Definition TCanvas.h:23

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

RooFit::YVar
RooCmdArg YVar(const RooAbsRealLValue &var, const RooCmdArg &arg={})
Definition RooGlobalFunc.cxx:826

RooFit::Binning
RooCmdArg Binning(const RooAbsBinning &binning)
Definition RooGlobalFunc.cxx:308

sigma
const Double_t sigma
Definition h1analysisProxy.h:11

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

rf301_composition
Definition rf301_composition.py:1