Logo ROOT   6.10/09
Reference Guide
DebuggingSampler.h
Go to the documentation of this file.
1 // @(#)root/roostats:$Id$
2 // Author: Kyle Cranmer, Lorenzo Moneta, Gregory Schott, Wouter Verkerke
3 /*************************************************************************
4  * Copyright (C) 1995-2008, Rene Brun and Fons Rademakers. *
5  * All rights reserved. *
6  * *
7  * For the licensing terms see $ROOTSYS/LICENSE. *
8  * For the list of contributors see $ROOTSYS/README/CREDITS. *
9  *************************************************************************/
10 
11 #ifndef ROOSTATS_DebuggingSampler
12 #define ROOSTATS_DebuggingSampler
13 
14 /** \class DebuggingSampler
15  \ingroup Roostats
16 
17 DebuggingSampler is a simple implementation of the DistributionCreator interface used for debugging.
18 The sampling distribution is uniformly random between [0,1] and is INDEPENDENT of the data. So it is not useful
19 for true statistical tests, but it is useful for debugging.
20 
21 */
22 
23 #include "Rtypes.h"
24 
25 #include <vector>
26 
29 
30 #include "RooRealVar.h"
31 #include "TRandom.h"
32 
33 namespace RooStats {
34 
36 
37  public:
39  fTestStatistic = new RooRealVar("UniformTestStatistic","UniformTestStatistic",0,0,1);
40  fRand = new TRandom();
41  }
42  virtual ~DebuggingSampler() {
43  delete fRand;
44  delete fTestStatistic;
45  }
46 
47  // Main interface to get a ConfInterval, pure virtual
49  paramsOfInterest = paramsOfInterest; // avoid warning
50  // normally this method would be complex, but here it is simple for debugging
51  std::vector<Double_t> testStatVec;
52  for(Int_t i=0; i<1000; ++i){
53  testStatVec.push_back( fRand->Uniform() );
54  }
55  return new SamplingDistribution("UniformSamplingDist", "for debugging", testStatVec );
56  }
57 
58  // Main interface to evaluate the test statistic on a dataset
59  virtual Double_t EvaluateTestStatistic(RooAbsData& /*data*/, RooArgSet& /*paramsOfInterest*/) {
60  // data = data; // avoid warning
61  // paramsOfInterest = paramsOfInterest; // avoid warning
62  return fRand->Uniform();
63  }
64 
65  // Get the TestStatistic
66  virtual TestStatistic* GetTestStatistic() const {
67  // TODO change to Roo... notifications
68  std::cout << "GetTestStatistic() IS NOT IMPLEMENTED FOR THIS SAMPLER. Returning NULL." << std::endl;
69  return NULL; /*fTestStatistic;*/
70  }
71 
72  // Get the Confidence level for the test
73  virtual Double_t ConfidenceLevel() const {return 1.-fSize;}
74 
75  // Common Initialization
76  virtual void Initialize(RooAbsArg& /* testStatistic */, RooArgSet& /* paramsOfInterest */, RooArgSet& /* nuisanceParameters */ ) {
77  }
78 
79  // Set the Pdf, add to the the workspace if not already there
80  virtual void SetPdf(RooAbsPdf&) {}
81 
82  // specify the parameters of interest in the interval
83  virtual void SetParameters(RooArgSet&) {}
84  // specify the nuisance parameters (eg. the rest of the parameters)
85  virtual void SetNuisanceParameters(const RooArgSet&) {}
86  // specify the values of parameters used when evaluating test statistic
87  virtual void SetParametersForTestStat(const RooArgSet& ) {}
88  // specify the conditional observables
89  virtual void SetGlobalObservables(const RooArgSet& ) {}
90 
91 
92  // set the size of the test (rate of Type I error) ( Eg. 0.05 for a 95% Confidence Interval)
93  virtual void SetTestSize(Double_t size) {fSize = size;}
94  // set the confidence level for the interval (eg. 0.95 for a 95% Confidence Interval)
95  virtual void SetConfidenceLevel(Double_t cl) {fSize = 1.-cl;}
96 
97  // Set the TestStatistic (want the argument to be a function of the data & parameter points
98  virtual void SetTestStatistic(TestStatistic* /*testStatistic*/) {
99  // TODO change to Roo... notifications
100  std::cout << "SetTestStatistic(...) IS NOT IMPLEMENTED FOR THIS SAMPLER" << std::endl;
101  }
102 
103  private:
107 
108  protected:
109  ClassDef(DebuggingSampler,1) // A simple implementation of the DistributionCreator interface
110  };
111 }
112 
113 
114 #endif
virtual void SetNuisanceParameters(const RooArgSet &)
TestStatSampler is an interface class for a tools which produce RooStats SamplingDistributions.
int Int_t
Definition: RtypesCore.h:41
#define NULL
Definition: RtypesCore.h:88
virtual SamplingDistribution * GetSamplingDistribution(RooArgSet &paramsOfInterest)
virtual void SetConfidenceLevel(Double_t cl)
#define ClassDef(name, id)
Definition: Rtypes.h:297
This is the base class for the ROOT Random number generators.
Definition: TRandom.h:27
virtual TestStatistic * GetTestStatistic() const
virtual void SetTestSize(Double_t size)
RooRealVar represents a fundamental (non-derived) real valued object.
Definition: RooRealVar.h:36
virtual void SetGlobalObservables(const RooArgSet &)
virtual Double_t EvaluateTestStatistic(RooAbsData &, RooArgSet &)
virtual void SetPdf(RooAbsPdf &)
RooAbsData is the common abstract base class for binned and unbinned datasets.
Definition: RooAbsData.h:37
This class simply holds a sampling distribution of some test statistic.
Namespace for the RooStats classes.
Definition: Asimov.h:20
double Double_t
Definition: RtypesCore.h:55
virtual Double_t Uniform(Double_t x1=1)
Returns a uniform deviate on the interval (0, x1).
Definition: TRandom.cxx:606
virtual void SetParametersForTestStat(const RooArgSet &)
virtual void SetParameters(RooArgSet &)
RooAbsPdf is the abstract interface for all probability density functions The class provides hybrid a...
Definition: RooAbsPdf.h:41
virtual void SetTestStatistic(TestStatistic *)
virtual Double_t ConfidenceLevel() const
RooAbsArg is the common abstract base class for objects that represent a value (of arbitrary type) an...
Definition: RooAbsArg.h:66
virtual void Initialize(RooAbsArg &, RooArgSet &, RooArgSet &)
TestStatistic is an interface class to provide a facility for construction test statistics distributi...
Definition: TestStatistic.h:31