#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 ()
const std::string &	name () const
T *	Hist ()
void	finalise (double scale=100)
	actually calculate the efficiencies
double	findTotalEfficiency ()
	these 1D and 2D Fill versionms should never be called directly in fact, are they even needed at all ?
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.

                                         :
    m_name(n.empty() ? std::string(h->GetName())+"_eff" : n) {
 
    m_hnumer = (T*)h->Clone( (m_name+"_n").c_str() );
    m_hdenom = (T*)h->Clone( (m_name+"_d").c_str() );
    m_heff   = (T*)h->Clone( m_name.c_str() );
 
    m_hmissed = (T*)h->Clone( (m_name+"_missed").c_str() );
 
    m_hnumer->Reset();
    m_hdenom->Reset();
    m_heff->Reset();
    m_hmissed->Reset();
    
  } 

◆ 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 45 of file T_Efficiency.h.

                                                               :
    m_name(n+"_eff") {
 
    m_hnumer = (T*)hnum->Clone( (m_name+"_n").c_str() );
    m_hdenom = (T*)hden->Clone( (m_name+"_d").c_str() );
    m_heff   = (T*)hnum->Clone( m_name.c_str() );
 
    m_hmissed = (T*)hnum->Clone( (m_name+"_missed").c_str() );
 
    m_heff->Reset();
    m_hmissed->Reset();
  } 

◆ ~T_Efficiency()

template<typename T>

virtual T_Efficiency< T >::~T_Efficiency ( )

inlinevirtual

Definition at line 62 of file T_Efficiency.h.

62{ }

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 166 of file T_Efficiency.h.

                                                                              { 
 
    for ( int i=1 ; i<=hd->GetNbinsX() ; i++ ) {  // cppcheck-suppress [ctunullpointer, nullPointer]; false positive
      double y = hd->GetBinContent(i);
      if ( y==0 ) hd->SetBinContent(i, 1e-20);
    }
 
    TGraphAsymmErrors* tg = new TGraphAsymmErrors( hn, hd, "cl=0.683 b(1,1) mode" );
 
    double* x      = tg->GetX();
    double* y      = tg->GetY();
  
    int n = tg->GetN();
 
    for ( int i=0 ; i<n ; i++ ) { 
    
      y[i]      *= scale;
 
      double yeup   = tg->GetErrorYhigh(i);
      double yedown = tg->GetErrorYlow(i);
      
      yeup   *= scale;
      yedown *= scale;
   
      tg->SetPoint( i, x[i], y[i] ); 
      
      tg->SetPointEYhigh( i, yeup ); 
      tg->SetPointEYlow( i, yedown ); 
 
      tg->SetPointEXhigh( i, 0 ); 
      tg->SetPointEXlow( i, 0 ); 
      
    } 
    
    return tg;
 
  }

◆ finalise()

template<typename T>

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

inline

actually calculate the efficiencies

Definition at line 71 of file T_Efficiency.h.

                                  { 
    getibinvec();
    m_heff->Reset();
    for ( size_t i=0 ; i<m_ibin.size() ; i++ ) { 
      double n = m_hnumer->GetBinContent(m_ibin[i]);
      double d = m_hdenom->GetBinContent(m_ibin[i]);
      
      double e  = 0;
      double ee = 0;
      if ( d!=0 ) { 
    e  = n/d; 
    ee = e*(1-e)/d; 
      } 
      
      // need proper error calculation...
      m_heff->SetBinContent( m_ibin[i], scale*e );
      m_heff->SetBinError( m_ibin[i], scale*std::sqrt(ee) );
 
    }
  }

◆ 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 131 of file T_Efficiency.h.

                               {
 
    getibinvec();
 
    double n_tot = 0;
    double d_tot = 0;
    
    for ( size_t i=0 ; i<m_ibin.size() ; i++ ) { 
      double n = m_hnumer->GetBinContent(m_ibin[i]);
      double d = m_hdenom->GetBinContent(m_ibin[i]);
      
      n_tot += n;
      d_tot += d;
    }
 
    if ( d_tot!=0 ) {
      return n_tot / d_tot;
    }
 
    return 0.;
  }

template<typename T>

void T_Efficiency< T >::Write ( )

inline

Definition at line 153 of file T_Efficiency.h.

               { 
    m_hnumer->Write();
    m_hdenom->Write();
    m_hmissed->Write();
    m_heff->Write();
  }

Definition at line 224 of file T_Efficiency.h.

The documentation for this class was generated from the following file:

T_Efficiency.h

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ T_Efficiency() [1/2]

◆ T_Efficiency() [2/2]

◆ ~T_Efficiency()

Member Function Documentation

◆ BayesInternal()

◆ finalise()

◆ findTotalEfficiency()

◆ getibinvec()

◆ Hist()

◆ name()

◆ Write()

Member Data Documentation

◆ m_hdenom

◆ m_heff

◆ m_hmissed

◆ m_hnumer

◆ m_ibin

◆ m_name