Loading Geometry

Generally, a reader class is provided for each format. The reader is created, then told to load a file, and it creates a Registry object (see Registry) containing the model. The registry is the final object from a reader, and its top-most volume can be used for visualisation or other operations.

GDML

In directory pyg4ometry/test/gdmlG4examples/ChargeExchangeMC/

import pyg4ometry

r = pyg4ometry.gdml.Reader("lht.gdml")
l = r.getRegistry().getWorldVolume()
v = pyg4ometry.visualisation.VtkViewerColouredMaterial()
v.addLogicalVolume(l)
v.view()

Geant4 example GDML from ChargeExchangeMC example.

FLUKA

In directory pyg4ometry/test/flukaCompoundExamples

Note FLUKA geometry can be loaded but cannot be visualised directly without conversion to a Geant4 model. This is not necessary for loading, but shown here.

import pyg4ometry

r = pyg4ometry.fluka.Reader("corrector-dipole2.inp")
flukaRegistry = r.flukaregistry

geantRegistry = pyg4ometry.convert.fluka2Geant4(flukaRegistry)

l = geantRegistry.getWorldVolume()
v = pyg4ometry.visualisation.VtkViewerColouredMaterial()
v.addLogicalVolume(l)
v.view()

FLUKA example of a dipole magnet.

ROOT

In directory pyg4ometry/test/root2Gdml

import pyg4ometry

r = pyg4ometry.io.ROOTTGeo.Reader("example.root")
l = r.getRegistry().getWorldVolume()
v = pyg4ometry.visualisation.VtkViewerColouredMaterial()
v.addLogicalVolume(l)
v.view()

ROOT example of Geant4’s LHT geometry.

STL

In directory pyg4ometry/test/stl

STL files are typically used for a single watertight solid mesh. This mesh is converted to a TesselatedSolid and then a logical volume which can be placed in a geometry. In directory pyg4ometry/test/stl.

import pyg4ometry

reg = pyg4ometry.geant4.Registry()
r = pyg4ometry.stl.Reader("utahteapot.stl", reg)
s = r.getSolid()
copper = pyg4ometry.geant4.MaterialPredefined("G4_Cu", reg)
l = pyg4ometry.geant4.LogicalVolume(s, copper, "utahteapot_lv", reg)
v = pyg4ometry.visualisation.VtkViewerColouredMaterial()
v.addLogicalVolume(l)
v.view()

alt:

Example of STL loading in pyg4ometry

STEP

In directory pyg4ometry/test/freecad

import pyg4ometry

r = pyg4ometry.freecad.Reader("08_AshTray.step")
r.relabelModel()
r.convertFlat()
l = r.getRegistry().getWorldVolume()
v = pyg4ometry.visualisation.VtkViewer()
v.addLogicalVolume(l)
v.view()