TOF/AliTOFQADataMakerRec.cxx

   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 /*
  25 Modified by fbellini on 01/11/2011
  26 - added histograms for LTM monitoring
  27 - fix for coverity
  28
  29 Modified 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
  35 Modified 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
  40 Modified 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
  54 Modified by fbellini on 13/09/2010
  55   - Set TLines as private members
  56   - Set image flag for expert histos
  57
  58 Modified 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
 110 ClassImp(AliTOFQADataMakerRec)
 111
 112 Int_t AliTOFQADataMakerRec::fgNbinsMultiplicity=200; //number of bins in multiplicity plot
 113 Int_t AliTOFQADataMakerRec::fgRangeMinMultiplicity=0;//min range in multiplicity plot
 114 Int_t AliTOFQADataMakerRec::fgRangeMaxMultiplicity=200;//max range in multiplicity plot
 115 Int_t AliTOFQADataMakerRec::fgNbinsTime=250;//number of bins in time plot
 116 const Float_t AliTOFQADataMakerRec::fgkNbinsWidthTime=2.44;//width of bins in time plot
 117 Float_t AliTOFQADataMakerRec::fgRangeMinTime=0.0;//range min in time plot
 118 Float_t AliTOFQADataMakerRec::fgRangeMaxTime=620.0; //range max in time plot
 119 Int_t AliTOFQADataMakerRec::fgCutNmaxFiredMacropad=5;//cut on number of max fired macropad
 120 const 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 //____________________________________________________________________________
 161 AliTOFQADataMakerRec::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 //__________________________________________________________________
 196 AliTOFQADataMakerRec& 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 //----------------------------------------------------------------------------
 207 AliTOFQADataMakerRec::~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   }
 228 for (Int_t sm=0;sm<22;sm++){
 229     if (fLineLTMbitId[sm])
 230       delete fLineLTMbitId[sm];
 231   }
 232 }
 233 //----------------------------------------------------------------------------
 234 AliTOFChannelOnlineStatusArray* 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 //____________________________________________________________________________
 278 void 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 //____________________________________________________________________________
 476 void 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 //____________________________________________________________________________
 547 void 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 //____________________________________________________________________________
 594 void 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 //____________________________________________________________________________
 826 void 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 //____________________________________________________________________________
 917 void 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 //____________________________________________________________________________
 981 void AliTOFQADataMakerRec::StartOfDetectorCycle()
 982 {
 983   //
 984   //Detector specific actions at start of cycle
 985   fCalibData = GetCalibData();
 986   fIsSOC=kTRUE;
 987   return;
 988 }
 989
 990 //____________________________________________________________________________
 991 void 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 //____________________________________________________________________________
1094 void 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 //---------------------------------------------------------------
1141 Int_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 //---------------------------------------------------------------
1184 Bool_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 //---------------------------------------------------------------
1199 Bool_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 //---------------------------------------------------------------
1212 Bool_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 //---------------------------------------------------------------
1230 Bool_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 //-----------------------------------------------------------------------------
1249 void 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 //-----------------------------------------------------------------------------
1306 void 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 //-------------------------------------------------------------------------
1316 Int_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 //----------------------------------------------------------------
1340 void 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 //_____________________________________________________________________________
1370 void 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 //_____________________________________________________________________________
1450 void 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 //_____________________________________________________________________________
1463 void 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 //_____________________________________________________________________________
1477 void 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 //------------------------------------------------------------------
1492 void 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 }