#include <CSCHitAssociator.h>

Public Types
typedef edm::DetSetVector< StripDigiSimLink >	DigiSimLinks

typedef edm::DetSet< StripDigiSimLink >	LayerLinks

typedef std::pair< uint32_t, EncodedEventId >	SimHitIdpr

typedef edm::DetSetVector< StripDigiSimLink >	WireDigiSimLinks

Public Member Functions
std::vector< SimHitIdpr >	associateCSCHitId (const CSCRecHit2D *) const

std::vector< SimHitIdpr >	associateHitId (const TrackingRecHit &) const

	CSCHitAssociator (const edm::Event &, const edm::EventSetup &, const edm::ParameterSet &)

	CSCHitAssociator (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

void	initEvent (const edm::Event &, const edm::EventSetup &)

Private Attributes
const CSCGeometry *	cscgeom

edm::InputTag	linksTag

const DigiSimLinks *	theDigiSimLinks

Detailed Description

Definition at line 24 of file CSCHitAssociator.h.

Member Typedef Documentation

typedef edm::DetSetVector<StripDigiSimLink> CSCHitAssociator::DigiSimLinks

Definition at line 27 of file CSCHitAssociator.h.

typedef edm::DetSet<StripDigiSimLink> CSCHitAssociator::LayerLinks

Definition at line 29 of file CSCHitAssociator.h.

typedef std::pair<uint32_t, EncodedEventId> CSCHitAssociator::SimHitIdpr

Definition at line 30 of file CSCHitAssociator.h.

typedef edm::DetSetVector<StripDigiSimLink> CSCHitAssociator::WireDigiSimLinks

Definition at line 28 of file CSCHitAssociator.h.

Constructor & Destructor Documentation

CSCHitAssociator::CSCHitAssociator	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const edm::ParameterSet &	conf
	)

Definition at line 8 of file CSCHitAssociator.cc.

References initEvent().

                                                                                                                 : 
   theDigiSimLinks(nullptr),
   linksTag(conf.getParameter<edm::InputTag>("CSClinksTag"))
 {
   initEvent(event,setup);
 }

CSCHitAssociator::CSCHitAssociator	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 15 of file CSCHitAssociator.cc.

References linksTag.

                                                                                            : 
   theDigiSimLinks(nullptr),
   linksTag(conf.getParameter<edm::InputTag>("CSClinksTag"))
 {
   iC.consumes<DigiSimLinks>(linksTag);
 }

Member Function Documentation

std::vector< CSCHitAssociator::SimHitIdpr > CSCHitAssociator::associateCSCHitId ( const CSCRecHit2D * cscrechit ) const

Definition at line 35 of file CSCHitAssociator.cc.

References CSCLayerGeometry::channel(), CSCRecHit2D::channels(), CSCRecHit2D::cscDetId(), cscgeom, edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), geometry, CSCGeometry::layer(), LogTrace, CSCRecHit2D::nStrips(), DetId::rawId(), and theDigiSimLinks.

Referenced by MuonAssociatorByHitsHelper::getMatchedIds().

                                                                                                              {
   std::vector<SimHitIdpr> simtrackids;
   
   unsigned int detId = cscrechit->geographicalId().rawId();
   int nchannels = cscrechit->nStrips();
   const CSCLayerGeometry * laygeom = cscgeom->layer(cscrechit->cscDetId())->geometry();
   
   DigiSimLinks::const_iterator layerLinks = theDigiSimLinks->find(detId);    
   
   if (layerLinks != theDigiSimLinks->end()) {
     
     for(int idigi = 0; idigi < nchannels; ++idigi) {
       // strip and readout channel numbers may differ in ME1/1A
       int istrip = cscrechit->channels(idigi);
       int channel = laygeom->channel(istrip);
       
       for (LayerLinks::const_iterator link=layerLinks->begin(); link!=layerLinks->end(); ++link) {
     int ch = static_cast<int>(link->channel());
     if (ch == channel) {
       SimHitIdpr currentId(link->SimTrackId(), link->eventId());
       if (find(simtrackids.begin(), simtrackids.end(), currentId) == simtrackids.end())
         simtrackids.push_back(currentId);
     }
       }
     }
     
   } else LogTrace("CSCHitAssociator")
     <<"*** WARNING in CSCHitAssociator::associateCSCHitId - CSC layer "<<detId<<" has no DigiSimLinks !"<<std::endl;   
   
   return simtrackids;
 }

std::vector< CSCHitAssociator::SimHitIdpr > CSCHitAssociator::associateHitId ( const TrackingRecHit & hit ) const

Definition at line 68 of file CSCHitAssociator.cc.

References CSCLayerGeometry::channel(), CSCRecHit2D::channels(), CSCRecHit2D::cscDetId(), cscgeom, edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), geometry, CSCGeometry::layer(), LogTrace, CSCRecHit2D::nStrips(), DetId::rawId(), and theDigiSimLinks.

 {
   std::vector<SimHitIdpr> simtrackids;
   
   const TrackingRecHit * hitp = &hit;
   const CSCRecHit2D * cscrechit = dynamic_cast<const CSCRecHit2D *>(hitp);
 
   if (cscrechit) {
     
     unsigned int detId = cscrechit->geographicalId().rawId();
     int nchannels = cscrechit->nStrips();
     const CSCLayerGeometry * laygeom = cscgeom->layer(cscrechit->cscDetId())->geometry();
 
     DigiSimLinks::const_iterator layerLinks = theDigiSimLinks->find(detId);    
 
     if (layerLinks != theDigiSimLinks->end()) {
       
       for(int idigi = 0; idigi < nchannels; ++idigi) {
     // strip and readout channel numbers may differ in ME1/1A
     int istrip = cscrechit->channels(idigi);
     int channel = laygeom->channel(istrip);
     
     for (LayerLinks::const_iterator link=layerLinks->begin(); link!=layerLinks->end(); ++link) {
       int ch = static_cast<int>(link->channel());
       if (ch == channel) {
         SimHitIdpr currentId(link->SimTrackId(), link->eventId());
         if (find(simtrackids.begin(), simtrackids.end(), currentId) == simtrackids.end())
           simtrackids.push_back(currentId);
       }
     }
       }
       
     } else LogTrace("CSCHitAssociator")
       <<"*** WARNING in CSCHitAssociator::associateHitId - CSC layer "<<detId<<" has no DigiSimLinks !"<<std::endl;   
     
   } else LogTrace("CSCHitAssociator")<<"*** WARNING in CSCHitAssociator::associateHitId, null dynamic_cast !";
   
   return simtrackids;
 }

void CSCHitAssociator::initEvent	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

Definition at line 22 of file CSCHitAssociator.cc.

References cscgeom, TTClusterAssociation_cfi::digiSimLinks, edm::EventSetup::get(), linksTag, LogTrace, edm::Handle< T >::product(), and theDigiSimLinks.

Referenced by CSCHitAssociator().

                                                                                   {
 
   edm::Handle<DigiSimLinks> digiSimLinks;
   LogTrace("CSCHitAssociator") <<"getting CSC Strip DigiSimLink collection - "<<linksTag;
   event.getByLabel(linksTag, digiSimLinks);
   theDigiSimLinks = digiSimLinks.product();
 
   // get CSC Geometry to use CSCLayer methods
   edm::ESHandle<CSCGeometry> mugeom;
   setup.get<MuonGeometryRecord>().get( mugeom );
   cscgeom = &*mugeom;
 }