#include <CSCSimHitMatcher.h>

Inheritance diagram for CSCSimHitMatcher:

Public Member Functions
void	camberIdsToString (const std::set< unsigned int > &) const

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

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

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

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

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

void	fitHitsInChamber (unsigned int detid, float &mean, float &slope) const

bool	hitStation (int, int) const

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

std::set< int >	hitWiregroupsInDetId (unsigned int, int margin_n_wg=0) const

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

float	LocalBendingInChamber (unsigned int detid) const

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

int	nCoincidenceChambers (int min_n_layers=4) const

int	nLayersWithHitsInChamber (unsigned int) const

int	nStations (int nl=4) const

GlobalPoint	simHitPositionKeyLayer (unsigned int chamberid) const

GlobalPoint	simHitsMeanPositionStation (int n) const

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

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

	~CSCSimHitMatcher ()

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	clear ()

void	matchSimHitsToSimTrack ()

Private Attributes
edm::ESGetToken< CSCGeometry, 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 CSC SimHit to SimTrack

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

Definition at line 16 of file CSCSimHitMatcher.h.

Constructor & Destructor Documentation

◆ CSCSimHitMatcher()

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

Definition at line 8 of file CSCSimHitMatcher.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("cscSimHit");
   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<CSCGeometry, MuonGeometryRecord>();
 
   //In order to make fitting ROOT histograms thread safe
   // one must call this undocumented function
   TMinuitMinimizer::UseStaticMinuit(false);
 }

◆ ~CSCSimHitMatcher()

CSCSimHitMatcher::~CSCSimHitMatcher ( )

inline

Definition at line 22 of file CSCSimHitMatcher.h.

22 {}

Member Function Documentation

◆ camberIdsToString()

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

◆ chamberIds()

std::set< unsigned int > CSCSimHitMatcher::chamberIds ( int type = MuonHitHelper::CSC_ALL ) const

Definition at line 100 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, CSCDetId::ring(), CSCDetId::station(), and MuonHitHelper::toCSCType().

Referenced by chamberIdsStation(), hitStation(), and nCoincidenceChambers().

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

◆ chamberIdsStation()

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

Definition at line 349 of file CSCSimHitMatcher.cc.

References chamberIds(), MuonHitHelper::CSC_ME12, MuonHitHelper::CSC_ME13, MuonHitHelper::CSC_ME1a, MuonHitHelper::CSC_ME1b, MuonHitHelper::CSC_ME21, MuonHitHelper::CSC_ME22, MuonHitHelper::CSC_ME31, MuonHitHelper::CSC_ME32, MuonHitHelper::CSC_ME41, MuonHitHelper::CSC_ME42, LaserDQM_cfg::p1, SiStripOfflineCRack_cfg::p2, chargedHadronTrackResolutionFilter_cfi::p3, mps_fire::result, and relativeConstraints::station.

                                                                           {
   set<unsigned int> result;
   switch (station) {
     case 1: {
       const auto& p1(chamberIds(MuonHitHelper::CSC_ME1a));
       const auto& p2(chamberIds(MuonHitHelper::CSC_ME1b));
       const auto& p3(chamberIds(MuonHitHelper::CSC_ME12));
       const auto& p4(chamberIds(MuonHitHelper::CSC_ME13));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       result.insert(p3.begin(), p3.end());
       result.insert(p4.begin(), p4.end());
       break;
     }
     case 2: {
       const auto& p1(chamberIds(MuonHitHelper::CSC_ME21));
       const auto& p2(chamberIds(MuonHitHelper::CSC_ME22));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       break;
     }
     case 3: {
       const auto& p1(chamberIds(MuonHitHelper::CSC_ME31));
       const auto& p2(chamberIds(MuonHitHelper::CSC_ME32));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       break;
     }
     case 4: {
       const auto& p1(chamberIds(MuonHitHelper::CSC_ME41));
       const auto& p2(chamberIds(MuonHitHelper::CSC_ME42));
       result.insert(p1.begin(), p1.end());
       result.insert(p2.begin(), p2.end());
       break;
     }
   };
   return result;
 }

◆ chamberIdsToString()

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

Definition at line 342 of file CSCSimHitMatcher.cc.

References AlCaHLTBitMon_ParallelJobs::p.

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

◆ clear()

void CSCSimHitMatcher::clear ( void )

private

Definition at line 388 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::clear().

Referenced by match().

388 { MuonSimHitMatcher::clear(); }

MuonSimHitMatcher::clear

void clear()

Definition: MuonSimHitMatcher.cc:139

◆ detIds()

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

Definition at line 86 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, CSCDetId::ring(), CSCDetId::station(), and MuonHitHelper::toCSCType().

Referenced by match().

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

◆ fitHitsInChamber()

void CSCSimHitMatcher::fitHitsInChamber	(	unsigned int	detid,
		float &	mean,
		float &	slope
	)		const

Definition at line 215 of file CSCSimHitMatcher.cc.

References ztail::d, l1tTowerCalibrationProducer_cfi::fit, MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInChamber(), phase1PixelTopology::layer, alignCSCRings::s, slope, x, and y.

                                                                                                 {
   const CSCDetId cscid(detid);
 
   const auto& sim_hits = hitsInChamber(detid);
 
   if (sim_hits.empty())
     return;
 
   vector<float> x;
   vector<float> y;
   vector<float> xe;
   vector<float> ye;
 
   const float HALF_STRIP_ERROR = 0.288675;
 
   for (const auto& h : sim_hits) {
     const LocalPoint& lp = h.entryPoint();
     const auto& d = h.detUnitId();
     float s = dynamic_cast<const CSCGeometry*>(geometry_)->layer(d)->geometry()->strip(lp);
     // shift to key half strip layer (layer 3)
     x.push_back(CSCDetId(d).layer() - 3);
     y.push_back(s);
     xe.push_back(float(0));
     ye.push_back(2 * HALF_STRIP_ERROR);
   }
   if (x.size() < 2)
     return;
 
   std::unique_ptr<TGraphErrors> gr(new TGraphErrors(x.size(), &x[0], &y[0], &xe[0], &ye[0]));
   TF1 fit("fit", "pol1", -3, 4);
   gr->Fit(&fit, "EMQN");
 
   intercept = fit.GetParameter(0);
   slope = fit.GetParameter(1);
 }

◆ hitStation()

bool CSCSimHitMatcher::hitStation	(	int	st,
		int	nlayers
	)		const

Definition at line 123 of file CSCSimHitMatcher.cc.

References relativeConstraints::chamber, chamberIds(), makeMuonMisalignmentScenario::endcap, l1ctLayer2EG_cff::id, nlayers, nLayersWithHitsInChamber(), DetId::rawId(), relativeConstraints::ring, and relativeConstraints::station.

Referenced by nStations().

                                                            {
   int nst = 0;
   for (const auto& ddt : chamberIds()) {
     const CSCDetId id(ddt);
     int ri(id.ring());
     if (id.station() != st)
       continue;
 
     // ME1/a case - check the ME1/b chamber
     if (st == 1 and ri == 4) {
       CSCDetId idME1b(id.endcap(), id.station(), 1, id.chamber(), 0);
       const int nl1a(nLayersWithHitsInChamber(id.rawId()));
       const int nl1b(nLayersWithHitsInChamber(idME1b.rawId()));
       if (nl1a + nl1b < nlayers)
         continue;
       ++nst;
     }
     // ME1/b case - check the ME1/a chamber
     else if (st == 1 and ri == 1) {
       CSCDetId idME1a(id.endcap(), id.station(), 4, id.chamber(), 0);
       const int nl1a(nLayersWithHitsInChamber(idME1a.rawId()));
       const int nl1b(nLayersWithHitsInChamber(id.rawId()));
       if (nl1a + nl1b < nlayers)
         continue;
       ++nst;
     }
     // default case
     else {
       const int nl(nLayersWithHitsInChamber(id.rawId()));
       if (nl < nlayers)
         continue;
       ++nst;
     }
   }
   return nst;
 }

◆ hitStripsInDetId()

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

Definition at line 295 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), l1ctLayer2EG_cff::id, phase1PixelTopology::layer, mps_fire::result, HLTBitAnalyser_cfi::simhits, and contentValuesCheck::ss.

Referenced by match().

                                                                                             {
   set<int> result;
   const auto& simhits = MuonSimHitMatcher::hitsInDetId(detid);
   CSCDetId id(detid);
   int max_nstrips = dynamic_cast<const CSCGeometry*>(geometry_)->layer(id)->geometry()->numberOfStrips();
   for (const auto& h : simhits) {
     const LocalPoint& lp = h.entryPoint();
     int central_strip = dynamic_cast<const CSCGeometry*>(geometry_)->layer(id)->geometry()->nearestStrip(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;
 }

◆ hitWiregroupsInDetId()

std::set< int > CSCSimHitMatcher::hitWiregroupsInDetId	(	unsigned int	detid,
		int	margin_n_wg = `0`
	)		const

Definition at line 313 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), l1ctLayer2EG_cff::id, phase1PixelTopology::layer, mps_fire::result, and HLTBitAnalyser_cfi::simhits.

                                                                                             {
   set<int> result;
   const auto& simhits = MuonSimHitMatcher::hitsInDetId(detid);
   CSCDetId id(detid);
   int max_n_wg = dynamic_cast<const CSCGeometry*>(geometry_)->layer(id)->geometry()->numberOfWireGroups();
   for (const auto& h : simhits) {
     const LocalPoint& lp = h.entryPoint();
     const auto& layer_geo = dynamic_cast<const CSCGeometry*>(geometry_)->layer(id)->geometry();
     int central_wg = layer_geo->wireGroup(layer_geo->nearestWire(lp));
     int wg_min = central_wg - margin_n_wg;
     wg_min = (wg_min > 0) ? wg_min : 1;
     int wg_max = central_wg + margin_n_wg;
     wg_max = (wg_max <= max_n_wg) ? wg_max : max_n_wg;
     for (int wg = wg_min; wg <= wg_max; ++wg)
       result.insert(wg);
   }
   return result;
 }

◆ init()

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

initialize the event

Definition at line 24 of file CSCSimHitMatcher.cc.

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

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

◆ LocalBendingInChamber()

float CSCSimHitMatcher::LocalBendingInChamber ( unsigned int detid ) const

Definition at line 164 of file CSCSimHitMatcher.cc.

References CSCDetId::chamber(), SiPixelRawToDigiRegional_cfi::deltaPhi, CSCDetId::endcap(), MuonSimHitMatcher::hitsInDetId(), nLayersWithHitsInChamber(), or, PV3DBase< T, PVType, FrameType >::phi(), CSCDetId::ring(), MuonSimHitMatcher::simHitsMeanPosition(), and CSCDetId::station().

                                                                       {
   const CSCDetId cscid(detid);
   if (nLayersWithHitsInChamber(detid) < 6)
     return -100;
   float phi_layer1 = -10;
   float phi_layer6 = 10;
 
   if (cscid.station() == 1 and (cscid.ring() == 1 or cscid.ring() == 4)) {
     // phi in layer 1
     const CSCDetId cscid1a(cscid.endcap(), cscid.station(), 4, cscid.chamber(), 1);
     const CSCDetId cscid1b(cscid.endcap(), cscid.station(), 1, cscid.chamber(), 1);
     const edm::PSimHitContainer& hits1a = MuonSimHitMatcher::hitsInDetId(cscid1a.rawId());
     const edm::PSimHitContainer& hits1b = MuonSimHitMatcher::hitsInDetId(cscid1b.rawId());
     const GlobalPoint& gp1a = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1a.rawId()));
     const GlobalPoint& gp1b = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1b.rawId()));
     if (!hits1a.empty() and !hits1b.empty())
       phi_layer1 = (gp1a.phi() + gp1b.phi()) / 2.0;
     else if (!hits1a.empty())
       phi_layer1 = gp1a.phi();
     else if (!hits1b.empty())
       phi_layer1 = gp1b.phi();
 
     // phi in layer 6
     const CSCDetId cscid6a(cscid.endcap(), cscid.station(), 4, cscid.chamber(), 6);
     const CSCDetId cscid6b(cscid.endcap(), cscid.station(), 1, cscid.chamber(), 6);
     const edm::PSimHitContainer& hits6a = MuonSimHitMatcher::hitsInDetId(cscid6a.rawId());
     const edm::PSimHitContainer& hits6b = MuonSimHitMatcher::hitsInDetId(cscid6b.rawId());
     const GlobalPoint& gp6a = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6a.rawId()));
     const GlobalPoint& gp6b = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6b.rawId()));
     if (!hits6a.empty() and !hits6b.empty())
       phi_layer6 = (gp6a.phi() + gp6b.phi()) / 2.0;
     else if (!hits6a.empty())
       phi_layer6 = gp6a.phi();
     else if (!hits6b.empty())
       phi_layer6 = gp6b.phi();
 
   } else {
     // phi in layer 1
     const CSCDetId cscid1(cscid.endcap(), cscid.station(), cscid.ring(), cscid.chamber(), 1);
     const GlobalPoint& gp1 = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1.rawId()));
     phi_layer1 = gp1.phi();
 
     // phi in layer 6
     const CSCDetId cscid6(cscid.endcap(), cscid.station(), cscid.ring(), cscid.chamber(), 6);
     const GlobalPoint& gp6 = simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6.rawId()));
     phi_layer6 = gp6.phi();
   }
   return deltaPhi(phi_layer6, phi_layer1);
 }

◆ match()

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

do the matching

Definition at line 30 of file CSCSimHitMatcher.cc.

References funct::abs(), clear(), MuonSimHitMatcher::detid_to_hits_, detIds(), MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), hitStripsInDetId(), l1ctLayer2EG_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), or, AlCaHLTBitMon_ParallelJobs::p, CSCDetId::ring(), HLTBitAnalyser_cfi::simhits, MuonSimHitMatcher::simHitsMeanPosition(), CSCDetId::station(), DigiDM_cff::strips, HLT_2022v14_cff::track, MuonSimHitMatcher::track_ids_, MuonSimHitMatcher::verbose_, and bphysicsOniaDQM_cfi::vertex.

                                                                            {
   clear();
 
   // instantiates the track ids and simhits
   MuonSimHitMatcher::match(track, vertex);
 
   // hard cut on non-CSC muons
   if (std::abs(track.momentum().eta()) < 0.9)
     return;
   if (std::abs(track.momentum().eta()) > 2.45)
     return;
 
   matchSimHitsToSimTrack();
 
   if (verbose_) {
     edm::LogInfo("CSCSimHitMatcher") << "nTrackIds " << track_ids_.size() << " nSelectedCSCSimHits " << hits_.size();
     edm::LogInfo("CSCSimHitMatcher") << "detids CSC " << detIds(0).size();
 
     for (const auto& id : detIds(0)) {
       const auto& simhits = hitsInDetId(id);
       const auto& simhits_gp = simHitsMeanPosition(simhits);
       const auto& strips = hitStripsInDetId(id);
       CSCDetId cscid(id);
       if (cscid.station() == 1 and (cscid.ring() == 1 or cscid.ring() == 4)) {
         edm::LogInfo("CSCSimHitMatcher") << "cscdetid " << CSCDetId(id) << ": " << simhits.size() << " "
                                          << simhits_gp.phi() << " " << detid_to_hits_[id].size();
         edm::LogInfo("CSCSimHitMatcher") << "nStrip " << strips.size();
         edm::LogInfo("CSCSimHitMatcher") << "strips : ";
         for (const auto& p : strips) {
           edm::LogInfo("CSCSimHitMatcher") << p;
         }
       }
     }
   }
 }

◆ matchSimHitsToSimTrack()

void CSCSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 66 of file CSCSimHitMatcher.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 CSC chambers
       if (discardEleHits_ && std::abs(pdgid) == 11)
         continue;
 
       const CSCDetId& layer_id(h.detUnitId());
       hits_.push_back(h);
       detid_to_hits_[h.detUnitId()].push_back(h);
       chamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
     }
   }
 }

◆ nCoincidenceChambers()

int CSCSimHitMatcher::nCoincidenceChambers ( int min_n_layers = 4 ) const

Definition at line 332 of file CSCSimHitMatcher.cc.

References chamberIds(), nLayersWithHitsInChamber(), and mps_fire::result.

                                                                  {
   int result = 0;
   const auto& chamber_ids = chamberIds(0);
   for (const auto& id : chamber_ids) {
     if (nLayersWithHitsInChamber(id) >= min_n_layers)
       result += 1;
   }
   return result;
 }

◆ nLayersWithHitsInChamber()

int CSCSimHitMatcher::nLayersWithHitsInChamber ( unsigned int detid ) const

Definition at line 114 of file CSCSimHitMatcher.cc.

References h, and MuonSimHitMatcher::hitsInChamber().

Referenced by hitStation(), LocalBendingInChamber(), and nCoincidenceChambers().

                                                                        {
   set<int> layers_with_hits;
   for (const auto& h : MuonSimHitMatcher::hitsInChamber(detid)) {
     const CSCDetId& idd(h.detUnitId());
     layers_with_hits.insert(idd.layer());
   }
   return layers_with_hits.size();
 }

◆ nStations()

int CSCSimHitMatcher::nStations ( int nl = 4 ) const

Definition at line 160 of file CSCSimHitMatcher.cc.

References hitStation(), and nlayers.

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

◆ simHitPositionKeyLayer()

GlobalPoint CSCSimHitMatcher::simHitPositionKeyLayer ( unsigned int chamberid ) const

◆ simHitsMeanPositionStation()

GlobalPoint CSCSimHitMatcher::simHitsMeanPositionStation ( int n ) const

◆ simHitsMeanStrip()

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

Definition at line 251 of file CSCSimHitMatcher.cc.

References ztail::d, f, MuonSimHitMatcher::geometry_, h, phase1PixelTopology::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();
     s = dynamic_cast<const CSCGeometry*>(geometry_)->layer(d)->geometry()->strip(lp);
     // convert to half-strip:
     s *= 2.;
     sums += s;
     ++n;
   }
   if (n == 0)
     return -1.f;
   return sums / n;
 }

◆ simHitsMeanWG()

float CSCSimHitMatcher::simHitsMeanWG ( const edm::PSimHitContainer & sim_hits ) const

Definition at line 272 of file CSCSimHitMatcher.cc.

References ztail::d, f, MuonSimHitMatcher::geometry(), MuonSimHitMatcher::geometry_, h, phase1PixelTopology::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();
     // find nearest wire
     const auto& layerG(dynamic_cast<const CSCGeometry*>(geometry_)->layer(d)->geometry());
     int nearestWire(layerG->nearestWire(lp));
     // then find the corresponding wire group
     s = layerG->wireGroup(nearestWire);
     sums += s;
     ++n;
   }
   if (n == 0)
     return -1.f;
   return sums / n;
 }

Member Data Documentation

◆ geomToken_

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> CSCSimHitMatcher::geomToken_

private

Definition at line 81 of file CSCSimHitMatcher.h.

Referenced by CSCSimHitMatcher(), and init().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CSCSimHitMatcher()

◆ ~CSCSimHitMatcher()

Member Function Documentation

◆ camberIdsToString()

◆ chamberIds()

◆ chamberIdsStation()

◆ chamberIdsToString()

◆ clear()

◆ detIds()

◆ fitHitsInChamber()

◆ hitStation()

◆ hitStripsInDetId()

◆ hitWiregroupsInDetId()

◆ init()

◆ LocalBendingInChamber()

◆ match()

◆ matchSimHitsToSimTrack()

◆ nCoincidenceChambers()

◆ nLayersWithHitsInChamber()

◆ nStations()

◆ simHitPositionKeyLayer()

◆ simHitsMeanPositionStation()

◆ simHitsMeanStrip()

◆ simHitsMeanWG()

Member Data Documentation

◆ geomToken_