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

set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

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`

kill_child_processes_at_exit

attach getModule -a $DETECTOR

typeset -Z 4 STRING
typeset -Z 2 FLOOR

for STRING in $STRINGS[*]; do

    set_variable  FIRST_FLOOR  1

    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}.$FORMAT    $BATCH
            
    done

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

done

detach