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-II MODEL ACCORDING TO S.S. OSTAPCHENKO, IEKP, KARLSRUHE AND 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) D. HECK, J. KNAPP, REPORT FZKA 6097 (1998) D. HECK, REPORT FZKA 7254 (2006) 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) 721-608-24075 PHONE: (49) 721-608-28134 OR : Dr. Dieter Heck e-mail: dieter.heck@kit.edu FAX: (49) 721-608-24075 PHONE: (49) 721-608-23777 AND SEND YOUR LIST-FILE BY E-MAIL NUMBER OF VERSION : 7.4001 DATE OF VERSION : JULY 22 2014 VERSION GENERATED FOR UNIX OR COMPATIBLE SYSTEMS ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ (RECL IS DEFINED IN BYTES) WITH TIMERC DATE ROUTINE ZENITH ANGLE DEPENDENCE FOR FLAT DETECTOR ARRAY COASTUSERLIB IS USED ==================== CONEX VERSION WITH CASCADE EQUATIONS ==================================== WARNING: DESPITE A LOT OF SUCCESSFUL TESTS, THE AUTHORS COULD NOT TEST ALL POSSIBLE COMBINATIONS OF MODELS AND OBSERVABLE TYPES. FOR THE MOMENT THE PRECISION EXPECTED BY THE USE OF CONEX OPTION SHOULD BE BETTER THAN 10% COMPARED TO USUAL CORSIKA FOR A FACTOR OF 10 GAIN IN CPU TIME BUT IT MAY DEPEND ON THE INITIAL CONDITIONS AND OBSERVABLES. AS A CONSEQUENCE THEY DECLINE ANY RESPONSABILITY IF AN ANALYSIS IS BIASED BY THE USE OF CONEX OPTION. THE USER SHOULD VALIDATE HIS ANALYSIS BY THE USE OF CORSIKA WITHOUT CONEX FOR PARTICLE BASED ANALYSIS. 1D SIMULATIONS BASED ON CONEX ARE ALREADY USED INTENSIVELY SINCE MANY YEARS. ==================================== CURVED VERSION WITH SLIDING PLANAR ATMOSPHERE ============================================= SLANT DEPTH FOR LONGITUDINAL DISTRIBUTIONS ========================================== THINNING IS ACTIVE ================== UPWARD VERSION FOR UPWARD GOING PARTICLES ========================================= ========== USERS RUN DIRECTIVES FOR THIS SIMULATION =========================== DATA CARDS FOR RUN STEERING ARE EXPECTED FROM STANDARD INPUT RUNNR 1 EVTNR 1 SEED 1 0 0 SEED 2 0 0 SEED 3 0 0 PRMPAR 14 ERANGE 1.000E+8 1.000E+8 ESLOPE 0.000E+00 THETAP 0. 0. PHIP 0.000E+00 0.000E+00 ECUTS 3.000E-01 3.000E-01 4.010E-04 4.010E-04 ELMFLG T T THIN 1.000E-06 1.000E+02 0.000E+00 THINH 1.000E+00 1.000E+02 NSHOW 1 USER huege HOST iklxds69 DIRECT './' OBSLEV 140000.0 ECTMAP 1.000E+05 STEPFC 1.000E+00 MUMULT T MUADDI T PAROUT F F MAXPRT 1 MAGNET 19.71 -14.18 LONGI T 5. T T RADNKG 5.000E+05 DATBAS F EXIT END OF DATACARD INPUT PARTICLES REGISTERED ON CURVED OBSLEV (X,Y) ARE NOT CARTESIAN COORDINATES! PLEASE READ THE USERS GUIDE SEE KEYWORD: CURVOUT/FLATOUT ===================================== ========== START OF RUN ======================================================= PRESENT TIME : 08.10.2014 07:24:12 UTC RANDOM NUMBER GENERATOR AT BEGIN OF RUN : SEQUENCE = 1 SEED = 1 CALLS = 0 BILLIONS = 0 SEQUENCE = 2 SEED = 2 CALLS = 0 BILLIONS = 0 SEQUENCE = 3 SEED = 3 CALLS = 0 BILLIONS = 0 SEQUENCE = 4 SEED = 4 CALLS = 0 BILLIONS = 0 SEQUENCE = 5 SEED = 5 CALLS = 0 BILLIONS = 0 SEQUENCE = 6 SEED = 6 CALLS = 0 BILLIONS = 0 SEQUENCE = 7 SEED = 3 CALLS = 0 BILLIONS = 0 SEQUENCE = 8 SEED = 5268784 CALLS = 0 BILLIONS = 0 ******************************************************* ** ** You are using the CoREAS V1 radio simulation code ** ** THIN/CURVED/SLANT/STACKIN/PRESHOWER: 1/1/1/0/0 ** ** Please cite the following article: ** T. Huege, M. Ludwig, C.W. James, AIP Conf. Proc. ** 1535, 128-132 (2013), doi:10.1063/1.4807534 ** ******************************************************* LONGITUDINAL OUTPUT TO FILE: ./DAT000001.long CONEX IS USED WITH conex_qgsjetII_gheisha.param ========== INTERACTION MODELS ================================================= QGSJET TREATS HIGH ENERGY HADRONIC INTERACTIONS QGSJET CROSS-SECTIONS ARE TAKEN ==================================================== | | | QUARK GLUON STRING JET -II MODEL | | | | HADRONIC INTERACTION MONTE CARLO | | BY | | S. OSTAPCHENKO | | | | e-mail: sergei@tf.phys.ntnu.no | | | | Version II-04 | | | | Publication to be cited when using this program: | | S.Ostapchenko, PRD 83 (2011) 014018 | | | | last modification: 09.04.2013 | | | | Any modification has to be approved by the author| ==================================================== qgaini: cross sections readout from the file: qgsdat-II-04 done qgaini: nuclear cross sections readout from the file sectnu-II-04 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 ###################################################### # SIMULATION WITH NKG NOT POSSIBLE IN CURVED VERSION # ###################################################### START EGS4 AIR SHOWER SUBROUTINE VERSION ELECTROMAGNETIC COMPONENT SIMULATED WITH EGS4 DATASET EGSDAT6_.4 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.00040 GEV KINETIC ENERGY GAMMAS CAN BE CREATED AND ANY GAMMA FOLLOWED DOWN TO 0.00040 GEV ENERGY ELECTRON HISTORIES ARE TERMINATED AT 0.00091 GEV GAMMA HISTORIES ARE TERMINATED AT 0.00040 GEV ========== CONSTANTS AND PARAMETERS =========================================== PHYSICAL CONSTANTS C(1) TO C(50) 6.37131500E+08 6.00000000E+05 2.00000000E+06 -1.35122092E+03 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 3.14159265E+00 -1.00000000E+00 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 CSTRBA(1) TO CSTRBA(11) 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.48200000E-01 5.16300000E-01 0.00000000E+00 0.00000000E+00 0.00000000E+00 6.78000000E-01 9.14000000E-01 ========== ATMOSPHERE ========================================================= ( US STANDARD ATMOSPHERE PARAMETERIZED BY LINSLEY ) HEIGHT H IN KM GIVES THICKNESS OF ATMOSPHERE T IN G/CM**2 H = 0.0... 4.0 KM ---> T = -1.86556E+02 + 1.2227E+03 * EXP( -H / 9.9419E+00) H = 4.0... 10.0 KM ---> T = -9.49190E+01 + 1.1449E+03 * EXP( -H / 8.7815E+00) H = 10.0... 40.0 KM ---> T = 6.12890E-01 + 1.3056E+03 * EXP( -H / 6.3614E+00) H = 40.0...100.0 KM ---> T = 0.00000E+00 + 5.4018E+02 * EXP( -H / 7.7217E+00) H = 100.0...112.8 KM ---> T = 1.12829E-02 - H / 1.0000E+04 ========== SHOWER PARAMETERS ================================================== PRIMARY PARTICLE IDENTIFICATION IS 14 PRIMARY ENERGY IS FIXED AT 1.000E+08 GEV THETA OF INCIDENCE IS FIXED TO 0.00 DEGREES PHI OF INCIDENCE IS FIXED TO 0.00 DEGREES HEIGHT OF FIRST INTERACTION IS CHOSEN RANDOMLY TARGET OF FIRST INTERACTION IS AT RANDOM STARTING ALTITUDE AT 11282920.00 CM (= 0.00 G/CM**2) WHICH IS AT TOP OF ATMOSPHERE OBSERVATION LEVEL # IN CM AND IN G/CM**2 1 1.40000000E+05 8.75500211E+02 ENERGY FRACTION FOR EM THINNING = 1.0000E-06 AND FOR HADRONIC THINNING = 1.0000E-06 HAD.WEIGHTS ARE LIMITED TO MAX. 1. EM WEIGHTS ARE LIMITED TO MAX. 100. PARTICLES WITH LORENTZ FACTOR LARGER THAN 1.0000E+05 ARE PRINTED OUT SHOWER PARTICLES ENERGY CUT : FOR HADRONS : 3.0000E-01 GEV FOR MUONS : 3.0000E-01 GEV FOR ELECTRONS : 4.0100E-04 GEV FOR GAMMAS : 4.0100E-04 GEV EARTH MAGNETIC FIELD STRENGTH IS 24.2807846 MICROTESLA WITH INCLINATION ANGLE -35.7324409 DEGREES ADDITIONAL INFORMATION ON MUON ORIGIN IS WRITTEN TO PARTICLE DATA FILE NUMBER OF SHOWERS TO GENERATE = 1 initialize MC ... model IIqgsjet lemodel gheisha !------------------------------ fconex wle /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.wghe20 ! gheisha fconex whe /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.wIIq20 ! qgsjet01 fconex dkz /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dkz20 ! decays fconex dkl /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dkl20 fconex dks /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dks20 fconex dkm /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dkm20 fconex dke /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dke20 fconex dkn /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dkn20 fconex dkg /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.dkg20 fconex p2le /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2ghe20 fconex p2he /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2IIq20 fconex p2d /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2dec20 !fopen ! to use output file (check, histo ...) !------------------------------ fegs dat /home/grubio/BUILD/CORSIKA-74001/build/conex/egs4.dat fegs out /dev/null !------------------------------ set ixmax 0 ! do not fit profile with G.H. set hacut1 0.1 ! cut for hadrons in GeV (not less than 0.1) set emcut1 0.001 ! cut for leptons in GeV (not less than 0.001) set enymin 0.1 ! minimum energy for gheisha initialize EGS4 ... CALL HATCH TO GET CROSS-SECTION DATA DUNIT REQUESTED&USED ARE: -2.00000E+00 8.16327E+08(CM.) EGS SUCCESSFULLY 'HATCHED' FOR ONE MEDIUM. ANY ELECTRON CAN BE FOLLOWED DOWN TO 0.200 MeV KINETIC ENERGY ANY PHOTON CAN BE FOLLOWED DOWN TO 0.200 MeV Full mode => minEhad : 14 0.000199526231 initialize QGSJET-II ... initialize Gheisha ... read LE model table from /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.wghe20 read LE pt2 table from /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2ghe20 cross section tables read HE model table from /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.wIIq20 read HE pt2 table from /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2IIq20 cross section tables read decay pt2 table from /home/grubio/BUILD/CORSIKA-74001/build/conex/conex.p2dec20 CONEXINI: GIVEN THSTEP NOT COMPATIBLE WITH CONEX USE 10. INSTEAD