]>
Commit | Line | Data |
---|---|---|
bc943889 | 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 | /* | |
50e28e1c | 17 | $Log: AliT0Preprocessor.cxx,v $ |
18 | Revision 1.8 2007/12/07 15:22:51 alla | |
19 | bug fixed by Alberto | |
271b5bd3 | 20 | |
67b67126 | 21 | Revision 1.7 2007/12/06 16:35:24 alla |
22 | new bugs fixed by Tomek | |
23 | ||
02cf5377 | 24 | Revision 1.5 2007/11/23 19:28:52 alla |
25 | bug fixed | |
26 | ||
bc943889 | 27 | Version 2.1 2007/11/21 |
28 | Preprocessor storing data to OCDB (T.Malkiewicz) | |
29 | ||
30 | Version 1.1 2006/10 | |
31 | Preliminary test version (T.Malkiewicz) | |
32 | */ | |
387638f7 | 33 | // T0 preprocessor: |
34 | // 1) takes data from DCS and passes it to the class AliTOFDataDCS | |
35 | // for processing and writes the result to the Reference DB. | |
36 | // 2) takes data form DAQ (both from Laser Calibration and Physics runs), | |
37 | // processes it, and stores either to OCDB or to Reference DB. | |
38 | ||
bc943889 | 39 | |
dc7ca31d | 40 | #include "AliT0Preprocessor.h" |
bc943889 | 41 | #include "AliT0DataDCS.h" |
42 | #include "AliT0CalibWalk.h" | |
43 | #include "AliT0CalibTimeEq.h" | |
dc7ca31d | 44 | |
45 | #include "AliCDBMetaData.h" | |
46 | #include "AliDCSValue.h" | |
8f1790c4 | 47 | #include "AliCDBEntry.h" |
dc7ca31d | 48 | #include "AliLog.h" |
dc7ca31d | 49 | |
50 | #include <TTimeStamp.h> | |
51 | #include <TFile.h> | |
5221c818 | 52 | #include <TObjString.h> |
dc7ca31d | 53 | #include <TNamed.h> |
54 | #include "AliT0Dqclass.h" | |
6a9d154f | 55 | #include "TClass.h" |
dc7ca31d | 56 | |
bc943889 | 57 | |
dc7ca31d | 58 | ClassImp(AliT0Preprocessor) |
59 | ||
60 | //____________________________________________________ | |
387638f7 | 61 | AliT0Preprocessor::AliT0Preprocessor(AliShuttleInterface* shuttle) : |
62 | AliPreprocessor("T00", shuttle), | |
63 | fData(0) | |
dc7ca31d | 64 | { |
bc943889 | 65 | //constructor |
2a8e3409 | 66 | AddRunType("PHYSICS"); |
67 | AddRunType("STANDALONE"); | |
e98e34e9 | 68 | // AddRunType("LASER"); |
dc7ca31d | 69 | } |
bc943889 | 70 | //____________________________________________________ |
dc7ca31d | 71 | |
72 | AliT0Preprocessor::~AliT0Preprocessor() | |
73 | { | |
50e28e1c | 74 | //destructor |
bc943889 | 75 | delete fData; |
76 | fData = 0; | |
77 | } | |
78 | //____________________________________________________ | |
dc7ca31d | 79 | |
bc943889 | 80 | void AliT0Preprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) |
81 | { | |
50e28e1c | 82 | // Creates AliT0DataDCS object |
bc943889 | 83 | AliPreprocessor::Initialize(run, startTime, endTime); |
84 | AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp(startTime).AsString(), TTimeStamp(endTime).AsString())); | |
8eb4ac75 | 85 | fData = new AliT0DataDCS(fRun, fStartTime, fEndTime, GetStartTimeDCSQuery(), GetEndTimeDCSQuery()); |
dc7ca31d | 86 | } |
bc943889 | 87 | //____________________________________________________ |
dc7ca31d | 88 | |
49ab1618 | 89 | Bool_t AliT0Preprocessor::ProcessDCS(){ |
90 | // Check whether DCS should be processed or not... | |
91 | TString runType = GetRunType(); | |
92 | Log(Form("ProcessDCS - RunType: %s",runType.Data())); | |
93 | ||
4b5d6a0a | 94 | if((runType == "STANDALONE")|| |
02888358 | 95 | (runType == "PHYSICS") ) { |
96 | // || (runType == "LASER")){ | |
97 | ||
98 | // return kFALSE; | |
99 | return kTRUE; | |
49ab1618 | 100 | }else{ |
101 | return kFALSE; | |
102 | } | |
103 | } | |
104 | //____________________________________________________ | |
bc943889 | 105 | |
49ab1618 | 106 | UInt_t AliT0Preprocessor::ProcessDCSDataPoints(TMap* dcsAliasMap){ |
107 | // Fills data into AliT0DataDCS object | |
108 | Log("Processing DCS DP"); | |
bc943889 | 109 | Bool_t resultDCSMap=kFALSE; |
110 | Bool_t resultDCSStore=kFALSE; | |
49ab1618 | 111 | |
8bfd9a3e | 112 | if(!dcsAliasMap) |
113 | { | |
114 | Log("No DCS input data"); | |
49ab1618 | 115 | return 1; |
8bfd9a3e | 116 | } |
117 | else | |
118 | { | |
119 | resultDCSMap=fData->ProcessData(*dcsAliasMap); | |
120 | if(!resultDCSMap) | |
121 | { | |
122 | Log("Error when processing DCS data"); | |
bc943889 | 123 | return 2;// return error Code for processed DCS data not stored |
124 | } | |
125 | else | |
126 | { | |
127 | AliCDBMetaData metaDataDCS; | |
128 | metaDataDCS.SetBeamPeriod(0); | |
129 | metaDataDCS.SetResponsible("Tomasz Malkiewicz"); | |
130 | metaDataDCS.SetComment("This preprocessor fills an AliTODataDCS object."); | |
131 | AliInfo("Storing DCS Data"); | |
fe3c9b71 | 132 | resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS); |
bc943889 | 133 | if (!resultDCSStore) |
8bfd9a3e | 134 | { |
bc943889 | 135 | Log("Some problems occurred while storing DCS data results in ReferenceDB"); |
387638f7 | 136 | return 2;// return error Code for processed DCS data not stored |
bc943889 | 137 | } |
8bfd9a3e | 138 | } |
139 | } | |
49ab1618 | 140 | return 0; |
141 | } | |
142 | //____________________________________________________ | |
02888358 | 143 | |
49ab1618 | 144 | UInt_t AliT0Preprocessor::ProcessLaser(){ |
145 | // Processing data from DAQ Standalone run | |
8f1790c4 | 146 | Log("Processing Laser calibration - Walk Correction"); |
147 | ||
148 | //Retrieve the last T0 calibration object | |
149 | ||
150 | Float_t parqtcold[24][2], parledold[24][2],parqtcnew[24][2], parlednew[24][2] , goodled[24][2], goodqtc[24][2]; | |
151 | ||
6a9d154f | 152 | // std::cout<<"sizeof "<<sizeof(parqtcold)<<std::endl; |
153 | memset(parqtcold, 0, sizeof(parqtcold)); | |
154 | memset(parqtcnew, 0, sizeof(parqtcnew)); | |
155 | memset(parledold, 0, sizeof(parledold)); | |
156 | memset(parlednew, 0, sizeof(parlednew)); | |
157 | Int_t iStore=0; | |
02888358 | 158 | |
8f1790c4 | 159 | |
160 | AliT0CalibWalk* clb=0; | |
161 | AliCDBEntry* entryCalib = GetFromOCDB("Calib", "Slewing_Walk"); | |
162 | if(!entryCalib) | |
163 | Log(Form("Cannot find any AliCDBEntry for [Calib, SlewingWalk]!")); | |
164 | else { | |
6a9d154f | 165 | clb =dynamic_cast<AliT0CalibWalk*>(entryCalib->GetObject()); |
0c2f5b5a | 166 | for(Int_t i=0; i<24; i++) |
8f1790c4 | 167 | { |
0c2f5b5a | 168 | for(Int_t ipar=0; ipar<2; ipar++) |
8f1790c4 | 169 | { |
0c2f5b5a | 170 | // std::cout<<"parqtcold "<<parqtcold[i][ipar]<<std::endl; |
171 | parqtcold[i][ipar] = clb->GetQTCpar(i,ipar); | |
172 | parledold[i][ipar] = clb->GetLEDpar(i, ipar); | |
173 | goodqtc[i][ipar] = 999; | |
174 | goodled[i][ipar] = 999; | |
02888358 | 175 | // std:: cout<<" old "<<i<<" "<<ipar<<" qtc "<< parqtcold[i][ipar]<<" led "<<parledold[i][ipar]<< std::endl; |
8f1790c4 | 176 | } |
177 | } | |
0c2f5b5a | 178 | } |
179 | ||
6a9d154f | 180 | |
181 | Bool_t resultLaser=kFALSE; | |
182 | //processing DAQ | |
183 | TList* list = GetFileSources(kDAQ, "LASER"); | |
184 | if (list) | |
185 | { | |
186 | TIter iter(list); | |
187 | TObjString *source; | |
188 | while ((source = dynamic_cast<TObjString *> (iter.Next()))) | |
189 | { | |
190 | const char *laserFile = GetFile(kDAQ, "LASER", source->GetName()); | |
191 | if (laserFile) | |
bc943889 | 192 | { |
02cf5377 | 193 | Log(Form("File with Id LASER found in source %s!", source->GetName())); |
8bfd9a3e | 194 | AliT0CalibWalk *laser = new AliT0CalibWalk(); |
8bfd9a3e | 195 | laser->MakeWalkCorrGraph(laserFile); |
8f1790c4 | 196 | //check difference with what was before |
0c2f5b5a | 197 | if(laser && clb ){ |
6a9d154f | 198 | iStore = 1; |
8f1790c4 | 199 | for(Int_t i=0; i<24; i++) |
200 | { | |
201 | for(Int_t ifit=0; ifit<2; ifit++) | |
202 | { | |
203 | parqtcnew[i][ifit] = laser->GetQTCpar(i,ifit); | |
204 | if( parqtcold[i][ifit] != 0 && parqtcnew[i][ifit] !=0) | |
205 | { | |
206 | goodqtc[i][ifit] = | |
207 | (parqtcnew[i][ifit] - parqtcold[i][ifit])/parqtcold[i][ifit]; | |
02888358 | 208 | // std::cout<<"qtc "<<i<<" "<<ifit<<" "<< goodqtc[i][ifit]<< std::endl; |
8f1790c4 | 209 | } |
210 | parlednew[i][ifit] = laser->GetLEDpar(i,ifit); | |
6a9d154f | 211 | if(parledold[i][ifit] != 0 && parlednew[i][ifit]!= 0 ) |
8f1790c4 | 212 | { |
213 | goodled[i][ifit]= | |
214 | (parlednew[i][ifit] - parledold[i][ifit])/parledold[i][ifit]; | |
02888358 | 215 | // std::cout<<"led "<<i<<" "<<ifit<<" "<< goodled[i][ifit]<< std::endl; |
8f1790c4 | 216 | } |
02888358 | 217 | // if(TMath::Abs(goodqtc[i][ifit])>0.1 || |
218 | // TMath::Abs(goodled[i][ifit])>0.1) | |
219 | // iStore = 0; | |
8f1790c4 | 220 | } |
221 | } | |
222 | } | |
6a9d154f | 223 | |
8f1790c4 | 224 | AliCDBMetaData metaData; |
225 | metaData.SetBeamPeriod(0); | |
bc943889 | 226 | metaData.SetResponsible("Tomek&Michal"); |
227 | metaData.SetComment("Walk correction from laser runs."); | |
8f1790c4 | 228 | if( iStore>0) |
229 | resultLaser=Store("Calib","Slewing_Walk", laser, &metaData, 0, 1); | |
bc943889 | 230 | delete laser; |
271b5bd3 | 231 | Log(Form("resultLaser = %d",resultLaser)); |
8bfd9a3e | 232 | } |
bc943889 | 233 | else |
234 | { | |
8bfd9a3e | 235 | Log(Form("Could not find file with Id LASER in source %s!", source->GetName())); |
bc943889 | 236 | return 1; |
237 | } | |
238 | } | |
239 | if (!resultLaser) | |
8bfd9a3e | 240 | { |
bc943889 | 241 | Log("No Laser Data stored"); |
242 | return 3;//return error code for failure in storing Laser Data | |
243 | } | |
67b67126 | 244 | } else { |
245 | Log("No sources found for id LASER!"); | |
246 | return 1; | |
247 | } | |
49ab1618 | 248 | return 0; |
249 | } | |
02888358 | 250 | |
49ab1618 | 251 | //____________________________________________________ |
252 | ||
253 | UInt_t AliT0Preprocessor::ProcessPhysics(){ | |
254 | //Processing data from DAQ Physics run | |
255 | Log("Processing Physics"); | |
02888358 | 256 | |
271b5bd3 | 257 | Bool_t resultOnline=kFALSE; |
49ab1618 | 258 | //processing DAQ |
259 | TList* listPhys = GetFileSources(kDAQ, "PHYSICS"); | |
260 | if (listPhys) | |
8bfd9a3e | 261 | { |
262 | TIter iter(listPhys); | |
263 | TObjString *sourcePhys; | |
264 | while ((sourcePhys = dynamic_cast<TObjString *> (iter.Next()))) | |
265 | { | |
266 | const char *filePhys = GetFile(kDAQ, "PHYSICS", sourcePhys->GetName()); | |
267 | if (filePhys) | |
268 | { | |
269 | AliT0CalibTimeEq *online = new AliT0CalibTimeEq(); | |
270 | online->Reset(); | |
02888358 | 271 | Bool_t writeok = online->ComputeOnlineParams(filePhys); |
8bfd9a3e | 272 | AliCDBMetaData metaData; |
273 | metaData.SetBeamPeriod(0); | |
02888358 | 274 | metaData.SetResponsible("Alla Maevskaya"); |
8bfd9a3e | 275 | metaData.SetComment("Time equalizing result."); |
02888358 | 276 | |
277 | if (writeok) resultOnline = Store("Calib","TimeDelay", online, &metaData, 0, 1); | |
278 | else { | |
279 | ||
280 | Log(Form("writeok = %d not enough data for equalizing",resultOnline)); | |
281 | return 0; | |
282 | } | |
271b5bd3 | 283 | Log(Form("resultOnline = %d",resultOnline)); |
8bfd9a3e | 284 | delete online; |
285 | } | |
271b5bd3 | 286 | else |
8bfd9a3e | 287 | { |
288 | Log(Form("Could not find file with Id PHYSICS in source %s!", sourcePhys->GetName())); | |
289 | return 1; | |
290 | } | |
271b5bd3 | 291 | |
8bfd9a3e | 292 | } |
293 | if (!resultOnline) | |
bc943889 | 294 | { |
271b5bd3 | 295 | Log("No Data stored"); |
bc943889 | 296 | return 4;//return error code for failure in storing OCDB Data |
297 | } | |
67b67126 | 298 | } else { |
299 | Log("No sources found for id PHYSICS!"); | |
300 | return 1; | |
301 | } | |
49ab1618 | 302 | return 0; |
303 | } | |
304 | //____________________________________________________ | |
305 | ||
306 | UInt_t AliT0Preprocessor::Process(TMap* dcsAliasMap ) | |
307 | { | |
308 | // T0 preprocessor return codes: | |
309 | // return=0 : all ok | |
310 | // return=1 : no DCS input data | |
311 | // return=2 : failed to store DCS data | |
312 | // return=3 : no Laser data (Walk correction) | |
313 | // return=4 : failed to store OCDB time equalized data | |
314 | // return=5 : no DAQ input for OCDB | |
315 | // return=6 : failed to retrieve DAQ data from OCDB | |
316 | // return=7 : failed to store T0 OCDB data | |
02888358 | 317 | // return=8 : not enough data for equalizing |
318 | Bool_t dcsDP = ProcessDCS(); | |
319 | Log(Form("dcsDP = %d",dcsDP)); | |
320 | TString runType = GetRunType(); | |
321 | Log(Form("RunType: %s",runType.Data())); | |
322 | //processing | |
323 | if(runType == "STANDALONE"){ | |
324 | if(dcsDP==1){ | |
325 | Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap); | |
326 | return iresultDCS; | |
327 | } | |
328 | } | |
329 | /* | |
330 | if(runType == "LASER"){ | |
331 | Int_t iresultLaser = ProcessLaser(); | |
332 | if(dcsDP==1){ | |
333 | Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap); | |
334 | return iresultDCS; | |
335 | } | |
336 | ||
337 | Log(Form("iresultLaser = %d",iresultLaser)); | |
338 | return iresultLaser; | |
339 | } | |
340 | */ | |
341 | else if(runType == "PHYSICS"){ | |
342 | Int_t iresultPhysics = ProcessPhysics(); | |
343 | if(dcsDP==1){ | |
344 | Int_t iresultDCS = ProcessDCSDataPoints(dcsAliasMap); | |
345 | return iresultDCS; | |
346 | } | |
347 | Log(Form("iresultPhysics = %d",iresultPhysics)); | |
79e2e9c6 | 348 | return iresultPhysics; |
02888358 | 349 | } |
350 | ||
351 | ||
79e2e9c6 | 352 | |
353 | return 0; | |
dc7ca31d | 354 | } |