Jpp  19.1.0-rc.1
the software that should make you happy
pre-calibration_D0ARCA009.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 # --------------------------------------------------------------------------------------------
5 #
6 # Steering script for the global-fit-of-global-fits procedure of the D0ARCA009 (94) detector.
7 #
8 # --------------------------------------------------------------------------------------------
9 
10 if [ -z $JPP_DIR ]; then
11  echo "Variable JPP_DIR undefined."
12  exit
13 fi
14 
15 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
16 
17 set_variable DETECTOR_ID 94
18 set_array RUNS 10360 10361 10363 10370 10390 10391 10411 10412
19 set_variable: MUL ACOUSTICS_MUL -1.21e-2
20 set_variable: THREADS ACOUSTICS_THREADS 5
21 set_variable SINGULARITY -
22 set_variable WORKDIR `pwd`
23 set_variable HOMEDIR $WORKDIR
24 
25 if do_usage $*; then
26  usage "$script <source directory> [user directory]"\
27  "\nThe source directory corresponds to a pre-calibrated D0ARCA006 detector."\
28  "\nThe user directory should be specified when this script is submitted in a batch queue."
29 fi
30 
31 case $# in
32  2) set_variable HOMEDIR $2;&
33  1) set_variable DIR ${1:a};;
34  *) fatal "Wrong number of arguments."
35 esac
36 
37 source JAcousticsToolkit.sh
38 
39 set_variable DETECTOR $DIR/detector.datx
40 set_variable TRIPOD $DIR/tripod.txt
41 set_variable DETECTOR_INITIAL detector_initial.datx
42 set_variable TRIPOD_INITIAL tripod_initial.txt
43 expand_array RUNS
44 
45 
46 # Run in HOMEDIR
47 
48 mkdir -p $HOMEDIR
49 cd $HOMEDIR
50 
51 JDetectorDB \
52  -D $DETECTOR_ID \
53  -@ "tcal = A03209429" \
54  -@ "pcal = A03209525" \
55  -@ "rcal = A03209511" \
56  -V "" \
57  -WW \
58  -o $DETECTOR_INITIAL
59 
60 cat>$TRIPOD_INITIAL<<EOF
61 # Very preliminary acoustic emitter positions from survey
62 3 +587848.700 +4016749.700 -3450.467
63 4 +587657.170 +4017113.940 -3459.600
64 5 +587906.500 +4017094.300 -3459.600
65 6 +586939.550 +4017018.130 -3447.020
66 7 +587198.628 +4016228.693 -3433.306
67 EOF
68 
69 eval `JPrintDetector -a $DETECTOR_INITIAL -O SUMMARY`; FINAL_STRINGS=($STRINGS[*])
70 eval `JPrintDetector -a $DETECTOR -O SUMMARY`; FIXED_STRINGS=($STRINGS[*])
71 
72 # Apply modifications.
73 
74 JEditDetector \
75  -a $DETECTOR_INITIAL \
76  -s "-1 mul $MUL" \
77  -o detector.datx --!
78 
79 eval `JPrintDetector -a detector.datx -O SUMMARY`
80 
81 for STRING in ${STRINGS[*]}; do
82 
83  set_variable MODULE `getModule -a detector.datx -L "$STRING 0"`
84 
85  JEditDetector \
86  -a detector.datx \
87  -M "$MODULE setz -11.7" \
88  -o detector.datx --!
89 done
90 
91 JEditDetector -a detector.datx -r "$FIXED_STRINGS[*]" -o detector.datx --!
92 JMergeDetector -a detector.datx -a $DETECTOR -o detector.datx --!
93 JEditDetector -a detector.datx -k "$FINAL_STRINGS[*]" -o detector.datx --!
94 
95 typeset -A TRIPODS
96 
97 get_tripods $TRIPOD_INITIAL TRIPODS; FINAL_TRIPODS=(${(k)TRIPODS[*]})
98 get_tripods $TRIPOD TRIPODS; FIXED_TRIPODS=(${(k)TRIPODS[*]})
99 
100 cp -p $TRIPOD_INITIAL tripod.txt
101 
102 JEditTripod -f tripod.txt -r "${FIXED_TRIPODS}" -o tripod.txt
103 JMergeTripod -f tripod.txt -f $TRIPOD -o tripod.txt
104 JEditTripod -f tripod.txt -k "${FINAL_TRIPODS}" -o tripod.txt
105 
106 cat>acoustics_trigger_parameters.txt<<EOF
107 Q = 0.0;
108 TMax_s = 0.020;
109 numberOfHits = 120;
110 EOF
111 
112 cat>disable.txt<<EOF
113 # disable transmissions
114 # $f4cf86ad-c521-44f1-aaab-9a8474b47d15$
115 # GIT 18.3.0-305-g4e82b04f4-D
116 # ROOT 6.30/02
117 # application JEditDisable
118 # command /data1/dejong/km3net/GIT/Jpp/out//Linux/bin//JEditDisable -f disable.log -o disable.log -a detector.datx -q -d 2 --!
119 # namespace KM3NET
120 # 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
121 2 808470373
122 2 808961261
123 2 808964925
124 2 808976233
125 3 808430571
126 3 808447094
127 3 808451782
128 3 808470373
129 3 808474602
130 3 808961261
131 3 808964925
132 3 808971936
133 3 808976233
134 3 808978646
135 3 808981369
136 3 808985803
137 4 808474602
138 4 808961261
139 4 808964925
140 5 808447094
141 5 808451782
142 5 808976207
143 5 808976233
144 6 808447094
145 6 808964808
146 7 808961261
147 7 808976319
148 EOF
149 
150 cat>transmitter.txt<<EOF
151 # $script
152 EOF
153 
154 cat>hydrophone.txt<<EOF
155 # Data from D0ARCA009_hydrophones.csv
156  9 -1 -0.57 -0.42 0.57
157  10 0 -0.67 0.39 0.62
158  11 0 -0.67 0.39 0.62
159  13 0 -0.68 -0.37 0.62
160  14 0 -0.20 -0.75 0.62
161  18 0 -0.70 -0.34 0.62
162  12 0 -0.12 0.76 0.62
163  15 0 -0.72 -0.27 0.62
164  19 0 -0.16 -0.76 0.62
165 EOF
166 
167 JAcousticsEventBuilder.sh detector.datx $RUNS[*]
168 
169 INPUT_FILES=(`ls KM3NeT_${(l:8::0::0:)DETECTOR_ID}_0*${^RUNS}_event.root`)
170 
171 echo $INPUT_FILES[*]
172 
173 
174 # Run in WORKDIR
175 
176 cd $WORKDIR
177 
178 if [ ! $HOMEDIR -ef $WORKDIR ]; then
179  cp -p $HOMEDIR/$ACOUSTICS_DETECTOR $WORKDIR
180  cp -p $HOMEDIR/${^ACOUSTICS_KEYS}.txt $WORKDIR
181  cp -p $HOMEDIR/${^ACOUSTICS_AUXS}.txt $WORKDIR
182  cp -p $HOMEDIR/${^INPUT_FILES} $WORKDIR
183 fi
184 
185 cat>script.txt<<EOF
186 
187 fix string $FIXED_STRINGS[*]
188 fix tripod $FIXED_TRIPODS[*]
189 
190 stage 0 0 2 250.0e-6 10.0 20 0.5
191 stage 1A 0 2 250.0e-6 10.0 20 0.5 0.5
192 stage 1B 0 2 250.0e-6 10.0 0 0.0005 0.2
193 stage 1C ${SINGULARITY}2 0 250.0e-6 10.0 20 0.3
194 stage 2A ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
195 stage 2C ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
196 
197 initialise
198 
199 stage 2a ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2 0.2
200 stage 2c ${SINGULARITY}2 0 100.0e-6 10.0 20 0.2
201 stage 3a ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2 0.2
202 stage 3c ${SINGULARITY}2 0 50.0e-6 10.0 20 0.2
203 EOF
204 
205 JSydney.sh detector.datx $INPUT_FILES[*] script.txt
206 
207 if [ ! $HOMEDIR -ef $WORKDIR ]; then
208  cp -vru * $HOMEDIR
209 fi
Generated by doxygen 1.9.1