plot-Ballarat.sh

Go to the documentation of this file.

#!/usr/bin/env 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       DYNAMICS_DEBUG        2
set_variable:  WORKDIR     DYNAMICS_WORKDIR      ./
set_variable:  FLOORS      DYNAMICS_FLOORS       1-18
set_variable:  HISTOGRAM   DYNAMICS_HISTOGRAM    H
set_variable:  TMAX_S      DYNAMICS_TMAX_S       100
set_variable:  TIMESTAMP   DYNAMICS_TIMESTAMP    UTC
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 parameters."
fi
 
set_variable     DETECTOR                $argv[1]
set_variable     INPUT_FILE              $argv[2,-1]
 
eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
for STRING in $STRINGS[*]; do
 
    JEditDetector -a $DETECTOR -k $STRING -o ${TMPDIR:-/tmp}/detector.datx
    
 
    $JPP_DIR/examples/JDynamics/JBallarat \
        -a ${TMPDIR:-/tmp}/detector.datx   \
        -f "$INPUT_FILE[*]"                \
        -o ${TMPDIR:-/tmp}/ballarat.root   \
        -T $TMAX_S                         \
        -d $DEBUG
 
    MODULES=()
 
    echo "Creating graphics for string $STRING"
 
    for (( FLOOR = ${FLOORS%%-*}; $FLOOR <= ${FLOORS##*-}; FLOOR += 1 )); do
 
        eval MODULE=\${STRING_${(l:4::0::0:)STRING}\[$(($FLOOR))\]:--1}
 
        JPlot1D \
            -f ${TMPDIR:-/tmp}/ballarat.root:"U\["${MODULE}"\].twist"      \
            -Y                                                             \
            -> "#Delta#phi [deg]"                                         \
            -s 1110                                                        \
            -T "${STRING}.${FLOOR}"                                        \
            -o ${TMPDIR:-/tmp}/U${STRING}_${FLOOR}.$FORMAT                 $BATCH
 
        MODULES+=($MODULE)
    done
 
 
    montage \
        -tile 6x3                                                          \
        -geometry +0+0                                                     \
        ${TMPDIR:-/tmp}/U${STRING}_{01..18}.$FORMAT                        \
        $WORKDIR/U${STRING}.$FORMAT                                        >& /dev/null
 
    JPlot1D \
        -w 1200x600                                                        \
        -f${TMPDIR:-/tmp}/ballarat.root:"${HISTOGRAM}\["${^MODULES}"\]"    \
        -y "-3.15 +3.15"                                                   \
        -\^ "twist [rad]"                                                  \
        -N "X 505"                                                         \
        -t "$TIMESTAMP"                                                    \
        -S 0.3                                                             \
        -T "${STRING}"                                                     \
        -O "]["                                                            \
        -o $WORKDIR/${STRING}.$FORMAT                                      $BATCH
done
 
montage \
    -tile 3x40 \
    -geometry +0+0 \
    $WORKDIR/0*.$FORMAT \
    $WORKDIR/twist.$FORMAT >& /dev/null