TCircle Class Reference

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

#include <TCircle.h>

Inheritance diagram for TCircle:

Public Member Functions
	TCircle (const AList< TMLink > &links)
	Constructor.
virtual	~TCircle ()
	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.
const HepPoint3D &	center (void) const
	returns position of center.
double	radius (void) const
	returns radius.
double	pt (void) const
	returns Pt.
double	impact (void) const
	returns impact parameter to the origin.
double	charge (void) const
	returns charge.
double	weight (const TMLink &l) const
	returns weight of TMLink in order to fit and make a circle.
int	fitForCurl (int ipConst=0)
	fits itself. Error was happened if return value is not zero.
void	property (double charge, double radius, HepPoint3D center)
	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).

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

Detailed Description

A class to represent a circle in tracking.

Definition at line 42 of file TCircle.h.

Constructor & Destructor Documentation

TCircle()

TCircle::TCircle

(

const AList< TMLink > &

links

)

Constructor.

Definition at line 22 of file TCircle.cxx.

: TTrackBase(a), _radius(0.), _charge(0.) {
 
    //...Set a defualt fitter...
    fitter(& TCircle::_fitter);
}

~TCircle()

TCircle::~TCircle ( )

virtual

Destructor.

Definition at line 29 of file TCircle.cxx.

                  {
}

Member Function Documentation

center()

const HepPoint3D & TCircle::center ( void  ) const

inline

returns position of center.

Definition at line 108 of file TCircle.h.

                          {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::center !!! fit not performed" << std::endl;
#endif
    return _center;
}

Referenced by TTrack::TTrack().

charge()

double TCircle::charge ( void  ) const

inline

returns charge.

Definition at line 147 of file TCircle.h.

                          {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::charge !!! fit not performed" << std::endl;
#endif
    return _charge;
 
}

dump()

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

Definition at line 33 of file TCircle.cxx.

                                                                {
    bool def = false;
    if (msg == "") def = true;
 
    if (def || msg.find("circle") != std::string::npos || msg.find("detail") != std::string::npos) {
    std::cout << pre;
    std::cout << "#links=" << _links.length();
    if (_fitted) {
        std::cout << ",charge=" << _charge;
        std::cout << ",center=" << _center;
        std::cout << ",radius=" << _radius;
        std::cout << std::endl << pre;
        std::cout << "pt=" << pt();
        std::cout << ",impact=" << impact();
        std::cout << std::endl;
    }
    else {
        std::cout << ", not fitted yet" << std::endl;
    }
    }
    if (! def) TTrackBase::dump(msg, pre);
}

Referenced by TBuilder::buildRphi().

fitForCurl()

int TCircle::fitForCurl

(

int

ipConst = 0

)

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

fits itself only for curl finder. Error was happened if return value is not zero.

Definition at line 136 of file TCircle.cxx.

                               {
    unsigned n = _links.length();
    if (n < 3) return -1;
 
    //IP check
    unsigned flagIP = 1;
    unsigned layerID = _links[0]->hit()->wire()->layerId();
    for(unsigned i = 0; i < n; i++){
      if(layerID != _links[i]->hit()->wire()->layerId()){
        flagIP = 0;
        break;
      }
    }
    if(ipConst != 0)flagIP = 1;
    if(flagIP == 1){
      _circle.add_point(0., 0., 0.5);
      //++_nPointsForFit;
    }
    for (unsigned i = 0; i < n; i++) {
        TMLink * l = _links[i];
        if (l == 0) continue;
 
        const TMDCWireHit * h = l->hit();
        if (h == 0) continue;
 
        //...Check next hit...
        HepPoint3D point;
        point = h->xyPosition();
        
        double weight = 1.0;
        weight = this->weight(* l);
 
        _circle.add_point(point.x(), point.y(), weight);
        //++_nPointsForFit;
    }
 
    //if (_nPointsForFit < 3) return -1;
    if(_circle.fit()<0.0 || _circle.kappa()==0.0) return -1;
    HepVector v(_circle.center());
    _center.setX(v(1));
    _center.setY(v(2));
    _radius = _circle.radius();
 
    //...Determine charge...Better way???
    int qSum = 0;
    for (unsigned i = 0; i < n; i++) {
        TMLink * l = _links[i];
        if (l == 0) continue;
 
        const TMDCWireHit * h = l->hit();
        if (h == 0) continue;
 
        float q = (_center.cross(h->xyPosition())).z();
        if (q > 0.) qSum += 1;
        else        qSum -= 1;
    }
    if (qSum >= 0) _charge = +1.;
    else           _charge = -1.;
    _radius *= _charge;
 
    _fitted = true;
    return 0;
}

impact()

double TCircle::impact ( void  ) const

inline

returns impact parameter to the origin.

Definition at line 136 of file TCircle.h.

                          {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::impact !!! fit not performed" << std::endl;
#endif
    return fabs(sqrt(_center.x() * _center.x() +
             _center.y() * _center.y())
        - fabs(_radius));
}

Referenced by dump(), and TMSelector::select().

objectType()

unsigned TCircle::objectType ( void  ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 165 of file TCircle.h.

                              {
    return Circle;
}

property()

void TCircle::property ( double  charge, double  radius, HepPoint3D  center  )

inline

sets circle properties.

Definition at line 157 of file TCircle.h.

                                                  {
    _charge = c;
    _radius = r;
    _center = e;
}

pt()

double TCircle::pt ( void  ) const

inline

returns Pt.

Definition at line 126 of file TCircle.h.

                      {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::pt !!! fit not performed" << std::endl;
#endif
   // return fabs(_radius) / Helix::ConstantAlpha;
    return fabs(_radius) / 333.564095;
}

Referenced by dump(), and TMSelector::select().

radius()

double TCircle::radius ( void  ) const

inline

returns radius.

Definition at line 117 of file TCircle.h.

                          {
#ifdef TRKRECO_DEBUG
    if (! _fitted) std::cout << "TCircle::radius !!! fit not performed" << std::endl;
#endif
    return _radius;
}

Referenced by TTrack::TTrack().

weight()

double TCircle::weight

(

const TMLink &

l

)

const

returns weight of TMLink in order to fit and make a circle.

Definition at line 111 of file TCircle.cxx.

                                      {
 
  //...Axial Wires
  int maxLink = 0;
  int localID[6];
  int layerID[6];
  int LayerID = l.hit()->wire()->layerId();
  for(int i=0;i<6;i++){
    if(l.neighbor(i)){
      maxLink = i;
      localID[i] = l.neighbor(i)->hit()->wire()->localId();
      layerID[i] = l.neighbor(i)->hit()->wire()->layerId();
    }else{
      break;
    }
  }
  if(maxLink != 1)return 1.0;
  if(layerID[0] == LayerID &&
     layerID[1] == LayerID)return 0.5;
  if(layerID[0]+layerID[1] != LayerID*2)return 1.0;
//Liuqg  if(localID[0] != localID[1])return 1.5;//1.0 or 2.0 ??
  return 1.0;
}

---

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

• TCircle.h
• TCircle.cxx