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;
  };
}
#endif /* SIMG4COMMON_MAPFIELD3DREGULAR_H */