JRateK40.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
 
set_variable:  DEBUG       PHYSICS_DEBUG          3
set_variable:  WORKDIR     PHYSICS_WORKDIR        ${TMPDIR:-/tmp}
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script"
fi
 
typeset -A FILTER
 
FILTER[350]=15
FILTER[370]=10
FILTER[400]=8
FILTER[420]=8
FILTER[450]=8
FILTER[470]=10
FILTER[500]=20
 
FORMULA="exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"
 
KEYS=($(echo ${(k)FILTER} | tr ' ' '\n' | sort -g))
 
for KEY in $KEYS; do
 
    echo "$KEY [nm]"
 
    $JPP_DIR/examples/JPhysics/JRateK40 \
        -o $WORKDIR/R\[$KEY\].root                           \
        -F "${${FORMULA//\[1\]/$KEY}//\[2\]/$FILTER[$KEY]}"  \
        -d $DEBUG
done
 
$JPP_DIR/examples/JPhysics/JRateK40 \
    -o $WORKDIR/R\[all\].root                          \
    -d $DEBUG
 
KEYS=(all ${(k)FILTER})
    
JPlot1D \
    -@ "titleOffset = 1.2"                             \
    -f$WORKDIR/R\[${^KEYS}\].root:h0                   \
    -x "280 600"                                       \
    -> "#lambda [nm]"                                 \
    -\^ "#partialR/#partial#lambda [kHz/nm]"           \
    -G X                                               \
    -G Y                                               \
    -T ""                                              \
    -L TR                                              \
    -o k40.$FORMAT                                     $BATCH
 
rm -f $WORKDIR/R\[${^KEYS}\].root