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