detector-XY:fit1d.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      ./
set_variable   FORMULA     "[0] + [1] * log(1.0 + 0.5*(x-[2])*(x-[2])/([3]*[3]))"  
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
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}.$FORMAT  $BATCH
 
    rm -f  $DETECTOR_ROOT
 
done