
Measuring the twopoint correlation function of the galaxies in the Universe
gives access to the underlying dark matter distribution, which is related to
cosmological parameters and to the physics of the primordial Universe. The
estimation of the correlation function for current galaxy surveys makes use of
the LandySzalay estimator, which is supposed to reach minimal variance. This
is only true, however, for a vanishing correlation function. We study the
LandySzalay estimator when these conditions are not fulfilled and propose a
new estimator that provides the smallest variance for a given survey geometry.
Our estimator is a linear combination of ratios between paircounts of data
and/or random catalogues (DD, RR and DR). The optimal combination for a given
geometry is determined by using lognormal mock catalogues. The resulting
estimator is biased in a modeldependent way, but we propose a simple iterative
procedure for obtaining an unbiased model independent estimator.Our method can
be easily applied to any dataset and requires few extra mock catalogues
compared to the standard LandySzalay analysis. Using various sets of simulated
data (lognormal, secondorder LPT and NBody), we obtain a 2025% gain on the
error bars on the twopoint correlation function for the SDSS geometry and
$\Lambda$CDM correlation function. When applied to SDSS data (DR7 and DR9), we
achieve a similar gain on the correlation functions, which translates into a
1015% improvement over the estimation of the densities of matter $\Omega_m$
and dark energy $\Omega_\Lambda$ in an open $\Lambda$CDM model. The constraints
derived from DR7 data with our estimator are similar to those obtained with the
DR9 data and the LandySzalay estimator, which covers a volume twice as large
and has a density that is three times higher.

Largescale structures in the Universe are a powerful tool to test
cosmological models and constrain cosmological parameters. A particular feature
of interest comes from Baryon Acoustic Oscillations (BAOs), which are sound
waves traveling in the hot plasma of the early Universe that stopped at the
recombination time. This feature can be observed as a localized bump in the
correlation function at the scale of the sound horizon $r_s$. As such, it
provides a standard ruler and a lot of constraining power in the correlation
function analysis of galaxy surveys. Moreover the detection of BAOs at the
expected scale gives a strong support to cosmological models. Both of these
studies (BAO detection and parameter constraints) rely on a statistical
modeling of the measured correlation function $\hat{\xi}$. Usually $\hat{\xi}$
is assumed to be gaussian, with a mean $\xi_\theta$ depending on the
cosmological model and a covariance matrix $C$ generally approximated as a
constant (i.e. independent of the model). In this article we study whether a
realistic modeldependent $C_\theta$ changes the results of cosmological
parameter constraints compared to the approximation of a constant covariance
matrix $C$. For this purpose, we use a new procedure to generate lognormal
realizations of the Luminous Red Galaxies sample of the Sloan Digital Sky
Survey Data Release 7 to obtain a modeldependent $C_\theta$ in a reasonable
time. The approximation of $C_\theta$ as a constant creates small changes in
the cosmological parameter constraints on our sample. We quantify this modeling
error using a lot of simulations and find that it only has a marginal influence
on cosmological parameter constraints for current and nextgeneration galaxy
surveys. It can be approximately taken into account by extending the $1\sigma$
intervals by a factor $\approx 1.3$.

We obtain constraints on cosmological parameters from the spherically
averaged redshiftspace correlation function of the CMASS Data Release 9 (DR9)
sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this
information with additional data from recent CMB, SN and BAO measurements. Our
results show no significant evidence of deviations from the standard
flatLambda CDM model, whose basic parameters can be specified by Omega_m =
0.285 + 0.009, 100 Omega_b = 4.59 + 0.09, n_s = 0.96 + 0.009, H_0 = 69.4 +
0.8 km/s/Mpc and sigma_8 = 0.80 + 0.02. The CMB+CMASS combination sets tight
constraints on the curvature of the Universe, with Omega_k = 0.0043 + 0.0049,
and the tensortoscalar amplitude ratio, for which we find r < 0.16 at the 95
per cent confidence level (CL). These data show a clear signature of a
deviation from scaleinvariance also in the presence of tensor modes, with n_s
<1 at the 99.7 per cent CL. We derive constraints on the fraction of massive
neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51
eV. We find no signature of a deviation from a cosmological constant from the
combination of all datasets, with a constraint of w_DE = 1.033 + 0.073 when
this parameter is assumed timeindependent, and no evidence of a departure from
this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1  a). The
achieved accuracy on our cosmological constraints is a clear demonstration of
the constraining power of current cosmological observations.

Baryon Acoustic Oscillations (BAO) are a feature imprinted in the density
field by acoustic waves travelling in the plasma of the early universe. Their
fixed scale can be used as a standard ruler to study the geometry of the
universe. BAO have been previously detected using correlation functions and
power spectra of the galaxy distribution. In this work, we present a new method
for the detection of the realspace structures associated with this feature.
These baryon acoustic structures are spherical shells with a relatively small
density contrast, surrounding high density central regions. We design a
specific wavelet adapted to the search for shells, and exploit the physics of
the process by making use of two different mass tracers, introducing a specific
statistic to detect the BAO features. We show the effect of the BAO signal in
this new statistic when applied to the Lambda  Cold Dark Matter (LCDM) model,
using an analytical approximation to the transfer function. We confirm the
reliability and stability of our method by using cosmological Nbody
simulations from the MareNostrum Institut de Ci\`encies de l'Espai (MICE). We
apply our method to the detection of BAO in a galaxy sample drawn from the
Sloan Digital Sky Survey (SDSS). We use the `Main' catalogue to trace the
shells, and the Luminous Red Galaxies (LRG) as tracers of the high density
central regions. Using this new method, we detect, with a high significance,
that the LRGs in our sample are preferentially located close to the centres of
shelllike structures in the density field, with characteristics similar to
those expected from BAOs. We show that stacking selected shells, we can find
their characteristic density profile. We have delineated a new feature of the
cosmic web, the BAO shells. As these are real spatial structures, the BAO
phenomenon can be studied in detail by examining those shells.

Baryon Acoustic Oscillations are a feature imprinted in the galaxy
distribution by acoustic waves traveling in the plasma of the early universe.
Their detection at the expected scale in largescale structures strongly
supports current cosmological models with a nearly linear evolution from
redshift approximately 1000, and the existence of dark energy. Besides, BAOs
provide a standard ruler for studying cosmic expansion. In this paper we focus
on methods for BAO detection using the correlation function measurement. For
each method, we want to understand the tested hypothesis (the hypothesis H0 to
be rejected) and the underlying assumptions. We first present wavelet methods
which are mildly modeldependent and mostly sensitive to the BAO feature. Then
we turn to fully modeldependent methods. We present the most often used method
based on the chi^2 statistic, but we find it has limitations. In general the
assumptions of the chi^2 method are not verified, and it only gives a rough
estimate of the significance. The estimate can become very wrong when
considering more realistic hypotheses, where the covariance matrix of the
measurement depends on cosmological parameters. Instead we propose to use a new
method based on two modifications: we modify the procedure for computing the
significance and make it rigorous, and we modify the statistic to obtain better
results in the case of varying covariance matrix. We verify with simulations
that correct significances are different from the ones obtained using the
classical chi^2 procedure. We also test a simple example of varying covariance
matrix. In this case we find that our modified statistic outperforms the
classical chi^2 statistic when both significances are correctly computed.
Finally we find that taking into account variations of the covariance matrix
can change both BAO detection levels and cosmological parameter constraints.