ATLAS Offline Software
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
T_Efficiency< T > Class Template Referenceabstract

#include <T_Efficiency.h>

Collaboration diagram for T_Efficiency< T >:

Public Member Functions

 T_Efficiency (T *h, const std::string &n)
 
 T_Efficiency (T *hnum, T *hden, const std::string &n, double)
 
virtual ~T_Efficiency ()
 
std::string name () const
 
T * Hist ()
 
void finalise (double scale=100)
 actually calculate the efficiencies More...
 
double findTotalEfficiency ()
 these 1D and 2D Fill versionms should never be called directly in fact, are they even needed at all ?
More...
 
void Write ()
 

Protected Member Functions

virtual void getibinvec (bool force=false)=0
 
TGraphAsymmErrors * BayesInternal (TH1 *hn, TH1 *hd, double scale=100) const
 

Protected Attributes

std::string m_name
 
T * m_hnumer
 
T * m_hdenom
 
T * m_hmissed
 
T * m_heff
 
std::vector< int > m_ibin
 

Detailed Description

template<typename T>
class T_Efficiency< T >

Definition at line 24 of file T_Efficiency.h.

Constructor & Destructor Documentation

◆ T_Efficiency() [1/2]

template<typename T >
T_Efficiency< T >::T_Efficiency ( T *  h,
const std::string &  n 
)
inline

Definition at line 28 of file T_Efficiency.h.

28  {
29 
30  std::string effname = n;
31  if ( effname=="" ) {
32  effname = std::string(h->GetName())+"_eff";
33  }
34 
35  m_name = effname;
36 
37  m_hnumer = (T*)h->Clone( (effname+"_n").c_str() );
38  m_hdenom = (T*)h->Clone( (effname+"_d").c_str() );
39  m_heff = (T*)h->Clone( effname.c_str() );
40 
41  m_hmissed = (T*)h->Clone( (effname+"_missed").c_str() );
42 
43  m_hnumer->Reset();
44  m_hdenom->Reset();
45  m_heff->Reset();
46  m_hmissed->Reset();
47 
48  }

◆ T_Efficiency() [2/2]

template<typename T >
T_Efficiency< T >::T_Efficiency ( T *  hnum,
T *  hden,
const std::string &  n,
double   
)
inline

can't call finialise here - need to call it in the derived class this->finalise(scale);

Definition at line 51 of file T_Efficiency.h.

51  {
52 
53  std::string effname = n+"_eff";
54 
55  m_name = effname;
56 
57  m_hnumer = (T*)hnum->Clone( (effname+"_n").c_str() );
58  m_hdenom = (T*)hden->Clone( (effname+"_d").c_str() );
59  m_heff = (T*)hnum->Clone( effname.c_str() );
60 
61  m_hmissed = (T*)hnum->Clone( (effname+"_missed").c_str() );
62 
63  m_heff->Reset();
64  m_hmissed->Reset();
65 
69  }

◆ ~T_Efficiency()

template<typename T >
virtual T_Efficiency< T >::~T_Efficiency ( )
inlinevirtual

Definition at line 71 of file T_Efficiency.h.

71 { }

Member Function Documentation

◆ BayesInternal()

template<typename T >
TGraphAsymmErrors* T_Efficiency< T >::BayesInternal ( TH1 *  hn,
TH1 *  hd,
double  scale = 100 
) const
inlineprotected

stupid root, told to divide, it skips bins where the nukber of entries is 0 (ok) but then complains that "number of points is not the same as the number of bins" now that would be ok, if these were user input values, but is stupid if this is some root policy. : root decides to do X and then prints a warning so instead, set the bin contents, for these bins to something really, really, really tiny ...

root is just such a pain - can't just get/set specific y values I mean - why some functions to get an x or y value and others only get all the x values, but the functions to return the errors only get ONE AT A TIME ? why isn't there a simple ScaleX() function? All this functionality no one cares about, and basic functionality missing

Definition at line 175 of file T_Efficiency.h.

175  {
176 
185 
186  for ( int i=1 ; i<=hd->GetNbinsX() ; i++ ) {
187  double y = hd->GetBinContent(i);
188  if ( y==0 ) hd->SetBinContent(i, 1e-20);
189  }
190 
191  TGraphAsymmErrors* tg = new TGraphAsymmErrors( hn, hd, "cl=0.683 b(1,1) mode" );
192 
193 
194  double* x = tg->GetX();
195  double* y = tg->GetY();
196 
197  int n = tg->GetN();
198 
199  for ( int i=0 ; i<n ; i++ ) {
200 
201  y[i] *= scale;
202 
203  double yeup = tg->GetErrorYhigh(i);
204  double yedown = tg->GetErrorYlow(i);
205 
206  yeup *= scale;
207  yedown *= scale;
208 
215 
216  tg->SetPoint( i, x[i], y[i] );
217 
218  tg->SetPointEYhigh( i, yeup );
219  tg->SetPointEYlow( i, yedown );
220 
221  tg->SetPointEXhigh( i, 0 );
222  tg->SetPointEXlow( i, 0 );
223 
224  }
225 
226  return tg;
227 
228  }

◆ finalise()

template<typename T >
void T_Efficiency< T >::finalise ( double  scale = 100)
inline

actually calculate the efficiencies

Definition at line 80 of file T_Efficiency.h.

80  {
81  getibinvec();
82  m_heff->Reset();
83  for ( size_t i=0 ; i<m_ibin.size() ; i++ ) {
84  double n = m_hnumer->GetBinContent(m_ibin[i]);
85  double d = m_hdenom->GetBinContent(m_ibin[i]);
86 
87  double e = 0;
88  double ee = 0;
89  if ( d!=0 ) {
90  e = n/d;
91  ee = e*(1-e)/d;
92  }
93 
94  // need proper error calculation...
95  m_heff->SetBinContent( m_ibin[i], scale*e );
96  m_heff->SetBinError( m_ibin[i], scale*std::sqrt(ee) );
97 
98  }
99  }

◆ findTotalEfficiency()

template<typename T >
double T_Efficiency< T >::findTotalEfficiency ( )
inline

these 1D and 2D Fill versionms should never be called directly in fact, are they even needed at all ?

calculate the efficiencies ...

Definition at line 140 of file T_Efficiency.h.

140  {
141 
142  getibinvec();
143 
144  double n_tot = 0;
145  double d_tot = 0;
146 
147  for ( size_t i=0 ; i<m_ibin.size() ; i++ ) {
148  double n = m_hnumer->GetBinContent(m_ibin[i]);
149  double d = m_hdenom->GetBinContent(m_ibin[i]);
150 
151  n_tot += n;
152  d_tot += d;
153  }
154 
155  if ( d_tot!=0 ) {
156  return n_tot / d_tot;
157  }
158 
159  return 0.;
160  }

◆ getibinvec()

template<typename T >
virtual void T_Efficiency< T >::getibinvec ( bool  force = false)
protectedpure virtual

Implemented in Efficiency2D, and Efficiency1D.

◆ Hist()

template<typename T >
T* T_Efficiency< T >::Hist ( )
inline

Definition at line 75 of file T_Efficiency.h.

75 { return m_heff; }

◆ name()

template<typename T >
std::string T_Efficiency< T >::name ( ) const
inline

Definition at line 73 of file T_Efficiency.h.

73 { return m_name; }

◆ Write()

template<typename T >
void T_Efficiency< T >::Write ( )
inline

Definition at line 162 of file T_Efficiency.h.

162  {
163  m_hnumer->Write();
164  m_hdenom->Write();
165  m_hmissed->Write();
166  m_heff->Write();
167  }

Member Data Documentation

◆ m_hdenom

template<typename T >
T* T_Efficiency< T >::m_hdenom
protected

Definition at line 236 of file T_Efficiency.h.

◆ m_heff

template<typename T >
T* T_Efficiency< T >::m_heff
protected

Definition at line 240 of file T_Efficiency.h.

◆ m_hmissed

template<typename T >
T* T_Efficiency< T >::m_hmissed
protected

Definition at line 238 of file T_Efficiency.h.

◆ m_hnumer

template<typename T >
T* T_Efficiency< T >::m_hnumer
protected

Definition at line 235 of file T_Efficiency.h.

◆ m_ibin

template<typename T >
std::vector<int> T_Efficiency< T >::m_ibin
protected

Definition at line 242 of file T_Efficiency.h.

◆ m_name

template<typename T >
std::string T_Efficiency< T >::m_name
protected

Definition at line 233 of file T_Efficiency.h.

The documentation for this class was generated from the following file:
T_Efficiency::getibinvec
virtual void getibinvec(bool force=false)=0
hist_file_dump.d
d
Definition: hist_file_dump.py:137
T_Efficiency::m_name
std::string m_name
Definition: T_Efficiency.h:233
yodamerge_tmp.scale
scale
Definition: yodamerge_tmp.py:138
x
#define x
T_Efficiency::m_hmissed
T * m_hmissed
Definition: T_Efficiency.h:238
TH1::SetBinContent
void SetBinContent(int, double)
Definition: rootspy.cxx:301
lumiFormat.i
int i
Definition: lumiFormat.py:92
beamspotman.n
n
Definition: beamspotman.py:731
perfmonmt-refit.n_tot
n_tot
Definition: perfmonmt-refit.py:104
T_Efficiency::m_heff
T * m_heff
Definition: T_Efficiency.h:240
TH1::GetBinContent
double GetBinContent(int) const
Definition: rootspy.cxx:298
DiTauMassTools::MaxHistStrategyV2::e
e
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:26
y
#define y
h
T_Efficiency::m_ibin
std::vector< int > m_ibin
Definition: T_Efficiency.h:242
T_Efficiency::m_hdenom
T * m_hdenom
Definition: T_Efficiency.h:236
TSU::T
unsigned long long T
Definition: L1TopoDataTypes.h:35
T_Efficiency::m_hnumer
T * m_hnumer
Definition: T_Efficiency.h:235
Generated on Tue Jul 2 2024 21:58:59 for ATLAS Offline Software by   doxygen 1.8.18