Jpp  16.0.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JAcousticsMonitor.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 if [ -z $JPP_DIR ]; then
5  echo "Variable JPP_DIR undefined."
6  exit
7 fi
8 
9 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10 
11 set_variable: DEBUG ACOUSTICS_DEBUG 2
12 set_variable: WORKDIR ACOUSTICS_WORKDIR ./
13 set_variable: FORMAT GRAPHICS_FORMAT gif
14 set_variable+ BATCH GRAPHICS_BATCH -B
15 
16 if do_usage $*; then
17  usage "$script <detector file> (run[-run])+"
18 fi
19 
20 if (( $# == 2 )); then
21  set_variable DETECTOR $argv[1]
22  set_array RUNS $argv[2]
23 else
24  fatal "Wrong number of arguments."
25 fi
26 
27 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
28 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
29 
30 source JAcoustics.sh $DETECTOR_ID
31 
32 CHECK_EXIT_CODE
33 
34 for RANGE in $RUNS[*]; do
35 
36  for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
37 
38  set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toashort.root
39  set_variable OUTPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
40  set_variable MONITOR_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
41 
42  if [[ ! -f $INPUT_FILE ]]; then
43 
44  JConvertDB \
45  -q toashort \
46  -@ "detid = `getDetector -D $DETECTOR_ID`" \
47  -@ "minrun = $RUN" \
48  -@ "maxrun = $RUN" \
49  -o $INPUT_FILE \
50  -d $DEBUG --!
51  fi
52 
53  if [[ ! -f $OUTPUT_FILE ]]; then
54 
55  JAcousticsEventBuilder \
56  -a $DETECTOR \
57  -f $INPUT_FILE \
58  -T $WORKDIR/tripod.txt \
59  -V $WORKDIR/sound_velocity.txt \
60  -W $WORKDIR/waveform.txt \
61  -@ $WORKDIR/acoustics_trigger_parameters.txt \
62  -H $WORKDIR/hydrophone.txt \
63  -o $OUTPUT_FILE \
64  -d $DEBUG --!
65  fi
66 
67  $JPP_DIR/examples/JAcoustics/JAcousticsMonitor \
68  -a $DETECTOR \
69  -f $OUTPUT_FILE \
70  -o $MONITOR_FILE \
71  -d $DEBUG --!
72 
73  JPlot1D \
74  -f ${MONITOR_FILE}:H\.\*size \
75  -L TL \
76  -\> "number of sensors" \
77  -\^ "number of events" \
78  -T "" \
79  -o $WORKDIR/size_${RUN}.$FORMAT $BATCH
80 
81  JPlot1D \
82  -f ${MONITOR_FILE}:H\.\*multiplicity \
83  -L TR \
84  -\> "multiplicity" \
85  -\^ "number of events" \
86  -Y \
87  -T "" \
88  -o $WORKDIR/multiplicity_${RUN}.$FORMAT $BATCH
89 
90  JPlot1D \
91  -f ${MONITOR_FILE}:H\.\*time \
92  -L TR \
93  -\> "#Deltat [s]" \
94  -\^ "number of events" \
95  -XX \
96  -T "" \
97  -o $WORKDIR/time_${RUN}.$FORMAT $BATCH
98 
99  JPlot1D \
100  -f ${MONITOR_FILE}:H\.\*quality \
101  -L TL \
102  -\> "10log(quality)" \
103  -\^ "number of events" \
104  -XX \
105  -T "" \
106  -o $WORKDIR/quality_${RUN}.$FORMAT $BATCH
107 
108  typeset -A TRIPODS
109 
110  get_tripods $WORKDIR/tripod.txt TRIPODS
111 
112  for EMITTER in ${(@k)TRIPODS}; do
113 
114  JPlot2D \
115  -f "${MONITOR_FILE}:H\[${EMITTER}\].event" \
116  -\> "string" \
117  -\< "floor" \
118  -O COLZ \
119  -T "[${EMITTER}]" \
120  -o $WORKDIR/event\[${EMITTER}\]_${RUN}.$FORMAT $BATCH
121 
122  JPlot2D \
123  -f "${MONITOR_FILE}:H\[${EMITTER}\].doubles" \
124  -\> "string" \
125  -\< "floor" \
126  -O COLZ \
127  -T "[${EMITTER}]" \
128  -o $WORKDIR/doubles\[${EMITTER}\]_${RUN}.$FORMAT $BATCH
129  done
130 
131  JPlot1D \
132  -w 1200x600 \
133  -f "${MONITOR_FILE}:G\[.*\].size" \
134  -y "0 150" \
135  -\^ "size" \
136  -t "$TIMESTAMP" \
137  -L BR \
138  -T "$RUN" \
139  -o $WORKDIR/size-${RUN}.$FORMAT $BATCH
140 
141  JPlot1D \
142  -w 1200x600 \
143  -f "${MONITOR_FILE}:G\[.*\].rms" \
144  -y "0 0.025" \
145  -\^ "rms [s]" \
146  -t "$TIMESTAMP" \
147  -L TR \
148  -T "$RUN" \
149  -o $WORKDIR/rms-${RUN}.$FORMAT $BATCH
150 
151  JPlot1D \
152  -w 1200x600 \
153  -f "${MONITOR_FILE}:G\[.*\].deviation" \
154  -y "0 0.050" \
155  -\^ "deviation [s]" \
156  -t "$TIMESTAMP" \
157  -L TR \
158  -T "$RUN" \
159  -o $WORKDIR/deviation-${RUN}.$FORMAT $BATCH
160 
161  JPlot1D \
162  -w 1200x600 \
163  -f "${MONITOR_FILE}:G\[.*\].quality" \
164  -y "0 6" \
165  -\^ "10log(quality)" \
166  -t "$TIMESTAMP" \
167  -L BR \
168  -T "$RUN" \
169  -o $WORKDIR/quality-${RUN}.$FORMAT $BATCH
170 
171  done
172 done
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:757
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
$argv
*set_variable DETECTOR_ID $argv[1] set_array RUNS $argv[2,-2]
Definition: JDataQuality.sh:24
f
o $QUALITY_ROOT d $DEBUG!JPlot1D f
Definition: JDataQuality.sh:66
JPHYSICS::H
static const double H
Planck constant [eV s].
Definition: JPhysics/JConstants.hh:28
exit
exit
Definition: JPizza.sh:36
Y
then fatal Wrong number of arguments fi set_variable STRING $argv[1] set_variable DETECTORXY_TXT $WORKDIR $DETECTORXY_TXT tail read X Y CHI2 RMS printf optimum n $X $Y $CHI2 $RMS awk v Y
Definition: detector-XY:fit1d.sh:33
V
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
echo
then echo
Definition: JDOMDAQDriver.sh:56
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
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
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
XX
set_variable LOG XX
Definition: JMuonPostfit.sh:71
o
then JConvertDetectorFormat a $DETECTOR[1] o
Definition: JCompareDetector.sh:45
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
OUTPUT_FILE
&set_variable OUTPUT_FILE
Definition: JARCACalibrateMuon.sh:29
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
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
DEBUG
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
Definition: JCalibratePMT.sh:21
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
KM3NeT_
do set_variable INPUT_FILE $WORKDIR KM3NeT_
Definition: JAcousticsEventBuilder.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
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