We compute the integrated doublon production rate in response to a lattice modulation for two-component fermions in an optical lattice. We derive a general formula for the integrated intensity, valid in the presence of inhomogeneous potentials such as the trap, which gives the integrated intensity in terms of equal time correlation functions only. Such a formula is thus well suited for direct numerical calculations. We show that, in the limit of large repulsion for commensurate fillings, or for temperature ranges for which the hopping is incoherent, the integrated doublon spectrum is directly related to the nearest-neighbor spin-spin correlation function. We compute its temperature dependence in this regime using finite temperature quantum Monte Carlo calculation