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

▶
mpexs benchmark results
SMART_READER_LITE
LIVE PREVIEW

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

Dec 03, 2023 104 likes •295 views

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

slide-1
SLIDE 1

MPEXS benchmark results

  • phase space data -

Akinori Kimura 14 February 2017

slide-2
SLIDE 2

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

slide-3
SLIDE 3

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

slide-4
SLIDE 4

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

  • X-rays source

4

(a) (b) (c) (d) (a) (b) (c) (d) 6 MeV 43,976,682 8 1.1 GB 15 MeV 91,321,327 8 2.4 GB 9 MeV 61,634,988 8 1.6 GB 18 MeV 114,334,660 7 3.0 GB 12 MeV 89,084,463 8 2.3 GB 21 MeV 167,747,990 8 4.4 GB (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

slide-5
SLIDE 5

14 February, 2017

Geometry Configuration of Simulation

5

size: 20 cm x 20 cm x 20 cm

water phantom

# 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

SSD: 100 cm

d

(source surface distance)

slide-6
SLIDE 6

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

slide-7
SLIDE 7

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

slide-8
SLIDE 8

14 February, 2017

Depth Dose at Center of Phantom (Primus electrons)

  • 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

8

0.0E+00 2.0E-12 4.0E-12 6.0E-12 8.0E-12 1.0E-11 1.2E-11 1.4E-11 1 2 3 4 5 6 7 8 9 10 11 12

dose per source particle [Gy] depth [cm] depth dose at center of phantom (4mm x4mm)

MPEXS 6MeV MPEXS 9MeV MPEXS 12MeV MPEXS 15MeV MPEXS 18MeV MPEXS 21MeV G4 6MeV G4 9MeV G4 12MeV G4 15MeV G4 18MeV G4 21MeV

4 mm 4 mm

phantom

preliminary result

slide-9
SLIDE 9

14 February, 2017

X-axis Profiles (Primus electrons at 1.5 cm depth)

X-axis

0.00E+00 2.00E-12 4.00E-12 6.00E-12 8.00E-12 1.00E-11 1.20E-11 1.40E-11

  • 10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10

dose per source particle [Gy]

X [cm]

Profiles of Primus electrons at 1.5 cm depth

6 MeV 9 MeV 12 MeV 15 MeV 18MeV 21 MeV EGSnrc 6 MeV EGSnrc 9 MeV EGSnrc 12 MeV EGSnrc 15 MeV EGSnrc 18MeV EGSnrc 21 MeV

  • MPEXS (circles) vs Geant4 (lines)
  • voxel of 8 mm out of plane and 2 mm x 2mm thick
  • electrons: 6, 9, 12, 18 and 21 MeV

9

2 mm 4 mm

phantom

0.0E+00 2.0E-12 4.0E-12 6.0E-12 8.0E-12 1.0E-11 1.2E-11 1.4E-11 1 2 3 4 5 6 7 8 9 10 11 12

dose per source particle [Gy] depth [cm] depth dose at center of phantom (4mm x4mm)

MPEXS 6MeV MPEXS 9MeV MPEXS 12MeV MPEXS 15MeV MPEXS 18MeV MPEXS 21MeV G4 6MeV G4 9MeV G4 12MeV G4 15MeV G4 18MeV G4 21MeV

1.5 cm

preliminary result

slide-10
SLIDE 10

14 February, 2017

Y-axis Profiles (Primus electrons at 1.5 cm depth)

Y

  • axis
  • MPEXS (circles) vs Geant4 (lines)
  • voxel of 8 mm out of plane and 2 mm x 2mm thick
  • electrons: 6, 9, 12, 18 and 21 MeV

10

2 mm 4 mm

phantom

0.00E+00 2.00E-12 4.00E-12 6.00E-12 8.00E-12 1.00E-11 1.20E-11 1.40E-11

  • 10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10

dose per source particle [Gy]

Y [cm]

Profiles of Primus electrons at 1.5 cm depth

6 MeV 9 MeV 12 MeV 15 MeV 18MeV 21 MeV EGSnrc 6 MeV EGSnrc 9 MeV EGSnrc 12 MeV EGSnrc 15 MeV EGSnrc 18MeV EGSnrc 21 MeV

0.0E+00 2.0E-12 4.0E-12 6.0E-12 8.0E-12 1.0E-11 1.2E-11 1.4E-11 1 2 3 4 5 6 7 8 9 10 11 12

dose per source particle [Gy] depth [cm] depth dose at center of phantom (4mm x4mm)

MPEXS 6MeV MPEXS 9MeV MPEXS 12MeV MPEXS 15MeV MPEXS 18MeV MPEXS 21MeV G4 6MeV G4 9MeV G4 12MeV G4 15MeV G4 18MeV G4 21MeV

1.5 cm

preliminary result

slide-11
SLIDE 11

14 February, 2017

Depth Dose at Center of Phantom (Primus X-rays)

  • MPEXS (circles) vs Geant4 (lines)
  • central axis voxel of 4 mm x 4 mm x 2 mm thick
  • X-rays: 6 and 18 MeV

11

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13 3.0E-13 3.5E-13 4.0E-13 4.5E-13 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

dose per source particle [Gy]

depth [cm] MPEXS 6MV MPEXS 18MV G4 6MV G4 18MV

4 mm 4 mm

phantom

preliminary result

slide-12
SLIDE 12

14 February, 2017

X-axis profiles (Primus X-rays)

  • MPEXS (circles) vs Geant4 (lines)
  • voxel of 8 mm out of plane and 2 mm x 2mm thick 


at 3.3 cm, 10.1 cm, or 18.1 cm depths

  • X-rays: 6 and 18 MeV

12

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13

  • 10
  • 9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10 dose per source particle [cGy]

X axis [cm]

X-axis profiles (siemens primus 6 MV X-ray)

depth: 3.2 to 3.4 cm depth: 10 to 10.2 cm depth: 18 to 18.2 cm EGSnrc: 3.2 to 3.4 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13 3.0E-13 3.5E-13 4.0E-13 4.5E-13

  • 10
  • 9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10 dose per source particle [cGy] X axis [cm]

X-axis profiles (siemens primus 18 MV X-ray)

depth: 3.2 to 3.4 cm depth: 10 to 10.2 cm depth: 18 to 18.2 cm EGSnrc 3.2 to 3.4 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13 3.0E-13 3.5E-13 4.0E-13 4.5E-13 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

dose per source particle [Gy]

depth [cm] MPEXS 6MV MPEXS 18MV G4 6MV G4 18MV

3.3 cm 10.1 cm 18.1 cm

preliminary result

slide-13
SLIDE 13

14 February, 2017

Y-axis profiles (Primus X-rays)

  • MPEXS (circles) vs Geant4 (lines)
  • voxel of 8 mm out of plane and 2 mm x 2mm thick 


at 3.3 cm, 10.1 cm, or 18.1 cm depths

  • X-rays: 6 and 18 MeV

13

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13 3.0E-13 3.5E-13 4.0E-13 4.5E-13 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

dose per source particle [Gy]

depth [cm] MPEXS 6MV MPEXS 18MV G4 6MV G4 18MV

3.3 cm 10.1 cm 18.1 cm

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13

  • 10
  • 9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10 dose per source particle [cGy]

Y axis [cm]

Y-axis profiles (siemens primus 6 MV X-ray)

depth: 3.2 to 3.4 cm depth: 10 to 10.2 cm depth: 18 to 18.2 cm EGSnrc: 3.2 to 3.4 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm

0.0E+00 5.0E-14 1.0E-13 1.5E-13 2.0E-13 2.5E-13 3.0E-13 3.5E-13 4.0E-13 4.5E-13

  • 10
  • 9 -8 -7 -6 -5 -4 -3 -2 -1 0

1 2 3 4 5 6 7 8 9 10 dose per source particle [cGy] Y axis [cm]

Y-axis profiles (siemens primus 18 MV X-ray)

depth: 3.2 to 3.4 cm depth: 10 to 10.2 cm depth: 18 to 18.2 cm EGSnrc 3.2 to 3.4 cm EGSnrc 10 to 10.2 cm EGSnrc 18 to 18.2 cm

preliminary result

slide-14
SLIDE 14

14 February, 2017

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:
  • source particles: Primus x-rays
  • total energy: 6MV
  • # of source particles: 10,394,741
  • Phase space data 2:
  • source particles: Primus electron
  • total energy: 6MeV
  • # of source particles: 5,498,221

14 user time system time total time MPEXS 17.98 sec. 7.64 sec. 25.62 sec. Geant4 1272.17 sec. 0.12 sec. 1274.30 sec.

x 47.4 speed up

user time system time total time MPEXS 10.75 sec. 3.76 sec. 14.61 sec. Geant4 692.63 sec. 0.07 sec. 692.63 sec.

x 49.7 speed up

slide-15
SLIDE 15

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

slide-16
SLIDE 16

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

slide-17
SLIDE 17

14 February, 2017

Another Work

  • MPEXSapp: a desktop application software of MPEXS
  • GUI: Qt 5.7.0

17