TRunge Class Reference

A class to represent a track in tracking. More...

#include <TRunge.h>

Inheritance diagram for TRunge:

Public Member Functions
	TRunge ()
	Constructors.
	TRunge (const TTrack &)
	TRunge (const Helix &)
	TRunge (const TRunge &)
	TRunge (const RecMdcTrack &)
	~TRunge ()
	Destructor.
unsigned	objectType (void) const
	returns object type
double	dr (void) const
	Track parameters (at pivot).
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
const HepPoint3D &	pivot (void) const
	pivot position
const Vector &	a (void) const
	returns helix parameter
const SymMatrix &	Ea (void) const
	returns error matrix
Helix	helix (void) const
	returns helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
int	BfieldID (void) const
	returns B field ID
double	StepSize (void) const
	returns step size
const double *	Yscal (void) const
	return error parameters for fitting with step size control
double	Eps (void) const
double	StepSizeMax (void) const
double	StepSizeMin (void) const
float	Mass (void) const
	return mass
double	MaxFlightLength (void) const
	return max flight length
int	approach (TMLink &, const RkFitMaterial, bool sagCorrection=true)
int	approach (TMLink &, float &tof, HepVector3D &p, const RkFitMaterial, bool sagCorrection=true)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest points between a line and this track
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material, unsigned &stepNum)
int	approach_point (TMLink &, const HepPoint3D &, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
const Vector &	a (const Vector &)
	set helix parameter
const SymMatrix &	Ea (const SymMatrix &)
	set helix error matrix
int	BfieldID (int)
	set B field map ID
double	StepSize (double)
	set step size to calc. trajectory
const double *	Yscal (const double *)
	set error parameters for fitting with step size control
double	Eps (double)
double	StepSizeMax (double)
double	StepSizeMin (double)
float	Mass (float)
	set mass for tof calc.
double	MaxFlightLength (double)
double	DisToWire (TMLink &, double, HepPoint3D, HepPoint3D &)
double	dEpath (double, double, double) const
void	eloss (double path, const RkFitMaterial *material, double mass, double y[6], int index) const
unsigned	Fly (const RkFitMaterial material)
	make the trajectory in cache, return the number of step
unsigned	Fly_SC (void)
double	intersect_cylinder (double r) const
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	Propagate (double y[6], const double &step, const RkFitMaterial material)
	propagate the track using 4th-order Runge-Kutta method
void	Function (const double y[6], double f[6])
void	Propagate1 (const double y[6], const double dydx[6], const double &step, double yout[6])
void	Propagate_QC (double y[6], double dydx[6], const double &steptry, const double &eps, const double yscal[6], double &stepdid, double &stepnext)
void	SetFirst (double y[6]) const
	set first point (position, momentum) s=0, phi=0
unsigned	Nstep (void) const
	access to trajectory cache
int	GetXP (unsigned stepNum, double y[6]) const
int	GetStep (unsigned stepNum, double &step) const
double	SetFlightLength (void)
	set flight length using wire hits
	TRunge ()
	Constructors.
	TRunge (const TTrack &)
	TRunge (const Helix &)
	TRunge (const TRunge &)
	TRunge (const RecMdcTrack &)
	~TRunge ()
	Destructor.
unsigned	objectType (void) const
	returns object type
double	dr (void) const
	Track parameters (at pivot).
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
const HepPoint3D &	pivot (void) const
	pivot position
const Vector &	a (void) const
	returns helix parameter
const SymMatrix &	Ea (void) const
	returns error matrix
Helix	helix (void) const
	returns helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
int	BfieldID (void) const
	returns B field ID
double	StepSize (void) const
	returns step size
const double *	Yscal (void) const
	return error parameters for fitting with step size control
double	Eps (void) const
double	StepSizeMax (void) const
double	StepSizeMin (void) const
float	Mass (void) const
	return mass
double	MaxFlightLength (void) const
	return max flight length
int	approach (TMLink &, const RkFitMaterial, bool sagCorrection=true)
int	approach (TMLink &, float &tof, HepVector3D &p, const RkFitMaterial, bool sagCorrection=true)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest points between a line and this track
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material, unsigned &stepNum)
int	approach_point (TMLink &, const HepPoint3D &, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
const Vector &	a (const Vector &)
	set helix parameter
const SymMatrix &	Ea (const SymMatrix &)
	set helix error matrix
int	BfieldID (int)
	set B field map ID
double	StepSize (double)
	set step size to calc. trajectory
const double *	Yscal (const double *)
	set error parameters for fitting with step size control
double	Eps (double)
double	StepSizeMax (double)
double	StepSizeMin (double)
float	Mass (float)
	set mass for tof calc.
double	MaxFlightLength (double)
double	DisToWire (TMLink &, double, HepPoint3D, HepPoint3D &)
double	dEpath (double, double, double) const
void	eloss (double path, const RkFitMaterial *material, double mass, double y[6], int index) const
unsigned	Fly (const RkFitMaterial material)
	make the trajectory in cache, return the number of step
unsigned	Fly_SC (void)
double	intersect_cylinder (double r) const
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	Propagate (double y[6], const double &step, const RkFitMaterial material)
	propagate the track using 4th-order Runge-Kutta method
void	Function (const double y[6], double f[6])
void	Propagate1 (const double y[6], const double dydx[6], const double &step, double yout[6])
void	Propagate_QC (double y[6], double dydx[6], const double &steptry, const double &eps, const double yscal[6], double &stepdid, double &stepnext)
void	SetFirst (double y[6]) const
	set first point (position, momentum) s=0, phi=0
unsigned	Nstep (void) const
	access to trajectory cache
int	GetXP (unsigned stepNum, double y[6]) const
int	GetStep (unsigned stepNum, double &step) const
double	SetFlightLength (void)
	set flight length using wire hits
	TRunge ()
	Constructors.
	TRunge (const TTrack &)
	TRunge (const Helix &)
	TRunge (const TRunge &)
	TRunge (const RecMdcTrack &)
	~TRunge ()
	Destructor.
unsigned	objectType (void) const
	returns object type
double	dr (void) const
	Track parameters (at pivot).
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
const HepPoint3D &	pivot (void) const
	pivot position
const Vector &	a (void) const
	returns helix parameter
const SymMatrix &	Ea (void) const
	returns error matrix
Helix	helix (void) const
	returns helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
int	BfieldID (void) const
	returns B field ID
double	StepSize (void) const
	returns step size
const double *	Yscal (void) const
	return error parameters for fitting with step size control
double	Eps (void) const
double	StepSizeMax (void) const
double	StepSizeMin (void) const
float	Mass (void) const
	return mass
double	MaxFlightLength (void) const
	return max flight length
int	approach (TMLink &, const RkFitMaterial, bool sagCorrection=true)
int	approach (TMLink &, float &tof, HepVector3D &p, const RkFitMaterial, bool sagCorrection=true)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest points between a line and this track
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material, unsigned &stepNum)
int	approach_point (TMLink &, const HepPoint3D &, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
const Vector &	a (const Vector &)
	set helix parameter
const SymMatrix &	Ea (const SymMatrix &)
	set helix error matrix
int	BfieldID (int)
	set B field map ID
double	StepSize (double)
	set step size to calc. trajectory
const double *	Yscal (const double *)
	set error parameters for fitting with step size control
double	Eps (double)
double	StepSizeMax (double)
double	StepSizeMin (double)
float	Mass (float)
	set mass for tof calc.
double	MaxFlightLength (double)
double	DisToWire (TMLink &, double, HepPoint3D, HepPoint3D &)
double	dEpath (double, double, double) const
void	eloss (double path, const RkFitMaterial *material, double mass, double y[6], int index) const
unsigned	Fly (const RkFitMaterial material)
	make the trajectory in cache, return the number of step
unsigned	Fly_SC (void)
double	intersect_cylinder (double r) const
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	Propagate (double y[6], const double &step, const RkFitMaterial material)
	propagate the track using 4th-order Runge-Kutta method
void	Function (const double y[6], double f[6])
void	Propagate1 (const double y[6], const double dydx[6], const double &step, double yout[6])
void	Propagate_QC (double y[6], double dydx[6], const double &steptry, const double &eps, const double yscal[6], double &stepdid, double &stepnext)
void	SetFirst (double y[6]) const
	set first point (position, momentum) s=0, phi=0
unsigned	Nstep (void) const
	access to trajectory cache
int	GetXP (unsigned stepNum, double y[6]) const
int	GetStep (unsigned stepNum, double &step) const
double	SetFlightLength (void)
	set flight length using wire hits
Public Member Functions inherited from TTrackBase
	TTrackBase ()
	Constructor.
	TTrackBase (const AList< TMLink > &links)
	Constructor.
virtual	~TTrackBase ()
	Destructor.
virtual unsigned	type (void) const
	returns type. Definition is depending on an object class.
virtual void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const AList< TMLink > &	links (unsigned mask=0) const
unsigned	nLinks (unsigned mask=0) const
	returns # of masked TMLinks assigned to this track object.
const AList< TMLink > &	cores (unsigned mask=0) const
	returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
unsigned	nCores (unsigned mask=0) const
	returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
void	update (void) const
	update cache.
void	append (TMLink &)
	appends a TMLink.
void	append (const AList< TMLink > &)
	appends TMLinks.
void	appendByApproach (AList< TMLink > &list, double maxSigma)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
virtual void	refine (double maxSigma)
	removes bad points by pull. The bad points are masked not to be used in fit.
virtual int	DropWorst ()
virtual void	removeLinks (void)
virtual double	distance (const TMLink &) const
	returns distance to a position of TMLink in TMLink space.
virtual int	approach (TMLink &) const
unsigned	testByApproach (const TMLink &list, double sigma) const
	returns # of good hits to be appended.
unsigned	testByApproach (const AList< TMLink > &list, double sigma) const
virtual int	fit (void)
	fits itself by a default fitter. Error was happened if return value is not zero.
const TMFitter *const	fitter (void) const
	returns a pointer to a default fitter.
const TMFitter *const	fitter (const TMFitter *)
	sets a default fitter.
void	falseFit ()
	false Fit
TMLink *	operator[] (unsigned i) const
const TTrackHEP *const	hep (void) const
	returns TTrackHEP.
unsigned	nHeps (void) const
	returns # of contributed TTrackHEP tracks.
const TTrackMC *const	mc (void) const
	returns a pointer to TTrackMC.
bool	fitted (void) const
	returns true if fitted.
bool	fittedWithCathode (void) const
	returns true if fitted with cathode hits(TEMPORARY).
	TTrackBase ()
	Constructor.
	TTrackBase (const AList< TMLink > &links)
	Constructor.
virtual	~TTrackBase ()
	Destructor.
virtual unsigned	type (void) const
	returns type. Definition is depending on an object class.
virtual void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const AList< TMLink > &	links (unsigned mask=0) const
unsigned	nLinks (unsigned mask=0) const
	returns # of masked TMLinks assigned to this track object.
const AList< TMLink > &	cores (unsigned mask=0) const
	returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
unsigned	nCores (unsigned mask=0) const
	returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
void	update (void) const
	update cache.
void	append (TMLink &)
	appends a TMLink.
void	append (const AList< TMLink > &)
	appends TMLinks.
void	appendByApproach (AList< TMLink > &list, double maxSigma)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
virtual void	refine (double maxSigma)
	removes bad points by pull. The bad points are masked not to be used in fit.
virtual int	DropWorst ()
virtual void	removeLinks (void)
virtual double	distance (const TMLink &) const
	returns distance to a position of TMLink in TMLink space.
virtual int	approach (TMLink &) const
unsigned	testByApproach (const TMLink &list, double sigma) const
	returns # of good hits to be appended.
unsigned	testByApproach (const AList< TMLink > &list, double sigma) const
virtual int	fit (void)
	fits itself by a default fitter. Error was happened if return value is not zero.
const TMFitter *const	fitter (void) const
	returns a pointer to a default fitter.
const TMFitter *const	fitter (const TMFitter *)
	sets a default fitter.
void	falseFit ()
	false Fit
TMLink *	operator[] (unsigned i) const
const TTrackHEP *const	hep (void) const
	returns TTrackHEP.
unsigned	nHeps (void) const
	returns # of contributed TTrackHEP tracks.
const TTrackMC *const	mc (void) const
	returns a pointer to TTrackMC.
bool	fitted (void) const
	returns true if fitted.
bool	fittedWithCathode (void) const
	returns true if fitted with cathode hits(TEMPORARY).
	TTrackBase ()
	Constructor.
	TTrackBase (const AList< TMLink > &links)
	Constructor.
virtual	~TTrackBase ()
	Destructor.
virtual unsigned	type (void) const
	returns type. Definition is depending on an object class.
virtual void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
const AList< TMLink > &	links (unsigned mask=0) const
unsigned	nLinks (unsigned mask=0) const
	returns # of masked TMLinks assigned to this track object.
const AList< TMLink > &	cores (unsigned mask=0) const
	returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
unsigned	nCores (unsigned mask=0) const
	returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0.
void	update (void) const
	update cache.
void	append (TMLink &)
	appends a TMLink.
void	append (const AList< TMLink > &)
	appends TMLinks.
void	appendByApproach (AList< TMLink > &list, double maxSigma)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
virtual void	refine (double maxSigma)
	removes bad points by pull. The bad points are masked not to be used in fit.
virtual int	DropWorst ()
virtual void	removeLinks (void)
virtual double	distance (const TMLink &) const
	returns distance to a position of TMLink in TMLink space.
virtual int	approach (TMLink &) const
unsigned	testByApproach (const TMLink &list, double sigma) const
	returns # of good hits to be appended.
unsigned	testByApproach (const AList< TMLink > &list, double sigma) const
virtual int	fit (void)
	fits itself by a default fitter. Error was happened if return value is not zero.
const TMFitter *const	fitter (void) const
	returns a pointer to a default fitter.
const TMFitter *const	fitter (const TMFitter *)
	sets a default fitter.
void	falseFit ()
	false Fit
TMLink *	operator[] (unsigned i) const
const TTrackHEP *const	hep (void) const
	returns TTrackHEP.
unsigned	nHeps (void) const
	returns # of contributed TTrackHEP tracks.
const TTrackMC *const	mc (void) const
	returns a pointer to TTrackMC.
bool	fitted (void) const
	returns true if fitted.
bool	fittedWithCathode (void) const
	returns true if fitted with cathode hits(TEMPORARY).

Friends
class	TRungeFitter

Additional Inherited Members
Protected Attributes inherited from TTrackBase
AList< TMLink >	_links
bool	_fitted
bool	_fittedWithCathode
TTrackMC *	_mc

Detailed Description

A class to represent a track in tracking.

Definition at line 65 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TRunge.h.

Constructor & Destructor Documentation

TRunge() [1/15]

TRunge::TRunge

(

)

Constructors.

Definition at line 46 of file TRunge.cxx.

    : TTrackBase()
    , _pivot( ORIGIN )
    , _a( Vector( 5, 0 ) )
    , _Ea( SymMatrix( 5, 0 ) )
    , _chi2( 0 )
    , _ndf( 0 )
    , _charge( 1 )
    , _bfieldID( 21 )
    , _Nstep( 0 )
    , _stepSize( default_stepSize )
    , _mass( 0.140 )
    , _eps( EPS )
    , _stepSizeMax( default_stepSizeMax )
    , _stepSizeMin( default_stepSizeMin )
    , _maxflightlength( default_maxflightlength )
    , m_pmgnIMF( nullptr ) {
 
  //...Set a default fitter...
  fitter( &TRunge::_fitter );
 
  _fitted            = false;
  _fittedWithCathode = false;
 
  _bfield = Bfield::getBfield( _bfieldID );
 
  _mass2 = _mass * _mass;
 
  _yscal[0] = 0.1;   // 1mm  -> error = 1mm*EPS
  _yscal[1] = 0.1;   // 1mm
  _yscal[2] = 0.1;   // 1mm
  _yscal[3] = 0.001; // 1MeV
  _yscal[4] = 0.001; // 1MeV
  _yscal[5] = 0.001; // 1MeV
}

Referenced by TRunge().

TRunge() [2/15]

TRunge::TRunge

(

const TTrack &

a

)

Definition at line 122 of file TRunge.cxx.

    : TTrackBase( a )
    , _pivot( a.helix().pivot() )
    , _a( a.helix().a() )
    , _Ea( a.helix().Ea() )
    , _chi2( 0 )
    , _ndf( 0 )
    , _bfieldID( 21 )
    , _Nstep( 0 )
    , _stepSize( default_stepSize )
    , _mass( 0.140 )
    , _eps( EPS )
    , _stepSizeMax( default_stepSizeMax )
    , _stepSizeMin( default_stepSizeMin )
    , _maxflightlength( default_maxflightlength )
    , m_pmgnIMF( nullptr ) {
  //  _a[2]=-_a[2];
  //...Set a default fitter...
  fitter( &TRunge::_fitter );
 
  _fitted            = false;
  _fittedWithCathode = false;
 
  if ( _a[2] < 0 ) _charge = -1;
  else _charge = 1;
 
  _bfield = Bfield::getBfield( _bfieldID );
 
  _mass2 = _mass * _mass;
 
  _yscal[0] = 0.1;   // 1mm  -> error = 1mm*EPS
  _yscal[1] = 0.1;   // 1mm
  _yscal[2] = 0.1;   // 1mm
  _yscal[3] = 0.001; // 1MeV
  _yscal[4] = 0.001; // 1MeV
  _yscal[5] = 0.001; // 1MeV
 
  // SetFlightLength();
  // cout<<"TR:: _maxflightlength="<<_maxflightlength<<endl;
}

TRunge() [3/15]

TRunge::TRunge

(

const Helix &

h

)

Definition at line 163 of file TRunge.cxx.

    : TTrackBase()
    , _pivot( h.pivot() )
    , _a( h.a() )
    , _Ea( h.Ea() )
    , _chi2( 0 )
    , _ndf( 0 )
    , _bfieldID( 21 )
    , _Nstep( 0 )
    , _stepSize( default_stepSize )
    , _mass( 0.140 )
    , _eps( EPS )
    , _stepSizeMax( default_stepSizeMax )
    , _stepSizeMin( default_stepSizeMin )
    , _maxflightlength( default_maxflightlength )
    , m_pmgnIMF( nullptr ) {
 
  //...Set a default fitter...
  fitter( &TRunge::_fitter );
 
  _fitted            = false;
  _fittedWithCathode = false;
 
  if ( _a[2] < 0 ) _charge = -1;
  else _charge = 1;
 
  _bfield = Bfield::getBfield( _bfieldID );
 
  _mass2 = _mass * _mass;
 
  _yscal[0] = 0.1;   // 1mm  -> error = 1mm*EPS
  _yscal[1] = 0.1;   // 1mm
  _yscal[2] = 0.1;   // 1mm
  _yscal[3] = 0.001; // 1MeV
  _yscal[4] = 0.001; // 1MeV
  _yscal[5] = 0.001; // 1MeV
}

TRunge() [4/15]

TRunge::TRunge

(

const TRunge &

a

)

Definition at line 201 of file TRunge.cxx.

    : TTrackBase( (TTrackBase&)a )
    , _pivot( a.pivot() )
    , _a( a.a() )
    , _Ea( a.Ea() )
    , _chi2( a.chi2() )
    , _ndf( a.ndf() )
    , _bfieldID( a.BfieldID() )
    , _Nstep( 0 )
    , _stepSize( a.StepSize() )
    , _mass( a.Mass() )
    , _eps( a.Eps() )
    , _stepSizeMax( a.StepSizeMax() )
    , _stepSizeMin( a.StepSizeMin() )
    , _maxflightlength( a.MaxFlightLength() )
    , m_pmgnIMF( nullptr ) {
 
  //...Set a default fitter...
  fitter( &TRunge::_fitter );
 
  _fitted            = false;
  _fittedWithCathode = false;
 
  if ( _a[2] < 0 ) _charge = -1;
  else _charge = 1;
 
  _bfield = Bfield::getBfield( _bfieldID );
 
  _mass2 = _mass * _mass;
 
  for ( unsigned i = 0; i < 6; i++ ) _yscal[i] = a.Yscal()[i];
}

TRunge() [5/15]

TRunge::TRunge

(

const RecMdcTrack &

a

)

Definition at line 82 of file TRunge.cxx.

    : TTrackBase()
    , _pivot( a.getPivot() )
    , _a( a.helix() )
    , _Ea( a.err() )
    , _chi2( a.chi2() )
    , _ndf( a.ndof() )
    , _bfieldID( 21 )
    , _Nstep( 0 )
    , _stepSize( default_stepSize )
    , _mass( 0.140 )
    , _eps( EPS )
    , _stepSizeMax( default_stepSizeMax )
    , _stepSizeMin( default_stepSizeMin )
    , _maxflightlength( default_maxflightlength )
    , m_pmgnIMF( nullptr ) {
  //    _a[2]=-_a[2];
  //...Set a default fitter...
  fitter( &TRunge::_fitter );
 
  _fitted            = false;
  _fittedWithCathode = false;
 
  if ( _a[2] < 0 ) _charge = -1;
  else _charge = 1;
 
  _bfield = Bfield::getBfield( _bfieldID );
 
  _mass2 = _mass * _mass;
 
  _yscal[0] = 0.1;   // 1mm  -> error = 1mm*EPS
  _yscal[1] = 0.1;   // 1mm
  _yscal[2] = 0.1;   // 1mm
  _yscal[3] = 0.001; // 1MeV
  _yscal[4] = 0.001; // 1MeV
  _yscal[5] = 0.001; // 1MeV
  // SetFlightLength();
  // cout<<"TR:: _maxflightlength="<<_maxflightlength<<endl;
}

~TRunge() [1/3]

TRunge::~TRunge

(

)

Destructor.

Definition at line 235 of file TRunge.cxx.

{}

TRunge() [6/15]

TRunge::TRunge

(

)

Constructors.

TRunge() [7/15]

TRunge::TRunge

(

const TTrack &

)

TRunge() [8/15]

TRunge::TRunge

(

const Helix &

)

TRunge() [9/15]

TRunge::TRunge

(

const TRunge &

)

TRunge() [10/15]

TRunge::TRunge

(

const RecMdcTrack &

)

~TRunge() [2/3]

TRunge::~TRunge

(

)

Destructor.

TRunge() [11/15]

TRunge::TRunge

(

)

Constructors.

TRunge() [12/15]

TRunge::TRunge

(

const TTrack &

)

TRunge() [13/15]

TRunge::TRunge

(

const Helix &

)

TRunge() [14/15]

TRunge::TRunge

(

const TRunge &

)

TRunge() [15/15]

TRunge::TRunge

(

const RecMdcTrack &

)

~TRunge() [3/3]

TRunge::~TRunge

(

)

Destructor.

Member Function Documentation

a() [1/6]

const Vector & TRunge::a

(

const Vector &

ta

)

set helix parameter

Definition at line 293 of file TRunge.cxx.

                                          {
  _a = ta;
  if ( _a[2] < 0 ) _charge = -1;
  else _charge = 1;
  _Nstep = 0;
  return ( _a );
}

a() [2/6]

const Vector & TRunge::a

(

const Vector &

)

set helix parameter

a() [3/6]

const Vector & TRunge::a

(

const Vector &

)

set helix parameter

a() [4/6]

const Vector & TRunge::a

(

void

)

const

returns helix parameter

Definition at line 249 of file TRunge.cxx.

{ return ( _a ); }

Referenced by dEpath(), DisToWire(), TRunge(), TRunge(), TRunge(), and TRunge().

a() [5/6]

const Vector & TRunge::a

(

void

)

const

returns helix parameter

a() [6/6]

const Vector & TRunge::a

(

void

)

const

returns helix parameter

approach() [1/6]

int TRunge::approach	(	TMLink &	l,
		const RkFitMaterial	material,
		bool	sagCorrection = true )

calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.

Definition at line 354 of file TRunge.cxx.

                                                                                    {
  float       tof;
  HepVector3D p;
  return ( approach( l, tof, p, material, doSagCorrection ) );
}

Referenced by approach().

approach() [2/6]

int TRunge::approach	(	TMLink &	,
		const RkFitMaterial	,
		bool	sagCorrection = true )

calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.

approach() [3/6]

int TRunge::approach	(	TMLink &	,
		const RkFitMaterial	,
		bool	sagCorrection = true )

calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.

approach() [4/6]

int TRunge::approach	(	TMLink &	l,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material,
		bool	sagCorrection = true )

Definition at line 360 of file TRunge.cxx.

                                             {
  HepPoint3D      xw;
  HepPoint3D      wireBackwardPosition;
  HepVector3D     v;
  HepPoint3D      onWire, onTrack;
  double          onWire_x, onWire_y, onWire_z, zwf, zwb;
  const TMDCWire& w    = *l.wire();
  xw                   = w.xyPosition();
  wireBackwardPosition = w.backwardPosition();
  v                    = w.direction();
 
  unsigned stepNum = 0;
  if ( approach_line( l, wireBackwardPosition, v, onWire, onTrack, tof, p, material,
                      stepNum ) < 0 )
  {
    // cout<<"TR::error approach_line"<<endl;
    return ( -1 );
  }
  zwf = w.forwardPosition().z();
  zwb = w.backwardPosition().z();
  // protection for strange wire hit info. (Only 'onWire' is corrected)
  if ( onWire.z() > zwf ) w.wirePosition( zwf, onWire, wireBackwardPosition, (HepVector3D&)v );
  else if ( onWire.z() < zwb )
    w.wirePosition( zwb, onWire, wireBackwardPosition, (HepVector3D&)v );
 
  // onWire,onTrack filled
 
  if ( !doSagCorrection )
  {
    l.positionOnWire( onWire );
    l.positionOnTrack( onTrack );
    double phipoint = atan( ( onTrack.y() - _pivot.y() ) / onTrack.x() - _pivot.x() );
    //    cout<<"dphi track:"<<l.dPhi()<<endl;
    //   cout<<"dphi rk:"<<phipoint-_a[0]<<endl;
    //    l.dPhi(phipoint-_a[0]);
    return ( 0 ); // no sag correction
  }
  // Sag correction
  //   loop for sag correction
  onWire_y = onWire.y();
  onWire_z = onWire.z();
 
  unsigned nTrial = 1;
  while ( nTrial < 100 )
  {
    // cout<<"TR: nTrial "<<nTrial<<endl;
    w.wirePosition( onWire_z, xw, wireBackwardPosition, (HepVector3D&)v );
    if ( approach_line( l, wireBackwardPosition, v, onWire, onTrack, tof, p, material,
                        stepNum ) < 0 )
      return ( -1 );
    if ( fabs( onWire_y - onWire.y() ) < 0.0001 ) break; // |dy|< 1 micron
    onWire_y = onWire.y();
    onWire_z += ( onWire.z() - onWire_z ) / 2;
    // protection for strange wire hit info.
    if ( onWire_z > zwf ) onWire_z = zwf;
    else if ( onWire_z < zwb ) onWire_z = zwb;
 
    nTrial++;
  }
  //  cout<<"TR  nTrial="<<nTrial<<endl;
  l.positionOnWire( onWire );
  l.positionOnTrack( onTrack );
  return ( nTrial );
}

approach() [5/6]

int TRunge::approach	(	TMLink &	,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	,
		bool	sagCorrection = true )

approach() [6/6]

int TRunge::approach	(	TMLink &	,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	,
		bool	sagCorrection = true )

approach_line() [1/9]

int TRunge::approach_line	(	TMLink &	l,
		const HepPoint3D &	w0,
		const HepVector3D &	v,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest points between a line and this track

Definition at line 426 of file TRunge.cxx.

                                                          {
  float       tof;
  HepVector3D p;
  return ( approach_line( l, w0, v, onLine, onTrack, tof, p, material ) );
}

Referenced by approach(), approach_line(), and approach_line().

approach_line() [2/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest points between a line and this track

approach_line() [3/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest points between a line and this track

approach_line() [4/9]

int TRunge::approach_line	(	TMLink &	l,
		const HepPoint3D &	w0,
		const HepVector3D &	v,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material )

Definition at line 434 of file TRunge.cxx.

                                                          {
  unsigned stepNum = 0;
  return ( approach_line( l, w0, v, onLine, onTrack, tof, p, material, stepNum ) );
}

approach_line() [5/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material )

approach_line() [6/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material )

approach_line() [7/9]

int TRunge::approach_line	(	TMLink &	l,
		const HepPoint3D &	w0,
		const HepVector3D &	v,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material,
		unsigned &	stepNum )

Definition at line 441 of file TRunge.cxx.

                                                                             {
  //  line = [w0] + t * [v]    -> [onLine]
  if ( _Nstep == 0 )
  {
    if ( _stepSize == 0 ) Fly_SC();
    else Fly( material );
  }
  // cout<<"TR::approach_line stepNum="<<stepNum<<endl;
 
  const double w0x = w0.x();
  const double w0y = w0.y();
  const double w0z = w0.z();
  // cout<<"TR::line w0="<<w0x<<","<<w0y<<","<<w0z<<endl;
  const double vx  = v.x();
  const double vy  = v.y();
  const double vz  = v.z();
  const double v_2 = vx * vx + vy * vy + vz * vz;
 
  const float clight = 29.9792458; //[cm/ns]
  // const float M2=_mass*_mass;
  // const float& M2=_mass2;
  const float p2         = _y[0][3] * _y[0][3] + _y[0][4] * _y[0][4] + _y[0][5] * _y[0][5];
  const float tof_factor = 1. / clight * sqrt( 1 + _mass2 / p2 );
 
  // search for the closest point in cache   (point - line)
  //  unsigned stepNum;
  double l2_old = DBL_MAX;
  if ( stepNum > _Nstep ) stepNum = 0;
  unsigned stepNumLo;
  unsigned stepNumHi;
  if ( stepNum == 0 && _stepSize != 0 )
  { // skip
    const double dx = _y[0][0] - w0x;
    const double dy = _y[0][1] - w0y;
    stepNum = (unsigned)( sqrt( ( dx * dx + dy * dy ) * ( 1 + _a[4] * _a[4] ) ) / _stepSize );
  }
  unsigned mergin;
  if ( _stepSize == 0 )
  {
    mergin = 10; // 10 step back
  }
  else
  {
    mergin = (unsigned)( 1.0 / _stepSize ); // 1mm back
  }
  if ( stepNum > mergin ) stepNum -= mergin;
  else stepNum = 0;
  if ( stepNum >= _Nstep ) stepNum = _Nstep - 1;
  // hunt
  //  unsigned inc=1;
  unsigned inc = ( mergin >> 1 ) + 1;
  stepNumLo    = stepNum;
  stepNumHi    = stepNum;
  for ( ;; )
  {
    const double dx = _y[stepNumHi][0] - w0x;
    const double dy = _y[stepNumHi][1] - w0y;
    const double dz = _y[stepNumHi][2] - w0z;
    const double t  = ( dx * vx + dy * vy + dz * vz ) / v_2;
    //    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zvec;
    zvec.setX( 0 );
    zvec.setY( 0 );
    zvec.setZ( 1 );
    HepVector3D hitv = t * v;
    HepPoint3D  pp, onw;
    pp.setX( _y[stepNumHi][0] );
    pp.setY( _y[stepNumHi][1] );
    pp.setZ( _y[stepNumHi][2] );
    double       zwire = hitv.dot( zvec ) + w0z;
    double       dtl   = DisToWire( l, zwire, pp, onw );
    const double l2    = dtl * dtl;
    if ( l2 > l2_old ) break;
    l2_old    = l2;
    stepNumLo = stepNumHi;
    //    inc+=inc;
    stepNumHi += inc;
    if ( stepNumHi >= _Nstep )
    {
      stepNumHi = _Nstep;
      break;
    }
  }
  // locate (2-bun-hou, bisection method)
  while ( stepNumHi - stepNumLo > 1 )
  {
    unsigned     j  = ( stepNumHi + stepNumLo ) >> 1;
    const double dx = _y[j][0] - w0x;
    const double dy = _y[j][1] - w0y;
    const double dz = _y[j][2] - w0z;
    const double t  = ( dx * vx + dy * vy + dz * vz ) / v_2;
    //    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zv;
    zv.setX( 0 );
    zv.setY( 0 );
    zv.setZ( 1 );
    HepVector3D hitv = t * v;
    HepPoint3D  point, onwir;
    point.setX( _y[j][0] );
    point.setY( _y[j][1] );
    point.setZ( _y[j][2] );
    double       zwir = hitv.dot( zv ) + w0z;
    double       dtll = DisToWire( l, zwir, point, onwir );
    const double l2   = dtll * dtll;
    if ( l2 > l2_old ) { stepNumHi = j; }
    else
    {
      l2_old    = l2;
      stepNumLo = j;
    }
  }
  // stepNum=stepNumHi;
  stepNum = stepNumLo;
  /*
  for(;stepNum<_Nstep;stepNum++){
    const double dx = _y[stepNum][0]-w0x;
    const double dy = _y[stepNum][1]-w0y;
    const double dz = _y[stepNum][2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    if(l2 > l2_old) break;
    l2_old=l2;
    //    const float p2 = _y[stepNum][3]*_y[stepNum][3]+
    //                      _y[stepNum][4]*_y[stepNum][4]+
    //                      _y[stepNum][5]*_y[stepNum][5];
    //    tof+=_stepSize/clight*sqrt(1+M2/p2);
  }
  */
  //  cout<<"TR  stepNum="<<stepNum<<endl;
  // if(stepNum>=_Nstep) return(-1);    // not found
  // stepNum--;
  if ( _stepSize == 0 )
  {
    double dstep = 0;
    for ( unsigned i = 0; i < stepNum; i++ ) dstep += _h[i];
    tof = dstep * tof_factor;
  }
  else { tof = stepNum * _stepSize * tof_factor; }
 
  // propagate the track and search for the closest point
  // 2-bun-hou (bisection method)
  /*
  double y[6],y_old[6];
  for(unsigned i=0;i<6;i++) y[i]=_y[stepNum][i];
  double step=_stepSize;
  double step2=0;
  for(;;){
    for(unsigned i=0;i<6;i++) y_old[i]=y[i];
    Propagate(y,step);
    const double dx = y[0]-w0x;
    const double dy = y[1]-w0y;
    const double dz = y[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    if(l2 > l2_old){  // back
      for(unsigned i=0;i<6;i++) y[i]=y_old[i];
    }else{        // propagate
      l2_old=l2;
      //      const float p2=y[3]*y[3]+y[4]*y[4]+y[5]*y[5];
      //      tof+=step/clight*sqrt(1+M2/p2);
      step2+=step;
    }
    step/=2;
    if(step < 0.0001) break;    // step < 1 [um]
  }
  */
  // Hasamiuchi-Hou (false position method)
  /*
  double y[6],y1[6],y2[6];
  for(unsigned i=0;i<6;i++) y1[i]=_y[stepNum][i];
  for(unsigned i=0;i<6;i++) y2[i]=_y[stepNum+2][i];
  double minStep=0;
  double maxStep=_stepSize*2;
  double step2=0;
  const double A[3][3]={{1-vx*vx/v_2, -vx*vy/v_2, -vx*vz/v_2},
                        { -vy*vx/v_2,1-vy*vy/v_2, -vy*vz/v_2},
                        { -vz*vx/v_2, -vz*vy/v_2,1-vz*vz/v_2}};
 
  double Y1=0;
  {
    const double dx = y1[0]-w0x;
    const double dy = y1[1]-w0y;
    const double dz = y1[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
    const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
    const double pmag=sqrt(y1[3]*y1[3]+y1[4]*y1[4]+y1[5]*y1[5]);
    for(int j=0;j<3;j++){
      double g=0;
      for(int k=0;k<3;k++) g+=A[j][k]*d[k];
      Y1+=y1[j+3]*g;
    }
    Y1=Y1/pmag/l;
  }
  double Y2=0;
  {
    const double dx = y2[0]-w0x;
    const double dy = y2[1]-w0y;
    const double dz = y2[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
    const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
    const double pmag=sqrt(y2[3]*y2[3]+y2[4]*y2[4]+y2[5]*y2[5]);
    for(int j=0;j<3;j++){
      double g=0;
      for(int k=0;k<3;k++) g+=A[j][k]*d[k];
      Y2+=y2[j+3]*g;
    }
    Y2=Y2/pmag/l;
  }
  if(Y1*Y2>=0) cout<<"TR  fatal error!"<<endl;
  double step_old= DBL_MAX;
  for(;;){
    const double step=(minStep*Y2-maxStep*Y1)/(Y2-Y1);
    if(fabs(step-step_old)<0.0001) break;
    step_old=step;
    for(unsigned i=0;i<6;i++) y[i]=y1[i];
    Propagate(y,step-minStep);
    double Y=0;
    {
      const double dx = y[0]-w0x;
      const double dy = y[1]-w0y;
      const double dz = y[2]-w0z;
      const double t = (dx*vx+dy*vy+dz*vz)/v_2;
      const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
      const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
      const double pmag=sqrt(y[3]*y[3]+y[4]*y[4]+y[5]*y[5]);
      for(int j=0;j<3;j++){
        double g=0;
        for(int k=0;k<3;k++) g+=A[j][k]*d[k];
        Y+=y[j+3]*g;
      }
      Y=Y/pmag/l;
    }
    if(Y1*Y<0){ //Y->Y2
      Y2=Y;
      for(unsigned i=0;i<6;i++) y2[i]=y[i];
      maxStep=step;
    }else{  //Y->Y1
      Y1=Y;
      for(unsigned i=0;i<6;i++) y1[i]=y[i];
      minStep=step;
    }
    if(Y1==Y2) break;
  }
  step2=step_old;
  */
  // Newton method
  double y[6];
  for ( unsigned i = 0; i < 6; i++ ) y[i] = _y[stepNum][i];
  // memcpy(y,_y[stepNum],sizeof(double)*6);
  double       step2   = 0;
  const double A[3][3] = { { 1 - vx * vx / v_2, -vx * vy / v_2, -vx * vz / v_2 },
                           { -vy * vx / v_2, 1 - vy * vy / v_2, -vy * vz / v_2 },
                           { -vz * vx / v_2, -vz * vy / v_2, 1 - vz * vz / v_2 } };
  double       factor  = 1;
  double       g[3];
  double       f[6];
  double       Af[3], Af2[3];
  unsigned     i, j, k;
  HepPoint3D   point, onwir;
  for ( i = 0; i < 10; i++ )
  {
    // cout<<"TR::line "<<y[0]<<","<<y[1]<<","<<y[2]<<","<<y[3]<<","<<y[4]<<","<<y[5]<<endl;
    const double dx   = y[0] - w0x;
    const double dy   = y[1] - w0y;
    const double dz   = y[2] - w0z;
    const double t    = ( dx * vx + dy * vy + dz * vz ) / v_2;
    const double d[3] = { dx - t * vx, dy - t * vy, dz - t * vz };
    // const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zv;
    zv.setX( 0 );
    zv.setY( 0 );
    zv.setZ( 1 );
    HepVector3D hitv = t * v;
    point.setX( y[0] );
    point.setY( y[1] );
    point.setZ( y[2] );
    double       zwir = hitv.dot( zv ) + w0z;
    double       dtll = DisToWire( l, zwir, point, onwir );
    const double l2   = dtll * dtll;
    for ( j = 0; j < 3; j++ )
    {
      g[j] = 0;
      for ( k = 0; k < 3; k++ ) g[j] += A[j][k] * d[k];
    }
    // cout<<"g="<<g[0]<<","<<g[1]<<","<<g[2]<<endl;
    Function( y, f );
    // cout<<"f="<<f[0]<<","<<f[1]<<","<<f[2]<<","<<f[3]<<","<<f[4]<<","<<f[5]<<endl;
    // cout<<"A="<<A[0][0]<<","<<A[0][1]<<","<<A[0][2]<<endl;
    // cout<<"A="<<A[1][0]<<","<<A[1][1]<<","<<A[1][2]<<endl;
    // cout<<"A="<<A[2][0]<<","<<A[2][1]<<","<<A[2][2]<<endl;
    double Y = 0;
    for ( j = 0; j < 3; j++ ) Y += y[j + 3] * g[j];
    double dYds = 0;
    for ( j = 0; j < 3; j++ ) dYds += f[j + 3] * g[j];
    // cout<<"dYds="<<dYds<<endl;
    for ( j = 0; j < 3; j++ )
    {
      Af[j]  = 0;
      Af2[j] = 0;
    }
    for ( j = 0; j < 3; j++ )
    {
      // Af[j]=0;
      // Af2[j]=0;
      for ( k = 0; k < 3; k++ ) Af[j] += ( A[j][k] * f[k] );
      for ( k = 0; k < 3; k++ ) Af2[j] += ( A[j][k] * Af[k] );
      dYds += y[j + 3] * Af2[j];
      // cout<<j<<" dYds="<<dYds<<endl;
    }
    // cout<<"dYds="<<dYds<<endl;
    const double step = -Y / dYds * factor;
    // cout<<"TR  step="<<step<<" i="<<i<<endl;
    if ( fabs( step ) < 0.00001 ) break; // step < 1 [um]
    if ( l2 > l2_old ) factor /= 2;
    l2_old = l2;
    Propagate( y, step, material );
    step2 += step;
  }
 
  tof += step2 * tof_factor;
 
  onTrack.setX( y[0] );
  onTrack.setY( y[1] );
  onTrack.setZ( y[2] );
  p.setX( y[3] );
  p.setY( y[4] );
  p.setZ( y[5] );
 
  const double dx = y[0] - w0x;
  const double dy = y[1] - w0y;
  const double dz = y[2] - w0z;
  const double s  = ( dx * vx + dy * vy + dz * vz ) / v_2;
  // onLine.setX(w0x+s*vx);
  // onLine.setY(w0y+s*vy);
  // onLine.setZ(w0z+s*vz);
  onLine.setX( onwir.x() );
  onLine.setY( onwir.y() );
  onLine.setZ( onwir.z() );
  return ( 0 );
}

approach_line() [8/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material,
		unsigned &	stepNum )

approach_line() [9/9]

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material,
		unsigned &	stepNum )

approach_point() [1/3]

int TRunge::approach_point	(	TMLink &	l,
		const HepPoint3D &	p0,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest point between a point and this track

Definition at line 786 of file TRunge.cxx.

                                                           {
  if ( _Nstep == 0 )
  {
    if ( _stepSize == 0 ) Fly_SC();
    else Fly( material );
  }
 
  double x0 = p0.x();
  double y0 = p0.y();
  double z0 = p0.z();
 
  // tof=0;
  // const float clight=29.9792458; //[cm/ns]
  // const float M2=_mass*_mass;
 
  // search for the closest point in cache
  unsigned stepNum;
  double   l2_old = DBL_MAX;
  for ( stepNum = 0; stepNum < _Nstep; stepNum++ )
  {
    double l2 = ( _y[stepNum][0] - x0 ) * ( _y[stepNum][0] - x0 ) +
                ( _y[stepNum][1] - y0 ) * ( _y[stepNum][1] - y0 ) +
                ( _y[stepNum][2] - z0 ) * ( _y[stepNum][2] - z0 );
    if ( l2 > l2_old ) break;
    l2_old = l2;
    // const double p2 = _y[stepNum][3]*_y[stepNum][3]+
    //                   _y[stepNum][4]*_y[stepNum][4]+
    //                   _y[stepNum][5]*_y[stepNum][5];
    // tof+=_stepSize/clight*sqrt(1+M2/p2);
  }
  if ( stepNum >= _Nstep ) return ( -1 ); // not found
  stepNum--;
 
  // propagate the track and search for the closest point
  double y[6], y_old[6];
  for ( unsigned i = 0; i < 6; i++ ) y[i] = _y[stepNum][i];
  double step = _stepSize;
  for ( ;; )
  {
    for ( unsigned i = 0; i < 6; i++ ) y_old[i] = y[i];
    Propagate( y, step, material );
    double l2 = ( y[0] - x0 ) * ( y[0] - x0 ) + ( y[1] - y0 ) * ( y[1] - y0 ) +
                ( y[2] - z0 ) * ( y[2] - z0 );
    if ( l2 > l2_old )
    { // back
      for ( unsigned i = 0; i < 6; i++ ) y[i] = y_old[i];
    }
    else
    { // propagate
      l2_old = l2;
      // const double p2=y[3]*y[3]+y[4]*y[4]+y[5]*y[5];
      // tof+=step/clight*sqrt(1+M2/p2);
    }
    step /= 2;
    if ( step < 0.0001 ) break; // step < 1 [um]
  }
 
  onTrack.setX( y[0] );
  onTrack.setY( y[1] );
  onTrack.setZ( y[2] );
  // p.setX(y[3]);
  // p.setY(y[4]);
  // p.setZ(y[5]);
  return ( 0 );
}

approach_point() [2/3]

int TRunge::approach_point	(	TMLink &	,
		const HepPoint3D &	,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest point between a point and this track

approach_point() [3/3]

int TRunge::approach_point	(	TMLink &	,
		const HepPoint3D &	,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material )

caluculate closest point between a point and this track

BfieldID() [1/6]

int TRunge::BfieldID

(

int

id

)

set B field map ID

Definition at line 307 of file TRunge.cxx.

                             {
  _bfieldID = id;
  _bfield   = Bfield::getBfield( _bfieldID );
  _Nstep    = 0;
  return ( _bfieldID );
}

BfieldID() [2/6]

int TRunge::BfieldID

(

int

)

set B field map ID

BfieldID() [3/6]

int TRunge::BfieldID

(

int

)

set B field map ID

BfieldID() [4/6]

int TRunge::BfieldID

(

void

)

const

returns B field ID

Definition at line 268 of file TRunge.cxx.

{ return ( _bfieldID ); }

Referenced by TRunge().

BfieldID() [5/6]

int TRunge::BfieldID

(

void

)

const

returns B field ID

BfieldID() [6/6]

int TRunge::BfieldID

(

void

)

const

returns B field ID

chi2() [1/3]

double TRunge::chi2

(

void

)

const

returns chi2.

Definition at line 257 of file TRunge.cxx.

{ return _chi2; }

Referenced by TRunge(), and TRunge().

chi2() [2/3]

double TRunge::chi2

(

void

)

const

returns chi2.

chi2() [3/3]

double TRunge::chi2

(

void

)

const

returns chi2.

dEpath() [1/3]

double TRunge::dEpath	(	double	mass,
		double	path,
		double	p ) const

Definition at line 1265 of file TRunge.cxx.

                                                                {
  double z       = 4.72234;
  double a       = 9.29212;
  double Density = 0.00085144;
  double coeff   = Density * z / a;
  double i       = 51.4709;
  double isq     = i * i * 1e-18;
  //       double sigma0_2 = 0.1569*coeff*path;
  //
  //  if (sigma0_2<0) return 0;
  //
  //  double psq = p * p;
  //  double bsq = psq / (psq + mass * mass);
  //
  //  // Correction for relativistic particles :
  //  double sigma_2 = sigma0_2*(1-0.5*bsq)/(1-bsq);
  //
  //  if (sigma_2<0) return 0;
  //
  //  // Return sigma in GeV !!
  //  return sqrt(sigma_2)*0.001;
  const double Me  = 0.000510999;
  double       psq = p * p;
  double       bsq = psq / ( psq + mass * mass );
  double       esq = psq / ( mass * mass );
 
  double s = Me / mass;
  double w = ( 4 * Me * esq / ( 1 + 2 * s * sqrt( 1 + esq ) + s * s ) );
 
  // Density correction :
  double cc, x0;
  cc = 1 + 2 * log( sqrt( isq ) / ( 28.8E-09 * sqrt( coeff ) ) );
  if ( cc < 5.215 ) x0 = 0.2;
  else x0 = 0.326 * cc - 1.5;
  double x1( 3 ), xa( cc / 4.606 ), aa;
  aa = 4.606 * ( xa - x0 ) / ( ( x1 - x0 ) * ( x1 - x0 ) * ( x1 - x0 ) );
  double delta( 0 );
  double x( log10( sqrt( esq ) ) );
  if ( x > x0 )
  {
    delta = 4.606 * x - cc;
    if ( x < x1 ) delta = delta + aa * ( x1 - x ) * ( x1 - x ) * ( x1 - x );
  }
 
  // Shell correction :
  float f1, f2, f3, f4, f5, ce;
  f1 = 1 / esq;
  f2 = f1 * f1;
  f3 = f1 * f2;
  f4 = ( f1 * 0.42237 + f2 * 0.0304 - f3 * 0.00038 ) * 1E12;
  f5 = ( f1 * 3.858 - f2 * 0.1668 + f3 * 0.00158 ) * 1E18;
  ce = f4 * isq + f5 * isq * sqrt( isq );
 
  return ( 0.0001535 * coeff / bsq *
           ( log( Me * esq * w / isq ) - 2 * bsq - delta - 2.0 * ce / z ) ) *
         path;
}

dEpath() [2/3]

double TRunge::dEpath	(	double	,
		double	,
		double	) const

dEpath() [3/3]

double TRunge::dEpath	(	double	,
		double	,
		double	) const

DisToWire() [1/3]

double TRunge::DisToWire	(	TMLink &	l,
		double	z,
		HepPoint3D	a,
		HepPoint3D &	b )

Definition at line 1236 of file TRunge.cxx.

                                          { // by liucy 2011/7/28
  double          zinter[400];
  double          ddist  = 999;
  double          finalz = 0;
  const TMDCWire& w      = *l.wire();
  HepPoint3D      onwire;
  //   for(int i=0;i<20;i++){
  //   zinter[i]=z-0.0002*i+0.002;
  HepPoint3D point = w.xyPosition( z );
  //     HepPoint3D point=w.xyPosition(zinter[i]);
  onwire.setX( point.x() );
  onwire.setY( point.y() );
  //  onwire.setZ(zinter[i]);
  onwire.setZ( z );
  //     double
  //     dist=sqrt((point.x()-a.x())*(point.x()-a.x())+(point.y()-a.y())*(point.y()-a.y())+(zinter[i]-a.z())*(zinter[i]-a.z()));
  double dist =
      sqrt( ( point.x() - a.x() ) * ( point.x() - a.x() ) +
            ( point.y() - a.y() ) * ( point.y() - a.y() ) + ( z - a.z() ) * ( z - a.z() ) );
  if ( dist < ddist )
  {
    ddist = dist;
    //       finalz=zinter[i];
    b = onwire;
  }
  //   }
  return ddist;
}

Referenced by approach_line().

DisToWire() [2/3]

double TRunge::DisToWire	(	TMLink &	,
		double	,
		HepPoint3D	,
		HepPoint3D &	)

DisToWire() [3/3]

double TRunge::DisToWire	(	TMLink &	,
		double	,
		HepPoint3D	,
		HepPoint3D &	)

dr() [1/3]

double TRunge::dr

(

void

)

const

Track parameters (at pivot).

Definition at line 237 of file TRunge.cxx.

{ return ( _a[0] ); }

dr() [2/3]

double TRunge::dr

(

void

)

const

Track parameters (at pivot).

dr() [3/3]

double TRunge::dr

(

void

)

const

Track parameters (at pivot).

dz() [1/3]

double TRunge::dz

(

void

)

const

Definition at line 243 of file TRunge.cxx.

{ return ( _a[3] ); }

Referenced by approach_line(), and intersect_xy_plane().

dz() [2/3]

double TRunge::dz

(

void

)

const

dz() [3/3]

double TRunge::dz

(

void

)

const

Ea() [1/6]

const SymMatrix & TRunge::Ea

(

const SymMatrix &

tEa

)

set helix error matrix

Definition at line 301 of file TRunge.cxx.

                                                  {
  _Ea    = tEa;
  _Nstep = 0;
  return ( _Ea );
}

Ea() [2/6]

const SymMatrix & TRunge::Ea

(

const SymMatrix &

)

set helix error matrix

Ea() [3/6]

const SymMatrix & TRunge::Ea

(

const SymMatrix &

)

set helix error matrix

Ea() [4/6]

const SymMatrix & TRunge::Ea

(

void

)

const

returns error matrix

Definition at line 251 of file TRunge.cxx.

{ return ( _Ea ); }

Referenced by TRunge(), TRunge(), and TRunge().

Ea() [5/6]

const SymMatrix & TRunge::Ea

(

void

)

const

returns error matrix

Ea() [6/6]

const SymMatrix & TRunge::Ea

(

void

)

const

returns error matrix

eloss() [1/3]

void TRunge::eloss	(	double	path,
		const RkFitMaterial *	material,
		double	mass,
		double	y[6],
		int	index ) const

Definition at line 1322 of file TRunge.cxx.

                                      {
  double psq = y[5] * y[5] + y[3] * y[3] + y[4] * y[4];
 
  double p = sqrt( psq );
 
  double dE = materiral->dE( mass, path, p );
  if ( index > 0 ) psq += dE * ( dE + 2 * sqrt( mass * mass + psq ) );
  else psq += dE * ( dE - 2 * sqrt( mass * mass + psq ) );
  double pnew  = sqrt( psq );
  double coeff = pnew / p;
 
  y[3] = y[3] * coeff;
  y[4] = y[4] * coeff;
  y[5] = y[5] * coeff;
}

Referenced by RkFitCylinder::updateTrack().

eloss() [2/3]

void TRunge::eloss	(	double	path,
		const RkFitMaterial *	material,
		double	mass,
		double	y[6],
		int	index ) const

eloss() [3/3]

void TRunge::eloss	(	double	path,
		const RkFitMaterial *	material,
		double	mass,
		double	y[6],
		int	index ) const

Eps() [1/6]

double TRunge::Eps

(

double

eps

)

Definition at line 325 of file TRunge.cxx.

                               {
  _eps   = eps;
  _Nstep = 0;
  return ( _eps );
}

Eps() [2/6]

double TRunge::Eps

(

double

)

Eps() [3/6]

double TRunge::Eps

(

double

)

Eps() [4/6]

double TRunge::Eps

(

void

)

const

Definition at line 273 of file TRunge.cxx.

{ return ( _eps ); }

Referenced by TRunge().

Eps() [5/6]

double TRunge::Eps

(

void

)

const

Eps() [6/6]

double TRunge::Eps

(

void

)

const

Fly() [1/3]

unsigned TRunge::Fly

(

const RkFitMaterial

material

)

make the trajectory in cache, return the number of step

Definition at line 853 of file TRunge.cxx.

                                                   {
  double   y[6];
  unsigned Nstep;
  double   flightlength;
 
  flightlength = 0;
  SetFirst( y );
  for ( Nstep = 0; Nstep < TRunge_MAXstep; Nstep++ )
  {
    for ( unsigned j = 0; j < 6; j++ ) _y[Nstep][j] = y[j];
    // memcpy(_y[Nstep],y,sizeof(double)*6);
 
    Propagate( y, _stepSize, material );
 
    if ( y[2] > 160 || y[2] < -80 || ( y[0] * y[0] + y[1] * y[1] ) > 8100 ) break;
    // if position is out side of MDC, stop to fly
    //   R>90cm, z<-80cm,160cm<z
 
    flightlength += _stepSize;
    if ( flightlength > _maxflightlength ) break;
 
    _h[Nstep] = _stepSize;
  }
  _Nstep = Nstep + 1;
  return ( _Nstep );
}

Referenced by approach_line(), and approach_point().

Fly() [2/3]

unsigned TRunge::Fly

(

const RkFitMaterial

material

)

make the trajectory in cache, return the number of step

Fly() [3/3]

unsigned TRunge::Fly

(

const RkFitMaterial

material

)

make the trajectory in cache, return the number of step

Fly_SC() [1/3]

unsigned TRunge::Fly_SC

(

void

)

Definition at line 1114 of file TRunge.cxx.

                              { // fly the particle track with stepsize control
  double   y[6], dydx[6], yscal[6];
  unsigned Nstep;
  double   step, stepdid, stepnext;
  unsigned i;
  double   flightlength;
 
  // yscal[0]=0.1; //1mm  -> error = 1mm*EPS
  // yscal[1]=0.1; //1mm
  // yscal[2]=0.1; //1mm
  // yscal[3]=0.001; //1MeV
  // yscal[4]=0.001; //1MeV
  // yscal[5]=0.001; //1MeV
 
  step         = default_stepSize0;
  flightlength = 0;
  SetFirst( y );
  for ( Nstep = 0; Nstep < TRunge_MAXstep; Nstep++ )
  {
    for ( i = 0; i < 6; i++ ) _y[Nstep][i] = y[i]; // save each step
    // memcpy(_y[Nstep],y,sizeof(double)*6);
 
    Function( y, dydx );
    // for(i=0;i<6;i++) yscal[i]=abs(y[i])+abs(dydx[i]*step)+TINY;
 
    Propagate_QC( y, dydx, step, _eps, _yscal, stepdid, stepnext );
 
    if ( y[2] > 160 || y[2] < -80 || ( y[0] * y[0] + y[1] * y[1] ) > 8100 ) break;
    // if position is out side of MDC, stop to fly
    //   R>90cm, z<-80cm,160cm<z
 
    flightlength += stepdid;
    if ( flightlength > _maxflightlength ) break;
 
    _h[Nstep] = stepdid;
    // cout<<"TR:"<<Nstep<<":step="<<stepdid<<endl;
    if ( stepnext < _stepSizeMin ) step = _stepSizeMin;
    else if ( stepnext > _stepSizeMax ) step = _stepSizeMax;
    else step = stepnext;
  }
  _Nstep = Nstep + 1;
  // cout<<"TR: Nstep="<<_Nstep<<endl;
  return ( _Nstep );
}

Referenced by approach_line(), and approach_point().

Fly_SC() [2/3]

unsigned TRunge::Fly_SC

(

void

)

Fly_SC() [3/3]

unsigned TRunge::Fly_SC

(

void

)

Function() [1/3]

void TRunge::Function	(	const double	y[6],
		double	f[6] )

Definition at line 934 of file TRunge.cxx.

                                                      {
  // return the value of formula
  //  y[6] = (x,y,z,px,py,pz)
  //  f[6] = ( dx[3]/ds, dp[3]/ds )
  //  dx/ds = p/|p|
  //  dp/ds = e/|e| 1/alpha (p x B)/|p||B|
  //    alpha = 1/cB = 333.5640952/B[Tesla]  [cm/(GeV/c)]
  //  const float Bx = _bfield->bx(y[0],y[1],y[2]);
  //  const float By = _bfield->by(y[0],y[1],y[2]);
  //  const float Bz = _bfield->bz(y[0],y[1],y[2]);  //[kGauss]
  double      pmag;
  double      factor;
  HepVector3D B;
  HepPoint3D  pos;
  pos.setX( (float)y[0] );
  pos.setY( (float)y[1] );
  pos.setZ( (float)y[2] );
  // cout<<"TR::pos="<<pos[0]<<","<<pos[1]<<","<<pos[2]<<endl;
  //  _bfield->fieldMap(pos,B);
  // cout<<"TR::B="<<B[0]<<","<<B[1]<<","<<B[2]<<endl;
  //  const double Bmag = sqrt(Bx*Bx+By*By+Bz*Bz);
  getPmgnIMF()->fieldVector( 10 * pos, B );
  pmag = sqrt( y[3] * y[3] + y[4] * y[4] + y[5] * y[5] );
  f[0] = y[3] / pmag; // p/|p|
  f[1] = y[4] / pmag;
  f[2] = y[5] / pmag;
 
  //  const factor = _charge/(alpha2/Bmag)/pmag/Bmag;
  // factor = ((double)_charge)/alpha2/10.; //[(GeV/c)/(cm kG)]
  // factor = ((double)_charge)/(333.564095/(-1000*(getPmgnIMF()->getReferField())))/10.;
  // //[(GeV/c)/(cm kG)]
  double Bmag = sqrt( B.x() * B.x() + B.y() * B.y() + B.z() * B.z() );
  factor      = ( (double)_charge ) / ( 0.333564095 );
  //  f[3] = factor*(f[1]*Bz-f[2]*By);
  //  f[4] = factor*(f[2]*Bx-f[0]*Bz);
  //  f[5] = factor*(f[0]*By-f[1]*Bx);
  f[3] = factor * ( f[1] * B.z() - f[2] * B.y() );
  f[4] = factor * ( f[2] * B.x() - f[0] * B.z() );
  f[5] = factor * ( f[0] * B.y() - f[1] * B.x() );
}

Referenced by approach_line(), Fly_SC(), Propagate(), Propagate1(), and Propagate_QC().

Function() [2/3]

void TRunge::Function	(	const double	y[6],
		double	f[6] )

Function() [3/3]

void TRunge::Function	(	const double	y[6],
		double	f[6] )

GetStep() [1/3]

int TRunge::GetStep	(	unsigned	stepNum,
		double &	step ) const

Definition at line 1000 of file TRunge.cxx.

                                                          {
  if ( stepNum >= _Nstep || stepNum >= TRunge_MAXstep ) return ( -1 );
  step = _h[stepNum];
  return ( 0 );
}

GetStep() [2/3]

int TRunge::GetStep	(	unsigned	stepNum,
		double &	step ) const

GetStep() [3/3]

int TRunge::GetStep	(	unsigned	stepNum,
		double &	step ) const

GetXP() [1/3]

int TRunge::GetXP	(	unsigned	stepNum,
		double	y[6] ) const

Definition at line 994 of file TRunge.cxx.

                                                       {
  if ( stepNum >= _Nstep || stepNum >= TRunge_MAXstep ) return ( -1 );
 
  for ( unsigned i = 0; i < 6; i++ ) y[i] = _y[stepNum][i];
  return ( 0 );
}

GetXP() [2/3]

int TRunge::GetXP	(	unsigned	stepNum,
		double	y[6] ) const

GetXP() [3/3]

int TRunge::GetXP	(	unsigned	stepNum,
		double	y[6] ) const

helix() [1/3]

Helix TRunge::helix

(

void

)

const

returns helix class

Definition at line 253 of file TRunge.cxx.

{ return ( Helix( _pivot, _a, _Ea ) ); }

Referenced by RkFitCylinder::intersect(), RkFitCylinder::intersect(), intersect_cylinder(), intersect_xy_plane(), intersect_yz_plane(), intersect_zx_plane(), pivot(), SetFlightLength(), TRunge(), and TRunge().

helix() [2/3]

Helix TRunge::helix

(

void

)

const

returns helix class

helix() [3/3]

Helix TRunge::helix

(

void

)

const

returns helix class

intersect_cylinder() [1/3]

double TRunge::intersect_cylinder

(

double

r

)

const

Definition at line 1338 of file TRunge.cxx.

                                                  {
  double m_rad  = helix().radius();
  double l      = helix().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 = helix().center().phi() - acos( cosPhi ) - helix().phi0();
  if ( dPhi < -M_PI ) dPhi += 2 * M_PI;
 
  return dPhi;
}

Referenced by RkFitCylinder::intersect(), RkFitCylinder::intersect(), intersect_yz_plane(), and intersect_zx_plane().

intersect_cylinder() [2/3]

double TRunge::intersect_cylinder

(

double

r

)

const

intersect_cylinder() [3/3]

double TRunge::intersect_cylinder

(

double

r

)

const

intersect_xy_plane() [1/3]

double TRunge::intersect_xy_plane

(

double

z

)

const

Definition at line 1380 of file TRunge.cxx.

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

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

intersect_xy_plane() [2/3]

double TRunge::intersect_xy_plane

(

double

z

)

const

intersect_xy_plane() [3/3]

double TRunge::intersect_xy_plane

(

double

z

)

const

intersect_yz_plane() [1/3]

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

Definition at line 1365 of file TRunge.cxx.

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

intersect_yz_plane() [2/3]

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

intersect_yz_plane() [3/3]

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

intersect_zx_plane() [1/3]

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

Definition at line 1350 of file TRunge.cxx.

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

intersect_zx_plane() [2/3]

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

intersect_zx_plane() [3/3]

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

kappa() [1/3]

double TRunge::kappa

(

void

)

const

Definition at line 241 of file TRunge.cxx.

{ return ( _a[2] ); }

kappa() [2/3]

double TRunge::kappa

(

void

)

const

kappa() [3/3]

double TRunge::kappa

(

void

)

const

Mass() [1/6]

float TRunge::Mass

(

float

mass

)

set mass for tof calc.

Definition at line 341 of file TRunge.cxx.

                               {
  _mass  = mass;
  _mass2 = _mass * _mass;
  return ( _mass );
}

Mass() [2/6]

float TRunge::Mass

(

float

)

set mass for tof calc.

Mass() [3/6]

float TRunge::Mass

(

float

)

set mass for tof calc.

Mass() [4/6]

float TRunge::Mass

(

void

)

const

return mass

Definition at line 277 of file TRunge.cxx.

{ return ( _mass ); }

Referenced by TRunge().

Mass() [5/6]

float TRunge::Mass

(

void

)

const

return mass

Mass() [6/6]

float TRunge::Mass

(

void

)

const

return mass

MaxFlightLength() [1/6]

double TRunge::MaxFlightLength

(

double

length

)

Definition at line 347 of file TRunge.cxx.

                                              {
  if ( length > 0 ) _maxflightlength = length;
  else _maxflightlength = default_maxflightlength;
  _Nstep = 0;
  return ( _maxflightlength );
}

MaxFlightLength() [2/6]

double TRunge::MaxFlightLength

(

double

)

MaxFlightLength() [3/6]

double TRunge::MaxFlightLength

(

double

)

MaxFlightLength() [4/6]

double TRunge::MaxFlightLength

(

void

)

const

return max flight length

Definition at line 279 of file TRunge.cxx.

{ return ( _maxflightlength ); }

Referenced by TRunge().

MaxFlightLength() [5/6]

double TRunge::MaxFlightLength

(

void

)

const

return max flight length

MaxFlightLength() [6/6]

double TRunge::MaxFlightLength

(

void

)

const

return max flight length

ndf() [1/3]

unsigned TRunge::ndf

(

void

)

const

returns NDF

Definition at line 255 of file TRunge.cxx.

{ return _ndf; }

Referenced by TRunge().

ndf() [2/3]

unsigned TRunge::ndf

(

void

)

const

returns NDF

ndf() [3/3]

unsigned TRunge::ndf

(

void

)

const

returns NDF

Nstep() [1/3]

unsigned TRunge::Nstep

(

void

)

const

access to trajectory cache

Definition at line 992 of file TRunge.cxx.

{ return ( _Nstep ); }

Referenced by Fly(), and Fly_SC().

Nstep() [2/3]

unsigned TRunge::Nstep

(

void

)

const

access to trajectory cache

Nstep() [3/3]

unsigned TRunge::Nstep

(

void

)

const

access to trajectory cache

objectType() [1/3]

unsigned TRunge::objectType ( void ) const

inlinevirtual

returns object type

Reimplemented from TTrackBase.

Definition at line 252 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TRunge.h.

{ return Runge; }

objectType() [2/3]

unsigned TRunge::objectType ( void ) const

virtual

returns object type

Reimplemented from TTrackBase.

objectType() [3/3]

unsigned TRunge::objectType ( void ) const

virtual

returns object type

Reimplemented from TTrackBase.

phi0() [1/3]

double TRunge::phi0

(

void

)

const

Definition at line 239 of file TRunge.cxx.

{ return ( _a[1] ); }

phi0() [2/3]

double TRunge::phi0

(

void

)

const

phi0() [3/3]

double TRunge::phi0

(

void

)

const

pivot() [1/6]

const HepPoint3D & TRunge::pivot

(

const HepPoint3D &

newpivot

)

set new pivot

!!!!! under construction !!!!! track parameter should be extracted after track propagation.

Definition at line 281 of file TRunge.cxx.

                                                            {
  /// !!!!! under construction !!!!!
  ///   track parameter should be extracted after track propagation.
  Helix tHelix( helix() );
  tHelix.pivot( newpivot );
  _pivot = newpivot;
  _a     = tHelix.a();
  _Ea    = tHelix.Ea();
  _Nstep = 0;
  return ( _pivot );
}

pivot() [2/6]

const HepPoint3D & TRunge::pivot

(

const HepPoint3D &

)

set new pivot

pivot() [3/6]

const HepPoint3D & TRunge::pivot

(

const HepPoint3D &

)

set new pivot

pivot() [4/6]

const HepPoint3D & TRunge::pivot

(

void

)

const

pivot position

Definition at line 247 of file TRunge.cxx.

{ return ( _pivot ); }

Referenced by intersect_xy_plane(), TRunge(), TRunge(), and TRunge().

pivot() [5/6]

const HepPoint3D & TRunge::pivot

(

void

)

const

pivot position

pivot() [6/6]

const HepPoint3D & TRunge::pivot

(

void

)

const

pivot position

Propagate() [1/3]

void TRunge::Propagate	(	double	y[6],
		const double &	step,
		const RkFitMaterial	material )

propagate the track using 4th-order Runge-Kutta method

Definition at line 880 of file TRunge.cxx.

                                                                                      {
  //  y[6] = (x,y,z,px,py,pz)
  double       f1[6], f2[6], f3[6], f4[6], yt[6];
  double       hh;
  double       h6;
  unsigned     i;
  const double mpi = 0.13957;
  double       me  = 0.000511;
  // eloss(step,&material, me,y,0);
  hh = step * 0.5;
  // cout<<"TR:Pro hh="<<hh<<endl;
  Function( y, f1 );
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + hh * f1[i];
  //{
  //  register double *a=y;
  //  register double *b=f1;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function( yt, f2 );
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + hh * f2[i];
  //{
  //  register double *a=y;
  //  register double *b=f2;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function( yt, f3 );
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + step * f3[i];
  //{
  //  register double *a=y;
  //  register double *b=f3;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+step*(*b);
  //}
  Function( yt, f4 );
 
  h6 = step / 6;
  for ( i = 0; i < 6; i++ ) y[i] += h6 * ( f1[i] + 2 * f2[i] + 2 * f3[i] + f4[i] );
  //{
  //  register double *a=f1;
  //  register double *b=f2;
  //  register double *c=f3;
  //  register double *d=f4;
  //  register double *t=y;
  //  register double *e=y+6;
  //  for(;t<e;a++,b++,c++,d++,t++)
  //    (*t)+=h6*((*a)+2*(*b)+2*(*c)+(*d));
  //}
}

Referenced by approach_line(), approach_point(), and Fly().

Propagate() [2/3]

void TRunge::Propagate	(	double	y[6],
		const double &	step,
		const RkFitMaterial	material )

propagate the track using 4th-order Runge-Kutta method

Propagate() [3/3]

void TRunge::Propagate	(	double	y[6],
		const double &	step,
		const RkFitMaterial	material )

propagate the track using 4th-order Runge-Kutta method

Propagate1() [1/3]

void TRunge::Propagate1	(	const double	y[6],
		const double	dydx[6],
		const double &	step,
		double	yout[6] )

Definition at line 1006 of file TRunge.cxx.

                                          {
  //  y[6] = (x,y,z,px,py,pz)
  double   f2[6], f3[6], f4[6], yt[6];
  double   hh;
  double   h6;
  unsigned i;
 
  hh = step * 0.5;
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + hh * dydx[i]; // 1st step
  //{
  //  const register double *a=y;
  //  const register double *b=dydx;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function( yt, f2 ); // 2nd step
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + hh * f2[i];
  //{
  //  const register double *a=y;
  //  const register double *b=f2;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function( yt, f3 ); // 3rd step
  for ( i = 0; i < 6; i++ ) yt[i] = y[i] + step * f3[i];
  //{
  //  const register double *a=y;
  //  const register double *b=f3;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+step*(*b);
  //}
  Function( yt, f4 ); // 4th step
 
  h6 = step / 6;
  for ( i = 0; i < 6; i++ ) yout[i] = y[i] + h6 * ( dydx[i] + 2 * f2[i] + 2 * f3[i] + f4[i] );
  //{
  //  const register double *a=dydx;
  //  const register double *b=f2;
  //  const register double *c=f3;
  //  const register double *d=f4;
  //  const register double *e=y;
  //  register double *t=yout;
  //  const register double *s=yout+6;
  //  for(;t<s;a++,b++,c++,d++,e++,t++)
  //    (*t)=(*e)+h6*((*a)+2*(*b)+2*(*c)+(*d));
  //}
}

Referenced by Propagate_QC().

Propagate1() [2/3]

void TRunge::Propagate1	(	const double	y[6],
		const double	dydx[6],
		const double &	step,
		double	yout[6] )

Propagate1() [3/3]

void TRunge::Propagate1	(	const double	y[6],
		const double	dydx[6],
		const double &	step,
		double	yout[6] )

Propagate_QC() [1/3]

void TRunge::Propagate_QC	(	double	y[6],
		double	dydx[6],
		const double &	steptry,
		const double &	eps,
		const double	yscal[6],
		double &	stepdid,
		double &	stepnext )

Definition at line 1063 of file TRunge.cxx.

                                              {
  // propagate with quality check, step will be scaled automatically
  double   ysav[6], dysav[6], ytemp[6];
  double   h, hh, errmax, temp;
  unsigned i;
 
  for ( i = 0; i < 6; i++ ) ysav[i] = y[i];
  // memcpy(ysav,y,sizeof(double)*6);
  for ( i = 0; i < 6; i++ ) dysav[i] = dydx[i];
  // memcpy(dysav,dydx,sizeof(double)*6);
 
  h = steptry;
  for ( ;; )
  {
    hh = h * 0.5; // 2 half step
    Propagate1( ysav, dysav, hh, ytemp );
    Function( ytemp, dydx );
    Propagate1( ytemp, dydx, hh, y );
    // if  check step size
    Propagate1( ysav, dysav, h, ytemp ); // 1 full step
    // error calculation
    errmax = 0.0;
    for ( i = 0; i < 6; i++ )
    {
      ytemp[i] = y[i] - ytemp[i];
      temp     = fabs( ytemp[i] / yscal[i] );
      if ( errmax < temp ) errmax = temp;
    }
    // cout<<"TR: errmax="<<errmax<<endl;
    errmax /= eps;
    if ( errmax <= 1.0 )
    { // step O.K.  calc. for next step
      stepdid  = h;
      stepnext = ( errmax > ERRCON ? SAFETY * h * exp( PGROW * log( errmax ) ) : 4.0 * h );
      break;
    }
    h = SAFETY * h * exp( PSHRNK * log( errmax ) );
  }
  for ( i = 0; i < 6; i++ ) y[i] += ytemp[i] * FCOR;
  //{
  //  register double *a=ytemp;
  //  register double *t=y;
  //  register double *e=y+6;
  //  for(;t<e;a++,t++) (*t)+=(*a)*FCOR;
  //}
}

Referenced by Fly_SC().

Propagate_QC() [2/3]

void TRunge::Propagate_QC	(	double	y[6],
		double	dydx[6],
		const double &	steptry,
		const double &	eps,
		const double	yscal[6],
		double &	stepdid,
		double &	stepnext )

Propagate_QC() [3/3]

void TRunge::Propagate_QC	(	double	y[6],
		double	dydx[6],
		const double &	steptry,
		const double &	eps,
		const double	yscal[6],
		double &	stepdid,
		double &	stepnext )

reducedchi2() [1/3]

double TRunge::reducedchi2

(

void

)

const

returns reduced-chi2

Definition at line 259 of file TRunge.cxx.

                                       {
  if ( _ndf == 0 )
  {
    std::cout << "error at TRunge::reducedchi2  ndf=0" << std::endl;
    return 0;
  }
  return ( _chi2 / _ndf );
}

reducedchi2() [2/3]

double TRunge::reducedchi2

(

void

)

const

returns reduced-chi2

reducedchi2() [3/3]

double TRunge::reducedchi2

(

void

)

const

returns reduced-chi2

SetFirst() [1/3]

void TRunge::SetFirst

(

double

y[6]

)

const

set first point (position, momentum) s=0, phi=0

Definition at line 975 of file TRunge.cxx.

                                         {
  //  y[6] = (x,y,z,px,py,pz)
  const double cosPhi0  = cos( _a[1] );
  const double sinPhi0  = sin( _a[1] );
  const double invKappa = 1 / abs( _a[2] );
 
  y[0] = _pivot.x() + _a[0] * cosPhi0; // x
  y[1] = _pivot.y() + _a[0] * sinPhi0; // y
  y[2] = _pivot.z() + _a[3];           // z
  y[3] = -sinPhi0 * invKappa;          // px
  y[4] = cosPhi0 * invKappa;           // py
  y[5] = _a[4] * invKappa;             // pz
}

Referenced by Fly(), and Fly_SC().

SetFirst() [2/3]

void TRunge::SetFirst

(

double

y[6]

)

const

set first point (position, momentum) s=0, phi=0

SetFirst() [3/3]

void TRunge::SetFirst

(

double

y[6]

)

const

set first point (position, momentum) s=0, phi=0

SetFlightLength() [1/3]

double TRunge::SetFlightLength

(

void

)

set flight length using wire hits

Definition at line 1159 of file TRunge.cxx.

                                     {
 
  // get a list of links to a wire hit
  const AList<TMLink>& cores = this->cores();
  // cout<<"TR:: cores.length="<<cores.length()<<endl;
 
  const Helix hel( this->helix() );
  double      tanl = hel.tanl();
  // double curv=hel.curv();
  // double rho = Helix::ConstantAlpha / hel.kappa();
  // double rho = 333.564095 / hel.kappa();
  const double Bz  = 1000 * getPmgnIMF()->getReferField();
  double       rho = 333.564095 / ( hel.kappa() * Bz );
 
  // search max phi
  double phi_max = -DBL_MAX;
  double phi_min = DBL_MAX;
  // loop over link
  for ( int j = 0; j < cores.length(); j++ )
  {
    TMLink& l = *cores[j];
    // fitting valid?
    const TMDCWire& wire = *l.wire();
    double          Wz   = 0;
    // double mindist;
    HepPoint3D  onWire;
    HepPoint3D  dummy_bP;
    HepVector3D dummy_dir;
    Helix       th( hel );
    for ( int i = 0; i < 10; i++ )
    {
      wire.wirePosition( Wz, onWire, dummy_bP, dummy_dir );
      th.pivot( onWire );
      if ( abs( th.dz() ) < 0.1 )
      { //<1mm
        // mindist=abs(hel.dr());
        break;
      }
      Wz += th.dz();
    }
    // onWire is closest point , th.x() is closest point on trk
    // Wz is z position of closest point
    // Z->dphi
    /*
    // Calculate position (x,y,z) along helix.
    // x = x0 + dr * cos(phi0) + (alpha / kappa) * (cos(phi0) - cos(phi0+phi))
    // y = y0 + dr * sin(phi0) + (alpha / kappa) * (sin(phi0) - sin(phi0+phi))
    // z = z0 + dz             - (alpha / kappa) * tan(lambda) * phi
    cout<<"TR:: Wz="<<Wz<<" tanl="<<tanl<<endl;
    double phi=( hel.pivot().z()+hel.dz()-Wz )/( curv*tanl );
    */
    //...Cal. dPhi to rotate...
    const HepPoint3D& xc = hel.center();
    const HepPoint3D& xt = hel.x();
    HepVector3D       v0, v1;
    v0                = xt - xc;
    v1                = th.x() - xc;
    const double vCrs = v0.x() * v1.y() - v0.y() * v1.x();
    const double vDot = v0.x() * v1.x() + v0.y() * v1.y();
    double       phi  = atan2( vCrs, vDot );
 
    double tz = hel.x( phi ).z();
    // cout<<"TR::  Wz="<<Wz<<" tz="<<tz<<endl;
 
    if ( phi > phi_max ) phi_max = phi;
    if ( phi < phi_min ) phi_min = phi;
    // cout<<"TR:: phi="<<phi<<endl;
  }
  // cout<<"TR:: phi_max="<<phi_max<<" phi_min="<<phi_min
  //     <<" rho="<<rho<<" tanl="<<tanl<<endl;
  // phi->length, set max filght length
  // tanl*=1.2; // for mergin
  _maxflightlength =
      abs( rho * ( phi_max - phi_min ) * sqrt( 1 + tanl * tanl ) ) * 1.2; // x1.1 mergin
 
  return ( _maxflightlength );
}

SetFlightLength() [2/3]

double TRunge::SetFlightLength

(

void

)

set flight length using wire hits

SetFlightLength() [3/3]

double TRunge::SetFlightLength

(

void

)

set flight length using wire hits

StepSize() [1/6]

double TRunge::StepSize

(

double

step

)

set step size to calc. trajectory

Definition at line 314 of file TRunge.cxx.

                                     {
  _stepSize = step;
  _Nstep    = 0;
  return ( _stepSize );
}

StepSize() [2/6]

double TRunge::StepSize

(

double

)

set step size to calc. trajectory

StepSize() [3/6]

double TRunge::StepSize

(

double

)

set step size to calc. trajectory

StepSize() [4/6]

double TRunge::StepSize

(

void

)

const

returns step size

Definition at line 270 of file TRunge.cxx.

{ return ( _stepSize ); }

Referenced by TRunge().

StepSize() [5/6]

double TRunge::StepSize

(

void

)

const

returns step size

StepSize() [6/6]

double TRunge::StepSize

(

void

)

const

returns step size

StepSizeMax() [1/6]

double TRunge::StepSizeMax

(

double

step

)

Definition at line 330 of file TRunge.cxx.

                                        {
  _stepSizeMax = step;
  _Nstep       = 0;
  return ( _stepSizeMax );
}

StepSizeMax() [2/6]

double TRunge::StepSizeMax

(

double

)

StepSizeMax() [3/6]

double TRunge::StepSizeMax

(

double

)

StepSizeMax() [4/6]

double TRunge::StepSizeMax

(

void

)

const

Definition at line 274 of file TRunge.cxx.

{ return ( _stepSizeMax ); }

Referenced by TRunge().

StepSizeMax() [5/6]

double TRunge::StepSizeMax

(

void

)

const

StepSizeMax() [6/6]

double TRunge::StepSizeMax

(

void

)

const

StepSizeMin() [1/6]

double TRunge::StepSizeMin

(

double

step

)

Definition at line 335 of file TRunge.cxx.

                                        {
  _stepSizeMin = step;
  _Nstep       = 0;
  return ( _stepSizeMin );
}

StepSizeMin() [2/6]

double TRunge::StepSizeMin

(

double

)

StepSizeMin() [3/6]

double TRunge::StepSizeMin

(

double

)

StepSizeMin() [4/6]

double TRunge::StepSizeMin

(

void

)

const

Definition at line 275 of file TRunge.cxx.

{ return ( _stepSizeMin ); }

Referenced by TRunge().

StepSizeMin() [5/6]

double TRunge::StepSizeMin

(

void

)

const

StepSizeMin() [6/6]

double TRunge::StepSizeMin

(

void

)

const

tanl() [1/3]

double TRunge::tanl

(

void

)

const

Definition at line 245 of file TRunge.cxx.

{ return ( _a[4] ); }

Referenced by intersect_xy_plane(), and SetFlightLength().

tanl() [2/3]

double TRunge::tanl

(

void

)

const

tanl() [3/3]

double TRunge::tanl

(

void

)

const

Yscal() [1/6]

const double * TRunge::Yscal

(

const double *

)

set error parameters for fitting with step size control

Yscal() [2/6]

const double * TRunge::Yscal

(

const double *

)

set error parameters for fitting with step size control

Yscal() [3/6]

const double * TRunge::Yscal

(

const double *

)

set error parameters for fitting with step size control

Yscal() [4/6]

const double * TRunge::Yscal

(

void

)

const

return error parameters for fitting with step size control

Definition at line 272 of file TRunge.cxx.

{ return ( _yscal ); }

Yscal() [5/6]

const double * TRunge::Yscal

(

void

)

const

return error parameters for fitting with step size control

Yscal() [6/6]

const double * TRunge::Yscal

(

void

)

const

return error parameters for fitting with step size control

Friends And Related Symbol Documentation

TRungeFitter

TRungeFitter

friend

Definition at line 222 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TRunge.h.

Referenced by SetFlightLength(), TRungeFitter, and Yscal().

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TRunge.h
• InstallArea/x86_64-el9-gcc13-opt/include/TrkReco/TRunge.h
• Reconstruction/TrkReco/include/TrkReco/TRunge.h
• TRunge.cxx