#include <GEMHitAssociator.h>

Public Types
typedef edm::DetSetVector < StripDigiSimLink >	DigiSimLinks

typedef edm::DetSet < StripDigiSimLink >	LayerLinks

typedef std::pair< uint32_t, EncodedEventId >	SimHitIdpr

Public Member Functions
std::vector< SimHitIdpr >	associateRecHit (const TrackingRecHit &hit) const

	GEMHitAssociator (const edm::ParameterSet &, edm::ConsumesCollector &&ic)

	GEMHitAssociator (const edm::Event &e, const edm::EventSetup &eventSetup, const edm::ParameterSet &conf)

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

	~GEMHitAssociator ()

Private Attributes
std::map< unsigned int, edm::PSimHitContainer >	_SimHitMap

bool	crossingframe

edm::InputTag	GEMdigisimlinkTag

edm::EDGetTokenT < edm::DetSetVector < StripDigiSimLink > >	GEMdigisimlinkToken_

edm::InputTag	GEMsimhitsTag

edm::EDGetTokenT < edm::PSimHitContainer >	GEMsimhitsToken_

edm::InputTag	GEMsimhitsXFTag

edm::EDGetTokenT < CrossingFrame< PSimHit > >	GEMsimhitsXFToken_

const DigiSimLinks *	theDigiSimLinks

bool	useGEMs_

Detailed Description

Definition at line 32 of file GEMHitAssociator.h.

Member Typedef Documentation

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

Definition at line 36 of file GEMHitAssociator.h.

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

Definition at line 37 of file GEMHitAssociator.h.

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

Definition at line 38 of file GEMHitAssociator.h.

Constructor & Destructor Documentation

GEMHitAssociator::GEMHitAssociator	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	ic
	)

Definition at line 6 of file GEMHitAssociator.cc.

References crossingframe, GEMdigisimlinkTag, GEMdigisimlinkToken_, GEMsimhitsTag, GEMsimhitsToken_, GEMsimhitsXFTag, GEMsimhitsXFToken_, and edm::InputTag::label().

                                                :
   GEMdigisimlinkTag(conf.getParameter<edm::InputTag>("GEMdigisimlinkTag")),
   // CrossingFrame used or not ?
   crossingframe(conf.getParameter<bool>("crossingframe")),
   useGEMs_(conf.getParameter<bool>("useGEMs")),
   GEMsimhitsTag(conf.getParameter<edm::InputTag>("GEMsimhitsTag")),
   GEMsimhitsXFTag(conf.getParameter<edm::InputTag>("GEMsimhitsXFTag"))
 {
   if (crossingframe){
     GEMsimhitsXFToken_=iC.consumes<CrossingFrame<PSimHit> >(GEMsimhitsXFTag);
   } else if (!GEMsimhitsTag.label().empty()) {
     GEMsimhitsToken_=iC.consumes<edm::PSimHitContainer>(GEMsimhitsTag);
   }
 
   GEMdigisimlinkToken_=iC.consumes< edm::DetSetVector<StripDigiSimLink> >(GEMdigisimlinkTag);
 }

GEMHitAssociator::GEMHitAssociator	(	const edm::Event &	e,
		const edm::EventSetup &	eventSetup,
		const edm::ParameterSet &	conf
	)

Definition at line 24 of file GEMHitAssociator.cc.

References initEvent().

                                                                                                                   :
   GEMdigisimlinkTag(conf.getParameter<edm::InputTag>("GEMdigisimlinkTag")),
   // CrossingFrame used or not ?
   crossingframe(conf.getParameter<bool>("crossingframe")),
   useGEMs_(conf.getParameter<bool>("useGEMs")),
   GEMsimhitsTag(conf.getParameter<edm::InputTag>("GEMsimhitsTag")),
   GEMsimhitsXFTag(conf.getParameter<edm::InputTag>("GEMsimhitsXFTag"))
 {
   initEvent(e,eventSetup);
 }

GEMHitAssociator::~GEMHitAssociator ( )

inline

Definition at line 47 of file GEMHitAssociator.h.

47 {}

Member Function Documentation

std::vector< GEMHitAssociator::SimHitIdpr > GEMHitAssociator::associateRecHit ( const TrackingRecHit & hit ) const

Definition at line 82 of file GEMHitAssociator.cc.

References GEMRecHit::clusterSize(), edm::DetSetVector< T >::end(), spr::find(), edm::DetSetVector< T >::find(), GEMRecHit::firstClusterStrip(), GEMRecHit::gemId(), i, theDigiSimLinks, and useGEMs_.

Referenced by MuonAssociatorByHitsHelper::getMatchedIds().

                                                                                                         {
   
   std::vector<SimHitIdpr> matched;
 
   if(useGEMs_){
     
       //std::cout<<"gemboo "<<useGEMs_<<std::endl;
 
       const TrackingRecHit * hitp = &hit;
       const GEMRecHit * gemrechit = dynamic_cast<const GEMRecHit *>(hitp);
 
       if (gemrechit) {
         
         GEMDetId gemDetId = gemrechit->gemId();
         int fstrip = gemrechit->firstClusterStrip();
         int cls = gemrechit->clusterSize();
         //int bx = gemrechit->BunchX();
 
         DigiSimLinks::const_iterator layerLinks = theDigiSimLinks->find(gemDetId);
 
         if (layerLinks != theDigiSimLinks->end()) {
         
         for(int i = fstrip; i < (fstrip+cls); ++i) {
                   
             for(LayerLinks::const_iterator itlink = layerLinks->begin(); itlink != layerLinks->end(); ++itlink) {
 
                 int ch = static_cast<int>(itlink->channel());
                 if(ch != i) continue;
 
                 SimHitIdpr currentId(itlink->SimTrackId(), itlink->eventId());
                 if(find(matched.begin(),matched.end(),currentId ) == matched.end())
                     matched.push_back(currentId);
 
             }
 
         }
 
         }else edm::LogWarning("GEMHitAssociator")
           <<"*** WARNING in GEMHitAssociator: GEM layer "<<gemDetId<<" has no DigiSimLinks !"<<std::endl;
         
       } else edm::LogWarning("GEMHitAssociator")<<"*** WARNING in GEMHitAssociator::associateRecHit, null dynamic_cast !";
 
   }
   
   return  matched;
 
 }

void GEMHitAssociator::initEvent	(	const edm::Event &	e,
		const edm::EventSetup &	eventSetup
	)

Definition at line 35 of file GEMHitAssociator.cc.

References _SimHitMap, crossingframe, GEMdigisimlinkTag, GEMsimhitsTag, GEMsimhitsXFTag, edm::Event::getByLabel(), edm::InputTag::label(), LogTrace, edm::Handle< T >::product(), theDigiSimLinks, and useGEMs_.

Referenced by GEMHitAssociator().

 {
 
   if(useGEMs_){
 
       if (crossingframe) {
         
         edm::Handle<CrossingFrame<PSimHit> > cf;
         LogTrace("GEMHitAssociator") <<"getting CrossingFrame<PSimHit> collection - "<<GEMsimhitsXFTag;
         e.getByLabel(GEMsimhitsXFTag, cf);
         
         std::auto_ptr<MixCollection<PSimHit> > 
           GEMsimhits( new MixCollection<PSimHit>(cf.product()) );
         LogTrace("GEMHitAssociator") <<"... size = "<<GEMsimhits->size();
 
         //   MixCollection<PSimHit> & simHits = *hits;
         
         for(MixCollection<PSimHit>::MixItr hitItr = GEMsimhits->begin();
         hitItr != GEMsimhits->end(); ++hitItr) 
           {
         _SimHitMap[hitItr->detUnitId()].push_back(*hitItr);
           }
         
       } else if (!GEMsimhitsTag.label().empty()) {
         edm::Handle<edm::PSimHitContainer> GEMsimhits;
         LogTrace("GEMHitAssociator") <<"getting PSimHit collection - "<<GEMsimhitsTag;
         e.getByLabel(GEMsimhitsTag, GEMsimhits);    
         LogTrace("GEMHitAssociator") <<"... size = "<<GEMsimhits->size();
         
         // arrange the hits by detUnit
         for(edm::PSimHitContainer::const_iterator hitItr = GEMsimhits->begin();
         hitItr != GEMsimhits->end(); ++hitItr)
           {
         _SimHitMap[hitItr->detUnitId()].push_back(*hitItr);
           }
       }
 
       edm::Handle<DigiSimLinks> digiSimLinks;
       LogTrace("GEMHitAssociator") <<"getting GEM Strip DigiSimLink collection - "<<GEMdigisimlinkTag;
       e.getByLabel(GEMdigisimlinkTag, digiSimLinks);
       theDigiSimLinks = digiSimLinks.product();
 
   }
 
 }