G4GeometryWorkspace Class Reference

#include <G4GeometryWorkspace.hh>

Public Types
using	pool_type = G4TWorkspacePool< G4GeometryWorkspace >

Public Member Functions
	G4GeometryWorkspace ()
	~G4GeometryWorkspace ()
void	UseWorkspace ()
void	ReleaseWorkspace ()
void	DestroyWorkspace ()
void	InitialiseWorkspace ()
void	SetVerbose (G4bool v)
G4bool	GetVerbose ()

Static Public Member Functions
static pool_type *	GetPool ()

Protected Member Functions
void	InitialisePhysicalVolumes ()
G4bool	CloneParameterisedSolids (G4PVParameterised *paramVol)
G4bool	CloneReplicaSolid (G4PVReplica *)

Detailed Description

Definition at line 58 of file G4GeometryWorkspace.hh.

Member Typedef Documentation

pool_type

using G4GeometryWorkspace::pool_type = G4TWorkspacePool<G4GeometryWorkspace>

Definition at line 62 of file G4GeometryWorkspace.hh.

Constructor & Destructor Documentation

G4GeometryWorkspace()

G4GeometryWorkspace::G4GeometryWorkspace

(

)

Definition at line 60 of file G4GeometryWorkspace.cc.

{
  fpLogicalVolumeSIM=
      &const_cast<G4LVManager&>(G4LogicalVolume::GetSubInstanceManager());
  fpPhysicalVolumeSIM=
      &const_cast<G4PVManager&>( G4VPhysicalVolume::GetSubInstanceManager() );
  fpReplicaSIM=
      &const_cast<G4PVRManager&>(G4PVReplica::GetSubInstanceManager());
  fpRegionSIM=
      &const_cast<G4RegionManager&>(G4Region::GetSubInstanceManager());
 
  // Create a work area for Logical Volumes in this thread
  // then capture its address
  InitialiseWorkspace();
  
  fLogicalVolumeOffset = fpLogicalVolumeSIM->GetOffset();
 
  fPhysicalVolumeOffset = fpPhysicalVolumeSIM->GetOffset();
 
  fReplicaOffset = fpReplicaSIM->GetOffset();
 
  fRegionOffset = fpRegionSIM->GetOffset();
}

~G4GeometryWorkspace()

G4GeometryWorkspace::~G4GeometryWorkspace

(

)

Definition at line 86 of file G4GeometryWorkspace.cc.

{  
}

Member Function Documentation

CloneParameterisedSolids()

G4bool G4GeometryWorkspace::CloneParameterisedSolids ( G4PVParameterised *  paramVol )

protected

Definition at line 231 of file G4GeometryWorkspace.cc.

{
  // Check whether it is a simple parameterisation or not
  //
  // G4VPVParameterisation *param= paramVol->GetParameterisation();
  // unsigned int numCopies= paramVol->GetMultiplicity();
  // unsigned int numDifferent= 0;
 
  G4LogicalVolume* logicalV= paramVol->GetLogicalVolume();
  G4VSolid* solid= logicalV->GetSolid();
  
  //  for( unsigned int i=0; i< numCopies; ++i)
  //  {
  //    G4VSolid *solidChk= param->ComputeSolid(i, paramVol);
  //    if( solidChk != solid)
  //    {
  //      ++numDifferent;
  //    }
  //  }
  //  if( numDifferent>0 )
  //  {
  //    G4ExceptionDescription ed;
  //    ed << "ERROR - Parameterisation using several instances of Solids \n"
  //       << "potentially to support different types of solids. \n"
  //       << "Geant4-MT currently does not support this type of \n"
  //       << "parameterisation, sorry !";
  //    G4Exception("G4GeometryWorkspace::CloneParameterisedVolume()",
  //                "GeomVol0001", FatalException, ed);
  //  }
  
  // Threads may attempt to clone a solids simultaneously.
  // Those cloned solids will be registered into a shared solid
  // store (C++ container). Need a lock to guarantee thread safety
  //
  G4AutoLock aLock(&solidclone);
  G4VSolid *workerSolid = solid->Clone();
  aLock.unlock();
  if( workerSolid != nullptr )
  {
    logicalV->InitialiseWorker(logicalV,workerSolid,0);
  }
  else
  {
    // In the case that not all solids support(ed) the Clone()
    // method, we do similar thing here to dynamically cast
    // and then get the clone method
    //
    G4ExceptionDescription ed;
    ed << "ERROR - Unable to initialise geometry for worker node. \n"
       << "A solid lacks the Clone() method - or Clone() failed. \n"
       << "   Type of solid: " << solid->GetEntityType() << "\n"
       << "   Parameters: " << *solid;
    G4Exception("G4GeometryWorkspace::CloneParameterisedVolume()",
                "GeomVol0003", FatalException, ed);
  }
  return true; // It Worked
}

Referenced by InitialisePhysicalVolumes().

CloneReplicaSolid()

G4bool G4GeometryWorkspace::CloneReplicaSolid ( G4PVReplica *  replicaPV )

protected

Definition at line 189 of file G4GeometryWorkspace.cc.

{
  // The solid Ptr is in the Logical Volume
  //
  G4LogicalVolume* logicalV = replicaPV->GetLogicalVolume();
  G4VSolid* solid = logicalV->GetSolid();
 
  G4AutoLock aLock(&solidclone);
  G4VSolid* workerSolid = solid->Clone();
  aLock.unlock();
 
  if( workerSolid != nullptr )
  {
    logicalV->InitialiseWorker(logicalV,workerSolid,0);
  }
  else
  {
    // In the case that not all solids support(ed) the Clone()
    // method, we do similar thing here to dynamically cast
    // and then get the clone method.
    //
    G4ExceptionDescription ed;
    ed << "ERROR - Unable to initialise geometry for worker node." << "\n"
       << "A solid lacks the Clone() method - or Clone() failed." << "\n"
       << "   Type of solid: " << solid->GetEntityType() << "\n"
       << "   Parameters: " << *solid;
    G4Exception("G4GeometryWorkspace::CloneParameterisedVolume()",
                "GeomVol0003", FatalException, ed);
    return false; 
  }
  return true; // It Worked
}

Referenced by InitialisePhysicalVolumes().

DestroyWorkspace()

void G4GeometryWorkspace::DestroyWorkspace

(

)

Definition at line 323 of file G4GeometryWorkspace.cc.

{
  G4PhysicalVolumeStore* physVolStore = G4PhysicalVolumeStore::GetInstance();
  for (size_t ip=0; ip<physVolStore->size(); ++ip)
  {
    G4VPhysicalVolume* physVol = (*physVolStore)[ip];
    G4LogicalVolume* logicalVol = physVol->GetLogicalVolume();
    G4PVReplica* g4PVReplica = nullptr;
    g4PVReplica =  dynamic_cast<G4PVReplica*>(physVol);
    if (g4PVReplica != nullptr)
    {
      g4PVReplica->TerminateWorker(g4PVReplica);
      G4PVParameterised* paramVol = nullptr;
      paramVol = dynamic_cast<G4PVParameterised*>(physVol);
      if (paramVol != nullptr)
      {
        // G4VSolid* solid = logicalVol->fSolid;
        logicalVol->TerminateWorker(logicalVol);
        // if( solid->IsClone() ) delete solid;
      }
      else
      {
        logicalVol->TerminateWorker(logicalVol);
      }
    }
    else
    {
      logicalVol->TerminateWorker(logicalVol);
    }
  }
  fpLogicalVolumeSIM->FreeSlave();
  fpPhysicalVolumeSIM->FreeSlave();
  fpReplicaSIM->FreeSlave();
  fpRegionSIM->FreeSlave();
}

Referenced by G4WorkerThread::UpdateGeometryAndPhysicsVectorFromMaster().

GetPool()

G4GeometryWorkspace::pool_type * G4GeometryWorkspace::GetPool ( )

static

Definition at line 53 of file G4GeometryWorkspace.cc.

{
  return &thePool;
}

Referenced by G4WorkerThread::BuildGeometryAndPhysicsVector(), G4WorkerThread::DestroyGeometryAndPhysicsVector(), and G4WorkerThread::UpdateGeometryAndPhysicsVectorFromMaster().

GetVerbose()

G4bool G4GeometryWorkspace::GetVerbose ( )

inline

Definition at line 75 of file G4GeometryWorkspace.hh.

{ return fVerbose; } 

InitialisePhysicalVolumes()

void G4GeometryWorkspace::InitialisePhysicalVolumes ( )

protected

Definition at line 136 of file G4GeometryWorkspace.cc.

{
  G4PhysicalVolumeStore* physVolStore = G4PhysicalVolumeStore::GetInstance();
  for (size_t ip=0; ip<physVolStore->size(); ++ip)
  {
    G4VPhysicalVolume* physVol = (*physVolStore)[ip];
    G4LogicalVolume *logicalVol = physVol->GetLogicalVolume();
 
    // Use shadow pointer
    //
    G4VSolid* solid = logicalVol->GetMasterSolid();
    G4PVReplica* g4PVReplica = nullptr;
    g4PVReplica = dynamic_cast<G4PVReplica*>(physVol);
    if (g4PVReplica == nullptr)
    {
      // Placement volume
      logicalVol->InitialiseWorker(logicalVol,solid,0);
    }
    else
    {          
      g4PVReplica->InitialiseWorker(g4PVReplica);
      if( !g4PVReplica->IsParameterised() )
      {
        logicalVol->InitialiseWorker(logicalVol,solid,0);
 
        // If the replica's solid (in LV) is changed during navigation,
        // it must be thread-private
        //
        CloneReplicaSolid( g4PVReplica );
      }
      else
      {
        G4PVParameterised *paramVol = dynamic_cast<G4PVParameterised*>(physVol);
        if (paramVol == nullptr)
        {
          G4Exception("G4GeometryWorkspace::CreateAndUseWorkspace()",
                      "GeomVol0003", FatalException,
                      "Cannot find Parameterisation for parameterised volume.");
        }
        CloneParameterisedSolids( paramVol );
      }
    }
  }
  if( fVerbose )
  {
    G4cout << "G4GeometryWorkspace::InitialisePhysicalVolumes: "
           << "Copying geometry - Done!" << G4endl;
  }
}

Referenced by InitialiseWorkspace().

InitialiseWorkspace()

void G4GeometryWorkspace::InitialiseWorkspace

(

)

Definition at line 292 of file G4GeometryWorkspace.cc.

{
  if( fVerbose )
  {
    G4cout << "G4GeometryWorkspace::InitialiseWorkspace():"
           << " Copying geometry - Start " << G4endl;
  }
  
  // Implementation originally in:
  //   G4WorkerThread::BuildGeometryAndPhysicsVector()
  // and improved for G4PVParamaterised
  
  // Geometry related, split classes mechanism:
  // Do *NOT* instantiate sub-instance for this thread, just copy the contents!
  //
  fpLogicalVolumeSIM->SlaveCopySubInstanceArray();
  fpPhysicalVolumeSIM->SlaveCopySubInstanceArray();
  fpReplicaSIM->SlaveCopySubInstanceArray();
  fpRegionSIM->SlaveInitializeSubInstance();
 
  InitialisePhysicalVolumes();
  
  if( fVerbose )
  {
    G4cout << "G4GeometryWorkspace::InitialiseWorkspace: "
           << "Copying geometry - Done!" << G4endl;
  }
}

Referenced by G4GeometryWorkspace(), and G4WorkerThread::UpdateGeometryAndPhysicsVectorFromMaster().

ReleaseWorkspace()

void G4GeometryWorkspace::ReleaseWorkspace

(

)

Definition at line 125 of file G4GeometryWorkspace.cc.

{
  fpLogicalVolumeSIM->UseWorkArea(0);
  fpPhysicalVolumeSIM->UseWorkArea(0);
  
  fpReplicaSIM->UseWorkArea(0);
  fpRegionSIM->UseWorkArea(0);
}

SetVerbose()

void G4GeometryWorkspace::SetVerbose ( G4bool  v )

inline

Definition at line 74 of file G4GeometryWorkspace.hh.

{ fVerbose = v; } 

UseWorkspace()

void G4GeometryWorkspace::UseWorkspace

(

)

Definition at line 93 of file G4GeometryWorkspace.cc.

{
  if( fVerbose ) 
  { 
    G4cout << "G4GeometryWorkspace::UseWorkspace: Start " << G4endl;
  }
 
  // Implementation originally in:
  //   G4WorkerThread::BuildGeometryAndPhysicsVector()
  // and improved for G4PVParamaterised
  
  // Geometry related, split classes mechanism: instantiate sub-instance
  // for this thread
  fpLogicalVolumeSIM->UseWorkArea(fLogicalVolumeOffset);
  fpPhysicalVolumeSIM->UseWorkArea(fPhysicalVolumeOffset);
  
  fpReplicaSIM->UseWorkArea(fReplicaOffset);
  fpRegionSIM->UseWorkArea(fRegionOffset);
 
  // When recycling a workspace 
  //   - it must be a lightweight operation, to reuse a valid work area
  //   - so it must NOT Initialise anything!
  // Do not call InitialisePhysicalVolumes();
 
  if( fVerbose )
  { 
     G4cout << "G4GeometryWorkspace::UseWorkspace:  End " << G4endl;
  }
}

---

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

• G4GeometryWorkspace.hh
• G4GeometryWorkspace.cc