1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | #include <stdint.h> |
9 | #include <vector> |
10 | |
11 | #include "JEventProcessor_DAQ_online.h" |
12 | #include <JANA/JApplication.h> |
13 | #include <JANA/JFactory.h> |
14 | |
15 | using namespace std; |
16 | using namespace jana; |
17 | |
18 | #include <DAQ/DF1TDCHit.h> |
19 | #include <DAQ/Df250PulseIntegral.h> |
20 | #include <DAQ/JEventSource_EVIO.h> |
21 | #include <TTAB/DTranslationTable.h> |
22 | |
23 | #include <TDirectory.h> |
24 | #include <TH2.h> |
25 | #include <TH1.h> |
26 | #include <TProfile.h> |
27 | #include <TProfile2D.h> |
28 | #include <TROOT.h> |
29 | |
30 | |
31 | static const int highcratenum=100; |
32 | |
33 | static TH2I *daq_occ_crates[highcratenum]; |
34 | static TProfile2D *daq_ped_crates[highcratenum]; |
35 | static TProfile2D *daq_TDClocked_crates[highcratenum]; |
36 | static TProfile2D *daq_TDCovr_crates[highcratenum]; |
37 | static TProfile *daq_hits_per_event; |
38 | static TProfile *daq_words_per_event; |
39 | static TH1D *daq_event_size; |
40 | static TH1D *daq_event_tdiff; |
41 | static TH1D *daq_words_by_type; |
42 | static bool ttab_labels_set = false; |
43 | |
44 | |
45 | extern "C"{ |
46 | void InitPlugin(JApplication *app){ |
47 | InitJANAPlugin(app); |
48 | app->AddProcessor(new JEventProcessor_DAQ_online()); |
49 | } |
50 | } |
51 | |
52 | |
53 | |
54 | |
55 | |
56 | JEventProcessor_DAQ_online::JEventProcessor_DAQ_online() |
57 | { |
58 | |
59 | } |
60 | |
61 | |
62 | |
63 | |
64 | JEventProcessor_DAQ_online::~JEventProcessor_DAQ_online() |
65 | { |
66 | |
67 | } |
68 | |
69 | |
70 | |
71 | |
72 | jerror_t JEventProcessor_DAQ_online::init(void) |
73 | { |
74 | printf("JEventProcessor_DAQ_online::init()\n"); |
75 | |
76 | |
77 | japp->RootWriteLock(); |
78 | |
79 | |
80 | maindir = gDirectory(TDirectory::CurrentDirectory()); |
81 | daqdir = maindir->mkdir("DAQ"); |
82 | daqdir->cd(); |
83 | |
84 | |
85 | for (int i=0; i<highcratenum; i++) { |
86 | daq_occ_crates[i] = NULL__null; |
87 | daq_ped_crates[i] = NULL__null; |
88 | daq_TDClocked_crates[i] = NULL__null; |
89 | daq_TDCovr_crates[i] = NULL__null; |
90 | } |
91 | |
92 | daq_hits_per_event = new TProfile("daq_hits_per_event", "Hits/event vs. rocid", 100, 0.5, 100.5); |
93 | daq_words_per_event = new TProfile("daq_words_per_event", "words/event vs. rocid", 100, 0.5, 100.5); |
94 | daq_event_size = new TH1D("daq_event_size", "Event size in kB", 1000, 0.0, 1.0E3); |
95 | daq_event_tdiff = new TH1D("daq_event_tdiff", "Time between events", 10000, 0.0, 1.0E2); |
96 | daq_words_by_type = new TH1D("daq_words_by_type", "Number of words in EVIO file by type", kNEVIOWordTypes, 0, (double)kNEVIOWordTypes); |
97 | |
98 | daq_words_per_event->GetXaxis()->SetBinLabel(1 ,"Trigger Bank"); |
99 | daq_words_per_event->GetXaxis()->SetBinLabel(99 ,"Residual"); |
100 | |
101 | daq_event_size->SetXTitle("Total event size (kB)"); |
102 | daq_event_tdiff->SetXTitle("#deltat between events (ms)"); |
103 | |
104 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kUnknown, "unknown"); |
105 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kEVIOEventNumber, "Event Number Word"); |
106 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kEVIOTimestamp, "Timestamp"); |
107 | |
108 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250BlockHeader, "f250 Block Header"); |
109 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250BlockTrailer, "f250 Block Trailer"); |
110 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250EventHeader, "f250 Event Header"); |
111 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250TriggerTime, "f250 Trigger Time"); |
112 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250WindowRawData, "f250 Window Raw Data"); |
113 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250WindowSum, "f250 Window Sum"); |
114 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250PulseRawData, "f250 Pulse Raw Data"); |
115 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250PulseIntegral, "f250 Pulse Integral"); |
116 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250PulseTime, "f250 Pulse Time"); |
117 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250PulsePedestal, "f250 Pulse Pedestal"); |
118 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250EventTrailer, "f250 Event Trailer"); |
119 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250DataNotValid, "f250 Data Not Valid"); |
120 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf250Filler, "f250 Filler Word"); |
121 | |
122 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125BlockHeader, "f125 Block Header"); |
123 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125BlockTrailer, "f125 Block Trailer"); |
124 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125EventHeader, "f125 Event Header"); |
125 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125TriggerTime, "f125 Trigger Time"); |
126 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125WindowRawData, "f125 Window Raw Data"); |
127 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125CDCPulse, "f125 CDC Pulse"); |
128 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125FDCPulse6, "f125 FDC Pulse (integral)"); |
129 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125FDCPulse9, "f125 FDC Pulse (peak)"); |
130 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125PulseIntegral, "f125 Pulse Integral"); |
131 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125PulseTime, "f125 Pulse Time"); |
132 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125PulsePedestal, "f125 Pulse Pedestal"); |
133 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125EventTrailer, "f125 Event Trailer"); |
134 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125DataNotValid, "f125 Data Not Valid"); |
135 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kf125Filler, "f125 Filler Word"); |
136 | |
137 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2BlockHeader, "F1v2 Block Header"); |
138 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2BLockTrailer, "F1v2 Block Trailer"); |
139 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2EventHeader, "F1v2 Event Header"); |
140 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2TriggerTime, "F1v2 Trigger Time"); |
141 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2ChipHeader, "F1v2 Chip Header"); |
142 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2Data, "F1v2 Data"); |
143 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2Filler, "F1v2 Filler"); |
144 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v2BreakWord, "F1v2 Break Word"); |
145 | |
146 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3BlockHeader, "F1v3 Block Header"); |
147 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3BLockTrailer, "F1v3 Block Trailer"); |
148 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3EventHeader, "F1v3 Event Header"); |
149 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3TriggerTime, "F1v3 Trigger Time"); |
150 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3ChipHeader, "F1v3 Chip Header"); |
151 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3Data, "F1v3 Data"); |
152 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3Filler, "F1v3 Filler"); |
153 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF1v3BreakWord, "F1v3 Break Word"); |
154 | |
155 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190GlobalHeader, "CAEN1190 GLobal Header"); |
156 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190GlobalTrailer, "CAEN1190 Global Trailer"); |
157 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190GlobalTriggerTime, "CAEN1190 Trigger Time"); |
158 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190TDCHeader, "CAEN1190 TDC Header"); |
159 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190TDCData, "CAEN1190 TDC Data"); |
160 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190TDCError, "CAEN1190 TDC Error"); |
161 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190TDCTrailer, "CAEN1190 TDC Trailer"); |
162 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kCAEN1190Filler, "CAEN1190 Filler"); |
163 | |
164 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kConfig, "DAQ Config"); |
165 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kConfigf250, "DAQ Config f250"); |
166 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kConfigf125, "DAQ Config f125"); |
167 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kConfigF1, "DAQ Config F1"); |
168 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kConfigCAEN1190, "DAQ Config CAEN1190"); |
169 | |
170 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kEPICSheader, "EPICS header"); |
171 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kEPICSdata, "EPICS data"); |
172 | |
173 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kF800FAFA, "0xf800fafa"); |
174 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kD00DD00D, "0xd00dd00d"); |
175 | |
176 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kTotWords, "Total words in all events"); |
177 | daq_words_by_type->GetXaxis()->SetBinLabel(1 + kNevents, "Number of events"); |
178 | |
179 | |
180 | maindir->cd(); |
181 | |
182 | |
183 | japp->RootUnLock(); |
184 | |
185 | return NOERROR; |
186 | } |
187 | |
188 | |
189 | |
190 | |
191 | void JEventProcessor_DAQ_online::AddROCIDLabels(JEventLoop *loop) |
192 | { |
193 | |
194 | |
195 | |
196 | |
197 | const DTranslationTable *ttab = NULL__null; |
198 | loop->GetSingle(ttab); |
199 | |
200 | japp->RootWriteLock(); |
201 | |
202 | |
203 | for(uint32_t rocid=2; rocid<99; rocid++){ |
204 | |
205 | |
206 | bool found_chan = false; |
207 | daq_hits_per_event->GetXaxis()->SetBinLabel(rocid, ""); |
208 | daq_words_per_event->GetXaxis()->SetBinLabel(rocid, ""); |
209 | for(uint32_t slot=2; slot<24; slot++){ |
210 | for(uint32_t channel=0; channel<3; channel++){ |
211 | try{ |
212 | DTranslationTable::csc_t csc = {rocid, slot, channel}; |
213 | const DTranslationTable::DChannelInfo &chinfo = ttab->GetDetectorIndex(csc); |
214 | daq_hits_per_event->GetXaxis()->SetBinLabel(rocid, ttab->DetectorName(chinfo.det_sys).c_str()); |
215 | daq_words_per_event->GetXaxis()->SetBinLabel(rocid, ttab->DetectorName(chinfo.det_sys).c_str()); |
216 | found_chan = true; |
217 | break; |
218 | }catch(JException &e){ |
219 | |
220 | } |
221 | } |
222 | if(found_chan) break; |
223 | } |
224 | } |
225 | |
226 | japp->RootUnLock(); |
227 | } |
228 | |
229 | |
230 | |
231 | |
232 | jerror_t JEventProcessor_DAQ_online::brun(JEventLoop *eventLoop, int32_t runnumber) |
233 | { |
234 | |
235 | return NOERROR; |
236 | } |
237 | |
238 | |
239 | |
240 | |
241 | jerror_t JEventProcessor_DAQ_online::evnt(JEventLoop *loop, uint64_t eventnumber) |
242 | { |
243 | |
244 | |
245 | |
246 | |
247 | |
248 | |
249 | vector<const DF1TDCHit*> f1tdchits; |
250 | vector<const Df250PulseIntegral*> f250PIs; |
251 | vector<const Df125PulseIntegral*> f125PIs; |
252 | vector<const Df125CDCPulse*> f125CDCs; |
253 | vector<const Df125FDCPulse*> f125FDCs; |
254 | vector<const DCAEN1290TDCHit*> caen1290hits; |
255 | |
256 | loop->Get(f1tdchits); |
257 | loop->Get(f250PIs); |
258 | loop->Get(f125PIs); |
259 | loop->Get(f125CDCs); |
260 | loop->Get(f125FDCs); |
261 | loop->Get(caen1290hits); |
262 | |
263 | ParseEventSize(loop->GetJEvent()); |
264 | |
265 | |
266 | if(!ttab_labels_set){ |
267 | ttab_labels_set = true; |
268 | AddROCIDLabels(loop); |
269 | } |
270 | |
271 | |
272 | uint32_t Nhits_rocid[101]; |
273 | for(uint32_t rocid=0; rocid<101; rocid++) Nhits_rocid[rocid] = 0; |
274 | |
275 | |
276 | japp->RootWriteLock(); |
277 | |
278 | if (daqdir!=NULL__null) daqdir->cd(); |
279 | |
280 | |
281 | |
282 | for(unsigned int i=0; i<f1tdchits.size(); i++) { |
283 | const DF1TDCHit *hit = f1tdchits[i]; |
284 | int rocid = hit->rocid; |
285 | int slot = hit->slot; |
286 | int channel = hit->channel; |
287 | int data_word = hit->data_word; |
288 | |
289 | if(rocid>=0 && rocid<=100) Nhits_rocid[rocid]++; |
290 | |
291 | if (daq_occ_crates[rocid]==NULL__null) { |
292 | printf("JEventProcessor_DAQ_online::evnt creating occupancy histogram for crate %i\n",rocid); |
293 | char cratename[255],title[255]; |
294 | sprintf(cratename,"daq_occ_crate%i",rocid); |
295 | sprintf(title,"Crate %i occupancy (TDC);Slot;Channel",rocid); |
296 | daq_occ_crates[rocid] = new TH2I(cratename,title,21,0.5,21.5,32,-0.5,31.5); |
297 | daq_occ_crates[rocid]->SetStats(0); |
298 | sprintf(cratename,"daq_TDClocked_crate%i",rocid); |
299 | sprintf(title,"Crate %i TDC lock status (TDC);Slot;Channel",rocid); |
300 | daq_TDClocked_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,32,-0.5,31.5); |
301 | daq_TDClocked_crates[rocid]->SetStats(0); |
302 | sprintf(cratename,"daq_TDCovr_crate%i",rocid); |
303 | sprintf(title,"Crate %i TDC overflow status (TDC);Slot;Channel",rocid); |
304 | daq_TDCovr_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,32,-0.5,31.5); |
305 | daq_TDCovr_crates[rocid]->SetStats(0); |
306 | |
307 | } |
308 | daq_occ_crates[rocid]->Fill(slot,channel); |
309 | daq_TDClocked_crates[rocid]->Fill(slot,channel,(data_word>>26)&(1)); |
310 | daq_TDCovr_crates[rocid]->Fill(slot,channel,(data_word>>25)&(1)); |
311 | daq_TDCovr_crates[rocid]->Fill(slot,channel,(data_word>>24)&(1)); |
312 | } |
313 | |
314 | |
315 | for(unsigned int i=0; i<f250PIs.size(); i++) { |
316 | const Df250PulseIntegral *hit = f250PIs[i]; |
317 | int rocid = hit->rocid; |
318 | int slot = hit->slot; |
319 | int channel = hit->channel; |
320 | |
321 | if(rocid>=0 && rocid<=100) { |
322 | Nhits_rocid[rocid]++; |
323 | |
324 | if (daq_occ_crates[rocid]==NULL__null) { |
325 | printf("JEventProcessor_DAQ_online::evnt creating occupancy histogram for crate %i\n",rocid); |
326 | char cratename[255],title[255]; |
327 | sprintf(cratename,"daq_occ_crate%i",rocid); |
328 | sprintf(title,"Crate %i occupancy (F250);Slot;Channel",rocid); |
329 | daq_occ_crates[rocid] = new TH2I(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
330 | daq_occ_crates[rocid]->SetStats(0); |
331 | } |
332 | daq_occ_crates[rocid]->Fill(slot,channel); |
333 | |
334 | if (daq_ped_crates[rocid]==NULL__null) { |
335 | printf("JEventProcessor_DAQ_online::evnt creating pedestal histogram for crate %i\n",rocid); |
336 | char cratename[255],title[255]; |
337 | sprintf(cratename,"daq_ped_crate%i",rocid); |
338 | sprintf(title,"Crate %i Average Pedestal (F250);Slot;Channel",rocid); |
339 | daq_ped_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
340 | daq_ped_crates[rocid]->SetStats(0); |
341 | } |
342 | if (hit->pedestal > 0) { |
343 | daq_ped_crates[rocid]->Fill(slot,channel,hit->pedestal); |
344 | } |
345 | } |
346 | } |
347 | |
348 | |
349 | for(unsigned int i=0; i<f125PIs.size(); i++) { |
350 | const Df125PulseIntegral *hit = f125PIs[i]; |
351 | int rocid = hit->rocid; |
352 | int slot = hit->slot; |
353 | int channel = hit->channel; |
354 | |
355 | if(hit->emulated) continue; |
356 | |
357 | if(rocid>=0 && rocid<=100) { |
358 | Nhits_rocid[rocid]++; |
359 | |
360 | if (daq_occ_crates[rocid]==NULL__null) { |
361 | printf("JEventProcessor_DAQ_online::evnt creating occupancy histogram for crate %i\n",rocid); |
362 | char cratename[255],title[255]; |
363 | sprintf(cratename,"daq_occ_crate%i",rocid); |
364 | sprintf(title,"Crate %i occupancy (F125);Slot;Channel",rocid); |
365 | daq_occ_crates[rocid] = new TH2I(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
366 | daq_occ_crates[rocid]->SetStats(0); |
367 | } |
368 | daq_occ_crates[rocid]->Fill(slot,channel); |
369 | |
370 | if (daq_ped_crates[rocid]==NULL__null) { |
371 | printf("JEventProcessor_DAQ_online::evnt creating pedestal histogram for crate %i\n",rocid); |
372 | char cratename[255],title[255]; |
373 | sprintf(cratename,"daq_ped_crate%i",rocid); |
374 | sprintf(title,"Crate %i Average Pedestal (F125);Slot;Channel",rocid); |
375 | daq_ped_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
376 | daq_ped_crates[rocid]->SetStats(0); |
377 | } |
378 | if (hit->pedestal > 0) { |
379 | daq_ped_crates[rocid]->Fill(slot,channel,hit->pedestal); |
380 | } |
381 | } |
382 | } |
383 | |
384 | |
385 | for(unsigned int i=0; i<f125CDCs.size(); i++) { |
386 | const Df125CDCPulse *hit = f125CDCs[i]; |
387 | int rocid = hit->rocid; |
388 | int slot = hit->slot; |
389 | int channel = hit->channel; |
390 | |
391 | if(rocid>=0 && rocid<=100) { |
392 | Nhits_rocid[rocid]++; |
393 | |
394 | if (daq_occ_crates[rocid]==NULL__null) { |
395 | printf("JEventProcessor_DAQ_online::evnt creating occupancy histogram for crate %i\n",rocid); |
396 | char cratename[255],title[255]; |
397 | sprintf(cratename,"daq_occ_crate%i",rocid); |
398 | sprintf(title,"Crate %i occupancy (F125);Slot;Channel",rocid); |
399 | daq_occ_crates[rocid] = new TH2I(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
400 | daq_occ_crates[rocid]->SetStats(0); |
401 | } |
402 | daq_occ_crates[rocid]->Fill(slot,channel); |
403 | |
404 | if (daq_ped_crates[rocid]==NULL__null) { |
405 | printf("JEventProcessor_DAQ_online::evnt creating pedestal histogram for crate %i\n",rocid); |
406 | char cratename[255],title[255]; |
407 | sprintf(cratename,"daq_ped_crate%i",rocid); |
408 | sprintf(title,"Crate %i Average Pedestal (F125);Slot;Channel",rocid); |
409 | daq_ped_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
410 | daq_ped_crates[rocid]->SetStats(0); |
411 | } |
412 | if (hit->pedestal > 0) { |
413 | daq_ped_crates[rocid]->Fill(slot,channel,hit->pedestal); |
414 | } |
415 | } |
416 | } |
417 | |
418 | |
419 | for(unsigned int i=0; i<f125FDCs.size(); i++) { |
420 | const Df125FDCPulse *hit = f125FDCs[i]; |
421 | int rocid = hit->rocid; |
422 | int slot = hit->slot; |
423 | int channel = hit->channel; |
424 | |
425 | if(rocid>=0 && rocid<=100) { |
426 | Nhits_rocid[rocid]++; |
427 | |
428 | if (daq_occ_crates[rocid]==NULL__null) { |
429 | printf("JEventProcessor_DAQ_online::evnt creating occupancy histogram for crate %i\n",rocid); |
430 | char cratename[255],title[255]; |
431 | sprintf(cratename,"daq_occ_crate%i",rocid); |
432 | sprintf(title,"Crate %i occupancy (F125);Slot;Channel",rocid); |
433 | daq_occ_crates[rocid] = new TH2I(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
434 | daq_occ_crates[rocid]->SetStats(0); |
435 | } |
436 | daq_occ_crates[rocid]->Fill(slot,channel); |
437 | |
438 | if (daq_ped_crates[rocid]==NULL__null) { |
439 | printf("JEventProcessor_DAQ_online::evnt creating pedestal histogram for crate %i\n",rocid); |
440 | char cratename[255],title[255]; |
441 | sprintf(cratename,"daq_ped_crate%i",rocid); |
442 | sprintf(title,"Crate %i Average Pedestal (F125);Slot;Channel",rocid); |
443 | daq_ped_crates[rocid] = new TProfile2D(cratename,title,21,0.5,21.5,16,-0.5,15.5); |
444 | daq_ped_crates[rocid]->SetStats(0); |
445 | } |
446 | if (hit->pedestal > 0) { |
447 | daq_ped_crates[rocid]->Fill(slot,channel,hit->pedestal); |
448 | } |
449 | } |
450 | } |
451 | |
452 | |
453 | for(unsigned int i=0; i<caen1290hits.size(); i++) { |
454 | const DCAEN1290TDCHit *hit = caen1290hits[i]; |
455 | int rocid = hit->rocid; |
456 | |
457 | |
458 | |
459 | if(rocid>=0 && rocid<=100) Nhits_rocid[rocid]++; |
460 | } |
461 | |
462 | |
463 | for(uint32_t rocid=0; rocid<101; rocid++) daq_hits_per_event->Fill(rocid, Nhits_rocid[rocid]); |
464 | |
465 | maindir->cd(); |
466 | |
467 | japp->RootUnLock(); |
468 | |
469 | return NOERROR; |
470 | } |
471 | |
472 | |
473 | |
474 | |
475 | void JEventProcessor_DAQ_online::ParseEventSize(JEvent &event) |
476 | { |
477 | |
478 | |
479 | |
480 | if(event.GetJEventSource()->className() != string("JEventSource_EVIO")){ |
481 | static bool warned = false; |
482 | if(!warned){ |
483 | cout << "WARNING: This is not an event source of type JEventSource_EVIO!" << endl; |
484 | cout << " Event size statistics filling unavailable!" << endl; |
485 | warned = true; |
486 | } |
487 | return; |
488 | } |
489 | |
490 | void *ref = event.GetRef(); |
491 | if(!ref) return; |
492 | uint32_t *istart = JEventSource_EVIO::GetEVIOBufferFromRef(ref); |
493 | uint32_t evio_buffsize = JEventSource_EVIO::GetEVIOBufferSizeFromRef(ref); |
494 | uint32_t evio_buffwords = evio_buffsize/sizeof(uint32_t); |
495 | uint32_t *iend = &istart[evio_buffwords]; |
496 | |
497 | if( istart==NULL__null ) return; |
498 | if( (evio_buffwords>=10) && (istart[7]==0xc0da0100) ){ |
499 | |
500 | istart= &istart[8]; |
501 | evio_buffsize -= 8*sizeof(uint32_t); |
| Value stored to 'evio_buffsize' is never read |
502 | evio_buffwords -= 8; |
503 | } |
504 | |
505 | |
506 | if( evio_buffwords >= 4 ){ |
507 | if( istart[1] == (0x60<<16) + (0xD<<8) + (0x1<<0) ){ |
508 | if( istart[2] == (0x61<<24) + (0x1<<16) + (0x1<<0) ){ |
509 | |
510 | japp->RootWriteLock(); |
511 | daq_words_by_type->Fill(kEPICSheader, 3.0); |
512 | daq_words_by_type->Fill(kEPICSdata, istart[0]/sizeof(uint32_t) - 3); |
513 | japp->RootUnLock(); |
514 | return; |
515 | } |
516 | } |
517 | } |
518 | |
519 | |
520 | uint32_t physics_event_len = istart[0]; |
521 | if( (istart[1] & 0xFF001000) != 0xFF001000 ) return; |
522 | |
523 | |
524 | uint32_t trigger_bank_len = istart[2]; |
525 | if( (istart[3] & 0xFF202000) != 0xFF202000 ) return; |
526 | uint64_t tlo = istart[2+5]; |
527 | uint64_t thi = istart[2+6]; |
528 | uint64_t timestamp = (thi<<32) + (tlo<<0); |
529 | |
530 | |
531 | uint32_t Nwords[100]; |
532 | uint32_t word_stats[kNEVIOWordTypes]; |
533 | for(uint32_t rocid=0; rocid<100; rocid++) Nwords[rocid] = 0; |
534 | for(uint32_t i=0; i<kNEVIOWordTypes; i++) word_stats[i] = 0; |
535 | |
536 | word_stats[kNevents]++; |
537 | word_stats[kTotWords] += evio_buffwords; |
538 | |
539 | |
540 | uint32_t *iptr = &istart[3+trigger_bank_len]; |
541 | while(iptr < iend){ |
542 | |
543 | uint32_t len = *iptr; |
544 | uint32_t rocid = (iptr[1]>>16) & 0XFF; |
545 | |
546 | if(rocid<100) Nwords[rocid] += len+1; |
547 | |
548 | uint32_t *imyend = &iptr[len+1]; |
549 | if(imyend > iend) imyend = iend; |
550 | |
551 | DataWordStats(iptr, imyend, word_stats); |
552 | |
553 | iptr = &iptr[len +1]; |
554 | } |
555 | |
556 | |
557 | japp->RootWriteLock(); |
558 | |
559 | |
560 | |
561 | |
562 | |
563 | |
564 | |
565 | uint32_t event_num = event.GetEventNumber(); |
566 | static uint32_t recent_event_nums[128]; |
567 | static uint64_t recent_timestamps[128]; |
568 | static uint32_t ievent = 0; |
569 | for(uint32_t i=0; i<ievent; i++){ |
570 | if(i>=128) break; |
571 | if(recent_event_nums[i] == (event_num-1)){ |
572 | double tdiff = (double)(timestamp - recent_timestamps[i])/250.0E6; |
573 | daq_event_tdiff->Fill(tdiff*1000.0); |
574 | break; |
575 | } |
576 | } |
577 | |
578 | |
579 | uint32_t idx = ievent%128; |
580 | recent_event_nums[idx] = event_num; |
581 | recent_timestamps[idx] = timestamp; |
582 | ievent++; |
583 | |
584 | |
585 | double physics_event_len_kB = (double)((physics_event_len+1)*sizeof(uint32_t))/1024.0; |
586 | daq_event_size->Fill(physics_event_len_kB); |
587 | uint32_t TotalWords = 0; |
588 | for(uint32_t rocid=0; rocid<100; rocid++){ |
589 | daq_words_per_event->Fill(rocid, Nwords[rocid]); |
590 | TotalWords += Nwords[rocid]; |
591 | } |
592 | |
593 | daq_words_per_event->Fill(1, trigger_bank_len+1); |
594 | daq_words_per_event->Fill(99, physics_event_len - trigger_bank_len - TotalWords); |
595 | |
596 | for(uint32_t i=0; i<kNEVIOWordTypes; i++){ |
597 | daq_words_by_type->Fill(i, (double)word_stats[i]); |
598 | } |
599 | |
600 | japp->RootUnLock(); |
601 | |
602 | } |
603 | |
604 | |
605 | |
606 | |
607 | void JEventProcessor_DAQ_online::DataWordStats(uint32_t *iptr, uint32_t *iend, uint32_t *word_stats) |
608 | { |
609 | |
610 | |
611 | |
612 | |
613 | iptr++; |
614 | uint32_t rocid = (*iptr++)>>16 & 0x0FFF; |
615 | while(iptr < iend){ |
616 | uint32_t data_block_bank_len = *iptr++; |
617 | uint32_t *iendbank = &iptr[data_block_bank_len]; |
618 | uint32_t det_id = ((*iptr) >> 16) & 0x0FFF; |
619 | iptr++; |
620 | |
621 | uint32_t Ntoprocess = data_block_bank_len - 1; |
622 | |
623 | #if 0 // I don't know if these words are actually implmented ?? |
624 | word_stats[kEVIOEventNumber]++; |
625 | word_stats[kEVIOTimestamp] += 2; |
626 | iptr++; |
627 | iptr++; |
628 | iptr++; |
629 | Ntoprocess -= 3; |
630 | #endif |
631 | uint32_t *irawdata = iptr; |
632 | |
633 | switch(det_id){ |
634 | case 0: |
635 | case 1: |
636 | case 3: |
637 | case 6: |
638 | case 16: |
639 | case 26: |
640 | ParseJLabModuleData(rocid, iptr, iendbank, word_stats); |
641 | break; |
642 | |
643 | case 20: |
644 | ParseCAEN1190(rocid, iptr, iendbank, word_stats); |
645 | break; |
646 | |
647 | case 0x55: |
648 | ParseModuleConfiguration(rocid, iptr, iendbank, word_stats); |
649 | break; |
650 | default: |
651 | break; |
652 | } |
653 | |
654 | uint32_t Nprocessed = (uint32_t)((uint64_t)iptr - (uint64_t)irawdata)/sizeof(uint32_t); |
655 | if(Nprocessed < Ntoprocess) word_stats[kUnknown] += Ntoprocess - Nprocessed; |
656 | iptr = iendbank; |
657 | } |
658 | |
659 | |
660 | } |
661 | |
662 | |
663 | |
664 | |
665 | void JEventProcessor_DAQ_online::ParseJLabModuleData(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
666 | { |
667 | while(iptr < iend){ |
668 | if(*iptr != 0xf800fafa) break; |
669 | word_stats[kF800FAFA]++; |
670 | iptr++; |
671 | } |
672 | |
673 | uint32_t mod_id = ((*iptr) >> 18) & 0x000F; |
674 | switch(mod_id){ |
675 | case DModuleType::FADC250: Parsef250Bank(rocid, iptr, iend, word_stats); break; |
676 | case DModuleType::FADC125: Parsef125Bank(rocid, iptr, iend, word_stats); break; |
677 | case DModuleType::F1TDC32: ParseF1v2TDCBank(rocid, iptr, iend, word_stats); break; |
678 | case DModuleType::F1TDC48: ParseF1v3TDCBank(rocid, iptr, iend, word_stats); break; |
679 | |
680 | |
681 | } |
682 | } |
683 | |
684 | |
685 | |
686 | |
687 | void JEventProcessor_DAQ_online::Parsef250Bank(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
688 | { |
689 | while(iptr<iend){ |
690 | |
691 | if(((*iptr>>31) & 0x1) == 0) { word_stats[kUnknown]++ ; iptr++; continue;} |
692 | |
693 | uint32_t window_width; |
694 | uint32_t window_words; |
695 | uint32_t data_type = (*iptr>>27) & 0x0F; |
696 | switch(data_type){ |
697 | case 0: word_stats[kf250BlockHeader]++; iptr++; break; |
698 | case 1: word_stats[kf250BlockTrailer]++; iptr++; break; |
699 | case 2: word_stats[kf250EventHeader]++; iptr++; break; |
700 | case 3: |
701 | word_stats[kf250TriggerTime]++; |
702 | iptr++; |
703 | if(((*iptr>>31) & 0x1) == 0){ word_stats[kf250TriggerTime]++; iptr++; } |
704 | break; |
705 | case 4: |
706 | window_width = (*iptr>>0) & 0x0FFF; |
707 | window_words = 1 + ((window_width+1)/2); |
708 | word_stats[kf250WindowRawData] += window_words; |
709 | iptr = &iptr[window_words]; |
710 | break; |
711 | case 7: word_stats[kf250PulseIntegral]++; iptr++; break; |
712 | case 8: word_stats[kf250PulseTime]++; iptr++; break; |
713 | case 10: word_stats[kf250PulsePedestal]++; iptr++; break; |
714 | case 13: word_stats[kf250EventTrailer]++; iptr++; break; |
715 | case 14: word_stats[kf250DataNotValid]++; iptr++; break; |
716 | case 15: word_stats[kf250Filler]++; iptr++; break; |
717 | |
718 | default: word_stats[kUnknown]++; iptr++; break; |
719 | } |
720 | } |
721 | } |
722 | |
723 | |
724 | |
725 | |
726 | void JEventProcessor_DAQ_online::Parsef125Bank(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
727 | { |
728 | while(iptr<iend){ |
729 | |
730 | if(((*iptr>>31) & 0x1) == 0) { word_stats[kUnknown]++ ; iptr++; continue;} |
731 | |
732 | uint32_t window_width; |
733 | uint32_t window_words; |
734 | uint32_t data_type = (*iptr>>27) & 0x0F; |
735 | switch(data_type){ |
736 | case 0: word_stats[kf125BlockHeader]++; iptr++; break; |
737 | case 1: word_stats[kf125BlockTrailer]++; iptr++; break; |
738 | case 2: word_stats[kf125EventHeader]++; iptr++; break; |
739 | case 3: |
740 | word_stats[kf125TriggerTime]++; |
741 | iptr++; |
742 | if(((*iptr>>31) & 0x1) == 0){ word_stats[kf125TriggerTime]++; iptr++; } |
743 | break; |
744 | case 4: |
745 | window_width = (*iptr>>0) & 0x0FFF; |
746 | window_words = 1 + ((window_width+1)/2); |
747 | word_stats[kf125WindowRawData] += window_words; |
748 | iptr = &iptr[window_words]; |
749 | break; |
750 | case 5: word_stats[kf125CDCPulse]++; |
751 | iptr++; |
752 | if(((*iptr>>31) & 0x1) == 0){ word_stats[kf125CDCPulse]++; iptr++; } |
753 | break; |
754 | case 6: word_stats[kf125FDCPulse6]++; |
755 | iptr++; |
756 | if(((*iptr>>31) & 0x1) == 0){ word_stats[kf125FDCPulse6]++; iptr++; } |
757 | break; |
758 | case 7: word_stats[kf125PulseIntegral]++; iptr++; break; |
759 | case 8: word_stats[kf125PulseTime]++; iptr++; break; |
760 | case 9: word_stats[kf125FDCPulse9]++; |
761 | iptr++; |
762 | if(((*iptr>>31) & 0x1) == 0){ word_stats[kf125FDCPulse9]++; iptr++; } |
763 | break; |
764 | case 10: word_stats[kf125PulsePedestal]++; iptr++; break; |
765 | case 13: word_stats[kf125EventTrailer]++; iptr++; break; |
766 | case 14: word_stats[kf125DataNotValid]++; iptr++; break; |
767 | case 15: word_stats[kf125Filler]++; iptr++; break; |
768 | |
769 | default: word_stats[kUnknown]++; iptr++; break; |
770 | } |
771 | } |
772 | } |
773 | |
774 | |
775 | |
776 | |
777 | void JEventProcessor_DAQ_online::ParseF1v2TDCBank(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
778 | { |
779 | while(iptr<iend){ |
780 | switch( (*iptr++) & 0xF8000000 ){ |
781 | case 0xC0000000: word_stats[kF1v2ChipHeader]++; break; |
782 | case 0xB8000000: word_stats[kF1v2Data]++; break; |
783 | case 0xF8000000: word_stats[kF1v2Filler]++; break; |
784 | case 0x80000000: word_stats[kF1v2BlockHeader]++; break; |
785 | case 0x88000000: word_stats[kF1v2BLockTrailer]++; break; |
786 | case 0x90000000: word_stats[kF1v2EventHeader]++; break; |
787 | case 0x98000000: word_stats[kF1v2TriggerTime]++; break; |
788 | case 0xF0000000: word_stats[kF1v2BreakWord]++; break; |
789 | default: word_stats[kUnknown]++; break; |
790 | } |
791 | } |
792 | } |
793 | |
794 | |
795 | |
796 | |
797 | void JEventProcessor_DAQ_online::ParseF1v3TDCBank(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
798 | { |
799 | while(iptr<iend){ |
800 | switch( (*iptr++) & 0xF8000000 ){ |
801 | case 0xC0000000: word_stats[kF1v3ChipHeader]++; break; |
802 | case 0xB8000000: word_stats[kF1v3Data]++; break; |
803 | case 0xF8000000: word_stats[kF1v3Filler]++; break; |
804 | case 0x80000000: word_stats[kF1v3BlockHeader]++; break; |
805 | case 0x88000000: word_stats[kF1v3BLockTrailer]++; break; |
806 | case 0x90000000: word_stats[kF1v3EventHeader]++; break; |
807 | case 0x98000000: word_stats[kF1v3TriggerTime]++; break; |
808 | case 0xF0000000: word_stats[kF1v3BreakWord]++; break; |
809 | default: word_stats[kUnknown]++; break; |
810 | } |
811 | } |
812 | } |
813 | |
814 | |
815 | |
816 | |
817 | void JEventProcessor_DAQ_online::ParseCAEN1190(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
818 | { |
819 | while(iptr<iend){ |
820 | |
821 | |
822 | |
823 | if(*iptr == 0xd00dd00d) { |
824 | word_stats[kD00DD00D]++; |
825 | iptr++; |
826 | continue; |
827 | } |
828 | |
829 | uint32_t type = (*iptr++) >> 27; |
830 | switch(type){ |
831 | case 0b01000: word_stats[kCAEN1190GlobalHeader]++; break; |
832 | case 0b10000: word_stats[kCAEN1190GlobalTrailer]++; break; |
833 | case 0b10001: word_stats[kCAEN1190GlobalTriggerTime]++; break; |
834 | case 0b00001: word_stats[kCAEN1190TDCHeader]++; break; |
835 | case 0b00000: word_stats[kCAEN1190TDCData]++; break; |
836 | case 0b00100: word_stats[kCAEN1190TDCError]++; break; |
837 | case 0b00011: word_stats[kCAEN1190TDCTrailer]++; break; |
838 | case 0b11000: word_stats[kCAEN1190Filler]++; break; |
839 | default: word_stats[kUnknown]++; break; |
840 | } |
841 | } |
842 | } |
843 | |
844 | |
845 | |
846 | |
847 | void JEventProcessor_DAQ_online::ParseModuleConfiguration(uint32_t rocid, uint32_t *&iptr, uint32_t *iend, uint32_t *word_stats) |
848 | { |
849 | while(iptr < iend){ |
850 | |
851 | word_stats[kConfig]++; |
852 | uint32_t Nvals = ((*iptr++) >> 24) & 0xFF; |
853 | |
854 | |
855 | for(uint32_t i=0; i< Nvals; i++){ |
856 | |
857 | switch((*iptr++)>>24){ |
858 | case 0x05: word_stats[kConfigf250]++; break; |
859 | case 0x0F: word_stats[kConfigf125]++; break; |
860 | case 0x06: word_stats[kConfigF1]++; break; |
861 | case 0x10: word_stats[kConfigCAEN1190]++; break; |
862 | default: word_stats[kConfig]++; break; |
863 | } |
864 | } |
865 | } |
866 | } |
867 | |
868 | |
869 | |
870 | |
871 | jerror_t JEventProcessor_DAQ_online::erun(void) |
872 | { |
873 | |
874 | |
875 | |
876 | return NOERROR; |
877 | } |
878 | |
879 | |
880 | |
881 | |
882 | jerror_t JEventProcessor_DAQ_online::fini(void) |
883 | { |
884 | |
885 | return NOERROR; |
886 | } |
887 | |