Cylinder.h

Go to the documentation of this file.

#ifndef Geom_Cylinder_H
#define Geom_Cylinder_H

#include "DataFormats/GeometrySurface/interface/Surface.h"
#include "DataFormats/GeometrySurface/interface/Plane.h"

class Cylinder : public virtual Surface {
public:
  typedef ReferenceCountingPointer<Cylinder> CylinderPointer;
  typedef ConstReferenceCountingPointer<Cylinder> ConstCylinderPointer;


  static CylinderPointer build(const PositionType& pos,
                   const RotationType& rot,
                   Scalar radius,
                   MediumProperties* mp=0) {
    return CylinderPointer(new Cylinder(pos, rot, radius, mp));
  }


  // -- DEPRECATED CONSTRUCTORS

  Cylinder( const PositionType& pos, const RotationType& rot, Scalar radius) :
    Surface( pos, rot), theRadius(radius) {}

  Cylinder( const PositionType& pos, const RotationType& rot, Scalar radius,
     MediumProperties* mp) : 
    Surface( pos, rot, mp), theRadius(radius) {}

  // -- Extension of Surface interface for cylinder

  Scalar radius() const {return theRadius;}

  // -- Implementation of Surface interface    

  using Surface::side;
  virtual Side side( const LocalPoint& p, Scalar toler) const;

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

  Plane fastTangent(const GlobalPoint& aPoint) const{
    GlobalVector yPlane(rotation().z());
    GlobalVector xPlane(yPlane.cross(aPoint-position()));
    return Plane(aPoint,RotationType(xPlane, yPlane));
  }

  Plane fastTangent(const LocalPoint& aPoint) const {
    return fastTangent(toGlobal(aPoint));
  } 

private:

  Scalar theRadius;

};

#endif