CDCTrack Class Reference

Class representing a sequence of three dimensional reconstructed hits. More...

#include <CDCTrack.h>

Inheritance diagram for CDCTrack:

Public Member Functions
	CDCTrack ()=default
	Default constructor for ROOT compatibility.
	CDCTrack (const std::vector< CDCRecoHit3D > &recoHits3D)
	Constructor from a series of hits.
	CDCTrack (const CDCSegment2D &segment)
	Constructor from a two dimensional segment filling the third dimension with 0 values.
bool	operator< (const CDCTrack &track) const
	Comparison of track - no particular order has been defined so far, all tracks are equivalent.
std::vector< CDCSegment3D >	splitIntoSegments () const
	Splits the track into segments.
ISuperLayer	getStartISuperLayer () const
	Getter for the superlayer id the track starts from.
ISuperLayer	getEndISuperLayer () const
	Getter for the superlayer id the track ends in.
const Vector3D &	getStartRecoPos3D () const
	Getter for the position of the first reconstructed hit.
const Vector3D &	getEndRecoPos3D () const
	Getter for the position of the last reconstructed hit.
void	setStartTrajectory3D (const CDCTrajectory3D &startTrajectory3D)
	Setter for the two dimensional trajectory.
void	setEndTrajectory3D (const CDCTrajectory3D &endTrajectory3D)
	Setter for the three dimensional trajectory.
const CDCTrajectory3D &	getStartTrajectory3D () const
	Getter for the two dimensional trajectory.
const CDCTrajectory3D &	getEndTrajectory3D () const
	Getter for the three dimensional trajectory.
AutomatonCell &	getAutomatonCell () const
	Mutable getter for the automaton cell.
AutomatonCell *	operator-> () const
	Indirection to the automaton cell for easier access to the flags.
void	unsetAndForwardMaskedFlag () const
	Unset the masked flag of the automaton cell of this segment and of all contained wire hits.
void	setAndForwardMaskedFlag () const
	Set the masked flag of the automaton cell of this segment and forward the masked flag to all contained wire hits.
void	receiveMaskedFlag () const
	Check all contained wire hits if one has the masked flag.
void	forwardTakenFlag (bool takenFlag=true) const
	Set the taken flag of all hits belonging to this track to the given value (default true), but do not touch the flag of the track itself.
void	sortByArcLength2D ()
	Sort the recoHits according to their perpS information.
void	shiftToPositiveArcLengths2D (bool doForAllTracks=false)
	Set all arcLengths to have positive values by shifting them by pi*radius if they are negative.
void	reverse ()
	Reverse the track inplace.
CDCTrack	reversed () const
	Return a reversed copy of the track.
MayBePtr< const CDCRecoHit3D >	find (const CDCWireHit &wireHit) const
	Finds the first CDCRecoHit3D that is based on the given wire hit - nullptr if none.
void	setHasMatchingSegment (bool hasMatchingSegment=true)
	Set the flag which indicates that the track has a matching segment (probably only used in the SegmentTrackCombiner).
bool	getHasMatchingSegment () const
	Get a flag which indicates that the track has a matching segment (probably set in the SegmentTrackCombiner).
float	getQualityIndicator () const
	Get the multivariate quality indicator in [0, 1] attached the CDCTrack.
void	setQualityIndicator (const float qualityIndicator)
	Set the multivariate quality indicator in [0, 1] attached the CDCTrack.

Static Public Member Functions
static CDCTrack	condense (const Path< const CDCTrack > &trackPath)
	concatenates several tracks from a path
static CDCTrack	condense (const Path< const CDCSegmentTriple > &segmentTriplePath)
	Reconstructs the hit content of the segment triple track to a CDCTrack averaging overlapping parts.
static CDCTrack	condense (const Path< const CDCSegmentPair > &segmentPairPath)
	Reconstructs the hit content of the axial stereo segment pair path to a CDCTrack averaging overlapping parts.

Public Attributes
T	elements
	STL member.

Private Attributes
AutomatonCell	m_automatonCell
	Memory for the automaton cell.
CDCTrajectory3D	m_startTrajectory3D
	Memory for the three dimensional trajectory at the start of the track.
CDCTrajectory3D	m_endTrajectory3D
	Memory for the three dimensional trajectory at the end of the track.
bool	m_hasMatchingSegment = false
	Flag which indicates that the track had a matching segment (can be used for filter decisions)
float	m_qualityIndicator = NAN
	Multivariate classifier output in [0, 1] correlated with probability that track is good match.

Detailed Description

Class representing a sequence of three dimensional reconstructed hits.

Definition at line 41 of file CDCTrack.h.

Constructor & Destructor Documentation

CDCTrack() [1/2]

CDCTrack ( const std::vector< CDCRecoHit3D > &  recoHits3D )

explicit

Constructor from a series of hits.

Definition at line 133 of file CDCTrack.cc.

  : std::vector<CDCRecoHit3D>(recoHits3D)
{
}

CDCTrack() [2/2]

CDCTrack ( const CDCSegment2D &  segment )

explicit

Constructor from a two dimensional segment filling the third dimension with 0 values.

Definition at line 138 of file CDCTrack.cc.

                                              :
  m_startTrajectory3D(segment.getTrajectory2D()),
  m_endTrajectory3D(segment.getTrajectory2D())
{
  if (segment.empty()) return;
 
  // Adjust the start point
  const CDCRecoHit2D& startRecoHit2D = segment.front();
  const CDCRecoHit2D& endRecoHit2D = segment.back();
 
  Vector3D startPos3D(startRecoHit2D.getRecoPos2D(), 0.0);
  Vector3D endPos3D(endRecoHit2D.getRecoPos2D(), 0.0);
 
  m_startTrajectory3D.setLocalOrigin(startPos3D);
  m_endTrajectory3D.setLocalOrigin(endPos3D);
 
  for (const CDCRecoHit2D& recoHit2D : segment) {
    const CDCRLWireHit& rlWireHit = recoHit2D.getRLWireHit();
    Vector3D recoPos3D(recoHit2D.getRecoPos2D(), 0.0);
    double perpS = m_startTrajectory3D.calcArcLength2D(recoPos3D);
    push_back(CDCRecoHit3D(rlWireHit, recoPos3D, perpS));
  }
 
  // TODO: Maybe enhance the estimation of the z coordinate with the superlayer slopes.
}

Member Function Documentation

condense() [1/3]

CDCTrack condense ( const Path< const CDCSegmentPair > &  segmentPairPath )

static

Reconstructs the hit content of the axial stereo segment pair path to a CDCTrack averaging overlapping parts.

Definition at line 263 of file CDCTrack.cc.

{
  CDCTrack track;
 
  //B2DEBUG(200,"Length of segmentTripleTrack is " << segmentTripleTrack.size() );
  if (segmentPairPath.empty()) return track;
 
  Path<const CDCSegmentPair>::const_iterator itSegmentPair = segmentPairPath.begin();
  const CDCSegmentPair* firstSegmentPair = *itSegmentPair++;
 
 
  // Keep the fit of the first segment pair to set it as the fit at the start of the track
  CDCTrajectory3D startTrajectory3D = firstSegmentPair->getTrajectory3D();
 
  double perpSOffset = 0.0;
  appendReconstructed(firstSegmentPair->getFromSegment(),
                      firstSegmentPair->getTrajectory3D(),
                      perpSOffset, track);
 
  while (itSegmentPair != segmentPairPath.end()) {
 
    const CDCSegmentPair* secondSegmentPair = *itSegmentPair++;
 
    B2ASSERT("Two segment pairs do not overlap in their segments",
             firstSegmentPair->getToSegment() == secondSegmentPair->getFromSegment());
 
    perpSOffset = appendReconstructedAverage(firstSegmentPair->getToSegment(),
                                             firstSegmentPair->getTrajectory3D(),
                                             perpSOffset,
                                             secondSegmentPair->getTrajectory3D(),
                                             track);
 
    firstSegmentPair = secondSegmentPair;
  }
 
  const CDCSegmentPair* lastSegmentPair = firstSegmentPair;
  appendReconstructed(lastSegmentPair->getToSegment(),
                      lastSegmentPair->getTrajectory3D(),
                      perpSOffset, track);
 
  // Keep the fit of the last segment pair to set it as the fit at the end of the track
  CDCTrajectory3D endTrajectory3D = lastSegmentPair->getTrajectory3D();
 
  // Move the reference point of the start fit to the first observed position
  double resetPerpSOffset = startTrajectory3D.setLocalOrigin(track.front().getRecoPos3D());
  track.setStartTrajectory3D(startTrajectory3D);
 
  // Move the reference point of the end fit to the last observed position
  endTrajectory3D.setLocalOrigin(track.back().getRecoPos3D());
  track.setEndTrajectory3D(endTrajectory3D);
 
  for (CDCRecoHit3D& recoHit3D : track) {
    recoHit3D.shiftArcLength2D(-resetPerpSOffset);
  }
 
  return track;
}

condense() [2/3]

CDCTrack condense ( const Path< const CDCSegmentTriple > &  segmentTriplePath )

static

Reconstructs the hit content of the segment triple track to a CDCTrack averaging overlapping parts.

Definition at line 197 of file CDCTrack.cc.

{
  CDCTrack track;
  // B2DEBUG(200,"Length of segmentTripleTrack is " << segmentTripleTrack.size() );
  if (segmentTriplePath.empty()) return track;
 
  Path<const CDCSegmentTriple>::const_iterator itSegmentTriple = segmentTriplePath.begin();
  const CDCSegmentTriple* firstSegmentTriple = *itSegmentTriple++;
 
  // Set the start fits of the track to the ones of the first segment
  CDCTrajectory3D startTrajectory3D = firstSegmentTriple->getTrajectory3D();
 
 
  double perpSOffset = 0.0;
  appendReconstructed(firstSegmentTriple->getStartSegment(),
                      firstSegmentTriple->getTrajectory3D(),
                      perpSOffset,
                      track);
 
  appendReconstructed(firstSegmentTriple->getMiddleSegment(),
                      firstSegmentTriple->getTrajectory3D(),
                      perpSOffset, track);
 
  while (itSegmentTriple != segmentTriplePath.end()) {
 
    const CDCSegmentTriple* secondSegmentTriple = *itSegmentTriple++;
    B2ASSERT("Two segment triples do not overlap in their axial segments",
             firstSegmentTriple->getEndSegment() == secondSegmentTriple->getStartSegment());
 
    perpSOffset = appendReconstructedAverage(firstSegmentTriple->getEndSegment(),
                                             firstSegmentTriple->getTrajectory3D(),
                                             perpSOffset,
                                             secondSegmentTriple->getTrajectory3D(),
                                             track);
 
    appendReconstructed(secondSegmentTriple->getMiddleSegment(),
                        secondSegmentTriple->getTrajectory3D(),
                        perpSOffset, track);
 
    firstSegmentTriple = secondSegmentTriple;
 
  }
 
  const CDCSegmentTriple* lastSegmentTriple = firstSegmentTriple;
 
  appendReconstructed(lastSegmentTriple->getEndSegment(),
                      lastSegmentTriple->getTrajectory3D(),
                      perpSOffset, track);
 
  // Set the end fits of the track to the ones of the last segment
  CDCTrajectory3D endTrajectory3D = lastSegmentTriple->getTrajectory3D();
 
  // Set the reference point on the trajectories to the last reconstructed hit
  double resetPerpSOffset = startTrajectory3D.setLocalOrigin(track.front().getRecoPos3D());
  track.setStartTrajectory3D(startTrajectory3D);
 
  endTrajectory3D.setLocalOrigin(track.back().getRecoPos3D());
  track.setEndTrajectory3D(endTrajectory3D);
 
  for (CDCRecoHit3D& recoHit3D : track) {
    recoHit3D.shiftArcLength2D(-resetPerpSOffset);
  }
 
  return track;
}

condense() [3/3]

CDCTrack condense ( const Path< const CDCTrack > &  trackPath )

static

concatenates several tracks from a path

FIXME : arc lengths are not set properly

Definition at line 164 of file CDCTrack.cc.

{
  if (trackPath.empty()) {
    return CDCTrack();
  } else if (trackPath.size() == 1) {
    return CDCTrack(*(trackPath[0]));
  } else {
    CDCTrack result;
    for (const CDCTrack* track :  trackPath) {
      for (const CDCRecoHit3D& recoHit3D : *track) {
        result.push_back(recoHit3D);
      }
    }
 
    CDCTrajectory3D startTrajectory3D = trackPath.front()->getStartTrajectory3D();
    CDCTrajectory3D endTrajectory3D = trackPath.back()->getStartTrajectory3D();
 
    double resetPerpSOffset =
      startTrajectory3D.setLocalOrigin(result.front().getRecoPos3D());
    result.setStartTrajectory3D(startTrajectory3D);
 
    endTrajectory3D.setLocalOrigin(result.back().getRecoPos3D());
    result.setEndTrajectory3D(endTrajectory3D);
 
    for (CDCRecoHit3D& recoHit3D : result) {
      recoHit3D.shiftArcLength2D(-resetPerpSOffset);
    }
 
    return result;
  }
}

find()

MayBePtr< const CDCRecoHit3D > find

(

const CDCWireHit &

wireHit

)

const

Finds the first CDCRecoHit3D that is based on the given wire hit - nullptr if none.

Definition at line 367 of file CDCTrack.cc.

{
  auto hasWireHit = [&wireHit](const CDCRecoHit3D & recoHit3D) {
    return recoHit3D.hasWireHit(wireHit);
  };
  auto itRecoHit3D = std::find_if(this->begin(), this->end(), hasWireHit);
  return itRecoHit3D == this->end() ? nullptr : &*itRecoHit3D;
}

forwardTakenFlag()

void forwardTakenFlag

(

bool

takenFlag = true

)

const

Set the taken flag of all hits belonging to this track to the given value (default true), but do not touch the flag of the track itself.

Definition at line 405 of file CDCTrack.cc.

{
  for (const CDCRecoHit3D& recoHit3D : *this) {
    recoHit3D.getWireHit().getAutomatonCell().setTakenFlag(takenFlag);
  }
}

getAutomatonCell()

AutomatonCell & getAutomatonCell ( ) const

inline

Mutable getter for the automaton cell.

Definition at line 125 of file CDCTrack.h.

      {
        return m_automatonCell;
      }

getEndISuperLayer()

ISuperLayer getEndISuperLayer ( ) const

inline

Getter for the superlayer id the track ends in.

Definition at line 79 of file CDCTrack.h.

      {
        return back().getISuperLayer();
      }

getEndRecoPos3D()

const Vector3D & getEndRecoPos3D ( ) const

inline

Getter for the position of the last reconstructed hit.

Definition at line 91 of file CDCTrack.h.

      {
        return back().getRecoPos3D();
      }

getEndTrajectory3D()

const CDCTrajectory3D & getEndTrajectory3D ( ) const

inline

Getter for the three dimensional trajectory.

The trajectory should start at the END of the track and follow its direction.

Definition at line 119 of file CDCTrack.h.

      {
        return m_endTrajectory3D;
      }

getHasMatchingSegment()

bool getHasMatchingSegment ( ) const

inline

Get a flag which indicates that the track has a matching segment (probably set in the SegmentTrackCombiner).

This flag can be used for filter decisions (e.g. if the track is fake).

Definition at line 182 of file CDCTrack.h.

      {
        return m_hasMatchingSegment;
      }

getQualityIndicator()

float getQualityIndicator ( ) const

inline

Get the multivariate quality indicator in [0, 1] attached the CDCTrack.

Definition at line 188 of file CDCTrack.h.

      {
        return m_qualityIndicator;
      }

getStartISuperLayer()

ISuperLayer getStartISuperLayer ( ) const

inline

Getter for the superlayer id the track starts from.

Definition at line 73 of file CDCTrack.h.

      {
        return front().getISuperLayer();
      }

getStartRecoPos3D()

const Vector3D & getStartRecoPos3D ( ) const

inline

Getter for the position of the first reconstructed hit.

Definition at line 85 of file CDCTrack.h.

      {
        return front().getRecoPos3D();
      }

getStartTrajectory3D()

const CDCTrajectory3D & getStartTrajectory3D ( ) const

inline

Getter for the two dimensional trajectory.

The trajectory should start at the start of the track and follow its direction.

Definition at line 112 of file CDCTrack.h.

      {
        return m_startTrajectory3D;
      }

operator->()

AutomatonCell * operator-> ( ) const

inline

Indirection to the automaton cell for easier access to the flags.

Definition at line 131 of file CDCTrack.h.

      {
        return &m_automatonCell;
      }

operator<()

bool operator< ( const CDCTrack &  track ) const

inline

Comparison of track - no particular order has been defined so far, all tracks are equivalent.

Definition at line 63 of file CDCTrack.h.

      {
        return false;
      }

receiveMaskedFlag()

void receiveMaskedFlag

(

)

const

Check all contained wire hits if one has the masked flag.

Set the masked flag of this segment in case at least one of the contained wire hits is flagged as masked.

Definition at line 394 of file CDCTrack.cc.

{
  for (const CDCRecoHit3D& recoHit3D : *this) {
    const CDCWireHit& wireHit = recoHit3D.getWireHit();
    if (wireHit.getAutomatonCell().hasMaskedFlag()) {
      getAutomatonCell().setMaskedFlag();
      return;
    }
  }
}

reverse()

void reverse

(

)

Reverse the track inplace.

Definition at line 336 of file CDCTrack.cc.

{
  if (empty()) return;
 
  // Exchange the forward and backward trajectory and reverse them
  std::swap(m_startTrajectory3D, m_endTrajectory3D);
  m_startTrajectory3D.reverse();
  m_endTrajectory3D.reverse();
 
  const CDCRecoHit3D& lastRecoHit3D = back();
  double lastPerpS = lastRecoHit3D.getArcLength2D();
  double newLastPerpS = m_startTrajectory3D.calcArcLength2D(lastRecoHit3D.getRecoPos3D());
 
  // Reverse the left right passage hypotheses and reverse the measured travel distance
  for (CDCRecoHit3D& recoHit3D : *this) {
    recoHit3D.reverse();
    double perpS = recoHit3D.getArcLength2D();
    recoHit3D.setArcLength2D(newLastPerpS + lastPerpS - perpS);
  }
 
  // Reverse the arrangement of hits.
  std::reverse(begin(), end());
}

reversed()

CDCTrack reversed

(

)

const

Return a reversed copy of the track.

Definition at line 360 of file CDCTrack.cc.

{
  CDCTrack reversedTrack(*this);
  reversedTrack.reverse();
  return reversedTrack;
}

setAndForwardMaskedFlag()

void setAndForwardMaskedFlag

(

)

const

Set the masked flag of the automaton cell of this segment and forward the masked flag to all contained wire hits.

Definition at line 385 of file CDCTrack.cc.

{
  getAutomatonCell().setMaskedFlag();
  for (const CDCRecoHit3D& recoHit3D : *this) {
    const CDCWireHit& wireHit = recoHit3D.getWireHit();
    wireHit.getAutomatonCell().setMaskedFlag();
  }
}

setEndTrajectory3D()

void setEndTrajectory3D ( const CDCTrajectory3D &  endTrajectory3D )

inline

Setter for the three dimensional trajectory.

The trajectory should start at the END of the track and follow its direction.

Definition at line 105 of file CDCTrack.h.

      {
        m_endTrajectory3D = endTrajectory3D;
      }

setHasMatchingSegment()

void setHasMatchingSegment ( bool  hasMatchingSegment = true )

inline

Set the flag which indicates that the track has a matching segment (probably only used in the SegmentTrackCombiner).

Definition at line 173 of file CDCTrack.h.

      {
        m_hasMatchingSegment = hasMatchingSegment;
      }

setQualityIndicator()

void setQualityIndicator ( const float  qualityIndicator )

inline

Set the multivariate quality indicator in [0, 1] attached the CDCTrack.

Definition at line 194 of file CDCTrack.h.

      {
        m_qualityIndicator = qualityIndicator;
      }

setStartTrajectory3D()

void setStartTrajectory3D ( const CDCTrajectory3D &  startTrajectory3D )

inline

Setter for the two dimensional trajectory.

The trajectory should start at the start of the track and follow its direction.

Definition at line 98 of file CDCTrack.h.

      {
        m_startTrajectory3D = startTrajectory3D;
      }

shiftToPositiveArcLengths2D()

void shiftToPositiveArcLengths2D

(

bool

doForAllTracks = false

)

Set all arcLengths to have positive values by shifting them by pi*radius if they are negative.

This can only be done if the radius is not infinity (for example cosmics). The flag can be used to do this for all tracks (default is to do this only for curlers)

Definition at line 423 of file CDCTrack.cc.

{
  const CDCTrajectory2D& startTrajectory2D = getStartTrajectory3D().getTrajectory2D();
  if (doForAllTracks or startTrajectory2D.isCurler(1.1)) {
    const double shiftValue = startTrajectory2D.getLocalCircle()->arcLengthPeriod();
    if (std::isfinite(shiftValue)) {
      for (CDCRecoHit3D& recoHit : *this) {
        if (recoHit.getArcLength2D() < 0) {
          recoHit.shiftArcLength2D(shiftValue);
        }
      }
    }
  }
}

sortByArcLength2D()

void sortByArcLength2D

(

)

Sort the recoHits according to their perpS information.

Definition at line 412 of file CDCTrack.cc.

{
  std::stable_sort(begin(),
                   end(),
  [](const CDCRecoHit3D & recoHit, const CDCRecoHit3D & otherRecoHit) {
    double arcLength = recoHit.getArcLength2D();
    double otherArcLength = otherRecoHit.getArcLength2D();
    return lessFloatHighNaN(arcLength, otherArcLength);
  });
}

splitIntoSegments()

std::vector< CDCSegment3D > splitIntoSegments

(

)

const

Splits the track into segments.

Note : No trajectory information is copied

Definition at line 321 of file CDCTrack.cc.

{
  vector<CDCSegment3D> result;
  ISuperLayer lastISuperLayer = -1;
  for (const CDCRecoHit3D& recoHit3D : *this) {
    ISuperLayer iSuperLayer = recoHit3D.getISuperLayer();
    if (result.empty() or lastISuperLayer != iSuperLayer) {
      result.emplace_back();
    }
    result.back().push_back(recoHit3D);
    lastISuperLayer = iSuperLayer;
  }
  return result;
}

unsetAndForwardMaskedFlag()

void unsetAndForwardMaskedFlag

(

)

const

Unset the masked flag of the automaton cell of this segment and of all contained wire hits.

Definition at line 376 of file CDCTrack.cc.

{
  getAutomatonCell().unsetMaskedFlag();
  for (const CDCRecoHit3D& recoHit3D : *this) {
    const CDCWireHit& wireHit = recoHit3D.getWireHit();
    wireHit.getAutomatonCell().unsetMaskedFlag();
  }
}

Member Data Documentation

m_automatonCell

AutomatonCell m_automatonCell

mutableprivate

Memory for the automaton cell.

Definition at line 201 of file CDCTrack.h.

m_endTrajectory3D

CDCTrajectory3D m_endTrajectory3D

private

Memory for the three dimensional trajectory at the end of the track.

Definition at line 207 of file CDCTrack.h.

m_hasMatchingSegment

bool m_hasMatchingSegment = false

private

Flag which indicates that the track had a matching segment (can be used for filter decisions)

Definition at line 210 of file CDCTrack.h.

m_qualityIndicator

float m_qualityIndicator = NAN

private

Multivariate classifier output in [0, 1] correlated with probability that track is good match.

Used to reject fakes and clones. Implemented here to be forwarded to RecoTrack.

Definition at line 214 of file CDCTrack.h.

m_startTrajectory3D

CDCTrajectory3D m_startTrajectory3D

private

Memory for the three dimensional trajectory at the start of the track.

Definition at line 204 of file CDCTrack.h.

---

The documentation for this class was generated from the following files:

• tracking/trackFindingCDC/eventdata/tracks/include/CDCTrack.h
• tracking/trackFindingCDC/eventdata/tracks/src/CDCTrack.cc