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.
	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.
	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	type (void) const
	returns type. Definition is depending on an object class.
const AList< TMLink > &	links (unsigned mask=0) const
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)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
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
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	type (void) const
	returns type. Definition is depending on an object class.
const AList< TMLink > &	links (unsigned mask=0) const
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)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
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
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	type (void) const
	returns type. Definition is depending on an object class.
const AList< TMLink > &	links (unsigned mask=0) const
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)
void	appendByDistance (AList< TMLink > &list, double maxDistance)
void	remove (TMLink &a)
	removes a TMLink.
void	remove (const AList< TMLink > &)
	removes TMLinks.
virtual void	refine (AList< TMLink > &list, double maxSigma)
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
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 InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

Constructor & Destructor Documentation

TCircle() [1/3]

TCircle::TCircle

(

const AList< TMLink > &

links

)

Constructor.

Definition at line 21 of file TCircle.cxx.

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

~TCircle() [1/3]

TCircle::~TCircle ( )

virtual

Destructor.

Definition at line 27 of file TCircle.cxx.

{}

TCircle() [2/3]

TCircle::TCircle

(

const AList< TMLink > &

links

)

Constructor.

~TCircle() [2/3]

virtual TCircle::~TCircle ( )

virtual

Destructor.

TCircle() [3/3]

TCircle::TCircle

(

const AList< TMLink > &

links

)

Constructor.

~TCircle() [3/3]

virtual TCircle::~TCircle ( )

virtual

Destructor.

Member Function Documentation

center() [1/3]

const HepPoint3D & TCircle::center ( void ) const

inline

returns position of center.

Definition at line 106 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

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

Referenced by property(), and TTrack::TTrack().

center() [2/3]

const HepPoint3D & TCircle::center

(

void

)

const

returns position of center.

center() [3/3]

const HepPoint3D & TCircle::center

(

void

)

const

returns position of center.

charge() [1/3]

double TCircle::charge ( void ) const

inline

returns charge.

Definition at line 136 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

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

Referenced by property().

charge() [2/3]

double TCircle::charge

(

void

)

const

returns charge.

charge() [3/3]

double TCircle::charge

(

void

)

const

returns charge.

dump() [1/3]

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 29 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().

dump() [2/3]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

dump() [3/3]

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

virtual

dumps debug information.

Reimplemented from TTrackBase.

fitForCurl() [1/3]

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

fitForCurl() [2/3]

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.

fitForCurl() [3/3]

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.

impact() [1/3]

double TCircle::impact ( void ) const

inline

returns impact parameter to the origin.

Definition at line 128 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/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().

impact() [2/3]

double TCircle::impact

(

void

)

const

returns impact parameter to the origin.

impact() [3/3]

double TCircle::impact

(

void

)

const

returns impact parameter to the origin.

objectType() [1/3]

unsigned TCircle::objectType ( void ) const

inlinevirtual

returns type.

Reimplemented from TTrackBase.

Definition at line 149 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

{ return Circle; }

objectType() [2/3]

virtual unsigned TCircle::objectType ( void ) const

virtual

returns type.

Reimplemented from TTrackBase.

objectType() [3/3]

virtual unsigned TCircle::objectType ( void ) const

virtual

returns type.

Reimplemented from TTrackBase.

property() [1/3]

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

inline

sets circle properties.

Definition at line 143 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

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

property() [2/3]

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

sets circle properties.

property() [3/3]

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

sets circle properties.

pt() [1/3]

double TCircle::pt ( void ) const

inline

returns Pt.

Definition at line 120 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/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().

pt() [2/3]

double TCircle::pt

(

void

)

const

returns Pt.

pt() [3/3]

double TCircle::pt

(

void

)

const

returns Pt.

radius() [1/3]

double TCircle::radius ( void ) const

inline

returns radius.

Definition at line 113 of file InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h.

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

Referenced by property(), and TTrack::TTrack().

radius() [2/3]

double TCircle::radius

(

void

)

const

returns radius.

radius() [3/3]

double TCircle::radius

(

void

)

const

returns radius.

weight() [1/3]

double TCircle::weight

(

const TMLink &

l

)

const

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

Definition at line 107 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;
}

Referenced by fitForCurl().

weight() [2/3]

double TCircle::weight

(

const TMLink &

l

)

const

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

weight() [3/3]

double TCircle::weight

(

const TMLink &

l

)

const

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

---

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

• InstallArea/x86_64-el9-gcc13-dbg/include/TrkReco/TCircle.h
• InstallArea/x86_64-el9-gcc13-opt/include/TrkReco/TCircle.h
• Reconstruction/TrkReco/include/TrkReco/TCircle.h
• TCircle.cxx