pyrats/physics.py
def N_merger_until_z(ximin, M0, z, z0):
    """
    Return the number of merger with a mass ratio greater than xmin
    for a halo of mass M0 at z0
    between z and z0 < z
    """
    intxi, residual = integrate.quad(
        lambda x: x**-1.995 * np.exp((x / 0.00972)**0.263), ximin, 1)
    result, residual = integrate.quad(
        lambda zz: (
            0.0104 * (1 + zz)**0.0993 * intxi * MofZ(M0, zz, z0)**0.133),
        z0, z)

    return result

    Return the number of merger with a mass ratio greater than xmin
    for a halo of mass M0 at z0
    between z and z0 < z
    
Return the number of merger with a mass ratio greater than xmin
for a halo of mass M0 at z0
between z and z0 < z
Return the number of merger with a mass ratio greater than xmin for a halo of mass M0 at z0 between z and z0 < z
def N_merger_until_z(ximin, M0, z, z0):
    intxi, residual = integrate.quad(
        lambda x: x**-1.995 * np.exp((x / 0.00972)**0.263), ximin, 1)
    result, residual = integrate.quad(
        lambda zz: (
            0.0104 * (1 + zz)**0.0993 * intxi * MofZ(M0, zz, z0)**0.133),
        z0, z)
    return result
def N_merger_until_z ( ximin , M0 , z , z0 ) : intxi , residual = integrate . quad ( lambda x : x ** - 1.995 * np . exp ( ( x / 0.00972 ) ** 0.263 ) , ximin , 1 ) result , residual = integrate . quad ( lambda zz : ( 0.0104 * ( 1 + zz ) ** 0.0993 * intxi * MofZ ( M0 , zz , z0 ) ** 0.133 ) , z0 , z ) return result
Return the number of merger with a mass ratio greater than xmin
for a halo of mass M0 at z0
between z and z0 < z
Return the number of merger with a mass ratio greater than xmin for a halo of mass M0 at z0 between z and z0 < z
"""
    Return the number of merger with a mass ratio greater than xmin
    for a halo of mass M0 at z0
    between z and z0 < z
    """
ximin
M0
z
z0