[guided]We now turn tail decay into a moment estimate. The correct tool is the layer-cake formula, applied to the non-negative measurable function $Z^q$ on the probability space $(\mathbb{R}^n,\mathcal{B}(\mathbb{R}^n),\mu)$. It gives
\begin{align*}
\int_{\mathbb{R}^n}Z(x)^q\,d\mu(x)
&=q\int_0^\infty t^{q-1}\mu(\{x\in\mathbb{R}^n:Z(x)>t\})\,d\mathcal{L}^1(t).
\end{align*}
The tail estimate is only stated for $t\ge 2$, so we split the one-dimensional [Lebesgue integral](/page/Lebesgue%20Integral) into the intervals $(0,2)$ and $[2,\infty)$. On $(0,2)$ we use only that $\mu$ is a probability measure:
\begin{align*}
q\int_0^2 t^{q-1}\mu(\{x\in\mathbb{R}^n:Z(x)>t\})\,d\mathcal{L}^1(t)
\le q\int_0^2 t^{q-1}\,d\mathcal{L}^1(t)=2^q.
\end{align*}
On $[2,\infty)$ we use Borell's tail bound:
\begin{align*}
q\int_2^\infty t^{q-1}\mu(\{x\in\mathbb{R}^n:Z(x)>t\})\,d\mathcal{L}^1(t)
\le q\int_2^\infty t^{q-1}\exp\left(-\frac{t}{2e}\right)\,d\mathcal{L}^1(t).
\end{align*}
We enlarge the integration domain from $[2,\infty)$ to $[0,\infty)$ because the integrand is non-negative. This gives
\begin{align*}
q\int_2^\infty t^{q-1}\exp\left(-\frac{t}{2e}\right)\,d\mathcal{L}^1(t)
\le q\int_0^\infty t^{q-1}\exp\left(-\frac{t}{2e}\right)\,d\mathcal{L}^1(t).
\end{align*}
Now substitute $u=t/(2e)$. Under this substitution $t=2eu$ and $d\mathcal{L}^1(t)=2e\,d\mathcal{L}^1(u)$, while the domain $[0,\infty)$ is mapped onto $[0,\infty)$. Hence
\begin{align*}
q\int_0^\infty t^{q-1}\exp\left(-\frac{t}{2e}\right)\,d\mathcal{L}^1(t)
&=q(2e)^q\int_0^\infty u^{q-1}e^{-u}\,d\mathcal{L}^1(u)\\
&=q(2e)^q\Gamma(q).
\end{align*}
Combining the small-tail and large-tail pieces yields
\begin{align*}
\int_{\mathbb{R}^n}Z(x)^q\,d\mu(x)\le 2^q+q(2e)^q\Gamma(q).
\end{align*}[/guided]