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         
set_variable:  FORMAT      GRAPHICS_FORMAT       gif
set_variable+  BATCH       GRAPHICS_BATCH        -B

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

case $# in 
    4) set_variable  WORKDIR     $4;&
    3) set_variable  INPUT_FILE  $3;&
    2) set_variable  TRIGGER     $2;&
    1) set_variable  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.$FORMAT    $BATCH