Now I have a Heaviside function $H(K-\sqrt{k^2+l^2})$ in a 2D $\hat k$ space, where $k$ and $l$ are two variables in that space. In a paper, it is said that the inverse Fourier transform of this Heaviside function is: $$ \frac{K}{\sqrt{x^2+y^2}} J_1(K\sqrt{x^2+y^2}) $$ where $x$ and $y$ are variables in 2D real space and $J_1$ is the first order Bessel function.
How to derive this result?
Attempts:
I try to perform this transform by: $$ \int_{-\infty}^{+\infty} d k \int_{-\infty}^{+\infty} d l \ H(K-\sqrt{k^2+l^2}) e^{i2\pi(kx+ly)} $$ Using the properity of Heaviside function, the integral becomes: $$ \iint_{k^2+l^2 \leq K^2} dk \ dl \ e^{i2\pi(kx+ly)} $$ Switch to polar coordinates(set $k=r\cos\theta$ and $l=r\sin\theta$): $$ \int_0^{K} dr \int_0^{2\pi} d\theta \ r e^{i2\pi r(x\cos\theta+y\sin\theta)} $$ After integral over $r$, we have: $$ \int_0^{2\pi} d\theta \ \frac{-1+[1-i2\pi K(x\cos\theta+y\sin\theta)]e^{i2\pi K(x\cos\theta+y\sin\theta)}} {4\pi^2(x\cos\theta+y\sin\theta)^2} $$ I don't know how to deal with this integral and how to transform it into Bessel function.