Jpp  16.0.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
acoustics-fit_00000100.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 100
18 set_array RUNS 11342 11345 11347 11374 11376
19 #set_variable: STAGES ACOUSTICS_STAGES 1 2 3
20 
21 if do_usage $*; then
22  usage "$script [<detector file> <tripod file>] [user directory]"\
23  "\nThe user directory should be specified when this script is submitted in a batch queue."\
24  "\nThe optional detector file and tripod file correspond to a pre-calibrated sub-detector."
25 fi
26 
27 set_variable WORKDIR `pwd`
28 set_variable HOMEDIR $WORKDIR
29 
30 case $# in
31  3) set_variable HOMEDIR $3;&
32  2) set_variable TRIPOD $2;
33  set_variable DETECTOR $1;;
34  1) set_variable HOMEDIR $1;;
35  0) ;;
36  *) fatal "Wrong number of arguments."
37 esac
38 
39 set_variable TRIPOD_INITIAL tripod_initial.txt
40 set_variable DETECTOR_INITIAL detector_initial.datx
41 
42 # Run in HOMEDIR
43 
44 mkdir -p $HOMEDIR
45 
46 cd $HOMEDIR
47 
48 if [[ ! -f $DETECTOR_INITIAL ]]; then
49 
50  JDetectorDB \
51  -D $DETECTOR_ID \
52  -r $RUNS[1] \
53  -WW \
54  -o $DETECTOR_INITIAL
55 fi
56 
57 if [[ ! -f $TRIPOD_INITIAL ]]; then
58 
59  cat>$TRIPOD_INITIAL<<EOF
60 # Very preliminary acoustic emitter positions from survey
61 2 +256815.500 +4743395.000 -2436.677
62 3 +257096.200 +4743636.000 -2439.354
63 4 +256919.238 +4743354.178 -2439.354
64 EOF
65 fi
66 
67 
68 # Apply modifications to the detector and tripod file, if any.
69 
70 JEditDetector \
71  -a $DETECTOR_INITIAL \
72  -o ${TMPDIR:-/tmp}/detector.datx
73 
74 if [[ -n "$DETECTOR" ]]; then
75 
76  eval `JPrintDetector -a $DETECTOR -O SUMMARY`
77 
78  JEditDetector -a ${TMPDIR:-/tmp}/detector.datx -r "$STRINGS[*]" -o ${TMPDIR:-/tmp}/detector.datx
79 
80  JMergeDetector -a ${TMPDIR:-/tmp}/detector.datx -a $DETECTOR -o detector.datx
81 
82  set_variable ACOUSTICS_STRINGS $STRINGS[*] # fix strings in acoustics-fit.sh
83 else
84 
85  JConvertDetectorFormat -a ${TMPDIR:-/tmp}/detector.datx -o detector.datx
86 fi
87 
88 if [[ -n "$TRIPOD" ]]; then
89 
90  source JAcoustics.sh --
91 
92  typeset -A TRIPODS
93 
94  get_tripods $TRIPOD TRIPODS
95 
96  typeset -A TRIPODS_INITIAL
97 
98  get_tripods $TRIPOD_INITIAL TRIPODS_INITIAL
99 
100  cp -p $TRIPOD_INITIAL tripod.txt
101 
102  for ID in ${(@k)TRIPODS_INITIAL}; do
103  if [[ -n "$TRIPODS[$ID]" ]]; then
104  JEditTripod -f tripod.txt -T "$ID set $TRIPODS[$ID]"
105  fi
106  done
107 
108  set_variable ACOUSTICS_TRIPODS ${(@k)TRIPODS[*]} # fix tripods in acoustics-fit.sh
109 else
110 
111  cp -p $TRIPOD_INITIAL tripod.txt
112 fi
113 
114 set_variable DETECTOR detector.datx
115 set_variable TRIPOD tripod.txt
116 
117 
118 JAcoustics.sh $DETECTOR_ID
119 
120 
121 JAcousticsEventBuilder.sh $DETECTOR $RUNS[*]
122 
123 # Run in WORKDIR
124 
125 cd $WORKDIR
126 
127 if [ ! $HOMEDIR -ef $WORKDIR ]; then
128 
129  cp -p $HOMEDIR/$TRIPOD $WORKDIR
130  cp -p $HOMEDIR/$DETECTOR $WORKDIR
131 
132  for RANGE in $RUNS[*]; do
133  for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
134  cp -p $HOMEDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root $WORKDIR
135  done
136  done
137 
138  JAcoustics.sh $DETECTOR_ID
139 
140 fi
141 
142 for STAGE in 1; do
143 
144  $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*event.root > $HOMEDIR/stage-${STAGE}.log
145 
146  if [ ! $HOMEDIR -ef $WORKDIR ]; then
147  cp -vru * $HOMEDIR
148  fi
149 done
150 
151 unset_variable ACOUSTICS_STRINGS
152 unset_variable ACOUSTICS_TRIPODS
153 
154 for STAGE in 2 3; do
155 
156  $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_*event.root > $HOMEDIR/stage-${STAGE}.log
157 
158  if [ ! $HOMEDIR -ef $WORKDIR ]; then
159  cp -vru * $HOMEDIR
160  fi
161 done
162 
k
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Definition: hydrophone-phi:fit.sh:35
RUNS
set_array RUNS
Definition: JOperationsRunMonitor.sh:24
file
then set_variable PMT_FILE set_variable DAQ_FILE set_variable OUTPUT_FILE set_variable DETECTOR else fatal Wrong number of arguments fi set_variable RUNBYRUN file
Definition: JRandomEventRunByRun.sh:46
cat
then JLigier sh continue fi cat
Definition: JDAQDriver.sh:51
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
exit
exit
Definition: JPizza.sh:36
HOMEDIR
then usage $script[< detector file >< tripod file >][user directory] nThe user directory should be specified when this script is submitted in a batch queue nThe optional detector file and tripod file correspond to a pre calibrated sub detector fi set_variable WORKDIR pwd set_variable HOMEDIR $WORKDIR case set_variable HOMEDIR
Definition: acoustics-fit_00000094.sh:32
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:758
echo
then echo
Definition: JDOMDAQDriver.sh:56
$STRINGS
do JCanberra a $DETECTOR f $INPUT_FILE o $WORKDIR canberra[${EMITTER}\] root T $WORKDIR tripod txt V $WORKDIR sound_velocity txt M $WORKDIR mechanics txt H $WORKDIR hydrophone txt E $EMITTER $DISABLE d $DEBUG!done kill_child_processes_at_exit attach getModule a $DETECTOR typeset Z STRING typeset Z FLOOR for STRING in $STRINGS[*]
Definition: JCanberra.sh:64
JTOOLS::n
const int n
Definition: JPolint.hh:676
for
do for((RUN=${RANGE%%-*};$RUN<=${RANGE##*-};RUN+=1))
$STAGE
*fatal invalid stage $STAGE(possible stages:1, 2, 3, A, D)"
DETECTOR_ID
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null if do_usage *then usage $script< detector identifier > fi case set_variable DETECTOR_ID
Definition: JTDC.sh:11
$RUNS
then fatal Detector file should be binary format fi eval JPrintDetector a $DETECTOR O IDENTIFIER typeset a INPUT_FILES for RANGE in $RUNS[*]
Definition: software/JCompass/JCompass.sh:47
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
TMPDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR if do_usage *then usage $script[(input file)+] fi set_variable DEBUG set_variable WORKDIR TMPDIR
Definition: JEnergyCorrection.sh:21
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
RUN
then rm i $OUTPUT_FILE fi let RUN
Definition: software/JDataWriter/JDataWriter.sh:103
strings
then usage $script< detector specific acoustics-fit script >< option > nAuxiliary script to make scan of pre stretching of detector strings(see JEditDetector)." "\nPossible options
Definition: acoustics-Zmul.sh:20
TRIPOD
&set_variable TRIPOD
Definition: acoustics-fit_00000094.sh:33
done
set_variable NUMBER_OF_ITERATIONS set_variable EPSILON cat acoustics_fit_parameters txt<< EOF $CONFIGURATION[*]Nmin=3;sigma_s=100.0e-6;stdev=10.0;mestimator=0;fixStrings=0;EOF for STRING in $STRINGS[*];do#fit stretching and(z) position of given string set_variable DETECTOR_TMP ${TMPDIR:-/tmp}/detector_A.datx JEditDetector-a $DETECTOR-o $DETECTOR_TMP-r $STRING JEditDetector-a $DETECTOR-o $DETECTOR-k $STRING for MUL in 0.005 0.001;do DX_M=0.2 for((N=0;$N< $NUMBER_OF_ITERATIONS;++N));do CHI2[3]=$CHI2[1] fitPositionOfString $STRING Z $DX_M fitStretchingOfString $STRING $MUL if(($CHI2[3]-$CHI2[1]< $EPSILON));then break fi done if(($N >=$NUMBER_OF_ITERATIONS));then printf"warning: reached maximum number of iterations %d - converenge %7.3f\n"$N $(($CHI2[3]-$CHI2[1])) fi done JMergeDetector-a $DETECTOR-a $DETECTOR_TMP-o $DETECTOR rm-f $DETECTOR_TMP JConvertDetectorFormat-a $DETECTOR-o $DETECTOR-r-d 0 > &dev null done
Definition: acoustics-fit.sh:563
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
usage
* usage
Definition: JLigier-local.sh:34
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
$DETECTOR
then set_variable MODULE getModule a $DETECTOR L $STRING $FLOOR JEditDetector a $DETECTOR M $MODULE add $X o $DETECTOR else echo No update of detector $DETECTOR
Definition: module-Z:fit.sh:118
Run
then error Missing CDF files error Run
Definition: JSirene.sh:58
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
cp
then cp
Definition: darkroomCalibration.sh:41
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
KM3NeT_
do set_variable INPUT_FILE $WORKDIR KM3NeT_
Definition: JAcousticsEventBuilder.sh:41
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
tripods
then fatal Not enough tripods
Definition: hydrophone-phi:run.sh:41
in
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:42
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
script
script
Definition: JAcoustics.sh:2
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
stage
then cat $TRIPOD_INITIAL<< EOF1 256877.5 4743716.7-2438.42 256815.5 4743395.0-2435.53 257096.2 4743636.0-2439.5EOFfiJEditDetector-a $DETECTOR_INITIAL-@"z = -2441.5"-s"-1 mul $MUL"-o $DETECTORcp-p $TRIPOD_INITIAL $TRIPODJAcoustics.sh $DETECTOR_IDJAcousticsEventBuilder.sh $DETECTOR $RUNS[*]cd $WORKDIRif[!$HOMEDIR-ef $WORKDIR];then cp-p $HOMEDIR/$TRIPOD $WORKDIR cp-p $HOMEDIR/$DETECTOR $WORKDIR for RANGE in $RUNS[*];do for((RUN=${RANGE%%-*};$RUN<=${RANGE##*-};RUN+=1));do cp-p $HOMEDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ ${(l:8::0::0:) RUN}_event.root $WORKDIR done done JAcoustics.sh $DETECTOR_IDfifor STAGE in`echo $STAGES`;do $JPP_DIR/examples/JAcoustics/acoustics-fit.sh $DETECTOR $TRIPOD $STAGE $WORKDIR/KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_ *event.root > $HOMEDIR stage
Definition: acoustics-fit_00000049.sh:102
$TRIPOD
then fatal Invalid tripod $TRIPOD
Definition: tripod-XY:run.sh:38
null
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
Definition: examples/JLogger/JLogger.sh:33
Generated by   doxygen 1.8.5