Context. One important step in the statistical analysis of the Ly-alpha forest data is the study of their second order properties. Usually, this is accomplished by means of the two-point correlation function or, alternatively, the K-function. In the computation of these functions it is necessary to take into account the presence of strong metal line complexes and strong Ly-alpha lines that can hidden part of the Ly-alpha forest and represent a non negligible source of bias. Aims. In this work, we show quantitatively what are the effects of the gaps introduced in the spectrum by the strong lines if they are not properly accounted for in the computation of the correlation properties. We propose a geometric method which is able to solve this problem and is computationally more efficient than the Monte Carlo (MC) technique that is typically adopted in Cosmology studies. The method is implemented in two different algorithms. The first one permits to obtain exact results, whereas the second one provides approximated results but is computationally very efficient. The proposed approach can be easily extended to deal with the case of two or more lists of lines that have to be analyzed at the same time. Methods. Numerical experiments are presented that illustrate the consequences to neglect the effects due to the strong lines and the excellent performances of the proposed approach. Results. The proposed method is able to remarkably improve the estimates of both the two-point correlation function and the K-function.

We present a systematic comparison of some usual estimators of the 2--point correlation function, some of them currently used in Cosmology, others extensively employed in the field of the statistical analysis of point processes. At small scales, it is known that the correlation function follows reasonably well a power--law expression $\xi(r) \propto r^{-\gamma}$. The accurate determination of the exponent $\gamma$ (the order of the pole) depends on the estimator used for $\xi(r)$; on the other hand, its behavior at large scale gives information on a possible trend to homogeneity. We study the concept, the possible bias, the dependence on random samples and the errors of each estimator. Errors are computed by means of artificial catalogues of Cox processes for which the analytical expression of the correlation function is known. We also introduce a new method for extracting simulated galaxy samples from cosmological simulations.