JDAQTimesliceSelector.sh

Go to the documentation of this file.

#!/usr/bin/env 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      ${TMPDIR:-/tmp}/                
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