CMS 3D CMS Logo
PFTrackAlgoTools.cc
Go to the documentation of this file.
1 #include "FWCore/MessageLogger/interface/MessageLogger.h"
2 #include "RecoParticleFlow/PFTracking/interface/PFTrackAlgoTools.h"
3 
4 namespace PFTrackAlgoTools {
5 
6  double dPtCut(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<double>& cuts,bool hltIterativeTracking = true){
7  switch (algo) {
8  case reco::TrackBase::ctf:
9  case reco::TrackBase::duplicateMerge:
10  case reco::TrackBase::initialStep:
11  case reco::TrackBase::highPtTripletStep:
12  case reco::TrackBase::lowPtQuadStep:
13  case reco::TrackBase::lowPtTripletStep:
14  case reco::TrackBase::pixelPairStep:
15  case reco::TrackBase::jetCoreRegionalStep:
16  return cuts[0];
17  case reco::TrackBase::detachedQuadStep:
18  case reco::TrackBase::detachedTripletStep:
19  return cuts[1];
20  case reco::TrackBase::mixedTripletStep:
21  return cuts[2];
22  case reco::TrackBase::pixelLessStep:
23  return cuts[3];
24  case reco::TrackBase::tobTecStep:
25  return cuts[4];
26  case reco::TrackBase::muonSeededStepInOut:
27  case reco::TrackBase::muonSeededStepOutIn:
28  return cuts.at(5);
29  case reco::TrackBase::hltIter0:
30  case reco::TrackBase::hltIter1:
31  case reco::TrackBase::hltIter2:
32  return cuts[0];
33  case reco::TrackBase::hltIter3:
34  return hltIterativeTracking ? cuts[1] : cuts[0];
35  case reco::TrackBase::hltIter4:
36  return hltIterativeTracking ? cuts[2] : cuts[0];
37  case reco::TrackBase::hltIterX:
38  return cuts[0];
39  default:
40  return hltIterativeTracking ? cuts.at(5):cuts[0];
41 
42  }
43  }
44 
45 
46 
47  unsigned int nHitCut(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<unsigned int>& cuts,bool hltIterativeTracking = true){
48  switch (algo) {
49  case reco::TrackBase::ctf:
50  case reco::TrackBase::duplicateMerge:
51  case reco::TrackBase::initialStep:
52  case reco::TrackBase::highPtTripletStep:
53  case reco::TrackBase::lowPtQuadStep:
54  case reco::TrackBase::lowPtTripletStep:
55  case reco::TrackBase::pixelPairStep:
56  case reco::TrackBase::jetCoreRegionalStep:
57  return cuts[0];
58  case reco::TrackBase::detachedQuadStep:
59  case reco::TrackBase::detachedTripletStep:
60  return cuts[1];
61  case reco::TrackBase::mixedTripletStep:
62  return cuts[2];
63  case reco::TrackBase::pixelLessStep:
64  return cuts[3];
65  case reco::TrackBase::tobTecStep:
66  return cuts[4];
67  case reco::TrackBase::muonSeededStepInOut:
68  case reco::TrackBase::muonSeededStepOutIn:
69  return cuts.at(5);
70  case reco::TrackBase::hltIter0:
71  case reco::TrackBase::hltIter1:
72  case reco::TrackBase::hltIter2:
73  return cuts[0];
74  case reco::TrackBase::hltIter3:
75  return hltIterativeTracking ? cuts[1] : cuts[0];
76  case reco::TrackBase::hltIter4:
77  return hltIterativeTracking ? cuts[2] : cuts[0];
78  case reco::TrackBase::hltIterX:
79  return cuts[0];
80  default:
81  return hltIterativeTracking ? cuts.at(5):cuts[0];
82 
83  }
84  }
85 
86 
87 
88  double errorScale(const reco::TrackBase::TrackAlgorithm& algo,const std::vector<double>& errorScale){
89  switch (algo) {
90  case reco::TrackBase::ctf:
91  case reco::TrackBase::duplicateMerge:
92  case reco::TrackBase::initialStep:
93  case reco::TrackBase::highPtTripletStep:
94  case reco::TrackBase::lowPtQuadStep:
95  case reco::TrackBase::lowPtTripletStep:
96  case reco::TrackBase::pixelPairStep:
97  case reco::TrackBase::jetCoreRegionalStep:
98  case reco::TrackBase::muonSeededStepInOut:
99  case reco::TrackBase::muonSeededStepOutIn:
100  case reco::TrackBase::detachedQuadStep:
101  case reco::TrackBase::detachedTripletStep:
102  case reco::TrackBase::mixedTripletStep:
103  case reco::TrackBase::hltIter0:
104  case reco::TrackBase::hltIter1:
105  case reco::TrackBase::hltIter2:
106  case reco::TrackBase::hltIter3:
107  case reco::TrackBase::hltIter4:
108  case reco::TrackBase::hltIterX:
109  return 1.0;
110  case reco::TrackBase::pixelLessStep:
111  return errorScale[0];
112  case reco::TrackBase::tobTecStep:
113  return errorScale[1];
114  default:
115  return 1E9;
116  }
117  }
118 
119 
120 bool isGoodForEGM(const reco::TrackBase::TrackAlgorithm& algo){
121 
122 
123  switch (algo) {
124  case reco::TrackBase::ctf:
125  case reco::TrackBase::duplicateMerge:
126  case reco::TrackBase::initialStep:
127  case reco::TrackBase::highPtTripletStep:
128  case reco::TrackBase::lowPtQuadStep:
129  case reco::TrackBase::lowPtTripletStep:
130  case reco::TrackBase::pixelPairStep:
131  case reco::TrackBase::jetCoreRegionalStep:
132  case reco::TrackBase::detachedQuadStep:
133  case reco::TrackBase::detachedTripletStep:
134  case reco::TrackBase::mixedTripletStep:
135  case reco::TrackBase::muonSeededStepInOut:
136  case reco::TrackBase::muonSeededStepOutIn:
137  case reco::TrackBase::hltIter0:
138  case reco::TrackBase::hltIter1:
139  case reco::TrackBase::hltIter2:
140  case reco::TrackBase::hltIter3:
141  case reco::TrackBase::hltIter4:
142  case reco::TrackBase::hltIterX:
143  return true;
144  default:
145  return false;
146  }
147 
148 }
149 
150 bool isGoodForEGMPrimary(const reco::TrackBase::TrackAlgorithm& algo){
151  switch (algo) {
152  case reco::TrackBase::ctf:
153  case reco::TrackBase::duplicateMerge:
154  case reco::TrackBase::cosmics:
155  case reco::TrackBase::initialStep:
156  case reco::TrackBase::highPtTripletStep:
157  case reco::TrackBase::lowPtQuadStep:
158  case reco::TrackBase::lowPtTripletStep:
159  case reco::TrackBase::pixelPairStep:
160  case reco::TrackBase::detachedQuadStep:
161  case reco::TrackBase::detachedTripletStep:
162  case reco::TrackBase::mixedTripletStep:
163  return true;
164  default:
165  return false;
166  }
167 
168 }
169 
170 bool isFifthStep(const reco::TrackBase::TrackAlgorithm& algo){
171  switch (algo) {
172  case reco::TrackBase::ctf:
173  case reco::TrackBase::duplicateMerge:
174  case reco::TrackBase::initialStep:
175  case reco::TrackBase::highPtTripletStep:
176  case reco::TrackBase::lowPtQuadStep:
177  case reco::TrackBase::lowPtTripletStep:
178  case reco::TrackBase::pixelPairStep:
179  case reco::TrackBase::jetCoreRegionalStep:
180  case reco::TrackBase::detachedQuadStep:
181  case reco::TrackBase::detachedTripletStep:
182  case reco::TrackBase::mixedTripletStep:
183  case reco::TrackBase::pixelLessStep:
184  case reco::TrackBase::muonSeededStepInOut:
185  case reco::TrackBase::muonSeededStepOutIn:
186  case reco::TrackBase::hltIter0:
187  case reco::TrackBase::hltIter1:
188  case reco::TrackBase::hltIter2:
189  case reco::TrackBase::hltIter3:
190  case reco::TrackBase::hltIter4:
191  case reco::TrackBase::hltIterX:
192  return false;
193  case reco::TrackBase::tobTecStep:
194  return true;
195  default:
196  return true;
197  }
198 
199 }
200 
201 
202 bool highQuality(const reco::TrackBase::TrackAlgorithm& algo){
203  switch (algo) {
204  case reco::TrackBase::initialStep:
205  case reco::TrackBase::highPtTripletStep:
206  case reco::TrackBase::lowPtQuadStep:
207  case reco::TrackBase::lowPtTripletStep:
208  case reco::TrackBase::pixelPairStep:
209  case reco::TrackBase::detachedQuadStep:
210  case reco::TrackBase::detachedTripletStep:
211  case reco::TrackBase::duplicateMerge:
212  case reco::TrackBase::jetCoreRegionalStep:
213  return true;
214  default:
215  return false;
216  }
217 
218 }
219 
220 
221 bool nonIterative(const reco::TrackBase::TrackAlgorithm& algo){
222  switch (algo) {
223  case reco::TrackBase::undefAlgorithm:
224  case reco::TrackBase::ctf:
225  case reco::TrackBase::cosmics:
226  return true;
227  default:
228  return false;
229 
230  }
231 
232 }
233 
234 
235 bool step45(const reco::TrackBase::TrackAlgorithm& algo){
236  switch (algo) {
237  case reco::TrackBase::mixedTripletStep:
238  case reco::TrackBase::pixelLessStep:
239  case reco::TrackBase::tobTecStep:
240  return true;
241  default:
242  return false;
243  }
244 
245 }
246 
247 
248 bool step5(const reco::TrackBase::TrackAlgorithm& algo){
249  return (algo==reco::TrackBase::tobTecStep||algo==reco::TrackBase::pixelLessStep);
250 }
251 
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 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
egammaForCoreTracking_cff.cuts
cuts
Definition: egammaForCoreTracking_cff.py:40
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
PFTrackAlgoTools::errorScale
double errorScale(const reco::TrackBase::TrackAlgorithm &, const std::vector< double > &)
Definition: PFTrackAlgoTools.cc:88
reco::TrackBase::TrackAlgorithm
TrackAlgorithm
track algorithm
Definition: TrackBase.h:99
PFTrackAlgoTools::nHitCut
unsigned int nHitCut(const reco::TrackBase::TrackAlgorithm &, const std::vector< unsigned int > &, bool)
Definition: PFTrackAlgoTools.cc:47
PFTrackAlgoTools::isGoodForEGM
bool isGoodForEGM(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:120
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
PFTrackAlgoTools::nonIterative
bool nonIterative(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:221
PFTrackAlgoTools::step45
bool step45(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:235
PFTrackAlgoTools::highQuality
bool highQuality(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:202
electrons_cff.lostHits
lostHits
Definition: electrons_cff.py:384
debug
#define debug
Definition: HDRShower.cc:19
PFTrackAlgoTools::isFifthStep
bool isFifthStep(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:170
PFTrackAlgoTools::dPtCut
double dPtCut(const reco::TrackBase::TrackAlgorithm &, const std::vector< double > &, bool)
Definition: PFTrackAlgoTools.cc:6
PFTrackAlgoTools::isGoodForEGMPrimary
bool isGoodForEGMPrimary(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:150
PFTrackAlgoTools::step5
bool step5(const reco::TrackBase::TrackAlgorithm &)
Definition: PFTrackAlgoTools.cc:248
patPFMETCorrections_cff.algo
algo
Definition: patPFMETCorrections_cff.py:119
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