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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
19 // EMCAL Preprocessor class. It runs by Shuttle at the end of the run,
20 // calculates stuff to be posted in OCDB
22 // Author: Boris Polichtchouk, 4 October 2006
23 // Adapted for EMCAL by Gustavo Conesa Balbastre, October 2006
24 // Updated by David Silvermyr May 2008, based on TPC code
25 ///////////////////////////////////////////////////////////////////////////////
31 #include "TParameter.h"
33 #include <TTimeStamp.h>
36 #include "AliShuttleInterface.h"
37 #include "AliEMCALPreprocessor.h"
39 #include "AliCDBMetaData.h"
40 #include "AliCaloCalibPedestal.h"
41 #include "AliCaloCalibSignal.h"
42 #include "AliEMCALSensorTempArray.h"
44 const Int_t kValCutTemp = 100; // discard temperatures > 100 degrees
45 const Int_t kDiffCutTemp = 5; // discard temperature differences > 5 degrees
46 const TString kPedestalRunType = "PEDESTAL"; // pedestal run identifier
47 const TString kPhysicsRunType = "PHYSICS"; // physics run identifier
48 const TString kStandAloneRunType = "STANDALONE"; // standalone run identifier
49 const TString kAmandaTemp = "PT_%02d.Temperature"; // Amanda string for temperature entries
50 //const Double_t kFitFraction = 0.7; // Fraction of DCS sensor fits required
51 const Double_t kFitFraction = -1.0; // Don't require minimum number of fits during commissioning
53 const TString kMetaResponsible = "David Silvermyr";
54 //legacy comments and return codes from TPC
55 const TString kMetaComment = "Preprocessor AliEMCAL data base entries.";
56 const int kReturnCodeNoInfo = 9;
57 const int kReturnCodeNoObject = 2;
58 const int kReturnCodeNoEntries = 1;
60 ClassImp(AliEMCALPreprocessor)
62 //_______________________________________________________________________________________
63 AliEMCALPreprocessor::AliEMCALPreprocessor() :
64 AliPreprocessor("EMC",0),
72 //_______________________________________________________________________________________
73 AliEMCALPreprocessor::AliEMCALPreprocessor(AliShuttleInterface* shuttle):
74 AliPreprocessor("EMC",shuttle),
79 // Constructor AddRunType(kPedestalRunType);
81 // define run types to be processed
82 AddRunType(kPedestalRunType);
83 AddRunType(kPhysicsRunType);
84 AddRunType(kStandAloneRunType);
87 //______________________________________________________________________________________________
88 AliEMCALPreprocessor::AliEMCALPreprocessor(const AliEMCALPreprocessor& ) :
89 AliPreprocessor("EMCAL",0),
90 fConfEnv(0), fTemp(0), fConfigOK(kTRUE)
92 Fatal("AliEMCALPreprocessor", "copy constructor not implemented");
95 // assignment operator; use copy ctor to make life easy.
96 //______________________________________________________________________________________________
97 AliEMCALPreprocessor& AliEMCALPreprocessor::operator = (const AliEMCALPreprocessor &source )
99 // assignment operator; use copy ctor
100 if (&source == this) return *this;
102 new (this) AliEMCALPreprocessor(source);
106 //____________________________________________________________________________
107 AliEMCALPreprocessor::~AliEMCALPreprocessor()
110 if (fTemp) delete fTemp;
113 //______________________________________________________________________________________________
114 void AliEMCALPreprocessor::Initialize(Int_t run, UInt_t startTime,
117 // Creates AliTestDataDCS object -- start maps half an hour beforre actual run start
118 UInt_t startTimeLocal = startTime-1800;
119 AliPreprocessor::Initialize(run, startTimeLocal, endTime);
121 AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
122 TTimeStamp((time_t)startTime,0).AsString(),
123 TTimeStamp((time_t)endTime,0).AsString()));
125 // Preprocessor configuration
126 AliCDBEntry* entry = GetFromOCDB("Config", "Preprocessor");
127 if (entry) fConfEnv = (TEnv*) entry->GetObject();
129 Log("AliEMCALPreprocessor: Preprocessor Config OCDB entry missing.\n");
134 // Temperature sensors
137 TString tempConf = fConfEnv->GetValue("Temperature","ON");
139 if (tempConf != "OFF" ) {
140 entry = GetFromOCDB("Config", "Temperature");
141 if (entry) confTree = (TTree*) entry->GetObject();
143 Log("AliEMCALPreprocessor: Temperature Config OCDB entry missing.\n");
147 fTemp = new AliEMCALSensorTempArray(startTimeLocal, fEndTime, confTree, kAmandaTemp);
148 fTemp->SetValCut(kValCutTemp);
149 fTemp->SetDiffCut(kDiffCutTemp);
155 //______________________________________________________________________________________________
156 UInt_t AliEMCALPreprocessor::Process(TMap* dcsAliasMap)
158 // Fills data into EMCAL calibrations objects
159 // Amanda servers provide information directly through dcsAliasMap
161 if (!fConfigOK) return kReturnCodeNoInfo;
163 TObjArray *resultArray = new TObjArray();
164 TString errorHandling = fConfEnv->GetValue("ErrorHandling","ON");
165 errorHandling.ToUpper();
169 if (errorHandling == "OFF" ) {
170 if (!dcsAliasMap) dcsResult = kReturnCodeNoEntries;
171 if (dcsAliasMap->GetEntries() == 0 ) dcsResult = kReturnCodeNoEntries;
172 status = new TParameter<int>("dcsResult",dcsResult);
173 resultArray->Add(status);
176 if (!dcsAliasMap) return kReturnCodeNoInfo;
177 if (dcsAliasMap->GetEntries() == 0 ) return kReturnCodeNoInfo;
180 TString runType = GetRunType();
182 // Temperature sensors are processed by AliEMCALCalTemp
183 TString tempConf = fConfEnv->GetValue("Temperature","ON");
185 if (tempConf != "OFF" ) {
186 UInt_t tempResult = MapTemperature(dcsAliasMap);
188 status = new TParameter<int>("tempResult",tempResult);
189 resultArray->Add(status);
192 // Other calibration information will be retrieved through FXS files
194 // TList* fileSourcesDAQ = GetFile(AliShuttleInterface::kDAQ, "pedestals");
195 // const char* fileNamePed = GetFile(AliShuttleInterface::kDAQ, "pedestals", "LDC1");
197 // TList* fileSourcesHLT = GetFile(AliShuttleInterface::kHLT, "calib");
198 // const char* fileNameHLT = GetFile(AliShuttleInterface::kHLT, "calib", "LDC1");
202 if(runType == kPedestalRunType) {
203 Int_t numSources = 1;
204 Int_t pedestalSource[2] = {AliShuttleInterface::kDAQ, AliShuttleInterface::kHLT} ;
205 TString source = fConfEnv->GetValue("Pedestal","DAQ");
207 if (source != "OFF" ) {
208 if ( source == "HLT") pedestalSource[0] = AliShuttleInterface::kHLT;
209 if (!GetHLTStatus()) pedestalSource[0] = AliShuttleInterface::kDAQ;
210 if (source == "HLTDAQ" ) {
212 pedestalSource[0] = AliShuttleInterface::kHLT;
213 pedestalSource[1] = AliShuttleInterface::kDAQ;
215 if (source == "DAQHLT" ) numSources=2;
216 UInt_t pedestalResult=0;
217 for (Int_t i=0; i<numSources; i++ ) {
218 pedestalResult = ExtractPedestals(pedestalSource[i]);
219 if ( pedestalResult == 0 ) break;
221 result += pedestalResult;
222 status = new TParameter<int>("pedestalResult",pedestalResult);
223 resultArray->Add(status);
227 // SIGNAL/LED ENTRIES:
229 if( runType == kPhysicsRunType || runType == kStandAloneRunType ) {
230 Int_t numSources = 1;
231 Int_t signalSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ;
232 TString source = fConfEnv->GetValue("Signal","DAQ");
234 if ( source != "OFF") {
235 if ( source == "HLT") signalSource[0] = AliShuttleInterface::kHLT;
236 if (!GetHLTStatus()) signalSource[0] = AliShuttleInterface::kDAQ;
237 if (source == "HLTDAQ" ) {
239 signalSource[0] = AliShuttleInterface::kHLT;
240 signalSource[1] = AliShuttleInterface::kDAQ;
242 if (source == "DAQHLT" ) numSources=2;
243 UInt_t signalResult=0;
244 for (Int_t i=0; i<numSources; i++ ) {
245 signalResult = ExtractSignal(signalSource[i]);
246 if ( signalResult == 0 ) break;
248 result += signalResult;
249 status = new TParameter<int>("signalResult",signalResult);
250 resultArray->Add(status);
255 // overall status at the end
256 if (errorHandling == "OFF" ) {
257 AliCDBMetaData metaData;
258 metaData.SetBeamPeriod(0);
259 metaData.SetResponsible(kMetaResponsible);
260 metaData.SetComment("Preprocessor AliEMCAL status.");
261 Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE);
262 resultArray->Delete();
270 //______________________________________________________________________________________________
271 UInt_t AliEMCALPreprocessor::MapTemperature(TMap* dcsAliasMap)
273 // extract DCS temperature maps. Perform fits to save space
276 TMap *map = fTemp->ExtractDCS(dcsAliasMap);
278 fTemp->MakeSplineFit(map);
279 Double_t fitFraction = 1.0*fTemp->NumFits()/fTemp->NumSensors();
280 if (fitFraction > kFitFraction ) {
281 AliInfo(Form("Temperature values extracted, fits performed.\n"));
284 Log ("Too few temperature maps fitted. \n");
285 result = kReturnCodeNoInfo;
289 Log("No temperature map extracted. \n");
290 result = kReturnCodeNoInfo;
293 // Now store the final CDB file
295 if ( result == 0 ) { // some info was found
296 AliCDBMetaData metaData;
297 metaData.SetBeamPeriod(0);
298 metaData.SetResponsible(kMetaResponsible);
299 metaData.SetComment(kMetaComment);
301 Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE);
302 if ( !storeOK ) result=1;
308 //______________________________________________________________________________________________
309 UInt_t AliEMCALPreprocessor::ExtractPedestals(Int_t sourceFXS)
313 // Read pedestal file from file exchange server
314 // Keep original entry from OCDB in case no new pedestals are available
316 AliCaloCalibPedestal *calibPed=0;
317 AliCDBEntry* entry = GetFromOCDB("Calib", "Pedestals");
318 if (entry) calibPed = (AliCaloCalibPedestal*)entry->GetObject();
319 if ( calibPed==NULL ) {
320 Log("AliEMCALPreprocessor: No previous EMCAL pedestal entry available.\n");
321 calibPed = new AliCaloCalibPedestal(AliCaloCalibPedestal::kEmCal);
324 TList* list = GetFileSources(sourceFXS,"pedestals");
325 if (list && list->GetEntries()>0) {
327 calibPed->Reset(); // let's make a fresh start before possibly adding stuff below
329 // loop through all files from LDCs
333 while (list->At(index)!=NULL) {
334 TObjString* fileNameEntry = (TObjString*) list->At(index);
335 if (fileNameEntry!=NULL) {
336 TString fileName = GetFile(sourceFXS, "pedestals",
337 fileNameEntry->GetString().Data());
338 TFile *f = TFile::Open(fileName);
340 Log ("Error opening pedestal file.");
341 result = kReturnCodeNoObject;
344 AliCaloCalibPedestal *calPed;
345 f->GetObject("emcCalibPedestal",calPed);
347 Log ("No pedestal calibration object in file.");
348 result = kReturnCodeNoObject;
351 if ( calPed->GetNEvents()>0 && calPed->GetNChanFills()>0 ) {
352 // add info for the modules available in the present file
353 Bool_t status = calibPed->AddInfo(calPed);
354 if (status) { changes++; }
364 // Store updated pedestal entry to OCDB
367 AliCDBMetaData metaData;
368 metaData.SetBeamPeriod(0);
369 metaData.SetResponsible(kMetaResponsible);
370 metaData.SetComment(kMetaComment);
372 Bool_t storeOK = StoreReferenceData("Calib", "Pedestals", calibPed, &metaData);
373 if ( !storeOK ) result++;
377 Log ("Error: no entries in input file list!");
378 result = kReturnCodeNoEntries;
384 //______________________________________________________________________________________________
385 UInt_t AliEMCALPreprocessor::ExtractSignal(Int_t sourceFXS)
389 // Read signal file from file exchange server
390 // Keep original entry from OCDB in case no new signal are available
392 AliCaloCalibSignal *calibSig=0;
393 AliCDBEntry* entry = GetFromOCDB("Calib", "LED");
394 if (entry) calibSig = (AliCaloCalibSignal*)entry->GetObject();
395 if ( calibSig==NULL ) {
396 Log("AliEMCALPreprocessor: No previous EMCAL signal entry available.\n");
397 calibSig = new AliCaloCalibSignal(AliCaloCalibSignal::kEmCal);
400 TList* list = GetFileSources(sourceFXS,"signal");
401 if (list && list->GetEntries()>0) {
403 /* DS: 17 June 2008 - commented out this reset to avoid crash in shuttle.
404 Not sure why it occurs, in standalone tests a Reset() seemed to work OK..
406 calibSig->Reset(); // let's make a fresh start before possibly adding stuff below
409 // loop through all files from LDCs
413 while (list->At(index)!=NULL) {
414 TObjString* fileNameEntry = (TObjString*) list->At(index);
415 if (fileNameEntry!=NULL) {
416 TString fileName = GetFile(sourceFXS, "signal",
417 fileNameEntry->GetString().Data());
418 TFile *f = TFile::Open(fileName);
420 Log ("Error opening signal file.");
421 result = kReturnCodeNoObject;
424 AliCaloCalibSignal *calSig;
425 f->GetObject("emcCalibSignal",calSig);
427 Log ("No signal calibration object in file.");
428 result = kReturnCodeNoObject;
431 if ( calSig->GetNEvents()>0 ) {
432 // add info for the modules available in the present file
433 Bool_t status = calibSig->AddInfo(calSig);
434 if (status) { changes++; }
444 // Store updated signal entry to OCDB
447 AliCDBMetaData metaData;
448 metaData.SetBeamPeriod(0);
449 metaData.SetResponsible(kMetaResponsible);
450 metaData.SetComment(kMetaComment);
452 Bool_t storeOK = Store("Calib", "LED", calibSig, &metaData, 0, kFALSE);
453 if ( !storeOK ) result++;
457 Log ("Error: no entries in input file list!");
458 result = kReturnCodeNoEntries;