پرتوهای کیهانی ذراتی پر انرژی هستند که از منابع گوناگونی سرچشمه می گیرند. کشف منابع این ذرات از مهترین مسائل موجود در فیزیک نجوم می باشد. این ذرات هنگام ورود به اتمسفر زمین به علت اندرکنش هایی که با هسته های هوا دارند، بهمنی از ذرات ثانویه را تشکیل می دهند. ذرات ثانویه تولید شده می توانند توسط آرایه های عظیمی از آشکارسازها که در روی زمین برپا شده آشکار سازی شوند. در کنار این آشکارسازی تمام اندرکنش ها و واپاشی هایی که در طول توسعه بهمن های گسترده هوایی (eas) رخ می دهد، می توانند توسط مدل های اندرکنش هادرونی که در کد corsika اجرا می شوند شبیه سازی شوند. این مدل ها دارای ویژگی های مختلفی هستند و نتایج متفاوتی را ارائه می دهند. با توجه به نتایج این شبیه سازی ها و فراوانی های بدست آمده از آشکارسازی ها می توان ترکیب جرمی و انرژی اولیه پرتوهای کیهانی را پیش بینی کرد. حالت مطلوب رسیدن به یک توافق سازگار بین مدل ها و نتایج تجربی می باشد و به همین منظور مدل ها روز به روز بهبود می یابند. در این کار ما بهمن های تولید شده توسط ذره اولیه پروتون که به وسیله سه مدل اندرکنش هادرونی انرژی پایین (gheisha, urqmd, fluka) شبیه سازی شده اند را مورد بررسی دقیق قرار دادیم.