CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Attributes
G4SimEvent Class Reference

#include <G4SimEvent.h>

Public Member Functions

void add (G4SimTrack *t)
 
void add (G4SimVertex *v)
 
void collisionPoint (const math::XYZTLorentzVectorD &v)
 
const math::XYZTLorentzVectorDcollisionPoint () const
 
 G4SimEvent ()
 
const G4SimTrackg4track (int i) const
 
const G4SimVertexg4vertex (int i) const
 
void hepEvent (const HepMC::GenEvent *r)
 
const HepMC::GenEventhepEvent () const
 
void load (edm::SimTrackContainer &c) const
 
void load (edm::SimVertexContainer &c) const
 
unsigned int nGenParts () const
 
void nparam (int n)
 
const int nparam () const
 
unsigned int nTracks () const
 
unsigned int nVertices () const
 
void param (const std::vector< float > &p)
 
const std::vector< float > & param () const
 
void weight (float w)
 
float weight () const
 
virtual ~G4SimEvent ()
 

Protected Attributes

math::XYZTLorentzVectorD collisionPoint_
 
std::vector< G4SimTrack * > g4tracks
 
std::vector< G4SimVertex * > g4vertices
 
const HepMC::GenEventhepMCEvent
 
int nparam_
 
std::vector< float > param_
 
float weight_
 

Detailed Description

Definition at line 14 of file G4SimEvent.h.

Constructor & Destructor Documentation

G4SimEvent::G4SimEvent ( )

Definition at line 11 of file G4SimEvent.cc.

12  : hepMCEvent(nullptr),
13  weight_(0),
15  nparam_(0),
16  param_(0) {}
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:14
int nparam_
Definition: G4SimEvent.h:42
float weight_
Definition: G4SimEvent.h:40
const HepMC::GenEvent * hepMCEvent
Definition: G4SimEvent.h:39
std::vector< float > param_
Definition: G4SimEvent.h:43
math::XYZTLorentzVectorD collisionPoint_
Definition: G4SimEvent.h:41
G4SimEvent::~G4SimEvent ( )
virtual

Definition at line 18 of file G4SimEvent.cc.

References g4tracks, g4vertices, and mps_fire::i.

18  {
19  // per suggestion by Chris Jones, it's faster
20  // that delete back() and pop_back()
21  //
22  unsigned int i = 0;
23 
24  for (i = 0; i < g4tracks.size(); i++) {
25  delete g4tracks[i];
26  g4tracks[i] = nullptr;
27  }
28  g4tracks.clear();
29 
30  for (i = 0; i < g4vertices.size(); i++) {
31  delete g4vertices[i];
32  g4vertices[i] = nullptr;
33  }
34  g4vertices.clear();
35 }
std::vector< G4SimTrack * > g4tracks
Definition: G4SimEvent.h:44
std::vector< G4SimVertex * > g4vertices
Definition: G4SimEvent.h:45

Member Function Documentation

void G4SimEvent::add ( G4SimTrack t)
inline

Definition at line 33 of file G4SimEvent.h.

References g4tracks.

Referenced by SimTrackManager::getOrCreateVertex(), SimTrackManager::reallyStoreTracks(), and counter.Counter::register().

33 { g4tracks.push_back(t); }
std::vector< G4SimTrack * > g4tracks
Definition: G4SimEvent.h:44
void G4SimEvent::add ( G4SimVertex v)
inline

Definition at line 34 of file G4SimEvent.h.

References g4vertices.

Referenced by counter.Counter::register().

34 { g4vertices.push_back(v); }
std::vector< G4SimVertex * > g4vertices
Definition: G4SimEvent.h:45
void G4SimEvent::collisionPoint ( const math::XYZTLorentzVectorD v)
inline

Definition at line 27 of file G4SimEvent.h.

References collisionPoint_, and findQualityFiles::v.

Referenced by RunManager::produce().

27 { collisionPoint_ = v; }
math::XYZTLorentzVectorD collisionPoint_
Definition: G4SimEvent.h:41
const math::XYZTLorentzVectorD& G4SimEvent::collisionPoint ( ) const
inline

Definition at line 28 of file G4SimEvent.h.

References collisionPoint_.

28 { return collisionPoint_; }
math::XYZTLorentzVectorD collisionPoint_
Definition: G4SimEvent.h:41
const G4SimTrack& G4SimEvent::g4track ( int  i) const
inline

Definition at line 35 of file G4SimEvent.h.

References g4tracks.

35 { return *g4tracks[i - 1]; }
std::vector< G4SimTrack * > g4tracks
Definition: G4SimEvent.h:44
const G4SimVertex& G4SimEvent::g4vertex ( int  i) const
inline

Definition at line 36 of file G4SimEvent.h.

References g4vertices.

36 { return *g4vertices[i - 1]; }
std::vector< G4SimVertex * > g4vertices
Definition: G4SimEvent.h:45
void G4SimEvent::hepEvent ( const HepMC::GenEvent r)
inline

Definition at line 23 of file G4SimEvent.h.

References hepMCEvent, and alignCSCRings::r.

Referenced by RunManagerMTWorker::produce(), and RunManager::produce().

23 { hepMCEvent = r; }
const HepMC::GenEvent * hepMCEvent
Definition: G4SimEvent.h:39
const HepMC::GenEvent* G4SimEvent::hepEvent ( ) const
inline

Definition at line 24 of file G4SimEvent.h.

References hepMCEvent.

24 { return hepMCEvent; }
const HepMC::GenEvent * hepMCEvent
Definition: G4SimEvent.h:39
void G4SimEvent::load ( edm::SimTrackContainer c) const

Definition at line 37 of file G4SimEvent.cc.

References G4SimTrack::energy(), g4tracks, GeV, mps_fire::i, G4SimTrack::id(), G4SimTrack::igenpart(), G4SimTrack::ivert(), G4SimTrack::momentum(), AlCaHLTBitMon_ParallelJobs::p, G4SimTrack::part(), mcMuonSeeds_cfi::SimTrack, protons_cff::t, G4SimTrack::trackerSurfaceMomentum(), and G4SimTrack::trackerSurfacePosition().

Referenced by OscarProducer::produce().

37  {
38  for (unsigned int i = 0; i < g4tracks.size(); i++) {
39  G4SimTrack* trk = g4tracks[i];
40  int ip = trk->part();
42  trk->momentum().x() / GeV, trk->momentum().y() / GeV, trk->momentum().z() / GeV, trk->energy() / GeV);
43  int iv = trk->ivert();
44  int ig = trk->igenpart();
45  int id = trk->id();
46  math::XYZVectorD tkpos(trk->trackerSurfacePosition().x() / cm,
47  trk->trackerSurfacePosition().y() / cm,
48  trk->trackerSurfacePosition().z() / cm);
50  trk->trackerSurfaceMomentum().y() / GeV,
51  trk->trackerSurfaceMomentum().z() / GeV,
52  trk->trackerSurfaceMomentum().e() / GeV);
53  // ip = particle ID as PDG
54  // pp = 4-momentum
55  // iv = corresponding G4SimVertex index
56  // ig = corresponding GenParticle index
57  SimTrack t = SimTrack(ip, p, iv, ig, tkpos, tkmom);
58  t.setTrackId(id);
59  t.setEventId(EncodedEventId(0));
60  c.push_back(t);
61  }
62  std::stable_sort(c.begin(), c.end(), IdSort());
63 }
double energy() const
Definition: G4SimTrack.h:60
const double GeV
Definition: MathUtil.h:16
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:14
int part() const
Definition: G4SimTrack.h:58
std::vector< G4SimTrack * > g4tracks
Definition: G4SimEvent.h:44
int ivert() const
Definition: G4SimTrack.h:61
const math::XYZVectorD & trackerSurfacePosition() const
Definition: G4SimTrack.h:66
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > XYZVectorD
spatial vector with cartesian internal representation
Definition: Vector3D.h:8
int igenpart() const
Definition: G4SimTrack.h:62
int id() const
Definition: G4SimTrack.h:57
const math::XYZTLorentzVectorD & trackerSurfaceMomentum() const
Definition: G4SimTrack.h:67
const math::XYZVectorD & momentum() const
Definition: G4SimTrack.h:59
void G4SimEvent::load ( edm::SimVertexContainer c) const

Definition at line 65 of file G4SimEvent.cc.

References g4vertices, mps_fire::i, G4SimVertex::parentIndex(), G4SimVertex::processType(), edm::second(), CoreSimVertex::setEventId(), SimVertex::setProcessType(), mcMuonSeeds_cfi::SimVertex, protons_cff::t, findQualityFiles::v, G4SimVertex::vertexGlobalTime(), G4SimVertex::vertexPosition(), and extraflags_cff::vtx.

65  {
66  for (unsigned int i = 0; i < g4vertices.size(); i++) {
68  //
69  // starting 1_1_0_pre3, SimVertex stores in cm !!!
70  //
71  math::XYZVectorD v3(vtx->vertexPosition().x() / cm, vtx->vertexPosition().y() / cm, vtx->vertexPosition().z() / cm);
72  float t = vtx->vertexGlobalTime() / second;
73  int iv = vtx->parentIndex();
74  // vv = position
75  // t = global time
76  // iv = index of the parent in the SimEvent SimTrack container (-1 if no parent)
77  SimVertex v = SimVertex(v3, t, iv, i);
78  v.setProcessType(vtx->processType());
80  c.push_back(v);
81  }
82 }
int parentIndex() const
Definition: G4SimVertex.h:17
void setProcessType(unsigned int ty)
Definition: SimVertex.h:35
U second(std::pair< T, U > const &p)
unsigned int processType() const
Definition: G4SimVertex.h:18
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > XYZVectorD
spatial vector with cartesian internal representation
Definition: Vector3D.h:8
const math::XYZVectorD & vertexPosition() const
index of the parent (-1 if no parent)
Definition: G4SimVertex.h:15
std::vector< G4SimVertex * > g4vertices
Definition: G4SimEvent.h:45
double vertexGlobalTime() const
Definition: G4SimVertex.h:16
void setEventId(EncodedEventId e)
Definition: CoreSimVertex.h:23
unsigned int G4SimEvent::nGenParts ( ) const
inline

Definition at line 22 of file G4SimEvent.h.

References hepMCEvent.

Referenced by RunManager::produce().

22 { return hepMCEvent->particles_size(); }
const HepMC::GenEvent * hepMCEvent
Definition: G4SimEvent.h:39
void G4SimEvent::nparam ( int  n)
inline

Definition at line 29 of file G4SimEvent.h.

References gen::n, and nparam_.

29 { nparam_ = n; }
int nparam_
Definition: G4SimEvent.h:42
const int G4SimEvent::nparam ( ) const
inline

Definition at line 30 of file G4SimEvent.h.

References nparam_.

30 { return nparam_; }
int nparam_
Definition: G4SimEvent.h:42
unsigned int G4SimEvent::nTracks ( ) const
inline

Definition at line 20 of file G4SimEvent.h.

References g4tracks.

Referenced by RunManager::produce().

20 { return g4tracks.size(); }
std::vector< G4SimTrack * > g4tracks
Definition: G4SimEvent.h:44
unsigned int G4SimEvent::nVertices ( ) const
inline

Definition at line 21 of file G4SimEvent.h.

References g4vertices.

Referenced by RunManager::produce().

21 { return g4vertices.size(); }
std::vector< G4SimVertex * > g4vertices
Definition: G4SimEvent.h:45
void G4SimEvent::param ( const std::vector< float > &  p)
inline

Definition at line 31 of file G4SimEvent.h.

References AlCaHLTBitMon_ParallelJobs::p, and param_.

const std::vector<float>& G4SimEvent::param ( ) const
inline

Definition at line 32 of file G4SimEvent.h.

References param_.

32 { return param_; }
std::vector< float > param_
Definition: G4SimEvent.h:43
void G4SimEvent::weight ( float  w)
inline

Definition at line 25 of file G4SimEvent.h.

References w, and weight_.

Referenced by RunManager::produce().

25 { weight_ = w; }
const double w
Definition: UKUtility.cc:23
float weight_
Definition: G4SimEvent.h:40
float G4SimEvent::weight ( ) const
inline

Definition at line 26 of file G4SimEvent.h.

References weight_.

26 { return weight_; }
float weight_
Definition: G4SimEvent.h:40

Member Data Documentation

math::XYZTLorentzVectorD G4SimEvent::collisionPoint_
protected

Definition at line 41 of file G4SimEvent.h.

Referenced by collisionPoint().

std::vector<G4SimTrack*> G4SimEvent::g4tracks
protected

Definition at line 44 of file G4SimEvent.h.

Referenced by add(), g4track(), load(), nTracks(), and ~G4SimEvent().

std::vector<G4SimVertex*> G4SimEvent::g4vertices
protected

Definition at line 45 of file G4SimEvent.h.

Referenced by add(), g4vertex(), load(), nVertices(), and ~G4SimEvent().

const HepMC::GenEvent* G4SimEvent::hepMCEvent
protected

Definition at line 39 of file G4SimEvent.h.

Referenced by hepEvent(), and nGenParts().

int G4SimEvent::nparam_
protected

Definition at line 42 of file G4SimEvent.h.

Referenced by nparam().

std::vector<float> G4SimEvent::param_
protected

Definition at line 43 of file G4SimEvent.h.

Referenced by param().

float G4SimEvent::weight_
protected

Definition at line 40 of file G4SimEvent.h.

Referenced by weight().