The Galileon model is a tensor-scalar theory of gravity which explains the late acceleration of the Universe expansion with no instabilities and recovers General Relativity in the strong field limit. Most constraints obtained so far on Galileon model parameters from cosmological data were derived for the limited subset of tracker solutions and reported tensions between the model and data. This paper presents the first exploration of the general solution of the Galileon model, which is confronted against recent cosmological data for both background observables and linear perturbations, using Monte-Carlo Markov chains. As representative scenarios of the Galileon models, we study the full Galileon model with disformal coupling to matter and the uncoupled cubic Galileon model. We find that the general solution of the full Galileon model provides a good fit to CMB spectra, while the cubic Galileon model does not. When extending the comparison to BAO and SNIa data, even the general solution of the full Galileon model fails at providing a good fit to all datasets simultaneously. Tensions remain if the models are extended with an additional free parameter, such as the sum of active neutrino masses or the normalization of the CMB lensing spectrum. Finally, the multi-messenger observation of GW170817 is also discussed in the framework of the scenarios considered. The time delay between the gravitational signal and its electromagnetic counterpart was computed \textit{a posteriori} in every scenario of the full Galileon model cosmological fit chains and found to be ruled out by this observation.

Clustering analyses of spectroscopic surveys are based upon density fluctuations, which are estimated by comparing the observed tracer density field to a selection function accounting for the survey density and geometry. However, this survey selection function is commonly partly inferred from the observed data itself, leading to so-called integral constraints, for which we propose a complete derivation. We discuss the normalisation of the introduced window functions, the shot noise contribution to the integral constraint corrections and wide-angle effects. Using this formalism, we review the well-known global integral constraint, arising when the expected mean galaxy density is taken to be the measured one. Another, stronger, constraint is imposed when the radial selection function itself is estimated from the data redshift distribution, as is often the case in the literature. We find that the impact of such a radial integral constraint can be as significant as the window function effect at large scales, depending on the survey geometry. Equations for this radial integral constraint are derived within our general formalism. We assess the validity of our approach by performing a Redshift Space Distortions (RSD) analysis on mock catalogues and emphasise that our results may be even more useful for analyses focusing on larger scales. Finally, as a further application, we show that unknown angular systematics can be mitigated by nulling the density fluctuations on a chosen angular scale. The induced loss of clustering is modelled by an angular integral constraint which can be combined with the radial one.