The halo mass function is a very important probe for testing cosmological models and our theory of structure formation. However, the only analytical deduction of the halo mass function (apart from fitting formulas) is based on the linear evolution of the cosmic density field and a linear collapse threshold derived from spherical collapse. Since both of these assumptions are not ideal, I will present an analytical estimate of the HMF from the non-linearly evolved density field. Our approach is based on excursion set statistics with correlated steps, which we combine with a PDF model of the cosmic density field. The parameters of this PDF model are fixed by the Kinetic Field Theory (KFT) of cosmic structure formation. This yields a closed-form expression of the HMF that only depends on the overdensity threshold. Treating this threshold as a free parameter, we find very good agreement with measurements from numerical simulations over the mass range from 10^10 to 10^16 solar masses.