MagneticField.h

#ifndef MAGNETICFIELD_H
#define MAGNETICFIELD_H 1



// MarlinTPC
#include "Field.h"

// CLHEP
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/EulerAngles.h"


namespace marlintpc {

class MagneticField: public Field {

public:
        MagneticField( std::string field_typename );

        virtual ~MagneticField();

        virtual Field *new_field( void ) = 0;

        virtual bool construct( EVENT::LCCollection * ) = 0;

        virtual CLHEP::Hep3Vector get_field_value( const CLHEP::Hep3Vector & );

        virtual CLHEP::Hep3Vector get_field_error( );

        void set_translation ( const CLHEP::Hep3Vector & );

        void set_rotation ( const double &, const double & );

        void set_scale_factor ( const double & );


protected:

        virtual CLHEP::Hep3Vector get_field_value_coil_coordinates( const CLHEP::Hep3Vector & ) = 0;

private:

        CLHEP::Hep3Vector _translation;
        double _x_rotation;
        double _y_rotation;
        double _scale_factor;
};

} // namespace marlintpc

#endif // MAGNETICFIELD_H