code: def fap_davies(Z, fmax, t, y, dy, normalization='standard'):
      
      N = len(t)
      fap_s = fap_single(Z, N, normalization=normalization)
      tau = tau_davies(Z, fmax, t, y, dy, normalization=normalization)
      return fap_s + tau
docs: Davies upper-bound to the false alarm probability
  
      (Eqn 5 of Baluev 2008)