JSirene.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                            Auxiliary script to run JSirene.
#
# ------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

set_variable:  DEBUG            SIRENE_DEBUG                             1
set_variable+  KEEP             SIRENE_KEEP                              -k
set_variable+  WRITE_EMSHOWERS  SIRENE_WRITE_EMSHOWERS                   -s
set_variable:  NUMBER_OF_HITS   SIRENE_NUMBER_OF_HITS                    1
set_variable:  SEED             SIRENE_SEED                              0
set_variable:  FACTOR           SIRENE_FACTOR                            1.0
set_variable:  CDF              SIRENE_CDF                               $JPP_DATA/I%p.dat

if do_usage $*; then
    usage "$script <detector file> (input file)+ <output file>"\
          "\nAuxiliary script for simulation of detector response."
fi

if (( $# < 3 )); then
    fatal "Wrong number of arguments."
fi

set_variable DETECTOR     $argv[1]
set_variable INPUT_FILE   $argv[2,-2]
set_variable OUTPUT_FILE  $argv[-1]
            
# Sanity checks

for TYPE in 1 2 5 6 13 14; do
    if [[ ! -f ${CDF/\%/${TYPE}} ]]; then
        error "Missing CDF files."
        error "Run: JMakePDF.sh"
        JMakePDF.sh -h  
        exit
    fi
done

timer_start

JSirene \
    -F $CDF              \
    -a $DETECTOR         \
    -f $INPUT_FILE       \
    -o $OUTPUT_FILE      \
    -@ "Tmax_ns=2.0;"    \
    -U $FACTOR           \
    $KEEP                \
    $WRITE_EMSHOWERS     \
    -N $NUMBER_OF_HITS   \
    -S $SEED             \
    -d $DEBUG --!

timer_stop
timer_print

JPrintMeta \
    -f $OUTPUT_FILE

JPrintSirene \
    -f $OUTPUT_FILE

JPrintTree \
    -f $OUTPUT_FILE