/
CompositePCRmagnet.py
36 lines (29 loc) · 1.41 KB
/
CompositePCRmagnet.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
# Mantid Repository : https://github.com/mantidproject/mantid
#
# Copyright © 2019 ISIS Rutherford Appleton Laboratory UKRI,
# NScD Oak Ridge National Laboratory, European Spallation Source,
# Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
# SPDX - License - Identifier: GPL - 3.0 +
# pylint: disable=invalid-name, anomalous-backslash-in-string, attribute-defined-outside-init
from mantid.api import IFunction1D, FunctionFactory
import numpy as np
class CompositePCRmagnet(IFunction1D):
def category(self):
return "Muon\\MuonSpecific"
def init(self):
self.declareParameter("A0", 0.2, "Amplitude")
self.declareParameter("KTdelta", 0.1, "KT Relaxation rate")
self.declareParameter("B0", 10, "Local Magnetic Field (G)")
self.declareParameter("Gauss", 0.2, "Gaussian relaxation rate")
def function1D(self, x):
A0 = self.getParameterValue("A0")
KTdelta = self.getParameterValue("KTdelta")
B0 = self.getParameterValue("B0")
Gaus = self.getParameterValue("Gauss")
gmu = 0.01355342 * 2 * np.pi
if B0 > 0:
output = 1.0 / 3 + 2.0 / 3 * np.cos(gmu * B0 * x) * np.exp(-((Gaus * x) ** 2) + 2.0 / 5 * (KTdelta * x) ** 2)
else:
output = 1.0 / 3 + 2.0 / 3 * (1 - (KTdelta * x) ** 2) * np.exp(-((KTdelta * x) ** 2) / 2)
return A0 * output
FunctionFactory.subscribe(CompositePCRmagnet)