/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/HiggsAnalysis/CombinedLimit/interface/SimplerLikelihoodRatioTestStatExt.h

Go to the documentation of this file.

#ifndef HiggsAnalysis_CombinedLimit_SimplerLikelihoodRatioTestStatExt_h
#define HiggsAnalysis_CombinedLimit_SimplerLikelihoodRatioTestStatExt_h

#include <memory>
#include <stdexcept>
#include <RooAbsPdf.h>
#include <RooAbsData.h>
#include <RooSimultaneous.h>
#include <RooArgSet.h>
#include <RooStats/TestStatistic.h>


class SimplerLikelihoodRatioTestStatOpt : public RooStats::TestStatistic {
    public:
        SimplerLikelihoodRatioTestStatOpt(const RooArgSet &obs, RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt, const RooArgSet & paramsNull = RooArgSet(), const RooArgSet & paramsAlt = RooArgSet(), bool factorize=true) ;

        virtual ~SimplerLikelihoodRatioTestStatOpt() ;

        virtual Double_t Evaluate(RooAbsData& data, RooArgSet& nullPOI) ;

        virtual const TString GetVarName() const {
            return TString::Format("-log(%s/%s)", pdfNull_->GetName(), pdfAlt_->GetName()); 
        }
    private:
        const RooArgSet *obs_; 
        RooAbsPdf *pdfNull_, *pdfAlt_;
        RooArgSet pdfCompNull_, pdfCompAlt_;
        RooArgList pdfDepObs_;
        RooArgSet snapNull_, snapAlt_; 
        std::auto_ptr<RooArgSet> paramsNull_, paramsAlt_;
        std::auto_ptr<RooAbsPdf> pdfNullOwned_, pdfAltOwned_;
        RooSimultaneous *simPdfNull_, *simPdfAlt_;
        std::vector<RooAbsPdf *> simPdfComponentsNull_, simPdfComponentsAlt_;

        double evalSimNLL(RooAbsData &data,  RooSimultaneous *pdf, std::vector<RooAbsPdf *> &components);
        double evalSimpleNLL(RooAbsData &data,  RooAbsPdf *pdf);
        void unrollSimPdf(RooSimultaneous *pdf, std::vector<RooAbsPdf *> &out);

}; // 

// ===== This below is identical to the RooStats::SimpleLikelihoodRatioTestStat also in implementation
//       I've made a copy here just to be able to put some debug hooks inside.
#if 0
class SimplerLikelihoodRatioTestStatExt : public RooStats::TestStatistic {
    public:
        SimplerLikelihoodRatioTestStatExt(const RooArgSet &obs, RooAbsPdf &pdfNull, RooAbsPdf &pdfAlt, const RooArgSet & paramsNull = RooArgSet(), const RooArgSet & paramsAlt = RooArgSet()) ;

        virtual ~SimplerLikelihoodRatioTestStatExt() ;

        virtual Double_t Evaluate(RooAbsData& data, RooArgSet& nullPOI) ;

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

    private:
        RooAbsPdf *pdfNull_, *pdfAlt_;
        RooArgSet snapNull_, snapAlt_; 
        std::auto_ptr<RooArgSet> paramsNull_, paramsAlt_;
        std::auto_ptr<RooAbsPdf> pdfNullOwned_, pdfAltOwned_;
}; // TestSimpleStatistics
#endif

#endif