CSCRecHitMatcher Class Reference

#include <CSCRecHitMatcher.h>

Public Member Functions
bool	areCSCRecHit2DsSame (const CSCRecHit2D &, const CSCRecHit2D &) const
bool	areCSCSegmentsSame (const CSCSegment &, const CSCSegment &) const
CSCSegment	bestCSCSegment (unsigned int)
std::set< unsigned int >	chamberIdsCSCRecHit2D () const
std::set< unsigned int >	chamberIdsCSCSegment () const
bool	cscRecHit2DInContainer (const CSCRecHit2D &, const CSCRecHit2DContainer &) const
const CSCRecHit2DContainer	cscRecHit2Ds () const
const CSCRecHit2DContainer &	cscRecHit2DsInChamber (unsigned int) const
const CSCRecHit2DContainer &	cscRecHit2DsInLayer (unsigned int) const
	CSCRecHitMatcher (edm::ParameterSet const &iPS, edm::ConsumesCollector &&iC)
bool	cscSegmentInContainer (const CSCSegment &, const CSCSegmentContainer &) const
const CSCSegmentContainer	cscSegments () const
const CSCSegmentContainer &	cscSegmentsInChamber (unsigned int) const
GlobalPoint	globalPoint (const CSCSegment &) const
void	init (const edm::Event &e, const edm::EventSetup &eventSetup)
bool	isCSCRecHit2DMatched (const CSCRecHit2D &) const
bool	isCSCSegmentMatched (const CSCSegment &) const
std::set< unsigned int >	layerIdsCSCRecHit2D () const
void	match (const SimTrack &t, const SimVertex &v)
	do the matching More...
int	nCSCRecHit2Ds () const
int	nCSCRecHit2DsInChamber (unsigned int) const
int	nCSCRecHit2DsInLayer (unsigned int) const
int	nCSCSegments () const
int	nCSCSegmentsInChamber (unsigned int) const
	~CSCRecHitMatcher ()

Private Member Functions
void	matchCSCRecHit2DsToSimTrack (const CSCRecHit2DCollection &)
void	matchCSCSegmentsToSimTrack (const CSCSegmentCollection &)

Private Attributes
std::map< unsigned int, CSCRecHit2DContainer >	chamber_to_cscRecHit2D_
std::map< unsigned int, CSCSegmentContainer >	chamber_to_cscSegment_
edm::ESHandle< CSCGeometry >	csc_geom_
std::unique_ptr< CSCDigiMatcher >	cscDigiMatcher_
const CSCGeometry *	cscGeometry_
edm::Handle < CSCRecHit2DCollection >	cscRecHit2DH_
edm::EDGetTokenT < CSCRecHit2DCollection >	cscRecHit2DToken_
edm::Handle< CSCSegmentCollection >	cscSegmentH_
edm::EDGetTokenT < CSCSegmentCollection >	cscSegmentToken_
std::map< unsigned int, CSCRecHit2DContainer >	layer_to_cscRecHit2D_
int	maxBXCSCRecHit2D_
int	maxBXCSCSegment_
int	minBXCSCRecHit2D_
int	minBXCSCSegment_
CSCRecHit2DContainer	no_cscRecHit2Ds_
CSCSegmentContainer	no_cscSegments_
int	verboseCSCRecHit2D_
int	verboseCSCSegment_

Detailed Description

Definition at line 31 of file CSCRecHitMatcher.h.

Constructor & Destructor Documentation

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

Definition at line 9 of file CSCRecHitMatcher.cc.

References edm::ParameterSet::getParameter(), eostools::move(), and TrackValidation_cff::pset.

                                                                                         {
  const auto& cscRecHit2D = pset.getParameter<edm::ParameterSet>("cscRecHit");
  maxBXCSCRecHit2D_ = cscRecHit2D.getParameter<int>("maxBX");
  minBXCSCRecHit2D_ = cscRecHit2D.getParameter<int>("minBX");
  verboseCSCRecHit2D_ = cscRecHit2D.getParameter<int>("verbose");

  const auto& cscSegment = pset.getParameter<edm::ParameterSet>("cscSegment");
  maxBXCSCSegment_ = cscSegment.getParameter<int>("maxBX");
  minBXCSCSegment_ = cscSegment.getParameter<int>("minBX");
  verboseCSCSegment_ = cscSegment.getParameter<int>("verbose");

  // make a new digi matcher
  cscDigiMatcher_ = std::make_unique<CSCDigiMatcher>(pset, std::move(iC));

  cscRecHit2DToken_ = iC.consumes<CSCRecHit2DCollection>(cscRecHit2D.getParameter<edm::InputTag>("inputTag"));
  cscSegmentToken_ = iC.consumes<CSCSegmentCollection>(cscSegment.getParameter<edm::InputTag>("inputTag"));
}

CSCRecHitMatcher::~CSCRecHitMatcher ( )

inline

Definition at line 37 of file CSCRecHitMatcher.h.

{}

Member Function Documentation

bool CSCRecHitMatcher::areCSCRecHit2DsSame	(	const CSCRecHit2D &	l,
		const CSCRecHit2D &	r
	)		const

Definition at line 271 of file CSCRecHitMatcher.cc.

References CSCRecHit2D::localPosition().

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

bool CSCRecHitMatcher::areCSCSegmentsSame	(	const CSCSegment &	l,
		const CSCSegment &	r
	)		const

Definition at line 275 of file CSCRecHitMatcher.cc.

References CSCSegment::localDirection(), and CSCSegment::localPosition().

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

CSCSegment CSCRecHitMatcher::bestCSCSegment

(

unsigned int

id

)

Definition at line 279 of file CSCRecHitMatcher.cc.

References gpuClustering::id.

                                                           {
  CSCSegment emptySegment;
  double chi2overNdf = 99;
  int index = 0;
  int foundIndex = -99;

  for (const auto& seg : chamber_to_cscSegment_[id]) {
    double newChi2overNdf(seg.chi2() / seg.degreesOfFreedom());
    if (newChi2overNdf < chi2overNdf) {
      chi2overNdf = newChi2overNdf;
      foundIndex = index;
    }
    ++index;
  }
  if (foundIndex != -99)
    return chamber_to_cscSegment_[id][foundIndex];
  return emptySegment;
}

std::set< unsigned int > CSCRecHitMatcher::chamberIdsCSCRecHit2D

(

)

const

Definition at line 175 of file CSCRecHitMatcher.cc.

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

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

std::set< unsigned int > CSCRecHitMatcher::chamberIdsCSCSegment

(

)

const

Definition at line 182 of file CSCRecHitMatcher.cc.

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

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

bool CSCRecHitMatcher::cscRecHit2DInContainer	(	const CSCRecHit2D &	sg,
		const CSCRecHit2DContainer &	c
	)		const

Definition at line 231 of file CSCRecHitMatcher.cc.

                                                                                                        {
  bool isSame = false;
  for (const auto& segment : c)
    if (areCSCRecHit2DsSame(sg, segment))
      isSame = true;
  return isSame;
}

const CSCRecHit2DContainer CSCRecHitMatcher::cscRecHit2Ds

(

)

const

Definition at line 213 of file CSCRecHitMatcher.cc.

References mps_fire::result.

                                                                {
  CSCRecHit2DContainer result;
  for (const auto& id : chamberIdsCSCRecHit2D()) {
    const auto& segmentsInChamber(cscRecHit2DsInChamber(id));
    result.insert(result.end(), segmentsInChamber.begin(), segmentsInChamber.end());
  }
  return result;
}

const CSCRecHit2DContainer & CSCRecHitMatcher::cscRecHit2DsInChamber

(

unsigned int

detid

)

const

Definition at line 195 of file CSCRecHitMatcher.cc.

                                                                                            {
  if (chamber_to_cscRecHit2D_.find(detid) == chamber_to_cscRecHit2D_.end())
    return no_cscRecHit2Ds_;
  return chamber_to_cscRecHit2D_.at(detid);
}

const CSCRecHit2DContainer & CSCRecHitMatcher::cscRecHit2DsInLayer

(

unsigned int

detid

)

const

Definition at line 189 of file CSCRecHitMatcher.cc.

                                                                                          {
  if (layer_to_cscRecHit2D_.find(detid) == layer_to_cscRecHit2D_.end())
    return no_cscRecHit2Ds_;
  return layer_to_cscRecHit2D_.at(detid);
}

bool CSCRecHitMatcher::cscSegmentInContainer	(	const CSCSegment &	sg,
		const CSCSegmentContainer &	c
	)		const

Definition at line 239 of file CSCRecHitMatcher.cc.

                                                                                                     {
  bool isSame = false;
  for (const auto& segment : c)
    if (areCSCSegmentsSame(sg, segment))
      isSame = true;
  return isSame;
}

const CSCSegmentContainer CSCRecHitMatcher::cscSegments

(

)

const

Definition at line 222 of file CSCRecHitMatcher.cc.

References mps_fire::result.

                                                              {
  CSCSegmentContainer result;
  for (const auto& id : chamberIdsCSCSegment()) {
    const auto& segmentsInChamber(cscSegmentsInChamber(id));
    result.insert(result.end(), segmentsInChamber.begin(), segmentsInChamber.end());
  }
  return result;
}

const CSCSegmentContainer & CSCRecHitMatcher::cscSegmentsInChamber

(

unsigned int

detid

)

const

Definition at line 201 of file CSCRecHitMatcher.cc.

                                                                                          {
  if (chamber_to_cscSegment_.find(detid) == chamber_to_cscSegment_.end())
    return no_cscSegments_;
  return chamber_to_cscSegment_.at(detid);
}

GlobalPoint CSCRecHitMatcher::globalPoint

(

const CSCSegment &

c

)

const

Definition at line 298 of file CSCRecHitMatcher.cc.

References CSCSegment::cscDetId(), and CSCSegment::localPosition().

                                                                   {
  return cscGeometry_->idToDet(c.cscDetId())->surface().toGlobal(c.localPosition());
}

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

Definition at line 27 of file CSCRecHitMatcher.cc.

References edm::EventSetup::get(), and edm::Event::getByToken().

                                                                               {
  cscDigiMatcher_->init(iEvent, iSetup);

  iEvent.getByToken(cscRecHit2DToken_, cscRecHit2DH_);
  iEvent.getByToken(cscSegmentToken_, cscSegmentH_);

  iSetup.get<MuonGeometryRecord>().get(csc_geom_);
  if (csc_geom_.isValid()) {
    cscGeometry_ = &*csc_geom_;
  } else {
    edm::LogWarning("CSCSimHitMatcher") << "+++ Info: CSC geometry is unavailable. +++\n";
  }
}

bool CSCRecHitMatcher::isCSCRecHit2DMatched

(

const CSCRecHit2D &

thisSg

)

const

Definition at line 247 of file CSCRecHitMatcher.cc.

                                                                           {
  return cscRecHit2DInContainer(thisSg, cscRecHit2Ds());
}

bool CSCRecHitMatcher::isCSCSegmentMatched

(

const CSCSegment &

thisSg

)

const

Definition at line 251 of file CSCRecHitMatcher.cc.

References cscSegments_cfi::cscSegments.

                                                                         {
  return cscSegmentInContainer(thisSg, cscSegments());
}

std::set< unsigned int > CSCRecHitMatcher::layerIdsCSCRecHit2D

(

)

const

Definition at line 168 of file CSCRecHitMatcher.cc.

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

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

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

do the matching

Definition at line 42 of file CSCRecHitMatcher.cc.

References cscSegments_cfi::cscSegments.

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

  // get the rechit collection
  const CSCRecHit2DCollection& cscRecHit2Ds = *cscRecHit2DH_.product();
  const CSCSegmentCollection& cscSegments = *cscSegmentH_.product();

  // now match the rechits
  matchCSCRecHit2DsToSimTrack(cscRecHit2Ds);
  matchCSCSegmentsToSimTrack(cscSegments);
}

void CSCRecHitMatcher::matchCSCRecHit2DsToSimTrack ( const CSCRecHit2DCollection &  rechits )

private

Definition at line 55 of file CSCRecHitMatcher.cc.

References CSCDetId::chamberId(), ztail::d, spr::find(), gpuClustering::id, AlCaHLTBitMon_ParallelJobs::p, and DetId::rawId().

                                                                                       {
  if (verboseCSCRecHit2D_)
    edm::LogInfo("CSCRecHitMatcher") << "Matching simtrack to CSC rechits";

  // fetch all layerIds with digis
  const auto& strip_ids = cscDigiMatcher_->detIdsStrip();
  const auto& wire_ids = cscDigiMatcher_->detIdsWire();

  // merge the two collections
  std::set<unsigned int> layer_ids;
  layer_ids.insert(strip_ids.begin(), strip_ids.end());
  layer_ids.insert(wire_ids.begin(), wire_ids.end());

  if (verboseCSCRecHit2D_)
    edm::LogInfo("CSCRecHitMatcher") << "Number of matched csc layer_ids " << layer_ids.size();

  for (const auto& id : layer_ids) {
    CSCDetId p_id(id);

    // print all the wires in the CSCChamber
    const auto& hit_wg(cscDigiMatcher_->wiregroupsInDetId(id));
    if (verboseCSCRecHit2D_) {
      edm::LogInfo("CSCRecHitMatcher") << "hit wg csc from simhit" << endl;
      for (const auto& p : hit_wg) {
        edm::LogInfo("CSCRecHitMatcher") << p;
      }
    }

    // print all the strips in the CSCChamber
    const auto& hit_strips(cscDigiMatcher_->stripsInDetId(id));
    if (verboseCSCRecHit2D_) {
      edm::LogInfo("CSCRecHitMatcher") << "hit strip csc from simhit" << endl;
      for (const auto& p : hit_strips) {
        edm::LogInfo("CSCRecHitMatcher") << p;
      }
    }

    // get the rechits
    const auto& rechits_in_det = rechits.get(p_id);
    for (auto d = rechits_in_det.first; d != rechits_in_det.second; ++d) {
      if (verboseCSCRecHit2D_)
        edm::LogInfo("CSCRecHitMatcher") << "rechit " << p_id << " " << *d;

      // does the wire number match?
      const bool wireMatch(std::find(hit_wg.begin(), hit_wg.end(), d->hitWire()) != hit_wg.end());

      // does the strip number match?
      bool stripMatch(false);
      for (size_t iS = 0; iS < d->nStrips(); ++iS) {
        if (std::find(hit_strips.begin(), hit_strips.end(), d->channels(iS)) != hit_strips.end())
          stripMatch = true;
      }

      // this rechit was matched to a matching simhit
      if (wireMatch and stripMatch) {
        if (verboseCSCRecHit2D_)
          edm::LogInfo("CSCRecHitMatcher") << "\t...was matched!";
        layer_to_cscRecHit2D_[id].push_back(*d);
        chamber_to_cscRecHit2D_[p_id.chamberId().rawId()].push_back(*d);
      }
    }
  }
}

void CSCRecHitMatcher::matchCSCSegmentsToSimTrack ( const CSCSegmentCollection &  cscSegments )

private

Definition at line 119 of file CSCRecHitMatcher.cc.

References ztail::d, AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), and HLT_FULL_cff::recHits.

                                                                                         {
  if (verboseCSCSegment_)
    edm::LogInfo("CSCRecHitMatcher") << "Matching simtrack to segments";
  // fetch all chamberIds with 2D rechits

  const auto& chamber_ids = chamberIdsCSCRecHit2D();
  if (verboseCSCSegment_)
    edm::LogInfo("CSCRecHitMatcher") << "Number of matched csc segments " << chamber_ids.size();
  for (const auto& id : chamber_ids) {
    CSCDetId p_id(id);

    // print all CSCRecHit2D in the CSCChamber
    const auto& csc_rechits(cscRecHit2DsInChamber(id));
    if (verboseCSCSegment_) {
      edm::LogInfo("CSCRecHitMatcher") << "hit csc rechits" << endl;
      for (const auto& p : csc_rechits) {
        edm::LogInfo("CSCRecHitMatcher") << p;
      }
    }

    // get the segments
    const auto& segments_in_det = cscSegments.get(p_id);
    for (auto d = segments_in_det.first; d != segments_in_det.second; ++d) {
      if (verboseCSCSegment_)
        edm::LogInfo("CSCRecHitMatcher") << "segment " << p_id << " " << *d << endl;

      //access the rechits
      const auto& recHits(d->recHits());

      int rechitsFound = 0;
      if (verboseCSCSegment_)
        edm::LogInfo("CSCRecHitMatcher") << recHits.size() << " csc rechits from segment " << endl;
      for (const auto& rh : recHits) {
        const CSCRecHit2D* cscrh(dynamic_cast<const CSCRecHit2D*>(rh));
        if (isCSCRecHit2DMatched(*cscrh))
          ++rechitsFound;
      }
      if (rechitsFound == 0)
        continue;
      if (verboseCSCSegment_) {
        edm::LogInfo("CSCRecHitMatcher") << "Found " << rechitsFound << " rechits out of "
                                         << cscRecHit2DsInChamber(id).size();
        edm::LogInfo("CSCRecHitMatcher") << "\t...was matched!";
      }
      chamber_to_cscSegment_[p_id.rawId()].push_back(*d);
    }
  }
}

int CSCRecHitMatcher::nCSCRecHit2Ds

(

)

const

Definition at line 255 of file CSCRecHitMatcher.cc.

References dqmiodumpmetadata::n.

                                          {
  int n = 0;
  const auto& ids = chamberIdsCSCRecHit2D();
  for (const auto& id : ids)
    n += cscRecHit2DsInChamber(id).size();
  return n;
}

int CSCRecHitMatcher::nCSCRecHit2DsInChamber

(

unsigned int

detid

)

const

Definition at line 209 of file CSCRecHitMatcher.cc.

{ return cscRecHit2DsInChamber(detid).size(); }

int CSCRecHitMatcher::nCSCRecHit2DsInLayer

(

unsigned int

detid

)

const

Definition at line 207 of file CSCRecHitMatcher.cc.

{ return cscRecHit2DsInLayer(detid).size(); }

int CSCRecHitMatcher::nCSCSegments

(

)

const

Definition at line 263 of file CSCRecHitMatcher.cc.

References dqmiodumpmetadata::n.

                                         {
  int n = 0;
  const auto& ids = chamberIdsCSCSegment();
  for (const auto& id : ids)
    n += cscSegmentsInChamber(id).size();
  return n;
}

int CSCRecHitMatcher::nCSCSegmentsInChamber

(

unsigned int

detid

)

const

Definition at line 211 of file CSCRecHitMatcher.cc.

{ return cscSegmentsInChamber(detid).size(); }

Member Data Documentation

std::map<unsigned int, CSCRecHit2DContainer> CSCRecHitMatcher::chamber_to_cscRecHit2D_

private

Definition at line 106 of file CSCRecHitMatcher.h.

std::map<unsigned int, CSCSegmentContainer> CSCRecHitMatcher::chamber_to_cscSegment_

private

Definition at line 107 of file CSCRecHitMatcher.h.

edm::ESHandle<CSCGeometry> CSCRecHitMatcher::csc_geom_

private

Definition at line 91 of file CSCRecHitMatcher.h.

std::unique_ptr<CSCDigiMatcher> CSCRecHitMatcher::cscDigiMatcher_

private

Definition at line 83 of file CSCRecHitMatcher.h.

const CSCGeometry* CSCRecHitMatcher::cscGeometry_

private

Definition at line 92 of file CSCRecHitMatcher.h.

edm::Handle<CSCRecHit2DCollection> CSCRecHitMatcher::cscRecHit2DH_

private

Definition at line 88 of file CSCRecHitMatcher.h.

edm::EDGetTokenT<CSCRecHit2DCollection> CSCRecHitMatcher::cscRecHit2DToken_

private

Definition at line 85 of file CSCRecHitMatcher.h.

edm::Handle<CSCSegmentCollection> CSCRecHitMatcher::cscSegmentH_

private

Definition at line 89 of file CSCRecHitMatcher.h.

edm::EDGetTokenT<CSCSegmentCollection> CSCRecHitMatcher::cscSegmentToken_

private

Definition at line 86 of file CSCRecHitMatcher.h.

std::map<unsigned int, CSCRecHit2DContainer> CSCRecHitMatcher::layer_to_cscRecHit2D_

private

Definition at line 105 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::maxBXCSCRecHit2D_

private

Definition at line 98 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::maxBXCSCSegment_

private

Definition at line 102 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::minBXCSCRecHit2D_

private

Definition at line 99 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::minBXCSCSegment_

private

Definition at line 103 of file CSCRecHitMatcher.h.

CSCRecHit2DContainer CSCRecHitMatcher::no_cscRecHit2Ds_

private

Definition at line 109 of file CSCRecHitMatcher.h.

CSCSegmentContainer CSCRecHitMatcher::no_cscSegments_

private

Definition at line 110 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::verboseCSCRecHit2D_

private

Definition at line 97 of file CSCRecHitMatcher.h.

int CSCRecHitMatcher::verboseCSCSegment_

private

Definition at line 101 of file CSCRecHitMatcher.h.