|
|
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.
|
|
|
|
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.
|
|
|
|
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.
|
|
|
|
Vijande, J., Granero, D., Perez-Calatayud, J., & Ballester, F. (2013). Monte Carlo dosimetric study of the medium dose rate CSM40 source. Appl. Radiat. Isot., 82, 283–288.
Abstract: The Cs-137 medium dose rate (MDR) CSM40 source model (Eckert & Ziegler BEBIG, Germany) is in clinical use but no dosimetric dataset has been published. This study aims to obtain dosimetric data for the CSM40 source for its use in clinical practice as required by the American Association of Physicists in Medicine (AAPM) and the European Society for Radiotherapy and Oncology (ESTRO). Penelope2008 and Geant4 Monte Carlo codes were used to characterize this source dosimetrically. It was located in an unbounded water phantom with composition and mass density as recommended by AAPM and ESTRO. Due to the low photon energies of Cs-137, absorbed dose was approximated by collisional kerma. Additional simulations were performed to obtain the air-kerma strength, sic. Mass-energy absorption coefficients in water and air were consistently derived and used to calculate collisional kerma. Results performed with both radiation transport codes showed agreement typically within 0.05%. Dose rate constant, radial dose function and anisotropy function are provided for the CSM40 and compared with published data for other commercially available Cs-137 sources. An uncertainty analysis has been performed. The data provided by this study can be used as input data and verification in the treatment planning systems. (C) 2013 Elsevier Ltd. All rights reserved.
|
|
|
|
Granero, D., Perez-Calatayud, J., Vijande, J., Ballester, F., & Rivard, M. J. (2014). Limitations of the TG-43 formalism for skin high-dose-rate brachytherapy dose calculations. Med. Phys., 41(2), 021703–8pp.
Abstract: Purpose: In skin high-dose-rate (HDR) brachytherapy, sources are located outside, in contact with, or implanted at some depth below the skin surface. Most treatment planning systems use the TG-43 formalism, which is based on single-source dose superposition within an infinite water medium without accounting for the true geometry in which conditions for scattered radiation are altered by the presence of air. The purpose of this study is to evaluate the dosimetric limitations of the TG-43 formalism in HDR skin brachytherapy and the potential clinical impact. Methods: Dose rate distributions of typical configurations used in skin brachytherapy were obtained: a 5 cm x 5 cm superficial mould; a source inside a catheter located at the skin surface with and without backscatter bolus; and a typical interstitial implant consisting of an HDR source in a catheter located at a depth of 0.5 cm. Commercially available HDR Co-60 and Ir-192 sources and a hypothetical Yb-169 source were considered. The Geant4Monte Carlo radiation transport code was used to estimate dose rate distributions for the configurations considered. These results were then compared to those obtained with the TG-43 dose calculation formalism. In particular, the influence of adding bolus material over the implant was studied. Results: For a 5 cm x 5 cm Ir-192 superficial mould and 0.5 cm prescription depth, dose differences in comparison to the TG-43 method were about -3%. When the source was positioned at the skin surface, dose differences were smaller than -1% for Co-60 and Ir-192, yet -3% for Yb-169. For the interstitial implant, dose differences at the skin surface were -7% for Co-60, -0.6% for Ir-192, and -2.5% for Yb-169. Conclusions: This study indicates the following: (i) for the superficial mould, no bolus is needed; (ii) when the source is in contact with the skin surface, no bolus is needed for either Co-60 and Ir-192. For lower energy radionuclides like Yb-169, bolus may be needed; and (iii) for the interstitial case, at least a 0.1 cm bolus is advised for Co-60 to avoid underdosing superficial target layers. For Ir-192 and Yb-169, no bolus is needed. For those cases where no bolus is needed, its use might be detrimental as the lack of radiation scatter may be beneficial to the patient, although the 2% tolerance for dose calculation accuracy recommended in the AAPM TG-56 report is not fulfilled.
|
|
|
|
Brown, J. M. C., Dimmock, M. R., Gillam, J. E., & Paganin, D. M. (2014). A low energy bound atomic electron Compton scattering model for Geant4. Nucl. Instrum. Methods Phys. Res. B, 338, 77–88.
Abstract: A two-body fully relativistic three-dimensional scattering framework has been utilised to develop an alternative Compton scattering computational model to those adapted from Ribberfors' work for Monte Carlo modelling of Compton scattering. Using a theoretical foundation that ensures the conservation of energy and momentum in the relativistic impulse approximation, this new model, the Monash University Compton scattering model, develops energy and directional algorithms for both the scattered photon and ejected Compton electron from first principles. The Monash University Compton scattering model was developed to address the limitation of the Compton electron directionality algorithms of other computational models adapted from Ribberfors' work. Here the development of the Monash University Compton scattering model, including its implementation in a Geant4 low energy electromagnetic physics class, G4LowEPComptonModel, is outlined. Assessment of the performance of G4LowEPComptonModel was undertaken in two steps: (1) comparison with respect to the two standard Compton scattering classes of Geant4 version 9.5, G4LivermoreComptonModel and G4PenelopeComptonModel, and (2) experimental comparison with respect to Compton electron kinetic energy spectra obtained from the Compton scattering of 662 key photons off the K-shell of gold. Both studies illustrate that the Monash University Compton scattering model, and in turn G4LowEPComptonModel, is a viable replacement for the majority of computational models that have been adapted from Ribberfors' work. It was also shown that the Monash University Compton scattering model is able to reproduce the Compton scattering triply differential cross-section Compton electron kinetic energy spectra of 662 keV photons K-shell scattering off of gold to within experimental uncertainty.
|
|
|
|
Campanario, F., & Kubocz, M. (2014). Higgs boson CP-properties of the gluonic contributions in Higgs plus three jet production via gluon fusion at the LHC. J. High Energy Phys., 10(10), 173–16pp.
Abstract: in high energy hadronic collisions, a general CP-violating Higgs boson Phi with accompanying jets can be efficiently produced via gluon fusion, which is mediated by heavy quark loops. In this article, we study the dominant sub-channel gg -> ggg Phi of the gluon fusion production process with triple real emission corrections at order alpha(5)(s). We go beyond the heavy top-quark approximation and include the full mass dependence of the top- and bottom-quark contributions. Furthermore, in a specific model we demonstrate the features of our program and show the impact of bottom-quark loop contributions in combination with large values of tan beta on differential distributions sensitive to CP-rneasurements of the Higgs boson.
|
|
|
|
Pujades, M. C., Granero, D., Vijande, J., Ballester, F., Perez-Calatayud, J., Papagiannis, P., et al. (2014). Air-kerma evaluation at the maze entrance of HDR brachytherapy facilities. J. Radiol. Prot., 34(4), 741–753.
Abstract: In the absence of procedures for evaluating the design of brachytherapy (BT) facilities for radiation protection purposes, the methodology used for external beam radiotherapy facilities is often adapted. The purpose of this study is to adapt the NCRP 151 methodology for estimating the air-kerma rate at the door in BT facilities. Such methodology was checked against Monte Carlo (MC) techniques using the code Geant4. Five different facility designs were studied for Ir-192 and Co-60 HDR applications to account for several different bunker layouts. For the estimation of the lead thickness needed at the door, the use of transmission data for the real spectra at the door instead of the ones emitted by Ir-192 and Co-60 will reduce the lead thickness by a factor of five for Ir-192 and ten for Co-60. This will significantly lighten the door and hence simplify construction and operating requirements for all bunkers. The adaptation proposed in this study to estimate the air-kerma rate at the door depends on the complexity of the maze: it provides good results for bunkers with a maze (i.e. similar to those used for linacs for which the NCRP 151 methodology was developed) but fails for less conventional designs. For those facilities, a specific Monte Carlo study is in order for reasons of safety and cost-effectiveness.
|
|
|
|
Ortega, P. G., Torres-Espallardo, I., Cerutti, F., Ferrari, A., Gillam, J. E., Lacasta, C., et al. (2015). Noise evaluation of Compton camera imaging for proton therapy. Phys. Med. Biol., 60(5), 1845–1863.
Abstract: Compton Cameras emerged as an alternative for real-time dose monitoring techniques for Particle Therapy (PT), based on the detection of prompt-gammas. As a consequence of the Compton scattering process, the gamma origin point can be restricted onto the surface of a cone (Compton cone). Through image reconstruction techniques, the distribution of the gamma emitters can be estimated, using cone-surfaces backprojections of the Compton cones through the image space, along with more sophisticated statistical methods to improve the image quality. To calculate the Compton cone required for image reconstruction, either two interactions, the last being photoelectric absorption, or three scatter interactions are needed. Because of the high energy of the photons in PT the first option might not be adequate, as the photon is not absorbed in general. However, the second option is less efficient. That is the reason to resort to spectral reconstructions, where the incoming. energy is considered as a variable in the reconstruction inverse problem. Jointly with prompt gamma, secondary neutrons and scattered photons, not strongly correlated with the dose map, can also reach the imaging detector and produce false events. These events deteriorate the image quality. Also, high intensity beams can produce particle accumulation in the camera, which lead to an increase of random coincidences, meaning events which gather measurements from different incoming particles. The noise scenario is expected to be different if double or triple events are used, and consequently, the reconstructed images can be affected differently by spurious data. The aim of the present work is to study the effect of false events in the reconstructed image, evaluating their impact in the determination of the beam particle ranges. A simulation study that includes misidentified events (neutrons and random coincidences) in the final image of a Compton Telescope for PT monitoring is presented. The complete chain of detection, from the beam particle entering a phantom to the event classification, is simulated using FLUKA. The range determination is later estimated from the reconstructed image obtained from a two and three-event algorithm based on Maximum Likelihood Expectation Maximization. The neutron background and random coincidences due to a therapeutic-like time structure are analyzed for mono-energetic proton beams. The time structure of the beam is included in the simulations, which will affect the rate of particles entering the detector.
|
|
|
|
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.
|
|
