JMonteCarloFileWriter.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
# utility script for Sirene Monte Carlo.
version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                        Example script to test ASCII I/O of Monte Carlo files.
#
# ------------------------------------------------------------------------------------------


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


source $JPP_DIR/setenv.sh $JPP_DIR

set_variable  DEBUG             1
set_variable  WORKDIR           /tmp/
set_variable  INPUT_FILE        $JPP_DATA/genhen.km3net_wpd_V2_0.evt
set_variable  DETECTOR          $JPP_DATA/km3net_reference.detx
set_variable  NUMBER_OF_EVENTS  100
set_variable  SEED              12345678
    

if ( do_usage $* ); then
    usage "$script [detector file [input file [working directory]]]"
fi

case $# in 
    3) WORKDIR=$3;&
    2) INPUT_FILE=$2;&
    1) DETECTOR=$1;;
esac


if (( 1 )); then

    if [[ -f $WORKDIR/sirene.root ]]; then
   rm -i $WORKDIR/sirene.root
    fi
    
    if [[ ! -f $WORKDIR/sirene.root ]]; then

   JSirene \
       -F $JPP_DATA/I%p.dat     \
       -a $DETECTOR             \
       -f $INPUT_FILE           \
       -o $WORKDIR/sirene.root  \
       -n $NUMBER_OF_EVENTS     \
       -S $SEED                 \
       -d $DEBUG
    fi
fi

if (( 1 )); then

    if [[ -f $WORKDIR/sirene.evt ]]; then
   rm -i $WORKDIR/sirene.evt
    fi
    
    if [[ ! -f $WORKDIR/sirene.evt ]]; then

   JSirene \
       -F $JPP_DATA/I%p.dat     \
       -a $DETECTOR             \
       -f $INPUT_FILE           \
       -o $WORKDIR/sirene.evt   \
       -n $NUMBER_OF_EVENTS     \
       -S $SEED                 \
       -d $DEBUG            
    fi
fi


rm -f $WORKDIR/sirene_root.txt
rm -f $WORKDIR/sirene_evt.txt

echo "Print $WORKDIR/sirene.root"
JPrint \
    -f $WORKDIR/sirene.root          \
    -n $NUMBER_OF_EVENTS             \
    > $WORKDIR/sirene_root.txt

echo "Print $WORKDIR/sirene.evt"
JPrint \
    -f $WORKDIR/sirene.evt           \
    -n $NUMBER_OF_EVENTS             \
    > $WORKDIR/sirene_evt.txt


echo "Comparing $WORKDIR/sirene_root.txt and $WORKDIR/sirene_evt.txt"
diff \
    $WORKDIR/sirene_root.txt         \
    $WORKDIR/sirene_evt.txt