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_ {0}
uint32_t	particleId_ {0}
float	simhit_energy_ {0.f}
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

typedef int SimCluster::Charge

electric charge type

Definition at line 33 of file SimCluster.h.

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

Definition at line 41 of file SimCluster.h.

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

reference to reco::GenParticle

Definition at line 40 of file SimCluster.h.

typedef math::XYZTLorentzVectorD SimCluster::LorentzVector

Lorentz vector.

Definition at line 34 of file SimCluster.h.

typedef math::XYZPointD SimCluster::Point

point in the space

Definition at line 36 of file SimCluster.h.

typedef math::PtEtaPhiMLorentzVector SimCluster::PolarLorentzVector

Lorentz vector.

Definition at line 35 of file SimCluster.h.

typedef math::XYZVectorD SimCluster::Vector

point in the space

Definition at line 37 of file SimCluster.h.

Constructor & Destructor Documentation

SimCluster::SimCluster

(

)

Definition at line 11 of file SimCluster.cc.

                       {
  // No operation
}

SimCluster::SimCluster

(

const SimTrack &

simtrk

)

Definition at line 15 of file SimCluster.cc.

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

                                             {
  addG4Track(simtrk);
  event_ = simtrk.eventId();
  particleId_ = simtrk.trackId();

  theMomentum_.SetPxPyPzE(
      simtrk.momentum().px(), simtrk.momentum().py(), simtrk.momentum().pz(), simtrk.momentum().E());
}

SimCluster::SimCluster	(	EncodedEventId	eventID,
		uint32_t	particleID
	)

Definition at line 24 of file SimCluster.cc.

References event_, and particleId_.

                                                                  {
  event_ = eventID;
  particleId_ = particleID;
}

SimCluster::~SimCluster

(

)

Definition at line 29 of file SimCluster.cc.

{}

Member Function Documentation

void SimCluster::addG4Track ( const SimTrack &  t )

inline

Definition at line 72 of file SimCluster.h.

References g4Tracks_.

Referenced by SimCluster().

{ g4Tracks_.push_back(t); }

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

inline

Definition at line 71 of file SimCluster.h.

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

{ genParticles_.push_back(ref); }

void SimCluster::addHitEnergy ( float  energy )

inline

add rechit energy

Definition at line 181 of file SimCluster.h.

References energies_.

{ energies_.emplace_back(energy); }

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

inline

add rechit with fraction

Definition at line 175 of file SimCluster.h.

References fractions_, and hits_.

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

void SimCluster::addSimHit ( const PCaloHit &  hit )

inline

add simhit's energy to cluster

Definition at line 216 of file SimCluster.h.

References PCaloHit::energy(), nsimhits_, and simhit_energy_.

                                      {
    simhit_energy_ += hit.energy();
    ++nsimhits_;
  }

math::XYZVectorF SimCluster::boostToCM ( ) const

inline

Vector to boost to the particle centre of mass frame.

Definition at line 97 of file SimCluster.h.

References p4().

{ return p4().BoostToCM(); }

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.

References g4Tracks_.

Referenced by operator<<(), and threeCharge().

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

void SimCluster::clearHitsAndFractions ( )

inline

clear the hits and fractions list

Definition at line 204 of file SimCluster.h.

References fractions_, hits_, and edm::swap().

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

void SimCluster::clearHitsEnergy ( )

inline

clear the energies list

Definition at line 210 of file SimCluster.h.

References energies_, and edm::swap().

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

float SimCluster::energy ( void  ) const

inline

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

Definition at line 104 of file SimCluster.h.

References p4().

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

{ return p4().E(); }

float SimCluster::et ( ) const

inline

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

Definition at line 107 of file SimCluster.h.

References p4().

{ return p4().Et(); }

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.

References p4().

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

{ return p4().Eta(); }

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.

References event_.

Referenced by ntupleDataFormat.Event::eventIdStr(), and operator<<().

{ return event_; }

g4t_iterator SimCluster::g4Track_begin ( ) const

inline

Definition at line 76 of file SimCluster.h.

References g4Tracks_.

Referenced by operator<<().

{ return g4Tracks_.begin(); }

g4t_iterator SimCluster::g4Track_end ( ) const

inline

Definition at line 77 of file SimCluster.h.

References g4Tracks_.

Referenced by operator<<().

{ return g4Tracks_.end(); }

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

inline

Definition at line 82 of file SimCluster.h.

References g4Tracks_.

{ return g4Tracks_; }

genp_iterator SimCluster::genParticle_begin ( ) const

inline

iterators

Definition at line 74 of file SimCluster.h.

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

Referenced by operator<<().

{ return genParticles_.begin(); }

genp_iterator SimCluster::genParticle_end ( ) const

inline

Definition at line 75 of file SimCluster.h.

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

Referenced by operator<<().

{ return genParticles_.end(); }

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

inline

Definition at line 80 of file SimCluster.h.

References genParticles_.

{ return genParticles_; }

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.

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

                                                                {
    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;
  }

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.

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

Referenced by HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), HGVHistoProducerAlgo::layerClusters_to_CaloParticles(), MultiClusterAssociatorByEnergyScoreImpl::makeConnections(), LCToCPAssociatorByEnergyScoreImpl::makeConnections(), and HGVHistoProducerAlgo::tracksters_to_SimTracksters().

                                                                 {
    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;
  }

bool SimCluster::longLived ( ) const

inline

is long lived?

Definition at line 166 of file SimCluster.h.

References longLivedTag, and status().

{ return status() & longLivedTag; }

float SimCluster::mass ( ) const

inline

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

Definition at line 110 of file SimCluster.h.

References p4().

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

{ return p4().M(); }

float SimCluster::massSqr ( ) const

inline

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

Definition at line 113 of file SimCluster.h.

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

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

math::XYZVectorF SimCluster::momentum ( ) const

inline

spatial momentum vector

Definition at line 94 of file SimCluster.h.

References p4().

{ return p4().Vect(); }

float SimCluster::mt ( ) const

inline

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

Definition at line 116 of file SimCluster.h.

References p4().

{ return p4().Mt(); }

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.

References p4().

{ return p4().Mt2(); }

int SimCluster::numberOfRecHits ( void  ) const

inline

Gives the total number of SimHits, in the cluster.

Definition at line 172 of file SimCluster.h.

References hits_.

{ return hits_.size(); }

int SimCluster::numberOfSimHits ( ) const

inline

Gives the total number of SimHits, in the cluster.

Definition at line 169 of file SimCluster.h.

References nsimhits_.

{ return nsimhits_; }

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.

References p4().

Referenced by Electron.Electron::ptErr().

{ return p4().P(); }

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.

References theMomentum_.

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

{ return theMomentum_; }

uint64_t SimCluster::particleId ( ) const

inline

Definition at line 68 of file SimCluster.h.

References particleId_.

{ return particleId_; }

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.

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

Referenced by Particle.Particle::__str__(), and operator<<().

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

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.

References p4().

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

{ return p4().Phi(); }

float SimCluster::pt ( ) const

inline

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

Definition at line 136 of file SimCluster.h.

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

{ return p4().Pt(); }

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.

References p4().

{ return p4().Px(); }

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.

References p4().

{ return p4().Py(); }

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.

References p4().

{ return p4().Pz(); }

float SimCluster::rapidity ( ) const

inline

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

Definition at line 152 of file SimCluster.h.

References p4().

Referenced by y().

{ return p4().Rapidity(); }

float SimCluster::simEnergy ( ) const

inline

returns the accumulated sim energy in the cluster

Definition at line 213 of file SimCluster.h.

References simhit_energy_.

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

{ return simhit_energy_; }

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.

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

Referenced by longLived(), and status().

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

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.

References p4().

Referenced by Tau.Tau::zImpact().

{ return p4().Theta(); }

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.

References charge(), and validate-o2o-wbm::f.

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

float SimCluster::y ( ) const

inline

Same as rapidity().

Definition at line 155 of file SimCluster.h.

References rapidity().

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

{ return rapidity(); }

Friends And Related Function Documentation

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

std::vector<float> SimCluster::energies_

private

Definition at line 229 of file SimCluster.h.

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

EncodedEventId SimCluster::event_

private

Definition at line 223 of file SimCluster.h.

Referenced by eventId(), and SimCluster().

std::vector<float> SimCluster::fractions_

private

Definition at line 228 of file SimCluster.h.

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

std::vector<SimTrack> SimCluster::g4Tracks_

private

references to G4 and reco::GenParticle tracks

Definition at line 234 of file SimCluster.h.

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

reco::GenParticleRefVector SimCluster::genParticles_

private

Definition at line 235 of file SimCluster.h.

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

std::vector<uint32_t> SimCluster::hits_

private

Definition at line 227 of file SimCluster.h.

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

const unsigned int SimCluster::longLivedTag = 65536

static

long lived flag

Definition at line 163 of file SimCluster.h.

Referenced by longLived().

uint64_t SimCluster::nsimhits_ {0}

private

Definition at line 222 of file SimCluster.h.

Referenced by addSimHit(), and numberOfSimHits().

uint32_t SimCluster::particleId_ {0}

private

Definition at line 225 of file SimCluster.h.

Referenced by particleId(), and SimCluster().

float SimCluster::simhit_energy_ {0.f}

private

Definition at line 226 of file SimCluster.h.

Referenced by addSimHit(), and simEnergy().

math::XYZTLorentzVectorF SimCluster::theMomentum_

private

Definition at line 231 of file SimCluster.h.

Referenced by p4(), and SimCluster().