CbmFieldMapSym3.cxx

Go to the documentation of this file.

/* Copyright (C) 2004-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Mohammad Al-Turany, Denis Bertini [committer] */
 
// -------------------------------------------------------------------------
// -----                    CbmFieldMapSym3 source file                -----
// -----         Created 12/01/04  by M. Al/Turany (CbmField.cxx)      -----
// -----                Redesign 13/02/06  by V. Friese                -----
// -------------------------------------------------------------------------
#include "CbmFieldMapSym3.h"
 
#include <TArrayF.h>  // for TArrayF
 
// -------------   Default constructor  ----------------------------------
CbmFieldMapSym3::CbmFieldMapSym3() : CbmFieldMap() { fType = 3; }
// ------------------------------------------------------------------------
 
 
// -------------   Standard constructor   ---------------------------------
CbmFieldMapSym3::CbmFieldMapSym3(const char* mapName, const char* fileType) : CbmFieldMap(mapName, fileType)
{
  fType = 3;
}
// ------------------------------------------------------------------------
 
 
// -------------   Constructor from CbmFieldPar   -------------------------
CbmFieldMapSym3::CbmFieldMapSym3(CbmFieldPar* fieldPar) : CbmFieldMap(fieldPar) { fType = 3; }
// ------------------------------------------------------------------------
 
 
// ------------   Destructor   --------------------------------------------
CbmFieldMapSym3::~CbmFieldMapSym3() {}
// ------------------------------------------------------------------------
 
 
// -----------   Get x component of the field   ---------------------------
Double_t CbmFieldMapSym3::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 CbmFieldMapSym3::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 CbmFieldMapSym3::GetBz(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBz, Local2Grid(local)) * local.hemiY * local.hemiZ;
}
// ------------------------------------------------------------------------
 
 
// -----------   Get the field value at a point   ------------------------
void CbmFieldMapSym3::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 * local.hemiZ;
}
// ------------------------------------------------------------------------