we shall use leading term for energy dependence of the mid rapidity multiplicity distribution (\[dNdy\]) which gives $$\frac{d^{2}N_{\alpha}}{dydp_{\text{T}}^{2}}=\frac{\gamma_{\alpha} }
{2T_{\alpha}^{2}}b_{\alpha}\left(  \frac{W}{Q_{0}}\right)  ^{2\lambda
/(2+\lambda)}\left[  1+\frac{Q_{0}\left(  \frac{W}{Q_{0}}\right)
^{\lambda/(2+\lambda)}\tau^{1/(2+\lambda)}}{n_{\alpha}\,T_{\alpha}}\right]
^{-n_{\alpha}}.\label{dNetc}%$$ The right hand side of Eq.(\[dNetc\]) should be an energy independent function of scaling variable $\tau$ only. Within approximations used so far there exists a simple solution to this requirement:$$T_{\alpha}=\kappa_{\alpha} \,Q_{0}\left(  \frac{W}{Q_{0}}\right)
^{\lambda/(2+\lambda)}\label{Tsol}%$$ where $\kappa_{\alpha}$ is a constant. Therefore GS predictions for Tsallis parameters are that $T_{\alpha}$ depends on energy as a power[^1], whereas $n_{\alpha}$ is a constant. A complete fit to high energy data for charged (unidentified) particles from NA49 energies up to the LHC [@Rybczynski:2012vj] shows rather small variation (of the order of 10%) of parameter $q(W)$ related to $n$ from Eq. (\[Tsallis\]) in the following way: $$n(W)=\frac{1}{q(W)-1}$$ which, however, translates into rather strong energy dependence of $n(W)$, especially for smaller energies where $q(W)$ is only slightly bigger than 1. Similar conclusion – as far as the energy dependence of the Tsallis parameters for identified particles is concerned – has been found in Ref. [@Cleymans:2013rfq] with temperature hardly depending on energy. One should note, however, that the multiplicity distribution used in Ref. [@Cleymans:2013rfq] slightly differs from the one of Eq. (\[Tsallis\]).

