JDAQSummaryslice.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script for JDAQSummaryslice.
#
# ------------------------------------------------------------------------------------------


if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable   WORKDIR      /tmp/         
set_variable   OUTPUT_FILE  summary.root  

if do_usage $*; then
    usage "$script <input file> [output file]"
fi

case $# in
    2) set_variable OUTPUT_FILE  $2;&
    1) set_variable INPUT_FILE   $1;;
    *) fatal "Wrong number of arguments."
esac

if (( 1 )); then

    print_variable     INPUT_FILE OUTPUT_FILE
    check_input_file   $INPUT_FILE

    timer_start
     
    $JPP_DIR/examples/JDAQ/JDAQSummaryslice \
        -f $INPUT_FILE   \
        -o $OUTPUT_FILE  \
        -d2
    
    timer_stop
    timer_print

fi

if (( 1 )); then 

    JPlot1D \
        -f ${OUTPUT_FILE}:h0            \
        -x "0 50"                       \
        -\> "Rate [kHz]"                \
        -\^ "number of events / kHz [a.u.]"   

    JPlot1D \
        -f ${OUTPUT_FILE}:h1            \
        -\> "high-rate veto PMT"        \
        -\^ "fraction of events [a.u.]"   

    JPlot1D \
        -f ${OUTPUT_FILE}:h2            \
        -\> "FIFO"                      \
        -\^ "fraction of events [a.u.]"   

    JPlot2D \
        -f ${OUTPUT_FILE}:hu            \
        -\> "sequence number"           \
        -\< "number of packets"         \
        -Z                              \
        -O COLZ

fi