doc/master/rf705__linearmorph_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussian.h"

#include "RooPolynomial.h"

#include "RooIntegralMorph.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "RooPlot.h"

#include "TH1.h"

using namespace RooFit;


void rf705_linearmorph()

{

   // C r e a t e   e n d   p o i n t   p d f   s h a p e s

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


   // Observable

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


   // Lower end point shape: a Gaussian

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

   RooGaussian g1("g1", "g1", x, g1mean, 2.0);


   // Upper end point shape: a Polynomial

   RooPolynomial g2("g2", "g2", x, RooArgSet(-0.03, -0.001));


   // C r e a t e   i n t e r p o l a t i n g   p d f

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


   // Create interpolation variable

   RooRealVar alpha("alpha", "alpha", 0, 1.0);


   // Specify sampling density on observable and interpolation variable

   x.setBins(1000, "cache");

   alpha.setBins(50, "cache");


   // Construct interpolating pdf in (x,a) represent g1(x) at a=a_min

   // and g2(x) at a=a_max

   RooIntegralMorph lmorph("lmorph", "lmorph", g1, g2, x, alpha);


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

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


   // Show end points as blue curves

   RooPlot *frame1 = x.frame();

   g1.plotOn(frame1);

   g2.plotOn(frame1);


   // Show interpolated shapes in red

   alpha.setVal(0.125);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.25);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.375);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.50);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.625);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.75);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.875);

   lmorph.plotOn(frame1, LineColor(kRed));

   alpha.setVal(0.95);

   lmorph.plotOn(frame1, LineColor(kRed));


   // S h o w   2 D   d i s t r i b u t i o n   o f   p d f ( x , a l p h a )

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


   // Create 2D histogram

   TH1 *hh = lmorph.createHistogram("hh", x, Binning(40), YVar(alpha, Binning(40)));

   hh->SetLineColor(kBlue);


   // F i t   p d f   t o   d a t a s e t   w i t h   a l p h a = 0 . 8

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


   // Generate a toy dataset at alpha = 0.8

   alpha = 0.8;

   std::unique_ptr<RooDataSet> data{lmorph.generate(x, 1000)};


   // Fit pdf to toy data

   lmorph.setCacheAlpha(true);

   lmorph.fitTo(*data, Verbose(true), PrintLevel(-1));


   // Plot fitted pdf and data overlaid

   RooPlot *frame2 = x.frame(Bins(100));

   data->plotOn(frame2);

   lmorph.plotOn(frame2);


   // S c a n   - l o g ( L )   v s   a l p h a

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


   // Show scan -log(L) of dataset w.r.t alpha

   RooPlot *frame3 = alpha.frame(Bins(100), Range(0.1, 0.9));


   // Make 2D pdf of histogram

   std::unique_ptr<RooAbsReal> nll{lmorph.createNLL(*data)};

   nll->plotOn(frame3, ShiftToZero());


   lmorph.setCacheAlpha(false);


   TCanvas *c = new TCanvas("rf705_linearmorph", "rf705_linearmorph", 800, 800);

   c->Divide(2, 2);

   c->cd(1);

   gPad->SetLeftMargin(0.15);

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

   frame1->Draw();

   c->cd(2);

   gPad->SetLeftMargin(0.20);

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

   hh->Draw("surf");

   c->cd(3);

   gPad->SetLeftMargin(0.15);

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

   frame3->Draw();

   c->cd(4);

   gPad->SetLeftMargin(0.15);

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

   frame2->Draw();


   return;

}

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

RooDataSet.h

RooGaussian.h

RooIntegralMorph.h

RooPlot.h

RooPolynomial.h

RooRealVar.h

kRed
@ kRed
Definition Rtypes.h:67

kBlue
@ kBlue
Definition Rtypes.h:67

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

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

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

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

RooIntegralMorph
Class RooIntegralMorph is an implementation of the histogram interpolation technique described by Ale...
Definition RooIntegralMorph.h:26

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

RooPolynomial
RooPolynomial implements a polynomial p.d.f of the form.
Definition RooPolynomial.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:109

RooFit::ShiftToZero
RooCmdArg ShiftToZero()
Definition RooGlobalFunc.cxx:205

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

RooFit::Detail::MathFuncs::nll
double nll(double pdf, double weight, int binnedL, int doBinOffset)
Definition MathFuncs.h:452

RooFit
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition CodegenImpl.h:72

rf705_linearmorph
Definition rf705_linearmorph.py:1

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