KfMeasurementU.h

Go to the documentation of this file.

/* Copyright (C) 2007-2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergey Gorbunov [committer], Igor Kulakov, Sergei Zharko */

 
#pragma once  // include this header only once per compilation unit
 
 
#include "KfDefs.h"
#include "KfUtils.h"
 
#include <boost/serialization/access.hpp>
 
#include <string>
 
namespace cbm::algo::kf
{
  template<typename DataT>
  class alignas(VcMemAlign) MeasurementU {
 
   public:
    friend class boost::serialization::access;

    MeasurementU() = default;

    MeasurementU(DataT cosPhi, DataT sinPhi, DataT u, DataT du2, DataT ndf)
      : fCosPhi(cosPhi)
      , fSinPhi(sinPhi)
      , fU(u)
      , fDu2(du2)
      , fNdf(ndf)
    {
    }

 
    void SetCosPhi(DataT cosPhi) { fCosPhi = cosPhi; }
    void SetSinPhi(DataT sinPhi) { fSinPhi = sinPhi; }
    void SetU(DataT u) { fU = u; }
    void SetDu2(DataT du2) { fDu2 = du2; }
    void SetNdf(DataT ndf) { fNdf = ndf; }
 
    DataT CosPhi() const { return fCosPhi; }
    DataT SinPhi() const { return fSinPhi; }
    DataT U() const { return fU; }
    DataT Du2() const { return fDu2; }
    DataT Ndf() const { return fNdf; }


    std::string ToString(int indentLevel = 0) const;

    bool IsFinite() const
    {
      return (utils::IsFinite(CosPhi()) && utils::IsFinite(SinPhi()) && utils::IsFinite(U()) && utils::IsFinite(Du2())
              && utils::IsFinite(Ndf()));
    }

    bool IsUndefined() const
    {
      return (utils::IsUndefined(CosPhi()) || utils::IsUndefined(SinPhi()) || utils::IsUndefined(U())
              || utils::IsUndefined(Du2()) || utils::IsUndefined(Ndf()));
    }

    friend class boost::serialization::access;
    template<class Archive>
    void serialize(Archive& ar, const unsigned int)
    {
      ar& fCosPhi;
      ar& fSinPhi;
      ar& fU;
      ar& fDu2;
      ar& fNdf;
    }
 
   private:
 
    DataT fCosPhi = defs::Undef<DataT>;  
    DataT fSinPhi = defs::Undef<DataT>;  
    DataT fU      = defs::Undef<DataT>;  
    DataT fDu2    = defs::Undef<DataT>;  

    DataT fNdf = defs::Undef<DataT>;
  } _fvecalignment;
 
}  // namespace cbm::algo::kf