plot-hobbit.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

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

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable     DEBUG                  ${CALIBRATE_DEBUG:-1}
set_variable     WORKDIR                /tmp/
    
if ( do_usage $* ); then
    usage "$script (input file)+"
fi

if (( $# < 1 )); then
    fatal "Not enough arguments."
fi

INPUT_FILES=($argv[1,-1])

if (( 1 )); then
    
    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:hq/g'`

    JPlot1D \
        -f "$HISTOGRAMS"              \
   -w 1000x500                   \
   -y "0 2"                      \
   -> "module index"            \
   -\^ "status"                  \
   -L TR -T ""                   \
   -o q.gif

fi

if (( 1 )); then
    
    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:h0/g'`

    JPlot1D \
        -f "$HISTOGRAMS"              \
   -w 1000x500                   \
        -y "-5 +5"                    \
        -> "module index"            \
        -\^ "#Deltat0 [ns]"           \
        -L TR -T ""                   \
   -o t0.gif

fi

if (( 1 )); then

    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:h0/g'`

    JVariance1D \
        -f "$HISTOGRAMS"              \
   -o variance.root              \
   -x "101 -10 +10"              \
   -F 0

    JPlot1D \
        -f variance.root:h0           \
        -> "#Deltat0 [ns]"           \
        -\^ "number of modules"       \
        -T "" -s 1110                 \
   -o t0s.gif

fi

if (( 1 )); then
    
    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:hc/g'`

    JPlot1D \
        -f "$HISTOGRAMS"              \
   -w 1000x500                   \
   -y "0 5"                      \
   -> "module index"            \
   -\^ "#chi^{2}/NDF"            \
   -L TR -T ""                   \
   -o chi2.gif

fi
    
if (( 1 )); then

    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:h2/g'`

    JPlot1D \
   -f "$HISTOGRAMS"              \
   -P Y                          \
   -y "-50 +50"                  \
   -> "#Deltat"                 \
   -\^ "number of hits"          \
   -T ""                         \
   -o d1t.gif

fi

if (( 1 )); then

    HISTOGRAMS=`echo ${INPUT_FILES} | sed 's/\.root/\.root:h2/g'`

    JPlot2D \
   -f "$HISTOGRAMS"              \
   -O COLZ                       \
   -> "module index"            \
   -< "#Deltat"                 \
   -T ""                         \
   -o d2t.gif
fi