TTrack Class Reference

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

#include <TTrack.h>

Inheritance diagram for TTrack:

List of all members.

Public Member Functions
void	append (const AList< TMLink > &)
	appends TMLinks.
void	append (TMLink &)
	appends a TMLink.
void	append (const AList< TMLink > &)
	appends TMLinks.
void	append (TMLink &)
	appends a TMLink.
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	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	appendByDistance (AList< TMLink > &list, double maxDistance)
	appends TMLinks by distance. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.
int	approach (TMLink &, bool sagCorrection) 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	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	approach2D (TMLink &) const
int	approach2D (TMLink &) const
void	assign (unsigned maskForWireHit)
	assigns wire hits to this track.
void	assign (unsigned maskForWireHit)
	assigns wire hits to this track.
const AList< TMLink > &	associateHits (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)
const AList< TMLink > &	associateHits (void) const
	sets/returns a list of associated TMLink which are used for table output.
TPoint2D	center (void) const
	returns position of helix center.
TPoint2D	center (void) const
	returns position of helix center.
double	charge (double)
	sets charge.
double	charge (void) const
	returns charge.
double	charge (double)
	sets charge.
double	charge (void) const
	returns charge.
double	chi2 (void) const
	returns chi2.
double	chi2 (void) const
	returns chi2.
double	confidenceLevel (void) const
	returns confidence level.
double	confidenceLevel (void) const
	returns confidence level.
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.
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.
TTrack *	daughter (TTrack *)
TTrack *	daughter (void) const
TTrack *	daughter (TTrack *)
TTrack *	daughter (void) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2, AList< HepPoint3D > &l3) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2, AList< HepPoint3D > &l3) const
void	deleteListForCurl (AList< HepPoint3D > &l1, AList< HepPoint3D > &l2) const
virtual double	distance (const TMLink &) const
	returns distance to a position of TMLink in TMLink space.
virtual double	distance (const TMLink &) const
	returns distance to a position of TMLink in TMLink space.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
void	falseFit ()
	false Fit
void	falseFit ()
	false Fit
const AList< TMLink > &	finalHits (const AList< TMLink > &hits)
const AList< TMLink > &	finalHits (void) const
	sets/returns a list of TMLink which are used for table output.
const AList< TMLink > &	finalHits (const AList< TMLink > &hits)
const AList< TMLink > &	finalHits (void) const
	sets/returns a list of TMLink which are used for table output.
unsigned	finder (unsigned finderMask)
unsigned	finder (void) const
	sets/returns finder.
unsigned	finder (unsigned finderMask)
unsigned	finder (void) const
	sets/returns finder.
virtual int	fit (void)
	fits itself by a default fitter. Error was happened if return value is not zero.
virtual int	fit (void)
	fits itself by a default fitter. Error was happened if return value is not zero.
int	fit2D (unsigned=0, double=0.1, double=0.015)
	fits itself with r-phi view. Error was happened if return value is not zero.
int	fit2D (unsigned=0, double=0.1, double=0.015)
	fits itself with r-phi view. Error was happened if return value is not zero.
bool	fitted (void) const
	returns true if fitted.
bool	fitted (void) const
	returns true if fitted.
bool	fittedWithCathode (void) const
	returns true if fitted with cathode hits(TEMPORARY).
bool	fittedWithCathode (void) const
	returns true if fitted with cathode hits(TEMPORARY).
const TMFitter *const	fitter (const TMFitter *)
	sets a default fitter.
const TMFitter *const	fitter (void) const
	returns a pointer to a default fitter.
const TMFitter *const	fitter (const TMFitter *)
	sets a default fitter.
const TMFitter *const	fitter (void) const
	returns a pointer to a default fitter.
unsigned	fitting (unsigned fitMask)
unsigned	fitting (void) const
	sets/returns fitting status.
unsigned	fitting (unsigned fitMask)
unsigned	fitting (void) const
	sets/returns fitting status.
int	HelCyl (double rhole, double rcyl, double zb, double zf, double epsl, double &phi, HepPoint3D &xp) const
	calculates an intersection of this track and a cylinder.
int	HelCyl (double rhole, double rcyl, double zb, double zf, double epsl, double &phi, HepPoint3D &xp) const
	calculates an intersection of this track and a cylinder.
const Helix &	helix (void) const
	returns helix parameter.
const Helix &	helix (void) const
	returns helix parameter.
const TTrackHEP *const	hep (void) const
	returns TTrackHEP.
const TTrackHEP *const	hep (void) const
	returns TTrackHEP.
double	impact (void) const
	returns signed impact parameter to the origin.
double	impact (void) const
	returns signed impact parameter to the origin.
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.
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.
const TTrackMC *const	mc (void) const
	returns a pointer to TTrackMC.
const TTrackMC *const	mc (void) const
	returns a pointer to TTrackMC.
TTrack *	mother (TTrack *)
TTrack *	mother (void) const
	sets/returns mother/daughter.
TTrack *	mother (TTrack *)
TTrack *	mother (void) const
	sets/returns mother/daughter.
void	movePivot (void)
	moves pivot to the inner most hit.
void	movePivot (void)
	moves pivot to the inner most hit.
const std::string &	name (const std::string &newName)
const std::string &	name (void) const
	returns/sets name.
const std::string &	name (const std::string &newName)
const std::string &	name (void) const
	returns/sets name.
unsigned	nCores (unsigned mask=0) const
	returns # 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.
unsigned	ndf (void) const
	returns NDF.
unsigned	ndf (void) const
	returns NDF.
unsigned	nHeps (void) const
	returns # of contributed TTrackHEP tracks.
unsigned	nHeps (void) const
	returns # of contributed TTrackHEP tracks.
unsigned	nLinks (unsigned mask=0) const
	returns # of masked TMLinks assigned to this track object.
unsigned	nLinks (unsigned mask=0) const
	returns # of masked TMLinks assigned to this track object.
unsigned	objectType (void) const
	returns type.
unsigned	objectType (void) const
	returns type.
TMLink *	operator[] (unsigned i) const
TMLink *	operator[] (unsigned i) const
Hep3Vector	p (void) const
	returns momentum.
Hep3Vector	p (void) const
	returns momentum.
double	pt (void) const
	returns Pt.
double	pt (void) const
	returns Pt.
double	ptot (void) const
	returns magnitude of momentum.
double	ptot (void) const
	returns magnitude of momentum.
double	pz (void) const
	returns Pz.
double	pz (void) const
	returns Pz.
unsigned	quality (unsigned qualityMask)
unsigned	quality (void) const
	sets/returns quality.
unsigned	quality (unsigned qualityMask)
unsigned	quality (void) const
	sets/returns quality.
double	radius (void) const
	returns signed radius.
double	radius (void) const
	returns signed radius.
virtual void	refine (double maxSigma)
	removes bad points by pull. The bad points are masked not to be used in fit.
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 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'.
void	refine2D (AList< TMLink > &list, float maxSigma)
	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))
void	refine2D (AList< TMLink > &list, float maxSigma)
	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))
void	remove (const AList< TMLink > &)
	removes TMLinks.
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
void	remove (TMLink &a)
	removes a TMLink.
virtual void	removeLinks (void)
virtual void	removeLinks (void)
const AList< TSegment > &	segments (void) const
AList< TSegment > &	segments (void)
	returns AList<TSegment>.
const AList< TSegment > &	segments (void) const
AList< TSegment > &	segments (void)
	returns AList<TSegment>.
unsigned	state (void) const
	returns/sets internal state.(for bank output)
unsigned	state (void) const
	returns/sets internal state.(for bank output)
int	stereoHitForCurl (AList< TMLink > &) const
	calculates arc length and z for a stereo hit. uses these functions for curl tracks(included svd version).
int	stereoHitForCurl (TMLink &link, TMLink &link1, TMLink &link2) const
int	stereoHitForCurl (TMLink &link, TMLink &link1) const
int	stereoHitForCurl (TMLink &link, AList< HepPoint3D > &arcZList) const
int	stereoHitForCurl (AList< TMLink > &) const
	calculates arc length and z for a stereo hit. uses these functions for curl tracks(included svd version).
int	stereoHitForCurl (TMLink &link, TMLink &link1, TMLink &link2) const
int	stereoHitForCurl (TMLink &link, TMLink &link1) const
int	stereoHitForCurl (TMLink &link, AList< HepPoint3D > &arcZList) const
int	szPosition (const HepPoint3D &p, HepPoint3D &szPosition) const
	calculates arc length for a point.
int	szPosition (const TSegment &segment, TMLink &link) const
	calculates arc length and z for a segment. Results are stored in TMLink.
int	szPosition (TMLink &link) const
	calculates arc length and z for a stereo hit.
int	szPosition (const HepPoint3D &p, HepPoint3D &szPosition) const
	calculates arc length for a point.
int	szPosition (const TSegment &segment, TMLink &link) const
	calculates arc length and z for a segment. Results are stored in TMLink.
int	szPosition (TMLink &link) const
	calculates arc length and z for a stereo hit.
unsigned	testByApproach (const AList< TMLink > &list, double sigma) 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
unsigned	testByApproach (const TMLink &list, double sigma) const
	returns # of good hits to be appended.
	TTrack (const Helix &)
	Constructor.
	TTrack (const T3DLine &)
	Constructor.
	TTrack (const TTrack &)
	Copy constructor.
	TTrack (const TCircle &)
	Constructor.
	TTrack ()
	Default constructor.
	TTrack (const Helix &)
	Constructor.
	TTrack (const T3DLine &)
	Constructor.
	TTrack (const TTrack &)
	Copy constructor.
	TTrack (const TCircle &)
	Constructor.
	TTrack ()
	Default constructor.
unsigned	type (void) const
	returns type. Definition is depending on an object type.
unsigned	type (void) const
	returns type. Definition is depending on an object type.
void	update (void) const
	update cache.
void	update (void) const
	update cache.
virtual	~TTrack ()
	Destructor.
virtual	~TTrack ()
	Destructor.
Protected Attributes
bool	_fitted
bool	_fittedWithCathode
AList< TMLink >	_links
AList< TMLink >	_links
TTrackMC *	_mc
TTrackMC *	_mc
Private Member Functions
unsigned	defineType (void) const
unsigned	defineType (void) const
int	dxda (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda, Vector &dzda) const
	calculates dXda. 'link' and 'dPhi' are inputs. Others are outputs.
int	dxda (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda, Vector &dzda) const
	calculates dXda. 'link' and 'dPhi' are inputs. Others are outputs.
int	dxda2D (double dPhi, Vector &dxda, Vector &dyda) const
int	dxda2D (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda) const
int	dxda2D (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda, Vector &dzda) const
int	dxda2D (double dPhi, Vector &dxda, Vector &dyda) const
int	dxda2D (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda) const
int	dxda2D (const TMLink &link, double dPhi, Vector &dxda, Vector &dyda, Vector &dzda) const
Private Attributes
AList< TMLink >	_associateHits
AList< TMLink >	_associateHits
double	_charge
double	_chi2
TTrack *	_daughter
TTrack *	_daughter
AList< TMLink >	_finalHits
AList< TMLink >	_finalHits
Helix *const	_helix
Helix *const	_helix
TTrack *	_mother
TTrack *	_mother
std::string	_name
unsigned	_ndf
AList< TSegment >	_segments
AList< TSegment >	_segments
unsigned	_state
unsigned	_type
IMagneticFieldSvc *	m_pmgnIMF
IMagneticFieldSvc *	m_pmgnIMF
Static Private Attributes
const THelixFitter	_fitter = THelixFitter("TTrack Default Helix Fitter")
Friends
class	Refit
class	TBuilder
class	TBuilder0
class	TBuilderCosmic
class	TBuilderCurl
class	TCosmicFitter
class	TCurlFinder
class	THelixFitter
class	TPMCurlFinder
class	TrkReco
class	TTrackManager

---

Detailed Description

A class to represent a track in tracking.

---

Constructor & Destructor Documentation

TTrack::TTrack (   )

Default constructor. : 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<< "Unable to open Magnetic field service"<<std::endl; } }

TTrack::TTrack (  const TCircle &   )

Constructor. : 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<< "Unable to open Magnetic field service"<<std::endl; } //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(); a[2] = 333.564095/(-1000*(m_pmgnIMF->getReferField())) / c.radius(); 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<< "Unable to open Magnetic field service"<<std::endl; } */ }

TTrack::TTrack (  const TTrack &   )

Copy constructor. : 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<< "Unable to open Magnetic field service"<<std::endl; } }

TTrack::TTrack (  const T3DLine &   )

Constructor.

TTrack::TTrack (  const Helix &   )

Constructor. : 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<< "Unable to open Magnetic field service"<<std::endl; } }

TTrack::~TTrack (   )  [virtual]

Destructor. { delete _helix; }

TTrack::TTrack (   )

Default constructor.

TTrack::TTrack (  const TCircle &   )

Constructor.

TTrack::TTrack (  const TTrack &   )

Copy constructor.

TTrack::TTrack (  const T3DLine &   )

Constructor.

TTrack::TTrack (  const Helix &   )

Constructor.

virtual TTrack::~TTrack (   )  [virtual]

Destructor.

---

Member Function Documentation

void TTrackBase::append (  const AList< TMLink > &   )  [inherited]

appends TMLinks.

void TTrackBase::append (  TMLink &   )  [inherited]

appends a TMLink.

void TTrackBase::append (  const AList< TMLink > &   )  [inherited]

appends TMLinks. { AList<TMLink> tmp; for (unsigned i = 0; i < a.length(); i++) { if ((_links.hasMember(a[i])) || (a[i]->hit()->state() & WireHitUsed)) { #ifdef TRKRECO_DEBUG_DETAIL std::cout << " TTrackBase::append(list) !!! "; std::cout << a[i]->wire()->name(); std::cout << " Hey!, this is already used! Don't mess with me."; std::cout << std::endl; #endif continue; } else { tmp.append(a[i]); } } _links.append(tmp); _updated = false; _fitted = false; _fittedWithCathode = false; // added by matsu ( 1999/05/24 ) }

void TTrackBase::append (  TMLink &   )  [inherited]

appends a TMLink. { if ((a.hit()->state() & WireHitUsed)) { #ifdef TRKRECO_DEBUG_DETAIL std::cout << "TTrackBase::append !!! " << a.wire()->name() << " this is already used by another track!" << std::endl; #endif return; } if (_links.hasMember(a)) { #ifdef TRKRECO_DEBUG_DETAIL std::cout << "TTrackBase::append !!! " << a.wire()->name() << " this is already included in this track!" << std::endl; #endif return; } _links.append(a); _updated = false; _fitted = false; _fittedWithCathode = false; // added by matsu ( 1999/05/24 ) }

void TTrackBase::appendByApproach (  AList< TMLink > &  list, double  maxSigma )  [inherited]

appends TMLinks by approach. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.

void TTrackBase::appendByApproach (  AList< TMLink > &  list, double  maxSigma )  [inherited]

appends TMLinks by approach. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned. { #ifdef TRKRECO_DEBUG_DETAIL std::cout << " TTrackBase::appendByApproach ... sigma=" << maxSigma << std::endl; #endif AList<TMLink> unused; unsigned n = list.length(); for (unsigned i = 0; i < n; i++) { TMLink & l = * list[i]; if ((_links.hasMember(l)) || (l.hit()->state() & WireHitUsed)) continue; //...Calculate closest approach... int err = approach(l); if (err < 0) { unused.append(l); continue; } //...Calculate sigma... float distance = (l.positionOnWire() - l.positionOnTrack()).mag(); float diff = fabs(distance - l.drift()); float sigma = diff / l.dDrift(); //...For debug... #ifdef TRKRECO_DEBUG_DETAIL std::cout << " sigma=" << sigma; std::cout << ",dist=" << distance; std::cout << ",diff=" << diff; std::cout << ",err=" << l.hit()->dDrift() << ","; if (sigma < maxSigma) std::cout << "ok,"; else std::cout << "X,"; l.dump("mc"); #endif //...Make sigma cut... if (sigma > maxSigma) { unused.append(l); continue; } //...OK... _links.append(l); _updated = false; _fitted = false; } list.remove(unused); }

void TTrackBase::appendByDistance (  AList< TMLink > &  list, double  maxDistance )  [inherited]

appends TMLinks by distance. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.

void TTrackBase::appendByDistance (  AList< TMLink > &  list, double  maxDistance )  [inherited]

appends TMLinks by distance. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned. { std::cout << "TTrackBase::appendByDistance !!! not implemented" << std::endl; list.removeAll(); }

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

int TTrack::approach (  TMLink &   )  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.

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

{ //...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(); //...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; double rho = 333.564095/(-1000*(m_pmgnIMF->getReferField())) / kappa; 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; }

int TTrack::approach (  TMLink &   )  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. { return approach(l, false); }

int TTrack::approach2D (  TMLink &   )  const

int TTrack::approach2D (  TMLink &   )  const

{ 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; }

void TTrack::assign (  unsigned  maskForWireHit  )

assigns wire hits to this track.

void TTrack::assign (  unsigned  maskForWireHit  )

assigns wire hits to this track. { 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); } }

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

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

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

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

{ _associateHits = list; return _associateHits; }

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

sets/returns a list of associated TMLink which are used for table output. { return _associateHits; }

TPoint2D TTrack::center (  void   )  const

returns position of helix center.

TPoint2D TTrack::center (  void   )  const [inline]

returns position of helix center. { return TPoint2D(_helix->center()); }

double TTrack::charge (  double   )

sets charge.

double TTrack::charge (  void   )  const

returns charge.

double TTrack::charge (  double   )  [inline]

sets charge. { return _charge = a; }

double TTrack::charge (  void   )  const [inline]

returns charge. { return _charge; }

double TTrack::chi2 (  void   )  const

returns chi2.

double TTrack::chi2 (  void   )  const [inline]

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

double TTrack::confidenceLevel (  void   )  const

returns confidence level.

double TTrack::confidenceLevel (  void   )  const [inline]

returns confidence level. { return chisq2confLevel((int) _ndf, _chi2); }

const AList<TMLink>& TTrackBase::cores (  unsigned  mask = 0  )  const [inherited]

returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0.

const AList< TMLink > & TTrackBase::cores (  unsigned  mask = 0  )  const [inherited]

returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0. { if (mask) std::cout << "TTrackBase::cores !!! mask is not supported" << std::endl; if (! _updated) update(); return _cores; }

TTrack* TTrack::daughter (  TTrack *   )

TTrack* TTrack::daughter (  void   )  const

TTrack * TTrack::daughter (  TTrack *   )  [inline]

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

TTrack * TTrack::daughter (  void   )  const [inline]

{ return _daughter; }

unsigned TTrack::defineType (  void   )  const [private]

unsigned TTrack::defineType (  void   )  const [private]

{ bool highPt = true; if (pt() < 0.2) highPt = false; //Bes bool fromIP = true; if (fabs(impact()) > 8.3) fromIP = false; if (fromIP && highPt) return _type = TrackTypeNormal; else if (fromIP && (! highPt)) return _type = TrackTypeCurl; if ((fabs(radius()) * 2. + fabs(impact())) < 81.) //Bes return _type = TrackTypeCircle; return _type = TrackTypeCosmic; }

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

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

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

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

virtual double TTrackBase::distance (  const TMLink &   )  const [virtual, inherited]

returns distance to a position of TMLink in TMLink space. Reimplemented in TLine0, TMLine, TLine0, and TMLine.

double TTrackBase::distance (  const TMLink &   )  const [virtual, inherited]

returns distance to a position of TMLink in TMLink space. Reimplemented in TLine0, TMLine, TLine0, and TMLine. { std::cout << "TTrackBase::distance !!! not implemented" << std::endl; return 0.; }

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

dumps debug information. Reimplemented from TTrackBase.

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

dumps debug information. Reimplemented from TTrackBase. { 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); }

int TTrack::dxda (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda, Vector &  dzda )  const [private]

calculates dXda. 'link' and 'dPhi' are inputs. Others are outputs.

int TTrack::dxda (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda, Vector &  dzda )  const [private]

calculates dXda. 'link' and 'dPhi' are inputs. Others are outputs. { //...Setup... const TMDCWire & w = * link.wire(); Vector a = _helix->a(); double dRho = a[0]; double phi0 = a[1]; double kappa = a[2]; // double rho = Helix::ConstantAlpha / kappa; // double rho = 333.564095 / kappa; double rho = 333.564095/(-1000*(m_pmgnIMF->getReferField())) / kappa; double tanLambda = a[4]; double sinPhi0 = sin(phi0); double cosPhi0 = cos(phi0); double sinPhi0dPhi = sin(phi0 + dPhi); double cosPhi0dPhi = cos(phi0 + dPhi); Vector dphida(5); //...Axial case... if (w.axial()) { HepPoint3D d = _helix->center() - w.xyPosition(); double dmag2 = d.mag2(); dphida[0] = (sinPhi0 * d.x() - cosPhi0 * d.y()) / dmag2; dphida[1] = (dRho + rho) * (cosPhi0 * d.x() + sinPhi0 * d.y()) / dmag2 - 1.; dphida[2] = (- rho / kappa) * (sinPhi0 * d.x() - cosPhi0 * d.y()) / dmag2; dphida[3] = 0.; dphida[4] = 0.; } //...Stereo case... else { HepPoint3D onTrack = _helix->x(dPhi); HepVector3D v = w.direction(); Vector c(3); c = HepPoint3D(w.backwardPosition() - (v * w.backwardPosition()) * v); Vector x(3); x = onTrack; Vector dxdphi(3); dxdphi[0] = rho * sinPhi0dPhi; dxdphi[1] = - rho * cosPhi0dPhi; dxdphi[2] = - rho * tanLambda; Vector d2xdphi2(3); d2xdphi2[0] = rho * cosPhi0dPhi; d2xdphi2[1] = rho * sinPhi0dPhi; d2xdphi2[2] = 0.; double dxdphi_dot_v = dxdphi[0]*v.x() + dxdphi[1]*v.y() + dxdphi[2]*v.z(); double x_dot_v = x[0]*v.x() + x[1]*v.y() + x[2]*v.z(); double dfdphi = - (dxdphi[0] - dxdphi_dot_v*v.x()) * dxdphi[0] - (dxdphi[1] - dxdphi_dot_v*v.y()) * dxdphi[1] - (dxdphi[2] - dxdphi_dot_v*v.z()) * dxdphi[2] - (x[0] - c[0] - x_dot_v*v.x()) * d2xdphi2[0] - (x[1] - c[1] - x_dot_v*v.y()) * d2xdphi2[1]; /* - (x[2] - c[2] - x_dot_v*v.z()) * d2xdphi2[2]; = 0. */ //dxda_phi, dyda_phi, dzda_phi : phi is fixed Vector dxda_phi(5); dxda_phi[0] = cosPhi0; dxda_phi[1] = -(dRho + rho)*sinPhi0 + rho*sinPhi0dPhi; dxda_phi[2] = -(rho/kappa)*( cosPhi0 - cosPhi0dPhi ); dxda_phi[3] = 0.; dxda_phi[4] = 0.; Vector dyda_phi(5); dyda_phi[0] = sinPhi0; dyda_phi[1] = (dRho + rho)*cosPhi0 - rho*cosPhi0dPhi; dyda_phi[2] = -(rho/kappa)*( sinPhi0 - sinPhi0dPhi ); dyda_phi[3] = 0.; dyda_phi[4] = 0.; Vector dzda_phi(5); dzda_phi[0] = 0.; dzda_phi[1] = 0.; dzda_phi[2] = (rho/kappa)*tanLambda*dPhi; dzda_phi[3] = 1.; dzda_phi[4] = -rho*dPhi; Vector d2xdphida(5); d2xdphida[0] = 0.; d2xdphida[1] = rho*cosPhi0dPhi; d2xdphida[2] = -(rho/kappa)*sinPhi0dPhi; d2xdphida[3] = 0.; d2xdphida[4] = 0.; Vector d2ydphida(5); d2ydphida[0] = 0.; d2ydphida[1] = rho*sinPhi0dPhi; d2ydphida[2] = (rho/kappa)*cosPhi0dPhi; d2ydphida[3] = 0.; d2ydphida[4] = 0.; Vector d2zdphida(5); d2zdphida[0] = 0.; d2zdphida[1] = 0.; d2zdphida[2] = (rho/kappa)*tanLambda; d2zdphida[3] = 0.; d2zdphida[4] = -rho; Vector dfda(5); for( int i = 0; i < 5; i++ ){ double d_dot_v = v.x()*dxda_phi[i] + v.y()*dyda_phi[i] + v.z()*dzda_phi[i]; dfda[i] = - (dxda_phi[i] - d_dot_v*v.x()) * dxdphi[0] - (dyda_phi[i] - d_dot_v*v.y()) * dxdphi[1] - (dzda_phi[i] - d_dot_v*v.z()) * dxdphi[2] - (x[0] - c[0] - x_dot_v*v.x()) * d2xdphida[i] - (x[1] - c[1] - x_dot_v*v.y()) * d2ydphida[i] - (x[2] - c[2] - x_dot_v*v.z()) * d2zdphida[i]; dphida[i] = - dfda[i] /dfdphi; } } dxda[0] = cosPhi0 + rho * sinPhi0dPhi * dphida[0]; dxda[1] = - (dRho + rho) * sinPhi0 + rho * sinPhi0dPhi * (1. + dphida[1]); dxda[2] = - rho / kappa * (cosPhi0 - cosPhi0dPhi) + rho * sinPhi0dPhi * dphida[2]; dxda[3] = rho * sinPhi0dPhi * dphida[3]; dxda[4] = rho * sinPhi0dPhi * dphida[4]; dyda[0] = sinPhi0 - rho * cosPhi0dPhi * dphida[0]; dyda[1] = (dRho + rho) * cosPhi0 - rho * cosPhi0dPhi * (1. + dphida[1]); dyda[2] = - rho / kappa * (sinPhi0 - sinPhi0dPhi) - rho * cosPhi0dPhi * dphida[2]; dyda[3] = - rho * cosPhi0dPhi * dphida[3]; dyda[4] = - rho * cosPhi0dPhi * dphida[4]; dzda[0] = - rho * tanLambda * dphida[0]; dzda[1] = - rho * tanLambda * dphida[1]; dzda[2] = rho / kappa * tanLambda * dPhi - rho * tanLambda * dphida[2]; dzda[3] = 1. - rho * tanLambda * dphida[3]; dzda[4] = - rho * dPhi - rho * tanLambda * dphida[4]; return 0; }

int TTrack::dxda2D (  double  dPhi, Vector &  dxda, Vector &  dyda )  const [private]

int TTrack::dxda2D (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda )  const [private]

int TTrack::dxda2D (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda, Vector &  dzda )  const [private]

int TTrack::dxda2D (  double  dPhi, Vector &  dxda, Vector &  dyda )  const [private]

{ //...Setup... double kappa = (_helix->a())[2]; if(kappa == 0.){ std::cout << "Error(?) : kappa == 0 in dxda2D of TrkReco." << std::endl; return 1; } double dRho = (_helix->a())[0]; double phi0 = (_helix->a())[1]; // double rho = Helix::ConstantAlpha / kappa; // double rho = 333.564095 / kappa; double rho = 333.564095/(-1000*(m_pmgnIMF->getReferField())) / kappa; double sinPhi0 = sin(phi0); double cosPhi0 = cos(phi0); double sinPhi0dPhi = sin(phi0 + dPhi); double cosPhi0dPhi = cos(phi0 + dPhi); Vector dphida(3); HepPoint3D d = _helix->center(); // d = center - (0,0,0) d.setZ(0.); double dmag2 = d.mag2(); if(dmag2 == 0.){ std::cout << "Error(?) : Distance(center-xyPosition) == 0 in dxda2D of TrkReco." << std::endl; return 1; } dphida[0] = (sinPhi0*d.x()-cosPhi0*d.y())/dmag2; dphida[1] = (dRho+rho)*(cosPhi0*d.x()+sinPhi0 * d.y())/dmag2 - 1.; dphida[2] = (-rho/kappa)*(sinPhi0*d.x()-cosPhi0*d.y())/dmag2; dxda[0] = cosPhi0+rho*sinPhi0dPhi*dphida[0]; dxda[1] = -(dRho+rho)*sinPhi0+rho*sinPhi0dPhi*(1.+dphida[1]); dxda[2] = -rho/kappa*(cosPhi0-cosPhi0dPhi)+rho*sinPhi0dPhi*dphida[2]; dyda[0] = sinPhi0-rho*cosPhi0dPhi*dphida[0]; dyda[1] = (dRho+rho)*cosPhi0-rho*cosPhi0dPhi*(1.+dphida[1]); dyda[2] = -rho/kappa*(sinPhi0-sinPhi0dPhi)-rho*cosPhi0dPhi*dphida[2]; return 0; }

int TTrack::dxda2D (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda )  const [private]

{ //...Setup... double kappa = (_helix->a())[2]; if(kappa == 0.){ std::cout << "Error(?) : kappa == 0 in dxda2D of TrkReco." << std::endl; return 1; } const TMDCWire &w = *link.wire(); double dRho = (_helix->a())[0]; double phi0 = (_helix->a())[1]; // double rho = Helix::ConstantAlpha / kappa; // double rho = 333.564095 / kappa; double rho = 333.564095/(-1000*(m_pmgnIMF->getReferField())) / kappa; double sinPhi0 = sin(phi0); double cosPhi0 = cos(phi0); double sinPhi0dPhi = sin(phi0 + dPhi); double cosPhi0dPhi = cos(phi0 + dPhi); Vector dphida(3); HepPoint3D d = _helix->center() - w.xyPosition(); d.setZ(0.); double dmag2 = d.mag2(); if(dmag2 == 0.){ std::cout << "Error(?) : Distance(center-xyPosition) == 0 in dxda2D of TrkReco." << std::endl; return 1; } dphida[0] = (sinPhi0*d.x()-cosPhi0*d.y())/dmag2; dphida[1] = (dRho+rho)*(cosPhi0*d.x()+sinPhi0 * d.y())/dmag2 - 1.; dphida[2] = (-rho/kappa)*(sinPhi0*d.x()-cosPhi0*d.y())/dmag2; dxda[0] = cosPhi0+rho*sinPhi0dPhi*dphida[0]; dxda[1] = -(dRho+rho)*sinPhi0+rho*sinPhi0dPhi*(1.+dphida[1]); dxda[2] = -rho/kappa*(cosPhi0-cosPhi0dPhi)+rho*sinPhi0dPhi*dphida[2]; dyda[0] = sinPhi0-rho*cosPhi0dPhi*dphida[0]; dyda[1] = (dRho+rho)*cosPhi0-rho*cosPhi0dPhi*(1.+dphida[1]); dyda[2] = -rho/kappa*(sinPhi0-sinPhi0dPhi)-rho*cosPhi0dPhi*dphida[2]; return 0; }

int TTrack::dxda2D (  const TMLink &  link, double  dPhi, Vector &  dxda, Vector &  dyda, Vector &  dzda )  const [private]

void TTrackBase::falseFit (   )  [inherited]

false Fit

void TTrackBase::falseFit (   )  [inline, inherited]

false Fit { _fitted = false; _fittedWithCathode = false; }

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

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

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

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

{ _finalHits = list; return _finalHits; }

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

sets/returns a list of TMLink which are used for table output. { return _finalHits; }

unsigned TTrack::finder (  unsigned  finderMask  )

unsigned TTrack::finder (  void   )  const

sets/returns finder.

unsigned TTrack::finder (  unsigned  finderMask  )  [inline]

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

unsigned TTrack::finder (  void   )  const [inline]

sets/returns finder. { return _state & TrackFinderMask; }

virtual int TTrackBase::fit (  void   )  [virtual, inherited]

fits itself by a default fitter. Error was happened if return value is not zero.

int TTrackBase::fit (  void   )  [virtual, inherited]

fits itself by a default fitter. Error was happened if return value is not zero. { return _fitter->fit(* this); }

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

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

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

fits itself with r-phi view. Error was happened if return value is not zero. { #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; }

bool TTrackBase::fitted (  void   )  const [inherited]

returns true if fitted.

bool TTrackBase::fitted (  void   )  const [inline, inherited]

returns true if fitted. { return _fitted; }

bool TTrackBase::fittedWithCathode (  void   )  const [inherited]

returns true if fitted with cathode hits(TEMPORARY).

bool TTrackBase::fittedWithCathode (  void   )  const [inline, inherited]

returns true if fitted with cathode hits(TEMPORARY). { return _fittedWithCathode; }

const TMFitter* const TTrackBase::fitter (  const TMFitter *   )  [inherited]

sets a default fitter.

const TMFitter* const TTrackBase::fitter (  void   )  const [inherited]

returns a pointer to a default fitter.

const TMFitter *const TTrackBase::fitter (  const TMFitter *   )  [inline, inherited]

sets a default fitter. { _fitted = false; return _fitter = a; }

const TMFitter *const TTrackBase::fitter (  void   )  const [inline, inherited]

returns a pointer to a default fitter. { return _fitter; }

unsigned TTrack::fitting (  unsigned  fitMask  )

unsigned TTrack::fitting (  void   )  const

sets/returns fitting status.

unsigned TTrack::fitting (  unsigned  fitMask  )  [inline]

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

unsigned TTrack::fitting (  void   )  const [inline]

sets/returns fitting status. { return (_state >> TrackFitShift) & TrackFitMask; }

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

calculates an intersection of this track and a cylinder.

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

calculates an intersection of this track and a cylinder. { 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 ; }

const Helix& TTrack::helix (  void   )  const

returns helix parameter.

const Helix & TTrack::helix (  void   )  const [inline]

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

const TTrackHEP* const TTrackBase::hep (  void   )  const [inherited]

returns TTrackHEP.

const TTrackHEP *const TTrackBase::hep (  void   )  const [inherited]

returns TTrackHEP. { unsigned n = _links.length(); CAList<TTrackHEP> hepList; CList<unsigned> hepCounter; for (unsigned i = 0; i < n; i++) { const TTrackHEP * hep = _links[i]->hit()->mc()->hep(); unsigned nH = hepList.length(); bool found = false; for (unsigned j = 0; j < nH; j++) { if (hepList[j] == hep) { found = true; ++(* hepCounter[j]); } } if (! found) { hepList.append(hep); unsigned c = 0; hepCounter.append(c); } } _nHeps = hepList.length(); _hep = 0; unsigned max = 0; for (unsigned i = 0; i < _nHeps; i++) { if ((* hepCounter[i]) > max) { max = (* hepCounter[i]); _hep = hepList[i]; } } return _hep; }

double TTrack::impact (  void   )  const

returns signed impact parameter to the origin.

double TTrack::impact (  void   )  const [inline]

returns signed impact parameter to the origin. { return fabs(_helix->radius()) - _helix->center().mag(); }

const AList<TMLink>& TTrackBase::links (  unsigned  mask = 0  )  const [inherited]

returns a list of masked TMLinks assigned to this track. 'mask' will be applied if mask is not 0.

const AList< TMLink > & TTrackBase::links (  unsigned  mask = 0  )  const [inherited]

returns a list of masked TMLinks assigned to this track. 'mask' will be applied if mask is not 0. { if (mask == 0) return _links; std::cout << "TTrackBase::links !!! mask is not supportted yet" << std::endl; return _links; }

const TTrackMC* const TTrackBase::mc (  void   )  const [inherited]

returns a pointer to TTrackMC.

const TTrackMC *const TTrackBase::mc (  void   )  const [inline, inherited]

returns a pointer to TTrackMC. { return _mc; }

TTrack* TTrack::mother (  TTrack *   )

TTrack* TTrack::mother (  void   )  const

sets/returns mother/daughter.

TTrack * TTrack::mother (  TTrack *   )  [inline]

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

TTrack * TTrack::mother (  void   )  const [inline]

sets/returns mother/daughter. { return _mother; }

void TTrack::movePivot (  void   )

moves pivot to the inner most hit.

void TTrack::movePivot (  void   )

moves pivot to the inner most hit. { 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); }

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

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

returns/sets name.

const std::string & TTrack::name (  const std::string &  newName  )  [inline]

{ return _name = a; }

const std::string & TTrack::name (  void   )  const [inline]

returns/sets name. { return _name; }

unsigned TTrackBase::nCores (  unsigned  mask = 0  )  const [inherited]

returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0.

unsigned TTrackBase::nCores (  unsigned  mask = 0  )  const [inherited]

returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0. { if (mask) std::cout << "TTrackBase::nCores !!! mask is not supported" << std::endl; if (! _updated) update(); return _cores.length(); }

unsigned TTrack::ndf (  void   )  const

returns NDF.

unsigned TTrack::ndf (  void   )  const [inline]

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

unsigned TTrackBase::nHeps (  void   )  const [inherited]

returns # of contributed TTrackHEP tracks.

unsigned TTrackBase::nHeps (  void   )  const [inherited]

returns # of contributed TTrackHEP tracks. { hep(); return _nHeps; }

unsigned TTrackBase::nLinks (  unsigned  mask = 0  )  const [inherited]

returns # of masked TMLinks assigned to this track object.

unsigned TTrackBase::nLinks (  unsigned  mask = 0  )  const [inherited]

returns # of masked TMLinks assigned to this track object. { unsigned n = _links.length(); if (mask == 0) return n; unsigned nn = 0; for (unsigned i = 0; i < n; i++) { const TMDCWireHit & h = * _links[i]->hit(); if (h.state() & mask) ++nn; } return nn; }

unsigned TTrack::objectType (  void   )  const [virtual]

returns type. Reimplemented from TTrackBase.

unsigned TTrack::objectType (  void   )  const [inline, virtual]

returns type. Reimplemented from TTrackBase. { return Track; }

TMLink* TTrackBase::operator[] (  unsigned  i  )  const [inherited]

TMLink * TTrackBase::operator[] (  unsigned  i  )  const [inline, inherited]

{ return _links[i]; }

Hep3Vector TTrack::p (  void   )  const

returns momentum.

Hep3Vector TTrack::p (  void   )  const [inline]

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

double TTrack::pt (  void   )  const

returns Pt.

double TTrack::pt (  void   )  const [inline]

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

double TTrack::ptot (  void   )  const

returns magnitude of momentum.

double TTrack::ptot (  void   )  const [inline]

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

double TTrack::pz (  void   )  const

returns Pz.

double TTrack::pz (  void   )  const [inline]

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

unsigned TTrack::quality (  unsigned  qualityMask  )

unsigned TTrack::quality (  void   )  const

sets/returns quality.

unsigned TTrack::quality (  unsigned  qualityMask  )  [inline]

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

unsigned TTrack::quality (  void   )  const [inline]

sets/returns quality. { return (_state >> TrackQualityShift) & TrackQualityMask; }

double TTrack::radius (  void   )  const

returns signed radius.

double TTrack::radius (  void   )  const [inline]

returns signed radius. { return _helix->radius(); }

virtual void TTrackBase::refine (  double  maxSigma  )  [virtual, inherited]

removes bad points by pull. The bad points are masked not to be used in fit.

virtual void TTrackBase::refine (  AList< TMLink > &  list, double  maxSigma )  [virtual, inherited]

removes bad points by pull. The bad points are removed from the track, and are returned in 'list'.

void TTrackBase::refine (  double  maxSigma  )  [virtual, inherited]

removes bad points by pull. The bad points are masked not to be used in fit. { AList<TMLink> bad = refineMain(sigma); // for (unsigned i = 0; i < bad.length(); i++) { // const TMDCWireHit * hit = bad[i]->hit(); // hit->state(hit->state() | WireHitInvalidForFit); // } #ifdef TRKRECO_DEBUG_DETAIL std::cout << " refine ... sigma=" << sigma << std::endl; Dump(bad, "detail sort", " "); #endif if (bad.length()) { _fitted = false; _updated = false; } }

void TTrackBase::refine (  AList< TMLink > &  list, double  maxSigma )  [virtual, inherited]

removes bad points by pull. The bad points are removed from the track, and are returned in 'list'. { AList<TMLink> bad = refineMain(sigma); #ifdef TRKRECO_DEBUG std::cout << " refine ... sigma=" << sigma << ", # of rejected hits="; std::cout << bad.length() << std::endl; #endif #ifdef TRKRECO_DEBUG Dump(bad, "sort pull mc", " "); #endif if (bad.length()) { _links.remove(bad); list.append(bad); _fitted = false; _updated = false; } }

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

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))

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

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)) { 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); }

void TTrackBase::remove (  const AList< TMLink > &   )  [inherited]

removes TMLinks.

void TTrackBase::remove (  TMLink &  a  )  [inherited]

removes a TMLink.

void TTrackBase::remove (  const AList< TMLink > &   )  [inline, inherited]

removes TMLinks. { _links.remove(a); _updated = false; _fitted = false; _fittedWithCathode = false; // mod. by matsu ( 1999/05/24 ) }

void TTrackBase::remove (  TMLink &  a  )  [inline, inherited]

removes a TMLink. { _links.remove(a); _updated = false; _fitted = false; _fittedWithCathode = false; // mod. by matsu ( 1999/05/24 ) }

virtual void TTrackBase::removeLinks (  void   )  [virtual, inherited]

void TTrackBase::removeLinks (  void   )  [virtual, inherited]

{ _links.removeAll(); }

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

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

returns AList<TSegment>.

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

{ return _segments; }

AList< TSegment > & TTrack::segments (  void   )  [inline]

returns AList<TSegment>. { return _segments; }

unsigned TTrack::state (  void   )  const

returns/sets internal state.(for bank output)

unsigned TTrack::state (  void   )  const [inline]

returns/sets internal state.(for bank output) { return _state; }

int TTrack::stereoHitForCurl (  AList< TMLink > &   )  const

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

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

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

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

int TTrack::stereoHitForCurl (  AList< TMLink > &   )  const

calculates arc length and z for a stereo hit. uses these functions for curl tracks(included svd version). { 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; }

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

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

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

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

calculates arc length for a point.

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

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

int TTrack::szPosition (  TMLink &  link  )  const

calculates arc length and z for a stereo hit.

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

calculates arc length for a point. { //...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; }

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

calculates arc length and z for a segment. Results are stored in TMLink. { //...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(); // std::cout<<"minDist szPosition TTrack:"<<minDist<<endl; for (unsigned i = 1; i < n; i++) { if (links[i]->hit()->drift() < minDist) { minDist = links[i]->drift(); // std::cout<<"minDist szPosition TTrack:"<<minDist<<endl; minL = links[i]; } } //...sz calculation... a.position(minL->position()); a.leftRight(2); a.hit(minL->hit()); int err = szPosition(a); // std::cout<<"err of szPosition TTrack:"<<err<<endl; if (err) return -2; return 0; }

int TTrack::szPosition (  TMLink &  link  )  const

calculates arc length and z for a stereo hit. { 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.); //...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; 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(); // std::cout<<"TTrack::wire direction:"<<h.wire()->direction()<<endl; 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); // std::cout<<"sz,Track::vmag:"<<vmag<<", helix_r:"<<r<<", wv:"<<wv<<", d:"<<sqrt(d2)<<endl; //...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(); // std::cout<<"X.z():"<<X.z()<<endl; // std::cout<<"szPosition::z(0),z(1),and leftRight:"<<z[0]<<","<<z[1]<<","<<link.leftRight()<<endl; // std::cout<<"szPosition::wire backwardPosition and forwardPosition:"<< h.wire()->backwardPosition().z()<<","<<h.wire()->forwardPosition().z()<<endl; #ifdef TRKRECO_DEBUG_DETAIL // std::cout << " l0, l1 = " << l[0] << ", " << l[1] << std::endl; // std::cout << " z0, z1 = " << z[0] << ", " << z[1] << std::endl; #endif //...Check z position... //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; }

unsigned TTrackBase::testByApproach (  const AList< TMLink > &  list, double  sigma )  const [inherited]

unsigned TTrackBase::testByApproach (  const TMLink &  list, double  sigma )  const [inherited]

returns # of good hits to be appended.

unsigned TTrackBase::testByApproach (  const AList< TMLink > &  list, double  sigma )  const [inherited]

{ #ifdef TRKRECO_DEBUG_DETAIL std::cout << " TTrackBase::testByApproach ... sigma=" << maxSigma << std::endl; #endif unsigned nOK = 0; unsigned n = list.length(); for (unsigned i = 0; i < n; i++) { TMLink & l = * list[i]; nOK += testByApproach(l, maxSigma); } return nOK; }

unsigned TTrackBase::testByApproach (  const TMLink &  list, double  sigma )  const [inherited]

returns # of good hits to be appended. { #ifdef TRKRECO_DEBUG_DETAIL std::cout << " TTrackBase::testByApproach ... sigma=" << maxSigma << std::endl; #endif TMLink l = link; //...Calculate closest approach... int err = approach(l); if (err < 0) return 0; //...Calculate sigma... float distance = l.distance(); float diff = fabs(distance - l.hit()->drift()); float sigma = diff / l.hit()->dDrift(); l.pull(sigma * sigma); //...For debug... #ifdef TRKRECO_DEBUG_DETAIL std::cout << " sigma=" << sigma; std::cout << ",dist=" << distance; std::cout << ",diff=" << diff << ","; if (sigma < maxSigma) std::cout << "ok,"; else std::cout << "X,"; l.dump("mc"); #endif //...Make sigma cut... if (sigma < maxSigma) return 1; return 0; }

unsigned TTrack::type (  void   )  const [virtual]

returns type. Definition is depending on an object type. Reimplemented from TTrackBase.

unsigned TTrack::type (  void   )  const [inline, virtual]

returns type. Definition is depending on an object type. Reimplemented from TTrackBase. { return defineType(); }

void TTrackBase::update (  void   )  const [inherited]

update cache. Reimplemented in TSegment, TSegment0, TSegment, and TSegment0.

void TTrackBase::update (  void   )  const [inherited]

update cache. Reimplemented in TSegment, TSegment0, TSegment, and TSegment0. { _cores.removeAll(); unsigned n = _links.length(); for (unsigned i = 0; i < n; i++) { TMLink * l = _links[i]; const TMDCWireHit & h = * l->hit(); if (h.state() & WireHitInvalidForFit) continue; if (! (h.state() & WireHitFittingValid)) continue; _cores.append(l); } _updated = true; }

---

Friends And Related Function Documentation

Refit [friend]

TBuilder [friend]

TBuilder0 [friend]

TBuilderCosmic [friend]

TBuilderCurl [friend]

TCosmicFitter [friend]

Reimplemented from TTrackBase.

TCurlFinder [friend]

THelixFitter [friend]

Reimplemented from TTrackBase.

TPMCurlFinder [friend]

TrkReco [friend]

TTrackManager [friend]

---

Member Data Documentation

AList<TMLink> TTrack::_associateHits [private]

AList<TMLink> TTrack::_associateHits [private]

double TTrack::_charge [private]

double TTrack::_chi2 [private]

TTrack* TTrack::_daughter [private]

TTrack* TTrack::_daughter [private]

AList<TMLink> TTrack::_finalHits [private]

AList<TMLink> TTrack::_finalHits [private]

bool TTrackBase::_fitted [mutable, protected, inherited]

bool TTrackBase::_fittedWithCathode [mutable, protected, inherited]

const THelixFitter TTrack::_fitter = THelixFitter("TTrack Default Helix Fitter") [static, private]

Reimplemented from TTrackBase.

Helix* const TTrack::_helix [private]

Helix* const TTrack::_helix [private]

AList<TMLink> TTrackBase::_links [protected, inherited]

AList<TMLink> TTrackBase::_links [protected, inherited]

TTrackMC* TTrackBase::_mc [protected, inherited]

TTrackMC* TTrackBase::_mc [protected, inherited]

TTrack* TTrack::_mother [private]

TTrack* TTrack::_mother [private]

std::string TTrack::_name [private]

unsigned TTrack::_ndf [private]

AList<TSegment> TTrack::_segments [private]

AList<TSegment> TTrack::_segments [private]

unsigned TTrack::_state [private]

unsigned TTrack::_type [mutable, private]

IMagneticFieldSvc* TTrack::m_pmgnIMF [private]

IMagneticFieldSvc* TTrack::m_pmgnIMF [private]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/TrkReco/TrkReco/TTrack.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Reconstruction/TrkReco/TrkReco-00-06-12/TrkReco/TTrack.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Reconstruction/TrkReco/TrkReco-00-06-12/src/TTrack.cxx

---