code: def _tilt(omega, b, psi):
      
      return np.cosh(psi/2.0)**b * np.exp(-psi**2/2.0 * omega)
docs: Compute the tilt of the PG density for value omega
      and tilt psi.
  
      :param omega: point at which to evaluate the density
      :param psi: tilt parameter