JDAQTimesliceSelector.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script for JDAQTimeslice.
#
# ------------------------------------------------------------------------------------------


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  histogram\[%\].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

check_input_file $INPUT_FILE

for TIMESLICE in JDAQTimeslice JDAQTimesliceL0 JDAQTimesliceL1 JDAQTimesliceL2 JDAQTimesliceSN; do

    echo Processing $INPUT_FILE with selection $TIMESLICE.

    $JPP_DIR/examples/JDAQ/JDAQTimesliceSelector \
        -f ${INPUT_FILE}                     \
        -o ${OUTPUT_FILE/\%/${TIMESLICE}}    \
        -C $TIMESLICE                        \
        -d $DEBUG
done

JPlot1D \
    -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h0    \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h0  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h0  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h0  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h0  \
    -\> "PMT"                                \
    -\^ "number of events [a.u.]"            \
    -T "" -L TR

JPlot1D \
    -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h1    \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h1  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h1  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h1  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h1  \
    -\> "time-over-threshold [ns]"           \
    -\^ "number of events [a.u.]"            \
    -T "" -L TR

JPlot1D \
    -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h2    \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h2  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h2  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h2  \
    -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h2  \
    -\> "time [ns]"                          \
    -\^ "number of events [a.u.]"            \
    -w 800x400                               \
    -T "" -L TR