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"
#include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"

class Cylinder final  : public Surface {
public:

  template<typename... Args>
    Cylinder(Scalar radius, Args&& ... args) :
    Surface(std::forward<Args>(args)...), theRadius(radius){}

  Cylinder(const PositionType& pos, const RotationType& rot, SimpleCylinderBounds const & bounds) :
    Surface(pos,rot, bounds.clone()),  theRadius(computeRadius(bounds)){}
  
  // average Rmin Rmax...
  static float computeRadius(Bounds const & bounds) {
    return  0.5f*(bounds.width() - bounds.thickness());
  }


  typedef ReferenceCountingPointer<Cylinder> CylinderPointer;
  typedef ConstReferenceCountingPointer<Cylinder> ConstCylinderPointer;
  typedef ReferenceCountingPointer<Cylinder> BoundCylinderPointer;
  typedef ConstReferenceCountingPointer<Cylinder> ConstBoundCylinderPointer;


  /*
  template<typename... Args>
  static CylinderPointer build(Args&& ... args) {
    return CylinderPointer(new Cylinder(std::forward<Args>(args)...));
  }
  */

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

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

  ~Cylinder(){}


  // -- 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 ConstReferenceCountingPointer<TangentPlane> tangentPlane (const GlobalPoint&) const;
  virtual ConstReferenceCountingPointer<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