cbmroot/KfTrackParam_8h_source.html

/* Copyright (C) 2007-2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt

   SPDX-License-Identifier: GPL-3.0-only

   Authors: Igor Kulakov [committer], Maksym Zyzak, Sergei Zharko */


#ifndef CA_CORE_CaTrackParam_h

#define CA_CORE_CaTrackParam_h 1


#include "KfDefs.h"

#include "KfSimd.h"

#include "KfUtils.h"


#include <boost/serialization/access.hpp>


#include <string>


namespace cbm::algo::kf

{

  template<typename T>

  class alignas(VcMemAlign) TrackParamBase {


   public:

    friend class boost::serialization::access;


    static constexpr int kNtrackParam{7};

    static constexpr int kNcovParam{(kNtrackParam) * (kNtrackParam + 1) / 2};


    typedef std::array<T, kNcovParam> CovMatrix_t;


    TrackParamBase() = default;


    template<typename T1>


    TrackParamBase(const TrackParamBase<T1>& tr)

    {

      Set(tr);

    }


    template<typename T1>


    void Set(const TrackParamBase<T1>& Tb)

    {

      CopyBase<T1, true, true>(0, Tb, 0);

    }


    void Set(const TrackParamBase<fvec>& Tb, const int ib) { CopyBase<fvec, true, false>(0, Tb, ib); }


    template<typename T1>


    void SetOneEntry(const int ia, const TrackParamBase<T1>& Tb)

    {

      CopyBase<T1, false, true>(ia, Tb, 0);

    }


    void SetOneEntry(const int ia, const TrackParamBase<fvec>& Tb, const int ib)

    {

      CopyBase<fvec, false, false>(ia, Tb, ib);

    }


    T Z() const { return fZ; }


    T X() const { return fX; }


    T Y() const { return fY; }


    T Tx() const { return fTx; }


    T Ty() const { return fTy; }


    T Qp() const { return fQp; }


    T Time() const { return fT; }


    T Vi() const { return fVi; }


    T ChiSq() const { return fChiSq; }


    T Ndf() const { return fNdf; }


    T ChiSqTime() const { return fChiSqTime; }


    T NdfTime() const { return fNdfTime; }


    T C(int i, int j) const

    {

      int ind = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      return fCovMatrix[ind];

    }


    template<int i, int j>


    T C() const

    {

      constexpr int ind = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      return fCovMatrix[ind];

    }


    T C00() const { return C<0, 0>(); }

    T C01() const { return C<0, 1>(); }

    T C02() const { return C<0, 2>(); }

    T C03() const { return C<0, 3>(); }

    T C04() const { return C<0, 4>(); }

    T C05() const { return C<0, 5>(); }

    T C06() const { return C<0, 6>(); }


    T C10() const { return C<1, 0>(); }

    T C11() const { return C<1, 1>(); }

    T C12() const { return C<1, 2>(); }

    T C13() const { return C<1, 3>(); }

    T C14() const { return C<1, 4>(); }

    T C15() const { return C<1, 5>(); }

    T C16() const { return C<1, 6>(); }


    T C20() const { return C<2, 0>(); }

    T C21() const { return C<2, 1>(); }

    T C22() const { return C<2, 2>(); }

    T C23() const { return C<2, 3>(); }

    T C24() const { return C<2, 4>(); }

    T C25() const { return C<2, 5>(); }

    T C26() const { return C<2, 6>(); }


    T C30() const { return C<3, 0>(); }

    T C31() const { return C<3, 1>(); }

    T C32() const { return C<3, 2>(); }

    T C33() const { return C<3, 3>(); }

    T C34() const { return C<3, 4>(); }

    T C35() const { return C<3, 5>(); }

    T C36() const { return C<3, 6>(); }


    T C40() const { return C<4, 0>(); }

    T C41() const { return C<4, 1>(); }

    T C42() const { return C<4, 2>(); }

    T C43() const { return C<4, 3>(); }

    T C44() const { return C<4, 4>(); }

    T C45() const { return C<4, 5>(); }

    T C46() const { return C<4, 6>(); }


    T C50() const { return C<5, 0>(); }

    T C51() const { return C<5, 1>(); }

    T C52() const { return C<5, 2>(); }

    T C53() const { return C<5, 3>(); }

    T C54() const { return C<5, 4>(); }

    T C55() const { return C<5, 5>(); }

    T C56() const { return C<5, 6>(); }


    T C60() const { return C<6, 0>(); }

    T C61() const { return C<6, 1>(); }

    T C62() const { return C<6, 2>(); }

    T C63() const { return C<6, 3>(); }

    T C64() const { return C<6, 4>(); }

    T C65() const { return C<6, 5>(); }

    T C66() const { return C<6, 6>(); }


    T GetZ() const { return fZ; }


    T GetX() const { return fX; }


    T GetXError() const { return sqrt(C00()); }


    T GetY() const { return fY; }


    T GetYError() const { return sqrt(C11()); }


    T GetTx() const { return fTx; }


    T GetTxError() const { return sqrt(C22()); }


    T GetTy() const { return fTy; }


    T GetTyError() const { return sqrt(C33()); }


    T GetQp() const { return fQp; }


    T GetQpError() const { return sqrt(C44()); }


    T GetTime() const { return fT; }


    T GetTimeError() const { return sqrt(C55()); }


    T GetVi() const { return fVi; }


    T GetViError() const { return sqrt(C66()); }


    const CovMatrix_t& GetCovMatrix() const { return fCovMatrix; }


    T GetCovariance(int i, int j) const { return C(i, j); }


    T GetChiSq() const { return fChiSq; }


    T GetNdf() const { return fNdf; }


    T GetChiSqTime() const { return fChiSqTime; }


    T GetNdfTime() const { return fNdfTime; }


    T GetCharge() const { return utils::iif(GetQp() > T(0.), T(1.), T(-1.)); }


    T GetPhi() const { return atan2(GetTy(), GetTx()); }


    T GetPhiError() const;


    T GetTheta() const { return atan(sqrt(GetTx() * GetTx() + GetTy() * GetTy())); }


    T GetThetaError() const;


    T GetP() const { return utils::iif(utils::fabs(GetQp()) > T(1.e-4), T(1.) / utils::fabs(GetQp()), T(1.e4)); }


    T GetPz() const { return GetP() / sqrt(T(1.) + GetTx() * GetTx() + GetTy() * GetTy()); }


    T GetPx() const { return GetPz() * GetTx(); }


    T GetPy() const { return GetPz() * GetTy(); }


    T GetPt() const

    {

      T t2 = GetTx() * GetTx() + GetTy() * GetTy();

      return GetP() * sqrt(t2 / (T(1.) + t2));

    }


    void SetZ(T v) { fZ = v; }


    void SetX(T v) { fX = v; }


    void SetY(T v) { fY = v; }


    void SetTx(T v) { fTx = v; }


    void SetTy(T v) { fTy = v; }


    void SetQp(T v) { fQp = v; }


    void SetTime(T v) { fT = v; }


    void SetVi(T v) { fVi = v; }


    void SetChiSq(T v) { fChiSq = v; }


    void SetNdf(T v) { fNdf = v; }


    void SetChiSqTime(T v) { fChiSqTime = v; }


    void SetNdfTime(T v) { fNdfTime = v; }


    void SetCovMatrix(const CovMatrix_t& val) { fCovMatrix = val; }


    void SetCovariance(int i, int j, T val)

    {

      int ind         = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      fCovMatrix[ind] = val;

    }


    template<int i, int j>


    void SetCovariance(T val)

    {

      constexpr int ind = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      fCovMatrix[ind]   = val;

    }


    void SetC00(T val) { SetCovariance<0, 0>(val); }

    void SetC10(T val) { SetCovariance<1, 0>(val); }

    void SetC11(T val) { SetCovariance<1, 1>(val); }

    void SetC20(T val) { SetCovariance<2, 0>(val); }

    void SetC21(T val) { SetCovariance<2, 1>(val); }

    void SetC22(T val) { SetCovariance<2, 2>(val); }

    void SetC30(T val) { SetCovariance<3, 0>(val); }

    void SetC31(T val) { SetCovariance<3, 1>(val); }

    void SetC32(T val) { SetCovariance<3, 2>(val); }

    void SetC33(T val) { SetCovariance<3, 3>(val); }

    void SetC40(T val) { SetCovariance<4, 0>(val); }

    void SetC41(T val) { SetCovariance<4, 1>(val); }

    void SetC42(T val) { SetCovariance<4, 2>(val); }

    void SetC43(T val) { SetCovariance<4, 3>(val); }

    void SetC44(T val) { SetCovariance<4, 4>(val); }

    void SetC50(T val) { SetCovariance<5, 0>(val); }

    void SetC51(T val) { SetCovariance<5, 1>(val); }

    void SetC52(T val) { SetCovariance<5, 2>(val); }

    void SetC53(T val) { SetCovariance<5, 3>(val); }

    void SetC54(T val) { SetCovariance<5, 4>(val); }

    void SetC55(T val) { SetCovariance<5, 5>(val); }

    void SetC60(T val) { SetCovariance<6, 0>(val); }

    void SetC61(T val) { SetCovariance<6, 1>(val); }

    void SetC62(T val) { SetCovariance<6, 2>(val); }

    void SetC63(T val) { SetCovariance<6, 3>(val); }

    void SetC64(T val) { SetCovariance<6, 4>(val); }

    void SetC65(T val) { SetCovariance<6, 5>(val); }

    void SetC66(T val) { SetCovariance<6, 6>(val); }


    T& Z() { return fZ; }


    T& X() { return fX; }


    T& Y() { return fY; }


    T& Tx() { return fTx; }


    T& Ty() { return fTy; }


    T& Qp() { return fQp; }


    T& Time() { return fT; }


    T& Vi() { return fVi; }


    CovMatrix_t& CovMatrix() { return fCovMatrix; }


    T& ChiSq() { return fChiSq; }


    T& Ndf() { return fNdf; }


    T& ChiSqTime() { return fChiSqTime; }


    T& NdfTime() { return fNdfTime; }


    T& C(int i, int j)

    {

      int ind = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      return fCovMatrix[ind];

    }


    template<int i, int j>


    T& C()

    {

      constexpr int ind = (j <= i) ? i * (1 + i) / 2 + j : j * (1 + j) / 2 + i;

      return fCovMatrix[ind];

    }


    T& C00() { return C<0, 0>(); }

    T& C01() { return C<0, 1>(); }

    T& C02() { return C<0, 2>(); }

    T& C03() { return C<0, 3>(); }

    T& C04() { return C<0, 4>(); }

    T& C05() { return C<0, 5>(); }

    T& C06() { return C<0, 6>(); }


    T& C10() { return C<1, 0>(); }

    T& C11() { return C<1, 1>(); }

    T& C12() { return C<1, 2>(); }

    T& C13() { return C<1, 3>(); }

    T& C14() { return C<1, 4>(); }

    T& C15() { return C<1, 5>(); }

    T& C16() { return C<1, 6>(); }


    T& C20() { return C<2, 0>(); }

    T& C21() { return C<2, 1>(); }

    T& C22() { return C<2, 2>(); }

    T& C23() { return C<2, 3>(); }

    T& C24() { return C<2, 4>(); }

    T& C25() { return C<2, 5>(); }

    T& C26() { return C<2, 6>(); }


    T& C30() { return C<3, 0>(); }

    T& C31() { return C<3, 1>(); }

    T& C32() { return C<3, 2>(); }

    T& C33() { return C<3, 3>(); }

    T& C34() { return C<3, 4>(); }

    T& C35() { return C<3, 5>(); }

    T& C36() { return C<3, 6>(); }


    T& C40() { return C<4, 0>(); }

    T& C41() { return C<4, 1>(); }

    T& C42() { return C<4, 2>(); }

    T& C43() { return C<4, 3>(); }

    T& C44() { return C<4, 4>(); }

    T& C45() { return C<4, 5>(); }

    T& C46() { return C<4, 6>(); }


    T& C50() { return C<5, 0>(); }

    T& C51() { return C<5, 1>(); }

    T& C52() { return C<5, 2>(); }

    T& C53() { return C<5, 3>(); }

    T& C54() { return C<5, 4>(); }

    T& C55() { return C<5, 5>(); }

    T& C56() { return C<5, 6>(); }


    T& C60() { return C<6, 0>(); }

    T& C61() { return C<6, 1>(); }

    T& C62() { return C<6, 2>(); }

    T& C63() { return C<6, 3>(); }

    T& C64() { return C<6, 4>(); }

    T& C65() { return C<6, 5>(); }

    T& C66() { return C<6, 6>(); }


    void ResetErrors(T c00, T c11, T c22, T c33, T c44, T c55, T c66);


    std::string ToString(int i = -1) const;


    std::string ToStringCorrelations(int i = -1) const;


    bool IsFinite(bool printWhenWrong) const;


    bool IsEntryConsistent(bool printWhenWrong, int i) const;


    bool IsConsistent(bool printWhenWrong, int nFilled) const;


    void InitVelocityRange(fscal minP);


    template<class Archive>


    void serialize(Archive& ar, const unsigned int /*version*/)

    {

      ar& fX;

      ar& fY;

      ar& fTx;

      ar& fTy;

      ar& fQp;

      ar& fZ;

      ar& fT;

      ar& fVi;

      ar& fCovMatrix;

      ar& fChiSq;

      ar& fNdf;

      ar& fChiSqTime;

      ar& fNdfTime;

    }


   private:

    template<typename T1, bool TDoAllA, bool TDoAllB>

    void CopyBase(const int ia, const TrackParamBase<T1>& Tb, const int ib);


   private:


    // Initializing parameters with NANs spoils the track fit where

    // the masked-out SIMD entries are suppressed by a multication by 0.

    // Therefore, we initialize the data members here with finite numbers.

    // For the numerical safety, with some reasonable numbers.


    // clang-format off


    CovMatrix_t fCovMatrix {1.,

                            0., 1.,

                            0., 0., 1.,

                            0., 0., 0., 1.,

                            0., 0., 0., 0., 1.,

                            0., 0., 0., 0., 0., 1.,

                            0., 0., 0., 0., 0., 0., 1.};


    // clang-format on


    T fX{0.};

    T fY{0.};

    T fZ{0.};

    T fTx{0.};

    T fTy{0.};

    T fQp{0.};

    T fT{0.};

    T fVi{0.};


    T fChiSq{0.};

    T fNdf{0.};

    T fChiSqTime{0.};

    T fNdfTime{0.};


  };  // class TrackParamBase


  template<typename T>


  class TrackParamBaseScalar : public TrackParamBase<T> {

   public:

    using TrackParamBase<T>::TrackParamBase;


    T GetEta() const { return -log(tan(this->GetTheta() * T(0.5))); }

  };


  // ---------------------------------------------------------------------------------------------------------------------


  template<typename T>

  class TrackParam : public TrackParamBaseScalar<T> {

   public:

    using TrackParamBaseScalar<T>::TrackParamBaseScalar;

  };


  template<>


  class TrackParam<fvec> : public TrackParamBase<fvec> {

   public:

    using TrackParamBase<fvec>::TrackParamBase;

  };


  typedef TrackParam<fvec> TrackParamV;

  typedef TrackParam<fscal> TrackParamS;

  typedef TrackParam<double> TrackParamD;


  template<typename T>


  inline T TrackParamBase<T>::GetPhiError() const

  {

    // phi = atan( tx / ty ); }


    T phiDdenom = GetTx() * GetTx() + GetTy() * GetTy();

    T phiDTx    = -GetTy() / phiDdenom;  // partial derivative of phi over Tx

    T phiDTy    = +GetTx() / phiDdenom;  // partial derivative of phi over Ty


    T varTx   = C22();  // variance of Tx

    T varTy   = C33();  // variance of Ty

    T covTxTy = C32();  // covariance of Tx and Ty


    T varPhi = phiDTx * phiDTx * varTx + phiDTy * phiDTy * varTy + T(2.) * phiDTx * phiDTy * covTxTy;

    return sqrt(varPhi);

  }


  template<typename T>


  inline T TrackParamBase<T>::GetThetaError() const

  {

    // theta = atan(sqrt( tx * tx + ty * ty) )


    T sumSqSlopes = GetTx() * GetTx() + GetTy() * GetTy();

    T thetaDdenom = sqrt(sumSqSlopes) * (T(1.) + sumSqSlopes);

    T thetaDTx    = GetTx() / thetaDdenom;

    T thetaDTy    = GetTy() / thetaDdenom;


    T varTx   = C22();  // variance of Tx

    T varTy   = C33();  // variance of Ty

    T covTxTy = C32();  // covariance of Tx and Ty


    T varTheta = thetaDTx * thetaDTx * varTx + thetaDTy * thetaDTy * varTy + T(2.) * thetaDTx * thetaDTy * covTxTy;


    return sqrt(varTheta);

  }


  // ---------------------------------------------------------------------------------------------------------------------

  //

  template<typename TdataA>

  template<typename TdataB, bool TDoAllA, bool TDoAllB>


  inline void TrackParamBase<TdataA>::CopyBase(const int ia, const TrackParamBase<TdataB>& Tb, const int ib)

  {

    auto copy = [&](TdataA& a, const TdataB& b) {

      utils::VecCopy<TdataA, TdataB, TDoAllA, TDoAllB>::CopyEntries(a, ia, b, ib);

    };


    copy(fX, Tb.GetX());

    copy(fY, Tb.GetY());

    copy(fTx, Tb.GetTx());

    copy(fTy, Tb.GetTy());

    copy(fQp, Tb.GetQp());

    copy(fZ, Tb.GetZ());

    copy(fT, Tb.GetTime());

    copy(fVi, Tb.GetVi());


    for (int i = 0; i < kNcovParam; ++i) {

      copy(fCovMatrix[i], Tb.GetCovMatrix()[i]);

    }


    copy(fChiSq, Tb.GetChiSq());

    copy(fNdf, Tb.GetNdf());

    copy(fChiSqTime, Tb.GetChiSqTime());

    copy(fNdfTime, Tb.GetNdfTime());


  }  // CopyBase


  // ---------------------------------------------------------------------------------------------------------------------

  //

  template<typename T>


  inline void TrackParamBase<T>::ResetErrors(T c00, T c11, T c22, T c33, T c44, T c55, T c66)

  {

    fCovMatrix.fill(0.);


    SetC00(c00);

    SetC11(c11);

    SetC22(c22);

    SetC33(c33);

    SetC44(c44);

    SetC55(c55);

    SetC66(c66);


    SetChiSq(0.);

    SetNdf(-5.);

    SetChiSqTime(0.);

    SetNdfTime(-2.);

  }


  // ---------------------------------------------------------------------------------------------------------------------

  //

  template<typename T>


  inline void TrackParamBase<T>::InitVelocityRange(fscal minP)

  {

    // initialise the velocity range with respect to the minimal momentum minP {GeV/c}

    using defs::ProtonMass;

    using defs::SpeedOfLightInv;

    fscal maxVi = sqrt(1. + (ProtonMass<double> / minP) * (ProtonMass<double> / minP)) * SpeedOfLightInv<double>;

    fscal minVi = SpeedOfLightInv<fscal>;

    fscal vmean = minVi + 0.4 * (maxVi - minVi);

    fscal dvi   = (maxVi - vmean) / 3.;

    SetVi(vmean);

    SetC66(dvi * dvi);

  }


}  // namespace cbm::algo::kf


#endif  // CA_CORE_CaTrackParam_h


ToString
std::string ToString(ECbmModuleId modId)
Definition CbmDefs.cxx:70

KfDefs.h
Common constant definitions for the Kalman Filter library.

sqrt
friend fvec sqrt(const fvec &a)
Definition KfSimdPseudo.h:118

v
fscal v[fmask::Size]
Definition KfSimdPseudo.h:4

log
friend fvec log(const fvec &a)
Definition KfSimdPseudo.h:122

KfSimd.h
Implementation selection for the SIMD utilities (VS or pseudo)

KfUtils.h

TrackParamBaseScalar
Scalar version of TrackParamBase.

cbm::algo::kf::TrackParamBaseScalar
Definition KfTrackParam.h:630

cbm::algo::kf::TrackParamBaseScalar::GetEta
T GetEta() const
Gets pseudo-rapidity.
Definition KfTrackParam.h:636

cbm::algo::kf::TrackParamBase
It is a technical base class of kf::TrackParam.

cbm::algo::kf::TrackParamBase::C06
T C06() const
Definition KfTrackParam.h:149

cbm::algo::kf::TrackParamBase::C62
T & C62()
Definition KfTrackParam.h:515

cbm::algo::kf::TrackParamBase::C05
T & C05()
Definition KfTrackParam.h:470

cbm::algo::kf::TrackParamBase::GetPx
T GetPx() const
Gets x-component of the momentum [GeV/ec].
Definition KfTrackParam.h:291

cbm::algo::kf::TrackParamBase::C00
T C00() const
Individual getters for covariance matrix elements.
Definition KfTrackParam.h:143

cbm::algo::kf::TrackParamBase::C10
T & C10()
Definition KfTrackParam.h:473

cbm::algo::kf::TrackParamBase::SetC20
void SetC20(T val)
Definition KfTrackParam.h:372

cbm::algo::kf::TrackParamBase::C63
T C63() const
Definition KfTrackParam.h:194

cbm::algo::kf::TrackParamBase::SetOneEntry
void SetOneEntry(const int ia, const TrackParamBase< fvec > &Tb, const int ib)
Definition KfTrackParam.h:74

cbm::algo::kf::TrackParamBase::C35
T C35() const
Definition KfTrackParam.h:172

cbm::algo::kf::TrackParamBase::C65
T C65() const
Definition KfTrackParam.h:196

cbm::algo::kf::TrackParamBase::C20
T & C20()
Definition KfTrackParam.h:481

cbm::algo::kf::TrackParamBase::SetChiSqTime
void SetChiSqTime(T v)
Sets Chi-square of time measurements.
Definition KfTrackParam.h:338

cbm::algo::kf::TrackParamBase::C51
T & C51()
Definition KfTrackParam.h:506

cbm::algo::kf::TrackParamBase::Ndf
T Ndf() const
Gets NDF of track fit model.
Definition KfTrackParam.h:111

cbm::algo::kf::TrackParamBase::C02
T C02() const
Definition KfTrackParam.h:145

cbm::algo::kf::TrackParamBase::Vi
T Vi() const
Gets inverse velocity [ns/cm].
Definition KfTrackParam.h:105

cbm::algo::kf::TrackParamBase::SetC55
void SetC55(T val)
Definition KfTrackParam.h:389

cbm::algo::kf::TrackParamBase::GetPhiError
T GetPhiError() const
Gets azimuthal angle error [rad].
Definition KfTrackParam.h:671

cbm::algo::kf::TrackParamBase::SetC64
void SetC64(T val)
Definition KfTrackParam.h:394

cbm::algo::kf::TrackParamBase::C41
T & C41()
Definition KfTrackParam.h:498

cbm::algo::kf::TrackParamBase::GetTxError
T GetTxError() const
Gets error of slope along x-axis.
Definition KfTrackParam.h:222

cbm::algo::kf::TrackParamBase::GetVi
T GetVi() const
Gets inverse velocity [ns/cm] in downstream direction.
Definition KfTrackParam.h:243

cbm::algo::kf::TrackParamBase::SetC11
void SetC11(T val)
Definition KfTrackParam.h:371

cbm::algo::kf::TrackParamBase::SetCovariance
void SetCovariance(int i, int j, T val)
Get covariance matrix element.
Definition KfTrackParam.h:350

cbm::algo::kf::TrackParamBase::SetTx
void SetTx(T v)
Sets slope along x-axis.
Definition KfTrackParam.h:317

cbm::algo::kf::TrackParamBase::Ndf
T & Ndf()
Reference to NDF of track fit model.
Definition KfTrackParam.h:434

cbm::algo::kf::TrackParamBase::SetC63
void SetC63(T val)
Definition KfTrackParam.h:393

cbm::algo::kf::TrackParamBase::GetXError
T GetXError() const
Gets x position error [cm].
Definition KfTrackParam.h:210

cbm::algo::kf::TrackParamBase::GetPy
T GetPy() const
Gets y-component of the momentum [GeV/ec].
Definition KfTrackParam.h:294

cbm::algo::kf::TrackParamBase::Set
void Set(const TrackParamBase< fvec > &Tb, const int ib)
Definition KfTrackParam.h:62

cbm::algo::kf::TrackParamBase::C62
T C62() const
Definition KfTrackParam.h:193

cbm::algo::kf::TrackParamBase::C15
T C15() const
Definition KfTrackParam.h:156

cbm::algo::kf::TrackParamBase::SetC51
void SetC51(T val)
Definition KfTrackParam.h:385

cbm::algo::kf::TrackParamBase::C43
T C43() const
Definition KfTrackParam.h:178

cbm::algo::kf::TrackParamBase::SetC62
void SetC62(T val)
Definition KfTrackParam.h:392

cbm::algo::kf::TrackParamBase::C32
T C32() const
Definition KfTrackParam.h:169

cbm::algo::kf::TrackParamBase::C61
T C61() const
Definition KfTrackParam.h:192

cbm::algo::kf::TrackParamBase::NdfTime
T NdfTime() const
Gets NDF of time measurements.
Definition KfTrackParam.h:117

cbm::algo::kf::TrackParamBase::GetYError
T GetYError() const
Gets y position error [cm].
Definition KfTrackParam.h:216

cbm::algo::kf::TrackParamBase::C64
T & C64()
Definition KfTrackParam.h:517

cbm::algo::kf::TrackParamBase::SetC44
void SetC44(T val)
Definition KfTrackParam.h:383

cbm::algo::kf::TrackParamBase::GetCovMatrix
const CovMatrix_t & GetCovMatrix() const
Gets covariance matrix.
Definition KfTrackParam.h:249

cbm::algo::kf::TrackParamBase::GetTy
T GetTy() const
Gets slope along y-axis.
Definition KfTrackParam.h:225

cbm::algo::kf::TrackParamBase::C60
T C60() const
Definition KfTrackParam.h:191

cbm::algo::kf::TrackParamBase::C54
T & C54()
Definition KfTrackParam.h:509

cbm::algo::kf::TrackParamBase::C40
T C40() const
Definition KfTrackParam.h:175

cbm::algo::kf::TrackParamBase::C21
T C21() const
Definition KfTrackParam.h:160

cbm::algo::kf::TrackParamBase::C42
T & C42()
Definition KfTrackParam.h:499

cbm::algo::kf::TrackParamBase::C
T C() const
Get covariance matrix element when indices are known at compile time.
Definition KfTrackParam.h:134

cbm::algo::kf::TrackParamBase::C02
T & C02()
Definition KfTrackParam.h:467

cbm::algo::kf::TrackParamBase::C33
T C33() const
Definition KfTrackParam.h:170

cbm::algo::kf::TrackParamBase::SetC66
void SetC66(T val)
Definition KfTrackParam.h:396

cbm::algo::kf::TrackParamBase::C36
T C36() const
Definition KfTrackParam.h:173

cbm::algo::kf::TrackParamBase::C16
T C16() const
Definition KfTrackParam.h:157

cbm::algo::kf::TrackParamBase::SetC30
void SetC30(T val)
Definition KfTrackParam.h:375

cbm::algo::kf::TrackParamBase::SetC32
void SetC32(T val)
Definition KfTrackParam.h:377

cbm::algo::kf::TrackParamBase::GetTimeError
T GetTimeError() const
Gets time error [ns].
Definition KfTrackParam.h:240

cbm::algo::kf::TrackParamBase::C55
T & C55()
Definition KfTrackParam.h:510

cbm::algo::kf::TrackParamBase::Tx
T & Tx()
Reference to slope along x-axis.
Definition KfTrackParam.h:413

cbm::algo::kf::TrackParamBase::C20
T C20() const
Definition KfTrackParam.h:159

cbm::algo::kf::TrackParamBase::GetCharge
T GetCharge() const
Gets charge.
Definition KfTrackParam.h:270

cbm::algo::kf::TrackParamBase::SetCovMatrix
void SetCovMatrix(const CovMatrix_t &val)
Sets covariance matrix.
Definition KfTrackParam.h:344

cbm::algo::kf::TrackParamBase::C40
T & C40()
Definition KfTrackParam.h:497

cbm::algo::kf::TrackParamBase::C52
T & C52()
Definition KfTrackParam.h:507

cbm::algo::kf::TrackParamBase::SetZ
void SetZ(T v)
Sets z position [cm].
Definition KfTrackParam.h:308

cbm::algo::kf::TrackParamBase::GetThetaError
T GetThetaError() const
Gets polar angle error [rad].
Definition KfTrackParam.h:691

cbm::algo::kf::TrackParamBase::C45
T & C45()
Definition KfTrackParam.h:502

cbm::algo::kf::TrackParamBase::C06
T & C06()
Definition KfTrackParam.h:471

cbm::algo::kf::TrackParamBase::SetC60
void SetC60(T val)
Definition KfTrackParam.h:390

cbm::algo::kf::TrackParamBase::Time
T & Time()
Reference to time [ns].
Definition KfTrackParam.h:422

cbm::algo::kf::TrackParamBase::C35
T & C35()
Definition KfTrackParam.h:494

cbm::algo::kf::TrackParamBase::C46
T C46() const
Definition KfTrackParam.h:181

cbm::algo::kf::TrackParamBase::GetPz
T GetPz() const
Gets z-component of the momentum [GeV/ec].
Definition KfTrackParam.h:288

cbm::algo::kf::TrackParamBase::C26
T & C26()
Definition KfTrackParam.h:487

cbm::algo::kf::TrackParamBase::C31
T & C31()
Definition KfTrackParam.h:490

cbm::algo::kf::TrackParamBase::GetZ
T GetZ() const
Gets z position [cm].
Definition KfTrackParam.h:204

cbm::algo::kf::TrackParamBase::C24
T & C24()
Definition KfTrackParam.h:485

cbm::algo::kf::TrackParamBase::C
T & C(int i, int j)
Get a reference to the covariance matrix element.
Definition KfTrackParam.h:446

cbm::algo::kf::TrackParamBase::SetC61
void SetC61(T val)
Definition KfTrackParam.h:391

cbm::algo::kf::TrackParamBase::C04
T C04() const
Definition KfTrackParam.h:147

cbm::algo::kf::TrackParamBase::SetC21
void SetC21(T val)
Definition KfTrackParam.h:373

cbm::algo::kf::TrackParamBase::GetP
T GetP() const
Gets momentum [GeV/ec]. For the straight tracks returns 1.e4 [GeV/ec].
Definition KfTrackParam.h:285

cbm::algo::kf::TrackParamBase::ResetErrors
void ResetErrors(T c00, T c11, T c22, T c33, T c44, T c55, T c66)
Resets variances of track parameters and chi2, ndf values.
Definition KfTrackParam.h:744

cbm::algo::kf::TrackParamBase::C24
T C24() const
Definition KfTrackParam.h:163

cbm::algo::kf::TrackParamBase::ChiSq
T ChiSq() const
Gets Chi-square of track fit model.
Definition KfTrackParam.h:108

cbm::algo::kf::TrackParamBase::Qp
T & Qp()
Reference to charge over momentum [ec/GeV].
Definition KfTrackParam.h:419

cbm::algo::kf::TrackParamBase::C05
T C05() const
Definition KfTrackParam.h:148

cbm::algo::kf::TrackParamBase::C30
T & C30()
Definition KfTrackParam.h:489

cbm::algo::kf::TrackParamBase::GetPhi
T GetPhi() const
Gets azimuthal angle [rad].
Definition KfTrackParam.h:273

cbm::algo::kf::TrackParamBase::SetC53
void SetC53(T val)
Definition KfTrackParam.h:387

cbm::algo::kf::TrackParamBase::C11
T C11() const
Definition KfTrackParam.h:152

cbm::algo::kf::TrackParamBase::Tx
T Tx() const
Gets slope along x-axis.
Definition KfTrackParam.h:93

cbm::algo::kf::TrackParamBase::GetTime
T GetTime() const
Gets time [ns].
Definition KfTrackParam.h:237

cbm::algo::kf::TrackParamBase::SetC00
void SetC00(T val)
Individual setters for covariance matrix elements.
Definition KfTrackParam.h:369

cbm::algo::kf::TrackParamBase::SetC33
void SetC33(T val)
Definition KfTrackParam.h:378

cbm::algo::kf::TrackParamBase::SetC42
void SetC42(T val)
Definition KfTrackParam.h:381

cbm::algo::kf::TrackParamBase::SetC65
void SetC65(T val)
Definition KfTrackParam.h:395

cbm::algo::kf::TrackParamBase::C03
T C03() const
Definition KfTrackParam.h:146

cbm::algo::kf::TrackParamBase::GetTyError
T GetTyError() const
Gets error of slope along y-axis.
Definition KfTrackParam.h:228

cbm::algo::kf::TrackParamBase::SetY
void SetY(T v)
Sets y position [cm].
Definition KfTrackParam.h:314

cbm::algo::kf::TrackParamBase::CopyBase
void CopyBase(const int ia, const TrackParamBase< T1 > &Tb, const int ib)
Copies all/one entries from the other class.

cbm::algo::kf::TrackParamBase::C31
T C31() const
Definition KfTrackParam.h:168

cbm::algo::kf::TrackParamBase::C50
T C50() const
Definition KfTrackParam.h:183

cbm::algo::kf::TrackParamBase::SetChiSq
void SetChiSq(T v)
Sets Chi-square of track fit model.
Definition KfTrackParam.h:332

cbm::algo::kf::TrackParamBase::X
T X() const
Gets x position [cm].
Definition KfTrackParam.h:87

cbm::algo::kf::TrackParamBase::SetC43
void SetC43(T val)
Definition KfTrackParam.h:382

cbm::algo::kf::TrackParamBase::Z
T & Z()
Reference to z position [cm].
Definition KfTrackParam.h:404

cbm::algo::kf::TrackParamBase::C64
T C64() const
Definition KfTrackParam.h:195

cbm::algo::kf::TrackParamBase::GetTheta
T GetTheta() const
Gets polar angle [rad].
Definition KfTrackParam.h:279

cbm::algo::kf::TrackParamBase::C22
T & C22()
Definition KfTrackParam.h:483

cbm::algo::kf::TrackParamBase::serialize
void serialize(Archive &ar, const unsigned int)
Definition KfTrackParam.h:560

cbm::algo::kf::TrackParamBase::SetVi
void SetVi(T v)
Sets inverse velocity [ns/cm].
Definition KfTrackParam.h:329

cbm::algo::kf::TrackParamBase::ChiSqTime
T ChiSqTime() const
Gets Chi-square of time measurements.
Definition KfTrackParam.h:114

cbm::algo::kf::TrackParamBase::GetTx
T GetTx() const
Gets slope along x-axis.
Definition KfTrackParam.h:219

cbm::algo::kf::TrackParamBase::C12
T & C12()
Definition KfTrackParam.h:475

cbm::algo::kf::TrackParamBase::ChiSq
T & ChiSq()
Reference to Chi-square of track fit model.
Definition KfTrackParam.h:431

cbm::algo::kf::TrackParamBase::Ty
T & Ty()
Reference to slope along y-axis.
Definition KfTrackParam.h:416

cbm::algo::kf::TrackParamBase::NdfTime
T & NdfTime()
Reference to NDF of time measurements.
Definition KfTrackParam.h:440

cbm::algo::kf::TrackParamBase::C13
T & C13()
Definition KfTrackParam.h:476

cbm::algo::kf::TrackParamBase::C14
T C14() const
Definition KfTrackParam.h:155

cbm::algo::kf::TrackParamBase::C11
T & C11()
Definition KfTrackParam.h:474

cbm::algo::kf::TrackParamBase::C
T & C()
Get a reference to the covariance matrix element when indices are known at compile time.
Definition KfTrackParam.h:457

cbm::algo::kf::TrackParamBase::GetViError
T GetViError() const
Gets inverse velocity error [ns/cm].
Definition KfTrackParam.h:246

cbm::algo::kf::TrackParamBase::SetNdf
void SetNdf(T v)
Sets NDF of track fit model.
Definition KfTrackParam.h:335

cbm::algo::kf::TrackParamBase::SetC41
void SetC41(T val)
Definition KfTrackParam.h:380

cbm::algo::kf::TrackParamBase::C65
T & C65()
Definition KfTrackParam.h:518

cbm::algo::kf::TrackParamBase::SetC31
void SetC31(T val)
Definition KfTrackParam.h:376

cbm::algo::kf::TrackParamBase::C13
T C13() const
Definition KfTrackParam.h:154

cbm::algo::kf::TrackParamBase::C41
T C41() const
Definition KfTrackParam.h:176

cbm::algo::kf::TrackParamBase::C10
T C10() const
Definition KfTrackParam.h:151

cbm::algo::kf::TrackParamBase::Y
T & Y()
Reference to y position [cm].
Definition KfTrackParam.h:410

cbm::algo::kf::TrackParamBase::C03
T & C03()
Definition KfTrackParam.h:468

cbm::algo::kf::TrackParamBase::C36
T & C36()
Definition KfTrackParam.h:495

cbm::algo::kf::TrackParamBase::C43
T & C43()
Definition KfTrackParam.h:500

cbm::algo::kf::TrackParamBase::GetNdfTime
T GetNdfTime() const
Gets NDF of time measurements.
Definition KfTrackParam.h:267

cbm::algo::kf::TrackParamBase::C51
T C51() const
Definition KfTrackParam.h:184

cbm::algo::kf::TrackParamBase::GetQp
T GetQp() const
Gets charge over momentum [ec/GeV].
Definition KfTrackParam.h:231

cbm::algo::kf::TrackParamBase::C56
T & C56()
Definition KfTrackParam.h:511

cbm::algo::kf::TrackParamBase::Y
T Y() const
Gets y position [cm].
Definition KfTrackParam.h:90

cbm::algo::kf::TrackParamBase::SetOneEntry
void SetOneEntry(const int ia, const TrackParamBase< T1 > &Tb)
Definition KfTrackParam.h:67

cbm::algo::kf::TrackParamBase::C53
T C53() const
Definition KfTrackParam.h:186

cbm::algo::kf::TrackParamBase::C22
T C22() const
Definition KfTrackParam.h:161

cbm::algo::kf::TrackParamBase::C52
T C52() const
Definition KfTrackParam.h:185

cbm::algo::kf::TrackParamBase::GetChiSqTime
T GetChiSqTime() const
Gets Chi-square of time measurements.
Definition KfTrackParam.h:264

cbm::algo::kf::TrackParamBase::C56
T C56() const
Definition KfTrackParam.h:189

cbm::algo::kf::TrackParamBase::Qp
T Qp() const
Gets charge over momentum [ec/GeV].
Definition KfTrackParam.h:99

cbm::algo::kf::TrackParamBase::Z
T Z() const
Gets z position [cm].
Definition KfTrackParam.h:84

cbm::algo::kf::TrackParamBase::InitVelocityRange
void InitVelocityRange(fscal minP)
Initializes inverse velocity range.
Definition KfTrackParam.h:765

cbm::algo::kf::TrackParamBase::GetPt
T GetPt() const
Gets transverse momentum.
Definition KfTrackParam.h:297

cbm::algo::kf::TrackParamBase::C60
T & C60()
Definition KfTrackParam.h:513

cbm::algo::kf::TrackParamBase::GetNdf
T GetNdf() const
Gets NDF of track fit model.
Definition KfTrackParam.h:261

cbm::algo::kf::TrackParamBase::C01
T & C01()
Definition KfTrackParam.h:466

cbm::algo::kf::TrackParamBase::GetQpError
T GetQpError() const
Gets error of charge over momentum [ec/GeV].
Definition KfTrackParam.h:234

cbm::algo::kf::TrackParamBase::Time
T Time() const
Gets time [ns].
Definition KfTrackParam.h:102

cbm::algo::kf::TrackParamBase::GetY
T GetY() const
Gets y position [cm].
Definition KfTrackParam.h:213

cbm::algo::kf::TrackParamBase::C34
T & C34()
Definition KfTrackParam.h:493

cbm::algo::kf::TrackParamBase::Vi
T & Vi()
Reference to inverse velocity [ns/cm].
Definition KfTrackParam.h:425

cbm::algo::kf::TrackParamBase::C46
T & C46()
Definition KfTrackParam.h:503

cbm::algo::kf::TrackParamBase::C23
T & C23()
Definition KfTrackParam.h:484

cbm::algo::kf::TrackParamBase::C
T C(int i, int j) const
Get covariance matrix element.
Definition KfTrackParam.h:123

cbm::algo::kf::TrackParamBase::C53
T & C53()
Definition KfTrackParam.h:508

cbm::algo::kf::TrackParamBase::SetX
void SetX(T v)
Sets x position [cm].
Definition KfTrackParam.h:311

cbm::algo::kf::TrackParamBase::SetNdfTime
void SetNdfTime(T v)
Sets NDF of time measurements.
Definition KfTrackParam.h:341

cbm::algo::kf::TrackParamBase::C45
T C45() const
Definition KfTrackParam.h:180

cbm::algo::kf::TrackParamBase::SetTy
void SetTy(T v)
Sets slope along y-axis.
Definition KfTrackParam.h:320

cbm::algo::kf::TrackParamBase::C04
T & C04()
Definition KfTrackParam.h:469

cbm::algo::kf::TrackParamBase::C25
T C25() const
Definition KfTrackParam.h:164

cbm::algo::kf::TrackParamBase::SetCovariance
void SetCovariance(T val)
Get covariance matrix element when indices are known at compile time.
Definition KfTrackParam.h:361

cbm::algo::kf::TrackParamBase::SetC50
void SetC50(T val)
Definition KfTrackParam.h:384

cbm::algo::kf::TrackParamBase::C66
T C66() const
Definition KfTrackParam.h:197

cbm::algo::kf::TrackParamBase::ChiSqTime
T & ChiSqTime()
Reference to Chi-square of time measurements.
Definition KfTrackParam.h:437

cbm::algo::kf::TrackParamBase::CovMatrix
CovMatrix_t & CovMatrix()
Reference to covariance matrix.
Definition KfTrackParam.h:428

cbm::algo::kf::TrackParamBase::C23
T C23() const
Definition KfTrackParam.h:162

cbm::algo::kf::TrackParamBase::C54
T C54() const
Definition KfTrackParam.h:187

cbm::algo::kf::TrackParamBase::C66
T & C66()
Definition KfTrackParam.h:519

cbm::algo::kf::TrackParamBase::SetC10
void SetC10(T val)
Definition KfTrackParam.h:370

cbm::algo::kf::TrackParamBase::CovMatrix_t
std::array< T, kNcovParam > CovMatrix_t
covariance matrix type
Definition KfTrackParam.h:41

cbm::algo::kf::TrackParamBase::C25
T & C25()
Definition KfTrackParam.h:486

cbm::algo::kf::TrackParamBase::C34
T C34() const
Definition KfTrackParam.h:171

cbm::algo::kf::TrackParamBase::C55
T C55() const
Definition KfTrackParam.h:188

cbm::algo::kf::TrackParamBase::C33
T & C33()
Definition KfTrackParam.h:492

cbm::algo::kf::TrackParamBase::SetC52
void SetC52(T val)
Definition KfTrackParam.h:386

cbm::algo::kf::TrackParamBase::TrackParamBase
TrackParamBase()=default

cbm::algo::kf::TrackParamBase::GetCovariance
T GetCovariance(int i, int j) const
Get covariance matrix element.
Definition KfTrackParam.h:255

cbm::algo::kf::TrackParamBase::C42
T C42() const
Definition KfTrackParam.h:177

cbm::algo::kf::TrackParamBase::C26
T C26() const
Definition KfTrackParam.h:165

cbm::algo::kf::TrackParamBase::C16
T & C16()
Definition KfTrackParam.h:479

cbm::algo::kf::TrackParamBase::C63
T & C63()
Definition KfTrackParam.h:516

cbm::algo::kf::TrackParamBase::C44
T & C44()
Definition KfTrackParam.h:501

cbm::algo::kf::TrackParamBase::C21
T & C21()
Definition KfTrackParam.h:482

cbm::algo::kf::TrackParamBase::C32
T & C32()
Definition KfTrackParam.h:491

cbm::algo::kf::TrackParamBase::Set
void Set(const TrackParamBase< T1 > &Tb)
Definition KfTrackParam.h:55

cbm::algo::kf::TrackParamBase::C61
T & C61()
Definition KfTrackParam.h:514

cbm::algo::kf::TrackParamBase::TrackParamBase
TrackParamBase(const TrackParamBase< T1 > &tr)
Definition KfTrackParam.h:46

cbm::algo::kf::TrackParamBase::GetChiSq
T GetChiSq() const
Gets Chi-square of track fit model.
Definition KfTrackParam.h:258

cbm::algo::kf::TrackParamBase::C00
T & C00()
Individual references to covariance matrix elements.
Definition KfTrackParam.h:465

cbm::algo::kf::TrackParamBase::C30
T C30() const
Definition KfTrackParam.h:167

cbm::algo::kf::TrackParamBase::Ty
T Ty() const
Gets slope along y-axis.
Definition KfTrackParam.h:96

cbm::algo::kf::TrackParamBase::C01
T C01() const
Definition KfTrackParam.h:144

cbm::algo::kf::TrackParamBase::C15
T & C15()
Definition KfTrackParam.h:478

cbm::algo::kf::TrackParamBase::SetC22
void SetC22(T val)
Definition KfTrackParam.h:374

cbm::algo::kf::TrackParamBase::C14
T & C14()
Definition KfTrackParam.h:477

cbm::algo::kf::TrackParamBase::C44
T C44() const
Definition KfTrackParam.h:179

cbm::algo::kf::TrackParamBase::SetC54
void SetC54(T val)
Definition KfTrackParam.h:388

cbm::algo::kf::TrackParamBase::SetQp
void SetQp(T v)
Sets charge over momentum [ec/GeV].
Definition KfTrackParam.h:323

cbm::algo::kf::TrackParamBase::X
T & X()
Reference to x position [cm].
Definition KfTrackParam.h:407

cbm::algo::kf::TrackParamBase::SetC40
void SetC40(T val)
Definition KfTrackParam.h:379

cbm::algo::kf::TrackParamBase::C50
T & C50()
Definition KfTrackParam.h:505

cbm::algo::kf::TrackParamBase::SetTime
void SetTime(T v)
Sets time [ns].
Definition KfTrackParam.h:326

cbm::algo::kf::TrackParamBase::C12
T C12() const
Definition KfTrackParam.h:153

cbm::algo::kf::TrackParamBase::GetX
T GetX() const
Gets x position [cm].
Definition KfTrackParam.h:207

cbm::algo::kf::TrackParam< fvec >
Definition KfTrackParam.h:653

cbm::algo::kf::TrackParam
TrackParam classes of different types.

cbm::algo::kf::fvec
Definition KfSimdPseudo.h:120

cbm::algo::kf::utils::VecCopy
Copies all/one SIMD entries from one class to the other class.
Definition KfUtils.h:179

cbm::algo::kf
Definition KfMeasurementTime.cxx:14

cbm::algo::kf::fscal
float fscal
Definition KfSimdPseudo.h:16