UPDATE: The result below can be understood as an almost sure stochastic equicontinuity condition. I don't know of any result establishing primitives of almost sure stochastic equicontinuity. If you know, this would very likely help me answer the question.

I need the following almost sure convergence result:

Let $||\cdot||$ denote a norm on a functional space (could be $||\cdot||_{\infty}$ or $||\cdot ||_2$ in $L^2$, for example). $\mathcal{G}=\{g:\mathcal{Z}\rightarrow \mathbb{R}\}$ of measurable functions. Let $Z_1,Z_2,\dots$ be $i.i.d$ random variables with $E[g(Z_i)]=0$ for any $g\in \mathcal{G}$. Let $\{\delta_n\}$ be a sequence of positive numbers such that $\delta_n\rightarrow 0$.
Then, under **ADDITIONAL ASSUMPTIONS**,
\begin{equation*}
\sup_{g\in\mathcal{G} \text{ s.t. }||g||\leq \delta_n}\left| \frac{1}{\sqrt{n}}\sum_{i=1}^n g(Z_i)\right|\rightarrow_{a.s.}0.
\end{equation*}

Note that all random variables are defined in the same probability space and are independent draws of the same distribution. The number $n$ here denotes the size of the sample, and the almost sure convergence is over the distribution of the $Z_i$ as $n\rightarrow \infty$. The space $\mathcal{G}$ does not depend on $n$, and $g$ is not random.

I need to find which reasonable assumptions could be made to establish this result. $\mathcal{G}$ is $P$-Donsker? Lipschitz functions? Bounds on moments? etc.

**Attending requests for clarification of the background:** this result is needed to establish a stochastic equicontinuity condition on a linear regression estimator with a generated regressor (it means that one of the regressors is estimated). This is a high level result, in the sense that it is a result for a generic estimator for the regressor, we want to give conditions that this estimator must satisfy, but not specify the estimator. I am using results in a paper by Chen, Linton and Keilegom (2003) which establish the consistency of the bootstrap for estimators which are based on the optimization of a function of the data, the parameter of interest and a nuisance infinite dimension parameter. So, suppose the model is $E[Y|X,W]=\beta X+\gamma h_0(W)$, but $h_0$ is not known. The term is thus estimated and the regression is done on $\hat{h}(W)$ instead. If it helps, you can think of $h_0(W)=E[V|V\leq 0, W]$, for another variable $V$ which is observed in the data, for example.(FIY: this is a huge simplification of the true situation, so don't give me suggestions about how to estimate the model I just described. In the actual model $V$ is not observed and must be predicted out of sample with machine learning and optimization methods.)

In order to apply one of the results in the paper I mentioned, I need to prove that some quantities such as the one above are $o_{a.s.}(1)$. For example, I give you a simplified version of one of them: \begin{equation*} \sup_{h,h_0\in\mathcal{H} \text{ s.t. }||h-h_0||\leq \delta_n}\left| \frac{1}{\sqrt{n}}\sum_{i=1}^n \left[X_i(h(W_i)-h_0(W_i))-E[X_i(h(W_i)-h_0(W_i))\right]\right| \end{equation*} I need to show that several objects that look more or less like the one above are $o_{a.s.}(1)$ as $n\rightarrow \infty$ for any sequence of positive $\delta_n\rightarrow 0$. Note that I can make assumptions about $\mathcal{H}$, $X$ and $W$. For example, I can say that $\mathcal{H}$ is the space of measurable negative Lipschitz functions, and I can say that $E[|X|^4]$ and $E[|W|^4]$ are finite. Basically I need conditions that would allow me to establish this result.

At the essence this is an empirical process. Define $\mathbb{G}_nf=\frac{1}{\sqrt{n}}\sum_{i=1}^n f(Z_i)-E[f(Z_i)]$ and let $\mathbb{G}$ be the Brownian Bridge, then if $\mathcal{G}$ is $P$-Donsker, $\mathbb{G}_ng\rightsquigarrow\mathbb{G}g$ for every fixed $g$. This is a convergence in distribution in a functional space. I need to show that $\sup_{||g||_{\mathcal{G}}\leq \delta_n} |\mathbb{G}_n g|\rightarrow_{a.s.} 0$.

I hope this is clearer and someone can help me. I ran out of ideas of things to try and need a fresh take.