#include <SiPixelDigitizer.h>

Inheritance diagram for cms::SiPixelDigitizer:

Public Member Functions
void	accumulate (edm::Event const &e, edm::EventSetup const &c) override

void	accumulate (PileUpEventPrincipal const &e, edm::EventSetup const &c, edm::StreamID const &) override

virtual void	beginJob ()

void	finalizeEvent (edm::Event &e, edm::EventSetup const &c) override

PileupMixingContent *	getEventPileupInfo () override

void	initializeEvent (edm::Event const &e, edm::EventSetup const &c) override

	SiPixelDigitizer (const edm::ParameterSet &conf, edm::ProducesCollector, edm::ConsumesCollector &iC)

void	StorePileupInformation (std::vector< int > &numInteractionList, std::vector< int > &bunchCrossingList, std::vector< float > &TrueInteractionList, std::vector< edm::EventID > &eventInfoList, int bunchSpacing) override

	~SiPixelDigitizer () override

Public Member Functions inherited from DigiAccumulatorMixMod
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)

virtual void	beginRun (edm::Run const &run, edm::EventSetup const &setup)

	DigiAccumulatorMixMod ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)

virtual void	endRun (edm::Run const &run, edm::EventSetup const &setup)

virtual void	finalizeBunchCrossing (edm::Event &event, edm::EventSetup const &setup, int bunchCrossing)

virtual void	initializeBunchCrossing (edm::Event const &event, edm::EventSetup const &setup, int bunchCrossing)

virtual	~DigiAccumulatorMixMod ()

Private Types
typedef std::vector< std::string >	vstring

Private Member Functions
void	accumulatePixelHits (edm::Handle< std::vector< PSimHit > >, size_t globalSimHitIndex, const unsigned int tofBin, edm::EventSetup const &c)

Private Attributes
std::unique_ptr< SiPixelDigitizerAlgorithm >	_pixeldigialgo

std::map< std::string, size_t >	crossingSimHitIndexOffset_
	Offset to add to the index of each sim hit to account for which crossing it's in. More...

std::map< unsigned int, PixelGeomDetUnit const * >	detectorUnits

bool	firstFinalizeEvent_

bool	firstInitializeEvent_

const std::string	geometryType

const std::string	hitsProducer

const int	NumberOfEndcapDisks

edm::ESHandle< TrackerGeometry >	pDD

std::unique_ptr< PileupMixingContent >	PileupInfo_

const bool	pilotBlades

edm::ESHandle< MagneticField >	pSetup

CLHEP::HepRandomEngine *	randomEngine_ = nullptr

const vstring	trackerContainers

Detailed Description

Definition at line 48 of file SiPixelDigitizer.h.

Member Typedef Documentation

◆ vstring

typedef std::vector<std::string> cms::SiPixelDigitizer::vstring

private

Definition at line 91 of file SiPixelDigitizer.h.

Constructor & Destructor Documentation

◆ SiPixelDigitizer()

SiPixelDigitizer::SiPixelDigitizer	(	const edm::ParameterSet &	conf,
		edm::ProducesCollector	producesCollector,
		edm::ConsumesCollector &	iC
	)

explicit

Definition at line 85 of file SiPixelDigitizer.cc.

       : firstInitializeEvent_(true),
         firstFinalizeEvent_(true),
         _pixeldigialgo(),
         hitsProducer(iConfig.getParameter<std::string>("hitsProducer")),
         trackerContainers(iConfig.getParameter<std::vector<std::string> >("RoutList")),
         geometryType(iConfig.getParameter<std::string>("PixGeometryType")),
         pilotBlades(iConfig.exists("enablePilotBlades") ? iConfig.getParameter<bool>("enablePilotBlades") : false),
         NumberOfEndcapDisks(iConfig.exists("NumPixelEndcap") ? iConfig.getParameter<int>("NumPixelEndcap") : 2) {
     edm::LogInfo("PixelDigitizer ") << "Enter the Pixel Digitizer";
  
     const std::string alias("simSiPixelDigis");
  
     producesCollector.produces<edm::DetSetVector<PixelDigi> >().setBranchAlias(alias);
     producesCollector.produces<edm::DetSetVector<PixelDigiSimLink> >().setBranchAlias(alias + "siPixelDigiSimLink");
  
     for (auto const& trackerContainer : trackerContainers) {
       edm::InputTag tag(hitsProducer, trackerContainer);
       iC.consumes<std::vector<PSimHit> >(edm::InputTag(hitsProducer, trackerContainer));
     }
     edm::Service<edm::RandomNumberGenerator> rng;
     if (!rng.isAvailable()) {
       throw cms::Exception("Configuration")
           << "SiPixelDigitizer requires the RandomNumberGeneratorService\n"
              "which is not present in the configuration file.  You must add the service\n"
              "in the configuration file or remove the modules that require it.";
     }
  
     _pixeldigialgo = std::make_unique<SiPixelDigitizerAlgorithm>(iConfig);
     if (NumberOfEndcapDisks != 2)
       producesCollector.produces<PixelFEDChannelCollection>();
   }

References _pixeldigialgo, SiStripOfflineCRack_cfg::alias, edm::ConsumesCollector::consumes(), Exception, hitsProducer, HLT_FULL_cff::InputTag, edm::Service< T >::isAvailable(), NumberOfEndcapDisks, edm::ProducesCollector::produces(), AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

◆ ~SiPixelDigitizer()

SiPixelDigitizer::~SiPixelDigitizer ( )

override

Definition at line 120 of file SiPixelDigitizer.cc.

120 { edm::LogInfo("PixelDigitizer ") << "Destruct the Pixel Digitizer"; }

Member Function Documentation

◆ accumulate() [1/2]

void SiPixelDigitizer::accumulate	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 207 of file SiPixelDigitizer.cc.

                                                                                        {
     // Step A: Get Inputs
     for (vstring::const_iterator i = trackerContainers.begin(), iEnd = trackerContainers.end(); i != iEnd; ++i) {
       edm::Handle<std::vector<PSimHit> > simHits;
       edm::InputTag tag(hitsProducer, *i);
  
       iEvent.getByLabel(tag, simHits);
       unsigned int tofBin = PixelDigiSimLink::LowTof;
       if ((*i).find(std::string("HighTof")) != std::string::npos)
         tofBin = PixelDigiSimLink::HighTof;
       accumulatePixelHits(simHits, crossingSimHitIndexOffset_[tag.encode()], tofBin, iSetup);
       // Now that the hits have been processed, I'll add the amount of hits in this crossing on to
       // the global counter. Next time accumulateStripHits() is called it will count the sim hits
       // as though they were on the end of this collection.
       // Note that this is only used for creating digi-sim links (if configured to do so).
       //       std::cout << "index offset, current hit count = " << crossingSimHitIndexOffset_[tag.encode()] << ", " << simHits->size() << std::endl;
       if (simHits.isValid())
         crossingSimHitIndexOffset_[tag.encode()] += simHits->size();
     }
   }

References accumulatePixelHits(), crossingSimHitIndexOffset_, PixelDigiSimLink::HighTof, hitsProducer, mps_fire::i, iEvent, PixelDigiSimLink::LowTof, FastTrackerRecHitCombiner_cfi::simHits, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

◆ accumulate() [2/2]

void SiPixelDigitizer::accumulate	(	PileUpEventPrincipal const &	e,
		edm::EventSetup const &	c,
		edm::StreamID const &	streamID
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 228 of file SiPixelDigitizer.cc.

                                                                  {
     // Step A: Get Inputs
     for (vstring::const_iterator i = trackerContainers.begin(), iEnd = trackerContainers.end(); i != iEnd; ++i) {
       edm::Handle<std::vector<PSimHit> > simHits;
       edm::InputTag tag(hitsProducer, *i);
  
       iEvent.getByLabel(tag, simHits);
       unsigned int tofBin = PixelDigiSimLink::LowTof;
       if ((*i).find(std::string("HighTof")) != std::string::npos)
         tofBin = PixelDigiSimLink::HighTof;
       accumulatePixelHits(simHits, crossingSimHitIndexOffset_[tag.encode()], tofBin, iSetup);
       // Now that the hits have been processed, I'll add the amount of hits in this crossing on to
       // the global counter. Next time accumulateStripHits() is called it will count the sim hits
       // as though they were on the end of this collection.
       // Note that this is only used for creating digi-sim links (if configured to do so).
       //       std::cout << "index offset, current hit count = " << crossingSimHitIndexOffset_[tag.encode()] << ", " << simHits->size() << std::endl;
       if (simHits.isValid())
         crossingSimHitIndexOffset_[tag.encode()] += simHits->size();
     }
   }

References accumulatePixelHits(), crossingSimHitIndexOffset_, PixelDigiSimLink::HighTof, hitsProducer, mps_fire::i, iEvent, PixelDigiSimLink::LowTof, FastTrackerRecHitCombiner_cfi::simHits, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

◆ accumulatePixelHits()

void SiPixelDigitizer::accumulatePixelHits	(	edm::Handle< std::vector< PSimHit > >	hSimHits,
		size_t	globalSimHitIndex,
		const unsigned int	tofBin,
		edm::EventSetup const &	c
	)

private

Definition at line 126 of file SiPixelDigitizer.cc.

                                                                           {
     if (hSimHits.isValid()) {
       std::set<unsigned int> detIds;
       std::vector<PSimHit> const& simHits = *hSimHits.product();
       edm::ESHandle<TrackerTopology> tTopoHand;
       iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
       const TrackerTopology* tTopo = tTopoHand.product();
       for (std::vector<PSimHit>::const_iterator it = simHits.begin(), itEnd = simHits.end(); it != itEnd;
            ++it, ++globalSimHitIndex) {
         unsigned int detId = (*it).detUnitId();
         if (detIds.insert(detId).second) {
           // The insert succeeded, so this detector element has not yet been processed.
           assert(detectorUnits[detId]);
           if (detectorUnits[detId] &&
               detectorUnits[detId]
                   ->type()
                   .isTrackerPixel()) {  // this test could be avoided and changed into a check of pixdet!=0
             std::map<unsigned int, PixelGeomDetUnit const*>::iterator itDet = detectorUnits.find(detId);
             if (itDet == detectorUnits.end())
               continue;
             auto pixdet = itDet->second;
             assert(pixdet != nullptr);
             //access to magnetic field in global coordinates
             GlobalVector bfield = pSetup->inTesla(pixdet->surface().position());
             LogDebug("PixelDigitizer ") << "B-field(T) at " << pixdet->surface().position()
                                         << "(cm): " << pSetup->inTesla(pixdet->surface().position());
             _pixeldigialgo->accumulateSimHits(
                 it, itEnd, globalSimHitIndex, tofBin, pixdet, bfield, tTopo, randomEngine_);
           }
         }
       }
     }
   }

References _pixeldigialgo, cms::cuda::assert(), detectorUnits, edm::EventSetup::get(), get, MagneticField::inTesla(), GeomDetEnumerators::isTrackerPixel(), LogDebug, edm::ESHandle< T >::product(), pSetup, randomEngine_, and FastTrackerRecHitCombiner_cfi::simHits.

Referenced by accumulate().

◆ beginJob()

virtual void cms::SiPixelDigitizer::beginJob ( void )

inlinevirtual

Definition at line 59 of file SiPixelDigitizer.h.

59 {}

◆ finalizeEvent()

void SiPixelDigitizer::finalizeEvent	(	edm::Event &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 252 of file SiPixelDigitizer.cc.

                                                                                     {
     edm::ESHandle<TrackerTopology> tTopoHand;
     iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
     const TrackerTopology* tTopo = tTopoHand.product();
  
     std::vector<edm::DetSet<PixelDigi> > theDigiVector;
     std::vector<edm::DetSet<PixelDigiSimLink> > theDigiLinkVector;
  
     if (firstFinalizeEvent_) {
       const unsigned int bunchspace = PileupInfo_->getMix_bunchSpacing();
       _pixeldigialgo->init_DynIneffDB(iSetup, bunchspace);
       firstFinalizeEvent_ = false;
     }
     _pixeldigialgo->calculateInstlumiFactor(PileupInfo_.get());
  
     if (_pixeldigialgo->killBadFEDChannels()) {
       std::unique_ptr<PixelFEDChannelCollection> PixelFEDChannelCollection_ =
           _pixeldigialgo->chooseScenario(PileupInfo_.get(), randomEngine_);
       if (PixelFEDChannelCollection_ == nullptr) {
         throw cms::Exception("NullPointerError") << "PixelFEDChannelCollection not set in chooseScenario function.\n";
       }
       iEvent.put(std::move(PixelFEDChannelCollection_));
     }
  
     for (const auto& iu : pDD->detUnits()) {
       if (iu->type().isTrackerPixel()) {
         //
  
         edm::DetSet<PixelDigi> collector(iu->geographicalId().rawId());
         edm::DetSet<PixelDigiSimLink> linkcollector(iu->geographicalId().rawId());
  
         _pixeldigialgo->digitize(
             dynamic_cast<const PixelGeomDetUnit*>(iu), collector.data, linkcollector.data, tTopo, randomEngine_);
         if (!collector.data.empty()) {
           theDigiVector.push_back(std::move(collector));
         }
         if (!linkcollector.data.empty()) {
           theDigiLinkVector.push_back(std::move(linkcollector));
         }
       }
     }
  
     // Step C: create collection with the cache vector of DetSet
     std::unique_ptr<edm::DetSetVector<PixelDigi> > output(new edm::DetSetVector<PixelDigi>(theDigiVector));
     std::unique_ptr<edm::DetSetVector<PixelDigiSimLink> > outputlink(
         new edm::DetSetVector<PixelDigiSimLink>(theDigiLinkVector));
  
     // Step D: write output to file
     iEvent.put(std::move(output));
     iEvent.put(std::move(outputlink));
  
     randomEngine_ = nullptr;  // to prevent access outside event
   }

References _pixeldigialgo, MixingModule_cfi::bunchspace, TrackerGeometry::detUnits(), Exception, firstFinalizeEvent_, edm::EventSetup::get(), get, iEvent, eostools::move(), convertSQLitetoXML_cfg::output, pDD, PileupInfo_, edm::ESHandle< T >::product(), and randomEngine_.

◆ getEventPileupInfo()

PileupMixingContent* cms::SiPixelDigitizer::getEventPileupInfo ( )

inlineoverridevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 70 of file SiPixelDigitizer.h.

70 { return PileupInfo_.get(); }

References PileupInfo_.

◆ initializeEvent()

void SiPixelDigitizer::initializeEvent	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 163 of file SiPixelDigitizer.cc.

                                                                                        {
     if (firstInitializeEvent_) {
       _pixeldigialgo->init(iSetup);
       firstInitializeEvent_ = false;
     }
  
     // Make sure that the first crossing processed starts indexing the sim hits from zero.
     // This variable is used so that the sim hits from all crossing frames have sequential
     // indices used to create the digi-sim link (if configured to do so) rather than starting
     // from zero for each crossing.
     crossingSimHitIndexOffset_.clear();
  
     // Cache random number engine
     edm::Service<edm::RandomNumberGenerator> rng;
     randomEngine_ = &rng->getEngine(e.streamID());
  
     _pixeldigialgo->initializeEvent();
     iSetup.get<TrackerDigiGeometryRecord>().get(geometryType, pDD);
     iSetup.get<IdealMagneticFieldRecord>().get(pSetup);
     edm::ESHandle<TrackerTopology> tTopoHand;
     iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
     const TrackerTopology* tTopo = tTopoHand.product();
  
     // FIX THIS! We only need to clear and (re)fill this map when the geometry type IOV changes.  Use ESWatcher to determine this.
     if (true) {  // Replace with ESWatcher
       detectorUnits.clear();
       for (const auto& iu : pDD->detUnits()) {
         unsigned int detId = iu->geographicalId().rawId();
         if (iu->type().isTrackerPixel()) {
           auto pixdet = dynamic_cast<const PixelGeomDetUnit*>(iu);
           assert(pixdet != nullptr);
           if (iu->subDetector() ==
               GeomDetEnumerators::SubDetector::PixelEndcap) {  // true ONLY for the phase 0 pixel deetctor
             unsigned int disk = tTopo->layer(detId);           // using the generic layer method
             //if using pilot blades, then allowing it for current detector only
             if ((disk == 3) && ((!pilotBlades) && (NumberOfEndcapDisks == 2)))
               continue;
           }
           detectorUnits.insert(std::make_pair(detId, pixdet));
         }
       }
     }
   }

References _pixeldigialgo, cms::cuda::assert(), crossingSimHitIndexOffset_, detectorUnits, TrackerGeometry::detUnits(), MillePedeFileConverter_cfg::e, firstInitializeEvent_, geometryType, edm::EventSetup::get(), get, edm::RandomNumberGenerator::getEngine(), TrackerTopology::layer(), NumberOfEndcapDisks, pDD, pilotBlades, GeomDetEnumerators::PixelEndcap, edm::ESHandle< T >::product(), pSetup, and randomEngine_.

◆ StorePileupInformation()

void cms::SiPixelDigitizer::StorePileupInformation	(	std::vector< int > &	numInteractionList,
		std::vector< int > &	bunchCrossingList,
		std::vector< float > &	TrueInteractionList,
		std::vector< edm::EventID > &	eventInfoList,
		int	bunchSpacing
	)

inlineoverridevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 61 of file SiPixelDigitizer.h.

                                                            {
       PileupInfo_ = std::make_unique<PileupMixingContent>(
           numInteractionList, bunchCrossingList, TrueInteractionList, eventInfoList, bunchSpacing);
     }

References genPUProtons_cfi::bunchCrossingList, HLTEgPhaseIITestSequence_cff::bunchSpacing, and PileupInfo_.

Member Data Documentation

◆ _pixeldigialgo

std::unique_ptr<SiPixelDigitizerAlgorithm> cms::SiPixelDigitizer::_pixeldigialgo

private

Definition at line 80 of file SiPixelDigitizer.h.

Referenced by accumulatePixelHits(), finalizeEvent(), initializeEvent(), and SiPixelDigitizer().

◆ crossingSimHitIndexOffset_

std::map<std::string, size_t> cms::SiPixelDigitizer::crossingSimHitIndexOffset_

private

Offset to add to the index of each sim hit to account for which crossing it's in.

I need to know what each sim hit index will be when the hits from all crossing frames are merged into one collection (assuming the MixingModule is configured to create the crossing frame for all sim hits). To do this I'll record how many hits were in each crossing, and then add that on to the index for a given hit in a given crossing. This assumes that the crossings are processed in the same order here as they are put into the crossing frame, which I'm pretty sure is true.
The key is the name of the sim hit collection.

Definition at line 89 of file SiPixelDigitizer.h.

Referenced by accumulate(), and initializeEvent().

◆ detectorUnits

std::map<unsigned int, PixelGeomDetUnit const*> cms::SiPixelDigitizer::detectorUnits

private

Definition at line 97 of file SiPixelDigitizer.h.

Referenced by accumulatePixelHits(), and initializeEvent().

◆ firstFinalizeEvent_

bool cms::SiPixelDigitizer::firstFinalizeEvent_

private

Definition at line 79 of file SiPixelDigitizer.h.

Referenced by finalizeEvent().

◆ firstInitializeEvent_

bool cms::SiPixelDigitizer::firstInitializeEvent_

private

Definition at line 78 of file SiPixelDigitizer.h.

Referenced by initializeEvent().

◆ geometryType

const std::string cms::SiPixelDigitizer::geometryType

private

Definition at line 94 of file SiPixelDigitizer.h.

Referenced by initializeEvent().

◆ hitsProducer

const std::string cms::SiPixelDigitizer::hitsProducer

private

Definition at line 92 of file SiPixelDigitizer.h.

Referenced by accumulate(), and SiPixelDigitizer().

◆ NumberOfEndcapDisks

const int cms::SiPixelDigitizer::NumberOfEndcapDisks

private

Definition at line 103 of file SiPixelDigitizer.h.

Referenced by initializeEvent(), and SiPixelDigitizer().

◆ pDD

edm::ESHandle<TrackerGeometry> cms::SiPixelDigitizer::pDD

private

Definition at line 95 of file SiPixelDigitizer.h.

Referenced by finalizeEvent(), and initializeEvent().

◆ PileupInfo_

std::unique_ptr<PileupMixingContent> cms::SiPixelDigitizer::PileupInfo_

private

Definition at line 100 of file SiPixelDigitizer.h.

Referenced by finalizeEvent(), getEventPileupInfo(), and StorePileupInformation().

◆ pilotBlades

const bool cms::SiPixelDigitizer::pilotBlades

private

Definition at line 102 of file SiPixelDigitizer.h.

Referenced by initializeEvent().

◆ pSetup

edm::ESHandle<MagneticField> cms::SiPixelDigitizer::pSetup

private

Definition at line 96 of file SiPixelDigitizer.h.

Referenced by accumulatePixelHits(), and initializeEvent().

◆ randomEngine_

CLHEP::HepRandomEngine* cms::SiPixelDigitizer::randomEngine_ = nullptr

private

Definition at line 98 of file SiPixelDigitizer.h.

Referenced by accumulatePixelHits(), finalizeEvent(), and initializeEvent().

◆ trackerContainers

const vstring cms::SiPixelDigitizer::trackerContainers

private

Definition at line 93 of file SiPixelDigitizer.h.

Referenced by accumulate(), and SiPixelDigitizer().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ vstring

Constructor & Destructor Documentation

◆ SiPixelDigitizer()

◆ ~SiPixelDigitizer()

Member Function Documentation

◆ accumulate() [1/2]

◆ accumulate() [2/2]

◆ accumulatePixelHits()

◆ beginJob()

◆ finalizeEvent()

◆ getEventPileupInfo()

◆ initializeEvent()

◆ StorePileupInformation()

Member Data Documentation

◆ _pixeldigialgo

◆ crossingSimHitIndexOffset_

◆ detectorUnits

◆ firstFinalizeEvent_

◆ firstInitializeEvent_

◆ geometryType

◆ hitsProducer

◆ NumberOfEndcapDisks

◆ pDD

◆ PileupInfo_

◆ pilotBlades

◆ pSetup

◆ randomEngine_

◆ trackerContainers