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