G4GeomTestVolume Class Reference

#include <G4GeomTestVolume.hh>

Public Member Functions
	G4GeomTestVolume (G4VPhysicalVolume *theTarget, G4double theTolerance=0.0, G4int numberOfPoints=10000, G4bool theVerbosity=true)
	~G4GeomTestVolume ()
G4double	GetTolerance () const
void	SetTolerance (G4double tolerance)
G4int	GetResolution () const
void	SetResolution (G4int points)
G4bool	GetVerbosity () const
void	SetVerbosity (G4bool verbosity)
G4int	GetErrorsThreshold () const
void	SetErrorsThreshold (G4int max)
void	TestOverlapInTree () const
void	TestRecursiveOverlap (G4int sLevel=0, G4int depth=-1)

Detailed Description

Definition at line 43 of file G4GeomTestVolume.hh.

Constructor & Destructor Documentation

G4GeomTestVolume()

G4GeomTestVolume::G4GeomTestVolume	(	G4VPhysicalVolume *	theTarget,
		G4double	theTolerance = 0.0,
		G4int	numberOfPoints = 10000,
		G4bool	theVerbosity = true
	)

Definition at line 43 of file G4GeomTestVolume.cc.

  : target(theTarget), tolerance(theTolerance),
    resolution(numberOfPoints), verbosity(theVerbosity)
{;}

~G4GeomTestVolume()

G4GeomTestVolume::~G4GeomTestVolume

(

)

Definition at line 54 of file G4GeomTestVolume.cc.

{;}

Member Function Documentation

GetErrorsThreshold()

G4int G4GeomTestVolume::GetErrorsThreshold

(

)

const

Definition at line 107 of file G4GeomTestVolume.cc.

{
  return maxErr;
}

GetResolution()

G4int G4GeomTestVolume::GetResolution

(

)

const

Definition at line 75 of file G4GeomTestVolume.cc.

{
  return resolution;
}

GetTolerance()

G4double G4GeomTestVolume::GetTolerance

(

)

const

Definition at line 59 of file G4GeomTestVolume.cc.

{
  return tolerance;
}

GetVerbosity()

G4bool G4GeomTestVolume::GetVerbosity

(

)

const

Definition at line 91 of file G4GeomTestVolume.cc.

{
  return verbosity;
}

SetErrorsThreshold()

void G4GeomTestVolume::SetErrorsThreshold

(

G4int

max

)

Definition at line 115 of file G4GeomTestVolume.cc.

{
  maxErr = max;
}

Referenced by G4GeometryMessenger::SetNewValue(), and TestRecursiveOverlap().

SetResolution()

void G4GeomTestVolume::SetResolution

(

G4int

points

)

Definition at line 83 of file G4GeomTestVolume.cc.

{
  resolution = np;
}

Referenced by G4GeometryMessenger::SetNewValue().

SetTolerance()

void G4GeomTestVolume::SetTolerance

(

G4double

tolerance

)

Definition at line 67 of file G4GeomTestVolume.cc.

{
  tolerance = tol;
}

Referenced by G4GeometryMessenger::SetNewValue().

SetVerbosity()

void G4GeomTestVolume::SetVerbosity

(

G4bool

verbosity

)

Definition at line 99 of file G4GeomTestVolume.cc.

{
  verbosity = verb;
}

Referenced by G4GeometryMessenger::SetNewValue().

TestOverlapInTree()

void G4GeomTestVolume::TestOverlapInTree

(

)

const

Definition at line 123 of file G4GeomTestVolume.cc.

{
  std::queue<G4VPhysicalVolume*> volumes;
  std::set<G4LogicalVolume*> checked;
 
  volumes.push(target);
  while (!volumes.empty())
  {
    G4VPhysicalVolume* current = volumes.front();
    volumes.pop();
 
    // check overlaps for daughters
    G4LogicalVolume* logical = current->GetLogicalVolume();
    G4int ndaughters = logical->GetNoDaughters();
    for (G4int i=0; i<ndaughters; ++i)
    {
      G4VPhysicalVolume* daughter = logical->GetDaughter(i);
      daughter->CheckOverlaps(resolution, tolerance, verbosity, maxErr);
    }
 
    // append the queue of volumes
    G4LogicalVolume* previousLogical = nullptr;
    for (G4int i=0; i<ndaughters; ++i)
    {
      G4VPhysicalVolume* daughter = logical->GetDaughter(i);
      G4LogicalVolume* daughterLogical = daughter->GetLogicalVolume();
      if (daughterLogical->GetNoDaughters() == 0) continue;
      G4bool found = (daughterLogical == previousLogical);
      if (!found) found = (checked.find(daughterLogical) != checked.cend());
      if (!found)
      {
        checked.emplace(daughterLogical);
        previousLogical = daughterLogical;
        volumes.push(daughter);
      }
      else
      {
        if (verbosity)
          G4cout << "Checking overlaps in tree of volume " << daughter->GetName()
                 << " (" << daughterLogical->GetSolid()->GetEntityType() << ")"
                 << " is omitted, to avoid duplication" << G4endl;
      }
    }
  }
}

TestRecursiveOverlap()

void G4GeomTestVolume::TestRecursiveOverlap	(	G4int	sLevel = 0,
		G4int	depth = -1
	)

Definition at line 172 of file G4GeomTestVolume.cc.

{
  // If reached requested level of depth (i.e. set to 0), exit.
  // If not depth specified (i.e. set to -1), visit the whole tree.
  // If requested initial level of depth is not zero, visit from beginning
  //
  if (depth == 0) return;
  if (depth != -1) depth--;
  if (slevel != 0) slevel--;
 
  //
  // As long as we reached the requested
  // initial level of depth, test ourselves
  //
  if ( slevel==0 )
  {
    target->CheckOverlaps(resolution, tolerance, verbosity, maxErr);
  }
 
  //
  // Loop over unique daughters
  //
  std::set<const G4LogicalVolume *> tested;
 
  const G4LogicalVolume *logical = target->GetLogicalVolume();
  G4int nDaughter = logical->GetNoDaughters();
  for( auto iDaughter=0; iDaughter<nDaughter; ++iDaughter )
  {
    G4VPhysicalVolume *daughter = logical->GetDaughter(iDaughter);
    
    // Tested already?
    //
    // const G4LogicalVolume *daughterLogical =
    //      daughter->GetLogicalVolume();
    // std::pair<std::set<const G4LogicalVolume *>::iterator, G4bool>
    //       there = tested.insert(daughterLogical);
    // if (!there.second) continue;
 
    //
    // Recurse
    //
    G4GeomTestVolume vTest( daughter, tolerance, resolution, verbosity );
    vTest.SetErrorsThreshold(maxErr);
    vTest.TestRecursiveOverlap( slevel,depth );
  }
}

Referenced by TestRecursiveOverlap().

---

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

• G4GeomTestVolume.hh
• G4GeomTestVolume.cc