Analysis_Efficiency Class Reference

#include <Analysis_Efficiency.h>

Inheritance diagram for Analysis_Efficiency:

Collaboration diagram for Analysis_Efficiency:

Public Member Functions
	Analysis_Efficiency (const std::string &name, double pTCut, double etaCut, double d0Cut, double z0Cut)
virtual void	initialise ()
	standard operation interface
virtual void	execute (const std::vector< TIDA::Track * > &referenceTracks, const std::vector< TIDA::Track * > &testTracks, TrackAssociator *associator)
virtual void	finalise ()
virtual void	execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, TrigObjectMatcher *)
virtual void	execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *)
virtual void	execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher, const TIDA::Event *event, double *)
virtual void	execute (const std::vector< TIDA::Track * > &, const std::vector< TIDA::Track * > &, TrackAssociator *, const TIDA::Event *, double *, TIDARoiDescriptor *)
const std::string &	name () const
	return identifier
virtual void	execute_vtx (const std::vector< TIDA::Vertex * > &, const std::vector< TIDA::Vertex * > &, const TIDA::Event *=0)
void	addHistogram (TH1 *h)
const std::map< std::string, TH1 * > &	getHistograms () const
	access the histograms
std::map< std::string, TH1 * >::const_iterator	begin () const
std::map< std::string, TH1 * >::const_iterator	end () const
void	setBeamRef (double x, double y, double z=0)
	set the beamline positions
void	setBeamRef (const std::vector< double > &v)
void	setBeamTest (double x, double y, double z=0)
void	setBeamTest (const std::vector< double > &v)
double	beamRefx () const
double	beamRefy () const
double	beamRefz () const
double	beamTestx () const
double	beamTesty () const
double	beamTestz () const
TIDA::FeatureStore &	store ()
const TIDA::Event *	event () const
void	setevent (TIDA::Event *e)
const TIDARoiDescriptor *	roi () const
void	setroi (TIDARoiDescriptor *r)

Protected Attributes
std::string	m_name
	identifier of the of the analysis - also used for the root directory into which the histograms are put
std::map< std::string, TH1 * >	m_histos
	lookup table for the histograms by name - does this need to be in the base class?
double	m_xBeamReference
	beamline positions reference sample
double	m_yBeamReference
double	m_zBeamReference
double	m_xBeamTest
	test sample
double	m_yBeamTest
double	m_zBeamTest
TIDA::FeatureStore	m_store
TIDA::Event *	m_event
TIDARoiDescriptor *	m_roi

Private Member Functions
double	phi (double p)

Private Attributes
TH1 *	m_h_integralD
TH1 *	m_h_integralN
TH1 *	m_h_integral
TH1 *	m_h_etaD
TH1 *	m_h_phiD
TH1 *	m_h_z0D
TH1 *	m_h_d0D
TH1 *	m_h_pTD
TH1 *	m_h_etaN
TH1 *	m_h_phiN
TH1 *	m_h_z0N
TH1 *	m_h_d0N
TH1 *	m_h_pTN
TH1 *	m_h_eta
TH1 *	m_h_phi
TH1 *	m_h_z0
TH1 *	m_h_d0
TH1 *	m_h_pT
double	m_pTCut
double	m_etaCut
double	m_d0Cut
double	m_z0Cut

Detailed Description

Definition at line 23 of file Analysis_Efficiency.h.

Constructor & Destructor Documentation

Analysis_Efficiency()

Analysis_Efficiency::Analysis_Efficiency	(	const std::string &	name,
		double	pTCut,
		double	etaCut,
		double	d0Cut,
		double	z0Cut )

Definition at line 18 of file Analysis_Efficiency.cxx.

                                                                                                                       :
  TrackAnalysis(name), m_pTCut(pTCut), m_etaCut(etaCut), m_d0Cut(d0Cut), m_z0Cut(z0Cut) {                                   
 
  // Create integral efficiency histograms
  m_h_integralD = new TH1D(std::string(m_name+"-denominator-integral").c_str(), std::string(m_name+" integral efficiency denominator").c_str(), 1, 0.0, 1.0);
  m_h_integralN = new TH1D(std::string(m_name+"-numerator-integral").c_str(),   std::string(m_name+" integral efficiency numerator").c_str(),   1, 0.0, 1.0);
  m_h_integral  = new TH1D(std::string(m_name+"-integral").c_str(),             std::string(m_name+" integral efficiency").c_str(),             1, 0.0, 1.0);
  addHistogram(m_h_integralD);
  addHistogram(m_h_integralN);
  addHistogram(m_h_integral);
    
  // Create differential efficiency histograms
  m_h_etaD = new TH1D(std::string(m_name+"-denominator-eta").c_str(), std::string(m_name+" eta efficiency denominator").c_str(), 100, -2.5, 2.5);
  m_h_phiD = new TH1D(std::string(m_name+"-denominator-phi").c_str(), std::string(m_name+" phi efficiency denominator").c_str(), 100, -3.142, 3.142);
  m_h_z0D  = new TH1D(std::string(m_name+"-denominator-z0").c_str(),  std::string(m_name+" z0 efficiency denominator").c_str(),  100, -200, 200);
  m_h_d0D  = new TH1D(std::string(m_name+"-denominator-d0").c_str(),  std::string(m_name+" d0 efficiency denominator").c_str(),  100, -20, 20);
  m_h_pTD  = new TH1D(std::string(m_name+"-denominator-pT").c_str(),  std::string(m_name+" pT efficiency denominator").c_str(),  100, -10000, 10000);
  addHistogram(m_h_etaD);
  addHistogram(m_h_phiD);
  addHistogram(m_h_z0D);
  addHistogram(m_h_d0D);
  addHistogram(m_h_pTD);
  
  m_h_etaN = new TH1D(std::string(m_name+"-numerator-eta").c_str(), std::string(m_name+" eta efficiency numerator").c_str(), 100, -2.5, 2.5);
  m_h_phiN = new TH1D(std::string(m_name+"-numerator-phi").c_str(), std::string(m_name+" phi efficiency numerator").c_str(), 100, -3.142, 3.142);
  m_h_z0N  = new TH1D(std::string(m_name+"-numerator-z0").c_str(),  std::string(m_name+" z0 efficiency numerator").c_str(),  100, -200, 200);
  m_h_d0N  = new TH1D(std::string(m_name+"-numerator-d0").c_str(),  std::string(m_name+" d0 efficiency numerator").c_str(),  100, -20, 20);
  m_h_pTN  = new TH1D(std::string(m_name+"-numerator-pT").c_str(),  std::string(m_name+" pT efficiency numerator").c_str(),  100, -10000, 10000);
  addHistogram(m_h_etaN);
  addHistogram(m_h_phiN);
  addHistogram(m_h_z0N);
  addHistogram(m_h_d0N);
  addHistogram(m_h_pTN);
  
  m_h_eta = new TH1D(std::string(m_name+"-eta").c_str(), std::string(m_name+" eta efficiency").c_str(), 100, -2.5, 2.5);
  m_h_phi = new TH1D(std::string(m_name+"-phi").c_str(), std::string(m_name+" phi efficiency").c_str(), 100, -3.142, 3.142);
  m_h_z0  = new TH1D(std::string(m_name+"-z0").c_str(),  std::string(m_name+" z0 efficiency").c_str(),  100, -200, 200);
  m_h_d0  = new TH1D(std::string(m_name+"-d0").c_str(),  std::string(m_name+" d0 efficiency").c_str(),  100, -20, 20);
  m_h_pT  = new TH1D(std::string(m_name+"-pT").c_str(),  std::string(m_name+" pT efficiency").c_str(),  100, -10000, 10000);
  addHistogram(m_h_eta);
  addHistogram(m_h_phi);
  addHistogram(m_h_z0);
  addHistogram(m_h_d0);
  addHistogram(m_h_pT);
  
}  

Member Function Documentation

addHistogram()

void TrackAnalysis::addHistogram ( TH1 * h )

inlineinherited

Definition at line 97 of file TrackAnalysis.h.

                              { 
    std::string name = h->GetName();
    m_histos.insert( std::map<std::string, TH1*>::value_type( name, h) );
  }

beamRefx()

double TrackAnalysis::beamRefx ( ) const

inlineinherited

Definition at line 124 of file TrackAnalysis.h.

{ return m_xBeamReference; }

beamRefy()

double TrackAnalysis::beamRefy ( ) const

inlineinherited

Definition at line 125 of file TrackAnalysis.h.

{ return m_yBeamReference; }

beamRefz()

double TrackAnalysis::beamRefz ( ) const

inlineinherited

Definition at line 126 of file TrackAnalysis.h.

{ return m_zBeamReference; }

beamTestx()

double TrackAnalysis::beamTestx ( ) const

inlineinherited

Definition at line 128 of file TrackAnalysis.h.

{ return  m_xBeamTest; }

beamTesty()

double TrackAnalysis::beamTesty ( ) const

inlineinherited

Definition at line 129 of file TrackAnalysis.h.

{ return  m_yBeamTest; }

beamTestz()

double TrackAnalysis::beamTestz ( ) const

inlineinherited

Definition at line 130 of file TrackAnalysis.h.

{ return  m_zBeamTest; }

begin()

std::map< std::string, TH1 * >::const_iterator TrackAnalysis::begin ( ) const

inlineinherited

Definition at line 105 of file TrackAnalysis.h.

{ return m_histos.begin(); }

end()

std::map< std::string, TH1 * >::const_iterator TrackAnalysis::end ( ) const

inlineinherited

Definition at line 106 of file TrackAnalysis.h.

{ return m_histos.end(); }

event()

const TIDA::Event * TrackAnalysis::event ( ) const

inlineinherited

Definition at line 135 of file TrackAnalysis.h.

{ return m_event; }

execute() [1/5]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > & , const std::vector< TIDA::Track * > & , TrackAssociator * , const TIDA::Event *  )

inlinevirtual

Reimplemented from TrackAnalysis.

Definition at line 68 of file TrackAnalysis.h.

                               { }

execute() [2/5]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > & , const std::vector< TIDA::Track * > & , TrackAssociator * , const TIDA::Event * , double * , TIDARoiDescriptor *  )

inlinevirtual

Reimplemented from TrackAnalysis.

Definition at line 81 of file TrackAnalysis.h.

                                 { }

execute() [3/5]

void Analysis_Efficiency::execute ( const std::vector< TIDA::Track * > & referenceTracks, const std::vector< TIDA::Track * > & testTracks, TrackAssociator * associator )

virtual

Implements TrackAnalysis.

Definition at line 73 of file Analysis_Efficiency.cxx.

                                               {
  
  // Loop over reference tracks
  std::vector<TIDA::Track*>::const_iterator reference, referenceEnd=referenceTracks.end();
  for(reference=referenceTracks.begin(); reference!=referenceEnd; ++reference) {
    
    // Get reference parameters
    double referenceEta = (*reference)->eta();
    double referencePhi = phi((*reference)->phi());
    double referenceZ0  = (*reference)->z0();
    double referenceD0  = (*reference)->a0();
    double referencePT  = (*reference)->pT();
    
    // Fill denominator
    bool passesEta = (fabs(referenceEta)<m_etaCut);
    bool passesPT  = (fabs(referencePT)>m_pTCut);
    bool passesD0  = (fabs(referenceD0)<m_d0Cut);
    bool passesZ0  = (fabs(referenceZ0)<m_z0Cut);
    if(passesPT&&passesD0&&passesZ0)            m_h_etaD->Fill(referenceEta);
    if(passesPT&&passesD0&&passesZ0&&passesEta) m_h_phiD->Fill(referencePhi);
    if(passesPT&&passesD0&&passesEta)           m_h_z0D->Fill(referenceZ0);
    if(passesPT&&passesZ0&&passesEta)           m_h_d0D->Fill(referenceD0);
    if(passesD0&&passesZ0&&passesEta)           m_h_pTD->Fill(referencePT);
    if(passesPT&&passesD0&&passesZ0&&passesEta) m_h_integralD->Fill(0.5);
    
    // Find matched tracks
    const TIDA::Track* test = associator->matched(*reference);
 
    // Fill histograms
    if(test ) {
        
      // Fill denominator
      if(passesPT&&passesD0&&passesZ0)            m_h_etaN->Fill(referenceEta);
      if(passesPT&&passesD0&&passesZ0&&passesEta) m_h_phiN->Fill(referencePhi);
      if(passesPT&&passesD0&&passesEta)           m_h_z0N->Fill(referenceZ0);
      if(passesPT&&passesZ0&&passesEta)           m_h_d0N->Fill(referenceD0);
      if(passesD0&&passesZ0&&passesEta)           m_h_pTN->Fill(referencePT);
      if(passesPT&&passesD0&&passesZ0&&passesEta) m_h_integralN->Fill(0.5);
      
    }
  }
}

execute() [4/5]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > & tracks1, const std::vector< TIDA::Track * > & tracks2, TrackAssociator * matcher, const TIDA::Event * event, double *  )

inlinevirtual

Reimplemented from TrackAnalysis.

Definition at line 73 of file TrackAnalysis.h.

                      { 
    execute( tracks1, tracks2, matcher, event );
  }

execute() [5/5]

virtual void TrackAnalysis::execute ( const std::vector< TIDA::Track * > & tracks1, const std::vector< TIDA::Track * > & tracks2, TrackAssociator * matcher, TrigObjectMatcher *  )

inlinevirtual

Reimplemented from TrackAnalysis.

Definition at line 61 of file TrackAnalysis.h.

                                   {
    execute( tracks1, tracks2, matcher );
  }

execute_vtx()

virtual void TrackAnalysis::execute_vtx ( const std::vector< TIDA::Vertex * > & , const std::vector< TIDA::Vertex * > & , const TIDA::Event * = 0 )

inlinevirtualinherited

Reimplemented in Analysis_Tier0, and AnalysisR3_Tier0.

Definition at line 89 of file TrackAnalysis.h.

                                      { }  

finalise()

void Analysis_Efficiency::finalise ( )

virtual

Implements TrackAnalysis.

Definition at line 120 of file Analysis_Efficiency.cxx.

                                   {
 
  m_h_phiD->Sumw2(); m_h_phiN->Sumw2(); m_h_phi->Sumw2(); m_h_phi->Divide(m_h_phiN, m_h_phiD, 1.0, 1.0, "B");
  m_h_etaD->Sumw2(); m_h_etaN->Sumw2(); m_h_eta->Sumw2(); m_h_eta->Divide(m_h_etaN, m_h_etaD, 1.0, 1.0, "B");
  m_h_z0D->Sumw2();  m_h_z0N->Sumw2();  m_h_z0->Sumw2();  m_h_z0->Divide(m_h_z0N, m_h_z0D, 1.0, 1.0, "B");
  m_h_d0D->Sumw2();  m_h_d0N->Sumw2();  m_h_d0->Sumw2();  m_h_d0->Divide(m_h_d0N, m_h_d0D, 1.0, 1.0, "B");
  m_h_pTD->Sumw2();  m_h_pTN->Sumw2();  m_h_pT->Sumw2();  m_h_pT->Divide(m_h_pTN, m_h_pTD, 1.0, 1.0, "B");
  m_h_integralD->Sumw2();  m_h_integralN->Sumw2();  m_h_integral->Sumw2();  m_h_integral->Divide(m_h_integralN, m_h_integralD, 1.0, 1.0, "B");
 
} 

getHistograms()

const std::map< std::string, TH1 * > & TrackAnalysis::getHistograms ( ) const

inlineinherited

access the histograms

Definition at line 103 of file TrackAnalysis.h.

{ return m_histos; }

initialise()

void Analysis_Efficiency::initialise ( )

virtual

standard operation interface

Implements TrackAnalysis.

Definition at line 67 of file Analysis_Efficiency.cxx.

                                     {
  
}

name()

const std::string & TrackAnalysis::name ( ) const

inlineinherited

return identifier

Definition at line 52 of file TrackAnalysis.h.

{ return m_name; }

phi()

double Analysis_Efficiency::phi ( double p )

private

Definition at line 133 of file Analysis_Efficiency.cxx.

                                        {
  if(p < -M_PI) p += 2*M_PI;
  if(p >  M_PI) p -= 2*M_PI;
  return p;
}

roi()

const TIDARoiDescriptor * TrackAnalysis::roi ( ) const

inlineinherited

Definition at line 138 of file TrackAnalysis.h.

{ return m_roi; }

setBeamRef() [1/2]

void TrackAnalysis::setBeamRef ( const std::vector< double > & v )

inlineinherited

Definition at line 112 of file TrackAnalysis.h.

                                                 {  
    if ( v.size()>0 ) m_xBeamReference = v[0]; 
    if ( v.size()>1 ) m_yBeamReference = v[1]; 
    if ( v.size()>2 ) m_zBeamReference = v[2]; 
  }

setBeamRef() [2/2]

void TrackAnalysis::setBeamRef ( double x, double y, double z = 0 )

inlineinherited

set the beamline positions

Definition at line 109 of file TrackAnalysis.h.

{  m_xBeamReference = x;  m_yBeamReference = y; m_zBeamReference = z; } 

setBeamTest() [1/2]

void TrackAnalysis::setBeamTest ( const std::vector< double > & v )

inlineinherited

Definition at line 118 of file TrackAnalysis.h.

                                                  { 
    if ( v.size()>0 ) m_xBeamTest = v[0]; 
    if ( v.size()>1 ) m_yBeamTest = v[1]; 
    if ( v.size()>2 ) m_zBeamTest = v[2]; 
  }

setBeamTest() [2/2]

void TrackAnalysis::setBeamTest ( double x, double y, double z = 0 )

inlineinherited

Definition at line 110 of file TrackAnalysis.h.

{  m_xBeamTest      = x;       m_yBeamTest = y;      m_zBeamTest = z; }

setevent()

void TrackAnalysis::setevent ( TIDA::Event * e )

inlineinherited

Definition at line 136 of file TrackAnalysis.h.

{ m_event=e; }

setroi()

void TrackAnalysis::setroi ( TIDARoiDescriptor * r )

inlineinherited

Definition at line 139 of file TrackAnalysis.h.

{ m_roi=r; }

store()

TIDA::FeatureStore & TrackAnalysis::store ( )

inlineinherited

Definition at line 133 of file TrackAnalysis.h.

{ return m_store; }

Member Data Documentation

m_d0Cut

double Analysis_Efficiency::m_d0Cut

private

Definition at line 72 of file Analysis_Efficiency.h.

m_etaCut

double Analysis_Efficiency::m_etaCut

private

Definition at line 71 of file Analysis_Efficiency.h.

m_event

TIDA::Event* TrackAnalysis::m_event

protectedinherited

Definition at line 164 of file TrackAnalysis.h.

m_h_d0

TH1* Analysis_Efficiency::m_h_d0

private

Definition at line 66 of file Analysis_Efficiency.h.

m_h_d0D

TH1* Analysis_Efficiency::m_h_d0D

private

Definition at line 56 of file Analysis_Efficiency.h.

m_h_d0N

TH1* Analysis_Efficiency::m_h_d0N

private

Definition at line 61 of file Analysis_Efficiency.h.

m_h_eta

TH1* Analysis_Efficiency::m_h_eta

private

Definition at line 63 of file Analysis_Efficiency.h.

m_h_etaD

TH1* Analysis_Efficiency::m_h_etaD

private

Definition at line 53 of file Analysis_Efficiency.h.

m_h_etaN

TH1* Analysis_Efficiency::m_h_etaN

private

Definition at line 58 of file Analysis_Efficiency.h.

m_h_integral

TH1* Analysis_Efficiency::m_h_integral

private

Definition at line 50 of file Analysis_Efficiency.h.

m_h_integralD

TH1* Analysis_Efficiency::m_h_integralD

private

Definition at line 48 of file Analysis_Efficiency.h.

m_h_integralN

TH1* Analysis_Efficiency::m_h_integralN

private

Definition at line 49 of file Analysis_Efficiency.h.

m_h_phi

TH1* Analysis_Efficiency::m_h_phi

private

Definition at line 64 of file Analysis_Efficiency.h.

m_h_phiD

TH1* Analysis_Efficiency::m_h_phiD

private

Definition at line 54 of file Analysis_Efficiency.h.

m_h_phiN

TH1* Analysis_Efficiency::m_h_phiN

private

Definition at line 59 of file Analysis_Efficiency.h.

m_h_pT

TH1* Analysis_Efficiency::m_h_pT

private

Definition at line 67 of file Analysis_Efficiency.h.

m_h_pTD

TH1* Analysis_Efficiency::m_h_pTD

private

Definition at line 57 of file Analysis_Efficiency.h.

m_h_pTN

TH1* Analysis_Efficiency::m_h_pTN

private

Definition at line 62 of file Analysis_Efficiency.h.

m_h_z0

TH1* Analysis_Efficiency::m_h_z0

private

Definition at line 65 of file Analysis_Efficiency.h.

m_h_z0D

TH1* Analysis_Efficiency::m_h_z0D

private

Definition at line 55 of file Analysis_Efficiency.h.

m_h_z0N

TH1* Analysis_Efficiency::m_h_z0N

private

Definition at line 60 of file Analysis_Efficiency.h.

m_histos

std::map<std::string, TH1*> TrackAnalysis::m_histos

protectedinherited

lookup table for the histograms by name - does this need to be in the base class?

Definition at line 149 of file TrackAnalysis.h.

m_name

std::string TrackAnalysis::m_name

protectedinherited

identifier of the of the analysis - also used for the root directory into which the histograms are put

Definition at line 145 of file TrackAnalysis.h.

m_pTCut

double Analysis_Efficiency::m_pTCut

private

Definition at line 70 of file Analysis_Efficiency.h.

m_roi

TIDARoiDescriptor* TrackAnalysis::m_roi

protectedinherited

Definition at line 165 of file TrackAnalysis.h.

m_store

TIDA::FeatureStore TrackAnalysis::m_store

protectedinherited

Definition at line 162 of file TrackAnalysis.h.

m_xBeamReference

double TrackAnalysis::m_xBeamReference

protectedinherited

beamline positions reference sample

Definition at line 153 of file TrackAnalysis.h.

m_xBeamTest

double TrackAnalysis::m_xBeamTest

protectedinherited

test sample

Definition at line 158 of file TrackAnalysis.h.

m_yBeamReference

double TrackAnalysis::m_yBeamReference

protectedinherited

Definition at line 154 of file TrackAnalysis.h.

m_yBeamTest

double TrackAnalysis::m_yBeamTest

protectedinherited

Definition at line 159 of file TrackAnalysis.h.

m_z0Cut

double Analysis_Efficiency::m_z0Cut

private

Definition at line 73 of file Analysis_Efficiency.h.

m_zBeamReference

double TrackAnalysis::m_zBeamReference

protectedinherited

Definition at line 155 of file TrackAnalysis.h.

m_zBeamTest

double TrackAnalysis::m_zBeamTest

protectedinherited

Definition at line 160 of file TrackAnalysis.h.

---

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

• Analysis_Efficiency.h
• Analysis_Efficiency.cxx