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

References excludeSingleEdgeHits_, CTPPSDiamondRecHit::getOOTIndex(), CTPPSDiamondRecHit::getToT(), and CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::hitVectorMap_.

Referenced by CTPPSDiamondLocalTrackFitter::produce().

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

void CTPPSDiamondTrackRecognition::clear ( void  )

overridevirtual

Reset internal state of a class instance.

Reimplemented from CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >.

Definition at line 23 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 43 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::setT(), CTPPSTimingLocalTrack::setTSigma(), CTPPSTimingLocalTrack::setValid(), CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::timeEval(), CTPPSTimingTrackRecognition< CTPPSDiamondLocalTrack, CTPPSDiamondRecHit >::tolerance_, HiIsolationCommonParameters_cff::track, and x.

Referenced by CTPPSDiamondLocalTrackFitter::produce().

{
  int numberOfTracks = 0;
  DimensionParameters param;

  auto getX = []( const CTPPSDiamondRecHit& hit ){ return hit.getX(); };
  auto getXWidth = []( const CTPPSDiamondRecHit& hit ){ return hit.getXWidth(); };
  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.getX0(), yRangeCenter, zRangeCenter );
      math::XYZPoint positionSigma( xTrack.getX0Sigma(), 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.getToT() > 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.setT( mean_time );
      newTrack.setTSigma( time_sigma );

      tracks.push_back( newTrack );
    }
  }

  return numberOfTracks;
}

Member Data Documentation

bool CTPPSDiamondTrackRecognition::excludeSingleEdgeHits_

private

Definition at line 39 of file CTPPSDiamondTrackRecognition.h.

Referenced by addHit(), and produceTracks().

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

private

Definition at line 38 of file CTPPSDiamondTrackRecognition.h.

Referenced by clear(), and produceTracks().