CMS 3D CMS Logo
HitPattern.h
Go to the documentation of this file.
1 // -*- C++ -*-
2 
3 #ifndef TrackReco_HitPattern_h
4 #define TrackReco_HitPattern_h
5 
6 //
7 // File: DataFormats/TrackReco/interface/HitPattern.h
8 //
9 // Marcel Vos, INFN Pisa
10 // v1.10 2007/05/08 bellan
11 // Zongru Wan, Kansas State University
12 // Jean-Roch Vlimant
13 // Kevin Burkett
14 // Boris Mangano
15 //
16 // Hit pattern is the summary information of the hits associated to track in
17 // AOD. When RecHits are no longer available, the compact hit pattern should
18 // allow basic track selection based on the hits in various subdetectors.
19 // The hits of a track are saved in unit16_t hitPattern[MaxHits].
20 //
21 // uint16_t
22 // +------------+---------------+---------------------------+-----------------+----------------+
23 // | tk/mu/mtd | sub-structure | sub-sub-structure | stereo | hit type |
24 // +------------+---------------+---------------------------+-----------------+----------------+
25 // | 11-10 | 9 8 7 | 6 5 4 3 | 2 | 1 0 | bit
26 // +------------+---------------+---------------------------+-----------------+----------------|
27 // | tk = 1 | PXB = 1 | layer = 1-3 | | hit type = 0-3 |
28 // | tk = 1 | PXF = 2 | disk = 1-2 | | hit type = 0-3 |
29 // | tk = 1 | TIB = 3 | layer = 1-4 | 0=rphi,1=stereo | hit type = 0-3 |
30 // | tk = 1 | TID = 4 | wheel = 1-3 | 0=rphi,1=stereo | hit type = 0-3 |
31 // | tk = 1 | TOB = 5 | layer = 1-6 | 0=rphi,1=stereo | hit type = 0-3 |
32 // | tk = 1 | TEC = 6 | wheel = 1-9 | 0=rphi,1=stereo | hit type = 0-3 |
33 // | mu = 0 | DT = 1 | 4*(stat-1)+superlayer | | hit type = 0-3 |
34 // | mu = 0 | CSC = 2 | 4*(stat-1)+(ring-1) | | hit type = 0-3 |
35 // | mu = 0 | RPC = 3 | 4*(stat-1)+2*layer+region | | hit type = 0-3 |
36 // | mu = 0 | GEM = 4 | 2*(stat-1)+2*(layer-1) | | hit type = 0-3 |
37 // | mu = 0 | ME0 = 5 | roll | | hit type = 0-3 |
38 // | mtd = 2 | BTL = 1 | moduleType = 1-3 | | hit type = 0-3 |
39 // | mtd = 2 | ETL = 2 | ring = 1-12 | | hit type = 0-3 |
40 // +------------+---------------+---------------------------+-----------------+----------------+
41 //
42 // hit type, see DataFormats/TrackingRecHit/interface/TrackingRecHit.h
43 // VALID = valid hit = 0
44 // MISSING = detector is good, but no rec hit found = 1
45 // INACTIVE = detector is off, so there was no hope = 2
46 // BAD = there were many bad strips within the ellipse = 3
47 //
48 // It had been shown by Zongru using a 100 GeV muon sample with 5000 events
49 // uniform in eta and phi, the average (maximum) number of tracker hits is
50 // 13 (17) and the average (maximum) number of muon detector hits is about
51 // 26 (50). If the number of hits of a track is larger than 80 then the extra
52 // hits are ignored by hit pattern. The static hit pattern array might be
53 // improved to a dynamic one in the future.
54 //
55 // Because of tracking with/without overlaps and with/without hit-splitting,
56 // the final number of hits per track is pretty "variable". Compared with the
57 // number of valid hits, the number of crossed layers with measurement should
58 // be more robust to discriminate between good and fake track.
59 //
60 // Since 4-bit for sub-sub-structure is not enough to specify a muon layer,
61 // the layer case counting methods are implemented for tracker only. This is
62 // different from the hit counting methods which are implemented for both
63 // tracker and muon detector.
64 //
65 // Given a tracker layer, specified by sub-structure and layer, the method
66 // getTrackerLayerCase(substr, layer) groups all of the hits in the hit pattern
67 // array for the layer together and returns one of the four cases
68 //
69 // crossed
70 // layer case 0: VALID + (MISSING, OFF, BAD) ==> with measurement
71 // layer case 1: MISSING + (OFF, BAD) ==> without measurement
72 // layer case 2: OFF, BAD ==> totally off or bad, cannot say much
73 // not crossed
74 // layer case NULL_RETURN: track outside acceptance or in gap ==> null
75 //
76 // Given a tracker layer, specified by sub-structure and layer, the method
77 // getTrackerMonoStereo(substr, layer) groups all of the valid hits in the hit
78 // pattern array for the layer together and returns
79 //
80 // 0: neither a valid mono nor a valid stereo hit
81 // MONO: valid mono hit
82 // STEREO: valid stereo hit
83 // MONO | STEREO: both
84 //
85 //
86 // Given a track, here is an example usage of hit pattern
87 //
88 // // hit pattern of the track
89 // const reco::HitPattern &p = track->hitPattern();
90 //
91 // // loop over the hits of the track.
92 // for (int i = 0; i < p.numberOfAllHits(HitPattern::TRACK_HITS); i++) {
93 // uint32_t hit = p.getHitPattern(HitPattern::TRACK_HITS, i);
94 //
95 // // if the hit is valid and in pixel barrel, print out the layer
96 // if (p.validHitFilter(hit) && p.pixelBarrelHitFilter(hit)){
97 // cout << "valid hit found in pixel barrel layer "
98 // << p.getLayer(hit)
99 // << endl;
100 // }
101 //
102 // // expert level: printout the hit in 11-bit binary format
103 // cout << "hit in 11-bit binary format = ";
104 // for (int j = 10; j >= 0; j--){
105 // int bit = (hit >> j) & 0x1;
106 // cout << bit;
107 // }
108 // cout << endl;
109 // }
110 //
111 // //count the number of valid pixel barrel *** hits ***
112 // cout << "number of of valid pixel barrel hits is "
113 // << p.numberOfValidPixelBarrelHits()
114 // << endl;
115 //
116 // //count the number of pixel barrel *** layers *** with measurement
117 // cout << "number of of pixel barrel layers with measurement is "
118 // << p.pixelBarrelLayersWithMeasurement()
119 // << endl;
120 //
121 
122 #include "DataFormats/DetId/interface/DetId.h"
123 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
124 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
125 #include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
126 #include "DataFormats/ForwardDetId/interface/MTDDetId.h"
127 #include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
128 #include "DataFormats/TrackingRecHit/interface/TrackingRecHitFwd.h"
129 #include "FWCore/Utilities/interface/Likely.h"
130 
131 #include <utility>
132 #include <algorithm>
133 #include <iostream>
134 #include <ostream>
135 #include <memory>
136 
137 class TrackerTopology;
138 
139 namespace test {
140  namespace TestHitPattern {
141  int test();
142  }
143 } // namespace test
144 
145 namespace reco {
146 
147  class HitPattern {
148  public:
149  enum { MONO = 1, STEREO = 2 };
150 
151  enum HIT_DETECTOR_TYPE { MUON_HIT = 0, TRACKER_HIT = 1, MTD_HIT = 2 };
152 
153  enum HIT_TYPE { VALID = 0, MISSING = 1, INACTIVE = 2, BAD = 3 };
154 
155  enum HitCategory { TRACK_HITS = 0, MISSING_INNER_HITS = 1, MISSING_OUTER_HITS = 2 };
156  const static unsigned short ARRAY_LENGTH = 57;
157  const static unsigned short HIT_LENGTH = 12;
158  const static unsigned short MaxHits = (8 * sizeof(uint16_t) * ARRAY_LENGTH) / HIT_LENGTH;
159 
160  static const uint32_t NULL_RETURN = 999999;
161  static const uint16_t EMPTY_PATTERN = 0x0;
162 
163  static bool trackerHitFilter(uint16_t pattern);
164  static bool muonHitFilter(uint16_t pattern);
165  static bool timingHitFilter(uint16_t pattern);
166 
167  static bool validHitFilter(uint16_t pattern);
168  static bool missingHitFilter(uint16_t pattern);
169  static bool inactiveHitFilter(uint16_t pattern);
170  static bool badHitFilter(uint16_t pattern);
171 
172  static bool pixelHitFilter(uint16_t pattern);
173  static bool pixelBarrelHitFilter(uint16_t pattern);
174  static bool pixelEndcapHitFilter(uint16_t pattern);
175  static bool stripHitFilter(uint16_t pattern);
176  static bool stripTIBHitFilter(uint16_t pattern);
177  static bool stripTIDHitFilter(uint16_t pattern);
178  static bool stripTOBHitFilter(uint16_t pattern);
179  static bool stripTECHitFilter(uint16_t pattern);
180 
181  static bool muonDTHitFilter(uint16_t pattern);
182  static bool muonCSCHitFilter(uint16_t pattern);
183  static bool muonRPCHitFilter(uint16_t pattern);
184  static bool muonGEMHitFilter(uint16_t pattern);
185  static bool muonME0HitFilter(uint16_t pattern);
186 
187  static bool timingBTLHitFilter(uint16_t pattern);
188  static bool timingETLHitFilter(uint16_t pattern);
189 
190  static uint32_t getHitType(uint16_t pattern);
191 
192  // mono (0) or stereo (1)
193  static uint32_t getSide(uint16_t pattern);
194  static uint32_t getLayer(uint16_t pattern);
195  static uint32_t getSubSubStructure(uint16_t pattern);
196  static uint32_t getSubStructure(uint16_t pattern);
197  static uint32_t getSubDetector(uint16_t pattern);
198 
200  static uint16_t getMuonStation(uint16_t pattern);
201 
203  static uint16_t getDTSuperLayer(uint16_t pattern); // only for DT patterns
204 
206  static uint16_t getCSCRing(uint16_t pattern);
207 
209  static uint16_t getRPCLayer(uint16_t pattern);
210 
212  static uint16_t getRPCregion(uint16_t pattern);
213 
215  static uint16_t getGEMStation(uint16_t pattern);
216 
218  static uint16_t getGEMLayer(uint16_t pattern);
219 
221  static uint16_t getBTLModType(uint16_t pattern);
222 
224  static uint16_t getETLRing(uint16_t pattern);
225 
226  HitPattern();
227 
228  ~HitPattern();
229 
230  HitPattern(const HitPattern &other);
231 
232  HitPattern &operator=(const HitPattern &other);
233 
234  template <typename I>
235  bool appendHits(const I &begin, const I &end, const TrackerTopology &ttopo);
236  bool appendHit(const TrackingRecHit &hit, const TrackerTopology &ttopo);
237  bool appendHit(const TrackingRecHitRef &ref, const TrackerTopology &ttopo);
238  bool appendHit(const DetId &id, TrackingRecHit::Type hitType, const TrackerTopology &ttopo);
239  bool appendHit(const uint16_t pattern, TrackingRecHit::Type hitType);
240 
247  bool appendTrackerHit(uint16_t subdet, uint16_t layer, uint16_t stereo, TrackingRecHit::Type hitType);
248 
255  bool appendMuonHit(const DetId &id, TrackingRecHit::Type hitType);
256 
257  // get the pattern of the position-th hit
258  uint16_t getHitPattern(HitCategory category, int position) const;
259 
260  void clear();
261 
262  // print the pattern of the position-th hit
263  void printHitPattern(HitCategory category, int position, std::ostream &stream) const;
264  void print(HitCategory category, std::ostream &stream = std::cout) const;
265 
266  // has valid hit in PXB/PXF layer x
267  bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const;
268 
269  int numberOfAllHits(HitCategory category) const; // not-null
270  int numberOfValidHits() const; // not-null, valid
271 
272  int numberOfAllTrackerHits(HitCategory category) const; // not-null, tracker
273  int numberOfValidTrackerHits() const; // not-null, valid, tracker
274  int numberOfValidPixelHits() const; // not-null, valid, pixel
275  int numberOfValidPixelBarrelHits() const; // not-null, valid, pixel PXB
276  int numberOfValidPixelEndcapHits() const; // not-null, valid, pixel PXF
277  int numberOfValidStripHits() const; // not-null, valid, strip
278  int numberOfValidStripTIBHits() const; // not-null, valid, strip TIB
279  int numberOfValidStripTIDHits() const; // not-null, valid, strip TID
280  int numberOfValidStripTOBHits() const; // not-null, valid, strip TOB
281  int numberOfValidStripTECHits() const; // not-null, valid, strip TEC
282 
283  int numberOfLostHits(HitCategory category) const; // not-null, not valid
284  int numberOfLostTrackerHits(HitCategory category) const; // not-null, not valid, tracker
285  int numberOfLostPixelHits(HitCategory category) const; // not-null, not valid, pixel
286  int numberOfLostPixelBarrelHits(HitCategory category) const; // not-null, not valid, pixel PXB
287  int numberOfLostPixelEndcapHits(HitCategory category) const; // not-null, not valid, pixel PXF
288  int numberOfLostStripHits(HitCategory category) const; // not-null, not valid, strip
289  int numberOfLostStripTIBHits(HitCategory category) const; // not-null, not valid, strip TIB
290  int numberOfLostStripTIDHits(HitCategory category) const; // not-null, not valid, strip TID
291  int numberOfLostStripTOBHits(HitCategory category) const; // not-null, not valid, strip TOB
292  int numberOfLostStripTECHits(HitCategory category) const; // not-null, not valid, strip TEC
293 
294  int numberOfTimingHits() const; // not-null timing
295  int numberOfValidTimingHits() const; // not-null, valid, timing
296  int numberOfValidTimingBTLHits() const; // not-null, valid, timing BTL
297  int numberOfValidTimingETLHits() const; // not-null, valid, timing ETL
298 
299  int numberOfLostTimingHits() const; // not-null, not valid, timing
300  int numberOfLostTimingBTLHits() const; // not-null, not valid, timing BTL
301  int numberOfLostTimingETLHits() const; // not-null, not valid, timing ETL
302 
303  int numberOfMuonHits() const; // not-null, muon
304  int numberOfValidMuonHits() const; // not-null, valid, muon
305  int numberOfValidMuonDTHits() const; // not-null, valid, muon DT
306  int numberOfValidMuonCSCHits() const; // not-null, valid, muon CSC
307  int numberOfValidMuonRPCHits() const; // not-null, valid, muon RPC
308  int numberOfValidMuonGEMHits() const; // not-null, valid, muon GEM
309  int numberOfValidMuonME0Hits() const; // not-null, valid, muon ME0
310 
311  int numberOfLostMuonHits() const; // not-null, not valid, muon
312  int numberOfLostMuonDTHits() const; // not-null, not valid, muon DT
313  int numberOfLostMuonCSCHits() const; // not-null, not valid, muon CSC
314  int numberOfLostMuonRPCHits() const; // not-null, not valid, muon RPC
315  int numberOfLostMuonGEMHits() const; // not-null, not valid, muon GEM
316  int numberOfLostMuonME0Hits() const; // not-null, not valid, muon ME0
317 
318  int numberOfBadHits() const; // not-null, bad (only used in Muon Ch.)
319  int numberOfBadMuonHits() const; // not-null, bad, muon
320  int numberOfBadMuonDTHits() const; // not-null, bad, muon DT
321  int numberOfBadMuonCSCHits() const; // not-null, bad, muon CSC
322  int numberOfBadMuonRPCHits() const; // not-null, bad, muon RPC
323  int numberOfBadMuonGEMHits() const; // not-null, bad, muon GEM
324  int numberOfBadMuonME0Hits() const; // not-null, bad, muon ME0
325 
326  int numberOfInactiveHits() const; // not-null, inactive
327  int numberOfInactiveTrackerHits() const; // not-null, inactive, tracker
328  int numberOfInactiveTimingHits() const; // not-null, inactive, timing
329 
330  // count strip layers that have non-null, valid mono and stereo hits
331  int numberOfValidStripLayersWithMonoAndStereo(uint16_t stripdet, uint16_t layer) const;
332  int numberOfValidStripLayersWithMonoAndStereo() const;
333  int numberOfValidTOBLayersWithMonoAndStereo(uint32_t layer = 0) const;
334  int numberOfValidTIBLayersWithMonoAndStereo(uint32_t layer = 0) const;
335  int numberOfValidTIDLayersWithMonoAndStereo(uint32_t layer = 0) const;
336  int numberOfValidTECLayersWithMonoAndStereo(uint32_t layer = 0) const;
337 
338  uint32_t getTrackerLayerCase(HitCategory category, uint16_t substr, uint16_t layer) const;
339  uint16_t getTrackerMonoStereo(HitCategory category, uint16_t substr, uint16_t layer) const;
340 
341  int trackerLayersWithMeasurementOld() const; // case 0: tracker
342  int trackerLayersWithMeasurement() const; // case 0: tracker
343  int pixelLayersWithMeasurementOld() const; // case 0: pixel
344  int pixelLayersWithMeasurement() const; // case 0: pixel
345  int stripLayersWithMeasurement() const; // case 0: strip
346  int pixelBarrelLayersWithMeasurement() const; // case 0: pixel PXB
347  int pixelEndcapLayersWithMeasurement() const; // case 0: pixel PXF
348  int stripTIBLayersWithMeasurement() const; // case 0: strip TIB
349  int stripTIDLayersWithMeasurement() const; // case 0: strip TID
350  int stripTOBLayersWithMeasurement() const; // case 0: strip TOB
351  int stripTECLayersWithMeasurement() const; // case 0: strip TEC
352 
353  int trackerLayersWithoutMeasurement(HitCategory category) const; // case 1: tracker
354  int trackerLayersWithoutMeasurementOld(HitCategory category) const; // case 1: tracker
355  int pixelLayersWithoutMeasurement(HitCategory category) const; // case 1: pixel
356  int stripLayersWithoutMeasurement(HitCategory category) const; // case 1: strip
357  int pixelBarrelLayersWithoutMeasurement(HitCategory category) const; // case 1: pixel PXB
358  int pixelEndcapLayersWithoutMeasurement(HitCategory category) const; // case 1: pixel PXF
359  int stripTIBLayersWithoutMeasurement(HitCategory category) const; // case 1: strip TIB
360  int stripTIDLayersWithoutMeasurement(HitCategory category) const; // case 1: strip TID
361  int stripTOBLayersWithoutMeasurement(HitCategory category) const; // case 1: strip TOB
362  int stripTECLayersWithoutMeasurement(HitCategory category) const; // case 1: strip TEC
363 
364  int trackerLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: tracker
365  int pixelLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: pixel
366  int stripLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: strip
367  int pixelBarrelLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: pixel PXB
368  int pixelEndcapLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: pixel PXF
369  int stripTIBLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: strip TIB
370  int stripTIDLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: strip TID
371  int stripTOBLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: strip TOB
372  int stripTECLayersTotallyOffOrBad(HitCategory category = TRACK_HITS) const; // case 2: strip TEC
373 
374  int trackerLayersNull() const; // case NULL_RETURN: tracker
375  int pixelLayersNull() const; // case NULL_RETURN: pixel
376  int stripLayersNull() const; // case NULL_RETURN: strip
377  int pixelBarrelLayersNull() const; // case NULL_RETURN: pixel PXB
378  int pixelEndcapLayersNull() const; // case NULL_RETURN: pixel PXF
379  int stripTIBLayersNull() const; // case NULL_RETURN: strip TIB
380  int stripTIDLayersNull() const; // case NULL_RETURN: strip TID
381  int stripTOBLayersNull() const; // case NULL_RETURN: strip TOB
382  int stripTECLayersNull() const; // case NULL_RETURN: strip TEC
383 
385  int muonStations(int subdet, int hitType) const;
386 
387  int muonStationsWithValidHits() const;
388  int muonStationsWithBadHits() const;
389  int muonStationsWithAnyHits() const;
390 
391  int dtStationsWithValidHits() const;
392  int dtStationsWithBadHits() const;
393  int dtStationsWithAnyHits() const;
394 
395  int cscStationsWithValidHits() const;
396  int cscStationsWithBadHits() const;
397  int cscStationsWithAnyHits() const;
398 
399  int rpcStationsWithValidHits() const;
400  int rpcStationsWithBadHits() const;
401  int rpcStationsWithAnyHits() const;
402 
403  int gemStationsWithValidHits() const;
404  int gemStationsWithBadHits() const;
405  int gemStationsWithAnyHits() const;
406 
407  int me0StationsWithValidHits() const;
408  int me0StationsWithBadHits() const;
409  int me0StationsWithAnyHits() const;
410 
412  int innermostMuonStationWithHits(int hitType) const;
413  int innermostMuonStationWithValidHits() const;
414  int innermostMuonStationWithBadHits() const;
415  int innermostMuonStationWithAnyHits() const;
416 
418  int outermostMuonStationWithHits(int hitType) const;
419  int outermostMuonStationWithValidHits() const;
420  int outermostMuonStationWithBadHits() const;
421  int outermostMuonStationWithAnyHits() const;
422 
423  int numberOfDTStationsWithRPhiView() const;
424  int numberOfDTStationsWithRZView() const;
425  int numberOfDTStationsWithBothViews() const;
426 
427  //only used by ROOT IO rule to read v12 HitPatterns
428  static bool fillNewHitPatternWithOldHitPattern_v12(const uint16_t oldHitPattern[],
429  uint8_t hitCount,
430  uint8_t beginTrackHits,
431  uint8_t endTrackHits,
432  uint8_t beginInner,
433  uint8_t endInner,
434  uint8_t beginOuter,
435  uint8_t endOuter,
436  reco::HitPattern *newObj);
437 
438  private:
439  // 3 bits for hit type
440  const static unsigned short HitTypeOffset = 0;
441  const static unsigned short HitTypeMask = 0x3;
442 
443  // 1 bit to identify the side in double-sided detectors
444  const static unsigned short SideOffset = 2;
445  const static unsigned short SideMask = 0x1;
446 
447  // 4 bits to identify the layer/disk/wheel within the substructure
448  const static unsigned short LayerOffset = 3;
449  const static unsigned short LayerMask = 0xF;
450 
451  // 3 bits to identify the tracker/muon detector substructure
452  const static unsigned short SubstrOffset = 7;
453  const static unsigned short SubstrMask = 0x7;
454 
455  // 2 bits to distinguish tracker, muon, mtd subsystems
456  const static unsigned short SubDetectorOffset = 10;
457  const static unsigned short SubDetectorMask = 0x3;
458 
459  const static unsigned short minTrackerWord = 1 << SubDetectorOffset;
460  const static unsigned short maxTrackerWord = (2 << SubDetectorOffset) - 1;
461  const static unsigned short minPixelWord = minTrackerWord | (1 << SubstrOffset);
462  const static unsigned short minStripWord = minTrackerWord | (3 << SubstrOffset);
463 
464  // detector side for tracker modules (mono/stereo)
465  static uint16_t isStereo(DetId i, const TrackerTopology &ttopo);
466  static bool stripSubdetectorHitFilter(uint16_t pattern, StripSubdetector::SubDetector substructure);
467 
468  static uint16_t encode(const TrackingRecHit &hit, const TrackerTopology &ttopo);
469  static uint16_t encode(const DetId &id, TrackingRecHit::Type hitType, const TrackerTopology &ttopo);
470  static uint16_t encode(uint16_t det, uint16_t subdet, uint16_t layer, uint16_t side, TrackingRecHit::Type hitType);
471 
472  // generic count methods
473  typedef bool filterType(uint16_t);
474 
475  template <typename F>
476  void call(HitCategory category, filterType typeFilter, F f) const;
477 
478  int countHits(HitCategory category, filterType filter) const;
479  int countTypedHits(HitCategory category, filterType typeFilter, filterType filter) const;
480 
481  bool insertTrackHit(const uint16_t pattern);
482  bool insertExpectedInnerHit(const uint16_t pattern);
483  bool insertExpectedOuterHit(const uint16_t pattern);
484  void insertHit(const uint16_t pattern);
485 
486  uint16_t getHitPatternByAbsoluteIndex(int position) const;
487 
488  std::pair<uint8_t, uint8_t> getCategoryIndexRange(HitCategory category) const;
489 
490  uint16_t hitPattern[ARRAY_LENGTH];
491  uint8_t hitCount;
492 
493  uint8_t beginTrackHits;
494  uint8_t endTrackHits;
495  uint8_t beginInner;
496  uint8_t endInner;
497  uint8_t beginOuter;
498  uint8_t endOuter;
499 
500  friend int ::test::TestHitPattern::test();
501 
502  template <int N>
503  friend struct PatternSet;
504  };
505 
506  inline std::pair<uint8_t, uint8_t> HitPattern::getCategoryIndexRange(HitCategory category) const {
507  switch (category) {
508  case TRACK_HITS:
509  return std::pair<uint8_t, uint8_t>(beginTrackHits, endTrackHits);
510  break;
511  case MISSING_INNER_HITS:
512  return std::pair<uint8_t, uint8_t>(beginInner, endInner);
513  break;
514  case MISSING_OUTER_HITS:
515  return std::pair<uint8_t, uint8_t>(beginOuter, endOuter);
516  break;
517  }
518  return std::pair<uint8_t, uint8_t>(-1, -1);
519  }
520 
521  template <typename I>
522  bool HitPattern::appendHits(const I &begin, const I &end, const TrackerTopology &ttopo) {
523  for (I hit = begin; hit != end; hit++) {
524  if UNLIKELY ((!appendHit(*hit, ttopo))) {
525  return false;
526  }
527  }
528  return true;
529  }
530 
531  inline uint16_t HitPattern::getHitPattern(HitCategory category, int position) const {
532  std::pair<uint8_t, uint8_t> range = getCategoryIndexRange(category);
533  if UNLIKELY ((position < 0 || (position + range.first) >= range.second)) {
534  return HitPattern::EMPTY_PATTERN;
535  }
536 
537  return getHitPatternByAbsoluteIndex(range.first + position);
538  }
539 
540  inline int HitPattern::countHits(HitCategory category, filterType filter) const {
541  int count = 0;
542  std::pair<uint8_t, uint8_t> range = getCategoryIndexRange(category);
543  for (int i = range.first; i < range.second; ++i) {
544  if (filter(getHitPatternByAbsoluteIndex(i))) {
545  ++count;
546  }
547  }
548  return count;
549  }
550 
551  template <typename F>
552  void HitPattern::call(HitCategory category, filterType typeFilter, F f) const {
553  std::pair<uint8_t, uint8_t> range = getCategoryIndexRange(category);
554  for (int i = range.first; i < range.second; i++) {
555  uint16_t pattern = getHitPatternByAbsoluteIndex(i);
556  // f() return false to ask to stop looping
557  if (typeFilter(pattern) && !f(pattern)) {
558  break;
559  }
560  }
561  }
562 
563  inline int HitPattern::countTypedHits(HitCategory category, filterType typeFilter, filterType filter) const {
564  int count = 0;
565  std::pair<uint8_t, uint8_t> range = getCategoryIndexRange(category);
566  for (int i = range.first; i < range.second; ++i) {
567  uint16_t pattern = getHitPatternByAbsoluteIndex(i);
568  if (typeFilter(pattern) && filter(pattern)) {
569  ++count;
570  }
571  }
572  return count;
573  }
574 
575  inline bool HitPattern::pixelHitFilter(uint16_t pattern) {
576  if UNLIKELY (!trackerHitFilter(pattern)) {
577  return false;
578  }
579 
580  uint32_t substructure = getSubStructure(pattern);
581  return (substructure == PixelSubdetector::PixelBarrel || substructure == PixelSubdetector::PixelEndcap);
582  }
583 
584  inline bool HitPattern::pixelBarrelHitFilter(uint16_t pattern) {
585  if UNLIKELY (!trackerHitFilter(pattern)) {
586  return false;
587  }
588 
589  uint32_t substructure = getSubStructure(pattern);
590  return (substructure == PixelSubdetector::PixelBarrel);
591  }
592 
593  inline bool HitPattern::pixelEndcapHitFilter(uint16_t pattern) {
594  if UNLIKELY (!trackerHitFilter(pattern)) {
595  return false;
596  }
597 
598  uint32_t substructure = getSubStructure(pattern);
599  return (substructure == PixelSubdetector::PixelEndcap);
600  }
601 
602  inline bool HitPattern::stripHitFilter(uint16_t pattern) {
603  return pattern > minStripWord && pattern <= maxTrackerWord;
604  }
605 
606  inline bool HitPattern::stripSubdetectorHitFilter(uint16_t pattern, StripSubdetector::SubDetector substructure) {
607  if UNLIKELY (!trackerHitFilter(pattern)) {
608  return false;
609  }
610 
611  return substructure == getSubStructure(pattern);
612  }
613 
614  inline bool HitPattern::stripTIBHitFilter(uint16_t pattern) {
615  return stripSubdetectorHitFilter(pattern, StripSubdetector::TIB);
616  }
617 
618  inline bool HitPattern::stripTIDHitFilter(uint16_t pattern) {
619  return stripSubdetectorHitFilter(pattern, StripSubdetector::TID);
620  }
621 
622  inline bool HitPattern::stripTOBHitFilter(uint16_t pattern) {
623  return stripSubdetectorHitFilter(pattern, StripSubdetector::TOB);
624  }
625 
626  inline bool HitPattern::stripTECHitFilter(uint16_t pattern) {
627  return stripSubdetectorHitFilter(pattern, StripSubdetector::TEC);
628  }
629 
630  inline bool HitPattern::muonDTHitFilter(uint16_t pattern) {
631  if UNLIKELY (!muonHitFilter(pattern)) {
632  return false;
633  }
634 
635  uint32_t substructure = getSubStructure(pattern);
636  return (substructure == (uint32_t)MuonSubdetId::DT);
637  }
638 
639  inline bool HitPattern::muonCSCHitFilter(uint16_t pattern) {
640  if UNLIKELY (!muonHitFilter(pattern)) {
641  return false;
642  }
643 
644  uint32_t substructure = getSubStructure(pattern);
645  return (substructure == (uint32_t)MuonSubdetId::CSC);
646  }
647 
648  inline bool HitPattern::muonRPCHitFilter(uint16_t pattern) {
649  if UNLIKELY (!muonHitFilter(pattern)) {
650  return false;
651  }
652 
653  uint32_t substructure = getSubStructure(pattern);
654  return (substructure == (uint32_t)MuonSubdetId::RPC);
655  }
656 
657  inline bool HitPattern::muonGEMHitFilter(uint16_t pattern) {
658  if UNLIKELY (!muonHitFilter(pattern)) {
659  return false;
660  }
661 
662  uint32_t substructure = getSubStructure(pattern);
663  return (substructure == (uint32_t)MuonSubdetId::GEM);
664  }
665 
666  inline bool HitPattern::muonME0HitFilter(uint16_t pattern) {
667  if UNLIKELY (!muonHitFilter(pattern))
668  return false;
669  uint16_t substructure = getSubStructure(pattern);
670  return (substructure == (uint16_t)MuonSubdetId::ME0);
671  }
672 
673  inline bool HitPattern::trackerHitFilter(uint16_t pattern) {
674  return pattern > minTrackerWord && pattern <= maxTrackerWord;
675  }
676 
677  inline bool HitPattern::muonHitFilter(uint16_t pattern) {
678  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
679  return false;
680  }
681 
682  return (((pattern >> SubDetectorOffset) & SubDetectorMask) == 0);
683  }
684 
685  inline bool HitPattern::timingBTLHitFilter(uint16_t pattern) {
686  if UNLIKELY (!timingHitFilter(pattern))
687  return false;
688  uint16_t substructure = getSubStructure(pattern);
689  return (substructure == (uint16_t)MTDDetId::BTL);
690  }
691 
692  inline bool HitPattern::timingETLHitFilter(uint16_t pattern) {
693  if UNLIKELY (!timingHitFilter(pattern))
694  return false;
695  uint16_t substructure = getSubStructure(pattern);
696  return (substructure == (uint16_t)MTDDetId::ETL);
697  }
698 
699  inline bool HitPattern::timingHitFilter(uint16_t pattern) {
700  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
701  return false;
702  }
703 
704  return (((pattern >> SubDetectorOffset) & SubDetectorMask) == 2);
705  }
706 
707  inline uint32_t HitPattern::getSubStructure(uint16_t pattern) {
708  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
709  return NULL_RETURN;
710  }
711 
712  return ((pattern >> SubstrOffset) & SubstrMask);
713  }
714 
715  inline uint32_t HitPattern::getLayer(uint16_t pattern) { return HitPattern::getSubSubStructure(pattern); }
716 
717  inline uint32_t HitPattern::getSubSubStructure(uint16_t pattern) {
718  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
719  return NULL_RETURN;
720  }
721 
722  return ((pattern >> LayerOffset) & LayerMask);
723  }
724 
725  inline uint32_t HitPattern::getSubDetector(uint16_t pattern) {
726  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
727  return NULL_RETURN;
728  }
729 
730  return ((pattern >> SubDetectorOffset) & SubDetectorMask);
731  }
732 
733  inline uint32_t HitPattern::getSide(uint16_t pattern) {
734  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
735  return NULL_RETURN;
736  }
737 
738  return (pattern >> SideOffset) & SideMask;
739  }
740 
741  inline uint32_t HitPattern::getHitType(uint16_t pattern) {
742  if UNLIKELY (pattern == HitPattern::EMPTY_PATTERN) {
743  return NULL_RETURN;
744  }
745 
746  return ((pattern >> HitTypeOffset) & HitTypeMask);
747  }
748 
749  inline uint16_t HitPattern::getMuonStation(uint16_t pattern) { return (getSubSubStructure(pattern) >> 2) + 1; }
750 
751  inline uint16_t HitPattern::getDTSuperLayer(uint16_t pattern) { return (getSubSubStructure(pattern) & 3); }
752 
753  inline uint16_t HitPattern::getCSCRing(uint16_t pattern) { return (getSubSubStructure(pattern) & 3) + 1; }
754 
755  inline uint16_t HitPattern::getRPCLayer(uint16_t pattern) {
756  uint16_t subSubStructure = getSubSubStructure(pattern);
757  uint16_t stat = subSubStructure >> 2;
758 
759  if LIKELY (stat <= 1) {
760  return ((subSubStructure >> 1) & 1) + 1;
761  }
762 
763  return 0;
764  }
765 
766  inline uint16_t HitPattern::getRPCregion(uint16_t pattern) { return getSubSubStructure(pattern) & 1; }
767 
769  inline uint16_t HitPattern::getGEMStation(uint16_t pattern)
770 
771  {
772  uint16_t sss = getSubSubStructure(pattern), stat = sss >> 1;
773  return stat + 1;
774  }
775 
777  inline uint16_t HitPattern::getBTLModType(uint16_t pattern) { return getSubSubStructure(pattern); }
778 
779  inline uint16_t HitPattern::getETLRing(uint16_t pattern) { return getSubSubStructure(pattern); }
780 
781  inline uint16_t HitPattern::getGEMLayer(uint16_t pattern) { return (getSubSubStructure(pattern) & 1) + 1; }
782 
783  inline bool HitPattern::validHitFilter(uint16_t pattern) { return getHitType(pattern) == HitPattern::VALID; }
784 
785  inline bool HitPattern::missingHitFilter(uint16_t pattern) { return getHitType(pattern) == HitPattern::MISSING; }
786 
787  inline bool HitPattern::inactiveHitFilter(uint16_t pattern) { return getHitType(pattern) == HitPattern::INACTIVE; }
788 
789  inline bool HitPattern::badHitFilter(uint16_t pattern) { return getHitType(pattern) == HitPattern::BAD; }
790 
791  inline int HitPattern::numberOfAllHits(HitCategory category) const {
792  std::pair<uint8_t, uint8_t> range = getCategoryIndexRange(category);
793  return range.second - range.first;
794  }
795 
796  inline int HitPattern::numberOfAllTrackerHits(HitCategory category) const {
797  return countHits(category, trackerHitFilter);
798  }
799 
800  inline int HitPattern::numberOfMuonHits() const { return countHits(TRACK_HITS, muonHitFilter); }
801 
802  inline int HitPattern::numberOfTimingHits() const { return countHits(TRACK_HITS, timingHitFilter); }
803 
804  inline int HitPattern::numberOfValidHits() const { return countHits(TRACK_HITS, validHitFilter); }
805 
806  inline int HitPattern::numberOfValidTrackerHits() const {
807  return countTypedHits(TRACK_HITS, validHitFilter, trackerHitFilter);
808  }
809 
810  inline int HitPattern::numberOfValidMuonHits() const {
811  return countTypedHits(TRACK_HITS, validHitFilter, muonHitFilter);
812  }
813 
814  inline int HitPattern::numberOfValidTimingHits() const {
815  return countTypedHits(TRACK_HITS, validHitFilter, timingHitFilter);
816  }
817 
818  inline int HitPattern::numberOfValidPixelHits() const {
819  return countTypedHits(TRACK_HITS, validHitFilter, pixelHitFilter);
820  }
821 
822  inline int HitPattern::numberOfValidPixelBarrelHits() const {
823  return countTypedHits(TRACK_HITS, validHitFilter, pixelBarrelHitFilter);
824  }
825 
826  inline int HitPattern::numberOfValidPixelEndcapHits() const {
827  return countTypedHits(TRACK_HITS, validHitFilter, pixelEndcapHitFilter);
828  }
829 
830  inline int HitPattern::numberOfValidStripHits() const {
831  return countTypedHits(TRACK_HITS, validHitFilter, stripHitFilter);
832  }
833 
834  inline int HitPattern::numberOfValidStripTIBHits() const {
835  return countTypedHits(TRACK_HITS, validHitFilter, stripTIBHitFilter);
836  }
837 
838  inline int HitPattern::numberOfValidStripTIDHits() const {
839  return countTypedHits(TRACK_HITS, validHitFilter, stripTIDHitFilter);
840  }
841 
842  inline int HitPattern::numberOfValidStripTOBHits() const {
843  return countTypedHits(TRACK_HITS, validHitFilter, stripTOBHitFilter);
844  }
845 
846  inline int HitPattern::numberOfValidStripTECHits() const {
847  return countTypedHits(TRACK_HITS, validHitFilter, stripTECHitFilter);
848  }
849 
850  inline int HitPattern::numberOfValidMuonDTHits() const {
851  return countTypedHits(TRACK_HITS, validHitFilter, muonDTHitFilter);
852  }
853 
854  inline int HitPattern::numberOfValidMuonCSCHits() const {
855  return countTypedHits(TRACK_HITS, validHitFilter, muonCSCHitFilter);
856  }
857 
858  inline int HitPattern::numberOfValidMuonRPCHits() const {
859  return countTypedHits(TRACK_HITS, validHitFilter, muonRPCHitFilter);
860  }
861 
862  inline int HitPattern::numberOfValidMuonGEMHits() const {
863  return countTypedHits(TRACK_HITS, validHitFilter, muonGEMHitFilter);
864  }
865 
866  inline int HitPattern::numberOfValidMuonME0Hits() const {
867  return countTypedHits(TRACK_HITS, validHitFilter, muonME0HitFilter);
868  }
869 
870  inline int HitPattern::numberOfValidTimingBTLHits() const {
871  return countTypedHits(TRACK_HITS, validHitFilter, timingBTLHitFilter);
872  }
873 
874  inline int HitPattern::numberOfValidTimingETLHits() const {
875  return countTypedHits(TRACK_HITS, validHitFilter, timingETLHitFilter);
876  }
877 
878  inline int HitPattern::numberOfLostHits(HitCategory category) const { return countHits(category, missingHitFilter); }
879 
880  inline int HitPattern::numberOfLostTrackerHits(HitCategory category) const {
881  return countTypedHits(category, missingHitFilter, trackerHitFilter);
882  }
883 
884  inline int HitPattern::numberOfLostMuonHits() const {
885  return countTypedHits(TRACK_HITS, missingHitFilter, muonHitFilter);
886  }
887 
888  inline int HitPattern::numberOfLostTimingHits() const {
889  return countTypedHits(TRACK_HITS, missingHitFilter, timingHitFilter);
890  }
891 
892  inline int HitPattern::numberOfLostTimingBTLHits() const {
893  return countTypedHits(TRACK_HITS, missingHitFilter, timingBTLHitFilter);
894  }
895 
896  inline int HitPattern::numberOfLostTimingETLHits() const {
897  return countTypedHits(TRACK_HITS, missingHitFilter, timingETLHitFilter);
898  }
899 
900  inline int HitPattern::numberOfLostPixelHits(HitCategory category) const {
901  return countTypedHits(category, missingHitFilter, pixelHitFilter);
902  }
903 
904  inline int HitPattern::numberOfLostPixelBarrelHits(HitCategory category) const {
905  return countTypedHits(category, missingHitFilter, pixelBarrelHitFilter);
906  }
907 
908  inline int HitPattern::numberOfLostPixelEndcapHits(HitCategory category) const {
909  return countTypedHits(category, missingHitFilter, pixelEndcapHitFilter);
910  }
911 
912  inline int HitPattern::numberOfLostStripHits(HitCategory category) const {
913  return countTypedHits(category, missingHitFilter, stripHitFilter);
914  }
915 
916  inline int HitPattern::numberOfLostStripTIBHits(HitCategory category) const {
917  return countTypedHits(category, missingHitFilter, stripTIBHitFilter);
918  }
919 
920  inline int HitPattern::numberOfLostStripTIDHits(HitCategory category) const {
921  return countTypedHits(category, missingHitFilter, stripTIDHitFilter);
922  }
923 
924  inline int HitPattern::numberOfLostStripTOBHits(HitCategory category) const {
925  return countTypedHits(category, missingHitFilter, stripTOBHitFilter);
926  }
927 
928  inline int HitPattern::numberOfLostStripTECHits(HitCategory category) const {
929  return countTypedHits(category, missingHitFilter, stripTECHitFilter);
930  }
931 
932  inline int HitPattern::numberOfLostMuonDTHits() const {
933  return countTypedHits(TRACK_HITS, missingHitFilter, muonDTHitFilter);
934  }
935 
936  inline int HitPattern::numberOfLostMuonCSCHits() const {
937  return countTypedHits(TRACK_HITS, missingHitFilter, muonCSCHitFilter);
938  }
939 
940  inline int HitPattern::numberOfLostMuonRPCHits() const {
941  return countTypedHits(TRACK_HITS, missingHitFilter, muonRPCHitFilter);
942  }
943 
944  inline int HitPattern::numberOfLostMuonGEMHits() const {
945  return countTypedHits(TRACK_HITS, missingHitFilter, muonGEMHitFilter);
946  }
947 
948  inline int HitPattern::numberOfLostMuonME0Hits() const {
949  return countTypedHits(TRACK_HITS, missingHitFilter, muonME0HitFilter);
950  }
951 
952  inline int HitPattern::numberOfBadHits() const { return countHits(TRACK_HITS, badHitFilter); }
953 
954  inline int HitPattern::numberOfBadMuonHits() const {
955  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonHitFilter);
956  }
957 
958  inline int HitPattern::numberOfBadMuonDTHits() const {
959  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonDTHitFilter);
960  }
961 
962  inline int HitPattern::numberOfBadMuonCSCHits() const {
963  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonCSCHitFilter);
964  }
965 
966  inline int HitPattern::numberOfBadMuonRPCHits() const {
967  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonRPCHitFilter);
968  }
969 
970  inline int HitPattern::numberOfBadMuonGEMHits() const {
971  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonGEMHitFilter);
972  }
973 
974  inline int HitPattern::numberOfBadMuonME0Hits() const {
975  return countTypedHits(TRACK_HITS, inactiveHitFilter, muonME0HitFilter);
976  }
977 
978  inline int HitPattern::numberOfInactiveHits() const { return countHits(TRACK_HITS, inactiveHitFilter); }
979 
980  inline int HitPattern::numberOfInactiveTrackerHits() const {
981  return countTypedHits(TRACK_HITS, inactiveHitFilter, trackerHitFilter);
982  }
983 
984  inline int HitPattern::trackerLayersWithMeasurementOld() const {
985  return pixelLayersWithMeasurement() + stripLayersWithMeasurement();
986  }
987 
988  inline int HitPattern::pixelLayersWithMeasurementOld() const {
989  return pixelBarrelLayersWithMeasurement() + pixelEndcapLayersWithMeasurement();
990  }
991 
992  inline int HitPattern::stripLayersWithMeasurement() const {
993  return stripTIBLayersWithMeasurement() + stripTIDLayersWithMeasurement() + stripTOBLayersWithMeasurement() +
994  stripTECLayersWithMeasurement();
995  }
996 
997  inline int HitPattern::trackerLayersWithoutMeasurementOld(HitCategory category) const {
998  return pixelLayersWithoutMeasurement(category) + stripLayersWithoutMeasurement(category);
999  }
1000 
1001  inline int HitPattern::pixelLayersWithoutMeasurement(HitCategory category) const {
1002  return pixelBarrelLayersWithoutMeasurement(category) + pixelEndcapLayersWithoutMeasurement(category);
1003  }
1004 
1005  inline int HitPattern::stripLayersWithoutMeasurement(HitCategory category) const {
1006  return stripTIBLayersWithoutMeasurement(category) + stripTIDLayersWithoutMeasurement(category) +
1007  stripTOBLayersWithoutMeasurement(category) + stripTECLayersWithoutMeasurement(category);
1008  }
1009 
1010  inline int HitPattern::trackerLayersTotallyOffOrBad(HitCategory category) const {
1011  return pixelLayersTotallyOffOrBad(category) + stripLayersTotallyOffOrBad(category);
1012  }
1013 
1014  inline int HitPattern::pixelLayersTotallyOffOrBad(HitCategory category) const {
1015  return pixelBarrelLayersTotallyOffOrBad(category) + pixelEndcapLayersTotallyOffOrBad(category);
1016  }
1017 
1018  inline int HitPattern::stripLayersTotallyOffOrBad(HitCategory category) const {
1019  return stripTIBLayersTotallyOffOrBad(category) + stripTIDLayersTotallyOffOrBad(category) +
1020  stripTOBLayersTotallyOffOrBad(category) + stripTECLayersTotallyOffOrBad(category);
1021  }
1022 
1023  inline int HitPattern::trackerLayersNull() const { return pixelLayersNull() + stripLayersNull(); }
1024 
1025  inline int HitPattern::pixelLayersNull() const { return pixelBarrelLayersNull() + pixelEndcapLayersNull(); }
1026 
1027  inline int HitPattern::stripLayersNull() const {
1028  return stripTIBLayersNull() + stripTIDLayersNull() + stripTOBLayersNull() + stripTECLayersNull();
1029  }
1030 
1031  inline int HitPattern::muonStationsWithValidHits() const { return muonStations(0, 0); }
1032 
1033  inline int HitPattern::muonStationsWithBadHits() const { return muonStations(0, 3); }
1034 
1035  inline int HitPattern::muonStationsWithAnyHits() const { return muonStations(0, -1); }
1036 
1037  inline int HitPattern::dtStationsWithValidHits() const { return muonStations(1, 0); }
1038 
1039  inline int HitPattern::dtStationsWithBadHits() const { return muonStations(1, 3); }
1040 
1041  inline int HitPattern::dtStationsWithAnyHits() const { return muonStations(1, -1); }
1042 
1043  inline int HitPattern::cscStationsWithValidHits() const { return muonStations(2, 0); }
1044 
1045  inline int HitPattern::cscStationsWithBadHits() const { return muonStations(2, 3); }
1046 
1047  inline int HitPattern::cscStationsWithAnyHits() const { return muonStations(2, -1); }
1048 
1049  inline int HitPattern::rpcStationsWithValidHits() const { return muonStations(3, 0); }
1050 
1051  inline int HitPattern::rpcStationsWithBadHits() const { return muonStations(3, 3); }
1052 
1053  inline int HitPattern::rpcStationsWithAnyHits() const { return muonStations(3, -1); }
1054 
1055  inline int HitPattern::gemStationsWithValidHits() const { return muonStations(4, 0); }
1056 
1057  inline int HitPattern::gemStationsWithBadHits() const { return muonStations(4, 3); }
1058 
1059  inline int HitPattern::gemStationsWithAnyHits() const { return muonStations(4, -1); }
1060 
1061  inline int HitPattern::me0StationsWithValidHits() const { return muonStations(5, 0); }
1062 
1063  inline int HitPattern::me0StationsWithBadHits() const { return muonStations(5, 3); }
1064 
1065  inline int HitPattern::me0StationsWithAnyHits() const { return muonStations(5, -1); }
1066 
1067  inline int HitPattern::innermostMuonStationWithValidHits() const { return innermostMuonStationWithHits(0); }
1068 
1069  inline int HitPattern::innermostMuonStationWithBadHits() const { return innermostMuonStationWithHits(3); }
1070 
1071  inline int HitPattern::innermostMuonStationWithAnyHits() const { return innermostMuonStationWithHits(-1); }
1072 
1073  inline int HitPattern::outermostMuonStationWithValidHits() const { return outermostMuonStationWithHits(0); }
1074 
1075  inline int HitPattern::outermostMuonStationWithBadHits() const { return outermostMuonStationWithHits(3); }
1076 
1077  inline int HitPattern::outermostMuonStationWithAnyHits() const { return outermostMuonStationWithHits(-1); }
1078 
1079  template <int N = HitPattern::MaxHits>
1080  struct PatternSet {
1081  static constexpr int MaxHits = N;
1082  unsigned char hit[N];
1083  unsigned char nhit;
1084 
1085  unsigned char const *begin() const { return hit; }
1086 
1087  unsigned char const *end() const { return hit + nhit; }
1088 
1089  unsigned char *begin() { return hit; }
1090 
1091  unsigned char *end() { return hit + nhit; }
1092 
1093  int size() const { return nhit; }
1094 
1095  unsigned char operator[](int i) const { return hit[i]; }
1096 
1097  PatternSet() : nhit(0) {}
1098 
1099  PatternSet(HitPattern::HitCategory category, HitPattern const &hp) { fill(category, hp); }
1100 
1101  void fill(HitPattern::HitCategory category, HitPattern const &hp) {
1102  int lhit = 0;
1103  auto unpack = [&lhit, this](uint16_t pattern) -> bool {
1104  unsigned char p = 255 & (pattern >> 3);
1105  hit[lhit++] = p;
1106 
1107  // bouble sort
1108  if (lhit > 1) {
1109  for (auto h = hit + lhit - 1; h != hit; --h) {
1110  if ((*(h - 1)) <= p) {
1111  break;
1112  }
1113  (*h) = *(h - 1);
1114  *(h - 1) = p;
1115  }
1116  }
1117  return lhit < MaxHits;
1118  };
1119 
1120  hp.call(category, HitPattern::validHitFilter, unpack);
1121  nhit = lhit;
1122  }
1123  };
1124 
1125  template <int N>
1126  inline PatternSet<N> commonHits(PatternSet<N> const &p1, PatternSet<N> const &p2) {
1127  PatternSet<N> comm;
1128  comm.nhit = std::set_intersection(p1.begin(), p1.end(), p2.begin(), p2.end(), comm.begin()) - comm.begin();
1129  return comm;
1130  }
1131 
1132 } // namespace reco
1133 
1134 #endif
reco::HitPattern::appendHit
bool appendHit(const TrackingRecHit &hit, const TrackerTopology &ttopo)
Definition: HitPattern.cc:201
reco::HitPattern::operator=
HitPattern & operator=(const HitPattern &other)
Definition: HitPattern.cc:37
MuonSubdetId::GEM
static constexpr int GEM
Definition: MuonSubdetId.h:14
Likely.h
reco::HitPattern::numberOfBadMuonDTHits
int numberOfBadMuonDTHits() const
Definition: HitPattern.h:958
reco::HitPattern::MISSING_OUTER_HITS
Definition: HitPattern.h:155
reco::HitPattern::getHitPattern
uint16_t getHitPattern(HitCategory category, int position) const
Definition: HitPattern.h:531
reco::HitPattern::getSubStructure
static uint32_t getSubStructure(uint16_t pattern)
Definition: HitPattern.h:707
reco::HitPattern::cscStationsWithAnyHits
int cscStationsWithAnyHits() const
Definition: HitPattern.h:1047
FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
reco::HitPattern::numberOfValidStripTECHits
int numberOfValidStripTECHits() const
Definition: HitPattern.h:846
taus_updatedMVAIds_cff.category
category
Definition: taus_updatedMVAIds_cff.py:31
reco::HitPattern::inactiveHitFilter
static bool inactiveHitFilter(uint16_t pattern)
Definition: HitPattern.h:787
MuonSubdetId::CSC
static constexpr int CSC
Definition: MuonSubdetId.h:12
reco::HitPattern::HIT_LENGTH
const static unsigned short HIT_LENGTH
Definition: HitPattern.h:157
mps_fire.i
i
Definition: mps_fire.py:428
PixelSubdetector.h
reco::HitPattern::LayerOffset
const static unsigned short LayerOffset
Definition: HitPattern.h:448
reco::HitPattern::numberOfLostMuonDTHits
int numberOfLostMuonDTHits() const
Definition: HitPattern.h:932
reco::HitPattern::numberOfValidTimingBTLHits
int numberOfValidTimingBTLHits() const
Definition: HitPattern.h:870
reco::HitPattern::getSubDetector
static uint32_t getSubDetector(uint16_t pattern)
Definition: HitPattern.h:725
PixelSubdetector::PixelEndcap
Definition: PixelSubdetector.h:11
reco::HitPattern::beginTrackHits
uint8_t beginTrackHits
Definition: HitPattern.h:493
reco::HitPattern::pixelEndcapLayersWithMeasurement
int pixelEndcapLayersWithMeasurement() const
Definition: HitPattern.cc:563
reco::HitPattern::filterType
bool filterType(uint16_t)
Definition: HitPattern.h:473
reco::HitPattern::innermostMuonStationWithBadHits
int innermostMuonStationWithBadHits() const
Definition: HitPattern.h:1069
PixelSubdetector::PixelBarrel
Definition: PixelSubdetector.h:11
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:594
reco::HitPattern::numberOfAllTrackerHits
int numberOfAllTrackerHits(HitCategory category) const
Definition: HitPattern.h:796
reco::HitPattern::trackerLayersWithoutMeasurement
int trackerLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:529
reco::HitPattern::numberOfTimingHits
int numberOfTimingHits() const
Definition: HitPattern.h:802
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
reco::HitPattern::numberOfLostStripTIBHits
int numberOfLostStripTIBHits(HitCategory category) const
Definition: HitPattern.h:916
reco::HitPattern::stripTIBLayersWithoutMeasurement
int stripTIBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:636
reco::HitPattern::beginInner
uint8_t beginInner
Definition: HitPattern.h:495
reco::HitPattern::muonRPCHitFilter
static bool muonRPCHitFilter(uint16_t pattern)
Definition: HitPattern.h:648
reco::HitPattern::numberOfValidTOBLayersWithMonoAndStereo
int numberOfValidTOBLayersWithMonoAndStereo(uint32_t layer=0) const
Definition: HitPattern.cc:411
reco::PatternSet::PatternSet
PatternSet(HitPattern::HitCategory category, HitPattern const &hp)
Definition: HitPattern.h:1099
reco::HitPattern::SideMask
const static unsigned short SideMask
Definition: HitPattern.h:445
TrackerTopology
Definition: TrackerTopology.h:16
reco::PatternSet::fill
void fill(HitPattern::HitCategory category, HitPattern const &hp)
Definition: HitPattern.h:1101
bJpsiMuMuTrigSettings_cff.filterType
filterType
Definition: bJpsiMuMuTrigSettings_cff.py:3
reco::HitPattern::validHitFilter
static bool validHitFilter(uint16_t pattern)
Definition: HitPattern.h:783
reco::HitPattern::numberOfLostTimingBTLHits
int numberOfLostTimingBTLHits() const
Definition: HitPattern.h:892
gather_cfg.cout
cout
Definition: gather_cfg.py:144
reco::HitPattern::HitTypeOffset
const static unsigned short HitTypeOffset
Definition: HitPattern.h:440
reco::HitPattern::clear
void clear()
Definition: HitPattern.cc:58
reco::HitPattern::numberOfLostStripTECHits
int numberOfLostStripTECHits(HitCategory category) const
Definition: HitPattern.h:928
reco::HitPattern::muonHitFilter
static bool muonHitFilter(uint16_t pattern)
Definition: HitPattern.h:677
reco::HitPattern::trackerHitFilter
static bool trackerHitFilter(uint16_t pattern)
Definition: HitPattern.h:673
reco::HitPattern::dtStationsWithValidHits
int dtStationsWithValidHits() const
Definition: HitPattern.h:1037
cms::cuda::stream
uint32_t const T *__restrict__ const uint32_t *__restrict__ int32_t int Histo::index_type cudaStream_t stream
Definition: HistoContainer.h:51
reco::HitPattern::numberOfBadMuonGEMHits
int numberOfBadMuonGEMHits() const
Definition: HitPattern.h:970
reco::HitPattern::numberOfMuonHits
int numberOfMuonHits() const
Definition: HitPattern.h:800
reco::HitPattern::gemStationsWithAnyHits
int gemStationsWithAnyHits() const
Definition: HitPattern.h:1059
reco::HitPattern::stripTIBLayersNull
int stripTIBLayersNull() const
Definition: HitPattern.cc:760
reco::HitPattern::getGEMLayer
static uint16_t getGEMLayer(uint16_t pattern)
GEM layer: 1,2. Only valid for muon GEM patterns, of course.
Definition: HitPattern.h:781
reco::HitPattern::MUON_HIT
Definition: HitPattern.h:151
reco::HitPattern::numberOfLostTimingHits
int numberOfLostTimingHits() const
Definition: HitPattern.h:888
reco::HitPattern::muonGEMHitFilter
static bool muonGEMHitFilter(uint16_t pattern)
Definition: HitPattern.h:657
reco::HitPattern::stripLayersWithoutMeasurement
int stripLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.h:1005
reco::HitPattern::HitPattern
HitPattern()
Definition: HitPattern.cc:19
reco::HitPattern::beginOuter
uint8_t beginOuter
Definition: HitPattern.h:497
reco::HitPattern::trackerLayersWithoutMeasurementOld
int trackerLayersWithoutMeasurementOld(HitCategory category) const
Definition: HitPattern.h:997
TrackingRecHitFwd.h
reco::PatternSet::hit
unsigned char hit[N]
Definition: HitPattern.h:1082
reco::HitPattern::getCSCRing
static uint16_t getCSCRing(uint16_t pattern)
CSC ring (1-4). Only valid for muon CSC patterns, of course.
Definition: HitPattern.h:753
reco::HitPattern::muonCSCHitFilter
static bool muonCSCHitFilter(uint16_t pattern)
Definition: HitPattern.h:639
reco::HitPattern::dtStationsWithBadHits
int dtStationsWithBadHits() const
Definition: HitPattern.h:1039
reco::HitPattern::stripTOBHitFilter
static bool stripTOBHitFilter(uint16_t pattern)
Definition: HitPattern.h:622
reco::HitPattern::fillNewHitPatternWithOldHitPattern_v12
static bool fillNewHitPatternWithOldHitPattern_v12(const uint16_t oldHitPattern[], uint8_t hitCount, uint8_t beginTrackHits, uint8_t endTrackHits, uint8_t beginInner, uint8_t endInner, uint8_t beginOuter, uint8_t endOuter, reco::HitPattern *newObj)
Definition: rootio_HitPattern.cc:60
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46
reco::HitPattern::rpcStationsWithBadHits
int rpcStationsWithBadHits() const
Definition: HitPattern.h:1051
reco::PatternSet::PatternSet
PatternSet()
Definition: HitPattern.h:1097
h
The Signals That Services Can Subscribe To This is based on ActivityRegistry h
Helper function to determine trigger accepts.
Definition: Activities.doc:4
reco::HitPattern::insertExpectedInnerHit
bool insertExpectedInnerHit(const uint16_t pattern)
Definition: HitPattern.cc:981
reco::HitPattern::numberOfLostHits
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:878
reco::HitPattern::numberOfValidPixelEndcapHits
int numberOfValidPixelEndcapHits() const
Definition: HitPattern.h:826
reco::HitPattern::stripTIBHitFilter
static bool stripTIBHitFilter(uint16_t pattern)
Definition: HitPattern.h:614
reco::HitPattern::trackerLayersTotallyOffOrBad
int trackerLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.h:1010
reco::HitPattern::maxTrackerWord
const static unsigned short maxTrackerWord
Definition: HitPattern.h:460
reco::HitPattern::getRPCLayer
static uint16_t getRPCLayer(uint16_t pattern)
RPC layer: for station 1 and 2, layer = 1(inner) or 2(outer); for station 3, 4 layer is always 0....
Definition: HitPattern.h:755
reco::HitPattern::me0StationsWithBadHits
int me0StationsWithBadHits() const
Definition: HitPattern.h:1063
reco::HitPattern::getHitType
static uint32_t getHitType(uint16_t pattern)
Definition: HitPattern.h:741
reco::HitPattern::getSide
static uint32_t getSide(uint16_t pattern)
Definition: HitPattern.h:733
reco::HitPattern::numberOfLostPixelEndcapHits
int numberOfLostPixelEndcapHits(HitCategory category) const
Definition: HitPattern.h:908
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:489
reco::HitPattern::stripTOBLayersWithoutMeasurement
int stripTOBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:656
reco::HitPattern::STEREO
Definition: HitPattern.h:149
edm::Ref
Definition: AssociativeIterator.h:58
reco::HitPattern::numberOfInactiveTrackerHits
int numberOfInactiveTrackerHits() const
Definition: HitPattern.h:980
F
static uInt32 F(BLOWFISH_CTX *ctx, uInt32 x)
Definition: blowfish.cc:163
reco::HitPattern::pixelBarrelLayersWithMeasurement
int pixelBarrelLayersWithMeasurement() const
Definition: HitPattern.cc:552
reco::HitPattern::numberOfBadHits
int numberOfBadHits() const
Definition: HitPattern.h:952
reco::HitPattern::numberOfLostMuonHits
int numberOfLostMuonHits() const
Definition: HitPattern.h:884
DetId
Definition: DetId.h:17
UNLIKELY
#define UNLIKELY(x)
Definition: Likely.h:21
reco::HitPattern::minTrackerWord
const static unsigned short minTrackerWord
Definition: HitPattern.h:459
reco::HitPattern::outermostMuonStationWithBadHits
int outermostMuonStationWithBadHits() const
Definition: HitPattern.h:1075
reco::HitPattern::numberOfValidPixelBarrelHits
int numberOfValidPixelBarrelHits() const
Definition: HitPattern.h:822
reco::HitPattern
Definition: HitPattern.h:147
TrackingRecHit.h
test
Definition: SmallWORMDict.h:13
reco::HitPattern::pixelLayersNull
int pixelLayersNull() const
Definition: HitPattern.h:1025
reco::HitPattern::pixelBarrelLayersNull
int pixelBarrelLayersNull() const
Definition: HitPattern.cc:738
reco::HitPattern::trackerLayersWithMeasurementOld
int trackerLayersWithMeasurementOld() const
Definition: HitPattern.h:984
reco::HitPattern::stripLayersWithMeasurement
int stripLayersWithMeasurement() const
Definition: HitPattern.h:992
reco::HitPattern::TRACKER_HIT
Definition: HitPattern.h:151
reco::HitPattern::numberOfLostStripTIDHits
int numberOfLostStripTIDHits(HitCategory category) const
Definition: HitPattern.h:920
reco::HitPattern::endOuter
uint8_t endOuter
Definition: HitPattern.h:498
trackingPlots.hp
hp
Definition: trackingPlots.py:1231
reco::HitPattern::encode
static uint16_t encode(const TrackingRecHit &hit, const TrackerTopology &ttopo)
Definition: HitPattern.cc:74
reco::HitPattern::numberOfValidMuonGEMHits
int numberOfValidMuonGEMHits() const
Definition: HitPattern.h:862
reco::HitPattern::muonStationsWithValidHits
int muonStationsWithValidHits() const
Definition: HitPattern.h:1031
reco::HitPattern::innermostMuonStationWithValidHits
int innermostMuonStationWithValidHits() const
Definition: HitPattern.h:1067
Exhume::I
const std::complex< double > I
Definition: I.h:8
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:510
reco::HitPattern::stripTIDHitFilter
static bool stripTIDHitFilter(uint16_t pattern)
Definition: HitPattern.h:618
reco::HitPattern::minPixelWord
const static unsigned short minPixelWord
Definition: HitPattern.h:461
reco::HitPattern::getTrackerLayerCase
uint32_t getTrackerLayerCase(HitCategory category, uint16_t substr, uint16_t layer) const
Definition: HitPattern.cc:427
reco::HitPattern::timingETLHitFilter
static bool timingETLHitFilter(uint16_t pattern)
Definition: HitPattern.h:692
reco::HitPattern::numberOfLostStripTOBHits
int numberOfLostStripTOBHits(HitCategory category) const
Definition: HitPattern.h:924
reco::HitPattern::getSubSubStructure
static uint32_t getSubSubStructure(uint16_t pattern)
Definition: HitPattern.h:717
TrackingRecHit::Type
Type
Definition: TrackingRecHit.h:45
reco::HitPattern::stripTIDLayersTotallyOffOrBad
int stripTIDLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:708
mps_fire.end
end
Definition: mps_fire.py:242
reco::HitPattern::print
void print(HitCategory category, std::ostream &stream=std::cout) const
Definition: HitPattern.cc:836
p2
double p2[4]
Definition: TauolaWrapper.h:90
reco::PatternSet::operator[]
unsigned char operator[](int i) const
Definition: HitPattern.h:1095
N
#define N
Definition: blowfish.cc:9
reco::HitPattern::minStripWord
const static unsigned short minStripWord
Definition: HitPattern.h:462
reco::HitPattern::numberOfLostTimingETLHits
int numberOfLostTimingETLHits() const
Definition: HitPattern.h:896
reco::HitPattern::pixelEndcapLayersTotallyOffOrBad
int pixelEndcapLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:687
reco::HitPattern::getGEMStation
static uint16_t getGEMStation(uint16_t pattern)
GEM station: 1,2. Only valid for muon GEM patterns, of course.
Definition: HitPattern.h:769
reco::HitPattern::getCategoryIndexRange
std::pair< uint8_t, uint8_t > getCategoryIndexRange(HitCategory category) const
Definition: HitPattern.h:506
submitPVResolutionJobs.count
count
Definition: submitPVResolutionJobs.py:352
reco::HitPattern::hasValidHitInPixelLayer
bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const
Definition: HitPattern.cc:334
StripSubdetector::TIB
static constexpr auto TIB
Definition: StripSubdetector.h:16
trackingPlots.other
other
Definition: trackingPlots.py:1464
reco::HitPattern::missingHitFilter
static bool missingHitFilter(uint16_t pattern)
Definition: HitPattern.h:785
reco::HitPattern::numberOfInactiveTimingHits
int numberOfInactiveTimingHits() const
reco::HitPattern::HIT_DETECTOR_TYPE
HIT_DETECTOR_TYPE
Definition: HitPattern.h:151
reco::HitPattern::muonStationsWithAnyHits
int muonStationsWithAnyHits() const
Definition: HitPattern.h:1035
MuonSubdetId::DT
static constexpr int DT
Definition: MuonSubdetId.h:11
h
reco::HitPattern::pixelBarrelLayersWithoutMeasurement
int pixelBarrelLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:614
reco::HitPattern::insertHit
void insertHit(const uint16_t pattern)
Definition: HitPattern.cc:952
ALCARECOTkAlBeamHalo_cff.filter
filter
Definition: ALCARECOTkAlBeamHalo_cff.py:27
reco::HitPattern::getRPCregion
static uint16_t getRPCregion(uint16_t pattern)
RPC region: 0 = barrel, 1 = endcap. Only valid for muon RPC patterns, of course.
Definition: HitPattern.h:766
reco::HitPattern::VALID
Definition: HitPattern.h:153
reco::HitPattern::numberOfValidMuonCSCHits
int numberOfValidMuonCSCHits() const
Definition: HitPattern.h:854
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
reco::HitPattern::SideOffset
const static unsigned short SideOffset
Definition: HitPattern.h:444
PixelSubdetector::SubDetector
SubDetector
Definition: PixelSubdetector.h:11
reco::HitPattern::pixelHitFilter
static bool pixelHitFilter(uint16_t pattern)
Definition: HitPattern.h:575
reco::HitPattern::stripTIBLayersWithMeasurement
int stripTIBLayersWithMeasurement() const
Definition: HitPattern.cc:574
reco::HitPattern::numberOfDTStationsWithBothViews
int numberOfDTStationsWithBothViews() const
Definition: HitPattern.cc:939
reco::HitPattern::getLayer
static uint32_t getLayer(uint16_t pattern)
Definition: HitPattern.h:715
reco::PatternSet::end
unsigned const char * end() const
Definition: HitPattern.h:1087
reco::HitPattern::printHitPattern
void printHitPattern(HitCategory category, int position, std::ostream &stream) const
Definition: HitPattern.cc:800
reco::HitPattern::numberOfInactiveHits
int numberOfInactiveHits() const
Definition: HitPattern.h:978
reco::HitPattern::cscStationsWithBadHits
int cscStationsWithBadHits() const
Definition: HitPattern.h:1045
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
reco::HitPattern::pixelEndcapHitFilter
static bool pixelEndcapHitFilter(uint16_t pattern)
Definition: HitPattern.h:593
reco::HitPattern::stripTECLayersTotallyOffOrBad
int stripTECLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:728
reco::HitPattern::appendHits
bool appendHits(const I &begin, const I &end, const TrackerTopology &ttopo)
Definition: HitPattern.h:522
reco::HitPattern::innermostMuonStationWithAnyHits
int innermostMuonStationWithAnyHits() const
Definition: HitPattern.h:1071
MTDDetId.h
reco::HitPattern::numberOfLostPixelHits
int numberOfLostPixelHits(HitCategory category) const
Definition: HitPattern.h:900
reco::HitPattern::outermostMuonStationWithValidHits
int outermostMuonStationWithValidHits() const
Definition: HitPattern.h:1073
reco::HitPattern::SubstrMask
const static unsigned short SubstrMask
Definition: HitPattern.h:453
reco::HitPattern::outermostMuonStationWithAnyHits
int outermostMuonStationWithAnyHits() const
Definition: HitPattern.h:1077
reco::HitPattern::pixelLayersTotallyOffOrBad
int pixelLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.h:1014
reco::HitPattern::numberOfBadMuonME0Hits
int numberOfBadMuonME0Hits() const
Definition: HitPattern.h:974
reco::HitPattern::timingHitFilter
static bool timingHitFilter(uint16_t pattern)
Definition: HitPattern.h:699
reco::PatternSet::nhit
unsigned char nhit
Definition: HitPattern.h:1083
position
static int position[264][3]
Definition: ReadPGInfo.cc:289
reco::HitPattern::numberOfValidTECLayersWithMonoAndStereo
int numberOfValidTECLayersWithMonoAndStereo(uint32_t layer=0) const
Definition: HitPattern.cc:423
reco::HitPattern::HIT_TYPE
HIT_TYPE
Definition: HitPattern.h:153
reco::HitPattern::EMPTY_PATTERN
static const uint16_t EMPTY_PATTERN
Definition: HitPattern.h:161
reco::PatternSet::begin
unsigned char * begin()
Definition: HitPattern.h:1089
reco::HitPattern::numberOfBadMuonRPCHits
int numberOfBadMuonRPCHits() const
Definition: HitPattern.h:966
topSingleLeptonDQM_PU_cfi.pattern
pattern
Definition: topSingleLeptonDQM_PU_cfi.py:39
reco::HitPattern::getETLRing
static uint16_t getETLRing(uint16_t pattern)
ETL Ring: 1-12. Only valid for ETL patterns of course.
Definition: HitPattern.h:779
reco::HitPattern::stripTIDLayersWithoutMeasurement
int stripTIDLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:646
reco::HitPattern::pixelLayersWithMeasurementOld
int pixelLayersWithMeasurementOld() const
Definition: HitPattern.h:988
p1
double p1[4]
Definition: TauolaWrapper.h:89
reco::HitPattern::rpcStationsWithValidHits
int rpcStationsWithValidHits() const
Definition: HitPattern.h:1049
reco::HitPattern::numberOfValidStripTIBHits
int numberOfValidStripTIBHits() const
Definition: HitPattern.h:834
reco::HitPattern::getMuonStation
static uint16_t getMuonStation(uint16_t pattern)
Muon station (1-4). Only valid for muon patterns, of course. only for patterns from muon,...
Definition: HitPattern.h:749
reco::HitPattern::numberOfValidTimingETLHits
int numberOfValidTimingETLHits() const
Definition: HitPattern.h:874
MuonSubdetId::ME0
static constexpr int ME0
Definition: MuonSubdetId.h:15
reco::HitPattern::muonStationsWithBadHits
int muonStationsWithBadHits() const
Definition: HitPattern.h:1033
reco::HitPattern::countTypedHits
int countTypedHits(HitCategory category, filterType typeFilter, filterType filter) const
Definition: HitPattern.h:563
reco::HitPattern::numberOfLostMuonCSCHits
int numberOfLostMuonCSCHits() const
Definition: HitPattern.h:936
reco::HitPattern::getHitPatternByAbsoluteIndex
uint16_t getHitPatternByAbsoluteIndex(int position) const
Definition: HitPattern.cc:293
reco::HitPattern::SubDetectorMask
const static unsigned short SubDetectorMask
Definition: HitPattern.h:457
reco::HitPattern::numberOfLostMuonGEMHits
int numberOfLostMuonGEMHits() const
Definition: HitPattern.h:944
MuonSubdetId.h
reco::HitPattern::LayerMask
const static unsigned short LayerMask
Definition: HitPattern.h:449
reco::HitPattern::isStereo
static uint16_t isStereo(DetId i, const TrackerTopology &ttopo)
Definition: HitPattern.cc:852
reco::HitPattern::timingBTLHitFilter
static bool timingBTLHitFilter(uint16_t pattern)
Definition: HitPattern.h:685
reco::HitPattern::badHitFilter
static bool badHitFilter(uint16_t pattern)
Definition: HitPattern.h:789
reco::HitPattern::pixelBarrelLayersTotallyOffOrBad
int pixelBarrelLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:676
reco::HitPattern::TRACK_HITS
Definition: HitPattern.h:155
reco::HitPattern::numberOfValidMuonDTHits
int numberOfValidMuonDTHits() const
Definition: HitPattern.h:850
TrackingRecHit
Definition: TrackingRecHit.h:21
reco::HitPattern::ARRAY_LENGTH
const static unsigned short ARRAY_LENGTH
Definition: HitPattern.h:156
reco::HitPattern::~HitPattern
~HitPattern()
Definition: HitPattern.cc:35
reco::HitPattern::pixelEndcapLayersNull
int pixelEndcapLayersNull() const
Definition: HitPattern.cc:749
reco::HitPattern::getBTLModType
static uint16_t getBTLModType(uint16_t pattern)
BTL Module type: 1,2,3. Only valid for BTL patterns of course.
Definition: HitPattern.h:777
reco::HitPattern::numberOfLostPixelBarrelHits
int numberOfLostPixelBarrelHits(HitCategory category) const
Definition: HitPattern.h:904
reco::HitPattern::numberOfValidStripTIDHits
int numberOfValidStripTIDHits() const
Definition: HitPattern.h:838
reco::HitPattern::stripTOBLayersTotallyOffOrBad
int stripTOBLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:718
MTDDetId::BTL
Definition: MTDDetId.h:27
reco::PatternSet::end
unsigned char * end()
Definition: HitPattern.h:1091
reco::PatternSet::size
int size() const
Definition: HitPattern.h:1093
reco::HitPattern::stripTIDLayersWithMeasurement
int stripTIDLayersWithMeasurement() const
Definition: HitPattern.cc:584
reco::HitPattern::numberOfAllHits
int numberOfAllHits(HitCategory category) const
Definition: HitPattern.h:791
reco::HitPattern::HitCategory
HitCategory
Definition: HitPattern.h:155
StripSubdetector::TEC
static constexpr auto TEC
Definition: StripSubdetector.h:19
reco::HitPattern::trackerLayersNull
int trackerLayersNull() const
Definition: HitPattern.h:1023
reco::HitPattern::gemStationsWithBadHits
int gemStationsWithBadHits() const
Definition: HitPattern.h:1057
reco::HitPattern::rpcStationsWithAnyHits
int rpcStationsWithAnyHits() const
Definition: HitPattern.h:1053
reco::HitPattern::muonStations
int muonStations(int subdet, int hitType) const
subdet = 0(all), 1(DT), 2(CSC), 3(RPC) 4(GEM); hitType=-1(all), 0=valid, 3=bad
Definition: HitPattern.cc:874
reco::PatternSet::begin
unsigned const char * begin() const
Definition: HitPattern.h:1085
reco::HitPattern::me0StationsWithAnyHits
int me0StationsWithAnyHits() const
Definition: HitPattern.h:1065
reco::HitPattern::numberOfValidStripHits
int numberOfValidStripHits() const
Definition: HitPattern.h:830
reco::HitPattern::numberOfValidTrackerHits
int numberOfValidTrackerHits() const
Definition: HitPattern.h:806
reco::HitPattern::numberOfValidStripTOBHits
int numberOfValidStripTOBHits() const
Definition: HitPattern.h:842
reco::HitPattern::appendMuonHit
bool appendMuonHit(const DetId &id, TrackingRecHit::Type hitType)
Definition: HitPattern.cc:277
DetId.h
reco::HitPattern::INACTIVE
Definition: HitPattern.h:153
MuonSubdetId::RPC
static constexpr int RPC
Definition: MuonSubdetId.h:13
reco::HitPattern::stripLayersNull
int stripLayersNull() const
Definition: HitPattern.h:1027
reco::HitPattern::innermostMuonStationWithHits
int innermostMuonStationWithHits(int hitType) const
hitType=-1(all), 0=valid, 3=bad; 0 = no stations at all
Definition: HitPattern.cc:887
LIKELY
#define LIKELY(x)
Definition: Likely.h:20
MTDDetId::ETL
Definition: MTDDetId.h:27
reco::HitPattern::stripSubdetectorHitFilter
static bool stripSubdetectorHitFilter(uint16_t pattern, StripSubdetector::SubDetector substructure)
Definition: HitPattern.h:606
reco::HitPattern::appendTrackerHit
bool appendTrackerHit(uint16_t subdet, uint16_t layer, uint16_t stereo, TrackingRecHit::Type hitType)
Definition: HitPattern.cc:273
reco::HitPattern::numberOfValidMuonHits
int numberOfValidMuonHits() const
Definition: HitPattern.h:810
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
reco::HitPattern::cscStationsWithValidHits
int cscStationsWithValidHits() const
Definition: HitPattern.h:1043
reco::HitPattern::stripLayersTotallyOffOrBad
int stripLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.h:1018
reco::HitPattern::BAD
Definition: HitPattern.h:153
StripSubdetector::TOB
static constexpr auto TOB
Definition: StripSubdetector.h:18
BeamSpotPI::unpack
std::pair< unsigned int, unsigned int > unpack(cond::Time_t since)
Definition: BeamSpotPayloadInspectorHelper.h:23
reco::HitPattern::HitTypeMask
const static unsigned short HitTypeMask
Definition: HitPattern.h:441
reco::HitPattern::stripTECHitFilter
static bool stripTECHitFilter(uint16_t pattern)
Definition: HitPattern.h:626
reco::HitPattern::countHits
int countHits(HitCategory category, filterType filter) const
Definition: HitPattern.h:540
reco::HitPattern::numberOfBadMuonCSCHits
int numberOfBadMuonCSCHits() const
Definition: HitPattern.h:962
reco::commonHits
PatternSet< N > commonHits(PatternSet< N > const &p1, PatternSet< N > const &p2)
Definition: HitPattern.h:1126
reco::HitPattern::hitCount
uint8_t hitCount
Definition: HitPattern.h:491
reco::HitPattern::numberOfBadMuonHits
int numberOfBadMuonHits() const
Definition: HitPattern.h:954
reco::HitPattern::pixelEndcapLayersWithoutMeasurement
int pixelEndcapLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:625
reco::HitPattern::numberOfValidMuonRPCHits
int numberOfValidMuonRPCHits() const
Definition: HitPattern.h:858
reco::HitPattern::numberOfLostTrackerHits
int numberOfLostTrackerHits(HitCategory category) const
Definition: HitPattern.h:880
reco::HitPattern::numberOfValidPixelHits
int numberOfValidPixelHits() const
Definition: HitPattern.h:818
reco::HitPattern::numberOfLostStripHits
int numberOfLostStripHits(HitCategory category) const
Definition: HitPattern.h:912
reco::HitPattern::MTD_HIT
Definition: HitPattern.h:151
test::TestHitPattern::test
int test()
reco::HitPattern::numberOfDTStationsWithRPhiView
int numberOfDTStationsWithRPhiView() const
Definition: HitPattern.cc:916
reco::HitPattern::NULL_RETURN
static const uint32_t NULL_RETURN
Definition: HitPattern.h:160
reco::HitPattern::endTrackHits
uint8_t endTrackHits
Definition: HitPattern.h:494
reco::HitPattern::SubstrOffset
const static unsigned short SubstrOffset
Definition: HitPattern.h:452
reco::HitPattern::numberOfLostMuonME0Hits
int numberOfLostMuonME0Hits() const
Definition: HitPattern.h:948
reco::HitPattern::MONO
Definition: HitPattern.h:149
reco::HitPattern::pixelBarrelHitFilter
static bool pixelBarrelHitFilter(uint16_t pattern)
Definition: HitPattern.h:584
reco::HitPattern::pixelLayersWithoutMeasurement
int pixelLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.h:1001
reco::HitPattern::numberOfValidStripLayersWithMonoAndStereo
int numberOfValidStripLayersWithMonoAndStereo() const
Definition: HitPattern.cc:390
reco::HitPattern::numberOfValidMuonME0Hits
int numberOfValidMuonME0Hits() const
Definition: HitPattern.h:866
reco::HitPattern::numberOfValidTIDLayersWithMonoAndStereo
int numberOfValidTIDLayersWithMonoAndStereo(uint32_t layer=0) const
Definition: HitPattern.cc:419
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:604
reco::HitPattern::numberOfValidTimingHits
int numberOfValidTimingHits() const
Definition: HitPattern.h:814
reco::HitPattern::outermostMuonStationWithHits
int outermostMuonStationWithHits(int hitType) const
hitType=-1(all), 0=valid, 3=bad; 0 = no stations at all
Definition: HitPattern.cc:902
reco::HitPattern::numberOfValidHits
int numberOfValidHits() const
Definition: HitPattern.h:804
reco::HitPattern::insertExpectedOuterHit
bool insertExpectedOuterHit(const uint16_t pattern)
Definition: HitPattern.cc:993
edm_modernize_messagelogger.stat
stat
Definition: edm_modernize_messagelogger.py:27
reco::HitPattern::stripTIBLayersTotallyOffOrBad
int stripTIBLayersTotallyOffOrBad(HitCategory category=TRACK_HITS) const
Definition: HitPattern.cc:698
SiStripSubdetector::Subdetector
Subdetector
Definition: SiStripEnums.h:5
reco::HitPattern::SubDetectorOffset
const static unsigned short SubDetectorOffset
Definition: HitPattern.h:456
reco::HitPattern::stripHitFilter
static bool stripHitFilter(uint16_t pattern)
Definition: HitPattern.h:602
reco::HitPattern::MaxHits
const static unsigned short MaxHits
Definition: HitPattern.h:158
reco::HitPattern::me0StationsWithValidHits
int me0StationsWithValidHits() const
Definition: HitPattern.h:1061
reco::HitPattern::gemStationsWithValidHits
int gemStationsWithValidHits() const
Definition: HitPattern.h:1055
reco::HitPattern::muonDTHitFilter
static bool muonDTHitFilter(uint16_t pattern)
Definition: HitPattern.h:630
reco::HitPattern::getTrackerMonoStereo
uint16_t getTrackerMonoStereo(HitCategory category, uint16_t substr, uint16_t layer) const
Definition: HitPattern.cc:455
reco::HitPattern::numberOfDTStationsWithRZView
int numberOfDTStationsWithRZView() const
Definition: HitPattern.cc:928
reco::HitPattern::getDTSuperLayer
static uint16_t getDTSuperLayer(uint16_t pattern)
DT superlayer (1-3). Where the "hit" was a DT segment, superlayer is 0. Only valid for muon DT patter...
Definition: HitPattern.h:751
reco::HitPattern::MISSING
Definition: HitPattern.h:153
reco::PatternSet
Definition: HitPattern.h:1080
StripSubdetector.h
reco::HitPattern::endInner
uint8_t endInner
Definition: HitPattern.h:496
reco::HitPattern::hitPattern
uint16_t hitPattern[ARRAY_LENGTH]
Definition: HitPattern.h:490
reco::HitPattern::numberOfLostMuonRPCHits
int numberOfLostMuonRPCHits() const
Definition: HitPattern.h:940
reco::HitPattern::numberOfValidTIBLayersWithMonoAndStereo
int numberOfValidTIBLayersWithMonoAndStereo(uint32_t layer=0) const
Definition: HitPattern.cc:415
reco::HitPattern::muonME0HitFilter
static bool muonME0HitFilter(uint16_t pattern)
Definition: HitPattern.h:666
StripSubdetector::TID
static constexpr auto TID
Definition: StripSubdetector.h:17
reco::HitPattern::stripTECLayersNull
int stripTECLayersNull() const
Definition: HitPattern.cc:790
hit
Definition: SiStripHitEffFromCalibTree.cc:88
reco::HitPattern::stripTECLayersWithoutMeasurement
int stripTECLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:666
reco::HitPattern::call
void call(HitCategory category, filterType typeFilter, F f) const
Definition: HitPattern.h:552
reco::PatternSet::MaxHits
static constexpr int MaxHits
Definition: HitPattern.h:1081
reco::HitPattern::insertTrackHit
bool insertTrackHit(const uint16_t pattern)
Definition: HitPattern.cc:963
reco::HitPattern::stripTOBLayersNull
int stripTOBLayersNull() const
Definition: HitPattern.cc:780
reco::HitPattern::stripTIDLayersNull
int stripTIDLayersNull() const
Definition: HitPattern.cc:770
reco::HitPattern::dtStationsWithAnyHits
int dtStationsWithAnyHits() const
Definition: HitPattern.h:1041
Generated for CMSSW Reference Manual by doxygen 1.8.16