/data/refman/pasoursint/CMSSW_5_3_0/src/DataFormats/GeometryVector/interface/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