Removing osolete variable. Adding title to trigger configuration obj also
[u/mrichter/AliRoot.git] / STEER / AliGRPPreprocessor.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 //                          Class AliGRPPreprocessor
20 //                  Global Run Parameters (GRP) preprocessor
21 //    Origin: Panos Christakoglou, UOA-CERN, Panos.Christakoglou@cern.ch
22 //    Modified: Ernesto.Lopez.Torres@cern.ch  CEADEN-CERN
23 //-------------------------------------------------------------------------
24
25 #include <TChain.h>
26 #include <TList.h>
27 #include <TMap.h>
28 #include <TObjString.h>
29 #include <TTimeStamp.h>
30 #include <TSystem.h>
31 #include <TFile.h>
32
33 #include "AliGRPPreprocessor.h"
34 #include "AliGRPDCS.h"
35 #include "AliDCSSensorArray.h"
36
37 #include "AliTriggerConfiguration.h"
38 #include "AliTriggerRunScalers.h"
39
40 #include "AliCDBMetaData.h"
41 #include "AliLog.h"
42
43 class AliDCSValue;
44 class AliShuttleInterface;
45
46 #include <TH1.h>
47
48 // needed for ReceivePromptRecoParameters
49 #include <TSQLServer.h>
50 #include <TSQLResult.h>
51 #include <TSQLRow.h>
52 #include <AliCDBManager.h>
53 #include <AliCDBMetaData.h>
54 #include <AliCDBId.h>
55 #include <AliTriggerConfiguration.h>
56
57 const Double_t kFitFraction = 0.7;                 // Fraction of DCS sensor fits required
58
59 ClassImp(AliGRPPreprocessor)
60
61 //_______________________________________________________________
62   const Int_t AliGRPPreprocessor::fgknDAQLbPar = 8; // num parameters in the logbook
63   const Int_t AliGRPPreprocessor::fgknDCSDP = 11;   // number of dcs dps
64   const char* AliGRPPreprocessor::fgkDCSDataPoints[AliGRPPreprocessor::fgknDCSDP] = {
65                    "LHCState",              // missing in DCS
66                    "L3Polarity",
67                    "DipolePolarity",
68                    "LHCLuminosity",         // missing in DCS
69                    "BeamIntensity",         // missing in DCS
70                    "L3Current",
71                    "DipoleCurrent",
72                    "CavernTemperature",
73                    "CavernAtmosPressure",
74                    "gva_cr5AtmosphericPressure",   // missing in DCS
75                    "gva_meyrinAtmosphericPressure" // missing in DCS
76                  };
77                  
78   const Short_t kSensors = 9; // start index position of sensor in DCS DPs
79
80   const char* AliGRPPreprocessor::fgkLHCState[20] = {
81                    "P", "PREPARE",
82                    "J", "PREINJECTION",
83                    "I", "INJECTION",
84                    "F", "FILLING",
85                    "A", "ADJUST",
86                    "U", "UNSTABLE BEAMS",
87                    "S", "STABLE BEAMS",
88                    "D", "BEAM DUMP",
89                    "R", "RECOVER",
90                    "C", "PRECYCLE"
91                  };
92
93   const char* kppError[] = {
94                    "",
95                    "(DAQ logbook ERROR)",
96                    "(DAQ FXS ERROR)",
97                    "(DCS FXS ERROR)",
98                    "(DCS data points ERROR)",
99                    "(Trigger Configuration ERROR)"
100   };
101
102 //_______________________________________________________________
103 AliGRPPreprocessor::AliGRPPreprocessor(AliShuttleInterface* shuttle):
104         AliPreprocessor("GRP",shuttle),  fPressure(0)
105 {
106   // constructor - shuttle must be instantiated!
107
108   AddRunType("PHYSICS");
109 }
110
111 //_______________________________________________________________
112 AliGRPPreprocessor::~AliGRPPreprocessor()
113 {
114   //destructor
115
116   delete fPressure;
117 }
118
119 //_______________________________________________________________
120 void AliGRPPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime)
121 {
122   // Initialize preprocessor
123
124   AliPreprocessor::Initialize(run, startTime, endTime);
125
126   AliInfo("Initialization of the GRP preprocessor.");
127
128   TClonesArray * array = new TClonesArray("AliDCSSensor",2);
129   for(Int_t j = 0; j < 2; j++) {
130     AliDCSSensor * sens = new ((*array)[j])AliDCSSensor;
131     sens->SetStringID(fgkDCSDataPoints[j+kSensors]);
132   }
133   AliInfo(Form("Pressure Entries: %d",array->GetEntries()));
134
135   fPressure = new AliDCSSensorArray(fStartTime, fEndTime, array);
136 }
137
138 //_______________________________________________________________
139 UInt_t AliGRPPreprocessor::Process(TMap* valueMap)
140 {
141   // process data retrieved by the Shuttle
142   
143   //=================//
144   // DAQ logbook     //
145   //=================//
146   UInt_t error = 0;
147   
148   TMap *grpmap = ProcessDaqLB();
149   if( grpmap->GetEntries() == fgknDAQLbPar ) {
150     Log(Form("DAQ logbook, successful!"));
151   } else {
152     Log(Form("DAQ logbook, missing parameters!!!"));
153     error |= 1;
154   }
155   //=================//
156   // DAQ FXS         //
157   //=================//
158   UInt_t iDaqFxs = ProcessDaqFxs();
159   if( iDaqFxs == 0 ) {
160     Log(Form("DAQ FXS, successful!"));
161   } else {
162     Log(Form("DAQ FXS, could not store run raw tag file!!!"));
163     error |= 2;
164   }
165   
166   //=================//
167   // DCS FXS         //
168   //=================//
169   UInt_t iDcsFxs = ProcessDcsFxs();
170   if( iDcsFxs == 0 ) {
171      Log(Form("DCS FXS, successful!"));
172   } else {
173      Log(Form("DCS FXS, Could not store CTP run configuration and scalers!!!"));
174     error |= 4;
175   }
176   
177   //=================//
178   // DCS data points //
179   //=================//
180   Int_t entries = ProcessDcsDPs( valueMap, grpmap );
181   if( entries < fgknDCSDP-5 ) { // FIXME (!= ) LHState and pressure map are not working yet...
182     Log(Form("Problem with the DCS data points!!!"));
183     error |= 8;
184   } else  Log(Form("DCS data points, successful!"));
185
186   //=======================//
187   // Trigger Configuration //
188   //=======================//
189   // either from DAQ logbook.....
190   const char * triggerConf = GetTriggerConfiguration();
191   if (triggerConf!= NULL) {
192     Log("Found trigger configuration in DAQ logbook");
193     AliTriggerConfiguration *runcfg = AliTriggerConfiguration::LoadConfigurationFromString(triggerConf);          
194     if (!runcfg) {
195       Log("Bad CTP run configuration file from DAQ logbook! The corresponding CDB entry will not be filled!");
196       error |= 16;
197     }
198     else {
199       TString titleCTPcfg = Form("CTP cfg for run %i from DAQ",fRun);
200       runcfg->SetTitle(titleCTPcfg);
201       AliCDBMetaData metaData;
202       metaData.SetBeamPeriod(0);
203       metaData.SetResponsible("Roman Lietava");
204       metaData.SetComment("CTP run configuration from DAQ logbook");
205       if (!Store("CTP","Config", runcfg, &metaData, 0, 0)) {
206         Log("Unable to store the CTP run configuration object to OCDB!");
207         error |= 16;
208       }
209     }
210   }
211   // ...or from DCS FXS
212   else{
213      Log("No trigger configuration found in the DAQ logbook!! Trying reading from DCS FXS...");
214      TString runcfgfile = GetFile(kDCS, "CTP_runconfig", "");
215      if (runcfgfile.IsNull()) {
216        Log("No CTP runconfig files has been found in DCS FXS!");
217        error |= 16;
218      }
219      else {
220        Log(Form("File with Id CTP_runconfig found! Copied to %s",runcfgfile.Data()));
221        AliTriggerConfiguration *runcfg = AliTriggerConfiguration::LoadConfiguration(runcfgfile);
222        if (!runcfg) {
223          Log("Bad CTP run configuration file from DCS FXS! The corresponding CDB entry will not be filled!");
224          error |= 16;;
225        }
226        else {
227          TString titleCTPcfg = Form("CTP cfg for run %i from DCS",fRun);
228          runcfg->SetTitle(titleCTPcfg);
229          AliCDBMetaData metaData;
230          metaData.SetBeamPeriod(0);
231          metaData.SetResponsible("Roman Lietava");
232          metaData.SetComment("CTP run configuration from DCS FXS");
233          if (!Store("CTP","Config", runcfg, &metaData, 0, 0)) {
234            Log("Unable to store the CTP run configuration object to OCDB!");
235            error |= 16;
236          }
237        }
238      }
239   }
240
241   grpmap->SetOwner(1);
242   AliInfo(Form("Final list entries: %d",grpmap->GetEntries()));
243   
244   AliCDBMetaData md;
245   md.SetResponsible("Ernesto Lopez Torres");
246   md.SetComment("Output parameters from the GRP preprocessor.");
247   
248   Bool_t result = Store("GRP", "Data", grpmap, &md);
249   
250   delete grpmap;
251   
252   if (result && !error ) {
253     Log("GRP Preprocessor Success");
254     return 0;
255   } else {
256     Log( Form("GRP Preprocessor FAILS!!! %s%s%s%s%s",
257                                  kppError[(error&1)?1:0],
258                                  kppError[(error&2)?2:0],
259                                  kppError[(error&4)?3:0],
260                                  kppError[(error&8)?4:0],
261                                  kppError[(error&16)?5:0]
262                                   ));
263     return error;
264   }
265 }
266
267 //_______________________________________________________________
268 TMap *AliGRPPreprocessor::ProcessDaqLB()
269 {
270   //Getting the DAQ lb information
271   
272   const char* timeStart         = GetRunParameter("time_start");
273   const char* timeEnd           = GetRunParameter("time_end");
274   const char* beamEnergy        = GetRunParameter("beamEnergy");
275   const char* beamType          = GetRunParameter("beamType");
276   const char* numberOfDetectors = GetRunParameter("numberOfDetectors");
277   const char* detectorMask      = GetRunParameter("detectorMask");
278   const char* lhcPeriod         = GetRunParameter("LHCperiod");
279   
280   TMap *mapDAQ = new TMap();
281   
282   if (timeStart) {
283     Log(Form("Start time for run %d: %s",fRun, timeStart));
284   } else {
285     Log(Form("Start time not put in logbook!"));
286   }
287   mapDAQ->Add(new TObjString("fAliceStartTime"), new TObjString(timeStart));
288
289   if (timeEnd) {
290     Log(Form("End time for run %d: %s",fRun, timeEnd));
291   } else {
292     Log(Form("End time not put in logbook!"));
293   }
294   mapDAQ->Add(new TObjString("fAliceStopTime"), new TObjString(timeEnd));
295
296   if (beamEnergy) {
297     Log(Form("Beam energy for run %d: %s",fRun, beamEnergy));
298   } else {
299     Log(Form("Beam energy not put in logbook!"));
300   }
301   mapDAQ->Add(new TObjString("fAliceBeamEnergy"), new TObjString(beamEnergy));
302
303   if (beamType) {
304     Log(Form("Beam type for run %d: %s",fRun, beamType));
305   } else {
306     Log(Form("Beam type not put in logbook!"));
307   }
308   mapDAQ->Add(new TObjString("fAliceBeamType"), new TObjString(beamType));
309
310   if (numberOfDetectors) {
311     Log(Form("Number of active detectors for run %d: %s",fRun, numberOfDetectors));
312   } else {
313     Log(Form("Number of active detectors not put in logbook!"));
314   }
315   mapDAQ->Add(new TObjString("fNumberOfDetectors"), new TObjString(numberOfDetectors));
316
317   if (detectorMask) {
318     Log(Form("Detector mask for run %d: %s",fRun, detectorMask));
319   } else {
320     Log(Form("Detector mask not put in logbook!"));
321   }
322   mapDAQ->Add(new TObjString("fDetectorMask"), new TObjString(detectorMask));
323
324   if (lhcPeriod) {
325     Log(Form("LHC period (DAQ) for run %d: %s",fRun, lhcPeriod));
326   } else {
327     Log(Form("LHCperiod not put in logbook!"));
328   }
329   mapDAQ->Add(new TObjString("fLHCPeriod"), new TObjString(lhcPeriod));
330   
331   mapDAQ->Add(new TObjString("fRunType"), new TObjString(GetRunType()));
332   Log( Form("Retrived %d parameters from logbook", mapDAQ->GetEntries() ) );
333
334   return mapDAQ;
335 }
336
337 //_______________________________________________________________
338 UInt_t AliGRPPreprocessor::ProcessDaqFxs()
339 {
340   //======DAQ FXS======//
341
342   TList* list = GetFileSources(kDAQ);  
343   if (!list) {
344     Log("No raw data tag list: connection problems with DAQ FXS logbook!");
345     return 1;
346   }
347
348   if (list->GetEntries() == 0) {
349     Log("no raw data tags in this run: nothing to merge!");
350     delete  list; list=0;
351     return 0;
352   }
353
354   TChain *fRawTagChain = new TChain("T");
355   Int_t nFiles=0;
356   TIterator* iter = list->MakeIterator();
357   TObject* obj = 0;
358   while ((obj = iter->Next())) {
359     TObjString* objStr = dynamic_cast<TObjString*> (obj);
360     if (objStr) {
361       Log(Form("Found source %s", objStr->String().Data()));
362       TList* list2 = GetFileIDs(kDAQ, objStr->String());
363       if (!list2) {
364         Log("No list with ids from DAQ was found: connection problems with DAQ FXS logbook!");
365         delete fRawTagChain; fRawTagChain=0;
366         return 1;
367       }
368       Log(Form("Number of ids: %d",list2->GetEntries()));
369       for(Int_t i = 0; i < list2->GetEntries(); i++) {
370         TObjString *idStr = (TObjString *)list2->At(i);
371         TString fileName = GetFile(kDAQ,idStr->String().Data(),objStr->String().Data());
372         if (fileName.Length() > 0) {
373            Log(Form("Adding file in the chain: %s",fileName.Data()));
374            fRawTagChain->Add(fileName.Data());
375            nFiles++;
376         } else {
377            Log(Form("Could not retrieve file with id %s from source %s: "
378                     "connection problems with DAQ FXS!",
379                      idStr->String().Data(), objStr->String().Data()));
380            delete list; list=0;
381            delete list2; list2=0;
382            delete fRawTagChain; fRawTagChain=0;
383            return 2;
384         }
385       }
386       delete list2;
387     }
388   }
389   
390   TString fRawDataFileName = "GRP_Merged.tag.root";
391   Log(Form("Merging %d raw data tags into file: %s", nFiles, fRawDataFileName.Data()));
392   if( fRawTagChain->Merge(fRawDataFileName) < 1 ) {
393     Log("Error merging raw data files!!!");
394     return 3;
395   }
396   
397   TString outputfile = Form("Run%d.Merged.RAW.tag.root", fRun);
398   Bool_t result = StoreRunMetadataFile(fRawDataFileName.Data(),outputfile.Data());
399   
400   if (!result) {
401     Log("Problem storing raw data tags in local file!!!");
402   } else {
403     Log("Raw data tags merged successfully!!");
404   }
405   
406   delete iter;
407   delete list;
408   delete fRawTagChain; fRawTagChain=0;
409
410   if (result == kFALSE) {
411     return 4;
412   }
413
414   return 0;
415
416 }
417
418 //_______________________________________________________________
419 UInt_t AliGRPPreprocessor::ProcessDcsFxs()
420 {
421   //======DCS FXS======//
422   // Get the CTP run configuration
423   // and scalers from DCS FXS
424
425 /*
426   {
427
428     // Get the CTP run configuration
429     TList* list = GetFileSources(kDCS,"CTP_runconfig");  
430     if (!list) {
431       Log("No CTP runconfig file: connection problems with DCS FXS logbook!");
432       return 1;
433     }
434   
435     if (list->GetEntries() == 0) {
436       Log("No CTP runconfig file to be processed!");
437       return 1;
438     }
439     else {
440       TIter iter(list);
441       TObjString *source;
442       while ((source = dynamic_cast<TObjString *> (iter.Next()))) {
443         TString runcfgfile = GetFile(kDCS, "CTP_runconfig", source->GetName());
444         if (runcfgfile.IsNull()) {
445           Log("No CTP runconfig files has been found: empty source!");
446         }
447         else {
448           Log(Form("File with Id CTP_runconfig found in source %s! Copied to %s",source->GetName(),runcfgfile.Data()));
449           AliTriggerConfiguration *runcfg = AliTriggerConfiguration::LoadConfiguration(runcfgfile);
450           if (!runcfg) {
451             Log("Bad CTP run configuration file! The corresponding CDB entry will not be filled!");
452             return 1;
453           }
454           else {
455             AliCDBMetaData metaData;
456             metaData.SetBeamPeriod(0);
457             metaData.SetResponsible("Roman Lietava");
458             metaData.SetComment("CTP run configuration");
459             if (!Store("CTP","Config", runcfg, &metaData, 0, 0)) {
460               Log("Unable to store the CTP run configuration object to OCDB!");
461             }
462           }
463         }
464       }
465     }
466     delete list;
467   }
468 */
469   {
470     // Get the CTP counters information
471     TList* list = GetFileSources(kDCS,"CTP_xcounters");  
472     if (!list) {
473       Log("No CTP counters file: connection problems with DAQ FXS logbook!");
474       return 1;
475     }
476   
477     if (list->GetEntries() == 0) {
478       Log("No CTP counters file to be processed!");
479       return 1;
480     }
481     else {
482       TIter iter(list);
483       TObjString *source;
484       while ((source = dynamic_cast<TObjString *> (iter.Next()))) {
485         TString countersfile = GetFile(kDCS, "CTP_xcounters", source->GetName());
486         if (countersfile.IsNull()) {
487           Log("No CTP counters files has been found: empty source!");
488         }
489         else {
490           Log(Form("File with Id CTP_xcounters found in source %s! Copied to %s",source->GetName(),countersfile.Data()));
491           AliTriggerRunScalers *scalers = AliTriggerRunScalers::ReadScalers(countersfile);
492           if (!scalers) {
493             Log("Bad CTP counters file! The corresponding CDB entry will not be filled!");
494             return 1;
495           }
496           else {
497             AliCDBMetaData metaData;
498             metaData.SetBeamPeriod(0);
499             metaData.SetResponsible("Roman Lietava");
500             metaData.SetComment("CTP scalers");
501             if (!Store("CTP","Scalers", scalers, &metaData, 0, 0)) {
502               Log("Unable to store the CTP scalers object to OCDB!");
503             }
504           }
505         }
506       }
507     }
508     delete list;
509   }
510
511   return 0;
512 }
513
514 //_______________________________________________________________
515 Int_t AliGRPPreprocessor::ProcessDcsDPs(TMap* valueMap, TMap* mapDCS)
516 {
517   //Getting the DCS dps
518   //===========//
519
520   //DCS data points
521   //===========//
522   
523   Int_t entries = 0;
524
525   AliInfo(Form("==========LHCState==========="));
526   TObjArray *aliasLHCState = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[0]);
527   if(!aliasLHCState) {
528     Log(Form("LHCState not found!!!"));
529   } else {
530     AliGRPDCS *dcs1 = new AliGRPDCS(aliasLHCState,fStartTime,fEndTime);
531     TString sLHCState = dcs1->ProcessDCS(2);
532     if (sLHCState) {
533       for( Int_t i=0; i<20; i+=2 ) {
534          if( sLHCState.CompareTo(fgkLHCState[i]) == 0 ) {
535             sLHCState = fgkLHCState[i+1];
536             break;
537          }
538       }
539       Log(Form("<LHCState> for run %d: %s",fRun, sLHCState.Data()));
540     } else {
541       Log("LHCState not put in TMap!");
542     }
543     mapDCS->Add(new TObjString("fLHCState"),new TObjString(sLHCState));
544     ++entries;
545   }
546
547   AliInfo(Form("==========L3Polarity==========="));
548   TObjArray *aliasL3Polarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[1]);
549   if(!aliasL3Polarity) {
550     Log(Form("L3Polarity not found!!!"));
551   } else {
552     AliGRPDCS *dcs6 = new AliGRPDCS(aliasL3Polarity,fStartTime,fEndTime);
553     TString sL3Polarity = dcs6->ProcessDCS(1);  
554     if (sL3Polarity) {
555       Log(Form("<L3Polarity> for run %d: %s",fRun, sL3Polarity.Data()));
556     } else {
557       Log("L3Polarity not put in TMap!");
558     }
559     mapDCS->Add(new TObjString("fL3Polarity"),new TObjString(sL3Polarity));
560     ++entries;
561   }
562
563   AliInfo(Form("==========DipolePolarity==========="));
564   TObjArray *aliasDipolePolarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[2]);
565   if(!aliasDipolePolarity) {
566     Log(Form("DipolePolarity not found!!!"));
567   } else {
568     AliGRPDCS *dcs8 = new AliGRPDCS(aliasDipolePolarity,fStartTime,fEndTime);
569     TString sDipolePolarity = dcs8->ProcessDCS(1);  
570     if (sDipolePolarity) {
571       Log(Form("<DipolePolarity> for run %d: %s",fRun, sDipolePolarity.Data()));
572     } else {
573       Log("DipolePolarity not put in TMap!");
574     }
575     mapDCS->Add(new TObjString("fDipolePolarity"),new TObjString(sDipolePolarity));
576     ++entries;
577   }
578
579   AliInfo(Form("==========LHCLuminosity==========="));
580   TObjArray *aliasLHCLuminosity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[3]);
581   if(!aliasLHCLuminosity) {
582     Log(Form("LHCLuminosity not found!!!"));
583   } else {
584     AliGRPDCS *dcs3 = new AliGRPDCS(aliasLHCLuminosity,fStartTime,fEndTime);
585     TString sMeanLHCLuminosity = dcs3->ProcessDCS(5);  
586     if (sMeanLHCLuminosity) {
587       Log(Form("<LHCLuminosity> for run %d: %s",fRun, sMeanLHCLuminosity.Data()));
588     } else {
589       Log("LHCLuminosity not put in TMap!");
590     }
591     mapDCS->Add(new TObjString("fLHCLuminosity"), new TObjString(sMeanLHCLuminosity));
592     ++entries;
593   }
594
595   AliInfo(Form("==========BeamIntensity==========="));
596   TObjArray *aliasBeamIntensity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[4]);
597   if(!aliasBeamIntensity) {
598     Log(Form("BeamIntensity not found!!!"));
599   } else {
600     AliGRPDCS *dcs4 = new AliGRPDCS(aliasBeamIntensity,fStartTime,fEndTime);
601     TString sMeanBeamIntensity = dcs4->ProcessDCS(5);  
602     if (sMeanBeamIntensity) {
603       Log(Form("<BeamIntensity> for run %d: %s",fRun, sMeanBeamIntensity.Data()));
604     } else {
605       Log("BeamIntensity not put in TMap!");
606     }
607     mapDCS->Add(new TObjString("fBeamIntensity"),new TObjString(sMeanBeamIntensity));
608     ++entries;
609   }
610
611   AliInfo(Form("==========L3Current==========="));
612   TObjArray *aliasL3Current = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[5]);
613   if(!aliasL3Current) {
614     Log(Form("L3Current not found!!!"));
615   } else {
616     AliGRPDCS *dcs5 = new AliGRPDCS(aliasL3Current,fStartTime,fEndTime);
617     TString sMeanL3Current = dcs5->ProcessDCS(5);  
618     if (sMeanL3Current) {
619       Log(Form("<L3Current> for run %d: %s",fRun, sMeanL3Current.Data()));
620     } else {
621       Log("L3Current not put in TMap!");
622     }
623     mapDCS->Add(new TObjString("fL3Current"),new TObjString(sMeanL3Current));
624     ++entries;
625   }
626
627
628   AliInfo(Form("==========DipoleCurrent==========="));
629   TObjArray *aliasDipoleCurrent = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[6]);
630   if(!aliasDipoleCurrent) {
631     Log(Form("DipoleCurrent not found!!!"));
632   }  else {
633     AliGRPDCS *dcs7 = new AliGRPDCS(aliasDipoleCurrent,fStartTime,fEndTime);
634     TString sMeanDipoleCurrent = dcs7->ProcessDCS(5);  
635     if (sMeanDipoleCurrent) {
636       Log(Form("<DipoleCurrent> for run %d: %s",fRun, sMeanDipoleCurrent.Data()));
637     } else {
638       Log("DipoleCurrent not put in TMap!");
639     }
640     mapDCS->Add(new TObjString("fDipoleCurrent"),new TObjString(sMeanDipoleCurrent));
641     ++entries;
642   }
643
644   AliInfo(Form("==========CavernTemperature==========="));
645   TObjArray *aliasCavernTemperature = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[7]);
646   if(!aliasCavernTemperature) {
647     Log(Form("CavernTemperature not found!!!"));
648   }  else {
649     AliGRPDCS *dcs9 = new AliGRPDCS(aliasCavernTemperature,fStartTime,fEndTime);
650     TString sMeanCavernTemperature = dcs9->ProcessDCS(5);  
651     if (sMeanCavernTemperature) {
652       Log(Form("<CavernTemperature> for run %d: %s",fRun, sMeanCavernTemperature.Data()));
653     } else {
654       Log("CavernTemperature not put in TMap!");
655     }
656     mapDCS->Add(new TObjString("fCavernTemperature"),new TObjString(sMeanCavernTemperature));
657     ++entries;
658   }
659
660   AliInfo(Form("==========CavernPressure==========="));
661   TObjArray *aliasCavernPressure = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[8]);
662   if(!aliasCavernPressure) {
663     Log("CavernPressure not found!!!");
664   } else {
665     AliGRPDCS *dcs10 = new AliGRPDCS(aliasCavernPressure,fStartTime,fEndTime);
666     TString sMeanCavernPressure = dcs10->ProcessDCS(5);  
667     if (sMeanCavernPressure) {
668       Log(Form("<CavernPressure> for run %d: %s",fRun, sMeanCavernPressure.Data()));
669     } else {
670       Log("CavernPressure not put in TMap!");
671     }
672     mapDCS->Add(new TObjString("fCavernPressure"),new TObjString(sMeanCavernPressure));
673     ++entries;
674   }
675
676
677    // NEEDS TO BE REVISED, CONFIRMED
678    AliInfo(Form("==========GenevaPressureMaps==========="));
679    AliDCSSensorArray *dcsSensorArray = GetPressureMap(valueMap);
680    if( fPressure->NumFits()==0 ) {
681      Log("Problem with the pressure sensor values!!!");
682    } else {
683       AliDCSSensor* sensorCr5 = dcsSensorArray->GetSensor(fgkDCSDataPoints[9]);
684       if( sensorCr5->GetFit() ) {
685         Log(Form("<GvaCr5Pressure> for run %d: Sensor Fit found",fRun));
686         mapDCS->Add( new TObjString("fCr5Pressure"), sensorCr5 );
687         ++entries;
688       } else {
689         Log(Form("ERROR Sensor Fit for %s not found: ", fgkDCSDataPoints[9] ));
690       }
691       
692       AliDCSSensor* sensorMeyrin = dcsSensorArray->GetSensor(fgkDCSDataPoints[10]);
693       if( sensorMeyrin->GetFit() ) {
694         Log(Form("<MeyrinPressure> for run %d: Sensor Fit found",fRun));
695         mapDCS->Add( new TObjString("fMeyrinPressure"), sensorMeyrin );
696         ++entries;
697       } else {
698         Log(Form("ERROR Sensor Fit for %s not found: ", fgkDCSDataPoints[10] ));
699       }
700
701    }
702
703   return entries;
704 }
705
706 //_______________________________________________________________
707 AliDCSSensorArray *AliGRPPreprocessor::GetPressureMap(TMap* dcsAliasMap)
708 {
709   // extract DCS pressure maps. Perform fits to save space
710   
711   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
712   if (map) {
713     fPressure->MakeSplineFit(map);
714     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors(); 
715     if (fitFraction > kFitFraction ) {
716       AliInfo(Form("Pressure values extracted, %d fits performed.", fPressure->NumFits()));
717     } else { 
718       AliInfo("Too few pressure maps fitted!!!");
719     }
720   } else {
721     AliInfo("no atmospheric pressure map extracted!!!");
722   }
723   delete map;
724  
725   return fPressure;
726 }
727
728 //_______________________________________________________________
729 /*UInt_t AliGRPPreprocessor::MapPressure(TMap* dcsAliasMap) {
730   // extract DCS pressure maps. Perform fits to save space
731   
732   UInt_t result=0;
733   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
734   if (map) {
735     fPressure->MakeSplineFit(map);
736     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors(); 
737     if (fitFraction > kFitFraction ) {
738       AliInfo(Form("Pressure values extracted, fits performed.\n"));
739     } else { 
740       Log ("Too few pressure maps fitted. \n");
741       result = 9;
742     }
743   } else {
744     Log("AliTPCPreprocsessor: no atmospheric pressure map extracted. \n");
745     result=9;
746   }
747   delete map;
748   // Now store the final CDB file
749   
750   if ( result == 0 ) {
751     AliCDBMetaData metaData;
752     metaData.SetBeamPeriod(0);
753     metaData.SetResponsible("Panos Christakoglou");
754     metaData.SetComment("Preprocessor AliGRP data base pressure entries.");
755     
756     Bool_t storeOK = Store("Calib", "Pressure", fPressure, &metaData, 0, 0);
757     if ( !storeOK ) result=1; 
758   }
759   
760   return result; 
761   }*/
762
763   
764 //_______________________________________________________________
765 Int_t AliGRPPreprocessor::ReceivePromptRecoParameters(UInt_t run, const char* dbHost, Int_t dbPort, const char* dbName, const char* user, const char* password, const char *cdbRoot)
766 {
767    //
768    // Retrieves logbook and trigger information from the online logbook
769    // This information is needed for prompt reconstruction
770    //
771    // Parameters are:
772    // Run number
773    // DAQ params: dbHost, dbPort, dbName, user, password, logbookTable, triggerTable
774    // cdbRoot
775    //
776    // returns:
777    //         positive on success: the return code is the run number of last run processed of the same run type already processed by the SHUTTLE
778    //         0 on success and no run was found
779    //         negative on error
780    //
781    // This function is NOT called during the preprocessor run in the Shuttle!
782    //
783
784    // defaults
785    if (dbPort == 0)
786      dbPort = 3306;
787
788    // CDB connection
789    AliCDBManager* cdb = AliCDBManager::Instance();
790    cdb->SetDefaultStorage(cdbRoot);
791
792    // SQL connection
793    TSQLServer* server = TSQLServer::Connect(Form("mysql://%s:%d/%s", dbHost, dbPort, dbName), user, password);
794
795    if (!server)
796    {
797      Printf("ERROR: Could not connect to DAQ LB");
798      return -1;
799    }
800
801    // main logbook
802    TString sqlQuery;
803    sqlQuery.Form("SELECT DAQ_time_start, run_type, detectorMask FROM logbook WHERE run = %d", run);
804    TSQLResult* result = server->Query(sqlQuery);
805    if (!result)
806    {
807      Printf("ERROR: Can't execute query <%s>!", sqlQuery.Data());
808      return -2;
809    }
810
811    if (result->GetRowCount() == 0)
812    {
813      Printf("ERROR: Run %d not found", run);
814      delete result;
815      return -3;
816    }
817
818    TSQLRow* row = result->Next();
819    if (!row)
820    {
821      Printf("ERROR: Could not receive data from run %d", run);
822      delete result;
823      return -4;
824    }
825
826    TString runType(row->GetField(1));
827
828    TMap grpData;
829    grpData.Add(new TObjString("DAQ_time_start"), new TObjString(row->GetField(0)));
830    grpData.Add(new TObjString("run_type"), new TObjString(runType));
831    grpData.Add(new TObjString("detectorMask"), new TObjString(row->GetField(2)));
832
833    delete row;
834    row = 0;
835
836    delete result;
837    result = 0;
838
839    Printf("Storing GRP/GRP/Data object with the following content");
840    grpData.Print();
841
842    AliCDBMetaData metadata;
843    metadata.SetResponsible("Jan Fiete Grosse-Oetringhaus");
844    metadata.SetComment("GRP Output parameters received during online running");
845
846    AliCDBId id("GRP/GRP/Data", run, run);
847    Bool_t success = cdb->Put(&grpData, id, &metadata);
848
849    grpData.DeleteAll();
850
851    if (!success)
852    {
853      Printf("ERROR: Could not store GRP/GRP/Data into OCDB");
854      return -5;
855    }
856
857    // Receive trigger information
858    sqlQuery.Form("SELECT configFile FROM logbook_trigger_config WHERE run = %d", run);
859    result = server->Query(sqlQuery);
860    if (!result)
861    {
862      Printf("ERROR: Can't execute query <%s>!", sqlQuery.Data());
863      return -11;
864    }
865
866    if (result->GetRowCount() == 0)
867    {
868      Printf("ERROR: Run %d not found in logbook_trigger_config", run);
869      delete result;
870      return -12;
871    }
872
873    row = result->Next();
874    if (!row)
875    {
876      Printf("ERROR: Could not receive logbook_trigger_config data from run %d", run);
877      delete result;
878      return -13;
879    }
880
881    TString triggerConfig(row->GetField(0));
882
883    delete row;
884    row = 0;
885
886    delete result;
887    result = 0;
888
889    Printf("Found trigger configuration: %s", triggerConfig.Data());
890
891    AliTriggerConfiguration *runcfg = AliTriggerConfiguration::LoadConfigurationFromString(triggerConfig);
892    if (!runcfg)
893    {
894      Printf("ERROR: Could not create CTP configuration object");
895      return -14;
896    }
897
898    metadata.SetComment("CTP run configuration received during online running");
899
900    AliCDBId id2("GRP/CTP/Config", run, run);
901    success = cdb->Put(runcfg, id2, &metadata);
902
903    delete runcfg;
904    runcfg = 0;
905
906    if (!success)
907    {
908      Printf("ERROR: Could not store GRP/CTP/Config into OCDB");
909      return -15;
910    }
911
912    // get last run with same run type that was already processed by the SHUTTLE
913
914    sqlQuery.Form("SELECT max(logbook.run) FROM logbook LEFT JOIN logbook_shuttle ON logbook_shuttle.run = logbook.run WHERE run_type = '%s' AND shuttle_done = 1", runType.Data());
915    result = server->Query(sqlQuery);
916    if (!result)
917    {
918      Printf("ERROR: Can't execute query <%s>!", sqlQuery.Data());
919      return -21;
920    }
921
922    if (result->GetRowCount() == 0)
923    {
924      Printf("ERROR: No result with query <%s>", sqlQuery.Data());
925      delete result;
926      return -22;
927    }
928
929    row = result->Next();
930    if (!row)
931    {
932      Printf("ERROR: Could not receive data for query <%s>", sqlQuery.Data());
933      delete result;
934      return -23;
935    }
936
937    TString lastRunStr(row->GetField(0));
938    Int_t lastRun = lastRunStr.Atoi();
939
940    Printf("Last run with same run type %s is %d", runType.Data(), lastRun);
941
942    delete row;
943    row = 0;
944
945    delete result;
946    result = 0;
947
948    server->Close();
949    delete server;
950    server = 0;
951
952    return lastRun;
953 }