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