#include <G4VCSGfaceted.hh>

Inheritance diagram for G4VCSGfaceted:

Public Member Functions
	G4VCSGfaceted (const G4String &name)

virtual	~G4VCSGfaceted ()

	G4VCSGfaceted (const G4VCSGfaceted &source)

const G4VCSGfaceted &	operator= (const G4VCSGfaceted &source)

virtual G4bool	CalculateExtent (const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pmin, G4double &pmax) const

virtual EInside	Inside (const G4ThreeVector &p) const

virtual G4ThreeVector	SurfaceNormal (const G4ThreeVector &p) const

virtual G4double	DistanceToIn (const G4ThreeVector &p, const G4ThreeVector &v) const

virtual G4double	DistanceToIn (const G4ThreeVector &p) const

virtual G4double	DistanceToOut (const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcNorm=false, G4bool validNorm=0, G4ThreeVector n=0) const

virtual G4double	DistanceToOut (const G4ThreeVector &p) const

virtual G4GeometryType	GetEntityType () const

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

virtual G4Polyhedron *	CreatePolyhedron () const =0

virtual void	DescribeYourselfTo (G4VGraphicsScene &scene) const

virtual G4VisExtent	GetExtent () const

virtual G4Polyhedron *	GetPolyhedron () const

G4int	GetCubVolStatistics () const

G4double	GetCubVolEpsilon () const

void	SetCubVolStatistics (G4int st)

void	SetCubVolEpsilon (G4double ep)

G4int	GetAreaStatistics () const

G4double	GetAreaAccuracy () const

void	SetAreaStatistics (G4int st)

void	SetAreaAccuracy (G4double ep)

virtual G4double	GetCubicVolume ()

virtual G4double	GetSurfaceArea ()

	G4VCSGfaceted (__void__ &)

Public Member Functions inherited from G4VSolid
	G4VSolid (const G4String &name)

virtual	~G4VSolid ()

G4bool	operator== (const G4VSolid &s) const

G4String	GetName () const

void	SetName (const G4String &name)

G4double	GetTolerance () const

virtual void	ComputeDimensions (G4VPVParameterisation p, const G4int n, const G4VPhysicalVolume pRep)

virtual G4ThreeVector	GetPointOnSurface () const

virtual G4VSolid *	Clone () const

void	DumpInfo () const

virtual const G4VSolid *	GetConstituentSolid (G4int no) const

virtual G4VSolid *	GetConstituentSolid (G4int no)

virtual const G4DisplacedSolid *	GetDisplacedSolidPtr () const

virtual G4DisplacedSolid *	GetDisplacedSolidPtr ()

	G4VSolid (__void__ &)

	G4VSolid (const G4VSolid &rhs)

G4VSolid &	operator= (const G4VSolid &rhs)

Protected Member Functions
virtual G4double	DistanceTo (const G4ThreeVector &p, const G4bool outgoing) const

G4ThreeVector	GetPointOnSurfaceGeneric () const

void	CopyStuff (const G4VCSGfaceted &source)

void	DeleteStuff ()

Protected Member Functions inherited from G4VSolid
void	CalculateClippedPolygonExtent (G4ThreeVectorList &pPolygon, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipCrossSection (G4ThreeVectorList *pVertices, const G4int pSectionIndex, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipBetweenSections (G4ThreeVectorList *pVertices, const G4int pSectionIndex, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis, G4double &pMin, G4double &pMax) const

void	ClipPolygon (G4ThreeVectorList &pPolygon, const G4VoxelLimits &pVoxelLimit, const EAxis pAxis) const

G4double	EstimateCubicVolume (G4int nStat, G4double epsilon) const

G4double	EstimateSurfaceArea (G4int nStat, G4double ell) const

Protected Attributes
G4int	numFace

G4VCSGface **	faces

G4double	fCubicVolume

G4double	fSurfaceArea

G4Polyhedron *	fpPolyhedron

Protected Attributes inherited from G4VSolid
G4double	kCarTolerance

Detailed Description

Definition at line 53 of file G4VCSGfaceted.hh.

Constructor & Destructor Documentation

G4VCSGfaceted::G4VCSGfaceted ( const G4String & name )

Definition at line 62 of file G4VCSGfaceted.cc.

   : G4VSolid(name),
     numFace(0), faces(0), fCubicVolume(0.), fSurfaceArea(0.), fpPolyhedron(0),
     fStatistics(1000000), fCubVolEpsilon(0.001), fAreaAccuracy(-1.)
 {
 }

G4VCSGfaceted::~G4VCSGfaceted ( )

virtual

Definition at line 84 of file G4VCSGfaceted.cc.

References DeleteStuff(), and fpPolyhedron.

 {
   DeleteStuff();
   delete fpPolyhedron;
 }

G4VCSGfaceted::G4VCSGfaceted ( const G4VCSGfaceted & source )

Definition at line 94 of file G4VCSGfaceted.cc.

References CopyStuff().

   : G4VSolid( source )
 {
   fStatistics = source.fStatistics;
   fCubVolEpsilon = source.fCubVolEpsilon;
   fAreaAccuracy = source.fAreaAccuracy;
 
   CopyStuff( source );
 }

G4VCSGfaceted::G4VCSGfaceted ( __void__ & a )

Definition at line 74 of file G4VCSGfaceted.cc.

   : G4VSolid(a),
     numFace(0), faces(0), fCubicVolume(0.), fSurfaceArea(0.), fpPolyhedron(0),
     fStatistics(1000000), fCubVolEpsilon(0.001), fAreaAccuracy(-1.)
 {
 }

Member Function Documentation

G4bool G4VCSGfaceted::CalculateExtent	(	const EAxis	pAxis,
		const G4VoxelLimits &	pVoxelLimit,
		const G4AffineTransform &	pTransform,
		G4double &	pmin,
		G4double &	pmax
	)		const

virtual

Implements G4VSolid.

Definition at line 176 of file G4VCSGfaceted.cc.

References G4VCSGface::CalculateExtent(), faces, G4SolidExtentList::GetExtent(), and numFace.

 {
   G4SolidExtentList  extentList( axis, voxelLimit );
 
   //
   // Loop over all faces, checking min/max extent as we go.
   //
   G4VCSGface **face = faces;
   do
   {
     (*face)->CalculateExtent( axis, voxelLimit, transform, extentList );
   } while( ++face < faces + numFace );
   
   //
   // Return min/max value
   //
   return extentList.GetExtent( min, max );
 }

void G4VCSGfaceted::CopyStuff ( const G4VCSGfaceted & source )

protected

Definition at line 134 of file G4VCSGfaceted.cc.

References G4VCSGface::Clone(), faces, fCubicVolume, fpPolyhedron, fSurfaceArea, and numFace.

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

 {
   numFace = source.numFace;
   if (numFace == 0) { return; }    // odd, but permissable?
   
   faces = new G4VCSGface*[numFace];
   
   G4VCSGface **face = faces,
        **sourceFace = source.faces;
   do
   {
     *face = (*sourceFace)->Clone();
   } while( ++sourceFace, ++face < faces+numFace );
   fCubicVolume = source.fCubicVolume;
   fSurfaceArea = source.fSurfaceArea;
   fpPolyhedron = 0;
 }

virtual G4Polyhedron* G4VCSGfaceted::CreatePolyhedron ( ) const

pure virtual

Reimplemented from G4VSolid.

Implemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Referenced by GetPolyhedron().

void G4VCSGfaceted::DeleteStuff ( )

protected

Definition at line 158 of file G4VCSGfaceted.cc.

References faces, and numFace.

Referenced by operator=(), G4Polycone::Reset(), G4Polyhedra::Reset(), and ~G4VCSGfaceted().

 {
   if (numFace)
   {
     G4VCSGface **face = faces;
     do
     {
       delete *face;
     } while( ++face < faces + numFace );
 
     delete [] faces;
   }
 }

void G4VCSGfaceted::DescribeYourselfTo ( G4VGraphicsScene & scene ) const

virtual

Implements G4VSolid.

Definition at line 403 of file G4VCSGfaceted.cc.

References G4VGraphicsScene::AddSolid().

 {
    scene.AddSolid( *this );
 }

G4double G4VCSGfaceted::DistanceTo	(	const G4ThreeVector &	p,
		const G4bool	outgoing
	)		const

protectedvirtual

Definition at line 385 of file G4VCSGfaceted.cc.

References G4VCSGface::Distance(), faces, G4VSolid::kCarTolerance, and numFace.

Referenced by DistanceToIn(), and DistanceToOut().

 {
   G4VCSGface **face = faces;
   G4double best = kInfinity;
   do
   {
     G4double distance = (*face)->Distance( p, outgoing );
     if (distance < best)  { best = distance; }
   } while( ++face < faces + numFace );
 
   return (best < 0.5*kCarTolerance) ? 0 : best;
 }

G4double G4VCSGfaceted::DistanceToIn	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v
	)		const

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 255 of file G4VCSGfaceted.cc.

References G4VCSGface::Distance(), faces, G4VSolid::kCarTolerance, and numFace.

Referenced by G4GenericPolycone::DistanceToIn(), G4Polyhedra::DistanceToIn(), and G4Polycone::DistanceToIn().

 {
   G4double distance = kInfinity;
   G4double distFromSurface = kInfinity;
   G4VCSGface **face = faces;
   G4VCSGface *bestFace = *face;
   do
   {
     G4double   faceDistance,
                faceDistFromSurface;
     G4ThreeVector   faceNormal;
     G4bool    faceAllBehind;
     if ((*face)->Intersect( p, v, false, kCarTolerance/2,
                 faceDistance, faceDistFromSurface,
                 faceNormal, faceAllBehind ) )
     {
       //
       // Intersecting face
       //
       if (faceDistance < distance)
       {
         distance = faceDistance;
         distFromSurface = faceDistFromSurface;
         bestFace = *face;
         if (distFromSurface <= 0) { return 0; }
       }
     }
   } while( ++face < faces + numFace );
   
   if (distance < kInfinity && distFromSurface<kCarTolerance/2)
   {
     if (bestFace->Distance(p,false) < kCarTolerance/2)  { distance = 0; }
   }
 
   return distance;
 }

G4double G4VCSGfaceted::DistanceToIn ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Reimplemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Definition at line 297 of file G4VCSGfaceted.cc.

References DistanceTo().

 {
   return DistanceTo( p, false );
 }

G4double G4VCSGfaceted::DistanceToOut	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v,
		const G4bool	calcNorm = `false`,
		G4bool *	validNorm = `0`,
		G4ThreeVector *	n = `0`
	)		const

virtual

Implements G4VSolid.

Definition at line 306 of file G4VCSGfaceted.cc.

References G4VCSGface::Distance(), faces, Inside(), G4VSolid::kCarTolerance, kSurface, and numFace.

 {
   G4bool allBehind = true;
   G4double distance = kInfinity;
   G4double distFromSurface = kInfinity;
   G4ThreeVector normal;
   
   G4VCSGface **face = faces;
   G4VCSGface *bestFace = *face;
   do
   {
     G4double  faceDistance,
               faceDistFromSurface;
     G4ThreeVector  faceNormal;
     G4bool    faceAllBehind;
     if ((*face)->Intersect( p, v, true, kCarTolerance/2,
                 faceDistance, faceDistFromSurface,
                 faceNormal, faceAllBehind ) )
     {
       //
       // Intersecting face
       //
       if ( (distance < kInfinity) || (!faceAllBehind) )  { allBehind = false; }
       if (faceDistance < distance)
       {
         distance = faceDistance;
         distFromSurface = faceDistFromSurface;
         normal = faceNormal;
         bestFace = *face;
         if (distFromSurface <= 0)  { break; }
       }
     }
   } while( ++face < faces + numFace );
   
   if (distance < kInfinity)
   {
     if (distFromSurface <= 0)
     {
       distance = 0;
     }
     else if (distFromSurface<kCarTolerance/2)
     {
       if (bestFace->Distance(p,true) < kCarTolerance/2)  { distance = 0; }
     }
 
     if (calcNorm)
     {
       *validNorm = allBehind;
       *n = normal;
     }
   }
   else
   { 
     if (Inside(p) == kSurface)  { distance = 0; }
     if (calcNorm)  { *validNorm = false; }
   }
 
   return distance;
 }

G4double G4VCSGfaceted::DistanceToOut ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 374 of file G4VCSGfaceted.cc.

References DistanceTo().

 {
   return DistanceTo( p, true );
 }

G4double G4VCSGfaceted::GetAreaAccuracy ( ) const

Definition at line 521 of file G4VCSGfaceted.cc.

 {
   return fAreaAccuracy;
 }

G4int G4VCSGfaceted::GetAreaStatistics ( ) const

Definition at line 512 of file G4VCSGfaceted.cc.

 {
   return fStatistics;
 }

G4double G4VCSGfaceted::GetCubicVolume ( )

virtual

Reimplemented from G4VSolid.

Definition at line 550 of file G4VCSGfaceted.cc.

References G4VSolid::EstimateCubicVolume(), and fCubicVolume.

 {
   if(fCubicVolume != 0.) {;}
   else   { fCubicVolume = EstimateCubicVolume(fStatistics,fCubVolEpsilon); }
   return fCubicVolume;
 }

G4double G4VCSGfaceted::GetCubVolEpsilon ( ) const

Definition at line 483 of file G4VCSGfaceted.cc.

 {
   return fCubVolEpsilon;
 }

G4int G4VCSGfaceted::GetCubVolStatistics ( ) const

Definition at line 474 of file G4VCSGfaceted.cc.

 {
   return fStatistics;
 }

G4GeometryType G4VCSGfaceted::GetEntityType ( ) const

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 448 of file G4VCSGfaceted.cc.

 {
   return G4String("G4CSGfaceted");
 }

G4VisExtent G4VCSGfaceted::GetExtent ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 414 of file G4VCSGfaceted.cc.

References faces, and numFace.

 {
   static const G4ThreeVector xMax(1,0,0), xMin(-1,0,0),
                              yMax(0,1,0), yMin(0,-1,0),
                              zMax(0,0,1), zMin(0,0,-1);
   static const G4ThreeVector *axes[6] =
      { &xMin, &xMax, &yMin, &yMax, &zMin, &zMax };
   
   G4double answers[6] =
      {-kInfinity, -kInfinity, -kInfinity, -kInfinity, -kInfinity, -kInfinity};
 
   G4VCSGface **face = faces;
   do
   {    
     const G4ThreeVector **axis = axes+5 ;
     G4double *answer = answers+5;
     do
     {
       G4double testFace = (*face)->Extent( **axis );
       if (testFace > *answer)  { *answer = testFace; }
     }
     while( --axis, --answer >= answers );
     
   } while( ++face < faces + numFace );
   
     return G4VisExtent( -answers[0], answers[1], 
                         -answers[2], answers[3],
                         -answers[4], answers[5]  );
 }

G4ThreeVector G4VCSGfaceted::GetPointOnSurfaceGeneric ( ) const

protected

Definition at line 589 of file G4VCSGfaceted.cc.

References faces, G4UniformRand, numFace, and G4VCSGface::SurfaceArea().

Referenced by G4GenericPolycone::GetPointOnSurface(), and G4Polyhedra::GetPointOnSurface().

 {
   // Preparing variables
   //
   G4ThreeVector answer=G4ThreeVector(0.,0.,0.);
   G4VCSGface **face = faces;
   G4double area = 0;
   G4int i;
   std::vector<G4double> areas; 
 
   // First step: calculate surface areas
   //
   do
   {
     G4double result = (*face)->SurfaceArea( );
     areas.push_back(result);
     area=area+result;
   } while( ++face < faces + numFace );
 
   // Second Step: choose randomly one surface
   //
   G4VCSGface **face1 = faces;
   G4double chose = area*G4UniformRand();
   G4double Achose1, Achose2;
   Achose1=0; Achose2=0.; 
   i=0;
 
   do
   {
     Achose2+=areas[i];
     if(chose>=Achose1 && chose<Achose2)
     {
       G4ThreeVector point;
       point= (*face1)->GetPointOnFace();
       return point;
     }
     i++;
     Achose1=Achose2;
   } while( ++face1 < faces + numFace );
 
   return answer;
 }

G4Polyhedron * G4VCSGfaceted::GetPolyhedron ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 572 of file G4VCSGfaceted.cc.

References CreatePolyhedron(), fpPolyhedron, HepPolyhedron::GetNumberOfRotationSteps(), and G4Polyhedron::GetNumberOfRotationStepsAtTimeOfCreation().

 {
   if (!fpPolyhedron ||
       fpPolyhedron->GetNumberOfRotationStepsAtTimeOfCreation() !=
       fpPolyhedron->GetNumberOfRotationSteps())
   {
     delete fpPolyhedron;
     fpPolyhedron = CreatePolyhedron();
   }
   return fpPolyhedron;
 }

G4double G4VCSGfaceted::GetSurfaceArea ( )

virtual

Reimplemented from G4VSolid.

Definition at line 561 of file G4VCSGfaceted.cc.

References G4VSolid::EstimateSurfaceArea(), and fSurfaceArea.

 {
   if(fSurfaceArea != 0.) {;}
   else   { fSurfaceArea = EstimateSurfaceArea(fStatistics,fAreaAccuracy); }
   return fSurfaceArea;
 }

EInside G4VCSGfaceted::Inside ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 208 of file G4VCSGfaceted.cc.

References faces, G4VCSGface::Inside(), G4VSolid::kCarTolerance, kOutside, kSurface, and numFace.

Referenced by DistanceToOut(), G4GenericPolycone::Inside(), G4Polyhedra::Inside(), and G4Polycone::Inside().

 {
   EInside answer=kOutside;
   G4VCSGface **face = faces;
   G4double best = kInfinity;
   do
   {
     G4double distance;
     EInside result = (*face)->Inside( p, kCarTolerance/2, &distance );
     if (result == kSurface) { return kSurface; }
     if (distance < best)
     {
       best = distance;
       answer = result;
     }
   } while( ++face < faces + numFace );
 
   return answer;
 }

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

Definition at line 108 of file G4VCSGfaceted.cc.

References CopyStuff(), DeleteStuff(), and G4VSolid::operator=().

Referenced by G4GenericPolycone::operator=(), G4Polycone::operator=(), and G4Polyhedra::operator=().

 {
   if (&source == this) { return *this; }
   
   // Copy base class data
   //
   G4VSolid::operator=(source);
 
   // Copy data
   //
   fStatistics = source.fStatistics;
   fCubVolEpsilon = source.fCubVolEpsilon;
   fAreaAccuracy = source.fAreaAccuracy;
 
   DeleteStuff();
   CopyStuff( source );
   
   return *this;
 }

void G4VCSGfaceted::SetAreaAccuracy ( G4double ep )

Definition at line 540 of file G4VCSGfaceted.cc.

References fSurfaceArea.

 {
   fSurfaceArea=0.;
   fAreaAccuracy=ep;
 }

void G4VCSGfaceted::SetAreaStatistics ( G4int st )

Definition at line 530 of file G4VCSGfaceted.cc.

References fSurfaceArea.

 {
   fSurfaceArea=0.;
   fStatistics=st;
 }

void G4VCSGfaceted::SetCubVolEpsilon ( G4double ep )

Definition at line 502 of file G4VCSGfaceted.cc.

References fCubicVolume.

 {
   fCubicVolume=0.;
   fCubVolEpsilon=ep;
 }

void G4VCSGfaceted::SetCubVolStatistics ( G4int st )

Definition at line 492 of file G4VCSGfaceted.cc.

References fCubicVolume.

 {
   fCubicVolume=0.;
   fStatistics=st;
 }

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

virtual

Implements G4VSolid.

Reimplemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Definition at line 457 of file G4VCSGfaceted.cc.

References G4VSolid::GetName(), and numFace.

 {
   os << "-----------------------------------------------------------\n"
      << "    *** Dump for solid - " << GetName() << " ***\n"
      << "    ===================================================\n"
      << " Solid type: G4VCSGfaceted\n"
      << " Parameters: \n"
      << "    number of faces: " << numFace << "\n"
      << "-----------------------------------------------------------\n";
 
   return os;
 }

G4ThreeVector G4VCSGfaceted::SurfaceNormal ( const G4ThreeVector & p ) const

virtual

Implements G4VSolid.

Definition at line 232 of file G4VCSGfaceted.cc.

References faces, and numFace.

 {
   G4ThreeVector answer;
   G4VCSGface **face = faces;
   G4double best = kInfinity;
   do
   {
     G4double distance;
     G4ThreeVector normal = (*face)->Normal( p, &distance );
     if (distance < best)
     {
       best = distance;
       answer = normal;
     }
   } while( ++face < faces + numFace );
 
   return answer;
 }