KFParticle.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//---------------------------------------------------------------------------------
// The KFParticle class
// .
// @author  S.Gorbunov, I.Kisel
// @version 1.0
// @since   13.05.07
//
// Class to reconstruct and store the decayed particle parameters.
// The method is described in CBM-SOFT note 2007-003,
// ``Reconstruction of decayed particles based on the Kalman filter'',
// http://www.gsi.de/documents/DOC-2007-May-14-1.pdf
//
// This class is ALICE interface to general mathematics in KFParticleBase
//
//  -= Copyright &copy ALICE HLT and CBM L1 Groups =-
//_________________________________________________________________________________

//#define NonhomogeneousField
#define HomogeneousField

#ifndef ALIKFPARTICLE_H
#define ALIKFPARTICLE_H

#include "KFParticleBase.h"
// #include "TMath.h"
#include <cmath>

#ifdef HomogeneousField
class KFPTrack;
class KFPVertex;
#endif

#ifdef NonhomogeneousField
class CbmKFTrackInterface;
class CbmKFVertexInterface;
#endif

class KFParticle : public KFParticleBase {

 public:
  //*
  //*  INITIALIZATION
  //*

  //* Set magnetic field for all particles
#ifdef HomogeneousField
  static void SetField(Double_t Bz);
#endif
  //* Constructor (empty)

  KFParticle() : KFParticleBase() { ; }

  //* Destructor (empty)

  ~KFParticle() { ; }

  //* Construction of mother particle by its 2-3-4 daughters

  KFParticle(const KFParticle &d1, const KFParticle &d2, Bool_t gamma = false);

  KFParticle(const KFParticle &d1, const KFParticle &d2, const KFParticle &d3);

  KFParticle(const KFParticle &d1, const KFParticle &d2, const KFParticle &d3, const KFParticle &d4);

  //* Initialisation from "cartesian" coordinates ( X Y Z Px Py Pz )
  //* Parameters, covariance matrix, charge and PID hypothesis should be provided

  void Create(const Double_t Param[], const Double_t Cov[], Int_t Charge, Double_t mass /*Int_t PID*/);

#ifdef HomogeneousField
  //* Initialisation from ALICE track, PID hypothesis shoould be provided

  KFParticle(const KFPTrack &track, Int_t PID);

  //* Initialisation from VVertex

  KFParticle(const KFPVertex &vertex);
#endif

#ifdef NonhomogeneousField
  KFParticle(CbmKFTrackInterface *Track, Double_t *z0 = 0, Int_t *qHypo = 0, Int_t *PID = 0);
  KFParticle(CbmKFVertexInterface &vertex);
#endif

  void CleanDaughtersId() { fDaughtersIds.clear(); }
  void SetNDaughters(int n) { fDaughtersIds.reserve(n); }
  void AddDaughter(int id) { fDaughtersIds.push_back(id); }

  //* Initialise covariance matrix and set current parameters to 0.0

  void Initialize();

  //* Set decay vertex parameters for linearisation

  void SetVtxGuess(Double_t x, Double_t y, Double_t z);

  //*
  //*  ACCESSORS
  //*

  //* Simple accessors

  Double_t GetX() const;    //* x of current position
  Double_t GetY() const;    //* y of current position
  Double_t GetZ() const;    //* z of current position
  Double_t GetPx() const;   //* x-compoment of 3-momentum
  Double_t GetPy() const;   //* y-compoment of 3-momentum
  Double_t GetPz() const;   //* z-compoment of 3-momentum
  Double_t GetE() const;    //* energy
  Double_t GetS() const;    //* decay length / momentum
  Int_t GetQ() const;       //* charge
  Double_t GetChi2() const; //* chi^2
  Int_t GetNDF() const;     //* Number of Degrees of Freedom

  Bool_t GetAtProductionVertex() const { return fAtProductionVertex; }
  void SetAtProductionVertex(Bool_t b) { fAtProductionVertex = b; }

#ifdef NonhomogeneousField
  float *GetFieldCoeff() { return fieldRegion; }
  void SetFieldCoeff(float c, int i) { fieldRegion[i] = c; }
#endif

  const Double_t &X() const { return fP[0]; }
  const Double_t &Y() const { return fP[1]; }
  const Double_t &Z() const { return fP[2]; }
  const Double_t &Px() const { return fP[3]; }
  const Double_t &Py() const { return fP[4]; }
  const Double_t &Pz() const { return fP[5]; }
  const Double_t &E() const { return fP[6]; }
  const Double_t &S() const { return fP[7]; }
  const Int_t &Q() const { return fQ; }
  const Double_t &Chi2() const { return fChi2; }
  const Int_t &NDF() const { return fNDF; }

  Double_t GetParameter(int i) const;
  Double_t GetCovariance(int i) const;
  Double_t GetCovariance(int i, int j) const;

  //* Accessors with calculations, value returned w/o error flag

  Double_t GetP() const;             //* momentum
  Double_t GetPt() const;            //* transverse momentum
  Double_t GetEta() const;           //* pseudorapidity
  Double_t GetPhi() const;           //* phi
  Double_t GetMomentum() const;      //* momentum (same as GetP() )
  Double_t GetMass() const;          //* mass
  Double_t GetDecayLength() const;   //* decay length
  Double_t GetDecayLengthXY() const; //* decay length in XY
  Double_t GetLifeTime() const;      //* life time
  Double_t GetR() const;             //* distance to the origin

  //* Accessors to estimated errors

  Double_t GetErrX() const;             //* x of current position
  Double_t GetErrY() const;             //* y of current position
  Double_t GetErrZ() const;             //* z of current position
  Double_t GetErrPx() const;            //* x-compoment of 3-momentum
  Double_t GetErrPy() const;            //* y-compoment of 3-momentum
  Double_t GetErrPz() const;            //* z-compoment of 3-momentum
  Double_t GetErrE() const;             //* energy
  Double_t GetErrS() const;             //* decay length / momentum
  Double_t GetErrP() const;             //* momentum
  Double_t GetErrPt() const;            //* transverse momentum
  Double_t GetErrEta() const;           //* pseudorapidity
  Double_t GetErrPhi() const;           //* phi
  Double_t GetErrMomentum() const;      //* momentum
  Double_t GetErrMass() const;          //* mass
  Double_t GetErrDecayLength() const;   //* decay length
  Double_t GetErrDecayLengthXY() const; //* decay length in XY
  Double_t GetErrLifeTime() const;      //* life time
  Double_t GetErrR() const;             //* distance to the origin

  //* Accessors with calculations( &value, &estimated sigma )
  //* error flag returned (0 means no error during calculations)

  int GetP(Double_t &P, Double_t &SigmaP) const;             //* momentum
  int GetPt(Double_t &Pt, Double_t &SigmaPt) const;          //* transverse momentum
  int GetEta(Double_t &Eta, Double_t &SigmaEta) const;       //* pseudorapidity
  int GetPhi(Double_t &Phi, Double_t &SigmaPhi) const;       //* phi
  int GetMomentum(Double_t &P, Double_t &SigmaP) const;      //* momentum
  int GetMass(Double_t &M, Double_t &SigmaM) const;          //* mass
  int GetDecayLength(Double_t &L, Double_t &SigmaL) const;   //* decay length
  int GetDecayLengthXY(Double_t &L, Double_t &SigmaL) const; //* decay length in XY
  int GetLifeTime(Double_t &T, Double_t &SigmaT) const;      //* life time
  int GetR(Double_t &R, Double_t &SigmaR) const;             //* R
  Double_t GetRapidity() const { return 0.5 * log((fP[6] + fP[5]) / (fP[6] - fP[5])); }
  Double_t GetTheta() const { return atan2(GetPt(), fP[5]); }

  //*
  //*  MODIFIERS
  //*

  Double_t &X();
  Double_t &Y();
  Double_t &Z();
  Double_t &Px();
  Double_t &Py();
  Double_t &Pz();
  Double_t &E();
  Double_t &S();
  Int_t &Q();
  Double_t &Chi2();
  Int_t &NDF();

  Double_t &Parameter(int i);
  Double_t &Covariance(int i);
  Double_t &Covariance(int i, int j);
  Double_t *Parameters();
  Double_t *CovarianceMatrix();

  //*
  //* CONSTRUCTION OF THE PARTICLE BY ITS DAUGHTERS AND MOTHER
  //* USING THE KALMAN FILTER METHOD
  //*

  //* Add daughter to the particle

  void AddDaughter(const KFParticle &Daughter);

  //* Add daughter via += operator: ex.{ D0; D0+=Pion; D0+= Kaon; }

  void operator+=(const KFParticle &Daughter);

  //* Set production vertex

  void SetProductionVertex(const KFParticle &Vtx);

  //* Set mass constraint

  void SetMassConstraint(Double_t Mass, Double_t SigmaMass = 0);

  //* Set no decay length for resonances

  void SetNoDecayLength();

  //* Everything in one go

  void Construct(const KFParticle *vDaughters[], int nDaughters, const KFParticle *ProdVtx = nullptr, Double_t Mass = -1, Bool_t IsConstrained = 0);

  //*
  //*                   TRANSPORT
  //*
  //*  ( main transportation parameter is S = SignedPath/Momentum )
  //*  ( parameters of decay & production vertices are stored locally )
  //*

  //* Transport the particle to its decay vertex

  void TransportToDecayVertex();

  //* Transport the particle to its production vertex

  void TransportToProductionVertex();

  //* Transport the particle close to xyz[] point

  void TransportToPoint(const Double_t xyz[]);

  //* Transport the particle close to VVertex
#ifdef HomogeneousField
  void TransportToVertex(const KFPVertex &v);
#endif
  //* Transport the particle close to another particle p

  void TransportToParticle(const KFParticle &p);

  //* Transport the particle on dS parameter (SignedPath/Momentum)

  void TransportToDS(Double_t dS);

  //* Get dS to a certain space point

  Double_t GetDStoPoint(const Double_t xyz[]) const;

  //* Get dS to other particle p (dSp for particle p also returned)

  void GetDStoParticle(const KFParticle &p, Double_t &DS, Double_t &DSp) const;

  //* Get dS to other particle p in XY-plane

  void GetDStoParticleXY(const KFParticleBase &p, Double_t &DS, Double_t &DSp) const;

  //*
  //* OTHER UTILITIES
  //*

  //* Calculate distance from another object [cm]

  Double_t GetDistanceFromVertex(const Double_t vtx[]) const;
  Double_t GetDistanceFromVertex(const KFParticle &Vtx) const;
#ifdef HomogeneousField
  Double_t GetDistanceFromVertex(const KFPVertex &Vtx) const;
#endif
  Double_t GetDistanceFromParticle(const KFParticle &p) const;

  //* Calculate sqrt(Chi2/ndf) deviation from another object
  //* ( v = [xyz]-vertex, Cv=[Cxx,Cxy,Cyy,Cxz,Cyz,Czz]-covariance matrix )

  Double_t GetDeviationFromVertex(const Double_t v[], const Double_t Cv[] = 0) const;
  Double_t GetDeviationFromVertex(const KFParticle &Vtx) const;
#ifdef HomogeneousField
  Double_t GetDeviationFromVertex(const KFPVertex &Vtx) const;
#endif
  Double_t GetDeviationFromParticle(const KFParticle &p) const;

  //* Calculate distance from another object [cm] in XY-plane

  Bool_t GetDistanceFromVertexXY(const Double_t vtx[], Double_t &val, Double_t &err) const;
  Bool_t GetDistanceFromVertexXY(const Double_t vtx[], const Double_t Cv[], Double_t &val, Double_t &err) const;
  Bool_t GetDistanceFromVertexXY(const KFParticle &Vtx, Double_t &val, Double_t &err) const;
#ifdef HomogeneousField
  Bool_t GetDistanceFromVertexXY(const KFPVertex &Vtx, Double_t &val, Double_t &err) const;
#endif

  Double_t GetDistanceFromVertexXY(const Double_t vtx[]) const;
  Double_t GetDistanceFromVertexXY(const KFParticle &Vtx) const;
#ifdef HomogeneousField
  Double_t GetDistanceFromVertexXY(const KFPVertex &Vtx) const;
#endif
  Double_t GetDistanceFromParticleXY(const KFParticle &p) const;

  //* Calculate sqrt(Chi2/ndf) deviation from another object in XY plane
  //* ( v = [xyz]-vertex, Cv=[Cxx,Cxy,Cyy,Cxz,Cyz,Czz]-covariance matrix )

  Double_t GetDeviationFromVertexXY(const Double_t v[], const Double_t Cv[] = 0) const;
  Double_t GetDeviationFromVertexXY(const KFParticle &Vtx) const;
#ifdef HomogeneousField
  Double_t GetDeviationFromVertexXY(const KFPVertex &Vtx) const;
#endif
  Double_t GetDeviationFromParticleXY(const KFParticle &p) const;

  //* Calculate opennig angle between two particles

  Double_t GetAngle(const KFParticle &p) const;
  Double_t GetAngleXY(const KFParticle &p) const;
  Double_t GetAngleRZ(const KFParticle &p) const;

  //* Subtract the particle from the vertex

  void SubtractFromVertex(KFParticle &v) const;

  //* Special method for creating gammas

  void ConstructGamma(const KFParticle &daughter1, const KFParticle &daughter2);

  // * Pseudo Proper Time of decay = (r*pt) / |pt| * M/|pt|
  // @primVertex - primary vertex
  // @mass - mass of the mother particle (in the case of "Hb -> JPsi" it would be JPsi mass)
  // @*timeErr2 - squared error of the decay time. If timeErr2 = 0 it isn't calculated
  Double_t GetPseudoProperDecayTime(const KFParticle &primVertex, const Double_t &mass, Double_t *timeErr2 = 0) const;

  void GetFieldValue(const Double_t xyz[], Double_t B[]) const;

 protected:
  //*
  //*  INTERNAL STUFF
  //*

  //* Method to access ALICE field
#ifdef HomogeneousField
  static Double_t GetFieldAlice();
#endif
  //* Other methods required by the abstract KFParticleBase class

  void GetDStoParticle(const KFParticleBase &p, Double_t &DS, Double_t &DSp) const;
  void Transport(Double_t dS, Double_t P[], Double_t C[]) const;
  static void GetExternalTrackParam(const KFParticleBase &p, Double_t &X, Double_t &Alpha, Double_t P[5]);

  // void GetDStoParticleALICE( const KFParticleBase &p, Double_t &DS, Double_t &DS1 ) const;

 private:
#ifdef HomogeneousField
  static Double_t fgBz; //* Bz compoment of the magnetic field
#endif
#ifdef NonhomogeneousField
  float fieldRegion[10];
#endif

  //  ClassDef( KFParticle, 1 );
};

//---------------------------------------------------------------------
//
//     Inline implementation of the KFParticle methods
//
//---------------------------------------------------------------------

#ifdef HomogeneousField
inline void KFParticle::SetField(Double_t Bz)
{
  fgBz = Bz;
}
#endif

inline KFParticle::KFParticle(const KFParticle &d1, const KFParticle &d2, const KFParticle &d3)
{
  KFParticle mother;
  mother += d1;
  mother += d2;
  mother += d3;
  *this = mother;
}

inline KFParticle::KFParticle(const KFParticle &d1, const KFParticle &d2, const KFParticle &d3, const KFParticle &d4)
{
  KFParticle mother;
  mother += d1;
  mother += d2;
  mother += d3;
  mother += d4;
  *this = mother;
}

inline void KFParticle::Initialize()
{
  KFParticleBase::Initialize();
}

inline void KFParticle::SetVtxGuess(Double_t x, Double_t y, Double_t z)
{
  KFParticleBase::SetVtxGuess(x, y, z);
}

inline Double_t KFParticle::GetX() const
{
  return KFParticleBase::GetX();
}

inline Double_t KFParticle::GetY() const
{
  return KFParticleBase::GetY();
}

inline Double_t KFParticle::GetZ() const
{
  return KFParticleBase::GetZ();
}

inline Double_t KFParticle::GetPx() const
{
  return KFParticleBase::GetPx();
}

inline Double_t KFParticle::GetPy() const
{
  return KFParticleBase::GetPy();
}

inline Double_t KFParticle::GetPz() const
{
  return KFParticleBase::GetPz();
}

inline Double_t KFParticle::GetE() const
{
  return KFParticleBase::GetE();
}

inline Double_t KFParticle::GetS() const
{
  return KFParticleBase::GetS();
}

inline Int_t KFParticle::GetQ() const
{
  return KFParticleBase::GetQ();
}

inline Double_t KFParticle::GetChi2() const
{
  return KFParticleBase::GetChi2();
}

inline Int_t KFParticle::GetNDF() const
{
  return KFParticleBase::GetNDF();
}

inline Double_t KFParticle::GetParameter(int i) const
{
  return KFParticleBase::GetParameter(i);
}

inline Double_t KFParticle::GetCovariance(int i) const
{
  return KFParticleBase::GetCovariance(i);
}

inline Double_t KFParticle::GetCovariance(int i, int j) const
{
  return KFParticleBase::GetCovariance(i, j);
}

inline Double_t KFParticle::GetP() const
{
  Double_t par, err;
  if (KFParticleBase::GetMomentum(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetPt() const
{
  Double_t par, err;
  if (KFParticleBase::GetPt(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetEta() const
{
  Double_t par, err;
  if (KFParticleBase::GetEta(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetPhi() const
{
  Double_t par, err;
  if (KFParticleBase::GetPhi(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetMomentum() const
{
  Double_t par, err;
  if (KFParticleBase::GetMomentum(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetMass() const
{
  Double_t par, err;
  if (KFParticleBase::GetMass(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetDecayLength() const
{
  Double_t par, err;
  if (KFParticleBase::GetDecayLength(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetDecayLengthXY() const
{
  Double_t par, err;
  if (KFParticleBase::GetDecayLengthXY(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetLifeTime() const
{
  Double_t par, err;
  if (KFParticleBase::GetLifeTime(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetR() const
{
  Double_t par, err;
  if (KFParticleBase::GetR(par, err))
    return 0;
  else
    return par;
}

inline Double_t KFParticle::GetErrX() const
{
  return sqrt(fabs(GetCovariance(0, 0)));
}

inline Double_t KFParticle::GetErrY() const
{
  return sqrt(fabs(GetCovariance(1, 1)));
}

inline Double_t KFParticle::GetErrZ() const
{
  return sqrt(fabs(GetCovariance(2, 2)));
}

inline Double_t KFParticle::GetErrPx() const
{
  return sqrt(fabs(GetCovariance(3, 3)));
}

inline Double_t KFParticle::GetErrPy() const
{
  return sqrt(fabs(GetCovariance(4, 4)));
}

inline Double_t KFParticle::GetErrPz() const
{
  return sqrt(fabs(GetCovariance(5, 5)));
}

inline Double_t KFParticle::GetErrE() const
{
  return sqrt(fabs(GetCovariance(6, 6)));
}

inline Double_t KFParticle::GetErrS() const
{
  return sqrt(fabs(GetCovariance(7, 7)));
}

inline Double_t KFParticle::GetErrP() const
{
  Double_t par, err;
  if (KFParticleBase::GetMomentum(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrPt() const
{
  Double_t par, err;
  if (KFParticleBase::GetPt(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrEta() const
{
  Double_t par, err;
  if (KFParticleBase::GetEta(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrPhi() const
{
  Double_t par, err;
  if (KFParticleBase::GetPhi(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrMomentum() const
{
  Double_t par, err;
  if (KFParticleBase::GetMomentum(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrMass() const
{
  Double_t par, err;
  if (KFParticleBase::GetMass(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrDecayLength() const
{
  Double_t par, err;
  if (KFParticleBase::GetDecayLength(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrDecayLengthXY() const
{
  Double_t par, err;
  if (KFParticleBase::GetDecayLengthXY(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrLifeTime() const
{
  Double_t par, err;
  if (KFParticleBase::GetLifeTime(par, err))
    return 1.e10;
  else
    return err;
}

inline Double_t KFParticle::GetErrR() const
{
  Double_t par, err;
  if (KFParticleBase::GetR(par, err))
    return 1.e10;
  else
    return err;
}

inline int KFParticle::GetP(Double_t &P, Double_t &SigmaP) const
{
  return KFParticleBase::GetMomentum(P, SigmaP);
}

inline int KFParticle::GetPt(Double_t &Pt, Double_t &SigmaPt) const
{
  return KFParticleBase::GetPt(Pt, SigmaPt);
}

inline int KFParticle::GetEta(Double_t &Eta, Double_t &SigmaEta) const
{
  return KFParticleBase::GetEta(Eta, SigmaEta);
}

inline int KFParticle::GetPhi(Double_t &Phi, Double_t &SigmaPhi) const
{
  return KFParticleBase::GetPhi(Phi, SigmaPhi);
}

inline int KFParticle::GetMomentum(Double_t &P, Double_t &SigmaP) const
{
  return KFParticleBase::GetMomentum(P, SigmaP);
}

inline int KFParticle::GetMass(Double_t &M, Double_t &SigmaM) const
{
  return KFParticleBase::GetMass(M, SigmaM);
}

inline int KFParticle::GetDecayLength(Double_t &L, Double_t &SigmaL) const
{
  return KFParticleBase::GetDecayLength(L, SigmaL);
}

inline int KFParticle::GetDecayLengthXY(Double_t &L, Double_t &SigmaL) const
{
  return KFParticleBase::GetDecayLengthXY(L, SigmaL);
}

inline int KFParticle::GetLifeTime(Double_t &T, Double_t &SigmaT) const
{
  return KFParticleBase::GetLifeTime(T, SigmaT);
}

inline int KFParticle::GetR(Double_t &R, Double_t &SigmaR) const
{
  return KFParticleBase::GetR(R, SigmaR);
}

inline Double_t &KFParticle::X()
{
  return KFParticleBase::X();
}

inline Double_t &KFParticle::Y()
{
  return KFParticleBase::Y();
}

inline Double_t &KFParticle::Z()
{
  return KFParticleBase::Z();
}

inline Double_t &KFParticle::Px()
{
  return KFParticleBase::Px();
}

inline Double_t &KFParticle::Py()
{
  return KFParticleBase::Py();
}

inline Double_t &KFParticle::Pz()
{
  return KFParticleBase::Pz();
}

inline Double_t &KFParticle::E()
{
  return KFParticleBase::E();
}

inline Double_t &KFParticle::S()
{
  return KFParticleBase::S();
}

inline Int_t &KFParticle::Q()
{
  return KFParticleBase::Q();
}

inline Double_t &KFParticle::Chi2()
{
  return KFParticleBase::Chi2();
}

inline Int_t &KFParticle::NDF()
{
  return KFParticleBase::NDF();
}

inline Double_t &KFParticle::Parameter(int i)
{
  return KFParticleBase::Parameter(i);
}

inline Double_t &KFParticle::Covariance(int i)
{
  return KFParticleBase::Covariance(i);
}

inline Double_t &KFParticle::Covariance(int i, int j)
{
  return KFParticleBase::Covariance(i, j);
}

inline Double_t *KFParticle::Parameters()
{
  return fP;
}

inline Double_t *KFParticle::CovarianceMatrix()
{
  return fC;
}

inline void KFParticle::operator+=(const KFParticle &Daughter)
{
  KFParticleBase::operator+=(Daughter);
}

inline void KFParticle::AddDaughter(const KFParticle &Daughter)
{
  KFParticleBase::AddDaughter(Daughter);
}

inline void KFParticle::SetProductionVertex(const KFParticle &Vtx)
{
  KFParticleBase::SetProductionVertex(Vtx);
}

inline void KFParticle::SetMassConstraint(Double_t Mass, Double_t SigmaMass)
{
  KFParticleBase::SetMassConstraint(Mass, SigmaMass);
}

inline void KFParticle::SetNoDecayLength()
{
  KFParticleBase::SetNoDecayLength();
}

inline void KFParticle::Construct(const KFParticle *vDaughters[], int nDaughters, const KFParticle *ProdVtx, Double_t Mass, Bool_t IsConstrained)
{
  KFParticleBase::Construct((const KFParticleBase **)vDaughters, nDaughters, (const KFParticleBase *)ProdVtx, Mass, IsConstrained);
}

inline void KFParticle::TransportToDecayVertex()
{
  KFParticleBase::TransportToDecayVertex();
}

inline void KFParticle::TransportToProductionVertex()
{
  KFParticleBase::TransportToProductionVertex();
}

inline void KFParticle::TransportToPoint(const Double_t xyz[])
{
  TransportToDS(GetDStoPoint(xyz));
}
#ifdef HomogeneousField
inline void KFParticle::TransportToVertex(const KFPVertex &v)
{
  TransportToPoint(KFParticle(v).fP);
}
#endif
inline void KFParticle::TransportToParticle(const KFParticle &p)
{
  Double_t dS, dSp;
  GetDStoParticle(p, dS, dSp);
  TransportToDS(dS);
}

inline void KFParticle::TransportToDS(Double_t dS)
{
  KFParticleBase::TransportToDS(dS);
}

inline Double_t KFParticle::GetDStoPoint(const Double_t xyz[]) const
{
#ifdef HomogeneousField
  return KFParticleBase::GetDStoPointBz(GetFieldAlice(), xyz);
#endif
#ifdef NonhomogeneousField
  return KFParticleBase::GetDStoPointCBM(xyz);
#endif
}

inline void KFParticle::GetDStoParticle(const KFParticle &p, Double_t &DS, Double_t &DSp) const
{
  GetDStoParticleXY(p, DS, DSp);
}

inline Double_t KFParticle::GetDistanceFromVertex(const Double_t vtx[]) const
{
  return KFParticleBase::GetDistanceFromVertex(vtx);
}

inline Double_t KFParticle::GetDeviationFromVertex(const Double_t v[], const Double_t Cv[]) const
{
  return KFParticleBase::GetDeviationFromVertex(v, Cv);
}

inline Double_t KFParticle::GetDistanceFromVertex(const KFParticle &Vtx) const
{
  return KFParticleBase::GetDistanceFromVertex(Vtx);
}

inline Double_t KFParticle::GetDeviationFromVertex(const KFParticle &Vtx) const
{
  return KFParticleBase::GetDeviationFromVertex(Vtx);
}
#ifdef HomogeneousField
inline Double_t KFParticle::GetDistanceFromVertex(const KFPVertex &Vtx) const
{
  return GetDistanceFromVertex(KFParticle(Vtx));
}

inline Double_t KFParticle::GetDeviationFromVertex(const KFPVertex &Vtx) const
{
  return GetDeviationFromVertex(KFParticle(Vtx));
}
#endif
inline Double_t KFParticle::GetDistanceFromParticle(const KFParticle &p) const
{
  return KFParticleBase::GetDistanceFromParticle(p);
}

inline Double_t KFParticle::GetDeviationFromParticle(const KFParticle &p) const
{
  return KFParticleBase::GetDeviationFromParticle(p);
}

inline void KFParticle::SubtractFromVertex(KFParticle &v) const
{
  KFParticleBase::SubtractFromVertex(v);
}

#ifdef HomogeneousField
inline Double_t KFParticle::GetFieldAlice()
{
  return fgBz;
}
#endif

#ifdef HomogeneousField
inline void KFParticle::GetFieldValue(const Double_t * /*xyz*/, Double_t B[]) const
{
  B[0] = B[1] = 0;
  B[2] = GetFieldAlice();
}
#endif

#ifdef NonhomogeneousField
inline void KFParticle::GetFieldValue(const Double_t xyz[], Double_t B[]) const
{
  //   FairField *MF = CbmKF::Instance()->GetMagneticField();
  //   MF->GetFieldValue( xyz, B );
}
#endif

inline void KFParticle::GetDStoParticle(const KFParticleBase &p, Double_t &DS, Double_t &DSp) const
{
  GetDStoParticleXY(p, DS, DSp);
}

inline void KFParticle::GetDStoParticleXY(const KFParticleBase &p, Double_t &DS, Double_t &DSp) const
{
#ifdef HomogeneousField
  KFParticleBase::GetDStoParticleBz(GetFieldAlice(), p, DS, DSp);
#endif
#ifdef NonhomogeneousField
  KFParticleBase::GetDStoParticleCBM(p, DS, DSp);
#endif
  // GetDStoParticleALICE( p, DS, DSp ) ;
}

inline void KFParticle::Transport(Double_t dS, Double_t P[], Double_t C[]) const
{
#ifdef HomogeneousField
  KFParticleBase::TransportBz(GetFieldAlice(), dS, P, C);
#endif
#ifdef NonhomogeneousField
  KFParticleBase::TransportCBM(dS, P, C);
#endif
}

inline void KFParticle::ConstructGamma(const KFParticle &daughter1, const KFParticle &daughter2)
{
#ifdef HomogeneousField
  KFParticleBase::ConstructGammaBz(daughter1, daughter2, GetFieldAlice());
#endif
}

#endif