MPEXS benchmark results - phase space data - Akinori Kimura 14 - PowerPoint PPT Presentation

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 are evaluated for the validation. • To compare the calculation speeds of the simulations between MPEXS and Geant4 with phase space data 2 14 February, 2017

Contents of Phase Space Data header • the total number of particles, the total number of photons, the max. kinetic • energy of the particles, the min. kinetic energy of the particles, etc. particle information of the total number of particles • LATCH : the particle’s history, the particle charge • E : the particle total energy • X : the particle X-position [cm], Y : the particle Y -position [cm] on a scoring plane • U : the X-direction cosine, V : the Y -direction cosine at the particle position • WT : the particle’s weight: the sign of the Z-direction cosine, 
 • W : = SQRT(1-(U^2+V^2)) 3 14 February, 2017

Dataset of Phase Space Data • SIEMENS Primus electrons or X-rays (a) total energy, (b) total number of source particles, (c) number of files, 
 • and (d) total file sizes • electron source 
 (a) (b) (c) (d) (a) (b) (c) (d) 6 MeV 15 MeV 43,976,682 8 1.1 GB 91,321,327 8 2.4 GB 9 MeV 18 MeV 61,634,988 8 1.6 GB 114,334,660 7 3.0 GB 12 MeV 21 MeV 89,084,463 8 2.3 GB 167,747,990 8 4.4 GB • X-rays source (a) (b) (c) (d) (a) (b) (c) (d) 6 MeV 83,171,454 8 2.2 GB 18 MeV 60,284,412 8 1.6 GB 4 14 February, 2017

Geometry Configuration of Simulation water phantom SSD: 100 cm size: 20 cm x 20 cm x 20 cm (source surface distance) d # of voxels: 100 x 100 x 100 phase space plane size: 10 cm x 10 cm distance to surface ( d ): electron phase space source: d =1 cm X-ray phase space source: d = 50 cm 5 14 February, 2017

Geant4 simulation • version: 10.2 p2 • Physics processes: • for electrons and positrons • G4eMultipleScattering, G4eIonisation, G4eBremsstrahlung, G4eplusAnnihilation • for X-rays • G4PhotoElectricEffect, G4ComptonScattering, G4GammaConversion • range cut value: 1 mm (default) 6 14 February, 2017

MPEXS • Physics configuration is the same as one of Geant4: • Physics processes: • for electrons and positrons • multiple scattering, ionization, Bremsstrahlung, positron annihilation • for X-rays • photoelectric effect, Compton scattering, gamma conversion • range cut value: 1 mm (default) 7 14 February, 2017

Depth Dose at Center of Phantom (Primus electrons) phantom MPEXS (circles) vs Geant4 (lines) • central axis voxel of 4 mm x 4 mm x 2 mm thick • electrons: 6, 9, 12, 18 and 21 MeV • 4 mm depth dose at center of phantom (4mm x4mm) 1 .4E-11 MPEXS 6MeV MPEXS 9MeV 1 .2E-11 MPEXS 12MeV 4 mm MPEXS 15MeV dose per source particle [Gy] 1 .0E-11 MPEXS 18MeV MPEXS 21MeV G4 6MeV 8 .0E-12 G4 9MeV G4 12MeV G4 15MeV 6 .0E-12 G4 18MeV G4 21MeV 4 .0E-12 2 .0E-12 0 .0E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 depth [cm] preliminary result 8 14 February, 2017

dose per source particle [Gy] Profiles of Primus electrons at 1.5 cm depth X [cm] X-axis Profiles (Primus electrons at 1.5 cm depth) phantom MPEXS (circles) vs Geant4 (lines) • voxel of 8 mm out of plane and 2 mm x 2mm thick • 4 mm electrons: 6, 9, 12, 18 and 21 MeV • X-axis 2 mm depth dose at center of phantom (4mm x4mm) 1 .40E-11 1 .4E-11 6 MeV 9 MeV MPEXS 6MeV 1 .20E-11 12 MeV MPEXS 9MeV 1 .2E-11 MPEXS 12MeV 15 MeV MPEXS 15MeV 18MeV 1 .00E-11 dose per source particle [Gy] 1 .0E-11 MPEXS 18MeV 21 MeV MPEXS 21MeV EGSnrc 6 MeV G4 6MeV 8 .00E-12 EGSnrc 9 MeV 8 .0E-12 G4 9MeV EGSnrc 12 MeV G4 12MeV EGSnrc 15 MeV 6 .00E-12 G4 15MeV 6 .0E-12 EGSnrc 18MeV G4 18MeV EGSnrc 21 MeV G4 21MeV 4 .00E-12 4 .0E-12 2 .00E-12 2 .0E-12 0 .0E+00 0 .00E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 depth [cm] 1.5 cm preliminary result 9 14 February, 2017

Profiles of Primus electrons at 1.5 cm depth dose per source particle [Gy] Y [cm] Y-axis Profiles (Primus electrons at 1.5 cm depth) phantom MPEXS (circles) vs Geant4 (lines) • voxel of 8 mm out of plane and 2 mm x 2mm thick • 4 mm electrons: 6, 9, 12, 18 and 21 MeV • Y -axis 2 mm depth dose at center of phantom (4mm x4mm) 1 .40E-11 6 MeV 1 .4E-11 9 MeV MPEXS 6MeV 12 MeV 1 .20E-11 MPEXS 9MeV 1 .2E-11 15 MeV MPEXS 12MeV 18MeV MPEXS 15MeV 1 .00E-11 21 MeV dose per source particle [Gy] 1 .0E-11 MPEXS 18MeV EGSnrc 6 MeV MPEXS 21MeV EGSnrc 9 MeV G4 6MeV 8 .00E-12 8 .0E-12 EGSnrc 12 MeV G4 9MeV EGSnrc 15 MeV G4 12MeV 6 .00E-12 EGSnrc 18MeV G4 15MeV 6 .0E-12 EGSnrc 21 MeV G4 18MeV G4 21MeV 4 .00E-12 4 .0E-12 2 .00E-12 2 .0E-12 0 .0E+00 0 .00E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 depth [cm] 1.5 cm preliminary result 10 14 February, 2017

Depth Dose at Center of Phantom (Primus X-rays) phantom MPEXS (circles) vs Geant4 (lines) • central axis voxel of 4 mm x 4 mm x 2 mm thick • X-rays: 6 and 18 MeV • 4 mm 4 .5E-13 MPEXS 6MV 4 .0E-13 MPEXS 18MV G4 6MV 3 .5E-13 4 mm G4 18MV dose per source particle [Gy] 3 .0E-13 2 .5E-13 2 .0E-13 1 .5E-13 1 .0E-13 5 .0E-14 0 .0E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 depth [cm] preliminary result 11 14 February, 2017

X-axis profiles (Primus X-rays) 4 .5E-13 MPEXS 6MV 4 .0E-13 MPEXS 18MV MPEXS (circles) vs Geant4 (lines) • G4 6MV 3 .5E-13 G4 18MV dose per source particle [Gy] 3 .0E-13 voxel of 8 mm out of plane and 2 mm x 2mm thick 
 • 2 .5E-13 at 3.3 cm, 10.1 cm, or 18.1 cm depths 2 .0E-13 X-rays: 6 and 18 MeV • 1 .5E-13 1 .0E-13 5 .0E-14 0 .0E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 depth [cm] 3.3 cm 10.1 cm 18.1 cm X-axis profiles (siemens primus 6 MV X-ray) X-axis profiles (siemens primus 18 MV X-ray) 4 .5E-13 depth: 3.2 to 3.4 cm depth: 3.2 to 3.4 cm 2 .5E-13 depth: 10 to 10.2 cm depth: 10 to 10.2 cm 4 .0E-13 depth: 18 to 18.2 cm depth: 18 to 18.2 cm dose per source particle [cGy] 3 .5E-13 2 .0E-13 EGSnrc 3.2 to 3.4 cm dose per source particle [cGy] EGSnrc: 3.2 to 3.4 cm 3 .0E-13 EGSnrc 10 to 10.2 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm EGSnrc 18 to 18.2 cm 1 .5E-13 2 .5E-13 2 .0E-13 1 .0E-13 1 .5E-13 1 .0E-13 5 .0E-14 5 .0E-14 0 .0E+00 0 .0E+00 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 preliminary result 12 14 February, 2017 X axis [cm] X axis [cm]

Y-axis profiles (Primus X-rays) 4 .5E-13 MPEXS 6MV 4 .0E-13 MPEXS 18MV MPEXS (circles) vs Geant4 (lines) • G4 6MV 3 .5E-13 G4 18MV dose per source particle [Gy] 3 .0E-13 voxel of 8 mm out of plane and 2 mm x 2mm thick 
 • 2 .5E-13 at 3.3 cm, 10.1 cm, or 18.1 cm depths 2 .0E-13 X-rays: 6 and 18 MeV • 1 .5E-13 1 .0E-13 5 .0E-14 0 .0E+00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 depth [cm] 3.3 cm 10.1 cm 18.1 cm Y-axis profiles (siemens primus 6 MV X-ray) Y-axis profiles (siemens primus 18 MV X-ray) 4 .5E-13 depth: 3.2 to 3.4 cm depth: 3.2 to 3.4 cm 2 .5E-13 depth: 10 to 10.2 cm depth: 10 to 10.2 cm 4 .0E-13 depth: 18 to 18.2 cm depth: 18 to 18.2 cm dose per source particle [cGy] 3 .5E-13 2 .0E-13 EGSnrc 3.2 to 3.4 cm dose per source particle [cGy] EGSnrc: 3.2 to 3.4 cm 3 .0E-13 EGSnrc 10 to 10.2 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm EGSnrc 18 to 18.2 cm 1 .5E-13 2 .5E-13 2 .0E-13 1 .0E-13 1 .5E-13 1 .0E-13 5 .0E-14 5 .0E-14 0 .0E+00 0 .0E+00 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 10 preliminary result 13 14 February, 2017 Y axis [cm] Y axis [cm]

Calculation Speeds MPEXS • GPU: NVIDIA GTX1080 • no. of processes: # threads: 256, # blocks: 256 • OS: CentOS 7, CUDA: 8.0 • Geant4 10.2p2 • CPU: intel Core i7-6700 3.4GHz • no. of processes: single thread • OS: CentOS 7 • Phase space data 1: • user time system time total time source particles: Primus x-rays • x 49.7 MPEXS 17.98 sec. 7.64 sec. 25.62 sec. total energy: 6MV • Geant4 1272.17 sec. 0.12 sec. 1274.30 sec. speed up # of source particles: 10,394,741 • Phase space data 2: • user time system time total time source particles: Primus electron • x 47.4 total energy: 6MeV • MPEXS 10.75 sec. 3.76 sec. 14.61 sec. # of source particles: 5,498,221 • Geant4 692.63 sec. 0.07 sec. 692.63 sec. speed up 14 14 February, 2017

Plans • to measure the calculation speed of MPEXS with phase space data on the in2p3 machine, cca.in2p3.fr • I have submitted a job as following command, but the job has been kept in the queue and has never run. • $ qsub -l os=cl7,GPU=1 -q mc_gpu_medium -pe multicores 1 exe • 1 GPU (Tesla K80), 1 CPU core • batch queue: mc_gpu_medium • exe: the script to run MPEXS • to optimize the code of MPEXS with phase space files to speed up 15 14 February, 2017

Summary • The preliminary results show the valid results between MPEXS and Geant4 with the phase space data of SIMENS Primus electrons or X-rays. • Calculation speeds of MPEXS with the phase space data is about 50 times faster than Geant4. • MPEXS: NVIDIA GTX1080 • Geant4: single thread of intel Core i7-6700 3.4GHz 16 14 February, 2017

Another Work • MPEXSapp : a desktop application software of MPEXS • GUI: Qt 5.7.0 17 14 February, 2017