#include <BestFitSigmaTestStat.h>

Inheritance diagram for BestFitSigmaTestStat:

Public Member Functions
	BestFitSigmaTestStat (const RooArgSet &observables, RooAbsPdf &pdf, const RooArgSet *nuisances, const RooArgSet &params, int verbosity=0)

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

virtual const TString	GetVarName () const

void	setPrintLevel (Int_t level)

Private Member Functions
bool	createNLL (RooAbsPdf &pdf, RooAbsData &data)

double	minNLL (bool constrained, RooRealVar *r=0)

Private Attributes
std::auto_ptr< RooAbsReal >	nll_

RooArgSet	nuisances_

std::auto_ptr< RooArgSet >	params_

RooAbsPdf *	pdf_

RooArgSet	poi_

RooArgSet	snap_

Int_t	verbosity_

Detailed Description

Definition at line 15 of file BestFitSigmaTestStat.h.

Constructor & Destructor Documentation

BestFitSigmaTestStat::BestFitSigmaTestStat	(	const RooArgSet &	observables,
		RooAbsPdf &	pdf,
		const RooArgSet *	nuisances,
		const RooArgSet &	params,
		int	verbosity = `0`
	)

Definition at line 42 of file BestFitSigmaTestStat.cc.

References nuisances_, params_, pdf_, poi_, and snap_.

 :
     pdf_(&pdf),
     verbosity_(verbosity)
 {
     params.snapshot(snap_,false);
     ((RooRealVar*)snap_.find(params.first()->GetName()))->setConstant(false);
     poi_.add(snap_);
     if (nuisances) { nuisances_.add(*nuisances); snap_.addClone(*nuisances, /*silent=*/true); }
     params_.reset(pdf_->getParameters(observables));
 }

Member Function Documentation

bool BestFitSigmaTestStat::createNLL	(	RooAbsPdf &	pdf,
		RooAbsData &	data
	)

private

Definition at line 93 of file BestFitSigmaTestStat.cc.

References nll_, and 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 BestFitSigmaTestStat::Evaluate	(	RooAbsData &	data,
		RooArgSet &	nullPOI
	)

virtual

Definition at line 60 of file BestFitSigmaTestStat.cc.

References gather_cfg::cout, createNLL(), minNLL(), nll_, params_, pdf_, poi_, alignCSCRings::r, and snap_.

 {
     // Take snapshot of initial state, to restore it at the end 
     RooArgSet initialState; params_->snapshot(initialState);
 
    // Initialize parameters
     *params_ = snap_;
 
     // Initialize signal strength
     RooRealVar *rIn = (RooRealVar *) poi_.first();
     RooRealVar *r   = (RooRealVar *) params_->find(rIn->GetName());
     bool canKeepNLL = createNLL(*pdf_, data);
 
     double initialR = rIn->getVal();
 
     // Perform unconstrained minimization
     r->setVal(initialR == 0 ? 0.5 : 0.5*initialR); //best guess
     r->setConstant(false);
 
     std::cout << "Doing a fit for with " << r->GetName() << " in range [ " << r->getMin() << " , " << r->getMax() << "]" << std::endl;
     std::cout << "Starting point is " << r->GetName() << " = " << r->getVal() << std::endl;
     minNLL(/*constrained=*/false, r);
     double bestFitR = r->getVal();
     std::cout << "Fit result was " << r->GetName() << " = " << r->getVal() << std::endl;
 
     //Restore initial state, to avoid issues with ToyMCSampler
     *params_ = initialState;
 
     if (!canKeepNLL) nll_.reset();
 
     return bestFitR;
 }

virtual const TString BestFitSigmaTestStat::GetVarName ( ) const

inlinevirtual

Definition at line 24 of file BestFitSigmaTestStat.h.

24 { return "mu-hat`"; }

double BestFitSigmaTestStat::minNLL	(	bool	constrained,
		RooRealVar *	r = `0`
	)

private

Definition at line 105 of file BestFitSigmaTestStat.cc.

References CascadeMinimizer::Constrained, CascadeMinimizer::minimize(), alignBH_cfg::mode, nll_, CascadeMinimizer::Unconstrained, and verbosity_.

Referenced by Evaluate().

 {
     CascadeMinimizer::Mode mode(constrained ? CascadeMinimizer::Constrained : CascadeMinimizer::Unconstrained);
     CascadeMinimizer minim(*nll_, mode, r);
     minim.minimize(verbosity_-2);
     return nll_->getVal();
 }