Let $(\Omega,\mathcal F,P)$ be a probability space and let $([0,1],\mathcal B[0,1],\lambda$) be the unit interval with Lebesgue measure on the Borel subsets of $[0,1]$. Let $f:[0,1]\times \Omega \to \mathbb R$ be $\mathcal B[0,1]\otimes \mathcal F$ measurable and bounded. I want to show the existence of a subsequence $(n_k)$ such that
$$\frac{1}{n_k}\sum_{i=1}^{n_k} f(t_i,\omega) \overset{L^1(\Omega,\mathcal F,P)}\to \int f(t,\omega)d\lambda \quad \text{as } k\to \infty $$
for $\lambda^\infty$-almost every sequence $(t_i)\subset [0,1]$, where $\lambda^\infty$ denotes the product measure on the $\sigma$-algebra $\mathcal B^\infty:=\bigotimes_{n=1}^\infty \mathcal B[0,1]$.
I tried the following:
Fix $\omega\in\Omega$. For each $j\geq1$, let $((t_i),\omega)\mapsto f(t_j,\omega)$ denote the composition of the maps $((t_i),\omega)\mapsto (t_j,\omega)$ and $(t,\omega) \mapsto f(t,\omega)$. Since the $j$th projection map $(t_i)\mapsto t_j$ on $[0,1]^{\infty}:=\prod_{i=1}^{\infty} [0,1]$ is $\mathcal B^\infty$ measurable, we see that $f(t_j,\omega)$ is $\mathcal B^\infty\otimes \mathcal F$ measurable. Moreover, by construction of the product measure $\lambda^\infty$, the projection maps are i.i.d. random variables on $([0,1]^{\infty},\mathcal B^{\infty},\lambda^{\infty})$, and so by the SLLN we get
$$\frac{1}{n}\sum_{i=1}^n f(t_i,\omega) \overset{\lambda^\infty\text{-a.s}}\to \int f(t,\omega)d\lambda \quad \text{as } n\to \infty $$
for each $\omega\in \Omega$. For each $n\geq 1$, define the map $$((t_i),\omega)\mapsto S^n((t_i),\omega):=\frac{1}{n}\sum_{i=1}^n f(t_i,\omega) -\int f(t,\omega)d\lambda.$$ By Fubini's theorem $\omega\to \int f(t,\omega)d\lambda$ is $\mathcal F$-measurable, and so $S^n((t_i),\omega)$ is seen to be $\mathcal B^\infty\otimes \mathcal F$ measurable. Also $S^n((t_i),\omega)$ is bounded since $f$ is bounded. Therefore the DCT gives
$$\int |S^n((t_i),\omega)| d\lambda^{\infty} \to 0 \quad \text{as } n\to \infty$$
for each $\omega\in \Omega$. By Fubini's theorem, these integrals are $\mathcal F$-measurable functions of $\omega$, and they are bounded, and so another application of the DCT gives
$$\int \int |S^n((t_i),\omega)| d\lambda^{\infty}dP \to 0 \quad \text{as } n\to \infty$$
We can interchange the order of integration by Fubini's theorem to obtain
$$\int \int |S^n((t_i),\omega)| dPd\lambda^{\infty} \to 0 \quad \text{as } n\to \infty$$
This says that $\int |S^n((t_i),\omega)| dP \to 0 $ as $n\to \infty $ w.r.t. the $L^1([0,1]^{\infty} ,\mathcal B^{\infty},\lambda^{\infty})$ norm. We can extract a subsequence $(n_k)$ converging $\lambda^{\infty}$-almost surely, i.e. such that
$$\frac{1}{n_k}\sum_{i=1}^{n_k} f(t_i,\omega) \overset{L^1(\Omega,\mathcal F,P)}\to \int f(t,\omega)d\lambda \quad \text{as } k\to \infty $$ for $\lambda^\infty$-almost every sequence $(t_i)\subset [0,1]$.
Is this reasoning correct? Thanks a lot for your help.