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 //-------------------------------------------------------------------------
  23
  24 #include <TChain.h>
  25 #include <TList.h>
  26 #include <TMap.h>
  27 #include <TObjString.h>
  28 #include <TTimeStamp.h>
  29 #include <TSystem.h>
  30 #include <TFile.h>
  31
  32 #include "AliGRPPreprocessor.h"
  33 #include "AliGRPDCS.h"
  34 #include "AliDCSSensorArray.h"
  35
  36 #include "AliTriggerConfiguration.h"
  37 #include "AliTriggerRunScalers.h"
  38
  39 #include "AliCDBMetaData.h"
  40 #include "AliLog.h"
  41
  42 class AliDCSValue;
  43 class AliShuttleInterface;
  44
  45 #include <TH1.h>
  46
  47 const Double_t kFitFraction = 0.7;                 // Fraction of DCS sensor fits required
  48
  49 ClassImp(AliGRPPreprocessor)
  50
  51 //_______________________________________________________________
  52   const char* AliGRPPreprocessor::fgkDCSDataPoints[12] = {"LHCState","LHCPeriod","LHCLuminosity","BeamIntensity","L3Current","L3Polarity","DipoleCurrent","DipolePolarity","CavernTemperature","CavernAtmosPressure","gva_cr5AtmosphericPressure","gva_meyrinAtmosphericPressure"};
  53
  54 //_______________________________________________________________
  55 AliGRPPreprocessor::AliGRPPreprocessor():
  56   AliPreprocessor("GRP",0), fPressure(0) {
  57   // default constructor - Don't use this!
  58
  59 }
  60
  61 //_______________________________________________________________
  62 AliGRPPreprocessor::AliGRPPreprocessor(AliShuttleInterface* shuttle):
  63   AliPreprocessor("GRP",shuttle), fPressure(0) {
  64   // constructor - shuttle must be instantiated!
  65
  66 }
  67
  68 //_______________________________________________________________
  69 AliGRPPreprocessor::~AliGRPPreprocessor() {
  70   //destructor
  71   delete fPressure;
  72 }
  73
  74 //_______________________________________________________________
  75 void AliGRPPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) {
  76   // Initialize preprocessor
  77
  78   AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp(startTime).AsString(), TTimeStamp(endTime).AsString()));
  79
  80   fRun = run;
  81   fStartTime = startTime;
  82   fEndTime = endTime;
  83   AliInfo("Initialization of the GRP preprocessor.");
  84
  85   TClonesArray * array = new TClonesArray("AliDCSSensor",2);
  86   for(Int_t j = 0; j < 2; j++) {
  87     AliDCSSensor * sens = new ((*array)[j])AliDCSSensor;
  88     sens->SetStringID(fgkDCSDataPoints[j+10]);
  89   }
  90   AliInfo(Form("Pressure Entries: %d",array->GetEntries()));
  91
  92   fPressure = new AliDCSSensorArray(fStartTime, fEndTime, array);
  93 }
  94
  95 //_______________________________________________________________
  96 UInt_t AliGRPPreprocessor::Process(TMap* valueMap) {
  97   // process data retrieved by the Shuttle
  98
  99   //=================//
 100   // DAQ logbook     //
 101   //=================//
 102
 103   TList *daqlblist = ProcessDaqLB();
 104   if(!daqlblist) {
 105     Log(Form("Problem with the DAQ logbook parameters!!!"));
 106     return 1;
 107   }
 108   TMap *m1 = (TMap *)daqlblist->At(0);
 109   TObjString *s1 = (TObjString *)m1->GetValue("fAliceStartTime");
 110   UInt_t iStartTime = atoi(s1->String().Data());
 111   TMap *m2 = (TMap *)daqlblist->At(1);
 112   TObjString *s2 = (TObjString *)m2->GetValue("fAliceStopTime");
 113   UInt_t iStopTime = atoi(s2->String().Data());
 114
 115   //=================//
 116   // DAQ FXS         //
 117   //=================//
 118   UInt_t iDaqFxs = ProcessDaqFxs();
 119   if(iDaqFxs == 0) {
 120         Log(Form("ProcessDaqFxs successful!"));
 121   } else {
 122         Log(Form("Could not store run raw tag file!"));
 123         return 1;
 124   }
 125
 126   //=================//
 127   // DCS FXS         //
 128   //=================//
 129   UInt_t iDcsFxs = ProcessDcsFxs();
 130   if(iDcsFxs == 0) {
 131         Log(Form("ProcessDcsFxs successful!"));
 132   } else {
 133         Log(Form("Could not store CTP run configuration and scalers!"));
 134         return 1;
 135   }
 136
 137   //=================//
 138   // DCS data points //
 139   //=================//
 140   TList *dcsdplist = ProcessDcsDPs(valueMap, iStartTime, iStopTime);
 141   if(!dcsdplist) {
 142     Log(Form("Problem with the DCS data points!!!"));
 143     return 1;
 144   }
 145   if(dcsdplist->GetEntries() != 10) {
 146     Log(Form("Problem with the DCS data points!!!"));
 147     // return 1; // TODO:COMMENTED FOR TESTING PURPOSES!
 148   }
 149   //NEEDS TO BE REVISED - BREAKS!!!
 150 //   AliDCSSensorArray *dcsSensorArray = GetPressureMap(valueMap,fPressure);
 151 //   if(!dcsSensorArray) {
 152 //     Log(Form("Problem with the pressure sensor values!!!"));
 153 //     return 0;
 154 //   }
 155
 156   TList * list = new TList();
 157   list = GetGlobalList(daqlblist,dcsdplist);
 158   list->SetOwner(1);
 159   AliInfo(Form("Final list entries: %d",list->GetEntries()));
 160
 161   AliCDBMetaData md;
 162   md.SetResponsible("Panos Christakoglou");
 163   md.SetComment("Output parameters from the GRP preprocessor.");
 164
 165   Bool_t result = Store("GRP", "Data", list, &md);
 166
 167   delete list;
 168
 169   if (result)
 170     return 0;
 171   else
 172     return 1;
 173 }
 174
 175 //_______________________________________________________________
 176 TList *AliGRPPreprocessor::GetGlobalList(TList *l1, TList *l2) {
 177   //Getting the global output TList
 178   TList *list = new TList();
 179   TMap *map = new TMap();
 180   for(Int_t i = 0; i < l1->GetEntries(); i++)
 181     list->AddLast(map = (TMap *)l1->At(i));
 182   for(Int_t i = 0; i < l2->GetEntries(); i++)
 183     list->AddLast(map = (TMap *)l2->At(i));
 184
 185   return list;
 186 }
 187
 188 //_______________________________________________________________
 189 TList *AliGRPPreprocessor::ProcessDaqLB() {
 190   //Getting the DAQ lb informnation
 191   const char* timeStart = GetRunParameter("time_start");
 192   const char* timeEnd = GetRunParameter("time_end");
 193   const char* beamEnergy = GetRunParameter("beamEnergy");
 194   const char* beamType = GetRunParameter("beamType");
 195   const char* numberOfDetectors = GetRunParameter("numberOfDetectors");
 196   const char* detectorMask = GetRunParameter("detectorMask");
 197   const char* lhcPeriod = GetRunParameter("LHCperiod");
 198
 199   if (timeStart) {
 200     Log(Form("Start time for run %d: %s",fRun, timeStart));
 201   } else {
 202     Log(Form("Start time not put in logbook!"));
 203   }
 204   TMap *mapDAQ1 = new TMap();
 205   mapDAQ1->Add(new TObjString("fAliceStartTime"),new TObjString(timeStart));
 206
 207   if (timeEnd) {
 208     Log(Form("End time for run %d: %s",fRun, timeEnd));
 209   } else {
 210     Log(Form("End time not put in logbook!"));
 211   }
 212   TMap *mapDAQ2 = new TMap();
 213   mapDAQ2->Add(new TObjString("fAliceStopTime"),new TObjString(timeEnd));
 214
 215   if (beamEnergy) {
 216     Log(Form("Beam energy for run %d: %s",fRun, beamEnergy));
 217   } else {
 218     Log(Form("Beam energy not put in logbook!"));
 219   }
 220   TMap *mapDAQ3 = new TMap();
 221   mapDAQ3->Add(new TObjString("fAliceBeamEnergy"),new TObjString(beamEnergy));
 222
 223   if (beamType) {
 224     Log(Form("Beam type for run %d: %s",fRun, beamType));
 225   } else {
 226     Log(Form("Beam type not put in logbook!"));
 227   }
 228   TMap *mapDAQ4 = new TMap();
 229   mapDAQ4->Add(new TObjString("fAliceBeamType"),new TObjString(beamType));
 230
 231   if (numberOfDetectors) {
 232     Log(Form("Number of active detectors for run %d: %s",fRun, numberOfDetectors));
 233   } else {
 234     Log(Form("Number of active detectors not put in logbook!"));
 235   }
 236   TMap *mapDAQ5 = new TMap();
 237   mapDAQ5->Add(new TObjString("fNumberOfDetectors"),new TObjString(numberOfDetectors));
 238
 239   if (detectorMask) {
 240     Log(Form("Detector mask for run %d: %s",fRun, detectorMask));
 241   } else {
 242     Log(Form("Detector mask not put in logbook!"));
 243   }
 244   TMap *mapDAQ6 = new TMap();
 245   mapDAQ6->Add(new TObjString("fDetectorMask"),new TObjString(detectorMask));
 246
 247   if (lhcPeriod) {
 248     Log(Form("LHC period (DAQ) for run %d: %s",fRun, lhcPeriod));
 249   } else {
 250     Log(Form("LHCperiod not put in logbook!"));
 251   }
 252   TMap *mapDAQ7 = new TMap();
 253   mapDAQ7->Add(new TObjString("fLHCPeriod"),new TObjString(lhcPeriod));
 254
 255   TList *list = new TList();
 256   list->Add(mapDAQ1); list->Add(mapDAQ2);
 257   list->Add(mapDAQ3); list->Add(mapDAQ4);
 258   list->Add(mapDAQ5); list->Add(mapDAQ6);
 259   list->Add(mapDAQ7);
 260
 261   return list;
 262 }
 263
 264 //_______________________________________________________________
 265 UInt_t AliGRPPreprocessor::ProcessDaqFxs() {
 266   //======DAQ FXS======//
 267
 268   TList* list = GetFileSources(kDAQ);
 269   if (!list) {
 270     Log("No raw data tag list: connection problems with DAQ FXS logbook!");
 271     return 1;
 272   }
 273
 274   if (list->GetEntries() == 0)
 275   {
 276         Log("no raw data tags in this run: nothing to merge!");
 277         delete  list; list=0;
 278         return 0;
 279   }
 280
 281   TChain *fRawTagChain = new TChain("T");
 282   Int_t nFiles=0;
 283   TIterator* iter = list->MakeIterator();
 284   TObject* obj = 0;
 285   while ((obj = iter->Next())) {
 286     TObjString* objStr = dynamic_cast<TObjString*> (obj);
 287     if (objStr) {
 288       Log(Form("Found source %s", objStr->String().Data()));
 289       TList* list2 = GetFileIDs(kDAQ, objStr->String());
 290       if (!list2) {
 291         Log("No list with ids from DAQ was found: connection problems with DAQ FXS logbook!");
 292         delete fRawTagChain; fRawTagChain=0;
 293         return 1;
 294       }
 295       Log(Form("Number of ids: %d",list2->GetEntries()));
 296       for(Int_t i = 0; i < list2->GetEntries(); i++) {
 297         TObjString *idStr = (TObjString *)list2->At(i);
 298         TString fileName = GetFile(kDAQ,idStr->String().Data(),objStr->String().Data());
 299         if (fileName.Length() > 0)
 300         {
 301                 Log(Form("Adding file in the chain: %s",fileName.Data()));
 302                 fRawTagChain->Add(fileName.Data());
 303                 nFiles++;
 304         } else {
 305                 Log(Form("Could not retrieve file with id %s from source %s: "
 306                         "connection problems with DAQ FXS!",
 307                                 idStr->String().Data(),objStr->String().Data()));
 308                 delete list; list=0;
 309                 delete list2; list2=0;
 310                 delete fRawTagChain; fRawTagChain=0;
 311                 return 2;
 312         }
 313       }
 314       delete list2;
 315     }
 316   }
 317
 318   TString fRawDataFileName = "GRP_Merged.tag.root";
 319   Log(Form("Merging %d raw data tags into file: %s", nFiles, fRawDataFileName.Data()));
 320   fRawTagChain->Merge(fRawDataFileName);
 321
 322   TString outputfile = Form("Run%d.Merged.RAW.tag.root", fRun);
 323   Bool_t result = StoreRunMetadataFile(fRawDataFileName.Data(),outputfile.Data());
 324
 325   if (!result)
 326   {
 327         Log("Problem storing raw data tags in local file!!");
 328   } else {
 329         Log("Raw data tags merged successfully!!");
 330   }
 331
 332   delete iter;
 333   delete list;
 334   delete fRawTagChain; fRawTagChain=0;
 335
 336   if (result == kFALSE)
 337   {
 338         return 3;
 339   }
 340
 341   return 0;
 342
 343 }
 344
 345 //_______________________________________________________________
 346 UInt_t AliGRPPreprocessor::ProcessDcsFxs() {
 347   //======DCS FXS======//
 348   // Get the CTP run configuration
 349   // and scalers from DCS FXS
 350
 351   {
 352     // Get the CTP run configuration
 353     TList* list = GetFileSources(kDCS,"CTP_runconfig");
 354     if (!list) {
 355       Log("No CTP runconfig file: connection problems with DCS FXS logbook!");
 356       return 1;
 357     }
 358
 359     if (list->GetEntries() == 0) {
 360       Log("No CTP runconfig file to be processed!");
 361     }
 362     else {
 363       TIter iter(list);
 364       TObjString *source;
 365       while ((source = dynamic_cast<TObjString *> (iter.Next()))) {
 366         TString runcfgfile = GetFile(kDCS, "CTP_runconfig", source->GetName());
 367         if (runcfgfile.IsNull()) {
 368           Log("No CTP runconfig files has been found: empty source!");
 369         }
 370         else {
 371           Log(Form("File with Id CTP_runconfig found in source %s! Copied to %s",source->GetName(),runcfgfile.Data()));
 372           AliTriggerConfiguration *runcfg = AliTriggerConfiguration::LoadConfiguration(runcfgfile);
 373           if (!runcfg) {
 374             Log("Bad CTP run configuration file! The corresponding CDB entry will not be filled!");
 375           }
 376           else {
 377             AliCDBMetaData metaData;
 378             metaData.SetBeamPeriod(0);
 379             metaData.SetResponsible("Roman Lietava");
 380             metaData.SetComment("CTP run configuration");
 381             if (!Store("CTP","Config", runcfg, &metaData, 0, 0)) {
 382               Log("Unable to store the CTP run configuration object to OCDB!");
 383             }
 384           }
 385         }
 386       }
 387     }
 388     delete list;
 389   }
 390
 391   {
 392     // Get the CTP counters information
 393     TList* list = GetFileSources(kDCS,"CTP_xcounters");
 394     if (!list) {
 395       Log("No CTP counters file: connection problems with DAQ FXS logbook!");
 396       return 1;
 397     }
 398
 399     if (list->GetEntries() == 0) {
 400       Log("No CTP counters file to be processed!");
 401     }
 402     else {
 403       TIter iter(list);
 404       TObjString *source;
 405       while ((source = dynamic_cast<TObjString *> (iter.Next()))) {
 406         TString countersfile = GetFile(kDCS, "CTP_xcounters", source->GetName());
 407         if (countersfile.IsNull()) {
 408           Log("No CTP counters files has been found: empty source!");
 409         }
 410         else {
 411           Log(Form("File with Id CTP_xcounters found in source %s! Copied to %s",source->GetName(),countersfile.Data()));
 412           AliTriggerRunScalers *scalers = AliTriggerRunScalers::ReadScalers(countersfile);
 413           if (!scalers) {
 414             Log("Bad CTP counters file! The corresponding CDB entry will not be filled!");
 415           }
 416           else {
 417             AliCDBMetaData metaData;
 418             metaData.SetBeamPeriod(0);
 419             metaData.SetResponsible("Roman Lietava");
 420             metaData.SetComment("CTP scalers");
 421             if (!Store("CTP","Scalers", scalers, &metaData, 0, 0)) {
 422               Log("Unable to store the CTP scalers object to OCDB!");
 423             }
 424           }
 425         }
 426       }
 427     }
 428     delete list;
 429   }
 430
 431   return 0;
 432 }
 433
 434 //_______________________________________________________________
 435 TList *AliGRPPreprocessor::ProcessDcsDPs(TMap* valueMap, UInt_t iStartTime, UInt_t iStopTime) {
 436   //Getting the DCS dps
 437   //===========//
 438
 439   TList *list = new TList();
 440
 441   //DCS data points
 442   //===========//
 443   AliInfo(Form("==========LHCState==========="));
 444   TObjArray *aliasLHCState = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[0]);
 445   if(!aliasLHCState) {
 446     Log(Form("LHCState not found!!!"));
 447     return list;
 448   }
 449   AliGRPDCS *dcs1 = new AliGRPDCS(aliasLHCState,iStartTime,iStopTime);
 450   TString sLHCState = dcs1->ProcessDCS(3);
 451   if (sLHCState) {
 452     Log(Form("<LHCState> for run %d: %s",fRun, sLHCState.Data()));
 453   } else {
 454     Log(Form("LHCState not put in TMap!"));
 455   }
 456   TMap *mapDCS1 = new TMap();
 457   mapDCS1->Add(new TObjString("fLHCState"),new TObjString(sLHCState));
 458   list->Add(mapDCS1);
 459
 460   AliInfo(Form("==========LHCPeriod==========="));
 461   TObjArray *aliasLHCPeriod = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[1]);
 462   if(!aliasLHCPeriod) {
 463     Log(Form("LHCPeriod not found!!!"));
 464     return list;
 465   }
 466   AliGRPDCS *dcs2 = new AliGRPDCS(aliasLHCPeriod,iStartTime,iStopTime);
 467   TString sLHCPeriod = dcs2->ProcessDCS(3);
 468   if (sLHCPeriod) {
 469     Log(Form("<LHCPeriod> for run %d: %s",fRun, sLHCPeriod.Data()));
 470   } else {
 471     Log(Form("LHCPeriod not put in TMap!"));
 472   }
 473   TMap *mapDCS2 = new TMap();
 474   mapDCS2->Add(new TObjString("fLHCCondition"),new TObjString(sLHCPeriod));
 475   list->Add(mapDCS2);
 476
 477   AliInfo(Form("==========LHCLuminosity==========="));
 478   TObjArray *aliasLHCLuminosity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[2]);
 479   if(!aliasLHCLuminosity) {
 480     Log(Form("LHCLuminosity not found!!!"));
 481     return list;
 482   }
 483   AliGRPDCS *dcs3 = new AliGRPDCS(aliasLHCLuminosity,iStartTime,iStopTime);
 484   TString sMeanLHCLuminosity = dcs3->ProcessDCS(2);
 485   if (sMeanLHCLuminosity) {
 486     Log(Form("<LHCLuminosity> for run %d: %s",fRun, sMeanLHCLuminosity.Data()));
 487   } else {
 488     Log(Form("LHCLuminosity not put in TMap!"));
 489   }
 490   TMap *mapDCS3 = new TMap();
 491   mapDCS3->Add(new TObjString("fLHCLuminosity"),new TObjString(sMeanLHCLuminosity));
 492   list->Add(mapDCS3);
 493
 494   AliInfo(Form("==========BeamIntensity==========="));
 495   TObjArray *aliasBeamIntensity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[3]);
 496   if(!aliasBeamIntensity) {
 497     Log(Form("BeamIntensity not found!!!"));
 498     return list;
 499   }
 500   AliGRPDCS *dcs4 = new AliGRPDCS(aliasBeamIntensity,iStartTime,iStopTime);
 501   TString sMeanBeamIntensity = dcs4->ProcessDCS(2);
 502   if (sMeanBeamIntensity) {
 503     Log(Form("<BeamIntensity> for run %d: %s",fRun, sMeanBeamIntensity.Data()));
 504   } else {
 505     Log(Form("BeamIntensity not put in TMap!"));
 506   }
 507   TMap *mapDCS4 = new TMap();
 508   mapDCS4->Add(new TObjString("fBeamIntensity"),new TObjString(sMeanBeamIntensity));
 509   list->Add(mapDCS4);
 510
 511   AliInfo(Form("==========L3Current==========="));
 512   TObjArray *aliasL3Current = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[4]);
 513   if(!aliasL3Current) {
 514     Log(Form("L3Current not found!!!"));
 515     return list;
 516   }
 517   AliGRPDCS *dcs5 = new AliGRPDCS(aliasL3Current,iStartTime,iStopTime);
 518   TString sMeanL3Current = dcs5->ProcessDCS(2);
 519   if (sMeanL3Current) {
 520     Log(Form("<L3Current> for run %d: %s",fRun, sMeanL3Current.Data()));
 521   } else {
 522     Log(Form("L3Current not put in TMap!"));
 523   }
 524   TMap *mapDCS5 = new TMap();
 525   mapDCS5->Add(new TObjString("fL3Current"),new TObjString(sMeanL3Current));
 526   list->Add(mapDCS5);
 527
 528   AliInfo(Form("==========L3Polarity==========="));
 529   TObjArray *aliasL3Polarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[5]);
 530   if(!aliasL3Polarity) {
 531     Log(Form("L3Polarity not found!!!"));
 532     return list;
 533   }
 534   AliGRPDCS *dcs6 = new AliGRPDCS(aliasL3Polarity,iStartTime,iStopTime);
 535   TString sL3Polarity = dcs6->ProcessDCS(4);
 536   if (sL3Polarity) {
 537     Log(Form("<L3Polarity> for run %d: %s",fRun, sL3Polarity.Data()));
 538   } else {
 539     Log(Form("L3Polarity not put in TMap!"));
 540   }
 541   TMap *mapDCS6 = new TMap();
 542   mapDCS6->Add(new TObjString("fL3Polarity"),new TObjString(sL3Polarity));
 543   list->Add(mapDCS6);
 544
 545   AliInfo(Form("==========DipoleCurrent==========="));
 546   TObjArray *aliasDipoleCurrent = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[6]);
 547   if(!aliasDipoleCurrent) {
 548     Log(Form("DipoleCurrent not found!!!"));
 549     return list;
 550   }
 551   AliGRPDCS *dcs7 = new AliGRPDCS(aliasDipoleCurrent,iStartTime,iStopTime);
 552   TString sMeanDipoleCurrent = dcs7->ProcessDCS(2);
 553   if (sMeanDipoleCurrent) {
 554     Log(Form("<DipoleCurrent> for run %d: %s",fRun, sMeanDipoleCurrent.Data()));
 555   } else {
 556     Log(Form("DipoleCurrent not put in TMap!"));
 557   }
 558   TMap *mapDCS7 = new TMap();
 559   mapDCS7->Add(new TObjString("fDipoleCurrent"),new TObjString(sMeanDipoleCurrent));
 560   list->Add(mapDCS7);
 561
 562   AliInfo(Form("==========DipolePolarity==========="));
 563   TObjArray *aliasDipolePolarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[7]);
 564   if(!aliasDipolePolarity) {
 565     Log(Form("DipolePolarity not found!!!"));
 566     return list;
 567   }
 568   AliGRPDCS *dcs8 = new AliGRPDCS(aliasDipolePolarity,iStartTime,iStopTime);
 569   TString sDipolePolarity = dcs8->ProcessDCS(4);
 570   if (sDipolePolarity) {
 571     Log(Form("<DipolePolarity> for run %d: %s",fRun, sDipolePolarity.Data()));
 572   } else {
 573     Log(Form("DipolePolarity not put in TMap!"));
 574   }
 575   TMap *mapDCS8 = new TMap();
 576   mapDCS8->Add(new TObjString("fDipolePolarity"),new TObjString(sDipolePolarity));
 577   list->Add(mapDCS8);
 578
 579   AliInfo(Form("==========CavernTemperature==========="));
 580   TObjArray *aliasCavernTemperature = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[8]);
 581   if(!aliasCavernTemperature) {
 582     Log(Form("CavernTemperature not found!!!"));
 583     return list;
 584   }
 585   AliGRPDCS *dcs9 = new AliGRPDCS(aliasCavernTemperature,iStartTime,iStopTime);
 586   TString sMeanCavernTemperature = dcs9->ProcessDCS(2);
 587   if (sMeanCavernTemperature) {
 588     Log(Form("<CavernTemperature> for run %d: %s",fRun, sMeanCavernTemperature.Data()));
 589   } else {
 590     Log(Form("CavernTemperature not put in TMap!"));
 591   }
 592   TMap *mapDCS9 = new TMap();
 593   mapDCS9->Add(new TObjString("fCavernTemperature"),new TObjString(sMeanCavernTemperature));
 594   list->Add(mapDCS9);
 595
 596   AliInfo(Form("==========CavernPressure==========="));
 597   TObjArray *aliasCavernPressure = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[9]);
 598   if(!aliasCavernPressure) {
 599     Log(Form("CavernPressure not found!!!"));
 600     return list;
 601   }
 602   AliGRPDCS *dcs10 = new AliGRPDCS(aliasCavernPressure,iStartTime,iStopTime);
 603   TString sMeanCavernPressure = dcs10->ProcessDCS(2);
 604   if (sMeanCavernPressure) {
 605     Log(Form("<CavernPressure> for run %d: %s",fRun, sMeanCavernPressure.Data()));
 606   } else {
 607     Log(Form("CavernPressure not put in TMap!"));
 608   }
 609   TMap *mapDCS10 = new TMap();
 610   mapDCS10->Add(new TObjString("fCavernPressure"),new TObjString(sMeanCavernPressure));
 611   list->Add(mapDCS10);
 612
 613   return list;
 614 }
 615
 616 //_______________________________________________________________
 617 AliDCSSensorArray *AliGRPPreprocessor::GetPressureMap(TMap* dcsAliasMap, AliDCSSensorArray *fPressure) {
 618   // extract DCS pressure maps. Perform fits to save space
 619
 620   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
 621   if (map) {
 622     fPressure->MakeSplineFit(map);
 623     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors();
 624     if (fitFraction > kFitFraction ) {
 625       AliInfo(Form("Pressure values extracted, fits performed.\n"));
 626     } else {
 627       AliInfo("Too few pressure maps fitted. \n");
 628     }
 629   } else {
 630     AliInfo("AliGRPDCS: no atmospheric pressure map extracted. \n");
 631   }
 632   delete map;
 633
 634   return fPressure;
 635 }
 636
 637 //_______________________________________________________________
 638 /*UInt_t AliGRPPreprocessor::MapPressure(TMap* dcsAliasMap) {
 639   // extract DCS pressure maps. Perform fits to save space
 640
 641   UInt_t result=0;
 642   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
 643   if (map) {
 644     fPressure->MakeSplineFit(map);
 645     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors();
 646     if (fitFraction > kFitFraction ) {
 647       AliInfo(Form("Pressure values extracted, fits performed.\n"));
 648     } else {
 649       Log ("Too few pressure maps fitted. \n");
 650       result = 9;
 651     }
 652   } else {
 653     Log("AliTPCPreprocsessor: no atmospheric pressure map extracted. \n");
 654     result=9;
 655   }
 656   delete map;
 657   // Now store the final CDB file
 658
 659   if ( result == 0 ) {
 660     AliCDBMetaData metaData;
 661     metaData.SetBeamPeriod(0);
 662     metaData.SetResponsible("Panos Christakoglou");
 663     metaData.SetComment("Preprocessor AliGRP data base pressure entries.");
 664
 665     Bool_t storeOK = Store("Calib", "Pressure", fPressure, &metaData, 0, 0);
 666     if ( !storeOK ) result=1;
 667   }
 668
 669   return result;
 670   }*/