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