|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
#include "FastSimulation/Event/interface/FBaseSimEvent.h"
#include "FastSimulation/Event/interface/FSimVertex.h"
inline const FSimVertex FSimTrack::vertex() const {
return (!mom_ ? FSimVertex(vertex_.position(), vertex_.parentIndex(), vertex_.vertexId(), nullptr)
: mom_->vertex(vertIndex()));
}
inline const FSimVertex& FSimTrack::endVertex() const {
if (endv_ < 0)
throw cms::Exception("FastSim")
<< "FSimTrack::endVertex() called for FSimTrack w/o end vertex, please contact FastSim developers" << std::endl;
return mom_->vertex(endv_);
}
inline const FSimTrack& FSimTrack::mother() const {
if (noMother()) {
throw cms::Exception("FastSim")
<< "FSimTrack::mother() called for FSimTrack w/o mother, please contact FastSim developers" << std::endl;
}
return vertex().parent();
}
inline const FSimTrack& FSimTrack::daughter(int i) const {
if (i < 0 || i >= nDaughters())
throw cms::Exception("FastSim")
<< "FSimTrack::daughter(int index) index out of range, please contact FastSim developers" << std::endl;
if (abs(type()) == 11 || abs(type()) == 13)
return mom_->track(daugh_[i]);
else
return endVertex().daughter(i);
}
inline int FSimTrack::nDaughters() const {
if (abs(type()) == 11 || abs(type()) == 13)
return daugh_.size();
else {
if (noEndVertex())
return 0;
else
return endVertex().nDaughters();
}
}
inline const std::vector<int>& FSimTrack::daughters() const {
if (abs(type()) == 11 || noEndVertex())
return daugh_;
else {
//if(noEndVertex()){
//throw cms::Exception("FastSim") << "FSimTrack::daughters() called for FSimTrack w/o end vertex, please contact FastSim developers" << std::endl;
//}
return endVertex().daughters();
}
}
inline bool FSimTrack::noEndVertex() const {
// The particle either has no end vertex index
if (endv_ < 0)
return true;
// or it's an electron/muon that has just Brem'ed, but continues its way
// ... but not those intermediate e/mu PYTHIA entries with prompt Brem
bool bremOutOfPipe = true;
if ((mom_->vertex(endv_)).position().Perp2() < 1.0)
bremOutOfPipe = false;
return ((abs(type()) == 11 || abs(type()) == 13) && bremOutOfPipe && endVertex().nDaughters() > 0 &&
endVertex().daughter(endVertex().nDaughters() - 1).type() == 22);
}
inline bool FSimTrack::noMother() const { return noVertex() || vertex().noParent(); }
inline bool FSimTrack::noDaughter() const { return noEndVertex() || !nDaughters(); }
inline const HepMC::GenParticle* FSimTrack::genParticle() const { return mom_->embdGenpart(genpartIndex()); }
|