[[TOC]] == 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 }}} {{{#!comment === Patching === To make the `TELESCOPE` keyword work properly, in `TELOUT` (search for this around line 75380) {{{#!fortran #if __CERENKOV__ && !__IACT__ *-- Author : R. Ulrich IKP KIT Karlsruhe 22/10/16 C======================================================================= INTEGER FUNCTION TELOUT(XPOSCER,YPOSCER,UDIRCER,VDIRCER) C----------------------------------------------------------------------- C ROUTINE TO CHECK WETHER CHERENKOV PHOTON HITS TELESCOPE C THIS FUNCTION IS CALLED FROM CERENK. C ARGUMENTS: C XPOSCER = X-POSITION OF PHOTON AT TELESCOPE C YPOSCER = Y-POSITION OF PHOTON AT TELESCOPE C UDIRCER = DIRECTION OF PHOTON RELATIVE TO X-DIRECTION C VDIRCER = DIRECTION OF PHOTON RELATIVE TO Y-DIRECTION C----------------------------------------------------------------------- }}} you have to replace {{{#!fortran TELZ = -OBSLEV(NOBSLV) + CERTELZ(J) }}} in `TELOUT` in `src/corsika.F` by {{{#!fortran TELZ = CERTELZ(J) }}} otherwise there is a z-distance between the sphere defining which photons are written to the output file and the telescope position. For inclined photons this leads to photons recorded far away from the telescope position. }}} === 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 {{{#!Spoiler {{{#!bash 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 {{{#!fortran 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) {{{#!fortran #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. {{{#!comment === The secret of wavelength 999 and 9999 in eventio: === > **IACTEXT** The interface to theTELOUTfunction is extended by parameters describing theemitting particle. This extended information is stored as an additional photon bunch(after the normal one) with mass, charge, energy, and emission time replacing thecx,cy,photons, andzemfields, respectively, and are identified by a wavelength of 9999.The compact output format is disabled for making that possible. In addition, all particlesarriving at the CORSIKA observation level are included in the eventio format output file,in a photon-bunch like block identified by array and detector numbers 999. == The 700nm problem == This is were Corsika checks the limits {{{#!fortran #if __IACT__ && !__CEFFIC__ IF ( WAVLGL .LT. 100.D0 .OR. WAVLGU .GT. 2000.D0 * .OR. WAVLGL .GE. WAVLGU ) THEN #else IF ( WAVLGL .LT. 100.D0 .OR. WAVLGU .GT. 700.D0 * .OR. WAVLGL .GE. WAVLGU ) THEN #endif }}} This is were the data is read from the files {{{#!fortran C======================================================================= SUBROUTINE TPDINI( OBLECE ) C----------------------------------------------------------------------- [...] C READ THE TABLE OF ATMOSPHERIC EXTINCTION AND STORE IT IN THE C ATMABS MATRIX. [...] DO I = 1, 105 READ(MCERABS,21) WLT 21 FORMAT(I4) READ(MCERABS,*,ERR=995,END=995) (ATMABS(I,J), J=0,50) ENDDO [...] C CALCULATE THE ATM. EXTINCTION FOR OBS. LEVEL INTERPOLATING THE TABLE DO IWL = 1, 105 FX0 = ATMABS(IWL,X0) FX1 = ATMABS(IWL,X1) ATEOBS(IWL) = LINEAR(X,X0,X1,FX0,FX1) ENDDO [...] C READ THE QUANTUM EFFECIENCY OF THE PMT OPEN(UNIT=MCERQEF,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)// * 'quanteff.dat',STATUS='OLD',FORM='FORMATTED',ERR=996) READ(MCERQEF,20) TEXTQEF READ(MCERQEF,23,ERR=997,END=997) (QUAEFF(I),I = 1,105) [...] C READ THE MIRROR REFLECTIVITY OPEN(UNIT=MCERMIR,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)// * 'mirreff.dat',STATUS='OLD',FORM='FORMATTED',ERR=998) READ(MCERMIR,20) TEXTREF READ(MCERMIR,23,END=999,ERR=999) (MIRREF(I), I=1,105) }}} This is where the index or the y-value and the corresponding x-values are calculated for mirror and pde efficiency. {{{#!fortran C===================================================================== SUBROUTINE TELEFF( ABSORB ) C----------------------------------------------------------------------- [...] C CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS RIWL = 1 + INT( (WL-180.D0)/5.D0 ) WLI0 = RIWL*5 + 175 WLI1 = RIWL*5 + 175 + 5 }}} This is the atmospheric absorption {{{#!fortran *-- Author : V. de Souza Filho, Uni. Campinas 22/06/1999 C======================================================================= SUBROUTINE ATABSO( ABSORB ) C----------------------------------------------------------------------- [...] C CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS RIWL = 1 + INT( (WL-180.D0)/5.D0 ) WLI0 = RIWL*5 + 175 WLI1 = RIWL*5 + 175 + 5 }}} }}} == 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! {{{#!Spoiler {{{ 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: {{{#!Spoiler {{{ 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 }}} }}} Run time between 5min and 10min.