JNarrabri.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       DYNAMICS_DEBUG                2
set_variable:  WORKDIR     DYNAMICS_WORKDIR              ./
set_variable:  TMAX_S      DYNAMICS_TMAX_S               1000.0                        
set_variable:  YMIN        DYNAMICS_YMIN                 0.00
set_variable:  YMAX        DYNAMICS_YMAX                 0.10
set_variable   DIR         $JPP_DIR/examples/JDynamics/  
set_variable:  TIMESTAMP   GRAPHICS_TIMESTAMP     utc
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector file> (input file)+"\
          "\nInput files correspond to the output of JKatoomba[.sh]."
fi
 
if (( $# < 2 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  DETECTOR            $argv[1]
set_variable  INPUT_FILE          $argv[2,-1]
 
eval `JPrintDetector -a $DETECTOR -O SUMMARY`
 
getMechanics.sh  $DETECTOR
 
$DIR/JNarrabri \
    -a $DETECTOR                              \
    -f "$INPUT_FILE[*]"                       \
    -o $WORKDIR/narrabri.root                 \
    -T $TMAX_S                                \
    -d $DEBUG --!
 
 
typeset -Z 4 STRING
 
for STRING in $STRINGS[*]; do
 
    JPrintResult \
        -f  "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt"    \
        -F "GetStdDev(1)"                                    \
        -F "GetStdDev(2)"                                    \
        -F "GetEntries - Integral"                           \
        -F "GetEntries"                                      | read -A RESULT
 
        printf "%04d %4.2f %4.2f %5.0f / %8.0f\n" $STRING $RESULT[*]
 
done
 
echo "Creating graphics output in directory $WORKDIR"
 
JPlot1D \
    -f $WORKDIR/narrabri.root:time                    \
    -> "#Deltat [s]"                                 \
    -\^ "number of events [a.u.]"                     \
    -XX                                               \
    -L TL                                             \
    -T ""                                             \
    -o $WORKDIR/time.$FORMAT                          $BATCH
 
JPlot1D \
    -w 1200x600                                       \
    -f $WORKDIR/narrabri.root:G\.\*orientation        \
    -\^ "orientation [rad]"                           \
    -y "-3.20 +3.20"                                  \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L BR                                             \
    -T ""                                             \
    -o $WORKDIR/orientation.$FORMAT                   $BATCH
 
JPlot1D \
    -w 1200x600                                       \
    -f $WORKDIR/narrabri.root:G\.\*amplitude          \
    -\^ "amplitude [rad]"                             \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L TR                                             \
    -T ""                                             \
    -o $WORKDIR/amplitude.$FORMAT                     $BATCH
 
JPlot1D \
    -w 1200x600                                       \
    -f $WORKDIR/narrabri.root:H\.\*orientation        \
    -\^ "orientation [rad]"                           \
    -y "-3.20 +3.20"                                  \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L BR                                             \
    -T ""                                             \
    -o $WORKDIR/orientation-L.$FORMAT                 $BATCH
 
JPlot1D \
    -w 1200x600                                       \
    -f $WORKDIR/narrabri.root:H\.\*amplitude          \
    -\^ "amplitude [rad]"                             \
    -N "X 505"                                        \
    -t "$TIMESTAMP"                                   \
    -L TR                                             \
    -T ""                                             \
    -o $WORKDIR/amplitude-L.$FORMAT                   $BATCH
 
for STRING in $STRINGS[*]; do
 
    JPlot1D \
        -w 1200x600                                               \
        -f  "$WORKDIR/narrabri.root:G\[${STRING}\]\.orientation"  \
        -f  "$WORKDIR/narrabri.root:H\[${STRING}\]\.orientation"  \
        -\^ "orientation [rad]"                                   \
        -y "-3.20 +3.20"                                          \
        -N "X 505"                                                \
        -t "$TIMESTAMP"                                           \
        -T "[${STRING}]"                                          \
        -o $WORKDIR/orientation_${STRING}.$FORMAT                 $BATCH
 
    JPlot1D \
        -w 1200x600                                               \
        -f  "$WORKDIR/narrabri.root:G\[${STRING}\]\.amplitude"    \
        -f  "$WORKDIR/narrabri.root:H\[${STRING}\]\.amplitude"    \
        -\^ "amplitude [rad]"                                     \
        -y "$YMIN $YMAX"                                          \
        -N "X 505"                                                \
        -t "$TIMESTAMP"                                           \
        -T "[${STRING}]"                                          \
        -o $WORKDIR/amplitude_${STRING}.$FORMAT                   $BATCH
 
    JPlot2D \
        -f  "$WORKDIR/narrabri.root:^\[${STRING}\]\.tilt"         \
        -> "dTx/dt [mrad/10 min]"                                \
        -< "dTy/dt [mrad/10 min]"                                \
        -O COLZ                                                   \
        -N "X 505"                                                \
        -N "Y 505"                                                \
        -s 101010                                                 \
        -T "[${STRING}]"                                          \
        -o $WORKDIR/tilt_${STRING}.$FORMAT                        $BATCH
done