larcorealg/larcorealg/GeoAlgo/GeoCone.cxx

Go to the documentation of this file.

#include "larcorealg/GeoAlgo/GeoAlgoException.h"  // for GeoAlgoException
#include "larcorealg/GeoAlgo/GeoCone.h"

#include <math.h>
#include <sstream>

namespace geoalgo {

  Cone::Cone() : HalfLine()
  {
    _length = 1;
    _radius = 1;
    _angle = atan(_radius/_length);
  }

  Cone::Cone(const double x,    const double y,    const double z,
             const double dirx, const double diry, const double dirz,
             const double length, const double radius)
    : HalfLine(x, y, z, dirx, diry, dirz)
  {
    if (length == 0){
      std::ostringstream msg;
      msg << "<<" << __FUNCTION__ << ">>"
          << " Cone Length cannot be 0." << std::endl;
      throw GeoAlgoException(msg.str());
    }
    _length = length;
    _radius = radius;
    _angle = atan(_radius/_length);
  }

  Cone::Cone(const Point_t& start, const Vector_t& dir,
             const double length, const double radius)
    : HalfLine( start,  dir )
  {
    if (length == 0){
      std::ostringstream msg;
      msg << "<<" << __FUNCTION__ << ">>"
          << " Cone Length cannot be 0." << std::endl;
      throw GeoAlgoException(msg.str());
    }
    _length = length;
    _radius = radius;
    _angle  = atan(_radius/_length);
  }

  double Cone::Length() const { return _length; }

  double Cone::Radius() const { return _radius; }

  double Cone::Angle() const { return _angle; }

  void Cone::Length(const double l)
  {
    if (l == 0){
      std::ostringstream msg;
      msg << "<<" << __FUNCTION__ << ">>"
          << " Cone Length cannot be 0." << std::endl;
      throw GeoAlgoException(msg.str());
    }
    _length = l;
    _angle = atan(_radius/_length);
  }

  void Cone::Radius(const double r) { _radius = r; _angle = atan(_radius/_length); }
}