CMS 3D CMS Logo

PixelPSimHitSelector.cc
Go to the documentation of this file.
1 
4 
6 
9 
11 
12 
14 {
15  // Look for psimhit collection associated o the tracker
16  PSimHitCollectionMap::const_iterator pSimHitCollections = pSimHitCollectionMap_.find("pixel");
17 
18  // Check that there are psimhit collections defined for the tracker
19  if (pSimHitCollections == pSimHitCollectionMap_.end()) return;
20 
21  // Grab all the PSimHit from the different sencitive volumes
23  std::vector<const CrossingFrame<PSimHit> *> cfPSimHitProductPointers;
24 
25  // Collect the product pointers to the different psimhit collection
26  for (std::size_t i = 0; i < pSimHitCollections->second.size(); ++i)
27  {
28  event.getByLabel("mix", pSimHitCollections->second[i], cfPSimHits);
29  cfPSimHitProductPointers.push_back(cfPSimHits.product());
30  }
31 
32  // Create a mix collection from the different psimhit collections
33  std::unique_ptr<MixCollection<PSimHit> > pSimHits( new MixCollection<PSimHit>(cfPSimHitProductPointers) );
34 
35  // Accessing dead pixel modules from DB:
36  edm::ESHandle<SiPixelQuality> siPixelBadModule;
37  setup.get<SiPixelQualityRcd>().get(siPixelBadModule);
38 
39  // Reading the DB information
40  std::vector<SiPixelQuality::disabledModuleType> badModules( siPixelBadModule->getBadComponentList() );
41  SiPixelQuality pixelQuality(badModules);
42 
43  // Select only psimhits from alive modules
44  for (MixCollection<PSimHit>::MixItr pSimHit = pSimHits->begin(); pSimHit != pSimHits->end(); ++pSimHit)
45  {
46  if ( !pixelQuality.IsModuleBad(pSimHit->detUnitId()) )
47  selection.push_back(*pSimHit);
48  }
49 }
selection
main part
Definition: corrVsCorr.py:98
def setup(process, global_tag, zero_tesla=False)
Definition: GeneralSetup.py:1
void select(PSimHitCollection &, edm::Event const &, edm::EventSetup const &) const override
Pre-process event information.
const std::vector< disabledModuleType > getBadComponentList() const
std::vector< PSimHit > PSimHitCollection
T const * product() const
Definition: Handle.h:81
const T & get() const
Definition: EventSetup.h:55
PSimHitCollectionMap pSimHitCollectionMap_
Definition: event.py:1