pre-calibration_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  DEBUG                  2
set_variable  DETECTOR_ID            75
set_array     RUNS                   9700 9777 9884 9984 10000 10113 10200 10263
set_variable: MUL                    ACOUSTICS_MUL    -1.21e-2

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       \
        -@ "tcal = A03166911" \
        -@ "pcal = A02980895" \
        -@ "rcal = A02980908" \
        -V ""                 \
        -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

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

    # Apply modifications.

    JEditDetector \
        -a $DETECTOR_INITIAL     \
        -s "-1 mul $MUL"         \
        -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 -11.7"  \
            -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;
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/${^ACOUSTICS_KEYS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}         $WORKDIR
fi

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

timer_start

initialise

stage_a  0  2  250.0e-6  0.5   1.0  20  >&  stage-1a.log;  backup $HOMEDIR/1a;  timer_lap
stage_b  0  2  250.0e-6  0.002 0.2   0  >&  stage-1b.log;  backup $HOMEDIR/1b;  timer_lap
stage_c  0  2  250.0e-6  0.2        20  >&  stage-1c.log;  backup $HOMEDIR/1c;  timer_lap
stage_a  2  0  100.0e-6  0.2   0.4  20  >&  stage-2a.log;  backup $HOMEDIR/2a;  timer_lap
stage_c  2  0  100.0e-6  0.2        20  >&  stage-2c.log;  backup $HOMEDIR/2c;  timer_lap
stage_a  2  0   50.0e-6  0.2   0.4  20  >&  stage-3a.log;  backup $HOMEDIR/3a;  timer_lap
stage_c  2  0   50.0e-6  0.2        20  >&  stage-3c.log;  backup $HOMEDIR/3c;  timer_lap

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

clean