JTriggerEfficiency1D.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Utility script to test JEfficiency.
#
# ------------------------------------------------------------------------------------------


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        ./
set_variable  INPUT_FILE     $JPP_DATA/sirene.root
set_variable  DETECTOR       $JPP_DATA/km3net_reference.detx
set_variable  TRIGGER        $JPP_DATA/trigger_parameters_arca.txt
set_variable  BACKGROUND_HZ  0


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

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


set_variable  OUTPUT_FILE    $WORKDIR/trigger_efficiency.root

if (( 1 )); then

    if [[ -f $OUTPUT_FILE ]] then
   rm -i $OUTPUT_FILE
    fi

    if [[ ! -f $OUTPUT_FILE ]] then

   JTriggerEfficiency \
       -a ${DETECTOR}                \
       -f ${INPUT_FILE}              \
       -o ${OUTPUT_FILE}             \
       -@ ${TRIGGER}                 \
       -d ${DEBUG}                   \
       -B "$BACKGROUND_HZ"           \
       --!
    fi
fi

set_variable  INPUT_FILE     $OUTPUT_FILE
set_variable  OUTPUT_FILE    $WORKDIR/efficiency.root

if (( 1 )); then 

    $JPP_DIR/examples/JTrigger/JTriggerEfficiency1D \
   -f $INPUT_FILE   \
   -o $OUTPUT_FILE  \
   -d $DEBUG
    
fi


JPlot1D \
    -f ${OUTPUT_FILE}:h1    \
    -> "number of hits"    \
    -\^ "number of events"  \
    -o counts.gif -T ""


JPlot1D \
    -f ${OUTPUT_FILE}:h2    \
    -y "0 1.1"              \
    -> "number of hits"    \
    -\^ "effieciency"       \
    -o efficiency.gif -T ""