#include <GEMSim.h>

Inheritance diagram for GEMSim:

Public Types
typedef edm::DetSet < StripDigiSimLink >	StripDigiSimLinks

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

GEMSimSetUp *	getGEMSimSetUp ()

void	setGEMSimSetUp (GEMSimSetUp *setup)

virtual void	setRandomEngine (CLHEP::HepRandomEngine &eng)=0

virtual void	simulate (const GEMEtaPartition *roll, const edm::PSimHitContainer &rpcHits)=0

virtual void	simulateNoise (const GEMEtaPartition *roll)=0

const StripDigiSimLinks &	stripDigiSimLinks () const

virtual	~GEMSim ()

Protected Types
typedef std::multimap < std::pair< unsigned int, int > , const PSimHit *, std::less < std::pair< unsigned int, int > > >	DetectorHitMap

Protected Member Functions
virtual void	addLinks (unsigned int strip, int bx)

	GEMSim (const edm::ParameterSet &config)

virtual void	init ()=0

Protected Attributes
DetectorHitMap	detectorHitMap_

GEMSimSetUp *	simSetUp_

StripDigiSimLinks	stripDigiSimLinks_

std::set< std::pair< int, int > >	strips_

Detailed Description

Base Class for the GEM strip response simulation

Author: Vadim Khotilovich

Definition at line 26 of file GEMSim.h.

Member Typedef Documentation

typedef std::multimap< std::pair<unsigned int, int>, const PSimHit*, std::less<std::pair<unsigned int, int> > > GEMSim::DetectorHitMap

protected

Definition at line 65 of file GEMSim.h.

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

Definition at line 30 of file GEMSim.h.

Constructor & Destructor Documentation

virtual GEMSim::~GEMSim ( )

inlinevirtual

Definition at line 32 of file GEMSim.h.

32 {}

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

inlineprotected

Definition at line 50 of file GEMSim.h.

50 {}

Member Function Documentation

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

protectedvirtual

creates links from Digi to SimTrack disabled for now

Definition at line 21 of file GEMSim.cc.

References DeDxDiscriminatorTools::charge(), detectorHitMap_, PSimHit::eventId(), 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(auto &charge: simTrackChargeMap)
   {
     int simTrackId = charge.first;
     stripDigiSimLinks_.push_back( 
       StripDigiSimLink(strip, simTrackId, eventIdMap[simTrackId], charge.second/totalCharge ));
   }
 }

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

virtual

Definition at line 6 of file GEMSim.cc.

References addLinks(), and strips_.

Referenced by GEMDigitizer::digitize().

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