GeneralBrokenLines V03-00-00
gblpy
Classes | Functions
gblpy.gbltst Namespace Reference

Classes

class  gblCurviSystem
 Curvilinear system as local system. More...
 
class  gblMeasSystem
 Measurement system as local system. More...
 

Functions

def example1 ()
 Create points on initial trajectory, create trajectory from points, fit and write trajectory to MP-II binary file, get track parameter corrections and covariance matrix at points. More...
 
def example2 ()
 Read trajectory from MP-II binary file and refit. More...
 
def gblSimpleJacobian (ds, cosl, bfac)
 Simple jacobian. More...
 

Detailed Description

Simple Test Program for General Broken Lines.

Created on Jul 27, 2011

@author: kleinwrt

Function Documentation

◆ example1()

def gblpy.gbltst.example1 ( )

Create points on initial trajectory, create trajectory from points, fit and write trajectory to MP-II binary file, get track parameter corrections and covariance matrix at points.

Equidistant measurement layers and thin scatterers, propagation with simple jacobian (quadratic in arc length differences). Curvilinear system (T,U,V) or measurement system (I,J,K) as local coordinate system and (q/p, slopes, offsets) as local track parameters.

This example simulates and refits tracks in a system of planar detectors with 2D measurements in a constant magnet field in Z direction using the curvilinear system as local system. The true track parameters are randomly smeared with respect to a (constant and straight) reference trajectory with direction (lambda, phi) and are used (only) for the on-the-fly simulation of the measurements and scatterers. The predictions from the reference trajectory are therefore always zero and the residuals needed (by addMeasurement) are equal to the measurements.

Definition at line 55 of file gbltst.py.

References gblpy.gbltst.gblSimpleJacobian().

◆ example2()

def gblpy.gbltst.example2 ( )

Read trajectory from MP-II binary file and refit.

Definition at line 281 of file gbltst.py.

Referenced by gblpy.gbltst.gblCurviSystem.getScatPrecision().

◆ gblSimpleJacobian()

def gblpy.gbltst.gblSimpleJacobian (   ds,
  cosl,
  bfac 
)

Simple jacobian.

Simple jacobian for (q/p, slopes, offsets) in curvilinear system quadratic in arc length difference.

Parameters
dsarc length difference; float
coslcos(lambda); float
bfacBz*c; float
Returns
jacobian to move by 'ds' on trajectory matrix(float)

Definition at line 326 of file gbltst.py.

Referenced by gblpy.gbltst.example1().