PV3DBase.h

Go to the documentation of this file.

#ifndef GeometryVector_PV3DBase_h
#define GeometryVector_PV3DBase_h

#include "DataFormats/GeometryVector/interface/Basic3DVector.h"
#include <iosfwd>

template <class T, class PVType, class FrameType>
class PV3DBase {
public:

  typedef T                                         ScalarType;
  typedef Basic3DVector<T>                          BasicVectorType;
  typedef typename BasicVectorType::Cylindrical     Cylindrical;
  typedef typename BasicVectorType::Spherical       Spherical;
  typedef typename BasicVectorType::Polar           Polar;
  typedef typename BasicVectorType::MathVector  MathVector;


  PV3DBase() : theVector() {}

  PV3DBase(const T & x, const T & y, const T & z) : theVector(x, y, z) {}

  PV3DBase( const Cylindrical& set) : theVector( set) {}

  PV3DBase( const Polar& set) : theVector( set) {}

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

  template <class U>
  explicit PV3DBase( const Basic3DVector<U>& v) : theVector(v) {}

  const BasicVectorType& basicVector() const { return theVector;}
#ifndef __REFLEX__
  MathVector const & mathVector() const { return theVector.v;}
  MathVector & mathVector() { return theVector.v;}
#endif  

  T x() const { return basicVector().x();}
  T y() const { return basicVector().y();}
  T z() const { return basicVector().z();}

  T mag2() const { return basicVector().mag2();}
  T mag() const  { return basicVector().mag();}
  T barePhi() const  { return basicVector().barePhi();}
  Geom::Phi<T> phi() const  { return basicVector().phi();}

  T perp2() const { return basicVector().perp2();}
  T perp() const  { return basicVector().perp();}
  T transverse() const  { return basicVector().transverse();}
  T bareTheta() const { return basicVector().bareTheta();}
  Geom::Theta<T> theta() const { return basicVector().theta();}
  T eta() const   { return basicVector().eta();}

protected:
  BasicVectorType theVector;
};

template <class T, class PV, class F>
inline std::ostream & operator<<( std::ostream& s, const PV3DBase<T,PV,F>& v) {
  return s << v.basicVector();
} 

#endif // GeometryVector_PV3DBase_h