Measurements of hadron production in 90 GeV/c proton-carbon interactions

This paper presents the multiplicity of neutral and charged hadrons produced in 90 GeV/c proton-carbon interactions from a dataset taken by the NA61/SHINE experiment in 2017. Particle identification via dE/dx was performed for the charged hadrons $\pi^{\pm}, K^{\pm}$ and $p/\bar{p}$; the neutral hadrons $K^{0}_{S}, \Lambda$, and $\bar{\Lambda}$ were identified via an invariant mass analysis of their decays to charged hadrons. Double-differential multiplicity results as a function of laboratory momentum and polar angle are presented for each particle species; these results provide vital constraints on the predicted neutrino beam flux for current and future long-baseline neutrino oscillation experiments.