code: def fap_baluev(Z, fmax, t, y, dy, normalization='standard'):
      
      cdf = cdf_single(Z, len(t), normalization)
      tau = tau_davies(Z, fmax, t, y, dy, normalization=normalization)
      return 1 - cdf * np.exp(-tau)
docs: Alias-free approximation to false alarm probability
  
      (Eqn 6 of Baluev 2008)