  In this article, we obtain hessian estimates for Kolmogorov-Fokker-Planck
operators in non-divergence form in several Banach function spaces. Our
approach relies on a representation formula and newly developed sparse
domination techniques in Harmonic Analysis. Our result when restricted to
weighted Lebesgue spaces yields sharp quantitative hessian estimates for the
Kolmogorov-Fokker-Planck operators.
