#include <GEMDigiMatcher.h>

Public Member Functions
const std::map< unsigned int, GEMPadDigiContainer >	allPads () const

std::set< unsigned int >	chamberIdsCluster (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	chamberIdsDigi (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	chamberIdsPad (int gem_type=MuonHitHelper::GEM_ALL) const

const GEMPadDigiClusterContainer &	clustersInChamber (unsigned int) const

const GEMPadDigiClusterContainer &	clustersInDetId (unsigned int) const

const GEMPadDigiClusterContainer &	clustersInSuperChamber (unsigned int) const

const GEMCoPadDigiContainer &	coPadsInSuperChamber (unsigned int) const

std::set< unsigned int >	detIdsCluster (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	detIdsDigi (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	detIdsPad (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	detIdsSimLink (int gem_type=MuonHitHelper::GEM_ALL) const

const GEMDigiContainer &	digisInChamber (unsigned int) const

const GEMDigiContainer &	digisInDetId (unsigned int) const

const GEMDigiContainer &	digisInSuperChamber (unsigned int) const

int	extrapolateHsfromGEMPad (unsigned int, int) const

int	extrapolateHsfromGEMStrip (unsigned int, int) const

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

GlobalPoint	getGlobalPointDigi (unsigned int rawId, const GEMDigi &d) const

GlobalPoint	getGlobalPointPad (unsigned int rawId, const GEMPadDigi &tp) const

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

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

std::shared_ptr< GEMSimHitMatcher >	muonSimHitMatcher ()

int	nCoPads () const
	How many coincidence pads in GEM did this simtrack get in total? More...

int	nLayersWithClustersInSuperChamber (unsigned int) const

int	nLayersWithDigisInSuperChamber (unsigned int) const

int	nLayersWithPadsInSuperChamber (unsigned int) const

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

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

const GEMPadDigiContainer &	padsInChamber (unsigned int) const

const GEMPadDigiContainer &	padsInDetId (unsigned int) const

const GEMPadDigiContainer &	padsInSuperChamber (unsigned int) const

std::set< int >	partitionNumbers () const

std::set< int >	partitionNumbersWithCoPads () const

const GEMDigiSimLinkContainer &	simLinksInDetId (unsigned int) const

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

std::set< unsigned int >	superChamberIdsCluster (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	superChamberIdsCoPad (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	superChamberIdsDigi (int gem_type=MuonHitHelper::GEM_ALL) const

std::set< unsigned int >	superChamberIdsPad (int gem_type=MuonHitHelper::GEM_ALL) const

	~GEMDigiMatcher ()

Private Member Functions
void	clear ()

void	matchClustersToSimTrack (const GEMPadDigiClusterCollection &)

void	matchCoPadsToSimTrack (const GEMCoPadDigiCollection &)

void	matchDigisSLToSimTrack (const edm::DetSetVector< GEMDigiSimLink > &)

void	matchDigisToSimTrack (const GEMDigiCollection &)

void	matchPadsToSimTrack (const GEMPadDigiCollection &)

template<class T >
std::set< unsigned int >	selectDetIds (const T &, int) const

Private Attributes
std::map< unsigned int, GEMPadDigiClusterContainer >	chamber_to_clusters_

std::map< unsigned int, GEMDigiContainer >	chamber_to_digis_

std::map< unsigned int, GEMPadDigiContainer >	chamber_to_pads_

std::map< unsigned int, GEMPadDigiClusterContainer >	detid_to_clusters_

std::map< unsigned int, GEMDigiContainer >	detid_to_digis_

std::map< unsigned int, GEMPadDigiContainer >	detid_to_pads_

std::map< unsigned int, GEMDigiSimLinkContainer >	detid_to_simLinks_

bool	discardEleHits_

edm::Handle< GEMPadDigiClusterCollection >	gemClustersH_

edm::EDGetTokenT< GEMPadDigiClusterCollection >	gemClusterToken_

edm::Handle< GEMCoPadDigiCollection >	gemCoPadsH_

edm::EDGetTokenT< GEMCoPadDigiCollection >	gemCoPadToken_

edm::Handle< GEMDigiCollection >	gemDigisH_

edm::Handle< edm::DetSetVector< GEMDigiSimLink > >	gemDigisSLH_

edm::EDGetTokenT< GEMDigiCollection >	gemDigiToken_

const GEMGeometry *	gemGeometry_

edm::Handle< GEMPadDigiCollection >	gemPadsH_

edm::EDGetTokenT< GEMPadDigiCollection >	gemPadToken_

edm::EDGetTokenT< edm::DetSetVector< GEMDigiSimLink > >	gemSimLinkToken_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_

int	matchDeltaStrip_

bool	matchToSimLink_

int	maxBXCluster_

int	maxBXCoPad_

int	maxBXDigi_

int	maxBXPad_

int	minBXCluster_

int	minBXCoPad_

int	minBXDigi_

int	minBXPad_

std::shared_ptr< GEMSimHitMatcher >	muonSimHitMatcher_

GEMPadDigiClusterContainer	no_gem_clusters_

GEMCoPadDigiContainer	no_gem_copads_

GEMDigiContainer	no_gem_digis_

GEMPadDigiContainer	no_gem_pads_

GEMDigiSimLinkContainer	no_gem_simLinks_

bool	simMuOnly_

std::map< unsigned int, GEMPadDigiClusterContainer >	superchamber_to_clusters_

std::map< unsigned int, GEMCoPadDigiContainer >	superchamber_to_copads_

std::map< unsigned int, GEMDigiContainer >	superchamber_to_digis_

std::map< unsigned int, GEMPadDigiContainer >	superchamber_to_pads_

bool	verboseCluster_

bool	verboseCoPad_

bool	verboseDigi_

bool	verbosePad_

bool	verboseSimLink_

Detailed Description

Definition at line 31 of file GEMDigiMatcher.h.

Constructor & Destructor Documentation

◆ GEMDigiMatcher()

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

Definition at line 5 of file GEMDigiMatcher.cc.

                                                                                      {
   const auto& gemSimLink = pset.getParameterSet("gemSimLink");
   simMuOnly_ = gemSimLink.getParameter<bool>("simMuOnly");
   discardEleHits_ = gemSimLink.getParameter<bool>("discardEleHits");
   verboseSimLink_ = gemSimLink.getParameter<int>("verbose");
  
   const auto& gemDigi = pset.getParameterSet("gemStripDigi");
   minBXDigi_ = gemDigi.getParameter<int>("minBX");
   maxBXDigi_ = gemDigi.getParameter<int>("maxBX");
   matchDeltaStrip_ = gemDigi.getParameter<int>("matchDeltaStrip");
   verboseDigi_ = gemDigi.getParameter<int>("verbose");
   matchToSimLink_ = gemDigi.getParameter<bool>("matchToSimLink");
  
   const auto& gemPad = pset.getParameterSet("gemPadDigi");
   minBXPad_ = gemPad.getParameter<int>("minBX");
   maxBXPad_ = gemPad.getParameter<int>("maxBX");
   verbosePad_ = gemPad.getParameter<int>("verbose");
  
   const auto& gemCluster = pset.getParameterSet("gemPadCluster");
   minBXCluster_ = gemCluster.getParameter<int>("minBX");
   maxBXCluster_ = gemCluster.getParameter<int>("maxBX");
   verboseCluster_ = gemCluster.getParameter<int>("verbose");
  
   const auto& gemCoPad = pset.getParameterSet("gemCoPadDigi");
   minBXCoPad_ = gemCoPad.getParameter<int>("minBX");
   maxBXCoPad_ = gemCoPad.getParameter<int>("maxBX");
   verboseCoPad_ = gemCoPad.getParameter<int>("verbose");
  
   // make a new simhits matcher
   muonSimHitMatcher_.reset(new GEMSimHitMatcher(pset, std::move(iC)));
  
   if (matchToSimLink_)
     gemSimLinkToken_ =
         iC.consumes<edm::DetSetVector<GEMDigiSimLink>>(gemSimLink.getParameter<edm::InputTag>("inputTag"));
   gemDigiToken_ = iC.consumes<GEMDigiCollection>(gemDigi.getParameter<edm::InputTag>("inputTag"));
   gemPadToken_ = iC.consumes<GEMPadDigiCollection>(gemPad.getParameter<edm::InputTag>("inputTag"));
   gemClusterToken_ = iC.consumes<GEMPadDigiClusterCollection>(gemCluster.getParameter<edm::InputTag>("inputTag"));
   gemCoPadToken_ = iC.consumes<GEMCoPadDigiCollection>(gemCoPad.getParameter<edm::InputTag>("inputTag"));
  
   geomToken_ = iC.esConsumes<GEMGeometry, MuonGeometryRecord>();
 }

References muonGEMDigiPSet::gemSimLink, eostools::move(), and muonDTDigis_cfi::pset.

◆ ~GEMDigiMatcher()

GEMDigiMatcher::~GEMDigiMatcher ( )

inline

Definition at line 37 of file GEMDigiMatcher.h.

37 {}

Member Function Documentation

◆ allPads()

const std::map<unsigned int, GEMPadDigiContainer> GEMDigiMatcher::allPads ( ) const

inline

Definition at line 83 of file GEMDigiMatcher.h.

83 { return detid_to_pads_; }

References detid_to_pads_.

◆ chamberIdsCluster()

std::set< unsigned int > GEMDigiMatcher::chamberIdsCluster ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 318 of file GEMDigiMatcher.cc.

                                                                          {
   return selectDetIds(chamber_to_clusters_, gem_type);
 }

◆ chamberIdsDigi()

std::set< unsigned int > GEMDigiMatcher::chamberIdsDigi ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 310 of file GEMDigiMatcher.cc.

                                                                       {
   return selectDetIds(chamber_to_digis_, gem_type);
 }

◆ chamberIdsPad()

std::set< unsigned int > GEMDigiMatcher::chamberIdsPad ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 314 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(chamber_to_pads_, gem_type);
 }

◆ clear()

void GEMDigiMatcher::clear ( void )

private

Definition at line 515 of file GEMDigiMatcher.cc.

                            {
   detid_to_simLinks_.clear();
  
   detid_to_digis_.clear();
   chamber_to_digis_.clear();
   superchamber_to_digis_.clear();
  
   detid_to_pads_.clear();
   chamber_to_pads_.clear();
   superchamber_to_pads_.clear();
  
   detid_to_clusters_.clear();
   chamber_to_clusters_.clear();
   superchamber_to_clusters_.clear();
  
   superchamber_to_copads_.clear();
 }

◆ clustersInChamber()

const GEMPadDigiClusterContainer & GEMDigiMatcher::clustersInChamber ( unsigned int detid ) const

Definition at line 380 of file GEMDigiMatcher.cc.

                                                                                             {
   if (chamber_to_clusters_.find(detid) == chamber_to_clusters_.end())
     return no_gem_clusters_;
   return chamber_to_clusters_.at(detid);
 }

◆ clustersInDetId()

const GEMPadDigiClusterContainer & GEMDigiMatcher::clustersInDetId ( unsigned int detid ) const

Definition at line 374 of file GEMDigiMatcher.cc.

                                                                                           {
   if (detid_to_clusters_.find(detid) == detid_to_clusters_.end())
     return no_gem_clusters_;
   return detid_to_clusters_.at(detid);
 }

◆ clustersInSuperChamber()

const GEMPadDigiClusterContainer & GEMDigiMatcher::clustersInSuperChamber ( unsigned int detid ) const

Definition at line 386 of file GEMDigiMatcher.cc.

                                                                                                  {
   if (superchamber_to_clusters_.find(detid) == superchamber_to_clusters_.end())
     return no_gem_clusters_;
   return superchamber_to_clusters_.at(detid);
 }

◆ coPadsInSuperChamber()

const GEMCoPadDigiContainer & GEMDigiMatcher::coPadsInSuperChamber ( unsigned int detid ) const

Definition at line 392 of file GEMDigiMatcher.cc.

                                                                                           {
   if (superchamber_to_copads_.find(detid) == superchamber_to_copads_.end())
     return no_gem_copads_;
   return superchamber_to_copads_.at(detid);
 }

◆ detIdsCluster()

std::set< unsigned int > GEMDigiMatcher::detIdsCluster ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 306 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(detid_to_clusters_, gem_type);
 }

◆ detIdsDigi()

std::set< unsigned int > GEMDigiMatcher::detIdsDigi ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 300 of file GEMDigiMatcher.cc.

                                                                   {
   return selectDetIds(detid_to_digis_, gem_type);
 }

◆ detIdsPad()

std::set< unsigned int > GEMDigiMatcher::detIdsPad ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 304 of file GEMDigiMatcher.cc.

304 { return selectDetIds(detid_to_pads_, gem_type); }

◆ detIdsSimLink()

std::set< unsigned int > GEMDigiMatcher::detIdsSimLink ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 296 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(detid_to_simLinks_, gem_type);
 }

◆ digisInChamber()

const GEMDigiContainer & GEMDigiMatcher::digisInChamber ( unsigned int detid ) const

Definition at line 344 of file GEMDigiMatcher.cc.

                                                                                {
   if (chamber_to_digis_.find(detid) == chamber_to_digis_.end())
     return no_gem_digis_;
   return chamber_to_digis_.at(detid);
 }

◆ digisInDetId()

const GEMDigiContainer & GEMDigiMatcher::digisInDetId ( unsigned int detid ) const

Definition at line 338 of file GEMDigiMatcher.cc.

                                                                              {
   if (detid_to_digis_.find(detid) == detid_to_digis_.end())
     return no_gem_digis_;
   return detid_to_digis_.at(detid);
 }

◆ digisInSuperChamber()

const GEMDigiContainer & GEMDigiMatcher::digisInSuperChamber ( unsigned int detid ) const

Definition at line 350 of file GEMDigiMatcher.cc.

                                                                                     {
   if (superchamber_to_digis_.find(detid) == superchamber_to_digis_.end())
     return no_gem_digis_;
   return superchamber_to_digis_.at(detid);
 }

◆ extrapolateHsfromGEMPad()

int GEMDigiMatcher::extrapolateHsfromGEMPad	(	unsigned int	,
		int
	)		const

◆ extrapolateHsfromGEMStrip()

int GEMDigiMatcher::extrapolateHsfromGEMStrip	(	unsigned int	,
		int
	)		const

◆ getGlobalPointDigi()

GlobalPoint GEMDigiMatcher::getGlobalPointDigi	(	unsigned int	rawId,
		const GEMDigi &	d
	)		const

Definition at line 501 of file GEMDigiMatcher.cc.

                                                                                          {
   GEMDetId gem_id(rawId);
   const LocalPoint& gem_lp = gemGeometry_->etaPartition(gem_id)->centreOfStrip(d.strip());
   const GlobalPoint& gem_gp = gemGeometry_->idToDet(gem_id)->surface().toGlobal(gem_lp);
   return gem_gp;
 }

References ztail::d.

◆ getGlobalPointPad()

GlobalPoint GEMDigiMatcher::getGlobalPointPad	(	unsigned int	rawId,
		const GEMPadDigi &	tp
	)		const

Definition at line 508 of file GEMDigiMatcher.cc.

                                                                                             {
   GEMDetId gem_id(rawId);
   const LocalPoint& gem_lp = gemGeometry_->etaPartition(gem_id)->centreOfPad(tp.pad());
   const GlobalPoint& gem_gp = gemGeometry_->idToDet(gem_id)->surface().toGlobal(gem_lp);
   return gem_gp;
 }

References cmsswSequenceInfo::tp.

◆ init()

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

Definition at line 47 of file GEMDigiMatcher.cc.

                                                                              {
   muonSimHitMatcher_->init(iEvent, iSetup);
  
   if (matchToSimLink_)
     iEvent.getByToken(gemSimLinkToken_, gemDigisSLH_);
   iEvent.getByToken(gemDigiToken_, gemDigisH_);
   iEvent.getByToken(gemPadToken_, gemPadsH_);
   iEvent.getByToken(gemClusterToken_, gemClustersH_);
   iEvent.getByToken(gemCoPadToken_, gemCoPadsH_);
  
   gemGeometry_ = &iSetup.getData(geomToken_);
 }

References edm::EventSetup::getData(), and iEvent.

◆ match()

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

do the matching

Definition at line 61 of file GEMDigiMatcher.cc.

                                                                 {
   // match simhits first
   muonSimHitMatcher_->match(t, v);
  
   // get the digi collections
   const GEMDigiCollection& gemDigis = *gemDigisH_.product();
   const GEMPadDigiCollection& gemPads = *gemPadsH_.product();
   const GEMPadDigiClusterCollection& gemClusters = *gemClustersH_.product();
   const GEMCoPadDigiCollection& gemCoPads = *gemCoPadsH_.product();
  
   clear();
  
   // hard cut on non-GEM muons
   if (std::abs(t.momentum().eta()) < 1.55)
     return;
  
   // now match the digis
   if (matchToSimLink_) {
     const edm::DetSetVector<GEMDigiSimLink>& gemDigisSL = *gemDigisSLH_.product();
     matchDigisSLToSimTrack(gemDigisSL);
   }
   matchDigisToSimTrack(gemDigis);
   matchPadsToSimTrack(gemPads);
   matchClustersToSimTrack(gemClusters);
   matchCoPadsToSimTrack(gemCoPads);
 }

References funct::abs(), clear(), submitPVValidationJobs::t, and findQualityFiles::v.

◆ matchClustersToSimTrack()

void GEMDigiMatcher::matchClustersToSimTrack ( const GEMPadDigiClusterCollection & clusters )

private

Definition at line 215 of file GEMDigiMatcher.cc.

                                                                                         {
   const auto& det_ids = muonSimHitMatcher_->detIds();
   for (auto id : det_ids) {
     GEMDetId p_id(id);
  
     auto clusters_in_det = clusters.get(p_id);
  
     for (auto cluster = clusters_in_det.first; cluster != clusters_in_det.second; ++cluster) {
       bool isMatched;
  
       // ignore 16-partition GE2/1 pads
       if (p_id.isGE21() and cluster->nPartitions() == GEMPadDigiCluster::GE21SplitStrip)
         continue;
  
       // check that the cluster BX is within the range
       if (cluster->bx() < minBXCluster_ || cluster->bx() > maxBXCluster_)
         continue;
  
       if (verboseCluster_)
         edm::LogInfo("GEMDigiMatcher") << "GEMCluster " << p_id << " " << *cluster << endl;
  
       // check that at least one pad was hit by the track
       for (const auto& p : cluster->pads()) {
         for (auto pad : detid_to_pads_[id]) {
           if (pad.pad() == p) {
             isMatched = true;
           }
         }
       }
       if (isMatched) {
         detid_to_clusters_[id].push_back(*cluster);
         chamber_to_clusters_[p_id.chamberId().rawId()].push_back(*cluster);
         superchamber_to_clusters_[p_id.superChamberId().rawId()].push_back(*cluster);
         if (verboseCluster_)
           edm::LogInfo("GEMDigiMatcher") << "...was matched!" << endl;
       }
     }
   }
 }

References GEMDetId::chamberId(), bsc_activity_cfg::clusters, GEMPadDigiCluster::GE21SplitStrip, triggerObjects_cff::id, GEMDetId::isGE21(), trackerHitRTTI::isMatched(), AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), and GEMDetId::superChamberId().

◆ matchCoPadsToSimTrack()

void GEMDigiMatcher::matchCoPadsToSimTrack ( const GEMCoPadDigiCollection & co_pads )

private

Definition at line 255 of file GEMDigiMatcher.cc.

                                                                                 {
   // loop on the GEM detids
   for (auto d : superChamberIdsPad()) {
     GEMDetId id(d);
  
     const auto& co_pads_in_det = co_pads.get(id);
     for (auto copad = co_pads_in_det.first; copad != co_pads_in_det.second; ++copad) {
       // check that the cluster BX is within the range
       if (copad->bx(1) < minBXCoPad_ || copad->bx(1) > maxBXCoPad_)
         continue;
  
       if (verboseCoPad_)
         edm::LogInfo("GEMDigiMatcher") << "GEMCoPadDigi: " << id << " " << *copad << endl;
  
       bool isMatchedL1 = false;
       bool isMatchedL2 = false;
       GEMDetId gemL1_id(id.region(), 1, id.station(), 1, id.chamber(), copad->roll());
       GEMDetId gemL2_id(id.region(), 1, id.station(), 2, id.chamber(), 0);
  
       // first pad is tightly matched
       for (const auto& p : padsInDetId(gemL1_id.rawId())) {
         if (p == copad->first()) {
           isMatchedL1 = true;
         }
       }
  
       // second pad can only be loosely matched
       for (const auto& p : padsInChamber(gemL2_id.rawId())) {
         if (p == copad->second()) {
           isMatchedL2 = true;
         }
       }
       if (isMatchedL1 and isMatchedL2) {
         superchamber_to_copads_[id.rawId()].push_back(*copad);
         if (verboseCoPad_)
           edm::LogInfo("GEMDigiMatcher") << "...was matched! " << endl;
       }
     }
   }
 }

References relativeConstraints::chamber, ztail::d, triggerObjects_cff::id, AlCaHLTBitMon_ParallelJobs::p, HLT_FULL_cff::region, and relativeConstraints::station.

◆ matchDigisSLToSimTrack()

void GEMDigiMatcher::matchDigisSLToSimTrack ( const edm::DetSetVector< GEMDigiSimLink > & digisSL )

private

Definition at line 88 of file GEMDigiMatcher.cc.

                                                                                           {
   if (verboseSimLink_)
     edm::LogInfo("GEMDigiMatcher") << "Matching simtrack to GEM simlinks" << endl;
  
   // loop on the simlinks
   for (auto itsimlink = digisSL.begin(); itsimlink != digisSL.end(); itsimlink++) {
     for (auto sl = itsimlink->data.begin(); sl != itsimlink->data.end(); ++sl) {
       GEMDetId p_id(sl->getDetUnitId());
  
       // ignore simlinks in non-matched chambers
       const auto& detids(muonSimHitMatcher_->detIds());
       if (detids.find(p_id.rawId()) == detids.end())
         continue;
  
       // no simhits in this chamber!
       if (muonSimHitMatcher_->hitsInDetId(p_id.rawId()).empty())
         continue;
  
       if (verboseSimLink_)
         edm::LogInfo("GEMDigiMatcher") << "GEMDigiSimLink " << p_id << " " << sl->getStrip() << " " << sl->getBx()
                                        << " " << sl->getEnergyLoss() << " " << sl->getTimeOfFlight() << " "
                                        << sl->getParticleType() << std::endl;
  
       // consider only the muon hits
       if (simMuOnly_ && std::abs(sl->getParticleType()) != 13)
         continue;
  
       // discard electron hits in the GEM chambers
       if (discardEleHits_ && std::abs(sl->getParticleType()) == 11)
         continue;
  
       // loop on the matched simhits
       for (const auto& simhit : muonSimHitMatcher_->hitsInDetId(p_id.rawId())) {
         // check if the simhit properties agree
         if (simhit.particleType() == sl->getParticleType() and simhit.trackId() == sl->getTrackId() and
             std::abs(simhit.energyLoss() - sl->getEnergyLoss()) < 0.001 and
             std::abs(simhit.timeOfFlight() - sl->getTimeOfFlight()) < 0.001 and
             simhit.entryPoint() == sl->getEntryPoint() and simhit.momentumAtEntry() == sl->getMomentumAtEntry()) {
           detid_to_simLinks_[p_id.rawId()].push_back(*sl);
           if (verboseSimLink_)
             edm::LogInfo("GEMDigiMatcher") << "...was matched!" << endl;
           break;
         }
       }
     }
   }
 }

References funct::abs(), edm::DetSetVector< T >::begin(), relativeConstraints::empty, and edm::DetSetVector< T >::end().

◆ matchDigisToSimTrack()

void GEMDigiMatcher::matchDigisToSimTrack ( const GEMDigiCollection & digis )

private

Definition at line 136 of file GEMDigiMatcher.cc.

                                                                         {
   if (verboseDigi_)
     edm::LogInfo("GEMDigiMatcher") << "Matching simtrack to GEM digis" << endl;
   for (auto id : muonSimHitMatcher_->detIds()) {
     GEMDetId p_id(id);
     const auto& hit_strips = muonSimHitMatcher_->hitStripsInDetId(id, matchDeltaStrip_);
     const auto& digis_in_det = digis.get(p_id);
  
     for (auto d = digis_in_det.first; d != digis_in_det.second; ++d) {
       bool isMatched = false;
  
       // check that the digi is within BX range
       if (d->bx() < minBXDigi_ || d->bx() > maxBXDigi_)
         continue;
  
       if (verboseDigi_)
         edm::LogInfo("GEMDigiMatcher") << "GEMDigi " << p_id << " " << *d << endl;
  
       // GEN-SIM-DIGI-RAW monte carlo
       if (matchToSimLink_) {
         // check that the digi matches to at least one GEMDigiSimLink
         for (const auto& sl : detid_to_simLinks_[p_id.rawId()]) {
           if (sl.getStrip() == d->strip() and sl.getBx() == d->bx()) {
             isMatched = true;
             break;
           }
         }
       }
       // GEN-SIM-RAW monte carlo
       else {
         // check that it matches a strip that was hit by SimHits from our track
         if (hit_strips.find(d->strip()) != hit_strips.end()) {
           isMatched = true;
         }
       }
       if (isMatched) {
         detid_to_digis_[p_id.rawId()].push_back(*d);
         chamber_to_digis_[p_id.chamberId().rawId()].push_back(*d);
         superchamber_to_digis_[p_id.superChamberId().rawId()].push_back(*d);
         if (verboseDigi_)
           edm::LogInfo("GEMDigiMatcher") << "...was matched!" << endl;
       }
     }
   }
 }

References GEMDetId::chamberId(), ztail::d, trackerHitRTTI::isMatched(), DetId::rawId(), and GEMDetId::superChamberId().

◆ matchPadsToSimTrack()

void GEMDigiMatcher::matchPadsToSimTrack ( const GEMPadDigiCollection & pads )

private

Definition at line 181 of file GEMDigiMatcher.cc.

                                                                          {
   const auto& det_ids = muonSimHitMatcher_->detIds();
   for (const auto& id : det_ids) {
     GEMDetId p_id(id);
  
     const auto& pads_in_det = pads.get(p_id);
  
     for (auto pad = pads_in_det.first; pad != pads_in_det.second; ++pad) {
       // ignore 16-partition GE2/1 pads
       if (p_id.isGE21() and pad->nPartitions() == GEMPadDigi::GE21SplitStrip)
         continue;
  
       // check that the pad BX is within the range
       if (pad->bx() < minBXPad_ || pad->bx() > maxBXPad_)
         continue;
  
       if (verbosePad_)
         edm::LogInfo("GEMDigiMatcher") << "GEMPad " << p_id << " " << *pad << endl;
  
       // check that it matches a pad that was hit by SimHits from our track
       for (auto digi : detid_to_digis_[p_id.rawId()]) {
         if (digi.strip() / 2 == pad->pad()) {
           detid_to_pads_[p_id.rawId()].push_back(*pad);
           chamber_to_pads_[p_id.chamberId().rawId()].push_back(*pad);
           superchamber_to_pads_[p_id.superChamberId().rawId()].push_back(*pad);
           if (verbosePad_)
             edm::LogInfo("GEMDigiMatcher") << "...was matched!" << endl;
           break;
         }
       }
     }
   }
 }

References GEMDetId::chamberId(), GEMPadDigi::GE21SplitStrip, GEMDetId::isGE21(), DetId::rawId(), and GEMDetId::superChamberId().

◆ muonSimHitMatcher()

std::shared_ptr<GEMSimHitMatcher> GEMDigiMatcher::muonSimHitMatcher ( )

inline

Definition at line 110 of file GEMDigiMatcher.h.

110 { return muonSimHitMatcher_; }

References muonSimHitMatcher_.

◆ nCoPads()

int GEMDigiMatcher::nCoPads ( ) const

How many coincidence pads in GEM did this simtrack get in total?

Definition at line 452 of file GEMDigiMatcher.cc.

                                   {
   int n = 0;
   const auto& ids = superChamberIdsCoPad();
   for (const auto& id : ids) {
     n += coPadsInSuperChamber(id).size();
   }
   return n;
 }

References dqmiodumpmetadata::n.

◆ nLayersWithClustersInSuperChamber()

int GEMDigiMatcher::nLayersWithClustersInSuperChamber ( unsigned int detid ) const

Definition at line 428 of file GEMDigiMatcher.cc.

                                                                               {
   set<int> layers;
   GEMDetId sch_id(detid);
   for (int iLayer = 1; iLayer <= 2; iLayer++) {
     GEMDetId ch_id(sch_id.region(), sch_id.ring(), sch_id.station(), iLayer, sch_id.chamber(), 0);
     // get the pads in this chamber
     const auto& clusters = clustersInChamber(ch_id.rawId());
     // at least one digi in this layer!
     if (!clusters.empty()) {
       layers.insert(iLayer);
     }
   }
   return layers.size();
 }

References GEMDetId::chamber(), bsc_activity_cfg::clusters, hgcalTopologyTester_cfi::layers, GEMDetId::region(), GEMDetId::ring(), and GEMDetId::station().

◆ nLayersWithDigisInSuperChamber()

int GEMDigiMatcher::nLayersWithDigisInSuperChamber ( unsigned int detid ) const

Definition at line 398 of file GEMDigiMatcher.cc.

                                                                            {
   set<int> layers;
   GEMDetId sch_id(detid);
   for (int iLayer = 1; iLayer <= 2; iLayer++) {
     GEMDetId ch_id(sch_id.region(), sch_id.ring(), sch_id.station(), iLayer, sch_id.chamber(), 0);
     // get the digis in this chamber
     const auto& digis = digisInChamber(ch_id.rawId());
     // at least one digi in this layer!
     if (!digis.empty()) {
       layers.insert(iLayer);
     }
   }
   return layers.size();
 }

References GEMDetId::chamber(), hgcalTopologyTester_cfi::layers, GEMDetId::region(), GEMDetId::ring(), and GEMDetId::station().

◆ nLayersWithPadsInSuperChamber()

int GEMDigiMatcher::nLayersWithPadsInSuperChamber ( unsigned int detid ) const

Definition at line 413 of file GEMDigiMatcher.cc.

                                                                           {
   set<int> layers;
   GEMDetId sch_id(detid);
   for (int iLayer = 1; iLayer <= 2; iLayer++) {
     GEMDetId ch_id(sch_id.region(), sch_id.ring(), sch_id.station(), iLayer, sch_id.chamber(), 0);
     // get the pads in this chamber
     const auto& pads = padsInChamber(ch_id.rawId());
     // at least one digi in this layer!
     if (!pads.empty()) {
       layers.insert(iLayer);
     }
   }
   return layers.size();
 }

References GEMDetId::chamber(), hgcalTopologyTester_cfi::layers, GEMDetId::region(), GEMDetId::ring(), and GEMDetId::station().

◆ nPads()

int GEMDigiMatcher::nPads ( ) const

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

Definition at line 443 of file GEMDigiMatcher.cc.

                                 {
   int n = 0;
   const auto& ids = superChamberIdsPad();
   for (const auto& id : ids) {
     n += padsInSuperChamber(id).size();
   }
   return n;
 }

References dqmiodumpmetadata::n.

◆ padNumbersInDetId()

std::set< int > GEMDigiMatcher::padNumbersInDetId ( unsigned int detid ) const

Definition at line 470 of file GEMDigiMatcher.cc.

                                                                       {
   set<int> result;
   const auto& digis = padsInDetId(detid);
   for (const auto& d : digis) {
     result.insert(d.pad());
   }
   return result;
 }

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

◆ padsInChamber()

const GEMPadDigiContainer & GEMDigiMatcher::padsInChamber ( unsigned int detid ) const

Definition at line 362 of file GEMDigiMatcher.cc.

                                                                                  {
   if (chamber_to_pads_.find(detid) == chamber_to_pads_.end())
     return no_gem_pads_;
   return chamber_to_pads_.at(detid);
 }

◆ padsInDetId()

const GEMPadDigiContainer & GEMDigiMatcher::padsInDetId ( unsigned int detid ) const

Definition at line 356 of file GEMDigiMatcher.cc.

                                                                                {
   if (detid_to_pads_.find(detid) == detid_to_pads_.end())
     return no_gem_pads_;
   return detid_to_pads_.at(detid);
 }

◆ padsInSuperChamber()

const GEMPadDigiContainer & GEMDigiMatcher::padsInSuperChamber ( unsigned int detid ) const

Definition at line 368 of file GEMDigiMatcher.cc.

                                                                                       {
   if (superchamber_to_pads_.find(detid) == superchamber_to_pads_.end())
     return no_gem_pads_;
   return superchamber_to_pads_.at(detid);
 }

◆ partitionNumbers()

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

Definition at line 479 of file GEMDigiMatcher.cc.

                                                    {
   std::set<int> result;
  
   const auto& detids = detIdsDigi();
   for (const auto& id : detids) {
     const GEMDetId& idd(id);
     result.insert(idd.roll());
   }
   return result;
 }

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

◆ partitionNumbersWithCoPads()

std::set< int > GEMDigiMatcher::partitionNumbersWithCoPads ( ) const

Definition at line 490 of file GEMDigiMatcher.cc.

                                                              {
   std::set<int> result;
  
   const auto& detids = superChamberIdsCoPad();
   for (const auto& id : detids) {
     const GEMDetId& idd(id);
     result.insert(idd.roll());
   }
   return result;
 }

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

◆ selectDetIds()

template<class T >

std::set< unsigned int > GEMDigiMatcher::selectDetIds	(	const T &	digis,
		int	gem_type
	)		const

private

Definition at line 183 of file GEMDigiMatcher.h.

                                                                                     {
   std::set<unsigned int> result;
   for (const auto& p : digis) {
     const auto& id = p.first;
     if (gem_type > 0) {
       GEMDetId detId(id);
       if (MuonHitHelper::toGEMType(detId.station(), detId.ring()) != gem_type)
         continue;
     }
     result.insert(p.first);
   }
   return result;
 }

References AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, GEMDetId::ring(), GEMDetId::station(), and MuonHitHelper::toGEMType().

◆ simLinksInDetId()

const GEMDigiSimLinkContainer& GEMDigiMatcher::simLinksInDetId ( unsigned int ) const

◆ stripNumbersInDetId()

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

Definition at line 461 of file GEMDigiMatcher.cc.

                                                                         {
   set<int> result;
   const auto& digis = digisInDetId(detid);
   for (const auto& d : digis) {
     result.insert(d.strip());
   }
   return result;
 }

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

◆ superChamberIdsCluster()

std::set< unsigned int > GEMDigiMatcher::superChamberIdsCluster ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 330 of file GEMDigiMatcher.cc.

                                                                               {
   return selectDetIds(superchamber_to_clusters_, gem_type);
 }

◆ superChamberIdsCoPad()

std::set< unsigned int > GEMDigiMatcher::superChamberIdsCoPad ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 334 of file GEMDigiMatcher.cc.

                                                                             {
   return selectDetIds(superchamber_to_copads_, gem_type);
 }

◆ superChamberIdsDigi()

std::set< unsigned int > GEMDigiMatcher::superChamberIdsDigi ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 322 of file GEMDigiMatcher.cc.

                                                                            {
   return selectDetIds(superchamber_to_digis_, gem_type);
 }

◆ superChamberIdsPad()

std::set< unsigned int > GEMDigiMatcher::superChamberIdsPad ( int gem_type = MuonHitHelper::GEM_ALL ) const

Definition at line 326 of file GEMDigiMatcher.cc.

                                                                           {
   return selectDetIds(superchamber_to_pads_, gem_type);
 }

Member Data Documentation

◆ chamber_to_clusters_

std::map<unsigned int, GEMPadDigiClusterContainer> GEMDigiMatcher::chamber_to_clusters_

private

Definition at line 170 of file GEMDigiMatcher.h.

◆ chamber_to_digis_

std::map<unsigned int, GEMDigiContainer> GEMDigiMatcher::chamber_to_digis_

private

Definition at line 162 of file GEMDigiMatcher.h.

◆ chamber_to_pads_

std::map<unsigned int, GEMPadDigiContainer> GEMDigiMatcher::chamber_to_pads_

private

Definition at line 166 of file GEMDigiMatcher.h.

◆ detid_to_clusters_

std::map<unsigned int, GEMPadDigiClusterContainer> GEMDigiMatcher::detid_to_clusters_

private

Definition at line 169 of file GEMDigiMatcher.h.

◆ detid_to_digis_

std::map<unsigned int, GEMDigiContainer> GEMDigiMatcher::detid_to_digis_

private

Definition at line 161 of file GEMDigiMatcher.h.

◆ detid_to_pads_

std::map<unsigned int, GEMPadDigiContainer> GEMDigiMatcher::detid_to_pads_

private

Definition at line 165 of file GEMDigiMatcher.h.

Referenced by allPads().

◆ detid_to_simLinks_

std::map<unsigned int, GEMDigiSimLinkContainer> GEMDigiMatcher::detid_to_simLinks_

private

Definition at line 159 of file GEMDigiMatcher.h.

◆ discardEleHits_

bool GEMDigiMatcher::discardEleHits_

private

Definition at line 142 of file GEMDigiMatcher.h.

◆ gemClustersH_

edm::Handle<GEMPadDigiClusterCollection> GEMDigiMatcher::gemClustersH_

private

Definition at line 130 of file GEMDigiMatcher.h.

◆ gemClusterToken_

edm::EDGetTokenT<GEMPadDigiClusterCollection> GEMDigiMatcher::gemClusterToken_

private

Definition at line 124 of file GEMDigiMatcher.h.

◆ gemCoPadsH_

edm::Handle<GEMCoPadDigiCollection> GEMDigiMatcher::gemCoPadsH_

private

Definition at line 131 of file GEMDigiMatcher.h.

◆ gemCoPadToken_

edm::EDGetTokenT<GEMCoPadDigiCollection> GEMDigiMatcher::gemCoPadToken_

private

Definition at line 125 of file GEMDigiMatcher.h.

◆ gemDigisH_

edm::Handle<GEMDigiCollection> GEMDigiMatcher::gemDigisH_

private

Definition at line 128 of file GEMDigiMatcher.h.

◆ gemDigisSLH_

edm::Handle<edm::DetSetVector<GEMDigiSimLink> > GEMDigiMatcher::gemDigisSLH_

private

Definition at line 127 of file GEMDigiMatcher.h.

◆ gemDigiToken_

edm::EDGetTokenT<GEMDigiCollection> GEMDigiMatcher::gemDigiToken_

private

Definition at line 122 of file GEMDigiMatcher.h.

◆ gemGeometry_

const GEMGeometry* GEMDigiMatcher::gemGeometry_

private

Definition at line 136 of file GEMDigiMatcher.h.

◆ gemPadsH_

edm::Handle<GEMPadDigiCollection> GEMDigiMatcher::gemPadsH_

private

Definition at line 129 of file GEMDigiMatcher.h.

◆ gemPadToken_

edm::EDGetTokenT<GEMPadDigiCollection> GEMDigiMatcher::gemPadToken_

private

Definition at line 123 of file GEMDigiMatcher.h.

◆ gemSimLinkToken_

edm::EDGetTokenT<edm::DetSetVector<GEMDigiSimLink> > GEMDigiMatcher::gemSimLinkToken_

private

Definition at line 121 of file GEMDigiMatcher.h.

◆ geomToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMDigiMatcher::geomToken_

private

Definition at line 135 of file GEMDigiMatcher.h.

◆ matchDeltaStrip_

int GEMDigiMatcher::matchDeltaStrip_

private

Definition at line 149 of file GEMDigiMatcher.h.

◆ matchToSimLink_

bool GEMDigiMatcher::matchToSimLink_

private

Definition at line 157 of file GEMDigiMatcher.h.

◆ maxBXCluster_

int GEMDigiMatcher::maxBXCluster_

private

Definition at line 146 of file GEMDigiMatcher.h.

◆ maxBXCoPad_

int GEMDigiMatcher::maxBXCoPad_

private

Definition at line 147 of file GEMDigiMatcher.h.

◆ maxBXDigi_

int GEMDigiMatcher::maxBXDigi_

private

Definition at line 144 of file GEMDigiMatcher.h.

◆ maxBXPad_

int GEMDigiMatcher::maxBXPad_

private

Definition at line 145 of file GEMDigiMatcher.h.

◆ minBXCluster_

int GEMDigiMatcher::minBXCluster_

private

Definition at line 146 of file GEMDigiMatcher.h.

◆ minBXCoPad_

int GEMDigiMatcher::minBXCoPad_

private

Definition at line 147 of file GEMDigiMatcher.h.

◆ minBXDigi_

int GEMDigiMatcher::minBXDigi_

private

Definition at line 144 of file GEMDigiMatcher.h.

◆ minBXPad_

int GEMDigiMatcher::minBXPad_

private

Definition at line 145 of file GEMDigiMatcher.h.

◆ muonSimHitMatcher_

std::shared_ptr<GEMSimHitMatcher> GEMDigiMatcher::muonSimHitMatcher_

private

Definition at line 133 of file GEMDigiMatcher.h.

Referenced by muonSimHitMatcher().

◆ no_gem_clusters_

GEMPadDigiClusterContainer GEMDigiMatcher::no_gem_clusters_

private

Definition at line 178 of file GEMDigiMatcher.h.

◆ no_gem_copads_

GEMCoPadDigiContainer GEMDigiMatcher::no_gem_copads_

private

Definition at line 179 of file GEMDigiMatcher.h.

◆ no_gem_digis_

GEMDigiContainer GEMDigiMatcher::no_gem_digis_

private

Definition at line 176 of file GEMDigiMatcher.h.

◆ no_gem_pads_

GEMPadDigiContainer GEMDigiMatcher::no_gem_pads_

private

Definition at line 177 of file GEMDigiMatcher.h.

◆ no_gem_simLinks_

GEMDigiSimLinkContainer GEMDigiMatcher::no_gem_simLinks_

private

Definition at line 175 of file GEMDigiMatcher.h.

◆ simMuOnly_

bool GEMDigiMatcher::simMuOnly_

private

Definition at line 141 of file GEMDigiMatcher.h.

◆ superchamber_to_clusters_

std::map<unsigned int, GEMPadDigiClusterContainer> GEMDigiMatcher::superchamber_to_clusters_

private

Definition at line 171 of file GEMDigiMatcher.h.

◆ superchamber_to_copads_

std::map<unsigned int, GEMCoPadDigiContainer> GEMDigiMatcher::superchamber_to_copads_

private

Definition at line 173 of file GEMDigiMatcher.h.

◆ superchamber_to_digis_

std::map<unsigned int, GEMDigiContainer> GEMDigiMatcher::superchamber_to_digis_

private

Definition at line 163 of file GEMDigiMatcher.h.

◆ superchamber_to_pads_

std::map<unsigned int, GEMPadDigiContainer> GEMDigiMatcher::superchamber_to_pads_

private

Definition at line 167 of file GEMDigiMatcher.h.

◆ verboseCluster_

bool GEMDigiMatcher::verboseCluster_

private

Definition at line 154 of file GEMDigiMatcher.h.

◆ verboseCoPad_

bool GEMDigiMatcher::verboseCoPad_

private

Definition at line 155 of file GEMDigiMatcher.h.

◆ verboseDigi_

bool GEMDigiMatcher::verboseDigi_

private

Definition at line 152 of file GEMDigiMatcher.h.

◆ verbosePad_

bool GEMDigiMatcher::verbosePad_

private

Definition at line 153 of file GEMDigiMatcher.h.

◆ verboseSimLink_

bool GEMDigiMatcher::verboseSimLink_

private

Definition at line 151 of file GEMDigiMatcher.h.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GEMDigiMatcher()

◆ ~GEMDigiMatcher()

Member Function Documentation

◆ allPads()

◆ chamberIdsCluster()

◆ chamberIdsDigi()

◆ chamberIdsPad()

◆ clear()

◆ clustersInChamber()

◆ clustersInDetId()

◆ clustersInSuperChamber()

◆ coPadsInSuperChamber()

◆ detIdsCluster()

◆ detIdsDigi()

◆ detIdsPad()

◆ detIdsSimLink()

◆ digisInChamber()

◆ digisInDetId()

◆ digisInSuperChamber()

◆ extrapolateHsfromGEMPad()

◆ extrapolateHsfromGEMStrip()

◆ getGlobalPointDigi()

◆ getGlobalPointPad()

◆ init()

◆ match()

◆ matchClustersToSimTrack()

◆ matchCoPadsToSimTrack()

◆ matchDigisSLToSimTrack()

◆ matchDigisToSimTrack()

◆ matchPadsToSimTrack()

◆ muonSimHitMatcher()

◆ nCoPads()

◆ nLayersWithClustersInSuperChamber()

◆ nLayersWithDigisInSuperChamber()

◆ nLayersWithPadsInSuperChamber()

◆ nPads()

◆ padNumbersInDetId()

◆ padsInChamber()

◆ padsInDetId()

◆ padsInSuperChamber()

◆ partitionNumbers()

◆ partitionNumbersWithCoPads()

◆ selectDetIds()

◆ simLinksInDetId()

◆ stripNumbersInDetId()

◆ superChamberIdsCluster()

◆ superChamberIdsCoPad()

◆ superChamberIdsDigi()

◆ superChamberIdsPad()

Member Data Documentation

◆ chamber_to_clusters_

◆ chamber_to_digis_

◆ chamber_to_pads_

◆ detid_to_clusters_

◆ detid_to_digis_

◆ detid_to_pads_

◆ detid_to_simLinks_

◆ discardEleHits_

◆ gemClustersH_

◆ gemClusterToken_

◆ gemCoPadsH_

◆ gemCoPadToken_

◆ gemDigisH_

◆ gemDigisSLH_

◆ gemDigiToken_

◆ gemGeometry_

◆ gemPadsH_

◆ gemPadToken_

◆ gemSimLinkToken_

◆ geomToken_

◆ matchDeltaStrip_

◆ matchToSimLink_

◆ maxBXCluster_

◆ maxBXCoPad_

◆ maxBXDigi_