Changes between Version 40 and Version 41 of InstallingCorsika


Ignore:
Timestamp:
Dec 21, 2018, 2:29:44 PM (3 months ago)
Author:
tbretz
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • InstallingCorsika

    v40 v41  
    55=== Extraction ===
    66
    7 Get the latest version of Corsika from here (right now this is at least 7.64):
    8 https://www.ikp.kit.edu/corsika/
     7Get 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.
    98
    109Extract Corsika
    1110{{{
    12 [0] tar xvf corsika-76400.tar.bz2
    13 }}}
    14 
     11[0] tar xvf corsika-76900.tar.bz2
     12}}}
     13
     14{{{#!comment
    1515=== Patching ===
    1616
     
    4545
    4646otherwise 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.
     47}}}
    4748
    4849=== Setup and Compilation ===
     
    5051Now, call coconut
    5152{{{
    52 [1] cd corsika-76400
     53[1] cd corsika-76900
    5354[2] ./coconut
    5455}}}
     
    5657You can use the defaults (just press return) for
    5758
    58 * Compile in 32 or 64bit mode? [64 bit]
    59 * Which high energy hadronic interaction model do you want to use? [QGSJET 01C]
    60 * Which low energy hadronic interaction model do you want to use? [GHEISHA 2002d]
    61 * Which detector geometry do you have? [horizontal flat]
     59* Compile in 32 or 64bit mode? [2: 64 bit]
     60* Which high energy hadronic interaction model do you want to use? [4: QGSJET 01C]
     61* Which low energy hadronic interaction model do you want to use? [1: GHEISHA 2002d]
     62* Which detector geometry do you have? [1: horizontal flat]
    6263
    6364Note that so far we were always using FLUKA: WHY??? Shell we use FLUKA as well?
     
    7677
    7778Now you will get asked more question, answer them as follows
    78 * Cherenkov light vertical (longitudinal) distribution option? [No Cherenkov light distribution at all]
    79 * Do you want Cherenkov light emission angle wavelength dependence? [depending on wavelength] {This corresponds to the CERWLEN option}
    80 
    81 Now, coconut can create the input files, configure the make system and make. To proceed, select  `*** Finish selection ***`.
     79* Cherenkov light vertical (longitudinal) distribution option? [3: No Cherenkov light distribution at all]
     80* Do you want Cherenkov light emission angle wavelength dependence? [2: depending on wavelength] {This corresponds to the CERWLEN option}
     81
     82Now, coconut can create the input files, configure the make system and make. To proceed, select  `*** Finish selection ***` (z).
    8283
    8384You should see
    8485{{{
    8586  Are you sure you want to continue with these current option selection:
    86     CERENKOV NOCLONG THIN TRAJECT VIEWCONE ATMEXT BERNLOHRDIR CERWLEN CEFFIC
     87     TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
    8788     yes or no ? (default: yes) >
    8889}}}
     
    9495{{{#!Spoiler
    9596{{{#!bash
     97tbretz@d4m:~/SW/corsika-76900> coconut
    9698 
    9799
     
    174176
    175177---------------------------------------------------------------------------
    176   options:   GHEISHA TIMEAUTO QGSJET01 HORIZONTAL
     178  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO
    177179
    178180 Which additional CORSIKA program options do you need ?
     
    264266    SELECTED         : VIEWCONE
    265267    NOT COMPATIBLE TO: VOLUMECORR INTTEST
    266     SELECTED path to installation: "/home/tbretz/SW/corsika-764xx/bernlohr"
     268    SELECTED path to installation: "/home/tbretz/SW/corsika-76900/bernlohr"
    267269    SELECTED         : ATMEXT
    268270    NOT COMPATIBLE TO: INTTEST
     
    270272
    271273---------------------------------------------------------------------------
    272   options:   GHEISHA TIMEAUTO QGSJET01 HORIZONTAL
    273   selection: BERNLOHRDIR CEFFIC THIN VIEWCONE TRAJECT NOCLONG CERWLEN CERENKOV ATMEXT
     274  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO
     275  selection: TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
    274276
    275277 Which additional CORSIKA program options do you need ?
     
    324326    x - exit make
    325327
    326     (multiple selections accepted, leading '-' removes option):
     328    (multiple selections accepted, leading '-' removes option): z
    327329  Are you sure you want to continue with these current option selection:
    328     BERNLOHRDIR CEFFIC THIN VIEWCONE TRAJECT NOCLONG CERWLEN CERENKOV ATMEXT
     330    TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
    329331    yes or no ? (default: yes) >
    330332
     
    333335  Your final selection to build CORSIKA is:
    334336
    335   options:   GHEISHA TIMEAUTO QGSJET01 HORIZONTAL
    336   selection: BERNLOHRDIR CEFFIC THIN VIEWCONE TRAJECT NOCLONG CERWLEN CERENKOV ATMEXT
     337  options:   HORIZONTAL QGSJET01 GHEISHA TIMEAUTO
     338  selection: TRAJECT VIEWCONE THIN CERENKOV CEFFIC BERNLOHRDIR CERWLEN ATMEXT NOCLONG
    337339
    338340
     
    385387checking how to run the C preprocessor... cpp
    386388checking how to get verbose linking output from gfortran... -v
    387 checking for Fortran 77 libraries of gfortran...  -L/usr/lib/gcc/x86_64-linux-gnu/5 -L/usr/lib/gcc/x86_64-linux-gnu/5/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/5/../../../../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/5/../../.. -lgfortran -lm -lquadmath
     389checking 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
    388390checking how to print strings... printf
    389391checking for a sed that does not truncate output... /bin/sed
     
    476478checking for thinning... (cached) yes
    477479checking for Cherenkov grid detector... (cached) yes
    478 checking for bernlohr... /home/tbretz/SW/corsika-764xx/bernlohr
     480checking for bernlohr... /home/tbretz/SW/corsika-76900/bernlohr
    479481checking for Cherenkov telescope detector... (cached) no
    480482checking new Cherenkov counting... (cached) no
     
    575577
    576578
    577 Compile CORSIKA in "/home/tbretz/SW/corsika-764xx/lib/unknown" and copy executable in "/home/tbretz/SW/corsika-764xx/run"
     579Compile CORSIKA in "/home/tbretz/SW/corsika-76900/lib/unknown" and copy executable in "/home/tbretz/SW/corsika-76900/run"
    578580
    579581Making install in bernlohr
    580 make[1]: Entering directory '/home/tbretz/SW/corsika-764xx/bernlohr'
     582make[1]: Entering directory '/home/tbretz/SW/corsika-76900/bernlohr'
    581583cc -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
    582584mv -f .deps/libbern_a-fileopen.Tpo .deps/libbern_a-fileopen.Po
     
    597599ar: `u' modifier ignored since `D' is the default (see `U')
    598600ranlib libatmext.a
    599 make[2]: Entering directory '/home/tbretz/SW/corsika-764xx/bernlohr'
    600  /bin/mkdir -p '/home/tbretz/SW/corsika-764xx/lib/unknown'
    601  /usr/bin/install -c -m 644  libbern.a libatmext.a '/home/tbretz/SW/corsika-764xx/lib/unknown'
    602  ( cd '/home/tbretz/SW/corsika-764xx/lib/unknown' && ranlib libbern.a )
    603  ( cd '/home/tbretz/SW/corsika-764xx/lib/unknown' && ranlib libatmext.a )
     601make[2]: Entering directory '/home/tbretz/SW/corsika-76900/bernlohr'
     602 /bin/mkdir -p '/home/tbretz/SW/corsika-76900/lib/unknown'
     603 /usr/bin/install -c -m 644  libbern.a libatmext.a '/home/tbretz/SW/corsika-76900/lib/unknown'
     604 ( cd '/home/tbretz/SW/corsika-76900/lib/unknown' && ranlib libbern.a )
     605 ( cd '/home/tbretz/SW/corsika-76900/lib/unknown' && ranlib libatmext.a )
    604606make[2]: Nothing to be done for 'install-data-am'.
    605 make[2]: Leaving directory '/home/tbretz/SW/corsika-764xx/bernlohr'
    606 make[1]: Leaving directory '/home/tbretz/SW/corsika-764xx/bernlohr'
     607make[2]: Leaving directory '/home/tbretz/SW/corsika-76900/bernlohr'
     608make[1]: Leaving directory '/home/tbretz/SW/corsika-76900/bernlohr'
    607609Making install in src
    608 make[1]: Entering directory '/home/tbretz/SW/corsika-764xx/src'
     610make[1]: Entering directory '/home/tbretz/SW/corsika-76900/src'
    609611cpp  -DHAVE_CONFIG_H -I../include -traditional-cpp corsika.F > corsikacompilefile.f
    610612gfortran -O0  -g -std=legacy -c -o corsika-corsikacompilefile.o `test -f 'corsikacompilefile.f' || echo './'`corsikacompilefile.f
     
    617619cc -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 &&\
    618620mv -f $depbase.Tpo $depbase.Po
    619 /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-764xx/lib/unknown -lbern -latmext     
    620 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-764xx/lib/unknown -lbern -latmext
    621 make[2]: Entering directory '/home/tbretz/SW/corsika-764xx/src'
    622  /bin/mkdir -p '/home/tbretz/SW/corsika-764xx/run'
    623   /bin/bash ../libtool   --mode=install /usr/bin/install -c corsika '/home/tbretz/SW/corsika-764xx/run'
    624 libtool: install: /usr/bin/install -c corsika /home/tbretz/SW/corsika-764xx/run/corsika
     621/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     
     622libtool: 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
     623make[2]: Entering directory '/home/tbretz/SW/corsika-76900/src'
     624 /bin/mkdir -p '/home/tbretz/SW/corsika-76900/run'
     625  /bin/bash ../libtool   --mode=install /usr/bin/install -c corsika '/home/tbretz/SW/corsika-76900/run'
     626libtool: install: /usr/bin/install -c corsika /home/tbretz/SW/corsika-76900/run/corsika
    625627make[2]: Nothing to be done for 'install-data-am'.
    626 make[2]: Leaving directory '/home/tbretz/SW/corsika-764xx/src'
    627 make[1]: Leaving directory '/home/tbretz/SW/corsika-764xx/src'
     628make[2]: Leaving directory '/home/tbretz/SW/corsika-76900/src'
     629make[1]: Leaving directory '/home/tbretz/SW/corsika-76900/src'
    628630Making install in .
    629 make[1]: Entering directory '/home/tbretz/SW/corsika-764xx'
    630 make[2]: Entering directory '/home/tbretz/SW/corsika-764xx'
     631make[1]: Entering directory '/home/tbretz/SW/corsika-76900'
     632make[2]: Entering directory '/home/tbretz/SW/corsika-76900'
    631633make  install-exec-hook
    632 make[3]: Entering directory '/home/tbretz/SW/corsika-764xx'
     634make[3]: Entering directory '/home/tbretz/SW/corsika-76900'
    633635 
    634636 
    635637 
    636 --> "corsika76400Linux_QGSJET_gheisha_thin" successfully installed in :
    637     /home/tbretz/SW/corsika-764xx/run/
     638--> "corsika76900Linux_QGSJET_gheisha_thin" successfully installed in :
     639    /home/tbretz/SW/corsika-76900/run/
    638640 
    639 --> You can run CORSIKA in /home/tbretz/SW/corsika-764xx/run/ using for instance :
    640     ./corsika76400Linux_QGSJET_gheisha_thin < all-inputs > output.txt
     641--> You can run CORSIKA in /home/tbretz/SW/corsika-76900/run/ using for instance :
     642    ./corsika76900Linux_QGSJET_gheisha_thin < all-inputs > output.txt
    641643 
    642644 
    643645 
    644 make[3]: Leaving directory '/home/tbretz/SW/corsika-764xx'
     646make[3]: Leaving directory '/home/tbretz/SW/corsika-76900'
    645647make[2]: Nothing to be done for 'install-data-am'.
    646 make[2]: Leaving directory '/home/tbretz/SW/corsika-764xx'
    647 make[1]: Leaving directory '/home/tbretz/SW/corsika-764xx'
     648make[2]: Leaving directory '/home/tbretz/SW/corsika-76900'
     649make[1]: Leaving directory '/home/tbretz/SW/corsika-76900'
    648650}}}
    649651}}}
     
    717719
    718720=== A problem with the IACT option might be this: ===
    719 > m)  Starting with version 1.38, the dynamic range in a telescope is basically unlimited dueto 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.
    720 
     721> 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.
     722
     723{{{#!comment
    721724=== The secret of wavelength 999 and 9999 in eventio: ===
    722725> **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.
     
    801804      WLI1 = RIWL*5 + 175 + 5
    802805}}}
     806}}}
    803807
    804808== MMCS ==
     
    815819Enabling 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.
    816820
    817 The disadvantage of the IACT option is that teh 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.
     821The 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.
    818822
    819823Using IACT means that we have to understand the automatic thinning and implement it and to make sure eventio works.