TrigElectron.cxx File Reference

#include "CaloGeoHelpers/CaloSampling.h"
#include "TrigParticle/TrigElectron.h"
#include "TruthUtils/ParticleConstants.h"

Include dependency graph for TrigElectron.cxx:

Go to the source code of this file.

Functions
bool	operator== (const TrigElectron &te1, const TrigElectron &te2)
	Two TrigElectron are the same if they have the same cluster and the same track (and if valid, i.e.
void	diff (const TrigElectron &te1, const TrigElectron &te2, std::map< std::string, double > &v_diff)
	Function compares two objects and returns "semi verbose" output in the form of map v_diff where varaible names and relative differences between two objects are returned.
MsgStream &	operator<< (MsgStream &msg, const TrigElectron &te)

Function Documentation

diff()

void diff	(	const TrigElectron &	te1,
		const TrigElectron &	te2,
		std::map< std::string, double > &	v_diff )

Function compares two objects and returns "semi verbose" output in the form of map v_diff where varaible names and relative differences between two objects are returned.

Definition at line 328 of file TrigElectron.cxx.

                                                                                                   {
 
  v_diff["RoI"]            = te1.roiId() - te2.roiId();
  v_diff["Charge"]         = std::abs( te1.charge() - te2.charge() );
  v_diff["Eta"]            = fabs( te1.eta() - te2.eta() );
  double d_phi             = fabs( te1.phi() - te2.phi() );
  v_diff["Phi"]            = (d_phi < M_PI ? d_phi : 2*M_PI - d_phi);
  v_diff["Pt"]             = 2 * fabs( (fabs(te1.pt()) - fabs(te2.pt())) / (fabs(te1.pt())+fabs(te2.pt())) );
  v_diff["Zvtx"]           = fabs( te1.Zvtx() - te2.Zvtx() );
  v_diff["trkEtaEtCalo"]   = fabs( te1.trkEtaAtCalo() - te2.trkEtaAtCalo() );
  double d_phi_extrapol    = fabs( te1.trkPhiAtCalo() - te2.trkPhiAtCalo() );
  v_diff["trkPhiEtCalo"]   = (d_phi_extrapol < M_PI ? d_phi_extrapol : 2*M_PI - d_phi_extrapol);
  v_diff["TRThits"]        = abs(te1.nTRTHits() - te2.nTRTHits());
  v_diff["TRThitsHiThres"] = abs(te1.nTRTHiThresholdHits() - te2.nTRTHiThresholdHits());
  //  v_diff["E"]              = fabs(te1.E() - te2.E());
  v_diff["F0"]             = fabs(te1.F0() - te2.F0());
  v_diff["F1"]             = fabs(te1.F1() - te2.F1());
  v_diff["F2"]             = fabs(te1.F2() - te2.F2());
  v_diff["F3"]             = fabs(te1.F3() - te2.F3());
 
  return;
}

operator<<()

MsgStream & operator<<	(	MsgStream &	msg,
		const TrigElectron &	te )

Definition at line 353 of file TrigElectron.cxx.

                                                                 {
  return ( msg << "TrigElectron " << (te.isValid() ? "(valid)" : "(not valid)") 
       << ": RoI = " << te.roiId() << "; Q = " << te.charge() 
       << "; Et/Pt = " << te.EtOverPt() << "; Pt   = " << te.pt()
       << "; TRT hits = " << te.nTRTHits() << "; TRT HiThreshold hits = " << te.nTRTHiThresholdHits() 
       << "; trkPt = " << te.trkPt() << " +- " << te.err_Pt()
       << "; eta = " << te.eta()   << " +- " << te.err_eta()
       << "; phi = " << te.phi()   << " +- " << te.err_phi()
       << "; Zvt = " << te.Zvtx()  << " +- " << te.err_Zvtx()
       << "; DeltaEta1 = " << te.DeltaEta1() << "; DeltaPhi2 = " << te.DeltaPhi2() 
       //      << "; E = "<<te.E()
       <<"; F0 = "<<te.F0()<<"; F1 = "<<te.F1()<<"; F2 = "<<te.F2()<<"; F3 = "<<te.F3() );
}

operator==()

bool operator==	(	const TrigElectron &	te1,
		const TrigElectron &	te2 )

Two TrigElectron are the same if they have the same cluster and the same track (and if valid, i.e.

well built)

Definition at line 317 of file TrigElectron.cxx.

                                                                     {
  return ( te1.isValid() && te2.isValid() 
       && te1.cluster() == te2.cluster() && te1.track() == te2.track());
}