I am computing a two-point correlator in 4D Euclidean space and I am struggling with one particular term. I have found that in momentum space my correlator goes as $$\langle \mathcal{O}(p)\mathcal{O}(q) \rangle \sim \delta(p+q) p^4\ln p.$$ I now want to Fourier transform this to get the correlator in position space. In other words I need to compute the integral $$\langle \mathcal{O}(x)\mathcal{O}(y) \rangle\sim\int d^4 p \; p^4 \ln(p) e^{-ip\cdot(x-y)}.$$
Some hints I have gotten, but haven't brought me far are the following:
transform into polar coordinates, $$\int d p \, d\theta \, d\psi \, d\phi \; p^7 \sin \theta \sin^2 \psi \ln(p) e^{-i|p|(|x|-|y|)\cos \theta}.$$ I'm not sure if this is correct because we want to identify the $\theta$ angle as the one between the vectors $p$ and $x$.
use a regulator function of the form $e^{-\epsilon p}$ for the 'radial' integral. I am confused here because I have never worked explicitly with these regulators before. I don't know how to compute an integral of the type $$\int dp \; p^7 \ln p e^{-ip x}.$$
consider the integral representations of the Bessel functions. I think this means use $$J_0(x) = \frac{1}{2\pi} \int^\pi_{-\pi} e^{ix\sin \phi}d\phi.$$ I don't know how to proceed once we have the Bessel function.