plot-Sirene.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}

# ------------------------------------------------------------------------------------------
#
#                         Utility script to plot JSirene results.
#
# ------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR

zmodload zsh/mathfunc

set_variable:  DEBUG       SIRENE_DEBUG          2
set_variable:  WORKDIR     SIRENE_WORKDIR        ./
set_variable:  FORMAT      GRAPHICS_FORMAT       gif
set_variable+  BATCH       GRAPHICS_BATCH        -B

if do_usage $*; then
    usage "$script (input file)+"
fi

if (( $# == 0 )); then
    fatal "No input files."
fi

INPUT_FILES=($argv[*])

typeset -A RADIATION

RADIATION[EErad]="10 19"
RADIATION[Brems]="20 29"
RADIATION[GNrad]="30 39"
RADIATION[DIS]="100 101"

for KEY VALUE in ${(@kv)RADIATION}; do

    JPlot1D \
        -f${^INPUT_FILES}:rad             \
        -XX -PX                           \
        -y "$VALUE" -Y                    \
        -\> "E [GeV]"                     \
        -\^ "Es [GeV]"                    \
        -T "$KEY"                         \
        -O "]["                           \
        -o $WORKDIR/rad_${KEY}.$FORMAT    $BATCH
        
done

RADIATION[MCS]="0 1"

for KEY VALUE in ${(@kv)RADIATION}; do

    JPlot1D \
        -f${^INPUT_FILES}:rms             \
        -XX -PX                           \
        -y "$VALUE" -Y                    \
        -z "1e-9 1e-2"                    \
        -\> "E [GeV]"                     \
        -\^ "#theta^{2} [rad^{2}]"        \
        -T "$KEY"                         \
        -O "]["                           \
        -o $WORKDIR/rad_${KEY}.$FORMAT    --! $BATCH
        
done

JPlot1D \
    -f${^INPUT_FILES}:cpu             \
    -XX -Y                            \
    -\> "E [GeV]"                     \
    -\^ "cpu [ms]"                    \
    -T ""                             \
    -o $WORKDIR/cpu.$FORMAT           $BATCH