Multidimensional models: conditional pdf with per-event errors
{
RooRealVar dterr(
"dterr",
"per-event error on dt", 0.01, 10);
RooLandau pdfDtErr(
"pdfDtErr",
"pdfDtErr", dterr, RooConst(1), RooConst(0.25));
RooDataSet *expDataDterr = pdfDtErr.generate(dterr, 10000);
TH1 *hh_decay = decay_gm.createHistogram(
"hh_decay", dt, Binning(50), YVar(dterr, Binning(50)));
RooPlot *frame = dt.
frame(Title(
"Slices of decay(dt|dterr) at various dterr"));
for (
Int_t ibin = 0; ibin < 100; ibin += 20) {
dterr.setBin(ibin);
decay_gm.plotOn(frame, Normalization(5.));
}
RooPlot *frame2 = dt.
frame(Title(
"Projection of decay(dt|dterr) on dt"));
decay_gm.plotOn(frame2, ProjWData(*expDataDterr, true));
TCanvas *
c =
new TCanvas(
"rf306_condpereventerrors",
"rf306_condperventerrors", 1200, 400);
gPad->SetLeftMargin(0.20);
gPad->SetLeftMargin(0.15);
gPad->SetLeftMargin(0.15);
}
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void data
RooDataSet is a container class to hold unbinned data.
Single or double sided decay function that can be analytically convolved with any RooResolutionModel ...
Class RooGaussModel implements a RooResolutionModel that models a Gaussian distribution.
Landau distribution p.d.f.
A RooPlot is a plot frame and a container for graphics objects within that frame.
static RooPlot * frame(const RooAbsRealLValue &var, double xmin, double xmax, Int_t nBins)
Create a new frame for a given variable in x.
void Draw(Option_t *options=nullptr) override
Draw this plot and all of the elements it contains.
RooRealVar represents a variable that can be changed from the outside.
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title.
virtual void SetLineColor(Color_t lcolor)
Set the line color.
TH1 is the base class of all histogram classes in ROOT.
void Draw(Option_t *option="") override
Draw this histogram with options.
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
[#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization
[#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization
[#1] INFO:NumericIntegration -- RooRealIntegral::init(gm1_conv_exp(-abs(@0)/@1)_dt_tau_[decay_gm]_Int[dt,dterr]) using numeric integrator RooIntegrator1D to calculate Int(dterr)
[#1] INFO:Plotting -- RooAbsReal::plotOn(decay_gm) plot on dt averages using data variables (dterr)
[#1] INFO:Plotting -- RooDataWeightedAverage::ctor(decay_gmDataWgtAvg) constructing data weighted average of function decay_gm_Norm[dt] over 100 data points of (dterr) with a total weight of 10000
.........................................................................................................................................................................................................................
- Date
- July 2008
- Author
- Wouter Verkerke
Definition in file rf306_condpereventerrors.C.