|
Ballester, F., Tedgren, A. C., Granero, D., Haworth, A., Mourtada, F., Fonseca, G. P., et al. (2015). A generic high-dose rate Ir-192 brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism. Med. Phys., 42(6), 3048–3062.
Abstract: Purpose: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) Ir-192 source and a virtual water phantom were designed, which can be imported into a TPS. Methods: A hypothetical, generic HDR Ir-192 source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic Ir-192 source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra (R) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision AcuRos (TM)]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and pENELopE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR Ir-192 source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. Results: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ACE algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 +/- 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ACE at clinically relevant distances. Conclusions: A hypothetical, generic HDR Ir-192 source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.
|
|
|
Assam, I., Vijande, J., Ballester, F., Perez-Calatayud, J., Poppe, B., & Siebert, F. A. (2022). Evaluation of dosimetric effects of metallic artifact reduction and tissue assignment on Monte Carlo dose calculations for I-125 prostate implants. Med. Phys., 49, 6195–6208.
Abstract: Purpose Monte Carlo (MC) simulation studies, aimed at evaluating the magnitude of tissue heterogeneity in I-125 prostate permanent seed implant brachytherapy (BT), customarily use clinical post-implant CT images to generate a virtual representation of a realistic patient model (virtual patient model). Metallic artifact reduction (MAR) techniques and tissue assignment schemes (TAS) are implemented on the post-implant CT images to mollify metallic artifacts due to BT seeds and to assign tissue types to the voxels corresponding to the bright seed spots and streaking artifacts, respectively. The objective of this study is to assess the combined influence of MAR and TAS on MC absorbed dose calculations in post-implant CT-based phantoms. The virtual patient models used for I-125 prostate implant MC absorbed dose calculations in this study are derived from the CT images of an external radiotherapy prostate patient without BT seeds and prostatic calcifications, thus averting the need to implement MAR and TAS. Methods The geometry of the IsoSeed I25.S17plus source is validated by comparing the MC calculated results of the TG-43 parameters for the line source approximation with the TG-43U1S2 consensus data. Four MC absorbed dose calculations are performed in two virtual patient models using the egs_brachy MC code: (1) TG-43-based D-w,w-TG(43), (2) D-w,D-w-MBDC that accounts for interseed scattering and attenuation (ISA), (3) D-m,D-m that examines ISA and tissue heterogeneity by scoring absorbed dose in tissue, and (4) D-w,D-m that unlike D-m,D-m scores absorbed dose in water. The MC absorbed doses (1) and (2) are simulated in a TG-43 patient phantom derived by assigning the densities of every voxel to 1.00 g cm(-3) (water), whereas MC absorbed doses (3) and (4) are scored in the TG-186 patient phantom generated by mapping the mass density of each voxel to tissue according to a CT calibration curve. The MC absorbed doses calculated in this study are compared with VariSeed v8.0 calculated absorbed doses. To evaluate the dosimetric effect of MAR and TAS, the MC absorbed doses of this work (independent of MAR and TAS) are compared to the MC absorbed doses of different I-125 source models from previous studies that were calculated with different MC codes using post-implant CT-based phantoms generated by implementing MAR and TAS on post-implant CT images. Results The very good agreement of TG-43 parameters of this study and the published consensus data within 3% validates the geometry of the IsoSeed I25.S17plus source. For the clinical studies, the TG-43-based calculations show a D-90 overestimation of more than 4% compared to the more realistic MC methods due to ISA and tissue composition. The results of this work generally show few discrepancies with the post-implant CT-based dosimetry studies with respect to the D-90 absorbed dose metric parameter. These discrepancies are mainly Type B uncertainties due to the different I-125 source models and MC codes. Conclusions The implementation of MAR and TAS on post-implant CT images have no dosimetric effect on the I-125 prostate MC absorbed dose calculation in post-implant CT-based phantoms.
|
|
|
Asai, M., Cortes-Giraldo, M. A., Gimenez-Alventosa, V., Gimenez, V., & Salvat, F. (2021). The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4. Front. Physics, 9, 738735–20pp.
Abstract: A translation of the penelope physics subroutines to C++, designed as an extension of the Geant4 toolkit, is presented. The Fortran code system penelope performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. Penelope implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of Geant4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original penelope programs. The combined code, named PenG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of Geant4.
|
|
|
n_TOF Collaboration(Alcayne, V. et al), Balibrea-Correa, J., Domingo-Pardo, C., Lerendegui-Marco, J., Babiano-Suarez, V., & Ladarescu, I. (2024). A Segmented Total Energy Detector (sTED) optimized for (n,γ) cross-section measurements at n_TOF EAR2. Radiat. Phys. Chem., 217, 11pp.
Abstract: The neutron time-of-flight facility nTOF at CERN is a spallation source dedicated to measurements of neutroninduced reaction cross-sections of interest in nuclear technologies, astrophysics, and other applications. Since 2014, Experimental ARea 2 (EAR2) is operational and delivers a neutron fluence of similar to 4 center dot 10(7) neutrons per nominal proton pulse, which is similar to 50 times higher than the one of Experimental ARea 1 (EAR1) of similar to 8 center dot 10(5) neutrons per pulse. The high neutron flux at EAR2 results in high counting rates in the detectors that challenged the previously existing capture detection systems. For this reason, a Segmented Total Energy Detector (sTED) has been developed to overcome the limitations in the detector's response, by reducing the active volume per module and by using a photo-multiplier (PMT) optimized for high counting rates. This paper presents the main characteristics of the sTED, including energy and time resolution, response to gamma-rays, and provides as well details of the use of the Pulse Height Weighting Technique (PHWT) with this detector. The sTED has been validated to perform neutron-capture cross-section measurements in EAR2 in the neutron energy range from thermal up to at least 400 keV. The detector has already been successfully used in several measurements at nTOF EAR2.
|
|
|
ATLAS Collaboration(Aad, G. et al), Amoros, G., Cabrera Urban, S., Castillo Gimenez, V., Costa, M. J., Escobar, C., et al. (2011). Measurement of the transverse momentum distribution of Z/gamma* bosons in proton-proton collisions at sqrt(s)=7 TeV with the ATLAS detector. Phys. Lett. B, 705(5), 415–434.
Abstract: A measurement of the Z/gamma* transverse momentum (p(T)(Z)) distribution in proton-proton collisions at sqrt(s) = 7 TeV is presented using Z/gamma* -> e(+)e(-) and Z/gamma* -> mu(+)mu(-) decays collected with the ATLAS detector in data sets with integrated luminosities of 35 pb(-1) and 40 pb(-1), respectively. The normalized differential cross sections are measured separately for electron and muon decay channels as well as for their combination up to pi of 350 GeV for invariant dilepton masses 66 GeV < m(ll) < 116 GeV. The measurement is compared to predictions of perturbative QCD and various event generators. The prediction of resummed QCD combined with fixed order perturbative QCD is found to be in good agreement with the data.
|
|
|
KM3NeT Collaboration(Aiello, S. et al), Alves Garre, S., Calvo, D., Carretero, V., Colomer, M., Corredoira, I., et al. (2020). gSeaGen: The KM3NeT GENIE-based code for neutrino telescopes. Comput. Phys. Commun., 256, 107477–15pp.
Abstract: The gSeaGen code is a GENIE-based application developed to efficiently generate high statistics samples of events, induced by neutrino interactions, detectable in a neutrino telescope. The gSeaGen code is able to generate events induced by all neutrino flavours, considering topological differences between tracktype and shower-like events. Neutrino interactions are simulated taking into account the density and the composition of the media surrounding the detector. The main features of gSeaGen are presented together with some examples of its application within the KM3NeT project. Program summary Program Title: gSeaGen CPC Library link to program files: http://dx.doi.org/10.17632/ymgxvy2br4.1 Licensing provisions: GPLv3 Programming language: C++ External routines/libraries: GENIE [1] and its external dependencies. Linkable to MUSIC [2] and PROPOSAL [3]. Nature of problem: Development of a code to generate detectable events in neutrino telescopes, using modern and maintained neutrino interaction simulation libraries which include the state-of-the-art physics models. The default application is the simulation of neutrino interactions within KM3NeT [4]. Solution method: Neutrino interactions are simulated using GENIE, a modern framework for Monte Carlo event generators. The GENIE framework, used by nearly all modern neutrino experiments, is considered as a reference code within the neutrino community. Additional comments including restrictions and unusual features: The code was tested with GENIE version 2.12.10 and it is linkable with release series 3. Presently valid up to 5 TeV. This limitation is not intrinsic to the code but due to the present GENIE valid energy range. References: [1] C. Andreopoulos at al., Nucl. Instrum. Meth. A614 (2010) 87. [2] P. Antonioli et al., Astropart. Phys. 7 (1997) 357. [3] J. H. Koehne et al., Comput. Phys. Commun. 184 (2013) 2070. [4] S. Adrian-Martinez et al., J. Phys. G: Nucl. Part. Phys. 43 (2016) 084001.
|
|
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
|
|
ATLAS Collaboration(Aad, G. et al), Cabrera Urban, S., Castillo Gimenez, V., Costa, M. J., Fassi, F., Ferrer, A., et al. (2013). Measurement of angular correlations in Drell-Yan lepton pairs to probe Z/gamma* boson transverse momentum at root s=7 TeV with the ATLAS detector. Phys. Lett. B, 720(1-3), 32–51.
Abstract: A measurement of angular correlations in Drell-Yan lepton pairs via the phi(eta)* observable is presented. This variable probes the same physics as the Z/gamma* boson transverse momentum with a better experimental resolution. The Z/gamma* -> e(+)e(-) and Z/gamma* -> mu(+)mu(-) decays produced in proton-proton collisions at a centre-of-mass energy of root s = 7 TeV are used. The data were collected with the ATLAS detector at the LHC and correspond to an integrated luminosity of 4.6 fb(-1). Normalised differential cross sections as a function of phi(eta)* are measured separately for electron and muon decay channels. These channels are then combined for improved accuracy. The cross section is also measured double differentially as a function of phi(eta)* for three independent bins of the Z boson rapidity. The results are compared to QCD calculations and to predictions from different Monte Carlo event generators. The data are reasonably well described, in all measured Z boson rapidity regions, by resummed QCD predictions combined with fixed-order perturbative QCD calculations or by some Monte Carlo event generators. The measurement precision is typically better by one order of magnitude than present theoretical uncertainties.
|
|