#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 =0.0

G4double	ry1 =0.0

G4double	dx1 =0.0

G4double	dy1 =0.0

G4double	rx2 =0.0

G4double	ry2 =0.0

G4double	dx2 =0.0

G4double	dy2 =0.0

G4bool	concave

Detailed Description

Definition at line 45 of file G4EnclosingCylinder.hh.

Constructor & Destructor Documentation

◆ G4EnclosingCylinder() [1/2]

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

Definition at line 39 of file G4EnclosingCylinder.cc.

  : startPhi(theStartPhi), totalPhi(theTotalPhi),
    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 91 of file G4EnclosingCylinder.cc.

92{

93}

◆ G4EnclosingCylinder() [2/2]

G4EnclosingCylinder::G4EnclosingCylinder ( __void__ & )

Definition at line 83 of file G4EnclosingCylinder.cc.

: radius(0.), zLo(0.), zHi(0.), phiIsOpen(0.), startPhi(0.), totalPhi(0.),
  concave(false)
{
}

Member Function Documentation

◆ MustBeOutside()

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

Definition at line 101 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 G4GenericPolycone::Inside(), G4Polycone::Inside(), G4Polyhedra::Inside(), and ShouldMiss().

◆ ShouldMiss()

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

Definition at line 130 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;
}