Class implementing the fitter using Karimakis method. More...

#include <CDCKarimakiFitter.h>

Inheritance diagram for CDCKarimakiFitter:

Public Member Functions
	CDCKarimakiFitter ()
	Default constructor.

	~CDCKarimakiFitter ()
	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 indictor that lines should be fitted by this fitter.

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

Static Public Member Functions
static const CDCKarimakiFitter &	getFitter ()
	Static getter for a general fitter instance with Karimakis method.

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

static const CDCKarimakiFitter &	getNoDriftVarianceFitter ()
	Static getter for a general fitter that does not use the drift length variances.

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

CDCTrajectory2D	fitGeneric (const AStartHits &startHits, const AEndHits &endHits) const
	Fits together two collections of hit types which are convertable 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 convertable 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 convertable 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.

Detailed Description

Class implementing the fitter using Karimakis method.

Definition at line 23 of file CDCKarimakiFitter.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 convertable 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 convertable 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 214 of file KarimakisMethod.cc.

{
  // Matrix of weighted sums
  Eigen::Matrix< double, 4, 4> sNoL = getWXYRSumMatrix(observations2D);
 
  // 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();
 
  // Determine NDF : Circle fit eats up to 3 degrees of freedom debpending on the constraints
  size_t ndf = observations2D.size() - 2;
 
  if (not isLineConstrained()) {
    --ndf;
  }
 
  // Parameters to be fitted
  UncertainPerigeeCircle resultCircle = fitKarimaki(sNoL(iW), meansNoL, cNoL, isLineConstrained());
  double chi2 = calcChi2Karimaki(resultCircle, sNoL(iW), cNoL);
 
  PerigeePrecision perigeePrecision = calcPrecisionKarimaki(resultCircle, sNoL, isLineConstrained());
 
  // Use in pivotingin caset the matrix is not full rank as is for the constrained cases-
  PerigeeCovariance perigeeCovariance = PerigeeUtil::covarianceFromPrecision(perigeePrecision);
 
  resultCircle.setChi2(chi2);
  resultCircle.setNDF(ndf);
  resultCircle.setPerigeeCovariance(perigeeCovariance);
  return resultCircle;
}

◆ getFitter()

const CDCKarimakiFitter & getFitter ( )

static

Static getter for a general fitter instance with Karimakis method.

Definition at line 20 of file CDCKarimakiFitter.cc.

{
  static CDCKarimakiFitter fitter;
  return fitter;
}

◆ getLineFitter()

const CDCKarimakiFitter & getLineFitter ( )

static

Static getter for a line fitter.

Definition at line 26 of file CDCKarimakiFitter.cc.

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

◆ getNoDriftVarianceFitter()

const CDCKarimakiFitter & getNoDriftVarianceFitter ( )

static

Static getter for a general fitter that does not use the drift length variances.

Definition at line 33 of file CDCKarimakiFitter.cc.

{
  static CDCKarimakiFitter noDriftVarianceFitter;
  noDriftVarianceFitter.setFitVariance(EFitVariance::c_DriftLength);
  return noDriftVarianceFitter;
}

◆ isLineConstrained()

bool isLineConstrained ( ) const

inlineinherited

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

Definition at line 36 of file KarimakisMethod.h.

      {
        return m_lineConstrained;
      }

◆ 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 42 of file KarimakisMethod.h.

      {
        m_lineConstrained = 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 29 of file KarimakisMethod.cc.

{
  size_t nObservations = observations2D.size();
  trajectory2D.clear();
  if (not nObservations) return;
 
  Vector2D ref = observations2D.getCentralPoint();
  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.setLocalOrigin(ref);
  trajectory2D.setLocalCircle(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 convertable 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 convertable 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 49 of file KarimakisMethod.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/CDCKarimakiFitter.h
tracking/trackFindingCDC/fitting/src/CDCKarimakiFitter.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()

◆ getLineFitter()

◆ getNoDriftVarianceFitter()

◆ isLineConstrained()

◆ setFitVariance()

◆ setLineConstrained()

◆ 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_useOrientation

◆ m_usePosition