CORSIKA (COsmic Ray SImulations for KAscade) is a widely-used Monte Carlo simulation program for modeling the development of extensive air showers initiated by high-energy cosmic rays in the atmosphere. It tracks the full particle cascade — from the primary cosmic ray interaction down to secondary particles including muons, electrons, and photons — and supports multiple hadronic interaction models such as SIBYLL, QGSJET, and DPMJET.
In this work, we use CORSIKA to simulate the atmospheric muon flux, adopting the Papini parameterization for the primary proton spectrum. Simulated events are processed through a weighted flux extraction pipeline, in which each muon is assigned a Papini weight corresponding to its parent primary energy, and the results are normalized to physical units via an equivalent exposure time derived from the total simulation weight.
For near-vertical muons (zenith angle 0°–5°), the simulated differential energy spectrum shows good agreement with experimental measurements across the energy range of 0.1–10⁵ GeV, validating the simulation framework and normalization procedure.
Current efforts are focused on large zenith angle muons (70°–80°). Preliminary results indicate a systematic discrepancy between the simulated flux and experimental data in this angular regime. The origin of this discrepancy is under investigation, with particular attention to the Earth's curvature effect — which becomes non-negligible for near-horizontal trajectories — and the geometric configuration of the simulation setup.