KalFitTrack Class Reference

Description of a track class (<- Helix.cc). More...

#include <KalFitTrack.h>

Inheritance diagram for KalFitTrack:

Public Member Functions
	KalFitTrack (const HepPoint3D &pivot, const CLHEP::HepVector &a, const CLHEP::HepSymMatrix &Ea, unsigned int m, double chiSq, unsigned int nhits)
	constructor
	~KalFitTrack (void)
	destructor
double	intersect_cylinder (double r) const
	Intersection with different geometry.
double	intersect_yz_plane (const HepTransform3D &plane, double x) const
double	intersect_zx_plane (const HepTransform3D &plane, double y) const
double	intersect_xy_plane (double z) const
void	msgasmdc (double path, int index)
	Calculate multiple scattering angle.
void	ms (double path, const KalFitMaterial &m, int index)
void	eloss (double path, const KalFitMaterial &m, int index)
	Calculate total energy lost in material.
double	smoother_Mdc (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
	Kalman smoother for Mdc.
double	smoother_Mdc_csmalign (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
double	smoother_Mdc (KalFitHitMdc &HitMdc, CLHEP::Hep3Vector &meas, KalFitHelixSeg &seg, double &dchi2, int csmflag)
double	update_hits (KalFitHitMdc &HitMdc, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, double &dtrack, double &dtracknew, double &dtdc, int csmflag)
	Include the Mdc wire hits.
double	update_hits (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	update_hits_csmalign (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc)
void	update_last (void)
	Record the current parameters as ..._last information :
void	update_forMdc (void)
void	point_last (const HepPoint3D &point)
	set and give out the last point of the track
const HepPoint3D &	point_last (void)
const HepPoint3D &	pivot_last (void) const
	returns helix parameters
const CLHEP::HepVector &	a_last (void) const
const CLHEP::HepSymMatrix &	Ea_last (void) const
const HepPoint3D &	pivot_forMdc (void) const
const CLHEP::HepVector &	a_forMdc (void) const
const CLHEP::HepSymMatrix &	Ea_forMdc (void) const
void	path_add (double path)
	Update the path length estimation.
void	addPathSM (double path)
double	getPathSM (void)
void	addTofSM (double time)
double	getTofSM (void)
void	fiTerm (double fi)
double	getFiTerm (void)
void	tof (double path)
	Update the tof estimation.
double	filter (double v_m, const CLHEP::HepVector &m_H, double v_d, double m_V)
double	cor_tanldep (double *p, double er)
	Correct the error according the current tanl value :
void	update_bit (int i)
int	insist (void) const
	Extractor :
int	type (void) const
int	trasan_id (void) const
double	r0 (void) const
double	mass (void) const
double	chiSq (void) const
double	chiSq_back (void) const
int	ndf_back (void) const
double	pathip (void) const
double	path_rd (void) const
double	path_ab (void) const
double *	pathl (void)
CLHEP::Hep3Vector *	mom (void)
double	tof (void) const
double	tof_kaon (void) const
double	tof_proton (void) const
double	p_kaon (void) const
double	p_proton (void) const
double	dchi2_max (void) const
double	r_max (void) const
unsigned int	nchits (void) const
unsigned int	nster (void) const
unsigned int	ncath (void) const
int	pat1 (void) const
int	pat2 (void) const
int	nhit_r (void) const
int	nhit_z (void) const
void	type (int t)
	Reinitialize (modificator).
void	trasan_id (int t)
void	insist (int t)
void	pathip (double pl)
void	p_kaon (double pl)
void	p_proton (double pl)
void	chiSq (double c)
void	chiSq_back (double c)
void	ndf_back (int n)
void	nchits (int n)
void	nster (int n)
void	add_nhit_r (void)
void	add_nhit_z (void)
double	PathL (int layer)
	Function to calculate the path length in the layer.
void	appendHitsMdc (KalFitHitMdc h)
	Functions for Mdc hits list.
void	HitsMdc (vector< KalFitHitMdc > &lh)
vector< KalFitHitMdc > &	HitsMdc (void)
KalFitHitMdc &	HitMdc (int i)
void	appendHelixSegs (KalFitHelixSeg s)
void	HelixSegs (vector< KalFitHelixSeg > &vs)
vector< KalFitHelixSeg > &	HelixSegs (void)
KalFitHelixSeg &	HelixSeg (int i)
void	order_wirhit (int index)
void	order_hits (void)
void	number_wirhit (void)
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot)
	Sets pivot position in a given mag field.
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot, double &pathl)
double	radius_numf (void) const
	Estimation of the radius in a given mag field.
double	getSigma (int layerId, double driftDist) const
double	getSigma (KalFitHitMdc &hitmdc, double tanlam, int lr, double dist) const
double	getDriftDist (KalFitHitMdc &hitmdc, double drifttime, int lr) const
double	getDriftTime (KalFitHitMdc &hitmdc, double toftime) const
double	getT0 (void) const
HepSymMatrix	getInitMatrix (void) const
double	getDigi () const
void	chgmass (int i)
int	nLayerUsed ()
void	resetLayerUsed ()
void	useLayer (int iLay)
	KalFitTrack (const HepPoint3D &pivot, const CLHEP::HepVector &a, const CLHEP::HepSymMatrix &Ea, unsigned int m, double chiSq, unsigned int nhits)
	constructor
	~KalFitTrack (void)
	destructor
double	intersect_cylinder (double r) const
	Intersection with different geometry.
double	intersect_yz_plane (const HepTransform3D &plane, double x) const
double	intersect_zx_plane (const HepTransform3D &plane, double y) const
double	intersect_xy_plane (double z) const
void	msgasmdc (double path, int index)
	Calculate multiple scattering angle.
void	ms (double path, const KalFitMaterial &m, int index)
void	eloss (double path, const KalFitMaterial &m, int index)
	Calculate total energy lost in material.
double	smoother_Mdc (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
	Kalman smoother for Mdc.
double	smoother_Mdc_csmalign (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
double	smoother_Mdc (KalFitHitMdc &HitMdc, CLHEP::Hep3Vector &meas, KalFitHelixSeg &seg, double &dchi2, int csmflag)
double	update_hits (KalFitHitMdc &HitMdc, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, double &dtrack, double &dtracknew, double &dtdc, int csmflag)
	Include the Mdc wire hits.
double	update_hits (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	update_hits_csmalign (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc)
void	update_last (void)
	Record the current parameters as ..._last information :
void	update_forMdc (void)
void	point_last (const HepPoint3D &point)
	set and give out the last point of the track
const HepPoint3D &	point_last (void)
const HepPoint3D &	pivot_last (void) const
	returns helix parameters
const CLHEP::HepVector &	a_last (void) const
const CLHEP::HepSymMatrix &	Ea_last (void) const
const HepPoint3D &	pivot_forMdc (void) const
const CLHEP::HepVector &	a_forMdc (void) const
const CLHEP::HepSymMatrix &	Ea_forMdc (void) const
void	path_add (double path)
	Update the path length estimation.
void	addPathSM (double path)
double	getPathSM (void)
void	addTofSM (double time)
double	getTofSM (void)
void	fiTerm (double fi)
double	getFiTerm (void)
void	tof (double path)
	Update the tof estimation.
double	filter (double v_m, const CLHEP::HepVector &m_H, double v_d, double m_V)
double	cor_tanldep (double *p, double er)
	Correct the error according the current tanl value :
void	update_bit (int i)
int	insist (void) const
	Extractor :
int	type (void) const
int	trasan_id (void) const
double	r0 (void) const
double	mass (void) const
double	chiSq (void) const
double	chiSq_back (void) const
int	ndf_back (void) const
double	pathip (void) const
double	path_rd (void) const
double	path_ab (void) const
double *	pathl (void)
CLHEP::Hep3Vector *	mom (void)
double	tof (void) const
double	tof_kaon (void) const
double	tof_proton (void) const
double	p_kaon (void) const
double	p_proton (void) const
double	dchi2_max (void) const
double	r_max (void) const
unsigned int	nchits (void) const
unsigned int	nster (void) const
unsigned int	ncath (void) const
int	pat1 (void) const
int	pat2 (void) const
int	nhit_r (void) const
int	nhit_z (void) const
void	type (int t)
	Reinitialize (modificator).
void	trasan_id (int t)
void	insist (int t)
void	pathip (double pl)
void	p_kaon (double pl)
void	p_proton (double pl)
void	chiSq (double c)
void	chiSq_back (double c)
void	ndf_back (int n)
void	nchits (int n)
void	nster (int n)
void	add_nhit_r (void)
void	add_nhit_z (void)
double	PathL (int layer)
	Function to calculate the path length in the layer.
void	appendHitsMdc (KalFitHitMdc h)
	Functions for Mdc hits list.
void	HitsMdc (vector< KalFitHitMdc > &lh)
vector< KalFitHitMdc > &	HitsMdc (void)
KalFitHitMdc &	HitMdc (int i)
void	appendHelixSegs (KalFitHelixSeg s)
void	HelixSegs (vector< KalFitHelixSeg > &vs)
vector< KalFitHelixSeg > &	HelixSegs (void)
KalFitHelixSeg &	HelixSeg (int i)
void	order_wirhit (int index)
void	order_hits (void)
void	number_wirhit (void)
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot)
	Sets pivot position in a given mag field.
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot, double &pathl)
double	radius_numf (void) const
	Estimation of the radius in a given mag field.
double	getSigma (int layerId, double driftDist) const
double	getSigma (KalFitHitMdc &hitmdc, double tanlam, int lr, double dist) const
double	getDriftDist (KalFitHitMdc &hitmdc, double drifttime, int lr) const
double	getDriftTime (KalFitHitMdc &hitmdc, double toftime) const
double	getT0 (void) const
HepSymMatrix	getInitMatrix (void) const
double	getDigi () const
void	chgmass (int i)
int	nLayerUsed ()
void	resetLayerUsed ()
void	useLayer (int iLay)
	KalFitTrack (const HepPoint3D &pivot, const CLHEP::HepVector &a, const CLHEP::HepSymMatrix &Ea, unsigned int m, double chiSq, unsigned int nhits)
	constructor
	~KalFitTrack (void)
	destructor
double	intersect_cylinder (double r) const
	Intersection with different geometry.
double	intersect_yz_plane (const HepTransform3D &plane, double x) const
double	intersect_zx_plane (const HepTransform3D &plane, double y) const
double	intersect_xy_plane (double z) const
void	msgasmdc (double path, int index)
	Calculate multiple scattering angle.
void	ms (double path, const KalFitMaterial &m, int index)
void	eloss (double path, const KalFitMaterial &m, int index)
	Calculate total energy lost in material.
double	smoother_Mdc (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
	Kalman smoother for Mdc.
double	smoother_Mdc_csmalign (KalFitHelixSeg &seg, CLHEP::Hep3Vector &meas, int &flg, int csmflag)
double	smoother_Mdc (KalFitHitMdc &HitMdc, CLHEP::Hep3Vector &meas, KalFitHelixSeg &seg, double &dchi2, int csmflag)
double	update_hits (KalFitHitMdc &HitMdc, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, double &dtrack, double &dtracknew, double &dtdc, int csmflag)
	Include the Mdc wire hits.
double	update_hits (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	update_hits_csmalign (KalFitHelixSeg &HelixSeg, int inext, CLHEP::Hep3Vector &meas, int way, double &dchi2, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc, int csmflag)
double	chi2_next (Helix &H, KalFitHitMdc &HitMdc)
void	update_last (void)
	Record the current parameters as ..._last information :
void	update_forMdc (void)
void	point_last (const HepPoint3D &point)
	set and give out the last point of the track
const HepPoint3D &	point_last (void)
const HepPoint3D &	pivot_last (void) const
	returns helix parameters
const CLHEP::HepVector &	a_last (void) const
const CLHEP::HepSymMatrix &	Ea_last (void) const
const HepPoint3D &	pivot_forMdc (void) const
const CLHEP::HepVector &	a_forMdc (void) const
const CLHEP::HepSymMatrix &	Ea_forMdc (void) const
void	path_add (double path)
	Update the path length estimation.
void	addPathSM (double path)
double	getPathSM (void)
void	addTofSM (double time)
double	getTofSM (void)
void	fiTerm (double fi)
double	getFiTerm (void)
void	tof (double path)
	Update the tof estimation.
double	filter (double v_m, const CLHEP::HepVector &m_H, double v_d, double m_V)
double	cor_tanldep (double *p, double er)
	Correct the error according the current tanl value :
void	update_bit (int i)
int	insist (void) const
	Extractor :
int	type (void) const
int	trasan_id (void) const
double	r0 (void) const
double	mass (void) const
double	chiSq (void) const
double	chiSq_back (void) const
int	ndf_back (void) const
double	pathip (void) const
double	path_rd (void) const
double	path_ab (void) const
double *	pathl (void)
CLHEP::Hep3Vector *	mom (void)
double	tof (void) const
double	tof_kaon (void) const
double	tof_proton (void) const
double	p_kaon (void) const
double	p_proton (void) const
double	dchi2_max (void) const
double	r_max (void) const
unsigned int	nchits (void) const
unsigned int	nster (void) const
unsigned int	ncath (void) const
int	pat1 (void) const
int	pat2 (void) const
int	nhit_r (void) const
int	nhit_z (void) const
void	type (int t)
	Reinitialize (modificator).
void	trasan_id (int t)
void	insist (int t)
void	pathip (double pl)
void	p_kaon (double pl)
void	p_proton (double pl)
void	chiSq (double c)
void	chiSq_back (double c)
void	ndf_back (int n)
void	nchits (int n)
void	nster (int n)
void	add_nhit_r (void)
void	add_nhit_z (void)
double	PathL (int layer)
	Function to calculate the path length in the layer.
void	appendHitsMdc (KalFitHitMdc h)
	Functions for Mdc hits list.
void	HitsMdc (vector< KalFitHitMdc > &lh)
vector< KalFitHitMdc > &	HitsMdc (void)
KalFitHitMdc &	HitMdc (int i)
void	appendHelixSegs (KalFitHelixSeg s)
void	HelixSegs (vector< KalFitHelixSeg > &vs)
vector< KalFitHelixSeg > &	HelixSegs (void)
KalFitHelixSeg &	HelixSeg (int i)
void	order_wirhit (int index)
void	order_hits (void)
void	number_wirhit (void)
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot)
	Sets pivot position in a given mag field.
const HepPoint3D &	pivot_numf (const HepPoint3D &newPivot, double &pathl)
double	radius_numf (void) const
	Estimation of the radius in a given mag field.
double	getSigma (int layerId, double driftDist) const
double	getSigma (KalFitHitMdc &hitmdc, double tanlam, int lr, double dist) const
double	getDriftDist (KalFitHitMdc &hitmdc, double drifttime, int lr) const
double	getDriftTime (KalFitHitMdc &hitmdc, double toftime) const
double	getT0 (void) const
HepSymMatrix	getInitMatrix (void) const
double	getDigi () const
void	chgmass (int i)
int	nLayerUsed ()
void	resetLayerUsed ()
void	useLayer (int iLay)
Public Member Functions inherited from KalmanFit::Helix
	Helix (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	Constructor with pivot, helix parameter a, and its error matrix.
	Helix (const HepPoint3D &pivot, const HepVector &a)
	Constructor without error matrix.
	Helix (const HepPoint3D &position, const Hep3Vector &momentum, double charge)
	Constructor with position, momentum, and charge.
virtual	~Helix ()
	Destructor.
const HepPoint3D &	center (void) const
	returns position of helix center(z = 0.);
const HepPoint3D &	pivot (void) const
	returns pivot position.
double	radius (void) const
	returns radious of helix.
HepPoint3D	x (double dPhi=0.) const
	returns position after rotating angle dPhi in phi direction.
double *	x (double dPhi, double p[3]) const
HepPoint3D	x (double dPhi, HepSymMatrix &Ex) const
	returns position and convariance matrix(Ex) after rotation.
Hep3Vector	direction (double dPhi=0.) const
	returns direction vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi=0.) const
	returns momentum vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi, HepSymMatrix &Em) const
	returns momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepSymMatrix &Em) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepPoint3D &x, HepSymMatrix &Emx) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
double	dr (void) const
	returns an element of parameters.
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
double	curv (void) const
double	sinPhi0 (void) const
double	cosPhi0 (void) const
const HepVector &	a (void) const
	returns helix parameters.
const HepSymMatrix &	Ea (void) const
	returns error matrix.
double	approach (KalFitHitMdc &hit, bool doSagCorrection) const
double	approach (HepPoint3D pfwd, HepPoint3D pbwd, bool doSagCorrection) const
const HepVector &	a (const HepVector &newA)
	sets helix parameters.
const HepSymMatrix &	Ea (const HepSymMatrix &newdA)
	sets helix paramters and error matrix.
const HepPoint3D &	pivot (const HepPoint3D &newPivot)
	sets pivot position.
void	set (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	sets helix pivot position, parameters, and error matrix.
void	ignoreErrorMatrix (void)
double	bFieldZ (double)
	sets/returns z componet of the magnetic field.
double	bFieldZ (void) const
double	alpha (void) const
Helix &	operator= (const Helix &)
	Copy operator.
HepMatrix	delApDelA (const HepVector &ap) const
HepMatrix	delXDelA (double phi) const
HepMatrix	delMDelA (double phi) const
HepMatrix	del4MDelA (double phi, double mass) const
HepMatrix	del4MXDelA (double phi, double mass) const
	Helix (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	Constructor with pivot, helix parameter a, and its error matrix.
	Helix (const HepPoint3D &pivot, const HepVector &a)
	Constructor without error matrix.
	Helix (const HepPoint3D &position, const Hep3Vector &momentum, double charge)
	Constructor with position, momentum, and charge.
virtual	~Helix ()
	Destructor.
const HepPoint3D &	center (void) const
	returns position of helix center(z = 0.);
const HepPoint3D &	pivot (void) const
	returns pivot position.
double	radius (void) const
	returns radious of helix.
HepPoint3D	x (double dPhi=0.) const
	returns position after rotating angle dPhi in phi direction.
double *	x (double dPhi, double p[3]) const
HepPoint3D	x (double dPhi, HepSymMatrix &Ex) const
	returns position and convariance matrix(Ex) after rotation.
Hep3Vector	direction (double dPhi=0.) const
	returns direction vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi=0.) const
	returns momentum vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi, HepSymMatrix &Em) const
	returns momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepSymMatrix &Em) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepPoint3D &x, HepSymMatrix &Emx) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
double	dr (void) const
	returns an element of parameters.
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
double	curv (void) const
double	sinPhi0 (void) const
double	cosPhi0 (void) const
const HepVector &	a (void) const
	returns helix parameters.
const HepSymMatrix &	Ea (void) const
	returns error matrix.
double	approach (KalFitHitMdc &hit, bool doSagCorrection) const
double	approach (HepPoint3D pfwd, HepPoint3D pbwd, bool doSagCorrection) const
const HepVector &	a (const HepVector &newA)
	sets helix parameters.
const HepSymMatrix &	Ea (const HepSymMatrix &newdA)
	sets helix paramters and error matrix.
const HepPoint3D &	pivot (const HepPoint3D &newPivot)
	sets pivot position.
void	set (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	sets helix pivot position, parameters, and error matrix.
void	ignoreErrorMatrix (void)
double	bFieldZ (double)
	sets/returns z componet of the magnetic field.
double	bFieldZ (void) const
double	alpha (void) const
Helix &	operator= (const Helix &)
	Copy operator.
HepMatrix	delApDelA (const HepVector &ap) const
HepMatrix	delXDelA (double phi) const
HepMatrix	delMDelA (double phi) const
HepMatrix	del4MDelA (double phi, double mass) const
HepMatrix	del4MXDelA (double phi, double mass) const
	Helix (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	Constructor with pivot, helix parameter a, and its error matrix.
	Helix (const HepPoint3D &pivot, const HepVector &a)
	Constructor without error matrix.
	Helix (const HepPoint3D &position, const Hep3Vector &momentum, double charge)
	Constructor with position, momentum, and charge.
virtual	~Helix ()
	Destructor.
const HepPoint3D &	center (void) const
	returns position of helix center(z = 0.);
const HepPoint3D &	pivot (void) const
	returns pivot position.
double	radius (void) const
	returns radious of helix.
HepPoint3D	x (double dPhi=0.) const
	returns position after rotating angle dPhi in phi direction.
double *	x (double dPhi, double p[3]) const
HepPoint3D	x (double dPhi, HepSymMatrix &Ex) const
	returns position and convariance matrix(Ex) after rotation.
Hep3Vector	direction (double dPhi=0.) const
	returns direction vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi=0.) const
	returns momentum vector after rotating angle dPhi in phi direction.
Hep3Vector	momentum (double dPhi, HepSymMatrix &Em) const
	returns momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepSymMatrix &Em) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
HepLorentzVector	momentum (double dPhi, double mass, HepPoint3D &x, HepSymMatrix &Emx) const
	returns 4momentum vector after rotating angle dPhi in phi direction.
double	dr (void) const
	returns an element of parameters.
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
double	curv (void) const
double	sinPhi0 (void) const
double	cosPhi0 (void) const
const HepVector &	a (void) const
	returns helix parameters.
const HepSymMatrix &	Ea (void) const
	returns error matrix.
double	approach (KalFitHitMdc &hit, bool doSagCorrection) const
double	approach (HepPoint3D pfwd, HepPoint3D pbwd, bool doSagCorrection) const
const HepVector &	a (const HepVector &newA)
	sets helix parameters.
const HepSymMatrix &	Ea (const HepSymMatrix &newdA)
	sets helix paramters and error matrix.
const HepPoint3D &	pivot (const HepPoint3D &newPivot)
	sets pivot position.
void	set (const HepPoint3D &pivot, const HepVector &a, const HepSymMatrix &Ea)
	sets helix pivot position, parameters, and error matrix.
void	ignoreErrorMatrix (void)
double	bFieldZ (double)
	sets/returns z componet of the magnetic field.
double	bFieldZ (void) const
double	alpha (void) const
Helix &	operator= (const Helix &)
	Copy operator.
HepMatrix	delApDelA (const HepVector &ap) const
HepMatrix	delXDelA (double phi) const
HepMatrix	delMDelA (double phi) const
HepMatrix	del4MDelA (double phi, double mass) const
HepMatrix	del4MXDelA (double phi, double mass) const

Static Public Member Functions
static void	setT0 (double t0)
static void	setInitMatrix (HepSymMatrix m)
static void	setMdcCalibFunSvc (const IMdcCalibFunSvc *calibsvc)
static void	setMagneticFieldSvc (IBesMagFieldSvc *)
static void	setIMdcGeomSvc (IMdcGeomSvc *igeomsvc)
static void	setMdcDigiCol (MdcDigiCol *digicol)
static int	nmass (void)
static double	mass (int i)
static void	LR (int x)
static int	lead (void)
	Magnetic field map.
static void	lead (int i)
static int	back (void)
static void	back (int i)
static int	resol (void)
static void	resol (int i)
static int	numf (void)
static void	numf (int i)
static void	setT0 (double t0)
static void	setInitMatrix (HepSymMatrix m)
static void	setMdcCalibFunSvc (const IMdcCalibFunSvc *calibsvc)
static void	setMagneticFieldSvc (IBesMagFieldSvc *)
static void	setIMdcGeomSvc (IMdcGeomSvc *igeomsvc)
static void	setMdcDigiCol (MdcDigiCol *digicol)
static int	nmass (void)
static double	mass (int i)
static void	LR (int x)
static int	lead (void)
	Magnetic field map.
static void	lead (int i)
static int	back (void)
static void	back (int i)
static int	resol (void)
static void	resol (int i)
static int	numf (void)
static void	numf (int i)
static void	setT0 (double t0)
static void	setInitMatrix (HepSymMatrix m)
static void	setMdcCalibFunSvc (const IMdcCalibFunSvc *calibsvc)
static void	setMagneticFieldSvc (IBesMagFieldSvc *)
static void	setIMdcGeomSvc (IMdcGeomSvc *igeomsvc)
static void	setMdcDigiCol (MdcDigiCol *digicol)
static int	nmass (void)
static double	mass (int i)
static void	LR (int x)
static int	lead (void)
	Magnetic field map.
static void	lead (int i)
static int	back (void)
static void	back (int i)
static int	resol (void)
static void	resol (int i)
static int	numf (void)
static void	numf (int i)

Static Public Attributes
static double	mdcGasRadlen_ = 0.
static int	tprop_ = 1
	for signal propagation correction
static int	debug_ = 0
	for debug
static double	chi2_hitf_ = 1000
	Cut chi2 for each hit.
static double	chi2_hits_ = 1000
static int	numf_ = 0
	Flag for treatment of non-uniform mag field.
static int	inner_steps_ = 3
static int	outer_steps_ = 3
static double	dchi2cutf_anal [43] = { 0. }
static double	dchi2cuts_anal [43] = { 0. }
static double	dchi2cutf_calib [43] = { 0. }
static double	dchi2cuts_calib [43] = { 0. }
static int	numfcor_ = 1
	NUMF treatment improved.
static double	Bznom_ = 10
static int	steplev_ = 0
static int	Tof_correc_ = 1
	Flag for TOF correction.
static int	strag_ = 1
	Flag to take account of energy loss straggling :
static double	factor_strag_ = 0.4
	factor of energy loss straggling for electron
static int	nmdc_hit2_ = 500
	Cut chi2 for each hit.
static double	chi2mdc_hit2_
static int	tofall_ = 1
static int	resolflag_ = 0
	wire resoltion flag
static int	LR_ = 1
	Use L/R decision from MdcRecHit information :
static int	drifttime_choice_ = 0
	the drifttime choice
Static Public Attributes inherited from KalmanFit::Helix
static const double	ConstantAlpha = 333.564095
	Constant alpha for uniform field.

Detailed Description

Description of a track class (<- Helix.cc).

Definition at line 32 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Constructor & Destructor Documentation

KalFitTrack() [1/3]

KalFitTrack::KalFitTrack	(	const HepPoint3D &	pivot,
		const CLHEP::HepVector &	a,
		const CLHEP::HepSymMatrix &	Ea,
		unsigned int	m,
		double	chiSq,
		unsigned int	nhits )

constructor

~KalFitTrack() [1/3]

KalFitTrack::~KalFitTrack

(

void

)

destructor

Definition at line 111 of file KalFitTrack.cxx.

                                {
  // delete all objects
}

KalFitTrack() [2/3]

KalFitTrack::KalFitTrack	(	const HepPoint3D &	pivot,
		const CLHEP::HepVector &	a,
		const CLHEP::HepSymMatrix &	Ea,
		unsigned int	m,
		double	chiSq,
		unsigned int	nhits )

constructor

~KalFitTrack() [2/3]

KalFitTrack::~KalFitTrack

(

void

)

destructor

KalFitTrack() [3/3]

KalFitTrack::KalFitTrack	(	const HepPoint3D &	pivot,
		const CLHEP::HepVector &	a,
		const CLHEP::HepSymMatrix &	Ea,
		unsigned int	m,
		double	chiSq,
		unsigned int	nhits )

constructor

~KalFitTrack() [3/3]

KalFitTrack::~KalFitTrack

(

void

)

destructor

Member Function Documentation

a_forMdc() [1/3]

const CLHEP::HepVector & KalFitTrack::a_forMdc ( void ) const

inline

Definition at line 142 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return a_forMdc_; }

a_forMdc() [2/3]

const CLHEP::HepVector & KalFitTrack::a_forMdc ( void ) const

inline

Definition at line 142 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return a_forMdc_; }

a_forMdc() [3/3]

const CLHEP::HepVector & KalFitTrack::a_forMdc ( void ) const

inline

Definition at line 142 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return a_forMdc_; }

a_last() [1/3]

const CLHEP::HepVector & KalFitTrack::a_last ( void ) const

inline

Definition at line 138 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return a_last_; }

a_last() [2/3]

const CLHEP::HepVector & KalFitTrack::a_last ( void ) const

inline

Definition at line 138 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return a_last_; }

a_last() [3/3]

const CLHEP::HepVector & KalFitTrack::a_last ( void ) const

inline

Definition at line 138 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return a_last_; }

add_nhit_r() [1/3]

void KalFitTrack::add_nhit_r ( void )

inline

Definition at line 210 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ nhit_r_++; }

add_nhit_r() [2/3]

void KalFitTrack::add_nhit_r ( void )

inline

Definition at line 210 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ nhit_r_++; }

add_nhit_r() [3/3]

void KalFitTrack::add_nhit_r ( void )

inline

Definition at line 210 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ nhit_r_++; }

add_nhit_z() [1/3]

void KalFitTrack::add_nhit_z ( void )

inline

Definition at line 211 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ nhit_z_++; }

add_nhit_z() [2/3]

void KalFitTrack::add_nhit_z ( void )

inline

Definition at line 211 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ nhit_z_++; }

add_nhit_z() [3/3]

void KalFitTrack::add_nhit_z ( void )

inline

Definition at line 211 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ nhit_z_++; }

addPathSM() [1/3]

void KalFitTrack::addPathSM

(

double

path

)

Definition at line 1359 of file KalFitTrack.cxx.

{ pathSM_ += path; }

Referenced by KalFitAlg::smoother_anal().

addPathSM() [2/3]

void KalFitTrack::addPathSM

(

double

path

)

addPathSM() [3/3]

void KalFitTrack::addPathSM

(

double

path

)

addTofSM() [1/3]

void KalFitTrack::addTofSM

(

double

time

)

Definition at line 1361 of file KalFitTrack.cxx.

{ tof2_ += time; }

Referenced by KalFitAlg::smoother_anal().

addTofSM() [2/3]

void KalFitTrack::addTofSM

(

double

time

)

addTofSM() [3/3]

void KalFitTrack::addTofSM

(

double

time

)

appendHelixSegs() [1/3]

void KalFitTrack::appendHelixSegs

(

KalFitHelixSeg

s

)

appendHelixSegs() [2/3]

void KalFitTrack::appendHelixSegs

(

KalFitHelixSeg

s

)

appendHelixSegs() [3/3]

void KalFitTrack::appendHelixSegs

(

KalFitHelixSeg

s

)

appendHitsMdc() [1/3]

void KalFitTrack::appendHitsMdc

(

KalFitHitMdc

h

)

Functions for Mdc hits list.

Definition at line 1874 of file KalFitTrack.cxx.

{ HitsMdc_.push_back( h ); }

Referenced by KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), and KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew().

appendHitsMdc() [2/3]

void KalFitTrack::appendHitsMdc

(

KalFitHitMdc

h

)

Functions for Mdc hits list.

appendHitsMdc() [3/3]

void KalFitTrack::appendHitsMdc

(

KalFitHitMdc

h

)

Functions for Mdc hits list.

back() [1/6]

void KalFitTrack::back ( int i )

static

Definition at line 1866 of file KalFitTrack.cxx.

{ back_ = i; }

back() [2/6]

void KalFitTrack::back ( int i )

static

back() [3/6]

void KalFitTrack::back ( int i )

static

back() [4/6]

int KalFitTrack::back ( void )

static

Definition at line 1867 of file KalFitTrack.cxx.

{ return back_; }

Referenced by KalFitAlg::initialize().

back() [5/6]

int KalFitTrack::back ( void )

static

back() [6/6]

int KalFitTrack::back ( void )

static

chgmass() [1/3]

void KalFitTrack::chgmass

(

int

i

)

Definition at line 1859 of file KalFitTrack.cxx.

                                   {
  mass_   = MASS[i];
  l_mass_ = i;
}

Referenced by KalFitAlg::complete_track(), and KalFitAlg::complete_track().

chgmass() [2/3]

void KalFitTrack::chgmass

(

int

i

)

chgmass() [3/3]

void KalFitTrack::chgmass

(

int

i

)

chi2_next() [1/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc )

Definition at line 2915 of file KalFitTrack.cxx.

                                                              {
 
  double            lr        = HitMdc.LR();
  const KalFitWire& Wire      = HitMdc.wire();
  int               wire_ID   = Wire.geoID();
  int               layerid   = HitMdc.wire().layer().layerId();
  double            entrangle = HitMdc.rechitptr()->getEntra();
 
  HepPoint3D fwd( Wire.fwd() );
  HepPoint3D bck( Wire.bck() );
  Hep3Vector wire = (Hep3Vector)fwd - (Hep3Vector)bck;
  Helix      work = H;
  work.ignoreErrorMatrix();
  work.pivot( ( fwd + bck ) * .5 );
  HepPoint3D x0kal = ( work.x( 0 ).z() - bck.z() ) / wire.z() * wire + bck;
  H.pivot( x0kal );
 
  Hep3Vector meas = H.momentum( 0 ).cross( wire ).unit();
 
  if ( wire_ID < 0 || wire_ID > 6796 )
  { // bes
    std::cout << "KalFitTrack : wire_ID problem : " << wire_ID << std::endl;
    return DBL_MAX;
  }
 
  double x[3]    = { pivot().x(), pivot().y(), pivot().z() };
  double pmom[3] = { momentum().x(), momentum().y(), momentum().z() };
  double tofest( 0 );
  double phi  = fmod( phi0() + M_PI4, M_PI2 );
  double csf0 = cos( phi );
  double snf0 = ( 1. - csf0 ) * ( 1. + csf0 );
  snf0        = sqrt( ( snf0 > 0. ) ? snf0 : 0. );
  if ( phi > M_PI ) snf0 = -snf0;
 
  if ( Tof_correc_ )
  {
    Hep3Vector ip( 0, 0, 0 );
    Helix      work = *(Helix*)this;
    work.ignoreErrorMatrix();
    work.pivot( ip );
    double phi_ip = work.phi0();
    if ( fabs( phi - phi_ip ) > M_PI )
    {
      if ( phi > phi_ip ) phi -= 2 * M_PI;
      else phi_ip -= 2 * M_PI;
    }
    double t = tanl();
    double l = fabs( radius() * ( phi - phi_ip ) * sqrt( 1 + t * t ) );
    double pmag( sqrt( 1.0 + t * t ) / kappa() );
    double mass_over_p( mass_ / pmag );
    double beta( 1.0 / sqrt( 1.0 + mass_over_p * mass_over_p ) );
    tofest = l / ( 29.9792458 * beta );
    //  if(csmflag==1 && HitMdc.wire().y()>0.)  tofest= -1. * tofest;
  }
 
  const HepSymMatrix& ea  = H.Ea();
  const HepVector&    v_a = H.a();
  double              dchi2R( DBL_MAX ), dchi2L( DBL_MAX );
 
  HepVector v_H( 5, 0 );
  v_H[0]         = -csf0 * meas.x() - snf0 * meas.y();
  v_H[3]         = -meas.z();
  HepMatrix v_HT = v_H.T();
 
  double    estim        = ( v_HT * v_a )[0];
  HepVector ea_v_H       = ea * v_H;
  HepMatrix ea_v_HT      = ( ea_v_H ).T();
  HepVector v_H_T_ea_v_H = v_HT * ea_v_H;
 
  HepSymMatrix eaNewL( 5 ), eaNewR( 5 );
  HepVector    aNewL( 5 ), aNewR( 5 );
 
  // double time = HitMdc.tdc();
  // if (Tof_correc_)
  //  time = time - tofest;
  double drifttime = getDriftTime( HitMdc, tofest );
  double ddl       = getDriftDist( HitMdc, drifttime, 0 );
  double ddr       = getDriftDist( HitMdc, drifttime, 1 );
  double erddl     = getSigma( HitMdc, H.a()[4], 0, ddl );
  double erddr     = getSigma( HitMdc, H.a()[4], 1, ddr );
 
  double dmeas2[2]    = { 0.1 * ddl, 0.1 * ddr };
  double er_dmeas2[2] = { 0., 0. };
  if ( resolflag_ == 1 )
  {
    er_dmeas2[0] = 0.1 * erddl;
    er_dmeas2[1] = 0.1 * erddr;
  }
  else if ( resolflag_ == 0 )
  {
    //  int layid = HitMdc.wire().layer().layerId();
    //  double sigma = getSigma(layid, dd);
    //  er_dmeas2[0] = er_dmeas2[1] = sigma;
  }
 
  if ( ( LR_ == 0 && lr != 1.0 ) || ( LR_ == 1 && lr == -1.0 ) )
  {
 
    double er_dmeasL, dmeasL;
    if ( Tof_correc_ )
    {
      dmeasL    = ( -1.0 ) * fabs( dmeas2[0] );
      er_dmeasL = er_dmeas2[0];
    }
    else
    {
      dmeasL    = ( -1.0 ) * fabs( HitMdc.dist()[0] );
      er_dmeasL = HitMdc.erdist()[0];
    }
 
    double AL = 1 / ( ( v_H_T_ea_v_H )[0] + er_dmeasL * er_dmeasL );
    eaNewL.assign( ea - ea_v_H * AL * ea_v_HT );
    double RkL = 1 - ( v_H_T_ea_v_H )[0] * AL;
    if ( 0. == RkL ) RkL = 1.e-4;
 
    HepVector diffL = ea_v_H * AL * ( dmeasL - estim );
    aNewL           = v_a + diffL;
    double sigL     = dmeasL - ( v_HT * aNewL )[0];
    dchi2L          = ( sigL * sigL ) / ( RkL * er_dmeasL * er_dmeasL );
  }
 
  if ( ( LR_ == 0 && lr != -1.0 ) || ( LR_ == 1 && lr == 1.0 ) )
  {
 
    double er_dmeasR, dmeasR;
    if ( Tof_correc_ )
    {
      dmeasR    = dmeas2[1];
      er_dmeasR = er_dmeas2[1];
    }
    else
    {
      dmeasR    = fabs( HitMdc.dist()[1] );
      er_dmeasR = HitMdc.erdist()[1];
    }
 
    double AR = 1 / ( ( v_H_T_ea_v_H )[0] + er_dmeasR * er_dmeasR );
    eaNewR.assign( ea - ea_v_H * AR * ea_v_HT );
    double RkR = 1 - ( v_H_T_ea_v_H )[0] * AR;
    if ( 0. == RkR ) RkR = 1.e-4;
 
    HepVector diffR = ea_v_H * AR * ( dmeasR - estim );
    aNewR           = v_a + diffR;
    double sigR     = dmeasR - ( v_HT * aNewR )[0];
    dchi2R          = ( sigR * sigR ) / ( RkR * er_dmeasR * er_dmeasR );
  }
 
  if ( dchi2R < dchi2L )
  {
    H.a( aNewR );
    H.Ea( eaNewR );
  }
  else
  {
    H.a( aNewL );
    H.Ea( eaNewL );
  }
  return ( ( dchi2R < dchi2L ) ? dchi2R : dchi2L );
}

chi2_next() [2/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc )

chi2_next() [3/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc )

chi2_next() [4/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc,
		int	csmflag )

Definition at line 3075 of file KalFitTrack.cxx.

                                                                           {
 
  double            lr        = HitMdc.LR();
  const KalFitWire& Wire      = HitMdc.wire();
  int               wire_ID   = Wire.geoID();
  int               layerid   = HitMdc.wire().layer().layerId();
  double            entrangle = HitMdc.rechitptr()->getEntra();
 
  HepPoint3D fwd( Wire.fwd() );
  HepPoint3D bck( Wire.bck() );
  Hep3Vector wire = (Hep3Vector)fwd - (Hep3Vector)bck;
  Helix      work = H;
  work.ignoreErrorMatrix();
  work.pivot( ( fwd + bck ) * .5 );
  HepPoint3D x0kal = ( work.x( 0 ).z() - bck.z() ) / wire.z() * wire + bck;
  H.pivot( x0kal );
 
  Hep3Vector meas = H.momentum( 0 ).cross( wire ).unit();
 
  if ( wire_ID < 0 || wire_ID > 6796 )
  { // bes
    std::cout << "KalFitTrack : wire_ID problem : " << wire_ID << std::endl;
    return DBL_MAX;
  }
 
  double x[3]    = { pivot().x(), pivot().y(), pivot().z() };
  double pmom[3] = { momentum().x(), momentum().y(), momentum().z() };
  double tofest( 0 );
  double phi  = fmod( phi0() + M_PI4, M_PI2 );
  double csf0 = cos( phi );
  double snf0 = ( 1. - csf0 ) * ( 1. + csf0 );
  snf0        = sqrt( ( snf0 > 0. ) ? snf0 : 0. );
  if ( phi > M_PI ) snf0 = -snf0;
 
  if ( Tof_correc_ )
  {
    Hep3Vector ip( 0, 0, 0 );
    Helix      work = *(Helix*)this;
    work.ignoreErrorMatrix();
    work.pivot( ip );
    double phi_ip = work.phi0();
    if ( fabs( phi - phi_ip ) > M_PI )
    {
      if ( phi > phi_ip ) phi -= 2 * M_PI;
      else phi_ip -= 2 * M_PI;
    }
    double t = tanl();
    double l = fabs( radius() * ( phi - phi_ip ) * sqrt( 1 + t * t ) );
    double pmag( sqrt( 1.0 + t * t ) / kappa() );
    double mass_over_p( mass_ / pmag );
    double beta( 1.0 / sqrt( 1.0 + mass_over_p * mass_over_p ) );
    tofest = l / ( 29.9792458 * beta );
    if ( csmflag == 1 && HitMdc.wire().y() > 0. ) tofest = -1. * tofest;
  }
 
  const HepSymMatrix& ea  = H.Ea();
  const HepVector&    v_a = H.a();
  double              dchi2R( DBL_MAX ), dchi2L( DBL_MAX );
 
  HepVector v_H( 5, 0 );
  v_H[0]         = -csf0 * meas.x() - snf0 * meas.y();
  v_H[3]         = -meas.z();
  HepMatrix v_HT = v_H.T();
 
  double    estim        = ( v_HT * v_a )[0];
  HepVector ea_v_H       = ea * v_H;
  HepMatrix ea_v_HT      = ( ea_v_H ).T();
  HepVector v_H_T_ea_v_H = v_HT * ea_v_H;
 
  HepSymMatrix eaNewL( 5 ), eaNewR( 5 );
  HepVector    aNewL( 5 ), aNewR( 5 );
 
  // double time = HitMdc.tdc();
  // if (Tof_correc_)
  //  time = time - tofest;
  double drifttime = getDriftTime( HitMdc, tofest );
  double ddl       = getDriftDist( HitMdc, drifttime, 0 );
  double ddr       = getDriftDist( HitMdc, drifttime, 1 );
  double erddl     = getSigma( HitMdc, H.a()[4], 0, ddl );
  double erddr     = getSigma( HitMdc, H.a()[4], 1, ddr );
 
  double dmeas2[2]    = { 0.1 * ddl, 0.1 * ddr };
  double er_dmeas2[2] = { 0., 0. };
  if ( resolflag_ == 1 )
  {
    er_dmeas2[0] = 0.1 * erddl;
    er_dmeas2[1] = 0.1 * erddr;
  }
  else if ( resolflag_ == 0 )
  {
    //  int layid = HitMdc.wire().layer().layerId();
    //  double sigma = getSigma(layid, dd);
    //  er_dmeas2[0] = er_dmeas2[1] = sigma;
  }
 
  if ( ( LR_ == 0 && lr != 1.0 ) || ( LR_ == 1 && lr == -1.0 ) )
  {
 
    double er_dmeasL, dmeasL;
    if ( Tof_correc_ )
    {
      dmeasL    = ( -1.0 ) * fabs( dmeas2[0] );
      er_dmeasL = er_dmeas2[0];
    }
    else
    {
      dmeasL    = ( -1.0 ) * fabs( HitMdc.dist()[0] );
      er_dmeasL = HitMdc.erdist()[0];
    }
 
    double AL = 1 / ( ( v_H_T_ea_v_H )[0] + er_dmeasL * er_dmeasL );
    eaNewL.assign( ea - ea_v_H * AL * ea_v_HT );
    double RkL = 1 - ( v_H_T_ea_v_H )[0] * AL;
    if ( 0. == RkL ) RkL = 1.e-4;
 
    HepVector diffL = ea_v_H * AL * ( dmeasL - estim );
    aNewL           = v_a + diffL;
    double sigL     = dmeasL - ( v_HT * aNewL )[0];
    dchi2L          = ( sigL * sigL ) / ( RkL * er_dmeasL * er_dmeasL );
  }
 
  if ( ( LR_ == 0 && lr != -1.0 ) || ( LR_ == 1 && lr == 1.0 ) )
  {
 
    double er_dmeasR, dmeasR;
    if ( Tof_correc_ )
    {
      dmeasR    = dmeas2[1];
      er_dmeasR = er_dmeas2[1];
    }
    else
    {
      dmeasR    = fabs( HitMdc.dist()[1] );
      er_dmeasR = HitMdc.erdist()[1];
    }
 
    double AR = 1 / ( ( v_H_T_ea_v_H )[0] + er_dmeasR * er_dmeasR );
    eaNewR.assign( ea - ea_v_H * AR * ea_v_HT );
    double RkR = 1 - ( v_H_T_ea_v_H )[0] * AR;
    if ( 0. == RkR ) RkR = 1.e-4;
 
    HepVector diffR = ea_v_H * AR * ( dmeasR - estim );
    aNewR           = v_a + diffR;
    double sigR     = dmeasR - ( v_HT * aNewR )[0];
    dchi2R          = ( sigR * sigR ) / ( RkR * er_dmeasR * er_dmeasR );
  }
 
  if ( dchi2R < dchi2L )
  {
    H.a( aNewR );
    H.Ea( eaNewR );
  }
  else
  {
    H.a( aNewL );
    H.Ea( eaNewL );
  }
  return ( ( dchi2R < dchi2L ) ? dchi2R : dchi2L );
}

chi2_next() [5/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc,
		int	csmflag )

chi2_next() [6/6]

double KalFitTrack::chi2_next	(	Helix &	H,
		KalFitHitMdc &	HitMdc,
		int	csmflag )

chiSq() [1/6]

void KalFitTrack::chiSq ( double c )

inline

Definition at line 204 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ chiSq_ = c; }

chiSq() [2/6]

void KalFitTrack::chiSq ( double c )

inline

Definition at line 204 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ chiSq_ = c; }

chiSq() [3/6]

void KalFitTrack::chiSq ( double c )

inline

Definition at line 204 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ chiSq_ = c; }

chiSq() [4/6]

double KalFitTrack::chiSq ( void ) const

inline

Definition at line 173 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_; }

Referenced by KalFitAlg::fillTds(), KalFitAlg::fillTds_ip(), KalFitAlg::fillTds_lead(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitTrack(), and KalFitAlg::start_seed().

chiSq() [5/6]

double KalFitTrack::chiSq ( void ) const

inline

Definition at line 173 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_; }

chiSq() [6/6]

double KalFitTrack::chiSq ( void ) const

inline

Definition at line 173 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_; }

chiSq_back() [1/6]

void KalFitTrack::chiSq_back ( double c )

inline

Definition at line 205 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ chiSq_back_ = c; }

chiSq_back() [2/6]

void KalFitTrack::chiSq_back ( double c )

inline

Definition at line 205 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ chiSq_back_ = c; }

chiSq_back() [3/6]

void KalFitTrack::chiSq_back ( double c )

inline

Definition at line 205 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ chiSq_back_ = c; }

chiSq_back() [4/6]

double KalFitTrack::chiSq_back ( void ) const

inline

Definition at line 174 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_back_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::start_seed().

chiSq_back() [5/6]

double KalFitTrack::chiSq_back ( void ) const

inline

Definition at line 174 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_back_; }

chiSq_back() [6/6]

double KalFitTrack::chiSq_back ( void ) const

inline

Definition at line 174 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return chiSq_back_; }

cor_tanldep() [1/3]

double KalFitTrack::cor_tanldep	(	double *	p,
		double	er )

Correct the error according the current tanl value :

cor_tanldep() [2/3]

double KalFitTrack::cor_tanldep	(	double *	p,
		double	er )

Correct the error according the current tanl value :

cor_tanldep() [3/3]

double KalFitTrack::cor_tanldep	(	double *	p,
		double	er )

Correct the error according the current tanl value :

dchi2_max() [1/3]

double KalFitTrack::dchi2_max ( void ) const

inline

Definition at line 186 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return dchi2_max_; }

dchi2_max() [2/3]

double KalFitTrack::dchi2_max ( void ) const

inline

Definition at line 186 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return dchi2_max_; }

dchi2_max() [3/3]

double KalFitTrack::dchi2_max ( void ) const

inline

Definition at line 186 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return dchi2_max_; }

Ea_forMdc() [1/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_forMdc ( void ) const

inline

Definition at line 143 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return Ea_forMdc_; }

Ea_forMdc() [2/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_forMdc ( void ) const

inline

Definition at line 143 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return Ea_forMdc_; }

Ea_forMdc() [3/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_forMdc ( void ) const

inline

Definition at line 143 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return Ea_forMdc_; }

Ea_last() [1/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_last ( void ) const

inline

Definition at line 139 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return Ea_last_; }

Ea_last() [2/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_last ( void ) const

inline

Definition at line 139 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return Ea_last_; }

Ea_last() [3/3]

const CLHEP::HepSymMatrix & KalFitTrack::Ea_last ( void ) const

inline

Definition at line 139 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return Ea_last_; }

eloss() [1/3]

void KalFitTrack::eloss	(	double	path,
		const KalFitMaterial &	m,
		int	index )

Calculate total energy lost in material.

Definition at line 272 of file KalFitTrack.cxx.

                                                                                {
#ifdef YDEBUG
  cout << "eloss():ea before: " << Ea() << endl;
#endif
  HepVector v_a     = a();
  double    v_a_2_2 = v_a[2] * v_a[2];
  double    v_a_4_2 = v_a[4] * v_a[4];
  double    pmag    = 1 / fabs( v_a[2] ) * sqrt( 1 + v_a_4_2 );
  double    psq     = pmag * pmag;
  double    E       = sqrt( mass_ * mass_ + psq );
  double    dE      = material.dE( mass_, path, pmag );
  // std::cout<<" eloss(): dE: "<<dE<<std::endl;//wangll
 
  if ( index > 0 ) psq += dE * ( dE + 2 * sqrt( mass_ * mass_ + psq ) );
  else
  {
    double dE_max = E - mass_;
    if ( dE < dE_max ) psq += dE * ( dE - 2 * sqrt( mass_ * mass_ + psq ) );
    else psq = -1.0;
  }
 
  if ( tofall_ && index < 0 )
  {
    // Kaon case :
    if ( p_kaon_ > 0 )
    {
      double psq_kaon = p_kaon_ * p_kaon_;
      double dE_kaon  = material.dE( MASS[3], path, p_kaon_ );
      psq_kaon += dE_kaon * ( dE_kaon - 2 * sqrt( MASS[3] * MASS[3] + psq_kaon ) );
      if ( psq_kaon < 0 ) psq_kaon = 0;
      p_kaon_ = sqrt( psq_kaon );
    }
    // Proton case :
    if ( p_proton_ > 0 )
    {
      double psq_proton = p_proton_ * p_proton_;
      double dE_proton  = material.dE( MASS[4], path, p_proton_ );
      psq_proton += dE_proton * ( dE_proton - 2 * sqrt( MASS[4] * MASS[4] + psq_proton ) );
      if ( psq_proton < 0 ) psq_proton = 0;
      p_proton_ = sqrt( psq_proton );
    }
  }
 
  double dpt;
  // cout<<"eloss(): psq = "<<psq<<endl;//wangll
  if ( psq < 0 ) dpt = 9999;
  else dpt = v_a[2] * pmag / sqrt( psq );
    // cout<<"eloss():k:   "<<v_a[2]<<"  k'  "<<dpt<<endl;//wangll
#ifdef YDEBUG
  cout << "eloss():k:   " << v_a[2] << "  k'  " << dpt << endl;
#endif
  // attempt to take account of energy loss for error matrix
 
  HepSymMatrix ea       = Ea();
  double       r_E_prim = ( E + dE ) / E;
 
  // 1/ Straggling in the energy loss :
  if ( strag_ )
  {
    double del_E( 0 );
    if ( l_mass_ == 0 ) { del_E = dE * factor_strag_; }
    else { del_E = material.del_E( mass_, path, pmag ); }
 
    ea[2][2] += ( v_a[2] * v_a[2] ) * E * E * del_E * del_E / ( psq * psq );
  }
 
  // Effect of the change of curvature (just variables change):
  double dpt2 = dpt * dpt;
  double A    = dpt * dpt2 / ( v_a_2_2 * v_a[2] ) * r_E_prim;
  double B    = v_a[4] / ( 1 + v_a_4_2 ) * dpt * ( 1 - dpt2 / v_a_2_2 * r_E_prim );
 
  double ea_2_0 = A * ea[2][0] + B * ea[4][0];
  double ea_2_1 = A * ea[2][1] + B * ea[4][1];
  double ea_2_2 = A * A * ea[2][2] + 2 * A * B * ea[2][4] + B * B * ea[4][4];
  double ea_2_3 = A * ea[2][3] + B * ea[4][3];
  double ea_2_4 = A * ea[2][4] + B * ea[4][4];
 
  v_a[2] = dpt;
  a( v_a );
 
  ea[2][0] = ea_2_0;
  // std::cout<<"ea[2][0] is "<<ea[2][0]<<" ea(3,1) is "<<ea(3,1)<<std::endl;
  ea[2][1] = ea_2_1;
  // std::cout<<"ea[2][2] is "<<ea[2][2]<<" ea(3,3) is "<<ea(3,3)<<std::endl;
  ea[2][2] = ea_2_2;
  ea[2][3] = ea_2_3;
  ea[2][4] = ea_2_4;
 
  Ea( ea );
  // cout<<"eloss():dE:   "<<dE<<endl;
  // cout<<"eloss():A:   "<<A<<"   B:    "<<B<<endl;
  // cout<<"eloss():ea after: "<<Ea()<<endl;
  r0_ = fabs( center().perp() - fabs( radius() ) );
}

Referenced by KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), KalFitAlg::smoother_calib(), KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

eloss() [2/3]

void KalFitTrack::eloss	(	double	path,
		const KalFitMaterial &	m,
		int	index )

Calculate total energy lost in material.

eloss() [3/3]

void KalFitTrack::eloss	(	double	path,
		const KalFitMaterial &	m,
		int	index )

Calculate total energy lost in material.

filter() [1/3]

double KalFitTrack::filter	(	double	v_m,
		const CLHEP::HepVector &	m_H,
		double	v_d,
		double	m_V )

filter() [2/3]

double KalFitTrack::filter	(	double	v_m,
		const CLHEP::HepVector &	m_H,
		double	v_d,
		double	m_V )

filter() [3/3]

double KalFitTrack::filter	(	double	v_m,
		const CLHEP::HepVector &	m_H,
		double	v_d,
		double	m_V )

fiTerm() [1/3]

void KalFitTrack::fiTerm

(

double

fi

)

Definition at line 1363 of file KalFitTrack.cxx.

{ fiTerm_ = fi; }

Referenced by KalFitAlg::smoother_anal().

fiTerm() [2/3]

void KalFitTrack::fiTerm

(

double

fi

)

fiTerm() [3/3]

void KalFitTrack::fiTerm

(

double

fi

)

getDigi() [1/3]

double KalFitTrack::getDigi

(

)

const

getDigi() [2/3]

double KalFitTrack::getDigi

(

)

const

getDigi() [3/3]

double KalFitTrack::getDigi

(

)

const

getDriftDist() [1/3]

double KalFitTrack::getDriftDist	(	KalFitHitMdc &	hitmdc,
		double	drifttime,
		int	lr ) const

Definition at line 186 of file KalFitTrack2.cxx.

                                                                                       {
  int layerid = hitmdc.wire().layer().layerId();
  int cellid  = MdcID::wire( hitmdc.rechitptr()->getMdcId() );
  if ( debug_ == 4 ) { std::cout << "the cellid is .." << cellid << std::endl; }
  double entrangle = hitmdc.rechitptr()->getEntra();
 
  // std::cout<<" entrangle: "<<entrangle<<std::endl;
 
  return CalibFunSvc_->driftTimeToDist( drifttime, layerid, cellid, lr, entrangle );
}

Referenced by chi2_next(), and chi2_next().

getDriftDist() [2/3]

double KalFitTrack::getDriftDist	(	KalFitHitMdc &	hitmdc,
		double	drifttime,
		int	lr ) const

getDriftDist() [3/3]

double KalFitTrack::getDriftDist	(	KalFitHitMdc &	hitmdc,
		double	drifttime,
		int	lr ) const

getDriftTime() [1/3]

double KalFitTrack::getDriftTime	(	KalFitHitMdc &	hitmdc,
		double	toftime ) const

Definition at line 62 of file KalFitTrack2.cxx.

                                                                             {
  const double vinner = 220.0e8; // cm/s
  const double vouter = 240.0e8; // cm/s
 
  int               layerid = hitmdc.wire().layer().layerId();
  double            zhit    = ( hitmdc.rechitptr() )->getZhit();
  const KalFitWire* wire    = &( hitmdc.wire() );
  HepPoint3D        fPoint  = wire->fwd();
  HepPoint3D        bPoint  = wire->bck();
 
  // unit is centimeter
  double wireLen = ( fPoint - bPoint ).x() * ( fPoint - bPoint ).x() +
                   ( fPoint - bPoint ).y() * ( fPoint - bPoint ).y() +
                   ( fPoint - bPoint ).z() * ( fPoint - bPoint ).z();
  wireLen         = sqrt( wireLen );
  double wireZLen = fabs( fPoint.z() - bPoint.z() );
  double tp       = 0.;
  double vp       = 0.;
  // west readout
  if ( 0 == layerid % 2 )
  {
    // inner chamber
    if ( layerid < 8 ) { vp = vinner; }
    else { vp = vouter; }
    tp = wireLen * fabs( zhit - bPoint.z() ) / wireZLen / vp;
  }
 
  // east readout
  if ( 1 == layerid % 2 )
  {
    // inner chamber
    if ( layerid < 8 ) { vp = vinner; }
    else { vp = vouter; }
    tp = wireLen * fabs( zhit - fPoint.z() ) / wireZLen / vp;
  }
 
  // s to ns
  tp = 1.0e9 * tp;
 
  if ( 0 == tprop_ ) tp = 0.;
 
  // std::cout<<"propogation time: "<<tp<<std::endl;
 
  int    wireid = hitmdc.wire().geoID();
  double drifttime1( 0. );
  double drifttime2( 0. );
  double drifttime3( 0. );
 
  MdcDigiCol::iterator iter = mdcDigiCol_->begin();
  for ( ; iter != mdcDigiCol_->end(); iter++ )
  {
    if ( ( *iter )->identify() == ( hitmdc.rechitptr() )->getMdcId() ) break;
  }
  // double t0 = get_T0();
  //  see the code of wulh in /Mdc/MdcCalibFunSvc/MdcCalibFunSvc/MdcCalibFunSvc.h,
  //  double getT0(int wireid) const { return m_t0[wireid]; }
 
  // double getTimeWalk(int layid, double Q) const ;
  double Q        = RawDataUtil::MdcCharge( ( *iter )->getChargeChannel() );
  double timewalk = CalibFunSvc_->getTimeWalk( layerid, Q );
 
  if ( debug_ == 4 )
  { std::cout << "CalibFunSvc_->getTimeWalk, timewalk =" << timewalk << std::endl; }
 
  double timeoffset = CalibFunSvc_->getT0( wireid );
  if ( debug_ == 4 )
  { std::cout << "CalibFunSvc_->getT0, timeoffset =" << timeoffset << std::endl; }
 
  double eventt0 = getT0();
 
  if ( debug_ == 4 )
  {
    std::cout << "the Event T0 we get in the function getDriftTime(...) is " << eventt0
              << std::endl;
  }
 
  // this tdc value come from MdcRecHit assigned by zhangyao
  double tdctime1 = hitmdc.tdc();
  double tdctime2( 0. );
  double tdctime3( 0. );
 
  if ( debug_ == 4 ) { std::cout << "tdctime1 be here is .." << tdctime1 << std::endl; }
 
  // this tdc value come from MdcDigiCol time channel
  // attention, if we use the iter like this: for(MdcDigiCol::iterator iter =
  // mdcDigiCol_->begin(); iter != mdcDigiCol_->end(); iter++) it cannot pass the gmake , throw
  // an error !
  if ( debug_ == 4 )
  {
    //  std::cout<<"the size of the mdcDigiCol_ is "<<mdcDigiCol_.size()<<std::endl;
  }
  // MdcDigiCol::iterator iter = mdcDigiCol_->begin();
  // for(; iter != mdcDigiCol_->end(); iter++ ) {
  //       if((*iter)->identify() == (hitmdc.rechitptr())->getMdcId()) break;
  //  }
  if ( debug_ == 4 )
  {
    std::cout << "the time channel be here is .." << ( *iter )->getTimeChannel() << std::endl;
  }
  tdctime2   = RawDataUtil::MdcTime( ( *iter )->getTimeChannel() );
  tdctime3   = hitmdc.rechitptr()->getTdc();
  drifttime1 = tdctime1 - eventt0 - toftime - timewalk - timeoffset - tp;
  drifttime2 = tdctime2 - eventt0 - toftime - timewalk - timeoffset - tp;
  drifttime3 = tdctime3 - eventt0 - toftime - timewalk - timeoffset - tp;
  if ( debug_ == 4 )
  {
    std::cout << "we now compare the three kind of tdc , the tdc get from timeChannel() is "
              << tdctime2 << " the tdc get from KalFitHitMdc is " << tdctime1
              << " the tdc from MdcRecHit is " << tdctime3 << " the drifttime1 is ..."
              << drifttime1 << " the drifttime 2 is ..." << drifttime2
              << " the drifttime3 is ..." << drifttime3 << std::endl;
  }
  // return drifttime3;
  // return drifttime1;
  if ( drifttime_choice_ == 0 ) return drifttime2;
  if ( drifttime_choice_ == 1 )
    // use the driftT caluculated by track-finding
    return hitmdc.rechitptr()->getDriftT();
 
  std::cerr << __FILE__ << ":" << __LINE__ << ": should not reach here" << std::endl;
  exit( 1 );
}

Referenced by chi2_next(), and chi2_next().

getDriftTime() [2/3]

double KalFitTrack::getDriftTime	(	KalFitHitMdc &	hitmdc,
		double	toftime ) const

getDriftTime() [3/3]

double KalFitTrack::getDriftTime	(	KalFitHitMdc &	hitmdc,
		double	toftime ) const

getFiTerm() [1/3]

double KalFitTrack::getFiTerm ( void )

inline

Definition at line 155 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return fiTerm_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

getFiTerm() [2/3]

double KalFitTrack::getFiTerm ( void )

inline

Definition at line 155 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return fiTerm_; }

getFiTerm() [3/3]

double KalFitTrack::getFiTerm ( void )

inline

Definition at line 155 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return fiTerm_; }

getInitMatrix() [1/3]

HepSymMatrix KalFitTrack::getInitMatrix

(

void

)

const

Definition at line 37 of file KalFitTrack2.cxx.

{ return initMatrix_; }

Referenced by KalFitAlg::smoother_calib().

getInitMatrix() [2/3]

HepSymMatrix KalFitTrack::getInitMatrix

(

void

)

const

getInitMatrix() [3/3]

HepSymMatrix KalFitTrack::getInitMatrix

(

void

)

const

getPathSM() [1/3]

double KalFitTrack::getPathSM ( void )

inline

Definition at line 149 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pathSM_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

getPathSM() [2/3]

double KalFitTrack::getPathSM ( void )

inline

Definition at line 149 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pathSM_; }

getPathSM() [3/3]

double KalFitTrack::getPathSM ( void )

inline

Definition at line 149 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pathSM_; }

getSigma() [1/6]

double KalFitTrack::getSigma	(	int	layerId,
		double	driftDist ) const

Definition at line 3235 of file KalFitTrack.cxx.

                                                                  {
  double sigma1, sigma2, f;
  driftDist *= 10; // mm
  if ( layerId < 8 )
  {
    if ( driftDist < 0.5 )
    {
      sigma1 = 0.112784;
      sigma2 = 0.229274;
      f      = 0.666;
    }
    else if ( driftDist < 1. )
    {
      sigma1 = 0.103123;
      sigma2 = 0.269797;
      f      = 0.934;
    }
    else if ( driftDist < 1.5 )
    {
      sigma1 = 0.08276;
      sigma2 = 0.17493;
      f      = 0.89;
    }
    else if ( driftDist < 2. )
    {
      sigma1 = 0.070109;
      sigma2 = 0.149859;
      f      = 0.89;
    }
    else if ( driftDist < 2.5 )
    {
      sigma1 = 0.064453;
      sigma2 = 0.130149;
      f      = 0.886;
    }
    else if ( driftDist < 3. )
    {
      sigma1 = 0.062383;
      sigma2 = 0.138806;
      f      = 0.942;
    }
    else if ( driftDist < 3.5 )
    {
      sigma1 = 0.061873;
      sigma2 = 0.145696;
      f      = 0.946;
    }
    else if ( driftDist < 4. )
    {
      sigma1 = 0.061236;
      sigma2 = 0.119584;
      f      = 0.891;
    }
    else if ( driftDist < 4.5 )
    {
      sigma1 = 0.066292;
      sigma2 = 0.148426;
      f      = 0.917;
    }
    else if ( driftDist < 5. )
    {
      sigma1 = 0.078074;
      sigma2 = 0.188148;
      f      = 0.911;
    }
    else if ( driftDist < 5.5 )
    {
      sigma1 = 0.088657;
      sigma2 = 0.27548;
      f      = 0.838;
    }
    else
    {
      sigma1 = 0.093089;
      sigma2 = 0.115556;
      f      = 0.367;
    }
  }
  else
  {
    if ( driftDist < 0.5 )
    {
      sigma1 = 0.112433;
      sigma2 = 0.327548;
      f      = 0.645;
    }
    else if ( driftDist < 1. )
    {
      sigma1 = 0.096703;
      sigma2 = 0.305206;
      f      = 0.897;
    }
    else if ( driftDist < 1.5 )
    {
      sigma1 = 0.082518;
      sigma2 = 0.248913;
      f      = 0.934;
    }
    else if ( driftDist < 2. )
    {
      sigma1 = 0.072501;
      sigma2 = 0.153868;
      f      = 0.899;
    }
    else if ( driftDist < 2.5 )
    {
      sigma1 = 0.065535;
      sigma2 = 0.14246;
      f      = 0.914;
    }
    else if ( driftDist < 3. )
    {
      sigma1 = 0.060497;
      sigma2 = 0.126489;
      f      = 0.918;
    }
    else if ( driftDist < 3.5 )
    {
      sigma1 = 0.057643;
      sigma2 = 0.112927;
      f      = 0.892;
    }
    else if ( driftDist < 4. )
    {
      sigma1 = 0.055266;
      sigma2 = 0.094833;
      f      = 0.887;
    }
    else if ( driftDist < 4.5 )
    {
      sigma1 = 0.056263;
      sigma2 = 0.124419;
      f      = 0.932;
    }
    else if ( driftDist < 5. )
    {
      sigma1 = 0.056599;
      sigma2 = 0.124248;
      f      = 0.923;
    }
    else if ( driftDist < 5.5 )
    {
      sigma1 = 0.061377;
      sigma2 = 0.146147;
      f      = 0.964;
    }
    else if ( driftDist < 6. )
    {
      sigma1 = 0.063978;
      sigma2 = 0.150591;
      f      = 0.942;
    }
    else if ( driftDist < 6.5 )
    {
      sigma1 = 0.072951;
      sigma2 = 0.15685;
      f      = 0.913;
    }
    else if ( driftDist < 7. )
    {
      sigma1 = 0.085438;
      sigma2 = 0.255109;
      f      = 0.931;
    }
    else if ( driftDist < 7.5 )
    {
      sigma1 = 0.101635;
      sigma2 = 0.315529;
      f      = 0.878;
    }
    else
    {
      sigma1 = 0.149529;
      sigma2 = 0.374697;
      f      = 0.89;
    }
  }
  double sigmax = sqrt( f * sigma1 * sigma1 + ( 1 - f ) * sigma2 * sigma2 ) * 0.1;
  return sigmax; // cm
}

Referenced by chi2_next(), and chi2_next().

getSigma() [2/6]

double KalFitTrack::getSigma	(	int	layerId,
		double	driftDist ) const

getSigma() [3/6]

double KalFitTrack::getSigma	(	int	layerId,
		double	driftDist ) const

getSigma() [4/6]

double KalFitTrack::getSigma	(	KalFitHitMdc &	hitmdc,
		double	tanlam,
		int	lr,
		double	dist ) const

Definition at line 198 of file KalFitTrack2.cxx.

                                                  {
  int    layerid   = hitmdc.wire().layer().layerId();
  double entrangle = hitmdc.rechitptr()->getEntra();
  //  double tanlam = hitmdc.rechitptr()->getTanl();
  double z = hitmdc.rechitptr()->getZhit();
  double Q = hitmdc.rechitptr()->getAdc();
  // std::cout<<" the driftdist  before getsigma is "<<dist<<" the layer
  // is"<<layerid<<std::endl; cout<<"layerid, lr, dist, entrangle, tanlam, z , Q =
  // "<<layerid<<", "<<lr<<", "<<dist<<", "<<entrangle<<",
  // "<<tanlam<<", "<<z<<", "<<Q<<endl;//wangll
  double temp = CalibFunSvc_->getSigma( layerid, lr, dist, entrangle, tanlam, z, Q );
  // std::cout<<" the sigma is "<<temp<<std::endl;
  return temp;
}

getSigma() [5/6]

double KalFitTrack::getSigma	(	KalFitHitMdc &	hitmdc,
		double	tanlam,
		int	lr,
		double	dist ) const

getSigma() [6/6]

double KalFitTrack::getSigma	(	KalFitHitMdc &	hitmdc,
		double	tanlam,
		int	lr,
		double	dist ) const

getT0() [1/3]

double KalFitTrack::getT0

(

void

)

const

Definition at line 47 of file KalFitTrack2.cxx.

                                      {
  //------------------get  event  start time-----------
 
  if ( debug_ == 4 )
  { std::cout << "in function KalFitTrack::getT0 ( ), EventT0_ = " << EventT0_ << std::endl; }
  return EventT0_;
}

Referenced by getDriftTime().

getT0() [2/3]

double KalFitTrack::getT0

(

void

)

const

getT0() [3/3]

double KalFitTrack::getT0

(

void

)

const

getTofSM() [1/3]

double KalFitTrack::getTofSM ( void )

inline

Definition at line 152 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return tof2_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

getTofSM() [2/3]

double KalFitTrack::getTofSM ( void )

inline

Definition at line 152 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return tof2_; }

getTofSM() [3/3]

double KalFitTrack::getTofSM ( void )

inline

Definition at line 152 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return tof2_; }

HelixSeg() [1/3]

KalFitHelixSeg & KalFitTrack::HelixSeg ( int i )

inline

Definition at line 225 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_[i]; }

Referenced by KalFitAlg::smoother_calib(), update_hits(), and update_hits_csmalign().

HelixSeg() [2/3]

KalFitHelixSeg & KalFitTrack::HelixSeg ( int i )

inline

Definition at line 225 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_[i]; }

HelixSeg() [3/3]

KalFitHelixSeg & KalFitTrack::HelixSeg ( int i )

inline

Definition at line 225 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_[i]; }

HelixSegs() [1/6]

void KalFitTrack::HelixSegs

(

vector< KalFitHelixSeg > &

vs

)

Referenced by KalFitAlg::complete_track(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), and KalFitAlg::smoother_calib().

HelixSegs() [2/6]

void KalFitTrack::HelixSegs

(

vector< KalFitHelixSeg > &

vs

)

HelixSegs() [3/6]

void KalFitTrack::HelixSegs

(

vector< KalFitHelixSeg > &

vs

)

HelixSegs() [4/6]

vector< KalFitHelixSeg > & KalFitTrack::HelixSegs ( void )

inline

Definition at line 224 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_; }

HelixSegs() [5/6]

vector< KalFitHelixSeg > & KalFitTrack::HelixSegs ( void )

inline

Definition at line 224 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_; }

HelixSegs() [6/6]

vector< KalFitHelixSeg > & KalFitTrack::HelixSegs ( void )

inline

Definition at line 224 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return HelixSegs_; }

HitMdc() [1/3]

KalFitHitMdc & KalFitTrack::HitMdc ( int i )

inline

Definition at line 220 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_[i]; }

Referenced by chi2_next(), chi2_next(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), smoother_Mdc(), KalFitAlg::start_seed(), update_hits(), and update_hits_csmalign().

HitMdc() [2/3]

KalFitHitMdc & KalFitTrack::HitMdc ( int i )

inline

Definition at line 220 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_[i]; }

HitMdc() [3/3]

KalFitHitMdc & KalFitTrack::HitMdc ( int i )

inline

Definition at line 220 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_[i]; }

HitsMdc() [1/6]

void KalFitTrack::HitsMdc

(

vector< KalFitHitMdc > &

lh

)

Definition at line 1873 of file KalFitTrack.cxx.

{ HitsMdc_ = lh; }

Referenced by KalFitAlg::complete_track(), KalFitAlg::complete_track(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew(), KalFitAlg::smoother_anal(), and KalFitAlg::start_seed().

HitsMdc() [2/6]

void KalFitTrack::HitsMdc

(

vector< KalFitHitMdc > &

lh

)

HitsMdc() [3/6]

void KalFitTrack::HitsMdc

(

vector< KalFitHitMdc > &

lh

)

HitsMdc() [4/6]

vector< KalFitHitMdc > & KalFitTrack::HitsMdc ( void )

inline

Definition at line 219 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_; }

Referenced by order_hits().

HitsMdc() [5/6]

vector< KalFitHitMdc > & KalFitTrack::HitsMdc ( void )

inline

Definition at line 219 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_; }

HitsMdc() [6/6]

vector< KalFitHitMdc > & KalFitTrack::HitsMdc ( void )

inline

Definition at line 219 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return HitsMdc_; }

insist() [1/6]

void KalFitTrack::insist ( int t )

inline

Definition at line 199 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ insist_ = t; }

insist() [2/6]

void KalFitTrack::insist ( int t )

inline

Definition at line 199 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ insist_ = t; }

insist() [3/6]

void KalFitTrack::insist ( int t )

inline

Definition at line 199 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ insist_ = t; }

insist() [4/6]

int KalFitTrack::insist ( void ) const

inline

Extractor :

Definition at line 168 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return insist_; }

insist() [5/6]

int KalFitTrack::insist ( void ) const

inline

Extractor :

Definition at line 168 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return insist_; }

insist() [6/6]

int KalFitTrack::insist ( void ) const

inline

Extractor :

Definition at line 168 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return insist_; }

intersect_cylinder() [1/3]

double KalFitTrack::intersect_cylinder

(

double

r

)

const

Intersection with different geometry.

Definition at line 127 of file KalFitTrack.cxx.

                                                       {
  double m_rad = radius();
  double l     = center().perp();
 
  double cosPhi = ( m_rad * m_rad + l * l - r * r ) / ( 2 * m_rad * l );
 
  if ( cosPhi < -1 || cosPhi > 1 ) return 0;
 
  double dPhi = center().phi() - acos( cosPhi ) - phi0();
 
  if ( dPhi < -M_PI ) dPhi += 2 * M_PI;
 
  return dPhi;
}

Referenced by KalFitCylinder::intersect(), KalFitCylinder::intersect(), intersect_yz_plane(), intersect_zx_plane(), and VertexExtrapolate::KalFitExt().

intersect_cylinder() [2/3]

double KalFitTrack::intersect_cylinder

(

double

r

)

const

Intersection with different geometry.

intersect_cylinder() [3/3]

double KalFitTrack::intersect_cylinder

(

double

r

)

const

Intersection with different geometry.

intersect_xy_plane() [1/3]

double KalFitTrack::intersect_xy_plane

(

double

z

)

const

Definition at line 172 of file KalFitTrack.cxx.

                                                       {
  if ( tanl() != 0 && radius() != 0 )
    return ( pivot().z() + dz() - z ) / ( radius() * tanl() );
  else return 0;
}

Referenced by KalFitCylinder::intersect(), and KalFitCylinder::intersect().

intersect_xy_plane() [2/3]

double KalFitTrack::intersect_xy_plane

(

double

z

)

const

intersect_xy_plane() [3/3]

double KalFitTrack::intersect_xy_plane

(

double

z

)

const

intersect_yz_plane() [1/3]

double KalFitTrack::intersect_yz_plane	(	const HepTransform3D &	plane,
		double	x ) const

Definition at line 157 of file KalFitTrack.cxx.

                                                                                    {
  HepPoint3D xc = plane * center();
  double     r  = radius();
  double     d  = r * r - ( x - xc.x() ) * ( x - xc.x() );
  if ( d < 0 ) return 0;
 
  double yy = xc.y();
  if ( yy > 0 ) yy -= sqrt( d );
  else yy += sqrt( d );
 
  double l = ( plane.inverse() * HepPoint3D( x, yy, 0 ) ).perp();
 
  return intersect_cylinder( l );
}

intersect_yz_plane() [2/3]

double KalFitTrack::intersect_yz_plane	(	const HepTransform3D &	plane,
		double	x ) const

intersect_yz_plane() [3/3]

double KalFitTrack::intersect_yz_plane	(	const HepTransform3D &	plane,
		double	x ) const

intersect_zx_plane() [1/3]

double KalFitTrack::intersect_zx_plane	(	const HepTransform3D &	plane,
		double	y ) const

Definition at line 142 of file KalFitTrack.cxx.

                                                                                    {
  HepPoint3D xc = plane * center();
  double     r  = radius();
  double     d  = r * r - ( y - xc.y() ) * ( y - xc.y() );
  if ( d < 0 ) return 0;
 
  double xx = xc.x();
  if ( xx > 0 ) xx -= sqrt( d );
  else xx += sqrt( d );
 
  double l = ( plane.inverse() * HepPoint3D( xx, y, 0 ) ).perp();
 
  return intersect_cylinder( l );
}

intersect_zx_plane() [2/3]

double KalFitTrack::intersect_zx_plane	(	const HepTransform3D &	plane,
		double	y ) const

intersect_zx_plane() [3/3]

double KalFitTrack::intersect_zx_plane	(	const HepTransform3D &	plane,
		double	y ) const

lead() [1/6]

void KalFitTrack::lead ( int i )

static

Definition at line 1864 of file KalFitTrack.cxx.

{ lead_ = i; }

lead() [2/6]

void KalFitTrack::lead ( int i )

static

lead() [3/6]

void KalFitTrack::lead ( int i )

static

lead() [4/6]

int KalFitTrack::lead ( void )

static

Magnetic field map.

Definition at line 1865 of file KalFitTrack.cxx.

{ return lead_; }

Referenced by KalFitAlg::initialize().

lead() [5/6]

int KalFitTrack::lead ( void )

static

Magnetic field map.

lead() [6/6]

int KalFitTrack::lead ( void )

static

Magnetic field map.

LR() [1/3]

void KalFitTrack::LR ( int x )

static

Definition at line 1872 of file KalFitTrack.cxx.

{ LR_ = x; }

Referenced by KalFitAlg::initialize().

LR() [2/3]

void KalFitTrack::LR ( int x )

static

LR() [3/3]

void KalFitTrack::LR ( int x )

static

mass() [1/6]

double KalFitTrack::mass ( int i )

static

Definition at line 1858 of file KalFitTrack.cxx.

{ return MASS[i]; }

mass() [2/6]

double KalFitTrack::mass ( int i )

static

mass() [3/6]

double KalFitTrack::mass ( int i )

static

mass() [4/6]

double KalFitTrack::mass ( void ) const

inline

Definition at line 172 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return mass_; }

Referenced by KalFitAlg::fillTds(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_lead(), and KalFitAlg::smoother_anal().

mass() [5/6]

double KalFitTrack::mass ( void ) const

inline

Definition at line 172 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return mass_; }

mass() [6/6]

double KalFitTrack::mass ( void ) const

inline

Definition at line 172 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return mass_; }

mom() [1/3]

CLHEP::Hep3Vector * KalFitTrack::mom ( void )

inline

Definition at line 180 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return mom_; }

mom() [2/3]

CLHEP::Hep3Vector * KalFitTrack::mom ( void )

inline

Definition at line 180 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return mom_; }

mom() [3/3]

CLHEP::Hep3Vector * KalFitTrack::mom ( void )

inline

Definition at line 180 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return mom_; }

ms() [1/3]

void KalFitTrack::ms	(	double	path,
		const KalFitMaterial &	m,
		int	index )

Definition at line 178 of file KalFitTrack.cxx.

                                                                             {
  HepSymMatrix ea = Ea();
  // cout<<"ms():path  "<<path<<endl;
  // cout<<"ms():ea before: "<<ea<<endl;
  double k   = kappa();
  double t   = tanl();
  double t2  = t * t;
  double pt2 = 1 + t2;
  double k2  = k * k;
 
  double pmag    = 1 / fabs( k ) * sqrt( pt2 );
  double dth     = material.mcs_angle( mass_, path, pmag );
  double dth2    = dth * dth;
  double pt2dth2 = pt2 * dth2;
 
  ea[1][1] += pt2dth2;
  ea[2][2] += k2 * t2 * dth2;
  ea[2][4] += k * t * pt2dth2;
  ea[4][4] += pt2 * pt2dth2;
 
  ea[3][3] += dth2 * path * path / 3 / ( 1 + t2 );
  ea[3][4] += dth2 * path / 2 * sqrt( 1 + t2 );
  ea[3][2] += dth2 * t / sqrt( 1 + t2 ) * k * path / 2;
  ea[0][0] += dth2 * path * path / 3;
  ea[0][1] += dth2 * sqrt( 1 + t2 ) * path / 2;
 
  Ea( ea );
  // cout<<"ms():ms angle in this: "<<dth<<endl;
  // cout<<"ms():ea after: "<<Ea()<<endl;
  if ( index < 0 )
  {
    double x0 = material.X0();
    if ( x0 ) path_rd_ += path / x0;
  }
}

Referenced by KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

ms() [2/3]

void KalFitTrack::ms	(	double	path,
		const KalFitMaterial &	m,
		int	index )

ms() [3/3]

void KalFitTrack::ms	(	double	path,
		const KalFitMaterial &	m,
		int	index )

msgasmdc() [1/3]

void KalFitTrack::msgasmdc	(	double	path,
		int	index )

Calculate multiple scattering angle.

Definition at line 214 of file KalFitTrack.cxx.

                                                   {
  double k  = kappa();
  double t  = tanl();
  double t2 = t * t;
  double k2 = k * k;
 
  double pmag = ( 1 / fabs( k ) ) * sqrt( 1 + t2 );
  double psq  = pmag * pmag;
  /*
     double Zprims = 3/2*0.076 + 0.580/9.79*4.99*(4.99+1) +
     0.041/183.85*74*(74+1) + 0.302/26.98 * 13 * (13+1);
     double chicc = 0.00039612 * sqrt(Zprims * 0.001168);
     double dth = 2.557 * chicc * sqrt(path * (mass_*mass_ + psq)) / psq;
   */
 
  // std::cout<<" mdcGasRadlen: "<<mdcGasRadlen_<<std::endl;
  double pathx = path / mdcGasRadlen_;
  double dth =
      0.0136 * sqrt( pathx * ( mass_ * mass_ + psq ) ) / psq * ( 1 + 0.038 * log( pathx ) );
  ;
  HepSymMatrix ea = Ea();
#ifdef YDEBUG
  cout << "msgasmdc():path  " << path << "  pathx " << pathx << endl;
  cout << "msgasmdc():dth  " << dth << endl;
  cout << "msgasmdc():ea before: " << ea << endl;
#endif
  double dth2 = dth * dth;
 
  ea[1][1] += ( 1 + t2 ) * dth2;
  ea[2][2] += k2 * t2 * dth2;
  ea[2][4] += k * t * ( 1 + t2 ) * dth2;
  ea[4][4] += ( 1 + t2 ) * ( 1 + t2 ) * dth2;
 
  // additionnal terms (terms proportional to l and l^2)
 
  ea[3][3] += dth2 * path * path / 3 / ( 1 + t2 );
  ea[3][4] += dth2 * path / 2 * sqrt( 1 + t2 );
  ea[3][2] += dth2 * t / sqrt( 1 + t2 ) * k * path / 2;
  ea[0][0] += dth2 * path * path / 3;
  ea[0][1] += dth2 * sqrt( 1 + t2 ) * path / 2;
 
  Ea( ea );
#ifdef YDEBUG
  cout << "msgasmdc():ea after: " << Ea() << endl;
#endif
  if ( index < 0 )
  {
    pathip_ += path;
    // RMK : put by hand, take care !!
    double x0 = mdcGasRadlen_; // for the Mdc gas
    path_rd_ += path / x0;
    tof( path );
#ifdef YDEBUG
    cout << "ms...pathip..." << pathip_ << endl;
#endif
  }
}

Referenced by KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), and KalFitAlg::smoother_calib().

msgasmdc() [2/3]

void KalFitTrack::msgasmdc	(	double	path,
		int	index )

Calculate multiple scattering angle.

msgasmdc() [3/3]

void KalFitTrack::msgasmdc	(	double	path,
		int	index )

Calculate multiple scattering angle.

ncath() [1/3]

unsigned int KalFitTrack::ncath ( void ) const

inline

Definition at line 190 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return ncath_; }

ncath() [2/3]

unsigned int KalFitTrack::ncath ( void ) const

inline

Definition at line 190 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return ncath_; }

ncath() [3/3]

unsigned int KalFitTrack::ncath ( void ) const

inline

Definition at line 190 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return ncath_; }

nchits() [1/6]

void KalFitTrack::nchits ( int n )

inline

Definition at line 207 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ nchits_ = n; }

nchits() [2/6]

void KalFitTrack::nchits ( int n )

inline

Definition at line 207 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ nchits_ = n; }

nchits() [3/6]

void KalFitTrack::nchits ( int n )

inline

Definition at line 207 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ nchits_ = n; }

nchits() [4/6]

unsigned int KalFitTrack::nchits ( void ) const

inline

Definition at line 188 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return nchits_; }

Referenced by KalFitAlg::fillTds(), KalFitAlg::fillTds_ip(), KalFitAlg::fillTds_lead(), and KalFitAlg::start_seed().

nchits() [5/6]

unsigned int KalFitTrack::nchits ( void ) const

inline

Definition at line 188 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return nchits_; }

nchits() [6/6]

unsigned int KalFitTrack::nchits ( void ) const

inline

Definition at line 188 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return nchits_; }

ndf_back() [1/6]

void KalFitTrack::ndf_back ( int n )

inline

Definition at line 206 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ ndf_back_ = n; }

ndf_back() [2/6]

void KalFitTrack::ndf_back ( int n )

inline

Definition at line 206 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ ndf_back_ = n; }

ndf_back() [3/6]

void KalFitTrack::ndf_back ( int n )

inline

Definition at line 206 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ ndf_back_ = n; }

ndf_back() [4/6]

int KalFitTrack::ndf_back ( void ) const

inline

Definition at line 175 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return ndf_back_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::start_seed().

ndf_back() [5/6]

int KalFitTrack::ndf_back ( void ) const

inline

Definition at line 175 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return ndf_back_; }

ndf_back() [6/6]

int KalFitTrack::ndf_back ( void ) const

inline

Definition at line 175 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return ndf_back_; }

nhit_r() [1/3]

int KalFitTrack::nhit_r ( void ) const

inline

Definition at line 193 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return nhit_r_; }

nhit_r() [2/3]

int KalFitTrack::nhit_r ( void ) const

inline

Definition at line 193 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return nhit_r_; }

nhit_r() [3/3]

int KalFitTrack::nhit_r ( void ) const

inline

Definition at line 193 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return nhit_r_; }

nhit_z() [1/3]

int KalFitTrack::nhit_z ( void ) const

inline

Definition at line 194 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return nhit_z_; }

nhit_z() [2/3]

int KalFitTrack::nhit_z ( void ) const

inline

Definition at line 194 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return nhit_z_; }

nhit_z() [3/3]

int KalFitTrack::nhit_z ( void ) const

inline

Definition at line 194 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return nhit_z_; }

nLayerUsed() [1/3]

int KalFitTrack::nLayerUsed ( )

inline

Definition at line 326 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

                   {
    int n = 0;
    for ( int i = 0; i < 43; i++ ) n += myLayerUsed[i];
    return n;
  }

Referenced by KalFitAlg::fillTds_ip().

nLayerUsed() [2/3]

int KalFitTrack::nLayerUsed ( )

inline

Definition at line 326 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

                   {
    int n = 0;
    for ( int i = 0; i < 43; i++ ) n += myLayerUsed[i];
    return n;
  }

nLayerUsed() [3/3]

int KalFitTrack::nLayerUsed ( )

inline

Definition at line 326 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

                   {
    int n = 0;
    for ( int i = 0; i < 43; i++ ) n += myLayerUsed[i];
    return n;
  }

nmass() [1/3]

int KalFitTrack::nmass ( void )

static

Definition at line 1857 of file KalFitTrack.cxx.

{ return NMASS; }

Referenced by KalFitAlg::complete_track(), and KalFitAlg::complete_track().

nmass() [2/3]

int KalFitTrack::nmass ( void )

static

nmass() [3/3]

int KalFitTrack::nmass ( void )

static

nster() [1/6]

void KalFitTrack::nster ( int n )

inline

Definition at line 208 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ nster_ = n; }

nster() [2/6]

void KalFitTrack::nster ( int n )

inline

Definition at line 208 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ nster_ = n; }

nster() [3/6]

void KalFitTrack::nster ( int n )

inline

Definition at line 208 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ nster_ = n; }

nster() [4/6]

unsigned int KalFitTrack::nster ( void ) const

inline

Definition at line 189 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return nster_; }

Referenced by KalFitAlg::fillTds(), KalFitAlg::fillTds_ip(), KalFitAlg::fillTds_lead(), and KalFitAlg::start_seed().

nster() [5/6]

unsigned int KalFitTrack::nster ( void ) const

inline

Definition at line 189 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return nster_; }

nster() [6/6]

unsigned int KalFitTrack::nster ( void ) const

inline

Definition at line 189 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return nster_; }

number_wirhit() [1/3]

void KalFitTrack::number_wirhit

(

void

)

Definition at line 449 of file KalFitTrack.cxx.

                                      {
  unsigned int nhit    = HitsMdc_.size();
  int          Num[50] = { 0 };
  for ( unsigned i = 0; i < nhit; i++ ) Num[HitsMdc_[i].wire().layer().layerId()]++;
  for ( unsigned i = 0; i < nhit; i++ )
    if ( Num[HitsMdc_[i].wire().layer().layerId()] > 2 ) HitsMdc_[i].chi2( -2 );
}

number_wirhit() [2/3]

void KalFitTrack::number_wirhit

(

void

)

number_wirhit() [3/3]

void KalFitTrack::number_wirhit

(

void

)

numf() [1/6]

void KalFitTrack::numf ( int i )

static

Definition at line 1870 of file KalFitTrack.cxx.

{ numf_ = i; }

numf() [2/6]

void KalFitTrack::numf ( int i )

static

numf() [3/6]

void KalFitTrack::numf ( int i )

static

numf() [4/6]

int KalFitTrack::numf ( void )

static

Definition at line 1871 of file KalFitTrack.cxx.

{ return numf_; }

Referenced by KalFitAlg::initialize().

numf() [5/6]

int KalFitTrack::numf ( void )

static

numf() [6/6]

int KalFitTrack::numf ( void )

static

order_hits() [1/3]

void KalFitTrack::order_hits

(

void

)

Definition at line 434 of file KalFitTrack.cxx.

                                   {
  for ( int it = 0; it < HitsMdc().size() - 1; it++ )
  {
    if ( HitsMdc_[it].wire().layer().layerId() == HitsMdc_[it + 1].wire().layer().layerId() )
    {
      if ( ( kappa() < 0 ) &&
           ( HitsMdc_[it].wire().localId() > HitsMdc_[it + 1].wire().localId() ) )
      { std::swap( HitsMdc_[it], HitsMdc_[it + 1] ); }
      if ( ( kappa() > 0 ) &&
           ( HitsMdc_[it].wire().localId() < HitsMdc_[it + 1].wire().localId() ) )
      { std::swap( HitsMdc_[it], HitsMdc_[it + 1] ); }
    }
  }
}

order_hits() [2/3]

void KalFitTrack::order_hits

(

void

)

order_hits() [3/3]

void KalFitTrack::order_hits

(

void

)

order_wirhit() [1/3]

void KalFitTrack::order_wirhit

(

int

index

)

Modifier Order the wire hits for mdc track

Definition at line 367 of file KalFitTrack.cxx.

                                          {
  unsigned int nhit      = HitsMdc_.size();
  Helix        tracktest = *(Helix*)this;
  int          ind       = 0;
  double*      Rad       = new double[nhit];
  double*      Ypos      = new double[nhit];
  for ( unsigned i = 0; i < nhit; i++ )
  {
 
    HepPoint3D fwd( HitsMdc_[i].wire().fwd() );
    HepPoint3D bck( HitsMdc_[i].wire().bck() );
    Hep3Vector wire = (CLHEP::Hep3Vector)fwd - (CLHEP::Hep3Vector)bck;
 
    // Modification for stereo wires :
    Helix work = tracktest;
    work.ignoreErrorMatrix();
    work.pivot( ( fwd + bck ) * .5 );
    HepPoint3D x0 = ( work.x( 0 ).z() - bck.z() ) / wire.z() * wire + bck;
 
    tracktest.pivot( x0 );
    Rad[ind]  = tracktest.x( 0 ).perp();
    Ypos[ind] = x0.y();
    ind++;
    // cout<<"Ypos: "<<Ypos[ind-1]<<endl;
  }
 
  // Reorder...
  if ( index < 0 )
    for ( int j, k = nhit - 1; k >= 0; k = j )
    {
      j = -1;
      for ( int i = 1; i <= k; i++ )
        if ( Rad[i - 1] > Rad[i] )
        {
          j = i - 1;
          std::swap( Rad[i], Rad[j] );
          std::swap( HitsMdc_[i], HitsMdc_[j] );
        }
    }
  if ( index > 0 )
    for ( int j, k = nhit - 1; k >= 0; k = j )
    {
      j = -1;
      for ( int i = 1; i <= k; i++ )
        if ( Rad[i - 1] < Rad[i] )
        {
          j = i - 1;
          std::swap( Rad[i], Rad[j] );
          std::swap( HitsMdc_[i], HitsMdc_[j] );
        }
    }
  if ( index == 0 )
    for ( int j, k = nhit - 1; k >= 0; k = j )
    {
      j = -1;
      for ( int i = 1; i <= k; i++ )
        if ( Ypos[i - 1] > Ypos[i] )
        {
          j = i - 1;
          std::swap( Ypos[i], Ypos[j] );
          std::swap( HitsMdc_[i], HitsMdc_[j] );
        }
    }
  delete[] Rad;
  delete[] Ypos;
}

Referenced by KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), and KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew().

order_wirhit() [2/3]

void KalFitTrack::order_wirhit

(

int

index

)

Modifier Order the wire hits for mdc track

order_wirhit() [3/3]

void KalFitTrack::order_wirhit

(

int

index

)

Modifier Order the wire hits for mdc track

p_kaon() [1/6]

void KalFitTrack::p_kaon ( double pl )

inline

Definition at line 202 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ p_kaon_ = pl; }

p_kaon() [2/6]

void KalFitTrack::p_kaon ( double pl )

inline

Definition at line 202 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ p_kaon_ = pl; }

p_kaon() [3/6]

void KalFitTrack::p_kaon ( double pl )

inline

Definition at line 202 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ p_kaon_ = pl; }

p_kaon() [4/6]

double KalFitTrack::p_kaon ( void ) const

inline

Definition at line 184 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return p_kaon_; }

Referenced by KalFitAlg::complete_track(), and KalFitAlg::complete_track().

p_kaon() [5/6]

double KalFitTrack::p_kaon ( void ) const

inline

Definition at line 184 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return p_kaon_; }

p_kaon() [6/6]

double KalFitTrack::p_kaon ( void ) const

inline

Definition at line 184 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return p_kaon_; }

p_proton() [1/6]

void KalFitTrack::p_proton ( double pl )

inline

Definition at line 203 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ p_proton_ = pl; }

p_proton() [2/6]

void KalFitTrack::p_proton ( double pl )

inline

Definition at line 203 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ p_proton_ = pl; }

p_proton() [3/6]

void KalFitTrack::p_proton ( double pl )

inline

Definition at line 203 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ p_proton_ = pl; }

p_proton() [4/6]

double KalFitTrack::p_proton ( void ) const

inline

Definition at line 185 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return p_proton_; }

Referenced by KalFitAlg::complete_track(), and KalFitAlg::complete_track().

p_proton() [5/6]

double KalFitTrack::p_proton ( void ) const

inline

Definition at line 185 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return p_proton_; }

p_proton() [6/6]

double KalFitTrack::p_proton ( void ) const

inline

Definition at line 185 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return p_proton_; }

pat1() [1/3]

int KalFitTrack::pat1 ( void ) const

inline

Definition at line 191 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pat1_; }

pat1() [2/3]

int KalFitTrack::pat1 ( void ) const

inline

Definition at line 191 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pat1_; }

pat1() [3/3]

int KalFitTrack::pat1 ( void ) const

inline

Definition at line 191 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pat1_; }

pat2() [1/3]

int KalFitTrack::pat2 ( void ) const

inline

Definition at line 192 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pat2_; }

pat2() [2/3]

int KalFitTrack::pat2 ( void ) const

inline

Definition at line 192 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pat2_; }

pat2() [3/3]

int KalFitTrack::pat2 ( void ) const

inline

Definition at line 192 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pat2_; }

path_ab() [1/3]

double KalFitTrack::path_ab ( void ) const

inline

Definition at line 178 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return path_ab_; }

path_ab() [2/3]

double KalFitTrack::path_ab ( void ) const

inline

Definition at line 178 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return path_ab_; }

path_ab() [3/3]

double KalFitTrack::path_ab ( void ) const

inline

Definition at line 178 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return path_ab_; }

path_add() [1/3]

void KalFitTrack::path_add

(

double

path

)

Update the path length estimation.

Definition at line 1354 of file KalFitTrack.cxx.

                                        {
  pathip_ += path;
  tof( path );
}

path_add() [2/3]

void KalFitTrack::path_add

(

double

path

)

Update the path length estimation.

path_add() [3/3]

void KalFitTrack::path_add

(

double

path

)

Update the path length estimation.

path_rd() [1/3]

double KalFitTrack::path_rd ( void ) const

inline

Definition at line 177 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return path_rd_; }

path_rd() [2/3]

double KalFitTrack::path_rd ( void ) const

inline

Definition at line 177 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return path_rd_; }

path_rd() [3/3]

double KalFitTrack::path_rd ( void ) const

inline

Definition at line 177 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return path_rd_; }

pathip() [1/6]

void KalFitTrack::pathip ( double pl )

inline

Definition at line 201 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ pathip_ = pl; }

pathip() [2/6]

void KalFitTrack::pathip ( double pl )

inline

Definition at line 201 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ pathip_ = pl; }

pathip() [3/6]

void KalFitTrack::pathip ( double pl )

inline

Definition at line 201 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ pathip_ = pl; }

pathip() [4/6]

double KalFitTrack::pathip ( void ) const

inline

Definition at line 176 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pathip_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

pathip() [5/6]

double KalFitTrack::pathip ( void ) const

inline

Definition at line 176 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pathip_; }

pathip() [6/6]

double KalFitTrack::pathip ( void ) const

inline

Definition at line 176 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pathip_; }

PathL() [1/3]

double KalFitTrack::PathL

(

int

layer

)

Function to calculate the path length in the layer.

PathL() [2/3]

double KalFitTrack::PathL

(

int

layer

)

Function to calculate the path length in the layer.

PathL() [3/3]

double KalFitTrack::PathL

(

int

layer

)

Function to calculate the path length in the layer.

pathl() [1/3]

double * KalFitTrack::pathl ( void )

inline

Definition at line 179 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return PathL_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and pivot_numf().

pathl() [2/3]

double * KalFitTrack::pathl ( void )

inline

Definition at line 179 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return PathL_; }

pathl() [3/3]

double * KalFitTrack::pathl ( void )

inline

Definition at line 179 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return PathL_; }

pivot_forMdc() [1/3]

const HepPoint3D & KalFitTrack::pivot_forMdc ( void ) const

inline

Definition at line 141 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pivot_forMdc_; }

pivot_forMdc() [2/3]

const HepPoint3D & KalFitTrack::pivot_forMdc ( void ) const

inline

Definition at line 141 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pivot_forMdc_; }

pivot_forMdc() [3/3]

const HepPoint3D & KalFitTrack::pivot_forMdc ( void ) const

inline

Definition at line 141 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pivot_forMdc_; }

pivot_last() [1/3]

const HepPoint3D & KalFitTrack::pivot_last ( void ) const

inline

returns helix parameters

Definition at line 137 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return pivot_last_; }

pivot_last() [2/3]

const HepPoint3D & KalFitTrack::pivot_last ( void ) const

inline

returns helix parameters

Definition at line 137 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return pivot_last_; }

pivot_last() [3/3]

const HepPoint3D & KalFitTrack::pivot_last ( void ) const

inline

returns helix parameters

Definition at line 137 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return pivot_last_; }

pivot_numf() [1/6]

const HepPoint3D & KalFitTrack::pivot_numf

(

const HepPoint3D &

newPivot

)

Sets pivot position in a given mag field.

Definition at line 1428 of file KalFitTrack.cxx.

                                                                      {
 
  int        nstep( 1 );
  HepPoint3D delta_x( ( newPivot - pivot() ).x() / double( inner_steps_ ),
                      ( newPivot - pivot() ).y() / double( inner_steps_ ),
                      ( newPivot - pivot() ).z() / double( inner_steps_ ) );
  int        i = 1;
 
  while ( i <= inner_steps_ )
  {
    HepPoint3D   nextPivot( pivot() + delta_x );
    double       xnp( nextPivot.x() ), ynp( nextPivot.y() ), znp( nextPivot.z() );
    HepSymMatrix Ea_now = Ea();
    HepPoint3D   piv( pivot() );
    double       xp( piv.x() ), yp( piv.y() ), zp( piv.z() );
    double       dr    = a()[0];
    double       phi0  = a()[1];
    double       kappa = a()[2];
    double       dz    = a()[3];
    double       tanl  = a()[4];
    double       m_rad( 0 );
    if ( numfcor_ == 1 ) m_rad = radius_numf();
    else m_rad = radius();
    double rdr  = dr + m_rad;
    double phi  = fmod( phi0 + M_PI4, M_PI2 );
    double csf0 = cos( phi );
    double snf0 = ( 1. - csf0 ) * ( 1. + csf0 );
    snf0        = sqrt( ( snf0 > 0. ) ? snf0 : 0. );
    if ( phi > M_PI ) snf0 = -snf0;
 
    double xc   = xp + rdr * csf0;
    double yc   = yp + rdr * snf0;
    double csf  = ( xc - xnp ) / m_rad;
    double snf  = ( yc - ynp ) / m_rad;
    double anrm = sqrt( csf * csf + snf * snf );
    csf /= anrm;
    snf /= anrm;
    phi         = atan2( snf, csf );
    double phid = fmod( phi - phi0 + M_PI8, M_PI2 );
    if ( phid > M_PI ) phid = phid - M_PI2;
    double drp = ( xp + dr * csf0 + m_rad * ( csf0 - csf ) - xnp ) * csf +
                 ( yp + dr * snf0 + m_rad * ( snf0 - snf ) - ynp ) * snf;
    double dzp = zp + dz - m_rad * tanl * phid - znp;
 
    HepVector ap( 5 );
    ap[0] = drp;
    ap[1] = fmod( phi + M_PI4, M_PI2 );
    ap[2] = kappa;
    ap[3] = dzp;
    ap[4] = tanl;
 
    // Modification due to non uniform magnetic field :
    if ( numf_ > 10 )
    {
 
      Hep3Vector x1( xp + dr * csf0, yp + dr * snf0, zp + dz );
      double     csf0p = cos( ap[1] );
      double     snf0p = ( 1. - csf0p ) * ( 1. + csf0p );
      snf0p            = sqrt( ( snf0p > 0. ) ? snf0p : 0. );
      if ( ap[1] > M_PI ) snf0p = -snf0p;
 
      Hep3Vector  x2( xnp + drp * csf0p, ynp + drp * snf0p, znp + dzp );
      Hep3Vector  dist( ( x1 - x2 ).x() / 100.0, ( x1 - x2 ).y() / 100.0,
                        ( x1 - x2 ).z() / 100.0 );
      HepPoint3D  middlebis( ( x1.x() + x2.x() ) / 2, ( x1.y() + x2.y() ) / 2,
                             ( x1.z() + x2.z() ) / 2 );
      HepVector3D field;
      MFSvc_->fieldVector( 10. * middlebis, field );
      field = 1000. * field;
      Hep3Vector dB( field.x(), field.y(), ( field.z() - 0.1 * Bznom_ ) );
      if ( field.z() )
      {
        double    akappa( fabs( kappa ) );
        double    sign = kappa / akappa;
        HepVector dp( 3 );
        dp = 0.299792458 * sign * dB.cross( dist );
        HepVector dhp( 3 );
        dhp[0] = -akappa * ( dp[0] * csf0p + dp[1] * snf0p );
        dhp[1] = kappa * akappa * ( dp[0] * snf0p - dp[1] * csf0p );
        dhp[2] = dp[2] * akappa + dhp[1] * tanl / kappa;
        if ( numfcor_ == 0 ) { ap[1] += dhp[0]; }
        ap[2] += dhp[1];
        ap[4] += dhp[2];
      }
    }
    HepMatrix m_del = delApDelA( ap );
    Ea_now.assign( m_del * Ea_now * m_del.T() );
    pivot( nextPivot );
    a( ap );
    Ea( Ea_now );
    i++;
  }
  return newPivot;
}

Referenced by KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), KalFitAlg::smoother_calib(), KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

pivot_numf() [2/6]

const HepPoint3D & KalFitTrack::pivot_numf

(

const HepPoint3D &

newPivot

)

Sets pivot position in a given mag field.

pivot_numf() [3/6]

const HepPoint3D & KalFitTrack::pivot_numf

(

const HepPoint3D &

newPivot

)

Sets pivot position in a given mag field.

pivot_numf() [4/6]

const HepPoint3D & KalFitTrack::pivot_numf	(	const HepPoint3D &	newPivot,
		double &	pathl )

Definition at line 1523 of file KalFitTrack.cxx.

                                                                                     {
 
  Helix tracktest = *(Helix*)this;
  tracktest.ignoreErrorMatrix();
  double tl = a()[4];
  double th = 90.0 - 180.0 * M_1_PI * atan( tl );
  /*
     int nstep(1);
     if (steplev_ == 1)
     nstep = 3;
     else if (steplev_ == 2 && (th > 140 || th <45))
     if ((pivot()-newPivot).mag()<3.)
     nstep = 3;
     else
     nstep = 6;
   */
  Hep3Vector delta_x( ( newPivot - pivot() ).x() / double( outer_steps_ ),
                      ( newPivot - pivot() ).y() / double( outer_steps_ ),
                      ( newPivot - pivot() ).z() / double( outer_steps_ ) );
  int        i = 1;
 
  while ( i <= outer_steps_ )
  {
    HepPoint3D nextPivot( pivot() + delta_x );
    double     xnp( nextPivot.x() ), ynp( nextPivot.y() ), znp( nextPivot.z() );
 
    HepSymMatrix Ea_now = Ea();
    HepPoint3D   piv( pivot() );
    double       xp( piv.x() ), yp( piv.y() ), zp( piv.z() );
 
    double dr    = a()[0];
    double phi0  = a()[1];
    double kappa = a()[2];
    double dz    = a()[3];
    double tanl  = a()[4];
 
    double m_rad( 0 );
    m_rad = radius();
 
    double rdr  = dr + m_rad;
    double phi  = fmod( phi0 + M_PI4, M_PI2 );
    double csf0 = cos( phi );
    double snf0 = ( 1. - csf0 ) * ( 1. + csf0 );
    snf0        = sqrt( ( snf0 > 0. ) ? snf0 : 0. );
    if ( phi > M_PI ) snf0 = -snf0;
 
    double xc   = xp + rdr * csf0;
    double yc   = yp + rdr * snf0;
    double csf  = ( xc - xnp ) / m_rad;
    double snf  = ( yc - ynp ) / m_rad;
    double anrm = sqrt( csf * csf + snf * snf );
    csf /= anrm;
    snf /= anrm;
    phi         = atan2( snf, csf );
    double phid = fmod( phi - phi0 + M_PI8, M_PI2 );
    if ( phid > M_PI ) phid = phid - M_PI2;
    double drp = ( xp + dr * csf0 + m_rad * ( csf0 - csf ) - xnp ) * csf +
                 ( yp + dr * snf0 + m_rad * ( snf0 - snf ) - ynp ) * snf;
    double dzp = zp + dz - m_rad * tanl * phid - znp;
 
    HepVector ap( 5 );
    ap[0] = drp;
    ap[1] = fmod( phi + M_PI4, M_PI2 );
    ap[2] = kappa;
    ap[3] = dzp;
    ap[4] = tanl;
 
    // std::cout<<" numf_: "<<numf_<<std::endl;
 
    // Modification due to non uniform magnetic field :
    if ( numf_ > 10 )
    {
 
      Hep3Vector x1( xp + dr * csf0, yp + dr * snf0, zp + dz );
 
      double csf0p = cos( ap[1] );
      double snf0p = ( 1. - csf0p ) * ( 1. + csf0p );
      snf0p        = sqrt( ( snf0p > 0. ) ? snf0p : 0. );
      if ( ap[1] > M_PI ) snf0p = -snf0p;
 
      Hep3Vector x2( xnp + drp * csf0p, ynp + drp * snf0p, znp + dzp );
 
      Hep3Vector dist( ( x1 - x2 ).x() / 100.0, ( x1 - x2 ).y() / 100.0,
                       ( x1 - x2 ).z() / 100.0 );
 
      HepPoint3D middlebis( ( x1.x() + x2.x() ) / 2, ( x1.y() + x2.y() ) / 2,
                            ( x1.z() + x2.z() ) / 2 );
 
      HepVector3D field;
      MFSvc_->fieldVector( 10. * middlebis, field );
      field = 1000. * field;
 
      // std::cout<<"B field: "<<field<<std::endl;
      Hep3Vector dB( field.x(), field.y(), ( field.z() - 0.1 * Bznom_ ) );
 
      // std::cout<<" dB: "<<dB<<std::endl;
 
      if ( field.z() )
      {
        double    akappa( fabs( kappa ) );
        double    sign = kappa / akappa;
        HepVector dp( 3 );
        dp = 0.299792458 * sign * dB.cross( dist );
 
        // std::cout<<"dp: "<<dp<<std::endl;
 
        HepVector dhp( 3 );
        dhp[0] = -akappa * ( dp[0] * csf0p + dp[1] * snf0p );
        dhp[1] = kappa * akappa * ( dp[0] * snf0p - dp[1] * csf0p );
        dhp[2] = dp[2] * akappa + dhp[1] * tanl / kappa;
 
        // std::cout<<"dhp: "<<dhp<<std::endl;
 
        ap[1] += dhp[0];
        ap[2] += dhp[1];
        ap[4] += dhp[2];
      }
    }
    HepMatrix m_del = delApDelA( ap );
    Ea_now.assign( m_del * Ea_now * m_del.T() );
    pivot( nextPivot );
    a( ap );
    Ea( Ea_now );
    i++;
 
    // std::cout<<" a: "<<a()<<std::endl;
  }
 
  // Estimation of the path length :
  double tanl_0( tracktest.a()[4] );
  double phi0_0( tracktest.a()[1] );
  double radius_0( tracktest.radius() );
  tracktest.pivot( newPivot );
 
  double phi0_1 = tracktest.a()[1];
  if ( fabs( phi0_1 - phi0_0 ) > M_PI )
  {
    if ( phi0_1 > phi0_0 ) phi0_1 -= 2 * M_PI;
    else phi0_0 -= 2 * M_PI;
  }
  if ( phi0_1 == phi0_0 ) phi0_1 = phi0_0 + 1.e-10;
  pathl = fabs( radius_0 * ( phi0_1 - phi0_0 ) * sqrt( 1 + tanl_0 * tanl_0 ) );
  return newPivot;
}

pivot_numf() [5/6]

const HepPoint3D & KalFitTrack::pivot_numf	(	const HepPoint3D &	newPivot,
		double &	pathl )

pivot_numf() [6/6]

const HepPoint3D & KalFitTrack::pivot_numf	(	const HepPoint3D &	newPivot,
		double &	pathl )

point_last() [1/6]

void KalFitTrack::point_last ( const HepPoint3D & point )

inline

set and give out the last point of the track

Definition at line 133 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ point_last_ = point; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::smoother_anal().

point_last() [2/6]

void KalFitTrack::point_last ( const HepPoint3D & point )

inline

set and give out the last point of the track

Definition at line 133 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ point_last_ = point; }

point_last() [3/6]

void KalFitTrack::point_last ( const HepPoint3D & point )

inline

set and give out the last point of the track

Definition at line 133 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ point_last_ = point; }

point_last() [4/6]

const HepPoint3D & KalFitTrack::point_last ( void )

inline

Definition at line 134 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return point_last_; }

point_last() [5/6]

const HepPoint3D & KalFitTrack::point_last ( void )

inline

Definition at line 134 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return point_last_; }

point_last() [6/6]

const HepPoint3D & KalFitTrack::point_last ( void )

inline

Definition at line 134 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return point_last_; }

r0() [1/3]

double KalFitTrack::r0 ( void ) const

inline

Definition at line 171 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return r0_; }

r0() [2/3]

double KalFitTrack::r0 ( void ) const

inline

Definition at line 171 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return r0_; }

r0() [3/3]

double KalFitTrack::r0 ( void ) const

inline

Definition at line 171 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return r0_; }

r_max() [1/3]

double KalFitTrack::r_max ( void ) const

inline

Definition at line 187 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return r_max_; }

r_max() [2/3]

double KalFitTrack::r_max ( void ) const

inline

Definition at line 187 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return r_max_; }

r_max() [3/3]

double KalFitTrack::r_max ( void ) const

inline

Definition at line 187 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return r_max_; }

radius_numf() [1/3]

double KalFitTrack::radius_numf

(

void

)

const

Estimation of the radius in a given mag field.

Definition at line 1393 of file KalFitTrack.cxx.

                                            {
 
  double Bz( Bznom_ );
  // std::cout<<"Bz: "<<Bz<<std::endl;
  if ( numf_ > 10 )
  {
    double dr   = a()[0];
    double phi0 = a()[1];
    double dz   = a()[3];
    double phi  = fmod( phi0 + M_PI4, M_PI2 );
    double csf0 = cos( phi );
    double snf0 = ( 1. - csf0 ) * ( 1. + csf0 );
    snf0        = sqrt( ( snf0 > 0. ) ? snf0 : 0. );
    if ( phi > M_PI ) snf0 = -snf0;
    // XYZPoint
    HepPoint3D x0( ( pivot().x() + dr * csf0 ), ( pivot().y() + dr * snf0 ),
                   ( pivot().z() + dz ) );
 
    // XYZVector b;
    HepVector3D b;
 
    // std::cout<<"b: "<<b<<std::endl;
 
    MFSvc_->fieldVector( 10. * x0, b );
 
    // std::cout<<"b: "<<b<<std::endl;
 
    b  = 10000. * b;
    Bz = b.z();
  }
  if ( Bz == 0 ) Bz = Bznom_;
  double ALPHA_loc = 10000. / 2.99792458 / Bz;
  return ALPHA_loc / a()[2];
}

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and pivot_numf().

radius_numf() [2/3]

double KalFitTrack::radius_numf

(

void

)

const

Estimation of the radius in a given mag field.

radius_numf() [3/3]

double KalFitTrack::radius_numf

(

void

)

const

Estimation of the radius in a given mag field.

resetLayerUsed() [1/3]

void KalFitTrack::resetLayerUsed ( )

inline

Definition at line 332 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

                        {
    for ( int i = 0; i < 43; i++ ) myLayerUsed[i] = 0;
  }

resetLayerUsed() [2/3]

void KalFitTrack::resetLayerUsed ( )

inline

Definition at line 332 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

                        {
    for ( int i = 0; i < 43; i++ ) myLayerUsed[i] = 0;
  }

resetLayerUsed() [3/3]

void KalFitTrack::resetLayerUsed ( )

inline

Definition at line 332 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

                        {
    for ( int i = 0; i < 43; i++ ) myLayerUsed[i] = 0;
  }

resol() [1/6]

void KalFitTrack::resol ( int i )

static

Definition at line 1868 of file KalFitTrack.cxx.

{ resolflag_ = i; }

resol() [2/6]

void KalFitTrack::resol ( int i )

static

resol() [3/6]

void KalFitTrack::resol ( int i )

static

resol() [4/6]

int KalFitTrack::resol ( void )

static

Definition at line 1869 of file KalFitTrack.cxx.

{ return resolflag_; }

Referenced by KalFitAlg::initialize().

resol() [5/6]

int KalFitTrack::resol ( void )

static

resol() [6/6]

int KalFitTrack::resol ( void )

static

setIMdcGeomSvc() [1/3]

void KalFitTrack::setIMdcGeomSvc ( IMdcGeomSvc * igeomsvc )

static

Definition at line 228 of file KalFitTrack2.cxx.

{ iGeomSvc_ = igeomsvc; }

Referenced by KalFitAlg::setGeomSvc_init().

setIMdcGeomSvc() [2/3]

void KalFitTrack::setIMdcGeomSvc ( IMdcGeomSvc * igeomsvc )

static

setIMdcGeomSvc() [3/3]

void KalFitTrack::setIMdcGeomSvc ( IMdcGeomSvc * igeomsvc )

static

setInitMatrix() [1/3]

void KalFitTrack::setInitMatrix ( HepSymMatrix m )

static

Definition at line 35 of file KalFitTrack2.cxx.

{ initMatrix_ = m; }

Referenced by KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), and KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew().

setInitMatrix() [2/3]

void KalFitTrack::setInitMatrix ( HepSymMatrix m )

static

setInitMatrix() [3/3]

void KalFitTrack::setInitMatrix ( HepSymMatrix m )

static

setMagneticFieldSvc() [1/3]

void KalFitTrack::setMagneticFieldSvc ( IBesMagFieldSvc * mf )

static

Definition at line 218 of file KalFitTrack2.cxx.

                                                           {
  /*ISvcLocator* svcLocator = Gaudi::svcLocator();
    StatusCode sc = svcLocator->service("MagneticFieldSvc",MFSvc_);
    if (sc.isFailure()){
    std::cout << "Could not load MagneticFieldSvc!" << std::endl;
    }*/
  MFSvc_ = mf;
  if ( MFSvc_ == 0 ) cout << "KalFitTrack2:: Could not load MagneticFieldSvc!" << endl;
}

Referenced by KalFitAlg::initialize().

setMagneticFieldSvc() [2/3]

void KalFitTrack::setMagneticFieldSvc ( IBesMagFieldSvc * )

static

setMagneticFieldSvc() [3/3]

void KalFitTrack::setMagneticFieldSvc ( IBesMagFieldSvc * )

static

setMdcCalibFunSvc() [1/3]

void KalFitTrack::setMdcCalibFunSvc ( const IMdcCalibFunSvc * calibsvc )

static

Definition at line 214 of file KalFitTrack2.cxx.

                                                                     {
  CalibFunSvc_ = calibsvc;
}

Referenced by KalFitAlg::setCalibSvc_init().

setMdcCalibFunSvc() [2/3]

void KalFitTrack::setMdcCalibFunSvc ( const IMdcCalibFunSvc * calibsvc )

static

setMdcCalibFunSvc() [3/3]

void KalFitTrack::setMdcCalibFunSvc ( const IMdcCalibFunSvc * calibsvc )

static

setMdcDigiCol() [1/3]

void KalFitTrack::setMdcDigiCol ( MdcDigiCol * digicol )

static

Definition at line 230 of file KalFitTrack2.cxx.

{ mdcDigiCol_ = digicol; }

Referenced by KalFitAlg::execute().

setMdcDigiCol() [2/3]

void KalFitTrack::setMdcDigiCol ( MdcDigiCol * digicol )

static

setMdcDigiCol() [3/3]

void KalFitTrack::setMdcDigiCol ( MdcDigiCol * digicol )

static

setT0() [1/3]

void KalFitTrack::setT0 ( double t0 )

static

Definition at line 39 of file KalFitTrack2.cxx.

                                           {
  //------------------set  event  start time-----------
 
  EventT0_ = eventstart;
  if ( debug_ == 4 )
  { std::cout << "in function KalFitTrack::setT0(...), EventT0_ = " << EventT0_ << std::endl; }
}

Referenced by KalFitAlg::execute().

setT0() [2/3]

void KalFitTrack::setT0 ( double t0 )

static

setT0() [3/3]

void KalFitTrack::setT0 ( double t0 )

static

smoother_Mdc() [1/6]

double KalFitTrack::smoother_Mdc	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

Kalman smoother for Mdc.

Referenced by KalFitAlg::smoother_anal(), and KalFitAlg::smoother_calib().

smoother_Mdc() [2/6]

double KalFitTrack::smoother_Mdc	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

Kalman smoother for Mdc.

smoother_Mdc() [3/6]

double KalFitTrack::smoother_Mdc	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

Kalman smoother for Mdc.

smoother_Mdc() [4/6]

double KalFitTrack::smoother_Mdc	(	KalFitHitMdc &	HitMdc,
		CLHEP::Hep3Vector &	meas,
		KalFitHelixSeg &	seg,
		double &	dchi2,
		int	csmflag )

smoother_Mdc() [5/6]

double KalFitTrack::smoother_Mdc	(	KalFitHitMdc &	HitMdc,
		CLHEP::Hep3Vector &	meas,
		KalFitHelixSeg &	seg,
		double &	dchi2,
		int	csmflag )

smoother_Mdc() [6/6]

double KalFitTrack::smoother_Mdc	(	KalFitHitMdc &	HitMdc,
		CLHEP::Hep3Vector &	meas,
		KalFitHelixSeg &	seg,
		double &	dchi2,
		int	csmflag )

smoother_Mdc_csmalign() [1/3]

double KalFitTrack::smoother_Mdc_csmalign	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

Referenced by KalFitAlg::smoother_calib().

smoother_Mdc_csmalign() [2/3]

double KalFitTrack::smoother_Mdc_csmalign	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

smoother_Mdc_csmalign() [3/3]

double KalFitTrack::smoother_Mdc_csmalign	(	KalFitHelixSeg &	seg,
		CLHEP::Hep3Vector &	meas,
		int &	flg,
		int	csmflag )

tof() [1/6]

void KalFitTrack::tof

(

double

path

)

Update the tof estimation.

Definition at line 1365 of file KalFitTrack.cxx.

                                   {
  double light_speed( 29.9792458 ); // light speed in cm/nsec
  double t = tanl();
  double pmag( sqrt( 1.0 + t * t ) / kappa() );
  if ( pmag != 0 )
  {
    double mass_over_p( mass_ / pmag );
    double beta( 1.0 / sqrt( 1.0 + mass_over_p * mass_over_p ) );
    tof_ += path / ( light_speed * beta );
  }
 
  if ( tofall_ )
  {
    if ( p_kaon_ > 0 )
    {
      double massk_over_p( MASS[3] / p_kaon_ );
      double beta_kaon( 1.0 / sqrt( 1.0 + massk_over_p * massk_over_p ) );
      tof_kaon_ += path / ( light_speed * beta_kaon );
    }
    if ( p_proton_ > 0 )
    {
      double massp_over_p( MASS[4] / p_proton_ );
      double beta_proton( 1.0 / sqrt( 1.0 + massp_over_p * massp_over_p ) );
      tof_proton_ += path / ( light_speed * beta_proton );
    }
  }
}

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

tof() [2/6]

void KalFitTrack::tof

(

double

path

)

Update the tof estimation.

tof() [3/6]

void KalFitTrack::tof

(

double

path

)

Update the tof estimation.

tof() [4/6]

double KalFitTrack::tof ( void ) const

inline

Definition at line 181 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return tof_; }

Referenced by msgasmdc(), and path_add().

tof() [5/6]

double KalFitTrack::tof ( void ) const

inline

Definition at line 181 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return tof_; }

tof() [6/6]

double KalFitTrack::tof ( void ) const

inline

Definition at line 181 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return tof_; }

tof_kaon() [1/3]

double KalFitTrack::tof_kaon ( void ) const

inline

Definition at line 182 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return tof_kaon_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

tof_kaon() [2/3]

double KalFitTrack::tof_kaon ( void ) const

inline

Definition at line 182 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return tof_kaon_; }

tof_kaon() [3/3]

double KalFitTrack::tof_kaon ( void ) const

inline

Definition at line 182 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return tof_kaon_; }

tof_proton() [1/3]

double KalFitTrack::tof_proton ( void ) const

inline

Definition at line 183 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return tof_proton_; }

Referenced by KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), and KalFitAlg::fillTds_back().

tof_proton() [2/3]

double KalFitTrack::tof_proton ( void ) const

inline

Definition at line 183 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return tof_proton_; }

tof_proton() [3/3]

double KalFitTrack::tof_proton ( void ) const

inline

Definition at line 183 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return tof_proton_; }

trasan_id() [1/6]

void KalFitTrack::trasan_id ( int t )

inline

Definition at line 198 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ trasan_id_ = t; }

trasan_id() [2/6]

void KalFitTrack::trasan_id ( int t )

inline

Definition at line 198 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ trasan_id_ = t; }

trasan_id() [3/6]

void KalFitTrack::trasan_id ( int t )

inline

Definition at line 198 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ trasan_id_ = t; }

trasan_id() [4/6]

int KalFitTrack::trasan_id ( void ) const

inline

Definition at line 170 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return trasan_id_; }

trasan_id() [5/6]

int KalFitTrack::trasan_id ( void ) const

inline

Definition at line 170 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return trasan_id_; }

trasan_id() [6/6]

int KalFitTrack::trasan_id ( void ) const

inline

Definition at line 170 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return trasan_id_; }

type() [1/6]

void KalFitTrack::type ( int t )

inline

Reinitialize (modificator).

Definition at line 197 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ type_ = t; }

type() [2/6]

void KalFitTrack::type ( int t )

inline

Reinitialize (modificator).

Definition at line 197 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ type_ = t; }

type() [3/6]

void KalFitTrack::type ( int t )

inline

Reinitialize (modificator).

Definition at line 197 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ type_ = t; }

type() [4/6]

int KalFitTrack::type ( void ) const

inline

Definition at line 169 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

{ return type_; }

Referenced by KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), and KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew().

type() [5/6]

int KalFitTrack::type ( void ) const

inline

Definition at line 169 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

{ return type_; }

type() [6/6]

int KalFitTrack::type ( void ) const

inline

Definition at line 169 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

{ return type_; }

update_bit() [1/3]

void KalFitTrack::update_bit

(

int

i

)

Definition at line 1843 of file KalFitTrack.cxx.

                                    {
  int j( 0 );
  if ( i < 31 )
  {
    j = (int)pow( 2., i );
    if ( !( pat1_ & j ) ) pat1_ = pat1_ | j;
  }
  else if ( i < 50 )
  {
    j = (int)pow( 2., ( i - 31 ) );
    if ( !( pat2_ & j ) ) pat2_ = pat2_ | j;
  }
}

update_bit() [2/3]

void KalFitTrack::update_bit

(

int

i

)

update_bit() [3/3]

void KalFitTrack::update_bit

(

int

i

)

update_forMdc() [1/3]

void KalFitTrack::update_forMdc

(

void

)

Definition at line 121 of file KalFitTrack.cxx.

                                      {
  pivot_forMdc_ = pivot();
  a_forMdc_     = a();
  Ea_forMdc_    = Ea();
}

Referenced by KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), and KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew().

update_forMdc() [2/3]

void KalFitTrack::update_forMdc

(

void

)

update_forMdc() [3/3]

void KalFitTrack::update_forMdc

(

void

)

update_hits() [1/6]

double KalFitTrack::update_hits	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

update_hits() [2/6]

double KalFitTrack::update_hits	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

update_hits() [3/6]

double KalFitTrack::update_hits	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

update_hits() [4/6]

double KalFitTrack::update_hits	(	KalFitHitMdc &	HitMdc,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		double &	dtrack,
		double &	dtracknew,
		double &	dtdc,
		int	csmflag )

Include the Mdc wire hits.

Referenced by KalFitAlg::filter_fwd_anal(), and KalFitAlg::filter_fwd_calib().

update_hits() [5/6]

double KalFitTrack::update_hits	(	KalFitHitMdc &	HitMdc,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		double &	dtrack,
		double &	dtracknew,
		double &	dtdc,
		int	csmflag )

Include the Mdc wire hits.

update_hits() [6/6]

double KalFitTrack::update_hits	(	KalFitHitMdc &	HitMdc,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		double &	dtrack,
		double &	dtracknew,
		double &	dtdc,
		int	csmflag )

Include the Mdc wire hits.

update_hits_csmalign() [1/3]

double KalFitTrack::update_hits_csmalign	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

Referenced by KalFitAlg::filter_fwd_calib().

update_hits_csmalign() [2/3]

double KalFitTrack::update_hits_csmalign	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

update_hits_csmalign() [3/3]

double KalFitTrack::update_hits_csmalign	(	KalFitHelixSeg &	HelixSeg,
		int	inext,
		CLHEP::Hep3Vector &	meas,
		int	way,
		double &	dchi2,
		int	csmflag )

update_last() [1/3]

void KalFitTrack::update_last

(

void

)

Record the current parameters as ..._last information :

Definition at line 115 of file KalFitTrack.cxx.

                                    {
  pivot_last_ = pivot();
  a_last_     = a();
  Ea_last_    = Ea();
}

update_last() [2/3]

void KalFitTrack::update_last

(

void

)

Record the current parameters as ..._last information :

update_last() [3/3]

void KalFitTrack::update_last

(

void

)

Record the current parameters as ..._last information :

useLayer() [1/3]

void KalFitTrack::useLayer ( int iLay )

inline

Definition at line 336 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

                            {
    if ( iLay >= 0 && iLay <= 43 ) myLayerUsed[iLay] = 1;
  }

useLayer() [2/3]

void KalFitTrack::useLayer ( int iLay )

inline

Definition at line 336 of file InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h.

                            {
    if ( iLay >= 0 && iLay <= 43 ) myLayerUsed[iLay] = 1;
  }

useLayer() [3/3]

void KalFitTrack::useLayer ( int iLay )

inline

Definition at line 336 of file Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h.

                            {
    if ( iLay >= 0 && iLay <= 43 ) myLayerUsed[iLay] = 1;
  }

Member Data Documentation

Bznom_

double KalFitTrack::Bznom_ = 10

static

Definition at line 289 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::beginRun(), KalFitAlg::execute(), KalFitAlg::initialize(), KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew(), pivot_numf(), pivot_numf(), radius_numf(), and setMdcDigiCol().

chi2_hitf_

double KalFitTrack::chi2_hitf_ = 1000

static

Cut chi2 for each hit.

Definition at line 272 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), and setMdcDigiCol().

chi2_hits_

double KalFitTrack::chi2_hits_ = 1000

static

Definition at line 272 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), and setMdcDigiCol().

chi2mdc_hit2_

double KalFitTrack::chi2mdc_hit2_

static

Definition at line 303 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by setMdcDigiCol().

dchi2cutf_anal

double KalFitTrack::dchi2cutf_anal = { 0. }

static

Definition at line 56 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut(), and setMdcDigiCol().

dchi2cutf_calib

double KalFitTrack::dchi2cutf_calib = { 0. }

static

Definition at line 58 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut(), and setMdcDigiCol().

dchi2cuts_anal

double KalFitTrack::dchi2cuts_anal = { 0. }

static

Definition at line 57 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut(), and setMdcDigiCol().

dchi2cuts_calib

double KalFitTrack::dchi2cuts_calib = { 0. }

static

Definition at line 59 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut(), and setMdcDigiCol().

debug_

int KalFitTrack::debug_ = 0

static

for debug

Definition at line 270 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by getDriftDist(), getDriftTime(), getT0(), KalFitAlg::initialize(), setMdcDigiCol(), and setT0().

drifttime_choice_

int KalFitTrack::drifttime_choice_ = 0

static

the drifttime choice

Definition at line 311 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by getDriftTime(), KalFitAlg::initialize(), and setMdcDigiCol().

factor_strag_

double KalFitTrack::factor_strag_ = 0.4

static

factor of energy loss straggling for electron

Definition at line 300 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by eloss(), KalFitAlg::initialize(), and setMdcDigiCol().

inner_steps_

int KalFitTrack::inner_steps_ = 3

static

Definition at line 275 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), pivot_numf(), and setMdcDigiCol().

LR_

int KalFitTrack::LR_ = 1

static

Use L/R decision from MdcRecHit information :

Definition at line 309 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by chi2_next(), chi2_next(), KalFitAlg::initialize(), KalFitAlg::kalman_fitting_anal(), KalFitAlg::kalman_fitting_calib(), KalFitAlg::kalman_fitting_csmalign(), KalFitAlg::kalman_fitting_MdcxReco_Csmc_Sew(), LR(), and setMdcDigiCol().

mdcGasRadlen_

double KalFitTrack::mdcGasRadlen_ = 0.

static

Definition at line 264 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by msgasmdc(), KalFitAlg::setBesFromGdml(), and setMdcDigiCol().

nmdc_hit2_

int KalFitTrack::nmdc_hit2_ = 500

static

Cut chi2 for each hit.

Definition at line 302 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), and setMdcDigiCol().

numf_

int KalFitTrack::numf_ = 0

static

Flag for treatment of non-uniform mag field.

Definition at line 274 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), numf(), numf(), pivot_numf(), pivot_numf(), radius_numf(), and setMdcDigiCol().

numfcor_

int KalFitTrack::numfcor_ = 1

static

NUMF treatment improved.

Definition at line 287 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::beginRun(), KalFitAlg::execute(), KalFitAlg::initialize(), pivot_numf(), and setMdcDigiCol().

outer_steps_

int KalFitTrack::outer_steps_ = 3

static

Definition at line 276 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), pivot_numf(), and setMdcDigiCol().

resolflag_

int KalFitTrack::resolflag_ = 0

static

wire resoltion flag

Definition at line 306 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by chi2_next(), chi2_next(), resol(), resol(), and setMdcDigiCol().

steplev_

int KalFitTrack::steplev_ = 0

static

Level of precision (= 1 : 5steps for all tracks; = 2: 5 steps for very non uniform part)

Definition at line 292 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize(), and setMdcDigiCol().

strag_

int KalFitTrack::strag_ = 1

static

Flag to take account of energy loss straggling :

Definition at line 298 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by eloss(), KalFitAlg::initialize(), and setMdcDigiCol().

Tof_correc_

int KalFitTrack::Tof_correc_ = 1

static

Flag for TOF correction.

Definition at line 295 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by chi2_next(), chi2_next(), KalFitAlg::initialize(), and setMdcDigiCol().

tofall_

int KalFitTrack::tofall_ = 1

static

Definition at line 304 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::complete_track(), KalFitAlg::complete_track(), eloss(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::fillTds_back(), KalFitAlg::initialize(), setMdcDigiCol(), and tof().

tprop_

int KalFitTrack::tprop_ = 1

static

for signal propagation correction

Definition at line 267 of file InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h.

Referenced by getDriftTime(), KalFitAlg::initialize(), and setMdcDigiCol().

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/KalFitAlg/KalFitTrack.h
• InstallArea/x86_64-el9-gcc13-opt/include/KalFitAlg/KalFitTrack.h
• Reconstruction/KalFitAlg/include/KalFitAlg/KalFitTrack.h
• KalFitTrack.cxx
• KalFitTrack2.cxx