Classes Preprocessor and DataDCS added

[u/mrichter/AliRoot.git] / VZERO / AliVZEROPreprocessor.cxx

#include "AliVZEROPreprocessor.h"
#include "AliVZEROCalibData.h"
#include "AliCDBMetaData.h"
#include "AliCDBEntry.h"
#include "AliDCSValue.h"
#include "AliLog.h"
#include <TFile.h>
#include <TTimeStamp.h>
#include <TObjString.h>
#include <TSystem.h>
#include <TList.h>

//
// This class is  a simple preprocessor for V0.
//

ClassImp(AliVZEROPreprocessor)

//______________________________________________________________________________________________
AliVZEROPreprocessor::AliVZEROPreprocessor(AliShuttleInterface* shuttle) :
  AliPreprocessor("V00", shuttle),
  fData(0)
 
{
  // constructor
}

//______________________________________________________________________________________________
AliVZEROPreprocessor::~AliVZEROPreprocessor()
{
  // destructor
}

//______________________________________________________________________________________________
void AliVZEROPreprocessor::Initialize(Int_t run, UInt_t startTime,
	UInt_t endTime)
{
  // Creates AliZDCDataDCS object

  AliPreprocessor::Initialize(run, startTime, endTime);

	Log(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
		TTimeStamp(startTime).AsString(),
		TTimeStamp(endTime).AsString()));

 	fRun = run;
        fStartTime = startTime;
        fEndTime = endTime;

	fData = new AliVZERODataDCS(fRun, fStartTime, fEndTime);
}



//______________________________________________________________________________________________
UInt_t AliVZEROPreprocessor::Process(TMap* dcsAliasMap)
{
  // Fills data retrieved from DCS and DAQ into a AliVZEROCalibData object and 
  // stores it into CalibrationDB


  // *** GET RUN TYPE ***
  TString runType = GetRunType();


  // *** REFERENCE DATA *** 
  
  TString fileName; 
  AliVZEROCalibData *calibData = new AliVZEROCalibData();
  
  // *************** From DCS ******************
  // Fills data into a AliVZERODataDCS object
  if(!dcsAliasMap) return 1;

 	// The processing of the DCS input data is forwarded to AliVZERODataDCS

	fData->ProcessData(*dcsAliasMap);
	//fData->Draw(""); 		// Draw the HV values as a function of time
	//dcsAliasMap->Print("");	// Print out the HV values

	// Writing VZERO PMs HV values into calibration object
	calibData->SetMeanHV(fData->GetMeanHV());
	calibData->SetWidthHV(fData->GetWidthHV());
  
  
   // *************** From DAQ ******************
	TString SourcesId[3];
	SourcesId[0] = "PEDESTALS";
	SourcesId[1] = "GAINS";
	SourcesId[2] = "SIGMAS";

	for(int iSource=0;iSource<3;iSource++){
  		TList* sourceList = GetFileSources(kDAQ, SourcesId[iSource].Data());
  		if (!sourceList){
			AliError(Form("No sources found for id %s", SourcesId[iSource].Data()));      		
      		return 1;   
		}
		AliInfo(Form("The following sources produced files with the id %s",SourcesId[iSource].Data()));
	        sourceList->Print();    

  		TIter iter(sourceList);
  		TObjString *source = 0;
		
		while((source=dynamic_cast<TObjString*> (iter.Next()))){
  			fileName = GetFile(kDAQ, SourcesId[iSource].Data(), source->GetName());
  		  	if (fileName.Length() > 0)
    			AliInfo(Form("Got the file %s, now we can extract some values.", fileName.Data()));
		  	FILE *file;
			if((file = fopen(fileName.Data(),"r")) == NULL){
            	AliError(Form("Cannot open file %s",fileName.Data()));
	    	  	return 1;
			}
		  	Float_t Values[128];
		  	for(Int_t j=0; j<128; j++)fscanf(file,"%f",&Values[j]);
			fclose(file);
	    	switch(iSource){
				case 0:
					calibData->SetPedestal(Values);
					break;
				case 1:
	    			        calibData->SetGain(Values);
					break;
				case 2:
	    			        calibData->SetSigma(Values);
					break;
			        }
		}
		delete sourceList;

	}
      
  
  // Check that everything was properly transmitted

//   for(Int_t j=0; j<128; j++){printf("ADCPedestal[%d] -> %f \n",j,calibData->GetPedestal(j));}
//   for(Int_t j=0; j<128; j++){printf("ADCGain[%d] -> %f \n",j,calibData->GetGain(j));}
//   for(Int_t j=0; j<128; j++){printf("ADCSigma[%d] -> %f \n",j,calibData->GetSigma(j));}
//   for(Int_t j=0; j<64; j++){printf("MeanHV[%d] -> %f \n",j,calibData->GetMeanHV(j));}
//   for(Int_t j=0; j<64; j++){printf("WidthHV[%d] -> %f \n",j,calibData->GetWidthHV(j));}
//   
  // Now we  store the VZERO Calibration Object into CalibrationDB
  
  AliCDBMetaData metaData;
  metaData.SetBeamPeriod(0);
  metaData.SetResponsible("Brigitte Cheynis");
  metaData.SetComment("This preprocessor fills an AliVZEROCalibData object");

  Bool_t result = Store("Calib", "Data", calibData, &metaData, 0, kTRUE);

  delete calibData;

  if (!result) return 1;   
  return 0;
}