pre-calibration_D1ARCA030.sh

Go to the documentation of this file.

#!/usr/bin/env zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the global-fit-of-global-fits procedure of the D1ARCA030 (267) detector.
#
# --------------------------------------------------------------------------------------------
#
# The start geometry is from D0ARCA030 version 1.1.3. The different
# files are expected to be in the $DIR directory, provided as command
# line argument to the script. Then can be recovered from the git
# archive with the following commands:
#
#   git clone git@git.km3net.de:auxiliary_data/calibration.git
#   getFile.sh calibration/ 232 1 detector G_DOM_1.1.3 detector
#   getFile.sh calibration/ 232 1 tripod D_1.1.3 tripod
 
 
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            267
set_array     RUNS                   24523-24538
set_variable: MUL                    ACOUSTICS_MUL          -1.21e-2
set_variable: THREADS                ACOUSTICS_THREADS       16
set_variable  SINGULARITY            -
set_variable  WORKDIR                `pwd`
set_variable  HOMEDIR                $WORKDIR
 
if do_usage $*; then
    usage "$script <source directory> [user directory]"\
          "\nThe source directory corresponds to a pre-calibrated D0ARCA030 detector."\
          "\nThe user directory should be specified when this script is submitted in a batch queue."
fi
 
case $# in
    2)  set_variable HOMEDIR              $2;&
    1)  set_variable DIR                  ${1:a};;
    *)  fatal "Wrong number of arguments."
esac
 
source JAcousticsToolkit.sh
 
# post-calibration of input detector
 
mkdir -p $WORKDIR/post-calibration
pushd    $WORKDIR/post-calibration
 
post-calibration_D0ARCA030.sh  $DIR
 
popd
 
set_variable  DETECTOR               $WORKDIR/post-calibration/detector.datx
set_variable  TRIPOD                 $WORKDIR/post-calibration/tripod.txt
set_variable  DETECTOR_INITIAL       detector_initial.datx
set_variable  TRIPOD_ABSOLUTE        tripod_absolute.txt
expand_array  RUNS
 
# Run in HOMEDIR
 
mkdir -p $HOMEDIR
cd       $HOMEDIR
 
cat>$TRIPOD_ABSOLUTE<<EOF
# TAB6 and TAB1, with there absolute coordinates
#  values corresponding to acoustic emitter position, 
#  see https://git.km3net.de/working_groups/calibration/-/issues/124
22   +587756.91    +4016036.98   -3429.73
23   +586738.80    +4018233.81   -3452.69
EOF
 
# Apply modifications.
# ID of D1ARCA030 assigned to the calibrated D0ARCA030 detector 
JEditDetector   -a $DETECTOR -@ "id=$DETECTOR_ID" -o detector.datx --!
 
# Copy the initial position of tripods
cp -p  $TRIPOD  tripod.txt
# Remove tripod 20
JEditTripod -f tripod.txt -o tripod.txt -r 20
 
# Shift the input geometry 10m up
SHIFT="0.0 0.0 10.0"
 
JEditTripod -f tripod.txt -o tripod.txt -T "-1 add $SHIFT"
JEditDetector -a detector.datx -o detector.datx -M "-1 add $SHIFT"
 
 
# Add the tripods with absolute positioning
typeset -A TRIPODS
get_tripods  $TRIPOD_ABSOLUTE  TRIPODS;  FIXED_TRIPODS=(${(k)TRIPODS[*]})
 
JEditTripod   -f tripod.txt  -r "${FIXED_TRIPODS}" -o tripod.txt
JMergeTripod  -f tripod.txt -f $TRIPOD_ABSOLUTE    -o tripod.txt
 
get_tripods  tripod.txt  TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
 
print "FIXED STRINGS:" $FIXED_STRINGS
print "FIXED TRIPODS:" $FIXED_TRIPODS
print "FINAL TRIPODS:" $FINAL_TRIPODS
 
 
cat>waveform.txt<<EOF
# Data from D1ARCA030_waveforms.csv
 28   9
-29   9
 25  10
-26  10
 29  13
-30  13
 34  14
-35  14
 49  16
-50  16
 23  18
-24  18
 39  19
-40  19
 47  21
-48  21
 41  22
-42  22
 43  23
-44  23
 40  24
-41  24
 44  25
-45  25
EOF
 
cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.030;
numberOfHits  = 400;
EOF
 
cat>disable.txt<<EOF
# list of bad transmissions
9 806649098
9 806657372
9 808432849
9 808454978
9 808474602
9 808957378
9 808958576
9 808981206
9 808981369
9 808987098
9 816919071
9 817305565
9 817595373
9 817605768
9 817608140
9 817802210
10 806649098
10 806657372
10 808430036
10 808437870
10 808468365
10 808474602
10 808956926
10 808964876
10 808974966
10 808977563
10 808979220
10 808984610
10 808984662
10 817305565
10 817569715
10 817584003
10 817606460
10 817608135
10 817616429
10 817802231
10 819737799
10 819852785
13 806649098
13 808430036
13 808430571
13 808447081
13 808454978
13 808474243
13 808474602
13 808480499
13 808952499
13 808956926
13 808964925
13 808974966
13 808981369
13 808984662
13 817305565
13 817802231
18 808447081
18 808961256
18 808976232
18 808976252
18 808985071
18 809503404
18 809521474
18 817305565
18 817605238
19 808956926
19 808971766
19 808976207
19 816929954
19 817305565
19 817569715
21 806476519
21 808454978
21 808474243
21 808961256
21 808964322
21 808971305
21 808979692
21 808984586
21 809538805
21 817300799
21 817305565
21 817320316
21 817565173
21 817581894
22 808440691
22 808454978
22 808474578
22 808944706
22 808956797
22 808979692
22 808980070
22 808981087
22 808981369
22 808981913
22 808982006
22 808986771
22 816929999
22 817300799
22 817305565
22 817319884
22 817597135
22 817605775
22 817608140
22 819762610
EOF
 
cat>transmitter.txt<<EOF
# Data from D1ARCA030_transmitters.csv
14    5  0   0.74   0.23   0.48
16    8  0   0.44   0.63   0.48
EOF
 
cat>hydrophone.txt<<EOF
# Data from D0ARCA028_hydrophones.csv
   9 -1  -0.57  -0.42   0.57
  10  0  -0.67   0.39   0.62
  11  0  -0.67   0.39   0.62
  13  0  -0.68  -0.37   0.62
  14  0  -0.20  -0.75   0.62
  12  0  -0.12   0.76   0.62
  15  0  -0.72  -0.27   0.62
  19  0  -0.16  -0.76   0.62
  16  0   0.07   0.77   0.62
  20  0   0.09  -0.77   0.62
  23  0  -0.59  -0.49   0.62
  24  0  -0.68  -0.37   0.62
  21  0  -0.46   0.62   0.62
  22  0  -0.26   0.73   0.62
  25  0  -0.74  -0.22   0.62
  30  0  -0.08  -0.77   0.62
  26  0  -0.06   0.77   0.62
  27  0  -0.21   0.74   0.62
  28  0   0.12   0.76   0.62
  32  0   0.47   0.62   0.62
   5  0  -0.76  -0.13   0.62
   1  0  -0.20  -0.75   0.62
   3  0   0.12  -0.76   0.62
   4  0   0.62  -0.46   0.62
   8  0  -0.52  -0.57   0.62
  14  0  -0.20  -0.75   0.62
  15  0  -0.75  -0.20   0.62
  17  0   0.08   0.77   0.62
  18  0  -0.22  -0.74   0.62
  31  0   0.71   0.31   0.62
EOF
 
JAcousticsEventBuilder.sh  detector.datx  $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/$ACOUSTICS_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  tripod  $FIXED_TRIPODS[*]
 
stage  1a  0                2  250.0e-6  10.0  20  1.0    1.0
stage  2a  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2    0.2
stage  2c  ${SINGULARITY}2  0  100.0e-6  10.0  20  0.2       
stage  3a  ${SINGULARITY}2  0   50.0e-6  10.0  10  0.15    0.15
stage  3c  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.3
 
# initialise
 
# stage  4a  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.2    0.2
# stage  4c  ${SINGULARITY}2  0   50.0e-6  10.0  20  0.2
EOF
 
JSydney.sh  detector.datx  $INPUT_FILES[*]  script.txt
 
#forcing the F0-F1 distance to its nominal value (private comm. with Edward).
set_variable DETECTOR_FINAL detector.datx 
eval `JPrintDetector -a $DETECTOR_FINAL -O SUMMARY`;
for STRING in ${STRINGS[*]}; do
   JEditDetector \
      -a $DETECTOR_FINAL       \
      -F "$STRING  64.72"   \
      -o $DETECTOR_FINAL    --!
done
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi