FM2DMagnetoStatic.h

Go to the documentation of this file.

#ifndef CLASSIC_FIELDMAP2DMAGNETOSTATIC_HH
#define CLASSIC_FIELDMAP2DMAGNETOSTATIC_HH
 
#include "Fields/Fieldmap.h"
 
class FM2DMagnetoStatic: public Fieldmap {
 
public:
    virtual bool getFieldstrength(const Vector_t<double, 3> &R, Vector_t<double, 3> &E, Vector_t<double, 3> &B) const;
    virtual bool getFieldDerivative(const Vector_t<double, 3> &R, Vector_t<double, 3> &E, Vector_t<double, 3> &B, const DiffDirection &dir) const;
    virtual void getFieldDimensions(double &zBegin, double &zEnd) const;
    virtual void getFieldDimensions(double &xIni, double &xFinal, double &yIni, double &yFinal, double &zIni, double &zFinal) const;
    virtual void swap();
    virtual void getInfo(Inform *msg);
    virtual double getFrequency() const;
    virtual void setFrequency(double freq);
 
    virtual bool isInside(const Vector_t<double, 3> &r) const;
private:
    FM2DMagnetoStatic(std::string aFilename);
    ~FM2DMagnetoStatic();
 
    virtual void readMap();
    virtual void freeMap();
 
    double *FieldstrengthBz_m;    
    double *FieldstrengthBr_m;    
    double rbegin_m;
    double rend_m;
    double zbegin_m;
    double zend_m;
    double hz_m;                   
    double hr_m;                   
    int num_gridpr_m;              
    int num_gridpz_m;              
    bool swap_m;
    friend class Fieldmap;
};
 
inline bool FM2DMagnetoStatic::isInside(const Vector_t<double, 3> &r) const
{
    return r(2) >= zbegin_m && r(2) < zend_m && std::sqrt(r(0)*r(0) + r(1)*r(1)) < rend_m;
}
 
#endif