#include <GEMDigiModule.h>

Public Types
typedef edm::DetSet < GEMDigiSimLink >	GEMDigiSimLinks

typedef edm::DetSet < StripDigiSimLink >	StripDigiSimLinks

Public Member Functions
void	fillDigis (int rollDetId, GEMDigiCollection &)

	GEMDigiModule (const edm::ParameterSet &)

const GEMDigiSimLinks &	gemDigiSimLinks () const

void	setGeometry (const GEMGeometry *)

void	simulate (const GEMEtaPartition , const edm::PSimHitContainer &, CLHEP::HepRandomEngine )

const StripDigiSimLinks &	stripDigiSimLinks () const

	~GEMDigiModule ()

Private Member Functions
void	addLinks (unsigned int strip, int bx)
	creates links from Digi to SimTrack More...

void	addLinksWithPartId (unsigned int strip, int bx)

Private Attributes
DetectorHitMap	detectorHitMap_

const GEMGeometry *	geometry_

std::vector< std::unique_ptr < GEMDigiModel > >	models

StripDigiSimLinks	stripDigiSimLinks_

Strips	strips_

GEMDigiSimLinks	theGemDigiSimLinks_

Detailed Description

Base Class for the GEM strip response simulation

Author: Sven Dildick by Yechan Kang

Definition at line 35 of file GEMDigiModule.h.

Member Typedef Documentation

typedef edm::DetSet<GEMDigiSimLink> GEMDigiModule::GEMDigiSimLinks

Definition at line 42 of file GEMDigiModule.h.

typedef edm::DetSet<StripDigiSimLink> GEMDigiModule::StripDigiSimLinks

Definition at line 41 of file GEMDigiModule.h.

Constructor & Destructor Documentation

GEMDigiModule::GEMDigiModule ( const edm::ParameterSet & config )

Definition at line 9 of file GEMDigiModule.cc.

References edm::ParameterSet::getParameter(), and models.

                                                           {
   bool simulateBkgNoise_(config.getParameter<bool>("simulateBkgNoise"));
   bool simulateIntrinsicNoise_(config.getParameter<bool>("simulateIntrinsicNoise"));
   if (simulateIntrinsicNoise_) {
     models.push_back(std::make_unique<GEMNoiseModel>(config));
   }
   if (simulateBkgNoise_) {
     models.push_back(std::make_unique<GEMBkgModel>(config));
   }
   models.push_back(std::make_unique<GEMSignalModel>(config));
 }

GEMDigiModule::~GEMDigiModule ( )

default

Member Function Documentation

void GEMDigiModule::addLinks	(	unsigned int	strip,
		int	bx
	)

private

creates links from Digi to SimTrack

Definition at line 50 of file GEMDigiModule.cc.

References RecoTauCleanerPlugins::charge, detectorHitMap_, PSimHit::eventId(), validate-o2o-wbm::f, MainPageGenerator::link, edm::DetSet< T >::push_back(), stripDigiSimLinks_, and PSimHit::trackId().

Referenced by fillDigis().

                                                        {
   std::pair<unsigned int, int> digi(strip, bx);
   auto channelHitItr = detectorHitMap_.equal_range(digi);
 
   // find the fraction contribution for each SimTrack
   std::map<int, float> simTrackChargeMap;
   std::map<int, EncodedEventId> eventIdMap;
   float totalCharge(0.);
   for (auto hitItr = channelHitItr.first; hitItr != channelHitItr.second; ++hitItr) {
     const PSimHit* hit(hitItr->second);
     // might be zero for unit tests and such
     if (hit == nullptr)
       continue;
 
     int simTrackId(hit->trackId());
     //float charge = hit->getCharge();
     const float charge(1.f);
     auto chargeItr = simTrackChargeMap.find(simTrackId);
     if (chargeItr == simTrackChargeMap.end()) {
       simTrackChargeMap[simTrackId] = charge;
       eventIdMap[simTrackId] = hit->eventId();
     } else {
       chargeItr->second += charge;
     }
     totalCharge += charge;
   }
 
   for (const auto& charge : simTrackChargeMap) {
     const int simTrackId(charge.first);
     auto link(StripDigiSimLink(strip, simTrackId, eventIdMap[simTrackId], charge.second / totalCharge));
     stripDigiSimLinks_.push_back(link);
   }
 }

void GEMDigiModule::addLinksWithPartId	(	unsigned int	strip,
		int	bx
	)

private

Definition at line 84 of file GEMDigiModule.cc.

References detectorHitMap_, PSimHit::eventId(), PSimHit::particleType(), edm::DetSet< T >::push_back(), theGemDigiSimLinks_, and PSimHit::trackId().

Referenced by fillDigis().

                                                                  {
   std::pair<unsigned int, int> digi(strip, bx);
   std::pair<DetectorHitMap::iterator, DetectorHitMap::iterator> channelHitItr = detectorHitMap_.equal_range(digi);
 
   for (DetectorHitMap::iterator hitItr = channelHitItr.first; hitItr != channelHitItr.second; ++hitItr) {
     const PSimHit* hit = (hitItr->second);
     // might be zero for unit tests and such
     if (hit == nullptr)
       continue;
 
     theGemDigiSimLinks_.push_back(GEMDigiSimLink(strip, bx, hit->particleType(), hit->trackId(), hit->eventId()));
   }
 }

void GEMDigiModule::fillDigis	(	int	rollDetId,
		GEMDigiCollection &	digis
	)

Definition at line 37 of file GEMDigiModule.cc.

References addLinks(), addLinksWithPartId(), ztail::d, and strips_.

                                                                      {
   for (const auto& d : strips_) {
     if (d.second == -999)
       continue;
     // (strip, bx)
     GEMDigi digi(d.first, d.second);
     digis.insertDigi(GEMDetId(rollDetId), digi);
     addLinks(d.first, d.second);
     addLinksWithPartId(d.first, d.second);
   }
   strips_.clear();
 }

const GEMDigiSimLinks& GEMDigiModule::gemDigiSimLinks ( ) const

inline

Definition at line 51 of file GEMDigiModule.h.

References theGemDigiSimLinks_.

51 { return theGemDigiSimLinks_; }

GEMDigiModule::theGemDigiSimLinks_

GEMDigiSimLinks theGemDigiSimLinks_

Definition: GEMDigiModule.h:65

void GEMDigiModule::setGeometry ( const GEMGeometry * geom )

Definition at line 98 of file GEMDigiModule.cc.

References isotrackApplyRegressor::model, and models.

                                                        {
   for (auto&& model : models) {
     model->setGeometry(geom);
   }
 }

void GEMDigiModule::simulate	(	const GEMEtaPartition *	roll,
		const edm::PSimHitContainer &	simHits,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 23 of file GEMDigiModule.cc.

References edm::DetSet< T >::clear(), detectorHitMap_, GEMEtaPartition::id(), isotrackApplyRegressor::model, models, DetId::rawId(), stripDigiSimLinks_, strips_, and theGemDigiSimLinks_.

                                                            {
   stripDigiSimLinks_.clear();
   detectorHitMap_.clear();
   stripDigiSimLinks_ = StripDigiSimLinks(roll->id().rawId());
   theGemDigiSimLinks_.clear();
   theGemDigiSimLinks_ = GEMDigiSimLinks(roll->id().rawId());
   for (auto&& model : models) {
     model->simulate(roll, simHits, engine, strips_, detectorHitMap_);
   }
   return;
 }