#include <MuonSimHitMatcher.h>

Inheritance diagram for MuonSimHitMatcher:

Public Member Functions
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 ()

Protected Member Functions
void	clear ()

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

Protected Attributes
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 muon SimHit to SimTrack

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

Definition at line 32 of file MuonSimHitMatcher.h.

Constructor & Destructor Documentation

◆ MuonSimHitMatcher()

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

Definition at line 7 of file MuonSimHitMatcher.cc.

References edm::ParameterSet::getParameterSet(), cscDigiValidation_cfi::simTrack, and ecalDetailedTimeRecHit_cfi::simVertex.

                                                                                          {
   const auto& simVertex = ps.getParameterSet("simVertex");
   const auto& simTrack = ps.getParameterSet("simTrack");
   verboseSimTrack_ = simTrack.getParameter<int>("verbose");
 
   simVertexInput_ = iC.consumes<edm::SimVertexContainer>(simVertex.getParameter<edm::InputTag>("inputTag"));
   simTrackInput_ = iC.consumes<edm::SimTrackContainer>(simTrack.getParameter<edm::InputTag>("inputTag"));
 }

◆ ~MuonSimHitMatcher()

MuonSimHitMatcher::~MuonSimHitMatcher ( )

inline

Definition at line 38 of file MuonSimHitMatcher.h.

38 {}

Member Function Documentation

◆ chamberIds()

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

◆ clear()

void MuonSimHitMatcher::clear ( void )

protected

Definition at line 139 of file MuonSimHitMatcher.cc.

Referenced by CSCSimHitMatcher::clear(), and GEMSimHitMatcher::clear().

                               {
   track_ids_.clear();
   trkid_to_index_.clear();
   hits_.clear();
   detid_to_hits_.clear();
   chamber_to_hits_.clear();
 }

◆ detIds()

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

◆ geometry()

const TrackingGeometry* MuonSimHitMatcher::geometry ( )

inline

Definition at line 69 of file MuonSimHitMatcher.h.

References geometry_.

Referenced by CSCSimHitMatcher::simHitsMeanWG().

69 { return geometry_; }

MuonSimHitMatcher::geometry_

const TrackingGeometry * geometry_

Definition: MuonSimHitMatcher.h:84

◆ getIdsOfSimTrackShower()

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

protected

Definition at line 47 of file MuonSimHitMatcher.cc.

References SimTrack::noVertex(), mps_fire::result, submitPVValidationJobs::t, and SimTrack::vertIndex().

                                                                                                               {
   vector<unsigned int> result;
   result.push_back(initial_trk_id);
 
   if (!simMuOnly_)
     return result;
 
   for (const auto& t : simTracks_) {
     SimTrack last_trk = t;
     // if (std::abs(t.type()) != 13) continue;
     bool is_child = false;
     while (true) {
       if (last_trk.noVertex())
         break;
       if (simVertices_[last_trk.vertIndex()].noParent())
         break;
 
       unsigned parentId = simVertices_[last_trk.vertIndex()].parentIndex();
       if (parentId == initial_trk_id) {
         is_child = true;
         break;
       }
 
       const auto& association = trkid_to_index_.find(parentId);
       if (association == trkid_to_index_.end())
         break;
 
       last_trk = simTracks_[association->second];
     }
     if (is_child) {
       result.push_back(t.trackId());
     }
   }
   return result;
 }

◆ hitsInChamber()

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

Definition at line 93 of file MuonSimHitMatcher.cc.

Referenced by CSCSimHitMatcher::fitHitsInChamber(), RPCSimHitMatcher::match(), ME0SimHitMatcher::match(), CSCSimHitMatcher::nLayersWithHitsInChamber(), and DTSimHitMatcher::nSuperLayersWithHitsInChamber().

                                                                                     {
   if (chamber_to_hits_.find(detid) == chamber_to_hits_.end())
     return no_hits_;
   return chamber_to_hits_.at(detid);
 }

◆ hitsInDetId()

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

Definition at line 87 of file MuonSimHitMatcher.cc.

Referenced by RPCSimHitMatcher::hitStripsInDetId(), ME0SimHitMatcher::hitStripsInDetId(), GEMSimHitMatcher::hitStripsInDetId(), CSCSimHitMatcher::hitStripsInDetId(), CSCSimHitMatcher::hitWiregroupsInDetId(), DTSimHitMatcher::hitWiresInDTLayerId(), CSCSimHitMatcher::LocalBendingInChamber(), CSCSimHitMatcher::match(), DTSimHitMatcher::match(), GEMSimHitMatcher::match(), ME0SimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

                                                                                   {
   if (detid_to_hits_.find(detid) == detid_to_hits_.end())
     return no_hits_;
   return detid_to_hits_.at(detid);
 }

◆ init()

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

initialize the event

Definition at line 17 of file MuonSimHitMatcher.cc.

References iEvent.

Referenced by CSCSimHitMatcher::init(), RPCSimHitMatcher::init(), ME0SimHitMatcher::init(), GEMSimHitMatcher::init(), and DTSimHitMatcher::init().

                                                                                 {
   hasGeometry_ = true;
 
   iEvent.getByToken(simTrackInput_, simTracksH_);
   iEvent.getByToken(simVertexInput_, simVerticesH_);
   iEvent.getByToken(simHitInput_, simHitsH_);
 }

◆ match()

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

do the matching

Definition at line 26 of file MuonSimHitMatcher.cc.

References submitPVValidationJobs::t, and HLT_2023v12_cff::track.

Referenced by CSCSimHitMatcher::match(), RPCSimHitMatcher::match(), ME0SimHitMatcher::match(), GEMSimHitMatcher::match(), and DTSimHitMatcher::match().

                                                                             {
   simTracks_ = *simTracksH_.product();
   simVertices_ = *simVerticesH_.product();
   simHits_ = *simHitsH_.product();
 
   // fill trkId2Index association:
   int no = 0;
   trkid_to_index_.clear();
   for (const auto& t : simTracks_) {
     trkid_to_index_[t.trackId()] = no;
     no++;
   }
 
   track_ids_ = getIdsOfSimTrackShower(track.trackId(), simTracks_, simVertices_);
   if (verboseSimTrack_) {
     edm::LogInfo("MuonSimHitMatcher") << "Printing track_ids" << std::endl;
     for (const auto& id : track_ids_)
       edm::LogInfo("MuonSimHitMatcher") << "id: " << id << std::endl;
   }
 }

◆ simHits()

const edm::PSimHitContainer & MuonSimHitMatcher::simHits ( int sub ) const

Definition at line 85 of file MuonSimHitMatcher.cc.

Referenced by ntupleDataFormat.TrackingParticle::bestMatchingTrack().

85 { return hits_; }

MuonSimHitMatcher::hits_

edm::PSimHitContainer hits_

Definition: MuonSimHitMatcher.h:105

◆ simHitsMeanMomentum()

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

Definition at line 119 of file MuonSimHitMatcher.cc.

References h, dqmiodumpmetadata::n, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                              {
   if (sim_hits.empty())
     return GlobalVector();  // point "zero"
 
   float sumx, sumy, sumz;
   sumx = sumy = sumz = 0.f;
   size_t n = 0;
   for (const auto& h : sim_hits) {
     const LocalVector& lv = h.momentumAtEntry();
     const GlobalVector& gv = geometry_->idToDet(h.detUnitId())->surface().toGlobal(lv);
     sumx += gv.x();
     sumy += gv.y();
     sumz += gv.z();
     ++n;
   }
   if (n == 0)
     return GlobalVector();
   return GlobalVector(sumx / n, sumy / n, sumz / n);
 }

◆ simHitsMeanPosition()

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

Definition at line 99 of file MuonSimHitMatcher.cc.

References h, and dqmiodumpmetadata::n.

Referenced by CSCSimHitMatcher::LocalBendingInChamber(), RPCSimHitMatcher::match(), ME0SimHitMatcher::match(), GEMSimHitMatcher::match(), DTSimHitMatcher::match(), and CSCSimHitMatcher::match().

                                                                                             {
   if (sim_hits.empty())
     return GlobalPoint();  // point "zero"
 
   float sumx, sumy, sumz;
   sumx = sumy = sumz = 0.f;
   size_t n = 0;
   for (const auto& h : sim_hits) {
     const LocalPoint& lp = h.entryPoint();
     const GlobalPoint& gp = geometry_->idToDet(h.detUnitId())->surface().toGlobal(lp);
     sumx += gp.x();
     sumy += gp.y();
     sumz += gp.z();
     ++n;
   }
   if (n == 0)
     return GlobalPoint();
   return GlobalPoint(sumx / n, sumy / n, sumz / n);
 }

◆ simHitsMeanPositionStation()

GlobalPoint MuonSimHitMatcher::simHitsMeanPositionStation ( int n ) const

Member Data Documentation

◆ chamber_to_hits_

std::map<unsigned int, edm::PSimHitContainer> MuonSimHitMatcher::chamber_to_hits_

protected

Definition at line 107 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::chamberIds(), DTSimHitMatcher::chamberIds(), GEMSimHitMatcher::chamberIds(), ME0SimHitMatcher::chamberIds(), RPCSimHitMatcher::chamberIds(), DTSimHitMatcher::hitsInChamber(), GEMSimHitMatcher::match(), ME0SimHitMatcher::match(), DTSimHitMatcher::match(), RPCSimHitMatcher::match(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

◆ detid_to_hits_

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

protected

Definition at line 106 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::detIds(), DTSimHitMatcher::detIds(), ME0SimHitMatcher::detIds(), RPCSimHitMatcher::detIds(), GEMSimHitMatcher::detIds(), CSCSimHitMatcher::match(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

◆ discardEleHits_

bool MuonSimHitMatcher::discardEleHits_

protected

Definition at line 80 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::CSCSimHitMatcher(), DTSimHitMatcher::DTSimHitMatcher(), GEMSimHitMatcher::GEMSimHitMatcher(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), GEMSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::ME0SimHitMatcher(), and RPCSimHitMatcher::RPCSimHitMatcher().

◆ geometry_

const TrackingGeometry* MuonSimHitMatcher::geometry_

protected

Definition at line 84 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::fitHitsInChamber(), geometry(), RPCSimHitMatcher::hitStripsInDetId(), ME0SimHitMatcher::hitStripsInDetId(), GEMSimHitMatcher::hitStripsInDetId(), CSCSimHitMatcher::hitStripsInDetId(), CSCSimHitMatcher::hitWiregroupsInDetId(), DTSimHitMatcher::hitWiresInDTChamberId(), DTSimHitMatcher::hitWiresInDTLayerId(), DTSimHitMatcher::hitWiresInDTSuperLayerId(), DTSimHitMatcher::init(), ME0SimHitMatcher::init(), RPCSimHitMatcher::init(), CSCSimHitMatcher::init(), GEMSimHitMatcher::init(), ME0SimHitMatcher::matchSimHitsToSimTrack(), GEMSimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::nLayersWithHitsInChamber(), GEMSimHitMatcher::simHitsGEMCentralPosition(), RPCSimHitMatcher::simHitsMeanStrip(), ME0SimHitMatcher::simHitsMeanStrip(), CSCSimHitMatcher::simHitsMeanStrip(), GEMSimHitMatcher::simHitsMeanStrip(), CSCSimHitMatcher::simHitsMeanWG(), and DTSimHitMatcher::simHitsMeanWire().

◆ hasGeometry_

bool MuonSimHitMatcher::hasGeometry_

protected

Definition at line 82 of file MuonSimHitMatcher.h.

◆ hits_

edm::PSimHitContainer MuonSimHitMatcher::hits_

protected

Definition at line 105 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::match(), GEMSimHitMatcher::match(), DTSimHitMatcher::match(), ME0SimHitMatcher::match(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

◆ no_hits_

edm::PSimHitContainer MuonSimHitMatcher::no_hits_

protected

Definition at line 102 of file MuonSimHitMatcher.h.

Referenced by DTSimHitMatcher::hitsInChamber(), DTSimHitMatcher::hitsInLayer(), ME0SimHitMatcher::hitsInSuperChamber(), GEMSimHitMatcher::hitsInSuperChamber(), and DTSimHitMatcher::hitsInSuperLayer().

◆ simHitInput_

edm::EDGetTokenT<edm::PSimHitContainer> MuonSimHitMatcher::simHitInput_

protected

Definition at line 88 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::CSCSimHitMatcher(), DTSimHitMatcher::DTSimHitMatcher(), GEMSimHitMatcher::GEMSimHitMatcher(), ME0SimHitMatcher::ME0SimHitMatcher(), and RPCSimHitMatcher::RPCSimHitMatcher().

◆ simHitPSet_

edm::ParameterSet MuonSimHitMatcher::simHitPSet_

protected

Definition at line 109 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::CSCSimHitMatcher(), DTSimHitMatcher::DTSimHitMatcher(), GEMSimHitMatcher::GEMSimHitMatcher(), ME0SimHitMatcher::ME0SimHitMatcher(), and RPCSimHitMatcher::RPCSimHitMatcher().

◆ simHits_

edm::PSimHitContainer MuonSimHitMatcher::simHits_

protected

Definition at line 97 of file MuonSimHitMatcher.h.

Referenced by RPCSimHitMatcher::match(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

◆ simHitsH_

edm::Handle<edm::PSimHitContainer> MuonSimHitMatcher::simHitsH_

protected

Definition at line 92 of file MuonSimHitMatcher.h.

◆ simMuOnly_

bool MuonSimHitMatcher::simMuOnly_

protected

Definition at line 79 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::CSCSimHitMatcher(), DTSimHitMatcher::DTSimHitMatcher(), GEMSimHitMatcher::GEMSimHitMatcher(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), GEMSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::ME0SimHitMatcher(), and RPCSimHitMatcher::RPCSimHitMatcher().

◆ simTrackInput_

edm::EDGetTokenT<edm::SimTrackContainer> MuonSimHitMatcher::simTrackInput_

protected

Definition at line 87 of file MuonSimHitMatcher.h.

◆ simTracks_

edm::SimTrackContainer MuonSimHitMatcher::simTracks_

protected

Definition at line 94 of file MuonSimHitMatcher.h.

◆ simTracksH_

edm::Handle<edm::SimTrackContainer> MuonSimHitMatcher::simTracksH_

protected

Definition at line 90 of file MuonSimHitMatcher.h.

◆ simVertexInput_

edm::EDGetTokenT<edm::SimVertexContainer> MuonSimHitMatcher::simVertexInput_

protected

Definition at line 86 of file MuonSimHitMatcher.h.

◆ simVertices_

edm::SimVertexContainer MuonSimHitMatcher::simVertices_

protected

Definition at line 95 of file MuonSimHitMatcher.h.

◆ simVerticesH_

edm::Handle<edm::SimVertexContainer> MuonSimHitMatcher::simVerticesH_

protected

Definition at line 91 of file MuonSimHitMatcher.h.

◆ track_ids_

std::vector<unsigned> MuonSimHitMatcher::track_ids_

protected

Definition at line 99 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::match(), GEMSimHitMatcher::match(), RPCSimHitMatcher::match(), DTSimHitMatcher::match(), ME0SimHitMatcher::match(), RPCSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::matchSimHitsToSimTrack(), DTSimHitMatcher::matchSimHitsToSimTrack(), CSCSimHitMatcher::matchSimHitsToSimTrack(), and GEMSimHitMatcher::matchSimHitsToSimTrack().

◆ trkid_to_index_

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

protected

Definition at line 100 of file MuonSimHitMatcher.h.

◆ verbose_

bool MuonSimHitMatcher::verbose_

protected

Definition at line 81 of file MuonSimHitMatcher.h.

Referenced by CSCSimHitMatcher::CSCSimHitMatcher(), DTSimHitMatcher::DTSimHitMatcher(), GEMSimHitMatcher::GEMSimHitMatcher(), DTSimHitMatcher::hitWiresInDTChamberId(), DTSimHitMatcher::hitWiresInDTLayerId(), DTSimHitMatcher::hitWiresInDTSuperLayerId(), DTSimHitMatcher::match(), GEMSimHitMatcher::match(), RPCSimHitMatcher::match(), CSCSimHitMatcher::match(), ME0SimHitMatcher::match(), GEMSimHitMatcher::matchSimHitsToSimTrack(), ME0SimHitMatcher::ME0SimHitMatcher(), and RPCSimHitMatcher::RPCSimHitMatcher().

◆ verboseSimTrack_

bool MuonSimHitMatcher::verboseSimTrack_

protected

Definition at line 78 of file MuonSimHitMatcher.h.

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MuonSimHitMatcher()

◆ ~MuonSimHitMatcher()

Member Function Documentation

◆ chamberIds()

◆ clear()

◆ detIds()

◆ geometry()

◆ getIdsOfSimTrackShower()

◆ hitsInChamber()

◆ hitsInDetId()

◆ init()

◆ match()

◆ simHits()

◆ simHitsMeanMomentum()

◆ simHitsMeanPosition()

◆ simHitsMeanPositionStation()

Member Data Documentation

◆ chamber_to_hits_

◆ detid_to_hits_

◆ discardEleHits_

◆ geometry_

◆ hasGeometry_

◆ hits_

◆ no_hits_

◆ simHitInput_

◆ simHitPSet_

◆ simHits_

◆ simHitsH_

◆ simMuOnly_

◆ simTrackInput_

◆ simTracks_

◆ simTracksH_

◆ simVertexInput_

◆ simVertices_

◆ simVerticesH_

◆ track_ids_

◆ trkid_to_index_

◆ verbose_

◆ verboseSimTrack_