pre-calibration_00000116.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA020 (116) 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            116
set_array     RUNS                   12617-12622
set_variable  MUL                    -1.21e-2

if do_usage $*; then
    usage "$script <detector file> <tripod file> [user directory]"\
          "\nThe user directory should be specified when this script is submitted in a batch queue."\
          "\nThe detector file and tripod file correspond to a pre-calibrated D0ARCA009 detector."
fi

set_variable      WORKDIR            `pwd`
set_variable      HOMEDIR            $WORKDIR

case $# in
    3)  set_variable HOMEDIR              $3;&
    2)  set_variable TRIPOD               $2;
        set_variable DETECTOR             $1;;
    *)  fatal "Wrong number of arguments."
esac

set_variable  TRIPOD_INITIAL         tripod_initial.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx


# Run in HOMEDIR

mkdir -p $HOMEDIR

cd  $HOMEDIR

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

    JDetectorDB \
        -D $DETECTOR_ID       \
        -r $RUNS[1]           \
        -V ""                 \
        -WW                   \
        -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
7   +587198.628  +4016228.693  -3433.306
8   +588257.896  +4017034.442  -3457.240
9   +587646.419  +4017109.956  -3456.980
10  +587763.722  +4017253.398  -3453.894
EOF
fi

# Apply modifications.

JEditDetector \
    -a $DETECTOR_INITIAL     \
    -s "-1 mul $MUL"         \
    -o detector.datx

eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  STRINGS_INITIAL=($STRINGS[*])
eval `JPrintDetector -a $DETECTOR         -O SUMMARY` 

JEditDetector   -a detector.datx  -r "$STRINGS[*]"          -o detector.datx
JMergeDetector  -a detector.datx  -a $DETECTOR              -o detector.datx
JEditDetector   -a detector.datx  -k "$STRINGS_INITIAL[*]"  -o detector.datx

FIXED_STRINGS=($STRINGS[*])

cp -p  $TRIPOD_INITIAL  tripod.txt

source JAcousticsToolkit.sh

typeset -A TRIPODS_INITIAL
typeset -A TRIPODS

get_tripods  $TRIPOD_INITIAL  TRIPODS_INITIAL
get_tripods  $TRIPOD          TRIPODS

JEditTripod   -f tripod.txt  -D "${(k)TRIPODS}"
JMergeTripod  -f tripod.txt  -f $TRIPOD  -o tripod.txt
JEditTripod   -f tripod.txt  -K "${(k)TRIPODS_INITIAL}"

FIXED_TRIPODS=(${(k)TRIPODS[*]})

set_variable  DETECTOR               detector.datx
set_variable  TRIPOD                 tripod.txt

JAcoustics.sh $DETECTOR_ID

cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 120;
EOF

JAcousticsEventBuilder.sh  $DETECTOR  $RUNS[*]

INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)


# Run in WORKDIR

cd  $WORKDIR

if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -p $HOMEDIR/$DETECTOR        $WORKDIR
    cp -p $HOMEDIR/$TRIPOD          $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}  $WORKDIR
    cp -p $HOMEDIR/{acoustics_fit_parameters,acoustics_trigger_parameters,disable,hydrophone,mechanics,sound_velocity,tripod,waveform}.txt  $WORKDIR
fi

source $JPP_DIR/examples/JAcoustics/acoustics-fit-toolkit.sh

timer_start

initialise

fixStrings  $FIXED_STRINGS[*]
fixTripods  $FIXED_TRIPODS[*]

stage_1A >&  stage-1A.log; timer_lap

initialise

stage_2  >&  stage-2.log;  timer_lap
stage_3  >&  stage-3.log;  timer_lap

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

clean