General statistics notes on counting experiments

From Vacuum Ultra Violet
Revision as of 14:22, 15 December 2022 by Mnuland@nikhef.nl (talk | contribs) (Create page)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Since VULCAN will measure light signals, whether they are reflected, fluoresced or other, basically all measurements will look the same. In the simplest case we look at the signal from a single SiPM, and we want to calculate the rate. This is a simple counting experiment.

Counting experiments

In a counting experiment, you count the number of occurances of an event in a certain time interval. Currently, we record 16us waveforms, but the exact length of the time interval is irrelevant. In these 16us, we want to know how much light was seen. So, we count the number of photo electron (p.e.) peaks. These peaks have a very distinct shape and a well-defined height. It depends on the measurement which peaks are relevant, e.g. for dark count measurements only the single p.e. peaks are interesting. In the end, we'll have a number of events in a certain time interval. These counts are Poisson distributed. So, if we take multiple measurements, and our process it truly a Poisson process, we'll be able to fit a nice Poisson distribution to the counts we measured. The maximum likelihood estimator (MLE) for the Poisson parameter λ is then simply the sample mean. This MLE is unbiased, efficient, complete and sufficient. λ is the expected value, as well as the variance, so when we have λ, we are basically done.

The rates we expected are rather high, especially at room temperature. For example, the VUV sensitive SiPMs have a dark count rate of ~1.5-2 MHz at room temperature. Since we can see a Poisson(100) distribution as the sum of 100 Poisson(1) distributions, we can apply the central limit theorem and approximate the Poisson distribution as a Gaussian if we have more than ~20 counts per time interval (or above 1.25 MHz). We'll need to apply a continuity correction too. For lower rates, we can simply take longer measurements in order to still be able to use the Gaussian approximation, or use Poisson distributions.