Physics307L:People/Gooden/Notebook/071015

POISSON STATISTICS

Objective

 * Using seemingly random data, most likely of cosmic origin, we hope to show using the Multichannel Analyser (MCA) that the data best fits a Poisson distribution rather than the Binomial or Gaussian Distributions. We are taking data over various apertures of time; for instance, 256 bins of 100 milli-seconds each. The large amount of incidents being recorded over varying sized bins will give us a large variety of distributions.

Experiment
We have a setup that consists of a photomultiplier tube that is attached to a NaI scintillator, both are housed in a structure of lead bricks. The arrangement is wired to a high voltage power supply (1000 volts) and then run through some sort of bridge and to a data acquisition board on a computer. The computer is running a program called PCAIII, which handles the data acquisition process. The photomultiplier tube and scintillator were connected by way of coaxial cables to the power supply which we connected to the bridge, and from the bridge into the computer. There were some erroneous cables coming from the data acquisition board that we had no need to mess with.
 * Setup

Once all the cables are secured and the power supply and bridge are warmed up, we were able to start taking data. Using PCAIII we simply configured how many bins of data we were taking and how much time each bin would get ("dwell time").You want to vary your bins, in both number and in dwell time. For example we chose 512 bins(for 800s, 2s, and 10s) and 4096 bins (for 10ms, 100ms, and 10s).
 * Procedure

Theory
When collecting large amounts of data it is wise to look at the probability distributions for that data.The binomial distribution is the most true distribution and from it we can derive the Gaussian and Poisson distributions as limiting cases. In this experiment we are concerning ourselves with the Poisson Distribution and seeing how closely our random data fit to it.

The Binomial Distribution
When analyzing any randomly distributed situation a binomial distribution:

$$B(x)=\frac{N!}{x!(N-x)!}p^nq^{N-n}$$

with a standard deviation of

$$\sigma=\sqrt{pN(1-p)}$$

and a mean of

$$a=pN$$

is used. With $$N$$ = the number of counts, $$p$$ = the probability of counts occurring, and $$q$$ = the probability of counts not occurring. In all instances $$p+q=1$$, since something either happens or it doesn't, $$p$$ and $$q$$ must sum to 1. In context of our experiment, we have a very large $$N$$ with a very small $$p$$. Undergoing several manipulations we can approximate the binomial distribution to be the Poisson distribution. More information can be found here

The Gaussian Distribution
When analyzing a situation in which there is a high probability of occurrence (large $$p$$) we use the Gaussian (or normal) distribution. The Gaussian distribution is given by

$$G(x)=\frac{1}{\sqrt{2\pi\sigma^2}}e^{-\frac{\left(x-a\right)^2}{2\sigma^2}}$$,

with $$a$$ = the mean, $$\sigma$$ = the standard deviation. The Gaussian distribution is often used to model probabilities and is useful because if the standard deviation and mean are optimal then the actual mean and standard deviation values will match those given theoretically. A very good tool for understanding the Gaussian distribution can be found here

The Poisson Distribution
The Poisson is a "discrete probability distribution" for the probability of a number of events occuring in a fixed period of time, such that the events occur at a known average rate. When analyzing a random situation in which there is a very low probability of occurrence (large $$N$$ and small $$p$$)we use the Poisson distribution. The standard form is given by:

$$P(x)=e^{-a}\frac{a^x}{x!}$$

with a standard deviation of

$$\sigma=\sqrt{a}$$,

with $$a$$ = the mean. The Poisson distribution appears only around zero and, unlike the Gaussian or binomial distributions, can only reflect positive integers. You could think of a Gaussian distribution that has been normalized so it can only take on positive values with a mean greater than zero. A place to get comfortable with Poisson distributions can be accessed here 1

Data and Results
We took four sets of data corresponding to different combinations of bin number and time delay. The time delay represented for how long measurements/counts were taken from each bin. The number of bins represents basically how many mini experiments were done for each data set. Using the MCA we were able to repeat the experiment of making counts of background radiation many many times with relative ease. The four sets of data that we took can be broken down such that, for data set 1 we used 512 bins and a time delay of 800ms, data set 2 512 bins and 100ms, data set 3 256 bins and 10s, and finaly for data set four we used 4096 bins with a time delay of 40s which took approximately 45 hours to run.

This is Data set 1, where mentioned above we used 512 bins and a time delay of 800ms. For this data set we found a mean value of $$a=10.53$$ and standard deviation of $$\sigma=3.55$$

This is Data set 2, where mentioned above we used 512 bins and a time delay of 100ms. For this data set we found a mean value of $$a=1.36$$ and standard deviation of $$\sigma=1.37$$

This is Data set 3, where mentioned above we used 256 bins and a time delay of 10s. For this data set we found a mean value of $$a=135.72$$ and standard deviation of $$\sigma=14.25$$

This is Data set 4, where mentioned above we used 4096 bins and a time delay of 40s. For this data set we found a mean value of $$a=543.60$$ and standard deviation of $$\sigma=34.63$$

Plot 2:This is a plot of the data, which was transformed into a probability by counting how many times it occured in the data set and dividing by the number of bins. This plot represents data set 2 for 100ms and plotted in green is the poisson, red is the data and in blue is the gaussian. Plot 3:This is a plot of the data, which was transformed into a probability by counting how many times it occured in the data set and dividing by the number of bins. This plot represents data set 3 for 10s and plotted in green is the poisson, red is the data and in blue is the gaussian. Plot 4:This is a plot of the data, which was transformed into a probability by counting how many times it occured in the data set and dividing by the number of bins. This plot represents data set 4 for 40s and plotted in green is the poisson, red is the data and in blue is the gaussian.
 * Plot 1:This is a plot of the data, which was transformed into a probability by counting how many times it occured in the data set and dividing by the number of bins. This plot represents data set 1 for 800ms and plotted in green is the poisson, red is the data and in blue is the gaussian.

Discussion
This experiment used a detector that picked up background radiation and by making counts of that background over selected periods of time many many times we produced the data given above. We were trying to compare the data taken, which in principle should be completely random, with the poisson distribution. As seen in plot 1 above,the data that we took fits well with the Poisson and also the Gaussian distribution. The poisson and gaussian do not differ by much with this data set. From plot 2 we can see that while most of the data that I plotted is not perfectly on the Poisson distrubtion is is much closer to fitting that distribution than the Gaussian. Plot 3 has the Poisson and Gaussian distributions being almost indistinguishable from eachother, while on the other hand the data for that set is very random and does not fit either distribution at all. I am not sure why there is so much error in that data set which is its difference from the expected probabilities. Then for plot 4, we can see that the data closely matches the plot of the Gaussian Distribution, however if you notice there is not plot of the Poisson. The reason for this is that because the counts were so large, ~600 or so, and because the Poisson uses of the counts and a factorial of them, Matlab was unable to compute the number and produced only infinities and NaN, which stands for Not a Number.