T3DLine Class Reference

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

#include <T3DLine.h>

Inheritance diagram for T3DLine:

Public Member Functions
	T3DLine ()
	Constructors.
	T3DLine (const TTrack &)
	T3DLine (const T3DLine &)
virtual	~T3DLine ()
	Destructor.
double	dr (void) const
	Track parameters.
double	phi0 (void) const
double	dz (void) const
double	tanl (void) const
double	cosPhi0 (void) const
double	sinPhi0 (void) const
Vector	a (void) const
const SymMatrix &	Ea (void) const
	returns error matrix
const HepPoint3D &	pivot (void) const
	pivot position
Helix	helix (void) const
	approximated helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
unsigned	objectType (void) const
	returns object type
HepPoint3D	x (double) const
	returns position on 3D line
HepPoint3D	x0 (void) const
	returns 3D line component x(t)=x0 + t * k
HepVector3D	k (void) const
int	approach (TMLink &, bool sagCorrection=true) const
	calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.
int	approach_line (const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack) const
	caluculate closest points between a line and this track
int	approach_point (const HepPoint3D &, HepPoint3D &onTrack) const
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
void	set (const HepPoint3D &, double t_dr, double t_phi0, double t_dz, double t_tanl)
	set track parameters,pivot
Vector	a (const Vector &)
	set track parameters
const SymMatrix &	Ea (const SymMatrix &)
	set error matrix
Public Member Functions inherited from TTrackBase
	TTrackBase ()
	Constructor.
	TTrackBase (const AList< TMLink > &links)
	Constructor.
virtual	~TTrackBase ()
	Destructor.
virtual unsigned	objectType (void) const
	returns object type.
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
	returns a list of masked TMLinks assigned to this track. 'mask' will be applied if mask is not 0.
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)
	appends TMLinks by approach. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.
void	appendByDistance (AList< TMLink > &list, double maxDistance)
	appends TMLinks by distance. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
	removes bad points by pull. The bad points are removed from the track, and are returned in 'list'.
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
	calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.
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	T3DLineFitter

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 50 of file T3DLine.h.

Constructor & Destructor Documentation

T3DLine() [1/3]

T3DLine::T3DLine

(

)

Constructors.

Definition at line 21 of file T3DLine.cxx.

  : TTrackBase(),
    _dr(0),_phi0(0),_dz(0),_tanl(0),
    _chi2(0),_ndf(0),_Ea(SymMatrix(4,0)),
    _pivot(ORIGIN){
 
  _cos_phi0 = 1;
  _sin_phi0 = 0;
 
  //...Set a default fitter...
  fitter(& T3DLine::_fitter);
 
  _fitted = false;
  _fittedWithCathode = false;
}

T3DLine() [2/3]

T3DLine::T3DLine

(

const TTrack &

a

)

Definition at line 37 of file T3DLine.cxx.

  : TTrackBase((TTrackBase &) a),
    _dr(a.helix().dr()),
    _phi0(a.helix().phi0()),
    _dz(a.helix().dz()),
    _tanl(a.helix().tanl()),
    _chi2(0),_ndf(0),_Ea(SymMatrix(4,0)),
    _pivot(a.helix().pivot()){
 
  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);
 
  //...Set a default fitter...
  fitter(& T3DLine::_fitter);
 
  _fitted = false;
  _fittedWithCathode = false;
}

T3DLine() [3/3]

T3DLine::T3DLine

(

const T3DLine &

a

)

Definition at line 59 of file T3DLine.cxx.

  : TTrackBase((TTrackBase &) a),
    _dr(a.dr()),
    _phi0(a.phi0()),
    _dz(a.dz()),
    _tanl(a.tanl()),
    _chi2(a.chi2()),_ndf(a.ndf()),
    _Ea(a.Ea()),
    _pivot(a.pivot()){
 
  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);
 
  //...Set a default fitter...
  fitter(& T3DLine::_fitter);
 
  _fitted = false;
  _fittedWithCathode = false;
}

~T3DLine()

T3DLine::~T3DLine ( )

virtual

Destructor.

Definition at line 56 of file T3DLine.cxx.

                 {
}

Member Function Documentation

a() [1/2]

Vector T3DLine::a

(

const Vector &

ta

)

set track parameters

Definition at line 187 of file T3DLine.cxx.

                                 {
  _dr = ta[0];
  _phi0 = ta[1];
  _dz = ta[2];
  _tanl = ta[3];
  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);
  return(ta);
}

a() [2/2]

Vector T3DLine::a

(

void

)

const

Definition at line 107 of file T3DLine.cxx.

                           {
  Vector ta(4);
  ta[0] = _dr;
  ta[1] = _phi0;
  ta[2] = _dz;
  ta[3] = _tanl;
  return(ta);
}

Referenced by helix().

approach()

int T3DLine::approach	(	TMLink &	l,
		bool	sagCorrection = true
	)		const

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 202 of file T3DLine.cxx.

                                                         {
 
  const TMDCWire& w=*l.wire();
  HepPoint3D xw = w.xyPosition();
  HepPoint3D wireBackwardPosition = w.backwardPosition();
  HepVector3D v = w.direction();
 
  HepPoint3D onWire,onTrack;
 
  if(approach_line(wireBackwardPosition,v,onWire,onTrack)<0)
    return(-1);
 
  // onWire,onTrack filled
  
  if(!doSagCorrection){
    l.positionOnWire(onWire);
    l.positionOnTrack(onTrack);
    return(0);              // no sag correction
  }
  // Sag correction
  //   loop for sag correction
  double onWire_y = onWire.y();
  double onWire_z = onWire.z();
 
  unsigned nTrial = 1;
  while(nTrial<100){
    w.wirePosition(onWire_z,xw,wireBackwardPosition,(HepVector3D&)v);
    if(approach_line(wireBackwardPosition,v,onWire,onTrack)<0)
      return(-1);
    if(fabs(onWire_y - onWire.y())<0.0001) break;   // |dy|< 1 micron
    onWire_y = onWire.y();
    onWire_z = onWire.z();
 
    nTrial++;
  }
 
  l.positionOnWire(onWire);
  l.positionOnTrack(onTrack);
  return(nTrial);
}

approach_line()

int T3DLine::approach_line	(	const HepPoint3D &	w0,
		const HepVector3D &	v,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack
	)		const

caluculate closest points between a line and this track

Definition at line 243 of file T3DLine.cxx.

                                                            {
  //  line = [w0] + s * [v]    -> [onLine]
  //  trk  = [x0] + t * [k]    -> [onTrack]
  //  if [v]//[k] then return(-1) error
 
  const HepVector3D k = this->k();
  const double v_k = v.dot(k);
  const double v_2 = v.mag2();
  const double k_2 = k.mag2();
  const double tk = v_k*v_k - v_2*k_2;
  if(tk==0) return(-1);
 
  const HepPoint3D x0 = this->x0();
  const HepPoint3D dx = x0 - w0;
  const double t = dx.dot( v_2 * k - v_k * v)/tk;
  const double s = dx.dot( v_k * k - k_2 * v)/tk;
 
  onLine = w0 + s * v;
  onTrack = x0 + t * k;
  return(0);
}

Referenced by approach().

approach_point()

int T3DLine::approach_point	(	const HepPoint3D &	p0,
		HepPoint3D &	onTrack
	)		const

caluculate closest point between a point and this track

Definition at line 266 of file T3DLine.cxx.

                                                                         {
  //  trk  = [x0] + t * [k]    -> [onTrack]
  //  if [v]//[k] then return(-1) error
 
  const HepVector3D k = this->k();
  const HepPoint3D x0 = this->x0();
  const HepPoint3D dx = p0 - x0;
  const double t = dx.dot(k)/k.mag2();
 
  onTrack = x0 + t * k;
  return(0);
}

chi2()

double T3DLine::chi2

(

void

)

const

returns chi2.

Definition at line 131 of file T3DLine.cxx.

                              {
  return _chi2;
}

cosPhi0()

double T3DLine::cosPhi0

(

void

)

const

Definition at line 95 of file T3DLine.cxx.

                                 {
  return _cos_phi0;
}

dr()

double T3DLine::dr

(

void

)

const

Track parameters.

Definition at line 79 of file T3DLine.cxx.

                            {
  return _dr;
}

Referenced by pivot().

dz()

double T3DLine::dz

(

void

)

const

Definition at line 87 of file T3DLine.cxx.

                            {
  return _dz;
}

Referenced by pivot().

Ea() [1/2]

const SymMatrix & T3DLine::Ea

(

const SymMatrix &

tEa

)

set error matrix

Definition at line 197 of file T3DLine.cxx.

                                                {
  _Ea = tEa;
  return(_Ea);
}

Ea() [2/2]

const SymMatrix & T3DLine::Ea

(

void

)

const

returns error matrix

Definition at line 116 of file T3DLine.cxx.

                                      {
  return(_Ea);
}

helix()

Helix T3DLine::helix

(

void

)

const

approximated helix class

Definition at line 120 of file T3DLine.cxx.

                              {
  HepVector a(5);
  a[0]=_dr;a[1]=_phi0;a[2]=1e-10;a[3]=_dz;a[4]=_tanl;
  Helix _helix(_pivot,a);
  return _helix;
}

Referenced by TBuilderCosmic::buildStereo().

k()

HepVector3D T3DLine::k

(

void

)

const

Definition at line 158 of file T3DLine.cxx.

                                {
  HepPoint3D p(-_sin_phi0,_cos_phi0,_tanl);
  return p;
}

Referenced by approach_line(), and approach_point().

ndf()

unsigned T3DLine::ndf

(

void

)

const

returns NDF

Definition at line 127 of file T3DLine.cxx.

                               {
  return _ndf;
}

objectType()

unsigned T3DLine::objectType ( void  ) const

inlinevirtual

returns object type

Reimplemented from TTrackBase.

Definition at line 150 of file T3DLine.h.

                                             {
  return Line3D;
}

phi0()

double T3DLine::phi0

(

void

)

const

Definition at line 83 of file T3DLine.cxx.

                              {
  return _phi0;
}

pivot() [1/2]

const HepPoint3D & T3DLine::pivot

(

const HepPoint3D &

newpivot

)

set new pivot

Definition at line 163 of file T3DLine.cxx.

                                                          {
  double dr=(_pivot.x()-newpivot.x())*_cos_phi0
           +(_pivot.y()-newpivot.y())*_sin_phi0 + _dr;
  double dz=_pivot.z()-newpivot.z()+_dz
           +_tanl*((_pivot.x()-newpivot.x())*_sin_phi0
          +(newpivot.y()-_pivot.y())*_cos_phi0);
  _dr=dr;
  _dz=dz;
  _pivot=newpivot;
  return _pivot;
}

pivot() [2/2]

const HepPoint3D & T3DLine::pivot

(

void

)

const

pivot position

Definition at line 103 of file T3DLine.cxx.

                                          {
  return _pivot;
}

reducedchi2()

double T3DLine::reducedchi2

(

void

)

const

returns reduced-chi2

Definition at line 135 of file T3DLine.cxx.

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

set()

void T3DLine::set	(	const HepPoint3D &	t_pivot,
		double	t_dr,
		double	t_phi0,
		double	t_dz,
		double	t_tanl
	)

set track parameters,pivot

Definition at line 175 of file T3DLine.cxx.

                                                              {
 
  _pivot = t_pivot;
  _dr = t_dr;
  _phi0 = t_phi0;
  _dz = t_dz;
  _tanl = t_tanl;
  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);
}

sinPhi0()

double T3DLine::sinPhi0

(

void

)

const

Definition at line 99 of file T3DLine.cxx.

                                 {
  return _sin_phi0;
}

tanl()

double T3DLine::tanl

(

void

)

const

Definition at line 91 of file T3DLine.cxx.

                              {
  return _tanl;
}

x()

HepPoint3D T3DLine::x

(

double

t

)

const

returns position on 3D line

Definition at line 143 of file T3DLine.cxx.

                                   {
  double tx= _pivot.x() + _dr * _cos_phi0 - t * _sin_phi0;
  double ty= _pivot.y() + _dr * _sin_phi0 + t * _cos_phi0;
  double tz= _pivot.z() + _dz             + t * _tanl;
  HepPoint3D p(tx,ty,tz);
  return p;
}

x0()

HepPoint3D T3DLine::x0

(

void

)

const

returns 3D line component x(t)=x0 + t * k

Definition at line 151 of file T3DLine.cxx.

                                {
  double tx= _pivot.x() + _dr * _cos_phi0;
  double ty= _pivot.y() + _dr * _sin_phi0;
  double tz= _pivot.z() + _dz;
  HepPoint3D p(tx,ty,tz);
  return p;
}

Referenced by approach_line(), and approach_point().

Friends And Related Function Documentation

T3DLineFitter

friend class T3DLineFitter

friend

Definition at line 147 of file T3DLine.h.

---

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

• T3DLine.h
• T3DLine.cxx