JTriggerReprocessor.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script to re-process trigger.
#
# ------------------------------------------------------------------------------------------


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        /tmp/
set_variable  INPUT_FILE     $JPP_DATA/trigger_efficiency+background.root
set_variable  DETECTOR       $JPP_DATA/km3net_reference.detx
set_variable  TRIGGER        $JPP_DATA/trigger_parameters_arca.txt


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

case $# in 
    4) WORKDIR=$4;&
    3) INPUT_FILE=$3;&
    2) TRIGGER=$2;&
    1) DETECTOR=$1;;
esac


print_variable     DETECTOR TRIGGER INPUT_FILE OUTPUT_FILE
check_input_file   $DETECTOR $TRIGGER
check_output_file  $OUTPUT_FILE
check_variable     JPP_DIR JPP_BIN JPP_LIB ROOTSYS
check_process      JTriggerProcessor


if ( ! reuse_file ${WORKDIR}/trigger_reprocessor.root ); then

   timer_start
   
   JTriggerReprocessor \
       -a ${DETECTOR}                           \
       -f ${INPUT_FILE}                         \
       -o ${WORKDIR}/trigger_reprocessor.root   \
       -@ $TRIGGER                              \
       -@"trigger3DShower.enabled=0"            \
       -@"trigger3DMuon.numberOfHits=15"        \
       -d ${DEBUG}                              \
       --!

   timer_stop
   timer_print
fi

   
$JPP_DIR/examples/JTrigger/JVolume1D \
    -f $INPUT_FILE                           \
    -o ${WORKDIR}/volume\[original\].root    \
    -X                                       \
    -d $DEBUG


$JPP_DIR/examples/JTrigger/JVolume1D \
    -f ${WORKDIR}/trigger_reprocessor.root   \
    -o ${WORKDIR}/volume\[reprocessed\].root \
    -X                                       \
    -d $DEBUG



JPlot1D \
    -f ${WORKDIR}/volume\[original\].root:hV    \
    -f ${WORKDIR}/volume\[reprocessed\].root:hV \
    -y "1e-2 20"            \
    -XXY                    \
    -> "E [GeV]"           \
    -\^ "V [km^{3}]"        \
    -T "" -L TL             \
    -o reprocess.gif