Using GBL with a step-wise helix. More...

#include <StepWiseHelixGeneralBrokenLineInterfaceProcessor.h>

Inheritance diagram for marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor:

Public Member Functions
virtual Processor *	newProcessor ()

	StepWiseHelixGeneralBrokenLineInterfaceProcessor ()

virtual void	init ()
	Initialize processor. More...

virtual void	processRunHeader (EVENT::LCRunHeader *run)

virtual void	processEvent (EVENT::LCEvent *evt)
	Process event. More...

virtual void	check (EVENT::LCEvent *evt)

virtual void	end ()

Protected Attributes
std::string	_inputCollectionName
	Name of the input collection – track seeds. More...

std::string	_outputTrackCollectionName
	Name of the output track collection. More...

double	_bfieldScaleFactor
	scale factor for magnetic field (default: 1.0) More...

bool	_outputIsPersistent
	whether the output is to be stored or not (default: true) More...

std::string	_fitOptions
	list of iterations with 'h' Huber, 't' Tukey or 'c' Cauchy down-weighting More...

bool	_writeMillepedeOut
	selects output to Millepede-II binary file More...

std::string	_fileNameMillepedeFile
	name of Millepede-II binary file More...

int	_milleCalcMethod
	defines which parameter should be calculated for Millepede-II binary file More...

int	_milleMinHits
	defines minimum number of hits on track for output to Millepede binary file More...

bool	_encodedModuleID
	Module ID is encoded in CellID0. More...

double	_maxStep
	maximal step size (default: 10.) More...

double	_x0PerUnitLength
	X0 per unit length (TPC gas as homogeneous scatterer) (default: 0.) More...

bool	_thickScatterer
	selects thick (otherwise thin) scatterers for multiple scattering (default: false) More...

double	_defaultMomentum
	default momentum for multiple scattering (for Bfield off, GeV) (default: 10.) More...

Private Member Functions
void	calcLineSeed (const std::vector< TrackerHit >, const double , double &, double &, double &, double &, double &) const
	Simple seed for straight line. More...

Private Attributes
gbl::GblTrajectory *	_trajectory
	GBL trajectory. More...

gbl::MilleBinary *	_milleBinary
	Millepede-II binary file. More...

bool	_curvature
	flag for curved track (helix, else straight line) More...

std::vector< unsigned int >	_theLabels
	Labels of (global) MP-II parameters. More...

TVectorD *	_correctionVector
	correction vector from GBL fit (for track parameters) More...

TMatrixDSym *	_covarianceMatrix
	covariance matrix from GBL More...

int	_ndf
	number of degrees of freedom in GBL fit More...

double	_chisquare
	chi2 from GBL fit More...

double	_lostweight
	weight lost by down-weighting More...

double	_refPointS
	arc-length of reference point More...

unsigned int	_refPointIndex
	label of reference point More...

double	_Xcenter
	TPC center, X coordinate. More...

double	_Ycenter
	TPC center, Y coordinate. More...

Detailed Description

Using GBL with a step-wise helix.

Processor interface to step-wise helix (or line) fitting with General Broken Lines, written by Claus Kleinwort

Author: C. Kleinwort, DESY

This processor interfaces the re-fitting functionality of the general broken line package by Claus Kleinwort. For information on the General Broken Lines, please refer to the publication and/or the wiki pages:

C. Kleinwort, General Broken Lines as advanced track fitting method, NIM A, 673 (2012), 107-110, doi:10.1016/j.nima.2012.01.024
http://www.wiki.terascale.de/index.php/GeneralBrokenLines

The standard behaviour is that a single, straightforward fit without any outlier rejection is performed. For multiple turns, including some outlier rejection; extend the fit options in the according steering parameter. The calculation of the measurement (pad) direction in XY uses "module.getLocalPadLayout().getCoordinateType()" to select polar or cartesian pad row geometry.

Inputs and requirements

Input to the processor are the found tracks with their associated hits (with the reference point inside the magnetic field!).
The curvilinear system is used as local system at each hit.
The curvilinear track parameters are used as local parameters at each hit.
A form of hit cleaning should be performed before calling the fit – before using the fit results.

Output

Output are the fitted tracks with the best estimation of track parameters and their uncertainties.
The specific energy loss is estimated from the median (and the MAD) of the 1/sqrt(E_scaled) distribution.
Optional output is a Millepede-II binary file, which can be inspected

Parameters

"InputSeedTracks":string	The name of the input collection of track candidates (default: "TPCSeedTracks")
"OutputTracks":string	The name of the output collection with the fitted tracks(default: "TPCTracks")
"WriteOutputToStorage":bool	To set whether the output should be written to storage
"BFieldScaleFactor":double	`Optional` parameter, scales magnetic field (map), use 1.0 (default) for field ON or 0.0 for field OFF
"FitOptions":string	`Optional` parameter, defines if multiple passes with "h" huber, "t" tukey or "c" cauchy downweighting should be performed
"WriteMillepedeBinary":bool	`Optional` parameter, set to "true" if you want a Millepede binary file (default:"false")
"MillepedeFilename":string	`Optional` parameter to set the name of the Millepede binary file (default:"")
"MillePedeCalcMethod":int	`Optional` parameter to define set of (global) parameters (default:0 = module alignment)
"MillePedeMinHits":int	`Optional` parameter to define the minimum number of hits on track for output to Millepede binary file (default: 42)
"EncodedModuleID":bool	`Optional` parameter, flag for encoding of module ID in CellID0 (default true)
"MaxStep":double	`Optional` parameter, maximal step size (with constant magnetic field) (default 10.0)
"X0PerUnitLength":bool	`Optional` parameter, X0 per unit length (TPC gas as homogeneous scatterer) (default: 0.)
"ThickScatterer":bool	`Optional` parameter, selects thick (otherwise thin) scatterers for multiple scattering (default: false)
"DefaultMomentum":bool	`Optional` parameter, default momentum for multiple scattering (for Bfield off, GeV) (default: 10.)

WORK IN PROGRESS: currently only works for modules with concentric pad rows (using center from GEAR)

Todo:

implement additional output of pull distributions, residuals, track paramters for each measurement point

A detailed description is available as lcnote LC-TOOL-2014-008. Changes with respect to this note:

Additional parameter "MillePedeMinHits" for help class "trackerAlcaSelector"
160223 CHK: Implement TPC gas as homogeneous scatterer (between first and last hit) modelled by thin or thick scatterers (per hit)

Definition at line 82 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Constructor & Destructor Documentation

◆ StepWiseHelixGeneralBrokenLineInterfaceProcessor()

marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::StepWiseHelixGeneralBrokenLineInterfaceProcessor ( )

Definition at line 47 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

References _bfieldScaleFactor, _defaultMomentum, _encodedModuleID, _fileNameMillepedeFile, _fitOptions, _inputCollectionName, _maxStep, _milleCalcMethod, _milleMinHits, _outputIsPersistent, _outputTrackCollectionName, _thickScatterer, _writeMillepedeOut, and _x0PerUnitLength.

Referenced by newProcessor().

Member Function Documentation

◆ calcLineSeed()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::calcLineSeed	(	const std::vector< TrackerHit *>	hitList,
		const double *	refPos,
		double &	omega,
		double &	phi,
		double &	d0,
		double &	tanLambda,
		double &	z0
	)		const

private

Simple seed for straight line.

Calculate seed from straight line defined by (radially) first and last hit

Parameters

[in]	hitList	list of hits
[in]	refPos	reference position
[out]	omega	curvature, = 0.
[out]	phi	XY slope
[out]	d0	XY offset
[out]	tanLambda	ZS slope
[out]	z0	ZS offset

Definition at line 412 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

Referenced by processEvent().

◆ check()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::check ( EVENT::LCEvent * evt )

virtual

Definition at line 392 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

Referenced by newProcessor().

◆ end()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::end ( )

virtual

Definition at line 395 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

References _milleBinary.

Referenced by newProcessor().

◆ init()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::init ( )

virtual

Initialize processor.

Definition at line 97 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

References _fileNameMillepedeFile, _milleBinary, _writeMillepedeOut, _Xcenter, and _Ycenter.

Referenced by newProcessor().

◆ newProcessor()

virtual Processor* marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::newProcessor ( )

inlinevirtual

Definition at line 85 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

References check(), end(), init(), processEvent(), processRunHeader(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ processEvent()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::processEvent ( EVENT::LCEvent * evt )

virtual

Process event.

Order hits by arc-length

Multiple scattering

Create vector of GBL points. from helper hits.

Create trajectory from list of GBL points.

Fit trajectory.

Fill output collection.

Definition at line 113 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

Referenced by newProcessor().

◆ processRunHeader()

void marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::processRunHeader ( EVENT::LCRunHeader * run )

virtual

Definition at line 109 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.cc.

Referenced by newProcessor().

Member Data Documentation

◆ _bfieldScaleFactor

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_bfieldScaleFactor

protected

scale factor for magnetic field (default: 1.0)

Definition at line 103 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _chisquare

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_chisquare

private

chi2 from GBL fit

Definition at line 124 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _correctionVector

TVectorD* marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_correctionVector

private

correction vector from GBL fit (for track parameters)

Definition at line 121 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

◆ _covarianceMatrix

TMatrixDSym* marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_covarianceMatrix

private

covariance matrix from GBL

Definition at line 122 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

◆ _curvature

bool marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_curvature

private

flag for curved track (helix, else straight line)

Definition at line 119 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _defaultMomentum

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_defaultMomentum

protected

default momentum for multiple scattering (for Bfield off, GeV) (default: 10.)

Definition at line 114 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _encodedModuleID

bool marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_encodedModuleID

protected

Module ID is encoded in CellID0.

Definition at line 110 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _fileNameMillepedeFile

std::string marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_fileNameMillepedeFile

protected

name of Millepede-II binary file

Definition at line 107 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by init(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _fitOptions

std::string marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_fitOptions

protected

list of iterations with 'h' Huber, 't' Tukey or 'c' Cauchy down-weighting

Definition at line 105 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _inputCollectionName

std::string marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_inputCollectionName

protected

Name of the input collection – track seeds.

Definition at line 101 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _lostweight

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_lostweight

private

weight lost by down-weighting

Definition at line 125 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _maxStep

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_maxStep

protected

maximal step size (default: 10.)

Definition at line 111 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _milleBinary

gbl::MilleBinary* marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_milleBinary

private

Millepede-II binary file.

Definition at line 118 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by end(), init(), and processEvent().

◆ _milleCalcMethod

int marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_milleCalcMethod

protected

defines which parameter should be calculated for Millepede-II binary file

Definition at line 108 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _milleMinHits

int marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_milleMinHits

protected

defines minimum number of hits on track for output to Millepede binary file

Definition at line 109 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _ndf

int marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_ndf

private

number of degrees of freedom in GBL fit

Definition at line 123 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _outputIsPersistent

bool marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_outputIsPersistent

protected

whether the output is to be stored or not (default: true)

Definition at line 104 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _outputTrackCollectionName

std::string marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_outputTrackCollectionName

protected

Name of the output track collection.

Definition at line 102 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _refPointIndex

unsigned int marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_refPointIndex

private

label of reference point

Definition at line 127 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _refPointS

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_refPointS

private

arc-length of reference point

Definition at line 126 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

◆ _theLabels

std::vector<unsigned int> marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_theLabels

private

Labels of (global) MP-II parameters.

Definition at line 120 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

◆ _thickScatterer

bool marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_thickScatterer

protected

selects thick (otherwise thin) scatterers for multiple scattering (default: false)

Definition at line 113 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _trajectory

gbl::GblTrajectory* marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_trajectory

private

GBL trajectory.

Definition at line 117 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent().

◆ _writeMillepedeOut

bool marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_writeMillepedeOut

protected

selects output to Millepede-II binary file

Definition at line 106 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by init(), processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _x0PerUnitLength

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_x0PerUnitLength

protected

X0 per unit length (TPC gas as homogeneous scatterer) (default: 0.)

Definition at line 112 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by processEvent(), and StepWiseHelixGeneralBrokenLineInterfaceProcessor().

◆ _Xcenter

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_Xcenter

private

TPC center, X coordinate.

Definition at line 128 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by init(), and processEvent().

◆ _Ycenter

double marlintpc::StepWiseHelixGeneralBrokenLineInterfaceProcessor::_Ycenter

private

TPC center, Y coordinate.

Definition at line 129 of file StepWiseHelixGeneralBrokenLineInterfaceProcessor.h.

Referenced by init(), and processEvent().

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

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Inputs and requirements

Output

Constructor & Destructor Documentation

◆ StepWiseHelixGeneralBrokenLineInterfaceProcessor()

Member Function Documentation

◆ calcLineSeed()

◆ check()

◆ end()

◆ init()

◆ newProcessor()

◆ processEvent()

◆ processRunHeader()

Member Data Documentation

◆ _bfieldScaleFactor

◆ _chisquare

◆ _correctionVector

◆ _covarianceMatrix

◆ _curvature

◆ _defaultMomentum

◆ _encodedModuleID

◆ _fileNameMillepedeFile

◆ _fitOptions

◆ _inputCollectionName

◆ _lostweight

◆ _maxStep

◆ _milleBinary

◆ _milleCalcMethod

◆ _milleMinHits

◆ _ndf

◆ _outputIsPersistent

◆ _outputTrackCollectionName

◆ _refPointIndex

◆ _refPointS

◆ _theLabels

◆ _thickScatterer

◆ _trajectory

◆ _writeMillepedeOut

◆ _x0PerUnitLength

◆ _Xcenter

◆ _Ycenter