Point on trajectory. More...

#include <GblPoint.h>

Public Member Functions
	GblPoint (const TMatrixD &aJacobian)
	Create a point.
virtual	~GblPoint ()
void	addMeasurement (const TMatrixD &aProjection, const TVectorD &aResiduals, const TVectorD &aPrecision)
	Add a mesurement to a point.
void	addMeasurement (const TVectorD &aResiduals, const TMatrixDSym &aPrecision)
	Add a mesurement to a point.
unsigned int	hasMeasurement () const
	Check for measurement at a point.
void	getMeasurement (SMatrix55 &aProjection, SVector5 &aResiduals, SVector5 &aPrecision) const
	Retrieve measurement of a point.
void	addScatterer (const TVectorD &aResiduals, const TVectorD &aPrecision)
	Add a (thin) scatterer to a point.
bool	hasScatterer () const
	Check for scatterer at a point.
void	getScatterer (SVector2 &aResiduals, SVector2 &aPrecision) const
	Retrieve scatterer of a point.
void	addLocals (const TMatrixD &aDerivatives)
	Add local derivatives to a point.
unsigned int	getNumLocals () const
	Retrieve number of local derivatives from a point.
TMatrixD	getLocalDerivatives () const
	Retrieve local derivatives from a point.
void	addGlobals (const std::vector< int > &aLabels, const TMatrixD &aDerivatives)
	Add global derivatives to a point.
unsigned int	getNumGlobals () const
	Retrieve number of global derivatives from a point.
std::vector< int >	getGlobalLabels () const
	Retrieve global derivatives labels from a point.
TMatrixD	getGlobalDerivatives () const
	Retrieve global derivatives from a point.
void	setLabel (unsigned int aLabel)
	Define label of point.
unsigned int	getLabel () const
	Retrieve label of point.
void	setOffset (int anOffset)
	Define offset for point.
int	getOffset () const
	Retrieve offset for point.
SMatrix55	getP2pJacobian () const
	Retrieve point-to-(previous)point jacobian.
void	addPrevJacobian (const SMatrix55 aJac)
	Define jacobian to previous scatterer.
void	addNextJacobian (const SMatrix55 aJac)
	Define jacobian to next scatterer.
void	getDerivatives (int aDirection, SMatrix22 &matW, SMatrix22 &matWJ, SVector2 &vecWd) const
	Retrieve derivatives of local track model.
Private Attributes
unsigned int	theLabel
	Label identifying point.
int	theOffset
	Offset number at point if not negative (else interpolation needed)
SMatrix55	p2pJacobian
	Point-to-point jacobian from previous point.
SMatrix55	prevJacobian
	Jacobian to previous scatterer (or first measurement)
SMatrix55	nextJacobian
	Jacobian to next scatterer (or last measurement)
unsigned int	measDim
	Dimension of measurement (1-5), 0 indicates absence of measurement.
SMatrix55	measProjection
	Projection from measurement to local system.
SVector5	measResiduals
	Measurement residuals.
SVector5	measPrecision
	Measurement precision (diagonal of inverse covariance matrix)
bool	transFlag
	Transformation exists?
TMatrixD	measTransformation
	Transformation of diagonalization (of precision matrix)
bool	scatFlag
	Scatterer present?
SVector2	scatResiduals
	Scattering residuals (initial kinks if iterating)
SVector2	scatPrecision
	Scattering precision (diagonal of inverse covariance matrix)
TMatrixD	localDerivatives
	Derivatives of measurement vs additional local (fit) parameters.
std::vector< int >	globalLabels
	Labels of global (MP-II) derivatives.
TMatrixD	globalDerivatives
	Derivatives of measurement vs additional global (MP-II) parameters.

Detailed Description

Point on trajectory.

User supplied point on (initial) trajectory.

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

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

Definition at line 43 of file GblPoint.h.

Constructor & Destructor Documentation

GblPoint::GblPoint ( const TMatrixD & aJacobian )

Create a point.

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

Parameters:

[in] aJacobian Transformation jacobian from previous point

Definition at line 15 of file GblPoint.cpp.

References p2pJacobian.

GblPoint::~GblPoint ( ) [virtual]

Definition at line 26 of file GblPoint.cpp.

Member Function Documentation

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

Add global derivatives to a point.

Point needs to have a measurement.

Parameters:

[in]	aLabels	Global derivatives labels
[in]	aDerivatives	Global derivatives (matrix)

Definition at line 161 of file GblPoint.cpp.

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

void GblPoint::addLocals ( const TMatrixD & aDerivatives )

Add local derivatives to a point.

Point needs to have a measurement.

Parameters:

[in] aDerivatives Local derivatives (matrix)

Definition at line 134 of file GblPoint.cpp.

References localDerivatives, measDim, measTransformation, and transFlag.

void GblPoint::addMeasurement	(	const TMatrixD &	aProjection,
		const TVectorD &	aResiduals,
		const TVectorD &	aPrecision
	)

Add a mesurement to a point.

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

Parameters:

[in]	aProjection	Projection from measurement to local system
[in]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (diagonal)

Definition at line 38 of file GblPoint.cpp.

References measDim, measPrecision, measProjection, and measResiduals.

Referenced by example1().

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

Add a mesurement to a point.

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

Parameters:

[in]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (matrix)

Definition at line 59 of file GblPoint.cpp.

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

void GblPoint::addNextJacobian ( const SMatrix55 aJac )

Define jacobian to next scatterer.

Parameters:

[in] aJac Jacobian

Definition at line 242 of file GblPoint.cpp.

References nextJacobian.

void GblPoint::addPrevJacobian ( const SMatrix55 aJac )

Define jacobian to previous scatterer.

Parameters:

[in] aJac Jacobian

Definition at line 225 of file GblPoint.cpp.

References prevJacobian.

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

Add a (thin) scatterer to a point.

Parameters:

[in]	aResiduals	Scatterer residuals
[in]	aPrecision	Scatterer precision (diagonal of inverse covariance matrix)

Definition at line 105 of file GblPoint.cpp.

References scatFlag, scatPrecision, and scatResiduals.

Referenced by example1().

void 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]	aDirection	Propagation direction (>0 forward, else backward)
[out]	matW	W
[out]	matWJ	W*J
[out]	vecWd	W*d

Definition at line 255 of file GblPoint.cpp.

References nextJacobian, and prevJacobian.

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

TMatrixD GblPoint::getGlobalDerivatives ( ) const

Retrieve global derivatives from a point.

Definition at line 186 of file GblPoint.cpp.

References globalDerivatives.

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

Retrieve global derivatives labels from a point.

Definition at line 181 of file GblPoint.cpp.

References globalLabels.

unsigned int GblPoint::getLabel ( ) const

Retrieve label of point.

Definition at line 199 of file GblPoint.cpp.

References theLabel.

TMatrixD GblPoint::getLocalDerivatives ( ) const

Retrieve local derivatives from a point.

Definition at line 151 of file GblPoint.cpp.

References localDerivatives.

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

Retrieve measurement of a point.

Parameters:

[out]	aProjection	Projection from (diagonalized) measurement to local system
[out]	aResiduals	Measurement residuals
[out]	aPrecision	Measurement precision (diagonal)

Definition at line 93 of file GblPoint.cpp.

References measPrecision, measProjection, and measResiduals.

unsigned int GblPoint::getNumGlobals ( ) const

Retrieve number of global derivatives from a point.

Definition at line 176 of file GblPoint.cpp.

References globalDerivatives.

unsigned int GblPoint::getNumLocals ( ) const

Retrieve number of local derivatives from a point.

Definition at line 146 of file GblPoint.cpp.

References localDerivatives.

Referenced by GblTrajectory::addPoint().

int GblPoint::getOffset ( ) const