TMLine Class Reference

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

#include <TMLine.h>

Inheritance diagram for TMLine:

Public Member Functions
	TMLine ()
	Constructor.
	TMLine (const AList< TMLink > &)
	Constructor.
virtual	~TMLine ()
	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)
double	distance (const HepPoint3D &) const
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.
	TMLine ()
	Constructor.
	TMLine (const AList< TMLink > &)
	Constructor.
virtual	~TMLine ()
	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)
double	distance (const HepPoint3D &) const
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 40 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

Constructor & Destructor Documentation

TMLine() [1/4]

TMLine::TMLine

(

)

Constructor.

Definition at line 23 of file TMLine.cxx.

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

TMLine() [2/4]

TMLine::TMLine

(

const AList< TMLink > &

a

)

Constructor.

Definition at line 36 of file TMLine.cxx.

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

~TMLine() [1/2]

TMLine::~TMLine ( )

virtual

Destructor.

Definition at line 49 of file TMLine.cxx.

                {
}

TMLine() [3/4]

TMLine::TMLine

(

)

Constructor.

TMLine() [4/4]

TMLine::TMLine

(

const AList< TMLink > &

)

Constructor.

~TMLine() [2/2]

virtual TMLine::~TMLine ( )

virtual

Destructor.

Member Function Documentation

a() [1/2]

double TMLine::a ( void  ) const

inline

returns coefficient a.

Definition at line 147 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

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

Referenced by TBuilder::build(), TBuilder::buildStereo(), TBuilderCurl::buildStereo(), and property().

a() [2/2]

double TMLine::a

(

void

)

const

returns coefficient a.

appendByszdistance() [1/2]

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

Definition at line 568 of file TMLine.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;
    }
}

appendByszdistance() [2/2]

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

appendSLY() [1/2]

void TMLine::appendSLY

(

AList< TMLink > &

list

)

Definition at line 561 of file TMLine.cxx.

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

appendSLY() [2/2]

void TMLine::appendSLY

(

AList< TMLink > &

list

)

b() [1/2]

double TMLine::b ( void  ) const

inline

returns coefficient b.

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

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

Referenced by TBuilder::build(), TBuilder::buildStereo(), TBuilderCurl::buildStereo(), and property().

b() [2/2]

double TMLine::b

(

void

)

const

returns coefficient b.

chi2() [1/2]

double TMLine::chi2

(

void

)

const

returns chi2.

Definition at line 113 of file TMLine.cxx.

                       {
#ifdef TRKRECO_DEBUG
    if (! _fitted)
    std::cout << "TMLine::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 reducedChi2().

chi2() [2/2]

double TMLine::chi2

(

void

)

const

returns chi2.

distance() [1/4]

double TMLine::distance ( const HepPoint3D &  p ) const

inline

Definition at line 176 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

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

distance() [2/4]

double TMLine::distance

(

const HepPoint3D &

)

const

distance() [3/4]

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

inlinevirtual

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

Reimplemented from TTrackBase.

Definition at line 165 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

                                       {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TMLine::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(), TBuilder::buildStereo(), refine(), TBuilder::removeFarSegment(), TBuilder::searchLine(), and TBuilder::selectStereoSegment().

distance() [4/4]

double TMLine::distance ( const TMLink &  ) const

virtual

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

Reimplemented from TTrackBase.

dump() [1/2]

void TMLine::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 TMLine.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 TMLine::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 TMLine::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 174 of file TMLine.cxx.

             {
  std::cout<<"   "<<__FILE__<<"   "<<__LINE__<<" ERROR : nsl"<<std::endl;
  //    if (_fitted) return 0;
 
    unsigned n = _links.length();
    int mask[100] = {0};
    int nsl[11] = {64,80,96,128,144,160,192,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 TMLine::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 TMLine::fit2p

(

)

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

Definition at line 403 of file TMLine.cxx.

              {
  //    if (_fitted) return 0;
 
  std::cout<<"   "<<__FILE__<<"   "<<__LINE__<<" ERROR : nsl"<<std::endl;
    unsigned n = _links.length();
    int mask[100] = {0};
    int nsl[11] = {64,80,96,128,144,160,192,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;
}

fit2p() [2/2]

int TMLine::fit2p

(

)

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

fit2s() [1/2]

int TMLine::fit2s

(

)

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

Definition at line 265 of file TMLine.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 TMLine::fit2s

(

)

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

fit2sp() [1/2]

int TMLine::fit2sp

(

)

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

Definition at line 325 of file TMLine.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;
}

fit2sp() [2/2]

int TMLine::fit2sp

(

)

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

objectType() [1/2]

unsigned TMLine::objectType ( void  ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 195 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

                             {
    return Line;
}

objectType() [2/2]

virtual unsigned TMLine::objectType ( void  ) const

virtual

returns type.

Reimplemented from TTrackBase.

property() [1/2]

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

inline

sets circle properties.

Definition at line 187 of file InstallArea/include/TrkReco/TrkReco/TMLine.h.

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

property() [2/2]

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

sets circle properties.

reducedChi2() [1/2]

double TMLine::reducedChi2

(

void

)

const

returns reduced-chi2.

Definition at line 607 of file TMLine.cxx.

                              {
#ifdef TRKRECO_DEBUG
    if (! _fitted)
    std::cout << "TMLine::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 = 1.;
    if (l.hit()) err = scale * l.hit()->dDrift();
    chi2 += c * c / err / err;
    }
 
    _reducedChi2 = chi2/(n-2);
    _fittedUpdated = true;
    return _reducedChi2;
}

reducedChi2() [2/2]

double TMLine::reducedChi2

(

void

)

const

returns reduced-chi2.

refine() [1/2]

void TMLine::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 135 of file TMLine.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 << "    TMLine::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;
    }
}

refine() [2/2]

void TMLine::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 TMLine::removeChits

(

)

remove extremly bad points.

Definition at line 511 of file TMLine.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;
}

removeChits() [2/2]

void TMLine::removeChits

(

)

remove extremly bad points.

removeSLY() [1/2]

void TMLine::removeSLY

(

AList< TMLink > &

list

)

Definition at line 554 of file TMLine.cxx.

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

removeSLY() [2/2]

void TMLine::removeSLY

(

AList< TMLink > &

list

)

---

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

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