G4PhysicalVolumeModel Class Reference

#include <G4PhysicalVolumeModel.hh>

Inheritance diagram for G4PhysicalVolumeModel:

Classes
class	G4PhysicalVolumeModelTouchable
class	G4PhysicalVolumeNodeID

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)
virtual	~G4PhysicalVolumeModel ()
void	DescribeYourselfTo (G4VGraphicsScene &)
G4String	GetCurrentDescription () const
G4String	GetCurrentTag () const
G4VPhysicalVolume *	GetTopPhysicalVolume () const
G4int	GetRequestedDepth () const
const G4VSolid *	GetClippingSolid () const
G4int	GetCurrentDepth () const
G4VPhysicalVolume *	GetCurrentPV () const
G4LogicalVolume *	GetCurrentLV () const
G4Material *	GetCurrentMaterial () 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
void	SetBaseFullPVPath (const std::vector< G4PhysicalVolumeNodeID > &baseFullPVPath)
void	SetRequestedDepth (G4int requestedDepth)
void	SetClippingSolid (G4VSolid *pClippingSolid)
void	SetClippingMode (ClippingMode mode)
G4bool	Validate (G4bool warn)
void	CurtailDescent ()
Public Member Functions inherited from G4VModel
	G4VModel (const G4Transform3D &modelTransformation=G4Transform3D(), const G4ModelingParameters *=0)
virtual	~G4VModel ()
virtual void	DescribeYourselfTo (G4VGraphicsScene &)=0
const G4ModelingParameters *	GetModelingParameters () const
const G4String &	GetType () const
virtual G4String	GetCurrentDescription () const
virtual G4String	GetCurrentTag () const
const G4VisExtent &	GetExtent () const
const G4String &	GetGlobalDescription () const
const G4String &	GetGlobalTag () const
const G4Transform3D &	GetTransformation () const
void	SetModelingParameters (const G4ModelingParameters *)
void	SetExtent (const G4VisExtent &)
void	SetType (const G4String &)
void	SetGlobalDescription (const G4String &)
void	SetGlobalTag (const G4String &)
void	SetTransformation (const G4Transform3D &)
virtual G4bool	Validate (G4bool warn=true)

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)
void	CalculateExtent ()

Protected Attributes
G4VPhysicalVolume *	fpTopPV
G4String	fTopPVName
G4int	fTopPVCopyNo
G4int	fRequestedDepth
G4bool	fUseFullExtent
G4int	fCurrentDepth
G4VPhysicalVolume *	fpCurrentPV
G4LogicalVolume *	fpCurrentLV
G4Material *	fpCurrentMaterial
G4Transform3D *	fpCurrentTransform
std::vector< G4PhysicalVolumeNodeID >	fBaseFullPVPath
std::vector< G4PhysicalVolumeNodeID >	fFullPVPath
std::vector< G4PhysicalVolumeNodeID >	fDrawnPVPath
G4bool	fCurtailDescent
G4VSolid *	fpClippingSolid
ClippingMode	fClippingMode
Protected Attributes inherited from G4VModel
G4String	fType
G4String	fGlobalTag
G4String	fGlobalDescription
G4VisExtent	fExtent
G4Transform3D	fTransform
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
	)

Definition at line 58 of file G4PhysicalVolumeModel.cc.

                      :
  G4VModel        (modelTransformation, pMP),
  fpTopPV         (pVPV),
  fRequestedDepth (requestedDepth),
  fUseFullExtent  (useFullExtent),
  fCurrentDepth   (0),
  fpCurrentPV     (0),
  fpCurrentLV     (0),
  fpCurrentMaterial (0),
  fpCurrentTransform (0),
  fCurtailDescent (false),
  fpClippingSolid (0),
  fClippingMode   (subtraction)
{
  if (fpTopPV) {
    
    fType = "G4PhysicalVolumeModel";
    fTopPVName = fpTopPV -> GetName ();
    fTopPVCopyNo = fpTopPV -> GetCopyNo ();
    std::ostringstream o;
    o << fpTopPV -> GetCopyNo ();
    fGlobalTag = fpTopPV -> GetName () + "." + o.str();
    fGlobalDescription = "G4PhysicalVolumeModel " + fGlobalTag;
 
    CalculateExtent ();
  }
}

~G4PhysicalVolumeModel()

G4PhysicalVolumeModel::~G4PhysicalVolumeModel ( )

virtual

Definition at line 91 of file G4PhysicalVolumeModel.cc.

{
  delete fpClippingSolid;
}

Member Function Documentation

CalculateExtent()

void G4PhysicalVolumeModel::CalculateExtent ( )

protected

Definition at line 96 of file G4PhysicalVolumeModel.cc.

{
  if (fUseFullExtent) {
    fExtent = fpTopPV -> GetLogicalVolume () -> GetSolid () -> GetExtent ();
  }
  else {
    G4BoundingSphereScene bsScene(this);
    const G4int tempRequestedDepth = fRequestedDepth;
    fRequestedDepth = -1;  // Always search to all depths to define extent.
    const G4ModelingParameters* tempMP = fpMP;
    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).
    fpMP = &mParams;
    DescribeYourselfTo (bsScene);
    G4double radius = bsScene.GetRadius();
    if (radius < 0.) {  // Nothing in the scene.
      fExtent = fpTopPV -> GetLogicalVolume () -> GetSolid () -> GetExtent ();
    } else {
      // Transform back to coordinates relative to the top
      // transformation, which is in G4VModel::fTransform.  This makes
      // it conform to all models, which are defined by a
      // transformation and an extent relative to that
      // transformation...
      G4Point3D centre = bsScene.GetCentre();
      centre.transform(fTransform.inverse());
      fExtent = G4VisExtent(centre, radius);
    }
    fpMP = tempMP;
    fRequestedDepth = tempRequestedDepth;
  }
}

Referenced by G4PhysicalVolumeModel(), and Validate().

CreateCurrentAttValues()

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

(

)

const

Definition at line 829 of file G4PhysicalVolumeModel.cc.

{
  std::vector<G4AttValue>* values = new std::vector<G4AttValue>;
  std::ostringstream oss;
  for (size_t i = 0; i < fFullPVPath.size(); ++i) {
    oss << fFullPVPath[i].GetPhysicalVolume()->GetName()
    << ':' << fFullPVPath[i].GetCopyNo();
    if (i != fFullPVPath.size() - 1) oss << '/';
  }
 
  if (!fpCurrentLV) {
     G4Exception
        ("G4PhysicalVolumeModel::CreateCurrentAttValues",
         "modeling0004",
         JustWarning,
         "Current logical volume not defined.");
     return values;
  }
 
  values->push_back(G4AttValue("PVPath", 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(),""));
  oss.str(""); oss << '\n' << *fpCurrentTransform;
  values->push_back(G4AttValue("Trans", 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(), and G4XXXStoredSceneHandler::PreAddSolid().

CurtailDescent()

void G4PhysicalVolumeModel::CurtailDescent ( )

inline

Definition at line 222 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 334 of file G4PhysicalVolumeModel.cc.

{
  // Maintain useful data members...
  fpCurrentPV = pVPV;
  fpCurrentLV = pLV;
  fpCurrentMaterial = pMaterial;
 
  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;
  fpCurrentTransform = &theNewAT;
 
  const G4VisAttributes* pVisAttribs = pLV->GetVisAttributes();
  if (!pVisAttribs) pVisAttribs = fpMP->GetDefaultVisAttributes();
  // Beware - pVisAttribs might still be zero - create a temporary default one...
  G4bool visAttsCreated = false;
  if (!pVisAttribs) {
    pVisAttribs = new G4VisAttributes;
    visAttsCreated = true;
  }
  // From here, can assume pVisAttribs is a valid pointer.
  
  // Make decision to draw...
  G4bool thisToBeDrawn = true;
  
  // Update full path of physical volumes...
  G4int copyNo = fpCurrentPV->GetCopyNo();
  fFullPVPath.push_back
  (G4PhysicalVolumeNodeID
   (fpCurrentPV,copyNo,fCurrentDepth,*fpCurrentTransform));
  
  // In case we need to copy the vis atts for modification...
  G4bool copyForVAM = false;
  const G4VisAttributes* pUnmodifiedVisAtts = pVisAttribs;
  G4VisAttributes* pModifiedVisAtts = 0;
  
  // Check if vis attributes are to be modified by a /vis/touchable/set/ command.
  const std::vector<G4ModelingParameters::VisAttributesModifier>& vams =
    fpMP->GetVisAttributesModifiers();
  std::vector<G4ModelingParameters::VisAttributesModifier>::const_iterator
    iModifier;
  for (iModifier = vams.begin();
       iModifier != vams.end();
       ++iModifier) {
    const G4ModelingParameters::PVNameCopyNoPath& vamPath =
      iModifier->GetPVNameCopyNoPath();
    if (vamPath.size() == fFullPVPath.size()) {
      // OK - there's a size match.  Check it out.
//      G4cout << "Size match" << G4endl;
      G4ModelingParameters::PVNameCopyNoPathConstIterator iVAMNameCopyNo;
      std::vector<G4PhysicalVolumeNodeID>::const_iterator iPVNodeId;
      for (iVAMNameCopyNo = vamPath.begin(), iPVNodeId = fFullPVPath.begin();
           iVAMNameCopyNo != vamPath.end();
           ++iVAMNameCopyNo, ++iPVNodeId) {
//        G4cout
//        << iVAMNameCopyNo->fName
//        << ',' << iVAMNameCopyNo->fCopyNo
//        << "; " << iPVNodeId->GetPhysicalVolume()->GetName()
//        << ','  << iPVNodeId->GetPhysicalVolume()->GetCopyNo()
//        << G4endl;
        if (!(
              iVAMNameCopyNo->GetName() ==
              iPVNodeId->GetPhysicalVolume()->GetName() &&
              iVAMNameCopyNo->GetCopyNo() ==
              iPVNodeId->GetPhysicalVolume()->GetCopyNo()
              )) {
          break;
        }
      }
      if (iVAMNameCopyNo == vamPath.end()) {
//        G4cout << "Match found" << G4endl;
        if (!copyForVAM) {
          pModifiedVisAtts = new G4VisAttributes(*pUnmodifiedVisAtts);
          pVisAttribs = pModifiedVisAtts;
          copyForVAM = true;
        }
        const G4VisAttributes& transVisAtts = iModifier->GetVisAttributes();
        switch (iModifier->GetVisAttributesSignifier()) {
          case G4ModelingParameters::VASVisibility:
            pModifiedVisAtts->SetVisibility(transVisAtts.IsVisible());
            break;
          case G4ModelingParameters::VASDaughtersInvisible:
            pModifiedVisAtts->SetDaughtersInvisible
            (transVisAtts.IsDaughtersInvisible());
            break;
          case G4ModelingParameters::VASColour:
            pModifiedVisAtts->SetColour(transVisAtts.GetColour());
            break;
          case G4ModelingParameters::VASLineStyle:
            pModifiedVisAtts->SetLineStyle(transVisAtts.GetLineStyle());
            break;
          case G4ModelingParameters::VASLineWidth:
            pModifiedVisAtts->SetLineWidth(transVisAtts.GetLineWidth());
            break;
          case G4ModelingParameters::VASForceWireframe:
            if (transVisAtts.GetForcedDrawingStyle() ==
                G4VisAttributes::wireframe) {
              pModifiedVisAtts->SetForceWireframe
              (transVisAtts.IsForceDrawingStyle());
            }
            break;
          case G4ModelingParameters::VASForceSolid:
            if (transVisAtts.GetForcedDrawingStyle() ==
                G4VisAttributes::solid) {
              pModifiedVisAtts->SetForceSolid
              (transVisAtts.IsForceDrawingStyle());
            }
            break;
          case G4ModelingParameters::VASForceAuxEdgeVisible:
            pModifiedVisAtts->SetForceAuxEdgeVisible
            (transVisAtts.IsForceAuxEdgeVisible());
            break;
          case G4ModelingParameters::VASForceLineSegmentsPerCircle:
            pModifiedVisAtts->SetForceLineSegmentsPerCircle
            (transVisAtts.GetForcedLineSegmentsPerCircle());
            break;
        }
      }
    }
  }
 
  // There are various reasons why this volume
  // might not be drawn...
  G4bool culling = fpMP->IsCulling();
  G4bool cullingInvisible = fpMP->IsCullingInvisible();
  G4bool markedVisible = pVisAttribs->IsVisible();
  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;
    }
  }
  
  // Record thisToBeDrawn in path...
  fFullPVPath.back().SetDrawn(thisToBeDrawn);
 
  if (thisToBeDrawn) {
 
    // Update path of drawn physical volumes...
    fDrawnPVPath.push_back
      (G4PhysicalVolumeNodeID
       (fpCurrentPV,copyNo,fCurrentDepth,*fpCurrentTransform,thisToBeDrawn));
 
    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 = 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 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.;
    // 4) Global culling is on....
    if (culling) {
      // 5) ..and culling of invisible volumes is on...
      if (cullingInvisible) {
    // 6) ...and the mother requests daughters invisible
    if (daughtersInvisible) daughtersToBeDrawn = false;
      }
      // 7) Or culling of covered daughters is requested...
      if (cullingCovered) {
    // 8) ...and surface drawing is operating...
    if (surfaceDrawing) {
      // 9) ...but only if mother is visible...
      if (thisToBeDrawn) {
        // 10) ...and opaque...
          if (opaque) daughtersToBeDrawn = false;
      }
    }
      }
    }
  }
 
  // Delete modified vis atts if created...
  if (copyForVAM) {
    delete pModifiedVisAtts;
    pVisAttribs = pUnmodifiedVisAtts;
  }
  
  // Vis atts for this volume no longer needed if created...
  if (visAttsCreated) delete pVisAttribs;
 
  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--;
    }
  }
 
  // 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 608 of file G4PhysicalVolumeModel.cc.

{
  sceneHandler.PreAddSolid (theAT, *pVisAttribs);
 
  G4VSolid* pSectionSolid = fpMP->GetSectionSolid();
  G4VSolid* pCutawaySolid = fpMP->GetCutawaySolid();
 
  if (!fpClippingSolid && !pSectionSolid && !pCutawaySolid) {
 
    pSol -> DescribeYourselfTo (sceneHandler);  // Standard treatment.
 
  } else {
 
    // Clipping, etc., performed by Boolean operations.
 
    // First, get polyhedron for current solid...
    if (pVisAttribs->IsForceLineSegmentsPerCircle())
      G4Polyhedron::SetNumberOfRotationSteps
    (pVisAttribs->GetForcedLineSegmentsPerCircle());
    else
      G4Polyhedron::SetNumberOfRotationSteps(fpMP->GetNoOfSides());
    const G4Polyhedron* pOriginal = pSol->GetPolyhedron();
    G4Polyhedron::ResetNumberOfRotationSteps();
 
    if (!pOriginal) {
 
      if (fpMP->IsWarning())
    G4cout <<
      "WARNING: G4PhysicalVolumeModel::DescribeSolid: solid\n  \""
           << pSol->GetName() <<
      "\" has no polyhedron.  Cannot by clipped."
           << G4endl;
      pSol -> DescribeYourselfTo (sceneHandler);  // Standard treatment.
 
    } else {
 
      G4Polyhedron resultant(*pOriginal);
      G4VisAttributes resultantVisAttribs(*pVisAttribs);
      G4VSolid* resultantSolid = 0;
 
      if (fpClippingSolid) {
    switch (fClippingMode) {
    default:
    case subtraction:
      resultantSolid = new G4SubtractionSolid
        ("resultant_solid", pSol, fpClippingSolid, theAT.inverse());
      break;
    case intersection:
      resultantSolid = new G4IntersectionSolid
        ("resultant_solid", pSol, fpClippingSolid, theAT.inverse());
      break;
    }
      }
 
      if (pSectionSolid) {
    resultantSolid = new G4IntersectionSolid
      ("sectioned_solid", pSol, pSectionSolid, theAT.inverse());
      }
 
      if (pCutawaySolid) {
    resultantSolid = new G4SubtractionSolid
      ("cutaway_solid", pSol, pCutawaySolid, theAT.inverse());
      }
 
      G4Polyhedron* tmpResultant = resultantSolid->GetPolyhedron();
      if (tmpResultant) resultant = *tmpResultant;
      else {
    if (fpMP->IsWarning())
      G4cout <<
        "WARNING: G4PhysicalVolumeModel::DescribeSolid: resultant polyhedron for"
        "\n  solid \"" << pSol->GetName() <<
        "\" not defined due to error during Boolean processing."
        "\n  Original will be drawn in red."
         << G4endl;
    resultantVisAttribs.SetColour(G4Colour::Red());
      }
 
      delete resultantSolid;
 
      // Finally, force polyhedron drawing...
      resultant.SetVisAttributes(resultantVisAttribs);
      sceneHandler.BeginPrimitives(theAT);
      sceneHandler.AddPrimitive(resultant);
      sceneHandler.EndPrimitives();
    }
  }
  sceneHandler.PostAddSolid ();
}

Referenced by DescribeAndDescend().

DescribeYourselfTo()

void G4PhysicalVolumeModel::DescribeYourselfTo ( G4VGraphicsScene &  sceneHandler )

virtual

Implements G4VModel.

Definition at line 135 of file G4PhysicalVolumeModel.cc.

{
  if (!fpTopPV) G4Exception
    ("G4PhysicalVolumeModel::DescribeYourselfTo",
     "modeling0012", FatalException, "No model.");
 
  if (!fpMP) G4Exception
    ("G4PhysicalVolumeModel::DescribeYourselfTo",
     "modeling0003", FatalException, "No modeling parameters.");
 
  // For safety...
  fCurrentDepth = 0;
 
  G4Transform3D startingTransformation = fTransform;
 
  VisitGeometryAndGetVisReps
    (fpTopPV,
     fRequestedDepth,
     startingTransformation,
     sceneHandler);
 
  // Clear data...
  fCurrentDepth     = 0;
  fpCurrentPV       = 0;
  fpCurrentLV       = 0;
  fpCurrentMaterial = 0;
  if (fFullPVPath.size() != fBaseFullPVPath.size()) {
    // They should be equal if pushing and popping is happening properly.
    G4Exception
    ("G4PhysicalVolumeModel::DescribeYourselfTo",
     "modeling0013",
     FatalException,
     "Path at start of modeling not equal to base path.  Something badly"
     "\nwrong.  Please contact visualisation coordinator.");
  }
  fDrawnPVPath.clear();
}

Referenced by CalculateExtent(), DescribeSolid(), G4LogicalVolumeModel::DescribeYourselfTo(), G4ASCIITreeSceneHandler::EndModeling(), G4VisCommandSceneAddVolume::SetNewValue(), and Validate().

GetAttDefs()

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

(

)

const

Definition at line 754 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)["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)["Trans"] =
    G4AttDef("Trans","Transformation 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(), G4XXXStoredSceneHandler::PreAddSolid(), and G4VisCommandList::SetNewValue().

GetClippingSolid()

const G4VSolid * G4PhysicalVolumeModel::GetClippingSolid ( ) const

inline

Definition at line 158 of file G4PhysicalVolumeModel.hh.

  {return fpClippingSolid;}

GetCurrentDepth()

G4int G4PhysicalVolumeModel::GetCurrentDepth ( ) const

inline

Definition at line 161 of file G4PhysicalVolumeModel.hh.

{return fCurrentDepth;}

Referenced by G4HepRepSceneHandler::AddPrimitive(), G4GMocrenFileSceneHandler::AddSolid(), and G4HepRepSceneHandler::AddSolid().

GetCurrentDescription()

G4String G4PhysicalVolumeModel::GetCurrentDescription ( ) const

virtual

Reimplemented from G4VModel.

Definition at line 186 of file G4PhysicalVolumeModel.cc.

{
  return "G4PhysicalVolumeModel " + GetCurrentTag ();
}

GetCurrentLV()

G4LogicalVolume * G4PhysicalVolumeModel::GetCurrentLV ( ) const

inline

Definition at line 167 of file G4PhysicalVolumeModel.hh.

{return fpCurrentLV;}

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

GetCurrentMaterial()

G4Material * G4PhysicalVolumeModel::GetCurrentMaterial ( ) const

inline

Definition at line 170 of file G4PhysicalVolumeModel.hh.

{return fpCurrentMaterial;}

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

GetCurrentPV()

G4VPhysicalVolume * G4PhysicalVolumeModel::GetCurrentPV ( ) const

inline

Definition at line 164 of file G4PhysicalVolumeModel.hh.

{return fpCurrentPV;}

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

GetCurrentTag()

G4String G4PhysicalVolumeModel::GetCurrentTag ( ) const

virtual

Reimplemented from G4VModel.

Definition at line 174 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().

GetDrawnPVPath()

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

inline

Definition at line 180 of file G4PhysicalVolumeModel.hh.

  {return fDrawnPVPath;}

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

GetFullPVPath()

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

inline

Definition at line 173 of file G4PhysicalVolumeModel.hh.

  {return fFullPVPath;}

GetRequestedDepth()

G4int G4PhysicalVolumeModel::GetRequestedDepth ( ) const

inline

Definition at line 156 of file G4PhysicalVolumeModel.hh.

{return fRequestedDepth;}

Referenced by G4ASCIITreeSceneHandler::EndModeling().

GetTopPhysicalVolume()

G4VPhysicalVolume * G4PhysicalVolumeModel::GetTopPhysicalVolume ( ) const

inline

Definition at line 154 of file G4PhysicalVolumeModel.hh.

{return fpTopPV;}

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

SetBaseFullPVPath()

void G4PhysicalVolumeModel::SetBaseFullPVPath ( const std::vector< G4PhysicalVolumeNodeID > &  baseFullPVPath )

inline

Definition at line 200 of file G4PhysicalVolumeModel.hh.

                     {
      fBaseFullPVPath = baseFullPVPath;
      fFullPVPath = baseFullPVPath;
  }

Referenced by G4VisCommandSceneAddVolume::SetNewValue().

SetClippingMode()

void G4PhysicalVolumeModel::SetClippingMode ( ClippingMode  mode )

inline

Definition at line 215 of file G4PhysicalVolumeModel.hh.

                                           {
    fClippingMode = mode;
  }

SetClippingSolid()

void G4PhysicalVolumeModel::SetClippingSolid ( G4VSolid *  pClippingSolid )

inline

Definition at line 211 of file G4PhysicalVolumeModel.hh.

                                                   {
    fpClippingSolid = pClippingSolid;
  }

SetRequestedDepth()

void G4PhysicalVolumeModel::SetRequestedDepth ( G4int  requestedDepth )

inline

Definition at line 207 of file G4PhysicalVolumeModel.hh.

                                                {
    fRequestedDepth = requestedDepth;
  }

Validate()

G4bool G4PhysicalVolumeModel::Validate ( G4bool  warn )

virtual

Reimplemented from G4VModel.

Definition at line 701 of file G4PhysicalVolumeModel.cc.

{
  G4TransportationManager* transportationManager =
    G4TransportationManager::GetTransportationManager ();
 
  size_t nWorlds = transportationManager->GetNoWorlds();
 
  G4bool found = false;
 
  std::vector<G4VPhysicalVolume*>::iterator iterWorld =
    transportationManager->GetWorldsIterator();
  for (size_t i = 0; i < nWorlds; ++i, ++iterWorld) {
    G4VPhysicalVolume* world = (*iterWorld);
    // The idea now is to seek a PV with the same name and copy no
    // in the hope it's the same one!!
    G4PhysicalVolumeModel searchModel (world);
    G4PhysicalVolumeSearchScene searchScene
      (&searchModel, fTopPVName, fTopPVCopyNo);
    G4ModelingParameters mp;  // Default modeling parameters for this search.
    mp.SetDefaultVisAttributes(fpMP? fpMP->GetDefaultVisAttributes(): 0);
    searchModel.SetModelingParameters (&mp);
    searchModel.DescribeYourselfTo (searchScene);
    G4VPhysicalVolume* foundVolume = searchScene.GetFoundVolume ();
    if (foundVolume) {
      if (foundVolume != fpTopPV && warn) {
    G4cout <<
      "G4PhysicalVolumeModel::Validate(): A volume of the same name and"
      "\n  copy number (\""
           << fTopPVName << "\", copy " << fTopPVCopyNo
           << ") still exists and is being used."
      "\n  But it is not the same volume you originally specified"
      "\n  in /vis/scene/add/."
           << G4endl;
      }
      fpTopPV = foundVolume;
      CalculateExtent ();
      found = true;
    }
  }
  if (found) return true;
  else {
    if (warn) {
      G4cout <<
    "G4PhysicalVolumeModel::Validate(): No volume of name and"
    "\n  copy number (\""
         << fTopPVName << "\", copy " << fTopPVCopyNo
         << ") exists."
         << G4endl;
    }
    return false;
  }
}

VisitGeometryAndGetVisReps()

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

protected

Definition at line 191 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;
  G4Material* pMaterial;
 
  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);
    if (fCurrentDepth == 0) nReplicas = 1;  // Just draw first
    G4VPVParameterisation* pP = pVPV -> GetParameterisation ();
    if (pP) {  // Parametrised volume.
      for (int n = 0; n < nReplicas; n++) {
    pSol = pP -> ComputeSolid (n, pVPV);
    pP -> ComputeTransformation (n, pVPV);
    pSol -> ComputeDimensions (pP, n, pVPV);
    pVPV -> SetCopyNo (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 = 0; n < nReplicas; n++) {
    G4ThreeVector translation;  // Null.
    G4RotationMatrix rotation;  // Null - life long 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())
          G4cout <<
        "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);
    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);
      }
    }
  }
 
  return;
}

Referenced by DescribeAndDescend(), and DescribeYourselfTo().

Member Data Documentation

fBaseFullPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fBaseFullPVPath

protected

Definition at line 260 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeYourselfTo(), and SetBaseFullPVPath().

fClippingMode

ClippingMode G4PhysicalVolumeModel::fClippingMode

protected

Definition at line 265 of file G4PhysicalVolumeModel.hh.

Referenced by DescribeSolid(), and SetClippingMode().

fCurrentDepth

G4int G4PhysicalVolumeModel::fCurrentDepth

protected

Definition at line 255 of file G4PhysicalVolumeModel.hh.

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

fCurtailDescent

G4bool G4PhysicalVolumeModel::fCurtailDescent

protected

Definition at line 263 of file G4PhysicalVolumeModel.hh.

Referenced by CurtailDescent(), and DescribeAndDescend().

fDrawnPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fDrawnPVPath

protected

Definition at line 262 of file G4PhysicalVolumeModel.hh.

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

fFullPVPath

std::vector<G4PhysicalVolumeNodeID> G4PhysicalVolumeModel::fFullPVPath

protected

Definition at line 261 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), DescribeAndDescend(), DescribeYourselfTo(), GetFullPVPath(), G4PhysicalVolumeModel::G4PhysicalVolumeModelTouchable::GetReplicaNumber(), G4PhysicalVolumeModel::G4PhysicalVolumeModelTouchable::GetRotation(), G4PhysicalVolumeModel::G4PhysicalVolumeModelTouchable::GetSolid(), G4PhysicalVolumeModel::G4PhysicalVolumeModelTouchable::GetTranslation(), G4PhysicalVolumeModel::G4PhysicalVolumeModelTouchable::GetVolume(), SetBaseFullPVPath(), and VisitGeometryAndGetVisReps().

fpClippingSolid

G4VSolid* G4PhysicalVolumeModel::fpClippingSolid

protected

Definition at line 264 of file G4PhysicalVolumeModel.hh.

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

fpCurrentLV

G4LogicalVolume* G4PhysicalVolumeModel::fpCurrentLV

protected

Definition at line 257 of file G4PhysicalVolumeModel.hh.

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

fpCurrentMaterial

G4Material* G4PhysicalVolumeModel::fpCurrentMaterial

protected

Definition at line 258 of file G4PhysicalVolumeModel.hh.

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

fpCurrentPV

G4VPhysicalVolume* G4PhysicalVolumeModel::fpCurrentPV

protected

Definition at line 256 of file G4PhysicalVolumeModel.hh.

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

fpCurrentTransform

G4Transform3D* G4PhysicalVolumeModel::fpCurrentTransform

protected

Definition at line 259 of file G4PhysicalVolumeModel.hh.

Referenced by CreateCurrentAttValues(), and DescribeAndDescend().

fpTopPV

G4VPhysicalVolume* G4PhysicalVolumeModel::fpTopPV

protected

Definition at line 249 of file G4PhysicalVolumeModel.hh.

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

fRequestedDepth

G4int G4PhysicalVolumeModel::fRequestedDepth

protected

Definition at line 252 of file G4PhysicalVolumeModel.hh.

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

fTopPVCopyNo

G4int G4PhysicalVolumeModel::fTopPVCopyNo

protected

Definition at line 251 of file G4PhysicalVolumeModel.hh.

Referenced by G4PhysicalVolumeModel(), and Validate().

fTopPVName

G4String G4PhysicalVolumeModel::fTopPVName

protected

Definition at line 250 of file G4PhysicalVolumeModel.hh.

Referenced by G4PhysicalVolumeModel(), and Validate().

fUseFullExtent

G4bool G4PhysicalVolumeModel::fUseFullExtent

protected

Definition at line 254 of file G4PhysicalVolumeModel.hh.

Referenced by CalculateExtent().

---

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

• G4PhysicalVolumeModel.hh
• G4PhysicalVolumeModel.cc