CUDA- based Geant4 Monte Carlo Simula8on for Radia8on Therapy. N. Henderson & K. Murakami GTC 2013
|
|
- Oliver Cook
- 5 years ago
- Views:
Transcription
1 CUDA- based Geant4 Monte Carlo Simula8on for Radia8on Therapy N. Henderson & K. Murakami GTC
2 The collabora8on Makoto Asai, SLAC Joseph Perl, SLAC Koichi Murakami, KEK- SLAC Takashi Sasaki, KEK Margot Gerritsen, ICME Nick Henderson, ICME Special thanks to the CUDA Center of Excellence Program 2
3 20,000,000 radiotherapy treatments per year in US 3
4 Image: capitalmedical.com 4
5 liver kidney Image: James M. Balter and Marc L. Kessler, JCO March 10, 2007 vol. 25 no
6 Simula8on methods Analy8c 8me: seconds to minutes accurate within 3-5% used in treatment planning Monte Carlo 8me: several hours to days of CPU 8me accurate within 1-2% used to verify treatment plans in certain cases 6
7 Geant4 Toolkit Enables Monte Carlo simula8on of par8cles travelling through and interac8ng with mayer Allows modeling of complex geometries Covers all elementary par8cles and nuclei for a wide energy range 7
8 Geant4 Applica8ons High Energy Physics Space & Radia8on Medical Physics ATLAS LISA gmocren Images from: Geant4 gallery and gmocren 8
9 Geant4 101 Geant4 simulates par8cles travelling through and interac8ng with mayer Example: photoelectric effect Image: CCSA by Wolfmankurd on Wikipedia 9
10 Paralleliza8on challenges in Geant4 Large and complex code base Sophis8cated geometry framework Elaborate physics models Branching, look- up tables, single- thread op8miza8ons 10
11 Thank goodness The simula8on is embarrassingly parallel! (the par8cles are independent) 11
12 Requirements for X- ray Radiotherapy Geometry is a voxelized box Physics is limited to low- energy electromagne8cs Material is modeled as water with different densi8es 12
13 Low energy electromagne8cs Gamma Compton scayering Photoelectric effect Gamma conversion Electron/positron Ioniza8on Bremsstrahlung Positron annihila8on Todo: Mul8ple scayering Image: CCSA by Man8corp on Wikipedia Image: public domain 13
14 What makes a physics process? Sample for angular and energy distribu8ons of secondary par8cles Deposit energy to material Produce secondary par8cles that must be tracked at a later point 14
15 Tracking algorithm Par8cles are tracked through space Each discrete move is called a step Physics process may occur along step, ager the step, or both step 1 par8cle secondary par8cle step 2 15
16 How are processes selected? Each process has an interac8on length (dx) Process with shortest dx is selected Ager the step dx is decreased or resampled 16
17 Energy deposi8on Happens along the step for the ioniza8on process May happen at the end of the step Secondary par8cles with too low energy are not generated, but treated as point- like energy deposi8ons step 1 par8cle secondary par8cle x step 2 17
18 G4CU: CUDA- base MC for RT Important data structures Algorithm summary Details Parallel stack Look- up tables 18
19 Struct- of- arrays data payern struct ParticleArray { // length of arrays int length; // kind of particle ParticleKind *kind; // position float *x, *y, *z; // direction float *dx, *dy, *dz; // particle energy float *energy; // voxel index int *vx, *vy, *vz, *vid; }; Common payern in CUDA to allow for coalesced memory access Experiments with transport showed this to be 3-4x faster than AOS 19
20 G4CU: algorithm query all processes to select step size reduce dose stacks to global array main loop apply all continuous processes decrease interaction lengths pop secondaries from stacks generate primary particles management apply limiting discrete process check termination conditions 20
21 Some details A process does a few things: 1. Changes direc8on and momentum of primary par8cle 2. Generates secondary par8cles 3. Deposits energy to the material We use thread local stacks to handle 2 and 3 21
22 Parallel stacks start end data pos template <typename T> struct pstack { // number of stacks int stack_num; // size of each stack int stack_size; // starts and ends of stacks int *start, *end; // stack positions int *pos; // stack data array T *data; }; 22
23 Dose reduc8on & storage Energy dose is stored in a thread local variable When par8cle moves to new voxel, dose is pushed onto stack Dose reduc8on is performed periodically We use Thrust index dose dose stack i join and sort dose stack j reduce by index
24 Used for interac8on length computa8on 40 bins, log spaced Linear and spline interpola8on Bremsstrahlung also uses 2D interpola8on Look- up tables 24
25 Benchmarks Configura8on: Geometry: 512 x 512 x 256 voxels Dose reduc8on frequency: every 200 itera8ons 128 blocks with 256 threads Primary par8cle is a 6 MeV gamma Tesla C
26 Simula8on 8me 100 million primary par8cles Time: 72 minutes ~ 23.1 primary par8cles per ms ~ 50-60x speedup over Geant4 on 1 CPU 26
27 Profile Component Percentage of overall 1me Physics processes 50 Energy dose reduc8on 30 Interac8on length 18 Run management 2 27
28 Physics process breakdown Process Par1cle Percentage of overall 1me Bremsstrahlung e- / e+ 23 Pair produc8on γ 7.5 Transport all 7 Photoelectric effect γ 7 Ioniza8on e- / e+ 3 Compton scayering γ 1.5 Positron annihila8on e+ 1 28
29 Dose distribu8on 29
30 Acknowledgements Geant4 Collabora8on, see Geant4 a simula1on toolkit, Nuclear Instruments and Methods in Physics Research A 506 (2003) Geant4 developments and Applica1ons, IEEE TransacGons on Nuclear Science 53 No. 1 (2006) NVIDIA & CCOE Program People: Koichi Murakami, KEK Makoto Asai, SLAC Joseph Perl, SLAC Takashi Sasaki, KEK Margot Gerritsen, ICME 30
Op#miza#on of CUDA- based Monte Carlo Simula#on for Radia#on Therapy. GTC 2014 N. Henderson & K. Murakami
Op#miza#on of CUDA- based Monte Carlo Simula#on for Radia#on Therapy GTC 2014 N. Henderson & K. Murakami The collabora#on Geant4 @ Special thanks to the CUDA Center of Excellence Program Makoto Asai, SLAC
More informationgpmc: GPU-Based Monte Carlo Dose Calculation for Proton Radiotherapy Xun Jia 8/7/2013
gpmc: GPU-Based Monte Carlo Dose Calculation for Proton Radiotherapy Xun Jia xunjia@ucsd.edu 8/7/2013 gpmc project Proton therapy dose calculation Pencil beam method Monte Carlo method gpmc project Started
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Monte Carlo Methods 1 Monte Carlo! Most accurate at predicting dose distributions! Based on
More informationThe Geant4-DNA project at the Physics-Biology frontier. JFY2015(BIO_02) reports and JFY2016(APP_01) plan Sébastien Incerti, CENBG Takashi Sasaki, KEK
The Geant4-DNA project at the Physics-Biology frontier JFY2015(BIO_02) reports and JFY2016(APP_01) plan Sébastien Incerti, CENBG Takashi Sasaki, KEK 2016 Project members Leader Sébastien Incerti Dr. CENBG/IN2P3
More informationMPEXS benchmark results
MPEXS benchmark results - phase space data - Akinori Kimura 14 February 2017 Aim To validate results of MPEXS with phase space data by comparing with Geant4 results Depth dose and lateral dose distributions
More informationA fast and accurate GPU-based proton transport Monte Carlo simulation for validating proton therapy treatment plans
A fast and accurate GPU-based proton transport Monte Carlo simulation for validating proton therapy treatment plans H. Wan Chan Tseung 1 J. Ma C. Beltran PTCOG 2014 13 June, Shanghai 1 wanchantseung.hok@mayo.edu
More informationElectron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations
Electron Dose Kernels (EDK) for Secondary Particle Transport in Deterministic Simulations A. Al-Basheer, G. Sjoden, M. Ghita Computational Medical Physics Team Nuclear & Radiological Engineering University
More informationImprovements in Monte Carlo simulation of large electron fields
SLAC-PUB-12908 October 2007 Improvements in Monte Carlo simulation of large electron fields Bruce A Faddegon 1, Joseph Perl 2 and Makoto Asai 2 1 University of California San Francisco Comprehensive Cancer
More informationELECTRON DOSE KERNELS TO ACCOUNT FOR SECONDARY PARTICLE TRANSPORT IN DETERMINISTIC SIMULATIONS
Computational Medical Physics Working Group Workshop II, Sep 30 Oct 3, 2007 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) ELECTRON DOSE
More informationMonte Carlo methods in proton beam radiation therapy. Harald Paganetti
Monte Carlo methods in proton beam radiation therapy Harald Paganetti Introduction: Proton Physics Electromagnetic energy loss of protons Distal distribution Dose [%] 120 100 80 60 40 p e p Ionization
More informationGPU-based high-performance computing for radiotherapy applications
GPU-based high-performance computing for radiotherapy applications Julien Bert, PhD CHRU de Brest LaTIM INSERM UMR1101 Radiotherapy Irradiation of the tumor: Maximum dose to the tumor Healthy surrounding
More informationSupercomputing the Cascade Processes of Radiation Transport
19 th World Conference on Non-Destructive Testing 2016 Supercomputing the Cascade Processes of Radiation Transport Mikhail ZHUKOVSKIY 1, Mikhail MARKOV 1, Sergey PODOLYAKO 1, Roman USKOV 1, Carsten BELLON
More informationRepresenting Range Compensators in the TOPAS Monte Carlo System
Representing Range Compensators in the TOPAS Monte Carlo System, Jan Schuemann, Jungwook Shin, Bruce Faddegon, Harald Paganetti, and Joseph Perl SLAC National Accelerator Laboratory University of Illinois
More informationStatus and plan for the hadron therapy simulation project in Japan
Status and plan for the hadron therapy simulation project in Japan Takashi Sasaki KEK and CREST, JST Takashi.Sasaki@kek.jp The Project The Development of Software Framework for Simulation in Radiotherapy
More informationGeant p01. Primary Par(cle. (based on slides by Makoto Asai) Geant4 Tutorial at Stanford 4 March 2014 Dennis Wright (SLAC)
Geant4 10.0 p01 Primary Par(cle (based on slides by Makoto Asai) Geant4 Tutorial at Stanford 4 March 2014 Dennis Wright (SLAC) Outline Primary Par(cle Genera(on Built- in primary par(cle generators par(cle
More informationOp#miza#on of a Fast CT System for More Accurate Radiotherapy. Erica Cherry, Meng Wu, and Rebecca Fahrig
Op#miza#on of a Fast CT System for More Accurate Radiotherapy Erica Cherry, Meng Wu, and Rebecca Fahrig Problem: mo#on during treatment Two overlapping thorax views during normal breathing liver lung spine
More informationGATE-RT Applications in Radiation Therapy
GATE-RT Applications in Radiation Therapy GATE Training, Orsay, April 2011 David SARRUT Creatis - CNRS Lyon - France GATE-RT: Radiation Therapy applications Related to dose distribution Radiotherapy: therapy
More informationDocumentation, examples and user support
Documentation, examples and user support Makoto Asai (SLAC PPA/SCA) November 12th, 2014 Geant4 tutorial @ ANS Winter Meeting 2014 Geant4 User Documenta0on Four manuals distributed at every Geant4 release
More informationPor$ng Monte Carlo Algorithms to the GPU. Ryan Bergmann UC Berkeley Serpent Users Group Mee$ng 9/20/2012 Madrid, Spain
Por$ng Monte Carlo Algorithms to the GPU Ryan Bergmann UC Berkeley Serpent Users Group Mee$ng 9/20/2012 Madrid, Spain 1 Outline Introduc$on to GPUs Why they are interes$ng How they operate Pros and cons
More informationAccelerating koblinger's method of compton scattering on GPU
Available online at www.sciencedirect.com Procedia Engineering 24 (211) 242 246 211 International Conference on Advances in Engineering Accelerating koblingers method of compton scattering on GPU Jing
More informationSoftware Suite for Particle Therapy Simulation
Software Suite for Particle Therapy Simulation Takashi Sasaki KEK Computing Research Center and CREST/JST Takashi.Sasaki@kek.jp The Project The Development of Software Framework for Simulation in Radiotherapy
More informationFrom Biological Cells to Populations of Individuals: Complex Systems Simulations with CUDA (S5133)
From Biological Cells to Populations of Individuals: Complex Systems Simulations with CUDA (S5133) Dr Paul Richmond Research Fellow University of Sheffield (NVIDIA CUDA Research Centre) Overview Complex
More informationMonte Carlo simulations. Lesson FYSKJM4710 Eirik Malinen
Monte Carlo simulations Lesson FYSKJM4710 Eirik Malinen MC simulations 1 Simulations of stochastic processes Interactions are stochastic: the path of a single ionizing particle may not be predicted Interactions
More informationMeasurement of depth-dose of linear accelerator and simulation by use of Geant4 computer code
reports of practical oncology and radiotherapy 1 5 (2 0 1 0) 64 68 available at www.sciencedirect.com journal homepage: http://www.rpor.eu/ Original article Measurement of depth-dose of linear accelerator
More informationARTCOLL PACKAGE FOR GAMMA ART6000 TM ROTATING GAMMA SYSTEM EMISSION SPECTRA CALCULATION
ARTCOLL PACKAGE FOR GAMMA ART6000 TM ROTATING GAMMA SYSTEM EMISSION SPECTRA CALCULATION Radovan D Ili}. PhD Institute of Nuclear Sciences Vin~a Belgrade, Serbia and Tomasz K. Helenowski, M.D. 936 Burnham
More informationGPU applications in Cancer Radiation Therapy at UCSD. Steve Jiang, UCSD Radiation Oncology Amit Majumdar, SDSC Dongju (DJ) Choi, SDSC
GPU applications in Cancer Radiation Therapy at UCSD Steve Jiang, UCSD Radiation Oncology Amit Majumdar, SDSC Dongju (DJ) Choi, SDSC Conventional Radiotherapy SIMULATION: Construciton, Dij Days PLANNING:
More informationOverview: Physically- Based Simulation of Single-Event Effects
Overview: Physically- Based Simulation of Single-Event Effects Robert A. Weller Institute for Space & Defense Electronics Vanderbilt University R. A. Reed M. H. Mendenhall D. R. Ball J. A. Pellish K. M.
More informationMonte Carlo simulations
MC simulations Monte Carlo simulations Eirik Malinen Simulations of stochastic processes Interactions are stochastic: the path of a single ioniing particle may not be predicted Interactions are quantified
More informationEvaluation of RayXpert for shielding design of medical facilities
Evaluation of Raypert for shielding design of medical facilities Sylvie Derreumaux 1,*, Sophie Vecchiola 1, Thomas Geoffray 2, and Cécile Etard 1 1 Institut for radiation protection and nuclear safety,
More informationProfiling & Tuning Applica1ons. CUDA Course July István Reguly
Profiling & Tuning Applica1ons CUDA Course July 21-25 István Reguly Introduc1on Why is my applica1on running slow? Work it out on paper Instrument code Profile it NVIDIA Visual Profiler Works with CUDA,
More informationFaster Navigation in Voxel Geometries (DICOM) Joseph Perl (SLAC/SCCS) G4NAMU AAPM Houston 27 July 2008
Faster Navigation in Voxel Geometries (DICOM) Joseph Perl (SLAC/SCCS) G4NAMU meeting @ AAPM Houston 27 July 2008 Technical Challenge #2 for Geant4 Med Apps: Patient Geometry HEP: constructive solid geometry
More informationimass: Computational NRF Spectra Signal from Geant4
imass: Computational NRF Spectra Signal from Geant4 John Perry, Shanjie Xiao, Tatjana Jevremovic School of Nuclear Engineering, Purdue University, West Lafayette, IN 47907 {joperry, xiaosj, tatjanaj} @purdue.edu
More informationLimitations in the PHOTON Monte Carlo gamma transport code
Nuclear Instruments and Methods in Physics Research A 480 (2002) 729 733 Limitations in the PHOTON Monte Carlo gamma transport code I. Orion a, L. Wielopolski b, * a St. Luke s/roosevelt Hospital, Columbia
More informationMathematical methods and simulations tools useful in medical radiation physics
Mathematical methods and simulations tools useful in medical radiation physics Michael Ljungberg, professor Department of Medical Radiation Physics Lund University SE-221 85 Lund, Sweden Major topic 1:
More informationMichael Speiser, Ph.D.
IMPROVED CT-BASED VOXEL PHANTOM GENERATION FOR MCNP MONTE CARLO Michael Speiser, Ph.D. Department of Radiation Oncology UT Southwestern Medical Center Dallas, TX September 1 st, 2012 CMPWG Workshop Medical
More informationHidenobu Tachibana The Cancer Institute Hospital of JFCR, Radiology Dept. The Cancer Institute of JFCR, Physics Dept.
2-D D Dose-CT Mapping in Geant4 Hidenobu Tachibana The Cancer Institute Hospital of JFCR, Radiology Dept. The Cancer Institute of JFCR, Physics Dept. Table of Contents Background & Purpose Materials Methods
More informationDosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System
Dosimetry Simulations with the UF-B Series Phantoms using the PENTRAN-MP Code System A. Al-Basheer, M. Ghita, G. Sjoden, W. Bolch, C. Lee, and the ALRADS Group Computational Medical Physics Team Nuclear
More information84 Luo heng-ming et al Vol. 12 Under the condition of one-dimensional slab geometry the transport equation for photons reads ±()N(x; μ; ) du d N(
Vol 12 No 7, July 23 cfl 23 Chin. Phys. Soc. 19-1963/23/12(7)/83-6 Chinese Physics and IOP Publishing Ltd The penetration, diffusion and energy deposition of high-energy photon * Luo heng-ming(λπξ) a),
More informationComplex Systems Simulations on the GPU
Complex Systems Simulations on the GPU Dr Paul Richmond Talk delivered by Peter Heywood University of Sheffield EMIT2015 Overview Complex Systems A Framework for Modelling Agents Benchmarking and Application
More informationOp#miza#on of an On- Board Imaging System for Rapid Radiotherapy. Erica Kemmerling, Meng Wu, He Yang, and Rebecca Fahrig
Op#miza#on of an On- Board Imaging System for Rapid Radiotherapy Erica Kemmerling, Meng Wu, He Yang, and Rebecca Fahrig Problem: mo#on during treatment Two overlapping thorax views during normal breathing
More informationGeant4 v9.5. Kernel III. Makoto Asai (SLAC) Geant4 Tutorial Course
Geant4 v9.5 Kernel III Makoto Asai (SLAC) Geant4 Tutorial Course Contents Fast simulation (Shower parameterization) Multi-threading Computing performance Kernel III - M.Asai (SLAC) 2 Fast simulation (shower
More informationAnalysis of Radiation Transport through Multileaf Collimators Using BEAMnrc Code
American Journal of Biomedical Engineering 216, 6(4): 124-131 DOI: 1.5923/j.ajbe.21664.3 Analysis of Radiation Transport through Multileaf Collimators Using BEAMnrc Code Ankit Kajaria 1,*, Neeraj Sharma
More informationGPU Based Convolution/Superposition Dose Calculation
GPU Based Convolution/Superposition Dose Calculation Todd McNutt 1, Robert Jacques 1,2, Stefan Pencea 2, Sharon Dye 2, Michel Moreau 2 1) Johns Hopkins University 2) Elekta Research Funded by and Licensed
More informationA dedicated tool for PET scanner simulations using FLUKA
A dedicated tool for PET scanner simulations using FLUKA P. G. Ortega FLUKA meeting June 2013 1 Need for in-vivo treatment monitoring Particles: The good thing is that they stop... Tumour Normal tissue/organ
More informationImplementation of the EGSnrc / BEAMnrc Monte Carlo code - Application to medical accelerator SATURNE43
International Journal of Innovation and Applied Studies ISSN 2028-9324 Vol. 6 No. 3 July 2014, pp. 635-641 2014 Innovative Space of Scientific Research Journals http://www.ijias.issr-journals.org/ Implementation
More informationGPU-based fast Monte Carlo simulation for radiotherapy dose calculation
GPU-based fast Monte Carlo simulation for radiotherapy dose calculation Xun Jia, Xuejun Gu, Yan Jiang Graves, Michael Folkerts, and Steve B. Jiang Center for Advanced Radiotherapy Technologies and Department
More informationIndrin Chetty Henry Ford Hospital Detroit, MI. AAPM Annual Meeting Houston 7:30-8:25 Mon 08/07/28 1/30
Review of TG105: Issues associated with clinical implementation of Monte Carlo-based photon and electron external beam treatment planning D. W. O. Rogers, Carleton Laboratory for Radiotherapy Physics,
More informationThe OpenGATE Collaboration
The OpenGATE Collaboration GATE developments and future releases GATE Users meeting, IEEE MIC 2015, San Diego Sébastien JAN CEA - IMIV sebastien.jan@cea.fr Outline Theranostics modeling Radiobiology Optics
More informationGeant4 Package in SPENVIS
N. Messios Belgian Institute for Space Aeronomy (BIRA-IASB) Outline Introduction Overview of the Geant4 models in SPENVIS Material definition GDML geometry definition Particle source GDML source geometry
More informationDevelopment of a Radiation Shielding Monte Carlo Code: RShieldMC
Development of a Radiation Shielding Monte Carlo Code: RShieldMC Shenshen GAO 1,2, Zhen WU 1,3, Xin WANG 1,2, Rui QIU 1,2, Chunyan LI 1,3, Wei LU 1,2, Junli LI 1,2*, 1.Department of Physics Engineering,
More informationDanesh TaNi & Amit Amritkar
GenIDLEST Co- Design Danesh TaNi & Amit Amritkar Collaborators Wu- chun Feng, Paul Sathre, Kaixi Hou, Sriram Chivukula, Hao Wang, Eric de Sturler, Kasia Swirydowicz Virginia Tech AFOSR- BRI Workshop Feb
More informationOutline. Monte Carlo Radiation Transport Modeling Overview (MCNP5/6) Monte Carlo technique: Example. Monte Carlo technique: Introduction
Monte Carlo Radiation Transport Modeling Overview () Lecture 7 Special Topics: Device Modeling Outline Principles of Monte Carlo modeling Radiation transport modeling with Utilizing Visual Editor (VisEd)
More informationLecture notes: Photon Monte Carlo simulation. Alex F Bielajew. Institute for National Measurement Standards. National Research Council of Canada
Lecture notes: Photon Monte Carlo simulation Alex F Bielajew Institute for National Measurement Standards National Research Council of Canada Ottawa, Canada K1A 0R6 Tel: 613-993-2715 FAX: 613-952-9865
More informationArtifact Mitigation in High Energy CT via Monte Carlo Simulation
PIERS ONLINE, VOL. 7, NO. 8, 11 791 Artifact Mitigation in High Energy CT via Monte Carlo Simulation Xuemin Jin and Robert Y. Levine Spectral Sciences, Inc., USA Abstract The high energy (< 15 MeV) incident
More informationInvestigation of tilted dose kernels for portal dose prediction in a-si electronic portal imagers
Investigation of tilted dose kernels for portal dose prediction in a-si electronic portal imagers Krista Chytyk MSc student Supervisor: Dr. Boyd McCurdy Introduction The objective of cancer radiotherapy
More informationMedical Physics Research Center, Mashhad University of Medical Sciences, Mashhad, Iran.
DXRaySMCS First User Friendly Interface Developed for Prediction of Diagnostic Radiology X-Ray Spectra Produced by Monte Carlo (MCNP-4C) Simulation in Iran M.T. Bahreyni Toosi a*, H. Moradi b, H. Zare
More informationMassively Parallel GPU-friendly Algorithms for PET. Szirmay-Kalos László, Budapest, University of Technology and Economics
Massively Parallel GPU-friendly Algorithms for PET Szirmay-Kalos László, http://cg.iit.bme.hu, Budapest, University of Technology and Economics (GP)GPU: CUDA (OpenCL) Multiprocessor N Multiprocessor 2
More informationTHE SIMULATION OF THE 4 MV VARIAN LINAC WITH EXPERIMENTAL VALIDATION
2007 International Nuclear Atlantic Conference - INAC 2007 Santos, SP, Brazil, September 30 to October 5, 2007 ASSOCIAÇÃO BRASILEIRA DE ENERGIA NUCLEAR - ABEN ISBN: 978-85-99141-02-1 THE SIMULATION OF
More informationDetermination of primary electron beam parameters in a Siemens Primus Linac using Monte Carlo simulation
Determination of primary electron beam parameters in a Siemens Primus Linac using Monte Carlo simulation Danial Seifi Makrani 1, Hadi Hasanzadeh 2*, Tayyeb Allahverdi Pourfallah 3, Arash Ghasemi 4, Majid
More informationArezoo Modiri Department of Radiation Oncology University of Maryland, Baltimore
Photon Optimization with GPU and Multi- Core CPU; What are the issues?, PhD Parallelization CPUs/Clusters/Cloud/GPUs Data management Outline Computation-Intensive Applications in Photon Radiotherapy Dose
More informationCut per region. Marc Verderi GEANT4 collaboration meeting 01/10/2002
Cut per region Marc Verderi GEANT4 collaboration meeting 01/10/2002 Introduction Cut here = «production threshold»; Not tracking cut; GEANT4 originally designed to allow a unique cut in range; Unique cut
More informationAbsorbed Fractions for Multi-Region Models of the Kidneys in ICRP/ICRU Voxel Phantoms
Absorbed Fractions for Multi-Region Models of the Kidneys in ICRP/ICRU Voxel Phantoms Sakae KINASE *, Masanori KIMURA, Shogo TAKAHARA and Toshimitsu HOMMA * 2-4 Shirane, Shirakata, Tokai-mura, Naka-gun,
More informationA Geant4 based treatment plan verification tool: commissioning of the LINAC model and DICOM-RT interface
A Geant4 based treatment plan verification tool: commissioning of the LINAC model and DICOM-RT interface Iwan Cornelius, Chris Poole. Stephen Warwick, Christian Langton Discipline of Physics, Faculty of
More informationValidation of GEANT4 Monte Carlo Simulation Code for 6 MV Varian Linac Photon Beam
Validation of GEANT4 Monte Carlo Code for 6 MV Varian Linac Photon Beam E. Salama ab*, A.S. Ali c, N. Emad d and A. Radi a a Physics Department, Faculty of Science, Ain Shams University, Cairo, Egypt;
More informationOBJECT ORIENTED DESIGN OF ANTHROPOMORPHIC PHANTOMS AND GEANT4-BASED IMPLEMENTATIONS ABSTRACT
International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 2009) Saratoga Springs, New York, May 3-7, 2009, on CD-ROM, American Nuclear Socety, LaGrange Park, IL (2009) OBJECT
More informationMCNP4C3-BASED SIMULATION OF A MEDICAL LINEAR ACCELERATOR
Computational Medical Physics Working Group Workshop II, Sep 3 Oct 3, 7 University of Florida (UF), Gainesville, Florida USA on CD-ROM, American Nuclear Society, LaGrange Park, IL (7) MCNP4C3-BASED SIMULATION
More informationATLAS Tracking Detector Upgrade studies using the Fast Simulation Engine
Journal of Physics: Conference Series PAPER OPEN ACCESS ATLAS Tracking Detector Upgrade studies using the Fast Simulation Engine To cite this article: Noemi Calace et al 2015 J. Phys.: Conf. Ser. 664 072005
More informationA Geometrical Modeller for HEP
A Geometrical Modeller for HEP R. Brun, A. Gheata CERN, CH 1211, Geneva 23, Switzerland M. Gheata ISS, RO 76900, Bucharest MG23, Romania For ALICE off-line collaboration Geometrical modelling generally
More informationAsynchronous and Fault-Tolerant Recursive Datalog Evalua9on in Shared-Nothing Engines
Asynchronous and Fault-Tolerant Recursive Datalog Evalua9on in Shared-Nothing Engines Jingjing Wang, Magdalena Balazinska, Daniel Halperin University of Washington Modern Analy>cs Requires Itera>on Graph
More informationAn approach to calculate and visualize intraoperative scattered radiation exposure
Peter L. Reicertz Institut für Medizinische Informatik An approach to calculate and visualize intraoperative scattered radiation exposure Markus Wagner University of Braunschweig Institute of Technology
More informationTopics for the TKR Software Review Tracy Usher, Leon Rochester
Topics for the TKR Software Review Tracy Usher, Leon Rochester Progress in reconstruction Reconstruction short-term plans Simulation Calibration issues Balloon-specific support Personnel and Schedule TKR
More informationPhysis: An Implicitly Parallel Framework for Stencil Computa;ons
Physis: An Implicitly Parallel Framework for Stencil Computa;ons Naoya Maruyama RIKEN AICS (Formerly at Tokyo Tech) GTC12, May 2012 1 è Good performance with low programmer produc;vity Mul;- GPU Applica;on
More informationOutline. In Situ Data Triage and Visualiza8on
In Situ Data Triage and Visualiza8on Kwan- Liu Ma University of California at Davis Outline In situ data triage and visualiza8on: Issues and strategies Case study: An earthquake simula8on Case study: A
More informationA new method of The Monte Carlo Simulation based on Hit stream for the LHAASO
A of The Monte Carlo Simulation based on Hit stream for the LHAASO Institute of High Energy Physics, CAS, Beijing, China. E-mail: wuhr@ihep.ac.cn Mingjun Chen Institute of High Energy Physics, CAS, Beijing,
More informationGPU-based fast gamma index calcuation
1 GPU-based fast gamma index calcuation 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 Xuejun Gu, Xun Jia, and Steve B. Jiang Center for Advanced Radiotherapy Technologies
More informationSoft GPGPUs for Embedded FPGAS: An Architectural Evaluation
Soft GPGPUs for Embedded FPGAS: An Architectural Evaluation 2nd International Workshop on Overlay Architectures for FPGAs (OLAF) 2016 Kevin Andryc, Tedy Thomas and Russell Tessier University of Massachusetts
More informationREAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT
REAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT Anand P Santhanam Assistant Professor, Department of Radiation Oncology OUTLINE Adaptive radiotherapy for head and
More informationMonte Carlo evaluation of the dosimetric uncertainty in matched 6 MV Elekta and Varian linear accelerators
Monte Carlo evaluation of the dosimetric uncertainty in matched 6 MV Elekta and Varian linear accelerators Dr Jessica E Lye PhD, BSc Hons. A thesis submitted in partial fulfilment of the requirements for
More informationOp#miza#on of D3Q19 La2ce Boltzmann Kernels for Recent Mul#- and Many-cores Intel Based Systems
Op#miza#on of D3Q19 La2ce Boltzmann Kernels for Recent Mul#- and Many-cores Intel Based Systems Ivan Giro*o 1235, Sebas4ano Fabio Schifano 24 and Federico Toschi 5 1 International Centre of Theoretical
More informationUnstructured Finite Volume Code on a Cluster with Mul6ple GPUs per Node
Unstructured Finite Volume Code on a Cluster with Mul6ple GPUs per Node Keith Obenschain & Andrew Corrigan Laboratory for Computa;onal Physics and Fluid Dynamics Naval Research Laboratory Washington DC,
More informationComparison of measured and Monte Carlo calculated electron beam central axis depth dose in water
Original article UDC: 681.324:519.245 Archive of Oncology 2001;9(2):83-7. Comparison of measured and Monte Carlo calculated electron beam central axis depth dose in water Darko LALIÆ 1 Radovan D. ILIÆ
More informationCOMPARISON OF DOSE CALCULATION ALGORITHMS FOR LEKSELL GAMMA KNIFE PERFEXION USING MONTE CARLO VOXEL PHANTOMS
COMPARISON OF DOSE CALCULATION ALGORITHMS FOR LEKSELL GAMMA KNIFE PERFEXION USING MONTE CARLO VOXEL PHANTOMS Jan Pipek 1, Josef Novotný Jr. 1,2,3, Josef Novotný 1, Petra Kozubíková 1 1 Faculty of Nuclear
More informationStarchart*: GPU Program Power/Performance Op7miza7on Using Regression Trees
Starchart*: GPU Program Power/Performance Op7miza7on Using Regression Trees Wenhao Jia, Princeton University Kelly A. Shaw, University of Richmond Margaret Martonosi, Princeton University *Sta7s7cal Tuning
More informationEVALUATION OF SPEEDUP OF MONTE CARLO CALCULATIONS OF TWO SIMPLE REACTOR PHYSICS PROBLEMS CODED FOR THE GPU/CUDA ENVIRONMENT
International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011) Rio de Janeiro, RJ, Brazil, May 8-12, 2011, on CD-ROM, Latin American Section (LAS)
More informationParallel Graph Coloring For Many- core Architectures
Parallel Graph Coloring For Many- core Architectures Mehmet Deveci, Erik Boman, Siva Rajamanickam Sandia Na;onal Laboratories Sandia National Laboratories is a multi-program laboratory managed and operated
More informationSource Model Tuning for a 6 MV Photon Beam used in Radiotherapy
Journal of Physics: Conference Series Source Model Tuning for a 6 MV Photon Beam used in Radiotherapy To cite this article: Lukas A Hirschi et al 2007 J. Phys.: Conf. Ser. 74 021008 View the article online
More informationIntroduction to Geant4 Visualization Geant4 v10.0.p3
Geant4 v10.0.p3 HepRep/HepRApp Makoto Asai (SLAC) Slides created by Joseph Perl (SLAC) So many options, it needs three title pages DAWN OpenGL 1 RayTracer OpenInventor DAWN HepRep/FRED 2 Qt gmocren DAWN
More informationSuitability Study of MCNP Monte Carlo Program for Use in Medical Physics
Nuclear Energy in Central Europe '98 Terme Catez, September 7 to 10, 1998 SI0100092 Suitability Study of MCNP Monte Carlo Program for Use in Medical Physics R. Jeraj Reactor Physics Division, Jozef Stefan
More informationClick to edit Master title style
Introduction to Serpent Code Fusion neutronics workshop, Cambridge, UK, June 11-12, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline style Serpent overview
More informationThe Cut and Thrust of CUDA
The Cut and Thrust of CUDA Luke Hodkinson Center for Astrophysics and Supercomputing Swinburne University of Technology Melbourne, Hawthorn 32000, Australia May 16, 2013 Luke Hodkinson The Cut and Thrust
More informationAssesing multileaf collimator effect on the build-up region using Monte Carlo method
Pol J Med Phys Eng. 2008;14(3):163-182. PL ISSN 1425-4689 doi: 10.2478/v10013-008-0014-0 website: http://www.pjmpe.waw.pl M. Zarza Moreno 1, 2, N. Teixeira 3, 4, A. P. Jesus 1, 2, G. Mora 1 Assesing multileaf
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Dose calculation algorithms! Correction based! Model based 1 Dose calculation algorithms!
More informationA Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,
More informationGeant4 Computing Performance Benchmarking and Monitoring
Journal of Physics: Conference Series PAPER OPEN ACCESS Geant4 Computing Performance Benchmarking and Monitoring To cite this article: Andrea Dotti et al 2015 J. Phys.: Conf. Ser. 664 062021 View the article
More informationStudy of scattered photons from the collimator system of Leksell Gamma Knife using the EGS4 Monte Carlo Code
Study of scattered photons from the collimator system of Leksell Gamma Knife using the EGS4 Monte Carlo Code Joel Y. C. Cheung Gamma Knife Centre, Canossa Hospital, 1 Old Peak Road, Hong Kong K. N. Yu
More informationUse of Monte Carlo modelling in radiotherapy linac design. David Roberts, PhD Senior Physicist Elekta
Use of Monte Carlo modelling in radiotherapy linac design David Roberts, PhD Senior Physicist Elekta Contents Overview of Elekta What we do Where we use Monte Carlo Codes and resources Example : Agility
More informationParallel Programming Pa,erns
Parallel Programming Pa,erns Bryan Mills, PhD Spring 2017 What is a programming pa,erns? Repeatable solu@on to commonly occurring problem It isn t a solu@on that you can t simply apply, the engineer has
More informationFusion PIC Code Performance Analysis on the Cori KNL System. T. Koskela*, J. Deslippe*,! K. Raman**, B. Friesen*! *NERSC! ** Intel!
Fusion PIC Code Performance Analysis on the Cori KNL System T. Koskela*, J. Deslippe*,! K. Raman**, B. Friesen*! *NERSC! ** Intel! tkoskela@lbl.gov May 18, 2017-1- Outline Introduc3on to magne3c fusion
More informationGPGPU in Film Production. Laurence Emms Pixar Animation Studios
GPGPU in Film Production Laurence Emms Pixar Animation Studios Outline GPU computing at Pixar Demo overview Simulation on the GPU Future work GPU Computing at Pixar GPUs have been used for real-time preview
More informationCRAY User Group Mee'ng May 2010
Applica'on Accelera'on on Current and Future Cray Pla4orms Alice Koniges, NERSC, Berkeley Lab David Eder, Lawrence Livermore Na'onal Laboratory (speakers) Robert Preissl, Jihan Kim (NERSC LBL), Aaron Fisher,
More information