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	(	const edm::ParameterSet &	iPS,
		edm::ConsumesCollector &&	iC
	)

Definition at line 5 of file DTSimHitMatcher.cc.

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

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

DTSimHitMatcher::~DTSimHitMatcher ( )

inline

Definition at line 20 of file DTSimHitMatcher.h.

References chamberIds(), chamberIdsStation(), detIds(), MuonHitHelper::DT_ALL, dtChamberIdsToString(), MillePedeFileConverter_cfg::e, hitsInChamber(), hitsInLayer(), hitsInSuperLayer(), hitStation(), hitWiresInDTChamberId(), hitWiresInDTLayerId(), hitWiresInDTSuperLayerId(), init(), layerIds(), match(), matchSimHitsToSimTrack(), gen::n, nCellsWithHitsInLayer(), nLayersWithHitsInChamber(), nLayersWithHitsInSuperLayer(), nStations(), nSuperLayersWithHitsInChamber(), simHitsMeanPositionStation(), simHitsMeanWire(), relativeConstraints::station, superlayerIds(), protons_cff::t, and findQualityFiles::v.

{}

Member Function Documentation

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

(

int

type = MuonHitHelper::DT_ALL

)

const

Definition at line 92 of file DTSimHitMatcher.cc.

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

Referenced by chamberIdsStation(), hitStation(), and ~DTSimHitMatcher().

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

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

(

int

station

)

const

Definition at line 317 of file DTSimHitMatcher.cc.

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, and mps_fire::result.

Referenced by ~DTSimHitMatcher().

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

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

(

int

type = MuonHitHelper::DT_ALL

)

const

Definition at line 78 of file DTSimHitMatcher.cc.

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

Referenced by match(), and ~DTSimHitMatcher().

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

void DTSimHitMatcher::dtChamberIdsToString

(

const std::set< unsigned int > &

set

)

const

Definition at line 309 of file DTSimHitMatcher.cc.

References AlCaHLTBitMon_ParallelJobs::p.

Referenced by ~DTSimHitMatcher().

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

const edm::PSimHitContainer & DTSimHitMatcher::hitsInChamber

(

unsigned int

detid

)

const

Definition at line 139 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

const edm::PSimHitContainer & DTSimHitMatcher::hitsInLayer

(

unsigned int

detid

)

const

Definition at line 118 of file DTSimHitMatcher.cc.

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

Referenced by nCellsWithHitsInLayer(), and ~DTSimHitMatcher().

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

const edm::PSimHitContainer & DTSimHitMatcher::hitsInSuperLayer

(

unsigned int

detid

)

const

Definition at line 128 of file DTSimHitMatcher.cc.

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

Referenced by nLayersWithHitsInSuperLayer(), and ~DTSimHitMatcher().

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

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

Definition at line 149 of file DTSimHitMatcher.cc.

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

Referenced by nStations(), and ~DTSimHitMatcher().

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

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

Definition at line 293 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

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

Definition at line 245 of file DTSimHitMatcher.cc.

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

Referenced by hitWiresInDTSuperLayerId(), and ~DTSimHitMatcher().

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

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

Definition at line 277 of file DTSimHitMatcher.cc.

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

Referenced by hitWiresInDTChamberId(), and ~DTSimHitMatcher().

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

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

initialize the event

Definition at line 18 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

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

(

)

const

Definition at line 106 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

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

do the matching

Definition at line 32 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

void DTSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 59 of file DTSimHitMatcher.cc.

References funct::abs(), MuonSimHitMatcher::chamber_to_hits_, MuonSimHitMatcher::detid_to_hits_, MuonSimHitMatcher::discardEleHits_, h, MuonSimHitMatcher::hits_, layer_to_hits_, BPhysicsValidation_cfi::pdgid, MuonSimHitMatcher::simHits_, MuonSimHitMatcher::simMuOnly_, superlayer_to_hits_, and MuonSimHitMatcher::track_ids_.

Referenced by match(), and ~DTSimHitMatcher().

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

int DTSimHitMatcher::nCellsWithHitsInLayer

(

unsigned int

detid

)

const

Definition at line 174 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

int DTSimHitMatcher::nLayersWithHitsInChamber

(

unsigned int

detid

)

const

Definition at line 210 of file DTSimHitMatcher.cc.

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

Referenced by hitStation(), and ~DTSimHitMatcher().

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

int DTSimHitMatcher::nLayersWithHitsInSuperLayer

(

unsigned int

detid

)

const

Definition at line 186 of file DTSimHitMatcher.cc.

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

Referenced by nLayersWithHitsInChamber(), and ~DTSimHitMatcher().

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

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

Definition at line 167 of file DTSimHitMatcher.cc.

References hitStation().

Referenced by ~DTSimHitMatcher().

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

int DTSimHitMatcher::nSuperLayersWithHitsInChamber

(

unsigned int

detid

)

const

Definition at line 198 of file DTSimHitMatcher.cc.

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

Referenced by hitStation(), and ~DTSimHitMatcher().

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

GlobalPoint DTSimHitMatcher::simHitsMeanPositionStation

(

int

n

)

const

Referenced by ~DTSimHitMatcher().

float DTSimHitMatcher::simHitsMeanWire

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 220 of file DTSimHitMatcher.cc.

References edmIntegrityCheck::d, f, MuonSimHitMatcher::geometry_, h, MuonHitHelper::isDT(), gen::n, and alignCSCRings::s.

Referenced by ~DTSimHitMatcher().

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

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

(

)

const

Definition at line 112 of file DTSimHitMatcher.cc.

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

Referenced by ~DTSimHitMatcher().

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

Member Data Documentation

edm::ESHandle<DTGeometry> DTSimHitMatcher::dt_geom_

private

Definition at line 79 of file DTSimHitMatcher.h.

Referenced by init().

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

private

Definition at line 76 of file DTSimHitMatcher.h.

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

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

private

Definition at line 77 of file DTSimHitMatcher.h.

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