Jpp
 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 
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>"
17 fi
18 
19 if (( $# > 1 )); then
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 set_variable INPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_toashort.root
32 set_variable OUTPUT_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_event.root
33 set_variable MONITOR_FILE $WORKDIR/KM3NeT_${(l:8::0::0:)DETECTOR_ID}_${(l:8::0::0:)RUN}_monitor.root
34 
35 if [[ ! -f $INPUT_FILE ]]; then
36 
37  JConvertDB \
38  -q toashort \
39  -@ "detid = `getDetector -D $DETECTOR_ID`" \
40  -@ "minrun = $RUN" \
41  -@ "maxrun = $RUN" \
42  -o $INPUT_FILE \
43  -d $DEBUG --!
44 fi
45 
46 if [[ ! -f $OUTPUT_FILE ]]; then
47 
48  eval JAcousticsEventBuilder \
49  -a $DETECTOR \
50  -f $INPUT_FILE \
51  -@ `make_equation TRIGGER` \
52  -T `print_tripods` \
53  -V `print_velocity` \
54  -H $WORKDIR/hydrophone.txt \
55  -o $OUTPUT_FILE \
56  -d $DEBUG --!
57 fi
58 
59 if [[ ! -f $MONITOR_FILE ]]; then
60 
61  $JPP_DIR/examples/JAcoustics/JAcousticsMonitor \
62  -a $DETECTOR \
63  -f $OUTPUT_FILE \
64  -o $MONITOR_FILE \
65  -d $DEBUG --!
66 fi
67 
68 JPlot1D \
69  -f ${MONITOR_FILE}:H\.\*size \
70  -L TL \
71  -\> "number of sensors" \
72  -\^ "number of events" \
73  -T "" \
74  -o $WORKDIR/size_${RUN}.gif $BACKGROUND
75 
76 JPlot1D \
77  -f ${MONITOR_FILE}:H\.\*multiplicity \
78  -L TR \
79  -\> "multiplicity" \
80  -\^ "number of events" \
81  -Y \
82  -T "" \
83  -o $WORKDIR/multiplicity_${RUN}.gif $BACKGROUND
84 
85 JPlot1D \
86  -f ${MONITOR_FILE}:H\.\*time \
87  -L TR \
88  -\> "#Deltat [s]" \
89  -\^ "number of events" \
90  -XX \
91  -T "" \
92  -o $WORKDIR/time_${RUN}.gif $BACKGROUND
93 
94 JPlot1D \
95  -f ${MONITOR_FILE}:H\.\*quality \
96  -L TL \
97  -\> "10log(quality)" \
98  -\^ "number of events" \
99  -XX \
100  -T "" \
101  -o $WORKDIR/quality_${RUN}.gif $BACKGROUND
102 
103 typeset -A TRIPODS
104 
105 get_tripods $WORKDIR/tripod.txt TRIPODS
106 
107 for EMITTER in ${(@k)TRIPODS}; do
108 
109  JPlot2D \
110  -f "${MONITOR_FILE}:H\[${EMITTER}\].event" \
111  -\> "string" \
112  -\< "floor" \
113  -O COLZ \
114  -T "[${EMITTER}]" \
115  -o $WORKDIR/event\[${EMITTER}\]_${RUN}.gif $BACKGROUND
116 
117  JPlot2D \
118  -f "${MONITOR_FILE}:H\[${EMITTER}\].doubles"\
119  -\> "string" \
120  -\< "floor" \
121  -O COLZ \
122  -T "[${EMITTER}]" \
123  -o $WORKDIR/doubles\[${EMITTER}\]_${RUN}.gif $BACKGROUND
124 done
125 
126 JPlot1D \
127  -w 1200x600 \
128  -f "${MONITOR_FILE}:G\[.*\].size" \
129  -y "0 150" \
130  -\^ "size" \
131  -t "$TIMESTAMP" \
132  -L BR \
133  -T "$RUN" \
134  -o $WORKDIR/size-${RUN}.gif $BACKGROUND
135 
136 JPlot1D \
137  -w 1200x600 \
138  -f "${MONITOR_FILE}:G\[.*\].rms" \
139  -y "0 0.025" \
140  -\^ "rms [s]" \
141  -t "$TIMESTAMP" \
142  -L TR \
143  -T "$RUN" \
144  -o $WORKDIR/rms-${RUN}.gif $BACKGROUND
145 
146 JPlot1D \
147  -w 1200x600 \
148  -f "${MONITOR_FILE}:G\[.*\].quality" \
149  -y "0 6" \
150  -\^ "10log(quality)" \
151  -t "$TIMESTAMP" \
152  -L BR \
153  -T "$RUN" \
154  -o $WORKDIR/quality-${RUN}.gif $BACKGROUND
function get_tripods()
Definition: JAcoustics.sh:100
set_variable DETECTOR
data_type w[N+1][M+1]
Definition: JPolint.hh:741
do set_variable INPUT_FILE $WORKDIR KM3NeT_
void set_variable(const std::string &name, const std::string &value)
Set environment variable.
then fatal No sound hydrophone file $HYDROPHONE_TXT fi JGraph f $HYDROPHONE_TXT o $HYDROPHONE_ROOT sort gr k
function print_tripods()
Definition: JAcoustics.sh:89
static const double H
Planck constant [eV s].
exit
Definition: JPizza.sh:36
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:74
*fatal Invalid argument *esac fi function print_velocity()
Definition: JAcoustics.sh:81
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 echo
do set_variable OUTPUT_DIRECTORY $WORKDIR T
then rm i $OUTPUT_FILE fi let RUN
&set_variable OUTPUT_FILE
do $DIR JTransitTime o
Definition: JTransitTime.sh:44
* usage
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable ACOUSTICS_DEBUG
Definition: JAcoustics.sh:17
*set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2,-2] set_variable OUTPUT_FILE $argv[-1]
Definition: JDomino.sh:39
TRIGGER
Definition: JAcoustics.sh:47
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
then JCalibrateToT a
Definition: JTuneHV.sh:103
then display $WORKDIR
Definition: plot-Domino.sh:127
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 typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:36
script
Definition: JAcoustics.sh:2
set_variable INPUT_FILE
then usage $script< detector file >< inputfile > fi case set_variable WORKDIR
Definition: JLegolas.sh:28
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
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62