multi-plotk40.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                         Example script to plot JFitK40 results.
#
# ------------------------------------------------------------------------------------------

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

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

if ( do_usage $* ); then
    usage "$script <detector file> <input file>"
fi

case $# in
    2)  set_variable  DETECTOR     $1
        set_variable  INPUT_FILE   $2;;
    *)  fatal "Invalid number of arguments"
esac

set_variable  TYPE    F     # histogram type

eval `JPrintDetector -a $DETECTOR -O SUMMARY`

attach getModule -a $DETECTOR

typeset -Z 3 STRING
typeset -Z 2 FLOOR

for STRING in $STRINGS[*]; do

    for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do

        print -u $FD_O "$STRING $FLOOR"

        read MODULE <& $FD_I

        echo "($STRING,$FLOOR) $MODULE"
            
        JPlot2D \
            -f ${INPUT_FILE}:${MODULE}.2${TYPE}    \
            -z "-0.1 1.5"                          \
            -\> "PMT pair"                         \
            -\< "#Deltat [ns]"                     \
            -\^ "Rate [Hz]"                        \
            -T "($STRING,$FLOOR)"                  \
            -O COLZ                                \
            -o /tmp/${TYPE}_${STRING}_${FLOOR}.gif \
            -B
            
    done

    montage \
        -tile 6x3                       \
        -geometry +0+0                  \
        /tmp/${TYPE}_${STRING}_*.gif    \
        ${TYPE}_${STRING}.gif

done

detach