acoustics-fit_00000075.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA006 (75) 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            75
#set_array     RUNS                   9345-9356
#set_variable: MUL                    ACOUSTICS_MUL      -1.17e-2
set_array     RUNS                   9700 9777 9884 9984 10000 10113 10200 10263
set_variable: MUL                    ACOUSTICS_MUL      -1.21e-2
set_variable: STAGES                 ACOUSTICS_STAGES   A 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       \
        -t A03166911          \
        -p A02980895          \
        -q A02980908          \
        -WW                   \
        -o $DETECTOR_INITIAL
fi

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

    cat>$TRIPOD_INITIAL<<EOF
# See https://git.km3net.de/calibration/input_tables for the origin(s) of the values
2  +587441.700 +4017014.700 -3459.600
3  +587848.700 +4016749.700 -3459.600
4  +587657.170 +4017113.940 -3459.600
5  +587906.500 +4017094.300 -3459.600
EOF
fi

# Apply modifications to the detector, if any.

JEditDetector \
    -a $DETECTOR_INITIAL     \
    -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