Jpp - 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+ BACKGROUND ACOUSTICS_BACKGROUND: -B
14 
15 if do_usage $*; then
16  usage "$script <detector file> (run[-run])+"
17 fi
18 
19 if (( $# == 2 )); then
21  set_array RUNS $argv[2]
22 else
23  fatal "Wrong number of arguments."
24 fi
25 
26 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
27 eval `JPrintDetector -a $DETECTOR -O SUMMARY`
28 
29 source JAcoustics.sh $DETECTOR_ID
30 
31 CHECK_EXIT_CODE
32 
33 for RANGE in $RUNS[*]; do
34 
35  for (( RUN = ${RANGE%%-*}; $RUN <= ${RANGE##*-}; RUN += 1 )); do
36 
37  set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toashort.root
38  set_variable OUTPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
39  set_variable MONITOR_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
40 
41  if [[ ! -f $INPUT_FILE ]]; then
42 
43  JConvertDB \
44  -q toashort \
45  -@ "detid = `getDetector -D $DETECTOR_ID`" \
46  -@ "minrun = $RUN" \
47  -@ "maxrun = $RUN" \
48  -o $INPUT_FILE \
49  -d $DEBUG --!
50  fi
51 
52  if [[ ! -f $OUTPUT_FILE ]]; then
53 
54  JAcousticsEventBuilder \
55  -a $DETECTOR \
56  -f $INPUT_FILE \
57  -T $WORKDIR/tripod.txt \
58  -V $WORKDIR/sound_velocity.txt \
59  -W $WORKDIR/waveform.txt \
60  -@ $WORKDIR/acoustics_trigger_parameters.txt \
61  -H $WORKDIR/hydrophone.txt \
62  -o $OUTPUT_FILE \
63  -d $DEBUG --!
64  fi
65 
66  if [[ ! -f $MONITOR_FILE ]]; then
67 
68  $JPP_DIR/examples/JAcoustics/JAcousticsMonitor \
69  -a $DETECTOR \
70  -f $OUTPUT_FILE \
71  -o $MONITOR_FILE \
72  -d $DEBUG --!
73  fi
74 
75  JPlot1D \
76  -f ${MONITOR_FILE}:H\.\*size \
77  -L TL \
78  -\> "number of sensors" \
79  -\^ "number of events" \
80  -T "" \
81  -o $WORKDIR/size_${RUN}.gif $BACKGROUND
82 
83  JPlot1D \
84  -f ${MONITOR_FILE}:H\.\*multiplicity \
85  -L TR \
86  -\> "multiplicity" \
87  -\^ "number of events" \
88  -Y \
89  -T "" \
90  -o $WORKDIR/multiplicity_${RUN}.gif $BACKGROUND
91 
92  JPlot1D \
93  -f ${MONITOR_FILE}:H\.\*time \
94  -L TR \
95  -\> "#Deltat [s]" \
96  -\^ "number of events" \
97  -XX \
98  -T "" \
99  -o $WORKDIR/time_${RUN}.gif $BACKGROUND
100 
101  JPlot1D \
102  -f ${MONITOR_FILE}:H\.\*quality \
103  -L TL \
104  -\> "10log(quality)" \
105  -\^ "number of events" \
106  -XX \
107  -T "" \
108  -o $WORKDIR/quality_${RUN}.gif $BACKGROUND
109 
110  typeset -A TRIPODS
111 
112  get_tripods $WORKDIR/tripod.txt TRIPODS
113 
114  for EMITTER in ${(@k)TRIPODS}; do
115 
116  JPlot2D \
117  -f "${MONITOR_FILE}:H\[${EMITTER}\].event" \
118  -\> "string" \
119  -\< "floor" \
120  -O COLZ \
121  -T "[${EMITTER}]" \
122  -o $WORKDIR/event\[${EMITTER}\]_${RUN}.gif $BACKGROUND
123 
124  JPlot2D \
125  -f "${MONITOR_FILE}:H\[${EMITTER}\].doubles"\
126  -\> "string" \
127  -\< "floor" \
128  -O COLZ \
129  -T "[${EMITTER}]" \
130  -o $WORKDIR/doubles\[${EMITTER}\]_${RUN}.gif $BACKGROUND
131  done
132 
133  JPlot1D \
134  -w 1200x600 \
135  -f "${MONITOR_FILE}:G\[.*\].size" \
136  -y "0 150" \
137  -\^ "size" \
138  -t "$TIMESTAMP" \
139  -L BR \
140  -T "$RUN" \
141  -o $WORKDIR/size-${RUN}.gif $BACKGROUND
142 
143  JPlot1D \
144  -w 1200x600 \
145  -f "${MONITOR_FILE}:G\[.*\].rms" \
146  -y "0 0.025" \
147  -\^ "rms [s]" \
148  -t "$TIMESTAMP" \
149  -L TR \
150  -T "$RUN" \
151  -o $WORKDIR/rms-${RUN}.gif $BACKGROUND
152 
153  JPlot1D \
154  -w 1200x600 \
155  -f "${MONITOR_FILE}:G\[.*\].quality" \
156  -y "0 6" \
157  -\^ "10log(quality)" \
158  -t "$TIMESTAMP" \
159  -L BR \
160  -T "$RUN" \
161  -o $WORKDIR/quality-${RUN}.gif $BACKGROUND
162 
163  done
164 done
data_type w[N+1][M+1]
Definition: JPolint.hh:741
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
then usage else fatal Wrong number of arguments fi JCookie sh eval JPrintDetector a $DETECTOR O IDENTIFIER source JAcoustics sh $DETECTOR_ID CHECK_EXIT_CODE for RANGE in $RUNS[*]
set_variable INPUT_FILE
static const double H
Planck constant [eV s].
exit
Definition: JPizza.sh:36
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
then JPizza f
Definition: JPizza.sh:46
then echo
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:20
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then rm i $OUTPUT_FILE fi let RUN
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
* usage
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
&set_variable OUTPUT_FILE
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then JCalibrateToT a
Definition: JTuneHV.sh:116
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
then usage $script< input_file >< detector_file > fi set_variable OUTPUT_DIR set_variable SELECTOR JDAQTimesliceL1 set_variable DEBUG case set_variable DEBUG
then display $WORKDIR
Definition: plot-Domino.sh:127
set_array RUNS
Definition: JDataQuality.sh:23
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR tmp set_variable JDAQ_TIMESLICE JDAQTimeslice if do_usage *then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:23
set_variable DETECTOR
do set_variable INPUT_FILE $WORKDIR KM3NeT_
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:38
script
Definition: JAcoustics.sh:2
source $JPP_DIR setenv csh $JPP_DIR eval JShellParser o a A
esac $JPP_BIN JLogger sh $LOGGER until pgrep JGetMessage</dev/null > dev null
esac done
Definition: JAddHDE.sh:21