Maria Grazia Pia

New Developments of the Goodness-of-Fit Statistical Toolkit

The Statistical Toolkit is a project for the development of open source software tools for statistical data analysis in experimental particle and nuclear physics. The second development cycle encompassed an extension of the software functionality and new tools to facilitate its usage in experimental environments. The new developments include additional goodness-of-fit tests, new implementations of existing tests to improve their statistical precision or computational performance, a new component to extend the usability of the toolkit with other data analysis systems, and new tools for an easier configuration and build of the system in the users computing environment. The computational performance of all the algorithms implemented has been studied

Validation of Geant4 Simulation of Electron Energy Deposition

Geant4-based simulations of the energy deposited by electrons in various materials are quantitatively compared to high-precision calorimetric measurements taken at Sandia Laboratories. The experimental data concern electron beams of energy between a few tens of kilolectron volt and 1 MeV at various incidence angles. Two experimental scenarios are evaluated: the longitudinal energy deposition pattern in a finely segmented detector, and the total energy deposited in a larger size calorimeter. The simulations are produced with Geant4 versions from 9.1 to 9.6; they involve models of electron-photon interactions in the standard and low energy electromagnetic packages, and various implementations of electron multiple scattering. Significant differences in compatibility with experimental data are observed in the longitudinal energy deposition patterns produced by the examined Geant4 versions, while the total deposited energy exhibits smaller variations across the various Geant4 versions, with the exception Geant4 9.4. The validation analysis, based on statistical methods, shows that the best compatibility between simulation and experimental energy deposition profiles is achieved using electromagnetic models based on the EEDL and EPDL evaluated data libraries with Geant4 9.1. The results document the accuracy achievable in the simulation of the energy deposited by low energy electrons with Geant4; they provide guidance for application in similar experimental scenarios and for improving Geant4.

Validation Test of Geant4 Simulation of Electron Backscattering

Backscattering is a sensitive probe of the accuracy of electron scattering algorithms implemented in Monte Carlo codes. The capability of the Geant4 toolkit to describe realistically the fraction of electrons backscattered from a target volume is extensively and quantitatively evaluated in comparison with experimental data retrieved from the literature. The validation test covers the energy range between approximately 100 eV and 20 MeV, and concerns a wide set of target elements. Multiple and single electron scattering models implemented in Geant4, as well as preassembled selections of physics models distributed within Geant4, are analyzed with statistical methods. The evaluations concern Geant4 versions from 9.1 to 10.1. Significant evolutions are observed over the range of Geant4 versions, not always in the direction of better compatibility with experiment. Goodness-of-fit tests complemented by categorical analysis tests identify a configuration based on Geant4 Urban multiple scattering model in Geant4 version 9.1 and a configuration based on single Coulomb scattering in Geant4 10.0 as the physics options best reproducing experimental data above a few tens of keV. At lower energies only single scattering demonstrates some capability to reproduce data down to a few keV. Recommended preassembled physics configurations appear incapable of describing electron backscattering compatible with experiment. With the support of statistical methods, a correlation is established between the validation of Geant4-based simulation of backscattering and of energy deposition.

Investigation of Geant4 Simulation of Electron Backscattering

A test of Geant4 simulation of electron backscattering recently published in this journal prompted further investigation into the causes of the observed behaviour. An interplay between features of geometry and physics algorithms implemented in Geant4 is found to significantly affect the accuracy of backscattering simulation in some physics configurations.

An exact framework for uncertainty quantification in Monte Carlo simulation

In the context of Monte Carlo (MC) simulation of particle transport Uncertainty Quantification (UQ) addresses the issue of predicting non statistical errors affecting the physical results, i.e. errors deriving mainly from uncertainties in the physics data and/or in the model they embed. In the case of a single uncertainty a simple analytical relation exists among its the Probability Density Function (PDF) and the corresponding PDF for the output of the simulation: this allows a complete statistical analysis of the results of the simulation. We examine the extension of this result to the multi-variate case, when more than one of the physical input parameters are affected by uncertainties: a typical scenario is the prediction of the dependence of the simulation on input cross section tabulations.

Algorithms and parameters for improved accuracy in physics data libraries

Recent efforts for the improvement of the accuracy of physics data libraries used in particle transport are summarized. Results are reported about a large scale validation analysis of atomic parameters used by major Monte Carlo systems (Geant4, EGS, MCNP, Penelope etc.); their contribution to the accuracy of simulation observables is documented. The results of this study motivated the development of a new atomic data management software package, which optimizes the provision of state-of-the-art atomic parameters to physics models. The effect of atomic parameters on the simulation of radioactive decay is illustrated. Ideas and methods to deal with physics models applicable to different energy ranges in the production of data libraries, rather than at runtime, are discussed.

Geant4 maintainability assessed with respect to software engineering references

We report a methodology developed to quantitatively assess the maintainability of Geant4 with respect to software engineering references. The level of maintainability is determined by combining a set of metrics values whose references are documented in literature.

Testable physics by design

The ability to test scientific software needs to be supported by adequate software design. Legacy software systems are often characterized by the difficulty to test parts of the software, mainly due to existing dependencies on other parts. Methods to improve the testability of physics software are discussed, along with open issues specific to physics software for Monte Carlo particle transport. The discussion is supported by examples drawn from the experience with validating Geant4 physics.

Editorial Conference Comments by the Editors

The Symposium on Radiation Measurements and Applications (SORMA) convened for the third time on the West Coast, May 22–26, 2016, at the Clark Kerr Campus of the University of California, Berkeley, CA, USA. With radiation detectors increasing in number, variety, and societal importance, we are alternating between SORMA (in Ann Arbor, MI, USA) and SORMA West so that the forum will be available every two years.