Point on trajectory. More...

#include <GblPoint.h>

Public Member Functions
void	addGlobals (const std::vector< int > &aLabels, const TMatrixD &aDerivatives)
	Add global derivatives to a point. More...

void	addLocals (const TMatrixD &aDerivatives)
	Add local derivatives to a point. More...

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...

void	addNextJacobian (const SMatrix55 &aJac)
	Define jacobian to next scatterer (by GBLTrajectory constructor) More...

void	addPrevJacobian (const SMatrix55 &aJac)
	Define jacobian to previous scatterer (by GBLTrajectory constructor) 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...

	GblPoint (const TMatrixD &aJacobian)
	Create a point. More...

	GblPoint (const SMatrix55 &aJacobian)

void	getDerivatives (int aDirection, SMatrix22 &matW, SMatrix22 &matWJ, SVector2 &vecWd) const
	Retrieve derivatives of local track model. More...

const TMatrixD &	getGlobalDerivatives () const
	Retrieve global derivatives from a point. More...

std::vector< int >	getGlobalLabels () const
	Retrieve global derivatives labels from a point. More...

unsigned int	getLabel () const
	Retrieve label of point. More...

const TMatrixD &	getLocalDerivatives () const
	Retrieve local derivatives from a point. More...

void	getMeasTransformation (TMatrixD &aTransformation) const
	Get measurement transformation (from diagonalization). More...

void	getMeasurement (SMatrix55 &aProjection, SVector5 &aResiduals, SVector5 &aPrecision) const
	Retrieve measurement of a point. More...

unsigned int	getNumGlobals () const
	Retrieve number of global derivatives from a point. More...

unsigned int	getNumLocals () const
	Retrieve number of local derivatives from a point. More...

int	getOffset () const
	Retrieve offset for point. More...

const SMatrix55 &	getP2pJacobian () const
	Retrieve point-to-(previous)point jacobian. 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...

unsigned int	hasMeasurement () const
	Check for measurement at a point. More...

bool	hasScatterer () const
	Check for scatterer at a point. More...

void	printPoint (unsigned int level=0) const
	Print GblPoint. More...

void	setLabel (unsigned int aLabel)
	Define label of point (by GBLTrajectory constructor) More...

void	setOffset (int anOffset)
	Define offset for point (by GBLTrajectory constructor) More...

virtual	~GblPoint ()

Private Attributes
TMatrixD	globalDerivatives
	Derivatives of measurement vs additional global (MP-II) parameters. More...

std::vector< int >	globalLabels
	Labels of global (MP-II) derivatives. More...

TMatrixD	localDerivatives
	Derivatives of measurement vs additional local (fit) parameters. More...

unsigned int	measDim
	Dimension of measurement (1-5), 0 indicates absence of measurement. More...

SVector5	measPrecision
	Measurement precision (diagonal of inverse covariance matrix) More...

SMatrix55	measProjection
	Projection from measurement to local system. More...

SVector5	measResiduals
	Measurement residuals. More...

TMatrixD	measTransformation
	Transformation of diagonalization (of meas. precision matrix) More...

SMatrix55	nextJacobian
	Jacobian to next scatterer (or last measurement) More...

SMatrix55	p2pJacobian
	Point-to-point jacobian from previous point. More...

SMatrix55	prevJacobian
	Jacobian to previous scatterer (or first measurement) More...

bool	scatFlag
	Scatterer present? More...

SVector2	scatPrecision
	Scattering precision (diagonal of inverse covariance matrix) More...

SVector2	scatResiduals
	Scattering residuals (initial kinks if iterating) More...

SMatrix22	scatTransformation
	Transformation of diagonalization (of scat. precision matrix) More...

unsigned int	theLabel
	Label identifying point. More...

int	theOffset
	Offset number at point if not negative (else interpolation needed) More...

bool	transFlag
	Transformation exists? More...

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 46 of file GblPoint.h.

Constructor & Destructor Documentation

gbl::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 18 of file GblPoint.cc.

References i, j, and p2pJacobian.

                                             :
         theLabel(0), theOffset(0), measDim(0), transFlag(false), measTransformation(), scatFlag(
                 false), localDerivatives(), globalLabels(), globalDerivatives() {
 
     for (unsigned int i = 0; i < 5; ++i) {
         for (unsigned int j = 0; j < 5; ++j) {
             p2pJacobian(i, j) = aJacobian(i, j);
         }
     }
 }

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

Definition at line 29 of file GblPoint.cc.

                                              :
         theLabel(0), theOffset(0), p2pJacobian(aJacobian), measDim(0), transFlag(
                 false), measTransformation(), scatFlag(false), localDerivatives(), globalLabels(), globalDerivatives() {
 
 }

gbl::GblPoint::~GblPoint ( )

virtual

Definition at line 35 of file GblPoint.cc.

35 {

36 }

Member Function Documentation

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]	aLabels	Global derivatives labels
[in]	aDerivatives	Global derivatives (matrix)

Definition at line 298 of file GblPoint.cc.

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

Referenced by MillePedeAlignmentAlgorithm::addGlobalData().

                                       {
     if (measDim) {
         globalLabels = aLabels;
         globalDerivatives.ResizeTo(aDerivatives);
         if (transFlag) {
             globalDerivatives = measTransformation * aDerivatives;
         } else {
             globalDerivatives = aDerivatives;
         }
 
     }
 }

void gbl::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 271 of file GblPoint.cc.

References localDerivatives, measDim, measTransformation, and transFlag.

                                                      {
     if (measDim) {
         localDerivatives.ResizeTo(aDerivatives);
         if (transFlag) {
             localDerivatives = measTransformation * aDerivatives;
         } else {
             localDerivatives = aDerivatives;
         }
     }
 }

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]	aProjection	Projection from local to measurement system
[in]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (diagonal)
[in]	minPrecision	Minimal precision to accept measurement

Definition at line 47 of file GblPoint.cc.

References i, j, measDim, measPrecision, measProjection, and measResiduals.

Referenced by ReferenceTrajectory::addMaterialEffectsCurvlinGbl(), and ReferenceTrajectory::addMaterialEffectsLocalGbl().

                              {
     measDim = aResiduals.GetNrows();
     unsigned int iOff = 5 - measDim;
     for (unsigned int i = 0; i < measDim; ++i) {
         measResiduals(iOff + i) = aResiduals[i];
         measPrecision(iOff + i) = (
                 aPrecision[i] >= minPrecision ? aPrecision[i] : 0.);
         for (unsigned int j = 0; j < measDim; ++j) {
             measProjection(iOff + i, iOff + j) = aProjection(i, j);
         }
     }
 }

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]	aProjection	Projection from local to measurement system
[in]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (matrix)
[in]	minPrecision	Minimal precision to accept measurement

Definition at line 72 of file GblPoint.cc.

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

                              {
     measDim = aResiduals.GetNrows();
     TMatrixDSymEigen measEigen(aPrecision);
     measTransformation.ResizeTo(measDim, measDim);
     measTransformation = measEigen.GetEigenVectors();
     measTransformation.T();
     transFlag = true;
     TVectorD transResiduals = measTransformation * aResiduals;
     TVectorD transPrecision = measEigen.GetEigenValues();
     TMatrixD transProjection = measTransformation * aProjection;
     unsigned int iOff = 5 - measDim;
     for (unsigned int i = 0; i < measDim; ++i) {
         measResiduals(iOff + i) = transResiduals[i];
         measPrecision(iOff + i) = (
                 transPrecision[i] >= minPrecision ? transPrecision[i] : 0.);
         for (unsigned int j = 0; j < measDim; ++j) {
             measProjection(iOff + i, iOff + j) = transProjection(i, j);
         }
     }
 }

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]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (diagonal)
[in]	minPrecision	Minimal precision to accept measurement

Definition at line 103 of file GblPoint.cc.

References i, measDim, measPrecision, measProjection, and measResiduals.

                                                          {
     measDim = aResiduals.GetNrows();
     unsigned int iOff = 5 - measDim;
     for (unsigned int i = 0; i < measDim; ++i) {
         measResiduals(iOff + i) = aResiduals[i];
         measPrecision(iOff + i) = (
                 aPrecision[i] >= minPrecision ? aPrecision[i] : 0.);
     }
     measProjection = ROOT::Math::SMatrixIdentity();
 }

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]	aResiduals	Measurement residuals
[in]	aPrecision	Measurement precision (matrix)
[in]	minPrecision	Minimal precision to accept measurement

Definition at line 124 of file GblPoint.cc.

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

                                                             {
     measDim = aResiduals.GetNrows();
     TMatrixDSymEigen measEigen(aPrecision);
     measTransformation.ResizeTo(measDim, measDim);
     measTransformation = measEigen.GetEigenVectors();
     measTransformation.T();
     transFlag = true;
     TVectorD transResiduals = measTransformation * aResiduals;
     TVectorD transPrecision = measEigen.GetEigenValues();
     unsigned int iOff = 5 - measDim;
     for (unsigned int i = 0; i < measDim; ++i) {
         measResiduals(iOff + i) = transResiduals[i];
         measPrecision(iOff + i) = (
                 transPrecision[i] >= minPrecision ? transPrecision[i] : 0.);
         for (unsigned int j = 0; j < measDim; ++j) {
             measProjection(iOff + i, iOff + j) = measTransformation(i, j);
         }
     }
 }

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

Define jacobian to next scatterer (by GBLTrajectory constructor)

Parameters

[in] aJac Jacobian

Definition at line 379 of file GblPoint.cc.

References nextJacobian.

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

                                                     {
     nextJacobian = aJac;
 }

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

Define jacobian to previous scatterer (by GBLTrajectory constructor)

Parameters

[in] aJac Jacobian

Definition at line 362 of file GblPoint.cc.

References prevJacobian.

                                                     {
     int ifail = 0;
 // to optimize: need only two last rows of inverse
 //  prevJacobian = aJac.InverseFast(ifail);
 //  block matrix algebra
     SMatrix23 CA = aJac.Sub<SMatrix23>(3, 0)
             * aJac.Sub<SMatrix33>(0, 0).InverseFast(ifail); // C*A^-1
     SMatrix22 DCAB = aJac.Sub<SMatrix22>(3, 3) - CA * aJac.Sub<SMatrix32>(0, 3); // D - C*A^-1 *B
     DCAB.InvertFast();
     prevJacobian.Place_at(DCAB, 3, 3);
     prevJacobian.Place_at(-DCAB * CA, 3, 0);
 }

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]	aResiduals	Scatterer residuals
[in]	aPrecision	Scatterer precision (diagonal of inverse covariance matrix)

Definition at line 188 of file GblPoint.cc.

References scatFlag, scatPrecision, scatResiduals, and scatTransformation.

Referenced by ReferenceTrajectory::addMaterialEffectsCurvlinGbl(), and ReferenceTrajectory::addMaterialEffectsLocalGbl().

                                     {
     scatFlag = true;
     scatResiduals(0) = aResiduals[0];
     scatResiduals(1) = aResiduals[1];
     scatPrecision(0) = aPrecision[0];
     scatPrecision(1) = aPrecision[1];
     scatTransformation = ROOT::Math::SMatrixIdentity();
 }

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]	aResiduals	Scatterer residuals
[in]	aPrecision	Scatterer precision (matrix)

Definition at line 214 of file GblPoint.cc.

References i, j, scatFlag, scatPrecision, scatResiduals, and scatTransformation.

                                        {
     scatFlag = true;
     TMatrixDSymEigen scatEigen(aPrecision);
     TMatrixD aTransformation = scatEigen.GetEigenVectors();
     aTransformation.T();
     TVectorD transResiduals = aTransformation * aResiduals;
     TVectorD transPrecision = scatEigen.GetEigenValues();
     for (unsigned int i = 0; i < 2; ++i) {
         scatResiduals(i) = transResiduals[i];
         scatPrecision(i) = transPrecision[i];
         for (unsigned int j = 0; j < 2; ++j) {
             scatTransformation(i, j) = aTransformation(i, j);
         }
     }
 }

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

Exceptions

std::overflow_error : matrix S is singular.

Definition at line 393 of file GblPoint.cc.

References gather_cfg::cout, nextJacobian, and prevJacobian.

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

                                {
 
     SMatrix22 matJ;
     SVector2 vecd;
     if (aDirection < 1) {
         matJ = prevJacobian.Sub<SMatrix22>(3, 3);
         matW = -prevJacobian.Sub<SMatrix22>(3, 1);
         vecd = prevJacobian.SubCol<SVector2>(0, 3);
     } else {
         matJ = nextJacobian.Sub<SMatrix22>(3, 3);
         matW = nextJacobian.Sub<SMatrix22>(3, 1);
         vecd = nextJacobian.SubCol<SVector2>(0, 3);
     }
 
     if (!matW.InvertFast()) {
         std::cout << " GblPoint::getDerivatives failed to invert matrix: "
                 << matW << std::endl;
         std::cout
                 << " Possible reason for singular matrix: multiple GblPoints at same arc-length"
                 << std::endl;
         throw std::overflow_error("Singular matrix inversion exception");
     }
     matWJ = matW * matJ;
     vecWd = matW * vecd;
 
 }

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

Retrieve global derivatives from a point.

Definition at line 323 of file GblPoint.cc.

References globalDerivatives.

                                                      {
     return globalDerivatives;
 }

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

Retrieve global derivatives labels from a point.

Definition at line 318 of file GblPoint.cc.

References globalLabels.

                                                {
     return globalLabels;
 }

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

Retrieve label of point.

Definition at line 336 of file GblPoint.cc.

References theLabel.

                                       {
     return theLabel;
 }

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

Retrieve local derivatives from a point.

Definition at line 288 of file GblPoint.cc.

References localDerivatives.

                                                     {
     return localDerivatives;
 }

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

Get measurement transformation (from diagonalization).

Parameters

[out] aTransformation Transformation matrix

Definition at line 171 of file GblPoint.cc.

References measDim, measTransformation, and transFlag.

                                                                     {
     aTransformation.ResizeTo(measDim, measDim);
     if (transFlag) {
         aTransformation = measTransformation;
     } else {
         aTransformation.UnitMatrix();
     }
 }

void gbl::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 160 of file GblPoint.cc.

References measPrecision, measProjection, and measResiduals.

                                     {
     aProjection = measProjection;
     aResiduals = measResiduals;
     aPrecision = measPrecision;
 }

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

Retrieve number of global derivatives from a point.

Definition at line 313 of file GblPoint.cc.

References globalDerivatives.

                                            {
     return globalDerivatives.GetNcols();
 }

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

Retrieve number of local derivatives from a point.

Definition at line 283 of file GblPoint.cc.

References localDerivatives.

                                           {
     return localDerivatives.GetNcols();
 }

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().

                               {
     return theOffset;
 }

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().

                                                 {
     return p2pJacobian;
 }

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

Retrieve scatterer of a point.

Parameters

[out]	aTransformation	Scatterer transformation from diagonalization
[out]	aResiduals	Scatterer residuals
[out]	aPrecision	Scatterer precision (diagonal)

Definition at line 242 of file GblPoint.cc.

References scatPrecision, scatResiduals, and scatTransformation.

                                     {
     aTransformation = scatTransformation;
     aResiduals = scatResiduals;
     aPrecision = scatPrecision;
 }

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

Get scatterer transformation (from diagonalization).

Parameters

[out] aTransformation Transformation matrix

Definition at line 253 of file GblPoint.cc.

References i, j, scatFlag, and scatTransformation.

                                                                     {
     aTransformation.ResizeTo(2, 2);
     if (scatFlag) {
         for (unsigned int i = 0; i < 2; ++i) {
             for (unsigned int j = 0; j < 2; ++j) {
                 aTransformation(i, j) = scatTransformation(i, j);
             }
         }
     } else {
         aTransformation.UnitMatrix();
     }
 }

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.

                                             {
     return measDim;
 }

bool gbl::GblPoint::hasScatterer ( ) const

Check for scatterer at a point.

Definition at line 232 of file GblPoint.cc.

References scatFlag.

                                   {
     return scatFlag;
 }

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

Print GblPoint.

Parameters

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

Definition at line 425 of file GblPoint.cc.

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

                                                   {
     std::cout << " GblPoint";
     if (theLabel) {
         std::cout << ", label " << theLabel;
         if (theOffset >= 0) {
             std::cout << ", offset " << theOffset;
         }
     }
     if (measDim) {
         std::cout << ", " << measDim << " measurements";
     }
     if (scatFlag) {
         std::cout << ", scatterer";
     }
     if (transFlag) {
         std::cout << ", diagonalized";
     }
     if (localDerivatives.GetNcols()) {
         std::cout << ", " << localDerivatives.GetNcols()
                 << " local derivatives";
     }
     if (globalDerivatives.GetNcols()) {
         std::cout << ", " << globalDerivatives.GetNcols()
                 << " global derivatives";
     }
     std::cout << std::endl;
     if (level > 0) {
         if (measDim) {
             std::cout << "  Measurement" << std::endl;
             std::cout << "   Projection: " << std::endl << measProjection
                     << std::endl;
             std::cout << "   Residuals: " << measResiduals << std::endl;
             std::cout << "   Precision: " << measPrecision << std::endl;
         }
         if (scatFlag) {
             std::cout << "  Scatterer" << std::endl;
             std::cout << "   Residuals: " << scatResiduals << std::endl;
             std::cout << "   Precision: " << scatPrecision << std::endl;
         }
         if (localDerivatives.GetNcols()) {
             std::cout << "  Local Derivatives:" << std::endl;
             localDerivatives.Print();
         }
         if (globalDerivatives.GetNcols()) {
             std::cout << "  Global Labels:";
             for (unsigned int i = 0; i < globalLabels.size(); ++i) {
                 std::cout << " " << globalLabels[i];
             }
             std::cout << std::endl;
             std::cout << "  Global Derivatives:" << std::endl;
             globalDerivatives.Print();
         }
         std::cout << "  Jacobian " << std::endl;
         std::cout << "   Point-to-point " << std::endl << p2pJacobian
                 << std::endl;
         if (theLabel) {
             std::cout << "   To previous offset " << std::endl << prevJacobian
                     << std::endl;
             std::cout << "   To next offset " << std::endl << nextJacobian
                     << std::endl;
         }
     }
 }

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

Define label of point (by GBLTrajectory constructor)

Parameters

[in] aLabel Label identifying point

Definition at line 331 of file GblPoint.cc.

References theLabel.

Referenced by SequenceTypes.DummyModule::__init__(), SequenceTypes.DummyBooleanModule::__init__(), and Vispa.Gui.FindDialog.FindDialog::reset().

                                            {
     theLabel = aLabel;
 }

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

Define offset for point (by GBLTrajectory constructor)

Parameters

[in] anOffset Offset number

Definition at line 344 of file GblPoint.cc.

References theOffset.

                                      {
     theOffset = anOffset;
 }

Member Data Documentation

TMatrixD gbl::GblPoint::globalDerivatives

private

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

Definition at line 107 of file GblPoint.h.

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

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

private

Labels of global (MP-II) derivatives.

Definition at line 106 of file GblPoint.h.

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

TMatrixD gbl::GblPoint::localDerivatives

private

Derivatives of measurement vs additional local (fit) parameters.

Definition at line 105 of file GblPoint.h.

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

unsigned int gbl::GblPoint::measDim

private

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

Definition at line 95 of file GblPoint.h.

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

SVector5 gbl::GblPoint::measPrecision

private

Measurement precision (diagonal of inverse covariance matrix)

Definition at line 98 of file GblPoint.h.

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

SMatrix55 gbl::GblPoint::measProjection

private

Projection from measurement to local system.

Definition at line 96 of file GblPoint.h.

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

SVector5 gbl::GblPoint::measResiduals

private

Measurement residuals.

Definition at line 97 of file GblPoint.h.

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

TMatrixD gbl::GblPoint::measTransformation

private

Transformation of diagonalization (of meas. precision matrix)

Definition at line 100 of file GblPoint.h.

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

SMatrix55 gbl::GblPoint::nextJacobian

private

Jacobian to next scatterer (or last measurement)

Definition at line 94 of file GblPoint.h.

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

SMatrix55 gbl::GblPoint::p2pJacobian

private

Point-to-point jacobian from previous point.

Definition at line 92 of file GblPoint.h.

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

SMatrix55 gbl::GblPoint::prevJacobian

private

Jacobian to previous scatterer (or first measurement)

Definition at line 93 of file GblPoint.h.

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

bool gbl::GblPoint::scatFlag

private

Scatterer present?

Definition at line 101 of file GblPoint.h.

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

SVector2 gbl::GblPoint::scatPrecision

private

Scattering precision (diagonal of inverse covariance matrix)

Definition at line 104 of file GblPoint.h.

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

SVector2 gbl::GblPoint::scatResiduals

private

Scattering residuals (initial kinks if iterating)

Definition at line 103 of file GblPoint.h.

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

SMatrix22 gbl::GblPoint::scatTransformation

private

Transformation of diagonalization (of scat. precision matrix)

Definition at line 102 of file GblPoint.h.

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

unsigned int gbl::GblPoint::theLabel

private

Label identifying point.

Definition at line 90 of file GblPoint.h.

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

int gbl::GblPoint::theOffset

private

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

Definition at line 91 of file GblPoint.h.

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

bool gbl::GblPoint::transFlag

private

Transformation exists?

Definition at line 99 of file GblPoint.h.

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

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation