plot-Sirene.sh

Go to the documentation of this file.

#!/usr/bin/env 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