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