SimCluster Class Reference

Monte Carlo truth information used for tracking validation. More...

#include <SimCluster.h>

Public Types
typedef int	Charge
	electric charge type More...
typedef std::vector< SimTrack >::const_iterator	g4t_iterator
typedef reco::GenParticleRefVector::iterator	genp_iterator
	reference to reco::GenParticle More...
typedef math::XYZTLorentzVectorD	LorentzVector
	Lorentz vector. More...
typedef math::XYZPointD	Point
	point in the space More...
typedef math::PtEtaPhiMLorentzVector	PolarLorentzVector
	Lorentz vector. More...
typedef math::XYZVectorD	Vector
	point in the space More...

Public Member Functions
void	addG4Track (const SimTrack &t)
void	addGenParticle (const reco::GenParticleRef &ref)
void	addHitEnergy (float energy)
	add rechit energy More...
void	addRecHitAndFraction (uint32_t hit, float fraction)
	add rechit with fraction More...
void	addSimHit (const PCaloHit &hit)
	add simhit's energy to cluster More...
math::XYZVectorF	boostToCM () const
	Vector to boost to the particle centre of mass frame. More...
float	charge () const
	Electric charge. Note this is taken from the first SimTrack only. More...
void	clearHitsAndFractions ()
	clear the hits and fractions list More...
void	clearHitsEnergy ()
	clear the energies list More...
float	energy () const
	Energy. Note this is taken from the first SimTrack only. More...
float	et () const
	Transverse energy. Note this is taken from the first SimTrack only. More...
float	eta () const
	Momentum pseudorapidity. Note this is taken from the simtrack before the calorimeter. More...
EncodedEventId	eventId () const
	Signal source, crossing number. More...
g4t_iterator	g4Track_begin () const
g4t_iterator	g4Track_end () const
const std::vector< SimTrack > &	g4Tracks () const
genp_iterator	genParticle_begin () const
	iterators More...
genp_iterator	genParticle_end () const
const reco::GenParticleRefVector &	genParticles () const
std::vector< std::pair< uint32_t, float > >	hits_and_energies () const
	Returns list of rechit IDs and energies for this SimCluster. More...
std::vector< std::pair< uint32_t, float > >	hits_and_fractions () const
	Returns list of rechit IDs and fractions for this SimCluster. More...
bool	longLived () const
	is long lived? More...
float	mass () const
	Mass. Note this is taken from the first SimTrack only. More...
float	massSqr () const
	Mass squared. Note this is taken from the first SimTrack only. More...
math::XYZVectorF	momentum () const
	spatial momentum vector More...
float	mt () const
	Transverse mass. Note this is taken from the first SimTrack only. More...
float	mtSqr () const
	Transverse mass squared. Note this is taken from the first SimTrack only. More...
int	numberOfRecHits () const
	Gives the total number of SimHits, in the cluster. More...
int	numberOfSimHits () const
	Gives the total number of SimHits, in the cluster. More...
float	p () const
	Magnitude of momentum vector. Note this is taken from the first SimTrack only. More...
const math::XYZTLorentzVectorF &	p4 () const
	Four-momentum Lorentz vector. Note this is taken from the first SimTrack only. More...
uint64_t	particleId () const
int	pdgId () const
	PDG ID. More...
float	phi () const
	Momentum azimuthal angle. Note this is taken from the first SimTrack only. More...
float	pt () const
	Transverse momentum. Note this is taken from the first SimTrack only. More...
float	px () const
	x coordinate of momentum vector. Note this is taken from the first SimTrack only. More...
float	py () const
	y coordinate of momentum vector. Note this is taken from the first SimTrack only. More...
float	pz () const
	z coordinate of momentum vector. Note this is taken from the first SimTrack only. More...
float	rapidity () const
	Rapidity. Note this is taken from the simtrack before the calorimeter. More...
	SimCluster ()
	SimCluster (const SimTrack &simtrk)
	SimCluster (EncodedEventId eventID, uint32_t particleID)
float	simEnergy () const
	returns the accumulated sim energy in the cluster More...
int	status () const
	Status word. More...
float	theta () const
	Momentum polar angle. Note this is taken from the first SimTrack only. More...
int	threeCharge () const
	Gives charge in unit of quark charge (should be 3 times "charge()") More...
float	y () const
	Same as rapidity(). More...
	~SimCluster ()

Static Public Attributes
static const unsigned int	longLivedTag = 65536
	long lived flag More...

Private Attributes
std::vector< float >	energies_
EncodedEventId	event_
std::vector< float >	fractions_
std::vector< SimTrack >	g4Tracks_
	references to G4 and reco::GenParticle tracks More...
reco::GenParticleRefVector	genParticles_
std::vector< uint32_t >	hits_
uint64_t	nsimhits_
uint32_t	particleId_
float	simhit_energy_
math::XYZTLorentzVectorF	theMomentum_

Friends
std::ostream &	operator<< (std::ostream &s, SimCluster const &tp)

Detailed Description

Monte Carlo truth information used for tracking validation.

Object with references to the original SimTrack and parent and daughter TrackingVertices. Simulation with high (~100) pileup was taking too much memory so the class was slimmed down and copies of the SimHits were removed.

Author

original author unknown, re-engineering and slimming by Subir Sarkar (subir.nosp@m..sar.nosp@m.kar@c.nosp@m.ern..nosp@m.ch), some tweaking and documentation by Mark Grimes (mark..nosp@m.grim.nosp@m.es@br.nosp@m.isto.nosp@m.l.ac..nosp@m.uk).

Date

original date unknown, re-engineering Jan-May 2013

Definition at line 29 of file SimCluster.h.

Member Typedef Documentation

Charge

typedef int SimCluster::Charge

electric charge type

Definition at line 33 of file SimCluster.h.

g4t_iterator

typedef std::vector<SimTrack>::const_iterator SimCluster::g4t_iterator

Definition at line 41 of file SimCluster.h.

genp_iterator

typedef reco::GenParticleRefVector::iterator SimCluster::genp_iterator

reference to reco::GenParticle

Definition at line 40 of file SimCluster.h.

LorentzVector

typedef math::XYZTLorentzVectorD SimCluster::LorentzVector

Lorentz vector.

Definition at line 34 of file SimCluster.h.

Point

typedef math::XYZPointD SimCluster::Point

point in the space

Definition at line 36 of file SimCluster.h.

PolarLorentzVector

typedef math::PtEtaPhiMLorentzVector SimCluster::PolarLorentzVector

Lorentz vector.

Definition at line 35 of file SimCluster.h.

Vector

typedef math::XYZVectorD SimCluster::Vector

point in the space

Definition at line 37 of file SimCluster.h.

Constructor & Destructor Documentation

SimCluster() [1/3]

SimCluster::SimCluster

(

)

Definition at line 11 of file SimCluster.cc.

                       {
  // No operation
}

SimCluster() [2/3]

SimCluster::SimCluster

(

const SimTrack &

simtrk

)

Definition at line 15 of file SimCluster.cc.

                                             {
  addG4Track(simtrk);
  event_ = simtrk.eventId();
  particleId_ = simtrk.trackId();
 
  theMomentum_.SetPxPyPzE(
      simtrk.momentum().px(), simtrk.momentum().py(), simtrk.momentum().pz(), simtrk.momentum().E());
}

References addG4Track(), event_, CoreSimTrack::eventId(), CoreSimTrack::momentum(), particleId_, theMomentum_, and CoreSimTrack::trackId().

SimCluster() [3/3]

SimCluster::SimCluster	(	EncodedEventId	eventID,
		uint32_t	particleID
	)

Definition at line 24 of file SimCluster.cc.

                                                                  {
  event_ = eventID;
  particleId_ = particleID;
}

References event_, EgammaObjectsElectrons_cfi::particleID, and particleId_.

~SimCluster()

SimCluster::~SimCluster

(

)

Definition at line 29 of file SimCluster.cc.

{}

Member Function Documentation

addG4Track()

void SimCluster::addG4Track ( const SimTrack &  t )

inline

Definition at line 72 of file SimCluster.h.

{ g4Tracks_.push_back(t); }

References g4Tracks_, and submitPVValidationJobs::t.

Referenced by SimCluster().

addGenParticle()

void SimCluster::addGenParticle ( const reco::GenParticleRef &  ref )

inline

Definition at line 71 of file SimCluster.h.

{ genParticles_.push_back(ref); }

References genParticles_, and edm::RefVector< C, T, F >::push_back().

addHitEnergy()

void SimCluster::addHitEnergy ( float  energy )

inline

add rechit energy

Definition at line 181 of file SimCluster.h.

{ energies_.emplace_back(energy); }

References energies_, and energy().

addRecHitAndFraction()

void SimCluster::addRecHitAndFraction ( uint32_t  hit, float  fraction  )

inline

add rechit with fraction

Definition at line 175 of file SimCluster.h.

                                                          {
    hits_.emplace_back(hit);
    fractions_.emplace_back(fraction);
  }

References HLT_FULL_cff::fraction, fractions_, and hits_.

addSimHit()

void SimCluster::addSimHit ( const PCaloHit &  hit )

inline

add simhit's energy to cluster

Definition at line 216 of file SimCluster.h.

{ simhit_energy_ += hit.energy(); }

References simhit_energy_.

boostToCM()

math::XYZVectorF SimCluster::boostToCM ( ) const

inline

Vector to boost to the particle centre of mass frame.

Definition at line 97 of file SimCluster.h.

{ return p4().BoostToCM(); }

References p4().

charge()

float SimCluster::charge ( void  ) const

inline

Electric charge. Note this is taken from the first SimTrack only.

Definition at line 85 of file SimCluster.h.

{ return g4Tracks_[0].charge(); }

References g4Tracks_.

Referenced by threeCharge().

clearHitsAndFractions()

void SimCluster::clearHitsAndFractions ( )

inline

clear the hits and fractions list

Definition at line 204 of file SimCluster.h.

                               {
    std::vector<uint32_t>().swap(hits_);
    std::vector<float>().swap(fractions_);
  }

References fractions_, and hits_.

clearHitsEnergy()

void SimCluster::clearHitsEnergy ( )

inline

clear the energies list

Definition at line 210 of file SimCluster.h.

{ std::vector<float>().swap(energies_); }

References energies_.

energy()

float SimCluster::energy ( void  ) const

inline

Energy. Note this is taken from the first SimTrack only.

Definition at line 104 of file SimCluster.h.

{ return p4().E(); }

References p4().

Referenced by addHitEnergy(), HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), and Jet.Jet::rawEnergy().

et()

float SimCluster::et ( ) const

inline

Transverse energy. Note this is taken from the first SimTrack only.

Definition at line 107 of file SimCluster.h.

{ return p4().Et(); }

References p4().

eta()

float SimCluster::eta ( void  ) const

inline

Momentum pseudorapidity. Note this is taken from the simtrack before the calorimeter.

Definition at line 148 of file SimCluster.h.

{ return p4().Eta(); }

References p4().

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

eventId()

EncodedEventId SimCluster::eventId ( ) const

inline

Signal source, crossing number.

Note this is taken from the first SimTrack only, but there shouldn't be any SimTracks from different crossings in the SimCluster.

Definition at line 66 of file SimCluster.h.

{ return event_; }

References event_.

Referenced by ntupleDataFormat.Event::eventIdStr().

g4Track_begin()

g4t_iterator SimCluster::g4Track_begin ( ) const

inline

Definition at line 76 of file SimCluster.h.

{ return g4Tracks_.begin(); }

References g4Tracks_.

g4Track_end()

g4t_iterator SimCluster::g4Track_end ( ) const

inline

Definition at line 77 of file SimCluster.h.

{ return g4Tracks_.end(); }

References g4Tracks_.

g4Tracks()

const std::vector<SimTrack>& SimCluster::g4Tracks ( ) const

inline

Definition at line 82 of file SimCluster.h.

{ return g4Tracks_; }

References g4Tracks_.

genParticle_begin()

genp_iterator SimCluster::genParticle_begin ( ) const

inline

iterators

Definition at line 74 of file SimCluster.h.

{ return genParticles_.begin(); }

References edm::RefVector< C, T, F >::begin(), and genParticles_.

genParticle_end()

genp_iterator SimCluster::genParticle_end ( ) const

inline

Definition at line 75 of file SimCluster.h.

{ return genParticles_.end(); }

References edm::RefVector< C, T, F >::end(), and genParticles_.

genParticles()

const reco::GenParticleRefVector& SimCluster::genParticles ( ) const

inline

Definition at line 80 of file SimCluster.h.

{ return genParticles_; }

References genParticles_.

hits_and_energies()

std::vector<std::pair<uint32_t, float> > SimCluster::hits_and_energies ( ) const

inline

Returns list of rechit IDs and energies for this SimCluster.

Definition at line 193 of file SimCluster.h.

                                                                {
    assert(hits_.size() == energies_.size());
    std::vector<std::pair<uint32_t, float>> result;
    result.reserve(hits_.size());
    for (size_t i = 0; i < hits_.size(); ++i) {
      result.emplace_back(hits_[i], energies_[i]);
    }
    return result;
  }

References cms::cuda::assert(), energies_, hits_, mps_fire::i, and mps_fire::result.

hits_and_fractions()

std::vector<std::pair<uint32_t, float> > SimCluster::hits_and_fractions ( ) const

inline

Returns list of rechit IDs and fractions for this SimCluster.

Definition at line 184 of file SimCluster.h.

                                                                 {
    std::vector<std::pair<uint32_t, float>> result;
    for (size_t i = 0; i < hits_.size(); ++i) {
      result.emplace_back(hits_[i], fractions_[i]);
    }
    return result;
  }

References fractions_, hits_, mps_fire::i, and mps_fire::result.

Referenced by HGCalHitCalibration::analyze(), HGCalShowerSeparation::analyze(), HGVHistoProducerAlgo::layerClusters_to_CaloParticles(), LayerClusterAssociatorByEnergyScoreImpl::makeConnections(), and HGVHistoProducerAlgo::multiClusters_to_CaloParticles().

longLived()

bool SimCluster::longLived ( ) const

inline

is long lived?

Definition at line 166 of file SimCluster.h.

{ return status() & longLivedTag; }

References longLivedTag, and status().

mass()

float SimCluster::mass ( ) const

inline

Mass. Note this is taken from the first SimTrack only.

Definition at line 110 of file SimCluster.h.

{ return p4().M(); }

References p4().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), and massSqr().

massSqr()

float SimCluster::massSqr ( ) const

inline

Mass squared. Note this is taken from the first SimTrack only.

Definition at line 113 of file SimCluster.h.

{ return pow(mass(), 2); }

References mass(), and funct::pow().

momentum()

math::XYZVectorF SimCluster::momentum ( ) const

inline

spatial momentum vector

Definition at line 94 of file SimCluster.h.

{ return p4().Vect(); }

References p4().

mt()

float SimCluster::mt ( ) const

inline

Transverse mass. Note this is taken from the first SimTrack only.

Definition at line 116 of file SimCluster.h.

{ return p4().Mt(); }

References p4().

mtSqr()

float SimCluster::mtSqr ( ) const

inline

Transverse mass squared. Note this is taken from the first SimTrack only.

Definition at line 120 of file SimCluster.h.

{ return p4().Mt2(); }

References p4().

numberOfRecHits()

int SimCluster::numberOfRecHits ( void  ) const

inline

Gives the total number of SimHits, in the cluster.

Definition at line 172 of file SimCluster.h.

{ return hits_.size(); }

References hits_.

numberOfSimHits()

int SimCluster::numberOfSimHits ( ) const

inline

Gives the total number of SimHits, in the cluster.

Definition at line 169 of file SimCluster.h.

{ return nsimhits_; }

References nsimhits_.

p()

float SimCluster::p ( ) const

inline

Magnitude of momentum vector. Note this is taken from the first SimTrack only.

Definition at line 101 of file SimCluster.h.

{ return p4().P(); }

References p4().

Referenced by Electron.Electron::ptErr().

p4()

const math::XYZTLorentzVectorF& SimCluster::p4 ( ) const

inline

Four-momentum Lorentz vector. Note this is taken from the first SimTrack only.

Definition at line 91 of file SimCluster.h.

{ return theMomentum_; }

References theMomentum_.

Referenced by boostToCM(), Tau.Tau::dxy_approx(), Tau.Tau::dz(), energy(), et(), eta(), mass(), momentum(), mt(), mtSqr(), p(), Lepton.Lepton::p4WithFSR(), phi(), pt(), px(), py(), pz(), rapidity(), and theta().

particleId()

uint64_t SimCluster::particleId ( ) const

inline

Definition at line 68 of file SimCluster.h.

{ return particleId_; }

References particleId_.

pdgId()

int SimCluster::pdgId ( ) const

inline

PDG ID.

Returns the PDG ID of the first associated gen particle. If there are no gen particles associated then it returns type() from the first SimTrack.

Definition at line 55 of file SimCluster.h.

                    {
    if (genParticles_.empty())
      return g4Tracks_[0].type();
    else
      return (*genParticles_.begin())->pdgId();
  }

References edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::empty(), g4Tracks_, and genParticles_.

Referenced by Particle.Particle::__str__().

phi()

float SimCluster::phi ( void  ) const

inline

Momentum azimuthal angle. Note this is taken from the first SimTrack only.

Definition at line 140 of file SimCluster.h.

{ return p4().Phi(); }

References p4().

Referenced by Particle.Particle::__str__(), and ntupleDataFormat.Track::phiPull().

pt()

float SimCluster::pt ( ) const

inline

Transverse momentum. Note this is taken from the first SimTrack only.

Definition at line 136 of file SimCluster.h.

{ return p4().Pt(); }

References p4().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

px()

float SimCluster::px ( ) const

inline

x coordinate of momentum vector. Note this is taken from the first SimTrack only.

Definition at line 124 of file SimCluster.h.

{ return p4().Px(); }

References p4().

py()

float SimCluster::py ( ) const

inline

y coordinate of momentum vector. Note this is taken from the first SimTrack only.

Definition at line 128 of file SimCluster.h.

{ return p4().Py(); }

References p4().

pz()

float SimCluster::pz ( ) const

inline

z coordinate of momentum vector. Note this is taken from the first SimTrack only.

Definition at line 132 of file SimCluster.h.

{ return p4().Pz(); }

References p4().

rapidity()

float SimCluster::rapidity ( ) const

inline

Rapidity. Note this is taken from the simtrack before the calorimeter.

Definition at line 152 of file SimCluster.h.

{ return p4().Rapidity(); }

References p4().

Referenced by y().

simEnergy()

float SimCluster::simEnergy ( ) const

inline

returns the accumulated sim energy in the cluster

Definition at line 213 of file SimCluster.h.

{ return simhit_energy_; }

References simhit_energy_.

Referenced by HGCalHitCalibration::analyze(), and HGCalShowerSeparation::analyze().

status()

int SimCluster::status ( void  ) const

inline

Status word.

Returns status() from the first gen particle, or -99 if there are no gen particles attached.

Definition at line 161 of file SimCluster.h.

{ return genParticles_.empty() ? -99 : (*genParticles_[0]).status(); }

References edm::RefVector< C, T, F >::empty(), genParticles_, and status().

Referenced by longLived(), and status().

theta()

float SimCluster::theta ( void  ) const

inline

Momentum polar angle. Note this is taken from the first SimTrack only.

Definition at line 144 of file SimCluster.h.

{ return p4().Theta(); }

References p4().

Referenced by Tau.Tau::zImpact().

threeCharge()

int SimCluster::threeCharge ( ) const

inline

Gives charge in unit of quark charge (should be 3 times "charge()")

Definition at line 87 of file SimCluster.h.

{ return lrintf(3.f * charge()); }

References charge(), and f.

y()

float SimCluster::y ( ) const

inline

Same as rapidity().

Definition at line 155 of file SimCluster.h.

{ return rapidity(); }

References rapidity().

Referenced by svgfig.Ellipse::__repr__(), geometryXMLparser.Alignable::pos(), ntupleDataFormat._HitObject::r(), and ntupleDataFormat._HitObject::r3D().

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  s, SimCluster const &  tp  )

friend

Definition at line 31 of file SimCluster.cc.

                                                            {
  s << "CP momentum, q, ID, & Event #: " << tp.p4() << " " << tp.charge() << " " << tp.pdgId() << " "
    << tp.eventId().bunchCrossing() << "." << tp.eventId().event() << std::endl;
 
  for (SimCluster::genp_iterator hepT = tp.genParticle_begin(); hepT != tp.genParticle_end(); ++hepT) {
    s << " HepMC Track Momentum " << (*hepT)->momentum().rho() << std::endl;
  }
 
  for (SimCluster::g4t_iterator g4T = tp.g4Track_begin(); g4T != tp.g4Track_end(); ++g4T) {
    s << " Geant Track Momentum  " << g4T->momentum() << std::endl;
    s << " Geant Track ID & type " << g4T->trackId() << " " << g4T->type() << std::endl;
    if (g4T->type() != tp.pdgId()) {
      s << " Mismatch b/t SimCluster and Geant types" << std::endl;
    }
  }
  s << " # of cells = " << tp.hits_.size()
    << ", effective cells = " << std::accumulate(tp.fractions_.begin(), tp.fractions_.end(), 0.f) << std::endl;
  return s;
}

Member Data Documentation

energies_

std::vector<float> SimCluster::energies_

private

Definition at line 226 of file SimCluster.h.

Referenced by addHitEnergy(), clearHitsEnergy(), and hits_and_energies().

event_

EncodedEventId SimCluster::event_

private

Definition at line 220 of file SimCluster.h.

Referenced by eventId(), and SimCluster().

fractions_

std::vector<float> SimCluster::fractions_

private

Definition at line 225 of file SimCluster.h.

Referenced by addRecHitAndFraction(), clearHitsAndFractions(), and hits_and_fractions().

g4Tracks_

std::vector<SimTrack> SimCluster::g4Tracks_

private

references to G4 and reco::GenParticle tracks

Definition at line 231 of file SimCluster.h.

Referenced by addG4Track(), charge(), g4Track_begin(), g4Track_end(), g4Tracks(), and pdgId().

genParticles_

reco::GenParticleRefVector SimCluster::genParticles_

private

Definition at line 232 of file SimCluster.h.

Referenced by addGenParticle(), genParticle_begin(), genParticle_end(), genParticles(), pdgId(), and status().

hits_

std::vector<uint32_t> SimCluster::hits_

private

Definition at line 224 of file SimCluster.h.

Referenced by addRecHitAndFraction(), clearHitsAndFractions(), hits_and_energies(), hits_and_fractions(), and numberOfRecHits().

longLivedTag

const unsigned int SimCluster::longLivedTag = 65536

static

long lived flag

Definition at line 163 of file SimCluster.h.

Referenced by longLived().

nsimhits_

uint64_t SimCluster::nsimhits_

private

Definition at line 219 of file SimCluster.h.

Referenced by numberOfSimHits().

particleId_

uint32_t SimCluster::particleId_

private

Definition at line 222 of file SimCluster.h.

Referenced by particleId(), and SimCluster().

simhit_energy_

float SimCluster::simhit_energy_

private

Definition at line 223 of file SimCluster.h.

Referenced by addSimHit(), and simEnergy().

theMomentum_

math::XYZTLorentzVectorF SimCluster::theMomentum_

private

Definition at line 228 of file SimCluster.h.

Referenced by p4(), and SimCluster().