CMS 3D CMS Logo
PFTrackAlgoTools.cc
Go to the documentation of this file.
1 #include "RecoParticleFlow/PFTracking/interface/PFTrackAlgoTools.h"
2 
3 namespace PFTrackAlgoTools {
4 
5  double dPtCut(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<double>& cuts,bool hltIterativeTracking = true){
6  switch (algo) {
7  case reco::TrackBase::ctf:
8  case reco::TrackBase::duplicateMerge:
9  case reco::TrackBase::initialStep:
10  case reco::TrackBase::highPtTripletStep:
11  case reco::TrackBase::lowPtQuadStep:
12  case reco::TrackBase::lowPtTripletStep:
13  case reco::TrackBase::pixelPairStep:
14  case reco::TrackBase::jetCoreRegionalStep:
15  return cuts[0];
16  case reco::TrackBase::detachedQuadStep:
17  case reco::TrackBase::detachedTripletStep:
18  return cuts[1];
19  case reco::TrackBase::mixedTripletStep:
20  return cuts[2];
21  case reco::TrackBase::pixelLessStep:
22  return cuts[3];
23  case reco::TrackBase::tobTecStep:
24  return cuts[4];
25  case reco::TrackBase::muonSeededStepInOut:
26  case reco::TrackBase::muonSeededStepOutIn:
27  return cuts.at(5);
28  case reco::TrackBase::hltIter0:
29  case reco::TrackBase::hltIter1:
30  case reco::TrackBase::hltIter2:
31  return cuts[0];
32  case reco::TrackBase::hltIter3:
33  return hltIterativeTracking ? cuts[1] : cuts[0];
34  case reco::TrackBase::hltIter4:
35  return hltIterativeTracking ? cuts[2] : cuts[0];
36  case reco::TrackBase::hltIterX:
37  return cuts[0];
38  default:
39  return hltIterativeTracking ? cuts.at(5):cuts[0];
40 
41  }
42  }
43 
44 
45 
46  unsigned int nHitCut(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<unsigned int>& cuts,bool hltIterativeTracking = true){
47  switch (algo) {
48  case reco::TrackBase::ctf:
49  case reco::TrackBase::duplicateMerge:
50  case reco::TrackBase::initialStep:
51  case reco::TrackBase::highPtTripletStep:
52  case reco::TrackBase::lowPtQuadStep:
53  case reco::TrackBase::lowPtTripletStep:
54  case reco::TrackBase::pixelPairStep:
55  case reco::TrackBase::jetCoreRegionalStep:
56  return cuts[0];
57  case reco::TrackBase::detachedQuadStep:
58  case reco::TrackBase::detachedTripletStep:
59  return cuts[1];
60  case reco::TrackBase::mixedTripletStep:
61  return cuts[2];
62  case reco::TrackBase::pixelLessStep:
63  return cuts[3];
64  case reco::TrackBase::tobTecStep:
65  return cuts[4];
66  case reco::TrackBase::muonSeededStepInOut:
67  case reco::TrackBase::muonSeededStepOutIn:
68  return cuts.at(5);
69  case reco::TrackBase::hltIter0:
70  case reco::TrackBase::hltIter1:
71  case reco::TrackBase::hltIter2:
72  return cuts[0];
73  case reco::TrackBase::hltIter3:
74  return hltIterativeTracking ? cuts[1] : cuts[0];
75  case reco::TrackBase::hltIter4:
76  return hltIterativeTracking ? cuts[2] : cuts[0];
77  case reco::TrackBase::hltIterX:
78  return cuts[0];
79  default:
80  return hltIterativeTracking ? cuts.at(5):cuts[0];
81 
82  }
83  }
84 
85 
86 
87  double errorScale(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<double>& errorScale){
88  switch (algo) {
89  case reco::TrackBase::ctf:
90  case reco::TrackBase::duplicateMerge:
91  case reco::TrackBase::initialStep:
92  case reco::TrackBase::highPtTripletStep:
93  case reco::TrackBase::lowPtQuadStep:
94  case reco::TrackBase::lowPtTripletStep:
95  case reco::TrackBase::pixelPairStep:
96  case reco::TrackBase::jetCoreRegionalStep:
97  case reco::TrackBase::muonSeededStepInOut:
98  case reco::TrackBase::muonSeededStepOutIn:
99  case reco::TrackBase::detachedQuadStep:
100  case reco::TrackBase::detachedTripletStep:
101  case reco::TrackBase::mixedTripletStep:
102  case reco::TrackBase::hltIter0:
103  case reco::TrackBase::hltIter1:
104  case reco::TrackBase::hltIter2:
105  case reco::TrackBase::hltIter3:
106  case reco::TrackBase::hltIter4:
107  case reco::TrackBase::hltIterX:
108  return 1.0;
109  case reco::TrackBase::pixelLessStep:
110  return errorScale[0];
111  case reco::TrackBase::tobTecStep:
112  return errorScale[1];
113  default:
114  return 1E9;
115  }
116  }
117 
118 
119 bool isGoodForEGM(const reco::TrackBase::TrackAlgorithm& algo){
120 
121 
122  switch (algo) {
123  case reco::TrackBase::ctf:
124  case reco::TrackBase::duplicateMerge:
125  case reco::TrackBase::initialStep:
126  case reco::TrackBase::highPtTripletStep:
127  case reco::TrackBase::lowPtQuadStep:
128  case reco::TrackBase::lowPtTripletStep:
129  case reco::TrackBase::pixelPairStep:
130  case reco::TrackBase::jetCoreRegionalStep:
131  case reco::TrackBase::detachedQuadStep:
132  case reco::TrackBase::detachedTripletStep:
133  case reco::TrackBase::mixedTripletStep:
134  case reco::TrackBase::muonSeededStepInOut:
135  case reco::TrackBase::muonSeededStepOutIn:
136  case reco::TrackBase::hltIter0:
137  case reco::TrackBase::hltIter1:
138  case reco::TrackBase::hltIter2:
139  case reco::TrackBase::hltIter3:
140  case reco::TrackBase::hltIter4:
141  case reco::TrackBase::hltIterX:
142  return true;
143  default:
144  return false;
145  }
146 
147 }
148 
149 bool isGoodForEGMPrimary(const reco::TrackBase::TrackAlgorithm& algo){
150  switch (algo) {
151  case reco::TrackBase::ctf:
152  case reco::TrackBase::duplicateMerge:
153  case reco::TrackBase::cosmics:
154  case reco::TrackBase::initialStep:
155  case reco::TrackBase::highPtTripletStep:
156  case reco::TrackBase::lowPtQuadStep:
157  case reco::TrackBase::lowPtTripletStep:
158  case reco::TrackBase::pixelPairStep:
159  case reco::TrackBase::detachedQuadStep:
160  case reco::TrackBase::detachedTripletStep:
161  case reco::TrackBase::mixedTripletStep:
162  return true;
163  default:
164  return false;
165  }
166 
167 }
168 
169 bool isFifthStep(const reco::TrackBase::TrackAlgorithm& algo){
170  switch (algo) {
171  case reco::TrackBase::ctf:
172  case reco::TrackBase::duplicateMerge:
173  case reco::TrackBase::initialStep:
174  case reco::TrackBase::highPtTripletStep:
175  case reco::TrackBase::lowPtQuadStep:
176  case reco::TrackBase::lowPtTripletStep:
177  case reco::TrackBase::pixelPairStep:
178  case reco::TrackBase::jetCoreRegionalStep:
179  case reco::TrackBase::detachedQuadStep:
180  case reco::TrackBase::detachedTripletStep:
181  case reco::TrackBase::mixedTripletStep:
182  case reco::TrackBase::pixelLessStep:
183  case reco::TrackBase::muonSeededStepInOut:
184  case reco::TrackBase::muonSeededStepOutIn:
185  case reco::TrackBase::hltIter0:
186  case reco::TrackBase::hltIter1:
187  case reco::TrackBase::hltIter2:
188  case reco::TrackBase::hltIter3:
189  case reco::TrackBase::hltIter4:
190  case reco::TrackBase::hltIterX:
191  return false;
192  case reco::TrackBase::tobTecStep:
193  return true;
194  default:
195  return true;
196  }
197 
198 }
199 
200 
201 bool highQuality(const reco::TrackBase::TrackAlgorithm& algo){
202  switch (algo) {
203  case reco::TrackBase::initialStep:
204  case reco::TrackBase::highPtTripletStep:
205  case reco::TrackBase::lowPtQuadStep:
206  case reco::TrackBase::lowPtTripletStep:
207  case reco::TrackBase::pixelPairStep:
208  case reco::TrackBase::detachedQuadStep:
209  case reco::TrackBase::detachedTripletStep:
210  case reco::TrackBase::duplicateMerge:
211  case reco::TrackBase::jetCoreRegionalStep:
212  return true;
213  default:
214  return false;
215  }
216 
217 }
218 
219 
220 bool nonIterative(const reco::TrackBase::TrackAlgorithm& algo){
221  switch (algo) {
222  case reco::TrackBase::undefAlgorithm:
223  case reco::TrackBase::ctf:
224  case reco::TrackBase::cosmics:
225  return true;
226  default:
227  return false;
228 
229  }
230 
231 }
232 
233 
234 bool step45(const reco::TrackBase::TrackAlgorithm& algo){
235  switch (algo) {
236  case reco::TrackBase::mixedTripletStep:
237  case reco::TrackBase::pixelLessStep:
238  case reco::TrackBase::tobTecStep:
239  return true;
240  default:
241  return false;
242  }
243 
244 }
245 
246 
247 bool step5(const reco::TrackBase::TrackAlgorithm& algo){
248  return (algo==reco::TrackBase::tobTecStep||algo==reco::TrackBase::pixelLessStep);
249 }
250 
251 #include "FWCore/MessageLogger/interface/MessageLogger.h"
252 bool goodPtResolution(const reco::TrackRef& trackref,
253  const std::vector<double>& DPtovPtCut,
254  const std::vector<unsigned>& NHitCut,
255  bool useIterTracking,
256  bool debug) {
257  //recheck that the track is high purity!
258  if (!trackref->quality(reco::TrackBase::highPurity))
259  return false;
260 
261  const double p = trackref->p();
262  const double pT = trackref->pt();
263  const double dpT = trackref->ptError();
264  const unsigned int nHit =
265  trackref->hitPattern().trackerLayersWithMeasurement();
266  const unsigned int nLostHit =
267  trackref->hitPattern().trackerLayersWithoutMeasurement(reco::HitPattern::TRACK_HITS);
268  const unsigned int lostHits = trackref->numberOfLostHits();
269  const double sigmaHad = std::sqrt(1.20*1.20/p+0.06*0.06) / (1.+lostHits);
270 
271  // Protection against 0 momentum tracks
272  if ( p < 0.05 ) return false;
273 
274  LogDebug("goodPtResolution") << " PFBlockAlgo: PFrecTrack->Track Pt= "
275  << pT << " DPt = " << dpT << std::endl;
276 
277 
278  double dptCut = dPtCut(trackref->algo(),DPtovPtCut,useIterTracking);
279  unsigned int nhitCut = nHitCut(trackref->algo(),NHitCut,useIterTracking);
280 
281  if ( ( dptCut > 0. &&
282  dpT/pT > dptCut*sigmaHad ) ||
283  nHit < nhitCut ) {
284  LogDebug("goodPtResolution") << " PFBlockAlgo: skip badly measured track"
285  << ", P = " << p
286  << ", Pt = " << pT
287  << " DPt = " << dpT
288  << ", N(hits) = " << nHit
289  << " (Lost : " << lostHits << "/" << nLostHit << ")"
290  << ", Algo = " << trackref->algo()
291  << std::endl;
292  LogDebug("goodPtResolution") << " cut is DPt/Pt < " << dptCut * sigmaHad << std::endl;
293  LogDebug("goodPtResolution") << " cut is NHit >= " << nhitCut << std::endl;
294  return false;
295  }
296 
297  return true;
298 }
299 
300 }
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
PFTrackAlgoTools::errorScale
double errorScale(const reco::TrackBase::TrackAlgorithm &, const std::vector< double > &)
Definition: PFTrackAlgoTools.cc:87
particleFlowClusterECALTimeSelected_cfi.cuts
cuts
Definition: particleFlowClusterECALTimeSelected_cfi.py:6
cutsCKF_cfi.lostHits
lostHits
Definition: cutsCKF_cfi.py:34
reco::TrackBase::TrackAlgorithm
TrackAlgorithm
track algorithm
Definition: TrackBase.h:99
LogDebug
#define LogDebug(id)
PFTrackAlgoTools::nHitCut
unsigned int nHitCut(const reco::TrackBase::TrackAlgorithm &, const std::vector< unsigned int > &, bool)
Definition: PFTrackAlgoTools.cc:46
PFTrackAlgoTools::isGoodForEGM
bool isGoodForEGM(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:119
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
PFTrackAlgoTools::nonIterative
bool nonIterative(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:220
PFTrackAlgoTools::step45
bool step45(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:234
PFTrackAlgoTools::highQuality
bool highQuality(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:201
debug
#define debug
Definition: HDRShower.cc:19
PFTrackAlgoTools::isFifthStep
bool isFifthStep(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:169
PFTrackAlgoTools::dPtCut
double dPtCut(const reco::TrackBase::TrackAlgorithm &, const std::vector< double > &, bool)
Definition: PFTrackAlgoTools.cc:5
PFTrackAlgoTools::isGoodForEGMPrimary
bool isGoodForEGMPrimary(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:149
PFTrackAlgoTools::step5
bool step5(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:247
patPFMETCorrections_cff.algo
algo
Definition: patPFMETCorrections_cff.py:123
PFTrackAlgoTools::goodPtResolution
bool goodPtResolution(const reco::TrackRef &, const std::vector< double > &DPtovPtCut, const std::vector< unsigned > &NHitCut, bool useIterTracking, bool debug=false)
Definition: PFTrackAlgoTools.cc:252
Generated for CMSSW Reference Manual by doxygen 1.8.11