Representing Range Compensators with Computational Geometry in TOPAS

1. Representing Range Compensators with Computational Geometry in TOPAS Forrest Iandola1,2 and Joseph Perl1 1 SLAC National Accelerator Laboratory 2 University of Illinois at Urbana-Champaign

2. Overview Medical Physics in 30 Seconds Introduction to TOPAS What is a range compensator? Subtraction Solid geometry Modeling compensators with Subtraction Solids Approximation with polyhedrons for (potentially) faster performance Forrest Iandola Computational Geometry Compensators

3. Medical Physics in 30 Seconds • Goal: kill cancer with radiation • Deliver radiation with protons, photons, other particles, or ions • Monte Carlo simulation of proton therapy beams is an up-and-coming field Forrest Iandola Computational Geometry Compensators

4. Introduction to TOPAS TOPAS (Tool for Particle Simulation) TOPAS = Monte Carlo simulation of radiation therapy beamlines User can easily customize beamline for specific treatment facilities Uses Geant4 for the “real” physics Forrest Iandola Computational Geometry Compensators

5. What is a Range Compensator? A radiation therapy beamline collimates the beam and produces a specific energy spread Range compensator produces an energy spread Construction: drill a number of holes out of a cylinder of lucite Each drill hole may have a unique depth “The thickness of the Lucite [plastic] will proportionally reduce the depth [energy] of the protons”1 1 http://neurosurgery.mgh.harvard.edu/protonbeam/ Forrest Iandola Computational Geometry Compensators

6. What is a Range Compensator? Forrest Iandola Computational Geometry Compensators

7. Subtraction Solids Geant4 supports boolean solid combinatorial geometry Subtraction solids Union solids It’s as simple as newSolid = Solid1 - Solid2 Overlap among subtracted solids is acceptable Solids can be recursively subtracted Forrest Iandola Computational Geometry Compensators

8. Compensator with Subtraction Solids Compensator comprised of a bigCylinder with n holes subtracted: newSolid1 = bigCylinder - smallCylinder1 newSolid2 = newSolid1 - smallCylinder2 … Compensator = newSolid(n-1) - smallCylinder(n) Forrest Iandola Computational Geometry Compensators

9. Compensator with Subtraction Solids Forrest Iandola Computational Geometry Compensators

10. Approximation for Performance Gains Approximate the drill holes with a collection of hexagons Lack of overlap among hexagons allows us to model all hexagons as a single polyhedron Future work: evaluate performance benefits (and accuracy reduction) with polyhedron method Forrest Iandola Computational Geometry Compensators

11. Approximation for Performance Gains Subtraction Solid Polyhedron Forrest Iandola Computational Geometry Compensators

12. Approximation for Performance Gains Future work: evaluate performance benefits (and accuracy reduction) with polyhedron method Forrest Iandola Computational Geometry Compensators

13. Acknowledgements Harald Paganetti (Massachusetts General Hospital and Harvard University) Jan Schuemann (Massachusetts General Hospital and Harvard University) Jungwook Shin (UC San Francisco) Bruce Faddegon (UC San Francisco) DOE and NIH for generous support Contact: forrest@slac.stanford.edu Forrest Iandola Computational Geometry Compensators