TTrack Class Reference

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

#include <TTrack.h>

Inheritance diagram for TTrack:

Public Member Functions
	TTrack ()
	Default constructor.
	TTrack (const TCircle &)
	Constructor.
	TTrack (const TTrack &)
	Copy constructor.
	TTrack (const T3DLine &)
	Constructor.
	TTrack (const Helix &)
	Constructor.
	TTrack (const TRunge &)
virtual	~TTrack ()
	Destructor.
const std::string &	name (void) const
	returns/sets name.
const std::string &	name (const std::string &newName)
TTrack *	mother (void) const
	sets/returns mother/daughter.
TTrack *	mother (TTrack *)
TTrack *	daughter (void) const
TTrack *	daughter (TTrack *)
unsigned	objectType (void) const
	returns type.
unsigned	type (void) const
	returns type. Definition is depending on an object type.
unsigned	finder (void) const
	sets/returns finder.
unsigned	finder (unsigned finderMask)
void	setFinderType (unsigned)
unsigned	getFinderType (void) const
unsigned	quality (void) const
	sets/returns quality.
unsigned	quality (unsigned qualityMask)
unsigned	fitting (void) const
	sets/returns fitting status.
unsigned	fitting (unsigned fitMask)
unsigned	state (void) const
	returns/sets internal state.(for bank output)
void	assign (unsigned maskForWireHit)
	assigns wire hits to this track.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
double	charge (void) const
	returns charge.
const Helix &	helix (void) const
	returns helix parameter.
Hep3Vector	p (void) const
	returns momentum.
double	ptot (void) const
	returns magnitude of momentum.
double	pt (void) const
	returns Pt.
double	pz (void) const
	returns Pz.
double	impact (void) const
	returns signed impact parameter to the origin.
TPoint2D	center (void) const
	returns position of helix center.
double	radius (void) const
	returns signed radius.
unsigned	ndf (void) const
	returns NDF.
double	chi2 (void) const
	returns chi2.
double	confidenceLevel (void) const
	returns confidence level.
AList< TSegment > &	segments (void)
	returns AList<TSegment>.
const AList< TSegment > &	segments (void) const
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.
int	approach (TMLink &, bool sagCorrection) const
int	approach2D (TMLink &) const
int	szPosition (TMLink &link) const
	calculates arc length and z for a stereo hit.
int	szPosition (const TSegment &segment, TMLink &link) const
	calculates arc length and z for a segment. Results are stored in TMLink.
int	szPosition (const HepPoint3D &p, HepPoint3D &szPosition) const
	calculates arc length for a point.
int	stereoHitForCurl (TMLink &link, AList< HepPoint3D > &arcZList) const
int	stereoHitForCurl (TMLink &link, TMLink &link1) const
int	stereoHitForCurl (TMLink &link, TMLink &link1, TMLink &link2) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2, AList< HepPoint3D > &l3) const
int	stereoHitForCurl (AList< TMLink > &) const
	calculates arc length and z for a stereo hit. uses these functions for curl tracks(included svd version).
double	charge (double)
	sets charge.
int	fit2D (unsigned=0, double=0.1, double=0.015)
	fits itself. Error was happened if return value is not zero.
void	refine2D (AList< TMLink > &list, float maxSigma)
	fits itself with cathode hits.
void	movePivot (void)
	moves pivot to the inner most hit.
int	HelCyl (double rhole, double rcyl, double zb, double zf, double epsl, double &phi, HepPoint3D &xp) const
	returns a cathode hit list.
const AList< TMLink > &	finalHits (void) const
	finds cathode hits associated to this track.
const AList< TMLink > &	finalHits (const AList< TMLink > &hits)
const AList< TMLink > &	associateHits (void) const
	sets/returns a list of associated TMLink which are used for table output.
const AList< TMLink > &	associateHits (const AList< TMLink > &hits)
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	TTrackManager
class	TrkReco
class	THelixFitter
class	TCosmicFitter
class	Refit
class	TBuilder
class	TBuilder0
class	TBuilderCosmic
class	TBuilderCurl
class	TPMCurlFinder
class	TCurlFinder

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 129 of file TTrack.h.

Constructor & Destructor Documentation

TTrack() [1/6]

TTrack::TTrack

(

)

Default constructor.

Definition at line 207 of file TTrack.cxx.

: TTrackBase(),
  _charge(1.),
  _helix(new Helix(ORIGIN, Vector(5, 0), SymMatrix(5, 0))),
  _ndf(0),
  _chi2(0.),
  _name("empty track"),
  _state(0),
  _mother(0),
  _daughter(0) {
  //jialk
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
  if(scmgn!=StatusCode::SUCCESS) { 
    std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
  }
}

TTrack() [2/6]

TTrack::TTrack

(

const TCircle &

c

)

Constructor.

Definition at line 81 of file TTrack.cxx.

: TTrackBase(c.links()),
  _helix(new Helix(ORIGIN, Vector(5, 0), SymMatrix(5, 0))), 
  _charge(c.charge()),
  _ndf(0),
  _chi2(0.),
  _name("none"),
  _type(0),
  _state(0),
  _mother(0),
  _daughter(0) {
 
    //jialk
    StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
    if(scmgn!=StatusCode::SUCCESS) { 
      std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
    }
 
    //_fitter.setMagneticFieldPointer(m_pmgnIMF);//yzhang add 2012-05-04 
    //cout<<"TTrack: "<<m_pmgnIMF->getReferField()<<endl;
    //...Set a defualt fitter...
    fitter(& TTrack::_fitter);
 
    //...Calculate helix parameters...
    Vector a(5);
    a[1] = fmod(atan2(_charge * (c.center().y() - ORIGIN.y()),
              _charge * (c.center().x() - ORIGIN.x()))
        + 4. * M_PI,
        2. * M_PI);
   // a[2] = Helix::ConstantAlpha / c.radius();
   // a[2] = 333.564095 / c.radius();
    const double Bz = -1000*m_pmgnIMF->getReferField();
    a[2] = 333.564095/ (c.radius() * Bz);
    a[0] = (c.center().x() - ORIGIN.x()) / cos(a[1]) - c.radius();
    a[3] = 0.;
    a[4] = 0.;
    _helix->a(a);
 
    //...Update links...
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++)
    _links[i]->track(this);
 
    _fitted = false;
    _fittedWithCathode = false;
/*
    //jialk
    StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
    if(scmgn!=StatusCode::SUCCESS) { 
      std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
    }
*/
}

TTrack() [3/6]

TTrack::TTrack

(

const TTrack &

a

)

Copy constructor.

Definition at line 135 of file TTrack.cxx.

: TTrackBase((TTrackBase &) a),
  _charge(a._charge),
  _segments(a._segments),
  _helix(new Helix(* a._helix)),
  _ndf(a._ndf),
  _chi2(a._chi2),
  _name("copy of " + a._name),
  _type(a._type),
//  _catHits(a._catHits),
  _state(a._state),
  _mother(a._mother),
  _daughter(a._daughter) { 
  //jialk
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
  if(scmgn!=StatusCode::SUCCESS) { 
    std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
  }
}

TTrack() [4/6]

TTrack::TTrack

(

const T3DLine &

)

Constructor.

TTrack() [5/6]

TTrack::TTrack

(

const Helix &

h

)

Constructor.

Definition at line 181 of file TTrack.cxx.

: TTrackBase(),
  _helix(new Helix(h)), 
  _ndf(0),
  _chi2(0.),
  _name("none"),
  _type(0),
  _state(0),
  _mother(0),
  _daughter(0) {
 
    //...Set a defualt fitter...
    fitter(& TTrack::_fitter);
 
    if(_helix->kappa() > 0.)_charge = 1.;
    else _charge = -1.;
 
    _fitted = false;
    _fittedWithCathode = false;
    //jialk
    StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
    if(scmgn!=StatusCode::SUCCESS) { 
      std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
    }
}

TTrack() [6/6]

TTrack::TTrack

(

const TRunge &

a

)

Definition at line 168 of file TTrack.cxx.

: TTrackBase((TTrackBase &) a),
  _helix(new Helix( a.helix())),
  _ndf(a.ndf()),
  _chi2(a.chi2()),
  _name("no"),
  _type(0),
  _state(0),
  _mother(0),
  _daughter(0) {
    if(_helix->kappa() > 0.)_charge = 1.;
        else _charge = -1.;
}

~TTrack()

TTrack::~TTrack ( )

virtual

Destructor.

Definition at line 224 of file TTrack.cxx.

                {
    delete _helix;
}

Member Function Documentation

approach() [1/2]

int TTrack::approach ( TMLink &  l ) const

virtual

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

Reimplemented from TTrackBase.

Definition at line 296 of file TTrack.cxx.

                                 {
    return approach(l, false);
}

Referenced by approach(), and TTrackManager::closest().

approach() [2/2]

int TTrack::approach	(	TMLink &	link,
		bool	sagCorrection
	)		const

Definition at line 3059 of file TTrack.cxx.

                                                          {
    //...Cal. dPhi to rotate...
    const TMDCWire & w = * link.wire();
    double wp[3]; w.xyPosition(wp);
    double wb[3]; w.backwardPosition(wb);
    double v[3];
    v[0] = w.direction().x();
    v[1] = w.direction().y();
    v[2] = w.direction().z();
    //...Sag correction...
    if (doSagCorrection) {
    HepVector3D dir = w.direction();
    HepPoint3D xw(wp[0], wp[1], wp[2]);
    HepPoint3D wireBackwardPosition(wb[0], wb[1], wb[2]);
    w.wirePosition(link.positionOnTrack().z(),
               xw,
               wireBackwardPosition,
               dir);
    v[0] = dir.x();
    v[1] = dir.y();
    v[2] = dir.z();
    wp[0] = xw.x();
    wp[1] = xw.y();
    wp[2] = xw.z();
    wb[0] = wireBackwardPosition.x();
    wb[1] = wireBackwardPosition.y();
    wb[2] = wireBackwardPosition.z();
    }
    //...Cal. dPhi to rotate...
    const HepPoint3D & xc = _helix->center();
    double xt[3]; _helix->x(0., xt);
    double x0 = - xc.x();
    double y0 = - xc.y();
    double x1 = wp[0] + x0;
    double y1 = wp[1] + y0;
    x0 += xt[0];
    y0 += xt[1];
    double dPhi = atan2(x0 * y1 - y0 * x1, x0 * x1 + y0 * y1);
    //...Setup...
    double kappa = _helix->kappa();
    double phi0 = _helix->phi0();
//yzhang 2012-08-30
    //...Axial case...
//    if (w.axial()) {
//  link.positionOnTrack(_helix->x(dPhi));
//  HepPoint3D x(wp[0], wp[1], wp[2]);
//  x.setZ(link.positionOnTrack().z());
//  link.positionOnWire(x);
//  link.dPhi(dPhi);
//  return 0;
//    }
#ifdef TRKRECO_DEBUG
    double firstdfdphi = 0.;
    static bool first = true;
    if (first) {
//  extern BelleTupleManager * BASF_Histogram;
//  BelleTupleManager * m = BASF_Histogram;
//  h_nTrial = m->histogram("TTrack::approach nTrial", 100, 0., 100.);
    }
#endif
 
    //...Stereo case...
    // double rho = Helix::ConstantAlpha / kappa;
    // double rho = 333.564095 / kappa;
    // yzhang 2012-05-03 delete
    //double rho = 333.564095/  kappa;
    //yzhang add 
    Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF);
    if(m_pmgnIMF==NULL) {
      std::cout<< __FILE__<<" "<<__LINE__<<" Unable to open Magnetic field service"<<std::endl;
    }
    const double Bz = -1000*m_pmgnIMF->getReferField();
    double rho   = 333.564095/(kappa * Bz); 
    //zhangy
    double tanLambda = _helix->tanl();
    static Vector x(3);
    double t_x[3];
    double t_dXdPhi[3];
    const double convergence = 1.0e-5;
    double l;
    unsigned nTrial = 0;
    while (nTrial < 100) {
 
    x = link.positionOnTrack(_helix->x(dPhi));
    t_x[0] = x[0];
    t_x[1] = x[1];
    t_x[2] = x[2];
 
        l = v[0] * t_x[0] + v[1] * t_x[1] + v[2] * t_x[2]
        - v[0] * wb[0] - v[1] * wb[1] - v[2] * wb[2];
 
    double rcosPhi = rho * cos(phi0 + dPhi);
    double rsinPhi = rho * sin(phi0 + dPhi);
    t_dXdPhi[0] =   rsinPhi;
    t_dXdPhi[1] = - rcosPhi;
    t_dXdPhi[2] = - rho * tanLambda;
 
    //...f = d(Distance) / d phi...
    double t_d2Xd2Phi[2];
    t_d2Xd2Phi[0] = rcosPhi;
    t_d2Xd2Phi[1] = rsinPhi;
 
    //...iw new...
    double n[3];
    n[0] = t_x[0] - wb[0];
    n[1] = t_x[1] - wb[1];
    n[2] = t_x[2] - wb[2];
    
    double a[3];
    a[0] = n[0] - l * v[0];
    a[1] = n[1] - l * v[1];
    a[2] = n[2] - l * v[2];
    double dfdphi = a[0] * t_dXdPhi[0]
        + a[1] * t_dXdPhi[1]
        + a[2] * t_dXdPhi[2];
 
#ifdef TRKRECO_DEBUG
    if (nTrial == 0) {
//      break;
        firstdfdphi = dfdphi;
    }
 
    //...Check bad case...
    if (nTrial > 3) {
        std::cout << "TTrack::approach ... " << w.name() << " "
              << "dfdphi(0)=" << firstdfdphi
              << ",(" << nTrial << ")=" << dfdphi << endl;
    }
#endif
 
    //...Is it converged?...
    if (fabs(dfdphi) < convergence)
        break;
 
    double dv = v[0] * t_dXdPhi[0]
        + v[1] * t_dXdPhi[1]
        + v[2] * t_dXdPhi[2];
    double t0 = t_dXdPhi[0] * t_dXdPhi[0]
        + t_dXdPhi[1] * t_dXdPhi[1]
        + t_dXdPhi[2] * t_dXdPhi[2];
    double d2fd2phi = t0 - dv * dv
        + a[0] * t_d2Xd2Phi[0]
        + a[1] * t_d2Xd2Phi[1];
//      + a[2] * t_d2Xd2Phi[2];
 
    dPhi -= dfdphi / d2fd2phi;
 
//      cout << "nTrial=" << nTrial << endl;
//      cout << "iw f,df,dphi=" << dfdphi << "," << d2fd2phi << "," << dPhi << endl;
 
    ++nTrial;
    }
 
    //...Cal. positions...
    link.positionOnWire(HepPoint3D(wb[0] + l * v[0],
                wb[1] + l * v[1],
                wb[2] + l * v[2]));
    link.dPhi(dPhi);
 
#ifdef TRKRECO_DEBUG
//    h_nTrial->accumulate((float) nTrial + .5);
#endif
 
    return nTrial;
}

approach2D()

int TTrack::approach2D

(

TMLink &

l

)

const

Definition at line 2891 of file TTrack.cxx.

                                  {
 
  const TMDCWire &w = *l.wire();
  double kappa = _helix->a()[2];
  double phi0  = _helix->a()[1];
  HepPoint3D xc(_helix->center());
  HepPoint3D xw(w.xyPosition());
  HepPoint3D xt(_helix->x());
  HepVector3D v0(xt-xc);
  HepVector3D v1(xw-xc);
 
  double vCrs = v0.x() * v1.y() - v0.y() * v1.x();
  double vDot = v0.x() * v1.x() + v0.y() * v1.y();
  double dPhi = atan2(vCrs, vDot);
 
  xt = _helix->x(dPhi);
  xt.setZ(0.);
  l.positionOnTrack(xt);
  xw.setZ(0.);
  l.positionOnWire(xw);
  l.dPhi(dPhi);
  return 0;
}

Referenced by fit2D().

assign()

void TTrack::assign

(

unsigned

maskForWireHit

)

assigns wire hits to this track.

Definition at line 3635 of file TTrack.cxx.

                               {
    hitMask |= WireHitUsed;
 
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink * l = _links[i];
    const TMDCWireHit * h = l->hit();
#ifdef TRKRECO_DEBUG
    if (h->track()) {
        std::cout << "TTrack::assign !!! hit(" << h->wire()->name();
        std::cout << ") already assigned" << std::endl;
    }
#endif
 
    //...This function will be moved to TTrackManager...
    h->track((TTrack *) this);
    h->state(h->state() | hitMask);
    }
}

associateHits() [1/2]

const AList< TMLink > & TTrack::associateHits ( const AList< TMLink > &  hits )

inline

Definition at line 717 of file TTrack.h.

                                                {
    _associateHits = list;
    return _associateHits;
}

associateHits() [2/2]

const AList< TMLink > & TTrack::associateHits ( void  ) const

inline

sets/returns a list of associated TMLink which are used for table output.

Definition at line 724 of file TTrack.h.

                                {
    return _associateHits;
}

center()

TPoint2D TTrack::center ( void  ) const

inline

returns position of helix center.

Definition at line 595 of file TTrack.h.

                         {
    return TPoint2D(_helix->center());
}

Referenced by stereoHitForCurl(), and szPosition().

charge() [1/2]

double TTrack::charge ( double  a )

inline

sets charge.

Definition at line 510 of file TTrack.h.

                       {
    return _charge = a;
}

charge() [2/2]

double TTrack::charge ( void  ) const

inline

returns charge.

Definition at line 504 of file TTrack.h.

                         {
    return _charge;
}

Referenced by TBuilder::buildStereo(), TBuilder0::buildStereo(), TBuilderCosmic::buildStereo(), TBuilderCurl::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), calVirtualCircle(), moveLR(), and TTrackManager::setCurlerFlags().

chi2()

double TTrack::chi2 ( void  ) const

inline

returns chi2.

Definition at line 495 of file TTrack.h.

                       {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TTrack::chi2 !!! chi2 not updated" << std::endl;
#endif
    return _chi2;
}

Referenced by fit2D(), and TTrackManager::sortTracksByQuality().

confidenceLevel()

double TTrack::confidenceLevel ( void  ) const

inline

returns confidence level.

Definition at line 547 of file TTrack.h.

                                  {
    return chisq2confLevel((int) _ndf, _chi2);
}

daughter() [1/2]

TTrack * TTrack::daughter ( TTrack *  a )

inline

Definition at line 663 of file TTrack.h.

                           {
    if (a) _state |= (TrackHasDaughter << TrackRelationShift);
    else   _state &= (~(TrackHasDaughter << TrackRelationShift));
    return _daughter = a;
}

daughter() [2/2]

TTrack * TTrack::daughter ( void  ) const

inline

Definition at line 657 of file TTrack.h.

                           {
    return _daughter;
}

Referenced by TTrackManager::setCurlerFlags().

deleteListForCurl() [1/2]

void TTrack::deleteListForCurl	(	AList< HepPoint3D > &	l1,
		AList< HepPoint3D > &	l2
	)		const

deleteListForCurl() [2/2]

void TTrack::deleteListForCurl	(	AList< HepPoint3D > &	l1,
		AList< HepPoint3D > &	l2,
		AList< HepPoint3D > &	l3
	)		const

dump()

void TTrack::dump ( const std::string &  message = std::string(""), const std::string &  prefix = std::string("")  ) const

virtual

dumps debug information.

Reimplemented from TTrackBase.

Definition at line 229 of file TTrack.cxx.

                                                                {
    bool def = false;
    if (msg == "") def = true;
    std::string pre = pre0;
    std::string tab;
    for (unsigned i = 0; i < pre.length(); i++)
    tab += " ";
    if (def || msg.find("track") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre
          << p() << "(pt=" << pt() << ")" << std::endl
          << tab
          << "links=" << _links.length()
          << "(cores=" << nCores()
          << "),chrg=" << _charge
          << ",ndf=" << _ndf
          << ",chi2=" << _chi2
          << std::endl;
    pre = tab;
    }
    if (msg.find("helix") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre
          << "pivot=" << _helix->pivot()
          << ",center=" << _helix->center() << std::endl
          << pre
          << "helix=(" << _helix->a()[0] << "," << _helix->a()[1]<< ","
          << _helix->a()[2] << "," << _helix->a()[3] << ","
          << _helix->a()[4] << ")" << std::endl;
    pre = tab;
    }
    if (! def) TTrackBase::dump(msg, pre);
}

finalHits() [1/2]

const AList< TMLink > & TTrack::finalHits ( const AList< TMLink > &  hits )

inline

Definition at line 704 of file TTrack.h.

                                            {
    _finalHits = list;
    return _finalHits;
}

finalHits() [2/2]

const AList< TMLink > & TTrack::finalHits ( void  ) const

inline

finds cathode hits associated to this track.

relation between cluster and wire relation between cluster and layer appends Svd Hits(TSvdHit) returns Svd Hits List sets/returns a list of TMLink which are used for table output.

Definition at line 711 of file TTrack.h.

                            {
    return _finalHits;
}

finder() [1/2]

unsigned TTrack::finder ( unsigned  finderMask )

inline

Definition at line 607 of file TTrack.h.

                         {
    _state |= (a & TrackFinderMask);
    return finder();
}

finder() [2/2]

unsigned TTrack::finder ( void  ) const

inline

sets/returns finder.

Definition at line 601 of file TTrack.h.

                         {
    return _state & TrackFinderMask;
}

Referenced by finder(), and TTrackManager::salvage().

fit2D()

int TTrack::fit2D	(	unsigned	ipFlag = 0,
		double	ipDistance = 0.1,
		double	ipError = 0.015
	)

fits itself. Error was happened if return value is not zero.

fits itself with r-phi view. Error was happened if return value is not zero.

Definition at line 2542 of file TTrack.cxx.

                                                                {
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    TTrack::fit2D(r-phi) ..." << std::endl;
#endif
  //if(_fitted)return 0;
 
  //...Check # of hits...
 
  //std::cout << "# = " << _links.length() << std::endl;
  //...Setup...
  unsigned nTrial(0);
  Vector a(_helix->a());
  double chi2;
  double chi2Old = 1.0e+99;
  Vector dchi2da(3);
  SymMatrix d2chi2d2a(3,0);
  Vector da(5), dxda(3), dyda(3);
  Vector dDda(3);
  const double convergence = 1.0e-5;
  Vector f(3);
  int err = 0;
  double factor = 1.0;
  unsigned usedWireNumber = 0;  
 
  //...Fitting loop...
  while(nTrial < 100){
    //...Set up...
    chi2 = 0.;
    for (unsigned j=0;j<3;++j) dchi2da[j] = 0.;
    d2chi2d2a = SymMatrix(3, 0);
    usedWireNumber = 0;
 
    //...Loop with hits...
    unsigned i = 0;
    while (TMLink * l = _links[i++]) {
      if(l->fit2D() == 0)continue;
      const TMDCWireHit &h = *l->hit();
      
      //...Cal. closest points...
      if(approach2D(*l) != 0)continue;
      double dPhi = l->dPhi();
      const HepPoint3D & onTrack = l->positionOnTrack();
      const HepPoint3D & onWire  = l->positionOnWire();
      HepPoint3D onTrack2(onTrack.x(),onTrack.y(),0.);
      HepPoint3D onWire2(onWire.x(),onWire.y(),0.);
 
      //...Obtain drift distance and its error...
      unsigned leftRight = WireHitRight;
      if (onWire2.cross(onTrack2).z() < 0.) leftRight = WireHitLeft;
      double distance   = h.drift(leftRight);
      double eDistance  = h.dDrift(leftRight);
      double eDistance2 = eDistance * eDistance;
      if(eDistance == 0.){
    std::cout << "Error(?) : Drift Distance Error = 0 in fit2D of TrkReco." << std::endl;
    continue;
      }
 
      //...Residual...
      HepVector3D v = onTrack2 - onWire2;
      double vmag = v.mag();
      double dDistance = vmag - distance;
 
      //...dxda...
      if(this->dxda2D(*l, dPhi, dxda, dyda) != 0)continue;
 
      //...Chi2 related...
      //Vector3 vw(0.,0.,1.);
      dDda = (vmag > 0.)
    ? (v.x()*dxda+v.y()*dyda)/vmag
    : Vector(3,0);
      if(vmag<=0.0){
    std::cout << "    in fit2D " << onTrack << ", " << onWire;
    h.dump();
    continue;
      }
      dchi2da     += (dDistance/eDistance2)*dDda;
      d2chi2d2a   += vT_times_v(dDda)/eDistance2;
      double pChi2 = dDistance*dDistance/eDistance2;
      chi2        += pChi2;
      
      //...Store results...
      l->update(onTrack2, onWire2, leftRight, pChi2);
      ++usedWireNumber;
    }
    if(ipFlag != 0){
      double kappa = _helix->a()[2];
      double phi0  = _helix->a()[1];
      HepPoint3D xc(_helix->center());
      HepPoint3D onWire(0.,0.,0.);
      HepPoint3D xt(_helix->x());
      HepVector3D v0(xt-xc);
      HepVector3D v1(onWire-xc);
      double vCrs = v0.x() * v1.y() - v0.y() * v1.x();
      double vDot = v0.x() * v1.x() + v0.y() * v1.y();
      double dPhi = atan2(vCrs, vDot);      
      HepPoint3D onTrack(_helix->x(dPhi).x(),_helix->x(dPhi).y(),0.);
      double distance   = ipDistance;
      double eDistance  = ipError;
      double eDistance2 = eDistance * eDistance;
 
      HepVector3D v = onTrack - onWire;
      double vmag = v.mag();
      double dDistance = vmag - distance;
 
      if(this->dxda2D(dPhi, dxda, dyda) != 0)goto ipOff;
      
      dDda = (vmag > 0.)
    ? (v.x()*dxda+v.y()*dyda)/vmag
    : Vector(3,0);
      if(vmag<=0.0){
    goto ipOff;
      }
      dchi2da     += (dDistance/eDistance2)*dDda;
      d2chi2d2a   += vT_times_v(dDda)/eDistance2;
      double pChi2 = dDistance*dDistance/eDistance2;
      chi2        += pChi2;
      
      ++usedWireNumber;
    }
  ipOff:
    if(usedWireNumber < 4){
      err = -2;
      break;
    }
 
    //...Check condition...
    double change = chi2Old - chi2;
    if(fabs(change) < convergence)break;
    if(change < 0.){
#ifdef TRKRECO_DEBUG_DETAIL
      std::cout << "chi2Old, chi2=" << chi2Old <<" "<< chi2 << std::endl;
#endif
      //change to the old value.
      a += factor*da;
      _helix->a(a);
 
      chi2 = 0.;
      for (unsigned j=0;j<3;++j) dchi2da[j] = 0.;
      d2chi2d2a = SymMatrix(3,0);
      usedWireNumber = 0;
   
      //...Loop with hits...
      unsigned i = 0;
      while (TMLink *l = _links[i++]) {
    if(l->fit2D() == 0)continue;
    const TMDCWireHit & h = * l->hit();
    
    //...Cal. closest points...
    if(approach2D(*l) != 0)continue;
    double dPhi = l->dPhi();
    const HepPoint3D & onTrack = l->positionOnTrack();
    const HepPoint3D & onWire  = l->positionOnWire();
    HepPoint3D onTrack2(onTrack.x(),onTrack.y(),0.);
    HepPoint3D onWire2(onWire.x(),onWire.y(),0.);
    
    //...Obtain drift distance and its error...
    unsigned leftRight = WireHitRight;
    if (onWire2.cross(onTrack2).z() < 0.) leftRight = WireHitLeft;
    double distance   = h.drift(leftRight);
    double eDistance  = h.dDrift(leftRight);
    double eDistance2 = eDistance * eDistance;
            
    //...Residual...
    HepVector3D v = onTrack2 - onWire2;
    double vmag = v.mag();
    double dDistance = vmag - distance;
 
    //...dxda...
    if(this->dxda2D(*l, dPhi, dxda, dyda) != 0)continue;
 
    //...Chi2 related...
    dDda = (vmag>0.)
      ? (v.x()*dxda + v.y()*dyda)/vmag
      : Vector(3,0);
    if(vmag<=0.0) {
      std::cout << "    in fit2D " << onTrack << ", " << onWire;
      h.dump();
      continue;
    }
    dchi2da     += (dDistance/eDistance2)*dDda;
    d2chi2d2a   += vT_times_v(dDda)/eDistance2;
    double pChi2 = dDistance*dDistance/eDistance2;
    chi2        += pChi2;
    
    //...Store results...
    l->update(onTrack2, onWire2, leftRight, pChi2);
    ++usedWireNumber;
      }
      if(ipFlag != 0){
    double kappa = _helix->a()[2];
    double phi0  = _helix->a()[1];
    HepPoint3D xc(_helix->center());
    HepPoint3D onWire(0.,0.,0.);
    HepPoint3D xt(_helix->x());
    HepVector3D v0(xt-xc);
    HepVector3D v1(onWire-xc);
    double vCrs = v0.x() * v1.y() - v0.y() * v1.x();
    double vDot = v0.x() * v1.x() + v0.y() * v1.y();
    double dPhi = atan2(vCrs, vDot);      
    HepPoint3D onTrack(_helix->x(dPhi).x(),_helix->x(dPhi).y(),0.);
    double distance   = ipDistance;
    double eDistance  = ipError;
    double eDistance2 = eDistance * eDistance;
    
    HepVector3D v = onTrack - onWire;
    double vmag = v.mag();
    double dDistance = vmag - distance;
    
    if(this->dxda2D(dPhi, dxda, dyda) != 0)goto ipOff2;
    
    dDda = (vmag > 0.)
      ? (v.x()*dxda+v.y()*dyda)/vmag
      : Vector(3,0);
    if(vmag<=0.0){
      goto ipOff2;
    }
    dchi2da     += (dDistance/eDistance2)*dDda;
    d2chi2d2a   += vT_times_v(dDda)/eDistance2;
    double pChi2 = dDistance*dDistance/eDistance2;
    chi2        += pChi2;
    
    ++usedWireNumber;
      }
    ipOff2:
      if(usedWireNumber < 4){
    err = -2;
    break;
      }
      //break;
      factor *= 0.75;
#ifdef TRKRECO_DEBUG_DETAIL 
      std::cout << "factor = " << factor << std::endl;
      std::cout << "chi2 = " << chi2 << std::endl;
#endif
      if(factor < 0.01)break;
    }
 
    chi2Old = chi2;
 
    //...Cal. helix parameters for next loop...
    f = solve(d2chi2d2a, dchi2da);
    da[0] = f[0];
    da[1] = f[1];
    da[2] = f[2];
    da[3] = 0.;
    da[4] = 0.;
    
    a -= factor*da;
    _helix->a(a);      
    ++nTrial;
  }
  if(err){
    return err;
  }
 
  //...Cal. error matrix...
  SymMatrix Ea(5,0);
  unsigned dim = 3;
  SymMatrix Eb = d2chi2d2a;
  SymMatrix Ec = Eb.inverse(err);
  Ea[0][0] = Ec[0][0];
  Ea[0][1] = Ec[0][1];
  Ea[0][2] = Ec[0][2];
  Ea[1][1] = Ec[1][1];
  Ea[1][2] = Ec[1][2];
  Ea[2][2] = Ec[2][2];
 
  //...Store information...
  if(!err){
    _helix->a(a);
    _helix->Ea(Ea);
  }else{
    err = -2;
  }
    
  _ndf  = usedWireNumber-dim;
  _chi2 = chi2;
    
  //_fitted = true;
 
#define JJJTEST 0
#if JJJTEST
  double gmaxX = -9999. ,gminX = 9999.;
  double gmaxY = -9999. ,gminY = 9999.;
  FILE *gnuplot, *data;
  double step = 200.;
  double dStep = 2.*M_PI/step;
  for(int i=0,size = _links.length();i<size;++i){
    TMLink * l = _links[i];
    double drift = l->hit()->distance(0);
    char name[100] = "dat";
    char counter[10] = "";
    sprintf(counter,"%02d",i);
    strcat(name,counter);
    if((data = fopen(name,"w")) != NULL){
      for(int ii=0;ii<step;++ii){
    double X = l->wire()->xyPosition().x() + drift*cos(dStep*static_cast<double>(ii));
    double Y = l->wire()->xyPosition().y() + drift*sin(dStep*static_cast<double>(ii));
    fprintf(data,"%lf, %lf\n",X,Y);
    if(gmaxX < X)gmaxX = X;
    if(gminX > X)gminX = X;
    if(gmaxY < Y)gmaxY = Y;
    if(gminY > Y)gminY = Y;
      }
      fclose(data);
    }
  }
  step = 300.;
  dStep = 2.*M_PI/step;
  if((data = fopen("datc","w")) != NULL){
    for(int ii=0;ii<step;++ii){
      double X = _helix->center().x() + _helix->radius()*cos(dStep*static_cast<double>(ii));
      double Y = _helix->center().y() + _helix->radius()*sin(dStep*static_cast<double>(ii));
      fprintf(data,"%lf, %lf\n",X,Y);
    }
    fclose(data);
  }
  if((gnuplot = popen("gnuplot","w")) != NULL){
    fprintf(gnuplot,"set size 0.721,1.0 \n");
    fprintf(gnuplot,"set xrange [%f:%f] \n",gminX,gmaxX);
    fprintf(gnuplot,"set yrange [%f:%f] \n",gminY,gmaxY);
    if(_links.length() == 4){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line \n");
    }else if(_links.length() == 5){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line \n");
    }else if(_links.length() == 6){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line \n");
    }else if(_links.length() == 7){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line, \"dat06\" with line \n");
    }else if(_links.length() == 8){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line, \"dat06\" with line, \"dat07\" with line \n");
    }else if(_links.length() == 9){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line, \"dat06\" with line, \"dat07\" with line, \"dat08\" with line \n");
    }else if(_links.length() == 10){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line, \"dat06\" with line, \"dat07\" with line, \"dat08\" with line, \"dat09\" with line \n");
    }else if(_links.length() >= 11){
      fprintf(gnuplot,"plot \"datc\" with line, \"dat00\" with line, \"dat01\" with line, \"dat02\" with line, \"dat03\" with line, \"dat04\" with line, \"dat05\" with line, \"dat06\" with line, \"dat07\" with line, \"dat08\" with line, \"dat09\" with line, \"dat10\" with line \n");
    }
    fflush(gnuplot);
    char tmp[8];
    gets(tmp);
    pclose(gnuplot);
  }
#endif //JJJTEST
 
  return err;
}

Referenced by refine2D().

fitting() [1/2]

unsigned TTrack::fitting ( unsigned  fitMask )

inline

Definition at line 636 of file TTrack.h.

                          {
    _state |= ((a << TrackFitShift) & TrackFitMask);
    return fitting();
}

fitting() [2/2]

unsigned TTrack::fitting ( void  ) const

inline

sets/returns fitting status.

Definition at line 630 of file TTrack.h.

                          {
    return (_state >> TrackFitShift) & TrackFitMask;
}

Referenced by fitting().

getFinderType()

unsigned TTrack::getFinderType ( void  ) const

inline

Definition at line 680 of file TTrack.h.

                                       {
        return _findertype;
}

HelCyl()

int TTrack::HelCyl	(	double	rhole,
		double	rcyl,
		double	zb,
		double	zf,
		double	epsl,
		double &	phi,
		HepPoint3D &	xp
	)		const

returns a cathode hit list.

calculates an intersection of this track and a cylinder.

Definition at line 316 of file TTrack.cxx.

                                   {
 
 
  int status(0);    // return value
  //---------------------------------------------------------------------
  //  value | ext |   status 
  //---------------------------------------------------------------------
  //     1. |  OK |
  //    -1. |  NO |  charge = 0
  //     0. |  NO |  | tanl | < 0.1  ( neglect | lamda | < 5.7 deg. )
  //        |     |  or | dPhi | > 2 pi ( neglect curly track )
  //        |     |  or cannot reach to r=rhole at z = zb or zf.
  //     2. |  OK |   backward , ext point set on z = zb
  //     3. |  OK |   forward  , ext point set on z = zf
  //---------------------------------------------------------------------
  // * when value = 0,2,3 , ext(z) <= zb  or ext(z) >= zf
 
  //--- debug
  //  std::cout << "   "  << std::endl;
  //  std::cout << "HelCyl called ..  rhole=" << rhole << " rCyl=" << rCyl ; 
  //  std::cout << " zb=" << zb << " zf=" << zf << " epsl=" << epsl << std::endl;
  //--- debug end
 
  // Check of Charge
 
  if ( int(_charge) == 0 ) {
    std::cout << "HelCyl gets a straight line !!" << std::endl;
    return -1 ;
  }
 
  // parameters
 
  HepPoint3D CenterCyl( 0., 0., 0. );
  HepPoint3D CenterTrk = _helix->center();
  double  rTrk = fabs( _helix->radius() );
 
  double rPivot = fabs( _helix->pivot().perp() );
 
  double phi0 = _helix->a()[1];
  double tanl = _helix->a()[4];
  //  double zdz = _helix->pivot().z() + _helix->a()[3];
  double dPhi;
  double zee;
 
 
  // Calculate intersections between cylinder and track
  // if return value = 2 track hitting barrel part
 
  HepPoint3D Cross1, Cross2;
 
  if (intersection( CenterTrk,_charge * rTrk ,CenterCyl,rCyl,
           epsl, 
           Cross1,Cross2)
      == 2 ) {
 
    double phiCyl = atan2( _charge * ( CenterTrk.y() - Cross1.y() ),
               _charge * ( CenterTrk.x() - Cross1.x() ) );
    phiCyl = ( phiCyl > 0. ) ? phiCyl : phiCyl + 2. * M_PI;
 
    dPhi = phiCyl - phi0;
 
    // dPhi region ( at cylinder )
    //      -pi <= dPhi < pi
 
 
    double PhiYobun = 1./ fabs( _helix->radius() );
    double zero = 0.00001;
 
    if( _charge >=0. ){
      if( dPhi > PhiYobun ) dPhi -= 2. * M_PI;
      if( -2. * M_PI  - zero <= dPhi <= ( -2.* M_PI + PhiYobun ) ) 
    dPhi += 2. * M_PI;
    }
 
    if( _charge < 0. ){
      if( dPhi < -PhiYobun ) dPhi += 2. * M_PI;
      if( 2. * M_PI + zero >= dPhi >= ( 2. * M_PI - PhiYobun ) )
    dPhi -= 2. * M_PI;
    }
 
    if( dPhi < - M_PI ) dPhi += 2. * M_PI;
    if( dPhi >= M_PI ) dPhi -= 2. * M_PI;
 
      //--debug 
      //  std::cout << "dPhi = " << dPhi << std::endl;
      //--debug end 
 
    xp.setX( Cross1.x() );
    xp.setY( Cross1.y() );
    xp.setZ( _helix->x(dPhi).z() );
    //    xp.setZ( zdz - _charge * rTrk * tanl * dPhi );
 
    if ( xp.z() > zb && xp.z() < zf ) {
      phi = dPhi;
 //--- debug ---
 //      std::cout << "return1 ( ext success )" <<  std::endl;
 //      std::cout << " xp:= " << xp.x() << ", " << xp.y() << ", " << xp.z() << std::endl;
  //--- debug  ----
      return 1 ;
    }
  }
 
 
  // tracks hitting endcaps
 
  if ( fabs(tanl) < 0.1 ) {
    //--- debug ---
    //    std::cout << "return0 ( ext failed , |tanl| < 0.1 )" <<  std::endl;
    //    std::cout << " xp:= " << xp.x() << ", " << xp.y() << ", " << xp.z() << std::endl;
    //--- debug ---
    return 0;
  }
 
  if ( tanl > 0. ) {
    zee = zf ;
    status = 3 ;
  } 
  else {
    zee = zb ;
    status = 2 ;
  }
 
  dPhi = _charge * ( _helix->x(0.).z() - zee )/rTrk/tanl;
  //  dPhi = _charge * ( zdz - zee )/rTrk/tanl;
 
  // Requre dPhi < 2*pi 
 
  if ( fabs(dPhi) > 2. * M_PI ) {
    //--- debug ---
    //    std::cout << " return0 ( ext failed , dPhi > 2pi )" << std::endl;
    //    std::cout << " xp:= " << xp.x() << ", " << xp.y() << ", " << xp.z() << std::endl;
    //--- debug  ---
    return 0 ;
  }
 
  xp.setX( _helix->x(dPhi).x() );
  xp.setY( _helix->x(dPhi).y() );
  xp.setZ( zee );
 
  double test = xp.perp2()  - rhole*rhole ;
  if ( test < 0. )  {
    //--- debug ---
    //    std::cout << "return0 ( cannot reach to rhole at z=edge )" << std::endl;
    //    std::cout << " xp:= " << xp.x() << ", " << xp.y() << ", " << xp.z() << std::endl;
    //--- debug ---
    return 0 ;
  }
 
  phi = dPhi ;
  //--- debug ---
  //  std::cout << "return" << status << std::endl;
  //  std::cout << " xp:= " << xp.x() << ", " << xp.y() << ", " << xp.z() << std::endl;
  //--- debug ---
 
  return status ;
 
}

helix()

const Helix & TTrack::helix ( void  ) const

inline

returns helix parameter.

Definition at line 477 of file TTrack.h.

                        {
#ifdef TRKRECO_DEBUG
    // if (! _fitted) std::cout << "TTrack::helix !!! helix not updated" << std::endl;
#endif
    return * _helix;
}

Referenced by TBuilder::buildStereo(), TBuilder0::buildStereo(), TBuilderCosmic::buildStereo(), TBuilderCurl::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), calVirtualCircle(), and TTrackManager::merge().

impact()

double TTrack::impact ( void  ) const

inline

returns signed impact parameter to the origin.

Definition at line 571 of file TTrack.h.

                         {
    return fabs(_helix->radius()) - _helix->center().mag();
}

Referenced by TBuilderCurl::buildStereo(), TBuilderCurl::buildStereoMC(), and TMSelector::select().

mother() [1/2]

TTrack * TTrack::mother ( TTrack *  a )

inline

Definition at line 649 of file TTrack.h.

                         {
    if (a) _state |= (TrackHasMother << TrackRelationShift);
    else   _state &= (~(TrackHasMother << TrackRelationShift));
    return _mother = a;
}

mother() [2/2]

TTrack * TTrack::mother ( void  ) const

inline

sets/returns mother/daughter.

Definition at line 643 of file TTrack.h.

                         {
    return _mother;
}

movePivot()

void TTrack::movePivot

(

void

)

moves pivot to the inner most hit.

Definition at line 262 of file TTrack.cxx.

                      {
     unsigned n = _links.length();
    if (! n) {
    std::cout << "TTrack::movePivot !!! can't move a pivot"
          << " because of no link";
    std::cout << std::endl;
    return;
    }
 
    //...Check cores...
     const AList<TMLink> & cores = TTrackBase::cores();
    const AList<TMLink> * links = & cores;
    if (cores.length() == 0) links = & _links;
 
    //...Hit loop...
    unsigned innerMost = 0;
    unsigned innerMostLayer = (* links)[0]->wire()->layerId();
    n = links->length();
    for (unsigned i = 1; i < n; i++) {
    TMLink * l = (* links)[i];
    if (l->wire()->layerId() < innerMostLayer) {
        innerMost = i;
        innerMostLayer = l->wire()->layerId();
    }
    }
 
    //...Move pivot...
    HepPoint3D newPivot = (* links)[innerMost]->positionOnWire();
    if (quality() & TrackQuality2D)
    newPivot.setZ(0.);
    _helix->pivot(newPivot);
}

name() [1/2]

const std::string & TTrack::name ( const std::string &  newName )

inline

Definition at line 522 of file TTrack.h.

                                {
    return _name = a;
}

name() [2/2]

const std::string & TTrack::name ( void  ) const

inline

returns/sets name.

Definition at line 516 of file TTrack.h.

                       {
    return _name;
}

Referenced by TTrackMC::dump(), fit2D(), and TTrackManager::salvageAssociateHits().

ndf()

unsigned TTrack::ndf ( void  ) const

inline

returns NDF.

Definition at line 486 of file TTrack.h.

                      {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TTrack::ndf !!! ndf not updated" << std::endl;
#endif
    return _ndf;
}

Referenced by TTrackManager::sortTracksByQuality().

objectType()

unsigned TTrack::objectType ( void  ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 559 of file TTrack.h.

                             {
    return Track;
}

p()

Hep3Vector TTrack::p ( void  ) const

inline

returns momentum.

Definition at line 553 of file TTrack.h.

                    {
    return _helix->momentum(0.);
}

Referenced by dump(), stereoHitForCurl(), and szPosition().

pt()

double TTrack::pt ( void  ) const

inline

returns Pt.

Definition at line 528 of file TTrack.h.

                     {
    return 1. / fabs(_helix->a()[2]);
}

Referenced by TBuilderCurl::buildStereo(), TBuilderCurl::buildStereoMC(), dump(), TTrackManager::maskNormal(), TTrackManager::merge(), TMSelector::select(), and TTrackManager::sortTracksByPt().

ptot()

double TTrack::ptot ( void  ) const

inline

returns magnitude of momentum.

Definition at line 540 of file TTrack.h.

                       {
    return (1. / fabs(_helix->a()[2])) * 
    sqrt(1. + _helix->a()[4] * _helix->a()[4]);
}

pz()

double TTrack::pz ( void  ) const

inline

returns Pz.

Definition at line 534 of file TTrack.h.

                     {
    return (1. / fabs(_helix->a()[2])) * _helix->a()[4];
}

Referenced by TBuilderCurl::buildStereo().

quality() [1/2]

unsigned TTrack::quality ( unsigned  qualityMask )

inline

Definition at line 620 of file TTrack.h.

                          {
//      _state = ((a << TrackQualityShift) & TrackQualityMask) |
//      (_state & (~ (TrackQualityMask << TrackQualityShift)));
    _state = ((a & TrackQualityMask) << TrackQualityShift) |
    (_state & (~ (TrackQualityMask << TrackQualityShift)));
    return quality();
}

quality() [2/2]

unsigned TTrack::quality ( void  ) const

inline

sets/returns quality.

Definition at line 614 of file TTrack.h.

                          {
    return (_state >> TrackQualityShift) & TrackQualityMask;
}

Referenced by movePivot(), and quality().

radius()

double TTrack::radius ( void  ) const

inline

returns signed radius.

Definition at line 577 of file TTrack.h.

                         {
    return _helix->radius();
}

Referenced by calVirtualCircle().

refine2D()

void TTrack::refine2D	(	AList< TMLink > &	list,
		float	maxSigma
	)

fits itself with cathode hits.

remove bad points by chi2. Bad points are returned in a 'list'. fit() should be called before calling this function. (using stereo wire as axial wire(z=0))

Definition at line 301 of file TTrack.cxx.

                                                     {
    unsigned n = _links.length();
    AList<TMLink> bad;
    for (unsigned i = 0; i < n; ++i) {
        if (_links[i]->pull() > maxSigma) bad.append(_links[i]);
    }
    _links.remove(bad);
    if (bad.length()){
      _fitted = false;
      fit2D();
    }
    list.append(bad);
}

segments() [1/2]

AList< TSegment > & TTrack::segments ( void  )

inline

returns AList<TSegment>.

Definition at line 583 of file TTrack.h.

                     {
    return _segments;
}

Referenced by TBuilder::buildStereoNew().

segments() [2/2]

const AList< TSegment > & TTrack::segments ( void  ) const

inline

Definition at line 589 of file TTrack.h.

                           {
    return _segments;
}

setFinderType()

void TTrack::setFinderType ( unsigned  type )

inline

Definition at line 676 of file TTrack.h.

                                       {
        _findertype=type;
}

state()

unsigned TTrack::state ( void  ) const

inline

returns/sets internal state.(for bank output)

Definition at line 671 of file TTrack.h.

                        {
    return _state;
}

stereoHitForCurl() [1/4]

int TTrack::stereoHitForCurl

(

AList< TMLink > &

list

)

const

calculates arc length and z for a stereo hit. uses these functions for curl tracks(included svd version).

Definition at line 2009 of file TTrack.cxx.

{
  if(list.length() == 0)return -1;
 
  HepPoint3D center = _helix->center();
  HepPoint3D tmp(-999., -999., 0.);
  double r  = fabs(_helix->curv());
  for(unsigned i = 0, size = list.length(); i < size; ++i){
    TMDCWireHit &h = *const_cast<TMDCWireHit*>(list[i]->hit());
    HepVector3D X = 0.5*(h.wire()->forwardPosition() +
              h.wire()->backwardPosition());
    HepVector3D x     = HepVector3D(X.x(), X.y(), 0.);
    HepVector3D w     = x - center;
    HepVector3D V     = h.wire()->direction();
    HepVector3D v     = HepVector3D(V.x(), V.y(), 0.);
    double   vmag2 = v.mag2();
    double   vmag  = sqrt(vmag2);
    //...temporary
    for(unsigned j = 0; j < 4; ++j)
      list[i]->arcZ(tmp,j);
    
    //...stereo?
    if (vmag == 0.) continue;
   
    double drift = h.drift(WireHitLeft);
    double R[2] = {r + drift, r - drift};
    double wv = w.dot(v);
    double d2[2];
    d2[0] = vmag2*R[0]*R[0] + (wv*wv - vmag2*w.mag2()); //...= v^2*(r^2 - w^2*sin()^2)...outer
    d2[1] = vmag2*R[1]*R[1] + (wv*wv - vmag2*w.mag2()); //...= v^2*(r^2 - w^2*sin()^2)...inner
    
    //...No crossing in R/Phi plane...
    if (d2[0] < 0. && d2[1] < 0.) continue;
 
    bool ok_inner(true);
    bool ok_outer(true);
    double d[2] = {-1., -1.};
    //...outer
    if(d2[0] >= 0.){
      d[0] = sqrt(d2[0]);
    }else{
      ok_outer = false;
    }
    if(d2[1] >= 0.){
      d[1] = sqrt(d2[1]);
    }else{
      ok_inner = false;
    }
    
    //...Cal. length and z to crossing points...
    double l[2][2];
    double z[2][2];
    //...outer
    if(ok_outer){
      l[0][0] = (- wv + d[0]) / vmag2; //...= (-wvcos()+d)/v/v = (-wcos() + (r^2 - w^2*sin()^2)^0.5)/v
      l[1][0] = (- wv - d[0]) / vmag2; //...= (-wvcos()+d)/v/v = (-wcos() - (r^2 - w^2*sin()^2)^0.5)/v
      z[0][0] = X.z() + l[0][0]*V.z();
      z[1][0] = X.z() + l[1][0]*V.z();
    }
    //...inner
    if(ok_inner){
      l[0][1] = (- wv + d[1]) / vmag2; //...= (-wvcos()+d)/v/v = (-wcos() + (r^2 - w^2*sin()^2)^0.5)/v
      l[1][1] = (- wv - d[1]) / vmag2; //...= (-wvcos()+d)/v/v = (-wcos() - (r^2 - w^2*sin()^2)^0.5)/v
      z[0][1] = X.z() + l[0][1]*V.z();
      z[1][1] = X.z() + l[1][1]*V.z();
    }
    
    //...Cal. xy position of crossing points...
    HepVector3D p[2][2];
#if 1
    HepVector3D tp[2][2];
#endif
    if(ok_outer){
      p[0][0] = x + l[0][0] * v;
      p[1][0] = x + l[1][0] * v;
#if 1
      tp[0][0] = p[0][0];
      tp[1][0] = p[1][0];
#endif
      HepVector3D tmp_pc = p[0][0] - center;
      HepVector3D pc0 = HepVector3D(tmp_pc.x(), tmp_pc.y(), 0.);
      p[0][0] -= drift/pc0.mag()*pc0;
      tmp_pc = p[1][0] - center;
      HepVector3D pc1 = HepVector3D(tmp_pc.x(), tmp_pc.y(), 0.);
      p[1][0] -= drift/pc1.mag()*pc1;
    }
#define JJTEST 0
    if(ok_inner){
      p[0][1] = x + l[0][1] * v;
      p[1][1] = x + l[1][1] * v;
#if JJTEST
      tp[0][1] = p[0][1];
      tp[1][1] = p[1][1];
#endif
      HepVector3D tmp_pc = p[0][1] - center;
      HepVector3D pc0 = HepVector3D(tmp_pc.x(), tmp_pc.y(), 0.);
      p[0][1] += drift/pc0.mag()*pc0;
      tmp_pc = p[1][1] - center;
      HepVector3D pc1 = HepVector3D(tmp_pc.x(), tmp_pc.y(), 0.);
      p[1][1] += drift/pc1.mag()*pc1;
    }
 
    //...boolean
    bool ok_xy[2][2];
    if(ok_outer){
      ok_xy[0][0] = true;
      ok_xy[1][0] = true;
    }else{
      ok_xy[0][0] = false;
      ok_xy[1][0] = false;
    }
    if(ok_inner){
      ok_xy[0][1] = true;
      ok_xy[1][1] = true;
    }else{
      ok_xy[0][1] = false;
      ok_xy[1][1] = false;
    }
    if(ok_outer){
      if (_charge * (center.x() * p[0][0].y() - center.y() * p[0][0].x())  < 0.)
    ok_xy[0][0] = false;
      if (_charge * (center.x() * p[1][0].y() - center.y() * p[1][0].x())  < 0.)
    ok_xy[1][0] = false;
    }
    if(ok_inner){
      if (_charge * (center.x() * p[0][1].y() - center.y() * p[0][1].x())  < 0.)
    ok_xy[0][1] = false;
      if (_charge * (center.x() * p[1][1].y() - center.y() * p[1][1].x())  < 0.)
    ok_xy[1][1] = false;
    }
    if(!ok_inner && ok_outer && (!ok_xy[0][0]) && (!ok_xy[1][0])){
      continue;
    }
    if(ok_inner && !ok_outer && (!ok_xy[0][1]) && (!ok_xy[1][1])){
      continue;
    }
    
    //...Check z position...
    if(ok_xy[0][0]){
      if (z[0][0] < h.wire()->backwardPosition().z() || 
      z[0][0] > h.wire()->forwardPosition().z()) ok_xy[0][0] = false;
    }
    if(ok_xy[1][0]){
      if (z[1][0] < h.wire()->backwardPosition().z() || 
      z[1][0] > h.wire()->forwardPosition().z()) ok_xy[1][0] = false;
    }
    if(ok_xy[0][1]){
      if (z[0][1] < h.wire()->backwardPosition().z() || 
      z[0][1] > h.wire()->forwardPosition().z()) ok_xy[0][1] = false;
    }
    if(ok_xy[1][1]){
      if (z[1][1] < h.wire()->backwardPosition().z() || 
      z[1][1] > h.wire()->forwardPosition().z()) ok_xy[1][1] = false;
    }
    if ((!ok_xy[0][0]) && (!ok_xy[1][0]) &&
    (!ok_xy[0][1]) && (!ok_xy[1][1])){
      continue;
    }
    double cosdPhi, dPhi;
    unsigned index;
    index = 0;
    if(ok_xy[0][0]){
      //...cal. arc length...
      cosdPhi = - center.dot((p[0][0] - center).unit()) / center.mag();
      if(fabs(cosdPhi) < 1.0){
    dPhi = acos(cosdPhi);
      }else if(cosdPhi >= 1.0){
    dPhi = 0.;
      }else{
    dPhi = M_PI;
      }
      list[i]->arcZ(HepPoint3D(r*dPhi, z[0][0], 0.), index);
      //std::cout << r*dPhi << ", " << z[0][0] << std::endl;
      ++index;
    }
    if(ok_xy[1][0]){
      //...cal. arc length...
      cosdPhi = - center.dot((p[1][0] - center).unit()) / center.mag();
      if(fabs(cosdPhi) < 1.0){
    dPhi = acos(cosdPhi);
      }else if(cosdPhi >= 1.0){
    dPhi = 0.;
      }else{
    dPhi = M_PI;
      }
      list[i]->arcZ(HepPoint3D(r*dPhi, z[1][0], 0.), index);
      //std::cout << r*dPhi << ", " << z[1][0] << std::endl;
      ++index;
    }
    if(ok_xy[0][1]){
      //...cal. arc length...
      cosdPhi = - center.dot((p[0][1] - center).unit()) / center.mag();
      if(fabs(cosdPhi) < 1.0){
    dPhi = acos(cosdPhi);
      }else if(cosdPhi >= 1.0){
    dPhi = 0.;
      }else{
    dPhi = M_PI;
      }
      list[i]->arcZ(HepPoint3D(r*dPhi, z[0][1], 0.), index);
      //std::cout << r*dPhi << ", " << z[0][1] << std::endl;
      ++index;
    }
    if(ok_xy[1][1]){
      //...cal. arc length...
      cosdPhi = - center.dot((p[1][1] - center).unit()) / center.mag();
      if(fabs(cosdPhi) < 1.0){
    dPhi = acos(cosdPhi);
      }else if(cosdPhi >= 1.0){
    dPhi = 0.;
      }else{
    dPhi = M_PI;
      }
      list[i]->arcZ(HepPoint3D(r*dPhi, z[1][1], 0.), index);
      //std::cout << r*dPhi << ", " << z[1][1] << std::endl;
      ++index;
    }
 
#if JJTEST
    double gmaxX = -9999. ,gminX = 9999.;
    double gmaxY = -9999. ,gminY = 9999.;
    FILE *gnuplot, *data;
    double step = 100.;
    double dStep = 2.*M_PI/step;
    if((data = fopen("dat1","w")) != NULL){
      if(ok_xy[0][0]){
    for(int ii=0;ii<step;++ii){
      double X = tp[0][0].x() + drift*cos(dStep*static_cast<double>(ii));
      double Y = tp[0][0].y() + drift*sin(dStep*static_cast<double>(ii));
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
    }
      }
      fclose(data);
    }
    if((data = fopen("dat2","w")) != NULL){
      if(ok_xy[1][0]){
    for(int ii=0;ii<step;++ii){
      double X = tp[1][0].x() + drift*cos(dStep*static_cast<double>(ii));
      double Y = tp[1][0].y() + drift*sin(dStep*static_cast<double>(ii));
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
    }
      }
      fclose(data);
    }
    if((data = fopen("dat3","w")) != NULL){
      if(ok_xy[0][1]){
    for(int ii=0;ii<step;++ii){
      double X = tp[0][1].x() + drift*cos(dStep*static_cast<double>(ii));
      double Y = tp[0][1].y() + drift*sin(dStep*static_cast<double>(ii));
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
    }
      }
      fclose(data);
    }
    if((data = fopen("dat4","w")) != NULL){
      if(ok_xy[1][1]){
    for(int ii=0;ii<step;++ii){
      double X = tp[1][1].x() + drift*cos(dStep*static_cast<double>(ii));
      double Y = tp[1][1].y() + drift*sin(dStep*static_cast<double>(ii));
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
    }
      }
      fclose(data);
    }
    HepVector3D tX = h.wire()->forwardPosition()-h.wire()->backwardPosition();
    HepVector3D tDist(tX.x(), tX.y(), 0.);
    double tD = tDist.mag();
    double vvvM = 1./ v.mag();
    HepVector3D tDire = vvvM*v;
    step = 2.;
    dStep = tD/step;
    if((data = fopen("dat5","w")) != NULL){
      for(int ii=0;ii<step+1;++ii){
    double X = h.wire()->backwardPosition().x()+dStep*static_cast<double>(ii)*tDire.x();
    double Y = h.wire()->backwardPosition().y()+dStep*static_cast<double>(ii)*tDire.y();
    fprintf(data,"%lf, %lf\n",X,Y);
    if(gmaxX < X)gmaxX = X;
    if(gminX > X)gminX = X;
    if(gmaxY < Y)gmaxY = Y;
    if(gminY > Y)gminY = Y;
      }
      fclose(data);
    }
    if((data = fopen("dat6","w")) != NULL){
      double X = h.wire()->backwardPosition().x();
      double Y = h.wire()->backwardPosition().y();
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
      fclose(data);    
    }
    if((data = fopen("dat7","w")) != NULL){
      double X = h.wire()->forwardPosition().x();
      double Y = h.wire()->forwardPosition().y();
      fprintf(data,"%lf, %lf\n",X,Y);
      if(gmaxX < X)gmaxX = X;
      if(gminX > X)gminX = X;
      if(gmaxY < Y)gmaxY = Y;
      if(gminY > Y)gminY = Y;
      fclose(data);  
    }
    step = 300.;
    dStep = 2.*M_PI/step;
    if((data = fopen("dat8","w")) != NULL){
      for(int ii=0;ii<step;++ii){
    double X = center.x() + r*cos(dStep*static_cast<double>(ii));
    double Y = center.y() + r*sin(dStep*static_cast<double>(ii));
    fprintf(data,"%lf, %lf\n",X,Y);
      }
      fclose(data);
    }
    if((data = fopen("dat9","w")) != NULL){
      if(ok_xy[0][0]){
    double X = p[0][0].x();
    double Y = p[0][0].y();
    fprintf(data,"%lf, %lf\n",X,Y);
      }
      fclose(data);  
    }
    if((data = fopen("dat10","w")) != NULL){
      if(ok_xy[1][0]){
    double X = p[1][0].x();
    double Y = p[1][0].y();
    fprintf(data,"%lf, %lf\n",X,Y);
      }
      fclose(data);  
    }
    if((data = fopen("dat11","w")) != NULL){
      if(ok_xy[0][1]){
    double X = p[0][1].x();
    double Y = p[0][1].y();
    fprintf(data,"%lf, %lf\n",X,Y);
      }
      fclose(data);  
    }
    if((data = fopen("dat12","w")) != NULL){
      if(ok_xy[1][1]){
    double X = p[1][1].x();
    double Y = p[1][1].y();
    fprintf(data,"%lf, %lf\n",X,Y);
      }
      fclose(data);  
    }
    std::cout << "Drift Distance = " << drift << ", xy1cm -> z" << V.z()/v.mag() << "cm, xyDist(cm) = " << tD << std::endl;
    if(tX.z()<0.)std::cout << "ERROR : F < B" << std::endl;
    if((gnuplot = popen("gnuplot","w")) != NULL){
      fprintf(gnuplot,"set size 0.721,1.0 \n");
      fprintf(gnuplot,"set xrange [%f:%f] \n",gminX,gmaxX);
      fprintf(gnuplot,"set yrange [%f:%f] \n",gminY,gmaxY);
      fprintf(gnuplot,"plot \"dat1\" with line, \"dat2\" with line, \"dat3\" with line, \"dat4\" with line, \"dat5\" with line, \"dat6\", \"dat7\", \"dat8\" with line, \"dat9\", \"dat10\", \"dat11\", \"dat12\" \n");
      fflush(gnuplot);
      char tmp[8];
      gets(tmp);
      pclose(gnuplot);
    }
#endif
  }
  return 0;
}

stereoHitForCurl() [2/4]

int TTrack::stereoHitForCurl	(	TMLink &	link,
		AList< HepPoint3D > &	arcZList
	)		const

stereoHitForCurl() [3/4]

int TTrack::stereoHitForCurl	(	TMLink &	link,
		TMLink &	link1
	)		const

stereoHitForCurl() [4/4]

int TTrack::stereoHitForCurl	(	TMLink &	link,
		TMLink &	link1,
		TMLink &	link2
	)		const

szPosition() [1/3]

int TTrack::szPosition	(	const HepPoint3D &	p,
		HepPoint3D &	szPosition
	)		const

calculates arc length for a point.

Definition at line 3608 of file TTrack.cxx.

                                                              {
 
    //...Cal. arc length...
    HepPoint3D center = _helix->center();
    HepPoint3D xy = p;
    xy.setZ(0.);
    double cosdPhi = - center.dot((xy - center).unit()) / center.mag();
    double dPhi;
    if (fabs(cosdPhi) <= 1.0) {
    dPhi = acos(cosdPhi);
    }
    else if (cosdPhi > 1.0) {
    dPhi = 0.0;
    }
    else {
    dPhi = M_PI;
    }
 
    //...Finish...
    sz.setX(_helix->curv() * dPhi);
    sz.setY(p.z());
    sz.setZ(0.);
 
    return 0;
}

szPosition() [2/3]

int TTrack::szPosition	(	const TSegment &	segment,
		TMLink &	link
	)		const

calculates arc length and z for a segment. Results are stored in TMLink.

Definition at line 3570 of file TTrack.cxx.

                                                             {
 
    //...Pick up a wire which represents segment position...
    AList<TMLink> links = segment.cores();
    unsigned n = links.length();
//    std::cout<<" szPosition TTrack:: segment.cores().length():"<<n<<endl;
//    for (unsigned i = 1; i < n; i++) {  
//   std::cout<<"drift distance"<<links[i]->hit()->drift()<<endl;
//    }  
    if (! n) return -1;
    TMLink * minL = links[0];
    float minDist = links[0]->drift();
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"minDist szPosition TTrack:"<<minDist<<endl;
#endif
    for (unsigned i = 1; i < n; i++) {
    if (links[i]->hit()->drift() < minDist) {
        minDist = links[i]->drift();
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout<<"minDist szPosition TTrack:"<<minDist<<endl;   
#endif
        minL = links[i];
    }
    }
 
    //...sz calculation...
    a.position(minL->position());
    a.leftRight(2);
    a.hit(minL->hit());
    int err = szPosition(a);
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"err of szPosition TTrack:"<<err<<endl;      
#endif
    if (err) return -2;
    return 0;
}

szPosition() [3/3]

int TTrack::szPosition

(

TMLink &

link

)

const

calculates arc length and z for a stereo hit.

Definition at line 3372 of file TTrack.cxx.

                                      {
    const TMDCWireHit & h = * link.hit();
    HepVector3D X = 0.5 * (h.wire()->forwardPosition()
            + h.wire()->backwardPosition());
    //    double theta = atan2(X.y(), X.x());
    //    HepVector3D lr(h.distance(WireHitLeft) * sin(theta),
    //      - h.distance(WireHitLeft) * cos(theta),
    //      0.);
 
    HepVector3D xx = HepVector3D(X.x(), X.y(), 0.);
    HepPoint3D center = _helix->center();
    HepVector3D yy = center - xx;
    HepVector3D ww = HepVector3D(yy.x(), yy.y(), 0.);
    double wwmag2 = ww.mag2();
    double wwmag = sqrt(wwmag2);
    HepVector3D lr(h.drift(WireHitLeft)/wwmag * ww.x(),
        h.drift(WireHitLeft)/wwmag * ww.y(),
        0.);
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"old lr "<<lr<<endl;
    std::cout<<"old X "<<X<<endl;
    std::cout<<"link.leftRight "<<link.leftRight()<<endl;
#endif
    //...Check left or right...
    //   // change to bes3 .. test..
    //if (link.leftRight() == WireHitRight) {
    //lr = - lr;
    //}
    //if (link.leftRight() == WireHitLeft) lr = - lr; 
    //else if (link.leftRight() == 2) lr = ORIGIN;
 
    //yzhang 2012-05-07 
    const double Bz = -1000*m_pmgnIMF->getReferField();
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"charge "<<_charge<<" Bz "<<Bz<<endl;
#endif
    if (_charge*Bz > 0){ //yzhang 2012-05-07 
      if (link.leftRight() == WireHitRight){
        lr = - lr;//right
      }
    }else{
      if (link.leftRight() == WireHitLeft){
        lr = - lr;//left
      }
    }
    if (link.leftRight() == 2) lr = ORIGIN;
    //zhangy
 
    X += lr;
 
    //...Prepare vectors...
    //    HepPoint3D center = _helix->center();
    HepPoint3D tmp(-9999., -9999., 0.);
    HepVector3D x = HepVector3D(X.x(), X.y(), 0.);
    HepVector3D w = x - center;
    // //modified the next sentence because the direction are different from belle.
    HepVector3D V = h.wire()->direction();
    //    //    to bes3
    // //    HepVector3D V = - h.wire()->direction();   
 
    HepVector3D v = HepVector3D(V.x(), V.y(), 0.);
    double vmag2 = v.mag2();
    double vmag = sqrt(vmag2);
    
    double r = _helix->curv();
    double wv = w.dot(v);
    //    //zsl for bes3
    //    wv = abs(wv);
    double d2 = wv * wv - vmag2 * (w.mag2() - r * r);
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"lr "<<lr<<endl;
    std::cout<<"forwardPosition "<<h.wire()->forwardPosition()<<endl;
    std::cout<<"backwardPosition "<<h.wire()->backwardPosition()<<endl;
    std::cout<<"X "<<X<<endl;
    std::cout<<"center "<<center<<endl;
    std::cout<<"xx "<<xx<<endl;
    std::cout<<"ww "<<ww<<endl;
    std::cout<<"TTrack::wire direction:"<<h.wire()->direction()<<endl;
    std::cout<<"x "<<x<<endl;
    std::cout<<"w "<<w<<endl;
    std::cout<<"sz,Track::vmag:"<<vmag<<", helix_r:"<<r<<", wv:"<<wv<<", d:"<<sqrt(d2)<<endl;     
#endif
 
     //...No crossing in R/Phi plane... This is too tight...
 
    if (d2 < 0.) {
    link.position(tmp);
 
#ifdef TRKRECO_DEBUG
    std::cout << "TTrack !!! stereo: 0. > d2 = " << d2 << " "
          << link.leftRight() << std::endl;
#endif
    return -1;
    }
    double d = sqrt(d2);
 
    //...Cal. length to crossing points...
    double l[2];
    l[0] = (- wv + d) / vmag2;
    l[1] = (- wv - d) / vmag2;
 
    //...Cal. z of crossing points...
    bool ok[2];
    ok[0] = true;
    ok[1] = true;
    double z[2];
    z[0] = X.z() + l[0] * V.z();
    z[1] = X.z() + l[1] * V.z();
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout<<"X.z():"<<X.z()<<endl;
    std::cout<<"szPosition::z(0) "<<z[0]<<" z(1)"<<z[1]<<" leftRight "<<link.leftRight()<<endl;
    std::cout<<"szPosition::wire backwardPosition and forwardPosition:"<< h.wire()->backwardPosition().z()<<","<<h.wire()->forwardPosition().z()<<endl;
    std::cout << "    l0, l1 = " << l[0] << ", " << l[1] << std::endl;
    std::cout << "    z0, z1 = " << z[0] << ", " << z[1] << std::endl;
    std::cout << "    backward = " << h.wire()->backwardPosition().z() << std::endl;
    std::cout << "    forward = " << h.wire()->forwardPosition().z() << std::endl;
#endif
 
    //...Check z position... //yzhang change 2012-05-03 
    //modified because Belle backward and forward are different from BESIII
    
    /*
    if (link.leftRight() == 2) {
    if (z[0] > h.wire()->backwardPosition().z()+20. 
        || z[0] < h.wire()->forwardPosition().z()-20.) ok[0] = false;
    if (z[1] > h.wire()->backwardPosition().z()+20.
        || z[1] < h.wire()->forwardPosition().z()-20.) ok[1] = false;
    }
    else {
    if (z[0] > h.wire()->backwardPosition().z()
        || z[0] < h.wire()->forwardPosition().z() ) ok[0] = false;
    if (z[1] > h.wire()->backwardPosition().z() 
        || z[1] < h.wire()->forwardPosition().z() ) ok[1] = false;
    }
    if ((! ok[0]) && (! ok[1])) {
    link.position(tmp);
    return -2;
    }*/
    //belle...
    if (link.leftRight() == 2) {
        if (z[0] < h.wire()->backwardPosition().z() - 20.
            || z[0] > h.wire()->forwardPosition().z() + 20.) ok[0] = false;
        if (z[1] < h.wire()->backwardPosition().z()-20.
            || z[1] > h.wire()->forwardPosition().z()+20.) ok[1] = false;
    }
    else {
        if (z[0] < h.wire()->backwardPosition().z()
            || z[0] > h.wire()->forwardPosition().z()) ok[0] = false;
        if (z[1] < h.wire()->backwardPosition().z()
            || z[1] > h.wire()->forwardPosition().z()) ok[1] = false;
    }
    if ((! ok[0]) && (! ok[1])) {
        link.position(tmp);
        return -2;
    }
    
 
    //...Cal. xy position of crossing points...
    HepVector3D p[2];
    p[0] = x + l[0] * v;
    p[1] = x + l[1] * v;
/*    if (_charge * (center.x() * p[0].y() - center.y() * p[0].x())  < 0.)  //liuqg, cosmic...
    ok[0] = false;
    if (_charge * (center.x() * p[1].y() - center.y() * p[1].x())  < 0.)
    ok[1] = false;
    if ((! ok[0]) && (! ok[1])){
      //        double tmp1 = _charge * (center.x() * p[0].y() - center.y() * p[0].x());
      //        double tmp2 = _charge * (center.x() * p[1].y() - center.y() * p[1].x()) ;
      //        if (link.leftRight() == 2) std::cout<<tmp1<<" "<<tmp2<<std::endl;
    link.position(tmp);
    return -3;
    }
*/
    //...Which one is the best?... Study needed...
    unsigned best = 0;
    if (ok[1]) best = 1;
 
    //...Cal. arc length...
    double cosdPhi = - center.dot((p[best] - center).unit()) / center.mag();
    double dPhi;
    if(fabs(cosdPhi)<=1.0) {
    dPhi = acos(cosdPhi);
    } else if (cosdPhi>1.0) {
    dPhi = 0.0;
    } else {
    dPhi = M_PI;
    }
 
    //...Finish...
    tmp.setX(r * dPhi);
    tmp.setY(z[best]);
    link.position(tmp);
 
    return 0;
}

Referenced by TBuilder::buildStereo(), TBuilder0::buildStereo(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), and szPosition().

type()

unsigned TTrack::type ( void  ) const

inlinevirtual

returns type. Definition is depending on an object type.

Reimplemented from TTrackBase.

Definition at line 565 of file TTrack.h.

                       {
    return defineType();
}

Referenced by TTrackManager::setCurlerFlags(), and setFinderType().

Friends And Related Function Documentation

Refit

friend class Refit

friend

Definition at line 373 of file TTrack.h.

TBuilder

friend class TBuilder

friend

Definition at line 374 of file TTrack.h.

TBuilder0

friend class TBuilder0

friend

Definition at line 375 of file TTrack.h.

TBuilderCosmic

friend class TBuilderCosmic

friend

Definition at line 376 of file TTrack.h.

TBuilderCurl

friend class TBuilderCurl

friend

Definition at line 377 of file TTrack.h.

TCosmicFitter

friend class TCosmicFitter

friend

Definition at line 372 of file TTrack.h.

TCurlFinder

friend class TCurlFinder

friend

Definition at line 379 of file TTrack.h.

THelixFitter

friend class THelixFitter

friend

Definition at line 371 of file TTrack.h.

TPMCurlFinder

friend class TPMCurlFinder

friend

Definition at line 378 of file TTrack.h.

TrkReco

friend class TrkReco

friend

Definition at line 370 of file TTrack.h.

TTrackManager

friend class TTrackManager

friend

Definition at line 369 of file TTrack.h.

---

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

• TTrack.h
• TTrack.cxx