GEMRecHitMatcher Class Reference

#include <GEMRecHitMatcher.h>

Public Member Functions
bool	areGEMRecHitSame (const GEMRecHit &l, const GEMRecHit &r) const
std::set< unsigned int >	chamberIds () const
std::set< unsigned int >	detIds () const
std::shared_ptr< GEMDigiMatcher >	gemDigiMatcher () const
	GEMRecHitMatcher (edm::ParameterSet const &iPS, edm::ConsumesCollector &&iC)
void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
bool	isGEMRecHitMatched (const GEMRecHit &thisRh) const
void	match (const SimTrack &t, const SimVertex &v)
	do the matching More...
int	nGEMRecHits () const
	How many recHits in GEM did this simtrack get in total? More...
int	nLayersWithRecHitsInSuperChamber (unsigned int) const
std::set< int >	partitionNumbers () const
bool	recHitInContainer (const GEMRecHit &rh, const GEMRecHitContainer &c) const
GlobalPoint	recHitMeanPosition (const GEMRecHitContainer &rechits) const
GlobalPoint	recHitPosition (const GEMRecHit &rechit) const
const GEMRecHitContainer &	recHits () const
const GEMRecHitContainer &	recHitsInChamber (unsigned int) const
const GEMRecHitContainer &	recHitsInDetId (unsigned int) const
const GEMRecHitContainer &	recHitsInSuperChamber (unsigned int) const
std::set< int >	stripNumbersInDetId (unsigned int) const
std::set< unsigned int >	superChamberIds () const
bool	verbose () const
	~GEMRecHitMatcher ()

Private Member Functions
void	matchRecHitsToSimTrack (const GEMRecHitCollection &recHits)

Private Attributes
std::map< unsigned int, GEMRecHitContainer >	chamber_to_recHits_
std::map< unsigned int, GEMRecHitContainer >	detid_to_recHits_
std::shared_ptr< GEMDigiMatcher >	gemDigiMatcher_
const GEMGeometry *	gemGeometry_
edm::Handle< GEMRecHitCollection >	gemRecHitH_
edm::EDGetTokenT< GEMRecHitCollection >	gemRecHitToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_
int	maxBX_
int	minBX_
const GEMRecHitContainer	no_recHits_
GEMRecHitContainer	recHits_
std::map< unsigned int, GEMRecHitContainer >	superchamber_to_recHits_
bool	verbose_

Detailed Description

Description: Matching of RecHits for SimTrack in GEM

Original Author : "Vadim Khotilovich" Contibuting Author : "Claudio Caputo"

Definition at line 30 of file GEMRecHitMatcher.h.

Constructor & Destructor Documentation

GEMRecHitMatcher()

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

Definition at line 7 of file GEMRecHitMatcher.cc.

References muonGEMRecHitPSet::gemRecHit, eostools::move(), and muonDTDigis_cfi::pset.

                                                                                         {
  const auto& gemRecHit = pset.getParameterSet("gemRecHit");
  minBX_ = gemRecHit.getParameter<int>("minBX");
  maxBX_ = gemRecHit.getParameter<int>("maxBX");
  verbose_ = gemRecHit.getParameter<int>("verbose");

  // make a new digi matcher
  gemDigiMatcher_.reset(new GEMDigiMatcher(pset, std::move(iC)));

  gemRecHitToken_ = iC.consumes<GEMRecHitCollection>(gemRecHit.getParameter<edm::InputTag>("inputTag"));
  geomToken_ = iC.esConsumes<GEMGeometry, MuonGeometryRecord>();
}

~GEMRecHitMatcher()

GEMRecHitMatcher::~GEMRecHitMatcher ( )

inline

Definition at line 36 of file GEMRecHitMatcher.h.

{}

Member Function Documentation

areGEMRecHitSame()

bool GEMRecHitMatcher::areGEMRecHitSame	(	const GEMRecHit &	l,
		const GEMRecHit &	r
	)		const

Definition at line 198 of file GEMRecHitMatcher.cc.

References MainPageGenerator::l.

                                                                                    {
  return l.localPosition() == r.localPosition() and l.BunchX() == r.BunchX();
}

chamberIds()

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

(

)

const

Definition at line 96 of file GEMRecHitMatcher.cc.

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

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

detIds()

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

(

void

)

const

Definition at line 89 of file GEMRecHitMatcher.cc.

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

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

gemDigiMatcher()

std::shared_ptr<GEMDigiMatcher> GEMRecHitMatcher::gemDigiMatcher ( ) const

inline

Definition at line 79 of file GEMRecHitMatcher.h.

References gemDigiMatcher_.

{ return gemDigiMatcher_; }

init()

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

Definition at line 20 of file GEMRecHitMatcher.cc.

References edm::EventSetup::getData(), and iEvent.

                                                                               {
  gemDigiMatcher_->init(iEvent, iSetup);

  iEvent.getByToken(gemRecHitToken_, gemRecHitH_);

  gemGeometry_ = &iSetup.getData(geomToken_);
}

isGEMRecHitMatched()

bool GEMRecHitMatcher::isGEMRecHitMatched

(

const GEMRecHit &

thisRh

)

const

Definition at line 194 of file GEMRecHitMatcher.cc.

References FastTrackerRecHitMaskProducer_cfi::recHits.

                                                                       {
  return recHitInContainer(thisRh, recHits());
}

match()

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

do the matching

Definition at line 29 of file GEMRecHitMatcher.cc.

References me0TriggerPseudoDigis_cff::gemRecHits, submitPVValidationJobs::t, and findQualityFiles::v.

                                                                  {
  // match digis first
  gemDigiMatcher_->match(t, v);

  // get the rechit collection
  const GEMRecHitCollection& gemRecHits = *gemRecHitH_.product();

  // now match the rechits
  matchRecHitsToSimTrack(gemRecHits);
}

matchRecHitsToSimTrack()

void GEMRecHitMatcher::matchRecHitsToSimTrack ( const GEMRecHitCollection &  recHits )

private

Definition at line 40 of file GEMRecHitMatcher.cc.

References GEMDetId::chamberId(), filterCSVwithJSON::copy, gather_cfg::cout, ztail::d, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, DetId::rawId(), GEMDetId::superChamberId(), and verbose.

                                                                                {
  // get the matched ids with digis
  const auto& det_ids = gemDigiMatcher_->detIdsDigi();

  // loop on those ids
  for (auto id : det_ids) {
    // now check the digis in this detid
    const auto& hit_digis = gemDigiMatcher_->stripNumbersInDetId(id);
    if (verbose()) {
      cout << "hit_digis_fat ";
      copy(hit_digis.begin(), hit_digis.end(), ostream_iterator<int>(cout, " "));
      cout << endl;
    }

    GEMDetId p_id(id);
    const auto& rechits_in_det = rechits.get(p_id);

    for (auto d = rechits_in_det.first; d != rechits_in_det.second; ++d) {
      if (verbose())
        edm::LogInfo("GEMDigiMatcher") << "recHit " << p_id << " " << *d << endl;

      // check that the rechit is within BX range
      if (d->BunchX() < minBX_ || d->BunchX() > maxBX_)
        continue;

      int firstStrip = d->firstClusterStrip();
      int cls = d->clusterSize();
      bool stripFound = false;

      // check that it matches a strip that was hit by digis from our track
      for (int i = firstStrip; i < (firstStrip + cls); i++) {
        if (hit_digis.find(i) != hit_digis.end())
          stripFound = true;
      }

      // this rechit did not correspond with any previously matched digi
      if (!stripFound)
        continue;
      if (verbose())
        edm::LogInfo("GEMDigiMatcher") << "oki" << endl;

      recHits_.push_back(*d);
      detid_to_recHits_[id].push_back(*d);
      chamber_to_recHits_[p_id.chamberId().rawId()].push_back(*d);
      superchamber_to_recHits_[p_id.superChamberId().rawId()].push_back(*d);
    }
  }
}

nGEMRecHits()

int GEMRecHitMatcher::nGEMRecHits

(

)

const

How many recHits in GEM did this simtrack get in total?

nLayersWithRecHitsInSuperChamber()

int GEMRecHitMatcher::nLayersWithRecHitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 128 of file GEMRecHitMatcher.cc.

References ztail::d, ALCARECOPPSCalTrackBasedSel_cff::detid, and hgcalTBTopologyTester_cfi::layers.

                                                                               {
  set<int> layers;
  for (const auto& d : recHitsInSuperChamber(detid)) {
    layers.insert(d.gemId().layer());
  }
  return layers.size();
}

partitionNumbers()

std::set< int > GEMRecHitMatcher::partitionNumbers

(

)

const

Definition at line 147 of file GEMRecHitMatcher.cc.

References mps_fire::result, and GEMDetId::roll().

                                                     {
  std::set<int> result;

  for (auto id : detIds()) {
    GEMDetId idd(id);
    result.insert(idd.roll());
  }
  return result;
}

recHitInContainer()

bool GEMRecHitMatcher::recHitInContainer	(	const GEMRecHit &	rh,
		const GEMRecHitContainer &	c
	)		const

Definition at line 186 of file GEMRecHitMatcher.cc.

References DummyCfis::c.

                                                                                               {
  bool isSame = false;
  for (const auto& thisRH : c)
    if (areGEMRecHitSame(thisRH, rh))
      isSame = true;
  return isSame;
}

recHitMeanPosition()

GlobalPoint GEMRecHitMatcher::recHitMeanPosition

(

const GEMRecHitContainer &

rechits

)

const

Definition at line 163 of file GEMRecHitMatcher.cc.

References ztail::d, and create_idmaps::n.

                                                                                       {
  GlobalPoint point_zero;
  if (rechit.empty())
    return point_zero;  // point "zero"

  float sumx, sumy, sumz;
  sumx = sumy = sumz = 0.f;
  size_t n = 0;
  for (const auto& d : rechit) {
    GlobalPoint gp = recHitPosition(d);
    if (gp == point_zero)
      continue;

    sumx += gp.x();
    sumy += gp.y();
    sumz += gp.z();
    ++n;
  }
  if (n == 0)
    return GlobalPoint();
  return GlobalPoint(sumx / n, sumy / n, sumz / n);
}

recHitPosition()

GlobalPoint GEMRecHitMatcher::recHitPosition

(

const GEMRecHit &

rechit

)

const

Definition at line 157 of file GEMRecHitMatcher.cc.

References GEMRecHit::gemId(), and GEMRecHit::localPosition().

                                                                          {
  const GEMDetId& idd = rechit.gemId();
  const LocalPoint& lp = rechit.localPosition();
  return gemGeometry_->idToDet(idd)->surface().toGlobal(lp);
}

recHits()

const GEMRecHitContainer& GEMRecHitMatcher::recHits ( void  ) const

inline

Definition at line 54 of file GEMRecHitMatcher.h.

References recHits_.

{ return recHits_; }

recHitsInChamber()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInChamber

(

unsigned int

detid

)

const

Definition at line 116 of file GEMRecHitMatcher.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid.

                                                                                     {
  if (chamber_to_recHits_.find(detid) == chamber_to_recHits_.end())
    return no_recHits_;
  return chamber_to_recHits_.at(detid);
}

recHitsInDetId()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInDetId

(

unsigned int

detid

)

const

Definition at line 110 of file GEMRecHitMatcher.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid.

                                                                                   {
  if (detid_to_recHits_.find(detid) == detid_to_recHits_.end())
    return no_recHits_;
  return detid_to_recHits_.at(detid);
}

recHitsInSuperChamber()

const GEMRecHitContainer & GEMRecHitMatcher::recHitsInSuperChamber

(

unsigned int

detid

)

const

Definition at line 122 of file GEMRecHitMatcher.cc.

References ALCARECOPPSCalTrackBasedSel_cff::detid.

                                                                                          {
  if (superchamber_to_recHits_.find(detid) == superchamber_to_recHits_.end())
    return no_recHits_;
  return superchamber_to_recHits_.at(detid);
}

stripNumbersInDetId()

std::set< int > GEMRecHitMatcher::stripNumbersInDetId

(

unsigned int

detid

)

const

Definition at line 136 of file GEMRecHitMatcher.cc.

References ztail::d, ALCARECOPPSCalTrackBasedSel_cff::detid, and mps_fire::result.

                                                                          {
  set<int> result;
  for (const auto& d : recHitsInDetId(detid)) {
    // loop on all strips hit in this rechit
    for (int iStrip = d.firstClusterStrip(); iStrip < d.firstClusterStrip() + d.clusterSize(); ++iStrip) {
      result.insert(iStrip);
    }
  }
  return result;
}

superChamberIds()

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

(

)

const

Definition at line 103 of file GEMRecHitMatcher.cc.

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

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

verbose()

bool GEMRecHitMatcher::verbose ( ) const

inline

Definition at line 71 of file GEMRecHitMatcher.h.

References verbose_.

{ return verbose_; }

Member Data Documentation

chamber_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::chamber_to_recHits_

private

Definition at line 100 of file GEMRecHitMatcher.h.

detid_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::detid_to_recHits_

private

Definition at line 99 of file GEMRecHitMatcher.h.

gemDigiMatcher_

std::shared_ptr<GEMDigiMatcher> GEMRecHitMatcher::gemDigiMatcher_

private

Definition at line 106 of file GEMRecHitMatcher.h.

Referenced by gemDigiMatcher().

gemGeometry_

const GEMGeometry* GEMRecHitMatcher::gemGeometry_

private

Definition at line 94 of file GEMRecHitMatcher.h.

gemRecHitH_

edm::Handle<GEMRecHitCollection> GEMRecHitMatcher::gemRecHitH_

private

Definition at line 91 of file GEMRecHitMatcher.h.

gemRecHitToken_

edm::EDGetTokenT<GEMRecHitCollection> GEMRecHitMatcher::gemRecHitToken_

private

Definition at line 90 of file GEMRecHitMatcher.h.

geomToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMRecHitMatcher::geomToken_

private

Definition at line 93 of file GEMRecHitMatcher.h.

maxBX_

int GEMRecHitMatcher::maxBX_

private

Definition at line 96 of file GEMRecHitMatcher.h.

minBX_

int GEMRecHitMatcher::minBX_

private

Definition at line 96 of file GEMRecHitMatcher.h.

no_recHits_

const GEMRecHitContainer GEMRecHitMatcher::no_recHits_

private

Definition at line 103 of file GEMRecHitMatcher.h.

recHits_

GEMRecHitContainer GEMRecHitMatcher::recHits_

private

Definition at line 104 of file GEMRecHitMatcher.h.

Referenced by recHits().

superchamber_to_recHits_

std::map<unsigned int, GEMRecHitContainer> GEMRecHitMatcher::superchamber_to_recHits_

private

Definition at line 101 of file GEMRecHitMatcher.h.

verbose_

bool GEMRecHitMatcher::verbose_

private

Definition at line 97 of file GEMRecHitMatcher.h.

Referenced by verbose().