Table of Contents
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 (excerpt):
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
[...]
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)
[...]
CONSTANTS FOR KAONS CKA(1) TO CKA(80)
[...]
CONSTANTS FOR ETAS CETA(1) TO CETA(5)
[...]
CONSTANTS FOR STRANGE BARYONS CSTRBA(1) TO CSTRBA(11)
[...]
========== 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
[...]
WORDS WRITTEN TO PARTICLE DATA FILE UP TO NOW = 144144
PARTICLES AT DETECTOR LEVEL:
FOR LEVEL 1
HEIGHT IN CM 4.100000E+05
HEIGHT IN G/CM**2 6.443448E+02
NO OF PROTONS = 0.
NO OF ANTIPROTONS = 0.
NO OF NEUTRONS = 0.
NO OF ANTINEUTRONS = 0.
NO OF GAMMAS = 1.596500E+04
NO OF POSITRONS = 9.870000E+02
NO OF ELECTRONS = 1.623000E+03
NO OF MU + = 2.
NO OF MU - = 2.
NO OF PI 0 = 0.
NO OF PI + = 0.
NO OF PI - = 0.
NO OF K0L = 0.
NO OF K0S = 0.
NO OF K + = 0.
NO OF K - = 0.
NO OF STR. BARYONS = 0.
NO OF DEUTERONS = 0.
NO OF TRITONS = 0.
NO OF 3HELIUM = 0.
NO OF ALPHAS = 0.
NO OF OTHER PARTIC.= 0.
NO OF DECAYED MUONS= 9.
NO OF ELIMIN. MUONS= 4.
NO OF NUCLEONS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF PIONS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF ETAS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF KAONS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF S.BARYONS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF RHO MESNS PRODUCED IN FIRST HADR. INTERACTION = 0
NO OF OTH. HADR.PRODUCED IN FIRST HADR. INTERACTION = 0
TOTAL MULTIPLICITY OF FIRST HADR. INTERACTION = 0
ELASTICITY OF FIRST HADR. INTERACTION = 0.0000
VERTICAL DEPTH (G/CM**2) OF FIRST HADR. INTERACTION = 0.0000
VERTICAL ALTITUDE (CM) OF FIRST HADR. INTERACTION =10503299.1
CROSS-SECTION MILLIBARN OF FIRST HADR. INTERACTION = 0.0000
TARGET MASS NUMBER OF FIRST HADR. INTERACTION = 0.0000
---------- NKG - OUTPUT OF SHOWER NO 1 --------------------------------
ELECTRON/GAMMA THRESHOLD AT 0.00300 / 0.00300 GEV
LEVEL THICKNESS HEIGHT ELECT. NUMBER PSEUDO-
NO. G/CM**2 CM AGE
1 100. 1681722. 155.354 0.314
2 200. 1254344. 2060.215 0.567
3 300. 984788. 6541.308 0.778
4 400. 778965. 9595.098 0.954
5 500. 610681. 8695.883 1.104
6 600. 467251. 5685.975 1.235
7 644. 410000. 4347.731 1.288
RADIAL BIN DISTANCE(CM) LOCAL AGE AT LEVEL NO. 1 AT HEIGHT: 410000. CM
1 229. 1.185
2 389. 1.154
3 661. 1.105
4 1123. 1.043
5 1908. 0.965
6 3241. 0.896
7 5506. 0.858
8 9353. 0.870
9 15887. 0.927
LATERAL ELECTRON DENSITY (/CM**2) AT LEVEL NO. 1 AT HEIGHT: 410000. CM
-----------------------------------------------------------------------------
DIST (CM) CZX CZXY CZY CZYX
-20000. 1.90878E-07 1.91545E-07 1.90859E-07 1.90948E-07
[...]
20000. 1.92439E-07 1.91024E-07 1.90956E-07 1.91253E-07
CERCNT = 235973.109
CHERENKOV PH. FROM ELECTRONS = 2.3597310E+05 CHERENKOV PH. FROM HADRONS = 0.0000000E+00
END OF SHOWER NO 1
PRESENT TIME : 05.06.2019 10:47:37 UTC
AND RANDOM NUMBER GENERATOR AT END OF EVENT : 1
SEQUENCE = 1 SEED = 62748 CALLS = 2020 BILLIONS = 0
SEQUENCE = 2 SEED = 62749 CALLS = 39537543 BILLIONS = 0
SEQUENCE = 3 SEED = 62750 CALLS = 612258645 BILLIONS = 2
SEQUENCE = 4 SEED = 62751 CALLS = 0 BILLIONS = 0
SEQUENCE = 5 SEED = 62752 CALLS = 0 BILLIONS = 0
========== RUN SUMMARY ========================================================
NUMBER OF GENERATED EVENTS = 1
TOTAL TIME USED = 385. SEC
TIME PER EVENT = 385.00 SEC
TOTAL SPACE ON MPATAP USED = 150696 WORDS
SPACE PER EVENT ON MPATAP = 150696 WORDS
AVERAGE HEIGHT OF 1ST INT. = 26.757 +- 0.000 G/CM**2
PRESENT TIME : 05.06.2019 10:47:37 UTC
========== END OF RUN ================================================
Run time between 5min and 10min.
Last modified
6 years ago
Last modified on 06/05/19 13:17:27
Note:
See TracWiki
for help on using the wiki.
