Logo ROOT   6.14/05
Reference Guide
ROCCalc.h
Go to the documentation of this file.
1 #ifndef ROOT_TMVA_ROCCalc
2 #define ROOT_TMVA_ROCCalc
3 
4 #include "Rtypes.h"
5 
6 #include <vector>
7 #include <sstream>
8 #include <iostream>
9 #include <iomanip>
10 
11 class TList;
12 class TTree;
13 class TString;
14 class TH1;
15 class TH1D;
16 class TH2;
17 class TH2F;
18 class TSpline;
19 
20 namespace TMVA {
21 
22  class MsgLogger;
23 
24 
25  class ROCCalc {
26 
27  public:
28  ROCCalc(TH1* mvaS, TH1* mvaB);
29 
30  ~ROCCalc();
31 
32 
33  TH1D* GetROC();
34  // return the signal eff for a given backgr. efficiency
35  Double_t GetEffSForEffBof(Double_t effBref, Double_t &effSerr);
36  // return the cut value
38  // return the area under the ROC curve
40  // return the statistical significance as function of the mva cut value
41  TH1* GetSignificance( Int_t nStot, Int_t nBtot);
42  TH1* GetPurity(Int_t nStot, Int_t nBtot);
43 
44  void ApplySignalAndBackgroundStyle( TH1* sig, TH1* bkg, TH1* any = 0 );
45 
46  TH1* GetMvaSpdf(){return fmvaSpdf;}
47  TH1* GetMvaBpdf(){return fmvaBpdf;}
48 
49  //false if is found some error in mvaS or mvaB
52 
53  private:
56  Int_t fMaxIter; // maximum number of iterations
57  Double_t fAbsTol; // absolute tolerance deviation
58 
59  Bool_t fStatus; //false if is found some error in mvaS or mvaB
60 
63 
64  TH1* fmvaS, *fmvaB; // the input mva distributions
65  TH1* fmvaSpdf, *fmvaBpdf; // the normalized (and rebinned) input mva distributions
66  Float_t fXmin, fXmax; // min and max of the mva distribution
67  Double_t fNevtS; // number of signal events (used in error calculation)
68  Int_t fCutOrientation; //+1 if larger mva value means more signal like, -1 otherwise
70  TSpline* fSplmvaCumS, *fSplmvaCumB; // spline of cumulated mva distributions
79 
80  Double_t fSignalCut; // MVA cut value for last demanded background rejection or signal efficiency
81 
82  mutable MsgLogger* fLogger; //! message logger
83  MsgLogger& Log() const { return *fLogger; }
84 
85  };
86 }
87 #endif
TSpline * fSplB
Definition: ROCCalc.h:69
TSpline * fSplS
Definition: ROCCalc.h:69
Int_t fnStot
Definition: ROCCalc.h:73
float Float_t
Definition: RtypesCore.h:53
TH1 * fmvaSpdf
Definition: ROCCalc.h:65
Int_t fMaxIter
Definition: ROCCalc.h:56
Base class for spline implementation containing the Draw/Paint methods.
Definition: TSpline.h:20
ROCCalc(TH1 *mvaS, TH1 *mvaB)
Definition: ROCCalc.cxx:48
TH1 * fmvaBpdf
Definition: ROCCalc.h:65
UInt_t fNbins
Definition: ROCCalc.h:61
TH1D * inveffBvsS
Definition: ROCCalc.h:78
Basic string class.
Definition: TString.h:131
TH1 * GetMvaSpdf()
Definition: ROCCalc.h:46
int Int_t
Definition: RtypesCore.h:41
bool Bool_t
Definition: RtypesCore.h:59
TSpline * fSplmvaCumB
Definition: ROCCalc.h:70
TH1 * fmvaS
Definition: ROCCalc.h:64
void ResetStatus()
Definition: ROCCalc.h:51
TSpline * fSplmvaCumS
Definition: ROCCalc.h:70
TH1 * fSignificance
Definition: ROCCalc.h:74
Double_t fAbsTol
Definition: ROCCalc.h:57
TH1D * GetROC()
get the ROC curve
Definition: ROCCalc.cxx:184
TH1 * fmvaB
Definition: ROCCalc.h:64
Double_t GetEffSForEffBof(Double_t effBref, Double_t &effSerr)
get the signal efficiency for a particular background efficiency that will be the value of the effici...
Definition: ROCCalc.cxx:298
Double_t Root(Double_t)
Root finding using Brents algorithm; taken from CERNLIB function RZERO.
Definition: ROCCalc.cxx:354
Bool_t fUseSplines
Definition: ROCCalc.h:62
MsgLogger * fLogger
Definition: ROCCalc.h:82
Double_t GetEffForRoot(Double_t theCut)
returns efficiency as function of cut
Definition: ROCCalc.cxx:329
Float_t fXmin
Definition: ROCCalc.h:66
TH1D * rejBvsS
Definition: ROCCalc.h:77
TH1D * effBvsS
Definition: ROCCalc.h:76
TH1 * GetPurity(Int_t nStot, Int_t nBtot)
Definition: ROCCalc.cxx:435
TH1 * fmvaScumul
Definition: ROCCalc.h:72
A doubly linked list.
Definition: TList.h:44
Double_t GetSignalReferenceCut()
Definition: ROCCalc.h:37
Bool_t GetStatus()
Definition: ROCCalc.h:50
Service class for 2-Dim histogram classes.
Definition: TH2.h:30
Float_t fXmax
Definition: ROCCalc.h:66
2-D histogram with a float per channel (see TH1 documentation)}
Definition: TH2.h:250
Bool_t fStatus
Definition: ROCCalc.h:59
TH1 * fPurity
Definition: ROCCalc.h:75
TH1 * GetSignificance(Int_t nStot, Int_t nBtot)
Definition: ROCCalc.cxx:446
unsigned int UInt_t
Definition: RtypesCore.h:42
MsgLogger & Log() const
message logger
Definition: ROCCalc.h:83
1-D histogram with a double per channel (see TH1 documentation)}
Definition: TH1.h:610
Double_t fSignalCut
Definition: ROCCalc.h:80
TH1 * fmvaBcumul
Definition: ROCCalc.h:72
double Double_t
Definition: RtypesCore.h:55
Int_t fCutOrientation
Definition: ROCCalc.h:68
The TH1 histogram class.
Definition: TH1.h:56
void ApplySignalAndBackgroundStyle(TH1 *sig, TH1 *bkg, TH1 *any=0)
Definition: ROCCalc.cxx:119
ostringstream derivative to redirect and format output
Definition: MsgLogger.h:59
Abstract ClassifierFactory template that handles arbitrary types.
TH1 * GetMvaBpdf()
Definition: ROCCalc.h:47
Double_t fNevtS
Definition: ROCCalc.h:67
Int_t fnBtot
Definition: ROCCalc.h:73
A TTree object has a header with a name and a title.
Definition: TTree.h:70
Double_t GetROCIntegral()
code to compute the area under the ROC ( rej-vs-eff ) curve
Definition: ROCCalc.cxx:274
TSpline * fSpleffBvsS
Definition: ROCCalc.h:71
const Bool_t kTRUE
Definition: RtypesCore.h:87
~ROCCalc()
destructor
Definition: ROCCalc.cxx:166