#include <TrackEfficiency.h>

Inheritance diagram for TrackEfficiency:

Collaboration diagram for TrackEfficiency:

Public Member Functions
	TrackEfficiency (const std::string &name)
virtual void	initialise ()
	standard operation interface
virtual void	execute (const std::vector< TIDA::Track * > &tracks1, const std::vector< TIDA::Track * > &tracks2, TrackAssociator *matcher)
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

Detailed Description

Definition at line 25 of file TrackEfficiency.h.

Constructor & Destructor Documentation

◆ TrackEfficiency()

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

inline

Definition at line 33 of file TrackEfficiency.h.

                                           : 
    TrackAnalysis( name )
  { }  

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.

124{ return m_xBeamReference; }

TrackAnalysis::m_xBeamReference

double m_xBeamReference

beamline positions reference sample

Definition TrackAnalysis.h:153

◆ beamRefy()

double TrackAnalysis::beamRefy ( ) const

inlineinherited

Definition at line 125 of file TrackAnalysis.h.

125{ return m_yBeamReference; }

TrackAnalysis::m_yBeamReference

double m_yBeamReference

Definition TrackAnalysis.h:154

◆ beamRefz()

double TrackAnalysis::beamRefz ( ) const

inlineinherited

Definition at line 126 of file TrackAnalysis.h.

126{ return m_zBeamReference; }

TrackAnalysis::m_zBeamReference

double m_zBeamReference

Definition TrackAnalysis.h:155

◆ beamTestx()

double TrackAnalysis::beamTestx ( ) const

inlineinherited

Definition at line 128 of file TrackAnalysis.h.

128{ return m_xBeamTest; }

TrackAnalysis::m_xBeamTest

double m_xBeamTest

test sample

Definition TrackAnalysis.h:158

◆ beamTesty()

double TrackAnalysis::beamTesty ( ) const

inlineinherited

Definition at line 129 of file TrackAnalysis.h.

129{ return m_yBeamTest; }

TrackAnalysis::m_yBeamTest

double m_yBeamTest

Definition TrackAnalysis.h:159

◆ beamTestz()

double TrackAnalysis::beamTestz ( ) const

inlineinherited

Definition at line 130 of file TrackAnalysis.h.

130{ return m_zBeamTest; }

TrackAnalysis::m_zBeamTest

double m_zBeamTest

Definition TrackAnalysis.h:160

◆ begin()

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

inlineinherited

Definition at line 105 of file TrackAnalysis.h.

105{ return m_histos.begin(); }

◆ end()

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

inlineinherited

Definition at line 106 of file TrackAnalysis.h.

106{ return m_histos.end(); }

◆ event()

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

inlineinherited

Definition at line 135 of file TrackAnalysis.h.

135{ return m_event; }

TrackAnalysis::m_event

TIDA::Event * m_event

Definition TrackAnalysis.h:164

◆ 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.

71 { }

◆ 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.

86 { }

◆ execute() [3/5]

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

virtual

Implements TrackAnalysis.

Definition at line 78 of file TrackEfficiency.cxx.

{ 
  std::cout << "execute() " << name() << std::endl;
 
  //    for ( int i=0 ; i<reftracks.size() ; i++ )  std::cout << i << " " << *reftracks[i]  << " ref"  << std::endl;
  //    for ( int i=0 ; i<testtracks.size() ; i++ ) std::cout << i << " " << *testtracks[i] << " tst"  << std::endl;
  
  // for an efficiency, we want at most 1 matching reco track
  // for each true track
  // for a purity, we want at most 1 true track for each 
  // reconstructed track
  // for a resolution we only want to consider matched pairs
  // of tracks
  
  // this is an efficiency
  for ( unsigned i=reftracks.size() ; i-- ;  ) { 
    
     // get the matched reco track for this true track if
     // there is one - returns NULL if no matching track
    
     //      std::cout << i << " " << *reftracks[i] << std::endl;
    
     const TIDA::Track* matchedreco = matcher->matched(reftracks[i]); 
      
     std::cout << "\t\tSUTT Analysis " << name() << "\t" << i << " " << *reftracks[i] << " -> ";
  
     if ( matchedreco )  {
 
    std::cout << *matchedreco << std::endl;
    
    // convert pt in MeV to GeV
    double pTt = reftracks[i]->pT()/1000; // convert to GeV
    double pTr = matchedreco->pT()/1000;   // convert to GeV
 
#if 0
    effvsPT->Fill(  reftracks[i]->pT() );
    effvseta->Fill( reftracks[i]->eta() );
    effvsphi->Fill( reftracks[i]->phi() );
    effvsd0->Fill( reftracks[i]->a0() );
 
 
    // fill the plots
    iptvipt->Fill( 1/pTt, 1/pTt-1/pTr );
    ptvpt->Fill( reftracks[i]->pT(), pTt-pTr );
 
    phivpt->Fill( reftracks[i]->pT(), reftracks[i]->phi()-matchedreco->phi() );
    etavpt->Fill( reftracks[i]->pT(), reftracks[i]->eta()-matchedreco->eta() );
    z0vpt->Fill( reftracks[i]->pT(),  reftracks[i]->z0()-matchedreco->z0() );
    d0vpt->Fill( reftracks[i]->pT(), reftracks[i]->a0()-matchedreco->a0() );
 
    ptveta->Fill( reftracks[i]->eta(), pTt-pTr); 
    iptveta->Fill( reftracks[i]->eta(), 1/pTt-1/pTr); 
    phiveta->Fill( reftracks[i]->eta(), reftracks[i]->phi()-matchedreco->phi() );
    etaveta->Fill( reftracks[i]->eta(), reftracks[i]->eta()-matchedreco->eta() );
    z0veta->Fill( reftracks[i]->eta(), reftracks[i]->z0()-matchedreco->z0() );
 
    z0vz0->Fill( reftracks[i]->z0(), reftracks[i]->z0()-matchedreco->z0() );
    d0vd0->Fill( reftracks[i]->a0(), reftracks[i]->a0()-matchedreco->a0() );
#endif
 
    std::cout << "SUTT pt res " 
          << "\tpt="    << reftracks[i]->pT() 
          << "\tpTr="   << pTt 
          << "\t1/pTt=" << 1/pTt 
          << "\tDpT="   << pTt-pTr  
          << "\tD1/pT=" << 1/pTt-1/pTr
          << std::endl;
 
     }
     else { 
    std::cout << " NULL" << std::endl;
#if 0
 
        effvsPT->FillTrueOnly(  reftracks[i]->pT() );
    effvseta->FillTrueOnly( reftracks[i]->eta() );
    effvsphi->FillTrueOnly( reftracks[i]->phi() );
    effvsd0->FillTrueOnly( reftracks[i]->a0() );
#endif
 
     }      
  }
 
 
 
}

◆ 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.

91 { }

◆ finalise()

void TrackEfficiency::finalise ( )

virtual

Implements TrackAnalysis.

Definition at line 168 of file TrackEfficiency.cxx.

                               {
  // sort out all the histograms and write them out
  std::cout << "SUTT TrackEfficiency::finalise() " << name() << std::endl;
 
#if 0
  gDirectory->pwd();
  Directory d(name());
  d.push();
  effvsPT->Finalise();
  effvseta->Finalise(Effplot::FitQuad);
  effvsphi->Finalise(Effplot::FitQuad);
  effvsd0->Finalise(Effplot::FitQuad);
  effvsPT->Write();
  effvseta->Write();
  effvsphi->Write();
  effvsd0->Write();
 
  iptvipt->Finalise();
  ptvpt->Finalise();
 
  phivpt->Finalise();
  etavpt->Finalise();
  z0vpt->Finalise();
  d0vpt->Finalise();
  
  ptveta->Finalise();
  iptveta->Finalise();
  phiveta->Finalise();
  etaveta->Finalise();
  z0veta->Finalise();
  
  z0vz0->Finalise();
  d0vd0->Finalise();
 
  // maybe put a Write() call into Finalise() ???
  iptvipt->Write();
  ptvpt->Write();
 
  phivpt->Write();
  etavpt->Write();
  z0vpt->Write();
  d0vpt->Write();
  
  ptveta->Write();
  iptveta->Write();
  phiveta->Write();
  etaveta->Write();
  z0veta->Write();
  
  z0vz0->Write();
  d0vd0->Write();
 
 
  d.pop();
 
#endif
} 

◆ getHistograms()

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

inlineinherited

access the histograms

Definition at line 103 of file TrackAnalysis.h.

103{ return m_histos; }

◆ initialise()

void TrackEfficiency::initialise ( )

virtual

standard operation interface

Implements TrackAnalysis.

Definition at line 20 of file TrackEfficiency.cxx.

                                 {
    // book the histograms
 
#if 0
  double limits[23] = { 0,     500,   1000,  1500,  2000,  2500, 
            3000,  4000,  5000,  6000,  7000,  8000, 
            10000, 12000, 14000, 16000, 20000, 25000, 
            30000, 35000, 40000, 50000, 60000 };
#endif
 
 
#if 0
  double limits[13] = { 0,     1000,  2000,  3000,  5000,  7000,  
            10000, 14000, 20000, 30000, 40000, 60000, 100000 };
 
  double limits2[24] = { 0,     1000,  2000,  3000,  5000,  7000, 8000, 8500, 9000, 9500,
             10000, 10500, 11000, 11500, 12000, 12500, 13000, 13500, 14000, 
             20000, 30000, 40000, 60000, 100000 };
#endif
 
  double limits2[18] = { 0,     1000,  2000,  3000,  5000,  7000,
             8000,  9000, 10000, 11000, 12000, 13000, 
             14000, 20000, 30000, 40000, 60000, 100000 };
 
  std::vector<double> limits;
  for ( int i=0 ; i<18 ; i++ ) limits.push_back(limits2[i]);
 
#if 0
  // efficiencies
  effvsPT  = new Effplot( "efficiency vs pt",  limits );
  effvseta = new Effplot( "efficiency vs eta", 8, -3,  3);
  effvsphi = new Effplot( "efficiency vs phi", 8, -M_PI, M_PI );
  //  effvsd0  = new Effplot( "efficiency vs d0", 20, -100, 100 );
  effvsd0  = new Effplot( "efficiency vs d0", 200, -200, 200 );
 
  // resolutions
  iptvipt  = new Resplot( "resolution ipt vs ipt",  20, -0.2, 0.2,   21, -0.01, 0.01);
 
  ptvpt   = new Resplot( "resolution pt vs pt",   limits,   41, -10,   10   );
  phivpt  = new Resplot( "resolution phi vs pt",  limits,   41, -0.01, 0.01 );
  etavpt  = new Resplot( "resolution eta vs pt",  limits,   41, -0.01,  0.01  );
  z0vpt   = new Resplot( "resolution z0 vs pt",   limits,   41, -1,    1  );
  d0vpt   = new Resplot( "resolution d0 vs pt",   limits,   61, -0.3,  0.3  );
 
  ptveta   = new Resplot( "resolution pt vs eta",   12, -3, 3,   41, -10, 10);
  iptveta  = new Resplot( "resolution ipt vs eta",  12, -3, 3,   41, -0.01, 0.01);
  phiveta  = new Resplot( "resolution phi vs eta",  12, -3, 3,   41, -0.01, 0.01);
  etaveta  = new Resplot( "resolution eta vs eta",  12, -3, 3,   41, -0.01, 0.01);
  z0veta   = new Resplot( "resolution z vs eta",    12, -3, 3,   41, -1, 1);
 
  z0vz0  = new Resplot( "resolution z vs z",    40, -200, 200,  41, -1, 1);
  d0vd0  = new Resplot( "resolution d0 vs d0",  20, -5,  5,     61, -0.3, 0.3); 
#endif
 
} 

◆ name()

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

inlineinherited

return identifier

Definition at line 52 of file TrackAnalysis.h.

52{ return m_name; }

TrackAnalysis::m_name

std::string m_name

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

Definition TrackAnalysis.h:145

◆ roi()

const TIDARoiDescriptor * TrackAnalysis::roi ( ) const

inlineinherited

Definition at line 138 of file TrackAnalysis.h.

138{ return m_roi; }

TrackAnalysis::m_roi

TIDARoiDescriptor * m_roi

Definition TrackAnalysis.h:165

◆ 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.

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

y

#define y

x

#define x

z

#define 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.

110{ 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.

136{ m_event=e; }

AllowedVariables::e

e

Definition AsgElectronSelectorTool.cxx:37

◆ setroi()

void TrackAnalysis::setroi ( TIDARoiDescriptor * r )

inlineinherited

Definition at line 139 of file TrackAnalysis.h.

139{ m_roi=r; }

r

int r

Definition globals.cxx:22

◆ store()

TIDA::FeatureStore & TrackAnalysis::store ( )

inlineinherited

Definition at line 133 of file TrackAnalysis.h.

133{ return m_store; }

TrackAnalysis::m_store

TIDA::FeatureStore m_store

Definition TrackAnalysis.h:162

Member Data Documentation

◆ m_event

TIDA::Event* TrackAnalysis::m_event

protectedinherited

Definition at line 164 of file TrackAnalysis.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_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_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:

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrackEfficiency()

Member Function Documentation

◆ addHistogram()

◆ beamRefx()

◆ beamRefy()

◆ beamRefz()

◆ beamTestx()

◆ beamTesty()

◆ beamTestz()

◆ begin()

◆ end()

◆ event()

◆ execute() [1/5]

◆ execute() [2/5]

◆ execute() [3/5]

◆ execute() [4/5]

◆ execute() [5/5]

◆ execute_vtx()

◆ finalise()

◆ getHistograms()

◆ initialise()

◆ name()

◆ roi()

◆ setBeamRef() [1/2]

◆ setBeamRef() [2/2]

◆ setBeamTest() [1/2]

◆ setBeamTest() [2/2]

◆ setevent()

◆ setroi()

◆ store()

Member Data Documentation

◆ m_event

◆ m_histos

◆ m_name

◆ m_roi

◆ m_store

◆ m_xBeamReference

◆ m_xBeamTest

◆ m_yBeamReference

◆ m_yBeamTest

◆ m_zBeamReference

◆ m_zBeamTest