acoustics-fit_00000049.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D_ORCA006 (49) detector.
#
# --------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable  DETECTOR_ID            49
set_array     RUNS                   7600 7700 7800 7900 8000 8100 8200
set_variable: MUL                    ACOUSTICS_MUL      +0.0e-2
set_variable: STAGES                 ACOUSTICS_STAGES   1 2 3

if do_usage $*; then
    usage "$script [user directory]"\
          "\nThe user directory should be specified when this script is submitted in a batch queue."
fi

set_variable      WORKDIR            `pwd`

if (( $# == 1 )); then
    set_variable  HOMEDIR            $1
else
    set_variable  HOMEDIR            $WORKDIR
fi

set_variable  TRIPOD_INITIAL         tripod_initial.txt
set_variable  TRIPOD                 tripod.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx
set_variable  DETECTOR               detector.datx


# Run in HOMEDIR

mkdir -p $HOMEDIR

cd  $HOMEDIR

if [[ ! -f $DETECTOR_INITIAL ]]; then

    JDetectorDB \
        -D $DETECTOR_ID       \
        -r $RUNS[1]           \
        -WW                   \
        -o $DETECTOR_INITIAL
fi

if [[ ! -f $TRIPOD_INITIAL ]]; then

    cat>$TRIPOD_INITIAL<<EOF
# Data are from V. Bertin, private communications.
1  256877.5  4743716.7  -2438.4
2  256815.5  4743395.0  -2435.5
3  257096.2  4743636.0  -2439.5
EOF
fi

# Apply modifications to the detector, if any.

JEditDetector \
    -a $DETECTOR_INITIAL     \
    -@ "z = -2441.5"         \
    -s "-1 mul $MUL"         \
    -o $DETECTOR

cp -p  $TRIPOD_INITIAL  $TRIPOD

JAcoustics.sh $DETECTOR_ID

JAcousticsEventBuilder.sh  $DETECTOR  $RUNS[*]


# Run in WORKDIR

cd  $WORKDIR

if [ ! $HOMEDIR -ef $WORKDIR ]; then

    cp -p $HOMEDIR/$TRIPOD        $WORKDIR
    cp -p $HOMEDIR/$DETECTOR      $WORKDIR

    for RANGE in $RUNS[*]; do
        for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
            cp -p $HOMEDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root  $WORKDIR
        done
    done

    JAcoustics.sh $DETECTOR_ID

fi

for STAGE in `echo $STAGES`; do

    $JPP_DIR/examples/JAcoustics/acoustics-fit.sh  $DETECTOR  $TRIPOD  $STAGE  $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*event.root  >  $HOMEDIR/stage-${STAGE}.log

    if [ ! $HOMEDIR -ef $WORKDIR ]; then
        cp -vru * $HOMEDIR
    fi
done