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 = nullptr) {
    return CylinderPointer(new Cylinder(radius, pos, rot, bounds));
  }
 
  static CylinderPointer build(Scalar radius,
                               const PositionType& pos,
                               const RotationType& rot,
                               Bounds* bounds = nullptr) {
    return CylinderPointer(new Cylinder(radius, pos, rot, bounds));
  }
 
  ~Cylinder() override {}
 
  // -- Extension of Surface interface for cylinder
 
  Scalar radius() const { return theRadius; }
 
  // -- Implementation of Surface interface
 
  using Surface::side;
  Side side(const LocalPoint& p, Scalar toler) const override;
 
  ConstReferenceCountingPointer<TangentPlane> tangentPlane(const GlobalPoint&) const override;
  ConstReferenceCountingPointer<TangentPlane> tangentPlane(const LocalPoint&) const override;
 
  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