#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) override

virtual void	beginJob ()

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

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

	SiPixelDigitizer (const edm::ParameterSet &conf, edm::EDProducer &mixMod)

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 > >)

Private Attributes
std::unique_ptr < SiPixelDigitizerAlgorithm >	_pixeldigialgo

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

bool	first

const std::string	geometryType

const std::string	hitsProducer

edm::ESHandle< TrackerGeometry >	pDD

edm::ESHandle< MagneticField >	pSetup

CLHEP::HepRandomEngine *	rndEngine

const vstring	trackerContainers

Detailed Description

Definition at line 44 of file SiPixelDigitizer.h.

Member Typedef Documentation

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

private

Definition at line 61 of file SiPixelDigitizer.h.

Constructor & Destructor Documentation

SiPixelDigitizer::SiPixelDigitizer	(	const edm::ParameterSet &	conf,
		edm::EDProducer &	mixMod
	)

explicit

Definition at line 89 of file SiPixelDigitizer.cc.

References _pixeldigialgo, edm::hlt::Exception, edm::RandomNumberGenerator::getEngine(), edm::Service< T >::isAvailable(), rndEngine, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                          :
     first(true),
     _pixeldigialgo(),
     hitsProducer(iConfig.getParameter<std::string>("hitsProducer")),
     trackerContainers(iConfig.getParameter<std::vector<std::string> >("ROUList")),
     geometryType(iConfig.getParameter<std::string>("GeometryType"))
   {
     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");
     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.";
     }
   
     rndEngine       = &(rng->getEngine());
     _pixeldigialgo.reset(new SiPixelDigitizerAlgorithm(iConfig,(*rndEngine)));
 
   }

SiPixelDigitizer::~SiPixelDigitizer ( )

virtual

Definition at line 115 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 174 of file SiPixelDigitizer.cc.

References accumulatePixelHits(), edm::Event::getByLabel(), hitsProducer, i, 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);
     }
   }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 186 of file SiPixelDigitizer.cc.

References accumulatePixelHits(), PileUpEventPrincipal::getByLabel(), hitsProducer, i, 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);
     }
   }

void SiPixelDigitizer::accumulatePixelHits ( edm::Handle< std::vector< PSimHit > > hSimHits )

private

Definition at line 125 of file SiPixelDigitizer.cc.

References _pixeldigialgo, detectorUnits, LogDebug, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, GloballyPositioned< T >::position(), pSetup, trackerHits::simHits, DetId::subdetId(), and GeomDet::surface().

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)) {
              PixelGeomDetUnit* pixdet = detectorUnits[detId];
              //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);
            }
          }
        }
     }
   }

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

inlinevirtual

Definition at line 56 of file SiPixelDigitizer.h.

56 {}

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 199 of file SiPixelDigitizer.cc.

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

                                                                                {
 
     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;
 
     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<PixelGeomDetUnit*>((*iu)),
                                  collector.data,
                                  linkcollector.data,
                  tTopo);
         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);
   }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 148 of file SiPixelDigitizer.cc.

References _pixeldigialgo, detectorUnits, first, geometryType, edm::EventSetup::get(), pDD, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, pSetup, and DetId::subdetId().

                                                                                   {
     if(first){
       _pixeldigialgo->init(iSetup);
       first = false;
     }
     _pixeldigialgo->initializeEvent();
     iSetup.get<TrackerDigiGeometryRecord>().get(geometryType, pDD);
     iSetup.get<IdealMagneticFieldRecord>().get(pSetup);
 
     // 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)) {  
           PixelGeomDetUnit* pixdet = dynamic_cast<PixelGeomDetUnit*>((*iu));
           assert(pixdet != 0);
           detectorUnits.insert(std::make_pair(detId, pixdet));
         }
       }
     }
   }