Cone.cc

Go to the documentation of this file.

 
 
#include "DataFormats/GeometrySurface/interface/Cone.h"
#include "DataFormats/GeometrySurface/interface/TangentPlane.h"
 
#include <iostream>
 
ConstReferenceCountingPointer<TangentPlane> Cone::tangentPlane(const GlobalPoint&) const {
  // FIXME: to be implemented...
  std::cout << "*** WARNING: Cone::tangentPlane not implemented." << std::endl;
  abort();
  return ConstReferenceCountingPointer<TangentPlane>();
}
 
ConstReferenceCountingPointer<TangentPlane> Cone::tangentPlane(const LocalPoint&) const {
  // FIXME: to be implemented...
  std::cout << "*** WARNING: Cone::tangentPlane not implemented." << std::endl;
  abort();
  return ConstReferenceCountingPointer<TangentPlane>();
}
 
Surface::Side Cone::side(const GlobalPoint& p, Scalar tolerance) const {
  // FIXME: should be done in local coordinates as this is not correct in the case the verstex is not on the (global) Z axis!!!!
 
  // tolerance is interpreted as max distance from cone surface.
  // FIXME: check case when vertex().z()==inf.
  GlobalPoint p1(p.x(), p.y(), p.z() - vertex().z());
 
  // handle the singularity of p=vertex (i.e. p1.mag() undefined)
  if (p1.mag() < tolerance)
    return SurfaceOrientation::onSurface;
  double delta = double(p1.theta()) - double(openingAngle());
  if (fabs(delta) < tolerance / p1.mag())
    return SurfaceOrientation::onSurface;
 
  if (p1.theta() < Geom::pi() / 2.) {
    return (delta > 0. ? SurfaceOrientation::positiveSide : SurfaceOrientation::negativeSide);
  } else {
    return (delta > 0. ? SurfaceOrientation::negativeSide : SurfaceOrientation::positiveSide);
  }
}