Bayesian luminosity function estimation in multi-depth datasets with selection effects: A case study for $3<z<5$ Lyman $\alpha$ emitters

TL;DR

This paper introduces a hierarchical Bayesian method to accurately estimate luminosity functions from multiple surveys with different depths and coverage, explicitly accounting for selection effects and uncertainties, demonstrated on high-redshift Lyman alpha emitters.

Contribution

The paper presents a novel Bayesian framework that combines data from diverse surveys to infer luminosity functions while explicitly modeling selection effects and measurement uncertainties.

Findings

Deep data are crucial for reducing biases at the faint end.

Wide-area data help constrain the bright end of the luminosity function.

The method provides consistent parameter estimates with previous lensing studies.

Abstract

We present a hierarchical Bayesian framework designed to infer the luminosity function of any class of object by jointly modelling data from multiple surveys with varying depth, completeness, and sky coverage. Our method explicitly accounts for selection effects and measurement uncertainties (e.g. in luminosity) and can be generalized to any extensive quantity, such as mass. We validated the model using mock catalogues; from this we determined that deep data reaching $\gtrsim 1.5$ dex below a characteristic luminosity ($\tilde{L}^\star$) are essential to reducing biases at the faint end ($\lesssim 0.1$ dex) and that wide-area data help constrain the bright end. As a proof of concept, we considered a combined sample of 1176 Lyman $\alpha$ emitters at redshift $3 < z < 5$ drawn from several MUSE surveys, ranging from ultra-deep ($\gtrsim 90$ hr) and narrow ($\lesssim 1$ arcmin$^2$) fields to shallow ($\lesssim 5$ hr) and wide ($\gtrsim 20$ arcmin$^2$) fields. With this complete sample, we constrain the luminosity function parameters $\log(\Phi^\star/\mathrm{Mpc^{-3}}) = -2.86^{+0.15}_{-0.17}$, $\log(L^\star/\mathrm{erg\,s^{-1}}) = 42.72^{+0.10}_{-0.09}$, and $\alpha = -1.81^{+0.09}_{-0.09}$, where the uncertainties represent the $90\%$ credible intervals. These values are in agreement with the results of studies based on gravitational lensing that reach $\log(L/\mathrm{erg\,s^{-1}}) \approx 41$, although differences in the faint-end slope underscore how systematic errors are starting to dominate. In contrast, wide-area surveys represent the natural extension needed to constrain the brightest Lyman $\alpha$ emitters [$\log(L/\mathrm{erg\,s^{-1}}) \gtrsim 43$], where statistical uncertainties still dominate.