Added marking SHUTTLE=DONE for invalid runs
[u/mrichter/AliRoot.git] / TPC / AliTPCPreprocessor.cxx
CommitLineData
54472e4f 1/**************************************************************************
2 * Copyright(c) 2007, 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
72df5829 16
54472e4f 17#include "AliTPCPreprocessor.h"
a7dce0bc 18#include "AliShuttleInterface.h"
54472e4f 19
20#include "AliCDBMetaData.h"
21#include "AliDCSValue.h"
22#include "AliLog.h"
23#include "AliTPCSensorTempArray.h"
a7dce0bc 24#include "AliTPCROC.h"
25#include "AliTPCCalROC.h"
26#include "AliTPCCalPad.h"
27#include "AliTPCCalibPedestal.h"
9f016d99 28#include "AliTPCCalibPulser.h"
29#include "AliTPCCalibCE.h"
9e66a8dd 30#include "TFile.h"
7b54c8e9 31#include "TTree.h"
32#include "TEnv.h"
54472e4f 33
34#include <TTimeStamp.h>
35
a7dce0bc 36const Int_t kValCutTemp = 100; // discard temperatures > 100 degrees
37const Int_t kDiffCutTemp = 5; // discard temperature differences > 5 degrees
38const TString kPedestalRunType = "PEDESTAL_RUN"; // pedestal run identifier
9f016d99 39const TString kPulserRunType = "PULSER_RUN"; // pulser run identifier
bdffc5fb 40const TString kPhysicsRunType = "PHYSICS"; // physics run identifier
67a165ed 41const TString kAmandaTemp = "tpc_PT_%d.Temperature"; // Amanda string for temperature entries
42const Double_t kFitFraction = 0.7; // Fraction of DCS sensor fits required
72df5829 43
54472e4f 44//
45// This class is the SHUTTLE preprocessor for the TPC detector.
f7f602cc 46// It contains several components, this far the part containing
54472e4f 47// temperatures is implemented
48//
49
50ClassImp(AliTPCPreprocessor)
51
52//______________________________________________________________________________________________
6d07bf74 53AliTPCPreprocessor::AliTPCPreprocessor(AliShuttleInterface* shuttle) :
54 AliPreprocessor("TPC",shuttle),
bdffc5fb 55 fConfEnv(0), fTemp(0), fHighVoltage(0), fConfigOK(kTRUE), fROC(0)
54472e4f 56{
57 // constructor
a7dce0bc 58 fROC = AliTPCROC::Instance();
54472e4f 59}
72df5829 60//______________________________________________________________________________________________
61// AliTPCPreprocessor::AliTPCPreprocessor(const AliTPCPreprocessor& org) :
62// AliPreprocessor(org),
bdffc5fb 63// fConfEnv(0), fTemp(0), fHighVoltage(0), fConfigOK(kTRUE)
72df5829 64// {
65// // copy constructor not implemented
f7f602cc 66// // -- missing underlying copy constructor in AliPreprocessor
67//
72df5829 68// Fatal("AliTPCPreprocessor", "copy constructor not implemented");
f7f602cc 69//
70// // fTemp = new AliTPCSensorTempArray(*(org.fTemp));
72df5829 71// }
54472e4f 72
73//______________________________________________________________________________________________
74AliTPCPreprocessor::~AliTPCPreprocessor()
75{
76 // destructor
a7dce0bc 77
54472e4f 78 delete fTemp;
bdffc5fb 79 delete fHighVoltage;
54472e4f 80}
72df5829 81//______________________________________________________________________________________________
82AliTPCPreprocessor& AliTPCPreprocessor::operator = (const AliTPCPreprocessor& )
83{
84 Fatal("operator =", "assignment operator not implemented");
85 return *this;
86}
87
54472e4f 88
89//______________________________________________________________________________________________
90void AliTPCPreprocessor::Initialize(Int_t run, UInt_t startTime,
91 UInt_t endTime)
92{
93 // Creates AliTestDataDCS object
94
95 AliPreprocessor::Initialize(run, startTime, endTime);
96
97 AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
98 TTimeStamp(startTime).AsString(),
99 TTimeStamp(endTime).AsString()));
100
7b54c8e9 101 // Preprocessor configuration
102
37092b2b 103 AliCDBEntry* entry = GetFromOCDB("Config", "Preprocessor");
7b54c8e9 104 if (entry) fConfEnv = (TEnv*) entry->GetObject();
105 if ( fConfEnv==0 ) {
7b54c8e9 106 Log("AliTPCPreprocsessor: Preprocessor Config OCDB entry missing.\n");
37092b2b 107 fConfigOK = kFALSE;
108 return;
7b54c8e9 109 }
110
18eade96 111 // Temperature sensors
112
f7f602cc 113 TTree *confTree = 0;
7b54c8e9 114 entry = GetFromOCDB("Config", "Temperature");
f7f602cc 115 if (entry) confTree = (TTree*) entry->GetObject();
18eade96 116 if ( confTree==0 ) {
18eade96 117 Log("AliTPCPreprocsessor: Temperature Config OCDB entry missing.\n");
118 fConfigOK = kFALSE;
f7f602cc 119 return;
18eade96 120 }
67a165ed 121 fTemp = new AliTPCSensorTempArray(fStartTime, fEndTime, confTree, kAmandaTemp);
18eade96 122 fTemp->SetValCut(kValCutTemp);
123 fTemp->SetDiffCut(kDiffCutTemp);
f7f602cc 124
18eade96 125
bdffc5fb 126 // High voltage measurements
127
128 if (false) { // high voltage maps not yet implemented...
129 confTree=0;
130 entry=0;
131 entry = GetFromOCDB("Config", "HighVoltage");
132 if (entry) confTree = (TTree*) entry->GetObject();
133 if ( confTree==0 ) {
134 Log("AliTPCPreprocsessor: High Voltage Config OCDB entry missing.\n");
135 fConfigOK = kFALSE;
136 return;
137 }
138 fHighVoltage = new AliDCSSensorArray(fStartTime, fEndTime, confTree);
139 }
54472e4f 140}
141
142//______________________________________________________________________________________________
143UInt_t AliTPCPreprocessor::Process(TMap* dcsAliasMap)
144{
145 // Fills data into TPC calibrations objects
146
54472e4f 147 // Amanda servers provide information directly through dcsAliasMap
148
a7dce0bc 149 if (!dcsAliasMap) return 9;
7b54c8e9 150 if (dcsAliasMap->GetEntries() == 0 ) return 9;
a7dce0bc 151 if (!fConfigOK) return 9;
152
153 TString runType = GetRunType();
154
54472e4f 155 // Temperature sensors are processed by AliTPCCalTemp
72df5829 156
18eade96 157
54472e4f 158 UInt_t tempResult = MapTemperature(dcsAliasMap);
159 UInt_t result=tempResult;
18eade96 160
bdffc5fb 161 // High Voltage recordings
162
163
164 if (false) { // High Voltage maps not yet implemented..
165
166 UInt_t hvResult = MapHighVoltage(dcsAliasMap);
167 result+=hvResult;
168
169 }
f7f602cc 170
54472e4f 171 // Other calibration information will be retrieved through FXS files
a7dce0bc 172 // examples:
54472e4f 173 // TList* fileSourcesDAQ = GetFile(AliShuttleInterface::kDAQ, "pedestals");
174 // const char* fileNamePed = GetFile(AliShuttleInterface::kDAQ, "pedestals", "LDC1");
175 //
176 // TList* fileSourcesHLT = GetFile(AliShuttleInterface::kHLT, "calib");
177 // const char* fileNameHLT = GetFile(AliShuttleInterface::kHLT, "calib", "LDC1");
178
9f016d99 179 // pedestal entries
54472e4f 180
a7dce0bc 181 if(runType == kPedestalRunType) {
7b54c8e9 182 Int_t pedestalSource = AliShuttleInterface::kDAQ;
183 TString source = fConfEnv->GetValue("Pedestal","DAQ");
184 source.ToUpper();
185 if ( source == "HLT" ) pedestalSource = AliShuttleInterface::kHLT;
9f016d99 186 if (!GetHLTStatus()) pedestalSource = AliShuttleInterface::kDAQ;
7b54c8e9 187 UInt_t pedestalResult = ExtractPedestals(pedestalSource);
a7dce0bc 188 result += pedestalResult;
7b54c8e9 189
a7dce0bc 190 }
191
9f016d99 192 // pulser trigger processing
193
194 if( runType == kPulserRunType ) {
195 Int_t pulserSource = AliShuttleInterface::kDAQ;
196 TString source = fConfEnv->GetValue("Pulser","DAQ");
197 source.ToUpper();
198 if ( source == "HLT" ) pulserSource = AliShuttleInterface::kHLT;
199 if (!GetHLTStatus()) pulserSource = AliShuttleInterface::kDAQ;
200 UInt_t pulserResult = ExtractPulser(pulserSource);
201 result += pulserResult;
202
203 }
204
205 // Central Electrode processing
206
bdffc5fb 207 if( runType == kPhysicsRunType ) { // CE input file not generated yet
9f016d99 208 Int_t ceSource = AliShuttleInterface::kDAQ;
209 TString source = fConfEnv->GetValue("CE","DAQ");
210 source.ToUpper();
211 if ( source == "HLT" ) ceSource = AliShuttleInterface::kHLT;
212 if (!GetHLTStatus()) ceSource = AliShuttleInterface::kDAQ;
213 UInt_t ceResult = ExtractCE(ceSource);
214 result += ceResult;
215
216 }
a7dce0bc 217
54472e4f 218 return result;
219}
220//______________________________________________________________________________________________
221UInt_t AliTPCPreprocessor::MapTemperature(TMap* dcsAliasMap)
222{
223
224 // extract DCS temperature maps. Perform fits to save space
225
72df5829 226 UInt_t result=0;
54472e4f 227 TMap *map = fTemp->ExtractDCS(dcsAliasMap);
228 if (map) {
229 fTemp->MakeSplineFit(map);
67a165ed 230 Double_t fitFraction = fTemp->NumFits()/fTemp->NumSensors();
231 if (fitFraction > kFitFraction ) {
232 AliInfo(Form("Temperature values extracted, fits performed.\n"));
233 } else {
234 Log ("Too few temperature maps fitted. \n");
235 result = 9;
236 }
54472e4f 237 } else {
67a165ed 238 Log("No temperature map extracted. \n");
72df5829 239 result=9;
54472e4f 240 }
241 delete map;
242 // Now store the final CDB file
a7dce0bc 243
244 if ( result == 0 ) {
72df5829 245 AliCDBMetaData metaData;
54472e4f 246 metaData.SetBeamPeriod(0);
247 metaData.SetResponsible("Haavard Helstrup");
248 metaData.SetComment("Preprocessor AliTPC data base entries.");
249
a7dce0bc 250 Bool_t storeOK = Store("Calib", "Temperature", fTemp, &metaData, 0, kFALSE);
251 if ( !storeOK ) result=1;
252
72df5829 253 }
54472e4f 254
255 return result;
7b54c8e9 256
54472e4f 257}
7b54c8e9 258
bdffc5fb 259//______________________________________________________________________________________________
260UInt_t AliTPCPreprocessor::MapHighVoltage(TMap* dcsAliasMap)
261{
262
263 // extract DCS HV maps. Perform fits to save space
264
265 UInt_t result=0;
266 TMap *map = fHighVoltage->ExtractDCS(dcsAliasMap);
267 if (map) {
268 fHighVoltage->MakeSplineFit(map);
269 Double_t fitFraction = fHighVoltage->NumFits()/fHighVoltage->NumSensors();
270 if (fitFraction > kFitFraction ) {
271 AliInfo(Form("High volatge recordings extracted, fits performed.\n"));
272 } else {
273 Log ("Too few high voltage recordings fitted. \n");
274 result = 9;
275 }
276 } else {
277 Log("No high voltage recordings extracted. \n");
278 result=9;
279 }
280 delete map;
281 // Now store the final CDB file
282
283 if ( result == 0 ) {
284 AliCDBMetaData metaData;
285 metaData.SetBeamPeriod(0);
286 metaData.SetResponsible("Haavard Helstrup");
287 metaData.SetComment("Preprocessor AliTPC data base entries.");
288
289 Bool_t storeOK = Store("Calib", "HighVoltage", fHighVoltage, &metaData, 0, kFALSE);
290 if ( !storeOK ) result=1;
291
292 }
293
294 return result;
295
296}
297
7b54c8e9 298
a7dce0bc 299//______________________________________________________________________________________________
7b54c8e9 300
301UInt_t AliTPCPreprocessor::ExtractPedestals(Int_t sourceFXS)
a7dce0bc 302{
303 //
304 // Read pedestal file from file exchage server
305 // Keep original entry from OCDB in case no new pedestals are available
306 //
307 AliTPCCalPad *calPadPed=0;
308 AliCDBEntry* entry = GetFromOCDB("Calib", "Pedestals");
309 if (entry) calPadPed = (AliTPCCalPad*)entry->GetObject();
310 if ( calPadPed==NULL ) {
a7dce0bc 311 Log("AliTPCPreprocsessor: No previous TPC pedestal entry available.\n");
7b54c8e9 312 calPadPed = new AliTPCCalPad("PedestalsMean","PedestalsMean");
313 }
314
315 AliTPCCalPad *calPadRMS=0;
316 entry = GetFromOCDB("Calib", "Noise");
317 if (entry) calPadRMS = (AliTPCCalPad*)entry->GetObject();
318 if ( calPadRMS==NULL ) {
7b54c8e9 319 Log("AliTPCPreprocsessor: No previous TPC noise entry available.\n");
320 calPadRMS = new AliTPCCalPad("PedestalsRMS","PedestalsRMS");
a7dce0bc 321 }
322
7b54c8e9 323
a7dce0bc 324 UInt_t result=0;
325
326 Int_t nSectors = fROC->GetNSectors();
7b54c8e9 327 TList* list = GetFileSources(sourceFXS,"pedestals");
9f016d99 328
329 if (list && list->GetEntries()>0) {
a7dce0bc 330
331// loop through all files from LDCs
332
bdffc5fb 333 Bool_t changed=false;
a7dce0bc 334 UInt_t index = 0;
335 while (list->At(index)!=NULL) {
336 TObjString* fileNameEntry = (TObjString*) list->At(index);
337 if (fileNameEntry!=NULL) {
7b54c8e9 338 TString fileName = GetFile(sourceFXS, "pedestals",
a7dce0bc 339 fileNameEntry->GetString().Data());
340 TFile *f = TFile::Open(fileName);
9f016d99 341 if (!f) {
342 Log ("Error opening pedestal file.");
343 result =2;
344 break;
345 }
a7dce0bc 346 AliTPCCalibPedestal *calPed;
347 f->GetObject("AliTPCCalibPedestal",calPed);
67a165ed 348 if ( !calPed ) {
349 Log ("No pedestal calibration object in file.");
350 result = 2;
351 break;
352 }
a7dce0bc 353
354 // replace entries for the sectors available in the present file
355
bdffc5fb 356 changed=true;
c9c23d80 357 for (Int_t sector=0; sector<nSectors; sector++) {
7b54c8e9 358 AliTPCCalROC *rocPed=calPed->GetCalRocPedestal(sector, kFALSE);
359 if ( rocPed ) calPadPed->SetCalROC(rocPed,sector);
360 AliTPCCalROC *rocRMS=calPed->GetCalRocRMS(sector, kFALSE);
361 if ( rocRMS ) calPadRMS->SetCalROC(rocRMS,sector);
a7dce0bc 362 }
363 }
364 ++index;
365 } // while(list)
366//
367// Store updated pedestal entry to OCDB
368//
bdffc5fb 369 if (changed) {
370 AliCDBMetaData metaData;
371 metaData.SetBeamPeriod(0);
372 metaData.SetResponsible("Haavard Helstrup");
373 metaData.SetComment("Preprocessor AliTPC data base entries.");
374
375 Bool_t storeOK = Store("Calib", "Pedestals", calPadPed, &metaData, 0, kTRUE);
376 if ( !storeOK ) ++result;
377 storeOK = Store("Calib", "PadNoise", calPadRMS, &metaData, 0, kTRUE);
378 if ( !storeOK ) ++result;
379 }
9f016d99 380 } else {
67a165ed 381 Log ("Error: no entries in input file list!");
9f016d99 382 result = 1;
7b54c8e9 383 }
a7dce0bc 384
a7dce0bc 385 return result;
386}
7b54c8e9 387
388//______________________________________________________________________________________________
389
390
9f016d99 391UInt_t AliTPCPreprocessor::ExtractPulser(Int_t sourceFXS)
392{
393 //
394 // Read pulser calibration file from file exchage server
395 // Keep original entry from OCDB in case no new pulser calibration is available
396 //
397 TObjArray *pulserObjects=0;
398 AliTPCCalPad *pulserTmean=0;
399 AliTPCCalPad *pulserTrms=0;
400 AliTPCCalPad *pulserQmean=0;
401 AliCDBEntry* entry = GetFromOCDB("Calib", "Pulser");
402 if (entry) pulserObjects = (TObjArray*)entry->GetObject();
403 if ( pulserObjects==NULL ) {
404 Log("AliTPCPreprocsessor: No previous TPC pulser entry available.\n");
405 pulserObjects = new TObjArray;
406 }
407
408 pulserTmean = (AliTPCCalPad*)pulserObjects->FindObject("PulserTmean");
409 if ( !pulserTmean ) {
410 pulserTmean = new AliTPCCalPad("PulserTmean","PulserTmean");
411 pulserObjects->Add(pulserTmean);
412 }
413 pulserTrms = (AliTPCCalPad*)pulserObjects->FindObject("PulserTrms");
414 if ( !pulserTrms ) {
415 pulserTrms = new AliTPCCalPad("PulserTrms","PulserTrms");
416 pulserObjects->Add(pulserTrms);
417 }
418 pulserQmean = (AliTPCCalPad*)pulserObjects->FindObject("PulserQmean");
419 if ( !pulserQmean ) {
420 pulserQmean = new AliTPCCalPad("PulserQmean","PulserQmean");
421 pulserObjects->Add(pulserQmean);
422 }
423
424
425 UInt_t result=0;
426
427 Int_t nSectors = fROC->GetNSectors();
428 TList* list = GetFileSources(sourceFXS,"pulser");
429
430 if (list && list->GetEntries()>0) {
431
432// loop through all files from LDCs
433
bdffc5fb 434 Bool_t changed=false;
9f016d99 435 UInt_t index = 0;
436 while (list->At(index)!=NULL) {
437 TObjString* fileNameEntry = (TObjString*) list->At(index);
438 if (fileNameEntry!=NULL) {
439 TString fileName = GetFile(sourceFXS, "pulser",
440 fileNameEntry->GetString().Data());
441 TFile *f = TFile::Open(fileName);
442 if (!f) {
443 Log ("Error opening pulser file.");
444 result =2;
445 break;
446 }
447 AliTPCCalibPulser *calPulser;
448 f->GetObject("AliTPCCalibPulser",calPulser);
67a165ed 449 if ( !calPulser ) {
450 Log ("No pulser calibration object in file.");
451 result = 2;
452 break;
453 }
9f016d99 454
455 // replace entries for the sectors available in the present file
456
bdffc5fb 457 changed=true;
9f016d99 458 for (Int_t sector=0; sector<nSectors; sector++) {
459 AliTPCCalROC *rocTmean=calPulser->GetCalRocT0(sector);
460 if ( rocTmean ) pulserTmean->SetCalROC(rocTmean,sector);
461 AliTPCCalROC *rocTrms=calPulser->GetCalRocRMS(sector);
462 if ( rocTrms ) pulserTrms->SetCalROC(rocTrms,sector);
463 AliTPCCalROC *rocQmean=calPulser->GetCalRocQ(sector);
464 if ( rocQmean ) pulserQmean->SetCalROC(rocQmean,sector);
465 }
466 }
467 ++index;
468 } // while(list)
469//
470// Store updated pedestal entry to OCDB
471//
bdffc5fb 472 if (changed) {
473 AliCDBMetaData metaData;
474 metaData.SetBeamPeriod(0);
475 metaData.SetResponsible("Haavard Helstrup");
476 metaData.SetComment("Preprocessor AliTPC data base entries.");
477
478 Bool_t storeOK = Store("Calib", "Pulser", pulserObjects, &metaData, 0, kTRUE);
479 if ( !storeOK ) ++result;
480 }
9f016d99 481 } else {
67a165ed 482 Log ("Error: no entries in input file list!");
9f016d99 483 result = 1;
484 }
485
486 return result;
487}
488
489UInt_t AliTPCPreprocessor::ExtractCE(Int_t sourceFXS)
490{
491 //
492 // Read Central Electrode file from file exchage server
493 // Keep original entry from OCDB in case no new CE calibration is available
494 //
495 TObjArray *ceObjects=0;
496 AliTPCCalPad *ceTmean=0;
497 AliTPCCalPad *ceTrms=0;
498 AliTPCCalPad *ceQmean=0;
499 AliCDBEntry* entry = GetFromOCDB("Calib", "CE");
500 if (entry) ceObjects = (TObjArray*)entry->GetObject();
501 if ( ceObjects==NULL ) {
502 Log("AliTPCPreprocsessor: No previous TPC central electrode entry available.\n");
503 ceObjects = new TObjArray;
504 }
505
506 ceTmean = (AliTPCCalPad*)ceObjects->FindObject("CETmean");
507 if ( !ceTmean ) {
508 ceTmean = new AliTPCCalPad("CETmean","CETmean");
509 ceObjects->Add(ceTmean);
510 }
511 ceTrms = (AliTPCCalPad*)ceObjects->FindObject("CETrms");
512 if ( !ceTrms ) {
513 ceTrms = new AliTPCCalPad("CETrms","CETrms");
514 ceObjects->Add(ceTrms);
515 }
516 ceQmean = (AliTPCCalPad*)ceObjects->FindObject("CEQmean");
517 if ( !ceQmean ) {
518 ceQmean = new AliTPCCalPad("CEQmean","CEQmean");
519 ceObjects->Add(ceQmean);
520 }
521
522
523 UInt_t result=0;
524
525 Int_t nSectors = fROC->GetNSectors();
526 TList* list = GetFileSources(sourceFXS,"CE");
527
528 if (list && list->GetEntries()>0) {
529
530// loop through all files from LDCs
531
532 UInt_t index = 0;
533 while (list->At(index)!=NULL) {
534 TObjString* fileNameEntry = (TObjString*) list->At(index);
535 if (fileNameEntry!=NULL) {
536 TString fileName = GetFile(sourceFXS, "CE",
537 fileNameEntry->GetString().Data());
538 TFile *f = TFile::Open(fileName);
539 if (!f) {
540 Log ("Error opening central electrode file.");
541 result =2;
542 break;
543 }
544 AliTPCCalibCE *calCE;
545 f->GetObject("AliTPCCalibCE",calCE);
546
547 // replace entries for the sectors available in the present file
548
549 for (Int_t sector=0; sector<nSectors; sector++) {
550 AliTPCCalROC *rocTmean=calCE->GetCalRocT0(sector);
551 if ( rocTmean ) ceTmean->SetCalROC(rocTmean,sector);
552 AliTPCCalROC *rocTrms=calCE->GetCalRocRMS(sector);
553 if ( rocTrms ) ceTrms->SetCalROC(rocTrms,sector);
554 AliTPCCalROC *rocQmean=calCE->GetCalRocQ(sector);
555 if ( rocQmean ) ceQmean->SetCalROC(rocQmean,sector);
556 }
557 }
558 ++index;
559 } // while(list)
560//
561// Store updated pedestal entry to OCDB
562//
563 AliCDBMetaData metaData;
564 metaData.SetBeamPeriod(0);
565 metaData.SetResponsible("Haavard Helstrup");
566 metaData.SetComment("Preprocessor AliTPC data base entries.");
567
568 Bool_t storeOK = Store("Calib", "CE", ceObjects, &metaData, 0, kTRUE);
569 if ( !storeOK ) ++result;
570
571 } else {
572 Log ("Error: no entries!");
573 result = 1;
574 }
575
576 return result;
577}