#include <FeldmanCousins.h>

Inheritance diagram for FeldmanCousins:

Public Member Functions
virtual void	applyOptions (const boost::program_options::variables_map &vm)

	FeldmanCousins ()

virtual const std::string &	name () const

virtual bool	run (RooWorkspace w, RooStats::ModelConfig mc_s, RooStats::ModelConfig mc_b, RooAbsData &data, double &limit, double &limitErr, const double hint)

Public Member Functions inherited from LimitAlgo
virtual void	applyDefaultOptions ()

	LimitAlgo ()

	LimitAlgo (const char *desc)

const boost::program_options::options_description &	options () const

virtual void	setNToys (const int)

virtual void	setToyNumber (const int)

virtual	~LimitAlgo ()

Static Private Attributes
static float	rAbsAccuracy_ = 0.1

static float	rRelAccuracy_ = 0.02

static float	toysFactor_ = 1

Additional Inherited Members
Protected Attributes inherited from LimitAlgo
boost::program_options::options_description	options_

Detailed Description

Compute limit using FeldmanCousins++

Author: Giovanni Petrucciani (UCSD)

Definition at line 13 of file FeldmanCousins.h.

Constructor & Destructor Documentation

FeldmanCousins::FeldmanCousins ( )

Definition at line 17 of file FeldmanCousins.cc.

References LimitAlgo::options_, rAbsAccuracy_, rRelAccuracy_, and toysFactor_.

                                :
     LimitAlgo("FeldmanCousins specific options") {
     options_.add_options()
         ("rAbsAcc", boost::program_options::value<float>(&rAbsAccuracy_)->default_value(rAbsAccuracy_), "Absolute accuracy on r to reach to terminate the scan")
         ("rRelAcc", boost::program_options::value<float>(&rRelAccuracy_)->default_value(rRelAccuracy_), "Relative accuracy on r to reach to terminate the scan")
         ("toysFactor", boost::program_options::value<float>(&toysFactor_)->default_value(toysFactor_),   "Increase the toys per point by this factor w.r.t. the minimum from adaptive sampling")
     ;
 }

Member Function Documentation

void FeldmanCousins::applyOptions ( const boost::program_options::variables_map & vm )

virtual

Reimplemented from LimitAlgo.

Definition at line 26 of file FeldmanCousins.cc.

27 {

28 }

virtual const std::string& FeldmanCousins::name ( void ) const

inlinevirtual

Implements LimitAlgo.

Definition at line 19 of file FeldmanCousins.h.

                                          {
     static const std::string name("FeldmanCousins");
     return name;
   }

bool FeldmanCousins::run	(	RooWorkspace *	w,
		RooStats::ModelConfig *	mc_s,
		RooStats::ModelConfig *	mc_b,
		RooAbsData &	data,
		double &	limit,
		double &	limitErr,
		const double *	hint
	)

virtual

Implements LimitAlgo.

Definition at line 30 of file FeldmanCousins.cc.

References cl, gather_cfg::cout, benchmark_cfg::fc, newFWLiteAna::found, i, lowerLimit_, max(), min, alignCSCRings::r, rAbsAccuracy_, rRelAccuracy_, and toysFactor_.

                                                                                                                                                                      {
   RooArgSet  poi(*mc_s->GetParametersOfInterest());
   RooRealVar *r = dynamic_cast<RooRealVar *>(poi.first());
 
   if ((hint != 0) && (*hint > r->getMin())) {
     r->setMax(std::min<double>(3*(*hint), r->getMax()));
   }
 
   RooStats::ModelConfig modelConfig(*mc_s);
   modelConfig.SetSnapshot(poi);
 
   RooStats::FeldmanCousins fc(data, modelConfig);
   fc.FluctuateNumDataEntries(mc_s->GetPdf()->canBeExtended());
   fc.UseAdaptiveSampling(true);
   fc.SetConfidenceLevel(cl);
   fc.AdditionalNToysFactor(toysFactor_);
 
   fc.SetNBins(10);
   do { 
       if (verbose > 1) std::cout << "scan in range [" << r->getMin() << ", " << r->getMax() << "]" << std::endl;
       std::auto_ptr<RooStats::PointSetInterval> fcInterval((RooStats::PointSetInterval *)fc.GetInterval());
       if (fcInterval.get() == 0) return false;
       if (verbose > 1) fcInterval->GetParameterPoints()->Print("V"); 
       RooDataHist* parameterScan = (RooDataHist*) fc.GetPointsToScan();
       int found = -1; 
       if (lowerLimit_) {
           for(int i=0; i<parameterScan->numEntries(); ++i){
               const RooArgSet * tmpPoint = parameterScan->get(i);
               if (!fcInterval->IsInInterval(*tmpPoint)) found = i;
               else break;
           }
       } else {
           for(int i=0; i<parameterScan->numEntries(); ++i){
               const RooArgSet * tmpPoint = parameterScan->get(i);
               bool inside = fcInterval->IsInInterval(*tmpPoint);
               if (inside) found = i;
               else if (found != -1) break;
           }
           if (found == -1) {
             if (verbose) std::cout << "Points are either all inside or all outside the bound." << std::endl;
             return false;
           }
       }
       double fcBefore = (found > -1 ? parameterScan->get(found)->getRealValue(r->GetName()) : r->getMin());
       double fcAfter  = (found < parameterScan->numEntries()-1 ? 
               parameterScan->get(found+1)->getRealValue(r->GetName()) : r->getMax());
       limit = 0.5*(fcAfter+fcBefore);
       limitErr = 0.5*(fcAfter-fcBefore);
       if (verbose > 0) std::cout << "  would be " << r->GetName() << " < " << limit << " +/- "<<limitErr << std::endl;
       r->setMin(std::max(r->getMin(), limit-3*limitErr)); 
       r->setMax(std::min(r->getMax(), limit+3*limitErr));
       if (limitErr < 4*std::max<float>(rAbsAccuracy_, rRelAccuracy_ * limit)) { // make last scan more precise
           fc.AdditionalNToysFactor(4*toysFactor_);
       }
   } while (limitErr > std::max<float>(rAbsAccuracy_, rRelAccuracy_ * limit));
 
   if (verbose > -1) {
       std::cout << "\n -- FeldmanCousins++ -- \n";
       std::cout << "Limit: " << r->GetName() << (lowerLimit_ ? "> " : "< ") << limit << " +/- " << limitErr << " @ " << cl * 100 << "% CL" << std::endl;
   }
   return true;
 }