Example #1
0
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
Example #2
0
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
Example #3
0
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)
Example #4
0
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
Example #5
0
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 N(x, l, p=p):
    '''Incomplete!!'''
    return 2. * mp.sinh(sqrt((p['mJ'] + l)/p['DJ']) * p['L1']) \
            * (mp.exp(mp.sqrt((p['mA'] + l)/p['DA']) * x) + \
            mp.exp(mp.sqrt((p['mA'] + l)/p['DA']) * (2*p['L2'] - x)))
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']