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