MapField3DRegular.h

Go to the documentation of this file.

#ifndef SIMG4COMMON_MAPFIELD3DREGULAR_H
#define SIMG4COMMON_MAPFIELD3DREGULAR_H
 
// Geant 4
#include "G4MagneticField.hh"
#include <vector>
 
namespace sim {
class MapField3DRegular : public G4MagneticField {
public:
  // Constructor
  explicit MapField3DRegular(const std::vector<double>& bX, const std::vector<double>& bY,
                             const std::vector<double>& bZ, const std::vector<double>& posX,
                             const std::vector<double>& posY, const std::vector<double>& posZ);
  // Destructor
  virtual ~MapField3DRegular() {}
 
  virtual void GetFieldValue(const G4double point[4], double* bField) const final;
 
private:
  std::vector<std::vector<std::vector<double>>> m_fieldX;
  std::vector<std::vector<std::vector<double>>> m_fieldY;
  std::vector<std::vector<std::vector<double>>> m_fieldZ;
  double m_minX, m_maxX, m_widthX;
  double m_minY, m_maxY, m_widthY;
  double m_minZ, m_maxZ, m_widthZ;
  size_t m_nX, m_nY, m_nZ;
};
} // namespace sim
#endif /* SIMG4COMMON_MAPFIELD3DREGULAR_H */