Investigation of the Background γ Radiation and Poisson Statistics
Jason Gross MIT - Department of Physics Jason Gross (8.13) Poisson Statistics October 7, 2011 1 / 26Investigation of the Background Radiation and Poisson Statistics - - PowerPoint PPT Presentation
Investigation of the Background Radiation and Poisson Statistics - - PowerPoint PPT Presentation
Investigation of the Background Radiation and Poisson Statistics Jason Gross MIT - Department of Physics Jason Gross (8.13) Poisson Statistics October 7, 2011 1 / 26 Random Variables Question: What does it mean for something to be
Random Variables
Question: What does it mean for something to be “random”?
Jason Gross (8.13) Poisson Statistics October 7, 2011 2 / 26Random Variables
XKCD’s Answer:
Jason Gross (8.13) Poisson Statistics October 7, 2011 3 / 26Random Variables
Wikipedia’s Answer: A random variable is a function from a probability space to a measurable space, typically the real numbers.
Jason Gross (8.13) Poisson Statistics October 7, 2011 4 / 26Random Variables
Wikipedia’s Answer: A random variable is a function from a probability space to a measurable space, typically the real numbers.
Jason Gross (8.13) Poisson Statistics October 7, 2011 4 / 26Random Variables
Wikipedia’s Answer: A random variable is a function from a probability space to a measurable space, typically the real numbers.
Jason Gross (8.13) Poisson Statistics October 7, 2011 4 / 26Random Variables
Wikipedia’s Answer: A random variable is a function from a probability space to a measurable space, typically the real numbers.
Jason Gross (8.13) Poisson Statistics October 7, 2011 4 / 26Random Variables
Wikipedia’s Answer: A random variable is a function from a probability space to a measurable space, typically the real numbers.
Jason Gross (8.13) Poisson Statistics October 7, 2011 4 / 26- f:
Random Variables
Yikes!
Jason Gross (8.13) Poisson Statistics October 7, 2011 6 / 26Random Variables
Let’s try again.
Jason Gross (8.13) Poisson Statistics October 7, 2011 7 / 26Random Variables
Question: What does it mean for something to be “random”? Answer:
Jason Gross (8.13) Poisson Statistics October 7, 2011 7 / 26Random Variables
Question: What does it mean for something to be “random”? Answer:
Jason Gross (8.13) Poisson Statistics October 7, 2011 7 / 26Random Variables
Question: What does it mean for something to be “random”? Answer: Why do we care?
Jason Gross (8.13) Poisson Statistics October 7, 2011 7 / 26Random Variables
Question: Why do we care?
Jason Gross (8.13) Poisson Statistics October 7, 2011 8 / 26Random Variables
Question: Why do we care? Answer: We want to be able to talk quantitatively about the relationship between measurements and theory.
Jason Gross (8.13) Poisson Statistics October 7, 2011 8 / 26Random Variables
Question: What does it mean for something to be “random”?
Jason Gross (8.13) Poisson Statistics October 7, 2011 9 / 26Random Variables
Law of Large Numbers: The average of a series of identical
statistically independent random
variables almost always converges to the true mean.
Jason Gross (8.13) Poisson Statistics October 7, 2011 10 / 26Random Variables
Law of Large Numbers: The average of a series of identical
statistically independent random
variables almost always converges to the true mean.
Jason Gross (8.13) Poisson Statistics October 7, 2011 10 / 26Random Variables
Law of Large Numbers: lim
n→∞ P
- |Xn − µ| > ε
- = 0
(weak law) P
- lim
n→∞ Xn = µ
- = 1 (strong
law)
Jason Gross (8.13) Poisson Statistics October 7, 2011 10 / 26Random Variables
Law of Large Numbers: lim
n→∞ P
- |Xn − µ| > ε
- = 0
(weak law) P
- lim
n→∞ Xn = µ
- = 1 (strong
law)
Jason Gross (8.13) Poisson Statistics October 7, 2011 10 / 26Random Variables
Law of Large Numbers: The average of a series of identical
statistically independent random
variables almost always converges to the true mean.
Jason Gross (8.13) Poisson Statistics October 7, 2011 10 / 26Random Variables
Question: Why is the law of large numbers enough? Answer:
Jason Gross (8.13) Poisson Statistics October 7, 2011 11 / 26Random Variables
Question: Why is the law of large numbers enough? Answer: Because everything is a random variable!
Jason Gross (8.13) Poisson Statistics October 7, 2011 11 / 26Binomial Distribution
Consider an event with two
- utcomes.
Binomial Distribution
P(x; n, p) = n x
- pxqn−x
= n! x!(n − x)!pxqn−x
Jason Gross (8.13) Poisson Statistics October 7, 2011 13 / 26Binomial Distribution
P(x; n, p) = n x
- pxqn−x
= n! x!(n − x)!pxqn−x
Jason Gross (8.13) Poisson Statistics October 7, 2011 13 / 26Binomial Distribution
P(x; n, p) = n x
- pxqn−x
= n! x!(n − x)!pxqn−x
Jason Gross (8.13) Poisson Statistics October 7, 2011 13 / 26Binomial Distribution
P(x; n, p) = n x
- pxqn−x
= n! x!(n − x)!pxqn−x
Jason Gross (8.13) Poisson Statistics October 7, 2011 13 / 26Binomial Distribution
P(x; n, p) = n x
- pxqn−x
= n! x!(n − x)!pxqn−x
Jason Gross (8.13) Poisson Statistics October 7, 2011 13 / 26Binomial Distribution
n 10 p 0.3 2 4 6 8 10 x 0.05 0.10 0.15 0.20 0.25 PxPx n x px 1 pnx
Jason Gross (8.13) Poisson Statistics October 7, 2011 14 / 26Poisson Distribution
Consider the limit of p ≪ 1, with µ = np fixed. P(x; µ) = µx x! e−µ
Jason Gross (8.13) Poisson Statistics October 7, 2011 15 / 26Poisson Distribution
Consider the limit of p ≪ 1, with µ = np fixed. P(x; µ) = µx x! e−µ
Jason Gross (8.13) Poisson Statistics October 7, 2011 15 / 26Poisson Distribution
Μ 3 2 4 6 8 10 x 0.05 0.10 0.15 0.20 PxPx Μx x Μ
Jason Gross (8.13) Poisson Statistics October 7, 2011 16 / 26Poisson Distribution
Question: Why do we care about this limit? Answer:
Jason Gross (8.13) Poisson Statistics October 7, 2011 17 / 26Poisson Distribution
Question: Why do we care about this limit? Answer: It’s everywhere! (when we count statistically independent events)
Jason Gross (8.13) Poisson Statistics October 7, 2011 17 / 26Gaussian Distribution
Consider the limit for µ large. P(x; µ) = 1 µ √ 2π e−1
2
x−µ
µ
2
Jason Gross (8.13) Poisson Statistics October 7, 2011 18 / 26Gaussian Distribution
Consider the limit for µ large. P(x; µ) = 1 µ √ 2π e−1
2
x−µ
µ
2
Jason Gross (8.13) Poisson Statistics October 7, 2011 18 / 26Gaussian Distribution
Μ 100 Σ 10 80 100 120 140 x 0.01 0.02 0.03 0.04 PxPx 1 2 Π Σ
- 1
Methodology
+HV Power Supply Scintillator NaI Photo Multiplier Tube Voltage Divider Preamplifier Amplifier Counter Jason Gross (8.13) Poisson Statistics October 7, 2011 20 / 26Methodology
Jason Gross (8.13) Poisson Statistics October 7, 2011 20 / 26 +HV Power Supply Scintillator NaI Photo Multiplier Tube Voltage Divider Preamplifier Amplifier CounterMethodology
Jason Gross (8.13) Poisson Statistics October 7, 2011 20 / 26 Scintillator NaI Photo Multiplier Tube Voltage ⑦ lightMethodology
Jason Gross (8.13) Poisson Statistics October 7, 2011 20 / 26 Scintillator NaI Photo Multiplier Tube Voltage ⑥ ③ e−Results
Number of Counts 1 2 3 4 5 6 10 20 30 40 50 Μ 1 s x 0.96 0.98 Χ2 PValue Curve Poisson 3.4 0.33 Normal 3.5102 5.51067 Frequency Bin Jason Gross (8.13) Poisson Statistics October 7, 2011 21 / 26Results
Number of Counts 2 4 6 8 10 12 14 5 10 15 20 25 Μ 4 s x 5.5 2.3 Χ2 PValue Curve Poisson 6.0 0.65 Normal- 84. 6.81013
Results
Number of Counts 5 10 15 20 5 10 15 20 Μ 10 s x 9.9 3.1 Χ2 PValue Curve Poisson 5.7 0.77 Normal- 28. 0.0052
Results
Number of Counts 80 100 120 140 160 2 4 6 8 10 12 14 Μ 100 s x 115. 11. Χ2 PValue Curve Poisson 11. 0.55 Normal- 10. 0.59
Results
Poisson Normal µ χ2 P-value χ2 P-value µ ≈ 1 3.4 0.32 350 5.5 · 10−67 µ ≈ 4 6.0 0.65 84 6.8 · 10−13 µ ≈ 10 5.7 0.77 28 0.0052 µ ≈ 100 11 0.55 10. 0.59
Jason Gross (8.13) Poisson Statistics October 7, 2011 25 / 26Thank You! Any questions?
Jason Gross (8.13) Poisson Statistics October 7, 2011 26 / 26