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        ${TMPDIR:-/tmp}/
set_variable   DETECTOR       $WORKDIR/km3net_reference.detx               
set_variable   NUMBER_OF_RUNS 10
set_variable:  FORMAT         GRAPHICS_FORMAT        gif
set_variable+  BATCH          GRAPHICS_BATCH         -B
set_array      ANGLES_DEG     0 1 2 3 4 5

if do_usage $*; then
    usage "$script <option>"\
          "\nPossible options: run, plot, clean."

fi

if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi

set_variable  OPTION              $argv[1]

if [[ $OPTION != "run" && $OPTION != "plot" && $OPTION != "clean" ]]; then
    fatal "Invalid option $OPTION."
fi

if [[ $OPTION == "run" ]]; then

    if [[ ! -f $DETECTOR ]] then
        JDetector.sh $DETECTOR
    fi
    
cat>$WORKDIR/trigger_parameters.txt<<EOF
trigger3DMuon.enabled        =  1;
trigger3DMuon.numberOfHits   =  5;
trigger3DMuon.gridAngle_deg  =  1;
ctMin           =  0.0;
TMaxLocal_ns    = 15.0;
EOF

    set_variable TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY
        
    INPUT_FILES=()

    for (( i = 1; $i <= $NUMBER_OF_RUNS; ++i)); do

        JSirene.sh             $DETECTOR  $JPP_DATA/genhen.km3net_wpd_V2_0.evt.gz  $WORKDIR/sirene_${i}.root
        JTriggerEfficiency.sh  $DETECTOR  $DETECTOR  $WORKDIR/sirene_${i}.root  $WORKDIR/trigger_efficiency_${i}.root  $WORKDIR/trigger_parameters.txt  $JPP_DATA/PMT_parameters.txt

        INPUT_FILES+=($WORKDIR/trigger_efficiency_${i}.root)
    done

    for ANGLE_DEG in $ANGLES_DEG[*]; do

        set_variable SIGMA_NS     3.0
        set_variable OUTLIERS     3
    
        set_variable OUTPUT_FILE  $WORKDIR/matrix\[${ANGLE_DEG}\ deg\].root

        $JPP_DIR/examples/JReconstruction//JMatrixNZ \
            -a $DETECTOR            \
            -f"$INPUT_FILES[*]"     \
            -o $OUTPUT_FILE         \
            -S ${SIGMA_NS}          \
            -A ${ANGLE_DEG}         \
            -O ${OUTLIERS}          \
            -d ${DEBUG} --!     
#           -U                      \
    done
fi

if [[ $OPTION == "plot" ]]; then

    if (( 0 )); then

        for H1 in h0 h1; do

            JPlot1D \
                -f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"      \
                -y "1 2e3" -Y                                             \
                -L TR                                                     \
                -T ""                                                     \
                -\^ "number of events [a.u.]"                             \
                -\> "#chi^{2}/NDF"                                        \
                -o chi2-${H1}.$FORMAT                                     $BATCH
        done
    fi

    if (( 1 )); then

        for H1 in p0 p1; do

            JPlot1D \
                -f"$WORKDIR/matrix["${^ANGLES_DEG}" deg].root:${H1}"      \
                -y "1e0 1e5" -Y                                           \
                -L TR                                                     \
                -T ""                                                     \
                -\^ "number of events [a.u.]"                             \
                -\> "probability"                                         \
                -o prob-${H1}.$FORMAT --!                                   $BATCH
        done
    fi
fi

if [[ $OPTION == "clean" ]]; then
    rm -f $WORKDIR/sirene_*.root
    rm -f $WORKDIR/trigger_efficiency_*.root
    rm -f $WORKDIR/matrix\[*\].root
fi