|
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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
#ifndef DataFormats_ParticleFlowReco_PFTrajectoryPoint_h
#define DataFormats_ParticleFlowReco_PFTrajectoryPoint_h
/**
*/
#include <vector>
#include <map>
#include <iosfwd>
#include "DataFormats/Math/interface/Point3D.h"
#include "Rtypes.h"
#include "DataFormats/Math/interface/LorentzVector.h"
#include "Math/GenVector/PositionVector3D.h"
namespace reco {
/**\class PFTrajectoryPoint
\brief A PFTrack holds several trajectory points, which basically
contain the position and momentum of a track at a given position.
\todo detId_, layer_, isTrackerLayer_ seem to be redundant
\todo deal with origin and end vertices of PFSimParticles
\todo remove HCAL exit
\author Renaud Bruneliere
\date July 2006
*/
class PFTrajectoryPoint {
public:
// Next typedef uses double in ROOT 6 rather than Double32_t due to a bug in ROOT 5,
// which otherwise would make ROOT5 files unreadable in ROOT6. This does not increase
// the size on disk, because due to the bug, double was actually stored on disk in ROOT 5.
typedef ROOT::Math::PositionVector3D<ROOT::Math::CylindricalEta3D<double> > REPPoint;
/// Define the different layers where the track can be propagated
enum LayerType {
/// Point of closest approach from beam axis (initial point in the case of PFSimParticle)
ClosestApproach = 0,
BeamPipeOrEndVertex = 1,
/// Preshower layer 1
PS1 = 2,
/// Preshower layer 2
PS2 = 3,
/// ECAL front face
ECALEntrance = 4,
/// expected maximum of the shower in ECAL, for an e/gamma particle
/// \todo: add an ECALShowerMaxHadrons
ECALShowerMax = 5,
/// HCAL front face
HCALEntrance = 6,
/// HCAL exit
HCALExit = 7,
/// HO layer
HOLayer = 8,
/// VFcal(HF) front face
VFcalEntrance = 9,
// Number of valid layers
NLayers = 10,
// Unknown
Unknown = -1
};
static const std::array<std::string, NLayers> layerTypeNames;
static LayerType layerTypeByName(const std::string& name);
/// default constructor. Set variables at default dummy values
PFTrajectoryPoint();
/// \brief constructor from values.
/// set detId to -1 if this point is not from a tracker layer
PFTrajectoryPoint(int detId, int layer, const math::XYZPoint& posxyz, const math::XYZTLorentzVector& momentum);
/// measurement detId
int detId() const { return detId_; }
/// trajectory point layer
int layer() const { return layer_; }
/// is this point valid ?
bool isValid() const {
if (layer_ == -1 && detId_ == -1)
return false;
else
return true;
}
/// is this point corresponding to an intersection with a tracker layer ?
bool isTrackerLayer() const {
if (detId_ >= 0)
return true;
else
return false;
}
/// cartesian position (x, y, z)
const math::XYZPoint& position() const { return posxyz_; }
/// trajectory position in (rho, eta, phi) base
const REPPoint& positionREP() const { return posrep_; }
/// calculate posrep_ once and for all
void calculatePositionREP() { posrep_.SetCoordinates(posxyz_.Rho(), posxyz_.Eta(), posxyz_.Phi()); }
/// 4-momenta quadrivector
const math::XYZTLorentzVector& momentum() const { return momentum_; }
bool operator==(const reco::PFTrajectoryPoint& other) const;
friend std::ostream& operator<<(std::ostream& out, const reco::PFTrajectoryPoint& trajPoint);
private:
/// \brief Is the measurement corresponding to a tracker layer?
/// or was it obtained by propagating the track to a certain position?
bool isTrackerLayer_;
/// detid if measurement is corresponding to a tracker layer
int detId_;
/// propagated layer
int layer_;
/// cartesian position (x, y, z)
math::XYZPoint posxyz_;
/// position in (rho, eta, phi) base (transient)
REPPoint posrep_;
/// momentum quadrivector
math::XYZTLorentzVector momentum_;
};
std::ostream& operator<<(std::ostream& out, const reco::PFTrajectoryPoint& trajPoint);
} // namespace reco
#endif
|