1 #include "AliVZEROPreprocessor.h"
2 #include "AliVZEROCalibData.h"
3 #include "AliVZEROTriggerData.h"
4 #include "AliCDBMetaData.h"
5 #include "AliCDBEntry.h"
6 #include "AliDCSValue.h"
8 #include "AliShuttleInterface.h"
9 #include "AliVZERODataFEE.h"
10 #include "AliVZERODataDCS.h"
13 #include <TTimeStamp.h>
14 #include <TObjString.h>
22 // This class is a simple preprocessor for VZERO detector.
24 // It gets High Voltage values for a given run from DCS and Pedestal values from DAQ
25 // and writes them as Calibration MetaData into OCDB/VZERO/Calib/Data
26 // It also retrieves FEE parameters from DCS archive
27 // and writes them as Trigger MetaData into OCDB/VZERO/Trigger/Data
28 // (to be used for trigger simulation)
31 ClassImp(AliVZEROPreprocessor)
33 //______________________________________________________________________________________________
34 AliVZEROPreprocessor::AliVZEROPreprocessor(AliShuttleInterface* shuttle) :
35 AliPreprocessor("V00", shuttle),
42 AddRunType("STANDALONE_PULSER");
43 AddRunType("STANDALONE_BC");
44 AddRunType("PHYSICS");
48 //______________________________________________________________________________________________
49 AliVZEROPreprocessor::~AliVZEROPreprocessor()
57 //______________________________________________________________________________________________
58 void AliVZEROPreprocessor::Initialize(Int_t run, UInt_t startTime,
61 // Creates AliVZERODataDCS object
63 AliPreprocessor::Initialize(run, startTime, endTime);
65 Log(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
66 TTimeStamp(startTime).AsString(),
67 TTimeStamp(endTime).AsString()));
70 // fStartTime = startTime;
71 // fEndTime = endTime;
72 fStartTime = GetStartTimeDCSQuery ();
73 fEndTime = GetEndTimeDCSQuery ();
74 time_t daqStart = (time_t) (((TString)GetRunParameter("DAQ_time_start")).Atoi());
75 time_t daqEnd = (time_t) (((TString)GetRunParameter("DAQ_time_end")).Atoi());
76 time_t ctpStart = (time_t) (((TString)GetRunParameter("TRGTimeStart")).Atoi());
77 time_t ctpEnd = (time_t) (((TString)GetRunParameter("TRGTimeEnd")).Atoi());
79 fData = new AliVZERODataDCS(fRun, fStartTime, fEndTime,(UInt_t)daqStart, (UInt_t)daqEnd,(UInt_t)ctpStart, (UInt_t)ctpEnd);
80 fFEEData = new AliVZERODataFEE(fRun, fStartTime, fEndTime);
84 //______________________________________________________________________________________________
85 UInt_t AliVZEROPreprocessor::Process(TMap* dcsAliasMap)
87 // Fills data retrieved from DCS and DAQ into a AliVZEROCalibData object and
88 // stores it into CalibrationDB
91 // *** GET RUN TYPE ***
92 TString runType = GetRunType();
95 // *** REFERENCE DATA ***
98 AliVZEROCalibData *calibData = new AliVZEROCalibData();
100 // *************** HV From DCS ******************
101 // Fills data into a AliVZERODataDCS object
102 if(!dcsAliasMap) return 1;
104 // The Processing of the DCS input data is forwarded to AliVZERODataDCS
105 if (!fData->ProcessData(*dcsAliasMap)) return 1;
107 // Writes VZERO PMs HV values into VZERO calibration object and Timing resolution parameters
108 calibData->FillDCSData(fData);
110 // *************** From DAQ ******************
112 TString sourcesId = "V00da_results";
114 TList* sourceList = GetFileSources(kDAQ, sourcesId.Data());
116 Log(Form("No sources found for id %s", sourcesId.Data()));
118 Log(Form("The following sources produced files with the id %s",sourcesId.Data()));
121 TIter iter(sourceList);
124 while((source=dynamic_cast<TObjString*> (iter.Next()))){
125 fileName = GetFile(kDAQ, sourcesId.Data(), source->GetName());
126 if (fileName.Length() > 0)
127 Log(Form("Got the file %s, now we can extract some values.", fileName.Data()));
129 if((file = fopen(fileName.Data(),"r")) == NULL){
130 Log(Form("Cannot open file %s",fileName.Data()));
132 Float_t pedMean[128], pedSigma[128], adcMean[128], adcSigma[128] ;
133 for(Int_t j=0; j<128; j++) {
134 Int_t resScan = fscanf(file,"%f %f %f %f",
135 &pedMean[j], &pedSigma[j], &adcMean[j], &adcSigma[j]);
136 if (resScan != 4) Log(Form("Bad data in file %s !",fileName.Data()));
140 calibData->SetPedestal(pedMean);
141 calibData->SetSigma(pedSigma);
142 calibData->SetADCsigma(adcSigma);
147 // Check that everything was properly transmitted
149 // for(Int_t j=0; j<128; j++){printf("Pedestal[%d] -> %f \n",j,calibData->GetPedestal(j));}
150 // for(Int_t j=0; j<128; j++){printf("pedSigma[%d] -> %f \n",j,calibData->GetSigma(j));}
151 // for(Int_t j=0; j<128; j++){printf("Gain[%d] -> %f \n",j,calibData->GetGain(j));}
152 // for(Int_t j=0; j<128; j++){printf("adcSigma[%d] -> %f \n",j,calibData->GetADCsigma(j));}
153 // for(Int_t j=0; j<64; j++){printf("MeanHV[%d] -> %f \n",j,calibData->GetMeanHV(j));}
154 // for(Int_t j=0; j<64; j++){printf("WidthHV[%d] -> %f \n",j,calibData->GetWidthHV(j));}
156 // Now we store the VZERO Calibration Object into CalibrationDB
158 Bool_t resECal=kTRUE;
161 // if(sourceList && sourceList->GetEntries()>0)
163 AliCDBMetaData metaData;
164 metaData.SetBeamPeriod(0);
165 metaData.SetResponsible("Brigitte Cheynis");
166 metaData.SetComment("This preprocessor fills an AliVZEROCalibData object");
168 resECal = Store("Calib", "Data", calibData, &metaData, 0, kTRUE);
170 if(resECal==kFALSE ) result = 1;
176 // -----------------------------------------------------------------------
177 // Retrieves Front End Electronics Parameters from DCS archive
178 // -----------------------------------------------------------------------
179 AliVZEROTriggerData *triggerData = new AliVZEROTriggerData();
181 // The processing of the DCS input data is forwarded to AliVZERODataFEE
182 fFEEData->ProcessData(*dcsAliasMap);
184 // Writes VZERO FEE parameters values into VZERO Trigger parametrization object
185 triggerData->FillData(fFEEData);
187 // Stores the VZERO Trigger Object into TriggerDB
191 metaData.SetBeamPeriod(0);
192 metaData.SetResponsible("Brigitte Cheynis");
193 metaData.SetComment("This preprocessor fills an AliVZEROTriggerData object");
195 resECal = Store("Trigger", "Data", triggerData, &metaData, 0, kTRUE);
196 if(resECal==kFALSE ) result = 1;
199 // *************** From DAQ DA - Equalization factors ******************
201 TH1F *eqFactors = new TH1F("VZEROEqualizationFactors","VZERO Equalization Factors for Pb-Pb",64,-0.5,63.5);
202 sourcesId = "V00DAEqualFactors";
204 TList* sourceList2 = GetFileSources(kDAQ, sourcesId.Data());
206 Log(Form("No sources found for id %s", sourcesId.Data()));
208 Log(Form("The following sources produced files with the id %s",sourcesId.Data()));
209 sourceList2->Print();
211 TIter iter2(sourceList2);
214 while((source2=dynamic_cast<TObjString*> (iter2.Next()))){
215 fileName = GetFile(kDAQ, sourcesId.Data(), source2->GetName());
216 if (fileName.Length() > 0)
217 Log(Form("Got the file %s, now we can extract some values.", fileName.Data()));
219 if((file2 = fopen(fileName.Data(),"r")) == NULL){
220 Log(Form("Cannot open file %s",fileName.Data()));
224 alpha[0] = alpha[65] = 0;
227 for(Int_t j=0; j<64; ++j) {
228 Int_t resScan = fscanf(file2,"%d %f", &tempCh, &tempAlpha);
229 if (resScan != 2) Log(Form("Bad data in file %s !",fileName.Data()));
230 alpha[tempCh+1] = tempAlpha;
234 // Check that everything was properly transmitted
235 printf("Equalization factors (0->64): ");
236 for(Int_t j=0; j<64; ++j) printf("%.5f ",alpha[j+1]);
239 eqFactors->SetContent(alpha);
244 // Now we store the VZERO Equalization Factors Object into OCDB
248 AliCDBMetaData metaData2;
249 metaData2.SetBeamPeriod(0);
250 metaData2.SetResponsible("Brigitte Cheynis");
251 metaData2.SetComment("VZERO Equalization Factors object filled by VZERO preprocessor");
253 resECal = Store("Calib", "EqualizationFactors", eqFactors, &metaData2, 0, kTRUE);
255 if(resECal==kFALSE ) result = 1;