There is an excellent explanation of this in Chapter 4 of L. N. Trefethen's Approximation Theory and Approximation Practice (henceforth ATAP; the first 6 chapters are available for free online). I will try to summarize it in your notation, but I recommend reading his explanation.
Chebyshev interpolation
Interpolation at the Chebyshev points gives your $\Pi^{CG}_n f$, which can also be written in terms of the first $n$ Chebyshev polynomials:
$$\Pi^{CG}_n f = \sum_{k=0}^n c_k T_k(x).$$
The question then becomes:
How are the coefficients $c_k$ of the degree-$n$ Chebyshev interpolant related to the coefficients $\hat{f}_k$ of the infinite Chebyshev series?
Aliasing
It turns out that on the Chebyshev grid with $n+1$ points, the following Chebyshev polynomials all take exactly the same values:
$$T_{2jn\pm m}(x), \ \ \ \ j = 0, 1, 2 \dots$$
The above statement is valid as long as $0\le m \le n$ and where the subscript is positive (of course). See Thm. 4.1 of ATAP. Thus each $T_k$ for $k > n$ is "aliased" to appear like some corresponding $T_k$ with $k\le n$.
Relation between the Chebyshev interpolant and the Chebyshev series
Once you know about aliasing, the connection between the coefficients $c_k$ and $\hat{f}_k$ is immediate. Since the Chebyshev interpolant and the infinite Chebyshev series must be equal at the Chebyshev points, we have
$$
\sum_{k=0}^n c_k T_k(x_j) = \sum_{k=0}^\infty \hat{f}_k T_k(x_j),
$$
where $x_j$ is any of the Chebyshev points.
But from the aliasing result, we know that each term with $k>n$ in the right-hand-side sum can be rewritten in terms of $T_k(x_j)$ for some $k\le n$. The result is Thm. 4.2 of ATAP, which gives the formulae (for the degree-$n$ interpolant coefficients)
\begin{align}
c_0 - f_0 & = \sum_{j=1}^\infty f_{2jn} \\
c_n - f_n & = \sum_{j=1}^\infty f_{(2j+1)n} \\
c_k - f_k & = \sum_{j=1}^\infty (f_{2jn+k} + f_{2jn-k}) & 1\le k \le n-1.
\end{align}
So the difference between the coefficients of the Chebyshev interpolant and the Chebyshev projection is given by certain sums of coefficients of the remainder of the Chebyshev series.
One of the major themes of ATAP is that the Chebyshev interpolant is almost as good as the Chebyshev projection for approximation purposes; see Theorems 7.2, 8.2, and 16.1 therein. The great virtue of the interpolant is that it is much easier to compute.
Best Answer
There is a missing absolute sign and it should be max not integral, i.e. $$ \max_{x\in[-1,1]}\lvert (x-x_0)\dots(x-x_n)\rvert=\max_{x\in[-1,1]}\lvert 2^{-n}T_{n+1}(x)\rvert=2^{-n}. $$
Edit: correcting not integral. The formula for error of polynomial interpolation with nodes $x_i$, $i=0,1,\dots,n$ is $$ f(x)-P_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}\prod_{i=0}^n (x-x_i) $$ so the maximum error is never more than $\sup\frac1{(n+1)!}\lvert f^{(n+1)}\rvert$ when using Chebyshev nodes.