|
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.
|
|
|
Barrientos, L., Borja-Lloret, M., Etxebeste, A., Muñoz, E., Oliver, J. F., Ros, A., et al. (2021). Performance evaluation of MACACO II Compton camera. Nucl. Instrum. Methods Phys. Res. A, 1014, 165702–7pp.
Abstract: The IRIS group at IFIC-Valencia has developed a second version of a Compton camera prototype for hadron therapy treatment monitoring, with the aim of improving the performance with respect to its predecessor. The system is composed of three Lanthanum (III) bromide (LaBr3) crystals coupled to silicon photomultipliers (SiPMs). The detector energy resolution has been improved to 5.6% FWHM at 511 keV and an angular resolution of 8.0 degrees has been obtained. Images of a Na-22 point-like source have been reconstructed selecting two and three interaction events. Moreover, the experimental data have been reproduced with Monte Carlo simulations using a Compton camera module (CCMod) in GATE v8.2 obtaining a good correlation.
|
|
|
Wieduwilt, P., Paschen, B., Schreeck, H., Schwenker, B., Soltau, J., Ahlburg, P., et al. (2021). Performance of production modules of the Belle II pixel detector in a high-energy particle beam. Nucl. Instrum. Methods Phys. Res. A, 991, 164978–15pp.
Abstract: The Belle II experiment at the Super B factory SuperKEKB, an asymmetric e(+) e(-) collider located in Tsukuba, Japan, is tailored to perform precision B physics measurements. The centre of mass energy of the collisions is equal to the rest mass of the gamma (4S) resonance of m(gamma(4S)) = 10.58 GeV. A high vertex resolution is essential for measuring the decay vertices of B mesons. Typical momenta of the decay products are ranging from a few tens of MeV to a few GeV and multiple scattering has a significant impact on the vertex resolution. The VerteX Detector (VXD) for Belle II is therefore designed to have as little material as possible inside the acceptance region. Especially the innermost two layers, populated by the PiXel Detector (PXD), have to be ultra-thin. The PXD is based on DEpleted P-channel Field Effect Transistors (DEPFETs) with a thickness of only 75 μm. Spatial resolution and hit efficiency of production detector modules were studied in beam tests performed at the DESY test beam facility. The spatial resolution was investigated as a function of the incidence angle and improvements due to charge sharing are demonstrated. The measured module performance is compatible with the requirements for Belle II.
|
|
|
ATLAS Collaboration(Aad, G. et al), Aparisi Pozo, J. A., Bailey, A. J., Cabrera Urban, S., Cardillo, F., Castillo, F. L., et al. (2021). Performance of the ATLAS RPC detector and Level-1 muon barrel trigger at root s=13 TeV. J. Instrum., 16(7), P07029–64pp.
Abstract: The ATLAS experiment at the Large Hadron Collider (LHC) employs a trigger system consisting of a first-level hardware trigger (L1) and a software-based high-level trigger. The L1 muon trigger system selects muon candidates, assigns them to the correct LHC bunch crossing and classifies them into one of six transverse-momentum threshold classes. The L1 muon trigger system uses resistive-plate chambers (RPCs) to generate the muon-induced trigger signals in the central (barrel) region of the ATLAS detector. The ATLAS RPCs are arranged in six concentric layers and operate in a toroidal magnetic field with a bending power of 1.5 to 5.5 Tm. The RPC detector consists of about 3700 gas volumes with a total surface area of more than 4000 m(2). This paper reports on the performance of the RPC detector and L1 muon barrel trigger using 60.8 fb(-1) of proton-proton collision data recorded by the ATLAS experiment in 2018 at a centre-of-mass energy of 13 TeV. Detector and trigger performance are studied using Z boson decays into a muon pair. Measurements of the RPC detector response, efficiency, and time resolution are reported. Measurements of the L1 muon barrel trigger efficiencies and rates are presented, along with measurements of the properties of the selected sample of muon candidates. Measurements of the RPC currents, counting rates and mean avalanche charge are performed using zero-bias collisions. Finally, RPC detector response and efficiency are studied at different high voltage and front-end discriminator threshold settings in order to extrapolate detector response to the higher luminosity expected for the High Luminosity LHC.
|
|
|
Escrig, S. et al, & Morales, A. I. (2021). Persistence of the Z=28 shell gap in A=75 isobars: Identification of a possible (1/2(-)) μs isomer in Co-75 and beta decay to Ni-75. Phys. Rev. C, 103(6), 064328–12pp.
Abstract: Background: The evolution of shell structure around doubly magic exotic nuclei is of great interest in nuclear physics and astrophysics. In the 'southwest' region of Ni-78, the development of deformation might trigger a major shift in our understanding of explosive nucleosynthesis. To this end, new spectroscopic information on key close-lying nuclei is very valuable. Purpose: We intend to measure the isomeric and beta decay of Co-75, with one-proton and two-neutron holes relative to Ni-78, to access new nuclear structure information in Co-75 and its beta-decay daughters Ni-75 and Ni-74. Methods: The nucleus Co-75 is produced in relativistic in-flight fission reactions of U-238 at the Radioactive Ion Beam Factory in the RIKEN Nishina Center. Its isomeric and f decay are studied exploiting the BigRIPS and EURICA setups. Results: We obtain partial beta-decay spectra for Ni-75 and Ni-74, and report a new isomeric transition in Co-75. The energy [E-gamma = 1914(2) keV] and half-life [t(1/2) = 13(6) μs] of the delayed gamma ray lend support for the existence of aJ(pi) = (1/2(-)) isomeric state at 1914(2) keV. A comparison with PFSDG-U shell-model calculations provides a good account for the observed states in Ni-75, but the first calculated 1/2(-) level in Co-75, a prolate K = 1/2 state, is predicted about 1 MeV below the observed (1/2(-)) level. Conclusions: The spherical-like structure of the lowest-lying excited states in Ni-75 is proved. In the case of Co-75, the results suggest that the dominance of the spherical configurations over the deformed ones might be stronger than expected below Ni-78. Further experimental efforts to discern the nature of the J(pi) = (1/2(-)) isomer are necessary.
|
|