Back to Home

SIMUPERL User's Guide

Last update: June 21, 2005

  1. Installation
  2. Steering parameters
  3. Booking F2K files
  4. Dumping data
  5. Data storage
  6. Event viewer
  7. References
    1. 1. Installation

      The first step to install SIMUPERL is to download install_simu.sh. Exceute the file in the directory where you want to install. The script will automatically download all necessary files and install them. Download Magnesium, Aluminum release of AMASIM.

      2. Steering parameters

      To simulate events, run master.pl located at 'simu2002/simuperl/'. For example,

      ./master.pl -jobs 50 -max 1000 -priority 100 -start 1 -generator nusim

      This will submit 50 jobs at a time and -max refers to the total number of jobs. Just type master.pl and enter. It will show you more options. Several global variables defined in master.pl can be overwritten in the command line. One may need to change following variables in master.pl:

      $ENV{CLUSTER}
      $ENV{USER}
      $ENV{DIR}
      $ENV{DET}
      $ENV{GENERATOR}
      $ENV{YEAR}
      $ENV{OUTDIR}
      (See Sec. 5)
      ...

      3. Booking F2K files

      ./book-ic.pl muon.ama.mmc.nusim.numu.000100.simu2005.icecube.gz muon.ama.mmc.nusim.numu.000100.simu2005.icecube.hbook

      book-ic.pl can be found at 'simu2002/simuperl/'. One may need to change $muntsrc in book-ic.pl if it is not working. Once booked, do chain hbook files by using a simple script (an example). One can also use queuejobs.pl for booking on parallel machines. But if you get an error message like "bash: /bin/ls: Argument list too long", then use the modified code (squeuejobs.pl). Once scripts are generated, submit the jobs by using foreach:

      $ foreach y (script_name.*)
      $ foreach? /usr/pbs/bin/qsub $y
      $ foreach? end

      4. Dumping data

      The size of data is huge, so it is not easy to handle the data. Thus, we want to dump some useful variables and play with them. One can write a kumac file to dump data (an example). Then the ASCII data need to be booked (an example originally written by Gary Hill). For neutrino events, I found that atmswght is sometimes too small (~10-40), and got an error message saying "expected real number". So, I made a cut (atmswght=0.or.atmswght>1.e-35) in the example kumac file. Then, I faced another problem - the cut results in blank lines in the dumped data files. One can easily remove the blank lines by

      $ awk 'NF>0' < input.dat > output.dat

      5. Data storage

      See this page for the proposed structure of simulation data storage. In order to set output directory, do 'ln -s $DATAHOME/sim data', and then set $ENV{OUTDIR} in master.pl to "data/IceCube/2005/amasim/IC4" or so instead of "results".

      6. Event viewer

      eview is 3-D event viewer for AMANDA/ICECUBE.

      7. References

      * If you have any questions about content of this page, email me at csong@icecube.wisc.edu.