hydrophone-t0: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               ${ACOUSTICS_DEBUG:-2}
set_variable  WORKDIR             ${ACOUSTICS_WORKDIR:-./}
set_variable  BACKGROUND          ${ACOUSTICS_BACKGROUND:+-B}

if ( do_usage $* ); then
    usage "$script [detector file]"
fi

if (( ${#} > 1 )); then
    fatal "Wrong number of arguments."
fi

set_variable  FORMULA             "[0] * exp(-0.5 * (x-[1])*(x-[1]) / ([2]*[2]))"

JFit \
    -f "$WORKDIR/hydrophone.root:\[%\]"          \
    -o "$WORKDIR/fit.root"                       \
    -F "$FORMULA"                                \
    -@"p0 = GetMaximum * 0.8"                    \
    -@"p1 = GetMean"                             \
    -@"p2 = GetRMS * 0.3"                        \
    -d $DEBUG

JPlot1D \
    -f "$WORKDIR/fit.root:\[%\]"                 \
    -\> "#Deltat [s]"                            \
    -\^ "number of events [a.u.]"                \
    -N "X 505"                                   \
    -T ""                                        \
    -d $DEBUG                                    \
    -o hydrophonet0.gif                          $BACKGROUND

if (( ${#} == 1 )); then

    set_variable  DETECTOR            $argv[1]

    eval `JPrintDetector -a $DETECTOR -O SUMMARY`

    set_array P1  `JPrintFit -f "$WORKDIR/fit.root:\[%\]" -@ p1`

    let TIME_OFFSET_NS="$P1[1] * 1.0e9"

    printf "Add time offset %12.6d to each hydrophone.\n" $TIME_OFFSET_NS
    
    for STRING in $STRINGS[*]; do

        set_variable  MODULE   `getModule -a $DETECTOR -L "$STRING 0"`

        JEditDetector \
            -a $DETECTOR                         \
            -o $DETECTOR                         \
            -M "$MODULE SUB $TIME_OFFSET_NS"     \
            -d $DEBUG
    done
fi