RPCSimHitMatcher Class Reference

#include <RPCSimHitMatcher.h>

Inheritance diagram for RPCSimHitMatcher:

Public Member Functions
std::set< unsigned int >	chamberIds (int type=MuonHitHelper::RPC_ALL) const
std::set< unsigned int >	detIds (int type=MuonHitHelper::RPC_ALL) const
bool	hitStation (int st) 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...
int	nStations () const
	RPCSimHitMatcher (const edm::ParameterSet &iPS, edm::ConsumesCollector &&iC)
float	simHitsMeanStrip (const edm::PSimHitContainer &sim_hits) const
	~RPCSimHitMatcher ()
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
edm::ESHandle< RPCGeometry >	rpc_geom_

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

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

Definition at line 14 of file RPCSimHitMatcher.h.

Constructor & Destructor Documentation

RPCSimHitMatcher()

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

Definition at line 6 of file RPCSimHitMatcher.cc.

    : MuonSimHitMatcher(ps, std::move(iC)) {
  simHitPSet_ = ps.getParameterSet("rpcSimHit");
  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_.

~RPCSimHitMatcher()

RPCSimHitMatcher::~RPCSimHitMatcher ( )

inline

Definition at line 20 of file RPCSimHitMatcher.h.

{}

Member Function Documentation

chamberIds()

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

(

int

type = MuonHitHelper::RPC_ALL

)

const

Definition at line 92 of file RPCSimHitMatcher.cc.

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

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, RPCDetId::region(), mps_fire::result, RPCDetId::ring(), RPCDetId::station(), and MuonHitHelper::toRPCType().

Referenced by hitStation(), and match().

detIds()

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

(

int

type = MuonHitHelper::RPC_ALL

)

const

Definition at line 78 of file RPCSimHitMatcher.cc.

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

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, RPCDetId::region(), mps_fire::result, RPCDetId::ring(), RPCDetId::station(), and MuonHitHelper::toRPCType().

Referenced by match().

hitStation()

bool RPCSimHitMatcher::hitStation

(

int

st

)

const

Definition at line 106 of file RPCSimHitMatcher.cc.

                                              {
  int nst = 0;
  for (const auto& ddt : chamberIds(0)) {
    const RPCDetId id(ddt);
    if (id.station() != st)
      continue;
    ++nst;
  }
  return nst;
}

References chamberIds(), triggerObjects_cff::id, and relativeConstraints::station.

Referenced by nStations().

hitStripsInDetId()

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

Definition at line 136 of file RPCSimHitMatcher.cc.

                                                                                            {
  set<int> result;
  RPCDetId id(detid);
  for (const auto& roll : dynamic_cast<const RPCGeometry*>(geometry_)->chamber(id)->rolls()) {
    int max_nstrips = roll->nstrips();
    for (const auto& h : MuonSimHitMatcher::hitsInDetId(roll->id().rawId())) {
      const LocalPoint& lp = h.entryPoint();
      int central_strip = static_cast<int>(roll->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 relativeConstraints::chamber, MuonSimHitMatcher::geometry_, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, mps_fire::result, and contentValuesCheck::ss.

Referenced by match().

init()

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

initialize the event

Definition at line 17 of file RPCSimHitMatcher.cc.

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

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

match()

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

do the matching

Definition at line 29 of file RPCSimHitMatcher.cc.

                                                                           {
  // instantiates the track ids and simhits
  MuonSimHitMatcher::match(track, vertex);
 
  if (hasGeometry_) {
    matchSimHitsToSimTrack();
 
    if (verbose_) {
      edm::LogInfo("RPCSimHitMatcher") << "nSimHits " << simHits_.size() << " nTrackIds " << track_ids_.size() << endl;
      edm::LogInfo("RPCSimHitMatcher") << "detids RPC " << detIds().size() << endl;
 
      const auto& ch_ids = chamberIds();
      for (const auto& id : ch_ids) {
        const auto& simhits = MuonSimHitMatcher::hitsInChamber(id);
        const auto& simhits_gp = simHitsMeanPosition(simhits);
        edm::LogInfo("RPCSimHitMatcher") << "RPCDetId " << RPCDetId(id) << ": nHits " << simhits.size() << " eta "
                                         << simhits_gp.eta() << " phi " << simhits_gp.phi() << " nCh "
                                         << chamber_to_hits_[id].size() << endl;
        const auto& strips = hitStripsInDetId(id);
        edm::LogInfo("RPCSimHitMatcher") << "nStrips " << strips.size() << endl;
        edm::LogInfo("RPCSimHitMatcher") << "strips : ";
        for (const auto& p : strips) {
          edm::LogInfo("RPCSimHitMatcher") << p;
        }
      }
    }
  }
}

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

matchSimHitsToSimTrack()

void RPCSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 58 of file RPCSimHitMatcher.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 RPC chambers
      if (discardEleHits_ && pdgid == 11)
        continue;
 
      const RPCDetId& layer_id(h.detUnitId());
      detid_to_hits_[h.detUnitId()].push_back(h);
      hits_.push_back(h);
      chamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
    }
  }
}

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, MuonSimHitMatcher::detid_to_hits_, MuonSimHitMatcher::discardEleHits_, MuonSimHitMatcher::hits_, EgammaValidation_cff::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, and MuonSimHitMatcher::track_ids_.

Referenced by match().

nStations()

int RPCSimHitMatcher::nStations

(

)

const

Definition at line 117 of file RPCSimHitMatcher.cc.

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

References hitStation().

simHitsMeanStrip()

float RPCSimHitMatcher::simHitsMeanStrip

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 119 of file RPCSimHitMatcher.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();
    const auto& d = h.detUnitId();
    sums += dynamic_cast<const RPCGeometry*>(geometry_)->roll(d)->strip(lp);
    ++n;
  }
  if (n == 0)
    return -1.f;
  return sums / n;
}

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

Member Data Documentation

rpc_geom_

edm::ESHandle<RPCGeometry> RPCSimHitMatcher::rpc_geom_

private

Definition at line 47 of file RPCSimHitMatcher.h.

Referenced by init().