AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA OOO OOO OOOO OOOO OO O O O O O O O O O O O OO O O O O O O O O O O OO O O O O O O O O O OOOO OO OO O O O O O OOOO O OO O O OOOOOOO O O O O O O O O OO O O O O OOO OOO O O OOOO OO O O O O COSMIC RAY SIMULATION FOR KASCADE A PROGRAM TO SIMULATE EXTENSIVE AIR SHOWERS IN ATMOSPHERE BASED ON A PROGRAM OF P.K.F. GRIEDER, UNIVERSITY BERN, SWITZERLAND QGSJET MODEL ACCORDING TO N.N. KALMYKOV AND S.S. OSTAPCHENKO, MSU, MOSCOW, RUSSIA HDPM MODEL ACCORDING TO J.N. CAPDEVIELLE, COLLEGE DE FRANCE, PARIS, FRANCE GHEISHA ROUTINES ACCORDING TO H. FESEFELDT, RWTH AACHEN, GERMANY EGS4 ACCORDING TO W.R. NELSON, H. HIRAYAMA, D.W.O. ROGERS, SLAC, STANFORD, USA NKG FORMULAS FOR FAST SIMULATION OF EL.MAG. PARTICLES REFERENCES: D. HECK, J.KNAPP, J.N. CAPDEVIELLE, G. SCHATZ, T. THOUW, REPORT FZKA 6019 (1998) SEE ALSO WEB PAGE http://www-ik.fzk.de/corsika/ KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) INSTITUT FUER KERNPHYSIK POSTFACH 3640 D-76021 KARLSRUHE GERMANY IN CASE OF PROBLEMS CONTACT: Dr. Tanguy Pierog e-mail: tanguy.pierog@kit.edu FAX: (49) 7247-82-4075 PHONE: (49) 7247-82-8134 OR : Dr. Dieter Heck e-mail: dieter.heck@kit.edu FAX: (49) 7247-82-4075 PHONE: (49) 7247-82-3777 AND SEND YOUR LIST-FILE BY E-MAIL NUMBER OF VERSION : 6.971 DATE OF VERSION : APRIL 6 2011 VERSION GENERATED FOR UNIX OR COMPATIBLE SYSTEMS ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ (RECL IS DEFINED IN BYTES) WITH NEW DATE_AND_TIME ROUTINE ZENITH ANGLE DEPENDENCE FOR FLAT DETECTOR ARRAY ========== USERS RUN DIRECTIVES FOR THIS SIMULATION =========================== DATA CARDS FOR RUN STEERING ARE EXPECTED FROM STANDARD INPUT RUNNR 1 run number EVTNR 1 number of first shower event NSHOW 100 number of showers to generate PRMPAR 1 particle type of prim. particle ESLOPE -2.7 slope of primary energy spectrum ERANGE 1.E5 1.E5 energy range of primary particle THETAP 0. 0. range of zenith angle (degree) PHIP 0. 0. range of azimuth angle (degree) SEED 1 0 0 seed for 1. random number sequence SEED 2 0 0 seed for 2. random number sequence OBSLEV 0.E2 observation level (in cm) OBSLEV 110.E2 OBSLEV 5000.E2 OBSLEV 1800.E2 FIXCHI 0. starting altitude (g/cm**2) MAGNET 20.0 42.8 magnetic field centr. Europe HADFLG 0 0 0 0 0 2 flags hadr.interact.&fragmentation ECUTS 0.3 0.3 0.003 0.003 energy cuts for particles MUADDI F additional info for muons MUMULT T muon multiple scattering angle ELMFLG T T em. interaction flags (NKG,EGS) STEPFC 1.0 mult. scattering step length fact. RADNKG 200.E2 outer radius for NKG lat.dens.distr. ARRANG 0. rotation of array to north LONGI T 20. T T longit.distr. & step size & fit & out ECTMAP 1.E3 cut on gamma factor for printout MAXPRT 100 max. number of printed events DIRECT ./ output directory DATBAS T write .dbase file PAROUT T F write DAT file USER you user DEBUG F 6 F 1000000 debug flag and log.unit for out EXIT terminates input END OF DATACARD INPUT ========== START OF RUN ======================================================= PRESENT TIME : 03.10.2012 13:51:04 RANDOM NUMBER GENERATOR AT BEGIN OF RUN : SEQUENCE = 1 SEED = 1 CALLS = 0 BILLIONS = 0 SEQUENCE = 2 SEED = 2 CALLS = 0 BILLIONS = 0 PARTICLE OUTPUT TO FILE : ./DAT000001 LONGITUDINAL OUTPUT TO FILE: ./DAT000001.long DBASE OUTPUT TO FILE : ./DAT000001.dbase ========== INTERACTION MODELS ================================================= QGSJET TREATS HIGH ENERGY HADRONIC INTERACTIONS QGSJET CROSS-SECTIONS ARE TAKEN ==================================================== | | | QUARK GLUON STRING JET MODEL | | | | HADRONIC INTERACTION MONTE CARLO | | BY | | N.N. KALMYKOV AND S.S. OSTAPCHENKO | | | | e-mail: serg@eas.npi.msu.su | | | | Publication to be cited when using this program: | | N.N. Kalmykov & S.S. Ostapchenko, A.I. Pavlov | | Nucl. Phys. B (Proc. Suppl.) 52B (1997) 17 | | | | last modification: June 15, 2005 by D. Heck | | (version qgsjet01c.f) | ==================================================== NUCLEUS PROJECTILES FRAGMENT REALISTICALLY IN THE FIRST INTERACTION FRAGMENTATION WITH EVAPORATION (PT AFTER JACEE) GHEISHA TREATS LOW ENERGY HADRONIC INTERACTIONS HIGH ENERGY INTERACTION MODEL USED ABOVE 80.000 GEV LAB ENERGY OR 12.324 GEV CM ENERGY MUON/TAU MULTIPLE SCATTERING AFTER MOLIERE ELECTROMAGNETIC COMPONENT SIMULATED WITH NKG START EGS4 AIR SHOWER SUBROUTINE VERSION ELECTROMAGNETIC COMPONENT SIMULATED WITH EGS4 DATASET EGSDAT6_3. AVAILABLE DATA FOR MEDIUM # 1, WHICH IS: MEDIUM=AIR-NTP ,STERNCID=AIR-NTP MIXT,RHO= 6.1250E-09,NE= 4,COMPOSITION IS : ASYM=N ,Z= 7.,A= 14.007,PZ= 5.39228E-02,RHOZ= 7.55280E-01 ASYM=O ,Z= 8.,A= 15.999,PZ= 1.44793E-02,RHOZ= 2.31660E-01 ASYM=AR,Z=18.,A= 39.948,PZ= 3.20917E-04,RHOZ= 1.28200E-02 ASYM=C ,Z= 6.,A= 12.011,PZ= 3.33024E-06,RHOZ= 4.00000E-05 EGS SUCCESSFULLY 'HATCHED' FOR ONE MEDIUM. ELECTRONS CAN BE CREATED AND ANY ELECTRON FOLLOWED DOWN TO 0.00300 GEV KINETIC ENERGY GAMMAS CAN BE CREATED AND ANY GAMMA FOLLOWED DOWN TO 0.00300 GEV ENERGY ELECTRON HISTORIES ARE TERMINATED AT 0.00351 GEV GAMMA HISTORIES ARE TERMINATED AT 0.00300 GEV ========== CONSTANTS AND PARAMETERS =========================================== PHYSICAL CONSTANTS C(1) TO C(50) 6.37131500E+08 6.00000000E+05 2.00000000E+06 0.00000000E+00 0.00000000E+00 4.58059628E-02 5.73089600E-01 5.28304202E-02 2.50000000E+00 2.07000000E+00 8.20000000E+00 1.00000000E-01 0.00000000E+00 0.00000000E+00 1.00002339E+00 9.67266341E-03 1.00000008E+00 5.75128639E-04 0.00000000E+00 0.00000000E+00 3.77000000E+01 1.53287300E-04 9.38641700E+00 2.00000000E-03 2.99792458E+10 1.00000000E+00 5.40302306E-01 1.57000000E+00 7.96326711E-04 2.10000000E-02 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.00000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.61909592E-01 8.99834237E-01 0.00000000E+00 1.03899199E+00 2.71383397E-01 1.37036000E+02 CONSTANTS FOR KAONS CKA(1) TO CKA(80) 0.00000000E+00 1.00000000E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 2.50000000E-01 5.00000000E-01 7.50000000E-01 1.00000000E+00 5.00000000E-01 2.00000000E-01 9.26406642E-01 1.12407138E+00 1.49600000E+02 1.49600000E+02 2.35531825E-01 2.06000000E-01 1.35000000E-01 2.22000000E-01 5.00000000E-01 0.00000000E+00 6.34300000E-01 6.89500000E-01 8.73700000E-01 6.62400000E-01 3.89500000E-01 2.44320613E+00 9.12400647E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.00000000E+00 1.00000000E+05 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 8.45600000E-01 9.01400000E-01 9.50100000E-01 9.82800000E-01 -2.15400000E-01 1.20000000E-02 -1.01000000E-02 1.27000000E+00 6.38000000E-01 5.70000000E-02 0.00000000E+00 1.84000000E+00 0.00000000E+00 1.00000000E+00 6.78000000E-01 7.60000000E-02 9.90000000E-03 2.22000000E+00 2.77000000E-02 0.00000000E+00 1.28800000E-02 2.78000000E-02 1.74000000E-02 1.19400000E-02 2.91000000E-02 0.00000000E+00 1.31000000E-02 3.30000000E-02 2.70000000E-02 1.24100000E-02 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 CONSTANTS FOR ETAS CETA(1) TO CETA(5) 3.94600000E-01 7.20000000E-01 9.49600000E-01 -1.07000000E+00 2.07000000E+00 CONSTANTS FOR STRANGE BARYONS...
inferno2302