KfToolsField.h

Go to the documentation of this file.

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

 
#pragma once  // include this header only once per compilation unit
 
#include "FairField.h"
#include "FairRunAna.h"
#include "KfDefs.h"
#include "KfFieldRegion.h"
 
namespace cbm::kf::tools
{
  inline void SetOriginalCbmField()
  {
    static FairField* fairField{nullptr};
 
    auto* run = FairRunAna::Instance();
 
    if (!run) {
      LOG(fatal) << "FairRunAna::Instance() is not initialized. Cannot set the original field function.";
      return;
    }
 
    fairField = run->GetField();
 
    if (nullptr == fairField) {
      LOG(debug) << "Kf Setup: No Cbm Field found";
      cbm::algo::kf::GlobalField::SetFieldFunction(cbm::algo::kf::EFieldType::Null, cbm::algo::kf::defs::ZeroFieldFn);
    }
    else {
 
      LOG(debug) << "Kf Setup: Cbm Field " << fairField->ClassName() << " " << fairField->GetName();
      LOG(debug) << "Field parameters at the origin: " << fairField->GetBx(0., 0., 0.) << ", "
                 << fairField->GetBy(0., 0., 0.) << ", " << fairField->GetBz(0., 0., 0.) << ", type "
                 << fairField->GetType();
 
      static auto fld = [&](double x, double y, double z) {
        double pos[3] = {x, y, z};
        double B[3]   = {0., 0., 0.};
        fairField->GetFieldValue(pos, B);
        return std::tuple(B[0], B[1], B[2]);
      };
      cbm::algo::kf::GlobalField::SetFieldFunction(cbm::algo::kf::EFieldType::Normal, fld);
    }
  }
 
}  // namespace cbm::kf::tools