detector-XY: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

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])/([4]*[4]) + 0.5*(y-[3])*(y-[3])/([5]*[5]))"

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

if (( $# == 0 || $# > 2 )); then
    fatal "Wrong number of arguments."
fi

set_variable  STRING              $argv[1]

set_variable  DETECTOR_TXT        $WORKDIR/detectorxy_${STRING}.txt
set_variable  DETECTOR_ROOT       $WORKDIR/detectorxy_${STRING}.root

if [[ ! -f $DETECTOR_TXT ]] then
    fatal "No detector data file $DETECTOR_TXT."
fi

JGraph2D \
    -f $DETECTOR_TXT                          \
    -o $DETECTOR_ROOT

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

let XMIN="$X - 0.5"
let XMAX="$X + 0.5"
let YMIN="$Y - 0.5"
let YMAX="$Y + 0.5"

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

set_array P2 `JPrintFit -f $WORKDIR/fit.root:user -@p2`
set_array P3 `JPrintFit -f $WORKDIR/fit.root:user -@p3`

printf "nominal %4d %7.3f %7.3f %7.3f\n" $STRING   0.0    0.0  `JPrintFit -f $WORKDIR/fit.root:user -x "  0.0    0.0 "`
printf "optimal %4d %7.3f %7.3f %7.3f\n" $STRING $P2[1] $P3[1] `JPrintFit -f $WORKDIR/fit.root:user -x "$P2[1] $P3[1]"`

if (( $# == 2 )); then

    set_variable  DETECTOR        $argv[2]

    JEditDetector -a $DETECTOR -S "$STRING add $P2[1] $P3[1] 0.0" -o $DETECTOR
fi

mv $WORKDIR/fit.root $DETECTOR_ROOT

typeset -Z 4 STRING

JPlot2D \
    -f ${DETECTOR_ROOT}:\.\*                  \
    -\> "#Deltax [m]"                         \
    -\< "#Deltay [m]"                         \
    -\^ "#chi^{2}/NDF"                        \
    -N "X 505"                                \
    -N "Y 505"                                \
    -O CONT3                                  \
    -T "[$STRING]"                            \
    -o $WORKDIR/detectorxy_${STRING}.gif      $BACKGROUND

rm -f  $DETECTOR_ROOT