#include <UGenericPolycone.hh>

Inheritance diagram for UGenericPolycone:

Public Member Functions
	UGenericPolycone (const std::string &name, double phiStart, double phiTotal, int numZPlanes, const double zPlane[], const double rInner[], const double rOuter[])

	UGenericPolycone (const std::string &name, double phiStart, double phiTotal, int numRZ, const double r[], const double z[])

virtual	~UGenericPolycone ()

VUSolid::EnumInside	Inside (const UVector3 &p) const

double	DistanceToIn (const UVector3 &p, const UVector3 &v, double aPstep=UUtils::kInfinity) const

UVector3	GetPointOnSurface () const

UGeometryType	GetEntityType () const

VUSolid *	Clone () const

std::ostream &	StreamInfo (std::ostream &os) const

bool	Reset ()

double	GetStartPhi () const

double	GetEndPhi () const

bool	IsOpen () const

int	GetNumRZCorner () const

UPolyconeSideRZ	GetCorner (int index) const

	UGenericPolycone (const UGenericPolycone &source)

UGenericPolycone &	operator= (const UGenericPolycone &source)

Public Member Functions inherited from UVCSGfaceted
	UVCSGfaceted (const std::string &name)

virtual	~UVCSGfaceted ()

	UVCSGfaceted (const UVCSGfaceted &source)

UVCSGfaceted &	operator= (const UVCSGfaceted &source)

VUSolid::EnumInside	InsideNoVoxels (const UVector3 &p) const

virtual bool	Normal (const UVector3 &p, UVector3 &n) const

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

virtual double	SafetyFromOutside (const UVector3 &aPoint, bool aAccurate=false) const

double	DistanceTo (const UVector3 &p, const bool outgoing) const

double	DistanceToOutNoVoxels (const UVector3 &p, const UVector3 &v, UVector3 &n, bool &aConvex) const

virtual double	DistanceToOut (const UVector3 &p, const UVector3 &v, UVector3 &n, bool &aConvex, double aPstep=UUtils::kInfinity) const

virtual double	SafetyFromInside (const UVector3 &aPoint, bool aAccurate=false) const

virtual double	SafetyFromInsideNoVoxels (const UVector3 &aPoint, bool aAccurate=false) const

int	GetCubVolStatistics () const

double	GetCubVolEpsilon () const

void	SetCubVolStatistics (int st)

void	SetCubVolEpsilon (double ep)

int	GetAreaStatistics () const

double	GetAreaAccuracy () const

void	SetAreaStatistics (int st)

void	SetAreaAccuracy (double ep)

virtual double	Capacity ()

virtual double	SurfaceArea ()

Public Member Functions inherited from VUSolid
	VUSolid ()

	VUSolid (const std::string &name)

virtual	~VUSolid ()

double	GetCarTolerance () const

double	GetRadTolerance () const

double	GetAngTolerance () const

void	SetCarTolerance (double eps)

void	SetRadTolerance (double eps)

void	SetAngTolerance (double eps)

virtual void	ExtentAxis (EAxisType aAxis, double &aMin, double &aMax) const

const std::string &	GetName () const

void	SetName (const std::string &aName)

virtual void	SamplePointsInside (int, UVector3 *) const

virtual void	SamplePointsOnSurface (int, UVector3 *) const

virtual void	SamplePointsOnEdge (int, UVector3 *) const

double	EstimateCubicVolume (int nStat, double epsilon) const

double	EstimateSurfaceArea (int nStat, double ell) const

Protected Member Functions
void	Create (double phiStart, double phiTotal, UReduciblePolygon *rz)

void	CopyStuff (const UGenericPolycone &source)

void	GetParametersList (int, double *) const

void	ComputeBBox (UBBox *, bool)

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

Protected Member Functions inherited from UVCSGfaceted
double	SafetyFromInsideSection (int index, const UVector3 &p, UBits &bits) const

int	GetSection (double z) const

UVector3	GetPointOnSurfaceGeneric () const

void	CopyStuff (const UVCSGfaceted &source)

void	DeleteStuff ()

void	FindCandidates (double z, std::vector< int > &candidates, bool sides=false)

void	InitVoxels (UReduciblePolygon &z, double radius)

Protected Attributes
double	startPhi

double	endPhi

bool	phiIsOpen

int	numCorner

UPolyconeSideRZ *	corners

UEnclosingCylinder *	enclosingCylinder

Protected Attributes inherited from UVCSGfaceted
int	numFace

UVCSGface **	faces

double	fCubicVolume

double	fSurfaceArea

std::vector< double >	fZs

std::vector< std::vector< int > >	fCandidates

int	fMaxSection

UBox	fBox

double	fBoxShift

bool	fNoVoxels

Additional Inherited Members
Public Types inherited from VUSolid
enum	EnumInside { eInside =0, eSurface =1, eOutside =2 }

enum	EAxisType { eXaxis =0, eYaxis =1, eZaxis =2 }

Static Public Member Functions inherited from VUSolid
static double	Tolerance ()

Static Protected Attributes inherited from VUSolid
static double	fgTolerance = 1.0E-9

static double	frTolerance = 1.0E-9

static double	faTolerance = 1.0E-9

Detailed Description

Definition at line 39 of file UGenericPolycone.hh.

Constructor & Destructor Documentation

UGenericPolycone::UGenericPolycone	(	const std::string &	name,
		double	phiStart,
		double	phiTotal,
		int	numZPlanes,
		const double	zPlane[],
		const double	rInner[],
		const double	rOuter[]
	)

Referenced by Clone().

UGenericPolycone::UGenericPolycone	(	const std::string &	name,
		double	phiStart,
		double	phiTotal,
		int	numRZ,
		const double	r[],
		const double	z[]
	)

Definition at line 37 of file UGenericPolycone.cc.

References Create().

   : UVCSGfaceted(name)
 {
   UReduciblePolygon* rz = new UReduciblePolygon(r, z, numRZ);
 
   Create(phiStart, phiTotal, rz);
 
   delete rz;
 }

UGenericPolycone::~UGenericPolycone ( )

virtual

Definition at line 247 of file UGenericPolycone.cc.

References corners, and enclosingCylinder.

 {
   delete [] corners;
   delete enclosingCylinder;
 }

UGenericPolycone::UGenericPolycone ( const UGenericPolycone & source )

Definition at line 257 of file UGenericPolycone.cc.

References CopyStuff().

   : UVCSGfaceted(source)
 {
   CopyStuff(source);
 }

Member Function Documentation

VUSolid * UGenericPolycone::Clone ( ) const

virtual

Implements VUSolid.

Definition at line 384 of file UGenericPolycone.cc.

References UGenericPolycone().

 {
   return new UGenericPolycone(*this);
 }

void UGenericPolycone::ComputeBBox	(	UBBox *	,
		bool
	)

inlineprotectedvirtual

Implements VUSolid.

Definition at line 121 of file UGenericPolycone.hh.

     {
       // Computes bounding box.
       std::cout << "ComputeBBox - Not implemented" << std::endl;
     }

void UGenericPolycone::CopyStuff ( const UGenericPolycone & source )

protected

Definition at line 286 of file UGenericPolycone.cc.

References corners, enclosingCylinder, endPhi, numCorner, phiIsOpen, and startPhi.

Referenced by operator=(), and UGenericPolycone().

 {
   //
   // Simple stuff
   //
   startPhi  = source.startPhi;
   endPhi    = source.endPhi;
   phiIsOpen = source.phiIsOpen;
   numCorner = source.numCorner;
 
   //
   // The corner array
   //
   corners = new UPolyconeSideRZ[numCorner];
 
   UPolyconeSideRZ* corn = corners,
                    *sourceCorn = source.corners;
   do
   {
     *corn = *sourceCorn;
   }
   while (++sourceCorn, ++corn < corners + numCorner);
 
   //
   // Enclosing cylinder
   //
   enclosingCylinder = new UEnclosingCylinder(*source.enclosingCylinder);
 }

void UGenericPolycone::Create	(	double	phiStart,
		double	phiTotal,
		UReduciblePolygon *	rz
	)

protected

Definition at line 59 of file UGenericPolycone.cc.

Referenced by UGenericPolycone().

 {
   //
   // Perform checks of rz values
   //
   if (rz->Amin() < 0.0)
   {
     std::ostringstream message;
     message << "Illegal input parameters - " << GetName() << std::endl
             << "                All R values must be >= 0 !";
     UUtils::Exception("UGenericPolycone::Create()", "GeomSolids0002",
                       FatalErrorInArguments, 1, message.str().c_str());
   }
 
   double rzArea = rz->Area();
   if (rzArea < -VUSolid::Tolerance())
     rz->ReverseOrder();
 
   else if (rzArea < -VUSolid::Tolerance())
   {
     std::ostringstream message;
     message << "Illegal input parameters - " << GetName() << std::endl
             << "                R/Z Cross section is zero or near zero: " << rzArea;
     UUtils::Exception("UGenericPolycone::Create()", "GeomSolids0002",
                       FatalErrorInArguments, 1, message.str().c_str());
   }
 
   if ((!rz->RemoveDuplicateVertices(VUSolid::Tolerance()))
       || (!rz->RemoveRedundantVertices(VUSolid::Tolerance())))
   {
     std::ostringstream message;
     message << "Illegal input parameters - " << GetName() << std::endl
             << "                Too few unique R/Z values !";
     UUtils::Exception("UGenericPolycone::Create()", "GeomSolids0002",
                       FatalErrorInArguments, 1, message.str().c_str());
   }
 
   if (rz->CrossesItself(1 / UUtils::kInfinity))
   {
     std::ostringstream message;
     message << "Illegal input parameters - " << GetName() << std::endl
             << "                R/Z segments Cross !";
     UUtils::Exception("UGenericPolycone::Create()", "GeomSolids0002",
                       FatalErrorInArguments, 1, message.str().c_str());
   }
 
   numCorner = rz->NumVertices();
 
   //
   // Phi opening? Account for some possible roundoff, and interpret
   // nonsense value as representing no phi opening
   //
   if (phiTotal <= 0 || phiTotal > 2 * UUtils::kPi - 1E-10)
   {
     phiIsOpen = false;
     startPhi = 0;
     endPhi = 2 * UUtils::kPi;
   }
   else
   {
     phiIsOpen = true;
 
     //
     // Convert phi into our convention
     //
     startPhi = phiStart;
     while (startPhi < 0) startPhi += 2 * UUtils::kPi;
 
     endPhi = phiStart + phiTotal;
     while (endPhi < startPhi) endPhi += 2 * UUtils::kPi;
   }
 
   //
   // Allocate corner array.
   //
   corners = new UPolyconeSideRZ[numCorner];
 
   //
   // Copy corners
   //
   UReduciblePolygonIterator iterRZ(rz);
 
   UPolyconeSideRZ* next = corners;
   iterRZ.Begin();
   do
   {
     next->r = iterRZ.GetA();
     next->z = iterRZ.GetB();
   }
   while (++next, iterRZ.Next());
 
   //
   // Allocate face pointer array
   //
   numFace = phiIsOpen ? numCorner + 2 : numCorner;
   faces = new UVCSGface*[numFace];
 
   //
   // Construct conical faces
   //
   // But! Don't construct a face if both points are at zero radius!
   //
   UPolyconeSideRZ* corner = corners,
                    *prev = corners + numCorner - 1,
                     *nextNext;
   UVCSGface** face = faces;
   do
   {
     next = corner + 1;
     if (next >= corners + numCorner) next = corners;
     nextNext = next + 1;
     if (nextNext >= corners + numCorner) nextNext = corners;
 
     if (corner->r < 1 / UUtils::kInfinity && next->r < 1 / UUtils::kInfinity) continue;
 
     //
     // We must decide here if we can dare declare one of our faces
     // as having a "valid" normal (i.e. allBehind = true). This
     // is never possible if the face faces "inward" in r.
     //
     bool allBehind;
     if (corner->z > next->z)
     {
       allBehind = false;
     }
     else
     {
       //
       // Otherwise, it is only true if the line passing
       // through the two points of the segment do not
       // split the r/z Cross section
       //
       allBehind = !rz->BisectedBy(corner->r, corner->z,
                                   next->r, next->z, VUSolid::Tolerance());
     }
 
     *face++ = new UPolyconeSide(prev, corner, next, nextNext,
                                 startPhi, endPhi - startPhi, phiIsOpen, allBehind);
   }
   while (prev = corner, corner = next, corner > corners);
 
   if (phiIsOpen)
   {
     //
     // Construct phi open edges
     //
     *face++ = new UPolyPhiFace(rz, startPhi, 0, endPhi);
     *face++ = new UPolyPhiFace(rz, endPhi,   0, startPhi);
   }
 
   //
   // We might have dropped a face or two: recalculate numFace
   //
   numFace = face - faces;
 
   //
   // Make enclosingCylinder
   //
   enclosingCylinder =
     new UEnclosingCylinder(rz->Amax(), rz->Bmax(), rz->Bmin(), phiIsOpen, phiStart, phiTotal);
 
   InitVoxels(*rz, enclosingCylinder->radius);
 
   fNoVoxels = fMaxSection < 2;
 
 }

double UGenericPolycone::DistanceToIn	(	const UVector3 &	p,
		const UVector3 &	v,
		double	aPstep = `UUtils::kInfinity`
	)		const

virtual

Reimplemented from UVCSGfaceted.

Definition at line 358 of file UGenericPolycone.cc.

References UVCSGfaceted::DistanceToIn(), enclosingCylinder, and UEnclosingCylinder::ShouldMiss().

 {
   //
   // Quick test
   //
   if (enclosingCylinder->ShouldMiss(p, v))
     return UUtils::kInfinity;
 
   //
   // Long answer
   //
   return UVCSGfaceted::DistanceToIn(p, v, aPstep);
 }

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

protectedvirtual

Implements VUSolid.

Definition at line 427 of file UGenericPolycone.cc.

References enclosingCylinder, and UEnclosingCylinder::Extent().

 {
   enclosingCylinder->Extent(aMin, aMax);
 }

UPolyconeSideRZ UGenericPolycone::GetCorner ( int index ) const

inline

Referenced by G4UGenericPolycone::GetCorner().

double UGenericPolycone::GetEndPhi ( ) const

inline

Referenced by G4UGenericPolycone::GetEndPhi().

UGeometryType UGenericPolycone::GetEntityType ( ) const

virtual

Reimplemented from UVCSGfaceted.

Definition at line 376 of file UGenericPolycone.cc.

 {
   return std::string("GenericPolycone");
 }

int UGenericPolycone::GetNumRZCorner ( ) const

inline

Referenced by G4UGenericPolycone::GetNumRZCorner().

void UGenericPolycone::GetParametersList	(	int	,
		double *
	)		const

inlineprotectedvirtual

Implements VUSolid.

Definition at line 119 of file UGenericPolycone.hh.

119 {}

UVector3 UGenericPolycone::GetPointOnSurface ( ) const

virtual

Implements VUSolid.

Definition at line 421 of file UGenericPolycone.cc.

References UVCSGfaceted::GetPointOnSurfaceGeneric().

 {
   return GetPointOnSurfaceGeneric();
 
 }

double UGenericPolycone::GetStartPhi ( ) const

inline

Referenced by G4UGenericPolycone::GetStartPhi().

VUSolid::EnumInside UGenericPolycone::Inside ( const UVector3 & p ) const

virtual

Reimplemented from UVCSGfaceted.

Definition at line 338 of file UGenericPolycone.cc.

References enclosingCylinder, VUSolid::eOutside, UVCSGfaceted::Inside(), and UEnclosingCylinder::MustBeOutside().

 {
   //
   // Quick test
   //
   if (enclosingCylinder->MustBeOutside(p)) return eOutside;
 
   //
   // Long answer
   //
   return UVCSGfaceted::Inside(p);
 }

bool UGenericPolycone::IsOpen ( ) const

inline

Referenced by G4UGenericPolycone::IsOpen().

UGenericPolycone & UGenericPolycone::operator= ( const UGenericPolycone & source )

Definition at line 267 of file UGenericPolycone.cc.

References CopyStuff(), corners, enclosingCylinder, and UVCSGfaceted::operator=().

 {
   if (this == &source) return *this;
 
   UVCSGfaceted::operator=(source);
 
   delete [] corners;
 
   delete enclosingCylinder;
 
   CopyStuff(source);
 
   return *this;
 }

bool UGenericPolycone::Reset ( )

Definition at line 319 of file UGenericPolycone.cc.

References VUSolid::GetName().

 {
 
   std::ostringstream message;
   message << "Solid " << GetName() << " built using generic construct."
           << std::endl << "Not applicable to the generic construct !";
   // UException("UGenericPolycone::Reset(,,)", "GeomSolids1001",
   //            JustWarning, message, "Parameters NOT resetted.");
   return 1;
 
 }

std::ostream & UGenericPolycone::StreamInfo ( std::ostream & os ) const

virtual

Reimplemented from UVCSGfaceted.

Definition at line 392 of file UGenericPolycone.cc.

References corners, endPhi, VUSolid::GetName(), numCorner, UPolyconeSideRZ::r, startPhi, and UPolyconeSideRZ::z.

 {
   int oldprc = os.precision(16);
   os << "-----------------------------------------------------------\n"
      << "       *** Dump for solid - " << GetName() << " ***\n"
      << "       ===================================================\n"
      << " Solid type: UGenericPolycone\n"
      << " Parameters: \n"
      << "       starting phi angle : " << startPhi / (UUtils::kPi / 180.0) << " degrees \n"
      << "       ending phi angle     : " << endPhi / (UUtils::kPi / 180.0) << " degrees \n";
   int i = 0;
 
   os << "    number of RZ points: " << numCorner << "\n"
      << "              RZ values (corners): \n";
   for (i = 0; i < numCorner; i++)
   {
     os << "                         "
        << corners[i].r << ", " << corners[i].z << "\n";
   }
 
   os << "-----------------------------------------------------------\n";
   os.precision(oldprc);
 
   return os;
 }