/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/DataFormats/GeometrySurface/interface/Plane.h

Go to the documentation of this file.

#ifndef Geom_Plane_H
#define Geom_Plane_H

#include "DataFormats/GeometrySurface/interface/Surface.h"
#include "boost/intrusive_ptr.hpp" 

class Plane : public virtual Surface {
public:
  typedef ReferenceCountingPointer<Plane> PlanePointer;
  typedef ConstReferenceCountingPointer<Plane> ConstPlanePointer;

  static PlanePointer build(const PositionType& pos,
                                    const RotationType& rot,
                                    MediumProperties* mp=0) {
    return PlanePointer(new Plane(pos, rot, mp));
  }

  // -- DEPRECATED CONSTRUCTORS

  Plane( const PositionType& pos, const RotationType& rot) :
    Surface( pos, rot) {}

  Plane( const PositionType& pos, const RotationType& rot, MediumProperties* mp) : 
    Surface( pos, rot, mp) {}

// extension of Surface interface for planes

  GlobalVector normalVector() const {
    return GlobalVector( rotation().zx(), rotation().zy(), rotation().zz());
  }

  float localZ (const GlobalPoint& gp) const {
    return normalVector().dot(gp-position());
  }

  float localZ (const GlobalVector& gv) const {
    return normalVector().dot(gv);
  }

// implementation of Surface interface    

  virtual SurfaceOrientation::Side side( const LocalPoint& p, Scalar toler) const {
    return (std::abs(p.z())<toler) ? SurfaceOrientation::onSurface : 
        (p.z()>0 ? SurfaceOrientation::positiveSide : SurfaceOrientation::negativeSide);
  }

  virtual SurfaceOrientation::Side side( const GlobalPoint& p, Scalar toler) const {
    Scalar lz = localZ(p);
    return (std::abs(lz)<toler ? SurfaceOrientation::onSurface : 
            (lz>0 ? SurfaceOrientation::positiveSide : SurfaceOrientation::negativeSide));
  }

  virtual ReferenceCountingPointer<TangentPlane> tangentPlane (const GlobalPoint&) const;

  virtual ReferenceCountingPointer<TangentPlane> tangentPlane (const LocalPoint&) const;




};

#endif