TSegment Class Reference

A class to relate TMDCWireHit and TTrack objects. More...

#include <TSegment.h>

Inheritance diagram for TSegment:

Public Member Functions
	TSegment ()
	Constructor.
	TSegment (const AList< TMLink > &)
virtual	~TSegment ()
	Destructor.
virtual unsigned	objectType (void) const
	returns type.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
unsigned	superLayerId () const
	returns super layer id.
const HepPoint3D &	position (void) const
	returns position.
const HepPoint3D &	outerPosition (void) const
const HepPoint3D &	lineTsf (void) const
	return line of Tsf for pos and dir
const HepPoint3D &	lineTsf (const HepPoint3D &)
const HepVector3D &	direction (void) const
	returns direction.
double	distance (const TSegment &) const
	calculates distance between two clusters. Smaller value indicates closer.
double	distance (const HepPoint3D &, const HepVector3D &) const
Range	rangeX (double min, double max) const
	returns Range of x-coordinate of TMLinks.
const AList< TMLink > &	inners (void) const
	returns TMLinks in inner/outer-most layer.
const AList< TMLink > &	outers (void) const
const TMLink &	center (void) const
	returns a TMLink which is the closest to the center.
unsigned	width (void) const
	returns width.
unsigned	innerWidth (void) const
	returns inner width.
unsigned	outerWidth (void) const
	returns outer width.
unsigned	innerMostLayer (void) const
	returns inner most layer.
unsigned	outerMostLayer (void) const
	returns outer most layer.
unsigned	clusterType (void) const
	returns cluster type. 0:empty, 1:short line, 2:long line, 3:V shage(from outside), 4:A shape, 5:X shape, 6:parallel, 7:complicated.
AList< TSegment >	split (void) const
	returns a list of sub TSegments in this cluster. If cluster type is 1, 2, or 7, no cluster is returned.
AList< TSegment >	splitTsf (AList< TMLink > &)
int	solveDualHits (void)
double	duality (void) const
void	solveLR (void)
	solve LR of hit in TSF.
void	solveThreeHits (void)
AList< TTrack > &	tracks (void)
AList< TSegment > &	innerLinks (void)
const AList< TSegment > &	innerLinks (void) const
AList< TSegment > &	outerLinks (void)
const AList< TSegment > &	outerLinks (void) const
unsigned	state (void) const
unsigned	state (unsigned)
	TSegment ()
	Constructor.
	TSegment (const AList< TMLink > &)
virtual	~TSegment ()
	Destructor.
virtual unsigned	objectType (void) const
	returns type.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
unsigned	superLayerId () const
	returns super layer id.
const HepPoint3D &	position (void) const
	returns position.
const HepPoint3D &	outerPosition (void) const
const HepPoint3D &	lineTsf (void) const
	return line of Tsf for pos and dir
const HepPoint3D &	lineTsf (const HepPoint3D &)
const HepVector3D &	direction (void) const
	returns direction.
double	distance (const TSegment &) const
	calculates distance between two clusters. Smaller value indicates closer.
double	distance (const HepPoint3D &, const HepVector3D &) const
Range	rangeX (double min, double max) const
	returns Range of x-coordinate of TMLinks.
const AList< TMLink > &	inners (void) const
	returns TMLinks in inner/outer-most layer.
const AList< TMLink > &	outers (void) const
const TMLink &	center (void) const
	returns a TMLink which is the closest to the center.
unsigned	width (void) const
	returns width.
unsigned	innerWidth (void) const
	returns inner width.
unsigned	outerWidth (void) const
	returns outer width.
unsigned	innerMostLayer (void) const
	returns inner most layer.
unsigned	outerMostLayer (void) const
	returns outer most layer.
unsigned	clusterType (void) const
	returns cluster type. 0:empty, 1:short line, 2:long line, 3:V shage(from outside), 4:A shape, 5:X shape, 6:parallel, 7:complicated.
AList< TSegment >	split (void) const
	returns a list of sub TSegments in this cluster. If cluster type is 1, 2, or 7, no cluster is returned.
AList< TSegment >	splitTsf (AList< TMLink > &)
int	solveDualHits (void)
double	duality (void) const
void	solveLR (void)
	solve LR of hit in TSF.
void	solveThreeHits (void)
AList< TTrack > &	tracks (void)
AList< TSegment > &	innerLinks (void)
const AList< TSegment > &	innerLinks (void) const
AList< TSegment > &	outerLinks (void)
const AList< TSegment > &	outerLinks (void) const
unsigned	state (void) const
unsigned	state (unsigned)
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).
	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).

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

Detailed Description

A class to relate TMDCWireHit and TTrack objects.

Definition at line 43 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

Constructor & Destructor Documentation

TSegment() [1/4]

TSegment::TSegment

(

)

Constructor.

Definition at line 30 of file TSegment.cxx.

: TTrackBase(),
  _innerWidth(0),
  _outerWidth(0),
  _nLayer(0),
  _clusterType(0),
  _duality(0.),
  _nDual(0),
  _angle(0.),
  _state(0),
  _lineTsf(0.,0.,0.) {
//    _times = 4.0;
//    if (_links[0]->wire()->stereo()) _times = 7.0;
    _times[0] = 7;
    _times[1] = 7;
    _times[2] = 4;
    _times[3] = 4;
    _times[4] = 4;
    _times[5] = 5;
    _times[6] = 5;
    _times[7] = 5;
    _times[8] = 5;
    _times[9] = 4;
    _times[10] = 4;
    _fitted = false;
}

TSegment() [2/4]

TSegment::TSegment

(

const AList< TMLink > &

a

)

Definition at line 57 of file TSegment.cxx.

: TTrackBase(a),
  _innerWidth(0),
  _outerWidth(0),
  _nLayer(0),
  _clusterType(0),
  _duality(0.),
  _nDual(0),
  _angle(0.),
  _state(0),
  _lineTsf(0.,0.,0.) {
    _links.sort(SortByWireId);
//    _times = 4.0;
//    if (_links[0]->wire()->stereo()) _times = 7.0;
    _times[0] = 7;
    _times[1] = 7;
    _times[2] = 4;
    _times[3] = 4;
    _times[4] = 4;
    _times[5] = 5;
    _times[6] = 5;
    _times[7] = 5;
    _times[8] = 5;
    _times[9] = 4;
    _times[10] = 4;
    _fitted = false;
}  //the above two innitial function are temporary for tsf!!!

~TSegment() [1/2]

TSegment::~TSegment ( )

virtual

Destructor.

Definition at line 112 of file TSegment.cxx.

                    {
}

TSegment() [3/4]

TSegment::TSegment

(

)

Constructor.

TSegment() [4/4]

TSegment::TSegment

(

const AList< TMLink > &

)

~TSegment() [2/2]

virtual TSegment::~TSegment ( )

virtual

Destructor.

Member Function Documentation

center() [1/2]

const TMLink & TSegment::center

(

void

)

const

returns a TMLink which is the closest to the center.

Referenced by rangeX(), and solveLR().

center() [2/2]

const TMLink & TSegment::center

(

void

)

const

returns a TMLink which is the closest to the center.

clusterType() [1/2]

unsigned TSegment::clusterType ( void  ) const

inline

returns cluster type. 0:empty, 1:short line, 2:long line, 3:V shage(from outside), 4:A shape, 5:X shape, 6:parallel, 7:complicated.

Definition at line 319 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                                {
    if (! nLinks()) return 0;
    if (_clusterType == 0) updateType();
    return _clusterType;
}

Referenced by dump(), and split().

clusterType() [2/2]

unsigned TSegment::clusterType

(

void

)

const

returns cluster type. 0:empty, 1:short line, 2:long line, 3:V shage(from outside), 4:A shape, 5:X shape, 6:parallel, 7:complicated.

direction() [1/2]

const HepVector3D & TSegment::direction ( void  ) const

inline

returns direction.

Definition at line 284 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                              {
    if (! _fitted) update();
    return _direction;
}

direction() [2/2]

const HepVector3D & TSegment::direction

(

void

)

const

returns direction.

distance() [1/4]

double TSegment::distance	(	const HepPoint3D &	p,
		const HepVector3D &	v
	)		const

Definition at line 249 of file TSegment.cxx.

                                                                    {
    HepVector3D dir = _position - p;
    if (dir.x() > M_PI) dir.setX(dir.x() - 2. * M_PI);
    else if (dir.x() < - M_PI) dir.setX(2. * M_PI + dir.x());
 
    float radial = fabs(v.unit().dot(dir));
    float radial2 = radial * radial;
 
    return (dir.mag2() - radial2) > 0.0 ? sqrt(dir.mag2() - radial2) : 0.;
}

distance() [2/4]

double TSegment::distance	(	const HepPoint3D &	,
		const HepVector3D &
	)		const

distance() [3/4]

double TSegment::distance

(

const TSegment &

c

)

const

calculates distance between two clusters. Smaller value indicates closer.

Definition at line 237 of file TSegment.cxx.

                                           {
    HepVector3D dir = c.position() - _position;
    if (dir.x() > M_PI) dir.setX(dir.x() - 2. * M_PI);
    else if (dir.x() < - M_PI) dir.setX(2. * M_PI + dir.x());
 
    float radial = fabs(_direction.dot(dir));
    float radial2 = radial * radial;
 
    return (dir.mag2() - radial2) > 0.0 ? sqrt(dir.mag2() - radial2) : 0.;
}

Referenced by solveDualHits().

distance() [4/4]

double TSegment::distance

(

const TSegment &

)

const

calculates distance between two clusters. Smaller value indicates closer.

duality() [1/2]

double TSegment::duality ( void  ) const

inline

Definition at line 327 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                            {
    return _duality;
}

duality() [2/2]

double TSegment::duality

(

void

)

const

dump() [1/2]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

Definition at line 116 of file TSegment.cxx.

                                                                 {
    if (! _fitted) update();
    bool def = false;
    if (msg == "") def = true;
 
    if (def || msg.find("cluster") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "#links=" << _links.length();
    std::cout << "(" << _innerWidth << "," << _outerWidth << ":";
    std::cout << clusterType() << ")," << _nDual << "," << _duality << ",";
    std::cout << _angle << std::endl;
    }
    if (def || msg.find("vector") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "pos" << _position << "," << "dir" << _direction;
    std::cout << std::endl;
    }
    if (def || msg.find("state") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "state=" << _state << std::cout << std::endl;
    }
    if (def || msg.find("link") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "unique links=" << NUniqueLinks(* this);
    std::cout << ",major links=" << NMajorLinks(* this);
    std::cout << ",branches=" << NLinkBranches(* this);
    std::cout << std::endl;
    }
    if (! def) TTrackBase::dump(msg, pre);
}

Referenced by THistogram::segments(), solveDualHits(), and superLayerId().

dump() [2/2]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

innerLinks() [1/4]

AList< TSegment > & TSegment::innerLinks ( void  )

inline

Definition at line 366 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                         {
    return _innerLinks;
}

innerLinks() [2/4]

AList< TSegment > & TSegment::innerLinks

(

void

)

innerLinks() [3/4]

const AList< TSegment > & TSegment::innerLinks ( void  ) const

inline

Definition at line 372 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                               {
    return _innerLinks;
}

innerLinks() [4/4]

const AList< TSegment > & TSegment::innerLinks

(

void

)

const

innerMostLayer() [1/2]

unsigned TSegment::innerMostLayer ( void  ) const

inline

returns inner most layer.

Definition at line 305 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                                   {
    if (! _fitted) update();
    return _innerMostLayer;
}

Referenced by width().

innerMostLayer() [2/2]

unsigned TSegment::innerMostLayer

(

void

)

const

returns inner most layer.

inners() [1/2]

const AList< TMLink > & TSegment::inners ( void  ) const

inline

returns TMLinks in inner/outer-most layer.

Definition at line 263 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                           {
    if (! _fitted) update();
    return _inners;
}

inners() [2/2]

const AList< TMLink > & TSegment::inners

(

void

)

const

returns TMLinks in inner/outer-most layer.

innerWidth() [1/2]

unsigned TSegment::innerWidth ( void  ) const

inline

returns inner width.

Definition at line 291 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                               {
    if (! _fitted) update();
    return _innerWidth;
}

innerWidth() [2/2]

unsigned TSegment::innerWidth

(

void

)

const

returns inner width.

lineTsf() [1/4]

const HepPoint3D & TSegment::lineTsf ( const HepPoint3D &  a )

inline

Definition at line 409 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                                     {
    return _lineTsf = a;
}

lineTsf() [2/4]

const HepPoint3D & TSegment::lineTsf

(

const HepPoint3D &

)

lineTsf() [3/4]

const HepPoint3D & TSegment::lineTsf ( void  ) const

inline

return line of Tsf for pos and dir

Definition at line 403 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                            {
    return _lineTsf;
}

lineTsf() [4/4]

const HepPoint3D & TSegment::lineTsf

(

void

)

const

return line of Tsf for pos and dir

objectType() [1/2]

unsigned TSegment::objectType ( void  ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 333 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                               {
    return Segment;
}

objectType() [2/2]

virtual unsigned TSegment::objectType ( void  ) const

virtual

returns type.

Reimplemented from TTrackBase.

outerLinks() [1/4]

AList< TSegment > & TSegment::outerLinks ( void  )

inline

Definition at line 378 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                         {
    return _outerLinks;
}

Referenced by OuterMostUniqueLink().

outerLinks() [2/4]

AList< TSegment > & TSegment::outerLinks

(

void

)

outerLinks() [3/4]

const AList< TSegment > & TSegment::outerLinks ( void  ) const

inline

Definition at line 384 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                               {
    return _outerLinks;
}

outerLinks() [4/4]

const AList< TSegment > & TSegment::outerLinks

(

void

)

const

outerMostLayer() [1/2]

unsigned TSegment::outerMostLayer ( void  ) const

inline

returns outer most layer.

Definition at line 312 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                                   {
    if (! _fitted) update();
    return _outerMostLayer;
}

Referenced by width().

outerMostLayer() [2/2]

unsigned TSegment::outerMostLayer

(

void

)

const

returns outer most layer.

outerPosition() [1/2]

const HepPoint3D & TSegment::outerPosition

(

void

)

const

outerPosition() [2/2]

const HepPoint3D & TSegment::outerPosition

(

void

)

const

outers() [1/2]

const AList< TMLink > & TSegment::outers ( void  ) const

inline

Definition at line 270 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                           {
    if (! _fitted) update();
    return _outers;
}

Referenced by solveDualHits().

outers() [2/2]

const AList< TMLink > & TSegment::outers

(

void

)

const

outerWidth() [1/2]

unsigned TSegment::outerWidth ( void  ) const

inline

returns outer width.

Definition at line 298 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                               {
    if (! _fitted) update();
    return _outerWidth;
}

outerWidth() [2/2]

unsigned TSegment::outerWidth

(

void

)

const

returns outer width.

position() [1/2]

const HepPoint3D & TSegment::position ( void  ) const

inline

returns position.

Definition at line 277 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                             {
    if (! _fitted) update();
    return _position;
}

Referenced by distance().

position() [2/2]

const HepPoint3D & TSegment::position

(

void

)

const

returns position.

rangeX() [1/2]

Range TSegment::rangeX	(	double	min,
		double	max
	)		const

returns Range of x-coordinate of TMLinks.

Definition at line 261 of file TSegment.cxx.

                                             {
#ifdef TRKRECO_DEBUG_DETAIL
    if (min > max) {
    std::cout << "TSegment::range !!! bad arguments:min,max=";
    std::cout << min << "," << max << std::endl;
    }
#endif
 
    unsigned n = _links.length();
    if (n == 0) return Range(0., 0.);
 
    //...Search for a center...
    bool found = false;
    double center;
    for (unsigned i = 0; i < n; i++) {
    double x = _links[i]->position().x();
    if (x < min || x > max) continue;
    center = x;
    found = true;
    break;
    }
    if (! found) return Range(0., 0.);
 
#ifdef TRKRECO_DEBUG_DETAIL
//     std::cout << "    center=" << center << std::endl;
#endif
 
    double distanceR = 0.;
    double distanceL = 0.;
    double distanceMax = max - min;
    for (unsigned i = 0; i < n; i++) {
    double x = _links[i]->position().x();
    if (x < min || x > max) continue;
 
    double distance0, distance1;
    if (x > center) {
        distance0 = x - center;
        distance1 = distanceMax - distance0;
    }
    else {
        distance1 = center - x;
        distance0 = distanceMax - distance1;
    }
 
    if (distance0 < distance1) {
        if (distance0 > distanceR) distanceR = distance0;
    }
    else {
        if (distance1 > distanceL) distanceL = distance1;
    }
 
#ifdef TRKRECO_DEBUG_DETAIL
//  std::cout << "    ";
//  std::cout << _links[i]->wire()->layerId() << "-";
//  std::cout << _links[i]->wire()->localId() << ",";
//  std::cout << _links[i]->position().x();
//  std::cout << ",0,1=" << distance0 << "," << distance1;
//  std::cout << ",l,r=" << distanceL << "," << distanceR;
//  std::cout << std::endl;
#endif
    }
 
    double right = center + distanceR;
    double left = center - distanceL;
 
    return Range(left, right);
}

rangeX() [2/2]

Range TSegment::rangeX	(	double	min,
		double	max
	)		const

returns Range of x-coordinate of TMLinks.

solveDualHits() [1/2]

int TSegment::solveDualHits

(

void

)

Definition at line 736 of file TSegment.cxx.

                            {
    updateType();
    if (_clusterType == 0) return 0;
    if (_nDual == 0) return 0;
    update();
    return 0;
 
    updateType();
    if (_clusterType == 0) return 0;
    if (_nDual == 0) return 0;
 
    AList<TMLink> seeds;
    AList<TMLink> duals;
    for (unsigned i = _innerMostLayer; i <= _outerMostLayer; i++) {
    AList<TMLink> list = SameLayer(_links, i);
 
    if (list.length() == 1) {
        seeds.append(list[0]);
    }
    else if (list.length() == 2) {
        if (Width(list) > 1) {
        const TMDCWireHit * h0 = list[0]->hit();
        const TMDCWireHit * h1 = list[1]->hit();
        double distance = (h0->xyPosition() - h1->xyPosition()).mag();
        distance = fabs(distance - 
                list[0]->drift() - 
                list[1]->drift());
        if (distance > 0.5) duals.append(list);
#ifdef TRKRECO_DEBUG_DETAIL
        h0->dump("", "        ");
        h1->dump("", "        ");
        std::cout << "        lyr=" << i;
        std::cout << ", duality distance = " << distance << std::endl;
#endif
        }
    }
    else if (list.length() == 0) {
        continue;
    }
#ifdef TRKRECO_DEBUG_DETAIL
    else {
        std::cout << "TSegment::solveDualHits !!! this is not expected 2";
        std::cout << std::endl;
        this->dump("cluster hits mc", "    ");
    }
#endif
    }
 
    //...Solve them...
    if (seeds.length() < 2) return -1;
    AList<TMLink> outers = InOut(seeds);
    const HepPoint3D & p = outers[0]->xyPosition();
    HepVector3D v = (outers[1]->xyPosition() - p).unit();
    unsigned n = duals.length() / 2;
    for (unsigned i = 0; i < n; i++) {
    TMLink * t0 = duals[i * 2 + 0];
    TMLink * t1 = duals[i * 2 + 1];
    HepVector3D x0 = t0->xyPosition() - p;
    HepVector3D x1 = t1->xyPosition() - p;
    double d0 = (x0 - (x0.dot(v) * v)).mag();
    double d1 = (x1 - (x1.dot(v) * v)).mag();
    if (d0 < d1) _links.remove(t1);
    else         _links.remove(t0);
    }
    update();
    return n;
}

Referenced by THistogram::segments().

solveDualHits() [2/2]

int TSegment::solveDualHits

(

void

)

solveLR() [1/2]

void TSegment::solveLR

(

void

)

solve LR of hit in TSF.

Definition at line 1604 of file TSegment.cxx.

                      {
    unsigned n = _links.length();
    //initial...
    if (n <= 3) {
      solveThreeHits();
      return;
    }
    for (int i = 0; i < n; ++i) {
    TMLink *l = _links[i];
    if (l->hit()->state() & WireHitPatternRight){
        unsigned newState = l->hit()->state()&(~WireHitPatternRight);
        l->hit()->state(newState);
    }
    if (l->hit()->state() & WireHitPatternLeft){
        unsigned newState = l->hit()->state()&(~WireHitPatternLeft);
        l->hit()->state(newState);
    }
    }
 
    HepPoint3D originCon(0., 0., 0.);
    HepPoint3D dirZ(0., 0., 1.);
    HepPoint3D center = closestPoint(originCon, _lineTsf);
    HepPoint3D v1 = (center - originCon).unit();
    HepPoint3D v2 = dirZ.cross(v1);
    for (int i = 0; i < n; ++i) {
//  cout<<"layerId: "<<_links[i]->wire()->layerId()<<"  localId: "<<_links[i]->wire()->localId()<<endl;
    const HepPoint3D & p = _links[i]->positionD();
    unsigned state = _links[i]->hit()->state();
 
    double cosA = (p - center).unit().dot(v1.unit());
    double cosB = (p - center).unit().dot(v2.unit());
    if (cosA*cosB < 0) state |= WireHitPatternLeft;
    else state |= WireHitPatternRight;
 
    _links[i]->hit()->state(state);
    }
}

solveLR() [2/2]

void TSegment::solveLR

(

void

)

solve LR of hit in TSF.

solveThreeHits() [1/2]

void TSegment::solveThreeHits

(

void

)

Definition at line 1658 of file TSegment.cxx.

                             {
    unsigned n = _links.length();
 
    for (unsigned i = 0; i < n; i++) {
    const TMDCWireHit * h = _links[i]->hit();
    const TMDCWire * w = h->wire();
    unsigned state = h->state();
 
    //...Cache pointers to a neighbor...
    const TMDCWire * neighbor[6];
    for (unsigned j = 0; j < 6; j++) neighbor[j] = w->neighbor(j);
 
    //...Decide hit pattern...
    unsigned pattern = 0;
    for (unsigned j = 0; j < 6; j++) {
        if (neighbor[j])
        if (neighbor[j]->hit())
            pattern += (1 << j);
    }
    state |= (pattern << WireHitNeighborHit);
 
    //...Solve LR locally...
    if ((pattern == 34) || (pattern == 42) ||
        (pattern == 40) || (pattern == 10) ||
        (pattern == 35) || (pattern == 50))
        state |= WireHitPatternRight;
    else if ((pattern == 17) || (pattern == 21) ||
         (pattern == 20) || (pattern ==  5) ||
         (pattern == 19) || (pattern == 49))
        state |= WireHitPatternLeft;
 
    //...Store it...
    h->state(state);
    }
}

Referenced by solveLR().

solveThreeHits() [2/2]

void TSegment::solveThreeHits

(

void

)

split() [1/2]

AList< TSegment > TSegment::split

(

void

)

const

returns a list of sub TSegments in this cluster. If cluster type is 1, 2, or 7, no cluster is returned.

Definition at line 388 of file TSegment.cxx.

                          {
    AList<TSegment> list;
 
    //...Do not split if cluster type is 1, 2, or 7...
    unsigned t = clusterType();
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    ... splitting : type=" << t << std::endl;
#endif 
    if (t == 0) return list;
    else if (t == 2) {
    // beta 5
    //if (_nDual > 2 && _duality > 0.7 && _angle > 0.7) return splitDual();
 
    // 1.67g
    // if (_nDual > 2 && _duality > 0.7) return splitDual();
 
    if (_nDual > 2 && _duality > 0.7 && _angle > 0.7) return splitDual();
    return list;
    }
    else if (t == 1) return list;
    else if (t == 7) return list;
 
    //...Parallel...
    else if (t == 6) return splitParallel();
    // else if (t == 6) return list;
 
    //...Avoid splitting of X or parallel...(future implementation)...
    else if (t > 4) return splitComplicated();
 
    //...A or V...
    return splitAV();
}

Referenced by THistogram::segments().

split() [2/2]

AList< TSegment > TSegment::split

(

void

)

const

returns a list of sub TSegments in this cluster. If cluster type is 1, 2, or 7, no cluster is returned.

splitTsf() [1/2]

AList< TSegment > TSegment::splitTsf

(

AList< TMLink > &

seedNeighbors

)

Definition at line 1134 of file TSegment.cxx.

                                                {
    //get links in each layer.
    AList<TSegment> list;
    AList<TMLink> links[4];   //links in each local layer.
    AList<TMLink> seeds2;     //links in the outer layer.
    AList<TMLink> exTsf;      //links in other local layers.
    TMLink *seed;             //link in the mostinner layer.
 
    for (unsigned i = 0; i < _links.length(); ++i) {
    TMLink * l = _links[i];
    links[l->wire()->localLayerId()].append(l);
    }
 
    //prepare for segment finding.
    if (links[0].length() == 0) {  //find in 234.
    if (links[3].length() == 0) return list;
    if (links[1].length() != 1) {
        cout<<"wrong in tsf, TSegment::splitTsf ...1"<<endl;
        return list;
    }
    seed = links[1][0];
    seeds2.append(links[3]);
    exTsf.append(links[2]);
    }
    else if (links[0].length() == 1) { //find in 1234, 124, 134, 123.
    seed = links[0][0];
    if (links[3].length() > 0) {     //1..4
        seeds2.append(links[3]);
        exTsf.append(links[1]);
        exTsf.append(links[2]);
    }
    else if (links[2].length() > 0) { //1..3
        seeds2.append(links[2]);
        exTsf.append(links[1]);
    }
    else return list;
    }
    else cout<<"wrong in tsf, TSegment::splitTsf ...2"<<endl;
    exTsf.append(seeds2);  //add the outer layer...
 
    //find segments
    for (unsigned i = 0; i < seeds2.length(); ++i) {
    if (seed->tsfTag() > 0 && seeds2[i]->tsfTag() > 0) continue;
    AList<TSegment> segments;
    HepPoint3D line[4];
    fitLine(seed, seeds2[i], line);
    segments = appendByLine(seed, seeds2[i], seedNeighbors, exTsf, line);
    list.append(segments);
    segments.removeAll();
    }
    if (seed->wire()->localLayerId() == 0) {   //for 123
    exTsf.removeAll();
    seeds2.removeAll();
    exTsf.append(links[1]);
    exTsf.append(links[2]);   //add the outer layer...
    for (int k = 0; k < links[2].length(); ++k){
        if (links[2][k]->tsfTag() == 0) seeds2.append(links[2][k]);
    }
    for (unsigned i = 0; i < seeds2.length(); ++i) {
        if (seed->tsfTag() > 0 && seeds2[i]->tsfTag() > 0) continue;
        AList<TSegment> segments2;
        HepPoint3D line2[4];
        fitLine(seed, seeds2[i], line2);
        segments2 = appendByLine(seed, seeds2[i], seedNeighbors, exTsf, line2);
        list.append(segments2);
        segments2.removeAll();
    }
    }
 
    return list;
}

Referenced by TMDCTsf::segments().

splitTsf() [2/2]

AList< TSegment > TSegment::splitTsf

(

AList< TMLink > &

)

state() [1/4]

unsigned TSegment::state ( unsigned  a )

inline

Definition at line 396 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                          {
    return _state = a;
}

state() [2/4]

unsigned TSegment::state

(

unsigned

)

state() [3/4]

unsigned TSegment::state ( void  ) const

inline

Definition at line 390 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                          {
    return _state;
}

Referenced by solveLR(), and solveThreeHits().

state() [4/4]

unsigned TSegment::state

(

void

)

const

superLayerId() [1/2]

unsigned TSegment::superLayerId ( void  ) const

inline

returns super layer id.

Definition at line 339 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                                 {
    unsigned id = (links())[0]->wire()->superLayerId();
#ifdef TRKRECO_DEBUG
    {
    const AList<TMLink> & list = links();
    unsigned n = list.length();
    for (unsigned i = 1; i < n; i++) {
        if (list[i]->hit()->wire()->superLayerId() != id) {
        std::cout << "TSegment::superLayerId !!! strange segment found";
        std::cout << std::endl;
        dump();
        break;
        }
    }
    }
#endif
    return id;
}

Referenced by superLayerId().

superLayerId() [2/2]

unsigned TSegment::superLayerId

(

)

const

returns super layer id.

tracks() [1/2]

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

inline

Definition at line 360 of file InstallArea/include/TrkReco/TrkReco/TSegment.h.

                     {
    return _tracks;
}

Referenced by TBuilder::buildRphi(), TBuilder::buildStereo(), and TBuilder::buildStereoNew().

tracks() [2/2]

AList< TTrack > & TSegment::tracks

(

void

)

width() [1/2]

unsigned TSegment::width

(

void

)

const

returns width.

Definition at line 1101 of file TSegment.cxx.

                          {
    unsigned maxWidth = 0;
    for (unsigned i = innerMostLayer(); i <= outerMostLayer(); i++) {
    AList<TMLink> tmp = SameLayer(links(), i);
    unsigned w = Width(tmp);
    if (w > maxWidth) maxWidth = w;
    }
    return maxWidth;
}

width() [2/2]

unsigned TSegment::width

(

void

)

const

returns width.

---

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

• InstallArea/include/TrkReco/TrkReco/TSegment.h
• Reconstruction/TrkReco/TrkReco-00-08-59-patch4-slc6tag/TrkReco/TSegment.h
• TSegment.cxx