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

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{Mp}{\mathrm{c}}^{-3})=-2.86_{-0.17}^{+0.15}$, $\log (L^{\star }/\mathrm{erg}{\mathrm{s}}^{-1})=42.72_{-0.09}^{+0.10}$, 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}{\mathrm{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}{\mathrm{s}}^{-1})\gtrsim 43$], where statistical uncertainties still dominate.

11 pages, 3 figures, 3 tables (in the main text). Published in A&A

References in corpus (38)

• Measuring Reddening with SDSS Stellar Spectra and Recalibrating SFD
• The Galaxy Luminosity Function and Luminosity Density at Redshift z=0.1
• The Evolution of the Stellar Mass Functions of Star-Forming and Quiescent Galaxies to z = 4 from the COSMOS/UltraVISTA Survey
• The Hubble Ultra Deep Field
• The MUSE second-generation VLT instrument
• The SDSS Quasar Survey: Quasar Luminosity Function from Data Release Three
• The 2dF QSO Redshift Survey - I. The Optical QSO Luminosity Function
• New Determinations of the UV Luminosity Functions from z~9 to z~2 show a remarkable consistency with halo growth and a constant star formation efficiency
• Extracting distribution parameters from multiple uncertain observations with selection biases
• High-Redshift Quasars Found in Sloan Digital Sky Survey Commissioning Data IV: Luminosity Function from the Fall Equatorial Stripe Sampl
• Ubiquitous giant Ly $\alpha$ nebulae around the brightest quasars at $z\sim 3.5$ revealed with MUSE
• Extended Lyman alpha haloes around individual high-redshift galaxies revealed by MUSE
• The Luminosity Function of Fermi-detected Flat-Spectrum Radio Quasars
• Observations of the Lyman-$\alpha$ Universe
• QSO MUSEUM I: A sample of 61 extended Ly$\alpha$-emission nebulae surrounding $z\sim 3$ quasars
• Soft X-ray AGN Luminosity Function from ROSAT Surveys II. Table of the binned Soft X-ray Luminosity Function
• The MUSE Hubble Ultra Deep Field surveys: Data release II
• Accounting for Source Uncertainties in Analyses of Astronomical Survey Data
• The MUSE Hubble Ultra Deep Field Survey VI: The Faint-End of the Lya Luminosity Function at 2.91 < z < 6.64 and Implications for Reionisation
• Mining Gravitational-wave Catalogs To Understand Binary Stellar Evolution: A New Hierarchical Bayesian Framework
• The Las Campanas Infrared Survey. IV. The Photometric Redshift Survey and the Rest-frame R-band Galaxy Luminosity Function at 0.5 <= z <= 1.5
• The MUSE Hubble Ultra Deep Field Survey XIII. Spatially resolved spectral properties of Lyman alpha haloes around star-forming galaxies at z > 3
• MUSE Analysis of Gas around Galaxies (MAGG) -- III: The gas and galaxy environment of z = 3-4.5 quasars
• Shedding Light on the Galaxy Luminosity Function
• The MUSE Extremely Deep Field: the Cosmic Web in Emission at High Redshift
• The MUSE-Wide Survey: A determination of the Lyman $\alpha$ emitter luminosity function at $3<z<6$
• Mass Functions of Supermassive Black Holes Across Cosmic Time
• LSDCat: Detection and cataloguing of emission-line sources in integral-field spectroscopy datacubes
• MUSE Analysis of Gas around Galaxies (MAGG) -- I: Survey design and the environment of a near pristine gas cloud at z~3.5
• The MUSE Ultra Deep Field (MUDF). II. Survey design and the gaseous properties of galaxy groups at 0.5 < z < 1.5
• A Flexible Method of Estimating Luminosity Functions
• Predicting the redshift 2 Halpha luminosity function using [OIII] emission line galaxies
• The MUSE Ultra Deep Field (MUDF). I. Discovery of a group of Ly$\alpha$ nebulae associated with a bright $z\approx 3.23$ quasar pair
• Faint end of the $z\sim 3-7$ luminosity function of Lyman-alpha emitters behind lensing clusters observed with MUSE
• J-PLUS: Unveiling the brightest-end of the Lyα luminosity function at 2.0<z<3.3 over 1000 deg^2
• First HETDEX Spectroscopic Determinations of Ly$\alpha$ and UV Luminosity Functions at $z=2-3$: Bridging a Gap Between Faint AGN and Bright Galaxies
• Probing the faint end Luminosity Function of Lyman Alpha Emitters at 3<z<7 behind 17 MUSE lensing clusters
• ODIN: Clustering Analysis of 14,000 Lyα Emitting Galaxies at z=2.4, 3.1, and 4.5