CbmFieldMapSym2.cxx

Go to the documentation of this file.

/* Copyright (C) 2005-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Elena Litvinenko, Denis Bertini [committer] */
 
// -------------------------------------------------------------------------
// -----                  CbmFieldMapSym2 source file               -----
// -----      Created 09/06/05  by E. Litvinenko (CbmFieldMapLIT.cxx)  -----
// -----                Redesign 20/02/06  by V. Friese                -----
// -------------------------------------------------------------------------
#include "CbmFieldMapSym2.h"
 
#include <TArrayF.h>  // for TArrayF
 
// -------------   Default constructor  ----------------------------------
CbmFieldMapSym2::CbmFieldMapSym2() : CbmFieldMap() { fType = 2; }
// ------------------------------------------------------------------------
 
 
// -------------   Standard constructor   ---------------------------------
CbmFieldMapSym2::CbmFieldMapSym2(const char* mapName, const char* fileType) : CbmFieldMap(mapName, fileType)
{
  fType = 2;
}
// ------------------------------------------------------------------------
 
 
// ------------   Constructor from CbmFieldPar   --------------------------
CbmFieldMapSym2::CbmFieldMapSym2(CbmFieldPar* fieldPar) : CbmFieldMap(fieldPar) { fType = 2; }
// ------------------------------------------------------------------------
 
 
// ------------   Destructor   --------------------------------------------
CbmFieldMapSym2::~CbmFieldMapSym2() {}
// ------------------------------------------------------------------------
 
 
// -----------   Get x component of the field   ---------------------------
Double_t CbmFieldMapSym2::GetBx(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBx, Local2Grid(local)) * local.hemiX * local.hemiY;
}
// ------------------------------------------------------------------------
 
 
// -----------   Get y component of the field   ---------------------------
Double_t CbmFieldMapSym2::GetBy(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBy, Local2Grid(local));
}
// ------------------------------------------------------------------------
 
 
// -----------   Get z component of the field   ---------------------------
Double_t CbmFieldMapSym2::GetBz(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBz, Local2Grid(local)) * local.hemiY;
}
// ------------------------------------------------------------------------
 
 
// -----------   Get the field value at a point   ------------------------
void CbmFieldMapSym2::GetFieldValue(const Double_t point[3], Double_t* bField)
{
  auto local = Global2LocalReflected(point[0], point[1], point[2]);
  auto grid  = Local2Grid(local);
  bField[0]  = Interpolate(fBx, grid) * local.hemiX * local.hemiY;
  bField[1]  = Interpolate(fBy, grid);
  bField[2]  = Interpolate(fBz, grid) * local.hemiY;
}
// ------------------------------------------------------------------------