acoustics-fit_00000094.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA009 (94) 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            94
set_array     RUNS                   10360 10361 10363 10370 10390 10391 10411 10412
set_variable: MUL                    ACOUSTICS_MUL      -1.2e-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 A03209429          \
        -p A03209525          \
        -q A03209511          \
        -W                    \
        -o $DETECTOR_INITIAL
fi

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

    cat>$TRIPOD_INITIAL<<EOF
# Very preliminary acoustic emitter positions from survey
3  +587848.700 +4016749.700 -3450.467
4  +587657.170 +4017113.940 -3459.600
5  +587906.500 +4017094.300 -3459.600
6  +586939.550 +4017018.130 -3447.020
7  +587198.628 +4016228.693 -3433.306
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 RUN in $RUNS[*]; do
        cp -p $HOMEDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root  $WORKDIR
    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