#include <GEMDigiProducer.h>

Inheritance diagram for GEMDigiProducer:

Public Types
typedef edm::DetSetVector< GEMDigiSimLink >	GEMDigiSimLinks

typedef edm::DetSetVector< StripDigiSimLink >	StripDigiSimLinks

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

	GEMDigiProducer (const edm::ParameterSet &ps)

void	produce (edm::Event &, const edm::EventSetup &) override

	~GEMDigiProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
edm::EDGetTokenT< CrossingFrame< PSimHit > >	cf_token

std::string	digiModelString_

GEMDigiModel *	gemDigiModel_

Detailed Description

Definition at line 20 of file GEMDigiProducer.h.

Member Typedef Documentation

typedef edm::DetSetVector<GEMDigiSimLink> GEMDigiProducer::GEMDigiSimLinks

Definition at line 26 of file GEMDigiProducer.h.

typedef edm::DetSetVector<StripDigiSimLink> GEMDigiProducer::StripDigiSimLinks

Definition at line 24 of file GEMDigiProducer.h.

Constructor & Destructor Documentation

GEMDigiProducer::GEMDigiProducer ( const edm::ParameterSet & ps )

explicit

Definition at line 29 of file GEMDigiProducer.cc.

References cf_token, digiModelString_, Exception, gemDigiModel_, reco::get(), edm::ParameterSet::getParameter(), LogDebug, and AlCaHLTBitMon_QueryRunRegistry::string.

   : digiModelString_(ps.getParameter<std::string>("digiModelString"))
 {
   produces<GEMDigiCollection>();
   produces<StripDigiSimLinks>("GEM");
   produces<GEMDigiSimLinks>("GEM");
 
   edm::Service<edm::RandomNumberGenerator> rng;
   if (!rng.isAvailable()){
     throw cms::Exception("Configuration")
       << "GEMDigiProducer::GEMDigiProducer() - RandomNumberGeneratorService is not present in configuration file.\n"
       << "Add the service in the configuration file or remove the modules that require it.";
   }
   gemDigiModel_ = GEMDigiModelFactory::get()->create("GEM" + digiModelString_ + "Model", ps);
   LogDebug("GEMDigiProducer") << "Using GEM" + digiModelString_ + "Model";
 
   std::string mix_(ps.getParameter<std::string>("mixLabel"));
   std::string collection_(ps.getParameter<std::string>("inputCollection"));
 
   cf_token = consumes<CrossingFrame<PSimHit> >(edm::InputTag(mix_, collection_));
 }

GEMDigiProducer::~GEMDigiProducer ( )

override

Definition at line 52 of file GEMDigiProducer.cc.

References gemDigiModel_.

 {
   delete gemDigiModel_;
 }

Member Function Documentation

void GEMDigiProducer::beginRun	(	const edm::Run &	,
		const edm::EventSetup &	eventSetup
	)

override

Definition at line 58 of file GEMDigiProducer.cc.

References gemDigiModel_, edm::EventSetup::get(), GEMDigiModel::setGeometry(), and GEMDigiModel::setup().

 {
   edm::ESHandle<GEMGeometry> hGeom;
   eventSetup.get<MuonGeometryRecord>().get(hGeom);
   gemDigiModel_->setGeometry(&*hGeom);
   gemDigiModel_->setup();
 }

void GEMDigiProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	eventSetup
	)

override

Definition at line 67 of file GEMDigiProducer.cc.

References cf_token, GEMGeometry::etaPartitions(), GEMDigiModel::fillDigis(), gemDigiModel_, GEMDigiModel::gemDigiSimLinks(), edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), GEMDigiModel::getGeometry(), hfClusterShapes_cfi::hits, LogDebug, eostools::move(), edm::Handle< T >::product(), edm::Event::put(), trackerHits::simHits, GEMDigiModel::simulateNoise(), GEMDigiModel::simulateSignal(), edm::Event::streamID(), and GEMDigiModel::stripDigiSimLinks().

 {
   edm::Service<edm::RandomNumberGenerator> rng;
   CLHEP::HepRandomEngine* engine = &rng->getEngine(e.streamID());
 
   edm::Handle<CrossingFrame<PSimHit> > cf;
   e.getByToken(cf_token, cf);
 
   std::unique_ptr<MixCollection<PSimHit> > hits(new MixCollection<PSimHit>(cf.product()));
 
   // Create empty output
   std::unique_ptr<GEMDigiCollection> digis(new GEMDigiCollection());
   std::unique_ptr<StripDigiSimLinks> stripDigiSimLinks(new StripDigiSimLinks() );
   std::unique_ptr<GEMDigiSimLinks> gemDigiSimLinks(new GEMDigiSimLinks() );
 
   // arrange the hits by eta partition
   std::map<uint32_t, edm::PSimHitContainer> hitMap;
   for (const auto& hit: *hits){
     hitMap[hit.detUnitId()].emplace_back(hit);
   }
 
   // simulate signal and noise for each eta partition
   const auto & etaPartitions(gemDigiModel_->getGeometry()->etaPartitions());
 
   for (const auto& roll: etaPartitions){
     const GEMDetId detId(roll->id());
     const uint32_t rawId(detId.rawId());
     const auto & simHits(hitMap[rawId]);
 
     LogDebug("GEMDigiProducer")
       << "GEMDigiProducer: found " << simHits.size() << " hit(s) in eta partition" << rawId;
 
     gemDigiModel_->simulateSignal(roll, simHits, engine);
     gemDigiModel_->simulateNoise(roll, engine);
     gemDigiModel_->fillDigis(rawId, *digis);
     (*stripDigiSimLinks).insert(gemDigiModel_->stripDigiSimLinks());
     (*gemDigiSimLinks).insert(gemDigiModel_->gemDigiSimLinks());
   }
 
   // store them in the event
   e.put(std::move(digis));
   e.put(std::move(stripDigiSimLinks),"GEM");
   e.put(std::move(gemDigiSimLinks),"GEM");
 }