|
|
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.
|
|
|
|
Amaldi, U., Bonomi, R., Braccini, S., Crescenti, M., Degiovanni, A., Garlasche, M., et al. (2010). Accelerators for hadrontherapy: From Lawrence cyclotrons to linacs. Nucl. Instrum. Methods Phys. Res. A, 620(2-3), 563–577.
Abstract: Hadrontherapy with protons and carbon ions is a fast developing methodology in radiation oncology. The accelerators used and planned for this purpose are reviewed starting from the cyclotrons used in the thirties. As discussed in the first part of this paper, normal and superconducting cyclotrons are still employed, together with synchrotrons, for proton therapy while for carbon ion therapy synchrotrons have been till now the only option. The latest developments concern a superconducting cyclotron for carbon ion therapy, fast-cycling high frequency linacs and 'single room' proton therapy facilities. These issues are discussed in the second part of the paper by underlining the present challenges, in particular the treatment of moving organs.
|
|
|
|
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.
|
|
|
|
Albiol, F., Corbi, A., & Albiol, A. (2017). 3D measurements in conventional X-ray imaging with RGB-D sensors. Med. Eng. Phys., 42, 73–79.
Abstract: A method for deriving 3D internal information in conventional X-ray settings is presented. It is based on the combination of a pair of radiographs from a patient and it avoids the use of X-ray-opaque fiducials and external reference structures. To achieve this goal, we augment an ordinary X-ray device with a consumer RGB-D camera. The patient' s rotation around the craniocaudal axis is tracked relative to this camera thanks to the depth information provided and the application of a modern surface-mapping algorithm. The measured spatial information is then translated to the reference frame of the X-ray imaging system. By using the intrinsic parameters of the diagnostic equipment, epipolar geometry, and X-ray images of the patient at different angles, 3D internal positions can be obtained. Both the RGB-D and Xray instruments are first geometrically calibrated to find their joint spatial transformation. The proposed method is applied to three rotating phantoms. The first two consist of an anthropomorphic head and a torso, which are filled with spherical lead bearings at precise locations. The third one is made of simple foam and has metal needles of several known lengths embedded in it. The results show that it is possible to resolve anatomical positions and lengths with a millimetric level of precision. With the proposed approach, internal 3D reconstructed coordinates and distances can be provided to the physician. It also contributes to reducing the invasiveness of ordinary X-ray environments and can replace other types of clinical explorations that are mainly aimed at measuring or geometrically relating elements that are present inside the patient's body.
|
|
|
|
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.
|
|
|
|
Pi, M., Barranco, M., Navarro, J., & Ancilotto, F. (2012). Nucleation and cavitation in parahydrogen. Chem. Phys., 399, 213–217.
Abstract: We have used a density functional approach to investigate thermal homogeneous nucleation and cavitation in parahydrogen. The effect of electrons as seeds of heterogeneous cavitation in liquid parahydrogen is also discussed within the capillary model.
|
|
|
|
Plaza, J., Martinez, T., Becares, V., Cano-Ott, D., Villamarin, D., de Rada, A. P., et al. (2023). Thermal neutron background at Laboratorio Subterraneo de Canfranc (LSC). Astropart Phys., 146, 102793–9pp.
Abstract: The thermal neutron background at Laboratorio Subterraneo de Canfranc (LSC) has been determined using several He-3 proportional counter detectors. Bare and Cd shielded counters were used in a series of long measurements. Pulse shape discrimination techniques were applied to discriminate between neutron and gamma signals as well as other intrinsic contributions. Montecarlo simulations allowed us to estimate the sensitivity of the detectors and calculate values for the background flux of thermal neutrons inside Hall-A of LSC. The obtained value is (3.5 +/- 0.8)x10(-6) n/cm(2)s, and is within an order of magnitude compared to similar facilities.
|
|
|
|
Jordan, D., Tain, J. L., Algora, A., Agramunt, J., Domingo-Pardo, C., Gomez-Hornillos, M. B., et al. (2013). Measurement of the neutron background at the Canfranc Underground Laboratory LSC. Astropart Phys., 42, 1–6.
Abstract: The energy distribution of the neutron background was measured for the first time at Hall A of the Canfranc Underground Laboratory. For this purpose we used a novel approach based on the combination of the information obtained with six large high-pressure He-3 proportional counters embedded in individual polyethylene blocks of different size. In this way not only the integral value but also the flux distribution as a function of neutron energy was determined in the range from 1 eV to 10 MeV. This information is of importance because different underground experiments show different neutron background energy dependence. The high sensitivity of the setup allowed to measure a neutron flux level which is about four orders of magnitude smaller that the neutron background at sea level. The integral value obtained is Phi(Hall A) = (3.44 +/- 0.35) x 10(-6) cm(-2) s(-1).
|
|
|
|
Menjo, H. et al, Faus-Golfe, A., & Velasco, J. (2011). Monte Carlo study of forward pi(0) production spectra to be measured by the LHCf experiment for the purpose of benchmarking hadron interaction models at 10(17) eV. Astropart Phys., 34(7), 513–520.
Abstract: The LHCf experiment aims to improve knowledge of forward neutral particle production spectra at the LHC energy which is relevant for the interpretation of air shower development of high energy cosmic rays. Two detectors, each composed of a pair of sampling and imaging calorimeters, have been installed at the forward region of IP1 to measure pi(0) energy spectra above 600 GeV. In this paper, we present a Monte Carlo study of the pi(0) measurements to be performed with one of the LHCf detectors for proton-proton collisions at root s = 14 TeV. In approximately 40 min of operation at luminosity 0.8 x 10(29) cm(-2) s(-1) during the beam commissioning phase of LHC, about 1.5 x 10(4) pi(0) events are expected to be obtained at two transverse positions of the detector. The backgrounds from interactions of secondary particles with beam pipes and interactions of beam particles with residual gas in the beam pipes are expected to be less than 0.1% of the signal from pi(0)s. We also discuss the capability of LHCf measurements to discriminate between the various hadron interaction models that are used for simulation of high energy air showers, such as DPMJET3.03, QGSJETII-03, SIBYLL2.1 and EPOS1.99.
|
|
