JToT.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#     Utility script to plot time slewing from JDETECTOR::JPMTDefaultSignaProcessor.
#
# ------------------------------------------------------------------------------------------


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


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable   DEBUG        3                  
set_variable   WORKDIR      /tmp/              
set_variable   OUTPUT_FILE  $WORKDIR/tot.root  
set_variable   NPE          1                  
set_variable:  FORMAT       GRAPHICS_FORMAT        gif
set_variable+  BATCH        GRAPHICS_BATCH         -B


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

case $# in
    1) set_variable  OUTPUT_FILE  $1;;
esac


if (( 1 )); then

    $JPP_DIR/examples/JDetector/JToT \
        -o $OUTPUT_FILE    \
        -n 999999          \
        -N $NPE            \
        -P %.gain=1.0      \
        -d $DEBUG

fi

if (( 1 )); then
    
    JPlot1D \
        -f ${OUTPUT_FILE}:h0            \
        -f ${OUTPUT_FILE}:h1            \
        -f ${OUTPUT_FILE}:h2            \
        -x "0 50"                       \
        -\> "time over threshold [ns]"  \
        -\^ "probability"               \
        -S 0.4                          \
        -y "1e-3 0.2" -Y                \
        -o $WORKDIR/tot.$FORMAT         $BATCH

fi

if (( 1 )); then
    
    JSum1D \
        -f ${OUTPUT_FILE}:h0            \
        -f ${OUTPUT_FILE}:h1            \
        -f ${OUTPUT_FILE}:h2            \
        -o $WORKDIR/sum.root

    JPlot1D \
        -f $WORKDIR/sum.root:h0         \
        -f $WORKDIR/sum.root:h1         \
        -f $WORKDIR/sum.root:h2         \
        -x "0 50"                       \
        -\> "time over threshold [ns]"  \
        -\^ "integrated probability"    \
        -S 0.4                          \
        -G XY                           \
        -y "1e-3 1.1" -Y                \
        -o $WORKDIR/tot-sum.$FORMAT     $BATCH

fi

if (( 1 )); then
    
    JPlot1D \
        -f ${OUTPUT_FILE}:h3            \
        -\> "transition time [ns]"      \
        -\^ "probability [a.u.]"        \
        -o $WORKDIR/t_ns.$FORMAT        $BATCH

fi