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_;
  };
}  // namespace Geom
 
#endif