def N_ref_macs(*, theta, mu):
    """
    Calculate N0 reference size for a macs command.
    Macs has a different meaning for the '-t theta' argument.
    Theta is the 'mutation rate per site per 4N generations'.
    theta = mu / (4 N)
    """
    return theta / (4 * mu)