#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 GEMRecHit *gemrechit) const

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

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

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 33 of file GEMHitAssociator.h.

Member Typedef Documentation

◆ DigiSimLinks

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

Definition at line 35 of file GEMHitAssociator.h.

◆ LayerLinks

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

Definition at line 36 of file GEMHitAssociator.h.

◆ SimHitIdpr

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

Definition at line 37 of file GEMHitAssociator.h.

Constructor & Destructor Documentation

◆ GEMHitAssociator() [1/2]

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

Definition at line 6 of file GEMHitAssociator.cc.

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

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

◆ GEMHitAssociator() [2/2]

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

Definition at line 22 of file GEMHitAssociator.cc.

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

References MillePedeFileConverter_cfg::e, and initEvent().

◆ ~GEMHitAssociator()

GEMHitAssociator::~GEMHitAssociator ( )

inline

Definition at line 46 of file GEMHitAssociator.h.

46 {}

Member Function Documentation

◆ associateRecHit()

std::vector< GEMHitAssociator::SimHitIdpr > GEMHitAssociator::associateRecHit ( const GEMRecHit * gemrechit ) const

Definition at line 70 of file GEMHitAssociator.cc.

                                                                                                         {
   std::vector<SimHitIdpr> matched;
  
   if (useGEMs_) {
     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;
 }

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

Referenced by MuonAssociatorByHitsHelper::getMatchedIds().

◆ initEvent()

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

Definition at line 34 of file GEMHitAssociator.cc.

                                                                                    {
   if (useGEMs_) {
     if (crossingframe) {
       edm::Handle<CrossingFrame<PSimHit>> cf;
       LogTrace("GEMHitAssociator") << "getting CrossingFrame<PSimHit> collection - " << GEMsimhitsXFTag;
       e.getByLabel(GEMsimhitsXFTag, cf);
  
       std::unique_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();
   }
 }