software/JAcoustics/plot-canberra.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:  XMIN        CANBERRA_XMIN         -1.0e-3
set_variable:  XMAX        CANBERRA_XMAX         +1.0e-3
set_variable:  ZMIN        CANBERRA_ZMIN         -0.5
set_variable:  ZMAX        CANBERRA_ZMAX         +0.5
set_variable+  DISABLE     CANBERRA_DISABLE      -! $WORKDIR/disable.txt
set_variable+  REVERT      CANBERRA_REVERT       -r
set_variable:  OPTION      CANBERRA_OPTION       -O HIST
set_variable:  FORMAT      GRAPHICS_FORMAT       gif
set_variable+  BATCH       GRAPHICS_BATCH        -B
 
if do_usage $*; then
    usage "$script <detector file> (input file)+"\
          "\nInput file corresponds to the output of JKatoomba[.sh]."
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
 
JCanberra.sh  $DETECTOR  $INPUT_FILES[*]  $WORKDIR/canberra\[%\].root
 
source JAcousticsToolkit.sh
 
CHECK_EXIT_CODE
 
typeset -A EMITTERS
 
get_tripods+       $WORKDIR/tripod.txt       EMITTERS
get_transmitters+  $WORKDIR/transmitter.txt  EMITTERS
 
 
typeset -Z 4 STRING
typeset -Z 2 FLOOR
 
for STRING in $STRINGS[*]; do
 
    echo -n "Creating graphics for string $STRING.."
 
    for (( FLOOR = $FIRST_FLOOR; $FLOOR <= $LAST_FLOOR; FLOOR += 1 )); do
 
        eval MODULE=\${STRING_${(l:4::0::0:)STRING}\[$(($FLOOR))\]:--1}
        
        INPUT_FILES=(`ls $WORKDIR/canberra\[*\].root`)
 
        JPlot1D \
            -f${^INPUT_FILES}:${MODULE}                         \
            -> "#Deltat [s]"                                   \
            -\^ "number of events [a.u.]"                       \
            -x "$XMIN $XMAX"                                    \
            $OPTION                                             \
            -N "X 505"                                          \
            -L TR                                               \
            -T "(${STRING},${FLOOR})"                           \
            -o ${TMPDIR:-/tmp}/toa_${STRING}_${FLOOR}.$FORMAT   $BATCH
    done
 
    montage \
        -tile 6x3                                               \
        -geometry +0+0                                          \
        ${TMPDIR:-/tmp}/toa_${STRING}_{01..18}.$FORMAT          \
        $WORKDIR/toa_${STRING}${REVERT}.$FORMAT                 >& /dev/null
    
    mv  ${TMPDIR:-/tmp}/toa_${STRING}_00.$FORMAT                \
        $WORKDIR/toa_${STRING}_00${REVERT}.$FORMAT
 
    echo " -> $WORKDIR/toa_${STRING}${REVERT}.$FORMAT"
 
    rm -f ${TMPDIR:-/tmp}/toa_${STRING}_*.$FORMAT
done
 
 
for EMITTER in ${(@k)EMITTERS}; do
 
    JPlot2D \
        -f $WORKDIR/canberra\[${EMITTER}\].root:mean      \
        -> "string"                                      \
        -< "floor"                                       \
        -\^ "#Deltat [ms]"                                \
        -z "$ZMIN $ZMAX"                                  \
        -N "Y 10"                                         \
        -T "[${EMITTER}]"                                 \
        -O COLZ                                           \
        -p 57                                             \
        -o $WORKDIR/mean_${EMITTER}${REVERT}.$FORMAT      $BATCH
 
    JPlot2D \
        -f $WORKDIR/canberra\[${EMITTER}\].root:sigma     \
        -> "string"                                      \
        -< "floor"                                       \
        -\^ "#sigma [ms]"                                 \
        -z "0.0 0.1"                                      \
        -N "Y 10"                                         \
        -T "[${EMITTER}]"                                 \
        -O COLZ                                           \
        -p 57                                             \
        -o $WORKDIR/sigma_${EMITTER}${REVERT}.$FORMAT     $BATCH
done
 
typeset -i NX
typeset -i NY
 
let "NY = sqrt(${#EMITTERS})"
let "NX = (${#EMITTERS} + $NY - 1)/ $NY"
 
for KEY in mean sigma; do
 
    montage \
        -tile ${NX}x${NY}                                 \
        -geometry +0+0                                    \
        $WORKDIR/${KEY}_*${REVERT}.$FORMAT                \
        $WORKDIR/${KEY}${REVERT}.$FORMAT                  >& /dev/null
done
 
                                         
rm -f $WORKDIR/canberra\[*\].root