1 #ifndef HiggsAnalysis_CombinedLimit_SimplerLikelihoodRatioTestStat_h
2 #define HiggsAnalysis_CombinedLimit_SimplerLikelihoodRatioTestStat_h
7 #include <RooAbsData.h>
9 #include <RooStats/TestStatistic.h>
13 SimplerLikelihoodRatioTestStat(RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt,
const RooArgSet & paramsNull = RooArgSet(),
const RooArgSet & paramsAlt = RooArgSet()) :
24 if (data.numEntries() != 1)
throw std::invalid_argument(
"HybridNew::TestSimpleStatistics: dataset doesn't have exactly 1 entry.");
25 const RooArgSet *
entry = data.get(0);
33 double altNLL =
pdfAlt_->getVal();
35 return -
log(nullNLL/altNLL);
39 return TString::Format(
"-log(%s/%s)",
pdfNull_->GetName(),
pdfAlt_->GetName());
virtual const TString GetVarName() const
std::auto_ptr< RooArgSet > paramsAlt_
std::auto_ptr< RooArgSet > paramsNull_
std::pair< std::string, MonitorElement * > entry
virtual Double_t Evaluate(RooAbsData &data, RooArgSet &nullPOI)
SimplerLikelihoodRatioTestStat(RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt, const RooArgSet ¶msNull=RooArgSet(), const RooArgSet ¶msAlt=RooArgSet())
char data[epos_bytes_allocation]