DTSimHitMatcher Class Reference

#include <DTSimHitMatcher.h>

Inheritance diagram for DTSimHitMatcher:

Public Member Functions
std::set< unsigned int >	chamberIds (int type=MuonHitHelper::DT_ALL) const
std::set< unsigned int >	chamberIdsStation (int station) const
std::set< unsigned int >	detIds (int type=MuonHitHelper::DT_ALL) const
void	dtChamberIdsToString (const std::set< unsigned int > &) const
	DTSimHitMatcher (const edm::ParameterSet &iPS, edm::ConsumesCollector &&iC)
const edm::PSimHitContainer &	hitsInChamber (unsigned int) const
const edm::PSimHitContainer &	hitsInLayer (unsigned int) const
const edm::PSimHitContainer &	hitsInSuperLayer (unsigned int) const
bool	hitStation (int, int, int) const
std::set< unsigned int >	hitWiresInDTChamberId (unsigned int, int margin_n_wires=0) const
std::set< unsigned int >	hitWiresInDTLayerId (unsigned int, int margin_n_wires=0) const
std::set< unsigned int >	hitWiresInDTSuperLayerId (unsigned int, int margin_n_wires=0) const
void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
	initialize the event More...
std::set< unsigned int >	layerIds () const
void	match (const SimTrack &t, const SimVertex &v)
	do the matching More...
int	nCellsWithHitsInLayer (unsigned int) const
int	nLayersWithHitsInChamber (unsigned int) const
int	nLayersWithHitsInSuperLayer (unsigned int) const
int	nStations (int nsl=1, int nl=3) const
int	nSuperLayersWithHitsInChamber (unsigned int) const
GlobalPoint	simHitsMeanPositionStation (int n) const
float	simHitsMeanWire (const edm::PSimHitContainer &sim_hits) const
std::set< unsigned int >	superlayerIds () const
	~DTSimHitMatcher ()
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< DTGeometry >	dt_geom_
std::map< unsigned int, edm::PSimHitContainer >	layer_to_hits_
std::map< unsigned int, edm::PSimHitContainer >	superlayer_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 DT SimHit to SimTrack

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

Definition at line 14 of file DTSimHitMatcher.h.

Constructor & Destructor Documentation

DTSimHitMatcher()

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

Definition at line 5 of file DTSimHitMatcher.cc.

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

~DTSimHitMatcher()

DTSimHitMatcher::~DTSimHitMatcher ( )

inline

Definition at line 20 of file DTSimHitMatcher.h.

{}

Member Function Documentation

chamberIds()

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

(

int

type = MuonHitHelper::DT_ALL

)

const

Definition at line 88 of file DTSimHitMatcher.cc.

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

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, DTChamberId::station(), MuonHitHelper::toDTType(), and DTChamberId::wheel().

Referenced by chamberIdsStation(), and hitStation().

chamberIdsStation()

std::set< unsigned int > DTSimHitMatcher::chamberIdsStation

(

int

station

)

const

Definition at line 296 of file DTSimHitMatcher.cc.

                                                                         {
  set<unsigned int> result;
  switch (station) {
    case 1: {
      const auto& p1(chamberIds(MuonHitHelper::DT_MB21p));
      const auto& p2(chamberIds(MuonHitHelper::DT_MB11p));
      const auto& p3(chamberIds(MuonHitHelper::DT_MB01));
      const auto& p4(chamberIds(MuonHitHelper::DT_MB11n));
      const auto& p5(chamberIds(MuonHitHelper::DT_MB21n));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      result.insert(p3.begin(), p3.end());
      result.insert(p4.begin(), p4.end());
      result.insert(p5.begin(), p5.end());
      break;
    }
    case 2: {
      const auto& p1(chamberIds(MuonHitHelper::DT_MB22p));
      const auto& p2(chamberIds(MuonHitHelper::DT_MB12p));
      const auto& p3(chamberIds(MuonHitHelper::DT_MB02));
      const auto& p4(chamberIds(MuonHitHelper::DT_MB12n));
      const auto& p5(chamberIds(MuonHitHelper::DT_MB22n));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      result.insert(p3.begin(), p3.end());
      result.insert(p4.begin(), p4.end());
      result.insert(p5.begin(), p5.end());
      break;
    }
    case 3: {
      const auto& p1(chamberIds(MuonHitHelper::DT_MB23p));
      const auto& p2(chamberIds(MuonHitHelper::DT_MB13p));
      const auto& p3(chamberIds(MuonHitHelper::DT_MB03));
      const auto& p4(chamberIds(MuonHitHelper::DT_MB13n));
      const auto& p5(chamberIds(MuonHitHelper::DT_MB23n));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      result.insert(p3.begin(), p3.end());
      result.insert(p4.begin(), p4.end());
      result.insert(p5.begin(), p5.end());
      break;
    }
    case 4: {
      const auto& p1(chamberIds(MuonHitHelper::DT_MB24p));
      const auto& p2(chamberIds(MuonHitHelper::DT_MB14p));
      const auto& p3(chamberIds(MuonHitHelper::DT_MB04));
      const auto& p4(chamberIds(MuonHitHelper::DT_MB14n));
      const auto& p5(chamberIds(MuonHitHelper::DT_MB24n));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      result.insert(p3.begin(), p3.end());
      result.insert(p4.begin(), p4.end());
      result.insert(p5.begin(), p5.end());
      break;
    }
  };
  return result;
}

References chamberIds(), MuonHitHelper::DT_MB01, MuonHitHelper::DT_MB02, MuonHitHelper::DT_MB03, MuonHitHelper::DT_MB04, MuonHitHelper::DT_MB11n, MuonHitHelper::DT_MB11p, MuonHitHelper::DT_MB12n, MuonHitHelper::DT_MB12p, MuonHitHelper::DT_MB13n, MuonHitHelper::DT_MB13p, MuonHitHelper::DT_MB14n, MuonHitHelper::DT_MB14p, MuonHitHelper::DT_MB21n, MuonHitHelper::DT_MB21p, MuonHitHelper::DT_MB22n, MuonHitHelper::DT_MB22p, MuonHitHelper::DT_MB23n, MuonHitHelper::DT_MB23p, MuonHitHelper::DT_MB24n, MuonHitHelper::DT_MB24p, p1, p2, p3, p4, mps_fire::result, and relativeConstraints::station.

detIds()

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

(

int

type = MuonHitHelper::DT_ALL

)

const

Definition at line 74 of file DTSimHitMatcher.cc.

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

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, DTChamberId::station(), MuonHitHelper::toDTType(), and DTChamberId::wheel().

Referenced by match().

dtChamberIdsToString()

void DTSimHitMatcher::dtChamberIdsToString

(

const std::set< unsigned int > &

set

)

const

Definition at line 289 of file DTSimHitMatcher.cc.

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

References AlCaHLTBitMon_ParallelJobs::p.

hitsInChamber()

const edm::PSimHitContainer & DTSimHitMatcher::hitsInChamber

(

unsigned int

detid

)

const

Definition at line 136 of file DTSimHitMatcher.cc.

                                                                                  {
  if (!MuonHitHelper::isDT(detid))
    return no_hits_;
 
  const DTWireId id(detid);
  if (chamber_to_hits_.find(id.chamberId().rawId()) == chamber_to_hits_.end())
    return no_hits_;
  return chamber_to_hits_.at(id.chamberId().rawId());
}

References MuonSimHitMatcher::chamber_to_hits_, triggerObjects_cff::id, MuonHitHelper::isDT(), and MuonSimHitMatcher::no_hits_.

hitsInLayer()

const edm::PSimHitContainer & DTSimHitMatcher::hitsInLayer

(

unsigned int

detid

)

const

Definition at line 116 of file DTSimHitMatcher.cc.

                                                                                {
  if (!MuonHitHelper::isDT(detid))
    return no_hits_;
 
  const DTWireId id(detid);
  if (layer_to_hits_.find(id.layerId().rawId()) == layer_to_hits_.end())
    return no_hits_;
  return layer_to_hits_.at(id.layerId().rawId());
}

References triggerObjects_cff::id, MuonHitHelper::isDT(), layer_to_hits_, and MuonSimHitMatcher::no_hits_.

Referenced by nCellsWithHitsInLayer().

hitsInSuperLayer()

const edm::PSimHitContainer & DTSimHitMatcher::hitsInSuperLayer

(

unsigned int

detid

)

const

Definition at line 126 of file DTSimHitMatcher.cc.

                                                                                     {
  if (!MuonHitHelper::isDT(detid))
    return no_hits_;
 
  const DTWireId id(detid);
  if (superlayer_to_hits_.find(id.superlayerId().rawId()) == superlayer_to_hits_.end())
    return no_hits_;
  return superlayer_to_hits_.at(id.superlayerId().rawId());
}

References triggerObjects_cff::id, MuonHitHelper::isDT(), MuonSimHitMatcher::no_hits_, and superlayer_to_hits_.

Referenced by nLayersWithHitsInSuperLayer().

hitStation()

bool DTSimHitMatcher::hitStation	(	int	st,
		int	nsuperlayers,
		int	nlayers
	)		const

Definition at line 146 of file DTSimHitMatcher.cc.

                                                                            {
  int nst = 0;
  for (const auto& ddt : chamberIds()) {
    const DTChamberId id(ddt);
    if (id.station() != st)
      continue;
 
    // require at least 1 superlayer
    const int nsl(nSuperLayersWithHitsInChamber(id.rawId()));
    if (nsl < nsuperlayers)
      continue;
 
    // require at least 3 layers hit per chamber
    const int nl(nLayersWithHitsInChamber(id.rawId()));
    if (nl < nlayers)
      continue;
    ++nst;
  }
  return nst;
}

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

Referenced by nStations().

hitWiresInDTChamberId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTChamberId	(	unsigned int	detid,
		int	margin_n_wires = 0
	)		const

Definition at line 277 of file DTSimHitMatcher.cc.

                                                                                                        {
  set<unsigned int> result;
  const auto& superLayers(dynamic_cast<const DTGeometry*>(geometry_)->chamber(DTChamberId(detid))->superLayers());
  for (const auto& sl : superLayers) {
    if (verbose_)
      edm::LogInfo("DTSimHitMatcher") << "hitWiresInDTChamberId::sl id " << sl->id() << endl;
    const auto& p(hitWiresInDTSuperLayerId(sl->id().rawId(), margin_n_wires));
    result.insert(p.begin(), p.end());
  }
  return result;
}

References relativeConstraints::chamber, MuonSimHitMatcher::geometry_, hitWiresInDTSuperLayerId(), AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonSimHitMatcher::verbose_.

hitWiresInDTLayerId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTLayerId	(	unsigned int	detid,
		int	margin_n_wires = 0
	)		const

Definition at line 239 of file DTSimHitMatcher.cc.

                                                                                                      {
  set<unsigned int> result;
  if (MuonHitHelper::isDT(detid)) {
    DTLayerId id(detid);
    int max_nwires = dynamic_cast<const DTGeometry*>(geometry_)->layer(id)->specificTopology().channels();
    for (int wn = 0; wn <= max_nwires; ++wn) {
      DTWireId wid(id, wn);
      for (const auto& h : MuonSimHitMatcher::hitsInDetId(wid.rawId())) {
        if (verbose_)
          edm::LogInfo("DTSimHitMatcher") << "central DTWireId " << wid << " simhit " << h << endl;
        int smin = wn - margin_n_wires;
        smin = (smin > 0) ? smin : 1;
        int smax = wn + margin_n_wires;
        smax = (smax <= max_nwires) ? smax : max_nwires;
        for (int ss = smin; ss <= smax; ++ss) {
          DTWireId widd(id, ss);
          if (verbose_)
            edm::LogInfo("DTSimHitMatcher") << "\tadding DTWireId to collection " << widd << endl;
          result.insert(widd.rawId());
        }
      }
    }
  }
  return result;
}

References MuonSimHitMatcher::geometry_, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, MuonHitHelper::isDT(), DetId::rawId(), mps_fire::result, contentValuesCheck::ss, and MuonSimHitMatcher::verbose_.

Referenced by hitWiresInDTSuperLayerId().

hitWiresInDTSuperLayerId()

std::set< unsigned int > DTSimHitMatcher::hitWiresInDTSuperLayerId	(	unsigned int	detid,
		int	margin_n_wires = 0
	)		const

Definition at line 265 of file DTSimHitMatcher.cc.

                                                                                                           {
  set<unsigned int> result;
  const auto& layers(dynamic_cast<const DTGeometry*>(geometry_)->superLayer(DTSuperLayerId(detid))->layers());
  for (const auto& l : layers) {
    if (verbose_)
      edm::LogInfo("DTSimHitMatcher") << "hitWiresInDTSuperLayerId::l id " << l->id() << endl;
    const auto& p(hitWiresInDTLayerId(l->id().rawId(), margin_n_wires));
    result.insert(p.begin(), p.end());
  }
  return result;
}

References MuonSimHitMatcher::geometry_, hitWiresInDTLayerId(), cmsLHEtoEOSManager::l, hgcalTopologyTester_cfi::layers, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, and MuonSimHitMatcher::verbose_.

Referenced by hitWiresInDTChamberId().

init()

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

initialize the event

Definition at line 16 of file DTSimHitMatcher.cc.

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

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

layerIds()

std::set< unsigned int > DTSimHitMatcher::layerIds

(

)

const

Definition at line 102 of file DTSimHitMatcher.cc.

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

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

match()

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

do the matching

Definition at line 28 of file DTSimHitMatcher.cc.

                                                                          {
  // instantiates the track ids and simhits
  MuonSimHitMatcher::match(track, vertex);
 
  if (hasGeometry_) {
    matchSimHitsToSimTrack();
 
    if (verbose_) {
      edm::LogInfo("DTSimHitMatcher") << "nTrackIds " << track_ids_.size() << " nSelectedDTSimHits " << hits_.size()
                                      << endl;
      edm::LogInfo("DTSimHitMatcher") << "detids DT " << detIds(0).size() << endl;
 
      const auto& dt_det_ids = detIds(0);
      for (const auto& id : dt_det_ids) {
        const auto& dt_simhits = MuonSimHitMatcher::hitsInDetId(id);
        const auto& dt_simhits_gp = simHitsMeanPosition(dt_simhits);
        edm::LogInfo("DTSimHitMatcher") << "DTWireId " << DTWireId(id) << ": nHits " << dt_simhits.size() << " eta "
                                        << dt_simhits_gp.eta() << " phi " << dt_simhits_gp.phi() << " nCh "
                                        << chamber_to_hits_[id].size() << endl;
      }
    }
  }
}

References MuonSimHitMatcher::chamber_to_hits_, detIds(), MuonSimHitMatcher::hasGeometry_, MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), MuonSimHitMatcher::simHitsMeanPosition(), HLT_2018_cff::track, MuonSimHitMatcher::track_ids_, MuonSimHitMatcher::verbose_, and bphysicsOniaDQM_cfi::vertex.

matchSimHitsToSimTrack()

void DTSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 52 of file DTSimHitMatcher.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 DT chambers
      if (discardEleHits_ && pdgid == 11)
        continue;
 
      const DTWireId layer_id(h.detUnitId());
      detid_to_hits_[h.detUnitId()].push_back(h);
      hits_.push_back(h);
      layer_to_hits_[layer_id.layerId().rawId()].push_back(h);
      superlayer_to_hits_[layer_id.superlayerId().rawId()].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_, layer_to_hits_, EgammaValidation_cff::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, superlayer_to_hits_, and MuonSimHitMatcher::track_ids_.

Referenced by match().

nCellsWithHitsInLayer()

int DTSimHitMatcher::nCellsWithHitsInLayer

(

unsigned int

detid

)

const

Definition at line 172 of file DTSimHitMatcher.cc.

                                                                   {
  set<int> layers_with_hits;
  const auto& hits = hitsInLayer(detid);
  for (const auto& h : hits) {
    if (MuonHitHelper::isDT(detid)) {
      const DTWireId idd(h.detUnitId());
      layers_with_hits.insert(idd.wire());
    }
  }
  return layers_with_hits.size();
}

References hfClusterShapes_cfi::hits, hitsInLayer(), and MuonHitHelper::isDT().

nLayersWithHitsInChamber()

int DTSimHitMatcher::nLayersWithHitsInChamber

(

unsigned int

detid

)

const

Definition at line 208 of file DTSimHitMatcher.cc.

                                                                      {
  int nLayers = 0;
  const auto& superLayers(dynamic_cast<const DTGeometry*>(geometry_)->chamber(DTChamberId(detid))->superLayers());
  for (const auto& sl : superLayers) {
    nLayers += nLayersWithHitsInSuperLayer(sl->id().rawId());
  }
  return nLayers;
}

References relativeConstraints::chamber, MuonSimHitMatcher::geometry_, MuonTCMETValueMapProducer_cff::nLayers, and nLayersWithHitsInSuperLayer().

Referenced by hitStation().

nLayersWithHitsInSuperLayer()

int DTSimHitMatcher::nLayersWithHitsInSuperLayer

(

unsigned int

detid

)

const

Definition at line 184 of file DTSimHitMatcher.cc.

                                                                         {
  set<int> layers_with_hits;
  const auto& hits = hitsInSuperLayer(detid);
  for (const auto& h : hits) {
    if (MuonHitHelper::isDT(detid)) {
      const DTLayerId idd(h.detUnitId());
      layers_with_hits.insert(idd.layer());
    }
  }
  return layers_with_hits.size();
}

References hfClusterShapes_cfi::hits, hitsInSuperLayer(), and MuonHitHelper::isDT().

Referenced by nLayersWithHitsInChamber().

nStations()

int DTSimHitMatcher::nStations	(	int	nsl = 1,
		int	nl = 3
	)		const

Definition at line 167 of file DTSimHitMatcher.cc.

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

References hitStation(), and nlayers.

nSuperLayersWithHitsInChamber()

int DTSimHitMatcher::nSuperLayersWithHitsInChamber

(

unsigned int

detid

)

const

Definition at line 196 of file DTSimHitMatcher.cc.

                                                                           {
  set<int> sl_with_hits;
  const auto& hits = MuonSimHitMatcher::hitsInChamber(detid);
  for (const auto& h : hits) {
    if (MuonHitHelper::isDT(detid)) {
      const DTSuperLayerId idd(h.detUnitId());
      sl_with_hits.insert(idd.superLayer());
    }
  }
  return sl_with_hits.size();
}

References hfClusterShapes_cfi::hits, MuonSimHitMatcher::hitsInChamber(), and MuonHitHelper::isDT().

Referenced by hitStation().

simHitsMeanPositionStation()

GlobalPoint DTSimHitMatcher::simHitsMeanPositionStation

(

int

n

)

const

simHitsMeanWire()

float DTSimHitMatcher::simHitsMeanWire

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 216 of file DTSimHitMatcher.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();
    if (MuonHitHelper::isDT(d)) {
      // find nearest wire
      s = dynamic_cast<const DTGeometry*>(geometry_)->layer(DTLayerId(d))->specificTopology().channel(lp);
    } else
      continue;
    sums += s;
    ++n;
  }
  if (n == 0)
    return -1.f;
  return sums / n;
}

References ztail::d, f, MuonSimHitMatcher::geometry_, MuonHitHelper::isDT(), dqmiodumpmetadata::n, and alignCSCRings::s.

superlayerIds()

std::set< unsigned int > DTSimHitMatcher::superlayerIds

(

)

const

Definition at line 109 of file DTSimHitMatcher.cc.

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

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

Member Data Documentation

dt_geom_

edm::ESHandle<DTGeometry> DTSimHitMatcher::dt_geom_

private

Definition at line 76 of file DTSimHitMatcher.h.

Referenced by init().

layer_to_hits_

std::map<unsigned int, edm::PSimHitContainer> DTSimHitMatcher::layer_to_hits_

private

Definition at line 73 of file DTSimHitMatcher.h.

Referenced by hitsInLayer(), layerIds(), and matchSimHitsToSimTrack().

superlayer_to_hits_

std::map<unsigned int, edm::PSimHitContainer> DTSimHitMatcher::superlayer_to_hits_

private

Definition at line 74 of file DTSimHitMatcher.h.

Referenced by hitsInSuperLayer(), matchSimHitsToSimTrack(), and superlayerIds().