| Peer-Reviewed

Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code

Received: 23 October 2016     Accepted: 3 November 2016     Published: 23 November 2016
Views:       Downloads:
Abstract

This study aimed to develop a new phase-space tool for Geant4 code using the HDF5 C++ scientific data library. The tool can be easy incorporated into existing Geant4 applications and provides particle recycling and rotational splitting capabilities which can be useful for faster modeling symmetric systems such as medical linear accelerator. The validation of this phase-space I/O routines has been performed in a very basic geometry. Thus, taking into account a homogeneous water phantom, the depth dose curve of a 20 MeV electron beam hitting a small target made of tungsten has been calculated and compared to ones produced in the simulation without phase-space technique and simulation with IAEA phase-space I/O routines. This study shown an excellent agreement found between different calculated depth dose curves, allowing us to validate our new phase-space I/O routines. Moreover, the CPU time spent by simulation without variance reduction technique can be reduced by 27% when this method was applied which is the same factor obtained with IAEA phase-space I/O routines. The phase-space tool describing in this paper, have been implemented for Geant4 code by using HDF5 C++ data managing library, the associated classes are few and easy to incorporate into an existing Geant4 user code, and it is freely available on GitHub ( https://github.com/EL-Bakkali-Jaafar/G4PhpH5/).

Published in American Journal of Mathematical and Computer Modelling (Volume 1, Issue 1)
DOI 10.11648/j.ajmcm.20160101.14
Page(s) 25-31
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2016. Published by Science Publishing Group

Keywords

Geant4, Phase-Space, Hdf5, C++, Variance Reduction Technique, Monte Carlo

References
[1] S. Agostinelli, J. Allison, K. Amako, J. Apostolakis, H. Araujo, P. Arce et al., “Geant4 - A Simulation Toolkit,” Nuclear Instruments and Methods A 506, 2003, pp. 250-303.
[2] Allison J et al., “Geant4 Developments and Applications,” IEEE Transactions on Nuclear Science 53, 2006, pp. 270–278.
[3] The HDF5 Group, HDF5 Dara Format. [Online] Available: http://www.hdfgroup.org/HDF5/release/obtainsrc.html.
[4] W. Wenjie, “Comparison of Portable Binary Data Formats within a Cosmological Simulation,” M. Sc. in High Performance Computing, The University of Edinburgh, 2005.
[5] B. Caccia, C. Andenna, and G. A. P Cirrone, “MedLinac2: a GEANT4 based software package for radiotherapy,” Annali dell’Istituto superiore di sanita 46, 2010, pp. 173–177.
[6] M. A Cortés-Giraldo, J. M Quesada Molina, M. I Gallardo, R. Capote, “Geant4 Interface to Work with IAEA Phase-Space Files”, 2009.
[7] M. A Cortés-Giraldo, J. M Quesada Molina, M. I Gallardo, R. Capote, “An implementation to read and write IAEA phase-space files in GEANT4-based simulations,” Int. J. Radiat. Biol. 88, 2012, pp.200-208.
[8] R. Capote R and I. Kawrakow, “Read/write routines implementing the IAEA phsp format, version of December 2009”. [Online] Available: http://www-nds.iaea.org/phsp/software/iaea phsp Dec2009.zip#phsp rw
[9] J. EL Bakkali, T. EL Bardouni, S. Safavi, M. Mohammed, S. Mroan, “Behaviors of the percentage depth dose curves along the beam axis of a phantom filled with different clinical PTO Objects, a Monte Carlo Geant4 study,” Radiation Physics and Chemistry 125, 2016, pp.199-204.
[10] J. EL Bakkali, T. EL Bardouni, “Validation of Monte Carlo Geant4 code for a 6 MV Varian linac,” Journal of King Saud University–Science, 2016, in press.
Cite This Article
  • APA Style

    Jaafar EL Bakkali, Abderrahim Doudouh, Khalid Bouyakhlef, Laila Baddouh, Keltoum Dahmani, et al. (2016). Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code. American Journal of Mathematical and Computer Modelling, 1(1), 25-31. https://doi.org/10.11648/j.ajmcm.20160101.14

    Copy | Download

    ACS Style

    Jaafar EL Bakkali; Abderrahim Doudouh; Khalid Bouyakhlef; Laila Baddouh; Keltoum Dahmani, et al. Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code. Am. J. Math. Comput. Model. 2016, 1(1), 25-31. doi: 10.11648/j.ajmcm.20160101.14

    Copy | Download

    AMA Style

    Jaafar EL Bakkali, Abderrahim Doudouh, Khalid Bouyakhlef, Laila Baddouh, Keltoum Dahmani, et al. Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code. Am J Math Comput Model. 2016;1(1):25-31. doi: 10.11648/j.ajmcm.20160101.14

    Copy | Download

  • @article{10.11648/j.ajmcm.20160101.14,
      author = {Jaafar EL Bakkali and Abderrahim Doudouh and Khalid Bouyakhlef and Laila Baddouh and Keltoum Dahmani and Hamid Mansouri},
      title = {Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code},
      journal = {American Journal of Mathematical and Computer Modelling},
      volume = {1},
      number = {1},
      pages = {25-31},
      doi = {10.11648/j.ajmcm.20160101.14},
      url = {https://doi.org/10.11648/j.ajmcm.20160101.14},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.ajmcm.20160101.14},
      abstract = {This study aimed to develop a new phase-space tool for Geant4 code using the HDF5 C++ scientific data library. The tool can be easy incorporated into existing Geant4 applications and provides particle recycling and rotational splitting capabilities which can be useful for faster modeling symmetric systems such as medical linear accelerator. The validation of this phase-space I/O routines has been performed in a very basic geometry. Thus, taking into account a homogeneous water phantom, the depth dose curve of a 20 MeV electron beam hitting a small target made of tungsten has been calculated and compared to ones produced in the simulation without phase-space technique and simulation with IAEA phase-space I/O routines. This study shown an excellent agreement found between different calculated depth dose curves, allowing us to validate our new phase-space I/O routines. Moreover, the CPU time spent by simulation without variance reduction technique can be reduced by 27% when this method was applied which is the same factor obtained with IAEA phase-space I/O routines. The phase-space tool describing in this paper, have been implemented for Geant4 code by using HDF5 C++ data managing library, the associated classes are few and easy to incorporate into an existing Geant4 user code, and it is freely available on GitHub ( https://github.com/EL-Bakkali-Jaafar/G4PhpH5/).},
     year = {2016}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Investigation of the HDF5 C++ Library in Development of New Phase-Space I/O for Radiotherapy Simulation Using Monte Carlo Geant4 Code
    AU  - Jaafar EL Bakkali
    AU  - Abderrahim Doudouh
    AU  - Khalid Bouyakhlef
    AU  - Laila Baddouh
    AU  - Keltoum Dahmani
    AU  - Hamid Mansouri
    Y1  - 2016/11/23
    PY  - 2016
    N1  - https://doi.org/10.11648/j.ajmcm.20160101.14
    DO  - 10.11648/j.ajmcm.20160101.14
    T2  - American Journal of Mathematical and Computer Modelling
    JF  - American Journal of Mathematical and Computer Modelling
    JO  - American Journal of Mathematical and Computer Modelling
    SP  - 25
    EP  - 31
    PB  - Science Publishing Group
    SN  - 2578-8280
    UR  - https://doi.org/10.11648/j.ajmcm.20160101.14
    AB  - This study aimed to develop a new phase-space tool for Geant4 code using the HDF5 C++ scientific data library. The tool can be easy incorporated into existing Geant4 applications and provides particle recycling and rotational splitting capabilities which can be useful for faster modeling symmetric systems such as medical linear accelerator. The validation of this phase-space I/O routines has been performed in a very basic geometry. Thus, taking into account a homogeneous water phantom, the depth dose curve of a 20 MeV electron beam hitting a small target made of tungsten has been calculated and compared to ones produced in the simulation without phase-space technique and simulation with IAEA phase-space I/O routines. This study shown an excellent agreement found between different calculated depth dose curves, allowing us to validate our new phase-space I/O routines. Moreover, the CPU time spent by simulation without variance reduction technique can be reduced by 27% when this method was applied which is the same factor obtained with IAEA phase-space I/O routines. The phase-space tool describing in this paper, have been implemented for Geant4 code by using HDF5 C++ data managing library, the associated classes are few and easy to incorporate into an existing Geant4 user code, and it is freely available on GitHub ( https://github.com/EL-Bakkali-Jaafar/G4PhpH5/).
    VL  - 1
    IS  - 1
    ER  - 

    Copy | Download

Author Information
  • Nuclear Medicine Department, Military Hospital Mohammed V, Rabat, Morocco

  • Nuclear Medicine Department, Military Hospital Mohammed V, Rabat, Morocco

  • University Mohammed V, Souissi, Faculty of Medicine and Pharmacy, Rabat, Morocco

  • Radiotherapy Department, Military Hospital Mohammed V, Rabat, Morocco

  • Radiotherapy Department, Military Hospital Mohammed V, Rabat, Morocco

  • Radiotherapy Department, Military Hospital Mohammed V, Rabat, Morocco

  • Sections