HbbAnalysis/src/EffUtility.cc

#include <cmath>


#include "UserCode/HbbAnalysis/interface/EffUtility.hh"

namespace HbbAnalysis {//namespace

  EffUtility::EffUtility():
    isInitialised_(false),
    nBins_(0),
    xMin_(0),
    xMax_(0),
    nSelected_(0),
    nTotal_(0)
  {}

  void EffUtility::initialise(const unsigned int aNBins, const double aXmin, const double aXmax){
    nBins_ = aNBins;
    xMin_ = aXmin;
    xMax_ = aXmax;
    nSelected_ = new unsigned int[nBins_];
    nTotal_    = new unsigned int[nBins_];
    
    for (unsigned int lBin(0); lBin<nBins_; lBin++){
      nSelected_[lBin] = 0;
      nTotal_[lBin] = 0;
    }

    isInitialised_ = true;

  }

  EffUtility::~EffUtility(){
    if (isInitialised_){
      delete[] nSelected_;
      delete[] nTotal_;
    }
  }

  const unsigned int EffUtility::numberOfBins(){
    return nBins_;
  }


  const double EffUtility::stepSize(){
    return (xMax_-xMin_)/nBins_;
  }
  

  const double EffUtility::xMin(){
    return xMin_;
  }

  const double EffUtility::xMax(){
    return xMax_;
  }

  void EffUtility::incrementPass(unsigned int aBin){
    nSelected_[aBin]++;
  }

  void EffUtility::incrementTotal(unsigned int aBin){
    nTotal_[aBin]++;
  }

  const float EffUtility::getRatio(unsigned int aBin){
    if (nTotal_[aBin]>0) return static_cast<float>(nSelected_[aBin])/nTotal_[aBin];
    return 0;
  }

  const float EffUtility::getRatioError(unsigned int aBin){
    if (nTotal_[aBin]>0) return sqrt(getRatio(aBin)*(1-getRatio(aBin))/getTotal(aBin));
    return 0;
  }

  const unsigned int EffUtility::getPassed(unsigned int aBin){
    return  nSelected_[aBin];
  }

  const unsigned int EffUtility::getTotal(unsigned int aBin){
    return nTotal_[aBin];
  }


}//namespace


Revision:	1.1.1.1 (vendor branch)
Committed:	Wed Apr 15 10:24:16 2009 UTC (16 years ago) by amagnan
Content type:	text/plain
Branch:	v00-00-00, MAIN
CVS Tags:	v01-00-00, beforeMETHacks, v00-05-03, v00-05-02, v00-05-01, v00-05-00, HbbAnaFor35X, v00-04-02, v00-04-01, v00-04-00, v00-03-00, v00-02-01, v00-02-00, v00-01-00, HbbAnaFor22X, HEAD
Changes since 1.1:	+0 -0 lines
Log Message:	Imported sources
#	User	Rev	Content
1	amagnan	1.1	#include <cmath>
2
3
4			#include "UserCode/HbbAnalysis/interface/EffUtility.hh"
5
6			namespace HbbAnalysis {//namespace
7
8			EffUtility::EffUtility():
9			isInitialised_(false),
10			nBins_(0),
11			xMin_(0),
12			xMax_(0),
13			nSelected_(0),
14			nTotal_(0)
15			{}
16
17			void EffUtility::initialise(const unsigned int aNBins, const double aXmin, const double aXmax){
18			nBins_ = aNBins;
19			xMin_ = aXmin;
20			xMax_ = aXmax;
21			nSelected_ = new unsigned int[nBins_];
22			nTotal_ = new unsigned int[nBins_];
23
24			for (unsigned int lBin(0); lBin<nBins_; lBin++){
25			nSelected_[lBin] = 0;
26			nTotal_[lBin] = 0;
27			}
28
29			isInitialised_ = true;
30
31			}
32
33			EffUtility::~EffUtility(){
34			if (isInitialised_){
35			delete[] nSelected_;
36			delete[] nTotal_;
37			}
38			}
39
40			const unsigned int EffUtility::numberOfBins(){
41			return nBins_;
42			}
43
44
45			const double EffUtility::stepSize(){
46			return (xMax_-xMin_)/nBins_;
47			}
48
49
50			const double EffUtility::xMin(){
51			return xMin_;
52			}
53
54			const double EffUtility::xMax(){
55			return xMax_;
56			}
57
58			void EffUtility::incrementPass(unsigned int aBin){
59			nSelected_[aBin]++;
60			}
61
62			void EffUtility::incrementTotal(unsigned int aBin){
63			nTotal_[aBin]++;
64			}
65
66			const float EffUtility::getRatio(unsigned int aBin){
67			if (nTotal_[aBin]>0) return static_cast<float>(nSelected_[aBin])/nTotal_[aBin];
68			return 0;
69			}
70
71			const float EffUtility::getRatioError(unsigned int aBin){
72			if (nTotal_[aBin]>0) return sqrt(getRatio(aBin)*(1-getRatio(aBin))/getTotal(aBin));
73			return 0;
74			}
75
76			const unsigned int EffUtility::getPassed(unsigned int aBin){
77			return nSelected_[aBin];
78			}
79
80			const unsigned int EffUtility::getTotal(unsigned int aBin){
81			return nTotal_[aBin];
82			}
83
84
85			}//namespace
86
87