plot-katoomba.sh

Go to the documentation of this file.

#!/usr/bin/env 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:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -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`
 
JAcoustics.sh   $DETECTOR_ID
 
source JAcousticsToolkit.sh
 
CHECK_EXIT_CODE
 
 
JPlot1D \
    -f${^INPUT_FILES}:chi2                       \
    -x "0 10.0"                                  \
    -> "#chi^{2}/NDF"                           \
    -\^ "number of events [a.u.]"                \
    -T ""                                        \
    -o chi2.$FORMAT                              $BATCH
 
JPlot1D \
    -f${^INPUT_FILES}:hn                         \
    -XX                                          \
    -> "number of iterations"                   \
    -\^ "number of events [a.u.]"                \
    -T ""                                        \
    -o step.$FORMAT                              $BATCH
 
 
typeset -A EMITTERS
 
get_tripods+     $WORKDIR/tripod.txt       EMITTERS
get_transmitters+ $WORKDIR/transmitter.txt EMITTERS
 
for EMITTER in ${(@k)EMITTERS}; do
 
    for (( i = 1; $i <= ${#INPUT_FILES}; ++i )); do
    
        JOpera2D \
            -f "${INPUT_FILES[$i]}:ha\[${EMITTER}\]"     \
            -f "${INPUT_FILES[$i]}:hb\[${EMITTER}\]"     \
            -u Subtract                                  \
            -o $WORKDIR/subtract.root                    >& /dev/null
 
        JOpera2D \
            -f "subtract.root:Subtract"                  \
            -f "${INPUT_FILES[$i]}:ha\[${EMITTER}\]"     \
            -u Divide                                    \
            -o $WORKDIR/ratio.root                       >& /dev/null
        
        JPlot2D \
            -f $WORKDIR/ratio.root:\.\*                  \
            -> "string"                                 \
            -< "floor"                                  \
            -\^ "inefficiency"                           \
            -z "1.0e-4 1.05" -Z                          \
            -O COLZ                                      \
            -T "$EMITTER"                                \
            -o ratio\[${EMITTER}\]_${i}.$FORMAT          $BATCH       
        
        rm -f $WORKDIR/subtract.root
        rm -f $WORKDIR/ratio.root
 
    done
done