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.
Abstract-Multiple time-to-digital converters coupled with silicon photomultipliers allow to timestamp several light photons generated by a scintillation event. Multichannel digital silicon photomultipliers opened the possibility to estimate a gammaphoton time mark by using several photoelectrons timestamps.We studied the already-existing statistics models of photoeletron time-stamping generation, while extending the current models by adding the skipping effect. Which accounts for the inability of the system to timestamp a continuous set of photoelectrons.In addition, we proposed two multiple photoelectron timemark estimators based on the best linear unbiased and the maximum likelihood estimation methods. We calculated the Cramér Rao lower bound for several system parameter and compared it to the proposed estimators' performance. We concluded that under certain system configurations the proposed estimators are efficient.Moreover, we investigated the effect of the dark count rate on the timing performance. Also, we introduced a filtering method that is based on measuring the time distance between adjacent timestamps. We performed a full Monte Carlo simulation to evaluate the proposed filter efficiency.Finally, we performed a full Monte Carlo simulation to compare the timemark estimators' performance. We concluded that the best linear unbiased estimator is as efficient as the maximum likelihood estimator. In addition, it was verified that the multichannel digital silicon photomultipliers have a stronger tolerance to dark counts in comparison current digital silicon photomultiplier architectures. I. CRAMÉR RAO LOWER BOUNDThe Cramér Rao Lower Bound (CRLB) determines a maximum estimation performance for problems that satisfy certain regularity conditions. The CRLB for the particular case of gamma-photon timemmark estimation was analyzed byIn single-photoelectron (q th ) timestamp estimation there is only one possible unbiased estimator, which is given bŷwhere t q represents the time-of-registration of the q th photoelectron and A is expressed asThe root mean square error (root-MSE) of the singlephotoelectron estimation method is equal to the square root of the variance of the PDF (probability density function) of the time-of-registration of the q th photoelectron [2].The maximum timing performance in the singlephotoelectron estimation case is determined by Eq. (3). We calculated the performance as a function of the photoelectron order for two different system parameter configurations (see Figs. 1 and 2) called datasets I and II. Dataset I is characterized by 100 ps FWHM jitter and 300 photoelectrons; dataset II by 700 ps FWHM jitter and 3800 photoelectrons. The scintillation decay constants were the same for both cases (LSO with properties [1]). We also calculated the CRLB (intrinsic limit) for each data set [2], [1]. II. MAXIMUM LIKELIHOOD ESTIMATIONThe likelihood function that corresponds to the t 1:Q timestamps of the first Q photoelectrons, when estimating the location parameter (T 0 ) is defined bywhere R repres...
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.