textbook:nrctextbook:chapter14

This is an old revision of the document!


14. Statistical uncertainties in radioactivity measurements

Chapter 14 from BASICS OF NUCLEAR PHYSICS AND OF RADIATION DETECTION AND MEASUREMENT - An open-access textbook for nuclear and radiochemistry students by Jukka Lehto

14.1. Count rate – activity

When measuring the activity ($\text{A}_x$) of a radioactive source the primary result is the total (gross) count rate ($\text{R}_g$) obtained from the measurement system (detector, amplifier and pulse counter)

$$R_g = \frac{X_g}{t}$$

[XIV.I]

where $\text{X}_g$ = number of collected total pulses and $\text{t}$ = measurement time. The unit of count rate is pulses per unit time: counts per second (s-1, cps) or counts per minute (cpm).

The observed gross count rate ($\text{R}_g$) includes, in addition to pulses resulting from the radioactive source (net pulses $\text{X}_n$), also pulses from background ($\text{X}_{bg}$) originating from various sources other than the actual source, such as cosmic radiation, presence of natural or pollution radionuclides and electric noise of the measurement system. These background pulses need to be counted separately in the absence of the radioactive source and the background count rate must be subtracted from the gross count rate to obtain the net count rate ($\text{R}_n$) originating from the radioactive source.

$$R_n = \frac{X_g}{t_g} - \frac{X_{bg}}{t_{bg}}$$

[XIV.II]

Activity of the source ($\text{A}_{x}$) is calculated either by comparing the net count rate of the source ($\text{R}_{x}$) to that obtained by measuring a standard source ($\text{R}_{st}$) with a known activity ($\text{A}_{st}$) in identical conditions as the unknown source

$$A_x = A_{st} \times \left( \frac{R_x}{R_{st}} \right)$$

[XIV.III]

or if the counting efficiency ($\text{E(%)}$) of the counting system is known dividing the net count rate with the counting efficiency.

$$A_x = R_x \left( \frac{E}{100} \right)$$

[XIV.IV]

What kind of uncertainties are involved here and how they are calculated are discussed below.

14.1. Systematic and random errors

In every measurement, including radioactivity measurement, there are two types of errors resulting in an uncertainty in the measurement result:

  1. Systematic errors arise from erroneous measurement system and they always function into same direction from the right result. If, for example, the activity of the standard is not what it is supposed to be or the settings of the measurement system, such as amplification of pulses, change during the measurement, this causes error to the observed activity. Even if the parallel results were close to each other, i.e. reproducible and precise, the results would not be accurate since they systematically deviate from the real value to certain direction in case of a systematic error.
  2. Random errors arise from the fact that the measurement system or the phenomenon measured, orboth, are intrinsically non-deterministic (stochastic):
  • The measurement systems are always non-ideal and do not always give the same response even though the measured quantity would have a constant value. For example, alpha particles for a certain transition of a radionuclide have always the same energy but we never obtain a perfect line peak spectrum, but the peak has a broadness depending on the system's limited preciseness in transforming alpha particle energies to electric pulses.
  • Some measured phenomena, such as radioactivity, are intrinsically stochastic. We can knowthe probability of nuclear transformations, decays, at certain time interval but it is impossible to find out their exact number since they vary in a stochastic manner.An important feature of a stochastic error is that the reproducibility and preciseness increases with the number of measurements.

Below we discuss in more detail the uncertainties arising from the stochastic nature of radioactive decay.

Figure XIV.1. Effect of systematic and random error on observed results. Left side: high precision but low accuracy. Right side: low precision but high accuracy.

14.2. Poisson and normal distribution

The variation of radioactive decay events and other stochastic processes with low and constant probabilities are mathematically described with Poisson distribution probability function (Equation XIV.V). The variation of radioactive decays (or particle/photon flux) is a fundamental physical characteristic of the radionuclide. If we consider a large enough number of radionuclides the number of decay rate varies with time following the equation:

$$P_x = \frac{m^x}{x!} \cdot e^{-m}$$

[XIV.V]

where $\text{P}_x$ is the probability for $\text{x}$ number of events occurring in unit time and $\text{m}$ is the most probable number of events. Since the number of decay events can have only integer values the graphical representation of Poisson distribution is a histogram. Poisson distribution is also not symmetrical but is slightly bended to lower values. To make treatment of results simpler the symmetric normal (Gaussian) distribution is used as an approximation to Poisson distribution. Figure XIV.2 shows the difference of Poisson and normal distributions for the probability to observe decay events in a number of identical time intervals. For a large number (<30) of events (m) Poisson and normal distributions are more or less identical.

Figure XIV.2. Poisson and normal distribution functions for the probability (P) to observe radioactive decay events (m) in a number of identical time intervals (http://nau.edu/cefns/labs/electron-microprobe/glg-510-class-notes/statistics/).

The mathematical formulation of normal distribution is as follows:

$$P_x = \frac{1}{{\sigma \sqrt{2\pi}}} \times e^{-\frac{{(x - m)^2}}{{2\sigma^2}}}$$

[XIV.VI]

where $\text{P}_x$ is the appearance probability of a stochastic event, $\text{m}$ is the real value of the events, $\sigma$ is the standard deviation of events at various time intervals.

14.3. Standard deviation

### To present the variation for a number of decay events/pulses the quantity used is standard deviation, s.Its mathematical expression is

$$\sigma_x = \sqrt{\frac{\sum (x_i - \bar{x})^2}{n-1}}$$

textbook/nrctextbook/chapter14.1737593726.txt.gz · Last modified: 2025-01-23 01:55 by Merja Herzig