wiki:InstallingCorsika

Version 49 (modified by tbretz, 4 months ago) (diff)

--

Installation

Extraction

Get the latest version of Corsik from https://www.ikp.kit.edu/corsika/ . Due to missing fetaures and a bug in earlier versions, this must be at least v7.69.

Extract Corsika

[0] tar xvf corsika-76900.tar.bz2

Setup and Compilation

Now, call coconut

[1] cd corsika-76900
[2] ./coconut

You can use the defaults (just press return) for

  • Compile in 32 or 64bit mode? [2: 64 bit]
  • Which high energy hadronic interaction model do you want to use? [4: QGSJET 01C]
  • Which low energy hadronic interaction model do you want to use? [1: GHEISHA 2002d]
  • Which detector geometry do you have? [1: horizontal flat]

Note that so far we were always using FLUKA: WHY??? Shell we use FLUKA as well?

Now you should get a long list of possible options

Which additional CORSIKA program options do you need?

Enter

(multiple selections accepted, leading '-' removes option): 1a 1c 1e 2a 7c 9 

These options correspond to: CERENKOV, IACT, CEFFIC, TRAJECT, VIEWCONE and ATMEXT.

Now you will get asked more question, answer them as follows

  • Cherenkov light vertical (longitudinal) distribution option? [3: No Cherenkov light distribution at all]
  • Do you want Cherenkov light emission angle wavelength dependence? [2: depending on wavelength] {This corresponds to the CERWLEN option}

Now, coconut can create the input files, configure the make system and make. To proceed, select *** Finish selection *** (z).

You should see

  Are you sure you want to continue with these current option selection:
     TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
     yes or no ? (default: yes) >

Which you acknowledge. Now, go ahead!

For completeness, here is the full session with Corsika 7.64

tbretz@d4m:~/SW/corsika-76900> coconut
 

|===========================================================================|
|                          Welcome to COCONUT (v3.1)                        |
|                  -- the CORSIKA CONfiguration UTility --                  |
|===========================================================================|
|                                                                           |
|                                                                           |
|           create an executable of a specific CORSIKA version              |
|                                                                           |
|                                                                           |
|        Please read the documentation for a detailed description           |
|                   of the options and how to use it.                       |
|                                                                           |
|             Try './coconut -h' to get some help about COCONUT             |
|    Use './coconut --expert' to enable additional configuration steps.     |
|                                                                           |
| (press 'Enter' to select an option followed by "[DEFAULT]" or "[CACHED]") |
|                                                                           |
|===========================================================================|




---------------------------------------------------------------------------
 Compile in 32 or 64bit mode ? 
    1 - Force 32bit mode
    2 - Use compiler default ('-m64' on a 64bit machine) [DEFAULT]

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 
    SELECTED         : NOM32 


---------------------------------------------------------------------------
 Which high energy hadronic interaction model do you want to use ?
    1 - DPMJET-III (2017.1) with PHOJET 1.20.0
    2 - EPOS LHC
    3 - NEXUS 3.97
    4 - QGSJET 01C (enlarged commons) [DEFAULT]
    5 - QGSJETII-04
    6 - SIBYLL 2.3c
    7 - VENUS 4.12

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 
    SELECTED         : QGSJET01 


---------------------------------------------------------------------------
 Which low energy hadronic interaction model do you want to use ?
    1 - GHEISHA 2002d (double precision) [DEFAULT]
    2 - FLUKA 
    3 - URQMD 1.3cr

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 
    SELECTED         : GHEISHA 


---------------------------------------------------------------------------
 Which detector geometry do you have ?
    1 - horizontal flat detector array [DEFAULT]
    2 - non-flat (volume) detector geometry
    3 - vertical string detector geometry

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 
    SELECTED         : HORIZONTAL 


---------------------------------------------------------------------------
  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO 

 Which additional CORSIKA program options do you need ?
    1a - Cherenkov version
    1b - Cherenkov version using Bernlohr IACT routines (for telescopes)
    1c - apply atm. absorption, mirror reflectivity & quantum eff.
    1d - Auger Cherenkov longitudinal distribution
    1e - TRAJECTory version to follow motion of source on the sky
    2 - LPM-effect without thinning
    2a - THINning version (includes LPM)
    2b - MULTIple THINning version (includes LPM)
    3 - PRESHOWER version for EeV gammas
    4 - NEUTRINO version
    4a - NUPRIM primary neutrino version with HERWIG
    4b - ICECUBE1 FIFO version
    4c - ICECUBE2 gzip/pipe output
    5 - STACK INput of secondaries, no primary particle
    6 - CHARMed particle/tau lepton version with PYTHIA
    6a - TAU LEPton version with PYTHIA
    7 - SLANT depth instead of vertical depth for longi-distribution
    7a - CURVED atmosphere version
    7b - UPWARD particles version
    7c - VIEWCONE version
    8a - shower PLOT version (PLOTSH) (only for single events)
    8b - shower PLOT(C) version (PLOTSH2) (only for single events)
    8c - ANAlysis HISTos & THIN (instead of particle file)
    8d - Auger-histo file & THIN
    8e - MUON-histo file
    9 - external atmosphere functions (table interpolation)
         (using bernlohr C-routines)
    9a - EFIELD version for electrical field in atmosphere
    9b - RIGIDITY Ooty version rejecting low-energy primaries entering Earth-magnetic field
    10a - DYNamic intermediate particle STACK
    10b - Remote Control for Corsika
    a - CONEX for high energy MC and cascade equations
    b - PARALLEL treatment of subshowers (includes LPM)
    c - CoREAS Radio Simulations
    d1 - Inclined observation plane
    d2 - ROOT particle OUTput file
    e - interaction test version (only for 1st interaction)
    f - Auger-info file instead of dbase file
    g - COMPACT particle output file
    h - MUPROD to write decaying muons
    h2 - prEHISTORY of muons: mother and grandmother
    k - annitest cross-section version (obsolete)
    l - hit Auger detector (steered by AUGSCT)
      - ------------------------------------------
    y - *** Reset selection ***
    z - *** Finish selection *** [DEFAULT]

    r - restart (reset all options to cached values)
    x - exit make

    (multiple selections accepted, leading '-' removes option): 1a 1c 1e 2a 7c 9


---------------------------------------------------------------------------
 Cherenkov light vertical (longitudinal) distribution option ? 
    1 - Photons counted only in the step where emitted [DEFAULT]
    2 - Photons counted in every step down to the observation level
        (compatible with old versions but inefficient)
    3 - No Cherenkov light distribution at all

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 3
    SELECTED         : NOCLONG 


---------------------------------------------------------------------------
 Do you want Cherenkov light emission angle wavelength dependence ? 
    1 - Emission angle is wavelength independent [DEFAULT]
    2 - Emission angle depending on wavelength

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 2
    SELECTED         : CERWLEN 
    SELECTED         : CERENKOV 
    NOT COMPATIBLE TO: COMPACT VOLUMECORR INTTEST ANAHIST AUGERHIST MUONHIST AUGCERLONG ICECUBE1 ICECUBE2 
    SELECTED         : CEFFIC 
    NOT COMPATIBLE TO: INTTEST CURVED AUGCERLONG 
    SELECTED         : TRAJECT 
    NOT COMPATIBLE TO: INTTEST UPWARD 
    SELECTED         : THIN 
    NOT COMPATIBLE TO: ICECUBE1 INTTEST LPM 
    SELECTED         : VIEWCONE 
    NOT COMPATIBLE TO: VOLUMECORR INTTEST 
    SELECTED path to installation: "/home/tbretz/SW/corsika-76900/bernlohr" 
    SELECTED         : ATMEXT 
    NOT COMPATIBLE TO: INTTEST 


---------------------------------------------------------------------------
  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO 
  selection: TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG

 Which additional CORSIKA program options do you need ?
    1a - Cherenkov version
    1b - Cherenkov version using Bernlohr IACT routines (for telescopes)
    1c - apply atm. absorption, mirror reflectivity & quantum eff.
    1d - Auger Cherenkov longitudinal distribution
    1e - TRAJECTory version to follow motion of source on the sky
    2 - LPM-effect without thinning
    2a - THINning version (includes LPM)
    2b - MULTIple THINning version (includes LPM)
    3 - PRESHOWER version for EeV gammas
    4 - NEUTRINO version
    4a - NUPRIM primary neutrino version with HERWIG
    4b - ICECUBE1 FIFO version
    4c - ICECUBE2 gzip/pipe output
    5 - STACK INput of secondaries, no primary particle
    6 - CHARMed particle/tau lepton version with PYTHIA
    6a - TAU LEPton version with PYTHIA
    7 - SLANT depth instead of vertical depth for longi-distribution
    7a - CURVED atmosphere version
    7b - UPWARD particles version
    7c - VIEWCONE version
    8a - shower PLOT version (PLOTSH) (only for single events)
    8b - shower PLOT(C) version (PLOTSH2) (only for single events)
    8c - ANAlysis HISTos & THIN (instead of particle file)
    8d - Auger-histo file & THIN
    8e - MUON-histo file
    9 - external atmosphere functions (table interpolation)
         (using bernlohr C-routines)
    9a - EFIELD version for electrical field in atmosphere
    9b - RIGIDITY Ooty version rejecting low-energy primaries entering Earth-magnetic field
    10a - DYNamic intermediate particle STACK
    10b - Remote Control for Corsika
    a - CONEX for high energy MC and cascade equations
    b - PARALLEL treatment of subshowers (includes LPM)
    c - CoREAS Radio Simulations
    d1 - Inclined observation plane
    d2 - ROOT particle OUTput file
    e - interaction test version (only for 1st interaction)
    f - Auger-info file instead of dbase file
    g - COMPACT particle output file
    h - MUPROD to write decaying muons
    h2 - prEHISTORY of muons: mother and grandmother
    k - annitest cross-section version (obsolete)
    l - hit Auger detector (steered by AUGSCT)
      - ------------------------------------------
    y - *** Reset selection ***
    z - *** Finish selection *** [DEFAULT]

    r - restart (reset all options to cached values)
    x - exit make

    (multiple selections accepted, leading '-' removes option): z
  Are you sure you want to continue with these current option selection:
    TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
    yes or no ? (default: yes) > 


-------------------------------------------------------------------------
  Your final selection to build CORSIKA is:

  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO 
  selection: TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG



---------------------------------------------------------------------------
 Configuration is finished. How do you want to proceed ? 
    f - Compiling and remove temporary files [DEFAULT]
    k - Compile and keep extracted CORSIKA source code
    n - Just extract source code. Do not compile!

    r - restart (reset all options to cached values)
    x - exit make

    (only one choice possible): 
    SELECTED         : COMPILE 
checking whether to enable maintainer-specific portions of Makefiles... no
checking build system type... x86_64-unknown-linux-gnu
checking host system type... x86_64-unknown-linux-gnu
checking for a BSD-compatible install... /usr/bin/install -c
checking whether build environment is sane... yes
checking for a thread-safe mkdir -p... /bin/mkdir -p
checking for gawk... gawk
checking whether make sets $(MAKE)... yes
checking whether make supports nested variables... yes
checking to compile without optimisation and system flags... (cached) no
checking whether to generate debug... (cached) yes
checking for pgf77... no
checking for ifc... no
checking for ifort... no
checking for gfortran... gfortran
checking whether the Fortran 77 compiler works... yes
checking for Fortran 77 compiler default output file name... a.out
checking for suffix of executables... 
checking whether we are cross compiling... no
checking for suffix of object files... o
checking whether we are using the GNU Fortran 77 compiler... yes
checking whether gfortran accepts -g... yes
checking for cc... cc
checking whether we are using the GNU C compiler... yes
checking whether cc accepts -g... yes
checking for cc option to accept ISO C89... none needed
checking whether cc understands -c and -o together... yes
checking for style of include used by make... GNU
checking dependency style of cc... gcc3
checking for g++... g++
checking whether we are using the GNU C++ compiler... yes
checking whether g++ accepts -g... yes
checking dependency style of g++... gcc3
checking for cpp... cpp
checking how to run the C preprocessor... cpp
checking how to get verbose linking output from gfortran... -v
checking for Fortran 77 libraries of gfortran...  -L/usr/lib/gcc/x86_64-linux-gnu/7 -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../../lib -L/lib/x86_64-linux-gnu -L/lib/../lib -L/usr/lib/x86_64-linux-gnu -L/usr/lib/../lib -L/usr/lib/gcc/x86_64-linux-gnu/7/../../.. -lgfortran -lm -lquadmath
checking how to print strings... printf
checking for a sed that does not truncate output... /bin/sed
checking for grep that handles long lines and -e... /bin/grep
checking for egrep... /bin/grep -E
checking for fgrep... /bin/grep -F
checking for ld used by cc... /usr/bin/ld
checking if the linker (/usr/bin/ld) is GNU ld... yes
checking for BSD- or MS-compatible name lister (nm)... /usr/bin/nm -B
checking the name lister (/usr/bin/nm -B) interface... BSD nm
checking whether ln -s works... yes
checking the maximum length of command line arguments... 1572864
checking whether the shell understands some XSI constructs... yes
checking whether the shell understands "+="... yes
checking how to convert x86_64-unknown-linux-gnu file names to x86_64-unknown-linux-gnu format... func_convert_file_noop
checking how to convert x86_64-unknown-linux-gnu file names to toolchain format... func_convert_file_noop
checking for /usr/bin/ld option to reload object files... -r
checking for objdump... objdump
checking how to recognize dependent libraries... pass_all
checking for dlltool... no
checking how to associate runtime and link libraries... printf %s\n
checking for ar... ar
checking for archiver @FILE support... @
checking for strip... strip
checking for ranlib... ranlib
checking command to parse /usr/bin/nm -B output from cc object... ok
checking for sysroot... no
checking for mt... mt
checking if mt is a manifest tool... no
checking for ANSI C header files... yes
checking for sys/types.h... yes
checking for sys/stat.h... yes
checking for stdlib.h... yes
checking for string.h... yes
checking for memory.h... yes
checking for strings.h... yes
checking for inttypes.h... yes
checking for stdint.h... yes
checking for unistd.h... yes
checking for dlfcn.h... yes
checking for objdir... .libs
checking if cc supports -fno-rtti -fno-exceptions... no
checking for cc option to produce PIC... -fPIC -DPIC
checking if cc PIC flag -fPIC -DPIC works... yes
checking if cc static flag -static works... yes
checking if cc supports -c -o file.o... yes
checking if cc supports -c -o file.o... (cached) yes
checking whether the cc linker (/usr/bin/ld -m elf_x86_64) supports shared libraries... yes
checking whether -lc should be explicitly linked in... no
checking dynamic linker characteristics... GNU/Linux ld.so
checking how to hardcode library paths into programs... immediate
checking whether stripping libraries is possible... yes
checking if libtool supports shared libraries... yes
checking whether to build shared libraries... yes
checking whether to build static libraries... yes
checking how to run the C++ preprocessor... g++ -E
checking for ld used by g++... /usr/bin/ld -m elf_x86_64
checking if the linker (/usr/bin/ld -m elf_x86_64) is GNU ld... yes
checking whether the g++ linker (/usr/bin/ld -m elf_x86_64) supports shared libraries... yes
checking for g++ option to produce PIC... -fPIC -DPIC
checking if g++ PIC flag -fPIC -DPIC works... yes
checking if g++ static flag -static works... yes
checking if g++ supports -c -o file.o... yes
checking if g++ supports -c -o file.o... (cached) yes
checking whether the g++ linker (/usr/bin/ld -m elf_x86_64) supports shared libraries... yes
checking dynamic linker characteristics... (cached) GNU/Linux ld.so
checking how to hardcode library paths into programs... immediate
checking if libtool supports shared libraries... yes
checking whether to build shared libraries... yes
checking whether to build static libraries... yes
checking for gfortran option to produce PIC... -fPIC
checking if gfortran PIC flag -fPIC works... yes
checking if gfortran static flag -static works... yes
checking if gfortran supports -c -o file.o... yes
checking if gfortran supports -c -o file.o... (cached) yes
checking whether the gfortran linker (/usr/bin/ld -m elf_x86_64) supports shared libraries... yes
checking dynamic linker characteristics... (cached) GNU/Linux ld.so
checking how to hardcode library paths into programs... immediate
./configure: line 19111: BOOST_PROGRAM_OPTIONS: command not found
checking to make dist... (cached) no
checking 32 bits compatibility... no
checking for stdint.h... (cached) yes
checking explicitly use the -m32 flag... (cached) no
checking explicitly use the -frecord-marker=4 flag... (cached) no
checking detector geometry... horizontal
checking time library... rc
checking high energy model... qgsjet01
checking for fluka... no
checking low energy model... gheisha
checking for thinning... (cached) yes
checking for Cherenkov grid detector... (cached) yes
checking for bernlohr... /home/tbretz/SW/corsika-76900/bernlohr
checking for Cherenkov telescope detector... (cached) no
checking new Cherenkov counting... (cached) no
checking no cherenkov longitudinla profile... (cached) yes
checking Cherenkov wavelength dependent... (cached) yes
checking simulate telescope detector... (cached) yes
checking for external IACT Cherenkov file... (cached) no
checking external atmosphere... (cached) yes
checking for multiple thinning... (cached) no
checking neutrino primaries... (cached) no
checking ICECUBE FiFo... (cached) no
checking ICECUBE gzip/pipe output... (cached) no
checking BAACK modules with C++0x option... (cached) no
checking BAACK modules with C++11 option... (cached) no
checking BAACK modules with C++14 option... (cached) no
checking BAACK modules with C++1z option... (cached) no
checking for BAACK... no
checking DYNamic particle STACK... (cached) no
checking emulation of normal stack... (cached) no
checking stores every saved particle to file... (cached) no
checking reimplementation of the ICECUBE1 module... (cached) no
checking user defined intermediate particle stack... (cached) no
checking Remote acces to corsika... (cached) no
checking plot shower... (cached) no
checking plot shower C... (cached) no
checking for cernlib... no
checking interaction test... (cached) no
checking slant profile... (cached) no
checking curved atmosphere... (cached) no
checking upward particles... (cached) no
checking viewcone... (cached) yes
checking cerenkov longitudinal profile... (cached) no
checking photon preshower... (cached) no
checking MUPROD... (cached) no
checking compact output... (cached) no
checking annitest... (cached) no
checking LPM effect... (cached) no
checking stack input... (cached) no
checking for pythia... no
checking pythia charm production... (cached) no
checking for herwig... no
checking herwig neutrino production (NUPRIM)... (cached) no
checking augerhit... (cached) no
checking pythia tau lepton production... (cached) no
checking particle history... (cached) no
checking trajectory... (cached) yes
checking electrical field... (cached) no
checking rigidity... (cached) no
checking for root... no
checking particle output in root file... (cached) no
checking machine independent output... (cached) no
checking for COASTUSERLIB... no
checking External COAST user library... (cached) no
checking CoREAS radio simulations... (cached) no
checking Inclined observation level... (cached) no
checking for conex... no
checking CONEX cascade equation (CONEX)... (cached) no
checking produce analysis histograms... (cached) no
checking augerinfo... (cached) no
checking augerhist... (cached) no
checking muonhist... (cached) no
checking parallel computation... (cached) no
checking for mpirunner_lib... no
checking parallel computation with MPI... (cached) no
checking do not compile binaries, just extract CORSIKA compilefile... (cached) no
checking to keep the CORSIKA compilefile... (cached) no
checking that generated files are newer than configure... done
configure: creating ./config.status
config.status: creating Makefile
config.status: creating baack/Makefile
config.status: creating bernlohr/Makefile
config.status: creating conex/Makefile
config.status: creating dpmjet/Makefile
config.status: creating epos/Makefile
config.status: creating pythia/Makefile
config.status: creating herwig/Makefile
config.status: creating nexus/Makefile
config.status: creating urqmd/Makefile
config.status: creating src/Makefile
config.status: creating run/Makefile
config.status: creating doc/Makefile
config.status: creating lib/Makefile
config.status: creating coast/Makefile
config.status: creating coast/Documentation/Makefile
config.status: creating coast/CorsikaOptions/rootout/Makefile
config.status: creating coast/CorsikaOptions/CoREAS/Makefile
config.status: creating coast/CorsikaOptions/Makefile
config.status: creating coast/CorsikaOptions/InclinedPlane/Makefile
config.status: creating coast/CorsikaFileIO/Makefile
config.status: creating coast/CorsikaInterface/Makefile
config.status: creating coast/CorsikaToROOT/Makefile
config.status: creating coast/CorsikaROOT/Makefile
config.status: creating coast/CorsikaIntern/Makefile
config.status: creating include/config.h
config.status: executing depfiles commands
config.status: executing libtool commands
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -


Compile CORSIKA in "/home/tbretz/SW/corsika-76900/lib/unknown" and copy executable in "/home/tbretz/SW/corsika-76900/run" 

Making install in bernlohr
make[1]: Entering directory '/home/tbretz/SW/corsika-76900/bernlohr'
cc -DHAVE_CONFIG_H -I. -I../include     -g -D_FILE_OFFSET_BITS=64 -MT libbern_a-fileopen.o -MD -MP -MF .deps/libbern_a-fileopen.Tpo -c -o libbern_a-fileopen.o `test -f 'fileopen.c' || echo './'`fileopen.c
mv -f .deps/libbern_a-fileopen.Tpo .deps/libbern_a-fileopen.Po
cc -DHAVE_CONFIG_H -I. -I../include     -g -D_FILE_OFFSET_BITS=64 -MT libbern_a-straux.o -MD -MP -MF .deps/libbern_a-straux.Tpo -c -o libbern_a-straux.o `test -f 'straux.c' || echo './'`straux.c
mv -f .deps/libbern_a-straux.Tpo .deps/libbern_a-straux.Po
rm -f libbern.a
ar cru libbern.a libbern_a-fileopen.o libbern_a-straux.o 
ar: `u' modifier ignored since `D' is the default (see `U')
ranlib libbern.a
cc -DHAVE_CONFIG_H -I. -I../include     -g -D_FILE_OFFSET_BITS=64 -MT libatmext_a-atmo.o -MD -MP -MF .deps/libatmext_a-atmo.Tpo -c -o libatmext_a-atmo.o `test -f 'atmo.c' || echo './'`atmo.c
mv -f .deps/libatmext_a-atmo.Tpo .deps/libatmext_a-atmo.Po
cc -DHAVE_CONFIG_H -I. -I../include     -g -D_FILE_OFFSET_BITS=64 -MT libatmext_a-fileopen.o -MD -MP -MF .deps/libatmext_a-fileopen.Tpo -c -o libatmext_a-fileopen.o `test -f 'fileopen.c' || echo './'`fileopen.c
mv -f .deps/libatmext_a-fileopen.Tpo .deps/libatmext_a-fileopen.Po
cc -DHAVE_CONFIG_H -I. -I../include     -g -D_FILE_OFFSET_BITS=64 -MT libatmext_a-straux.o -MD -MP -MF .deps/libatmext_a-straux.Tpo -c -o libatmext_a-straux.o `test -f 'straux.c' || echo './'`straux.c
mv -f .deps/libatmext_a-straux.Tpo .deps/libatmext_a-straux.Po
rm -f libatmext.a
ar cru libatmext.a libatmext_a-atmo.o libatmext_a-fileopen.o libatmext_a-straux.o 
ar: `u' modifier ignored since `D' is the default (see `U')
ranlib libatmext.a
make[2]: Entering directory '/home/tbretz/SW/corsika-76900/bernlohr'
 /bin/mkdir -p '/home/tbretz/SW/corsika-76900/lib/unknown'
 /usr/bin/install -c -m 644  libbern.a libatmext.a '/home/tbretz/SW/corsika-76900/lib/unknown'
 ( cd '/home/tbretz/SW/corsika-76900/lib/unknown' && ranlib libbern.a )
 ( cd '/home/tbretz/SW/corsika-76900/lib/unknown' && ranlib libatmext.a )
make[2]: Nothing to be done for 'install-data-am'.
make[2]: Leaving directory '/home/tbretz/SW/corsika-76900/bernlohr'
make[1]: Leaving directory '/home/tbretz/SW/corsika-76900/bernlohr'
Making install in src
make[1]: Entering directory '/home/tbretz/SW/corsika-76900/src'
cpp  -DHAVE_CONFIG_H -I../include -traditional-cpp corsika.F > corsikacompilefile.f
gfortran -O0  -g -std=legacy -c -o corsika-corsikacompilefile.o `test -f 'corsikacompilefile.f' || echo './'`corsikacompilefile.f
depbase=`echo tobuf.o | sed 's|[^/]*$|.deps/&|;s|\.o$||'`;\
cc -DHAVE_CONFIG_H -I. -I../include        -g -D_FILE_OFFSET_BITS=64 -MT tobuf.o -MD -MP -MF $depbase.Tpo -c -o tobuf.o tobuf.c &&\
mv -f $depbase.Tpo $depbase.Po
gfortran -O0  -g -std=legacy -c -o corsika-qgsjet01d.o `test -f 'qgsjet01d.f' || echo './'`qgsjet01d.f
gfortran -O0  -g -std=legacy -c -o corsika-gheisha_2002d.o `test -f 'gheisha_2002d.f' || echo './'`gheisha_2002d.f
depbase=`echo timerc.o | sed 's|[^/]*$|.deps/&|;s|\.o$||'`;\
cc -DHAVE_CONFIG_H -I. -I../include        -g -D_FILE_OFFSET_BITS=64 -MT timerc.o -MD -MP -MF $depbase.Tpo -c -o timerc.o timerc.c &&\
mv -f $depbase.Tpo $depbase.Po
/bin/bash ../libtool  --tag=F77   --mode=link gfortran -O0  -g -std=legacy  -D_FILE_OFFSET_BITS=64 -o corsika corsika-corsikacompilefile.o tobuf.o  corsika-qgsjet01d.o    corsika-gheisha_2002d.o   timerc.o -L/home/tbretz/SW/corsika-76900/lib/unknown -lbern -latmext      
libtool: link: gfortran -O0 -g -std=legacy -D_FILE_OFFSET_BITS=64 -o corsika corsika-corsikacompilefile.o tobuf.o corsika-qgsjet01d.o corsika-gheisha_2002d.o timerc.o  -L/home/tbretz/SW/corsika-76900/lib/unknown -lbern -latmext
make[2]: Entering directory '/home/tbretz/SW/corsika-76900/src'
 /bin/mkdir -p '/home/tbretz/SW/corsika-76900/run'
  /bin/bash ../libtool   --mode=install /usr/bin/install -c corsika '/home/tbretz/SW/corsika-76900/run'
libtool: install: /usr/bin/install -c corsika /home/tbretz/SW/corsika-76900/run/corsika
make[2]: Nothing to be done for 'install-data-am'.
make[2]: Leaving directory '/home/tbretz/SW/corsika-76900/src'
make[1]: Leaving directory '/home/tbretz/SW/corsika-76900/src'
Making install in .
make[1]: Entering directory '/home/tbretz/SW/corsika-76900'
make[2]: Entering directory '/home/tbretz/SW/corsika-76900'
make  install-exec-hook
make[3]: Entering directory '/home/tbretz/SW/corsika-76900'
 
 
 
--> "corsika76900Linux_QGSJET_gheisha_thin" successfully installed in :
    /home/tbretz/SW/corsika-76900/run/
 
--> You can run CORSIKA in /home/tbretz/SW/corsika-76900/run/ using for instance :
    ./corsika76900Linux_QGSJET_gheisha_thin < all-inputs > output.txt
 
 
 
make[3]: Leaving directory '/home/tbretz/SW/corsika-76900'
make[2]: Nothing to be done for 'install-data-am'.
make[2]: Leaving directory '/home/tbretz/SW/corsika-76900'
make[1]: Leaving directory '/home/tbretz/SW/corsika-76900'

Options description

CERENKOV[1a]: Enables production of Cherenkov photons

IACT[1b]: Redirects photons to eventio (telescope.dat) file (CER files exists, but only for one telescope and it is empty), extends wavelength range to 2000nm.

CERWLEN: Calculates a wavelength dependent refractive index for each step and replaces the production height (8-th) in the photon output by the wavelength.

IACTEXT: Write also particles to eventio file

CEFFIC[1c]: All three look-up-tables (QE, Atmabs, Mirror) are limited to 105 values between 180nm and 700nm -> This turns the extension to 2000nm in the IACT option off

TRAJECT[1e]: To be able to simulate along a trajetory (zd, az, orientation to magnetic field) with "TRAFLG T" (Turn off with "TRAFLG F")

THIN[2a]: Thinning is mainly interesting for >PeV showers and calculates average results for high densities. With CEFFIC and/or CERWLEN, the wavelength is written 8th value to the Cherenkov bunch into the CER-file instead of the weight. To turn the thinning off for the production use "THIN 0 1 0" (The production should behave like with no THIN option). To turn thinning only off for electromagnetic-particles use "THINEM 0 1"

VIEWCONE[7c]: Logically allows to simulate a cone rather than a rectangle in Alt/Az

ATMEXT[9]: To use an atmosphere from an input file rather than a built-in atmosphere

This is from the tracking of the electrons through the Earth's magnetic field

C  LIMITING FACTOR FOR STEP SIZE OF ELECTRON IN MAGNETIC FIELD
#if __CERENKOV__ && __IACT__
C  LIMIT IN DEFLECTION ANGLE IS 2.5 MILLIRADIAN = 0.143 DEG
C  WE USE A LIMIT OF ABOUT 0.05 DEG (APPROX. 1 MILLIRAD)
      BLIMIT   = 0.001D0 / BNORM
#else
C  WE USE A LIMIT OF ABOUT 11.4 DEG (0.2 RAD)
      BLIMIT   = 0.2D0 / BNORM
#endif

This is from the muon tracking (MUTRAC)

#if __CERENKOV__ && __IACT__
C  SCATTERING ANGLES OF MUONS SHOULD BE SMALLER THAN THE PIXEL SIZE.
*     AUX = MIN( 1.D0, 0.015D0*GAMMA )
C  THE SAME SHOULD HOLD FOR DEFLECTION IN THE GEOMAGNETIC FIELD.
C  HERE USING A MAXIMUM RMS SCATTERING / DEFLECTION ANGLE OF 0.05 DEG
C  AND APPROXIMATE ALL BETA*GAMMA TERMS BY GAMMA.
Cxx      Write(*,*) 'mu step old-style step=',MIN( 1.D0,0.015D0*GAMMA )
C  FOR A MEAN SCATTERING ANGLE THETA WE HAVE A STEP LENGTH OF ABOUT
C  (THETA / (13.6 MEV/(BETA*C*P))**2 RADIATION LENGTHS (PDG),
C  NOT TAKING INTO ACCOUNT THE NON-GAUSSIAN PART OF THE DISTIBUTION.
C  FOR THE MOMENT DON''T CARE ABOUT THE DIFFERENCE BETWEEN THE
C  'COULOMB SCATTERING LENGTH' 37.7 G/CM**2 (=C(21)) AND THE
C  RADIATION LENGTH OF 36.66 OR 36.62 G/CM**2 IN AIR.
C  NOTE: PI/180/(13.6 MEV/(BETA*C*P)) APPROX 0.136*GAMMA FOR MUONS.
      AUX = MIN( 1.D0, ((0.05*0.136)*GAMMA)**2 )
      IF ( BNORMC .GT. 0.D0 ) THEN
C  NOTE: PI/180*PAMA(5)*BETA*GAMMA APPROX 0.00185*GAMMA
         AUX = MIN( AUX, (0.05*0.00185)*GAMMA*RHOF(H)/(BNORMC*C(21)) )
      ENDIF
Cxx      Write(*,*) 'mu step new-style step=',AUX
#else
      AUX = MIN( 10.D0, 0.015D0*GAMMA )
#endif

For me this seems only relevant if someone wants to backtrack a single muon, i.e. determin its arrival direction.

An interesting possibility of the IACT option is:

j) Starting with version 1.25, the package has been prepared for importance sampling ofcore position offsets. This would mean that actual core offsets can be generated in anon-uniform distribution and can extend to different distances, depending on primaryenergy, primary type, zenith angle and so on. This package, however, does not provide areal implementation of importance sampling (other than for testing that the later stagesof the processing properly get the weights for each event). If you do nothing about it,you will get uniformly distributed core offsets as before. If you plan to make use ofimportance sampling, you have to replace the file ’sampling.c’ with an implementationof your choice.

A problem with the IACT option might be this:

m) Starting with version 1.38, the dynamic range in a telescope is basically unlimited due to automatic thinning of bunches. When a detector sphere is hit by more than the givenmaximum number of bunches, the actual number of bunches is reduced by increasingpowers of two, by discarding every second bunch and increasing the bunch size of theremaining bunches by factors of two. Very large eventio buffers are now possible on64-bit machines. They were formerly limited to less than 1 GiB per telescope array andevent and now to 2 GiB on 32-bit machines but can be increased to 4 Terabytes on 64-bitmachines.

MMCS

The current version of MARS is NOT compatible anymore with MMCS output. Although the code is still there and can be enables, it is turned off by default (see MPhotonData.cc)

Conclusion

If we want to do any absorption in Corsika (CEFFIC), we have to patch Corsika to allow wavelength larger than 700nm or accept that we will miss photons above 700nm.

To be able to define telescope positions ("TELESCOPE"), the IACT option is not required. Plain Corsika can nowadays do that.

Enabling IACT might make sense, due to the code snippets shown above, the automatic thinning (faster) and the work on improving runtime by faster interpolation algorithms.

The disadvantage of the IACT option is that the data of all telescopes end up in the same file and that *eventio* format has to be used. Although, the iACT options writes a temporary file for each individual telescope, they are sorted together into one file later.

Using IACT means that we have to understand the automatic thinning and implement it and to make sure eventio works.

I have the impression that eventio files are larger than the corsika output (which is a bit strange)... this has to be checked with long runs and high statistics!

Examples

Here is an example input card for gammas (not all options have yet reasonable values!) Make sure to adapt the DATDIR properly!

RUNNR   100                      run number                                                                                        
EVTNR   0                        number of first shower event                                                                      
NSHOW   1                        number of showers to be generated in one run                                                      

SEED    62748   0       0        seed for hadronic part                                                                            
SEED    62749   0       0        seed for EGS4 part                                                                                
SEED    62750   0       0        seed for Cherenkov part                                                                           
SEED    62751   0       0        seed for 4. random number sequence                                                                
SEED    62752   0       0        seed for 5. random number sequence                                                                

PRMPAR  1                        primary particle code: 1 = photon, 3 = e-, 14 = p+
FIXCHI  0                        starting altitude of primary particle?
TSTART T 			 taken from IACT option (start emission at border of atmosphere)

ERANGE  10000    10000           energy range of primary particle in GeV
ESLOPE  -2                       slope of energy spectrum

PHIP    0.      0.               range of azimuth angle in deg
THETAP	0.	0.               range of zenith angle in deg
VIEWCONE 0.0    0.0              range of viewcone in deg

TRAFLG F                         TRAJECTORY option
* SRCPOS ra dec                  source position
* TRATM  y m d h m s duration    observation time
* TLAT   d m s N                 observatory latitude
* TLONG  d m s N                 observatory longitude

CWAVLG   250.  700.              simulation range of wavelength band for cherenkov radiation in nm
CWAVRANG 250.  700               input range for all three CERQEFF files
* CERFIL  1                      cherenkov output file, last item of bunch gives wavelength (part of CERWLEN)
CERSIZ  1                        maximal bunch size of CHphotons treated together
CERQEF F T F                     quantum efficiency of PMT, atmospheric absorption of Cherenkov photons, mirror reflectivity of telescope

* ATMOD   22                     atmospheric parametrisation, here: US STD atm parametrised by Keilhauer
* MAGNET
OBSLEV  4100.E2                  observation level above sea level in cm

ATMOSPHERE 7 T

TELESCOPE  0 0 0 2.0e2 0         they must not overlap...

CSCAT 1 0 0                      number of uses of each Cherenkov event, maximum scattering of core location in cm

* ECUTS 0.3 0.3 0.03 0.03        low energy cutoffs, here are std. values in GeV

* PAROUT T F
DATBAS T                         write header data in human redable form (ideal to fill that into our database)

DIRECT ./                        defines name of output directory
DATDIR corsika-76400/run         defines name of input directory
USER   tbretz

THIN 0 0.1 0                     turn off thinning

Here is the coresponding output of CORSIKA:

 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)
             D. HECK, J. KNAPP, REPORT FZKA 6097 (1998)
             SEE ALSO WEB PAGE https://www.ikp.kit.edu/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@partner.kit.edu
                         FAX:   (49) 721-608-24075 
                       PHONE:   (49) 721-608-23777 
   AND SEND YOUR LIST-FILE BY E-MAIL

 NUMBER OF VERSION :  7.6900
 DATE   OF VERSION : DECEMBER  20, 2018

 VERSION GENERATED FOR UNIX OR COMPATIBLE SYSTEMS
 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  (RECL IS DEFINED IN BYTES)
  WITH TIMERC DATE ROUTINE
 ZENITH ANGLE DEPENDENCE FOR FLAT DETECTOR ARRAY
  
 INTERFACE FOR EXTERNAL ATMOSPHERIC PROFILES ENABLED
 ===================================================

 CHERENKOV RADIATION IS GENERATED
 ================================

 WITH WAVELENGTH DEPENDENCE
 ATMOSPHERIC ABSORPTION, MIRROR REFLECTIVITY AND QUANTUM EFFICIENCY MAY BE RESPECTED

 TRAJECTORY VERSION TO FOLLOW EMITTING SOURCE
 ============================================
  
 THINNING IS ACTIVE
 ==================

 PRIMARY DIRECTION IS SELECTED FROM VIEWING CONE
 ===============================================

  
 ========== USERS RUN DIRECTIVES FOR THIS SIMULATION ===========================

 DATA CARDS FOR RUN STEERING ARE EXPECTED FROM STANDARD INPUT
   
RUNNR   103                      run number
EVTNR   0                        number of first shower event
NSHOW   1                        number of showers to be generated in one run

SEED    62748   0       0        seed for hadronic part
SEED    62749   0       0        seed for EGS4 part
SEED    62750   0       0        seed for Cherenkov part
SEED    62751   0       0        seed for 4. random number sequence
SEED    62752   0       0        seed for 5. random number sequence

PRMPAR  1                        primary particle code: 1 = photon, 3 = e-, 14 = p+
FIXCHI  0                        starting altitude of primary particle?
TSTART T     taken from IACT option (start emission at border of atmosphere)

ERANGE  10000    10000           energy range of primary particle in GeV
ESLOPE  -2                       slope of energy spectrum

PHIP    0.      0.               range of azimuth angle in deg
THETAP 0. 0.              range of zenith angle in deg
VIEWCONE 0.0    0.0              range of viewcone in deg

TRAFLG F                         TRAJECTORY option
* SRCPOS ra dec                  source position
* TRATM  y m d h m s duration    observation time
* TLAT   d m s N                 observatory latitude
* TLONG  d m s N                 observatory longitude

 Atmospheric profile 7 with 50 levels read from file atmprof7.dat


 Results of the atmosphere fit:
 Layer 1:   0.00 km < h <   3.00 km: a =     -144.838, b =      1192.34, c =       994186
 Layer 2:   3.00 km < h <  10.00 km: a =     -124.071, b =      1173.98, c =       967530
 Layer 3:  10.00 km < h <  50.00 km: a =     0.360027, b =      1412.08, c =       636143
 Layer 4:  50.00 km < h < 100.00 km: a = -0.000824761, b =      810.682, c =       735640
 Layer 5: 100.00 km < h < 109.23 km: a =   0.00221589, b =            1, c =  4.92961e+09


 Altitude [km]    rho(table)     rho(fit)       thick(table)  thick(fit)
        0.0      1.19820e-03   1.19931e-03      1.04750e+03   1.04750e+03
        1.0      1.08430e-03   1.08455e-03      9.33510e+02   9.33403e+02
        2.0      9.81420e-04   9.80765e-04      8.30320e+02   8.30225e+02
        3.0      8.88120e-04   8.89886e-04      7.36920e+02   7.36919e+02
        4.0      8.03260e-04   8.02504e-04      6.52420e+02   6.52375e+02
        5.0      7.25830e-04   7.23703e-04      5.76020e+02   5.76133e+02
        6.0      6.54990e-04   6.52639e-04      5.07040e+02   5.07377e+02
        7.0      5.90060e-04   5.88554e-04      4.44830e+02   4.45372e+02
        8.0      5.30430e-04   5.30762e-04      3.88840e+02   3.89456e+02
        9.0      4.75610e-04   4.78644e-04      3.38570e+02   3.39031e+02
       10.0      4.25210e-04   4.60898e-04      2.93560e+02   2.93557e+02
       11.0      3.78710e-04   3.93854e-04      2.53390e+02   2.50907e+02
       12.0      3.35320e-04   3.36562e-04      2.17710e+02   2.14461e+02
       13.0      2.94600e-04   2.87604e-04      1.86240e+02   1.83317e+02
       14.0      2.56450e-04   2.45768e-04      1.58720e+02   1.56703e+02
       15.0      2.21010e-04   2.10017e-04      1.34870e+02   1.33961e+02
       16.0      1.88580e-04   1.79467e-04      1.14420e+02   1.14527e+02
       17.0      1.59740e-04   1.53361e-04      9.70320e+01   9.79196e+01
       18.0      1.34720e-04   1.31052e-04      8.23400e+01   8.37281e+01
       19.0      1.13400e-04   1.11989e-04      6.99620e+01   7.16010e+01
       20.0      9.54910e-05   9.56985e-05      5.95440e+01   6.12380e+01
       21.0      8.05660e-05   8.17778e-05      5.07640e+01   5.23824e+01
       22.0      6.81200e-05   6.98820e-05      4.33480e+01   4.48150e+01
       23.0      5.77090e-05   5.97167e-05      3.70720e+01   3.83484e+01
       24.0      4.89770e-05   5.10300e-05      3.17510e+01   3.28224e+01
       25.0      4.16330e-05   4.36069e-05      2.72310e+01   2.81003e+01
       27.5      2.79110e-05   2.94361e-05      1.86580e+01   1.90856e+01
       30.0      1.88520e-05   1.98704e-05      1.28900e+01   1.30004e+01
       32.5      1.28100e-05   1.34132e-05      8.98200e+00   8.89271e+00
       35.0      8.75820e-06   9.05433e-06      6.31960e+00   6.11987e+00
       37.5      6.04210e-06   6.11197e-06      4.49210e+00   4.24812e+00
       40.0      4.21740e-06   4.12578e-06      3.22480e+00   2.98462e+00
       42.5      2.98450e-06   2.78504e-06      2.33450e+00   2.13171e+00
       45.0      2.14370e-06   1.88000e-06      1.70000e+00   1.55597e+00
       47.5      1.56200e-06   1.26906e-06      1.24110e+00   1.16733e+00
       50.0      1.14930e-06   1.23132e-06      9.04990e-01   9.04984e-01
       55.0      6.27730e-07   6.24007e-07      4.73870e-01   4.58220e-01
       60.0      3.37430e-07   3.16234e-07      2.39460e-01   2.31809e-01
       65.0      1.74280e-07   1.60261e-07      1.15470e-01   1.17069e-01
       70.0      8.54410e-08   8.12167e-08      5.28480e-02   5.89215e-02
       75.0      3.98420e-08   4.11589e-08      2.28330e-02   2.94534e-02
       80.0      1.58630e-08   2.08585e-08      9.11640e-03   1.45196e-02
       85.0      6.77050e-09   1.05706e-08      3.77200e-03   6.95142e-03
       90.0      2.86410e-09   5.35698e-09      1.49980e-03   3.11604e-03
       95.0      1.15620e-09   2.71480e-09      5.50290e-04   1.17236e-03
      100.0      4.05740e-10   2.02856e-10      1.87330e-04   1.87336e-04
      105.0      1.26610e-10   2.02856e-10      6.71190e-05   8.59079e-05
      109.2      1.26610e-10   2.02856e-10      3.72008e-44   0.00000e+00


CWAVLG   250.  700.              simulation range of wavelength band for cherenkov radiation in nm
CWAVRANG 250.  700               input range for all three CERQEFF files
* CERFIL  1                      cherenkov output file, last item of bunch gives wavelength (part of CERWLEN)
CERSIZ  1                      maximal bunch size of CHphotons treated together
CERQEF F T F                     quantum efficiency of PMT, atmospheric absorption of Cherenkov photons, mirror reflectivity of telescope

* ATMOD   22                     atmospheric parametrisation, here: US STD atm parametrised by Keilhauer
* MAGNET
OBSLEV  4100.E2                  observation level above sea level in cm

ATMOSPHERE 7 T

TELESCOPE  0 0 0 2.0e2 0         they must not overlap...

CSCAT 1 0 0                      number of uses of each Cherenkov event, maximum scattering of core location in cm

* ECUTS 0.3 0.3 0.03 0.03        low energy cutoffs, here are std. values in GeV

* PAROUT T F
DATBAS T

DIRECT ./                        defines name of output directory
DATDIR ./run/                    defines name of input directory
USER   tbretz

THIN 0 0.1 0                     turn off thinning
 *** NO MORE DIRECTIVES FOUND ***

 ATMOSPHERIC REFRACTION IS TAKEN INTO ACCOUNT
 ============================================



 FITTED ATMOSPHERIC PARAMETERS:
 HLAY =   0.0000000000000000        300000.00000000000        1000000.0000000000        5000000.0000000000        10000000.000000000     
 AATM =  -144.83813187979536       -124.07131244989489       0.36002723273742276       -8.2476119271945869E-004   2.2158928172410280E-003
 BATM =   1192.3381318797951        1173.9762361493856        1412.0788076142780        810.68229015979523        1.0000000000000000     
 CATM =   994186.38000000000        967529.82871759660        636143.04000000004        735639.91469823441        4929612368.8889332     

 ========== START OF RUN =======================================================

 PRESENT TIME : 05.06.2019  10:41:12 UTC
 RANDOM NUMBER GENERATOR AT BEGIN OF RUN :
 SEQUENCE =  1  SEED =     62748  CALLS =         0  BILLIONS =         0
 SEQUENCE =  2  SEED =     62749  CALLS =         0  BILLIONS =         0
 SEQUENCE =  3  SEED =     62750  CALLS =         0  BILLIONS =         0
 SEQUENCE =  4  SEED =     62751  CALLS =         0  BILLIONS =         0
 SEQUENCE =  5  SEED =     62752  CALLS =         0  BILLIONS =         0

 PARTICLE OUTPUT TO FILE: ./DAT000103                                                                                                                                                                                                                                                    

 CHERENKOV OUTPUT TO FILE: ./CER000103                                                                                                                                                                                                                                                    

 DBASE OUTPUT TO FILE: ./DAT000103.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:  Jan. 30, 2013  by T. Pierog  |
 |               (version qgsjet01d.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.58059654E-02  5.73089633E-01  5.28304217E-02  2.50000000E+00  2.07000000E+00
 8.20000000E+00  1.00000000E-01  0.00000000E+00  0.00000000E+00  1.00002339E+00
 9.67266333E-03  1.00000008E+00  5.75170728E-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.61909598E-01
 8.99834275E-01  0.00000000E+00  1.03899198E+00  2.71383371E-01  1.37035999E+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.26406662E-01  1.12407142E+00  1.49600000E+02
 1.49600000E+02  2.35531824E-01  2.06000000E-01  1.35000000E-01  2.22000000E-01
 5.00000000E-01  0.00000000E+00  6.35600000E-01  6.92762000E-01  8.74200000E-01
 6.78300000E-01  4.06900000E-01  2.44320607E+00  9.12400642E-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.42300000E-01  8.98200000E-01  9.48900000E-01  9.82400000E-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.98000000E-02
 0.00000000E+00  1.28800000E-02  2.96000000E-02  1.96000000E-02  1.19400000E-02
 2.82000000E-02  0.00000000E+00  1.31000000E-02  2.82000000E-02  1.38000000E-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.97200000E-01  7.26500000E-01  9.57500000E-01  9.00000000E-04  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.40900000E-01
 5.16300000E-01  0.00000000E+00  0.00000000E+00  0.00000000E+00  6.78000000E-01
 9.14000000E-01


 ========== ATMOSPHERE =========================================================

  ( EXTERNAL ATMOSPHERE FROM TABLE  7 AS FITTED )

  HEIGHT H IN KM GIVES THICKNESS OF ATMOSPHERE T IN G/CM**2
 H =   0.0...  3.0 KM ---> T = -1.44838E+02 + 1.1923E+03 * EXP( -H / 9.9419E+00)
 H =   3.0... 10.0 KM ---> T = -1.24071E+02 + 1.1740E+03 * EXP( -H / 9.6753E+00)
 H =  10.0... 50.0 KM ---> T =  3.60027E-01 + 1.4121E+03 * EXP( -H / 6.3614E+00)
 H =  50.0...100.0 KM ---> T = -8.24761E-04 + 8.1068E+02 * EXP( -H / 7.3564E+00)
 H = 100.0...109.2 KM ---> T =  2.21589E-03 - H / 4.9296E+04


 ========== SHOWER PARAMETERS ==================================================
 PRIMARY PARTICLE IDENTIFICATION IS            1
 PRIMARY ENERGY IS FIXED AT            1.0000E+04 GEV
 HEIGHT OF FIRST INTERACTION IS CHOSEN RANDOMLY
 TARGET OF FIRST INTERACTION IS AT RANDOM
 STARTING ALTITUDE AT     10923492.64 CM (=      0.00 G/CM**2)
                  WHICH IS AT TOP OF ATMOSPHERE

 OBSERVATION LEVEL # IN  CM    AND IN   G/CM**2 
          1       4.10000000E+05       6.44344752E+02
 THRESHOLD FOR THINNING SMALLER THAN ELCUT
  
 PARTICLES WITH LORENTZ FACTOR LARGER THAN     1.0000E+04 ARE PRINTED OUT
 SHOWER PARTICLES ENERGY CUT :
      FOR HADRONS   :      3.0000E-01 GEV
      FOR MUONS     :      3.0000E-01 GEV
      FOR ELECTRONS :      3.0000E-03 GEV
      FOR GAMMAS    :      3.0000E-03 GEV


 EARTH MAGNETIC FIELD STRENGTH IS    47.8015976      MICROTESLA
      WITH INCLINATION ANGLE         64.7376328      DEGREES
 MULTI-CORE #           1 :    0.0000000000000000        0.0000000000000000     
 CHERENKOV TELESCOPE #           1 , POS=   0.0000000000000000        0.0000000000000000        0.0000000000000000      , R=   200.00000000000000      , ID=           0
 CHERENKOV WAVELENGTH BAND FROM    250.000000      TO    700.000000      NANOMETER
 CHERENKOV BUNCH SIZE IS SET TO    1.0000000000000000     
 CHERENKOV PHOTONS ARE WRITTEN TO SEPARATE FILE WITH WAVELENGTH INSTEAD OF WEIGHT
 ATMOSPHERIC EXTINCTION COEFF. FOR CERENKOV PHOTONS, 180-700nm, in STEPS of 1 km  IS READ

 NUMBER OF SHOWERS TO GENERATE =           1



 ========== SHOWER NO          1 ===============================================

 PRESENT TIME : 05.06.2019  10:41:12 UTC
 AND RANDOM NUMBER GENERATOR AT BEGIN OF EVENT :       1
 SEQUENCE =  1  SEED =     62748  CALLS =         0  BILLIONS =         0
 SEQUENCE =  2  SEED =     62749  CALLS =         0  BILLIONS =         0
 SEQUENCE =  3  SEED =     62750  CALLS =         0  BILLIONS =         0
 SEQUENCE =  4  SEED =     62751  CALLS =         0  BILLIONS =         0
 SEQUENCE =  5  SEED =     62752  CALLS =         0  BILLIONS =         0

 TRACKING STARTS AT MARGIN OF ATMOSPHERE

 PRIMARY PARTICLE IS ELECTROMAGNETIC

                    TYPE       GAMMA    COSTHETA      PHIX      PHIY      HEIGHT      TIME       X-CM       Y-CM       GEN/CHI    WEIGHT   ALEVEL  E ON STACK

  FIRST INTERACTION AT    25116.097500000000       M ALTITUDE

 NO MORE SECONDARIES FOUND ON STACK
         218 PARTICLES WRITTEN TO STACK
         218 PARTICLES READ FROM STACK
           0 BLOCKS MAXIMUM STACK FILE SIZE
           0 SHIFTS TO EXTERNAL STACK
       18579 PARTICLES WRITTEN TO MPATAP
      253062 CHERENKOV BUNCHES WRITTEN TO MCETAB


 ENERGY - MULTIPLICITY MATRIX OF SHOWER NO          1
 ENERGY RUNS VERTICALLY, MULTIPLICITY HORIZONTALLY

Run time between 5min and 10min.