Changed names of static libraries, targets and updated README_CMake
[u/mrichter/AliRoot.git] / TRD / AliTRDqaBlackEvents.cxx
CommitLineData
48ccaa07 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * withount fee, provided thats the abov copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id: AliTRDqaBlackEvents.cxx 23387 2008-01-17 17:25:16Z cblume $ */
17
18////////////////////////////////////////////////////////////////////////////
19// //
20// QA of black events //
21// //
22// Author: //
23// Sylwester Radomski (radomski@physi.uni-heidelberg.de) //
24// //
25////////////////////////////////////////////////////////////////////////////
26
27#include "TH1D.h"
28#include "TH2D.h"
29#include "TH2S.h"
48ccaa07 30#include "TF1.h"
31#include "TFile.h"
32#include "TCanvas.h"
33#include "TPad.h"
34#include "TLatex.h"
35#include "TStyle.h"
36#include "TGraph.h"
37
c87495af 38#include "AliRawReader.h"
48ccaa07 39
ef1860eb 40#include "AliTRDrawStreamOld.h"
48ccaa07 41#include "AliTRDqaBlackEvents.h"
42
48ccaa07 43ClassImp(AliTRDqaBlackEvents)
44
45///////////////////////////////////////////////////////////////////////////////////////////////////
46
47AliTRDqaBlackEvents::AliTRDqaBlackEvents()
48 :TObject()
49 ,fnEvents(0)
50 ,fCreateFull(0)
51 ,fThresh(0)
52 ,fCount(0)
53 ,fRefEv(0)
54 ,fOccupancy(0)
55 ,fDetRob(0)
56 ,fTBEvent(0)
57 ,fRefHistPed(0)
58 ,fRefHistNoise(0)
59 ,fErrorHC(0)
60 ,fErrorMCM(0)
61 ,fErrorADC(0)
62 ,fErrorSMHC(0)
63 ,fErrorSMMCM(0)
64 ,fErrorSMADC(0)
65 ,fErrorGraphHC(0)
66 ,fErrorGraphMCM(0)
67 ,fErrorGraphADC(0)
68 ,fGraphMCM(0)
69 ,fMcmTracks(0)
70 ,fMapMCM(0)
71 ,fFracMCM(0)
72 ,fSMHCped(0)
73 ,fSMHCerr(0)
74 ,fNoiseTotal(0)
75 ,fPP(0)
76 ,fMinNoise(0.5)
77 ,fMaxNoise(2)
78 ,fFitType(0)
48ccaa07 79{
80 //
81 // Constructor
82 // to create the histograms call Init()
83 //
84
8d3a4626 85 for (Int_t i = 0; i < kDET; i++) {
86 fPed[i] = 0x0;
87 fNoise[i] = 0x0;
88 fChPP[i] = 0x0;
89 fNPointDist[i] = 0x0;
90 fChPed[i] = 0x0;
91 fChNoise[i] = 0x0;
92 fNPoint[i] = 0x0;
93 fData[i] = 0x0;
94 fSignal[i] = 0x0;
95 fnEntriesRM[i] = 0x0;
96 fnEntriesRMDist[i] = 0x0;
97 fChPedRes[i] = 0x0;
98 fChNoiseRes[i] = 0x0;
99 fErrorLocHC[i] = 0x0;
100 fErrorLocMCM[i] = 0x0;
101 fErrorLocADC[i] = 0x0;
102 }
103 for (Int_t i = 0; i < 3; i++) {
104 fGraphPP[i] = 0x0;
105 fSMLink[i] = 0x0;
106 fGrLink[i] = 0x0;
107 fppThresh[i] = 0;
108 fnPP[i] = 0;
109 fnLink[i] = 0;
110 }
111 for (Int_t i = 0; i < 2; i++) {
112 fnErrorHC[i] = 0;
113 fnErrorMCM[i] = 0;
114 fnErrorADC[i] = 0;
115 }
116 for (Int_t i = 0; i < kSM; i++) {
117 fSmNoiseRms[i] = 0x0;
118 fSmNoiseFit[i] = 0x0;
119 fSmPP[i] = 0x0;
120 }
121 for (Int_t i = 0; i < kSM+1; i++) {
122 fNumberADC[i] = 0x0;
123 fnADCinSM[i] = 0;
124 }
125 for (Int_t i = 0; i < 1000; i++) {
126 fEvNoDist[i] = 0;
127 }
2f9bdd85 128 strncpy(fRefFileName,"",256);
8d3a4626 129
48ccaa07 130}
131
132///////////////////////////////////////////////////////////////////////////////////////////////////
133
134AliTRDqaBlackEvents::AliTRDqaBlackEvents(const AliTRDqaBlackEvents &qa)
135 :TObject(qa)
136 ,fnEvents(0)
137 ,fCreateFull(0)
138 ,fThresh(0)
139 ,fCount(0)
140 ,fRefEv(0)
141 ,fOccupancy(0)
142 ,fDetRob(0)
143 ,fTBEvent(0)
144 ,fRefHistPed(0)
145 ,fRefHistNoise(0)
146 ,fErrorHC(0)
147 ,fErrorMCM(0)
148 ,fErrorADC(0)
149 ,fErrorSMHC(0)
150 ,fErrorSMMCM(0)
151 ,fErrorSMADC(0)
152 ,fErrorGraphHC(0)
153 ,fErrorGraphMCM(0)
154 ,fErrorGraphADC(0)
155 ,fGraphMCM(0)
156 ,fMcmTracks(0)
157 ,fMapMCM(0)
158 ,fFracMCM(0)
159 ,fSMHCped(0)
160 ,fSMHCerr(0)
161 ,fNoiseTotal(0)
162 ,fPP(0)
163 ,fMinNoise(0.5)
164 ,fMaxNoise(2)
165 ,fFitType(0)
48ccaa07 166{
167 //
168 // Copy constructor
169 // to create the histograms call Init()
170 //
171
8d3a4626 172 for (Int_t i = 0; i < kDET; i++) {
173 fPed[i] = 0x0;
174 fNoise[i] = 0x0;
175 fChPP[i] = 0x0;
176 fNPointDist[i] = 0x0;
177 fChPed[i] = 0x0;
178 fChNoise[i] = 0x0;
179 fNPoint[i] = 0x0;
180 fData[i] = 0x0;
181 fSignal[i] = 0x0;
182 fnEntriesRM[i] = 0x0;
183 fnEntriesRMDist[i] = 0x0;
184 fChPedRes[i] = 0x0;
185 fChNoiseRes[i] = 0x0;
186 fErrorLocHC[i] = 0x0;
187 fErrorLocMCM[i] = 0x0;
188 fErrorLocADC[i] = 0x0;
189 }
190 for (Int_t i = 0; i < 3; i++) {
191 fGraphPP[i] = 0x0;
192 fSMLink[i] = 0x0;
193 fGrLink[i] = 0x0;
194 fppThresh[i] = 0;
195 fnPP[i] = 0;
196 fnLink[i] = 0;
197 }
198 for (Int_t i = 0; i < 2; i++) {
199 fnErrorHC[i] = 0;
200 fnErrorMCM[i] = 0;
201 fnErrorADC[i] = 0;
202 }
203 for (Int_t i = 0; i < kSM; i++) {
204 fSmNoiseRms[i] = 0x0;
205 fSmNoiseFit[i] = 0x0;
206 fSmPP[i] = 0x0;
207 }
208 for (Int_t i = 0; i < kSM+1; i++) {
209 fNumberADC[i] = 0x0;
210 fnADCinSM[i] = 0;
211 }
212 for (Int_t i = 0; i < 1000; i++) {
213 fEvNoDist[i] = 0;
214 }
215 strncpy(fRefFileName,"",256);
216
48ccaa07 217}
218
219///////////////////////////////////////////////////////////////////////////////////////////////////
220
221void AliTRDqaBlackEvents::Init()
222{
223 //
224 // creates histograms
225 //
226
227 //TFile *file = new
228 //Info("Init", "Statring");
229
230 fnEvents = 0;
231
232 // histograms for chambers
233 for(Int_t det=0; det<kDET; det++) {
234
235 fNPoint[det] = new TH2D(Form("entries_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
236 //fData[det] = new TH3F(Form("data_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5, 50, -0.5, 49.5);
237
238 // pedestal noise maps using RMS and Fit
239 fChPed[det] = new TH2D(Form("ped_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
240 fChNoise[det] = new TH2D(Form("noise_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
241
242 //fChPed[det] = new TH2D(Form("ped_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
243 //fChNoise[det] = new TH2D(Form("noise_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
244
245 // distribution per detector
246 fPed[det] = new TH1D(Form("pedDist_%d", det), ";pedestals (ADC counts)", 100, 5, 15);
247 fNoise[det] = new TH1D(Form("noiseDist_%d", det), ";noise (ADC counts)", 100, 0, 5);
248 fSignal[det] = new TH1D(Form("signal_%d", det), ";signal (ADC counts)", 100, -0.5, 99.5);
249 fChPP[det] = new TH1D(Form("pp_%d", det), ";pp (ADC)", 200, -0.5, 199.5);
250
251 fnEntriesRM[det] = new TH2D(Form("entriesRM_%d", det), ";ROB,MCM", 8, -0.5, 7.5, 16, -0.5, 15.5);
252
253 // histograms after reference subtraction
254 fChPedRes[det] = new TH2D(Form("pedRef_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
255 fChNoiseRes[det] = new TH2D(Form("noiseRef_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
256
257
258 // error codes
259 fErrorLocMCM[det] = new TH2D(Form("errorLocMCM_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
260 fErrorLocADC[det] = new TH2D(Form("errorLocADC_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
261 fErrorLocHC[det] = new TH2D(Form("errorLocHC_%d", det), "", 16, -0.5, 15.5, 144, -0.5, 143.5);
262 }
263
264 // histogram for each MCM
265 for(Int_t i=0; i < kDET * kROB * kMCM; i++)
266 fFullCounter[i] = 0;
267
268 // histograms from the whole detector
269 fOccupancy = new TH1D("occupancy", "", 20, -0.5, 19.5);
270 fDetRob = new TH2D("DetRob", ";detector;ROB", kDET, -0.5, 539.5, 8, -0.5, 7.5);
271 fTBEvent = new TH2D("tbEvent", ";event ID;time bin", 100, -0.5, 99.5, 30, -0.5, 29.5);
272
273 // errors statistics and location
274 fErrorHC = new TH1D("errorHC", ";error ID;", 18, -3.5, 14.5);
275 fErrorMCM = new TH1D("errorMCM", ";error ID;", 18, -3.5, 14.5);
276 fErrorADC = new TH1D("errorADC", ";error ID;", 18, -3.5, 14.5);
277
278 fErrorSMHC = new TH1D("errorSM_HC", ";SM id", 18, -0.5, 17.5);
279 fErrorSMMCM = new TH1D("errorSM_MCM", ";SM id", 18, -0.5, 17.5);
280 fErrorSMADC = new TH1D("errorSM_ADC", ";SM id", 18, -0.5, 17.5);
281
282
283 fErrorGraphHC = new TGraph();
284 fErrorGraphMCM = new TGraph();
285 fErrorGraphADC = new TGraph();
286
287 fGraphMCM = new TGraph();
288
289 for(Int_t i=0; i<3; i++) {
290 fGraphPP[i] = new TGraph();
291 }
292
293
294 fMapMCM = new TH2D("mapMCM", ";det;mcm", 540, -0.5, 539.5, kROB*kMCM, -0.5, kROB*kMCM-0.5);
295 fFracMCM = new TH1D("fracMCM", ";frequency", 100, 0, 1);
296
297
298 fErrorGraphHC->GetHistogram()->SetTitle("fraction of events with HC error;event number");
299 fErrorGraphMCM->GetHistogram()->SetTitle("fraction of events with MCM error;event number;");
300 fErrorGraphADC->GetHistogram()->SetTitle("fraction of events with ADC error;event number;");
301
302
303 fSMHCped = new TH2D("smHcPed", ";super module;half chamber", 18, -0.5, 17.5, 60, -0.5, 59.5);
304
305 // link monitor
306 const char *linkName[3] = {"smLink", "smBeaf", "smData"};
307 const char *linkGrName[3] = {"grSmLink", "grSmBeaf", "grSmData"};
308 for(Int_t i=0; i<3; i++) {
309 fSMLink[i] = new TH2D(linkName[i], ";super module;link", 18, -0.5, 17.5, 60, -0.5, 59.5);
310 fGrLink[i] = new TGraph();
311 fGrLink[i]->SetName(linkGrName[i]);
312 }
313
314 //fZSsize = new TH1D("zssizeSingle", ";threshold;nADC", 40, -0.5, 39.5);
315
316
317 //Info("Init", "Done");
318
319 // number of ADC channels fired per SM and in total
320 for(Int_t sm=0; sm<kSM+1; sm++)
321 fNumberADC[sm] = new TGraph();
322
323 //
324 fNoiseTotal = new TH1D("noiseTotal", "noise (ADC)", 250, 0, 10);
325 fPP = new TH1D("peakPeak", "p-p (ADC)", 200, -0.5, 199.5);
326
327 for(Int_t sm=0; sm<kSM; sm++) {
328 fSmNoiseRms[sm] = new TH1D(Form("noiseRms_sm%d", sm), ";noise from RMS (ADC)", 100, 0, 10);
329 fSmNoiseFit[sm] = new TH1D(Form("noiseFit_sm%d", sm), ";noise frim Fit (ADC)", 100, 0, 10);
330 fSmPP[sm] = new TH1D(Form("peakPeak_sm%d", sm), ";peak-peak (ADC)", 200, -0.5, 199.5);
331 }
332
333 // event number consistancy
334 for(Int_t i=0; i<1000; i++) {
335 fEvNoDist[i] = new TH1D(Form("mcmEvDist_%d", i), ";#Delta Events", 201, -100.5, 100.5);
336 }
337 fRefEv = -1;
338
339 fMcmTracks = new TObjArray();
340
341 // clean data direct
342 for(Int_t det=0; det<kDET; det++)
343 for(Int_t row=0; row<kROW; row++)
344 for(Int_t pad=0; pad<kPAD; pad++)
345 for(Int_t ch=0; ch<kCH; ch++) {
346 fDataDirect[det][row][pad][ch] = 0;
347 fSignalDirect[det][ch] = 0; // overdone
348 }
349}
350
351
352///////////////////////////////////////////////////////////////////////////////////////////////////
353
354void AliTRDqaBlackEvents::Reset()
355{
356 //
357 // Resets the histograms
358 //
359
360 for(Int_t i=0; i<kDET; i++) {
361 //fData[i]->Reset();
362 fChPed[i]->Reset();
363 fChNoise[i]->Reset();
364 }
365}
366
367
368///////////////////////////////////////////////////////////////////////////////////////////////////
369
370void AliTRDqaBlackEvents::SetRefFile(const char *filename) {
371
8d3a4626 372 strncpy(fRefFileName,filename,256);
373
48ccaa07 374}
375
376///////////////////////////////////////////////////////////////////////////////////////////////////
377
378void AliTRDqaBlackEvents::ReadRefHists(Int_t det) {
c87495af 379 //
380 // Read the reference histograms
381 //
382
48ccaa07 383 fRefHistPed = 0;
384 fRefHistNoise = 0;
385
386 TFile *file = 0;
387 if (fRefFileName) TFile::Open(fRefFileName);
388 if (!file) return;
389
390 fRefHistPed = (TH2D*)file->Get(Form("ped_%d",det));
391 fRefHistNoise = (TH2D*)file->Get(Form("noise_%d", det));
392
393 if (file) file->Close();
394}
395
396///////////////////////////////////////////////////////////////////////////////////////////////////
397
398void AliTRDqaBlackEvents::StartEvent()
399{
400 //
401 // start an event
402 //
403
404 // clear the mcm data
405 for(Int_t i=0; i < kDET * kROB * kMCM; i++) {
406 if (fFullSignal[i]) fFullSignal[i]->Reset();
407 fFullCounter[i] = 0;
408 }
409
410 for(Int_t i=0; i<2; i++) {
411 fnErrorHC[i] = 0;
412 fnErrorMCM[i] = 0;
413 fnErrorADC[i] = 0;
414 }
415
416
417 Int_t ppThresh[3] = {10, 20, 40};
418 for(Int_t i=0; i<3; i++) {
419 fppThresh[i] = ppThresh[i];
420 fnPP[i] = 0;
421 fnLink[i] = 0;
422 }
423
424 for(Int_t sm=0; sm<kSM+1; sm++) fnADCinSM[sm] = 0;
425
426 if (fRefEv > 0) fRefEv++;
427 fEvNoDist[999]->Reset(); // keep only the last event
428
429}
430
431///////////////////////////////////////////////////////////////////////////////////////////////////
432
ef1860eb 433void AliTRDqaBlackEvents::AddBuffer(AliTRDrawStreamOld *data, AliRawReader * const reader)
48ccaa07 434{
435
436 //printf ("try to read data\n");
437 Int_t nextBuff = data->NextBuffer();
438 //printf("done ...\n");
439
440 if (nextBuff == 0) return;
441
442 Int_t sm = reader->GetEquipmentId() - 1024;
443 //printf("reading SM %d\n", sm);
444 AliInfo(Form("reading SM %d", sm));
445
446 if (sm < 0 || sm > 18) return;
447
448 // lopp over stacks, links ...
449
450 for (Int_t istack = 0; istack < 5; istack++) {
451 for (Int_t ilink = 0; ilink < 12; ilink++) {
452
453 //printf("HC = %d %d\n", istack, ilink);
454
455 Int_t det = sm * 30 + istack * 6 + ilink/2;
456
457 // check if data delivered
458 if (!(data->IsLinkActiveInStack(istack, ilink))) continue;
459 fSMLink[0]->Fill(sm, istack * 12 + ilink);
460 fnLink[0]++;
461
462 // check if beaf-beaf
463 if (data->GetLinkMonitorError(istack, ilink)) {
464 fSMLink[1]->Fill(sm, istack * 12 + ilink);
465 fnLink[1]++;
466 continue;
467 }
468
469 // fill histogram with HC header errors
470 Int_t nErrHc = 0;
471 Int_t nErrHcTot = 0;
472
473 nErrHc = FillBits(fErrorHC, data->GetH0ErrorCode(istack, ilink), 0);
474 if (!nErrHc) fErrorHC->Fill(-3);
475 nErrHcTot += nErrHc;
476
477 nErrHc = FillBits(fErrorHC, data->GetH1ErrorCode(istack, ilink), 2);
478 if (!nErrHc) fErrorHC->Fill(-2);
479 nErrHcTot += nErrHc;
480
481 nErrHc = FillBits(fErrorHC, data->GetHCErrorCode(istack, ilink), 4);
482 if (!nErrHc) fErrorHC->Fill(-1);
483 nErrHcTot += nErrHc;
484
485 // trending
486 fnErrorHC[0]++;
487 if (nErrHcTot > 0) {
488 fnErrorHC[1]++;
489 fErrorSMHC->Fill(sm);
490 }
491
492 // data integrity protection
493
494 //if (data->GetHCErrorCode(istack, ilink) > 0) continue;
495 if (data->GetH0ErrorCode(istack, ilink) > 0) continue;
496 if (data->GetH1ErrorCode(istack, ilink) > 0) continue;
497
498 fSMLink[2]->Fill(sm, istack * 12 + ilink);
499 fnLink[2]++;
500
501
502 for (Int_t imcm = 0; imcm < data->GetHCMCMmax(istack, ilink); imcm++ ){
503
504 //printf("mcm = %d %d %d\n", istack, ilink, imcm);
505
506 // fill MCM error code
507
508 Int_t nErrMcm = 0;
509 Int_t nErrMcmTot = 0;
510
511 nErrMcm = FillBits(fErrorMCM, data->GetMCMhdErrorCode(istack, ilink, imcm), 0);
512 if (!nErrMcm) fErrorMCM->Fill(-3);
513 nErrMcmTot += nErrMcm;
514
515 nErrMcm = FillBits(fErrorMCM, data->GetMCMADCMaskErrorCode(istack, ilink, imcm), 5);
516 if (!nErrMcm) fErrorMCM->Fill(-2);
517 nErrMcmTot += nErrMcm;
518
519 nErrMcm = FillBits(fErrorMCM, data->GetMCMErrorCode(istack, ilink, imcm), 10);
520 if (!nErrMcm) fErrorMCM->Fill(-1);
521 nErrMcmTot += nErrMcm;
522
523 // trending
524 fnErrorMCM[0]++;
525 if (nErrMcmTot > 0) {
526 fnErrorMCM[1]++;
527 fErrorSMMCM->Fill(sm);
528 }
529
530 // MCM protection
531 if ( (data->GetMCMhdErrorCode(istack,ilink,imcm)) & 2 ) continue;
532 //if ((data->GetMCMADCMaskErrorCode(istack,ilink,imcm))) continue;
533 //if ((data->GetMCMErrorCode(istack,ilink,imcm))) continue;
534
535 Int_t mcmEvent = data->GetEventNumber(istack, ilink, imcm);
536
537 // set the reference event number
538 if (fRefEv < 0) {
539 fRefEv = mcmEvent;
540 printf("Reference Event Number = %d (%d %d %d)\n", fRefEv, istack, ilink, imcm);
541 }
542
543 // fill event distribution
544 if (!(fnEvents%10)) {
545 fEvNoDist[fnEvents/10]->Fill(mcmEvent - fRefEv);
546 }
547
548 fEvNoDist[999]->Fill(mcmEvent - fRefEv);
549
550 Int_t mcm = data->GetMCM(istack, ilink, imcm);
551 Int_t rob = data->GetROB(istack, ilink, imcm);
552
553 // create a structure for an MCM if needed
554 Int_t mcmIndex = det * (kMCM * kROB) + rob * kMCM + mcm;
555 if (fCreateFull && !fFullSignal[mcmIndex])
556 fFullSignal[mcmIndex] =
557 new TH2S(Form("mcm_%d_%d_%d_%d_%d", sm, istack, ilink/2, rob, mcm),
558 Form("mcm-%d-%d-%d-%d-%d;ADC;time bin", sm, istack, ilink/2, rob, mcm),
559 21, -0.5, 20.5, 30, -0.5, 29.5);
560
561
562 //Int_t zsADC[21][40];
563 /*
564 for(Int_t ina=0; ina<21; ina++)
565 for(Int_t th=0; th<40; th++)
566 zsADC[ina][th] = 0;
567 */
568
569 // first loop over ADC chanels
570
571 for (Int_t iadc=0; iadc < data->GetADCcount(istack, ilink, imcm); iadc++) {
572
573 //printf("ADC = %d\n", iadc);
574
575
576 // fill ADC error bits
577 Int_t nErrAdc = FillBits(fErrorADC, data->GetADCErrorCode(), 0);
578 if (!nErrAdc) fErrorADC->Fill(-1);
579
580 fnErrorADC[0]++;
581 if (nErrAdc > 0) {
582 fnErrorADC[1]++;
583 fErrorSMADC->Fill(sm);
584 }
585
586 // ADC protection
587 if ((data->GetADCErrorCode(istack,ilink,imcm,iadc))) continue;
588
589 Int_t minV = 1024;
590 Int_t maxV = 0;
591
592 Int_t *sig = data->GetSignalDirect(istack, ilink, imcm, iadc);
593
594 //Int_t adc = data->GetADCnumber(istack, ilink, imcm, iadc);
595 Int_t row = data->GetRow(istack, ilink, imcm);
596 Int_t col = data->GetCol(istack, ilink, imcm, iadc);
597
598 // loop over Time Bins and fill histograms
599 for(Int_t k=0; k < data->GetNumberOfTimeBins(istack, ilink); k++) {
600
601 //fSignal[det]->Fill(sig[k]);
602 //fData[det]->Fill(row, col, sig[k]); // slow
603
604 if (sig[k] < kCH) {
605 fSignalDirect[det][sig[k]]++;
606 fDataDirect[det][row][col][sig[k]]++; // direct data
607 }
608
609 // peak-peak
610 minV = (minV < sig[k]) ? minV : sig[k];
611 maxV = (maxV > sig[k]) ? maxV : sig[k];
612
613 // check for active MCMs
614 if (fCreateFull && fFullSignal[mcmIndex]) {
615 if (sig[k] > fThresh || sig[k] < 0) fFullCounter[mcmIndex]++;
616 //if (sm == 0 && istack == 0 && ilink/2 == 1 && rob == 1 && mcm == 15) fFullCounter[mcmIndex]++; // special
617 //fFullSignal[mcmIndex]->Fill(adc, k, sig[k]); // slow
618 }
619
620 // zero suppresion tests
621 /*
622 for(Int_t th=0; th<40; th++) {
623 if (sig[k] > th) {
624 zsADC[iadc][th] = 1;
625 if (iadc > 0) zsADC[iadc-1][th] = 1;
626 if (iadc < 10) zsADC[iadc+1][th] = 1;
627 }
628 }
629 */
630
631 } // tb
632
633 if (maxV > 0) {
634 fnADCinSM[sm]++;
635 fnADCinSM[kSM]++;
636 }
637
638 Int_t adcPP = maxV - minV;
639 //if (adcPP == 100) fFullCounter[mcmIndex] += 10;
640
641 fPP->Fill(adcPP);
642 fChPP[det]->Fill(adcPP);
643 fSmPP[sm]->Fill(adcPP);
644
645 for(Int_t i=0; i<3; i++) {
646 if ((adcPP) > fppThresh[i]) fnPP[i]++;
647 }
648
649
650 } // adc
651
652 // fill ZS histos
653 /*
654 for(Int_t th=0; th<40; th++) {
655 Int_t nnADC = 0;
656 for(Int_t ins=0; ins<21; ins++)
657 nnADC += zsADC[ins][th];
658 fZSsize->Fill(th, nnADC);
659 }
660 */
661
662 // fill active MCMs
663
664 if (fCreateFull && fFullSignal[mcmIndex] && (fFullCounter[mcmIndex] > fCount)) {
665
666 for (Int_t iadc=0; iadc < data->GetADCcount(istack, ilink, imcm); iadc++) {
667
668 // ADC protection
669 if ((data->GetADCErrorCode(istack,ilink,imcm,iadc))) continue;
670
671 //Int_t row = data->GetRow(istack, ilink, imcm);
672 //Int_t col = data->GetCol(istack, ilink, imcm, iadc);
673 Int_t adc = data->GetADCnumber(istack, ilink, imcm, iadc);
674
675 Int_t *sig = data->GetSignalDirect(istack, ilink, imcm, iadc);
676
677 // loop over Time Bins and fill histograms
678 for(Int_t k=0; k < data->GetNumberOfTimeBins(istack, ilink); k++) {
679 fFullSignal[mcmIndex]->Fill(adc, k, sig[k]); // slow
680 }
681 } // tb
682 }
683
684 } // mcm
685 } // link
686 } // stack
687
688 // printf("end of loops\n");
689}
690
691
692///////////////////////////////////////////////////////////////////////////////////////////////////
693
694void AliTRDqaBlackEvents::FinishEvent()
695{
c87495af 696 //
697 // Processing at the end of the current event
698 //
48ccaa07 699
700 for(Int_t i=0; i<3; i++) {
701 fGraphPP[i]->SetPoint(fnEvents, fnEvents, fnPP[i]);
702 }
703
704 // trend of the number of links
705 for(Int_t i=0; i<3; i++) {
706 fGrLink[i]->SetPoint(fnEvents, fnEvents, fnLink[i]);
707 }
708
709 // save interesting histos
710 Int_t mcmTrackCandidate = 0;
711 for(Int_t i = 0; i < kDET * kROB * kMCM; i++) {
712 if ((fFullCounter[i] > fCount) && fFullSignal[i] && CheckMCM(i) ) {
713
714 fMcmTracks->AddLast(fFullSignal[i]->Clone(Form("event_%d_%s", fnEvents, fFullSignal[i]->GetName())));
715 mcmTrackCandidate++;
716
717 Int_t mcmTrackletDet = i/(kROB * kMCM);
718 Int_t mcmTrackletMcm = i%(kROB * kMCM);
719 fMapMCM->Fill(mcmTrackletDet, mcmTrackletMcm);
720 }
721 }
722
723 fGraphMCM->SetPoint(fnEvents, fnEvents, mcmTrackCandidate);
724 AliInfo(Form("Number of MCM track candidates = %d\n", mcmTrackCandidate));
725
726
727 // update fraction of error graphs
728 Double_t err;
729
730 err = (fnErrorHC[0] > 0)? 100.*fnErrorHC[1]/fnErrorHC[0] : -1;
731 fErrorGraphHC->SetPoint(fnEvents, fnEvents, err);
732
733 err = (fnErrorMCM[0] > 0)? 100.*fnErrorMCM[1]/fnErrorMCM[0] : -1;
734 fErrorGraphMCM->SetPoint(fnEvents, fnEvents, err);
735
736 err = (fnErrorADC[0] > 0)? 100.*fnErrorADC[1]/fnErrorADC[0] : -1;
737 fErrorGraphADC->SetPoint(fnEvents, fnEvents, err);
738
739 // number of fired ADC per SM
740 for(Int_t sm=0; sm<kSM+1; sm++)
741 fNumberADC[sm]->SetPoint(fnEvents, fnEvents, fnADCinSM[sm]);
742
743 fnEvents++;
744}
745
746///////////////////////////////////////////////////////////////////////////////////////////////////
747
748void AliTRDqaBlackEvents::Process(const char *filename)
749{
750 //
751 // Process something
752 //
753
754 char fn[256];
8d3a4626 755 strncpy(fn,filename,256);
48ccaa07 756
757 //AliInfo(Form("FILENAME = %s (%s)\n", filename, fn));
758
759 Int_t map[kDET];
760
761 TH1D *hist = new TH1D("fitSignal", "", 50, -0.5, 49.5);
762 TF1 *fit = new TF1("fit", "gaus(0)", 0, 20);
763 fit->SetParameters(1e3, 10, 1);
764
765 for(Int_t det=0; det<kDET; det++) {
766
767 //AliInfo(Form("processing chamber %d\n", det));
768
769 map[det] = 0;
770 //if (fData[det]->GetSum() < 10) continue;
771 //if (fDataDirect[det][10][10][10] < 20) continue;
772 //map[det] = 1;
773
774
775 // rewrite signal-direct
776 for(Int_t ch=0; ch<kCH; ch++) {
777 fSignal[det]->Fill(ch, fSignalDirect[det][ch]);
778 }
779
780 // read reference distributions
781 ReadRefHists(det);
782
783 //for(Int_t row=0; row<fData[det]->GetXaxis()->GetNbins(); row++) {
784 //for(Int_t pad=0; pad<fData[det]->GetYaxis()->GetNbins(); pad++) {
785
786 for(Int_t row=0; row<kROW; row++) {
787 for(Int_t pad=0; pad<kPAD; pad++) {
788
789 // project the histogramm
790 hist->Reset();
791 //for(Int_t bb=0; bb<50; bb++) {
792 for(Int_t bb=0; bb<kCH; bb++) {
793 //Int_t dataBin = fData[det]->FindBin(row, pad, bb);
794 //Double_t v = fData[det]->GetBinContent(dataBin);
795 hist->SetBinContent(bb+1, fDataDirect[det][row][pad][bb]);
796 }
797
798 Int_t bin = fChPed[det]->FindBin(row, pad);
799
800 if (hist->GetSum() > 1) {
801
802 map[det] = 1;
803 Double_t ped = 0, noise = 0;
804
805 if (fFitType == 0) {
806 fit->SetParameters(1e3, 10, 1);
807 hist->Fit(fit, "q0", "goff", 0, 20);
808 TF1 *f = hist->GetFunction("fit");
809 ped = TMath::Abs(f->GetParameter(1));
810 noise = TMath::Abs(f->GetParameter(2));
811 fSmNoiseFit[det/30]->Fill(noise);
812 } else {
813 ped = hist->GetMean();
814 noise = hist->GetRMS();
815 fSmNoiseRms[det/30]->Fill(noise);
816 //if (pad == 0)
817 // AliInfo(Form("data %f %f %f\n", hist->GetSum(), ped, noise));
818 }
819
820 fChPed[det]->SetBinContent(bin, ped);
821 fChNoise[det]->SetBinContent(bin, noise);
822 fNoiseTotal->Fill(noise);
823
824 // subtract reference values
825 Double_t refped = 0;
826 Double_t refnoise = 0;
827
828 if (fRefHistPed) refped = fRefHistPed->GetBinContent(bin);
829 if (fRefHistNoise) refnoise = fRefHistPed->GetBinContent(bin);
830
831 fChPedRes[det]->SetBinContent(bin, ped-refped);
832 fChNoiseRes[det]->SetBinContent(bin, noise-refnoise);
833
834 fPed[det]->Fill(ped);
835 fNoise[det]->Fill(noise);
836
837 // fill SM-HC plot
838 Int_t sm = det / 30;
839 Int_t hc = (pad < kPAD/2) ? 2* (det % 30) : 2* (det % 30) + 1;
840 if (ped > 9. && ped < 11) fSMHCped->Fill(sm, hc, 1./1152.); // number of pads in HC
841
842 } else {
843
844 // not enought data found
845 fChPed[det]->SetBinContent(bin, 0);
846 fChNoise[det]->SetBinContent(bin, 0);
847 fChPedRes[det]->SetBinContent(bin, 0);
848 fChNoiseRes[det]->SetBinContent(bin, 0);
849 }
850
851 //delete hist;
852 }
853 }
854 }
855
856
857 //AliInfo(Form("Number of events = %d\n", fnEvents));
858
859 // normalize number of entries histos
860 Int_t max = 0;
861 for(Int_t i=0; i<kDET; i++) {
862 if (!map[i]) continue;
863 for(Int_t j=0; j<fNPoint[i]->GetXaxis()->GetNbins(); j++) {
864 for(Int_t k=0; k<fNPoint[i]->GetYaxis()->GetNbins(); k++) {
865 Int_t dataBin = fNPoint[i]->FindBin(j, k);
866 Double_t v = fNPoint[i]->GetBinContent(dataBin);
867 if (v > max) max = (Int_t)v;
868 }
869 }
870 }
871
872 char entriesDistName[100];
873
874 for(Int_t i=0; i<kDET; i++) {
875
876 if (!map[i]) continue;
877
878 sprintf(entriesDistName, "entriesDist_%d", i);
879 fNPointDist[i] = new TH1D(entriesDistName, ";number of events", max+2, -0.5, max+1.5);
880
881 for(Int_t j=0; j<fNPoint[i]->GetXaxis()->GetNbins(); j++) {
882 for(Int_t k=0; k<fNPoint[i]->GetYaxis()->GetNbins(); k++) {
883 Int_t dataBin = fNPoint[i]->FindBin(j, k);
884 Double_t v = fNPoint[i]->GetBinContent(dataBin);
885 //if (v > fnEvents) AliInfo(Form("N = %d V = %lf\n", fnEvents, v));
886 fNPointDist[i]->Fill(v);
887 }
888 }
889
890 fNPoint[i]->Scale(1./fnEvents);
891 }
892
893
894 for(Int_t i=0; i<kDET; i++) {
895 fnEntriesRM[i]->SetMaximum(fnEvents * 1.5);
896 }
897
898 // save histograms
899
900 //AliInfo(Form("FILENAME 2 = %s (%d)\n", fn, fn));
901 TFile *file = new TFile(fn, "recreate");
902 for(Int_t det = 0; det < kDET; det++) {
903 if (!map[det]) continue;
904 fChPed[det]->Write();
905 fChNoise[det]->Write();
906 fNPoint[det]->Write();
907 fNPointDist[det]->Write();
908 fPed[det]->Write();
909 fNoise[det]->Write();
910 fSignal[det]->Write();
911 fnEntriesRM[det]->Write();
912 fChPP[det]->Write();
913
914 fChPedRes[det]->Write();
915 fChNoiseRes[det]->Write();
916
917 // save error hists
918 fErrorLocMCM[det]->SetMinimum(0);
919 fErrorLocMCM[det]->SetMaximum(fnEvents);
920 fErrorLocMCM[det]->Write();
921
922 fErrorLocADC[det]->SetMinimum(0);
923 fErrorLocADC[det]->SetMaximum(fnEvents);
924 fErrorLocADC[det]->Write();
925 }
926
927 for(Int_t sm=0; sm<kSM; sm++) {
928 fSmNoiseRms[sm]->Write();
929 fSmNoiseFit[sm]->Write();
930 fSmPP[sm]->Write();
931 }
932
933
934
935 Int_t nMcm = 0;
936 for(Int_t i=0; i < kDET * kROB * kMCM; i++) {
937 if (fFullSignal[i] && fFullCounter[i] > fCount) {
938 fFullSignal[i]->Write();
939 nMcm++;
940 }
941 }
942
943 AliInfo(Form("Number of saved MCMs = %d\n", nMcm));
944
945 fMcmTracks->Write();
946 AliInfo(Form("Number of tracks = %d\n", fMcmTracks->GetEntries()));
947
948 // permanently problematic MCMs
949 for(Int_t det=0; det<kDET; det++) {
950 for(Int_t mcm=0; mcm<kROB*kMCM; mcm++) {
951
952 Int_t mRob = mcm / kMCM;
953 Int_t mMcm = mcm % kMCM;
954 Int_t bin = fMapMCM->FindBin(det, mcm);
955 Double_t frac = 1. * fMapMCM->GetBinContent(bin) / fnEvents;
956 fFracMCM->Fill(frac);
957
958 if (frac > 0.7) {
959 AliInfo(Form("{%d, %d, %d, %f}, \n", det, mRob, mMcm, frac));
960 }
961 }
962 }
963
964
965 fOccupancy->Write();
966 fDetRob->Write();
967 fTBEvent->Write();
968
969 // error hists
970 fErrorHC->Write();
971 fErrorMCM->Write();
972 fErrorADC->Write();
973
974 fErrorSMHC->Write();
975 fErrorSMMCM->Write();
976 fErrorSMADC->Write();
977
978 // write graphs
979 fErrorGraphHC->Write("trendErrorHC");
980 fErrorGraphMCM->Write("trendErrorMCM");
981 fErrorGraphADC->Write("trendErrorADC");
982
983 fGraphMCM->Write("trendMCM");
984
985 for(Int_t i=0; i<3; i++) {
986 fGraphPP[i]->Write(Form("fracPP_%d", i));
987 }
988
989 //fZSsize->Scale(1./fnEvents);
990 //fZSsize->Write();
991
992 fMapMCM->SetMaximum(fnEvents);
993 fMapMCM->Write();
994 fFracMCM->Write();
995
996 fSMHCped->Write();
997
998 for(Int_t i=0; i<3; i++ ) {
999 fSMLink[i]->Write();
1000 fGrLink[i]->Write();
1001 }
1002
1003 for(Int_t sm=0; sm<kSM; sm++)
1004 fNumberADC[sm]->Write(Form("nADCinSM%d",sm));
1005
1006 fNumberADC[kSM]->Write("nADCinEvent");
1007
1008 fNoiseTotal->Write();
1009 fPP->Write();
1010
1011 for(Int_t i=0; i<1000; i++) {
1012 if (fEvNoDist[i]->GetSum() > 0) fEvNoDist[i]->Write();
1013 }
1014
1015
1016 file->Close();
1017 delete file;
1018}
1019
1020///////////////////////////////////////////////////////////////////////////////////////////////////
1021
c87495af 1022Int_t AliTRDqaBlackEvents::CheckMCM(Int_t index) const {
1023 //
1024 // Checks a single MCM
1025 //
1026
48ccaa07 1027 return 1;
1028
1029 static Int_t data[21][3] = {
1030 {1, 0, 1},
1031 {242, 0, 0},
1032 {242, 0, 1},
1033 {242, 0, 2},
1034 {242, 0, 4},
1035 {242, 0, 5},
1036 {242, 0, 6},
1037 {242, 0, 8},
1038 {242, 0, 12},
1039 {251, 7, 7},
1040 {254, 3, 11},
1041 {259, 3, 14},
1042 {260, 1, 9},
1043 {260, 3, 15},
1044 {273, 1, 7},
1045 {273, 1, 15},
1046 {276, 5, 11},
1047 {280, 6, 2},
1048 {299, 6, 4},
1049 {511, 2, 9},
1050 {517, 7, 15}
1051 };
1052
1053 for(Int_t i=0; i<21; i++) {
1054 Int_t wIndex = data[i][0] * kROB*kMCM + data[i][1] * kMCM + data[i][2];
1055 if (index == wIndex) return 0;
1056 }
1057
1058 return 1;
1059}
1060
1061///////////////////////////////////////////////////////////////////////////////////////////////////
1062
1063
1064
1065
1066
1067
1068
1069void AliTRDqaBlackEvents::DrawChamber(const char *filename, Int_t det, Int_t w, Int_t h)
1070{
1071 //
1072 // Draw raport for one chamber:
1073 // pedestal map, noise map, distribution of pedestal and noise
1074 //
1075 // input:
1076 // name of the file with histograms (created with Process())
1077 // detector Id (0 - 539)
1078 //
1079
1080 // setup global style
1081 gStyle->SetPalette(1);
1082 gStyle->SetOptStat(0);
1083 gStyle->SetPadTopMargin(0.02);
1084 gStyle->SetPadBottomMargin(0.05);
1085
1086 TFile *file = new TFile(filename, "READ");
1087
1088 TCanvas *c = new TCanvas("blackEvents",Form("blackEvents %d",det), w, h);
1089 c->SetVertical(kFALSE);
1090 c->Divide(3,1, 0.01, 0.01);
1091 c->cd(3);
1092
1093 TPad *mPad = (TPad*) gPad;
1094 mPad->Divide(1,2,0.01,0.01);
1095
1096 c->cd(1);
1097 TH2D *h2 = (TH2D*)file->Get(Form("ped_%d",det));
1098 h2->SetMinimum(5);
1099 h2->SetMaximum(15);
1100 h2->SetTitle(";Z direction;#phi direction");
1101 h2->Draw("colz");
1102
1103 c->cd(2);
1104 h2 = (TH2D*)file->Get(Form("noise_%d",det));
1105 h2->SetMinimum(fMinNoise);
1106 h2->SetMaximum(fMaxNoise);
1107 h2->SetTitle(";Z direction;#phi direction");
1108 h2->Draw("colz");
1109
1110 mPad->cd(1);
1111 //gPad->SetLogy();
1112 TH1D *h1 = (TH1D*)file->Get(Form("pedDist_%d", det));
1113 h1->Draw();
1114
1115 mPad->cd(2);
1116 gPad->SetLogy();
1117 h1 = (TH1D*)file->Get(Form("noiseDist_%d", det));
1118 h1->Draw();
1119
1120 h1->Fit("gaus");
1121 TF1 *f = h1->GetFunction("gaus");
1122 const char *tt = Form("#mu = %.2f #sigma = %0.2f ", f->GetParameter(1),f->GetParameter(2));
1123 TLatex *ll = new TLatex(2, 100, tt);
1124 ll->SetTextSize(0.06);
1125 ll->Draw();
1126}
1127
1128///////////////////////////////////////////////////////////////////////////////////////////////////
1129
1130void AliTRDqaBlackEvents::DrawSm(const char *filename, Int_t sm, Int_t w, Int_t h)
1131{
1132 //
1133 // ????????????
1134 //
1135
1136 gStyle->SetPalette(1);
1137 gStyle->SetOptStat(0);
1138
1139 gStyle->SetPadTopMargin(0.02);
1140 //gStyle->SetPadBottomMargin(0.05);
1141 //gStyle->SetPadLeftMargin(0.02);
1142 //gStyle->SetPadRightMargin(0.02);
1143
1144 TFile *file = new TFile(filename, "READ");
1145
1146 TCanvas *c = new TCanvas("blackEventsSM",Form("blackEvents SM %d",sm), w, h);
1147 c->SetVertical(kFALSE);
1148 c->Divide(5, 6, 0.001, 0.01);
1149
1150 for(Int_t i=0; i<30; i++) {
1151
1152 TH2D *h2 = (TH2D*)file->Get(Form("noise_%d",i+30*sm));
1153 if (!h2) continue;
1154 h2->SetMinimum(fMinNoise);
1155 h2->SetMaximum(fMaxNoise);
1156
1157 // to be replaced by the official calculation
1158 Int_t stack = i/6;
1159 Int_t layer = i%6;
1160 Int_t index = (5-layer)*5 + stack + 1;
1161 //AliInfo(Form("%d %d %d %d\n", i, stack, layer, index));
1162 c->cd(index);
1163 gPad->SetBottomMargin(0.02);
1164 gPad->SetTopMargin(0.02);
1165
1166 h2->Draw("col");
1167 }
1168}
1169
1170///////////////////////////////////////////////////////////////////////////////////////////////////
1171
1172Int_t AliTRDqaBlackEvents::FillBits(TH1D *hist, Int_t code, Int_t offset) {
c87495af 1173 //
1174 // Fill bits
1175 //
48ccaa07 1176
1177 Int_t nb = 0;
1178 UInt_t test = 1;
1179 for(Int_t i=0; i<8; i++) {
1180 if (code & test) {
1181 hist->Fill(i+offset);
1182 nb++;
1183 }
1184 test *= 2;
1185 }
1186
1187 return nb;
1188}
1189
1190///////////////////////////////////////////////////////////////////////////////////////////////////