#include <GEMRecHitMatcher.h>

Public Member Functions
bool	areGEMRecHitSame (const GEMRecHit &l, const GEMRecHit &r) const

std::set< unsigned int >	chamberIds () const

std::set< unsigned int >	detIds () const

std::shared_ptr< GEMDigiMatcher >	gemDigiMatcher () const

	GEMRecHitMatcher (edm::ParameterSet const &iPS, edm::ConsumesCollector &&iC)

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

bool	isGEMRecHitMatched (const GEMRecHit &thisRh) const

void	match (const SimTrack &t, const SimVertex &v)
	do the matching More...

int	nGEMRecHits () const
	How many recHits in GEM did this simtrack get in total? More...

int	nLayersWithRecHitsInSuperChamber (unsigned int) const

std::set< int >	partitionNumbers () const

bool	recHitInContainer (const GEMRecHit &rh, const GEMRecHitContainer &c) const

GlobalPoint	recHitMeanPosition (const GEMRecHitContainer &rechits) const

GlobalPoint	recHitPosition (const GEMRecHit &rechit) const

const GEMRecHitContainer &	recHits () const

const GEMRecHitContainer &	recHitsInChamber (unsigned int) const

const GEMRecHitContainer &	recHitsInDetId (unsigned int) const

const GEMRecHitContainer &	recHitsInSuperChamber (unsigned int) const

std::set< int >	stripNumbersInDetId (unsigned int) const

std::set< unsigned int >	superChamberIds () const

bool	verbose () const

	~GEMRecHitMatcher ()

Private Member Functions
void	matchRecHitsToSimTrack (const GEMRecHitCollection &recHits)

Private Attributes
std::map< unsigned int, GEMRecHitContainer >	chamber_to_recHits_

std::map< unsigned int, GEMRecHitContainer >	detid_to_recHits_

edm::ESHandle< GEMGeometry >	gem_geom_

std::shared_ptr< GEMDigiMatcher >	gemDigiMatcher_

const GEMGeometry *	gemGeometry_

edm::Handle< GEMRecHitCollection >	gemRecHitH_

edm::EDGetTokenT< GEMRecHitCollection >	gemRecHitToken_

int	maxBX_

int	minBX_

const GEMRecHitContainer	no_recHits_

GEMRecHitContainer	recHits_

std::map< unsigned int, GEMRecHitContainer >	superchamber_to_recHits_

bool	verbose_

Detailed Description

Description: Matching of RecHits for SimTrack in GEM

Original Author : "Vadim Khotilovich" Contibuting Author : "Claudio Caputo"

Definition at line 31 of file GEMRecHitMatcher.h.

Constructor & Destructor Documentation

◆ GEMRecHitMatcher()

GEMRecHitMatcher::GEMRecHitMatcher	(	edm::ParameterSet const &	iPS,
		edm::ConsumesCollector &&	iC
	)

Definition at line 7 of file GEMRecHitMatcher.cc.

                                                                                          {
   const auto& gemRecHit = pset.getParameterSet("gemRecHit");
   minBX_ = gemRecHit.getParameter<int>("minBX");
   maxBX_ = gemRecHit.getParameter<int>("maxBX");
   verbose_ = gemRecHit.getParameter<int>("verbose");
  
   // make a new digi matcher
   gemDigiMatcher_.reset(new GEMDigiMatcher(pset, std::move(iC)));
  
   gemRecHitToken_ = iC.consumes<GEMRecHitCollection>(gemRecHit.getParameter<edm::InputTag>("inputTag"));
 }

References muonGEMRecHitPSet::gemRecHit, eostools::move(), and muonDTDigis_cfi::pset.

◆ ~GEMRecHitMatcher()

GEMRecHitMatcher::~GEMRecHitMatcher ( )

inline

Definition at line 37 of file GEMRecHitMatcher.h.

37 {}

Member Function Documentation

◆ areGEMRecHitSame()

bool GEMRecHitMatcher::areGEMRecHitSame	(	const GEMRecHit &	l,
		const GEMRecHit &	r
	)		const

Definition at line 202 of file GEMRecHitMatcher.cc.

                                                                                     {
   return l.localPosition() == r.localPosition() and l.BunchX() == r.BunchX();
 }

References cmsLHEtoEOSManager::l, and alignCSCRings::r.

◆ chamberIds()

std::set< unsigned int > GEMRecHitMatcher::chamberIds ( ) const

Definition at line 100 of file GEMRecHitMatcher.cc.

                                                         {
   std::set<unsigned int> result;
   for (const auto& p : chamber_to_recHits_)
     result.insert(p.first);
   return result;
 }

References AlCaHLTBitMon_ParallelJobs::p, and mps_fire::result.

◆ detIds()

std::set< unsigned int > GEMRecHitMatcher::detIds ( void ) const

Definition at line 93 of file GEMRecHitMatcher.cc.

                                                     {
   std::set<unsigned int> result;
   for (const auto& p : detid_to_recHits_)
     result.insert(p.first);
   return result;
 }

References AlCaHLTBitMon_ParallelJobs::p, and mps_fire::result.

◆ gemDigiMatcher()

std::shared_ptr<GEMDigiMatcher> GEMRecHitMatcher::gemDigiMatcher ( ) const

inline

Definition at line 80 of file GEMRecHitMatcher.h.

80 { return gemDigiMatcher_; }

References gemDigiMatcher_.

◆ init()

void GEMRecHitMatcher::init	(	const edm::Event &	e,
		const edm::EventSetup &	eventSetup
	)

Definition at line 19 of file GEMRecHitMatcher.cc.

                                                                                {
   gemDigiMatcher_->init(iEvent, iSetup);
  
   iEvent.getByToken(gemRecHitToken_, gemRecHitH_);
  
   iSetup.get<MuonGeometryRecord>().get(gem_geom_);
   if (gem_geom_.isValid()) {
     gemGeometry_ = &*gem_geom_;
   } else {
     std::cout << "+++ Info: GEM geometry is unavailable. +++\n";
   }
 }

References gather_cfg::cout, edm::EventSetup::get(), get, and iEvent.

◆ isGEMRecHitMatched()

bool GEMRecHitMatcher::isGEMRecHitMatched ( const GEMRecHit & thisRh ) const

Definition at line 198 of file GEMRecHitMatcher.cc.

                                                                        {
   return recHitInContainer(thisRh, recHits());
 }

References FastTrackerRecHitMaskProducer_cfi::recHits.

◆ match()

void GEMRecHitMatcher::match	(	const SimTrack &	t,
		const SimVertex &	v
	)

do the matching

Definition at line 33 of file GEMRecHitMatcher.cc.

                                                                   {
   // match digis first
   gemDigiMatcher_->match(t, v);
  
   // get the rechit collection
   const GEMRecHitCollection& gemRecHits = *gemRecHitH_.product();
  
   // now match the rechits
   matchRecHitsToSimTrack(gemRecHits);
 }

References gemRecHits_cfi::gemRecHits, OrderedSet::t, and findQualityFiles::v.

◆ matchRecHitsToSimTrack()

void GEMRecHitMatcher::matchRecHitsToSimTrack ( const GEMRecHitCollection & recHits )

private

Definition at line 44 of file GEMRecHitMatcher.cc.

                                                                                 {
   // get the matched ids with digis
   const auto& det_ids = gemDigiMatcher_->detIdsDigi();
  
   // loop on those ids
   for (auto id : det_ids) {
     // now check the digis in this detid
     const auto& hit_digis = gemDigiMatcher_->stripNumbersInDetId(id);
     if (verbose()) {
       cout << "hit_digis_fat ";
       copy(hit_digis.begin(), hit_digis.end(), ostream_iterator<int>(cout, " "));
       cout << endl;
     }
  
     GEMDetId p_id(id);
     const auto& rechits_in_det = rechits.get(p_id);
  
     for (auto d = rechits_in_det.first; d != rechits_in_det.second; ++d) {
       if (verbose())
         cout << "recHit " << p_id << " " << *d << endl;
  
       // check that the rechit is within BX range
       if (d->BunchX() < minBX_ || d->BunchX() > maxBX_)
         continue;
  
       int firstStrip = d->firstClusterStrip();
       int cls = d->clusterSize();
       bool stripFound = false;
  
       // check that it matches a strip that was hit by digis from our track
       for (int i = firstStrip; i < (firstStrip + cls); i++) {
         if (hit_digis.find(i) != hit_digis.end())
           stripFound = true;
       }
  
       // this rechit did not correspond with any previously matched digi
       if (!stripFound)
         continue;
       if (verbose())
         cout << "oki" << endl;
  
       recHits_.push_back(*d);
       detid_to_recHits_[id].push_back(*d);
       chamber_to_recHits_[p_id.chamberId().rawId()].push_back(*d);
       superchamber_to_recHits_[p_id.superChamberId().rawId()].push_back(*d);
     }
   }
 }

References GEMDetId::chamberId(), filterCSVwithJSON::copy, gather_cfg::cout, ztail::d, mps_fire::i, triggerObjects_cff::id, DetId::rawId(), TrackInfoProducer_cfi::rechits, GEMDetId::superChamberId(), and verbose.

◆ nGEMRecHits()

int GEMRecHitMatcher::nGEMRecHits ( ) const

How many recHits in GEM did this simtrack get in total?

◆ nLayersWithRecHitsInSuperChamber()

int GEMRecHitMatcher::nLayersWithRecHitsInSuperChamber ( unsigned int detid ) const

Definition at line 132 of file GEMRecHitMatcher.cc.

                                                                                {
   set<int> layers;
   for (const auto& d : recHitsInSuperChamber(detid)) {
     layers.insert(d.gemId().layer());
   }
   return layers.size();
 }

References ztail::d, and hgcalTopologyTester_cfi::layers.

◆ partitionNumbers()

std::set< int > GEMRecHitMatcher::partitionNumbers ( ) const

Definition at line 151 of file GEMRecHitMatcher.cc.

                                                      {
   std::set<int> result;
  
   for (auto id : detIds()) {
     GEMDetId idd(id);
     result.insert(idd.roll());
   }
   return result;
 }

References mps_fire::result, and GEMDetId::roll().

◆ recHitInContainer()

bool GEMRecHitMatcher::recHitInContainer	(	const GEMRecHit &	rh,
		const GEMRecHitContainer &	c
	)		const

Definition at line 190 of file GEMRecHitMatcher.cc.

                                                                                                {
   bool isSame = false;
   for (const auto& thisRH : c)
     if (areGEMRecHitSame(thisRH, rh))
       isSame = true;
   return isSame;
 }

References HltBtagPostValidation_cff::c.

◆ recHitMeanPosition()

GlobalPoint GEMRecHitMatcher::recHitMeanPosition ( const GEMRecHitContainer & rechits ) const

Definition at line 167 of file GEMRecHitMatcher.cc.

                                                                                        {
   GlobalPoint point_zero;
   if (rechit.empty())
     return point_zero;  // point "zero"
  
   float sumx, sumy, sumz;
   sumx = sumy = sumz = 0.f;
   size_t n = 0;
   for (const auto& d : rechit) {
     GlobalPoint gp = recHitPosition(d);
     if (gp == point_zero)
       continue;
  
     sumx += gp.x();
     sumy += gp.y();
     sumz += gp.z();
     ++n;
   }
   if (n == 0)
     return GlobalPoint();
   return GlobalPoint(sumx / n, sumy / n, sumz / n);
 }

References ztail::d, runTauDisplay::gp, and dqmiodumpmetadata::n.

◆ recHitPosition()

GlobalPoint GEMRecHitMatcher::recHitPosition ( const GEMRecHit & rechit ) const

Definition at line 161 of file GEMRecHitMatcher.cc.

                                                                           {
   const GEMDetId& idd = rechit.gemId();
   const LocalPoint& lp = rechit.localPosition();
   return gemGeometry_->idToDet(idd)->surface().toGlobal(lp);
 }

References GEMRecHit::gemId(), and GEMRecHit::localPosition().

◆ recHits()

const GEMRecHitContainer& GEMRecHitMatcher::recHits ( void ) const

inline

Definition at line 55 of file GEMRecHitMatcher.h.

55 { return recHits_; }

References recHits_.

◆ recHitsInChamber()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInChamber ( unsigned int detid ) const

Definition at line 120 of file GEMRecHitMatcher.cc.

                                                                                      {
   if (chamber_to_recHits_.find(detid) == chamber_to_recHits_.end())
     return no_recHits_;
   return chamber_to_recHits_.at(detid);
 }

◆ recHitsInDetId()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInDetId ( unsigned int detid ) const

Definition at line 114 of file GEMRecHitMatcher.cc.

                                                                                    {
   if (detid_to_recHits_.find(detid) == detid_to_recHits_.end())
     return no_recHits_;
   return detid_to_recHits_.at(detid);
 }

◆ recHitsInSuperChamber()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInSuperChamber ( unsigned int detid ) const

Definition at line 126 of file GEMRecHitMatcher.cc.

                                                                                           {
   if (superchamber_to_recHits_.find(detid) == superchamber_to_recHits_.end())
     return no_recHits_;
   return superchamber_to_recHits_.at(detid);
 }

◆ stripNumbersInDetId()

std::set< int > GEMRecHitMatcher::stripNumbersInDetId ( unsigned int detid ) const

Definition at line 140 of file GEMRecHitMatcher.cc.

                                                                           {
   set<int> result;
   for (const auto& d : recHitsInDetId(detid)) {
     // loop on all strips hit in this rechit
     for (int iStrip = d.firstClusterStrip(); iStrip < d.firstClusterStrip() + d.clusterSize(); ++iStrip) {
       result.insert(iStrip);
     }
   }
   return result;
 }

References ztail::d, and mps_fire::result.

◆ superChamberIds()

std::set< unsigned int > GEMRecHitMatcher::superChamberIds ( ) const

Definition at line 107 of file GEMRecHitMatcher.cc.

                                                              {
   std::set<unsigned int> result;
   for (const auto& p : superchamber_to_recHits_)
     result.insert(p.first);
   return result;
 }

References AlCaHLTBitMon_ParallelJobs::p, and mps_fire::result.

◆ verbose()

bool GEMRecHitMatcher::verbose ( ) const

inline

Definition at line 72 of file GEMRecHitMatcher.h.

72 { return verbose_; }

References verbose_.

Member Data Documentation

◆ chamber_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::chamber_to_recHits_

private

Definition at line 101 of file GEMRecHitMatcher.h.

◆ detid_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::detid_to_recHits_

private

Definition at line 100 of file GEMRecHitMatcher.h.

◆ gem_geom_

edm::ESHandle<GEMGeometry> GEMRecHitMatcher::gem_geom_

private

Definition at line 94 of file GEMRecHitMatcher.h.

◆ gemDigiMatcher_

std::shared_ptr<GEMDigiMatcher> GEMRecHitMatcher::gemDigiMatcher_

private

Definition at line 107 of file GEMRecHitMatcher.h.

Referenced by gemDigiMatcher().

◆ gemGeometry_

const GEMGeometry* GEMRecHitMatcher::gemGeometry_

private

Definition at line 95 of file GEMRecHitMatcher.h.

◆ gemRecHitH_

edm::Handle<GEMRecHitCollection> GEMRecHitMatcher::gemRecHitH_

private

Definition at line 92 of file GEMRecHitMatcher.h.

◆ gemRecHitToken_

edm::EDGetTokenT<GEMRecHitCollection> GEMRecHitMatcher::gemRecHitToken_

private

Definition at line 91 of file GEMRecHitMatcher.h.

◆ maxBX_

int GEMRecHitMatcher::maxBX_

private

Definition at line 97 of file GEMRecHitMatcher.h.

◆ minBX_

int GEMRecHitMatcher::minBX_

private

Definition at line 97 of file GEMRecHitMatcher.h.

◆ no_recHits_

const GEMRecHitContainer GEMRecHitMatcher::no_recHits_

private

Definition at line 104 of file GEMRecHitMatcher.h.

◆ recHits_

GEMRecHitContainer GEMRecHitMatcher::recHits_

private

Definition at line 105 of file GEMRecHitMatcher.h.

Referenced by recHits().

◆ superchamber_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::superchamber_to_recHits_

private

Definition at line 102 of file GEMRecHitMatcher.h.

◆ verbose_

bool GEMRecHitMatcher::verbose_

private

Definition at line 98 of file GEMRecHitMatcher.h.

Referenced by verbose().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GEMRecHitMatcher()

◆ ~GEMRecHitMatcher()

Member Function Documentation

◆ areGEMRecHitSame()

◆ chamberIds()

◆ detIds()

◆ gemDigiMatcher()

◆ init()

◆ isGEMRecHitMatched()

◆ match()

◆ matchRecHitsToSimTrack()

◆ nGEMRecHits()

◆ nLayersWithRecHitsInSuperChamber()

◆ partitionNumbers()

◆ recHitInContainer()

◆ recHitMeanPosition()

◆ recHitPosition()

◆ recHits()

◆ recHitsInChamber()

◆ recHitsInDetId()

◆ recHitsInSuperChamber()

◆ stripNumbersInDetId()

◆ superChamberIds()

◆ verbose()

Member Data Documentation

◆ chamber_to_recHits_

◆ detid_to_recHits_

◆ gem_geom_

◆ gemDigiMatcher_

◆ gemGeometry_

◆ gemRecHitH_

◆ gemRecHitToken_

◆ maxBX_

◆ minBX_

◆ no_recHits_

◆ recHits_

◆ superchamber_to_recHits_

◆ verbose_