]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - TOF/AliTOFQADataMakerRec.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (non-incremental) | history | HEAD
Coverity fixes.
[u/mrichter/AliRoot.git] / TOF / AliTOFQADataMakerRec.cxx
... / ...
CommitLineData
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 * without fee, provided that the above 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///////////////////////////////////////////////////////////////////////
17// //
18// Produces the data needed to calculate the TOF quality assurance. //
19// QA objects are 1 & 2 Dimensional histograms. //
20// author: S.Arcelli //
21// //
22///////////////////////////////////////////////////////////////////////
23
24/*
25Modified by fbellini on 01/11/2011
26- added histograms for LTM monitoring
27- fix for coverity
28
29Modified by fbellini on 17/10/2011
30- fix for memory leak in constructor
31- added methods to read histos ranges from config file in DQM
32- added CTTM maps + relative methods to retrieve CTTM numbering
33- removed obslete comments
34
35Modified by fbellini & rshanoian on 06/07/2011
36- changes for trigger classes implementation
37- fRunNumber added as private member
38- added time vs BCID plot
39
40Modified by fbellini on 18/01/2011
41- reduced histo binning to reduce size
42- added decoding errors plot
43- added channel maps and options for DQM shifters
44- new list of recPoints and ESDs plots
45- removed hTrmChannels035 and hTrmChannels3671
46
47 Modified by bvonhall on 03/11/2010
48 - modified declaration of hTrmChannels035 and hTrmChannels3671 in EndOfDetectorCycle()
49 to prevent memory corruption
50
51 Modified by adecaro on 18/10/2010
52 - fTOFRawStream object set as private member
53
54Modified by fbellini on 13/09/2010
55 - Set TLines as private members
56 - Set image flag for expert histos
57
58Modified by fbellini on 14/06/2010
59 - Updated plots
60 - use LoadRawDataBuffersV2()
61
62 Modified by fbellini on 10/05/2010
63 - Fixed EndOfDetectorCycle() memory corruption bug
64
65 Modified by fbellini on 22/04/2010
66 - Added filter for physics events
67
68 Modified by fbellini on 16/04/2010
69 - Added EnableDqmShifterOpt()
70 - Modified EndOfDetectorCycle() with options for DQM
71 - Updated ESDs QA
72
73 Modified by fbellini on 30/03/2010
74 - Changed raws time histos range to 610ns
75 - Added FilterLTMData() and FilterSpare() methods
76 - Added check on enabled channels for raw data
77 - Updated RecPoints QA
78
79 Modified by fbellini on 02/03/2010
80 - Fixed raw data decoding methods (use AliTOFRawStream::LoadRawDataBuffer())
81 - Added filter for noisy channels and read map from OCDB
82 - Added GetCalibData() method
83 - Added CheckVolumeID() and CheckEquipID() methods
84 - Updated Raw QA
85*/
86#include <iostream>
87#include <fstream>
88
89#include <TClonesArray.h>
90#include <TH1F.h>
91#include <TH2F.h>
92#include <TLine.h>
93#include <TPaveText.h>
94
95#include "AliCDBManager.h"
96#include "AliCDBEntry.h"
97#include "AliESDEvent.h"
98#include "AliESDtrack.h"
99#include "AliQAChecker.h"
100#include "AliRawReader.h"
101#include "AliTOFRawStream.h"
102#include "AliTOFcluster.h"
103#include "AliTOFQADataMakerRec.h"
104#include "AliTOFrawData.h"
105#include "AliTOFGeometry.h"
106#include "AliTOFChannelOnlineStatusArray.h"
107#include "AliTOFDecoderSummaryData.h"
108#include "AliTOFDecoderV2.h"
109
110ClassImp(AliTOFQADataMakerRec)
111
112Int_t AliTOFQADataMakerRec::fgNbinsMultiplicity=200; //number of bins in multiplicity plot
113Int_t AliTOFQADataMakerRec::fgRangeMinMultiplicity=0;//min range in multiplicity plot
114Int_t AliTOFQADataMakerRec::fgRangeMaxMultiplicity=200;//max range in multiplicity plot
115Int_t AliTOFQADataMakerRec::fgNbinsTime=250;//number of bins in time plot
116const Float_t AliTOFQADataMakerRec::fgkNbinsWidthTime=2.44;//width of bins in time plot
117Float_t AliTOFQADataMakerRec::fgRangeMinTime=0.0;//range min in time plot
118Float_t AliTOFQADataMakerRec::fgRangeMaxTime=620.0; //range max in time plot
119Int_t AliTOFQADataMakerRec::fgCutNmaxFiredMacropad=5;//cut on number of max fired macropad
120const Int_t AliTOFQADataMakerRec::fgkFiredMacropadLimit=50;//cut on number of max fired macropad
121
122
123//____________________________________________________________________________
124 AliTOFQADataMakerRec::AliTOFQADataMakerRec() :
125 AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kTOF), "TOF Quality Assurance Data Maker"),
126 fCalibData(0x0),
127 fEnableNoiseFiltering(kFALSE),
128 fEnableDqmShifterOpt(kFALSE),
129 fIsSOC(kFALSE),
130 fLineExpTimeMin(0x0),
131 fLineExpTimeMax(0x0),
132 fLineExpTotMin(0x0),
133 fLineExpTotMax(0x0),
134 fTOFRawStream(AliTOFRawStream()),
135 fDecoderSummary(0),
136 fRunNumber(-1),
137 fCalib(AliTOFcalib())
138{
139 //
140 // ctor
141 //
142 // fLineExpTimeMin = new TLine(200., 0., 200., 0.);
143 // fLineExpTimeMax = new TLine(250., 0., 250., 0.);
144 // fLineExpTotMin = new TLine(5., 0., 5., 0.);
145 // fLineExpTotMax = new TLine(20., 0., 20., 0.);
146 for (Int_t sm=0;sm<17;sm++){
147 fLineSMid[sm] = new TLine( sm+1, 0., sm+1, 91.);
148 }
149
150 for (Int_t sm=0;sm<71;sm++){
151 fLineLTMid[sm] = new TLine( sm+1, 0., sm+1, 23.);
152 }
153
154 for (Int_t sm=0;sm<22;sm++){
155 fLineLTMbitId[sm] = new TLine( 0., sm+1, 72. ,sm+1);
156 }
157
158}
159
160//____________________________________________________________________________
161AliTOFQADataMakerRec::AliTOFQADataMakerRec(const AliTOFQADataMakerRec& qadm) :
162 AliQADataMakerRec(),
163 fCalibData(qadm.fCalibData),
164 fEnableNoiseFiltering(qadm.fEnableNoiseFiltering),
165 fEnableDqmShifterOpt(qadm.fEnableDqmShifterOpt),
166 fIsSOC(qadm.fIsSOC),
167 fLineExpTimeMin(qadm.fLineExpTimeMin),
168 fLineExpTimeMax(qadm.fLineExpTimeMax),
169 fLineExpTotMin(qadm.fLineExpTotMin),
170 fLineExpTotMax(qadm.fLineExpTotMax),
171 fTOFRawStream(qadm.fTOFRawStream),
172 fDecoderSummary(qadm.fDecoderSummary),
173 fRunNumber(qadm.fRunNumber),
174 fCalib(qadm.fCalib)
175{
176 //
177 //copy ctor
178 //
179 SetName((const char*)qadm.GetName()) ;
180 SetTitle((const char*)qadm.GetTitle());
181
182 for (Int_t sm=0;sm<17;sm++){
183 fLineSMid[sm]=qadm.fLineSMid[sm];
184 }
185
186 for (Int_t sm=0;sm<71;sm++){
187 fLineLTMid[sm] = qadm.fLineLTMid[sm];
188 }
189
190 for (Int_t sm=0;sm<22;sm++){
191 fLineLTMbitId[sm] = qadm.fLineLTMbitId[sm];
192 }
193}
194
195//__________________________________________________________________
196AliTOFQADataMakerRec& AliTOFQADataMakerRec::operator = (const AliTOFQADataMakerRec& qadm )
197{
198 //
199 // assignment operator.
200 //
201 this->~AliTOFQADataMakerRec();
202 new(this) AliTOFQADataMakerRec(qadm);
203 return *this;
204}
205
206//----------------------------------------------------------------------------
207AliTOFQADataMakerRec::~AliTOFQADataMakerRec()
208{
209 //destructor
210 fTOFRawStream.Clear();
211 fCalib.Clear();
212 if (fLineExpTimeMin)
213 delete fLineExpTimeMin;
214 if (fLineExpTimeMax)
215 delete fLineExpTimeMax;
216 if (fLineExpTotMin)
217 delete fLineExpTotMin;
218 if (fLineExpTotMax)
219 delete fLineExpTotMax;
220 for (Int_t sm=0;sm<17;sm++){
221 if (fLineSMid[sm])
222 delete fLineSMid[sm];
223 }
224 for (Int_t sm=0;sm<71;sm++){
225 if (fLineLTMid[sm])
226 delete fLineLTMid[sm];
227 }
228for (Int_t sm=0;sm<22;sm++){
229 if (fLineLTMbitId[sm])
230 delete fLineLTMbitId[sm];
231 }
232}
233//----------------------------------------------------------------------------
234AliTOFChannelOnlineStatusArray* AliTOFQADataMakerRec::GetCalibData()
235{
236 //
237 // Retrive TOF calib objects from OCDB
238 //
239 AliCDBManager *man = AliCDBManager::Instance();
240 AliCDBEntry *cdbe=0;
241
242 if (fRun<=0) fRunNumber=145288; //reference run from LHC11a
243 else fRunNumber=fRun;
244
245 if (man->GetRun()!=fRunNumber){
246 fRunNumber=man->GetRun();
247 AliWarning(Form("Run number mismatch found: setting it to value from current AliCDBManager instance = %i", fRunNumber));
248 }
249 cdbe = man->Get("TOF/Calib/Status",fRunNumber);
250
251 if(!cdbe){
252 // for DQM online
253 AliWarning("Load of calibration data from default (alien://) storage failed!");
254 printf("Calibration data will be loaded from local storage - ok if on DQM station!");
255 man->SetDefaultStorage("local:///local/cdb/");
256 cdbe = man->Get("TOF/Calib/Status",fRun);
257
258 if(!cdbe){
259 AliWarning("Load of calibration data from local DQM machine storage failed!");
260 AliWarning("Calibration data will be loaded from local ($ALICE_ROOT) storage ");
261 man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
262 cdbe = man->Get("TOF/Calib/Status",fRunNumber);
263 }
264 }
265 // Retrieval of data in directory TOF/Calib/Data:
266 AliTOFChannelOnlineStatusArray * array = 0;
267 if (cdbe) {
268 printf("======= OCDB object for TOF retrieved from run %i in %s\n",fRunNumber,cdbe->GetName());
269 array = (AliTOFChannelOnlineStatusArray *)cdbe->GetObject();
270 }
271 if (!array) AliFatal("No calibration data from calibration database !");
272
273 fCalib.Init(fRunNumber);
274 return array;
275}
276
277//____________________________________________________________________________
278void AliTOFQADataMakerRec::InitRaws()
279{
280 //
281 // create Raws histograms in Raws subdir
282 //
283 ReadHistogramRangeFromFile(gSystem->Getenv("TOFDQMHISTO_CONFIGFILE"));
284
285 const Bool_t expert = kTRUE ;
286 const Bool_t saveCorr = kTRUE ;
287 const Bool_t image = kTRUE ;
288
289 TH1I * h0 = new TH1I("hTOFRaws","TOF raw hit multiplicity; TOF raw hits number; Events ",fgNbinsMultiplicity, fgRangeMinMultiplicity, fgRangeMaxMultiplicity);
290 TH1I * h1 = new TH1I("hTOFRawsIA","TOF raw hit multiplicity - I/A side; TOF raw hits number;Events ",fgNbinsMultiplicity, fgRangeMinMultiplicity, fgRangeMaxMultiplicity);
291 TH1I * h2 = new TH1I("hTOFRawsOA","TOF raw hit multiplicity - O/A side; TOF raw hits number;Events ",fgNbinsMultiplicity, fgRangeMinMultiplicity, fgRangeMaxMultiplicity);
292 TH1I * h3 = new TH1I("hTOFRawsIC","TOF raw hit multiplicity - I/C side; TOF raw hits number;Events ",fgNbinsMultiplicity, fgRangeMinMultiplicity, fgRangeMaxMultiplicity);
293 TH1I * h4 = new TH1I("hTOFRawsOC","TOF raw hit multiplicity - O/C side; TOF raw hits number;Events ",fgNbinsMultiplicity, fgRangeMinMultiplicity, fgRangeMaxMultiplicity);
294
295 TH1F * h5 = new TH1F("hTOFRawsTime", "TOF Raws - Hit time (ns);Measured Hit time [ns];Hits", fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
296 TH1F * h6 = new TH1F("hTOFRawsTimeIA", "TOF Raws - Hit time (ns) - I/A side;Measured Hit time [ns];Hits", fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
297 TH1F * h7 = new TH1F("hTOFRawsTimeOA", "TOF Raws - Hit time (ns) - O/A side;Measured Hit time [ns];Hits", fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
298 TH1F * h8 = new TH1F("hTOFRawsTimeIC", "TOF Raws - Hit time (ns) - I/C side;Measured Hit time [ns];Hits", fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
299 TH1F * h9 = new TH1F("hTOFRawsTimeOC", "TOF Raws - Hit time (ns) - O/C side;Measured Hit time [ns];Hits", fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
300
301 TH1F * h10 = new TH1F("hTOFRawsToT", "TOF Raws - Hit ToT (ns);Measured Hit ToT (ns);Hits", 100, 0., 48.8);
302
303 TH1F * h11 = new TH1F("hTOFRawsToTIA", "TOF Raws - Hit ToT (ns) - I/A side;Measured Hit ToT (ns);Hits", 100, 0., 48.8);
304 TH1F * h12 = new TH1F("hTOFRawsToTOA", "TOF Raws - Hit ToT (ns) - O/A side;Measured Hit ToT (ns);Hits", 100, 0., 48.8);
305 TH1F * h13 = new TH1F("hTOFRawsToTIC", "TOF Raws - Hit ToT (ns) - I/C side;Measured Hit ToT (ns);Hits", 100, 0., 48.8);
306 TH1F * h14 = new TH1F("hTOFRawsToTOC", "TOF Raws - Hit ToT (ns) - O/C side;Measured Hit ToT (ns);Hits", 100, 0., 48.8);
307
308 TH1F * h15 = new TH1F("hTOFRawsLTMHits", "LTMs OR signals; Crate; Counts", 72, 0., 72.);
309 TH2F * h16 = new TH2F("hTOFrefMap", "TOF enabled channel reference map;sector;strip", 72, 0., 18., 91, 0., 91.);
310 TH2F * h17 = new TH2F("hTOFRawHitMap","TOF raw hit map;sector;strip", 72, 0., 18., 91, 0., 91.);
311 TH2I * h18 = new TH2I("hTOFDecodingErrors","Decoding error monitoring; DDL; Error ", 72, 0, 72, 13,1,14);
312
313 h18->GetYaxis()->SetBinLabel(1,"DRM ");
314 h18->GetYaxis()->SetBinLabel(2,"LTM ");
315 h18->GetYaxis()->SetBinLabel(3,"TRM 3 ");
316 h18->GetYaxis()->SetBinLabel(4,"TRM 4");
317 h18->GetYaxis()->SetBinLabel(5,"TRM 5");
318 h18->GetYaxis()->SetBinLabel(6,"TRM 6");
319 h18->GetYaxis()->SetBinLabel(7,"TRM 7");
320 h18->GetYaxis()->SetBinLabel(8,"TRM 8");
321 h18->GetYaxis()->SetBinLabel(9,"TRM 9");
322 h18->GetYaxis()->SetBinLabel(10,"TRM 10");
323 h18->GetYaxis()->SetBinLabel(11,"TRM 11");
324 h18->GetYaxis()->SetBinLabel(12,"TRM 12");
325 h18->GetYaxis()->SetBinLabel(13,"recovered");
326
327 TH1F * h19 = new TH1F("hTOFOrphansTime", "TOF Raws - Orphans time (ns);Measured Hit time [ns];Hits",fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
328 TH2F * h20 = new TH2F("hTOFRawTimeVsTRM035", "TOF raws - Hit time vs TRM - crates 0 to 35; TRM index = DDL*10+TRM(0-9);TOF raw time [ns]", 361, 0., 361.,fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
329 TH2F * h21 = new TH2F("hTOFRawTimeVsTRM3671", "TOF raws - Hit time vs TRM - crates 36 to 72; TRM index = DDL**10+TRM(0-9);TOF raw time [ns]", 361, 360., 721.,fgNbinsTime, fgRangeMinTime,fgRangeMaxTime);
330 TH2F * h22 = new TH2F("hTOFTimeVsStrip","TOF raw hit time vs. MRPC (along z axis); MRPC index along z axis; Raws TOF time (ns) ", 91,0.,91,fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
331 TH2F * h23 = new TH2F("hTOFtimeVsBCID","TOF time vs BCID; BCID; time (ns) ", 3564, 0., 3564.,fgNbinsTime,fgRangeMinTime,fgRangeMaxTime);
332 TH2F * h24 = new TH2F("hTOFchannelEfficiencyMap","TOF channels (HWok && efficient && !noisy && !problematic);sector;strip", 72, 0., 18., 91, 0., 91.);
333 TH2F * h25 = new TH2F("hTOFhitsCTTM","Map of hit pads according to CTTM numbering;LTM index;bit index", 72, 0., 72., 23, 0., 23.);
334 TH2F * h26 = new TH2F("hTOFmacropadCTTM","Map of hit macropads according to CTTM numbering;LTM index; bit index", 72, 0., 72., 23, 0., 23.);
335 TH2F * h27 = new TH2F("hTOFmacropadDeltaPhiTime","#Deltat vs #Delta#Phi of hit macropads;#Delta#Phi (degrees);#DeltaBX", 18, 0., 180., 20, 0., 20.0);
336 TH2I *h28 = new TH2I("hBXVsCttmBit","BX ID in TOF matching window vs trg channel; trg channel; BX", 1728, 0, 1728, 24, 0, 24);
337 TH2F *h29 = new TH2F("hTimeVsCttmBit","TOF raw time vs trg channel; trg channel; raw time (ns)", 1728, 0., 1728., fgNbinsTime, fgRangeMinTime, fgRangeMaxTime);
338
339 h25->GetYaxis()->SetTickLength(-0.02);
340 h26->GetYaxis()->SetTickLength(-0.02);
341 h25->GetYaxis()->SetNdivisions(210);
342 h26->GetYaxis()->SetNdivisions(210);
343 h25->GetXaxis()->SetTickLength(-0.02);
344 h26->GetXaxis()->SetTickLength(-0.02);
345 h25->GetXaxis()->SetLabelOffset(0.015);
346 h26->GetXaxis()->SetLabelOffset(0.015);
347 h25->GetXaxis()->SetNdivisions(515);
348 h26->GetXaxis()->SetNdivisions(515);
349
350 h0->Sumw2() ;
351 h1->Sumw2() ;
352 h2->Sumw2() ;
353 h3->Sumw2() ;
354 h4->Sumw2() ;
355 // h5->Sumw2() ;
356 h6->Sumw2() ;
357 h7->Sumw2() ;
358 h8->Sumw2() ;
359 h9->Sumw2() ;
360 // h10->Sumw2() ;
361 h11->Sumw2() ;
362 h12->Sumw2() ;
363 h13->Sumw2() ;
364 h14->Sumw2() ;
365 h15->Sumw2() ;
366 h16->Sumw2() ;
367 h17->Sumw2() ;
368 h18->Sumw2() ;
369 h19->Sumw2() ;
370 h20->Sumw2() ;
371 h21->Sumw2() ;
372 h22->Sumw2() ;
373 h23->Sumw2() ;
374 h24->Sumw2() ;
375 h25->Sumw2() ;
376 h26->Sumw2() ;
377 h27->Sumw2() ;
378 h28->Sumw2() ;
379 h29->Sumw2() ;
380
381 //add lines for DQM shifter
382 fLineExpTimeMin = new TLine(200., 0., 200., 0.);
383 fLineExpTimeMax = new TLine(250., 0., 250., 0.);
384 fLineExpTotMin = new TLine(5., 0., 5., 0.);
385 fLineExpTotMax = new TLine(20., 0., 20., 0.);
386
387 fLineExpTimeMin->SetLineColor(kGreen);
388 fLineExpTimeMin->SetLineWidth(2);
389
390 fLineExpTimeMax->SetLineColor(kGreen);
391 fLineExpTimeMax->SetLineWidth(2);
392
393 fLineExpTotMin->SetLineColor(kGreen);
394 fLineExpTotMin->SetLineWidth(2);
395
396 fLineExpTotMax->SetLineColor(kGreen);
397 fLineExpTotMax->SetLineWidth(2);
398
399 for (Int_t sm=0;sm<17;sm++){
400 fLineSMid[sm]->SetLineColor(kMagenta);
401 fLineSMid[sm]->SetLineWidth(2);
402 }
403
404 h5->GetListOfFunctions()->Add(fLineExpTimeMin);
405 h5->GetListOfFunctions()->Add(fLineExpTimeMax);
406 h10->GetListOfFunctions()->Add(fLineExpTotMin);
407 h10->GetListOfFunctions()->Add(fLineExpTotMax);
408
409 for (Int_t sm=0;sm<17;sm++){
410 h16->GetListOfFunctions()->Add(fLineSMid[sm]);
411 h17->GetListOfFunctions()->Add(fLineSMid[sm]);
412 }
413
414 for (Int_t sm=0;sm<71;sm++){
415 fLineLTMid[sm]->SetLineColor(kBlack);
416 fLineLTMid[sm]->SetLineWidth(1);
417 h26->GetListOfFunctions()->Add(fLineLTMid[sm]);
418 h25->GetListOfFunctions()->Add(fLineLTMid[sm]);
419 }
420 for (Int_t sm=0;sm<22;sm++){
421 fLineLTMbitId[sm]->SetLineColor(kBlack);
422 fLineLTMbitId[sm]->SetLineWidth(1);
423 h26->GetListOfFunctions()->Add(fLineLTMbitId[sm]);
424 h25->GetListOfFunctions()->Add(fLineLTMbitId[sm]);
425 }
426
427 TPaveText *phosHoleBox=new TPaveText(13,38,16,53,"b");
428 phosHoleBox->SetFillStyle(0);
429 phosHoleBox->SetFillColor(kWhite);
430 phosHoleBox->SetLineColor(kMagenta);
431 phosHoleBox->SetLineWidth(2);
432 phosHoleBox->AddText("PHOS");
433 h16->GetListOfFunctions()->Add(phosHoleBox);
434 h17->GetListOfFunctions()->Add(phosHoleBox);
435
436 // h0->SetDrawOption("logy");
437 // h5->SetDrawOption("logy");
438 // h10->SetDrawOption("logy");
439
440 Add2RawsList(h0, 0, !expert, image, !saveCorr) ;
441 Add2RawsList(h1, 1, expert, !image, !saveCorr) ;
442 Add2RawsList(h2, 2, expert, !image, !saveCorr) ;
443 Add2RawsList(h3, 3, expert, !image, !saveCorr) ;
444 Add2RawsList(h4, 4, expert, !image, !saveCorr) ;
445 Add2RawsList(h5, 5, !expert, image, !saveCorr) ;
446 Add2RawsList(h6, 6, expert, !image, !saveCorr) ;
447 Add2RawsList(h7, 7, expert, !image, !saveCorr) ;
448 Add2RawsList(h8, 8, expert, !image, !saveCorr) ;
449 Add2RawsList(h9, 9, expert, !image, !saveCorr) ;
450 Add2RawsList(h10, 10, !expert, image, !saveCorr) ;
451 Add2RawsList(h11, 11, expert, !image, !saveCorr) ;
452 Add2RawsList(h12, 12, expert, !image, !saveCorr) ;
453 Add2RawsList(h13, 13, expert, !image, !saveCorr) ;
454 Add2RawsList(h14, 14, expert, !image, !saveCorr) ;
455 Add2RawsList(h15, 15, expert, !image, !saveCorr) ;
456 Add2RawsList(h16, 16, !expert, image, !saveCorr) ;
457 Add2RawsList(h17, 17, !expert, image, !saveCorr) ;
458 Add2RawsList(h18, 18, expert, !image, !saveCorr) ;
459 Add2RawsList(h19, 19, expert, !image, !saveCorr) ;
460 Add2RawsList(h20, 20, expert, !image, !saveCorr) ;
461 Add2RawsList(h21, 21, expert, !image, !saveCorr) ;
462 Add2RawsList(h22, 22, !expert, image, !saveCorr) ;
463 Add2RawsList(h23, 23, !expert, !image, !saveCorr) ;
464 Add2RawsList(h24, 24, !expert, !image, !saveCorr) ;
465 Add2RawsList(h25, 25, !expert, !image, !saveCorr) ;
466 Add2RawsList(h26, 26, !expert, image, !saveCorr) ;
467 Add2RawsList(h27, 27, !expert, image, !saveCorr) ;
468 Add2RawsList(h28, 28, !expert, !image, !saveCorr) ;
469 Add2RawsList(h29, 29, !expert, !image, !saveCorr) ;
470
471//
472 ClonePerTrigClass(AliQAv1::kRAWS); // this should be the last line
473}
474
475//____________________________________________________________________________
476void AliTOFQADataMakerRec::InitRecPoints()
477{
478 //
479 // create RecPoints histograms in RecPoints subdir
480 //
481
482 const Bool_t expert = kTRUE ;
483 const Bool_t image = kTRUE ;
484
485 TH1I * h0 = new TH1I("hTOFRecPoints", "TOF RecPoints multiplicity ; TOF RecPoints number;Events",200, 0, 200) ;
486
487 TH1F * h1 = new TH1F("hTOFRecPointsTimeIA", "RecPoints Time Spectrum in TOF (ns)- I/A side; Calibrated TOF time [ns];Events", 250, 0., 610.) ;
488 TH1F * h2 = new TH1F("hTOFRecPointsTimeOA", "RecPoints Time Spectrum in TOF (ns)- O/A side; Calibrated TOF time [ns];Events", 250, 0., 610.) ;
489 TH1F * h3 = new TH1F("hTOFRecPointsTimeIC", "RecPoints Time Spectrum in TOF (ns)- I/C side; Calibrated TOF time [ns];Events", 250, 0., 610.) ;
490 TH1F * h4 = new TH1F("hTOFRecPointsTimeOC", "RecPoints Time Spectrum in TOF (ns)- O/C side; Calibrated TOF time [ns];Events", 250, 0., 610.) ;
491
492 TH1F * h5 = new TH1F("hTOFRecPointsRawTimeIA", "RecPoints raw Time Spectrum in TOF (ns)-I/A side; Measured TOF time [ns];Hits", 250, 0., 610.) ;
493 TH1F * h6 = new TH1F("hTOFRecPointsRawTimeOA", "RecPoints raw Time Spectrum in TOF (ns)-O/A side; Measured TOF time [ns];Hits", 250, 0., 610.) ;
494 TH1F * h7 = new TH1F("hTOFRecPointsRawTimeIC", "RecPoints raw Time Spectrum in TOF (ns)-I/C side; Measured TOF time [ns];Hits", 250, 0., 610.) ;
495 TH1F * h8 = new TH1F("hTOFRecPointsRawTimeOC", "RecPoints raw Time Spectrum in TOF (ns)-O/C side; Measured TOF time [ns];Hits", 250, 0., 610.) ;
496
497 TH1F * h9 = new TH1F("hTOFRecPointsToTIA", "RecPoints ToT Spectrum in TOF (ns)-I/A side; Measured TOT [ns];Hits", 100, 0., 48.8 ) ;
498 TH1F * h10 = new TH1F("hTOFRecPointsToTOA", "RecPoints ToT Spectrum in TOF (ns)-O/A side; Measured TOT [ns];Hits", 100, 0., 48.8 ) ;
499 TH1F * h11 = new TH1F("hTOFRecPointsToTIC", "RecPoints ToT Spectrum in TOF (ns)-I/C side; Measured TOT [ns];Hits", 100, 0., 48.8 ) ;
500 TH1F * h12 = new TH1F("hTOFRecPointsToTOC", "RecPoints ToT Spectrum in TOF (ns)-O/C side; Measured TOT [ns];Hits", 100, 0., 48.8 ) ;
501
502 TH2F * h13 = new TH2F("hTOFRecPointsClusMap","RecPoints map; sector ;strip", 72, 0., 18., 91, 0., 91.) ;
503 TH2F * h14 = new TH2F("hTOFRecPointsTimeVsStrip","RecPoints TOF time vs. strip (theta); Strip index; RecPoints TOF time (ns) ",91, 0., 91., 250, 0., 610.) ;
504 TH2F * h15 = new TH2F("hTOFRecPointsTimeVsTRM035","TOF RecPoints time vs TRM - crates 0 to 35; TRM index = DDL*10+TRM(0-9);TOF time [ns]", 361, 0., 361., 250, 0., 610.0) ;
505 TH2F * h16 = new TH2F("hTOFRecPointsTimeVsTRM3671","TOF RecPoints time vs TRM - crates 36 to 72; TRM index = DDL**10+TRM(0-9);TOF time [ns]", 361, 360., 721., 250, 0., 610.0) ;
506
507 h0->Sumw2() ;
508 h1->Sumw2() ;
509 h2->Sumw2() ;
510 h3->Sumw2() ;
511 h4->Sumw2() ;
512 h5->Sumw2() ;
513 h6->Sumw2() ;
514 h7->Sumw2() ;
515 h8->Sumw2() ;
516 h9->Sumw2() ;
517 h10->Sumw2() ;
518 h11->Sumw2() ;
519 h12->Sumw2() ;
520 h13->Sumw2() ;
521 h14->Sumw2() ;
522 h15->Sumw2() ;
523 h16->Sumw2() ;
524
525 Add2RecPointsList(h0, 0, !expert, image) ;
526 Add2RecPointsList(h1, 1, !expert, image) ;
527 Add2RecPointsList(h2, 2, !expert, image) ;
528 Add2RecPointsList(h3, 3, !expert, image) ;
529 Add2RecPointsList(h4, 4, !expert, image) ;
530 Add2RecPointsList(h5, 5, expert, !image) ;
531 Add2RecPointsList(h6, 6, expert, !image) ;
532 Add2RecPointsList(h7, 7, expert, !image) ;
533 Add2RecPointsList(h8, 8, expert, !image) ;
534 Add2RecPointsList(h9, 9, !expert, !image) ;
535 Add2RecPointsList(h10, 10, !expert, !image) ;
536 Add2RecPointsList(h11, 11, !expert, !image) ;
537 Add2RecPointsList(h12, 12, !expert, !image) ;
538 Add2RecPointsList(h13, 13, expert, !image) ;
539 Add2RecPointsList(h14, 14, expert, image) ;
540 Add2RecPointsList(h15, 15, expert, !image) ;
541 Add2RecPointsList(h16, 16, expert, !image) ;
542 //
543 ClonePerTrigClass(AliQAv1::kRECPOINTS); // this should be the last line
544}
545
546//____________________________________________________________________________
547void AliTOFQADataMakerRec::InitESDs()
548{
549 //
550 //create ESDs histograms in ESDs subdir
551 //
552
553 const Bool_t expert = kTRUE ;
554 const Bool_t image = kTRUE ;
555
556 TH1I * h0 = new TH1I("hTOFESDs", "Number of matched TOF tracks per event;Number of TOF matched ESD tracks;Counts", 200, 0, 200) ;
557 TH1F * h1 = new TH1F("hTOFESDsTime", "Matched ESDs tracks: TOF Time spectrum; Calibrated TOF time [ns];Counts", 250, 0., 610. ) ;
558 TH1F * h2 = new TH1F("hTOFESDsRawTime", "Matched ESDs tracks: TOF raw Time spectrum;Measured TOF time [ns];Counts", 250, 0., 610.) ;
559 TH1F * h3 = new TH1F("hTOFESDsToT", "Matched ESDs tracks: TOF ToT spectrum; ESDs ToT [ns];Counts",100, 0., 48.8) ;
560 TH1F * h4 = new TH1F("hTOFESDskTOFOUT", "p_{T} distribution of tracks with kTOFout; p_{T} (GeV/c);Counts", 50, 0.20, 5.00) ;
561 TH1F * h5 = new TH1F("hTOFESDskTIME", "p_{T} distribution of tracks with kTOFout && kTIME; p_{T} (GeV/c);Counts", 50, 0.20, 5.00) ;
562 TH1F * h6 = new TH1F("hTOFESDsMatched", "p_{T} distribution of tracks with kTOFout && TOFtime>0; p_{T} (GeV/c);Counts", 50, 0.20, 5.00) ;
563 TH1F * h7 = new TH1F("hTOFESDsMatchingProb", "TPC-TOF track-matching probability;TOF matching probability (%) ;Counts",101, -1.0, 100) ;
564 TH1F * h8 = new TH1F("hTOFESDsDiffTime", "ESDs t_{TOF}-t_{exp,pi} spectrum in TOF (ps); t_{TOF}-t_{exp,pi} [ps];Counts", 200, -2440., 2440.) ;
565 TH1F * h9 = new TH1F("hTOFHitsLength", "Matched ESDs tracks: Length Spectrum; Track length [cm];Counts", 800, 0., 800) ;
566
567 h0->Sumw2() ;
568 h1->Sumw2() ;
569 h2->Sumw2() ;
570 h3->Sumw2() ;
571 h4->Sumw2() ;
572 h5->Sumw2() ;
573 h6->Sumw2() ;
574 h7->Sumw2() ;
575 h8->Sumw2() ;
576 h9->Sumw2() ;
577
578 Add2ESDsList(h0, 0, !expert, image) ;
579 Add2ESDsList(h1, 1, !expert, image) ;
580 Add2ESDsList(h2, 2, expert, !image) ;
581 Add2ESDsList(h3, 3, !expert, !image) ;
582 Add2ESDsList(h4, 4, expert, image) ;
583 Add2ESDsList(h5, 5, expert, image) ;
584 Add2ESDsList(h6, 6, expert, image) ;
585 Add2ESDsList(h7, 7, expert, image) ;
586 Add2ESDsList(h8, 8, expert, !image) ;
587 Add2ESDsList(h9, 9, !expert, !image) ;
588 //
589 ClonePerTrigClass(AliQAv1::kESDS); // this should be the last line
590}
591
592
593//____________________________________________________________________________
594void AliTOFQADataMakerRec::MakeRaws(AliRawReader* rawReader)
595{
596 //
597 // makes data from Raws
598 //
599 if (rawReader->GetType()==7) {
600
601 Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;//in ns
602 Double_t tot2ns=AliTOFGeometry::ToTBinWidth()*1E-3;
603 Int_t ntof[5]; /* 0=tot, 1=IA, 2=OA, 3=IC, 4=OC*/
604 for (Int_t j=0;j<5;j++){ ntof[j]=0;}
605 Int_t equipmentID[5]; //(ddl, trm, chain,tdc,channel)
606 Int_t volumeID[5]; //(sector,plate,strip,padX,padZ)
607 Int_t volumeID2[5]; //(sector,plate,strip,padZ,padX) to use AliTOFGeometry::GetIndex()
608 Int_t chIndex=-1;
609 Int_t indexCTTM[2]={-1,-1};
610 Int_t indexGeo2CTTM[2]={-1,-1};
611 Float_t macropadPhiTimeUPC[fgkFiredMacropadLimit][2];
612 for (Int_t ii=0;ii<2;ii++){
613 for (Int_t jj=0;jj<fgkFiredMacropadLimit;jj++){
614 macropadPhiTimeUPC[jj][ii]=-999.0;
615 }
616 }
617
618 TClonesArray * clonesRawData;
619 fTOFRawStream.SetRawReader(rawReader);
620 Int_t BCID=rawReader->GetBCID();
621
622 Int_t nFiredMacropad=0,
623 iFiredMacropad=-1;
624 nFiredMacropad=GetNumberOfFiredMacropad(rawReader);
625
626 //uncomment if needed to apply DeltaBC correction
627 //fTOFRawStream.ApplyBCCorrections(kTRUE);
628
629 if (fDecoderSummary){
630 fDecoderSummary->Reset();
631 }
632 for (Int_t iDDL = 0; iDDL < AliTOFGeometry::NDDL()*AliTOFGeometry::NSectors(); iDDL++){
633 rawReader->Reset();
634 fTOFRawStream.LoadRawDataBuffersV2(iDDL);
635
636 //* get decoding error counters
637 fDecoderSummary = ( (AliTOFDecoderV2*) fTOFRawStream.GetDecoderV2() )->GetDecoderSummaryData();
638 if ( (fDecoderSummary) && (fDecoderSummary ->GetErrorDetected()) ) {
639 Int_t errorSlotID=(Int_t) fDecoderSummary->GetErrorSlotID();
640 FillRawsData(18,iDDL,errorSlotID);
641 if (fDecoderSummary -> GetRecoverError() )
642 FillRawsData(18,iDDL,13);
643 }
644
645 clonesRawData = (TClonesArray*)fTOFRawStream.GetRawData();
646 for (Int_t iRawData = 0; iRawData<clonesRawData->GetEntriesFast(); iRawData++) {
647 AliTOFrawData *tofRawDatum = (AliTOFrawData*)clonesRawData->UncheckedAt(iRawData);
648 Float_t tofRawTime=tofRawDatum->GetTOF()*tdc2ns;
649
650 if (tofRawDatum->GetTOF()){
651 equipmentID[0]=iDDL;
652 equipmentID[1]=tofRawDatum->GetTRM();
653 equipmentID[2]=tofRawDatum->GetTRMchain();
654 equipmentID[3]=tofRawDatum->GetTDC();
655 equipmentID[4]=tofRawDatum->GetTDCchannel();
656
657 if (CheckEquipID(equipmentID)){
658 fTOFRawStream.EquipmentId2VolumeId(iDDL,
659 tofRawDatum->GetTRM(),
660 tofRawDatum->GetTRMchain(),
661 tofRawDatum->GetTDC(),
662 tofRawDatum->GetTDCchannel(),
663 volumeID);
664 //LTM data
665 if (FilterLTMData(equipmentID)) { //counts LTM hits
666 if (equipmentID[2]==1) //crate left, A-side or C-side
667 FillRawsData(15,equipmentID[0]);
668 else
669 FillRawsData(15,equipmentID[0]-1);
670
671 //retrieve CTTM index (Ltm, bit)
672 GetCTTMIndex(equipmentID, indexCTTM);
673
674 //get BX index within TOF-matching window
675 Int_t indexBC=-1;
676 indexBC= TMath::Nint(tofRawTime/24.4);
677
678 Int_t indexCttmChannel=indexCTTM[0]*24+indexCTTM[1];
679 FillRawsData(28,indexCttmChannel,indexBC);
680 FillRawsData(29,indexCttmChannel,tofRawTime);
681
682 //fired macropad map (from LTM hits) - only for low multi evts (UPC)
683 if ((nFiredMacropad<=fgCutNmaxFiredMacropad)){
684 iFiredMacropad++;
685 AliInfo(Form("Event found with %i fired macropads in BCID = %i!",nFiredMacropad,BCID));
686 FillRawsData(26,indexCTTM[0],indexCTTM[1]);
687 Float_t halfSMphi=-999.0;
688 if (indexCTTM[0]<36)
689 halfSMphi=indexCTTM[0]*10.+5.;
690 else halfSMphi=(indexCTTM[0]-36)*10.+5.;
691 macropadPhiTimeUPC[iFiredMacropad][0]=halfSMphi;
692 macropadPhiTimeUPC[iFiredMacropad][1]=indexBC;
693 }
694 }
695
696 //TRM data
697 if (CheckVolumeID(volumeID)){
698 volumeID2[0]=volumeID[0];
699 volumeID2[1]=volumeID[1];
700 volumeID2[2]=volumeID[2];
701 volumeID2[3]=volumeID[4];
702 volumeID2[4]=volumeID[3];
703 chIndex=AliTOFGeometry::GetIndex(volumeID2);
704
705 //fill hit map according to CTTM numbering
706 GetGeo2CTTMIndex(volumeID2, indexGeo2CTTM);
707 if ((nFiredMacropad<=fgCutNmaxFiredMacropad)){
708 FillRawsData(25,indexGeo2CTTM[0],indexGeo2CTTM[1]);
709 }
710 //hits selection
711 if (tofRawDatum->GetTOT()){
712 if (!(fCalibData->GetNoiseStatus(chIndex)==AliTOFChannelOnlineStatusArray::kTOFNoiseBad)
713 && (fCalibData->GetHWStatus(chIndex) == AliTOFChannelOnlineStatusArray::kTOFHWOk)) {//noise and enabled filter
714 ntof[0]++; //counter for tof hits
715
716 //fill global spectra for DQM plots
717 FillRawsData(5, tofRawTime) ;//in ns
718 FillRawsData(10, tofRawDatum->GetTOT()*tot2ns) ;//in ns
719 FillRawsData(23, BCID, tofRawTime) ;//in ns
720
721 //fill side-related spectra for experts plots
722 Int_t ddlACside=iDDL/36; // 0 or 1
723 Int_t ddlPerSm=iDDL%4;
724
725 if (volumeID2[0]>4 && volumeID2[0]<14){ //O side
726 if (ddlPerSm<2){ //A side
727 ntof[1]++;
728 FillRawsData(6, tofRawTime) ;
729 FillRawsData(11, tofRawDatum->GetTOT()*tot2ns) ;
730 } else { //C side
731 ntof[3]++;
732 FillRawsData(8, tofRawTime) ;
733 FillRawsData(13, tofRawDatum->GetTOT()*tot2ns) ;
734 }
735 } else {
736 if (volumeID2[0]<5 || volumeID2[0]>13){ //I side
737 if (ddlPerSm<2){ //A side
738 ntof[2]++;
739 FillRawsData(7, tofRawTime) ;
740 FillRawsData(12, tofRawDatum->GetTOT()*tot2ns) ;
741 } else {//C side
742 ntof[4]++;
743 FillRawsData(9, tofRawTime) ;
744 FillRawsData(14, tofRawDatum->GetTOT()*tot2ns) ;
745 }
746 }
747 }
748
749 //compute TRM offset
750 Int_t trm= iDDL*10+(equipmentID[1]-3);
751 FillRawsData(20+ddlACside,trm,tofRawTime);
752 FillRawsData(22,GetStripIndex(volumeID),tofRawTime) ;
753 Short_t fea = volumeID2[4]/12;
754 Float_t hitmapx = volumeID2[0] + ((Double_t)(3 - fea) + 0.5) *0.25;
755 FillRawsData(17,hitmapx,GetStripIndex(volumeID2));
756 if (fCalib.IsChannelEnabled(chIndex,kTRUE,kTRUE))//checks also if efficient and if problematic
757 FillRawsData(24,hitmapx,GetStripIndex(volumeID2));
758 }//noise filter
759 }//end hit selection
760 else { //orphans
761 if (!(fCalibData->GetNoiseStatus(chIndex) == AliTOFChannelOnlineStatusArray::kTOFNoiseBad)
762 && (fCalibData->GetHWStatus(chIndex) == AliTOFChannelOnlineStatusArray::kTOFHWOk))
763 FillRawsData(19, tofRawTime) ;//in ns
764 }//end orphans
765 }//end volumeID check
766 }//end equipID check
767 }//end tof check
768 }//loop on raw data
769 clonesRawData->Clear();
770 } // DDL Loop
771
772 for (Int_t j=0;j<5;j++) FillRawsData(j,ntof[j]);
773 fTOFRawStream.Clear();
774
775 if ((nFiredMacropad<=fgCutNmaxFiredMacropad)){
776 Float_t deltaPhiMacropad=-999.;
777 Float_t deltaTimeMacropad=-999.;
778 for (Int_t j=0;j<fgCutNmaxFiredMacropad+1; j++){
779 for (Int_t k=j+1;k<fgCutNmaxFiredMacropad+1; k++){
780 if ((macropadPhiTimeUPC[j][0]>0.0)&&(macropadPhiTimeUPC[k][0]>0.0)){
781 deltaPhiMacropad=TMath::Abs(macropadPhiTimeUPC[j][0]-macropadPhiTimeUPC[k][0]);
782 deltaTimeMacropad=TMath::Abs(macropadPhiTimeUPC[j][1]-macropadPhiTimeUPC[k][1]);
783 if (deltaPhiMacropad<=180.)
784 FillRawsData(27, deltaPhiMacropad,deltaTimeMacropad);
785 else
786 FillRawsData(27, TMath::Abs(360.0-deltaPhiMacropad),deltaTimeMacropad);
787 }
788 }
789 }
790 }//end cut on number of fired macropad
791 } else {
792 AliDebug(1,Form("Event of type %d found. Skipping non-physics event for QA.\n", rawReader->GetType()));
793 }
794
795 //fill reference map for DQM shifter only once in a detector cycle
796 if (fIsSOC) {
797 Int_t geoId[5]={-1,-1,-1,-1,-1};// pgeoId=(sm, mod, strip, padZ, padX)
798 Int_t detId[5]={-1,-1,-1,-1,-1};//detID=(ddl,trm,tdc, chain,channel)
799 Int_t trmIndex=-1;
800 for (Int_t ch = 0; ch < fCalibData->GetSize(); ch++) {
801 AliTOFGeometry::GetVolumeIndices(ch,geoId);
802 AliTOFRawStream::Geant2EquipmentId(geoId,detId);
803 if ((detId[1]<0)||(detId[0]<0)) continue;
804 trmIndex=(detId[1]-3)+detId[0]*10;
805
806 if ( (!(fCalibData->GetNoiseStatus(ch)==AliTOFChannelOnlineStatusArray::kTOFNoiseBad))
807 && (fCalibData->GetHWStatus(ch) == AliTOFChannelOnlineStatusArray::kTOFHWOk) ){
808 //fill reference map with info from OCDB
809 Short_t fea = geoId[4]/12;
810 Float_t hitmapx = geoId[0] + ((Double_t)(3 - fea) + 0.5)*0.25;
811 FillRawsData(16,hitmapx, GetStripIndex(geoId));
812 }
813 }
814 fIsSOC=kFALSE;
815 }
816
817 //enable options for DQM shifter
818 EnableDqmShifterOpt(kTRUE);
819 //
820 IncEvCountCycleRaws();
821 IncEvCountTotalRaws();
822 //
823}
824
825//____________________________________________________________________________
826void AliTOFQADataMakerRec::MakeRecPoints(TTree * clustersTree)
827{
828 //
829 // Make data from Clusters
830 //
831
832 Double_t tdc2ns=AliTOFGeometry::TdcBinWidth()*1E-3;
833 Double_t tot2ns=AliTOFGeometry::ToTBinWidth()*1E-3;
834
835 Int_t volumeID2[5];//(sm, plate,strip, padZ,padX)
836 //Int_t volumeID[5];//(sm, plate,strip, padX,padZ)
837
838 TBranch *branch=clustersTree->GetBranch("TOF");
839 if (!branch) {
840 AliError("can't get the branch with the TOF clusters !");
841 return;
842 }
843
844 static TClonesArray dummy("AliTOFcluster",10000);
845 dummy.Clear();
846 TClonesArray *clusters=&dummy;
847 branch->SetAddress(&clusters);
848
849 // Import the tree
850 clustersTree->GetEvent(0);
851
852 FillRecPointsData(0,(Int_t)clusters->GetEntriesFast()) ;
853
854 TIter next(clusters) ;
855 AliTOFcluster * c ;
856 while ( (c = dynamic_cast<AliTOFcluster *>(next())) ) {
857
858 // volumeID2[0] = c->GetDetInd(0);
859 // volumeID2[1] = c->GetDetInd(1);
860 // volumeID2[2] = c->GetDetInd(2);
861 // volumeID2[3] = c->GetDetInd(4); //padX
862 // volumeID2[4] = c->GetDetInd(3); //padZ
863
864 for (Int_t i=0;i<5;i++){
865 volumeID2[i]=c->GetDetInd(i); //X and Z indeces inverted in RecPoints
866 }
867 //Int_t chIndex=AliTOFGeometry::GetIndex(volumeID2);
868 Int_t iDDL=AliTOFRawStream::Geant2DDL(volumeID2);
869 Int_t iTRM=AliTOFRawStream::Geant2TRM(volumeID2);
870 Short_t fea = volumeID2[4]/12;
871 Float_t hitmapx = volumeID2[0] + ((Double_t)(3 - fea) + 0.5) *0.25;
872 Int_t ddlACside=iDDL/36; // 0 or 1
873 Int_t ddlPerSm=iDDL%4;
874
875 if ((c->GetTDCRAW()) && (c->GetTDC()) && (c->GetToT())){
876 if (volumeID2[0]>4 && volumeID2[0]<14){ //I side
877 if (ddlPerSm<2){ //A side
878 FillRecPointsData(1, c->GetTDC()*tdc2ns) ;//in ns
879 FillRecPointsData(5, c->GetTDCRAW()*tdc2ns) ;//in ns
880 FillRecPointsData(9, c->GetToT()*tot2ns) ;//in ns
881 } else {//C side
882 FillRecPointsData(3, c->GetTDC()*tdc2ns) ;//in ns
883 FillRecPointsData(7, c->GetTDCRAW()*tdc2ns) ;//in ns
884 FillRecPointsData(11, c->GetToT()*tot2ns) ;//in ns
885 }
886 } else {
887 if (volumeID2[0]<5 || volumeID2[0]>13){ //O side
888 if (ddlPerSm<2){ //A side
889 FillRecPointsData(2, c->GetTDC()*tdc2ns) ;//in ns
890 FillRecPointsData(6, c->GetTDCRAW()*tdc2ns) ;//in ns
891 FillRecPointsData(10, c->GetToT()*tot2ns) ;//in ns
892 } else { //C side
893 FillRecPointsData(4, c->GetTDC()*tdc2ns) ;//in ns
894 FillRecPointsData(8, c->GetTDCRAW()*tdc2ns) ;//in ns
895 FillRecPointsData(12, c->GetToT()*tot2ns) ;//in ns
896 }
897 }
898 }
899 FillRecPointsData(13,hitmapx,GetStripIndex(volumeID2));
900 FillRecPointsData(14,GetStripIndex(volumeID2), c->GetTDC()*tdc2ns) ;
901 Int_t trm= iDDL*10+(iTRM-3);
902 if (ddlACside==0) { //A side
903 FillRecPointsData(15,trm,c->GetTDC()*tdc2ns);
904 } else {//C side
905 FillRecPointsData(16,trm,c->GetTDC()*tdc2ns);
906 }
907 }//hit selection
908 }//end while
909 EnableDqmShifterOpt(kFALSE);
910 //
911 IncEvCountCycleRecPoints();
912 IncEvCountTotalRecPoints();
913 //
914}
915
916//____________________________________________________________________________
917void AliTOFQADataMakerRec::MakeESDs(AliESDEvent * esd)
918{
919 //
920 // make QA data from ESDs
921 //
922 const Double_t speedOfLight = TMath::C()*1E2*1E-12; // cm/ps
923 const Double_t pionMass = 0.13957018; //GeV/c^2
924
925 Int_t ntrk = esd->GetNumberOfTracks() ;
926 Int_t ntpc=0;
927 Int_t ntofout=0;
928
929 while (ntrk--) {
930 AliESDtrack *track=esd->GetTrack(ntrk);
931 Double_t tofTime=track->GetTOFsignal();//in ps
932 Double_t tofTimeRaw=track->GetTOFsignalRaw();//in ps
933 Double_t tofToT=track->GetTOFsignalToT(); //in ps
934
935 UInt_t status=track->GetStatus();
936 if (track->IsOn(AliESDtrack::kTPCrefit)) {
937 ntpc++;
938 Double_t y=track->Eta();
939 if (TMath::Abs(y)<0.9) { //select TOF acceptance
940 if ((status&AliESDtrack::kTOFout)!=0) { //define matching
941 ntofout++;
942 FillESDsData(1,tofTime*1E-3);
943 FillESDsData(2,tofTimeRaw*1E-3);
944 FillESDsData(3,tofToT*1E-3);
945 FillESDsData(4,track->Pt());
946
947 Double_t length =track->GetIntegratedLength();
948 Double_t mom2=(track->Pt()*track->Pt())+(track->Pz()*track->Pz());
949 Double_t piTexp = TMath::Sqrt(1+(pionMass*pionMass/mom2))*length/speedOfLight; //in ps
950 FillESDsData(8,tofTime-piTexp);
951 FillESDsData(9,length);
952
953 if ((status&AliESDtrack::kTIME)!=0)
954 FillESDsData(5,track->Pt());
955
956 if (tofTime>0)
957 FillESDsData(6,track->Pt());
958 } //end check on matched tracks
959 }
960 }//end check on TPCrefit
961 }
962
963 FillESDsData(0,ntofout) ;
964 if(ntpc>0){
965 Float_t ratio = (Float_t)ntofout/(Float_t)ntpc*100.; //matching probability
966 FillESDsData(7,ratio) ;
967 }
968
969 if(ntofout>0) {
970 Float_t ratio = (Float_t)ntofout/(Float_t)ntpc*100; //matched over propagated to TOF outer radius
971 FillESDsData(8,ratio) ;
972 }
973 EnableDqmShifterOpt(kFALSE);
974 //
975 IncEvCountCycleESDs();
976 IncEvCountTotalESDs();
977 //
978}
979
980//____________________________________________________________________________
981void AliTOFQADataMakerRec::StartOfDetectorCycle()
982{
983 //
984 //Detector specific actions at start of cycle
985 fCalibData = GetCalibData();
986 fIsSOC=kTRUE;
987 return;
988}
989
990//____________________________________________________________________________
991void AliTOFQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
992{
993 //Detector specific actions at end of cycle
994 // do the QA checking
995 ResetEventTrigClasses();
996 //
997 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
998 if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) continue ;
999 SetEventSpecie(AliRecoParam::ConvertIndex(specie));
1000
1001 for (int itc=-1;itc<GetNTrigClasses();itc++) { // RS: loop over eventual clones per trigger class
1002
1003 if (fEnableDqmShifterOpt) {
1004 // RS: fetch the histograms for given trigger class
1005 TObjArray& arrRW = *GetRawsDataOfTrigClass(itc);
1006
1007 // Help make the raw qa histogram easier to interpret for the DQM shifter
1008 if (!arrRW[ 0] || !arrRW[ 5] || !arrRW[10] || !arrRW[15] || !arrRW[16] || !arrRW[17]) continue;
1009
1010 printf("=========>Processed %i physics raw of specie %s with TrigGlass %d\n",
1011 GetEvCountCycleRaws(itc),AliRecoParam::GetEventSpecieName(specie), itc);
1012
1013 //Double_t monitorPeriodLength=fProcessedRawEventN*600*1E-9;//in s
1014
1015 if (fCalibData){
1016 //set minima and maxima to allow log scale
1017 Double_t yTimeMax = ((TH1*)arrRW[5])->GetMaximum()*1.05;
1018 Double_t yTotMax = ((TH1*)arrRW[10])->GetMaximum()*1.05;
1019 fLineExpTimeMin->SetY2(yTimeMax);
1020 fLineExpTimeMax->SetY2(yTimeMax);
1021 fLineExpTotMin->SetY2(yTotMax);
1022 fLineExpTotMax->SetY2(yTotMax);
1023 //
1024 for (Int_t j=0;j<18;j++){
1025 if ((j==0)||(j==5)||(j==10)||(j==15)||(j==16)||(j==17)) {
1026 TH1* htmp = (TH1*)arrRW[j];
1027 htmp->GetXaxis()->SetLabelOffset(0.005);
1028 htmp->GetXaxis()->SetLabelSize(0.05);
1029 htmp->GetXaxis()->SetTitleOffset(0.8);
1030 htmp->GetXaxis()->SetTitleSize(0.05);
1031 htmp->GetYaxis()->SetLabelOffset(0.005);
1032 htmp->GetYaxis()->SetLabelSize(0.06);
1033 htmp->GetYaxis()->SetTitleOffset(0.8);
1034 htmp->GetYaxis()->SetTitleSize(0.06);
1035 }
1036 }
1037 //make up for all histos
1038 for(Int_t j=0;j<5;j++) {
1039 TH1* htmp = (TH1*)arrRW[j];
1040 if (!htmp) continue;
1041 htmp->SetMarkerColor(kBlue);
1042 htmp->SetMarkerStyle(8);
1043 htmp->SetMarkerSize(0.7);
1044 }
1045 for(Int_t j=5;j<15;j++) {
1046 TH1* htmp = (TH1*)arrRW[j];
1047 if (!htmp) continue;
1048 htmp->SetLineColor(kBlue);
1049 htmp->SetLineWidth(1);
1050 htmp->SetMarkerColor(kBlue);
1051 }
1052
1053 TH1* htmp = (TH1*)arrRW[15];
1054 htmp->SetLineColor(kBlue);
1055 htmp->SetLineWidth(1);
1056 htmp->SetMarkerStyle(8);
1057 htmp->SetMarkerSize(0.7);
1058 htmp->SetMarkerColor(kBlue);//Option("bar");
1059 //
1060 TString title25 = Form("Map of hit pads according to CTTM numbering (Max Fired Macropad = %i)",fgCutNmaxFiredMacropad);
1061 TString title26 = Form("Map of hit macropads according to CTTM numbering (Max Fired Macropad = %i)",fgCutNmaxFiredMacropad);
1062 TString title27 = Form("#Deltat vs #Delta#Phi of hit macropads (Max Fired Macropad = %i)",fgCutNmaxFiredMacropad);
1063
1064 if ( (htmp=(TH1*)arrRW[16]) ) htmp->SetOption("colz");
1065 if ( (htmp=(TH1*)arrRW[17]) ) htmp->SetOption("colz");
1066 if ( (htmp=(TH1*)arrRW[18]) ) htmp->SetOption("colz");
1067 if ( (htmp=(TH1*)arrRW[22]) ) htmp->SetOption("colz");
1068 if ( (htmp=(TH1*)arrRW[23]) ) htmp->SetOption("colz");
1069 if ( (htmp=(TH1*)arrRW[24]) ) htmp->SetOption("colz");
1070 if ( (htmp=(TH1*)arrRW[28]) ) htmp->SetOption("colz");
1071 if ( (htmp=(TH1*)arrRW[29]) ) htmp->SetOption("colz");
1072
1073 if ( (htmp=(TH1*)arrRW[25]) ) {
1074 htmp->SetOption("colz");
1075 htmp->SetTitle(title25.Data());
1076 }
1077 if ( (htmp=(TH1*)arrRW[26]) ) {
1078 htmp->SetOption("colz");
1079 htmp->SetTitle(title26.Data());
1080 }
1081 if ( (htmp=(TH1*)arrRW[27]) ){
1082 htmp->SetOption("colz");
1083 htmp->SetTitle(title27.Data());
1084 }
1085 }
1086 }//END ENABLE DQM SHIFTER OPT
1087 } // RS: loop over trigger classes
1088 } //end for
1089 //
1090 AliQAChecker::Instance()->Run(AliQAv1::kTOF, task, list) ;
1091 //
1092}
1093//____________________________________________________________________________
1094void AliTOFQADataMakerRec::GetMapIndeces(const Int_t* const in , Int_t* out)
1095{
1096 //
1097 //return appropriate indeces for the theta-phi map
1098 //
1099
1100 Int_t npadX = AliTOFGeometry::NpadX();
1101 Int_t npadZ = AliTOFGeometry::NpadZ();
1102 Int_t nStripA = AliTOFGeometry::NStripA();
1103 Int_t nStripB = AliTOFGeometry::NStripB();
1104 Int_t nStripC = AliTOFGeometry::NStripC();
1105
1106 Int_t isector = in[0];
1107 Int_t iplate = in[1];
1108 Int_t istrip = in[2];
1109 Int_t ipadX = in[3];
1110 Int_t ipadZ = in[4];
1111
1112 Int_t stripOffset = 0;
1113 switch (iplate) {
1114 case 0:
1115 stripOffset = 0;
1116 break;
1117 case 1:
1118 stripOffset = nStripC;
1119 break;
1120 case 2:
1121 stripOffset = nStripC+nStripB;
1122 break;
1123 case 3:
1124 stripOffset = nStripC+nStripB+nStripA;
1125 break;
1126 case 4:
1127 stripOffset = nStripC+nStripB+nStripA+nStripB;
1128 break;
1129 default:
1130 AliDebug(1,Form("Wrong plate number in TOF (%d) !",iplate));
1131 break;
1132 };
1133 Int_t zindex=npadZ*(istrip+stripOffset)+(ipadZ+1);
1134 Int_t phiindex=npadX*isector+ipadX+1;
1135 out[0]=zindex;
1136 out[1]=phiindex;
1137
1138}
1139
1140//---------------------------------------------------------------
1141Int_t AliTOFQADataMakerRec::GetStripIndex(const Int_t * const in)
1142{
1143 /* return tof strip index between 0 and 91 */
1144
1145 Int_t nStripA = AliTOFGeometry::NStripA();
1146 Int_t nStripB = AliTOFGeometry::NStripB();
1147 Int_t nStripC = AliTOFGeometry::NStripC();
1148
1149 // Int_t isector = in[0];
1150 Int_t iplate = in[1];
1151 Int_t istrip = in[2];
1152 //Int_t ipadX = in[3];
1153 //Int_t ipadZ = in[4];
1154
1155 Int_t stripOffset = 0;
1156 switch (iplate) {
1157 case 0:
1158 stripOffset = 0;
1159 break;
1160 case 1:
1161 stripOffset = nStripC;
1162 break;
1163 case 2:
1164 stripOffset = nStripC+nStripB;
1165 break;
1166 case 3:
1167 stripOffset = nStripC+nStripB+nStripA;
1168 break;
1169 case 4:
1170 stripOffset = nStripC+nStripB+nStripA+nStripB;
1171 break;
1172 default:
1173 AliDebug(1,Form("Wrong plate number in TOF (%d) !",iplate));
1174 stripOffset=-1;
1175 break;
1176 };
1177
1178 if (stripOffset<0 || stripOffset>92) return -1;
1179 else
1180 return (stripOffset+istrip);
1181
1182}
1183//---------------------------------------------------------------
1184Bool_t AliTOFQADataMakerRec::CheckVolumeID(const Int_t * const volumeID)
1185{
1186 //
1187 //Checks volume ID validity
1188 //
1189 for (Int_t j=0;j<5;j++){
1190 if (volumeID[j]<0) {
1191 AliDebug(1,Form("Invalid detector volume index for volumeID[%i]",j));
1192 return kFALSE;
1193 }
1194 }
1195 return kTRUE;
1196}
1197
1198//---------------------------------------------------------------
1199Bool_t AliTOFQADataMakerRec::CheckEquipID(const Int_t * const equipmentID)
1200{
1201 //
1202 //Checks equipment ID validity
1203 for (Int_t j=0;j<5;j++){
1204 if (equipmentID[j]<0) {
1205 AliDebug(1,Form("Invalid equipment volume index for equipmentID[%i]",j));
1206 return kFALSE;
1207 }
1208 }
1209 return kTRUE;
1210}
1211//---------------------------------------------------------------
1212Bool_t AliTOFQADataMakerRec::FilterLTMData(const Int_t * const equipmentID) const
1213{
1214 /*It returns kTRUE if data come from LTM.
1215 It thus filters trigger-related signals */
1216
1217 Int_t ddl, trm, tdc;
1218 //if (!CheckEquipID(equipmentID)) return kFALSE;
1219 ddl = equipmentID[0];
1220 trm = equipmentID[1];
1221 tdc = equipmentID[3];
1222
1223 if ((ddl%2==1) && (trm==3) && (tdc>11 && tdc<15))
1224 return kTRUE;
1225 else
1226 return kFALSE;
1227
1228}
1229//---------------------------------------------------------------
1230Bool_t AliTOFQADataMakerRec::FilterSpare(const Int_t * const equipmentID) const
1231{
1232 /*It returns kTRUE if data come from spare
1233 equipment ID.
1234 So far only check on TRM 3 crate left is implemented */
1235
1236 Int_t ddl, trm, tdc;
1237 //if (!CheckEquipID(equipmentID)) return kFALSE;
1238 ddl = equipmentID[0];
1239 trm = equipmentID[1];
1240 tdc = equipmentID[3];
1241
1242 if ((ddl%2==1) && (trm==3) && (tdc>2 && tdc<12))
1243 return kTRUE;
1244 else
1245 return kFALSE;
1246}
1247
1248//-----------------------------------------------------------------------------
1249void AliTOFQADataMakerRec::GetGeo2LTMIndex(const Int_t * const detind, Int_t *indexLTM) {
1250 //
1251 // getting LTMmatrix indexes for current digit
1252 //
1253 Int_t stripId=GetStripIndex(detind);
1254
1255 if (detind[1]==0 || detind[1]==1 || (detind[1]==2 && detind[2]<=7)) { //A side
1256 if (detind[4]<24){ //R
1257 indexLTM[0] = detind[0]*2;
1258 } else { //L
1259 indexLTM[0] = detind[0]*2+1;
1260 }
1261 indexLTM[1]=stripId;
1262
1263 } else { //C side
1264 if (detind[4]<24){
1265 indexLTM[0] = detind[0]*2+36;
1266 } else {
1267 indexLTM[0] = (detind[0]*2+1)+36;
1268 }
1269 indexLTM[1]=90-stripId;
1270 }
1271
1272 // if (indexLTM[0]<36) { //A side
1273 // if (detind[1] ==0){
1274 // indexLTM[1] = detind[2];
1275 // }
1276 // else if (detind[1] ==1){
1277 // indexLTM[1] = detind[2]+nStripB;
1278 // }
1279 // else if (detind[1] ==2){
1280 // indexLTM[1] = detind[2]+19*2;
1281 // }
1282 // else{
1283 // AliError("Smth Wrong!!!");
1284 // }
1285 // }
1286 // else { //C side
1287 // if (detind[1]==2){
1288 // if (detind[4]<24)
1289 // indexLTM[1] = (nStripAL-detind[2])+nStripC+nStripB;
1290 // else
1291 // indexLTM[1] = (nStripAR-detind[2])+nStripC+nStripB;
1292 // }
1293 // else if (detind[1] ==3){
1294 // indexLTM[1] = (nStripB-detind[2])+nStripC;
1295 // }
1296 // else if (detind[1] ==4){
1297 // indexLTM[1] = nStripC-detind[2];
1298 // }
1299 // else{
1300 // AliError("Smth Wrong!!!");
1301 // }
1302 // }
1303}
1304
1305//-----------------------------------------------------------------------------
1306void AliTOFQADataMakerRec::GetGeo2CTTMIndex(Int_t *detind, Int_t *indexCTTM) {
1307 //
1308 // Returns CTTM index corresponding to the detector element detind
1309 //
1310 GetGeo2LTMIndex(detind,indexCTTM);
1311 indexCTTM[1]/=2;
1312 return;
1313}
1314
1315//-------------------------------------------------------------------------
1316Int_t AliTOFQADataMakerRec::GetNumberOfFiredMacropad(AliRawReader * rawReader){
1317
1318 Int_t nFired=0;
1319 TClonesArray * clonesRawData;
1320 if (!rawReader) return 0;
1321 for (Int_t iDDL = 0; iDDL < AliTOFGeometry::NDDL()*AliTOFGeometry::NSectors(); iDDL++){
1322 rawReader->Reset();
1323 fTOFRawStream.LoadRawDataBuffersV2(iDDL);
1324 clonesRawData = (TClonesArray*)fTOFRawStream.GetRawData();
1325 for (Int_t iRawData = 0; iRawData<clonesRawData->GetEntriesFast(); iRawData++) {
1326 AliTOFrawData *tofRawDatum = (AliTOFrawData*)clonesRawData->UncheckedAt(iRawData);
1327 if (tofRawDatum->GetTOF()){
1328 if ( (tofRawDatum->GetTRM()==3)&&
1329 (tofRawDatum->GetTDC()>11)&&
1330 (tofRawDatum->GetTDC()<15)) {
1331 nFired+=1;
1332 }
1333 }
1334 }
1335 }//loop over DDLs
1336 return nFired;
1337}
1338
1339//----------------------------------------------------------------
1340void AliTOFQADataMakerRec::GetCTTMIndex(Int_t *equipid, Int_t *indexCTTM) {
1341 //
1342 // Returns CTTM index corresponding to the equipment id equipid, only for LTM hits
1343 // equipid = (crate, trm, chain, tdc, channel)
1344
1345 if ( (equipid[1]!=3)||(equipid[3]<12) ){
1346 indexCTTM[0]=-1;
1347 indexCTTM[1]=-1;
1348 return;
1349 }
1350 Int_t modStrip2LTM[3][8]={ { 0, 1, 2, 3, 4, 5, 6, 7},
1351 { 8, 9, 10, 11, 12, 13, 14, 15},
1352 {16, 17, 18, 19, 20, 21, 22, 23}
1353 };
1354
1355 Int_t DDL2LTMmatrix[72]={0,1,37,36,2,3,39,38,4,5,41,40,6,7,43,42,8,9,45,44,10,11,47,46,12,13,49,48,14,15,51,50,16,17,53,52,18,19,
1356 55,54,20,21,57,56,22,23,59,58,24,25,61,60,26,27,63,62,28,29,65,64,30,31,67,66,32,33,69,68,34,35,71,70};
1357
1358 Int_t itdc=equipid[3]%12;
1359 Int_t crate=-1;
1360 if (equipid[2]==0)
1361 crate=equipid[0]-1;
1362 else crate=equipid[0];
1363
1364 indexCTTM[0]=DDL2LTMmatrix[crate];
1365 indexCTTM[1]=modStrip2LTM[itdc][equipid[4]];
1366 return;
1367}
1368
1369//_____________________________________________________________________________
1370void AliTOFQADataMakerRec::ReadHistogramRangeFromFile(const Char_t * filename)
1371{
1372 //
1373 // read histogram ranges from configuration file
1374 //
1375 if (!filename) {
1376 AliInfo("Config file with histograms ranges not found or invalid -> use default values.");
1377 SetDefaultHistogramRange();
1378 SetDefaultCutNmaxFiredMacropad();
1379 return;
1380 }
1381
1382 std::fstream configFile;
1383 configFile.open(filename, std::fstream::in);
1384 if (!configFile.is_open()){
1385 AliInfo("Cannot open config file with histograms ranges -> use default values.");
1386 SetDefaultHistogramRange();
1387 return;
1388 }
1389
1390 //check file size
1391 Int_t begin = configFile.tellg();
1392 configFile.seekg(0, std::fstream::end); /* end */
1393 Int_t end = configFile.tellg();
1394 Int_t size = end - begin;
1395 configFile.seekg(0, std::fstream::beg); /* rewind file */
1396 if (size <= 0){
1397 AliInfo(Form("Unexpected EOF of config file with histograms ranges. File size: %d -> use default values", size));
1398 SetDefaultHistogramRange();
1399 return;
1400 }
1401
1402 Int_t minMulti=9999, maxMulti=-9999;
1403 Int_t nbinsMulti=0,nbinsTime=0;
1404 Float_t minTime=9999.0, maxTime=-9999.0;
1405 Int_t cutFiredMacropad=0;
1406 TString endoflist;
1407 while (!configFile.eof()) {
1408 configFile >> cutFiredMacropad >> minMulti >> maxMulti >> minTime >> maxTime;
1409 configFile >> endoflist;
1410 if (endoflist.Contains("end")) break;
1411 }
1412
1413 //set multiplicity histo ranges
1414 if (minMulti>maxMulti){
1415 AliInfo("Invalid range for multiplicity histogram set. Changing to default values.");
1416 SetDefaultMultiHistogramRange();
1417 } else {
1418 nbinsMulti = maxMulti-minMulti;
1419 SetNbinsMultiplicityHisto(nbinsMulti);
1420 SetMultiplicityHistoRange(minMulti,maxMulti);
1421 //AliInfo(Form("Setting multiplicity histogram ranges to: multMin = %i - multMax = %i - nMultBins = %i", fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
1422 }
1423
1424 //set time histo ranges
1425 if (minTime>maxTime){
1426 AliInfo("Invalid range for time histogram set. Changing to default values.");
1427 SetDefaultTimeHistogramRange();
1428 } else {
1429 nbinsTime = TMath::Nint((maxTime - minTime)/fgkNbinsWidthTime);//ns
1430 maxTime=minTime+nbinsTime*fgkNbinsWidthTime;//ns
1431 SetNbinsTimeHisto(nbinsTime);
1432 SetTimeHistoRange(minTime,maxTime);
1433 //AliInfo(Form("Setting time histogram ranges to: timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i", fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
1434 }
1435
1436 if ((cutFiredMacropad>0)&&(cutFiredMacropad<fgkFiredMacropadLimit)){
1437 AliInfo("Invalid value for cut on fired macropad. Changing to default values.");
1438 SetDefaultCutNmaxFiredMacropad();
1439 } else {
1440 SetCutNmaxFiredMacropad(cutFiredMacropad);
1441 //AliInfo(Form("Setting cut on fired macropad to: = %i",cutFiredMacropad));
1442 }
1443 AliInfo(Form("Setting: multMin = %i - multMax = %i - nMultBins = %i, timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i, cutMaxFiredMacropad = %i",
1444 fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity, fgRangeMinTime, fgRangeMaxTime,fgNbinsTime, cutFiredMacropad));
1445 configFile.close();
1446 return;
1447}
1448
1449//_____________________________________________________________________________
1450void AliTOFQADataMakerRec::SetDefaultHistogramRange()
1451{
1452 //
1453 // set default histogram ranges (tuned on 2011 pp collisions)
1454 //
1455 AliInfo("Setting all histogram ranges to default values.");
1456 SetDefaultMultiHistogramRange();
1457 SetDefaultTimeHistogramRange();
1458 SetDefaultCutNmaxFiredMacropad();
1459 return;
1460}
1461
1462//_____________________________________________________________________________
1463void AliTOFQADataMakerRec::SetDefaultMultiHistogramRange()
1464{
1465 //
1466 // set default histogram ranges (tuned on 2011 pp collisions)
1467 //
1468 SetMultiplicityHistoRange (0, 200);
1469 SetNbinsMultiplicityHisto(200);
1470 AliInfo("Setting Multiplicity histogram ranges to default values.");
1471 AliInfo(Form("multMin = %i - multMax = %i - nMultBins = %i",
1472 fgRangeMinMultiplicity, fgRangeMaxMultiplicity, fgNbinsMultiplicity));
1473 return;
1474}
1475
1476//_____________________________________________________________________________
1477void AliTOFQADataMakerRec::SetDefaultTimeHistogramRange()
1478{
1479 //
1480 // set default histogram ranges (tuned on 2011 pp collisions)
1481 //
1482 SetNbinsTimeHisto(250);
1483 SetTimeHistoRange (0.0,610.);
1484
1485 AliInfo("Setting Time histogram ranges to default values:");
1486 AliInfo(Form("timeMin = %5.2f ns - timeMax = %5.2f ns - nTimeBins = %i",
1487 fgRangeMinTime, fgRangeMaxTime,fgNbinsTime));
1488 return;
1489}
1490
1491//------------------------------------------------------------------
1492void AliTOFQADataMakerRec::SetDefaultCutNmaxFiredMacropad()
1493{
1494 //
1495 // set default cut on fired macropad
1496 //
1497 SetCutNmaxFiredMacropad(5);
1498 AliInfo(Form("Setting cut on fired macropad to default values: NfiredMacropad = %i", fgCutNmaxFiredMacropad));
1499 return;
1500}