| 1 |
// @(#)root/hist:$Id: TLimit.h 20882 2007-11-19 11:31:26Z rdm $
|
| 2 |
// Author: Christophe.Delaere@cern.ch 21/08/2002
|
| 3 |
|
| 4 |
#ifndef ROOT_TLimit
|
| 5 |
#define ROOT_TLimit
|
| 6 |
|
| 7 |
#ifndef ROOT_TObject
|
| 8 |
#include "TObject.h"
|
| 9 |
#endif
|
| 10 |
|
| 11 |
#include "TVectorDfwd.h"
|
| 12 |
|
| 13 |
class TConfidenceLevel;
|
| 14 |
class TRandom;
|
| 15 |
class TLimitDataSource;
|
| 16 |
class TArrayD;
|
| 17 |
class TOrdCollection;
|
| 18 |
class TH1;
|
| 19 |
|
| 20 |
//____________________________________________________________________
|
| 21 |
//
|
| 22 |
// TLimit
|
| 23 |
//
|
| 24 |
// This class computes 95% Confidence Levels.
|
| 25 |
//
|
| 26 |
// Implemented by C. Delaere from the mclimit code written by Tom Junk.
|
| 27 |
// reference: HEP-EX/9902006
|
| 28 |
// See http://cern.ch/thomasj/searchlimits/ecl.html for more details.
|
| 29 |
//____________________________________________________________________
|
| 30 |
|
| 31 |
class TLimit {
|
| 32 |
protected:
|
| 33 |
static bool Fluctuate(TLimitDataSource * input, TLimitDataSource * output, bool init,TRandom *, bool stat=false);
|
| 34 |
static Double_t LogLikelihood(Double_t s, Double_t b, Double_t b2, Double_t d);
|
| 35 |
|
| 36 |
public:
|
| 37 |
TLimit() {}
|
| 38 |
virtual ~TLimit() {}
|
| 39 |
static TConfidenceLevel *ComputeLimit(TLimitDataSource * data,
|
| 40 |
Int_t nmc =50000,
|
| 41 |
bool stat = false,
|
| 42 |
TRandom * generator = 0);
|
| 43 |
static TConfidenceLevel *ComputeLimit(Double_t s, Double_t b, Int_t d,
|
| 44 |
Int_t nmc =50000,
|
| 45 |
bool stat = false,
|
| 46 |
TRandom * generator = 0);
|
| 47 |
static TConfidenceLevel *ComputeLimit(Double_t s, Double_t b, Int_t d,
|
| 48 |
TVectorD* se, TVectorD* be, TObjArray*,
|
| 49 |
Int_t nmc =50000,
|
| 50 |
bool stat = false,
|
| 51 |
TRandom * generator = 0);
|
| 52 |
static TConfidenceLevel *ComputeLimit(TH1* s, TH1* b, TH1* d,
|
| 53 |
Int_t nmc =50000,
|
| 54 |
bool stat = false,
|
| 55 |
TRandom * generator = 0);
|
| 56 |
static TConfidenceLevel *ComputeLimit(TH1* s, TH1* b, TH1* d,
|
| 57 |
TVectorD* se, TVectorD* be, TObjArray*,
|
| 58 |
Int_t nmc =50000,
|
| 59 |
bool stat = false,
|
| 60 |
TRandom * generator = 0);
|
| 61 |
private:
|
| 62 |
static TArrayD *fgTable; // a log table... just to speed up calculation
|
| 63 |
static TOrdCollection *fgSystNames; // Collection of systematics names
|
| 64 |
ClassDef(TLimit, 2) // Class to compute 95% CL limits
|
| 65 |
};
|
| 66 |
|
| 67 |
#endif
|
| 68 |
|
