General Broken Lines Trajectory. More...

Inheritance diagram for gblpy.gblfit.GblTrajectory:

Public Member Functions
def	__init__ (self, hasCurv=True, aDim=[0, 1])
	Create new trajectory. More...

def	addPoint (self, point)
	Add point to trajectory. More...

def	getNumPoints (self)
	Get number of points on trajectory. More...

def	getNumbers (self)
	Get (parameter) number of trajectory. More...

def	addExternalSeed (self, aLabel, aSeed)
	Add external seed to trajectory. More...

def	printPoints (self)
	Print points of trajectory. More...

def	printData (self)
	Print data of trajectory. More...

def	getData (self)
	Get data of trajectory. More...

def	milleOut (self, aFile, doublePrec=False)
	Write (data blocks of) trajectory to MP (binary) file (as float or double values). More...

def	milleIn (self, aFile)
	Read (data blocks of) trajectory from MP (binary) file. More...

def	getResults (self, aLabel)
	Get results (corrections, covarinace matrix) at point in forward or backward direction. More...

def	getMeasResults (self, aLabel)
	Get residuals at point from measurement. More...

def	getScatResults (self, aLabel)
	Get residuals at point from scatterer. More...

def	getResidual (self, iData)
	Get residuals from data of trajectory. More...

def	fit (self, optionList="", aLabel=0)
	Perform fit of trajectory. More...

Private Member Functions
def	__getJacobian (self, aLabel)
	Get jacobian for transformation from fit to track parameters at point. More...

def	__getFitToLocalJacobian (self, aPoint, measDim, nJacobian=1)
	Get (part of) jacobian for transformation from (trajectory) fit to track parameters at point. More...

def	__getFitToKinkJacobian (self, aPoint)
	Get jacobian for transformation from (trajectory) fit to kink parameters at point. More...

def	__getFitToStepJacobian (self, aPoint)
	Get jacobian for transformation from (trajectory) fit to step parameters at point. More...

def	__getFitToKinkAndStepJacobian (self, aPoint)
	Get jacobian for transformation from (trajectory) fit to kink and step parameters at point. More...

def	__getResAndErr (self, aData, used=True)

Private Attributes
	__numPoints
	number of points on trajectory; int More...

	__numOffsets
	number of (points with) offsets on trajectory; int More...

	__numCurvature
	'curvature' is fit parameter (=1); int More...

	__numParameters
	number fit parameters; int More...

	__numLocals
	number of local parameters; int More...

	__externalPoint
	label of point with external seed; int More...

	__dimensions
	active components of offsets (both ([0,1]) or single ([0] or [1]); list(int) More...

	__points
	points on trajectory; list(GblPoint) More...

	__data
	data (blocks) of trajectory; list(GblData) More...

	__externalSeed
	external seed (for local, fit parameters); matrix(float) More...

	__measDataIndex
	mapping points to data blocks from measurements; list(int) More...

	__scatDataIndex
	mapping points to data blocks from scatterers; list(int) More...

	__skippedMeasLabel
	label of point with measurements skipped in fit (for unbiased residuals) More...

	__matrix
	Calculate Jacobians to previous/next scatterer from point to point ones. More...

	__vector

Detailed Description

General Broken Lines Trajectory.

Definition at line 768 of file gblfit.py.

Constructor & Destructor Documentation

◆ init()

def gblpy.gblfit.GblTrajectory.__init__	(	self,
		hasCurv = `True`,
		aDim = `[0, 1]`
	)

Create new trajectory.

Parameters

hasCurv	flag for curvature; bool
aDim	active offset components (0:u_1, 1:u_2); list(int)

Definition at line 775 of file gblfit.py.

Member Function Documentation

◆ __getFitToKinkAndStepJacobian()

def gblpy.gblfit.GblTrajectory.__getFitToKinkAndStepJacobian	(	self,
		aPoint
	)

private

Get jacobian for transformation from (trajectory) fit to kink and step parameters at point.

Jacobian broken lines (q/p,..,u_i-1,u_i,u_i+1..) to kink (du') and step (du) parameters.

Thick sactterers only.

Parameters

aPoint point to use; GblPoint

Returns: labels for fit parameters with non zero derivatives, corresponding transformation matrix; list(vector(int), matrix(float))

Definition at line 1112 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__numCurvature, and gblpy.gblfit.GblTrajectory.__numLocals.

Referenced by gblpy.gblfit.GblTrajectory.fit().

◆ __getFitToKinkJacobian()

def gblpy.gblfit.GblTrajectory.__getFitToKinkJacobian	(	self,
		aPoint
	)

private

Get jacobian for transformation from (trajectory) fit to kink parameters at point.

Jacobian broken lines (q/p,..,u_i-1,u_i,u_i+1..) to kink (du') parameters.

Parameters

aPoint point to use; GblPoint

Returns: labels for fit parameters with non zero derivatives, corresponding transformation matrix; list(vector(int), matrix(float))

Definition at line 1042 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__numCurvature, and gblpy.gblfit.GblTrajectory.__numLocals.

Referenced by gblpy.gblfit.GblTrajectory.fit().

◆ __getFitToLocalJacobian()

def gblpy.gblfit.GblTrajectory.__getFitToLocalJacobian	(	self,
		aPoint,
		measDim,
		nJacobian = `1`
	)

private

Get (part of) jacobian for transformation from (trajectory) fit to track parameters at point.

Jacobian broken lines (q/p,..,u_i,u_i+1..) to local (q/p,u',u) parameters.

Parameters

aPoint	point to use; GblPoint
measDim	dimension of 'measurement' (2, 4 or 5); int
nJacobian	direction (0: to previous offset, 1: to next offset); int

Returns: labels for fit parameters with non zero derivatives, corresponding transformation matrix; list(vector(int), matrix(float))

Definition at line 952 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__numCurvature, and gblpy.gblfit.GblTrajectory.__numLocals.

Referenced by gblpy.gblfit.GblTrajectory.__getJacobian(), and gblpy.gblfit.GblTrajectory.fit().

◆ __getFitToStepJacobian()

def gblpy.gblfit.GblTrajectory.__getFitToStepJacobian	(	self,
		aPoint
	)

private

Get jacobian for transformation from (trajectory) fit to step parameters at point.

Jacobian broken lines (q/p,..,u_i-1,u_i,u_i+1..) to kink (du') and step (du) parameters.

Thick sactterers only.

Parameters

aPoint point to use; GblPoint

Returns: labels for fit parameters with non zero derivatives, corresponding transformation matrix; list(vector(int), matrix(float))

Definition at line 1081 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__numCurvature, and gblpy.gblfit.GblTrajectory.__numLocals.

Referenced by gblpy.gblfit.GblTrajectory.fit().

◆ __getJacobian()

def gblpy.gblfit.GblTrajectory.__getJacobian	(	self,
		aLabel
	)

private

Get jacobian for transformation from fit to track parameters at point.

Parameters

aLabel (signed) label of point; int

Returns: labels of required fit parameters and jacobian; list

Definition at line 901 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__getFitToLocalJacobian(), gblpy.gblfit.GblTrajectory.__numCurvature, gblpy.gblfit.GblTrajectory.__numLocals, gblpy.gblfit.GblTrajectory.__numPoints, and gblpy.gblfit.GblTrajectory.__points.

Referenced by gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.getResults().

◆ __getResAndErr()

def gblpy.gblfit.GblTrajectory.__getResAndErr	(	self,
		aData,
		used = `True`
	)

private

Definition at line 1145 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__data, gblpy.gblfit.GblTrajectory.__matrix, and gblpy.gblfit.GblTrajectory.getResidual().

Referenced by gblpy.gblfit.GblTrajectory.getMeasResults(), gblpy.gblfit.GblTrajectory.getResidual(), and gblpy.gblfit.GblTrajectory.getScatResults().

◆ addExternalSeed()

def gblpy.gblfit.GblTrajectory.addExternalSeed	(	self,
		aLabel,
		aSeed
	)

Add external seed to trajectory.

Parameters

aLabel	label of point with external seed; int
aSeed	seed (precision matrix of track parameters at point); matrix(float)

Definition at line 836 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__externalPoint, and gblpy.gblfit.GblTrajectory.__externalSeed.

◆ addPoint()

def gblpy.gblfit.GblTrajectory.addPoint	(	self,
		point
	)

Add point to trajectory.

Points have to be ordered in arc length.

Parameters

point point to be added; GblPoint

Returns: label of added point (1..number(points)); int

Definition at line 808 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__numLocals, gblpy.gblfit.GblTrajectory.__numPoints, and gblpy.gblfit.GblTrajectory.__points.

◆ fit()

def gblpy.gblfit.GblTrajectory.fit	(	self,
		optionList = `""`,
		aLabel = `0`
	)

Perform fit of trajectory.

Parameters

optionList	M-estimators to be used (one iteration per character); string
aLabel	label of point where to skip measurements (for unbiased residuals)

Returns: Chi2, Ndf, loss of weight from fit ([0., -1, 0.] if fit failed); list

Definition at line 1233 of file gblfit.py.

◆ getData()

def gblpy.gblfit.GblTrajectory.getData ( self )

Get data of trajectory.

Returns: data blocks; list

Definition at line 856 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__data.

◆ getMeasResults()

def gblpy.gblfit.GblTrajectory.getMeasResults	(	self,
		aLabel
	)

Get residuals at point from measurement.

Get (diagonalized) residual, error of measurement and residual and down-weighting factor for measurement at point

Parameters

aLabel label of point; int

Returns: vectors with residuals, errors of maeansurements and residuals, down-weighting factors; list

Definition at line 1185 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__getResAndErr(), gblpy.gblfit.GblTrajectory.__measDataIndex, and gblpy.gblfit.GblTrajectory.__skippedMeasLabel.

◆ getNumbers()

def gblpy.gblfit.GblTrajectory.getNumbers ( self )

Get (parameter) number of trajectory.

Returns: number of (fit) parameters; int

Definition at line 828 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__numCurvature, gblpy.gblfit.GblTrajectory.__numLocals, and gblpy.gblfit.GblTrajectory.__numParameters.

◆ getNumPoints()

def gblpy.gblfit.GblTrajectory.getNumPoints ( self )

Get number of points on trajectory.

Returns: number of points; int

Definition at line 821 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__numPoints.

◆ getResidual()

def gblpy.gblfit.GblTrajectory.getResidual	(	self,
		iData
	)

Get residuals from data of trajectory.

Parameters

iData index of data block; int

Definition at line 1224 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__getResAndErr().

Referenced by gblpy.gblfit.GblTrajectory.__getResAndErr().

◆ getResults()

def gblpy.gblfit.GblTrajectory.getResults	(	self,
		aLabel
	)

Get results (corrections, covarinace matrix) at point in forward or backward direction.

The point is identified by its label (1..number(points)), the sign distinguishes the backward (facing previous point) and forward 'side' (facing next point). For scatterers the track direction may change in between.

Parameters

aLabel signed label of point (<0 backward, >0 forward); int

Returns: correction vector, covarinace matrix for track parameters; list

Definition at line 1167 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.__matrix, and gblpy.gblfit.GblTrajectory.__vector.

◆ getScatResults()

def gblpy.gblfit.GblTrajectory.getScatResults	(	self,
		aLabel
	)

Get residuals at point from scatterer.

Get (diagonalized) residual, error of measurement and residual and down-weighting factor for scatterering kinks at point

Parameters

aLabel label of point; int

Returns: vectors with residuals, errors of maeansurements and residuals, down-weighting factors; list

Definition at line 1206 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__getResAndErr(), and gblpy.gblfit.GblTrajectory.__scatDataIndex.

◆ milleIn()

def gblpy.gblfit.GblTrajectory.milleIn	(	self,
		aFile
	)

Read (data blocks of) trajectory from MP (binary) file.

Parameters

aFile MP file

Definition at line 876 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__data, gblpy.gblfit.GblTrajectory.__dimensions, gblpy.gblfit.GblTrajectory.__numCurvature, gblpy.gblfit.GblTrajectory.__numLocals, and gblpy.gblfit.GblTrajectory.__numParameters.

◆ milleOut()

def gblpy.gblfit.GblTrajectory.milleOut	(	self,
		aFile,
		doublePrec = `False`
	)

Write (data blocks of) trajectory to MP (binary) file (as float or double values).

Parameters

aFile	MP file
doublePrec	flag for storage in as double values

Definition at line 864 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__data.

◆ printData()

def gblpy.gblfit.GblTrajectory.printData ( self )

Print data of trajectory.

Definition at line 847 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__data.

◆ printPoints()

def gblpy.gblfit.GblTrajectory.printPoints ( self )

Print points of trajectory.

Definition at line 841 of file gblfit.py.

References gblpy.gblfit.GblTrajectory.__points.

Member Data Documentation

◆ __data

gblpy.gblfit.GblTrajectory.__data

private

data (blocks) of trajectory; list(GblData)

Down weight (data) outliers.

Parameters

aMethod M-estimator; int

Returns: loss of weight (sum(1-down_weighting)); float

Calculate predictions.

Definition at line 793 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getResAndErr(), gblpy.gblfit.GblTrajectory.fit(), gblpy.gblfit.GblTrajectory.getData(), gblpy.gblfit.GblTrajectory.milleIn(), gblpy.gblfit.GblTrajectory.milleOut(), and gblpy.gblfit.GblTrajectory.printData().

◆ __dimensions

gblpy.gblfit.GblTrajectory.__dimensions

private

active components of offsets (both ([0,1]) or single ([0] or [1]); list(int)

Definition at line 789 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getFitToKinkAndStepJacobian(), gblpy.gblfit.GblTrajectory.__getFitToKinkJacobian(), gblpy.gblfit.GblTrajectory.__getFitToLocalJacobian(), gblpy.gblfit.GblTrajectory.__getFitToStepJacobian(), gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.milleIn().

◆ __externalPoint

gblpy.gblfit.GblTrajectory.__externalPoint

private

label of point with external seed; int

Definition at line 787 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.addExternalSeed(), and gblpy.gblfit.GblTrajectory.fit().

◆ __externalSeed

gblpy.gblfit.GblTrajectory.__externalSeed

private

external seed (for local, fit parameters); matrix(float)

Definition at line 795 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.addExternalSeed(), and gblpy.gblfit.GblTrajectory.fit().

◆ __matrix

gblpy.gblfit.GblTrajectory.__matrix

private

Calculate Jacobians to previous/next scatterer from point to point ones.

Prepare fit; generate data from points.

Build linear equation system from data.

Definition at line 1375 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getResAndErr(), and gblpy.gblfit.GblTrajectory.getResults().

◆ __measDataIndex

gblpy.gblfit.GblTrajectory.__measDataIndex

private

mapping points to data blocks from measurements; list(int)

Definition at line 797 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.getMeasResults().

◆ __numCurvature

gblpy.gblfit.GblTrajectory.__numCurvature

private

'curvature' is fit parameter (=1); int

Definition at line 781 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getFitToKinkAndStepJacobian(), gblpy.gblfit.GblTrajectory.__getFitToKinkJacobian(), gblpy.gblfit.GblTrajectory.__getFitToLocalJacobian(), gblpy.gblfit.GblTrajectory.__getFitToStepJacobian(), gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.fit(), gblpy.gblfit.GblTrajectory.getNumbers(), and gblpy.gblfit.GblTrajectory.milleIn().

◆ __numLocals

gblpy.gblfit.GblTrajectory.__numLocals

private

number of local parameters; int

Definition at line 785 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getFitToKinkAndStepJacobian(), gblpy.gblfit.GblTrajectory.__getFitToKinkJacobian(), gblpy.gblfit.GblTrajectory.__getFitToLocalJacobian(), gblpy.gblfit.GblTrajectory.__getFitToStepJacobian(), gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.addPoint(), gblpy.gblfit.GblTrajectory.fit(), gblpy.gblfit.GblTrajectory.getNumbers(), and gblpy.gblfit.GblTrajectory.milleIn().

◆ __numOffsets

gblpy.gblfit.GblTrajectory.__numOffsets

private

number of (points with) offsets on trajectory; int

Define offsets from list of points.

Definition at line 779 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.fit().

◆ __numParameters

gblpy.gblfit.GblTrajectory.__numParameters

private

number fit parameters; int

Definition at line 783 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.fit(), gblpy.gblfit.GblTrajectory.getNumbers(), and gblpy.gblfit.GblTrajectory.milleIn().

◆ __numPoints

gblpy.gblfit.GblTrajectory.__numPoints

private

number of points on trajectory; int

Definition at line 777 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.addPoint(), gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.getNumPoints().

◆ __points

gblpy.gblfit.GblTrajectory.__points

private

points on trajectory; list(GblPoint)

Definition at line 791 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.__getJacobian(), gblpy.gblfit.GblTrajectory.addPoint(), gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.printPoints().

◆ __scatDataIndex

gblpy.gblfit.GblTrajectory.__scatDataIndex

private

mapping points to data blocks from scatterers; list(int)

Definition at line 799 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.fit(), and gblpy.gblfit.GblTrajectory.getScatResults().

◆ __skippedMeasLabel

gblpy.gblfit.GblTrajectory.__skippedMeasLabel

private

label of point with measurements skipped in fit (for unbiased residuals)

Definition at line 801 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.getMeasResults().

◆ __vector

gblpy.gblfit.GblTrajectory.__vector

private

Definition at line 1376 of file gblfit.py.

Referenced by gblpy.gblfit.GblTrajectory.getResults().

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

gblfit.py

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ __init__()

Member Function Documentation

◆ __getFitToKinkAndStepJacobian()

◆ __getFitToKinkJacobian()

◆ __getFitToLocalJacobian()

◆ __getFitToStepJacobian()

◆ __getJacobian()

◆ __getResAndErr()

◆ addExternalSeed()

◆ addPoint()

◆ fit()

◆ getData()

◆ getMeasResults()

◆ getNumbers()

◆ getNumPoints()

◆ getResidual()

◆ getResults()

◆ getScatResults()

◆ milleIn()

◆ milleOut()

◆ printData()

◆ printPoints()

Member Data Documentation

◆ __data

◆ __dimensions

◆ __externalPoint

◆ __externalSeed

◆ __matrix

◆ __measDataIndex

◆ __numCurvature

◆ __numLocals

◆ __numOffsets

◆ __numParameters

◆ __numPoints

◆ __points

◆ __scatDataIndex

◆ __skippedMeasLabel

◆ __vector

◆ init()