#include <ProfiledLikelihoodRatioTestStat.h>

Inheritance diagram for ProfiledLikelihoodRatioTestStat:

Public Member Functions
virtual Double_t	Evaluate (RooAbsData &data, RooArgSet &nullPOI)

virtual const TString	GetVarName () const

	ProfiledLikelihoodRatioTestStat (RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt, const RooArgSet *nuisances, const RooArgSet &paramsNull=RooArgSet(), const RooArgSet &paramsAlt=RooArgSet())

Private Attributes
RooArgSet	nuisances_

std::auto_ptr< RooArgSet >	paramsAlt_

std::auto_ptr< RooArgSet >	paramsNull_

RooAbsPdf *	pdfAlt_

RooAbsPdf *	pdfNull_

RooArgSet	snapAlt_

RooArgSet	snapNull_

Detailed Description

Definition at line 10 of file ProfiledLikelihoodRatioTestStat.h.

Constructor & Destructor Documentation

ProfiledLikelihoodRatioTestStat::ProfiledLikelihoodRatioTestStat	(	RooAbsPdf &	pdfNull,
		RooAbsPdf &	pdfAlt,
		const RooArgSet *	nuisances,
		const RooArgSet &	paramsNull = `RooArgSet()`,
		const RooArgSet &	paramsAlt = `RooArgSet()`
	)

inline

Definition at line 12 of file ProfiledLikelihoodRatioTestStat.h.

References nuisances_, snapAlt_, and snapNull_.

                                                                                                                                                                                                   : 
             pdfNull_(&pdfNull), pdfAlt_(&pdfAlt),
             paramsNull_(pdfNull_->getVariables()), 
             paramsAlt_(pdfAlt_->getVariables()) 
         {
             snapNull_.addClone(paramsNull);
             snapAlt_.addClone(paramsAlt);
             if (nuisances) nuisances_.addClone(*nuisances);
         }

Member Function Documentation

Double_t ProfiledLikelihoodRatioTestStat::Evaluate	(	RooAbsData &	data,
		RooArgSet &	nullPOI
	)

virtual

Definition at line 6 of file ProfiledLikelihoodRatioTestStat.cc.

References create_public_lumi_plots::log, nuisances_, paramsAlt_, paramsNull_, pdfAlt_, pdfNull_, snapAlt_, and snapNull_.

 {
     if (data.numEntries() != 1) throw std::invalid_argument("HybridNew::TestSimpleStatistics: dataset doesn't have exactly 1 entry.");
     CloseCoutSentry(true);
 
     const RooArgSet *entry = data.get(0);
     *paramsNull_ = *entry;
     *paramsNull_ = nuisances_;
     *paramsNull_ = snapNull_;
     *paramsNull_ = nullPOI;
 
     pdfNull_->fitTo(data, RooFit::Constrain(nuisances_), RooFit::Hesse(0), RooFit::PrintLevel(-1), RooFit::PrintEvalErrors(-1));
     double nullNLL = pdfNull_->getVal();
 
     *paramsAlt_ = *entry;
     *paramsAlt_ = nuisances_;
     *paramsAlt_ = snapAlt_;
 
     pdfAlt_->fitTo(data, RooFit::Constrain(nuisances_), RooFit::Hesse(0), RooFit::PrintLevel(-1), RooFit::PrintEvalErrors(-1));
 
     double altNLL = pdfAlt_->getVal();
     return -log(nullNLL/altNLL);
 }

virtual const TString ProfiledLikelihoodRatioTestStat::GetVarName ( ) const

inlinevirtual

Definition at line 24 of file ProfiledLikelihoodRatioTestStat.h.

References pdfAlt_, and pdfNull_.

                                                  {
             return TString::Format("-log(%s/%s)", pdfNull_->GetName(), pdfAlt_->GetName()); 
         }