The purpose of this analysis is to conduct a simulation study and evaluate the performance of a few hadronic interaction models in the production of $K^{\pm}$ charged particles measured by the ALICE experiment in proton-proton collisions at $\sqrt s$ = 900 GeV using modified Hagedorn function
with embedded transverse flow velocity. Simulation models in this analysis are widely used to describe the propagation of cosmic air showers in the atmosphere. Among the employed hadronic interaction models, we present results from the EPOS-LHC, Pythia, QGSJETII-04, and Sibyll models, as they are widely used and able to achieve good agreement with experimental data. Up to about 1.3 GeV/c, EPOS-LHC and Sibyll show a good prediction of the data, while Pythia, QGSJETII-04 underestimate the data over the entire energy range.