pre-calibration_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

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
expand_array  RUNS


# 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

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

    # Apply modifications.

    JEditDetector \
        -a $DETECTOR_INITIAL     \
        -s "-1 addz -6.9"        \
        -o $DETECTOR --!

    eval `JPrintDetector -a $DETECTOR -O SUMMARY`

    for STRING in ${STRINGS[*]}; do

        set_variable  MODULE  `getModule -a $DETECTOR -L "$STRING 0"`

        JEditDetector \
            -a $DETECTOR             \
            -M "$MODULE setz -2.9"   \
            -o $DETECTOR --!
    done
fi

if [[ ! -f $TRIPOD ]]; then
    cp -p  $TRIPOD_INITIAL  $TRIPOD
fi

JAcoustics.sh $DETECTOR_ID

cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.020;
quantile      = 0.9;
numberOfHits  = 90;
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

JConvertDetectorFormat -a $DETECTOR -o $HOMEDIR/detector_backup.datx
JDetachPMTs -a $DETECTOR -o $DETECTOR

timer_start

initialise
    
stage_1B  0.002 0.1   0  >&  stage-1B.log;  backup $HOMEDIR/1B;  timer_lap
stage_2a  0.2   0.5  20  >&  stage-2a.log;  backup $HOMEDIR/2a;  timer_lap
stage_2b  0.2        20  >&  stage-2b.log;  backup $HOMEDIR/2b;  timer_lap
stage_3a  0.2   0.5  20  >&  stage-3a.log;  backup $HOMEDIR/3a;  timer_lap
stage_3b  0.2        20  >&  stage-3b.log;  backup $HOMEDIR/3b;  timer_lap

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

clean

JAttachPMTs -a $DETECTOR -f $HOMEDIR/detector_backup.datx -o $DETECTOR

rm -f $HOMEDIR/detector_backup.datx