6. Conclusions
We introduced two different methods to reconstruct the lateral distribution function of air shower muons: the profile and the integrated likelihoods. Both likelihoods extend a previous approach by considering the detector timing. Although we applied the likelihoods to a specific cosmic ray detector, they can be used for any kind of segmented particle counters with time resolution. We found an optimal distance of 450 m to measure the shower size parameter in a triangular array with 750 m between detectors. The new likelihoods improve the reconstruction in two aspects. Firstly, by raising the number of muons a detector can handle before saturating, more events can be reconstructed. The recovery is more significant close to 10 EeV, the upper limit of the considered energy range, a region where events are usually scarce. Secondly, we reduced the statistical fluctuations of the parameter that measures the shower size from 1 EeV upwards. This decrease allows for a more powerful discrimination between different primary masses based on the number of muons. By comparing to an ideal muon counter, we established that the resolutions achieved with the new likelihoods are close to the lower bound given the detector size and spacing. We also showed that the approximations introduced for the profile and integrated likelihoods do not bias the reconstructed shower size parameter and kept the coverage of its 1σ confidence interval close to the expected Gaussian nominal value. The shower size parameters reconstructed with the integrated and the profile likelihoods are very similar. Nevertheless the profile likelihood is the preferred reconstruction method given the much shorter time it takes to process the data. The correspondence between the profile and the integrated likelihood results, shows the robustness of these techniques to reconstruct the muon lateral distribution with an array of segmented counters.