#include <Efficiency2D.h>

Inheritance diagram for Efficiency2D:

Collaboration diagram for Efficiency2D:

Public Member Functions
	Efficiency2D (TH2F *h, const std::string &n="")

	Efficiency2D (TH2F hnum, TH2F hden, const std::string &n, double scale=100)

	~Efficiency2D ()

void	Fill (double x, double y, double w=1)
	fill methods ... More...

void	FillDenom (double x, double y, float w=1)

TGraphAsymmErrors *	BayesX (int slice, double scale=100)
	evaluate the uncertainties correctly ... More...

TGraphAsymmErrors *	BayesY (int slice, double scale=100)

TH1D *	sliceX (int i)

TH1D *	sliceY (int i)

int	slicesX () const

int	slicesY () const

const std::string &	name () const

TH2F *	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)

std::string	slicename (const std::string &s, int i) const

TGraphAsymmErrors *	BayesInternal (TH1 hn, TH1 hd, double scale=100) const

Protected Attributes
std::string	m_name

TH2F *	m_hnumer

TH2F *	m_hdenom

TH2F *	m_hmissed

TH2F *	m_heff

std::vector< int >	m_ibin

Detailed Description

Definition at line 19 of file Efficiency2D.h.

Constructor & Destructor Documentation

◆ Efficiency2D() [1/2]

Efficiency2D::Efficiency2D	(	TH2F *	h,
		const std::string &	n = `""`
	)

inline

Definition at line 23 of file Efficiency2D.h.

23 :

24 T_Efficiency<TH2F>( h, n ) { }

◆ Efficiency2D() [2/2]

Efficiency2D::Efficiency2D	(	TH2F *	hnum,
		TH2F *	hden,
		const std::string &	n,
		double	scale = `100`
	)

inline

Definition at line 26 of file Efficiency2D.h.

                                                                              :
     T_Efficiency<TH2F>( hnum, hden, n, scale ) { 
     getibinvec(true);
     finalise( scale );
   }  

◆ ~Efficiency2D()

Efficiency2D::~Efficiency2D ( )

inline

Definition at line 32 of file Efficiency2D.h.

32 { }

Member Function Documentation

◆ BayesInternal()

TGraphAsymmErrors* T_Efficiency< TH2F >::BayesInternal	(	TH1 *	hn,
		TH1 *	hd,
		double	scale = `100`
	)		const

inlineprotectedinherited

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;
  
   }

◆ BayesX()

TGraphAsymmErrors* Efficiency2D::BayesX	(	int	slice,
		double	scale = `100`
	)

inline

evaluate the uncertainties correctly ...

Definition at line 49 of file Efficiency2D.h.

                                                          { 
     
     if ( slice<0 || slice>=slicesX() ) return 0; 
     
     TH1D* hn = 0;
     TH1D* hd = 0;
     
     if ( m_hnumer ) { 
       // get some slice from the numerator histo
       hn = m_hnumer->ProjectionY( slicename( name()+"_yslice", slice).c_str(), slice+1, slice+1, "e");
     }
     
     if ( m_hdenom ) { 
       // get some slice from the denominator histo
       hd = m_hdenom->ProjectionY( slicename( name()+"_yslice", slice).c_str(), slice+1, slice+1, "e");
     }
     
     return BayesInternal( hn, hd, scale );
     
   }

◆ BayesY()

TGraphAsymmErrors* Efficiency2D::BayesY	(	int	slice,
		double	scale = `100`
	)

inline

Definition at line 72 of file Efficiency2D.h.

                                                          { 
  
     if ( slice<0 || slice>=slicesY() ) return 0; 
  
     TH1D* hn = 0;
     TH1D* hd = 0;
  
     if ( m_hnumer ) { 
       // get some slice from the numerator histo
       hn = m_hnumer->ProjectionX( slicename( name()+"_yslice", slice).c_str(), slice+1, slice+1, "e");
     }
  
     if ( m_hdenom ) { 
       // get some slice from the denominator histo
       hd = m_hdenom->ProjectionX( slicename( name()+"_yslice", slice).c_str(), slice+1, slice+1, "e");
     }
  
     return BayesInternal( hn, hd, scale );
      
   }

◆ Fill()

void Efficiency2D::Fill	(	double	x,
		double	y,
		double	w = `1`
	)

inline

fill methods ...

Definition at line 36 of file Efficiency2D.h.

                                              {
     m_hnumer->Fill( float(x), float(y), float(w) );
     m_hdenom->Fill( float(x), float(y), float(w) );
   }

◆ FillDenom()

void Efficiency2D::FillDenom	(	double	x,
		double	y,
		float	w = `1`
	)

inline

Definition at line 41 of file Efficiency2D.h.

                                                  { 
     m_hmissed->Fill(float(x), float(y), float(w));
     m_hdenom->Fill(float(x), float(y), float(w));
   }

◆ finalise()

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

inlineinherited

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()

double T_Efficiency< TH2F >::findTotalEfficiency

inlineinherited

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.;
   }

◆ getibinvec()

virtual void Efficiency2D::getibinvec ( bool force = false )

inlineprotectedvirtual

Implements T_Efficiency< TH2F >.

Definition at line 118 of file Efficiency2D.h.

                                             { 
     if ( !force && !m_ibin.empty() ) return;
     for ( int i=1 ; i<=m_hdenom->GetNbinsX() ; i++ ) { 
       for ( int j=1 ; j<=m_hdenom->GetNbinsY() ; j++ ) { 
     m_ibin.push_back( m_hdenom->GetBin(i,j) );
       }
     }
   }

◆ Hist()

TH2F * T_Efficiency< TH2F >::Hist

inlineinherited

Definition at line 66 of file T_Efficiency.h.

66 { return m_heff; }

◆ name()

const std::string& T_Efficiency< TH2F >::name

inlineinherited

Definition at line 64 of file T_Efficiency.h.

64 { return m_name; }

◆ slicename()

std::string Efficiency2D::slicename	(	const std::string &	s,
		int	i
	)		const

inlineprotected

Definition at line 128 of file Efficiency2D.h.

                                                          { 
     if ( i<10 ) return s + "0" + std::to_string(i);
     else        return s + std::to_string(i);
   }

◆ slicesX()

int Efficiency2D::slicesX ( ) const

inline

Definition at line 105 of file Efficiency2D.h.

                       { 
     if ( m_hdenom ) return m_hdenom->GetXaxis()->GetNbins();
     return 0;
   }

◆ slicesY()

int Efficiency2D::slicesY ( ) const

inline

Definition at line 111 of file Efficiency2D.h.

                       { 
     if ( m_hdenom ) return m_hdenom->GetYaxis()->GetNbins();
     return 0;
   }

◆ sliceX()

TH1D* Efficiency2D::sliceX ( int i )

inline

Definition at line 94 of file Efficiency2D.h.

                         { 
     if ( i<0 || i>=slicesX() ) return 0; 
     return m_heff->ProjectionY( slicename( name()+"_xslice", i).c_str(), i+1, i+1, "e");
   }

◆ sliceY()

TH1D* Efficiency2D::sliceY ( int i )

inline

Definition at line 99 of file Efficiency2D.h.

                         { 
     if ( i<0 || i>=slicesY() ) return 0; 
     return m_heff->ProjectionX( slicename( name()+"_yslice", i).c_str(), i+1, i+1, "e");
   }

◆ Write()

void T_Efficiency< TH2F >::Write

inlineinherited

Definition at line 153 of file T_Efficiency.h.

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

Member Data Documentation

◆ m_hdenom

TH2F * T_Efficiency< TH2F >::m_hdenom

protectedinherited

Definition at line 227 of file T_Efficiency.h.

◆ m_heff

TH2F * T_Efficiency< TH2F >::m_heff

protectedinherited

Definition at line 231 of file T_Efficiency.h.

◆ m_hmissed

TH2F * T_Efficiency< TH2F >::m_hmissed

protectedinherited

Definition at line 229 of file T_Efficiency.h.

◆ m_hnumer

TH2F * T_Efficiency< TH2F >::m_hnumer

protectedinherited

Definition at line 226 of file T_Efficiency.h.

◆ m_ibin

std::vector<int> T_Efficiency< TH2F >::m_ibin

protectedinherited

Definition at line 233 of file T_Efficiency.h.

◆ m_name

std::string T_Efficiency< TH2F >::m_name

protectedinherited

Definition at line 224 of file T_Efficiency.h.

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

Efficiency2D.h

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Efficiency2D() [1/2]

◆ Efficiency2D() [2/2]

◆ ~Efficiency2D()

Member Function Documentation

◆ BayesInternal()

◆ BayesX()

◆ BayesY()

◆ Fill()

◆ FillDenom()

◆ finalise()

◆ findTotalEfficiency()

◆ getibinvec()

◆ Hist()

◆ name()

◆ slicename()

◆ slicesX()

◆ slicesY()

◆ sliceX()

◆ sliceY()

◆ Write()

Member Data Documentation

◆ m_hdenom

◆ m_heff

◆ m_hmissed

◆ m_hnumer

◆ m_ibin

◆ m_name