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::CTPPSDiamondTrackRecognition ( const edm::ParameterSet & iConfig )

Definition at line 14 of file CTPPSDiamondTrackRecognition.cc.

16 : CTPPSTimingTrackRecognition<CTPPSDiamondLocalTrack, CTPPSDiamondRecHit>(iConfig),

Member Function Documentation

◆ addHit()

void CTPPSDiamondTrackRecognition::addHit ( const CTPPSDiamondRecHit & recHit )

overridevirtual

Feed a new hit to the tracks recognition algorithm.

Implements CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 27 of file CTPPSDiamondTrackRecognition.cc.

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

References excludeSingleEdgeHits_, CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::hitVectorMap_, and rpcPointValidation_cfi::recHit.

Referenced by CTPPSDiamondLocalTrackFitter::produce().

◆ clear()

void CTPPSDiamondTrackRecognition::clear ( void )

overridevirtual

Reset internal state of a class instance.

Reimplemented from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 20 of file CTPPSDiamondTrackRecognition.cc.

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

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

Referenced by CTPPSDiamondLocalTrackFitter::produce().

◆ produceTracks()

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 36 of file CTPPSDiamondTrackRecognition.cc.

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