module-Z: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])/([3]*[3]))"
set_variable  MINIMAL_DZ          0.01

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

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

set_variable  STRING              $argv[1]
set_variable  FLOORS              $argv[2]

for (( FLOOR = ${FLOORS%%-*}; $FLOOR <= ${FLOORS##*-}; FLOOR += 1 )); do

    typeset -Z 0 STRING
    typeset -Z 0 FLOOR

    set_variable  MODULE_TXT      $WORKDIR/modulez_${STRING}_${FLOOR}.txt
    set_variable  MODULE_ROOT     $WORKDIR/modulez_${STRING}_${FLOOR}.root

    if [[ ! -f $MODULE_TXT ]] then
        fatal "No file $MODULE_TXT."
    fi

    sort -g  -k 2,2 $MODULE_TXT | tail -1 | read X ZMAX RMS
    sort -gr -k 2,2 $MODULE_TXT | tail -1 | read X ZMIN RMS

    printf "minimum %4d %2d %7.3f %7.3f\n" $STRING $FLOOR   $X   $ZMIN

    let DZ="0.5*($ZMAX - $ZMIN)"

    if (( $ZMAX - $ZMIN < $MINIMAL_DZ )); then
        set_variable OPTION "-= p1 = 0.0; p2 = 0.0; p3 = 1.0"
    else
        set_variable OPTION 
    fi

    JGraph \
        -f $MODULE_TXT                            \
        -o $MODULE_ROOT

    JFit \
        -f ${MODULE_ROOT}:\.\*                    \
        -o $WORKDIR/fit.root                      \
        -F "$FORMULA"                             \
        -@ "p0 =  $ZMIN"                          \
        -@ "p1 =  $DZ"                            \
        -@ "p2 =  $X"                             \
        -@ "p3 =  0.15"                           \
        "$OPTION"                                 \
        -d $DEBUG

    set_array    P2   `JPrintFit -f $WORKDIR/fit.root:\.\* -@p2`

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

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

        set_variable  DETECTOR        $argv[3]

        if (( $ZMAX - $ZMIN >= $MINIMAL_DZ )); then
            
            set_variable  MODULE          `getModule -a $DETECTOR -L "$STRING $FLOOR"`

            JEditDetector -a $DETECTOR -M "$MODULE add 0.0 0.0 $P2[1]" -o $DETECTOR

        else
            echo "No update of detector $DETECTOR; z-range ($ZMAX-$ZMIN) < $MINIMAL_DZ."
        fi
    fi

    mv $WORKDIR/fit.root $MODULE_ROOT

    typeset -Z 4 STRING
    typeset -Z 2 FLOOR

    JPlot1D \
        -f ${MODULE_ROOT}:\.\*                         \
        -\> "#Deltaz [m]"                              \
        -\^ "#chi^{2}/NDF"                             \
        -N "X 505"                                     \
        -T "($STRING,$FLOOR)"                          \
        -o $WORKDIR/modulez_${STRING}_${FLOOR}.gif     $BACKGROUND

    rm -f  $MODULE_ROOT

done

if [[ "$FLOORS" == "1-18" ]]; then

    echo -n "Creating graphics for string $STRING.."

    montage \
        -tile 6x3                           \
        -geometry +0+0                      \
	$WORKDIR/modulez_${STRING}_*.gif    \
        $WORKDIR/modulez_${STRING}.gif      >& /dev/null

    echo " -> $WORKDIR/modulez_${STRING}.gif"

    rm -f $WORKDIR/modulez_${STRING}_*.gif

fi

---

Generated by   1.8.5