#include <RPCSimHitMatcher.h>

Inheritance diagram for RPCSimHitMatcher:

Public Member Functions
std::set< unsigned int >	chamberIds (int type=MuonHitHelper::RPC_ALL) const

std::set< unsigned int >	detIds (int type=MuonHitHelper::RPC_ALL) const

bool	hitStation (int st) const

std::set< int >	hitStripsInDetId (unsigned int, int margin_n_strips=0) const

void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
	initialize the event More...

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

int	nStations () const

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

float	simHitsMeanStrip (const edm::PSimHitContainer &sim_hits) const

	~RPCSimHitMatcher ()

Public Member Functions inherited from MuonSimHitMatcher
std::set< unsigned int >	chamberIds (int type=0) const

std::set< unsigned int >	detIds (int type=0) const

const TrackingGeometry *	geometry ()

const edm::PSimHitContainer &	hitsInChamber (unsigned int) const

const edm::PSimHitContainer &	hitsInDetId (unsigned int) const

void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
	initialize the event More...

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

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

const edm::PSimHitContainer &	simHits (int) const

GlobalVector	simHitsMeanMomentum (const edm::PSimHitContainer &sim_hits) const

GlobalPoint	simHitsMeanPosition (const edm::PSimHitContainer &sim_hits) const

GlobalPoint	simHitsMeanPositionStation (int n) const

	~MuonSimHitMatcher ()

Private Member Functions
void	matchSimHitsToSimTrack ()

Private Attributes
edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	geomToken_

Additional Inherited Members
Protected Member Functions inherited from MuonSimHitMatcher
void	clear ()

std::vector< unsigned int >	getIdsOfSimTrackShower (unsigned trk_id, const edm::SimTrackContainer &simTracks, const edm::SimVertexContainer &simVertices)

Protected Attributes inherited from MuonSimHitMatcher
std::map< unsigned int, edm::PSimHitContainer >	chamber_to_hits_

std::map< unsigned int, edm::PSimHitContainer >	detid_to_hits_

bool	discardEleHits_

const TrackingGeometry *	geometry_

bool	hasGeometry_

edm::PSimHitContainer	hits_

edm::PSimHitContainer	no_hits_

edm::EDGetTokenT< edm::PSimHitContainer >	simHitInput_

edm::ParameterSet	simHitPSet_

edm::PSimHitContainer	simHits_

edm::Handle< edm::PSimHitContainer >	simHitsH_

bool	simMuOnly_

edm::EDGetTokenT< edm::SimTrackContainer >	simTrackInput_

edm::SimTrackContainer	simTracks_

edm::Handle< edm::SimTrackContainer >	simTracksH_

edm::EDGetTokenT< edm::SimVertexContainer >	simVertexInput_

edm::SimVertexContainer	simVertices_

edm::Handle< edm::SimVertexContainer >	simVerticesH_

std::vector< unsigned >	track_ids_

std::map< unsigned int, unsigned int >	trkid_to_index_

bool	verbose_

bool	verboseSimTrack_

Detailed Description

Description: Matching of RPC SimHit to SimTrack

Author: Sven Dildick (TAMU), Tao Huang (TAMU)

Definition at line 16 of file RPCSimHitMatcher.h.

Constructor & Destructor Documentation

◆ RPCSimHitMatcher()

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

Definition at line 6 of file RPCSimHitMatcher.cc.

References MuonSimHitMatcher::discardEleHits_, geomToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), MuonSimHitMatcher::simHitInput_, MuonSimHitMatcher::simHitPSet_, MuonSimHitMatcher::simMuOnly_, and MuonSimHitMatcher::verbose_.

     : MuonSimHitMatcher(ps, std::move(iC)) {
   simHitPSet_ = ps.getParameterSet("rpcSimHit");
   verbose_ = simHitPSet_.getParameter<int>("verbose");
   simMuOnly_ = simHitPSet_.getParameter<bool>("simMuOnly");
   discardEleHits_ = simHitPSet_.getParameter<bool>("discardEleHits");
 
   simHitInput_ = iC.consumes<edm::PSimHitContainer>(simHitPSet_.getParameter<edm::InputTag>("inputTag"));
   geomToken_ = iC.esConsumes<RPCGeometry, MuonGeometryRecord>();
 }

◆ ~RPCSimHitMatcher()

RPCSimHitMatcher::~RPCSimHitMatcher ( )

inline

Definition at line 22 of file RPCSimHitMatcher.h.

22 {}

Member Function Documentation

◆ chamberIds()

std::set< unsigned int > RPCSimHitMatcher::chamberIds ( int type = MuonHitHelper::RPC_ALL ) const

Definition at line 85 of file RPCSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, hcalRecHitTable_cff::detId, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonHitHelper::toRPCType().

Referenced by hitStation(), and match().

                                                                 {
   std::set<unsigned int> result;
   for (const auto& p : chamber_to_hits_) {
     const auto& id = p.first;
     if (type > 0) {
       RPCDetId detId(id);
       if (MuonHitHelper::toRPCType(detId.region(), detId.station(), detId.ring()) != type)
         continue;
     }
     result.insert(id);
   }
   return result;
 }

◆ detIds()

std::set< unsigned int > RPCSimHitMatcher::detIds ( int type = MuonHitHelper::RPC_ALL ) const

Definition at line 71 of file RPCSimHitMatcher.cc.

References hcalRecHitTable_cff::detId, MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonHitHelper::toRPCType().

Referenced by match().

                                                             {
   std::set<unsigned int> result;
   for (const auto& p : detid_to_hits_) {
     const auto& id = p.first;
     if (type > 0) {
       RPCDetId detId(id);
       if (MuonHitHelper::toRPCType(detId.region(), detId.station(), detId.ring()) != type)
         continue;
     }
     result.insert(id);
   }
   return result;
 }

◆ hitStation()

bool RPCSimHitMatcher::hitStation ( int st ) const

Definition at line 99 of file RPCSimHitMatcher.cc.

References chamberIds(), l1ctLayer2EG_cff::id, and relativeConstraints::station.

Referenced by nStations().

                                               {
   int nst = 0;
   for (const auto& ddt : chamberIds(0)) {
     const RPCDetId id(ddt);
     if (id.station() != st)
       continue;
     ++nst;
   }
   return nst;
 }

◆ hitStripsInDetId()

std::set< int > RPCSimHitMatcher::hitStripsInDetId	(	unsigned int	detid,
		int	margin_n_strips = `0`
	)		const

Definition at line 129 of file RPCSimHitMatcher.cc.

References relativeConstraints::chamber, MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), l1ctLayer2EG_cff::id, mps_fire::result, nano_mu_digi_cff::roll, and contentValuesCheck::ss.

Referenced by match().

                                                                                             {
   set<int> result;
   RPCDetId id(detid);
   for (const auto& roll : dynamic_cast<const RPCGeometry*>(geometry_)->chamber(id)->rolls()) {
     int max_nstrips = roll->nstrips();
     for (const auto& h : MuonSimHitMatcher::hitsInDetId(roll->id().rawId())) {
       const LocalPoint& lp = h.entryPoint();
       int central_strip = static_cast<int>(roll->topology().channel(lp));
       int smin = central_strip - margin_n_strips;
       smin = (smin > 0) ? smin : 1;
       int smax = central_strip + margin_n_strips;
       smax = (smax <= max_nstrips) ? smax : max_nstrips;
       for (int ss = smin; ss <= smax; ++ss)
         result.insert(ss);
     }
   }
   return result;
 }

◆ init()

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

initialize the event

Definition at line 18 of file RPCSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, geomToken_, edm::EventSetup::getData(), iEvent, and MuonSimHitMatcher::init().

                                                                                {
   geometry_ = &iSetup.getData(geomToken_);
   MuonSimHitMatcher::init(iEvent, iSetup);
 }

◆ match()

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

do the matching

Definition at line 24 of file RPCSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, chamberIds(), detIds(), MuonSimHitMatcher::hitsInChamber(), hitStripsInDetId(), l1ctLayer2EG_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), AlCaHLTBitMon_ParallelJobs::p, HLTBitAnalyser_cfi::simhits, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simHitsMeanPosition(), DigiDM_cff::strips, HLT_2024v10_cff::track, MuonSimHitMatcher::track_ids_, MuonSimHitMatcher::verbose_, and bphysicsOniaDQM_cfi::vertex.

                                                                            {
   // instantiates the track ids and simhits
   MuonSimHitMatcher::match(track, vertex);
 
   matchSimHitsToSimTrack();
 
   if (verbose_) {
     edm::LogInfo("RPCSimHitMatcher") << "nSimHits " << simHits_.size() << " nTrackIds " << track_ids_.size();
     edm::LogInfo("RPCSimHitMatcher") << "detids RPC " << detIds().size();
 
     const auto& ch_ids = chamberIds();
     for (const auto& id : ch_ids) {
       const auto& simhits = MuonSimHitMatcher::hitsInChamber(id);
       const auto& simhits_gp = simHitsMeanPosition(simhits);
       edm::LogInfo("RPCSimHitMatcher") << "RPCDetId " << RPCDetId(id) << ": nHits " << simhits.size() << " eta "
                                        << simhits_gp.eta() << " phi " << simhits_gp.phi() << " nCh "
                                        << chamber_to_hits_[id].size();
       const auto& strips = hitStripsInDetId(id);
       edm::LogInfo("RPCSimHitMatcher") << "nStrips " << strips.size();
       edm::LogInfo("RPCSimHitMatcher") << "strips : ";
       for (const auto& p : strips) {
         edm::LogInfo("RPCSimHitMatcher") << p;
       }
     }
   }
 }

◆ matchSimHitsToSimTrack()

void RPCSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 51 of file RPCSimHitMatcher.cc.

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, MuonSimHitMatcher::detid_to_hits_, MuonSimHitMatcher::discardEleHits_, h, MuonSimHitMatcher::hits_, EgammaValidation_cff::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, and MuonSimHitMatcher::track_ids_.

Referenced by match().

                                               {
   for (const auto& track_id : track_ids_) {
     for (const auto& h : simHits_) {
       if (h.trackId() != track_id)
         continue;
       int pdgid = h.particleType();
       if (simMuOnly_ && std::abs(pdgid) != 13)
         continue;
       // discard electron hits in the RPC chambers
       if (discardEleHits_ && std::abs(pdgid) == 11)
         continue;
 
       const RPCDetId& layer_id(h.detUnitId());
       detid_to_hits_[h.detUnitId()].push_back(h);
       hits_.push_back(h);
       chamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
     }
   }
 }

◆ nStations()

int RPCSimHitMatcher::nStations ( ) const

Definition at line 110 of file RPCSimHitMatcher.cc.

References hitStation().

110 { return (hitStation(1) + hitStation(2) + hitStation(3) + hitStation(4)); }

RPCSimHitMatcher::hitStation

bool hitStation(int st) const

Definition: RPCSimHitMatcher.cc:99

◆ simHitsMeanStrip()

float RPCSimHitMatcher::simHitsMeanStrip ( const edm::PSimHitContainer & sim_hits ) const

Definition at line 112 of file RPCSimHitMatcher.cc.

References ztail::d, f, MuonSimHitMatcher::geometry_, h, dqmiodumpmetadata::n, and nano_mu_digi_cff::roll.

                                                                                   {
   if (sim_hits.empty())
     return -1.f;
 
   float sums = 0.f;
   size_t n = 0;
   for (const auto& h : sim_hits) {
     const LocalPoint& lp = h.entryPoint();
     const auto& d = h.detUnitId();
     sums += dynamic_cast<const RPCGeometry*>(geometry_)->roll(d)->strip(lp);
     ++n;
   }
   if (n == 0)
     return -1.f;
   return sums / n;
 }

Member Data Documentation

◆ geomToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> RPCSimHitMatcher::geomToken_

private

Definition at line 49 of file RPCSimHitMatcher.h.

Referenced by init(), and RPCSimHitMatcher().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RPCSimHitMatcher()

◆ ~RPCSimHitMatcher()

Member Function Documentation

◆ chamberIds()

◆ detIds()

◆ hitStation()

◆ hitStripsInDetId()

◆ init()

◆ match()

◆ matchSimHitsToSimTrack()

◆ nStations()

◆ simHitsMeanStrip()

Member Data Documentation

◆ geomToken_