We present a full kinematic analysis of neutrino-nucleus charged current quasielastic interactions based on the Local Fermi Gas model and the Random Phase Approximation. The model was implemented in the NEUT Monte Carlo framework, which allows us to investigate potentially measurable observables, including hadron distributions. We compare the predictions simultaneously to the most recent T2K and MINERvA charged current (CC) inclusive, CC0 π and transverse kinematic-imbalance variable results. We pursuit a microscopic interpretation of the relevant reaction mechanisms, with the aim to achieving in neutrino oscillation experiments a correct reconstruction of the incoming neutrino kinematics, free of conceptual biasses. Such study is of the utmost importance for the ambitious experimental program which is underway to precisely determine neutrino properties, test the three-generation paradigm, establish the order of mass eigenstates and investigate leptonic CP violation.