PndCATrackParamVector Class Reference

#include <PndCATrackParamVector.h>

Classes
struct	PndCATrackFitParam

Public Member Functions
	PndCATrackParamVector ()
void	Reset ()
void	InitCovMatrix (float_v d2QMom=0.f)
void	InitByTarget (const PndCATarget &target)
void	InitByHit (const PndCAHitV &hit, const PndCAParam &param, const float_v &dQMom)
void	InitDirection (float_v r0, float_v r1, float_v r2)
float_v	X () const
float_v	Y () const
float_v	Z () const
float_v	SinPhi () const
float_v	DzDs () const
float_v	QPt () const
float_v	Angle () const
int_v	ISec () const
float_v	X0 () const
float_v	X1 () const
float_v	X2 () const
float_v	Tx1 () const
float_v	Tx2 () const
float_v	QMomentum () const
float_v	SignCosPhi () const
float_v	Chi2 () const
int_v	NDF () const
float_v	Err2Y () const
float_v	Err2Z () const
float_v	Err2SinPhi () const
float_v	Err2DzDs () const
float_v	Err2QPt () const
float_v	GetX () const
float_v	GetY () const
float_v	GetZ () const
float_v	GetSinPhi () const
float_v	GetDzDs () const
float_v	GetQPt () const
float_v	GetSignCosPhi () const
float_v	GetChi2 () const
int_v	GetNDF () const
float_v	GetKappa (const float_v &Bz) const
float_v	GetCosPhiPositive () const
float_v	GetCosPhi () const
float_v	Err2X1 () const
float_v	Err2X2 () const
float_v	Err2QMomentum () const
const float_v &	Par (int i) const
const float_v &	Cov (int i) const
void	SetTrackParam (const PndCATrackParamVector &param, const float_m &m=float_m(true))
void	SetTrackParam (const PndCATrackParamVector &p, int k, int m)
void	SetPar (int i, const float_v &v)
void	SetPar (int i, const float_v &v, const float_m &m)
void	SetCov (int i, const float_v &v)
void	SetCov (int i, const float_v &v, const float_m &m)
void	SetX (const float_v &v)
void	SetY (const float_v &v)
void	SetZ (const float_v &v)
void	SetX (const float_v &v, const float_m &m)
void	SetY (const float_v &v, const float_m &m)
void	SetZ (const float_v &v, const float_m &m)
void	SetSinPhi (const float_v &v)
void	SetSinPhi (const float_v &v, const float_m &m)
void	SetDzDs (const float_v &v)
void	SetDzDs (const float_v &v, const float_m &m)
void	SetQPt (const float_v &v)
void	SetQMomentum (const float_v &v)
void	SetQPt (const float_v &v, const float_m &m)
void	SetSignCosPhi (const float_v &v)
void	SetSignCosPhi (const float_v &v, const float_m &m)
void	SetChi2 (const float_v &v)
void	SetChi2 (const float_v &v, const float_m &m)
void	SetNDF (int v)
void	SetNDF (const int_v &v)
void	SetNDF (const int_v &v, const int_m &m)
void	SetAngle (const float_v &v)
void	SetAngle (const float_v &v, const float_m &m)
void	SetISec (const int_v &v)
void	SetISec (const int_v &v, const int_m &m)
void	SetErr2Y (float_v v)
void	SetErr2Z (float_v v)
void	SetErr2QPt (float_v v)
float_v	GetDist2 (const PndCATrackParamVector &t) const
float_v	GetDistXZ2 (const PndCATrackParamVector &t) const
float_v	GetS (const float_v &x, const float_v &y, const float_v &Bz) const
void	GetDCAPoint (const float_v &x, const float_v &y, const float_v &z, float_v *px, float_v *py, float_v *pz, const float_v &Bz) const
float_m	Transport0 (const int_v &ista, const PndCAParam &param, const float_m &mask=float_m(true))
float_m	Transport0 (const PndCAHitV &hit, const PndCAParam &param, const float_m &mask=float_m(true))
float_m	Transport0 (const PndCAHit &hit, const PndCAParam &p, const float_m &mask=float_m(true))
float_m	TransportToXWithMaterial (const float_v &x, const float_v &XOverX0, const float_v &XThimesRho, const float_v &Bz, const float maxSinPhi=.999f)
float_m	TransportToX (const float_v &x, const float_v &Bz, const float maxSinPhi=.999f, const float_m &mask=float_m(true))
float_m	TransportToX (const float_v &x, PndCATrackLinearisationVector &t0, const float_v &Bz, const float maxSinPhi=.999f, float_v *DL=0, const float_m &mask=float_m(true))
float_m	Transport0ToX (const float_v &x, const float_v &Bz, const float_m &mask)
float_m	TransportToX (const float_v &x, const float_v &sinPhi0, const float_v &Bz, const float_v maxSinPhi=.999f, const float_m &mask=float_m(true))
float_m	TransportToXWithMaterial (const float_v &x, PndCATrackLinearisationVector &t0, PndCATrackFitParam &par, const float_v &XOverX0, const float_v &XThimesRho, const float_v &Bz, const float maxSinPhi=.999f, const float_m &mask=float_m(true))
float_m	TransportToXWithMaterial (const float_v &x, PndCATrackFitParam &par, const float_v &XOverX0, const float_v &XThimesRho, const float_v &Bz, const float maxSinPhi=.999f)
float_m	Rotate (const float_v &alpha, PndCATrackLinearisationVector &t0, const float maxSinPhi=.999f, const float_m &mask=float_m(true))
float_m	Rotate (const float_v &alpha, const float maxSinPhi=.999f, const float_m &mask=float_m(true))
void	RotateXY (float_v alpha, float_v &x, float_v &y, float_v &sin, const float_m &mask=float_m(true)) const
float_m	FilterWithMaterial (const float_v &y, const float_v &z, float_v err2Y, float_v errYZ, float_v err2Z, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const int_v &hitNDF=int_v(2), const float_v &chi2Cut=10e10f)
float_m	FilterWithMaterial (const float_v &y, const float_v &z, const PndCAStripInfo &info, float_v err2, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
void	CalculateFitParameters (PndCATrackFitParam &par, const float_v &mass=0.13957f)
float_m	CorrectForMeanMaterial (const float_v &xOverX0, const float_v &xTimesRho, const PndCATrackFitParam &par, const float_m &_mask)
float_m	FilterVtx (const float_v &xV, const float_v &yV, const PndCAX1X2MeasurementInfo &info, float_v &extrDx, float_v &extrDy, float_v J[], const float_m &active=float_m(true))
float_m	TransportJ0ToX0 (const float_v &x0, const float_v &cBz, float_v &extrDy, float_v &extrDz, float_v J[], const float_m &active=float_m(true))
float_m	Transport (const int_v &ista, const PndCAParam &param, const float_m &mask=float_m(true))
float_m	Transport (const PndCAHitV &hit, const PndCAParam &p, const float_m &mask=float_m(true))
float_m	Filter (const PndCAHitV &hit, const PndCAParam &param, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
float_m	Transport (const PndCAHit &hit, const PndCAParam &p, const float_m &mask=float_m(true))
float_m	Filter (const PndCAHit &hit, const PndCAParam &param, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
float_m	Accept (const PndCAHit &hit, const PndCAParam &param, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f) const
float_m	AcceptWithMaterial (const float_v &y, const float_v &z, float_v err2Y, float_v errYZ, float_v err2Z, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const int_v &hitNDF=int_v(2), const float_v &chi2Cut=10e10f) const
float_m	AcceptWithMaterial (const float_v &y, const float_v &z, const PndCAStripInfo &info, float_v err2, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f) const
float_m	AddTarget (const PndCATarget &target, const float_m &mask=float_m(true))

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

Public Attributes
float_v	fX
float_v	fSignCosPhi
float_v	fP [5]
float_v	fC [15]
float_v	fChi2
int_v	fNDF
float_v	fAlpha
int_v	fISec

Private Member Functions
const float_v *	Par () const
const float_v *	Cov () const
float_v *	Par ()
float_v *	Cov ()

Friends
class	PndCATrackParam

Detailed Description

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

Definition at line 38 of file PndCATrackParamVector.h.

Constructor & Destructor Documentation

PndCATrackParamVector::PndCATrackParamVector ( )

inline

Definition at line 43 of file PndCATrackParamVector.h.

References fC, fP, and i.

      : fX( Vc::Zero ),
      fSignCosPhi( Vc::Zero ),
      fChi2( Vc::Zero ),
        fNDF( Vc::Zero ),
        fISec(Vc::Zero)
    {
      for ( int i = 0; i <  5; ++i ) fP[i].setZero();
      for ( int i = 0; i < 15; ++i ) fC[i].setZero();
    }

Member Function Documentation

float_m PndCATrackParamVector::Accept	(	const PndCAHit &	hit,
		const PndCAParam &	param,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)		const

Definition at line 1026 of file PndCATrackParamVector.cxx.

References AcceptWithMaterial(), PndCAHit::Err2R(), PndCAHit::Err2X1(), PndCAHit::Err2X2(), PndCAHit::ErrX12(), f, PndCAHit::IStation(), NDF(), PndCAParam::Station(), PndCAHit::X1(), PndCAHit::X1Corrected(), and PndCAHit::X2().

Referenced by PndCAGBTracker::PickUpHits().

{
  if ( param.Station( hit.IStation() ).NDF == 1 ) {
    return AcceptWithMaterial( hit.X1Corrected(this->SinPhi()), hit.X2(), param.Station( hit.IStation() ).f, hit.Err2R(), 0.999f, mask, chi2Cut );
  }
  return AcceptWithMaterial( hit.X1(), hit.X2(), hit.Err2X1(), hit.ErrX12(), hit.Err2X2(), 0.999f, mask, int_v(param.Station( hit.IStation() ).NDF), chi2Cut );
}

float_m PndCATrackParamVector::AcceptWithMaterial	(	const float_v &	y,
		const float_v &	z,
		float_v	err2Y,
		float_v	errYZ,
		float_v	err2Z,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const int_v &	hitNDF = int_v(2),
		const float_v &	chi2Cut = 10e10f
	)		const

Definition at line 1034 of file PndCATrackParamVector.cxx.

References CAMath::Abs(), c10, c11, c20, c21, f, fC, fChi2, fP, mrfdata_8b_error::success, and z0.

Referenced by Accept().

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

  const float_v c00 = fC[0];
  const float_v c10 = fC[1];
  const float_v c11 = fC[2];
  const float_v c20 = fC[3];
  const float_v c21 = fC[4];
  const float_v
    z0 = y - fP[0],
    z1 = z - fP[1];
    
  err2Y += c00;
  err2Z += c11;
  errYZ += c10;
  float_v d1 = ( err2Y*err2Z - errYZ*errYZ );

  float_m success = mask;//  if ( ISUNLIKELY( err2Y < 1.e-8f ) || ISUNLIKELY( err2Z < 1.e-8f ) ) return 0;
  success &= (err2Y > 1.e-8f ) && ( err2Z > 1.e-8f );

  const float_v mS0 = err2Z;
  const float_v mS1 = -errYZ;
  const float_v mS2 = err2Y;

  // K = CHtS

  const float_v k20 = (c20 * mS0 + c21*mS1),   k21 = (c20 * mS1 + c21*mS2);
  const float_v sinPhi = d1*fP[2] + ( k20 * z0 + k21 * z1 );
  success &= CAMath::Abs( sinPhi ) < maxSinPhi*d1;
  float_v chi2 = ( mS0 * z0 * z0 + 2 * z0 * z1 * mS1 + mS2 * z1 * z1 );
  success &=  chi2 < d1*(chi2Cut-fChi2);
  return success;
}

float_m PndCATrackParamVector::AcceptWithMaterial	(	const float_v &	y,
		const float_v &	z,
		const PndCAStripInfo &	info,
		float_v	err2,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)		const

Definition at line 1071 of file PndCATrackParamVector.cxx.

References c10, c11, f, mrfdata_8b_error::success, y, and z.

{
  assert( maxSinPhi > 0.f );
  float_m success = mask;
  success &= ( err2 > 1.e-8f );

  const float_v& c00 = fC[0];
  const float_v& c10 = fC[1];
  const float_v& c11 = fC[2];

  const float_v& sb = info.sin;
  const float_v& cb = info.cos;

  const float_v u = cb*y + sb*z;
  const float_v zeta = cb*fP[0] + sb*fP[1] - u;

  // F = CH'
  const float_v  F0 = cb*c00 + sb*c10;
  const float_v  F1 = cb*c10 + sb*c11;
        
  const float_v  HCH = ( F0*cb + F1*sb );

  const float_v  w = ( err2 + HCH );
  success &= zeta*zeta < w*(chi2Cut - fChi2);
  return success;
}

float_m PndCATrackParamVector::AddTarget	(	const PndCATarget &	target,
		const float_m &	mask = float_m( true )
	)

Definition at line 870 of file PndCATrackParamVector.cxx.

References PndCATarget::B(), PndCATarget::Err2X1(), PndCATarget::Err2X2(), FilterVtx(), fNDF, PndCATarget::NDF(), TransportJ0ToX0(), PndCATarget::X0(), PndCATarget::X1(), and PndCATarget::X2().

{
  float_m active = mask;

  float_v eY, eZ;
  //float_v dz = target.X0() - X0();
  float_v J[6];
    // H = 1 0 J[0] J[1] J[2]
    //     0 1 J[3] J[4] J[5]
  active &= TransportJ0ToX0( target.X0(), static_cast<float_v>(target.B()), eY, eZ, J, active );
  active &= FilterVtx( target.X1(), target.X2(), PndCAX1X2MeasurementInfo( target.Err2X1(), 0.f, target.Err2X2() ), eY, eZ, J, active );
  fNDF (static_cast<int_m>(active)) += target.NDF();

  return active;
}

float_v PndCATrackParamVector::Angle ( ) const

inline

Definition at line 96 of file PndCATrackParamVector.h.

References fAlpha.

Referenced by PndCAGBTracker::PickUpHits().

{ return fAlpha; }

float_v PndCATrackParamVector::ApproximateBetheBloch ( const float_v &  beta2 )

static

Definition at line 471 of file PndCATrackParamVector.cxx.

References f, CAMath::Log(), and One.

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

  const float_v &beta2_1subBeta2 = beta2 / ( float_v( Vc::One ) - beta2 ); // beta2 * CAMath::Reciprocal( float_v( Vc::One ) - beta2 );
  const float_v &_0p000153_beta2 = 0.153e-3f / beta2;
  const float log_3p5mul5940 = 9.942227380852058f; // log( 3.5 * 5940 )
  const float log_5940 = 8.68946441235669f; // log( 5940 )
  const float_v log_beta2_1subBeta2 = CAMath::Log( beta2_1subBeta2 );
    
  float_v ret = _0p000153_beta2 * ( log_5940 + log_beta2_1subBeta2 - beta2 );
  ret( beta2_1subBeta2 > 3.5f*3.5f ) =
    _0p000153_beta2 * ( log_3p5mul5940 + 0.5f * log_beta2_1subBeta2 - beta2 );
  ret.setZero( beta2 >= float_v( Vc::One ) );
  return ret;
}

float_v PndCATrackParamVector::BetheBlochGas ( const float_v &  bg )

static

Definition at line 450 of file PndCATrackParamVector.cxx.

References BetheBlochGeant(), f, 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_v rho = 0.9e-3f;
  const float_v x0  = 2.f;
  const float_v x1  = 4.f;
  const float_v mI  = 140.e-9f;
  const float_v mZA = 0.49555f;

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

float_v PndCATrackParamVector::BetheBlochGeant ( const float_v &  bg, const float_v &  kp0 = 2.33f, const float_v &  kp1 = 0.20f, const float_v &  kp2 = 3.00f, const float_v &  kp3 = 173e-9f, const float_v &  kp4 = 0.49848f  )

static

Definition at line 397 of file PndCATrackParamVector.cxx.

References f, CAMath::Log(), mK, One, r, CAMath::Sqrt(), x, x0, and Zero.

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-3f; // [GeV*cm^2/g]
  const float _2me  = 1.022e-3f;    // [GeV/c^2]
  const float_v &rho = kp0;
  const float_v &x0  = kp1 * 2.303f;
  const float_v &x1  = kp2 * 2.303f;
  const float_v &mI  = kp3;
  const float_v &mZA = kp4;
  const float_v &maxT = _2me * bg2;    // neglecting the electron mass

  //*** Density effect
  float_v d2( Vc::Zero );
  const float_v x = 0.5f * CAMath::Log( bg2 );
  const float_v lhwI = CAMath::Log( 28.816f * 1e-9f * CAMath::Sqrt( rho * mZA ) / mI );
  d2( x > x1 ) = lhwI + x - 0.5f;
  const float_v &r = ( x1 - x ) / ( x1 - x0 );
  d2( x > x0 && x <= x1 ) = lhwI + x - 0.5f + ( 0.5f - lhwI - x0 ) * r * r * r;

  return mK*mZA*( float_v( Vc::One ) + bg2 ) / bg2*( 0.5f*CAMath::Log( maxT * maxT / ( mI*mI ) ) - bg2 / ( float_v( Vc::One ) + bg2 ) - d2 );
}

float_v PndCATrackParamVector::BetheBlochSolid ( const float_v &  bg )

static

Definition at line 438 of file PndCATrackParamVector.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 PndCATrackParamVector::CalculateFitParameters	(	PndCATrackFitParam &	par,
		const float_v &	mass = 0.13957f
	)

Definition at line 493 of file PndCATrackParamVector.cxx.

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

Referenced by TransportToXWithMaterial().

{
  //*!

  const float_v p2 = ( float_v( Vc::One ) + fP[3] * fP[3] );
  const float_v k2 = fP[4] * fP[4];
  const float_v mass2 = mass * mass;
  const float_v beta2 = p2 / ( p2 + mass2 * k2 );

  float_v pp2 = 10000.f; pp2( k2 > 1.e-8f ) = p2 / k2; // impuls 2

  //par.fBethe = BetheBlochGas( pp2/mass2);
  // par.fBethe = ApproximateBetheBloch( pp2 / mass2 ); Why pp2 / mass2?, should be beta2.
  par.fBethe = ApproximateBetheBloch( beta2 );
  par.fE = CAMath::Sqrt( pp2 + mass2 );
  par.fTheta2 = 14.1f * 14.1f / ( beta2 * pp2 * 1.e6f );
  par.fEP2 = par.fE / pp2;

  // Approximate energy loss fluctuation (M.Ivanov)

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

  par.fK22 = ( float_v( Vc::One ) + fP[3] * fP[3] );
  par.fK33 = par.fK22 * par.fK22;
  par.fK43 = fP[3] * fP[4] * par.fK22;
  par.fK44 = fP[3] * fP[3] * fP[4] * fP[4];
}

float_v PndCATrackParamVector::Chi2 ( ) const

inline

Definition at line 111 of file PndCATrackParamVector.h.

References fChi2.

{ return fChi2; }

float_m PndCATrackParamVector::CorrectForMeanMaterial	(	const float_v &	xOverX0,
		const float_v &	xTimesRho,
		const PndCATrackFitParam &	par,
		const float_m &	_mask
	)

Definition at line 524 of file PndCATrackParamVector.cxx.

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

Referenced by TransportToXWithMaterial().

{
  //------------------------------------------------------------------
  // 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_v &fC22 = fC[5];
  float_v &fC33 = fC[9];
  float_v &fC40 = fC[10];
  float_v &fC41 = fC[11];
  float_v &fC42 = fC[12];
  float_v &fC43 = fC[13];
  float_v &fC44 = fC[14];

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

  const float_v &dE = par.fBethe * xTimesRho;
  assert( dE >= 0 );
  float_m mask = _mask && dE <= 0.3f * par.fE; //30% energy loss is too much!
  const float_v &corr = ( float_v( Vc::One ) - par.fEP2 * dE );
  mask &= corr >= 0.3f && corr <= 1.3f;

  fP[4]( mask ) *= corr;
  fC40 ( mask ) *= corr;
  fC41 ( mask ) *= corr;
  fC42 ( mask ) *= corr;
  fC43 ( mask ) *= corr;
  fC44 ( mask ) *= corr * corr;
  fC44 ( mask ) += par.fSigmadE2 * CAMath::Abs( dE );

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

  assert( xOverX0 >= 0 );
  const float_v &theta2 = par.fTheta2 * xOverX0;
  fC22( mask ) += theta2 * par.fK22 * ( float_v( Vc::One ) - fP[2] * fP[2] );
  fC33( mask ) += theta2 * par.fK33;
  fC43( mask ) += theta2 * par.fK43;
  fC44( mask ) += theta2 * par.fK44;

  return mask;
}

const float_v& PndCATrackParamVector::Cov ( int  i ) const

inline

Definition at line 141 of file PndCATrackParamVector.h.

References fC, and i.

Referenced by PndCATrackParam::PndCATrackParam(), PndFTSCATrackParam::PndFTSCATrackParam(), and SetTrackParam().

{ return fC[i]; }

const float_v* PndCATrackParamVector::Cov ( ) const

inlineprivate

Definition at line 146 of file PndCATrackParamVector.h.

References fC.

{ return fC; }

float_v* PndCATrackParamVector::Cov ( )

inlineprivate

Definition at line 148 of file PndCATrackParamVector.h.

References fC.

{ return fC; }

float_v PndCATrackParamVector::DzDs ( ) const

inline

Definition at line 93 of file PndCATrackParamVector.h.

References fP.

Referenced by PndCAGBTracker::PickUpHits(), Transport0ToX(), TransportJ0ToX0(), TransportToX(), and Tx2().

{ return fP[3]; }

float_v PndCATrackParamVector::Err2DzDs ( ) const

inline

Definition at line 117 of file PndCATrackParamVector.h.

References fC.

{ return fC[9]; }

float_v PndCATrackParamVector::Err2QMomentum ( ) const

inline

Definition at line 138 of file PndCATrackParamVector.h.

References fC.

{ return fC[14]; }

float_v PndCATrackParamVector::Err2QPt ( ) const

inline

Definition at line 118 of file PndCATrackParamVector.h.

References fC.

{ return fC[14]; }

float_v PndCATrackParamVector::Err2SinPhi ( ) const

inline

Definition at line 116 of file PndCATrackParamVector.h.

References fC.

{ return fC[5]; }

float_v PndCATrackParamVector::Err2X1 ( ) const

inline

Definition at line 134 of file PndCATrackParamVector.h.

References fC.

Referenced by PndCAGBTracker::PickUpHits().

{ return fC[0]; }

float_v PndCATrackParamVector::Err2X2 ( ) const

inline

Definition at line 135 of file PndCATrackParamVector.h.

References fC.

Referenced by PndCAGBTracker::PickUpHits().

{ return fC[2]; }

float_v PndCATrackParamVector::Err2Y ( ) const

inline

Definition at line 114 of file PndCATrackParamVector.h.

References fC.

{ return fC[0]; }

float_v PndCATrackParamVector::Err2Z ( ) const

inline

Definition at line 115 of file PndCATrackParamVector.h.

References fC.

{ return fC[2]; }

float_m PndCATrackParamVector::Filter	(	const PndCAHitV &	hit,
		const PndCAParam &	param,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 832 of file PndCATrackParamVector.cxx.

References PndCAHitV::Err2R(), PndCAHitV::Err2X1(), PndCAHitV::Err2X2(), PndCAHitV::ErrX12(), f, PndCAStation::f, FilterWithMaterial(), PndCAHitV::IStations(), PndCAHitV::IsValid(), m, PndCAStation::NDF, NDF(), PndCAParam::Station(), PndCAHitV::X1(), PndCAHitV::X1Corrected(), and PndCAHitV::X2().

Referenced by PndCAGBTracker::Create1Plets(), PndCAFunctionality::FitIteration(), and PndCAGBTracker::PickUpHits().

{
  const int iSta = hit.IStations()[hit.IsValid().firstOne()];
  if ( param.Station( iSta ).NDF == 1 ) {
    float_m m = FilterWithMaterial( hit.X1Corrected(*this), hit.X2(), param.Station( iSta ).f, hit.Err2R(), 0.999f, mask, chi2Cut );
    return m;
  }
  //SG!!
  float_m m = FilterWithMaterial( hit.X1(), hit.X2(), float_v(3.f)*hit.Err2X1(), hit.ErrX12(), float_v(3.f)*hit.Err2X2(), 0.999f, mask, int_v(param.Station( hit.IStations()[0] ).NDF), chi2Cut );
  return m;
}

float_m PndCATrackParamVector::Filter	(	const PndCAHit &	hit,
		const PndCAParam &	param,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 862 of file PndCATrackParamVector.cxx.

References PndCAHit::Err2R(), PndCAHit::Err2X1(), PndCAHit::Err2X2(), PndCAHit::ErrX12(), f, FilterWithMaterial(), PndCAHit::IStation(), NDF(), PndCAParam::Station(), PndCAHit::X1(), PndCAHit::X1Corrected(), and PndCAHit::X2().

{
  if ( param.Station( hit.IStation() ).NDF == 1 ) {
    return FilterWithMaterial( hit.X1Corrected(this->SinPhi()), hit.X2(), param.Station( hit.IStation() ).f, hit.Err2R(), 0.999f, mask, chi2Cut );
  }
  return FilterWithMaterial( hit.X1(), hit.X2(), hit.Err2X1(), hit.ErrX12(), hit.Err2X2(), 0.999f, mask, int_v(param.Station( hit.IStation() ).NDF), chi2Cut );
}

float_m PndCATrackParamVector::FilterVtx ( const float_v &  xV, const float_v &  yV, const PndCAX1X2MeasurementInfo &  info, float_v &  extrDx, float_v &  extrDy, float_v  J[], const float_m &  active = float_m( true )  )

inline

Definition at line 456 of file PndCATrackParamVector.h.

References PndCAX1X2MeasurementInfo::C00(), PndCAX1X2MeasurementInfo::C10(), c10, PndCAX1X2MeasurementInfo::C11(), c11, c20, c21, c22, fC, fChi2, fP, Y(), and Z().

Referenced by AddTarget().

{
  float_v zeta0, zeta1, S00, S10, S11, si;
  float_v F00, F10, F20, F30, F40, F01, F11, F21, F31, F41 ;
  float_v K00, K10, K20, K30, K40, K01, K11, K21, K31, K41;
  float_v& c00 = fC[0];
  float_v& c10 = fC[1];
  float_v& c11 = fC[2];
  float_v& c20 = fC[3];
  float_v& c21 = fC[4];
  float_v& c22 = fC[5];
  float_v& c30 = fC[6];
  float_v& c31 = fC[7];
  float_v& c32 = fC[8];
  float_v& c33 = fC[9];
  float_v& c40 = fC[10];
  float_v& c41 = fC[11];
  float_v& c42 = fC[12];
  float_v& c43 = fC[13];
  float_v& c44 = fC[14];
  
  zeta0 = Y() + extrDy - yV;
  zeta1 = Z() + extrDz - zV;

    // H = 1 0 J[0] J[1] J[2]
    //     0 1 J[3] J[4] J[5]
  
  // F = CH'
  F00 = c00;       F01 = c10;    
  F10 = c10;       F11 = c11;
  F20 = J[0]*c22;  F21 = J[3]*c22;
  F30 = J[1]*c33;  F31 = J[4]*c33;
  F40 = J[2]*c44;  F41 = J[5]*c44;

  S00 = info.C00() + F00 + J[0]*F20 + J[1]*F30 + J[2]*F40;
  S10 = info.C10() + F10 + J[3]*F20 + J[4]*F30 + J[5]*F40;
  S11 = info.C11() + F11 + J[3]*F21 + J[4]*F31 + J[5]*F41;

  si = 1.f/(S00*S11 - S10*S10);
  float_v S00tmp = S00;
  S00 =  si*S11;
  S10 = -si*S10;
  S11 =  si*S00tmp;

  fChi2(active) +=  zeta0*zeta0*S00 + 2.f*zeta0*zeta1*S10 + zeta1*zeta1*S11;

  K00 = F00*S00 + F01*S10;  K01 = F00*S10 + F01*S11;
  K10 = F10*S00 + F11*S10;  K11 = F10*S10 + F11*S11;
  K20 = F20*S00 + F21*S10;  K21 = F20*S10 + F21*S11;
  K30 = F30*S00 + F31*S10;  K31 = F30*S10 + F31*S11;
  K40 = F40*S00 + F41*S10;  K41 = F40*S10 + F41*S11;

  fP[0](active)  -= K00*zeta0 + K01*zeta1;
  fP[1](active)  -= K10*zeta0 + K11*zeta1;
  fP[2](active)  -= K20*zeta0 + K21*zeta1;
  fP[3](active)  -= K30*zeta0 + K31*zeta1;
  fP[4](active)  -= K40*zeta0 + K41*zeta1;

  c00(active) -=  ( K00*F00 + K01*F01 );
  c10(active) -=  ( K10*F00 + K11*F01 );
  c11(active) -=  ( K10*F10 + K11*F11 );
  c20(active)  = -( K20*F00 + K21*F01 );
  c21(active)  = -( K20*F10 + K21*F11 );
  c22(active) -=  ( K20*F20 + K21*F21 );
  c30(active)  = -( K30*F00 + K31*F01 );
  c31(active)  = -( K30*F10 + K31*F11 );
  c32(active)  = -( K30*F20 + K31*F21 );
  c33(active) -=  ( K30*F30 + K31*F31 );
  c40(active)  = -( K40*F00 + K41*F01 );
  c41(active)  = -( K40*F10 + K41*F11 );
  c42(active)  = -( K40*F20 + K41*F21 );
  c43(active)  = -( K40*F30 + K41*F31 );
  c44(active) -=  ( K40*F40 + K41*F41 );
  
  return active;
}

float_m PndCATrackParamVector::FilterWithMaterial	(	const float_v &	y,
		const float_v &	z,
		float_v	err2Y,
		float_v	errYZ,
		float_v	err2Z,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const int_v &	hitNDF = int_v(2),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 620 of file PndCATrackParamVector.cxx.

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

Referenced by Filter().

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

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


  const float_v
    z0 = y - fP[0],
    z1 = z - fP[1];
  
    // foreach_bit( int i, mask ) {
    //   cout << i << " e = " << err2Y[i] << " " << errYZ[i] << " " << err2Z[i] << " c = " << c00[i] << " " << c10[i] << " " << c11[i] <<  " z = " << z0[i] <<  " " << z1[i]  << " chi2 = " << fChi2[i]<< endl;
    // }
  
  err2Y += c00;
  err2Z += c11;
  errYZ += c10;
  float_v d = float_v( Vc::One ) / ( err2Y*err2Z - errYZ*errYZ );


  float_m success = mask;//  if ( ISUNLIKELY( err2Y < 1.e-8f ) || ISUNLIKELY( err2Z < 1.e-8f ) ) return 0;
  success &= (err2Y > 1.e-8f ) && ( err2Z > 1.e-8f );

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

    // K = CHtS

  const float_v
    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_v sinPhi = fP[2] + k20 * z0 + k21 * z1;

  success &= CAMath::Abs( sinPhi ) < maxSinPhi;

  fNDF( static_cast<int_m>(success) )  += hitNDF;
  fChi2(success) += mS0 * z0 * z0 + 2 * z0 * z1 * mS1 + mS2 * z1 * z1;
  success &= fChi2 < chi2Cut;
  if ( success.isEmpty() ) return success; // TODO move upper

    // foreach_bit( int i, success ) {
    //   cout << i << " " << mS0[i] << " " << mS1[i] << " " << mS2[i] << " " << (mS0 * z0 * z0)[i] << " " << (2 * z0 * z1 * mS1)[i] << " " << (mS2 * z1 * z1)[i] << " " << fChi2[i] << endl;
    // }
    // std::cout << success << std::endl;
  fP[ 0](success) += k00 * z0 + k01 * z1;
  fP[ 1](success) += k10 * z0 + k11 * z1;
  fP[ 2](success) = sinPhi ;
  fP[ 3](success) += k30 * z0 + k31 * z1;
  fP[ 4](success) += k40 * z0 + k41 * z1;
  fC[ 0](success) -= (k00 * c00 + k01 * c10); //c00
  
  fC[ 1](success) -= (k10 * c00 + k11 * c10); //c10
  fC[ 2](success) -= (k10 * c10 + k11 * c11); //c11
  
  fC[ 3](success) -= (k20 * c00 + k21 * c10); //c20
  fC[ 4](success) -= (k20 * c10 + k21 * c11); //c21
  fC[ 5](success) -= (k20 * c20 + k21 * c21); //c22
  
  fC[ 6](success) -= (k30 * c00 + k31 * c10); //c30
  fC[ 7](success) -= (k30 * c10 + k31 * c11); //c31
  fC[ 8](success) -= (k30 * c20 + k31 * c21); //c32
  fC[ 9](success) -= (k30 * c30 + k31 * c31); //c33

  fC[10](success) -= (k40 * c00 + k41 * c10); //c40
  fC[11](success) -= (k40 * c10 + k41 * c11); //c41
// std::cout << fC[11] << " " << fC[10] <<  std::endl;

  fC[12](success) -= (k40 * c20 + k41 * c21); //c42
  fC[13](success) -= (k40 * c30 + k41 * c31); //c43
  fC[14](success) -= (k40 * c40 + k41 * c41); //c44
  
  return success;
}

float_m PndCATrackParamVector::FilterWithMaterial	(	const float_v &	y,
		const float_v &	z,
		const PndCAStripInfo &	info,
		float_v	err2,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 715 of file PndCATrackParamVector.cxx.

References CAMath::Abs(), c10, c11, c20, c21, PndCAStripInfo::cos, f, fC, fChi2, fNDF, fP, PndCAStripInfo::sin, mrfdata_8b_error::success, and z.

{
  assert( maxSinPhi > 0.f );
  float_m success = mask;
  success &= ( err2 > 1.e-8f );

  const float_v& c00 = fC[0];
  const float_v& c10 = fC[1];
  const float_v& c11 = fC[2];

  const float_v& sb = info.sin;
  const float_v& cb = info.cos;

  const float_v u = cb*y + sb*z;
  const float_v zeta = cb*fP[0] + sb*fP[1] - u;

    // F = CH'
  const float_v  F0 = cb*c00 + sb*c10;
  const float_v  F1 = cb*c10 + sb*c11;
        
  const float_v  HCH = ( F0*cb + F1*sb );

  const float_v  wi = 1/( err2 + HCH );
  const float_v  zetawi = zeta *wi;
  
  fChi2(success) += zeta * zetawi;
  success &= fChi2 < chi2Cut;
  if ( success.isEmpty() ) return success;

  const float_v& c20 = fC[3];
  const float_v& c21 = fC[4];
  const float_v& c30 = fC[6];
  const float_v& c31 = fC[7];
  const float_v& c40 = fC[10];
  const float_v& c41 = fC[11];
  
  const float_v  F2 = cb*c20 + sb*c21;
  const float_v  F3 = cb*c30 + sb*c31;
  const float_v  F4 = cb*c40 + sb*c41;

        
  const float_v  K1 = F1*wi;
  const float_v  K2 = F2*wi;
  const float_v  K3 = F3*wi;
  const float_v  K4 = F4*wi;
 

  const float_v sinPhi = fP[2] - F2*zetawi;

  success &= CAMath::Abs( sinPhi ) < maxSinPhi;

  fNDF( static_cast<int_m>(success) )  += 1; // TODO

  fP[ 0](success) -= F0*zetawi;
  fP[ 1](success) -= F1*zetawi;
  fP[ 2](success) = sinPhi ;
  fP[ 3](success) -= F3*zetawi;
  fP[ 4](success) -= F4*zetawi;

  fC[ 0](success) -= F0*F0*wi;
  
  fC[ 1](success) -= K1*F0;
  fC[ 2](success) -= K1*F1;
  
  fC[ 3](success) -= K2*F0;
  fC[ 4](success) -= K2*F1;
  fC[ 5](success) -= K2*F2;
  
  fC[ 6](success) -= K3*F0;
  fC[ 7](success) -= K3*F1;
  fC[ 8](success) -= K3*F2;
  fC[ 9](success) -= K3*F3;

  fC[10](success) -= K4*F0;
  fC[11](success) -= K4*F1;

  fC[12](success) -= K4*F2;
  fC[13](success) -= K4*F3;
  fC[14](success) -= K4*F4;
  
  return success;
}

float_v PndCATrackParamVector::GetChi2 ( ) const

inline

Definition at line 127 of file PndCATrackParamVector.h.

References fChi2.

Referenced by SetTrackParam().

{ return fChi2; }

float_v PndCATrackParamVector::GetCosPhi ( ) const

inline

Definition at line 132 of file PndCATrackParamVector.h.

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

Referenced by GetDCAPoint(), GetS(), Rotate(), and RotateXY().

{ return fSignCosPhi*CAMath::Sqrt( float_v( Vc::One ) - SinPhi()*SinPhi() ); }

float_v PndCATrackParamVector::GetCosPhiPositive ( ) const

inline

Definition at line 131 of file PndCATrackParamVector.h.

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

{ return CAMath::Sqrt( float_v( Vc::One ) - SinPhi()*SinPhi() ); }

void PndCATrackParamVector::GetDCAPoint	(	const float_v &	x,
		const float_v &	y,
		const float_v &	z,
		float_v *	px,
		float_v *	py,
		float_v *	pz,
		const float_v &	Bz
	)		const

Definition at line 190 of file PndCATrackParamVector.cxx.

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

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

  const float_v &x0 = GetX();
  const float_v &y0 = GetY();
  const float_v &k  = GetKappa( Bz );
  const float_v &ex = GetCosPhi();
  const float_v &ey = GetSinPhi();
  const float_v &dx = x - x0;
  const float_v &dy = y - y0;
  const float_v &ax = dx * k + ey;
  const float_v &ay = dy * k - ex;
  const float_v &a = sqrt( ax * ax + ay * ay );
  *xp = x0 + ( dx - ey * ( ( dx * dx + dy * dy ) * k - 2.f * ( -dx * ey + dy * ex ) ) / ( a + 1.f ) ) / a;
  *yp = y0 + ( dy + ex * ( ( dx * dx + dy * dy ) * k - 2.f * ( -dx * ey + dy * ex ) ) / ( a + 1.f ) ) / a;
  const float_v s = GetS( x, y, Bz );
  *zp = GetZ() + GetDzDs() * s;
  const float_v dZ = CAMath::Abs( GetDzDs() * CAMath::TwoPi() / k );
  const float_m mask = CAMath::Abs( k ) > 1.e-2f && dZ > .1f;
  ( *zp )( mask ) += CAMath::Round( ( z - *zp ) / dZ ) * dZ;
}

float_v PndCATrackParamVector::GetDist2

(

const PndCATrackParamVector &

t

)

const

Definition at line 153 of file PndCATrackParamVector.cxx.

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

{
  // get squared distance between tracks

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

float_v PndCATrackParamVector::GetDistXZ2

(

const PndCATrackParamVector &

t

)

const

Definition at line 163 of file PndCATrackParamVector.cxx.

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

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

  const float_v &dx = GetX() - t.GetX();
  const float_v &dz = GetZ() - t.GetZ();
  return dx*dx + dz*dz;
}

float_v PndCATrackParamVector::GetDzDs ( ) const

inline

Definition at line 124 of file PndCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint().

{ return fP[3]; }

float_v PndCATrackParamVector::GetKappa ( const float_v &  Bz ) const

inline

Definition at line 130 of file PndCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint(), and GetS().

{ return fP[4]*Bz; }

int_v PndCATrackParamVector::GetNDF ( ) const

inline

Definition at line 128 of file PndCATrackParamVector.h.

References fNDF.

Referenced by SetTrackParam().

{ return fNDF; }

float_v PndCATrackParamVector::GetQPt ( ) const

inline

Definition at line 125 of file PndCATrackParamVector.h.

References fP.

{ return fP[4]; }

float_v PndCATrackParamVector::GetS	(	const float_v &	x,
		const float_v &	y,
		const float_v &	Bz
	)		const

Definition at line 173 of file PndCATrackParamVector.cxx.

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

Referenced by GetDCAPoint().

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

  const float_v &k  = GetKappa( Bz );
  const float_v &ex = GetCosPhi();
  const float_v &ey = GetSinPhi();
  const float_v &x = _x - GetX();
  const float_v &y = _y - GetY();
  float_v dS = x * ex + y * ey;
  const float_m &mask = CAMath::Abs( k ) > 1.e-4f;
  if ( !mask.isEmpty() ) {
    dS( mask ) = CAMath::ATan2( k * dS, 1 + k * ( x * ey - y * ex ) ) / k;
  }
  return dS;
}

float_v PndCATrackParamVector::GetSignCosPhi ( ) const

inline

Definition at line 126 of file PndCATrackParamVector.h.

References fSignCosPhi.

{ return fSignCosPhi; }

float_v PndCATrackParamVector::GetSinPhi ( ) const

inline

Definition at line 123 of file PndCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint(), and GetS().

{ return fP[2]; }

float_v PndCATrackParamVector::GetX ( ) const

inline

Definition at line 120 of file PndCATrackParamVector.h.

References fX.

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

{ return fX; }

float_v PndCATrackParamVector::GetY ( ) const

inline

Definition at line 121 of file PndCATrackParamVector.h.

References fP.

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

{ return fP[0]; }

float_v PndCATrackParamVector::GetZ ( ) const

inline

Definition at line 122 of file PndCATrackParamVector.h.

References fP.

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

{ return fP[1]; }

void PndCATrackParamVector::InitByHit	(	const PndCAHitV &	hit,
		const PndCAParam &	param,
		const float_v &	dQMom
	)

Definition at line 105 of file PndCATrackParamVector.cxx.

References PndCAHitV::Angle(), PndCAHitV::Err2X1(), PndCAHitV::Err2X2(), f, SetAngle(), SetChi2(), SetCov(), SetDzDs(), SetNDF(), SetQPt(), SetSignCosPhi(), SetSinPhi(), SetX(), SetY(), SetZ(), PndCAHitV::X0(), PndCAHitV::X1(), and PndCAHitV::X2().

{
  SetX( hit.X0() );
  SetY( hit.X1() );
  SetZ( hit.X2() );
  SetSinPhi( 0.f );
  SetSignCosPhi( 1.f );
  SetDzDs( 0.f );  
  SetQPt( 0.f );
  SetAngle( hit.Angle() );
  
  SetCov( 0, hit.Err2X1() );
  SetCov( 1, 0.f );
  SetCov( 2, hit.Err2X2() );
  SetCov( 3, 0.f );
  SetCov( 4, 0.f );
  SetCov( 5, 1.f );
  SetCov( 6, 0.f );
  SetCov( 7, 0.f );
  SetCov( 8, 0.f );
  SetCov( 9, 1.f );
  SetCov(10, 0.f );
  SetCov(11, 0.f );
  SetCov(12, 0.f );
  SetCov(13, 0.f );
  SetCov(14, dQMom*dQMom );

  SetChi2( 0.f );
    
  SetNDF( -3 );

  UNUSED_PARAM1(param);
  // const int ista0 = hit.IStation();
  // const PndCAStation sta0 = param.Station(ista0);
  
  // L1FieldValue B;
  // sta0.fieldSlice.GetFieldValue( hit.X1(), hit.X2(), B );
  // SetField( 1, B, sta0.z );

  // if ( ista0 > 0 ) {
  //   const L1Station staC = param.Station(ista0/2); // The field will be used to propaget to the target. Therefore take station in a middle between target and hit.
  //   const float_v tx1 = hit.X1()/hit.X0(); // suppose target is at 0
  //   const float_v tx2 = hit.X2()/hit.X0(); // suppose target is at 0
  //   staC.fieldSlice.GetFieldValue( tx1*staC.z, tx2*staC.z, B );
  //   SetField( 0, B, staC.z );
  // }
}

void PndCATrackParamVector::InitByTarget

(

const PndCATarget &

target

)

Definition at line 72 of file PndCATrackParamVector.cxx.

References PndCATarget::Err2QMom(), PndCATarget::Err2X1(), PndCATarget::Err2X2(), f, PndCATarget::NDF(), SetChi2(), SetCov(), SetDzDs(), SetNDF(), SetQPt(), SetSignCosPhi(), SetSinPhi(), SetX(), SetY(), SetZ(), PndCATarget::X0(), PndCATarget::X1(), and PndCATarget::X2().

Referenced by PndCAGBTracker::Create1Plets(), and PndCAFunctionality::FitIteration().

{
  SetX( t.X0() );
  SetY( t.X1() );
  SetZ( t.X2() );
  SetSinPhi( 1.f );
  SetSignCosPhi( 0.f );
  SetDzDs( 0.f );  
  SetQPt( 0.f );
  SetCov( 0.f, t.Err2X1() );
  SetCov( 1, 0.f );
  SetCov( 2, t.Err2X2() );
  SetCov( 3, 0.f );
  SetCov( 4, 0.f );
  SetCov( 5, 1.f );
  SetCov( 6, 0.f );
  SetCov( 7, 0.f );
  SetCov( 8, 0.f );
  SetCov( 9, 1.f );
  SetCov(10, 0.f );
  SetCov(11, 0.f );
  SetCov(12, 0.f );
  SetCov(13, 0.f );
  SetCov(14, t.Err2QMom() );

  SetChi2( 0.f );
    
  SetNDF( -5 + t.NDF() );

  // SetField( 0, t.B(), t.X0() );
}

void PndCATrackParamVector::InitCovMatrix

(

float_v

d2QMom = 0.f

)

Definition at line 49 of file PndCATrackParamVector.cxx.

References f, SetChi2(), SetCov(), and SetNDF().

Referenced by PndCAFunctionality::FitIteration().

{
  SetCov( 0.f, 100.f );
  SetCov( 1, 0.f );
  SetCov( 2, 100.f );
  SetCov( 3, 0.f );
  SetCov( 4, 0.f );
  SetCov( 5, 1.f );
  SetCov( 6, 0.f );
  SetCov( 7, 0.f );
  SetCov( 8, 0.f );
  SetCov( 9, 1.f );
  SetCov(10, 0.f );
  SetCov(11, 0.f );
  SetCov(12, 0.f );
  SetCov(13, 0.f );
  SetCov(14, d2QMom );

  SetChi2( 0.f );
    
  SetNDF( -5 );
}

void PndCATrackParamVector::InitDirection ( float_v  r0, float_v  r1, float_v  r2  )

inline

Definition at line 69 of file PndCATrackParamVector.h.

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

Referenced by PndCAGBTracker::Create1Plets(), and PndCAFunctionality::FitIteration().

    {
      const float_v r = sqrt( r0*r0+r1*r1 );
      SetSinPhi( r1/r );
      SetSignCosPhi( 1.f);// r0/abs(r0) );
      SetDzDs( r2/r );
    }

int_v PndCATrackParamVector::ISec ( ) const

inline

Definition at line 97 of file PndCATrackParamVector.h.

References fISec.

Referenced by PndCAGBTracker::PickUpHits().

{ return fISec; }

int_v PndCATrackParamVector::NDF ( ) const

inline

Definition at line 112 of file PndCATrackParamVector.h.

References fNDF.

Referenced by Accept(), and Filter().

{ return fNDF; }

const float_v& PndCATrackParamVector::Par ( int  i ) const

inline

Definition at line 140 of file PndCATrackParamVector.h.

References fP, and i.

Referenced by PndCATrackParam::PndCATrackParam(), PndFTSCATrackParam::PndFTSCATrackParam(), and SetTrackParam().

{ return fP[i]; }

const float_v* PndCATrackParamVector::Par ( ) const

inlineprivate

Definition at line 145 of file PndCATrackParamVector.h.

References fP.

{ return fP; }

float_v* PndCATrackParamVector::Par ( )

inlineprivate

Definition at line 147 of file PndCATrackParamVector.h.

References fP.

{ return fP; }

float_v PndCATrackParamVector::QMomentum ( ) const

inline

Definition at line 104 of file PndCATrackParamVector.h.

References QPt().

Referenced by PndCAFunctionality::Fit(), and PndCAFunctionality::FitIteration().

{ return QPt(); } // used for triplets comparison

float_v PndCATrackParamVector::QPt ( ) const

inline

Definition at line 94 of file PndCATrackParamVector.h.

References fP.

Referenced by QMomentum(), and TransportToX().

{ return fP[4]; }

void PndCATrackParamVector::Reset ( )

inline

Definition at line 54 of file PndCATrackParamVector.h.

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

              {
    fX.setZero();
    fSignCosPhi = float_v(1.f);
    fChi2.setZero();
    fNDF.setZero();
    //fISec.setZero();
    {
      for ( int i = 0; i <  5; ++i ) fP[i].setZero();
      for ( int i = 0; i < 15; ++i ) fC[i].setZero();
    }
  }

float_m PndCATrackParamVector::Rotate	(	const float_v &	alpha,
		PndCATrackLinearisationVector &	t0,
		const float	maxSinPhi = .999f,
		const float_m &	mask = float_m( true )
	)

Definition at line 568 of file PndCATrackParamVector.cxx.

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

Referenced by PndCAGBTracker::PickUpHits(), Transport(), and Transport0().

{
  //* Rotate the coordinate system in XY on the angle alpha
  if ( (CAMath::Abs(alpha) < 1e-6f || !mask).isFull() ) return mask;
  
  const float_v cA = CAMath::Cos( alpha );
  const float_v sA = CAMath::Sin( alpha );
  const float_v x0 = X(), y0 = Y(), sP = t0.SinPhi(), cP = t0.CosPhi();
  const float_v cosPhi = cP * cA + sP * sA;
  const float_v sinPhi = -cP * sA + sP * cA;

  float_m ReturnMask(mask);
  ReturnMask &= (!( CAMath::Abs( sinPhi ) > maxSinPhi || CAMath::Abs( cosPhi ) < 1.e-4f || CAMath::Abs( cP ) < 1.e-4f  ));

  float_v tmp = alpha*0.15915494f;// 1/(2.f*3.1415f);
  ReturnMask &= abs(tmp - round(tmp)) < 0.167f; // allow turn by 60 degree only TODO scalar

  //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_v j0 = cP / cosPhi;
  const float_v j2 = cosPhi / cP;
  const float_v d = SinPhi() - sP;

  SetX( x0*cA +  y0*sA, ReturnMask );
  SetY(-x0*sA +  y0*cA, ReturnMask );
  t0.SetCosPhi( cosPhi );
  t0.SetSinPhi( sinPhi );

  SetSinPhi( sinPhi + j2*d, ReturnMask );
  
  fC[0](ReturnMask) *= j0 * j0;
  fC[1](ReturnMask) *= j0;
  fC[3](ReturnMask) *= j0;
  fC[6](ReturnMask) *= j0;
  fC[10](ReturnMask) *= j0;

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

  fAlpha(ReturnMask) += alpha;
  
  return ReturnMask;
}

float_m PndCATrackParamVector::Rotate ( const float_v &  alpha, const float  maxSinPhi = .999f, const float_m &  mask = float_m( true )  )

inline

Definition at line 392 of file PndCATrackParamVector.h.

References CAMath::Abs(), alpha, CAMath::Cos(), d, f, fAlpha, fC, GetCosPhi(), SetSinPhi(), SetX(), SetY(), CAMath::Sin(), SinPhi(), X(), Y(), and y0.

{
  // Rotate the coordinate system in XY on the angle alpha
  if ( (CAMath::Abs(alpha) < 1e-6f || !mask).isFull() ) return mask;
  
  const float_v cA = CAMath::Cos( alpha );
  const float_v sA = CAMath::Sin( alpha );
  const float_v x0 = X(), y0 = Y(), sP = SinPhi(), cP = GetCosPhi();
  const float_v cosPhi = cP * cA + sP * sA;
  const float_v sinPhi = -cP * sA + sP * cA;

  float_m ReturnMask(mask);
  ReturnMask &= (!( CAMath::Abs( sinPhi ) > maxSinPhi || CAMath::Abs( cosPhi ) < 1.e-4f || CAMath::Abs( cP ) < 1.e-4f  ));

  float_v tmp = alpha*0.15915494f;// 1/(2.f*3.1415f);
  ReturnMask &= abs(tmp - round(tmp)) < 0.167f; // allow turn by 60 degree only TODO scalar

  //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_v j0 = cP / cosPhi;
  const float_v j2 = cosPhi / cP;
  const float_v d = SinPhi() - sP;

  SetX( x0*cA +  y0*sA, ReturnMask );
  SetY(-x0*sA +  y0*cA, ReturnMask );  

  SetSinPhi( sinPhi + j2*d, ReturnMask );
  
  fC[0](ReturnMask) *= j0 * j0;
  fC[1](ReturnMask) *= j0;
  fC[3](ReturnMask) *= j0;
  fC[6](ReturnMask) *= j0;
  fC[10](ReturnMask) *= j0;

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

  fAlpha(ReturnMask) += alpha;
  
  return ReturnMask;
}

void PndCATrackParamVector::RotateXY ( float_v  alpha, float_v &  x, float_v &  y, float_v &  sin, const float_m &  mask = float_m( true )  ) const

inline

Definition at line 442 of file PndCATrackParamVector.h.

References CAMath::Cos(), f, GetCosPhi(), CAMath::Sin(), sin(), SinPhi(), X(), x, Y(), and y.

{
  //* Rotate the coordinate system in XY on the angle alpha
  if ( (abs(alpha) < 1e-6f || !mask).isFull() ) return;
  
  const float_v cA = CAMath::Cos( alpha );
  const float_v sA = CAMath::Sin( alpha );

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

void PndCATrackParamVector::SetAngle ( const float_v &  v )

inline

Definition at line 200 of file PndCATrackParamVector.h.

References fAlpha, and v.

Referenced by PndCAGBTracker::Create1Plets(), PndCAFunctionality::FitIteration(), and InitByHit().

{  fAlpha = v; }

void PndCATrackParamVector::SetAngle ( const float_v &  v, const float_m &  m  )

inline

Definition at line 201 of file PndCATrackParamVector.h.

References fAlpha, and v.

{  fAlpha(m) = v; }

void PndCATrackParamVector::SetChi2 ( const float_v &  v )

inline

Definition at line 194 of file PndCATrackParamVector.h.

References fChi2, and v.

Referenced by ConvertTrackParamToVector(), InitByHit(), InitByTarget(), and InitCovMatrix().

{  fChi2 = v; }

void PndCATrackParamVector::SetChi2 ( const float_v &  v, const float_m &  m  )

inline

Definition at line 195 of file PndCATrackParamVector.h.

References fChi2, and v.

{  fChi2(m) = v; }

void PndCATrackParamVector::SetCov ( int  i, const float_v &  v  )

inline

Definition at line 176 of file PndCATrackParamVector.h.

References fC, i, and v.

Referenced by ConvertTrackParamToVector(), InitByHit(), InitByTarget(), and InitCovMatrix().

{ fC[i] = v; }

void PndCATrackParamVector::SetCov ( int  i, const float_v &  v, const float_m &  m  )

inline

Definition at line 177 of file PndCATrackParamVector.h.

References fC, i, and m.

{ fC[i]( m ) = v; }

void PndCATrackParamVector::SetDzDs ( const float_v &  v )

inline

Definition at line 187 of file PndCATrackParamVector.h.

References fP, and v.

Referenced by InitByHit(), InitByTarget(), and InitDirection().

{  fP[3] = v; }

void PndCATrackParamVector::SetDzDs ( const float_v &  v, const float_m &  m  )

inline

Definition at line 188 of file PndCATrackParamVector.h.

References fP, and m.

{  fP[3]( m ) = v; }

void PndCATrackParamVector::SetErr2QPt ( float_v  v )

inline

Definition at line 208 of file PndCATrackParamVector.h.

References fC, and v.

{ fC[14] = v; }

void PndCATrackParamVector::SetErr2Y ( float_v  v )

inline

Definition at line 206 of file PndCATrackParamVector.h.

References fC, and v.

{ fC[0] = v; }

void PndCATrackParamVector::SetErr2Z ( float_v  v )

inline

Definition at line 207 of file PndCATrackParamVector.h.

References fC, and v.

{ fC[2] = v; }

void PndCATrackParamVector::SetISec ( const int_v &  v )

inline

Definition at line 203 of file PndCATrackParamVector.h.

References fISec, and v.

Referenced by PndCAGBTracker::Create1Plets().

{  fISec = v; }

void PndCATrackParamVector::SetISec ( const int_v &  v, const int_m &  m  )

inline

Definition at line 204 of file PndCATrackParamVector.h.

References fISec, and v.

{  fISec(m) = v; }

void PndCATrackParamVector::SetNDF ( int  v )

inline

Definition at line 196 of file PndCATrackParamVector.h.

References fNDF, and v.

Referenced by ConvertTrackParamToVector(), InitByHit(), InitByTarget(), and InitCovMatrix().

{ fNDF = v; }

void PndCATrackParamVector::SetNDF ( const int_v &  v )

inline

Definition at line 197 of file PndCATrackParamVector.h.

References fNDF, and v.

{ fNDF = v; }

void PndCATrackParamVector::SetNDF ( const int_v &  v, const int_m &  m  )

inline

Definition at line 198 of file PndCATrackParamVector.h.

References fNDF, and v.

{ fNDF(m) = v; }

void PndCATrackParamVector::SetPar ( int  i, const float_v &  v  )

inline

Definition at line 174 of file PndCATrackParamVector.h.

References fP, i, and v.

Referenced by ConvertTrackParamToVector().

{ fP[i] = v; }

void PndCATrackParamVector::SetPar ( int  i, const float_v &  v, const float_m &  m  )

inline

Definition at line 175 of file PndCATrackParamVector.h.

References fP, i, and m.

{ fP[i]( m ) = v; }

void PndCATrackParamVector::SetQMomentum ( const float_v &  v )

inline

Definition at line 190 of file PndCATrackParamVector.h.

References SetQPt().

Referenced by PndCAFunctionality::FitIteration().

{  SetQPt(v); }

void PndCATrackParamVector::SetQPt ( const float_v &  v )

inline

Definition at line 189 of file PndCATrackParamVector.h.

References fP, and v.

Referenced by InitByHit(), InitByTarget(), and SetQMomentum().

{  fP[4] = v; }

void PndCATrackParamVector::SetQPt ( const float_v &  v, const float_m &  m  )

inline

Definition at line 191 of file PndCATrackParamVector.h.

References fP, and m.

{  fP[4]( m ) = v; }

void PndCATrackParamVector::SetSignCosPhi ( const float_v &  v )

inline

Definition at line 192 of file PndCATrackParamVector.h.

References fSignCosPhi, and v.

Referenced by ConvertTrackParamToVector(), InitByHit(), InitByTarget(), and InitDirection().

{  fSignCosPhi = v; }

void PndCATrackParamVector::SetSignCosPhi ( const float_v &  v, const float_m &  m  )

inline

Definition at line 193 of file PndCATrackParamVector.h.

References fSignCosPhi, and v.

{  fSignCosPhi(m) = v; }

void PndCATrackParamVector::SetSinPhi ( const float_v &  v )

inline

Definition at line 185 of file PndCATrackParamVector.h.

References fP, and v.

Referenced by InitByHit(), InitByTarget(), InitDirection(), and Rotate().

{  fP[2] = v; }

void PndCATrackParamVector::SetSinPhi ( const float_v &  v, const float_m &  m  )

inline

Definition at line 186 of file PndCATrackParamVector.h.

References fP, and m.

{  fP[2]( m ) = v; }

void PndCATrackParamVector::SetTrackParam ( const PndCATrackParamVector &  param, const float_m &  m = float_m( true )  )

inline

Definition at line 151 of file PndCATrackParamVector.h.

References Cov(), fC, fChi2, fISec, fNDF, fP, fSignCosPhi, fX, GetChi2(), GetNDF(), i, m, Par(), SignCosPhi(), and X().

Referenced by PndCAGBTracker::PickUpHits().

  {
    for(int i=0; i<5; i++) fP[i](m) = param.Par()[i];
    for(int i=0; i<15; i++) fC[i](m) = param.Cov()[i];
    fX(m) = param.X();
    fSignCosPhi(m) = param.SignCosPhi();
    fChi2(m) = param.GetChi2();
    fNDF(static_cast<int_m>(m)) = param.GetNDF();
    fISec(m) = param.fISec;
  }

void PndCATrackParamVector::SetTrackParam ( const PndCATrackParamVector &  p, int  k, int  m  )

inline

Definition at line 162 of file PndCATrackParamVector.h.

References Cov(), f, fAlpha, fC, fChi2, fISec, fNDF, fP, fSignCosPhi, fX, i, m, and Par().

  {
    fX[k]=p.fX[m];
    fSignCosPhi = float_v(1.f);
    for(int i=0; i<5; i++) fP[i][k] = p.Par()[i][m];
    for(int i=0; i<15; i++) fC[i][k] = p.Cov()[i][m];
    fChi2[k] = p.fChi2[m];
    fNDF[k] = p.fNDF[m];
    fAlpha[k] = p.fAlpha[m];
    fISec[k] = p.fISec[m];
  }

void PndCATrackParamVector::SetX ( const float_v &  v )

inline

Definition at line 179 of file PndCATrackParamVector.h.

References fX, and v.

Referenced by ConvertTrackParamToVector(), InitByHit(), InitByTarget(), and Rotate().

{  fX = v;    }

void PndCATrackParamVector::SetX ( const float_v &  v, const float_m &  m  )

inline

Definition at line 182 of file PndCATrackParamVector.h.

References fX, and v.

{  fX( m ) = v;    }

void PndCATrackParamVector::SetY ( const float_v &  v )

inline

Definition at line 180 of file PndCATrackParamVector.h.

References fP, and v.

Referenced by InitByHit(), InitByTarget(), and Rotate().

{  fP[0] = v; }

void PndCATrackParamVector::SetY ( const float_v &  v, const float_m &  m  )

inline

Definition at line 183 of file PndCATrackParamVector.h.

References fP, and m.

{  fP[0]( m ) = v; }

void PndCATrackParamVector::SetZ ( const float_v &  v )

inline

Definition at line 181 of file PndCATrackParamVector.h.

References fP, and v.

Referenced by InitByHit(), and InitByTarget().

{  fP[1] = v; }

void PndCATrackParamVector::SetZ ( const float_v &  v, const float_m &  m  )

inline

Definition at line 184 of file PndCATrackParamVector.h.

References fP, and m.

{  fP[1]( m ) = v; }

float_v PndCATrackParamVector::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 110 of file PndCATrackParamVector.h.

References fSignCosPhi.

Referenced by PndCATrackLinearisationVector::PndCATrackLinearisationVector(), SetTrackParam(), Tx1(), and Tx2().

{ return fSignCosPhi; }

float_v PndCATrackParamVector::SinPhi ( ) const

inline

Definition at line 92 of file PndCATrackParamVector.h.

References fP.

Referenced by GetCosPhi(), GetCosPhiPositive(), PndCAGBTracker::PickUpHits(), Rotate(), RotateXY(), Transport0ToX(), TransportJ0ToX0(), TransportToX(), Tx1(), Tx2(), and PndCAHitV::X1Corrected().

{ return fP[2]; }

float_m PndCATrackParamVector::Transport	(	const int_v &	ista,
		const PndCAParam &	param,
		const float_m &	mask = float_m( true )
	)

Definition at line 800 of file PndCATrackParamVector.cxx.

References PndCAParam::cBz(), PndCAParam::GetR(), and TransportToX().

Referenced by PndCAFunctionality::FitIteration().

{
  float_m active = mask;
  PndCATrackLinearisationVector tE( *this );
  active &= TransportToX( param.GetR( ista, active ), tE, param.cBz(), 0.999f, NULL, active );
  //SG!!
  /*
  PndCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  PndCATrackLinearisationVector tE( *this );
  active &= TransportToXWithMaterial( param.GetR( ista, active ), tE, fitPar, param.GetXOverX0(int_v(ista),active), param.GetXTimesRho(int_v(ista),active), param.cBz(), 0.999f, active );
  */
  return active;
}

float_m PndCATrackParamVector::Transport	(	const PndCAHitV &	hit,
		const PndCAParam &	p,
		const float_m &	mask = float_m( true )
	)

Definition at line 815 of file PndCATrackParamVector.cxx.

References CAMath::Abs(), PndCAHitV::Angle(), PndCAParam::cBz(), f, fAlpha, fX, PndCAHitV::IsValid(), Rotate(), TransportToX(), and PndCAHitV::X0().

{
  if ( ((CAMath::Abs(-fAlpha + hit.Angle()) < 1e-6f && CAMath::Abs(-fX + hit.X0()) < 2e-4f) || !mask).isFull() ) return mask;
  float_m active = mask & hit.IsValid();
  PndCATrackLinearisationVector tR( *this );
  active &= Rotate(  -fAlpha + hit.Angle(), tR, .999f, active ); 
  active &= TransportToX( hit.X0(), tR,  param.cBz(), 0.999f, NULL, active );

  //SG!!
  /* 
  PndCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  active &= TransportToXWithMaterial( hit.X0(), tR, fitPar, param.GetXOverX0(hit.IStations(),active), param.GetXTimesRho(hit.IStations(),active), param.cBz(), 0.999f, active );
  */
  return active;
}

float_m PndCATrackParamVector::Transport	(	const PndCAHit &	hit,
		const PndCAParam &	p,
		const float_m &	mask = float_m( true )
	)

Definition at line 844 of file PndCATrackParamVector.cxx.

References CAMath::Abs(), PndCAHit::Angle(), PndCAParam::cBz(), f, fAlpha, fX, Rotate(), TransportToX(), and PndCAHit::X0().

{
  if ( ((CAMath::Abs(-fAlpha + hit.Angle()) < 1e-6f && CAMath::Abs(-fX + hit.X0()) < 2e-4f) || !mask).isFull() ) return mask;
  
  float_m active = mask;

  PndCATrackLinearisationVector tR( *this );
  active &= Rotate(  -fAlpha + hit.Angle(), tR, .999f, active );

  active &= TransportToX( hit.X0(), tR,  param.cBz(), 0.999f, NULL, active );
  /*
  PndCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  active &= TransportToXWithMaterial( hit.X0(), tR, fitPar, param.GetXOverX0(hit.IStation()), param.GetXTimesRho(hit.IStation()), param.cBz(), 0.999f, active );
  */
  return active;
}

float_m PndCATrackParamVector::Transport0 ( const int_v &  ista, const PndCAParam &  param, const float_m &  mask = float_m( true )  )

inline

Definition at line 565 of file PndCATrackParamVector.h.

References PndCAParam::cBz(), PndCAParam::GetR(), and Transport0ToX().

Referenced by PndCAGBTracker::PickUpHits().

{
  float_m active = mask;
  //PndCATrackLinearisationVector tE( *this );
  //active &= TransportToX( param.GetR( ista, active ), tE, param.cBz(), 0.999f, NULL, active );
  active &= Transport0ToX( param.GetR( ista, active ), param.cBz(), active);
  return active;
}

float_m PndCATrackParamVector::Transport0	(	const PndCAHitV &	hit,
		const PndCAParam &	param,
		const float_m &	mask = float_m( true )
	)

Definition at line 1014 of file PndCATrackParamVector.cxx.

References PndCAHitV::Angle(), PndCAParam::cBz(), fAlpha, PndCAHitV::IsValid(), Rotate(), Transport0ToX(), and PndCAHitV::X0().

{
  //if ( ((CAMath::Abs(-fAlpha + hit.Angle()) < 1e-6f && CAMath::Abs(-fX + hit.X0()) < 2e-4f) || !mask).isFull() ) return mask;
  float_m active = mask & hit.IsValid();
  active &= Rotate(  -fAlpha + hit.Angle(), .999f, active ); 
  active &= Transport0ToX( hit.X0(),  param.cBz(), active );
  return active;
}

float_m PndCATrackParamVector::Transport0 ( const PndCAHit &  hit, const PndCAParam &  p, const float_m &  mask = float_m( true )  )

inline

Definition at line 574 of file PndCATrackParamVector.h.

References CAMath::Abs(), PndCAHit::Angle(), PndCAParam::cBz(), f, fAlpha, fX, Rotate(), Transport0ToX(), and PndCAHit::X0().

{
  if ( ((CAMath::Abs(-fAlpha + hit.Angle()) < 1e-6f && CAMath::Abs(-fX + hit.X0()) < 2e-4f) || !mask).isFull() ) return mask;  
  float_m active = mask;
  //PndCATrackLinearisationVector tR( *this );
  active &= Rotate(  -fAlpha + hit.Angle(), .999f, active );
  //active &= Rotate(  -fAlpha + hit.Angle(), tR, .999f, active );  
  //active &= TransportToX( hit.X0(), tR,  param.cBz(), 0.999f, NULL, active );
  active &=Transport0ToX( hit.X0(),  param.cBz(), active );
  return active;
}

float_m PndCATrackParamVector::Transport0ToX	(	const float_v &	x,
		const float_v &	Bz,
		const float_m &	mask
	)

Definition at line 953 of file PndCATrackParamVector.cxx.

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

Referenced by PndCAGBTracker::Create1Plets(), and Transport0().

{
  //* 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_v ey = SinPhi();
  const float_v ex = CAMath::Sqrt( float_v( Vc::One ) - ey*ey );
  const float_v dx = x - X();
 
  float_m mask = _mask && CAMath::Abs( ey ) <= .999f;

  const float_v exi = 1.f / ex;
  const float_v dS = dx * exi;
  const float_v h2 = dx * exi * exi *exi;
  const float_v h4 = 0.5f*dx* Bz * h2;
  const float_v dxBz = dx * Bz;

  fX    += dx;
  fP[0] += dS*ey + h4 * fP[4];
  fP[1] += dS * DzDs();
  fP[2] += dxBz * fP[4];

  const float_v c00 = fC[0];
  const float_v c10 = fC[1];
  const float_v c11 = fC[2];
  const float_v c20 = fC[3];
  const float_v c21 = fC[4];
  const float_v c22 = fC[5];
  const float_v c30 = fC[6];
  const float_v c31 = fC[7];
  const float_v c32 = fC[8];
  const float_v c33 = fC[9];
  const float_v c40 = fC[10];
  const float_v c41 = fC[11];
  const float_v c42 = fC[12];
  const float_v c43 = fC[13];
  const float_v 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[10] = c40 + h2 * c42 + h4 * c44;
  fC[11] = c41 + dS * c43;
  fC[12] = c42 + dxBz * c44;
  return mask;
}

float_m PndCATrackParamVector::TransportJ0ToX0 ( const float_v &  x0, const float_v &  cBz, float_v &  extrDy, float_v &  extrDz, float_v  J[], const float_m &  active = float_m( true )  )

inline

Definition at line 533 of file PndCATrackParamVector.h.

References CAMath::Abs(), dx, DzDs(), h2, h4, One, CAMath::RSqrt(), SinPhi(), and X().

Referenced by AddTarget().

{
  const float_v &ey = SinPhi();
  const float_v &dx = x - X();
  const float_v &exi = CAMath::RSqrt( float_v( Vc::One ) - ey * ey ); // RSqrt

  const float_v &dxBz = dx * cBz;
  const float_v &dS = dx * exi;
  const float_v &h2 = dS * exi * exi;
  const float_v &h4 = .5f * h2 * dxBz;

  float_m mask = active && CAMath::Abs( exi ) <= 1.e4f;

  extrDy(mask) = dS*ey;
  extrDz(mask) = dS*DzDs();
  J[0](mask) = dS;
  J[1](mask) = 0;
  J[2](mask) = h4;
  J[3](mask) = dS;
  J[4](mask) = dS;
  J[5](mask) = 0;
  return mask;
}

float_m PndCATrackParamVector::TransportToX	(	const float_v &	x,
		const float_v &	Bz,
		const float	maxSinPhi = .999f,
		const float_m &	mask = float_m( true )
	)

Definition at line 333 of file PndCATrackParamVector.cxx.

References t0.

Referenced by Transport(), and TransportToXWithMaterial().

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

//  assert( ( x == 0 && mask ).isEmpty() );

  PndCATrackLinearisationVector t0( *this );

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

float_m PndCATrackParamVector::TransportToX	(	const float_v &	x,
		PndCATrackLinearisationVector &	t0,
		const float_v &	Bz,
		const float	maxSinPhi = .999f,
		float_v *	DL = 0,
		const float_m &	mask = float_m( true )
	)

Definition at line 220 of file PndCATrackParamVector.cxx.

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

{
  //* 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
  //*
// std::cout << t0.QPt()[0] << " " << fP[4][0] << std::endl;
  const float_v ex = t0.CosPhi();
  const float_v ey = t0.SinPhi();
  const float_v k  = t0.QPt() * Bz;
  const float_v dx = x - X();

  float_v ey1 = k * dx + ey;

  // check for intersection with X=x

  float_v ex1 = CAMath::Sqrt( 1.f - ey1 * ey1 );
  ex1( ex < Vc::Zero ) = -ex1;

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

  const float_v cci = 1.f / cc;
  const float_v exi = 1.f / ex;
  const float_v ex1i = 1.f / ex1;

  float_m mask = _mask && CAMath::Abs( ey1 ) <= maxSinPhi && CAMath::Abs( cc ) >= 1.e-4f && CAMath::Abs( ex ) >= 1.e-4f && CAMath::Abs( ex1 ) >= 1.e-4f;

  const float_v tg = ss * cci; // tan((phi1+phi)/2)

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

  dl( cc < Vc::Zero ) = -dl;
  const float_v dSin = CAMath::Max( float_v( -1.f ), CAMath::Min( float_v( Vc::One ), dl * k * 0.5f ) );
  const float_v dS = ( CAMath::Abs( k ) > 1.e-4f )  ? ( 2 * CAMath::ASin( dSin ) / k ) : dl;
  const float_v dz = dS * t0.DzDs();

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

  const float_v 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_v h2 = dx * ( 1.f + ey * ey1 + ex * ex1 ) * exi * ex1i * cci;
  const float_v h4 = dx2 * ( cc + ss * ey1 * ex1i ) * cci * cci * Bz;
  const float_v dxBz = dx * Bz;

  ex1( !mask ) = ex;
  ey1( !mask ) = ey;

// std::cout << "QPt  " << t0.QPt()[0] << " " << fP[4][0] << 
//              "  CosPhi " << t0.CosPhi()[0] << " " << ex1[0] << std::endl;

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

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

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

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

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

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

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

  fC[10]( mask ) = c40 + h2 * c42 + h4 * c44;
  fC[11]( mask ) = c41 + dS * c43;
  fC[12]( mask ) = c42 + dxBz * c44;
  fC[13]( mask ) = c43;
  fC[14]( mask ) = c44;
//std::cout << fC[10] << " " << fC[11] <<" " << h2 << " " << h4 << " " << c44 << " " << mask << " " << Bz << std::endl;
  //debugKF() << mask << "\n" << *this << std::endl;
  return mask;
}

float_m PndCATrackParamVector::TransportToX ( const float_v &  x, const float_v &  sinPhi0, const float_v &  Bz, const float_v  maxSinPhi = .999f, const float_m &  mask = float_m( true )  )

inline

mvz start 23.01.2010

mvz end 23.01.2010

Definition at line 313 of file PndCATrackParamVector.h.

References CAMath::Abs(), c20, c22, dx, DzDs(), f, fC, fP, fX, h2, h4, One, QPt(), CAMath::RSqrt(), 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
  //*

  //debugKF() << "Start TransportToX(" << x << ", " << _mask << ")\n" << *this << std::endl;

  const float_v &ey = sinPhi0;
  const float_v &dx = x - X();
  const float_v &exi = float_v( Vc::One ) * CAMath::RSqrt( float_v( Vc::One ) - ey * ey ); // RSqrt

  const float_v &dxBz = dx * Bz;
  const float_v &dS = dx * exi;
  const float_v &h2 = dS * exi * exi;
  const float_v &h4 = .5f * h2 * dxBz;
//  const float_v &sinPhi = SinPhi() * (float_v( Vc::One ) - 0.5f * dxBz * QPt() *dxBz * QPt()/ ( float_v( Vc::One ) - SinPhi()*SinPhi() )) + dxBz * QPt();
  const float_v &sinPhi = SinPhi() + dxBz * QPt();

  float_m mask = _mask && CAMath::Abs( exi ) <= 1.e4f;
  mask &= ( (CAMath::Abs( sinPhi ) <= maxSinPhi) || (maxSinPhi <= 0.f) );


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


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

  const float_v two( 2.f );

  fC[0] ( mask ) += h2 * h2 * c22 + h4 * h4 * c44
                     + two * ( h2 * c20 + h4 * ( c40 + h2 * c42 ) );

  //fC[1] ( mask ) += h2 * c21 + h4 * c41 + dS * ( c30 + h2 * c32 + h4 * c43 );
  fC[2] ( mask ) += dS * ( two * c31 + dS * c33 );

  fC[3] ( mask ) += h2 * c22 + h4 * c42 + dxBz * ( c40 + h2 * c42 + h4 * c44 );
  //fC[4] ( mask ) += dS * c32 + dxBz * ( c41 + dS * c43 );
  const float_v &dxBz_c44 = dxBz * c44;
  fC[12]( mask ) += dxBz_c44;
  fC[5] ( mask ) += dxBz * ( two * c42 + dxBz_c44 );

  //fC[6] ( mask ) += h2 * c32 + h4 * c43;
  fC[7] ( mask ) += dS * c33;
  //fC[8] ( mask ) += dxBz * c43;
  //fC[9] ( mask ) = c33;

  fC[10]( mask ) += h2 * c42 + h4 * c44;
  //fC[11]( mask ) += dS * c43;
  //fC[13]( mask ) = c43;
  //fC[14]( mask ) = c44;

  //debugKF() << mask << "\n" << *this << std::endl;
  return mask;
}

float_m PndCATrackParamVector::TransportToXWithMaterial	(	const float_v &	x,
		const float_v &	XOverX0,
		const float_v &	XThimesRho,
		const float_v &	Bz,
		const float	maxSinPhi = .999f
	)

Definition at line 381 of file PndCATrackParamVector.cxx.

References CalculateFitParameters(), and par.

Referenced by TransportToXWithMaterial().

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

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

float_m PndCATrackParamVector::TransportToXWithMaterial	(	const float_v &	x,
		PndCATrackLinearisationVector &	t0,
		PndCATrackFitParam &	par,
		const float_v &	XOverX0,
		const float_v &	XThimesRho,
		const float_v &	Bz,
		const float	maxSinPhi = .999f,
		const float_m &	mask = float_m( true )
	)

Definition at line 347 of file PndCATrackParamVector.cxx.

References CorrectForMeanMaterial(), PndCATrackLinearisationVector::DzDs(), PndCATrackLinearisationVector::SinPhi(), sqrt(), and TransportToX().

{
  //* Transport the track parameters to X=x  taking into account material budget
//   UNUSED_PARAM1(par);
  // const float kRho = 1.025e-3f;//0.9e-3;
//  const float kRadLen = 29.532f;//28.94;
//  const float kRhoOverRadLen = kRho / kRadLen;
  // const float kRhoOverRadLen = 7.68e-5;
  float_v dl;

  float_m mask = mask_ && TransportToX( x, t0, Bz,  maxSinPhi, &dl, mask_ );
//   if ( !mask.isEmpty() ) {
//     CorrectForMeanMaterial( dl * kRhoOverRadLen, dl * kRho, par, mask );
//   }

  const float_v dzds2 = t0.DzDs()*t0.DzDs();
  const float_v sphi2 = t0.SinPhi()*t0.SinPhi();
  float_v koeff = sqrt( (1 + dzds2)*( 1 + sphi2 ) ); // = dr/dx
  mask &= CorrectForMeanMaterial( XOverX0*koeff, XThimesRho*koeff, par, mask );

  return mask;
}

float_m PndCATrackParamVector::TransportToXWithMaterial	(	const float_v &	x,
		PndCATrackFitParam &	par,
		const float_v &	XOverX0,
		const float_v &	XThimesRho,
		const float_v &	Bz,
		const float	maxSinPhi = .999f
	)

Definition at line 372 of file PndCATrackParamVector.cxx.

References t0, and TransportToXWithMaterial().

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

  PndCATrackLinearisationVector t0( *this );
  return TransportToXWithMaterial( x, t0, par, XOverX0, XThimesRho, Bz, maxSinPhi );
}

float_v PndCATrackParamVector::Tx1 ( ) const

inline

Definition at line 102 of file PndCATrackParamVector.h.

References SignCosPhi(), SinPhi(), and sqrt().

{ return SinPhi()/(SignCosPhi()*sqrt( 1 - SinPhi()*SinPhi() )); } // CHECKME

float_v PndCATrackParamVector::Tx2 ( ) const

inline

Definition at line 103 of file PndCATrackParamVector.h.

References DzDs(), SignCosPhi(), SinPhi(), and sqrt().

{ return DzDs()/(SignCosPhi()*sqrt( 1 - SinPhi()*SinPhi() )); } // dx2/dx0 = dz/dx

float_v PndCATrackParamVector::X ( ) const

inline

Definition at line 89 of file PndCATrackParamVector.h.

References fX.

Referenced by Rotate(), RotateXY(), SetTrackParam(), Transport0ToX(), TransportJ0ToX0(), TransportToX(), and X0().

{ return fX;    }

float_v PndCATrackParamVector::X0 ( ) const

inline

Definition at line 99 of file PndCATrackParamVector.h.

References X().

{ return X(); }

float_v PndCATrackParamVector::X1 ( ) const

inline

Definition at line 100 of file PndCATrackParamVector.h.

References Y().

Referenced by PndCAGBTracker::PickUpHits().

{ return Y(); }

float_v PndCATrackParamVector::X2 ( ) const

inline

Definition at line 101 of file PndCATrackParamVector.h.

References Z().

Referenced by PndCAGBTracker::PickUpHits().

{ return Z(); }

float_v PndCATrackParamVector::Y ( ) const

inline

Definition at line 90 of file PndCATrackParamVector.h.

References fP.

Referenced by FilterVtx(), PndCAGBTracker::PickUpHits(), Rotate(), RotateXY(), TransportToX(), and X1().

{ return fP[0]; }

float_v PndCATrackParamVector::Z ( ) const

inline

Definition at line 91 of file PndCATrackParamVector.h.

References fP.

Referenced by FilterVtx(), PndCAGBTracker::PickUpHits(), TransportToX(), and X2().

{ return fP[1]; }

Friends And Related Function Documentation

friend class PndCATrackParam

friend

Definition at line 144 of file PndCATrackParamVector.h.

Member Data Documentation

float_v PndCATrackParamVector::fAlpha

Definition at line 307 of file PndCATrackParamVector.h.

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

float_v PndCATrackParamVector::fC[15]

Definition at line 303 of file PndCATrackParamVector.h.

Referenced by AcceptWithMaterial(), CorrectForMeanMaterial(), Cov(), Err2DzDs(), Err2QMomentum(), Err2QPt(), Err2SinPhi(), Err2X1(), Err2X2(), Err2Y(), Err2Z(), FilterVtx(), FilterWithMaterial(), PndCAGBTracker::PickUpHits(), PndCATrackParamVector(), Reset(), Rotate(), SetCov(), SetErr2QPt(), SetErr2Y(), SetErr2Z(), SetTrackParam(), Transport0ToX(), and TransportToX().

float_v PndCATrackParamVector::fChi2

Definition at line 304 of file PndCATrackParamVector.h.

Referenced by AcceptWithMaterial(), Chi2(), FilterVtx(), FilterWithMaterial(), GetChi2(), PndCAGBTracker::PickUpHits(), Reset(), SetChi2(), and SetTrackParam().

int_v PndCATrackParamVector::fISec

Definition at line 308 of file PndCATrackParamVector.h.

Referenced by ISec(), PndCAGBTracker::PickUpHits(), SetISec(), and SetTrackParam().

int_v PndCATrackParamVector::fNDF

Definition at line 305 of file PndCATrackParamVector.h.

Referenced by AddTarget(), FilterWithMaterial(), GetNDF(), NDF(), PndCAGBTracker::PickUpHits(), Reset(), SetNDF(), and SetTrackParam().

float_v PndCATrackParamVector::fP[5]

Definition at line 302 of file PndCATrackParamVector.h.

Referenced by AcceptWithMaterial(), CalculateFitParameters(), CorrectForMeanMaterial(), DzDs(), FilterVtx(), FilterWithMaterial(), GetDzDs(), GetKappa(), GetQPt(), GetSinPhi(), GetY(), GetZ(), Par(), PndCAGBTracker::PickUpHits(), PndCATrackParamVector(), QPt(), Reset(), SetDzDs(), SetPar(), SetQPt(), SetSinPhi(), SetTrackParam(), SetY(), SetZ(), SinPhi(), Transport0ToX(), TransportToX(), Y(), and Z().

float_v PndCATrackParamVector::fSignCosPhi

Definition at line 301 of file PndCATrackParamVector.h.

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

float_v PndCATrackParamVector::fX

Definition at line 300 of file PndCATrackParamVector.h.

Referenced by GetX(), Reset(), SetTrackParam(), SetX(), Transport(), Transport0(), Transport0ToX(), TransportToX(), and X().

---

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

• PndCATrackParamVector.h
• PndCATrackParamVector.cxx