GENFIT  Rev:NoNumberAvailable
Public Member Functions | Private Attributes | List of all members
gbl::GblPoint Class Reference

Point on trajectory. More...

#include <GblPoint.h>

Collaboration diagram for gbl::GblPoint:
Collaboration graph

Public Member Functions

 GblPoint (const TMatrixD &aJacobian)
 Create a point. More...
 
 GblPoint (const SMatrix55 &aJacobian)
 
virtual ~GblPoint ()
 
void addMeasurement (const TMatrixD &aProjection, const TVectorD &aResiduals, const TVectorD &aPrecision, double minPrecision=0.)
 Add a measurement to a point. More...
 
void addMeasurement (const TMatrixD &aProjection, const TVectorD &aResiduals, const TMatrixDSym &aPrecision, double minPrecision=0.)
 Add a measurement to a point. More...
 
void addMeasurement (const TVectorD &aResiduals, const TVectorD &aPrecision, double minPrecision=0.)
 Add a measurement to a point. More...
 
void addMeasurement (const TVectorD &aResiduals, const TMatrixDSym &aPrecision, double minPrecision=0.)
 Add a measurement to a point. More...
 
unsigned int hasMeasurement () const
 Check for measurement at a point. More...
 
void getMeasurement (SMatrix55 &aProjection, SVector5 &aResiduals, SVector5 &aPrecision) const
 Retrieve measurement of a point. More...
 
void getMeasTransformation (TMatrixD &aTransformation) const
 Get measurement transformation (from diagonalization). More...
 
void addScatterer (const TVectorD &aResiduals, const TVectorD &aPrecision)
 Add a (thin) scatterer to a point. More...
 
void addScatterer (const TVectorD &aResiduals, const TMatrixDSym &aPrecision)
 Add a (thin) scatterer to a point. More...
 
bool hasScatterer () const
 Check for scatterer at a point. More...
 
void getScatterer (SMatrix22 &aTransformation, SVector2 &aResiduals, SVector2 &aPrecision) const
 Retrieve scatterer of a point. More...
 
void getScatTransformation (TMatrixD &aTransformation) const
 Get scatterer transformation (from diagonalization). More...
 
void addLocals (const TMatrixD &aDerivatives)
 Add local derivatives to a point. More...
 
unsigned int getNumLocals () const
 Retrieve number of local derivatives from a point. More...
 
const TMatrixD & getLocalDerivatives () const
 Retrieve local derivatives from a point. More...
 
void addGlobals (const std::vector< int > &aLabels, const TMatrixD &aDerivatives)
 Add global derivatives to a point. More...
 
unsigned int getNumGlobals () const
 Retrieve number of global derivatives from a point. More...
 
std::vector< int > getGlobalLabels () const
 Retrieve global derivatives labels from a point. More...
 
const TMatrixD & getGlobalDerivatives () const
 Retrieve global derivatives from a point. More...
 
void setLabel (unsigned int aLabel)
 Define label of point (by GBLTrajectory constructor) More...
 
unsigned int getLabel () const
 Retrieve label of point. More...
 
void setOffset (int anOffset)
 Define offset for point (by GBLTrajectory constructor) More...
 
int getOffset () const
 Retrieve offset for point. More...
 
const SMatrix55getP2pJacobian () const
 Retrieve point-to-(previous)point jacobian. More...
 
void addPrevJacobian (const SMatrix55 &aJac)
 Define jacobian to previous scatterer (by GBLTrajectory constructor) More...
 
void addNextJacobian (const SMatrix55 &aJac)
 Define jacobian to next scatterer (by GBLTrajectory constructor) More...
 
void getDerivatives (int aDirection, SMatrix22 &matW, SMatrix22 &matWJ, SVector2 &vecWd) const
 Retrieve derivatives of local track model. More...
 
void printPoint (unsigned int level=0) const
 Print GblPoint. More...
 

Private Attributes

unsigned int theLabel
 Label identifying point. More...
 
int theOffset
 Offset number at point if not negative (else interpolation needed) More...
 
SMatrix55 p2pJacobian
 Point-to-point jacobian from previous point. More...
 
SMatrix55 prevJacobian
 Jacobian to previous scatterer (or first measurement) More...
 
SMatrix55 nextJacobian
 Jacobian to next scatterer (or last measurement) More...
 
unsigned int measDim
 Dimension of measurement (1-5), 0 indicates absence of measurement. More...
 
SMatrix55 measProjection
 Projection from measurement to local system. More...
 
SVector5 measResiduals
 Measurement residuals. More...
 
SVector5 measPrecision
 Measurement precision (diagonal of inverse covariance matrix) More...
 
bool transFlag
 Transformation exists? More...
 
TMatrixD measTransformation
 Transformation of diagonalization (of meas. precision matrix) More...
 
bool scatFlag
 Scatterer present? More...
 
SMatrix22 scatTransformation
 Transformation of diagonalization (of scat. precision matrix) More...
 
SVector2 scatResiduals
 Scattering residuals (initial kinks if iterating) More...
 
SVector2 scatPrecision
 Scattering precision (diagonal of inverse covariance matrix) More...
 
TMatrixD localDerivatives
 Derivatives of measurement vs additional local (fit) parameters. More...
 
std::vector< int > globalLabels
 Labels of global (MP-II) derivatives. More...
 
TMatrixD globalDerivatives
 Derivatives of measurement vs additional global (MP-II) parameters. More...
 

Detailed Description

Point on trajectory.

User supplied point on (initial) trajectory.

Must have jacobian for propagation from previous point. May have:

  1. Measurement (1D - 5D)
  2. Scatterer (thin, 2D kinks)
  3. Additional local parameters (with derivatives). Fitted together with track parameters.
  4. Additional global parameters (with labels and derivatives). Not fitted, only passed on to (binary) file for fitting with Millepede-II.

Definition at line 48 of file GblPoint.h.

Constructor & Destructor Documentation

◆ GblPoint() [1/2]

gbl::GblPoint::GblPoint ( const TMatrixD &  aJacobian)

Create a point.

Create point on (initial) trajectory. Needs transformation jacobian from previous point.

Parameters
[in]aJacobianTransformation jacobian from previous point

Definition at line 18 of file GblPoint.cc.

References p2pJacobian.

◆ GblPoint() [2/2]

gbl::GblPoint::GblPoint ( const SMatrix55 aJacobian)

Definition at line 29 of file GblPoint.cc.

◆ ~GblPoint()

gbl::GblPoint::~GblPoint ( )
virtual

Definition at line 35 of file GblPoint.cc.

Member Function Documentation

◆ addGlobals()

void gbl::GblPoint::addGlobals ( const std::vector< int > &  aLabels,
const TMatrixD &  aDerivatives 
)

Add global derivatives to a point.

Point needs to have a measurement.

Parameters
[in]aLabelsGlobal derivatives labels
[in]aDerivativesGlobal derivatives (matrix)

Definition at line 298 of file GblPoint.cc.

References globalDerivatives, globalLabels, measDim, measTransformation, and transFlag.

Referenced by genfit::GblFitterInfo::constructGblPoint(), and genfit::GFGbl::processTrackWithRep().

◆ addLocals()

void gbl::GblPoint::addLocals ( const TMatrixD &  aDerivatives)

Add local derivatives to a point.

Point needs to have a measurement.

Parameters
[in]aDerivativesLocal derivatives (matrix)

Definition at line 271 of file GblPoint.cc.

References localDerivatives, measDim, measTransformation, and transFlag.

Referenced by genfit::GblFitterInfo::constructGblPoint().

◆ addMeasurement() [1/4]

void gbl::GblPoint::addMeasurement ( const TMatrixD &  aProjection,
const TVectorD &  aResiduals,
const TMatrixDSym &  aPrecision,
double  minPrecision = 0. 
)

Add a measurement to a point.

Add measurement (in meas. system) with arbitrary precision (inverse covariance) matrix. Will be diagonalized. ((up to) 2D: position, 4D: slope+position, 5D: curvature+slope+position)

Parameters
[in]aProjectionProjection from local to measurement system
[in]aResidualsMeasurement residuals
[in]aPrecisionMeasurement precision (matrix)
[in]minPrecisionMinimal precision to accept measurement

Definition at line 72 of file GblPoint.cc.

References measDim, measPrecision, measProjection, measResiduals, measTransformation, and transFlag.

◆ addMeasurement() [2/4]

void gbl::GblPoint::addMeasurement ( const TMatrixD &  aProjection,
const TVectorD &  aResiduals,
const TVectorD &  aPrecision,
double  minPrecision = 0. 
)

Add a measurement to a point.

Add measurement (in meas. system) with diagonal precision (inverse covariance) matrix. ((up to) 2D: position, 4D: slope+position, 5D: curvature+slope+position)

Parameters
[in]aProjectionProjection from local to measurement system
[in]aResidualsMeasurement residuals
[in]aPrecisionMeasurement precision (diagonal)
[in]minPrecisionMinimal precision to accept measurement

Definition at line 47 of file GblPoint.cc.

References measDim, measPrecision, measProjection, and measResiduals.

Referenced by genfit::GblFitterInfo::constructGblPoint(), and genfit::GFGbl::processTrackWithRep().

◆ addMeasurement() [3/4]

void gbl::GblPoint::addMeasurement ( const TVectorD &  aResiduals,
const TMatrixDSym &  aPrecision,
double  minPrecision = 0. 
)

Add a measurement to a point.

Add measurement in local system with arbitrary precision (inverse covariance) matrix. Will be diagonalized. ((up to) 2D: position, 4D: slope+position, 5D: curvature+slope+position)

Parameters
[in]aResidualsMeasurement residuals
[in]aPrecisionMeasurement precision (matrix)
[in]minPrecisionMinimal precision to accept measurement

Definition at line 124 of file GblPoint.cc.

References measDim, measPrecision, measProjection, measResiduals, measTransformation, and transFlag.

◆ addMeasurement() [4/4]

void gbl::GblPoint::addMeasurement ( const TVectorD &  aResiduals,
const TVectorD &  aPrecision,
double  minPrecision = 0. 
)

Add a measurement to a point.

Add measurement in local system with diagonal precision (inverse covariance) matrix. ((up to) 2D: position, 4D: slope+position, 5D: curvature+slope+position)

Parameters
[in]aResidualsMeasurement residuals
[in]aPrecisionMeasurement precision (diagonal)
[in]minPrecisionMinimal precision to accept measurement

Definition at line 103 of file GblPoint.cc.

References measDim, measPrecision, measProjection, and measResiduals.

◆ addNextJacobian()

void gbl::GblPoint::addNextJacobian ( const SMatrix55 aJac)

Define jacobian to next scatterer (by GBLTrajectory constructor)

Parameters
[in]aJacJacobian

Definition at line 379 of file GblPoint.cc.

References nextJacobian.

Referenced by gbl::GblTrajectory::calcJacobians().

◆ addPrevJacobian()

void gbl::GblPoint::addPrevJacobian ( const SMatrix55 aJac)

Define jacobian to previous scatterer (by GBLTrajectory constructor)

Parameters
[in]aJacJacobian

Definition at line 362 of file GblPoint.cc.

References prevJacobian.

◆ addScatterer() [1/2]

void gbl::GblPoint::addScatterer ( const TVectorD &  aResiduals,
const TMatrixDSym &  aPrecision 
)

Add a (thin) scatterer to a point.

Add scatterer with arbitrary precision (inverse covariance) matrix. Will be diagonalized. Changes local track direction.

The precision matrix for the local slopes is defined by the angular scattering error theta_0 and the scalar products c_1, c_2 of the offset directions in the local frame with the track direction:

       (1 - c_1*c_1 - c_2*c_2)   |  1 - c_1*c_1     - c_1*c_2  |
  P =  ~~~~~~~~~~~~~~~~~~~~~~~ * |                             |
           theta_0*theta_0       |    - c_1*c_2   1 - c_2*c_2  |
Parameters
[in]aResidualsScatterer residuals
[in]aPrecisionScatterer precision (matrix)

Definition at line 214 of file GblPoint.cc.

References scatFlag, scatPrecision, scatResiduals, and scatTransformation.

◆ addScatterer() [2/2]

void gbl::GblPoint::addScatterer ( const TVectorD &  aResiduals,
const TVectorD &  aPrecision 
)

Add a (thin) scatterer to a point.

Add scatterer with diagonal precision (inverse covariance) matrix. Changes local track direction.

Parameters
[in]aResidualsScatterer residuals
[in]aPrecisionScatterer precision (diagonal of inverse covariance matrix)

Definition at line 188 of file GblPoint.cc.

References scatFlag, scatPrecision, scatResiduals, and scatTransformation.

Referenced by genfit::GblFitterInfo::constructGblPoint(), and genfit::GFGbl::processTrackWithRep().

◆ getDerivatives()

void gbl::GblPoint::getDerivatives ( int  aDirection,
SMatrix22 matW,
SMatrix22 matWJ,
SVector2 vecWd 
) const

Retrieve derivatives of local track model.

Linearized track model: F_u(q/p,u',u) = J*u + S*u' + d*q/p, W is inverse of S, negated for backward propagation.

Parameters
[in]aDirectionPropagation direction (>0 forward, else backward)
[out]matWW
[out]matWJW*J
[out]vecWdW*d
Exceptions
std::overflow_error: matrix S is singular.

Definition at line 393 of file GblPoint.cc.

References nextJacobian, and prevJacobian.

Referenced by gbl::GblTrajectory::getFitToKinkJacobian(), and gbl::GblTrajectory::getFitToLocalJacobian().

◆ getGlobalDerivatives()

const TMatrixD & gbl::GblPoint::getGlobalDerivatives ( ) const

Retrieve global derivatives from a point.

Definition at line 323 of file GblPoint.cc.

References globalDerivatives.

◆ getGlobalLabels()

std::vector< int > gbl::GblPoint::getGlobalLabels ( ) const

Retrieve global derivatives labels from a point.

Definition at line 318 of file GblPoint.cc.

References globalLabels.

◆ getLabel()

unsigned int gbl::GblPoint::getLabel ( ) const

Retrieve label of point.

Definition at line 336 of file GblPoint.cc.

References theLabel.

◆ getLocalDerivatives()

const TMatrixD & gbl::GblPoint::getLocalDerivatives ( ) const

Retrieve local derivatives from a point.

Definition at line 288 of file GblPoint.cc.

References localDerivatives.

◆ getMeasTransformation()

void gbl::GblPoint::getMeasTransformation ( TMatrixD &  aTransformation) const

Get measurement transformation (from diagonalization).

Parameters
[out]aTransformationTransformation matrix

Definition at line 171 of file GblPoint.cc.

References measDim, measTransformation, and transFlag.

◆ getMeasurement()

void gbl::GblPoint::getMeasurement ( SMatrix55 aProjection,
SVector5 aResiduals,
SVector5 aPrecision 
) const

Retrieve measurement of a point.

Parameters
[out]aProjectionProjection from (diagonalized) measurement to local system
[out]aResidualsMeasurement residuals
[out]aPrecisionMeasurement precision (diagonal)

Definition at line 160 of file GblPoint.cc.

References measPrecision, measProjection, and measResiduals.

◆ getNumGlobals()

unsigned int gbl::GblPoint::getNumGlobals ( ) const

Retrieve number of global derivatives from a point.

Definition at line 313 of file GblPoint.cc.

References globalDerivatives.

◆ getNumLocals()

unsigned int gbl::GblPoint::getNumLocals ( ) const

Retrieve number of local derivatives from a point.

Definition at line 283 of file GblPoint.cc.

References localDerivatives.

◆ getOffset()

int gbl::GblPoint::getOffset ( ) const

Retrieve offset for point.

Definition at line 349 of file GblPoint.cc.

References theOffset.

Referenced by gbl::GblTrajectory::getFitToKinkJacobian(), and gbl::GblTrajectory::getFitToLocalJacobian().

◆ getP2pJacobian()

const SMatrix55 & gbl::GblPoint::getP2pJacobian ( ) const

Retrieve point-to-(previous)point jacobian.

Definition at line 354 of file GblPoint.cc.

References p2pJacobian.

Referenced by gbl::GblTrajectory::calcJacobians().

◆ getScatterer()

void gbl::GblPoint::getScatterer ( SMatrix22 aTransformation,
SVector2 aResiduals,
SVector2 aPrecision 
) const

Retrieve scatterer of a point.

Parameters
[out]aTransformationScatterer transformation from diagonalization
[out]aResidualsScatterer residuals
[out]aPrecisionScatterer precision (diagonal)

Definition at line 242 of file GblPoint.cc.

References scatPrecision, scatResiduals, and scatTransformation.

◆ getScatTransformation()

void gbl::GblPoint::getScatTransformation ( TMatrixD &  aTransformation) const

Get scatterer transformation (from diagonalization).

Parameters
[out]aTransformationTransformation matrix

Definition at line 253 of file GblPoint.cc.

References scatFlag, and scatTransformation.

◆ hasMeasurement()

unsigned int gbl::GblPoint::hasMeasurement ( ) const

Check for measurement at a point.

Get dimension of measurement (0 = none).

Returns
measurement dimension

Definition at line 150 of file GblPoint.cc.

References measDim.

Referenced by genfit::GblFitter::processTrackWithRep().

◆ hasScatterer()

bool gbl::GblPoint::hasScatterer ( ) const

Check for scatterer at a point.

Definition at line 232 of file GblPoint.cc.

References scatFlag.

Referenced by genfit::GblFitter::processTrackWithRep().

◆ printPoint()

void gbl::GblPoint::printPoint ( unsigned int  level = 0) const

Print GblPoint.

Parameters
[in]levelprint level (0: minimum, >0: more)

Definition at line 425 of file GblPoint.cc.

References globalDerivatives, globalLabels, localDerivatives, measDim, measPrecision, measProjection, measResiduals, nextJacobian, p2pJacobian, prevJacobian, scatFlag, scatPrecision, scatResiduals, theLabel, theOffset, and transFlag.

◆ setLabel()

void gbl::GblPoint::setLabel ( unsigned int  aLabel)

Define label of point (by GBLTrajectory constructor)

Parameters
[in]aLabelLabel identifying point

Definition at line 331 of file GblPoint.cc.

References theLabel.

◆ setOffset()

void gbl::GblPoint::setOffset ( int  anOffset)

Define offset for point (by GBLTrajectory constructor)

Parameters
[in]anOffsetOffset number

Definition at line 344 of file GblPoint.cc.

References theOffset.

Member Data Documentation

◆ globalDerivatives

TMatrixD gbl::GblPoint::globalDerivatives
private

Derivatives of measurement vs additional global (MP-II) parameters.

Definition at line 109 of file GblPoint.h.

Referenced by addGlobals(), getGlobalDerivatives(), getNumGlobals(), and printPoint().

◆ globalLabels

std::vector<int> gbl::GblPoint::globalLabels
private

Labels of global (MP-II) derivatives.

Definition at line 108 of file GblPoint.h.

Referenced by addGlobals(), getGlobalLabels(), and printPoint().

◆ localDerivatives

TMatrixD gbl::GblPoint::localDerivatives
private

Derivatives of measurement vs additional local (fit) parameters.

Definition at line 107 of file GblPoint.h.

Referenced by addLocals(), getLocalDerivatives(), getNumLocals(), and printPoint().

◆ measDim

unsigned int gbl::GblPoint::measDim
private

Dimension of measurement (1-5), 0 indicates absence of measurement.

Definition at line 97 of file GblPoint.h.

Referenced by addGlobals(), addLocals(), addMeasurement(), getMeasTransformation(), hasMeasurement(), and printPoint().

◆ measPrecision

SVector5 gbl::GblPoint::measPrecision
private

Measurement precision (diagonal of inverse covariance matrix)

Definition at line 100 of file GblPoint.h.

Referenced by addMeasurement(), getMeasurement(), and printPoint().

◆ measProjection

SMatrix55 gbl::GblPoint::measProjection
private

Projection from measurement to local system.

Definition at line 98 of file GblPoint.h.

Referenced by addMeasurement(), getMeasurement(), and printPoint().

◆ measResiduals

SVector5 gbl::GblPoint::measResiduals
private

Measurement residuals.

Definition at line 99 of file GblPoint.h.

Referenced by addMeasurement(), getMeasurement(), and printPoint().

◆ measTransformation

TMatrixD gbl::GblPoint::measTransformation
private

Transformation of diagonalization (of meas. precision matrix)

Definition at line 102 of file GblPoint.h.

Referenced by addGlobals(), addLocals(), addMeasurement(), and getMeasTransformation().

◆ nextJacobian

SMatrix55 gbl::GblPoint::nextJacobian
private

Jacobian to next scatterer (or last measurement)

Definition at line 96 of file GblPoint.h.

Referenced by addNextJacobian(), getDerivatives(), and printPoint().

◆ p2pJacobian

SMatrix55 gbl::GblPoint::p2pJacobian
private

Point-to-point jacobian from previous point.

Definition at line 94 of file GblPoint.h.

Referenced by GblPoint(), getP2pJacobian(), and printPoint().

◆ prevJacobian

SMatrix55 gbl::GblPoint::prevJacobian
private

Jacobian to previous scatterer (or first measurement)

Definition at line 95 of file GblPoint.h.

Referenced by addPrevJacobian(), getDerivatives(), and printPoint().

◆ scatFlag

bool gbl::GblPoint::scatFlag
private

Scatterer present?

Definition at line 103 of file GblPoint.h.

Referenced by addScatterer(), getScatTransformation(), hasScatterer(), and printPoint().

◆ scatPrecision

SVector2 gbl::GblPoint::scatPrecision
private

Scattering precision (diagonal of inverse covariance matrix)

Definition at line 106 of file GblPoint.h.

Referenced by addScatterer(), getScatterer(), and printPoint().

◆ scatResiduals

SVector2 gbl::GblPoint::scatResiduals
private

Scattering residuals (initial kinks if iterating)

Definition at line 105 of file GblPoint.h.

Referenced by addScatterer(), getScatterer(), and printPoint().

◆ scatTransformation

SMatrix22 gbl::GblPoint::scatTransformation
private

Transformation of diagonalization (of scat. precision matrix)

Definition at line 104 of file GblPoint.h.

Referenced by addScatterer(), getScatterer(), and getScatTransformation().

◆ theLabel

unsigned int gbl::GblPoint::theLabel
private

Label identifying point.

Definition at line 92 of file GblPoint.h.

Referenced by getLabel(), printPoint(), and setLabel().

◆ theOffset

int gbl::GblPoint::theOffset
private

Offset number at point if not negative (else interpolation needed)

Definition at line 93 of file GblPoint.h.

Referenced by getOffset(), printPoint(), and setOffset().

◆ transFlag

bool gbl::GblPoint::transFlag
private

Transformation exists?

Definition at line 101 of file GblPoint.h.

Referenced by addGlobals(), addLocals(), addMeasurement(), getMeasTransformation(), and printPoint().


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