Belle II Software development
trgrawdataModule.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9#include <trg/trg/modules/trgrawdata/trgrawdataModule.h>
10#include <trg/gdl/dbobjects/TRGGDLDBUnpacker.h>
11#include <framework/dataobjects/EventMetaData.h>
12#include <framework/datastore/StoreArray.h>
13#include <rawdata/dataobjects/RawTRG.h>
14
15#include <string.h>
16
17using namespace std;
18using namespace Belle2;
19
20
21REG_MODULE(TRGRAWDATA);
22
24{
25
26 setDescription("TRG rawdata analyser");
28
29 addParam("histRecord", m_histRecord,
30 "Recording event by event timing distribution histogram or not",
31 false);
32 addParam("mydebug", m_mydebug, "debug flag", 0);
33 addParam("cc_check", m_cc_check, "cc mode", true);
34 addParam("print_cc", m_print_cc, "Print Disordered clocks", true);
35 addParam("print_clkcyc_err", m_print_clkcyc_err, "Print Errs of Disordered clocks", false);
36 addParam("print_dbmap", m_print_dbmap, "Print Database Bit Map", false);
37
38 addParam("hdr_nwd_2d0", m_hdr_nwd_2d0, "Number of word for 2D0 Header", 3);
39 addParam("hdr_nwd_2d1", m_hdr_nwd_2d1, "Number of word for 2D1 Header", 3);
40 addParam("hdr_nwd_2d2", m_hdr_nwd_2d2, "Number of word for 2D2 Header", 3);
41 addParam("hdr_nwd_2d3", m_hdr_nwd_2d3, "Number of word for 2D3 Header", 3);
42 addParam("hdr_nwd_3d0", m_hdr_nwd_3d0, "Number of word for 3D0 Header", 3);
43 addParam("hdr_nwd_3d1", m_hdr_nwd_3d1, "Number of word for 3D1 Header", 3);
44 addParam("hdr_nwd_3d2", m_hdr_nwd_3d2, "Number of word for 3D2 Header", 3);
45 addParam("hdr_nwd_3d3", m_hdr_nwd_3d3, "Number of word for 3D3 Header", 3);
46 addParam("hdr_nwd_nn0", m_hdr_nwd_nn0, "Number of word for NN0 Header", 3);
47 addParam("hdr_nwd_nn1", m_hdr_nwd_nn1, "Number of word for NN1 Header", 3);
48 addParam("hdr_nwd_nn2", m_hdr_nwd_nn2, "Number of word for NN2 Header", 3);
49 addParam("hdr_nwd_nn3", m_hdr_nwd_nn3, "Number of word for NN3 Header", 3);
50 addParam("hdr_nwd_sl0", m_hdr_nwd_sl0, "Number of word for TSF0 Header", 3);
51 addParam("hdr_nwd_sl1", m_hdr_nwd_sl1, "Number of word for TSF1 Header", 3);
52 addParam("hdr_nwd_sl2", m_hdr_nwd_sl2, "Number of word for TSF2 Header", 3);
53 addParam("hdr_nwd_sl3", m_hdr_nwd_sl3, "Number of word for TSF3 Header", 3);
54 addParam("hdr_nwd_sl4", m_hdr_nwd_sl4, "Number of word for TSF4 Header", 3);
55 addParam("hdr_nwd_sl5", m_hdr_nwd_sl5, "Number of word for TSF5 Header", 3);
56 addParam("hdr_nwd_sl6", m_hdr_nwd_sl6, "Number of word for TSF6 Header", 3);
57 addParam("hdr_nwd_sl8", m_hdr_nwd_sl8, "Number of word for TSF8 Header", 3);
58 addParam("hdr_nwd_etf", m_hdr_nwd_etf, "Number of word for ETF Header", 3);
59 addParam("hdr_nwd_grl", m_hdr_nwd_grl, "Number of word for GRL Header", 3);
60 addParam("hdr_nwd_top", m_hdr_nwd_top, "Number of word for TOP Header", 3);
61 addParam("hdr_nwd_gdl", m_hdr_nwd_gdl, "Number of word for GDL Header", 6);
62
63 addParam("cpr_2d0", m_cpr_2d0, "Copper number of 2D0", 0x11000001U);
64 addParam("cpr_2d1", m_cpr_2d1, "Copper number of 2D1", 0x11000001U);
65 addParam("cpr_2d2", m_cpr_2d2, "Copper number of 2D2", 0x11000002U);
66 addParam("cpr_2d3", m_cpr_2d3, "Copper number of 2D3", 0x11000002U);
67 addParam("cpr_3d0", m_cpr_3d0, "Copper number of 3D0", 0x11000003U);
68 addParam("cpr_3d1", m_cpr_3d1, "Copper number of 3D1", 0x11000003U);
69 addParam("cpr_3d2", m_cpr_3d2, "Copper number of 3D2", 0x11000004U);
70 addParam("cpr_3d3", m_cpr_3d3, "Copper number of 3D3", 0x11000004U);
71 addParam("cpr_nn0", m_cpr_nn0, "Copper number of NN0", 0x11000005U);
72 addParam("cpr_nn1", m_cpr_nn1, "Copper number of NN1", 0x11000005U);
73 addParam("cpr_nn2", m_cpr_nn2, "Copper number of NN2", 0x11000006U);
74 addParam("cpr_nn3", m_cpr_nn3, "Copper number of NN3", 0x11000006U);
75 addParam("cpr_sl0", m_cpr_sl0, "Copper number of SL0", 0x11000007U);
76 addParam("cpr_sl1", m_cpr_sl1, "Copper number of SL1", 0x11000009U);//*
77 addParam("cpr_sl2", m_cpr_sl2, "Copper number of SL2", 0x11000008U);
78 addParam("cpr_sl3", m_cpr_sl3, "Copper number of SL3", 0x11000008U);
79 addParam("cpr_sl4", m_cpr_sl4, "Copper number of SL4", 0x11000007U);//*
80 addParam("cpr_sl5", m_cpr_sl5, "Copper number of SL5", 0x11000009U);
81 addParam("cpr_sl6", m_cpr_sl6, "Copper number of SL6", 0x1100000aU);
82 addParam("cpr_sl7", m_cpr_sl7, "Copper number of SL7", 0x1100000a);
83 addParam("cpr_sl8", m_cpr_sl8, "Copper number of SL8", 0x1100000aU);
84 addParam("cpr_gdl", m_cpr_gdl, "Copper number of GDL", 0x15000001U);
85 addParam("cpr_etf", m_cpr_etf, "Copper number of ETF", 0x15000001U);
86 addParam("cpr_grl", m_cpr_grl, "Copper number of GRL", 0x15000002);
87 addParam("cpr_top", m_cpr_top, "Copper number of TOP", 0x12000001U);
88
89 addParam("on_2d0", m_on_2d0, "ON/OFF of 2D0", true);
90 addParam("on_2d1", m_on_2d1, "ON/OFF of 2D1", true);
91 addParam("on_2d2", m_on_2d2, "ON/OFF of 2D2", true);
92 addParam("on_2d3", m_on_2d3, "ON/OFF of 2D3", true);
93 addParam("on_3d0", m_on_3d0, "ON/OFF of 3D0", true);
94 addParam("on_3d1", m_on_3d1, "ON/OFF of 3D1", true);
95 addParam("on_3d2", m_on_3d2, "ON/OFF of 3D2", true);
96 addParam("on_3d3", m_on_3d3, "ON/OFF of 3D3", true);
97 addParam("on_nn0", m_on_nn0, "ON/OFF of NN0", true);
98 addParam("on_nn1", m_on_nn1, "ON/OFF of NN1", true);
99 addParam("on_nn2", m_on_nn2, "ON/OFF of NN2", true);
100 addParam("on_nn3", m_on_nn3, "ON/OFF of NN3", true);
101 addParam("on_sl0", m_on_sl0, "ON/OFF of SL0", true);
102 addParam("on_sl1", m_on_sl1, "ON/OFF of SL1", true);
103 addParam("on_sl2", m_on_sl2, "ON/OFF of SL2", true);
104 addParam("on_sl3", m_on_sl3, "ON/OFF of SL3", true);
105 addParam("on_sl4", m_on_sl4, "ON/OFF of SL4", true);
106 addParam("on_sl5", m_on_sl5, "ON/OFF of SL5", true);
107 addParam("on_sl6", m_on_sl6, "ON/OFF of SL6", true);
108 addParam("on_sl7", m_on_sl7, "ON/OFF of SL7", true);
109 addParam("on_sl8", m_on_sl8, "ON/OFF of SL8", true);
110 addParam("on_gdl", m_on_gdl, "ON/OFF of GDL", true);
111 addParam("on_etf", m_on_etf, "ON/OFF of ETF", true);
112 addParam("on_grl", m_on_grl, "ON/OFF of GRL", true);
113 addParam("on_top", m_on_top, "ON/OFF of TOP", true);
114
115 addParam("nwd_2ds", m_nwd_2ds, "N word of 2Ds", 3939); // 2624x48
116 addParam("nwd_3ds", m_nwd_3ds, "N word of 3Ds", 3939); // 2624x48
117 addParam("nwd_nns", m_nwd_nns, "N word of NNs", 0xc03); // 2048x48
118 addParam("nwd_sl0", m_nwd_sl0, "N word of SL0", 0xc03); // 2048x48
119 addParam("nwd_sl1", m_nwd_sl1, "N word of SL1", 0x1803); // 4096x48
120 addParam("nwd_sl2", m_nwd_sl2, "N word of SL2", 0x1803); // 4096x48
121 addParam("nwd_sl3", m_nwd_sl3, "N word of SL3", 0x1803); // 4096x48
122 addParam("nwd_sl4", m_nwd_sl4, "N word of SL4", 0xc03); // 2048x48
123 addParam("nwd_sl5", m_nwd_sl5, "N word of SL5", 0x1803); // 4096x48
124 addParam("nwd_sl6", m_nwd_sl6, "N word of SL6", 0xc03); // 2048x48
125 addParam("nwd_sl7", m_nwd_sl7, "N word of SL7", 0x603); // 1024x48
126 addParam("nwd_sl8", m_nwd_sl8, "N word of SL8", 0x603); // 1024x48
127 addParam("nwd_gdl", m_nwd_gdl, "N word of GDL", 646); // 640x32
128 addParam("nwd_etf", m_nwd_etf, "N word of ETF", 1539); // 1024x48
129 addParam("nwd_grl", m_nwd_grl, "N word of GRL", 0xc03); // 1024x48
130 addParam("nwd_top", m_nwd_top, "N word of TOP", 0x303); // 1024*24
131
132 addParam("nclk_2ds", m_nclk_2ds, "N clk of 2Ds", 48);
133 addParam("nclk_3ds", m_nclk_3ds, "N clk of 3Ds", 48);
134 addParam("nclk_nns", m_nclk_nns, "N clk of NNs", 48);
135 addParam("nclk_sl0", m_nclk_sl0, "N clk of SL0", 48);
136 addParam("nclk_sl1", m_nclk_sl1, "N clk of SL1", 48);
137 addParam("nclk_sl2", m_nclk_sl2, "N clk of SL2", 48);
138 addParam("nclk_sl3", m_nclk_sl3, "N clk of SL3", 48);
139 addParam("nclk_sl4", m_nclk_sl4, "N clk of SL4", 48);
140 addParam("nclk_sl5", m_nclk_sl5, "N clk of SL5", 48);
141 addParam("nclk_sl6", m_nclk_sl6, "N clk of SL6", 48);
142 addParam("nclk_sl7", m_nclk_sl7, "N clk of SL7", 48);
143 addParam("nclk_sl8", m_nclk_sl8, "N clk of SL8", 48);
144 addParam("nclk_gdl", m_nclk_gdl, "N clk of GDL", 32);
145 addParam("nclk_etf", m_nclk_etf, "N clk of ETF", 48);
146 addParam("nclk_grl", m_nclk_grl, "N clk of GRL", 48);
147 addParam("nclk_top", m_nclk_top, "N clk of TOP", 24);
148
149 addParam("nclk_cdctrg", m_nclk_cdctrg, "N clock of CDCTRG", 48);
150
151 addParam("scale_2ds", m_scale_2ds, "Scale factor of 2Ds", 256);
152 addParam("scale_3ds", m_scale_3ds, "Scale factor of 3Ds", 256);
153 addParam("scale_nns", m_scale_nns, "Scale factor of NNs", 256);
154 addParam("scale_sl0", m_scale_sl0, "Scale factor of SL0", 256);
155 addParam("scale_sl1", m_scale_sl1, "Scale factor of SL1", 256);
156 addParam("scale_sl2", m_scale_sl2, "Scale factor of SL2", 256);
157 addParam("scale_sl3", m_scale_sl3, "Scale factor of SL3", 256);
158 addParam("scale_sl4", m_scale_sl4, "Scale factor of SL4", 256);
159 addParam("scale_sl5", m_scale_sl5, "Scale factor of SL5", 256);
160 addParam("scale_sl6", m_scale_sl6, "Scale factor of SL6", 256);
161 addParam("scale_sl7", m_scale_sl7, "Scale factor of SL7", 256);
162 addParam("scale_sl8", m_scale_sl8, "Scale factor of SL8", 256);
163 addParam("scale_gdl", m_scale_gdl, "Scale factor of GDL", 1);
164 addParam("scale_etf", m_scale_etf, "Scale factor of ETF", 512);
165 addParam("scale_grl", m_scale_grl, "Scale factor of GRL", 1);
166 addParam("scale_top", m_scale_top, "Scale factor of TOP", 32);
167
168 addParam("fmid_2d0", m_fmid_2d0, "Firmware ID of 2D0", 0x32442020);
169 addParam("fmid_2d1", m_fmid_2d1, "Firmware ID of 2D1", 0x32442020);
170 addParam("fmid_2d2", m_fmid_2d2, "Firmware ID of 2D2", 0x32442020);
171 addParam("fmid_2d3", m_fmid_2d3, "Firmware ID of 2D3", 0x32442020);
172 addParam("fmid_3d0", m_fmid_3d0, "Firmware ID of 3D0", 0x33440001);
173 addParam("fmid_3d1", m_fmid_3d1, "Firmware ID of 3D1", 0x33440001);
174 addParam("fmid_3d2", m_fmid_3d2, "Firmware ID of 3D2", 0x33440001);
175 addParam("fmid_3d3", m_fmid_3d3, "Firmware ID of 3D3", 0x33440001);
176 addParam("fmid_nn0", m_fmid_nn0, "Firmware ID of NN0", 0x4e4e5402);
177 addParam("fmid_nn1", m_fmid_nn1, "Firmware ID of NN1", 0x4e4e5402);
178 addParam("fmid_nn2", m_fmid_nn2, "Firmware ID of NN2", 0x4e4e5402);
179 addParam("fmid_nn3", m_fmid_nn3, "Firmware ID of NN3", 0x4e4e5402);
180 addParam("fmid_sl0", m_fmid_sl0, "Firmware ID of SL0", 0x54534630);
181 addParam("fmid_sl1", m_fmid_sl1, "Firmware ID of SL1", 0x54534631);
182 addParam("fmid_sl2", m_fmid_sl2, "Firmware ID of SL2", 0x54534632);
183 addParam("fmid_sl3", m_fmid_sl3, "Firmware ID of SL3", 0x54534633);
184 addParam("fmid_sl4", m_fmid_sl4, "Firmware ID of SL4", 0x54534634);
185 addParam("fmid_sl5", m_fmid_sl5, "Firmware ID of SL5", 0x54534635);
186 addParam("fmid_sl6", m_fmid_sl6, "Firmware ID of SL6", 0x54534636);
187 addParam("fmid_sl7", m_fmid_sl7, "Firmware ID of SL7", 0x54534637);
188 addParam("fmid_sl8", m_fmid_sl8, "Firmware ID of SL8", 0x54534d38);
189 addParam("fmid_gdl", m_fmid_gdl, "Firmware ID of GDL", 0x47444c20U);
190 addParam("fmid_etf", m_fmid_etf, "Firmware ID of ETF", 0x45544620);
191 addParam("fmid_grl", m_fmid_grl, "Firmware ID of GRL", 0x47824c20);
192 addParam("fmid_top", m_fmid_top, "Firmware ID of TOP", 0);
193
194}
195
197{
198 B2INFO("trgrawdata: initialize() started.");
199 if (! m_unpacker) {
200 B2INFO("no database of gdl unpacker");
201 }
202 if (m_print_dbmap) {
203
204 //const int nword_header = m_unpacker->get_nword_header(); //not used for now
205 //const int n_clocks = m_unpacker->getnClks(); //not used for now
206 //const int nBits = m_unpacker->getnBits(); //not used for now
207 const int n_leafs = m_unpacker->getnLeafs();
208 const int n_leafsExtra = m_unpacker->getnLeafsExtra();
209 //const int conf = m_unpacker->getconf(); //not used for now
210
211 int aBitMap[320][2] = {{0}};
212 int aBitMap_extra[100][3] = {{ -1}};
213 for (int i = 0; i < n_leafsExtra; i++) {
214 aBitMap_extra[i][0] = m_unpacker->getBitMap_extra(i, 0);
215 aBitMap_extra[i][1] = m_unpacker->getBitMap_extra(i, 1);
216 aBitMap_extra[i][2] = m_unpacker->getBitMap_extra(i, 2);
217 }
218
219 char LeafNames[200][100] = {{0}};
220 for (int i = 0; i < 200; i++) {
221 std::cout << "LeafBitMap[" << i << "] = " << m_unpacker->getLeafMap(i) << std::endl;
222 strcpy(LeafNames[i], m_unpacker->getLeafnames(i));
223 std::cout << "LeafNames[" << i << "] = " << m_unpacker->getLeafnames(i) << std::endl;
224 }
225 for (int i = 0; i < n_leafs; i++) {
226 aBitMap[i][0] = m_unpacker->getBitMap(i, 0);
227 aBitMap[i][1] = m_unpacker->getBitMap(i, 1);
228 }
229 for (int i = 0; i < 200; i++) {
230 int bin = m_unpacker->getLeafMap(i) + 1;
231 if (0 < bin && bin <= n_leafs) {
232 std::cout << "leaf(" << i
233 << "), bin(" << bin
234 << "), LeafNames[leaf](" << LeafNames[i]
235 << "), BitMap[bin-1][0](" << aBitMap[bin - 1][0]
236 << "), BitMap[bin-1][1](" << aBitMap[bin - 1][1]
237 << ")" << std::endl;
238 }
239 }
240 // for leafsExtra
241 for (int i = 0; i < 200; i++) {
242 int bin = m_unpacker->getLeafMap(i) + 1;
243 int j = bin - n_leafs - 1;
244 if (n_leafs < bin && bin <= n_leafs + n_leafsExtra) {
245 std::cout << "i(" << i
246 << "), bin(" << bin
247 << "), LeafNames[leaf](" << LeafNames[i]
248 << "), BitMap_extra[j][0](buf[" << aBitMap_extra[j][0]
249 << "]), BitMap_extra[j][1](downto " << aBitMap_extra[j][1]
250 << "), BitMap_extra[j][2](" << aBitMap_extra[j][2]
251 << " bit length)" << std::endl;
252 }
253 }
254
255 }
256
257}
258
260{
261
262 B2INFO("trgrawdata: event() started.");
263 StoreArray<RawTRG> raw_trgarray;
264
266 int _exp = bevt->getExperiment();
267 int _run = bevt->getRun();
268 unsigned _eve = bevt->getEvent();
269 if (0)
270 std::cout << "**************_eve(" << _eve
271 << "), c_nModules(" << c_nModules
272 << "), m_nclk_cdctrg(" << m_nclk_cdctrg
273 << "), _exp(" << _exp
274 << "), _run(" << _run
275 << "), _eve(" << _eve
276 << "), _Form(" << Form("he%010d", (int)_eve)
277 << ")" << std::endl;
278
279 TH2D* h_0 = new TH2D(Form("he%010d", (int)_eve),
280 Form("e%02dr%04d ev %08d", _exp, _run, (int)_eve),
281 c_nModules * 4, 0, c_nModules * 4,
282 m_nclk_cdctrg, 0, m_nclk_cdctrg);
283
284 for (unsigned i = 0; i < c_nModules; ++i) {
285 h_0->GetXaxis()->SetBinLabel(4 * i + 1, Form("%s hdevt", moduleNames[i]));
286 h_0->GetXaxis()->SetBinLabel(4 * i + 2, Form("%s dtevt", moduleNames[i]));
287 h_0->GetXaxis()->SetBinLabel(4 * i + 3, Form("%s L1rvc", moduleNames[i]));
288 h_0->GetXaxis()->SetBinLabel(4 * i + 4, Form("%s cc", moduleNames[i]));
289 }
290
291 /*
292 GDL only. To get gdlrvc12
293 */
294 unsigned gdlrvc12 = -1;
295 for (int i = 0; i < raw_trgarray.getEntries(); i++) {
296 for (int j = 0; j < raw_trgarray[i]->GetNumEntries(); j++) {
297 int cprid = (int)raw_trgarray[i]->GetNodeID(j);
298 unsigned eve20 = (_eve & 0xfffff);
299 for (int hslb = 0; hslb < 2; hslb++) {
300 int nword = raw_trgarray[i]->GetDetectorNwords(j, hslb);
301 int* buf = raw_trgarray[i]->GetDetectorBuffer(j, hslb);
302 if (0x15000001 == cprid && hslb == 0) {
303
304 if (m_mydebug) {
305 printf("aa:GDL start: 0x%x%c exp(%d), run(%d), eve(%u), eve(0x%x), nword(%d)\n",
306 cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword);
307 }
308
309 unsigned buf0 = 0;
310 if (nword > 0) buf0 = (unsigned)buf[0];
311 unsigned buf1 = 0;
312 if (nword > 1) buf1 = (unsigned)buf[1];
313 unsigned buf2 = 0;
314 if (nword > 2) buf2 = (unsigned)buf[2];
315 //unsigned buf2cnttrg20 = (buf2 >> 12) & 0xfffff;
316 //unsigned buf2rvc12 = (buf2 & 0xfff);
317
318 if (nword < 3) {
319 m_cntr_nwe_badnwd[e_gdl]++;
320 printf("ab:GDL(0x%x%c) exp(%d), run(%d), eve(%u), eve(0x%x), nword(%d), ",
321 cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword);
322 printf("buf0(0x%x), buf1(0x%x), buf2(0x%x)\n", buf0, buf1, buf2);
323 continue;
324 }
325
326 if (nword < 8) {
327 printf("ad:GDL(0x%x%c) exp(%d), run(%d), eve(%u), eve(0x%x), nword(%d)\n",
328 cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword);
329 m_cntr_nwe_badnwd[e_gdl]++;
330 } else {
331 unsigned buf5 = (unsigned)buf[m_hdr_nwd_gdl - 1];
332 // 626 = 646 -
333 unsigned buf6 = (unsigned)buf[626]; // 1st wd in the last clock cycle
334 unsigned buf7 = (unsigned)buf[627]; // 2nd wd in the last clock cycle
335 unsigned lastdd = (buf6 >> 24) & 0xff;
336 unsigned last_rvc = (buf6 >> 4) & 0xfff;
337 unsigned cnttrg_data_15to12 = buf6 & (0xf);
338 unsigned cnttrg_data_12 = (buf7 >> 20) & (0xfff);
339 gdlrvc12 = (buf5 & 0xfff);
340 unsigned gdlhdrcnttrg20 = (buf5 >> 12) & 0xfffff;
341
342 //unsigned wdperclk = (m_nwd_gdl - m_hdr_nwd_gdl) / m_nclk_gdl;
343 //unsigned ibuf_lastclk_msw = m_nwd_gdl - wdperclk;
344 //unsigned buf_lastclk_msw = (unsigned)buf[ibuf_lastclk_msw];
345 //unsigned buf_lastclk_msw2 = (unsigned)buf[ibuf_lastclk_msw + 1];
346 //unsigned buf_lastclk_msw3 = (unsigned)buf[ibuf_lastclk_msw + 2];
347 unsigned cnttrg_data_16 = (cnttrg_data_15to12 << 12) + cnttrg_data_12;
348 unsigned eve16 = (_eve & 0xffff);
349
350 h_0->SetBinContent(4 * (int)e_gdl + 1, m_nclk_gdl, gdlhdrcnttrg20);
351 h_0->SetBinContent(4 * (int)e_gdl + 2, m_nclk_gdl, cnttrg_data_16);
352 h_0->SetBinContent(4 * (int)e_gdl + 3, m_nclk_gdl, gdlrvc12);
353
354 if (nword != m_nwd_gdl || buf0 != m_fmid_gdl || lastdd != 0xdd ||
355 gdlhdrcnttrg20 != eve20 ||
356 !(eve16 == cnttrg_data_16 + 1 || (eve16 == 0 && cnttrg_data_16 == 0xffff))
357 ) {
358
359 unsigned diag = 0;
360 diag |= (nword != m_nwd_gdl) ? 1 : 0;
361 diag |= (buf0 != m_fmid_gdl) ? 2 : 0;
362 diag |= (lastdd != 0xdd) ? 4 : 0;
363 diag |= (gdlhdrcnttrg20 != eve20) ? 8 : 0;
364 diag |= !(eve16 == cnttrg_data_16 + 1 || (eve16 == 0 && cnttrg_data_16 == 0xffff)) ? 16 : 0;
365 printf("ae:GDL(0x%x%c) exp(%d), run(%d), evedaq(%u=0x%x), gdlhdrcnttrg20(%u,0x%x), "
366 "cnttrg_data16(%u=0x%x), l1rvc(%u)-lastrvc(%u)=%u, nword(%d), diag(%u)\n",
367 cprid, 'a' + hslb, _exp, _run, _eve, _eve, gdlhdrcnttrg20, gdlhdrcnttrg20,
368 cnttrg_data_16, cnttrg_data_16, gdlrvc12, last_rvc, gdlrvc12 - last_rvc, nword, diag); //TODO can the difference be negative?
369 printf("\n");
370 m_cntr_nwn_badtrg[e_gdl]++;
371 } else {
372 m_cntr_nwn[e_gdl]++;
373 }
374 }
375 }
376 }
377 }
378 }
379
380 /* Other than GDL */
381 for (int i = 0; i < raw_trgarray.getEntries(); i++) {
382 for (int j = 0; j < raw_trgarray[i]->GetNumEntries(); j++) {
383 unsigned cprid = raw_trgarray[i]->GetNodeID(j);
384 //unsigned eve16 = (_eve & 0xffff);
385 unsigned eve20 = (_eve & 0xfffff);
386 for (int hslb = 0; hslb < 2; hslb++) {
387 unsigned e_mdl = -1;
388 unsigned _nwd = -1;
389 unsigned _scale = -1;
390 unsigned _nclk = -1;
391 unsigned _hdr_nwd = 3;
392 if (m_cpr_2d0 == cprid && hslb == 0 && m_on_2d0) {
393 e_mdl = e_2d0; _nwd = m_nwd_2ds; _scale = m_scale_2ds; _nclk = m_nclk_2ds; _hdr_nwd = m_hdr_nwd_2d0;
394 } else if (m_cpr_2d1 == cprid && hslb == 1 && m_on_2d1) {
395 e_mdl = e_2d1; _nwd = m_nwd_2ds; _scale = m_scale_2ds; _nclk = m_nclk_2ds; _hdr_nwd = m_hdr_nwd_2d1;
396 } else if (m_cpr_2d2 == cprid && hslb == 0 && m_on_2d2) {
397 e_mdl = e_2d2; _nwd = m_nwd_2ds; _scale = m_scale_2ds; _nclk = m_nclk_2ds; _hdr_nwd = m_hdr_nwd_2d2;
398 } else if (m_cpr_2d3 == cprid && hslb == 1 && m_on_2d3) {
399 e_mdl = e_2d3; _nwd = m_nwd_2ds; _scale = m_scale_2ds; _nclk = m_nclk_2ds; _hdr_nwd = m_hdr_nwd_2d3;
400 } else if (m_cpr_3d0 == cprid && hslb == 0 && m_on_3d0) {
401 e_mdl = e_3d0; _nwd = m_nwd_3ds; _scale = m_scale_3ds; _nclk = m_nclk_3ds; _hdr_nwd = m_hdr_nwd_3d0;
402 } else if (m_cpr_3d1 == cprid && hslb == 1 && m_on_3d1) {
403 e_mdl = e_3d1; _nwd = m_nwd_3ds; _scale = m_scale_3ds; _nclk = m_nclk_3ds; _hdr_nwd = m_hdr_nwd_3d1;
404 } else if (m_cpr_3d2 == cprid && hslb == 0 && m_on_3d2) {
405 e_mdl = e_3d2; _nwd = m_nwd_3ds; _scale = m_scale_3ds; _nclk = m_nclk_3ds; _hdr_nwd = m_hdr_nwd_3d2;
406 } else if (m_cpr_3d3 == cprid && hslb == 1 && m_on_3d3) {
407 e_mdl = e_3d3; _nwd = m_nwd_3ds; _scale = m_scale_3ds; _nclk = m_nclk_3ds; _hdr_nwd = m_hdr_nwd_3d3;
408 } else if (m_cpr_nn0 == cprid && hslb == 0 && m_on_nn0) {
409 e_mdl = e_nn0; _nwd = m_nwd_nns; _scale = m_scale_nns; _nclk = m_nclk_nns; _hdr_nwd = m_hdr_nwd_nn0;
410 } else if (m_cpr_nn1 == cprid && hslb == 1 && m_on_nn1) {
411 e_mdl = e_nn1; _nwd = m_nwd_nns; _scale = m_scale_nns; _nclk = m_nclk_nns; _hdr_nwd = m_hdr_nwd_nn1;
412 } else if (m_cpr_nn2 == cprid && hslb == 0 && m_on_nn2) {
413 e_mdl = e_nn2; _nwd = m_nwd_nns; _scale = m_scale_nns; _nclk = m_nclk_nns; _hdr_nwd = m_hdr_nwd_nn2;
414 } else if (m_cpr_nn3 == cprid && hslb == 1 && m_on_nn3) {
415 e_mdl = e_nn3; _nwd = m_nwd_nns; _scale = m_scale_nns; _nclk = m_nclk_nns; _hdr_nwd = m_hdr_nwd_nn3;
416 } else if (m_cpr_sl0 == cprid && hslb == 0 && m_on_sl0) {
417 e_mdl = e_sl0; _nwd = m_nwd_sl0; _scale = m_scale_sl0; _nclk = m_nclk_sl0; _hdr_nwd = m_hdr_nwd_sl0;
418 } else if (m_cpr_sl1 == cprid && hslb == 0 && m_on_sl1) { //ex w/ 1
419 e_mdl = e_sl1; _nwd = m_nwd_sl1; _scale = m_scale_sl1; _nclk = m_nclk_sl1; _hdr_nwd = m_hdr_nwd_sl1;
420 } else if (m_cpr_sl2 == cprid && hslb == 0 && m_on_sl2) {
421 e_mdl = e_sl2; _nwd = m_nwd_sl2; _scale = m_scale_sl2; _nclk = m_nclk_sl2; _hdr_nwd = m_hdr_nwd_sl2;
422 } else if (m_cpr_sl3 == cprid && hslb == 1 && m_on_sl3) {
423 e_mdl = e_sl3; _nwd = m_nwd_sl3; _scale = m_scale_sl3; _nclk = m_nclk_sl3; _hdr_nwd = m_hdr_nwd_sl3;
424 } else if (m_cpr_sl4 == cprid && hslb == 1 && m_on_sl4) { //ex w/ 4
425 e_mdl = e_sl4; _nwd = m_nwd_sl4; _scale = m_scale_sl4; _nclk = m_nclk_sl4; _hdr_nwd = m_hdr_nwd_sl4;
426 } else if (m_cpr_sl5 == cprid && hslb == 1 && m_on_sl5) {
427 e_mdl = e_sl5; _nwd = m_nwd_sl5; _scale = m_scale_sl5; _nclk = m_nclk_sl5; _hdr_nwd = m_hdr_nwd_sl5;
428 } else if (m_cpr_sl6 == cprid && hslb == 0 && m_on_sl6) {
429 e_mdl = e_sl6; _nwd = m_nwd_sl6; _scale = m_scale_sl6; _nclk = m_nclk_sl6; _hdr_nwd = m_hdr_nwd_sl6;
430 } else if (m_cpr_sl8 == cprid && hslb == 0 && m_on_sl8) {
431 e_mdl = e_sl8; _nwd = m_nwd_sl8; _scale = m_scale_sl8; _nclk = m_nclk_sl8; _hdr_nwd = m_hdr_nwd_sl8;
432 } else if (m_cpr_etf == cprid && hslb == 1 && m_on_etf) {
433 e_mdl = e_etf; _nwd = m_nwd_etf; _scale = m_scale_etf; _nclk = m_nclk_etf; _hdr_nwd = m_hdr_nwd_etf;
434 } else if (m_cpr_top == cprid && hslb == 0 && m_on_top) {
435 e_mdl = e_top; _nwd = m_nwd_top; _scale = m_scale_top; _nclk = m_nclk_top; _hdr_nwd = m_hdr_nwd_top;
436 } else {
437 continue;
438 }
439
440 unsigned nword = raw_trgarray[i]->GetDetectorNwords(j, hslb);
441 int* buf = raw_trgarray[i]->GetDetectorBuffer(j, hslb);
442 unsigned buf0 = 0;
443 if (nword > 0) buf0 = (unsigned)buf[0];
444 unsigned buf1 = 0;
445 if (nword > 1) buf1 = (unsigned)buf[1];
446 unsigned buf2 = 0;
447 if (nword > 2) buf2 = (unsigned)buf[2];
448 unsigned buf2cnttrg20 = (buf2 >> 12) & 0xfffff;
449 unsigned buf2rvc12 = (buf2 & 0xfff);
450
451 if (nword < 3) {
452 // err
453 m_cntr_nwe_badnwd[e_mdl]++;
454 printf("ah:%s(0x%x%c) exp(%d), run(%d), eve(%u), eve(0x%x), nword(%u), ",
455 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword);
456 printf("buf0(0x%x), buf1(0x%x), buf2(0x%x)\n", buf0, buf1, buf2);
457 continue;
458 } else if (nword == _hdr_nwd) {
459 // header-only event
460 if (gdlrvc12 != buf2rvc12) {
461 printf("ai:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), gdlrvc12(0x%x), buf2rvc12(0x%x)\n",
462 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, gdlrvc12, buf2rvc12);
463 m_cntr_nw3_badrvc[e_mdl]++;
464 } else if (eve20 != buf2cnttrg20) {
465 printf("aj:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), cnttrg20(0x%x)\n",
466 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, buf2cnttrg20);
467 m_cntr_nw3_badtrg[e_mdl]++;
468 } else if ((_eve % _scale) == 0) {
469 m_cntr_nw3_badvet[e_mdl]++;
470 } else {
471 m_cntr_nw3[e_mdl]++;
472 }
473 continue;
474 } else if (nword != _nwd) {
475 // err
476 printf("bo:wrong nword: %s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), buf0(0x%x), buf2(0x%x)\n",
477 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, buf0, buf2);
478 m_cntr_nwe_badnwd[e_mdl]++;
479 continue;
480 }
481
482 unsigned buf3 = (unsigned)buf[_hdr_nwd];
483 unsigned buf3dddd = ((buf3 >> 16) & 0xffff);
484 unsigned buf4 = (unsigned)buf[_hdr_nwd + 1];
485 unsigned buf5 = (unsigned)buf[_hdr_nwd + 2];
486 unsigned wdperclk = (_nwd - _hdr_nwd) / _nclk;
487 //unsigned ibuf_lsw = _nwd - 1;
488 unsigned ibuf_lastclk_msw = _nwd - wdperclk;
489 //unsigned buf_lastclk_msw = (unsigned)buf[ibuf_lastclk_msw];
490 unsigned buf_lastclk_msw2 = (unsigned)buf[ibuf_lastclk_msw + 1];
491 unsigned buf_lastclk_msw3 = (unsigned)buf[ibuf_lastclk_msw + 2];
492 //unsigned dddd_lastclk = ((buf_lastclk_msw >> 16) & 0xffff);
493 unsigned datacnttrg32 = 0;
494 unsigned buf4_msb16 = (buf4 & 0xffff);
495 unsigned buf5_lsb16 = ((buf5 >> 16) & 0xffff);
496
497 if (m_mydebug) {
498 printf("af:Debug0: %s(0x%x%c) exp(%d), run(%d), eve(%u), eve(0x%x), nword(%u)\n",
499 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword);
500 }
501
502 /*
503 *
504 * Module Start
505 *
506 *
507 */
508
509 if ((0x11000009 == cprid && hslb == 0) || // SL1
510 (0x11000008 == cprid && hslb == 1) || // SL1
511 (0x11000009 == cprid && hslb == 1) // SL5
512 ) {
513 /* SL1 SL3 SL5 4096x48
514 data_b2l_r(4095 downto 4048) <= x"dddd" & "00000" & revoclk & m_cntr125M(15 downto 0);--i=0..5
515 data_b2l_r(4031 downto 4000) <= cnttrg;--i=6..9
516 buf[3] = dddd * revoclk
517 buf[4] = m_cntr125M & x"0000"
518 buf[5] = cnttrg
519 */
520 datacnttrg32 = buf5;
521
522 } else if ((0x11000007 == cprid && hslb == 0) || // SL0
523 (0x11000007 == cprid && hslb == 1) || // SL4
524 (0x1100000a == cprid && hslb == 0) // SL6
525 ) {
526 /* SL0 SL4 SL6 2048*48
527 data_b2l_r(2047 downto 2000) <= x"dddd" &"00000"& revoclk & m_cntr125M(15 downto 0);
528 data_b2l_r(1999 downto 1968) <= cnttrg;
529
530 data_b2l(2047, 2016) <= buf[3]
531 data_b2l(2015, 1984) <= buf[4]
532 data_b2l(1983, 1952) <= buf[5]
533 data_b2l(1999, 1968) <= cnttrg;
534 */
535 datacnttrg32 = (buf4_msb16 << 16) + buf5_lsb16;
536
537 } else if ((m_cpr_2d0 == cprid && hslb == 0 && m_on_2d0) ||
538 (m_cpr_2d1 == cprid && hslb == 1 && m_on_2d1) ||
539 (m_cpr_2d2 == cprid && hslb == 0 && m_on_2d2) ||
540 (m_cpr_2d3 == cprid && hslb == 1 && m_on_2d3)) {
541 /*
542 2D2D2D: 2048 bit x 48 clock = 64 word x 48 clock
543 */
544
545 unsigned buf_lastclk_lsb = (unsigned)buf[_nwd - 1];
546 unsigned buf_lastclk_lsb2 = (unsigned)buf[_nwd - 2];
547 //unsigned nworking = (buf_lastclk_lsb2 >> 24) & 0xff;
548 //unsigned idbuf = (buf_lastclk_lsb2 >> 16) & 0xff;
549 //unsigned ihdr = (buf_lastclk_lsb2 >> 8) & 0xff;
550 datacnttrg32 = ((buf_lastclk_lsb2 & 0xff) << 24) + (buf_lastclk_lsb >> 8);
551
552 } else if ((m_cpr_3d0 == cprid && hslb == 0 && m_on_3d0) ||
553 (m_cpr_3d1 == cprid && hslb == 1 && m_on_3d1) ||
554 (m_cpr_3d2 == cprid && hslb == 0 && m_on_3d2) ||
555 (m_cpr_3d3 == cprid && hslb == 1 && m_on_3d3)) {
556 /*
557 3D3D3D: 2048 bit x 48 clock = 64 word x 48 clock
558 */
559
560 // data_b2l_r(2623 downto 2576) <= x"dddd" &"00000"& revoclk & m_cntr125M(15 downto 0);
561 // data_b2l_r(2575 downto 2544) <= cnttrg;
562 // [3] (2623, 2592)
563 // [4] (2591, 2560)
564 // [5] (2559, 2528)
565 unsigned datacnttrg_msb16 = (buf_lastclk_msw2 & 0xffff);
566 unsigned datacnttrg_lsb16 = ((buf_lastclk_msw3 >> 16) & 0xffff);
567 datacnttrg32 = (datacnttrg_msb16 << 16) + datacnttrg_lsb16;
568
569 } else if ((m_cpr_nn0 == cprid && hslb == 0 && m_on_nn0) ||
570 (m_cpr_nn1 == cprid && hslb == 1 && m_on_nn1) ||
571 (m_cpr_nn2 == cprid && hslb == 0 && m_on_nn2) ||
572 (m_cpr_nn3 == cprid && hslb == 1 && m_on_nn3)) {
573 /*
574 int e_mdl = -1;
575 if (m_cpr_nn0 == cprid && hslb == 0 && m_on_nn0) e_mdl = e_nn0;
576 if (m_cpr_nn1 == cprid && hslb == 1 && m_on_nn1) e_mdl = e_nn1;
577 if (m_cpr_nn2 == cprid && hslb == 0 && m_on_nn2) e_mdl = e_nn2;
578 if (m_cpr_nn3 == cprid && hslb == 1 && m_on_nn3) e_mdl = e_nn3;
579 */
580 /*
581 NNNNNN: 2048 bit x 48 clock = 64 word x 48 clock
582 data_b2l_r(1999 downto 1968) <= cnttrg;
583 */
584
585 // [3] 2047 2016
586 // [4] 2015 1984
587 // [5] 1983 1952
588 // [3 + 64] = [m_hdr_nwd + wdperclk]
589 // [3 + 65] = [m_hdr_nwd + wdperclk + 1]
590 // [3 + 66] = [m_hdr_nwd + wdperclk + 2]
591 // [3 + 64 * 47] = [m_hdr_nwd + wdperclk * (_nclk - 1)]
592 // [3 + 65 * 47] = [m_hdr_nwd + wdperclk * (_nclk - 1) + 1]
593 // [3 + 66 * 47] = [m_hdr_nwd + wdperclk * (_nclk - 1) + 2]
594 unsigned datacnttrg_msb16 = (buf_lastclk_msw2 & 0xffff);
595 unsigned datacnttrg_lsb16 = ((buf_lastclk_msw3 >> 16) & 0xffff);
596 datacnttrg32 = (datacnttrg_msb16 << 16) + datacnttrg_lsb16;
597 } else if (m_cpr_top == cprid && hslb == 0 && m_on_top) {
598 // TOPTRG
599 } // Modules
600
601 bool something_bad = false;
602 h_0->SetBinContent(4 * e_mdl + 1, _nclk, buf2cnttrg20);
603 h_0->SetBinContent(4 * e_mdl + 2, _nclk, datacnttrg32);
604 h_0->SetBinContent(4 * e_mdl + 3, _nclk, buf2rvc12);
605 if (gdlrvc12 != buf2rvc12) {
606 printf("ba:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), gdlrvc12(0x%x), buf2rvc12(0x%x)\n",
607 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, gdlrvc12, buf2rvc12);
608 m_cntr_nwn_badrvc[e_mdl]++;
609 something_bad = true;
610 }
611 if (buf3dddd == 0xbbbb) {
612 m_cntr_nwn_badbbb[e_mdl]++;
613 printf("bv:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), cnttrg20(0x%x), dddd(0x%x), datacnttrg32(0x%x)\n",
614 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, buf2cnttrg20, buf3dddd, datacnttrg32);
615 something_bad = true;
616 }
617 if (buf3dddd != 0xdddd) {
618 m_cntr_nwn_badddd[e_mdl]++;
619 printf("bb:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), cnttrg20(0x%x), dddd(0x%x), datacnttrg32(0x%x)\n",
620 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, buf2cnttrg20, buf3dddd, datacnttrg32);
621 something_bad = true;
622 }
623 if (datacnttrg32 + 1 != _eve) {
624 if (!(
625 (m_cpr_top == cprid) ||
626 (m_cpr_3d0 <= cprid && cprid <= m_cpr_3d3)
627 )
628 ) {
629 printf("bc:%s(0x%x%c) exp(%d), run(%d), eve(%u,0x%x), nword(%u), cnttrg20(0x%x), dddd(0x%x), datacnttrg32(0x%x)\n",
630 moduleNames[e_mdl], cprid, 'a' + hslb, _exp, _run, _eve, _eve, nword, buf2cnttrg20, buf3dddd, datacnttrg32);
631 m_cntr_nwn_badtrg[e_mdl]++;
632 something_bad = true;
633 }
634 }
635 if (! something_bad) m_cntr_nwn[e_mdl]++;
636
637 } // for hslb
638 }
639 }
640
641 if (m_cc_check) {
642 for (int i = 0; i < raw_trgarray.getEntries(); i++) {
643 for (int j = 0; j < raw_trgarray[i]->GetNumEntries(); j++) {
644 unsigned cprid = raw_trgarray[i]->GetNodeID(j);
645// int _exp = raw_trgarray[i]->GetExpNo(j);
646// int _run = raw_trgarray[i]->GetRunNo(j);
647 unsigned _eveRaw = raw_trgarray[i]->GetEveNo(j);
648 if ((0x11000001 <= cprid && cprid <= 0x1100000a) ||
649 (0x15000001 <= cprid && cprid <= 0x15000002) ||
650 (0x12000001 == cprid) // TOPTRG
651 ) {
652 } else {
653 continue;
654 }
655 for (int hslb = 0; hslb < 2; hslb++) {
656 int nclks = m_nclk_cdctrg;
657 unsigned _hdr_nwd = 3;
658 unsigned e_mdl = -1;
659 unsigned _nwd = -1;
660 if (m_cpr_top == cprid && hslb == 0 && m_on_top) {
661 e_mdl = e_top; _nwd = m_nwd_top;
662 nclks = m_nclk_top; _hdr_nwd = m_hdr_nwd_top;
663 } else if (m_cpr_2d0 == cprid && hslb == 0 && m_on_2d0) {
664 e_mdl = e_2d0; _nwd = m_nwd_2ds; _hdr_nwd = m_hdr_nwd_2d0;
665 } else if (m_cpr_2d1 == cprid && hslb == 1 && m_on_2d1) {
666 e_mdl = e_2d1; _nwd = m_nwd_2ds; _hdr_nwd = m_hdr_nwd_2d1;
667 } else if (m_cpr_2d2 == cprid && hslb == 0 && m_on_2d2) {
668 e_mdl = e_2d2; _nwd = m_nwd_2ds; _hdr_nwd = m_hdr_nwd_2d2;
669 } else if (m_cpr_2d3 == cprid && hslb == 1 && m_on_2d3) {
670 e_mdl = e_2d3; _nwd = m_nwd_2ds; _hdr_nwd = m_hdr_nwd_2d3;
671 } else if (m_cpr_3d0 == cprid && hslb == 0 && m_on_3d0) {
672 e_mdl = e_3d0; _nwd = m_nwd_3ds; _hdr_nwd = m_hdr_nwd_3d0;
673 } else if (m_cpr_3d1 == cprid && hslb == 1 && m_on_3d1) {
674 e_mdl = e_3d1; _nwd = m_nwd_3ds; _hdr_nwd = m_hdr_nwd_3d1;
675 } else if (m_cpr_3d2 == cprid && hslb == 0 && m_on_3d2) {
676 e_mdl = e_3d2; _nwd = m_nwd_3ds; _hdr_nwd = m_hdr_nwd_3d2;
677 } else if (m_cpr_3d3 == cprid && hslb == 1 && m_on_3d3) {
678 e_mdl = e_3d3; _nwd = m_nwd_3ds; _hdr_nwd = m_hdr_nwd_3d3;
679 } else if (m_cpr_nn0 == cprid && hslb == 0 && m_on_nn0) {
680 e_mdl = e_nn0; _nwd = m_nwd_nns; _hdr_nwd = m_hdr_nwd_nn0;
681 } else if (m_cpr_nn1 == cprid && hslb == 1 && m_on_nn1) {
682 e_mdl = e_nn1; _nwd = m_nwd_nns; _hdr_nwd = m_hdr_nwd_nn1;
683 } else if (m_cpr_nn2 == cprid && hslb == 0 && m_on_nn2) {
684 e_mdl = e_nn2; _nwd = m_nwd_nns; _hdr_nwd = m_hdr_nwd_nn2;
685 } else if (m_cpr_nn3 == cprid && hslb == 1 && m_on_nn3) {
686 e_mdl = e_nn3; _nwd = m_nwd_nns; _hdr_nwd = m_hdr_nwd_nn3;
687 } else if (m_cpr_sl0 == cprid && hslb == 0 && m_on_sl0) {
688 e_mdl = e_sl0; _nwd = m_nwd_sl0; _hdr_nwd = m_hdr_nwd_sl0;
689 } else if (m_cpr_sl1 == cprid && hslb == 0 && m_on_sl1) {//ex w/ 4
690 e_mdl = e_sl1; _nwd = m_nwd_sl1; _hdr_nwd = m_hdr_nwd_sl1;
691 } else if (m_cpr_sl2 == cprid && hslb == 0 && m_on_sl2) {
692 e_mdl = e_sl2; _nwd = m_nwd_sl2; _hdr_nwd = m_hdr_nwd_sl2;
693 } else if (m_cpr_sl3 == cprid && hslb == 1 && m_on_sl3) {
694 e_mdl = e_sl3; _nwd = m_nwd_sl3; _hdr_nwd = m_hdr_nwd_sl3;
695 } else if (m_cpr_sl4 == cprid && hslb == 1 && m_on_sl4) {//ex w/ 1
696 e_mdl = e_sl4; _nwd = m_nwd_sl4; _hdr_nwd = m_hdr_nwd_sl4;
697 } else if (m_cpr_sl5 == cprid && hslb == 1 && m_on_sl5) {
698 e_mdl = e_sl5; _nwd = m_nwd_sl5; _hdr_nwd = m_hdr_nwd_sl5;
699 } else if (m_cpr_sl6 == cprid && hslb == 0 && m_on_sl6) {
700 e_mdl = e_sl6; _nwd = m_nwd_sl6; _hdr_nwd = m_hdr_nwd_sl6;
701 } else if (m_cpr_sl8 == cprid && hslb == 0 && m_on_sl8) {
702 e_mdl = e_sl8; _nwd = m_nwd_sl8; _hdr_nwd = m_hdr_nwd_sl8;
703 } else if (m_cpr_etf == cprid && hslb == 1 && m_on_etf) {
704 e_mdl = e_etf; _nwd = m_nwd_etf; _hdr_nwd = m_hdr_nwd_etf;
705 } else if (m_cpr_gdl == cprid && hslb == 0 && m_on_gdl) {
706 e_mdl = e_gdl; _nwd = m_nwd_gdl;
707 nclks = m_nclk_gdl; _hdr_nwd = 6;
708 } else {
709 continue;
710 }
711 unsigned nword = raw_trgarray[i]->GetDetectorNwords(j, hslb);
712 if (nword == _hdr_nwd) continue;
713 if (nword != _nwd) {
714 printf("bu:Nword mismatch:nword(%u),expected(%u). eve(%u), %s(0x%x%c), nword(%u)\n",
715 nword, _nwd, _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword);
716 }
717 int* buf = raw_trgarray[i]->GetDetectorBuffer(j, hslb);
718 int i47 = _hdr_nwd + (nclks - 1) * (_nwd - _hdr_nwd) / nclks;
719 if (i47 > int(nword) - 1) {
720 if (m_print_clkcyc_err)
721 printf("bp:data truncation. eve(%u), %s(0x%x%c), nword(%u)\n",
722 _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword);
723 break;
724 }
725 unsigned buf47 = buf[i47];
726 bool break_this_module = false;
727 bool cc_disorder = false;
728
729 unsigned prev_cc = buf47 & 0xffff;
730 if (e_mdl == e_gdl) {
731 prev_cc = (buf47 >> 4) & 0xfff;
732 }
733 char ccc[400];
734 sprintf(ccc, "i47(%d),buf47(0x%x),%x", i47, buf47, prev_cc);
735 std::string ccs(ccc);
736 for (int clk = 0; clk < nclks; clk++) {
737 unsigned ibuf = _hdr_nwd + clk * (nword - _hdr_nwd) / nclks;
738 if (ibuf > nword - 1) {
739 if (m_print_clkcyc_err)
740 printf("bq:data truncation. eve(%u), %s(0x%x%c), nword(%u)\n",
741 _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword);
742 break_this_module = true;
743 m_cntr_bad_nwd[e_mdl]++;
744 break;
745 }
746 unsigned ddddcc = buf[ibuf];
747 unsigned dddd = (ddddcc >> 16) & 0xffff;
748 unsigned cc = ddddcc & 0xffff;
749 if (e_mdl == e_gdl) {
750 cc = (ddddcc >> 4) & 0xfff;
751 }
752 h_0->SetBinContent(4 * e_mdl + 4, clk + 1, cc);
753 sprintf(ccc, ",%x", cc);
754 if ((int)e_2d0 <= (int)e_mdl && (int)e_mdl <= (int)e_2d3 &&
755 m_nwd_2ds == 3939/* =2624bit */) {
756 // (downto) = wd i
757 // ( ) = 75 300
758 // ( ) = 76 304
759 // ( ) = 79 316
760 // (79, 72) = 79 318 addr_o
761 // (71, 64) = 79 319 addr_i
762 // (63, 56) = 80 320 nworking
763 // (55, 48) = 80 321 rennum
764 // (47, 40) = 80 322 itarget
765 // (39, 8) = 80 323-326 = cnttrg
766 // ( 7, 0) = 81
767 int i_addr = ibuf + 79; // addr_i and addr_o
768 int val_addr = buf[i_addr];
769 int addr_o = (val_addr >> 8) & 0xff;
770 int addr_i = val_addr & 0xff;
771 sprintf(ccc, ",%x(%d,%d)", cc, addr_i, addr_o);
772 }
773 ccs += ccc;
774 if (e_mdl == e_gdl) {
775 dddd = (ddddcc >> 24);
776 if (dddd != 0xdd) {
777 if (m_print_clkcyc_err)
778 printf("br:dddd not found. eve(%u), %s(0x%x%c), nword(%u)\n",
779 _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword);
780 break_this_module = true;
781 m_cntr_bad_ddd[e_mdl]++;
782 break;
783 }
784 } else if (dddd != 0xdddd) {
785 if (m_print_clkcyc_err)
786 printf("bs:dddd not found. eve(%u), %s(0x%x%c), nword(%u)\n",
787 _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword);
788 break_this_module = true;
789 m_cntr_bad_ddd[e_mdl]++;
790 break;
791 }
792 // 0x500 == 1280
793 if (!(cc == prev_cc + 4 || prev_cc == cc + 4 * (nclks - 1) || prev_cc + 0x500 == cc + 4 * (nclks - 1)
794 || cc + 0x500 == prev_cc + 4)) {
795 cc_disorder = true;
796 ccs += "*";
797 }
798 prev_cc = cc;
799 } // clk
800 if (break_this_module) continue; // bad dddd
801 if (cc_disorder) {
802 m_cntr_bad_odr[e_mdl]++;
803 if (m_print_cc) {
804 printf("bt:ccdisorder: eve(%u), %s(0x%x%c), nword(%u), %s\n",
805 _eveRaw, moduleNames[e_mdl], cprid, 'a' + hslb, nword, ccs.c_str());
806 }
807 } else {
808 m_cntr_good_odr[e_mdl]++;
809 }
810 } // hslb
811 } // j < raw_trgarray[i]->GetNumEntries()
812 } // i < raw_trgarray.getEntries()
813 }
814
815 if (! m_histRecord) {
816 h_0->Delete();
817 }
818}
819
HistoModule()
Constructor.
Definition HistoModule.h:32
void setDescription(const std::string &description)
Sets the description of the module.
Definition Module.cc:214
void setPropertyFlags(unsigned int propertyFlags)
Sets the flags for the module properties.
Definition Module.cc:208
@ c_ParallelProcessingCertified
This module can be run in parallel processing mode safely (All I/O must be done through the data stor...
Definition Module.h:80
Accessor to arrays stored in the data store.
Definition StoreArray.h:113
int getEntries() const
Get the number of objects in the array.
Definition StoreArray.h:216
Type-safe access to single objects in the data store.
Definition StoreObjPtr.h:96
static const unsigned int c_nModules
number of modules
virtual void initialize() override
initialize
virtual void event() override
Event.
virtual void defineHisto() override
Define Histogram.
void addParam(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
Definition Module.h:559
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Definition Module.h:649
Abstract base class for different kinds of events.
STL namespace.