PndFTSCATrackParamVector Class Reference

#include <PndFTSCATrackParamVector.h>

Classes
struct	PndFTSCATrackFitParam

Public Member Functions
	PndFTSCATrackParamVector ()
void	ConvertTrackParamToVector (PndFTSCATrackParam t0[float_v::Size], int nTracksV)
void	InitCovMatrix (float_v d2QMom=0.f)
void	InitByTarget (const FTSCATarget &target)
void	InitByHit (const FTSCAHitV &hit, const PndFTSCAParam &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
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 PndFTSCATrackParamVector &param, const float_m &m=float_m(true))
void	SetTrackParamOne (int iV, const PndFTSCATrackParamVector &param, int iVa)
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	SetErr2Y (float_v v)
void	SetErr2Z (float_v v)
void	SetErr2QPt (float_v v)
float_v	GetDist2 (const PndFTSCATrackParamVector &t) const
float_v	GetDistXZ2 (const PndFTSCATrackParamVector &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	TransportToX0WithMaterial (const float_v &x, const float_v &XOverX0, const float_v &XThimesRho, const float_v &Bz, const float maxSinPhi=.999f)
float_m	TransportToX0 (const float_v &x, const float_v &Bz, const float maxSinPhi=.999f, const float_m &mask=float_m(true))
float_m	TransportToX0 (const float_v &x, PndFTSCATrackLinearisationVector &t0, const float_v &Bz, const float maxSinPhi=.999f, float_v *DL=0, const float_m &mask=float_m(true))
float_m	TransportToX0 (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	TransportToX0WithMaterial (const float_v &x, PndFTSCATrackLinearisationVector &t0, PndFTSCATrackFitParam &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	TransportToX0WithMaterial (const float_v &x, PndFTSCATrackFitParam &par, const float_v &XOverX0, const float_v &XThimesRho, const float_v &Bz, const float maxSinPhi=.999f)
float_m	Rotate (const float_v &alpha, PndFTSCATrackLinearisationVector &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 FTSCAStripInfo &info, float_v err2, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
float_m	FilterWithMaterial (const float_v &y, const float_v &z, const float_v &r, const FTSCAStripInfo &info, float_v err2, float maxSinPhi=0.999f, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
void	CalculateFitParameters (PndFTSCATrackFitParam &par, const float_v &mass=0.13957f)
float_m	CorrectForMeanMaterial (const float_v &xOverX0, const float_v &xTimesRho, const PndFTSCATrackFitParam &par, const float_m &_mask)
float_m	FilterVtx (const float_v &xV, const float_v &yV, const CAX1X2MeasurementInfo &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 PndFTSCAParam &param, const float_m &mask=float_m(true))
float_m	Transport (const FTSCAHitV &hit, const PndFTSCAParam &p, const float_m &mask=float_m(true))
float_m	Filter (const FTSCAHitV &hit, const PndFTSCAParam &param, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
float_m	Transport (const FTSCAHit &hit, const PndFTSCAParam &p, const float_m &mask=float_m(true))
float_m	Filter (const FTSCAHit &hit, const PndFTSCAParam &param, const float_m &mask=float_m(true), const float_v &chi2Cut=10e10f)
float_m	AddTarget (const FTSCATarget &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)

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

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

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

Detailed Description

PndFTSCATrackParamVector class describes the track parametrisation which is used by the PndFTSCATracker slice tracker.

Definition at line 1269 of file PndFTSCATrackParamVector.h.

Constructor & Destructor Documentation

PndFTSCATrackParamVector::PndFTSCATrackParamVector ( )

inline

Definition at line 1274 of file PndFTSCATrackParamVector.h.

References fP, and i.

      : fX( Vc::Zero ),
      fSignCosPhi( Vc::Zero ),
      fChi2( Vc::Zero ),
      fNDF( 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 PndFTSCATrackParamVector::AddTarget	(	const FTSCATarget &	target,
		const float_m &	mask = float_m( true )
	)

Definition at line 1504 of file PndFTSCATrackParamVector.cxx.

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

Referenced by PndFTSCAGBTracker::FitTrackCA().

{
  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(), CAX1X2MeasurementInfo( target.Err2X1(), 0.f, target.Err2X2() ), eY, eZ, J, active );
  fNDF (static_cast<int_m>(active)) += target.NDF();

  return active;
}

float_v PndFTSCATrackParamVector::Angle ( ) const

inline

Definition at line 1317 of file PndFTSCATrackParamVector.h.

References fAlpha.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), SetTrackParam(), and SetTrackParamOne().

{ return fAlpha; }

float_v PndFTSCATrackParamVector::ApproximateBetheBloch ( const float_v &  beta2 )

static

Definition at line 1019 of file PndFTSCATrackParamVector.cxx.

References 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 PndFTSCATrackParamVector::BetheBlochGas ( const float_v &  bg )

static

Definition at line 998 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::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 945 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::BetheBlochSolid ( const float_v &  bg )

static

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

Definition at line 1041 of file PndFTSCATrackParamVector.cxx.

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

Referenced by PndFTSCAGBTracker::FitTrack(), Transport(), and TransportToX0WithMaterial().

{
  //*!

  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 PndFTSCATrackParamVector::Chi2 ( ) const

inline

Definition at line 1331 of file PndFTSCATrackParamVector.h.

References fChi2.

Referenced by ConvertTrackParamToVector(), operator<<(), and PndFTSCAGBTracker::Refit_1().

{ return fChi2; }

void PndFTSCATrackParamVector::ConvertTrackParamToVector	(	PndFTSCATrackParam	t0[float_v::Size],
		int	nTracksV
	)

Definition at line 557 of file PndFTSCATrackParamVector.cxx.

References Angle(), Chi2(), Cov(), NDF(), Par(), SetAngle(), SetChi2(), SetCov(), SetNDF(), SetPar(), SetSignCosPhi(), SetX(), SignCosPhi(), t0, and X().

Referenced by PndFTSCAGBTracker::FitTracks().

{
  float_v tmpVec;
  int_v tmpVecShort;
  float_v::Memory tmpFloat;
  int_v::Memory tmpShort;

  for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].X();
  tmpVec.load( tmpFloat );
  SetX(tmpVec);
  for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].SignCosPhi();
  tmpVec.load( tmpFloat );
  SetSignCosPhi(tmpVec);

  for(int iP=0; iP<5; iP++)
  {
    for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].Par()[iP];
    tmpVec.load( tmpFloat );
    SetPar(iP,tmpVec);
  }
  for(int iC=0; iC<15; iC++)
  {
    for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].Cov()[iC];
    tmpVec.load( tmpFloat );
    SetCov(iC,tmpVec);
  }
  for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].Chi2();
  tmpVec.load( tmpFloat );
  SetChi2(tmpVec);
  for(int iV=0; iV < nTracksV; iV++) tmpShort[iV] = t0[iV].NDF();
  tmpVecShort.load( tmpShort );
  SetNDF(tmpVecShort);

  for(int iV=0; iV < nTracksV; iV++) tmpFloat[iV] = t0[iV].Angle();
  tmpVec.load( tmpFloat );
  SetAngle(tmpVec);
}

float_m PndFTSCATrackParamVector::CorrectForMeanMaterial	(	const float_v &	xOverX0,
		const float_v &	xTimesRho,
		const PndFTSCATrackFitParam &	par,
		const float_m &	_mask
	)

Definition at line 1072 of file PndFTSCATrackParamVector.cxx.

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

Referenced by TransportToX0WithMaterial().

{
  //------------------------------------------------------------------
  // 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 >= float_v(Vc::Zero) || !_mask ).isFull() );
  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 >= float_v(Vc::Zero) || !mask ).isFull() );
  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& PndFTSCATrackParamVector::Cov ( int  i ) const

inline

Definition at line 1361 of file PndFTSCATrackParamVector.h.

References i.

Referenced by PndFTSCAGBTracker::FitTrack(), operator<<(), PndFTSCAGBTracker::PickUpHits(), SetTrackParam(), and SetTrackParamOne().

{ return fC[i]; }

const float_v* PndFTSCATrackParamVector::Cov ( ) const

inlineprivate

Definition at line 1366 of file PndFTSCATrackParamVector.h.

References fC.

Referenced by ConvertTrackParamToVector().

{ return fC; }

float_v* PndFTSCATrackParamVector::Cov ( )

inlineprivate

Definition at line 1368 of file PndFTSCATrackParamVector.h.

References fC.

{ return fC; }

float_v PndFTSCATrackParamVector::DzDs ( ) const

inline

Definition at line 1314 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by TransportJ0ToX0(), TransportToX0(), and Tx2().

{ return fP[3]; }

float_v PndFTSCATrackParamVector::Err2DzDs ( ) const

inline

Definition at line 1337 of file PndFTSCATrackParamVector.h.

{ return fC[9]; }

float_v PndFTSCATrackParamVector::Err2QMomentum ( ) const

inline

Definition at line 1358 of file PndFTSCATrackParamVector.h.

{ return fC[14]; }

float_v PndFTSCATrackParamVector::Err2QPt ( ) const

inline

Definition at line 1338 of file PndFTSCATrackParamVector.h.

{ return fC[14]; }

float_v PndFTSCATrackParamVector::Err2SinPhi ( ) const

inline

Definition at line 1336 of file PndFTSCATrackParamVector.h.

{ return fC[5]; }

float_v PndFTSCATrackParamVector::Err2X1 ( ) const

inline

Definition at line 1354 of file PndFTSCATrackParamVector.h.

Referenced by PndFTSCAGBTracker::PickUpHits().

{ return fC[0]; }

float_v PndFTSCATrackParamVector::Err2X2 ( ) const

inline

Definition at line 1355 of file PndFTSCATrackParamVector.h.

{ return fC[2]; }

float_v PndFTSCATrackParamVector::Err2Y ( ) const

inline

Definition at line 1334 of file PndFTSCATrackParamVector.h.

{ return fC[0]; }

float_v PndFTSCATrackParamVector::Err2Z ( ) const

inline

Definition at line 1335 of file PndFTSCATrackParamVector.h.

{ return fC[2]; }

float_m PndFTSCATrackParamVector::Filter	(	const FTSCAHitV &	hit,
		const PndFTSCAParam &	param,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 1464 of file PndFTSCATrackParamVector.cxx.

References FTSCAHitV::Err2X1(), FTSCAHitV::Err2X2(), FTSCAHitV::ErrX12(), FTSCAStation::f, FilterWithMaterial(), FTSCAHitV::IStations(), FTSCAHitV::IsValid(), FTSCAStation::NDF, NDF(), PndFTSCAParam::Station(), FTSCAHitV::X1(), and FTSCAHitV::X2().

Referenced by CAFunctionality::FitIteration(), PndFTSCAGBTracker::FitTrackCA(), PndFTSCAGBTracker::PickUpHits(), and PndFTSCAGBTracker::Refit().

{
  const int iSta = hit.IStations()[hit.IsValid().firstOne()];
  if ( param.Station( iSta ).NDF == 1 ) {
#ifdef DRIFT_TUBES
    return FilterWithMaterial( hit.XWire1(), hit.XWire2(), hit.RSigned(), param.Station( iSta ).f, hit.Err2R(), 0.999f, mask, chi2Cut );
#else
    assert(0);
#endif
  }
  return FilterWithMaterial( hit.X1(), hit.X2(), hit.Err2X1(), hit.ErrX12(), hit.Err2X2(), 0.999f, mask, int_v(param.Station( hit.IStations()[0] ).NDF), chi2Cut );
}

float_m PndFTSCATrackParamVector::Filter	(	const FTSCAHit &	hit,
		const PndFTSCAParam &	param,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 1492 of file PndFTSCATrackParamVector.cxx.

References FTSCAHit::Err2X1(), FTSCAHit::Err2X2(), FTSCAHit::ErrX12(), FilterWithMaterial(), FTSCAHit::IStation(), NDF(), PndFTSCAParam::Station(), FTSCAHit::X1(), and FTSCAHit::X2().

{
  if ( param.Station( hit.IStation() ).NDF == 1 ) {
#ifdef DRIFT_TUBES
    return FilterWithMaterial( hit.XWire1(), hit.XWire2(), hit.RSigned(), param.Station( hit.IStation() ).f, hit.Err2R(), 0.999f, mask, chi2Cut );
#else
    assert(0);
#endif
  }
  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 PndFTSCATrackParamVector::FilterVtx ( const float_v &  xV, const float_v &  yV, const CAX1X2MeasurementInfo &  info, float_v &  extrDx, float_v &  extrDy, float_v  J[], const float_m &  active = float_m( true )  )

inline

Definition at line 1917 of file PndFTSCATrackParamVector.h.

References CAX1X2MeasurementInfo::C00(), CAX1X2MeasurementInfo::C10(), c10, CAX1X2MeasurementInfo::C11(), c11, c20, c21, c22, 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 PndFTSCATrackParamVector::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 1168 of file PndFTSCATrackParamVector.cxx.

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

Referenced by Filter(), and PndFTSCAGBTracker::FitTrack().

{
  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 PndFTSCATrackParamVector::FilterWithMaterial	(	const float_v &	y,
		const float_v &	z,
		const FTSCAStripInfo &	info,
		float_v	err2,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 1263 of file PndFTSCATrackParamVector.cxx.

References CAMath::Abs(), c10, c11, c20, c21, FTSCAStripInfo::cos, f, fC, fChi2, fNDF, fP, FTSCAStripInfo::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_m PndFTSCATrackParamVector::FilterWithMaterial	(	const float_v &	y,
		const float_v &	z,
		const float_v &	r,
		const FTSCAStripInfo &	info,
		float_v	err2,
		float	maxSinPhi = 0.999f,
		const float_m &	mask = float_m( true ),
		const float_v &	chi2Cut = 10e10f
	)

Definition at line 1347 of file PndFTSCATrackParamVector.cxx.

References CAMath::Abs(), c10, c11, c20, c21, FTSCAStripInfo::cos, f, fC, fChi2, fNDF, fP, h1, r, rsqrt(), FTSCAStripInfo::sin, SinPhi(), sqrt(), mrfdata_8b_error::success, and y0.

{
    // linearize track in current point, which must be x == x0
    // distance between wire and track are : r = - ({x,y,z} - {x0,y0,z0}) * e_o / |e_o|, where ort e_o = |[e_t x e_s]|
  const float_v& etx = sqrt( 1 - SinPhi()*SinPhi() );
  const float_v& ety = SinPhi();
  const float_v& etz = 1.f;
  const float_v& eox = - ety * info.cos - etz * info.sin; //  e_o[x] = - e_t[y] * e_s[z] + e_t[z] * e_s[y]
  const float_v& eoy = etx * info.cos; // e_o[y] = - 0 + e_t[x] * e_s[z]
  const float_v& eoz = etx * info.sin; // e_o[z] = - e_t[x] * e_s[y] + 0

  const float_v& iEo = rsqrt( eox*eox + eoy*eoy + eoz*eoz );
  const float_v& h0 = eoy*iEo;
  const float_v& h1 = eoz*iEo;

  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 zeta = h0*(fP[0] - y0) + h1*(fP[1] - z0) - r;

    // F = CH'
  const float_v  F0 = h0*c00 + h1*c10;
  const float_v  F1 = h0*c10 + h1*c11;

  const float_v  HCH = ( F0*h0 + F1*h1 );

  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 = h0*c20 + h1*c21;
  const float_v  F3 = h0*c30 + h1*c31;
  const float_v  F4 = h0*c40 + h1*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 PndFTSCATrackParamVector::GetChi2 ( ) const

inline

Definition at line 1347 of file PndFTSCATrackParamVector.h.

References fChi2.

Referenced by SetTrackParam(), and SetTrackParamOne().

{ return fChi2; }

float_v PndFTSCATrackParamVector::GetCosPhi ( ) const

inline

Definition at line 1352 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::GetCosPhiPositive ( ) const

inline

Definition at line 1351 of file PndFTSCATrackParamVector.h.

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

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

void PndFTSCATrackParamVector::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 737 of file PndFTSCATrackParamVector.cxx.

References a, CAMath::Abs(), dx, dy, 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 PndFTSCATrackParamVector::GetDist2

(

const PndFTSCATrackParamVector &

t

)

const

Definition at line 700 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::GetDistXZ2

(

const PndFTSCATrackParamVector &

t

)

const

Definition at line 710 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::GetDzDs ( ) const

inline

Definition at line 1344 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint().

{ return fP[3]; }

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

inline

Definition at line 1350 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint(), and GetS().

{ return fP[4]*Bz; }

int_v PndFTSCATrackParamVector::GetNDF ( ) const

inline

Definition at line 1348 of file PndFTSCATrackParamVector.h.

References fNDF.

Referenced by SetTrackParam(), and SetTrackParamOne().

{ return fNDF; }

float_v PndFTSCATrackParamVector::GetQPt ( ) const

inline

Definition at line 1345 of file PndFTSCATrackParamVector.h.

References fP.

{ return fP[4]; }

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

Definition at line 720 of file PndFTSCATrackParamVector.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 PndFTSCATrackParamVector::GetSignCosPhi ( ) const

inline

Definition at line 1346 of file PndFTSCATrackParamVector.h.

References fSignCosPhi.

{ return fSignCosPhi; }

float_v PndFTSCATrackParamVector::GetSinPhi ( ) const

inline

Definition at line 1343 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by GetDCAPoint(), and GetS().

{ return fP[2]; }

float_v PndFTSCATrackParamVector::GetX ( ) const

inline

Definition at line 1340 of file PndFTSCATrackParamVector.h.

References fX.

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

{ return fX; }

float_v PndFTSCATrackParamVector::GetY ( ) const

inline

Definition at line 1341 of file PndFTSCATrackParamVector.h.

References fP.

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

{ return fP[0]; }

float_v PndFTSCATrackParamVector::GetZ ( ) const

inline

Definition at line 1342 of file PndFTSCATrackParamVector.h.

References fP.

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

{ return fP[1]; }

void PndFTSCATrackParamVector::InitByHit	(	const FTSCAHitV &	hit,
		const PndFTSCAParam &	param,
		const float_v &	dQMom
	)

Definition at line 652 of file PndFTSCATrackParamVector.cxx.

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

Referenced by PndFTSCAGBTracker::Create1Plets(), and PndFTSCAGBTracker::FitTrackCA().

{
  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 FTSCAStation sta0 = param.Station(ista0);

  // CAFieldValue B;
  // sta0.fieldSlice.GetFieldValue( hit.X1(), hit.X2(), B );
  // SetField( 1, B, sta0.z );

  // if ( ista0 > 0 ) {
  //   const FTSCAStation 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 PndFTSCATrackParamVector::InitByTarget

(

const FTSCATarget &

target

)

Definition at line 619 of file PndFTSCATrackParamVector.cxx.

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

Referenced by PndFTSCAGBTracker::Create1Plets(), CAFunctionality::FitIteration(), PndFTSCAGBTracker::FitTrackCA(), and PndFTSCAGBTracker::Refit().

{
  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 PndFTSCATrackParamVector::InitCovMatrix

(

float_v

d2QMom = 0.f

)

Definition at line 596 of file PndFTSCATrackParamVector.cxx.

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

{
  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 PndFTSCATrackParamVector::InitDirection ( float_v  r0, float_v  r1, float_v  r2  )

inline

Definition at line 1290 of file PndFTSCATrackParamVector.h.

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

Referenced by PndFTSCAGBTracker::Create1Plets(), CAFunctionality::FitIteration(), PndFTSCAGBTracker::FitTrackCA(), and PndFTSCAGBTracker::Refit().

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

int_v PndFTSCATrackParamVector::NDF ( ) const

inline

Definition at line 1332 of file PndFTSCATrackParamVector.h.

References fNDF.

Referenced by ConvertTrackParamToVector(), Filter(), operator<<(), and PndFTSCAGBTracker::Refit_1().

{ return fNDF; }

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

inline

Definition at line 1360 of file PndFTSCATrackParamVector.h.

References fP, and i.

Referenced by PndFTSCAGBTracker::FitTrack(), operator<<(), SetTrackParam(), and SetTrackParamOne().

{ return fP[i]; }

const float_v* PndFTSCATrackParamVector::Par ( ) const

inlineprivate

Definition at line 1365 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by ConvertTrackParamToVector().

{ return fP; }

float_v* PndFTSCATrackParamVector::Par ( )

inlineprivate

Definition at line 1367 of file PndFTSCATrackParamVector.h.

References fP.

{ return fP; }

float_v PndFTSCATrackParamVector::QMomentum ( ) const

inline

Definition at line 1324 of file PndFTSCATrackParamVector.h.

References QPt().

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

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

float_v PndFTSCATrackParamVector::QPt ( ) const

inline

Definition at line 1315 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), QMomentum(), and TransportToX0().

{ return fP[4]; }

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

Definition at line 1116 of file PndFTSCATrackParamVector.cxx.

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

Referenced by PndFTSCAGBTracker::FitTrack(), and Transport().

{
  //* 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 PndFTSCATrackParamVector::Rotate ( const float_v &  alpha, const float  maxSinPhi = .999f, const float_m &  mask = float_m( true )  )

inline

Definition at line 1855 of file PndFTSCATrackParamVector.h.

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

{
  //* Rotate the coordinate system in XY on the angle alpha
  if ( (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 x = X(), y = Y(), sP = SinPhi(), cP = GetCosPhi();
  const float_v cosPhi = cP * cA + sP * sA;
  const float_v sinPhi = -cP * sA + sP * cA;

  float_m mReturn = mask && (CAMath::Abs( sinPhi ) < maxSinPhi) && (CAMath::Abs( cosPhi ) > 1.e-2f) && (CAMath::Abs( cP ) > 1.e-2f);

  mReturn &= abs(alpha) < 3.1415f * 0.25f; // allow turn by 45 degree only

  const float_v j0 = cP / cosPhi;
  const float_v j2 = cosPhi / cP;

  SetX( x*cA +  y*sA, mReturn);
  SetY( -x*sA +  y*cA, mReturn );
  SetSignCosPhi( CAMath::Abs(cosPhi)/cosPhi, mReturn );
  SetSinPhi( sinPhi, mReturn );


  //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](mReturn) *= j0 * j0;
  fC[1](mReturn) *= j0;
  fC[3](mReturn) *= j0;
  fC[6](mReturn) *= j0;
  fC[10](mReturn) *= j0;

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

  fAlpha(mReturn) += alpha;
  //cout<<"      "<<fC[0]<<" "<<fC[1]<<" "<<fC[6]<<" "<<fC[10]<<" "<<fC[4]<<" "<<fC[5]<<" "<<fC[8]<<" "<<fC[12]<<endl;
  return mReturn;
}

void PndFTSCATrackParamVector::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 1903 of file PndFTSCATrackParamVector.h.

References CAMath::Cos(), 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 PndFTSCATrackParamVector::SetAngle ( const float_v &  v )

inline

Definition at line 1419 of file PndFTSCATrackParamVector.h.

References fAlpha, and v.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::Create1Plets(), CAFunctionality::FitIteration(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), InitByHit(), and PndFTSCAGBTracker::Refit().

{  fAlpha = v; }

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

inline

Definition at line 1420 of file PndFTSCATrackParamVector.h.

References fAlpha, and v.

{  fAlpha(m) = v; }

void PndFTSCATrackParamVector::SetChi2 ( const float_v &  v )

inline

Definition at line 1413 of file PndFTSCATrackParamVector.h.

References fChi2, and v.

Referenced by ConvertTrackParamToVector(), FTSCATrack::Fit2Times(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTracks(), InitByHit(), InitByTarget(), InitCovMatrix(), PndFTSCAGBTracker::Refit(), and PndFTSCAGBTracker::Refit_1().

{  fChi2 = v; }

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

inline

Definition at line 1414 of file PndFTSCATrackParamVector.h.

References fChi2, and v.

{  fChi2(m) = v; }

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

inline

Definition at line 1395 of file PndFTSCATrackParamVector.h.

References i, and v.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTracks(), InitByHit(), InitByTarget(), and InitCovMatrix().

{ fC[i] = v; }

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

inline

Definition at line 1396 of file PndFTSCATrackParamVector.h.

References i, and m.

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

void PndFTSCATrackParamVector::SetDzDs ( const float_v &  v )

inline

Definition at line 1406 of file PndFTSCATrackParamVector.h.

References fP, and v.

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

{  fP[3] = v; }

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

inline

Definition at line 1407 of file PndFTSCATrackParamVector.h.

References fP, and m.

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

void PndFTSCATrackParamVector::SetErr2QPt ( float_v  v )

inline

Definition at line 1424 of file PndFTSCATrackParamVector.h.

References v.

{ fC[14] = v; }

void PndFTSCATrackParamVector::SetErr2Y ( float_v  v )

inline

Definition at line 1422 of file PndFTSCATrackParamVector.h.

References v.

{ fC[0] = v; }

void PndFTSCATrackParamVector::SetErr2Z ( float_v  v )

inline

Definition at line 1423 of file PndFTSCATrackParamVector.h.

References v.

{ fC[2] = v; }

void PndFTSCATrackParamVector::SetNDF ( int  v )

inline

Definition at line 1415 of file PndFTSCATrackParamVector.h.

References fNDF, and v.

Referenced by ConvertTrackParamToVector(), FTSCATrack::Fit2Times(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTracks(), InitByHit(), InitByTarget(), InitCovMatrix(), and PndFTSCAGBTracker::Refit().

{ fNDF = v; }

void PndFTSCATrackParamVector::SetNDF ( const int_v &  v )

inline

Definition at line 1416 of file PndFTSCATrackParamVector.h.

References fNDF, and v.

{ fNDF = v; }

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

inline

Definition at line 1417 of file PndFTSCATrackParamVector.h.

References fNDF, and v.

{ fNDF(m) = v; }

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

inline

Definition at line 1393 of file PndFTSCATrackParamVector.h.

References fP, i, and v.

Referenced by ConvertTrackParamToVector().

{ fP[i] = v; }

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

inline

Definition at line 1394 of file PndFTSCATrackParamVector.h.

References fP, i, and m.

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

void PndFTSCATrackParamVector::SetQMomentum ( const float_v &  v )

inline

Definition at line 1409 of file PndFTSCATrackParamVector.h.

References SetQPt().

Referenced by CAFunctionality::FitIteration().

{  SetQPt(v); }

void PndFTSCATrackParamVector::SetQPt ( const float_v &  v )

inline

Definition at line 1408 of file PndFTSCATrackParamVector.h.

References fP, and v.

Referenced by PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), InitByHit(), InitByTarget(), and SetQMomentum().

{  fP[4] = v; }

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

inline

Definition at line 1410 of file PndFTSCATrackParamVector.h.

References fP, and m.

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

void PndFTSCATrackParamVector::SetSignCosPhi ( const float_v &  v )

inline

Definition at line 1411 of file PndFTSCATrackParamVector.h.

References fSignCosPhi, and v.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::FitTrack(), InitByHit(), InitByTarget(), InitDirection(), and Rotate().

{  fSignCosPhi = v; }

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

inline

Definition at line 1412 of file PndFTSCATrackParamVector.h.

References fSignCosPhi, and v.

{  fSignCosPhi(m) = v; }

void PndFTSCATrackParamVector::SetSinPhi ( const float_v &  v )

inline

Definition at line 1404 of file PndFTSCATrackParamVector.h.

References fP, and v.

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

{  fP[2] = v; }

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

inline

Definition at line 1405 of file PndFTSCATrackParamVector.h.

References fP, and m.

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

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

inline

Definition at line 1371 of file PndFTSCATrackParamVector.h.

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

Referenced by PndFTSCAGBTracker::FitTrackCA().

  {
    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();
    fAlpha(static_cast<int_m>(m)) = param.Angle();
  }

void PndFTSCATrackParamVector::SetTrackParamOne ( int  iV, const PndFTSCATrackParamVector &  param, int  iVa  )

inline

Definition at line 1382 of file PndFTSCATrackParamVector.h.

References Angle(), Cov(), fAlpha, fChi2, fNDF, fP, fSignCosPhi, fX, GetChi2(), GetNDF(), i, Par(), SignCosPhi(), and X().

Referenced by FTSCANPletV::CopyOne().

  {
    for(int i=0; i<5; i++) fP[i][iV] = param.Par()[i][iVa];
    for(int i=0; i<15; i++) fC[i][iV] = param.Cov()[i][iVa];
    fX[iV] = param.X()[iVa];
    fSignCosPhi[iV] = param.SignCosPhi()[iVa];
    fChi2[iV] = param.GetChi2()[iVa];
    fNDF[iV] = param.GetNDF()[iVa];
    fAlpha[iV] = param.Angle()[iVa];
  }

void PndFTSCATrackParamVector::SetX ( const float_v &  v )

inline

Definition at line 1398 of file PndFTSCATrackParamVector.h.

References fX, and v.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::FitTrackCA(), InitByHit(), InitByTarget(), PndFTSCAGBTracker::Refit_1(), and Rotate().

{  fX = v;    }

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

inline

Definition at line 1401 of file PndFTSCATrackParamVector.h.

References fX, and v.

{  fX( m ) = v;    }

void PndFTSCATrackParamVector::SetY ( const float_v &  v )

inline

Definition at line 1399 of file PndFTSCATrackParamVector.h.

References fP, and v.

Referenced by PndFTSCAGBTracker::FitTrackCA(), InitByHit(), InitByTarget(), PndFTSCAGBTracker::Refit_1(), and Rotate().

{  fP[0] = v; }

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

inline

Definition at line 1402 of file PndFTSCATrackParamVector.h.

References fP, and m.

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

void PndFTSCATrackParamVector::SetZ ( const float_v &  v )

inline

Definition at line 1400 of file PndFTSCATrackParamVector.h.

References fP, and v.

Referenced by PndFTSCAGBTracker::FitTrackCA(), InitByHit(), and InitByTarget().

{  fP[1] = v; }

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

inline

Definition at line 1403 of file PndFTSCATrackParamVector.h.

References fP, and m.

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

float_v PndFTSCATrackParamVector::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 1330 of file PndFTSCATrackParamVector.h.

References fSignCosPhi.

Referenced by ConvertTrackParamToVector(), operator<<(), PndFTSCATrackLinearisationVector::PndFTSCATrackLinearisationVector(), SetTrackParam(), SetTrackParamOne(), Tx1(), and Tx2().

{ return fSignCosPhi; }

float_v PndFTSCATrackParamVector::SinPhi ( ) const

inline

Definition at line 1313 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by FilterWithMaterial(), PndFTSCAGBTracker::FitTrack(), GetCosPhi(), GetCosPhiPositive(), Rotate(), RotateXY(), TransportJ0ToX0(), TransportToX0(), Tx1(), and Tx2().

{ return fP[2]; }

float_m PndFTSCATrackParamVector::Transport	(	const int_v &	ista,
		const PndFTSCAParam &	param,
		const float_m &	mask = float_m( true )
	)

Definition at line 1440 of file PndFTSCATrackParamVector.cxx.

References CalculateFitParameters(), PndFTSCAParam::cBz(), PndFTSCAParam::GetX0(), PndFTSCAParam::GetXOverX0(), PndFTSCAParam::GetXTimesRho(), and TransportToX0WithMaterial().

Referenced by PndFTSCAGBTracker::Create1Plets(), CAFunctionality::FitIteration(), PndFTSCAGBTracker::FitTrackCA(), PndFTSCAGBTracker::PickUpHits(), and PndFTSCAGBTracker::Refit().

{
  float_m active = mask;

  PndFTSCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  PndFTSCATrackLinearisationVector tE( *this );
  active &= TransportToX0WithMaterial( param.GetX0( 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 PndFTSCATrackParamVector::Transport	(	const FTSCAHitV &	hit,
		const PndFTSCAParam &	p,
		const float_m &	mask = float_m( true )
	)

Definition at line 1451 of file PndFTSCATrackParamVector.cxx.

References CAMath::Abs(), FTSCAHitV::Angle(), CalculateFitParameters(), PndFTSCAParam::cBz(), fAlpha, fX, PndFTSCAParam::GetXOverX0(), PndFTSCAParam::GetXTimesRho(), FTSCAHitV::IStations(), FTSCAHitV::IsValid(), Rotate(), TransportToX0WithMaterial(), and FTSCAHitV::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();
  PndFTSCATrackLinearisationVector tR( *this );
  active &= Rotate(  -fAlpha + hit.Angle(), tR, .999f, active );

  PndFTSCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  active &= TransportToX0WithMaterial( hit.X0(), tR, fitPar, param.GetXOverX0(hit.IStations(),active), param.GetXTimesRho(hit.IStations(),active), param.cBz(), 0.999f, active );
  return active;
}

float_m PndFTSCATrackParamVector::Transport	(	const FTSCAHit &	hit,
		const PndFTSCAParam &	p,
		const float_m &	mask = float_m( true )
	)

Definition at line 1477 of file PndFTSCATrackParamVector.cxx.

References CAMath::Abs(), FTSCAHit::Angle(), CalculateFitParameters(), PndFTSCAParam::cBz(), fAlpha, fX, PndFTSCAParam::GetXOverX0(), PndFTSCAParam::GetXTimesRho(), FTSCAHit::IStation(), Rotate(), TransportToX0WithMaterial(), and FTSCAHit::X0().

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

  float_m active = mask;

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

  PndFTSCATrackFitParam fitPar;
  CalculateFitParameters( fitPar );
  active &= TransportToX0WithMaterial( hit.X0(), tR, fitPar, param.GetXOverX0(hit.IStation()), param.GetXTimesRho(hit.IStation()), param.cBz(), 0.999f, active );
  return active;
}

float_m PndFTSCATrackParamVector::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 1994 of file PndFTSCATrackParamVector.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(active) = dS*ey;
  extrDz(active) = dS*DzDs();
  J[0](active) = dS;
  J[1](active) = 0;
  J[2](active) = h4;
  J[3](active) = dS;
  J[4](active) = dS;
  J[5](active) = 0;
  return active;
}

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

Definition at line 880 of file PndFTSCATrackParamVector.cxx.

References t0.

Referenced by TransportToX0WithMaterial().

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

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

  PndFTSCATrackLinearisationVector t0( *this );

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

float_m PndFTSCATrackParamVector::TransportToX0	(	const float_v &	x,
		PndFTSCATrackLinearisationVector &	t0,
		const float_v &	Bz,
		const float	maxSinPhi = .999f,
		float_v *	DL = 0,
		const float_m &	mask = float_m( true )
	)

Definition at line 767 of file PndFTSCATrackParamVector.cxx.

References CAMath::Abs(), CAMath::ASin(), c10, c11, c20, c21, c22, PndFTSCATrackLinearisationVector::CosPhi(), d, debugKF(), dx, dy, PndFTSCATrackLinearisationVector::DzDs(), f, fC, fP, fX, h2, h4, CAMath::Max(), CAMath::Min(), One, PndFTSCATrackLinearisationVector::QPt(), PndFTSCATrackLinearisationVector::SetCosPhi(), PndFTSCATrackLinearisationVector::SetSinPhi(), PndFTSCATrackLinearisationVector::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 PndFTSCATrackParamVector::TransportToX0 ( 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 1518 of file PndFTSCATrackParamVector.h.

References CAMath::Abs(), c20, c22, debugKF(), dx, DzDs(), 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 TransportToX0(" << 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;
//#define LOSE_DEBUG
#ifdef LOSE_DEBUG
  std::cout << " TrTo-sinPhi0 = " << sinPhi0 << std::endl;
#endif
//  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();
#ifdef LOSE_DEBUG
  std::cout << " TrTo-sinPhi = " << sinPhi << std::endl;
#endif
  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 PndFTSCATrackParamVector::TransportToX0WithMaterial	(	const float_v &	x,
		const float_v &	XOverX0,
		const float_v &	XThimesRho,
		const float_v &	Bz,
		const float	maxSinPhi = .999f
	)

Definition at line 929 of file PndFTSCATrackParamVector.cxx.

References CalculateFitParameters(), and par.

Referenced by PndFTSCAGBTracker::FitTrack(), Transport(), and TransportToX0WithMaterial().

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

  PndFTSCATrackFitParam par;
  CalculateFitParameters( par );
  return TransportToX0WithMaterial( x, par, XOverX0, XThimesRho, Bz, maxSinPhi );
}

float_m PndFTSCATrackParamVector::TransportToX0WithMaterial	(	const float_v &	x,
		PndFTSCATrackLinearisationVector &	t0,
		PndFTSCATrackFitParam &	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 894 of file PndFTSCATrackParamVector.cxx.

References CorrectForMeanMaterial(), PndFTSCATrackLinearisationVector::DzDs(), PndFTSCATrackLinearisationVector::SinPhi(), sqrt(), and TransportToX0().

{
  //* 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_ && TransportToX0( 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 );
  CorrectForMeanMaterial( XOverX0*koeff, XThimesRho*koeff, par, mask );

  return mask;
}

float_m PndFTSCATrackParamVector::TransportToX0WithMaterial	(	const float_v &	x,
		PndFTSCATrackFitParam &	par,
		const float_v &	XOverX0,
		const float_v &	XThimesRho,
		const float_v &	Bz,
		const float	maxSinPhi = .999f
	)

Definition at line 920 of file PndFTSCATrackParamVector.cxx.

References t0, and TransportToX0WithMaterial().

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

  PndFTSCATrackLinearisationVector t0( *this );
  return TransportToX0WithMaterial( x, t0, par, XOverX0, XThimesRho, Bz, maxSinPhi );
}

float_v PndFTSCATrackParamVector::Tx1 ( ) const

inline

Definition at line 1322 of file PndFTSCATrackParamVector.h.

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

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

float_v PndFTSCATrackParamVector::Tx2 ( ) const

inline

Definition at line 1323 of file PndFTSCATrackParamVector.h.

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

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

float_v PndFTSCATrackParamVector::X ( ) const

inline

Definition at line 1310 of file PndFTSCATrackParamVector.h.

References fX.

Referenced by ConvertTrackParamToVector(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), operator<<(), PndFTSCAGBTracker::PickUpHits(), PndFTSCAGBTracker::Refit_1(), Rotate(), RotateXY(), SetTrackParam(), SetTrackParamOne(), TransportJ0ToX0(), TransportToX0(), and X0().

{ return fX;    }

float_v PndFTSCATrackParamVector::X0 ( ) const

inline

Definition at line 1319 of file PndFTSCATrackParamVector.h.

References X().

Referenced by AddTarget().

{ return X(); }

float_v PndFTSCATrackParamVector::X1 ( ) const

inline

Definition at line 1320 of file PndFTSCATrackParamVector.h.

References Y().

Referenced by PndFTSCAGBTracker::PickUpHits().

{ return Y(); }

float_v PndFTSCATrackParamVector::X2 ( ) const

inline

Definition at line 1321 of file PndFTSCATrackParamVector.h.

References Z().

{ return Z(); }

float_v PndFTSCATrackParamVector::Y ( ) const

inline

Definition at line 1311 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by FilterVtx(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), PndFTSCAGBTracker::Refit_1(), Rotate(), RotateXY(), TransportToX0(), and X1().

{ return fP[0]; }

float_v PndFTSCATrackParamVector::Z ( ) const

inline

Definition at line 1312 of file PndFTSCATrackParamVector.h.

References fP.

Referenced by FilterVtx(), PndFTSCAGBTracker::FitTrack(), PndFTSCAGBTracker::FitTrackCA(), PndFTSCAGBTracker::PickUpHits(), PndFTSCAGBTracker::Refit_1(), TransportToX0(), and X2().

{ return fP[1]; }

Friends And Related Function Documentation

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

friend

Definition at line 1543 of file PndFTSCATrackParamVector.cxx.

{
  if ( out == std::cerr ) {
    out << "------------------------------ Track Param ------------------------------"
        << "\n             X: " << t.X()
        << "\n    SignCosPhi: " << t.SignCosPhi()
        << "\n          Chi2: " << t.Chi2()
        << "\n           NDF: " << t.NDF()
        << "\n             Y: " << t.Par()[0]
        << "\n             Z: " << t.Par()[1]
        << "\n        SinPhi: " << t.Par()[2]
        << "\n          DzDs: " << t.Par()[3]
        << "\n          q/Pt: " << t.Par()[4]
        << "\nCovariance Matrix\n";
    int i = 0;
    out << std::setprecision( 2 );
    for ( int step = 1; step <= 5; ++step ) {
      int end = i + step;
      for ( ; i < end; ++i ) {
        out << t.Cov()[i] << '\t';
      }
      out << "\n";
    }
    out << std::setprecision( 6 );
    return out << std::endl;
  }
  for ( int j = 0; j < uint_v::Size; ++j ) {
    out << t.X()[j] << " "
        << t.SignCosPhi()[j] << " "
        << t.Chi2()[j] << " "
        << t.NDF()[j]
        << std::endl;
    for ( int i = 0; i < 5; i++ ) out << t.Par()[i][j] << " ";
    out << std::endl;
    for ( int i = 0; i < 15; i++ ) out << t.Cov()[i][j] << " ";
    out << std::endl;
  }
  return out;
}

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

friend

Definition at line 1522 of file PndFTSCATrackParamVector.cxx.

{
  float_v::Memory x, s, p[5], c[15], chi2;
  int_v::Memory ndf;
  for ( int j = 0; j < uint_v::Size; ++j ) {
    in >> x[j];
    in >> s[j];
    for ( int i = 0; i < 5; i++ ) in >> p[i][j];
    for ( int i = 0; i < 15; i++ ) in >> c[i][j];
    in >> chi2[j];
    in >> ndf[j];
  }
  t.fX.load( x );
  t.fSignCosPhi.load( s );
  for ( int i = 0; i < 5; i++ ) t.fP[i].load( p[i] );
  for ( int i = 0; i < 5; i++ ) t.fC[i].load( c[i] );
  t.fChi2.load( chi2 );
  t.fNDF.load( ndf );
  return in;
}

friend class PndFTSCATrackParam

friend

Definition at line 1364 of file PndFTSCATrackParamVector.h.

Member Data Documentation

float_v PndFTSCATrackParamVector::fAlpha

private

Definition at line 1513 of file PndFTSCATrackParamVector.h.

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

float_v PndFTSCATrackParamVector::fC[15]

private

Definition at line 1509 of file PndFTSCATrackParamVector.h.

Referenced by CorrectForMeanMaterial(), Cov(), FilterWithMaterial(), operator>>(), Rotate(), and TransportToX0().

float_v PndFTSCATrackParamVector::fChi2

private

Definition at line 1510 of file PndFTSCATrackParamVector.h.

Referenced by Chi2(), FilterVtx(), FilterWithMaterial(), GetChi2(), operator>>(), SetChi2(), SetTrackParam(), and SetTrackParamOne().

int_v PndFTSCATrackParamVector::fNDF

private

Definition at line 1511 of file PndFTSCATrackParamVector.h.

Referenced by AddTarget(), FilterWithMaterial(), GetNDF(), NDF(), operator>>(), SetNDF(), SetTrackParam(), and SetTrackParamOne().

float_v PndFTSCATrackParamVector::fP[5]

private

Definition at line 1508 of file PndFTSCATrackParamVector.h.

Referenced by CalculateFitParameters(), CorrectForMeanMaterial(), DzDs(), FilterVtx(), FilterWithMaterial(), GetDzDs(), GetKappa(), GetQPt(), GetSinPhi(), GetY(), GetZ(), operator>>(), Par(), PndFTSCATrackParamVector(), QPt(), SetDzDs(), SetPar(), SetQPt(), SetSinPhi(), SetTrackParam(), SetTrackParamOne(), SetY(), SetZ(), SinPhi(), TransportToX0(), Y(), and Z().

float_v PndFTSCATrackParamVector::fSignCosPhi

private

Definition at line 1507 of file PndFTSCATrackParamVector.h.

Referenced by GetCosPhi(), GetSignCosPhi(), operator>>(), SetSignCosPhi(), SetTrackParam(), SetTrackParamOne(), and SignCosPhi().

float_v PndFTSCATrackParamVector::fX

private

Definition at line 1506 of file PndFTSCATrackParamVector.h.

Referenced by GetX(), operator>>(), SetTrackParam(), SetTrackParamOne(), SetX(), Transport(), TransportToX0(), and X().

---

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

• PndFTSCATrackParamVector.h
• PndFTSCATrackParamVector.cxx