CTPPSDiamondTrackRecognition Class Reference

Class performing smart reconstruction for PPS Diamond Detectors. More...

#include <CTPPSDiamondTrackRecognition.h>

Inheritance diagram for CTPPSDiamondTrackRecognition:

Public Member Functions
void	addHit (const CTPPSDiamondRecHit &recHit) override
	Feed a new hit to the tracks recognition algorithm. More...
void	clear () override
	Reset internal state of a class instance. More...
	CTPPSDiamondTrackRecognition (const edm::ParameterSet &iConfig)
int	produceTracks (edm::DetSet< CTPPSDiamondLocalTrack > &tracks) override
	Produce a collection of tracks for the current station, given its hits collection. More...
Public Member Functions inherited from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >
	CTPPSTimingTrackRecognition (const edm::ParameterSet &iConfig)
virtual	~CTPPSTimingTrackRecognition ()=default

Private Attributes
bool	excludeSingleEdgeHits_
std::unordered_map< int, int >	mhMap_

Additional Inherited Members
Protected Types inherited from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >
typedef std::vector< CTPPSDiamondRecHit >	HitVector
typedef std::unordered_map< int, HitVector >	HitVectorMap
typedef std::vector< CTPPSDiamondLocalTrack >	TrackVector
Protected Member Functions inherited from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >
SpatialRange	getHitSpatialRange (const HitVector &hits)
void	producePartialTracks (const HitVector &hits, const DimensionParameters &param, float(*getHitCenter)(const CTPPSDiamondRecHit &), float(*getHitRangeWidth)(const CTPPSDiamondRecHit &), void(*setTrackCenter)(CTPPSDiamondLocalTrack &, float), void(*setTrackSigma)(CTPPSDiamondLocalTrack &, float), TrackVector &result)
bool	timeEval (const HitVector &hits, float &meanTime, float &timeSigma) const
Protected Attributes inherited from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >
HitVectorMap	hitVectorMap_
	RecHit vectors that should be processed separately while reconstructing tracks. More...
reco::FormulaEvaluator	pixelEfficiencyFunction_
const float	resolution_
const float	sigma_
const float	threshold_
const float	thresholdFromMaximum_
const float	tolerance_

Detailed Description

Class performing smart reconstruction for PPS Diamond Detectors.

Date

Jan 2017

Definition at line 26 of file CTPPSDiamondTrackRecognition.h.

Constructor & Destructor Documentation

CTPPSDiamondTrackRecognition::CTPPSDiamondTrackRecognition

(

const edm::ParameterSet &

iConfig

)

Definition at line 15 of file CTPPSDiamondTrackRecognition.cc.

    : CTPPSTimingTrackRecognition<CTPPSDiamondLocalTrack, CTPPSDiamondRecHit>(iConfig),
      excludeSingleEdgeHits_(iConfig.getParameter<bool>("excludeSingleEdgeHits")) {}

Member Function Documentation

void CTPPSDiamondTrackRecognition::addHit ( const CTPPSDiamondRecHit &  recHit )

overridevirtual

Feed a new hit to the tracks recognition algorithm.

Implements CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 28 of file CTPPSDiamondTrackRecognition.cc.

References excludeSingleEdgeHits_, CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::hitVectorMap_, CTPPSDiamondRecHit::ootIndex(), and CTPPSDiamondRecHit::toT().

Referenced by CTPPSDiamondLocalTrackFitter::produce().

                                                                          {
  if (excludeSingleEdgeHits_ && recHit.toT() <= 0.)
    return;
  // store hit parameters
  hitVectorMap_[recHit.ootIndex()].emplace_back(recHit);
}

void CTPPSDiamondTrackRecognition::clear ( void  )

overridevirtual

Reset internal state of a class instance.

Reimplemented from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 21 of file CTPPSDiamondTrackRecognition.cc.

References CTPPSTimingTrackRecognition< TRACK_TYPE, HIT_TYPE >::clear(), and mhMap_.

Referenced by CTPPSDiamondLocalTrackFitter::produce().

                                         {
  CTPPSTimingTrackRecognition<CTPPSDiamondLocalTrack, CTPPSDiamondRecHit>::clear();
  mhMap_.clear();
}

int CTPPSDiamondTrackRecognition::produceTracks ( edm::DetSet< CTPPSDiamondLocalTrack > &  tracks )

overridevirtual

Produce a collection of tracks for the current station, given its hits collection.

Implements CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 37 of file CTPPSDiamondTrackRecognition.cc.

References CTPPSDiamondLocalTrack::containsHit(), excludeSingleEdgeHits_, f, CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::getHitSpatialRange(), hfClusterShapes_cfi::hits, CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::hitVectorMap_, mhMap_, position, CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::producePartialTracks(), edm::DetSet< T >::push_back(), CTPPSTimingLocalTrack::setPosition(), CTPPSTimingLocalTrack::setPositionSigma(), CTPPSTimingLocalTrack::setTime(), CTPPSTimingLocalTrack::setTimeSigma(), CTPPSTimingLocalTrack::setValid(), CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::timeEval(), CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::tolerance_, HLT_2018_cff::track, x, and hit::x.

Referenced by CTPPSDiamondLocalTrackFitter::produce().

                                                                                         {
  int numberOfTracks = 0;
  DimensionParameters param;

  auto getX = [](const CTPPSDiamondRecHit& hit) { return hit.x(); };
  auto getXWidth = [](const CTPPSDiamondRecHit& hit) { return hit.xWidth(); };
  auto setX = [](CTPPSDiamondLocalTrack& track, float x) { track.setPosition(math::XYZPoint(x, 0., 0.)); };
  auto setXSigma = [](CTPPSDiamondLocalTrack& track, float sigma) {
    track.setPositionSigma(math::XYZPoint(sigma, 0., 0.));
  };

  for (const auto& hitBatch : hitVectorMap_) {
    // separate the tracking for each bunch crossing
    const auto& oot = hitBatch.first;
    const auto& hits = hitBatch.second;

    auto hitRange = getHitSpatialRange(hits);

    TrackVector xPartTracks;

    // produce tracks in x dimension
    param.rangeBegin = hitRange.xBegin;
    param.rangeEnd = hitRange.xEnd;
    producePartialTracks(hits, param, getX, getXWidth, setX, setXSigma, xPartTracks);

    if (xPartTracks.empty())
      continue;

    const float yRangeCenter = 0.5f * (hitRange.yBegin + hitRange.yEnd);
    const float zRangeCenter = 0.5f * (hitRange.zBegin + hitRange.zEnd);
    const float ySigma = 0.5f * (hitRange.yEnd - hitRange.yBegin);
    const float zSigma = 0.5f * (hitRange.zEnd - hitRange.zBegin);

    for (const auto& xTrack : xPartTracks) {
      math::XYZPoint position(xTrack.x0(), yRangeCenter, zRangeCenter);
      math::XYZPoint positionSigma(xTrack.x0Sigma(), ySigma, zSigma);

      const int multipleHits = (mhMap_.find(oot) != mhMap_.end()) ? mhMap_[oot] : 0;
      CTPPSDiamondLocalTrack newTrack(position, positionSigma, 0.f, 0.f, oot, multipleHits);

      // find contributing hits
      HitVector componentHits;
      for (const auto& hit : hits)
        if (newTrack.containsHit(hit, tolerance_) && (!excludeSingleEdgeHits_ || hit.toT() > 0.))
          componentHits.emplace_back(hit);
      // compute timing information
      float mean_time = 0.f, time_sigma = 0.f;
      bool valid_hits = timeEval(componentHits, mean_time, time_sigma);
      newTrack.setValid(valid_hits);
      newTrack.setTime(mean_time);
      newTrack.setTimeSigma(time_sigma);

      tracks.push_back(newTrack);
    }
  }

  return numberOfTracks;
}

Member Data Documentation

bool CTPPSDiamondTrackRecognition::excludeSingleEdgeHits_

private

Definition at line 38 of file CTPPSDiamondTrackRecognition.h.

Referenced by addHit(), and produceTracks().

std::unordered_map<int, int> CTPPSDiamondTrackRecognition::mhMap_

private

Definition at line 37 of file CTPPSDiamondTrackRecognition.h.

Referenced by clear(), and produceTracks().