TotemTimingTrackRecognition Class Reference

#include <TotemTimingTrackRecognition.h>

Inheritance diagram for TotemTimingTrackRecognition:

Public Member Functions
void	addHit (const TotemTimingRecHit &recHit) override
	Add new hit to the set from which the tracks are reconstructed. More...
int	produceTracks (edm::DetSet< TotemTimingLocalTrack > &tracks) override
	Produces a collection of tracks for the current station, given its hits collection. More...
	TotemTimingTrackRecognition (const edm::ParameterSet &iConfig)
Public Member Functions inherited from CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >
virtual void	clear ()
	Reset internal state of a class instance. More...
	CTPPSTimingTrackRecognition (const edm::ParameterSet &iConfig)
virtual	~CTPPSTimingTrackRecognition ()=default

Additional Inherited Members
Protected Types inherited from CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >
typedef std::vector < TotemTimingRecHit >	HitVector
typedef std::unordered_map < int, HitVector >	HitVectorMap
typedef std::vector < TotemTimingLocalTrack >	TrackVector
Protected Member Functions inherited from CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >
SpatialRange	getHitSpatialRange (const HitVector &hits)
void	producePartialTracks (const HitVector &hits, const DimensionParameters &param, float(*getHitCenter)(const TotemTimingRecHit &), float(*getHitRangeWidth)(const TotemTimingRecHit &), void(*setTrackCenter)(TotemTimingLocalTrack &, float), void(*setTrackSigma)(TotemTimingLocalTrack &, float), TrackVector &result)
bool	timeEval (const HitVector &hits, float &meanTime, float &timeSigma) const
Protected Attributes inherited from CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >
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 intended to perform general CTPPS timing detectors track recognition, as well as construction of specialized classes (for now CTPPSDiamond and TotemTiming local tracks).

Definition at line 24 of file TotemTimingTrackRecognition.h.

Constructor & Destructor Documentation

TotemTimingTrackRecognition::TotemTimingTrackRecognition

(

const edm::ParameterSet &

iConfig

)

Definition at line 15 of file TotemTimingTrackRecognition.cc.

    : CTPPSTimingTrackRecognition<TotemTimingLocalTrack, TotemTimingRecHit>(iConfig) {}

Member Function Documentation

void TotemTimingTrackRecognition::addHit ( const TotemTimingRecHit &  recHit )

overridevirtual

Add new hit to the set from which the tracks are reconstructed.

Implements CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >.

Definition at line 20 of file TotemTimingTrackRecognition.cc.

References CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::hitVectorMap_, TotemTimingRecHit::NO_T_AVAILABLE, and CTPPSTimingRecHit::time().

                                                                        {
  if (recHit.time() != TotemTimingRecHit::NO_T_AVAILABLE)
    hitVectorMap_[0].emplace_back(recHit);
}

int TotemTimingTrackRecognition::produceTracks ( edm::DetSet< TotemTimingLocalTrack > &  tracks )

overridevirtual

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

Implements CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >.

Definition at line 27 of file TotemTimingTrackRecognition.cc.

References CTPPSTimingLocalTrack::containsHit(), CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::getHitSpatialRange(), CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::hitVectorMap_, position, CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::producePartialTracks(), edm::DetSet< T >::push_back(), CTPPSTimingLocalTrack::setPosition(), CTPPSTimingLocalTrack::setPositionSigma(), CTPPSTimingLocalTrack::setTime(), CTPPSTimingLocalTrack::setTimeSigma(), CTPPSTimingLocalTrack::setValid(), CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::threshold_, CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::timeEval(), CTPPSTimingTrackRecognition< TotemTimingLocalTrack, TotemTimingRecHit >::tolerance_, HLT_FULL_cff::track, x, hit::x, y, and hit::y.

                                                                                       {
  int numberOfTracks = 0;
  DimensionParameters param;

  auto getX = [](const TotemTimingRecHit& hit) { return hit.x(); };
  auto getXWidth = [](const TotemTimingRecHit& hit) { return hit.xWidth(); };
  auto setX = [](TotemTimingLocalTrack& track, float x) { track.setPosition(math::XYZPoint(x, 0., 0.)); };
  auto setXSigma = [](TotemTimingLocalTrack& track, float sigma) {
    track.setPositionSigma(math::XYZPoint(sigma, 0., 0.));
  };
  auto getY = [](const TotemTimingRecHit& hit) { return hit.y(); };
  auto getYWidth = [](const TotemTimingRecHit& hit) { return hit.yWidth(); };
  auto setY = [](TotemTimingLocalTrack& track, float y) { track.setPosition(math::XYZPoint(0., y, 0.)); };
  auto setYSigma = [](TotemTimingLocalTrack& track, float sigma) {
    track.setPositionSigma(math::XYZPoint(0., sigma, 0.));
  };

  for (const auto& hitBatch : hitVectorMap_) {
    const auto& hits = hitBatch.second;
    const auto& hitRange = getHitSpatialRange(hits);

    TrackVector xPartTracks, yPartTracks;

    param.rangeBegin = hitRange.xBegin;
    param.rangeEnd = hitRange.xEnd;
    producePartialTracks(hits, param, getX, getXWidth, setX, setXSigma, xPartTracks);

    param.rangeBegin = hitRange.yBegin;
    param.rangeEnd = hitRange.yEnd;
    producePartialTracks(hits, param, getY, getYWidth, setY, setYSigma, yPartTracks);

    if (xPartTracks.empty() && yPartTracks.empty())
      continue;

    unsigned int validHitsNumber = (unsigned int)threshold_ + 1;

    for (const auto& xTrack : xPartTracks) {
      for (const auto& yTrack : yPartTracks) {
        math::XYZPoint position(xTrack.x0(), yTrack.y0(), 0.5f * (hitRange.zBegin + hitRange.zEnd));
        math::XYZPoint positionSigma(xTrack.x0Sigma(), yTrack.y0Sigma(), 0.5f * (hitRange.zEnd - hitRange.zBegin));

        TotemTimingLocalTrack newTrack(position, positionSigma, 0., 0.);

        HitVector componentHits;
        for (const auto& hit : hits)
          if (newTrack.containsHit(hit, tolerance_))
            componentHits.emplace_back(hit);
        if (componentHits.size() < validHitsNumber)
          continue;

        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);
        // in a next iteration, we will be setting validity / numHits / numPlanes
        tracks.push_back(newTrack);
      }
    }
  }

  return numberOfTracks;
}