.. index:: MuH
Fitting function for use by Muon scientists defined by:
A(t)=\frac{A_0(t)}{6}e^{-\lambda t}e^{-\frac{(\sigma t)^2}{2}}\left(1+\cos(\omega_{D}t + \phi)+2\cos\left(\frac{\omega_D}{2}t+\phi\right)+2\cos\left(\frac{3\omega_D}{2}t+\phi\right)\right)
A_0 is the amplitude,
\lambda is the exponential decay rate,
\sigma is the gaussian decay rate,
\omega_D = 2 \pi \nu_D where \nu_D (MHz) is the oscillating frequency,
and \phi (rad) is the phase at time t=0.
.. plot::
from mantid.simpleapi import FunctionWrapper
import matplotlib.pyplot as plt
import numpy as np
x = np.arange(0.1,16,0.1)
y = FunctionWrapper("MuH", A0 = 0.5, NuD = 0.5, Lambda = 0.3, Sigma = 0.05, Phi = 0.0)
fig, ax=plt.subplots()
ax.plot(x, y(x))
ax.set_xlabel('t($\mu$s)')
ax.set_ylabel('A(t)')
.. attributes::
.. properties::
[1] T. Lancaster et al., J. Phys.: Condens. Matter 21 346004 (2009).
.. categories::
.. sourcelink::