Search

The applicability extent of hadron therapy for tumor treatment is currently limited by the lack of reliable online monitoring techniques. An active topic of investigation is the research of monitoring systems based on the detection of secondary radiation produced during treatment. MACACO, a multi-layer Compton camera based on LaBr3 scintillator crystals and SiPMs, is being developed at IFIC-Valencia for this purpose. This work reports the results obtained from measurements of a 150 MeV proton beam impinging on a PMMA target. A neural network trained on Monte Carlo simulations is used for event selection, increasing the signal to background ratio before image reconstruction. Images of the measured prompt gamma distributions are reconstructed by means of a spectral reconstruction code, through which the 4.439 MeV spectral line is resolved. Images of the emission distribution at this energy are reconstructed, allowing calculation of the distal fall-off and identification of target displacements of 3 mm.

The IRIS group of IFIC-Valencia is developing a Compton telescope for treatment monitoring in hadron therapy. The system is composed of three detector layers, each of them consisting of a monolithic LaBr 3 scintillator crystal coupled to Silicon Photomultiplier (SiPM) arrays. A first version of the prototype was developed within the European project ENVISION. Tests with point-like sources and in-beam were carried out with successful results, demonstrating the feasibility of the proposed technology. The limitations of the system were also identified. A new version of the prototype is under development with the aim of improving the performance of the telescope detector planes. To this end, the initial tests of a detector plane based on a new SiPM array have been carried out. The first characterization, coincidence and imaging tests with this new version of the detector plane are presented. In addition, the suitability of CeBr 3 scintillator crystals for this application is being investigated.

International audience; 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°. Moreover, with the aim of demonstrating the versatility of such a tool the performance of two different Compton camera designs was evaluated and compared.

The SAFIR development represents a novel Positron Emission Tomography (PET) detector, conceived for preclinical fast acquisitions inside the bore of a Magnetic Resonance Imaging (MRI) scanner. The goal is hybrid and simultaneous PET/MRI dynamic studies at unprecedented temporal resolutions of a few seconds. The detector relies on matrices of scintillating LSO-based crystals coupled one-to-one with SiPM arrays and readout by fast ASICs with excellent timing resolution and high rate capabilities. The paper describes the detector concept and the initial results in terms of simulations and characterisation measurements.

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.

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.

SAFIR is a high-rate PET insert for a preclinical 7 T MRI scanner, aiming at studying fast biological processes in small animals, using tracer activities up to 500 MBq. The full detector comprises 15360 channels, each composed of a LYSO crystal and a $S$ iPM photo-sensor. A dedicated readout system, operating inside the MRI magnet of 200 mm inner diameter acquires 40 kHz of hits per channel. PETA6 ASICs digitize amplitude and time of the hits. FPGAs read, process and transfer, via optical links, the digital data. We present the readout system in detail, together with first test results.

A three-layer Compton Telescope, capable of measuring two and three-interaction events, is being developed at IRIS group (Image Reconstruction, Instrumentation and Simulations for medical imaging applications) of IFIC-Valencia. It is based on continuous LaBr 3 crystals coupled to SiPMs; consequently, a large number of projections can be obtained (high granularity). In order to achieve an accurate reconstruction process, precise estimations of the System Matrix and the Sensitivity Matrix are needed. Given that Monte-Carlo calculations of the aforementioned quantities are computationally expensive due to the high granularity of the imaging device, an analytical model would be advisable. The group has already proposed and assessed a model for two-interaction events. In the present work we propose a similar model for three-interaction events. Whereas three-interaction events allow to obtain the energy of incident gamma-rays without full-absorption requirement, efficiency loss and richer diversity of signal/background events must be handled. In this work we present a preliminary validation of the proposed model, by comparing its performance and predictions with Monte-Carlo simulations. Images of simulated point-like monochromatic sources, as well as a preliminary background study, are also reported.

The suitability of monolithic Lutetium Fine Silicate (LFS) continuous scintillation crystals coupled to 64-channel SiPM-TSV matrices for small animal PET has been evaluated. Three crystals of size 25.8 × 25 . 8 × 15 mm 3 were tested: LYSO white coating, LFS white coating (LFS-W) and LFS specular coating (LFS-S). Electronically collimated measurements were carried out to compare the performance of LYSO and LFS crystals. Mean energy resolutions of 14%, 17% and 19% were measured for LYSO, LFS-S and LFS-W respectively. Similar intrinsic spatial resolutions were found for all three crystals in the central region. Different versions of the light model used to estimate the interaction position of the photons in the detector have been studied. Further measurements with two detector heads of each crystal type in coincidence have also been carried out with a moveable 22Na source. A ML-MLEM image reconstruction code is used to reconstruct the data. Final position resolutions for all three configurations show similar values. Taking into account that LYSO and LFS present similar performance levels (and considering that LFS has better timing resolution than LYSO) we conclude that LFS is a valid substitute for LYSO.

Accidental coincidences are one of the main sources of image degradation in positron emission tomography. It is possible to compensate for their degradation e ects, but an accurate method to estimate the randoms rate is required. Two conventional methods are used for random rate estimation: the singles rate method and the delayed window method. In this work we propose a mathematical model that describes the process of accidental coincidence formation. By using it, we are able to predict the correct number of randoms as well as the estimations provided the singles rate and delayed window methods. The model is used to propose a novel estimation method: the singles-prompt . The aim of this work is to assess the performance of the singles-prompt method and, specially, the model capabilities at several levels. The results agree with all the predictions of the model. In particular, the singles rate and delayed window estimations behave as described by the model and the singles-prompt method estimates the correct number of randoms regardless of the source distribution and total activities.

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.

In a dual PET ring detector, a high resolution detector ring is placed within a standard PET ring. For the outer ring detector with a spatial resolution significantly worse than the resolution of the inner ring, it has been shown that the combined detector performs significantly better than a ring only system. To test performance in a system with a high resolution ring detectors, a LYSO PET module was built with total size equal to the standard size of PET detector modules, and a resolution of 1.5 mm FWHM. This paper reports characterization and calibration of the LYSO PET module and results of the preliminary PET imaging test.

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 ([Formula: see text]) 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 [Formula: see text]Na and an [Formula: see text]Y radioactive sources have been successfully reconstructed. The full width half maximum of the reconstructed images is below 4 mm for a [Formula: see text]Na source at a distance of 5 cm.

A high-resolution silicon detector probe, in coincidence with a conventional PET scanner, is expected to provide images of higher quality than those achievable using the scanner alone. Spatial resolution should improve due to the finer pixelization of the probe detector, while increased sensitivity in the probe vicinity is expected to decrease noise. A PET-probe prototype is being developed utilizing this principle. The system includes a probe consisting of ten layers of silicon detectors, each a 80 × 52 array of 1 × 1 × 1 mm(3) pixels, to be operated in coincidence with a modern clinical PET scanner. Detailed simulation studies of this system have been performed to assess the effect of the additional probe information on the quality of the reconstructed images. A grid of point sources was simulated to study the contribution of the probe to the system resolution at different locations over the field of view (FOV). A resolution phantom was used to demonstrate the effect on image resolution for two probe positions. A homogeneous source distribution with hot and cold regions was used to demonstrate that the localized improvement in resolution does not come at the expense of the overall quality of the image. Since the improvement is constrained to an area close to the probe, breast imaging is proposed as a potential application for the novel geometry. In this sense, a simplified breast phantom, adjacent to heart and torso compartments, was simulated and the effect of the probe on lesion detectability, through measurements of the local contrast recovery coefficient-to-noise ratio (CNR), was observed. The list-mode ML-EM algorithm was used for image reconstruction in all cases. As expected, the point spread function of the PET-probe system was found to be non-isotropic and vary with position, offering improvement in specific regions. Increase in resolution, of factors of up to 2, was observed in the region close to the probe. Images of the resolution phantom showed visible improvement in resolution when including the probe in the simulations. The image quality study demonstrated that contrast and spill-over ratio in other areas of the FOV were not sacrificed for this enhancement. The CNR study performed on the breast phantom indicates increased lesion detectability provided by the probe.

The development of novel detection devices and systems such as the AX-positron emission tomography (PET) demonstrator often introduce or increase the measurement of atypical coincidence events such as inter-crystal scattering (ICS). In more standard systems, ICS events often go undetected and the small measured fraction may be ignored. As the measured quantity of such events in the data increases, so too does the importance of considering them during image reconstruction. Generally, treatment of ICS events will attempt to determine which of the possible candidate lines of response (LoRs) correctly determine the annihilation photon trajectory. However, methods of assessment often have low success rates or are computationally demanding. In this investigation alternative approaches are considered. Experimental data was taken using the AX-PET prototype and a NEMA phantom. Three methods of ICS treatment were assessed--each of which considered all possible candidate LoRs during image reconstruction. Maximum likelihood expectation maximization was used in conjunction with both standard (line-like) and novel (V-like in this investigation) detection responses modeled within the system matrix. The investigation assumed that no information other than interaction locations was available to distinguish between candidates, yet the methods assessed all provided means by which such information could be included. In all cases it was shown that the signal to noise ratio is increased using ICS events. However, only one method, which used full modeling of the ICS response in the system matrix--the V-like model--provided enhancement in all figures of merit assessed in this investigation. Finally, the optimal method of ICS incorporation was demonstrated using data from two small animals measured using the AX-PET demonstrator.

AX-PET stands for a new geometrical concept for a high resolution and high sensitivity PET scanner, based on an axial arrangement of long scintillating crystals in the tomograph, for a parallax free PET detector. Two identical AX-PET modules – consisting of matrices of LYSO crystals interleaved with WLS strips – have been built. They form the AX-PET Demonstrator, which has been extensively characterized and successfully used for the reconstruction of images of several phantoms. In this paper we report on the current status of the project, with emphasis on the most relevant results achieved both in terms of detector characterization and image reconstruction. We also discuss the recent preliminary results obtained with the digital SiPM from Philips (dSiPM), which are currently being tested as a possible alternative photodetector for AX-PET. With their very good intrinsic time resolution, dSiPM could add Time of Flight capability to the AX-PET concept.

AX-PET is a novel PET detector based on axially oriented crystals and orthogonal wavelength shifter (WLS) strips, both individually read out by silicon photo-multipliers. Its design decouples sensitivity and spatial resolution, by reducing the parallax error due to the layered arrangement of the crystals. Additionally the granularity of AX-PET enhances the capability to track photons within the detector yielding a large fraction of inter-crystal scatter events. These events, if properly processed, can be included in the reconstruction stage further increasing the sensitivity. Its unique features require dedicated Monte-Carlo simulations, enabling the development of the device, interpreting data and allowing the development of reconstruction codes. At the same time the non-conventional design of AX-PET poses several challenges to the simulation and modeling tasks, mostly related to the light transport and distribution within the crystals and WLS strips, as well as the electronics readout. In this work we present a hybrid simulation tool based on an analytical model and a Monte-Carlo based description of the AX-PET demonstrator. It was extensively validated against experimental data, providing excellent agreement.

The AX-PET collaboration has developed a novel concept for high resolution PET imaging to overcome some of the performance limitations of classical PET cameras, in particular the compromise between spatial resolution and sensitivity introduced by the parallax error. The detector consists of an arrangement of long LYSO scintillating crystals axially oriented around the field of view together with arrays of wave length shifter strips orthogonal to the crystals. This matrix allows a precise 3D measurement of the photon interaction point. This is valid both for photoelectric absorption at 511 keV and for Compton scattering down to deposited energies of about 100 keV. Crystals and WLS strips are individually read out using Geiger-mode Avalanche Photo Diodes (G-APDs). The sensitivity of such a detector can be adjusted by changing the number of layers and the resolution is defined by the crystal and strip dimensions. Two AX-PET modules were built and fully characterized in dedicated test set-ups at CERN, with point-like 22 Na sources. Their performance in terms of energy ( R energy ≈ 11.8 % (FWMH) at 511 keV) and spatial resolution was assessed ( σ axial ≈ 0.65 mm ), both individually and for the two modules in coincidence. Test campaigns at ETH Zurich and at the company AAA allowed the tomographic reconstructions of more complex phantoms validating the 3D reconstruction algorithms. The concept of the AX-PET modules will be presented together with some characterization results. We describe a count rate model which allows to optimize the planing of the tomographic scans.

The probe arrangement in positron emission tomography (PET) is a combination of a high resolution detector (probe) placed within a whole body PET scanner. An additional subset of high-resolution data is recorded that contains events where one of the annihilation photon pair interacted in the high resolution detector. A study was designed to evaluate impact of probe geometry and efficiency on the quality of obtained image. The performance was predicted by estimating mean and variance of efficient estimators such as MLEM based on imaging matrix evaluated at particular probe geometry. Results were compared to data measured at a high-resolution PET probe prototype constructed at University of Michigan. The prototype consists of two arcs of standard PET detector elements at 50 cm radius combined with silicon sensors segmented into 1 mm3 at 5.7 cm radius. A high-resolution Derenzo phantom filled with 18F-flourodeoxyglucose was rotated during a scan. Limited angle probe geometry was simulated by selection of a subset of phantom angular views. Data reconstruction was matched to algorithms evaluated in theoretical predictions. Two parameters, contrast recovery coefficient (CRC) and signal to noise ratio (SNR) were compared to predicted values. The data shows a 50 % increase in CRC and SNR at an event ratio of 10 % where improvement is maximized for a subset of angular positions. Splitting high-resolution data to two perpendicular angular positions improves CRC and SNR.

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 × 16 arrays of 1.4 mm × 1.4 mm pixels and the other in 40 × 26 arrays of 1.0 mm × 1.0 mm 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 × 6 arrays of 6 mm × 12 mm × 30 mm 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.5 mm 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.8 mm and 4.0 mm phantom features respectively, were observed, while previously unresolvable 3.2 mm 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.

The interest in using continuous monolithic crystals in positron emission tomography (PET) has grown in the last years. Coupled to silicon photomultipliers (SiPMs), the detector can combine high sensitivity and high resolution, the two main factors to be maximized in a positron emission tomograph. In this work, the position determination capability of a detector comprised of a [Formula: see text] mm(3) LYSO crystal coupled to an [Formula: see text]-pixel array of SiPMs is evaluated. The 3D interaction position of γ-rays is estimated using an analytical model of the light distribution including reflections on the facets of the crystal. Monte Carlo simulations have been performed to evaluate different crystal reflectors and geometries. The method has been characterized and applied to different cases. Intrinsic resolution obtained with the position estimation method used in this work, applied to experimental data, achieves sub-millimetre resolution values. Average resolution over the detector surface for 5 mm thick crystal is ∼0.9 mm FWHM and ∼1.2 mm FWHM for 10 mm thick crystal. Depth of interaction resolution is close to 2 mm FWHM in both cases, while the FWTM is ∼5.3 mm for 5 mm thick crystal and ∼9.6 mm for 10 mm thick crystal.

A Compton telescope for dose monitoring in hadron therapy is under development at IFIC. The system consists of three layers of LaBr3 crystals coupled to silicon photomultiplier arrays. 22Na sources have been successfully imaged reconstructing the data with an ML-EM code. Calibration and temperature stabilization are necessary for the prototype operation at low coincidence rates. A spatial resolution of 7.8 mm FWHM has been obtained in the first imaging tests., This work was supported in part by the European Commissions 7th Framework Programme through the ENVISION (G. A. num 241851) project and the Marie Curie ITN ENTERVISION (G. A. num 264552), by the Spanish Ministerio de Economia y Competitividad (FPA2010-14891, FIS2011-14585-E, and FPA2014-53599-R), the UVEG (UV-INV-PRECOMP12-80755), and the Generalitat Valenciana (GV/2013/133 and GVISIC/2012/020). Group members are supported by Ramón y Cajal, UVEG Atraccio de Talent, and Generalitat Valenciana contracts.

Random coincidences degrade the image in Positron Emission Tomography, PET. To compensate for their degradation effects, the rate of random coincidences should be estimated. Under certain circumstances, current estimation methods fail to provide accurate results. We propose a novel method, "Singles-Prompts" (SP), that includes the information conveyed by prompt coincidences and models the pile-up. The SP method has the same structure than the well-known "Singles Rate" (SR) approach. Hence, SP can straightforwardly replace SR. In this work, the SP method has been extensively assessed and compared to two conventional methods, SR and the delayed window (DW) method, in a preclinical PET scenario using Monte-Carlo simulations. SP offers accurate estimates for the randoms rates, while SR and DW tend to overestimate the rates (∼10%, and 5%, respectively). With pile-up, the SP method is more robust than SR (but less than DW). At the image level, the contrast is overestimated in SR-corrected images, +16%, while SP produces the correct value. Spill-over is slightly reduced using SP instead of SR. The DW images values are similar to those of SP except for low-statistic scenarios, where DW behaves as if randoms were not compensated for. In particular, the contrast is reduced, -16%. In general, the better estimations of SP translate into better image quality.

Axial PET is a novel geometrical concept for Positron Emission Tomography (PET), based on layers of long scintillating crystals axially aligned with the bore axis. The axial coordinate is obtained from arrays of wavelength shifting (WLS) plastic strips placed orthogonally to the crystals. This article describes the design, construction and performance evaluation of a demonstrator set-up which consists of two identical detector modules, used in coincidence. Each module comprises 48 LYSO crystals of 100 mm length and 156 WLS strips. Crystals and strips are readout by Geiger-mode Avalanche Photo Diodes (G-APDs). The signals from the two modules are processed by fully analog front-end electronics and recorded in coincidence by a VME-based data acquisition system. Measurements with point-like 22 Na sources, with the modules used both individually and in coincidence mode, allowed for a complete performance evaluation up to the focal plane reconstruction of point sources. The results obtained are in good agreement with expectations and proved the set-up to be ready for the next evaluation phase with PET phantoms filled with radiotracers.

The design and construction of a PET camera module with high sensitivity, full 3-D spatial reconstruction and very good energy resolution is presented. The basic principle consists of an axial arrangement of long scintillation crystals around the Field Of View (FOV), providing a measurement of the transverse coordinates of the interacting 511 keV gamma ray. On top of each layer of crystals, an array of Wave-Length Shifter (WLS) strips, which collect the light leaving the crystals sideways, is positioned orthogonal to the crystal direction. The signals in the WLS strips allow a precise measurement of the z (axial) co-ordinate of the 511 keV γ-ray gamma impact. The construction of two modules used for demonstration of the concept is described. First preliminary results on spatial and energy resolution from one full module will be shown.

Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.

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 22Na 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.

Hadron therapy offers the possibility of delivering a large amount of radiation dose to tumors with minimal absorption by the surrounding healthy tissue. In order to fully exploit the advantages of this technique, the use of real-time beam monitoring devices becomes mandatory. Compton imaging devices can be employed to map the distribution of prompt gamma emission during the treatment and thus assess its correct delivery. The Compton telescope prototype developed at IFIC-Valencia for this purpose is made of three layers of LaBr3 crystals coupled to silicon photomultipliers. The system has been tested in a 4.44 MeV gamma field at the 3 MV Tandetron accelerator at HZDR, Dresden. Images of the target with the system in three different positions separated by 10 mm were successfully reconstructed. This indicates the ability of MACACO for imaging the prompt gamma rays emitted at such energies.

We describe a PET device based on a novel method to extract the coordinates of the interaction point of the 511keV γ rays from 100 mm long and thin LYSO (Lutetium Yttrium OxyorthoSilicate) scintillator bars, positioned axially in the tomograph. The coordinate along the hit crystal is measured by using a hodoscope of Wave Length Shifting (WLS) plastic strips mounted perpendicularly to each plane of scintillators. As photodetectors, new Geiger mode Avalanche PhotoDetectors (G-APDs) with integrated electronics are being used to detect both the hit crystal in a block (x and y coordinates) and the interaction point in the crystal (z coordinate) through the light escaping from the crystal and transmitted to the WLS strips. In this way, the γ interaction point can be determined with a spatial resolution of few cubic millimeters down to a minimum deposited energy of about 50 keV, resulting in a volumetric precision very close to the limits imposed by the physics of the positron annihilation. The method allows to increase the detection efficiency without affecting the spatial resolution by adding scintillator planes in the radial direction. A demonstrator scanner, based on two matrices of 8 × 6 LYS crystals and 312 WLS strips, slotted in between the crystals, is under construction. Preliminary results from the feasibility studies of the various components will be presented.

The Small Animal Fast Insert for mRi (SAFIR) will be a PET insert for the Bruker BioSpin 70/30. It aims at applications where fast processes such as blood perfusion in the rodent brain are to be monitored comprehensively and non-invasively. Employing electronics originally designed for time-of-flight applications, coincidence resolving times of less than a nanosecond can be achieved allowing for short coincidence time windows. Consequently, random contributions to the coincidence events are suppressed, making short acquisition frames with very high tracer concentrations possible. This will allow collecting sufficient count statistics in a few seconds. Geant4-based Monte Carlo simulations were used to characterize the performance of the reference design consisting of polished LSO-like crystals, one-to-one coupled to Silicon Photomultipliers. The crystals are grouped into 8×8 matrices, which are arranged into 24 modules and 10 rings. Similar methods as described in the NEMA NU 4-2008 standard were employed. The simulation results on NECR, sensitivity, and spatial resolution will be presented. Most notably, the NECR at 500 MBq is almost nine times higher than at 50 MBq for the given scanner. Combined with a very high sensitivity, this allows for short acquisition times using these very high injected doses. In addition, essentially random-free measurements at standard activities below 50 MBq are possible.

A three layer Compton telescope for dose monitoring in hadron therapy has been developed at IFIC. It consists of three detector layers and each layer is made of a continuous LaBr3 crystal coupled to four silicon photomultiplier (SiPMs) arrays. A dedicated image reconstruction method has also been developed. The telescope has been characterized and tested with radioactive sources and in-beam. Simulations with GATE have been carried out in order to optimize the detector geometry and response. Ongoing work is mainly focused on performance optimization and on the simultaneous acquisition and processing of two- and three-layer interaction events.

Continuous scintillator crystals can offer several advantages over pixellated crystals. Apart from an increase of the active volume, spatial resolution is not limited to the crystal size and Depth of interaction (DoI) can be estimated from the light distribution. Continuous crystals coupled to finely segmented photodetectors can improve simultaneously spatial resolution and sensitivity, the main two factors to maximize in Positron Emission Tomography (PET). The main drawback is the difficulty in the position estimation within the crystal; thus, accurate position estimation algorithms are needed. However, if the standard cuboid crystals are used in the full ring geometry, there will be large wedge-shaped gaps between detector heads. For maximizing the sensitivity and avoiding data truncation, a scanner based on tapered continuous crystals is considered. In this work, a basic comparison between two scanners of the same dimensions based on tapered and cuboid crystals was performed. The study is focused on the quantification of the increase on sensitivity while first preliminary results were obtained in terms of image quality and spatial resolution, An increase of 14.08% on the average sensitivity was obtained for tapered crystals. Image reconstruction was implemented with Filtered Back-Projection (FBP) method using the ideal position interaction of the annihilation photon within the crystal given by the simulations. An improvement on the recovery of homogeneous regions was observed due to less data truncation using tapered crystal.

Hadron-therapy (HT) delivers high dose to the tumor, significantly sparing surrounding healthy tissues. Positron Emission Tomography (PET) can be used to image short-lived positron-emitters, produced by the hadron beam in the patient. PET is nowadays the most assessed imaging technique to monitor HT treatments. The number of positron-emitters produced is much lower than in conventional clinical scenarios, dedicated scanners and acquisition protocols are then required. Sensitivity becomes a crucial parameter since the duration of the acquisition is limited by the beam delivery during treatment (in-beam PET) as well as by the physiological wash-out of the radioisotopes after treatment (in-room PET). AX-PET, a novel PET concept with axially oriented crystals, provides a higher sensitivity almost fully decoupled from its 3D spatial resolution. The flexible design can overcome many of the geometrical constraints imposed by HT monitoring. This work presents a preliminary feasibility study of the application of AX-PET to the HT scenario in in-room modality.

A three layer Compton telescope for dose monitoring in hadron therapy is under development at IFIC - Valencia. It will consist of three detector layers and each layer will be made of a continuous LaBr 3 crystal coupled to four Silicon Photomultiplier (SiPM) arrays. A prototype made of two detector layers has already been assembled and tested. The detector layers consist of a 32×36×5 mm3 and 32×36×10 mm3 LaBr 3 crystal, each coupled to four SiPM arrays. The front-end electronics of each detector is based on the VATA64HDR16 application specific integrated circuit (ASIC). The gamma rays coming from a 22Na source have been detected with the two detector layers working in time coincidence, allowing energy and position estimation.

Hadron-therapy (HT) aims to treat tumors by maximizing the dose released to the target and sparing the dose to normal tissues. For a successful outcome it is very important to determine where the maximum dose is deposited; therefore range verification is necessary for treatment optimization and patient safety. Secondary positron emitting isotopes and prompt gamma radiation are produced after the hadron beam passage. This secondary radiation coming from tissue activation could be used for quality control of treatment, provided that it can be detected and employed to reconstruct the beam path using imaging techniques. This is one of the main goals of the ENVISION project: to evaluate and develop on-line monitoring devices for HT, like PET for detecting the annihilation photons from positrons and Compton Cameras (CCs) for prompt gamma radiation detection. In both technologies high sensitivity is required to increase the signal-to-noise ratio of the reconstructed image for a given therapeutic dose. This simulation study focuses on the sensitivity optimization of such devices, PET and CC, taking into account its use and constraints for on-line HT monitoring. Several configurations of both technologies have been investigated using sources generated from hadron beams. In the case of PET data, the time-of-flight (TOF) information has been included too. For individual hadron beams acquired after 5 minutes, differences in range of 3 mm are detected for all the PET configurations, except for the partial-ring of 60 cm diameter. In addition, patient data from a carbon ion treatment at GSI have been simulated and reconstructed. The system with higher sensitivity and angular sampling recovers more accurately areas with no activity (nasal cavity). In the case of the CC data, the quality of the reconstructed image when using 2-interaction events is notably improved when the detector layers are placed covering a larger solid angle.

A prototype of a high resolution small animal PET scanner with continuous LYSO crystals coupled to SiPM matrices has been developed. Images of point-like sources in different positions of the FOV have been obtained with a FWHM better than 1 mm. In the near future, a full-ring PET will be developed. Due to the increase in the number of detector heads, a new Data Acquisition (DAQ) system is needed. First tests with the VATA64HDR16 ASIC from Gamma Medica - IDEAS have been carried out. Also, the use of tapered crystals is being considered for reduction of the gaps between detectors and therefore increase the full-ring sensitivity.

Real-time monitoring techniques are receiving increased attention due to their great potential for improving the precision of treatment delivery for hadron-therapy (HT). The study of prompt-γ(PG) focuses on the energetic photon radiation exiting the patient few nanoseconds after the beam irradiation. Its high intensity compared to other similar technologies like in-beam Positron Emission Tomography (PET) and good correlation with the original dose deposition profile are promising features. On the other hand, the use of neutral particles to monitor the dose comes together with a high noise background, mainly due to neutrons and scattered γ, which are uncorrelated with the original dose map. The search of a robust detector that fully exploits the information from PG is still a matter of research.

The usage of long, axially oriented scintillator crystals in a PET scanner has been shown by the AX-PET Demonstrator as a possible solution for a high resolution and high sensitivity PET detector. In the AX-PET implementation, arrays of wavelength shifting (WLS) strips, placed orthogonally behind every crystal layer, are used to define the axial coordinate. After extensive characterization measurements, the AX-PET Demonstrator has been successfully used for the reconstruction of several phantoms and a few rodents. Possible extensions of the AX-PET concept towards Time Of Flight capabilities have been investigated, using Philips digital SiPMs as alternative photodetector. Promising CRT values equal to 211 ps have been demonstrated, using long crystals with dual sided readout and mean timing method. Finally, we report about an alternative way to reconstruct the axial coordinate: exploiting the light sharing of 100 mm long crystals with special surface treatment resulted in axial resolutions of the order of 4 mm FWHM.

A Compton telescope for dose delivery monitoring in hadron therapy is under development at IFIC Valencia within the European project ENVISION. The telescope will consist of three detector planes, each one composed of a LaBr3 continuous scintillator crystal coupled to four silicon photomultiplier (SiPM) arrays. After the development of a first prototype which served to assess the principle, a second prototype with larger crystals has been assembled and is being tested. The current version of the prototype consists of two detector layers, each one composed of a 32.5 × 35 mm2 crystal coupled to four SiPM arrays. The VATA64HDR16 ASIC has been employed as front-end electronics. The readout system consists of a custom made data acquisition board. Tests with point-like sources have been carried out in the laboratory, assessing the correct functioning of the device. The system optimization is ongoing.

In the development of prototype systems for positron emission tomography a valid and robust image reconstruction algorithm is required. However, prototypes often employ novel detector and system geometries which may change rapidly under optimization. In addition, developing systems generally produce highly granular, or possibly continuous detection domains which require some level of on-the-fly calculation for retention of measurement precision. In this investigation a new method of on-the-fly system matrix calculation is proposed that provides advantages in application to such list-mode systems in terms of flexibility in system modeling. The new method is easily adaptable to complicated system geometries and available computational resources. Detection uncertainty models are used as random number generators to produce ensembles of possible photon trajectories at image reconstruction time for each datum in the measurement list. However, the result of this approach is that the system matrix elements change at each iteration in a non-repetitive manner. The resulting algorithm is considered the simulation of a one-pass list (SOPL) which is generated and the list traversed during image reconstruction. SOPL alters the system matrix in use at each iteration and so behavior within the maximum likelihood-expectation maximization algorithm was investigated. A two-pixel system and a small two dimensional imaging model are used to illustrate the process and quantify aspects of the algorithm. The two-dimensional imaging system showed that, while incurring a penalty in image resolution, in comparison to a non-random equal-computation counterpart, SOPL provides much enhanced noise properties. In addition, enhancement in system matrix quality is straightforward (by increasing the number of samples in the ensemble) so that the resolution penalty can be recovered when desired while retaining improvement in noise properties. Finally the approach is tested and validated against a standard (highly accurate) system matrix using experimental data from a prototype system--the AX-PET.

Positron Emission Tomography (PET) is based on the detection in coincidence of the two photons created in a positron annihilation. In conventional PET, this coincidence identification is usually carried out through a coincidence electronic unit. An accidental coincidence occurs when two photons arising from different annihilations are classified as a coincidence. Accidental coincidences are one of the main sources of image degradation in PET. Some novel systems allow coincidences to be selected post-acquisition in software, or in real time through a digital coincidence engine in an FPGA. These approaches provide the user with extra flexibility in the sorting process and allow the application of alternative coincidence sorting procedures. In this work a novel sorting procedure based on Artificial Neural Network (ANN) techniques has been developed. It has been compared to a conventional coincidence sorting algorithm based on a time coincidence window. The data have been obtained from Monte-Carlo simulations. A small animal PET scanner has been implemented to this end. The efficiency (the ratio of correct identifications) can be selected for both methods. In one case by changing the actual value of the coincidence window used, and in the other by changing a threshold at the output of the neural network. At matched efficiencies, the ANN-based method always produces a sorted output with a smaller random fraction. In addition, two differential trends are found: the conventional method presents a maximum achievable efficiency, while the ANN-based method is able to increase the efficiency up to unity, the ideal value, at the cost of increasing the random fraction. Images reconstructed using ANN sorted data (no compensation for randoms) present better contrast, and those image features which are more affected by randoms are enhanced. For the image quality phantom used in the paper, the ANN method decreases the spillover ratio by a factor of 18%., This work was supported by RD12/0036/0020 Red Tematica de Investigacion Coperativa en Cancer (RTICC) 2012; Ministerio de Ciencia e Innovacion through a JAE-Doc contract associated to CSIC and FSE co-financed; also under projects TEC-2007-61047 and FPA-2010-14891.