back-propagation_D0ARCA021.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
# --------------------------------------------------------------------------------------------
#
# Steering script for the back-propagation procedure of the D0ARCA021 (133) detector from the
# D0ARCA028(160) 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                   13491-13498
set_variable: THREADS                ACOUSTICS_THREADS       8
set_variable: ARCHIVE                CALIBRATION_ARCHIVE     
set_variable  WORKDIR                `pwd`
set_variable  HOMEDIR                $WORKDIR
 
#===========NEEDS ADJUSTMENT===========
if do_usage $*; then
    usage "$script"\
          "\nRequires variable CALIBRATION_ARCHIVE to be set for the user."
fi
 
source JAcousticsToolkit.sh
 
set_variable  DETECTOR               detector.datx
set_variable  TRIPOD_INITIAL         tripod_initial.txt
set_variable  TRIPOD_ABSOLUTE        tripod_absolute.txtdo
set_variable  TRIPOD                 tripod.txt
expand_array  RUNS
 
# Run in HOMEDIR
mkdir -p $HOMEDIR
cd       $HOMEDIR
 
# retrieve D0ARCA021 input files from calibration archive
getFile.sh $ARCHIVE $DETECTOR_ID 13491 detector G_1.0.0 $WORKDIR/D0ARCA021_G_1.0.0
getFile.sh $ARCHIVE $DETECTOR_ID 13491 tripod $WORKDIR/tripod_initial
getFile.sh $ARCHIVE $DETECTOR_ID 13491 transmitter $WORKDIR/transmitter
getFile.sh $ARCHIVE $DETECTOR_ID 13491 hydrophone $WORKDIR/hydrophone
getFile.sh $ARCHIVE $DETECTOR_ID 13491 waveform $WORKDIR/waveform
 
# copy D0ARCA028 files from previous back-propagation step
cp -p /sps/km3net/users/lvoorend/calibration_workdir/backpropogate_D1ARCA030/reproduce_ARCA28/detector.datx $WORKDIR/D0ARCA028.datx
cp -p /sps/km3net/users/lvoorend/calibration_workdir/backpropogate_D1ARCA030/reproduce_ARCA28/tripod.txt $WORKDIR/tripod_D0ARCA028.txt
 
# ceate detector file for D0ARCA021 by copying D0ARCA028 and removing the modules that are not present in D0ARCA021
cp -p $WORKDIR/D0ARCA028.datx $WORKDIR/detector.datx
 
JEditDetector \
    -a $DETECTOR     \
    -r "1 3 4 8 14 15 17 18 31"         \
    -@ "id=$DETECTOR_ID"          \
    -o $DETECTOR
 
# Use JMergeDetector to merge ARCA28 pcal onto ARCA21 detector file (required to get correct modules on DU 14 and 15)
JMergeDetector \
    -a $WORKDIR/D0ARCA021_G_1.0.0.datx     \
    -@ "pcal=$DETECTOR"                     \
    -o $DETECTOR
 
# create TRIPOD for D0ARCA021 to be optimised in the back-propagation procedure.
typeset -A TRIPODS
 
JEditTripod \
    -f $WORKDIR/tripod_D0ARCA028.txt   \
    -o $TRIPOD_ABSOLUTE   # no emitters that have to be removed (all present in D0ARCA021)
 
get_tripods  $TRIPOD_INITIAL  TRIPODS;  FINAL_TRIPODS=(${(k)TRIPODS[*]})
get_tripods  $TRIPOD_ABSOLUTE TRIPODS;  FIXED_TRIPODS=(${(k)TRIPODS[*]})
 
cp -p  $TRIPOD_INITIAL  tripod.txt
 
JEditTripod   -f tripod.txt  -r "${FIXED_TRIPODS}" -o tripod.txt
JMergeTripod  -f tripod.txt  -f $TRIPOD_ABSOLUTE   -o tripod.txt
JEditTripod   -f tripod.txt  -k "${FINAL_TRIPODS}" -o tripod.txt
 
cat>acoustics_trigger_parameters.txt<<EOF
Q             = 0.0;
TMax_s        = 0.020;
numberOfHits  = 300;
EOF
 
cat>disable.txt<<EOF
# disable transmissions
# $e807fd96-3473-48a0-a5d8-5cbe0cc1ebc4$
# GIT 18.3.0-305-g4e82b04f4-D
# ROOT 6.30/02
# application JEditDisable
# command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
# namespace KM3NET
# system Linux verster.lorentz.leidenuniv.nl 6.6.8-200.fc39.x86_64 #1 SMP PREEMPT_DYNAMIC Thu Dec 21 04:01:49 UTC 2023 x86_64
2 808470373
2 808961261
2 808964925
2 808976233
3 808430571
3 808447094
3 808451782
3 808470373
3 808474602
3 808961261
3 808964925
3 808971936
3 808976233
3 808978646
3 808981369
3 808985803
4 808474602
4 808961261
4 808964925
5 808447094
5 808451782
5 808976207
5 808976233
6 808447094
6 808964808
7 806476519
7 808961261
7 808974966
7 808976319
8 808430571
8 808956926
8 808964757
8 808974966
8 808980070
8 817612559
9 808454978
9 808474602
9 808958576
9 808981206
9 808981369
9 808987098
9 816919071
9 817595373
9 817802210
10 808430036
10 808437870
10 808956926
10 808964876
10 808971766
10 808977563
10 817616429
10 817802231
12 808454978
12 808474243
12 808491794
12 808982550
12 808995481
12 808998966
13 808430036
13 808430571
13 808447081
13 808451689
13 808470373
13 808474243
13 808474602
13 808952499
13 808964925
13 808974966
13 808976233
13 808976352
13 808981369
13 808984662
13 817802231
EOF
 
cat>mechanics.txt<<EOF
-1 0.00094 294.291
EOF
 
cat>acoustics_fit_parameters.txt<<'EOF'
# $fd0a2502-e154-4566-b8fb-63f49e5ffffd$
# GIT test_elongated_shower_pde-123-ged3beb762-D
# ROOT 6.24/06
# application JEditProperties
# command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditProperties -f acoustics_fit_parameters.txt -o acoustics_fit_parameters.txt -@ \"mestimator  =   1;\" -@ \"option      =   3;\" -@ \"stdev       =  10.0e+5;\" -d 2
# namespace KM3NET
# system Linux doesburg.lorentz.leidenuniv.nl 5.15.13-200.fc35.x86_64 #1 SMP Wed Jan 5 16:39:13 UTC 2022 x86_64
Nmin         = 3;
Tmax_s       = 600.0;
chi2perNDF   = 100.0;
mestimator   = 1;
option       = 3;
sigma_s      = 50.0e-6;
stdev        = 10.0e+5;\
EOF
 
cat>sound_velocity.txt<<EOF
1541.00 -17e-3 -2000.00
EOF
 
# Build acoustic events for D0ARCA028 pre-calibration runs
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/$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
 
# Refit string 14, 15 and AB12 
fix  string  5 6 7 9 10 11 12 13 16 19 20 21 22 23 24 25 26 27 28 30 32
fix  tripod  7 8 9 10 13
 
stage  1A   -3  2  250.0e-6  10.0  20  0.5    0.5     
stage  2A   -3  2  100.0e-6  10.0  20  0.2    0.2
stage  3A   -3  2   50.0e-6  10.0  20  0.2    0.2
stage  4A    3  2   50.0e-6  10.0  20  0.2    0.2
 
EOF
 
JSydney.sh  $DETECTOR  $INPUT_FILES[*]  script.txt
 
# Apply creep correction to D0ARCA028 detector file
JKatoomba.sh  $DETECTOR  $INPUT_FILES[*]  katoomba.root
 
JCreep \
    -a $DETECTOR            \
    -f katoomba.root        \
    -A    
 
mkdir -p $WORKDIR/5JCreep
cp -p $DETECTOR $WORKDIR/5JCreep/$DETECTOR              
 
if [ ! $HOMEDIR -ef $WORKDIR ]; then
    cp -vru * $HOMEDIR
fi