1 |
#ifndef GUARD_StandardHypoTestInvDemo_h
|
2 |
#define GUARD_StandardHypoTestInvDemo_h
|
3 |
|
4 |
|
5 |
#include "TFile.h"
|
6 |
#include "RooWorkspace.h"
|
7 |
#include "RooAbsPdf.h"
|
8 |
#include "RooRealVar.h"
|
9 |
#include "RooDataSet.h"
|
10 |
#include "RooStats/ModelConfig.h"
|
11 |
#include "TGraphErrors.h"
|
12 |
#include "TGraphAsymmErrors.h"
|
13 |
#include "TCanvas.h"
|
14 |
#include "TLine.h"
|
15 |
|
16 |
#include "RooStats/HybridCalculator.h"
|
17 |
#include "RooStats/FrequentistCalculator.h"
|
18 |
#include "RooStats/ToyMCSampler.h"
|
19 |
#include "RooStats/HypoTestPlot.h"
|
20 |
|
21 |
#include "RooStats/NumEventsTestStat.h"
|
22 |
#include "RooStats/ProfileLikelihoodTestStat.h"
|
23 |
#include "RooStats/SimpleLikelihoodRatioTestStat.h"
|
24 |
#include "RooStats/RatioOfProfiledLikelihoodsTestStat.h"
|
25 |
#include "RooStats/MaxLikelihoodEstimateTestStat.h"
|
26 |
|
27 |
#include "RooStats/HypoTestInverter.h"
|
28 |
#include "RooStats/HypoTestInverterResult.h"
|
29 |
#include "RooStats/HypoTestInverterPlot.h"
|
30 |
|
31 |
// internal routine to run the inverter
|
32 |
RooStats::HypoTestInverterResult* RunInverter(RooWorkspace * w, const char * modelSBName, const char * modelBName, const char * dataName,
|
33 |
int type, int testStatType, int npoints, double poimin, double poimax, int ntoys,
|
34 |
bool useCls, bool useNumberCounting, const char * nuisPriorName);
|
35 |
|
36 |
void StandardHypoTestInvDemo(const char * fileName =0,
|
37 |
const char * wsName = "combined",
|
38 |
const char * modelSBName = "ModelConfig",
|
39 |
const char * modelBName = "",
|
40 |
const char * dataName = "obsData",
|
41 |
int calculatorType = 0,
|
42 |
int testStatType = 3,
|
43 |
bool useCls = true ,
|
44 |
int npoints = 5,
|
45 |
double poimin = 0,
|
46 |
double poimax = 5,
|
47 |
int ntoys=1000,
|
48 |
bool useNumberCounting = false,
|
49 |
const char * nuisPriorName = 0);
|
50 |
void ReadResult(const char * fileName, const char * resultName, bool useCLs);
|
51 |
|
52 |
|
53 |
#endif
|