G4PhysicalVolumeModel Class Reference

#include <G4PhysicalVolumeModel.hh>

Inheritance diagram for G4PhysicalVolumeModel:

Classes
class	G4PhysicalVolumeNodeID
class	G4PhysicalVolumeModelTouchable
struct	TouchableProperties

Public Types
enum	{ UNLIMITED = -1 }
enum	ClippingMode { subtraction , intersection }

Public Member Functions
	G4PhysicalVolumeModel (G4VPhysicalVolume *=0, G4int requestedDepth=UNLIMITED, const G4Transform3D &modelTransformation=G4Transform3D(), const G4ModelingParameters *=0, G4bool useFullExtent=false, const std::vector< G4PhysicalVolumeNodeID > &baseFullPVPath=std::vector< G4PhysicalVolumeNodeID >())
virtual	~G4PhysicalVolumeModel ()
void	DescribeYourselfTo (G4VGraphicsScene &)
G4String	GetCurrentDescription () const
G4String	GetCurrentTag () const
G4VPhysicalVolume *	GetTopPhysicalVolume () const
G4int	GetRequestedDepth () const
const G4VSolid *	GetClippingSolid () const
G4int	GetCurrentDepth () const
const G4Transform3D &	GetTransformation () const
G4VPhysicalVolume *	GetCurrentPV () const
G4int	GetCurrentPVCopyNo () const
G4LogicalVolume *	GetCurrentLV () const
G4Material *	GetCurrentMaterial () const
const G4Transform3D &	GetCurrentTransform () const
const std::vector< G4PhysicalVolumeNodeID > &	GetBaseFullPVPath () const
const std::vector< G4PhysicalVolumeNodeID > &	GetFullPVPath () const
const std::vector< G4PhysicalVolumeNodeID > &	GetDrawnPVPath () const
const std::map< G4String, G4AttDef > *	GetAttDefs () const
std::vector< G4AttValue > *	CreateCurrentAttValues () const
const std::map< G4int, G4int > &	GetMapOfDrawnTouchables () const
G4int	GetTotalDrawnTouchables ()
G4int	GetMaxFullDepth () const
const std::map< G4int, G4int > &	GetMapOfAllTouchables () const
G4int	GetTotalAllTouchables ()
void	SetRequestedDepth (G4int requestedDepth)
void	SetClippingSolid (G4VSolid *pClippingSolid)
void	SetClippingMode (ClippingMode mode)
G4bool	Validate (G4bool warn)
void	Abort () const
void	CurtailDescent () const
void	CalculateExtent ()
Public Member Functions inherited from G4VModel
	G4VModel (const G4ModelingParameters *=0)
virtual	~G4VModel ()
const G4ModelingParameters *	GetModelingParameters () const
const G4String &	GetType () const
const G4VisExtent &	GetExtent () const
const G4String &	GetGlobalDescription () const
const G4String &	GetGlobalTag () const
void	SetModelingParameters (const G4ModelingParameters *)
void	SetExtent (const G4VisExtent &)
void	SetType (const G4String &)
void	SetGlobalDescription (const G4String &)
void	SetGlobalTag (const G4String &)

Static Public Member Functions
static G4ModelingParameters::PVNameCopyNoPath	GetPVNameCopyNoPath (const std::vector< G4PhysicalVolumeNodeID > &)
static G4String	GetPVNamePathString (const std::vector< G4PhysicalVolumeNodeID > &)
Static Public Member Functions inherited from G4VModel
static const G4ModelingParameters *	GetCurrentModelingParameters ()
static void	SetCurrentModelingParameters (const G4ModelingParameters *)

Protected Member Functions
void	VisitGeometryAndGetVisReps (G4VPhysicalVolume *, G4int requestedDepth, const G4Transform3D &, G4VGraphicsScene &)
void	DescribeAndDescend (G4VPhysicalVolume *, G4int requestedDepth, G4LogicalVolume *, G4VSolid *, G4Material *, const G4Transform3D &, G4VGraphicsScene &)
virtual void	DescribeSolid (const G4Transform3D &theAT, G4VSolid *pSol, const G4VisAttributes *pVisAttribs, G4VGraphicsScene &sceneHandler)

Protected Attributes
G4VPhysicalVolume *	fpTopPV
G4String	fTopPVName
G4int	fTopPVCopyNo
G4int	fRequestedDepth
G4bool	fUseFullExtent
G4Transform3D	fTransform
G4int	fCurrentDepth
G4VPhysicalVolume *	fpCurrentPV
G4int	fCurrentPVCopyNo
G4LogicalVolume *	fpCurrentLV
G4Material *	fpCurrentMaterial
G4Transform3D	fCurrentTransform
std::vector< G4PhysicalVolumeNodeID >	fBaseFullPVPath
std::vector< G4PhysicalVolumeNodeID >	fFullPVPath
std::vector< G4PhysicalVolumeNodeID >	fDrawnPVPath
G4bool	fAbort
G4bool	fCurtailDescent
G4VSolid *	fpClippingSolid
ClippingMode	fClippingMode
G4int	fNClippers
std::map< G4int, G4int >	fMapDrawnTouchables
G4int	fTotalDrawnTouchables
std::map< G4int, G4int >	fMapAllTouchables
G4int	fTotalAllTouchables
G4int	fMaxFullDepth
Protected Attributes inherited from G4VModel
G4String	fType
G4String	fGlobalTag
G4String	fGlobalDescription
G4VisExtent	fExtent
const G4ModelingParameters *	fpMP

Detailed Description

Definition at line 82 of file G4PhysicalVolumeModel.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
UNLIMITED

Definition at line 86 of file G4PhysicalVolumeModel.hh.

{UNLIMITED = -1};

ClippingMode

enum G4PhysicalVolumeModel::ClippingMode

Enumerator
subtraction
intersection

Definition at line 88 of file G4PhysicalVolumeModel.hh.

{subtraction, intersection};

Constructor & Destructor Documentation

G4PhysicalVolumeModel()

G4PhysicalVolumeModel::G4PhysicalVolumeModel	(	G4VPhysicalVolume *	pVPV = 0,
		G4int	requestedDepth = UNLIMITED,
		const G4Transform3D &	modelTransformation = G4Transform3D(),
		const G4ModelingParameters *	pMP = 0,
		G4bool	useFullExtent = false,
		const std::vector< G4PhysicalVolumeNodeID > &	baseFullPVPath = std::vector<G4PhysicalVolumeNodeID>() )

Definition at line 62 of file G4PhysicalVolumeModel.cc.

: G4VModel           (pMP)
, fpTopPV            (pVPV)
, fTopPVCopyNo       (pVPV? pVPV->GetCopyNo(): 0)
, fRequestedDepth    (requestedDepth)
, fUseFullExtent     (useFullExtent)
, fTransform         (modelTransform)
, fCurrentDepth      (0)
, fpCurrentPV        (fpTopPV)
, fCurrentPVCopyNo   (fpTopPV? fpTopPV->GetCopyNo(): 0)
, fpCurrentLV        (fpTopPV? fpTopPV->GetLogicalVolume(): 0)
, fpCurrentMaterial  (fpCurrentLV? fpCurrentLV->GetMaterial(): 0)
, fCurrentTransform  (modelTransform)
, fBaseFullPVPath    (baseFullPVPath)
, fFullPVPath        (fBaseFullPVPath)
, fAbort             (false)
, fCurtailDescent    (false)
, fpClippingSolid    (0)
, fClippingMode      (subtraction)
, fNClippers         (0)
, fTotalDrawnTouchables (0)
, fTotalAllTouchables(0)
, fMaxFullDepth      (0)
{
  fType = "G4PhysicalVolumeModel";
 
  if (!fpTopPV) {
 
    // In some circumstances creating an "empty" G4PhysicalVolumeModel is
    // allowed, so I have supressed the G4Exception below.  If it proves to
    // be a problem we might have to re-instate it, but it is unlikley to
    // be used except by visualisation experts.  See, for example, /vis/list,
    // where it is used simply to get a list of G4AttDefs.
    //    G4Exception
    //    ("G4PhysicalVolumeModel::G4PhysicalVolumeModel",
    //     "modeling0010", FatalException, "Null G4PhysicalVolumeModel pointer.");
 
    fTopPVName = "NULL";
    fGlobalTag = "Empty";
    fGlobalDescription = "G4PhysicalVolumeModel " + fGlobalTag;
 
  } else {
 
    fTopPVName = fpTopPV -> GetName ();
    std::ostringstream oss;
    oss << fpTopPV->GetName() << ':' << fpTopPV->GetCopyNo()
    << " BasePath:" << fBaseFullPVPath;
    fGlobalTag = oss.str();
    fGlobalDescription = "G4PhysicalVolumeModel " + fGlobalTag;
    CalculateExtent ();
  }
}

Referenced by G4LogicalVolumeModel::DescribeYourselfTo(), and G4LogicalVolumeModel::G4LogicalVolumeModel().

~G4PhysicalVolumeModel()

G4PhysicalVolumeModel::~G4PhysicalVolumeModel ( )

virtual

Definition at line 121 of file G4PhysicalVolumeModel.cc.

{
  delete fpClippingSolid;
}

Member Function Documentation

Abort()

void G4PhysicalVolumeModel::Abort ( ) const

inline

Definition at line 275 of file G4PhysicalVolumeModel.hh.

{fAbort = true;}

CalculateExtent()

void G4PhysicalVolumeModel::CalculateExtent

(

)

Definition at line 148 of file G4PhysicalVolumeModel.cc.

{
  // To handle paramaterisations, set copy number and compute dimensions
  // to get extent right
  G4VPVParameterisation* pP = fpTopPV -> GetParameterisation ();
  if (pP) {
    fpTopPV -> SetCopyNo (fTopPVCopyNo);
    G4VSolid* solid = pP -> ComputeSolid (fTopPVCopyNo, fpTopPV);
    solid -> ComputeDimensions (pP, fTopPVCopyNo, fpTopPV);
  }
  if (fUseFullExtent) {
    fExtent = fpTopPV -> GetLogicalVolume () -> GetSolid () -> GetExtent ();
  } else {
    // Calculate extent of *drawn* volumes, i.e., ignoring culled, e.g.,
    // invisible volumes, by traversing the whole geometry hierarchy below
    // this physical volume.
    G4BoundingExtentScene beScene(this);
    const G4int tempRequestedDepth = fRequestedDepth;
    const G4Transform3D tempTransform = fTransform;
    const G4ModelingParameters* tempMP = fpMP;
    fRequestedDepth = -1;  // Always search to all depths to define extent.
    fTransform = G4Transform3D();  // Extent is in local cooridinates
    G4ModelingParameters mParams
      (0,      // No default vis attributes needed.
       G4ModelingParameters::wf,  // wireframe (not relevant for this).
       true,   // Global culling.
       true,   // Cull invisible volumes.
       false,  // Density culling.
       0.,     // Density (not relevant if density culling false).
       true,   // Cull daughters of opaque mothers.
       24);    // No of sides (not relevant for this operation).
    mParams.SetSpecialMeshRendering(true);  // Avoids traversing parameterisations
    fpMP = &mParams;
    DescribeYourselfTo (beScene);
    fExtent = beScene.GetBoundingExtent();
    fpMP = tempMP;
    fTransform = tempTransform;
    fRequestedDepth = tempRequestedDepth;
  }
  G4double radius = fExtent.GetExtentRadius();
  if (radius < 0.) {  // Nothing in the scene - revert to top extent
    fExtent = fpTopPV -> GetLogicalVolume () -> GetSolid () -> GetExtent ();
  }
  fExtent.Transform(fTransform);
}

Referenced by G4PhysicalVolumeModel(), and G4VVisCommandGeometrySet::Set().

CreateCurrentAttValues()

std::vector< G4AttValue > * G4PhysicalVolumeModel::CreateCurrentAttValues

(

)

const

Definition at line 1043 of file G4PhysicalVolumeModel.cc.

{
  std::vector<G4AttValue>* values = new std::vector<G4AttValue>;
 
  if (!fpCurrentLV) {
     G4Exception
        ("G4PhysicalVolumeModel::CreateCurrentAttValues",
         "modeling0004",
         JustWarning,
         "Current logical volume not defined.");
     return values;
  }
 
  std::ostringstream oss; oss << fFullPVPath;
  values->push_back(G4AttValue("PVPath", oss.str(),""));
 
  oss.str(""); oss << fBaseFullPVPath;
  values->push_back(G4AttValue("BasePVPath", oss.str(),""));
 
  values->push_back(G4AttValue("LVol", fpCurrentLV->GetName(),""));
  G4VSolid* pSol = fpCurrentLV->GetSolid();
 
  values->push_back(G4AttValue("Solid", pSol->GetName(),""));
 
  values->push_back(G4AttValue("EType", pSol->GetEntityType(),""));
 
  oss.str(""); oss << '\n' << *pSol;
  values->push_back(G4AttValue("DmpSol", oss.str(),""));
 
  const G4RotationMatrix localRotation = fpCurrentPV->GetObjectRotationValue();
  const G4ThreeVector& localTranslation = fpCurrentPV->GetTranslation();
  oss.str(""); oss << '\n' << G4Transform3D(localRotation,localTranslation);
  values->push_back(G4AttValue("LocalTrans", oss.str(),""));
 
  oss.str(""); oss << '\n' << pSol->GetExtent() << std::endl;
  values->push_back(G4AttValue("LocalExtent", oss.str(),""));
 
  oss.str(""); oss << '\n' << fCurrentTransform;
  values->push_back(G4AttValue("GlobalTrans", oss.str(),""));
 
  oss.str(""); oss << '\n' << (pSol->GetExtent()).Transform(fCurrentTransform) << std::endl;
  values->push_back(G4AttValue("GlobalExtent", oss.str(),""));
 
  G4String matName = fpCurrentMaterial? fpCurrentMaterial->GetName(): G4String("No material");
  values->push_back(G4AttValue("Material", matName,""));
 
  G4double matDensity = fpCurrentMaterial? fpCurrentMaterial->GetDensity(): 0.;
  values->push_back(G4AttValue("Density", G4BestUnit(matDensity,"Volumic Mass"),""));
 
  G4State matState = fpCurrentMaterial? fpCurrentMaterial->GetState(): kStateUndefined;
  oss.str(""); oss << matState;
  values->push_back(G4AttValue("State", oss.str(),""));
 
  G4double matRadlen = fpCurrentMaterial? fpCurrentMaterial->GetRadlen(): 0.;
  values->push_back(G4AttValue("Radlen", G4BestUnit(matRadlen,"Length"),""));
 
  G4Region* region = fpCurrentLV->GetRegion();
  G4String regionName = region? region->GetName(): G4String("No region");
  values->push_back(G4AttValue("Region", regionName,""));
 
  oss.str(""); oss << fpCurrentLV->IsRootRegion();
  values->push_back(G4AttValue("RootRegion", oss.str(),""));
 
  return values;
}

Referenced by G4VSceneHandler::LoadAtts(), G4ASCIITreeSceneHandler::RequestPrimitives(), and G4VisCommandsTouchable::SetNewValue().

CurtailDescent()

void G4PhysicalVolumeModel::CurtailDescent ( ) const

inline

Definition at line 278 of file G4PhysicalVolumeModel.hh.

{fCurtailDescent = true;}

Referenced by G4ASCIITreeSceneHandler::RequestPrimitives().

DescribeAndDescend()

void G4PhysicalVolumeModel::DescribeAndDescend ( G4VPhysicalVolume * pVPV, G4int requestedDepth, G4LogicalVolume * pLV, G4VSolid * pSol, G4Material * pMaterial, const G4Transform3D & theAT, G4VGraphicsScene & sceneHandler )

protected

Definition at line 423 of file G4PhysicalVolumeModel.cc.

{
  // Maintain useful data members...
  fpCurrentPV = pVPV;
  fCurrentPVCopyNo = pVPV->GetCopyNo();
  fpCurrentLV = pLV;
  fpCurrentMaterial = pMaterial;
 
  // Create a nodeID for use below - note the "drawn" flag is true
  G4int copyNo = fpCurrentPV->GetCopyNo();
  auto nodeID = G4PhysicalVolumeNodeID
  (fpCurrentPV,copyNo,fCurrentDepth,fCurrentTransform);
 
  // Update full path of physical volumes...
  fFullPVPath.push_back(nodeID);
 
  const G4RotationMatrix objectRotation = pVPV -> GetObjectRotationValue ();
  const G4ThreeVector&  translation     = pVPV -> GetTranslation ();
  G4Transform3D theLT (G4Transform3D (objectRotation, translation));
 
  // Compute the accumulated transformation...
  // Note that top volume's transformation relative to the world
  // coordinate system is specified in theAT == startingTransformation
  // = fTransform (see DescribeYourselfTo), so first time through the
  // volume's own transformation, which is only relative to its
  // mother, i.e., not relative to the world coordinate system, should
  // not be accumulated.
  G4Transform3D theNewAT (theAT);
  if (fCurrentDepth != 0) theNewAT = theAT * theLT;
  fCurrentTransform = theNewAT;
 
  const G4VisAttributes* pVisAttribs = pLV->GetVisAttributes();
  //  If the volume does not have any vis attributes, create it.
  G4VisAttributes* tempVisAtts = nullptr;
  if (!pVisAttribs) {
    if (fpMP->GetDefaultVisAttributes()) {
      tempVisAtts = new G4VisAttributes(*fpMP->GetDefaultVisAttributes());
    } else {
      tempVisAtts = new G4VisAttributes;
    }
    // The user may request /vis/viewer/set/colourByDensity.
    if (fpMP->GetCBDAlgorithmNumber() == 1) {
      // Algorithm 1: 3 parameters: Simple rainbow mapping.
      if (fpMP->GetCBDParameters().size() != 3) {
        G4Exception("G4PhysicalVolumeModelTouchable::DescribeAndDescend",
                    "modeling0014",
                    FatalErrorInArgument,
                    "Algorithm-parameter mismatch for Colour By Density");
      } else {
        const G4double d = pMaterial? pMaterial->GetDensity(): 0.;
        const G4double d0 = fpMP->GetCBDParameters()[0]; // Invisible d < d0.
        const G4double d1 = fpMP->GetCBDParameters()[1]; // Rainbow d0->d1->d2.
        const G4double d2 = fpMP->GetCBDParameters()[2]; // Blue d > d2.
        if (d < d0) { // Density < d0 is invisible.
          tempVisAtts->SetVisibility(false);
        } else { // Intermediate densities are on a spectrum.
          G4double red, green, blue;
          if (d < d1) {
            red = (d1-d)/(d1-d0); green = (d-d0)/(d1-d0); blue = 0.;
          } else if (d < d2) {
            red = 0.; green = (d2-d)/(d2-d1); blue = (d-d1)/(d2-d1);
          } else {  // Density >= d2 is blue.
            red = 0.; green = 0.; blue = 1.;
          }
          tempVisAtts->SetColour(G4Colour(red,green,blue));
        }
      }
    } else if (fpMP->GetCBDAlgorithmNumber() == 2) {
      // Algorithm 2
      // ...etc.
    }
    pVisAttribs = tempVisAtts;
  }
  // From here, can assume pVisAttribs is a valid pointer.  This is necessary
  // because PreAddSolid needs a vis attributes object.
 
  // Check if vis attributes are to be modified by a /vis/touchable/set/ command.
  const auto& vams = fpMP->GetVisAttributesModifiers();
  if (vams.size()) {
    // OK, we have some VAMs (Vis Attributes Modifiers).
    for (const auto& vam: vams) {
      const auto& vamPath = vam.GetPVNameCopyNoPath();
      if (vamPath.size() == fFullPVPath.size()) {
        // OK, we have a size match.
        // Check the volume name/copy number path.
        auto iVAMNameCopyNo = vamPath.begin();
        auto iPVNodeId = fFullPVPath.begin();
        for (; iVAMNameCopyNo != vamPath.end(); ++iVAMNameCopyNo, ++iPVNodeId) {
          if (!(
                iVAMNameCopyNo->GetName() ==
                iPVNodeId->GetPhysicalVolume()->GetName() &&
                iVAMNameCopyNo->GetCopyNo() ==
                iPVNodeId->GetPhysicalVolume()->GetCopyNo()
                )) {
            // This path element does NOT match.
            break;
          }
        }
        if (iVAMNameCopyNo == vamPath.end()) {
          // OK, the paths match (the above loop terminated normally).
          // Create a vis atts object for the modified vis atts.
          // It is static so that we may return a reliable pointer to it.
          static G4VisAttributes modifiedVisAtts;
          // Initialise it with the current operational vis atts and reset the pointer.
          modifiedVisAtts = *pVisAttribs;
          pVisAttribs = &modifiedVisAtts;
          const G4VisAttributes& transVisAtts = vam.GetVisAttributes();
          switch (vam.GetVisAttributesSignifier()) {
            case G4ModelingParameters::VASVisibility:
              modifiedVisAtts.SetVisibility(transVisAtts.IsVisible());
              break;
            case G4ModelingParameters::VASDaughtersInvisible:
              modifiedVisAtts.SetDaughtersInvisible
              (transVisAtts.IsDaughtersInvisible());
              break;
            case G4ModelingParameters::VASColour:
              modifiedVisAtts.SetColour(transVisAtts.GetColour());
              break;
            case G4ModelingParameters::VASLineStyle:
              modifiedVisAtts.SetLineStyle(transVisAtts.GetLineStyle());
              break;
            case G4ModelingParameters::VASLineWidth:
              modifiedVisAtts.SetLineWidth(transVisAtts.GetLineWidth());
              break;
            case G4ModelingParameters::VASForceWireframe:
              if (transVisAtts.IsForceDrawingStyle()) {
                if (transVisAtts.GetForcedDrawingStyle() ==
                    G4VisAttributes::wireframe) {
                  modifiedVisAtts.SetForceWireframe(true);
                }
              }
              break;
            case G4ModelingParameters::VASForceSolid:
              if (transVisAtts.IsForceDrawingStyle()) {
                if (transVisAtts.GetForcedDrawingStyle() ==
                    G4VisAttributes::solid) {
                  modifiedVisAtts.SetForceSolid(true);
                }
              }
              break;
            case G4ModelingParameters::VASForceCloud:
              if (transVisAtts.IsForceDrawingStyle()) {
                if (transVisAtts.GetForcedDrawingStyle() ==
                    G4VisAttributes::cloud) {
                  modifiedVisAtts.SetForceCloud(true);
                }
              }
              break;
            case G4ModelingParameters::VASForceNumberOfCloudPoints:
              modifiedVisAtts.SetForceNumberOfCloudPoints
              (transVisAtts.GetForcedNumberOfCloudPoints());
              break;
            case G4ModelingParameters::VASForceAuxEdgeVisible:
              if (transVisAtts.IsForceAuxEdgeVisible()) {
                modifiedVisAtts.SetForceAuxEdgeVisible
                (transVisAtts.IsForcedAuxEdgeVisible());
              }
              break;
            case G4ModelingParameters::VASForceLineSegmentsPerCircle:
              modifiedVisAtts.SetForceLineSegmentsPerCircle
              (transVisAtts.GetForcedLineSegmentsPerCircle());
              break;
          }
        }
      }
    }
  }
 
  // Check for special mesh rendering
  if (fpMP->IsSpecialMeshRendering()) {
    G4bool potentialG4Mesh = false;
    if (fpMP->GetSpecialMeshVolumes().empty()) {
      // No volumes specified - all are potentially possible
      potentialG4Mesh = true;
    } else {
      // Name and (optionally) copy number of container volume is specified
      for (const auto& pvNameCopyNo: fpMP->GetSpecialMeshVolumes()) {
        if (pVPV->GetName() == pvNameCopyNo.GetName()) {
          // We have a name match
          if (pvNameCopyNo.GetCopyNo() < 0) {
            // Any copy number is OK
            potentialG4Mesh = true;
          } else {
            if (pVPV->GetCopyNo() == pvNameCopyNo.GetCopyNo()) {
              // We have a name and copy number match
              potentialG4Mesh = true;
            }
          }
        }
      }
    }
    if (potentialG4Mesh) {
      // Create - or at least attempt to create - a mesh. If it cannot be created
      // out of this pVPV the type will be "invalid".
      G4Mesh mesh(pVPV,theNewAT);
      if (mesh.GetMeshType() != G4Mesh::invalid) {
        // Create "artificial" nodeID to represent the replaced volumes
        G4int artCopyNo = 0;
        auto artPV = mesh.GetParameterisedVolume();
        auto artDepth = fCurrentDepth + 1;
        auto artNodeID = G4PhysicalVolumeNodeID(artPV,artCopyNo,artDepth);
        fFullPVPath.push_back(artNodeID);
        fDrawnPVPath.push_back(artNodeID);
        sceneHandler.AddCompound(mesh);
        fFullPVPath.pop_back();
        fDrawnPVPath.pop_back();
        delete tempVisAtts;  // Needs cleaning up (Coverity warning!!)
        return;  // Mesh found and processed - nothing more to do.
      }  // else continue processing
    }
  }
 
  auto currentFullDepth = fCurrentDepth + (G4int)fBaseFullPVPath.size();
 
  // Respect transparency by depth set by user
  const auto& transparencyByDepthParameter = fpMP->GetTransparencyByDepth();
  const auto& maxDepth = sceneHandler.GetMaxGeometryDepth();
  if (transparencyByDepthParameter > 0 && maxDepth > 0) {
    const auto& option = fpMP->GetTransparencyByDepthOption();
    G4double alpha = 1.;  // Multiplies pre-existing opacity
    switch (option) {
      case 1: {  // Unwrap - simply make invisible by depth
        if (currentFullDepth <= transparencyByDepthParameter) alpha = 0.;
        break;
      }
      case 2: {  // Fade a layer at a time
        alpha = currentFullDepth - transparencyByDepthParameter;
        alpha = std::min(1.,std::max(0.,alpha));
        break;
      }
      case 3: {  // X-ray: progressive transparancy throughout depth
        // Algorithm by Andrea Barresi January 2025.
        // Linear function from 0 to 1 in the [0,maxdepth] range
        // Shift according to k value so that:
        // alpha = 1 for all volumes when k = 0,
        // alpha = 0 for all volumes when k = maxDepth.
        const auto& k = transparencyByDepthParameter;
        alpha = G4double(currentFullDepth) / maxDepth - 2 * k / maxDepth + 1;
        alpha = std::min(1., std::max(0., alpha));
      }
      default: {}
    }
    if (alpha < 1.) {
      // As above, create a vis atts object for the modified vis atts.
      // It is static so that we may return a reliable pointer to it.
      static G4VisAttributes transparencyByDepthVisAtts;
      // Initialise it with the current operational vis atts and reset the pointer.
      transparencyByDepthVisAtts = *pVisAttribs;
      pVisAttribs = &transparencyByDepthVisAtts;
      // Adjust the transparency (multiply any existing opacity)
      auto colour = pVisAttribs->GetColour();
      colour.SetAlpha(colour.GetAlpha()*alpha);
      transparencyByDepthVisAtts.SetColour(colour);
    }
  }
 
  // Make decision to draw...
  G4bool thisToBeDrawn = true;
 
  // There are various reasons why this volume
  // might not be drawn...
  G4bool culling = fpMP->IsCulling();
  G4bool cullingInvisible = fpMP->IsCullingInvisible();
  G4bool markedVisible
  = pVisAttribs->IsVisible() && pVisAttribs->GetColour().GetAlpha() > 0;
  G4bool cullingLowDensity = fpMP->IsDensityCulling();
  G4double density = pMaterial? pMaterial->GetDensity(): 0;
  G4double densityCut = fpMP -> GetVisibleDensity ();
 
  // 1) Global culling is on....
  if (culling) {
    // 2) Culling of invisible volumes is on...
    if (cullingInvisible) {
      // 3) ...and the volume is marked not visible...
      if (!markedVisible) thisToBeDrawn = false;
    }
    // 4) Or culling of low density volumes is on...
    if (cullingLowDensity) {
      // 5) ...and density is less than cut value...
      if (density < densityCut) thisToBeDrawn = false;
    }
  }
  // 6) The user has asked for all further traversing to be aborted...
  if (fAbort) thisToBeDrawn = false;
 
  // Set "drawn" flag (it was true by default) - thisToBeDrawn may be false
  nodeID.SetDrawn(thisToBeDrawn);
 
  fMapAllTouchables[currentFullDepth]++;  // Increment for every touchable at each depth
  if (fMaxFullDepth < currentFullDepth) fMaxFullDepth = currentFullDepth;
 
  if (thisToBeDrawn) {
 
    // Update path of drawn physical volumes...
    fDrawnPVPath.push_back(nodeID);
 
    fMapDrawnTouchables[currentFullDepth]++;  // Increment for every touchable drawn at each depth
 
    if (fpMP->IsExplode() && fDrawnPVPath.size() == 1) {
      // For top-level drawn volumes, explode along radius...
      G4Transform3D centering = G4Translate3D(fpMP->GetExplodeCentre());
      G4Transform3D centred = centering.inverse() * theNewAT;
      G4Scale3D oldScale;
      G4Rotate3D oldRotation;
      G4Translate3D oldTranslation;
      centred.getDecomposition(oldScale, oldRotation, oldTranslation);
      G4double explodeFactor = fpMP->GetExplodeFactor();
      G4Translate3D newTranslation =
      G4Translate3D(explodeFactor * oldTranslation.dx(),
                    explodeFactor * oldTranslation.dy(),
                    explodeFactor * oldTranslation.dz());
      theNewAT = centering * newTranslation * oldRotation * oldScale;
    }
 
    DescribeSolid (theNewAT, pSol, pVisAttribs, sceneHandler);
 
  }
 
  // Make decision to draw daughters, if any.  There are various
  // reasons why daughters might not be drawn...
 
  // First, reasons that do not depend on culling policy...
  G4int nDaughters = (G4int)pLV->GetNoDaughters();
  G4bool daughtersToBeDrawn = true;
  // 1) There are no daughters...
  if (!nDaughters) daughtersToBeDrawn = false;
  // 2) We are at the limit if requested depth...
  else if (requestedDepth == 0) daughtersToBeDrawn = false;
  // 3) The user has asked for all further traversing to be aborted...
  else if (fAbort) daughtersToBeDrawn = false;
  // 4) The user has asked that the descent be curtailed...
  else if (fCurtailDescent) daughtersToBeDrawn = false;
 
  // Now, reasons that depend on culling policy...
  else {
    G4bool daughtersInvisible = pVisAttribs->IsDaughtersInvisible();
    // Culling of covered daughters request.  This is computed in
    // G4VSceneHandler::CreateModelingParameters() depending on view
    // parameters...
    G4bool cullingCovered = fpMP->IsCullingCovered();
    G4bool surfaceDrawing =
    fpMP->GetDrawingStyle() == G4ModelingParameters::hsr ||
    fpMP->GetDrawingStyle() == G4ModelingParameters::hlhsr;
    if (pVisAttribs->IsForceDrawingStyle()) {
      switch (pVisAttribs->GetForcedDrawingStyle()) {
        default:
        case G4VisAttributes::wireframe: surfaceDrawing = false; break;
        case G4VisAttributes::solid: surfaceDrawing = true; break;
      }
    }
    G4bool opaque = pVisAttribs->GetColour().GetAlpha() >= 1.;
    // 5) Global culling is on....
    if (culling) {
      // 6) ..and culling of invisible volumes is on...
      if (cullingInvisible) {
        // 7) ...and the mother requests daughters invisible
        if (daughtersInvisible) daughtersToBeDrawn = false;
      }
      // 8) Or culling of covered daughters is requested...
      if (cullingCovered) {
        // 9) ...and surface drawing is operating...
        if (surfaceDrawing) {
          // 10) ...but only if mother is visible...
          if (thisToBeDrawn) {
            // 11) ...and opaque...
            if (opaque) daughtersToBeDrawn = false;
          }
        }
      }
    }
  }
 
  if (daughtersToBeDrawn) {
    for (G4int iDaughter = 0; iDaughter < nDaughters; iDaughter++) {
      // Store daughter pVPV in local variable ready for recursion...
      G4VPhysicalVolume* pDaughterVPV = pLV -> GetDaughter (iDaughter);
      // Descend the geometry structure recursively...
      fCurrentDepth++;
      VisitGeometryAndGetVisReps
      (pDaughterVPV, requestedDepth - 1, theNewAT, sceneHandler);
      fCurrentDepth--;
    }
  }
 
  // Clean up
  delete tempVisAtts;
 
  // Reset for normal descending of next volume at this level...
  fCurtailDescent = false;
 
  // Pop item from paths physical volumes...
  fFullPVPath.pop_back();
  if (thisToBeDrawn) {
    fDrawnPVPath.pop_back();
  }
}

Referenced by VisitGeometryAndGetVisReps().

DescribeSolid()

void G4PhysicalVolumeModel::DescribeSolid ( const G4Transform3D & theAT, G4VSolid * pSol, const G4VisAttributes * pVisAttribs, G4VGraphicsScene & sceneHandler )

protectedvirtual

Reimplemented in G4LogicalVolumeModel.

Definition at line 870 of file G4PhysicalVolumeModel.cc.

{
  G4DisplacedSolid* pSectionSolid = fpMP->GetSectionSolid();
  G4DisplacedSolid* pCutawaySolid = fpMP->GetCutawaySolid();
 
  if (fNClippers <= 0 || fNClippers > 1) {
 
    // Normal case - no clipping, etc. - or, illegally, more than one of those
    sceneHandler.PreAddSolid (theAT, *pVisAttribs);
    pSol -> DescribeYourselfTo (sceneHandler);  // Standard treatment.
    sceneHandler.PostAddSolid ();
 
  } else {  // fNClippers == 1
 
    G4VSolid*         pResultantSolid = nullptr;
    G4DisplacedSolid* pDisplacedSolid = nullptr;
 
    if (fpClippingSolid) {
      pDisplacedSolid = new G4DisplacedSolid("clipper", fpClippingSolid, theAT.inverse());
      switch (fClippingMode) {
        case subtraction:
          if (SubtractionBoundingLimits(pSol)) {
            pResultantSolid = new G4SubtractionSolid
            ("subtracted_clipped_solid", pSol, pDisplacedSolid);
          }
          break;
        case intersection:
          if (IntersectionBoundingLimits(pSol, pDisplacedSolid)) {
            pResultantSolid = new G4IntersectionSolid
            ("intersected_clipped_solid", pSol, pDisplacedSolid);
          }
          break;
      }
 
    } else if (pSectionSolid) {
      pDisplacedSolid = new G4DisplacedSolid("intersector", pSectionSolid, theAT.inverse());
      if (IntersectionBoundingLimits(pSol, pDisplacedSolid)) {
        pResultantSolid = new G4IntersectionSolid("sectioned_solid", pSol, pDisplacedSolid);
      }
 
    } else if (pCutawaySolid) {
      pDisplacedSolid = new G4DisplacedSolid("cutaway", pCutawaySolid, theAT.inverse());
      switch (fpMP->GetCutawayMode()) {
        case G4ModelingParameters::cutawayUnion:
          if (SubtractionBoundingLimits(pSol)) {
            pResultantSolid = new G4SubtractionSolid("cutaway_solid", pSol, pDisplacedSolid);
          }
          break;
        case G4ModelingParameters::cutawayIntersection:
          if (IntersectionBoundingLimits(pSol, pDisplacedSolid)) {
            pResultantSolid = new G4IntersectionSolid("cutaway_solid", pSol, pDisplacedSolid);
          }
          break;
      }
    }
 
    if (pResultantSolid) {
      sceneHandler.PreAddSolid (theAT, *pVisAttribs);
      pResultantSolid -> DescribeYourselfTo (sceneHandler);
      sceneHandler.PostAddSolid ();
    }
 
    delete pResultantSolid;
    delete pDisplacedSolid;
  }
}

Referenced by DescribeAndDescend().

DescribeYourselfTo()

void G4PhysicalVolumeModel::DescribeYourselfTo ( G4VGraphicsScene & sceneHandler )

virtual

Implements G4VModel.

Definition at line 194 of file G4PhysicalVolumeModel.cc.

{
  if (!fpTopPV) {
    G4Exception
    ("G4PhysicalVolumeModel::DescribeYourselfTo",
     "modeling0012", FatalException, "No model.");
    return;  // Should never reach here, but keeps Coverity happy.
  }
 
  if (!fpMP) {
    G4Exception
    ("G4PhysicalVolumeModel::DescribeYourselfTo",
     "modeling0013", FatalException, "No modeling parameters.");
    return;  // Should never reach here, but keeps Coverity happy.
  }
 
  fNClippers = 0;
  G4DisplacedSolid* pSectionSolid = fpMP->GetSectionSolid();
  G4DisplacedSolid* pCutawaySolid = fpMP->GetCutawaySolid();
  if (fpClippingSolid) fNClippers++;
  if (pSectionSolid)   fNClippers++;
  if (pCutawaySolid)   fNClippers++;
  if (fNClippers > 1) {
    G4ExceptionDescription ed;
    ed << "More than one solid cutter/clipper:";
    if (fpClippingSolid) ed << "\nclipper in force";
    if (pSectionSolid)   ed << "\nsectioner in force";
    if (pCutawaySolid)   ed << "\ncutaway in force";
    G4Exception("G4PhysicalVolumeModel::DescribeSolid", "modeling0016", JustWarning, ed);
  }
 
  G4Transform3D startingTransformation = fTransform;
 
  fMapDrawnTouchables.clear();  // Keeps count of touchables drawn at each depth
  fMapAllTouchables.clear();  // Keeps count of all touchables at each depth
 
  VisitGeometryAndGetVisReps
    (fpTopPV,
     fRequestedDepth,
     startingTransformation,
     sceneHandler);
 
  fTotalDrawnTouchables = 0;
  for (const auto& entry : fMapDrawnTouchables) {
    fTotalDrawnTouchables += entry.second;
  }
 
  fTotalAllTouchables = 0;
  for (const auto& entry : fMapAllTouchables) {
    fTotalAllTouchables += entry.second;
  }
 
  // Reset or clear data...
  fCurrentDepth     = 0;
  fpCurrentPV       = fpTopPV;
  fCurrentPVCopyNo  = fpTopPV->GetCopyNo();
  fpCurrentLV       = fpTopPV->GetLogicalVolume();
  fpCurrentMaterial = fpCurrentLV? fpCurrentLV->GetMaterial(): 0;
  fFullPVPath       = fBaseFullPVPath;
  fDrawnPVPath.clear();
  fAbort            = false;
  fCurtailDescent   = false;
}

Referenced by CalculateExtent(), DescribeSolid(), G4LogicalVolumeModel::DescribeYourselfTo(), G4VSceneHandler::Draw3DRectMeshAsDots(), G4VSceneHandler::Draw3DRectMeshAsSurfaces(), G4VisManager::DrawGeometry(), G4VSceneHandler::DrawTetMeshAsDots(), G4VSceneHandler::DrawTetMeshAsSurfaces(), G4ASCIITreeSceneHandler::EndModeling(), G4TouchableUtils::FindTouchableProperties(), G4VisCommandSceneAddLocalAxes::SetNewValue(), G4VisCommandSceneAddVolume::SetNewValue(), G4VisCommandSetTouchable::SetNewValue(), G4VisCommandSetVolumeForField::SetNewValue(), G4VisCommandsTouchable::SetNewValue(), G4VisCommandViewerCentreOn::SetNewValue(), and G4VtkUnstructuredGridPipeline::SetUnstructuredGridData().

GetAttDefs()

const std::map< G4String, G4AttDef > * G4PhysicalVolumeModel::GetAttDefs

(

)

const

Definition at line 962 of file G4PhysicalVolumeModel.cc.

{
    G4bool isNew;
    std::map<G4String,G4AttDef>* store
      = G4AttDefStore::GetInstance("G4PhysicalVolumeModel", isNew);
    if (isNew) {
      (*store)["PVPath"] =
      G4AttDef("PVPath","Physical Volume Path","Physics","","G4String");
      (*store)["BasePVPath"] =
      G4AttDef("BasePVPath","Base Physical Volume Path","Physics","","G4String");
      (*store)["LVol"] =
      G4AttDef("LVol","Logical Volume","Physics","","G4String");
      (*store)["Solid"] =
      G4AttDef("Solid","Solid Name","Physics","","G4String");
      (*store)["EType"] =
      G4AttDef("EType","Entity Type","Physics","","G4String");
      (*store)["DmpSol"] =
      G4AttDef("DmpSol","Dump of Solid properties","Physics","","G4String");
      (*store)["LocalTrans"] =
      G4AttDef("LocalTrans","Local transformation of volume","Physics","","G4String");
      (*store)["LocalExtent"] =
      G4AttDef("LocalExtent","Local extent of volume","Physics","","G4String");
      (*store)["GlobalTrans"] =
      G4AttDef("GlobalTrans","Global transformation of volume","Physics","","G4String");
      (*store)["GlobalExtent"] =
      G4AttDef("GlobalExtent","Global extent of volume","Physics","","G4String");
      (*store)["Material"] =
      G4AttDef("Material","Material Name","Physics","","G4String");
      (*store)["Density"] =
      G4AttDef("Density","Material Density","Physics","G4BestUnit","G4double");
      (*store)["State"] =
      G4AttDef("State","Material State (enum undefined,solid,liquid,gas)","Physics","","G4String");
      (*store)["Radlen"] =
      G4AttDef("Radlen","Material Radiation Length","Physics","G4BestUnit","G4double");
      (*store)["Region"] =
      G4AttDef("Region","Cuts Region","Physics","","G4String");
      (*store)["RootRegion"] =
      G4AttDef("RootRegion","Root Region (0/1 = false/true)","Physics","","G4bool");
    }
  return store;
}

Referenced by G4VSceneHandler::LoadAtts(), G4ASCIITreeSceneHandler::RequestPrimitives(), G4VisCommandList::SetNewValue(), and G4VisCommandsTouchable::SetNewValue().

GetBaseFullPVPath()

const std::vector< G4PhysicalVolumeNodeID > & G4PhysicalVolumeModel::GetBaseFullPVPath ( ) const

inline

Definition at line 203 of file G4PhysicalVolumeModel.hh.

  {return fBaseFullPVPath;}

GetClippingSolid()

const G4VSolid * G4PhysicalVolumeModel::GetClippingSolid ( ) const

inline

Definition at line 179 of file G4PhysicalVolumeModel.hh.

  {return fpClippingSolid;}

GetCurrentDepth()

G4int G4PhysicalVolumeModel::GetCurrentDepth ( ) const

inline

Definition at line 182 of file G4PhysicalVolumeModel.hh.

{return fCurrentDepth;}

Referenced by G4GMocrenFileSceneHandler::AddSolid().

GetCurrentDescription()

G4String G4PhysicalVolumeModel::GetCurrentDescription ( ) const

virtual

Reimplemented from G4VModel.

Definition at line 271 of file G4PhysicalVolumeModel.cc.

{
  return "G4PhysicalVolumeModel " + GetCurrentTag ();
}

GetCurrentLV()

G4LogicalVolume * G4PhysicalVolumeModel::GetCurrentLV ( ) const

inline

Definition at line 194 of file G4PhysicalVolumeModel.hh.

{return fpCurrentLV;}

Referenced by G4GMocrenFileSceneHandler::AddPrimitive(), G4OpenInventorSceneHandler::AddPrimitive(), and G4ASCIITreeSceneHandler::RequestPrimitives().

GetCurrentMaterial()

G4Material * G4PhysicalVolumeModel::GetCurrentMaterial ( ) const

inline

Definition at line 197 of file G4PhysicalVolumeModel.hh.

{return fpCurrentMaterial;}

Referenced by G4GMocrenFileSceneHandler::AddSolid(), G4GMocrenFileSceneHandler::AddSolid(), and G4ASCIITreeSceneHandler::RequestPrimitives().

GetCurrentPV()

G4VPhysicalVolume * G4PhysicalVolumeModel::GetCurrentPV ( ) const

inline

Definition at line 188 of file G4PhysicalVolumeModel.hh.

{return fpCurrentPV;}

Referenced by G4GMocrenFileSceneHandler::AddPrimitive(), G4GMocrenFileSceneHandler::AddSolid(), and G4ASCIITreeSceneHandler::RequestPrimitives().

GetCurrentPVCopyNo()

G4int G4PhysicalVolumeModel::GetCurrentPVCopyNo ( ) const

inline

Definition at line 191 of file G4PhysicalVolumeModel.hh.

{return fCurrentPVCopyNo;}

GetCurrentTag()

G4String G4PhysicalVolumeModel::GetCurrentTag ( ) const

virtual

Reimplemented from G4VModel.

Definition at line 259 of file G4PhysicalVolumeModel.cc.

{
  if (fpCurrentPV) {
    std::ostringstream o;
    o << fpCurrentPV -> GetCopyNo ();
    return fpCurrentPV -> GetName () + ":" + o.str();
  }
  else {
    return "WARNING: NO CURRENT VOLUME - global tag is " + fGlobalTag;
  }
}

Referenced by GetCurrentDescription().

GetCurrentTransform()

const G4Transform3D & G4PhysicalVolumeModel::GetCurrentTransform ( ) const

inline

Definition at line 200 of file G4PhysicalVolumeModel.hh.

{return fCurrentTransform;}

Referenced by G4ASCIITreeSceneHandler::RequestPrimitives(), and G4VisCommandsTouchable::SetNewValue().

GetDrawnPVPath()

const std::vector< G4PhysicalVolumeNodeID > & G4PhysicalVolumeModel::GetDrawnPVPath ( ) const

inline

Definition at line 217 of file G4PhysicalVolumeModel.hh.

  {return fDrawnPVPath;}

Referenced by G4GMocrenFileSceneHandler::AddSolid(), G4VTreeSceneHandler::PreAddSolid(), and G4ASCIITreeSceneHandler::RequestPrimitives().

GetFullPVPath()

const std::vector< G4PhysicalVolumeNodeID > & G4PhysicalVolumeModel::GetFullPVPath ( ) const

inline

Definition at line 210 of file G4PhysicalVolumeModel.hh.

  {return fFullPVPath;}

Referenced by G4OpenGLQtViewer::addPVSceneTreeElement(), and G4RayTracerSceneHandler::BuildVisAttsMap().

GetMapOfAllTouchables()

const std::map< G4int, G4int > & G4PhysicalVolumeModel::GetMapOfAllTouchables ( ) const

inline

Definition at line 254 of file G4PhysicalVolumeModel.hh.

{return fMapAllTouchables;}

GetMapOfDrawnTouchables()

const std::map< G4int, G4int > & G4PhysicalVolumeModel::GetMapOfDrawnTouchables ( ) const

inline

Definition at line 245 of file G4PhysicalVolumeModel.hh.

{return fMapDrawnTouchables;}

GetMaxFullDepth()

G4int G4PhysicalVolumeModel::GetMaxFullDepth ( ) const

inline

Definition at line 251 of file G4PhysicalVolumeModel.hh.

{return fMaxFullDepth;}

GetPVNameCopyNoPath()

G4ModelingParameters::PVNameCopyNoPath G4PhysicalVolumeModel::GetPVNameCopyNoPath ( const std::vector< G4PhysicalVolumeNodeID > & path )

static

Definition at line 126 of file G4PhysicalVolumeModel.cc.

{
  G4ModelingParameters::PVNameCopyNoPath PVNameCopyNoPath;
  for (const auto& node: path) {
    PVNameCopyNoPath.push_back
    (G4ModelingParameters::PVNameCopyNo
     (node.GetPhysicalVolume()->GetName(),node.GetCopyNo()));
  }
  return PVNameCopyNoPath;
}

Referenced by G4VViewer::TouchableSetColour(), G4VViewer::TouchableSetVisibility(), and G4VVisCommand::Twinkle().

GetPVNamePathString()

G4String G4PhysicalVolumeModel::GetPVNamePathString ( const std::vector< G4PhysicalVolumeNodeID > & path )

static

Definition at line 138 of file G4PhysicalVolumeModel.cc.

{
  std::ostringstream oss;
  oss << path;
  return oss.str();
}

Referenced by G4VViewer::TouchableSetColour(), and G4VViewer::TouchableSetVisibility().

GetRequestedDepth()

G4int G4PhysicalVolumeModel::GetRequestedDepth ( ) const

inline

Definition at line 177 of file G4PhysicalVolumeModel.hh.

{return fRequestedDepth;}

Referenced by G4ASCIITreeSceneHandler::EndModeling().

GetTopPhysicalVolume()

G4VPhysicalVolume * G4PhysicalVolumeModel::GetTopPhysicalVolume ( ) const

inline

Definition at line 175 of file G4PhysicalVolumeModel.hh.

{return fpTopPV;}

Referenced by G4GMocrenFileSceneHandler::AddSolid(), and G4ASCIITreeSceneHandler::EndModeling().

GetTotalAllTouchables()

G4int G4PhysicalVolumeModel::GetTotalAllTouchables ( )

inline

Definition at line 257 of file G4PhysicalVolumeModel.hh.

{return fTotalAllTouchables;}

GetTotalDrawnTouchables()

G4int G4PhysicalVolumeModel::GetTotalDrawnTouchables ( )

inline

Definition at line 248 of file G4PhysicalVolumeModel.hh.

{return fTotalDrawnTouchables;}

GetTransformation()

const G4Transform3D & G4PhysicalVolumeModel::GetTransformation ( ) const

inline

Definition at line 185 of file G4PhysicalVolumeModel.hh.

{return fTransform;}

SetClippingMode()

void G4PhysicalVolumeModel::SetClippingMode ( ClippingMode mode )

inline

Definition at line 268 of file G4PhysicalVolumeModel.hh.

                                           {
    fClippingMode = mode;
  }

Referenced by G4VisCommandSceneAddVolume::SetNewValue().

SetClippingSolid()

void G4PhysicalVolumeModel::SetClippingSolid ( G4VSolid * pClippingSolid )

inline

Definition at line 264 of file G4PhysicalVolumeModel.hh.

                                                   {
    fpClippingSolid = pClippingSolid;
  }

Referenced by G4VisCommandSceneAddVolume::SetNewValue().

SetRequestedDepth()

void G4PhysicalVolumeModel::SetRequestedDepth ( G4int requestedDepth )

inline

Definition at line 260 of file G4PhysicalVolumeModel.hh.

                                                {
    fRequestedDepth = requestedDepth;
  }

Validate()

G4bool G4PhysicalVolumeModel::Validate ( G4bool warn )

virtual

Reimplemented from G4VModel.

Definition at line 941 of file G4PhysicalVolumeModel.cc.

{
// Not easy to see how to validate this sort of model. Previously there was
// a check that a volume of the same name (fTopPVName) existed somewhere in
// the geometry tree but under some circumstances this consumed lots of CPU
// time. Instead, let us simply check that the volume (fpTopPV) exists in the
// physical volume store.
  const auto& pvStore = G4PhysicalVolumeStore::GetInstance();
  auto iterator = find(pvStore->begin(),pvStore->end(),fpTopPV);
  if (iterator == pvStore->end()) {
    if (warn) {
      G4ExceptionDescription ed;
      ed << "Attempt to validate a volume that is no longer in the physical volume store.";
      G4Exception("G4PhysicalVolumeModel::Validate", "modeling0015", JustWarning, ed);
    }
    return false;
  } else {
    return true;
  }
}

Referenced by G4VisCommandSceneAddVolume::SetNewValue().

VisitGeometryAndGetVisReps()

void G4PhysicalVolumeModel::VisitGeometryAndGetVisReps ( G4VPhysicalVolume * pVPV, G4int requestedDepth, const G4Transform3D & theAT, G4VGraphicsScene & sceneHandler )

protected

Definition at line 276 of file G4PhysicalVolumeModel.cc.

{
  // Visits geometry structure to a given depth (requestedDepth), starting
  //   at given physical volume with given starting transformation and
  //   describes volumes to the scene handler.
  // requestedDepth < 0 (default) implies full visit.
  // theAT is the Accumulated Transformation.
 
  // Find corresponding logical volume and (later) solid, storing in
  // local variables to preserve re-entrancy.
  G4LogicalVolume* pLV  = pVPV -> GetLogicalVolume ();
  G4VSolid* pSol = nullptr;
  G4Material* pMaterial = nullptr;
 
  if (!(pVPV -> IsReplicated ())) {
    // Non-replicated physical volume.
    pSol = pLV -> GetSolid ();
    pMaterial = pLV -> GetMaterial ();
    DescribeAndDescend (pVPV, requestedDepth, pLV, pSol, pMaterial,
            theAT, sceneHandler);
  }
  else {
    // Replicated or parametrised physical volume.
    EAxis axis;
    G4int nReplicas;
    G4double width;
    G4double offset;
    G4bool consuming;
    pVPV -> GetReplicationData (axis, nReplicas, width,  offset, consuming);
    G4int nBegin = 0;
    G4int nEnd = nReplicas;
    if (fCurrentDepth == 0) { // i.e., top volume
      nBegin = fTopPVCopyNo;  // Describe only one volume, namely the one
      nEnd = nBegin + 1;      // specified by the given copy number.
    }
    G4VPVParameterisation* pP = pVPV -> GetParameterisation ();
    if (pP) {  // Parametrised volume.
      for (int n = nBegin; n < nEnd; n++) {
        pSol = pP -> ComputeSolid (n, pVPV);
        pP -> ComputeTransformation (n, pVPV);
        pSol -> ComputeDimensions (pP, n, pVPV);
        pVPV -> SetCopyNo (n);
        fCurrentPVCopyNo = n;
        // Create a touchable of current parent for ComputeMaterial.
        // fFullPVPath has not been updated yet so at this point it
        // corresponds to the parent.
        G4PhysicalVolumeModelTouchable parentTouchable(fFullPVPath);
        pMaterial = pP -> ComputeMaterial (n, pVPV, &parentTouchable);
        DescribeAndDescend (pVPV, requestedDepth, pLV, pSol, pMaterial,
                            theAT, sceneHandler);
      }
    }
    else {  // Plain replicated volume.  From geometry_guide.txt...
      // The replica's positions are claculated by means of a linear formula.
      // Replication may occur along:
      // 
      // o Cartesian axes (kXAxis,kYAxis,kZAxis)
      // 
      //   The replications, of specified width have coordinates of
      //   form (-width*(nReplicas-1)*0.5+n*width,0,0) where n=0.. nReplicas-1
      //   for the case of kXAxis, and are unrotated.
      // 
      // o Radial axis (cylindrical polar) (kRho)
      // 
      //   The replications are cons/tubs sections, centred on the origin
      //   and are unrotated.
      //   They have radii of width*n+offset to width*(n+1)+offset
      //                      where n=0..nReplicas-1
      // 
      // o Phi axis (cylindrical polar) (kPhi)
      //   The replications are `phi sections' or wedges, and of cons/tubs form
      //   They have phi of offset+n*width to offset+(n+1)*width where
      //   n=0..nReplicas-1
      // 
      pSol = pLV -> GetSolid ();
      pMaterial = pLV -> GetMaterial ();
      G4ThreeVector originalTranslation = pVPV -> GetTranslation ();
      G4RotationMatrix* pOriginalRotation = pVPV -> GetRotation ();
      G4double originalRMin = 0., originalRMax = 0.;
      if (axis == kRho && pSol->GetEntityType() == "G4Tubs") {
        originalRMin = ((G4Tubs*)pSol)->GetInnerRadius();
        originalRMax = ((G4Tubs*)pSol)->GetOuterRadius();
      }
      G4bool visualisable = true;
      for (int n = nBegin; n < nEnd; n++) {
        G4ThreeVector translation;  // Identity.
        G4RotationMatrix rotation;  // Identity - life enough for visualizing.
        G4RotationMatrix* pRotation = 0;
        switch (axis) {
          default:
          case kXAxis:
            translation = G4ThreeVector (-width*(nReplicas-1)*0.5+n*width,0,0);
            break;
          case kYAxis:
            translation = G4ThreeVector (0,-width*(nReplicas-1)*0.5+n*width,0);
            break;
          case kZAxis:
            translation = G4ThreeVector (0,0,-width*(nReplicas-1)*0.5+n*width);
            break;
          case kRho:
            if (pSol->GetEntityType() == "G4Tubs") {
              ((G4Tubs*)pSol)->SetInnerRadius(width*n+offset);
              ((G4Tubs*)pSol)->SetOuterRadius(width*(n+1)+offset);
            } else {
              if (fpMP->IsWarning())
                G4warn <<
                "G4PhysicalVolumeModel::VisitGeometryAndGetVisReps: WARNING:"
                "\n  built-in replicated volumes replicated in radius for "
                << pSol->GetEntityType() <<
                "-type\n  solids (your solid \""
                << pSol->GetName() <<
                "\") are not visualisable."
                << G4endl;
              visualisable = false;
            }
            break;
          case kPhi:
            rotation.rotateZ (-(offset+(n+0.5)*width));
            // Minus Sign because for the physical volume we need the
            // coordinate system rotation.
            pRotation = &rotation;
            break;
        }
        pVPV -> SetTranslation (translation);
        pVPV -> SetRotation    (pRotation);
        pVPV -> SetCopyNo (n);
        fCurrentPVCopyNo = n;
        if (visualisable) {
          DescribeAndDescend (pVPV, requestedDepth, pLV, pSol, pMaterial,
                              theAT, sceneHandler);
        }
      }
      // Restore originals...
      pVPV -> SetTranslation (originalTranslation);
      pVPV -> SetRotation    (pOriginalRotation);
      if (axis == kRho && pSol->GetEntityType() == "G4Tubs") {
        ((G4Tubs*)pSol)->SetInnerRadius(originalRMin);
        ((G4Tubs*)pSol)->SetOuterRadius(originalRMax);
      }
    }
  }
}

Referenced by DescribeAndDescend(), and DescribeYourselfTo().

Member Data Documentation

fAbort

G4bool G4PhysicalVolumeModel::fAbort

mutableprotected

Definition at line 322 of file G4PhysicalVolumeModel.hh.

Referenced by Abort(), DescribeAndDescend(), DescribeYourselfTo(), and G4PhysicalVolumeModel().

fBaseFullPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fBaseFullPVPath

protected

Definition at line 319 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetBaseFullPVPath().

fClippingMode

ClippingMode G4PhysicalVolumeModel::fClippingMode

protected

Definition at line 325 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeSolid(), G4PhysicalVolumeModel(), and SetClippingMode().

fCurrentDepth

G4int G4PhysicalVolumeModel::fCurrentDepth

protected

Definition at line 313 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetCurrentDepth(), and VisitGeometryAndGetVisReps().

fCurrentPVCopyNo

G4int G4PhysicalVolumeModel::fCurrentPVCopyNo

protected

Definition at line 315 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetCurrentPVCopyNo(), and VisitGeometryAndGetVisReps().

fCurrentTransform

G4Transform3D G4PhysicalVolumeModel::fCurrentTransform

protected

Definition at line 318 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), G4PhysicalVolumeModel(), and GetCurrentTransform().

fCurtailDescent

G4bool G4PhysicalVolumeModel::fCurtailDescent

mutableprotected

Definition at line 323 of file G4PhysicalVolumeModel.hh.

Referenced by CurtailDescent(), DescribeAndDescend(), DescribeYourselfTo(), and G4PhysicalVolumeModel().

fDrawnPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fDrawnPVPath

protected

Definition at line 321 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), DescribeYourselfTo(), and GetDrawnPVPath().

fFullPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fFullPVPath

protected

Definition at line 320 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetFullPVPath(), and VisitGeometryAndGetVisReps().

fMapAllTouchables

std::map<G4int,G4int> G4PhysicalVolumeModel::fMapAllTouchables

protected

Definition at line 329 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), DescribeYourselfTo(), and GetMapOfAllTouchables().

fMapDrawnTouchables

std::map<G4int,G4int> G4PhysicalVolumeModel::fMapDrawnTouchables

protected

Definition at line 327 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), DescribeYourselfTo(), and GetMapOfDrawnTouchables().

fMaxFullDepth

G4int G4PhysicalVolumeModel::fMaxFullDepth

protected

Definition at line 331 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeAndDescend(), G4PhysicalVolumeModel(), and GetMaxFullDepth().

fNClippers

G4int G4PhysicalVolumeModel::fNClippers

protected

Definition at line 326 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeSolid(), DescribeYourselfTo(), and G4PhysicalVolumeModel().

fpClippingSolid

G4VSolid* G4PhysicalVolumeModel::fpClippingSolid

protected

Definition at line 324 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeSolid(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetClippingSolid(), SetClippingSolid(), and ~G4PhysicalVolumeModel().

fpCurrentLV

G4LogicalVolume* G4PhysicalVolumeModel::fpCurrentLV

protected

Definition at line 316 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetCurrentLV().

fpCurrentMaterial

G4Material* G4PhysicalVolumeModel::fpCurrentMaterial

protected

Definition at line 317 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetCurrentMaterial().

fpCurrentPV

G4VPhysicalVolume* G4PhysicalVolumeModel::fpCurrentPV

protected

Definition at line 314 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetCurrentPV(), and GetCurrentTag().

fpTopPV

G4VPhysicalVolume* G4PhysicalVolumeModel::fpTopPV

protected

Definition at line 306 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent(), G4LogicalVolumeModel::DescribeYourselfTo(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetTopPhysicalVolume(), and Validate().

fRequestedDepth

G4int G4PhysicalVolumeModel::fRequestedDepth

protected

Definition at line 309 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent(), DescribeYourselfTo(), G4PhysicalVolumeModel(), GetRequestedDepth(), and SetRequestedDepth().

fTopPVCopyNo

G4int G4PhysicalVolumeModel::fTopPVCopyNo

protected

Definition at line 308 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent(), G4PhysicalVolumeModel(), and VisitGeometryAndGetVisReps().

fTopPVName

G4String G4PhysicalVolumeModel::fTopPVName

protected

Definition at line 307 of file G4PhysicalVolumeModel.hh.

Referenced by G4PhysicalVolumeModel().

fTotalAllTouchables

G4int G4PhysicalVolumeModel::fTotalAllTouchables

protected

Definition at line 330 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetTotalAllTouchables().

fTotalDrawnTouchables

G4int G4PhysicalVolumeModel::fTotalDrawnTouchables

protected

Definition at line 328 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetTotalDrawnTouchables().

fTransform

G4Transform3D G4PhysicalVolumeModel::fTransform

protected

Definition at line 312 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent(), DescribeYourselfTo(), G4PhysicalVolumeModel(), and GetTransformation().

fUseFullExtent

G4bool G4PhysicalVolumeModel::fUseFullExtent

protected

Definition at line 311 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent(), and G4PhysicalVolumeModel().

---

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

• G4PhysicalVolumeModel.hh
• G4PhysicalVolumeModel.cc