pre-calibration_D0ORCA010.sh

Go to the documentation of this file.

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

# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D0ORCA010 (100) 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            100
set_array     RUNS                   11342 11345 11347 11374 11376
set_variable: TREADS                 ACOUSTICS_TREADS    5

if do_usage $*; then
    usage "$script <source directory> [user directory]"\
          "\nThe source directory corresponds to a pre-calibrated D_ORCA006 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


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]           \
        -WW                   \
        -o $DETECTOR_INITIAL
fi

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

    cat>$TRIPOD_INITIAL<<EOF
# Very preliminary acoustic emitter positions from survey
2  +256815.500 +4743395.000 -2436.677
3  +257096.200 +4743636.000 -2439.354
4  +256919.238 +4743354.178 -2439.354
EOF
fi

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

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

    # Apply modifications.

    JEditDetector \
        -a $DETECTOR_INITIAL     \
        -o detector.datx --!

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

    JEditDetector \
        -a detector.datx                 \
        -R "817295048 B1"                \
        -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[*]

cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.015;
numberOfHits  = 150;
EOF

cat>disable.txt<<EOF
# disable transmissions
# <emitter identifier> <module identifier>
# all emit. and ORCA.0003.16
# emitter 2 and ORCA.0009.12
# emitter 3 and ORCA.0018.06
# emitter 3 and ORCA.0001.12
# emitter 2 and ORCA.0011.13
# emitter 4 and ORCA.0011.02
-1 808966194
 2 808472260
 3 817301292
 3 809503416
 2 808978694
 4 808971812
EOF

cat>transmitter.txt<<EOF
# $script
EOF

cat>hydrophone.txt<<EOF
# Data from D0ORCA010_hydrophones.csv
   1 -1  -0.93   0.46   0.55
   2  0   0.02   1.04   0.55
   3  0   0.46   0.93   0.55
   9  0  -0.75  -0.72   0.55
  10 -1   0.11   1.04   0.55
  11 -1   0.54   0.89   0.55
  12  0   0.46   0.93   0.55
  17  0   0.61   0.84   0.55
  18  0   1.03   0.16   0.55
  19  0   0.11   1.04   0.55
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/${^ACOUSTICS_AUXS}.txt  $WORKDIR
    cp -p $HOMEDIR/${^INPUT_FILES}         $WORKDIR
fi

cat>script.txt<<EOF

fix  string  $FIXED_STRINGS[*]
fix  tripod  $FIXED_TRIPODS[*]

stage   0   1  0  100.0e-6  10.0  10  0.5       
stage  1B   1  0  100.0e-6  10.0   0  0.001  0.1
stage  2A   1  0  100.0e-6  10.0  20  0.2    0.5
stage  2C   1  0  100.0e-6  10.0  20  0.2       

initialise

stage  2a   1  0  100.0e-6  10.0  20  0.2    0.5
stage  2c   1  0  100.0e-6  10.0  20  0.2       
stage  3a   1  0   50.0e-6  10.0  20  0.2    0.5
stage  3c   1  0   50.0e-6  10.0  20  0.2       
EOF

JSydney.sh  $DETECTOR  $INPUT_FILES[*]  script.txt

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