Mean field approach to reconstructed neutrino energy distributions in accelerator-based experiments

TL;DR

This paper develops a mean field theoretical approach using HF-CRPA to model neutrino energy reconstruction in accelerator experiments, revealing significant differences from traditional models and impacting oscillation analysis.

Contribution

It introduces a mean field approach with HF-CRPA for neutrino energy reconstruction, highlighting effects of nuclear dynamics often neglected in simpler models.

Findings

HF-CRPA shows low reconstructed energy tails due to elastic distortion.

Nuclear dynamics reshape energy distributions beyond plane wave models.

RFG cannot reproduce HF-CRPA results with any binding energy value.

Abstract

The reconstruction of the neutrino energy is crucial in oscillation experiments that use interactions with nuclei to detect the neutrino. The common reconstruction procedure is based on the kinematics of the final-state lepton. The interpretation of the reconstructed energy in terms of the real neutrino energy must rely on a model for the neutrino-nucleus interaction. The Relativistic Fermi Gas (RFG) model is frequently used in these analyses. In the Hartree-Fock (HF) model for quasielastic nucleon knockout, the bound nucleon wave functions are obtained using an effective nucleon-nucleon force. The final-state wave function is constructed from continuum states in the same potential which have the correct asymptotic behavior. The Continuum Random Phase Approximation (CRPA) model extends the HF approach taking long range correlations into account in a self-consistent way. Considering only single-nucleon processes, the distributions of reconstructed neutrino energies obtained within the HF-CRPA approach are compared with the results of the RFG, an RPWIA calculation, and the RPA+np-nh model of Martini et al. We find that the distributions of reconstructed energies for a fixed incoming energy in the HF-CRPA display additional strength in the low reconstructed energy tails compared to models without elastic distortion of the outgoing nucleon. This asymmetry redistributes strength from higher to lower values of the reconstructed energy. The mean field description of the nuclear dynamics results in a reshaping of the reconstructed energy distribution that cannot be accounted for in a plane wave impulse approximation model, even by modifying ad hoc parameters such as the binding energy. In particular it is shown that in the RFG calculations there is no value of the binding energy which is able to reproduce the entire T2K $\nu_\mu$ oscillated spectrum as calculated in HF-CRPA.