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