JPrintK40Rates.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                     Utility script to evaluate rates due to K40 decays.
#
# ------------------------------------------------------------------------------------------


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   darkcount_Hz   700.0
set_variable   rateL0_Hz     4500.0

# https://elog.km3net.de/Analysis/597

set_array      rateL1_Hz     568  49.1  5.48  0.479  # 0.0613  0.000781  0.000133

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


function typeout()
{
    printf "(%4.0f, " $(( $rateL0_Hz + $darkcount_Hz ))
    printf "{ %6.1f" $rateL1_Hz[1]
    printf ", %6.2f" $rateL1_Hz[2,-1]
    printf "})\n"
}


# https://git.km3net.de/common/jpp/-/issues/351

set_variable  WORKDIR    ${TMPDIR:-/tmp}/.K40

mkdir -p     $WORKDIR

cp $DIR/JRateK40.cc         $WORKDIR/rateL0.cc 
cp $DIR/JMultiplicityK40.cc $WORKDIR/rateL1.cc 

cat>$WORKDIR/Makefile<<EOF
include \$(JPP_DIR)/make/JMakefile
EOF

for TARGET in KM3NeT KM3NeT_highQE; do
    if (( 1 )); then
    make -C $WORKDIR clean
    make -C $WORKDIR $TARGET

    $WORKDIR/rateL0 \
        -d 0         >& $WORKDIR/rateL0_${TARGET}.txt

    $WORKDIR/rateL1 \
        -n 100000000 \
        -F 4         \
        -G 2         \
        -D "0.216 5" \
        -A 1         \
        -U           \
        -o /dev/null \
        -d 0         >& $WORKDIR/rateL1_${TARGET}.txt
    fi
    eval set_variable  rateL0_${TARGET}_Hz  `awk '/Rate *1Dx1D/   { print $4 }' $WORKDIR/rateL0_${TARGET}.txt`
    eval set_array     rateL1_${TARGET}_Hz  `awk '/Rate\[.*\] *=/ { print $3 }' $WORKDIR/rateL1_${TARGET}.txt`

done

#rm -rf $WORKDIR


typeout

# correct rates

rateL0_Hz=$(( $rateL0_Hz * $rateL0_KM3NeT_highQE_Hz / $rateL0_KM3NeT_Hz ))

darkcount_Hz=$(( $darkcount_Hz * 3 ))

for (( i = 1; i <= ${#rateL1_Hz}; ++i )); do
    rateL1_Hz[$i]=$(( $rateL1_Hz[$i] * $rateL1_KM3NeT_highQE_Hz[$i] / $rateL1_KM3NeT_Hz[$i] ))
done

typeout