Use momentum conservation to justify
$$
\int_{0}^{1} d z z\left(q_{S}\left(z, Q^{2}\right)+g\left(z, Q^{2}\right)\right)=1
$$
Hence, determine the $\delta(1-z)$ term in $P_{g g}$ and verify
$$
P_{g g}(z)=6\left(\frac{1-z}{z}+\frac{z}{(1-z)_{+}}+z(1-z)\right)+\left(\frac{11}{2}-\frac{n_{f}}{3}\right) \delta(1-z),
$$
where $n_{f}$ is the number of quark flavors.