CbmBsField.h

Go to the documentation of this file.

/* Copyright (C) 2006-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Elena Litvinenko, Denis Bertini [committer], Florian Uhlig */
 
// -------------------------------------------------------------------------
// -----                      CbmBsField  header file                  -----
// -----                Created 12/05/06  by E.I.Litvinenko            -----
// -------------------------------------------------------------------------
 

 
#ifndef CBMBSF_H
#define CBMBSF_H
 
#include "CbmFieldMap.h"  // for CbmFieldMap
 
#include <Rtypes.h>      // for THashConsistencyHolder, ClassDef
#include <RtypesCore.h>  // for Double_t, Int_t, Float_t
#include <TString.h>     // for TString
 
#include <cstring>  // for strcpy
 
class CbmFieldPar;
class TArrayF;
 
class CbmBsField : public CbmFieldMap {
 public:
  CbmBsField(const char* pBsName, Int_t symType = 3);
  CbmBsField();
  CbmBsField(CbmFieldPar* fieldPar);
 
  virtual ~CbmBsField();

  virtual void Init() override;
 
  void writeBsRootfile(const char* name);  // Write Field Splined  to root file
  void readBsRootfile(const char* name);   // Read  Field Splined  from root file
  void readFortranAsciifiles(const char* MapFile1, const char* MapFile2,
                             const char* MapFile3);  // for the first use only
 
  void CalculateMapFromBs(Int_t pNx = 0, Int_t pNy = 0,
                          Int_t pNz = 0);  // calculate map pNx,pNy,Pnz dim (default - LL1-3,LL2-3,LL3-3)
 
 
  TArrayF* GetBsBx() { return fBsBx; }
  TArrayF* GetBsBy() { return fBsBy; }
  TArrayF* GetBsBz() { return fBsBz; }
  TArrayF* GetX() { return fX; }
  TArrayF* GetY() { return fY; }
  TArrayF* GetZ() { return fZ; }
 
  void GetBsName(char* pBsName) { strcpy(pBsName, fBsName.Data()); }
  void SetBsName(const char* pBsName)
  {
    fBsName = pBsName;
    this->SetName(pBsName);
  }

  virtual Double_t GetBx(Double_t x, Double_t y, Double_t z) override;
  virtual Double_t GetBy(Double_t x, Double_t y, Double_t z) override;
  virtual Double_t GetBz(Double_t x, Double_t y, Double_t z) override;

  virtual void GetFieldValue(const Double_t Point[3], Double_t* Bfield) override;
 
 
  // private:
 protected:
  Int_t NDIM, LL1, LL2, LL3, II1, II2, II3;
 
  TArrayF* fBsBx;  // NDIM coefficients for Bx
  TArrayF* fBsBy;  // NDIM coefficients for By
  TArrayF* fBsBz;  // NDIM coefficients for Bz
  TArrayF* fX;     // LL1 x coordinates of grid points
  TArrayF* fY;     // LL2 y coordinates of grid points
  TArrayF* fZ;     // LL3 z coordinates of grid points
 
  Float_t* UX1;  
  Float_t* UX2;  
  Float_t* UX3;  
  Float_t* F0;   
  Float_t* G0;   
  Float_t* U0;   
 
  TString fBsName;  // basic field name used in Bs root file
 protected:
  void PALC0(Double_t X, Double_t Y, Double_t Z, Double_t* BX, Double_t* BY,
             Double_t* BZ);  // from Fortran code
  Float_t SPL0(Double_t T, Double_t X0, Double_t X1, Double_t X2, Double_t X3,
               Double_t X4);  // from Fortran code
 
 private:
  CbmBsField(const CbmBsField&);
  CbmBsField& operator=(const CbmBsField&);
 
  ClassDefOverride(CbmBsField, 2)
};
 
#endif