Commutator of scalar fields

So, in the calculation of $D(t,r) = \left[ \phi(x) , \phi(y) \right]$, where $t= x^0 - y^0,~ \vec{r} = \vec{x} - \vec{y}$ you need to calculate the following integral $$D(t,r) = \frac{1}{2\pi^2 r} \int\limits_0^\infty dp \frac{ p \sin(p r) \sin \left[(p^2 + m^2)^{1/2} t \right]} { (p^2 + m^2 )^{1/2}}$$ For $m=0$, the integral is simple. We get $$D(t,r) = \frac{1}{4\pi r} \left[ \delta(t - r) - \delta(t + r) \right]$$ I even know what the answer for $m \neq 0$. I have no idea how to calculate it though. Any help?

Answer

Using Gradshteyn and Ryzhik (seventh edition) 3.876 (1) $$\int_0^\infty \frac{\sin{(p \sqrt{x^2+a^2})}}{\sqrt{x^2+a^2}} \cos(b~x)dx=\frac{\pi}{2} J_0(a\sqrt{p^2-b^2}) ~~[0