|
Jordan, D., Algora, A., & Tain, J. L. (2016). An event generator for simulations of complex beta-decay experiments. Nucl. Instrum. Methods Phys. Res. A, 828, 52–57.
Abstract: This article describes a Monte Carlo event generator for the design, optimization and performance characterization of beta decay spectroscopy experimental set-ups. The event generator has been developed within the Geant4 simulation architecture and provides new features and greater flexibility in comparison with the current available decay generator.
|
|
|
AGATA Collaboration, Domingo-Pardo, C., Bazzacco, D., Doornenbal, P., Farnea, E., Gadea, A., et al. (2012). Conceptual design and performance study for the first implementation of AGATA at the in-flight RIB facility of GSI. Nucl. Instrum. Methods Phys. Res. A, 694, 297–312.
Abstract: The main objective of the Advanced GAmma Tracking Array (AGATA) is the investigation of the structure of exotic nuclei at the new generation of RIB facilities. As part of the preparatory phase for FAIR-NUSTAR, AGATA is going to be installed at the FRS fragmentation facility of the GSI centre for an experimental campaign to be performed in 2012 and 2013. Owing to its gamma-ray tracking capabilities and the envisaged enhancement in resolving power, a series of in-flight gamma-ray spectroscopy experiments are being planned. The present work describes the conceptual design of this first implementation of AGATA at GSI-FRS, and provides information about the expected performance figures. According to the characteristics of each particular experiment, it is foreseen that the target-array distance is adjusted in order to achieve the optimum compromise between detection efficiency and energy resolution, or to cover an specific angular range of the emitted electromagnetic radiation. Thus, a comprehensive Monte Carlo study of the detection sensitivity in terms of photopeak efficiency, resolution and peak-to-total ratio, as a function of the target-array distance is presented. Several configurations have been investigated, and MC-calculations indicate that a remarkable enhancement in resolving power can be achieved when double-cluster AGATA detectors are developed and implemented. Several experimental effects are also investigated. This concerns the impact of passive materials between the target and the array, the angular distribution of the detection efficiency and the influence of target thickness effects and transition lifetimes in the attainable detection sensitivity. A short overview on half-life measurements via lineshape effects utilizing AGATA is also presented. (C) 2012 Elsevier B.V. All rights reserved.
|
|
|
Jaworski, G., Palacz, M., Nyberg, J., de Angelis, G., de France, G., Di Nitto, A., et al. (2012). Monte Carlo simulation of a single detector unit for the neutron detector array NEDA. Nucl. Instrum. Methods Phys. Res. A, 673, 64–72.
Abstract: A study of the dimensions and performance of a single detector of the future neutron detector array NEDA was performed by means of Monte Carlo simulations, using GEANT4. Two different liquid scintillators were evaluated: the hydrogen based BC501A and the deuterated BC537. The efficiency and the probability that one neutron will trigger a signal in more than one detector were investigated as a function of the detector size. The simulations were validated comparing the results to experimental measurements performed with two existing neutron detectors, with different geometries, based on the liquid scintillator BC501.
|
|
|
Guerrero, C., Cano-Ott, D., Mendoza, E., Tain, J. L., Algora, A., Berthoumieux, E., et al. (2012). Monte Carlo simulation of the n_TOF Total Absorption Calorimeter. Nucl. Instrum. Methods Phys. Res. A, 671, 108–117.
Abstract: The n_TOF Total Absorption Calorimeter (TAC) is a 4 pi BaF2 segmented detector used at CERN for measuring neutron capture cross-sections of importance for the design of advanced nuclear reactors. This work presents the simulation code that has been developed in GEANT4 for the accurate determination of the detection efficiency of the TAC for neutron capture events. The code allows to calculate the efficiency of the TAC for every neutron capture state, as a function of energy, crystal multiplicity, and counting rate. The code includes all instrumental effects such as the single crystal detection threshold and energy resolution, finite size of the coincidence time window, and signal pile-up. The results from the simulation have been validated with experimental data for a large set of electromagnetic de-excitation patterns: beta-decay of well known calibration sources, neutron capture reactions in light nuclei with well known level schemes like Ti-nat, reference samples used in (n,gamma) measurements like Au-197 and experimental data from an actinide sample like Pu-240. The systematic uncertainty in the determination of the detection efficiency has been estimated for all the cases. As a representative example, the accuracy reached for the case of Au-197(n,gamma) ranges between 0.5% and 2%, depending on the experimental and analysis conditions. Such a value matches the high accuracy required for the nuclear cross-section data needed in advanced reactor design.
|
|
|
Langer, C., Algora, A., Couture, A., Csatlos, M., Gulyas, J., Heil, M., et al. (2011). Simulations and developments of the Low Energy Neutron detector Array LENA. Nucl. Instrum. Methods Phys. Res. A, 659(1), 411–418.
Abstract: Prototypes of the Low Energy Neutron detector Array (LENA) have been tested and compared with detailed GEANT simulations. LENA will consist of plastic scintillation bars with the dimensions 1000 x 45 x 10 mm(3). The tests have been performed with gamma-ray sources and neutrons originating from the neutron-induced fission of (235)U. The simulations agreed very well with the measured response and were therefore used to simulate the response to mono-energetic neutrons with different detection thresholds. LENA will be used to detect low-energy neutrons from (p,n)-type reactions with low momentum transfer foreseen at the R(3)B and EXL setups at FAIR, Darmstadt.
|
|
|
AGATA Collaboration(Soderstrom, P. A. et al), & Gadea, A. (2011). Interaction position resolution simulations and in-beam measurements of the AGATA HPGe detectors. Nucl. Instrum. Methods Phys. Res. A, 638(1), 96–109.
Abstract: The interaction position resolution of the segmented HPGe detectors of an AGATA triple cluster detector has been studied through Monte Carlo simulations and in an in-beam experiment. A new method based on measuring the energy resolution of Doppler-corrected gamma-ray spectra at two different target to detector distances is described. This gives the two-dimensional position resolution in the plane perpendicular to the direction of the emitted gamma-ray. The gamma-ray tracking was used to determine the full energy of the gamma-rays and the first interaction point, which is needed for the Doppler correction. Five different heavy-ion induced fusion-evaporation reactions and a reference reaction were selected for the simulations. The results of the simulations show that the method works very well and gives a systematic deviation of <1 mm in the FVVHM of the interaction position resolution for the gamma-ray energy range from 60 keV to 5 MeV. The method was tested with real data from an in-beam measurement using a (30)5i beam at 64 MeV on a thin C-12 target. Pulse-shape analysis of the digitized detector waveforms and gamma-ray tracking was performed to determine the position of the first interaction point, which was used for the Doppler corrections. Results of the dependency of the interaction position resolution on the gamma-ray energy and on the energy, axial location and type of the first interaction point, are presented. The FVVHM of the interaction position resolution varies roughly linearly as a function of gamma-ray energy from 8.5 mm at 250 key to 4 mm at 1.5 MeV, and has an approximately constant value of about 4 mm in the gamma-ray energy range from 1.5 to 4 MeV.
|
|
|
AGATA Collaboration, Farnea, E., Recchia, F., Bazzacco, D., Kroll, T., Podolyak, Z., et al. (2010). Conceptual design and Monte Carlo simulations of the AGATA array. Nucl. Instrum. Methods Phys. Res. A, 621(1-3), 331–343.
Abstract: The aim of the Advanced GAmma Tracking Array (AGATA) project is the construction of an array based on the novel concepts of pulse shape analysis and gamma-ray tracking with highly segmented Ge semiconductor detectors. The conceptual design of AGATA and its performance evaluation under different experimental conditions has required the development of a suitable Monte Carlo code. In this article, the description of the code as well as simulation results relevant for AGATA, are presented.
|
|
|
ATLAS Collaboration(Adragna, P. et al), Castelo, J., Castillo Gimenez, V., Cuenca, C., Ferrer, A., Fullana, E., et al. (2010). Measurement of pion and proton response and longitudinal shower profiles up to 20 nuclear interaction lengths with the ATLAS Tile calorimeter. Nucl. Instrum. Methods Phys. Res. A, 615(2), 158–181.
Abstract: The response of pions and protons in the energy range of 20-180 GeV, produced at CERN's SPS H8 test-beam line in the ATLAS iron-scintillator Tile hadron calorimeter, has been measured. The test-beam configuration allowed the measurement of the longitudinal shower development for pions and protons up to 20 nuclear interaction lengths. It was found that pions penetrate deeper in the calorimeter than protons. However, protons induce showers that are wider laterally to the direction of the impinging particle. Including the measured total energy response, the pion-to-proton energy ratio and the resolution, all observations are consistent with a higher electromagnetic energy fraction in pion-induced showers. The data are compared with GEANT4 simulations using several hadronic physics lists. The measured longitudinal shower profiles are described by an analytical shower parametrization within an accuracy of 5-10%. The amount of energy leaking out behind the calorimeter is determined and parametrized as a function of the beam energy and the calorimeter depth. This allows for a leakage correction of test-beam results in the standard projective geometry.
|
|
|
Oliver, S., Rodriguez Bosca, S., & Gimenez-Alventosa, V. (2024). Enabling particle transport on CAD-based geometries for radiation simulations with penRed. Comput. Phys. Commun., 298, 109091–11pp.
Abstract: Geometry construction is a fundamental aspect of any radiation transport simulation, regardless of the Monte Carlo code being used. Typically, this process is tedious, time-consuming, and error-prone. The conventional approach involves defining geometries using mathematical objects or surfaces. However, this method comes with several limitations, especially when dealing with complex models, particularly those with organic shapes. Furthermore, since each code employs its own format and methodology for defining geometries, sharing and reproducing simulations among researchers becomes a challenging task. Consequently, many codes have implemented support for simulating over geometries constructed via Computer-Aided Design (CAD) tools. Unfortunately, this feature is lacking in penRed and other PENELOPE physics-based codes. Therefore, the objective of this work is to implement such support within the penRed framework. New version program summary Program Title: Parallel Engine for Radiation Energy Deposition (penRed) CPC Library link to program files: https://doi.org/10.17632/rkw6tvtngy.2 Developer's repository link: https://github.com/PenRed/PenRed Code Ocean capsule: https://codeocean.com/capsule/1041417/tree Licensing provisions: GNU Affero General Public License v3 Programming language: C++ standard 2011. Journal reference of previous version: V. Gimenez-Alventosa, V. Gimenez Gomez, S. Oliver, PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE, Computer Physics Communications 267 (2021) 108065. doi:https://doi.org/10.1016/j.cpc.2021.108065. Does the new version supersede the previous version?: Yes Reasons for the new version: Implements the capability to simulate on CAD constructed geometries, among many other features and fixes. Summary of revisions: All changes applied through the code versions are summarized in the file CHANGELOG.md in the repository package. Nature of problem: While Monte Carlo codes have proven valuable in simulating complex radiation scenarios, they rely heavily on accurate geometrical representations. In the same way as many other Monte Carlo codes, penRed employs simple geometric quadric surfaces like planes, spheres and cylinders to define geometries. However, since these geometric models offer a certain level of flexibility, these representations have limitations when it comes to simulating highly intricate and irregular shapes. Anatomic structures, for example, require detailed representations of organs, tissues and bones, which are difficult to achieve using basic geometric objects. Similarly, complex devices or intricate mechanical systems may have designs that cannot be accurately represented within the constraints of such geometric models. Moreover, when the complexity of the model increases, geometry construction process becomes more difficult, tedious, time-consuming and error-prone [2]. Also, as each Monte Carlo geometry library uses its own format and construction method, reproducing the same geometry among different codes is a challenging task. Solution method: To face the problems stated above, the objective of this work is to implement the capability to simulate using irregular and adaptable meshed geometries in the penRed framework. This kind of meshes can be constructed using Computer-Aided Design (CAD) tools, the use of which is very widespread and streamline the design process. This feature has been implemented in a new geometry module named “MESH_BODY” specific for this kind of geometries. This one is freely available and usable within the official penRed package1. It can be used since penRed version 1.9.3b and above.
|
|
|
Gimenez-Alventosa, V., Gimenez, V., & Oliver, S. (2021). PenRed: An extensible and parallel Monte-Carlo framework for radiation transport based on PENELOPE. Comput. Phys. Commun., 267, 108065–12pp.
Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.
|
|