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      /tmp/

if ( do_usage $* ); then
    usage "$script [working directory]"
fi

case $# in
    1) 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[*]

HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:pmt/g'`)

JPlot1D \
    -f "$HISTOGRAMS[*]"    \
    -\> "cos(#theta)"      \
    -\^ "acceptance"       \
    -y "0 1.6"             \
    -L TR                  \
    -T "PMT"               \
    -o pmt.gif

M=2

HISTOGRAMS=(`echo ${INPUT_FILES} | eval sed 's/\.root/\.root:M\\.${M}/g'`)

JPlot1D \
    -f "$HISTOGRAMS[*]"    \
    -\> "D [m]"            \
    -\^ "rate [Hz]"        \
    -x "0 5"               \
    -y "1 1e2" -Y          \
    -T "M = ${M}"          \
    -o k40.gif