TRunge Class Reference

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

#include <TRunge.h>

Inheritance diagram for TRunge:

List of all members.

Public Member Functions
	TRunge ()
	Constructors.
	TRunge (const TTrack &)
	TRunge (const Helix &)
	TRunge (const TRunge &)
	TRunge (const RecMdcTrack &)
	~TRunge ()
	Destructor.
unsigned	objectType (void) const
	returns object type
double	dr (void) const
	Track parameters (at pivot).
double	phi0 (void) const
double	kappa (void) const
double	dz (void) const
double	tanl (void) const
const HepPoint3D &	pivot (void) const
	pivot position
const Vector &	a (void) const
	returns helix parameter
const SymMatrix &	Ea (void) const
	returns error matrix
Helix	helix (void) const
	returns helix class
unsigned	ndf (void) const
	returns NDF
double	chi2 (void) const
	returns chi2.
double	reducedchi2 (void) const
	returns reduced-chi2
int	BfieldID (void) const
	returns B field ID
double	StepSize (void) const
	returns step size
const double *	Yscal (void) const
	return error parameters for fitting with step size control
double	Eps (void) const
double	StepSizeMax (void) const
double	StepSizeMin (void) const
float	Mass (void) const
	return mass
double	MaxFlightLength (void) const
	return max flight length
int	approach (TMLink &, const RkFitMaterial, bool sagCorrection=true)
int	approach (TMLink &, float &tof, HepVector3D &p, const RkFitMaterial, bool sagCorrection=true)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, const RkFitMaterial material)
	caluculate closest points between a line and this track
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material)
int	approach_line (TMLink &, const HepPoint3D &, const HepVector3D &, HepPoint3D &onLine, HepPoint3D &onTrack, float &tof, HepVector3D &p, const RkFitMaterial material, unsigned &stepNum)
int	approach_point (TMLink &, const HepPoint3D &, HepPoint3D &onTrack, const RkFitMaterial material) const
	caluculate closest point between a point and this track
const HepPoint3D &	pivot (const HepPoint3D &)
	set new pivot
const Vector &	a (const Vector &)
	set helix parameter
const SymMatrix &	Ea (const SymMatrix &)
	set helix error matrix
int	BfieldID (int)
	set B field map ID
double	StepSize (double)
	set step size to calc. trajectory
const double *	Yscal (const double *)
	set error parameters for fitting with step size control
double	Eps (double)
double	StepSizeMax (double)
double	StepSizeMin (double)
float	Mass (float)
	set mass for tof calc.
double	MaxFlightLength (double)
double	DisToWire (TMLink &, double, HepPoint3D, HepPoint3D &)
double	dEpath (double, double, double) const
void	eloss (double path, const RkFitMaterial *material, double mass, double y[6], int index) const
unsigned	Fly (const RkFitMaterial material) const
	make the trajectory in cache, return the number of step
unsigned	Fly_SC (void) const
double	intersect_cylinder (double r) const
double	intersect_yz_plane (const HepTransform3D &plane, double x) const
double	intersect_zx_plane (const HepTransform3D &plane, double y) const
double	intersect_xy_plane (double z) const
void	Propagate (double y[6], const double &step, const RkFitMaterial material) const
	propagate the track using 4th-order Runge-Kutta method
void	Function (const double y[6], double f[6]) const
void	Propagate1 (const double y[6], const double dydx[6], const double &step, double yout[6]) const
void	Propagate_QC (double y[6], double dydx[6], const double &steptry, const double &eps, const double yscal[6], double &stepdid, double &stepnext) const
void	SetFirst (double y[6]) const
	set first point (position, momentum) s=0, phi=0
unsigned	Nstep (void) const
	access to trajectory cache
int	GetXP (unsigned stepNum, double y[6]) const
int	GetStep (unsigned stepNum, double &step) const
double	SetFlightLength (void)
	set flight length using wire hits
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).
Protected Attributes
AList< TMLink >	_links
bool	_fitted
bool	_fittedWithCathode
TTrackMC *	_mc
Private Attributes
HepPoint3D	_pivot
Vector	_a
SymMatrix	_Ea
double	_chi2
unsigned	_ndf
Bfield *	_bfield
int	_bfieldID
double	_maxflightlength
double	_stepSize
double	_yscal [6]
double	_eps
double	_stepSizeMax
double	_stepSizeMin
float	_mass
const RkFitMaterial	FitMdcMaterial
float	_mass2
int	_charge
double	_y [TRunge_MAXstep][6]
double	_h [TRunge_MAXstep]
unsigned	_Nstep
IMagneticFieldSvc *	m_pmgnIMF
Static Private Attributes
static const TRungeFitter	_fitter = TRungeFitter("TRunge Default fitter")
Friends
class	TRungeFitter

---

Detailed Description

A class to represent a track in tracking.

Definition at line 65 of file TRunge.h.

---

Constructor & Destructor Documentation

TRunge::TRunge

(

)

Constructors.

Definition at line 47 of file TRunge.cxx.

References _bfield, _bfieldID, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _mass, _mass2, _yscal, TTrackBase::fitter(), Bfield::getBfield(), and m_pmgnIMF.

  : TTrackBase(),
    _pivot(ORIGIN),_a(Vector(5,0)),_Ea(SymMatrix(5,0)), 
    _chi2(0),_ndf(0),_charge(1),_bfieldID(21),_Nstep(0),
    _stepSize(default_stepSize),_mass(0.140),_eps(EPS),
    _stepSizeMax(default_stepSizeMax),_stepSizeMin(default_stepSizeMin),
    _maxflightlength(default_maxflightlength) {

  //...Set a default fitter...
  fitter(& TRunge::_fitter);

  _fitted = false;
  _fittedWithCathode = false;

  _bfield = Bfield::getBfield(_bfieldID);

  _mass2=_mass*_mass;

  _yscal[0]=0.1; //1mm  -> error = 1mm*EPS
  _yscal[1]=0.1; //1mm
  _yscal[2]=0.1; //1mm
  _yscal[3]=0.001; //1MeV
  _yscal[4]=0.001; //1MeV
  _yscal[5]=0.001; //1MeV
  //jialk
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF); 
  if(scmgn!=StatusCode::SUCCESS) { 
    std::cout<< __FILE__<<" Unable to open Magnetic field service"<<std::endl;
  }

}

TRunge::TRunge

(

const TTrack &

)

Definition at line 114 of file TRunge.cxx.

References _a, _bfield, _bfieldID, _charge, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _mass, _mass2, _yscal, TTrackBase::fitter(), Bfield::getBfield(), and m_pmgnIMF.

  : TTrackBase(a),
    _pivot(a.helix().pivot()),_a(a.helix().a()),_Ea(a.helix().Ea()),
    _chi2(0),_ndf(0),_bfieldID(21),_Nstep(0),
    _stepSize(default_stepSize),_mass(0.140),_eps(EPS),
    _stepSizeMax(default_stepSizeMax),_stepSizeMin(default_stepSizeMin),
    _maxflightlength(default_maxflightlength) {
//  _a[2]=-_a[2];
  //...Set a default fitter...
  fitter(& TRunge::_fitter);

  _fitted = false;
  _fittedWithCathode = false;

  if(_a[2]<0) _charge=-1;
  else        _charge=1;

  _bfield = Bfield::getBfield(_bfieldID);

  _mass2=_mass*_mass;

  _yscal[0]=0.1; //1mm  -> error = 1mm*EPS
  _yscal[1]=0.1; //1mm
  _yscal[2]=0.1; //1mm
  _yscal[3]=0.001; //1MeV
  _yscal[4]=0.001; //1MeV
  _yscal[5]=0.001; //1MeV
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF);
  if(scmgn!=StatusCode::SUCCESS) {
          std::cout<< "Unable to open Magnetic field service"<<std::endl;
  }
  //SetFlightLength();
  //cout<<"TR:: _maxflightlength="<<_maxflightlength<<endl;
}

TRunge::TRunge

(

const Helix &

)

Definition at line 149 of file TRunge.cxx.

References _a, _bfield, _bfieldID, _charge, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _mass, _mass2, _yscal, TTrackBase::fitter(), and Bfield::getBfield().

  : TTrackBase(),
    _pivot(h.pivot()),_a(h.a()),_Ea(h.Ea()),
    _chi2(0),_ndf(0),_bfieldID(21),_Nstep(0),
    _stepSize(default_stepSize),_mass(0.140),_eps(EPS),
    _stepSizeMax(default_stepSizeMax),_stepSizeMin(default_stepSizeMin),
    _maxflightlength(default_maxflightlength) {

  //...Set a default fitter...
  fitter(& TRunge::_fitter);

  _fitted = false;
  _fittedWithCathode = false;

  if(_a[2]<0) _charge=-1;
  else        _charge=1;

  _bfield = Bfield::getBfield(_bfieldID);

  _mass2=_mass*_mass;

  _yscal[0]=0.1; //1mm  -> error = 1mm*EPS
  _yscal[1]=0.1; //1mm
  _yscal[2]=0.1; //1mm
  _yscal[3]=0.001; //1MeV
  _yscal[4]=0.001; //1MeV
  _yscal[5]=0.001; //1MeV
}

TRunge::TRunge

(

const TRunge &

)

Definition at line 178 of file TRunge.cxx.

References _a, _bfield, _bfieldID, _charge, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _mass, _mass2, _yscal, a(), TTrackBase::fitter(), Bfield::getBfield(), and genRecEmupikp::i.

  : TTrackBase((TTrackBase &) a),
    _pivot(a.pivot()),_a(a.a()),_Ea(a.Ea()),
    _chi2(a.chi2()),_ndf(a.ndf()),_bfieldID(a.BfieldID()),_Nstep(0),
    _stepSize(a.StepSize()),_mass(a.Mass()),_eps(a.Eps()),
    _stepSizeMax(a.StepSizeMax()),_stepSizeMin(a.StepSizeMin()),
    _maxflightlength(a.MaxFlightLength()) {

  //...Set a default fitter...
  fitter(& TRunge::_fitter);

  _fitted = false;
  _fittedWithCathode = false;

  if(_a[2]<0) _charge=-1;
  else        _charge=1;

  _bfield = Bfield::getBfield(_bfieldID);

  _mass2=_mass*_mass;

  for(unsigned i=0;i<6;i++) _yscal[i]=a.Yscal()[i];
}

TRunge::TRunge

(

const RecMdcTrack &

)

Definition at line 79 of file TRunge.cxx.

References _a, _bfield, _bfieldID, _charge, TTrackBase::_fitted, TTrackBase::_fittedWithCathode, _fitter, _mass, _mass2, _yscal, TTrackBase::fitter(), Bfield::getBfield(), and m_pmgnIMF.

  : TTrackBase(),
  _pivot(a.getPivot()),_a(a.helix()),_Ea(a.err()),
  _chi2(a.chi2()),_ndf(a.ndof()),_bfieldID(21),_Nstep(0),
  _stepSize(default_stepSize),_mass(0.140),_eps(EPS),
  _stepSizeMax(default_stepSizeMax),_stepSizeMin(default_stepSizeMin),
  _maxflightlength(default_maxflightlength) {
//    _a[2]=-_a[2];
    //...Set a default fitter...
    fitter(& TRunge::_fitter);

    _fitted = false;
    _fittedWithCathode = false;

    if(_a[2]<0) _charge=-1;
    else        _charge=1;

    _bfield = Bfield::getBfield(_bfieldID);

    _mass2=_mass*_mass;

    _yscal[0]=0.1; //1mm  -> error = 1mm*EPS
    _yscal[1]=0.1; //1mm
    _yscal[2]=0.1; //1mm
    _yscal[3]=0.001; //1MeV
    _yscal[4]=0.001; //1MeV
    _yscal[5]=0.001; //1MeV
    StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF);
    if(scmgn!=StatusCode::SUCCESS) {
      std::cout<< "Unable to open Magnetic field service"<<std::endl;
    }
    //SetFlightLength();
    //cout<<"TR:: _maxflightlength="<<_maxflightlength<<endl;

  }

TRunge::~TRunge

(

)

Destructor.

Definition at line 203 of file TRunge.cxx.

                {
}

---

Member Function Documentation

const Vector & TRunge::a

(

const Vector &

)

set helix parameter

Definition at line 299 of file TRunge.cxx.

References _a, _charge, and _Nstep.

                                       {
  _a=ta;
  if(_a[2]<0) _charge=-1;
  else        _charge=1;
  _Nstep=0;
  return(_a);
}

const Vector & TRunge::a

(

void

)

const

returns helix parameter

Definition at line 230 of file TRunge.cxx.

References _a.

Referenced by dEpath(), and TRunge().

                                 {
  return(_a);
}

void TTrackBase::append

(

const AList< TMLink > &

)

[inherited]

appends TMLinks.

Definition at line 384 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, TTrackBase::_updated, genRecEmupikp::i, and WireHitUsed.

                                          {
    AList<TMLink> tmp;
    for (unsigned i = 0; i < a.length(); i++) {
        if ((_links.hasMember(a[i])) || (a[i]->hit()->state() & WireHitUsed)) {
#ifdef TRKRECO_DEBUG_DETAIL
            std::cout << "    TTrackBase::append(list) !!! ";
            std::cout << a[i]->wire()->name();
            std::cout << " Hey!, this is already used! Don't mess with me.";
            std::cout << std::endl;
#endif
            continue;
        }
        else {
            tmp.append(a[i]);
        }
    }
    _links.append(tmp);
    _updated = false;
    _fitted = false;
    _fittedWithCathode = false; // added by matsu ( 1999/05/24 )
}

void TTrackBase::append

(

TMLink &

)

[inherited]

appends a TMLink.

Definition at line 362 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, TTrackBase::_updated, TMLink::hit(), TMDCWire::name(), TMDCWireHit::state(), TMLink::wire(), and WireHitUsed.

Referenced by TBuilderCosmic::buildStereo(), TBuilder0::buildStereo(), TBuilder::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), TBuilder::buildStereoNew(), THistogram::clusters(), THistogram::clusters0(), TCurlFinder::dividing2DTrack(), TrkReco::execute(), TConformalFinder::link(), TCurlFinder::make2DTrack(), TTrackManager::merge(), TTrackManager::salvageAssociateHits(), TSegment0::splitAV(), TSegment::splitAV(), TSegment0::splitParallel(), TSegment::splitParallel(), and TConformalFinder0::standardFinding().

                             {
    if ((a.hit()->state() & WireHitUsed)) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "TTrackBase::append !!! " << a.wire()->name()
                  << " this is already used by another track!" << std::endl;
#endif
        return;
    }
    if (_links.hasMember(a)) {
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "TTrackBase::append !!! " << a.wire()->name()
                  << " this is already included in this track!" << std::endl;
#endif
        return;
    }
    _links.append(a);
    _updated = false;
    _fitted = false;
    _fittedWithCathode = false; // added by matsu ( 1999/05/24 )
}

void TTrackBase::appendByApproach	(	AList< TMLink > &	list,
		double	maxSigma
	)			[inherited]

appends TMLinks by approach. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.

Definition at line 101 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, TTrackBase::approach(), TTrackBase::distance(), showlog::err, genRecEmupikp::i, and WireHitUsed.

Referenced by TBuilder0::appendClusters(), and TTrackManager::salvage().

                                                                  {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    TTrackBase::appendByApproach ... sigma=" << maxSigma << std::endl;
#endif

    AList<TMLink> unused;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
        TMLink & l = * list[i];

        if ((_links.hasMember(l)) || (l.hit()->state() & WireHitUsed))
            continue;

        //...Calculate closest approach...
        int err = approach(l);
        if (err < 0) {
            unused.append(l);
            continue;
        }

        //...Calculate sigma...
        float distance = (l.positionOnWire() - l.positionOnTrack()).mag();
        float diff = fabs(distance - l.drift());
        float sigma = diff / l.dDrift();

        //...For debug...
#ifdef TRKRECO_DEBUG_DETAIL
        std::cout << "    sigma=" << sigma;
        std::cout << ",dist=" << distance;
        std::cout << ",diff=" << diff;
        std::cout << ",err=" << l.hit()->dDrift() << ",";
        if (sigma < maxSigma) std::cout << "ok,";
        else                  std::cout << "X,";
        l.dump("mc");
#endif
                
        //...Make sigma cut...
        if (sigma > maxSigma) {
            unused.append(l);
            continue;
        }

        //...OK...
        _links.append(l);
        _updated = false;
        _fitted = false;
    }
    list.remove(unused);
}

void TTrackBase::appendByDistance	(	AList< TMLink > &	list,
		double	maxDistance
	)			[inherited]

appends TMLinks by distance. 'list' is an input. Unappended TMLinks will be removed from 'list' when returned.

Definition at line 152 of file TTrackBase.cxx.

                                                               {
    std::cout << "TTrackBase::appendByDistance !!! not implemented" << std::endl;
    list.removeAll();
}

int TTrackBase::approach

(

TMLink &

)

const [virtual, inherited]

calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.

Reimplemented in TTrack.

Definition at line 95 of file TTrackBase.cxx.

Referenced by TTrackBase::appendByApproach(), and TTrackBase::testByApproach().

                                   {
    std::cout << "TTrackBase::approach !!! not implemented" << std::endl;
    return -1;
}

int TRunge::approach	(	TMLink &	,
		float &	tof,
		HepVector3D &	p,
		const	RkFitMaterial,
		bool	sagCorrection = true
	)

Definition at line 366 of file TRunge.cxx.

References _pivot, approach_line(), TMLink::positionOnTrack(), TMLink::positionOnWire(), v, w, and TMLink::wire().

                                                                         {
  HepPoint3D xw;
  HepPoint3D wireBackwardPosition;
  HepVector3D v;
  HepPoint3D onWire,onTrack;
  double onWire_x,onWire_y, onWire_z, zwf, zwb;
  const TMDCWire& w=*l.wire();
  xw = w.xyPosition();
  wireBackwardPosition = w.backwardPosition();
  v = w.direction();

  unsigned stepNum=0;
  if(approach_line(l,wireBackwardPosition,v,onWire,onTrack,tof,p,material,stepNum)<0){
    //cout<<"TR::error approach_line"<<endl;
    return(-1);
  }
  zwf = w.forwardPosition().z();
  zwb = w.backwardPosition().z();
  // protection for strange wire hit info. (Only 'onWire' is corrected)
  if(onWire.z() > zwf) 
    w.wirePosition(zwf,onWire,wireBackwardPosition,(HepVector3D&)v);
  else if(onWire.z() < zwb) 
    w.wirePosition(zwb,onWire,wireBackwardPosition,(HepVector3D&)v);

  // onWire,onTrack filled
  
  if(!doSagCorrection){
    l.positionOnWire(onWire);
    l.positionOnTrack(onTrack);
    double phipoint=atan((onTrack.y()-_pivot.y())/onTrack.x()-_pivot.x());
//    cout<<"dphi track:"<<l.dPhi()<<endl;
 //   cout<<"dphi rk:"<<phipoint-_a[0]<<endl;
//    l.dPhi(phipoint-_a[0]);
    return(0);                          // no sag correction
  }
  // Sag correction
  //   loop for sag correction
  onWire_y = onWire.y();
  onWire_z = onWire.z();

  unsigned nTrial = 1;
  while(nTrial<100){
    //cout<<"TR: nTrial "<<nTrial<<endl;
    w.wirePosition(onWire_z,xw,wireBackwardPosition,(HepVector3D&)v);
    if(approach_line(l,wireBackwardPosition,v,onWire,onTrack,tof,p,material,stepNum)<0)
      return(-1);
    if(fabs(onWire_y - onWire.y())<0.0001) break;       // |dy|< 1 micron
    onWire_y = onWire.y();
    onWire_z += (onWire.z()-onWire_z)/2;
    // protection for strange wire hit info.
    if(onWire_z > zwf)      onWire_z=zwf;
    else if(onWire_z < zwb) onWire_z=zwb;

    nTrial++;
  }
  //  cout<<"TR  nTrial="<<nTrial<<endl;
  l.positionOnWire(onWire);
  l.positionOnTrack(onTrack);
  return(nTrial);
}

int TRunge::approach	(	TMLink &	,
		const	RkFitMaterial,
		bool	sagCorrection = true
	)

calculates the closest approach to a wire in real space. Results are stored in TMLink. Return value is negative if error happened.

Definition at line 360 of file TRunge.cxx.

                                                                                 {
  float tof;
  HepVector3D p;
  return(approach(l,tof,p,material,doSagCorrection));
}

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material,
		unsigned &	stepNum
	)

Definition at line 442 of file TRunge.cxx.

References _a, _h, _mass2, _Nstep, _stepSize, _y, EvtCyclic3::A, DBL_MAX, DisToWire(), dz(), Fly(), Fly_SC(), Function(), genRecEmupikp::i, ganga-rec::j, Propagate(), s, and t().

                                                                                                      {
  //  line = [w0] + t * [v]    -> [onLine]
  if(_Nstep==0){
    if(_stepSize==0) Fly_SC();
    else Fly(material);

  }
  //cout<<"TR::approach_line stepNum="<<stepNum<<endl;

  const double w0x = w0.x();
  const double w0y = w0.y();
  const double w0z = w0.z();
  //cout<<"TR::line w0="<<w0x<<","<<w0y<<","<<w0z<<endl;
  const double vx = v.x();
  const double vy = v.y();
  const double vz = v.z();
  const double v_2 = vx*vx+vy*vy+vz*vz;

  const float clight=29.9792458; //[cm/ns]   
  //const float M2=_mass*_mass;
  //const float& M2=_mass2;
  const float p2 = _y[0][3]*_y[0][3]+_y[0][4]*_y[0][4]+_y[0][5]*_y[0][5];
  const float tof_factor=1./clight*sqrt(1+_mass2/p2);

  // search for the closest point in cache   (point - line)
  //  unsigned stepNum;
  double l2_old= DBL_MAX;
  if(stepNum>_Nstep) stepNum=0;
  unsigned stepNumLo;
  unsigned stepNumHi;
  if(stepNum==0 && _stepSize!=0){ // skip
    const double dx = _y[0][0]-w0x;
    const double dy = _y[0][1]-w0y;
    stepNum=(unsigned)
      (sqrt( (dx*dx+dy*dy)*(1+_a[4]*_a[4]) )/_stepSize);
  }
  unsigned mergin;
  if(_stepSize==0){
    mergin=10;  //10 step back
  }else{
    mergin=(unsigned)(1.0/_stepSize);   //1mm back
  }
  if(stepNum>mergin) stepNum-=mergin;
  else           stepNum=0;
  if(stepNum>=_Nstep) stepNum=_Nstep-1;
  // hunt
  //  unsigned inc=1;
  unsigned inc=(mergin>>1)+1;
  stepNumLo=stepNum;
  stepNumHi=stepNum;
  for(;;){
    const double dx = _y[stepNumHi][0]-w0x;
    const double dy = _y[stepNumHi][1]-w0y;
    const double dz = _y[stepNumHi][2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
//    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zvec;
    zvec.setX(0);
    zvec.setY(0);
    zvec.setZ(1);
    HepVector3D hitv=t*v;
    HepPoint3D pp,onw;
    pp.setX(_y[stepNumHi][0]);
    pp.setY(_y[stepNumHi][1]);
    pp.setZ(_y[stepNumHi][2]); 
    double zwire=hitv.dot(zvec)+w0z;
    double dtl=DisToWire(l,zwire,pp,onw);
    const double l2=dtl*dtl;      
    if(l2 > l2_old) break;
    l2_old=l2;
    stepNumLo=stepNumHi;
    //    inc+=inc;
    stepNumHi+=inc;
    if(stepNumHi>=_Nstep){
      stepNumHi=_Nstep;
      break;
    }
  }
  // locate (2-bun-hou, bisection method)
  while(stepNumHi-stepNumLo>1){
    unsigned j=(stepNumHi+stepNumLo)>>1;
    const double dx = _y[j][0]-w0x;
    const double dy = _y[j][1]-w0y;
    const double dz = _y[j][2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
//    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zv;
    zv.setX(0);
    zv.setY(0);
    zv.setZ(1);
    HepVector3D hitv=t*v;
    HepPoint3D point,onwir;
    point.setX(_y[j][0]);
    point.setY(_y[j][1]);
    point.setZ(_y[j][2]);
    double zwir=hitv.dot(zv)+w0z;
    double dtll=DisToWire(l,zwir,point,onwir);
    const double l2 =dtll*dtll;
    if(l2 > l2_old){
      stepNumHi=j;
    }else{
      l2_old=l2;
      stepNumLo=j;
    }
  }
  //stepNum=stepNumHi;
  stepNum=stepNumLo;
  /*
  for(;stepNum<_Nstep;stepNum++){
    const double dx = _y[stepNum][0]-w0x;
    const double dy = _y[stepNum][1]-w0y;
    const double dz = _y[stepNum][2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    if(l2 > l2_old) break;
    l2_old=l2;
    //    const float p2 = _y[stepNum][3]*_y[stepNum][3]+
    //                      _y[stepNum][4]*_y[stepNum][4]+
    //                      _y[stepNum][5]*_y[stepNum][5];
    //    tof+=_stepSize/clight*sqrt(1+M2/p2);
  }
  */
  //  cout<<"TR  stepNum="<<stepNum<<endl;
  //if(stepNum>=_Nstep) return(-1);     // not found
  //stepNum--;
  if(_stepSize==0){
    double dstep=0;
    for(unsigned i=0;i<stepNum;i++) dstep+=_h[i];
    tof=dstep*tof_factor;
  }else{
    tof=stepNum*_stepSize*tof_factor;
  }

  // propagate the track and search for the closest point
  //2-bun-hou (bisection method)
  /*
  double y[6],y_old[6];
  for(unsigned i=0;i<6;i++) y[i]=_y[stepNum][i];
  double step=_stepSize;
  double step2=0;
  for(;;){
    for(unsigned i=0;i<6;i++) y_old[i]=y[i];
    Propagate(y,step);
    const double dx = y[0]-w0x;
    const double dy = y[1]-w0y;
    const double dz = y[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    if(l2 > l2_old){  // back
      for(unsigned i=0;i<6;i++) y[i]=y_old[i];
    }else{            // propagate
      l2_old=l2;
      //      const float p2=y[3]*y[3]+y[4]*y[4]+y[5]*y[5];
      //      tof+=step/clight*sqrt(1+M2/p2);
      step2+=step;
    }
    step/=2;
    if(step < 0.0001) break;    // step < 1 [um]
  }
  */
  // Hasamiuchi-Hou (false position method)
  /*
  double y[6],y1[6],y2[6];
  for(unsigned i=0;i<6;i++) y1[i]=_y[stepNum][i];
  for(unsigned i=0;i<6;i++) y2[i]=_y[stepNum+2][i];
  double minStep=0;
  double maxStep=_stepSize*2;
  double step2=0;
  const double A[3][3]={{1-vx*vx/v_2, -vx*vy/v_2, -vx*vz/v_2},
                        { -vy*vx/v_2,1-vy*vy/v_2, -vy*vz/v_2},
                        { -vz*vx/v_2, -vz*vy/v_2,1-vz*vz/v_2}};

  double Y1=0;
  {
    const double dx = y1[0]-w0x;
    const double dy = y1[1]-w0y;
    const double dz = y1[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
    const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
    const double pmag=sqrt(y1[3]*y1[3]+y1[4]*y1[4]+y1[5]*y1[5]);
    for(int j=0;j<3;j++){
      double g=0;
      for(int k=0;k<3;k++) g+=A[j][k]*d[k];
      Y1+=y1[j+3]*g;
    }
    Y1=Y1/pmag/l;
  }
  double Y2=0;
  {
    const double dx = y2[0]-w0x;
    const double dy = y2[1]-w0y;
    const double dz = y2[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
    const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
    const double pmag=sqrt(y2[3]*y2[3]+y2[4]*y2[4]+y2[5]*y2[5]);
    for(int j=0;j<3;j++){
      double g=0;
      for(int k=0;k<3;k++) g+=A[j][k]*d[k];
      Y2+=y2[j+3]*g;
    }
    Y2=Y2/pmag/l;
  }
  if(Y1*Y2>=0) cout<<"TR  fatal error!"<<endl;
  double step_old= DBL_MAX;
  for(;;){
    const double step=(minStep*Y2-maxStep*Y1)/(Y2-Y1);
    if(fabs(step-step_old)<0.0001) break;
    step_old=step;
    for(unsigned i=0;i<6;i++) y[i]=y1[i];
    Propagate(y,step-minStep);
    double Y=0;
    {
      const double dx = y[0]-w0x;
      const double dy = y[1]-w0y;
      const double dz = y[2]-w0z;
      const double t = (dx*vx+dy*vy+dz*vz)/v_2;
      const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
      const double l = sqrt( (dx*dx+dy*dy+dz*dz)-(t*t)/v_2 );
      const double pmag=sqrt(y[3]*y[3]+y[4]*y[4]+y[5]*y[5]);
      for(int j=0;j<3;j++){
        double g=0;
        for(int k=0;k<3;k++) g+=A[j][k]*d[k];
        Y+=y[j+3]*g;
      }
      Y=Y/pmag/l;
    }
    if(Y1*Y<0){ //Y->Y2
      Y2=Y;
      for(unsigned i=0;i<6;i++) y2[i]=y[i];
      maxStep=step;
    }else{      //Y->Y1
      Y1=Y;
      for(unsigned i=0;i<6;i++) y1[i]=y[i];
      minStep=step;
    }
    if(Y1==Y2) break;
  }
  step2=step_old;
  */
  // Newton method
  double y[6];
  for(unsigned i=0;i<6;i++) y[i]=_y[stepNum][i];
  //memcpy(y,_y[stepNum],sizeof(double)*6);
  double step2=0;
  const double A[3][3]={{1-vx*vx/v_2, -vx*vy/v_2, -vx*vz/v_2},
                        { -vy*vx/v_2,1-vy*vy/v_2, -vy*vz/v_2},
                        { -vz*vx/v_2, -vz*vy/v_2,1-vz*vz/v_2}};
  double factor=1;
  double g[3];
  double f[6];
  double Af[3],Af2[3];
  unsigned i,j,k;
  HepPoint3D point,onwir;
  for(i=0;i<10;i++){
    //cout<<"TR::line "<<y[0]<<","<<y[1]<<","<<y[2]<<","<<y[3]<<","<<y[4]<<","<<y[5]<<endl;
    const double dx = y[0]-w0x;
    const double dy = y[1]-w0y;
    const double dz = y[2]-w0z;
    const double t = (dx*vx+dy*vy+dz*vz)/v_2;
    const double d[3]={dx-t*vx,dy-t*vy,dz-t*vz};
   // const double l2 = (dx*dx+dy*dy+dz*dz)-(t*t)/v_2;
    HepVector3D zv;
    zv.setX(0);
    zv.setY(0);
    zv.setZ(1);
    HepVector3D hitv=t*v;
    point.setX(y[0]);
    point.setY(y[1]);
    point.setZ(y[2]);
    double zwir=hitv.dot(zv)+w0z;
    double dtll=DisToWire(l,zwir,point,onwir);
    const double l2=dtll*dtll;
    for(j=0;j<3;j++){
      g[j]=0;
      for(k=0;k<3;k++) g[j]+=A[j][k]*d[k];
    }
    //cout<<"g="<<g[0]<<","<<g[1]<<","<<g[2]<<endl;
    Function(y,f);
    //cout<<"f="<<f[0]<<","<<f[1]<<","<<f[2]<<","<<f[3]<<","<<f[4]<<","<<f[5]<<endl;
    //cout<<"A="<<A[0][0]<<","<<A[0][1]<<","<<A[0][2]<<endl;
    //cout<<"A="<<A[1][0]<<","<<A[1][1]<<","<<A[1][2]<<endl;
    //cout<<"A="<<A[2][0]<<","<<A[2][1]<<","<<A[2][2]<<endl;
    double Y=0;
    for(j=0;j<3;j++) Y+=y[j+3]*g[j];
    double dYds=0;
    for(j=0;j<3;j++) dYds+=f[j+3]*g[j];
    //cout<<"dYds="<<dYds<<endl;
    for(j=0;j<3;j++){
      Af[j]=0;
      Af2[j]=0;
    }
    for(j=0;j<3;j++){
      //Af[j]=0;
      //Af2[j]=0;
      for(k=0;k<3;k++) Af[j]+=(A[j][k]*f[k]);
      for(k=0;k<3;k++) Af2[j]+=(A[j][k]*Af[k]);
      dYds+=y[j+3]*Af2[j];
      //cout<<j<<" dYds="<<dYds<<endl;
    }
    //cout<<"dYds="<<dYds<<endl;
    const double step=-Y/dYds*factor;
    //cout<<"TR  step="<<step<<" i="<<i<<endl;
    if(fabs(step) < 0.00001) break;     // step < 1 [um]
    if(l2 > l2_old) factor/=2;
    l2_old=l2;
    Propagate(y,step,material);
    step2+=step;
  }

  tof+=step2*tof_factor;

  onTrack.setX(y[0]);
  onTrack.setY(y[1]);
  onTrack.setZ(y[2]);
  p.setX(y[3]);
  p.setY(y[4]);
  p.setZ(y[5]);

  const double dx = y[0]-w0x;
  const double dy = y[1]-w0y;
  const double dz = y[2]-w0z;
  const double s = (dx*vx+dy*vy+dz*vz)/v_2;
 // onLine.setX(w0x+s*vx);
 // onLine.setY(w0y+s*vy);
 // onLine.setZ(w0z+s*vz);
  onLine.setX(onwir.x());
  onLine.setY(onwir.y());
  onLine.setZ(onwir.z());
  return(0);
}

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		float &	tof,
		HepVector3D &	p,
		const RkFitMaterial	material
	)

Definition at line 435 of file TRunge.cxx.

References approach_line().

                                                                                   {
  unsigned stepNum=0;
  return(approach_line(l,w0,v,onLine,onTrack,tof,p,material ,stepNum));
}

int TRunge::approach_line	(	TMLink &	,
		const HepPoint3D &	,
		const HepVector3D &	,
		HepPoint3D &	onLine,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material
	)

caluculate closest points between a line and this track

Definition at line 428 of file TRunge.cxx.

Referenced by approach(), and approach_line().

                                                                                                {
  float tof;
  HepVector3D p;
  return(approach_line(l,w0,v,onLine,onTrack,tof,p,material));
}

int TRunge::approach_point	(	TMLink &	,
		const HepPoint3D &	,
		HepPoint3D &	onTrack,
		const RkFitMaterial	material
	)			const

caluculate closest point between a point and this track

Definition at line 777 of file TRunge.cxx.

References _Nstep, _stepSize, _y, DBL_MAX, Fly(), Fly_SC(), genRecEmupikp::i, and Propagate().

                                                                                                                {
  if(_Nstep==0){
    if(_stepSize==0) Fly_SC();
    else Fly(material);
  }

  double x0=p0.x();
  double y0=p0.y();
  double z0=p0.z();

  //tof=0;
  //const float clight=29.9792458; //[cm/ns]   
  //const float M2=_mass*_mass;

  // search for the closest point in cache
  unsigned stepNum;
  double l2_old= DBL_MAX;
  for(stepNum=0;stepNum<_Nstep;stepNum++){
    double l2=(_y[stepNum][0]-x0)*(_y[stepNum][0]-x0)+
             (_y[stepNum][1]-y0)*(_y[stepNum][1]-y0)+
             (_y[stepNum][2]-z0)*(_y[stepNum][2]-z0);
    if(l2 > l2_old) break;
    l2_old=l2;
    //const double p2 = _y[stepNum][3]*_y[stepNum][3]+
    //                  _y[stepNum][4]*_y[stepNum][4]+
    //                  _y[stepNum][5]*_y[stepNum][5];
    //tof+=_stepSize/clight*sqrt(1+M2/p2);
  }
  if(stepNum>=_Nstep) return(-1);       // not found
  stepNum--;

  // propagate the track and search for the closest point
  double y[6],y_old[6];
  for(unsigned i=0;i<6;i++) y[i]=_y[stepNum][i];
  double step=_stepSize;
  for(;;){
    for(unsigned i=0;i<6;i++) y_old[i]=y[i];
    Propagate(y,step,material);
    double l2=(y[0]-x0)*(y[0]-x0)+(y[1]-y0)*(y[1]-y0)+(y[2]-z0)*(y[2]-z0);
    if(l2 > l2_old){  // back
      for(unsigned i=0;i<6;i++) y[i]=y_old[i];
    }else{            // propagate
      l2_old=l2;
      //const double p2=y[3]*y[3]+y[4]*y[4]+y[5]*y[5];
      //tof+=step/clight*sqrt(1+M2/p2);
    }
    step/=2;
    if(step < 0.0001) break;    // step < 1 [um]
  }

  onTrack.setX(y[0]);
  onTrack.setY(y[1]);
  onTrack.setZ(y[2]);
  //p.setX(y[3]);
  //p.setY(y[4]);
  //p.setZ(y[5]);
  return(0);
}

int TRunge::BfieldID

(

int

)

set B field map ID

Definition at line 313 of file TRunge.cxx.

References _bfield, _bfieldID, _Nstep, and Bfield::getBfield().

                          {
  _bfieldID=id;
  _bfield = Bfield::getBfield(_bfieldID);
  _Nstep=0;
  return(_bfieldID);
}

int TRunge::BfieldID

(

void

)

const

returns B field ID

Definition at line 258 of file TRunge.cxx.

References _bfieldID.

                              {
  return(_bfieldID);
}

double TRunge::chi2

(

void

)

const

returns chi2.

Definition at line 246 of file TRunge.cxx.

References _chi2.

                             {
  return _chi2;
}

const AList< TMLink > & TTrackBase::cores

(

unsigned

mask = 0

)

const [inherited]

returns a list of masked TMLinks for fit. 'mask' will be applied if mask is not 0.

Definition at line 317 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_updated, and TTrackBase::update().

Referenced by TBuilder::buildStereoNew(), TConformalFinder::expand(), TTrack::movePivot(), SetFlightLength(), and TTrack::szPosition().

                                     {
    if (mask)
        std::cout << "TTrackBase::cores !!! mask is not supported" << std::endl;
    if (! _updated) update();
    return _cores;
}

double TRunge::dEpath	(	double	,
		double	,
		double
	)			const

Definition at line 1239 of file TRunge.cxx.

References a(), f1, f2, genRecEmupikp::i, s, w, and x.

                                                            {
       double z=4.72234;
       double a=9.29212;
       double Density=0.00085144;
       double coeff=Density*z/a;
       double i=51.4709;
       double isq=i * i * 1e-18;
//       double sigma0_2 = 0.1569*coeff*path;
//  
//  if (sigma0_2<0) return 0;
// 
//  double psq = p * p;
//  double bsq = psq / (psq + mass * mass);
//  
//  // Correction for relativistic particles :
//  double sigma_2 = sigma0_2*(1-0.5*bsq)/(1-bsq);
//
//  if (sigma_2<0) return 0;
//
//  // Return sigma in GeV !!
//  return sqrt(sigma_2)*0.001;
const double Me = 0.000510999;
  double psq = p * p;
  double bsq = psq / (psq + mass * mass);
  double esq = psq / (mass * mass);

  double s = Me / mass;
  double w = (4 * Me * esq
          / (1 + 2 * s * sqrt(1 + esq)
         + s * s));

  // Density correction :
  double cc, x0;
  cc = 1+2*log(sqrt(isq)/(28.8E-09*sqrt(coeff)));
  if (cc < 5.215)
    x0 = 0.2;
  else
    x0 = 0.326*cc-1.5;
  double x1(3), xa(cc/4.606), aa;
  aa = 4.606*(xa-x0)/((x1-x0)*(x1-x0)*(x1-x0));
  double delta(0);
double x(log10(sqrt(esq)));
  if (x > x0){
    delta = 4.606*x-cc;
    if (x < x1) delta=delta+aa*(x1-x)*(x1-x)*(x1-x);
  }

 // Shell correction :  
  float f1, f2, f3, f4, f5, ce;
  f1 = 1/esq;
  f2 = f1*f1;
  f3 = f1*f2;
  f4 = (f1*0.42237+f2*0.0304-f3*0.00038)*1E12;
  f5 = (f1*3.858-f2*0.1668+f3*0.00158)*1E18;
  ce = f4*isq+f5*isq*sqrt(isq);

  return (0.0001535 * coeff / bsq
      * (log(Me * esq * w / isq)
         - 2 * bsq-delta-2.0*ce/z)) * path;
}

double TTrackBase::distance

(

const TMLink &

)

const [virtual, inherited]

returns distance to a position of TMLink in TMLink space.

Reimplemented in TLine0, and TMLine.

Definition at line 89 of file TTrackBase.cxx.

Referenced by TTrackBase::appendByApproach(), TTrack::fit2D(), and TTrackBase::testByApproach().

                                         {
    std::cout << "TTrackBase::distance !!! not implemented" << std::endl;
    return 0.;
}

double TRunge::DisToWire	(	TMLink &	,
		double	,
		HepPoint3D	,
		HepPoint3D &
	)

Definition at line 1215 of file TRunge.cxx.

References check_raw_filter::dist, w, and TMLink::wire().

Referenced by approach_line().

                                                                        {// by liucy 2011/7/28
   double zinter[400];
   double ddist=999;
   double finalz=0;
   const TMDCWire& w=*l.wire();
   HepPoint3D onwire;
//   for(int i=0;i<20;i++){
  //   zinter[i]=z-0.0002*i+0.002;
    HepPoint3D point=w.xyPosition(z);
//     HepPoint3D point=w.xyPosition(zinter[i]);
     onwire.setX(point.x());
     onwire.setY(point.y());
   //  onwire.setZ(zinter[i]);
     onwire.setZ(z);
//     double dist=sqrt((point.x()-a.x())*(point.x()-a.x())+(point.y()-a.y())*(point.y()-a.y())+(zinter[i]-a.z())*(zinter[i]-a.z()));
      double dist=sqrt((point.x()-a.x())*(point.x()-a.x())+(point.y()-a.y())*(point.y()-a.y())+(z-a.z())*(z-a.z()));
     if(dist<ddist){
       ddist=dist;
//       finalz=zinter[i];
       b=onwire;
     }
//   }   
   return ddist;
}

double TRunge::dr

(

void

)

const

Track parameters (at pivot).

Definition at line 206 of file TRunge.cxx.

References _a.

                           {
  return(_a[0]);
}

int TTrackBase::DropWorst

(

)

[virtual, inherited]

Definition at line 212 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, and genRecEmupikp::i.

Referenced by TrkReco::execute().

                      {
  AList<TMLink> bad;
  int jbad=-1;
  double sigma=0;
  unsigned n = _links.length();
  double chi2=0;
  for (unsigned i = 0; i < n; i++){
   // if(g_sigmalyr1)g_sigmalyr1->fill(sqrt(_links[i]->pull()),_links[i]->wire()->layerId());
    chi2+=_links[i]->pull();
    if (_links[i]->pull() > sigma){
      jbad=i; 
      sigma = _links[i]->pull();
    }
  }
        
  chi2=chi2/n;
  //if(sigma>9){bad.append(_links[jbad]);if(g_sigmalyr)g_sigmalyr->fill(sqrt(sigma),_links[jbad]->wire()->layerId());}
//  if(sigma>20.25&&chi2>2.08){bad.append(_links[jbad]);}
  if(sigma>9){bad.append(_links[jbad]);}
  if (bad.length()) {
    _links.remove(bad);
    _fitted = false;
    _updated = false;
  }
  return bad.length();
}

void TTrackBase::dump	(	const std::string &	message = std::string(""),
		const std::string &	prefix = std::string("")
	)			const [virtual, inherited]

dumps debug information.

Reimplemented in TCircle, TLine0, TMLine, TSegment, TSegment0, and TTrack.

Definition at line 58 of file TTrackBase.cxx.

References TTrackBase::_links, TTrackBase::_updated, Dump(), TTrackBase::mc(), RealDBUtil::npos, and TTrackBase::update().

Referenced by TTrack::dump(), TSegment0::dump(), TSegment::dump(), TMLine::dump(), TLine0::dump(), and TCircle::dump().

                                                                   {
    bool mc = (msg.find("mc") != std::string::npos);
    bool pull = (msg.find("pull") != std::string::npos);
    bool flag = (msg.find("flag") != std::string::npos);
    bool detail = (msg.find("detail") != std::string::npos);
    if (detail)
        mc = pull = flag = true;

    if (detail || (msg.find("layer") != std::string::npos)) {
        if (! _updated) update();
    }
    if (detail || (msg.find("hits") != std::string::npos)) {
        Dump(_links, msg, pre);
    }
}

double TRunge::dz

(

void

)

const

Definition at line 218 of file TRunge.cxx.

References _a.

Referenced by approach_line(), and intersect_xy_plane().

                           {
  return(_a[3]);
}

const SymMatrix & TRunge::Ea

(

const SymMatrix &

)

set helix error matrix

Definition at line 307 of file TRunge.cxx.

References _Ea, and _Nstep.

                                               {
  _Ea=tEa;
  _Nstep=0;
  return(_Ea);
}

const SymMatrix & TRunge::Ea

(

void

)

const

returns error matrix

Definition at line 234 of file TRunge.cxx.

References _Ea.

                                     {
  return(_Ea);
}

void TRunge::eloss	(	double	path,
		const RkFitMaterial *	material,
		double	mass,
		double	y[6],
		int	index
	)			const

Definition at line 1299 of file TRunge.cxx.

References RkFitMaterial::dE().

Referenced by RkFitCylinder::updateTrack().

                                                                                                      {
double psq=y[5]*y[5]+y[3]*y[3]+y[4]*y[4];

double p=sqrt(psq);

double dE=materiral->dE(mass,path,p);
if (index > 0)
  psq += dE * (dE + 2 * sqrt(mass * mass + psq));
else
  psq += dE * (dE - 2 * sqrt(mass * mass + psq));
double pnew=sqrt(psq);
double coeff=pnew/p;

y[3]=y[3]*coeff;
y[4]=y[4]*coeff;
y[5]=y[5]*coeff;
}

double TRunge::Eps

(

double

)

Definition at line 331 of file TRunge.cxx.

References _eps, and _Nstep.

                            {
  _eps=eps;
  _Nstep=0;
  return(_eps);
}

double TRunge::Eps

(

void

)

const

Definition at line 269 of file TRunge.cxx.

References _eps.

                            {
  return(_eps);
}

void TTrackBase::falseFit

(

)

[inline, inherited]

false Fit

Definition at line 229 of file TTrackBase.h.

References TTrackBase::_fitted, and TTrackBase::_fittedWithCathode.

                    {
    _fitted = false;
    _fittedWithCathode = false;
}

int TTrackBase::fit

(

void

)

[virtual, inherited]

fits itself by a default fitter. Error was happened if return value is not zero.

Definition at line 357 of file TTrackBase.cxx.

References TTrackBase::_fitter, and TMFitter::fit().

Referenced by TBuilder0::buildRphi(), TBuilder::buildRphi(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), and TConformalFinder::expand().

                    {
    return _fitter->fit(* this);
}

bool TTrackBase::fitted

(

void

)

const [inline, inherited]

returns true if fitted.

Definition at line 222 of file TTrackBase.h.

References TTrackBase::_fitted.

Referenced by TCosmicFitter::fit(), TBuilder::searchLine(), and TMSelector::select().

                             {
    return _fitted;
}

bool TTrackBase::fittedWithCathode

(

void

)

const [inline, inherited]

returns true if fitted with cathode hits(TEMPORARY).

Definition at line 243 of file TTrackBase.h.

References TTrackBase::_fittedWithCathode.

                                        {
    return _fittedWithCathode;
}

const TMFitter *const TTrackBase::fitter

(

const TMFitter *

)

[inline, inherited]

sets a default fitter.

Definition at line 261 of file TTrackBase.h.

References TTrackBase::_fitted, and TTrackBase::_fitter.

                                     {
    _fitted = false;
    return _fitter = a;
}

const TMFitter *const TTrackBase::fitter

(

void

)

const [inline, inherited]

returns a pointer to a default fitter.

Definition at line 255 of file TTrackBase.h.

References TTrackBase::_fitter.

Referenced by T3DLine::T3DLine(), TCircle::TCircle(), TLine0::TLine0(), TMLine::TMLine(), TRunge(), and TTrack::TTrack().

                             {
    return _fitter;
}

unsigned TRunge::Fly

(

const RkFitMaterial

material

)

const

make the trajectory in cache, return the number of step

Definition at line 836 of file TRunge.cxx.

References _h, _maxflightlength, _Nstep, _stepSize, _y, ganga-rec::j, Nstep(), Propagate(), SetFirst(), and TRunge_MAXstep.

Referenced by approach_line(), and approach_point().

                                                       {
  double y[6];
  unsigned Nstep;
  double flightlength;

  flightlength=0;
  SetFirst(y);
  for(Nstep=0;Nstep<TRunge_MAXstep;Nstep++){
    for(unsigned j=0;j<6;j++) _y[Nstep][j]=y[j];
    //memcpy(_y[Nstep],y,sizeof(double)*6);

    Propagate(y,_stepSize,material);

    if( y[2]>160 || y[2]<-80 || (y[0]*y[0]+y[1]*y[1])>8100 ) break;
    //if position is out side of MDC, stop to fly
    //  R>90cm, z<-80cm,160cm<z

    flightlength += _stepSize;
    if(flightlength>_maxflightlength) break;

    _h[Nstep]=_stepSize;
  }
  _Nstep=Nstep+1;
  return(_Nstep);
}

unsigned TRunge::Fly_SC

(

void

)

const

Definition at line 1097 of file TRunge.cxx.

References _eps, _h, _maxflightlength, _Nstep, _stepSizeMax, _stepSizeMin, _y, _yscal, default_stepSize0, Function(), genRecEmupikp::i, Nstep(), Propagate_QC(), SetFirst(), and TRunge_MAXstep.

Referenced by approach_line(), and approach_point().

                                 {//fly the particle track with stepsize control
  double y[6],dydx[6],yscal[6];
  unsigned Nstep;
  double step,stepdid,stepnext;
  unsigned i;
  double flightlength;

  //yscal[0]=0.1; //1mm  -> error = 1mm*EPS
  //yscal[1]=0.1; //1mm
  //yscal[2]=0.1; //1mm
  //yscal[3]=0.001; //1MeV
  //yscal[4]=0.001; //1MeV
  //yscal[5]=0.001; //1MeV

  step=default_stepSize0;
  flightlength=0;
  SetFirst(y);
  for(Nstep=0;Nstep<TRunge_MAXstep;Nstep++){
    for(i=0;i<6;i++) _y[Nstep][i]=y[i]; // save each step
    //memcpy(_y[Nstep],y,sizeof(double)*6);

    Function(y,dydx);
    //for(i=0;i<6;i++) yscal[i]=abs(y[i])+abs(dydx[i]*step)+TINY;

    Propagate_QC(y,dydx,step,_eps,_yscal,stepdid,stepnext);

    if( y[2]>160 || y[2]<-80 || (y[0]*y[0]+y[1]*y[1])>8100 ) break;
    //if position is out side of MDC, stop to fly
    //  R>90cm, z<-80cm,160cm<z

    flightlength += stepdid;
    if(flightlength>_maxflightlength) break;

    _h[Nstep]=stepdid;
    //cout<<"TR:"<<Nstep<<":step="<<stepdid<<endl;
    if(stepnext<_stepSizeMin)      step=_stepSizeMin;
    else if(stepnext>_stepSizeMax) step=_stepSizeMax;
    else step=stepnext;
  }
  _Nstep=Nstep+1;
  //cout<<"TR: Nstep="<<_Nstep<<endl;
  return(_Nstep);
}

void TRunge::Function	(	const double	y[6],
		double	f[6]
	)			const

Definition at line 916 of file TRunge.cxx.

References _charge, EvtCyclic3::B, IMagneticFieldSvc::fieldVector(), m_pmgnIMF, and boss::pos.

Referenced by approach_line(), Fly_SC(), Propagate(), Propagate1(), and Propagate_QC().

                                                        {
  // return the value of formula
  //  y[6] = (x,y,z,px,py,pz)
  //  f[6] = ( dx[3]/ds, dp[3]/ds )
  //  dx/ds = p/|p|
  //  dp/ds = e/|e| 1/alpha (p x B)/|p||B|
  //    alpha = 1/cB = 333.5640952/B[Tesla]  [cm/(GeV/c)]
  //  const float Bx = _bfield->bx(y[0],y[1],y[2]);
  //  const float By = _bfield->by(y[0],y[1],y[2]);
  //  const float Bz = _bfield->bz(y[0],y[1],y[2]);  //[kGauss]
  double pmag;
  double factor;
  HepVector3D B;
  HepPoint3D pos;
  pos.setX((float)y[0]);
  pos.setY((float)y[1]);
  pos.setZ((float)y[2]);
  //cout<<"TR::pos="<<pos[0]<<","<<pos[1]<<","<<pos[2]<<endl;
//  _bfield->fieldMap(pos,B);
  //cout<<"TR::B="<<B[0]<<","<<B[1]<<","<<B[2]<<endl;
  //  const double Bmag = sqrt(Bx*Bx+By*By+Bz*Bz);
  m_pmgnIMF->fieldVector(10*pos,B);
  pmag = sqrt(y[3]*y[3]+y[4]*y[4]+y[5]*y[5]);
  f[0] = y[3]/pmag;     // p/|p|
  f[1] = y[4]/pmag;
  f[2] = y[5]/pmag;

  //  const factor = _charge/(alpha2/Bmag)/pmag/Bmag;
 // factor = ((double)_charge)/alpha2/10.; //[(GeV/c)/(cm kG)]
 // factor = ((double)_charge)/(333.564095/(-1000*(m_pmgnIMF->getReferField())))/10.; //[(GeV/c)/(cm kG)]
  double Bmag=sqrt(B.x()*B.x()+B.y()*B.y()+B.z()*B.z());
  factor = ((double)_charge)/(0.333564095);
  //  f[3] = factor*(f[1]*Bz-f[2]*By);
  //  f[4] = factor*(f[2]*Bx-f[0]*Bz);
  //  f[5] = factor*(f[0]*By-f[1]*Bx);
  f[3] = factor*(f[1]*B.z()-f[2]*B.y());
  f[4] = factor*(f[2]*B.x()-f[0]*B.z());
  f[5] = factor*(f[0]*B.y()-f[1]*B.x());
}

int TRunge::GetStep	(	unsigned	stepNum,
		double &	step
	)			const

Definition at line 985 of file TRunge.cxx.

References _h, _Nstep, and TRunge_MAXstep.

                                                      {
  if(stepNum>=_Nstep || stepNum>=TRunge_MAXstep) return(-1);
  step=_h[stepNum];
  return(0);
}

int TRunge::GetXP	(	unsigned	stepNum,
		double	y[6]
	)			const

Definition at line 979 of file TRunge.cxx.

References _Nstep, _y, genRecEmupikp::i, and TRunge_MAXstep.

                                                   {
  if(stepNum>=_Nstep || stepNum>=TRunge_MAXstep) return(-1);

  for(unsigned i=0;i<6;i++) y[i]=_y[stepNum][i];
  return(0);
}

Helix TRunge::helix

(

void

)

const

returns helix class

Definition at line 238 of file TRunge.cxx.

References _a, _Ea, _pivot, and Helix.

Referenced by RkFitCylinder::intersect(), intersect_cylinder(), intersect_xy_plane(), intersect_yz_plane(), intersect_zx_plane(), and pivot().

                             {
  return(Helix(_pivot,_a,_Ea));
}

const TTrackHEP *const TTrackBase::hep

(

void

)

const [inherited]

returns TTrackHEP.

Definition at line 407 of file TTrackBase.cxx.

References TTrackBase::_hep, TTrackBase::_links, TTrackBase::_nHeps, genRecEmupikp::i, ganga-rec::j, max, and TTrackBase::mc().

Referenced by TTrackBase::nHeps().

                          {
    unsigned n = _links.length();
    CAList<TTrackHEP> hepList;
    CList<unsigned> hepCounter;
    for (unsigned i = 0; i < n; i++) {
        const TTrackHEP * hep = _links[i]->hit()->mc()->hep();
        unsigned nH = hepList.length();
        bool found = false;
        for (unsigned j = 0; j < nH; j++) {
            if (hepList[j] == hep) {
                found = true;
                ++(* hepCounter[j]);
            }
        }

        if (! found) {
            hepList.append(hep);
            unsigned c = 0;
            hepCounter.append(c);
        }
    }

    _nHeps = hepList.length();
    _hep = 0;
    unsigned max = 0;
    for (unsigned i = 0; i < _nHeps; i++) {
        if ((* hepCounter[i]) > max) {
            max = (* hepCounter[i]);
            _hep = hepList[i];
        }
    }

    return _hep;
}

double TRunge::intersect_cylinder

(

double

r

)

const

Definition at line 1316 of file TRunge.cxx.

References Helix::center(), helix(), M_PI, Helix::phi0(), and Helix::radius().

Referenced by RkFitCylinder::intersect(), intersect_yz_plane(), and intersect_zx_plane().

                                                {
    double m_rad = helix().radius();
    double l = helix().center().perp();
    double cosPhi = (m_rad * m_rad + l * l  - r * r) / (2 * m_rad * l);

    if(cosPhi < -1 || cosPhi > 1) return 0.;

    double dPhi = helix().center().phi() - acos(cosPhi) - helix().phi0();
    if(dPhi < -M_PI) dPhi += 2 * M_PI;

    return dPhi;
}

double TRunge::intersect_xy_plane

(

double

z

)

const

Definition at line 1362 of file TRunge.cxx.

References dz(), helix(), pivot(), Helix::radius(), Helix::tanl(), and tanl().

Referenced by RkFitCylinder::intersect().

                                                {
    if (helix().tanl() != 0 && helix().radius() != 0)
      return (helix().pivot().z() + helix().dz() - z) / (helix().radius() * helix().tanl());
    else return 0;
}

double TRunge::intersect_yz_plane	(	const HepTransform3D &	plane,
		double	x
	)			const

Definition at line 1345 of file TRunge.cxx.

References Helix::center(), helix(), intersect_cylinder(), and Helix::radius().

                                                                             {
    HepPoint3D xc = plane * helix().center();
    double r = helix().radius();
    double d = r * r - (x - xc.x()) * (x - xc.x());
    if(d < 0) return 0;

    double yy = xc.y();
    if(yy > 0) yy -= sqrt(d);
    else yy += sqrt(d);

    double l = (plane.inverse() *
                HepPoint3D(x, yy, 0)).perp();

    return intersect_cylinder(l);
}

double TRunge::intersect_zx_plane	(	const HepTransform3D &	plane,
		double	y
	)			const

Definition at line 1328 of file TRunge.cxx.

References Helix::center(), helix(), intersect_cylinder(), and Helix::radius().

                                                                             {
    HepPoint3D xc = plane * helix().center();
    double r = helix().radius();
    double d = r * r - (y - xc.y()) * (y - xc.y());
    if(d < 0) return 0;

    double xx = xc.x();
    if(xx > 0) xx -= sqrt(d);
    else xx += sqrt(d);

    double l = (plane.inverse() *
                HepPoint3D(xx, y, 0)).perp();

    return intersect_cylinder(l);
}

double TRunge::kappa

(

void

)

const

Definition at line 214 of file TRunge.cxx.

References _a.

                              {
  return(_a[2]);
}

const AList< TMLink > & TTrackBase::links

(

unsigned

mask = 0

)

const [inherited]

returns a list of masked TMLinks assigned to this track. 'mask' will be applied if mask is not 0.

Definition at line 297 of file TTrackBase.cxx.

References TTrackBase::_links.

Referenced by TConformalFinder0::appendCluster(), TBuilder::build(), TBuilderCurl::buildStereo(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo0(), TBuilderCurl::buildStereoMC(), TBuilder::buildStereoNew(), TBuilderCurl::check(), TCurlFinder::check2DCircle(), TCurlFinder::check3DTrack(), TCurlFinder::dividing2DTrack(), TCurlFinder::dumpType1(), TCurlFinder::dumpType2(), TrkReco::execute(), TConformalFinder0::findCloseClusters(), TFastFinder::findCloseHits(), TConformalFinder0::findCloseHits(), TConformalFinder0::findClusterLink(), TCosmicFitter::fit(), TCurlFinder::fitWDD(), TCurlFinder::make2DTrack(), TCurlFinder::make3DTrack(), TCurlFinder::makeWithMC(), TCurlFinder::mask3DTrack(), TTrackManager::merge(), TTrack::movePivot(), TCurlFinder::plotCircle(), TCurlFinder::plotTrack(), TMSelector::preSelect(), TCurlFinder::removeStereo(), TConformalFinder::salvage(), TCurlFinder::salvage3DTrack(), TBuilder::searchLine(), TMDCTsf::segments(), TBuilderCurl::setArcZ(), TConformalFinder::slowFinding2D(), TSegment0::splitAV(), TSegment::splitAV(), TSegment::splitTsf(), TConformalFinder::stereoSegments(), TSegment0::superLayerId(), TSegment::superLayerId(), TTrack::szPosition(), TCurlFinder::trace2DTrack(), TCurlFinder::trace3DTrack(), and TSegment::width().

                                     {
    if (mask == 0) return _links;

    std::cout << "TTrackBase::links !!! mask is not supportted yet" << std::endl;
    return _links;
}

float TRunge::Mass

(

float

)

set mass for tof calc.

Definition at line 347 of file TRunge.cxx.

References _mass, and _mass2.

                            {
  _mass=mass;
  _mass2=_mass*_mass;
  return(_mass);
}

float TRunge::Mass

(

void

)

const

return mass

Definition at line 279 of file TRunge.cxx.

References _mass.

                            {
  return(_mass);
}

double TRunge::MaxFlightLength

(

double

)

Definition at line 353 of file TRunge.cxx.

References _maxflightlength, _Nstep, and default_maxflightlength.

                                           {
  if(length>0) _maxflightlength=length;
  else _maxflightlength=default_maxflightlength;
  _Nstep=0;
  return(_maxflightlength);
}

double TRunge::MaxFlightLength

(

void

)

const

return max flight length

Definition at line 283 of file TRunge.cxx.

References _maxflightlength.

                                        {
  return(_maxflightlength);
}

const TTrackMC *const TTrackBase::mc

(

void

)

const [inline, inherited]

returns a pointer to TTrackMC.

Definition at line 249 of file TTrackBase.h.

References TTrackBase::_mc.

Referenced by TTrackBase::dump(), and TTrackBase::hep().

                         {
    return _mc;
}

unsigned TTrackBase::nCores

(

unsigned

mask = 0

)

const [inherited]

returns # of masked TMLinks for fit. 'mask' will be applied if mask is not 0.

Definition at line 325 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_updated, and TTrackBase::update().

Referenced by TTrack::dump(), and TConformalFinder::expand().

                                      {
    if (mask)
        std::cout << "TTrackBase::nCores !!! mask is not supported" << std::endl;
    if (! _updated) update();
    return _cores.length();
}

unsigned TRunge::ndf

(

void

)

const

returns NDF

Definition at line 242 of file TRunge.cxx.

References _ndf.

                              {
  return _ndf;
}

unsigned TTrackBase::nHeps

(

void

)

const [inherited]

returns # of contributed TTrackHEP tracks.

Definition at line 443 of file TTrackBase.cxx.

References TTrackBase::_nHeps, and TTrackBase::hep().

                            {
    hep();
    return _nHeps;
}

unsigned TTrackBase::nLinks

(

unsigned

mask = 0

)

const [inherited]

returns # of masked TMLinks assigned to this track object.

Definition at line 305 of file TTrackBase.cxx.

References TTrackBase::_links, and genRecEmupikp::i.

Referenced by TBuilderCurl::buildStereo(), TBuilderCurl::buildStereoMC(), TBuilderCurl::check(), TCurlFinder::check2DCircle(), TCurlFinder::check3DTrack(), TSegment0::clusterType(), TSegment::clusterType(), TCurlFinder::dividing2DTrack(), TCurlFinder::dumpType1(), TCurlFinder::dumpType2(), TCurlFinder::make2DTrack(), TCurlFinder::makeWithMC(), TCurlFinder::merge3DTrack(), TCurlFinder::plotCircle(), TCurlFinder::plotTrack(), TMSelector::preSelect(), TCurlFinder::salvage3DTrack(), TBuilderCurl::setArcZ(), TSegment0::splitComplicated(), TSegment::splitComplicated(), TSegment0::splitParallel(), TSegment::splitParallel(), TSegment0::updateType(), and TSegment::updateType().

                                      {
    unsigned n = _links.length();
    if (mask == 0) return n;
    unsigned nn = 0;
    for (unsigned i = 0; i < n; i++) {
        const TMDCWireHit & h = * _links[i]->hit();
        if (h.state() & mask) ++nn;
    }
    return nn;
}

unsigned TRunge::Nstep

(

void

)

const

access to trajectory cache

Definition at line 975 of file TRunge.cxx.

References _Nstep.

Referenced by Fly(), and Fly_SC().

                                {
  return(_Nstep);
}

unsigned TRunge::objectType

(

void

)

const [inline, virtual]

returns object type

Reimplemented from TTrackBase.

Definition at line 256 of file TRunge.h.

References Runge.

                                            {
  return Runge;
}

TMLink * TTrackBase::operator[]

(

unsigned

i

)

const [inline, inherited]

Definition at line 237 of file TTrackBase.h.

References TTrackBase::_links.

                                       {
    return _links[i];
}

double TRunge::phi0

(

void

)

const

Definition at line 210 of file TRunge.cxx.

References _a.

                             {
  return(_a[1]);
}

const HepPoint3D & TRunge::pivot

(

const HepPoint3D &

)

set new pivot

!!!!! under construction !!!!! track parameter should be extracted after track propagation.

Definition at line 287 of file TRunge.cxx.

References _a, _Ea, _Nstep, _pivot, Helix::a(), Helix::Ea(), helix(), and Helix::pivot().

                                                         {
  Helix tHelix(helix());
  tHelix.pivot(newpivot);
  _pivot=newpivot;
  _a=tHelix.a();
  _Ea=tHelix.Ea();
  _Nstep=0;
  return(_pivot);
}

const HepPoint3D & TRunge::pivot

(

void

)

const

pivot position

Definition at line 226 of file TRunge.cxx.

References _pivot.

Referenced by intersect_xy_plane().

                                         {
  return(_pivot);
}

void TRunge::Propagate	(	double	y[6],
		const double &	step,
		const RkFitMaterial	material
	)			const

propagate the track using 4th-order Runge-Kutta method

Definition at line 862 of file TRunge.cxx.

References f1, f2, Function(), genRecEmupikp::i, me, and mpi.

Referenced by approach_line(), approach_point(), and Fly().

                                                                                        {
  //  y[6] = (x,y,z,px,py,pz)
  double f1[6],f2[6],f3[6],f4[6],yt[6];
  double hh;
  double h6;
  unsigned i;
  const double mpi=0.13957;
  double me=0.000511;
  //eloss(step,&material, me,y,0);
  hh = step*0.5;
  //cout<<"TR:Pro hh="<<hh<<endl;
  Function(y,f1);
  for(i=0;i<6;i++) yt[i]=y[i]+hh*f1[i];
  //{
  //  register double *a=y;
  //  register double *b=f1;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function(yt,f2);
  for(i=0;i<6;i++) yt[i]=y[i]+hh*f2[i];
  //{
  //  register double *a=y;
  //  register double *b=f2;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function(yt,f3);
  for(i=0;i<6;i++) yt[i]=y[i]+step*f3[i];
  //{
  //  register double *a=y;
  //  register double *b=f3;
  //  register double *t=yt;
  //  register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+step*(*b);
  //}
  Function(yt,f4);

  h6 = step/6;
  for(i=0;i<6;i++) y[i]+=h6*(f1[i]+2*f2[i]+2*f3[i]+f4[i]);
  //{
  //  register double *a=f1;
  //  register double *b=f2;
  //  register double *c=f3;
  //  register double *d=f4;
  //  register double *t=y;
  //  register double *e=y+6;
  //  for(;t<e;a++,b++,c++,d++,t++)
  //    (*t)+=h6*((*a)+2*(*b)+2*(*c)+(*d));
  //}
}

void TRunge::Propagate1	(	const double	y[6],
		const double	dydx[6],
		const double &	step,
		double	yout[6]
	)			const

Definition at line 991 of file TRunge.cxx.

References f2, Function(), and genRecEmupikp::i.

Referenced by Propagate_QC().

                                                                 {
  //  y[6] = (x,y,z,px,py,pz)
  double f2[6],f3[6],f4[6],yt[6];
  double hh;
  double h6;
  unsigned i;

  hh = step*0.5;
  for(i=0;i<6;i++) yt[i]=y[i]+hh*dydx[i];       // 1st step
  //{
  //  const register double *a=y;
  //  const register double *b=dydx;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function(yt,f2);                              // 2nd step
  for(i=0;i<6;i++) yt[i]=y[i]+hh*f2[i];
  //{
  //  const register double *a=y;
  //  const register double *b=f2;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+hh*(*b);
  //}
  Function(yt,f3);                              // 3rd step
  for(i=0;i<6;i++) yt[i]=y[i]+step*f3[i];
  //{
  //  const register double *a=y;
  //  const register double *b=f3;
  //  register double *t=yt;
  //  const register double *e=y+6;
  //  for(;a<e;a++,b++,t++) (*t)=(*a)+step*(*b);
  //}
  Function(yt,f4);                              // 4th step

  h6 = step/6;
  for(i=0;i<6;i++) yout[i]=y[i]+h6*(dydx[i]+2*f2[i]+2*f3[i]+f4[i]);
  //{
  //  const register double *a=dydx;
  //  const register double *b=f2;
  //  const register double *c=f3;
  //  const register double *d=f4;
  //  const register double *e=y;
  //  register double *t=yout;
  //  const register double *s=yout+6;
  //  for(;t<s;a++,b++,c++,d++,e++,t++)
  //    (*t)=(*e)+h6*((*a)+2*(*b)+2*(*c)+(*d));
  //}
}

void TRunge::Propagate_QC	(	double	y[6],
		double	dydx[6],
		const double &	steptry,
		const double &	eps,
		const double	yscal[6],
		double &	stepdid,
		double &	stepnext
	)			const

Definition at line 1048 of file TRunge.cxx.

References ERRCON, exp(), FCOR, Function(), genRecEmupikp::i, PGROW, Propagate1(), PSHRNK, SAFETY, and subSeperate::temp.

Referenced by Fly_SC().

                                                                  {
  //propagate with quality check, step will be scaled automatically
  double ysav[6],dysav[6],ytemp[6];
  double h,hh,errmax,temp;
  unsigned i;

  for(i=0;i<6;i++) ysav[i]=y[i];
  //memcpy(ysav,y,sizeof(double)*6);
  for(i=0;i<6;i++) dysav[i]=dydx[i];
  //memcpy(dysav,dydx,sizeof(double)*6);

  h=steptry;
  for(;;){
    hh=h*0.5;                           // 2 half step
    Propagate1(ysav,dysav,hh,ytemp);
    Function(ytemp,dydx);
    Propagate1(ytemp,dydx,hh,y);
    //if  check step size
    Propagate1(ysav,dysav,h,ytemp);     // 1 full step
    // error calculation
    errmax=0.0;
    for(i=0;i<6;i++){
      ytemp[i]=y[i]-ytemp[i];
      temp=fabs(ytemp[i]/yscal[i]);
      if(errmax < temp) errmax=temp;
    }
    //cout<<"TR: errmax="<<errmax<<endl;
    errmax/=eps;
    if(errmax <= 1.0){                  // step O.K.  calc. for next step
      stepdid=h;
      stepnext=(errmax>ERRCON ? SAFETY*h*exp(PGROW*log(errmax)) : 4.0*h);
      break;
    }
    h=SAFETY*h*exp(PSHRNK*log(errmax));
  }
  for(i=0;i<6;i++) y[i]+=ytemp[i]*FCOR;
  //{
  //  register double *a=ytemp;
  //  register double *t=y;
  //  register double *e=y+6;
  //  for(;t<e;a++,t++) (*t)+=(*a)*FCOR;
  //}

}

double TRunge::reducedchi2

(

void

)

const

returns reduced-chi2

Definition at line 250 of file TRunge.cxx.

References _chi2, and _ndf.

                                    {
  if(_ndf==0){
    std::cout<<"error at TRunge::reducedchi2  ndf=0"<<std::endl;
    return 0;
  }
  return (_chi2/_ndf);
}

void TTrackBase::refine

(

double

maxSigma

)

[virtual, inherited]

removes bad points by pull. The bad points are masked not to be used in fit.

Definition at line 194 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_updated, Dump(), and TTrackBase::refineMain().

                               {
    AList<TMLink> bad = refineMain(sigma);
//      for (unsigned i = 0; i < bad.length(); i++) {
//      const TMDCWireHit * hit = bad[i]->hit();
//      hit->state(hit->state() | WireHitInvalidForFit);
//      }

#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    refine ... sigma=" << sigma << std::endl;
    Dump(bad, "detail sort", "        ");
#endif

    if (bad.length()) {
        _fitted = false;
        _updated = false;
    }
}

void TTrackBase::refine	(	AList< TMLink > &	list,
		double	maxSigma
	)			[virtual, inherited]

removes bad points by pull. The bad points are removed from the track, and are returned in 'list'.

Definition at line 170 of file TTrackBase.cxx.

References TTrackBase::_fitted, TTrackBase::_links, TTrackBase::_updated, Dump(), and TTrackBase::refineMain().

Referenced by TBuilder0::appendClusters(), TBuilderCurl::buildStereo(), TBuilderCosmic::buildStereo(), TBuilder0::buildStereo(), TBuilder::buildStereo(), TBuilder0::buildStereo0(), and TBuilderCurl::buildStereoMC().

                                                     {
    AList<TMLink> bad = refineMain(sigma);
#ifdef TRKRECO_DEBUG
    std::cout << "    refine ... sigma=" << sigma << ", # of rejected hits=";
    std::cout << bad.length() << std::endl;
#endif
#ifdef TRKRECO_DEBUG
    Dump(bad, "sort pull mc", "        ");
#endif

    if (bad.length()) {
        _links.remove(bad);
        list.append(bad);
        _fitted = false;
        _updated = false;
    }
}

void TTrackBase::remove

(

const AList< TMLink > &

)

[inline, inherited]

removes TMLinks.

Definition at line 213 of file TTrackBase.h.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, and TTrackBase::_updated.

                                          {
    _links.remove(a);
    _updated = false;
    _fitted = false;
    _fittedWithCathode = false; // mod. by matsu ( 1999/05/24 )
}

void TTrackBase::remove

(

TMLink &

a

)

[inline, inherited]

removes a TMLink.

Definition at line 204 of file TTrackBase.h.

References TTrackBase::_fitted, TTrackBase::_fittedWithCathode, TTrackBase::_links, and TTrackBase::_updated.

Referenced by TConformalFinder0::appendClusters2(), TBuilder::buildStereoNew(), TCurlFinder::dividing2DTrack(), TCurlFinder::makeWithMC(), TTrackManager::maskCurl(), TTrackManager::maskNormal(), TCurlFinder::removeStereo(), TBuilder::searchLine(), TConformalFinder0::specialFinding(), TSegment0::splitComplicated(), and TSegment::splitComplicated().

                             {
    _links.remove(a);
    _updated = false;
    _fitted = false;
    _fittedWithCathode = false; // mod. by matsu ( 1999/05/24 )
}

void TTrackBase::removeLinks

(

void

)

[virtual, inherited]

Definition at line 189 of file TTrackBase.cxx.

References TTrackBase::_links.

Referenced by TBuilderCosmic::buildStereo(), and TrkReco::execute().

                           {
  _links.removeAll();
}

void TRunge::SetFirst

(

double

y[6]

)

const

set first point (position, momentum) s=0, phi=0

Definition at line 956 of file TRunge.cxx.

References _a, _pivot, abs, cos(), and sin().

Referenced by Fly(), and Fly_SC().

                                      {
  //  y[6] = (x,y,z,px,py,pz)
  const double cosPhi0=cos(_a[1]);
  const double sinPhi0=sin(_a[1]);
  const double invKappa=1/abs(_a[2]);

  y[0]= _pivot.x()+_a[0]*cosPhi0;       //x
  y[1]= _pivot.y()+_a[0]*sinPhi0;       //y
  y[2]= _pivot.z()+_a[3];               //z
  y[3]= -sinPhi0*invKappa;              //px
  y[4]=  cosPhi0*invKappa;              //py
  y[5]= _a[4]*invKappa;                 //pz
}

double TRunge::SetFlightLength

(

void

)

set flight length using wire hits

Definition at line 1141 of file TRunge.cxx.

References _maxflightlength, abs, TTrackBase::cores(), DBL_MAX, Helix::dz(), IMagneticFieldSvc::getReferField(), genRecEmupikp::i, ganga-rec::j, m_pmgnIMF, Helix::pivot(), tanl(), TMDCWire::wirePosition(), and Helix::x().

                                  {

  // get a list of links to a wire hit
  const AList<TMLink>& cores=this->cores();
  //cout<<"TR:: cores.length="<<cores.length()<<endl;

  const Helix hel(this->helix());
  double tanl=hel.tanl();
  //double curv=hel.curv();
  //double rho = Helix::ConstantAlpha / hel.kappa();
 // double rho = 333.564095 / hel.kappa();  
  const double Bz = 1000*m_pmgnIMF->getReferField();
  double rho = 333.564095/(hel.kappa()*Bz);  

  // search max phi
  double phi_max=- DBL_MAX;
  double phi_min= DBL_MAX;
  // loop over link
  for(int j=0;j<cores.length();j++){
    TMLink& l=*cores[j];
    // fitting valid?
    const TMDCWire& wire=*l.wire();
    double Wz=0;
    //double mindist;
    HepPoint3D onWire;
    HepPoint3D dummy_bP;
    HepVector3D dummy_dir;
    Helix th(hel);
    for(int i=0;i<10;i++){
      wire.wirePosition(Wz,onWire,dummy_bP,dummy_dir);
      th.pivot(onWire);
      if(abs(th.dz())<0.1){     //<1mm
        //mindist=abs(hel.dr());
        break;
      }
      Wz+=th.dz();
    }
    //onWire is closest point , th.x() is closest point on trk
    //Wz is z position of closest point
    //Z->dphi
    /*
    // Calculate position (x,y,z) along helix.
    // x = x0 + dr * cos(phi0) + (alpha / kappa) * (cos(phi0) - cos(phi0+phi))
    // y = y0 + dr * sin(phi0) + (alpha / kappa) * (sin(phi0) - sin(phi0+phi))
    // z = z0 + dz             - (alpha / kappa) * tan(lambda) * phi
    cout<<"TR:: Wz="<<Wz<<" tanl="<<tanl<<endl;
    double phi=( hel.pivot().z()+hel.dz()-Wz )/( curv*tanl );
    */
    //...Cal. dPhi to rotate...
    const HepPoint3D & xc = hel.center();
    const HepPoint3D & xt = hel.x();
    HepVector3D v0, v1;
    v0 = xt - xc;
    v1 = th.x() - xc;
    const double vCrs = v0.x() * v1.y() - v0.y() * v1.x();
    const double vDot = v0.x() * v1.x() + v0.y() * v1.y();
    double phi = atan2(vCrs, vDot);

    double tz=hel.x(phi).z();
    //cout<<"TR::  Wz="<<Wz<<" tz="<<tz<<endl;
    
    if(phi>phi_max) phi_max=phi;
    if(phi<phi_min) phi_min=phi;
    //cout<<"TR:: phi="<<phi<<endl;
  }
  //cout<<"TR:: phi_max="<<phi_max<<" phi_min="<<phi_min
  //    <<" rho="<<rho<<" tanl="<<tanl<<endl;
  //phi->length, set max filght length
  //tanl*=1.2; // for mergin
  _maxflightlength=
    abs( rho*(phi_max-phi_min)*sqrt(1+tanl*tanl) )*1.2; //x1.1 mergin

  return(_maxflightlength);
}

double TRunge::StepSize

(

double

)

set step size to calc. trajectory

Definition at line 320 of file TRunge.cxx.

References _Nstep, and _stepSize.

                                  {
  _stepSize=step;
  _Nstep=0;
  return(_stepSize);
}

double TRunge::StepSize

(

void

)

const

returns step size

Definition at line 262 of file TRunge.cxx.

References _stepSize.

                                 {
  return(_stepSize);
}

double TRunge::StepSizeMax

(

double

)

Definition at line 336 of file TRunge.cxx.

References _Nstep, and _stepSizeMax.

                                     {
  _stepSizeMax=step;
  _Nstep=0;
  return(_stepSizeMax);
}

double TRunge::StepSizeMax

(

void

)

const

Definition at line 272 of file TRunge.cxx.

References _stepSizeMax.

                                    {
  return(_stepSizeMax);
}

double TRunge::StepSizeMin

(

double

)

Definition at line 341 of file TRunge.cxx.

References _Nstep, and _stepSizeMin.

                                     {
  _stepSizeMin=step;
  _Nstep=0;
  return(_stepSizeMin);
}

double TRunge::StepSizeMin

(

void

)

const

Definition at line 275 of file TRunge.cxx.

References _stepSizeMin.

                                    {
  return(_stepSizeMin);
}

double TRunge::tanl

(

void

)

const

Definition at line 222 of file TRunge.cxx.

References _a.

Referenced by intersect_xy_plane(), and SetFlightLength().

                             {
  return(_a[4]);
}

unsigned TTrackBase::testByApproach	(	const AList< TMLink > &	list,
		double	sigma
	)			const [inherited]

Definition at line 240 of file TTrackBase.cxx.

References genRecEmupikp::i, and TTrackBase::testByApproach().

                                                                            {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    TTrackBase::testByApproach ... sigma=" << maxSigma << std::endl;
#endif

    unsigned nOK = 0;
    unsigned n = list.length();
    for (unsigned i = 0; i < n; i++) {
        TMLink & l = * list[i];
        nOK += testByApproach(l, maxSigma);
    }
    return nOK;
}

unsigned TTrackBase::testByApproach	(	const TMLink &	list,
		double	sigma
	)			const [inherited]

returns # of good hits to be appended.

Definition at line 255 of file TTrackBase.cxx.

References TTrackBase::approach(), TMDCWireHit::dDrift(), TMLink::distance(), TTrackBase::distance(), TMDCWireHit::drift(), TMLink::dump(), showlog::err, TMLink::hit(), and TMLink::pull().

Referenced by TTrackBase::testByApproach().

                                                                     {
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    TTrackBase::testByApproach ... sigma=" << maxSigma << std::endl;
#endif
    TMLink l = link;

    //...Calculate closest approach...
    int err = approach(l);
    if (err < 0) return 0;
    //...Calculate sigma...
    float distance = l.distance();
    float diff = fabs(distance - l.hit()->drift());
    float sigma = diff / l.hit()->dDrift();
    l.pull(sigma * sigma);

    //...For debug...
#ifdef TRKRECO_DEBUG_DETAIL
    std::cout << "    sigma=" << sigma;
    std::cout << ",dist=" << distance;
    std::cout << ",diff=" << diff << ",";
    if (sigma < maxSigma) std::cout << "ok,";
    else                  std::cout << "X,";
    l.dump("mc");
#endif
                
    //...Make sigma cut...
    if (sigma < maxSigma) return 1;

    return 0;
}

unsigned TTrackBase::type

(

void

)

const [inline, virtual, inherited]

returns type. Definition is depending on an object class.

Reimplemented in TTrack.

Definition at line 274 of file TTrackBase.h.

                           {
    return 0;
}

void TTrackBase::update

(

void

)

const [inherited]

update cache.

Reimplemented in TSegment, and TSegment0.

Definition at line 75 of file TTrackBase.cxx.

References TTrackBase::_cores, TTrackBase::_links, TTrackBase::_updated, genRecEmupikp::i, TMDCWireHit::state(), WireHitFittingValid, and WireHitInvalidForFit.

Referenced by TTrackBase::cores(), TTrackBase::dump(), and TTrackBase::nCores().

                             {
    _cores.removeAll();
    unsigned n = _links.length();
    for (unsigned i = 0; i < n; i++) {
        TMLink * l = _links[i];
        const TMDCWireHit & h = * l->hit();
        if (h.state() & WireHitInvalidForFit) continue;
        if (! (h.state() & WireHitFittingValid)) continue;
        _cores.append(l);
    }
    _updated = true;
}

const double* TRunge::Yscal

(

const double *

)

set error parameters for fitting with step size control

const double * TRunge::Yscal

(

void

)

const

return error parameters for fitting with step size control

Definition at line 266 of file TRunge.cxx.

References _yscal.

                                     {
  return(_yscal);
}

---

Friends And Related Function Documentation

friend class TRungeFitter [friend]

Definition at line 224 of file TRunge.h.

---

Member Data Documentation

Vector TRunge::_a [private]

Definition at line 216 of file TRunge.h.

Referenced by a(), approach_line(), dr(), dz(), helix(), kappa(), phi0(), pivot(), SetFirst(), tanl(), and TRunge().

Bfield* TRunge::_bfield [private]

Definition at line 226 of file TRunge.h.

Referenced by BfieldID(), and TRunge().

int TRunge::_bfieldID [private]

Definition at line 227 of file TRunge.h.

Referenced by BfieldID(), and TRunge().

int TRunge::_charge [private]

Definition at line 245 of file TRunge.h.

Referenced by a(), Function(), and TRunge().

double TRunge::_chi2 [private]

Definition at line 219 of file TRunge.h.

Referenced by chi2(), and reducedchi2().

SymMatrix TRunge::_Ea [private]

Definition at line 217 of file TRunge.h.

Referenced by Ea(), helix(), and pivot().

double TRunge::_eps [private]

Definition at line 234 of file TRunge.h.

Referenced by Eps(), and Fly_SC().

bool TTrackBase::_fitted [mutable, protected, inherited]

Definition at line 162 of file TTrackBase.h.

Referenced by TMLine::a(), TLine0::a(), TTrackBase::append(), TTrackBase::appendByApproach(), TMLine::appendByszdistance(), TLine0::appendByszdistance(), TMLine::appendSLY(), TLine0::appendSLY(), TMLine::b(), TLine0::b(), TCircle::center(), TCircle::charge(), TTrack::chi2(), TMLine::chi2(), TLine0::chi2(), TSegment0::direction(), TSegment::direction(), TMLine::distance(), TLine0::distance(), TTrackBase::DropWorst(), TSegment0::dump(), TSegment::dump(), TCircle::dump(), TTrackBase::falseFit(), THelixFitter::fit(), TCosmicFitter::fit(), TMLine::fit2(), TLine0::fit2(), TMLine::fit2p(), TLine0::fit2p(), TMLine::fit2s(), TLine0::fit2s(), TMLine::fit2sp(), TLine0::fit2sp(), TCircle::fitForCurl(), TTrackBase::fitted(), TTrackBase::fitter(), TCircle::impact(), TSegment0::innerMostLayer(), TSegment::innerMostLayer(), TSegment::inners(), TSegment0::innerWidth(), TSegment::innerWidth(), TTrack::ndf(), TSegment0::outerMostLayer(), TSegment::outerMostLayer(), TSegment::outers(), TSegment0::outerWidth(), TSegment::outerWidth(), TSegment0::position(), TSegment::position(), TCircle::pt(), TCircle::radius(), TMLine::reducedChi2(), TLine0::reducedChi2(), TTrackBase::refine(), TMLine::refine(), TLine0::refine(), TTrack::refine2D(), TTrackBase::remove(), TMLine::removeChits(), TLine0::removeChits(), TMLine::removeSLY(), TLine0::removeSLY(), T3DLine::T3DLine(), TRunge(), TSegment::TSegment(), TSegment0::TSegment0(), TTrack::TTrack(), TSegment0::update(), TSegment::update(), TSegment0::updateType(), and TSegment::updateType().

bool TTrackBase::_fittedWithCathode [mutable, protected, inherited]

Definition at line 163 of file TTrackBase.h.

Referenced by TTrackBase::append(), TTrackBase::falseFit(), TTrackBase::fittedWithCathode(), TTrackBase::remove(), T3DLine::T3DLine(), TRunge(), and TTrack::TTrack().

const TRungeFitter TRunge::_fitter = TRungeFitter("TRunge Default fitter") [static, private]

Reimplemented from TTrackBase.

Definition at line 222 of file TRunge.h.

Referenced by TRunge().

double TRunge::_h[TRunge_MAXstep] [mutable, private]

Definition at line 248 of file TRunge.h.

Referenced by approach_line(), Fly(), Fly_SC(), and GetStep().

AList<TMLink> TTrackBase::_links [protected, inherited]

Definition at line 161 of file TTrackBase.h.

Referenced by TTrackBase::append(), TTrackBase::appendByApproach(), TMLine::appendByszdistance(), TLine0::appendByszdistance(), TMLine::appendSLY(), TLine0::appendSLY(), TTrack::assign(), TBuilderCurl::buildStereo(), TMLine::chi2(), TLine0::chi2(), TTrackBase::DropWorst(), TTrackBase::dump(), TTrack::dump(), TSegment0::dump(), TSegment::dump(), TMLine::dump(), TLine0::dump(), TCircle::dump(), TSegment::expandSeg(), TMLine::fit2(), TLine0::fit2(), TTrack::fit2D(), TMLine::fit2p(), TLine0::fit2p(), TMLine::fit2s(), TLine0::fit2s(), TMLine::fit2sp(), TLine0::fit2sp(), TCircle::fitForCurl(), TTrackBase::hep(), TTrackBase::links(), TTrack::movePivot(), TTrackBase::nLinks(), TTrackBase::operator[](), TSegment0::rangeX(), TSegment::rangeX(), TMLine::reducedChi2(), TLine0::reducedChi2(), TTrackBase::refine(), TMLine::refine(), TLine0::refine(), TTrack::refine2D(), TTrackBase::refineMain(), TTrackBase::remove(), TMLine::removeChits(), TLine0::removeChits(), TTrackBase::removeLinks(), TMLine::removeSLY(), TLine0::removeSLY(), TSegment::segSalvage(), TSegment0::solveDualHits(), TSegment::solveDualHits(), TSegment::solveLR(), TSegment::solveThreeHits(), TSegment0::splitAV(), TSegment::splitAV(), TSegment0::splitComplicated(), TSegment::splitComplicated(), TSegment0::splitDual(), TSegment::splitDual(), TSegment0::splitParallel(), TSegment::splitParallel(), TSegment::splitTsf(), TSegment::TSegment(), TSegment0::TSegment0(), TTrack::TTrack(), TTrackBase::update(), TSegment0::update(), TSegment::update(), TSegment0::updateDuality(), TSegment::updateDuality(), TSegment0::updateType(), and TSegment::updateType().

float TRunge::_mass [private]

Definition at line 238 of file TRunge.h.

Referenced by Mass(), and TRunge().

float TRunge::_mass2 [private]

Definition at line 244 of file TRunge.h.

Referenced by approach_line(), Mass(), and TRunge().

double TRunge::_maxflightlength [private]

Definition at line 229 of file TRunge.h.

Referenced by Fly(), Fly_SC(), MaxFlightLength(), and SetFlightLength().

TTrackMC* TTrackBase::_mc [protected, inherited]

Definition at line 165 of file TTrackBase.h.

Referenced by TTrackBase::mc().

unsigned TRunge::_ndf [private]

Definition at line 220 of file TRunge.h.

Referenced by ndf(), and reducedchi2().

unsigned TRunge::_Nstep [mutable, private]

Definition at line 249 of file TRunge.h.

Referenced by a(), approach_line(), approach_point(), BfieldID(), Ea(), Eps(), Fly(), Fly_SC(), GetStep(), GetXP(), MaxFlightLength(), Nstep(), pivot(), StepSize(), StepSizeMax(), and StepSizeMin().

HepPoint3D TRunge::_pivot [private]

Definition at line 215 of file TRunge.h.

Referenced by approach(), helix(), pivot(), and SetFirst().

double TRunge::_stepSize [private]

Definition at line 231 of file TRunge.h.

Referenced by approach_line(), approach_point(), Fly(), and StepSize().

double TRunge::_stepSizeMax [private]

Definition at line 235 of file TRunge.h.

Referenced by Fly_SC(), and StepSizeMax().

double TRunge::_stepSizeMin [private]

Definition at line 236 of file TRunge.h.

Referenced by Fly_SC(), and StepSizeMin().

double TRunge::_y[TRunge_MAXstep][6] [mutable, private]

Definition at line 247 of file TRunge.h.

Referenced by approach_line(), approach_point(), Fly(), Fly_SC(), and GetXP().

double TRunge::_yscal[6] [private]

Definition at line 233 of file TRunge.h.

Referenced by Fly_SC(), TRunge(), and Yscal().

const RkFitMaterial TRunge::FitMdcMaterial [private]

Definition at line 239 of file TRunge.h.

IMagneticFieldSvc* TRunge::m_pmgnIMF [private]

Definition at line 251 of file TRunge.h.

Referenced by Function(), SetFlightLength(), and TRunge().

---