plot-katoomba.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

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

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable:  DEBUG       ACOUSTICS_DEBUG        2
set_variable:  WORKDIR     ACOUSTICS_WORKDIR      ./
set_variable+  BACKGROUND  ACOUSTICS_BACKGROUND:  -B

if do_usage $*; then
    usage "$script <detector file> (input file)+>"
fi

if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi

set_variable  DETECTOR            $argv[1]
set_array     INPUT_FILES         $argv[2,-1]

eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
eval `JPrintDetector -a $DETECTOR -O SUMMARY`

source JAcoustics.sh   $DETECTOR_ID

CHECK_EXIT_CODE


JPlot1D \
    -f "${INPUT_FILES//\.root/.root:chi2}"       \
    -x "0 10.0"                                  \
    -\> "#chi^{2}/NDF"                           \
    -\^ "number of events [a.u.]"                \
    -T ""                                        \
    -o chi2.gif                                  $BACKGROUND


typeset -A TRIPODS

get_tripods $WORKDIR/tripod.txt TRIPODS

for EMITTER in ${(@k)TRIPODS}; do

    for (( i = 1; $i <= ${#INPUT_FILES}; ++i )); do
    
        JOpera2D \
            -f "${INPUT_FILES[$i]}:hb\[${EMITTER}\]"     \
            -f "${INPUT_FILES[$i]}:ha\[${EMITTER}\]"     \
            -u Divide                                    \
            -o $WORKDIR/ratio.root
        
        JPlot2D \
            -f $WORKDIR/ratio.root:\.\*                  \
            -\> "string"                                 \
            -\< "floor"                                  \
            -\^ "efficiency"                             \
            -z "0.0 1.05"                                \
            -O COLZ                                      \
            -T "$EMITTER"                               \
            -o ratio\[${EMITTER}\]_${i}.gif              $BACKGROUND       
        
        rm -f $WORKDIR/ratio.root

    done
done