Jpp  19.1.0
the software that should make you happy
software/JTrigger/JTriggerEfficiencyRunByRun.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 # \author mdejong
4 #
5 script=${0##*/}
6 
7 # ------------------------------------------------------------------------------------------
8 #
9 # Script to run JTriggerEfficiency in run-by-run mode.
10 #
11 # ------------------------------------------------------------------------------------------
12 
13 if [ -z $JPP_DIR ]; then
14  echo "Variable JPP_DIR undefined."
15  exit
16 fi
17 
18 source $JPP_DIR/setenv.sh $JPP_DIR
19 
20 set_variable: DEBUG TRIGGEREFFICIENCY_DEBUG 1
21 set_variable: WORKDIR TRIGGEREFFICIENCY_WORKDIR ./
22 set_variable: BACKGROUND_HZ TRIGGEREFFICIENCY_BACKGROUND_HZ `getK40Rates`
23 set_variable: SEED TRIGGEREFFICIENCY_SEED 0
24 set_variable: NUMBER_OF_TIMESLICES TRIGGEREFFICIENCY_NUMBER_OF_TIMESLICES 0
25 set_variable: FACTOR TRIGGEREFFICIENCY_FACTOR 1.0
26 set_variable: SAMPLER TRIGGEREFFICIENCY_SAMPLER "100 1000"
27 set_variable+ OPTION TRIGGEREFFICIENCY_TRIGGERED_EVENTS_ONLY -O
28 set_variable+ DISABLE_TURBOT TRIGGEREFFICIENCY_DISABLE_TURBOT Y
29 
30 if do_usage $*; then
31  usage "$script <detector file> <detector file> (input file)+ <output file> <DAQ file> <PMT parameters file>"\
32  "\nAuxiliary script to trigger Monte Carlo events."
33 fi
34 
35 if (( $# >= 6 )); then
36  set_variable PMT_FILE $argv[-1]
37  set_variable DAQ_FILE $argv[-2]
38  set_variable OUTPUT_FILE $argv[-3]
39  set_array INPUT_FILES $argv[3,-4]
40  set_variable DETECTORB $argv[2]
41  set_variable DETECTORA $argv[1]
42 else
43  fatal "Wrong number of arguments."
44 fi
45 
46 timer_start
47 
48 JPrintTree -f $DAQ_FILE -@ "type = JDAQSummaryslice" -k number_of_entries | read N
49 
50 set_variable RUNBYRUN "file=$DAQ_FILE; sampler=${SAMPLER%% *} $(($N / 2))"
51 
52 set_variable DETECTORA_2 `mktemp $WORKDIR/${${DETECTORA:t}:r}_XXXXXX.datx`
53 set_variable DETECTORB_2 `mktemp $WORKDIR/${${DETECTORB:t}:r}_XXXXXX.datx`
54 
55 JConvertDetectorFormat -a ${DETECTORA} -o ${DETECTORA_2}
56 JConvertDetectorFormat -a ${DETECTORB} -o ${DETECTORB_2}
57 
58 
59 # Procedure to determine optical modules with out-sync-problem.
60 # The given detector file is accordingly updated.
61 
62 if [[ "$DISABLE_TURBOT" != "Y" ]]; then
63 
64  JTurbot \
65  -f ${DAQ_FILE} \
66  -a ${DETECTORA_2} \
67  -A \
68  -d ${DEBUG} \
69  -N ${NUMBER_OF_TIMESLICES} \
70  --!
71 
72  timer_lap
73 fi
74 
75 
76 # Procedure to determine Jpp version used during data taking as default handling of UDP packet losses has changed.
77 # The given detector file is accordingly updated.
78 
79 JPrintMeta -f ${DAQ_FILE} -k GIT | read DAQ_VERSION
80 
81 echo ${(s/./)DAQ_VERSION} | read DAQ_MAJOR DAQ_MINOR DAQ_PATCH
82 
83 if [[ "$DAQ_MAJOR" != [0-9][0-9]* ]] || (( $DAQ_MAJOR <= 12 )); then # Jpp-version dependence
84 
85  notice "Set UDP counter and trailer disable."
86 
87  JEditDetector \
88  -a ${DETECTORB_2} \
89  -P "-1 -1 set UDP_COUNTER_DISABLE" \
90  -P "-1 -1 set UDP_TRAILER_DISABLE" \
91  -o ${DETECTORB_2}
92 fi
93 
94 
95 JTriggerEfficiency \
96  -a ${DETECTORA_2} \
97  -b ${DETECTORB_2} \
98  -f"${INPUT_FILES[*]}" \
99  -o ${OUTPUT_FILE} \
100  -P ${PMT_FILE} \
101  -P "QE=${FACTOR}" \
102  ${OPTION} \
103  -S ${SEED} \
104  -B "$BACKGROUND_HZ" \
105  -r "${RUNBYRUN}" \
106  -d ${DEBUG} \
107  --!
108 
109 timer_lap
110 
111 
112 JPrintMeta \
113  -f ${OUTPUT_FILE}
114 
115 JPrintTree \
116  -f ${OUTPUT_FILE}
117 
118 JTriggerMonitor \
119  -f ${OUTPUT_FILE} \
120  -o /dev/null
121 
122 
123 rm -f ${DETECTORA_2}
124 rm -f ${DETECTORB_2}