PndCATrackParam Class Reference

#include <PndCATrackParam.h>

Classes
struct	PndCATrackFitParam

Public Member Functions
	PndCATrackParam ()
	PndCATrackParam (const TrackParamVector &v, int i)
float	X () const
float	Y () const
float	Z () const
float	SinPhi () const
float	DzDs () const
float	QPt () const
float	QMomentum () const
float	SignCosPhi () const
float	Chi2 () const
int	NDF () const
float	Err2Y () const
float	Err2Z () const
float	Err2SinPhi () const
float	Err2DzDs () const
float	Err2QPt () const
float	Angle () const
int	ISec () const
float	GetX () const
float	GetY () const
float	GetZ () const
float	GetSinPhi () const
float	GetDzDs () const
float	GetQPt () const
float	GetSignCosPhi () const
float	GetChi2 () const
int	GetNDF () const
float	GetKappa (float Bz) const
float	GetCosPhiPositive () const
float	GetCosPhi () const
float	GetErr2Y () const
float	GetErr2Z () const
float	GetErr2SinPhi () const
float	GetErr2DzDs () const
float	GetErr2QPt () const
float	Err2X1 () const
float	Err2X2 () const
float	Err2QMomentum () const
const float *	Par () const
const float *	Cov () const
const float *	GetPar () const
const float *	GetCov () const
void	SetPar (int i, float v)
void	SetCov (int i, float v)
void	SetX (float v)
void	SetY (float v)
void	SetZ (float v)
void	SetSinPhi (float v)
void	SetDzDs (float v)
void	SetQPt (float v)
void	SetSignCosPhi (float v)
void	SetChi2 (float v)
void	SetNDF (int v)
void	SetAngle (float v)
void	SetISec (int v)
void	SetErr2QPt (float v)
void	InitDirection (float r0, float r1, float r2)
float	GetDist2 (const PndCATrackParam &t) const
float	GetDistXZ2 (const PndCATrackParam &t) const
float	GetS (float x, float y, float Bz) const
void	GetDCAPoint (float x, float y, float z, float &px, float &py, float &pz, float Bz) const
bool	TransportToX (float x, float Bz, float maxSinPhi=.999)
bool	TransportToXWithMaterial (float x, float Bz, float maxSinPhi=.999)
bool	TransportToX (float x, PndCATrackLinearisation &t0, float Bz, float maxSinPhi=.999, float *DL=0)
bool	TransportToX (float x, float sinPhi0, float cosPhi0, float Bz, float maxSinPhi=.999)
bool	TransportToXWithMaterial (float x, PndCATrackLinearisation &t0, PndCATrackFitParam &par, float Bz, float maxSinPhi=.999)
bool	TransportToXWithMaterial (float x, PndCATrackFitParam &par, float Bz, float maxSinPhi=.999)
void	CalculateFitParameters (PndCATrackFitParam &par, float mass=0.13957)
bool	CorrectForMeanMaterial (float xOverX0, float xTimesRho, const PndCATrackFitParam &par)
bool	Rotate (float alpha, float maxSinPhi=.999)
bool	Rotate (float alpha, PndCATrackLinearisation &t0, float maxSinPhi=.999)
void	RotateXY (float alpha, float &x, float &y, float &sin) const
bool	Filter (float y, float z, float err2Y, float errYZ, float err2Z, float maxSinPhi=.999)
void	Print () const
void	ResetCovMatrix ()
PndCATrackParam	GetGlobalParam (float alpha) const
void	Reset ()
bool	Transport (const PndCAHit &hit, float Bz)
bool	Filter (const PndCAHit &hit)
bool	IsValid () const
void	SetAsInvalid ()

Static Public Member Functions
static float	ApproximateBetheBloch (float beta2)
static float	BetheBlochGeant (float bg, float kp0=2.33, float kp1=0.20, float kp2=3.00, float kp3=173e-9, float kp4=0.49848)
static float	BetheBlochSolid (float bg)
static float	BetheBlochGas (float bg)

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

Detailed Description

PndCATrackParam class describes the track parametrisation which is used by the PndCATracker slice tracker.

Definition at line 28 of file PndCATrackParam.h.

Constructor & Destructor Documentation

PndCATrackParam::PndCATrackParam ( )

inline

Definition at line 34 of file PndCATrackParam.h.

References Reset().

{ Reset();}

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

inline

Definition at line 35 of file PndCATrackParam.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] ),
        fISec( v.ISec()[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 PndCATrackParam::Angle ( ) const

inline

Definition at line 81 of file PndCATrackParam.h.

References fAlpha.

Referenced by PndCATracking::CATrackParToFairTrackParP().

{ return fAlpha; }

float PndCATrackParam::ApproximateBetheBloch ( float  beta2 )

inlinestatic

Definition at line 318 of file PndCATrackParam.h.

References log().

Referenced by CalculateFitParameters().

{
  //------------------------------------------------------------------
  // This is an approximation of the Bethe-Bloch formula with
  // the density effect taken into account at beta*gamma > 3.5
  // (the approximation is reasonable only for solid materials)
  //------------------------------------------------------------------
  if ( beta2 >= 1 )
    return 0;
  else {
    const float beta2_beta21i = beta2 / ( 1 - beta2 );
    if ( beta2_beta21i > 12.25 ) // 3.5^2 = 12.25
      return 0.153e-3 / beta2 * ( 9.94223 + 0.5 * log( beta2_beta21i ) - beta2 ); // log( 3.5*5940 ) = 9.94223
  else
      return 0.153e-3 / beta2 * ( 8.6895 + log( beta2_beta21i ) - beta2 ); // log( 5940 ) = 8.6895
  }
}

float PndCATrackParam::BetheBlochGas ( float  bg )

static

Definition at line 262 of file PndCATrackParam.cxx.

References BetheBlochGeant(), and x0.

{
  //------------------------------------------------------------------
  // This is an approximation of the Bethe-Bloch formula,
  // reasonable for gas materials.
  // All the parameters are, in fact, for Ne.
  // The returned value is in [GeV]
  //------------------------------------------------------------------

  const float rho = 0.9e-3;
  const float x0  = 2.;
  const float x1  = 4.;
  const float mI  = 140.e-9;
  const float mZA = 0.49555;

  return BetheBlochGeant( bg, rho, x0, x1, mI, mZA );
}

float PndCATrackParam::BetheBlochGeant ( float  bg, float  kp0 = 2.33, float  kp1 = 0.20, float  kp2 = 3.00, float  kp3 = 173e-9, float  kp4 = 0.49848  )

static

Definition at line 206 of file PndCATrackParam.cxx.

References CAMath::Log(), me, mK, r, CAMath::Sqrt(), x, and x0.

Referenced by BetheBlochGas(), and BetheBlochSolid().

{
  //
  // This is the parameterization of the Bethe-Bloch formula inspired by Geant.
  //
  // bg2  - (beta*gamma)^2
  // kp0 - density [g/cm^3]
  // kp1 - density effect first junction point
  // kp2 - density effect second junction point
  // kp3 - mean excitation energy [GeV]
  // kp4 - mean Z/A
  //
  // The default values for the kp* parameters are for silicon.
  // The returned value is in [GeV/(g/cm^2)].
  //

  const float mK  = 0.307075e-3; // [GeV*cm^2/g]
  const float me  = 0.511e-3;    // [GeV/c^2]
  const float rho = kp0;
  const float x0  = kp1 * 2.303;
  const float x1  = kp2 * 2.303;
  const float mI  = kp3;
  const float mZA = kp4;
  const float maxT = 2 * me * bg2;    // neglecting the electron mass

  //*** Density effect
  float d2 = 0.;
  const float x = 0.5 * CAMath::Log( bg2 );
  const float lhwI = CAMath::Log( 28.816 * 1e-9 * CAMath::Sqrt( rho * mZA ) / mI );
  if ( x > x1 ) {
    d2 = lhwI + x - 0.5;
  } else if ( x > x0 ) {
    const float r = ( x1 - x ) / ( x1 - x0 );
    d2 = lhwI + x - 0.5 + ( 0.5 - lhwI - x0 ) * r * r * r;
  }

  return mK*mZA*( 1 + bg2 ) / bg2*( 0.5*CAMath::Log( 2*me*bg2*maxT / ( mI*mI ) ) - bg2 / ( 1 + bg2 ) - d2 );
}

float PndCATrackParam::BetheBlochSolid ( float  bg )

static

Definition at line 250 of file PndCATrackParam.cxx.

References BetheBlochGeant().

{
  //------------------------------------------------------------------
  // This is an approximation of the Bethe-Bloch formula,
  // reasonable for solid materials.
  // All the parameters are, in fact, for Si.
  // The returned value is in [GeV]
  //------------------------------------------------------------------

  return BetheBlochGeant( bg );
}

void PndCATrackParam::CalculateFitParameters ( PndCATrackFitParam &  par, float  mass = 0.13957  )

inline

Definition at line 336 of file PndCATrackParam.h.

References ApproximateBetheBloch(), PndCATrackParam::PndCATrackFitParam::fBethe, PndCATrackParam::PndCATrackFitParam::fE, PndCATrackParam::PndCATrackFitParam::fEP2, PndCATrackParam::PndCATrackFitParam::fK22, PndCATrackParam::PndCATrackFitParam::fK33, PndCATrackParam::PndCATrackFitParam::fK43, PndCATrackParam::PndCATrackFitParam::fK44, fP, PndCATrackParam::PndCATrackFitParam::fSigmadE2, PndCATrackParam::PndCATrackFitParam::fTheta2, p2, and CAMath::Sqrt().

Referenced by Transport(), and TransportToXWithMaterial().

{
  const float p2 = ( 1. + fP[3] * fP[3] );
  const float k2 = fP[4] * fP[4];
  const float mass2 = mass * mass;
  
  const float beta2 = p2 / ( p2 + mass2 * k2 );
        
  const float pp2 = ( k2 > 1.e-8 ) ? p2 / k2 : 10000; // impuls 2
        
    //par.fBethe = BetheBlochGas( pp2/mass2);
  par.fBethe = ApproximateBetheBloch( pp2 / mass2 );
  par.fE = CAMath::Sqrt( pp2 + mass2 );
  par.fTheta2 = 198.81e-6  / ( beta2 * pp2 ); // 14.1^2 * 1e-6 
  par.fEP2 = par.fE / pp2;                    // have tried reduce number of "/", but it was slower. (may be bacause of additional of constants = memory)

  // Approximate energy loss fluctuation (M.Ivanov)

  const float knst = 0.07; // To be tuned.
  par.fSigmadE2 = knst * par.fEP2 * fP[4];
  par.fSigmadE2 = par.fSigmadE2 * par.fSigmadE2;

  par.fK22 = p2;
  par.fK33 = par.fK22 * par.fK22;
  par.fK43 = fP[3] * fP[4] * par.fK22;
  par.fK44 = (p2 - 1.f) * k2;

}

float PndCATrackParam::Chi2 ( ) const

inline

Definition at line 72 of file PndCATrackParam.h.

References fChi2.

Referenced by PndCAGBTracker::CreateTracks(), and PndCATrack::Fit().

{ return fChi2; }

bool PndCATrackParam::CorrectForMeanMaterial ( float  xOverX0, float  xTimesRho, const PndCATrackFitParam &  par  )

inline

Definition at line 501 of file PndCATrackParam.h.

References CAMath::Abs(), corr, dE, PndCATrackParam::PndCATrackFitParam::fBethe, fC, PndCATrackParam::PndCATrackFitParam::fE, PndCATrackParam::PndCATrackFitParam::fEP2, PndCATrackParam::PndCATrackFitParam::fK22, PndCATrackParam::PndCATrackFitParam::fK33, PndCATrackParam::PndCATrackFitParam::fK43, PndCATrackParam::PndCATrackFitParam::fK44, fP, PndCATrackParam::PndCATrackFitParam::fSigmadE2, and PndCATrackParam::PndCATrackFitParam::fTheta2.

{
  //------------------------------------------------------------------
  // This function corrects the track parameters for the crossed material.
  // "xOverX0"   - X/X0, the thickness in units of the radiation length.
  // "xTimesRho" - is the product length*density (g/cm^2).
  //------------------------------------------------------------------
  // float &fC22 = fC[5];
  // float &fC33 = fC[9];
  // float &fC40 = fC[10];
  // float &fC41 = fC[11];
  // float &fC42 = fC[12];
  // float &fC43 = fC[13];
  // float &fC44 = fC[14];

  //Energy losses************************

  const float dE = par.fBethe * xTimesRho;
  if ( CAMath::Abs( dE ) > 0.3 * par.fE ) return 0; //30% energy loss is too much!
  const float corr = ( 1. - par.fEP2 * dE );
  if ( corr < 0.3 || corr > 1.3 ) return 0;

  fP[4] *= corr;
  fC[10] *= corr;
  fC[11] *= corr;
  fC[12] *= corr;
  fC[13] *= corr;
  fC[14] *= corr * corr;
  fC[14] += par.fSigmadE2 * CAMath::Abs( dE );
//  std::cout << "dE "<<dE<<"    corr "<<corr<<"  fBethe  " <<par.fBethe<<"  XxRo  "<<xTimesRho<<std::endl;

  //Multiple scattering******************

  const float theta2 = par.fTheta2 * CAMath::Abs( xOverX0 );
  fC[5] += theta2 * par.fK22 * ( 1. - fP[2] * fP[2] );
  fC[9] += theta2 * par.fK33;
  fC[13] += theta2 * par.fK43;
  fC[14] += theta2 * par.fK44;
  
  return 1;
}

const float* PndCATrackParam::Cov ( ) const

inline

Definition at line 109 of file PndCATrackParam.h.

References fC.

Referenced by PndCATracking::CATrackParToFairTrackParP(), and PndCAGBTracker::Merge().

{ return fC; }

float PndCATrackParam::DzDs ( ) const

inline

Definition at line 63 of file PndCATrackParam.h.

References fP.

Referenced by PndCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawTrackParam(), and TransportToX().

{ return fP[3]; }

float PndCATrackParam::Err2DzDs ( ) const

inline

Definition at line 78 of file PndCATrackParam.h.

References fC.

{ return fC[9]; }

float PndCATrackParam::Err2QMomentum ( ) const

inline

Definition at line 106 of file PndCATrackParam.h.

References fC.

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

{ return fC[14]; }

float PndCATrackParam::Err2QPt ( ) const

inline

Definition at line 79 of file PndCATrackParam.h.

References fC.

{ return fC[14]; }

float PndCATrackParam::Err2SinPhi ( ) const

inline

Definition at line 77 of file PndCATrackParam.h.

References fC.

{ return fC[5]; }

float PndCATrackParam::Err2X1 ( ) const

inline

Definition at line 104 of file PndCATrackParam.h.

References fC.

{ return fC[0]; }

float PndCATrackParam::Err2X2 ( ) const

inline

Definition at line 105 of file PndCATrackParam.h.

References fC.

{ return fC[2]; }

float PndCATrackParam::Err2Y ( ) const

inline

Definition at line 75 of file PndCATrackParam.h.

References fC.

{ return fC[0]; }

float PndCATrackParam::Err2Z ( ) const

inline

Definition at line 76 of file PndCATrackParam.h.

References fC.

{ return fC[2]; }

bool PndCATrackParam::Filter ( float  y, float  z, float  err2Y, float  errYZ, float  err2Z, float  maxSinPhi = .999  )

inline

Definition at line 237 of file PndCATrackParam.h.

References CAMath::Abs(), c10, c11, c20, c21, d, f, fC, fChi2, fNDF, fP, ISUNLIKELY, and z0.

Referenced by Filter().

{
  assert( maxSinPhi > 0.f );
  //* Add the y,z measurement with the Kalman filter

  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];
//  float c22 = fC[5];
  const float c30 = fC[6];
  const float c31 = fC[7];
//  float c32 = fC[8];
//  float c33 = fC[9];
  const float c40 = fC[10];
  const float c41 = fC[11];
//  float c42 = fC[12];
//  float c43 = fC[13];
//  float c44 = fC[14];
  
  float d = 1.f / ( err2Y*err2Z + err2Y*c11 + err2Z*c00 + c00*c11 - c10*c10 - 2*errYZ*c10 - errYZ*errYZ );
  err2Y += c00;
  err2Z += c11;
  errYZ += c10;
  
  const float
  z0 = y - fP[0],
  z1 = z - fP[1];

  if ( ISUNLIKELY( err2Y < 1.e-8f ) || ISUNLIKELY( err2Z < 1.e-8f ) ) return 0;

  const float mS0 = err2Z*d;
  const float mS1 = -errYZ*d;
  const float mS2 = err2Y*d;

  // K = CHtS

  const float 
  k00 = c00 * mS0 + c10*mS1,   k01 = c00 * mS1 + c10*mS2,
  k10 = c10 * mS0 + c11*mS1,   k11 = c10 * mS1 + c11*mS2,
  k20 = c20 * mS0 + c21*mS1,   k21 = c20 * mS1 + c21*mS2,
  k30 = c30 * mS0 + c31*mS1,   k31 = c30 * mS1 + c31*mS2,
  k40 = c40 * mS0 + c41*mS1,   k41 = c40 * mS1 + c41*mS2; 

  const float sinPhi = fP[2] + k20 * z0 + k21 * z1;

  if ( ISUNLIKELY( CAMath::Abs( sinPhi ) >= maxSinPhi ) ) return 0;

  fNDF  += 2;
  fChi2 += mS0 * z0 * z0 + mS2 * z1 * z1 + 2 * z0 * z1 * mS1;

  fP[ 0] += k00 * z0 + k01 * z1;
  fP[ 1] += k10 * z0 + k11 * z1;
  fP[ 2] = sinPhi ;
  fP[ 3] += k30 * z0 + k31 * z1;
  fP[ 4] += k40 * z0 + k41 * z1;

  fC[ 0] -= (k00 * c00 + k01 * c10); //c00
  
  fC[ 1] -= (k10 * c00 + k11 * c10); //c10
  fC[ 2] -= (k10 * c10 + k11 * c11); //c11
  
  fC[ 3] -= (k20 * c00 + k21 * c10); //c20
  fC[ 4] -= (k20 * c10 + k21 * c11); //c21
  fC[ 5] -= (k20 * c20 + k21 * c21); //c22
  
  fC[ 6] -= (k30 * c00 + k31 * c10); //c30
  fC[ 7] -= (k30 * c10 + k31 * c11); //c31
  fC[ 8] -= (k30 * c20 + k31 * c21); //c32
  fC[ 9] -= (k30 * c30 + k31 * c31); //c33

  fC[10] -= (k40 * c00 + k41 * c10); //c40
  fC[11] -= (k40 * c10 + k41 * c11); //c41
  fC[12] -= (k40 * c20 + k41 * c21); //c42
  fC[13] -= (k40 * c30 + k41 * c31); //c43
  fC[14] -= (k40 * c40 + k41 * c41); //c44
  
  return 1;
}

bool PndCATrackParam::Filter ( const PndCAHit &  hit )

inline

Definition at line 582 of file PndCATrackParam.h.

References PndCAHit::Err2X1(), PndCAHit::Err2X2(), PndCAHit::ErrX12(), Filter(), PndCAHit::X1(), and PndCAHit::X2().

{
  return Filter( hit.X1(), hit.X2(), hit.Err2X1(), hit.ErrX12(), hit.Err2X2(), 0.999f );
}

float PndCATrackParam::GetChi2 ( ) const

inline

Definition at line 91 of file PndCATrackParam.h.

References fChi2.

{ return fChi2; }

float PndCATrackParam::GetCosPhi ( ) const

inline

Definition at line 96 of file PndCATrackParam.h.

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

Referenced by PndCATracking::CATrackParToFairTrackParP(), GetDCAPoint(), GetS(), Rotate(), and RotateXY().

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

float PndCATrackParam::GetCosPhiPositive ( ) const

inline

Definition at line 95 of file PndCATrackParam.h.

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

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

const float* PndCATrackParam::GetCov ( ) const

inline

Definition at line 112 of file PndCATrackParam.h.

References fC.

{ return fC; }    

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

Definition at line 73 of file PndCATrackParam.cxx.

References a, CAMath::Abs(), dx, dy, GetCosPhi(), GetDzDs(), GetKappa(), GetS(), GetSinPhi(), GetX(), GetY(), GetZ(), CAMath::Nint(), s, sqrt(), CAMath::TwoPi(), x0, and y0.

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

  float x0 = GetX();
  float y0 = GetY();
  float k  = GetKappa( Bz );
  float ex = GetCosPhi();
  float ey = GetSinPhi();
  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 = GetS( x, y, Bz );
  zp = GetZ() + GetDzDs() * s;
  if ( CAMath::Abs( k ) > 1.e-2 ) {
    float dZ = CAMath::Abs( GetDzDs() * CAMath::TwoPi() / k );
    if ( dZ > .1 ) {
      zp += CAMath::Nint( ( z - zp ) / dZ ) * dZ;
    }
  }
}

float PndCATrackParam::GetDist2

(

const PndCATrackParam &

t

)

const

Definition at line 39 of file PndCATrackParam.cxx.

References dx, dy, dz, GetX(), GetY(), and GetZ().

{
  // get squared distance between tracks

  float dx = GetX() - t.GetX();
  float dy = GetY() - t.GetY();
  float dz = GetZ() - t.GetZ();
  return dx*dx + dy*dy + dz*dz;
}

float PndCATrackParam::GetDistXZ2

(

const PndCATrackParam &

t

)

const

Definition at line 49 of file PndCATrackParam.cxx.

References dx, dz, GetX(), and GetZ().

{
  // get squared distance between tracks in X&Z

  float dx = GetX() - t.GetX();
  float dz = GetZ() - t.GetZ();
  return dx*dx + dz*dz;
}

float PndCATrackParam::GetDzDs ( ) const

inline

Definition at line 88 of file PndCATrackParam.h.

References fP.

Referenced by GetDCAPoint(), and Print().

{ return fP[3]; }

float PndCATrackParam::GetErr2DzDs ( ) const

inline

Definition at line 101 of file PndCATrackParam.h.

References fC.

Referenced by Print().

{ return fC[9]; }

float PndCATrackParam::GetErr2QPt ( ) const

inline

Definition at line 102 of file PndCATrackParam.h.

References fC.

Referenced by Print().

{ return fC[14]; }

float PndCATrackParam::GetErr2SinPhi ( ) const

inline

Definition at line 100 of file PndCATrackParam.h.

References fC.

Referenced by Print().

{ return fC[5]; }

float PndCATrackParam::GetErr2Y ( ) const

inline

Definition at line 98 of file PndCATrackParam.h.

References fC.

Referenced by Print().

{ return fC[0]; }

float PndCATrackParam::GetErr2Z ( ) const

inline

Definition at line 99 of file PndCATrackParam.h.

References fC.

Referenced by Print().

{ return fC[2]; }

PndCATrackParam PndCATrackParam::GetGlobalParam

(

float

alpha

)

const

Definition at line 492 of file PndCATrackParam.cxx.

References r, and Rotate().

{
  PndCATrackParam r = *this;
  r.Rotate( -alpha );
  return r;
}

float PndCATrackParam::GetKappa ( float  Bz ) const

inline

Definition at line 94 of file PndCATrackParam.h.

References fP.

Referenced by GetDCAPoint(), and GetS().

{ return fP[4]*Bz; }

int PndCATrackParam::GetNDF ( ) const

inline

Definition at line 92 of file PndCATrackParam.h.

References fNDF.

{ return fNDF; }

const float* PndCATrackParam::GetPar ( ) const

inline

Definition at line 111 of file PndCATrackParam.h.

References fP.

{ return fP; }

float PndCATrackParam::GetQPt ( ) const

inline

Definition at line 89 of file PndCATrackParam.h.

References fP.

Referenced by Print().

{ return fP[4]; }

float PndCATrackParam::GetS	(	float	x,
		float	y,
		float	Bz
	)		const

Definition at line 59 of file PndCATrackParam.cxx.

References CAMath::Abs(), CAMath::ATan2(), GetCosPhi(), GetKappa(), GetSinPhi(), GetX(), and GetY().

Referenced by GetDCAPoint().

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

  float k  = GetKappa( Bz );
  float ex = GetCosPhi();
  float ey = GetSinPhi();
  x -= GetX();
  y -= GetY();
  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;
}

float PndCATrackParam::GetSignCosPhi ( ) const

inline

Definition at line 90 of file PndCATrackParam.h.

References fSignCosPhi.

Referenced by Print().

{ return fSignCosPhi; }

float PndCATrackParam::GetSinPhi ( ) const

inline

Definition at line 87 of file PndCATrackParam.h.

References fP.

Referenced by GetDCAPoint(), GetS(), and Print().

{ return fP[2]; }

float PndCATrackParam::GetX ( ) const

inline

Definition at line 84 of file PndCATrackParam.h.

References fX.

Referenced by GetDCAPoint(), GetDist2(), GetDistXZ2(), GetS(), and Print().

{ return fX; }

float PndCATrackParam::GetY ( ) const

inline

Definition at line 85 of file PndCATrackParam.h.

References fP.

Referenced by GetDCAPoint(), GetDist2(), GetS(), and Print().

{ return fP[0]; }

float PndCATrackParam::GetZ ( ) const

inline

Definition at line 86 of file PndCATrackParam.h.

References fP.

Referenced by GetDCAPoint(), GetDist2(), GetDistXZ2(), and Print().

{ return fP[1]; }

void PndCATrackParam::InitDirection ( float  r0, float  r1, float  r2  )

inline

Definition at line 132 of file PndCATrackParam.h.

References r, SetDzDs(), SetSignCosPhi(), SetSinPhi(), sign(), and sqrt().

  {
    const float r = sqrt( r0*r0+r1*r1 );
    SetSinPhi( r1/r );
    SetSignCosPhi( sign(r0) );
    SetDzDs( r2/r );
  }

int PndCATrackParam::ISec ( ) const

inline

Definition at line 82 of file PndCATrackParam.h.

References fISec.

{ return fISec; }

bool PndCATrackParam::IsValid ( ) const

inline

Definition at line 207 of file PndCATrackParam.h.

References fChi2.

{ return fChi2 != -1; }

int PndCATrackParam::NDF ( ) const

inline

Definition at line 73 of file PndCATrackParam.h.

References fNDF.

Referenced by PndCATrack::Fit().

{ return fNDF; }

const float* PndCATrackParam::Par ( ) const

inline

Definition at line 108 of file PndCATrackParam.h.

References fP.

Referenced by PndCAGBTracker::Merge().

{ return fP; }

void PndCATrackParam::Print

(

)

const

Definition at line 457 of file PndCATrackParam.cxx.

References GetDzDs(), GetErr2DzDs(), GetErr2QPt(), GetErr2SinPhi(), GetErr2Y(), GetErr2Z(), GetQPt(), GetSignCosPhi(), GetSinPhi(), GetX(), GetY(), and GetZ().

{
  //* print parameters

#if !defined(HLTCA_GPUCODE)
  std::cout << "track: x=" << GetX() << " c=" << GetSignCosPhi() << ", P= " << GetY() << " " << GetZ() << " " << GetSinPhi() << " " << GetDzDs() << " " << GetQPt() << std::endl;
  std::cout << "errs2: " << GetErr2Y() << " " << GetErr2Z() << " " << GetErr2SinPhi() << " " << GetErr2DzDs() << " " << GetErr2QPt() << std::endl;
#endif
}

float PndCATrackParam::QMomentum ( ) const

inline

Definition at line 65 of file PndCATrackParam.h.

References fP.

Referenced by PndCANPlet::QMomentum().

{ return fP[4]; }

float PndCATrackParam::QPt ( ) const

inline

Definition at line 64 of file PndCATrackParam.h.

References fP.

Referenced by PndCATracking::CATrackParToFairTrackParP(), and TransportToX().

{ return fP[4]; }

void PndCATrackParam::Reset ( )

inline

Definition at line 202 of file PndCATrackParam.h.

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

Referenced by PndCATrackParam().

{ 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; fISec=0;}

void PndCATrackParam::ResetCovMatrix ( )

inline

Definition at line 191 of file PndCATrackParam.h.

References fC.

    {
      fC[0] = 10.f; 
      fC[1] = 0.f;   fC[2] = 10.f;
      fC[3] = 0.f;   fC[4] = 0.f;  fC[5] = 1.f;
      fC[6] = 0.f;   fC[7] = 0.f;  fC[8] = 0.f;  fC[9] = 1.f;
      fC[10] = 0.f;  fC[11] = 0.f; fC[12] = 0.f; fC[13] = 0.f;  fC[14] = 10.f;
    }

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

Definition at line 290 of file PndCATrackParam.cxx.

References CAMath::Abs(), alpha, CAMath::Cos(), fabs(), fAlpha, fC, GetCosPhi(), SetSignCosPhi(), SetSinPhi(), SetX(), SetY(), CAMath::Sin(), SinPhi(), X(), x, Y(), and y.

Referenced by GetGlobalParam(), 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 = GetCosPhi();
  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;

  SetX( x*cA +  y*sA );
  SetY( -x*sA +  y*cA );
  SetSignCosPhi( CAMath::Abs(cosPhi)/cosPhi );
  SetSinPhi( 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
  //cout<<"alpha="<<alpha<<" "<<x<<" "<<y<<" "<<sP<<" "<<cP<<" "<<j0<<" "<<j2<<endl;
  //cout<<"      "<<fC[0]<<" "<<fC[1]<<" "<<fC[6]<<" "<<fC[10]<<" "<<fC[4]<<" "<<fC[5]<<" "<<fC[8]<<" "<<fC[12]<<endl;
  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;
  //cout<<"      "<<fC[0]<<" "<<fC[1]<<" "<<fC[6]<<" "<<fC[10]<<" "<<fC[4]<<" "<<fC[5]<<" "<<fC[8]<<" "<<fC[12]<<endl;

  fAlpha += alpha;
  
  return 1;
}

bool PndCATrackParam::Rotate	(	float	alpha,
		PndCATrackLinearisation &	t0,
		float	maxSinPhi = .999
	)

Definition at line 338 of file PndCATrackParam.cxx.

References CAMath::Abs(), alpha, CAMath::Cos(), PndCATrackLinearisation::CosPhi(), d, fabs(), fAlpha, fC, PndCATrackLinearisation::SetCosPhi(), PndCATrackLinearisation::SetSinPhi(), SetSinPhi(), SetX(), SetY(), CAMath::Sin(), PndCATrackLinearisation::SinPhi(), SinPhi(), X(), x0, Y(), and y0.

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

  const float cA = CAMath::Cos( alpha );
  const float sA = CAMath::Sin( alpha );
  const float x0 = X(), y0 = Y(), sP = t0.SinPhi(), cP = t0.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
  
  //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

  const float j0 = cP / cosPhi;
  const float j2 = cosPhi / cP;
  const float d[2] = {Y() - y0, SinPhi() - sP};

  SetX( x0*cA +  y0*sA );
  SetY( -x0*sA +  y0*cA + j0*d[0] );
  t0.SetCosPhi( cosPhi );
  t0.SetSinPhi( sinPhi );

  SetSinPhi( sinPhi + j2*d[1] );

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

void PndCATrackParam::RotateXY ( float  alpha, float &  x, float &  y, float &  sin  ) const

inline

Definition at line 224 of file PndCATrackParam.h.

References CAMath::Cos(), GetCosPhi(), CAMath::Sin(), SinPhi(), X(), and Y().

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

  const float cA = CAMath::Cos( alpha );
  const float sA = CAMath::Sin( alpha );

  x = ( X()*cA +  Y()*sA );
  y = ( -X()*sA +  Y()*cA );
  sin = -GetCosPhi() * sA + SinPhi() * cA;

}

void PndCATrackParam::SetAngle ( float  v )

inline

Definition at line 127 of file PndCATrackParam.h.

References fAlpha, and v.

{ fAlpha = v; }

void PndCATrackParam::SetAsInvalid ( )

inline

Definition at line 208 of file PndCATrackParam.h.

References fChi2.

{ fChi2 = -1; }

void PndCATrackParam::SetChi2 ( float  v )

inline

Definition at line 124 of file PndCATrackParam.h.

References fChi2, and v.

{  fChi2 = v; }

void PndCATrackParam::SetCov ( int  i, float  v  )

inline

Definition at line 115 of file PndCATrackParam.h.

References fC, i, and v.

{ fC[i] = v; }

void PndCATrackParam::SetDzDs ( float  v )

inline

Definition at line 121 of file PndCATrackParam.h.

References fP, and v.

Referenced by InitDirection().

{  fP[3] = v; }

void PndCATrackParam::SetErr2QPt ( float  v )

inline

Definition at line 130 of file PndCATrackParam.h.

References fC, and v.

{ fC[14] = v; }

void PndCATrackParam::SetISec ( int  v )

inline

Definition at line 128 of file PndCATrackParam.h.

References fISec, and v.

{ fISec = v; }

void PndCATrackParam::SetNDF ( int  v )

inline

Definition at line 125 of file PndCATrackParam.h.

References fNDF, and v.

{ fNDF = v; }

void PndCATrackParam::SetPar ( int  i, float  v  )

inline

Definition at line 114 of file PndCATrackParam.h.

References fP, i, and v.

{ fP[i] = v; }

void PndCATrackParam::SetQPt ( float  v )

inline

Definition at line 122 of file PndCATrackParam.h.

References fP, and v.

{  fP[4] = v; }

void PndCATrackParam::SetSignCosPhi ( float  v )

inline

Definition at line 123 of file PndCATrackParam.h.

References fSignCosPhi, and v.

Referenced by InitDirection(), and Rotate().

{  fSignCosPhi = v; }

void PndCATrackParam::SetSinPhi ( float  v )

inline

Definition at line 120 of file PndCATrackParam.h.

References fP, and v.

Referenced by InitDirection(), and Rotate().

{  fP[2] = v; }

void PndCATrackParam::SetX ( float  v )

inline

Definition at line 117 of file PndCATrackParam.h.

References fX, and v.

Referenced by Rotate().

{  fX = v;    }

void PndCATrackParam::SetY ( float  v )

inline

Definition at line 118 of file PndCATrackParam.h.

References fP, and v.

Referenced by Rotate().

{  fP[0] = v; }

void PndCATrackParam::SetZ ( float  v )

inline

Definition at line 119 of file PndCATrackParam.h.

References fP, and v.

{  fP[1] = v; }

float PndCATrackParam::SignCosPhi ( ) const

inline

The sign of cos phi is always positive in the slice tracker. Only after coordinate transformation can the sign change to negative.

Definition at line 71 of file PndCATrackParam.h.

References fSignCosPhi.

{ return fSignCosPhi; }

float PndCATrackParam::SinPhi ( ) const

inline

Definition at line 62 of file PndCATrackParam.h.

References fP.

Referenced by PndCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawTrackParam(), GetCosPhi(), GetCosPhiPositive(), Rotate(), RotateXY(), and TransportToX().

{ return fP[2]; }

bool PndCATrackParam::Transport ( const PndCAHit &  hit, float  Bz  )

inline

Definition at line 569 of file PndCATrackParam.h.

References PndCAHit::Angle(), CalculateFitParameters(), fAlpha, Rotate(), TransportToXWithMaterial(), and PndCAHit::X0().

Referenced by PndCAGBTracker::Merge().

{
  // TODO material. See Vector part
  PndCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  PndCATrackLinearisation tR( *this );
  const bool rotated = Rotate(  -fAlpha + hit.Angle(), tR, .999f );
  PndCATrackLinearisation tE( *this );
  const bool transported = TransportToXWithMaterial( hit.X0(), tE, fitPar, Bz, 0.999f );
  return rotated & transported;
}

bool PndCATrackParam::TransportToX	(	float	x,
		float	Bz,
		float	maxSinPhi = .999
	)

Definition at line 182 of file PndCATrackParam.cxx.

References t0.

Referenced by TransportToXWithMaterial().

{
  //* Transport the track parameters to X=x

  PndCATrackLinearisation t0( *this );

  return TransportToX( x, t0, Bz, maxSinPhi );
}

bool PndCATrackParam::TransportToX ( float  x, PndCATrackLinearisation &  t0, float  Bz, float  maxSinPhi = .999, float *  DL = 0  )

inline

Definition at line 369 of file PndCATrackParam.h.

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

{
  //* 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 = t0.CosPhi();
  const float ey = t0.SinPhi();
  const float k  = t0.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 * t0.DzDs();

  if ( DL ) *DL = -dS * CAMath::Sqrt( 1.f + t0.DzDs() * t0.DzDs() );


  const float d[3] = { fP[2] - t0.SinPhi(), fP[3] - t0.DzDs(), fP[4] - t0.QPt() };

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

  t0.SetCosPhi( ex1 );
  t0.SetSinPhi( ey1 );

  fX    = X() + dx;
  fP[0] = Y() + dy     + h2 * d[0]           +   h4 * d[2];
  fP[1] = Z() + dz               + dS * d[1];
  fP[2] = t0.SinPhi() +     d[0]           + dxBz * d[2];
  if(CAMath::Abs(fP[2]) > maxSinPhi) fP[2] = t0.SinPhi();

#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

//std::cout << fC[0] << "  "<<fC1[0]<<"       "<<fC[2] << "  "<<fC1[2]<<"       "<<fC[5] << "  "<<fC1[5]<<"       "<<fC[9] << "  "<<fC1[9]<<"       "<<fC[14] << "  "<<fC1[14]<<std::endl;
  return 1;
}

bool PndCATrackParam::TransportToX	(	float	x,
		float	sinPhi0,
		float	cosPhi0,
		float	Bz,
		float	maxSinPhi = .999
	)

Definition at line 106 of file PndCATrackParam.cxx.

References CAMath::Abs(), c10, c11, c20, c21, c22, dx, DzDs(), fC, fP, fX, h2, h4, QPt(), SinPhi(), and X().

{
  //* Transport the track parameters to X=x, using linearization at phi0 with 0 curvature,
  //* 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 = cosPhi0;
  const float ey = sinPhi0;
  const float dx = x - X();

  if ( CAMath::Abs( ex ) < 1.e-4 ) return 0;
  const float exi = 1. / ex;
  
  const float dxBz = dx * Bz;
  const float dS = dx * exi;
  const float h2 = dS * exi * exi;
  const float h4 = .5 * h2 * dxBz;

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

  const float sinPhi = SinPhi() + dxBz * QPt();
  if ( maxSinPhi > 0 && CAMath::Abs( sinPhi ) > maxSinPhi ) return 0;

  fX    = X() + dx;
  fP[0] += dS * ey + h2 * ( SinPhi() - ey )  +   h4 * QPt();
  fP[1] += dS * DzDs();
  fP[2] = sinPhi;

  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 * ( h2 * c20 + h4 * c40 + h2 * h4 * c42 )  );

  fC[1] = c10 + h2 * c21 + h4 * c41 + dS * ( c30 + h2 * c32 + h4 * c43 );
  fC[2] = c11 + 2 * 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 * 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;

  return 1;
}

bool PndCATrackParam::TransportToXWithMaterial	(	float	x,
		float	Bz,
		float	maxSinPhi = .999
	)

Definition at line 191 of file PndCATrackParam.cxx.

References CalculateFitParameters(), and par.

Referenced by Transport(), and TransportToXWithMaterial().

{
  //* Transport the track parameters to X=x taking into account material budget

  PndCATrackFitParam par;
  CalculateFitParameters( par );
  return TransportToXWithMaterial( x, par, Bz, maxSinPhi );
}

bool PndCATrackParam::TransportToXWithMaterial ( float  x, PndCATrackLinearisation &  t0, PndCATrackFitParam &  par, float  Bz, float  maxSinPhi = .999  )

inline

Definition at line 543 of file PndCATrackParam.h.

References TransportToX().

{
  //* Transport the track parameters to X=x  taking into account material budget

  const float kRho = 1.54e-3;//1.025e-3 ;//0.9e-3;
//  const float kRadLen = 29.532;//28.94;
  //const float kRhoOverRadLen = kRho / kRadLen;
  const float kRhoOverRadLen = 7.68e-5;
  float dl;

  if ( !TransportToX( x, t0, Bz,  maxSinPhi, &dl ) ) return 0;

  UNUSED_PARAM3(kRho, kRhoOverRadLen, par);//TODO
  // CorrectForMeanMaterial( dl*kRhoOverRadLen, dl*kRho, par );
  return 1;
}

bool PndCATrackParam::TransportToXWithMaterial ( float  x, PndCATrackFitParam &  par, float  Bz, float  maxSinPhi = .999  )

inline

Definition at line 560 of file PndCATrackParam.h.

References t0, and TransportToXWithMaterial().

{
  //* Transport the track parameters to X=x  taking into account material budget

  PndCATrackLinearisation t0( *this );
  return TransportToXWithMaterial( x, t0, par, Bz, maxSinPhi );
}

float PndCATrackParam::X ( ) const

inline

Definition at line 59 of file PndCATrackParam.h.

References fX.

Referenced by PndCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawTrackParam(), Rotate(), RotateXY(), and TransportToX().

{ return fX;    }

float PndCATrackParam::Y ( ) const

inline

Definition at line 60 of file PndCATrackParam.h.

References fP.

Referenced by PndCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawTrackParam(), Rotate(), RotateXY(), and TransportToX().

{ return fP[0]; }

float PndCATrackParam::Z ( ) const

inline

Definition at line 61 of file PndCATrackParam.h.

References fP.

Referenced by PndCATracking::CATrackParToFairTrackParP(), PndFTSCADisplay::DrawTrackParam(), and TransportToX().

{ return fP[1]; }

Member Data Documentation

float PndCATrackParam::fAlpha

private

Definition at line 220 of file PndCATrackParam.h.

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

float PndCATrackParam::fC[15]

private

Definition at line 216 of file PndCATrackParam.h.

Referenced by CorrectForMeanMaterial(), Cov(), Err2DzDs(), Err2QMomentum(), Err2QPt(), Err2SinPhi(), Err2X1(), Err2X2(), Err2Y(), Err2Z(), Filter(), GetCov(), GetErr2DzDs(), GetErr2QPt(), GetErr2SinPhi(), GetErr2Y(), GetErr2Z(), PndCATrackParam(), Reset(), ResetCovMatrix(), Rotate(), SetCov(), SetErr2QPt(), and TransportToX().

float PndCATrackParam::fChi2

private

Definition at line 217 of file PndCATrackParam.h.

Referenced by Chi2(), Filter(), GetChi2(), IsValid(), Reset(), SetAsInvalid(), and SetChi2().

int PndCATrackParam::fISec

private

Definition at line 221 of file PndCATrackParam.h.

Referenced by ISec(), Reset(), and SetISec().

int PndCATrackParam::fNDF

private

Definition at line 218 of file PndCATrackParam.h.

Referenced by Filter(), GetNDF(), NDF(), Reset(), and SetNDF().

float PndCATrackParam::fP[5]

private

Definition at line 215 of file PndCATrackParam.h.

Referenced by CalculateFitParameters(), CorrectForMeanMaterial(), DzDs(), Filter(), GetDzDs(), GetKappa(), GetPar(), GetQPt(), GetSinPhi(), GetY(), GetZ(), Par(), PndCATrackParam(), QMomentum(), QPt(), Reset(), SetDzDs(), SetPar(), SetQPt(), SetSinPhi(), SetY(), SetZ(), SinPhi(), TransportToX(), Y(), and Z().

float PndCATrackParam::fSignCosPhi

private

Definition at line 214 of file PndCATrackParam.h.

Referenced by GetCosPhi(), GetSignCosPhi(), Reset(), SetSignCosPhi(), and SignCosPhi().

float PndCATrackParam::fX

private

Definition at line 213 of file PndCATrackParam.h.

Referenced by GetX(), Reset(), SetX(), TransportToX(), and X().

---

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

• PndCATrackParam.h
• PndCATrackParam.cxx