The particle showers produced in the atmosphere due to the interactions of primary cosmic
particles require a thorough understanding in the backdrop of searches for rare interactions. In
this work, we made a comprehensive study of air shower simulations using various combinations of
hadronic models and particle transport code of the CORSIKA package. The primary proton and helium
distributions are taken as power law which are scaled to match the measured flux in balloon and
satellite-based experiments at the top of atmosphere. The shower simulation includes production,
transport and decays of secondaries up to the ground level. In this study, we focus on the bulk
of the spectra and particles which is computationally intensive and hence parallel processing of
events is done on computer cluster. We provide a way to normalize the simulation results to be
compared with the ground-based measurements namely, single and multiple muon yields and their
charge ratios as a function of zenith angle and momentum. This provides a basis for comparisons
among the six model combinations used in this study and the differences are outlined. Most of the
hadronic models in CORSIKA produce the bulk ground based measurements fairly well. We use one of
the best model combinations to quantitatively predict the absolute and relative yields of various
particles at ground level as well as their correlations with primaries and with each other. The
leptonic ratios are obtained as a function of energy and zenith angle which are important inputs
for the neutrino oscillation physics.