#include <G4VCSGfaceted.hh>

Inheritance diagram for G4VCSGfaceted:

Collaboration diagram for G4VCSGfaceted:

Public Member Functions
	G4VCSGfaceted (const G4String &name)

virtual	~G4VCSGfaceted ()

	G4VCSGfaceted (const G4VCSGfaceted &source)

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	Extent (G4ThreeVector &pMin, G4ThreeVector &pMax) 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)

G4double	EstimateCubicVolume (G4int nStat, G4double epsilon) const

G4double	EstimateSurfaceArea (G4int nStat, G4double ell) const

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

Protected Attributes
G4int	numFace

G4VCSGface **	faces

G4double	fCubicVolume

G4double	fSurfaceArea

G4bool	fRebuildPolyhedron

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 69 of file G4VCSGfaceted.cc.

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

G4VCSGfaceted::~G4VCSGfaceted ( )

virtual

Definition at line 93 of file G4VCSGfaceted.cc.

 {
   DeleteStuff();
   delete fpPolyhedron; fpPolyhedron = 0;
 }

Here is the call graph for this function:

G4VCSGfaceted::G4VCSGfaceted ( const G4VCSGfaceted & source )

Definition at line 103 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

G4VCSGfaceted::G4VCSGfaceted ( __void__ & a )

Definition at line 82 of file G4VCSGfaceted.cc.

   : G4VSolid(a),
     numFace(0), faces(0), fCubicVolume(0.), fSurfaceArea(0.),
     fRebuildPolyhedron(false), 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.

Reimplemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Definition at line 187 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

void G4VCSGfaceted::CopyStuff ( const G4VCSGfaceted & source )

protected

Definition at line 143 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4Polyhedron* G4VCSGfaceted::CreatePolyhedron ( ) const

pure virtual

Reimplemented from G4VSolid.

Implemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Here is the caller graph for this function:

void G4VCSGfaceted::DeleteStuff ( )

protected

Definition at line 168 of file G4VCSGfaceted.cc.

 {
   if (numFace)
   {
     G4VCSGface **face = faces;
     do    // Loop checking, 13.08.2015, G.Cosmo
     {
       delete *face;
     } while( ++face < faces + numFace );
 
     delete [] faces;
   }
   delete fpPolyhedron; fpPolyhedron = 0;
 }

Here is the caller graph for this function:

void G4VCSGfaceted::DescribeYourselfTo ( G4VGraphicsScene & scene ) const

virtual

Implements G4VSolid.

Definition at line 414 of file G4VCSGfaceted.cc.

 {
    scene.AddSolid( *this );
 }

Here is the call graph for this function:

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

protectedvirtual

Definition at line 396 of file G4VCSGfaceted.cc.

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

Here is the caller graph for this function:

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

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 266 of file G4VCSGfaceted.cc.

 {
   G4double distance = kInfinity;
   G4double distFromSurface = kInfinity;
   G4VCSGface **face = faces;
   G4VCSGface *bestFace = *face;
   do    // Loop checking, 13.08.2015, G.Cosmo
   {
     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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Implements G4VSolid.

Reimplemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Definition at line 308 of file G4VCSGfaceted.cc.

 {
   return DistanceTo( p, false );
 }

Here is the call graph for this function:

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 317 of file G4VCSGfaceted.cc.

 {
   G4bool allBehind = true;
   G4double distance = kInfinity;
   G4double distFromSurface = kInfinity;
   G4ThreeVector normal;
   
   G4VCSGface **face = faces;
   G4VCSGface *bestFace = *face;
   do    // Loop checking, 13.08.2015, G.Cosmo
   {
     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;
 }

Here is the call graph for this function:

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

virtual

Implements G4VSolid.

Definition at line 385 of file G4VCSGfaceted.cc.

 {
   return DistanceTo( p, true );
 }

Here is the call graph for this function:

G4double G4VCSGfaceted::GetAreaAccuracy ( ) const

Definition at line 532 of file G4VCSGfaceted.cc.

 {
   return fAreaAccuracy;
 }

G4int G4VCSGfaceted::GetAreaStatistics ( ) const

Definition at line 523 of file G4VCSGfaceted.cc.

 {
   return fStatistics;
 }

G4double G4VCSGfaceted::GetCubicVolume ( )

virtual

Reimplemented from G4VSolid.

Definition at line 561 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

G4double G4VCSGfaceted::GetCubVolEpsilon ( ) const

Definition at line 494 of file G4VCSGfaceted.cc.

 {
   return fCubVolEpsilon;
 }

G4int G4VCSGfaceted::GetCubVolStatistics ( ) const

Definition at line 485 of file G4VCSGfaceted.cc.

 {
   return fStatistics;
 }

G4GeometryType G4VCSGfaceted::GetEntityType ( ) const

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 459 of file G4VCSGfaceted.cc.

 {
   return G4String("G4CSGfaceted");
 }

G4VisExtent G4VCSGfaceted::GetExtent ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 425 of file G4VCSGfaceted.cc.

 {
   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    // Loop checking, 13.08.2015, G.Cosmo
   {    
     const G4ThreeVector **axis = axes+5 ;
     G4double *answer = answers+5;
     do    // Loop checking, 13.08.2015, G.Cosmo
     {
       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 604 of file G4VCSGfaceted.cc.

 {
   // 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    // Loop checking, 13.08.2015, G.Cosmo
   {
     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;
 }

Here is the caller graph for this function:

G4Polyhedron * G4VCSGfaceted::GetPolyhedron ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 583 of file G4VCSGfaceted.cc.

 {
   if (!fpPolyhedron ||
       fRebuildPolyhedron ||
       fpPolyhedron->GetNumberOfRotationStepsAtTimeOfCreation() !=
       fpPolyhedron->GetNumberOfRotationSteps())
   {
     G4AutoLock l(&polyhedronMutex);
     delete fpPolyhedron;
     fpPolyhedron = CreatePolyhedron();
     fRebuildPolyhedron = false;
     l.unlock();
   }
   return fpPolyhedron;
 }

Here is the call graph for this function:

G4double G4VCSGfaceted::GetSurfaceArea ( )

virtual

Reimplemented from G4VSolid.

Definition at line 572 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

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

virtual

Implements G4VSolid.

Reimplemented in G4Polycone, G4Polyhedra, and G4GenericPolycone.

Definition at line 219 of file G4VCSGfaceted.cc.

 {
   EInside answer=kOutside;
   G4VCSGface **face = faces;
   G4double best = kInfinity;
   do    // Loop checking, 13.08.2015, G.Cosmo
   {
     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;
 }

Here is the caller graph for this function:

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

Definition at line 117 of file G4VCSGfaceted.cc.

 {
   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;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4VCSGfaceted::SetAreaAccuracy ( G4double ep )

Definition at line 551 of file G4VCSGfaceted.cc.

 {
   fSurfaceArea=0.;
   fAreaAccuracy=ep;
 }

void G4VCSGfaceted::SetAreaStatistics ( G4int st )

Definition at line 541 of file G4VCSGfaceted.cc.

 {
   fSurfaceArea=0.;
   fStatistics=st;
 }

void G4VCSGfaceted::SetCubVolEpsilon ( G4double ep )

Definition at line 513 of file G4VCSGfaceted.cc.

 {
   fCubicVolume=0.;
   fCubVolEpsilon=ep;
 }

void G4VCSGfaceted::SetCubVolStatistics ( G4int st )

Definition at line 503 of file G4VCSGfaceted.cc.

 {
   fCubicVolume=0.;
   fStatistics=st;
 }

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

virtual

Implements G4VSolid.

Reimplemented in G4Polyhedra, G4Polycone, and G4GenericPolycone.

Definition at line 468 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

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

virtual

Implements G4VSolid.

Definition at line 243 of file G4VCSGfaceted.cc.

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

Here is the call graph for this function:

Member Data Documentation

G4VCSGface** G4VCSGfaceted::faces

protected

Definition at line 120 of file G4VCSGfaceted.hh.

G4double G4VCSGfaceted::fCubicVolume

protected

Definition at line 121 of file G4VCSGfaceted.hh.

G4Polyhedron* G4VCSGfaceted::fpPolyhedron

mutableprotected

Definition at line 124 of file G4VCSGfaceted.hh.

G4bool G4VCSGfaceted::fRebuildPolyhedron

mutableprotected

Definition at line 123 of file G4VCSGfaceted.hh.

G4double G4VCSGfaceted::fSurfaceArea

protected

Definition at line 122 of file G4VCSGfaceted.hh.

G4int G4VCSGfaceted::numFace

protected

Definition at line 119 of file G4VCSGfaceted.hh.

The documentation for this class was generated from the following files:

source/geant4.10.03.p03/source/geometry/solids/specific/include/G4VCSGfaceted.hh
source/geant4.10.03.p03/source/geometry/solids/specific/src/G4VCSGfaceted.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation