[guided]The concentration step must be done with the original class $\mathcal H$, not with $\mathcal H^\pm$, because the sharp Lipschitz constant uses the range $[0,1]$. Define $\Phi:\mathcal Z^n\to[0,\infty)$ by
\begin{align*}
\Phi(z_1,\dots,z_n):=\sup_{h\in\mathcal H}\left|\frac{1}{n}\sum_{i=1}^{n}h(z_i)-P h\right|.
\end{align*}
Then the random variable we want to control is exactly
\begin{align*}
S=\Phi(Z_1,\dots,Z_n).
\end{align*}
We verify the bounded-differences hypothesis. Fix two deterministic samples $(z_1,\dots,z_n)$ and $(z_1',\dots,z_n')$ that differ only at coordinate $j$. For a fixed $h\in\mathcal H$, the population term $P h$ is unchanged, and all empirical summands except the $j$-th one cancel. Hence
\begin{align*}
\left|\left(\frac{1}{n}\sum_{i=1}^{n}h(z_i)-P h\right)-\left(\frac{1}{n}\sum_{i=1}^{n}h(z_i')-P h\right)\right|=\frac{1}{n}|h(z_j)-h(z_j')|.
\end{align*}
Because every $h\in\mathcal H$ takes values in $[0,1]$, we have $|h(z_j)-h(z_j')|\le 1$, so the displayed difference is at most $1/n$. Passing from a fixed $h$ to the supremum cannot increase the coordinate sensitivity beyond this common bound, and therefore
\begin{align*}
|\Phi(z_1,\dots,z_n)-\Phi(z_1',\dots,z_n')|\le \frac{1}{n}.
\end{align*}
We now apply McDiarmid's bounded differences inequality, whose conclusion says that a [measurable function](/page/Measurable%20Function) of independent variables with coordinate sensitivities $c_1,\dots,c_n$ satisfies
\begin{align*}
\mathbb P\left(\Phi(Z_1,\dots,Z_n)-\mathbb E[\Phi(Z_1,\dots,Z_n)]\ge u\right)\le \exp\left(-\frac{2u^2}{\sum_{i=1}^{n}c_i^2}\right).
\end{align*}
Here $c_i=1/n$ for every $i$, so
\begin{align*}
\sum_{i=1}^{n}c_i^2=\sum_{i=1}^{n}\frac{1}{n^2}=\frac{1}{n}.
\end{align*}
Thus
\begin{align*}
\mathbb P\left(S-\mathbb E[S]\ge u\right)\le \exp(-2nu^2).
\end{align*}
Choosing
\begin{align*}
u=\sqrt{\frac{t}{2n}}
\end{align*}
makes the right-hand side equal to $e^{-t}$. Combining this concentration estimate with the symmetrization estimate
\begin{align*}
\mathbb E[S]\le 2\mathbb E\left[\mathfrak R_n(\mathcal H^\pm;Z_1,\dots,Z_n)\right]
\end{align*}
gives
\begin{align*}
\mathbb P\left(S\le 2\mathbb E\left[\mathfrak R_n(\mathcal H^\pm;Z_1,\dots,Z_n)\right]+\sqrt{\frac{t}{2n}}\right)\ge 1-e^{-t}.
\end{align*}
This is exactly the asserted high-probability uniform deviation bound.[/guided]