PndFTSCATrackParam Class Reference

#include <PndFTSCATrackParam.h>

Public Member Functions
	PndFTSCATrackParam ()
	PndFTSCATrackParam (const TrackParamVector &v, int i)
float	X0 () const
float	X1 () const
float	X2 () const
float	X () const
float	Y () const
float	Z () const
float	SinPhi () const
float	DzDs () const
float	QPt () const
float	QP () const
float	QMomentum () const
float	SignCosPhi () const
float	Chi2 () const
int	NDF () const
float	Err2X1 () const
float	Err2X2 () const
float	Err2Y () const
float	Err2Z () const
float	Err2SinPhi () const
float	Err2DzDs () const
float	Err2QPt () const
float	Angle () const
float	Kappa (float Bz) const
float	CosPhi () const
float	Err2QMomentum () const
const float *	Par () const
const float *	Cov () const
bool	GetXYZPxPyPzQ (float &x, float &y, float &z, float &px, float &py, float &pz, int &q, float cov[21]) const
void	SetSinPhi (float v)
void	GetDCAPoint (float x, float y, float z, float &px, float &py, float &pz, float Bz) const
bool	TransportToX0 (float x, float Bz, float maxSinPhi=.999)
bool	Rotate (float alpha, float maxSinPhi=.999)
bool	Transport (const FTSCAHit &hit, const PndFTSCAParam &param)
bool	IsValid () const
void	SetAsInvalid ()

Private Member Functions
void	Reset ()
float	S (float x, float y, float Bz) const

Private Attributes
float	fX
float	fSignCosPhi
float	fP [5]
float	fC [15]
float	fChi2
int	fNDF
float	fAlpha

Friends
std::istream &	operator>> (std::istream &, PndFTSCATrackParam &)
std::ostream &	operator<< (std::ostream &, const PndFTSCATrackParam &)

Detailed Description

Definition at line 245 of file PndFTSCATrackParam.h.

Constructor & Destructor Documentation

PndFTSCATrackParam::PndFTSCATrackParam ( )

inline

Definition at line 251 of file PndFTSCATrackParam.h.

References Reset().

{ Reset();}

PndFTSCATrackParam::PndFTSCATrackParam ( const TrackParamVector &  v, int  i  )

inline

Definition at line 252 of file PndFTSCATrackParam.h.

References PndCATrackParamVector::Cov(), fC, fP, i, and PndCATrackParamVector::Par().

    : fX( v.X()[i] ),
      fSignCosPhi( v.SignCosPhi()[i] ),
      fChi2( v.Chi2()[i] ),
      fNDF( v.NDF()[i] ),
      fAlpha( v.Angle()[i] )
  {
    for ( int j = 0; j <  5; ++j ) fP[j] = v.Par()[j][i];
    for ( int j = 0; j < 15; ++j ) fC[j] = v.Cov()[j][i];
  }

Member Function Documentation

float PndFTSCATrackParam::Angle ( ) const

inline

Definition at line 289 of file PndFTSCATrackParam.h.

References fAlpha.

Referenced by PndFtsCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawGBHits(), PndFTSCADisplay::DrawGBTrack(), PndFTSCADisplay::DrawRecoTrack(), PndFTSTopoReconstructor::GetChiToPrimVertex(), GetXYZPxPyPzQ(), and PndFTSTopoReconstructor::Init().

{ return fAlpha; }

float PndFTSCATrackParam::Chi2 ( ) const

inline

Definition at line 277 of file PndFTSCATrackParam.h.

References fChi2.

Referenced by PndFTSCAGBTracker::CreateTracks(), FTSCATrack::Fit(), FTSCATrack::Fit2Times(), PndFTSTopoReconstructor::Init(), and operator<<().

{ return fChi2; }

float PndFTSCATrackParam::CosPhi ( ) const

inline

Definition at line 292 of file PndFTSCATrackParam.h.

References fSignCosPhi, SinPhi(), and CAMath::Sqrt().

Referenced by PndFTSCADisplay::DrawGBTrackFast(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), Rotate(), S(), and TransportToX0().

{ return fSignCosPhi*CAMath::Sqrt( 1 - SinPhi()*SinPhi() ); }

const float* PndFTSCATrackParam::Cov ( ) const

inline

Definition at line 297 of file PndFTSCATrackParam.h.

References fC.

Referenced by PndFtsCATracking::CATrackParToFairTrackParP(), PndFTSTopoReconstructor::GetChiToPrimVertex(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), and operator<<().

{ return fC; }

float PndFTSCATrackParam::DzDs ( ) const

inline

Definition at line 271 of file PndFTSCATrackParam.h.

References fP.

Referenced by PndFTSCADisplay::DrawRecoTrack(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), QP(), and TransportToX0().

{ return fP[3]; }

float PndFTSCATrackParam::Err2DzDs ( ) const

inline

Definition at line 286 of file PndFTSCATrackParam.h.

References fC.

{ return fC[9]; }

float PndFTSCATrackParam::Err2QMomentum ( ) const

inline

Definition at line 294 of file PndFTSCATrackParam.h.

References fC.

Referenced by FTSCANPlet::QMomentumErr(), and FTSCANPlet::QMomentumErr2().

{ return fC[14]; }

float PndFTSCATrackParam::Err2QPt ( ) const

inline

Definition at line 287 of file PndFTSCATrackParam.h.

References fC.

{ return fC[14]; }

float PndFTSCATrackParam::Err2SinPhi ( ) const

inline

Definition at line 285 of file PndFTSCATrackParam.h.

References fC.

{ return fC[5]; }

float PndFTSCATrackParam::Err2X1 ( ) const

inline

Definition at line 280 of file PndFTSCATrackParam.h.

References Err2Y().

{ return Err2Y(); }

float PndFTSCATrackParam::Err2X2 ( ) const

inline

Definition at line 281 of file PndFTSCATrackParam.h.

References Err2Z().

{ return Err2Z(); }

float PndFTSCATrackParam::Err2Y ( ) const

inline

Definition at line 283 of file PndFTSCATrackParam.h.

References fC.

Referenced by Err2X1().

{ return fC[0]; }

float PndFTSCATrackParam::Err2Z ( ) const

inline

Definition at line 284 of file PndFTSCATrackParam.h.

References fC.

Referenced by Err2X2().

{ return fC[2]; }

void PndFTSCATrackParam::GetDCAPoint	(	float	x,
		float	y,
		float	z,
		float &	px,
		float &	py,
		float &	pz,
		float	Bz
	)		const

Definition at line 92 of file PndFTSCATrackParam.cxx.

References a, CAMath::Abs(), CosPhi(), dx, dy, DzDs(), Kappa(), CAMath::Nint(), s, S(), SinPhi(), sqrt(), CAMath::TwoPi(), X(), x0, Y(), y0, and Z().

Referenced by PndFTSCADisplay::DrawGBTrackFast().

{
  //* Get the track point closest to the (x,y,z)

  float x0 = X();
  float y0 = Y();
  float k  = Kappa( Bz );
  float ex = CosPhi();
  float ey = SinPhi();
  float dx = x - x0;
  float dy = y - y0;
  float ax = dx * k + ey;
  float ay = dy * k - ex;
  float a = sqrt( ax * ax + ay * ay );
  xp = x0 + ( dx - ey * ( ( dx * dx + dy * dy ) * k - 2 * ( -dx * ey + dy * ex ) ) / ( a + 1 ) ) / a;
  yp = y0 + ( dy + ex * ( ( dx * dx + dy * dy ) * k - 2 * ( -dx * ey + dy * ex ) ) / ( a + 1 ) ) / a;
  float s = S( x, y, Bz );
  zp = Z() + DzDs() * s;
  if ( CAMath::Abs( k ) > 1.e-2 ) {
    float dZ = CAMath::Abs( DzDs() * CAMath::TwoPi() / k );
    if ( dZ > .1 ) {
      zp += CAMath::Nint( ( z - zp ) / dZ ) * dZ;
    }
  }
}

bool PndFTSCATrackParam::GetXYZPxPyPzQ	(	float &	x,
		float &	y,
		float &	z,
		float &	px,
		float &	py,
		float &	pz,
		int &	q,
		float	cov[21]
	)		const

Definition at line 198 of file PndFTSCATrackParam.cxx.

References CAMath::Abs(), alpha, Angle(), ASSERT, cos(), CosPhi(), Cov(), DzDs(), f, fabs(), i, ok, pt(), QPt(), sin(), SinPhi(), x, X(), y, Y(), and Z().

{
  {  // -- convert parameters
  
    x = X();
    y = Y();
    z = Z();
    
    const float pt = CAMath::Abs( 1.f / QPt() );
    q = -( QPt() >=0 ? 1 : -1 );
    
    const float cosL = DzDs();
    px = pt * CosPhi();
    py = pt * SinPhi();
    pz = pt * cosL;

      // -- convert cov matrix
      // get jacobian
    float J[6][6];
    for (int i = 0; i < 6; i++)
      for (int j = 0; j < 6; j++)
        J[i][j] = 0;
    J[0][0] = 1; // x -> x
    J[1][1] = 1; // y -> y
    J[2][2] = 1; // z -> z
    J[3][3] = pt * CosPhi() / SinPhi();
    J[3][5] = - q * pt * pt * CosPhi(); // q/pt -> px
    J[4][3] = pt; // sinPhi -> py
    J[4][5] = - q* pt * pt * SinPhi(); // q/pt -> py
    J[5][4] = pt; // dz/ds -> pz
    J[5][5] = - q* pt * pt * cosL; // q/pt -> pz

    float CovIn[6][6]; // triangular -> symmetric matrix
    {
      CovIn[0][0] = .000001f*.00001f; // dx. From nowhere. TODO
      for (int i = 1; i < 6; i++) {
        CovIn[i][0] = 0;
        CovIn[0][i] = 0;
      }
      int k = 0;
      for (int i = 1; i < 6; i++) {
        for (int j = 1; j <= i; j++, k++) {
          CovIn[i][j] = Cov()[k];
          CovIn[j][i] = Cov()[k];
        }
      }
    }
    
    float CovInJ[6][6];      // CovInJ = CovIn * J^t
    for (int i = 0; i < 6; i++)
      for (int j = 0; j < 6; j++) {
        CovInJ[i][j] = 0;
        for (int k = 0; k < 6; k++) {
          CovInJ[i][j] += CovIn[i][k] * J[j][k];
        }
      }
    
    float CovOut[6][6];      // CovOut = J * CovInJ
    for (int i = 0; i < 6; i++)
      for (int j = 0; j < 6; j++) {
        CovOut[i][j] = 0;
        for (int k = 0; k < 6; k++) {
          CovOut[i][j] += J[i][k] * CovInJ[k][j];
        }
      }

      // symmetric matrix -> triangular
    {
      int k = 0;
      for (int i = 0; i < 6; i++) {
        for (int j = 0; j <= i; j++, k++) {
          cov[k] = CovOut[i][j];
          ASSERT( !finite(CovOut[i][j]) ||  CovOut[i][j] == 0 || fabs( 1. - CovOut[j][i]/CovOut[i][j] ) <= 0.05,
          "CovOut[" << i << "][" << j << "] == CovOut[" << j << "][" << i << "] : " << CovOut[i][j] << " == " << CovOut[j][i]);
        }
      }
    }
     
    {   // check cov matrix
      bool ok = true;
      int k = 0;
      for (int i = 0; i < 6; i++) {
        for (int j = 0; j <= i; j++, k++) {
          ok &= finite( cov[k] );
        }
        ok &= ( cov[k-1] > 0 );
      }
      if (!ok) return false;
    }
  }

    // Rotate in XY plane to convert to global
  { 
    const float alpha = Angle();

#if defined(PANDA_STT) || defined(PANDA_FTS)
    assert(0); // different coordinate transform for panda
#endif
    const float cA = cos( alpha );
    const float sA = sin( alpha );

      //float J[6][6] = { {  cA, sA, 0,  0,  0,  0 }, // X
      //                  { -sA, cA, 0,  0,  0,  0 }, // Y
      //                  {  0,  0,  1,  0,  0,  0 }, // Z
      //                  {  0,  0,  0,  cA, sA, 0 }, // Px
      //                  {  0,  0,  0, -sA, cA, 0 }, // Py
      //                  {  0,  0,  0,  0,  0,  1 } }; // Pz
  
    float J[2][2] = { {  cA, sA },
                      { -sA, cA } };

    { // convert x, y TODO optimize
      const float xt  = x,  yt = y;
      x = xt*cA  +  yt*sA;
      y = -xt*sA +  yt*cA;

      float Cov1[2][2]; // triangular -> symmetric matrix
      Cov1[0][0] = cov[0];
      Cov1[0][1] = cov[1];
      Cov1[1][0] = cov[1];
      Cov1[1][1] = cov[2];
    
      float Cov2[2][2];      // Cov2 = Cov1 * J^t
      for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++) {
          Cov2[i][j] = 0;
          for (int k = 0; k < 2; k++) {
            Cov2[i][j] += Cov1[i][k] * J[j][k];
          }
        }
    
        // Cov1 = J * Cov2
      for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++) {
          Cov1[i][j] = 0;
          for (int k = 0; k < 2; k++) {
            Cov1[i][j] += J[i][k] * Cov2[k][j];
          }
        }

        // symmetric matrix -> triangular
      cov[0] = Cov1[0][0];
      cov[1] = Cov1[0][1];
      cov[2] = Cov1[1][1];
    }

    { // convert Px, Py TODO optimize
      const float xt  = px,  yt = py;
      px = xt*cA  +  yt*sA;
      py = -xt*sA +  yt*cA;

      float Cov1[2][2]; // triangular -> symmetric matrix
      Cov1[0][0] = cov[9];
      Cov1[0][1] = cov[13];
      Cov1[1][0] = cov[13];
      Cov1[1][1] = cov[14];

    
      float Cov2[2][2];      // Cov2 = Cov1 * J^t
      for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++) {
          Cov2[i][j] = 0;
          for (int k = 0; k < 2; k++) {
            Cov2[i][j] += Cov1[i][k] * J[j][k];
          }
        }
    
        // Cov1 = J * Cov2
      for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++) {
          Cov1[i][j] = 0;
          for (int k = 0; k < 2; k++) {
            Cov1[i][j] += J[i][k] * Cov2[k][j];
          }
        }

        // symmetric matrix -> triangular

      cov[9]  = Cov1[0][0];
      cov[13] = Cov1[0][1];
      cov[14] = Cov1[1][1];
    }
  }

  return true;
}

bool PndFTSCATrackParam::IsValid ( ) const

inline

Definition at line 313 of file PndFTSCATrackParam.h.

References fChi2.

{ return fChi2 != -1; }

float PndFTSCATrackParam::Kappa ( float  Bz ) const

inline

Definition at line 291 of file PndFTSCATrackParam.h.

References fP.

Referenced by PndFTSCADisplay::DrawGBTrackFast(), GetDCAPoint(), and S().

{ return fP[4]*Bz; }

int PndFTSCATrackParam::NDF ( ) const

inline

Definition at line 278 of file PndFTSCATrackParam.h.

References fNDF.

Referenced by FTSCATrack::Fit(), FTSCATrack::Fit2Times(), PndFTSTopoReconstructor::Init(), and operator<<().

{ return fNDF; }

const float* PndFTSCATrackParam::Par ( ) const

inline

Definition at line 296 of file PndFTSCATrackParam.h.

References fP.

Referenced by operator<<().

{ return fP; }

float PndFTSCATrackParam::QMomentum ( ) const

inline

Definition at line 274 of file PndFTSCATrackParam.h.

References QPt().

Referenced by FTSCANPlet::QMomentum().

{ return QPt(); }

float PndFTSCATrackParam::QP ( ) const

inline

Definition at line 273 of file PndFTSCATrackParam.h.

References DzDs(), QPt(), and sqrt().

Referenced by PndFtsCATracking::CATrackParToFairTrackParP(), FTSCATrack::Fit(), and PndFTSCAGBTracker::Merge().

{ return QPt() / sqrt(DzDs()*DzDs() + 1); }

float PndFTSCATrackParam::QPt ( ) const

inline

Definition at line 272 of file PndFTSCATrackParam.h.

References fP.

Referenced by PndFTSCADisplay::DrawGBTrackFast(), PndFTSCADisplay::DrawRecoTrack(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), QMomentum(), QP(), and TransportToX0().

{ return fP[4]; }

void PndFTSCATrackParam::Reset ( )

inlineprivate

Definition at line 317 of file PndFTSCATrackParam.h.

References fC, fChi2, fNDF, fP, fSignCosPhi, and fX.

Referenced by PndFTSCATrackParam().

               {
    fX = 0;
    fSignCosPhi = 0;
    for(int i=0; i < 5; i++) fP[i] = 0;
    for(int i=0; i < 15; i++) fC[i] = 0;
    fChi2 = 0;
    fNDF = 0;
  }

bool PndFTSCATrackParam::Rotate	(	float	alpha,
		float	maxSinPhi = .999
	)

Definition at line 120 of file PndFTSCATrackParam.cxx.

References CAMath::Abs(), alpha, CAMath::Cos(), CosPhi(), fabs(), fAlpha, fC, fP, fSignCosPhi, fX, CAMath::Sin(), SinPhi(), x, X(), y, and Y().

Referenced by PndFTSCADisplay::DrawGBHits(), PndFTSCADisplay::DrawTrack(), and Transport().

{
  //* Rotate the coordinate system in XY on the angle alpha

  const float cA = CAMath::Cos( alpha );
  const float sA = CAMath::Sin( alpha );
  const float x = X(), y = Y(), sP = SinPhi(), cP = CosPhi();
  const float cosPhi = cP * cA + sP * sA;
  const float sinPhi = -cP * sA + sP * cA;

  if ( CAMath::Abs( sinPhi ) > maxSinPhi || CAMath::Abs( cosPhi ) < 1.e-2 || CAMath::Abs( cP ) < 1.e-2  ) return 0;

  if ( fabs(alpha) > 3.1415 * 0.25 ) return 0; // allow turn by 45 degree only
    
  const float j0 = cP / cosPhi;
  const float j2 = cosPhi / cP;

  fX = x*cA +  y*sA;
  fP[0] = -x*sA +  y*cA;
  fSignCosPhi = CAMath::Abs(cosPhi) / cosPhi;
  fP[2] = sinPhi;

  //  float J[5][5] = { { j0, 0, 0,  0,  0 }, // Y
  //                    {  0, 1, 0,  0,  0 }, // Z
  //                    {  0, 0, j2, 0,  0 }, // SinPhi
  //                    {  0, 0, 0,  1,  0 }, // DzDs
  //                    {  0, 0, 0,  0,  1 } }; // Kappa
  fC[0] *= j0 * j0;
  fC[1] *= j0;
  fC[3] *= j0;
  fC[6] *= j0;
  fC[10] *= j0;

  fC[3] *= j2;
  fC[4] *= j2;
  fC[5] *= j2 * j2;
  fC[8] *= j2;
  fC[12] *= j2;

  fAlpha += alpha;
  
  return 1;
}

float PndFTSCATrackParam::S ( float  x, float  y, float  Bz  ) const

private

Definition at line 78 of file PndFTSCATrackParam.cxx.

References CAMath::Abs(), CAMath::ATan2(), CosPhi(), Kappa(), SinPhi(), X(), and Y().

Referenced by GetDCAPoint().

{
  //* Get XY path length to the given point

  float k  = Kappa( Bz );
  float ex = CosPhi();
  float ey = SinPhi();
  x -= X();
  y -= Y();
  float dS = x * ex + y * ey;
  if ( CAMath::Abs( k ) > 1.e-4 ) dS = CAMath::ATan2( k * dS, 1 + k * ( x * ey - y * ex ) ) / k;
  return dS;
}

void PndFTSCATrackParam::SetAsInvalid ( )

inline

Definition at line 314 of file PndFTSCATrackParam.h.

References fChi2.

Referenced by PndFTSCAGBTracker::FitTracks().

{ fChi2 = -1; }

void PndFTSCATrackParam::SetSinPhi ( float  v )

inline

Definition at line 301 of file PndFTSCATrackParam.h.

References fP, and v.

Referenced by PndFTSCADisplay::DrawTrack().

{  fP[2] = v; }

float PndFTSCATrackParam::SignCosPhi ( ) const

inline

Definition at line 276 of file PndFTSCATrackParam.h.

References fSignCosPhi.

Referenced by PndFTSCADisplay::DrawRecoTrack(), and operator<<().

{ return fSignCosPhi; }

float PndFTSCATrackParam::SinPhi ( ) const

inline

Definition at line 270 of file PndFTSCATrackParam.h.

References fP.

Referenced by CosPhi(), PndFTSCADisplay::DrawGBTrackFast(), PndFTSCADisplay::DrawRecoTrack(), PndFTSCADisplay::DrawTrack(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), Rotate(), S(), and TransportToX0().

{ return fP[2]; }

bool PndFTSCATrackParam::Transport	(	const FTSCAHit &	hit,
		const PndFTSCAParam &	param
	)

Definition at line 165 of file PndFTSCATrackParam.cxx.

References FTSCAHit::Angle(), PndFTSCAParam::cBz(), fAlpha, Rotate(), TransportToX0(), and FTSCAHit::X0().

Referenced by PndFTSCAGBTracker::Merge().

{
  const bool rotated = Rotate( -fAlpha + hit.Angle() );
  const bool transported = TransportToX0( hit.X0(), param.cBz() );
  return rotated & transported;
}

bool PndFTSCATrackParam::TransportToX0 ( float  x, float  Bz, float  maxSinPhi = .999  )

inline

Definition at line 337 of file PndFTSCATrackParam.h.

References CAMath::Abs(), CAMath::ASin(), c10, c11, c20, c21, c22, CosPhi(), dx, dy, dz, DzDs(), f, fC, fP, fX, h2, h4, QPt(), SinPhi(), CAMath::Sqrt(), X(), Y(), and Z().

Referenced by PndFTSCADisplay::DrawGBHits(), PndFTSTopoReconstructor::GetChiToPrimVertex(), PndFTSTopoReconstructor::Init(), and Transport().

{
    //* Transport the track parameters to X=x, using linearization at t0, and the field value Bz
    //* maxSinPhi is the max. allowed value for |t0.SinPhi()|
    //* linearisation of trajectory t0 is also transported to X=x,
    //* returns 1 if OK
    //*

  const float ex = CosPhi();
  const float ey = SinPhi();
  const float k  = QPt() * Bz;
  const float dx = x - X();

  const float ey1 = k * dx + ey;

    // check for intersection with X=x

  if ( CAMath::Abs( ey1 ) > maxSinPhi ) return 0;

  float ex1 = CAMath::Sqrt( 1.f - ey1 * ey1 );
  if ( ex < 0 ) ex1 = -ex1;

  const float dx2 = dx * dx;
  const float ss = ey + ey1;
  const float cc = ex + ex1;

  if ( ( CAMath::Abs( cc ) < 1.e-4 || CAMath::Abs( ex ) < 1.e-4 || CAMath::Abs( ex1 ) < 1.e-4 ) ) return 0;
  
  const float cci = 1.f / cc;
  const float exi = 1.f / ex;
  const float ex1i = 1.f / ex1;
  
  const float tg = ss * cci; // tan((phi1+phi)/2)

  const float dy = dx * tg;
  float dl = dx * CAMath::Sqrt( 1.f + tg * tg );

  if ( cc < 0 ) dl = -dl;
  float dSin = dl * k * 0.5;
  if ( dSin > 1.f ) dSin = 1.f;
  if ( dSin < -1.f ) dSin = -1.f;
  const float dS = ( CAMath::Abs( k ) > 1.e-4 )  ? ( 2 * CAMath::ASin( dSin ) / k ) : dl;
  const float dz = dS * DzDs();


    //float H0[5] = { 1,0, h2,  0, h4 };
    //float H1[5] = { 0, 1, 0, dS,  0 };
    //float H2[5] = { 0, 0, 1,  0, dxBz };
    //float H3[5] = { 0, 0, 0,  1,  0 };
    //float H4[5] = { 0, 0, 0,  0,  1 };

  const float h2 = dx * ( 1.f + ey * ey1 + ex * ex1 ) * exi * ex1i * cci;
  const float h4 = dx2 * ( cc + ss * ey1 * ex1i ) * cci * cci * Bz;
  const float dxBz = dx * Bz;

  fX    = X() + dx;
  fP[0] = Y() + dy;
  fP[1] = Z() + dz;
  fP[2] = ey1;

#if 1
  const float c00 = fC[0];
  const float c10 = fC[1];
  const float c11 = fC[2];
  const float c20 = fC[3];
  const float c21 = fC[4];
  const float c22 = fC[5];
  const float c30 = fC[6];
  const float c31 = fC[7];
  const float c32 = fC[8];
  const float c33 = fC[9];
  const float c40 = fC[10];
  const float c41 = fC[11];
  const float c42 = fC[12];
  const float c43 = fC[13];
  const float c44 = fC[14];

  fC[0] = ( c00  + h2 * h2 * c22 + h4 * h4 * c44
          + 2.f * ( h2 * c20 + h4 * c40 + h2 * h4 * c42 )  );

  fC[1] = c10 + h2 * c21 + h4 * c41 + dS * ( c30 + h2 * c32 + h4 * c43 );
  fC[2] = c11 + 2.f * dS * c31 + dS * dS * c33;

  fC[3] = c20 + h2 * c22 + h4 * c42 + dxBz * ( c40 + h2 * c42 + h4 * c44 );
  fC[4] = c21 + dS * c32 + dxBz * ( c41 + dS * c43 );
  fC[5] = c22 + 2.f * dxBz * c42 + dxBz * dxBz * c44;

  fC[6] = c30 + h2 * c32 + h4 * c43;
  fC[7] = c31 + dS * c33;
  fC[8] = c32 + dxBz * c43;
  fC[9] = c33;

  fC[10] = c40 + h2 * c42 + h4 * c44;
  fC[11] = c41 + dS * c43;
  fC[12] = c42 + dxBz * c44;
  fC[13] = c43;
  fC[14] = c44;
#else
  fC[0] = ( fC[0]  + h2 * h2 * fC[5] + h4 * h4 * fC[14]
          + 2 * ( h2 * fC[3] + h4 * fC[10] + h2 * h4 * fC[12] )  );

  fC[1] = fC[1] + h2 * fC[4] + h4 * fC[11] + dS * ( fC[6] + h2 * fC[8] + h4 * fC[13] );
  fC[2] = fC[2] + 2 * dS * fC[7] + dS * dS * fC[9];

  fC[3] = fC[3] + h2 * fC[5] + h4 * fC[12] + dxBz * ( fC[10] + h2 * fC[12] + h4 * fC[14] );
  fC[4] = fC[4] + dS * fC[8] + dxBz * ( fC[11] + dS * fC[13] );
  fC[5] = fC[5] + 2 * dxBz * fC[12] + dxBz * dxBz * fC[14];

  fC[6] = fC[6] + h2 * fC[8] + h4 * fC[13];
  fC[7] = fC[7] + dS * fC[9];
  fC[8] = fC[8] + dxBz * fC[13];
  fC[9] = fC[9];

  fC[10] = fC[10] + h2 * fC[12] + h4 * fC[14];
  fC[11] = fC[11] + dS * fC[13];
  fC[12] = fC[12] + dxBz * fC[14];
  fC[13] = fC[13];
  fC[14] = fC[14];
#endif

  return 1;
}

float PndFTSCATrackParam::X ( ) const

inline

Definition at line 267 of file PndFTSCATrackParam.h.

References fX.

Referenced by PndFtsCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawGBTrackFast(), PndFTSCADisplay::DrawRecoTrack(), PndFTSCADisplay::DrawTrack(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), PndFTSCAGBTracker::Merge(), operator<<(), Rotate(), S(), and TransportToX0().

{ return fX;    }

float PndFTSCATrackParam::X0 ( ) const

inline

Definition at line 263 of file PndFTSCATrackParam.h.

References fX.

{ return fX; }

float PndFTSCATrackParam::X1 ( ) const

inline

Definition at line 264 of file PndFTSCATrackParam.h.

References Y().

{ return Y(); }

float PndFTSCATrackParam::X2 ( ) const

inline

Definition at line 265 of file PndFTSCATrackParam.h.

References Z().

{ return Z(); }

float PndFTSCATrackParam::Y ( ) const

inline

Definition at line 268 of file PndFTSCATrackParam.h.

References fP.

Referenced by PndFtsCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawGBHits(), PndFTSCADisplay::DrawGBTrackFast(), PndFTSCADisplay::DrawRecoTrack(), PndFTSCADisplay::DrawTrack(), PndFTSTopoReconstructor::GetChiToPrimVertex(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), PndFTSCAGBTracker::Merge(), Rotate(), S(), TransportToX0(), and X1().

{ return fP[0]; }

float PndFTSCATrackParam::Z ( ) const

inline

Definition at line 269 of file PndFTSCATrackParam.h.

References fP.

Referenced by PndFTSCADisplay::DrawGBHits(), PndFTSCADisplay::DrawRecoTrack(), PndFTSCADisplay::DrawTrack(), PndFTSTopoReconstructor::GetChiToPrimVertex(), GetDCAPoint(), GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), PndFTSCAGBTracker::Merge(), TransportToX0(), and X2().

{ return fP[1]; }

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  , const PndFTSCATrackParam &    )

friend

Definition at line 184 of file PndFTSCATrackParam.cxx.

{
  out << t.X() << " "
  << t.SignCosPhi() << " "
  << t.Chi2() << " "
  << t.NDF()
  << '\n';
  for ( int i = 0; i < 5; i++ ) out << t.Par()[i] << " ";
  out << '\n';
  for ( int i = 0; i < 15; i++ ) out << t.Cov()[i] << " ";
  out << '\n';
  return out;
}

std::istream& operator>> ( std::istream &  , PndFTSCATrackParam &    )

friend

Definition at line 173 of file PndFTSCATrackParam.cxx.

{
  in >> t.fX;
  in >> t.fSignCosPhi;
  in >> t.fChi2;
  in >> t.fNDF;
  for ( int i = 0; i < 5; i++ ) in >> t.fP[i];
  for ( int i = 0; i < 15; i++ ) in >> t.fC[i];
  return in;
}

Member Data Documentation

float PndFTSCATrackParam::fAlpha

private

Definition at line 334 of file PndFTSCATrackParam.h.

Referenced by Angle(), Rotate(), and Transport().

float PndFTSCATrackParam::fC[15]

private

Definition at line 330 of file PndFTSCATrackParam.h.

Referenced by Cov(), Err2DzDs(), Err2QMomentum(), Err2QPt(), Err2SinPhi(), Err2Y(), Err2Z(), operator>>(), PndFTSCATrackParam(), Reset(), Rotate(), and TransportToX0().

float PndFTSCATrackParam::fChi2

private

Definition at line 331 of file PndFTSCATrackParam.h.

Referenced by Chi2(), IsValid(), operator>>(), Reset(), and SetAsInvalid().

int PndFTSCATrackParam::fNDF

private

Definition at line 332 of file PndFTSCATrackParam.h.

Referenced by NDF(), operator>>(), and Reset().

float PndFTSCATrackParam::fP[5]

private

Definition at line 329 of file PndFTSCATrackParam.h.

Referenced by DzDs(), Kappa(), operator>>(), Par(), PndFTSCATrackParam(), QPt(), Reset(), Rotate(), SetSinPhi(), SinPhi(), TransportToX0(), Y(), and Z().

float PndFTSCATrackParam::fSignCosPhi

private

Definition at line 328 of file PndFTSCATrackParam.h.

Referenced by CosPhi(), operator>>(), Reset(), Rotate(), and SignCosPhi().

float PndFTSCATrackParam::fX

private

Definition at line 327 of file PndFTSCATrackParam.h.

Referenced by operator>>(), Reset(), Rotate(), TransportToX0(), X(), and X0().

---

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

• PndFTSCATrackParam.h
• PndFTSCATrackParam.cxx