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