#include <ProfiledLikelihoodRatioTestStatExt.h>

Public Member Functions
virtual Double_t	Evaluate (RooAbsData &data, RooArgSet &nullPOI)
virtual const TString	GetVarName () const
	ProfiledLikelihoodRatioTestStatOpt (const RooArgSet &obs, RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt, const RooArgSet *nuisances, const RooArgSet &paramsNull=RooArgSet(), const RooArgSet &paramsAlt=RooArgSet(), int verbosity=0)
void	setPrintLevel (Int_t level)
Private Member Functions
bool	createNLL (RooAbsPdf &pdf, RooAbsData &data, std::auto_ptr< RooAbsReal > &nll)
double	minNLL (std::auto_ptr< RooAbsReal > &nll)
Private Attributes
std::auto_ptr< RooAbsReal >	nllAlt_
std::auto_ptr< RooAbsReal >	nllNull_
RooArgSet	nuisances_
std::auto_ptr< RooArgSet >	paramsAlt_
std::auto_ptr< RooArgSet >	paramsNull_
RooAbsPdf *	pdfAlt_
RooAbsPdf *	pdfNull_
RooArgSet	snapAlt_
RooArgSet	snapNull_
Int_t	verbosity_

Detailed Description

Definition at line 19 of file ProfiledLikelihoodRatioTestStatExt.h.

Constructor & Destructor Documentation

ProfiledLikelihoodRatioTestStatOpt::ProfiledLikelihoodRatioTestStatOpt	(	const RooArgSet &	obs,
		RooAbsPdf &	pdfNull,
		RooAbsPdf &	pdfAlt,
		const RooArgSet *	nuisances,
		const RooArgSet &	paramsNull = `RooArgSet()`,
		const RooArgSet &	paramsAlt = `RooArgSet()`,
		int	verbosity = `0`
	)

Definition at line 40 of file ProfiledLikelihoodRatioTestStatExt.cc.

References gather_cfg::cout, DBG, nuisances_, pdfAlt_, pdfNull_, snapAlt_, and snapNull_.

: 
    pdfNull_(&pdfNull), pdfAlt_(&pdfAlt),
    paramsNull_(pdfNull_->getVariables()), 
    paramsAlt_(pdfAlt_->getVariables()), 
    verbosity_(verbosity)
{
    snapNull_.addClone(paramsNull);
    snapAlt_.addClone(paramsAlt);
    DBG(DBG_TestStat_params, (std::cout << "Null snapshot" << pdfNull_->GetName())) DBG(DBG_TestStat_params, (snapNull_.Print("V")))
    DBG(DBG_TestStat_params, (std::cout << "Alt  snapshot" << pdfAlt_->GetName()))  DBG(DBG_TestStat_params, (snapAlt_.Print("V")))
    if (nuisances) nuisances_.addClone(*nuisances);
}

Member Function Documentation

bool ProfiledLikelihoodRatioTestStatOpt::createNLL	(	RooAbsPdf &	pdf,
		RooAbsData &	data,
		std::auto_ptr< RooAbsReal > &	nll
	)		`[private]`

Definition at line 85 of file ProfiledLikelihoodRatioTestStatExt.cc.

References nuisances_.

Referenced by Evaluate().

{
    if (typeid(pdf) == typeid(RooSimultaneousOpt)) {
        if (nll_.get() == 0) nll_.reset(pdf.createNLL(data, RooFit::Constrain(nuisances_)));
        else ((cacheutils::CachingSimNLL&)(*nll_)).setData(data);
        return true;
    } else {
        nll_.reset(pdf.createNLL(data, RooFit::Constrain(nuisances_)));
        return false;
    }
}

Double_t ProfiledLikelihoodRatioTestStatOpt::Evaluate	(	RooAbsData &	data,
		RooArgSet &	nullPOI
	)		`[virtual]`

Definition at line 59 of file ProfiledLikelihoodRatioTestStatExt.cc.

References gather_cfg::cout, createNLL(), DBG, DBGV, minNLL(), nllAlt_, nllNull_, nuisances_, paramsAlt_, paramsNull_, pdfAlt_, pdfNull_, snapAlt_, and snapNull_.

{
    *paramsNull_ = nuisances_;
    *paramsNull_ = snapNull_;
    *paramsNull_ = nullPOI;
        
    DBGV(DBG_TestStat_params, (std::cout << "Parameters of null pdf (pre fit)" << pdfNull_->GetName()))
    DBGV(DBG_TestStat_params, (paramsNull_->Print("V")))

    bool canKeepNullNLL = createNLL(*pdfNull_, data, nllNull_);
    double nullNLL = minNLL(nllNull_);
    if (!canKeepNullNLL) nllNull_.reset();
    
    *paramsAlt_ = nuisances_;
    *paramsAlt_ = snapAlt_;
        
    DBGV(DBG_TestStat_params, (std::cout << "Parameters of alt pdf " << pdfAlt_->GetName()))
    DBGV(DBG_TestStat_params, (paramsAlt_->Print("V")))
    bool canKeepAltNLL = createNLL(*pdfAlt_, data, nllAlt_);
    double altNLL = minNLL(nllAlt_);
    if (!canKeepAltNLL) nllAlt_.reset();
    
    DBG(DBG_TestStat_params, (printf("Pln: null = %+8.4f, alt = %+8.4f\n", nullNLL, altNLL)))
    return nullNLL-altNLL;
}

virtual const TString ProfiledLikelihoodRatioTestStatOpt::GetVarName ( ) const [inline, virtual]

Definition at line 27 of file ProfiledLikelihoodRatioTestStatExt.h.

References pdfAlt_, and pdfNull_.

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

double ProfiledLikelihoodRatioTestStatOpt::minNLL ( std::auto_ptr< RooAbsReal > & nll ) [private]

Definition at line 98 of file ProfiledLikelihoodRatioTestStatExt.cc.

References gather_cfg::cout, data, nuisances_, and verbosity_.

Referenced by Evaluate().

{
#if defined(DBG_TestStat_NOFIT) && (DBG_TestStat_NOFIT > 0)
    if (verbosity_ > 0) std::cout << "Profiling likelihood for pdf " << pdf.GetName() << std::endl;
    std::auto_ptr<RooAbsReal> nll_(pdf.createNLL(data, RooFit::Constrain(nuisances_)));
    return nll_->getVal();
#endif
    RooMinimizer minim(*nll_);
    minim.setStrategy(0);
    minim.setPrintLevel(verbosity_-2);
    minim.minimize(ROOT::Math::MinimizerOptions::DefaultMinimizerType().c_str(), ROOT::Math::MinimizerOptions::DefaultMinimizerAlgo().c_str());
    if (verbosity_ > 1) {
        std::auto_ptr<RooFitResult> res(minim.save());
        res->Print("V");
    }
    return nll_->getVal();
}