Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
TestNonCentral.C
Go to the documentation of this file.
1/// \file
2/// \ingroup tutorial_roostats
3/// \notebook -js
4///
5/// \macro_image
6/// \macro_output
7/// \macro_code
8///
9/// \author Lorenzo Moneta
10
11using namespace RooFit;
12void TestNonCentral()
13{
14
15 RooWorkspace w("w");
16 // k <2, must use sum
17 w.factory("NonCentralChiSquare::nc(x[0,50],k[1.99,0,5],lambda[5])");
18 // kk > 2 can use bessel
19 w.factory("NonCentralChiSquare::ncc(x,kk[2.01,0,5],lambda)");
20 // kk > 2, force sum
21 w.factory("NonCentralChiSquare::nccc(x,kk,lambda)");
22 ((RooNonCentralChiSquare *)w.pdf("nccc"))->SetForceSum(true);
23
24 // a normal "central" chi-square for comparison when lambda->0
25 w.factory("ChiSquarePdf::cs(x,k)");
26
27 // w.var("kk")->setVal(4.); // test a large kk
28
29 std::unique_ptr<RooDataSet> ncdata{w.pdf("nc")->generate(*w.var("x"), 100)};
30 std::unique_ptr<RooDataSet> csdata{w.pdf("cs")->generate(*w.var("x"), 100)};
31 RooPlot *plot = w.var("x")->frame();
32 ncdata->plotOn(plot, MarkerColor(kRed));
33 csdata->plotOn(plot, MarkerColor(kBlue));
34 w.pdf("nc")->plotOn(plot, LineColor(kRed));
35 w.pdf("ncc")->plotOn(plot, LineColor(kGreen));
36 w.pdf("nccc")->plotOn(plot, LineColor(kYellow), LineStyle(kDashed));
37 w.pdf("cs")->plotOn(plot, LineColor(kBlue), LineStyle(kDotted));
38 plot->Draw();
39}
@ kRed
Definition Rtypes.h:66
@ kGreen
Definition Rtypes.h:66
@ kBlue
Definition Rtypes.h:66
@ kYellow
Definition Rtypes.h:66
@ kDashed
Definition TAttLine.h:48
@ kDotted
Definition TAttLine.h:48
winID h TVirtualViewer3D TVirtualGLPainter char TVirtualGLPainter plot
The PDF of the Non-Central Chi Square distribution for n degrees of freedom.
Plot frame and a container for graphics objects within that frame.
Definition RooPlot.h:43
Persistable container for RooFit projects.
virtual void Draw(Option_t *option="")
Default Draw method for all objects.
Definition TObject.cxx:274
RooCmdArg MarkerColor(Color_t color)
RooCmdArg LineColor(Color_t color)
RooCmdArg LineStyle(Style_t style)
The namespace RooFit contains mostly switches that change the behaviour of functions of PDFs (or othe...
Definition JSONIO.h:26