KfField.cxx

Go to the documentation of this file.

/* Copyright (C) 2024 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergei Zharko [committer] */

 
#include "KfField.h"
 
#include <mutex>
#include <sstream>
#include <tuple>
 
using cbm::algo::kf::EFieldMode;
using cbm::algo::kf::Field;
using cbm::algo::kf::FieldBuilder;
using cbm::algo::kf::detail::ConcreteField;
 
// ---------------------------------------------------------------------------------------------------------------------
//
template<typename T>
std::string ConcreteField<T, EFieldMode::Interpolated>::ToString(int indentLevel, int verbose) const
{
  constexpr char IndentChar = '\t';
  std::stringstream msg;
  std::string indent(indentLevel, IndentChar);
  msg << indent << "Interpolated magnetic field";
  msg << "\n" << indent << "Field type: " << static_cast<int>(fFieldType);
  msg << '\n' << indent << "Field in the primary vertex: " << fPrimVertexField.ToString();
  msg << '\n' << indent << "Field slices:";
  for (int iSlice = 0; iSlice < GetNofFieldSlices(); ++iSlice) {
    const auto& fldSlice = fvFieldSlices[iSlice];
    msg << '\n' << indent << iSlice << ") " << fldSlice.ToString(0, verbose);
  }
  return msg.str();
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
template<typename T>
std::string ConcreteField<T, EFieldMode::Original>::ToString(int indentLevel, int /*verbose*/) const
{
  constexpr char IndentChar = '\t';
  std::stringstream msg;
  std::string indent(indentLevel, IndentChar);
  msg << indent << "Original field function";
  msg << "\n" << indent << "Field type: " << static_cast<int>(fFieldType);
  msg << "\n" << indent << "Field in the primary vertex: " << fPrimVertexField.ToString();
  msg << "\n" << indent << "Field slice z-positions: ";
  for (int iSlice = 0; iSlice < GetNofFieldSlices(); ++iSlice) {
    msg << "\n " << indent << iSlice << ") " << fvFieldSliceZ[iSlice];
  }
  return msg.str();
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
template<typename T>
std::string Field<T>::ToString(int indentLevel, int verbose) const
{
  constexpr char IndentChar = '\t';
  std::stringstream msg;
  std::string indent(indentLevel, IndentChar);
  kf::visit([&msg, indentLevel, verbose](const auto& v) { msg << v.ToString(indentLevel, verbose); }, fFieldVariant);
  return msg.str();
}
 
namespace cbm::algo::kf
{
  template class Field<float>;
  template class Field<double>;
  template class Field<fvec>;
}  // namespace cbm::algo::kf
 
 
// ---------------------------------------------------------------------------------------------------------------------
//
void FieldBuilder::AddSliceReference(double halfSizeX, double halfSizeY, double zRef)
{
  if (!fSliceReferences.emplace(halfSizeX, halfSizeY, zRef).second) {
    std::stringstream msg;
    msg << "FieldBuilder::AddReference: attempt of adding another slice reference with zRef = " << zRef
        << "(halfSizeX = " << halfSizeX << ", halfSizeY = " << halfSizeY << ").\nThe next slice references were "
        << "added:";
    for (const auto& el : fSliceReferences) {
      msg << "\n\t- halfSizeX = " << el.fHalfSizeX << ", halfSizeY = " << el.fHalfSizeY << ", zRef = " << el.fRefZ;
    }
    throw std::logic_error(msg.str());
  }
}