\(\renewcommand\AA{\unicode{x212B}}\)

# Error Propagation¶

The purpose of this document is to explain how Mantid deals with error propagation and how it is used in its algorithms.

## Theory¶

In order to deal with error propagation, Mantid treats errors as Gaussian probabilities (also known as a bell curve or normal probabilities) and each observation as independent. Meaning that if \(X = 100 \pm 1\) then it is still possible for a value of \(102\) to occur, but less likely than \(101\) or \(99\), and a value of \(105\) is far less likely still than any of these values.

## Plus and Minus Algorithm¶

The Plus v1 algorithm adds two datasets together, propagating the uncertainties. Mantid calculates the result of \(X_1 + X_2\) as

\(X = X_1 + X_2\)

with uncertainty

\(\sigma_X = \sqrt{ \left( \sigma_{X_1} \right)^2 + \left( \sigma_{X_2} \right)^2 }\).

Consider the example where \(X_1 = 101 \pm 2\) and \(X_2 = 99 \pm 2\). Then for this algorithm:

\(X = X_1 + X_2 = 101 + 99 = 200\)

\(\sigma_X = \sqrt{ 2^2 + 2^2} = \sqrt{8} = 2.8284\)

Hence the result of Plus v1 can be summarised as \(X = 200 \pm \sqrt{8}\).

Mantid deals with the Minus v1 algorithm similarly: the result of \(X_1 - X_2\) is

\(X = X_1 - X_2\)

with error

\(\sigma_X = \sqrt{ \left( \sigma_{X_1} \right)^2 + \left( \sigma_{X_2} \right)^2 }\).

## Multiply and Divide Algorithm¶

The Multiply v1 and Divide v1 algorithms propagate the uncertainties according to (see also here):

\(\sigma_X = \left|X\right| \sqrt{ \left( \frac{\sigma_{X_1}}{X_1} \right)^2 + \left( \frac{\sigma_{X_2}}{X_2} \right)^2 }\),

where \(X\) is the result of the multiplication, \(X = X_1 \cdot X_2\), or the division, \(X = X_1 / X_2\).

Considering the example above where \(X_1 = 101 \pm 2\) and \(X_2 = 99 \pm 2\). Mantid would calculate the result of \(X_1 / X_2\) as \(X = 101 / 99 = 1.0202\), with uncertainty \(\sigma_X = 1.0202 \sqrt{ \left(2/101\right)^2 + \left(2/99\right)^2} = 0.0288\).

For Multiply v1, the result of \(X_1 \times X_2\) is \(X = 101 \times 99 = 9999\), with uncertainty \(\sigma_X = 9999 \sqrt{ \left(2/101\right)^2 + \left(2/99\right)^2} = 282.8568\).

**Category**: Concepts