CbmFieldMap.h

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], Volker Friese */
 
// -------------------------------------------------------------------------
// -----                      CbmFieldMap header file                  -----
// -----          Created 12/01/04  by M. Al/Turany (CbmField.h)       -----
// -----                Redesign 13/02/06  by V. Friese                -----
// -------------------------------------------------------------------------
 
 
#ifndef CBMFIELDMAP_H
#define CBMFIELDMAP_H 1
 
 
#include <FairField.h>  // for FairField
 
#include <Rtypes.h>      // for THashConsistencyHolder, ClassDef
#include <RtypesCore.h>  // for Double_t, Int_t, Bool_t, Option_t
#include <TString.h>     // for TString
 
class CbmFieldMapCreator;
class CbmFieldMapData;
class CbmFieldPar;
class TArrayF;
 
class CbmFieldMap : public FairField {
 
 
public:
  CbmFieldMap();
 
 
  CbmFieldMap(const char* mapName, const char* fileType = "R");
 
 
  CbmFieldMap(CbmFieldPar* fieldPar);
 
 
  CbmFieldMap(CbmFieldMapCreator* creator);
 
 
  virtual ~CbmFieldMap();
 
 
  virtual void Init();
 
 
  virtual void Init(Int_t nX, Double_t xMin, Double_t xMax, Int_t nY, Double_t yMin, Double_t yMax, Int_t nZ,
                    Double_t zMin, Double_t zMax, TArrayF* bx, TArrayF* by, TArrayF* bz);
 
 
  virtual Double_t GetBx(Double_t x, Double_t y, Double_t z);
  virtual Double_t GetBy(Double_t x, Double_t y, Double_t z);
  virtual Double_t GetBz(Double_t x, Double_t y, Double_t z);
 
 
  virtual Bool_t IsInside(Double_t x, Double_t y, Double_t z, Int_t& ix, Int_t& iy, Int_t& iz, Double_t& dx,
                          Double_t& dy, Double_t& dz);
 
 
  void WriteAsciiFile(const char* fileName);
 
 
  void WriteRootFile(const char* fileName, const char* mapName);
 
 
  virtual void SetPosition(Double_t x, Double_t y, Double_t z);
 
 
  virtual void SetScale(Double_t factor) { fScale = factor; }
 
 
  Double_t GetXmin() const { return fXmin; }
  Double_t GetYmin() const { return fYmin; }
  Double_t GetZmin() const { return fZmin; }
  Double_t GetXmax() const { return fXmax; }
  Double_t GetYmax() const { return fYmax; }
  Double_t GetZmax() const { return fZmax; }
  Double_t GetXstep() const { return fXstep; }
  Double_t GetYstep() const { return fYstep; }
  Double_t GetZstep() const { return fZstep; }
  Int_t GetNx() const { return fNx; }
  Int_t GetNy() const { return fNy; }
  Int_t GetNz() const { return fNz; }
 
 
  Double_t GetPositionX() const { return fPosX; }
  Double_t GetPositionY() const { return fPosY; }
  Double_t GetPositionZ() const { return fPosZ; }
 
 
  Double_t GetScale() const { return fScale; }
 
 
  TArrayF* GetBx() const { return fBx; }
  TArrayF* GetBy() const { return fBy; }
  TArrayF* GetBz() const { return fBz; }
 
 
  const char* GetFileName() { return fFileName.Data(); }
 
 
  virtual void Print(Option_t* = "") const;
 
 
protected:
  void Reset();
 
 
  void ReadAsciiFile(const char* fileName);
 
  void ReadAsciiFile2018(const char* fileName);
 
  void ReadRootFile(const char* fileName, const char* mapName);
 
 
  void SetField(const CbmFieldMapData* data);
 
 
  Double_t Interpolate(Double_t dx, Double_t dy, Double_t dz);
 
 
  TString fFileName;
 
 
  Double_t fScale;
 
 
  Double_t fPosX, fPosY, fPosZ;
 
 
  Double_t fXmin, fXmax, fXstep;
  Double_t fYmin, fYmax, fYstep;
  Double_t fZmin, fZmax, fZstep;
 
 
  Int_t fNx, fNy, fNz;
 
 
  TArrayF* fBx;
  TArrayF* fBy;
  TArrayF* fBz;
 
 
  Double_t fHa[2][2][2];  
  Double_t fHb[2][2];     
  Double_t fHc[2];        
 
  Double_t fBxOrigin;  
  Double_t fByOrigin;  
  Double_t fBzOrigin;  
 
private:
  CbmFieldMap(const CbmFieldMap&);
  CbmFieldMap& operator=(const CbmFieldMap&);
 
  ClassDef(CbmFieldMap, 2)
};
 
 
#endif