CbmFieldMapSym1.cxx

Go to the documentation of this file.

/* Copyright (C) 2008-2020 Justus-Liebig-Universitaet Giessen, Giessen
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Elena Litvinenko, Elena Lebedeva [committer] */
 
// -------------------------------------------------------------------------
// -----                  CbmFieldMapSym1 source file                  -----
// -----                 Created 22/09/08  by E.Litvinenko             -----
// -------------------------------------------------------------------------
#include "CbmFieldMapSym1.h"
 
#include <TArrayF.h>  // for TArrayF
 
// -------------   Default constructor  ----------------------------------
CbmFieldMapSym1::CbmFieldMapSym1() : CbmFieldMap() { fType = 5; }
// ------------------------------------------------------------------------
 
 
// -------------   Standard constructor   ---------------------------------
CbmFieldMapSym1::CbmFieldMapSym1(const char* mapName, const char* fileType) : CbmFieldMap(mapName, fileType)
{
  fType = 5;
}
// ------------------------------------------------------------------------
 
 
// ------------   Constructor from CbmFieldPar   --------------------------
CbmFieldMapSym1::CbmFieldMapSym1(CbmFieldPar* fieldPar) : CbmFieldMap(fieldPar) { fType = 5; }
// ------------------------------------------------------------------------
 
 
// ------------   Destructor   --------------------------------------------
CbmFieldMapSym1::~CbmFieldMapSym1() {}
// ------------------------------------------------------------------------
 
 
// -----------   Get x component of the field   ---------------------------
Double_t CbmFieldMapSym1::GetBx(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBx, Local2Grid(local)) * local.hemiX;
}
// ------------------------------------------------------------------------
 
 
// -----------   Get y component of the field   ---------------------------
Double_t CbmFieldMapSym1::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 CbmFieldMapSym1::GetBz(Double_t x, Double_t y, Double_t z)
{
  auto local = Global2LocalReflected(x, y, z);
  return Interpolate(fBz, Local2Grid(local));
}
// ------------------------------------------------------------------------
 
 
// -----------   Get the field value at a point   ------------------------
void CbmFieldMapSym1::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;
  bField[1]  = Interpolate(fBy, grid);
  bField[2]  = Interpolate(fBz, grid);
}
// ------------------------------------------------------------------------