Fixed AliFMDPreprocessor to use GetFromOCDB. It is still AliFMDParameters
[u/mrichter/AliRoot.git] / FMD / AliFMDPreprocessor.cxx
1 /**************************************************************************
2  * Copyright(c) 2004, 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 /* $Id$ */
16 /** @file    AliFMDPreprocessor.cxx
17     @author  Hans Hjersing Dalsgaard <canute@nbi.dk>
18     @date    Mon Mar 27 12:39:09 2006
19     @brief   Shuttle "preprocessor" for the FMD
20 */
21 //___________________________________________________________________
22 //
23 // The class processes data points from DCS (via Amanada), and DAQ DA
24 // files (via FXS) to make calibration data for the FMD. 
25 //
26 // Data points: 
27 //    *  Nothing yet. 
28 //
29 // DAQ FXS file:
30 //    * pedestals - a (ASCII) Comma Separated Values files with the
31 //                  fields 
32 //                       rcu     DDL number 
33 //                       board   FEC board number 
34 //                       chip    ALTRO chip number on FEC
35 //                       channel ALTRO channel number
36 //                       strip   VA1 strip number
37 //                       sample  Sample number
38 //                       ped     Mean of ADC spectra
39 //                       noise   Spread of ADC spectra
40 //                       mu      Mean of Gaussian fit to ADC spectra
41 //                       sigma   Variance of Gaussian fit to ADC spectra
42 //                       chi2    Chi^2 per degrees of freedom of fit
43 //    * Gains     - a (ASCII) Comma Separated Values files with the
44 //                  fields 
45 //                       rcu     DDL number 
46 //                       board   FEC board number 
47 //                       chip    ALTRO chip number on FEC
48 //                       channel ALTRO channel number
49 //                       strip   VA1 strip number
50 //                       gain    Slope of gain
51 //                       error   Error on gain
52 //                       chi2    Chi^2 per degrees of freedom of fit
53 //                  
54 // Latest changes by Christian Holm Christensen
55 //
56
57 // #include <iostream>
58
59 #include <fstream>
60 #include "AliFMDPreprocessor.h"
61 #include "AliFMDCalibPedestal.h"
62 #include "AliFMDCalibGain.h"
63 #include "AliFMDCalibStripRange.h"
64 #include "AliFMDCalibSampleRate.h"
65 #include "AliFMDParameters.h"
66 #include "AliCDBMetaData.h"
67 #include "AliCDBManager.h"
68 // #include "AliDCSValue.h"
69 #include "AliLog.h"
70 #include <TTimeStamp.h>
71 // #include <TFile.h>
72 #include <TObjString.h>
73 #include <TString.h>
74 #include <TNamed.h>
75
76
77 ClassImp(AliFMDPreprocessor)
78 #if 0 // Do not remove - here to make Emacs happy
79 ;
80 #endif 
81
82 //____________________________________________________
83 Bool_t AliFMDPreprocessor::GetAndCheckFileSources(TList*&     list,
84                                                   Int_t       system, 
85                                                   const char* id) 
86 {
87   // Convinience function 
88   // Parameters: 
89   //   list     On return, list of files. 
90   //   system   Alice system (DAQ, DCS, ...)
91   //   id       File id
92   // Return:
93   //   kTRUE on success. 
94   list = GetFileSources(system, id);
95   if (!list) { 
96     TString sys;
97     switch (system) { 
98     case kDAQ: sys = "DAQ";     break;
99     case kDCS: sys = "DCS";     break;
100     default:   sys = "unknown"; break;
101     }
102     Log(Form("Failed to get file sources for %s/%s", sys.Data(), system));
103     return kFALSE;
104   }
105   return kTRUE;
106 }
107
108 //____________________________________________________
109 AliCDBEntry* 
110 AliFMDPreprocessor::GetFromCDB(const char* second, const char* third)
111 {
112   return GetFromOCDB(second, third);
113 }
114
115
116 //____________________________________________________
117 UInt_t AliFMDPreprocessor::Process(TMap* /* dcsAliasMap */)
118 {
119   // Main member function. 
120   // Parameters: 
121   //    dcsAliassMap   Map of DCS data point aliases.
122   // Return 
123   //    ? 
124
125   // Do we need this ?
126   // if(!dcsAliasMap) return 1;
127   // 
128   // Invoking the cdb manager and the FMD parameters class
129   // AliCDBManager* cdb   = AliCDBManager::Instance();
130   // cdb->SetDefaultStorage("local://$ALICE_ROOT");
131   // cdb->SetRun(0);
132   AliFMDParameters* pars = AliFMDParameters::Instance();
133   pars->Init(this, false, AliFMDParameters::kAltroMap);
134
135   // This is if the SOR contains Fee parameters, and we run a DA to
136   // extract these parameters.   The same code could work if we get
137   // the information from DCS via the FXS 
138   TList* files = 0;
139   GetAndCheckFileSources(files, kDAQ,"info");
140   // if (!ifiles) return 1;
141   AliFMDCalibSampleRate*      calibRate  = 0;
142   AliFMDCalibStripRange*      calibRange = 0;
143   AliFMDCalibZeroSuppression* calibZero  = 0;
144   GetInfoCalibration(files, calibRate, calibRange, calibZero);
145   
146   // Gt the run type 
147   TString runType(GetRunType()); 
148
149   //Creating calibration objects
150   AliFMDCalibPedestal* calibPed  = 0;
151   AliFMDCalibGain*     calibGain = 0;
152   if (runType.Contains("PEDESTAL", TString::kIgnoreCase) && 
153       GetAndCheckFileSources(files, kDAQ, "pedestal"))  
154     calibPed = GetPedestalCalibration(files);
155   if (runType.Contains("PULSER", TString::kIgnoreCase) && 
156       GetAndCheckFileSources(files, kDAQ, "gain"))
157     calibGain = GetGainCalibration(files);
158
159   
160   
161   //Storing Calibration objects  
162   AliCDBMetaData metaData;
163   metaData.SetBeamPeriod(0);
164   metaData.SetResponsible("Hans H. Dalsgaard");
165   metaData.SetComment("Preprocessor stores pedestals and gains for the FMD.");
166   
167   Bool_t resultPed   = kFALSE;
168   Bool_t resultGain  = kFALSE;
169   Bool_t resultRange = kFALSE;
170   Bool_t resultRate  = kFALSE;
171   Bool_t resultZero  = kFALSE;
172   if(calibPed)  { 
173     resultPed  = Store("Calib","Pedestal", calibPed, &metaData);
174     delete calibPed;
175   }
176   if(calibGain) { 
177     resultGain = Store("Calib","PulseGain", calibGain, &metaData);
178     delete calibGain;
179   }
180   if(calibRange) { 
181     resultRange = Store("Calib","StripRange", calibRange, &metaData);
182     delete calibRange;
183   }
184   if(calibRate) { 
185     resultRate = Store("Calib","SampleRate", calibRate, &metaData);
186     delete calibRate;
187   }
188   if(calibZero) { 
189     resultZero = Store("Calib","ZeroSuppression", calibZero, &metaData);
190     delete calibZero;
191   }
192   
193   return (resultPed && resultGain ? 0 : 1);
194 #if 0
195   // Disabled until we implement GetInfoCalibration properly
196   return (resultPed   && 
197           resultGain  && 
198           resultRange && 
199           resultRate  &&
200           resultZero 
201           ? 0 : 1);
202 #endif
203 }
204
205 //____________________________________________________________________
206 Bool_t
207 AliFMDPreprocessor::GetInfoCalibration(TList* files, 
208                                        AliFMDCalibSampleRate*&      s,
209                                        AliFMDCalibStripRange*&      r, 
210                                        AliFMDCalibZeroSuppression*& z)
211 {
212   // Get info calibrations. 
213   // Parameters:
214   //     files List of files. 
215   //     s     On return, newly allocated object 
216   //     r     On return, newly allocated object 
217   //     z     On return, newly allocated object 
218   // Return: 
219   //     kTRUE on success
220   if (!files) return kTRUE; // Should really be false
221   if (files->GetEntries() <= 0) return kTRUE;
222   
223   s = new AliFMDCalibSampleRate();
224   r = new AliFMDCalibStripRange();
225   z = new AliFMDCalibZeroSuppression();
226   
227   // AliFMDParameters*    pars     = AliFMDParameters::Instance();
228   TIter                iter(files);
229   TObjString*          fileSource;
230
231   while((fileSource = dynamic_cast<TObjString*>(iter.Next()))) {
232     const Char_t* filename = GetFile(kDAQ, "info", fileSource->GetName());
233     std::ifstream in(filename);
234     if(!in) {
235       AliError(Form("File %s not found!", filename));
236       continue;
237     }
238   }
239   return kTRUE;
240 }
241
242   
243 //____________________________________________________________________
244 AliFMDCalibPedestal* 
245 AliFMDPreprocessor::GetPedestalCalibration(TList* pedFiles)
246 {
247   // Read DAQ DA produced CSV files of pedestals, and return a
248   // calibration object. 
249   // Parameters:
250   //   pedFiles     List of pedestal files 
251   // Return 
252   //   A pointer to a newly allocated AliFMDCalibPedestal object, or
253   //   null in case of errors. 
254   if(!pedFiles) return 0;
255
256   AliFMDCalibPedestal* calibPed = new AliFMDCalibPedestal();
257   AliFMDParameters*    pars     = AliFMDParameters::Instance();
258   TIter                iter(pedFiles);
259   TObjString*          fileSource;
260   
261   while((fileSource = dynamic_cast<TObjString*>(iter.Next()))) {
262     const Char_t* filename = GetFile(kDAQ, "pedestal", fileSource->GetName());
263     std::ifstream in(filename);
264     if(!in) {
265       AliError(Form("File %s not found!", filename));
266       continue;
267     }
268
269     // Get header (how long is it ?)
270     TString header;
271     header.ReadLine(in);
272     header.ToLower();
273     if(!header.Contains("pedestal")) {
274       AliError("File header is not from pedestal!");
275       continue;
276     }
277     Log("File contains data from pedestals");
278     
279     // Read columns line
280     int lineno = 2;
281     header.ReadLine(in);
282     
283     // Loop until EOF
284     while(!in.eof()) {
285       if(in.bad()) { 
286         AliError(Form("Bad read at line %d in %s", lineno, filename));
287         break;
288       }
289       UInt_t ddl=2, board, chip, channel, strip, tb;
290       Float_t ped, noise, mu, sigma, chi2ndf;
291       Char_t c[10];
292           
293       in // >> ddl      >> c[0] 
294          >> board    >> c[1]
295          >> chip     >> c[2]
296          >> channel  >> c[3]
297          >> strip    >> c[4]
298          >> tb       >> c[5]
299          >> ped      >> c[6]
300          >> noise    >> c[7]
301          >> mu       >> c[8]
302          >> sigma    >> c[9]
303          >> chi2ndf;
304       lineno++;
305       // Ignore trailing garbage 
306       if (strip > 127) continue;
307       
308       //Setting the pedestals via the hardware address
309       UShort_t det, sec, str;
310       Char_t ring;
311           
312       pars->Hardware2Detector(ddl,board,chip,channel,det,ring,sec,str);
313       strip += str;
314       calibPed->Set(det,ring,sec,strip,ped,noise);
315     }
316   }
317   return calibPed;
318 }       
319
320 //____________________________________________________________________
321 AliFMDCalibGain* 
322 AliFMDPreprocessor::GetGainCalibration(TList* gainFiles)
323 {
324   // Read DAQ DA produced CSV files of pedestals, and return a
325   // calibration object. 
326   // Parameters:
327   //   pedFiles     List of pedestal files 
328   // Return 
329   //   A pointer to a newly allocated AliFMDCalibPedestal object, or
330   //   null in case of errors. 
331   if(!gainFiles) return 0;
332   
333   AliFMDCalibGain*  calibGain  = new AliFMDCalibGain();
334   AliFMDParameters* pars       = AliFMDParameters::Instance();
335   TIter             iter(gainFiles);
336   TObjString*       fileSource;
337   while((fileSource = dynamic_cast<TObjString *>(iter.Next()))) {
338     const Char_t* filename = GetFile(kDAQ, "gain", fileSource->GetName());
339     std::ifstream in(filename);
340     if(!in) {
341       AliError(Form("File %s not found!", filename));
342       continue;
343     }
344
345     //Get header (how long is it ?)
346     TString header;
347     header.ReadLine(in);
348     header.ToLower();
349     if(!header.Contains("gain")) {
350       AliError("File header is not from gain!");
351       continue;
352     }
353     Log("File contains data from pulse gain");
354
355     // Read column headers
356     header.ReadLine(in);
357
358     int lineno  = 2;
359     // Read until EOF 
360     while(!in.eof()) {
361       if(in.bad()) { 
362         AliError(Form("Bad read at line %d in %s", lineno, filename));
363         break;
364       }
365       UInt_t ddl=2, board, chip, channel, strip;
366       Float_t gain,error,  chi2ndf;
367       Char_t c[7];
368               
369       in // >> ddl      >> c[0] 
370          >> board    >> c[1]
371          >> chip     >> c[2]
372          >> channel  >> c[3]
373          >> strip    >> c[4]
374          >> gain     >> c[5]
375          >> error    >> c[6]
376          >> chi2ndf;
377       lineno++;
378       // Ignore trailing garbage
379       if(strip > 127) continue;
380       
381       //Setting the pedestals via the hardware address
382       UShort_t det, sec, str;
383       Char_t ring;
384       pars->Hardware2Detector(ddl,board,chip,channel,det,ring,sec,str);
385
386       strip += str;
387       calibGain->Set(det,ring,sec,strip,gain);
388     }
389   }
390   return calibGain;
391 }
392
393 //____________________________________________________________________
394 //
395 // EOF
396 //