#include <ME0SimHitMatcher.h>

Inheritance diagram for ME0SimHitMatcher:

Public Member Functions
std::set< unsigned int >	chamberIds () const

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

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

std::set< int >	hitPartitions () const

const edm::PSimHitContainer &	hitsInSuperChamber (unsigned int) 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...

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

int	nCoincidenceChambers (int min_n_layers=4) const

int	nLayersWithHitsInSuperChamber (unsigned int) const

int	nPadsWithHits () const

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

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

std::set< unsigned int >	superChamberIdsCoincidences (int min_n_layers=4) const

	~ME0SimHitMatcher ()

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
std::map< unsigned int, std::set< int > >	detids_to_copads_

std::map< unsigned int, std::set< int > >	detids_to_pads_

edm::ESGetToken< ME0Geometry, MuonGeometryRecord >	geomToken_

std::map< unsigned int, edm::PSimHitContainer >	superChamber_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 ME0 SimHit to SimTrack

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

Definition at line 16 of file ME0SimHitMatcher.h.

Constructor & Destructor Documentation

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

Definition at line 5 of file ME0SimHitMatcher.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("me0SimHit");
   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<ME0Geometry, MuonGeometryRecord>();
 }

ME0SimHitMatcher::~ME0SimHitMatcher ( )

inline

Definition at line 22 of file ME0SimHitMatcher.h.

22 {}

Member Function Documentation

std::set< unsigned int > ME0SimHitMatcher::chamberIds ( ) const

Definition at line 92 of file ME0SimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, and mps_fire::result.

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

std::set< unsigned int > ME0SimHitMatcher::detIds ( void ) const

Definition at line 85 of file ME0SimHitMatcher.cc.

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, and mps_fire::result.

Referenced by hitPartitions(), match(), matchSimHitsToSimTrack(), and nPadsWithHits().

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

std::set< int > ME0SimHitMatcher::hitPadsInDetId ( unsigned int detid ) const

Definition at line 175 of file ME0SimHitMatcher.cc.

References detids_to_pads_, and none.

Referenced by nPadsWithHits().

                                                                      {
   set<int> none;
   if (detids_to_pads_.find(detid) == detids_to_pads_.end())
     return none;
   return detids_to_pads_.at(detid);
 }

std::set< int > ME0SimHitMatcher::hitPartitions ( ) const

Definition at line 182 of file ME0SimHitMatcher.cc.

References detIds(), mps_fire::result, and ME0DetId::roll().

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

const edm::PSimHitContainer & ME0SimHitMatcher::hitsInSuperChamber ( unsigned int detid ) const

Definition at line 106 of file ME0SimHitMatcher.cc.

References MuonSimHitMatcher::no_hits_, and superChamber_to_hits_.

Referenced by nLayersWithHitsInSuperChamber().

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

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

Definition at line 156 of file ME0SimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), gpuClustering::id, mps_fire::result, and contentValuesCheck::ss.

Referenced by match().

                                                                                             {
   set<int> result;
   const auto& simhits = hitsInDetId(detid);
   ME0DetId id(detid);
   int max_nstrips = dynamic_cast<const ME0Geometry*>(geometry_)->etaPartition(id)->nstrips();
   for (const auto& h : simhits) {
     const LocalPoint& lp = h.entryPoint();
     int central_strip =
         1 + static_cast<int>(dynamic_cast<const ME0Geometry*>(geometry_)->etaPartition(id)->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;
 }

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

initialize the event

Definition at line 17 of file ME0SimHitMatcher.cc.

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

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

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

do the matching

Definition at line 23 of file ME0SimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, detIds(), MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInChamber(), hitStripsInDetId(), gpuClustering::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), AlCaHLTBitMon_ParallelJobs::p, MuonSimHitMatcher::simHitsMeanPosition(), DigiDM_cff::strips, MuonSimHitMatcher::track_ids_, and MuonSimHitMatcher::verbose_.

                                                                            {
   // instantiates the track ids and simhits
   MuonSimHitMatcher::match(track, vertex);
 
   matchSimHitsToSimTrack();
 
   if (verbose_) {
     edm::LogInfo("ME0SimHitMatcher") << "nTrackIds " << track_ids_.size() << " nSelectedME0SimHits " << hits_.size();
     edm::LogInfo("ME0SimHitMatcher") << "detids ME0 " << detIds().size();
 
     const auto& me0_ch_ids = detIds();
     for (const auto& id : me0_ch_ids) {
       const auto& me0_simhits = MuonSimHitMatcher::hitsInChamber(id);
       const auto& me0_simhits_gp = simHitsMeanPosition(me0_simhits);
       edm::LogInfo("ME0SimHitMatcher") << "me0chid " << ME0DetId(id) << ": nHits " << me0_simhits.size() << " phi "
                                        << me0_simhits_gp.phi() << " nCh " << chamber_to_hits_[id].size();
       const auto& strips = hitStripsInDetId(id);
       edm::LogInfo("ME0SimHitMatcher") << "nStrip " << strips.size();
       edm::LogInfo("ME0SimHitMatcher") << "strips : ";
       for (const auto& p : strips) {
         edm::LogInfo("ME0SimHitMatcher") << p;
       }
     }
   }
 }

void ME0SimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 49 of file ME0SimHitMatcher.cc.

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, ztail::d, MuonSimHitMatcher::detid_to_hits_, detIds(), detids_to_pads_, MuonSimHitMatcher::discardEleHits_, MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), gpuClustering::id, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, superChamber_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 ME0 chambers
       if (discardEleHits_ && std::abs(pdgid) == 11)
         continue;
 
       const ME0DetId& layer_id(h.detUnitId());
       detid_to_hits_[h.detUnitId()].push_back(h);
       hits_.push_back(h);
       chamber_to_hits_[layer_id.layerId().rawId()].push_back(h);
       superChamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
     }
   }
 
   // find pads with hits
   const auto& detids = detIds();
   for (const auto& d : detids) {
     ME0DetId id(d);
     const auto& hits = hitsInDetId(d);
     const auto& roll = dynamic_cast<const ME0Geometry*>(geometry_)->etaPartition(id);
     // int max_npads = roll->npads();
     set<int> pads;
     for (const auto& h : hits) {
       const LocalPoint& lp = h.entryPoint();
       pads.insert(1 + static_cast<int>(roll->padTopology().channel(lp)));
     }
     detids_to_pads_[d] = pads;
   }
 }

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

Definition at line 133 of file ME0SimHitMatcher.cc.

References superChamberIdsCoincidences().

                                                                  {
   return superChamberIdsCoincidences(min_n_layers).size();
 }

int ME0SimHitMatcher::nLayersWithHitsInSuperChamber ( unsigned int detid ) const

Definition at line 112 of file ME0SimHitMatcher.cc.

References h, and hitsInSuperChamber().

Referenced by superChamberIdsCoincidences().

                                                                             {
   set<int> layers_with_hits;
   const auto& hits = hitsInSuperChamber(detid);
   for (const auto& h : hits) {
     const ME0DetId& idd(h.detUnitId());
     layers_with_hits.insert(idd.layer());
   }
   return layers_with_hits.size();
 }

int ME0SimHitMatcher::nPadsWithHits ( ) const

Definition at line 193 of file ME0SimHitMatcher.cc.

References detIds(), hitPadsInDetId(), and mps_fire::result.

                                           {
   int result = 0;
   const auto& pad_ids = detIds();
   for (const auto& id : pad_ids) {
     result += hitPadsInDetId(id).size();
   }
   return result;
 }

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

Definition at line 137 of file ME0SimHitMatcher.cc.

References ztail::d, validate-o2o-wbm::f, MuonSimHitMatcher::geometry_, h, 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 ME0Geometry*>(geometry_)->etaPartition(d)->strip(lp);
     sums += s;
     ++n;
   }
   if (n == 0)
     return -1.f;
   return sums / n;
 }

std::set< unsigned int > ME0SimHitMatcher::superChamberIds ( ) const

Definition at line 99 of file ME0SimHitMatcher.cc.

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

Referenced by superChamberIdsCoincidences().

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

std::set< unsigned int > ME0SimHitMatcher::superChamberIdsCoincidences ( int min_n_layers = 4 ) const

Definition at line 122 of file ME0SimHitMatcher.cc.

References nLayersWithHitsInSuperChamber(), AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and superChamberIds().

Referenced by nCoincidenceChambers().

                                                                                          {
   set<unsigned int> result;
   // loop over the super chamber Ids
   for (const auto& p : superChamberIds()) {
     if (nLayersWithHitsInSuperChamber(p) >= min_n_layers) {
       result.insert(p);
     }
   }
   return result;
 }