MagneticFieldMap3D.h

Go to the documentation of this file.

 
#ifndef SIMCORE_MAGNETICFIELDMAP3D_H_
#define SIMCORE_MAGNETICFIELDMAP3D_H_
 
// Geant4
#include "Framework/Logger.h"
#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>>> x_field_;
  vector<vector<vector<double>>> y_field_;
  vector<vector<vector<double>>> z_field_;
 
  /*
   * 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 x_offset_;
  double y_offset_;
  double z_offset_;
 
  /*
   * Flags for inverting dimensions.
   */
  bool invert_x_, invert_y_, invert_z_;
 
  enableLogging("MagneticFieldMap3D")
};
 
}  // namespace simcore
// namespace
 
#endif