JMatrixNZ.sh

Go to the documentation of this file.

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

# ------------------------------------------------------------------------------------------
#
#                 Utility script to test JMatrixNZ on Monte Carlo true muons.
#
# ------------------------------------------------------------------------------------------


if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi


source $JPP_DIR/setenv.sh $JPP_DIR


set_variable  DEBUG          2
set_variable  WORKDIR        $JPP_DATA
set_variable  INPUT_FILE     $WORKDIR/trigger_efficiency+background.root 
set_variable  DETECTOR       $WORKDIR/km3net_reference.detx
set_variable  JPP            $JPP_DIR/examples/JFit/

if ( do_usage $* ); then
    usage "$script [detector file [input file]]"
fi

case $# in
    2) INPUT_FILE=$2;&
    1) DETECTOR=$1;;
esac


if (( 1 )); then

    print_variable     DETECTOR INPUT_FILE
    check_input_file   $DETECTOR $INPUT_FILE

    for ALPHA_DEG in 0 1 2 3 4 5; do

   set_variable SIGMA_NS    3.0
   set_variable OUTLIERS    3

   set_variable OUTPUT_FILE matrix\[${ALPHA_DEG}\ deg\].root
   
   $JPP/JMatrixNZ \
       -a $DETECTOR            \
       -f $INPUT_FILE          \
       -o $OUTPUT_FILE         \
       -S ${SIGMA_NS}          \
       -A ${ALPHA_DEG}         \
       -O ${OUTLIERS}          \
       -d ${DEBUG} --!     
#      -U                      \

    done

fi


if (( 0 )); then

    for HISTOGRAM in h0 h1; do

   JPlot1D  \
       -f matrix\[0\ deg\].root:${HISTOGRAM}    \
       -f matrix\[1\ deg\].root:${HISTOGRAM}    \
       -f matrix\[2\ deg\].root:${HISTOGRAM}    \
       -f matrix\[3\ deg\].root:${HISTOGRAM}    \
       -f matrix\[4\ deg\].root:${HISTOGRAM}    \
       -f matrix\[5\ deg\].root:${HISTOGRAM}    \
       -y "1 2e3" -Y                            \
       -LL                                      \
       -\^ "number of events [a.u.]"            \
       -> "#chi^{2}/NDF"                       \
       -o chi2-${HISTOGRAM}.gif
       
    done

fi


if (( 1 )); then

    for HISTOGRAM in p0 p1; do

   JPlot1D  \
       -f matrix\[0\ deg\].root:${HISTOGRAM}    \
       -f matrix\[1\ deg\].root:${HISTOGRAM}    \
       -f matrix\[2\ deg\].root:${HISTOGRAM}    \
       -f matrix\[3\ deg\].root:${HISTOGRAM}    \
       -f matrix\[4\ deg\].root:${HISTOGRAM}    \
       -f matrix\[5\ deg\].root:${HISTOGRAM}    \
       -y "1e1 1e4" -Y                          \
       -LL                                      \
       -\^ "number of events [a.u.]"            \
       -> "probability"                        \
       -o prob-${HISTOGRAM}.gif
       
    done

fi