#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 ( ) [inline, virtual]

Definition at line 32 of file GEMSim.h.

{}

GEMSim::GEMSim ( const edm::ParameterSet & config ) [inline, protected]

Definition at line 50 of file GEMSim.h.

{}

Member Function Documentation

void GEMSim::addLinks	(	unsigned int	strip,
		int	bx
	)		`[protected, virtual]`

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