Framework for evaluating photon-counting detectors under pile-up conditions.

Purpose: While X-ray photon-counting detectors (PCDs) promise to revolutionize medical imaging, theoretical frameworks to evaluate them are commonly limited to incident fluence rates sufficiently low that the detector response can be considered linear. However, typical clinical operating conditions lead to a significant level of pile-up, invalidating this assumption of a linear response. Here, we present a framework that aims to evaluate PCDs, taking into account their non-linear behavior. Approach: We employ small-signal analysis to study the behavior of PCDs under pile-up conditions. The response is approximated as linear around a given operating point, determined by the incident spectrum and fluence rate. The detector response is subsequently described by the proposed perturbation point spread function (pPSF). We demonstrate this approach using Monte-Carlo simulations of idealized direct- and indirect-conversion PCDs. Results: The pPSFs of two PCDs are calculated. It is then shown how the pPSF allows to determine the sensitivity of the detector signal to an arbitrary lesion. This example illustrates the detrimental influence of pile-up, which may cause non-intuitive effects such as contrast/contrast-to-noise ratio inversion or cancellation between/within energy bins. Conclusions: The proposed framework permits quantifying the spectral and spatial performance of PCDs under clinically realistic conditions at a given operating point. The presented example illustrates why PCDs should not be analyzed assuming that they are linear systems. The framework can, for example, be used to guide the development of PCDs and PCD-based systems. Furthermore, it can be applied to adapt commonly used measures, such as the modulation transfer function, to non-linear PCDs.

Design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomic quantification in hybrid imaging: A feasibility study.

BACKGROUND: Accuracy and precision assessment in radiomic features is important for the determination of their potential to characterize cancer lesions. In this regard, simulation of different imaging conditions using specialized phantoms is increasingly being investigated. In this study, the design and evaluation of a modular multimodality imaging phantom to simulate heterogeneous uptake and enhancement patterns for radiomics quantification in hybrid imaging is presented. METHODS: A modular multimodality imaging phantom was constructed that could simulate different patterns of heterogeneous uptake and enhancement patterns in positron emission tomography (PET), single-photon emission computed tomography (SPECT), computed tomography (CT), and magnetic resonance (MR) imaging. The phantom was designed to be used as an insert in the standard NEMA-NU2 IEC body phantom casing. The entire phantom insert is composed of three segments, each containing three separately fillable compartments. The fillable compartments between segments had different sizes in order to simulate heterogeneous patterns at different spatial scales. The compartments were separately filled with different ratios of 99m Tc-pertechnetate, 18 F-fluorodeoxyglucose ([18 F]FDG), iodine- and gadolinium-based contrast agents for SPECT, PET, CT, and T1 -weighted MR imaging respectively. Image acquisition was performed using standard oncological protocols on all modalities and repeated five times for repeatability assessment. A total of 93 radiomic features were calculated. Variability was assessed by determining the coefficient of quartile variation (CQV) of the features. Comparison of feature repeatability at different modalities and spatial scales was performed using Kruskal-Wallis-, Mann-Whitney U-, one-way ANOVA- and independent t-tests. RESULTS: Heterogeneous uptake and enhancement could be simulated on all four imaging modalities. Radiomic features in SPECT were significantly less stable than in all other modalities. Features in PET were significantly less stable than in MR and CT. A total of 20 features, particularly in the gray-level co-occurrence matrix (GLCM) and gray-level run-length matrix (GLRLM) class, were found to be relatively stable in all four modalities for all three spatial scales of heterogeneous patterns (with CQV < 10%). CONCLUSION: The phantom was suitable for simulating heterogeneous uptake and enhancement patterns in [18 F]FDG-PET, 99m Tc-SPECT, CT, and T1 -weighted MR images. The results of this work indicate that the phantom might be useful for the further development and optimization of imaging protocols for radiomic quantification in hybrid imaging modalities.

Time of Flight in Perspective: Instrumental and Computational Aspects of Time Resolution in Positron Emission Tomography.

The first time-of-flight positron emission tomography (TOF-PET) scanners were developed as early as in the 1980s. However, the poor light output and low detection efficiency of TOF-capable detectors available at the time limited any gain in image quality achieved with these TOF-PET scanners over the traditional non-TOF PET scanners. The discovery of LSO and other Lu-based scintillators revived interest in TOF-PET and led to the development of a second generation of scanners with high sensitivity and spatial resolution in the mid-2000s. The introduction of the silicon photomultiplier (SiPM) has recently yielded a third generation of TOF-PET systems with unprecedented imaging performance. Parallel to these instrumentation developments, much progress has been made in the development of image reconstruction algorithms that better utilize the additional information provided by TOF. Overall, the benefits range from a reduction in image variance (SNR increase), through allowing joint estimation of activity and attenuation, to better reconstructing data from limited angle systems. In this work, we review these developments, focusing on three broad areas: 1) timing theory and factors affecting the time resolution of a TOF-PET system; 2) utilization of TOF information for improved image reconstruction; and 3) quantification of the benefits of TOF compared to non-TOF PET. Finally, we offer a brief outlook on the TOF-PET developments anticipated in the short and longer term. Throughout this work, we aim to maintain a clinically driven perspective, treating TOF as one of multiple (and sometimes competitive) factors that can aid in the optimization of PET imaging performance.

Silicon photomultiplier-based scintillation detectors for photon-counting CT: A feasibility study.

PURPOSE: The implementation of photon-counting detectors is widely expected to be the next breakthrough in X-ray computed tomography (CT) instrumentation. A small number of prototype scanners equipped with direct-conversion detectors based on room-temperature semiconductors, such as CdTe and CdZnTe (CZT), are currently installed at medical centers. Here, we investigate the feasibility of using silicon photomultiplier (SiPM)-based scintillation detectors in photon-counting computed tomography (PCCT) scanners, as a potential alternative to CdTe and CZT detectors. METHODS: We introduce a model that allows us to compute the expected energy resolution as well as the expected pulse shape and associated rate capability of SiPM-based PCCT detectors. The model takes into account SiPM saturation and optical crosstalk, because these phenomena may substantially affect the performance of SiPM-based PCCT detectors with sub-mm pixels. We present model validation experiments using a single-pixel detector consisting of a 0.9 × 0.9 × 1.0 mm3 LuAP:Ce scintillation crystal coupled to a 1 × 1 mm2 SiPM. We subsequently use the validated model to compute the expected performance of the fast scintillators LYSO:Ce, LuAP:Ce, and LaBr3 :Ce, coupled to currently available SiPMs, as well as to a more advanced SiPM prototype with improved dynamic range, for sub-mm pixel sizes. RESULTS: The model was found to be in good agreement with the validation experiments, both with respect to energy resolution and pulse shape. It shows how saturation progressively degrades the energy resolution of detectors equipped with currently available SiPMs as the pixel size decreases. Moreover, the expected pulse duration is relatively long (~200 ns) with these SiPMs. However, when LuAP:Ce and LaBr3 :Ce are coupled to the more advanced SiPM prototype, the pulse duration improves to less than 60 ns, which is in the same order of magnitude as pulses from CdTe and CZT detectors. It follows that sufficient rate capability can be achieved with pixel sizes of 400 µm or smaller. Moreover, LaBr3 :Ce detectors can provide an energy resolution of 11.5%-13.5% at 60 keV, comparable to CdTe and CZT detectors. CONCLUSIONS: This work provides first evidence that it may be feasible to develop SiPM-based scintillation detectors for PCCT that can compete with CdTe and CZT detectors in terms of energy resolution and rate capability.

Physics and technology of time-of-flight PET detectors.

The imaging performance of clinical positron emission tomography (PET) systems has evolved impressively during the last ∼15 years. A main driver of these improvements has been the introduction of time-of-flight (TOF) detectors with high spatial resolution and detection efficiency, initially based on photomultiplier tubes, later silicon photomultipliers. This review aims to offer insight into the challenges encountered, solutions developed, and lessons learned during this period. Detectors based on fast, bright, inorganic scintillators form the scope of this work, as these are used in essentially all clinical TOF-PET systems today. The improvement of the coincidence resolving time (CRT) requires the optimization of the entire detection chain and a sound understanding of the physics involved facilitates this effort greatly. Therefore, the theory of scintillation detector timing is reviewed first. Once the fundamentals have been set forth, the principal detector components are discussed: the scintillator and the photosensor. The parameters that influence the CRT are examined and the history, state-of-the-art, and ongoing developments are reviewed. Finally, the interplay between these components and the optimization of the overall detector design are considered. Based on the knowledge gained to date, it appears feasible to improve the CRT from the values of 200-400 ps achieved by current state-of-the-art TOF-PET systems to about 100 ps or less, even though this may require the implementation of advanced methods such as time resolution recovery. At the same time, it appears unlikely that a system-level CRT in the order of ∼10 ps can be reached with conventional scintillation detectors. Such a CRT could eliminate the need for conventional tomographic image reconstruction and a search for new approaches to timestamp annihilation photons with ultra-high precision is therefore warranted. While the focus of this review is on timing performance, it attempts to approach the topic from a clinically driven perspective, i.e. bearing in mind that the ultimate goal is to optimize the value of PET in research and (personalized) medicine.

Fast automated multi-criteria planning for HDR brachytherapy explored for prostate cancer.

We developed a fast and fully-automated, multi-criteria treatment planning workflow for high dose rate brachytherapy (HDR-BT). In this workflow, the patient-CT with catheter reconstructions and dwell positions are imported from the clinical TPS into a novel system for automated dwell time optimisation. The optimised dwell times are then imported into the clinical TPS. The aims of automation were (1) planner-independent, enhanced plan quality, (2) short optimisation times. Our in-house developed system for fully automated, multi-criteria external beam radiotherapy (EBRT) treatment planning (Erasmus-iCycle) was adapted for optimisation of HDR-BT dose distributions. The investigations were performed with planning CT scans with catheter reconstructions and delineations of twenty-five low- and intermediate-risk prostate cancer patients who were previously treated in our center with [Formula: see text] Gy HDR-BT. Automatically generated plans (autoplans) were compared to the corresponding clinical plans. All evaluations were performed in the clinical TPS. The requested 95% tumour coverage was obtained for all autoplans, while this was only observed in 23/25 clinical plans. All autoplans showed a consistent reduction of the [Formula: see text] for the highest prioritised OAR, the urethra. The average and maximum reductions were 6.3%-point and 12.1%-point of the prescribed dose, respectively. In addition, conformality of the autoplans was higher. The autoplans had slightly smaller delivery times. Autoplanning took on average 4.6 s, including computation of the dose kernels.

Correlations between the shifts in prompt gamma emission profiles and the changes in daily target coverage during simulated pencil beam scanning proton therapy.

The aim of this study was to investigate the feasibility of using prompt gamma (PG) ray emission profiles to monitor changes in dose to the planning target volume (PTV) during pencil beam scanning (PBS) proton therapy as a result of day-to-day variation in patient anatomy. For 11 prostate patients, we simulated treatment plan delivery using the patients' daily anatomy as observed in the planning CT and 7-9 control CT scans, including the detected PG profiles resulting from the 5%, 10%, and 20% most intense proton pencil beams. For each patient, we determined the changes in dosimetric parameters for the high- and low-dose PTVs between the simulations performed using the planning CT scan and the different control CT scans and correlated these to changes in the PG emission profiles. Changes in coverage of the high- and low-dose PTV correlated most strongly with the median and mean absolute PG emission profile shifts of the 5% most intense pencil beams, respectively. With a mean Pearson correlation coefficient of -0.76 (SD: 0.17) for the high-dose PTV and of -0.60 (SD: 0.51) for the low-dose PTV. We showed, as a proof of principle, that PG emission profiles obtained during PBS proton therapy could be used to detect changes in PTV coverage due to day-to-day anatomical variation.

Simulation of proton range monitoring in an anthropomorphic phantom using multi-slat collimators and time-of-flight detection of prompt-gamma quanta.

Prompt-gamma (PG) imaging has the potential for monitoring proton therapy in real time. Different approaches are investigated. We focus on developing multi-slat collimators to image PG quanta, aiming at optimizing collimator performance to detect deviations in treatment delivery. We investigated six different multi-slat configurations, which have either optimal (analytical) intrinsic spatial resolution at fixed efficiency, or otherwise; at different distances from the proton pencil-beam axis (15 cm-35 cm). We used Geant4 to simulate irradiations of the head (energy: 130 MeV) and pelvis (200 MeV) of an anthropomorphic phantom, with and without physiologic/morphologic or setup changes of clinical dosimetric relevance. The particles escaping the phantom were transported through each of these multi-slat configurations and the gamma counts profiles were recorded at the collimator exit. Median filtering was applied to the registered PG-profiles to mitigate the effects of septa shadowing and statistical fluctuations. Time-of-flight discrimination was used to enhance the signal-to-background ratio, which appeared crucial for 200 MeV irradiations. Visual detection of the artificially introduced changes was possible by comparing the PG to the depth-dose profiles. Moreover, 2 mm range shifts could be detected in the head irradiation case using a simple linear regression fit to the falloff of the PG-profile. The influence of changes in complex, patient-like dose distributions on the PG-profiles obtained with multi-slat collimation is first studied in this work, which further gives insight on collimator design optimization and highlights its potential and simplicity for detecting proton treatment deviations over a wide range of Bragg peak positions.

Sub-3 mm, near-200 ps TOF/DOI-PET imaging with monolithic scintillator detectors in a 70 cm diameter tomographic setup.

Recently, a monolithic scintillator detector for time-of-flight (TOF)/depth-of-interaction (DOI) positron emission tomography (PET) was developed. It has a detector spatial resolution of ~1.7 mm full-width-at-half-maximum (FWHM), a coincidence resolving time (CRT) of ~215 ps FWHM, and ~4.7 mm FWHM DOI resolution. Here, we demonstrate, for the first time, the imaging performance of this detector in a 70 cm diameter PET geometry. We built a tomographic setup representative of a whole-body clinical scanner, comprising two coaxially rotating arms, each carrying a detector module, and a central, rotating phantom table. The fully automated setup sequentially acquires all possible lines of response (LORs) of a complete detector ring, using a step-and-shoot acquisition approach. The modules contained 2 × 2 detectors, each detector consisting of a 32 mm × 32 mm × 22 mm LYSO crystal and a digital silicon photomultiplier (dSiPM) array. The system spatial resolution was assessed using a Na-22 point source at different radial distances in the field-of-view (FOV). Using 2D filtered back projection (2D FBP, non-TOF), tangential and radial spatial resolutions of ~2.9 mm FWHM were obtained at the center of the FOV. The use of DOI information resulted in almost uniform spatial resolution throughout the FOV up to a radial distance of 25 cm, where the radial and tangential resolution are ~3.3 mm FWHM and ~4.7 mm FWHM, respectively, whereas without DOI the resolution deteriorates to ~9 mm FWHM. Additional measurements were performed with a Na-22 filled Derenzo-like phantom at different locations within the FOV. Images reconstructed with a TOF maximum-likelihood expectation-maximization (TOF ML-EM) algorithm show that the system is able to clearly resolve 3 mm diameter hot rods up to 25 cm radial distance. The excellent and uniform spatial resolution, combined with an energy resolution of 10.2% FWHM and a CRT of ~212 ps FWHM, indicates a great potential for monolithic scintillators as practical high-performance detectors in TOF/DOI-PET systems.

Optimization of an ultralow-dose high-resolution pediatric PET scanner design based on monolithic scintillators with dual-sided digital SiPM readout: a simulation study.

The goal of this simulation study is the performance evaluation and comparison of six potential designs for a time-of-flight PET scanner for pediatric patients of up to about 12 years of age. It is designed to have a high sensitivity and provide high-contrast and high-resolution images. The simulated pediatric PET is a full ring scanner, consisting of 32 × 32 mm2 monolithic LYSO:Ce crystals coupled to digital silicon photomultiplier arrays. The six considered designs differ in axial lengths (27.2 cm, 54.4 cm and 102 cm) and crystal thicknesses (22 mm and 11 mm). The simulations are based on measured detector response data. We study two possible detector arrangements: 22 mm-thick crystals with dual-sided readout and 11 mm-thick crystals with back-sided readout. The six designs are simulated by means of the GEANT4 application for tomographic emission software, using the measured spatial, energy and time response of the monolithic scintillator detectors as input. The performance of the six designs is compared on the basis of four studies: (1) spatial resolution; (2) NEMA NU2-2012 sensitivity and scatter fraction (SF) tests; (3) non-prewhitening signal-to-noise ratio observer study; and (4) receiver operating characteristics analysis. Based on the results, two designs are identified as cost-effective solutions for fast and efficient imaging of children: one with 54.4 cm axial field-of-view (FOV) and 22 mm-thick crystals, and another one with 102 cm axial FOV and 11 cm-thick crystals. The first one has a higher center point sensitivity than the second one, but requires dual-sided readout. The second design has the advantage of allowing a whole-body scan in a single bed position acquisition. Both designs have the potential to provide an excellent spatial resolution (∼2 mm) and an ultra-high sensitivity (>100 cps [Formula: see text]).

Improved image quality using monolithic scintillator detectors with dual-sided readout in a whole-body TOF-PET ring: a simulation study.

We have recently built and characterized the performance of a monolithic scintillator detector based on a 32 mm × 32 mm × 22 mm LYSO:Ce crystal read out by digital silicon photomultiplier (dSiPM) arrays coupled to the crystal front and back surfaces in a dual-sided readout (DSR) configuration. The detector spatial resolution appeared to be markedly better than that of a detector consisting of the same crystal with conventional back-sided readout (BSR). Here, we aim to evaluate the influence of this difference in the detector spatial response on the quality of reconstructed images, so as to quantify the potential benefit of the DSR approach for high-resolution, whole-body time-of-flight (TOF) positron emission tomography (PET) applications. We perform Monte Carlo simulations of clinical PET systems based on BSR and DSR detectors, using the results of our detector characterization experiments to model the detector spatial responses. We subsequently quantify the improvement in image quality obtained with DSR compared to BSR, using clinically relevant metrics such as the contrast recovery coefficient (CRC) and the area under the localized receiver operating characteristic curve (ALROC). Finally, we compare the results with simulated rings of pixelated detectors with DOI capability. Our results show that the DSR detector produces significantly higher CRC and increased ALROC values than the BSR detector. The comparison with pixelated systems indicates that one would need to choose a crystal size of 3.2 mm with three DOI layers to match the performance of the BSR detector, while a pixel size of 1.3 mm with three DOI layers would be required to get on par with the DSR detector.

Towards monolithic scintillator based TOF-PET systems: practical methods for detector calibration and operation.

Gamma-ray detectors based on thick monolithic scintillator crystals can achieve spatial resolutions <2 mm full-width-at-half-maximum (FWHM) and coincidence resolving times (CRTs) better than 200 ps FWHM. Moreover, they provide high sensitivity and depth-of-interaction (DOI) information. While these are excellent characteristics for clinical time-of-flight (TOF) positron emission tomography (PET), the application of monolithic scintillators has so far been hampered by the lengthy and complex procedures needed for position- and time-of-interaction estimation. Here, the algorithms previously developed in our group are revised to make the calibration and operation of a large number of monolithic scintillator detectors in a TOF-PET system practical. In particular, the k-nearest neighbor (k-NN) classification method for x,y-position estimation is accelerated with an algorithm that quickly preselects only the most useful reference events, reducing the computation time for position estimation by a factor of ~200 compared to the previously published k-NN 1D method. Also, the procedures for estimating the DOI and time of interaction are revised to enable full detector calibration by means of fan-beam or flood irradiations only. Moreover, a new technique is presented to allow the use of events in which some of the photosensor pixel values and/or timestamps are missing (e.g. due to dead time), so as to further increase system sensitivity. The accelerated methods were tested on a monolithic scintillator detector specifically developed for clinical PET applications, consisting of a 32 mm × 32 mm × 22 mm LYSO : Ce crystal coupled to a digital photon counter (DPC) array. This resulted in a spatial resolution of 1.7 mm FWHM, an average DOI resolution of 3.7 mm FWHM, and a CRT of 214 ps. Moreover, the possibility of using events missing the information of up to 16 out of 64 photosensor pixels is shown. This results in only a small deterioration of the detector performance.

A 32 mm × 32 mm × 22 mm monolithic LYSO:Ce detector with dual-sided digital photon counter readout for ultrahigh-performance TOF-PET and TOF-PET/MRI.

New applications for positron emission tomography (PET) and combined PET/magnetic resonance imaging (MRI) are currently emerging, for example in the fields of neurological, breast, and pediatric imaging. Such applications require improved image quality, reduced dose, shorter scanning times, and more precise quantification. This can be achieved by means of dedicated scanners based on ultrahigh-performance detectors, which should provide excellent spatial resolution, precise depth-of-interaction (DOI) estimation, outstanding time-of-flight (TOF) capability, and high detection efficiency. Here, we introduce such an ultrahigh-performance TOF/DOI PET detector, based on a 32 mm × 32 mm × 22 mm monolithic LYSO:Ce crystal. The 32 mm × 32 mm front and back faces of the crystal are coupled to a digital photon counter (DPC) array, in so-called dual-sided readout (DSR) configuration. The fully digital detector offers a spatial resolution of ~1.1 mm full width at half maximum (FWHM)/~1.2 mm mean absolute error, together with a DOI resolution of ~2.4 mm FWHM, an energy resolution of 10.2% FWHM, and a coincidence resolving time of 147 ps FWHM. The time resolution closely approaches the best results (135 ps FWHM) obtained to date with small crystals made from the same material coupled to the same DPC arrays, illustrating the excellent correction for optical and electronic transit time spreads that can be achieved in monolithic scintillators using maximum-likelihood techniques for estimating the time of interaction. The performance barely degrades for events with missing data (up to 6 out of 32 DPC dies missing), permitting the use of almost all events registered under realistic acquisition conditions. Moreover, the calibration procedures and computational methods used for position and time estimation follow recently made improvements that make them fast and practical, opening up realistic perspectives for using DSR monolithic scintillator detectors in TOF-PET and TOF-PET/MRI systems.

The price of robustness; impact of worst-case optimization on organ-at-risk dose and complication probability in intensity-modulated proton therapy for oropharyngeal cancer patients.

PURPOSE: To quantify the impact of the degree of robustness against setup errors and range errors on organ-at-risk (OAR) dose and normal tissue complication probabilities (NTCPs) in intensity-modulated proton therapy for oropharyngeal cancer patients. MATERIAL AND METHODS: For 20 oropharyngeal cases (10 unilateral and 10 bilateral), robust treatment plans were generated using 'minimax' worst-case optimization. We varied the robustness against setup errors ('setup robustness') from 1 to 7mm and the robustness against range errors ('range robustness') from 1% to 7% (+1mm). We evaluated OAR doses and NTCP-values for xerostomia, dysphagia and larynx edema. RESULTS: Varying the degree of setup robustness was found to have a considerably larger impact than varying the range robustness. Increasing setup robustness from 1mm to 3, 5, and 7mm resulted in average NTCP-values to increase by 1.9, 4.4 and 7.5 percentage point, whereas they increased by only 0.4, 0.8 and 1.2 percentage point when increasing range robustness from 1% to 3%, 5% and 7%. The degree of setup robustness was observed to have a clinically significant impact in bilateral cases in particular. CONCLUSIONS: For oropharyngeal cancer patients, minimizing setup errors should be given a higher priority than minimizing range errors.

Comparative Characterization Study of a LaBr3(Ce) Scintillation Crystal in Two Surface Wrapping Scenarios: Absorptive and Reflective.

The properties of a 50 mm × 50 mm × 30 mm monolithic LaBr3:Ce scintillator crystal coupled to a position-sensitive multi-anode photomultiplier (PMT, Hamamatsu H9500), representing the absorbing detector of a Compton camera under study for online ion (proton) beam range verification in hadron therapy, was evaluated in combination with either absorptive or reflective crystal surface coating. This study covered an assessment of the energy and position-dependent energy resolution, exhibiting a factor of 2.5-3.5 improvement for the reflectively wrapped crystal at 662 keV. The spatial dependency was investigated using a collimated (137)Cs source, showing a steep degradation of the energy resolution at the edges and corners of the absorptively wrapped crystal. Furthermore, the time resolution was determined to be 273 ps (FWHM) and 536 ps (FWHM) with reflective and absorptive coating, respectively, using a (60)Co source. In contrast, the light spread function (LSF) of the light amplitude distribution on the PMT segments improved for the absorptively wrapped detector. Both wrapping modalities showed almost no differences in the energy-dependent photopeak detection efficiency.

Time-based position estimation in monolithic scintillator detectors.

Gamma-ray detectors based on bright monolithic scintillation crystals coupled to pixelated photodetectors are currently being considered for several applications in the medical imaging field. In a typical monolithic detector, both the light intensity and the time of arrival of the earliest scintillation photons can be recorded by each of the photosensor pixels every time a gamma interaction occurs. Generally, the time stamps are used to determine the gamma interaction time while the light intensities are used to estimate the 3D position of the interaction point. In this work we show that the spatio-temporal distribution of the time stamps also carries information on the location of the gamma interaction point and thus the time stamps can be used as explanatory variables for position estimation. We present a model for the spatial resolution obtainable when the interaction position is estimated using exclusively the time stamp of the first photon detected on each of the photosensor pixels. The model is shown to be in agreement with experimental measurements on a 16 mm × 16 mm × 10 mm LSO : Ce,0.2%Ca crystal coupled to a digital photon counter (DPC) array where a spatial resolution of 3 mm (root mean squared error) is obtained. Finally we discuss the effects of the main parameters such as scintillator rise and decay time, light output and photosensor single photon time resolution and pixel size.

Time-resolved imaging of prompt-gamma rays for proton range verification using a knife-edge slit camera based on digital photon counters.

Proton range monitoring may facilitate online adaptive proton therapy and improve treatment outcomes. Imaging of proton-induced prompt gamma (PG) rays using a knife-edge slit collimator is currently under investigation as a potential tool for real-time proton range monitoring. A major challenge in collimated PG imaging is the suppression of neutron-induced background counts. In this work, we present an initial performance test of two knife-edge slit camera prototypes based on arrays of digital photon counters (DPCs). PG profiles emitted from a PMMA target upon irradiation with a 160 MeV proton pencil beams (about 6.5 × 10(9) protons delivered in total) were measured using detector modules equipped with four DPC arrays coupled to BGO or LYSO : Ce crystal matrices. The knife-edge slit collimator and detector module were placed at 15 cm and 30 cm from the beam axis, respectively, in all cases. The use of LYSO : Ce enabled time-of-flight (TOF) rejection of background events, by synchronizing the DPC readout electronics with the 106 MHz radiofrequency signal of the cyclotron. The signal-to-background (S/B) ratio of 1.6 obtained with a 1.5 ns TOF window and a 3 MeV-7 MeV energy window was about 3 times higher than that obtained with the same detector module without TOF discrimination and 2 times higher than the S/B ratio obtained with the BGO module. Even 1 mm shifts of the Bragg peak position translated into clear and consistent shifts of the PG profile if TOF discrimination was applied, for a total number of protons as low as about 6.5 × 10(8) and a detector surface of 6.6 cm × 6.6 cm.

Time estimation with multichannel digital silicon photomultipliers.

Accuracy in timemark estimation is crucial for time-of-flight positron emission tomography, in order to ensure high quality images after reconstruction. Since the introduction of multichannel digital silicon photomultipliers, it is possible to acquire several photoelectron timestamps for each individual gamma event. We study several timemark estimators based on multiple photoelectron timestamps by means of a comprehensive statistical model. In addition, we calculate the MSE of the estimators in comparison to the Cramér-Rao lower bound as a function of the system design parameters. We investigate the effect of skipping some of the photoelectron timestamps, which is a direct consequence of the limited number of time-to-digital converters and we propose a technique to compensate for this effect. In addition, we carry out an extensive analysis to evaluate the influence of dark counts on the detector timing performance. Moreover, we investigate the improvement of the timing performance that can be obtained with dark count filtering and we propose an appropriate filtering method based on measuring the time difference between sorted timestamps. Finally, we perform a full Monte Carlo simulation to compare different timemark estimators by exploring several system design parameters. It is demonstrated that a simple weighted-average estimator can achieve a comparable performance as the more complex maximum likelihood estimator.