#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	matchPadsInDetId (const GEMDetId &simhitId, const GEMPadDigiCollection &)

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	(	edm::ParameterSet const &	iPS,
		edm::ConsumesCollector &&	iC
	)

Definition at line 5 of file GEMDigiMatcher.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), and eostools::move().

                                                                                      {
   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>();
 }

GEMDigiMatcher::~GEMDigiMatcher ( )

inline

Definition at line 37 of file GEMDigiMatcher.h.

37 {}

Member Function Documentation

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

inline

Definition at line 83 of file GEMDigiMatcher.h.

References detid_to_pads_.

83 { return detid_to_pads_; }

GEMDigiMatcher::detid_to_pads_

std::map< unsigned int, GEMPadDigiContainer > detid_to_pads_

Definition: GEMDigiMatcher.h:173

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

Definition at line 341 of file GEMDigiMatcher.cc.

                                                                          {
   return selectDetIds(chamber_to_clusters_, gem_type);
 }

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

Definition at line 333 of file GEMDigiMatcher.cc.

                                                                       {
   return selectDetIds(chamber_to_digis_, gem_type);
 }

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

Definition at line 337 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(chamber_to_pads_, gem_type);
 }

void GEMDigiMatcher::clear ( void )

private

Definition at line 538 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();
 }

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

Definition at line 403 of file GEMDigiMatcher.cc.

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

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

Definition at line 397 of file GEMDigiMatcher.cc.

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

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

Definition at line 409 of file GEMDigiMatcher.cc.

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

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

Definition at line 415 of file GEMDigiMatcher.cc.

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

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

Definition at line 329 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(detid_to_clusters_, gem_type);
 }

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

Definition at line 323 of file GEMDigiMatcher.cc.

                                                                   {
   return selectDetIds(detid_to_digis_, gem_type);
 }

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

Definition at line 327 of file GEMDigiMatcher.cc.

327 { return selectDetIds(detid_to_pads_, gem_type); }

GEMDigiMatcher::selectDetIds

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

Definition: GEMDigiMatcher.h:191

GEMDigiMatcher::detid_to_pads_

std::map< unsigned int, GEMPadDigiContainer > detid_to_pads_

Definition: GEMDigiMatcher.h:173

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

Definition at line 319 of file GEMDigiMatcher.cc.

                                                                      {
   return selectDetIds(detid_to_simLinks_, gem_type);
 }

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

Definition at line 367 of file GEMDigiMatcher.cc.

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

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

Definition at line 361 of file GEMDigiMatcher.cc.

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

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

Definition at line 373 of file GEMDigiMatcher.cc.

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

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

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

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

Definition at line 524 of file GEMDigiMatcher.cc.

References GEMDigi::strip().

                                                                                          {
   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;
 }

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

Definition at line 531 of file GEMDigiMatcher.cc.

References GEMPadDigi::pad().

                                                                                             {
   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;
 }

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

Definition at line 47 of file GEMDigiMatcher.cc.

References edm::Event::getByToken(), and edm::EventSetup::getHandle().

                                                                              {
   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_);
 
   const auto gemH = iSetup.getHandle(geomToken_);
   if (!gemH.isValid()) {
     gemGeometry_ = nullptr;
     edm::LogError("GEMDigiMatcher") << "Failed to initialize GEM geometry.";
   }
   gemGeometry_ = gemH.product();
 }

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

do the matching

Definition at line 66 of file GEMDigiMatcher.cc.

References funct::abs(), clear(), and CoreSimTrack::momentum().

                                                                 {
   // 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);
 }

void GEMDigiMatcher::matchClustersToSimTrack ( const GEMPadDigiClusterCollection & clusters )

private

Definition at line 245 of file GEMDigiMatcher.cc.

References GEMDetId::chamberId(), trackerHitRTTI::isMatched(), AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), and GEMDetId::superChamberId().

                                                                                         {
   for (auto it = clusters.begin(); it != clusters.end(); ++it) {
     const GEMDetId p_id = (*it).first;
     const auto clvec = (*it).second;
     for (auto cluster = clvec.first; cluster != clvec.second; ++cluster) {
       bool isMatched = false;
 
       // 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 (const auto& pad : detid_to_pads_[p_id.rawId()]) {
           if (pad.pad() == p) {
             isMatched = true;
           }
         }
       }
       if (isMatched) {
         detid_to_clusters_[p_id.rawId()].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;
       }
     }
   }
 }

void GEMDigiMatcher::matchCoPadsToSimTrack ( const GEMCoPadDigiCollection & co_pads )

private

Definition at line 278 of file GEMDigiMatcher.cc.

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

                                                                                 {
   // 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;
       }
     }
   }
 }

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

private

Definition at line 93 of file GEMDigiMatcher.cc.

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

                                                                                           {
   if (verboseSimLink_)
     edm::LogInfo("GEMDigiMatcher") << "Matching simtrack to GEM simlinks" << endl;
 
   // loop on the simlinks
   for (auto itsimlink = digisSL.begin(); itsimlink != digisSL.end(); itsimlink++) {
     GEMDetId p_id(itsimlink->id);
     for (auto sl = itsimlink->data.begin(); sl != itsimlink->data.end(); ++sl) {
       // 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->getTrackId() << 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.trackId() == sl->getTrackId() and simhit.particleType() == sl->getParticleType()) {
           detid_to_simLinks_[p_id.rawId()].push_back(*sl);
           if (verboseSimLink_)
             edm::LogInfo("GEMDigiMatcher") << "...was matched!" << endl;
           break;
         }
       }
     }
   }
 }

void GEMDigiMatcher::matchDigisToSimTrack ( const GEMDigiCollection & digis )

private

Definition at line 136 of file GEMDigiMatcher.cc.

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

                                                                         {
   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;
       }
     }
   }
 }

void GEMDigiMatcher::matchPadsInDetId	(	const GEMDetId &	simhitId,
		const GEMPadDigiCollection &
	)

private

void GEMDigiMatcher::matchPadsToSimTrack ( const GEMPadDigiCollection & pads )

private

Definition at line 181 of file GEMDigiMatcher.cc.

References GEMDetId::chamber(), GEMDetId::chamberId(), GEMPadDigi::GE21, GEMPadDigi::GE21SplitStrip, GEMDetId::isGE11(), GEMDetId::isGE21(), GEMDetId::isME0(), GEMDetId::layer(), or, DetId::rawId(), GEMDetId::region(), GEMDetId::ring(), GEMDetId::roll(), GEMDetId::station(), and GEMDetId::superChamberId().

                                                                          {
   for (auto it = pads.begin(); it != pads.end(); ++it) {
     const GEMDetId& p_id = (*it).first;
     const auto& padsvec = (*it).second;
 
     for (auto pad = padsvec.first; pad != padsvec.second; ++pad) {
       // 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;
 
       auto digivec = detid_to_digis_[p_id.rawId()];
       /*
         For GE1/1 and ME0 (and obsolete 8-partition GE2/1) geometries, the matching
         is pretty simple. Each simhit is converted to a digi. Two digis in neighboring
         strips are ed together into a pad. So for these geometries you just need
         to match pads to simtracks that are in the same eta partition (detid) as the
         simhit is in.  By convention, all pads in the 16-partition GE2/1 geometry are
         assigned to odd partition numbers. For the 16-partition GE2/1 geometry, you may
         have a track with simhits in eta partition 1 and 2 on the same strip and only
         produce 1 pad associated to eta partition 1. Therefore, for pads with odd
         partition number N, you need to consider the digis in the even partition number
         N+1.
       */
       if (p_id.roll() % 2 == 1 && p_id.isGE21() && pad->nPartitions() == GEMPadDigi::GE21SplitStrip) {
         // make the GEMDetId for the neighboring partition
         GEMDetId p_id2(p_id.region(), p_id.ring(), p_id.station(), p_id.layer(), p_id.chamber(), p_id.roll() + 1);
 
         // make a temporary container for its digis
         auto digivec2 = detid_to_digis_[p_id2.rawId()];
 
         // now add it to the container for the digis in the odd partition number
         digivec.insert(digivec.end(), digivec2.begin(), digivec2.end());
       }
       for (const auto& digi : digivec) {
         // for 8-partition geometries, the pad number equals the strip number divided by two
         const bool match8Partition(digi.strip() / 2 == pad->pad());
 
         // for 16-partition geometries, the pad number is the strip number itself
         const bool match16Partition(digi.strip() == pad->pad());
 
         // now consider the different cases separately
         const bool matchGE0(p_id.isME0() and match8Partition);
         const bool matchGE11(p_id.isGE11() and match8Partition);
         const bool matchGE21_8(p_id.isGE21() and pad->nPartitions() == GEMPadDigi::GE21 and match8Partition);
         const bool matchGE21_16(p_id.isGE21() and pad->nPartitions() == GEMPadDigi::GE21SplitStrip and
                                 match16Partition);
 
         // OR them together
         if (matchGE0 or matchGE11 or matchGE21_8 or matchGE21_16) {
           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;
         }
       }
     }
   }
 }

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

inline

Definition at line 110 of file GEMDigiMatcher.h.

References muonSimHitMatcher_.

110 { return muonSimHitMatcher_; }

GEMDigiMatcher::muonSimHitMatcher_

std::shared_ptr< GEMSimHitMatcher > muonSimHitMatcher_

Definition: GEMDigiMatcher.h:140

int GEMDigiMatcher::nCoPads ( ) const

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

Definition at line 475 of file GEMDigiMatcher.cc.

References dqmiodumpmetadata::n.

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

int GEMDigiMatcher::nLayersWithClustersInSuperChamber ( unsigned int detid ) const

Definition at line 451 of file GEMDigiMatcher.cc.

References GEMDetId::chamber(), HLT_FULL_cff::clusters, LayerTriplets::layers(), GEMDetId::region(), GEMDetId::ring(), and GEMDetId::station().

                                                                               {
   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();
 }

int GEMDigiMatcher::nLayersWithDigisInSuperChamber ( unsigned int detid ) const

Definition at line 421 of file GEMDigiMatcher.cc.

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

                                                                            {
   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();
 }

int GEMDigiMatcher::nLayersWithPadsInSuperChamber ( unsigned int detid ) const

Definition at line 436 of file GEMDigiMatcher.cc.

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

                                                                           {
   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();
 }

int GEMDigiMatcher::nPads ( ) const

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

Definition at line 466 of file GEMDigiMatcher.cc.

References dqmiodumpmetadata::n.

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

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

Definition at line 493 of file GEMDigiMatcher.cc.

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

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

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

Definition at line 385 of file GEMDigiMatcher.cc.

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

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

Definition at line 379 of file GEMDigiMatcher.cc.

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

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

Definition at line 391 of file GEMDigiMatcher.cc.

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

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

Definition at line 502 of file GEMDigiMatcher.cc.

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

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

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

Definition at line 513 of file GEMDigiMatcher.cc.

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

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

template<class T >

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

private

Definition at line 191 of file GEMDigiMatcher.h.

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

                                                                                     {
   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;
 }

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

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

Definition at line 484 of file GEMDigiMatcher.cc.

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

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

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

Definition at line 353 of file GEMDigiMatcher.cc.

                                                                               {
   return selectDetIds(superchamber_to_clusters_, gem_type);
 }

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

Definition at line 357 of file GEMDigiMatcher.cc.

                                                                             {
   return selectDetIds(superchamber_to_copads_, gem_type);
 }

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

Definition at line 345 of file GEMDigiMatcher.cc.

                                                                            {
   return selectDetIds(superchamber_to_digis_, gem_type);
 }

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

Definition at line 349 of file GEMDigiMatcher.cc.

                                                                           {
   return selectDetIds(superchamber_to_pads_, gem_type);
 }

Member Data Documentation

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

private

Definition at line 178 of file GEMDigiMatcher.h.

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

private

Definition at line 170 of file GEMDigiMatcher.h.

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

private

Definition at line 174 of file GEMDigiMatcher.h.

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

private

Definition at line 177 of file GEMDigiMatcher.h.

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

private

Definition at line 169 of file GEMDigiMatcher.h.

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

private

Definition at line 173 of file GEMDigiMatcher.h.

Referenced by allPads().

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

private

Definition at line 167 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::discardEleHits_

private

Definition at line 149 of file GEMDigiMatcher.h.

edm::Handle<GEMPadDigiClusterCollection> GEMDigiMatcher::gemClustersH_

private

Definition at line 137 of file GEMDigiMatcher.h.

edm::EDGetTokenT<GEMPadDigiClusterCollection> GEMDigiMatcher::gemClusterToken_

private

Definition at line 131 of file GEMDigiMatcher.h.

edm::Handle<GEMCoPadDigiCollection> GEMDigiMatcher::gemCoPadsH_

private

Definition at line 138 of file GEMDigiMatcher.h.

edm::EDGetTokenT<GEMCoPadDigiCollection> GEMDigiMatcher::gemCoPadToken_

private

Definition at line 132 of file GEMDigiMatcher.h.

edm::Handle<GEMDigiCollection> GEMDigiMatcher::gemDigisH_

private

Definition at line 135 of file GEMDigiMatcher.h.

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

private

Definition at line 134 of file GEMDigiMatcher.h.

edm::EDGetTokenT<GEMDigiCollection> GEMDigiMatcher::gemDigiToken_

private

Definition at line 129 of file GEMDigiMatcher.h.

const GEMGeometry* GEMDigiMatcher::gemGeometry_

private

Definition at line 143 of file GEMDigiMatcher.h.

edm::Handle<GEMPadDigiCollection> GEMDigiMatcher::gemPadsH_

private

Definition at line 136 of file GEMDigiMatcher.h.

edm::EDGetTokenT<GEMPadDigiCollection> GEMDigiMatcher::gemPadToken_

private

Definition at line 130 of file GEMDigiMatcher.h.

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

private

Definition at line 128 of file GEMDigiMatcher.h.

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

private

Definition at line 142 of file GEMDigiMatcher.h.

int GEMDigiMatcher::matchDeltaStrip_

private

Definition at line 156 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::matchToSimLink_

private

Definition at line 164 of file GEMDigiMatcher.h.

int GEMDigiMatcher::maxBXCluster_

private

Definition at line 153 of file GEMDigiMatcher.h.

int GEMDigiMatcher::maxBXCoPad_

private

Definition at line 154 of file GEMDigiMatcher.h.

int GEMDigiMatcher::maxBXDigi_

private

Definition at line 151 of file GEMDigiMatcher.h.

int GEMDigiMatcher::maxBXPad_

private

Definition at line 152 of file GEMDigiMatcher.h.

int GEMDigiMatcher::minBXCluster_

private

Definition at line 153 of file GEMDigiMatcher.h.

int GEMDigiMatcher::minBXCoPad_

private

Definition at line 154 of file GEMDigiMatcher.h.

int GEMDigiMatcher::minBXDigi_

private

Definition at line 151 of file GEMDigiMatcher.h.

int GEMDigiMatcher::minBXPad_

private

Definition at line 152 of file GEMDigiMatcher.h.

std::shared_ptr<GEMSimHitMatcher> GEMDigiMatcher::muonSimHitMatcher_

private

Definition at line 140 of file GEMDigiMatcher.h.

Referenced by muonSimHitMatcher().

GEMPadDigiClusterContainer GEMDigiMatcher::no_gem_clusters_

private

Definition at line 186 of file GEMDigiMatcher.h.

GEMCoPadDigiContainer GEMDigiMatcher::no_gem_copads_

private

Definition at line 187 of file GEMDigiMatcher.h.

GEMDigiContainer GEMDigiMatcher::no_gem_digis_

private

Definition at line 184 of file GEMDigiMatcher.h.

GEMPadDigiContainer GEMDigiMatcher::no_gem_pads_

private

Definition at line 185 of file GEMDigiMatcher.h.

GEMDigiSimLinkContainer GEMDigiMatcher::no_gem_simLinks_

private

Definition at line 183 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::simMuOnly_

private

Definition at line 148 of file GEMDigiMatcher.h.

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

private

Definition at line 179 of file GEMDigiMatcher.h.

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

private

Definition at line 181 of file GEMDigiMatcher.h.

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

private

Definition at line 171 of file GEMDigiMatcher.h.

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

private

Definition at line 175 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::verboseCluster_

private

Definition at line 161 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::verboseCoPad_

private

Definition at line 162 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::verboseDigi_

private

Definition at line 159 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::verbosePad_

private

Definition at line 160 of file GEMDigiMatcher.h.

bool GEMDigiMatcher::verboseSimLink_

private

Definition at line 158 of file GEMDigiMatcher.h.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation