14 #ifndef __CINT__ // problem including this file with CINT
32 using namespace RooFit;
33 using namespace RooStats;
54 RooExtendPdf bkg_ext_pdf(
"bkg_ext_pdf",
"",bkg_pdf,bkg_yield);
78 &nuisance_parameters, &bkg_yield_prior);
93 if (! myHybridResult) {
94 std::cerr <<
"\nError returned from Hypothesis test" << std::endl;
119 HybridPlot* myHybridPlot = myHybridResult->GetPlot(
"myHybridPlot",
"Plot of results with HybridCalculatorOriginal",100);
120 myHybridPlot->
Draw();
123 double clsb_data = myHybridResult->CLsplusb();
124 double clb_data = myHybridResult->CLb();
125 double cls_data = myHybridResult->CLs();
126 double data_significance = myHybridResult->Significance();
127 double min2lnQ_data = myHybridResult->GetTestStat_data();
130 double mean_sb_toys_test_stat = myHybridPlot->GetSBmean();
131 myHybridResult->SetDataTestStatistics(mean_sb_toys_test_stat);
132 double toys_significance = myHybridResult->Significance();
134 std::cout <<
"Completed HybridCalculatorOriginal example:\n";
135 std::cout <<
" - -2lnQ = " << min2lnQ_data << endl;
136 std::cout <<
" - CL_sb = " << clsb_data << std::endl;
137 std::cout <<
" - CL_b = " << clb_data << std::endl;
138 std::cout <<
" - CL_s = " << cls_data << std::endl;
139 std::cout <<
" - significance of data = " << data_significance << std::endl;
140 std::cout <<
" - mean significance of toys = " << toys_significance << std::endl;
HybridCalculatorOriginal class.
void Draw(const char *options="")
Draw on current pad.
virtual void SetSeed(UInt_t seed=0)
Set the random generator seed.
RooCmdArg Extended(Bool_t flag=kTRUE)
void HybridOriginalDemo(int ntoys=1000)
void SetNumberOfToys(unsigned int ntoys)
static TRandom * randomGenerator()
Return a pointer to a singleton random-number generator implementation.
Double_t getVal(const RooArgSet *set=0) const
This class provides the plots for the result of a study performed with the HybridCalculatorOriginal c...
virtual HybridResult * GetHypoTest() const
inherited methods from HypoTestCalculator interface
RooRealVar represents a fundamental (non-derived) real valued object.
void SetGenerateBinned(bool on=true)
void setConstant(Bool_t value=kTRUE)
void UseNuisance(bool on=true)
RooDataSet is a container class to hold unbinned data.
RooConstVar & RooConst(Double_t val)
RooDataSet * generate(const RooArgSet &whatVars, Int_t nEvents, const RooCmdArg &arg1, const RooCmdArg &arg2=RooCmdArg::none(), const RooCmdArg &arg3=RooCmdArg::none(), const RooCmdArg &arg4=RooCmdArg::none(), const RooCmdArg &arg5=RooCmdArg::none())
Generate a new dataset containing the specified variables with events sampled from our distribution...
RooPolynomial implements a polynomial p.d.f of the form By default coefficient a_0 is chosen to be 1...
void SetTestStatistic(int index)
set the desired test statistics: index=1 : 2 * log( L_sb / L_b ) (DEFAULT) index=2 : number of genera...