pre-calibration_00000133.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ARCA021 (133) 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            133
set_array     RUNS                   13290-13297
set_variable  MUL                    -1.21e-2

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

set_variable      WORKDIR            `pwd`
set_variable      HOMEDIR            $WORKDIR

case $# in
    2)  set_variable HOMEDIR              $2;&
    1)  set_variable DIR                  ${1:a};;
    *)  fatal "Wrong number of arguments."
esac


# post-calibration of input detector

mkdir -p $WORKDIR/post-calibration

pushd    $WORKDIR/post-calibration

$JPP_DIR/examples/JAcoustics/post-calibration_00000094.sh  $DIR

popd

set_variable  DIR                    $WORKDIR/post-calibration


set_variable  DETECTOR               $DIR/detector.datx
set_variable  TRIPOD                 $DIR/tripod.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx
set_variable  TRIPOD_INITIAL         tripod_initial.txt
expand_array  RUNS


# 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
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
12  +587600.000  +4017500.000  -3457.240
13  +587510.740  +4016869.160  -3451.700
EOF
fi

JEditDetector \
    -a $DETECTOR           \
    -r "19"                \
    -o detector_tmp.datx

eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`;  FINAL_STRINGS=($STRINGS[*])
eval `JPrintDetector -a detector_tmp.datx -O SUMMARY`;  FIXED_STRINGS=($STRINGS[*])

if [[ ! -f detector.datx ]]; then

    # Apply modifications.

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

    eval `JPrintDetector -a detector.datx -O SUMMARY`

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

        set_variable  MODULE  `getModule -a detector.datx -L "$STRING 0"`

        JEditDetector \
            -a detector.datx         \
            -M "$MODULE setz -11.7"  \
            -o detector.datx --!
    done

    JEditDetector   -a detector.datx  -r "$FIXED_STRINGS[*]"  -o detector.datx
    JMergeDetector  -a detector.datx  -a detector_tmp.datx    -o detector.datx
    JEditDetector   -a detector.datx  -k "$FINAL_STRINGS[*]"  -o detector.datx
fi

source JAcousticsToolkit.sh

typeset -A TRIPODS

get_tripods  $TRIPOD_INITIAL  TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
get_tripods  $TRIPOD          TRIPODS;  FIXED_TRIPODS=(${(k)TRIPODS[*]})

if [[ ! -f tripod.txt ]]; then

    cp -p  $TRIPOD_INITIAL  tripod.txt

    JEditTripod   -f tripod.txt  -r "${FIXED_TRIPODS}" -o tripod.txt
    JMergeTripod  -f tripod.txt  -f $TRIPOD            -o tripod.txt
    JEditTripod   -f tripod.txt  -k "${FINAL_TRIPODS}" -o tripod.txt
fi

set_variable  DETECTOR               detector.datx
set_variable  TRIPOD                 tripod.txt

echo Fixed strings: $FIXED_STRINGS[*]
echo Fixed tripods: $FIXED_TRIPODS[*]

JAcoustics.sh $DETECTOR_ID

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

JAcousticsEventBuilder.sh  $DETECTOR  $RUNS[*]

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

# Run in WORKDIR

cd  $WORKDIR

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

if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -p $HOMEDIR/$DETECTOR               $WORKDIR
    cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}         $WORKDIR
fi

timer_start

initialise

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

stage_0  0  2  250.0e-6  10.0  0.5        20  >&  stage-0.log;   backup $HOMEDIR/0;   timer_lap
stage_a  0  2  250.0e-6  10.0  0.5   1.0  20  >&  stage-1A.log;  backup $HOMEDIR/1A;  timer_lap
stage_b  0  2  250.0e-6  10.0  0.002 0.2   0  >&  stage-1B.log;  backup $HOMEDIR/1B;  timer_lap
stage_c  2  0  250.0e-6  10.0  0.3        20  >&  stage-1C.log;  backup $HOMEDIR/1C;  timer_lap
stage_a  2  0  100.0e-6  10.0  0.2   0.4  20  >&  stage-2A.log;  backup $HOMEDIR/2A;  timer_lap
stage_c  2  0  100.0e-6  10.0  0.2        20  >&  stage-2C.log;  backup $HOMEDIR/2C;  timer_lap

initialise

stage_a  2  0  100.0e-6  10.0  0.2   0.4  20  >&  stage-2a.log;  backup $HOMEDIR/2a;  timer_lap
stage_c  2  0  100.0e-6  10.0  0.2        20  >&  stage-2c.log;  backup $HOMEDIR/2c;  timer_lap
stage_a  2  0   50.0e-6  10.0  0.2   0.4  20  >&  stage-3a.log;  backup $HOMEDIR/3a;  timer_lap
stage_c  2  0   50.0e-6  10.0  0.2        20  >&  stage-3c.log;  backup $HOMEDIR/3c;  timer_lap

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

clean