Doubly differential transverse momentum spectra of strange mesons [Formula: see text] and [Formula: see text] within rapidity interval of [Formula: see text] are studied by using Monte Carlo event generators in pp, pPb and PbPb collisions at different center-of-mass energies. The predictions of simulations are compared to the ALICE data. The fitting of doubly differential transverse momentum spectra of ALICE data is performed by using Tsallis distribution function. To check the quality of the fitting, chemical potential is also set to zero value. The extracted values of the non-extensive parameter ([Formula: see text]), a chemical potential [Formula: see text] and effective Tsallis temperature [Formula: see text] for different particle species have been obtained from the thermodynamically consistent Tsallis distribution function. It is observed that the effective Tsallis temperature increases at different center-of-mass energies. The predictions of the MC event generators, i.e. PYTHIA8 CR and HIJING 2.0 are compared to the ALICE data, deviations are observed which may be due to the reshuffling of parton colors and string fragmentation.