TOF/AliTOFPreprocessor.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 /* $Id$ */
  17
  18
  19 //#include <Riostream.h>
  20 //#include <stdio.h>
  21 //#include <stdlib.h>
  22
  23 #include <TFile.h>
  24 #include <TH2S.h>
  25 #include <TH1F.h>
  26 #include <TCanvas.h>
  27 #include <TMath.h>
  28 #include <TObjArray.h>
  29 #include <TObjString.h>
  30 #include <TTimeStamp.h>
  31
  32 #include "AliCDBMetaData.h"
  33 #include "AliCDBEntry.h"
  34 #include "AliLog.h"
  35 #include "AliTOFChannelOnlineArray.h"
  36 #include "AliTOFChannelOnlineStatusArray.h"
  37 #include "AliTOFDataDCS.h"
  38 #include "AliTOFDCSmaps.h"
  39 #include "AliTOFLvHvDataPoints.h"
  40 #include "AliTOFGeometry.h"
  41 #include "AliTOFPreprocessor.h"
  42 #include "AliTOFFEEReader.h"
  43 #include "AliTOFRawStream.h"
  44 #include "AliTOFCableLengthMap.h"
  45 #include "AliTOFcalibHisto.h"
  46 #include "AliTOFFEEDump.h"
  47 #include "TChain.h"
  48 #include "AliTOFDeltaBCOffset.h"
  49 #include "AliTOFCTPLatency.h"
  50 #include "AliTOFRunParams.h"
  51 #include "AliTOFT0FillOnlineCalib.h"
  52 #include "AliTOFHitField.h"
  53 #include "AliTOFChannelOffline.h"
  54 #include "TF1.h"
  55 #include "TGeoManager.h"
  56 #include "AliGeomManager.h"
  57 #include "AliTOFReadoutInfo.h"
  58
  59 // TOF preprocessor class.
  60 // It takes data from DCS and passes them to the class AliTOFDataDCS, which
  61 // processes them. The result is then written to the CDB.
  62 // Analogously, it takes data form DAQ (both at Run level and inclusive -
  63 // of all the runs - level, processes them, and stores both Reference Data
  64 // and Online Calibration files in the CDB.
  65 // Processing of Pulser/Noise Run data and of TOF FEE DCS map included also.
  66
  67 // return codes:
  68 // return=0 : all ok
  69 // return=1 : no DCS input data Map
  70 // return=2 : no DCS input data processing
  71 // return=3 : no DCS processed data was stored in Ref Data
  72 // return=4 : no DAQ input for Ref Data
  73 // return=5 : failed to store DAQ Ref Data
  74 // return=6 : failed to retrieve DAQ data for calibration
  75 // return=7 : problems in processing histos in the input DAQ file
  76 // return=8 : failed to store Online Delays
  77 // return=9 : failed to store Reference Data for Pulser
  78 // return=10: failed to retrieve Pulser data
  79 // return=11: failed to store Pulser map in OCDB
  80 // return=12: failed to store Reference Data for Noise
  81 // return=13: failed to retrieve Noise data
  82 // return=14: failed to store Noise map in OCDB
  83 // return=15: failed to retrieve FEE data from FXS
  84 // return=16: failed to retrieve FEE data from OCDB
  85 // return=17: failed to store FEE data in OCDB
  86 // return=18: failed to store FEE reference data in OCDB
  87 // return=20: failed in retrieving status variable
  88 // return=100 : no DCS input data Map (HV and LV status)
  89 // return=200 : no DCS input data processing (HV and LV status)
  90 // return=300 : no DCS processed data was stored in Ref Data (HV and LV status)
  91
  92 ClassImp(AliTOFPreprocessor)
  93
  94 const Int_t    AliTOFPreprocessor::fgkBinRangeAve = 13;    // number of bins where to calculate the mean
  95 const Double_t AliTOFPreprocessor::fgkIntegralThr = 100;   // min number of entries to perform computation of delay per channel
  96 const Double_t AliTOFPreprocessor::fgkThrPar      = 0.013; // parameter used to trigger the calculation of the delay
  97
  98 //_____________________________________________________________________________
  99
 100 AliTOFPreprocessor::AliTOFPreprocessor(AliShuttleInterface* shuttle) :
 101   AliPreprocessor("TOF", shuttle),
 102   fData(0),
 103   fHVLVmaps(0),
 104   fCal(0),
 105   fNChannels(0),
 106   fStoreRefData(kTRUE),
 107   fFDRFlag(kFALSE),
 108   fStatus(0),
 109   fMatchingWindow(0),
 110   fLatencyWindow(0),
 111   fIsStatusMapChanged(0)
 112 {
 113   // constructor
 114   AddRunType("PHYSICS");
 115   AddRunType("PULSER");
 116   AddRunType("NOISE");
 117
 118 }
 119
 120 //_____________________________________________________________________________
 121
 122 AliTOFPreprocessor::~AliTOFPreprocessor()
 123 {
 124   // destructor
 125 }
 126
 127 //______________________________________________________________________________
 128 void AliTOFPreprocessor::Initialize(Int_t run, UInt_t startTime,
 129         UInt_t endTime)
 130 {
 131   // Creates AliTOFDataDCS object
 132
 133   AliPreprocessor::Initialize(run, startTime, endTime);
 134
 135         AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s \n\tStartTime DCS Query %s \n\tEndTime DCS Query %s", run,
 136                 TTimeStamp(startTime).AsString(),
 137                 TTimeStamp(endTime).AsString(), ((TTimeStamp)GetStartTimeDCSQuery()).AsString(), ((TTimeStamp)GetEndTimeDCSQuery()).AsString()));
 138
 139         fData = new AliTOFDataDCS(fRun, fStartTime, fEndTime, GetStartTimeDCSQuery(), GetEndTimeDCSQuery());
 140         fHVLVmaps = new AliTOFLvHvDataPoints(fRun, fStartTime, fEndTime, GetStartTimeDCSQuery(), GetEndTimeDCSQuery());
 141         fNChannels = AliTOFGeometry::NSectors()*(2*(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())*AliTOFGeometry::NpadZ()*AliTOFGeometry::NpadX();
 142 }
 143 //_____________________________________________________________________________
 144 Bool_t AliTOFPreprocessor::ProcessDCS(){
 145
 146   // check whether DCS should be processed or not...
 147
 148   TString runType = GetRunType();
 149   Log(Form("RunType %s",runType.Data()));
 150
 151   if (runType != "PHYSICS"){
 152     return kFALSE;
 153   }
 154
 155   return kTRUE;
 156 }
 157 //_____________________________________________________________________________
 158
 159 UInt_t AliTOFPreprocessor::ProcessDCSDataPoints(TMap *dcsAliasMap)
 160 {
 161   // Fills data into a AliTOFDataDCS object
 162
 163
 164   Log("Processing DCS DP");
 165   TH1::AddDirectory(0);
 166
 167   Bool_t resultDCSMap=kFALSE;
 168   Bool_t resultDCSStore=kFALSE;
 169
 170   // processing DCS
 171
 172   fData->SetFDRFlag(fFDRFlag);
 173
 174   if (!dcsAliasMap){
 175     Log("No DCS map found: TOF exiting from Shuttle");
 176     if (fData){
 177             delete fData;
 178             fData = 0;
 179     }
 180     return 1;// return error Code for DCS input data not found
 181   }
 182   else {
 183
 184   // The processing of the DCS input data is forwarded to AliTOFDataDCS
 185     resultDCSMap=fData->ProcessData(*dcsAliasMap);
 186     if(!resultDCSMap){
 187       Log("Some problems occurred while processing DCS data, TOF exiting from Shuttle");
 188       if (fData){
 189               delete fData;
 190               fData = 0;
 191       }
 192       return 2;// return error Code for processed DCS data not stored
 193     }
 194     else{
 195       AliCDBMetaData metaDataDCS;
 196       metaDataDCS.SetBeamPeriod(0);
 197       metaDataDCS.SetResponsible("Chiara Zampolli");
 198       metaDataDCS.SetComment("This preprocessor fills an AliTOFDataDCS object.");
 199       AliInfo("Storing DCS Data");
 200       resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS);
 201       if (!resultDCSStore){
 202         Log("Some problems occurred while storing DCS data results in Reference Data, TOF exiting from Shuttle");
 203         if (fData){
 204                 delete fData;
 205                 fData = 0;
 206         }
 207         return 3;// return error Code for processed DCS data not stored
 208                  // in reference data
 209
 210       }
 211     }
 212   }
 213   if (fData){
 214           delete fData;
 215           fData = 0;
 216   }
 217
 218   return 0;
 219 }
 220 //_____________________________________________________________________________
 221
 222 UInt_t AliTOFPreprocessor::ProcessHVandLVdps(TMap *dcsAliasMap)
 223 {
 224   //
 225   //Fills data into a AliTOFLvHvDataPoints object
 226   // Merges fStatus object with LV and HV status at SOR
 227   // Updates fStatus object with LV and HV status
 228   //    at EOR in case of correct end of run
 229   //    at last but two value in case of end-of-run caused by TOF detector.
 230   //
 231
 232   Log("Processing HV and LV DCS DPs");
 233   TH1::AddDirectory(0);
 234
 235   Bool_t resultDCSMap=kFALSE;
 236
 237   // processing DCS HV and LV data points
 238
 239   fHVLVmaps->SetFDRFlag(fFDRFlag);
 240
 241   if (!dcsAliasMap){
 242     Log("No DCS map found: TOF exiting from Shuttle");
 243     if (fHVLVmaps){
 244       delete fHVLVmaps;
 245       fHVLVmaps = 0;
 246     }
 247     return 100;// return error Code for DCS input data not found
 248   }
 249   else {
 250
 251     // The processing of the DCS input data is forwarded to AliTOFDataDCS
 252     //if (0) { // AdC
 253     resultDCSMap = fHVLVmaps->ProcessData(*dcsAliasMap);
 254     if (!resultDCSMap) {
 255       Log("Some problems occurred while processing DCS data, TOF exiting from Shuttle");
 256       if (fHVLVmaps) {
 257         delete fHVLVmaps;
 258         fHVLVmaps = 0;
 259       }
 260       return 200;// return error Code for processed DCS data not stored
 261     }
 262     else {
 263
 264       // check with plots. Start...
 265       /*
 266       TH1F *hROsor = new TH1F("hROsor","RO status map at SOR",91*96*18,-0.5,91*96*18-0.5);
 267       for (Int_t ii=1; ii<=91*96*18; ii++) hROsor->SetBinContent(ii,-1);
 268       for (Int_t ii=0; ii<91*96*18; ii++) {
 269         if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWBad)
 270           hROsor->SetBinContent(ii+1,0);
 271         else if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWOk)
 272           hROsor->SetBinContent(ii+1,1);
 273       }
 274
 275       TH1F *hROandHVandLVsor = new TH1F("hROandHVandLVsor","RO.and.HV.andLV status map at SOR",91*96*18,-0.5,91*96*18-0.5);
 276       for (Int_t ii=1; ii<=91*96*18; ii++) hROandHVandLVsor->SetBinContent(ii,-1);
 277       TH1F *hROandHVandLVeor = new TH1F("hROandHVandLVeor","RO.and.HV.andLV status map at EOR",91*96*18,-0.5,91*96*18-0.5);
 278       for (Int_t ii=1; ii<=91*96*18; ii++) hROandHVandLVeor->SetBinContent(ii,-1);
 279       */
 280
 281       AliTOFDCSmaps * lvANDhvMap = (AliTOFDCSmaps*)fHVLVmaps->GetHVandLVmapAtSOR(); // Get LV.and.HV status map at SOR
 282       for (Int_t index=0; index<fNChannels; index++) {
 283         if ( ( lvANDhvMap->GetCellValue(index)==0 &&
 284                fStatus->GetHWStatus(index) != AliTOFChannelOnlineStatusArray::kTOFHWBad ) ||
 285              ( lvANDhvMap->GetCellValue(index)==1 &&
 286                fStatus->GetHWStatus(index) != AliTOFChannelOnlineStatusArray::kTOFHWOk ) ) {
 287           fStatus->SetHWStatus(index, AliTOFChannelOnlineStatusArray::kTOFHWBad);
 288           fIsStatusMapChanged=kTRUE;
 289         }
 290       }
 291
 292       // check with plots. Start...
 293       /*
 294       for (Int_t ii=0; ii<91*96*18; ii++) {
 295         if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWBad)
 296           hROandHVandLVsor->SetBinContent(ii+1,0);
 297         else if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWOk)
 298           hROandHVandLVsor->SetBinContent(ii+1,1);
 299       }
 300       */
 301
 302       lvANDhvMap = (AliTOFDCSmaps*)fHVLVmaps->GetHVandLVmapAtEOR(); // Get LV.and.HV status map at EOR
 303       for (Int_t index=0; index<fNChannels; index++) {
 304         if ( ( lvANDhvMap->GetCellValue(index)==0 &&
 305                fStatus->GetHWStatus(index)!=AliTOFChannelOnlineStatusArray::kTOFHWBad ) ||
 306              ( lvANDhvMap->GetCellValue(index)==1 &&
 307                fStatus->GetHWStatus(index) != AliTOFChannelOnlineStatusArray::kTOFHWOk ) ) {
 308           fStatus->SetHWStatus(index, AliTOFChannelOnlineStatusArray::kTOFHWBad);
 309           fIsStatusMapChanged=kTRUE;
 310         }
 311       }
 312
 313       // check with plots. Start...
 314       /*
 315       for (Int_t ii=0; ii<91*96*18; ii++) {
 316         if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWBad)
 317           hROandHVandLVeor->SetBinContent(ii+1,0);
 318         else if (fStatus->GetHWStatus(ii)==AliTOFChannelOnlineStatusArray::kTOFHWOk)
 319           hROandHVandLVeor->SetBinContent(ii+1,1);
 320       }
 321
 322       TCanvas *canvas = new TCanvas("canvas","",10,10,1000,1000);
 323       canvas->SetFillColor(0);
 324       canvas->Divide(2,2);
 325       canvas->cd(1);
 326       hROsor->SetLineWidth(2);
 327       hROsor->Draw();
 328       canvas->cd(2);
 329       hROandHVandLVsor->SetLineWidth(2);
 330       hROandHVandLVsor->Draw();
 331       canvas->cd(3);
 332       hROandHVandLVeor->SetLineWidth(2);
 333       hROandHVandLVeor->Draw();
 334       canvas->cd();
 335       */
 336
 337     }
 338     //} // AdC
 339   }
 340
 341
 342   /* check whether we don't need to update OCDB.
 343    * in this case we can return without errors. */
 344
 345   if (!fIsStatusMapChanged) {
 346     AliInfo("TOF HW status config has not changed. Do not overwrite stored file.");
 347     return 0; // return ok
 348   }
 349
 350   TString runType = GetRunType();
 351   if (runType != "PHYSICS") {
 352     AliInfo(Form("Run Type = %s, waiting to store status map",GetRunType()));
 353     return 0; // return ok
 354   }
 355
 356   // update the OCDB with the current FEE.and.HV.and.LV
 357   // since even a little difference has been detected.
 358
 359   AliCDBMetaData metaData;
 360   metaData.SetBeamPeriod(0);
 361   metaData.SetResponsible("Roberto Preghenella");
 362   metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object from FEE.and.HV.and.LV data.");
 363   AliInfo("Storing Status data from current run. Collected RO.and.HV.and.LV infos @ EOR");
 364   // store TOF channel status
 365   if (!Store("Calib", "Status", fStatus, &metaData, 0, kTRUE)) {
 366     // failed
 367     Log("problems while storing RO.and.HV.and.LV Status data object");
 368     if (fStatus){
 369       delete fStatus;
 370       fStatus = 0;
 371     }
 372     if (fHVLVmaps) {
 373       delete fHVLVmaps;
 374       fHVLVmaps = 0;
 375     }
 376     return 17; // return error code for problems while TOF channel status
 377   }
 378
 379   // everything fine. return
 380
 381   if (fStatus){
 382     delete fStatus;
 383     fStatus = 0;
 384   }
 385
 386   if (fHVLVmaps) {
 387     delete fHVLVmaps;
 388     fHVLVmaps = 0;
 389   }
 390
 391   return 0;
 392 }
 393
 394 //_____________________________________________________________________________
 395
 396 UInt_t AliTOFPreprocessor::ProcessOnlineDelays()
 397 {
 398   // Processing data from DAQ for online calibration
 399
 400   Bool_t updateOCDB = kFALSE;
 401   Log("Processing DAQ delays");
 402
 403   // reading configuration map
 404   TString compDelays = "kFALSE";
 405   Int_t deltaStartingRun = fRun;
 406   Int_t startingRun = fRun-deltaStartingRun;
 407   Int_t binRangeAve = fgkBinRangeAve;
 408   Double_t integralThr = fgkIntegralThr;
 409   Double_t thrPar = fgkThrPar;
 410
 411   AliCDBEntry *cdbEntry = GetFromOCDB("Calib","Config");
 412   if (!cdbEntry) {
 413           Log(Form("No Configuration entry found in CDB, using default values: ComputingDelays = %s, StartingRun = %i",compDelays.Data(), startingRun));
 414   }
 415   else {
 416           TMap *configMap = (TMap*)cdbEntry->GetObject();
 417           if (!configMap){
 418                   Log(Form("No map found in Config entry in CDB, using default values: ComputingDelays = %s, StartingRun = %i",compDelays.Data(), startingRun));
 419           }
 420           else{
 421                   TObjString *strDelays = (TObjString*)configMap->GetValue("ComputingDelays");
 422                   if (strDelays) {
 423                           compDelays = (TString) strDelays->GetString();
 424                   }
 425                   else {
 426                           Log(Form("No ComputingDelays value found in Map from Config entry in CDB, using default value: ComputingDelays =  %s",compDelays.Data()));
 427                   }
 428                   TObjString *strRun = (TObjString*)configMap->GetValue("StartingRun");
 429                   if (strRun) {
 430                           TString tmpstr = strRun->GetString();
 431                           startingRun = tmpstr.Atoi();
 432                           deltaStartingRun = fRun - startingRun;
 433                   }
 434                   else {
 435                           Log(Form("No StartingRun value found in Map from Config entry in CDB, using default value: StartingRun = %i",startingRun));
 436                   }
 437                   TObjString *strBinRangeAve = (TObjString*)configMap->GetValue("BinRangeAve");
 438                   if (strBinRangeAve) {
 439                           TString tmpstr = strBinRangeAve->GetString();
 440                           binRangeAve = tmpstr.Atoi();
 441                   }
 442                   else {
 443                           Log(Form("No BinRangeAve value found in Map from Config entry in CDB, using default value: BinRangeAve = %i",binRangeAve));
 444                   }
 445                   TObjString *strIntegralThr = (TObjString*)configMap->GetValue("IntegralThr");
 446                   if (strIntegralThr) {
 447                           TString tmpstr = strIntegralThr->GetString();
 448                           integralThr = tmpstr.Atof();
 449                   }
 450                   else {
 451                           Log(Form("No IntegralThr value found in Map from Config entry in CDB, using default value: IntegralThr = %f",integralThr));
 452                   }
 453                   TObjString *strThrPar = (TObjString*)configMap->GetValue("ThrPar");
 454                   if (strThrPar) {
 455                           TString tmpstr = strThrPar->GetString();
 456                           thrPar = tmpstr.Atof();
 457                   }
 458                   else {
 459                           Log(Form("No ThrPar value found in Map from Config entry in CDB, using default value: ThrPar = %f",thrPar));
 460                   }
 461           }
 462   }
 463   if (compDelays == "kTRUE") fFDRFlag = kFALSE;
 464   else fFDRFlag = kTRUE;
 465
 466   delete cdbEntry;
 467   cdbEntry = 0x0;
 468
 469   Log(Form("ComputingDelays = %s, StartingRun = %i",compDelays.Data(),startingRun));
 470
 471   /* init array with current calibration, if any */
 472   fCal = new AliTOFChannelOnlineArray(fNChannels);
 473   AliTOFChannelOnlineArray *curCal = NULL;
 474
 475   AliCDBEntry *cdbEntry2 = GetFromOCDB("Calib","ParOnlineDelay");
 476   if (!cdbEntry2 || !cdbEntry2->GetObject()) {
 477     /* no CDB entry found. set update flag */
 478     Log("     ************ WARNING ************");
 479     Log("No CDB ParOnlineDelay entry found, creating a new one!");
 480     Log("     *********************************");
 481     updateOCDB = kTRUE;
 482   }
 483   else {
 484     Log("Found previous ParOnlineDelay entry. Using it to init calibration");
 485     curCal = (AliTOFChannelOnlineArray *)cdbEntry2->GetObject();
 486     for (Int_t i = 0; i < fNChannels; i++)
 487       fCal->SetDelay(i, curCal->GetDelay(i));
 488   }
 489
 490
 491   TH1::AddDirectory(0);
 492
 493   Bool_t resultDAQRef=kFALSE;
 494   Bool_t resultTOFPP=kFALSE;
 495   TH2S *h2 = 0x0;
 496   // processing DAQ
 497
 498   TFile * daqFile=0x0;
 499
 500   if(fStoreRefData){
 501     //retrieving data at Run level
 502           TList* list = GetFileSources(kDAQ, "RUNLevel");
 503           if (list !=0x0 && list->GetEntries()!=0)
 504                   {
 505                           AliInfo("The following sources produced files with the id RUNLevel");
 506                           list->Print();
 507                           for (Int_t jj=0;jj<list->GetEntries();jj++){
 508                                   TObjString * str = dynamic_cast<TObjString*> (list->At(jj));
 509                                   if (!str) {
 510                                     AliError("dynamic_cast returned NULL");
 511                                     return 4;
 512                                   }
 513                                   AliInfo(Form("found source %s", str->String().Data()));
 514                                   // file to be stored run per run
 515                                   TString fileNameRun = GetFile(kDAQ, "RUNLevel", str->GetName());
 516                                   if (fileNameRun.Length()>0){
 517                                           AliInfo(Form("Got the file %s, now we can store the Reference Data for the current Run.", fileNameRun.Data()));
 518                                           daqFile = new TFile(fileNameRun.Data(),"READ");
 519                                           h2 = (TH2S*) daqFile->Get("htof");
 520                                           AliCDBMetaData metaDataHisto;
 521                                           metaDataHisto.SetBeamPeriod(0);
 522                                           metaDataHisto.SetResponsible("Chiara Zampolli");
 523                                           metaDataHisto.SetComment("This preprocessor stores the array of histos object as Reference Data.");
 524                                           AliInfo("Storing Reference Data");
 525                                           resultDAQRef = StoreReferenceData("Calib","DAQData",h2, &metaDataHisto);
 526                                           if (!resultDAQRef){
 527                                                   Log("some problems occurred::No Reference Data stored, TOF exiting from Shuttle");
 528                                                   delete h2;
 529                                                   delete list;
 530                                                   delete fCal;
 531                                                   fCal=0x0;
 532                                                   return 5;//return error code for failure in storing Ref Data
 533                                           }
 534                                           daqFile->Close();
 535                                           delete daqFile;
 536                                   }
 537
 538                                   else{
 539                                           Log("The input data file from DAQ (run-level) was not found, TOF exiting from Shuttle ");
 540                                           delete list;
 541                                           delete fCal;
 542                                           fCal=0x0;
 543                                           return 4;//return error code for failure in retrieving Ref Data
 544                                   }
 545                           }
 546                           delete list;
 547                   }
 548           else{
 549                   Log("The input data file list from DAQ (run-level) was not found, TOF exiting from Shuttle ");
 550                   delete fCal;
 551                   fCal=0x0;
 552                   return 4;//return error code for failure in retrieving Ref Data
 553           }
 554   }
 555
 556
 557   //Total files, with cumulative histos
 558
 559   TList* listTot = GetFileSources(kDAQ, "DELAYS");
 560   if (listTot !=0x0 && listTot->GetEntries()!=0)
 561           {
 562                   AliInfo("The following sources produced files with the id DELAYS");
 563                   listTot->Print();
 564                   for (Int_t jj=0;jj<listTot->GetEntries();jj++){
 565                           TObjString * str = dynamic_cast<TObjString*> (listTot->At(jj));
 566                           AliInfo(Form("found source %s", str->String().Data()));
 567
 568                           // file with summed histos, to extract calib params
 569                           TString fileName = GetFile(kDAQ, "DELAYS", str->GetName());
 570                           if (fileName.Length()>0){
 571                                   AliInfo(Form("Got the file %s, now we can extract some values.", fileName.Data()));
 572
 573                                   daqFile = new TFile(fileName.Data(),"READ");
 574                                   if (h2) delete h2;
 575                                   h2 = (TH2S*) daqFile->Get("htoftot");
 576                                   if (!h2){
 577                                           Log("some problems occurred:: No histo retrieved, TOF exiting from Shuttle");
 578                                           delete listTot;
 579                                           delete daqFile;
 580                                           delete fCal;
 581                                           fCal=0x0;
 582                                           return 7; //return error code for histograms not existing/junky
 583                                   }
 584                                   else {
 585                                           static const Int_t kSize=h2->GetNbinsX();
 586                                           static const Int_t kNBins=h2->GetNbinsY();
 587                                           static const Double_t kXBinmin=h2->GetYaxis()->GetBinLowEdge(1);
 588                                           if (kSize != fNChannels){
 589                                                   Log(" number of bins along x different from number of pads, found only a subset of the histograms, TOF exiting from Shuttle");
 590                                                   delete listTot;
 591                                                   delete h2;
 592                                                   delete daqFile;
 593                                                   delete fCal;
 594                                                   fCal=0x0;
 595                                                   return 7; //return error code for histograms not existing/junky
 596                                           }
 597                                           Int_t nNotStatistics = 0; // number of channel with not enough statistics
 598
 599                                           /* FDR flag set. do not compute delays, use nominal cable delays */
 600                                           if (fFDRFlag) {
 601
 602                                             Log(" Not computing delays according to flag set in Config entry in OCDB!");
 603                                             FillWithCosmicCalibration(fCal);
 604
 605                                             /* check whether the new calibration is different from the previous one */
 606                                             if (curCal) { /* well, check also whether we have a previous calibration */
 607                                               for (Int_t i = 0; i < fNChannels; i++) {
 608                                                 if (fCal->GetDelay(i) != curCal->GetDelay(i)) {
 609                                                   updateOCDB = kTRUE;
 610                                                   break;
 611                                                 }
 612                                               }
 613                                             }
 614                                             else /* otherwise update OCDB */
 615                                               updateOCDB = kTRUE;
 616
 617                                           }
 618
 619                                           else {  // computing delays if not in FDR runs
 620
 621                                             updateOCDB = kTRUE; /* always update OCDB when computing delays */
 622
 623                                                   for (Int_t ich=0;ich<kSize;ich++){
 624                                                           /* check whether channel has been read out during current run.
 625                                                            * if the status is bad it means it has not been read out.
 626                                                            * in this case skip channel in order to not affect the mean */
 627                                                           if (fStatus->GetHWStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFHWBad){
 628                                                                   AliDebug(2,Form(" Channel %i found bad according to FEEmap, (HW status = %i), skipping from delay computing",ich, (Int_t)fStatus->GetHWStatus(ich)));
 629                                                                   continue;
 630                                                           }
 631                                                           AliDebug(2,Form(" Channel %i found ok according to FEEmap, starting delay computing",ich));
 632                                                           TH1S *h1 = new TH1S("h1","h1",kNBins,kXBinmin-0.5,kNBins*1.+kXBinmin-0.5);
 633                                                           for (Int_t ibin=0;ibin<kNBins;ibin++){
 634                                                                   h1->SetBinContent(ibin+1,h2->GetBinContent(ich+1,ibin+1));
 635                                                           }
 636                                                           if(h1->Integral()<integralThr) {
 637                                                                   nNotStatistics++;
 638                                                                   Log(Form(" Not enough statistics for bin %i, skipping this channel",ich));  // printing message only if not in FDR runs
 639                                                                   delete h1;
 640                                                                   h1=0x0;
 641                                                                   continue;
 642                                                           }
 643                                                           Bool_t found=kFALSE;
 644                                                           Float_t minContent=h1->Integral()*thrPar;
 645                                                           Int_t nbinsX = h1->GetNbinsX();
 646                                                           Int_t startBin=1;
 647                                                           for (Int_t j=1; j<=nbinsX; j++){
 648                                                                   if ((
 649                                                                        h1->GetBinContent(j) +
 650                                                                        h1->GetBinContent(j+1)+
 651                                                                        h1->GetBinContent(j+2)+
 652                                                                        h1->GetBinContent(j+3))>minContent){
 653                                                                           found=kTRUE;
 654                                                                           startBin=j;
 655                                                                           break;
 656                                                                   }
 657                                                           }
 658                                                           if(!found) AliInfo(Form("WARNING!!! no start of fit found for histo # %i",ich));
 659                                                           // Now calculate the mean over the interval.
 660                                                           Double_t mean = 0;
 661                                                           Double_t sumw2 = 0;
 662                                                           Double_t nent = 0;
 663                                                           for(Int_t k=0;k<binRangeAve;k++){
 664                                                                   mean=mean+h1->GetBinCenter(startBin+k)*h1->GetBinContent(startBin+k);
 665                                                                   nent=nent+h1->GetBinContent(startBin+k);
 666                                                                   sumw2=sumw2+(h1->GetBinCenter(startBin+k))*(h1->GetBinCenter(startBin+k))*(h1->GetBinContent(startBin+k));
 667                                                           }
 668                                                           mean= mean/nent; //<x>
 669                                                           sumw2=sumw2/nent; //<x^2>
 670                                                           Double_t rmsmean= 0;
 671                                                           rmsmean = TMath::Sqrt((sumw2-mean*mean)/nent);
 672                                                           if (ich<fNChannels) {
 673                                                                   Float_t delay = mean*AliTOFGeometry::TdcBinWidth()*1.E-3; // delay in ns
 674                                                                   fCal->SetDelay(ich,delay);  // delay in ns
 675                                                                   AliDebug(2,Form("Setting delay %f (ns) for channel %i",delay,ich));
 676                                                           }
 677                                                           delete h1;
 678                                                           h1=0x0;
 679                                                   }
 680                                           }
 681                                           if (nNotStatistics!=0) Log(Form("Too little statistics for %d channels!",nNotStatistics));
 682                                   }
 683                                   delete h2;
 684                                   daqFile->Close();
 685                                   delete daqFile;
 686                           }
 687                           else{
 688                                   Log("The Cumulative data file from DAQ does not exist, TOF exiting from Shuttle");
 689                                   delete listTot;
 690                                   delete fCal;
 691                                   fCal=0x0;
 692                                   return 6;//return error code for problems in retrieving DAQ data
 693                           }
 694                   }
 695                   delete listTot;
 696           }
 697   else{
 698     Log("Problem: no list for Cumulative data file from DAQ was found, TOF exiting from Shuttle");
 699     delete fCal;
 700     fCal=0x0;
 701     return 6; //return error code for problems in retrieving DAQ data
 702   }
 703
 704   /* check whether we don't need to update OCDB.
 705    * in this case we can return without errors and
 706    * the current FEE is stored in the fStatus object. */
 707   if (!updateOCDB) {
 708     AliInfo("update OCDB flag not set. Do not overwrite stored file.");
 709     return 0; /* return ok */
 710   }
 711
 712   daqFile=0;
 713   AliCDBMetaData metaData;
 714   metaData.SetBeamPeriod(0);
 715   metaData.SetResponsible("Chiara Zampolli");
 716   metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineArray object for online calibration - delays.");
 717   AliInfo("Storing Calibration Data");
 718   resultTOFPP = Store("Calib","ParOnlineDelay",fCal, &metaData,deltaStartingRun,kTRUE);
 719   if(!resultTOFPP){
 720     Log("Some problems occurred while storing online object resulting from DAQ data processing");
 721     delete fCal;
 722     fCal=0x0;
 723     return 8;//return error code for problems in storing DAQ data
 724   }
 725
 726   if (fCal){
 727     delete fCal;
 728     fCal = 0;
 729   }
 730
 731   return 0;
 732 }
 733
 734 //_____________________________________________________________________________
 735
 736 UInt_t
 737 AliTOFPreprocessor::ProcessT0Fill()
 738 {
 739   // Processing data from DAQ for T0-fill measurement
 740
 741   Log("Processing T0-fill");
 742
 743 #if 0
 744   /* instance and setup CDB manager */
 745   AliCDBManager *cdb = AliCDBManager::Instance();
 746   /* load geometry */
 747   if (!gGeoManager) AliGeomManager::LoadGeometry();
 748 #endif
 749
 750   /* get params from OCDB */
 751   AliCDBEntry *cdbe = NULL;
 752
 753   /*
 754    * at this stage status object is not on OCDB yet
 755    * since it will be stored later. nevertheless we
 756    * should have the array in memory since it has been
 757    * already setup by ProcessFEF.
 758    */
 759
 760   /* check status and latency window available */
 761   if (!fStatus || !fLatencyWindow){
 762     AliError("No valid fStatus or fLatencyWindow found, some errors must have occurred!!");
 763     return 21;
 764   }
 765
 766   /* get offline calibration from OCDB */
 767   cdbe = GetFromOCDB("Calib", "ParOffline");
 768   if (!cdbe) {
 769     Log("cannot get \"ParOffline\" entry from OCDB");
 770     return 21;
 771   }
 772   TObjArray *offlineArray = (TObjArray *)cdbe->GetObject();
 773   AliTOFChannelOffline *channelOffline;
 774   if (!offlineArray) {
 775     Log("cannot get \"ParOffline\" object from CDB entry");
 776     return 21;
 777   }
 778   Log("got \"ParOffline\" object");
 779
 780   /* get deltaBC offset from OCDB */
 781   cdbe = GetFromOCDB("Calib", "DeltaBCOffset");
 782   if (!cdbe) {
 783     Log("cannot get \"DeltaBCOffset\" entry from OCDB");
 784     return 21;
 785   }
 786   AliTOFDeltaBCOffset *deltaBCOffsetObject = (AliTOFDeltaBCOffset *)cdbe->GetObject();
 787   if (!deltaBCOffsetObject) {
 788     Log("cannot get \"DeltaBCOffset\" object from CDB entry");
 789     return 21;
 790   }
 791   Int_t deltaBCOffset = deltaBCOffsetObject->GetDeltaBCOffset();
 792   Log(Form("got \"DeltaBCOffset\" object: deltaBCOffset=%d (BC bins)", deltaBCOffset));
 793
 794   /* get CTP latency from OCDB */
 795   cdbe = GetFromOCDB("Calib", "CTPLatency");
 796   if (!cdbe) {
 797     Log("cannot get \"CTPLatency\" entry from OCDB");
 798     return 21;
 799   }
 800   AliTOFCTPLatency *ctpLatencyObject = (AliTOFCTPLatency *)cdbe->GetObject();
 801   if (!ctpLatencyObject) {
 802     Log("cannot get \"CTPLatency\" object from CDB entry");
 803     return 21;
 804   }
 805   Float_t ctpLatency = ctpLatencyObject->GetCTPLatency();
 806   Log(Form("got \"CTPLatency\" object: ctpLatency=%f (ps)", ctpLatency));
 807
 808   /* get file sources from FXS */
 809   TList *fileList = GetFileSources(kDAQ, "HITS");
 810   if (!fileList || fileList->GetEntries() == 0) {
 811     Log("cannot get DAQ source file list or empty list");
 812     return 21;
 813   }
 814   Log(Form("got DAQ source file list: %d files", fileList->GetEntries()));
 815   fileList->Print();
 816
 817   /* create tree chain using file sources */
 818   TChain chain("hitTree");
 819   for (Int_t ifile = 0; ifile < fileList->GetEntries(); ifile++) {
 820     TObjString *str = (TObjString *)fileList->At(ifile);
 821     TString filename = GetFile(kDAQ, "HITS", str->GetName());
 822     chain.Add(filename);
 823     Log(Form("file added to input chain: source=%s, filename=%s", str->String().Data(), filename.Data()));
 824   }
 825   Int_t nhits = chain.GetEntries();
 826   Log(Form("input chain ready: %d hits", nhits));
 827
 828   /* setup input chain */
 829   AliTOFHitField *hit = new AliTOFHitField();
 830   chain.SetBranchAddress("hit", &hit);
 831
 832   /* create calib histo and geometry */
 833   AliTOFcalibHisto calibHisto;
 834   calibHisto.LoadCalibHisto();
 835   AliTOFGeometry tofGeo;
 836
 837   /* constants */
 838   Float_t c = TMath::C() * 1.e2 / 1.e12; /* cm/ps */
 839   Float_t c_1 = 1. / c;
 840   /* variables */
 841   Int_t index, timebin, totbin, deltaBC, l0l1latency, det[5];
 842   Float_t timeps, totns, corrps, length, timeexp, timezero, pos[3], latencyWindow;
 843
 844   /* histos */
 845   TH1F *hT0Fill = new TH1F("hT0Fill", "T0 fill;t - t_{exp}^{(c)} (ps);", 2000, -24400., 24400.);
 846
 847   /* loop over hits */
 848   for (Int_t ihit = 0; ihit < nhits; ihit++) {
 849
 850     /* get entry */
 851    chain.GetEntry(ihit);
 852
 853     /* get hit info */
 854     index = hit->GetIndex();
 855     timebin = hit->GetTimeBin();
 856     totbin = hit->GetTOTBin();
 857     deltaBC = hit->GetDeltaBC();
 858     l0l1latency = hit->GetL0L1Latency();
 859     latencyWindow = fLatencyWindow[index] * 1.e3;
 860
 861     /* convert time in ps and tot in ns */
 862     timeps = timebin * AliTOFGeometry::TdcBinWidth();
 863     totns = totbin * AliTOFGeometry::ToTBinWidth() * 1.e-3;
 864     /* get calibration correction in ps */
 865
 866
 867     channelOffline = (AliTOFChannelOffline *)offlineArray->At(index);
 868     if (totns < AliTOFGeometry::SlewTOTMin()) totns = AliTOFGeometry::SlewTOTMin();
 869     if (totns > AliTOFGeometry::SlewTOTMax()) totns = AliTOFGeometry::SlewTOTMax();
 870     corrps = 0.;
 871     for (Int_t ipar = 0; ipar < 6; ipar++) corrps += channelOffline->GetSlewPar(ipar) * TMath::Power(totns, ipar);
 872     corrps *= 1.e3;
 873     /* perform time correction */
 874     //    timeps = timeps + (deltaBC - deltaBCOffset) * AliTOFGeometry::BunchCrossingBinWidth() + l0l1latency * AliTOFGeometry::BunchCrossingBinWidth() + ctpLatency - latencyWindow - corrps; /* deltaBC correction removed for the time being */
 875     timeps = timeps + l0l1latency * AliTOFGeometry::BunchCrossingBinWidth() + ctpLatency - latencyWindow - corrps;
 876     /* compute length and expected time */
 877     tofGeo.GetVolumeIndices(index, det);
 878     tofGeo.GetPosPar(det, pos);
 879     length = 0.;
 880     for (Int_t i = 0; i < 3; i++) length += pos[i] * pos[i];
 881     length = TMath::Sqrt(length);
 882     timeexp = length * c_1;
 883     /* compute time zero */
 884     timezero = timeps - timeexp;
 885
 886     /* fill histos */
 887     hT0Fill->Fill(timezero);
 888   }
 889
 890   /* rebin until maximum bin has required minimum entries */
 891   Int_t maxBin = hT0Fill->GetMaximumBin();
 892   Float_t maxBinContent = hT0Fill->GetBinContent(maxBin);
 893   Float_t binWidth = hT0Fill->GetBinWidth(maxBin);
 894   while (maxBinContent < 400 && binWidth < 90.) {
 895     hT0Fill->Rebin(2);
 896     maxBin = hT0Fill->GetMaximumBin();
 897     maxBinContent = hT0Fill->GetBinContent(maxBin);
 898     binWidth = hT0Fill->GetBinWidth(maxBin);
 899   }
 900   Float_t maxBinCenter = hT0Fill->GetBinCenter(maxBin);
 901
 902   /* rough fit of the edge */
 903   TF1 *gaus = (TF1 *)gROOT->GetFunction("gaus");
 904   gaus->SetParameter(1, maxBinCenter);
 905   Float_t fitMin = maxBinCenter - 1000.; /* fit from 1 ns before max */
 906   Float_t fitMax = maxBinCenter + 1000.; /* fit until 1 ns above max */
 907   hT0Fill->Fit("gaus", "q0", "", fitMin, fitMax);
 908   /* better fit of the edge */
 909   Float_t mean, sigma;
 910   for (Int_t istep = 0; istep < 10; istep++) {
 911     mean = gaus->GetParameter(1);
 912     sigma = gaus->GetParameter(2);
 913     fitMin = mean - 3. * sigma;
 914     fitMax = mean;
 915     hT0Fill->Fit("gaus", "q0", "", fitMin, fitMax);
 916   }
 917   /* print params */
 918   mean = gaus->GetParameter(1);
 919   sigma = gaus->GetParameter(2);
 920   Float_t meane = gaus->GetParError(1);
 921   Float_t sigmae = gaus->GetParError(2);
 922   Log(Form("edge fit: mean  = %f +- %f ps", mean, meane));
 923   Log(Form("edge fit: sigma = %f +- %f ps", sigma, sigmae));
 924   /* check error */
 925   if (meane > 300.) {
 926     Log("error on mean is large: store default T0-fill value (0 ps)");
 927     mean = 0.;
 928   }
 929   if (sigmae > 300.) {
 930     Log("error on sigma is large: store default TOFreso value (200 ps)");
 931     sigma = 200.;
 932   }
 933
 934   /* scratch values from the fit and use max bin center as t0-fill */
 935   mean = maxBinCenter;
 936   sigma = -1.;
 937   Log(Form("do not care about fitted value, just use max bin as t0-fill: %f ps", mean));
 938
 939   /* create RunParams object */
 940   UInt_t timestamp[1] = {0};
 941   Float_t t0[1] = {mean};
 942   Float_t tofReso[1] = {sigma};
 943   Float_t t0Spread[1] = {-1.};
 944   AliTOFRunParams *runParamsObject = new AliTOFRunParams(1);
 945   runParamsObject->SetTimestamp(timestamp);
 946   runParamsObject->SetT0(t0);
 947   runParamsObject->SetTOFResolution(tofReso);
 948   runParamsObject->SetT0Spread(t0Spread);
 949
 950   /* store reference data */
 951   if(fStoreRefData){
 952     AliCDBMetaData metaDataHisto;
 953     metaDataHisto.SetBeamPeriod(0);
 954     metaDataHisto.SetResponsible("Roberto Preghenella");
 955     metaDataHisto.SetComment("online T0-fill histogram");
 956     if (!StoreReferenceData("Calib","T0Fill", hT0Fill, &metaDataHisto)) {
 957       Log("error while storing reference data");
 958       delete hT0Fill;
 959       delete hit;
 960       delete runParamsObject;
 961       return 21;
 962     }
 963     Log("reference data successfully stored");
 964   }
 965
 966   AliCDBMetaData metaData;
 967   metaData.SetBeamPeriod(0);
 968   metaData.SetResponsible("Roberto Preghenella");
 969   metaData.SetComment("online RunParams measurement");
 970   if (!Store("Calib", "RunParams", runParamsObject, &metaData, 0, kFALSE)) {
 971     Log("error while storing RunParams object");
 972     delete hT0Fill;
 973     delete hit;
 974     delete runParamsObject;
 975     return 21;
 976   }
 977   Log("RunParams object successfully stored");
 978
 979   delete hT0Fill;
 980   delete hit;
 981   delete runParamsObject;
 982   return 0;
 983
 984 }
 985
 986 //_____________________________________________________________________________
 987
 988 UInt_t
 989 AliTOFPreprocessor::ProcessNoiseCalibTrg()
 990 {
 991   // Processing data from DAQ using calibration triggers for noise measurement
 992
 993   Log("Processing Noise (calibration trigger)");
 994
 995   /* check status and matching window available */
 996   if (!fStatus || !fMatchingWindow){
 997     AliError("No valid fStatus or fMatchingWindow found, some errors must have occurred!!");
 998     return 22;
 999   }
1000
1001   Float_t noiseThr = 1000.;   // setting default threshold for noise to 1000 Hz
1002   // reading config map
1003   AliCDBEntry *cdbEntry = GetFromOCDB("Calib","ConfigNoise");
1004   if (!cdbEntry) Log(Form("No Configuration entry found in CDB, using default values: NoiseThr = %f",noiseThr));
1005   else {
1006     TMap *configMap = (TMap*)cdbEntry->GetObject();
1007     if (!configMap) Log(Form("No map found in Config entry in CDB, using default values: NoiseThr = %f", noiseThr));
1008     else {
1009       TObjString *strNoiseThr = (TObjString*)configMap->GetValue("NoiseThr");
1010       if (strNoiseThr) {
1011         TString tmpstr = strNoiseThr->GetString();
1012         noiseThr = tmpstr.Atoi();
1013       }
1014       else Log(Form("No NoiseThr value found in Map from ConfigNoise entry in CDB, using default value: NoiseThr = %f",noiseThr));
1015     }
1016   }
1017
1018   /* get file sources from FXS */
1019   TList *fileList = GetFileSources(kDAQ, "CALIB");
1020   if (!fileList || fileList->GetEntries() == 0) {
1021     Log("cannot get DAQ source file list or empty list");
1022     return 22;
1023   }
1024   Log(Form("got DAQ source file list: %d files", fileList->GetEntries()));
1025   fileList->Print();
1026
1027   /* open input file (only one expected) */
1028   TObjString *str = (TObjString *)fileList->At(0);
1029   TString filename = GetFile(kDAQ, "CALIB", str->GetName());
1030   Log(Form("opening input file: source=%s, filename=%s", str->String().Data(), filename.Data()));
1031   TFile *filein = TFile::Open(filename.Data());
1032   if (!filein || !filein->IsOpen()) {
1033     Log("cannot open input file");
1034     return 22;
1035   }
1036
1037   /* get histo from input file */
1038   TH1F *hCalibHit = (TH1F *)filein->Get("hCalibHit");
1039   if (!hCalibHit) {
1040     Log("cannot get \"hCalibHit\" histo");
1041     return 22;
1042   }
1043
1044   /* create and set noise rate histo and check rate */
1045   TH1F *hNoiseRate = new TH1F("hNoiseRate", ";index;rate (Hz)", fNChannels, 0., fNChannels);
1046   Float_t rate, rate_err;
1047   for (Int_t ich = 0; ich < fNChannels; ich++) {
1048     /* check channel enabled */
1049     if (fStatus->GetHWStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFHWBad) continue;
1050     /* set noise rate histo */
1051     rate = hCalibHit->GetBinContent(ich + 1);
1052     rate_err = hCalibHit->GetBinError(ich + 1);
1053     rate /= fMatchingWindow[ich] * 1.e-9;
1054     rate_err /= fMatchingWindow[ich] * 1.e-9;
1055     hNoiseRate->SetBinContent(ich + 1, rate);
1056     hNoiseRate->SetBinError(ich + 1, rate_err);
1057     /* check error */
1058     if (rate_err == 0.) continue;
1059     /* check noise rate and set noise flags */
1060     if ((rate - 3. * rate_err) > noiseThr) {
1061       Log(Form("channel %d detected as noisy: rate = (%f +- %f) Hz", ich, rate, rate_err));
1062       if (fStatus->GetNoiseStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFNoiseOk) {
1063         Log(Form("channel %d noise status changed from Ok to Bad", ich));
1064         fStatus->SetNoiseStatus(ich, AliTOFChannelOnlineStatusArray::kTOFNoiseBad);
1065         fIsStatusMapChanged = kTRUE;
1066       }
1067       else Log(Form("channel %d noise status unchanged", ich));
1068     }
1069     else if ((rate + 3. * rate_err) < noiseThr) {
1070       if (fStatus->GetNoiseStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFNoiseBad) {
1071         Log(Form("channel %d noise status changed from Bad to Ok", ich));
1072         fStatus->SetNoiseStatus(ich, AliTOFChannelOnlineStatusArray::kTOFNoiseOk);
1073         fIsStatusMapChanged = kTRUE;
1074       }
1075     }
1076   }
1077
1078   /* store reference data */
1079   if(fStoreRefData){
1080     AliCDBMetaData metaDataHisto;
1081     metaDataHisto.SetBeamPeriod(0);
1082     metaDataHisto.SetResponsible("Roberto Preghenella");
1083     metaDataHisto.SetComment("calibration trigger noise rate histogram");
1084     if (!StoreReferenceData("Calib","CalibNoise", hNoiseRate, &metaDataHisto)) {
1085       Log("error while storing reference data");
1086       delete hNoiseRate;
1087       filein->Close();
1088       return 22;
1089     }
1090     Log("reference data successfully stored");
1091   }
1092
1093   delete hNoiseRate;
1094   filein->Close();
1095   return 0;
1096 }
1097
1098 //_____________________________________________________________________________
1099
1100 UInt_t
1101 AliTOFPreprocessor::ProcessReadout()
1102 {
1103   // Processing data from DAQ to compute reaodut efficiency
1104
1105   Log("Processing Readout");
1106
1107   /* get file sources from FXS */
1108   TList *fileList = GetFileSources(kDAQ, "READOUT");
1109   if (!fileList || fileList->GetEntries() == 0) {
1110     Log("cannot get DAQ source file list or empty list");
1111     return 22;
1112   }
1113   Log(Form("got DAQ source file list: %d files", fileList->GetEntries()));
1114   fileList->Print();
1115
1116   /* open input file (only one expected) */
1117   TObjString *str = (TObjString *)fileList->At(0);
1118   TString filename = GetFile(kDAQ, "READOUT", str->GetName());
1119   Log(Form("opening input file: source=%s, filename=%s", str->String().Data(), filename.Data()));
1120   TFile *filein = TFile::Open(filename.Data());
1121   if (!filein || !filein->IsOpen()) {
1122     Log("cannot open input file");
1123     return 23;
1124   }
1125
1126   /* get histo from input file */
1127   TH1F *hChainEfficiency = (TH1F *)filein->Get("hChainEfficiency");
1128   if (!hChainEfficiency) {
1129     Log("cannot get \"hChainEfficiency\" histo");
1130     return 23;
1131   }
1132
1133   /* fill channel efficiency histo */
1134
1135  /* temporarly disable warnings */
1136   AliLog::EType_t logLevel = (AliLog::EType_t)AliLog::GetGlobalLogLevel();
1137   AliLog::SetGlobalLogLevel(AliLog::kError);
1138
1139   TH1F *hChannelEfficiency = new TH1F("hChannelEfficiency", "Channel readout efficiency;index;efficiency", fNChannels, 0., fNChannels);
1140   Int_t chainIndex, det[5], dummy, index;
1141   Float_t effi, effi_err;
1142   /* loop over DDLs */
1143   for (Int_t iddl = 0; iddl < 72; iddl++) {
1144     /* loop over TRMs */
1145     for (Int_t itrm = 0; itrm < 10; itrm++) {
1146       /* loop over chains */
1147       for (Int_t ichain = 0; ichain < 2; ichain++) {
1148         chainIndex = ichain + 2 * itrm + 20 * iddl;
1149         effi = hChainEfficiency->GetBinContent(chainIndex + 1);
1150         effi_err = hChainEfficiency->GetBinError(chainIndex + 1);
1151         /* loop over TDCs */
1152         for (Int_t itdc = 0; itdc < 15; itdc++) {
1153           /* loop over channels */
1154           for (Int_t ichannel = 0; ichannel < 8; ichannel++) {
1155
1156             /* get channel index */
1157             AliTOFRawStream::EquipmentId2VolumeId(iddl, itrm + 3, ichain, itdc, ichannel, det);
1158             dummy = det[4];
1159             det[4] = det[3];
1160             det[3] = dummy;
1161             /* check valid index */
1162             if (det[0] < 0 || det[0] > 17 ||
1163                 det[1] < 0 || det[1] > 5 ||
1164                 det[2] < 0 || det[2] > 18 ||
1165                 det[3] < 0 || det[3] > 1 ||
1166                 det[4] < 0 || det[4] > 47) continue;
1167             index = AliTOFGeometry::GetIndex(det);
1168
1169             /* set channel efficiency */
1170             hChannelEfficiency->SetBinContent(index + 1, effi);
1171             hChannelEfficiency->SetBinError(index + 1, effi_err);
1172
1173           }
1174         }
1175       }
1176     }
1177   }
1178
1179   /* re-enable warnings */
1180   AliLog::SetGlobalLogLevel(logLevel);
1181
1182   /* store reference data */
1183   if(fStoreRefData){
1184     /* setup TOF readout info object */
1185     AliTOFReadoutInfo readoutInfo;
1186     readoutInfo.SetChainEfficiency((TH1F *)filein->Get("hChainEfficiency"));
1187     readoutInfo.SetTRMData((TH1F *)filein->Get("hTRMData"));
1188     readoutInfo.SetTRMEmptyEvent((TH1F *)filein->Get("hTRMEmptyEvent"));
1189     readoutInfo.SetTRMBadEventCounter((TH1F *)filein->Get("hTRMBadEventCounter"));
1190     readoutInfo.SetTRMBadCRC((TH1F *)filein->Get("hTRMBadCRC"));
1191     readoutInfo.SetChainData((TH1F *)filein->Get("hChainData"));
1192     readoutInfo.SetChainBadStatus((TH1F *)filein->Get("hChainBadStatus"));
1193     readoutInfo.SetChainBadEventCounter((TH1F *)filein->Get("hChainBadEventCounter"));
1194     readoutInfo.SetTDCError((TH1F *)filein->Get("hTDCError"));
1195     readoutInfo.SetTDCErrorFlags((TH2F *)filein->Get("hTDCErrorFlags"));
1196
1197     AliCDBMetaData metaDataHisto;
1198     metaDataHisto.SetBeamPeriod(0);
1199     metaDataHisto.SetResponsible("Roberto Preghenella");
1200     metaDataHisto.SetComment("readout info data");
1201     if (!StoreReferenceData("Calib","ReadoutInfo", &readoutInfo, &metaDataHisto)) {
1202       Log("error while storing reference data");
1203       filein->Close();
1204       return 23;
1205     }
1206     Log("reference data successfully stored");
1207   }
1208
1209   AliCDBMetaData metaData;
1210   metaData.SetBeamPeriod(0);
1211   metaData.SetResponsible("Roberto Preghenella");
1212   metaData.SetComment("online ReadoutEfficiency measurement");
1213   if (!Store("Calib", "ReadoutEfficiency", hChannelEfficiency, &metaData, 0, kFALSE)) {
1214     Log("error while storing ReadoutEfficiency object");
1215     delete hChannelEfficiency;
1216     filein->Close();
1217     return 23;
1218   }
1219   Log("ReadoutEfficiency object successfully stored");
1220
1221   delete hChannelEfficiency;
1222   filein->Close();
1223   return 0;
1224 }
1225
1226 //_____________________________________________________________________________
1227
1228 UInt_t AliTOFPreprocessor::ProcessPulserData()
1229 {
1230   // Processing Pulser Run data for TOF channel status
1231
1232   Log("Processing Pulser");
1233
1234   if (fStatus==0x0){
1235           AliError("No valid fStatus found, some errors must have occurred!!");
1236           return 20;
1237   }
1238
1239   TH1::AddDirectory(0);
1240
1241   Bool_t resultPulserRef=kFALSE;
1242   Bool_t resultPulser=kFALSE;
1243
1244   static const Int_t kSize = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
1245   TH1S * htofPulser = new TH1S("hTOFpulser","histo with signals on TOF during pulser", kSize,-0.5,kSize-0.5);
1246   for (Int_t ibin =1;ibin<=kSize;ibin++){
1247           htofPulser->SetBinContent(ibin,-1);
1248   }
1249
1250   // processing pulser
1251
1252   TFile * daqFile=0x0;
1253   TH1S *h1=0x0;
1254
1255   //retrieving Pulser data
1256   TList* listPulser = GetFileSources(kDAQ, "PULSER");
1257   if (listPulser !=0x0 && listPulser->GetEntries()!=0)
1258           {
1259                   AliInfo("The following sources produced files with the id PULSER");
1260                   listPulser->Print();
1261                   Int_t nPulser = 0;
1262                   for (Int_t jj=0;jj<listPulser->GetEntries();jj++){
1263                           Int_t nPulserSource = 0;
1264                           TObjString * str = dynamic_cast<TObjString*> (listPulser->At(jj));
1265                           if (!str) {
1266                             AliError("dynamic_cast returned NULL");
1267                             return 4;
1268                           }
1269                           AliInfo(Form("found source %s", str->String().Data()));
1270                           // file to be stored run per run
1271                           TString fileNamePulser = GetFile(kDAQ, "PULSER", str->GetName());
1272                           if (fileNamePulser.Length()>0){
1273                                   // storing refernce data
1274                                   AliInfo(Form("Got the file %s, now we can process pulser data.", fileNamePulser.Data()));
1275                                   daqFile = new TFile(fileNamePulser.Data(),"READ");
1276                                   h1 = (TH1S*) daqFile->Get("hTOFpulser");
1277                                   for (Int_t ibin=0;ibin<kSize;ibin++){
1278                                           if ((h1->GetBinContent(ibin+1))!=-1){
1279                                                   if ((htofPulser->GetBinContent(ibin+1))==-1){
1280                                                           htofPulser->SetBinContent(ibin+1,h1->GetBinContent(ibin+1));
1281                                                   }
1282                                                   else {
1283                                                           Log(Form("Something strange occurred during Pulser run, channel %i already read by another LDC, please check!",ibin));
1284                                                   }
1285                                           }
1286                                   }
1287
1288                                   // elaborating infos
1289                                   Double_t mean =0;
1290                                   Int_t nread=0;
1291                                   Int_t nreadNotEmpty=0;
1292                                   for (Int_t ientry=1;ientry<=h1->GetNbinsX();ientry++){
1293
1294                                           AliDebug(3,Form(" channel %i pulser status before pulser = %i, with global status = %i",ientry,(Int_t)fStatus->GetPulserStatus(ientry),(Int_t)fStatus->GetStatus(ientry)));
1295                                           /* check whether channel has been read out during current run.
1296                                            * if the status is bad it means it has not been read out.
1297                                            * in this case skip channel in order to not affect the mean */
1298                                           if (fStatus->GetHWStatus(ientry-1) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
1299                                                   continue;
1300                                           nPulser++;
1301                                           nPulserSource++;
1302                                           if (h1->GetBinContent(ientry)==-1) continue;
1303                                           else {
1304                                                   if (h1->GetBinContent(ientry)>0) {
1305                                                           nreadNotEmpty++;
1306                                                           AliDebug(2,Form(" channel %i is ok with entry = %f; so far %i channels added ",ientry-1,h1->GetBinContent(ientry),nreadNotEmpty));
1307                                                   }
1308                                                   mean+=h1->GetBinContent(ientry);
1309                                                   nread++;
1310                                           }
1311                                   }
1312                                   if (nread!=0) {
1313                                           mean/=nread;
1314                                           AliDebug(2,Form(" nread =  %i , nreadNotEmpty = %i, mean = %f",nread,nreadNotEmpty,mean));
1315                                           for (Int_t ich =0;ich<fNChannels;ich++){
1316                                                   if (h1->GetBinContent(ich+1)==-1) continue;
1317                                                   AliDebug(3,Form(" channel %i pulser status before pulser = %i",ich,(Int_t)fStatus->GetPulserStatus(ich)));
1318
1319                                                   /* check whether channel has been read out during current run.
1320                                                    * if the status is bad it means it has not been read out.
1321                                                    * in this case skip channel in order to leave its status
1322                                                    * unchanged */
1323                                                   if (fStatus->GetHWStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
1324                                                           continue;
1325
1326                                                   if (h1->GetBinContent(ich+1)<0.05*mean){
1327                                                           fStatus->SetPulserStatus(ich,AliTOFChannelOnlineStatusArray::kTOFPulserBad);  // bad status for pulser
1328                                                           AliDebug(2,Form( " channel %i pulser status after pulser = %i (bad, content = %f), with global status = %i",ich,(Int_t)fStatus->GetPulserStatus(ich),h1->GetBinContent(ich+1),(Int_t)fStatus->GetStatus(ich)));
1329                                                   }
1330                                                   else {
1331                                                           fStatus->SetPulserStatus(ich,AliTOFChannelOnlineStatusArray::kTOFPulserOk);  // good status for pulser
1332                                                           AliDebug(2,Form( " channel %i pulser status after pulser = %i (good), with global status = %i",ich,(Int_t)fStatus->GetPulserStatus(ich),(Int_t)fStatus->GetStatus(ich)));
1333                                                   }
1334                                           }
1335                                   }
1336                                   else {
1337                                           Log("No channels read!! No action taken, keeping old status");
1338                                   }
1339
1340                                   daqFile->Close();
1341                                   delete daqFile;
1342                                   delete h1;
1343                           }
1344
1345                           else{
1346                                   Log("The input data file from DAQ (pulser) was not found, TOF exiting from Shuttle ");
1347                                   delete listPulser;
1348                                   delete htofPulser;
1349                                   htofPulser = 0x0;
1350                                   if (fStatus){
1351                                           delete fStatus;
1352                                           fStatus = 0;
1353                                   }
1354                                   return 10;//return error code for failure in retrieving Ref Data
1355                           }
1356                           AliDebug(2,Form(" Number of channels processed during pulser run from source %i = %i",jj, nPulserSource));
1357                   }
1358                   AliDebug(2,Form(" Number of channels processed during pulser run = %i",nPulser));
1359                   delete listPulser;
1360           }
1361
1362   else{
1363           Log("The input data file list from DAQ (pulser) was not found, TOF exiting from Shuttle ");
1364           delete htofPulser;
1365           htofPulser = 0x0;
1366           if (fStatus){
1367                   delete fStatus;
1368                   fStatus = 0;
1369           }
1370           return 10;//return error code for failure in retrieving Ref Data
1371   }
1372
1373   //storing in OCDB
1374
1375   AliCDBMetaData metaData;
1376   metaData.SetBeamPeriod(0);
1377   metaData.SetResponsible("Chiara Zampolli");
1378   metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object after a Pulser run.");
1379   AliInfo("Storing Calibration Data from Pulser Run");
1380   resultPulser = Store("Calib","Status",fStatus, &metaData,0,kTRUE);
1381   if(!resultPulser){
1382     Log("Some problems occurred while storing online object resulting from Pulser data processing");
1383     delete htofPulser;
1384     htofPulser = 0x0;
1385     if (fStatus){
1386             delete fStatus;
1387             fStatus = 0;
1388     }
1389     return 11;//return error code for problems in storing Pulser data
1390   }
1391
1392   if(fStoreRefData){
1393
1394     AliCDBMetaData metaDataHisto;
1395     metaDataHisto.SetBeamPeriod(0);
1396     metaDataHisto.SetResponsible("Chiara Zampolli");
1397     metaDataHisto.SetComment("This preprocessor stores the Ref data from a pulser run.");
1398     AliInfo("Storing Reference Data");
1399     resultPulserRef = StoreReferenceData("Calib","PulserData",htofPulser, &metaDataHisto);
1400     if (!resultPulserRef){
1401       Log("some problems occurred::No Reference Data for pulser stored, TOF exiting from Shuttle");
1402       delete htofPulser;
1403       htofPulser = 0x0;
1404       if (fStatus){
1405               delete fStatus;
1406               fStatus = 0;
1407       }
1408       return 9;//return error code for failure in storing Ref Data
1409     }
1410   }
1411
1412   daqFile=0;
1413
1414   delete htofPulser;
1415   htofPulser = 0x0;
1416
1417   if (fStatus){
1418     delete fStatus;
1419     fStatus = 0;
1420   }
1421
1422   return 0;
1423 }
1424 //_____________________________________________________________________________
1425
1426 UInt_t AliTOFPreprocessor::ProcessNoiseData()
1427 {
1428
1429   // Processing Noise Run data for TOF channel status
1430
1431   Log("Processing Noise");
1432
1433   if (fStatus==0x0){
1434           AliError("No valid fStatus found, some errors must have occurred!!");
1435           return 20;
1436   }
1437
1438   Float_t noiseThr = 1;   // setting default threshold for noise to 1 Hz
1439   // reading config map
1440   AliCDBEntry *cdbEntry = GetFromOCDB("Calib","ConfigNoise");
1441   if (!cdbEntry) {
1442           Log(Form("No Configuration entry found in CDB, using default values: NoiseThr = %f",noiseThr));
1443   }
1444   else {
1445           TMap *configMap = (TMap*)cdbEntry->GetObject();
1446           if (!configMap){
1447                   Log(Form("No map found in Config entry in CDB, using default values: NoiseThr = %f", noiseThr));
1448           }
1449           else{
1450                   TObjString *strNoiseThr = (TObjString*)configMap->GetValue("NoiseThr");
1451                   if (strNoiseThr) {
1452                           TString tmpstr = strNoiseThr->GetString();
1453                           noiseThr = tmpstr.Atoi();
1454                   }
1455                   else {
1456                           Log(Form("No NoiseThr value found in Map from ConfigNoise entry in CDB, using default value: NoiseThr = %f",noiseThr));
1457                   }
1458           }
1459   }
1460
1461   delete cdbEntry;
1462   cdbEntry = 0x0;
1463
1464   TH1::AddDirectory(0);
1465
1466   Bool_t resultNoiseRef=kFALSE;
1467   Bool_t resultNoise=kFALSE;
1468
1469   static const Int_t kSize = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
1470   TH1F * htofNoise = new TH1F("hTOFnoise","histo with signals on TOF during noise", kSize,-0.5,kSize-0.5);
1471   for (Int_t ibin =1;ibin<=kSize;ibin++){
1472           htofNoise->SetBinContent(ibin,-1);
1473   }
1474
1475   // processing noise
1476
1477   TFile * daqFile=0x0;
1478   TH1F * h1=0x0;
1479
1480   // useful counters
1481   Int_t nNoise = 0;
1482   Int_t nNoisyChannels = 0;
1483   Int_t nNotNoisyChannels = 0;
1484   Int_t nChannelsFromDA = 0;
1485   Int_t nMatchingWindowNullNonZero = 0;
1486   Int_t nMatchingWindowNullEqualZero = 0;
1487
1488   // retrieving Noise data
1489   TList* listNoise = GetFileSources(kDAQ, "NOISE");
1490   if (listNoise !=0x0 && listNoise->GetEntries()!=0)
1491           {
1492                   AliInfo("The following sources produced files with the id NOISE");
1493                   listNoise->Print();
1494                   for (Int_t jj=0;jj<listNoise->GetEntries();jj++){
1495                           Int_t nNoiseSource = 0;
1496                           TObjString * str = dynamic_cast<TObjString*> (listNoise->At(jj));
1497                           if (!str) {
1498                             AliError("dynamic_cast returned NULL");
1499                             return 4;
1500                           }
1501                           AliInfo(Form("found source %s", str->String().Data()));
1502                           // file to be stored run per run
1503                           TString fileNameNoise = GetFile(kDAQ, "NOISE", str->GetName());
1504                           if (fileNameNoise.Length()>0){
1505                                   // storing reference data
1506                                   AliInfo(Form("Got the file %s, now we can process noise data.", fileNameNoise.Data()));
1507                                   daqFile = new TFile(fileNameNoise.Data(),"READ");
1508                                   h1 = (TH1F*) daqFile->Get("hTOFnoise");
1509                                   for (Int_t ibin=0;ibin<kSize;ibin++){
1510                                           if ((h1->GetBinContent(ibin+1))!=-1){
1511                                                   nNoiseSource++;
1512                                                   // checking the matching window for current channel
1513                                                   if (fMatchingWindow[ibin] == 0){
1514                                                           Log(Form("Matching window for channel %i null, but the channel was read by the LDC! skipping channel, BUT Please check!",ibin));
1515                                                           if ((h1->GetBinContent(ibin+1))!=0) nMatchingWindowNullNonZero++;
1516                                                           if ((h1->GetBinContent(ibin+1))==0) nMatchingWindowNullEqualZero++;
1517                                                           continue;
1518                                                   }
1519                                                   if ((htofNoise->GetBinContent(ibin+1))==-1){
1520                                                           htofNoise->SetBinContent(ibin+1,h1->GetBinContent(ibin+1)/(fMatchingWindow[ibin]*1.E-9));
1521                                                           if ((h1->GetBinContent(ibin+1))!= 0) AliDebug(2,Form("Channel = %i, Matching window = %i, Content = %f", ibin, fMatchingWindow[ibin], htofNoise->GetBinContent(ibin+1)));
1522                                                   }
1523                                                   else {
1524                                                           Log(Form("Something strange occurred during Noise run, channel %i already read by another LDC, please check!",ibin));
1525                                                   }
1526                                           }
1527                                   }
1528
1529                                   Log(Form(" Number of channels processed during noise run from source %i = %i",jj, nNoiseSource));
1530                                   daqFile->Close();
1531                                   delete daqFile;
1532                                   delete h1;
1533                                   daqFile = 0x0;
1534                                   h1 = 0x0;
1535
1536                           }
1537                           else{
1538                                   Log("The input data file from DAQ (noise) was not found, TOF exiting from Shuttle ");
1539                                   delete listNoise;
1540                                   listNoise = 0x0;
1541                                   delete htofNoise;
1542                                   htofNoise = 0x0;
1543                                   if (fStatus){
1544                                           delete fStatus;
1545                                           fStatus = 0;
1546                                   }
1547                                   if (fMatchingWindow){
1548                                           delete [] fMatchingWindow;
1549                                           fMatchingWindow = 0;
1550                                   }
1551                                   return 13;//return error code for failure in retrieving Ref Data
1552                           }
1553                   }
1554           }
1555
1556   else{
1557           Log("The input data file list from DAQ (noise) was not found, TOF exiting from Shuttle ");
1558           delete htofNoise;
1559           htofNoise = 0x0;
1560           if (fStatus){
1561                   delete fStatus;
1562                   fStatus = 0;
1563           }
1564           if (fMatchingWindow){
1565                   delete [] fMatchingWindow;
1566                   fMatchingWindow = 0;
1567           }
1568           return 13;//return error code for failure in retrieving Ref Data
1569   }
1570
1571   // elaborating infos to set NOISE status
1572   for (Int_t ich =0;ich<fNChannels;ich++){
1573           if (htofNoise->GetBinContent(ich+1)== -1) continue;
1574
1575           nChannelsFromDA++;
1576
1577           AliDebug(3,Form(" channel %i noise status before noise = %i, with global status = %i",ich,(Int_t)fStatus->GetNoiseStatus(ich),(Int_t)fStatus->GetStatus(ich)));
1578           //AliDebug(2,Form( " channel %i status before noise = %i",ich,(Int_t)fStatus->GetNoiseStatus(ich)));
1579
1580           /* check whether channel has been read out during current run.
1581            * if the status is bad it means it has not been read out.
1582            * in this case skip channel in order to leave its status
1583            * unchanged */
1584
1585           if ((fStatus->GetHWStatus(ich)) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
1586                   continue;
1587
1588           nNoise++;
1589           if (htofNoise->GetBinContent(ich+1) >= noiseThr){
1590                   fStatus->SetNoiseStatus(ich,AliTOFChannelOnlineStatusArray::kTOFNoiseBad); // bad status for noise
1591                   AliDebug(3,Form( " channel %i noise status after noise = %i, with global status = %i",ich,(Int_t)fStatus->GetNoiseStatus(ich),(Int_t)fStatus->GetStatus(ich)));
1592                   nNoisyChannels++;
1593           }
1594           else {
1595                   fStatus->SetNoiseStatus(ich,AliTOFChannelOnlineStatusArray::kTOFNoiseOk); // good status for noise
1596                   AliDebug(3,Form(" channel %i noise status after noise = %i, with global status = %i",ich,(Int_t)fStatus->GetNoiseStatus(ich),(Int_t)fStatus->GetStatus(ich)));
1597                   nNotNoisyChannels++;
1598           }
1599   }
1600
1601   Log(Form(" Number of channels processed by DA during noise run, independetly from TOFFEE = %i",nChannelsFromDA));
1602   Log(Form(" Number of channels processed during noise run (that were ON according to TOFFEE) = %i",nNoise));
1603   Log(Form(" Number of noisy channels found during noise run = %i",nNoisyChannels));
1604   Log(Form(" Number of not noisy channels found during noise run = %i",nNotNoisyChannels));
1605   Log(Form(" Number of channels with matching window NULL (so skipped), but Non Zero content = %i",nMatchingWindowNullNonZero));
1606   Log(Form(" Number of channels with matching window NULL (so skipped), and Zero content = %i",nMatchingWindowNullEqualZero));
1607
1608   delete listNoise;
1609
1610   //daqFile=0;
1611
1612   //storing in OCDB
1613
1614   AliCDBMetaData metaData;
1615   metaData.SetBeamPeriod(0);
1616   metaData.SetResponsible("Chiara Zampolli");
1617   metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object after a Noise run.");
1618   AliInfo("Storing Calibration Data from Noise Run");
1619   resultNoise = Store("Calib","Status",fStatus, &metaData,0,kTRUE);
1620   if(!resultNoise){
1621     Log("Some problems occurred while storing online object resulting from Noise data processing");
1622     delete htofNoise;
1623     htofNoise = 0x0;
1624     if (fStatus){
1625             delete fStatus;
1626             fStatus = 0;
1627     }
1628     if (fMatchingWindow){
1629             delete [] fMatchingWindow;
1630             fMatchingWindow = 0;
1631     }
1632     return 14;//return error code for problems in storing Noise data
1633   }
1634
1635   if(fStoreRefData){
1636
1637     AliCDBMetaData metaDataHisto;
1638     metaDataHisto.SetBeamPeriod(0);
1639     metaDataHisto.SetResponsible("Chiara Zampolli");
1640     metaDataHisto.SetComment("This preprocessor stores the Ref data from a noise run. ");
1641     AliInfo("Storing Reference Data");
1642     resultNoiseRef = StoreReferenceData("Calib","NoiseData",htofNoise, &metaDataHisto);
1643     if (!resultNoiseRef){
1644       Log("some problems occurred::No Reference Data for noise stored");
1645       delete htofNoise;
1646       htofNoise = 0x0;
1647       if (fStatus){
1648               delete fStatus;
1649               fStatus = 0;
1650       }
1651       if (fMatchingWindow){
1652               delete [] fMatchingWindow;
1653               fMatchingWindow = 0;
1654       }
1655       return 12;//return error code for failure in storing Ref Data
1656     }
1657   }
1658
1659   delete htofNoise;
1660   htofNoise = 0x0;
1661
1662   if (fStatus){
1663     delete fStatus;
1664     fStatus = 0;
1665   }
1666
1667   if (fMatchingWindow){
1668           delete [] fMatchingWindow;
1669           fMatchingWindow = 0;
1670   }
1671
1672   return 0;
1673 }
1674 //_____________________________________________________________________________
1675
1676 UInt_t AliTOFPreprocessor::ProcessFEEData()
1677 {
1678   // Processing Pulser Run data for TOF channel status
1679   // dummy for the time being
1680
1681   Log("Processing FEE");
1682
1683   //Bool_t updateOCDB = kFALSE;
1684   AliTOFFEEReader feeReader;
1685
1686   TH1C hCurrentFEE("hCurrentFEE","histo with current FEE channel status", fNChannels, 0, fNChannels);
1687
1688   /* load current TOF FEE(dump) from DCS FXS,
1689    * setup TOFFEEdump object */
1690
1691   const char * toffeeFileName = GetFile(kDCS,"TofFeeMap","");
1692   AliInfo(Form("toffee file name = %s", toffeeFileName));
1693   if (toffeeFileName == NULL) {
1694     return 15;
1695   }
1696   AliTOFFEEDump feedump;
1697   feedump.ReadFromFile(toffeeFileName);
1698
1699   /* load current TOF FEE(light) config from DCS FXS, parse,
1700    * fill current FEE histogram and set FEE status */
1701
1702   const char * nameFile = GetFile(kDCS,"TofFeeLightMap","");
1703   AliInfo(Form("toffeeLight file name = %s",nameFile));
1704   if (nameFile == NULL) {
1705           return 15;
1706   }
1707   feeReader.LoadFEElightConfig(nameFile);
1708   Int_t parseFee = feeReader.ParseFEElightConfig();
1709   AliDebug(2,Form("%i enabled channels found in FEElight configuration",parseFee));
1710   /* load stored TOF FEE from OCDB and compare it with current FEE.
1711    * if stored FEE is different from current FEE set update flag.
1712    * if there is no stored FEE in OCDB set update flag */
1713
1714   fMatchingWindow = new Int_t[fNChannels];
1715   fLatencyWindow = new Int_t[fNChannels];
1716
1717   AliCDBEntry *cdbEntry = GetFromOCDB("Calib","Status");
1718   if (!cdbEntry) {
1719           /* no CDB entry found. set update flag */
1720           Log("     ************ WARNING ************");
1721           Log("No CDB Status entry found, creating a new one!");
1722           Log("     *********************************");
1723           fStatus = new AliTOFChannelOnlineStatusArray(fNChannels);
1724           //updateOCDB = kTRUE;
1725           fIsStatusMapChanged = kTRUE;
1726   }
1727   else {
1728           if (cdbEntry) cdbEntry->SetOwner(kFALSE);
1729           /* CDB entry OK. loop over channels */
1730           fStatus = (AliTOFChannelOnlineStatusArray*) cdbEntry->GetObject();
1731           delete cdbEntry;
1732           cdbEntry = 0x0;
1733           /* cehck whether status object has latency window data */
1734           if (!fStatus->HasLatencyWindow()) {
1735             /* create new status object and update OCDB */
1736             Log("     ************ WARNING ************");
1737             Log("CDB Status entry found but has no latency window data, creating a new one!");
1738             Log("     *********************************");
1739             delete fStatus;
1740             fStatus = new AliTOFChannelOnlineStatusArray(fNChannels);
1741             //updateOCDB = kTRUE;
1742             fIsStatusMapChanged = kTRUE;
1743           }
1744   }
1745   for (Int_t iChannel = 0; iChannel < fNChannels; iChannel++){
1746           //AliDebug(2,Form("********** channel %i",iChannel));
1747           /* compare current FEE channel status with stored one
1748            * if different set update flag and break loop */
1749           //AliDebug(2,Form( " channel %i status before FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1750           fMatchingWindow[iChannel] = feeReader.GetMatchingWindow(iChannel);
1751           fLatencyWindow[iChannel] = feeReader.GetLatencyWindow(iChannel);
1752           if (feeReader.IsChannelEnabled(iChannel)) {
1753                   hCurrentFEE.SetBinContent(iChannel + 1, 1);
1754                   if (fStatus->GetHWStatus(iChannel)!=AliTOFChannelOnlineStatusArray::kTOFHWOk){
1755                           //updateOCDB = kTRUE;
1756                           fIsStatusMapChanged = kTRUE;
1757                           fStatus->SetHWStatus(iChannel,AliTOFChannelOnlineStatusArray::kTOFHWOk);
1758                           AliDebug(3,Form( " changed into enabled: channel %i status after FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1759                   }
1760                   if (fStatus->GetLatencyWindow(iChannel)!=fLatencyWindow[iChannel]){
1761                           //updateOCDB = kTRUE;
1762                           fIsStatusMapChanged = kTRUE;
1763                           fStatus->SetLatencyWindow(iChannel,fLatencyWindow[iChannel]);
1764                           AliDebug(3,Form( " changed latency window: channel %i latency window after FEE = %i",iChannel,fStatus->GetLatencyWindow(iChannel)));
1765                   }
1766           }
1767           else {
1768                   if (fStatus->GetHWStatus(iChannel)!=AliTOFChannelOnlineStatusArray::kTOFHWBad){
1769                           //updateOCDB = kTRUE;
1770                           fIsStatusMapChanged = kTRUE;
1771                           fStatus->SetHWStatus(iChannel,AliTOFChannelOnlineStatusArray::kTOFHWBad);
1772                           AliDebug(3,Form( " changed into disabled: channel %i status after FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1773                   }
1774           }
1775   }
1776
1777
1778   /* check whether we don't have to store reference data.
1779    * in this case we return without errors. */
1780   if (fStoreRefData) {
1781     /* store reference data */
1782     AliCDBMetaData metaDataHisto;
1783     metaDataHisto.SetBeamPeriod(0);
1784     metaDataHisto.SetResponsible("Roberto Preghenella");
1785     metaDataHisto.SetComment("This preprocessor stores the FEE Ref data of the current run.");
1786     AliInfo("Storing FEE reference data");
1787     /* store FEE reference data */
1788     if (!StoreReferenceData("Calib", "FEEData", &hCurrentFEE, &metaDataHisto)) {
1789       /* failed */
1790       Log("problems while storing FEE reference data");
1791       if (fStatus){
1792         delete fStatus;
1793         fStatus = 0;
1794       }
1795       return 18; /* error return code for problems while storing FEE reference data */
1796     }
1797
1798     /* store TOF FEE dump reference data */
1799     AliCDBMetaData metaDatadump;
1800     metaDatadump.SetBeamPeriod(0);
1801     metaDatadump.SetResponsible("Roberto Preghenella");
1802     metaDatadump.SetComment("This preprocessor stores the TOF FEE dump Ref data of the current run.");
1803     AliInfo("Storing TOF FEE dump reference data");
1804     /* store FEE reference data */
1805     if (!StoreReferenceData("Calib", "FEEDump", &feedump, &metaDatadump)) {
1806       /* failed */
1807       Log("problems while storing TOF FEE dump reference data");
1808       return 18; /* error return code for problems while storing FEE reference data */
1809     }
1810   }
1811
1812   return 0;
1813
1814 }
1815
1816 //_____________________________________________________________________________
1817
1818 UInt_t AliTOFPreprocessor::Process(TMap *dcsAliasMap)
1819 {
1820   //
1821   // Main AliTOFPreprocessor method called by SHUTTLE
1822   //
1823
1824   TString runType = GetRunType();
1825   Log(Form("RunType %s",runType.Data()));
1826
1827   // processing
1828
1829   /* always process FEE data */
1830   Int_t iresultFEE = ProcessFEEData();
1831   if (iresultFEE != 0)
1832     return iresultFEE;
1833
1834   if (runType == "PULSER") {
1835     Int_t iresultPulser = ProcessPulserData();
1836     return iresultPulser;
1837   }
1838
1839   if (runType == "NOISE") { // for the time being associating noise runs with pedestal runs; proper run type to be defined
1840     Int_t iresultNoise = ProcessNoiseData();
1841     return iresultNoise;
1842   }
1843
1844   if (runType == "PHYSICS") {
1845     //    Int_t iresultDAQ = ProcessOnlineDelays();
1846     Int_t iresultDAQ = ProcessT0Fill();
1847     Int_t iresultNoiseCalib = ProcessNoiseCalibTrg();
1848     Int_t iresultReadout = ProcessReadout();
1849     Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap);
1850     Int_t iResultHVandLVdps = ProcessHVandLVdps(dcsAliasMap);
1851     return iresultDAQ+iresultNoiseCalib+iresultDCS+iResultHVandLVdps+iresultReadout;
1852   }
1853
1854   // storing
1855   return 0;
1856 }
1857
1858
1859 //_____________________________________________________________________________
1860
1861 void
1862 AliTOFPreprocessor::FillWithCosmicCalibration(AliTOFChannelOnlineArray *cal)
1863 {
1864   /*
1865    * fill with cosmic calibration
1866    */
1867
1868   Log(" Using cosmic-ray calibration.");
1869
1870   AliTOFcalibHisto calibHisto;
1871   calibHisto.SetFullCorrectionFlag(AliTOFcalibHisto::kTimeSlewingCorr, kFALSE);
1872   Log(Form(" loading calibration histograms from %s", calibHisto.GetCalibHistoFileName()));
1873   Log(Form(" loading calibration parameters from %s", calibHisto.GetCalibParFileName()));
1874   calibHisto.LoadCalibPar();
1875
1876   /* loop over channel index */
1877   for (Int_t iIndex = 0; iIndex < fNChannels; iIndex++) {
1878     cal->SetDelay(iIndex, calibHisto.GetFullCorrection(iIndex));
1879   }
1880
1881 }
1882
1883 //_____________________________________________________________________________
1884
1885 void
1886 AliTOFPreprocessor::FillWithCableLengthMap(AliTOFChannelOnlineArray *cal)
1887 {
1888   /*
1889    * fill with cosmic calibration
1890    */
1891
1892   Log(" Using cable-length map.");
1893   AliTOFRawStream tofrs;
1894   Int_t det[5], dummy, index;
1895   Float_t cableTimeShift;
1896
1897   /* temporarly disable warnings */
1898   AliLog::EType_t logLevel = (AliLog::EType_t)AliLog::GetGlobalLogLevel();
1899   AliLog::SetGlobalLogLevel(AliLog::kError);
1900
1901   /* loop over EO indeces */
1902   for (Int_t iddl = 0; iddl < 72; iddl++)
1903     for (Int_t islot = 3; islot <= 12; islot++)
1904       for (Int_t ichain = 0; ichain < 2; ichain++)
1905         for (Int_t itdc = 0; itdc < 15; itdc++)
1906           for (Int_t ichannel = 0; ichannel < 8; ichannel++) {
1907
1908             /* get DO index */
1909             tofrs.EquipmentId2VolumeId(iddl, islot, ichain, itdc, ichannel, det);
1910
1911             /* swap det[3] and det[4] indeces (needed to obtain correct channel index) */
1912             dummy = det[3];
1913             det[3] = det[4];
1914             det[4] = dummy;
1915
1916             /* check DO index */
1917             if (det[0] < 0 || det[0] > 17 ||
1918                 det[1] < 0 || det[1] > 4 ||
1919                 det[2] < 0 || det[2] > 18 ||
1920                 det[3] < 0 || det[3] > 1 ||
1921                 det[4] < 0 || det[4] > 47)
1922               continue;
1923
1924             /* get channel index */
1925             index = AliTOFGeometry::GetIndex(det);
1926             if (index < 0) continue;
1927
1928             /* get cable time shift */
1929             cableTimeShift = AliTOFCableLengthMap::GetCableTimeShift(iddl, islot, ichain, itdc);
1930
1931             /* set delay */
1932             if (index<fNChannels) {
1933               cal->SetDelay(index,cableTimeShift);  // delay in ns
1934               AliDebug(2,Form("Setting delay %f (ns) for channel %i",cableTimeShift,index));
1935             }
1936
1937           } /* loop over EO indeces */
1938
1939   /* re-enable warnings */
1940   AliLog::SetGlobalLogLevel(logLevel);
1941
1942 }
1943
1944