Class implementing the Riemann fit for two dimensional trajectory circle. More...

#include <CDCRiemannFitter.h>

Inheritance diagram for CDCRiemannFitter:

Public Member Functions
	CDCRiemannFitter ()
	Default constructor.

	~CDCRiemannFitter ()
	Default destructor.

CDCTrajectory2D	fit (const CDCObservations2D &observations2D) const
	Fits a collection of observation drift circles.

CDCTrajectory2D	fit (CDCObservations2D &&observations2D) const
	Fits a collection of observation drift circles.

CDCTrajectory2D	fit (const CDCTrack &track) const
	Fits the track.

CDCTrajectory2D	fit (const CDCSegment3D &segment) const
	Fits the segment.

CDCTrajectory2D	fit (const CDCSegment2D &fromSegment, const CDCSegment2D &toSegment) const
	Fits to the two segments.

CDCTrajectory2D	fit (const CDCSegment2D &segment) const
	Fits the segment.

CDCTrajectory2D	fit (const std::vector< const CDCWireHit * > &wireHit) const
	Fits to the wire hit observations.

CDCTrajectory2D	fit (const std::vector< const CDCWire * > &wires) const
	Fits to the wire positions. Explicit specialisation to be used from python.

CDCTrajectory2D	fit (const CDCWireHitSegment &wireHits) const
	Fits to the wire positions. Explicit specialisation to be used from python.

void	update (CDCTrajectory2D &trajectory2D, const CDCObservations2D &observations2D) const
	Update the trajectory with a fit to the observations.

void	update (CDCTrajectory2D &trajectory2D, CDCObservations2D &&observations2D) const
	Update the trajectory with a fit to the observations.

void	update (CDCTrajectory2D &trajectory2D, const CDCSegment2D &segment) const
	Updates the given trajectory inplace from the given segment.

void	update (CDCTrajectory2D &trajectory2D, const CDCAxialSegmentPair &axialSegmentPair) const
	Updates the given trajectory inplace from the given segment pair.

void	update (CDCTrajectory2D &trajectory2D, CDCObservations2D &observations2D) const
	Executes the fit and updates the trajectory parameters This may render the information in the observation object.

void	useOnlyPosition ()
	Setup the fitter to use only the reconstructed positions of the hits.

void	useOnlyOrientation ()
	Setup the fitter to use only reference position and the drift length with right left orientation.

void	usePositionAndOrientation ()
	Setup the fitter to use both the reconstructed position and the reference position and the drift length with right left orientation.

void	setFitVariance (EFitVariance fitVariance)
	Setup the fitter to use the given variance measure by default.

bool	isLineConstrained () const
	Getter for the indicator that lines should be fitted by this fitter.

bool	isOriginConstrained () const
	Getter for the indicator that curves through the origin should be fitted by this fitter.

void	setLineConstrained (bool constrained=true)
	Indicator if this fitter is setup to fit lines.

void	setOriginConstrained (bool constrained=true)
	Indicator if this fitter is setup to fit curves through the origin.

Static Public Member Functions
static const CDCRiemannFitter &	getFitter ()
	Static getter for a general Riemann fitter.

static const CDCRiemannFitter &	getFitter (bool fromOrigin, bool line)
	Static getter for a given Riemann fitter.

static const CDCRiemannFitter &	getLineFitter ()
	Static getter for a line fitter.

static const CDCRiemannFitter &	getOriginCircleFitter ()
	Static getter for an origin circle fitter.

Private Member Functions
CDCTrajectory2D	fitGeneric (const AHits &hits) const
	Fits a collection of hit typs which are convertible to observation circles.

CDCTrajectory2D	fitGeneric (const AStartHits &startHits, const AEndHits &endHits) const
	Fits together two collections of hit types which are convertible to observation circles.

void	updateGeneric (CDCTrajectory2D &trajectory2D, const AHits &hits) const
	Updates a given trajectory with a fit to a collection of hits types, which are convertible to observation circles.

void	updateGeneric (CDCTrajectory2D &trajectory2D, const AStartHits &startHits, const AEndHits &endHits) const
	Updates a given trajectory with a fit to two collection of hit types, which are convertible to observation circles.

UncertainPerigeeCircle	fitInternal (CDCObservations2D &observations2D) const
	Internal method doing the heavy work.

Private Attributes
bool	m_usePosition
	Flag indicating the reconstructed position shall be used in the fit.

bool	m_useOrientation
	Flag indicating the reference position and drift length with right left orientation shall be used in the fit.

EFitVariance	m_fitVariance
	Default variance to be used in the fit.

bool	m_lineConstrained
	Memory for the flag indicating that lines should be fitted.

bool	m_originConstrained
	Memory for the flag indicating that curves through the origin shall be fitter.

Detailed Description

Class implementing the Riemann fit for two dimensional trajectory circle.

Definition at line 24 of file CDCRiemannFitter.h.

Member Function Documentation

◆ fit() [1/9]

CDCTrajectory2D fit ( CDCObservations2D && observations2D ) const

inherited

Fits a collection of observation drift circles.

Definition at line 46 of file CDCFitter2D.icc.h.

    {
      CDCTrajectory2D result;
      update(result, observations2D);
      return result;
    }

◆ fit() [2/9]

CDCTrajectory2D fit ( const CDCObservations2D & observations2D ) const

inherited

Fits a collection of observation drift circles.

Definition at line 43 of file CDCFitter2D.icc.h.

    {
      return fit(CDCObservations2D(observations2D));
    }

◆ fit() [3/9]

CDCTrajectory2D fit	(	const CDCSegment2D &	fromSegment,
		const CDCSegment2D &	toSegment
	)		const

inherited

Fits to the two segments.

Definition at line 61 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(fromSegment, toSegment);
    }

◆ fit() [4/9]

CDCTrajectory2D fit ( const CDCSegment2D & segment ) const

inherited

Fits the segment.

Definition at line 64 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(segment);
    }

◆ fit() [5/9]

CDCTrajectory2D fit ( const CDCSegment3D & segment ) const

inherited

Fits the segment.

Definition at line 58 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(segment);
    }

◆ fit() [6/9]

CDCTrajectory2D fit ( const CDCTrack & track ) const

inherited

Fits the track.

Definition at line 55 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(track);
    }

◆ fit() [7/9]

CDCTrajectory2D fit ( const CDCWireHitSegment & wireHits ) const

inherited

Fits to the wire positions. Explicit specialisation to be used from python.

Definition at line 73 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(wireHits);
    }

◆ fit() [8/9]

CDCTrajectory2D fit ( const std::vector< const CDCWire * > & wires ) const

inherited

Fits to the wire positions. Explicit specialisation to be used from python.

Definition at line 70 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(wires);
    }

◆ fit() [9/9]

CDCTrajectory2D fit ( const std::vector< const CDCWireHit * > & wireHit ) const

inherited

Fits to the wire hit observations.

Definition at line 67 of file CDCFitter2D.icc.h.

    {
      return fitGeneric(wireHits);
    }

◆ fitGeneric() [1/2]

CDCTrajectory2D fitGeneric ( const AHits & hits ) const

privateinherited

Fits a collection of hit typs which are convertible to observation circles.

Definition at line 84 of file CDCFitter2D.icc.h.

    {
      CDCTrajectory2D result;
      updateGeneric(result, hits);
      return result;
    }

◆ fitGeneric() [2/2]

CDCTrajectory2D fitGeneric	(	const AStartHits &	startHits,
		const AEndHits &	endHits
	)		const

privateinherited

Fits together two collections of hit types which are convertible to observation circles.

Definition at line 88 of file CDCFitter2D.icc.h.

    {
      CDCTrajectory2D result;
      updateGeneric(result, startHits, endHits);
      return result;
    }

◆ fitInternal()

UncertainPerigeeCircle fitInternal ( CDCObservations2D & observations2D ) const

privateinherited

Internal method doing the heavy work.

Definition at line 294 of file ExtendedRiemannsMethod.cc.

{
  using namespace NParabolicParameterIndices;
 
  // Matrix of weighted sums
  Eigen::Matrix< double, 5, 5 > s = getWXYRLSumMatrix(observations2D);
 
  // The same as above without drift lengths
  Eigen::Matrix<double, 4, 4> sNoL = s.block<4, 4>(0, 0);
 
  // Determine NDF : Circle fit eats up to 3 degrees of freedom depending on the constraints
  size_t ndf = observations2D.size() - 1;
 
  if (not isOriginConstrained()) {
    --ndf;
  }
 
  if (not isLineConstrained()) {
    --ndf;
  }
 
  // Parameters to be fitted
  UncertainPerigeeCircle resultCircle;
  double chi2 = 0;
 
  size_t nObservationsWithDriftRadius = observations2D.getNObservationsWithDriftRadius();
  if (nObservationsWithDriftRadius > 0) {
    resultCircle = UncertainPerigeeCircle(::fit(s, isLineConstrained(), isOriginConstrained()));
    chi2 = calcChi2(resultCircle, s);
  } else {
    if (not isOriginConstrained()) {
      // Alternative implementation for comparison
 
      // Matrix of averages
      Eigen::Matrix< double, 4, 4> aNoL = sNoL / sNoL(iW);
 
      // Measurement means
      Eigen::Matrix< double, 4, 1> meansNoL = aNoL.row(iW);
 
      // Covariance matrix
      Eigen::Matrix< double, 4, 4> cNoL = aNoL - meansNoL * meansNoL.transpose();
 
      resultCircle = UncertainPerigeeCircle(fitSeperateOffset(meansNoL, cNoL, isLineConstrained()));
 
    } else {
      resultCircle = UncertainPerigeeCircle(::fit(sNoL, isLineConstrained(), isOriginConstrained()));
    }
 
    chi2 = calcChi2(resultCircle, sNoL);
  }
 
  // Covariance calculation does not need the drift lengths, which is why we do not forward them.
  PerigeePrecision perigeePrecision =
    calcPrecision(resultCircle, sNoL, isLineConstrained(), isOriginConstrained());
 
  // Use in pivoting in case the matrix is not full rank as it is for the constrained cases
  PerigeeCovariance perigeeCovariance = PerigeeUtil::covarianceFromPrecision(perigeePrecision);
 
  resultCircle.setNDF(ndf);
  resultCircle.setChi2(chi2);
  resultCircle.setPerigeeCovariance(perigeeCovariance);
  return resultCircle;
}

◆ getFitter() [1/2]

const CDCRiemannFitter & getFitter ( )

static

Static getter for a general Riemann fitter.

Definition at line 20 of file CDCRiemannFitter.cc.

{
  static CDCRiemannFitter fitter;
  return fitter;
}

◆ getFitter() [2/2]

const CDCRiemannFitter & getFitter	(	bool	fromOrigin,
		bool	line
	)

static

Static getter for a given Riemann fitter.

Definition at line 26 of file CDCRiemannFitter.cc.

{
  static CDCRiemannFitter fitter;
  if (fromOrigin) fitter.setOriginConstrained();
  if (line) fitter.setLineConstrained();
  return fitter;
}

◆ getLineFitter()

const CDCRiemannFitter & getLineFitter ( )

static

Static getter for a line fitter.

Definition at line 34 of file CDCRiemannFitter.cc.

{
  static CDCRiemannFitter lineFitter;
  lineFitter.setLineConstrained();
  return lineFitter;
}

◆ getOriginCircleFitter()

const CDCRiemannFitter & getOriginCircleFitter ( )

static

Static getter for an origin circle fitter.

Definition at line 41 of file CDCRiemannFitter.cc.

{
  static CDCRiemannFitter originCircleFitter;
  originCircleFitter.setOriginConstrained();
  return originCircleFitter;
}

◆ isLineConstrained()

bool isLineConstrained ( ) const

inlineinherited

Getter for the indicator that lines should be fitted by this fitter.

Definition at line 39 of file ExtendedRiemannsMethod.h.

      {
        return m_lineConstrained;
      }

◆ isOriginConstrained()

bool isOriginConstrained ( ) const

inlineinherited

Getter for the indicator that curves through the origin should be fitted by this fitter.

Definition at line 45 of file ExtendedRiemannsMethod.h.

      {
        return m_originConstrained;
      }

◆ setFitVariance()

void setFitVariance ( EFitVariance fitVariance )

inherited

Setup the fitter to use the given variance measure by default.

Definition at line 120 of file CDCFitter2D.icc.h.

    {
      m_fitVariance = fitVariance;
    }

◆ setLineConstrained()

void setLineConstrained ( bool constrained = true )

inlineinherited

Indicator if this fitter is setup to fit lines.

Definition at line 51 of file ExtendedRiemannsMethod.h.

      {
        m_lineConstrained = constrained;
      }

◆ setOriginConstrained()

void setOriginConstrained ( bool constrained = true )

inlineinherited

Indicator if this fitter is setup to fit curves through the origin.

Definition at line 57 of file ExtendedRiemannsMethod.h.

      {
        m_originConstrained = constrained;
      }

◆ update() [1/5]

void update	(	CDCTrajectory2D &	trajectory2D,
		CDCObservations2D &&	observations2D
	)		const

inherited

Update the trajectory with a fit to the observations.

Definition at line 52 of file CDCFitter2D.icc.h.

    {
      AFitMethod::update(trajectory2D, observations2D);
    }

◆ update() [2/5]

void update	(	CDCTrajectory2D &	trajectory2D,
		CDCObservations2D &	observations2D
	)		const

inherited

Executes the fit and updates the trajectory parameters This may render the information in the observation object.

Definition at line 34 of file ExtendedRiemannsMethod.cc.

{
  size_t nObservations = observations2D.size();
  trajectory2D.clear();
  if (not nObservations) return;
 
  Vector2D origin = Vector2D(0.0, 0.0);
  Vector2D centralPoint = observations2D.getCentralPoint();
 
  const Vector2D& ref = isOriginConstrained() ? origin : centralPoint;
  observations2D.passiveMoveBy(ref);
 
  UncertainPerigeeCircle perigeeCircle = fitInternal(observations2D);
 
  double frontX = observations2D.getX(0);
  double frontY = observations2D.getY(0);
  Vector2D frontPos(frontX, frontY);
 
  double backX = observations2D.getX(nObservations - 1);
  double backY = observations2D.getY(nObservations - 1);
  Vector2D backPos(backX, backY);
 
  Vector2D overPos(0, 0);
  double totalPerps = (perigeeCircle->arcLengthBetween(frontPos, overPos) +
                       perigeeCircle->arcLengthBetween(overPos, backPos));
 
  if (totalPerps < 0) {
    perigeeCircle.reverse();
  }
 
  trajectory2D = CDCTrajectory2D(ref, perigeeCircle);
}

◆ update() [3/5]

void update	(	CDCTrajectory2D &	trajectory2D,
		const CDCAxialSegmentPair &	axialSegmentPair
	)		const

inherited

Updates the given trajectory inplace from the given segment pair.

Definition at line 79 of file CDCFitter2D.icc.h.

    {
      return updateGeneric(trajectory2D, axialSegmentPair);
    }

◆ update() [4/5]

void update	(	CDCTrajectory2D &	trajectory2D,
		const CDCObservations2D &	observations2D
	)		const

inherited

Update the trajectory with a fit to the observations.

Definition at line 49 of file CDCFitter2D.icc.h.

    {
      return update(trajectory2D, CDCObservations2D(observations2D));
    }

◆ update() [5/5]

void update	(	CDCTrajectory2D &	trajectory2D,
		const CDCSegment2D &	segment
	)		const

inherited

Updates the given trajectory inplace from the given segment.

Definition at line 76 of file CDCFitter2D.icc.h.

    {
      updateGeneric(trajectory2D, segment);
    }

◆ updateGeneric() [1/2]

void updateGeneric	(	CDCTrajectory2D &	trajectory2D,
		const AHits &	hits
	)		const

privateinherited

Updates a given trajectory with a fit to a collection of hits types, which are convertible to observation circles.

Definition at line 95 of file CDCFitter2D.icc.h.

    {
      CDCObservations2D observations2D;
      observations2D.setFitVariance(m_fitVariance);
 
      if (m_usePosition) {
        observations2D.setFitPos(EFitPos::c_RecoPos);
        observations2D.appendRange(hits);
      }
      if (m_useOrientation) {
        observations2D.setFitPos(EFitPos::c_RLDriftCircle);
        observations2D.appendRange(hits);
      }
 
      if (observations2D.size() < 4) {
        trajectory2D.clear();
      } else {
        AFitMethod::update(trajectory2D, observations2D);
      }
    }

◆ updateGeneric() [2/2]

void updateGeneric	(	CDCTrajectory2D &	trajectory2D,
		const AStartHits &	startHits,
		const AEndHits &	endHits
	)		const

privateinherited

Updates a given trajectory with a fit to two collection of hit types, which are convertible to observation circles.

Definition at line 102 of file CDCFitter2D.icc.h.

    {
      CDCObservations2D observations2D;
      observations2D.setFitVariance(m_fitVariance);
 
      if (m_usePosition) {
        observations2D.setFitPos(EFitPos::c_RecoPos);
        observations2D.appendRange(startHits);
      }
      if (m_useOrientation) {
        observations2D.setFitPos(EFitPos::c_RLDriftCircle);
        observations2D.appendRange(startHits);
      }
 
      if (m_usePosition) {
        observations2D.setFitPos(EFitPos::c_RecoPos);
        observations2D.appendRange(endHits);
      }
      if (m_useOrientation) {
        observations2D.setFitPos(EFitPos::c_RLDriftCircle);
        observations2D.appendRange(endHits);
      }
 
      if (observations2D.size() < 4) {
        trajectory2D.clear();
      } else {
        AFitMethod::update(trajectory2D, observations2D);
      }
    }

◆ useOnlyOrientation()

void useOnlyOrientation

inherited

Setup the fitter to use only reference position and the drift length with right left orientation.

Definition at line 114 of file CDCFitter2D.icc.h.

    {
      m_usePosition = false;
      m_useOrientation = true;
    }

◆ useOnlyPosition()

void useOnlyPosition

inherited

Setup the fitter to use only the reconstructed positions of the hits.

Definition at line 111 of file CDCFitter2D.icc.h.

    {
      m_usePosition = true;
      m_useOrientation = false;
    }

◆ usePositionAndOrientation()

void usePositionAndOrientation

inherited

Setup the fitter to use both the reconstructed position and the reference position and the drift length with right left orientation.

Definition at line 117 of file CDCFitter2D.icc.h.

    {
      m_usePosition = true;
      m_useOrientation = true;
    }

Member Data Documentation

◆ m_fitVariance

EFitVariance m_fitVariance

privateinherited

Default variance to be used in the fit.

Definition at line 130 of file CDCFitter2D.h.

◆ m_lineConstrained

bool m_lineConstrained

privateinherited

Memory for the flag indicating that lines should be fitted.

Definition at line 64 of file ExtendedRiemannsMethod.h.

◆ m_originConstrained

bool m_originConstrained

privateinherited

Memory for the flag indicating that curves through the origin shall be fitter.

Definition at line 67 of file ExtendedRiemannsMethod.h.

◆ m_useOrientation

bool m_useOrientation

privateinherited

Flag indicating the reference position and drift length with right left orientation shall be used in the fit.

Definition at line 127 of file CDCFitter2D.h.

◆ m_usePosition

bool m_usePosition

privateinherited

Flag indicating the reconstructed position shall be used in the fit.

Definition at line 124 of file CDCFitter2D.h.

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

tracking/trackFindingCDC/fitting/include/CDCRiemannFitter.h
tracking/trackFindingCDC/fitting/src/CDCRiemannFitter.cc

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Function Documentation

◆ fit() [1/9]

◆ fit() [2/9]

◆ fit() [3/9]

◆ fit() [4/9]

◆ fit() [5/9]

◆ fit() [6/9]

◆ fit() [7/9]

◆ fit() [8/9]

◆ fit() [9/9]

◆ fitGeneric() [1/2]

◆ fitGeneric() [2/2]

◆ fitInternal()

◆ getFitter() [1/2]

◆ getFitter() [2/2]

◆ getLineFitter()

◆ getOriginCircleFitter()

◆ isLineConstrained()

◆ isOriginConstrained()

◆ setFitVariance()

◆ setLineConstrained()

◆ setOriginConstrained()

◆ update() [1/5]

◆ update() [2/5]

◆ update() [3/5]

◆ update() [4/5]

◆ update() [5/5]

◆ updateGeneric() [1/2]

◆ updateGeneric() [2/2]

◆ useOnlyOrientation()

◆ useOnlyPosition()

◆ usePositionAndOrientation()

Member Data Documentation

◆ m_fitVariance

◆ m_lineConstrained

◆ m_originConstrained

◆ m_useOrientation

◆ m_usePosition