GEMDigiModel Class Reference

#include <GEMDigiModel.h>

Inheritance diagram for GEMDigiModel:

Public Types
typedef edm::DetSet< GEMDigiSimLink >	GEMDigiSimLinks
typedef edm::DetSet< StripDigiSimLink >	StripDigiSimLinks

Public Member Functions
void	fillDigis (int rollDetId, GEMDigiCollection &)
const GEMDigiSimLinks &	gemDigiSimLinks () const
const GEMGeometry *	getGeometry ()
void	setGeometry (const GEMGeometry *geom)
virtual void	setup ()=0
virtual std::vector< std::pair< int, int > >	simulateClustering (const GEMEtaPartition *, const PSimHit *, const int, CLHEP::HepRandomEngine *engine)=0
virtual void	simulateNoise (const GEMEtaPartition *, CLHEP::HepRandomEngine *engine)=0
virtual void	simulateSignal (const GEMEtaPartition *, const edm::PSimHitContainer &, CLHEP::HepRandomEngine *engine)=0
const StripDigiSimLinks &	stripDigiSimLinks () const
virtual	~GEMDigiModel ()

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

Protected Member Functions
void	addLinks (unsigned int strip, int bx)
	creates links from Digi to SimTrack More...
void	addLinksWithPartId (unsigned int strip, int bx)
	GEMDigiModel (const edm::ParameterSet &)

Protected Attributes
DetectorHitMap	detectorHitMap_
const GEMGeometry *	geometry_
StripDigiSimLinks	stripDigiSimLinks_
std::set< std::pair< int, int > >	strips_
GEMDigiSimLinks	theGemDigiSimLinks_

Detailed Description

Base Class for the GEM strip response simulation

Author

Sven Dildick

Definition at line 32 of file GEMDigiModel.h.

Member Typedef Documentation

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

protected

Definition at line 76 of file GEMDigiModel.h.

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

Definition at line 37 of file GEMDigiModel.h.

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

Definition at line 36 of file GEMDigiModel.h.

Constructor & Destructor Documentation

virtual GEMDigiModel::~GEMDigiModel ( )

inlinevirtual

Definition at line 39 of file GEMDigiModel.h.

{}

GEMDigiModel::GEMDigiModel ( const edm::ParameterSet &  )

inlineprotected

Definition at line 61 of file GEMDigiModel.h.

{}

Member Function Documentation

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

protected

creates links from Digi to SimTrack

Definition at line 21 of file GEMDigiModel.cc.

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

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

protected

Definition at line 60 of file GEMDigiModel.cc.

References detectorHitMap_, PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::eventId(), PSimHit::momentumAtEntry(), PSimHit::particleType(), PSimHit::processType(), edm::DetSet< T >::push_back(), theGemDigiSimLinks_, PSimHit::timeOfFlight(), 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(digi, hit->entryPoint(), hit->momentumAtEntry(), hit->timeOfFlight(), hit->energyLoss(),
                                                       hit->particleType(), hit->detUnitId(), hit->trackId(), hit->eventId(), hit->processType()));

  }
}

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

Definition at line 5 of file GEMDigiModel.cc.

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

Referenced by GEMDigiProducer::produce().

{
  for (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& GEMDigiModel::gemDigiSimLinks ( ) const

inline

Definition at line 57 of file GEMDigiModel.h.

Referenced by GEMDigiProducer::produce().

{return theGemDigiSimLinks_;}

const GEMGeometry* GEMDigiModel::getGeometry ( )

inline

Definition at line 43 of file GEMDigiModel.h.

References hcalnoiseinfoproducer_cfi::fillDigis, and GeneralSetup::setup().

Referenced by GEMDigiProducer::produce().

{return geometry_;}

void GEMDigiModel::setGeometry ( const GEMGeometry *  geom )

inline

Definition at line 41 of file GEMDigiModel.h.

References relativeConstraints::geom.

Referenced by GEMDigiProducer::beginRun().

{geometry_ = geom;}

virtual void GEMDigiModel::setup ( )

pure virtual

Implemented in GEMTrivialModel, and GEMSimpleModel.

Referenced by GEMDigiProducer::beginRun().

virtual std::vector<std::pair<int,int> > GEMDigiModel::simulateClustering ( const GEMEtaPartition *  , const PSimHit *  , const int  , CLHEP::HepRandomEngine *  engine  )

pure virtual

Implemented in GEMSimpleModel, and GEMTrivialModel.

virtual void GEMDigiModel::simulateNoise ( const GEMEtaPartition *  , CLHEP::HepRandomEngine *  engine  )

pure virtual

Implemented in GEMSimpleModel, and GEMTrivialModel.

Referenced by GEMDigiProducer::produce().

virtual void GEMDigiModel::simulateSignal ( const GEMEtaPartition *  , const edm::PSimHitContainer &  , CLHEP::HepRandomEngine *  engine  )

pure virtual

Implemented in GEMSimpleModel, and GEMTrivialModel.

Referenced by GEMDigiProducer::produce().

const StripDigiSimLinks& GEMDigiModel::stripDigiSimLinks ( ) const

inline

Definition at line 56 of file GEMDigiModel.h.

Referenced by GEMDigiProducer::produce().

{return stripDigiSimLinks_;}

Member Data Documentation

DetectorHitMap GEMDigiModel::detectorHitMap_

protected

Definition at line 78 of file GEMDigiModel.h.

Referenced by addLinks(), addLinksWithPartId(), GEMTrivialModel::simulateSignal(), and GEMSimpleModel::simulateSignal().

const GEMGeometry* GEMDigiModel::geometry_

protected

Definition at line 63 of file GEMDigiModel.h.

Referenced by GEMSimpleModel::getSimHitBx().

StripDigiSimLinks GEMDigiModel::stripDigiSimLinks_

protected

Definition at line 79 of file GEMDigiModel.h.

Referenced by addLinks(), GEMTrivialModel::simulateSignal(), and GEMSimpleModel::simulateSignal().

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

protected

Definition at line 65 of file GEMDigiModel.h.

Referenced by fillDigis(), GEMSimpleModel::simulateNoise(), GEMTrivialModel::simulateSignal(), and GEMSimpleModel::simulateSignal().

GEMDigiSimLinks GEMDigiModel::theGemDigiSimLinks_

protected

Definition at line 80 of file GEMDigiModel.h.

Referenced by addLinksWithPartId(), and GEMSimpleModel::simulateSignal().