1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 //#include <Riostream.h>
26 #include <TObjArray.h>
27 #include <TObjString.h>
28 #include <TTimeStamp.h>
30 #include "AliCDBMetaData.h"
31 #include "AliCDBEntry.h"
33 #include "AliTOFChannelOnlineArray.h"
34 #include "AliTOFChannelOnlineStatusArray.h"
35 #include "AliTOFDataDCS.h"
36 #include "AliTOFGeometry.h"
37 #include "AliTOFPreprocessor.h"
38 #include "AliTOFFEEReader.h"
39 #include "AliTOFRawStream.h"
40 #include "AliTOFCableLengthMap.h"
41 #include "AliTOFcalibHisto.h"
44 // TOF preprocessor class.
45 // It takes data from DCS and passes them to the class AliTOFDataDCS, which
46 // processes them. The result is then written to the CDB.
47 // Analogously, it takes data form DAQ (both at Run level and inclusive -
48 // of all the runs - level, processes them, and stores both Reference Data
49 // and Online Calibration files in the CDB.
50 // Processing of Pulser/Noise Run data and of TOF FEE DCS map included also.
54 // return=1 : no DCS input data Map
55 // return=2 : no DCS input data processing
56 // return=3 : no DCS processed data was stored in Ref Data
57 // return=4 : no DAQ input for Ref Data
58 // return=5 : failed to store DAQ Ref Data
59 // return=6 : failed to retrieve DAQ data for calibration
60 // return=7 : problems in processing histos in the input DAQ file
61 // return=8 : failed to store Online Delays
62 // return=9 : failed to store Reference Data for Pulser
63 // return=10: failed to retrieve Pulser data
64 // return=11: failed to store Pulser map in OCDB
65 // return=12: failed to store Reference Data for Noise
66 // return=13: failed to retrieve Noise data
67 // return=14: failed to store Noise map in OCDB
68 // return=15: failed to retrieve FEE data from FXS
69 // return=16: failed to retrieve FEE data from OCDB
70 // return=17: failed to store FEE data in OCDB
71 // return=18: failed to store FEE reference data in OCDB
72 // return=20: failed in retrieving status variable
74 ClassImp(AliTOFPreprocessor)
76 const Int_t AliTOFPreprocessor::fgkBinRangeAve = 13; // number of bins where to calculate the mean
77 const Double_t AliTOFPreprocessor::fgkIntegralThr = 100; // min number of entries to perform computation of delay per channel
78 const Double_t AliTOFPreprocessor::fgkThrPar = 0.013; // parameter used to trigger the calculation of the delay
80 //_____________________________________________________________________________
82 AliTOFPreprocessor::AliTOFPreprocessor(AliShuttleInterface* shuttle) :
83 AliPreprocessor("TOF", shuttle),
93 AddRunType("PHYSICS");
99 //_____________________________________________________________________________
101 AliTOFPreprocessor::~AliTOFPreprocessor()
106 //______________________________________________________________________________
107 void AliTOFPreprocessor::Initialize(Int_t run, UInt_t startTime,
110 // Creates AliTOFDataDCS object
112 AliPreprocessor::Initialize(run, startTime, endTime);
114 AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s \n\tStartTime DCS Query %s \n\tEndTime DCS Query %s", run,
115 TTimeStamp(startTime).AsString(),
116 TTimeStamp(endTime).AsString(), ((TTimeStamp)GetStartTimeDCSQuery()).AsString(), ((TTimeStamp)GetEndTimeDCSQuery()).AsString()));
118 fData = new AliTOFDataDCS(fRun, fStartTime, fEndTime, GetStartTimeDCSQuery(), GetEndTimeDCSQuery());
119 fNChannels = AliTOFGeometry::NSectors()*(2*(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())*AliTOFGeometry::NpadZ()*AliTOFGeometry::NpadX();
121 //_____________________________________________________________________________
122 Bool_t AliTOFPreprocessor::ProcessDCS(){
124 // check whether DCS should be processed or not...
126 TString runType = GetRunType();
127 Log(Form("RunType %s",runType.Data()));
129 if (runType != "PHYSICS"){
135 //_____________________________________________________________________________
137 UInt_t AliTOFPreprocessor::ProcessDCSDataPoints(TMap* dcsAliasMap)
139 // Fills data into a AliTOFDataDCS object
142 Log("Processing DCS DP");
143 TH1::AddDirectory(0);
145 Bool_t resultDCSMap=kFALSE;
146 Bool_t resultDCSStore=kFALSE;
150 fData->SetFDRFlag(fFDRFlag);
153 Log("No DCS map found: TOF exiting from Shuttle");
158 return 1;// return error Code for DCS input data not found
162 // The processing of the DCS input data is forwarded to AliTOFDataDCS
163 resultDCSMap=fData->ProcessData(*dcsAliasMap);
165 Log("Some problems occurred while processing DCS data, TOF exiting from Shuttle");
170 return 2;// return error Code for processed DCS data not stored
173 AliCDBMetaData metaDataDCS;
174 metaDataDCS.SetBeamPeriod(0);
175 metaDataDCS.SetResponsible("Chiara Zampolli");
176 metaDataDCS.SetComment("This preprocessor fills an AliTOFDataDCS object.");
177 AliInfo("Storing DCS Data");
178 resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS);
179 if (!resultDCSStore){
180 Log("Some problems occurred while storing DCS data results in Reference Data, TOF exiting from Shuttle");
185 return 3;// return error Code for processed DCS data not stored
198 //_____________________________________________________________________________
200 UInt_t AliTOFPreprocessor::ProcessOnlineDelays()
202 // Processing data from DAQ for online calibration
204 Bool_t updateOCDB = kFALSE;
205 Log("Processing DAQ delays");
207 // reading configuration map
208 TString compDelays = "kFALSE";
209 Int_t deltaStartingRun = fRun;
210 Int_t startingRun = fRun-deltaStartingRun;
211 Int_t binRangeAve = fgkBinRangeAve;
212 Double_t integralThr = fgkIntegralThr;
213 Double_t thrPar = fgkThrPar;
215 AliCDBEntry *cdbEntry = GetFromOCDB("Calib","Config");
217 Log(Form("No Configuration entry found in CDB, using default values: ComputingDelays = %s, StartingRun = %i",compDelays.Data(), startingRun));
220 TMap *configMap = (TMap*)cdbEntry->GetObject();
222 Log(Form("No map found in Config entry in CDB, using default values: ComputingDelays = %s, StartingRun = %i",compDelays.Data(), startingRun));
225 TObjString *strDelays = (TObjString*)configMap->GetValue("ComputingDelays");
227 compDelays = (TString) strDelays->GetString();
230 Log(Form("No ComputingDelays value found in Map from Config entry in CDB, using default value: ComputingDelays = %s",compDelays.Data()));
232 TObjString *strRun = (TObjString*)configMap->GetValue("StartingRun");
234 TString tmpstr = strRun->GetString();
235 startingRun = tmpstr.Atoi();
236 deltaStartingRun = fRun - startingRun;
239 Log(Form("No StartingRun value found in Map from Config entry in CDB, using default value: StartingRun = %i",startingRun));
241 TObjString *strBinRangeAve = (TObjString*)configMap->GetValue("BinRangeAve");
242 if (strBinRangeAve) {
243 TString tmpstr = strBinRangeAve->GetString();
244 binRangeAve = tmpstr.Atoi();
247 Log(Form("No BinRangeAve value found in Map from Config entry in CDB, using default value: BinRangeAve = %i",binRangeAve));
249 TObjString *strIntegralThr = (TObjString*)configMap->GetValue("IntegralThr");
250 if (strIntegralThr) {
251 TString tmpstr = strIntegralThr->GetString();
252 integralThr = tmpstr.Atof();
255 Log(Form("No IntegralThr value found in Map from Config entry in CDB, using default value: IntegralThr = %i",integralThr));
257 TObjString *strThrPar = (TObjString*)configMap->GetValue("ThrPar");
259 TString tmpstr = strThrPar->GetString();
260 thrPar = tmpstr.Atof();
263 Log(Form("No ThrPar value found in Map from Config entry in CDB, using default value: ThrPar = %i",thrPar));
267 if (compDelays == "kTRUE") fFDRFlag = kFALSE;
268 else fFDRFlag = kTRUE;
273 Log(Form("ComputingDelays = %s, StartingRun = %i",compDelays.Data(),startingRun));
275 /* init array with current calibration, if any */
276 fCal = new AliTOFChannelOnlineArray(fNChannels);
277 AliTOFChannelOnlineArray *curCal = NULL;
279 AliCDBEntry *cdbEntry2 = GetFromOCDB("Calib","ParOnlineDelay");
280 if (!cdbEntry2 || !cdbEntry2->GetObject()) {
281 /* no CDB entry found. set update flag */
282 Log(" ************ WARNING ************");
283 Log("No CDB ParOnlineDelay entry found, creating a new one!");
284 Log(" *********************************");
288 Log("Found previous ParOnlineDelay entry. Using it to init calibration");
289 curCal = (AliTOFChannelOnlineArray *)cdbEntry2->GetObject();
290 for (Int_t i = 0; i < fNChannels; i++)
291 fCal->SetDelay(i, curCal->GetDelay(i));
295 TH1::AddDirectory(0);
297 Bool_t resultDAQRef=kFALSE;
298 Bool_t resultTOFPP=kFALSE;
305 //retrieving data at Run level
306 TList* list = GetFileSources(kDAQ, "RUNLevel");
307 if (list !=0x0 && list->GetEntries()!=0)
309 AliInfo("The following sources produced files with the id RUNLevel");
311 for (Int_t jj=0;jj<list->GetEntries();jj++){
312 TObjString * str = dynamic_cast<TObjString*> (list->At(jj));
313 AliInfo(Form("found source %s", str->String().Data()));
314 // file to be stored run per run
315 TString fileNameRun = GetFile(kDAQ, "RUNLevel", str->GetName());
316 if (fileNameRun.Length()>0){
317 AliInfo(Form("Got the file %s, now we can store the Reference Data for the current Run.", fileNameRun.Data()));
318 daqFile = new TFile(fileNameRun.Data(),"READ");
319 h2 = (TH2S*) daqFile->Get("htof");
320 AliCDBMetaData metaDataHisto;
321 metaDataHisto.SetBeamPeriod(0);
322 metaDataHisto.SetResponsible("Chiara Zampolli");
323 metaDataHisto.SetComment("This preprocessor stores the array of histos object as Reference Data.");
324 AliInfo("Storing Reference Data");
325 resultDAQRef = StoreReferenceData("Calib","DAQData",h2, &metaDataHisto);
327 Log("some problems occurred::No Reference Data stored, TOF exiting from Shuttle");
332 return 5;//return error code for failure in storing Ref Data
339 Log("The input data file from DAQ (run-level) was not found, TOF exiting from Shuttle ");
343 return 4;//return error code for failure in retrieving Ref Data
349 Log("The input data file list from DAQ (run-level) was not found, TOF exiting from Shuttle ");
352 return 4;//return error code for failure in retrieving Ref Data
357 //Total files, with cumulative histos
359 TList* listTot = GetFileSources(kDAQ, "DELAYS");
360 if (listTot !=0x0 && listTot->GetEntries()!=0)
362 AliInfo("The following sources produced files with the id DELAYS");
364 for (Int_t jj=0;jj<listTot->GetEntries();jj++){
365 TObjString * str = dynamic_cast<TObjString*> (listTot->At(jj));
366 AliInfo(Form("found source %s", str->String().Data()));
368 // file with summed histos, to extract calib params
369 TString fileName = GetFile(kDAQ, "DELAYS", str->GetName());
370 if (fileName.Length()>0){
371 AliInfo(Form("Got the file %s, now we can extract some values.", fileName.Data()));
373 daqFile = new TFile(fileName.Data(),"READ");
375 h2 = (TH2S*) daqFile->Get("htoftot");
377 Log("some problems occurred:: No histo retrieved, TOF exiting from Shuttle");
382 return 7; //return error code for histograms not existing/junky
385 static const Int_t kSize=h2->GetNbinsX();
386 static const Int_t kNBins=h2->GetNbinsY();
387 static const Double_t kXBinmin=h2->GetYaxis()->GetBinLowEdge(1);
388 if (kSize != fNChannels){
389 Log(" number of bins along x different from number of pads, found only a subset of the histograms, TOF exiting from Shuttle");
395 return 7; //return error code for histograms not existing/junky
397 Int_t nNotStatistics = 0; // number of channel with not enough statistics
399 /* FDR flag set. do not compute delays, use nominal cable delays */
402 Log(" Not computing delays according to flag set in Config entry in OCDB!");
403 FillWithCosmicCalibration(fCal);
405 /* check whether the new calibration is different from the previous one */
406 if (curCal) { /* well, check also whether we have a previous calibration */
407 for (Int_t i = 0; i < fNChannels; i++) {
408 if (fCal->GetDelay(i) != curCal->GetDelay(i)) {
414 else /* otherwise update OCDB */
419 else { // computing delays if not in FDR runs
421 updateOCDB = kTRUE; /* always update OCDB when computing delays */
423 for (Int_t ich=0;ich<kSize;ich++){
424 /* check whether channel has been read out during current run.
425 * if the status is bad it means it has not been read out.
426 * in this case skip channel in order to not affect the mean */
427 if (fStatus->GetHWStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFHWBad){
428 AliDebug(2,Form(" Channel %i found bad according to FEEmap, (HW status = %i), skipping from delay computing",ich, (Int_t)fStatus->GetHWStatus(ich)));
431 AliDebug(2,Form(" Channel %i found ok according to FEEmap, starting delay computing",ich));
432 TH1S *h1 = new TH1S("h1","h1",kNBins,kXBinmin-0.5,kNBins*1.+kXBinmin-0.5);
433 for (Int_t ibin=0;ibin<kNBins;ibin++){
434 h1->SetBinContent(ibin+1,h2->GetBinContent(ich+1,ibin+1));
436 if(h1->Integral()<integralThr) {
438 Log(Form(" Not enough statistics for bin %i, skipping this channel",ich)); // printing message only if not in FDR runs
444 Float_t minContent=h1->Integral()*thrPar;
445 Int_t nbinsX = h1->GetNbinsX();
447 for (Int_t j=1; j<=nbinsX; j++){
449 h1->GetBinContent(j) +
450 h1->GetBinContent(j+1)+
451 h1->GetBinContent(j+2)+
452 h1->GetBinContent(j+3))>minContent){
458 if(!found) AliInfo(Form("WARNING!!! no start of fit found for histo # %i",ich));
459 // Now calculate the mean over the interval.
463 for(Int_t k=0;k<binRangeAve;k++){
464 mean=mean+h1->GetBinCenter(startBin+k)*h1->GetBinContent(startBin+k);
465 nent=nent+h1->GetBinContent(startBin+k);
466 sumw2=sumw2+(h1->GetBinCenter(startBin+k))*(h1->GetBinCenter(startBin+k))*(h1->GetBinContent(startBin+k));
468 mean= mean/nent; //<x>
469 sumw2=sumw2/nent; //<x^2>
471 rmsmean = TMath::Sqrt((sumw2-mean*mean)/nent);
472 if (ich<fNChannels) {
473 Float_t delay = mean*AliTOFGeometry::TdcBinWidth()*1.E-3; // delay in ns
474 fCal->SetDelay(ich,delay); // delay in ns
475 AliDebug(2,Form("Setting delay %f (ns) for channel %i",delay,ich));
481 if (nNotStatistics!=0) Log(Form("Too little statistics for %d channels!",nNotStatistics));
488 Log("The Cumulative data file from DAQ does not exist, TOF exiting from Shuttle");
492 return 6;//return error code for problems in retrieving DAQ data
498 Log("Problem: no list for Cumulative data file from DAQ was found, TOF exiting from Shuttle");
501 return 6; //return error code for problems in retrieving DAQ data
504 /* check whether we don't need to update OCDB.
505 * in this case we can return without errors and
506 * the current FEE is stored in the fStatus object. */
508 AliInfo("update OCDB flag not set. Do not overwrite stored file.");
509 return 0; /* return ok */
513 AliCDBMetaData metaData;
514 metaData.SetBeamPeriod(0);
515 metaData.SetResponsible("Chiara Zampolli");
516 metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineArray object for online calibration - delays.");
517 AliInfo("Storing Calibration Data");
518 resultTOFPP = Store("Calib","ParOnlineDelay",fCal, &metaData,deltaStartingRun,kTRUE);
520 Log("Some problems occurred while storing online object resulting from DAQ data processing");
523 return 8;//return error code for problems in storing DAQ data
533 //_____________________________________________________________________________
535 UInt_t AliTOFPreprocessor::ProcessPulserData()
537 // Processing Pulser Run data for TOF channel status
539 Log("Processing Pulser");
542 AliError("No valid fStatus found, some errors must have occurred!!");
546 TH1::AddDirectory(0);
548 Bool_t resultPulserRef=kFALSE;
549 Bool_t resultPulser=kFALSE;
551 static const Int_t kSize = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
552 TH1S * htofPulser = new TH1S("hTOFpulser","histo with signals on TOF during pulser", kSize,-0.5,kSize-0.5);
553 for (Int_t ibin =1;ibin<=kSize;ibin++){
554 htofPulser->SetBinContent(ibin,-1);
562 //retrieving Pulser data
563 TList* listPulser = GetFileSources(kDAQ, "PULSER");
564 if (listPulser !=0x0 && listPulser->GetEntries()!=0)
566 AliInfo("The following sources produced files with the id PULSER");
569 for (Int_t jj=0;jj<listPulser->GetEntries();jj++){
570 Int_t nPulserSource = 0;
571 TObjString * str = dynamic_cast<TObjString*> (listPulser->At(jj));
572 AliInfo(Form("found source %s", str->String().Data()));
573 // file to be stored run per run
574 TString fileNamePulser = GetFile(kDAQ, "PULSER", str->GetName());
575 if (fileNamePulser.Length()>0){
576 // storing refernce data
577 AliInfo(Form("Got the file %s, now we can process pulser data.", fileNamePulser.Data()));
578 daqFile = new TFile(fileNamePulser.Data(),"READ");
579 h1 = (TH1S*) daqFile->Get("hTOFpulser");
580 for (Int_t ibin=0;ibin<kSize;ibin++){
581 if ((h1->GetBinContent(ibin+1))!=-1){
582 if ((htofPulser->GetBinContent(ibin+1))==-1){
583 htofPulser->SetBinContent(ibin+1,h1->GetBinContent(ibin+1));
586 Log(Form("Something strange occurred during Pulser run, channel %i already read by another LDC, please check!",ibin));
594 Int_t nreadNotEmpty=0;
595 for (Int_t ientry=1;ientry<=h1->GetNbinsX();ientry++){
597 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)));
598 /* check whether channel has been read out during current run.
599 * if the status is bad it means it has not been read out.
600 * in this case skip channel in order to not affect the mean */
601 if (fStatus->GetHWStatus(ientry-1) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
605 if (h1->GetBinContent(ientry)==-1) continue;
607 if (h1->GetBinContent(ientry)>0) {
609 AliDebug(2,Form(" channel %i is ok with entry = %f; so far %i channels added ",ientry-1,h1->GetBinContent(ientry),nreadNotEmpty));
611 mean+=h1->GetBinContent(ientry);
617 AliDebug(2,Form(" nread = %i , nreadNotEmpty = %i, mean = %f",nread,nreadNotEmpty,mean));
618 for (Int_t ich =0;ich<fNChannels;ich++){
619 if (h1->GetBinContent(ich+1)==-1) continue;
620 AliDebug(3,Form(" channel %i pulser status before pulser = %i",ich,(Int_t)fStatus->GetPulserStatus(ich)));
622 /* check whether channel has been read out during current run.
623 * if the status is bad it means it has not been read out.
624 * in this case skip channel in order to leave its status
626 if (fStatus->GetHWStatus(ich) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
629 if (h1->GetBinContent(ich+1)<0.05*mean){
630 fStatus->SetPulserStatus(ich,AliTOFChannelOnlineStatusArray::kTOFPulserBad); // bad status for pulser
631 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)));
634 fStatus->SetPulserStatus(ich,AliTOFChannelOnlineStatusArray::kTOFPulserOk); // good status for pulser
635 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)));
640 Log("No channels read!! No action taken, keeping old status");
649 Log("The input data file from DAQ (pulser) was not found, TOF exiting from Shuttle ");
657 return 10;//return error code for failure in retrieving Ref Data
659 AliDebug(2,Form(" Number of channels processed during pulser run from source %i = %i",jj, nPulserSource));
661 AliDebug(2,Form(" Number of channels processed during pulser run = %i",nPulser));
666 Log("The input data file list from DAQ (pulser) was not found, TOF exiting from Shuttle ");
673 return 10;//return error code for failure in retrieving Ref Data
678 AliCDBMetaData metaData;
679 metaData.SetBeamPeriod(0);
680 metaData.SetResponsible("Chiara Zampolli");
681 metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object after a Pulser run.");
682 AliInfo("Storing Calibration Data from Pulser Run");
683 resultPulser = Store("Calib","Status",fStatus, &metaData,0,kTRUE);
685 Log("Some problems occurred while storing online object resulting from Pulser data processing");
692 return 11;//return error code for problems in storing Pulser data
697 AliCDBMetaData metaDataHisto;
698 metaDataHisto.SetBeamPeriod(0);
699 metaDataHisto.SetResponsible("Chiara Zampolli");
701 sprintf(comment,"This preprocessor stores the Ref data from a pulser run.");
702 metaDataHisto.SetComment(comment);
703 AliInfo("Storing Reference Data");
704 resultPulserRef = StoreReferenceData("Calib","PulserData",htofPulser, &metaDataHisto);
705 if (!resultPulserRef){
706 Log("some problems occurred::No Reference Data for pulser stored, TOF exiting from Shuttle");
713 return 9;//return error code for failure in storing Ref Data
729 //_____________________________________________________________________________
731 UInt_t AliTOFPreprocessor::ProcessNoiseData()
734 // Processing Noise Run data for TOF channel status
736 Log("Processing Noise");
739 AliError("No valid fStatus found, some errors must have occurred!!");
743 Float_t noiseThr = 1; // setting default threshold for noise to 1 Hz
744 // reading config map
745 AliCDBEntry *cdbEntry = GetFromOCDB("Calib","ConfigNoise");
747 Log(Form("No Configuration entry found in CDB, using default values: NoiseThr = %d",noiseThr));
750 TMap *configMap = (TMap*)cdbEntry->GetObject();
752 Log(Form("No map found in Config entry in CDB, using default values: NoiseThr = %d", noiseThr));
755 TObjString *strNoiseThr = (TObjString*)configMap->GetValue("NoiseThr");
757 TString tmpstr = strNoiseThr->GetString();
758 noiseThr = tmpstr.Atoi();
761 Log(Form("No NoiseThr value found in Map from ConfigNoise entry in CDB, using default value: NoiseThr = %i",noiseThr));
769 TH1::AddDirectory(0);
771 Bool_t resultNoiseRef=kFALSE;
772 Bool_t resultNoise=kFALSE;
774 static const Int_t kSize = AliTOFGeometry::NPadXSector()*AliTOFGeometry::NSectors();
775 TH1F * htofNoise = new TH1F("hTOFnoise","histo with signals on TOF during noise", kSize,-0.5,kSize-0.5);
776 for (Int_t ibin =1;ibin<=kSize;ibin++){
777 htofNoise->SetBinContent(ibin,-1);
787 Int_t nNoisyChannels = 0;
788 Int_t nNotNoisyChannels = 0;
789 Int_t nChannelsFromDA = 0;
790 Int_t nMatchingWindowNullNonZero = 0;
791 Int_t nMatchingWindowNullEqualZero = 0;
793 // retrieving Noise data
794 TList* listNoise = GetFileSources(kDAQ, "NOISE");
795 if (listNoise !=0x0 && listNoise->GetEntries()!=0)
797 AliInfo("The following sources produced files with the id NOISE");
799 for (Int_t jj=0;jj<listNoise->GetEntries();jj++){
800 Int_t nNoiseSource = 0;
801 TObjString * str = dynamic_cast<TObjString*> (listNoise->At(jj));
802 AliInfo(Form("found source %s", str->String().Data()));
803 // file to be stored run per run
804 TString fileNameNoise = GetFile(kDAQ, "NOISE", str->GetName());
805 if (fileNameNoise.Length()>0){
806 // storing reference data
807 AliInfo(Form("Got the file %s, now we can process noise data.", fileNameNoise.Data()));
808 daqFile = new TFile(fileNameNoise.Data(),"READ");
809 h1 = (TH1F*) daqFile->Get("hTOFnoise");
810 for (Int_t ibin=0;ibin<kSize;ibin++){
811 if ((h1->GetBinContent(ibin+1))!=-1){
813 // checking the matching window for current channel
814 if (fMatchingWindow[ibin] == 0){
815 Log(Form("Matching window for channel %i null, but the channel was read by the LDC! skipping channel, BUT Please check!",ibin));
816 if ((h1->GetBinContent(ibin+1))!=0) nMatchingWindowNullNonZero++;
817 if ((h1->GetBinContent(ibin+1))==0) nMatchingWindowNullEqualZero++;
820 if ((htofNoise->GetBinContent(ibin+1))==-1){
821 htofNoise->SetBinContent(ibin+1,h1->GetBinContent(ibin+1)/(fMatchingWindow[ibin]*1.E-9));
822 if ((h1->GetBinContent(ibin+1))!= 0) AliDebug(2,Form("Channel = %i, Matching window = %i, Content = %f", ibin, fMatchingWindow[ibin], htofNoise->GetBinContent(ibin+1)));
825 Log(Form("Something strange occurred during Noise run, channel %i already read by another LDC, please check!",ibin));
830 Log(Form(" Number of channels processed during noise run from source %i = %i",jj, nNoiseSource));
839 Log("The input data file from DAQ (noise) was not found, TOF exiting from Shuttle ");
848 if (fMatchingWindow){
849 delete [] fMatchingWindow;
852 return 13;//return error code for failure in retrieving Ref Data
858 Log("The input data file list from DAQ (noise) was not found, TOF exiting from Shuttle ");
865 if (fMatchingWindow){
866 delete [] fMatchingWindow;
869 return 13;//return error code for failure in retrieving Ref Data
872 // elaborating infos to set NOISE status
873 for (Int_t ich =0;ich<fNChannels;ich++){
874 if (htofNoise->GetBinContent(ich+1)== -1) continue;
878 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)));
879 //AliDebug(2,Form( " channel %i status before noise = %i",ich,(Int_t)fStatus->GetNoiseStatus(ich)));
881 /* check whether channel has been read out during current run.
882 * if the status is bad it means it has not been read out.
883 * in this case skip channel in order to leave its status
886 if ((fStatus->GetHWStatus(ich)) == AliTOFChannelOnlineStatusArray::kTOFHWBad)
890 if (htofNoise->GetBinContent(ich+1) >= noiseThr){
891 fStatus->SetNoiseStatus(ich,AliTOFChannelOnlineStatusArray::kTOFNoiseBad); // bad status for noise
892 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)));
896 fStatus->SetNoiseStatus(ich,AliTOFChannelOnlineStatusArray::kTOFNoiseOk); // good status for noise
897 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)));
902 Log(Form(" Number of channels processed by DA during noise run, independetly from TOFFEE = %i",nChannelsFromDA));
903 Log(Form(" Number of channels processed during noise run (that were ON according to TOFFEE) = %i",nNoise));
904 Log(Form(" Number of noisy channels found during noise run = %i",nNoisyChannels));
905 Log(Form(" Number of not noisy channels found during noise run = %i",nNotNoisyChannels));
906 Log(Form(" Number of channels with matching window NULL (so skipped), but Non Zero content = %i",nMatchingWindowNullNonZero));
907 Log(Form(" Number of channels with matching window NULL (so skipped), and Zero content = %i",nMatchingWindowNullEqualZero));
915 AliCDBMetaData metaData;
916 metaData.SetBeamPeriod(0);
917 metaData.SetResponsible("Chiara Zampolli");
918 metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object after a Noise run.");
919 AliInfo("Storing Calibration Data from Noise Run");
920 resultNoise = Store("Calib","Status",fStatus, &metaData,0,kTRUE);
922 Log("Some problems occurred while storing online object resulting from Noise data processing");
929 if (fMatchingWindow){
930 delete [] fMatchingWindow;
933 return 14;//return error code for problems in storing Noise data
938 AliCDBMetaData metaDataHisto;
939 metaDataHisto.SetBeamPeriod(0);
940 metaDataHisto.SetResponsible("Chiara Zampolli");
942 sprintf(comment,"This preprocessor stores the Ref data from a noise run. ");
943 metaDataHisto.SetComment(comment);
944 AliInfo("Storing Reference Data");
945 resultNoiseRef = StoreReferenceData("Calib","NoiseData",htofNoise, &metaDataHisto);
946 if (!resultNoiseRef){
947 Log("some problems occurred::No Reference Data for noise stored");
954 if (fMatchingWindow){
955 delete [] fMatchingWindow;
958 return 12;//return error code for failure in storing Ref Data
970 if (fMatchingWindow){
971 delete [] fMatchingWindow;
977 //_____________________________________________________________________________
979 UInt_t AliTOFPreprocessor::ProcessFEEData()
981 // Processing Pulser Run data for TOF channel status
982 // dummy for the time being
984 Log("Processing FEE");
986 Bool_t updateOCDB = kFALSE;
987 AliTOFFEEReader feeReader;
989 TH1C hCurrentFEE("hCurrentFEE","histo with current FEE channel status", fNChannels, 0, fNChannels);
991 /* load current TOF FEE config from DCS FXS, parse,
992 * fill current FEE histogram and set FEE status */
994 const char * nameFile = GetFile(kDCS,"TofFeeLightMap","");
995 AliInfo(Form("nameFile = %s",nameFile));
996 if (nameFile == NULL) {
999 feeReader.LoadFEElightConfig(nameFile);
1000 Int_t parseFee = feeReader.ParseFEElightConfig();
1001 AliDebug(2,Form("%i enabled channels found in FEElight configuration",parseFee));
1002 /* load stored TOF FEE from OCDB and compare it with current FEE.
1003 * if stored FEE is different from current FEE set update flag.
1004 * if there is no stored FEE in OCDB set update flag */
1006 fMatchingWindow = new Int_t[fNChannels];
1008 AliCDBEntry *cdbEntry = GetFromOCDB("Calib","Status");
1010 /* no CDB entry found. set update flag */
1011 Log(" ************ WARNING ************");
1012 Log("No CDB Status entry found, creating a new one!");
1013 Log(" *********************************");
1014 fStatus = new AliTOFChannelOnlineStatusArray(fNChannels);
1018 if (cdbEntry) cdbEntry->SetOwner(kFALSE);
1019 /* CDB entry OK. loop over channels */
1020 fStatus = (AliTOFChannelOnlineStatusArray*) cdbEntry->GetObject();
1024 for (Int_t iChannel = 0; iChannel < fNChannels; iChannel++){
1025 //AliDebug(2,Form("********** channel %i",iChannel));
1026 /* compare current FEE channel status with stored one
1027 * if different set update flag and break loop */
1028 //AliDebug(2,Form( " channel %i status before FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1029 fMatchingWindow[iChannel] = feeReader.GetMatchingWindow(iChannel);
1030 if (feeReader.IsChannelEnabled(iChannel)) {
1031 hCurrentFEE.SetBinContent(iChannel + 1, 1);
1032 if (fStatus->GetHWStatus(iChannel)!=AliTOFChannelOnlineStatusArray::kTOFHWOk){
1034 fStatus->SetHWStatus(iChannel,AliTOFChannelOnlineStatusArray::kTOFHWOk);
1035 AliDebug(3,Form( " changed into enabled: channel %i status after FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1039 if (fStatus->GetHWStatus(iChannel)!=AliTOFChannelOnlineStatusArray::kTOFHWBad){
1041 fStatus->SetHWStatus(iChannel,AliTOFChannelOnlineStatusArray::kTOFHWBad);
1042 AliDebug(3,Form( " changed into disabled: channel %i status after FEE = %i",iChannel,(Int_t)fStatus->GetHWStatus(iChannel)));
1048 /* check whether we don't have to store reference data.
1049 * in this case we return without errors. */
1050 if (fStoreRefData) {
1051 /* store reference data */
1052 AliCDBMetaData metaDataHisto;
1053 metaDataHisto.SetBeamPeriod(0);
1054 metaDataHisto.SetResponsible("Roberto Preghenella");
1055 metaDataHisto.SetComment("This preprocessor stores the FEE Ref data of the current run.");
1056 AliInfo("Storing FEE reference data");
1057 /* store FEE reference data */
1058 if (!StoreReferenceData("Calib", "FEEData", &hCurrentFEE, &metaDataHisto)) {
1060 Log("problems while storing FEE reference data");
1065 return 18; /* error return code for problems while storing FEE reference data */
1069 /* check whether we don't need to update OCDB.
1070 * in this case we can return without errors and
1071 * the current FEE is stored in the fStatus object. */
1073 AliInfo("TOF FEE config has not changed. Do not overwrite stored file.");
1074 return 0; /* return ok */
1077 TString runType = GetRunType();
1078 if (runType != "PHYSICS") {
1079 AliInfo(Form("Run Type = %s, waiting to store status map",GetRunType()));
1080 return 0; /* return ok */
1083 /* update the OCDB with the current FEE since even
1084 * a little difference has been detected. */
1086 AliCDBMetaData metaData;
1087 metaData.SetBeamPeriod(0);
1088 metaData.SetResponsible("Roberto Preghenella");
1089 metaData.SetComment("This preprocessor fills an AliTOFChannelOnlineStatusArray object from FEE data.");
1090 AliInfo("Storing Status data from current run after FEE parsing");
1091 /* store FEE data */
1092 if (!Store("Calib", "Status", fStatus, &metaData, 0, kTRUE)) {
1094 Log("problems while storing FEE data object");
1099 return 17; /* return error code for problems while storing FEE data */
1102 /* everything fine. return */
1113 //_____________________________________________________________________________
1115 UInt_t AliTOFPreprocessor::Process(TMap* dcsAliasMap)
1121 TString runType = GetRunType();
1122 Log(Form("RunType %s",runType.Data()));
1126 /* always process FEE data */
1127 Int_t iresultFEE = ProcessFEEData();
1128 if (iresultFEE != 0)
1131 if (runType == "PULSER") {
1132 Int_t iresultPulser = ProcessPulserData();
1133 return iresultPulser;
1136 if (runType == "NOISE") { // for the time being associating noise runs with pedestal runs; proper run type to be defined
1137 Int_t iresultNoise = ProcessNoiseData();
1138 return iresultNoise;
1141 if (runType == "PHYSICS") {
1142 Int_t iresultDAQ = ProcessOnlineDelays();
1143 if (iresultDAQ != 0) {
1147 Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap);
1157 //_____________________________________________________________________________
1160 AliTOFPreprocessor::FillWithCosmicCalibration(AliTOFChannelOnlineArray *cal)
1163 * fill with cosmic calibration
1166 Log(" Using cosmic-ray calibration.");
1168 AliTOFcalibHisto calibHisto;
1169 Log(Form(" loading calibration histograms from %s", calibHisto.GetCalibHistoFileName()));
1170 Log(Form(" loading calibration parameters from %s", calibHisto.GetCalibParFileName()));
1171 calibHisto.LoadCalibPar();
1173 /* loop over channel index */
1174 for (Int_t iIndex = 0; iIndex < fNChannels; iIndex++) {
1175 cal->SetDelay(iIndex, calibHisto.GetFullCorrection(iIndex));
1180 //_____________________________________________________________________________
1183 AliTOFPreprocessor::FillWithCableLengthMap(AliTOFChannelOnlineArray *cal)
1186 * fill with cosmic calibration
1189 Log(" Using cable-length map.");
1190 AliTOFRawStream tofrs;
1191 Int_t det[5], dummy, index;
1192 Float_t cableTimeShift;
1194 /* temporarly disable warnings */
1195 AliLog::EType_t logLevel = (AliLog::EType_t)AliLog::GetGlobalLogLevel();
1196 AliLog::SetGlobalLogLevel(AliLog::kError);
1198 /* loop over EO indeces */
1199 for (Int_t iddl = 0; iddl < 72; iddl++)
1200 for (Int_t islot = 3; islot <= 12; islot++)
1201 for (Int_t ichain = 0; ichain < 2; ichain++)
1202 for (Int_t itdc = 0; itdc < 15; itdc++)
1203 for (Int_t ichannel = 0; ichannel < 8; ichannel++) {
1206 tofrs.EquipmentId2VolumeId(iddl, islot, ichain, itdc, ichannel, det);
1208 /* swap det[3] and det[4] indeces (needed to obtain correct channel index) */
1213 /* check DO index */
1214 if (det[0] < 0 || det[0] > 17 ||
1215 det[1] < 0 || det[1] > 4 ||
1216 det[2] < 0 || det[2] > 18 ||
1217 det[3] < 0 || det[3] > 1 ||
1218 det[4] < 0 || det[4] > 47)
1221 /* get channel index */
1222 index = AliTOFGeometry::GetIndex(det);
1224 /* get cable time shift */
1225 cableTimeShift = AliTOFCableLengthMap::GetCableTimeShift(iddl, islot, ichain, itdc);
1228 if (index<fNChannels) {
1229 cal->SetDelay(index,cableTimeShift); // delay in ns
1230 AliDebug(2,Form("Setting delay %f (ns) for channel %i",cableTimeShift,index));
1233 } /* loop over EO indeces */
1235 /* re-enable warnings */
1236 AliLog::SetGlobalLogLevel(logLevel);