#include <G4TessellatedSolid.hh>

Inheritance diagram for G4TessellatedSolid:

Public Member Functions
	G4TessellatedSolid ()

virtual	~G4TessellatedSolid ()

	G4TessellatedSolid (const G4String &name)

	G4TessellatedSolid (__void__ &)

	G4TessellatedSolid (const G4TessellatedSolid &ts)

G4TessellatedSolid &	operator= (const G4TessellatedSolid &right)

G4TessellatedSolid &	operator+= (const G4TessellatedSolid &right)

G4bool	AddFacet (G4VFacet *aFacet)

G4VFacet *	GetFacet (G4int i) const

G4int	GetNumberOfFacets () 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

virtual G4double	DistanceToOut (const G4ThreeVector &p, const G4ThreeVector &v, const G4bool calcNorm, G4bool validNorm, G4ThreeVector norm) const

virtual G4bool	Normal (const G4ThreeVector &p, G4ThreeVector &n) const

virtual G4double	SafetyFromOutside (const G4ThreeVector &p, G4bool aAccurate=false) const

virtual G4double	SafetyFromInside (const G4ThreeVector &p, G4bool aAccurate=false) const

virtual G4GeometryType	GetEntityType () const

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

virtual G4VSolid *	Clone () const

virtual G4ThreeVector	GetPointOnSurface () const

virtual G4double	GetSurfaceArea ()

virtual G4double	GetCubicVolume ()

void	SetSolidClosed (const G4bool t)

G4bool	GetSolidClosed () const

void	SetMaxVoxels (G4int max)

G4SurfaceVoxelizer &	GetVoxels ()

virtual G4bool	CalculateExtent (const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pMin, G4double &pMax) const

G4double	GetMinXExtent () const

G4double	GetMaxXExtent () const

G4double	GetMinYExtent () const

G4double	GetMaxYExtent () const

G4double	GetMinZExtent () const

G4double	GetMaxZExtent () const

G4ThreeVectorList *	CreateRotatedVertices (const G4AffineTransform &pT) const

virtual G4Polyhedron *	CreatePolyhedron () const

virtual G4Polyhedron *	GetPolyhedron () const

virtual G4NURBS *	CreateNURBS () const

virtual void	DescribeYourselfTo (G4VGraphicsScene &scene) const

virtual G4VisExtent	GetExtent () const

G4int	AllocatedMemoryWithoutVoxels ()

G4int	AllocatedMemory ()

void	DisplayAllocatedMemory ()

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 G4bool	CalculateExtent (const EAxis pAxis, const G4VoxelLimits &pVoxelLimit, const G4AffineTransform &pTransform, G4double &pMin, G4double &pMax) const =0

virtual EInside	Inside (const G4ThreeVector &p) const =0

virtual G4ThreeVector	SurfaceNormal (const G4ThreeVector &p) const =0

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

virtual G4double	DistanceToIn (const G4ThreeVector &p) const =0

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

virtual G4double	DistanceToOut (const G4ThreeVector &p) const =0

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

virtual G4double	GetCubicVolume ()

virtual G4double	GetSurfaceArea ()

virtual G4GeometryType	GetEntityType () const =0

virtual G4ThreeVector	GetPointOnSurface () const

virtual G4VSolid *	Clone () const

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

void	DumpInfo () const

virtual void	DescribeYourselfTo (G4VGraphicsScene &scene) const =0

virtual G4VisExtent	GetExtent () const

virtual G4Polyhedron *	CreatePolyhedron () const

virtual G4NURBS *	CreateNURBS () const

virtual G4Polyhedron *	GetPolyhedron () 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)

Additional Inherited Members
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 inherited from G4VSolid
G4double	kCarTolerance

Detailed Description

Definition at line 128 of file G4TessellatedSolid.hh.

Constructor & Destructor Documentation

◆ G4TessellatedSolid() [1/4]

G4TessellatedSolid::G4TessellatedSolid ( )

Definition at line 101 of file G4TessellatedSolid.cc.

                                        : G4VSolid("dummy")
{
  Initialize();
}

Referenced by Clone().

◆ ~G4TessellatedSolid()

G4TessellatedSolid::~G4TessellatedSolid ( )

virtual

Definition at line 131 of file G4TessellatedSolid.cc.

{
  DeleteObjects ();
}

◆ G4TessellatedSolid() [2/4]

G4TessellatedSolid::G4TessellatedSolid ( const G4String & name )

Definition at line 111 of file G4TessellatedSolid.cc.

  : G4VSolid(name)
{
  Initialize();
}

◆ G4TessellatedSolid() [3/4]

G4TessellatedSolid::G4TessellatedSolid ( __void__ & a )

Definition at line 122 of file G4TessellatedSolid.cc.

                                                   : G4VSolid(a)
{
  Initialize();
  fMinExtent.set(0,0,0);
  fMaxExtent.set(0,0,0);
}

◆ G4TessellatedSolid() [4/4]

G4TessellatedSolid::G4TessellatedSolid ( const G4TessellatedSolid & ts )

Definition at line 140 of file G4TessellatedSolid.cc.

  : G4VSolid(ts), fpPolyhedron(0)
{
  Initialize();
 
  CopyObjects(ts);
}

Member Function Documentation

◆ AddFacet()

G4bool G4TessellatedSolid::AddFacet ( G4VFacet * aFacet )

Definition at line 220 of file G4TessellatedSolid.cc.

{
  // Add the facet to the vector.
  //
  if (fSolidClosed)
  {
    G4Exception("G4TessellatedSolid::AddFacet()", "GeomSolids1002",
                JustWarning, "Attempt to add facets when solid is closed.");
    return false;
  }
  else if (aFacet->IsDefined())
  {
    set<G4VertexInfo,G4VertexComparator>::iterator begin
      = fFacetList.begin(), end = fFacetList.end(), pos, it;
    G4ThreeVector p = aFacet->GetCircumcentre();
    G4VertexInfo value;
    value.id = fFacetList.size();
    value.mag2 = p.mag2();
 
    G4bool found = false;
    if (!OutsideOfExtent(p, kCarTolerance))
    {
      G4double kCarTolerance24 = kCarTolerance * kCarTolerance / 4.0;
      pos = fFacetList.lower_bound(value);
 
      it = pos;
      while (!found && it != end)
      {
        G4int id = (*it).id;
        G4ThreeVector q = fFacets[id]->GetCircumcentre();
        if ((found = (fFacets[id] == aFacet))) break;
        G4double dif = q.mag2() - value.mag2;
        if (dif > kCarTolerance24) break;
        it++;
      }
 
      if (fFacets.size() > 1)
      {
        it = pos;
        while (!found && it != begin)
        {
          --it;
          G4int id = (*it).id;
          G4ThreeVector q = fFacets[id]->GetCircumcentre();
          found = (fFacets[id] == aFacet);
          if (found) break;
          G4double dif = value.mag2 - q.mag2();
          if (dif > kCarTolerance24) break;
        }
      }
    }
 
    if (!found)
    {
      fFacets.push_back(aFacet);
      fFacetList.insert(value);
    }
 
    return true;
  }
  else
  {
    G4Exception("G4TessellatedSolid::AddFacet()", "GeomSolids1002",
                JustWarning, "Attempt to add facet not properly defined.");    
    aFacet->StreamInfo(G4cout);
    return false;
  }
}

Referenced by G4tgbVolume::FindOrConstructG4Solid(), operator+=(), and G4GDMLReadSolids::TessellatedRead().

◆ AllocatedMemory()

G4int G4TessellatedSolid::AllocatedMemory ( )

Definition at line 1991 of file G4TessellatedSolid.cc.

{
  G4int size = AllocatedMemoryWithoutVoxels();
  G4int sizeInsides = fInsides.GetNbytes();
  G4int sizeVoxels = fVoxels.AllocatedMemory();
  size += sizeInsides + sizeVoxels;
  return size;
}

Referenced by DisplayAllocatedMemory().

◆ AllocatedMemoryWithoutVoxels()

G4int G4TessellatedSolid::AllocatedMemoryWithoutVoxels ( )

Definition at line 1965 of file G4TessellatedSolid.cc.

{
  G4int base = sizeof(*this);
  base += fVertexList.capacity() * sizeof(G4ThreeVector);
  base += fRandir.capacity() * sizeof(G4ThreeVector);
 
  G4int limit = fFacets.size();
  for (G4int i = 0; i < limit; i++)
  {
    G4VFacet &facet = *fFacets[i];
    base += facet.AllocatedMemory();
  }
 
  std::set<G4VFacet *>::const_iterator beg, end, it;
  beg = fExtremeFacets.begin();
  end = fExtremeFacets.end();
  for (it = beg; it != end; it++)
  {
    G4VFacet &facet = *(*it);
    base += facet.AllocatedMemory();
  }
  return base;
}

Referenced by AllocatedMemory(), and DisplayAllocatedMemory().

◆ CalculateExtent()

G4bool G4TessellatedSolid::CalculateExtent	(	const EAxis	pAxis,
		const G4VoxelLimits &	pVoxelLimit,
		const G4AffineTransform &	pTransform,
		G4double &	pMin,
		G4double &	pMax
	)		const

virtual

Implements G4VSolid.

Definition at line 1742 of file G4TessellatedSolid.cc.

{
  G4ThreeVectorList transVertexList(fVertexList);
  G4int size = fVertexList.size();
 
  // Put solid into transformed frame
  for (G4int i=0; i < size; ++i)
  {
    pTransform.ApplyPointTransform(transVertexList[i]);
  }
 
  // Find min and max extent in each dimension
  G4ThreeVector minExtent(kInfinity, kInfinity, kInfinity);
  G4ThreeVector maxExtent(-kInfinity, -kInfinity, -kInfinity);
 
  size = transVertexList.size();
  for (G4int i=0; i< size; ++i)
  {
    for (G4int axis=G4ThreeVector::X; axis < G4ThreeVector::SIZE; ++axis)
    {
      G4double coordinate = transVertexList[i][axis];
      if (coordinate < minExtent[axis])
      { minExtent[axis] = coordinate; }
      if (coordinate > maxExtent[axis])
      { maxExtent[axis] = coordinate; }
    }
  }
 
  // Check for containment and clamp to voxel boundaries
  for (G4int axis=G4ThreeVector::X; axis < G4ThreeVector::SIZE; ++axis)
  {
    EAxis geomAxis = kXAxis; // U geom classes use different index type
    switch(axis)
    {
    case G4ThreeVector::X: geomAxis = kXAxis; break;
    case G4ThreeVector::Y: geomAxis = kYAxis; break;
    case G4ThreeVector::Z: geomAxis = kZAxis; break;
    }
    G4bool isLimited = pVoxelLimit.IsLimited(geomAxis);
    G4double voxelMinExtent = pVoxelLimit.GetMinExtent(geomAxis);
    G4double voxelMaxExtent = pVoxelLimit.GetMaxExtent(geomAxis);
 
    if (isLimited)
    {
      if ( minExtent[axis] > voxelMaxExtent+kCarTolerance ||
        maxExtent[axis] < voxelMinExtent-kCarTolerance    )
      {
        return false ;
      }
      else
      {
        if (minExtent[axis] < voxelMinExtent)
        {
          minExtent[axis] = voxelMinExtent ;
        }
        if (maxExtent[axis] > voxelMaxExtent)
        {
          maxExtent[axis] = voxelMaxExtent;
        }
      }
    }
  }
 
  // Convert pAxis into G4ThreeVector index
  G4int vecAxis=0;
  switch(pAxis)
  {
  case kXAxis: vecAxis = G4ThreeVector::X; break;
  case kYAxis: vecAxis = G4ThreeVector::Y; break;
  case kZAxis: vecAxis = G4ThreeVector::Z; break;
  default: break;
  } 
 
  pMin = minExtent[vecAxis] - kCarTolerance;
  pMax = maxExtent[vecAxis] + kCarTolerance;
 
  return true;
}

Referenced by G4GenericTrap::CalculateExtent().

◆ Clone()

G4VSolid * G4TessellatedSolid::Clone ( ) const

virtual

Reimplemented from G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1564 of file G4TessellatedSolid.cc.

{
  return new G4TessellatedSolid(*this);
}

◆ CreateNURBS()

G4NURBS * G4TessellatedSolid::CreateNURBS ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 1714 of file G4TessellatedSolid.cc.

{
  return 0;
}

Referenced by G4GenericTrap::CreateNURBS().

◆ CreatePolyhedron()

G4Polyhedron * G4TessellatedSolid::CreatePolyhedron ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 1680 of file G4TessellatedSolid.cc.

{
  G4int nVertices = fVertexList.size();
  G4int nFacets   = fFacets.size();
  G4PolyhedronArbitrary *polyhedron =
    new G4PolyhedronArbitrary (nVertices, nFacets);
  for (G4ThreeVectorList::const_iterator v= fVertexList.begin();
                                         v!=fVertexList.end(); ++v)
  {
    polyhedron->AddVertex(*v);
  }
 
  G4int size = fFacets.size();
  for (G4int i = 0; i < size; ++i)
  {
    G4VFacet &facet = *fFacets[i];
    G4int v[4];
    int n = facet.GetNumberOfVertices();
    if (n > 4) n = 4;
    else if (n == 3) v[3] = 0;
    for (G4int j=0; j<n; ++j)
    {
      G4int k = facet.GetVertexIndex(j);
      v[j] = k+1;
    }
    polyhedron->AddFacet(v[0],v[1],v[2],v[3]);
  }
  polyhedron->SetReferences();  
 
  return (G4Polyhedron*) polyhedron;
}

Referenced by G4GenericTrap::CreatePolyhedron(), and GetPolyhedron().

◆ CreateRotatedVertices()

G4ThreeVectorList * G4TessellatedSolid::CreateRotatedVertices ( const G4AffineTransform & pT ) const

◆ DescribeYourselfTo()

void G4TessellatedSolid::DescribeYourselfTo ( G4VGraphicsScene & scene ) const

virtual

Implements G4VSolid.

Definition at line 1673 of file G4TessellatedSolid.cc.

{
  scene.AddSolid (*this);
}

Referenced by G4GenericTrap::DescribeYourselfTo().

◆ DisplayAllocatedMemory()

void G4TessellatedSolid::DisplayAllocatedMemory ( )

Definition at line 549 of file G4TessellatedSolid.cc.

{
  G4int without = AllocatedMemoryWithoutVoxels();
  G4int with = AllocatedMemory();
  G4double ratio = (G4double) with / without;
  G4cout << "G4TessellatedSolid - Allocated memory without voxel overhead "
         << without << "; with " << with << "; ratio: " << ratio << endl;
}

Referenced by SetSolidClosed().

◆ DistanceToIn() [1/2]

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

virtual

Implements G4VSolid.

Definition at line 1610 of file G4TessellatedSolid.cc.

{
  return SafetyFromOutside(p,false);
}

◆ DistanceToIn() [2/2]

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

virtual

Implements G4VSolid.

Definition at line 1617 of file G4TessellatedSolid.cc.

{
  return DistanceToInCore(p,v,kInfinity);
}

Referenced by G4GenericTrap::DistanceToIn(), and SafetyFromInside().

◆ DistanceToOut() [1/2]

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

virtual

Implements G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1630 of file G4TessellatedSolid.cc.

{
  return SafetyFromInside(p,false);
}

Referenced by G4ExtrudedSolid::DistanceToOut(), G4GenericTrap::DistanceToOut(), and SafetyFromOutside().

◆ DistanceToOut() [2/2]

G4double G4TessellatedSolid::DistanceToOut	(	const G4ThreeVector &	p,
		const G4ThreeVector &	v,
		const G4bool	calcNorm,
		G4bool *	validNorm,
		G4ThreeVector *	norm
	)		const

virtual

Implements G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1653 of file G4TessellatedSolid.cc.

{
  G4ThreeVector n;
  G4bool valid;
 
  G4double dist = DistanceToOutCore(p, v, n, valid);
  if (calcNorm)
  {
    *norm = n;
    *validNorm = valid;
  }
  return dist;
}

◆ GetCubicVolume()

G4double G4TessellatedSolid::GetCubicVolume ( )

virtual

Reimplemented from G4VSolid.

Definition at line 1875 of file G4TessellatedSolid.cc.

{
  if(fCubicVolume != 0.) {;}
  else   { fCubicVolume = G4VSolid::GetCubicVolume(); }
  return fCubicVolume;
}

◆ GetEntityType()

G4GeometryType G4TessellatedSolid::GetEntityType ( ) const

virtual

Implements G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1535 of file G4TessellatedSolid.cc.

{
  return fGeometryType;
}

◆ GetExtent()

G4VisExtent G4TessellatedSolid::GetExtent ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 1868 of file G4TessellatedSolid.cc.

{
  return G4VisExtent (fMinExtent.x(), fMaxExtent.x(), fMinExtent.y(), fMaxExtent.y(), fMinExtent.z(), fMaxExtent.z());
}

Referenced by G4GenericTrap::GetExtent().

◆ GetFacet()

G4VFacet * G4TessellatedSolid::GetFacet ( G4int i ) const

inline

Definition at line 303 of file G4TessellatedSolid.hh.

{
  return fFacets[i];
}

Referenced by operator+=(), and G4GDMLWriteSolids::TessellatedWrite().

◆ GetMaxXExtent()

G4double G4TessellatedSolid::GetMaxXExtent ( ) const

Definition at line 1833 of file G4TessellatedSolid.cc.

{
  return fMaxExtent.x();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetMaxYExtent()

G4double G4TessellatedSolid::GetMaxYExtent ( ) const

Definition at line 1847 of file G4TessellatedSolid.cc.

{
  return fMaxExtent.y();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetMaxZExtent()

G4double G4TessellatedSolid::GetMaxZExtent ( ) const

Definition at line 1861 of file G4TessellatedSolid.cc.

{
  return fMaxExtent.z();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetMinXExtent()

G4double G4TessellatedSolid::GetMinXExtent ( ) const

Definition at line 1826 of file G4TessellatedSolid.cc.

{
  return fMinExtent.x();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetMinYExtent()

G4double G4TessellatedSolid::GetMinYExtent ( ) const

Definition at line 1840 of file G4TessellatedSolid.cc.

{
  return fMinExtent.y();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetMinZExtent()

G4double G4TessellatedSolid::GetMinZExtent ( ) const

Definition at line 1854 of file G4TessellatedSolid.cc.

{
  return fMinExtent.z();
}

Referenced by G4ExtrudedSolid::Inside().

◆ GetNumberOfFacets()

G4int G4TessellatedSolid::GetNumberOfFacets ( ) const

Definition at line 612 of file G4TessellatedSolid.cc.

{
  return fFacets.size();
}

Referenced by operator+=(), and G4GDMLWriteSolids::TessellatedWrite().

◆ GetPointOnSurface()

G4ThreeVector G4TessellatedSolid::GetPointOnSurface ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 1899 of file G4TessellatedSolid.cc.

{
  // Select randomly a facet and return a random point on it
 
  G4int i = (G4int) G4RandFlat::shoot(0., fFacets.size());
  return fFacets[i]->GetPointOnFace();
}

Referenced by G4GenericTrap::GetPointOnSurface().

◆ GetPolyhedron()

G4Polyhedron * G4TessellatedSolid::GetPolyhedron ( ) const

virtual

Reimplemented from G4VSolid.

Definition at line 1723 of file G4TessellatedSolid.cc.

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

Referenced by G4GenericTrap::GetPolyhedron().

◆ GetSolidClosed()

G4bool G4TessellatedSolid::GetSolidClosed ( ) const

Definition at line 584 of file G4TessellatedSolid.cc.

{
  return fSolidClosed;
}

◆ GetSurfaceArea()

G4double G4TessellatedSolid::GetSurfaceArea ( )

virtual

Reimplemented from G4VSolid.

Definition at line 1884 of file G4TessellatedSolid.cc.

{
  if (fSurfaceArea != 0.) return fSurfaceArea;
 
  G4int size = fFacets.size();
  for (G4int i = 0; i < size; ++i)
  {
    G4VFacet &facet = *fFacets[i];
    fSurfaceArea += facet.GetArea();
  }
  return fSurfaceArea;
}

◆ GetVoxels()

G4SurfaceVoxelizer & G4TessellatedSolid::GetVoxels ( )

inline

Definition at line 313 of file G4TessellatedSolid.hh.

{
  return fVoxels;
}

◆ Inside()

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

virtual

Implements G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1579 of file G4TessellatedSolid.cc.

{
  EInside location;
 
  if (fVoxels.GetCountOfVoxels() > 1)
  {
    location = InsideVoxels(aPoint);
  }
  else
  {
    location = InsideNoVoxels(aPoint);
  }
  return location;
}

Referenced by G4GenericTrap::Inside(), SafetyFromInside(), and SafetyFromOutside().

◆ Normal()

G4bool G4TessellatedSolid::Normal	(	const G4ThreeVector &	p,
		G4ThreeVector &	n
	)		const

virtual

Definition at line 956 of file G4TessellatedSolid.cc.

{
  G4double minDist;
  G4VFacet *facet = 0;
 
  if (fVoxels.GetCountOfVoxels() > 1)
  {
    vector<G4int> curVoxel(3);
    fVoxels.GetVoxel(curVoxel, p);
    const vector<G4int> &candidates = fVoxels.GetCandidates(curVoxel);
    // fVoxels.GetCandidatesVoxelArray(p, candidates, 0);
 
    if (G4int limit = candidates.size())
    {
      minDist = kInfinity;
      for(G4int i = 0 ; i < limit ; ++i)
      {      
        G4int candidate = candidates[i];
        G4VFacet &fct = *fFacets[candidate];
        G4double dist = fct.Distance(p,minDist);
        if (dist < minDist) minDist = dist;
        if (dist <= kCarToleranceHalf)
        {
          aNormal = fct.GetSurfaceNormal();
          return true;
        }
      }
    }
    minDist = MinDistanceFacet(p, true, facet);
  }
  else
  {
    minDist = kInfinity;
    G4int size = fFacets.size();
    for (G4int i = 0; i < size; ++i)
    {
      G4VFacet &f = *fFacets[i];
      G4double dist = f.Distance(p, minDist);
      if (dist < minDist)
      {
        minDist  = dist;
        facet = &f;
      }
    }
  }
 
  if (minDist != kInfinity)
  {
    if (facet)  { aNormal = facet->GetSurfaceNormal(); }
    return minDist <= kCarToleranceHalf;
  }
  else
  {
#ifdef G4VERBOSE
    std::ostringstream message;
    message << "Point p is not on surface !?" << endl
      << "          No facets found for point: " << p << " !" << endl
      << "          Returning approximated value for normal.";
 
    G4Exception("G4TessellatedSolid::SurfaceNormal(p)",
                "GeomSolids1002", JustWarning, message );
#endif
    aNormal = (p.z() > 0 ? G4ThreeVector(0,0,1) : G4ThreeVector(0,0,-1));
    return false;
  }
}

Referenced by SurfaceNormal().

◆ operator+=()

G4TessellatedSolid & G4TessellatedSolid::operator+= ( const G4TessellatedSolid & right )

Definition at line 599 of file G4TessellatedSolid.cc.

{
  G4int size = right.GetNumberOfFacets();
  for (G4int i = 0; i < size; ++i)
    AddFacet(right.GetFacet(i)->GetClone());
 
  return *this;
}

◆ operator=()

G4TessellatedSolid & G4TessellatedSolid::operator= ( const G4TessellatedSolid & right )

Definition at line 153 of file G4TessellatedSolid.cc.

{
  if (&ts == this) return *this;
 
  // Copy base class data
  G4VSolid::operator=(ts);
 
  DeleteObjects ();
 
  Initialize();
 
  CopyObjects (ts);
 
  return *this;
}

Referenced by G4ExtrudedSolid::operator=().

◆ SafetyFromInside()

G4double G4TessellatedSolid::SafetyFromInside	(	const G4ThreeVector &	p,
		G4bool	aAccurate = `false`
	)		const

virtual

Definition at line 1486 of file G4TessellatedSolid.cc.

{  
#if G4SPECSDEBUG
  if ( Inside(p) == kOutside )
  {
    std::ostringstream message;
    G4int oldprc = message.precision(16) ;
    message << "Point p is already outside!?" << endl
      << "Position:"  << endl << endl
      << "p.x() = "   << p.x()/mm << " mm" << endl
      << "p.y() = "   << p.y()/mm << " mm" << endl
      << "p.z() = "   << p.z()/mm << " mm" << endl
      << "DistanceToIn(p) == " << DistanceToIn(p);
    message.precision(oldprc) ;
    G4Exception("G4TriangularFacet::DistanceToOut(p)",
                "GeomSolids1002", JustWarning, message);
  }
#endif
 
  G4double minDist;
 
  if (OutsideOfExtent(p, kCarTolerance)) return 0.0;
 
  if (fVoxels.GetCountOfVoxels() > 1)
  {
    G4VFacet *facet;
    minDist = MinDistanceFacet(p, true, facet);
  }
  else
  {
    minDist = kInfinity;
    G4double dist = 0.0;
    G4int size = fFacets.size();
    for (G4int i = 0; i < size; ++i)
    {
      G4VFacet &facet = *fFacets[i];
      dist = facet.Distance(p,minDist);
      if (dist < minDist) minDist  = dist;
    }
  }
  return minDist;
}

Referenced by DistanceToOut().

◆ SafetyFromOutside()

G4double G4TessellatedSolid::SafetyFromOutside	(	const G4ThreeVector &	p,
		G4bool	aAccurate = `false`
	)		const

virtual

Definition at line 1427 of file G4TessellatedSolid.cc.

{
#if G4SPECSDEBUG
  if ( Inside(p) == kInside )
  {
    std::ostringstream message;
    G4int oldprc = message.precision(16) ;
    message << "Point p is already inside!?" << endl
      << "Position:"  << endl << endl
      << "p.x() = "   << p.x()/mm << " mm" << endl
      << "p.y() = "   << p.y()/mm << " mm" << endl
      << "p.z() = "   << p.z()/mm << " mm" << endl
      << "DistanceToOut(p) == " << DistanceToOut(p);
    message.precision(oldprc) ;
    G4Exception("G4TriangularFacet::DistanceToIn(p)",
                "GeomSolids1002", JustWarning, message);
  }
#endif
 
  G4double minDist;
 
  if (!aAccurate)
    return fVoxels.DistanceToBoundingBox(p);
 
  if (fVoxels.GetCountOfVoxels() > 1)
  {
    if (!OutsideOfExtent(p, kCarTolerance))
    {
      vector<G4int> startingVoxel(3);
      fVoxels.GetVoxel(startingVoxel, p);
      const vector<G4int> &candidates = fVoxels.GetCandidates(startingVoxel);
      if (candidates.size() == 0 && fInsides.GetNbits())
      {
        G4int index = fVoxels.GetPointIndex(p);
        if (fInsides[index]) return 0.;
      }
    }
 
    G4VFacet *facet;
    minDist = MinDistanceFacet(p, true, facet);
  }
  else
  {
    minDist = kInfinity;
    G4int size = fFacets.size();
    for (G4int i = 0; i < size; ++i)
    {
      G4VFacet &facet = *fFacets[i];
      G4double dist = facet.Distance(p,minDist);
      if (dist < minDist) minDist  = dist;
    }
  }
  return minDist;
}

Referenced by DistanceToIn().

◆ SetMaxVoxels()

void G4TessellatedSolid::SetMaxVoxels ( G4int max )

inline

Definition at line 308 of file G4TessellatedSolid.hh.

{
  fVoxels.SetMaxVoxels(max);
}

◆ SetSolidClosed()

void G4TessellatedSolid::SetSolidClosed ( const G4bool t )

Definition at line 560 of file G4TessellatedSolid.cc.

{
  if (t)
  {
    CreateVertexList();
 
    SetExtremeFacets();
 
    Voxelize();
 
#ifdef G4SPECSDEBUG
    DisplayAllocatedMemory();
#endif
 
  }  
  fSolidClosed = t;
}

Referenced by G4GDMLReadSolids::TessellatedRead().

◆ StreamInfo()

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

virtual

Implements G4VSolid.

Reimplemented in G4ExtrudedSolid.

Definition at line 1542 of file G4TessellatedSolid.cc.

{
  os << endl;
  os << "Geometry Type    = " << fGeometryType  << endl;
  os << "Number of facets = " << fFacets.size() << endl;
 
  G4int size = fFacets.size();
  for (G4int i = 0; i < size; ++i)
  {
    os << "FACET #          = " << i + 1 << endl;
    G4VFacet &facet = *fFacets[i];
    facet.StreamInfo(os);
  }
  os << endl;
 
  return os;
}

◆ SurfaceNormal()

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

virtual

Implements G4VSolid.

Definition at line 1596 of file G4TessellatedSolid.cc.

{
  G4ThreeVector n;
  Normal(p, n);
  return n;
}

Referenced by G4GenericTrap::SurfaceNormal().

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4TessellatedSolid() [1/4]

◆ ~G4TessellatedSolid()

◆ G4TessellatedSolid() [2/4]

◆ G4TessellatedSolid() [3/4]

◆ G4TessellatedSolid() [4/4]

Member Function Documentation

◆ AddFacet()

◆ AllocatedMemory()

◆ AllocatedMemoryWithoutVoxels()

◆ CalculateExtent()

◆ Clone()

◆ CreateNURBS()

◆ CreatePolyhedron()

◆ CreateRotatedVertices()

◆ DescribeYourselfTo()

◆ DisplayAllocatedMemory()

◆ DistanceToIn() [1/2]

◆ DistanceToIn() [2/2]

◆ DistanceToOut() [1/2]

◆ DistanceToOut() [2/2]

◆ GetCubicVolume()

◆ GetEntityType()

◆ GetExtent()

◆ GetFacet()

◆ GetMaxXExtent()

◆ GetMaxYExtent()

◆ GetMaxZExtent()

◆ GetMinXExtent()

◆ GetMinYExtent()

◆ GetMinZExtent()

◆ GetNumberOfFacets()

◆ GetPointOnSurface()

◆ GetPolyhedron()

◆ GetSolidClosed()

◆ GetSurfaceArea()

◆ GetVoxels()

◆ Inside()

◆ Normal()

◆ operator+=()

◆ operator=()

◆ SafetyFromInside()

◆ SafetyFromOutside()

◆ SetMaxVoxels()

◆ SetSolidClosed()

◆ StreamInfo()

◆ SurfaceNormal()