JMultiplicityK40.sh

Go to the documentation of this file.

#!/bin/zsh
#
# \author mdejong
#
version=1.0
script=${0##*/}
 
# ------------------------------------------------------------------------------------------
#
#                     Utility script to make coincidence rate plots.
#
# ------------------------------------------------------------------------------------------
 
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi
 
 
source $JPP_DIR/setenv.sh $JPP_DIR
 
 
set_variable   DIR         $JPP_DIR/examples/JPhysics/  
set_variable   DEBUG       2                            
set_variable   WORKDIR     ${TMPDIR:-/tmp}/                        
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script [working directory]"
fi
 
case $# in
    1) set_variable  WORKDIR  $1;;
esac
 
if [[ ! -f $WORKDIR/k40\[KM3NeT\].root ]]; then 
 
    $DIR/JMultiplicityK40 \
        -n 100000000 \
        -F 4         \
        -G 2         \
        -D "0.216 5" \
        -A 1         \
        -U           \
        -o $WORKDIR/k40\[KM3NeT\].root
fi
 
for (( a = 0.0; $a <= 1.0; a += 0.5 )); do
 
    set_variable LABEL `printf '%3.1f' $a`
 
    if [[ ! -f $WORKDIR/k40\[${LABEL}\].root ]]; then 
 
        $DIR/JMultiplicityK40 \
            -n 100000000 \
            -F 4         \
            -G 2         \
            -D "0.216 5" \
            -A 3 -a $a   \
            -U           \
            -o $WORKDIR/k40\[${LABEL}\].root
    fi
done
 
INPUT_FILES=(`ls $WORKDIR/k40\[*\].root`)
 
echo $INPUT_FILES[*]
 
JPlot1D \
    -f "${INPUT_FILES//\.root/.root:pmt}" \
    -> "cos(#theta)"      \
    -\^ "acceptance"       \
    -y "0 1.6"             \
    -L TR                  \
    -T "PMT"               \
    -o pmt.$FORMAT         $BATCH
 
M=2
 
JPlot1D \
    -f "${INPUT_FILES//\.root/.root:M.${M}}" \
    -> "D [m]"            \
    -\^ "rate [Hz]"        \
    -x "0 5"               \
    -y "1 1e2" -Y          \
    -L TR                  \
    -T "M = ${M}"          \
    -o k40.$FORMAT         $BATCH