MagneticFieldMap3D.h

Go to the documentation of this file.

 
#ifndef SIMCORE_MAGNETICFIELDMAP3D_H_
#define SIMCORE_MAGNETICFIELDMAP3D_H_
 
// Geant4
#include "G4MagneticField.hh"
 
// STL
#include <vector>
using std::vector;
 
namespace simcore {
 
class MagneticFieldMap3D : public G4MagneticField {
 public:
  MagneticFieldMap3D(const char* filename, double xOffset, double yOffset,
                     double zOffset);
 
  void GetFieldValue(const double point[4], double* bfield) const;
 
 private:
  /*
   * Storage space for the table.
   */
  vector<vector<vector<double>>> xField_;
  vector<vector<vector<double>>> yField_;
  vector<vector<vector<double>>> zField_;
 
  /*
   * The dimensions of the table.
   */
  int nx_, ny_, nz_;
 
  /*
   * The physical limits of the defined region.
   */
  double minx_, maxx_, miny_, maxy_, minz_, maxz_;
 
  /*
   * The physical extent of the defined region.
   */
  double dx_, dy_, dz_;
 
  /*
   * Offsets if field map is not in global coordinates
   */
  double xOffset_;
  double yOffset_;
  double zOffset_;
 
  /*
   * Flags for inverting dimensions.
   */
  bool invertX_, invertY_, invertZ_;
};
 
}  // namespace simcore
// namespace
 
#endif