compute the spectral function with the setting in colangelo2009