General Broken Lines are extending the scope of the original broken lines from a tracking setup with a solenoidal magnetic field, independent tracking in and perpendicular to the bending plane and multiple scattering in approximated thick scatterers to arbitrary propagation of the trajectory, arbitrary measurement planes and thick scatterers described by an equivalent pair of thin scatterers.
A trajectory based on General Broken Lines is a track refit to add the description of multiple scattering to an initial trajectory based on the propagation in a magnetic field (and average energy loss). It is constructed from a sequence of (pairs of) thin scatterers describing the multiple scattering in the material between adjacent measurement planes. Predictions for the measurements are obtained by interpolation between the enclosing scatterers and triplets of adjacent scatterers define kink angles with variance according to the material of the central scatterer. This requires the propagation (on the initial trajectory) from a measurement plane or scatterer to the previous and next scatterer. The propagation has to be provided by the user in form of the jacobians for the local track parameters (transformation matrices). In the case of an homogeneous magnetic field they can be calculated from NIM A566:687-698,2006 by A. Strandlie and W. Wittek. The initial trajectory can be the result of a fit of the measurements (internal seed) or a prediction from another detector part (external seed). A General Broken Lines fit with external seed and one (additional) measurement is equivalent to the filtering step of the track fit with a Kalman filter.
As the track refit can provide the complete covariance matrix of all track parameters General Broken Lines are well suited as track model for calibration and alignment with Millepede II.
Implementations in Fortran, Python and C++ are provided by DESY under the terms of the LGPLv2 license and are maintained in a publicly accessible subversion repository at DESY that can be browsed via viewvc or WebSVN.
The source code is publicly available via the subversion command line:
svn checkout http://svnsrv.desy.de/public/GeneralBrokenLines/tags/V02-03-01 GeneralBrokenLines
The further development is concentrating on the C++ version.
The usage of the software package is described in corresponding README files. The Makefile produces a test executable demonstrating the usage with a simplified track model. All versions can directly output the trajectory information to Millepede-II binary files.
Online documentation (doxygen) is available for:
The C++ versions before V02-00-00 have worked with ROOT (root.cern.ch) as linear algebra package. Since version V02-00-00 the EIGEN template library (eigen.tuxfamily.org) is in use.
Fixes for the geometric constraint with external measurement and for internal iterations (setZero for linear equation system added). The GblPoint constructor can optionally reserve now space for several measurements.
Multiple measurements can now be added to a point to implement ambiguities.
Minor fixes and updates.
For composed trajectories added geometric constraint (for common vertex) to full kinematic constraint.
Cleanup of the C++ examples.
The CMake file for the examples has been fixed.
Examples have been reworked and extended. A problem for a composed trajectory with external measurement in combination with additional local parameters has been fixed. Treatment of down-weighting for (un)biased residuals improved.
Another problem (in V02) with sometimes uninitialised values (in composed trajectories) has been fixed.
A problem (in V02) with sometimes uninitialised values has been fixed.
The χ2 calculation for M-estimators has been fixed (usage of proper objective functions).
Several GBL methods implemented as templates to allow more EIGEN compile time optimization.
Replaced ROOT by Eigen3 as linear algebra package. User input and output with ROOT matrices is still possible with the preprocessor directive
Unbiased residuals implemented for Python and C++.
Authorship, copyright and license defined, cleanups and C++11 compatibility.
Bug fixes and fixes for (C++) code analysis by covertity (coverity.com), diagonalization of scatterers implemented.
Adaptations for use by ILCSOFT (ilcsoft.desy.de).
First "complete" version building with CMAKE.
For questions, suggestions or help please contact: Claus(dot)Kleinwort(at)desy(dot)de