KFPTrack Class Reference

#include <KFPTrack.h>

Public Member Functions
	KFPTrack ()
	~KFPTrack ()
int	GetID () const
bool	GetXYZPxPyPz (float *p) const
bool	GetCovarianceXYZPxPyPz (float cv[21]) const
bool	GetCovarianceXYZPxPyPz (double cv[21]) const
void	GetXYZ (float *position) const
void	GetPxPyPz (float *position) const
void	XvYvZv (float *position) const
void	PxPyPz (float *position) const
void	XvYvZv (double *position) const
void	PxPyPz (double *position) const
float	GetX () const
float	GetY () const
float	GetZ () const
float	GetPx () const
float	GetPy () const
float	GetPz () const
float	GetPt () const
float	GetP () const
void	GetCovarianceMatrix (float *covmatrix)
float	GetParameter (int i) const
float	GetCovariance (int i) const
int	Charge () const
float	GetChi2perNDF () const
float	GetChi2 () const
int	GetNDF () const
const float *	GetTrack () const
const float *	GetCovMatrix () const
void	SetParameters (float *position)
void	SetParameters (float x, float y, float z, float px, float py, float pz)
void	SetXYZ (float x, float y, float z)
void	SetPxPyPz (float px, float py, float pz)
void	SetID (int id)
void	SetX (float x)
void	SetY (float y)
void	SetZ (float z)
void	SetPx (float px)
void	SetPy (float py)
void	SetPz (float pz)
void	SetCharge (int q)
void	SetChi2 (float chi)
void	SetNDF (int ndf)
void	SetCovarianceMatrix (const float *C)
void	RotateXY (float alpha)
int	Id () const
void	SetId (int id)

Private Attributes
int	ID
float	fP [6]
float	fC [21]
float	fChi2
int	fQ
int	fNDF
int	fId

Detailed Description

Definition at line 6 of file KFPTrack.h.

Constructor & Destructor Documentation

KFPTrack::KFPTrack ( )

inline

Definition at line 10 of file KFPTrack.h.

{ }

KFPTrack::~KFPTrack ( )

inline

Definition at line 11 of file KFPTrack.h.

{ }

Member Function Documentation

int KFPTrack::Charge ( ) const

inline

Definition at line 62 of file KFPTrack.h.

References fQ.

Referenced by KFParticleSIMD::Create(), KFParticleFinder::FindParticles(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

{ return fQ; }

float KFPTrack::GetChi2 ( ) const

inline

Definition at line 64 of file KFPTrack.h.

References fChi2.

Referenced by KFParticleSIMD::Create(), KFParticleFinder::FindParticles(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

{ return fChi2;      }

float KFPTrack::GetChi2perNDF ( ) const

inline

Definition at line 63 of file KFPTrack.h.

References fChi2, and fNDF.

{ return fChi2/fNDF; }

float KFPTrack::GetCovariance ( int  i ) const

inline

Definition at line 60 of file KFPTrack.h.

References fC, and i.

{ return fC[i]; }

void KFPTrack::GetCovarianceMatrix ( float *  covmatrix )

inline

Definition at line 54 of file KFPTrack.h.

References fC, and i.

         {
           for (int i=0; i<21; i++)
             covmatrix[i] = fC[i];
         }

bool KFPTrack::GetCovarianceXYZPxPyPz ( float  cv[21] ) const

inline

Definition at line 21 of file KFPTrack.h.

References fC, and i.

Referenced by KFParticleSIMD::Create(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

         {
           for (int i=0; i<21; i++)
             cv[i] = fC[i];
           return 1;
         }

bool KFPTrack::GetCovarianceXYZPxPyPz ( double  cv[21] ) const

inline

Definition at line 28 of file KFPTrack.h.

References fC, and i.

         {
           for (int i=0; i<21; i++)
             cv[i] = fC[i];
           return 1;
         }

const float* KFPTrack::GetCovMatrix ( ) const

inline

Definition at line 68 of file KFPTrack.h.

References fC.

Referenced by KFParticleFinder::FindParticles().

{ return fC; }

int KFPTrack::GetID ( ) const

inline

Definition at line 13 of file KFPTrack.h.

References ID.

{ return ID; }

int KFPTrack::GetNDF ( ) const

inline

Definition at line 65 of file KFPTrack.h.

References fNDF.

Referenced by KFParticleSIMD::Create(), KFParticleFinder::FindParticles(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

{ return fNDF; }

float KFPTrack::GetP ( ) const

inline

Definition at line 52 of file KFPTrack.h.

References fP, and sqrt().

Referenced by KFParticleFinder::Find2DaughterDecay().

{ return sqrt(fP[3]*fP[3]+fP[4]*fP[4]+fP[5]*fP[5]); }

float KFPTrack::GetParameter ( int  i ) const

inline

Definition at line 59 of file KFPTrack.h.

References fP, and i.

{ return fP[i]; }

float KFPTrack::GetPt ( ) const

inline

Definition at line 51 of file KFPTrack.h.

References fP, and sqrt().

Referenced by KFParticleFinder::Find2DaughterDecay(), and KFParticleFinder::FindDMesLambdac().

{ return sqrt(fP[3]*fP[3]+fP[4]*fP[4]); }

float KFPTrack::GetPx ( ) const

inline

Definition at line 47 of file KFPTrack.h.

References fP.

Referenced by RotateXY().

{ return fP[3]; }

void KFPTrack::GetPxPyPz ( float *  position ) const

inline

Definition at line 37 of file KFPTrack.h.

References fP.

{position[0] = fP[3]; position[1] = fP[4]; position[2] = fP[5];}

float KFPTrack::GetPy ( ) const

inline

Definition at line 48 of file KFPTrack.h.

References fP.

Referenced by RotateXY().

{ return fP[4]; }

float KFPTrack::GetPz ( ) const

inline

Definition at line 49 of file KFPTrack.h.

References fP.

{ return fP[5]; }

const float* KFPTrack::GetTrack ( ) const

inline

Definition at line 67 of file KFPTrack.h.

References fP.

Referenced by KFParticleFinder::FindParticles().

{ return fP; }

float KFPTrack::GetX ( ) const

inline

Definition at line 44 of file KFPTrack.h.

References fP.

Referenced by RotateXY().

{ return fP[0]; }

void KFPTrack::GetXYZ ( float *  position ) const

inline

Definition at line 36 of file KFPTrack.h.

References fP.

{position[0] = fP[0]; position[1] = fP[1]; position[2] = fP[2];}

bool KFPTrack::GetXYZPxPyPz ( float *  p ) const

inline

Definition at line 15 of file KFPTrack.h.

References fP, and i.

         {
           for(int i=0; i<6; i++)
             p[i] = fP[i];
           return 1;
         }

float KFPTrack::GetY ( ) const

inline

Definition at line 45 of file KFPTrack.h.

References fP.

Referenced by RotateXY().

{ return fP[1]; }

float KFPTrack::GetZ ( ) const

inline

Definition at line 46 of file KFPTrack.h.

References fP.

Referenced by RotateXY().

{ return fP[2]; }

int KFPTrack::Id ( ) const

inline

Definition at line 108 of file KFPTrack.h.

References fId.

Referenced by KFParticleFinder::FindParticles().

{ return fId; };

void KFPTrack::PxPyPz ( float *  position ) const

inline

Definition at line 40 of file KFPTrack.h.

References fP.

Referenced by KFParticleSIMD::Create(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

{position[0] = fP[3]; position[1] = fP[4]; position[2] = fP[5];}

void KFPTrack::PxPyPz ( double *  position ) const

inline

Definition at line 42 of file KFPTrack.h.

References fP.

{position[0] = fP[3]; position[1] = fP[4]; position[2] = fP[5];}

void KFPTrack::RotateXY

(

float

alpha

)

Definition at line 7 of file KFPTrack.cxx.

References cos(), fC, GetPx(), GetPy(), GetX(), GetY(), GetZ(), i, SetPx(), SetPy(), SetX(), SetY(), SetZ(), sin(), x, and y.

Referenced by PndFTSTopoReconstructor::Init().

{
  const float cA = cos( alpha );
  const float sA = sin( alpha );

#if 0
    //float J[6][6] = { {  cA, sA, 0,  0,  0,  0 }, // X
    //                  { -sA, cA, 0,  0,  0,  0 }, // Y
    //                  {  0,  0,  1,  0,  0,  0 }, // Z
    //                  {  0,  0,  0,  cA, sA, 0 }, // Px
    //                  {  0,  0,  0, -sA, cA, 0 }, // Py
    //                  {  0,  0,  0,  0,  0,  1 } }; // Pz
  
#if defined(PANDA_STT) || defined(PANDA_FTS)
  float J[2][2] = { { sA,  cA },
                    { cA, -sA } };
#else
  float J[2][2] = { {  cA, sA },
                    { -sA, cA } };
#endif

  { // convert x, y TODO optimize
    const float x  = GetX(),  y = GetY();

#if defined(PANDA_STT) || defined(PANDA_FTS)
    SetX( x*sA +  y*cA );
    SetY( x*cA -  y*sA );
#else
    SetX(  x*cA +  y*sA );
    SetY( -x*sA +  y*cA );
#endif

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

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

  { // convert Px, Py TODO optimize
    const float x  = GetPx(),  y = GetPy();
#if defined(PANDA_STT) || defined(PANDA_FTS)
    SetPx( -x*sA -  y*cA );
    SetPy( -x*cA +  y*sA );
    float J[2][2] = { { -sA, -cA },
                      { -cA, sA } };
#else
    SetPx( x*cA  +  y*sA );
    SetPy( -x*sA +  y*cA );
#endif

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

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

      // symmetric matrix -> triangular

    fC[9] = Cov1[0][0];
    fC[13] = Cov1[0][1];
    fC[14] = Cov1[1][1];
  }
#if defined(PANDA_STT) || defined(PANDA_FTS)
  //Change sign of px, py, z
//   SetPx(-GetPx());
//   SetPy(-GetPy());
  SetZ( -GetZ() );
//   fC[ 3] = -fC[ 3];
//   fC[ 4] = -fC[ 4];
//   fC[ 8] = -fC[ 8];
//   fC[12] = -fC[12];

  fC[ 3] = -fC[ 3];
  fC[ 4] = -fC[ 4];
  fC[ 6] = -fC[ 6];
  fC[ 7] = -fC[ 7];
  fC[10] = -fC[10];
  fC[11] = -fC[11];
  fC[17] = -fC[17];
  fC[18] = -fC[18];
  fC[19] = -fC[19];
#endif

#else
#if defined(PANDA_STT) || defined(PANDA_FTS)
    //float J[6][6] = { {  sA,  cA,  0,  0,    0,  0 }, // X
    //                  {  cA, -sA,  0,  0,    0,  0 }, // Y
    //                  {  0,   0,  -1,  0,    0,  0 }, // Z
    //                  {  0,   0,   0, -sA, -cA,  0 }, // Px
    //                  {  0,   0,   0, -cA,  sA,  0 }, // Py
    //                  {  0,   0,   0,  0,    0,  1 } }; // Pz

    const float x  = GetX(),  y = GetY();

    SetX( x*sA +  y*cA );
    SetY( x*cA -  y*sA );

    const float px  = GetPx(),  py = GetPy();

    SetPx( -px*sA -  py*cA );
    SetPy( -px*cA +  py*sA );

    SetZ( -GetZ() );

    float cov[21];
    for(int iC=0; iC<21; iC++)
      cov[iC] = fC[iC];

    fC[0] = cA*cA* cov[2] + 2*cA*cov[1]*sA + cov[0]*sA*sA;
    fC[1] = cA*cA*cov[1] + cA *(cov[0] - cov[2])* sA - cov[1]* sA*sA;
    fC[2] = cA*cA*cov[0] - 2* cA *cov[1] *sA + cov[2]* sA*sA,

    fC[3] = -cA* cov[4] - cov[3]* sA;
    fC[4] = -cA* cov[3] + cov[4]* sA;
    fC[5] = cov[5];

    fC[6] = -cA*cA* cov[11] - cA* (cov[10] + cov[7])* sA - cov[6]* sA*sA;
    fC[7] = -cA*cA* cov[10] + cA* (cov[11] - cov[6])* sA + cov[7]* sA*sA;
    fC[8] =  cA* cov[12] + cov[8]* sA;
    fC[9] = cA*cA* cov[14] + 2* cA* cov[13]* sA + cov[9]* sA*sA;

    fC[10] = -cA*cA* cov[7] + cA* (cov[11] - cov[6])* sA + cov[10]* sA*sA;
    fC[11] = -cA*cA* cov[6] + cA* (cov[10] + cov[7])* sA - cov[11]* sA*sA;
    fC[12] =  cA* cov[8] - cov[12]* sA;
    fC[13] = cA*cA* cov[13] + cA* (-cov[14] + cov[9])* sA - cov[13]* sA*sA;
    fC[14] = cA*cA* cov[9] - 2 *cA* cov[13]* sA + cov[14]* sA*sA;

    fC[15] = cA *cov[16] + cov[15]* sA;
    fC[16] = cA *cov[15] - cov[16]* sA;
    fC[17] = -cov[17];
    fC[18] =  -cA* cov[19] - cov[18]* sA;
    fC[19] = -cA * cov[18] + cov[19]* sA;
    fC[20] = cov[20];
#else
    //float J[6][6] = { {  cA, -sA,  0,    0,    0,   0 }, // X
    //                  {  sA,  cA,  0,    0,    0,   0 }, // Y
    //                  {  0,   0,   1,    0,    0,   0 }, // Z
    //                  {  0,   0,   0,   cA,  -sA,   0 }, // Px
    //                  {  0,   0,   0,   sA,   cA,   0 }, // Py
    //                  {  0,   0,   0,    0,    0,   1 } }; // Pz

    const float x  = GetX(),  y = GetY();

    SetX( x*cA -  y*sA );
    SetY( x*sA +  y*cA );

    const float px  = GetPx(),  py = GetPy();

    SetPx( px*cA - py*sA );
    SetPy( px*sA + py*cA );

    float cov[21];
    for(int iC=0; iC<21; iC++)
      cov[iC] = fC[iC];

    fC[0] = cA*cA* cov[0] - 2*cA*cov[1]*sA + cov[2]*sA*sA;

    fC[1] = cA*cA*cov[1] + cA *(cov[0] - cov[2])* sA - cov[1]* sA*sA;
    fC[2] = cA*cA*cov[2] + 2* cA *cov[1] *sA + cov[0]* sA*sA,

    fC[3] = cA* cov[3] - cov[4]* sA;
    fC[4] = cA* cov[4] + cov[3]* sA;
    fC[5] = cov[5];

    fC[6] = cA* (cA* cov[6] - cov[10]* sA) - sA* (cA* cov[7] - cov[11]* sA);
    fC[7] = sA* (cA* cov[6] - cov[10]* sA) + cA* (cA* cov[7] - cov[11]* sA);
    fC[8] = cA* cov[8] - cov[12]* sA;
    fC[9] = cA* (cA* cov[9] - cov[13]* sA) - sA *(cA* cov[13] - cov[14]* sA);

    fC[10] = cA *(cov[10]* cA + cov[6]* sA) - (cov[11]* cA + cov[7] *sA) *sA;
    fC[11] = cA *(cov[11]* cA + cov[7]* sA) + (cov[10]* cA + cov[6] *sA) *sA;
    fC[12] =  cov[12] *cA + cov[8]* sA;
    fC[13] = -sA* (cA* cov[14] + cov[13]* sA) + cA* (cA* cov[13] + cov[9]* sA),
    fC[14] = cA *(cA* cov[14] + cov[13]* sA) + sA* (cA* cov[13] + cov[9]* sA);

    fC[15] = cA *cov[15] - cov[16]* sA;
    fC[16] = cA *cov[16] + cov[15]* sA;
    fC[17] = cov[17];
    fC[18] = cA* cov[18] - cov[19]* sA;
    fC[19] = cA * cov[19] + cov[18]* sA;
    fC[20] = cov[20];
#endif
#endif

} // RotateXY

void KFPTrack::SetCharge ( int  q )

inline

Definition at line 96 of file KFPTrack.h.

References fQ.

Referenced by PndFTSTopoReconstructor::Init().

{ fQ = q; }

void KFPTrack::SetChi2 ( float  chi )

inline

Definition at line 97 of file KFPTrack.h.

References fChi2.

Referenced by PndFTSTopoReconstructor::Init().

{ fChi2 = chi; }

void KFPTrack::SetCovarianceMatrix ( const float *  C )

inline

Definition at line 100 of file KFPTrack.h.

References fC, and i.

Referenced by PndFTSTopoReconstructor::Init().

  {
    for (int i=0; i<21; i++)
      fC[i] = C[i];
  }

void KFPTrack::SetID ( int  id )

inline

Definition at line 88 of file KFPTrack.h.

References ID.

{ID = id;}

void KFPTrack::SetId ( int  id )

inline

Definition at line 109 of file KFPTrack.h.

References fId.

Referenced by PndFTSTopoReconstructor::Init().

{ fId = id; };

void KFPTrack::SetNDF ( int  ndf )

inline

Definition at line 98 of file KFPTrack.h.

References fNDF.

Referenced by PndFTSTopoReconstructor::Init().

{ fNDF = ndf; }

void KFPTrack::SetParameters ( float *  position )

inline

Definition at line 70 of file KFPTrack.h.

References fP, and i.

       { 
         for(int i=0; i<6; i++)
           fP[i] = position[i];
       }

void KFPTrack::SetParameters ( float  x, float  y, float  z, float  px, float  py, float  pz  )

inline

Definition at line 75 of file KFPTrack.h.

References fP, pz, x, y, and z.

       { 
         fP[0] = x;  fP[1] = y;  fP[2] = z;
         fP[3] = px; fP[4] = py; fP[5] = pz;
       }

void KFPTrack::SetPx ( float  px )

inline

Definition at line 93 of file KFPTrack.h.

References fP.

Referenced by PndFTSTopoReconstructor::Init(), and RotateXY().

{ fP[3] = px; }

void KFPTrack::SetPxPyPz ( float  px, float  py, float  pz  )

inline

Definition at line 84 of file KFPTrack.h.

References fP, and pz.

       { 
         fP[3] = px; fP[4] = py; fP[5] = pz;
       }

void KFPTrack::SetPy ( float  py )

inline

Definition at line 94 of file KFPTrack.h.

References fP.

Referenced by PndFTSTopoReconstructor::Init(), and RotateXY().

{ fP[4] = py; }

void KFPTrack::SetPz ( float  pz )

inline

Definition at line 95 of file KFPTrack.h.

References fP, and pz.

Referenced by PndFTSTopoReconstructor::Init().

{ fP[5] = pz; }

void KFPTrack::SetX ( float  x )

inline

Definition at line 90 of file KFPTrack.h.

References fP, and x.

Referenced by PndFTSTopoReconstructor::Init(), and RotateXY().

{ fP[0] = x; }

void KFPTrack::SetXYZ ( float  x, float  y, float  z  )

inline

Definition at line 80 of file KFPTrack.h.

References fP, x, y, and z.

       { 
         fP[0] = x;  fP[1] = y;  fP[2] = z;
       }

void KFPTrack::SetY ( float  y )

inline

Definition at line 91 of file KFPTrack.h.

References fP, and y.

Referenced by PndFTSTopoReconstructor::Init(), and RotateXY().

{ fP[1] = y; }

void KFPTrack::SetZ ( float  z )

inline

Definition at line 92 of file KFPTrack.h.

References fP, and z.

Referenced by PndFTSTopoReconstructor::Init(), and RotateXY().

{ fP[2] = z; }

void KFPTrack::XvYvZv ( float *  position ) const

inline

Definition at line 39 of file KFPTrack.h.

References fP.

Referenced by KFParticleSIMD::Create(), KFParticle::KFParticle(), and KFParticleSIMD::KFParticleSIMD().

{position[0] = fP[0]; position[1] = fP[1]; position[2] = fP[2];}

void KFPTrack::XvYvZv ( double *  position ) const

inline

Definition at line 41 of file KFPTrack.h.

References fP.

{position[0] = fP[0]; position[1] = fP[1]; position[2] = fP[2];}

Member Data Documentation

float KFPTrack::fC[21]

private

Definition at line 115 of file KFPTrack.h.

Referenced by GetCovariance(), GetCovarianceMatrix(), GetCovarianceXYZPxPyPz(), GetCovMatrix(), RotateXY(), and SetCovarianceMatrix().

float KFPTrack::fChi2

private

Definition at line 116 of file KFPTrack.h.

Referenced by GetChi2(), GetChi2perNDF(), and SetChi2().

int KFPTrack::fId

private

Definition at line 120 of file KFPTrack.h.

Referenced by Id(), and SetId().

int KFPTrack::fNDF

private

Definition at line 118 of file KFPTrack.h.

Referenced by GetChi2perNDF(), GetNDF(), and SetNDF().

float KFPTrack::fP[6]

private

Definition at line 114 of file KFPTrack.h.

Referenced by GetP(), GetParameter(), GetPt(), GetPx(), GetPxPyPz(), GetPy(), GetPz(), GetTrack(), GetX(), GetXYZ(), GetXYZPxPyPz(), GetY(), GetZ(), PxPyPz(), SetParameters(), SetPx(), SetPxPyPz(), SetPy(), SetPz(), SetX(), SetXYZ(), SetY(), SetZ(), and XvYvZv().

int KFPTrack::fQ

private

Definition at line 117 of file KFPTrack.h.

Referenced by Charge(), and SetCharge().

int KFPTrack::ID

private

Definition at line 109 of file KFPTrack.h.

Referenced by GetID(), and SetID().

---

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

• KFPTrack.h
• KFPTrack.cxx