plot-acoustics-fit.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

zmodload zsh/mathfunc

set_variable   DEBUG       2
set_variable   WORKDIR     ${TMPDIR:-/tmp}/
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

if do_usage $*; then
    usage "$script (input file)+"\
          "\nInput files correspond to stage files of a pre-calibration procedure."
fi

if (( $# == 0 )); then
    fatal "Wrong number of arguments."
fi

set_array     INPUT_FILES         $argv[1,-1]

let "YMIN = 0.0"
let "YMAX = 0.0"

OUTPUT_FILES=()

for (( i = 1; i <= ${#INPUT_FILES}; ++i )); do

    set_variable  INPUT_FILE  "$INPUT_FILES[${i}]"
    set_variable  TITLE       "${${INPUT_FILE##*/stage-}%%.log}"

    if ( `ls $INPUT_FILE >& /dev/null` ); then
        echo "Processing $INPUT_FILE"
    else
        fatal "Missing input file $INPUT_FILE"
    fi
    
    #awk '/chi2\[(1|3)\]/ { if ($3 < 1.0e10) { printf "%5d %9.5f\n", NR, $3 }}' $INPUT_FILE  >  $WORKDIR/g_${i}.txt
    awk '/chi2\[.\]/ { if ($3 < 1.0e10) { printf "%5d %9.5f\n", NR, $3 }}' $INPUT_FILE  >  $WORKDIR/g_${i}.txt

    sort -g -k 2 $WORKDIR/g_${i}.txt | tail -1 | read N Y

    if (( $Y > $YMAX )); then
        let "YMAX = $Y"
    fi
    
    JGraph \
        -f $WORKDIR/g_${i}.txt                   \
        -o $WORKDIR/g_${i}.root                  \
        -T $TITLE

    OUTPUT_FILES+=$WORKDIR/g_${i}.root
done

let "YMAX = floor(2.0 * ($YMAX + 0.1 * ($YMAX - $YMIN)) + 1.0) * 0.5"

JPlot1D \
    -f${^OUTPUT_FILES}":.*"                  \
    -\> "step"                               \
    -\^ "#chi^{2}/NDF"                       \
    -y "$YMIN $YMAX"                         \
    -N "X 505"                               \
    -S  0.5                                  \
    -T ""                                    \
    -L TR                                    \
    -o pre-calibration.$FORMAT               $BATCH

rm -f $WORKDIR/g_*.root