GEMDigiModule Class Reference

#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

GEMDigiSimLinks

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

Definition at line 42 of file GEMDigiModule.h.

StripDigiSimLinks

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

Definition at line 41 of file GEMDigiModule.h.

Constructor & Destructor Documentation

GEMDigiModule()

GEMDigiModule::GEMDigiModule

(

const edm::ParameterSet &

config

)

Definition at line 9 of file GEMDigiModule.cc.

                                                          {
  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::~GEMDigiModule ( )

default

Member Function Documentation

addLinks()

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

private

creates links from Digi to SimTrack

Definition at line 50 of file GEMDigiModule.cc.

References nano_mu_digi_cff::bx, ALCARECOTkAlJpsiMuMu_cff::charge, detectorHitMap_, f, edm::DetSet< T >::push_back(), nano_mu_digi_cff::strip, and stripDigiSimLinks_.

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

addLinksWithPartId()

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

private

Definition at line 84 of file GEMDigiModule.cc.

References nano_mu_digi_cff::bx, detectorHitMap_, edm::DetSet< T >::push_back(), nano_mu_digi_cff::strip, and theGemDigiSimLinks_.

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()));
  }
}

fillDigis()

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();
}

gemDigiSimLinks()

const GEMDigiSimLinks& GEMDigiModule::gemDigiSimLinks ( ) const

inline

Definition at line 51 of file GEMDigiModule.h.

References theGemDigiSimLinks_.

{ return theGemDigiSimLinks_; }

setGeometry()

void GEMDigiModule::setGeometry

(

const GEMGeometry *

geom

)

Definition at line 98 of file GEMDigiModule.cc.

References relativeConstraints::geom, and isotrackApplyRegressor::model.

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

simulate()

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_, isotrackApplyRegressor::model, nano_mu_digi_cff::roll, FastTrackerRecHitCombiner_cfi::simHits, 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;
}

stripDigiSimLinks()

const StripDigiSimLinks& GEMDigiModule::stripDigiSimLinks ( ) const

inline

Definition at line 50 of file GEMDigiModule.h.

References stripDigiSimLinks_.

{ return stripDigiSimLinks_; }

Member Data Documentation

detectorHitMap_

DetectorHitMap GEMDigiModule::detectorHitMap_

private

Definition at line 63 of file GEMDigiModule.h.

Referenced by addLinks(), addLinksWithPartId(), and simulate().

geometry_

const GEMGeometry* GEMDigiModule::geometry_

private

Definition at line 54 of file GEMDigiModule.h.

models

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

private

Definition at line 56 of file GEMDigiModule.h.

stripDigiSimLinks_

StripDigiSimLinks GEMDigiModule::stripDigiSimLinks_

private

Definition at line 64 of file GEMDigiModule.h.

Referenced by addLinks(), simulate(), and stripDigiSimLinks().

strips_

Strips GEMDigiModule::strips_

private

Definition at line 62 of file GEMDigiModule.h.

Referenced by fillDigis(), and simulate().

theGemDigiSimLinks_

GEMDigiSimLinks GEMDigiModule::theGemDigiSimLinks_

private

Definition at line 65 of file GEMDigiModule.h.

Referenced by addLinksWithPartId(), gemDigiSimLinks(), and simulate().