detector-XY:fit1d.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+  BACKGROUND  ACOUSTICS_BACKGROUND:  -B
set_variable   FORMULA     "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3]))"  

if do_usage $*; then
    usage "$script <string identifier>"
fi

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

set_variable  STRING              $argv[1]

set_variable  DETECTORXY_TXT      $WORKDIR/detectorxy_${STRING}.txt

sort -gr -k 4,4 $DETECTORXY_TXT | tail -1 | read X Y CHI2 RMS

printf "optimum %7.3f %7.3f %7.3f %7.3f\n" $X $Y $CHI2 $RMS

awk -v Y=$Y '{ if ( $2 == Y )  { printf "%7.3f %7.3f %7.3f\n", $1, $3, $4 } }' $DETECTORXY_TXT  >!  $WORKDIR/detectorx_${STRING}.txt
awk -v X=$X '{ if ( $1 == X )  { printf "%7.3f %7.3f %7.3f\n", $2, $3, $4 } }' $DETECTORXY_TXT  >!  $WORKDIR/detectory_${STRING}.txt

for var in x y; do

    set_variable  DETECTOR_TXT        $WORKDIR/detector${var}_${STRING}.txt
    set_variable  DETECTOR_ROOT       $WORKDIR/detector${var}_${STRING}.root

    JGraph \
        -f $DETECTOR_TXT                          \
        -o $DETECTOR_ROOT

    sort -gr -k 2,2 $DETECTOR_TXT | tail -1 | read X Z RMS

    let XMIN="$X - 0.5"
    let XMAX="$X + 0.5"

    JFit \
        -f ${DETECTOR_ROOT}:\.\*                  \
        -o $WORKDIR/fit.root                      \
        -F "$FORMULA"                             \
        -x "$XMIN $XMAX"                          \
        -@ "p0 =  $Z"                             \
        -@ "p1 =  2.0"                            \
        -@ "p2 =  $X"                             \
        -@ "p3 =  0.5"                            \
        -d $DEBUG

    mv $WORKDIR/fit.root $DETECTOR_ROOT

    JPlot1D \
        -f ${DETECTOR_ROOT}:\.\*                              \
        -\> "#Delta${var} [m]"                                \
        -\^ "#chi^{2}/NDF"                                    \
        -N "X 505"                                            \
        -T "[${(l:4::0::0:)STRING}].${var}"                   \
        -o $WORKDIR/detector${var}_${(l:4::0::0:)STRING}.gif  $BACKGROUND

    rm -f  $DETECTOR_ROOT

done