#include <G4EnclosingCylinder.hh>

Public Member Functions
	G4EnclosingCylinder (const G4ReduciblePolygon *rz, G4bool phiIsOpen, G4double startPhi, G4double totalPhi)

	~G4EnclosingCylinder ()

G4bool	MustBeOutside (const G4ThreeVector &p) const

G4bool	ShouldMiss (const G4ThreeVector &p, const G4ThreeVector &v) const

	G4EnclosingCylinder (__void__ &)

Protected Attributes
G4double	radius

G4double	zLo

G4double	zHi

G4bool	phiIsOpen

G4double	startPhi

G4double	totalPhi

G4double	rx1

G4double	ry1

G4double	dx1

G4double	dy1

G4double	rx2

G4double	ry2

G4double	dx2

G4double	dy2

G4bool	concave

Detailed Description

Definition at line 53 of file G4EnclosingCylinder.hh.

Constructor & Destructor Documentation

◆ G4EnclosingCylinder() [1/2]

G4EnclosingCylinder::G4EnclosingCylinder	(	const G4ReduciblePolygon *	rz,
		G4bool	phiIsOpen,
		G4double	startPhi,
		G4double	totalPhi
	)

Definition at line 48 of file G4EnclosingCylinder.cc.

  : startPhi(theStartPhi), totalPhi(theTotalPhi),
    rx1(0.), ry1(0.), dx1(0.), dy1(0.),
    rx2(0.), ry2(0.), dx2(0.), dy2(0.),     
    concave(theTotalPhi > pi)
{
  //
  // Obtain largest r and smallest and largest z
  //
  radius = rz->Amax();
  zHi = rz->Bmax();
  zLo = rz->Bmin();
 
  G4double kCarTolerance = G4GeometryTolerance::GetInstance()
                           ->GetSurfaceTolerance();
  //
  // Save phi info
  //
  phiIsOpen = thePhiIsOpen;
  if ( phiIsOpen )
  {    
    rx1 = std::cos(startPhi);
    ry1 = std::sin(startPhi);
    dx1 = +ry1*10*kCarTolerance;
    dy1 = -rx1*10*kCarTolerance;
    
    rx2 = std::cos(startPhi+totalPhi);
    ry2 = std::sin(startPhi+totalPhi);
    dx2 = -ry2*10*kCarTolerance;
    dy2 = +rx2*10*kCarTolerance;
  }
  
  //
  // Add safety
  //
  radius += 10*kCarTolerance;
  zLo    -= 10*kCarTolerance;
  zHi    += 10*kCarTolerance;
}

◆ ~G4EnclosingCylinder()

G4EnclosingCylinder::~G4EnclosingCylinder ( )

Definition at line 105 of file G4EnclosingCylinder.cc.

106{

107}

◆ G4EnclosingCylinder() [2/2]

G4EnclosingCylinder::G4EnclosingCylinder ( __void__ & )

Definition at line 95 of file G4EnclosingCylinder.cc.

: radius(0.), zLo(0.), zHi(0.), phiIsOpen(0.), startPhi(0.), totalPhi(0.),
  rx1(0.), ry1(0.), dx1(0.), dy1(0.), rx2(0.), ry2(0.), dx2(0.), dy2(0.),
  concave(false)
{
}

Member Function Documentation

◆ MustBeOutside()

G4bool G4EnclosingCylinder::MustBeOutside ( const G4ThreeVector & p ) const

Definition at line 117 of file G4EnclosingCylinder.cc.

{
  if (p.perp() > 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;
}

Referenced by G4Polycone::Inside(), G4Polyhedra::Inside(), and ShouldMiss().

◆ ShouldMiss()

G4bool G4EnclosingCylinder::ShouldMiss	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v
	)		const

Definition at line 148 of file G4EnclosingCylinder.cc.

{
  if (!MustBeOutside(p)) return false;
  
  G4double cross = p.x()*v.y() - p.y()*v.x();
  if (cross > radius) return true;
  
  if (p.perp() > radius)
  {
    G4double dot = p.x()*v.x() + p.y()*v.y();
    if (dot > 0) return true;
  }
 
  return false;
}