Jpp  master_rocky-40-g5f0272dcd
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  CHECK_EXIT_CODE
73 
74  timer_lap
75 fi
76 
77 
78 # Procedure to determine Jpp version used during data taking as default handling of UDP packet losses has changed.
79 # The given detector file is accordingly updated.
80 
81 JPrintMeta -f ${DAQ_FILE} -k GIT | read DAQ_VERSION
82 
83 echo ${(s/./)DAQ_VERSION} | read DAQ_MAJOR DAQ_MINOR DAQ_PATCH
84 
85 if [[ "$DAQ_MAJOR" != [0-9][0-9]* ]] || (( $DAQ_MAJOR <= 12 )); then # Jpp-version dependence
86 
87  notice "Set UDP counter and trailer disable."
88 
89  JEditDetector \
90  -a ${DETECTORB_2} \
91  -P "-1 -1 set UDP_COUNTER_DISABLE" \
92  -P "-1 -1 set UDP_TRAILER_DISABLE" \
93  -o ${DETECTORB_2}
94 
95  CHECK_EXIT_CODE
96 fi
97 
98 
99 JTriggerEfficiency \
100  -a ${DETECTORA_2} \
101  -b ${DETECTORB_2} \
102  -f"${INPUT_FILES[*]}" \
103  -o ${OUTPUT_FILE} \
104  -P ${PMT_FILE} \
105  -P "QE=${FACTOR}" \
106  ${OPTION} \
107  -S ${SEED} \
108  -B "$BACKGROUND_HZ" \
109  -r "${RUNBYRUN}" \
110  -d ${DEBUG} \
111  --!
112 
113 CHECK_EXIT_CODE
114 
115 timer_lap
116 
117 
118 JPrintMeta \
119  -f ${OUTPUT_FILE}
120 
121 JPrintTree \
122  -f ${OUTPUT_FILE}
123 
124 JTriggerMonitor \
125  -f ${OUTPUT_FILE} \
126  -o /dev/null
127 
128 
129 rm -f ${DETECTORA_2}
130 rm -f ${DETECTORB_2}