JFilter.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test trigger algorithms.
#
# ------------------------------------------------------------------------------------------


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


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable     DEBUG           2
set_variable     WORKDIR         $JPP_DATA
set_variable     INPUT_FILE      $WORKDIR/trigger_efficiency+background.root
set_variable     OUTPUT_FILE     $WORKDIR/filter\%.root
set_variable     DETECTOR        $WORKDIR/km3net_reference.detx
set_variable     HISTOGRAM       1


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

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


if (( 1 )); then

    check_input_file  $DETECTOR $INPUT_FILE

    for OPTION in A B C D E F; do
    
   FILE="${OUTPUT_FILE/\%/[${OPTION}]}"
   
   print_variable DETECTOR INPUT_FILE FILE
   
   $JPP_DIR/examples/JTrigger/JFilter \
       -a $DETECTOR     \
       -f $INPUT_FILE   \
       -o $FILE         \
       -c $OPTION       \
       -T 20            \
       -R 200           \
       -H $HISTOGRAM    \
       -d $DEBUG
   
    done

fi


if (( 1 )); then

    X_LABEL[1]="Number of hits"
    X_LABEL[2]="E_{#nu} [GeV]"

    X_OPTION[1]=""
    X_OPTION[2]="-XX"

    X_RANGE[1]="-x 1 50"
    X_RANGE[2]="-x 2  7"

    FILE_E=""
    FILE_P=""

    for OPTION in A B C D E F; do
   FILE_E+=" ${OUTPUT_FILE/\%/[${OPTION}]}:he"
   FILE_P+=" ${OUTPUT_FILE/\%/[${OPTION}]}:hp"
    done

    FILE_E=${FILE_E## }
    FILE_P=${FILE_P## }

    JPlot1D \
   -f $FILE_E               \
   -y "0 1.1"               \
   -\^ "efficiency"         \
   -> $X_LABEL[$HISTOGRAM] \
   $X_OPTION[$HISTOGRAM]    \
   $X_RANGE[$HISTOGRAM]     \
   -L TR


    JPlot1D \
   -f $FILE_P               \
   -y "0 1.1"               \
   -\^ "purity"             \
   -> $X_LABEL[$HISTOGRAM] \
   $X_OPTION[$HISTOGRAM]    \
   $X_RANGE[$HISTOGRAM]     \
   -L TR

fi