#include <Analysis_ResolutionCosmics.h>

Inheritance diagram for Analysis_ResolutionCosmics:

Collaboration diagram for Analysis_ResolutionCosmics:

Public Member Functions
	Analysis_ResolutionCosmics (const std::string &name)
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_res_eta
TH1 *	m_h_res_phi
TH1 *	m_h_res_d0
TH1 *	m_h_res_z0
TH1 *	m_h_res_invpT
TH1 *	m_h_pull_eta
TH1 *	m_h_pull_phi
TH1 *	m_h_pull_d0
TH1 *	m_h_pull_z0
TH1 *	m_h_pull_invpT

Detailed Description

Definition at line 23 of file Analysis_ResolutionCosmics.h.

Constructor & Destructor Documentation

◆ Analysis_ResolutionCosmics()

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

Definition at line 18 of file Analysis_ResolutionCosmics.cxx.

                                                                            : TrackAnalysis(name) {
 
  // Create resolution histograms
  m_h_res_eta   = new TH1D(std::string(m_name+"-Res-eta").c_str(), std::string(m_name+" eta resolution").c_str(), 50, -0.05, 0.05);
  m_h_res_phi   = new TH1D(std::string(m_name+"-Res-phi").c_str(), std::string(m_name+" phi resolution").c_str(), 50, -0.05, 0.05);
  m_h_res_z0    = new TH1D(std::string(m_name+"-Res-z0").c_str(),  std::string(m_name+" z0 resolution").c_str(), 50, -3.0, 3.0);
  m_h_res_d0    = new TH1D(std::string(m_name+"-Res-d0").c_str(),  std::string(m_name+" d0 resolution").c_str(), 50, -3.0, 3.0);
  m_h_res_invpT = new TH1D(std::string(m_name+"-Res-pT").c_str(),  std::string(m_name+" inv-pT resolution").c_str(), 100, -0.0004, 0.0004);
  addHistogram(m_h_res_eta);
  addHistogram(m_h_res_phi);
  addHistogram(m_h_res_z0);
  addHistogram(m_h_res_d0);
  addHistogram(m_h_res_invpT);
  
  // Create pull histograms
  m_h_pull_eta   = new TH1D(std::string(m_name+"-Pull-eta").c_str(), std::string(m_name+" eta pull").c_str(), 50, -20.0, 20.0);
  m_h_pull_phi   = new TH1D(std::string(m_name+"-Pull-phi").c_str(), std::string(m_name+" phi pull").c_str(), 50, -20.0, 20.0);
  m_h_pull_z0    = new TH1D(std::string(m_name+"-Pull-z0").c_str(),  std::string(m_name+" z0 pull").c_str(), 50, -20.0, 20.0);
  m_h_pull_d0    = new TH1D(std::string(m_name+"-Pull-d0").c_str(),  std::string(m_name+" d0 pull").c_str(), 50, -20.0, 20.0);
  m_h_pull_invpT = new TH1D(std::string(m_name+"-Pull-pT").c_str(),  std::string(m_name+" inv-pT pull").c_str(), 50, -20.0, 20.0);
  addHistogram(m_h_pull_eta);
  addHistogram(m_h_pull_phi);
  addHistogram(m_h_pull_z0);
  addHistogram(m_h_pull_d0);
  addHistogram(m_h_pull_invpT);
 
}

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 Analysis_ResolutionCosmics::execute	(	const std::vector< TIDA::Track * > &	referenceTracks,
		const std::vector< TIDA::Track * > &	testTracks,
		TrackAssociator *	associator )

virtual

Implements TrackAnalysis.

Definition at line 54 of file Analysis_ResolutionCosmics.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();
    
    // Find matched tracks
    const TIDA::Track* test = associator->matched(*reference);
 
    // Fill histograms
    if(test) {
 
      // Get test parameters
      double testEta = -test->eta();
      double testPhi = phi(test->phi()+M_PI);
      double testZ0  = test->z0();
      double testD0  = -test->a0();
      double testPT  = -test->pT();
    
      // Skip problematic tracks
      if(referencePT==0 || testPT==0) continue;
      
      // Get errors
      double eeta   = sqrt( test->deta()*test->deta() + (*reference)->deta()*(*reference)->deta() );
      double ephi   = sqrt( test->dphi()*test->dphi() + (*reference)->dphi()*(*reference)->dphi() );
      double ez0    = sqrt( test->dz0() *test->dz0()  + (*reference)->dz0() *(*reference)->dz0() );
      double ed0    = sqrt( test->da0() *test->da0()  + (*reference)->da0() *(*reference)->da0() );
      double einvpT = sqrt( test->dpT() *test->dpT()  + (*reference)->dpT() *(*reference)->dpT() );
      
      // Fill resolution plots
      m_h_res_eta->Fill(fabs(referenceEta)-fabs(testEta));
      m_h_res_phi->Fill(phi(referencePhi-testPhi));
      m_h_res_z0->Fill(fabs(referenceZ0)-fabs(testZ0));
      m_h_res_d0->Fill(fabs(referenceD0)-fabs(testD0));
      m_h_res_invpT->Fill(fabs(1.0/referencePT)-fabs(1.0/testPT));
      
      // Fill pull plots
      if(eeta!=0)   m_h_pull_eta->Fill((fabs(referenceEta)-fabs(testEta))/eeta);           else m_h_pull_eta->Fill(-1000.0);
      if(ephi!=0)   m_h_pull_phi->Fill((phi(referencePhi-testPhi))/ephi);                  else m_h_pull_phi->Fill(-1000.0);
      if(ez0!=0)    m_h_pull_z0->Fill((fabs(referenceZ0)-fabs(testZ0))/ez0);               else m_h_pull_z0->Fill(-1000.0);
      if(ed0!=0)    m_h_pull_d0->Fill((fabs(referenceD0)-fabs(testD0))/ed0);               else m_h_pull_d0->Fill(-1000.0);
      if(einvpT!=0) m_h_pull_invpT->Fill((fabs(1.0/referencePT)-fabs(1.0/testPT))/einvpT); else m_h_pull_invpT->Fill(-1000.0);
 
    }
  }
}

◆ 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 Analysis_ResolutionCosmics::finalise ( )

virtual

Implements TrackAnalysis.

Definition at line 112 of file Analysis_ResolutionCosmics.cxx.

                                          {
 
} 

◆ 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 Analysis_ResolutionCosmics::initialise ( )

virtual

standard operation interface

Implements TrackAnalysis.

Definition at line 48 of file Analysis_ResolutionCosmics.cxx.

                                            {
 
}

◆ name()

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

inlineinherited

return identifier

Definition at line 52 of file TrackAnalysis.h.

52{ return m_name; }

◆ phi()

double Analysis_ResolutionCosmics::phi ( double p )

private

Definition at line 118 of file Analysis_ResolutionCosmics.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.

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_h_pull_d0

TH1* Analysis_ResolutionCosmics::m_h_pull_d0

private

Definition at line 57 of file Analysis_ResolutionCosmics.h.

◆ m_h_pull_eta

TH1* Analysis_ResolutionCosmics::m_h_pull_eta

private

Definition at line 55 of file Analysis_ResolutionCosmics.h.

◆ m_h_pull_invpT

TH1* Analysis_ResolutionCosmics::m_h_pull_invpT

private

Definition at line 59 of file Analysis_ResolutionCosmics.h.

◆ m_h_pull_phi

TH1* Analysis_ResolutionCosmics::m_h_pull_phi

private

Definition at line 56 of file Analysis_ResolutionCosmics.h.

◆ m_h_pull_z0

TH1* Analysis_ResolutionCosmics::m_h_pull_z0

private

Definition at line 58 of file Analysis_ResolutionCosmics.h.

◆ m_h_res_d0

TH1* Analysis_ResolutionCosmics::m_h_res_d0

private

Definition at line 50 of file Analysis_ResolutionCosmics.h.

◆ m_h_res_eta

TH1* Analysis_ResolutionCosmics::m_h_res_eta

private

Definition at line 48 of file Analysis_ResolutionCosmics.h.

◆ m_h_res_invpT

TH1* Analysis_ResolutionCosmics::m_h_res_invpT

private

Definition at line 52 of file Analysis_ResolutionCosmics.h.

◆ m_h_res_phi

TH1* Analysis_ResolutionCosmics::m_h_res_phi

private

Definition at line 49 of file Analysis_ResolutionCosmics.h.

◆ m_h_res_z0

TH1* Analysis_ResolutionCosmics::m_h_res_z0

private

Definition at line 51 of file Analysis_ResolutionCosmics.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

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Analysis_ResolutionCosmics()

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

◆ phi()

◆ roi()

◆ setBeamRef() [1/2]

◆ setBeamRef() [2/2]

◆ setBeamTest() [1/2]

◆ setBeamTest() [2/2]

◆ setevent()

◆ setroi()

◆ store()

Member Data Documentation

◆ m_event

◆ m_h_pull_d0

◆ m_h_pull_eta

◆ m_h_pull_invpT

◆ m_h_pull_phi

◆ m_h_pull_z0

◆ m_h_res_d0

◆ m_h_res_eta

◆ m_h_res_invpT

◆ m_h_res_phi

◆ m_h_res_z0

◆ m_histos

◆ m_name

◆ m_roi

◆ m_store

◆ m_xBeamReference

◆ m_xBeamTest

◆ m_yBeamReference

◆ m_yBeamTest

◆ m_zBeamReference

◆ m_zBeamTest