T3DLine Class Reference

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

#include <T3DLine.h>

Inheritance diagram for T3DLine:

List of all members.

Public Member Functions
	T3DLine ()
	Constructors.
	T3DLine (const TTrack &)
	T3DLine (const T3DLine &)
virtual	~T3DLine ()
	Destructor.
double	dr (void) const
	Track parameters.
double	phi0 (void) const
double	dz (void) const
double	tanl (void) const
double	cosPhi0 (void) const
double	sinPhi0 (void) const
Vector	a (void) const
const SymMatrix &	Ea (void) const
	returns error matrix
const HepPoint3D &	pivot (void) const
	pivot position
Helix	helix (void) const
	approximated helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
unsigned	objectType (void) const
	returns object type
HepPoint3D	x (double) const
	returns position on 3D line
HepPoint3D	x0 (void) const
	returns 3D line component x(t)=x0 + t * k
HepVector3D	k (void) const
int	approach (TMLink &, bool sagCorrection=true) const
	calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.
int	approach_line (const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack) const
	caluculate closest points between a line and this track
int	approach_point (const HepPoint3D &, HepPoint3D &onTrack) const
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
void	set (const HepPoint3D &, double t_dr, double t_phi0, double t_dz, double t_tanl)
	set track parameters,pivot
Vector	a (const Vector &)
	set track parameters
const SymMatrix &	Ea (const SymMatrix &)
	set error matrix
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).
Protected Attributes
AList< TMLink >	_links
bool	_fitted
bool	_fittedWithCathode
TTrackMC *	_mc
Private Attributes
HepPoint3D	_pivot
double	_dr
double	_phi0
double	_dz
double	_tanl
SymMatrix	_Ea
double	_chi2
unsigned	_ndf
double	_cos_phi0
double	_sin_phi0
Static Private Attributes
static const T3DLineFitter	_fitter = T3DLineFitter("T3DLine Default fitter")
Friends
class	T3DLineFitter

---

Detailed Description

A class to represent a track in tracking.

Definition at line 50 of file T3DLine.h.

---

Constructor & Destructor Documentation

T3DLine::T3DLine

(

)

Constructors.

Definition at line 21 of file T3DLine.cxx.

References _cos_phi0, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _sin_phi0, and TTrackBase::fitter().

  : TTrackBase(),
    _dr(0),_phi0(0),_dz(0),_tanl(0),
    _chi2(0),_ndf(0),_Ea(SymMatrix(4,0)),
    _pivot(ORIGIN){

  _cos_phi0 = 1;
  _sin_phi0 = 0;

  //...Set a default fitter...
  fitter(& T3DLine::_fitter);

  _fitted = false;
  _fittedWithCathode = false;
}

T3DLine::T3DLine

(

const TTrack &

)

Definition at line 37 of file T3DLine.cxx.

References _cos_phi0, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _phi0, _sin_phi0, cos(), TTrackBase::fitter(), and sin().

  : TTrackBase((TTrackBase &) a),
    _dr(a.helix().dr()),
    _phi0(a.helix().phi0()),
    _dz(a.helix().dz()),
    _tanl(a.helix().tanl()),
    _chi2(0),_ndf(0),_Ea(SymMatrix(4,0)),
    _pivot(a.helix().pivot()){

  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);

  //...Set a default fitter...
  fitter(& T3DLine::_fitter);

  _fitted = false;
  _fittedWithCathode = false;
}

T3DLine::T3DLine

(

const T3DLine &

)

Definition at line 59 of file T3DLine.cxx.

References _cos_phi0, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _phi0, _sin_phi0, cos(), TTrackBase::fitter(), and sin().

  : TTrackBase((TTrackBase &) a),
    _dr(a.dr()),
    _phi0(a.phi0()),
    _dz(a.dz()),
    _tanl(a.tanl()),
    _chi2(a.chi2()),_ndf(a.ndf()),
    _Ea(a.Ea()),
    _pivot(a.pivot()){

  _cos_phi0 = cos(_phi0);
  _sin_phi0 = sin(_phi0);

  //...Set a default fitter...
  fitter(& T3DLine::_fitter);

  _fitted = false;
  _fittedWithCathode = false;
}

T3DLine::~T3DLine

(

)

[virtual]

Destructor.

Definition at line 56 of file T3DLine.cxx.

                 {
}

---

Member Function Documentation

Vector T3DLine::a

(

const Vector &

)

set track parameters

Definition at line 187 of file T3DLine.cxx.

References _cos_phi0, _dr, _dz, _phi0, _sin_phi0, _tanl, cos(), and sin().

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

Vector T3DLine::a

(

void

)

const

Definition at line 107 of file T3DLine.cxx.

References _dr, _dz, _phi0, and _tanl.

Referenced by helix().

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

void TTrackBase::append

(

const AList< TMLink > &

)

[inherited]

appends TMLinks.

Definition at line 384 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, TTrackBase::_updated, genRecEmupikp::i, and WireHitUsed.

                                          {
    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.

Definition at line 362 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, TTrackBase::_updated, TMLink::hit(), TMDCWire::name(), TMDCWireHit::state(), TMLink::wire(), and WireHitUsed.

Referenced by TBuilderCosmic::buildStereo(), TBuilder0::buildStereo(), TBuilder::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), TBuilder::buildStereoNew(), THistogram::clusters(), THistogram::clusters0(), TCurlFinder::dividing2DTrack(), TrkReco::execute(), TConformalFinder::link(), TCurlFinder::make2DTrack(), TTrackManager::merge(), TTrackManager::salvageAssociateHits(), TSegment0::splitAV(), TSegment::splitAV(), TSegment0::splitParallel(), TSegment::splitParallel(), and TConformalFinder0::standardFinding().

                             {
    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.

Definition at line 101 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, TTrackBase::approach(), TTrackBase::distance(), showlog::err, genRecEmupikp::i, and WireHitUsed.

Referenced by TBuilder0::appendClusters(), and TTrackManager::salvage().

                                                                  {
#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.

Definition at line 152 of file TTrackBase.cxx.

                                                               {
    std::cout << "TTrackBase::appendByDistance !!! not implemented" << std::endl;
    list.removeAll();
}

int TTrackBase::approach

(

TMLink &

)

const [virtual, inherited]

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

Reimplemented in TTrack.

Definition at line 95 of file TTrackBase.cxx.

Referenced by TTrackBase::appendByApproach(), and TTrackBase::testByApproach().

                                   {
    std::cout << "TTrackBase::approach !!! not implemented" << std::endl;
    return -1;
}

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

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

Definition at line 202 of file T3DLine.cxx.

References approach_line(), TMLink::positionOnTrack(), TMLink::positionOnWire(), v, w, and TMLink::wire().

                                                         {

  const TMDCWire& w=*l.wire();
  HepPoint3D xw = w.xyPosition();
  HepPoint3D wireBackwardPosition = w.backwardPosition();
  HepVector3D v = w.direction();

  HepPoint3D onWire,onTrack;

  if(approach_line(wireBackwardPosition,v,onWire,onTrack)<0)
    return(-1);

  // onWire,onTrack filled
  
  if(!doSagCorrection){
    l.positionOnWire(onWire);
    l.positionOnTrack(onTrack);
    return(0);                          // no sag correction
  }
  // Sag correction
  //   loop for sag correction
  double onWire_y = onWire.y();
  double onWire_z = onWire.z();

  unsigned nTrial = 1;
  while(nTrial<100){
    w.wirePosition(onWire_z,xw,wireBackwardPosition,(HepVector3D&)v);
    if(approach_line(wireBackwardPosition,v,onWire,onTrack)<0)
      return(-1);
    if(fabs(onWire_y - onWire.y())<0.0001) break;       // |dy|< 1 micron
    onWire_y = onWire.y();
    onWire_z = onWire.z();

    nTrial++;
  }

  l.positionOnWire(onWire);
  l.positionOnTrack(onTrack);
  return(nTrial);
}

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

caluculate closest points between a line and this track

Definition at line 243 of file T3DLine.cxx.

References k(), s, t(), and x0().

Referenced by approach().

                                                                        {
  //  line = [w0] + s * [v]    -> [onLine]
  //  trk  = [x0] + t * [k]    -> [onTrack]
  //  if [v]//[k] then return(-1) error

  const HepVector3D k = this->k();
  const double v_k = v.dot(k);
  const double v_2 = v.mag2();
  const double k_2 = k.mag2();
  const double tk = v_k*v_k - v_2*k_2;
  if(tk==0) return(-1);

  const HepPoint3D x0 = this->x0();
  const HepPoint3D dx = x0 - w0;
  const double t = dx.dot( v_2 * k - v_k * v)/tk;
  const double s = dx.dot( v_k * k - k_2 * v)/tk;

  onLine = w0 + s * v;
  onTrack = x0 + t * k;
  return(0);
}

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

caluculate closest point between a point and this track

Definition at line 266 of file T3DLine.cxx.

References k(), t(), and x0().

                                                                         {
  //  trk  = [x0] + t * [k]    -> [onTrack]
  //  if [v]//[k] then return(-1) error

  const HepVector3D k = this->k();
  const HepPoint3D x0 = this->x0();
  const HepPoint3D dx = p0 - x0;
  const double t = dx.dot(k)/k.mag2();

  onTrack = x0 + t * k;
  return(0);
}

double T3DLine::chi2

(

void

)

const

returns chi2.

Definition at line 131 of file T3DLine.cxx.

References _chi2.

                              {
  return _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.

Definition at line 317 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_updated, and TTrackBase::update().

Referenced by TBuilder::buildStereoNew(), TConformalFinder::expand(), TTrack::movePivot(), TRunge::SetFlightLength(), and TTrack::szPosition().

                                     {
    if (mask)
        std::cout << "TTrackBase::cores !!! mask is not supported" << std::endl;
    if (! _updated) update();
    return _cores;
}

double T3DLine::cosPhi0

(

void

)

const

Definition at line 95 of file T3DLine.cxx.

References _cos_phi0.

                                 {
  return _cos_phi0;
}

double TTrackBase::distance

(

const TMLink &

)

const [virtual, inherited]

returns distance to a position of TMLink in TMLink space.

Reimplemented in TLine0, and TMLine.

Definition at line 89 of file TTrackBase.cxx.

Referenced by TTrackBase::appendByApproach(), TTrack::fit2D(), and TTrackBase::testByApproach().

                                         {
    std::cout << "TTrackBase::distance !!! not implemented" << std::endl;
    return 0.;
}

double T3DLine::dr

(

void

)

const

Track parameters.

Definition at line 79 of file T3DLine.cxx.

References _dr.

Referenced by pivot().

                            {
  return _dr;
}

int TTrackBase::DropWorst

(

)

[virtual, inherited]

Definition at line 212 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, and genRecEmupikp::i.

Referenced by TrkReco::execute().

                      {
  AList<TMLink> bad;
  int jbad=-1;
  double sigma=0;
  unsigned n = _links.length();
  double chi2=0;
  for (unsigned i = 0; i < n; i++){
   // if(g_sigmalyr1)g_sigmalyr1->fill(sqrt(_links[i]->pull()),_links[i]->wire()->layerId());
    chi2+=_links[i]->pull();
    if (_links[i]->pull() > sigma){
      jbad=i; 
      sigma = _links[i]->pull();
    }
  }
        
  chi2=chi2/n;
  //if(sigma>9){bad.append(_links[jbad]);if(g_sigmalyr)g_sigmalyr->fill(sqrt(sigma),_links[jbad]->wire()->layerId());}
//  if(sigma>20.25&&chi2>2.08){bad.append(_links[jbad]);}
  if(sigma>9){bad.append(_links[jbad]);}
  if (bad.length()) {
    _links.remove(bad);
    _fitted = false;
    _updated = false;
  }
  return bad.length();
}

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

dumps debug information.

Reimplemented in TCircle, TLine0, TMLine, TSegment, TSegment0, and TTrack.

Definition at line 58 of file TTrackBase.cxx.

References TTrackBase::_links, TTrackBase::_updated, Dump(), TTrackBase::mc(), RealDBUtil::npos, and TTrackBase::update().

Referenced by TTrack::dump(), TSegment0::dump(), TSegment::dump(), TMLine::dump(), TLine0::dump(), and TCircle::dump().

                                                                   {
    bool mc = (msg.find("mc") != std::string::npos);
    bool pull = (msg.find("pull") != std::string::npos);
    bool flag = (msg.find("flag") != std::string::npos);
    bool detail = (msg.find("detail") != std::string::npos);
    if (detail)
        mc = pull = flag = true;

    if (detail || (msg.find("layer") != std::string::npos)) {
        if (! _updated) update();
    }
    if (detail || (msg.find("hits") != std::string::npos)) {
        Dump(_links, msg, pre);
    }
}

double T3DLine::dz

(

void

)

const

Definition at line 87 of file T3DLine.cxx.

References _dz.

Referenced by pivot().

                            {
  return _dz;
}

const SymMatrix & T3DLine::Ea

(

const SymMatrix &

)

set error matrix

Definition at line 197 of file T3DLine.cxx.

References _Ea.

                                                {
  _Ea = tEa;
  return(_Ea);
}

const SymMatrix & T3DLine::Ea

(

void

)

const

returns error matrix

Definition at line 116 of file T3DLine.cxx.

References _Ea.

                                      {
  return(_Ea);
}

void TTrackBase::falseFit

(

)

[inline, inherited]

false Fit

Definition at line 229 of file TTrackBase.h.

References TTrackBase::_fitted, and TTrackBase::_fittedWithCathode.

                    {
    _fitted = false;
    _fittedWithCathode = false;
}

int TTrackBase::fit

(

void

)

[virtual, inherited]

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

Definition at line 357 of file TTrackBase.cxx.

References TTrackBase::_fitter, and TMFitter::fit().

Referenced by TBuilder0::buildRphi(), TBuilder::buildRphi(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), and TConformalFinder::expand().

                    {
    return _fitter->fit(* this);
}

bool TTrackBase::fitted

(

void

)

const [inline, inherited]

returns true if fitted.

Definition at line 222 of file TTrackBase.h.

References TTrackBase::_fitted.

Referenced by TCosmicFitter::fit(), TBuilder::searchLine(), and TMSelector::select().

                             {
    return _fitted;
}

bool TTrackBase::fittedWithCathode

(

void

)

const [inline, inherited]

returns true if fitted with cathode hits(TEMPORARY).

Definition at line 243 of file TTrackBase.h.

References TTrackBase::_fittedWithCathode.

                                        {
    return _fittedWithCathode;
}

const TMFitter *const TTrackBase::fitter

(

const TMFitter *

)

[inline, inherited]

sets a default fitter.

Definition at line 261 of file TTrackBase.h.

References TTrackBase::_fitted, and TTrackBase::_fitter.

                                     {
    _fitted = false;
    return _fitter = a;
}

const TMFitter *const TTrackBase::fitter

(

void

)

const [inline, inherited]

returns a pointer to a default fitter.

Definition at line 255 of file TTrackBase.h.

References TTrackBase::_fitter.

Referenced by T3DLine(), TCircle::TCircle(), TLine0::TLine0(), TMLine::TMLine(), TRunge::TRunge(), and TTrack::TTrack().

                             {
    return _fitter;
}

Helix T3DLine::helix

(

void

)

const

approximated helix class

Definition at line 120 of file T3DLine.cxx.

References _dr, _dz, _phi0, _pivot, _tanl, and a().

Referenced by TBuilderCosmic::buildStereo().

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

const TTrackHEP *const TTrackBase::hep

(

void

)

const [inherited]

returns TTrackHEP.

Definition at line 407 of file TTrackBase.cxx.

References TTrackBase::_hep, TTrackBase::_links, TTrackBase::_nHeps, genRecEmupikp::i, ganga-rec::j, max, and TTrackBase::mc().

Referenced by TTrackBase::nHeps().

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

HepVector3D T3DLine::k

(

void

)

const

Definition at line 158 of file T3DLine.cxx.

References _cos_phi0, _sin_phi0, and _tanl.

Referenced by approach_line(), and approach_point().

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

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.

Definition at line 297 of file TTrackBase.cxx.

References TTrackBase::_links.

Referenced by TConformalFinder0::appendCluster(), TBuilder::build(), TBuilderCurl::buildStereo(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), TBuilder::buildStereoNew(), TBuilderCurl::check(), TCurlFinder::check2DCircle(), TCurlFinder::check3DTrack(), TCurlFinder::dividing2DTrack(), TCurlFinder::dumpType1(), TCurlFinder::dumpType2(), TrkReco::execute(), TConformalFinder0::findCloseClusters(), TFastFinder::findCloseHits(), TConformalFinder0::findCloseHits(), TConformalFinder0::findClusterLink(), TCosmicFitter::fit(), TCurlFinder::fitWDD(), TCurlFinder::make2DTrack(), TCurlFinder::make3DTrack(), TCurlFinder::makeWithMC(), TCurlFinder::mask3DTrack(), TTrackManager::merge(), TTrack::movePivot(), TCurlFinder::plotCircle(), TCurlFinder::plotTrack(), TMSelector::preSelect(), TCurlFinder::removeStereo(), TConformalFinder::salvage(), TCurlFinder::salvage3DTrack(), TBuilder::searchLine(), TMDCTsf::segments(), TBuilderCurl::setArcZ(), TConformalFinder::slowFinding2D(), TSegment0::splitAV(), TSegment::splitAV(), TSegment::splitTsf(), TConformalFinder::stereoSegments(), TSegment0::superLayerId(), TSegment::superLayerId(), TTrack::szPosition(), TCurlFinder::trace2DTrack(), TCurlFinder::trace3DTrack(), and TSegment::width().

                                     {
    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 [inline, inherited]

returns a pointer to TTrackMC.

Definition at line 249 of file TTrackBase.h.

References TTrackBase::_mc.

Referenced by TTrackBase::dump(), and TTrackBase::hep().

                         {
    return _mc;
}

unsigned TTrackBase::nCores

(

unsigned

mask = 0

)

const [inherited]

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

Definition at line 325 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_updated, and TTrackBase::update().

Referenced by TTrack::dump(), and TConformalFinder::expand().

                                      {
    if (mask)
        std::cout << "TTrackBase::nCores !!! mask is not supported" << std::endl;
    if (! _updated) update();
    return _cores.length();
}

unsigned T3DLine::ndf

(

void

)

const

returns NDF

Definition at line 127 of file T3DLine.cxx.

References _ndf.

                               {
  return _ndf;
}

unsigned TTrackBase::nHeps

(

void

)

const [inherited]

returns # of contributed TTrackHEP tracks.

Definition at line 443 of file TTrackBase.cxx.

References TTrackBase::_nHeps, and TTrackBase::hep().

                            {
    hep();
    return _nHeps;
}

unsigned TTrackBase::nLinks

(

unsigned

mask = 0

)

const [inherited]

returns # of masked TMLinks assigned to this track object.

Definition at line 305 of file TTrackBase.cxx.

References TTrackBase::_links, and genRecEmupikp::i.

Referenced by TBuilderCurl::buildStereo(), TBuilderCurl::buildStereoMC(), TBuilderCurl::check(), TCurlFinder::check2DCircle(), TCurlFinder::check3DTrack(), TSegment0::clusterType(), TSegment::clusterType(), TCurlFinder::dividing2DTrack(), TCurlFinder::dumpType1(), TCurlFinder::dumpType2(), TCurlFinder::make2DTrack(), TCurlFinder::makeWithMC(), TCurlFinder::merge3DTrack(), TCurlFinder::plotCircle(), TCurlFinder::plotTrack(), TMSelector::preSelect(), TCurlFinder::salvage3DTrack(), TBuilderCurl::setArcZ(), TSegment0::splitComplicated(), TSegment::splitComplicated(), TSegment0::splitParallel(), TSegment::splitParallel(), TSegment0::updateType(), and TSegment::updateType().

                                      {
    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 T3DLine::objectType

(

void

)

const [inline, virtual]

returns object type

Reimplemented from TTrackBase.

Definition at line 150 of file T3DLine.h.

References Line3D.

                                             {
  return Line3D;
}

TMLink * TTrackBase::operator[]

(

unsigned

i

)

const [inline, inherited]

Definition at line 237 of file TTrackBase.h.

References TTrackBase::_links.

                                       {
    return _links[i];
}

double T3DLine::phi0

(

void

)

const

Definition at line 83 of file T3DLine.cxx.

References _phi0.

                              {
  return _phi0;
}

const HepPoint3D & T3DLine::pivot

(

const HepPoint3D &

)

set new pivot

Definition at line 163 of file T3DLine.cxx.

References _cos_phi0, _dr, _dz, _pivot, _sin_phi0, _tanl, dr(), and dz().

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

const HepPoint3D & T3DLine::pivot

(

void

)

const

pivot position

Definition at line 103 of file T3DLine.cxx.

References _pivot.

                                          {
  return _pivot;
}

double T3DLine::reducedchi2

(

void

)

const

returns reduced-chi2

Definition at line 135 of file T3DLine.cxx.

References _chi2, and _ndf.

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

void TTrackBase::refine

(

double

maxSigma

)

[virtual, inherited]

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

Definition at line 194 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_updated, Dump(), and TTrackBase::refineMain().

                               {
    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'.

Definition at line 170 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, Dump(), and TTrackBase::refineMain().

Referenced by TBuilder0::appendClusters(), TBuilderCurl::buildStereo(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo(), TBuilder::buildStereo(), TBuilder0::buildStereo0(), and TBuilderCurl::buildStereoMC().

                                                     {
    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 TTrackBase::remove

(

const AList< TMLink > &

)

[inline, inherited]

removes TMLinks.

Definition at line 213 of file TTrackBase.h.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, and TTrackBase::_updated.

                                          {
    _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.

Definition at line 204 of file TTrackBase.h.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, and TTrackBase::_updated.

Referenced by TConformalFinder0::appendClusters2(), TBuilder::buildStereoNew(), TCurlFinder::dividing2DTrack(), TCurlFinder::makeWithMC(), TTrackManager::maskCurl(), TTrackManager::maskNormal(), TCurlFinder::removeStereo(), TBuilder::searchLine(), TConformalFinder0::specialFinding(), TSegment0::splitComplicated(), and TSegment::splitComplicated().

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

void TTrackBase::removeLinks

(

void

)

[virtual, inherited]

Definition at line 189 of file TTrackBase.cxx.

References TTrackBase::_links.

Referenced by TBuilderCosmic::buildStereo(), and TrkReco::execute().

                           {
  _links.removeAll();
}

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

set track parameters,pivot

Definition at line 175 of file T3DLine.cxx.

References _cos_phi0, _dr, _dz, _phi0, _pivot, _sin_phi0, _tanl, cos(), and sin().

                                                                      {

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

double T3DLine::sinPhi0

(

void

)

const

Definition at line 99 of file T3DLine.cxx.

References _sin_phi0.

                                 {
  return _sin_phi0;
}

double T3DLine::tanl

(

void

)

const

Definition at line 91 of file T3DLine.cxx.

References _tanl.

                              {
  return _tanl;
}

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

Definition at line 240 of file TTrackBase.cxx.

References genRecEmupikp::i, and TTrackBase::testByApproach().

                                                                            {
#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.

Definition at line 255 of file TTrackBase.cxx.

References TTrackBase::approach(), TMDCWireHit::dDrift(), TMLink::distance(), TTrackBase::distance(), TMDCWireHit::drift(), TMLink::dump(), showlog::err, TMLink::hit(), and TMLink::pull().

Referenced by TTrackBase::testByApproach().

                                                                     {
#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 TTrackBase::type

(

void

)

const [inline, virtual, inherited]

returns type. Definition is depending on an object class.

Reimplemented in TTrack.

Definition at line 274 of file TTrackBase.h.

                           {
    return 0;
}

void TTrackBase::update

(

void

)

const [inherited]

update cache.

Reimplemented in TSegment, and TSegment0.

Definition at line 75 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_links, TTrackBase::_updated, genRecEmupikp::i, TMDCWireHit::state(), WireHitFittingValid, and WireHitInvalidForFit.

Referenced by TTrackBase::cores(), TTrackBase::dump(), and TTrackBase::nCores().

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

HepPoint3D T3DLine::x

(

double

)

const

returns position on 3D line

Definition at line 143 of file T3DLine.cxx.

References _cos_phi0, _dr, _dz, _pivot, _sin_phi0, and _tanl.

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

HepPoint3D T3DLine::x0

(

void

)

const

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

Definition at line 151 of file T3DLine.cxx.

References _cos_phi0, _dr, _dz, _pivot, and _sin_phi0.

Referenced by approach_line(), and approach_point().

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

---

Friends And Related Function Documentation

friend class T3DLineFitter [friend]

Definition at line 147 of file T3DLine.h.

---

Member Data Documentation

double T3DLine::_chi2 [private]

Definition at line 138 of file T3DLine.h.

Referenced by chi2(), and reducedchi2().

double T3DLine::_cos_phi0 [private]

Definition at line 142 of file T3DLine.h.

Referenced by a(), cosPhi0(), k(), pivot(), set(), T3DLine(), x(), and x0().

double T3DLine::_dr [private]

Definition at line 131 of file T3DLine.h.

Referenced by a(), dr(), helix(), pivot(), set(), x(), and x0().

double T3DLine::_dz [private]

Definition at line 133 of file T3DLine.h.

Referenced by a(), dz(), helix(), pivot(), set(), x(), and x0().

SymMatrix T3DLine::_Ea [private]

Definition at line 136 of file T3DLine.h.

Referenced by Ea().

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

Definition at line 162 of file TTrackBase.h.

Referenced by TMLine::a(), TLine0::a(), TTrackBase::append(), TTrackBase::appendByApproach(), TMLine::appendByszdistance(), TLine0::appendByszdistance(), TMLine::appendSLY(), TLine0::appendSLY(), TMLine::b(), TLine0::b(), TCircle::center(), TCircle::charge(), TTrack::chi2(), TMLine::chi2(), TLine0::chi2(), TSegment0::direction(), TSegment::direction(), TMLine::distance(), TLine0::distance(), TTrackBase::DropWorst(), TSegment0::dump(), TSegment::dump(), TCircle::dump(), TTrackBase::falseFit(), THelixFitter::fit(), TCosmicFitter::fit(), TMLine::fit2(), TLine0::fit2(), TMLine::fit2p(), TLine0::fit2p(), TMLine::fit2s(), TLine0::fit2s(), TMLine::fit2sp(), TLine0::fit2sp(), TCircle::fitForCurl(), TTrackBase::fitted(), TTrackBase::fitter(), TCircle::impact(), TSegment0::innerMostLayer(), TSegment::innerMostLayer(), TSegment::inners(), TSegment0::innerWidth(), TSegment::innerWidth(), TTrack::ndf(), TSegment0::outerMostLayer(), TSegment::outerMostLayer(), TSegment::outers(), TSegment0::outerWidth(), TSegment::outerWidth(), TSegment0::position(), TSegment::position(), TCircle::pt(), TCircle::radius(), TMLine::reducedChi2(), TLine0::reducedChi2(), TTrackBase::refine(), TMLine::refine(), TLine0::refine(), TTrack::refine2D(), TTrackBase::remove(), TMLine::removeChits(), TLine0::removeChits(), TMLine::removeSLY(), TLine0::removeSLY(), T3DLine(), TRunge::TRunge(), TSegment::TSegment(), TSegment0::TSegment0(), TTrack::TTrack(), TSegment0::update(), TSegment::update(), TSegment0::updateType(), and TSegment::updateType().

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

Definition at line 163 of file TTrackBase.h.

Referenced by TTrackBase::append(), TTrackBase::falseFit(), TTrackBase::fittedWithCathode(), TTrackBase::remove(), T3DLine(), TRunge::TRunge(), and TTrack::TTrack().

const T3DLineFitter T3DLine::_fitter = T3DLineFitter("T3DLine Default fitter") [static, private]

Reimplemented from TTrackBase.

Definition at line 145 of file T3DLine.h.

Referenced by T3DLine().

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

Definition at line 161 of file TTrackBase.h.

Referenced by TTrackBase::append(), TTrackBase::appendByApproach(), TMLine::appendByszdistance(), TLine0::appendByszdistance(), TMLine::appendSLY(), TLine0::appendSLY(), TTrack::assign(), TBuilderCurl::buildStereo(), TMLine::chi2(), TLine0::chi2(), TTrackBase::DropWorst(), TTrackBase::dump(), TTrack::dump(), TSegment0::dump(), TSegment::dump(), TMLine::dump(), TLine0::dump(), TCircle::dump(), TSegment::expandSeg(), TMLine::fit2(), TLine0::fit2(), TTrack::fit2D(), TMLine::fit2p(), TLine0::fit2p(), TMLine::fit2s(), TLine0::fit2s(), TMLine::fit2sp(), TLine0::fit2sp(), TCircle::fitForCurl(), TTrackBase::hep(), TTrackBase::links(), TTrack::movePivot(), TTrackBase::nLinks(), TTrackBase::operator[](), TSegment0::rangeX(), TSegment::rangeX(), TMLine::reducedChi2(), TLine0::reducedChi2(), TTrackBase::refine(), TMLine::refine(), TLine0::refine(), TTrack::refine2D(), TTrackBase::refineMain(), TTrackBase::remove(), TMLine::removeChits(), TLine0::removeChits(), TTrackBase::removeLinks(), TMLine::removeSLY(), TLine0::removeSLY(), TSegment::segSalvage(), TSegment0::solveDualHits(), TSegment::solveDualHits(), TSegment::solveLR(), TSegment::solveThreeHits(), TSegment0::splitAV(), TSegment::splitAV(), TSegment0::splitComplicated(), TSegment::splitComplicated(), TSegment0::splitDual(), TSegment::splitDual(), TSegment0::splitParallel(), TSegment::splitParallel(), TSegment::splitTsf(), TSegment::TSegment(), TSegment0::TSegment0(), TTrack::TTrack(), TTrackBase::update(), TSegment0::update(), TSegment::update(), TSegment0::updateDuality(), TSegment::updateDuality(), TSegment0::updateType(), and TSegment::updateType().

TTrackMC* TTrackBase::_mc [protected, inherited]

Definition at line 165 of file TTrackBase.h.

Referenced by TTrackBase::mc().

unsigned T3DLine::_ndf [private]

Definition at line 139 of file T3DLine.h.

Referenced by ndf(), and reducedchi2().

double T3DLine::_phi0 [private]

Definition at line 132 of file T3DLine.h.

Referenced by a(), helix(), phi0(), set(), and T3DLine().

HepPoint3D T3DLine::_pivot [private]

Definition at line 129 of file T3DLine.h.

Referenced by helix(), pivot(), set(), x(), and x0().

double T3DLine::_sin_phi0 [private]

Definition at line 143 of file T3DLine.h.

Referenced by a(), k(), pivot(), set(), sinPhi0(), T3DLine(), x(), and x0().

double T3DLine::_tanl [private]

Definition at line 134 of file T3DLine.h.

Referenced by a(), helix(), k(), pivot(), set(), tanl(), and x().

---