TLine0 Class Reference

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

#include <TLine0.h>

Inheritance diagram for TLine0:

Public Member Functions
	TLine0 ()
	Constructor.
	TLine0 (const AList< TMLink > &)
	Constructor.
virtual	~TLine0 ()
	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.
double	a (void) const
	returns coefficient a.
double	b (void) const
	returns coefficient b.
double	chi2 (void) const
	returns chi2.
double	reducedChi2 (void) const
	returns reduced-chi2.
double	distance (const TMLink &) const
	returns distance to a position of TMLink itself. (not to a wire)
int	fit2 ()
	fits itself. Error was happened if return value is not zero.
int	fit2s ()
	fits itself using single hits in a wire-layer. Error was happened if return value is not zero.
int	fit2p ()
	fits itself using isolated hits. Error was happened if return value is not zero.
int	fit2sp ()
	fits itself using single hits in a wire-layer. Error was happened if return value is not zero.
void	removeChits ()
	remove extremly bad points.
void	refine (AList< TMLink > &list, float maxSigma)
	remove bad points by chi2. Bad points are returned in a 'list'. fit() should be called before calling this function.
void	removeSLY (AList< TMLink > &list)
void	appendSLY (AList< TMLink > &list)
void	appendByszdistance (AList< TMLink > &list, unsigned isl, float maxSigma)
void	property (double a, double b, double det)
	sets circle properties.
	TLine0 ()
	Constructor.
	TLine0 (const AList< TMLink > &)
	Constructor.
virtual	~TLine0 ()
	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.
double	a (void) const
	returns coefficient a.
double	b (void) const
	returns coefficient b.
double	chi2 (void) const
	returns chi2.
double	reducedChi2 (void) const
	returns reduced-chi2.
double	distance (const TMLink &) const
	returns distance to a position of TMLink itself. (not to a wire)
int	fit2 ()
	fits itself. Error was happened if return value is not zero.
int	fit2s ()
	fits itself using single hits in a wire-layer. Error was happened if return value is not zero.
int	fit2p ()
	fits itself using isolated hits. Error was happened if return value is not zero.
int	fit2sp ()
	fits itself using single hits in a wire-layer. Error was happened if return value is not zero.
void	removeChits ()
	remove extremly bad points.
void	refine (AList< TMLink > &list, float maxSigma)
	remove bad points by chi2. Bad points are returned in a 'list'. fit() should be called before calling this function.
void	removeSLY (AList< TMLink > &list)
void	appendSLY (AList< TMLink > &list)
void	appendByszdistance (AList< TMLink > &list, unsigned isl, float maxSigma)
void	property (double a, double b, double det)
	sets circle properties.
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 represent a track in tracking.

Definition at line 30 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

Constructor & Destructor Documentation

TLine0() [1/4]

TLine0::TLine0

(

)

Constructor.

Definition at line 23 of file TLine0.cxx.

: TTrackBase(),
  _a(0.),
  _b(0.),
  _det(0.),
  _fittedUpdated(false),
  _chi2(0.),
  _reducedChi2(0.) {
 
    //...Set a defualt fitter...
    fitter(& TLine0::_fitter);
}

TLine0() [2/4]

TLine0::TLine0

(

const AList< TMLink > &

a

)

Constructor.

Definition at line 36 of file TLine0.cxx.

: TTrackBase(a),
  _a(0.),
  _b(0.),
  _det(0.),
  _fittedUpdated(false),
  _chi2(0.),
  _reducedChi2(0.) {
 
    //...Set a defualt fitter...
    fitter(& TLine0::_fitter);
}

~TLine0() [1/2]

TLine0::~TLine0 ( )

virtual

Destructor.

Definition at line 49 of file TLine0.cxx.

                {
}

TLine0() [3/4]

TLine0::TLine0

(

)

Constructor.

TLine0() [4/4]

TLine0::TLine0

(

const AList< TMLink > &

)

Constructor.

~TLine0() [2/2]

virtual TLine0::~TLine0 ( )

virtual

Destructor.

Member Function Documentation

a() [1/2]

double TLine0::a ( void  ) const

inline

returns coefficient a.

Definition at line 127 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

                    {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::a !!! fit not performed" << std::endl;
#endif
    return _a;
}

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

a() [2/2]

double TLine0::a

(

void

)

const

returns coefficient a.

appendByszdistance() [1/2]

void TLine0::appendByszdistance	(	AList< TMLink > &	list,
		unsigned	isl,
		float	maxSigma
	)

Definition at line 565 of file TLine0.cxx.

                                                                             {
 
    //... intialize
    unsigned nb = _links.length();
    
    //....Select good hit
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * list[i];
        if(l.hit()->wire()->superLayerId() == isl){
      double dist = distance(l);
      if (dist < maxSigma) {
            _links.append(l);
      }
    }
    }
 
    unsigned na = _links.length();
    if(nb != na){
      AList<TMLink> bad;
      //... remove duplicated hits
      for(unsigned i = 0; i<na ;i++){
         TMLink & l = * _links[i];
         if(i < na - 1) {
           TMLink & lnext = * _links[i+1];
           if(l.hit()->wire()->layerId() == lnext.hit()->wire()->layerId() ){
             if(l.hit()->wire()->localId() == lnext.hit()->wire()->localId()){
               bad.append(l);
         }
       }
     }
      }
      if(bad.length() > 0) _links.remove(bad);
      _fitted = false;
      _fittedUpdated = false;
    }
}

Referenced by TBuilder0::buildStereo().

appendByszdistance() [2/2]

void TLine0::appendByszdistance	(	AList< TMLink > &	list,
		unsigned	isl,
		float	maxSigma
	)

appendSLY() [1/2]

void TLine0::appendSLY

(

AList< TMLink > &

list

)

Definition at line 558 of file TLine0.cxx.

                                      {
      _links.append(list);
      _fitted = false;
      _fittedUpdated = false;
}

appendSLY() [2/2]

void TLine0::appendSLY

(

AList< TMLink > &

list

)

b() [1/2]

double TLine0::b ( void  ) const

inline

returns coefficient b.

Definition at line 136 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

                    {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::b !!! fit not performed" << std::endl;
#endif
    return _b;
}

Referenced by TBuilder0::buildStereo(), TBuilderCosmic::buildStereo(), TBuilderCurl::buildStereoMC(), and property().

b() [2/2]

double TLine0::b

(

void

)

const

returns coefficient b.

chi2() [1/2]

double TLine0::chi2

(

void

)

const

returns chi2.

Definition at line 113 of file TLine0.cxx.

                       {
#ifdef TRKRECO_DEBUG_DETAIL
    if (! _fitted) std::cout << "TLine0::chi2 !!! fit not performed" << std::endl;
#endif
 
    if (_fittedUpdated) return _chi2;
    _chi2 = 0.;
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
    double c = y - _a * x - _b;
    _chi2 += c * c;
    }
    _fittedUpdated = true;
    return _chi2;
}

Referenced by TBuilder0::buildStereo(), and reducedChi2().

chi2() [2/2]

double TLine0::chi2

(

void

)

const

returns chi2.

distance() [1/2]

double TLine0::distance ( const TMLink &  l ) const

inlinevirtual

returns distance to a position of TMLink itself. (not to a wire)

Reimplemented from TTrackBase.

Definition at line 145 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

                                       {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TLine0::distance !!! fit not performed" << std::endl;
#endif
    double dy = fabs(_a * l.position().x() + _b - l.position().y());
    double invCos = sqrt(1. + _a * _a);
    return dy / invCos;
}

Referenced by appendByszdistance(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), and refine().

distance() [2/2]

double TLine0::distance ( const TMLink &  ) const

virtual

returns distance to a position of TMLink itself. (not to a wire)

Reimplemented from TTrackBase.

dump() [1/2]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

Definition at line 53 of file TLine0.cxx.

                                                               {
    bool def = false;
    if (msg == "") def = true;
 
    if (def || msg.find("line") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "#links=" << _links.length();
    std::cout << ",a=" << _a;
    std::cout << ",b=" << _b;
    std::cout << ",det=" << _det;
    std::cout << std::endl;
    }
    if (! def) TTrackBase::dump(msg, pre);
}

dump() [2/2]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

fit2() [1/2]

int TLine0::fit2

(

)

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

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

Definition at line 173 of file TLine0.cxx.

             {
  //    if (_fitted) return 0;
 
    unsigned n = _links.length();
    int mask[100] = {0};
    int nsl[11] = {64,80,96,128,144,160,176,208,240,256,288};
    int npos = 0, nneg = 0;
    for (unsigned i = 0; i < n -1 ; i++) {
    TMLink & l = * _links[i];
        for (unsigned j = i+1; j < n  ; j++) {
       TMLink & s = * _links[j];
           if(l.hit()->wire()->layerId() == s.hit()->wire()->layerId()){
         //... Check 3 consective hits
             if(i > 0 && (mask[i-1] == 1 && mask[j] == 1) ){
               TMLink & t = * _links[i-1];
               if(l.hit()->wire()->layerId() == t.hit()->wire()->layerId()){
             mask[i] = 1;
           }
         }
             int ilast = nsl[l.hit()->wire()->superLayerId()]-1;
             int ilocal = l.hit()->wire()->localId();
             int jlocal = s.hit()->wire()->localId();
             if(ilocal > 0 && ilocal < ilast){
               if(abs(jlocal-ilocal) > 1 ) {
                  mask[i] = 1;
                  mask[j] = 1;
           }
         }else if(ilocal == 0){
               if(jlocal > 1 && jlocal < ilast){
                 mask[i] = 1;
                 mask[j] = 1;
           }
         }else if(ilocal == ilast){
               if(jlocal > 0 && jlocal < ilast-1){
                 mask[i] = 1;
                 mask[j] = 1;
           }
         }
           }
        }
    //...
        if(mask[i] == 0){
          if(l.position().y() >= 0) npos += 1;
          if(l.position().y() <  0) nneg += 1;
    }
    }            
 
    //....
    double sumX = 0., sumY = 0., sumX2 = 0., sumXY = 0., sumY2 = 0.;
    int nused = 0;
    int lyid[2];
    for (unsigned i = 0; i < n; i++) {
 
        if(mask[i] == 1) continue;
 
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
        if(abs(npos-nneg) > 3){
          if(npos > nneg && y < 0 ) continue;
          if(npos < nneg && y > 0 ) continue;
    }
    sumX  += x;
    sumY  += y;
    sumX2 += x * x;
    sumXY += x * y;
    sumY2 += y * y;
        if(nused < 2){
          lyid[nused] = l.hit()->wire()->layerId();
    }
        nused += 1;
    }
 
     if(nused < 2 || (nused == 2 && lyid[0]==lyid[1])) {
       return -2;
     }
     double sum = double(nused);
    _det = sum * sumX2 - sumX * sumX;
    if(_det == 0.) {
    return -1;
    }
    _a = (sumXY * sum - sumX * sumY) / _det;
    _b = (sumX2 * sumY - sumX * sumXY) / _det;    
 
    _fitted = true;
    return 0;
}

fit2() [2/2]

int TLine0::fit2

(

)

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

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

fit2p() [1/2]

int TLine0::fit2p

(

)

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

Definition at line 401 of file TLine0.cxx.

              {
  //    if (_fitted) return 0;
 
    unsigned n = _links.length();
    int mask[100] = {0};
    int nsl[11] = {64,80,96,128,144,160,176,208,240,256,288};
    double phi_ave = 0.;
    int nphi = 0;
    double Crad = 180./3.141592;
    for (unsigned i = 0; i < n -1 ; i++) {
    TMLink & l = * _links[i];
        for (unsigned j = i+1; j < n  ; j++) {
       TMLink & s = * _links[j];
           if(l.hit()->wire()->layerId() == s.hit()->wire()->layerId()){
         //... Check 3 consective hits
             if(i > 0 && (mask[i-1] == 1 && mask[j] == 1) ){
               TMLink & t = * _links[i-1];
               if(l.hit()->wire()->layerId() == t.hit()->wire()->layerId()){
             mask[i] = 1;
           }
         }
             int ilast = nsl[l.hit()->wire()->superLayerId()]-1;
             int ilocal = l.hit()->wire()->localId();
             int jlocal = s.hit()->wire()->localId();
             if(ilocal > 0 && ilocal < ilast){
               if(abs(jlocal-ilocal) > 1 ) {
                  mask[i] = 1;
                  mask[j] = 1;
           }
         }else if(ilocal == 0){
               if(jlocal > 1 && jlocal < ilast){
                 mask[i] = 1;
                 mask[j] = 1;
           }
         }else if(ilocal == ilast){
               if(jlocal > 0 && jlocal < ilast-1){
                 mask[i] = 1;
                 mask[j] = 1;
           }
         }
           }
        }
    //...
    //...
        if(mask[i] != 1){
          double phi = Crad*atan2(l.position().y(), l.position().x());
          phi_ave += phi;
          nphi += 1;
    }
    }            
 
    //...
    if(mask[n-1] != 1){
      TMLink & l = * _links[n-1];
      double phi = Crad*atan2(l.position().y(), l.position().x());
      phi_ave += phi;
      nphi += 1;
    }
    double phi_max = 0.;
    double phi_min = 0.;
    if(nphi> 0){
      phi_max = phi_ave/n + 40;
      phi_min = phi_ave/n - 40;
    }
 
    //....
    double sumX = 0., sumY = 0., sumX2 = 0., sumXY = 0., sumY2 = 0.;
    int nused = 0;
    int lyid[2];
    for (unsigned i = 0; i < n; i++) {
 
        if(mask[i] == 1) continue;
 
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
        double phi = Crad*atan2(l.position().y(), l.position().x());
        if(phi > phi_max && phi<phi_min ) continue;
 
    sumX  += x;
    sumY  += y;
    sumX2 += x * x;
    sumXY += x * y;
    sumY2 += y * y;
        if(nused < 2){
          lyid[nused] = l.hit()->wire()->layerId();
    }
        nused += 1;
    }
 
     if(nused < 2 || (nused == 2 && lyid[0]==lyid[1])) {
       return -2;
     }
     double sum = double(nused);
    _det = sum * sumX2 - sumX * sumX;
    if(_det == 0.) {
    return -1;
    }
    _a = (sumXY * sum - sumX * sumY) / _det;
    _b = (sumX2 * sumY - sumX * sumXY) / _det;    
 
    _fitted = true;
    return 0;
}

Referenced by TBuilder0::buildStereo0().

fit2p() [2/2]

int TLine0::fit2p

(

)

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

fit2s() [1/2]

int TLine0::fit2s

(

)

fits itself using single hits in a wire-layer. Error was happened if return value is not zero.

Definition at line 263 of file TLine0.cxx.

              {
  //    if (_fitted) return 0;
 
    unsigned n = _links.length();
    int mask[100] = {0};
    int npos = 0, nneg = 0;
    for (unsigned i = 0; i < n -1 ; i++) {
    TMLink & l = * _links[i];
        for (unsigned j = i+1; j < n  ; j++) {
       TMLink & s = * _links[j];
           if(l.hit()->wire()->layerId() == s.hit()->wire()->layerId()){
              mask[i] = 1;
              mask[j] = 1;
           }
        }
    //...
        if(mask[i] == 0){
          if(l.position().y() >= 0) npos += 1;
          if(l.position().y() <  0) nneg += 1;
    }
    }            
 
    //....
    double sumX = 0., sumY = 0., sumX2 = 0., sumXY = 0., sumY2 = 0.;
    int nused = 0;
    int lyid[2];
    for (unsigned i = 0; i < n; i++) {
 
        if(mask[i] == 1) continue;
 
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
        if(npos > nneg && y < 0 ) continue;
        if(npos < nneg && y > 0 ) continue;
 
    sumX  += x;
    sumY  += y;
    sumX2 += x * x;
    sumXY += x * y;
    sumY2 += y * y;
        nused += 1;
    }
 
     if(nused < 4) {
       return -2;
     }
     double sum = double(nused);
    _det = sum * sumX2 - sumX * sumX;
    if(_det == 0.) {
    return -1;
    }
    _a = (sumXY * sum - sumX * sumY) / _det;
    _b = (sumX2 * sumY - sumX * sumXY) / _det;    
 
    _fitted = true;
    return 0;
}

fit2s() [2/2]

int TLine0::fit2s

(

)

fits itself using single hits in a wire-layer. Error was happened if return value is not zero.

fit2sp() [1/2]

int TLine0::fit2sp

(

)

fits itself using single hits in a wire-layer. Error was happened if return value is not zero.

Definition at line 323 of file TLine0.cxx.

               {
  //    if (_fitted) return 0;
 
    unsigned n = _links.length();
    int mask[100] = {0};
    double phi_ave = 0.;
    int nphi = 0;
    double Crad = 180./3.141592;
    for (unsigned i = 0; i < n -1 ; i++) {
    TMLink & l = * _links[i];
        for (unsigned j = i+1; j < n  ; j++) {
       TMLink & s = * _links[j];
           if(l.hit()->wire()->layerId() == s.hit()->wire()->layerId()){
              mask[i] = 1;
              mask[j] = 1;
           }
        }
    //...
        if(mask[i] != 1){
          double phi = Crad*atan2(l.position().y(), l.position().x());
          phi_ave += phi;
          nphi += 1;
    }
    }            
 
    //...
    if(mask[n-1] != 1){
      TMLink & l = * _links[n-1];
      double phi = Crad*atan2(l.position().y(), l.position().x());
      phi_ave += phi;
      nphi += 1;
    }
    double phi_max = 0.;
    double phi_min = 0.;
    if(nphi> 0){
      phi_max = phi_ave/n + 40;
      phi_min = phi_ave/n - 40;
    }
 
    //....
    double sumX = 0., sumY = 0., sumX2 = 0., sumXY = 0., sumY2 = 0.;
    int nused = 0;
    int lyid[2];
    for (unsigned i = 0; i < n; i++) {
 
        if(mask[i] == 1) continue;
 
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
        double phi = Crad*atan2(l.position().y(), l.position().x());
        if(phi > phi_max && phi<phi_min ) continue;
 
    sumX  += x;
    sumY  += y;
    sumX2 += x * x;
    sumXY += x * y;
    sumY2 += y * y;
        nused += 1;
    }
 
     if(nused < 4) {
       return -2;
     }
     double sum = double(nused);
    _det = sum * sumX2 - sumX * sumX;
    if(_det == 0.) {
    return -1;
    }
    _a = (sumXY * sum - sumX * sumY) / _det;
    _b = (sumX2 * sumY - sumX * sumXY) / _det;    
 
    _fitted = true;
    return 0;
}

Referenced by TBuilder0::buildStereo0().

fit2sp() [2/2]

int TLine0::fit2sp

(

)

fits itself using single hits in a wire-layer. Error was happened if return value is not zero.

objectType() [1/2]

unsigned TLine0::objectType ( void  ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 164 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

                             {
    return Line;
}

objectType() [2/2]

virtual unsigned TLine0::objectType ( void  ) const

virtual

returns type.

Reimplemented from TTrackBase.

property() [1/2]

void TLine0::property ( double  a, double  b, double  det  )

inline

sets circle properties.

Definition at line 156 of file InstallArea/include/TrkReco/TrkReco/TLine0.h.

                                               {
    _a = a;
    _b = b;
    _det = det;
}

property() [2/2]

void TLine0::property	(	double	a,
		double	b,
		double	det
	)

sets circle properties.

reducedChi2() [1/2]

double TLine0::reducedChi2

(

void

)

const

returns reduced-chi2.

Definition at line 604 of file TLine0.cxx.

                              {
#ifdef TRKRECO_DEBUG_DETAIL
    if (! _fitted) std::cout << "TLine0::reducedChi2 !!! fit not performed" << std::endl;
#endif
 
    if (_fittedUpdated) return _reducedChi2;
    double chi2 = 0.;
    double scale = 20.;
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * _links[i];
 
    double x = l.position().x();
    double y = l.position().y();
    double c = y - _a * x - _b;
        double err = scale * l.hit()->dDrift(); 
    chi2 += c * c / err / err;
    }
 
    _reducedChi2 = chi2/(n-2);
    _fittedUpdated = true;
    return _reducedChi2;
}

Referenced by TBuilder0::buildStereo().

reducedChi2() [2/2]

double TLine0::reducedChi2

(

void

)

const

returns reduced-chi2.

refine() [1/2]

void TLine0::refine	(	AList< TMLink > &	list,
		float	maxSigma
	)

remove bad points by chi2. Bad points are returned in a 'list'. fit() should be called before calling this function.

Definition at line 134 of file TLine0.cxx.

                                                   {
    AList<TMLink> bad;
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++) {
    TMLink & l = * _links[i];
    double dist = distance(l);
    if (dist > maxSigma) bad.append(l);
    }
 
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    TLine0::refine ... rejected hits:max distance=" << maxSigma;
    std::cout << std::endl;
    bad.sort(SortByWireId);
    for (unsigned i = 0; i < bad.length(); i++) {
    TMLink & l = * _links[i];
    std::cout << "        ";
    std::cout << l.wire()->layerId() << "-";
    std::cout << l.wire()->localId();
    std::cout << "(";
    if (l.hit()->mc()) {
        if (l.hit()->mc()->hep()) std::cout << l.hit()->mc()->hep()->id();
        else std::cout << "0";
    }
    std::cout << "),";
    std::cout << l.position() << "," << distance(l);
    if (distance(l) > maxSigma) std::cout << " X";
    std::cout << std::endl;
    }
#endif
 
    if (bad.length()) {
    _links.remove(bad);
    list.append(bad);
    _fitted = false;
    _fittedUpdated = false;
    }
}

Referenced by TBuilder0::buildStereo(), and TBuilder0::buildStereo0().

refine() [2/2]

void TLine0::refine	(	AList< TMLink > &	list,
		float	maxSigma
	)

remove bad points by chi2. Bad points are returned in a 'list'. fit() should be called before calling this function.

removeChits() [1/2]

void TLine0::removeChits

(

)

remove extremly bad points.

Definition at line 508 of file TLine0.cxx.

                    {
 
    unsigned n = _links.length();
    int nlyr[50] = {0};
    int nneg = 0, npos = 0;
    for (unsigned i = 0; i < n -1 ; i++) {
    TMLink & l = * _links[i];
        nlyr[l.hit()->wire()->layerId()] += 1;
        if(l.position().y() < 0.){
           nneg += 1;
    }else {
           npos += 1;
    }
    }            
 
    //...
    AList<TMLink> bad;
    for (unsigned i = 0; i < n; i++) {
 
    TMLink & l = * _links[i];
 
    //...if # of hits in a wire layer, don't use...
    if (nlyr[l.hit()->wire()->layerId()] > 3) {
            bad.append(l);
            continue;
    }
    //...remove extremely bad poinits
        if(abs(nneg - npos) > 3) {
          if(npos > nneg && l.position().y() < 0 ) bad.append(l);
          if(npos < nneg && l.position().y() > 0 ) bad.append(l);
    }
    }
    //...
    if (bad.length() > 0 && bad.length() < n) {
    _links.remove(bad);
    }
 
    //... For the next fit
   _fitted = false;
   _fittedUpdated = false;
}

Referenced by TBuilder0::buildStereo0().

removeChits() [2/2]

void TLine0::removeChits

(

)

remove extremly bad points.

removeSLY() [1/2]

void TLine0::removeSLY

(

AList< TMLink > &

list

)

Definition at line 551 of file TLine0.cxx.

                                      {
      _links.remove(list);
      _fitted = false;
      _fittedUpdated = false;
}

Referenced by TBuilder0::buildStereo().

removeSLY() [2/2]

void TLine0::removeSLY

(

AList< TMLink > &

list

)

---

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

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