A CAD Interface for GEANT4 - arXiv

1 1A CAD Interface for GEANT4 Christopher M Poole, Iwan Cornelius, Jamie V Trapp, Christian M LangtonAbstract Typically used as a tool for Monte Carlo simulationof high energy physics experiments, GEANT4 is increasinglybeing employed for the simulation of complex radiotherapytreatments. Often the specification of components within a clinicallinear accelerator treatment head is provided in a CAD fileformat. Direct import of these CAD files into GEANT4 may notbe possible, and complex components such as individual leaveswithin a multi-leaf collimator may be difficult to define via othermeans. Solutions that allow for users to work around the limitedsupport in the GEANT4 toolkit for loading predefined CADgeometries has been presented by others, however these solutionsrequire intermediate file format conversion using commercialsoftware.

2 Here within we describe a technique that allows forCAD models to be directly loaded as geometry without the needfor commercial software and intermediate file format of the Interface was tested using a set of CAD modelsof various complexity; for the models used in testing, no importerrors were reported and all geometry was found to be navigableby Terms Monte Carlo, GEANT4 , BEAMnrc, CAD, ge-ometry, radiotherapyI. INTRODUCTIONGEOMETRY and Tracking ( GEANT4 ) is a C++ toolkitspecifically designed to track particles traversing a ge-ometry whilst being subject to physical processes, it findsapplication in fields such as nuclear and particle physics andspace engineering, with increasing use in medical physics[2], [6], [10], [11], [14].

3 Numerous physical processes canbe modeled including photo-nuclear interactions, optical pro-cesses such as scintillation and Cherenkov radiation and otherparticle interactions over a wide energy range (250eVuptoT eVenergies); the full gamut of processes available tothe user is described by others [2]. Additionally, the toolkitprovides functionality for the inclusion and exclusion of thedesired processes as well as sufficient extensibility to includedcustom or user defined processes. Fast and effective geometrydefinition is also available to the user for relatively simpleobjects using constructs such asG4 Orbfor defining orbsor spheres,G4 Boxfor defining rectangular prisms and theconcept of boolean solids [2].

4 At the most fundamental level,the GEANT4 geometry hierarchy is divided into solids, logicalvolumes and physical volumes where solids described shape,logical volumes define material properties and mother daughterrelations, and physical volumes define placement within themother Description Mark-up Language (GDML) is acomprehensive geometry description format using ExtensibleCorrespondence can be directed authors are with the Discipline of Physics, Faculty of Science and Tech-nology, and the Institute of Health and Biomedical Innovation at QueenslandUniversity of Technology, Brisbane, work is funded by the Queensland Cancer Physics Collaborative, andCancer Australia (Department of Health and Ageing) Research Grant Language (XML) that allows for the persistence ofmany aspects of a geometry, including material properties andassemblies [7].

5 Functionality within the toolkit readily allowsfor geometries to be saved and reloaded as GDML, howeverthere is little support for converting pre-existing user definedComputer aided design (CAD) models to GDML or indeeddirectly loading these same CAD models as geometry [8]. Ofthe two methods available to the user for accomplishing thistask, both are reliant of intermediate file format conversionusing commercial software [5], [8].Many CAD packages export files compliant to the standardfor the exchange of product model data (STEP - ISO 10303),a standard designed to supersede the still widely used initialgraphics exchange specification (IGES) [4], [13]. As such, ISOstandard STEP has been the target format for loading CADgeometry into GEANT4 .

6 Persistence example, with identifier G02 , distributed with GEANT4 describes loading STEPT ools (STEP Tools Incorporated, New York) files directly, withthe intermediate conversion of STEP to STEP Tools usingthe commercial ST-Viewer program or ST-Developer libraries(STEP Tools Incorporated, New York), refer to (a). Thisprocess allows for assemblies of components to be loadeddirectly, however the STEP Tools programs and librariesmay be prohibitively expensive for some users. Constantineet al described the process of converting STEP to GDML using the commercial software FastRad (Tests & Radiation- Toulouse) [5], refer to (b). FastRad again may beconsidered prohibitively expensive for some users with therequirement of annual licensing, in addition to this, the trialversion limits the conversion process to no more than 20elements per , CAD file formats such as stereo-lithographyformat (STL) or Stanford polygon file format (PLY) de-scribing a 3-dimensional volume use tessellated polygonmeshes (usually triangular or quadrangular) to define itsclosed surface.

7 A cloud of points in 3-dimensional spacedefine the vertexes of the mesh and a collection of facesdefine the interconnects between the vertexes. The equiva-lent geometric construct in GEANT4 to this scheme is thetessellated solid (G4 TessellatedSolid); except for thetwo previously mentioned techniques, mappings between theG4 TessellatedSolidconstruct and CAD file formatsare not readily available [8]. Programming toolkits however,already exist for the creation and manipulation of triangularand quadrangular meshes; the kind exportable by many mod-ern CAD packages [1], [3]. Most notably the templated C++library VCGLIB (Visual Computing Laboratory, Italy) offersadvanced mesh manipulation routines and point cloud surfacereconstruction algorithms.

8 In addition to this, VCGLIB offersan import/export Interface for many common CAD file formats[1]. [ ] 5 May 20112 STEPCADSTEP ToolGEANT4exportST-viewerSTEP Tool Interface (a)STEPCADGDMLGEANT4exportFastR ADGDML Interface (b)PLY, OFF, Interface (c)Fig. 1: A diagrammatic comparison between two currently available CAD import techniques where (a) shows the techniquedescribed in persistence example G02, (b) shows the STEP to GDML conversion technique [8], and (c) shows the new directimport technique proposed here within we describe a very simple technique that allowsfor CAD models to be directly loaded as geometry in GEANT4using VCGLIB. This technique is capable of handling manyfile formates exportable by CAD programs whilst not relyingon intermediate file format METHODSA.

9 GEANT4 CAD InterfaceThree C++ classes were created describing a vertex, faceand a triangular mesh. The classCADV ertexof base typevcg::Vertexwas used to store a single vertex definedby three double precision values(x, y, z). For faces, theclassCADF aceof base typevcg::Facewas used to storea reference to all vertices that define the face(a, b, c).Finally, the classCADM eshtemplated on twostd::vectorcontainers was created. The first vector container was tem-plated on the custom typeCADV ertex, and the second onthe custom typeCADFace- in effect theCADM eshtypecontained two collections, one describing all the vertices ina mesh and another describing all of the interconnect betweenthe vertexes (the faces).

10 Through an instance of aCADMesh,the description of a triangular mesh with a cloud of verticesinterconnected by a group of faces could be in VCGLIB extend the VCGLIB core to pro-vide additional mesh manipulation functionality including im-porters and exporters. Using theOpenmethod in the wrappervcg::tri::io::ImporterPLYand parsing to it a PLYfile and a reference to an instance ofCADMesh, the instanceofCADM eshwas populated with a 3-dimensional triangularfacet mesh described the PLY file. Other triangular facetmesh importers include OBJ, OFF, and STL file formats [1].Random access to faces in the mesh was available usingtheCADMesh::FaceIteratoriterator where each faceprovided a pointer to each vertex(a, b, c)of the current faceand each vertex provided pointers to its coordinates(x, y, z).