Cylindrical3DMagneticField.h

#ifndef CYLINDRICAL3DMAGNETICFIELD_H
#define CYLINDRICAL3DMAGNETICFIELD_H 1

// stl
#include <vector>

// CLHEP
#include <CLHEP/Matrix/Vector.h>
#include <CLHEP/Vector/ThreeVector.h>

// MarlinTPC
#include "MagneticField.h"

namespace EVENT {

        class LCCollection;
}

namespace marlintpc {

class Cylindrical3DMagneticField: public MagneticField {

public:
        Cylindrical3DMagneticField();

        virtual ~Cylindrical3DMagneticField();

        virtual Field *new_field( void );

        virtual bool construct( EVENT::LCCollection * );

        struct FieldMapEntry {

                float _Br;
                float _Bphi;
                float _Bz;
        };

protected:

        CLHEP::Hep3Vector get_field_value_coil_coordinates( const CLHEP::Hep3Vector & );


private:

        double _r_min;
        double _r_max;
        double _phi_min;
        double _phi_max;
        double _z_min;
        double _z_max;

        int _r_bins;
        int _phi_bins;
        int _z_bins;

        double _dr;
        double _dphi;
        double _dz;

        std::vector<std::vector<std::vector<FieldMapEntry> > > _field_map;
};

} // namespace marlintpc

#endif // CYLINDRICAL3DMAGNETICFIELD_H