How to see that $$\int_0^{\infty} \sqrt{a^2+b^2} \exp \left[ -\frac{a^2+b^2}{c} \right] da= cb^2K_2 \left[ \frac{b}{c} \right],$$ where $$ K_\alpha(x) = \int_0^\infty e^{-x\cosh t} \cosh \alpha t \,dt ,$$ is the modified Bessel functions of the second kind.
EDIT
I realized that I made a stupid mistake in previous calculations. Now I get the following integral for which I should understand how to get the final result, sorry for my ignorance.
$$\int_0^{\infty} a \sqrt{a^2-b^2} \exp \left[ -\frac{a}{c} \right] da= cb^2K_2 \left[ \frac{b}{c} \right],$$
$$I=\int_0^{\infty} \sqrt{a^2+b^2} \exp \left[ -\frac{a^2+b^2}{c} \right] da=\frac{c\sqrt{\pi } }{2} e^{-\frac{b^2}{c}} \, U\left(-\frac{1}{2},0,\frac{b^2}{c}\right)$$ where appears the Tricomi hypergeometric function.
This can be also expressed as $$I=\frac{b^2}{4} e^{-\frac{b^2}{2 c}} \left(K_0\left(\frac{b^2}{2 c}\right)+K_1\left(\frac{b^2}{2 c}\right) \right)$$ but this does not seem to be the rhs (except if I am mistaken).