Spectral functions at finite temperature and two-loop order are investigated, for a medium consisting of massless particles. We consider them in the timelike and spacelike domains, allowing the propagating particles to be any valid combination of bosons and fermions. Divergences (if present) are analytically derived and set aside for the remaining finite part to be calculated numerically.To illustrate the utility of these 'master' functions, we consider transverse and longitudinal parts of the QCD vector channel spectral function.