plot-equalizer.sh

Go to the documentation of this file.

#!/bin/zsh
#
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
set_variable   DEBUG       2      
set_variable   WORKDIR     `pwd`  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> (input file)+"
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable     DETECTOR               $argv[1]
set_array        INPUT_FILES            $argv[2,-1]
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
kill_child_processes_at_exit
 
attach getModule -a $DETECTOR
 
for STRING in $STRINGS[*]; do
 
    rm -f ${TMPDIR:-/tmp}/H1-${(l:4::0::0:)STRING}_*.$FORMAT
 
    echo -n "Creating graphics for string ${(l:4::0::0:)STRING}:"
 
    set_variable  FIRST_FLOOR  1
 
    for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
        echo -n " ${FLOOR}"
 
        sput $STRING $FLOOR
        sget  MODULE
 
        JPlot1D \
            -f "${INPUT_FILES//root/root:H1\[${MODULE}\]}"                   \
            -y "0.0 10.0"                                                    \
            -> "PMT"                                                        \
            -\^ "rate [Hz]"                                                  \
            -T "(${(l:4::0::0:)STRING},${(l:2::0::0:)FLOOR})"                \
            -o ${TMPDIR:-/tmp}/H1-${(l:4::0::0:)STRING}_${(l:2::0::0:)FLOOR}.$FORMAT    $BATCH
    done
 
    montage \
        -tile 6x3                               \
        -geometry +0+0                          \
        ${TMPDIR:-/tmp}/H1-${(l:4::0::0:)STRING}_*.$FORMAT     \
        H1-${(l:4::0::0:)STRING}.$FORMAT            >& /dev/null
 
    echo
done
 
 
for STRING in $STRINGS[*]; do
 
    echo -n "Creating graphics for string ${(l:4::0::0:)STRING}"
 
    HISTOGRAMS=()
 
    for INPUT_FILE in $INPUT_FILES[*]; do
 
        set_variable  TOKEN         ${INPUT_FILE%%/*}
        set_variable  STRING_TXT    ${TMPDIR:-/tmp}/H0\[${TOKEN}\].txt
        set_variable  STRING_ROOT   ${TMPDIR:-/tmp}/H0\[${TOKEN}\].root
 
        rm -f $STRING_TXT $STRING_ROOT
 
        set_variable  FIRST_FLOOR  1
 
        for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
            sput $STRING $FLOOR
            sget  MODULE
 
            set_variable P1 `JPrintResult -f "${INPUT_FILE}:H0\[${MODULE}\]" -F "GetMean"`
            set_variable P2  0.0
            
            printf "%2d %7.3f %7.4f\n" $FLOOR $P1 $P2  >>  $STRING_TXT
 
        done
 
        JGraph -f $STRING_TXT -o $STRING_ROOT
 
        HISTOGRAMS+=("${STRING_ROOT}:.*")
 
    done
 
    JPlot1D \
        -f "$HISTOGRAMS[*]"                          \
        -x "0.5 18.5"                                \
        -y "6.0 8.0"                                 \
        -> "floor"                                  \
        -\^ "rate [Hz]"                              \
        -T "${(l:4::0::0:)STRING}"                   \
        -o H0-${(l:4::0::0:)STRING}.$FORMAT          $BATCH
    
    echo
 
done
 
detach