CoordinateSets.h

Go to the documentation of this file.

#ifndef GeometryVector_Geom_CoordinateSets_h
#define GeometryVector_Geom_CoordinateSets_h

#include <cmath>

namespace Geom {

    template <typename T>
    class Polar2Cartesian {
    public:
    Polar2Cartesian( const T& r, const T& phi) :
        r_(r), phi_(phi) {}

    const T& r() const   {return r_;}
    const T& phi() const {return phi_;}

    T x() const {return r_ * cos(phi_);}
    T y() const {return r_ * sin(phi_);}

    private:
    T r_;
    T phi_;
    };

    template <typename T>
    class Cylindrical2Cartesian {
    public:
    Cylindrical2Cartesian( const T& r, const T& phi, const T& z) :
        r_(r), phi_(phi), z_(z) {}

    const T& r() const   {return r_;}
    const T& phi() const {return phi_;}
    const T& z() const   {return z_;}

    T x() const {return r_ * cos(phi_);}
    T y() const {return r_ * sin(phi_);}

    private:
    T r_;
    T phi_;
    T z_;
    };

    template <typename T>
    class Spherical2Cartesian {
    public:
    Spherical2Cartesian( const T& theta, const T& phi, const T& mag) :
        theta_(theta), phi_(phi), r_(mag), 
        transv_( sin(theta)*mag) {}

    const T& theta() const   {return theta_;}
    const T& phi() const     {return phi_;}
    const T& r() const       {return r_;}

    T x() const {return transv_ * cos(phi());}
    T y() const {return transv_ * sin(phi());}
    T z() const {return cos(theta()) * r();}

    private:
    T theta_;
    T phi_;
    T r_;
    T transv_;
    };


    template <typename T>
    class Cartesian2Cartesian3D {
    public:
    Cartesian2Cartesian3D( const T& x, const T& y, const T& z) :
        x_(x), y_(y), z_(z) {}

    const T& x() const   {return x_;}
    const T& y() const   {return y_;}
    const T& z() const   {return z_;}
   private:
    T x_;
    T y_;
    T z_;
    };
}

#endif