JConvertDusj.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Auxiliary script to convert files to DUSJ format.
#
# ------------------------------------------------------------------------------------------

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

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable     DEBUG        2
set_variable     WORKDIR      /tmp/

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

case $# in
    4) set_variable WORKDIR      $4;&
    3) set_variable OUTPUT_FILE  $3;
       set_variable DETECTOR     $2;
       set_variable INPUT_FILE   $1;;
    *) fatal "Wrong number of arguments."
esac

if [[ ${OUTPUT_FILE##*\.} == "evt" ]]; then

    JConvertDusj \
        -f $INPUT_FILE         \
        -a $DETECTOR           \
        -o $WORKDIR/dusj.evt   \
        -d $DEBUG

    awk '
    /primarylepton:/ { $1="track_in:"; print $0; next }
    /weights:/       { if (NF > 1)     print $0; next }
    /w2list:/        { if (NF > 1)     print $0; next }
    /w3list:/        { if (NF > 1)     print $0; next }
    //               {                 print $0; next }
    ' $WORKDIR/dusj.evt > $OUTPUT_FILE

    rm -f $WORKDIR/dusj.evt

else

    JConvertDusj \
        -f $INPUT_FILE         \
        -a $DETECTOR           \
        -o $OUTPUT_FILE        \
        -d $DEBUG

fi