#include <DTSimHitMatcher.h>

Inheritance diagram for DTSimHitMatcher:

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

std::set< unsigned int >	chamberIdsStation (int station) const

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

void	dtChamberIdsToString (const std::set< unsigned int > &) const

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

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

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

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

bool	hitStation (int, int, int) const

std::set< unsigned int >	hitWiresInDTChamberId (unsigned int, int margin_n_wires=0) const

std::set< unsigned int >	hitWiresInDTLayerId (unsigned int, int margin_n_wires=0) const

std::set< unsigned int >	hitWiresInDTSuperLayerId (unsigned int, int margin_n_wires=0) const

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

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

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

int	nCellsWithHitsInLayer (unsigned int) const

int	nLayersWithHitsInChamber (unsigned int) const

int	nLayersWithHitsInSuperLayer (unsigned int) const

int	nStations (int nsl=1, int nl=3) const

int	nSuperLayersWithHitsInChamber (unsigned int) const

GlobalPoint	simHitsMeanPositionStation (int n) const

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

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

	~DTSimHitMatcher ()

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< DTGeometry, MuonGeometryRecord >	geomToken_

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

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

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 DT SimHit to SimTrack

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

Definition at line 16 of file DTSimHitMatcher.h.

Constructor & Destructor Documentation

◆ DTSimHitMatcher()

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

Definition at line 5 of file DTSimHitMatcher.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("dtSimHit");
   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<DTGeometry, MuonGeometryRecord>();
 }

◆ ~DTSimHitMatcher()

DTSimHitMatcher::~DTSimHitMatcher ( )

inline

Definition at line 22 of file DTSimHitMatcher.h.

22 {}

Member Function Documentation

◆ chamberIds()

std::set< unsigned int > DTSimHitMatcher::chamberIds ( int type = MuonHitHelper::DT_ALL ) const

Definition at line 80 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, DTChamberId::station(), MuonHitHelper::toDTType(), and DTChamberId::wheel().

Referenced by chamberIdsStation(), and hitStation().

                                                                   {
   std::set<unsigned int> result;
   for (const auto& p : chamber_to_hits_) {
     const auto& id = p.first;
     if (dt_type > 0) {
       DTChamberId detId(id);
       if (MuonHitHelper::toDTType(detId.wheel(), detId.station()) != dt_type)
         continue;
     }
     result.insert(id);
   }
   return result;
 }

◆ chamberIdsStation()

std::set< unsigned int > DTSimHitMatcher::chamberIdsStation ( int station ) const

Definition at line 288 of file DTSimHitMatcher.cc.

                                                                          {
   set<unsigned int> result;
   switch (station) {
     case 1: {
       const auto& p1(chamberIds(MuonHitHelper::DT_MB21p));
       const auto& p2(chamberIds(MuonHitHelper::DT_MB11p));
       const auto& p3(chamberIds(MuonHitHelper::DT_MB01));
       const auto& p4(chamberIds(MuonHitHelper::DT_MB11n));
       const auto& p5(chamberIds(MuonHitHelper::DT_MB21n));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       result.insert(p3.begin(), p3.end());
       result.insert(p4.begin(), p4.end());
       result.insert(p5.begin(), p5.end());
       break;
     }
     case 2: {
       const auto& p1(chamberIds(MuonHitHelper::DT_MB22p));
       const auto& p2(chamberIds(MuonHitHelper::DT_MB12p));
       const auto& p3(chamberIds(MuonHitHelper::DT_MB02));
       const auto& p4(chamberIds(MuonHitHelper::DT_MB12n));
       const auto& p5(chamberIds(MuonHitHelper::DT_MB22n));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       result.insert(p3.begin(), p3.end());
       result.insert(p4.begin(), p4.end());
       result.insert(p5.begin(), p5.end());
       break;
     }
     case 3: {
       const auto& p1(chamberIds(MuonHitHelper::DT_MB23p));
       const auto& p2(chamberIds(MuonHitHelper::DT_MB13p));
       const auto& p3(chamberIds(MuonHitHelper::DT_MB03));
       const auto& p4(chamberIds(MuonHitHelper::DT_MB13n));
       const auto& p5(chamberIds(MuonHitHelper::DT_MB23n));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       result.insert(p3.begin(), p3.end());
       result.insert(p4.begin(), p4.end());
       result.insert(p5.begin(), p5.end());
       break;
     }
     case 4: {
       const auto& p1(chamberIds(MuonHitHelper::DT_MB24p));
       const auto& p2(chamberIds(MuonHitHelper::DT_MB14p));
       const auto& p3(chamberIds(MuonHitHelper::DT_MB04));
       const auto& p4(chamberIds(MuonHitHelper::DT_MB14n));
       const auto& p5(chamberIds(MuonHitHelper::DT_MB24n));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       result.insert(p3.begin(), p3.end());
       result.insert(p4.begin(), p4.end());
       result.insert(p5.begin(), p5.end());
       break;
     }
   };
   return result;
 }

◆ detIds()

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

Definition at line 66 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, DTChamberId::station(), MuonHitHelper::toDTType(), and DTChamberId::wheel().

Referenced by match().

                                                               {
   std::set<unsigned int> result;
   for (const auto& p : detid_to_hits_) {
     const auto& id = p.first;
     if (dt_type > 0) {
       DTWireId detId(id);
       if (MuonHitHelper::toDTType(detId.wheel(), detId.station()) != dt_type)
         continue;
     }
     result.insert(id);
   }
   return result;
 }

◆ dtChamberIdsToString()

void DTSimHitMatcher::dtChamberIdsToString ( const std::set< unsigned int > & set ) const

Definition at line 281 of file DTSimHitMatcher.cc.

References AlCaHLTBitMon_ParallelJobs::p.

                                                                                 {
   for (const auto& p : set) {
     DTChamberId detId(p);
     edm::LogInfo("DTSimHitMatcher") << " " << detId << "\n";
   }
 }

◆ hitsInChamber()

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

Definition at line 128 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, l1ctLayer2EG_cff::id, MuonHitHelper::isDT(), and MuonSimHitMatcher::no_hits_.

                                                                                   {
   if (!MuonHitHelper::isDT(detid))
     return no_hits_;
 
   const DTWireId id(detid);
   if (chamber_to_hits_.find(id.chamberId().rawId()) == chamber_to_hits_.end())
     return no_hits_;
   return chamber_to_hits_.at(id.chamberId().rawId());
 }

◆ hitsInLayer()

const edm::PSimHitContainer & DTSimHitMatcher::hitsInLayer ( unsigned int detid ) const

Definition at line 108 of file DTSimHitMatcher.cc.

References l1ctLayer2EG_cff::id, MuonHitHelper::isDT(), layer_to_hits_, and MuonSimHitMatcher::no_hits_.

Referenced by nCellsWithHitsInLayer().

                                                                                 {
   if (!MuonHitHelper::isDT(detid))
     return no_hits_;
 
   const DTWireId id(detid);
   if (layer_to_hits_.find(id.layerId().rawId()) == layer_to_hits_.end())
     return no_hits_;
   return layer_to_hits_.at(id.layerId().rawId());
 }

◆ hitsInSuperLayer()

const edm::PSimHitContainer & DTSimHitMatcher::hitsInSuperLayer ( unsigned int detid ) const

Definition at line 118 of file DTSimHitMatcher.cc.

References l1ctLayer2EG_cff::id, MuonHitHelper::isDT(), MuonSimHitMatcher::no_hits_, and superlayer_to_hits_.

Referenced by nLayersWithHitsInSuperLayer().

                                                                                      {
   if (!MuonHitHelper::isDT(detid))
     return no_hits_;
 
   const DTWireId id(detid);
   if (superlayer_to_hits_.find(id.superlayerId().rawId()) == superlayer_to_hits_.end())
     return no_hits_;
   return superlayer_to_hits_.at(id.superlayerId().rawId());
 }

◆ hitStation()

bool DTSimHitMatcher::hitStation	(	int	st,
		int	nsuperlayers,
		int	nlayers
	)		const

Definition at line 138 of file DTSimHitMatcher.cc.

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

Referenced by nStations().

                                                                             {
   int nst = 0;
   for (const auto& ddt : chamberIds()) {
     const DTChamberId id(ddt);
     if (id.station() != st)
       continue;
 
     // require at least 1 superlayer
     const int nsl(nSuperLayersWithHitsInChamber(id.rawId()));
     if (nsl < nsuperlayers)
       continue;
 
     // require at least 3 layers hit per chamber
     const int nl(nLayersWithHitsInChamber(id.rawId()));
     if (nl < nlayers)
       continue;
     ++nst;
   }
   return nst;
 }

◆ hitWiresInDTChamberId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTChamberId	(	unsigned int	detid,
		int	margin_n_wires = `0`
	)		const

Definition at line 269 of file DTSimHitMatcher.cc.

References relativeConstraints::chamber, MuonSimHitMatcher::geometry_, hitWiresInDTSuperLayerId(), AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonSimHitMatcher::verbose_.

                                                                                                         {
   set<unsigned int> result;
   const auto& superLayers(dynamic_cast<const DTGeometry*>(geometry_)->chamber(DTChamberId(detid))->superLayers());
   for (const auto& sl : superLayers) {
     if (verbose_)
       edm::LogInfo("DTSimHitMatcher") << "hitWiresInDTChamberId::sl id " << sl->id();
     const auto& p(hitWiresInDTSuperLayerId(sl->id().rawId(), margin_n_wires));
     result.insert(p.begin(), p.end());
   }
   return result;
 }

◆ hitWiresInDTLayerId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTLayerId	(	unsigned int	detid,
		int	margin_n_wires = `0`
	)		const

Definition at line 231 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), l1ctLayer2EG_cff::id, MuonHitHelper::isDT(), pixelTopology::layer, DetId::rawId(), mps_fire::result, contentValuesCheck::ss, and MuonSimHitMatcher::verbose_.

Referenced by hitWiresInDTSuperLayerId().

                                                                                                       {
   set<unsigned int> result;
   if (MuonHitHelper::isDT(detid)) {
     DTLayerId id(detid);
     int max_nwires = dynamic_cast<const DTGeometry*>(geometry_)->layer(id)->specificTopology().channels();
     for (int wn = 0; wn <= max_nwires; ++wn) {
       DTWireId wid(id, wn);
       for (const auto& h : MuonSimHitMatcher::hitsInDetId(wid.rawId())) {
         if (verbose_)
           edm::LogInfo("DTSimHitMatcher") << "central DTWireId " << wid << " simhit " << h;
         int smin = wn - margin_n_wires;
         smin = (smin > 0) ? smin : 1;
         int smax = wn + margin_n_wires;
         smax = (smax <= max_nwires) ? smax : max_nwires;
         for (int ss = smin; ss <= smax; ++ss) {
           DTWireId widd(id, ss);
           if (verbose_)
             edm::LogInfo("DTSimHitMatcher") << "\tadding DTWireId to collection " << widd;
           result.insert(widd.rawId());
         }
       }
     }
   }
   return result;
 }

◆ hitWiresInDTSuperLayerId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTSuperLayerId	(	unsigned int	detid,
		int	margin_n_wires = `0`
	)		const

Definition at line 257 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, hitWiresInDTLayerId(), cmsLHEtoEOSManager::l, hgcalTBTopologyTester_cfi::layers, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonSimHitMatcher::verbose_.

Referenced by hitWiresInDTChamberId().

                                                                                                            {
   set<unsigned int> result;
   const auto& layers(dynamic_cast<const DTGeometry*>(geometry_)->superLayer(DTSuperLayerId(detid))->layers());
   for (const auto& l : layers) {
     if (verbose_)
       edm::LogInfo("DTSimHitMatcher") << "hitWiresInDTSuperLayerId::l id " << l->id();
     const auto& p(hitWiresInDTLayerId(l->id().rawId(), margin_n_wires));
     result.insert(p.begin(), p.end());
   }
   return result;
 }

◆ init()

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

initialize the event

Definition at line 17 of file DTSimHitMatcher.cc.

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

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

◆ layerIds()

std::set< unsigned int > DTSimHitMatcher::layerIds ( ) const

Definition at line 94 of file DTSimHitMatcher.cc.

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

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

◆ match()

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

do the matching

Definition at line 23 of file DTSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, detIds(), MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), l1ctLayer2EG_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), MuonSimHitMatcher::simHitsMeanPosition(), HLT_2022v15_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("DTSimHitMatcher") << "nTrackIds " << track_ids_.size() << " nSelectedDTSimHits " << hits_.size();
     edm::LogInfo("DTSimHitMatcher") << "detids DT " << detIds(0).size();
 
     const auto& dt_det_ids = detIds(0);
     for (const auto& id : dt_det_ids) {
       const auto& dt_simhits = MuonSimHitMatcher::hitsInDetId(id);
       const auto& dt_simhits_gp = simHitsMeanPosition(dt_simhits);
       edm::LogInfo("DTSimHitMatcher") << "DTWireId " << DTWireId(id) << ": nHits " << dt_simhits.size() << " eta "
                                       << dt_simhits_gp.eta() << " phi " << dt_simhits_gp.phi() << " nCh "
                                       << chamber_to_hits_[id].size();
     }
   }
 }

◆ matchSimHitsToSimTrack()

void DTSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 44 of file DTSimHitMatcher.cc.

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, MuonSimHitMatcher::detid_to_hits_, MuonSimHitMatcher::discardEleHits_, h, MuonSimHitMatcher::hits_, layer_to_hits_, EgammaValidation_cff::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, superlayer_to_hits_, 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 DT chambers
       if (discardEleHits_ && std::abs(pdgid) == 11)
         continue;
 
       const DTWireId layer_id(h.detUnitId());
       detid_to_hits_[h.detUnitId()].push_back(h);
       hits_.push_back(h);
       layer_to_hits_[layer_id.layerId().rawId()].push_back(h);
       superlayer_to_hits_[layer_id.superlayerId().rawId()].push_back(h);
       chamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
     }
   }
 }

◆ nCellsWithHitsInLayer()

int DTSimHitMatcher::nCellsWithHitsInLayer ( unsigned int detid ) const

Definition at line 164 of file DTSimHitMatcher.cc.

References h, hfClusterShapes_cfi::hits, hitsInLayer(), and MuonHitHelper::isDT().

                                                                    {
   set<int> layers_with_hits;
   const auto& hits = hitsInLayer(detid);
   for (const auto& h : hits) {
     if (MuonHitHelper::isDT(detid)) {
       const DTWireId idd(h.detUnitId());
       layers_with_hits.insert(idd.wire());
     }
   }
   return layers_with_hits.size();
 }

◆ nLayersWithHitsInChamber()

int DTSimHitMatcher::nLayersWithHitsInChamber ( unsigned int detid ) const

Definition at line 200 of file DTSimHitMatcher.cc.

References relativeConstraints::chamber, MuonSimHitMatcher::geometry_, MuonTCMETValueMapProducer_cff::nLayers, and nLayersWithHitsInSuperLayer().

Referenced by hitStation().

                                                                       {
   int nLayers = 0;
   const auto& superLayers(dynamic_cast<const DTGeometry*>(geometry_)->chamber(DTChamberId(detid))->superLayers());
   for (const auto& sl : superLayers) {
     nLayers += nLayersWithHitsInSuperLayer(sl->id().rawId());
   }
   return nLayers;
 }

◆ nLayersWithHitsInSuperLayer()

int DTSimHitMatcher::nLayersWithHitsInSuperLayer ( unsigned int detid ) const

Definition at line 176 of file DTSimHitMatcher.cc.

References h, hfClusterShapes_cfi::hits, hitsInSuperLayer(), and MuonHitHelper::isDT().

Referenced by nLayersWithHitsInChamber().

                                                                          {
   set<int> layers_with_hits;
   const auto& hits = hitsInSuperLayer(detid);
   for (const auto& h : hits) {
     if (MuonHitHelper::isDT(detid)) {
       const DTLayerId idd(h.detUnitId());
       layers_with_hits.insert(idd.layer());
     }
   }
   return layers_with_hits.size();
 }

◆ nStations()

int DTSimHitMatcher::nStations	(	int	nsl = `1`,
		int	nl = `3`
	)		const

Definition at line 159 of file DTSimHitMatcher.cc.

References hitStation(), and nlayers.

                                                                   {
   return (hitStation(1, nsuperlayers, nlayers) + hitStation(2, nsuperlayers, nlayers) +
           hitStation(3, nsuperlayers, nlayers) + hitStation(4, nsuperlayers, nlayers));
 }

◆ nSuperLayersWithHitsInChamber()

int DTSimHitMatcher::nSuperLayersWithHitsInChamber ( unsigned int detid ) const

Definition at line 188 of file DTSimHitMatcher.cc.

References h, hfClusterShapes_cfi::hits, MuonSimHitMatcher::hitsInChamber(), and MuonHitHelper::isDT().

Referenced by hitStation().

                                                                            {
   set<int> sl_with_hits;
   const auto& hits = MuonSimHitMatcher::hitsInChamber(detid);
   for (const auto& h : hits) {
     if (MuonHitHelper::isDT(detid)) {
       const DTSuperLayerId idd(h.detUnitId());
       sl_with_hits.insert(idd.superLayer());
     }
   }
   return sl_with_hits.size();
 }

◆ simHitsMeanPositionStation()

GlobalPoint DTSimHitMatcher::simHitsMeanPositionStation ( int n ) const

◆ simHitsMeanWire()

float DTSimHitMatcher::simHitsMeanWire ( const edm::PSimHitContainer & sim_hits ) const

Definition at line 208 of file DTSimHitMatcher.cc.

References ztail::d, f, MuonSimHitMatcher::geometry_, h, MuonHitHelper::isDT(), pixelTopology::layer, dqmiodumpmetadata::n, and alignCSCRings::s.

                                                                                 {
   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();
     float s;
     const auto& d = h.detUnitId();
     if (MuonHitHelper::isDT(d)) {
       // find nearest wire
       s = dynamic_cast<const DTGeometry*>(geometry_)->layer(DTLayerId(d))->specificTopology().channel(lp);
     } else
       continue;
     sums += s;
     ++n;
   }
   if (n == 0)
     return -1.f;
   return sums / n;
 }

◆ superlayerIds()

std::set< unsigned int > DTSimHitMatcher::superlayerIds ( ) const

Definition at line 101 of file DTSimHitMatcher.cc.

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

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

Member Data Documentation

◆ geomToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTSimHitMatcher::geomToken_

private

Definition at line 78 of file DTSimHitMatcher.h.

Referenced by DTSimHitMatcher(), and init().

◆ layer_to_hits_

std::map<unsigned int, edm::PSimHitContainer> DTSimHitMatcher::layer_to_hits_

private

Definition at line 75 of file DTSimHitMatcher.h.

Referenced by hitsInLayer(), layerIds(), and matchSimHitsToSimTrack().

◆ superlayer_to_hits_

std::map<unsigned int, edm::PSimHitContainer> DTSimHitMatcher::superlayer_to_hits_

private

Definition at line 76 of file DTSimHitMatcher.h.

Referenced by hitsInSuperLayer(), matchSimHitsToSimTrack(), and superlayerIds().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTSimHitMatcher()

◆ ~DTSimHitMatcher()

Member Function Documentation

◆ chamberIds()

◆ chamberIdsStation()

◆ detIds()

◆ dtChamberIdsToString()

◆ hitsInChamber()

◆ hitsInLayer()

◆ hitsInSuperLayer()

◆ hitStation()

◆ hitWiresInDTChamberId()

◆ hitWiresInDTLayerId()

◆ hitWiresInDTSuperLayerId()

◆ init()

◆ layerIds()

◆ match()

◆ matchSimHitsToSimTrack()

◆ nCellsWithHitsInLayer()

◆ nLayersWithHitsInChamber()

◆ nLayersWithHitsInSuperLayer()

◆ nStations()

◆ nSuperLayersWithHitsInChamber()

◆ simHitsMeanPositionStation()

◆ simHitsMeanWire()

◆ superlayerIds()

Member Data Documentation

◆ geomToken_

◆ layer_to_hits_

◆ superlayer_to_hits_