ME0SimHitMatcher Class Reference

#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 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::ESHandle< ME0Geometry >	me0_geom_
std::map< unsigned int, edm::PSimHitContainer >	superChamber_to_hits_

Additional Inherited Members
Protected Member Functions inherited from MuonSimHitMatcher
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 14 of file ME0SimHitMatcher.h.

Constructor & Destructor Documentation

ME0SimHitMatcher()

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

Definition at line 5 of file ME0SimHitMatcher.cc.

    : 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"));
}

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

~ME0SimHitMatcher()

ME0SimHitMatcher::~ME0SimHitMatcher ( )

inline

Definition at line 20 of file ME0SimHitMatcher.h.

{}

Member Function Documentation

chamberIds()

std::set< unsigned int > ME0SimHitMatcher::chamberIds

(

)

const

Definition at line 100 of file ME0SimHitMatcher.cc.

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

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

detIds()

std::set< unsigned int > ME0SimHitMatcher::detIds

(

void

)

const

Definition at line 93 of file ME0SimHitMatcher.cc.

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

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

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

hitPadsInDetId()

std::set< int > ME0SimHitMatcher::hitPadsInDetId

(

unsigned int

detid

)

const

Definition at line 183 of file ME0SimHitMatcher.cc.

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

References detids_to_pads_, and none.

Referenced by nPadsWithHits().

hitPartitions()

std::set< int > ME0SimHitMatcher::hitPartitions

(

)

const

Definition at line 190 of file ME0SimHitMatcher.cc.

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

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

hitsInSuperChamber()

const edm::PSimHitContainer & ME0SimHitMatcher::hitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 114 of file ME0SimHitMatcher.cc.

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

References MuonSimHitMatcher::no_hits_, and superChamber_to_hits_.

Referenced by nLayersWithHitsInSuperChamber().

hitStripsInDetId()

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

Definition at line 164 of file ME0SimHitMatcher.cc.

                                                                                            {
  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;
}

References MuonSimHitMatcher::geometry_, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, mps_fire::result, HLTBitAnalyser_cfi::simhits, contentValuesCheck::ss, and ecaldqm::topology().

Referenced by match().

init()

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

initialize the event

Definition at line 16 of file ME0SimHitMatcher.cc.

                                                                               {
  iSetup.get<MuonGeometryRecord>().get(me0_geom_);
  if (me0_geom_.isValid()) {
    geometry_ = &*me0_geom_;
  } else {
    hasGeometry_ = false;
    edm::LogWarning("ME0SimHitMatcher") << "+++ Info: ME0 geometry is unavailable. +++\n";
  }
  MuonSimHitMatcher::init(iEvent, iSetup);
}

References MuonSimHitMatcher::geometry_, edm::EventSetup::get(), get, MuonSimHitMatcher::hasGeometry_, iEvent, MuonSimHitMatcher::init(), edm::ESHandleBase::isValid(), and me0_geom_.

match()

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

do the matching

Definition at line 28 of file ME0SimHitMatcher.cc.

                                                                           {
  // instantiates the track ids and simhits
  MuonSimHitMatcher::match(track, vertex);
 
  if (hasGeometry_) {
    matchSimHitsToSimTrack();
 
    if (verbose_) {
      edm::LogInfo("ME0SimHitMatcher") << "nTrackIds " << track_ids_.size() << " nSelectedME0SimHits " << hits_.size()
                                       << endl;
      edm::LogInfo("ME0SimHitMatcher") << "detids ME0 " << detIds().size() << endl;
 
      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() << endl;
        const auto& strips = hitStripsInDetId(id);
        edm::LogInfo("ME0SimHitMatcher") << "nStrip " << strips.size() << endl;
        edm::LogInfo("ME0SimHitMatcher") << "strips : ";
        for (const auto& p : strips) {
          edm::LogInfo("ME0SimHitMatcher") << p;
        }
      }
    }
  }
}

References MuonSimHitMatcher::chamber_to_hits_, detIds(), MuonSimHitMatcher::hasGeometry_, MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInChamber(), hitStripsInDetId(), triggerObjects_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), AlCaHLTBitMon_ParallelJobs::p, MuonSimHitMatcher::simHitsMeanPosition(), DigiDM_cff::strips, HLT_2018_cff::track, MuonSimHitMatcher::track_ids_, MuonSimHitMatcher::verbose_, and bphysicsOniaDQM_cfi::vertex.

matchSimHitsToSimTrack()

void ME0SimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 57 of file ME0SimHitMatcher.cc.

                                              {
  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;
  }
}

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, ztail::d, MuonSimHitMatcher::detid_to_hits_, detIds(), detids_to_pads_, MuonSimHitMatcher::discardEleHits_, MuonSimHitMatcher::geometry_, hfClusterShapes_cfi::hits, MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, EgammaValidation_cff::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, superChamber_to_hits_, and MuonSimHitMatcher::track_ids_.

Referenced by match().

nCoincidenceChambers()

int ME0SimHitMatcher::nCoincidenceChambers

(

int

min_n_layers = 4

)

const

Definition at line 141 of file ME0SimHitMatcher.cc.

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

References superChamberIdsCoincidences().

nLayersWithHitsInSuperChamber()

int ME0SimHitMatcher::nLayersWithHitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 120 of file ME0SimHitMatcher.cc.

                                                                            {
  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();
}

References hfClusterShapes_cfi::hits, and hitsInSuperChamber().

Referenced by superChamberIdsCoincidences().

nPadsWithHits()

int ME0SimHitMatcher::nPadsWithHits

(

)

const

Definition at line 201 of file ME0SimHitMatcher.cc.

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

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

simHitsMeanStrip()

float ME0SimHitMatcher::simHitsMeanStrip

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 145 of file ME0SimHitMatcher.cc.

                                                                                  {
  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;
}

References ztail::d, f, MuonSimHitMatcher::geometry_, dqmiodumpmetadata::n, and alignCSCRings::s.

superChamberIds()

std::set< unsigned int > ME0SimHitMatcher::superChamberIds

(

)

const

Definition at line 107 of file ME0SimHitMatcher.cc.

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

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

Referenced by superChamberIdsCoincidences().

superChamberIdsCoincidences()

std::set< unsigned int > ME0SimHitMatcher::superChamberIdsCoincidences

(

int

min_n_layers = 4

)

const

Definition at line 130 of file ME0SimHitMatcher.cc.

                                                                                         {
  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;
}

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

Referenced by nCoincidenceChambers().

Member Data Documentation

detids_to_copads_

std::map<unsigned int, std::set<int> > ME0SimHitMatcher::detids_to_copads_

private

Definition at line 71 of file ME0SimHitMatcher.h.

detids_to_pads_

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

private

Definition at line 68 of file ME0SimHitMatcher.h.

Referenced by hitPadsInDetId(), and matchSimHitsToSimTrack().

me0_geom_

edm::ESHandle<ME0Geometry> ME0SimHitMatcher::me0_geom_

private

Definition at line 65 of file ME0SimHitMatcher.h.

Referenced by init().

superChamber_to_hits_

std::map<unsigned int, edm::PSimHitContainer> ME0SimHitMatcher::superChamber_to_hits_

private

Definition at line 73 of file ME0SimHitMatcher.h.

Referenced by hitsInSuperChamber(), matchSimHitsToSimTrack(), and superChamberIds().