Efficiency1D Class Reference

#include <Efficiency1D.h>

Inheritance diagram for Efficiency1D:

Collaboration diagram for Efficiency1D:

Public Member Functions
	Efficiency1D (TH1F *h, const std::string &n="")
	Efficiency1D (TH1F *hnum, TH1F *hden, const std::string &n, double scale=100)
	~Efficiency1D ()
virtual void	Fill (double x, double w=1)
virtual void	FillDenom (double x, float w=1)
TGraphAsymmErrors *	Bayes (double scale=100)
	evaluate the uncertainties correctly ... More...
std::string	name () const
TH1F *	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)
TGraphAsymmErrors *	BayesInternal (TH1 *hn, TH1 *hd, double scale=100) const

Protected Attributes
std::string	m_name
TH1F *	m_hnumer
TH1F *	m_hdenom
TH1F *	m_hmissed
TH1F *	m_heff
std::vector< int >	m_ibin

Detailed Description

Definition at line 19 of file Efficiency1D.h.

Constructor & Destructor Documentation

Efficiency1D() [1/2]

Efficiency1D::Efficiency1D ( TH1F *  h, const std::string &  n = ""  )

inline

Definition at line 23 of file Efficiency1D.h.

                                               :
    T_Efficiency<TH1F>( h, n ) { }  

Efficiency1D() [2/2]

Efficiency1D::Efficiency1D ( TH1F *  hnum, TH1F *  hden, const std::string &  n, double  scale = 100  )

inline

Definition at line 26 of file Efficiency1D.h.

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

~Efficiency1D()

Efficiency1D::~Efficiency1D ( )

inline

Definition at line 32 of file Efficiency1D.h.

{ } 

Member Function Documentation

Bayes()

TGraphAsymmErrors* Efficiency1D::Bayes ( double  scale = 100 )

inline

evaluate the uncertainties correctly ...

Definition at line 49 of file Efficiency1D.h.

                                             { 
    return BayesInternal( m_hnumer, m_hdenom, scale );
  }

BayesInternal()

TGraphAsymmErrors* T_Efficiency< TH1F >::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 175 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;
 
  }

Fill()

virtual void Efficiency1D::Fill ( double  x, double  w = 1  )

inlinevirtual

Definition at line 37 of file Efficiency1D.h.

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

FillDenom()

virtual void Efficiency1D::FillDenom ( double  x, float  w = 1  )

inlinevirtual

Definition at line 42 of file Efficiency1D.h.

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

finalise()

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

inlineinherited

actually calculate the efficiencies

Definition at line 80 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< TH1F >::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 140 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 Efficiency1D::getibinvec ( bool  force = false )

inlineprotectedvirtual

Implements T_Efficiency< TH1F >.

Definition at line 55 of file Efficiency1D.h.

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

Hist()

TH1F * T_Efficiency< TH1F >::Hist

inlineinherited

Definition at line 75 of file T_Efficiency.h.

{ return m_heff; }

name()

std::string T_Efficiency< TH1F >::name

inlineinherited

Definition at line 73 of file T_Efficiency.h.

{ return m_name; } 

Write()

void T_Efficiency< TH1F >::Write

inlineinherited

Definition at line 162 of file T_Efficiency.h.

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

Member Data Documentation

m_hdenom

TH1F * T_Efficiency< TH1F >::m_hdenom

protectedinherited

Definition at line 236 of file T_Efficiency.h.

m_heff

TH1F * T_Efficiency< TH1F >::m_heff

protectedinherited

Definition at line 240 of file T_Efficiency.h.

m_hmissed

TH1F * T_Efficiency< TH1F >::m_hmissed

protectedinherited

Definition at line 238 of file T_Efficiency.h.

m_hnumer

TH1F * T_Efficiency< TH1F >::m_hnumer

protectedinherited

Definition at line 235 of file T_Efficiency.h.

m_ibin

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

protectedinherited

Definition at line 242 of file T_Efficiency.h.

m_name

std::string T_Efficiency< TH1F >::m_name

protectedinherited

Definition at line 233 of file T_Efficiency.h.

---

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

• Efficiency1D.h