#include <SiPixelDigitizer.h>

Inheritance diagram for cms::SiPixelDigitizer:

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

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

virtual void	beginJob ()

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

virtual PileupMixingContent *	getEventPileupInfo ()

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

	SiPixelDigitizer (const edm::ParameterSet &conf, edm::one::EDProducerBase &mixMod, edm::ConsumesCollector &iC)

virtual void	StorePileupInformation (std::vector< int > &numInteractionList, std::vector< int > &bunchCrossingList, std::vector< float > &TrueInteractionList)

virtual	~SiPixelDigitizer ()

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 > >, CLHEP::HepRandomEngine *)

CLHEP::HepRandomEngine *	randomEngine (edm::StreamID const &streamID)

Private Attributes
std::unique_ptr < SiPixelDigitizerAlgorithm >	_pixeldigialgo

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

bool	first

const std::string	geometryType

const std::string	hitsProducer

const int	NumberOfEndcapDisks

edm::ESHandle< TrackerGeometry >	pDD

PileupMixingContent *	PileupInfo_

const bool	pilotBlades

edm::ESHandle< MagneticField >	pSetup

std::vector < CLHEP::HepRandomEngine * >	randomEngines_

const vstring	trackerContainers

Detailed Description

Definition at line 48 of file SiPixelDigitizer.h.

Member Typedef Documentation

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

private

Definition at line 76 of file SiPixelDigitizer.h.

Constructor & Destructor Documentation

SiPixelDigitizer::SiPixelDigitizer	(	const edm::ParameterSet &	conf,
		edm::one::EDProducerBase &	mixMod,
		edm::ConsumesCollector &	iC
	)

explicit

Definition at line 95 of file SiPixelDigitizer.cc.

References _pixeldigialgo, edm::ConsumesCollector::consumes(), edm::hlt::Exception, hitsProducer, edm::Service< T >::isAvailable(), AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

                                                                                                                           :
     first(true),
     _pixeldigialgo(),
     hitsProducer(iConfig.getParameter<std::string>("hitsProducer")),
     trackerContainers(iConfig.getParameter<std::vector<std::string> >("ROUList")),
     geometryType(iConfig.getParameter<std::string>("GeometryType")),
     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"); 
     
     mixMod.produces<edm::DetSetVector<PixelDigi> >().setBranchAlias(alias);
     mixMod.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.reset(new SiPixelDigitizerAlgorithm(iConfig));
   }

SiPixelDigitizer::~SiPixelDigitizer ( )

virtual

Definition at line 125 of file SiPixelDigitizer.cc.

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

Member Function Documentation

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 194 of file SiPixelDigitizer.cc.

References accumulatePixelHits(), edm::Event::getByLabel(), hitsProducer, i, randomEngine(), trackerHits::simHits, edm::Event::streamID(), GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

                                                                                   {
     // 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);
       accumulatePixelHits(simHits, randomEngine(iEvent.streamID()));
     }
   }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 206 of file SiPixelDigitizer.cc.

References accumulatePixelHits(), PileUpEventPrincipal::getByLabel(), hitsProducer, i, randomEngine(), trackerHits::simHits, GlobalPosition_Frontier_DevDB_cff::tag, and trackerContainers.

                                                                                                                            {
     // 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);
       accumulatePixelHits(simHits, randomEngine(streamID));
     }
   }

void SiPixelDigitizer::accumulatePixelHits	(	edm::Handle< std::vector< PSimHit > >	hSimHits,
		CLHEP::HepRandomEngine *	engine
	)

private

Definition at line 135 of file SiPixelDigitizer.cc.

References _pixeldigialgo, detectorUnits, LogDebug, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pSetup, trackerHits::simHits, and DetId::subdetId().

Referenced by accumulate().

                                                                                                              {
     if(hSimHits.isValid()) {
        std::set<unsigned int> detIds;
        std::vector<PSimHit> const& simHits = *hSimHits.product();
        for(std::vector<PSimHit>::const_iterator it = simHits.begin(), itEnd = simHits.end(); it != itEnd; ++it) {
          unsigned int detId = (*it).detUnitId();
          if(detIds.insert(detId).second) {
            // The insert succeeded, so this detector element has not yet been processed.
            unsigned int isub = DetId(detId).subdetId();
            if((isub == PixelSubdetector::PixelBarrel) || (isub == PixelSubdetector::PixelEndcap)) {
          std::map<unsigned int, PixelGeomDetUnit const *>::iterator itDet = detectorUnits.find(detId);       
          if (itDet == detectorUnits.end()) continue;
              auto pixdet = itDet->second;
              //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, pixdet, bfield, engine);
            }
          }
        }
     }
   }

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

inlinevirtual

Definition at line 60 of file SiPixelDigitizer.h.

60 {}

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 219 of file SiPixelDigitizer.cc.

References _pixeldigialgo, edm::EventSetup::get(), edm::RandomNumberGenerator::getEngine(), getEventPileupInfo(), convertSQLitetoXML_cfg::output, pDD, PileupInfo_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::ESHandle< class >::product(), edm::Event::put(), edm::Event::streamID(), and DetId::subdetId().

                                                                                {
 
     edm::Service<edm::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());
 
     edm::ESHandle<TrackerTopology> tTopoHand;
     iSetup.get<IdealGeometryRecord>().get(tTopoHand);
     const TrackerTopology *tTopo=tTopoHand.product();
 
     std::vector<edm::DetSet<PixelDigi> > theDigiVector;
     std::vector<edm::DetSet<PixelDigiSimLink> > theDigiLinkVector;
  
     PileupInfo_ = getEventPileupInfo();
     _pixeldigialgo->calculateInstlumiFactor(PileupInfo_);   
 
     for(TrackingGeometry::DetUnitContainer::const_iterator iu = pDD->detUnits().begin(); iu != pDD->detUnits().end(); iu ++){
       DetId idet=DetId((*iu)->geographicalId().rawId());
       unsigned int isub=idet.subdetId();
       
       if((isub == PixelSubdetector::PixelBarrel) || (isub == PixelSubdetector::PixelEndcap)) {
 
     //
 
         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,
                                  engine);
         if(collector.data.size() > 0) {
           theDigiVector.push_back(std::move(collector));
         }
         if(linkcollector.data.size() > 0) {
           theDigiLinkVector.push_back(std::move(linkcollector));
         }
       }
     }
     
     // Step C: create collection with the cache vector of DetSet 
     std::auto_ptr<edm::DetSetVector<PixelDigi> > 
       output(new edm::DetSetVector<PixelDigi>(theDigiVector) );
     std::auto_ptr<edm::DetSetVector<PixelDigiSimLink> > 
       outputlink(new edm::DetSetVector<PixelDigiSimLink>(theDigiLinkVector) );
 
     // Step D: write output to file 
     iEvent.put(output);
     iEvent.put(outputlink);
   }

virtual PileupMixingContent* cms::SiPixelDigitizer::getEventPileupInfo ( )

inlinevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 68 of file SiPixelDigitizer.h.

References PileupInfo_.

Referenced by finalizeEvent().

68 { return PileupInfo_; }

cms::SiPixelDigitizer::PileupInfo_

PileupMixingContent * PileupInfo_

Definition: SiPixelDigitizer.h:85

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 160 of file SiPixelDigitizer.cc.

References _pixeldigialgo, detectorUnits, first, geometryType, edm::EventSetup::get(), NumberOfEndcapDisks, pDD, pilotBlades, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::ESHandle< class >::product(), pSetup, TrackerTopology::pxfDisk(), and DetId::subdetId().

                                                                                   {
     if(first){
       _pixeldigialgo->init(iSetup);
       first = false;
     }
     _pixeldigialgo->initializeEvent();
     iSetup.get<TrackerDigiGeometryRecord>().get(geometryType, pDD);
     iSetup.get<IdealMagneticFieldRecord>().get(pSetup);
     edm::ESHandle<TrackerTopology> tTopoHand;
     iSetup.get<IdealGeometryRecord>().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(TrackingGeometry::DetUnitContainer::const_iterator iu = pDD->detUnits().begin(); iu != pDD->detUnits().end(); ++iu) {
         unsigned int detId = (*iu)->geographicalId().rawId();
         DetId idet=DetId(detId);
         unsigned int isub=idet.subdetId();
         if((isub == PixelSubdetector::PixelBarrel) || (isub == PixelSubdetector::PixelEndcap)) {  
           auto pixdet = dynamic_cast<const PixelGeomDetUnit*>((*iu));
           assert(pixdet != 0);
       if (isub==PixelSubdetector::PixelEndcap) {
         unsigned int disk = tTopo->pxfDisk(detId);
         //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));
         }
       }
     }
   }

CLHEP::HepRandomEngine * SiPixelDigitizer::randomEngine ( edm::StreamID const & streamID )

private

Definition at line 271 of file SiPixelDigitizer.cc.

References edm::RandomNumberGenerator::getEngine(), cmsHarvester::index, randomEngines_, and edm::StreamID::value().

Referenced by accumulate().

                                                                                   {
     unsigned int index = streamID.value();
     if(index >= randomEngines_.size()) {
       randomEngines_.resize(index + 1, nullptr);
     }
     CLHEP::HepRandomEngine* ptr = randomEngines_[index];
     if(!ptr) {
       edm::Service<edm::RandomNumberGenerator> rng;
       ptr = &rng->getEngine(streamID);
       randomEngines_[index] = ptr;
     }
     return ptr;
   }

virtual void cms::SiPixelDigitizer::StorePileupInformation	(	std::vector< int > &	numInteractionList,
		std::vector< int > &	bunchCrossingList,
		std::vector< float > &	TrueInteractionList
	)

inlinevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 62 of file SiPixelDigitizer.h.

References PileupInfo_.

                                                        {
       PileupInfo_ = new PileupMixingContent(numInteractionList, bunchCrossingList, TrueInteractionList);
     }