|
Gimenez-Alventosa, V., Gimenez, V., Ballester, F., Vijande, J., & Andreo, P. (2020). Monte Carlo calculation of beam quality correction factors for PTW cylindrical ionization chambers in photon beams. Phys. Med. Biol., 65(20), 205005–11pp.
Abstract: The beam quality correction factork(Q)for megavoltage photon beams has been calculated for eight PTW (Freiburg, Germany) ionization chambers (Farmer chambers PTW30010, PTW30011, PTW30012, and PTW30013, Semiflex 3D chambers PTW31021, PTW31010, and PTW31013, and the PinPoint 3D chamber PTW31016). Simulations performed on the widely used NE-2571 ionization chamber have been used to benchmark the results. The Monte Carlo code PENELOPE/penEasy was used to calculate the absorbed dose to a point in water and the absorbed dose to the active air volume of the chambers for photon beams in the range 4 to 24 MV. Of the nine ionization chambers analysed, only five are included in the current version of the International Code of Practice for dosimetry based on standards of absorbed dose to water (IAEA TRS 398). The values reported in this work agree with those in the literature within the uncertainty estimates and are to be included in the average values of the data obtained by different working groups for the forthcoming update of TRS 398.
|
|
|
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.
|
|
|
Gomez-Cadenas, J. J., Benlloch-Rodriguez, J. M., & Ferrario, P. (2017). Monte Carlo study of the coincidence resolving time of a liquid xenon PET scanner, using Cherenkov radiation. J. Instrum., 12, P08023–13pp.
Abstract: In this paper we use detailed Monte Carlo simulations to demonstrate that liquid xenon (LXe) can be used to build a Cherenkov-based TOF-PET, with an intrinsic coincidence resolving time (CRT) in the vicinity of 10 ps. This extraordinary performance is due to three facts: a) the abundant emission of Cherenkov photons by liquid xenon; b) the fact that LXe is transparent to Cherenkov light; and c) the fact that the fastest photons in LXe have wavelengths higher than 300 nm, therefore making it possible to separate the detection of scintillation and Cherenkov light. The CRT in a Cherenkov LXe TOF-PET detector is, therefore, dominated by the resolution (time jitter) introduced by the photosensors and the electronics. However, we show that for sufficiently fast photosensors (e.g, an overall 40 ps jitter, which can be achieved by current micro-channel plate photomultipliers) the overall CRT varies between 30 and 55 ps, depending on the detection efficiency. This is still one order of magnitude better than commercial CRT devices and improves by a factor 3 the best CRT obtained with small laboratory prototypes.
|
|
|
Guadilla, V. et al, Algora, A., Tain, J. L., Agramunt, J., Jordan, D., Monserrate, M., et al. (2017). Characterization of a cylindrical plastic beta-detector with Monte Carlo simulations of optical photons. Nucl. Instrum. Methods Phys. Res. A, 854, 134–138.
Abstract: In this work we report on the Monte Carlo study performed to understand and reproduce experimental measurements of a new plastic beta-detector with cylindrical geometry. Since energy deposition simulations differ from the experimental measurements for such a geometry, we show how the simulation of production and transport of optical photons does allow one to obtain the shapes of the experimental spectra. Moreover, taking into account the computational effort associated with this kind of simulation, we develop a method to convert the simulations of energy deposited into light collected, depending only on the interaction point in the detector. This method represents a useful solution when extensive simulations have to be done, as in the case of the calculation of the response function of the spectrometer in a total absorption gamma-ray spectroscopy analysis.
|
|
|
Guadilla, V., Algora, A., Estienne, M., Fallot, M., Gelletly, W., Porta, A., et al. (2024). First measurements with a new fl-electron detector for spectral shape studies. J. Instrum., 19(2), P02027–21pp.
Abstract: The shape of the electron spectrum emitted in /3 decay carries a wealth of information about nuclear structure and fundamental physics. In spite of that, few dedicated measurements have been made of /3 -spectrum shapes. In this work we present a newly developed detector for /3 electrons based on a telescope concept. A thick plastic scintillator is employed in coincidence with a thin silicon detector. The first measurements employing this detector have been carried out with mono -energetic electrons from the high-energy resolution electron -beam spectrometer at Bordeaux. Here we report on the good reproduction of the experimental spectra of mono -energetic electrons using Monte Carlo simulations. This is a crucial step for future experiments, where a detailed Monte Carlo characterization of the detector is needed to determine the shape of the /3 -electron spectra by deconvolution of the measured spectra with the response function of the detector. A chamber to contain two telescope assemblies has been designed for future /3 -decay experiments at the Ion Guide Isotope Separator On -Line facility in Jyvaskyla, aimed at improving our understanding of reactor antineutrino spectra.
|
|