Jpp  18.2.0-rc.1
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JRecalibrateK40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 script=${0##*/}
3 
4 # ------------------------------------------------------------------------------------------
5 #
6 # Script to run K40 time calibration following HV tuning.
7 #
8 # ------------------------------------------------------------------------------------------
9 
10 if [ -z $JPP_DIR ]; then
11  echo "Variable JPP_DIR undefined."
12  exit
13 fi
14 
15 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
16 
17 set_variable: DEBUG CALIBRATION_DEBUG 1
18 set_variable: WORKDIR CALIBRATION_WORKDIR `pwd`
19 set_variable: DAQ_TIMESLICE CALIBRATION_TIMESLICE JDAQTimesliceL1
20 set_variable+ DATABASE_OPTION CALIBRATION_USEPMTID -U
21 set_variable+ MODEL CALIBRATION_MODEL -D
22 
23 if do_usage $*; then
24  usage "$script <detector file> <input file> <input file> <input file>"\
25  "\nThe 1st and 3rd input file should correspond to a data taking run in which the HV settings are complete (e.g. \"vendor\" or \"tuned\")."\
26  "\nThe 2nd input file should correspond to a run with hybrid HV settings."\
27  "\nIn this, half of the PMTs should have the same HV settings as either the 1st or 3rd input file."
28 fi
29 
30 if (( $# == 4 )); then
31  set_variable DETECTOR $argv[1]
32  set_array INPUT_FILES $argv[2,-1]
33 else
34  fatal "Wrong number of arguments."
35 fi
36 
37 JCookie.sh
38 
39 #
40 # Get background option for given input file.
41 #
42 # \param 1 input file
43 # \return option
44 #
45 function get_option()
46 {
47  if (( `JPrintTree -f $1 | awk '/JDAQSummaryslice/ { print $3 } '` > 0 )); then
48  echo rates
49  else
50  echo counts
51  fi
52 }
53 
54 
55 #
56 # Standard calibration.
57 #
58 # \param 1 detector
59 # \param 2 input file
60 #
61 function calibrate()
62 {
63  JPrintDAQHeader -f $2 | read ID RUN FRAME_INDEX UTC
64 
65  JCalibrateK40 \
66  -a $1 \
67  -f $2 \
68  -C $DAQ_TIMESLICE \
69  -o $WORKDIR/monitor_${RUN}.root \
70  -b `get_option $2` \
71  -d $DEBUG --!
72 
73  JMergeCalibrateK40 \
74  -f $WORKDIR/monitor_${RUN}.root \
75  -o $WORKDIR/merge_${RUN}.root \
76  -@ "LIVETIME_S = 10" \
77  -d $DEBUG
78 
79  JConvertDetectorFormat \
80  -a $1 \
81  -o $WORKDIR/detector_${RUN}.datx
82 
83  JFitK40 \
84  -a $WORKDIR/detector_${RUN}.datx \
85  -f $WORKDIR/merge_${RUN}.root \
86  -o $WORKDIR/fit_${RUN}.root \ $MODEL \
87  -A \
88  -w \
89  -d $DEBUG
90 }
91 
92 
93 #
94 # Hybrid calibration.
95 #
96 # \param 1 detector
97 # \param 2 reference run
98 # \param 3 input file
99 #
100 function hybrid()
101 {
102  JPrintDAQHeader -f $3 | read ID RUN FRAME_INDEX UTC
103 
104  JTDC \
105  -D $DETECTOR_ID \
106  -r "$2 $RUN" \
107  $DATABASE_OPTION \
108  -o "$WORKDIR/TDC_${2}_${RUN}.txt"
109 
110  JCalibrateK40 \
111  -a $1 \
112  -f $3 \
113  -C $DAQ_TIMESLICE \
114  -o $WORKDIR/monitor_${RUN}.root \
115  -b `get_option $3` \
116  -d $DEBUG --!
117 
118  JMergeCalibrateK40 \
119  -f $WORKDIR/monitor_${RUN}.root \
120  -o $WORKDIR/merge_${RUN}.root \
121  -@ "LIVETIME_S = 10" \
122  -d $DEBUG
123 
124  JConvertDetectorFormat \
125  -a $1 \
126  -o $WORKDIR/detector_${RUN}.datx
127 
128  JFitK40 \
129  -a $WORKDIR/detector_${RUN}.datx \
130  -f $WORKDIR/merge_${RUN}.root \
131  -o $WORKDIR/fit_${RUN}.root \ $MODEL \
132  -A \
133  -w \
134  -! $WORKDIR/TDC_${2}_${RUN}.txt \
135  -d $DEBUG
136 }
137 
138 
139 # Get run numbers
140 
141 eval `JPrintDetector -a $DETECTOR -O IDENTIFIER`
142 
143 RUNS=()
144 
145 for INPUT_FILE in $INPUT_FILES[*]; do
146 
147  JPrintDAQHeader -f $INPUT_FILE | read ID RUN FRAME_INDEX UTC
148 
149  if (( $ID != $DETECTOR_ID )); then
150  fatal "Detector identifier from $INPUT_FILE $ID != $DETECTOR_ID"
151  fi
152 
153  RUNS+=($RUN)
154 
155 done
156 
157 
158 calibrate $DETECTOR ${INPUT_FILES[1]}
159 
160 hybrid $WORKDIR/detector_${RUNS[1]}.datx ${RUNS[1]} ${INPUT_FILES[2]}
161 
162 hybrid $WORKDIR/detector_${RUNS[2]}.datx ${RUNS[2]} ${INPUT_FILES[3]}
163 
DAQ_TIMESLICE
&set_variable DAQ_TIMESLICE
Definition: JLegolas.sh:32
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:778
$argv
then usage $script< detector identifier >< startdate\"YYYY-MM-DDHH:MM:SS\"><finaldate\"YYYY-MM-DDHH:MM:SS\"><QA/QCfile> fi case set_variable QAQC_TXT $argv[4]
Definition: JDataMonitor.sh:24
RUNS
set_array RUNS
Definition: JOperationsRunMonitor.sh:24
JCOMPAREHISTOGRAMS::read
std::istream & read(std::istream &in, JTestSummary &summary, const char delimiter= ' ')
Read test summary.
Definition: JTestSummary.hh:222
f
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
exit
exit
Definition: JPizza.sh:36
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
JPHYSICS::C
static const double C
Physics constants.
Definition: JPhysics/JConstants.hh:25
set_variable
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable
Definition: JAcoustics.sh:21
awk
do set_variable STRING_TXT awk
Definition: plot-canberra.sh:195
a
then JCalibrateToT a
Definition: JTuneHV.sh:113
RUN
then rm i $OUTPUT_FILE fi let RUN
Definition: software/JDataWriter/JDataWriter.sh:103
$INPUT_FILES
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR set_variable DEBUG set_variable WORKDIR if do_usage *then usage for INPUT_FILE in $INPUT_FILES[*]
Definition: JCheckTriggerParameters.sh:31
WORKDIR
then echo Variable JPP_DIR undefined exit fi source $JPP_DIR setenv sh $JPP_DIR &dev null set_variable DEBUG set_variable WORKDIR
Definition: JLegolas.sh:20
usage
* usage
Definition: JLigier-local.sh:34
$INPUT_FILE
then JFileTuna f $INPUT_FILE
Definition: file-Tuna.sh:66
INPUT_FILE
set_variable INPUT_FILE
Definition: JARCACalibrateMuon.sh:30
INPUT_FILES
do sput $STRING $FLOOR sget MODULE INPUT_FILES
Definition: JCanberra.sh:77
$ID
then ls $DIR $ID
Definition: getArchive.sh:85
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:47
$DETECTOR
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
Definition: module-Z:fit.sh:118
$WORKDIR
then display $WORKDIR
Definition: plot-Domino.sh:128
o
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
Definition: compare-detector.sh:31
_
then eval ls $DIR $ID $TYPE _
Definition: getArchive.sh:104
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
echo
then echo
Definition: JQAQC.sh:90
JLANG::null
static JNullStream null
Null I/O stream.
Definition: JNullStream.hh:51
in
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46
A
source $JPP_DIR setenv csh $JPP_DIR &dev null eval JShellParser o a A
Definition: JShellParser.csh:15
script
script
Definition: JAcoustics.sh:2
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
done
esac done
Definition: JAddHDE.sh:21
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5