if measure_with_energy:
        se_suc_prob = np.sqrt(
            np.count_nonzero(np.array(response.energies) <= ref_energy)
            / (len(response.energies) - 1)
        )