|
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
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
|
#ifndef __DataFormats_PatCandidates_IsolatedTrack_h__
#define __DataFormats_PatCandidates_IsolatedTrack_h__
/*
\class pat::IsolatedTrack IsolatedTrack.h "DataFormats/PatCandidates/interface/IsolatedTrack.h"
\brief Small class to store key info on isolated tracks
pat::IsolatedTrack stores important info on isolated tracks. Draws from
packedPFCandidates, lostTracks, and generalTracks.
\author Bennett Marsh
*/
#include "DataFormats/Candidate/interface/LeafCandidate.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackBase.h"
#include "DataFormats/TrackReco/interface/HitPattern.h"
#include "DataFormats/PatCandidates/interface/PackedCandidate.h"
#include "DataFormats/PatCandidates/interface/PFIsolation.h"
namespace pat {
class IsolatedTrack : public reco::LeafCandidate {
public:
IsolatedTrack()
: LeafCandidate(0, LorentzVector(0, 0, 0, 0)),
pfIsolationDR03_(pat::PFIsolation()),
miniIsolation_(pat::PFIsolation()),
matchedCaloJetEmEnergy_(0.),
matchedCaloJetHadEnergy_(0.),
pfLepOverlap_(false),
pfNeutralSum_(0.),
dz_(0.),
dxy_(0.),
dzError_(0.),
dxyError_(0.),
fromPV_(-1),
trackQuality_(0),
dEdxStrip_(0),
dEdxPixel_(0),
hitPattern_(reco::HitPattern()),
crossedEcalStatus_(std::vector<uint16_t>()),
crossedHcalStatus_(std::vector<uint32_t>()),
deltaEta_(0),
deltaPhi_(0),
packedCandRef_(PackedCandidateRef()),
nearestPFPackedCandRef_(PackedCandidateRef()),
nearestLostTrackPackedCandRef_(PackedCandidateRef()) {}
explicit IsolatedTrack(const PFIsolation& iso,
const PFIsolation& miniiso,
float caloJetEm,
float caloJetHad,
bool pfLepOverlap,
float pfNeutralSum,
const LorentzVector& p4,
int charge,
int id,
float dz,
float dxy,
float dzError,
float dxyError,
const reco::HitPattern& hp,
float dEdxS,
float dEdxP,
int fromPV,
int tkQual,
const std::vector<uint16_t>& ecalst,
const std::vector<uint32_t>& hcalst,
int dEta,
int dPhi,
const PackedCandidateRef& pcref,
const PackedCandidateRef& refToNearestPF,
const PackedCandidateRef& refToNearestLostTrack)
: LeafCandidate(charge, p4, Point(0., 0., 0.), id),
pfIsolationDR03_(iso),
miniIsolation_(miniiso),
matchedCaloJetEmEnergy_(caloJetEm),
matchedCaloJetHadEnergy_(caloJetHad),
pfLepOverlap_(pfLepOverlap),
pfNeutralSum_(pfNeutralSum),
dz_(dz),
dxy_(dxy),
dzError_(dzError),
dxyError_(dxyError),
fromPV_(fromPV),
trackQuality_(tkQual),
dEdxStrip_(dEdxS),
dEdxPixel_(dEdxP),
hitPattern_(hp),
crossedEcalStatus_(ecalst),
crossedHcalStatus_(hcalst),
deltaEta_(dEta),
deltaPhi_(dPhi),
packedCandRef_(pcref),
nearestPFPackedCandRef_(refToNearestPF),
nearestLostTrackPackedCandRef_(refToNearestLostTrack) {}
~IsolatedTrack() override {}
const PFIsolation& pfIsolationDR03() const { return pfIsolationDR03_; }
const PFIsolation& miniPFIsolation() const { return miniIsolation_; }
float matchedCaloJetEmEnergy() const { return matchedCaloJetEmEnergy_; }
float matchedCaloJetHadEnergy() const { return matchedCaloJetHadEnergy_; }
bool pfLepOverlap() const { return pfLepOverlap_; }
float pfNeutralSum() const { return pfNeutralSum_; }
float dz() const { return dz_; }
float dzError() const override { return dzError_; }
float dxy() const { return dxy_; }
float dxyError() const override { return dxyError_; }
int fromPV() const { return fromPV_; }
bool isHighPurityTrack() const {
return (trackQuality_ & (1 << reco::TrackBase::highPurity)) >> reco::TrackBase::highPurity;
}
bool isTightTrack() const { return (trackQuality_ & (1 << reco::TrackBase::tight)) >> reco::TrackBase::tight; }
bool isLooseTrack() const { return (trackQuality_ & (1 << reco::TrackBase::loose)) >> reco::TrackBase::loose; }
const reco::HitPattern& hitPattern() const { return hitPattern_; }
/// helper functions for string cut parser
int lostInnerLayers() const {
return hitPattern_.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS);
}
int lostLayers() const { return hitPattern_.trackerLayersWithoutMeasurement(reco::HitPattern::TRACK_HITS); }
int lostOuterLayers() const {
return hitPattern_.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS);
}
float dEdxStrip() const { return dEdxStrip_; }
float dEdxPixel() const { return dEdxPixel_; }
//! just the status code part of an EcalChannelStatusCode for all crossed Ecal cells
const std::vector<uint16_t>& crossedEcalStatus() const { return crossedEcalStatus_; }
//! just the status code part of an HcalChannelStatus for all crossed Hcal cells
const std::vector<uint32_t>& crossedHcalStatus() const { return crossedHcalStatus_; }
//! difference in eta/phi between initial traj and intersection w/ ecal
//! Values are between +-0.5 with a precision of 0.002
float deltaEta() const { return float(deltaEta_) / 500.f; }
float deltaPhi() const { return float(deltaPhi_) / 500.f; }
const PackedCandidateRef& packedCandRef() const { return packedCandRef_; }
const PackedCandidateRef& nearestPFPackedCandRef() const { return nearestPFPackedCandRef_; }
const PackedCandidateRef& nearestLostTrackPackedCandRef() const { return nearestPFPackedCandRef_; }
protected:
PFIsolation pfIsolationDR03_;
PFIsolation miniIsolation_;
float matchedCaloJetEmEnergy_; //energy of nearest calojet within a given dR;
float matchedCaloJetHadEnergy_;
bool pfLepOverlap_;
float pfNeutralSum_;
float dz_, dxy_, dzError_, dxyError_;
int fromPV_; //only stored for packedPFCandidates
int trackQuality_;
float dEdxStrip_, dEdxPixel_; //in MeV/mm
reco::HitPattern hitPattern_;
std::vector<uint16_t> crossedEcalStatus_;
std::vector<uint32_t> crossedHcalStatus_;
int deltaEta_, deltaPhi_;
PackedCandidateRef packedCandRef_; // stored only for packedPFCands/lostTracks. NULL for generalTracks
PackedCandidateRef nearestPFPackedCandRef_;
PackedCandidateRef nearestLostTrackPackedCandRef_;
};
typedef std::vector<IsolatedTrack> IsolatedTrackCollection;
} // namespace pat
#endif
|