CMS 3D CMS Logo
HITrackClusterRemover.cc
Go to the documentation of this file.
1 #include "FWCore/Framework/interface/Frameworkfwd.h"
2 #include "FWCore/Framework/interface/stream/EDProducer.h"
3 #include "FWCore/Framework/interface/Event.h"
4 #include "FWCore/Framework/interface/ESHandle.h"
5 #include "FWCore/ParameterSet/interface/ParameterSet.h"
6 #include "FWCore/Utilities/interface/InputTag.h"
7 
8 #include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
9 #include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
10 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
11 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit1D.h"
12 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
13 #include "DataFormats/TrackerRecHit2D/interface/ProjectedSiStripRecHit2D.h"
14 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
15 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
16 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
17 #include "DataFormats/Common/interface/Handle.h"
18 #include "DataFormats/SiStripDetId/interface/SiStripDetId.h"
19 #include "DataFormats/Common/interface/DetSetVector.h"
20 #include "DataFormats/Common/interface/ValueMap.h"
21 #include "DataFormats/Common/interface/DetSetVectorNew.h"
22 #include "DataFormats/Provenance/interface/ProductID.h"
23 #include "DataFormats/Common/interface/ContainerMask.h"
24 
25 #include "DataFormats/DetId/interface/DetId.h"
26 
27 #include "DataFormats/TrackReco/interface/TrackFwd.h"
28 #include "DataFormats/TrackReco/interface/Track.h"
29 #include "DataFormats/TrackerRecHit2D/interface/ClusterRemovalInfo.h"
30 
31 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
32 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
33 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
34 #include "Geometry/CommonDetUnit/interface/GeomDetType.h"
35 
36 #include "RecoTracker/TransientTrackingRecHit/interface/Traj2TrackHits.h"
37 
38 //
39 // class decleration
40 //
41 
42 class HITrackClusterRemover : public edm::stream::EDProducer<> {
43 public:
44  HITrackClusterRemover(const edm::ParameterSet &iConfig);
45  ~HITrackClusterRemover() override;
46  void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) override;
47 
48 private:
49  struct ParamBlock {
50  ParamBlock() : isSet_(false), usesCharge_(false) {}
51  ParamBlock(const edm::ParameterSet &iConfig)
52  : isSet_(true),
53  usesCharge_(iConfig.exists("maxCharge")),
54  usesSize_(iConfig.exists("maxSize")),
55  cutOnPixelCharge_(iConfig.exists("minGoodPixelCharge")),
56  cutOnStripCharge_(iConfig.exists("minGoodStripCharge")),
57  maxChi2_(iConfig.getParameter<double>("maxChi2")),
58  maxCharge_(usesCharge_ ? iConfig.getParameter<double>("maxCharge") : 0),
59  minGoodPixelCharge_(cutOnPixelCharge_ ? iConfig.getParameter<double>("minGoodPixelCharge") : 0),
60  minGoodStripCharge_(cutOnStripCharge_ ? iConfig.getParameter<double>("minGoodStripCharge") : 0),
61  maxSize_(usesSize_ ? iConfig.getParameter<uint32_t>("maxSize") : 0) {}
62  bool isSet_, usesCharge_, usesSize_, cutOnPixelCharge_, cutOnStripCharge_;
63  float maxChi2_, maxCharge_, minGoodPixelCharge_, minGoodStripCharge_;
64  size_t maxSize_;
65  };
66  static const unsigned int NumberOfParamBlocks = 6;
67 
68  bool doTracks_;
69  bool doStrip_, doPixel_;
70  bool mergeOld_;
71  typedef edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > PixelMaskContainer;
72  typedef edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > StripMaskContainer;
73  edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster> > pixelClusters_;
74  edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > stripClusters_;
75  edm::EDGetTokenT<reco::TrackCollection> tracks_;
76  edm::EDGetTokenT<reco::ClusterRemovalInfo> oldRemovalInfo_;
77  edm::EDGetTokenT<PixelMaskContainer> oldPxlMaskToken_;
78  edm::EDGetTokenT<StripMaskContainer> oldStrMaskToken_;
79  std::vector<edm::EDGetTokenT<edm::ValueMap<int> > > overrideTrkQuals_;
80  edm::EDGetTokenT<SiStripRecHit2DCollection> rphiRecHitToken_, stereoRecHitToken_;
81  // edm::EDGetTokenT<SiPixelRecHitCollection> pixelRecHitsToken_;
82 
83  ParamBlock pblocks_[NumberOfParamBlocks];
84  void readPSet(const edm::ParameterSet &iConfig,
85  const std::string &name,
86  int id1 = -1,
87  int id2 = -1,
88  int id3 = -1,
89  int id4 = -1,
90  int id5 = -1,
91  int id6 = -1);
92 
93  std::vector<uint8_t> pixels, strips; // avoid unneed alloc/dealloc of this
94  edm::ProductID pixelSourceProdID, stripSourceProdID; // ProdIDs refs must point to (for consistency tests)
95 
96  inline void process(const TrackingRecHit *hit, unsigned char chi2, const TrackerGeometry *tg);
97  inline void process(const OmniClusterRef &cluRef, SiStripDetId &detid, bool fromTrack);
98 
99  template <typename T>
100  std::unique_ptr<edmNew::DetSetVector<T> > cleanup(const edmNew::DetSetVector<T> &oldClusters,
101  const std::vector<uint8_t> &isGood,
102  reco::ClusterRemovalInfo::Indices &refs,
103  const reco::ClusterRemovalInfo::Indices *oldRefs);
104 
105  // Carries in full removal info about a given det from oldRefs
106  void mergeOld(reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices &oldRefs);
107 
108  bool clusterWasteSolution_, doStripChargeCheck_, doPixelChargeCheck_;
109  std::string stripRecHits_, pixelRecHits_;
110  bool filterTracks_;
111  int minNumberOfLayersWithMeasBeforeFiltering_;
112  reco::TrackBase::TrackQuality trackQuality_;
113  std::vector<bool> collectedStrips_;
114  std::vector<bool> collectedPixels_;
115 
116  float sensorThickness(const SiStripDetId &detid) const;
117 };
118 
119 using namespace std;
120 using namespace edm;
121 using namespace reco;
122 
123 void HITrackClusterRemover::readPSet(
124  const edm::ParameterSet &iConfig, const std::string &name, int id1, int id2, int id3, int id4, int id5, int id6) {
125  if (iConfig.exists(name)) {
126  ParamBlock pblock(iConfig.getParameter<ParameterSet>(name));
127  if (id1 == -1) {
128  fill(pblocks_, pblocks_ + NumberOfParamBlocks, pblock);
129  } else {
130  pblocks_[id1] = pblock;
131  if (id2 != -1)
132  pblocks_[id2] = pblock;
133  if (id3 != -1)
134  pblocks_[id3] = pblock;
135  if (id4 != -1)
136  pblocks_[id4] = pblock;
137  if (id5 != -1)
138  pblocks_[id5] = pblock;
139  if (id6 != -1)
140  pblocks_[id6] = pblock;
141  }
142  }
143 }
144 
145 HITrackClusterRemover::HITrackClusterRemover(const ParameterSet &iConfig)
146  : doTracks_(iConfig.exists("trajectories")),
147  doStrip_(iConfig.existsAs<bool>("doStrip") ? iConfig.getParameter<bool>("doStrip") : true),
148  doPixel_(iConfig.existsAs<bool>("doPixel") ? iConfig.getParameter<bool>("doPixel") : true),
149  mergeOld_(iConfig.exists("oldClusterRemovalInfo")),
150  clusterWasteSolution_(true),
151  doStripChargeCheck_(
152  iConfig.existsAs<bool>("doStripChargeCheck") ? iConfig.getParameter<bool>("doStripChargeCheck") : false),
153  doPixelChargeCheck_(
154  iConfig.existsAs<bool>("doPixelChargeCheck") ? iConfig.getParameter<bool>("doPixelChargeCheck") : false),
155  stripRecHits_(doStripChargeCheck_ ? iConfig.getParameter<std::string>("stripRecHits")
156  : std::string("siStripMatchedRecHits")),
157  pixelRecHits_(doPixelChargeCheck_ ? iConfig.getParameter<std::string>("pixelRecHits")
158  : std::string("siPixelRecHits")) {
159  mergeOld_ = mergeOld_ && !iConfig.getParameter<InputTag>("oldClusterRemovalInfo").label().empty();
160  if (iConfig.exists("overrideTrkQuals"))
161  overrideTrkQuals_.push_back(consumes<edm::ValueMap<int> >(iConfig.getParameter<InputTag>("overrideTrkQuals")));
162  if (iConfig.exists("clusterLessSolution"))
163  clusterWasteSolution_ = !iConfig.getParameter<bool>("clusterLessSolution");
164  if ((doPixelChargeCheck_ && !doPixel_) || (doStripChargeCheck_ && !doStrip_))
165  throw cms::Exception("Configuration Error")
166  << "HITrackClusterRemover: Charge check asked without cluster collection ";
167  if (doPixelChargeCheck_)
168  throw cms::Exception("Configuration Error")
169  << "HITrackClusterRemover: Pixel cluster charge check not yet implemented";
170 
171  if (doPixel_ && clusterWasteSolution_)
172  produces<edmNew::DetSetVector<SiPixelCluster> >();
173  if (doStrip_ && clusterWasteSolution_)
174  produces<edmNew::DetSetVector<SiStripCluster> >();
175  if (clusterWasteSolution_)
176  produces<ClusterRemovalInfo>();
177 
178  assert(!clusterWasteSolution_);
179 
180  fill(pblocks_, pblocks_ + NumberOfParamBlocks, ParamBlock());
181  readPSet(iConfig, "Common", -1);
182  if (doPixel_) {
183  readPSet(iConfig, "Pixel", 0, 1);
184  readPSet(iConfig, "PXB", 0);
185  readPSet(iConfig, "PXE", 1);
186  }
187  if (doStrip_) {
188  readPSet(iConfig, "Strip", 2, 3, 4, 5);
189  readPSet(iConfig, "StripInner", 2, 3);
190  readPSet(iConfig, "StripOuter", 4, 5);
191  readPSet(iConfig, "TIB", 2);
192  readPSet(iConfig, "TID", 3);
193  readPSet(iConfig, "TOB", 4);
194  readPSet(iConfig, "TEC", 5);
195  }
196 
197  bool usingCharge = false;
198  for (size_t i = 0; i < NumberOfParamBlocks; ++i) {
199  if (!pblocks_[i].isSet_)
200  throw cms::Exception("Configuration Error")
201  << "HITrackClusterRemover: Missing configuration for detector with subDetID = " << (i + 1);
202  if (pblocks_[i].usesCharge_ && !usingCharge) {
203  throw cms::Exception("Configuration Error") << "HITrackClusterRemover: Configuration for subDetID = " << (i + 1)
204  << " uses cluster charge, which is not enabled.";
205  }
206  }
207 
208  if (!clusterWasteSolution_) {
209  produces<edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > >();
210  produces<edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > >();
211  }
212  trackQuality_ = reco::TrackBase::undefQuality;
213  filterTracks_ = false;
214  if (iConfig.exists("TrackQuality")) {
215  filterTracks_ = true;
216  trackQuality_ = reco::TrackBase::qualityByName(iConfig.getParameter<std::string>("TrackQuality"));
217  minNumberOfLayersWithMeasBeforeFiltering_ =
218  iConfig.existsAs<int>("minNumberOfLayersWithMeasBeforeFiltering")
219  ? iConfig.getParameter<int>("minNumberOfLayersWithMeasBeforeFiltering")
220  : 0;
221  }
222 
223  if (doTracks_)
224  tracks_ = consumes<reco::TrackCollection>(iConfig.getParameter<InputTag>("trajectories"));
225  if (doPixel_)
226  pixelClusters_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(iConfig.getParameter<InputTag>("pixelClusters"));
227  if (doStrip_)
228  stripClusters_ = consumes<edmNew::DetSetVector<SiStripCluster> >(iConfig.getParameter<InputTag>("stripClusters"));
229  if (mergeOld_) {
230  oldRemovalInfo_ = consumes<ClusterRemovalInfo>(iConfig.getParameter<InputTag>("oldClusterRemovalInfo"));
231  oldPxlMaskToken_ = consumes<PixelMaskContainer>(iConfig.getParameter<InputTag>("oldClusterRemovalInfo"));
232  oldStrMaskToken_ = consumes<StripMaskContainer>(iConfig.getParameter<InputTag>("oldClusterRemovalInfo"));
233  }
234 
235  if (doStripChargeCheck_) {
236  rphiRecHitToken_ = consumes<SiStripRecHit2DCollection>(InputTag(stripRecHits_, "rphiRecHit"));
237  stereoRecHitToken_ = consumes<SiStripRecHit2DCollection>(InputTag(stripRecHits_, "stereoRecHit"));
238  }
239  // if(doPixelChargeCheck_) pixelRecHitsToken_ = consumes<SiPixelRecHitCollection>(InputTag(pixelRecHits_));
240 }
241 
242 HITrackClusterRemover::~HITrackClusterRemover() {}
243 
244 void HITrackClusterRemover::mergeOld(ClusterRemovalInfo::Indices &refs, const ClusterRemovalInfo::Indices &oldRefs) {
245  for (size_t i = 0, n = refs.size(); i < n; ++i) {
246  refs[i] = oldRefs[refs[i]];
247  }
248 }
249 
250 template <typename T>
251 std::unique_ptr<edmNew::DetSetVector<T> > HITrackClusterRemover::cleanup(
252  const edmNew::DetSetVector<T> &oldClusters,
253  const std::vector<uint8_t> &isGood,
254  reco::ClusterRemovalInfo::Indices &refs,
255  const reco::ClusterRemovalInfo::Indices *oldRefs) {
256  typedef typename edmNew::DetSetVector<T> DSV;
257  typedef typename edmNew::DetSetVector<T>::FastFiller DSF;
258  typedef typename edmNew::DetSet<T> DS;
259  auto output = std::make_unique<DSV>();
260  output->reserve(oldClusters.size(), oldClusters.dataSize());
261 
262  unsigned int countOld = 0;
263  unsigned int countNew = 0;
264 
265  // cluster removal loop
266  const T *firstOffset = &oldClusters.data().front();
267  for (typename DSV::const_iterator itdet = oldClusters.begin(), enddet = oldClusters.end(); itdet != enddet; ++itdet) {
268  DS oldDS = *itdet;
269 
270  if (oldDS.empty())
271  continue; // skip empty detsets
272 
273  uint32_t id = oldDS.detId();
274  DSF outds(*output, id);
275 
276  for (typename DS::const_iterator it = oldDS.begin(), ed = oldDS.end(); it != ed; ++it) {
277  uint32_t index = ((&*it) - firstOffset);
278  countOld++;
279  if (isGood[index]) {
280  outds.push_back(*it);
281  countNew++;
282  refs.push_back(index);
283  //std::cout << "HITrackClusterRemover::cleanup " << typeid(T).name() << " reference " << index << " to " << (refs.size() - 1) << std::endl;
284  }
285  }
286  if (outds.empty())
287  outds.abort(); // not write in an empty DSV
288  }
289 
290  if (oldRefs != nullptr)
291  mergeOld(refs, *oldRefs);
292  return output;
293 }
294 
295 float HITrackClusterRemover::sensorThickness(const SiStripDetId &detid) const {
296  if (detid.subdetId() >= SiStripDetId::TIB) {
297  if (detid.subdetId() == SiStripDetId::TOB)
298  return 0.047;
299  if (detid.moduleGeometry() == SiStripModuleGeometry::W5 || detid.moduleGeometry() == SiStripModuleGeometry::W6 ||
300  detid.moduleGeometry() == SiStripModuleGeometry::W7)
301  return 0.047;
302  return 0.029; // so it is TEC ring 1-4 or TIB or TOB;
303  } else if (detid.subdetId() == PixelSubdetector::PixelBarrel)
304  return 0.0285;
305  else
306  return 0.027;
307 }
308 
309 void HITrackClusterRemover::process(OmniClusterRef const &ocluster, SiStripDetId &detid, bool fromTrack) {
310  SiStripRecHit2D::ClusterRef cluster = ocluster.cluster_strip();
311  if (cluster.id() != stripSourceProdID)
312  throw cms::Exception("Inconsistent Data")
313  << "HITrackClusterRemover: strip cluster ref from Product ID = " << cluster.id()
314  << " does not match with source cluster collection (ID = " << stripSourceProdID << ")\n.";
315 
316  uint32_t subdet = detid.subdetId();
317  assert(cluster.id() == stripSourceProdID);
318  if (pblocks_[subdet - 1].usesSize_ && (cluster->amplitudes().size() > pblocks_[subdet - 1].maxSize_))
319  return;
320  if (!fromTrack) {
321  int clusCharge = 0;
322  for (std::vector<uint8_t>::const_iterator iAmp = cluster->amplitudes().begin(); iAmp != cluster->amplitudes().end();
323  ++iAmp) {
324  clusCharge += *iAmp;
325  }
326  if (pblocks_[subdet - 1].cutOnStripCharge_ &&
327  (clusCharge > (pblocks_[subdet - 1].minGoodStripCharge_ * sensorThickness(detid))))
328  return;
329  }
330 
331  if (collectedStrips_.size() <= cluster.key())
332  edm::LogError("BadCollectionSize") << collectedStrips_.size() << " is smaller than " << cluster.key();
333 
334  assert(collectedStrips_.size() > cluster.key());
335  strips[cluster.key()] = false;
336  if (!clusterWasteSolution_)
337  collectedStrips_[cluster.key()] = true;
338 }
339 
340 void HITrackClusterRemover::process(const TrackingRecHit *hit, unsigned char chi2, const TrackerGeometry *tg) {
341  SiStripDetId detid = hit->geographicalId();
342  uint32_t subdet = detid.subdetId();
343 
344  assert((subdet > 0) && (subdet <= NumberOfParamBlocks));
345 
346  // chi2 cut
347  if (chi2 > Traj2TrackHits::toChi2x5(pblocks_[subdet - 1].maxChi2_))
348  return;
349 
350  if (GeomDetEnumerators::isTrackerPixel(tg->geomDetSubDetector(subdet))) {
351  if (!doPixel_)
352  return;
353  // this is a pixel, and i *know* it is
354  const SiPixelRecHit *pixelHit = static_cast<const SiPixelRecHit *>(hit);
355 
356  SiPixelRecHit::ClusterRef cluster = pixelHit->cluster();
357 
358  if (cluster.id() != pixelSourceProdID)
359  throw cms::Exception("Inconsistent Data")
360  << "HITrackClusterRemover: pixel cluster ref from Product ID = " << cluster.id()
361  << " does not match with source cluster collection (ID = " << pixelSourceProdID << ")\n.";
362 
363  assert(cluster.id() == pixelSourceProdID);
364  //DBG// cout << "HIT NEW PIXEL DETID = " << detid.rawId() << ", Cluster [ " << cluster.key().first << " / " << cluster.key().second << " ] " << endl;
365 
366  // if requested, cut on cluster size
367  if (pblocks_[subdet - 1].usesSize_ && (cluster->pixels().size() > pblocks_[subdet - 1].maxSize_))
368  return;
369 
370  // mark as used
371  pixels[cluster.key()] = false;
372 
373  //if(!clusterWasteSolution_) collectedPixel[detid.rawId()].insert(cluster);
374  assert(collectedPixels_.size() > cluster.key());
375  //assert(detid.rawId() == cluster->geographicalId()); //This condition fails
376  if (!clusterWasteSolution_)
377  collectedPixels_[cluster.key()] = true;
378 
379  } else { // aka Strip
380  if (!doStrip_)
381  return;
382  const type_info &hitType = typeid(*hit);
383  if (hitType == typeid(SiStripRecHit2D)) {
384  const SiStripRecHit2D *stripHit = static_cast<const SiStripRecHit2D *>(hit);
385  //DBG// cout << "Plain RecHit 2D: " << endl;
386  process(stripHit->omniClusterRef(), detid, true);
387  } else if (hitType == typeid(SiStripRecHit1D)) {
388  const SiStripRecHit1D *hit1D = static_cast<const SiStripRecHit1D *>(hit);
389  process(hit1D->omniClusterRef(), detid, true);
390  } else if (hitType == typeid(SiStripMatchedRecHit2D)) {
391  const SiStripMatchedRecHit2D *matchHit = static_cast<const SiStripMatchedRecHit2D *>(hit);
392  //DBG// cout << "Matched RecHit 2D: " << endl;
393  process(matchHit->monoClusterRef(), detid, true);
394  process(matchHit->stereoClusterRef(), detid, true);
395  } else if (hitType == typeid(ProjectedSiStripRecHit2D)) {
396  const ProjectedSiStripRecHit2D *projHit = static_cast<const ProjectedSiStripRecHit2D *>(hit);
397  //DBG// cout << "Projected RecHit 2D: " << endl;
398  process(projHit->originalHit().omniClusterRef(), detid, true);
399  } else
400  throw cms::Exception("NOT IMPLEMENTED")
401  << "Don't know how to handle " << hitType.name() << " on detid " << detid.rawId() << "\n";
402  }
403 }
404 
405 /* Schematic picture of n-th step Iterative removal
406  * (that os removing clusters after n-th step tracking)
407  * clusters: [ C1 ] -> [ C2 ] -> ... -> [ Cn ] -> [ Cn + 1 ]
408  * ^ ^ ^--- OUTPUT "new" ID
409  * |-- before any removal |----- Source clusters
410  * |-- OUTPUT "old" ID |----- Hits in Traj. point here
411  * | \----- Old ClusterRemovalInfo "new" ID
412  * \-- Old ClusterRemovalInfo "old" ID
413  */
414 
415 void HITrackClusterRemover::produce(Event &iEvent, const EventSetup &iSetup) {
416  ProductID pixelOldProdID, stripOldProdID;
417 
418  edm::ESHandle<TrackerGeometry> tgh;
419  iSetup.get<TrackerDigiGeometryRecord>().get(tgh);
420 
421  Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusters;
422  if (doPixel_) {
423  iEvent.getByToken(pixelClusters_, pixelClusters);
424  pixelSourceProdID = pixelClusters.id();
425  }
426  //DBG// std::cout << "HITrackClusterRemover: Read pixel " << pixelClusters_.encode() << " = ID " << pixelSourceProdID << std::endl;
427 
428  Handle<edmNew::DetSetVector<SiStripCluster> > stripClusters;
429  if (doStrip_) {
430  iEvent.getByToken(stripClusters_, stripClusters);
431  stripSourceProdID = stripClusters.id();
432  }
433  //DBG// std::cout << "HITrackClusterRemover: Read strip " << stripClusters_.encode() << " = ID " << stripSourceProdID << std::endl;
434 
435  std::unique_ptr<ClusterRemovalInfo> cri;
436  if (clusterWasteSolution_) {
437  if (doStrip_ && doPixel_)
438  cri = std::make_unique<ClusterRemovalInfo>(pixelClusters, stripClusters);
439  else if (doStrip_)
440  cri = std::make_unique<ClusterRemovalInfo>(stripClusters);
441  else if (doPixel_)
442  cri = std::make_unique<ClusterRemovalInfo>(pixelClusters);
443  }
444 
445  Handle<ClusterRemovalInfo> oldRemovalInfo;
446  if (mergeOld_ && clusterWasteSolution_) {
447  iEvent.getByToken(oldRemovalInfo_, oldRemovalInfo);
448  // Check ProductIDs
449  if ((oldRemovalInfo->stripNewRefProd().id() == stripClusters.id()) &&
450  (oldRemovalInfo->pixelNewRefProd().id() == pixelClusters.id())) {
451  cri->getOldClustersFrom(*oldRemovalInfo);
452 
453  pixelOldProdID = oldRemovalInfo->pixelRefProd().id();
454  stripOldProdID = oldRemovalInfo->stripRefProd().id();
455 
456  } else {
457  edm::EDConsumerBase::Labels labels;
458  labelsForToken(oldRemovalInfo_, labels);
459  throw cms::Exception("Inconsistent Data")
460  << "HITrackClusterRemover: "
461  << "Input collection product IDs are [pixel: " << pixelClusters.id() << ", strip: " << stripClusters.id()
462  << "] \n"
463  << "\t but the *old* ClusterRemovalInfo " << labels.productInstance << " refers as 'new product ids' to "
464  << "[pixel: " << oldRemovalInfo->pixelNewRefProd().id()
465  << ", strip: " << oldRemovalInfo->stripNewRefProd().id() << "]\n"
466  << "NOTA BENE: when running HITrackClusterRemover with an old ClusterRemovalInfo the hits in the trajectory "
467  "MUST be already re-keyed.\n";
468  }
469  } else { // then Old == Source
470  pixelOldProdID = pixelSourceProdID;
471  stripOldProdID = stripSourceProdID;
472  }
473 
474  if (doStrip_) {
475  strips.resize(stripClusters->dataSize());
476  fill(strips.begin(), strips.end(), true);
477  }
478  if (doPixel_) {
479  pixels.resize(pixelClusters->dataSize());
480  fill(pixels.begin(), pixels.end(), true);
481  }
482  if (mergeOld_) {
483  edm::Handle<PixelMaskContainer> oldPxlMask;
484  edm::Handle<StripMaskContainer> oldStrMask;
485  iEvent.getByToken(oldPxlMaskToken_, oldPxlMask);
486  iEvent.getByToken(oldStrMaskToken_, oldStrMask);
487  LogDebug("HITrackClusterRemover") << "to merge in, " << oldStrMask->size() << " strp and " << oldPxlMask->size()
488  << " pxl";
489  oldStrMask->copyMaskTo(collectedStrips_);
490  oldPxlMask->copyMaskTo(collectedPixels_);
491  assert(stripClusters->dataSize() >= collectedStrips_.size());
492  collectedStrips_.resize(stripClusters->dataSize(), false);
493  } else {
494  collectedStrips_.resize(stripClusters->dataSize(), false);
495  collectedPixels_.resize(pixelClusters->dataSize(), false);
496  }
497 
498  if (doTracks_) {
499  Handle<reco::TrackCollection> tracks;
500  iEvent.getByToken(tracks_, tracks);
501 
502  std::vector<Handle<edm::ValueMap<int> > > quals;
503  if (!overrideTrkQuals_.empty()) {
504  quals.resize(1);
505  iEvent.getByToken(overrideTrkQuals_[0], quals[0]);
506  }
507  int it = 0;
508  for (const auto &track : *tracks) {
509  if (filterTracks_) {
510  bool goodTk = true;
511  if (!quals.empty()) {
512  int qual = (*(quals[0])).get(it++);
513  if (qual < 0) {
514  goodTk = false;
515  }
516  //note that this does not work for some trackquals (goodIterative or undefQuality)
517  else
518  goodTk = (qual & (1 << trackQuality_)) >> trackQuality_;
519  } else
520  goodTk = (track.quality(trackQuality_));
521  if (!goodTk)
522  continue;
523  if (track.hitPattern().trackerLayersWithMeasurement() < minNumberOfLayersWithMeasBeforeFiltering_)
524  continue;
525  }
526  auto const &chi2sX5 = track.extra()->chi2sX5();
527  assert(chi2sX5.size() == track.recHitsSize());
528  auto hb = track.recHitsBegin();
529  for (unsigned int h = 0; h < track.recHitsSize(); h++) {
530  auto hit = *(hb + h);
531  if (!hit->isValid())
532  continue;
533  process(hit, chi2sX5[h], tgh.product());
534  }
535  }
536  }
537 
538  if (doStripChargeCheck_) {
539  edm::Handle<SiStripRecHit2DCollection> rechitsrphi;
540  iEvent.getByToken(rphiRecHitToken_, rechitsrphi);
541  const SiStripRecHit2DCollection::DataContainer *rphiRecHits = &(rechitsrphi).product()->data();
542  for (SiStripRecHit2DCollection::DataContainer::const_iterator recHit = rphiRecHits->begin();
543  recHit != rphiRecHits->end();
544  recHit++) {
545  SiStripDetId detid = recHit->geographicalId();
546  process(recHit->omniClusterRef(), detid, false);
547  }
548  edm::Handle<SiStripRecHit2DCollection> rechitsstereo;
549  iEvent.getByToken(stereoRecHitToken_, rechitsstereo);
550  const SiStripRecHit2DCollection::DataContainer *stereoRecHits = &(rechitsstereo).product()->data();
551  for (SiStripRecHit2DCollection::DataContainer::const_iterator recHit = stereoRecHits->begin();
552  recHit != stereoRecHits->end();
553  recHit++) {
554  SiStripDetId detid = recHit->geographicalId();
555  process(recHit->omniClusterRef(), detid, false);
556  }
557  }
558  // if(doPixelChargeCheck_) {
559  // edm::Handle<SiPixelRecHitCollection> pixelrechits;
560  // iEvent.getByToken(pixelRecHitsToken_,pixelrechits);
561  // }
562 
563  if (doPixel_ && clusterWasteSolution_) {
564  OrphanHandle<edmNew::DetSetVector<SiPixelCluster> > newPixels = iEvent.put(
565  cleanup(*pixelClusters, pixels, cri->pixelIndices(), mergeOld_ ? &oldRemovalInfo->pixelIndices() : nullptr));
566  //DBG// std::cout << "HITrackClusterRemover: Wrote pixel " << newPixels.id() << " from " << pixelSourceProdID << std::endl;
567  cri->setNewPixelClusters(newPixels);
568  }
569  if (doStrip_ && clusterWasteSolution_) {
570  OrphanHandle<edmNew::DetSetVector<SiStripCluster> > newStrips = iEvent.put(
571  cleanup(*stripClusters, strips, cri->stripIndices(), mergeOld_ ? &oldRemovalInfo->stripIndices() : nullptr));
572  //DBG// std::cout << "HITrackClusterRemover: Wrote strip " << newStrips.id() << " from " << stripSourceProdID << std::endl;
573  cri->setNewStripClusters(newStrips);
574  }
575 
576  if (clusterWasteSolution_) {
577  // double fraction_pxl= cri->pixelIndices().size() / (double) pixels.size();
578  // double fraction_strp= cri->stripIndices().size() / (double) strips.size();
579  // edm::LogWarning("HITrackClusterRemover")<<" fraction: " << fraction_pxl <<" "<<fraction_strp;
580  iEvent.put(std::move(cri));
581  }
582 
583  pixels.clear();
584  strips.clear();
585 
586  if (!clusterWasteSolution_) {
587  //auto_ptr<edmNew::DetSetVector<SiPixelClusterRefNew> > removedPixelClsuterRefs(new edmNew::DetSetVector<SiPixelClusterRefNew>());
588  //auto_ptr<edmNew::DetSetVector<SiStripRecHit1D::ClusterRef> > removedStripClsuterRefs(new edmNew::DetSetVector<SiStripRecHit1D::ClusterRef>());
589 
590  LogDebug("HITrackClusterRemover") << "total strip to skip: "
591  << std::count(collectedStrips_.begin(), collectedStrips_.end(), true);
592  // std::cout << "HITrackClusterRemover " <<"total strip to skip: "<<std::count(collectedStrips_.begin(),collectedStrips_.end(),true) <<std::endl;
593  iEvent.put(std::make_unique<StripMaskContainer>(edm::RefProd<edmNew::DetSetVector<SiStripCluster> >(stripClusters),
594  collectedStrips_));
595 
596  LogDebug("HITrackClusterRemover") << "total pxl to skip: "
597  << std::count(collectedPixels_.begin(), collectedPixels_.end(), true);
598  iEvent.put(std::make_unique<PixelMaskContainer>(edm::RefProd<edmNew::DetSetVector<SiPixelCluster> >(pixelClusters),
599  collectedPixels_));
600  }
601  collectedStrips_.clear();
602  collectedPixels_.clear();
603 }
604 
605 #include "FWCore/PluginManager/interface/ModuleDef.h"
606 #include "FWCore/Framework/interface/MakerMacros.h"
607 DEFINE_FWK_MODULE(HITrackClusterRemover);
edm::RefProd
Definition: EDProductfwd.h:25
ClusterRemovalInfo.h
SummaryClient_cfi.labels
labels
Definition: SummaryClient_cfi.py:61
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
Handle.h
PDWG_EXOHSCP_cff.tracks
tracks
Definition: PDWG_EXOHSCP_cff.py:28
electrons_cff.bool
bool
Definition: electrons_cff.py:372
SiStripDetId::moduleGeometry
SiStripModuleGeometry moduleGeometry() const
Definition: SiStripDetId.h:109
mps_fire.i
i
Definition: mps_fire.py:355
HITrackClusterRemover::ParamBlock::cutOnPixelCharge_
bool cutOnPixelCharge_
Definition: HITrackClusterRemover.cc:62
HITrackClusterRemover::ParamBlock::maxChi2_
float maxChi2_
Definition: HITrackClusterRemover.cc:63
HITrackClusterRemover::doPixelChargeCheck_
bool doPixelChargeCheck_
Definition: HITrackClusterRemover.cc:108
HITrackClusterRemover::pixelClusters_
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > > pixelClusters_
Definition: HITrackClusterRemover.cc:73
SiStripModuleGeometry::W7
LumiMonitor_cff.pixelClusters
pixelClusters
Definition: LumiMonitor_cff.py:10
HITrackClusterRemover
Definition: HITrackClusterRemover.cc:42
funct::false
false
Definition: Factorize.h:34
ContainerMask.h
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
edm::RefProd::id
ProductID id() const
Accessor for product ID.
Definition: RefProd.h:124
TrackerGeometry.h
PixelSubdetector::PixelBarrel
Definition: PixelSubdetector.h:11
reco::ClusterRemovalInfo::Indices
std::vector< uint32_t > Indices
Definition: ClusterRemovalInfo.h:22
ESHandle.h
HITrackClusterRemover::oldPxlMaskToken_
edm::EDGetTokenT< PixelMaskContainer > oldPxlMaskToken_
Definition: HITrackClusterRemover.cc:77
HITrackClusterRemover::ParamBlock::ParamBlock
ParamBlock(const edm::ParameterSet &iConfig)
Definition: HITrackClusterRemover.cc:51
Traj2TrackHits::toChi2x5
static unsigned char toChi2x5(float chi2)
Definition: Traj2TrackHits.h:25
convertSQLitetoXML_cfg.output
output
Definition: convertSQLitetoXML_cfg.py:32
edmNew::DetSetVector::size
size_type size() const
Definition: DetSetVectorNew.h:596
edm::EDGetTokenT
Definition: EDGetToken.h:33
edm
HLT enums.
Definition: AlignableModifier.h:19
HITrackClusterRemover::pblocks_
ParamBlock pblocks_[NumberOfParamBlocks]
Definition: HITrackClusterRemover.cc:83
HITrackClusterRemover::ParamBlock::maxCharge_
float maxCharge_
Definition: HITrackClusterRemover.cc:63
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
reco::TrackBase::undefQuality
Definition: TrackBase.h:151
globals_cff.id1
id1
Definition: globals_cff.py:32
reco::ClusterRemovalInfo::stripRefProd
const SiStripClusterRefProd & stripRefProd() const
Definition: ClusterRemovalInfo.h:58
SiPixelCluster.h
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:231
GeomDetType.h
HITrackClusterRemover::process
void process(const TrackingRecHit *hit, unsigned char chi2, const TrackerGeometry *tg)
Definition: HITrackClusterRemover.cc:340
cms::cuda::assert
assert(be >=bs)
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
EDProducer.h
HITrackClusterRemover::StripMaskContainer
edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > > StripMaskContainer
Definition: HITrackClusterRemover.cc:72
SiStripDetId.h
edmNew::DetSetVector::begin
const_iterator begin(bool update=false) const
Definition: DetSetVectorNew.h:572
SiStripRecHit2D
Definition: SiStripRecHit2D.h:7
edm::ContainerMask
Definition: ContainerMask.h:36
HITrackClusterRemover::trackQuality_
reco::TrackBase::TrackQuality trackQuality_
Definition: HITrackClusterRemover.cc:112
edm::ParameterSet::existsAs
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:160
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
hltPixelTracks_cff.chi2
chi2
Definition: hltPixelTracks_cff.py:25
edm::Handle
Definition: AssociativeIterator.h:50
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:47
OmniClusterRef
Definition: OmniClusterRef.h:12
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
ProjectedSiStripRecHit2D
Definition: ProjectedSiStripRecHit2D.h:8
HITrackClusterRemover::pixelRecHits_
std::string pixelRecHits_
Definition: HITrackClusterRemover.cc:109
HITrackClusterRemover::ParamBlock::usesCharge_
bool usesCharge_
Definition: HITrackClusterRemover.cc:62
SiStripModuleGeometry::W6
HITrackClusterRemover::strips
std::vector< uint8_t > strips
Definition: HITrackClusterRemover.cc:93
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
HITrackClusterRemover::pixelSourceProdID
edm::ProductID pixelSourceProdID
Definition: HITrackClusterRemover.cc:94
edm::Ref
Definition: AssociativeIterator.h:58
reco::ClusterRemovalInfo::pixelIndices
Indices & pixelIndices()
Definition: ClusterRemovalInfo.h:49
HITrackClusterRemover::readPSet
void readPSet(const edm::ParameterSet &iConfig, const std::string &name, int id1=-1, int id2=-1, int id3=-1, int id4=-1, int id5=-1, int id6=-1)
Definition: HITrackClusterRemover.cc:123
HITrackClusterRemover::stripSourceProdID
edm::ProductID stripSourceProdID
Definition: HITrackClusterRemover.cc:94
HITrackClusterRemover::~HITrackClusterRemover
~HITrackClusterRemover() override
Definition: HITrackClusterRemover.cc:242
MakerMacros.h
h
OmniClusterRef::cluster_strip
ClusterStripRef cluster_strip() const
Definition: OmniClusterRef.h:44
Track.h
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
TrackFwd.h
SiStripModuleGeometry::W5
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
HITrackClusterRemover::ParamBlock::usesSize_
bool usesSize_
Definition: HITrackClusterRemover.cc:62
HITrackClusterRemover::collectedPixels_
std::vector< bool > collectedPixels_
Definition: HITrackClusterRemover.cc:114
HLT_2018_cff.stripClusters
stripClusters
Definition: HLT_2018_cff.py:8683
HITrackClusterRemover::ParamBlock::isSet_
bool isSet_
Definition: HITrackClusterRemover.cc:62
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
HITrackClusterRemover::collectedStrips_
std::vector< bool > collectedStrips_
Definition: HITrackClusterRemover.cc:113
HITrackClusterRemover::sensorThickness
float sensorThickness(const SiStripDetId &detid) const
Definition: HITrackClusterRemover.cc:295
ProductID.h
edm::ESHandle< TrackerGeometry >
edmNew::DetSet
Definition: DetSetNew.h:22
ProjectedSiStripRecHit2D::originalHit
SiStripRecHit2D originalHit() const
Definition: ProjectedSiStripRecHit2D.h:56
HITrackClusterRemover::minNumberOfLayersWithMeasBeforeFiltering_
int minNumberOfLayersWithMeasBeforeFiltering_
Definition: HITrackClusterRemover.cc:111
HITrackClusterRemover::HITrackClusterRemover
HITrackClusterRemover(const edm::ParameterSet &iConfig)
Definition: HITrackClusterRemover.cc:145
HITrackClusterRemover::filterTracks_
bool filterTracks_
Definition: HITrackClusterRemover.cc:110
HITrackClusterRemover::rphiRecHitToken_
edm::EDGetTokenT< SiStripRecHit2DCollection > rphiRecHitToken_
Definition: HITrackClusterRemover.cc:80
HITrackClusterRemover::stripClusters_
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > stripClusters_
Definition: HITrackClusterRemover.cc:74
Traj2TrackHits.h
SiPixelRecHit.h
edmNew::DetSetVector::dataSize
size_type dataSize() const
Definition: DetSetVectorNew.h:594
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
HITrackClusterRemover::NumberOfParamBlocks
static const unsigned int NumberOfParamBlocks
Definition: HITrackClusterRemover.cc:66
LaserSeedGenerator_cfi.stereoRecHits
stereoRecHits
Definition: LaserSeedGenerator_cfi.py:13
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HITrackClusterRemover::produce
void produce(edm::Event &iEvent, const edm::EventSetup &iSetup) override
Definition: HITrackClusterRemover.cc:415
HITrackClusterRemover::tracks_
edm::EDGetTokenT< reco::TrackCollection > tracks_
Definition: HITrackClusterRemover.cc:75
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:674
funct::true
true
Definition: Factorize.h:173
ntuplemaker.fill
fill
Definition: ntuplemaker.py:304
TrackerDigiGeometryRecord.h
SiPixelRecHitCollection.h
HITrackClusterRemover::ParamBlock::cutOnStripCharge_
bool cutOnStripCharge_
Definition: HITrackClusterRemover.cc:62
HITrackClusterRemover::overrideTrkQuals_
std::vector< edm::EDGetTokenT< edm::ValueMap< int > > > overrideTrkQuals_
Definition: HITrackClusterRemover.cc:79
HITrackClusterRemover::doTracks_
bool doTracks_
Definition: HITrackClusterRemover.cc:68
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:79016
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:670
TrackerGeometry::geomDetSubDetector
const GeomDetEnumerators::SubDetector geomDetSubDetector(int subdet) const
Definition: TrackerGeometry.cc:203
edm::ParameterSet
Definition: ParameterSet.h:36
edm::LogError
Definition: MessageLogger.h:183
SiStripCluster.h
Event.h
HITrackClusterRemover::doStripChargeCheck_
bool doStripChargeCheck_
Definition: HITrackClusterRemover.cc:108
SiStripMatchedRecHit2D::stereoClusterRef
OmniClusterRef const & stereoClusterRef() const
Definition: SiStripMatchedRecHit2D.h:34
HITrackClusterRemover::ParamBlock::minGoodStripCharge_
float minGoodStripCharge_
Definition: HITrackClusterRemover.cc:63
KineDebug3::count
void count()
Definition: KinematicConstrainedVertexUpdatorT.h:21
HITrackClusterRemover::oldRemovalInfo_
edm::EDGetTokenT< reco::ClusterRemovalInfo > oldRemovalInfo_
Definition: HITrackClusterRemover.cc:76
HITrackClusterRemover::stripRecHits_
std::string stripRecHits_
Definition: HITrackClusterRemover.cc:109
ModuleDef.h
GeomDetEnumerators::isTrackerPixel
bool isTrackerPixel(GeomDetEnumerators::SubDetector m)
Definition: GeomDetEnumerators.cc:68
iEvent
int iEvent
Definition: GenABIO.cc:224
SiStripRecHit2DCollection.h
edm::stream::EDProducer
Definition: EDProducer.h:38
SiStripDetId::TOB
static constexpr auto TOB
Definition: SiStripDetId.h:39
ProjectedSiStripRecHit2D.h
HITrackClusterRemover::ParamBlock::minGoodPixelCharge_
float minGoodPixelCharge_
Definition: HITrackClusterRemover.cc:63
edm::EventSetup
Definition: EventSetup.h:57
edmNew::DetSetVector::DataContainer
std::vector< data_type > DataContainer
Definition: DetSetVectorNew.h:204
DetSetVector.h
SiStripRecHit1D
Definition: SiStripRecHit1D.h:8
reco::ClusterRemovalInfo::stripNewRefProd
const SiStripClusterRefProd & stripNewRefProd() const
Definition: ClusterRemovalInfo.h:61
get
#define get
SiStripMatchedRecHit2D::monoClusterRef
OmniClusterRef const & monoClusterRef() const
Definition: SiStripMatchedRecHit2D.h:35
hcalSimParameters_cfi.hb
hb
Definition: hcalSimParameters_cfi.py:57
edm::ContainerMask::size
size_t size() const
Definition: ContainerMask.h:53
InputTag.h
TrackingRecHit
Definition: TrackingRecHit.h:21
SiStripRecHit1D.h
ValueMap.h
reco::ClusterRemovalInfo::pixelRefProd
const SiPixelClusterRefProd & pixelRefProd() const
Definition: ClusterRemovalInfo.h:57
HITrackClusterRemover::oldStrMaskToken_
edm::EDGetTokenT< StripMaskContainer > oldStrMaskToken_
Definition: HITrackClusterRemover.cc:78
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
HITrackClusterRemover::ParamBlock::ParamBlock
ParamBlock()
Definition: HITrackClusterRemover.cc:50
HITrackClusterRemover::mergeOld_
bool mergeOld_
Definition: HITrackClusterRemover.cc:70
HITrackClusterRemover::mergeOld
void mergeOld(reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices &oldRefs)
Definition: HITrackClusterRemover.cc:244
GeomDet.h
edmNew::DetSetVector
Definition: DetSetNew.h:13
edm::Ref::id
ProductID id() const
Accessor for product ID.
Definition: Ref.h:244
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
edm::OrphanHandle
Definition: EDProductfwd.h:39
edmNew::DetSetVector::data
data_type const * data(size_t cell) const
Definition: DetSetVectorNew.h:602
SiStripMatchedRecHit2D
Definition: SiStripMatchedRecHit2D.h:8
edm::ContainerMask::copyMaskTo
void copyMaskTo(std::vector< bool > &) const
Definition: ContainerMask.h:87
reco::ClusterRemovalInfo::pixelNewRefProd
const SiPixelClusterRefProd & pixelNewRefProd() const
Definition: ClusterRemovalInfo.h:60
DetId.h
Frameworkfwd.h
T
long double T
Definition: Basic3DVectorLD.h:48
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:577
edm::ValueMap< int >
Exception
Definition: hltDiff.cc:246
reco::ClusterRemovalInfo::stripIndices
Indices & stripIndices()
Definition: ClusterRemovalInfo.h:50
HITrackClusterRemover::cleanup
std::unique_ptr< edmNew::DetSetVector< T > > cleanup(const edmNew::DetSetVector< T > &oldClusters, const std::vector< uint8_t > &isGood, reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices *oldRefs)
Definition: HITrackClusterRemover.cc:251
SiStripDetId::TIB
static constexpr auto TIB
Definition: SiStripDetId.h:37
HITrackClusterRemover::pixels
std::vector< uint8_t > pixels
Definition: HITrackClusterRemover.cc:93
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
HITrackClusterRemover::PixelMaskContainer
edm::ContainerMask< edmNew::DetSetVector< SiPixelCluster > > PixelMaskContainer
Definition: HITrackClusterRemover.cc:71
HITrackClusterRemover::stereoRecHitToken_
edm::EDGetTokenT< SiStripRecHit2DCollection > stereoRecHitToken_
Definition: HITrackClusterRemover.cc:80
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
edmNew::DetSetVector::FastFiller
Definition: DetSetVectorNew.h:236
HITrackClusterRemover::doPixel_
bool doPixel_
Definition: HITrackClusterRemover.cc:69
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:250
HITrackClusterRemover::doStrip_
bool doStrip_
Definition: HITrackClusterRemover.cc:69
HLT_2018_cff.track
track
Definition: HLT_2018_cff.py:10352
cms::Exception
Definition: Exception.h:70
ParameterSet.h
SiStripDetId
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:18
SiStripMatchedRecHit2D.h
globals_cff.id2
id2
Definition: globals_cff.py:33
HITrackClusterRemover::clusterWasteSolution_
bool clusterWasteSolution_
Definition: HITrackClusterRemover.cc:108
edm::Event
Definition: Event.h:73
TrackerSingleRecHit::omniClusterRef
OmniClusterRef const & omniClusterRef() const
Definition: TrackerSingleRecHit.h:41
SiStripRecHit2D.h
DetSetVectorNew.h
edm::InputTag
Definition: InputTag.h:15
label
const char * label
Definition: PFTauDecayModeTools.cc:11
HITrackClusterRemover::ParamBlock::maxSize_
size_t maxSize_
Definition: HITrackClusterRemover.cc:64
HITrackClusterRemover::ParamBlock
Definition: HITrackClusterRemover.cc:49
edm::ProductID
Definition: ProductID.h:27
hit
Definition: SiStripHitEffFromCalibTree.cc:88
edm::ProductLabels
Definition: ProductLabels.h:4
TrackerGeometry
Definition: TrackerGeometry.h:14
LaserSeedGenerator_cfi.rphiRecHits
rphiRecHits
Definition: LaserSeedGenerator_cfi.py:9
Generated for CMSSW Reference Manual by doxygen 1.8.16