JTriggerEfficiency1D.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example histogram to histogram trigger efficiency.
#
# ------------------------------------------------------------------------------------------


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) set_variable  WORKDIR     $2;&
    1) set_variable  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 ""