The distribution of horizontal subspaces is always Lagrangian not only for Riemannian metrics, but also for the Ehresmann connection associated to Finsler metrics and the Levi-Civita connection of pseudo-Riemannian metrics. It may not be Lagrangian for general sprays though.
See my answer to this question Intuition for Levi-Civita connection via Hamiltonian flows