JConvertDusj.sh

Go to the documentation of this file.

#!/usr/bin/env 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  ${TMPDIR:-/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