1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | #include <stdint.h> |
9 | #include <vector> |
10 | |
11 | #include <TMath.h> |
12 | |
13 | |
14 | #include "JEventProcessor_CDC_expert.h" |
15 | #include <JANA/JApplication.h> |
16 | |
17 | |
18 | using namespace std; |
19 | using namespace jana; |
20 | |
21 | |
22 | #include "CDC/DCDCHit.h" |
23 | #include "CDC/DCDCDigiHit.h" |
24 | #include "DAQ/Df125PulseIntegral.h" |
25 | #include "DAQ/Df125PulsePedestal.h" |
26 | #include "DAQ/Df125WindowRawData.h" |
27 | #include "DAQ/Df125CDCPulse.h" |
28 | |
29 | #include <TDirectory.h> |
30 | #include <TH2.h> |
31 | #include <TH1.h> |
32 | |
33 | |
34 | |
35 | |
36 | |
37 | |
38 | static TH1D *cdc_e = NULL__null; |
39 | static TH2D *cdc_e_vs_n = NULL__null; |
40 | |
41 | static TH1D *cdc_t = NULL__null; |
42 | static TH2D *cdc_t_vs_n = NULL__null; |
43 | |
44 | |
45 | static TH2D *cdc_e_ring[29]; |
46 | static TH2D *cdc_t_ring[29]; |
47 | |
48 | |
49 | static TH2D *cdc_e_vs_t; |
50 | static TH2D *cdc_e_vs_t_ring[29]; |
51 | |
52 | static TH2I *cdc_raw_int_vs_t; |
53 | static TH2I *cdc_raw_int_vs_t_ring[29]; |
54 | |
55 | |
56 | static TH2I *cdc_o_badt; |
57 | static TH2I *cdc_o_overflow; |
58 | |
59 | static TH2I *cdc_ped_ring[29]; |
60 | static TH1I *cdc_ped_badt; |
61 | static TH1I *cdc_ped_overflow; |
62 | |
63 | static TH2I *cdc_windata_ped_ring[29]; |
64 | |
65 | static TH2I *cdc_windata_ped_roc25; |
66 | static TH2I *cdc_windata_ped_roc26; |
67 | static TH2I *cdc_windata_ped_roc27; |
68 | static TH2I *cdc_windata_ped_roc28; |
69 | |
70 | static TH2I *cdc_raw_t_ring[29]; |
71 | static TH1I *cdc_raw_t_badt; |
72 | static TH1I *cdc_raw_t_overflow; |
73 | |
74 | static TH2I *cdc_raw_amp_ring[29]; |
75 | static TH1I *cdc_raw_amp_badt; |
76 | static TH1I *cdc_raw_amp_overflow; |
77 | |
78 | |
79 | |
80 | static TH2I *cdc_raw_intpp_ring[29]; |
81 | static TH1I *cdc_raw_intpp_badt; |
82 | static TH1I *cdc_raw_intpp_overflow; |
83 | |
84 | |
85 | |
86 | static TH2I *cdc_raw_int_ring[29]; |
87 | static TH1I *cdc_raw_int_badt; |
88 | static TH1I *cdc_raw_int_overflow; |
89 | |
90 | |
91 | |
92 | |
93 | |
94 | bool hists3675 = false; |
95 | |
96 | |
97 | |
98 | |
99 | |
100 | |
101 | extern "C"{ |
102 | void InitPlugin(JApplication *app){ |
103 | InitJANAPlugin(app); |
104 | app->AddProcessor(new JEventProcessor_CDC_expert()); |
105 | } |
106 | } |
107 | |
108 | |
109 | |
110 | |
111 | |
112 | JEventProcessor_CDC_expert::JEventProcessor_CDC_expert() { |
113 | } |
114 | |
115 | |
116 | |
117 | |
118 | |
119 | JEventProcessor_CDC_expert::~JEventProcessor_CDC_expert() { |
120 | } |
121 | |
122 | |
123 | |
124 | |
125 | jerror_t JEventProcessor_CDC_expert::init(void) { |
126 | |
127 | |
128 | |
129 | |
130 | |
131 | |
132 | |
133 | |
134 | return NOERROR; |
135 | } |
136 | |
137 | |
138 | |
139 | |
140 | |
141 | jerror_t JEventProcessor_CDC_expert::brun(JEventLoop *eventLoop, int32_t runnumber) { |
142 | |
143 | |
144 | |
145 | |
146 | japp->RootWriteLock(); |
147 | |
148 | |
149 | |
150 | |
151 | Int_t IMAX = 400000; |
152 | Int_t PMAX = 512; |
153 | Int_t AMAX = 4096; |
154 | Int_t RTMAX = 12000; |
155 | Int_t RTVSNMAX = 8192; |
156 | |
157 | Char_t rtunits[8] = "0.125ns"; |
158 | |
159 | if (runnumber > 3675) { |
160 | |
161 | |
162 | |
163 | |
164 | |
165 | |
166 | |
167 | |
168 | |
169 | |
170 | |
171 | |
172 | IMAX = 16384; |
173 | PMAX = 256; |
174 | AMAX = 512; |
175 | RTMAX = 2048; |
176 | RTVSNMAX = 1024; |
177 | |
178 | sprintf(rtunits,"0.8ns"); |
179 | |
180 | } |
181 | |
182 | |
183 | const Int_t EMAX = 21000; |
184 | |
185 | |
186 | |
187 | const Int_t TMAX = 2000; |
188 | |
189 | |
190 | |
191 | const Int_t NSTRAWS = 3522; |
192 | const Float_t HALF = 0.5; |
193 | const Float_t NSTRAWSPH = 3522.5; |
194 | |
195 | |
196 | |
197 | TDirectory *main = gDirectory(TDirectory::CurrentDirectory()); |
198 | gDirectory(TDirectory::CurrentDirectory())->mkdir("CDC_expert")->cd(); |
199 | |
200 | |
201 | |
202 | |
203 | |
204 | |
205 | const Int_t straws[29] = {0,42,42,54,54,66,66,80,80,93,93,106,106,123,123,135,135,146,146,158,158,170,170,182,182,197,197,209,209}; |
206 | |
207 | |
208 | cdc_e = new TH1D("cdc_e","CDC charge (fC);charge (fC)",200,0,EMAX); |
209 | cdc_e_vs_n = new TH2D("cdc_e_vs_n","CDC charge (fC) vs straw number;straw;charge (fC)",NSTRAWS,HALF,NSTRAWSPH,100,0,EMAX); |
210 | |
211 | cdc_e_vs_t = new TH2D("cdc_e_vs_t","CDC charge (fC) vs time (ns);time (ns);charge (fC)",150,-250,TMAX,100,0,EMAX); |
212 | |
213 | |
214 | cdc_t = new TH1D("cdc_t","CDC time (ns);time (ns)",300,-250,TMAX); |
215 | cdc_t_vs_n = new TH2D("cdc_t_vs_n","CDC time (ns) vs straw number;straw;time (ns)",NSTRAWS,HALF,NSTRAWSPH,150,-250,TMAX); |
216 | |
217 | |
218 | |
219 | |
220 | cdc_raw_int_vs_t = new TH2I("cdc_raw_int_vs_t",Form("CDC integral (ADC units), pedestal subtracted, vs raw time (units of %s);time (%s);integral, pedestal subtracted (ADC units)",rtunits,rtunits),(Int_t)256,0,RTVSNMAX,100,0,IMAX); |
221 | |
222 | |
223 | |
224 | cdc_windata_ped_roc25 = new TH2I("cdc_windata_ped_roc25","CDC pedestal (ADC units) from raw window data vs slot*100+channel, ROC 25;slot*100 + channel;pedestal",1600,200+HALF,1800+HALF,(Int_t)PMAX/4,0,PMAX); |
225 | cdc_windata_ped_roc26 = new TH2I("cdc_windata_ped_roc26","CDC pedestal (ADC units) from raw window data vs slot*100+channel, ROC 26;slot*100 + channel;pedestal",1600,200+HALF,1800+HALF,(Int_t)PMAX/4,0,PMAX); |
226 | cdc_windata_ped_roc27 = new TH2I("cdc_windata_ped_roc27","CDC pedestal (ADC units) from raw window data vs slot*100+channel, ROC 27;slot*100 + channel;pedestal",1600,200+HALF,1800+HALF,(Int_t)PMAX/4,0,PMAX); |
227 | cdc_windata_ped_roc28 = new TH2I("cdc_windata_ped_roc28","CDC pedestal (ADC units) from raw window data vs slot*100+channel, ROC 28;slot*100 + channel;pedestal",1600,200+HALF,1800+HALF,(Int_t)PMAX/4,0,PMAX); |
228 | |
229 | |
230 | |
231 | if (runnumber > 3675) { |
232 | |
233 | gDirectory(TDirectory::CurrentDirectory())->mkdir("bad_t","CDC Bad time flagged")->cd(); |
234 | |
235 | hists3675 = true; |
236 | |
237 | cdc_o_badt = new TH2I("cdc_o_badt","CDC occupancy by straw,ring, events with bad time flagged;straw;ring",209,0.5,209.5,28,0.5,28.5); |
238 | cdc_ped_badt = new TH1I("cdc_ped_badt","CDC pedestal, events with bad time flagged;straw;pedestal",256,0,PMAX); |
239 | cdc_raw_t_badt = new TH1I("cdc_raw_t_badt",Form("CDC raw time (units of %s), events with bad time flagged;straw;raw time (%s)",rtunits,rtunits),256,0,RTMAX); |
240 | cdc_raw_amp_badt = new TH1I("cdc_raw_amp_badt","CDC amplitude (ADC units), events with bad time flagged;ADC units",256,0,AMAX); |
241 | cdc_raw_int_badt = new TH1I("cdc_raw_intpp_badt","CDC integral (ADC units), pedestal subtracted, events with bad time flagged;ADC units",100,0,IMAX); |
242 | cdc_raw_intpp_badt = new TH1I("cdc_raw_intpp_badt","CDC integral (ADC units), including pedestal, events with bad time flagged;ADC units",100,0,IMAX); |
243 | |
244 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
245 | gDirectory(TDirectory::CurrentDirectory())->mkdir("overflows","CDC overflow flagged")->cd(); |
246 | |
247 | cdc_o_overflow = new TH2I("cdc_o_overflow","CDC overflow occupancy by straw,ring;straw;ring",209,0.5,209.5,28,0.5,28.5); |
248 | cdc_ped_overflow = new TH1I("cdc_ped_overflow","CDC pedestal, events with ADC overflow;pedestal",256,0,PMAX); |
249 | cdc_raw_t_overflow = new TH1I("cdc_raw_t_overflow",Form("CDC raw time (units of %s), events with ADC overflow;raw time (%s)",rtunits,rtunits),256,0,RTMAX); |
250 | cdc_raw_amp_overflow = new TH1I("cdc_raw_amp_overflow","CDC amplitude (ADC units), events with ADC overflow;ADC units",256,0,AMAX); |
251 | cdc_raw_int_overflow = new TH1I("cdc_raw_intpp_overflow","CDC integral (ADC units), pedestal subtracted, events with ADC overflow;ADC units",100,0,IMAX); |
252 | cdc_raw_intpp_overflow = new TH1I("cdc_raw_intpp_overflow","CDC integral (ADC units), including pedestal, events with ADC overflow;ADC units",100,0,IMAX); |
253 | |
254 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
255 | |
256 | } |
257 | |
258 | Int_t i; |
259 | |
260 | |
261 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_e_vs_t","CDC rings: charge vs time")->cd(); |
262 | |
263 | for (i=1; i<29; i++) { |
264 | cdc_e_vs_t_ring[i] = new TH2D(Form("cdc_e_vs_t_ring[%i]",i),"CDC charge (fC) vs time (ns);time (ns);charge (fC)",150,0,TMAX,100,0,EMAX); |
265 | } |
266 | |
267 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
268 | |
269 | |
270 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_int_vs_raw_t","CDC rings: integral vs raw time (pedestal subtracted)")->cd(); |
271 | |
272 | for (i=1; i<29; i++) { |
273 | cdc_raw_int_vs_t_ring[i] = new TH2I(Form("cdc_raw_int_vs_t_ring[%i]",i),Form("CDC integral (ADC units), pedestal subtracted, vs raw time (%s);raw time (%s);integral, pedestal subtracted (ADC units)",rtunits,rtunits),256,0,RTVSNMAX,100,0,IMAX); |
274 | |
275 | } |
276 | |
277 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
278 | |
279 | |
280 | |
281 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_e","CDC rings: charge vs straw")->cd(); |
282 | |
283 | for (i=1; i<29; i++) { |
284 | cdc_e_ring[i] = new TH2D(Form("cdc_e_ring[%i]",i),Form("CDC charge (fC), ring %i;straw;charge (fC)",i),straws[i],HALF,straws[i]+HALF,100,0,EMAX); |
285 | } |
286 | |
287 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
288 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings__t","CDC rings: time vs straw")->cd(); |
289 | |
290 | for (i=1; i<29; i++) { |
291 | cdc_t_ring[i] = new TH2D(Form("cdc_t_ring[%i]",i),Form("CDC time (ns), ring %i;straw;time (ns)",i),straws[i],HALF,straws[i]+HALF,150,0,TMAX); |
292 | } |
293 | |
294 | |
295 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
296 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_pedestal","CDC rings: pedestal vs straw")->cd(); |
297 | |
298 | for (i=1; i<29; i++) { |
299 | cdc_ped_ring[i] = new TH2I(Form("cdc_ped_ring[%i]",i),Form("CDC pedestal (ADC units), ring %i;straw;pedestal",i),straws[i],HALF,straws[i]+HALF,(Int_t)PMAX/2,0,PMAX); |
300 | } |
301 | |
302 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
303 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_windata_pedestal","CDC rings: pedestal from raw window data vs straw")->cd(); |
304 | |
305 | for (i=1; i<29; i++) { |
306 | cdc_windata_ped_ring[i] = new TH2I(Form("cdc_windata_ped_ring[%i]",i),Form("CDC pedestal (ADC units) from raw window data, ring %i;straw;pedestal",i),straws[i],HALF,straws[i]+HALF,(Int_t)PMAX/2,0,PMAX); |
307 | } |
308 | |
309 | |
310 | |
311 | |
312 | |
313 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
314 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_raw_t","CDC rings: raw time vs straw")->cd(); |
315 | |
316 | for (i=1; i<29; i++) { |
317 | cdc_raw_t_ring[i] = new TH2I(Form("cdc_raw_t_ring[%i]",i),Form("CDC raw time (units of %s), ring %i;straw;raw time (%s)",rtunits,i,rtunits),straws[i],HALF,straws[i]+HALF,256,0,RTVSNMAX); |
318 | } |
319 | |
320 | |
321 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
322 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_raw_amp","CDC rings: amplitude")->cd(); |
323 | |
324 | for (i=1; i<29; i++) { |
325 | cdc_raw_amp_ring[i] = new TH2I(Form("cdc_raw_amp_ring[%i]",i),Form("CDC amplitude (ADC units), ring %i",i),straws[i],HALF,straws[i]+HALF,256,0,AMAX); |
326 | } |
327 | |
328 | |
329 | |
330 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
331 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_raw_integral","CDC rings: integral vs straw (pedestal subtracted)")->cd(); |
332 | |
333 | for (i=1; i<29; i++) { |
334 | cdc_raw_int_ring[i] = new TH2I(Form("cdc_raw_int_ring[%i]",i),Form("CDC integral (ADC units), pedestal subtracted, ring %i",i),straws[i],HALF,straws[i]+HALF,100,0,IMAX); |
335 | } |
336 | |
337 | |
338 | gDirectory(TDirectory::CurrentDirectory())->cd("../"); |
339 | gDirectory(TDirectory::CurrentDirectory())->mkdir("rings_raw_integral_incl_ped","CDC rings: integral vs straw (including pedestal)")->cd(); |
340 | |
341 | for (i=1; i<29; i++) { |
342 | cdc_raw_intpp_ring[i] = new TH2I(Form("cdc_raw_intpp_ring[%i]",i),Form("CDC integral (ADC units), including pedestal, ring %i",i),straws[i],HALF,straws[i]+HALF,100,0,IMAX); |
343 | } |
344 | |
345 | |
346 | |
347 | |
348 | main->cd(); |
349 | |
350 | japp->RootUnLock(); |
351 | |
352 | return NOERROR; |
353 | } |
354 | |
355 | |
356 | |
357 | |
358 | |
359 | jerror_t JEventProcessor_CDC_expert::evnt(JEventLoop *eventLoop, uint64_t eventnumber) { |
360 | |
361 | |
362 | |
363 | |
364 | |
365 | |
366 | float q,t; |
367 | |
368 | uint32_t qf,ocount; |
369 | uint32_t tr,p,a; |
370 | uint32_t integral; |
371 | uint32_t integ; |
372 | |
373 | uint16_t ring,straw; |
374 | uint16_t n; |
375 | |
376 | Bool_t PED_SUB; |
377 | |
378 | uint32_t total_ped; |
379 | uint32_t nsamples_integral; |
380 | uint32_t nsamples_pedestal; |
381 | |
382 | uint32_t rocid; |
383 | uint32_t slot; |
384 | uint32_t channel; |
385 | |
386 | |
387 | const uint16_t NPEDSAMPLES=16; |
388 | |
389 | |
390 | int straw_offset[29] = {0,0,42,84,138,192,258,324,404,484,577,670,776,882,1005,1128,1263,1398,1544,1690,1848,2006,2176,2346,2528,2710,2907,3104,3313}; |
391 | |
392 | |
393 | |
394 | |
395 | |
396 | |
397 | vector<const DCDCHit*> hits; |
398 | eventLoop->Get(hits); |
399 | |
400 | |
401 | vector<const DCDCDigiHit*> digihits; |
402 | eventLoop->Get(digihits); |
403 | |
404 | |
405 | vector<const Df125WindowRawData*> wrdvector; |
406 | eventLoop->Get(wrdvector); |
407 | |
408 | japp->RootWriteLock(); |
409 | |
410 | |
411 | for(uint32_t i=0; i<hits.size(); i++) { |
412 | |
413 | const DCDCHit *hit = hits[i]; |
414 | |
415 | if(hit->q>0.0) { |
416 | |
417 | q = hit->q; |
418 | t = hit->t; |
419 | ring = hit->ring; |
420 | straw = hit->straw; |
421 | |
422 | n = straw_offset[ring] + straw; |
423 | |
424 | if (q > 0.0) { |
425 | cdc_e->Fill(q); |
426 | cdc_e_vs_n->Fill(n,q); |
427 | } |
428 | |
429 | |
430 | cdc_t->Fill(t); |
431 | cdc_t_vs_n->Fill(n,t); |
432 | |
433 | |
434 | cdc_e_vs_t->Fill(t,q); |
435 | cdc_e_vs_t_ring[ring]->Fill(t,q); |
436 | |
437 | cdc_e_ring[ring]->Fill(straw,q); |
438 | cdc_t_ring[ring]->Fill(straw,t); |
439 | } |
440 | } |
441 | |
442 | |
443 | for(uint32_t i=0; i<digihits.size(); i++) { |
444 | |
445 | const DCDCDigiHit *digihit = digihits[i]; |
446 | |
447 | |
448 | const Df125PulseIntegral *pi = NULL__null; |
449 | const Df125PulsePedestal *pp = NULL__null; |
450 | const Df125WindowRawData *windat = NULL__null; |
451 | const Df125CDCPulse *cp = NULL__null; |
452 | |
453 | |
454 | vector<uint16_t> samples; |
455 | uint32_t winped=0; |
456 | |
457 | PED_SUB = kFALSE; |
458 | total_ped = 0; |
| Value stored to 'total_ped' is never read |
459 | |
460 | rocid = 0; |
461 | slot = 0; |
462 | channel = 0; |
463 | qf = 0; |
464 | ocount = 0; |
465 | a = 0; |
466 | |
467 | |
468 | digihit->GetSingle(pi); |
469 | if (pi) { |
470 | rocid = pi->rocid; |
471 | slot = pi->slot; |
472 | channel = pi->channel; |
473 | pi->GetSingle(windat); |
474 | } else if (i < (uint32_t)wrdvector.size()) { |
475 | windat = wrdvector[i]; |
476 | } |
477 | |
478 | nsamples_integral = pi ? pi->nsamples_integral : 0; |
479 | nsamples_pedestal = pi ? pi->nsamples_pedestal : 0; |
480 | |
481 | if ((nsamples_integral > 0) && (nsamples_pedestal > 0)) PED_SUB = kTRUE; |
482 | |
483 | digihit->GetSingle(pp); |
484 | if(pp) a = pp->pulse_peak; |
485 | |
486 | |
487 | digihit->GetSingle(cp); |
488 | if (cp) { |
489 | rocid = cp->rocid; |
490 | slot = cp->slot; |
491 | channel = cp->channel; |
492 | a = cp->first_max_amp; |
493 | qf = cp->time_quality_bit; |
494 | ocount = cp->overflow_count; |
495 | } |
496 | |
497 | |
498 | |
499 | if (windat) { |
500 | |
501 | if (windat->samples.size()>=NPEDSAMPLES) { |
502 | |
503 | winped = 0; |
504 | |
505 | for (uint16_t j=0; j<NPEDSAMPLES; j++) winped += (uint32_t)windat->samples[j]; |
506 | |
507 | winped = (uint32_t)winped/16.0; |
508 | |
509 | if (winped > 0) { |
510 | |
511 | if (rocid == 25) cdc_windata_ped_roc25->Fill(100*slot + channel,winped); |
512 | if (rocid == 26) cdc_windata_ped_roc26->Fill(100*slot + channel,winped); |
513 | if (rocid == 27) cdc_windata_ped_roc27->Fill(100*slot + channel,winped); |
514 | if (rocid == 28) cdc_windata_ped_roc28->Fill(100*slot + channel,winped); |
515 | |
516 | } |
517 | |
518 | } |
519 | } |
520 | |
521 | |
522 | |
523 | |
524 | if((digihit->pulse_integral>0)||(digihit->pulse_time>0)) { |
525 | |
526 | ring = digihit->ring; |
527 | straw = digihit->straw; |
528 | |
529 | p = digihit->pedestal; |
530 | tr = digihit->pulse_time; |
531 | integral = digihit->pulse_integral; |
532 | |
533 | integ = 0; |
534 | |
535 | |
536 | |
537 | if (PED_SUB) { |
538 | total_ped = p*nsamples_integral/nsamples_pedestal; |
539 | integ = integral - total_ped; |
540 | } |
541 | |
542 | |
543 | |
544 | if (PED_SUB) cdc_raw_int_vs_t->Fill(tr,integ); |
545 | if (PED_SUB) cdc_raw_int_vs_t_ring[ring]->Fill(tr,integ); |
546 | |
547 | cdc_ped_ring[ring]->Fill(straw,p); |
548 | cdc_raw_t_ring[ring]->Fill(straw,tr); |
549 | cdc_raw_amp_ring[ring]->Fill(straw,a); |
550 | if (PED_SUB) cdc_raw_int_ring[ring]->Fill(straw,integ); |
551 | cdc_raw_intpp_ring[ring]->Fill(straw,integral); |
552 | if (winped) cdc_windata_ped_ring[ring]->Fill(straw,winped); |
553 | |
554 | |
555 | if (cp && hists3675) { |
556 | if (qf==1) { |
557 | cdc_o_badt->Fill(straw,ring); |
558 | cdc_ped_badt->Fill(p); |
559 | cdc_raw_t_badt->Fill(tr); |
560 | cdc_raw_amp_badt->Fill(a); |
561 | if (PED_SUB) cdc_raw_int_badt->Fill(integ); |
562 | cdc_raw_intpp_badt->Fill(integral); |
563 | } |
564 | |
565 | |
566 | if (ocount>0) { |
567 | cdc_o_overflow->Fill(straw,ring); |
568 | cdc_ped_overflow->Fill(p); |
569 | cdc_raw_t_overflow->Fill(tr); |
570 | cdc_raw_amp_overflow->Fill(a); |
571 | if (PED_SUB) cdc_raw_int_overflow->Fill(integ); |
572 | cdc_raw_intpp_overflow->Fill(integral); |
573 | } |
574 | } |
575 | |
576 | |
577 | } |
578 | |
579 | } |
580 | |
581 | |
582 | japp->RootUnLock(); |
583 | |
584 | |
585 | return NOERROR; |
586 | } |
587 | |
588 | |
589 | |
590 | |
591 | |
592 | jerror_t JEventProcessor_CDC_expert::erun(void) { |
593 | |
594 | |
595 | |
596 | return NOERROR; |
597 | } |
598 | |
599 | |
600 | |
601 | |
602 | |
603 | jerror_t JEventProcessor_CDC_expert::fini(void) { |
604 | |
605 | return NOERROR; |
606 | } |
607 | |
608 | |
609 | |
610 | |