| 894 | | TOTAL SPACE ON MPATAP USED = 6552 WORDS |
| 895 | | SPACE PER EVENT ON MPATAP = 6552 WORDS |
| 896 | | AVERAGE HEIGHT OF 1ST INT. = 26.809 +- 0.000 G/CM**2 |
| 897 | | |
| | 894 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
| | 895 | |
| | 896 | OOO OOO OOOO OOOO OO O O O |
| | 897 | O O O O O O O O OO O O O O |
| | 898 | O O O O O O OO O O O O |
| | 899 | O O O O O OOOO OO OO O O |
| | 900 | O O O OOOO O OO O O OOOOOOO |
| | 901 | O O O O O O O O OO O O O O |
| | 902 | OOO OOO O O OOOO OO O O O O |
| | 903 | |
| | 904 | COSMIC RAY SIMULATION FOR KASCADE |
| | 905 | |
| | 906 | A PROGRAM TO SIMULATE EXTENSIVE AIR SHOWERS IN ATMOSPHERE |
| | 907 | |
| | 908 | BASED ON A PROGRAM OF P.K.F. GRIEDER, UNIVERSITY BERN, SWITZERLAND |
| | 909 | QGSJET MODEL ACCORDING TO N.N. KALMYKOV AND S.S. OSTAPCHENKO, MSU, MOSCOW, RUSSIA |
| | 910 | HDPM MODEL ACCORDING TO J.N. CAPDEVIELLE, COLLEGE DE FRANCE, PARIS, FRANCE |
| | 911 | GHEISHA ROUTINES ACCORDING TO H. FESEFELDT, RWTH AACHEN, GERMANY |
| | 912 | EGS4 ACCORDING TO W.R. NELSON, H. HIRAYAMA, D.W.O. ROGERS, SLAC, STANFORD, USA |
| | 913 | NKG FORMULAS FOR FAST SIMULATION OF EL.MAG. PARTICLES |
| | 914 | |
| | 915 | REFERENCES: D. HECK, J.KNAPP, J.N. CAPDEVIELLE, G. SCHATZ, T. THOUW, |
| | 916 | REPORT FZKA 6019 (1998) |
| | 917 | D. HECK, J. KNAPP, REPORT FZKA 6097 (1998) |
| | 918 | SEE ALSO WEB PAGE https://www.ikp.kit.edu/corsika/ |
| | 919 | KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) |
| | 920 | INSTITUT FUER KERNPHYSIK |
| | 921 | POSTFACH 3640 |
| | 922 | D-76021 KARLSRUHE |
| | 923 | GERMANY |
| | 924 | |
| | 925 | IN CASE OF PROBLEMS CONTACT: Dr. Tanguy Pierog |
| | 926 | e-mail: tanguy.pierog@kit.edu |
| | 927 | FAX: (49) 721-608-24075 |
| | 928 | PHONE: (49) 721-608-28134 |
| | 929 | OR : Dr. Dieter Heck |
| | 930 | e-mail: dieter.heck@partner.kit.edu |
| | 931 | FAX: (49) 721-608-24075 |
| | 932 | PHONE: (49) 721-608-23777 |
| | 933 | AND SEND YOUR LIST-FILE BY E-MAIL |
| | 934 | |
| | 935 | NUMBER OF VERSION : 7.6900 |
| | 936 | DATE OF VERSION : DECEMBER 20, 2018 |
| | 937 | |
| | 938 | VERSION GENERATED FOR UNIX OR COMPATIBLE SYSTEMS |
| | 939 | ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ |
| | 940 | (RECL IS DEFINED IN BYTES) |
| | 941 | WITH TIMERC DATE ROUTINE |
| | 942 | ZENITH ANGLE DEPENDENCE FOR FLAT DETECTOR ARRAY |
| 899 | | PRESENT TIME : 05.06.2019 10:39:26 UTC |
| 900 | | |
| 901 | | ========== END OF RUN ================================================ |
| | 944 | INTERFACE FOR EXTERNAL ATMOSPHERIC PROFILES ENABLED |
| | 945 | =================================================== |
| | 946 | |
| | 947 | CHERENKOV RADIATION IS GENERATED |
| | 948 | ================================ |
| | 949 | |
| | 950 | WITH WAVELENGTH DEPENDENCE |
| | 951 | ATMOSPHERIC ABSORPTION, MIRROR REFLECTIVITY AND QUANTUM EFFICIENCY MAY BE RESPECTED |
| | 952 | |
| | 953 | TRAJECTORY VERSION TO FOLLOW EMITTING SOURCE |
| | 954 | ============================================ |
| | 955 | |
| | 956 | THINNING IS ACTIVE |
| | 957 | ================== |
| | 958 | |
| | 959 | PRIMARY DIRECTION IS SELECTED FROM VIEWING CONE |
| | 960 | =============================================== |
| | 961 | |
| | 962 | |
| | 963 | ========== USERS RUN DIRECTIVES FOR THIS SIMULATION =========================== |
| | 964 | |
| | 965 | DATA CARDS FOR RUN STEERING ARE EXPECTED FROM STANDARD INPUT |
| | 966 | |
| | 967 | RUNNR 103 run number |
| | 968 | EVTNR 0 number of first shower event |
| | 969 | NSHOW 1 number of showers to be generated in one run |
| | 970 | |
| | 971 | SEED 62748 0 0 seed for hadronic part |
| | 972 | SEED 62749 0 0 seed for EGS4 part |
| | 973 | SEED 62750 0 0 seed for Cherenkov part |
| | 974 | SEED 62751 0 0 seed for 4. random number sequence |
| | 975 | SEED 62752 0 0 seed for 5. random number sequence |
| | 976 | |
| | 977 | PRMPAR 1 primary particle code: 1 = photon, 3 = e-, 14 = p+ |
| | 978 | FIXCHI 0 starting altitude of primary particle? |
| | 979 | TSTART T taken from IACT option (start emission at border of atmosphere) |
| | 980 | |
| | 981 | ERANGE 10000 10000 energy range of primary particle in GeV |
| | 982 | ESLOPE -2 slope of energy spectrum |
| | 983 | |
| | 984 | PHIP 0. 0. range of azimuth angle in deg |
| | 985 | THETAP 0. 0. range of zenith angle in deg |
| | 986 | VIEWCONE 0.0 0.0 range of viewcone in deg |
| | 987 | |
| | 988 | TRAFLG F TRAJECTORY option |
| | 989 | * SRCPOS ra dec source position |
| | 990 | * TRATM y m d h m s duration observation time |
| | 991 | * TLAT d m s N observatory latitude |
| | 992 | * TLONG d m s N observatory longitude |
| | 993 | |
| | 994 | Atmospheric profile 7 with 50 levels read from file atmprof7.dat |
| | 995 | |
| | 996 | |
| | 997 | Results of the atmosphere fit: |
| | 998 | Layer 1: 0.00 km < h < 3.00 km: a = -144.838, b = 1192.34, c = 994186 |
| | 999 | Layer 2: 3.00 km < h < 10.00 km: a = -124.071, b = 1173.98, c = 967530 |
| | 1000 | Layer 3: 10.00 km < h < 50.00 km: a = 0.360027, b = 1412.08, c = 636143 |
| | 1001 | Layer 4: 50.00 km < h < 100.00 km: a = -0.000824761, b = 810.682, c = 735640 |
| | 1002 | Layer 5: 100.00 km < h < 109.23 km: a = 0.00221589, b = 1, c = 4.92961e+09 |
| | 1003 | |
| | 1004 | |
| | 1005 | Altitude [km] rho(table) rho(fit) thick(table) thick(fit) |
| | 1006 | 0.0 1.19820e-03 1.19931e-03 1.04750e+03 1.04750e+03 |
| | 1007 | 1.0 1.08430e-03 1.08455e-03 9.33510e+02 9.33403e+02 |
| | 1008 | 2.0 9.81420e-04 9.80765e-04 8.30320e+02 8.30225e+02 |
| | 1009 | 3.0 8.88120e-04 8.89886e-04 7.36920e+02 7.36919e+02 |
| | 1010 | 4.0 8.03260e-04 8.02504e-04 6.52420e+02 6.52375e+02 |
| | 1011 | 5.0 7.25830e-04 7.23703e-04 5.76020e+02 5.76133e+02 |
| | 1012 | 6.0 6.54990e-04 6.52639e-04 5.07040e+02 5.07377e+02 |
| | 1013 | 7.0 5.90060e-04 5.88554e-04 4.44830e+02 4.45372e+02 |
| | 1014 | 8.0 5.30430e-04 5.30762e-04 3.88840e+02 3.89456e+02 |
| | 1015 | 9.0 4.75610e-04 4.78644e-04 3.38570e+02 3.39031e+02 |
| | 1016 | 10.0 4.25210e-04 4.60898e-04 2.93560e+02 2.93557e+02 |
| | 1017 | 11.0 3.78710e-04 3.93854e-04 2.53390e+02 2.50907e+02 |
| | 1018 | 12.0 3.35320e-04 3.36562e-04 2.17710e+02 2.14461e+02 |
| | 1019 | 13.0 2.94600e-04 2.87604e-04 1.86240e+02 1.83317e+02 |
| | 1020 | 14.0 2.56450e-04 2.45768e-04 1.58720e+02 1.56703e+02 |
| | 1021 | 15.0 2.21010e-04 2.10017e-04 1.34870e+02 1.33961e+02 |
| | 1022 | 16.0 1.88580e-04 1.79467e-04 1.14420e+02 1.14527e+02 |
| | 1023 | 17.0 1.59740e-04 1.53361e-04 9.70320e+01 9.79196e+01 |
| | 1024 | 18.0 1.34720e-04 1.31052e-04 8.23400e+01 8.37281e+01 |
| | 1025 | 19.0 1.13400e-04 1.11989e-04 6.99620e+01 7.16010e+01 |
| | 1026 | 20.0 9.54910e-05 9.56985e-05 5.95440e+01 6.12380e+01 |
| | 1027 | 21.0 8.05660e-05 8.17778e-05 5.07640e+01 5.23824e+01 |
| | 1028 | 22.0 6.81200e-05 6.98820e-05 4.33480e+01 4.48150e+01 |
| | 1029 | 23.0 5.77090e-05 5.97167e-05 3.70720e+01 3.83484e+01 |
| | 1030 | 24.0 4.89770e-05 5.10300e-05 3.17510e+01 3.28224e+01 |
| | 1031 | 25.0 4.16330e-05 4.36069e-05 2.72310e+01 2.81003e+01 |
| | 1032 | 27.5 2.79110e-05 2.94361e-05 1.86580e+01 1.90856e+01 |
| | 1033 | 30.0 1.88520e-05 1.98704e-05 1.28900e+01 1.30004e+01 |
| | 1034 | 32.5 1.28100e-05 1.34132e-05 8.98200e+00 8.89271e+00 |
| | 1035 | 35.0 8.75820e-06 9.05433e-06 6.31960e+00 6.11987e+00 |
| | 1036 | 37.5 6.04210e-06 6.11197e-06 4.49210e+00 4.24812e+00 |
| | 1037 | 40.0 4.21740e-06 4.12578e-06 3.22480e+00 2.98462e+00 |
| | 1038 | 42.5 2.98450e-06 2.78504e-06 2.33450e+00 2.13171e+00 |
| | 1039 | 45.0 2.14370e-06 1.88000e-06 1.70000e+00 1.55597e+00 |
| | 1040 | 47.5 1.56200e-06 1.26906e-06 1.24110e+00 1.16733e+00 |
| | 1041 | 50.0 1.14930e-06 1.23132e-06 9.04990e-01 9.04984e-01 |
| | 1042 | 55.0 6.27730e-07 6.24007e-07 4.73870e-01 4.58220e-01 |
| | 1043 | 60.0 3.37430e-07 3.16234e-07 2.39460e-01 2.31809e-01 |
| | 1044 | 65.0 1.74280e-07 1.60261e-07 1.15470e-01 1.17069e-01 |
| | 1045 | 70.0 8.54410e-08 8.12167e-08 5.28480e-02 5.89215e-02 |
| | 1046 | 75.0 3.98420e-08 4.11589e-08 2.28330e-02 2.94534e-02 |
| | 1047 | 80.0 1.58630e-08 2.08585e-08 9.11640e-03 1.45196e-02 |
| | 1048 | 85.0 6.77050e-09 1.05706e-08 3.77200e-03 6.95142e-03 |
| | 1049 | 90.0 2.86410e-09 5.35698e-09 1.49980e-03 3.11604e-03 |
| | 1050 | 95.0 1.15620e-09 2.71480e-09 5.50290e-04 1.17236e-03 |
| | 1051 | 100.0 4.05740e-10 2.02856e-10 1.87330e-04 1.87336e-04 |
| | 1052 | 105.0 1.26610e-10 2.02856e-10 6.71190e-05 8.59079e-05 |
| | 1053 | 109.2 1.26610e-10 2.02856e-10 3.72008e-44 0.00000e+00 |
| | 1054 | |
| | 1055 | |
| | 1056 | CWAVLG 250. 700. simulation range of wavelength band for cherenkov radiation in nm |
| | 1057 | CWAVRANG 250. 700 input range for all three CERQEFF files |
| | 1058 | * CERFIL 1 cherenkov output file, last item of bunch gives wavelength (part of CERWLEN) |
| | 1059 | CERSIZ 1 maximal bunch size of CHphotons treated together |
| | 1060 | CERQEF F T F quantum efficiency of PMT, atmospheric absorption of Cherenkov photons, mirror reflectivity of telescope |
| | 1061 | |
| | 1062 | * ATMOD 22 atmospheric parametrisation, here: US STD atm parametrised by Keilhauer |
| | 1063 | * MAGNET |
| | 1064 | OBSLEV 4100.E2 observation level above sea level in cm |
| | 1065 | |
| | 1066 | ATMOSPHERE 7 T |
| | 1067 | |
| | 1068 | TELESCOPE 0 0 0 2.0e2 0 they must not overlap... |
| | 1069 | |
| | 1070 | CSCAT 1 0 0 number of uses of each Cherenkov event, maximum scattering of core location in cm |
| | 1071 | |
| | 1072 | * ECUTS 0.3 0.3 0.03 0.03 low energy cutoffs, here are std. values in GeV |
| | 1073 | |
| | 1074 | * PAROUT T F |
| | 1075 | DATBAS T |
| | 1076 | |
| | 1077 | DIRECT ./ defines name of output directory |
| | 1078 | DATDIR ./run/ defines name of input directory |
| | 1079 | USER tbretz |
| | 1080 | |
| | 1081 | THIN 0 0.1 0 turn off thinning |
| | 1082 | *** NO MORE DIRECTIVES FOUND *** |
