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