changing return type from Bool_t to expected UInt_t
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALPreprocessor.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 // EMCAL Preprocessor class. It runs by Shuttle at the end of the run,
20 // calculates stuff to be posted in OCDB
21 //
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 ///////////////////////////////////////////////////////////////////////////////
26
27 //Root
28 #include "TFile.h"
29 #include "TTree.h"
30 #include "TEnv.h"
31 #include "TParameter.h"
32
33 #include <TTimeStamp.h>
34
35 //AliRoot
36 #include "AliShuttleInterface.h"
37 #include "AliEMCALPreprocessor.h"
38 #include "AliLog.h"
39 #include "AliDCSValue.h"
40 #include "AliCDBMetaData.h"
41 #include "AliEMCALTriggerTRUDCSConfig.h"
42 #include "AliEMCALTriggerDCSConfig.h"
43 #include "AliCaloCalibPedestal.h"
44 #include "AliCaloCalibSignal.h"
45 #include "AliEMCALSensorTempArray.h"
46
47 const Int_t kValCutTemp = 100;               // discard temperatures > 100 degrees
48 const Int_t kDiffCutTemp = 5;                // discard temperature differences > 5 degrees
49 const TString kPedestalRunType = "PEDESTAL";  // pedestal run identifier
50 const TString kPhysicsRunType = "PHYSICS";   // physics run identifier
51 const TString kStandAloneRunType = "STANDALONE_BC"; // standalone run identifier
52 const TString kAmandaTemp = "PT_%02d.Temperature"; // Amanda string for temperature entries
53 //const Double_t kFitFraction = 0.7;                 // Fraction of DCS sensor fits required 
54 const Double_t kFitFraction = -1.0;          // Don't require minimum number of fits during commissioning 
55
56 const TString kMetaResponsible = "David Silvermyr";
57 //legacy comments and return codes from TPC
58 const TString kMetaComment = "Preprocessor AliEMCAL data base entries.";
59 const int kReturnCodeNoInfo = 9;
60 const int kReturnCodeNoObject = 2;
61 const int kReturnCodeNoEntries = 1;
62
63 const int kNTRU = 30; // From 2011; 10 SuperModules (SM) * 3 TRU per SM
64
65 ClassImp(AliEMCALPreprocessor)
66   
67 //_______________________________________________________________________________________
68 AliEMCALPreprocessor::AliEMCALPreprocessor() :
69   AliPreprocessor("EMC",0),
70   fConfEnv(0), 
71   fTemp(0), 
72   fConfigOK(kTRUE)
73 {
74   //default constructor
75 }
76
77 //_______________________________________________________________________________________
78 AliEMCALPreprocessor::AliEMCALPreprocessor(AliShuttleInterface* shuttle):
79   AliPreprocessor("EMC",shuttle),
80   fConfEnv(0), 
81   fTemp(0), 
82   fConfigOK(kTRUE)
83 {
84   // Constructor AddRunType(kPedestalRunType);
85   
86   // define run types to be processed
87   AddRunType(kPedestalRunType);
88   AddRunType(kPhysicsRunType);
89 }
90
91 //______________________________________________________________________________________________
92 AliEMCALPreprocessor::AliEMCALPreprocessor(const AliEMCALPreprocessor&  ) :
93   AliPreprocessor("EMCAL",0),
94   fConfEnv(0), fTemp(0), fConfigOK(kTRUE)
95 {
96   Fatal("AliEMCALPreprocessor", "copy constructor not implemented");
97 }
98
99 // assignment operator; use copy ctor to make life easy.
100 //______________________________________________________________________________________________
101 AliEMCALPreprocessor& AliEMCALPreprocessor::operator = (const AliEMCALPreprocessor &source ) 
102 {
103   // assignment operator; use copy ctor
104   if (&source == this) return *this;
105   
106   new (this) AliEMCALPreprocessor(source);
107   return *this;
108 }
109
110 //____________________________________________________________________________
111 AliEMCALPreprocessor::~AliEMCALPreprocessor()
112 {
113   // destructor
114   if (fTemp) delete fTemp;
115 }
116
117 //______________________________________________________________________________________________
118 void AliEMCALPreprocessor::Initialize(Int_t run, UInt_t startTime,
119                                       UInt_t endTime)
120 {
121   // Creates AliTestDataDCS object -- start maps half an hour beforre actual run start
122   UInt_t startTimeLocal = startTime-1800;
123   AliPreprocessor::Initialize(run, startTimeLocal, endTime);
124   
125   AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
126                TTimeStamp((time_t)startTime,0).AsString(),
127                TTimeStamp((time_t)endTime,0).AsString()));
128   
129   // Preprocessor configuration
130   AliCDBEntry* entry = GetFromOCDB("Config", "Preprocessor");
131   if (entry) fConfEnv = (TEnv*) entry->GetObject();
132   if ( fConfEnv==0 ) {
133     Log("AliEMCALPreprocessor: Preprocessor Config OCDB entry missing.\n");
134     fConfigOK = kFALSE;
135     return;
136   }
137   
138   // Temperature sensors
139   TTree *confTree = 0;
140   
141   TString tempConf = fConfEnv->GetValue("Temperature","ON");
142   tempConf.ToUpper();
143   if (tempConf != "OFF" ) {
144     entry = GetFromOCDB("Config", "Temperature");
145     if (entry) confTree = (TTree*) entry->GetObject();
146     if ( confTree==0 ) {
147       Log("AliEMCALPreprocessor: Temperature Config OCDB entry missing.\n");
148       fConfigOK = kFALSE;
149       return;
150     }
151     fTemp = new AliEMCALSensorTempArray(startTimeLocal, fEndTime, confTree, kAmandaTemp);
152     fTemp->SetValCut(kValCutTemp);
153     fTemp->SetDiffCut(kDiffCutTemp);
154   }
155   
156   return;
157 }
158
159 //______________________________________________________________________________________________
160 UInt_t AliEMCALPreprocessor::Process(TMap* dcsAliasMap)
161 {
162   // Fills data into EMCAL calibrations objects
163   // Amanda servers provide information directly through dcsAliasMap
164   
165   if (!fConfigOK) return kReturnCodeNoInfo;
166   UInt_t result = 0;
167   TObjArray *resultArray = new TObjArray();
168   TString errorHandling = fConfEnv->GetValue("ErrorHandling","ON");
169   errorHandling.ToUpper();
170   TObject * status;
171   
172   UInt_t dcsResult=0;
173   if (errorHandling == "OFF" ) {
174     if (!dcsAliasMap) dcsResult = kReturnCodeNoEntries;
175     else if (dcsAliasMap->GetEntries() == 0 ) dcsResult = kReturnCodeNoEntries;  
176     status = new TParameter<int>("dcsResult",dcsResult);
177     resultArray->Add(status);
178   } 
179   else {
180     if (!dcsAliasMap) return kReturnCodeNoInfo;
181     else if (dcsAliasMap->GetEntries() == 0 ) return kReturnCodeNoInfo;
182   }
183   
184     
185   TString runType = GetRunType();
186   
187   // Temperature sensors are processed by AliEMCALCalTemp
188   TString tempConf = fConfEnv->GetValue("Temperature","ON");
189   tempConf.ToUpper();
190   if (tempConf != "OFF" && dcsAliasMap ) {
191     UInt_t tempResult = MapTemperature(dcsAliasMap);
192     result=tempResult;
193     status = new TParameter<int>("tempResult",tempResult);
194     resultArray->Add(status);
195   }
196   // Trigger configuration processing: only for Physics runs
197   TString triggerConf = fConfEnv->GetValue("Trigger","ON");
198   triggerConf.ToUpper();
199   if( runType == kPhysicsRunType || runType == kPedestalRunType ) {
200     //  if (triggerConf != "OFF" && dcsAliasMap ) {
201     if ( dcsAliasMap ) {
202       UInt_t triggerResult = MapTriggerConfig(dcsAliasMap);
203       AliInfo(Form("triggerConf %s\n", triggerConf.Data()));
204       result+=triggerResult;
205       status = new TParameter<int>("triggerResult",triggerResult);
206       resultArray->Add(status);
207     }
208   }
209   
210   // Other calibration information will be retrieved through FXS files
211   //  examples:
212   //    TList* fileSourcesDAQ = GetFile(AliShuttleInterface::kDAQ, "pedestals");
213   //    const char* fileNamePed = GetFile(AliShuttleInterface::kDAQ, "pedestals", "LDC1");
214   //
215   //    TList* fileSourcesHLT = GetFile(AliShuttleInterface::kHLT, "calib");
216   //    const char* fileNameHLT = GetFile(AliShuttleInterface::kHLT, "calib", "LDC1");
217   
218   // PEDESTAL ENTRIES:
219   
220   if ( runType == kPedestalRunType ) {
221     Int_t numSources = 1;
222     Int_t pedestalSource[2] = {AliShuttleInterface::kDAQ, AliShuttleInterface::kHLT} ;
223     TString source = fConfEnv->GetValue("Pedestal","DAQ");
224     source.ToUpper();
225     if (source != "OFF" ) { 
226       if ( source == "HLT") pedestalSource[0] = AliShuttleInterface::kHLT;
227       if (!GetHLTStatus()) pedestalSource[0] = AliShuttleInterface::kDAQ;
228       if (source == "HLTDAQ" ) {
229         numSources=2;
230         pedestalSource[0] = AliShuttleInterface::kHLT;
231         pedestalSource[1] = AliShuttleInterface::kDAQ;
232       }
233       if (source == "DAQHLT" ) numSources=2;
234       UInt_t pedestalResult=0;
235       for (Int_t i=0; i<numSources; i++ ) {     
236         pedestalResult = ExtractPedestals(pedestalSource[i]);
237         if ( pedestalResult == 0 ) break;
238       }
239       result += pedestalResult;
240       status = new TParameter<int>("pedestalResult",pedestalResult);
241       resultArray->Add(status);
242     }
243   }
244   
245   // SIGNAL/LED ENTRIES:
246   if( runType == kPhysicsRunType ) {
247     Int_t numSources = 1;
248     Int_t signalSource[2] = {AliShuttleInterface::kDAQ,AliShuttleInterface::kHLT} ;
249     TString source = fConfEnv->GetValue("Signal","DAQ");
250     source.ToUpper();
251     if ( source != "OFF") { 
252       if ( source == "HLT") signalSource[0] = AliShuttleInterface::kHLT;
253       if (!GetHLTStatus()) signalSource[0] = AliShuttleInterface::kDAQ;
254       if (source == "HLTDAQ" ) {
255         numSources=2;
256         signalSource[0] = AliShuttleInterface::kHLT;
257         signalSource[1] = AliShuttleInterface::kDAQ;
258       }
259       if (source == "DAQHLT" ) numSources=2;
260       UInt_t signalResult=0;
261       for (Int_t i=0; i<numSources; i++ ) {     
262         signalResult = ExtractSignal(signalSource[i]);
263         if ( signalResult == 0 ) break;
264       }
265       result += signalResult;
266       status = new TParameter<int>("signalResult",signalResult);
267       resultArray->Add(status);
268     }
269   }
270   
271   
272   // overall status at the end
273   if (errorHandling == "OFF" ) {
274     AliCDBMetaData metaData;
275     metaData.SetBeamPeriod(0);
276     metaData.SetResponsible(kMetaResponsible);
277     metaData.SetComment("Preprocessor AliEMCAL status.");
278     Bool_t storeOK = Store("Calib", "PreprocStatus", resultArray, &metaData, 0, kFALSE);
279     resultArray->Delete();
280     result = 0;
281     if ( !storeOK )  result=1;
282     return result;
283   } 
284   else { 
285     return result;
286   }
287   
288 }
289 //______________________________________________________________________________________________
290 UInt_t AliEMCALPreprocessor::MapTemperature(TMap* dcsAliasMap)
291 { // extract DCS temperature maps. Perform fits to save space
292   UInt_t result=0;
293
294   TMap *map = fTemp->ExtractDCS(dcsAliasMap);
295   if (map) {
296     fTemp->MakeSplineFit(map);
297     Double_t fitFraction = 1.0*fTemp->NumFits()/fTemp->NumSensors(); 
298     if (fitFraction > kFitFraction ) {
299       AliInfo(Form("Temperature values extracted, fits performed.\n"));
300     } 
301     else { 
302       Log ("Too few temperature maps fitted. \n");
303       result = kReturnCodeNoInfo;
304     }
305   } 
306   else {
307     Log("No temperature map extracted. \n");
308     result = kReturnCodeNoInfo;
309   }
310   delete map;
311   // Now store the final CDB file
312   
313   if ( result == 0 ) { // some info was found
314     AliCDBMetaData metaData;
315     metaData.SetBeamPeriod(0);
316     metaData.SetResponsible(kMetaResponsible);
317     metaData.SetComment(kMetaComment);
318     
319     Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE);
320     if ( !storeOK )  result=1;
321     AliInfo(Form("Temperature info stored. result %d\n", result));
322   }
323   
324   return result;
325 }
326
327 //______________________________________________________________________________________________
328 UInt_t AliEMCALPreprocessor::MapTriggerConfig(TMap* dcsAliasMap)
329 { // extract DCS trigger info
330   AliInfo(Form("Get TRU info from DCS DPs.\n"));
331   Int_t i, iTRU;
332   const Int_t bufsize = 1000;
333   char buf[bufsize];
334
335   AliDCSValue *dcsVal;
336   TObjArray *arrL0ALGSEL, *arrPEAKFINDER, *arrGLOBALTHRESH, *arrCOSMTHRESH;
337   TObjArray *arrMASK[6];
338
339   // overall object to hold STU and DCS config info
340   // DS comment: for now only holds TRU info, i.e. only partially filled
341   // (STU info only in raw data header; unfortunately not also picked up via DCS DPs)
342   AliEMCALTriggerDCSConfig *trigConfig = new AliEMCALTriggerDCSConfig();
343   // allocate space for TRU objects
344   TClonesArray *truArr = new TClonesArray("AliEMCALTriggerTRUDCSConfig", kNTRU);
345   for( iTRU = 0; iTRU < kNTRU; iTRU++){
346     new((*truArr)[iTRU]) AliEMCALTriggerTRUDCSConfig();
347   }
348   trigConfig->SetTRUArr(truArr);
349
350   // loop through all TRUs
351   bool debug = true; // debug flag for AliInfo printouts for each TRU
352   for( iTRU = 0; iTRU < kNTRU; iTRU++){
353     if (debug) AliInfo( Form("iTRU %d \n", iTRU) );
354     // get the shuttled values
355     snprintf( buf, bufsize, "EMC_TRU%02d_L0ALGSEL", iTRU );
356     arrL0ALGSEL = (TObjArray*) dcsAliasMap->GetValue( buf );
357     snprintf( buf, bufsize, "EMC_TRU%02d_PEAKFINDER", iTRU );
358     arrPEAKFINDER = (TObjArray*) dcsAliasMap->GetValue( buf );
359     snprintf( buf, bufsize, "EMC_TRU%02d_GLOBALTHRESH", iTRU );
360     arrGLOBALTHRESH = (TObjArray*) dcsAliasMap->GetValue( buf );
361     snprintf( buf, bufsize, "EMC_TRU%02d_COSMTHRESH", iTRU );
362     arrCOSMTHRESH = (TObjArray*) dcsAliasMap->GetValue( buf );
363     
364     for( i = 0; i < 6; i++ ){
365       snprintf( buf, bufsize, "EMC_TRU%02d_MASK%d", iTRU, i );
366       arrMASK[i] = (TObjArray*) dcsAliasMap->GetValue( buf );
367     }
368     
369     // fill the objects
370     AliEMCALTriggerTRUDCSConfig* truConfig = trigConfig->GetTRUDCSConfig(iTRU);
371     if( ! truConfig ){
372       AliWarning( Form("EMC DCS TRU%02d config not retrieved!\n", iTRU ));
373       continue;
374     }
375
376     // get last entries. fill the TRU object
377     if( ! arrL0ALGSEL ){
378       AliWarning( Form("EMC DCS TRU%02d L0ALGSEL alias not found!\n", iTRU ));
379     }
380     else{
381       if (debug) AliInfo( Form("arrL0ALGSEL has %d entries \n", arrL0ALGSEL->GetEntries()) );
382       if ( arrL0ALGSEL->GetEntries() > 0 ) {
383         dcsVal = (AliDCSValue *) arrL0ALGSEL->At( arrL0ALGSEL->GetEntries() - 1 );
384         if (dcsVal) truConfig->SetL0SEL( dcsVal->GetUInt() );
385       }
386     }
387
388     if( ! arrPEAKFINDER ){
389       AliWarning( Form("EMC DCS TRU%02d PEAKFINDER alias not found!\n", iTRU ));
390     }
391     else{
392       if (debug) AliInfo( Form("arrPEAKFINDER has %d entries \n", arrPEAKFINDER->GetEntries()) );
393       if ( arrPEAKFINDER->GetEntries() > 0 ) {
394         dcsVal = (AliDCSValue *) arrPEAKFINDER->At( arrPEAKFINDER->GetEntries() - 1 );
395         if (dcsVal) truConfig->SetSELPF( dcsVal->GetUInt() );
396       }
397     }
398
399     if( ! arrGLOBALTHRESH ){
400       AliWarning( Form("EMC DCS TRU%02d GLOBALTHRESH alias not found!\n", iTRU ));
401     }
402     else{
403       if (debug) AliInfo( Form("arrGLOBALTHRESH has %d entries \n", arrGLOBALTHRESH->GetEntries()) );
404       if ( arrGLOBALTHRESH->GetEntries() > 0 ) {
405         dcsVal = (AliDCSValue *) arrGLOBALTHRESH->At( arrGLOBALTHRESH->GetEntries() - 1 );
406         if (dcsVal) truConfig->SetGTHRL0( dcsVal->GetUInt() );
407       }
408     }
409
410     if( ! arrCOSMTHRESH ){
411       AliWarning( Form("EMC DCS TRU%02d COSMTHRESH alias not found!\n", iTRU ));
412     }
413     else{
414       if (debug) AliInfo( Form("arrCOSMTHRESH has %d entries \n", arrCOSMTHRESH->GetEntries()) );
415       if ( arrCOSMTHRESH->GetEntries() > 0 ) {
416         dcsVal = (AliDCSValue *) arrCOSMTHRESH->At( arrCOSMTHRESH->GetEntries() - 1 );
417         if (dcsVal) truConfig->SetL0COSM( dcsVal->GetUInt() );
418       }
419     }
420     
421     for( i = 0; i < 6; i++ ){
422       if( ! arrMASK[i] ){
423         AliWarning( Form("EMC DCS TRU%02d MASK%d alias not found!\n", iTRU, i ));
424       }
425       else{
426         if (debug) AliInfo( Form("arrMASK[%d] has %d entries \n", i, arrMASK[i]->GetEntries()) );
427         if ( arrMASK[i]->GetEntries() > 0 ) {
428           dcsVal = (AliDCSValue *) arrMASK[i]->At( arrMASK[i]->GetEntries() - 1 );
429           if (dcsVal) truConfig->SetMaskReg( dcsVal->GetUInt(), i );
430         }
431       }
432     }
433     
434   } // TRUs
435   AliInfo(Form("TRU info retrieved.\n"));
436   // save the objects
437   AliCDBMetaData metaData;
438   metaData.SetBeamPeriod(0);
439   metaData.SetResponsible(kMetaResponsible);
440   metaData.SetComment(kMetaComment); 
441       
442   UInt_t result=0;
443   Bool_t storeOK = Store("Calib", "Trigger", trigConfig, &metaData, 0, kFALSE);
444   if ( !storeOK )  result=1;
445   AliInfo(Form("TRU info stored. result %d\n", result));
446
447   return result;
448 }
449
450 //______________________________________________________________________________________________
451 UInt_t AliEMCALPreprocessor::ExtractPedestals(Int_t sourceFXS)
452 {
453   //  Read pedestal file from file exchange server
454   //  Only store if new pedestal info is available
455   //
456   UInt_t result=0;
457
458   AliCaloCalibPedestal *calibPed = new AliCaloCalibPedestal(AliCaloCalibPedestal::kEmCal);
459   calibPed->Init();
460
461   TList* list = GetFileSources(sourceFXS,"pedestals");
462   if (list && list->GetEntries()>0) {
463     
464     //  loop through all files from LDCs
465
466     int changes = 0;
467     UInt_t index = 0;
468     while (list->At(index)!=NULL) {
469       TObjString* fileNameEntry = (TObjString*) list->At(index);
470       if (fileNameEntry!=NULL) {
471         TString fileName = GetFile(sourceFXS, "pedestals",
472                                    fileNameEntry->GetString().Data());
473         TFile *f = TFile::Open(fileName);
474         if (!f) {
475           Log ("Error opening pedestal file.");
476           result = kReturnCodeNoObject;
477           break;
478         }
479         AliCaloCalibPedestal *calPed;
480         f->GetObject("emcCalibPedestal",calPed);
481         if ( !calPed ) {
482           Log ("No pedestal calibration object in file.");
483           result = kReturnCodeNoObject;
484           break;
485         }
486         if ( calPed->GetNEvents()>0 && calPed->GetNChanFills()>0 ) {
487           // add info for the modules available in the present file
488           Bool_t status = calibPed->AddInfo(calPed);
489           if (status) { changes++; }
490         }
491         
492         delete calPed; 
493         f->Close();
494       }
495       index++;
496     }  // while(list)
497     
498     //
499     //  Store updated pedestal entry to OCDB
500     //
501     if (changes>0) {
502       AliCDBMetaData metaData;
503       metaData.SetBeamPeriod(0);
504       metaData.SetResponsible(kMetaResponsible);
505       metaData.SetComment(kMetaComment); 
506       
507       Bool_t storeOK = StoreReferenceData("Calib", "Pedestals", calibPed, &metaData);
508       if ( !storeOK ) result++;
509     }
510   } 
511   else {
512     Log ("Error: no entries in input file list!");
513     result = kReturnCodeNoEntries;
514   }
515   
516   return result;
517 }
518
519 //______________________________________________________________________________________________
520 UInt_t AliEMCALPreprocessor::ExtractSignal(Int_t sourceFXS)
521 { //  Read signal file from file exchange server
522   //  Only store if new signal info is available
523   //
524   UInt_t result=0;
525   AliCaloCalibSignal *calibSig = new AliCaloCalibSignal(AliCaloCalibSignal::kEmCal); 
526   
527   TList* list = GetFileSources(sourceFXS,"signal");
528   if (list && list->GetEntries()>0) {
529
530     //  loop through all files from LDCs
531     
532     int changes = 0;
533     UInt_t index = 0;
534     while (list->At(index)!=NULL) {
535       TObjString* fileNameEntry = (TObjString*) list->At(index);
536       if (fileNameEntry!=NULL) {
537         TString fileName = GetFile(sourceFXS, "signal",
538                                    fileNameEntry->GetString().Data());
539         TFile *f = TFile::Open(fileName);
540         if (!f) {
541           Log ("Error opening signal file.");
542           result = kReturnCodeNoObject;
543           break;
544         }
545         AliCaloCalibSignal *calSig;
546         f->GetObject("emcCalibSignal",calSig);
547         if ( !calSig ) {
548           Log ("No signal calibration object in file.");
549           result = kReturnCodeNoObject;
550           break;
551         }
552         if ( calSig->GetNEvents()>0 ) {
553           // add info for the modules available in the present file
554           Bool_t status = calibSig->AddInfo(calSig);
555           if (status) { changes++; }
556         }
557         
558         delete calSig; 
559         f->Close();
560       }
561       index++;
562     }  // while(list)
563     
564     //
565     //  Store updated signal entry to OCDB
566     //
567     if (changes>0) {
568       AliCDBMetaData metaData;
569       metaData.SetBeamPeriod(0);
570       metaData.SetResponsible(kMetaResponsible);
571       metaData.SetComment(kMetaComment); 
572       
573       Bool_t storeOK = Store("Calib", "LED", calibSig, &metaData, 0, kFALSE);
574       if ( !storeOK ) result++;
575     }
576   } 
577   else {
578     Log ("Error: no entries in input file list!");
579     result = kReturnCodeNoEntries;
580   }
581
582   return result;
583 }
584
585