1 #include "AliVZEROPreprocessor.h"
2 #include "AliVZEROCalibData.h"
3 #include "AliCDBMetaData.h"
4 #include "AliCDBEntry.h"
5 #include "AliDCSValue.h"
8 #include <TTimeStamp.h>
9 #include <TObjString.h>
14 // This class is a simple preprocessor for V0.
17 ClassImp(AliVZEROPreprocessor)
19 //______________________________________________________________________________________________
20 AliVZEROPreprocessor::AliVZEROPreprocessor(AliShuttleInterface* shuttle) :
21 AliPreprocessor("V00", shuttle),
27 AddRunType("PHYSICS");
31 //______________________________________________________________________________________________
32 AliVZEROPreprocessor::~AliVZEROPreprocessor()
39 //______________________________________________________________________________________________
40 void AliVZEROPreprocessor::Initialize(Int_t run, UInt_t startTime,
43 // Creates AliZDCDataDCS object
45 AliPreprocessor::Initialize(run, startTime, endTime);
47 Log(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
48 TTimeStamp(startTime).AsString(),
49 TTimeStamp(endTime).AsString()));
52 fStartTime = startTime;
55 fData = new AliVZERODataDCS(fRun, fStartTime, fEndTime);
59 //______________________________________________________________________________________________
60 UInt_t AliVZEROPreprocessor::Process(TMap* dcsAliasMap)
62 // Fills data retrieved from DCS and DAQ into a AliVZEROCalibData object and
63 // stores it into CalibrationDB
66 // *** GET RUN TYPE ***
67 TString runType = GetRunType();
70 // *** REFERENCE DATA ***
73 AliVZEROCalibData *calibData = new AliVZEROCalibData();
75 // *************** From DCS ******************
76 // Fills data into a AliVZERODataDCS object
77 if(!dcsAliasMap) return 1;
79 // The processing of the DCS input data is forwarded to AliVZERODataDCS
81 fData->ProcessData(*dcsAliasMap);
82 //fData->Draw(""); // Draws the HV values as a function of time
83 //dcsAliasMap->Print(""); // Prints out the HV values
85 // Writes VZERO PMs HV values into VZERO calibration object
86 calibData->SetMeanHV(fData->GetMeanHV());
87 calibData->SetWidthHV(fData->GetWidthHV());
89 // *************** From DAQ ******************
91 TString SourcesId = "V00da_results";
93 TList* sourceList = GetFileSources(kDAQ, SourcesId.Data());
95 Log(Form("No sources found for id %s", SourcesId.Data()));
97 Log(Form("The following sources produced files with the id %s",SourcesId.Data()));
100 TIter iter(sourceList);
101 TObjString *source = 0;
103 while((source=dynamic_cast<TObjString*> (iter.Next()))){
104 fileName = GetFile(kDAQ, SourcesId.Data(), source->GetName());
105 if (fileName.Length() > 0)
106 Log(Form("Got the file %s, now we can extract some values.", fileName.Data()));
108 if((file = fopen(fileName.Data(),"r")) == NULL){
109 Log(Form("Cannot open file %s",fileName.Data()));
111 Float_t PEDmean[128], PEDsigma[128], ADCmean[128], ADCsigma[128] ;
112 for(Int_t j=0; j<128; j++)fscanf(file,"%f %f %f %f",
113 &PEDmean[j], &PEDsigma[j], &ADCmean[j], &ADCsigma[j]);
116 calibData->SetPedestal(PEDmean);
117 calibData->SetSigma(PEDsigma);
118 calibData->SetGain(ADCmean);
119 calibData->SetADCsigma(ADCsigma);
124 // Check that everything was properly transmitted
126 // for(Int_t j=0; j<128; j++){printf("Pedestal[%d] -> %f \n",j,calibData->GetPedestal(j));}
127 // for(Int_t j=0; j<128; j++){printf("PedSigma[%d] -> %f \n",j,calibData->GetSigma(j));}
128 // for(Int_t j=0; j<128; j++){printf("Gain[%d] -> %f \n",j,calibData->GetGain(j));}
129 // for(Int_t j=0; j<128; j++){printf("ADCsigma[%d] -> %f \n",j,calibData->GetADCsigma(j));}
130 // for(Int_t j=0; j<64; j++){printf("MeanHV[%d] -> %f \n",j,calibData->GetMeanHV(j));}
131 // for(Int_t j=0; j<64; j++){printf("WidthHV[%d] -> %f \n",j,calibData->GetWidthHV(j));}
133 // Now we store the VZERO Calibration Object into CalibrationDB
135 Bool_t result = kFALSE;
136 if(sourceList && sourceList->GetEntries()>0)
138 AliCDBMetaData metaData;
139 metaData.SetBeamPeriod(0);
140 metaData.SetResponsible("Brigitte Cheynis");
141 metaData.SetComment("This preprocessor fills an AliVZEROCalibData object");
143 result = Store("Calib", "Data", calibData, &metaData, 0, kTRUE);
149 if (!result) return 1;