|
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.
|
|
|
NEXT Collaboration(Simon, A. et al), Gomez-Cadenas, J. J., Alvarez, V., Benlloch-Rodriguez, J. M., Botas, A., Carcel, S., et al. (2017). Application and performance of an ML-EM algorithm in NEXT. J. Instrum., 12, P08009–22pp.
Abstract: The goal of the NEXT experiment is the observation of neutrinoless double beta decay in Xe-136 using a gaseous xenon TPC with electroluminescent amplification and specialized photodetector arrays for calorimetry and tracking. The NEXT Collaboration is exploring a number of reconstruction algorithms to exploit the full potential of the detector. This paper describes one of them: the Maximum Likelihood Expectation Maximization (ML-EM) method, a generic iterative algorithm to find maximum-likelihood estimates of parameters that has been applied to solve many different types of complex inverse problems. In particular, we discuss a bi-dimensional version of the method in which the photosensor signals integrated over time are used to reconstruct a transverse projection of the event. First results show that, when applied to detector simulation data, the algorithm achieves nearly optimal energy resolution (better than 0.5% FWHM at the Q value of 136Xe) for events distributed over the full active volume of the TPC.
|
|
|
Grkovski, M., Brzezinski, K., Cindro, V., Clinthorne, N. H., Kagan, H., Lacasta, C., et al. (2015). Evaluation of a high resolution silicon PET insert module. Nucl. Instrum. Methods Phys. Res. A, 788, 86–94.
Abstract: Conventional PET systems can be augmented with additional detectors placed in close proximity of the region of interest. We developed a high resolution PET insert module to evaluate the added benefit of such a combination. The insert module consists of two back-to-back 1 mm thick silicon sensors, each segmented into 1040 1 mm(2) pads arranged in a 40 by 26 array. A set of 16 VATAGP7.1 ASICs and a custom assembled data acquisition board were used to read out the signal from the insert module. Data were acquired in slice (20) geometry with a Jaszczak phantom (rod diameters of 12-4.8 mm) Filled with F-18-FDG and the images were reconstructed with ML-EM method. Both data with full and limited angular coverage from the insert module were considered and three types of coincidence events were combined. The ratio of high-resolution data that substantially improves quality of the reconstructed image for the region near the surface of the insert module was estimated to be about 4%. Results from our previous studies suggest that such ratio could be achieved at a moderate technological expense by using an equivalent of two insert modules (an effective sensor thickness of 4 mm).
|
|
|
Bouhova-Thacker, E., Kostyukhin, V., Koffas, T., Liebig, W., Limper, M., Piacquadio, G. N., et al. (2010). Expected Performance of Vertex Reconstruction in the ATLAS Experiment at the LHC. IEEE Trans. Nucl. Sci., 57(2), 760–767.
Abstract: In the harsh environment of the Large Hadron Collider at CERN (design luminosity of 10(34) cm(-2) s(-1)) efficient reconstruction of vertices is crucial for many physics analyses. Described in this paper is the expected performance of the vertex reconstruction used in the ATLAS experiment. The algorithms for the reconstruction of primary and secondary vertices as well as for finding photon conversions and vertex reconstruction in jets are described. The implementation of vertex algorithms which follows a very modular design based on object-oriented C++ is presented. A user-friendly concept allows event reconstruction and physics analyses to compare and optimize their choice among different vertex reconstruction strategies. The performance of implemented algorithms has been studied on a variety of Monte Carlo samples and results are presented.
|
|
|
Brzezinski, K., Oliver, J. F., Gillam, J., Rafecas, M., Studen, A., Grkovski, M., et al. (2016). Experimental evaluation of the resolution improvement provided by a silicon PET probe. J. Instrum., 11, P09016–13pp.
Abstract: A high-resolution PET system, which incorporates a silicon detector probe into a conventional PET scanner, has been proposed to obtain increased image quality in a limited region of interest. Detailed simulation studies have previously shown that the additional probe information improves the spatial resolution of the reconstructed image and increases lesion detectability, with no cost to other image quality measures. The current study expands on the previous work by using a laboratory prototype of the silicon PET-probe system to examine the resolution improvement in an experimental setting. Two different versions of the probe prototype were assessed, both consisting of a back-to-back pair of 1-mm thick silicon pad detectors, one arranged in 32 x 16 arrays of 1.4mm x 1.4mm pixels and the other in 40 x 26 arrays of 1.0mm x 1.0mm pixels. Each detector was read out by a set of VATAGP7 ASICs and a custom-designed data acquisition board which allowed trigger and data interfacing with the PET scanner, itself consisting of BGO block detectors segmented into 8 x 6 arrays of 6mm x 12mm x 30mm crystals. Limited-angle probe data was acquired from a group of Na-22 point-like sources in order to observe the maximum resolution achievable using the probe system. Data from a Derenzo-like resolution phantom was acquired, then scaled to obtain similar statistical quality as that of previous simulation studies. In this case, images were reconstructed using measurements of the PET ring alone and with the inclusion of the probe data. Images of the Na-22 source demonstrated a resolution of 1.5mm FWHM in the probe data, the PET ring resolution being approximately 6 mm. Profiles taken through the image of the Derenzo-like phantom showed a clear increase in spatial resolution. Improvements in peak-to-valley ratios of 50% and 38%, in the 4.8mm and 4.0mm phantom features respectively, were observed, while previously unresolvable 3.2mm features were brought to light by the addition of the probe. These results support the possibility of improving the image resolution of a clinical PET scanner using the silicon PET-probe.
|
|