CLHEP 2.0.4.7 Reference Documentation     CLHEP Home Page     CLHEP Documentation     CLHEP Bug Reports

Point3D.h

Go to the documentation of this file.

// -*- C++ -*-
// $Id: Point3D.h,v 1.3.4.1 2004/07/23 23:09:28 garren Exp $
// ---------------------------------------------------------------------------
//
// This file is a part of the CLHEP - a Class Library for High Energy Physics.
//
// History:
// 09.09.96 E.Chernyaev - initial version
// 12.06.01 E.Chernyaev - CLHEP-1.7: introduction of BasicVector3D to decouple
//                        the functionality from CLHEP::Hep3Vector
// 01.04.03 E.Chernyaev - CLHEP-1.9: template version
//

#ifndef HEP_POINT3D_H
#define HEP_POINT3D_H

#include <iosfwd>
#include "CLHEP/Geometry/defs.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Geometry/BasicVector3D.h"

namespace HepGeom {

  class Transform3D;

  template<class T>
  class Point3D : public BasicVector3D<T> {};

  template<>
  class Point3D<float> : public BasicVector3D<float> {
  public:
    Point3D() {}

    Point3D(float x, float y, float z) : BasicVector3D<float>(x,y,z) {}

    explicit Point3D(const float * a)
      : BasicVector3D<float>(a[0],a[1],a[2]) {}

    Point3D(const Point3D<float> & v) : BasicVector3D<float>(v) {}

    Point3D(const BasicVector3D<float> & v) : BasicVector3D<float>(v) {}

    ~Point3D() {}

    Point3D<float> & operator=(const Point3D<float> & v) {
      set(v.x(),v.y(),v.z()); return *this;
    }

    Point3D<float> & operator=(const BasicVector3D<float> & v) {
      set(v.x(),v.y(),v.z()); return *this;
    }

    float distance2() const { return mag2(); }

    float distance2(const Point3D<float> & p) const {
      float dx = p.x()-x(), dy = p.y()-y(), dz = p.z()-z();
      return dx*dx + dy*dy + dz*dz;
    }

    float distance() const { return sqrt(distance2()); }

    float distance(const Point3D<float> & p) const {
      return sqrt(distance2(p));
    }

    Point3D<float> & transform(const Transform3D & m);
  };

  Point3D<float>
  operator*(const Transform3D & m, const Point3D<float> & p);

  template<>
  class Point3D<double> : public BasicVector3D<double> {
  public:
    Point3D() {}

    Point3D(double x, double y, double z) : BasicVector3D<double>(x,y,z) {}

    explicit Point3D(const float * a)
      : BasicVector3D<double>(a[0],a[1],a[2]) {}

    explicit Point3D(const double * a)
      : BasicVector3D<double>(a[0],a[1],a[2]) {}

    Point3D(const Point3D<double> & v) : BasicVector3D<double>(v) {}

    Point3D(const BasicVector3D<float> & v) : BasicVector3D<double>(v) {}

    Point3D(const BasicVector3D<double> & v) : BasicVector3D<double>(v) {}

    ~Point3D() {}

    Point3D(const CLHEP::Hep3Vector & v)
      : BasicVector3D<double>(v.x(),v.y(),v.z()) {}

    operator CLHEP::Hep3Vector () const { return CLHEP::Hep3Vector(x(),y(),z()); }

    Point3D<double> & operator=(const Point3D<double> & v) {
      set(v.x(),v.y(),v.z()); return *this;
    }

    Point3D<double> & operator=(const BasicVector3D<float> & v) {
      set(v.x(),v.y(),v.z()); return *this;
    }

    Point3D<double> & operator=(const BasicVector3D<double> & v) {
      set(v.x(),v.y(),v.z()); return *this;
    }

    double distance2() const { return mag2(); }

    double distance2(const Point3D<double> & p) const {
      double dx = p.x()-x(), dy = p.y()-y(), dz = p.z()-z();
      return dx*dx + dy*dy + dz*dz;
    }

    double distance() const { return sqrt(distance2()); }

    double distance(const Point3D<double> & p) const {
      return sqrt(distance2(p));
    }

    Point3D<double> & transform(const Transform3D & m);
  };

  Point3D<double>
  operator*(const Transform3D & m, const Point3D<double> & p);

} /* namespace HepGeom */


#ifdef ENABLE_BACKWARDS_COMPATIBILITY
//  backwards compatibility will be enabled ONLY in CLHEP 1.9
#include "CLHEP/config/CLHEP.h"
#include "CLHEP/Geometry/Normal3D.h"
#include "CLHEP/Geometry/Transform3D.h"
typedef HepGeom::Point3D<double> HepPoint3D;
#endif

#endif /* HEP_POINT3D_H */

---