|
|
Gimenez-Alventosa, V., Gimenez, V., Ballester, F., Vijande, J., & Andreo, P. (2020). Monte Carlo calculation of beam quality correction factors for PTW cylindrical ionization chambers in photon beams. Phys. Med. Biol., 65(20), 205005–11pp.
Abstract: The beam quality correction factork(Q)for megavoltage photon beams has been calculated for eight PTW (Freiburg, Germany) ionization chambers (Farmer chambers PTW30010, PTW30011, PTW30012, and PTW30013, Semiflex 3D chambers PTW31021, PTW31010, and PTW31013, and the PinPoint 3D chamber PTW31016). Simulations performed on the widely used NE-2571 ionization chamber have been used to benchmark the results. The Monte Carlo code PENELOPE/penEasy was used to calculate the absorbed dose to a point in water and the absorbed dose to the active air volume of the chambers for photon beams in the range 4 to 24 MV. Of the nine ionization chambers analysed, only five are included in the current version of the International Code of Practice for dosimetry based on standards of absorbed dose to water (IAEA TRS 398). The values reported in this work agree with those in the literature within the uncertainty estimates and are to be included in the average values of the data obtained by different working groups for the forthcoming update of TRS 398.
|
|
|
|
Roser, J., Muñoz, E., Barrientos, L., Barrio, J., Bernabeu, J., Borja-Lloret, M., et al. (2020). Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events. Phys. Med. Biol., 65(14), 145005–17pp.
Abstract: Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enablein vivorange verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.
|
|
|
|
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.
|
|
|
|
Muñoz, E., Barrientos, L., Bernabeu, J., Borja-Lloret, M., Llosa, G., Ros, A., et al. (2020). A spectral reconstruction algorithm for two-plane Compton cameras. Phys. Med. Biol., 65(2), 025011–17pp.
Abstract: One factor limiting the current applicability extent of hadron therapy is the lack of a reliable method for real time treatment monitoring. The use of Compton imaging systems as monitors requires the correct reconstruction of the distribution of prompt gamma productions during patient irradiation. In order to extract the maximum information from all the measurable events, we implemented a spectral reconstruction method that assigns to all events a probability of being either partial or total energy depositions. The method, implemented in a list-mode maximum likelihood expectation maximization algorithm, generates a four dimensional image in the joint spatial-spectral domain, in which the voxels containing the emission positions and energies are obtained. The analytical model used for the system response function is also employed to derive an analytical expression for the sensitivity, which is calculated via Monte Carlo integration. The performance of the method is evaluated through reconstruction of various experimental and simulated sources with different spatial and energy distributions. The results show that the proposed method can recover the spectral and spatial information simultaneously, but only under the assumption of ideal measurements. The analysis of the Monte Carlo simulations has led to the identification of two important degradation sources: the mispositioning of the gamma interaction point and the missing energy recorded in the interaction. Both factors are related to the high energy transferred to the recoil electrons, which can travel far from the interaction point and even escape the detector. These effects prevent the direct application of the current method in more realistic scenarios. Nevertheless, experimental point-like sources have been accurately reconstructed and the spatial distributions and spectral emission of complex simulated phantoms can be identified.
|
|
|
|
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.
|
|
|
|
Muñoz, E., Barrio, J., Bernabeu, J., Etxebeste, A., Lacasta, C., Llosa, G., et al. (2018). Study and comparison of different sensitivity models for a two-plane Compton camera. Phys. Med. Biol., 63(13), 135004–19pp.
Abstract: Given the strong variations in the sensitivity of Compton cameras for the detection of events originating from different points in the field of view (FoV), sensitivity correction is often necessary in Compton image reconstruction. Several approaches for the calculation of the sensitivity matrix have been proposed in the literature. While most of these models are easily implemented and can be useful in many cases, they usually assume high angular coverage over the scattered photon, which is not the case for our prototype. In this work, we have derived an analytical model that allows us to calculate a detailed sensitivity matrix, which has been compared to other sensitivity models in the literature. Specifically, the proposed model describes the probability of measuring a useful event in a two-plane Compton camera, including the most relevant physical processes involved. The model has been used to obtain an expression for the system and sensitivity matrices for iterative image reconstruction. These matrices have been validated taking Monte Carlo simulations as a reference. In order to study the impact of the sensitivity, images reconstructed with our sensitivity model and with other models have been compared. Images have been reconstructed from several simulated sources, including point-like sources and extended distributions of activity, and also from experimental data measured with Na-22 sources. Results show that our sensitivity model is the best suited for our prototype. Although other models in the literature perform successfully in many scenarios, they are not applicable in all the geometrical configurations of interest for our system. In general, our model allows to effectively recover the intensity of point-like sources at different positions in the FoV and to reconstruct regions of homogeneous activity with minimal variance. Moreover, it can be employed for all Compton camera configurations, including those with low angular coverage over the scatterer.
|
|
|
|
Gimenez-Alventosa, V., Gimenez, V., Ballester, F., Vijande, J., & Andreo, P. (2018). Correction factors for ionization chamber measurements with the 'Valencia' and 'large field Valencia' brachytherapy applicators. Phys. Med. Biol., 63(12), 125004–10pp.
Abstract: Treatment of small skin lesions using HDR brachytherapy applicators is a widely used technique. The shielded applicators currently available in clinical practice are based on a tungsten-alloy cup that collimates the source-emitted radiation into a small region, hence protecting nearby tissues. The goal of this manuscript is to evaluate the correction factors required for dose measurements with a plane-parallel ionization chamber typically used in clinical brachytherapy for the 'Valencia' and 'large field Valencia' shielded applicators. Monte Carlo simulations have been performed using the PENELOPE-2014 system to determine the absorbed dose deposited in a water phantom and in the chamber active volume with a Type A uncertainty of the order of 0.1%. The average energies of the photon spectra arriving at the surface of the water phantom differ by approximately 10%, being 384 keV for the 'Valencia' and 343 keV for the 'large field Valencia'. The ionization chamber correction factors have been obtained for both applicators using three methods, their values depending on the applicator being considered. Using a depth-independent global chamber perturbation correction factor and no shift of the effective point of measurement yields depth-dose differences of up to 1% for the 'Valencia' applicator. Calculations using a depth-dependent global perturbation factor, or a shift of the effective point of measurement combined with a constant partial perturbation factor, result in differences of about 0.1% for both applicators. The results emphasize the relevance of carrying out detailed Monte Carlo studies for each shielded brachytherapy applicator and ionization chamber.
|
|
|
|
Muñoz, E., Barrio, J., Etxebeste, A., Ortega, P. G., Lacasta, C., Oliver, J. F., et al. (2017). Performance evaluation of MACACO: a multilayer Compton camera. Phys. Med. Biol., 62(18), 7321–7341.
Abstract: Compton imaging devices have been proposed and studied for a wide range of applications. We have developed a Compton camera prototype which can be operated with two or three detector layers based on monolithic lanthanum bromide (LaBr3) crystals coupled to silicon photomultipliers (SiPMs), to be used for proton range verification in hadron therapy. In this work, we present the results obtained with our prototype in laboratory tests with radioactive sources and in simulation studies. Images of a Na-22 and an Y-88 radioactive sources have been successfully reconstructed. The full width half maximum of the reconstructed images is below 4 mm for a Na-22 source at a distance of 5 cm.
|
|
|
|
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).
|
|
|
|
Garcia-Cases, F., Perez-Calatayud, J., Ballester, F., Vijande, J., & Granero, D. (2018). Peripheral dose around a mobile linac for intraoperative radiotherapy: radiation protection aspects. J. Radiol. Prot., 38(4), 1393–1411.
Abstract: The aim of this work is to analyse the scattered radiation produced by the mobile accelerator Mobetron 1000. To do so, detailed Monte Carlo simulations using two different codes, Penelope2008 and Geant4, were performed. Measurements were also done. To quantify the attenuation due to the internal structures, present in the accelerator head, on the scattered radiation produced, some of the main structural shielding in the Mobetron 1000 has been incorporated into the geometry simulation. Results are compared with measurements. Some discrepancies between the calculated and measured dose values were found. These differences can be traced back to the importance of the radiation component due to low energy scattered electrons. This encouraged us to perform additional calculations to separate the role played by this component. Ambient dose equivalent, H*(10), outside of the operating room (OR) has been evaluated using Geant4. H*(10) has been measured inside and outside the OR, being its values compatible with those reported in the literature once the low energy electron component is removed. With respect to the role played by neutrons, estimations of neutron H*(10) using Geant4 together with H*(10) measurements has been performed for the case of the 12 MeV electron beam. The values obtained agree with the experimental values existing in the literature, being much smaller than those registered in conventional accelerators. This study is a useful tool for the clinical user to investigate the radiation protection issues arising with the use of these accelerators in ORs without structural shielding. These results will also enable to better fix the maximum number of treatments that could be performed while insuring adequate radiological protection of workers and public in the hospital.
|
|
