TCircleFitter Class Reference

A class to fit a TTrackBase object to a circle. More...

#include <TCircleFitter.h>

Inheritance diagram for TCircleFitter:

Public Member Functions
	TCircleFitter (const std::string &name)
	Constructor.
virtual	~TCircleFitter ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	fit (TTrackBase &) const
	TCircleFitter (const std::string &name)
	Constructor.
virtual	~TCircleFitter ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	fit (TTrackBase &) const
	TCircleFitter (const std::string &name)
	Constructor.
virtual	~TCircleFitter ()
	Destructor.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	fit (TTrackBase &) const
Public Member Functions inherited from TMFitter
	TMFitter (const std::string &name)
	Constructor.
virtual	~TMFitter ()
	Destructor.
const std::string &	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	TMFitter (const std::string &name)
	Constructor.
virtual	~TMFitter ()
	Destructor.
const std::string &	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	TMFitter (const std::string &name)
	Constructor.
virtual	~TMFitter ()
	Destructor.
const std::string &	name (void) const
	returns name.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const

Additional Inherited Members
Protected Member Functions inherited from TMFitter
void	fitDone (TTrackBase &) const
	sets the fitted flag. (Bad implementation)
void	fitDone (TTrackBase &) const
	sets the fitted flag. (Bad implementation)
void	fitDone (TTrackBase &) const
	sets the fitted flag. (Bad implementation)

Detailed Description

A class to fit a TTrackBase object to a circle.

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

Constructor & Destructor Documentation

TCircleFitter() [1/3]

TCircleFitter::TCircleFitter

(

const std::string &

name

)

Constructor.

Definition at line 21 of file TCircleFitter.cxx.

    : TMFitter( name ), _charge( 0. ), _radius( 0. ), _center( HepPoint3D( 0., 0., 0. ) ) {}

~TCircleFitter() [1/3]

TCircleFitter::~TCircleFitter ( )

virtual

Destructor.

Definition at line 24 of file TCircleFitter.cxx.

{}

TCircleFitter() [2/3]

TCircleFitter::TCircleFitter

(

const std::string &

name

)

Constructor.

~TCircleFitter() [2/3]

virtual TCircleFitter::~TCircleFitter ( )

virtual

Destructor.

TCircleFitter() [3/3]

TCircleFitter::TCircleFitter

(

const std::string &

name

)

Constructor.

~TCircleFitter() [3/3]

virtual TCircleFitter::~TCircleFitter ( )

virtual

Destructor.

Member Function Documentation

dump() [1/3]

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

dumps debug information.

dump() [2/3]

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

dumps debug information.

dump() [3/3]

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

dumps debug information.

fit() [1/3]

int TCircleFitter::fit ( TTrackBase & t ) const

virtual

Implements TMFitter.

Definition at line 26 of file TCircleFitter.cxx.

                                            {
 
#ifdef TRKRECO_DEBUG_DETAIL
  std::cout << "    TCircleFitter::fit ..." << std::endl;
#endif
 
  //...Already fitted ?...
  if ( t.fitted() ) return TFitAlreadyFitted;
 
  //...Check # of hits...
  if ( t.links().length() < 3 ) return TFitErrorFewHits;
 
  //...Hit loop...
  Lpav circle;
  /*
      circle.add_point(-173.011, 94.3169);
      circle.add_point(-187.6, 101.528);
      circle.add_point(-202.018, 108.601);
      circle.add_point(-216.639, 115.691);
      circle.add_point(-231.152, 122.49);
      circle.add_point(-246.207, 129.496);
      circle.add_point(-261.179, 136.407);
      circle.add_point(-275.198, 142.764);
      circle.add_point(-290.519, 149.572);
      circle.add_point(-304.42, 155.7);
      circle.add_point(-319.947, 162.285);
      circle.add_point(-335.495, 168.731);
      circle.add_point(-608.394, 264.577);
      circle.add_point(-612.842, 265.869);
      circle.add_point(-627.175, 270.011);
      circle.add_point(-643.654, 274.719);
      circle.add_point(-657.809, 278.683);
      circle.add_point(-674.464, 283.209);
      circle.add_point(-690.7, 287.541);
      circle.add_point(-704.925, 291.091);
    */
  unsigned n = t.links().length();
 
  for ( unsigned i = 0; i < n; i++ )
  {
    TMLink*            l = t.links()[i];
    const TMDCWireHit* h = l->hit();
 
    //...Check next hit...
    HepPoint3D point;
    if ( h->state() & WireHitPatternLeft ) point = h->position( WireHitLeft );
    else if ( h->state() & WireHitPatternRight ) point = h->position( WireHitRight );
    else point = h->xyPosition();
    // float weight = 1. / (h->distance() * h->distance());
    // float weight = 1. / h->distance();
 
    //  if (h->state() & WireHitPatternLeft) cout<<"h: "<<h->wire()->layerId()<<" local: "
    //                      <<h->wire()->localId()<<" TCirclefitter:
    // L"<<endl;    else cout<<"h: "<<h->wire()->layerId()<<" local: "
    //      <<h->wire()->localId()<<" TCirclefitter: R"<<endl;
 
    circle.add_point( point.x(), point.y() ); //, weight);
                                              // yuany
 
    //  cout<<"TCirclefitter("<<sqrt(point.x()*point.x()+point.y()*point.y())<<","
    //  <<point.y()/point.x()<<") dis:"
    //  <<(h->drift(0))*10000/40<<" point "<<point<<" L "<<(h->state()&WireHitPatternLeft)
    //  <<" R "<<(h->state()&WireHitPatternRight)<<endl;
 
    //  cout<<"TCirclefitter "<<point.x()<<"    "<<point.y()<<"    L
    //"<<(h->state()&WireHitPatternLeft)<<" R "<<(h->state()&WireHitPatternRight)<<endl;
  }
 
  if ( circle.fit() < 0.0 || circle.kappa() == 0.0 ) return TFitFailed;
  // yuany
  //    cout<<"circle chi2:"<<circle.fit()<<" kappa:"<<circle.kappa()<<"
  // center:"<<circle.center()<<" radius:"<<circle.radius()<<endl;
 
  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 = t.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;
 
  if ( t.objectType() == Circle ) ( (TCircle&)t ).property( _charge, _radius, _center );
  fitDone( t );
  return 0;
}

fit() [2/3]

virtual int TCircleFitter::fit ( TTrackBase & ) const

virtual

Implements TMFitter.

fit() [3/3]

virtual int TCircleFitter::fit ( TTrackBase & ) const

virtual

Implements TMFitter.

---

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

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