acoustics-Zmul.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      `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 acoustics-fit 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
        set_variable  ACOUSTICS_STAGES   A

        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