Jpp
JDAQTimesliceSelector.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 version=1.0
6 script=${0##*/}
7 
8 # ------------------------------------------------------------------------------------------
9 #
10 # Example script for JDAQTimeslice.
11 #
12 # ------------------------------------------------------------------------------------------
13 
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 
21 source $JPP_DIR/setenv.sh $JPP_DIR
22 
23 
24 set_variable DEBUG 2
25 set_variable WORKDIR /tmp/
26 define_variable INPUT_FILE
27 set_variable OUTPUT_FILE histogram\[%\].root
28 
29 if ( do_usage $* ); then
30  usage "$script <input file> [output file]"
31 fi
32 
33 case $# in
34  2) set_variable OUTPUT_FILE $2;&
35  1) set_variable INPUT_FILE $1;;
36  *) fatal "Wrong number of arguments."
37 esac
38 
39 check_input_file $INPUT_FILE
40 
41 for TIMESLICE in JDAQTimeslice JDAQTimesliceL0 JDAQTimesliceL1 JDAQTimesliceL2 JDAQTimesliceSN; do
42 
43  echo Processing $INPUT_FILE with selection $TIMESLICE.
44 
45  $JPP_DIR/examples/JDAQ/JDAQTimesliceSelector \
46  -f ${INPUT_FILE} \
47  -o ${OUTPUT_FILE/\%/${TIMESLICE}} \
48  -C $TIMESLICE \
49  -d $DEBUG
50 done
51 
52 JPlot1D \
53  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h0 \
54  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h0 \
55  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h0 \
56  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h0 \
57  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h0 \
58  -> "PMT" \
59  -\^ "number of events [a.u.]" \
60  -T "" -L TR
61 
62 JPlot1D \
63  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h1 \
64  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h1 \
65  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h1 \
66  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h1 \
67  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h1 \
68  -> "time-over-threshold [ns]" \
69  -\^ "number of events [a.u.]" \
70  -T "" -L TR
71 
72 JPlot1D \
73  -f ${OUTPUT_FILE/\%/JDAQTimeslice}:h2 \
74  -f ${OUTPUT_FILE/\%/JDAQTimesliceL0}:h2 \
75  -f ${OUTPUT_FILE/\%/JDAQTimesliceL1}:h2 \
76  -f ${OUTPUT_FILE/\%/JDAQTimesliceL2}:h2 \
77  -f ${OUTPUT_FILE/\%/JDAQTimesliceSN}:h2 \
78  -> "time [ns]" \
79  -\^ "number of events [a.u.]" \
80  -w 800x400 \
81  -T "" -L TR