CSCSimHitMatcher Class Reference

#include <CSCSimHitMatcher.h>

Inheritance diagram for CSCSimHitMatcher:

Public Member Functions
void	camberIdsToString (const std::set< unsigned int > &) const
std::set< unsigned int >	chamberIds (int type=MuonHitHelper::CSC_ALL) const
std::set< unsigned int >	chamberIdsStation (int station) const
void	chamberIdsToString (const std::set< unsigned int > &set) const
	CSCSimHitMatcher (const edm::ParameterSet &iPS, edm::ConsumesCollector &&iC)
std::set< unsigned int >	detIds (int type=MuonHitHelper::CSC_ALL) const
bool	hitStation (int, int) const
std::set< int >	hitStripsInDetId (unsigned int, int margin_n_strips=0) const
std::set< int >	hitWiregroupsInDetId (unsigned int, int margin_n_wg=0) const
void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
	initialize the event More...
float	LocalBendingInChamber (unsigned int detid) const
void	match (const SimTrack &t, const SimVertex &v)
	do the matching More...
int	nCoincidenceChambers (int min_n_layers=4) const
int	nLayersWithHitsInChamber (unsigned int) const
int	nStations (int nl=4) const
GlobalPoint	simHitPositionKeyLayer (unsigned int chamberid) const
GlobalPoint	simHitsMeanPositionStation (int n) const
float	simHitsMeanStrip (const edm::PSimHitContainer &sim_hits) const
float	simHitsMeanWG (const edm::PSimHitContainer &sim_hits) const
	~CSCSimHitMatcher ()
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< CSCGeometry >	csc_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 CSC SimHit to SimTrack

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

Definition at line 14 of file CSCSimHitMatcher.h.

Constructor & Destructor Documentation

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

Definition at line 5 of file CSCSimHitMatcher.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("cscSimHit");
  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"));
}

CSCSimHitMatcher::~CSCSimHitMatcher ( )

inline

Definition at line 20 of file CSCSimHitMatcher.h.

References camberIdsToString(), chamberIds(), chamberIdsStation(), chamberIdsToString(), MuonHitHelper::CSC_ALL, detIds(), MillePedeFileConverter_cfg::e, hitStation(), hitStripsInDetId(), hitWiregroupsInDetId(), init(), LocalBendingInChamber(), match(), matchSimHitsToSimTrack(), gen::n, nCoincidenceChambers(), nLayersWithHitsInChamber(), nStations(), simHitPositionKeyLayer(), simHitsMeanPositionStation(), simHitsMeanStrip(), simHitsMeanWG(), relativeConstraints::station, lumiQTWidget::t, and findQualityFiles::v.

{}

Member Function Documentation

void CSCSimHitMatcher::camberIdsToString

(

const std::set< unsigned int > &

)

const

Referenced by ~CSCSimHitMatcher().

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

(

int

type = MuonHitHelper::CSC_ALL

)

const

Definition at line 98 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::chamber_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, CSCDetId::ring(), CSCDetId::station(), and MuonHitHelper::toCSCType().

Referenced by chamberIdsStation(), hitStation(), nCoincidenceChambers(), and ~CSCSimHitMatcher().

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

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

(

int

station

)

const

Definition at line 334 of file CSCSimHitMatcher.cc.

References chamberIds(), MuonHitHelper::CSC_ME12, MuonHitHelper::CSC_ME13, MuonHitHelper::CSC_ME1a, MuonHitHelper::CSC_ME1b, MuonHitHelper::CSC_ME21, MuonHitHelper::CSC_ME22, MuonHitHelper::CSC_ME31, MuonHitHelper::CSC_ME32, MuonHitHelper::CSC_ME41, MuonHitHelper::CSC_ME42, p1, p2, p3, p4, and mps_fire::result.

Referenced by ~CSCSimHitMatcher().

                                                                          {
  set<unsigned int> result;
  switch (station) {
    case 1: {
      const auto& p1(chamberIds(MuonHitHelper::CSC_ME1a));
      const auto& p2(chamberIds(MuonHitHelper::CSC_ME1b));
      const auto& p3(chamberIds(MuonHitHelper::CSC_ME12));
      const auto& p4(chamberIds(MuonHitHelper::CSC_ME13));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      result.insert(p3.begin(), p3.end());
      result.insert(p4.begin(), p4.end());
      break;
    }
    case 2: {
      const auto& p1(chamberIds(MuonHitHelper::CSC_ME21));
      const auto& p2(chamberIds(MuonHitHelper::CSC_ME22));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      break;
    }
    case 3: {
      const auto& p1(chamberIds(MuonHitHelper::CSC_ME31));
      const auto& p2(chamberIds(MuonHitHelper::CSC_ME32));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      break;
    }
    case 4: {
      const auto& p1(chamberIds(MuonHitHelper::CSC_ME41));
      const auto& p2(chamberIds(MuonHitHelper::CSC_ME42));
      result.insert(p1.begin(), p1.end());
      result.insert(p2.begin(), p2.end());
      break;
    }
  };
  return result;
}

void CSCSimHitMatcher::chamberIdsToString

(

const std::set< unsigned int > &

set

)

const

Definition at line 326 of file CSCSimHitMatcher.cc.

References AlCaHLTBitMon_ParallelJobs::p.

Referenced by ~CSCSimHitMatcher().

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

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

(

int

type = MuonHitHelper::CSC_ALL

)

const

Definition at line 84 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::detid_to_hits_, AlCaHLTBitMon_ParallelJobs::p, mps_fire::result, CSCDetId::ring(), CSCDetId::station(), and MuonHitHelper::toCSCType().

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

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

bool CSCSimHitMatcher::hitStation	(	int	st,
		int	nlayers
	)		const

Definition at line 121 of file CSCSimHitMatcher.cc.

References relativeConstraints::chamber, chamberIds(), makeMuonMisalignmentScenario::endcap, triggerObjects_cff::id, nLayersWithHitsInChamber(), DetId::rawId(), relativeConstraints::ring, and relativeConstraints::station.

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

                                                           {
  int nst = 0;
  for (const auto& ddt : chamberIds()) {
    const CSCDetId id(ddt);
    int ri(id.ring());
    if (id.station() != st) continue;

    // ME1/a case - check the ME1/b chamber
    if (st == 1 and ri == 4) {
      CSCDetId idME1b(id.endcap(), id.station(), 1, id.chamber(), 0);
      const int nl1a(nLayersWithHitsInChamber(id.rawId()));
      const int nl1b(nLayersWithHitsInChamber(idME1b.rawId()));
      if (nl1a + nl1b < nlayers) continue;
      ++nst;
    }
    // ME1/b case - check the ME1/a chamber
    else if (st == 1 and ri == 1) {
      CSCDetId idME1a(id.endcap(), id.station(), 4, id.chamber(), 0);
      const int nl1a(nLayersWithHitsInChamber(idME1a.rawId()));
      const int nl1b(nLayersWithHitsInChamber(id.rawId()));
      if (nl1a + nl1b < nlayers) continue;
      ++nst;
    }
    // default case
    else {
      const int nl(nLayersWithHitsInChamber(id.rawId()));
      if (nl < nlayers) continue;
      ++nst;
    }
  }
  return nst;
}

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

Definition at line 270 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, and mps_fire::result.

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

                                                                            {
  set<int> result;
  const auto& simhits = MuonSimHitMatcher::hitsInDetId(detid);
  CSCDetId id(detid);
  int max_nstrips = dynamic_cast<const CSCGeometry*>(geometry_)
                        ->layer(id)
                        ->geometry()
                        ->numberOfStrips();
  for (const auto& h : simhits) {
    const LocalPoint& lp = h.entryPoint();
    int central_strip = dynamic_cast<const CSCGeometry*>(geometry_)
                            ->layer(id)
                            ->geometry()
                            ->nearestStrip(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;
}

std::set< int > CSCSimHitMatcher::hitWiregroupsInDetId	(	unsigned int	detid,
		int	margin_n_wg = 0
	)		const

Definition at line 294 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::geometry_, h, MuonSimHitMatcher::hitsInDetId(), triggerObjects_cff::id, and mps_fire::result.

Referenced by ~CSCSimHitMatcher().

                                                                            {
  set<int> result;
  const auto& simhits = MuonSimHitMatcher::hitsInDetId(detid);
  CSCDetId id(detid);
  int max_n_wg = dynamic_cast<const CSCGeometry*>(geometry_)
                     ->layer(id)
                     ->geometry()
                     ->numberOfWireGroups();
  for (const auto& h : simhits) {
    const LocalPoint& lp = h.entryPoint();
    const auto& layer_geo =
        dynamic_cast<const CSCGeometry*>(geometry_)->layer(id)->geometry();
    int central_wg = layer_geo->wireGroup(layer_geo->nearestWire(lp));
    int wg_min = central_wg - margin_n_wg;
    wg_min = (wg_min > 0) ? wg_min : 1;
    int wg_max = central_wg + margin_n_wg;
    wg_max = (wg_max <= max_n_wg) ? wg_max : max_n_wg;
    for (int wg = wg_min; wg <= wg_max; ++wg) result.insert(wg);
  }
  return result;
}

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

initialize the event

Definition at line 18 of file CSCSimHitMatcher.cc.

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

Referenced by ~CSCSimHitMatcher().

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

float CSCSimHitMatcher::LocalBendingInChamber

(

unsigned int

detid

)

const

Definition at line 159 of file CSCSimHitMatcher.cc.

References CSCDetId::chamber(), hiPixelPairStep_cff::deltaPhi, CSCDetId::endcap(), MuonSimHitMatcher::hitsInDetId(), nLayersWithHitsInChamber(), or, PV3DBase< T, PVType, FrameType >::phi(), CSCDetId::ring(), MuonSimHitMatcher::simHitsMeanPosition(), and CSCDetId::station().

Referenced by ~CSCSimHitMatcher().

                                                                      {
  const CSCDetId cscid(detid);
  if (nLayersWithHitsInChamber(detid) < 6) return -100;
  float phi_layer1 = -10;
  float phi_layer6 = 10;

  if (cscid.station() == 1 and (cscid.ring() == 1 or cscid.ring() == 4)) {
    // phi in layer 1
    const CSCDetId cscid1a(cscid.endcap(), cscid.station(), 4, cscid.chamber(),
                           1);
    const CSCDetId cscid1b(cscid.endcap(), cscid.station(), 1, cscid.chamber(),
                           1);
    const edm::PSimHitContainer& hits1a =
        MuonSimHitMatcher::hitsInDetId(cscid1a.rawId());
    const edm::PSimHitContainer& hits1b =
        MuonSimHitMatcher::hitsInDetId(cscid1b.rawId());
    const GlobalPoint& gp1a =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1a.rawId()));
    const GlobalPoint& gp1b =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1b.rawId()));
    if (!hits1a.empty() and !hits1b.empty())
      phi_layer1 = (gp1a.phi() + gp1b.phi()) / 2.0;
    else if (!hits1a.empty())
      phi_layer1 = gp1a.phi();
    else if (!hits1b.empty())
      phi_layer1 = gp1b.phi();

    // phi in layer 6
    const CSCDetId cscid6a(cscid.endcap(), cscid.station(), 4, cscid.chamber(),
                           6);
    const CSCDetId cscid6b(cscid.endcap(), cscid.station(), 1, cscid.chamber(),
                           6);
    const edm::PSimHitContainer& hits6a =
        MuonSimHitMatcher::hitsInDetId(cscid6a.rawId());
    const edm::PSimHitContainer& hits6b =
        MuonSimHitMatcher::hitsInDetId(cscid6b.rawId());
    const GlobalPoint& gp6a =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6a.rawId()));
    const GlobalPoint& gp6b =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6b.rawId()));
    if (!hits6a.empty() and !hits6b.empty())
      phi_layer6 = (gp6a.phi() + gp6b.phi()) / 2.0;
    else if (!hits6a.empty())
      phi_layer6 = gp6a.phi();
    else if (!hits6b.empty())
      phi_layer6 = gp6b.phi();

  } else {
    // phi in layer 1
    const CSCDetId cscid1(cscid.endcap(), cscid.station(), cscid.ring(),
                          cscid.chamber(), 1);
    const GlobalPoint& gp1 =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid1.rawId()));
    phi_layer1 = gp1.phi();

    // phi in layer 6
    const CSCDetId cscid6(cscid.endcap(), cscid.station(), cscid.ring(),
                          cscid.chamber(), 6);
    const GlobalPoint& gp6 =
        simHitsMeanPosition(MuonSimHitMatcher::hitsInDetId(cscid6.rawId()));
    phi_layer6 = gp6.phi();
  }
  return deltaPhi(phi_layer6, phi_layer1);
}

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

do the matching

Definition at line 32 of file CSCSimHitMatcher.cc.

References MuonSimHitMatcher::detid_to_hits_, detIds(), MuonSimHitMatcher::hasGeometry_, MuonSimHitMatcher::hits_, MuonSimHitMatcher::hitsInDetId(), hitStripsInDetId(), triggerObjects_cff::id, MuonSimHitMatcher::match(), matchSimHitsToSimTrack(), or, AlCaHLTBitMon_ParallelJobs::p, CSCDetId::ring(), MuonSimHitMatcher::simHitsMeanPosition(), CSCDetId::station(), RecoTauPiZeroBuilderPlugins_cfi::strips, MuonSimHitMatcher::track_ids_, and MuonSimHitMatcher::verbose_.

Referenced by ~CSCSimHitMatcher().

                                                                           {
  // instantiates the track ids and simhits
  MuonSimHitMatcher::match(track, vertex);

  if (hasGeometry_) {
    matchSimHitsToSimTrack();

    if (verbose_) {
      edm::LogInfo("CSCSimHitMatcher")
          << "nTrackIds " << track_ids_.size() << " nSelectedCSCSimHits "
          << hits_.size() << endl;
      edm::LogInfo("CSCSimHitMatcher")
          << "detids CSC " << detIds(0).size() << endl;

      for (const auto& id : detIds(0)) {
        const auto& simhits = hitsInDetId(id);
        const auto& simhits_gp = simHitsMeanPosition(simhits);
        const auto& strips = hitStripsInDetId(id);
        CSCDetId cscid(id);
        if (cscid.station() == 1 and (cscid.ring() == 1 or cscid.ring() == 4)) {
          edm::LogInfo("CSCSimHitMatcher")
              << "cscdetid " << CSCDetId(id) << ": " << simhits.size() << " "
              << simhits_gp.phi() << " " << detid_to_hits_[id].size() << endl;
          edm::LogInfo("CSCSimHitMatcher")
              << "nStrip " << strips.size() << endl;
          edm::LogInfo("CSCSimHitMatcher") << "strips : ";
          for (const auto& p : strips) {
            edm::LogInfo("CSCSimHitMatcher") << p;
          }
        }
      }
    }
  }
}

void CSCSimHitMatcher::matchSimHitsToSimTrack ( )

private

Definition at line 67 of file CSCSimHitMatcher.cc.

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

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

                                              {
  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 CSC chambers
      if (discardEleHits_ && pdgid == 11) continue;

      const CSCDetId& layer_id(h.detUnitId());
      hits_.push_back(h);
      detid_to_hits_[h.detUnitId()].push_back(h);
      chamber_to_hits_[layer_id.chamberId().rawId()].push_back(h);
    }
  }
}

int CSCSimHitMatcher::nCoincidenceChambers

(

int

min_n_layers = 4

)

const

Definition at line 317 of file CSCSimHitMatcher.cc.

References chamberIds(), nLayersWithHitsInChamber(), and mps_fire::result.

Referenced by ~CSCSimHitMatcher().

                                                                 {
  int result = 0;
  const auto& chamber_ids = chamberIds(0);
  for (const auto& id : chamber_ids) {
    if (nLayersWithHitsInChamber(id) >= min_n_layers) result += 1;
  }
  return result;
}

int CSCSimHitMatcher::nLayersWithHitsInChamber

(

unsigned int

detid

)

const

Definition at line 112 of file CSCSimHitMatcher.cc.

References h, and MuonSimHitMatcher::hitsInChamber().

Referenced by hitStation(), LocalBendingInChamber(), nCoincidenceChambers(), and ~CSCSimHitMatcher().

                                                                       {
  set<int> layers_with_hits;
  for (const auto& h : MuonSimHitMatcher::hitsInChamber(detid)) {
    const CSCDetId& idd(h.detUnitId());
    layers_with_hits.insert(idd.layer());
  }
  return layers_with_hits.size();
}

int CSCSimHitMatcher::nStations

(

int

nl = 4

)

const

Definition at line 154 of file CSCSimHitMatcher.cc.

References hitStation().

Referenced by ~CSCSimHitMatcher().

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

GlobalPoint CSCSimHitMatcher::simHitPositionKeyLayer

(

unsigned int

chamberid

)

const

Referenced by ~CSCSimHitMatcher().

GlobalPoint CSCSimHitMatcher::simHitsMeanPositionStation

(

int

n

)

const

Referenced by ~CSCSimHitMatcher().

float CSCSimHitMatcher::simHitsMeanStrip

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 224 of file CSCSimHitMatcher.cc.

References edmIntegrityCheck::d, f, MuonSimHitMatcher::geometry_, h, gen::n, and alignCSCRings::s.

Referenced by ~CSCSimHitMatcher().

                                               {
  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 CSCGeometry*>(geometry_)
            ->layer(d)
            ->geometry()
            ->strip(lp);
    // convert to half-strip:
    s *= 2.;
    sums += s;
    ++n;
  }
  if (n == 0) return -1.f;
  return sums / n;
}

float CSCSimHitMatcher::simHitsMeanWG

(

const edm::PSimHitContainer &

sim_hits

)

const

Definition at line 247 of file CSCSimHitMatcher.cc.

References edmIntegrityCheck::d, f, geometry, MuonSimHitMatcher::geometry_, h, gen::n, and alignCSCRings::s.

Referenced by ~CSCSimHitMatcher().

                                               {
  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();
    // find nearest wire
    const auto& layerG(
        dynamic_cast<const CSCGeometry*>(geometry_)->layer(d)->geometry());
    int nearestWire(layerG->nearestWire(lp));
    // then find the corresponding wire group
    s = layerG->wireGroup(nearestWire);
    sums += s;
    ++n;
  }
  if (n == 0) return -1.f;
  return sums / n;
}

Member Data Documentation

edm::ESHandle<CSCGeometry> CSCSimHitMatcher::csc_geom_

private

Definition at line 76 of file CSCSimHitMatcher.h.

Referenced by init().