JParramatta.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:  ZMAX_M      ACOUSTICS_ZMAX_M       0.0
set_variable:  TIMESTAMP   GRAPHICS_TIMESTAMP     utc
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable:  SIZE        GRAPHICS_SIZE          0.5
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script (input file)+"\
          "\nInput files correspond to the output of JKatoomba[.sh]."
fi
 
if (( $# < 1 )); then
    fatal "Wrong number of arguments."
fi
 
set_variable  INPUT_FILE          $argv[1,-1]
 
JParramatta \
    -f "$INPUT_FILE[*]"                       \
    -o $WORKDIR/parramatta.root               \
    -Z $ZMAX_M                                \
    -d $DEBUG --!
 
echo "Creating graphics output in directory $WORKDIR"
 
JPlot1D \
    -w 1200x600                                               \
    -f $WORKDIR/parramatta.root:g1                            \
    -\^ "#chi^{2}/NDF"                                        \
    -N "X 505"                                                \
    -t "$TIMESTAMP"                                           \
    -S $SIZE                                                  \
    -T ""                                                     \
    -o $WORKDIR/chi2-time.$FORMAT                             $BATCH
 
JPlot1D \
    -w 1200x600                                               \
    -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].amplitude"  \
    -\^ "amplitude [mrad]"                                    \
    -N "X 505"                                                \
    -t "$TIMESTAMP"                                           \
    -S $SIZE                                                  \
    -T ""                                                     \
    -L TR                                                     \
    -o $WORKDIR/amplitude.$FORMAT                             $BATCH
 
JPlot1D \
    -w 1200x600                                               \
    -f "$WORKDIR/parramatta.root:G\[[0-9][0-9]*\].stretching" \
    -\^ "stretching [%]"                                      \
    -N "X 505"                                                \
    -t "$TIMESTAMP"                                           \
    -S $SIZE                                                  \
    -T ""                                                     \
    -L BR                                                     \
    -o $WORKDIR/stretching.$FORMAT                            $BATCH
 
JPlot1D \
    -f "$WORKDIR/parramatta.root:hx"                          \
    -f "$WORKDIR/parramatta.root:hy"                          \
    -\^ "#DeltaTx/y [mrad]"                                   \
    -> "string"                                              \
    -T ""                                                     \
    -o $WORKDIR/txy.$FORMAT                                   $BATCH
 
typeset -Z4  STRING
 
for H1 in `JPrintName -f "$WORKDIR/parramatta.root:^\[[0-9][0-9]*\]\.tiltdeviation"`; do
 
    set_variable  STRING  `echo $H1 | sed 's/.*\[\‍([0-9]*\‍)\].*/\1/'`
 
    JPlot2D \
        -f  "$WORKDIR/parramatta.root:^\[${STRING}\]\.tiltdeviation" \
        -< "Ty - <Ty> [mrad]"                                       \
        -> "Tx - <Tx> [mrad]"                                       \
        -O COLZ                                                      \
        -G X                                                         \
        -G Y                                                         \
        -N "X 505"                                                   \
        -N "Y 505"                                                   \
        -s 101110                                                    \
        -T "[${STRING}]"                                             \
        -o $WORKDIR/tiltdeviation_${STRING}.$FORMAT                  $BATCH
done
 
if (( 0 )); then
 
    set_variable FORMULA      "[0] + [1]*x + [2]*x*x"
 
    for H1 in `JPrintName -f "$WORKDIR/parramatta.root:stretching\[[0-9][0-9]*\]"`; do
 
        set_variable  STRING  `echo $H1 | sed 's/.*\[\‍([0-9]*\‍)\].*/\1/'`
 
        H1=${H1/\[/\\\[}
        H1=${H1/\]/\\\]}
 
        JFit \
            -f "$WORKDIR/parramatta.root:${H1}"                       \
            -o ${TMPDIR:-/tmp}/fit.root                               \
            -F "$FORMULA"                                             \
            -@ "p0 = 0.0"                                             \
            -@ "p1 = 1.0e-5"                                          \
            -@ "p2 = 1.0e-6"                                          \
            -d $DEBUG
 
        JPlot1D \
            -f ${TMPDIR:-/tmp}/fit.root:\.\*                          \
            -> "|T| [mrad]"                                          \
            -\^ "stretching [%]"                                      \
            -y "-0.25 +0.25"                                          \
            -N "X 505"                                                \
            -T "$STRING"                                              \
            -o $WORKDIR/profile_${STRING}.$FORMAT                     $BATCH    
    done
fi