multi-plotk40.sh

Go to the documentation of this file.

#!/usr/bin/env 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`
 
 
 
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 ${TMPDIR:-/tmp}/${TYPE}_${STRING}_${FLOOR}.$FORMAT    $BATCH
            
    done
 
    montage \
        -tile 6x3                          \
        -geometry +0+0                     \
        ${TMPDIR:-/tmp}/${TYPE}_${STRING}_*.$FORMAT   \
        ${TYPE}_${STRING}.$FORMAT
 
done