In the Figure the pseudo-Voigt fitting of the Celestine spectrum is given. The sum of the squares of the deviations is larger than that of the fit with two q-Gaussians given in the Figure . Let us note that the far-wing behavior is different. The choice of q=1.01 for the q-Gaussian as representative of a Gaussian function could be questionable. Therefore, we add the Fig. , which is showing the fit with q=1.001 and q=1.02. The minimized sum of the squares of deviations is the same for q=1.01 and q=1.001.