GeneralBrokenLines V03-01-02
gblpy
|
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... | |
Simple Test Program for General Broken Lines. Created on Jul 27, 2011 @author: kleinwrt
def 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 gblSimpleJacobian().
def gbltst.example2 | ( | ) |
Read trajectory from MP-II binary file and refit.
Definition at line 281 of file gbltst.py.
Referenced by gbltst.gblCurviSystem.getScatPrecision().
def gbltst.gblSimpleJacobian | ( | ds, | |
cosl, | |||
bfac | |||
) |
Simple jacobian.
Simple jacobian for (q/p, slopes, offsets) in curvilinear system quadratic in arc length difference.
ds | arc length difference; float |
cosl | cos(lambda); float |
bfac | Bz*c; float |
Definition at line 328 of file gbltst.py.
Referenced by example1().