[step:Apply Maker's theorem to the non-stationary ergodic average]
The functions $f_j$ and $f_\infty$ belong to $L^1(X,\mathcal B,\mu)$, and the previous step gives $f_j\to f_\infty$ in $L^1$ and almost everywhere. Define the invariant $\sigma$-algebra
\begin{align*}
\mathcal I_T:=\{A\in\mathcal B:T^{-1}A=A\text{ modulo }\mu\text{-null sets}\}.
\end{align*}
Let
\begin{align*}
\mathbb E_\mu[\cdot\mid\mathcal I_T]:L^1(X,\mathcal B,\mu)\to L^1(X,\mathcal I_T,\mu)
\end{align*}
denote conditional expectation with respect to $\mathcal I_T$. We use Maker's theorem in the tail-envelope form: if $g_j\to g_\infty$ almost everywhere and in $L^1$, and if the map $R_N^g:X\to[0,\infty]$ defined by
\begin{align*}
R_N^g(x):=\sup_{j\geq N}|g_j(x)-g_\infty(x)|
\end{align*}
belongs to $L^1(X,\mathcal B,\mu)$ with
\begin{align*}
\int_X R_N^g(x)\,d\mu(x)\to0,
\end{align*}
then
\begin{align*}
\frac{1}{n}\sum_{k=0}^{n-1}(g_k-g_\infty)(T^k x)\to0
\end{align*}
for $\mu$-a.e. $x$ and in $L^1(X,\mathcal B,\mu)$. We now verify the tail-envelope hypothesis using Breiman's maximal convergence theorem for finite conditional information. In the form needed here, the theorem says: if $\mathcal Q$ is a finite measurable partition of a probability space $(X,\mathcal B,\mu)$, if $(\mathcal A_j)_{j\geq0}$ is an increasing sequence of sub-$\sigma$-algebras, if
\begin{align*}
\mathcal A_\infty:=\bigvee_{j=0}^{\infty}\mathcal A_j,
\end{align*}
and if the maps $F_j:X\to[0,\infty]$ and $F_\infty:X\to[0,\infty]$ are defined by
\begin{align*}
F_j:=I_\mu(\mathcal Q\mid\mathcal A_j)
\end{align*}
and
\begin{align*}
F_\infty:=I_\mu(\mathcal Q\mid\mathcal A_\infty),
\end{align*}
then for each $N\in\mathbb N$ the maximal tail map $S_N:X\to[0,\infty]$ defined by
\begin{align*}
S_N(x):=\sup_{j\geq N}|F_j(x)-F_\infty(x)|
\end{align*}
belongs to $L^1(X,\mathcal B,\mu)$ and satisfies
\begin{align*}
\int_X S_N(x)\,d\mu(x)\to0.
\end{align*}
The hypotheses match our situation with $\mathcal Q=\mathcal P$, $F_j=f_j$, and $F_\infty=f_\infty$, because the previous step defined $\mathcal A_j=\mathcal G\vee\mathcal P_{[1,j]}$, $\mathcal A_\infty=\mathcal G\vee\bigvee_{i=1}^{\infty}T^i\mathcal P$, and $f_j=I_\mu(\mathcal P\mid\mathcal A_j)$. Therefore the map $R_N:X\to[0,\infty]$ defined by
\begin{align*}
R_N(x):=\sup_{j\geq N}|f_j(x)-f_\infty(x)|.
\end{align*}
belongs to $L^1(X,\mathcal B,\mu)$ and satisfies
\begin{align*}
\int_X R_N(x)\,d\mu(x)\to0.
\end{align*}
Therefore Maker's theorem applies to the sequence $(f_j)_{j\ge0}$ and gives
\begin{align*}
\frac{1}{n}\sum_{k=0}^{n-1} f_k(T^k x)\to \mathbb E_\mu[f_\infty\mid\mathcal I_T](x)
\end{align*}
for $\mu$-a.e. $x$ and in $L^1(X,\mathcal B,\mu)$. Since the system is ergodic, every member of $\mathcal I_T$ has $\mu$-measure $0$ or $1$, so
\begin{align*}
\mathbb E_\mu[f_\infty\mid\mathcal I_T]=\int_X f_\infty(y)\,d\mu(y)
\end{align*}
$\mu$-a.e. Combining this with the decomposition from the previous step gives
\begin{align*}
\frac{1}{n}I_n^{\mathcal P\mid\mathcal G}(x)\to \int_X f_\infty(y)\,d\mu(y)
\end{align*}
for $\mu$-a.e. $x$ and in $L^1(X,\mathcal B,\mu)$.
[/step]