#include <UEnclosingCylinder.hh>

Public Member Functions
	UEnclosingCylinder (double r, double lo, double hi, bool phiIsOpen, double startPhi, double totalPhi)

	~UEnclosingCylinder ()

bool	MustBeOutside (const UVector3 &p) const

bool	ShouldMiss (const UVector3 &p, const UVector3 &v) const

double	DistanceTo (const UVector3 &p, const UVector3 &v) const

double	SafetyFromOutside (const UVector3 &p) const

void	Extent (UVector3 &aMin, UVector3 &aMax) const

Data Fields
double	radius

Protected Attributes
double	zLo

double	zHi

bool	phiIsOpen

double	startPhi

double	totalPhi

double	rx1

double	ry1

double	dx1

double	dy1

double	rx2

double	ry2

double	dx2

double	dy2

bool	concave

UTubs *	tube

Detailed Description

Definition at line 31 of file UEnclosingCylinder.hh.

Constructor & Destructor Documentation

UEnclosingCylinder::UEnclosingCylinder	(	double	r,
		double	lo,
		double	hi,
		bool	phiIsOpen,
		double	startPhi,
		double	totalPhi
	)

Definition at line 31 of file UEnclosingCylinder.cc.

References dx1, dx2, dy1, dy2, phiIsOpen, radius, rx1, rx2, ry1, ry2, startPhi, VUSolid::Tolerance(), totalPhi, tube, zHi, and zLo.

   : startPhi(theStartPhi), totalPhi(theTotalPhi),
     rx1(0.), ry1(0.), dx1(0.), dy1(0.),
     rx2(0.), ry2(0.), dx2(0.), dy2(0.),
     concave(theTotalPhi > UUtils::kPi)
 {
   //
   // Obtain largest r and smallest and largest z
   //
 
   /*
   radius = rz->Amax();
   zHi = rz->Bmax();
   zLo = rz->Bmin();
   */
 
   radius = r;
   zHi = hi;
   zLo = lo;
 
   double fTolerance = VUSolid::Tolerance();
 
   tube = new UTubs("", 0, radius + fTolerance, zHi - zLo, theStartPhi, theTotalPhi);
 
   //
   // Save phi info
   //
   phiIsOpen = thePhiIsOpen;
   if (phiIsOpen)
   {
     rx1 = std::cos(startPhi);
     ry1 = std::sin(startPhi);
     dx1 = +ry1 * 10 * fTolerance;
     dy1 = -rx1 * 10 * fTolerance;
 
     rx2 = std::cos(startPhi + totalPhi);
     ry2 = std::sin(startPhi + totalPhi);
     dx2 = -ry2 * 10 * fTolerance;
     dy2 = +rx2 * 10 * fTolerance;
   }
 
   //
   // Add safety
   //
   radius += 10 * fTolerance;
   zLo   -= 10 * fTolerance;
   zHi   += 10 * fTolerance;
 }

UEnclosingCylinder::~UEnclosingCylinder ( )

Definition at line 87 of file UEnclosingCylinder.cc.

88 {

89 }

Member Function Documentation

double UEnclosingCylinder::DistanceTo	(	const UVector3 &	p,
		const UVector3 &	v
	)		const

Definition at line 179 of file UEnclosingCylinder.cc.

References UTubs::DistanceToIn(), and tube.

 {
   return tube->DistanceToIn(p, v);
 }

void UEnclosingCylinder::Extent	(	UVector3 &	aMin,
		UVector3 &	aMax
	)		const

Definition at line 173 of file UEnclosingCylinder.cc.

References radius, zHi, and zLo.

Referenced by UGenericPolycone::Extent(), and UPolyhedra::Extent().

 {
   aMin = UVector3(-radius, -radius, zLo);
   aMax = UVector3(radius, radius, zHi);
 }

bool UEnclosingCylinder::MustBeOutside ( const UVector3 & p ) const

Definition at line 99 of file UEnclosingCylinder.cc.

References concave, dx1, dx2, dy1, dy2, UVector3::Perp2(), phiIsOpen, radius, rx1, rx2, ry1, ry2, UVector3::x, UVector3::y, UVector3::z, zHi, and zLo.

Referenced by UGenericPolycone::Inside(), UPolyhedra::Inside(), and ShouldMiss().

 {
 //  if (p.Perp() > radius ) return true;
   if (p.Perp2() > radius * radius) return true;
   if (p.z < zLo) return true;
   if (p.z > zHi) return true;
 
   if (phiIsOpen)
   {
     if (concave)
     {
       if (((p.x - dx1)*ry1 - (p.y - dy1)*rx1) < 0) return false;
       if (((p.x - dx2)*ry2 - (p.y - dy2)*rx2) > 0) return false;
     }
     else
     {
       if (((p.x - dx1)*ry1 - (p.y - dy1)*rx1) > 0) return true;
       if (((p.x - dx2)*ry2 - (p.y - dy2)*rx2) < 0) return true;
     }
   }
 
   return false;
 }

double UEnclosingCylinder::SafetyFromOutside ( const UVector3 & p ) const

Definition at line 184 of file UEnclosingCylinder.cc.

References UTubs::SafetyFromOutside(), and tube.

Referenced by UPolycone::SafetyFromOutside().

 {
   return tube->SafetyFromOutside(p);
 }

bool UEnclosingCylinder::ShouldMiss	(	const UVector3 &	p,
		const UVector3 &	v
	)		const

Definition at line 131 of file UEnclosingCylinder.cc.

References MustBeOutside(), UVector3::Perp2(), radius, UVector3::x, and UVector3::y.

Referenced by UGenericPolycone::DistanceToIn(), and UPolyhedra::DistanceToIn().

 {
   if (!MustBeOutside(p)) return false;
 
 //  if (p.z < zLo - VUSolid::Tolerance() && v.z <= 0) return true;
 //  if (p.z > zHi + VUSolid::Tolerance() && v.z >= 0) return true;
 
   double cross = p.x * v.y - p.y * v.x;
   if (cross > radius) return true;
 
   double r2 = p.Perp2();
   if (r2 > radius * radius)
   {
     double dot = p.x * v.x + p.y * v.y;
     if (dot > 0) return true;
 
 //    double n = v.Perp2();
 //    if (Dot < std::sqrt(r2 - radius * radius) * n) return true;
   }
 
 
   /*
   if (phiIsOpen)
   {
     if (concave)
     {
       if ( ((p.x-dx1)*ry1 - (p.y-dy1)*rx1) < 0) return false;
       if ( ((p.x-dx2)*ry2 - (p.y-dy2)*rx2) > 0) return false;
     }
     else
     {
       if ( ((p.x-dx1)*ry1 - (p.y-dy1)*rx1) > 0) return true;
       if ( ((p.x-dx2)*ry2 - (p.y-dy2)*rx2) < 0) return true;
     }
     return false;
   }
   */
 
   return false;
 }