A++ » ROOSTATS » RooStats::AsymptoticCalculator

class RooStats::AsymptoticCalculator: public RooStats::HypoTestCalculatorGeneric

Function Members (Methods)

public:
virtual~AsymptoticCalculator()
RooStats::AsymptoticCalculatorAsymptoticCalculator(const RooStats::AsymptoticCalculator&)
RooStats::AsymptoticCalculatorAsymptoticCalculator(RooAbsData& data, const RooStats::ModelConfig& altModel, const RooStats::ModelConfig& nullModel, bool nominalAsimov = false)
static TClass*Class()
static RooAbsData*GenerateAsimovData(const RooAbsPdf& pdf, const RooArgSet& observables)
const RooStats::ModelConfig*RooStats::HypoTestCalculatorGeneric::GetAlternateModel() const
const RooArgSet&GetBestFitParams() const
const RooArgSet&GetBestFitPoi() const
const RooAbsData*RooStats::HypoTestCalculatorGeneric::GetData() const
static doubleGetExpectedPValues(double pnull, double palt, double nsigma, bool usecls, bool oneSided = true)
virtual const RooArgSet*RooStats::HypoTestCalculatorGeneric::GetFitInfo() const
virtual RooStats::HypoTestResult*GetHypoTest() const
const RooRealVar*GetMuHat() const
const RooStats::ModelConfig*RooStats::HypoTestCalculatorGeneric::GetNullModel() const
RooStats::TestStatSampler*RooStats::HypoTestCalculatorGeneric::GetTestStatSampler() const
RooStats::HypoTestCalculatorRooStats::HypoTestCalculator::HypoTestCalculator()
RooStats::HypoTestCalculatorRooStats::HypoTestCalculator::HypoTestCalculator(const RooStats::HypoTestCalculator&)
RooStats::HypoTestCalculatorGenericRooStats::HypoTestCalculatorGeneric::HypoTestCalculatorGeneric(const RooStats::HypoTestCalculatorGeneric&)
RooStats::HypoTestCalculatorGenericRooStats::HypoTestCalculatorGeneric::HypoTestCalculatorGeneric(const RooAbsData& data, const RooStats::ModelConfig& altModel, const RooStats::ModelConfig& nullModel, RooStats::TestStatSampler* sampler = 0)
boolInitialize() const
virtual TClass*IsA() const
boolIsOneSidedDiscovery() const
boolIsTwoSided() const
static RooAbsData*MakeAsimovData(const RooStats::ModelConfig& model, const RooArgSet& allParamValues, RooArgSet& globObs)
static RooAbsData*MakeAsimovData(RooAbsData& data, const RooStats::ModelConfig& model, const RooArgSet& poiValues, RooArgSet& globObs, const RooArgSet* genPoiValues = 0)
RooStats::AsymptoticCalculator&operator=(const RooStats::AsymptoticCalculator&)
virtual voidSetAlternateModel(const RooStats::ModelConfig& altModel)
virtual voidRooStats::HypoTestCalculator::SetCommonModel(const RooStats::ModelConfig& model)
virtual voidSetData(RooAbsData& data)
virtual voidSetNullModel(const RooStats::ModelConfig& nullModel)
voidSetOneSided(bool on)
voidSetOneSidedDiscovery(bool on)
static voidSetPrintLevel(int level)
voidSetQTilde(bool on)
voidSetTwoSided()
virtual voidShowMembers(TMemberInspector& insp) const
virtual voidStreamer(TBuffer&)
voidStreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
voidRooStats::HypoTestCalculatorGeneric::UseSameAltToys()
protected:
virtual intRooStats::HypoTestCalculatorGeneric::CheckHook() const
static doubleEvaluateNLL(RooAbsPdf& pdf, RooAbsData& data, const RooArgSet* condObs, const RooArgSet* poiSet = 0)
static voidFillBins(const RooAbsPdf& pdf, const RooArgList& obs, RooAbsData& data, int& index, double& binVolume, int& ibin)
static RooAbsData*GenerateAsimovDataSinglePdf(const RooAbsPdf& pdf, const RooArgSet& obs, const RooRealVar& weightVar, RooCategory* channelCat = 0)
static RooAbsData*GenerateCountingAsimovData(RooAbsPdf& pdf, const RooArgSet& obs, const RooRealVar& weightVar, RooCategory* channelCat = 0)
virtual voidRooStats::HypoTestCalculatorGeneric::PostHook() const
virtual intRooStats::HypoTestCalculatorGeneric::PreAltHook(RooArgSet*, double) const
virtual voidRooStats::HypoTestCalculatorGeneric::PreHook() const
virtual intRooStats::HypoTestCalculatorGeneric::PreNullHook(RooArgSet*, double) const
static boolSetObsToExpected(RooAbsPdf& pdf, const RooArgSet& obs)
static boolSetObsToExpected(RooProdPdf& prod, const RooArgSet& obs)

Data Members

private:
RooAbsData*fAsimovDataasimov data set
RooArgSetfAsimovGlobObssnapshot of Asimov global observables
RooArgSetfBestFitParamssnapshot of all best fitted Parameter values
RooArgSetfBestFitPoisnapshot of best fitted POI values
boolfIsInitialized! flag to check if calculator is initialized
doublefNLLAsimov
doublefNLLObs
boolfNominalAsimovmake Asimov at nominal parameter values
boolfOneSidedfor one sided PL test statistic (upper limits)
boolfOneSidedDiscoveryfor one sided PL test statistic (for discovery)
intfUseQTildeflag to indicate if using qtilde or not (-1 (default based on RooRealVar)), 0 false, 1 (true)
static intfgPrintLevelcontrol print level (0 minimal, 1 normal, 2 debug)

Class Charts

Inheritance Chart:
RooStats::HypoTestCalculator
RooStats::HypoTestCalculatorGeneric
RooStats::AsymptoticCalculator

Function documentation

AsymptoticCalculator(RooAbsData& data, const RooStats::ModelConfig& altModel, const RooStats::ModelConfig& nullModel, bool nominalAsimov = false)
    HypoTestCalculatorGeneric(data, altModel, nullModel, 0)
 {
 }
bool Initialize() const
 initialize the calculator by performin g a global fit and make the Asimov data set
HypoTestResult * GetHypoTest() const
 re-implement HypoTest computation using the asymptotic
RooAbsData * MakeAsimovData(RooAbsData& data, const RooStats::ModelConfig& model, const RooArgSet& poiValues, RooArgSet& globObs, const RooArgSet* genPoiValues = 0)
 make the asimov data from the ModelConfig and list of poi - return data set annd snapshoot of global obs
 poiValues is the snapshot of POI used for finding the best buisance parameter values (conditioned at these values)
 genPoiValues is optionally a different set of POI values used for generating. By default the same POI are used for generating and for finding the nuisance parameters
RooAbsData * MakeAsimovData(const RooStats::ModelConfig& model, const RooArgSet& allParamValues, RooArgSet& globObs)
 make a nominal asimov data from the ModelConfig and parameter values
 The parameter values (including the nunisance) could be given from a fit to data or be at the nominal values
RooAbsData * GenerateAsimovData(const RooAbsPdf& pdf, const RooArgSet& observables)
double GetExpectedPValues(double pnull, double palt, double nsigma, bool usecls, bool oneSided = true)
 function given the null and the alt p value - return the expected one given the N - sigma value
void SetOneSided(bool on)
 set test statistic for one sided (upper limits)
{ fOneSided = on; }
void SetTwoSided()
 set the test statistics for two sided (in case of upper limits
 for discovery does not make really sense)
{ fOneSided = false; fOneSidedDiscovery = false;}
void SetOneSidedDiscovery(bool on)
 set the test statistics for one-sided discovery
void SetNullModel(const RooStats::ModelConfig& nullModel)
 re-implementation of  setters since they are needed to re-initialize the calculator
void SetAlternateModel(const RooStats::ModelConfig& altModel)
void SetData(RooAbsData& data)
bool IsTwoSided() const
{ return (!fOneSided && !fOneSidedDiscovery); }
bool IsOneSidedDiscovery() const
{ return fOneSidedDiscovery; }
void SetQTilde(bool on)
 set using of qtilde, by default is controlled if RoORealVar is limited or not
{ fUseQTilde = on; }
const RooArgSet & GetBestFitPoi() const
 return snapshot of the best fit parameter
{ return fBestFitPoi; }
const RooRealVar * GetMuHat() const
 return best fit parameter (firs of poi)
{ return dynamic_cast<RooRealVar*>(fBestFitPoi.first()); }
const RooArgSet & GetBestFitParams() const
 return best fit value for all parameters
{ return fBestFitPoi; }
void SetPrintLevel(int level)
RooAbsData * GenerateAsimovDataSinglePdf(const RooAbsPdf& pdf, const RooArgSet& obs, const RooRealVar& weightVar, RooCategory* channelCat = 0)
 // configure TestStatSampler for the Null run
 int PreNullHook(RooArgSet *parameterPoint, double obsTestStat) const;
 // configure TestStatSampler for the Alt run
 int PreAltHook(RooArgSet *parameterPoint, double obsTestStat) const;
void FillBins(const RooAbsPdf& pdf, const RooArgList& obs, RooAbsData& data, int& index, double& binVolume, int& ibin)
bool SetObsToExpected(RooAbsPdf& pdf, const RooArgSet& obs)
bool SetObsToExpected(RooProdPdf& prod, const RooArgSet& obs)