TSegment0 Class Reference

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

#include <TSegment0.h>

Inheritance diagram for TSegment0:

Public Member Functions
	TSegment0 ()
	Constructor.
	TSegment0 (const AList< TMLink > &)
virtual	~TSegment0 ()
	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 HepVector3D &	direction (void) const
	returns direction.
double	distance (const TSegment0 &) 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.
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< TSegment0 >	split (void) const
	returns a list of sub TSegments in this cluster. If cluster type is 1, 2, or 7, no cluster is returned.
int	solveDualHits (void)
double	duality (void) const
AList< TTrack > &	tracks (void)
AList< TSegment0 > &	innerLinks (void)
Public Member Functions inherited from TTrackBase
	TTrackBase ()
	Constructor.
	TTrackBase (const AList< TMLink > &links)
	Constructor.
virtual	~TTrackBase ()
	Destructor.
virtual unsigned	type (void) const
	returns type. Definition is depending on an object class.
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 41 of file TSegment0.h.

Constructor & Destructor Documentation

TSegment0() [1/2]

TSegment0::TSegment0

(

)

Constructor.

Definition at line 22 of file TSegment0.cxx.

: TTrackBase(),
  _innerWidth(0),
  _outerWidth(0),
  _nLayer(0),
  _clusterType(0),
  _duality(0.),
  _nDual(0),
  _angle(0.) {
    _fitted = false;
}

TSegment0() [2/2]

TSegment0::TSegment0

(

const AList< TMLink > &

a

)

Definition at line 34 of file TSegment0.cxx.

: TTrackBase(a),
  _innerWidth(0),
  _outerWidth(0),
  _nLayer(0),
  _clusterType(0),
  _duality(0.),
  _nDual(0),
  _angle(0.) {
    _links.sort(SortByWireId);
    _fitted = false;
}

~TSegment0()

TSegment0::~TSegment0 ( )

virtual

Destructor.

Definition at line 47 of file TSegment0.cxx.

                      {
}

Member Function Documentation

clusterType()

unsigned TSegment0::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 214 of file TSegment0.h.

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

Referenced by dump(), TConformalFinder0::findClusters2(), and split().

direction()

const HepVector3D & TSegment0::direction ( void ) const

inline

returns direction.

Definition at line 179 of file TSegment0.h.

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

distance() [1/2]

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

Definition at line 143 of file TSegment0.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/2]

double TSegment0::distance

(

const TSegment0 &

c

)

const

calculates distance between two clusters. Smaller value indicates closer.

Definition at line 131 of file TSegment0.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().

duality()

double TSegment0::duality ( void ) const

inline

Definition at line 222 of file TSegment0.h.

                             {
    return _duality;
}

dump()

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

Definition at line 51 of file TSegment0.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) TTrackBase::dump(msg, pre);
}

Referenced by TConformalFinder0::findClusters(), TConformalFinder0::findClusters2(), THistogram::segments0(), solveDualHits(), and superLayerId().

innerLinks()

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

inline

Definition at line 261 of file TSegment0.h.

                          {
    return _innerLinks;
}

innerMostLayer()

unsigned TSegment0::innerMostLayer ( void ) const

inline

returns inner most layer.

Definition at line 200 of file TSegment0.h.

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

innerWidth()

unsigned TSegment0::innerWidth ( void ) const

inline

returns inner width.

Definition at line 186 of file TSegment0.h.

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

objectType()

unsigned TSegment0::objectType ( void ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 228 of file TSegment0.h.

                                {
    return Segment;
}

outerMostLayer()

unsigned TSegment0::outerMostLayer ( void ) const

inline

returns outer most layer.

Definition at line 207 of file TSegment0.h.

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

outerWidth()

unsigned TSegment0::outerWidth ( void ) const

inline

returns outer width.

Definition at line 193 of file TSegment0.h.

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

position()

const HepPoint3D & TSegment0::position ( void ) const

inline

returns position.

Definition at line 172 of file TSegment0.h.

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

Referenced by distance().

rangeX()

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

returns Range of x-coordinate of TMLinks.

Definition at line 155 of file TSegment0.cxx.

                                              {
#ifdef TRKRECO_DEBUG_DETAIL
    if (min > max) {
    std::cout << "TSegment0::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);
}

solveDualHits()

int TSegment0::solveDualHits

(

void

)

Definition at line 624 of file TSegment0.cxx.

                             {
    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 - h0->drift() - h1->drift());
        if (distance > 0.5) duals.append(list);
#ifdef TRKRECO_DEBUG_DETAIL
//      h0->dump();
//      h1->dump();
//      std::cout << "duality distance = " << distance << std::endl;
//      std::cout << "i = " << i << std::endl;
#endif
        }
    }
    else if (list.length() == 0) {
        continue;
    }
#ifdef TRKRECO_DEBUG_DETAIL
    else {
        std::cout << "TSegment0::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);
    }
    return n;
}

Referenced by TConformalFinder0::findClusters(), TConformalFinder0::findClusters2(), and THistogram::segments0().

split()

AList< TSegment0 > TSegment0::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 281 of file TSegment0.cxx.

                           {
    AList<TSegment0> 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();
    if (_nDual > 2 && _duality > 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();
 
    //...Avoid splitting of X or parallel...(future implementation)...
    else if (t > 4) return splitComplicated();
 
    //...A or V...
    return splitAV();
}

Referenced by TConformalFinder0::findClusters(), and THistogram::segments0().

superLayerId()

unsigned TSegment0::superLayerId ( void ) const

inline

returns super layer id.

Definition at line 234 of file TSegment0.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 << "TSegment0::superLayerId !!! strange segment found";
        std::cout << std::endl;
        dump();
        break;
        }
    }
    }
#endif
    return id;
}

Referenced by superLayerId().

tracks()

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

inline

Definition at line 255 of file TSegment0.h.

                      {
    return _tracks;
}

---

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

• TSegment0.h
• TSegment0.cxx