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