#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 *)

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

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

Detailed Description

Definition at line 27 of file FTTrack.h.

Constructor & Destructor Documentation

◆ FTTrack()

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

inline

constructor

Definition at line 160 of file 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()

FTTrack::~FTTrack ( )

inline

destructor

Definition at line 180 of file FTTrack.h.

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

Member Function Documentation

◆ append_stereo()

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 FTTrack.h.

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

◆ append_stereo_cache()

void FTTrack::append_stereo_cache ( FTSegment * src )

inline

append stereo segment to the cache

Definition at line 205 of file FTTrack.h.

{
  _stereo_segments_cache->append(src);
}

◆ axial_segments()

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

inline

returns axial segments

Definition at line 191 of file FTTrack.h.

{
  return _axial_segments;
}

Referenced by printout().

◆ chi2_kappa_tmp()

float FTTrack::chi2_kappa_tmp ( void ) const

inline

returns sigmaKappa at linking

Definition at line 230 of file FTTrack.h.

{
  return _chi2_kappa;
}

◆ d_z()

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

◆ get_nhits()

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

◆ helix()

Helix * FTTrack::helix ( void ) const

inline

returns helix parameters

Definition at line 251 of file FTTrack.h.

{
  return _helix;
}

◆ kappa_tmp()

float FTTrack::kappa_tmp ( void ) const

inline

returns kappa at linking

Definition at line 223 of file FTTrack.h.

{
  return _kappa;
}

◆ linkStereoSegments()

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

◆ lpav()

const Lpav & FTTrack::lpav ( void ) const

inline

returns lpav

Definition at line 237 of file FTTrack.h.

{
  return *_la;
}

Referenced by FTSegment::update3D().

◆ printout()

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();  
}

◆ r_phi2Fit()

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_phi3Fit()

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_phi4Fit()

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_phiFit()

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_phiReFit()

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

◆ s_zFit()

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

◆ setFTFinder()

void FTTrack::setFTFinder ( FTFinder * finder )

inline

Definition at line 323 of file FTTrack.h.

{
  m_ftFinder = finder;
}

◆ SigmaS()

float FTTrack::SigmaS ( float s )

inline

add s for culculation of tanLambda

Definition at line 265 of file FTTrack.h.

{
  return _SigmaS += s;
}

◆ SigmaSS()

float FTTrack::SigmaSS ( float ss )

inline

add s for culculation of dz, tanLambda

Definition at line 279 of file FTTrack.h.

{
  return _SigmaSS += ss;
}

◆ SigmaSZ()

float FTTrack::SigmaSZ ( float sz )

inline

add s for culculation of dz, tanLambda

Definition at line 272 of file FTTrack.h.

{
  return _SigmaSZ += sz;
}

◆ SigmaZ()

float FTTrack::SigmaZ ( float z )

inline

add z for culculation of tanLambda

Definition at line 258 of file FTTrack.h.

{
  return _SigmaZ += z;
}

◆ stereo_segments()

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

inline

returns stereo_segments

Definition at line 198 of file FTTrack.h.

{
  return *_stereo_segments;
}

Referenced by printout().

◆ updateSZ()

void FTTrack::updateSZ ( void )

inline

update s and z information for linking

Definition at line 300 of file 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);
  }
}

◆ Zav()

const zav & FTTrack::Zav ( void ) const

inline

returns zav

Definition at line 244 of file FTTrack.h.

{
  return *_za;
}

Member Data Documentation

◆ [union]

union { ... } FTTrack

◆ [union]

union { ... } FTTrack

◆ _axial_segments

FTList<FTSegment *>& FTTrack::_axial_segments

Definition at line 130 of file 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 FTTrack.h.

Referenced by chi2_kappa_tmp().

◆ _helix

Helix* FTTrack::_helix

Definition at line 129 of file FTTrack.h.

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

◆ _kappa

float FTTrack::_kappa

Definition at line 137 of file FTTrack.h.

Referenced by kappa_tmp(), and r_phiFit().

◆ _la

Lpav* FTTrack::_la

Definition at line 127 of file 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 FTTrack.h.

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

◆ _SigmaSS

float FTTrack::_SigmaSS

Definition at line 142 of file FTTrack.h.

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

◆ _SigmaSZ

float FTTrack::_SigmaSZ

Definition at line 145 of file FTTrack.h.

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

◆ _SigmaZ

float FTTrack::_SigmaZ

Definition at line 144 of file 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 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 FTTrack.h.

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

◆ _stereo_segments_cache

FTList<FTSegment *>* FTTrack::_stereo_segments_cache

Definition at line 132 of file FTTrack.h.

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

◆ _za

zav* FTTrack::_za

Definition at line 128 of file FTTrack.h.

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

◆ minDr

float FTTrack::minDr

Definition at line 113 of file FTTrack.h.

◆ minPt

float FTTrack::minPt

Definition at line 112 of file FTTrack.h.

◆ param

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

static

Definition at line 134 of file FTTrack.h.

Referenced by r_phiFit(), and s_zFit().

◆ Testime

float FTTrack::Testime

Definition at line 117 of file FTTrack.h.

◆ xtCoEff

float FTTrack::xtCoEff

Definition at line 116 of file FTTrack.h.

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

Public Member Functions

Public Attributes

Static Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FTTrack()

◆ ~FTTrack()

Member Function Documentation

◆ append_stereo()

◆ append_stereo_cache()

◆ axial_segments()

◆ chi2_kappa_tmp()

◆ d_z()

◆ get_nhits()

◆ helix()

◆ kappa_tmp()

◆ linkStereoSegments()

◆ lpav()

◆ printout()

◆ r_phi2Fit()

◆ r_phi3Fit()

◆ r_phi4Fit()

◆ r_phiFit()

◆ r_phiReFit()

◆ s_zFit()

◆ setFTFinder()

◆ SigmaS()

◆ SigmaSS()

◆ SigmaSZ()

◆ SigmaZ()

◆ stereo_segments()

◆ updateSZ()

◆ Zav()

Member Data Documentation

◆ [union]

◆ [union]

◆ _axial_segments

◆ _chi2_kappa

◆ _helix

◆ _kappa

◆ _la

◆ _SigmaS

◆ _SigmaSS

◆ _SigmaSZ

◆ _SigmaZ

◆ _stereo_segments

◆ _stereo_segments_by_superLayer

◆ _stereo_segments_cache

◆ _za

◆ minDr

◆ minPt

◆ param

◆ Testime

◆ xtCoEff