Jpp  18.0.0-rc.2
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JCalibrateK40.sh
Go to the documentation of this file.
1 #!/bin/zsh
2 #
3 #
4 # \author mdejong
5 #
6 version=1.0
7 script=${0##*/}
8 
9 # ------------------------------------------------------------------------------------------
10 #
11 # Utility script to test JCalibrateK40/JFitK40.
12 #
13 # ------------------------------------------------------------------------------------------
14 
15 if [ -z $JPP_DIR ]; then
16  echo "Variable JPP_DIR undefined."
17  exit
18 fi
19 
20 source $JPP_DIR/setenv.sh $JPP_DIR
21 
22 set_variable DEBUG 2
23 set_variable WORKDIR ${TMPDIR:-/tmp}/
24 set_variable: FORMAT GRAPHICS_FORMAT gif
25 set_variable+ BATCH GRAPHICS_BATCH -B
26 
27 if do_usage $*; then
28  usage "$script [working directory]"
29 fi
30 
31 case $# in
32  1) set_variable WORKDIR $1;;
33 esac
34 
35 set_variable DETECTOR $WORKDIR/detector.detx
36 set_variable BACKGROUND_HZ 5e3 440 50 5 0.5
37 set_variable QE 1.0
38 set_variable numberOfSlices 6000
39 
40 
41 if (( 1 )); then
42 
43  set_variable detectorType 1
44  set_variable numberOfLines 1
45  set_variable numberOfFloors 1
46 
47  JDetector \
48  -D 1 \
49  -@ "detectorType = $detectorType" \
50  -@ "numberOfLines = $numberOfLines" \
51  -@ "numberOfFloors = $numberOfFloors" \
52  -o $DETECTOR -V V2 \
53  -d $DEBUG
54 
55 fi
56 
57 if (( 1 )); then
58 
59  debug "Generating random background."
60 
61  if ( ! reuse_file $WORKDIR/timeslice.root ); then
62 
63  JRandomTimesliceWriter \
64  -a $DETECTOR \
65  -o $WORKDIR/timeslice.root \
66  -B "$BACKGROUND_HZ" \
67  -n $numberOfSlices \
68  -P "pmt= 101 0 QE=$QE" \
69  -P "pmt= 101 1 QE=$QE" \
70  -P "pmt= 101 2 QE=$QE" \
71  -P "pmt= 101 1 QE=$QE" \
72  -P "pmt= 101 3 QE=$QE" \
73  -P "pmt= 101 4 QE=$QE" \
74  -P "pmt= 101 5 QE=$QE" \
75  -P "pmt= 101 6 QE=$QE" \
76  -P "pmt= 101 7 QE=$QE" \
77  -P "pmt= 101 8 QE=$QE" \
78  -P "pmt= 101 9 QE=$QE" \
79  -P "pmt= 101 10 QE=$QE" \
80  -P "pmt= 101 11 QE=$QE" \
81  -P "pmt= 101 12 QE=$QE" \
82  -P "pmt= 101 11 QE=$QE" \
83  -P "pmt= 101 13 QE=$QE" \
84  -P "pmt= 101 14 QE=$QE" \
85  -P "pmt= 101 15 QE=$QE" \
86  -P "pmt= 101 16 QE=$QE" \
87  -P "pmt= 101 17 QE=$QE" \
88  -P "pmt= 101 18 QE=$QE" \
89  -P "pmt= 101 19 QE=$QE" \
90  -P "pmt= 101 20 QE=$QE" \
91  -P "pmt= 101 21 QE=$QE" \
92  -P "pmt= 101 22 QE=$QE" \
93  -P "pmt= 101 21 QE=$QE" \
94  -P "pmt= 101 23 QE=$QE" \
95  -P "pmt= 101 24 QE=$QE" \
96  -P "pmt= 101 25 QE=$QE" \
97  -P "pmt= 101 26 QE=$QE" \
98  -P "pmt= 101 27 QE=$QE" \
99  -P "pmt= 101 28 QE=$QE" \
100  -P "pmt= 101 29 QE=$QE" \
101  -P "pmt= 101 30 QE=$QE" \
102  -d $DEBUG --!
103 
104  fi
105 fi
106 
107 if (( 1 )); then
108 
109  debug "Processing data."
110 
111  if ( ! reuse_file $WORKDIR/trigger_processor.root ); then
112 
113  JTriggerProcessor \
114  -a $DETECTOR \
115  -f $WORKDIR/timeslice.root \
116  -o $WORKDIR/trigger_processor.root \
117  -@"ctMin=-1.0" \
118  -@"TMaxLocal_ns=50" \
119  -@"writeL1=1" \
120  -d $DEBUG
121 
122  fi
123 fi
124 
125 if (( 1 )); then
126 
127  debug "Monitoring data."
128 
129  JCalibrateK40 \
130  -f $WORKDIR/trigger_processor.root \
131  -C JDAQTimesliceL1 \
132  -a $DETECTOR \
133  -o $WORKDIR/monitor.root \
134  -d $DEBUG
135 
136  debug "Merging data."
137 
138  JMergeCalibrateK40 \
139  -f $WORKDIR/monitor.root \
140  -o $WORKDIR/merge.root \
141  -d $DEBUG
142 
143  debug "Fitting data."
144  debug "Detector: $DETECTOR"
145 
146  JFitK40 \
147  -a $DETECTOR \
148  -f $WORKDIR/merge.root \
149  -o $WORKDIR/fit.root \
150  -w \
151  -d $DEBUG
152 
153 fi
154 
155 if (( 1 )); then
156 
157  debug "Plotting results"
158 
159  JPlot2D \
160  -f $WORKDIR/monitor.root:\.\*2S \
161  -OLEGO \
162  -\> "PMT pair" \
163  -\< "#DeltaT [ns]" \
164  -T "" \
165  -o monitor.$FORMAT $BATCH
166 
167 fi
168 
169 if (( 1 )); then
170 
171  debug "Plotting standard deviations"
172 
173  JPlot2D \
174  -f $WORKDIR/fit.root:\.\*2F \
175  -OLEGO \
176  -\> "PMT pair" \
177  -\< "#DeltaT [ns]" \
178  -z "-3 +3" \
179  -T "standard deviation" \
180  -o stdev.$FORMAT $BATCH
181 
182 fi
183 
184 if (( 1 )); then
185 
186  debug "Plotting results for t0"
187 
188  JPlot1D \
189  -f $WORKDIR/fit.root:\.\*1t0 \
190  -y "-10 +10" \
191  -\> "PMT" \
192  -\^ "#DeltaT [ns]" \
193  -L TR \
194  -T "time calibration" \
195  -o t0.$FORMAT $BATCH
196 
197 fi
198 
199 if (( 1 )); then
200 
201  debug "Plotting results for sigma"
202 
203  JPlot1D \
204  -f $WORKDIR/fit.root:\.\*1sigma \
205  -y "0 5" \
206  -\> "PMT" \
207  -\^ "#sigma [ns]" \
208  -L TR \
209  -T "time resolution" \
210  -o sigma.$FORMAT $BATCH
211 
212 fi
213 
214 
215 if (( 1 )); then
216 
217  debug "Plotting results for QE"
218 
219  JPlot1D \
220  -f $WORKDIR/fit.root:\.\*1QE \
221  -y "0 2" \
222  -\> "PMT" \
223  -\^ "QE" \
224  -L TR \
225  -T "QE" \
226  -o QE.$FORMAT $BATCH
227 
228 fi
detectorType
remove to overwrite exit fi JDetector D $ID detectorType
Definition: software/JDataFilter/JDataFilter.sh:114
JTOOLS::w
data_type w[N+1][M+1]
Definition: JPolint.hh:778
version
version
Definition: JEditTuneHV.sh:5
JCALIBRATE::JFitK40
JFitK40_t< TF2 > JFitK40
Type definition for backward compatibility.
Definition: JFitK40.hh:807
f
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
TMPDIR
usr share Modules init zsh export TMPDIR
Definition: launch_pre-calibration_D1ORCA011.sh:27
exit
exit
Definition: JPizza.sh:36
V
V(JDAQEvent-JTriggerReprocessor)*1.0/(JDAQEvent+1.0e-10)
o
then fatal Wrong number of arguments fi JConvertDetectorFormat a o
Definition: compare-detector.sh:31
GAUSS_HERMITE::sigma
const double sigma[]
Definition: JQuadrature.cc:74
echo
then echo
Definition: JDOMDAQDriver.sh:56
JPHYSICS::C
static const double C
Physics constants.
Definition: JPhysics/JConstants.hh:25
JTOOLS::n
const int n
Definition: JPolint.hh:697
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
a
then JCalibrateToT a
Definition: JTuneHV.sh:116
T
do set_variable OUTPUT_DIRECTORY $WORKDIR T
Definition: JCalibrateHeight.sh:61
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
F
then awk F
Definition: getMIIseacurrent.sh:86
usage
* usage
Definition: JLigier-local.sh:34
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
monitor
static void monitor()
Definition: DataQueue.cc:32
DETECTOR
set_variable DETECTOR
Definition: JARCACalibrateMuon.sh:31
detector
do set_variable DETECTOR_TXT $WORKDIR detector
Definition: detector-XY:fit1d.sh:38
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
D
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
debug
int debug
debug level
Definition: archive-put-wiki-detectors.sh:92
P
then $DIR JPlotNPE PDG P
Definition: JPlotNPE-PDG.sh:62
DEBUG
#define DEBUG(A)
Message macros.
Definition: JMessage.hh:62
Generated by   doxygen 1.8.5