EDIT: This question has been edited thanks to a comment. One of my definitions was wrong, so I have rewritten the whole question.
I was reading this paper about $T \bar{T}$ deformations of $2d$-QFTs in the plane. All is fine until the beginning of section 3. There they start talking about the limit of a large number of degrees of freedom. This means large central charge $c$ in a CFT and something similar for general QFTs. They say something like:
In the large $c$ limit correlation functions of $T_{ij}$ factorize ($T_{ij}$ is the energy-momentum tensor in Euclidean coordinates).
The connected contribution to an $n$-point function of energy-momentum tensors is proportional to $c$ at large $c$, so that when we compute a general correlation function and look at the contribution to it which is a product of $k$ connected components, then this will scale as $c^k$.
I will assume connected correlation functions are defined in a similar way as connected parts of the $S$-matrix in Weinberg's QFT I. For instance, for a 6 point function of the energy-mometum tensor:
$$ \langle T_1 T_2 T_3 T_4 T_5 T_6 \rangle = \langle T_1 T_2 T_3 T_4 T_5 T_6 \rangle_c \\ \qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad + \langle T_1 T_2 T_3 \rangle_c \langle T_4 T_5 T_6 \rangle_c + \text{permutations} \\ \qquad \qquad \qquad \qquad \qquad \qquad + \langle T_1 T_2 \rangle_c \langle T_3 T_4 \rangle_c \langle T_5 T_6 \rangle_c + \text{permutations} $$
There is no term like $\langle T_1\rangle_c \langle T_2 \rangle_c \langle T_3 \rangle_c \langle T_4 \rangle_c \langle T_5 \rangle_c \langle T_6 \rangle_c$ because $ \langle T_{ij} \rangle_c = \langle T_{ij} \rangle = 0$ ( in the plane you can set this to zero). Note that this implies $\langle T_{ij} T_{kl} \rangle_c = \langle T_{ij} T_{kl} \rangle$.
QUESTION: Why do connected $n$-point functions scale as $c$ at large $c$? I mean, why is the following true for any $n$?
$$ \langle T_{i_1 j_1} T_{i_2 j_2} ... T_{i_n j_n} \rangle_c \underset{c \rightarrow \infty}{\propto} c .$$
The only thing that comes to my mind that relates the energy-momentum tensor correlation functions and $c$ is the OPE
$$T(z) T(w) \sim \frac{c/2}{(z-w)^4} + \frac{2 T(w)}{(z-w)^2} + \frac{\partial T(w)}{z-w}. $$
This works inside correlation functions, and I don't know how would it act inside connected correlation functions. Moreover, if this was the correct approach, it would give a factor of $c$ for every pair of $T$s inside the correlator. For instance
$$ \langle T(z_1) T(z_2) T(z_3) T(z_4) \rangle \sim \frac{1/4}{(z_1-z_2)^4(z_3-z_4)^4} c^2. $$