def intA(p=p, extra_par={}):
old_p = p.copy()
p.update(extra_par)
wJ = sqrt((p["mJ"]) / p["DJ"])
wA = sqrt((p["mA"]) / p["DA"])
result = (
sinh(wA * p["d"])
/ wA
/ cosh(wA * p["d"])
* (
1 / (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"]))
- 1.0 / (p["r"] * p["d"] * exp(-p["a"] * p["L1"]))
)
)
p.update(old_p)
return result
def J(x, p=p, extra_par={}):
old_p = p.copy()
p.update(extra_par)
assert x >= 0
assert x <= p["L1"]
wJ = sqrt((p["mJ"]) / p["DJ"])
wA = sqrt((p["mA"]) / p["DA"])
result = p["K"] * exp(wJ * (x - p["L1"])) / 2.0 * (1 - wA / wJ * tanh(wA * p["d"])) * (
1 / (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"]))
- 1.0 / (p["r"] * p["d"] * exp(-p["a"] * p["L1"]))
) + p["K"] * exp(-wJ * (x - p["L1"])) / 2.0 * (1 + wA / wJ * tanh(wA * p["d"])) * (
1 / (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"]))
- 1.0 / (p["r"] * p["d"] * exp(-p["a"] * p["L1"]))
)
p.update(old_p)
return result
def check_consistency(p=p, extra_par={}):
old_p = p.copy()
p.update(extra_par)
wJ = sqrt(p["mJ"] / p["DJ"])
wA = sqrt(p["mA"] / p["DA"])
assert p["r"] >= 1.0
assert 1.0 / p["r"] * (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"])) <= 1.0
p.update(old_p)
def A(x, p=p, extra_par={}):
old_p = p.copy()
p.update(extra_par)
assert x >= p["L1"]
assert x <= p["L1"] + p["d"]
wJ = sqrt((p["mJ"]) / p["DJ"])
wA = sqrt((p["mA"]) / p["DA"])
result = (
cosh(wA * (p["L1"] + p["d"] - x))
/ cosh(wA * p["d"])
* (
1 / (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"]))
- 1.0 / (p["r"] * p["d"] * exp(-p["a"] * p["L1"]))
)
)
p.update(old_p)
return result
def G(l, p=p, extra_par={}):
old_p = p.copy()
p.update(extra_par)
wJ = sqrt((p["mJ"] + l) / p["DJ"])
wA = sqrt((p["mA"] + l) / p["DA"])
result = exp(l * p["t1"]) * (cosh(wJ * p["L1"]) + wA / wJ * tanh(wA * p["d"]) * sinh(wJ * p["L1"]))
p.update(old_p)
return result
def F(x, p=p):
return (1 - p['b'] * mp.exp(-x * p['t2'])) * mp.exp(x * p['t1']) \
* (mp.cosh(mp.sqrt((p['mJ'] + x)/p['DJ']) * p['L1']) \
+ mp.sinh(mp.sqrt((p['mJ'] + x)/p['DJ']) * p['L1']) \
* mp.sqrt((p['mJ'] + x)/p['DJ']) / mp.sqrt((p['mA'] + x)/p['DA']) \
* mp.coth(mp.sqrt((p['mA'] + x)/p['DA']) * p['d'])) - p['a']
