TRungeFitter Class Reference

A class to fit a TTrackBase object to a 3D line. More...

#include <TRungeFitter.h>

Inheritance diagram for TRungeFitter:

List of all members.

Public Member Functions
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
void	dump (const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
	dumps debug information.
virtual int	fit (TTrackBase &, float t0Offset) const
virtual int	fit (TTrackBase &) const
virtual int	fit (TTrackBase &, float t0Offset) const
virtual int	fit (TTrackBase &) const
const std::string &	name (void) const
	returns name.
const std::string &	name (void) const
	returns name.
void	propagation (int)
void	propagation (int)
void	sag (bool)
void	sag (bool)
void	tof (bool)
void	tof (bool)
	TRungeFitter (const std::string &name, bool m_sag, int m_prop, bool m_tof)
	TRungeFitter (const std::string &name)
	Constructor.
	TRungeFitter (const std::string &name, bool m_sag, int m_prop, bool m_tof)
	TRungeFitter (const std::string &name)
	Constructor.
virtual	~TRungeFitter ()
	Destructor.
virtual	~TRungeFitter ()
	Destructor.
Protected Member Functions
void	fitDone (TTrackBase &) const
	sets the fitted flag. (Bad implementation)
void	fitDone (TTrackBase &) const
	sets the fitted flag. (Bad implementation)
Private Attributes
int	_propagation
bool	_sag
bool	_tof

---

Detailed Description

A class to fit a TTrackBase object to a 3D line.

---

Constructor & Destructor Documentation

TRungeFitter::TRungeFitter (  const std::string &  name  )

Constructor. : TMFitter(name), _sag(true),_propagation(1),_tof(false){ }

TRungeFitter::TRungeFitter (  const std::string &  name, bool  m_sag, int  m_prop, bool  m_tof )

: TMFitter(name), _sag(m_sag),_propagation(m_prop),_tof(m_tof){ }

TRungeFitter::~TRungeFitter (   )  [virtual]

Destructor. { }

TRungeFitter::TRungeFitter (  const std::string &  name  )

Constructor.

TRungeFitter::TRungeFitter (  const std::string &  name, bool  m_sag, int  m_prop, bool  m_tof )

virtual TRungeFitter::~TRungeFitter (   )  [virtual]

Destructor.

---

Member Function Documentation

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

dumps debug information. Reimplemented from TMFitter.

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

dumps debug information. Reimplemented from TMFitter.

virtual int TRungeFitter::fit (  TTrackBase &  , float  t0Offset )  const [virtual]

virtual int TRungeFitter::fit (  TTrackBase &   )  const [virtual]

Implements TMFitter.

int TRungeFitter::fit (  TTrackBase &  , float  t0Offset )  const [virtual]

{ //std::cout<<"TRungeFitter::fit start"<<std::endl; //...Type check... if(tb.objectType() != Runge) return TFitUnavailable; TRunge& t = (TRunge&) tb; //...Already fitted ?... if(t.fitted()) return TFitAlreadyFitted; //...Count # of hits... AList<TMLink> cores = t.cores(); unsigned nCores = cores.length(); unsigned nStereoCores = NStereoHits(cores); //...Check # of hits... if ((nStereoCores < 2) || (nCores - nStereoCores < 3)) return TFitErrorFewHits; //...Move pivot to ORIGIN... const HepPoint3D pivot_bak = t.pivot(); t.pivot(ORIGIN); //...Setup... Vector a(5),da(5); a = t.a(); Vector dDda(5); Vector dchi2da(5); SymMatrix d2chi2d2a(5,0); const SymMatrix zero5(5,0); double chi2; double chi2Old = DBL_MAX; int err = 0; double factor = 0.1; Vector beta(5); SymMatrix alpha(5,0); SymMatrix alpha2(5,0); double (*d)[5]=new double[nCores][5]; // double (*d2)[5]=new double[nCores][5]; Vector ea(5); float tof; HepVector3D p; unsigned i; double distance; double eDistance; // ea... init const double ea_init=0.000001; ea[0]=ea_init; //dr ea[1]=ea_init; //phi0 ea[2]=ea_init; //kappa ea[3]=ea_init; //dz ea[4]=ea_init; //tanl //std::cout<<"TRF ::"<<a[0]<<","<<a[1]<<","<<a[2]<<","<<a[3]<<","<<a[4]<<std::endl; //long int lclock0=clock()/1000; //long int lclock=lclock0; Vector def_a(t.a()); Vector ta(def_a); //...Fitting loop... unsigned nTrial = 0; while(nTrial < 100){ //...Set up... chi2 = 0; for (unsigned j=0;j<5;j++) dchi2da[j]=0; d2chi2d2a=zero5; def_a=t.a(); //#### loop for shifted helix parameter set #### for(unsigned j=0;j<5;j++){ ta=def_a; ta[j]+=ea[j]; t.a(ta); //...Loop with hits... i=0; //std::cout<<"TRF:: cores="<<cores.length()<<std::endl; while(TMLink* l = cores[i++]){ //...Cal. closest points... t.approach(*l,tof,p,_sag); const HepPoint3D& onTrack=l->positionOnTrack(); const HepPoint3D& onWire=l->positionOnWire(); d[i-1][j]=(onTrack-onWire).mag(); //std::cout<<"TRF:: i="<<i<<std::endl; }//end of loop with hits //lclock=clock()/1000; //std::cout<<"TRF clock="<<lclock-lclock0<<std::endl; //lclock0=lclock; } /* for(int j=0;j<5;j++){ ta=def_a; ta[j]-=ea[j]; t.a(ta); //...Loop with hits... float tof_dummy; Vector3 p_dummy; unsigned i=0; while(TMLink* l = cores[i++]){ //...Cal. closest points... t.approach(*l,tof_dummy,p_dummy,_sag); const HepPoint3D& onTrack=l->positionOnTrack(); const HepPoint3D& onWire=l->positionOnWire(); d2[i-1][j]=(onTrack-onWire).mag(); }//end of loop with hits } */ t.a(def_a); //#### original parameter set and calc. chi2 #### //...Loop with hits... i=0; while(TMLink* l = cores[i++]){ const TMDCWireHit& h = *l->hit(); //...Cal. closest points... if(t.approach(*l,tof,p,_sag)<0){ std::cout<<"TrkReco:TRF:: bad wire"<<std::endl; continue; } const HepPoint3D& onTrack=l->positionOnTrack(); const HepPoint3D& onWire=l->positionOnWire(); unsigned leftRight = WireHitRight; if (onWire.cross(onTrack).z() < 0.) leftRight = WireHitLeft; //...Obtain drift distance and its error... if ((t0Offset == 0.) && (_propagation==0) && (! _tof)) { //...No correction... distance = l->drift(leftRight); eDistance = h.dDrift(leftRight); }else{ //...T0 and propagation corrections... int wire = h.wire()->id(); int side = leftRight; if (side==0) side = -1; float tp[3] = {p.x(),p.y(),p.z()}; float x[3] = {onWire.x(), onWire.y(), onWire.z()}; float time = h.reccdc()->tdc + t0Offset - tof; float dist; float edist; int prop = _propagation; //zsl calcdc_driftdist_(& prop, // & wire, // & side, // tp, // x, // & time, // & dist, // & edist); distance = (double) dist; eDistance = (double) edist; } double eDistance2=eDistance*eDistance; //...Residual... const double d0 = (onTrack-onWire).mag(); //if(d0>2) std::cout<<"TRF:: strange dist. d0="<<d0<<" x="<<distance // <<" ex="<<eDistance<<std::endl; double dDistance = d0 - distance; //...dDda... for(int j=0;j<5;j++){ dDda[j]=(d[i-1][j]-d0)/ea[j]; //if(dDda[j]==0) std::cout<<"TRF:: dDda=0 j="<<j<<" ea="<<ea[j]<<std::endl; } // for(int j=0;j<5;j++) dDda[j]=(d[i-1][j]-d2[i-1][j])/ea[j]/2.; //...Chi2 related... dchi2da += (dDistance / eDistance2) * dDda; d2chi2d2a += vT_times_v(dDda) / eDistance2; double pChi2 = dDistance * dDistance / eDistance2; chi2 += pChi2; //if(!(pChi2<0 || pChi2>=0)){ // std::cout<<"TRF:: pChi2="<<pChi2<<" X="<<d0<<" fx="<<distance // <<" ex="<<eDistance<<std::endl; //} //...Store results... l->update(onTrack, onWire, leftRight, pChi2); }//end of loop with hits //...Check condition... double change = chi2Old - chi2; if (0 <= change && change < 0.01) break; if (change < 0.) { factor *= 100; a += da; //recover if(-0.01 < change){ d2chi2d2a=alpha; chi2=chi2Old; break; } }else if(change > 0.){ chi2Old = chi2; factor *= 0.1; alpha=d2chi2d2a; beta=dchi2da; }else{ std::cout<<"TrkReco:TRF:: bad chi2 = "<<chi2<<std::endl; break; // protection for nan } alpha2=alpha; for(unsigned j=0;j<5;j++) alpha2[j][j]*=(1+factor); //...Cal. helix parameters for next loop... da = solve(alpha2,beta); //lclock=clock()/1000; //std::cout<<"TRF "<<nTrial<<": " // <<"cl="<<lclock-lclock0<<": " // <<a[0]<<","<<a[1]<<","<<a[2]<<","<<a[3]<<","<<a[4]<<" " // <<chi2<<"/"<<nCores<<":"<<factor // <<" :da "<<da[0]<<","<<da[1]<<","<<da[2]<<","<<da[3]<<","<<da[4]<<std::endl; //lclock0=lclock; a -= da; t.a(a); //ea = 0.0001*da; for(i=0;i<5;i++){ ea[i]=0.0001*abs(da[i]); if(ea[i]>ea_init) ea[i]=ea_init; if(ea[i]<1.0e-10) ea[i]=1.0e-10; } ++nTrial; } //std::cout<<"TRungeFitter:: nTrial="<<nTrial<<std::endl; //...Cal. error matrix... SymMatrix Ea(5,0); unsigned dim; dim=5; Ea = d2chi2d2a.inverse(err); //...Store information... if(! err){ t.a(a); t.Ea(Ea); t._fitted = true; }else{ err = TFitFailed; } t._ndf = nCores - dim; t._chi2 = chi2; //...Recover pivot... t.pivot(pivot_bak); delete [] d; // delete [] d2; return err; }

int TRungeFitter::fit (  TTrackBase &   )  const [virtual]

Implements TMFitter. { return fit(tb,0); }

void TMFitter::fitDone (  TTrackBase &   )  const [protected, inherited]

sets the fitted flag. (Bad implementation)

void TMFitter::fitDone (  TTrackBase &   )  const [protected, inherited]

sets the fitted flag. (Bad implementation) { t._fitted = true; }

const std::string& TMFitter::name (  void   )  const [inherited]

returns name.

const std::string & TMFitter::name (  void   )  const [inline, inherited]

returns name. { return _name; }

void TRungeFitter::propagation (  int   )

void TRungeFitter::propagation (  int   )

{ _propagation = _in; }

void TRungeFitter::sag (  bool   )

void TRungeFitter::sag (  bool   )

{ _sag = _in; }

void TRungeFitter::tof (  bool   )

void TRungeFitter::tof (  bool   )

{ _tof = _in; }

---

Member Data Documentation

int TRungeFitter::_propagation [private]

bool TRungeFitter::_sag [private]

bool TRungeFitter::_tof [private]

---

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

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/TrkReco/TrkReco/TRungeFitter.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Reconstruction/TrkReco/TrkReco-00-06-12/TrkReco/TRungeFitter.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Reconstruction/TrkReco/TrkReco-00-06-12/src/TRungeFitter.cxx

---