| Home | << 1 2 3 4 5 6 >> |
|
Borja-Lloret, M., Barrientos, L., Bernabeu, J., Lacasta, C., Muñoz, E., Ros, A., et al. (2023). Influence of the background in Compton camera images for proton therapy treatment monitoring. Phys. Med. Biol., 68(14), 144001–16pp.
Abstract: Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system. Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities. Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns. Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.
|
|
Borys, D. et al, & Brzezinski, K. (2022). ProTheRaMon-a GATE simulation framework for proton therapy range monitoring using PET imaging. Phys. Med. Biol., 67(22), 224002–15pp.
Abstract: Objective. This paper reports on the implementation and shows examples of the use of the ProTheRaMon framework for simulating the delivery of proton therapy treatment plans and range monitoring using positron emission tomography (PET). ProTheRaMon offers complete processing of proton therapy treatment plans, patient CT geometries, and intra-treatment PET imaging, taking into account therapy and imaging coordinate systems and activity decay during the PET imaging protocol specific to a given proton therapy facility. We present the ProTheRaMon framework and illustrate its potential use case and data processing steps for a patient treated at the Cyclotron Centre Bronowice (CCB) proton therapy center in Krakow, Poland. Approach. The ProTheRaMon framework is based on GATE Monte Carlo software, the CASToR reconstruction package and in-house developed Python and bash scripts. The framework consists of five separated simulation and data processing steps, that can be further optimized according to the user's needs and specific settings of a given proton therapy facility and PET scanner design. Main results. ProTheRaMon is presented using example data from a patient treated at CCB and the J-PET scanner to demonstrate the application of the framework for proton therapy range monitoring. The output of each simulation and data processing stage is described and visualized. Significance. We demonstrate that the ProTheRaMon simulation platform is a high-performance tool, capable of running on a computational cluster and suitable for multi-parameter studies, with databases consisting of large number of patients, as well as different PET scanner geometries and settings for range monitoring in a clinical environment. Due to its modular structure, the ProTheRaMon framework can be adjusted for different proton therapy centers and/or different PET detector geometries. It is available to the community via github (Borys et al 2022).
Keywords: proton therapy; GATE; Monte Carlo simulations; J-PET; medical imaging
|
|
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.
|
|
Brzezinski, K. et al. (2023). Detection of range shifts in proton beam therapy using the J-PET scanner: a patient simulation study. Phys. Med. Biol., 68(14), 145016–17pp.
Abstract: Objective. The Jagiellonian positron emission tomography (J-PET) technology, based on plastic scintillators, has been proposed as a cost effective tool for detecting range deviations during proton therapy. This study investigates the feasibility of using J-PET for range monitoring by means of a detailed Monte Carlo simulation study of 95 patients who underwent proton therapy at the Cyclotron Centre Bronowice (CCB) in Krakow, Poland. Approach. Discrepancies between prescribed and delivered treatments were artificially introduced in the simulations by means of shifts in patient positioning and in the Hounsfield unit to the relative proton stopping power calibration curve. A dual-layer, cylindrical J-PET geometry was simulated in an in-room monitoring scenario and a triple-layer, dual-head geometry in an in-beam protocol. The distribution of range shifts in reconstructed PET activity was visualized in the beam's eye view. Linear prediction models were constructed from all patients in the cohort, using the mean shift in reconstructed PET activity as a predictor of the mean proton range deviation. Main results. Maps of deviations in the range of reconstructed PET distributions showed agreement with those of deviations in dose range in most patients. The linear prediction model showed a good fit, with coefficient of determination r (2) = 0.84 (in-room) and 0.75 (in-beam). Residual standard error was below 1 mm: 0.33 mm (in-room) and 0.23 mm (in-beam). Significance. The precision of the proposed prediction models shows the sensitivity of the proposed J-PET scanners to shifts in proton range for a wide range of clinical treatment plans. Furthermore, it motivates the use of such models as a tool for predicting proton range deviations and opens up new prospects for investigations into the use of intra-treatment PET images for predicting clinical metrics that aid in the assessment of the quality of delivered treatment.
|
|
Calatayud-Jordan, J., Candela-Juan, C., Palma, J. D., Pujades-Claumarchirant, M. C., Soriano, A., Gracia-Ochoa, M., et al. (2021). Influence of the simultaneous calibration of multiple ring dosimeters on the individual absorbed dose. J. Radiol. Prot., 41(2), 384–397.
Abstract: Ring dosimeters for personal dosimetry are calibrated in accredited laboratories following ISO 4037-3 guidelines. The simultaneous irradiation of multiple dosimeters would save time, but has to be carefully studied, since the scattering conditions could change and influence the absorbed dose in nearby dosimeters. Monte Carlo simulations using PENELOPE-2014 were performed to explore the need to increase the uncertainty of H-p (0.07) in the simultaneous irradiation of three and five DXT-RAD 707H-2 (Thermo Scientific) ring dosimeters with beam qualities: N-30, N-80 and N-300. Results show that the absorbed dose in each dosimeter is compatible with each of the others and with the reference simulation (a single dosimeter), with a coverage probability of 95% (k = 2). Comparison with experimental data yielded consistent results with the same coverage probability. Therefore, five ring dosimeters can be simultaneously irradiated with beam qualities ranging, at least, between N-30 and N-300 with a negligible impact on the uncertainty of H-p (0.07).
Keywords: ring dosimeters; personal dosimetry; calibration; Monte Carlo; ISO 4037
|
|
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.
Keywords: QCD Phenomenology; Monte Carlo Simulations
|
|
Candela-Juan, C., Niatsetski, Y., van der Laarse, R., Granero, D., Ballester, F., Perez-Calatayud, J., et al. (2016). Design and characterization of a new high-dose-rate brachytherapy Valencia applicator for larger skin lesions. Med. Phys., 43(4), 1639–1648.
Abstract: Purpose: The aims of this study were (i) to design a new high-dose-rate (HDR) brachytherapy applicator for treating surface lesions with planning target volumes larger than 3 cm in diameter and up to 5 cm in size, using the microSelectron-HDR or Flexitron afterloader (Elekta Brachytherapy) with a Ir-192 source; (ii) to calculate by means of the Monte Carlo (MC) method the dose distribution for the new applicator when it is placed against a water phantom; and (iii) to validate experimentally the dose distributions in water. Methods: The PENELOPE2008 MC code was used to optimize dwell positions and dwell times. Next, the dose distribution in a water phantom and the leakage dose distribution around the applicator were calculated. Finally, MC data were validated experimentally for a 192Ir mHDR-v2 source by measuring (i) dose distributions with radiochromic EBT3 films (ISP); (ii) percentage depth-dose (PDD) curve with the parallel-plate ionization chamber Advanced Markus (PTW); and (iii) absolute dose rate with EBT3 films and the PinPoint T31016 (PTW) ionization chamber. Results: The new applicator is made of tungsten alloy (Densimet) and consists of a set of interchangeable collimators. Three catheters are used to allocate the source at prefixed dwell positions with preset weights to produce a homogenous dose distribution at the typical prescription depth of 3 mm in water. The same plan is used for all available collimators. PDD, absolute dose rate per unit of air kerma strength, and off-axis profiles in a cylindrical water phantom are reported. These data can be used for treatment planning. Leakage around the applicator was also scored. The dose distributions, PDD, and absolute dose rate calculated agree within experimental uncertainties with the doses measured: differences of MC data with chamber measurements are up to 0.8% and with radiochromic films are up to 3.5%. Conclusions: The new applicator and the dosimetric data provided here will be a valuable tool in clinical practice, making treatment of large skin lesions simpler, faster, and safer. Also the dose to surrounding healthy tissues is minimal.
Keywords: skin applicator; Valencia applicator; HDR brachytherapy; dosimetry; Monte Carlo
|
|
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.
Keywords: gamma-Ray spectroscopy; Tracking; Monte Carlo
|
|
Etxebeste, A., Dauvergne, D., Fontana, M., Letang, J. M., Llosa, G., Muñoz, E., et al. (2020). CCMod: a GATE module for Compton camera imaging simulation. Phys. Med. Biol., 65(5), 055004–17pp.
Abstract: Compton cameras are gamma-ray imaging systems which have been proposed for a wide variety of applications such as medical imaging, nuclear decommissioning or homeland security. In the design and optimization of such a system Monte Carlo simulations play an essential role. In this work, we propose a generic module to perform Monte Carlo simulations and analyses of Compton Camera imaging which is included in the open-source GATE/Geant4 platform. Several digitization stages have been implemented within the module to mimic the performance of the most commonly employed detectors (e.g. monolithic blocks, pixelated scintillator crystals, strip detectors...). Time coincidence sorter and sequence coincidence reconstruction are also available in order to aim at providing modules to facilitate the comparison and reproduction of the data taken with different prototypes. All processing steps may be performed during the simulation (on-the-fly mode) or as a post-process of the output files (offline mode). The predictions of the module have been compared with experimental data in terms of energy spectra, angular resolution, efficiency and back-projection image reconstruction. Consistent results within a 3-sigma interval were obtained for the energy spectra except for low energies where small differences arise. The angular resolution measure for incident photons of 1275 keV was also in good agreement between both data sets with a value close to 13 degrees. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.
Keywords: Monte Carlo; simulation; gamma imaging; Compton camera
|
|
Etxebeste, A., Barrio, J., Bernabeu, J., Lacasta, C., Llosa, G., Muñoz, E., et al. (2019). Study of sensitivity and resolution for full ring PET prototypes based on continuous crystals and analytical modeling of the light distribution. Phys. Med. Biol., 64(3), 035015–17pp.
Abstract: Sensitivity and spatial resolution are the main parameters to maximize in the performance of a PET scanner. For this purpose, detectors consisting of a combination of continuous crystals optically coupled to segmented photodetectors have been employed. With the use of continuous crystals the sensitivity is increased with respect to the pixelated crystals. In addition, spatial resolution is no longer limited to the crystal size. The main drawback is the difficulty in determining the interaction position. In this work, we present the characterization of the performance of a full ring based on cuboid continuous crystals coupled to SiPMs. To this end, we have employed the simulations developed in a previous work for our experimental detector head. Sensitivity could be further enhanced by using tapered crystals. This enhancement is obtained by increasing the solid angle coverage, reducing the wedge-shaped gaps between contiguous detectors. The performance of the scanners based on both crystal geometries was characterized following NEMA NU 4-2008 standardized protocol in order to compare them. An average sensitivity gain over the entire axial field of view of 13.63% has been obtained with tapered geometry while similar performance of the spatial resolution has been proven with both scanners. The activity at which NECR and true peak occur is smaller and the peak value is greater for tapered crystals than for cuboid crystals. Moreover, a higher degree of homogeneity was obtained in the sensitivity map due to the tighter packing of the crystals, which reduces the gaps and results in a better recovery of homogeneous regions than for the cuboid configuration. Some of the results obtained, such as spatial resolution, depend on the interaction position estimation and may vary if other method is employed.
|