MapField2DRegular.h

Go to the documentation of this file.

#ifndef SIMG4COMMON_MAPFIELD2DREGULAR_H
#define SIMG4COMMON_MAPFIELD2DREGULAR_H
 
// Geant 4
#include "G4MagneticField.hh"
#include <vector>
 
namespace sim {
  class MapField2DRegular : public G4MagneticField {
    public:
    // Constructor
    explicit MapField2DRegular(const std::vector<double>& bR,
                               const std::vector<double>& bZ,
                               const std::vector<double>& posR,
                               const std::vector<double>& posZ);
    // Destructor
    virtual ~MapField2DRegular() {}
 
    virtual void GetFieldValue(const G4double point[4], double* bField) const final;
 
  private:
    std::vector<std::vector<double>> m_fieldR;
    std::vector<std::vector<double>> m_fieldZ;
    double m_minR, m_maxR, m_widthR;
    double m_minZ, m_maxZ, m_widthZ;
    size_t m_nR, m_nZ;
  };
}
 
#endif /* SIMG4COMMON_MAPFIELD2DREGULAR_H */