Cylinder.cc

Go to the documentation of this file.

#include "DataFormats/GeometrySurface/interface/Cylinder.h"
#include "DataFormats/GeometrySurface/interface/TangentPlane.h"
#include "DataFormats/GeometrySurface/interface/GeomExceptions.h"

#include <cfloat>

Surface::Side Cylinder::side( const LocalPoint& p, Scalar toler) const
{
  Scalar lz = p.perp() - radius();
  return (fabs(lz)<toler ? SurfaceOrientation::onSurface : 
      (lz>0 ? SurfaceOrientation::positiveSide : SurfaceOrientation::negativeSide));
}

ReferenceCountingPointer<TangentPlane> 
Cylinder::tangentPlane (const LocalPoint& aPoint) const
{
  return tangentPlane(toGlobal(aPoint));
}

ReferenceCountingPointer<TangentPlane> 
Cylinder::tangentPlane (const GlobalPoint& aPoint) const
{
  //
  // Tangent plane at specified point. In order to avoid
  // possible numerical problems currently no attempt is made 
  // to verify, if the point is actually on the cylinder.
  //
  // local y parallel to axis
  GlobalVector yPlane(rotation().z());
  // local x normal to y and a vector linking the specified
  // point with the axis
  GlobalVector xPlane(yPlane.cross(aPoint-position()));
  Scalar size = std::max(std::max(std::abs(xPlane.x()),std::abs(xPlane.y())), std::abs(xPlane.z()));
  if ( size<FLT_MIN ) 
    throw GeometryError("Attempt to construct TangentPlane on cylinder axis");
//   // local z defined by x and y (should point outwards from axis)
//   GlobalVector zPlane(xPlane.cross(yPlane));
  // rotation constructor will normalize...
  return ReferenceCountingPointer<TangentPlane>(new TangentPlane(aPoint,
                                 RotationType(xPlane,
                                          yPlane),
                                 this));
}