from aart import *

# Johnson’s SU distribution (e.g., Eq.12 2008.03879)
def profile(r,a,gammap,mup,sigmap):
    #(Eq. 58 P1)
    return np.exp(-(1/2)*(gammap+ np.arcsinh((r-mup)/sigmap))**2)/sqrt((r-mup)**2+sigmap**2)