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	update_hits (RecCgemCluster *Cluster, double recR, int way, 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	setClusterRefVec (ClusterRefVec clusterRefVec)

ClusterRefVec	getClusterRefVec ()

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	update_hits (RecCgemCluster *Cluster, double recR, int way, 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	setClusterRefVec (ClusterRefVec clusterRefVec)

ClusterRefVec	getClusterRefVec ()

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)
	unsets error matrix. Error calculations will be ignored after this function call until an error matrix be set again. 0 matrix will be return as a return value for error matrix when you call functions which returns an error matrix.

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)
	unsets error matrix. Error calculations will be ignored after this function call until an error matrix be set again. 0 matrix will be return as a return value for error matrix when you call functions which returns an error matrix.

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 MdcCalibFunSvc *calibsvc)

static void	setMagneticFieldSvc (IMagneticFieldSvc *)

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 MdcCalibFunSvc *calibsvc)

static void	setMagneticFieldSvc (IMagneticFieldSvc *)

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 35 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Constructor & Destructor Documentation

◆ KalFitTrack() [1/2]

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

constructor

◆ ~KalFitTrack() [1/2]

KalFitTrack::~KalFitTrack ( void )

destructor

Definition at line 108 of file KalFitTrack.cxx.

{
    // delete all objects
 
}

◆ KalFitTrack() [2/2]

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

constructor

◆ ~KalFitTrack() [2/2]

KalFitTrack::~KalFitTrack ( void )

destructor

Member Function Documentation

◆ a_forMdc() [1/2]

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

inline

Definition at line 156 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

156{ return a_forMdc_; }

◆ a_forMdc() [2/2]

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

inline

Definition at line 156 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

156{ return a_forMdc_; }

◆ a_last() [1/2]

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

inline

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

152{ return a_last_; }

◆ a_last() [2/2]

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

inline

Definition at line 152 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

152{ return a_last_; }

◆ add_nhit_r() [1/2]

void KalFitTrack::add_nhit_r ( void )

inline

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

225{ nhit_r_++;}

◆ add_nhit_r() [2/2]

void KalFitTrack::add_nhit_r ( void )

inline

Definition at line 225 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

225{ nhit_r_++;}

◆ add_nhit_z() [1/2]

void KalFitTrack::add_nhit_z ( void )

inline

Definition at line 226 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

226{ nhit_z_++;}

◆ add_nhit_z() [2/2]

void KalFitTrack::add_nhit_z ( void )

inline

Definition at line 226 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

226{ nhit_z_++;}

◆ addPathSM() [1/2]

void KalFitTrack::addPathSM ( double path )

Definition at line 1301 of file KalFitTrack.cxx.

                                      {
    pathSM_ += path;
}

Referenced by KalFitAlg::smoother_anal().

◆ addPathSM() [2/2]

void KalFitTrack::addPathSM ( double path )

◆ addTofSM() [1/2]

void KalFitTrack::addTofSM ( double time )

Definition at line 1306 of file KalFitTrack.cxx.

                                     {
    tof2_ += time;
} 

Referenced by KalFitAlg::smoother_anal().

◆ addTofSM() [2/2]

void KalFitTrack::addTofSM ( double time )

◆ appendHelixSegs() [1/2]

void KalFitTrack::appendHelixSegs ( KalFitHelixSeg s )

◆ appendHelixSegs() [2/2]

void KalFitTrack::appendHelixSegs ( KalFitHelixSeg s )

◆ appendHitsMdc() [1/2]

void KalFitTrack::appendHitsMdc ( KalFitHitMdc h )

Functions for Mdc hits list.

Definition at line 1842 of file KalFitTrack.cxx.

1842{ 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/2]

void KalFitTrack::appendHitsMdc ( KalFitHitMdc h )

Functions for Mdc hits list.

◆ back() [1/4]

void KalFitTrack::back ( int i )

static

Definition at line 1834 of file KalFitTrack.cxx.

1834{ back_ = i;}

◆ back() [2/4]

static void KalFitTrack::back ( int i )

static

◆ back() [3/4]

int KalFitTrack::back ( void )

static

Definition at line 1835 of file KalFitTrack.cxx.

1835{ return back_;}

Referenced by KalFitAlg::initialize().

◆ back() [4/4]

static int KalFitTrack::back ( void )

static

◆ chgmass() [1/2]

void KalFitTrack::chgmass ( int i )

Definition at line 1827 of file KalFitTrack.cxx.

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

Referenced by KalFitAlg::complete_track().

◆ chgmass() [2/2]

void KalFitTrack::chgmass ( int i )

◆ chi2_next() [1/4]

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

Definition at line 2931 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/4]

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

◆ chi2_next() [3/4]

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

Definition at line 3076 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() [4/4]

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

◆ chiSq() [1/4]

void KalFitTrack::chiSq ( double c )

inline

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

219{ chiSq_ = c; }

◆ chiSq() [2/4]

void KalFitTrack::chiSq ( double c )

inline

Definition at line 219 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

219{ chiSq_ = c; }

◆ chiSq() [3/4]

double KalFitTrack::chiSq ( void ) const

inline

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

188{ return chiSq_; }

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

◆ chiSq() [4/4]

double KalFitTrack::chiSq ( void ) const

inline

Definition at line 188 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

188{ return chiSq_; }

◆ chiSq_back() [1/4]

void KalFitTrack::chiSq_back ( double c )

inline

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

220{ chiSq_back_ = c; }

◆ chiSq_back() [2/4]

void KalFitTrack::chiSq_back ( double c )

inline

Definition at line 220 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

220{ chiSq_back_ = c; }

◆ chiSq_back() [3/4]

double KalFitTrack::chiSq_back ( void ) const

inline

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

189{ return chiSq_back_; }

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

◆ chiSq_back() [4/4]

double KalFitTrack::chiSq_back ( void ) const

inline

Definition at line 189 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

189{ return chiSq_back_; }

◆ cor_tanldep() [1/2]

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

Correct the error according the current tanl value :

◆ cor_tanldep() [2/2]

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

Correct the error according the current tanl value :

◆ dchi2_max() [1/2]

double KalFitTrack::dchi2_max ( void ) const

inline

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

201{ return dchi2_max_; }

◆ dchi2_max() [2/2]

double KalFitTrack::dchi2_max ( void ) const

inline

Definition at line 201 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

201{ return dchi2_max_; }

◆ Ea_forMdc() [1/2]

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

inline

Definition at line 157 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

157{ return Ea_forMdc_; }

◆ Ea_forMdc() [2/2]

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

inline

Definition at line 157 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

157{ return Ea_forMdc_; }

◆ Ea_last() [1/2]

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

inline

Definition at line 153 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

153{ return Ea_last_; }

◆ Ea_last() [2/2]

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

inline

Definition at line 153 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

153{ return Ea_last_; }

◆ eloss() [1/2]

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

Calculate total energy lost in material.

Definition at line 281 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::Cgem_filter_anal(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::fitGemHits(), KalFitAlg::fromFHitToInnerWall(), KalFitAlg::smoother_anal(), KalFitAlg::smoother_calib(), KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

◆ eloss() [2/2]

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

Calculate total energy lost in material.

◆ filter() [1/2]

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

◆ filter() [2/2]

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

◆ fiTerm() [1/2]

void KalFitTrack::fiTerm ( double fi )

Definition at line 1311 of file KalFitTrack.cxx.

                                 {
    fiTerm_ = fi;
}

Referenced by KalFitAlg::smoother_anal().

◆ fiTerm() [2/2]

void KalFitTrack::fiTerm ( double fi )

◆ getClusterRefVec() [1/2]

ClusterRefVec KalFitTrack::getClusterRefVec ( )

inline

Definition at line 243 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

243{ return ClusterRefVec_;} //added hy Huang Zhen

◆ getClusterRefVec() [2/2]

ClusterRefVec KalFitTrack::getClusterRefVec ( )

inline

Definition at line 243 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

243{ return ClusterRefVec_;} //added hy Huang Zhen

◆ getDigi() [1/2]

double KalFitTrack::getDigi ( ) const

◆ getDigi() [2/2]

double KalFitTrack::getDigi ( ) const

◆ getDriftDist() [1/2]

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

Definition at line 194 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().

◆ getDriftDist() [2/2]

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

◆ getDriftTime() [1/2]

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

Definition at line 74 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(); 
}

Referenced by chi2_next().

◆ getDriftTime() [2/2]

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

◆ getFiTerm() [1/2]

double KalFitTrack::getFiTerm ( void )

inline

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

169{ return fiTerm_;}

Referenced by KalFitAlg::fillTds_back().

◆ getFiTerm() [2/2]

double KalFitTrack::getFiTerm ( void )

inline

Definition at line 169 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

169{ return fiTerm_;}

◆ getInitMatrix() [1/2]

HepSymMatrix KalFitTrack::getInitMatrix ( void ) const

Definition at line 42 of file KalFitTrack2.cxx.

{
    return initMatrix_ ;
}

Referenced by KalFitAlg::smoother_calib().

◆ getInitMatrix() [2/2]

HepSymMatrix KalFitTrack::getInitMatrix ( void ) const

◆ getPathSM() [1/2]

double KalFitTrack::getPathSM ( void )

inline

Definition at line 163 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

163{ return pathSM_;}

Referenced by KalFitAlg::fillTds_back().

◆ getPathSM() [2/2]

double KalFitTrack::getPathSM ( void )

inline

Definition at line 163 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

163{ return pathSM_;}

◆ getSigma() [1/4]

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

Definition at line 3222 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().

◆ getSigma() [2/4]

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

◆ getSigma() [3/4]

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

Definition at line 209 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() [4/4]

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

◆ getT0() [1/2]

double KalFitTrack::getT0 ( void ) const

Definition at line 57 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/2]

double KalFitTrack::getT0 ( void ) const

◆ getTofSM() [1/2]

double KalFitTrack::getTofSM ( void )

inline

Definition at line 166 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

166{ return tof2_;}

Referenced by KalFitAlg::fillTds_back().

◆ getTofSM() [2/2]

double KalFitTrack::getTofSM ( void )

inline

Definition at line 166 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

166{ return tof2_;}

◆ HelixSeg() [1/2]

KalFitHelixSeg & KalFitTrack::HelixSeg ( int i )

inline

Definition at line 240 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

240{ return HelixSegs_[i];}

Referenced by KalFitAlg::smoother_calib().

◆ HelixSeg() [2/2]

KalFitHelixSeg & KalFitTrack::HelixSeg ( int i )

inline

Definition at line 240 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

240{ return HelixSegs_[i];}

◆ HelixSegs() [1/4]

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

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

◆ HelixSegs() [2/4]

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

◆ HelixSegs() [3/4]

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

inline

Definition at line 239 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

239{ return HelixSegs_;}

◆ HelixSegs() [4/4]

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

inline

Definition at line 239 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

239{ return HelixSegs_;}

◆ HitMdc() [1/2]

KalFitHitMdc & KalFitTrack::HitMdc ( int i )

inline

Definition at line 235 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

235{ return HitsMdc_[i];}

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

◆ HitMdc() [2/2]

KalFitHitMdc & KalFitTrack::HitMdc ( int i )

inline

Definition at line 235 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

235{ return HitsMdc_[i];}

◆ HitsMdc() [1/4]

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

Definition at line 1841 of file KalFitTrack.cxx.

1841{ HitsMdc_ = lh;}

Referenced by 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/4]

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

◆ HitsMdc() [3/4]

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

inline

Definition at line 234 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

234{ return HitsMdc_;}

Referenced by order_hits().

◆ HitsMdc() [4/4]

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

inline

Definition at line 234 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

234{ return HitsMdc_;}

◆ insist() [1/4]

void KalFitTrack::insist ( int t )

inline

Definition at line 214 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

214{ insist_ = t;}

◆ insist() [2/4]

void KalFitTrack::insist ( int t )

inline

Definition at line 214 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

214{ insist_ = t;}

◆ insist() [3/4]

int KalFitTrack::insist ( void ) const

inline

Extractor :

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

183{ return insist_; }

◆ insist() [4/4]

int KalFitTrack::insist ( void ) const

inline

Extractor :

Definition at line 183 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

183{ return insist_; }

◆ intersect_cylinder() [1/2]

double KalFitTrack::intersect_cylinder ( double r ) const

Intersection with different geometry.

Definition at line 128 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 KalFitAlg::Cgem_filter_anal(), KalFitAlg::fitGemHits(), KalFitAlg::fromFHitToInnerWall(), KalFitAlg::innerwall(), KalFitCylinder::intersect(), intersect_yz_plane(), intersect_zx_plane(), and update_hits().

◆ intersect_cylinder() [2/2]

double KalFitTrack::intersect_cylinder ( double r ) const

Intersection with different geometry.

◆ intersect_xy_plane() [1/2]

double KalFitTrack::intersect_xy_plane ( double z ) const

Definition at line 180 of file KalFitTrack.cxx.

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

Referenced by KalFitCylinder::intersect().

◆ intersect_xy_plane() [2/2]

double KalFitTrack::intersect_xy_plane ( double z ) const

◆ intersect_yz_plane() [1/2]

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

Definition at line 162 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/2]

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

◆ intersect_zx_plane() [1/2]

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

Definition at line 144 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/2]

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

◆ lead() [1/4]

void KalFitTrack::lead ( int i )

static

Definition at line 1832 of file KalFitTrack.cxx.

1832{ lead_ = i;}

◆ lead() [2/4]

static void KalFitTrack::lead ( int i )

static

◆ lead() [3/4]

int KalFitTrack::lead ( void )

static

Magnetic field map.

Definition at line 1833 of file KalFitTrack.cxx.

1833{ return lead_;}

Referenced by KalFitAlg::initialize().

◆ lead() [4/4]

static int KalFitTrack::lead ( void )

static

Magnetic field map.

◆ LR() [1/2]

void KalFitTrack::LR ( int x )

static

Definition at line 1840 of file KalFitTrack.cxx.

1840{ LR_ = x;}

Referenced by KalFitAlg::initialize().

◆ LR() [2/2]

static void KalFitTrack::LR ( int x )

static

◆ mass() [1/4]

double KalFitTrack::mass ( int i )

static

Definition at line 1826 of file KalFitTrack.cxx.

1826{ return MASS[i];}

◆ mass() [2/4]

static double KalFitTrack::mass ( int i )

static

◆ mass() [3/4]

double KalFitTrack::mass ( void ) const

inline

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

187{ return mass_; }

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

◆ mass() [4/4]

double KalFitTrack::mass ( void ) const

inline

Definition at line 187 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

187{ return mass_; }

◆ mom() [1/2]

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

inline

Definition at line 195 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

195{ return mom_; }

◆ mom() [2/2]

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

inline

Definition at line 195 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

195{ return mom_; }

◆ ms() [1/2]

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

Definition at line 187 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 KalFitAlg::fitGemHits(), KalFitAlg::fromFHitToInnerWall(), KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

◆ ms() [2/2]

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

◆ msgasmdc() [1/2]

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

Calculate multiple scattering angle.

Definition at line 224 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::Cgem_filter_anal(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::smoother_anal(), and KalFitAlg::smoother_calib().

◆ msgasmdc() [2/2]

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

Calculate multiple scattering angle.

◆ ncath() [1/2]

unsigned int KalFitTrack::ncath ( void ) const

inline

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

205{ return ncath_; }

◆ ncath() [2/2]

unsigned int KalFitTrack::ncath ( void ) const

inline

Definition at line 205 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

205{ return ncath_; }

◆ nchits() [1/4]

void KalFitTrack::nchits ( int n )

inline

Definition at line 222 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

222{ nchits_ = n; }

n

const Int_t n

Definition: DataBase/tau_mode.c:65

◆ nchits() [2/4]

void KalFitTrack::nchits ( int n )

inline

Definition at line 222 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

222{ nchits_ = n; }

◆ nchits() [3/4]

unsigned int KalFitTrack::nchits ( void ) const

inline

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

203{ return nchits_; }

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

◆ nchits() [4/4]

unsigned int KalFitTrack::nchits ( void ) const

inline

Definition at line 203 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

203{ return nchits_; }

◆ ndf_back() [1/4]

void KalFitTrack::ndf_back ( int n )

inline

Definition at line 221 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

221{ ndf_back_ = n; }

◆ ndf_back() [2/4]

void KalFitTrack::ndf_back ( int n )

inline

Definition at line 221 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

221{ ndf_back_ = n; }

◆ ndf_back() [3/4]

int KalFitTrack::ndf_back ( void ) const

inline

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

190{ return ndf_back_; }

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

◆ ndf_back() [4/4]

int KalFitTrack::ndf_back ( void ) const

inline

Definition at line 190 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

190{ return ndf_back_; }

◆ nhit_r() [1/2]

int KalFitTrack::nhit_r ( void ) const

inline

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

208{ return nhit_r_; }

◆ nhit_r() [2/2]

int KalFitTrack::nhit_r ( void ) const

inline

Definition at line 208 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

208{ return nhit_r_; }

◆ nhit_z() [1/2]

int KalFitTrack::nhit_z ( void ) const

inline

Definition at line 209 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

209{ return nhit_z_; }

◆ nhit_z() [2/2]

int KalFitTrack::nhit_z ( void ) const

inline

Definition at line 209 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

209{ return nhit_z_; }

◆ nLayerUsed() [1/2]

int KalFitTrack::nLayerUsed ( )

inline

Definition at line 343 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ nLayerUsed() [2/2]

int KalFitTrack::nLayerUsed ( )

inline

Definition at line 343 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

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

◆ nmass() [1/2]

int KalFitTrack::nmass ( void )

static

Definition at line 1825 of file KalFitTrack.cxx.

1825{ return NMASS;}

Referenced by KalFitAlg::complete_track().

◆ nmass() [2/2]

static int KalFitTrack::nmass ( void )

static

◆ nster() [1/4]

void KalFitTrack::nster ( int n )

inline

Definition at line 223 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

223{ nster_ = n; }

◆ nster() [2/4]

void KalFitTrack::nster ( int n )

inline

Definition at line 223 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

223{ nster_ = n; }

◆ nster() [3/4]

unsigned int KalFitTrack::nster ( void ) const

inline

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

204{ return nster_; }

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

◆ nster() [4/4]

unsigned int KalFitTrack::nster ( void ) const

inline

Definition at line 204 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

204{ return nster_; }

◆ number_wirhit() [1/2]

void KalFitTrack::number_wirhit ( void )

Definition at line 457 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/2]

void KalFitTrack::number_wirhit ( void )

◆ numf() [1/4]

void KalFitTrack::numf ( int i )

static

Definition at line 1838 of file KalFitTrack.cxx.

1838{ numf_ = i;}

KalFitTrack::numf_

static int numf_

Flag for treatment of non-uniform mag field.

Definition: InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h:291

◆ numf() [2/4]

static void KalFitTrack::numf ( int i )

static

◆ numf() [3/4]

int KalFitTrack::numf ( void )

static

Definition at line 1839 of file KalFitTrack.cxx.

1839{ return numf_;}

Referenced by KalFitAlg::initialize().

◆ numf() [4/4]

static int KalFitTrack::numf ( void )

static

◆ order_hits() [1/2]

void KalFitTrack::order_hits ( void )

Definition at line 442 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/2]

void KalFitTrack::order_hits ( void )

◆ order_wirhit() [1/2]

void KalFitTrack::order_wirhit ( int index )

Modifier Order the wire hits for mdc track

Definition at line 380 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/2]

void KalFitTrack::order_wirhit ( int index )

Modifier Order the wire hits for mdc track

◆ p_kaon() [1/4]

void KalFitTrack::p_kaon ( double pl )

inline

Definition at line 217 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

217{ p_kaon_ = pl;}

◆ p_kaon() [2/4]

void KalFitTrack::p_kaon ( double pl )

inline

Definition at line 217 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

217{ p_kaon_ = pl;}

◆ p_kaon() [3/4]

double KalFitTrack::p_kaon ( void ) const

inline

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

199{ return p_kaon_; }

Referenced by KalFitAlg::complete_track().

◆ p_kaon() [4/4]

double KalFitTrack::p_kaon ( void ) const

inline

Definition at line 199 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

199{ return p_kaon_; }

◆ p_proton() [1/4]

void KalFitTrack::p_proton ( double pl )

inline

Definition at line 218 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

218{ p_proton_ = pl;}

◆ p_proton() [2/4]

void KalFitTrack::p_proton ( double pl )

inline

Definition at line 218 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

218{ p_proton_ = pl;}

◆ p_proton() [3/4]

double KalFitTrack::p_proton ( void ) const

inline

Definition at line 200 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

200{ return p_proton_; }

Referenced by KalFitAlg::complete_track().

◆ p_proton() [4/4]

double KalFitTrack::p_proton ( void ) const

inline

Definition at line 200 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

200{ return p_proton_; }

◆ pat1() [1/2]

int KalFitTrack::pat1 ( void ) const

inline

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

206{ return pat1_; }

◆ pat1() [2/2]

int KalFitTrack::pat1 ( void ) const

inline

Definition at line 206 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

206{ return pat1_; }

◆ pat2() [1/2]

int KalFitTrack::pat2 ( void ) const

inline

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

207{ return pat2_; }

◆ pat2() [2/2]

int KalFitTrack::pat2 ( void ) const

inline

Definition at line 207 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

207{ return pat2_; }

◆ path_ab() [1/2]

double KalFitTrack::path_ab ( void ) const

inline

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

193{ return path_ab_; }

◆ path_ab() [2/2]

double KalFitTrack::path_ab ( void ) const

inline

Definition at line 193 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

193{ return path_ab_; }

◆ path_add() [1/2]

void KalFitTrack::path_add ( double path )

Update the path length estimation.

Definition at line 1294 of file KalFitTrack.cxx.

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

◆ path_add() [2/2]

void KalFitTrack::path_add ( double path )

Update the path length estimation.

◆ path_rd() [1/2]

double KalFitTrack::path_rd ( void ) const

inline

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

192{ return path_rd_; }

◆ path_rd() [2/2]

double KalFitTrack::path_rd ( void ) const

inline

Definition at line 192 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

192{ return path_rd_; }

◆ pathip() [1/4]

void KalFitTrack::pathip ( double pl )

inline

Definition at line 216 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

216{ pathip_ = pl;}

◆ pathip() [2/4]

void KalFitTrack::pathip ( double pl )

inline

Definition at line 216 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

216{ pathip_ = pl;}

◆ pathip() [3/4]

double KalFitTrack::pathip ( void ) const

inline

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

191{ return pathip_; }

Referenced by KalFitAlg::fillTds_back().

◆ pathip() [4/4]

double KalFitTrack::pathip ( void ) const

inline

Definition at line 191 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

191{ return pathip_; }

◆ PathL() [1/2]

double KalFitTrack::PathL ( int layer )

Function to calculate the path length in the layer.

◆ PathL() [2/2]

double KalFitTrack::PathL ( int layer )

Function to calculate the path length in the layer.

◆ pathl() [1/2]

double * KalFitTrack::pathl ( void )

inline

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

194{ return PathL_; }

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

◆ pathl() [2/2]

double * KalFitTrack::pathl ( void )

inline

Definition at line 194 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

194{ return PathL_; }

◆ pivot_forMdc() [1/2]

const HepPoint3D & KalFitTrack::pivot_forMdc ( void ) const

inline

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

155{ return pivot_forMdc_;}

◆ pivot_forMdc() [2/2]

const HepPoint3D & KalFitTrack::pivot_forMdc ( void ) const

inline

Definition at line 155 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

155{ return pivot_forMdc_;}

◆ pivot_last() [1/2]

const HepPoint3D & KalFitTrack::pivot_last ( void ) const

inline

returns helix parameters

Definition at line 151 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

151{ return pivot_last_; }

◆ pivot_last() [2/2]

const HepPoint3D & KalFitTrack::pivot_last ( void ) const

inline

returns helix parameters

Definition at line 151 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

151{ return pivot_last_; }

◆ pivot_numf() [1/4]

const HepPoint3D & KalFitTrack::pivot_numf ( const HepPoint3D & newPivot )

Sets pivot position in a given mag field.

Definition at line 1379 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::Cgem_filter_anal(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), KalFitAlg::fitGemHits(), KalFitAlg::fromFHitToInnerWall(), KalFitAlg::innerwall(), KalFitAlg::smoother_anal(), KalFitAlg::smoother_calib(), update_hits(), KalFitElement::updateTrack(), KalFitElement::updateTrack_alreadyfound(), and KalFitElement::updateTrack_rphi().

◆ pivot_numf() [2/4]

const HepPoint3D & KalFitTrack::pivot_numf ( const HepPoint3D & newPivot )

Sets pivot position in a given mag field.

◆ pivot_numf() [3/4]

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

Definition at line 1483 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() [4/4]

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

◆ point_last() [1/4]

void KalFitTrack::point_last ( const HepPoint3D & point )

inline

set and give out the last point of the track

Definition at line 146 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

146{ point_last_ = point;}

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

◆ point_last() [2/4]

void KalFitTrack::point_last ( const HepPoint3D & point )

inline

set and give out the last point of the track

Definition at line 146 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

146{ point_last_ = point;}

◆ point_last() [3/4]

const HepPoint3D & KalFitTrack::point_last ( void )

inline

Definition at line 147 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

147{ return point_last_;}

◆ point_last() [4/4]

const HepPoint3D & KalFitTrack::point_last ( void )

inline

Definition at line 147 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

147{ return point_last_;}

◆ r0() [1/2]

double KalFitTrack::r0 ( void ) const

inline

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

186{ return r0_; }

◆ r0() [2/2]

double KalFitTrack::r0 ( void ) const

inline

Definition at line 186 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

186{ return r0_; }

◆ r_max() [1/2]

double KalFitTrack::r_max ( void ) const

inline

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

202{ return r_max_; }

◆ r_max() [2/2]

double KalFitTrack::r_max ( void ) const

inline

Definition at line 202 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

202{ return r_max_; }

◆ radius_numf() [1/2]

double KalFitTrack::radius_numf ( void ) const

Estimation of the radius in a given mag field.

Definition at line 1341 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(), and pivot_numf().

◆ radius_numf() [2/2]

double KalFitTrack::radius_numf ( void ) const

Estimation of the radius in a given mag field.

◆ resetLayerUsed() [1/2]

void KalFitTrack::resetLayerUsed ( )

inline

Definition at line 350 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ resetLayerUsed() [2/2]

void KalFitTrack::resetLayerUsed ( )

inline

Definition at line 350 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

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

◆ resol() [1/4]

void KalFitTrack::resol ( int i )

static

Definition at line 1836 of file KalFitTrack.cxx.

1836{ resolflag_ = i;}

◆ resol() [2/4]

static void KalFitTrack::resol ( int i )

static

◆ resol() [3/4]

int KalFitTrack::resol ( void )

static

Definition at line 1837 of file KalFitTrack.cxx.

1837{ return resolflag_;}

Referenced by KalFitAlg::initialize().

◆ resol() [4/4]

static int KalFitTrack::resol ( void )

static

◆ setClusterRefVec() [1/2]

void KalFitTrack::setClusterRefVec ( ClusterRefVec clusterRefVec )

inline

Definition at line 242 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

242{ ClusterRefVec_ = clusterRefVec; } //added hy Huang Zhen

◆ setClusterRefVec() [2/2]

void KalFitTrack::setClusterRefVec ( ClusterRefVec clusterRefVec )

inline

Definition at line 242 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

242{ ClusterRefVec_ = clusterRefVec; } //added hy Huang Zhen

◆ setIMdcGeomSvc() [1/2]

void KalFitTrack::setIMdcGeomSvc ( IMdcGeomSvc * igeomsvc )

static

Definition at line 239 of file KalFitTrack2.cxx.

{
    iGeomSvc_ = igeomsvc;
}

Referenced by KalFitAlg::setGeomSvc_init().

◆ setIMdcGeomSvc() [2/2]

static void KalFitTrack::setIMdcGeomSvc ( IMdcGeomSvc * igeomsvc )

static

◆ setInitMatrix() [1/2]

void KalFitTrack::setInitMatrix ( HepSymMatrix m )

static

Definition at line 37 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/2]

static void KalFitTrack::setInitMatrix ( HepSymMatrix m )

static

◆ setMagneticFieldSvc() [1/2]

void KalFitTrack::setMagneticFieldSvc ( IMagneticFieldSvc * mf )

static

Definition at line 228 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/2]

static void KalFitTrack::setMagneticFieldSvc ( IMagneticFieldSvc * )

static

◆ setMdcCalibFunSvc() [1/2]

void KalFitTrack::setMdcCalibFunSvc ( const MdcCalibFunSvc * calibsvc )

static

Definition at line 223 of file KalFitTrack2.cxx.

{
    CalibFunSvc_ = calibsvc;
}

Referenced by KalFitAlg::setCalibSvc_init().

◆ setMdcCalibFunSvc() [2/2]

static void KalFitTrack::setMdcCalibFunSvc ( const MdcCalibFunSvc * calibsvc )

static

◆ setMdcDigiCol() [1/2]

void KalFitTrack::setMdcDigiCol ( MdcDigiCol * digicol )

static

Definition at line 244 of file KalFitTrack2.cxx.

{
    mdcDigiCol_ = digicol;
}

Referenced by KalFitAlg::execute().

◆ setMdcDigiCol() [2/2]

static void KalFitTrack::setMdcDigiCol ( MdcDigiCol * digicol )

static

◆ setT0() [1/2]

void KalFitTrack::setT0 ( double t0 )

static

Definition at line 47 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/2]

static void KalFitTrack::setT0 ( double t0 )

static

◆ smoother_Mdc() [1/4]

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/4]

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

Kalman smoother for Mdc.

◆ smoother_Mdc() [3/4]

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

◆ smoother_Mdc() [4/4]

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

◆ smoother_Mdc_csmalign() [1/2]

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

Referenced by KalFitAlg::smoother_calib().

◆ smoother_Mdc_csmalign() [2/2]

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

◆ tof() [1/4]

void KalFitTrack::tof ( double path )

Update the tof estimation.

Definition at line 1316 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().

◆ tof() [2/4]

void KalFitTrack::tof ( double path )

Update the tof estimation.

◆ tof() [3/4]

double KalFitTrack::tof ( void ) const

inline

Definition at line 196 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

196{ return tof_; }

Referenced by msgasmdc(), and path_add().

◆ tof() [4/4]

double KalFitTrack::tof ( void ) const

inline

Definition at line 196 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

196{ return tof_; }

◆ tof_kaon() [1/2]

double KalFitTrack::tof_kaon ( void ) const

inline

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

197{ return tof_kaon_; }

Referenced by KalFitAlg::fillTds_back().

◆ tof_kaon() [2/2]

double KalFitTrack::tof_kaon ( void ) const

inline

Definition at line 197 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

197{ return tof_kaon_; }

◆ tof_proton() [1/2]

double KalFitTrack::tof_proton ( void ) const

inline

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

198{ return tof_proton_; }

Referenced by KalFitAlg::fillTds_back().

◆ tof_proton() [2/2]

double KalFitTrack::tof_proton ( void ) const

inline

Definition at line 198 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

198{ return tof_proton_; }

◆ trasan_id() [1/4]

void KalFitTrack::trasan_id ( int t )

inline

Definition at line 213 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

213{ trasan_id_ = t;}

◆ trasan_id() [2/4]

void KalFitTrack::trasan_id ( int t )

inline

Definition at line 213 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

213{ trasan_id_ = t;}

◆ trasan_id() [3/4]

int KalFitTrack::trasan_id ( void ) const

inline

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

185{ return trasan_id_; }

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

◆ trasan_id() [4/4]

int KalFitTrack::trasan_id ( void ) const

inline

Definition at line 185 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

185{ return trasan_id_; }

◆ type() [1/4]

void KalFitTrack::type ( int t )

inline

Reinitialize (modificator)

Definition at line 212 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

212{ type_ = t;}

◆ type() [2/4]

void KalFitTrack::type ( int t )

inline

Reinitialize (modificator)

Definition at line 212 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

212{ type_ = t;}

◆ type() [3/4]

int KalFitTrack::type ( void ) const

inline

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

184{ return type_; }

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

◆ type() [4/4]

int KalFitTrack::type ( void ) const

inline

Definition at line 184 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

184{ return type_; }

◆ update_bit() [1/2]

void KalFitTrack::update_bit ( int i )

Definition at line 1812 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/2]

void KalFitTrack::update_bit ( int i )

◆ update_forMdc() [1/2]

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/2]

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	(	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::Cgem_filter_anal(), KalFitAlg::filter_fwd_anal(), and KalFitAlg::filter_fwd_calib().

◆ 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.

◆ update_hits() [5/6]

double KalFitTrack::update_hits	(	RecCgemCluster *	Cluster,
		double	recR,
		int	way,
		int	csmflag
	)

Definition at line 2813 of file KalFitTrack.cxx.

{
    //const double a_stero[3] = {(45.94*3.1415926/180),(31.10*3.1415926/180),(32.99*3.1415926/180)};
    //const double r_layer[3] = {87.5026,132.7686,175.2686};
    //const double x_reso[3]={0.1372,0.1476,0.1412};
    //const double v_reso[3]={0.1273,0.1326,0.1378};
 
    // --- measurement
    double recZ   = Cluster->getRecZ()/10;
    double recPhi = Cluster->getrecphi();
    while(recPhi > M_PI) recPhi-=2*M_PI;
    while(recPhi <-M_PI) recPhi+=2*M_PI;
    HepVector v_measu(2,0);
    v_measu(1) = recPhi;
    v_measu(2) = recZ;
 
    // --- estimation
    //double Phi = intersect_cylinder(recR);
    HepPoint3D x0kal = x(intersect_cylinder(recR));
    pivot_numf(x0kal);
    HepVector v_estim(2,0);
    v_estim(1) = x0kal.phi();
    v_estim(2) = x0kal.z();
 
    const HepSymMatrix& ea = Ea();
    HepVector v_a = a();
 
    // --- derivative matrix
    double dPhi=intersect_cylinder(recR);
    double dr = v_a(1);
    double phi0 = v_a(2);
    double kappa = v_a(3);
    double tanl = v_a(5);
    double x0 = x0kal.x();
    double y0 = x0kal.y();
    double Alpha = alpha();
    HepMatrix H(2, 5, 0);
    H(1,1) = -y0*cos(phi0)/(y0*y0+x0*x0)+x0*sin(phi0)/(x0*x0+y0*y0);
    H(1,2) = -y0/(y0*y0+x0*x0)*((-1)*dr*sin(phi0)+Alpha/kappa*(sin(phi0+dPhi)-sin(phi0)))+x0/(x0*x0+y0*y0)*(dr*cos(phi0)+Alpha/kappa*(cos(phi0)-cos(phi0+dPhi)));
    H(1,3) = -y0/(y0*y0+x0*x0)*(-1)*Alpha/(kappa*kappa)*(cos(phi0)-cos(phi0+dPhi))+x0/(x0*x0+y0*y0)*(-1)*Alpha/(kappa*kappa)*(sin(phi0)-sin(phi0+dPhi));
    H(2,3) = Alpha/(kappa*kappa)*tanl*dPhi;
    H(2,4) = 1.0;
    H(2,5) = -1*(Alpha/kappa)*dPhi;
 
    // --- error matrix of Cgem
    int layer  = Cluster->getlayerid();
    ISvcLocator* svcLocator = Gaudi::svcLocator();
    ICgemGeomSvc* ISvc;
    StatusCode Cgem_sc=svcLocator->service("CgemGeomSvc", ISvc);
    CgemGeomSvc* CgemGeomSvc_=dynamic_cast<CgemGeomSvc *>(ISvc);
    CgemGeoLayer* CgemLayer = CgemGeomSvc_->getCgemLayer(layer);
    double R_x = (CgemLayer->getInnerROfAnodeCu2() + CgemLayer->getOuterROfAnodeCu2())/2;
    //R_v = (CgemLayer->getInnerROfAnodeCu1() + CgemLayer->getOuterROfAnodeCu1())/2;
    double a_stero = (CgemLayer->getAngleOfStereo())*M_PI/180;
 
    ICgemCalibFunSvc* CgemCalibSvc;
    StatusCode sc = svcLocator->service ("CgemCalibFunSvc", CgemCalibSvc);
    int iView(0), mode(2);
    double Q(100), T(100);
    double Phi_momentum = momentum(dPhi).phi();
    double Phi_position = x0kal.phi();
    double delta_phi=Phi_momentum - Phi_position;
    while(delta_phi>M_PI)  delta_phi-=CLHEP::twopi;
    while(delta_phi<-M_PI) delta_phi+=CLHEP::twopi;
    double sigma_X = CgemCalibSvc->getSigma(layer,iView,mode,delta_phi,Q,T);// in mm
    double sigma_V = CgemCalibSvc->getSigma(layer,iView,mode,delta_phi,Q,T);// in mm
 
    HepSymMatrix V(2,0);
    //V(1,1) = pow(x_reso[layer]/r_layer[layer],2);
    //V(2,2) = pow(sqrt(0.01*v_reso[layer]*v_reso[layer]/(sin(a_stero[layer])*sin(a_stero[layer]))+0.01*x_reso[layer]*x_reso[layer]/(tan(a_stero[layer])*tan(a_stero[layer]))),2);
    V(1,1) = pow(sigma_X/R_x,2);
    V(2,2) = pow(sqrt(0.01*sigma_V*sigma_V/(sin(a_stero)*sin(a_stero))+0.01*sigma_X*sigma_X/(tan(a_stero)*tan(a_stero))),2);
    // --- gain matrix
    int ierr=-1;
    HepMatrix K = ea*H.T()*(V+H*ea*H.T()).inverse(ierr);
    if(ierr != 0)cout<<"errer in inverse operation of matrix!"<<endl;
 
    // --- new error matrix
    HepSymMatrix EaNew(5,0);
    EaNew.assign(ea-K*H*ea);
    Ea(EaNew);
 
    // --- diff
    HepVector v_diff = v_measu - v_estim;
    while(v_diff(1) > M_PI) v_diff(1)-=2*M_PI;
    while(v_diff(1) <- M_PI) v_diff(1)+=2*M_PI;
 
    // --- new parameters
    //HepVector aNew = v_a + K*v_diff;
    a(v_a + K*v_diff);
 
    // --- new v_estim
    HepPoint3D x0kal_new = x(intersect_cylinder(recR));
    HepVector v_estim_new(2,0);
    v_estim_new(1) = x0kal_new.phi();
    v_estim_new(2) = x0kal_new.z();
 
    // --- difference between measurement and updated estimation
    v_diff = v_measu - v_estim_new;
    while(v_diff(1) > M_PI) v_diff(1)-=2*M_PI;
    while(v_diff(1) <- M_PI) v_diff(1)+=2*M_PI;
 
    // --- new derivative matrix 
    //track_new.pivot_numf(x0kal_new,pathl);
    //HepVector a_new = track_new.a();
    //HepVector Ea_new = track_new.ea();
 
    // --- R matrix and Chisuqre 
    HepSymMatrix R(2,0);
    R.assign(V-H*EaNew*H.T());
    HepVector dChi2 = v_diff.T()*R.inverse(ierr)*v_diff;
    if(ierr != 0)cout<<"errer in inverse operation of matrix!"<<endl;
 
    // --- update Chisqure of track
    if(way>0) chiSq(chiSq()+dChi2(1));
    if(way<0) chiSq_back(chiSq_back()+dChi2(1));
    return dChi2(1);
}

◆ update_hits() [6/6]

double KalFitTrack::update_hits	(	RecCgemCluster *	Cluster,
		double	recR,
		int	way,
		int	csmflag
	)

◆ update_hits_csmalign() [1/2]

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/2]

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

◆ update_last() [1/2]

void KalFitTrack::update_last ( void )

Record the current parameters as ..._last information :

Definition at line 114 of file KalFitTrack.cxx.

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

◆ update_last() [2/2]

void KalFitTrack::update_last ( void )

Record the current parameters as ..._last information :

◆ useLayer() [1/2]

void KalFitTrack::useLayer ( int iLay )

inline

Definition at line 354 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ useLayer() [2/2]

void KalFitTrack::useLayer ( int iLay )

inline

Definition at line 354 of file Reconstruction/KalFitAlg/KalFitAlg-00-15-08/KalFitAlg/KalFitTrack.h.

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

Member Data Documentation

◆ Bznom_

double KalFitTrack::Bznom_ = 10

static

Definition at line 306 of file InstallArea/include/KalFitAlg/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(), and radius_numf().

◆ chi2_hitf_

double KalFitTrack::chi2_hitf_ = 1000

static

Cut chi2 for each hit.

Definition at line 289 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize().

◆ chi2_hits_

double KalFitTrack::chi2_hits_ = 1000

static

Definition at line 289 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize().

◆ chi2mdc_hit2_

static double KalFitTrack::chi2mdc_hit2_

static

Definition at line 320 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

◆ dchi2cutf_anal

double KalFitTrack::dchi2cutf_anal = {0.}

static

Definition at line 295 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut().

◆ dchi2cutf_calib

double KalFitTrack::dchi2cutf_calib = {0.}

static

Definition at line 299 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut().

◆ dchi2cuts_anal

double KalFitTrack::dchi2cuts_anal = {0.}

static

Definition at line 297 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut().

◆ dchi2cuts_calib

double KalFitTrack::dchi2cuts_calib = {0.}

static

Definition at line 301 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::setDchisqCut().

◆ debug_

int KalFitTrack::debug_ = 0

static

for debug

Definition at line 287 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ drifttime_choice_

int KalFitTrack::drifttime_choice_ = 0

static

the drifttime choice

Definition at line 328 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ factor_strag_

double KalFitTrack::factor_strag_ = 0.4

static

factor of energy loss straggling for electron

Definition at line 317 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ inner_steps_

int KalFitTrack::inner_steps_ = 3

static

Definition at line 292 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ LR_

int KalFitTrack::LR_ = 1

static

Use L/R decision from MdcRecHit information :

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

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

◆ mdcGasRadlen_

double KalFitTrack::mdcGasRadlen_ = 0.

static

Definition at line 281 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ nmdc_hit2_

int KalFitTrack::nmdc_hit2_ = 500

static

Cut chi2 for each hit.

Definition at line 319 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize().

◆ numf_

int KalFitTrack::numf_ = 0

static

Flag for treatment of non-uniform mag field.

Definition at line 291 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ numfcor_

int KalFitTrack::numfcor_ = 1

static

NUMF treatment improved.

Definition at line 304 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ outer_steps_

int KalFitTrack::outer_steps_ = 3

static

Definition at line 293 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ resolflag_

int KalFitTrack::resolflag_ = 0

static

wire resoltion flag

Definition at line 323 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by chi2_next(), and resol().

◆ 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 309 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

Referenced by KalFitAlg::initialize().

◆ strag_

int KalFitTrack::strag_ = 1

static

Flag to take account of energy loss straggling :

Definition at line 315 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ Tof_correc_

int KalFitTrack::Tof_correc_ = 1

static

Flag for TOF correction.

Definition at line 312 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ tofall_

int KalFitTrack::tofall_ = 1

static

Definition at line 321 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

◆ tprop_

int KalFitTrack::tprop_ = 1

static

for signal propagation correction

Definition at line 284 of file InstallArea/include/KalFitAlg/KalFitAlg/KalFitTrack.h.

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

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

Public Member Functions

Static Public Member Functions

Static Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ KalFitTrack() [1/2]

◆ ~KalFitTrack() [1/2]

◆ KalFitTrack() [2/2]

◆ ~KalFitTrack() [2/2]

Member Function Documentation

◆ a_forMdc() [1/2]

◆ a_forMdc() [2/2]

◆ a_last() [1/2]

◆ a_last() [2/2]

◆ add_nhit_r() [1/2]

◆ add_nhit_r() [2/2]

◆ add_nhit_z() [1/2]

◆ add_nhit_z() [2/2]

◆ addPathSM() [1/2]

◆ addPathSM() [2/2]

◆ addTofSM() [1/2]

◆ addTofSM() [2/2]

◆ appendHelixSegs() [1/2]

◆ appendHelixSegs() [2/2]

◆ appendHitsMdc() [1/2]

◆ appendHitsMdc() [2/2]

◆ back() [1/4]

◆ back() [2/4]

◆ back() [3/4]

◆ back() [4/4]

◆ chgmass() [1/2]

◆ chgmass() [2/2]

◆ chi2_next() [1/4]

◆ chi2_next() [2/4]

◆ chi2_next() [3/4]

◆ chi2_next() [4/4]

◆ chiSq() [1/4]

◆ chiSq() [2/4]

◆ chiSq() [3/4]

◆ chiSq() [4/4]

◆ chiSq_back() [1/4]

◆ chiSq_back() [2/4]

◆ chiSq_back() [3/4]

◆ chiSq_back() [4/4]

◆ cor_tanldep() [1/2]

◆ cor_tanldep() [2/2]

◆ dchi2_max() [1/2]

◆ dchi2_max() [2/2]

◆ Ea_forMdc() [1/2]

◆ Ea_forMdc() [2/2]

◆ Ea_last() [1/2]

◆ Ea_last() [2/2]

◆ eloss() [1/2]

◆ eloss() [2/2]

◆ filter() [1/2]

◆ filter() [2/2]

◆ fiTerm() [1/2]

◆ fiTerm() [2/2]

◆ getClusterRefVec() [1/2]

◆ getClusterRefVec() [2/2]

◆ getDigi() [1/2]

◆ getDigi() [2/2]

◆ getDriftDist() [1/2]

◆ getDriftDist() [2/2]

◆ getDriftTime() [1/2]

◆ getDriftTime() [2/2]

◆ getFiTerm() [1/2]

◆ getFiTerm() [2/2]

◆ getInitMatrix() [1/2]

◆ getInitMatrix() [2/2]

◆ getPathSM() [1/2]

◆ getPathSM() [2/2]

◆ getSigma() [1/4]

◆ getSigma() [2/4]

◆ getSigma() [3/4]

◆ getSigma() [4/4]

◆ getT0() [1/2]

◆ getT0() [2/2]

◆ getTofSM() [1/2]

◆ getTofSM() [2/2]