There is a good reason to use $z$ instead of $e^{sT}$. Before starting with any analysis, let me remind you that in analysis of signals and systems we are interested in analyzing the frequency spectrum of the signal, i.e. the Laplace transform on the imaginary line $s = jw$. And since your signal $x[n]$ is discrete, then its frequency spectrum its periodic, so its more general define $s=jwT$.
Now, let $X(z) = \mathcal{Z}\{x[n]\}$ of a causal or non-causal discrete signal $x[n]$, i.e.
$$ X(z) = \sum_{n=-\infty}^{\infty} x[n] z^{-n}. $$
Since $z\in\mathbb{C}$ we have $z = |z| e^{j\arg z}$. Without loss of generality we rewrite $|z| = r$ and $\arg z = wT$, i.e. $z=r e^{jwT}$ (note that not necessarily $r=1$). Then
$$ \begin{aligned}
X(z) &= \sum_{n=-\infty}^{\infty} x[n] z^{-n}\\
&= \sum_{n=-\infty}^{\infty} x[n] (r e^{jwT})^{-n}\\
% &= \sum_{n=-\infty}^{\infty} (x[n] r^{-n}) e^{-njwT}\\
&= \sum_{n=-\infty}^{\infty} (x[n] r^{-n}) (e^{jwT})^{-n},
\end{aligned} $$
which implies $X(z) = \left. \mathcal{L}\{x[n]r^{-n}\} \right\rvert_{s=jwT} = \mathcal{F}\{x[n]r^{-n}\}$. As a consequence, $X(z)$ is a Fourier transform more generic than the Fourier transform $X(e^{jwT}) = \mathcal{F}\{x[n]\}$ of our signal of interest.
So, if the convergence radius of $X(z)$ is less than unity then $X(e^{jwT})$ does not exist and therefore its Fourier transform does not either, which represents a problem because there are many signals with this problem of convergence, e.g. non-causal signals such as a digital image filter. Therefore, it is convenient (and even necessary in non-causal signals) to use the Z-transform.
Or informally, use $z$ instead of $\left. e^{sT} \right\rvert_{s=jwT} = e^{jwT}$ whenever you can.
We also recommend to see this link about radius convergence.
At this point, it is clear that the Z-transform has the same objective as the Laplace transform: ensure the convergence of the transform in some region of $\mathbb{C}$, where the Z-transform does it for discrete signals and Laplace transform for continuous signals.
If we use 3D rotational symmetry and define
$$x^{\mu}~=~(T,0,0,\underbrace{R}_{\geq 0})\quad\text{and}\quad x^{\prime \mu}~=~(t,r\sin\theta\cos\phi,r \sin\theta \sin\phi,r\cos\theta),$$
and
$$\rho~=~\sqrt{r^2+\alpha^2}\quad\text{and}\quad \Delta~=~R^2-T^2, $$
then OP's integral (1) becomes
$$ \begin{align}
I(R,T,\alpha)~=~& \int_{\mathbb{R}} \!\mathrm{d}t \int_{\mathbb{R}_+} \!r^2 \mathrm{d}r \int_0^{\pi} \!\sin\theta \mathrm{d}\theta\int_0^{2\pi} \!\mathrm{d}\phi~\cr
&\delta(\rho^2+R^2-2Rr\cos\theta-(T-t)^2)~\Theta(T-t) ~\delta(\rho^2-t^2)~\Theta(t)\cr
~\stackrel{R\neq 0}{=}& 2\pi \Theta(T) \int_0^T \!\mathrm{d}t \int_{\mathbb{R}_+} \!r^2 \mathrm{d}r \int_{-1}^1 \! \mathrm{d}\cos\theta~\cr
&\frac{1}{2Rr}\delta\left(\frac{\rho^2+R^2-(T-t)^2}{2Rr}-\cos\theta\right) ~\frac{\delta(\rho-t)}{\rho+t}\cr
~=~& 2\pi \Theta(T) \int_{\mathbb{R}_+} \!\frac{r^2 \mathrm{d}r}{2Rr~2\rho} ~\Theta\left(\left|\frac{\rho^2+R^2-(T-\rho)^2}{2Rr}\right|<1\right) ~\Theta(0<\rho<T)\cr
~=~& 2\pi \Theta(T) \int_{\alpha}^{\infty} \!\frac{ \mathrm{d}\rho}{4R} ~\Theta(|\Delta +2T\rho|< 2Rr) ~\Theta(\rho<T)\cr
~=~& \frac{\pi}{2R} \Theta(T-\alpha) \int_{\alpha}^T \! \mathrm{d}\rho ~\Theta(4R^2(\rho^2-\alpha^2)-(\Delta+2T\rho)^2)\cr
~=~& \frac{\pi}{2R} \Theta(T-\alpha) \int_{\alpha}^T \! \mathrm{d}\rho ~\Theta(A\rho^2+B\rho+C), \end{align}$$
where
$$A~=~4\Delta, \quad B~=~-4T\Delta \quad\text{and}\quad C~=~-\Delta^2 -4R^2\alpha^2.$$
The discriminant is
$$ D~=~B^2-4AC~=~16R^2\Delta(\Delta +\alpha^2).$$
The roots are
$$ \rho_{\pm}~=~\frac{-B\pm\sqrt{D}}{2A}~=~\frac{T\pm R\sqrt{1+\alpha^2/\Delta}}{2}. $$
Case $-\alpha^2\leq \Delta\leq 0$: Then $D\leq 0$ and the integral $I(R,T,\alpha)=0$ vanishes.
Case $\Delta> 0$: Then $D>0$ and the roots $\rho_{\pm}$ are outside the integration region $[\alpha,T]$, so the integral $I(R,T,\alpha)=0$ vanishes.
Case $\Delta \leq -\alpha^2$: Then $R\leq T$, and therefore $R\sqrt{1+\alpha^2/\Delta}\leq T$, so the root $\rho_+< T$.
$$I(R,T,\alpha)~=~\frac{\pi}{2R} \Theta(T-\alpha)~\Theta(-\Delta-\alpha^2)~\underbrace{m([\rho_-,\rho_+]\cap [\alpha,T])}_{= \text{ length of the interval}}
,$$
where $m$ denotes the Lebesgue measure.
Let us for the remainder of this answer assume that $\alpha=0$. Then the discriminant $D\geq 0$. The roots are
$$ \rho_{\pm}~=~\frac{-B\pm\sqrt{D}}{2A}~=~\frac{T\pm R}{2}. $$
It follow that the integral is supported in the future light-cone:
$$I(R,T,\alpha=0)~=~\frac{\pi}{2R} \Theta(T)~R\Theta(-\Delta)~=~\frac{\pi}{2}\Theta(T)\Theta(-\Delta).$$
(The case $R=0$ follows e.g. from continuity of the limit $R\to 0^+$.)
Best Answer
You have that $$ \int_0^{+\infty}{\delta(\cos x)e^{-x}dx}=\sum_{n=0}^{+\infty}{e^{-\left(\frac{\pi}{2}+n\pi\right)}}=\frac{e^{-\frac{\pi}{2}}}{1-e^{-\pi}}=\frac{1}{2\mathrm{sh}\left(\frac{\pi}{2}\right)} $$