[u/mrichter/AliRoot.git] / STEER / AliGRPPreprocessor.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

//-------------------------------------------------------------------------
//                          Class AliGRPPreprocessor
//                  Global Run Parameters (GRP) preprocessor
//    Origin: Panos Christakoglou, UOA-CERN, Panos.Christakoglou@cern.ch
//-------------------------------------------------------------------------

#include <TChain.h>
#include <TList.h>
#include <TMap.h>
#include <TObjString.h>
#include <TTimeStamp.h>

#include "AliGRPPreprocessor.h"
#include "AliGRPDCS.h"
#include "AliDCSSensorArray.h"

#include "AliCDBMetaData.h"
#include "AliLog.h"

class AliDCSValue;
class AliShuttleInterface;

#include <TH1.h>

const Double_t kFitFraction = 0.7;                 // Fraction of DCS sensor fits required

ClassImp(AliGRPPreprocessor)

//_______________________________________________________________
  const char* AliGRPPreprocessor::fgkDCSDataPoints[12] = {"LHCState","LHCPeriod","LHCLuminosity","BeamIntensity","L3Current","L3Polarity","DipoleCurrent","DipolePolarity","CavernTemperature","CavernAtmosPressure","gva_cr5AtmosphericPressure","gva_meyrinAtmosphericPressure"};

//_______________________________________________________________
AliGRPPreprocessor::AliGRPPreprocessor():
  AliPreprocessor("GRP",0), fPressure(0) {
  // default constructor - Don't use this!
  
}

//_______________________________________________________________
AliGRPPreprocessor::AliGRPPreprocessor(AliShuttleInterface* shuttle):
  AliPreprocessor("GRP",shuttle), fPressure(0) {
  // constructor - shuttle must be instantiated!
  
}

//_______________________________________________________________
AliGRPPreprocessor::~AliGRPPreprocessor() {
  //destructor
  delete fPressure;
}

//_______________________________________________________________
void AliGRPPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) {
  // Initialize preprocessor
  
  AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp(startTime).AsString(), TTimeStamp(endTime).AsString()));
  
  fRun = run;
  fStartTime = startTime;
  fEndTime = endTime;
  AliInfo("Initialization of the GRP preprocessor.");

  TClonesArray * array = new TClonesArray("AliDCSSensor",2);
  for(Int_t j = 0; j < 2; j++) {
    AliDCSSensor * sens = new ((*array)[j])AliDCSSensor;
    sens->SetStringID(fgkDCSDataPoints[j+10]);
  }
  AliInfo(Form("Pressure Entries: %d",array->GetEntries()));

  fPressure = new AliDCSSensorArray(fStartTime, fEndTime, array);
}

//_______________________________________________________________
UInt_t AliGRPPreprocessor::Process(TMap* valueMap) {
  // process data retrieved by the Shuttle
  
  //=================//
  // DAQ logbook     //
  //=================//
  TList *daqlblist = ProcessDaqLB();
  if(!daqlblist) {
    Log(Form("Problem with the DAQ logbook parameters!!!"));
    return 0;
  }
  TMap *m1 = (TMap *)daqlblist->At(0);
  TObjString *s1 = (TObjString *)m1->GetValue("fAliceStartTime");
  UInt_t iStartTime = atoi(s1->String().Data());
  TMap *m2 = (TMap *)daqlblist->At(1);
  TObjString *s2 = (TObjString *)m2->GetValue("fAliceStopTime");
  UInt_t iStopTime = atoi(s2->String().Data());
  TMap *m3 = (TMap *)daqlblist->At(6);
  TObjString *s3 = (TObjString *)m3->GetValue("fLHCPeriod");
  TString productionYear = "";

  //=================//
  // DAQ FXS         //
  //=================//
  UInt_t iDaqFxs = ProcessDaqFxs(s3->String(),productionYear);
  if(iDaqFxs == 0) Log(Form("Raw data merged tags copied succesfully in AliEn!!!"));
 
  //=================//
  // DCS data points //
  //=================//
  TList *dcsdplist = ProcessDcsDPs(valueMap, iStartTime, iStopTime);
  if(!dcsdplist) {
    Log(Form("Problem with the DCS data points!!!"));
    return 0;
  }    
  if(dcsdplist->GetEntries() != 10) {
    Log(Form("Problem with the DCS data points!!!"));
    return 0;
  }
  //NEEDS TO BE REVISED - BREAKS!!!
//   AliDCSSensorArray *dcsSensorArray = GetPressureMap(valueMap,fPressure);
//   if(!dcsSensorArray) {
//     Log(Form("Problem with the pressure sensor values!!!"));
//     return 0;
//   }

  TList * list = new TList();
  list = GetGlobalList(daqlblist,dcsdplist);
  list->SetOwner(1);
  AliInfo(Form("Final list entries: %d",list->GetEntries()));
  
  AliCDBMetaData md;
  md.SetResponsible("Panos Christakoglou");
  md.SetComment("Output parameters from the GRP preprocessor.");
  
  Bool_t result = Store("GRP", "Data", list, &md);
  
  delete list;
  
  if (result)
    return 0;
  else
    return 1;
}

//_______________________________________________________________
TList *AliGRPPreprocessor::GetGlobalList(TList *l1, TList *l2) {
  //Getting the global output TList
  TList *list = new TList();
  TMap *map = new TMap();
  for(Int_t i = 0; i < l1->GetEntries(); i++) 
    list->AddLast(map = (TMap *)l1->At(i));
  for(Int_t i = 0; i < l2->GetEntries(); i++) 
    list->AddLast(map = (TMap *)l2->At(i));

  return list;
}

//_______________________________________________________________
TList *AliGRPPreprocessor::ProcessDaqLB() {
  //Getting the DAQ lb informnation
  const char* timeStart = GetRunParameter("time_start");
  const char* timeEnd = GetRunParameter("time_end");
  const char* beamEnergy = GetRunParameter("beamEnergy");
  const char* beamType = GetRunParameter("beamType");
  const char* numberOfDetectors = GetRunParameter("numberOfDetectors");
  const char* detectorMask = GetRunParameter("detectorMask");
  const char* lhcPeriod = GetRunParameter("LHCperiod");

  if (timeStart) {
    Log(Form("Start time for run %d: %s",fRun, timeStart));
  } else {
    Log(Form("Start time not put in logbook!"));
  }
  TMap *mapDAQ1 = new TMap();
  mapDAQ1->Add(new TObjString("fAliceStartTime"),new TObjString(timeStart));

  if (timeEnd) {
    Log(Form("End time for run %d: %s",fRun, timeEnd));
  } else {
    Log(Form("End time not put in logbook!"));
  }
  TMap *mapDAQ2 = new TMap();
  mapDAQ2->Add(new TObjString("fAliceStopTime"),new TObjString(timeEnd));

  if (beamEnergy) {
    Log(Form("Beam energy for run %d: %s",fRun, beamEnergy));
  } else {
    Log(Form("Beam energy not put in logbook!"));
  }
  TMap *mapDAQ3 = new TMap();
  mapDAQ3->Add(new TObjString("fAliceBeamEnergy"),new TObjString(beamEnergy));

  if (beamType) {
    Log(Form("Beam type for run %d: %s",fRun, beamType));
  } else {
    Log(Form("Beam type not put in logbook!"));
  }
  TMap *mapDAQ4 = new TMap();
  mapDAQ4->Add(new TObjString("fAliceBeamType"),new TObjString(beamType));

  if (numberOfDetectors) {
    Log(Form("Number of active detectors for run %d: %s",fRun, numberOfDetectors));
  } else {
    Log(Form("Number of active detectors not put in logbook!"));
  }
  TMap *mapDAQ5 = new TMap();
  mapDAQ5->Add(new TObjString("fNumberOfDetectors"),new TObjString(numberOfDetectors));

  if (detectorMask) {
    Log(Form("Detector mask for run %d: %s",fRun, detectorMask));
  } else {
    Log(Form("Detector mask not put in logbook!"));
  }
  TMap *mapDAQ6 = new TMap();
  mapDAQ6->Add(new TObjString("fDetectorMask"),new TObjString(detectorMask));

  if (lhcPeriod) {
    Log(Form("LHC period (DAQ) for run %d: %s",fRun, lhcPeriod));
  } else {
    Log(Form("LHCperiod not put in logbook!"));
  }
  TMap *mapDAQ7 = new TMap();
  mapDAQ7->Add(new TObjString("fLHCPeriod"),new TObjString(lhcPeriod));

  TList *list = new TList();
  list->Add(mapDAQ1); list->Add(mapDAQ2);
  list->Add(mapDAQ3); list->Add(mapDAQ4);
  list->Add(mapDAQ5); list->Add(mapDAQ6);
  list->Add(mapDAQ7);

  return list;
}

//_______________________________________________________________
UInt_t AliGRPPreprocessor::ProcessDaqFxs(TString lhcperiod, TString productionYear) {
  //======DAQ FXS======//
  TChain *fRawTagChain = new TChain("T");
  TString fRawDataFileName;
  TList* list = GetFileSources(kDAQ);  
  if (!list) {
    Log("No raw data tag list found!!!");
    return 1;
  }
  TIterator* iter = list->MakeIterator();
  TObject* obj = 0;
  while ((obj = iter->Next())) {
    TObjString* objStr = dynamic_cast<TObjString*> (obj);
    if (objStr) {
      Log(Form("Found source %s", objStr->String().Data()));
      TList* list2 = GetFileIDs(kDAQ, objStr->String());
      if (!list2) {
	Log("No list with ids from DAQ was found!!!");
	return 2;
      }
      Log(Form("Number of ids: %d",list2->GetEntries()));
      for(Int_t i = 0; i < list2->GetEntries(); i++) {
	TObjString *idStr = (TObjString *)list2->At(i);
	//Log(Form("Filename1: %s",idStr->String().Data()));
	TString fileName = GetFile(kDAQ,idStr->String().Data(),objStr->String().Data());      
	Log(Form("Adding file in the chain: %s",fileName.Data()));
	fRawTagChain->Add(fileName.Data());
	fRawDataFileName = fileName(0,fileName.First("_"));
      }
      delete list2;
    }
  }
  delete iter;
  delete list;
  fRawDataFileName += "_GRP_Merged.tag.root";
  Log(Form("Merging raw data tags into file: %s",fRawDataFileName.Data()));

  TString outputfile = "alien:///alice/data/"; 
  outputfile += productionYear.Data(); outputfile += "/";
  outputfile += lhcperiod.Data(); outputfile += "/";
  outputfile += fRun; outputfile += "/raw/"; 
  //StoreTagFiles(fRawDataFileName.Data(),outputfile.Data());

  return 0;
}

//_______________________________________________________________
TList *AliGRPPreprocessor::ProcessDcsDPs(TMap* valueMap, UInt_t iStartTime, UInt_t iStopTime) {
  //Getting the DCS dps
  //===========//
  
  TList *list = new TList();

  //DCS data points
  //===========//
  AliInfo(Form("==========LHCState==========="));
  TObjArray *aliasLHCState = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[0]);
  if(!aliasLHCState) {
    Log(Form("LHCState not found!!!"));
    return list;
  }
  AliGRPDCS *dcs1 = new AliGRPDCS(aliasLHCState,iStartTime,iStopTime);
  TString sLHCState = dcs1->ProcessDCS(3);  
  if (sLHCState) {
    Log(Form("<LHCState> for run %d: %s",fRun, sLHCState.Data()));
  } else {
    Log(Form("LHCState not put in TMap!"));
  }
  TMap *mapDCS1 = new TMap();
  mapDCS1->Add(new TObjString("fLHCState"),new TObjString(sLHCState));
  list->Add(mapDCS1);

  AliInfo(Form("==========LHCPeriod==========="));
  TObjArray *aliasLHCPeriod = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[1]);
  if(!aliasLHCPeriod) {
    Log(Form("LHCPeriod not found!!!"));
    return list;
  }
  AliGRPDCS *dcs2 = new AliGRPDCS(aliasLHCPeriod,iStartTime,iStopTime);
  TString sLHCPeriod = dcs2->ProcessDCS(3);  
  if (sLHCPeriod) {
    Log(Form("<LHCPeriod> for run %d: %s",fRun, sLHCPeriod.Data()));
  } else {
    Log(Form("LHCPeriod not put in TMap!"));
  }
  TMap *mapDCS2 = new TMap();
  mapDCS2->Add(new TObjString("fLHCCondition"),new TObjString(sLHCPeriod));
  list->Add(mapDCS2);

  AliInfo(Form("==========LHCLuminosity==========="));
  TObjArray *aliasLHCLuminosity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[2]);
  if(!aliasLHCLuminosity) {
    Log(Form("LHCLuminosity not found!!!"));
    return list;
  }
  AliGRPDCS *dcs3 = new AliGRPDCS(aliasLHCLuminosity,iStartTime,iStopTime);
  TString sMeanLHCLuminosity = dcs3->ProcessDCS(2);  
  if (sMeanLHCLuminosity) {
    Log(Form("<LHCLuminosity> for run %d: %s",fRun, sMeanLHCLuminosity.Data()));
  } else {
    Log(Form("LHCLuminosity not put in TMap!"));
  }
  TMap *mapDCS3 = new TMap();
  mapDCS3->Add(new TObjString("fLHCLuminosity"),new TObjString(sMeanLHCLuminosity));
  list->Add(mapDCS3);

  AliInfo(Form("==========BeamIntensity==========="));
  TObjArray *aliasBeamIntensity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[3]);
  if(!aliasBeamIntensity) {
    Log(Form("BeamIntensity not found!!!"));
    return list;
  }
  AliGRPDCS *dcs4 = new AliGRPDCS(aliasBeamIntensity,iStartTime,iStopTime);
  TString sMeanBeamIntensity = dcs4->ProcessDCS(2);  
  if (sMeanBeamIntensity) {
    Log(Form("<BeamIntensity> for run %d: %s",fRun, sMeanBeamIntensity.Data()));
  } else {
    Log(Form("BeamIntensity not put in TMap!"));
  }
  TMap *mapDCS4 = new TMap();
  mapDCS4->Add(new TObjString("fBeamIntensity"),new TObjString(sMeanBeamIntensity));
  list->Add(mapDCS4);

  AliInfo(Form("==========L3Current==========="));
  TObjArray *aliasL3Current = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[4]);
  if(!aliasL3Current) {
    Log(Form("L3Current not found!!!"));
    return list;
  }
  AliGRPDCS *dcs5 = new AliGRPDCS(aliasL3Current,iStartTime,iStopTime);
  TString sMeanL3Current = dcs5->ProcessDCS(2);  
  if (sMeanL3Current) {
    Log(Form("<L3Current> for run %d: %s",fRun, sMeanL3Current.Data()));
  } else {
    Log(Form("L3Current not put in TMap!"));
  }
  TMap *mapDCS5 = new TMap();
  mapDCS5->Add(new TObjString("fL3Current"),new TObjString(sMeanL3Current));
  list->Add(mapDCS5);

  AliInfo(Form("==========L3Polarity==========="));
  TObjArray *aliasL3Polarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[5]);
  if(!aliasL3Polarity) {
    Log(Form("L3Polarity not found!!!"));
    return list;
  }
  AliGRPDCS *dcs6 = new AliGRPDCS(aliasL3Polarity,iStartTime,iStopTime);
  TString sL3Polarity = dcs6->ProcessDCS(4);  
  if (sL3Polarity) {
    Log(Form("<L3Polarity> for run %d: %s",fRun, sL3Polarity.Data()));
  } else {
    Log(Form("L3Polarity not put in TMap!"));
  }
  TMap *mapDCS6 = new TMap();
  mapDCS6->Add(new TObjString("fL3Polarity"),new TObjString(sL3Polarity));
  list->Add(mapDCS6);

  AliInfo(Form("==========DipoleCurrent==========="));
  TObjArray *aliasDipoleCurrent = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[6]);
  if(!aliasDipoleCurrent) {
    Log(Form("DipoleCurrent not found!!!"));
    return list;
  }
  AliGRPDCS *dcs7 = new AliGRPDCS(aliasDipoleCurrent,iStartTime,iStopTime);
  TString sMeanDipoleCurrent = dcs7->ProcessDCS(2);  
  if (sMeanDipoleCurrent) {
    Log(Form("<DipoleCurrent> for run %d: %s",fRun, sMeanDipoleCurrent.Data()));
  } else {
    Log(Form("DipoleCurrent not put in TMap!"));
  }
  TMap *mapDCS7 = new TMap();
  mapDCS7->Add(new TObjString("fDipoleCurrent"),new TObjString(sMeanDipoleCurrent));
  list->Add(mapDCS7);

  AliInfo(Form("==========DipolePolarity==========="));
  TObjArray *aliasDipolePolarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[7]);
  if(!aliasDipolePolarity) {
    Log(Form("DipolePolarity not found!!!"));
    return list;
  }
  AliGRPDCS *dcs8 = new AliGRPDCS(aliasDipolePolarity,iStartTime,iStopTime);
  TString sDipolePolarity = dcs8->ProcessDCS(4);  
  if (sDipolePolarity) {
    Log(Form("<DipolePolarity> for run %d: %s",fRun, sDipolePolarity.Data()));
  } else {
    Log(Form("DipolePolarity not put in TMap!"));
  }
  TMap *mapDCS8 = new TMap();
  mapDCS8->Add(new TObjString("fDipolePolarity"),new TObjString(sDipolePolarity));
  list->Add(mapDCS8);

  AliInfo(Form("==========CavernTemperature==========="));
  TObjArray *aliasCavernTemperature = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[8]);
  if(!aliasCavernTemperature) {
    Log(Form("CavernTemperature not found!!!"));
    return list;
  }
  AliGRPDCS *dcs9 = new AliGRPDCS(aliasCavernTemperature,iStartTime,iStopTime);
  TString sMeanCavernTemperature = dcs9->ProcessDCS(2);  
  if (sMeanCavernTemperature) {
    Log(Form("<CavernTemperature> for run %d: %s",fRun, sMeanCavernTemperature.Data()));
  } else {
    Log(Form("CavernTemperature not put in TMap!"));
  }
  TMap *mapDCS9 = new TMap();
  mapDCS9->Add(new TObjString("fCavernTemperature"),new TObjString(sMeanCavernTemperature));
  list->Add(mapDCS9);

  AliInfo(Form("==========CavernPressure==========="));
  TObjArray *aliasCavernPressure = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[9]);
  if(!aliasCavernPressure) {
    Log(Form("CavernPressure not found!!!"));
    return list;
  }
  AliGRPDCS *dcs10 = new AliGRPDCS(aliasCavernPressure,iStartTime,iStopTime);
  TString sMeanCavernPressure = dcs10->ProcessDCS(2);  
  if (sMeanCavernPressure) {
    Log(Form("<CavernPressure> for run %d: %s",fRun, sMeanCavernPressure.Data()));
  } else {
    Log(Form("CavernPressure not put in TMap!"));
  }
  TMap *mapDCS10 = new TMap();
  mapDCS10->Add(new TObjString("fCavernPressure"),new TObjString(sMeanCavernPressure));
  list->Add(mapDCS10);

  return list;
}

//_______________________________________________________________
AliDCSSensorArray *AliGRPPreprocessor::GetPressureMap(TMap* dcsAliasMap, AliDCSSensorArray *fPressure) {
  // extract DCS pressure maps. Perform fits to save space
  
  TMap *map = fPressure->ExtractDCS(dcsAliasMap);
  if (map) {
    fPressure->MakeSplineFit(map);
    Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors(); 
    if (fitFraction > kFitFraction ) {
      AliInfo(Form("Pressure values extracted, fits performed.\n"));
    } else { 
      AliInfo("Too few pressure maps fitted. \n");
    }
  } else {
    AliInfo("AliGRPDCS: no atmospheric pressure map extracted. \n");
  }
  delete map;
 
  return fPressure;
}

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