JDAQEvent.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script for JDAQEvent.
#
# ------------------------------------------------------------------------------------------


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/
define_variable  INPUT_FILE
set_variable     OUTPUT_FILE  $WORKDIR/histogram.root
#set_variable     TRIGGER_MASK    10    # 3D Shower
#set_variable     TRIGGER_MASK   100    # MX Shower
#set_variable     TRIGGER_MASK 10000    # 3D Muon
set_variable     TRIGGER_MASK 11111    # All

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/JDAQEvent \
        -f $INPUT_FILE   \
        -o $OUTPUT_FILE  \
        -T $TRIGGER_MASK \
        -d $DEBUG --!
    
    timer_stop
    timer_print
fi

if (( 1 )); then
    JPlot1D \
        -f ${OUTPUT_FILE}:hn            \
        -f ${OUTPUT_FILE}:hm            \
        -\> "hits"                      \
        -\^ "number of events [a.u.]"   

    JPlot1D \
        -f ${OUTPUT_FILE}:h0            \
        -\> "PMT"                       \
        -\^ "number of events [a.u.]"   
    
    JPlot1D \
        -f ${OUTPUT_FILE}:h1            \
        -\> "time-over-threshold [ns]"  \
        -\^ "number of events [a.u.]"   \
        
    JPlot1D \
        -f ${OUTPUT_FILE}:h2            \
        -\> "time [ns]"                 \
        -\^ "number of events [a.u.]"   \
        -w 800x400
fi