Mean and Variance of Photon Counting with Deadtime
|
|
- Sabina O’Neal’
- 5 years ago
- Views:
Transcription
1 Mean and Variance of Photon Counting with Deadtime Daniel F. Yu and Jeffrey A. Fessler 445 EECS Bldg., University of Michigan, Ann Arbor, MI Abstract The statistics of photon counting by systems affected by deadtime are potentially important for statistical image reconstruction methods. We present a new way of analyzing the moments of the counting process for a counter system affected by various models of deadtime related to PET and SPECT imaging. We derive simple and exact expressions for the first and second moments of the number of recorded events under various models. From our mean expression for a SPECT deadtime model, we derive a simple estimator for the actual intensity of the underlying Poisson process; simulations show that our estimator is unbiased even for extremely high count rates. I. INTRODUCTION Every photon counting system exhibits a characteristic called deadtime. Since the pulses produced by a detector have finite time duration, if a second pulse occurs before the first has disappeared, the two pulses will overlap to form a single distorted pulse []. Depending on the system, one or both arrivals will be lost. In PET or SPECT scanners, the length of pulse resolving time, often just called deadtime, denoted ç, is around 2çs. Counting systems are usually classified into two categories: nonparalyzable (type I) or paralyzable (type II). In a nonparalyzable system, each recorded photon produces a deadtime of length ç; if an arrival is recorded at t, then any arrival from t to t + ç will not be recorded. In a paralyzable system, each photon arrival, whether recorded or not, produces a deadtime of length ç; if there is an arrival at t,thenanyarrival from t to t + ç will not be recorded. In some SPECT systems [2], we encounter a third model that is similar to the paralyzable model: if two photons arrive within ç of each other, then neither photon will be recorded (e.g., due to pulse pile-up); we call this the type III model. The asymptotic moments of the nonparalyzable model are well known [3]. For the paralyzable model, the exact expression for the mean of the number of recorded events from time to t, denoted Y ètè, has been derived previously [4]. However, for the type III model, only an approximate expression for the mean number of recorded events has been derived [2]. In this paper, we derive the exact mean and variance expressions of Y ètè for both type II and type III models. II. STATISTICAL ANALYSIS OF DEADTIME We define a photon arrival to mean a photon interacting with the scintillator with sufficient deposited energy to trigger detection. The photon arrival process N ètè counts the number of arrivals during the time interval è;të, and the photon recording process Y ètè counts the number of recorded events. We assume that N ètè is a homogeneous Poisson process with rate ç (photon arrivals per unit time) which stays constant with time. We also assume, for the sake of simplicity, that ç is known and deterministic. A. Asymptotic Analysis via Renewal Theory The counting processes in all three types of systems discussed above are examples of renewal processes [3], and renewal theory has been the classical basis for deadtime analysis [5]. A renewal process involves recurrent patterns connected with repeated trials. Roughly speaking, if after each occurrence of a pattern E, the random process starts from scratch in the sense that the trials following an occurrence of E form a replica of the whole process, then the process qualifies as a renewal process. If we define E to be the state of the counter is ready to record the next photon arrival, then after each occurrence of E, the counting process is statistically equivalent. A very useful random variable to define is T E, the waiting time between one renewal and the next (renewal here means return to E). Note that in the context of photon counting system, with E defined as above, the number of renewals from to t is almost exactly the number of recorded events from to t. IfT E has ensemble mean ç E and variance çe 2, then the number of renewals from to t, ~Y ètè, is asymptotically Gaussian distributed [6] [3] with the following moments: Eë ~ Y ètèë ç tç E ; Varë ~ Y ètèë ç tç 2 Eç 3 E; () where ç indicates that the ratio of the two sides tends to unity as tç E!. Hence asymptotically, the mean and For type III deadtime, we define renewal as return to E after recording an event.
2 variance of the waiting time between renewals forms a sort of duality relationship with the mean and variance of the number of renewals. For the other two deadtime models, if we try to derive EëY ètèë from EëT E ë, it is much more difficult to obtain a simple closed form expression because the EëT E ë we get is probably an infinite sum and it is often not easy to obtain every term in this sum; the variance of T E is even more complicated. Therefore, in the following section, we describe a new approach for deriving the moments of counting processes. B. Exact Mean and Variance of Counting Processes We first consider a general counting process Y where Y èt ;t 2 è denotes the number of recorded events during the time interval èt ;t 2 ë and Y ètè is a shorthand for Y è;tè. We define the instantaneous rate æ : R! ë; è of the process Y ètè as: æèsè, EëY ès + æè, Y èsèëæ; (2) æ! and the instantaneous second moment æ : R! ë; è as: æèsè, EëèY ès + æè, Y èsèè 2 ëæ: (3) æ! We also define the correlation function æ : R 2! ë; è as: æès ;s 2 è, EëèY ès + æ è, Y ès èè æ æ ;æ 2! We assume 2 èy ès 2 + æ 2 è, Y ès 2 èèëèæ æ 2 è: (4) (i) æ and æ are well-defined ç-almost everywhere, and æ is well defined ç 2 -almost everywhere, and æ and æ are integrable with respect to ç and ç 2 over any finite interval and rectangle, respectively; (ii) EëY ès; s + æèëæ and EëY 2 ès; s + æèëæ are uniformly bounded for all s and æ 2 è; è; (iii) EëY ès ;s +æ èy ès 2 ;s 2 +æ 2 èëèæ æ 2 è is uniformly bounded for all s, s 2,andæ ;æ 2 2 è; è such that ès ;s + æ è ë ès 2 ;s 2 + æ 2 è ;. These Assumptions hold for a wide variety of counting processes, including any homogeneous Poisson process with finite intensity. 2 ç and ç2 denote Lebesgue measures on R and R 2, respectively. For analysis purposes, we artificially divide the time interval ë;të into n segments of length æ each, i.e., t næ. We have Y ètè EëY ètèë i i Y èiæ; èi +èæè; (5) EëY èiæ; èi +èæèë; (6) R f æèsèds; (7) where we define the following piecewise constant function: f æ èsè, é : EëY èjæ; èj +èæèëæ; if s 2 èjæ; èj +èæë; ç j ç n, ; otherwise. () Since æètè is well-defined almost everywhere in the interval ë;të and EëY ès; s + æèëæ is uniformly bounded, by the Lebesgue Dominated Convergence theorem (LDCT) [7], æ! R f æèsèdçèsè R f æ èsèdçèsè æ! æèsèds: (9) Hence, we have the following simple general expression for the mean of the counting process in terms of its instantaneous rate: EëY ètèë æèsèds: () We consider the second moment by a similar argument: EëY 2 ètèë Eëè i i Y èiæ; èi +èæèè 2 ë EëY 2 èiæ; èi +èæèèë + 2 n,2 i ji+ EëY èiæ; èi +èæèy èjæ; èj +èæèë () R g æèsèdçèsè +2 h æ ès ;s 2 èdç 2 ès ;s 2 è; (2) R 2 where we define the following piecewise constant functions: g æ èsè, é : EëY 2 èjæ; èj +èæèëæ; if s 2 èjæ; èj +èæë and ç j ç n, ; otherwise, (3)
3 and EëY èiæ; èi +èæèæ if s 2 èiæ; èi +èæë, é Y èjæ; s 2 2 èjæ; èj +èæë, h æ ès ;s 2 è, èj +èæèëæ 2 ; ç i ç n, 2, and i +ç j ç n, é: ; otherwise. (4) Since æ is well-defined almost everywhere in ë;tëæë;të and EëY ès ;s + æèy ès 2 ;s 2 + æèëæ 2 is uniformly bounded, by LDCT and Fubini s Theorem [7], æ! R 2 h æ ès ;s 2 èdç 2 ès ;s 2 è h æ ès ;s 2 èdç 2 ès ;s 2 è R 2 æ! æès ;s 2 èds ds 2 : (5) s 2 Similarly, one can show that æ! R g æèsèdçèsè æèsèds: (6) Thus using (2), (5), and (6), we have the following general expression for the second moment of Y ètè: EëY 2 ètèë æèsèds +2 æès ;s 2 èds 2 ds : (7) s In the context of counting processes with deadtime, which includes all random processes considered in this paper, the process satisfies this additional assumption: (iv) there exists a positive æ such that æ 2 è;æ è, Y ès; s + æè ç. This assumption greatly simplifies the derivations for the moments of counting processes affected by deadtime, since for æéæ éç, EëY 2 ès; s + æèë EëY ès; s + æèë () using 2 and 2,so æèsè æèsè: (9) Thus we obtain the following corollary of (7) for random processes satisfying assumptions (i) to (iv): EëY 2 ètèë EëY ètèë + 2 æès ;s 2 èds 2 ds : (2) s Furthermore, if Y ètè has stationary increments, then æèsè is constant and æès ;s 2 èæè;s 2, s è and we can further simplify the results () and (2) to the following: EëY ètèë æt (2) EëY 2 ètèë æt +2 èt, sèæè;sèds: (22) III. SINGLE PHOTON COUNTING A. Mean and Variance of Recorded Singles Counts, Model Type II First we consider the paralyzable model in which if the waiting time for a photon arrival is less than ç, then this photon is not recorded. We derive the mean and variance of Y ètè, the number of recorded events from time to time t. We observe that Y ètè inherits the stationary increment property of the arrival process N ètè. Wefirstderive EëY è;æèë, where we pick æéçsuch that the number of recorded events during è;æë is either or. Let T denote the time of the first photon arrival after time ; it is exponentially distributed. If there is an arrival at T s, é s é æ, and there is no arrival between s, ç and s (in fact, we only need to make sure there is no arrival between s, ç and, i.e., N èè, N ès, ç è,since the first arrival after occurs at s), then there will be a recorded event during the interval è;æë. Thus EëY è;æèë æ æ PëY è;æèjt sëf T èsèds PëN ès, ç; è jt sëf T èsèds e,çèç,sè çe,çs ds çæe,çç : (23) Hence by the definition given in (2), the instantaneous rate of Y ètè is æ çe,çç ; (24) and by (2), we easily obtain the following result (e.g., []), EëY ètèë çte,çç : (25) The variance of Y ètè for the type II model is (see Appendix A): VarëY ètèë çte,çç è, è2çç, çç 2 tèe,çç è: (26) Figure shows the mean and variance of the singles count for a detector affected by deadtime of type II. Since the mean and variance can differ greatly, Y ètè is not Poisson. B. Mean and Variance of Recorded Singles Counts, Model Type III Now we turn to the type of system described in [2], in which if the waiting time for a photon arrival is less than ç, then neither this photon nor the previous photon
4 . 2 x 5 Ideal mean: E[N(t)] Type II mean: E[Y(t)] Type II variance: Var[Y(t)]. 2 x 5 Ideal mean: E[N(t)] Type III mean: E[Y(t)] Type III variance: Var[Y(t)] E[N(t)], E[Y(t)], Var[Y(t)].2. E[N(t)], E[Y(t)], Var[Y(t)] x 5 Fig. Mean and variance for paralyzable (type II) systems, with t s;ç 2çs. will be recorded. We again observe that Y ètè inherits the stationary increment property of the arrival process N ètè. We first derive EëY è;æèë, where we pick æ é ç such that the number of recorded events during è;æë is still either or. Hence, EëY è;æèë æ æ æ PëY è;æèjt sëf T èsèds PëN ès, ç; è ëpëès; s + ç èëf T èsèds e,çèç,sè e,çç çe,çs ds çæe,ç2ç : (27) Hence for this system, the instantaneous rate as defined in (2) is æ çe,ç2ç ; (2) and by (2), the expected number of recorded events for a type III system is exactly: EëY ètèë çte,ç2ç : (29) The variance of Y ètè for the type III model is (we omit the derivation due to space constraints): VarëY ètèë çte,ç2ç +2e,3çç èçt, çç, è +e,4çç è4ç 2 ç 2, 4ç 2 tç +2, 2çt +4çç è: (3) Figure 2 shows the (exact) mean and variance of the singles count Y ètè for type III systems. Again Y ètè is not Poisson x 5 Fig. 2 Mean and variance for type III systems, with t s;ç 2çs. IV. COUNT RATE CORRECTION FOR SYSTEM TYPE III For a quantitatively accurate reconstruction, we must correct for the effect of deadtime to avoid underestimation of source activity. For type III systems, Engeland et al [2] proposed the following correction formula, ^ç Y t è + 2Y t ç + 6Y t ç 2 è; (3) 2 which they obtained by solving an approximate mean waiting time expression up to second order in ç by means of the expansion ç a + bç + cç 2. We propose to estimate the true count rate by solving numerically our exact expression (29), i.e., solve Y t ^çe,2^çç (32) for ^ç given Y and t. One could solve analytically the exact mean waiting time expression up to second order in ç, which yields exactly the same estimator as (3), but this estimator does not the mean waiting time expression exactly. Figure 3 compares our new estimator (32) and the estimator proposed in [2]. It shows that our new estimator is unbiased even at very high count rates. The error bars are not shown in the figure as they are smaller than the plotting symbols. When t is large, the standard deviation is very small when compared to the mean of Y ètè, thus these estimates have extremely small standard deviations. By solving (32) numerically, we obtain essentially perfect deadtime correction for a type III system. V. DISCUSSION We have analyzed the mean and variance of the recorded singles counts for two distinct models of deadtime. In both cases, the variance can be significantly
5 x our estimate of uncorrected estimate of Engelands estimate x 4 Fig. 3 2 realizations, with t s; ç 2çs. less than the mean, indicating that the counting statistics are not Poisson in the presence of deadtime. Deadtime losses can be significant in practical SPECT and PET systems, particularly in fully 3D PET imaging and in SPECT transmission measurements with a scanning line source. The count rates for a detector block (PET) or detector zone (SPECT) can be significant enough to yield non-poisson statistics for the total counts recorded by the block or zone. However 3, in the practical situations that we are aware of, the count rates for individual detector elements within the block or zone are usually not high enough to correspond to significant deadtime losses. Even though the variance of the counts recorded by a block can be significantly lower than the mean, the variance of the counts recorded by an individual detector element is nevertheless quite close to the mean and likely to be well approximated by a Poisson distribution. Furthermore, the correlation between individual detectors will be fairly small. Thus it appears that statistical image reconstruction based on Poisson models, while certainly not optimal, should be adequate in practice even under fairly large deadtime losses, provided the deadtime loss factor is included in the system matrix. We must add one caveat to this conclusion however. Although pairs of individual detectors have small correlation, the correlation coefficient between the sum of one group of detectors and the sum of all other detectors in a block may not be small in the presence of deadtime. The effect of such correlations on image reconstruction algorithms is unknown and may deserve further investigation. Another natural extension of this work would be to consider systems with random resolving times ç. As long as the minimum resolving time is greater than zero, assumption (iv) would still hold and the derivations would be similar. 3 Due to space constraints, we omit detailed analysis and only present our conclusions. VI. APPENDI A We derive the variance of Y ètè for deadtime model II, the paralyzable model. We first derive æè;sè. We consider two cases. CASE : éséç We pick æ such that é æ é s é s + æ é ç. Two recorded events cannot correspond to photons that arrived within ç of each other. Hence for é s é ç, EëY è;æèy ès; s + æèë, and by the definition given in (2): æè;sè. CASE 2: çését We pick æ such that æéçand s +æ étand æés,ç. For séç, Y è;æè and Y ès; s + æè are statistically independent, since the event there is an arrival during è;æë is statistically independent from the event there is an arrival during ès; s + æë, because they are at least ç apart in time. Hence by (23), EëY è;æèy ès; s + æèë E 2 ëy è;æèë èçæe,çç è 2 ; (33) and æè;sèèçe,çç è 2 : (34) Combining the above two cases and using (22) yields EëY 2 ètèë æt +2 ç èt, sèèçe,çç è 2 ds çte,çç +ëèt, ç èèçe,çç èë 2 : (35) Using VarëY ètèë EëY 2 ètèë, E 2 ëy ètèë, with (25) and (35), and simplifying yields (26). VII. REFERENCES [] J. A. Sorenson and M. E. Phelps, Physics in nuclear medicine, Saunders, Philadelphia, 2 edition, 97. [2] U. Engeland, T. Striker, and H. Luig, Count-rate statistics of the gamma camera, Phys. Med. Biol., vol. 43, no., pp , October 99. [3] W. Feller, An Introduction to Probability Theory and Its Applications, John Wiley and Sons, 96. [4] F. Carloni, A. Corberi, M. Marseguerra, and C. M. Porceddu, Asymptotic method of dead-time correction in Poissonian distribution, Nucl Instrum Methods, vol. 7, no., pp. 7 76, February 97. [5] J. Libert, Statistique de comptage avec distribution uniofmre de temps mort, Nucl.Instr.Meth., vol. 5, no. 3, pp , May 97. [6] D. R. Cox, Renewal Theory, Methuen and Co., 962. [7] A. M. Bruckner, J. B. Bruckner, and B. S. Thomson, Real Analysis, Prentice Hall, 997. This work was supported in part by NIH grants CA- 67 and CA
Explicit evaluation of the transmission factor T 1. Part I: For small dead-time ratios. by Jorg W. MUller
Rapport BIPM-87/5 Explicit evaluation of the transmission factor T (8,E) Part I: For small dead-time ratios by Jorg W. MUller Bureau International des Poids et Mesures, F-930 Sevres Abstract By a detailed
More informationObjective Functions for Tomographic Reconstruction from. Randoms-Precorrected PET Scans. gram separately, this process doubles the storage space for
Objective Functions for Tomographic Reconstruction from Randoms-Precorrected PET Scans Mehmet Yavuz and Jerey A. Fessler Dept. of EECS, University of Michigan Abstract In PET, usually the data are precorrected
More informationScan Time Optimization for Post-injection PET Scans
Presented at 998 IEEE Nuc. Sci. Symp. Med. Im. Conf. Scan Time Optimization for Post-injection PET Scans Hakan Erdoğan Jeffrey A. Fessler 445 EECS Bldg., University of Michigan, Ann Arbor, MI 4809-222
More informationRadiation Detection and Measurement
Radiation Detection and Measurement June 2008 Tom Lewellen Tkldog@u.washington.edu Types of radiation relevant to Nuclear Medicine Particle Symbol Mass (MeV/c 2 ) Charge Electron e-,! - 0.511-1 Positron
More informationWeek 4: Chap. 4 Basic Detector Principles
Week 4: Chap. 4 Basic Detector Principles General use of Statistical Distributions Basic Detector Principles -- Current Mode -- Pulse Mode --- resolution ---- Fano Factor --- efficiency --- dead time Ion
More informationSystem Modeling for Gamma Ray Imaging Systems
System Modeling for Gamma Ray Imaging Systems Daniel J. Lingenfelter and Jeffrey A. Fessler COMMUNICATIONS & SIGNAL PROCESSING LABORATORY Department of Electrical Engineering and Computer Science The University
More informationAn Introduction to Probability Theory and Its Applications
An Introduction to Probability Theory and Its Applications WILLIAM FELLER (1906-1970) Eugene Higgins Professor of Mathematics Princeton University VOLUME II SECOND EDITION JOHN WILEY & SONS Contents I
More informationSample Spectroscopy System Hardware
Semiconductor Detectors vs. Scintillator+PMT Detectors Semiconductors are emerging technology - Scint.PMT systems relatively unchanged in 50 years. NaI(Tl) excellent for single-photon, new scintillation
More informationOn background radiation and effective source strength. by Jorg,~. Müller. Bureau International des Poids et Mesures, F Sèvres
BIPM Working Party Note 226 On background radiation and effective source strength by Jorg,~. Müller Bureau International des Poids et Mesures, F-92310 Sèvres Any measurement of the activity of a radiation
More informationA conceptual interpretation of the renewal theorem with applications
Risk, Reliability and Societal Safety Aven & Vinnem (eds) 2007 Taylor & Francis Group, London, ISBN 978-0-415-44786-7 A conceptual interpretation of the renewal theorem with applications J.A.M. van der
More informationPositron Emission Tomography
Positron Emission Tomography Presenter: Difei Wang June,2018 Universität Bonn Contents 2 / 24 1 2 3 4 Positron emission Detected events Detectors and configuration Data acquisition Positron emission Positron
More informationLecture 1: Brief Review on Stochastic Processes
Lecture 1: Brief Review on Stochastic Processes A stochastic process is a collection of random variables {X t (s) : t T, s S}, where T is some index set and S is the common sample space of the random variables.
More informationDISCRETE STOCHASTIC PROCESSES Draft of 2nd Edition
DISCRETE STOCHASTIC PROCESSES Draft of 2nd Edition R. G. Gallager January 31, 2011 i ii Preface These notes are a draft of a major rewrite of a text [9] of the same name. The notes and the text are outgrowths
More informationList of Nuclear Medicine Radionuclides. Nuclear Medicine Imaging Systems: The Scintillation Camera. Crystal and light guide
Nuclear Medicine Imaging Systems: The Scintillation Camera List of Nuclear Medicine Radionuclides Tc99m 140.5 kev 6.03 hours I-131 364, 637 kev 8.06 days I-123 159 kev 13.0 hours I-125 35 kev 60.2 days
More informationChapter 5. Statistical Models in Simulations 5.1. Prof. Dr. Mesut Güneş Ch. 5 Statistical Models in Simulations
Chapter 5 Statistical Models in Simulations 5.1 Contents Basic Probability Theory Concepts Discrete Distributions Continuous Distributions Poisson Process Empirical Distributions Useful Statistical Models
More informationRESPONSE FUNCTION STUDY FOR ENERGY TO LIGHT CONVERSION IN ORGANIC LIQUID SCINTILLATORS
RESPONSE FUNCTION STUDY FOR ENERGY TO LIGHT CONVERSION IN ORGANIC LIQUID SCINTILLATORS S. Prasad *, A. Enqvist, S. D. Clarke, S. A. Pozzi, E. W. Larsen 1 Department of Nuclear Engineering and Radiological
More informationTHE COMPTON EFFECT Last Revised: January 5, 2007
B2-1 THE COMPTON EFFECT Last Revised: January 5, 2007 QUESTION TO BE INVESTIGATED: How does the energy of a scattered photon change after an interaction with an electron? INTRODUCTION: When a photon is
More informationWeek 7: Ch. 10 Spec. w/ Scintillation Ctrs. Photomultiplier Devices
Week 7: Ch. 0 Spec. w/ Scintillation Ctrs. multiplier Devices Spectroscopy with Scint. Counters -- gamma-ray interactions, reprise -- observed spectra --- spectral components, backscatter --- summing --
More informationA.V. SAVKIN AND I.R. PETERSEN uncertain systems in which the uncertainty satisæes a certain integral quadratic constraint; e.g., see ë5, 6, 7ë. The ad
Journal of Mathematical Systems, Estimation, and Control Vol. 6, No. 3, 1996, pp. 1í14 cæ 1996 Birkhíauser-Boston Robust H 1 Control of Uncertain Systems with Structured Uncertainty æ Andrey V. Savkin
More informationMaximum-Likelihood Deconvolution in the Spatial and Spatial-Energy Domain for Events With Any Number of Interactions
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 59, NO. 2, APRIL 2012 469 Maximum-Likelihood Deconvolution in the Spatial and Spatial-Energy Domain for Events With Any Number of Interactions Weiyi Wang, Member,
More informationRadionuclide Imaging MII Positron Emission Tomography (PET)
Radionuclide Imaging MII 3073 Positron Emission Tomography (PET) Positron (β + ) emission Positron is an electron with positive charge. Positron-emitting radionuclides are most commonly produced in cyclotron
More informationIntroduction to the Mathematics of Medical Imaging
Introduction to the Mathematics of Medical Imaging Second Edition Charles L. Epstein University of Pennsylvania Philadelphia, Pennsylvania EiaJTL Society for Industrial and Applied Mathematics Philadelphia
More informationASM Study Manual for Exam P, Second Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA Errata
ASM Study Manual for Exam P, Second Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA (krzysio@krzysio.net) Errata Effective July 5, 3, only the latest edition of this manual will have its errata
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationL n = l n (π n ) = length of a longest increasing subsequence of π n.
Longest increasing subsequences π n : permutation of 1,2,...,n. L n = l n (π n ) = length of a longest increasing subsequence of π n. Example: π n = (π n (1),..., π n (n)) = (7, 2, 8, 1, 3, 4, 10, 6, 9,
More informationEfficient Algorithms for Pulse Parameter Estimation, Pulse Peak Localization And Pileup Reduction in Gamma Ray Spectroscopy M.W.Raad 1, L.
Efficient Algorithms for Pulse Parameter Estimation, Pulse Peak Localization And Pileup Reduction in Gamma Ray Spectroscopy M.W.Raad 1, L. Cheded 2 1 Computer Engineering Department, 2 Systems Engineering
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More informationGEOMETRIC -discrete A discrete random variable R counts number of times needed before an event occurs
STATISTICS 4 Summary Notes. Geometric and Exponential Distributions GEOMETRIC -discrete A discrete random variable R counts number of times needed before an event occurs P(X = x) = ( p) x p x =,, 3,...
More informationRadiation Detection and Measurement
Radiation Detection and Measurement Third Edition Glenn F. Knoll Professor of Nuclear Engineering and Radiological Sciences University of Michigan Ann Arbor, Michigan John Wiley & Sons, Inc. New York/Chichester/Weinheim/Brisbane/Toronto/Singapore
More informationSynthesis of Reinforcement Learning, Neural Networks, and PI Control Applied to a Simulated Heating Coil
Synthesis of Reinforcement Learning, Neural Networks, and PI Control Applied to a Simulated Heating Coil Charles W. Anderson 1, Douglas C. Hittle 2, Alon D. Katz 2, and R. Matt Kretchmar 1 1 Department
More informationCommunication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi
Communication Engineering Prof. Surendra Prasad Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture - 41 Pulse Code Modulation (PCM) So, if you remember we have been talking
More informationBirth-Death Processes
Birth-Death Processes Birth-Death Processes: Transient Solution Poisson Process: State Distribution Poisson Process: Inter-arrival Times Dr Conor McArdle EE414 - Birth-Death Processes 1/17 Birth-Death
More informationCross-Selling in a Call Center with a Heterogeneous Customer Population. Technical Appendix
Cross-Selling in a Call Center with a Heterogeneous Customer opulation Technical Appendix Itay Gurvich Mor Armony Costis Maglaras This technical appendix is dedicated to the completion of the proof of
More informationRobust Maximum- Likelihood Position Estimation in Scintillation Cameras
Robust Maximum- Likelihood Position Estimation in Scintillation Cameras Jeffrey A. Fessler: W. Leslie Rogers, and Neal H. Clinthorne Division of Nuclear Medicine, University of Michigan ABSTRACT The classical
More informationRecap. Probability, stochastic processes, Markov chains. ELEC-C7210 Modeling and analysis of communication networks
Recap Probability, stochastic processes, Markov chains ELEC-C7210 Modeling and analysis of communication networks 1 Recap: Probability theory important distributions Discrete distributions Geometric distribution
More informationRadiation Detection and Measurement
Radiation Detection and Measurement Fourth Edition Glenn F. Knoll Professor Emeritus of Nuclear Engineering and Radiological Sciences University of Michigan Ann Arbor, Michigan WILEY John Wiley & Sons,
More informationin Cross-Section Data
Sensitivity of Photoneutron Production to Perturbations in Cross-Section Data S. D. Clarke Purdue University, West Lafayette, Indiana S. A. Pozzi University of Michigan, Ann Arbor, Michigan E. Padovani
More informationRandom Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg
Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws Symeon Chatzinotas February 11, 2013 Luxembourg Outline 1. Random Matrix Theory 1. Definition 2. Applications 3. Asymptotics 2. Ensembles
More informationMTH 202 : Probability and Statistics
MTH 202 : Probability and Statistics Lecture 5-8 : 15, 20, 21, 23 January, 2013 Random Variables and their Probability Distributions 3.1 : Random Variables Often while we need to deal with probability
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions A
Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationIEOR 6711: Stochastic Models I SOLUTIONS to the First Midterm Exam, October 7, 2008
IEOR 6711: Stochastic Models I SOLUTIONS to the First Midterm Exam, October 7, 2008 Justify your answers; show your work. 1. A sequence of Events. (10 points) Let {B n : n 1} be a sequence of events in
More informationSolution: The process is a compound Poisson Process with E[N (t)] = λt/p by Wald's equation.
Solutions Stochastic Processes and Simulation II, May 18, 217 Problem 1: Poisson Processes Let {N(t), t } be a homogeneous Poisson Process on (, ) with rate λ. Let {S i, i = 1, 2, } be the points of the
More informationFigure 10.1: Recording when the event E occurs
10 Poisson Processes Let T R be an interval. A family of random variables {X(t) ; t T} is called a continuous time stochastic process. We often consider T = [0, 1] and T = [0, ). As X(t) is a random variable
More informationPhysics in Nuclear Medicine
SIMON R. CHERRY, PH.D. Professor Department of Biomedical Engineering University of California-Davis Davis, California JAMES A. SORENSON, PH.D. Emeritus Professor of Medical Physics University of Wisconsin-Madison
More informationProbability and Stochastic Processes
Probability and Stochastic Processes A Friendly Introduction Electrical and Computer Engineers Third Edition Roy D. Yates Rutgers, The State University of New Jersey David J. Goodman New York University
More informationApplied Probability and Stochastic Processes
Applied Probability and Stochastic Processes In Engineering and Physical Sciences MICHEL K. OCHI University of Florida A Wiley-Interscience Publication JOHN WILEY & SONS New York - Chichester Brisbane
More informationLossless Compression of List-Mode 3D PET Data 1
Lossless Compression of List-Mode 3D PET Data 1 Evren Asma, David W. Shattuck and Richard M. Leahy Signal and Image Processing Institute, University of Southern California, Los Angeles, CA 989-56 Abstract
More informationRADIATION DETECTION AND MEASUREMENT
RADIATION DETECTION AND MEASUREMENT SECOND EDITION GLENN F. KNOLL Professor of Nuclear Engineering The University of Michigan Ann Arbor, Michigan WILEY JOHN WILEY & SONS New York Chichester Brisbane Toronto
More informationAUTOMATED TEMPLATE MATCHING METHOD FOR NMIS AT THE Y-12 NATIONAL SECURITY COMPLEX
AUTOMATED TEMPLATE MATCHING METHOD FOR NMIS AT THE Y-1 NATIONAL SECURITY COMPLEX J. A. Mullens, J. K. Mattingly, L. G. Chiang, R. B. Oberer, J. T. Mihalczo ABSTRACT This paper describes a template matching
More informationS1. Proofs of the Key Properties of CIATA-Ph
Article submitted to INFORMS Journal on Computing; manuscript no. (Please, provide the mansucript number!) S-1 Online Supplement to Modeling and Simulation of Nonstationary and Non-Poisson Processes by
More informationCOMPTON SCATTERING OF GAMMA RAYS
COMPTON SCATTERING OF GAMMA RAYS v2.7 Last revised: R. A. Schumacher, January 2017 I. INTRODUCTION Compton scattering is the name given to the scattering of high-energy gamma rays from electrons. The gamma
More informationA source model for ISDN packet data traffic *
1 A source model for ISDN packet data traffic * Kavitha Chandra and Charles Thompson Center for Advanced Computation University of Massachusetts Lowell, Lowell MA 01854 * Proceedings of the 28th Annual
More informationModeling Recurrent Events in Panel Data Using Mixed Poisson Models
Modeling Recurrent Events in Panel Data Using Mixed Poisson Models V. Savani and A. Zhigljavsky Abstract This paper reviews the applicability of the mixed Poisson process as a model for recurrent events
More informationSupplementary Figures
Supplementary Figures Supplementary Figure. X-ray diffraction pattern of CH 3 NH 3 PbI 3 film. Strong reflections of the () family of planes is characteristics of the preferred orientation of the perovskite
More informationCOUNTING ERRORS AND STATISTICS RCT STUDY GUIDE Identify the five general types of radiation measurement errors.
LEARNING OBJECTIVES: 2.03.01 Identify the five general types of radiation measurement errors. 2.03.02 Describe the effect of each source of error on radiation measurements. 2.03.03 State the two purposes
More informationRandom Processes. DS GA 1002 Probability and Statistics for Data Science.
Random Processes DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Aim Modeling quantities that evolve in time (or space)
More informationEDRP lecture 7. Poisson process. Pawe J. Szab owski
EDRP lecture 7. Poisson process. Pawe J. Szab owski 2007 Counting process Random process fn t ; t 0g is called a counting process, if N t is equal total number of events that have happened up to moment
More informationA Scaled Diæerence Chi-square Test Statistic. Albert Satorra. Universitat Pompeu Fabra. and. Peter M. Bentler. August 3, 1999
A Scaled Diæerence Chi-square Test Statistic for Moment Structure Analysis æ Albert Satorra Universitat Pompeu Fabra and Peter M. Bentler University of California, Los Angeles August 3, 1999 æ Research
More informationThe Statistics of the Number of Neutron Collisions Prior to Absorption
NUCLEAR SCIENCE AND ENGINEERING: 153, 60 68 ~2006! The Statistics of the Number of Neutron Collisions Prior to Absorption Sara A. Pozzi* Oak Ridge National Laboratory, P.O. Box 2008, MS 6010, Oak Ridge,
More informationA SEARCH FOR ASTROPHYSICAL POINT SOURCES OF 100 TEV GAMMA RAYS BY THE UMC COLLABORATION
A SEARCH FOR ASTROPHYSICAL POINT SOURCES OF 100 TEV GAMMA RAYS BY THE UMC COLLABORATION T.A.McKay, b A.Borione, b M.Catanese, a C.E.Covault, b J.W.Cronin, b B.E.Fick, b K.G.Gibbs, b K.D.Green, a S.Hauptfeld,
More informationCRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS
CRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS É. M. Zsolnay Department of Nuclear Techniques, Budapest University of Technology
More informationA new timing model for calculating the intrinsic timing resolution of a scintillator detector
INSTITUTE OF PHYSICS PUBLISHING Phys. Med. Biol. 5 (7) 3 7 PHYSICS IN MEDICINE AND BIOLOGY doi:.88/3-955/5/4/6 A new timing model for calculating the intrinsic timing resolution of a scintillator detector
More informationProblems 5: Continuous Markov process and the diffusion equation
Problems 5: Continuous Markov process and the diffusion equation Roman Belavkin Middlesex University Question Give a definition of Markov stochastic process. What is a continuous Markov process? Answer:
More informationSummer School in Statistics for Astronomers & Physicists June 5-10, 2005
Summer School in Statistics for Astronomers & Physicists June 5-10, 2005 Session on Statistical Inference for Astronomers Poisson Processes and Gaussian Processes Donald Richards Department of Statistics
More informationOn the static assignment to parallel servers
On the static assignment to parallel servers Ger Koole Vrije Universiteit Faculty of Mathematics and Computer Science De Boelelaan 1081a, 1081 HV Amsterdam The Netherlands Email: koole@cs.vu.nl, Url: www.cs.vu.nl/
More information22.56J Noninvasive Imaging in Biology and Medicine Instructor: Prof. Alan Jasanoff Fall 2005, TTh 1-2:30
22.56J Noninvasive Imaging in Biology and Medicine Instructor: Prof. Alan Jasanoff Fall 2005, TTh 1-2:30 Sample problems HW1 1. Look up (e.g. in the CRC Manual of Chemistry and Physics www.hbcpnetbase.com)
More informationLossless Compression of Dynamic PET Data
IEEE TRANSACTIONS ON NUCLEAR SCIENCE, VOL. 50, NO. 1, FEBRUARY 2003 9 Lossless Compression of Dynamic PET Data Evren Asma, Student Member, IEEE, David W. Shattuck, Member, IEEE, and Richard M. Leahy, Senior
More informationSTATISTICAL HANDLING OF RADIOACTIVITY MEASUREMENTS
STATISTICAL HANDLING OF RADIOACTIVITY MEASUREMENTS OBJECTIVES: 1. To learn how to obtain accurate and reliable measurements of radioactivity in samples, notwithstanding the randomness of radioactive decay.
More informationORTEC. Counters/Ratemeters/ Multichannel Scalers. Choosing the Right Counting Solution. The Basic Functions of Counting Systems
Choosing the Right Counting Solution ORTEC offers a variety of instruments to configure counting systems for many applications, from simple to complex. The following descriptions and selection charts will
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationRandom Processes Why we Care
Random Processes Why we Care I Random processes describe signals that change randomly over time. I Compare: deterministic signals can be described by a mathematical expression that describes the signal
More informationAttenuation of Radiation in Matter. Attenuation of gamma particles
Attenuation of Radiation in Matter In this experiment we will examine how radiation decreases in intensity as it passes through a substance. Since radiation interacts with matter, its intensity will decrease
More information11/10/2014. Chapter 1: Introduction to Medical Imaging. Projection (Transmission) vs. Emission Imaging. Emission Imaging
Chapter 1: Introduction to Medical Imaging Overview of Modalities Properties of an Image: Limitations on Information Content Contrast (both object & image): Brightness difference Sharpness (blur): Smallest
More informationWEIBULL RENEWAL PROCESSES
Ann. Inst. Statist. Math. Vol. 46, No. 4, 64-648 (994) WEIBULL RENEWAL PROCESSES NIKOS YANNAROS Department of Statistics, University of Stockholm, S-06 9 Stockholm, Sweden (Received April 9, 993; revised
More informationChemical Engineering 412
Chemical Engineering 412 Introductory Nuclear Engineering Lecture 26 Radiation Detection & Measurement II Spiritual Thought 2 I would not hold the position in the Church I hold today had I not followed
More informationUniform Quadratic Penalties Cause Nonuniform Spatial Resolution
Uniform Quadratic Penalties Cause Nonuniform Spatial Resolution Jeffrey A. Fessler and W. Leslie Rogers 3480 Kresge 111, Box 0552, University of Michigan, Ann Arbor, MI 48109-0552 ABSTRACT This paper examines
More informationPosition Mapping, Energy Calibration, and Flood Correction Improve the Performances of Small Gamma Camera Using PSPMT
Position Mapping, Energy Calibration, and Flood Correction Improve the Performances of Small Gamma Camera Using PSPMT Myung Hwan Jeong, Yong Choi, Member, IEEE, Yong Hyun Chung, Tae Yong Song, Jin Ho Jung,
More informationWooldridge, Introductory Econometrics, 4th ed. Appendix C: Fundamentals of mathematical statistics
Wooldridge, Introductory Econometrics, 4th ed. Appendix C: Fundamentals of mathematical statistics A short review of the principles of mathematical statistics (or, what you should have learned in EC 151).
More informationFinite State Machines. CS 447 Wireless Embedded Systems
Finite State Machines CS 447 Wireless Embedded Systems Outline Discrete systems Finite State Machines Transitions Timing Update functions Determinacy and Receptiveness 1 Discrete Systems Operates in sequence
More informationNuclear Medicine Intro & Physics from Medical Imaging Signals and Systems, Chapter 7, by Prince and Links
Nuclear Medicine Intro & Physics from Medical Imaging Signals and Systems, Chapter 7, by Prince and Links NM - introduction Relies on EMISSION of photons from body (versus transmission of photons through
More informationA8824: Statistics Notes David Weinberg, Astronomy 8824: Statistics Notes 6 Estimating Errors From Data
Where does the error bar go? Astronomy 8824: Statistics otes 6 Estimating Errors From Data Suppose you measure the average depression of flux in a quasar caused by absorption from the Lyman-alpha forest.
More informationSlides 8: Statistical Models in Simulation
Slides 8: Statistical Models in Simulation Purpose and Overview The world the model-builder sees is probabilistic rather than deterministic: Some statistical model might well describe the variations. An
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationProspects for achieving < 100 ps FWHM coincidence resolving time in time-of-flight PET
Prospects for achieving < 100 ps FWHM coincidence resolving time in time-of-flight PET, 28-Feb-2012, ICTR-PHE, Geneva, Switzerland 1 Time-of-flight PET Colon cancer, left upper quadrant peritoneal node
More informationContinuous-time Markov Chains
Continuous-time Markov Chains Gonzalo Mateos Dept. of ECE and Goergen Institute for Data Science University of Rochester gmateosb@ece.rochester.edu http://www.ece.rochester.edu/~gmateosb/ October 23, 2017
More informationPart I Stochastic variables and Markov chains
Part I Stochastic variables and Markov chains Random variables describe the behaviour of a phenomenon independent of any specific sample space Distribution function (cdf, cumulative distribution function)
More informationE. K. A. ADVANCED PHYSICS LABORATORY STATISTICS OF COUNTING WITH A GEIGER COUNTER ARTIFICIAL RADIOACTIVITY
E. K. A. ADVANCED PHYSICS LABORATORY STATISTICS OF COUNTING WITH A GEIGER COUNTER ARTIFICIAL RADIOACTIVITY 1. INTRODUCTION The Geiger Müller (GM tube) detector for ionizing particles operates on the principle
More informationStatistics of Accumulating Signal
Instrument Science Report NICMOS 98-008 Statistics of Accumulating Signal W.B. Sparks, STScI May 1998 ABSTRACT NICMOS detectors accumulate charge which can be non-destructively. The charge accumulation
More informationTHE POISSON WAITING TIME DISTRIBUTION: ANALYSIS AND APPLICATION
THE POISSON WAITING TIME DISTRIBUTION: ANALYSIS AND APPLICATION by Kyle Park ECED 6530 Project Report 12/12/2013 Table of Contents Summary ii Nomenclature ii List of Figures iii 1 The Waiting Time Distribution
More informationDesign of a Lanthanum Bromide Detector for TOF PET
Design of a Lanthanum Bromide Detector for TOF PET A. Kuhn, S. Surti, Member, IEEE, J. S. Karp, Senior Member, IEEE, P. S. Raby, K. S. Shah, A. E. Perkins, Member, IEEE, G. Muehllehner, Fellow Member,
More informationFoundations of Probability and Statistics
Foundations of Probability and Statistics William C. Rinaman Le Moyne College Syracuse, New York Saunders College Publishing Harcourt Brace College Publishers Fort Worth Philadelphia San Diego New York
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
More informationPoisson Processes. Stochastic Processes. Feb UC3M
Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written
More informationECE 565 Notes on Shot Noise, Photocurrent Statistics, and Integrate-and-Dump Receivers M. M. Hayat 3/17/05
ECE 565 Notes on Shot Noise, Photocurrent Statistics, and Integrate-and-Dump Receivers M. M. Hayat 3/17/5 Let P (t) represent time-varying optical power incident on a semiconductor photodetector with quantum
More informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationSimulation of stationary processes. Timo Tiihonen
Simulation of stationary processes Timo Tiihonen 2014 Tactical aspects of simulation Simulation has always a goal. How to organize simulation to reach the goal sufficiently well and without unneccessary
More informationChapter 6: Random Processes 1
Chapter 6: Random Processes 1 Yunghsiang S. Han Graduate Institute of Communication Engineering, National Taipei University Taiwan E-mail: yshan@mail.ntpu.edu.tw 1 Modified from the lecture notes by Prof.
More informationErrata for the ASM Study Manual for Exam P, Fourth Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA
Errata for the ASM Study Manual for Exam P, Fourth Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA (krzysio@krzysio.net) Effective July 5, 3, only the latest edition of this manual will have its
More informationAppendix A. The Particle in a Box: A Demonstration of Quantum Mechanical Principles for a Simple, One-Dimensional, One-Electron Model System
Appendix A The Particle in a Box: A Demonstration of Quantum Mechanical Principles for a Simple, One-Dimensional, One-Electron Model System Real quantum mechanical systems have the tendency to become mathematically
More informationLecture 3. Heat flow. 3.1 Conservation of energy
Lecture 3 Heat flow 3.1 Conservation of energy One of the recuring themes in this class and in all of physics is that of conservation of energy. You have all heard this before: energy can be neither created
More information