FTTrack Class Reference

#include <FTTrack.h>

Public Member Functions
	FTTrack (FTList< FTSegment * > &axial_segments, float kappa, float chi2_kappa)
	constructor
	~FTTrack ()
	destructor
int	r_phiFit (void)
	do r-phi circle fit
int	r_phiReFit (float vx, float xy, int vtx_flag)
	do r-phi refit
int	r_phi2Fit (float vx, float xy, int vtx_flag)
int	r_phi3Fit (int l, float vx, float xy, int vtx_flag)
int	r_phi4Fit (float vx, float xy, int vtx_flag)
int	s_zFit (void)
	do s-z linear fit
void	updateSZ (void)
	update s and z information for linking
int	linkStereoSegments (void)
	link stereo segments by tanLambda
Helix *	helix (void) const
	returns helix parameters
const Lpav &	lpav (void) const
	returns lpav
const zav &	Zav (void) const
	returns zav
FTList< FTSegment * > &	axial_segments (void) const
	returns axial segments
FTList< FTSegment * > &	stereo_segments (void) const
	returns stereo_segments
float	kappa_tmp (void) const
	returns kappa at linking
float	chi2_kappa_tmp (void) const
	returns sigmaKappa at linking
float	d_z (float s, float z) const
void	printout ()
float	SigmaZ (float z)
	add z for culculation of tanLambda
float	SigmaS (float s)
	add s for culculation of tanLambda
float	SigmaSS (float ss)
	add s for culculation of dz, tanLambda
float	SigmaSZ (float sz)
	add s for culculation of dz, tanLambda
void	append_stereo (FTSegment *, float s, float z)
	append stereo segment to the stereo segment list
void	append_stereo_cache (FTSegment *)
	append stereo segment to the cache
int	get_nhits (void)
	calculate the wire hits number
void	setFTFinder (FTFinder *)
	FTTrack (FTList< FTSegment * > &axial_segments, float kappa, float chi2_kappa)
	constructor
	~FTTrack ()
	destructor
int	r_phiFit (void)
	do r-phi circle fit
int	r_phiReFit (float vx, float xy, int vtx_flag)
	do r-phi refit
int	r_phi2Fit (float vx, float xy, int vtx_flag)
int	r_phi3Fit (int l, float vx, float xy, int vtx_flag)
int	r_phi4Fit (float vx, float xy, int vtx_flag)
int	s_zFit (void)
	do s-z linear fit
void	updateSZ (void)
	update s and z information for linking
int	linkStereoSegments (void)
	link stereo segments by tanLambda
Helix *	helix (void) const
	returns helix parameters
const Lpav &	lpav (void) const
	returns lpav
const zav &	Zav (void) const
	returns zav
FTList< FTSegment * > &	axial_segments (void) const
	returns axial segments
FTList< FTSegment * > &	stereo_segments (void) const
	returns stereo_segments
float	kappa_tmp (void) const
	returns kappa at linking
float	chi2_kappa_tmp (void) const
	returns sigmaKappa at linking
float	d_z (float s, float z) const
void	printout ()
float	SigmaZ (float z)
	add z for culculation of tanLambda
float	SigmaS (float s)
	add s for culculation of tanLambda
float	SigmaSS (float ss)
	add s for culculation of dz, tanLambda
float	SigmaSZ (float sz)
	add s for culculation of dz, tanLambda
void	append_stereo (FTSegment *, float s, float z)
	append stereo segment to the stereo segment list
void	append_stereo_cache (FTSegment *)
	append stereo segment to the cache
int	get_nhits (void)
	calculate the wire hits number
void	setFTFinder (FTFinder *)

Public Attributes
float	minPt
float	minDr
float	xtCoEff
float	Testime
Lpav *	_la
zav *	_za
Helix *	_helix
FTList< FTSegment * > &	_axial_segments
FTList< FTSegment * > *	_stereo_segments
FTList< FTSegment * > *	_stereo_segments_cache
FTList< FTList< FTSegment * > * > *	_stereo_segments_by_superLayer
union {
float   _kappa
float   _SigmaS
};
union {
float   _chi2_kappa
float   _SigmaSS
};
float	_SigmaZ
float	_SigmaSZ
union {
float   _kappa
float   _SigmaS
};
union {
float   _chi2_kappa
float   _SigmaSS
};

Static Public Attributes
static MdcParameter *	param = MdcParameter::instance()

Detailed Description

Definition at line 27 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Constructor & Destructor Documentation

FTTrack() [1/2]

FTTrack::FTTrack ( FTList< FTSegment * > &  axial_segments, float  kappa, float  chi2_kappa  )

inline

constructor

Definition at line 160 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

  : _la(NULL),
    _za(NULL),
    _helix(NULL),
    _axial_segments(axial_segments),
    _stereo_segments(NULL),
    _stereo_segments_cache(NULL),
    _stereo_segments_by_superLayer(NULL),
    //_wireHits(*(new FTList<FTWire *>(50))),
    _kappa(kappa),
    _chi2_kappa(chi2_kappa)
{
  StatusCode scmgn = Gaudi::svcLocator()->service("MagneticFieldSvc",m_pmgnIMF);
  if(scmgn!=StatusCode::SUCCESS) {
    std::cout<< "Unable to open Magnetic field service"<<std::endl;
  }
}

~FTTrack() [1/2]

FTTrack::~FTTrack ( )

inline

destructor

Definition at line 180 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  delete &_axial_segments;
  delete _la;
  delete _za;
  delete _stereo_segments;
  delete _helix;
}

FTTrack() [2/2]

FTTrack::FTTrack	(	FTList< FTSegment * > &	axial_segments,
		float	kappa,
		float	chi2_kappa
	)

constructor

~FTTrack() [2/2]

FTTrack::~FTTrack

(

)

destructor

Member Function Documentation

append_stereo() [1/2]

void FTTrack::append_stereo ( FTSegment *  src, float  s = 0, float  z = 0  )

inline

append stereo segment to the stereo segment list

Definition at line 212 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  _stereo_segments->append(src);
  _SigmaS += s;
  _SigmaZ += z;
  _SigmaSZ += s*z;
  _SigmaSS += s*s;
}

append_stereo() [2/2]

void FTTrack::append_stereo	(	FTSegment *	,
		float	s,
		float	z
	)

append stereo segment to the stereo segment list

append_stereo_cache() [1/2]

void FTTrack::append_stereo_cache ( FTSegment *  src )

inline

append stereo segment to the cache

Definition at line 205 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  _stereo_segments_cache->append(src);
}

append_stereo_cache() [2/2]

void FTTrack::append_stereo_cache

(

FTSegment *

)

append stereo segment to the cache

axial_segments() [1/2]

FTList< FTSegment * > & FTTrack::axial_segments ( void  ) const

inline

returns axial segments

Definition at line 191 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _axial_segments;
}

Referenced by printout().

axial_segments() [2/2]

FTList< FTSegment * > & FTTrack::axial_segments

(

void

)

const

returns axial segments

chi2_kappa_tmp() [1/2]

float FTTrack::chi2_kappa_tmp ( void  ) const

inline

returns sigmaKappa at linking

Definition at line 230 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _chi2_kappa;
}

chi2_kappa_tmp() [2/2]

float FTTrack::chi2_kappa_tmp

(

void

)

const

returns sigmaKappa at linking

d_z() [1/2]

float FTTrack::d_z ( float  s, float  z  ) const

inline

returns difference between z and estimated z at stereo segment linking if its valid

Definition at line 286 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  int n = _stereo_segments->length();
  if (!n) return 0;
  float Dz = (n==1)
    ? z - s*_SigmaZ/_SigmaS
    : z - (s*(_SigmaSZ-_SigmaS*_SigmaZ)+_SigmaSS*_SigmaZ-_SigmaS*_SigmaSZ)
      / (_SigmaSS - _SigmaS*_SigmaS);
  return (Dz < 10*(7-n)) ? Dz : 9999.;
}

Referenced by linkStereoSegments(), and updateSZ().

d_z() [2/2]

float FTTrack::d_z	(	float	s,
		float	z
	)		const

returns difference between z and estimated z at stereo segment linking if its valid

get_nhits() [1/2]

int FTTrack::get_nhits

(

void

)

calculate the wire hits number

Definition at line 484 of file FTTrack.cxx.

{
  int nhits = 0;
  int n = _axial_segments.length();
  for (int i = 0; i^n; i++) {
    nhits += _axial_segments[i]->wireHits().length();
  }
  n = _stereo_segments->length();
  for (int j = 0; j^n; j++) {
    nhits += (*_stereo_segments)[j]->wireHits().length();
  }
  return nhits;
}

get_nhits() [2/2]

int FTTrack::get_nhits

(

void

)

calculate the wire hits number

helix() [1/2]

Helix * FTTrack::helix ( void  ) const

inline

returns helix parameters

Definition at line 251 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _helix;
}

helix() [2/2]

Helix * FTTrack::helix

(

void

)

const

returns helix parameters

kappa_tmp() [1/2]

float FTTrack::kappa_tmp ( void  ) const

inline

returns kappa at linking

Definition at line 223 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _kappa;
}

kappa_tmp() [2/2]

float FTTrack::kappa_tmp

(

void

)

const

returns kappa at linking

linkStereoSegments() [1/2]

int FTTrack::linkStereoSegments

(

void

)

link stereo segments by tanLambda

Definition at line 342 of file FTTrack.cxx.

                               {
  delete _stereo_segments_cache;
  int n = _stereo_segments_by_superLayer->length();
  for (int i = 0; i^n; i++){
    FTList<FTSegment *> & segments = *(*_stereo_segments_by_superLayer)[i];
    int m = segments.length();
    float min_D_z = 9998.;
    float S = 0.;
    float Z = 0.;
    FTSegment * selected = NULL;
    for (int j = 0; j^m; j++){
      FTSegment * s = segments[j];
      float s_tmp = s->s();
      float z_tmp = s->z();
      double D_z = fabs(d_z(s_tmp,z_tmp));
      if (D_z < min_D_z){
    selected = s;
    min_D_z = D_z;
    S = s_tmp;
    Z = z_tmp;
      }
    }
    if (selected){
      _stereo_segments->append(selected);
      _SigmaS += S;
      _SigmaZ += Z;
      _SigmaSZ += S*Z;
      _SigmaSS += S*S;
    }
  }
  _stereo_segments_by_superLayer->deleteAll();
  delete _stereo_segments_by_superLayer;
  return 0;
}

linkStereoSegments() [2/2]

int FTTrack::linkStereoSegments

(

void

)

link stereo segments by tanLambda

lpav() [1/2]

const Lpav & FTTrack::lpav ( void  ) const

inline

returns lpav

Definition at line 237 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return *_la;
}

Referenced by FTSegment::update3D().

lpav() [2/2]

const Lpav & FTTrack::lpav

(

void

)

const

returns lpav

printout() [1/2]

void FTTrack::printout

(

)

printout _tracks infomation added by X.-R. Lu

Definition at line 499 of file FTTrack.cxx.

                      {
  int n = axial_segments().length();
  for (int i = 0; i^n; i++) axial_segments()[i]->printout();
  n = stereo_segments().length();
  for (int i = 0; i^n; i++) stereo_segments()[i]->printout();  
}

printout() [2/2]

void FTTrack::printout

(

)

printout _tracks infomation added by X.-R. Lu

r_phi2Fit() [1/2]

int FTTrack::r_phi2Fit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

Definition at line 191 of file FTTrack.cxx.

{
  if (vtx_flag) _la->fit(vx, vy, 20*_la->chisq()/_la->nc());
  int n = _axial_segments.length();
  _la->clear();
  int k=0,l=0;
  float temp=0;
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
     int m = hits.length();
     for (int j = 0; j^m; j++){
      k++;
      FTWire & h = * hits[j];
      const float x = h.x();
      const float y = h.y();
      double d0 = _la->d((double)x,(double)y);
      float cellsize = h.layer().csize();
      if (m_ftFinder->evtTiming){
    float time = h.time() + m_ftFinder->evtTiming ;
    if(time<-100) continue;
    h.distance(m_ftFinder->t2x(h.layer(),time));
    //  h.distance(time*40*0.0001);
      }
      float delta = h.distance()/h.layer().r();
      if (fabs(d0) > 0.5*cellsize) continue; // remove bad hits
      if (h.distance()<0.1 || h.distance()>0.6) continue;
      if (d0>0){            // left or right
    h.stateOR(FTWireHitRight);
    d0 = -d0;
      }else{
    h.stateOR(FTWireHitLeft);
    delta = -delta;
      }
#ifndef OnlineMode
      //fill histogram 
       //g_sigmaxy->fill(h.distance()+d0, 1.0);
#endif
      _la->add_point(x - delta*y, y + delta*x, 1);
      if(temp<fabs(h.distance()+d0)){
        temp=fabs(h.distance()+d0);
    l=k;
      }
    }
  }
  HepVector b = _la->Hpar(m_ftFinder->pivot);
  double  chi2 = _la->fit();            // refit using drift distance
#ifndef OnlineMode
  // g_chi2xy->fill(chi2/_la->nc(), 1.0);
#endif
  if(k>5){    
    return  l;
  }
  else{
    return -99;
  }
}

r_phi2Fit() [2/2]

int FTTrack::r_phi2Fit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

r_phi3Fit() [1/2]

int FTTrack::r_phi3Fit	(	int	l,
		float	vx,
		float	xy,
		int	vtx_flag
	)

Definition at line 249 of file FTTrack.cxx.

{
  if (vtx_flag) _la->fit(vx, vy, 20*_la->chisq()/_la->nc());
  int n = _axial_segments.length();
  _la->clear();
  int k=0;
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
    int m = hits.length();
    for (int j = 0; j^m; j++){
      k++;
      FTWire & h = * hits[j];
      const float x = h.x();
      const float y = h.y();
      double d0 = _la->d((double)x,(double)y);
      float cellsize = h.layer().csize();
      if (m_ftFinder->evtTiming){
    float time = h.time() + m_ftFinder->evtTiming ;
    if(time<-100) continue;
    h.distance(m_ftFinder->t2x(h.layer(),time));
    //h.distance(time*40*0.0001);
      }
      float delta = h.distance()/h.layer().r();
      if (fabs(d0) > 0.5*cellsize) continue; // remove bad hits
      if (d0>0){            // left or right
    h.stateOR(FTWireHitRight);
    d0 = -d0;
      }else{
    h.stateOR(FTWireHitLeft);
    delta = -delta;
      }
      if (k==l){
    //delete hits[j];
    m=hits.remove(j);
    continue;
      }
      _la->add_point(x - delta*y, y + delta*x, 1);
      // h.setChi2(h.distance()+d0);
    }
  }
  HepVector b = _la->Hpar(m_ftFinder->pivot);
  double chi2 = _la->fit();            // refit using drift distance
  return 1;
}

r_phi3Fit() [2/2]

int FTTrack::r_phi3Fit	(	int	l,
		float	vx,
		float	xy,
		int	vtx_flag
	)

r_phi4Fit() [1/2]

int FTTrack::r_phi4Fit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

Definition at line 295 of file FTTrack.cxx.

{
  if (vtx_flag) _la->fit(vx, vy, 20*_la->chisq()/_la->nc());
  int n = _axial_segments.length();
  _la->clear();
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
    int m = hits.length();
    for (int j = 0; j^m; j++){
      FTWire & h = * hits[j];
      const float x = h.x();
      const float y = h.y();
      double d0 = _la->d((double)x,(double)y);
      float cellsize = h.layer().csize();
      if (m_ftFinder->evtTiming){
    float time = h.time() + m_ftFinder->evtTiming ;
    if(time<-100) continue;
    h.distance(m_ftFinder->t2x(h.layer(),time));
    //h.distance(time*40*0.0001);
      }
      float delta = h.distance()/h.layer().r();
      if (fabs(d0) > 0.5*cellsize) continue; // remove bad hits
      if (d0>0){            // left or right
    h.stateOR(FTWireHitRight);
    d0 = -d0;
      }else{
    h.stateOR(FTWireHitLeft);
    delta = -delta;
      }
    
      _la->add_point(x - delta*y, y + delta*x, 1);
      h.setChi2(h.distance()+d0);
#ifndef OnlineMode
      //fill histogram 
      // g_sigmaxy->fill(h.distance()+d0, 1.0);
#endif
    }
  }
  HepVector b = _la->Hpar(m_ftFinder->pivot);
  double chi2 = _la->fit();            // refit using drift distance
#ifndef OnlineMode
 // g_chi2xy->fill(chi2/_la->nc(), 1.0);
#endif
  return 1;
}

r_phi4Fit() [2/2]

int FTTrack::r_phi4Fit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

r_phiFit() [1/2]

int FTTrack::r_phiFit

(

void

)

do r-phi circle fit

Definition at line 37 of file FTTrack.cxx.

{
  IMessageSvc *msgSvc;
  Gaudi::svcLocator()->service("MessageSvc", msgSvc);
 
  MsgStream log(msgSvc, "FTFinder");
 
  // static const float alpha(333.564095);
  if (fabs(_kappa) > 1.2/param->_minPt) return 0;
 
  _la = new Lpav;
  int n = _axial_segments.length();
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
    int m = hits.length();
    for (int j = 0; j^m; j++){
      FTWire & h = * hits[j];
      if (h.stateAND(FTWireFittingInvalid)){
    hits.remove(j);
    continue;
      }
      //      log << MSG::DEBUG << "layer: "<< h.layer().layerId() << " phi: " << h.phi() << " df_distance: "<< h.distance() << endreq;
      double par = h.distance()/(0.25*h.layer().csize());
      _la->add_point((double)h.x(),(double)h.y(),exp(-par*par));
    }
  }
 
  double chi2 = _la->fit();
  HepVector a = _la->Hpar(m_ftFinder->pivot);
  log << MSG::DEBUG << " chi2/_la cut(1): " << chi2 << " / " << _la->nc() 
      <<" a1(" << param->_minDr << "): " << a(1) 
      <<" a3(" << 1.05/param->_minPt << "):" << a(3) <<endreq;
  if (chi2/_la->nc() > 1.) return 0;
  if (fabs(a(3))>(1.05/param->_minPt)) return 0; 
  if (fabs(a(1))>param->_minDr) return 0;
  _la->clear();
  log << MSG::DEBUG << " passed chi2/_la, a(3) and a(1) cut" <<endreq;
  int layer0 = 0;
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
    if (!_axial_segments[i]->superLayer().superLayerId()) layer0 = 1;
    int m = hits.length();
    for (int j = 0; j^m; j++){
      FTWire & h = * hits[j];
      const float x = h.x();
      const float y = h.y();
      double d0 = _la->d((double)x,(double)y);
      float delta = h.distance()/h.layer().r();
      if (fabs(d0) > 0.7*h.layer().csize()) continue; // remove bad hits
      if (d0>0){            // left or right
    d0 = -d0;
      }else{
    delta = -delta;
      }
      _la->add_point(x - delta*y, y + delta*x, 1);
    }
  }
  if (!layer0){      // salvage hits from complecated segments in layer0
    FTList<FTSegment *> & salvage = 
      m_ftFinder->superLayer(0)->complecated_segments();
    HepVector center = _la->center();
    const float xc = center(1);
    const float yc = center(2);
    const float rc = a(1)+(-1. / 2.99792458 /m_pmgnIMF->getReferField())/a(3); // rho+drho(signed)
    int nn = salvage.length();
    for (int i = 0; i^nn; i++){
      int salvaged = 0;
      FTList<FTWire *> & hits = salvage[i]->wireHits();
      int m = hits.length();
      for (int j = 0; j^m; j++){
    FTWire & h = * hits[j];
    float x = h.x();
    float y = h.y();
    float r = h.layer().r();
    if ((y*xc-x*yc)/(r*rc)<0.707) break;
    double d0 = _la->d((double)x,(double)y);
    float delta = h.distance()/r;
    if (fabs(d0) > 0.7*h.layer().csize()) continue; // remove bad hits
    if (d0>0){          // left or right
      h.stateOR(FTWireHitRight);
      d0 = -d0;
    }else{
      h.stateOR(FTWireHitLeft);
      delta = -delta;
    }
    _la->add_point(x - delta*y, y + delta*x, 1);
    salvaged = 1;
      }
      if (salvaged){
    _axial_segments.append(salvage[i]);
    break;
      }
    }
  }
  chi2 = _la->fit();        // refit using drift distance
  _stereo_segments = new FTList<FTSegment *>(6);  
  _stereo_segments_cache = new FTList<FTSegment *>(3);
  _stereo_segments_by_superLayer = new FTList<FTList<FTSegment *> *>(6);
  _za = new zav;
  _SigmaS = 0.;
  _SigmaZ = 0.;
  _SigmaSS = 0.;
  _SigmaSZ = 0.;
  return 1;
}

r_phiFit() [2/2]

int FTTrack::r_phiFit

(

void

)

do r-phi circle fit

r_phiReFit() [1/2]

int FTTrack::r_phiReFit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

do r-phi refit

Definition at line 144 of file FTTrack.cxx.

{
  if (vtx_flag) _la->fit(vx, vy, 20*_la->chisq()/_la->nc());
  int n = _axial_segments.length();
  _la->clear();
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = _axial_segments[i]->wireHits();
    int m = hits.length();
    for (int j = 0; j^m; j++){
      FTWire & h = * hits[j];
      const float x = h.x();
      const float y = h.y();
      double d0 = _la->d((double)x,(double)y);
      float cellsize = h.layer().csize();
      if (m_ftFinder->evtTiming){
    float time = h.time() + m_ftFinder->evtTiming;
    if(time<-100) continue; // discard the bad TDC time
    h.distance(m_ftFinder->t2x(h.layer(),time));
    //h.distance(time*40*0.0001);
      }
      float delta = h.distance()/h.layer().r();
      if (fabs(d0) > 0.5*cellsize) continue; // remove bad hits
      if (d0>0){            // left or right
    h.stateOR(FTWireHitRight);
    d0 = -d0;
      }else{
    h.stateOR(FTWireHitLeft);
    delta = -delta;
      }
      h.setChi2(h.distance()+d0);
#ifndef OnlineMode
      //fill histogram 
      g_sigmaxy->fill(h.distance()+d0, 1.0);
#endif
      _la->add_point(x - delta*y, y + delta*x, 1);
    }
  }
  HepVector b = _la->Hpar(m_ftFinder->pivot);
  double  chi2 = _la->fit();            // refit using drift distance
#ifndef OnlineMode
  g_chi2xy->fill(chi2/_la->nc(), 1.0);
#endif
      
  return 1;
} 

r_phiReFit() [2/2]

int FTTrack::r_phiReFit	(	float	vx,
		float	xy,
		int	vtx_flag
	)

do r-phi refit

s_zFit() [1/2]

int FTTrack::s_zFit

(

void

)

do s-z linear fit

Definition at line 378 of file FTTrack.cxx.

                   {
  IMessageSvc *msgSvc;
  Gaudi::svcLocator()->service("MessageSvc", msgSvc);
 
  MsgStream log(msgSvc, "FTFinder");
 
  HepVector a = _la->Hpar(m_ftFinder->pivot);
  int n = _stereo_segments->length();
  log<<MSG::DEBUG<<"number of stereo segments: "<< n << endreq;
  if (n < (param->_nseg)){
    a(4) = -9999.;
    a(5) = -9999.;
    _helix = new Helix(m_ftFinder->pivot,a);
    log<<MSG::DEBUG<<"cut by _nseg" << param->_nseg << endreq;
    return 0;
  }
  FTList<double> zList(50);
  FTList<double> sList(50);
  FTList<FTWire *> hList(50);
  for (int i = 0; i^n; i++){
    FTList<FTWire *> & hits = (*_stereo_segments)[i]->wireHits();
    int m = hits.length();
    for (int j = 0; j^m; j++){
      FTWire & h = * hits[j];
      double z;
      if (!(h.z(*_la,z))) continue;
      double s = _la->s(h.layer().r());
      float cellsize = h.layer().csize();
      log<<MSG::DEBUG<<"cellsize: " << cellsize << endreq;
 
      if (m_ftFinder->evtTiming){
    float time = h.time() + m_ftFinder->evtTiming;
    if(time<-100) continue; //discard the bad TDC time
    double distance = m_ftFinder->t2x(h.layer(),time);
    h.distance(distance);
    log<<MSG::DEBUG<<"m_ftFinder->evtTiming: "<< m_ftFinder->evtTiming << " TDC time: " << h.time() << " distance: "<< distance << endreq;
      }
 
      double par = h.distance()/(0.25*cellsize);
      _za->add(s,z,exp(-par*par));
      sList.append(s);
      zList.append(z);
      hList.append(&h);
    }
  }
  if (hList.length() < (param->_nlength)){
    a(4) = -9999.;
    a(5) = -9999.;
    _helix = new Helix(m_ftFinder->pivot,a);
    log << MSG::DEBUG << "cut by _nlength: " << param->_nlength << endreq;
    return 0;
  }   
  double chi2 = _za->calculate();
  _za->clear();
  n = hList.length();
  for (int i = 0; i^n; i++){
    double d = _za->d(sList[i],zList[i]);
    float z_distance = hList[i]->distance_z();
    if (fabs(fabs(d)-z_distance) > (param->_z_cut1)){ 
      log<<MSG::DEBUG<<"cut by _z_cut1: "<< param->_z_cut1 << endreq;
      zList.remove2(i);
      sList.remove2(i);
      n = hList.remove(i);
      continue;
    }
    _za->add(sList[i], (d>0) ? zList[i]-z_distance : zList[i]+z_distance, 1.);
  }
 
  if (_za->nc() < (param->_nc)){
    a(4) = -9999.;
    a(5) = -9999.;
    _helix = new Helix(m_ftFinder->pivot,a);
    log<<MSG::DEBUG<<"cut by _nc: " << param->_nc << endreq;
    return 0;
  }    
  chi2 = _za->calculate();
  _za->clear();
  for (int i = 0; i^n; i++){
    double d = _za->d(sList[i],zList[i]);
    float z_distance = hList[i]->distance_z();
    hList[i]->setChi2(fabs(d)-z_distance);
 #ifndef OnlineMode
    //fill ntuple
    g_sigmaz->fill(fabs(d)-z_distance, 1.0);
#endif
    if (fabs(fabs(d)-z_distance) > (param->_z_cut2)) continue; 
    _za->add(sList[i], (d>0) ? zList[i]-z_distance : zList[i]+z_distance, 1.);
  }
 
  if (_za->nc() < (param->_nc)){
    a(4) = -9999.;
    a(5) = -9999.;
    _helix = new Helix(m_ftFinder->pivot,a);
    log<<MSG::DEBUG<<"cut by _nc" << param->_nc << endreq;
    return 0;
  }    
  chi2 = _za->calculate();
#ifndef OnlineMode
  g_chi2sz->fill(chi2/_za->nc(), 1.0);
#endif
  a(4) = _za->b();      // dz
  a(5) = _za->a();      // tanLambda
  _helix = new Helix(m_ftFinder->pivot,a);
  return 1;
}

s_zFit() [2/2]

int FTTrack::s_zFit

(

void

)

do s-z linear fit

setFTFinder() [1/2]

void FTTrack::setFTFinder ( FTFinder *  finder )

inline

Definition at line 323 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  m_ftFinder = finder;
}

setFTFinder() [2/2]

void FTTrack::setFTFinder

(

FTFinder *

)

SigmaS() [1/2]

float FTTrack::SigmaS ( float  s )

inline

add s for culculation of tanLambda

Definition at line 265 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _SigmaS += s;
}

SigmaS() [2/2]

float FTTrack::SigmaS

(

float

s

)

add s for culculation of tanLambda

SigmaSS() [1/2]

float FTTrack::SigmaSS ( float  ss )

inline

add s for culculation of dz, tanLambda

Definition at line 279 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _SigmaSS += ss;
}

SigmaSS() [2/2]

float FTTrack::SigmaSS

(

float

ss

)

add s for culculation of dz, tanLambda

SigmaSZ() [1/2]

float FTTrack::SigmaSZ ( float  sz )

inline

add s for culculation of dz, tanLambda

Definition at line 272 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _SigmaSZ += sz;
}

SigmaSZ() [2/2]

float FTTrack::SigmaSZ

(

float

sz

)

add s for culculation of dz, tanLambda

SigmaZ() [1/2]

float FTTrack::SigmaZ ( float  z )

inline

add z for culculation of tanLambda

Definition at line 258 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return _SigmaZ += z;
}

SigmaZ() [2/2]

float FTTrack::SigmaZ

(

float

z

)

add z for culculation of tanLambda

stereo_segments() [1/2]

FTList< FTSegment * > & FTTrack::stereo_segments ( void  ) const

inline

returns stereo_segments

Definition at line 198 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return *_stereo_segments;
}

Referenced by printout().

stereo_segments() [2/2]

FTList< FTSegment * > & FTTrack::stereo_segments

(

void

)

const

returns stereo_segments

updateSZ() [1/2]

void FTTrack::updateSZ ( void  )

inline

update s and z information for linking

Definition at line 300 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  int n = _stereo_segments_cache->length();
  if (n==1){
    FTSegment * s = _stereo_segments_cache->first();
    float s_tmp = s->s();
    float z_tmp = s->z();
    if (d_z(s_tmp,z_tmp) < 9998.){
      _SigmaS += s_tmp;
      _SigmaZ += z_tmp;
      _SigmaSS += s_tmp*s_tmp;
      _SigmaSZ += s_tmp*z_tmp;
      _stereo_segments->append(s);
    }
    _stereo_segments_cache->clear();
  }else if(n){
    _stereo_segments_by_superLayer->append(_stereo_segments_cache);
    _stereo_segments_cache = new FTList<FTSegment *>(3);
  }
}

updateSZ() [2/2]

void FTTrack::updateSZ

(

void

)

update s and z information for linking

Zav() [1/2]

const zav & FTTrack::Zav ( void  ) const

inline

returns zav

Definition at line 244 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

{
  return *_za;
}

Zav() [2/2]

const zav & FTTrack::Zav

(

void

)

const

returns zav

Member Data Documentation

union { ... } FTTrack::@36

union { ... } FTTrack::@38

union { ... } FTTrack::@68

union { ... } FTTrack::@70

_axial_segments

FTList< FTSegment * > & FTTrack::_axial_segments

Definition at line 130 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by axial_segments(), get_nhits(), r_phi2Fit(), r_phi3Fit(), r_phi4Fit(), r_phiFit(), r_phiReFit(), and ~FTTrack().

_chi2_kappa

float FTTrack::_chi2_kappa

Definition at line 141 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by chi2_kappa_tmp().

_helix

Helix * FTTrack::_helix

Definition at line 129 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by helix(), s_zFit(), and ~FTTrack().

_kappa

float FTTrack::_kappa

Definition at line 137 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by kappa_tmp(), and r_phiFit().

_la

Lpav * FTTrack::_la

Definition at line 127 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by lpav(), r_phi2Fit(), r_phi3Fit(), r_phi4Fit(), r_phiFit(), r_phiReFit(), s_zFit(), and ~FTTrack().

_SigmaS

float FTTrack::_SigmaS

Definition at line 138 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo(), d_z(), linkStereoSegments(), r_phiFit(), SigmaS(), and updateSZ().

_SigmaSS

float FTTrack::_SigmaSS

Definition at line 142 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo(), d_z(), linkStereoSegments(), r_phiFit(), SigmaSS(), and updateSZ().

_SigmaSZ

float FTTrack::_SigmaSZ

Definition at line 145 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo(), d_z(), linkStereoSegments(), r_phiFit(), SigmaSZ(), and updateSZ().

_SigmaZ

float FTTrack::_SigmaZ

Definition at line 144 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo(), d_z(), linkStereoSegments(), r_phiFit(), SigmaZ(), and updateSZ().

_stereo_segments

FTList< FTSegment * > * FTTrack::_stereo_segments

Definition at line 131 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo(), d_z(), get_nhits(), linkStereoSegments(), r_phiFit(), s_zFit(), stereo_segments(), updateSZ(), and ~FTTrack().

_stereo_segments_by_superLayer

FTList< FTList< FTSegment * > * > * FTTrack::_stereo_segments_by_superLayer

Definition at line 133 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by linkStereoSegments(), r_phiFit(), and updateSZ().

_stereo_segments_cache

FTList< FTSegment * > * FTTrack::_stereo_segments_cache

Definition at line 132 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by append_stereo_cache(), linkStereoSegments(), r_phiFit(), and updateSZ().

_za

zav * FTTrack::_za

Definition at line 128 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by r_phiFit(), s_zFit(), Zav(), and ~FTTrack().

minDr

float FTTrack::minDr

Definition at line 113 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

minPt

float FTTrack::minPt

Definition at line 112 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

param

MdcParameter * FTTrack::param = MdcParameter::instance()

static

Definition at line 134 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

Referenced by r_phiFit(), and s_zFit().

Testime

float FTTrack::Testime

Definition at line 117 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

xtCoEff

float FTTrack::xtCoEff

Definition at line 116 of file InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h.

---

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

• InstallArea/include/MdcFastTrkAlg/MdcFastTrkAlg/FTTrack.h
• Reconstruction/MdcFastTrkAlg/MdcFastTrkAlg-00-04-09/MdcFastTrkAlg/FTTrack.h
• FTTrack.cxx