Point3DBase.h

Go to the documentation of this file.

#ifndef GeometryVector_Point3DBase_h
#define GeometryVector_Point3DBase_h

#include "DataFormats/GeometryVector/interface/PointTag.h"
#include "DataFormats/GeometryVector/interface/PV3DBase.h"
#include "DataFormats/GeometryVector/interface/Point2DBase.h"
#include "DataFormats/GeometryVector/interface/Vector3DBase.h"

template <class T, class FrameTag>
class Point3DBase : public PV3DBase<T, PointTag, FrameTag> {
public:
  typedef PV3DBase<T, PointTag, FrameTag> BaseClass;
  typedef Vector3DBase<T, FrameTag> VectorType;
  typedef typename BaseClass::Cylindrical Cylindrical;
  typedef typename BaseClass::Spherical Spherical;
  typedef typename BaseClass::Polar Polar;
  typedef typename BaseClass::BasicVectorType BasicVectorType;

  Point3DBase() {}

  template <class U>
  Point3DBase(const Point3DBase<U, FrameTag>& p) : BaseClass(p.basicVector()) {}

  Point3DBase(const T& x, const T& y, const T& z) : BaseClass(x, y, z) {}

  explicit Point3DBase(const Cylindrical& set) : BaseClass(set) {}

  explicit Point3DBase(const Polar& set) : BaseClass(set) {}

  Point3DBase(const Geom::Theta<T>& th, const Geom::Phi<T>& ph, const T& r) : BaseClass(th, ph, r) {}

  Point3DBase(const T& x, const T& y) : BaseClass(x, y, 0) {}
  explicit Point3DBase(Point2DBase<T, FrameTag> p) : BaseClass(p.x(), p.y(), 0) {}

  template <class U>
  explicit Point3DBase(const Basic3DVector<U>& v) : BaseClass(v) {}

  // equality
  bool operator==(const Point3DBase& rh) const { return this->basicVector() == rh.basicVector(); }

  template <class U>
  Point3DBase& operator+=(const Vector3DBase<U, FrameTag>& v) {
    this->theVector += v.basicVector();
    return *this;
  }

  template <class U>
  Point3DBase& operator-=(const Vector3DBase<U, FrameTag>& v) {
    this->theVector -= v.basicVector();
    return *this;
  }
};

template <typename T, typename U, class Frame>
inline Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> operator+(const Point3DBase<T, Frame>& p,
                                                                           const Vector3DBase<U, Frame>& v) {
  typedef Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> RT;
  return RT(p.basicVector() + v.basicVector());
}

template <typename T, typename U, class Frame>
inline Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> operator+(const Vector3DBase<T, Frame>& p,
                                                                           const Point3DBase<U, Frame>& v) {
  typedef Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> RT;
  return RT(p.basicVector() + v.basicVector());
}

template <typename T, typename U, class Frame>
inline Vector3DBase<typename PreciseFloatType<T, U>::Type, Frame> operator-(const Point3DBase<T, Frame>& p1,
                                                                            const Point3DBase<U, Frame>& p2) {
  typedef Vector3DBase<typename PreciseFloatType<T, U>::Type, Frame> RT;
  return RT(p1.basicVector() - p2.basicVector());
}

template <typename T, typename U, class Frame>
inline Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> operator-(const Point3DBase<T, Frame>& p,
                                                                           const Vector3DBase<U, Frame>& v) {
  typedef Point3DBase<typename PreciseFloatType<T, U>::Type, Frame> RT;
  return RT(p.basicVector() - v.basicVector());
}
#endif  // GeometryVector_Point3DBase_h