|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
#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));
}
ConstReferenceCountingPointer<TangentPlane> Cylinder::tangentPlane(const LocalPoint& aPoint) const {
return tangentPlane(toGlobal(aPoint));
}
ConstReferenceCountingPointer<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 ConstReferenceCountingPointer<TangentPlane>(new TangentPlane(aPoint, RotationType(xPlane, yPlane)));
}
|