doc/v636/rf310__sliceplot_8C.html

#include "RooRealVar.h"

#include "RooDataSet.h"

#include "RooGaussModel.h"

#include "RooDecay.h"

#include "RooBMixDecay.h"

#include "RooCategory.h"

#include "TCanvas.h"

#include "TAxis.h"

#include "RooPlot.h"

using namespace RooFit;


void rf310_sliceplot()

{


   // C r e a t e   B   d e c a y   p d f   w it h   m i x i n g

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


   // Decay time observables

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


   // Discrete observables mixState (B0tag==B0reco?) and tagFlav (B0tag==B0(bar)?)

   RooCategory mixState("mixState", "B0/B0bar mixing state");

   RooCategory tagFlav("tagFlav", "Flavour of the tagged B0");


   // Define state labels of discrete observables

   mixState.defineType("mixed", -1);

   mixState.defineType("unmixed", 1);

   tagFlav.defineType("B0", 1);

   tagFlav.defineType("B0bar", -1);


   // Model parameters

   RooRealVar dm("dm", "delta m(B)", 0.472, 0., 1.0);

   RooRealVar tau("tau", "B0 decay time", 1.547, 1.0, 2.0);

   RooRealVar w("w", "Flavor Mistag rate", 0.03, 0.0, 1.0);

   RooRealVar dw("dw", "Flavor Mistag rate difference between B0 and B0bar", 0.01);


   // Build a gaussian resolution model

   RooRealVar bias1("bias1", "bias1", 0);

   RooRealVar sigma1("sigma1", "sigma1", 0.01);

   RooGaussModel gm1("gm1", "gauss model 1", dt, bias1, sigma1);


   // Construct a decay pdf, smeared with single gaussian resolution model

   RooBMixDecay bmix_gm1("bmix", "decay", dt, mixState, tagFlav, tau, dm, w, dw, gm1, RooBMixDecay::DoubleSided);


   // Generate BMixing data with above set of event errors

   std::unique_ptr<RooDataSet> data{bmix_gm1.generate({dt, tagFlav, mixState}, 20000)};


   // P l o t   f u l l   d e c a y   d i s t r i b u t i o n

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


   // Create frame, plot data and pdf projection (integrated over tagFlav and mixState)

   RooPlot *frame = dt.frame(Title("Inclusive decay distribution"));

   data->plotOn(frame);

   bmix_gm1.plotOn(frame);


   // P l o t   d e c a y   d i s t r .   f o r   m i x e d   a n d   u n m i x e d   s l i c e   o f   m i x S t a t e

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


   // Create frame, plot data (mixed only)

   RooPlot *frame2 = dt.frame(Title("Decay distribution of mixed events"));

   data->plotOn(frame2, Cut("mixState==mixState::mixed"));


   // Position slice in mixState at "mixed" and plot slice of pdf in mixstate over data (integrated over tagFlav)

   bmix_gm1.plotOn(frame2, Slice(mixState, "mixed"));


   // Create frame, plot data (unmixed only)

   RooPlot *frame3 = dt.frame(Title("Decay distribution of unmixed events"));

   data->plotOn(frame3, Cut("mixState==mixState::unmixed"));


   // Position slice in mixState at "unmixed" and plot slice of pdf in mixstate over data (integrated over tagFlav)

   bmix_gm1.plotOn(frame3, Slice(mixState, "unmixed"));


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

   c->Divide(3);

   c->cd(1);

   gPad->SetLeftMargin(0.15);

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

   gPad->SetLogy();

   frame->Draw();

   c->cd(2);

   gPad->SetLeftMargin(0.15);

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

   gPad->SetLogy();

   frame2->Draw();

   c->cd(3);

   gPad->SetLeftMargin(0.15);

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

   gPad->SetLogy();

   frame3->Draw();

}

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

RooBMixDecay.h

RooCategory.h

RooDataSet.h

RooDecay.h

RooGaussModel.h

RooPlot.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

w
winID w
Definition TGWin32VirtualGLProxy.cxx:39

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

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

RooBMixDecay
Class RooBMixDecay is a RooAbsAnaConvPdf implementation that describes the decay of B mesons with the...
Definition RooBMixDecay.h:23

RooBMixDecay::DoubleSided
@ DoubleSided
Definition RooBMixDecay.h:26

RooCategory
Object to represent discrete states.
Definition RooCategory.h:28

RooGaussModel
Class RooGaussModel implements a RooResolutionModel that models a Gaussian distribution.
Definition RooGaussModel.h:25

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

RooPlot::GetYaxis
TAxis * GetYaxis() const
Definition RooPlot.cxx:1223

RooPlot::Draw
void Draw(Option_t *options=nullptr) override
Draw this plot and all of the elements it contains.
Definition RooPlot.cxx:596

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

TAttAxis::SetTitleOffset
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
Definition TAttAxis.cxx:280

TCanvas
The Canvas class.
Definition TCanvas.h:23

RooFit::Cut
RooCmdArg Cut(const char *cutSpec)
Definition RooGlobalFunc.cxx:300

ROOT::Internal::Slice
T Slice(const T &histo, std::vector< Int_t > &args)
Creates a sliced copy of the given histogram.
Definition TH1.h:74

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

rf310_sliceplot
Definition rf310_sliceplot.py:1

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