Download LAB 5

Richard Kass
Physics 416
Fall 2005
LAB 5
The goal of this lab is to become familiar with the concept of energy resolution and to measure
the energy resolution of your NaI detector. In a previous lab, Lab 4, you measured the response
of the NaI detector to gamma rays and saw that a gamma ray shows up as a "peak" or "bump" in
the MCA (Multi Channel Analyzer) spectrum. The shape of the peak is to a very good
approximation described by a Gaussian function. The mean of this Gaussian corresponds to the
energy of the gamma ray while the width of the Gaussian tells us something about how precisely
the NaI detector measures the energy of the gamma ray.
The standard definition of energy resolution (R) is:
Full WidthHalf Maximum
Position of Peak
For a Gaussian distribution the position of the peak = mean, and Full Width Half Maximum
(FWHM) is related to the standard deviation () by FWHM=2.36.
R
(I) To get a feeling for the spread in gamma ray energy from a relatively long lived isotope
estimate the natural line width (E) of Cs137 using the Heisenberg Uncertainty Principle:
Et 
16
with  6.58 10 eV  sec , t taken to be the half-life of Cs137, 30 years.
Using your energy calibration curve from Lab 4 roughly how many voltage bins does this energy
spread (E) cover? Could you possibly ever see this energy spread using the equipment that we
have in this lab?
(II) Measure the energy resolution of your NaI detector using Co60, Na22, and Cs137. Take a
spectrum of each sample and save the data run to disk. Make sure that you accumulate enough
data for each sample by waiting (5-10 minutes) until the spectrum has "smoothed" itself out.
(a) Crude calculation of energy resolution:
Calculate the energy resolution using the voltages (x-axis values) from your MCA spectrum. For
each of the five peaks (2 from Na22, 2 from Co60, and 1 from Cs137) estimate the FWHM and
location of the peak in terms of the voltage. Calculate the energy resolution (R) for the five
gamma rays and make a plot of your energy resolution (y-axis) vs. the true peak energy (x-axis).
(b) A better calculation of the energy resolution:
WARNING: You will need to use Kaleidagraph for this section!!
Fit your data to a function that represents the measured gamma ray spectrum. The bins near the
gamma ray the spectrum can be thought to consist of two components, the gamma ray itself (our
signal), which we can model with a Gaussian, and background. If there were no background then
the estimated number of counts in bin i would be NiS :
S
i
N 
A

e

( xi   ) 2
2
2
Here xi is energy at the center of bin i. Use your energy calibration curve from Lab 4 to convert
your x-axis values (volts) into energy (MeV). We would like to determine the constants A, , and
 that best describe our data. Unfortunately, we also have background in our spectrum. Let's
assume that NiB , the number of background counts in bin i, (i.e. what we would have if we took
the radioactive source away) in the area of interest can be modeled by a straight line (is this
reasonable?):
B
Ni  C  Dxi
Again, the constants C and D are to be determined from our data. Now, since our data really
consists of background and signal the best representation of our data is given by the sum of the
signal and background functions:
(x   ) 2
A  i2  2
S B
Ni  e
 C  Dx i

The problem now is to determine the 5 constants (A, C, D, , ) from our data points. In
principle this is simple, we use a technique such as Maximum Likelihood or Least Squares to
find the 5 parameters. In practice this is difficult since the above function is non-linear in some of
the variables, thus we would have to solve a set of 5 coupled non-linear equations! Rather than
write a computer program to do this (you can if you want), we will use a feature of Kaleidagraph
that fits complicated functions to data. Review the handout “A Quick Guide to Kaleidagraph” or
read Chapter 9 of the Kaleidagraph Learning Guide to learn how to fit your data to a function
and then fit your 5 gamma-rays to the above function.
Note: Kaleidagraph will try to fit all of the data points in a file to the function of interest. In most cases however we only want to
fit a small part of the spectrum (e.g. 50 bins out of 250) to the function. Therefore it is best to copy the data in the region of
interest to a new file and work with this file while performing the fits.
For each of the five peaks calculate the energy resolution using the mean () and standard
deviation () obtained from your fit. Don't forget to convert standard deviation into FWHM (you
could fit for the FWHM instead of )! Plot the resolutions from part (b) along with the
resolutions from parts (a). Is the energy resolution constant as a function of energy?
Extra Credit:
Redo your Kaleidagraph fits using the optional weighting algorithm. Here you check the Weight
Data option under DEFINE when using Curve Fit. You should weight each data point using its
statistical error. What are the 2‘s of the fits? Are they reasonable?
Super Extra Credit:
Simulate, by writing a computer program, the performance of a NaI detector. In this part of the
lab we would like to get a feeling for what a Co60 gamma ray spectrum would look like if the
NaI detector had a percent energy resolution (100R) of 3%, 6%, or 25%. To do this you will
write a program that simulates the response of a detector with a known energy resolution. In this
exercise assume that the response of the detector is given by a Gaussian function with known
mean () and standard deviation (). When Co60 decays it gives off two distinct gamma rays
(not always, but 99% of the time), a lower energy gamma with E1 = 1.172 MeV and a higher
energy gamma ray with E2 = 1.333 MeV. Let's represent the response of the NaI detector as a
Gaussian with = E1 or E2 depending on which gamma ray we are detecting. Using the
definition of energy resolution (R), and defining E as the energy of the gamma ray the of the
Gaussian is given by:
E R

2.36
From a previous lab we know how to generate a Gaussian distribution using random numbers.
Recall that the Gaussian you generated had = 0 and = 1. However for this lab we want a
Gaussian with = E1 or E2 and  given by the energy resolution. To transform (apart from a
normalization constant) from a set of numbers {g1...gn} that are distributed according to a
Gaussian distribution with = 0 and = 1 to a set of numbers {G1...Gn} also Gaussian
distributed but with = m and = s use:
Gi  s  gi  m
Using the definition of a Gaussian distribution show that the above transformation really
does what it claims to do!
Use the following prescription to simulate on the computer the energy response of the detector to
Co60. The goal here is to make a histogram similar to the one displayed by the LABVIEW
program MCA when you measured the CO60 spectrum.
a) Choose an energy resolution (first do 3%, then 6%, then 25%).
b) Use your Gaussian number generator to obtain a set of gi's. (you need 20,000 g’s)
c) Transform each gi to the new Gaussian variable Gi. Remember s depends on the
chosen energy resolution and m is the true energy of the gamma ray. The G's now
represent the measured energy of the gamma rays. For each g you generate you need to
decide if represents the low energy or high energy gamma ray. For the low energy
gamma m = 1.172, for the high energy gamma m = 1.333 MeV.
d) Keep track of the number of measured gamma rays with energies in the interval
[E, E+E]. The LABVIEW program uses a bin size of E ~ 15 keV. I suggest using 25
bins with this bin size, the first bin starting at E = 1.1 MeV.
For each energy resolution generate 104 low energy gamma rays (E1) and 104 high energy gamma
rays (E1), this represents 104 decays of Co60. Make a histogram of your energy spectrum for
each of the three energy resolutions. Which of the three histograms looks most like your actual
Co60 spectrum? Comment on how easy/difficult it is to tell that Co60 decays to two gamma rays
if the energy resolution of your detector is ~25%.