Background
Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro‐ and nanodosimetry, imaging, radiation protection, and nuclear ...medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing.
Aims
To respond to these needs, we developed G4‐Med, a benchmarking and regression testing system of Geant4 for medical physics.
Materials and Methods
G4‐Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4‐Med are executed on the CERN computing infrastructure via the use of the geant‐val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes.
Results
This paper describes the tests included in G4‐Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data.
Discussion
Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics.
Conclusion
The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.
A model for the simulation of orientational effects in straight and bent periodic atomic structures is presented. The continuum potential approximation has been adopted. The model allows the ...manipulation of particle trajectories by means of straight and bent crystals and the scaling of the cross sections of hadronic and electromagnetic processes for channeled particles. Based on such a model, an extension of the Geant4 toolkit has been developed. The code has been validated against data from channeling experiments carried out at CERN.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
•Nuclear reactions are fundamental for Ion-therapy simulations.•Theoretical models running time are too long for medical applications.•We are proposing the idea of using a Deep Learning algorithm to ...emulate the model.•We obtained promising preliminary results.
A reliable model to simulate nuclear interactions is fundamental for Ion-therapy. We already showed how BLOB (“Boltzmann-Langevin One Body”), a model developed to simulate heavy ion interactions up to few hundreds of MeV/u, could simulate also 12C reactions in the same energy domain. However, its computation time is too long for any medical application. For this reason we present the possibility of emulating it with a Deep Learning algorithm.
The BLOB final state is a Probability Density Function (PDF) of finding a nucleon in a position of the phase space. We discretised this PDF and trained a Variational Auto-Encoder (VAE) to reproduce such a discrete PDF. As a proof of concept, we developed and trained a VAE to emulate BLOB in simulating the interactions of 12C with 12C at 62 MeV/u. To have more control on the generation, we forced the VAE latent space to be organised with respect to the impact parameter (b) training a classifier of b jointly with the VAE.
The distributions obtained from the VAE are similar to the input ones and the computation time needed to use the VAE as a generator is negligible.
We show that it is possible to use a Deep Learning approach to emulate a model developed to simulate nuclear reactions in the energy range of interest for Ion-therapy. We foresee the implementation of the generation part in C++ and to interface it with the most used Monte Carlo toolkit: Geant4.
The Geant4 collaboration regularly performs validation and regression tests. The results are stored in a central repository and can be easily accessed via a web application. In this article we ...describe the Geant4 physics validation repository which consists of a relational database storing experimental data and Geant4 test results, a java API and a web application. The functionality of these components and the technology choices we made are also described.
Simulation of orientational coherent effects via Geant4 beam manipulation of high-and very-high-energy particle beams is a hot topic in accelerator physics. Coherent effects of ultra-relativistic ...particles in bent crystals allow the steering of particle trajectories thanks to the strong electrical field generated between atomic planes. Recently, a collimation experiment with bent crystals was carried out at the CERN-LHC, paving the way to the usage of such technology in current and future accelerators. Geant4 is a widely used object-oriented tool-kit for the Monte Carlo simulation of the interaction of particles with matter in high-energy physics. Moreover, its areas of application include also nuclear and accelerator physics, as well as studies in medical and space science. We present the first Geant4 extension for the simulation of orientational effects in straight and bent crystals for high energy charged particles. The model allows the manipulation of particle trajectories by means of straight and bent crystals and the scaling of the cross sections of hadronic and electromagnetic processes for channeled particles. Based on such a model, an extension of the Geant4 toolkit has been developed. The code and the model have been validated by comparison with published experimental data regarding the deflection efficiency via channeling and the variation of the rate of inelastic nuclear interactions.
•Geant4 is widely used for developing Monte Carlo simulations in medical physics.•It has severe limitations in reproducing nuclear interactions below 100 MeV/n.•We interfaced two models dedicated to ...low energy nuclear interactions with Geant4.•We obtained good agreement with experimental data of 12C fragmentation at 62 MeV/n.
Monte Carlo (MC) simulations are widely used for medical applications and nuclear reaction models are fundamental for the simulation of the particle interactions with patients in ion therapy. Therefore, it is of utmost importance to have reliable models in MC simulations for such interactions. Geant4 is one of the most used toolkits for MC simulation. However, its models showed severe limitations in reproducing the yields measured in the interaction of ion beams below 100 MeV/u with thin targets. For this reason, we interfaced two models, SMF (“Stochastic Mean Field”) and BLOB (“Boltzmann-Langevin One Body”), dedicated to simulate such reactions, with Geant4.
Both SMF and BLOB are semi-classical, one-body approaches to solve the Boltzmann-Langevin equation. They include an identical treatment of the mean-field propagation, on the basis of the same effective interaction, but they differ in the way fluctuations are included. Furthermore, we tested a correction to the excitation energy calculated for the light fragments emerging from the simulations and a simple coalescence model.
While both SMF and BLOB have been developed to simulate heavy ion interactions, they show very good results in reproducing the experimental yields of light fragments, up to alpha particles, obtained in the interaction of 12C with a thin carbon target at 62 MeV/u.
BLOB in particular gives promising results and this stresses the importance of integrating it into the Geant4 toolkit.
ATLAS's current software framework, Gaudi Athena, has been very successful for the experiment in LHC Runs 1 and 2. However, its single-threaded design has been recognised for some time to be ...increasingly problematic as CPUs have increased core counts and decreased available memory per core. Even the multi-process version of Athena, AthenaMP, will not scale to the range of architectures we expect to use beyond Run2. ATLAS examined the requirements on an updated multi-threaded framework and laid out plans for a new framework, including better support for High Level Trigger use cases, in 2014. In this paper we report on our progress in developing the new multi-threaded task parallel extension of Athena, AthenaMT. Implementing AthenaMT has required many significant code changes. Progress has been made in updating key concepts of the framework, allowing different levels of thread safety in algorithmic code. Substantial advances have also been made in implementing a data flow centric design, as well as on the development of the new 'event views' infrastructure. These event views support partial event processing and are an essential component to support the High Level Trigger's processing of certain regions of interest. A major effort has also been invested to have an early version of AthenaMT that can run simulation on many core architectures, which has augmented the understanding gained from work on earlier ATLAS demonstrators.
As Moore's Law continues to deliver more and more transistors, the mainstream processor industry is preparing to expand its investments in areas other than simple core count. These new interests ...include deep integration of on-chip components, advanced vector units, memory, cache and interconnect technologies. We examine these moving trends with parallelized and vectorized High Energy Physics workloads in mind. In particular, we report on practical experience resulting from experiments with scalable HEP benchmarks on the Intel "Ivy Bridge-EP" and "Haswell" processor families. In addition, we examine the benefits of the new "Haswell" microarchitecture and its impact on multiple facets of HEP software. Finally, we report on the power efficiency of new systems.
This paper summarizes the five years of CERN openlab's efforts focused on the Intel Xeon Phi co-processor, from the time of its inception to public release. We consider the architecture of the device ...vis a vis the characteristics of HEP software and identify key opportunities for HEP processing, as well as scaling limitations. We report on improvements and speedups linked to parallelization and vectorization on benchmarks involving software frameworks such as Geant4 and ROOT. Finally, we extrapolate current software and hardware trends and project them onto accelerators of the future, with the specifics of offline and online HEP processing in mind.
We present an overview of recent improvements of hadronic models in Geant4 for the physics configurations (Physics Lists) relevant to applications in high energy experiments. During last year the ...improvements have concentrated on the study of unphysical discontinuities in calorimeter observables in the transition regions between the models used in Physics Lists. The microscopic origin of these have been investigated, and possible improvements of Geant4 code are currently under validation. In this paper we discuss the status of the latest version of Geant4 with emphasis on the most promising new developments, namely the Fritiof based and CHIPS Physics Lists.