32 void rf607_fitresult()
41 RooRealVar mean(
"mean",
"mean of gaussians",5,-10,10) ;
42 RooRealVar sigma1(
"sigma1",
"width of gaussians",0.5,0.1,10) ;
43 RooRealVar sigma2(
"sigma2",
"width of gaussians",1,0.1,10) ;
45 RooGaussian sig1(
"sig1",
"Signal component 1",
x,mean,sigma1) ;
46 RooGaussian sig2(
"sig2",
"Signal component 2",
x,mean,sigma2) ;
54 RooRealVar sig1frac(
"sig1frac",
"fraction of component 1 in signal",0.8,0.,1.) ;
58 RooRealVar bkgfrac(
"bkgfrac",
"fraction of background",0.5,0.,1.) ;
95 RooPlot* frame =
new RooPlot(sigma1,sig1frac,0.45,0.60,0.65,0.90) ;
96 frame->
SetTitle(
"Covariance between sigma1 and sig1frac") ;
97 r->
plotOn(frame,sigma1,sig1frac,
"ME12ABHV") ;
105 cout <<
"EDM = " << r->
edm() << endl ;
106 cout <<
"-log(L) at minimum = " << r->
minNll() << endl ;
109 cout <<
"final value of floating parameters" << endl ;
113 cout <<
"correlation between sig1frac and a0 is " << r->
correlation(sig1frac,a0) << endl ;
114 cout <<
"correlation between bkgfrac and mean is " << r->
correlation(
"bkgfrac",
"mean") << endl ;
121 cout <<
"correlation matrix" << endl ;
123 cout <<
"covariance matrix" << endl ;
131 TFile
f(
"rf607_fitresult.root",
"RECREATE") ;
139 TCanvas*
c =
new TCanvas(
"rf607_fitresult",
"rf607_fitresult",800,400) ;
const TMatrixDSym & correlationMatrix() const
Return correlation matrix ;.
virtual void SetTitleOffset(Float_t offset=1)
Set distance between the axis and the axis title Offset is a correction factor with respect to the "s...
virtual Int_t Write(const char *name=0, Int_t option=0, Int_t bufsize=0)
Write this object to the current directory.
RooPlot * plotOn(RooPlot *frame, const RooAbsArg &par1, const RooAbsArg &par2, const char *options="ME") const
RooAddPdf is an efficient implementation of a sum of PDFs of the form.
const RooArgList & floatParsFinal() const
R__EXTERN TStyle * gStyle
const TMatrixDSym & covarianceMatrix() const
Return covariance matrix.
TVirtualPad * cd(Int_t subpadnumber=0)
Set current canvas & pad.
void SetTitle(const char *name)
Set the title of the RooPlot to 'title'.
virtual void Print(Option_t *options=0) const
Print TNamed name and title.
virtual void Print(Option_t *options=0) const
This method must be overridden when a class wants to print itself.
RooRealVar represents a fundamental (non-derived) real valued object.
TH2 * correlationHist(const char *name="correlation_matrix") const
Return TH2D of correlation matrix.
Service class for 2-Dim histogram classes.
virtual void Draw(Option_t *option="")
Draw this histogram with options.
RooDataSet is a container class to hold unbinned data.
A RooPlot is a plot frame and a container for graphics objects within that frame. ...
void Print(Option_t *name="") const
Print the matrix as a table of elements.
RooCmdArg Save(Bool_t flag=kTRUE)
virtual void Divide(Int_t nx=1, Int_t ny=1, Float_t xmargin=0.01, Float_t ymargin=0.01, Int_t color=0)
Automatic pad generation by division.
void SetOptStat(Int_t stat=1)
The type of information printed in the histogram statistics box can be selected via the parameter mod...
Chebychev polynomial p.d.f.
Double_t correlation(const RooAbsArg &par1, const RooAbsArg &par2) const
virtual void Draw(Option_t *options=0)
Draw this plot and all of the elements it contains.