acoustics-Zmul.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      `pwd`
set_array      MULS        -2.0e-2 -1.7e-2 -1.6e-2 -1.5e-2 -1.4e-2 -1.3e-2 -1.2e-2 -1.1e-2 -1.0e-2 -0.9e-2 -0.9e-2 -0.8e-2 -0.5e-2 +0.0e-2 +0.5e-2
set_variable   FORMULA     "[0] + [1] * ((x-[2])*(x-[2])/([3]*[3]))"
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <detector specific pre-calibration script> <option>"\
          "\nAuxiliary script to make scan of pre-stretching of detector strings (see JEditDetector)."\
          "\nPossible options: run, plot or fit."
fi
 
if (( $# < 1 || $# > 2)); then
    fatal "Wrong number of arguments."
fi
 
# global variables
 
set_variable  SCRIPT              $argv[1]
set_variable  OPTION              $argv[-1]
 
if [[ $OPTION != "run" && $OPTION != "plot" && $OPTION != "fit" ]]; then
    fatal "Invalid option $OPTION."
fi
 
if [[ $OPTION == "run" ]]; then
 
    for MUL in $MULS[*]; do
 
        set_variable  DIR                $WORKDIR/mul=$MUL
 
        mkdir -p $DIR
        cd       $DIR
 
        # set varibales in script
 
        set_variable  ACOUSTICS_MUL      $MUL
 
        nohup $SCRIPT >& all.log &
    done
fi
 
if [[ $OPTION == "plot" || $OPTION == "fit" ]]; then
 
    set_variable  ACOUSTICS_ZMUL_TXT    ${TMPDIR:-/tmp}/acoustics-zmul.txt
    set_variable  ACOUSTICS_ZMUL_ROOT   ${TMPDIR:-/tmp}/acoustics-zmul.root
 
    rm -f $ACOUSTICS_ZMUL_TXT >& /dev/null
 
    for MUL in $MULS[*]; do
 
        set_variable  DIR                $WORKDIR/mul=$MUL
        set_variable  FILE               `ls $DIR/stage-(1|A).log`
 
        tail -1 $FILE | grep "Elapsed time" >& /dev/null
 
        if (( $? == 0 )); then
 
            set_array  RESULT  `grep -v "warning" $FILE | tail -2 | head -1`
 
            printf "%4.1e %6.3f\n" $MUL $RESULT[-1]  >>  $ACOUSTICS_ZMUL_TXT
        fi
    done
    
    JGraph \
        -f $ACOUSTICS_ZMUL_TXT                    \
        -o $ACOUSTICS_ZMUL_ROOT
 
    XMIN=`sort -k1g $ACOUSTICS_ZMUL_TXT | head -1 | awk '{ print $1 }'`
    XMAX=`sort -k1g $ACOUSTICS_ZMUL_TXT | tail -1 | awk '{ print $1 }'`
 
    let "XMIN = $XMIN - 1.0e-3"
    let "XMAX = $XMAX + 1.0e-3"
 
    if [[ $OPTION == "fit" ]]; then
 
        sort -gr -k2gr $ACOUSTICS_ZMUL_TXT | tail -1 | read X Y
 
        JFit \
            -f ${ACOUSTICS_ZMUL_ROOT}:\.\*            \
            -o ${TMPDIR:-/tmp}/fit.root                          \
            -F "$FORMULA"                             \
            -@ "p0 = $Y"                              \
            -@ "p1 = 0.1"                             \
            -@ "p2 = $X"                              \
            -@ "p3 = 0.002"                           \
            -d $DEBUG
 
        set_array P2 `JPrintFit -f ${TMPDIR:-/tmp}/fit.root:\.\* -@p2`
 
        printf "Optimal stretching factor %8.2e\n" $P2[1]
 
        mv ${TMPDIR:-/tmp}/fit.root $ACOUSTICS_ZMUL_ROOT
 
    fi
        
    JPlot1D \
        -f $ACOUSTICS_ZMUL_ROOT:\.\*        \
        -x "$XMIN $XMAX"                    \
        -> "M"                             \
        -\^ "#chi^{2}/NDF"                  \
        -T ""                               \
        -o acoustics-zmul.$FORMAT           $BATCH                   
fi