13 #ifndef RooStats_FeldmanCousins 17 #ifndef RooStats_RooStatsUtils 21 #ifndef RooStats_PointSetInterval 50 FeldmanCousins::FeldmanCousins() :
55 fAdaptiveSampling(false),
58 fDoProfileConstruction(true),
59 fSaveBeltToFile(false),
62 // default constructor
63 // fWS = new RooWorkspace("FeldmanCousinsWS");
64 // fOwnsWorkspace = true;
81 fAdaptiveSampling(false),
82 fAdditionalNToysFactor(1.),
85 fDoProfileConstruction(true),
86 fSaveBeltToFile(false),
165 ooccoutP(&
fModel,
Generation) <<
"FeldmanCousins: Model has nuisance parameters, will do profile construction" << endl;
170 while ((myarg2 = dynamic_cast<RooRealVar*>(it2.
Next()))) {
191 "profileConstruction",
197 *parameters = *parameterScan->
get(i);
199 profileConstructionPoints->
add(*parameters);
215 while ((myarg = dynamic_cast<RooRealVar*>(it.
Next()))) {
Double_t fAdditionalNToysFactor
virtual RooAbsReal * createNLL(RooAbsData &data, const RooLinkedList &cmdList)
Construct representation of -log(L) of PDFwith given dataset.
TIterator * createIterator(Bool_t dir=kIterForward) const
ModelConfig is a simple class that holds configuration information specifying how a model should be u...
ConfidenceBelt * GetConfidenceBelt()
get confidence belt
const RooArgSet * GetObservables() const
get RooArgSet for observables (return NULL if not existing)
virtual Bool_t add(const RooAbsCollection &col, Bool_t silent=kFALSE)
Add a collection of arguments to this collection by calling add() for each element in the source coll...
virtual const RooArgSet * get() const
virtual void SetTestSize(Double_t size)
set the size of the test (rate of Type I error) ( Eg. 0.05 for a 95% Confidence Interval) ...
RooCmdArg CloneData(Bool_t flag)
Double_t getVal(const RooArgSet *set=0) const
virtual PointSetInterval * GetInterval() const
Main interface to get a ConfInterval (will be a PointSetInterval)
TestStatSampler is an interface class for a tools which produce RooStats SamplingDistributions.
void GuessObsAndNuisance(const RooAbsData &data)
guesses Observables and ParametersOfInterest if not already set
RooFit::MsgLevel globalKillBelow() const
virtual void SetNEventsPerToy(const Int_t nevents)
Bool_t equals(const RooAbsCollection &otherColl) const
Check if this and other collection have identically named contents.
static RooMsgService & instance()
Return reference to singleton instance.
void SetParameterPointsToTest(RooAbsData &pointsToTest)
User-defined set of points to test.
virtual PointSetInterval * GetInterval() const
Main interface to get a ConfInterval (will be a PointSetInterval)
RooDataSet is a container class to hold N-dimensional binned data.
void AdditionalNToysFactor(double fact)
give user ability to ask for more toys
virtual ~FeldmanCousins()
destructor if(fOwnsWorkspace && fWS) delete fWS;
void setBins(Int_t nBins, const char *name=0)
RooRealVar represents a fundamental (non-derived) real valued object.
void CreateTestStatSampler() const
initializes fTestStatSampler data member (mutable)
Bool_t fDoProfileConstruction
ToyMCSampler * fTestStatSampler
void setGlobalKillBelow(RooFit::MsgLevel level)
virtual void SetData(RooAbsData &data)
Set the DataSet.
RooAbsData is the common abstract base class for binned and unbinned datasets.
ProfileLikelihoodTestStat is an implementation of the TestStatistic interface that calculates the pro...
ToyMCSampler is an implementation of the TestStatSampler interface.
RooDataSet is a container class to hold unbinned data.
virtual void add(const RooArgSet &row, Double_t weight=1.0, Double_t weightError=0)
Add a data point, with its coordinates specified in the 'data' argset, to the data set...
The FeldmanCousins class (like the Feldman-Cousins technique) is essentially a specific configuration...
void SetTestStatSampler(TestStatSampler &sampler)
in addition to interface we also need: Set the TestStatSampler (eg.
RooAbsData * fPointsToTest
void UseAdaptiveSampling(bool flag=true)
adaptive sampling algorithm to speed up interval caculation
Namespace for the RooStats classes.
RooAbsPdf * GetPdf() const
get model PDF (return NULL if pdf has not been specified or does not exist)
const RooArgSet * GetParametersOfInterest() const
get RooArgSet containing the parameter of interest (return NULL if not existing)
PointSetInterval is a concrete implementation of the ConfInterval interface.
virtual Int_t numEntries() const
Return the number of bins.
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
virtual void SetObservables(const RooArgSet &o)
void CreateParameterPoints() const
initializes fPointsToTest data member (mutable)
void CreateConfBelt(bool flag=true)
should create confidence belt
NeymanConstruction is a concrete implementation of the NeymanConstruction interface that...
RooAbsPdf is the abstract interface for all probability density functions The class provides hybrid a...
TestStatSampler * GetTestStatSampler() const
Returns instance of TestStatSampler.
const RooArgSet * GetNuisanceParameters() const
get RooArgSet containing the nuisance parameters (return NULL if not existing)
void SetLeftSideTailFraction(Double_t leftSideFraction=0.)
fLeftSideTailFraction*fSize defines lower edge of acceptance region.
virtual RooAbsReal * createProfile(const RooArgSet ¶msOfInterest)
Create a RooProfileLL object that eliminates all nuisance parameters in the present function...
void SaveBeltToFile(bool flag=true)
save teh confidence belt to a file
virtual void SetModel(const ModelConfig &)
set the model
ConfidenceBelt * fConfBelt
virtual Int_t numEntries() const