typo fixed
[u/mrichter/AliRoot.git] / TRD / AliTRDcalibDB.cxx
CommitLineData
3551db50 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/* $Id$ */
17
18///////////////////////////////////////////////////////////////////////////////
19// //
20// Class providing the calibration parameters by accessing the CDB //
21// //
2745a409 22// Request an instance with AliTRDcalibDB::Instance() //
3551db50 23// If a new event is processed set the event number with SetRun //
24// Then request the calibration data //
25// //
2745a409 26// Author: //
27// Jan Fiete Grosse-Oetringhaus (Jan.Fiete.Grosse-Oetringhaus@cern.ch) //
28// //
3551db50 29///////////////////////////////////////////////////////////////////////////////
30
a7ac01d2 31#include <TClonesArray.h>
e232d349 32#include <TObjArray.h>
3551db50 33
2745a409 34#include "AliCDBManager.h"
2745a409 35#include "AliCDBEntry.h"
36#include "AliLog.h"
3551db50 37
51a0ce25 38#include "AliTRDPIDReference.h"
3551db50 39#include "AliTRDcalibDB.h"
a9e22145 40#include "AliTRDtrapConfig.h"
41#include "AliTRDtrapConfigHandler.h"
eda91732 42#include "AliTRDCommonParam.h"
3551db50 43
7754cd1f 44#include "Cal/AliTRDCalROC.h"
7754cd1f 45#include "Cal/AliTRDCalPad.h"
46#include "Cal/AliTRDCalDet.h"
e232d349 47#include "Cal/AliTRDCalDCS.h"
d6a1ec13 48#include "Cal/AliTRDCalDCSv2.h"
ea3eaa08 49#include "Cal/AliTRDCalDCSFEEv2.h"
720a0a16 50#include "Cal/AliTRDCalPID.h"
7754cd1f 51#include "Cal/AliTRDCalMonitoring.h"
7754cd1f 52#include "Cal/AliTRDCalChamberStatus.h"
7754cd1f 53#include "Cal/AliTRDCalPadStatus.h"
54#include "Cal/AliTRDCalSingleChamberStatus.h"
9dcc64cc 55#include "Cal/AliTRDCalTrkAttach.h"
ea3eaa08 56#include "Cal/AliTRDCalOnlineGainTable.h"
3551db50 57
58ClassImp(AliTRDcalibDB)
59
2745a409 60AliTRDcalibDB *AliTRDcalibDB::fgInstance = 0;
61Bool_t AliTRDcalibDB::fgTerminated = kFALSE;
3551db50 62
63//_ singleton implementation __________________________________________________
64AliTRDcalibDB* AliTRDcalibDB::Instance()
65{
66 //
67 // Singleton implementation
68 // Returns an instance of this class, it is created if neccessary
63a700c6 69 //
3551db50 70
6d50f529 71 if (fgTerminated != kFALSE) {
3551db50 72 return 0;
6d50f529 73 }
3551db50 74
6d50f529 75 if (fgInstance == 0) {
3551db50 76 fgInstance = new AliTRDcalibDB();
6d50f529 77 }
e254cad1 78
3551db50 79 return fgInstance;
4e009ce4 80
3551db50 81}
82
2745a409 83//_____________________________________________________________________________
3551db50 84void AliTRDcalibDB::Terminate()
85{
86 //
87 // Singleton implementation
6bf5f0ce 88 // Deletes the instance of this class and sets the terminated flag,
89 // instances cannot be requested anymore
3551db50 90 // This function can be called several times.
91 //
92
93 fgTerminated = kTRUE;
94
6d50f529 95 if (fgInstance != 0) {
3551db50 96 delete fgInstance;
97 fgInstance = 0;
98 }
2745a409 99
3551db50 100}
101
102//_____________________________________________________________________________
103AliTRDcalibDB::AliTRDcalibDB()
2745a409 104 :TObject()
105 ,fRun(-1)
6d50f529 106 ,fPRFsmp(0)
107 ,fPRFbin(0)
108 ,fPRFlo(0)
109 ,fPRFhi(0)
110 ,fPRFwid(0)
111 ,fPRFpad(0)
9dcc64cc 112 ,fPIDResponse(NULL)
ea3eaa08 113 ,fOnlineGainTableID(0)
a9e22145 114 ,fTrapConfig(0x0)
f16844ef 115 ,fTrapConfigName("")
116 ,fTrapConfigVersion("")
3551db50 117{
118 //
2745a409 119 // Default constructor
3551db50 120 //
2745a409 121 // TODO Default runnumber is set to 0, this should be changed later
122 // to an invalid value (e.g. -1) to prevent
3551db50 123 // TODO invalid calibration data to be used.
2745a409 124 //
125
95867fd1 126 for (Int_t i = 0; i < kCDBCacheSize; ++i) {
2745a409 127 fCDBCache[i] = 0;
128 fCDBEntries[i] = 0;
129 }
3551db50 130
2745a409 131 // Create the sampled PRF
132 SamplePRF();
d4c6453d 133
2745a409 134}
135
136//_____________________________________________________________________________
137AliTRDcalibDB::AliTRDcalibDB(const AliTRDcalibDB &c)
138 :TObject(c)
6d50f529 139 ,fRun(-1)
140 ,fPRFsmp(0)
141 ,fPRFbin(0)
142 ,fPRFlo(0)
143 ,fPRFhi(0)
144 ,fPRFwid(0)
145 ,fPRFpad(0)
9dcc64cc 146 ,fPIDResponse(NULL)
ea3eaa08 147 ,fOnlineGainTableID(0)
a9e22145 148 ,fTrapConfig(0x0)
f16844ef 149 ,fTrapConfigName("")
150 ,fTrapConfigVersion("")
2745a409 151{
152 //
153 // Copy constructor (not that it make any sense for a singleton...)
154 //
155
95867fd1 156 for (Int_t i = 0; i < kCDBCacheSize; ++i) {
2745a409 157 fCDBCache[i] = 0;
3551db50 158 fCDBEntries[i] = 0;
159 }
6a739e92 160
161 // Create the sampled PRF
162 SamplePRF();
2745a409 163
164}
165
166//_____________________________________________________________________________
167AliTRDcalibDB &AliTRDcalibDB::operator=(const AliTRDcalibDB &c)
168{
169 //
170 // Assignment operator (same as above ...)
171 //
172
173 if (this != &c) {
174 AliFatal("No assignment operator defined");
175 }
4e009ce4 176
2745a409 177 return *this;
178
acba9bad 179}
3551db50 180
181//_____________________________________________________________________________
182AliTRDcalibDB::~AliTRDcalibDB()
183{
184 //
185 // destructor
186 //
187
6d50f529 188 if (fPRFsmp) {
189 delete [] fPRFsmp;
190 fPRFsmp = 0;
6a739e92 191 }
ea3eaa08 192
193 if (fPIDResponse) {
194 delete fPIDResponse;
195 fPIDResponse = 0x0;
196 }
6a739e92 197
3551db50 198 Invalidate();
2745a409 199
acba9bad 200}
3551db50 201
63a700c6 202//_caching functions____________________________________________________________
95867fd1 203const TObject *AliTRDcalibDB::GetCachedCDBObject(Int_t id)
63a700c6 204{
2745a409 205 //
206 // Retrieves a cdb object with the given id. The objects are cached as
207 // long as the run number is not changed.
208 //
209 // Put together the available objects here by using the lines
210 // a) For usual calibration objects:
211 // case kID<Name> :
212 // return CacheCDBEntry(kID<Name>,"TRD/Calib/<Path>");
213 // break;
214 // See function CacheCDBEntry for details.
215 // and
216 // b) For calibration data which depends on two objects: One containing
217 // a value per detector and one the local fluctuations per pad:
218 // case kID<Name> :
219 // return CacheMergeCDBEntry(kID<Name>,"TRD/Calib/<padPath>","TRD/Calib/<chamberPath>");
220 // break;
221 // See function CacheMergeCDBEntry for details.
222 //
63a700c6 223
6d50f529 224 switch (id) {
2745a409 225
226 // Parameters defined per pad and chamber
227 case kIDVdriftPad :
6d50f529 228 return CacheCDBEntry(kIDVdriftPad ,"TRD/Calib/LocalVdrift");
2745a409 229 break;
230 case kIDVdriftChamber :
6d50f529 231 return CacheCDBEntry(kIDVdriftChamber ,"TRD/Calib/ChamberVdrift");
2745a409 232 break;
eda91732 233 case kIDExBChamber :
234 return CacheCDBEntry(kIDExBChamber ,"TRD/Calib/ChamberExB");
235 break;
2745a409 236 case kIDT0Pad :
6d50f529 237 return CacheCDBEntry(kIDT0Pad ,"TRD/Calib/LocalT0");
2745a409 238 break;
239 case kIDT0Chamber :
6d50f529 240 return CacheCDBEntry(kIDT0Chamber ,"TRD/Calib/ChamberT0");
2745a409 241 break;
242 case kIDGainFactorPad :
6d50f529 243 return CacheCDBEntry(kIDGainFactorPad ,"TRD/Calib/LocalGainFactor");
2745a409 244 break;
245 case kIDGainFactorChamber :
6d50f529 246 return CacheCDBEntry(kIDGainFactorChamber ,"TRD/Calib/ChamberGainFactor");
2745a409 247 break;
ea3eaa08 248
249 case kIDOnlineGainFactor :
250 switch(GetOnlineGainTableID()) {
251 case 0:
252 // For testing purposes only !!!
253 AliInfo("No gain table name from OCDB. Use default table!");
254 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Krypton_2011-01");
255 break;
256 case 1:
257 // Online gain table ID 1
258 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Krypton_2011-01");
259 break;
412d07b9 260 case 2:
261 // Online gain table ID 2
262 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Gaintbl_Uniform_FGAN0_2011-01");
263 break;
264 case 3:
265 // Online gain table ID 3
266 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Gaintbl_Uniform_FGAN8_2011-01");
267 break;
3abf1c06 268 case 4:
269 // Online gain table ID 4
270 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Krypton_2011-02");
271 break;
4a05be5f 272 case 5:
273 // Online gain table ID 5
274 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Krypton_2011-03");
275 break;
276 case 6:
277 // Online gain table ID 6
278 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Gaintbl_Uniform_FGAN0_2012-01");
279 break;
280 case 7:
281 // Online gain table ID 7
282 return CacheCDBEntry(kIDOnlineGainFactor ,"TRD/Calib/Gaintbl_Uniform_FGAN8_2012-01");
283 break;
ea3eaa08 284 }
285 break;
286
df83a620 287 case kIDNoiseChamber :
ea3eaa08 288 return CacheCDBEntry(kIDNoiseChamber ,"TRD/Calib/DetNoise");
df83a620 289 break;
290 case kIDNoisePad :
ea3eaa08 291 return CacheCDBEntry(kIDNoisePad ,"TRD/Calib/PadNoise");
df83a620 292 break;
293
2745a409 294 // Parameters defined per pad
295 case kIDPRFWidth :
6d50f529 296 return CacheCDBEntry(kIDPRFWidth ,"TRD/Calib/PRFWidth");
2745a409 297 break;
298
299 // Status values
2745a409 300 case kIDChamberStatus :
6d50f529 301 return CacheCDBEntry(kIDChamberStatus ,"TRD/Calib/ChamberStatus");
2745a409 302 break;
2745a409 303 case kIDPadStatus :
6d50f529 304 return CacheCDBEntry(kIDPadStatus ,"TRD/Calib/PadStatus");
2745a409 305 break;
306
307 // Global parameters
308 case kIDMonitoringData :
6d50f529 309 return CacheCDBEntry(kIDMonitoringData ,"TRD/Calib/MonitoringData");
2745a409 310 break;
aa617684 311 case kIDFEE :
312 return CacheCDBEntry(kIDFEE ,"TRD/Calib/FEE");
2745a409 313 break;
a9e22145 314 case kIDTrapConfig :
315 return CacheCDBEntry(kIDFEE ,"TRD/Calib/TrapConfig");
316 break;
e232d349 317 case kIDDCS :
5f481838 318 return CacheCDBEntry(kIDDCS ,"TRD/Calib/DCS");
e232d349 319 break;
44dbae42 320 case kIDPIDNN :
321 return CacheCDBEntry(kIDPIDNN ,"TRD/Calib/PIDNN");
ea3eaa08 322 break;
2745a409 323 case kIDPIDLQ :
6d50f529 324 return CacheCDBEntry(kIDPIDLQ ,"TRD/Calib/PIDLQ");
2745a409 325 break;
9dcc64cc 326 case kIDPIDLQ1D:
327 return CacheCDBEntry(kIDPIDLQ1D ,"TRD/Calib/PIDLQ1D");
328 break;
a7ac01d2 329 case kIDRecoParam :
9dcc64cc 330 return CacheCDBEntry(kIDRecoParam ,"TRD/Calib/RecoParam");
331 break;
332 case kIDAttach :
333 return CacheCDBEntry(kIDAttach ,"TRD/Calib/TrkAttach");
a7ac01d2 334 break;
e2a1c98b 335 case kIDPHQ :
336 return CacheCDBEntry(kIDPHQ ,"TRD/Calib/PHQ");
337 break;
2745a409 338 }
339
340 return 0;
341
63a700c6 342}
343
344//_____________________________________________________________________________
95867fd1 345AliCDBEntry *AliTRDcalibDB::GetCDBEntry(const char *cdbPath)
63a700c6 346{
347 //
348 // Retrieves an entry with path <cdbPath> from the CDB.
349 //
350
95867fd1 351 AliCDBEntry *entry = AliCDBManager::Instance()->Get(cdbPath,fRun);
2745a409 352 if (!entry) {
a7ac01d2 353 AliError(Form("Failed to get entry: %s",cdbPath));
63a700c6 354 return 0;
355 }
356
63a700c6 357 return entry;
2745a409 358
63a700c6 359}
360
361//_____________________________________________________________________________
95867fd1 362const TObject *AliTRDcalibDB::CacheCDBEntry(Int_t id, const char *cdbPath)
63a700c6 363{
364 //
365 // Caches the entry <id> with cdb path <cdbPath>
366 //
d6a1ec13 367
2745a409 368 if (!fCDBCache[id]) {
63a700c6 369 fCDBEntries[id] = GetCDBEntry(cdbPath);
95867fd1 370 if (fCDBEntries[id]) {
63a700c6 371 fCDBCache[id] = fCDBEntries[id]->GetObject();
95867fd1 372 }
d6a1ec13 373 }
374
63a700c6 375 return fCDBCache[id];
ea3eaa08 376
63a700c6 377}
378
3551db50 379//_____________________________________________________________________________
380void AliTRDcalibDB::SetRun(Long64_t run)
381{
382 //
7754cd1f 383 // Sets current run number. Calibration data is read from the corresponding file.
3551db50 384 // When the run number changes the caching is invalidated.
385 //
7754cd1f 386
6d50f529 387 if (fRun == run) {
3551db50 388 return;
6d50f529 389 }
7754cd1f 390
3551db50 391 fRun = run;
95867fd1 392
3551db50 393 Invalidate();
2745a409 394
3551db50 395}
7754cd1f 396
3551db50 397//_____________________________________________________________________________
398void AliTRDcalibDB::Invalidate()
399{
400 //
401 // Invalidates cache (when run number is changed).
402 //
403
95867fd1 404 for (Int_t i = 0; i < kCDBCacheSize; ++i) {
2745a409 405 if (fCDBEntries[i]) {
406 if (AliCDBManager::Instance()->GetCacheFlag() == kFALSE) {
11f5dd43 407 if ((fCDBEntries[i]->IsOwner() == kFALSE) &&
408 (fCDBCache[i])) {
e254cad1 409 delete fCDBCache[i];
2745a409 410 }
e254cad1 411 delete fCDBEntries[i];
412 }
3551db50 413 fCDBEntries[i] = 0;
2745a409 414 fCDBCache[i] = 0;
3551db50 415 }
416 }
2745a409 417
ea3eaa08 418 fOnlineGainTableID = 0;
419
df83a620 420}
ea3eaa08 421
df83a620 422//_____________________________________________________________________________
423Float_t AliTRDcalibDB::GetNoise(Int_t det, Int_t col, Int_t row)
424{
425 //
426 // Returns the noise level in ADC counts for the given pad.
427 //
428
429 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
430 (GetCachedCDBObject(kIDNoisePad));
431 if (!calPad) {
432 return -1;
433 }
434
435 AliTRDCalROC *roc = calPad->GetCalROC(det);
436 if (!roc) {
437 return -1;
438 }
439
440 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
441 (GetCachedCDBObject(kIDNoiseChamber));
442 if (!calChamber) {
443 return -1;
444 }
445
446 return calChamber->GetValue(det) * roc->GetValue(col,row);
447
448}
449
450//_____________________________________________________________________________
451AliTRDCalROC *AliTRDcalibDB::GetNoiseROC(Int_t det)
452{
453 //
454 // Returns the Vdrift calibration object for a given ROC
455 // containing one number per pad
456 //
457
458 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
459 (GetCachedCDBObject(kIDNoisePad));
460 if (!calPad) {
461 return 0;
462 }
463
464 AliTRDCalROC *roc = calPad->GetCalROC(det);
465 if (!roc) {
466 return 0;
467 }
468 else {
469 return roc;
470 }
471
472}
473
474//_____________________________________________________________________________
475const AliTRDCalDet *AliTRDcalibDB::GetNoiseDet()
476{
477 //
478 // Returns the Vdrift calibration object
479 // containing one number per detector
480 //
481
482 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
483 (GetCachedCDBObject(kIDNoiseChamber));
484 if (!calChamber) {
485 return 0;
486 }
487 else {
488 return calChamber;
489 }
490
3551db50 491}
492
3551db50 493//_____________________________________________________________________________
494Float_t AliTRDcalibDB::GetVdrift(Int_t det, Int_t col, Int_t row)
495{
496 //
497 // Returns the drift velocity for the given pad.
498 //
7754cd1f 499
95867fd1 500 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
6d50f529 501 (GetCachedCDBObject(kIDVdriftPad));
502 if (!calPad) {
3551db50 503 return -1;
6d50f529 504 }
3551db50 505
95867fd1 506 AliTRDCalROC *roc = calPad->GetCalROC(det);
6d50f529 507 if (!roc) {
3551db50 508 return -1;
6d50f529 509 }
3551db50 510
95867fd1 511 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
6d50f529 512 (GetCachedCDBObject(kIDVdriftChamber));
513 if (!calChamber) {
7754cd1f 514 return -1;
6d50f529 515 }
7754cd1f 516
95867fd1 517 return calChamber->GetValue(det) * roc->GetValue(col,row);
2745a409 518
7754cd1f 519}
afb9f880 520
521//_____________________________________________________________________________
522AliTRDCalROC *AliTRDcalibDB::GetVdriftROC(Int_t det)
523{
524 //
525 // Returns the Vdrift calibration object for a given ROC
526 // containing one number per pad
527 //
528
529 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
530 (GetCachedCDBObject(kIDVdriftPad));
531 if (!calPad) {
532 return 0;
533 }
534
535 AliTRDCalROC *roc = calPad->GetCalROC(det);
536 if (!roc) {
537 return 0;
538 }
539 else {
540 return roc;
541 }
542
543}
544
545//_____________________________________________________________________________
546const AliTRDCalDet *AliTRDcalibDB::GetVdriftDet()
547{
548 //
549 // Returns the Vdrift calibration object
550 // containing one number per detector
551 //
552
553 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
554 (GetCachedCDBObject(kIDVdriftChamber));
555 if (!calChamber) {
556 return 0;
557 }
558 else {
559 return calChamber;
560 }
561
562}
7754cd1f 563
e2a1c98b 564//_____________________________________________________________________________
565TObjArray * AliTRDcalibDB::GetPHQ()
566{
567 //
568 //return PHQ calibration object
569 //
570 TObjArray *arr = (TObjArray *) (GetCachedCDBObject(kIDPHQ));
571 return arr;
572}
573
7754cd1f 574//_____________________________________________________________________________
575Float_t AliTRDcalibDB::GetVdriftAverage(Int_t det)
576{
577 //
578 // Returns the average drift velocity for the given detector
579 //
580
95867fd1 581 const AliTRDCalDet *calDet = dynamic_cast<const AliTRDCalDet *>
6d50f529 582 (GetCachedCDBObject(kIDVdriftChamber));
583 if (!calDet) {
7754cd1f 584 return -1;
6d50f529 585 }
7754cd1f 586
587 return calDet->GetValue(det);
2745a409 588
acba9bad 589}
eda91732 590//_____________________________________________________________________________
591const AliTRDCalDet *AliTRDcalibDB::GetExBDet()
592{
593 //
594 // Returns the exB calibration object
595 // containing one number per detector
596 //
597
598 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *> (GetCachedCDBObject(kIDExBChamber));
599
600 Double_t meanexb = 100.0;
601 if (calChamber) meanexb = calChamber->GetMean();
602 //printf("mean %f\n",meanexb);
3551db50 603
eda91732 604 if ((!calChamber) || (meanexb > 70.0)) {
605
606 const AliTRDCalDet *calChambervdrift = dynamic_cast<const AliTRDCalDet *>
607 (GetCachedCDBObject(kIDVdriftChamber));
608 if (!calChambervdrift) {
609 return 0;
610 }
611 else {
612 AliTRDCalDet *calDetExB = new AliTRDCalDet("lorentz angle tan","lorentz angle tan (detector value)");
613 for(Int_t k = 0; k < 540; k++){
614 calDetExB->SetValue(k,AliTRDCommonParam::Instance()->GetOmegaTau(calChambervdrift->GetValue(k)));
eda91732 615 }
616 return calDetExB;
617 }
618 }
619 else return calChamber;
620
621}
3551db50 622//_____________________________________________________________________________
623Float_t AliTRDcalibDB::GetT0(Int_t det, Int_t col, Int_t row)
624{
625 //
626 // Returns t0 for the given pad.
627 //
628
d4c6453d 629 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
630 (GetCachedCDBObject(kIDT0Pad));
6d50f529 631 if (!calPad) {
3551db50 632 return -1;
6d50f529 633 }
3551db50 634
d4c6453d 635 AliTRDCalROC *roc = calPad->GetCalROC(det);
6d50f529 636 if (!roc) {
3551db50 637 return -1;
6d50f529 638 }
3551db50 639
d4c6453d 640 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
641 (GetCachedCDBObject(kIDT0Chamber));
6d50f529 642 if (!calChamber) {
7754cd1f 643 return -1;
6d50f529 644 }
7754cd1f 645
8d786fb8 646 return calChamber->GetValue(det) + roc->GetValue(col,row);
2745a409 647
7754cd1f 648}
afb9f880 649
56178ff4 650//_____________________________________________________________________________
651AliTRDCalROC *AliTRDcalibDB::GetT0ROC(Int_t det)
652{
653 //
654 // Returns the t0 calibration object for a given ROC
655 // containing one number per pad
656 //
657
658 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
659 (GetCachedCDBObject(kIDT0Pad));
660 if (!calPad) {
661 return 0;
662 }
663
664 AliTRDCalROC *roc = calPad->GetCalROC(det);
665 if (!roc) {
666 return 0;
667 }
668 else {
669 return roc;
670 }
671
672}
673
674//_____________________________________________________________________________
675const AliTRDCalDet *AliTRDcalibDB::GetT0Det()
676{
677 //
678 // Returns the t0 calibration object
679 // containing one number per detector
680 //
681
682 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
683 (GetCachedCDBObject(kIDT0Chamber));
684 if (!calChamber) {
685 return 0;
686 }
687 else {
688 return calChamber;
689 }
690
691}
692
7754cd1f 693//_____________________________________________________________________________
694Float_t AliTRDcalibDB::GetT0Average(Int_t det)
695{
696 //
697 // Returns the average t0 for the given detector
698 //
699
5aba4f09 700 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
701 (GetCachedCDBObject(kIDT0Pad));
702 if (!calPad) {
703 return -1;
704 }
705
706 AliTRDCalROC *roc = calPad->GetCalROC(det);
707 if (!roc) {
708 return -1;
709 }
710
95867fd1 711 const AliTRDCalDet *calDet = dynamic_cast<const AliTRDCalDet *>
6d50f529 712 (GetCachedCDBObject(kIDT0Chamber));
713 if (!calDet) {
7754cd1f 714 return -1;
6d50f529 715 }
7754cd1f 716
4eb02ab1 717 Double_t sum = 0.0;
5aba4f09 718 for (Int_t channel = 0; channel < roc->GetNchannels(); ++channel) {
4eb02ab1 719 sum += roc->GetValue(channel);
5aba4f09 720 }
4eb02ab1 721 sum /= roc->GetNchannels();
722 sum += calDet->GetValue(det);
723 return sum;
2745a409 724
acba9bad 725}
3551db50 726
ea3eaa08 727//_____________________________________________________________________________
728AliTRDCalOnlineGainTableROC* AliTRDcalibDB::GetOnlineGainTableROC(Int_t det)
729{
730 //
731 // Returns the online gain factor table for a given ROC.
732 //
c36ddf52 733
734 if (!HasOnlineFilterGain()) {
735 return 0x0;
736 }
ea3eaa08 737
738 const AliTRDCalOnlineGainTable *calOnline
739 = dynamic_cast<const AliTRDCalOnlineGainTable *>
740 (GetCachedCDBObject(kIDOnlineGainFactor));
741 if (!calOnline) {
742 return 0x0;
743 }
744
745 return calOnline->GetGainTableROC(det);
746
747}
748
749//_____________________________________________________________________________
750Float_t AliTRDcalibDB::GetOnlineGainFactor(Int_t det, Int_t col, Int_t row)
751{
752 //
753 // Returns the online gain factor for the given pad.
754 //
c36ddf52 755
756 if (!HasOnlineFilterGain()) {
757 return 0x0;
758 }
ea3eaa08 759
760 const AliTRDCalOnlineGainTable *calOnline
761 = dynamic_cast<const AliTRDCalOnlineGainTable *>
762 (GetCachedCDBObject(kIDOnlineGainFactor));
763 if (!calOnline) {
764 return -1;
765 }
766
767 return calOnline->GetGainCorrectionFactor(det,row,col);
768
769}
770
56178ff4 771//_____________________________________________________________________________
772AliTRDCalROC *AliTRDcalibDB::GetGainFactorROC(Int_t det)
773{
774 //
775 // Returns the gain factor calibration object for a given ROC
56178ff4 776 //
777
778 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
779 (GetCachedCDBObject(kIDGainFactorPad));
780 if (!calPad) {
781 return 0;
782 }
783
784 AliTRDCalROC *roc = calPad->GetCalROC(det);
785 if (!roc) {
786 return 0;
787 }
788 else {
789 return roc;
790 }
791
792}
793
c36ddf52 794//_____________________________________________________________________________
795Float_t AliTRDcalibDB::GetGainFactor(Int_t det, Int_t col, Int_t row)
796{
797 //
798 // Returns the gain factor for the given pad.
799 //
800
801 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
802 (GetCachedCDBObject(kIDGainFactorPad));
803 if (!calPad) {
804 return -1;
805 }
806
807 AliTRDCalROC *roc = calPad->GetCalROC(det);
808 if (!roc) {
809 return -1;
810 }
811
812 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
813 (GetCachedCDBObject(kIDGainFactorChamber));
814 if (!calChamber) {
815 return -1;
816 }
817
818 return calChamber->GetValue(det) * roc->GetValue(col,row);
819
820}
821
56178ff4 822//_____________________________________________________________________________
823const AliTRDCalDet *AliTRDcalibDB::GetGainFactorDet()
824{
825 //
826 // Returns the gain factor calibration object
827 // containing one number per detector
828 //
829
830 const AliTRDCalDet *calChamber = dynamic_cast<const AliTRDCalDet *>
831 (GetCachedCDBObject(kIDGainFactorChamber));
832 if (!calChamber) {
833 return 0;
834 }
835 else {
836 return calChamber;
837 }
838
839}
840
7754cd1f 841//_____________________________________________________________________________
842Float_t AliTRDcalibDB::GetGainFactorAverage(Int_t det)
843{
844 //
845 // Returns the average gain factor for the given detector
846 //
847
95867fd1 848 const AliTRDCalDet *calDet = dynamic_cast<const AliTRDCalDet *>
6d50f529 849 (GetCachedCDBObject(kIDGainFactorChamber));
850 if (!calDet) {
7754cd1f 851 return -1;
6d50f529 852 }
7754cd1f 853
854 return calDet->GetValue(det);
2745a409 855
acba9bad 856}
6a739e92 857
524fc8fa 858//_____________________________________________________________________________
859AliTRDCalROC *AliTRDcalibDB::GetPRFROC(Int_t det)
860{
861 //
862 // Returns the PRF calibration object for a given ROC
863 // containing one number per pad
864 //
865
866 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
867 (GetCachedCDBObject(kIDPRFWidth));
868 if (!calPad) {
869 return 0;
870 }
871
872 AliTRDCalROC *roc = calPad->GetCalROC(det);
873 if (!roc) {
874 return 0;
875 }
876 else {
877 return roc;
878 }
879
880}
881
6a739e92 882//_____________________________________________________________________________
883Float_t AliTRDcalibDB::GetPRFWidth(Int_t det, Int_t col, Int_t row)
884{
885 //
886 // Returns the PRF width for the given pad.
887 //
888
95867fd1 889 const AliTRDCalPad *calPad = dynamic_cast<const AliTRDCalPad *>
6d50f529 890 (GetCachedCDBObject(kIDPRFWidth));
891 if (!calPad) {
6a739e92 892 return -1;
6d50f529 893 }
6a739e92 894
95867fd1 895 AliTRDCalROC *roc = calPad->GetCalROC(det);
6d50f529 896 if (!roc) {
6a739e92 897 return -1;
6d50f529 898 }
6a739e92 899
95867fd1 900 return roc->GetValue(col,row);
2745a409 901
acba9bad 902}
3551db50 903
904//_____________________________________________________________________________
f1dcad37 905Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
3551db50 906{
e232d349 907 //
a835b5d7 908 // Returns number of time bins from the DCS
e232d349 909 //
f1dcad37 910
a835b5d7 911 // default value - has not been set
912 Int_t nUndef = -1;
c6f7c6cb 913
a835b5d7 914 // Get the corresponding parameter
915 TString tbstr = "";
916 GetDCSConfigParOption(kTimebin, 0, tbstr);
917
918 // Check if there is any content in the string first
919 if (tbstr.Length() == 0) {
920 AliError("Parameter for number of timebins is empty!");
c6f7c6cb 921 return nUndef;
e232d349 922 }
c6f7c6cb 923
a835b5d7 924 // Check if we have the correct config parameter
925 TString tbident = "tb";
926 TString tbsubstr = tbstr(0,2);
927 if (!tbsubstr.EqualTo(tbident)) {
928 AliError(Form("Parameter for number of timebins is corrupted (%s)!", tbstr.Data()));
929 return nUndef;
930 }
c6f7c6cb 931
a835b5d7 932 tbstr.Remove(0,2);
933 // check if there is more than a number
934 if (!tbstr.IsDigit()) {
935 AliError(Form("Parameter for number of timebins is corrupted (%s)!", tbstr.Data()));
936 return nUndef;
e232d349 937 }
d6a1ec13 938
a835b5d7 939 Int_t ntb = tbstr.Atoi();
940 AliInfo(Form("Number of timebins from CDB: %d", ntb));
c6f7c6cb 941
a835b5d7 942 return ntb;
ea3eaa08 943
e232d349 944}
945
946//_____________________________________________________________________________
f1dcad37 947void AliTRDcalibDB::GetFilterType(TString &filterType)
948{
949 //
950 // Returns the filter type
951 //
952
fc358d0f 953 GetDCSConfigParOption(kFltrSet, 0, filterType);
d6a1ec13 954
af086c5e 955}
956
957//_____________________________________________________________________________
ea3eaa08 958Int_t AliTRDcalibDB::GetOnlineGainTableID()
959{
960 //
961 // Get the gain table ID from the DCS
962 //
963
964 if (fOnlineGainTableID > 0) {
965 return fOnlineGainTableID;
966 }
967
968 const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
969 if (!dcsArr){
970 return -1;
971 }
972
973 Int_t esor = 0; // Take SOR
974 Int_t calver = 0; // Check CalDCS version
975 if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
976 if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
977
978 if (calver == 1) {
979
980 // No data for old DCS object available, anyway
981 return -1;
982
983 }
984 else if (calver == 2) {
985
986 // DCSv2 object
987 const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
988 if(!calDCSv2){
989 return -1;
990 }
991
992 TString tableName = "";
993 for (Int_t i = 0; i < 540; i++) {
994 const AliTRDCalDCSFEEv2 *calDCSFEEv2 = calDCSv2->GetCalDCSFEEObj(0);
995 tableName = calDCSFEEv2->GetGainTableName();
996 if (tableName.Length() > 0) {
997 break;
998 }
999 }
412d07b9 1000 if (tableName.CompareTo("Krypton_2011-01") == 0) {
ea3eaa08 1001 fOnlineGainTableID = 1;
1002 return fOnlineGainTableID;
1003 }
412d07b9 1004 if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2011-01") == 0) {
1005 fOnlineGainTableID = 2;
1006 return fOnlineGainTableID;
1007 }
1008 if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2011-01") == 0) {
1009 fOnlineGainTableID = 3;
1010 return fOnlineGainTableID;
1011 }
3abf1c06 1012 if (tableName.CompareTo("Krypton_2011-02") == 0) {
1013 fOnlineGainTableID = 4;
1014 return fOnlineGainTableID;
1015 }
4a05be5f 1016 if (tableName.CompareTo("Krypton_2011-03") == 0) {
1017 fOnlineGainTableID = 5;
1018 return fOnlineGainTableID;
1019 }
1020 if (tableName.CompareTo("Gaintbl_Uniform_FGAN0_2012-01") == 0) {
1021 fOnlineGainTableID = 6;
1022 return fOnlineGainTableID;
1023 }
1024 if (tableName.CompareTo("Gaintbl_Uniform_FGAN8_2012-01") == 0) {
1025 fOnlineGainTableID = 7;
1026 return fOnlineGainTableID;
1027 }
ea3eaa08 1028
1029 }
1030 else {
1031
1032 AliError("NO DCS/DCSv2 OCDB entry found!");
1033 return -1;
1034
1035 }
1036
45b38fb0 1037 return -1;
1038
ea3eaa08 1039}
1040
1041//_____________________________________________________________________________
1042void AliTRDcalibDB::GetGlobalConfiguration(TString &config)
1043{
af086c5e 1044 //
1045 // Get Configuration from the DCS
1046 //
d6a1ec13 1047
af086c5e 1048 const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
1049 if(!dcsArr){
a835b5d7 1050 AliError("No DCS CDB Object available!");
af086c5e 1051 config = "";
1052 return;
1053 }
d6a1ec13 1054
a835b5d7 1055 Int_t idSOR = 0, idEOR=1; // The index of SOR and EOR
1056 Bool_t hasSOR = (dcsArr->At(idSOR));
1057 Bool_t hasEOR = (dcsArr->At(idEOR));
1058 TString cfgSOR = "", cfgEOR = ""; // The configuration at SOR/EOR
d6a1ec13 1059
a835b5d7 1060 // The SOR object is mandatory
1061 if (!hasSOR) {
1062 AliError("NO SOR object found in CDB file!");
1063 config = "";
1064 return;
1065 }
1066 if (!hasEOR) AliWarning("NO EOR object found in CDB file!");
ea3eaa08 1067
a835b5d7 1068 // Check CalDCS version
1069 Int_t calver = 0;
1070 if (!strcmp(dcsArr->At(idSOR)->ClassName(),"AliTRDCalDCS")) calver = 1;
1071 else if (!strcmp(dcsArr->At(idSOR)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
ea3eaa08 1072
a835b5d7 1073 // Get the configuration strings
1074 if (calver == 1) {
1075 // DCS object
1076 const AliTRDCalDCS *calSOR = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(idSOR));
1077 cfgSOR = calSOR->GetGlobalConfigName();
1078 if (hasEOR) {
1079 const AliTRDCalDCS *calEOR = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(idEOR));
1080 cfgEOR = calEOR->GetGlobalConfigName();
1081 }
ea3eaa08 1082 }
1083 else if (calver == 2) {
d6a1ec13 1084 // DCSv2 object
a835b5d7 1085 const AliTRDCalDCSv2 *calv2SOR = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(idSOR));
1086 cfgSOR = calv2SOR->GetGlobalConfigName();
1087 if (hasEOR) {
1088 const AliTRDCalDCSv2 *calv2EOR = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(idEOR));
1089 cfgEOR = calv2EOR->GetGlobalConfigName();
1090 }
ea3eaa08 1091 }
1092 else {
ea3eaa08 1093 AliError("NO DCS/DCSv2 OCDB entry found!");
a835b5d7 1094 config = "";
1095 return;
1096 }
ea3eaa08 1097
a835b5d7 1098 // If there is no EOR entry, return the SOR value
1099 if (!hasEOR || cfgEOR.Length()==0) {
1100 config = cfgSOR;
1101 return;
ea3eaa08 1102 }
d6a1ec13 1103
a835b5d7 1104 // Check if the configuration is the same for both
1105 if (cfgSOR.EqualTo(cfgEOR)) {
1106 config = cfgSOR;
1107 return;
1108 }
1109
1110 // When both SOR and EOR have an entry but are different, the config is not defined
1111 AliError("Inconsistent configuration at start and end of run found!");
1112 config = "";
1113 return;
1114
af086c5e 1115}
f1dcad37 1116
f8059b39 1117//_____________________________________________________________________________
1118void AliTRDcalibDB::GetGlobalConfigurationVersion(TString &version)
1119{
1120 //
1121 // Get Version of Configuration from the DCS
1122 //
1123
1124 const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
1125 if(!dcsArr){
1126 version = "";
1127 return;
1128 }
1129
1130 Int_t esor = 0; // Take SOR
1131 Int_t calver = 0; // Check CalDCS version
1132 if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
1133 if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
1134
1135 if (calver == 1) {
1136
1137 // DCS object
1138 const AliTRDCalDCS *calDCS = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(esor));
1139 if(!calDCS){
1140 version = "";
1141 return;
1142 }
1143 version = calDCS->GetGlobalConfigVersion();
1144
1145 }
1146 else if (calver == 2) {
1147
1148 // DCSv2 object
1149 const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
1150 if(!calDCSv2){
1151 version = "";
1152 return;
1153 }
1154 version = calDCSv2->GetGlobalConfigVersion();
1155
1156 }
1157 else {
1158
1159 AliError("NO DCS/DCSv2 OCDB entry found!");
1160
1161 }
1162
1163}
1164
af086c5e 1165//_____________________________________________________________________________
fc358d0f 1166void AliTRDcalibDB::GetDCSConfigParOption(Int_t cfgType, Int_t option, TString &cfgo)
af086c5e 1167{
1168 //
fc358d0f 1169 // Get a configuration (see enum in header file) or the options of a configuration
1170 // option == 0 returns the configuration itself
1171 // option > 0 returns the optional parameter Nr. (option) of the configuration (cfgType)
af086c5e 1172 //
ea3eaa08 1173
fc358d0f 1174 // define the delimiters
1175 TString cdelim = "_";
1176 TString odelim = "-";
ea3eaa08 1177
fc358d0f 1178 // get the full configuration name
1179 TString cname;
af086c5e 1180 GetGlobalConfiguration(cname);
fc358d0f 1181 TObjArray *carr = cname.Tokenize(cdelim);
1182 Int_t nconfig = carr->GetEntries();
1183
1184 // protect
1185 if (nconfig == 0) {
1186 AliError("Bad DCS configuration name!");
1187 cfgo = "";
1188 return;
1189 } else if ((nconfig-1) < cfgType) {
1190 AliError("Not enough DCS configuration parameters!");
1191 cfgo = "";
1192 return;
1193 }
1194
1195 TString fullcfg = ((TObjString*)carr->At(cfgType))->GetString();
1196
1197 if (fullcfg.Contains(odelim)) {
ea3eaa08 1198
fc358d0f 1199 TObjArray *oarr = fullcfg.Tokenize(odelim);
1200 Int_t noptions = oarr->GetEntries();
1201
1202 // protect
1203 if ((noptions-1) < option) {
1204 AliError("Not enough DCS configuration options defined!");
1205 cfgo = "";
1206 return;
1207 }
1208
1209 cfgo = ((TObjString*)oarr->At(option))->GetString();
1210 return;
1211
1212 }
1213 else {
1214
1215 if (option != 0) {
1216 AliError("Not enough DCS configuration options defined!");
1217 cfgo = "";
1218 return;
1219 }
1220 cfgo = fullcfg;
1221 return;
1222
1223 }
1224
1225}
1226
1227//_____________________________________________________________________________
1228Bool_t AliTRDcalibDB::HasOnlineFilterPedestal()
1229{
1230 //
1231 // Checks whether pedestal filter was applied online
1232 //
1233
1234 TString filterconfig;
1235 GetFilterType(filterconfig);
ea3eaa08 1236
af086c5e 1237 return filterconfig.Contains("p");
ea3eaa08 1238
af086c5e 1239}
1240
1241//_____________________________________________________________________________
ea3eaa08 1242Bool_t AliTRDcalibDB::HasOnlineFilterGain()
1243{
af086c5e 1244 //
1245 // Checks whether online gain filter was applied
1246 //
ea3eaa08 1247
fc358d0f 1248 TString filterconfig;
1249 GetFilterType(filterconfig);
ea3eaa08 1250
af086c5e 1251 return filterconfig.Contains("g");
ea3eaa08 1252
af086c5e 1253}
1254
1255//_____________________________________________________________________________
ea3eaa08 1256Bool_t AliTRDcalibDB::HasOnlineTailCancellation()
1257{
af086c5e 1258 //
1259 // Checks whether online tail cancellation was applied
1260 //
ea3eaa08 1261
fc358d0f 1262 TString filterconfig;
1263 GetFilterType(filterconfig);
ea3eaa08 1264
af086c5e 1265 return filterconfig.Contains("t");
ea3eaa08 1266
e232d349 1267}
1268
7754cd1f 1269//_____________________________________________________________________________
1270Char_t AliTRDcalibDB::GetPadStatus(Int_t det, Int_t col, Int_t row)
1271{
1272 //
1273 // Returns the status of the given pad
1274 //
1275
aa617684 1276 const AliTRDCalPadStatus *cal = dynamic_cast<const AliTRDCalPadStatus *>
1277 (GetCachedCDBObject(kIDPadStatus));
6d50f529 1278 if (!cal) {
7754cd1f 1279 return -1;
6d50f529 1280 }
7754cd1f 1281
95867fd1 1282 const AliTRDCalSingleChamberStatus *roc = cal->GetCalROC(det);
6d50f529 1283 if (!roc) {
7754cd1f 1284 return -1;
6d50f529 1285 }
7754cd1f 1286
95867fd1 1287 return roc->GetStatus(col,row);
2745a409 1288
7754cd1f 1289}
1290
0e09df31 1291//_____________________________________________________________________________
1292AliTRDCalSingleChamberStatus* AliTRDcalibDB::GetPadStatusROC(Int_t det)
1293{
1294 //
1295 // Returns the pad status calibration object for a given ROC
1296 //
1297
1298 const AliTRDCalPadStatus *cal = dynamic_cast<const AliTRDCalPadStatus *>
1299 (GetCachedCDBObject(kIDPadStatus));
1300 if (!cal) {
1301 return 0;
1302 }
1303
1304 AliTRDCalSingleChamberStatus *roc = cal->GetCalROC(det);
1305 if (!roc) {
1306 return 0;
1307 }
1308 else {
1309 return roc;
1310 }
1311
1312}
1313
7754cd1f 1314//_____________________________________________________________________________
1315Char_t AliTRDcalibDB::GetChamberStatus(Int_t det)
1316{
1317 //
1318 // Returns the status of the given chamber
1319 //
1320
95867fd1 1321 const AliTRDCalChamberStatus *cal = dynamic_cast<const AliTRDCalChamberStatus *>
6d50f529 1322 (GetCachedCDBObject(kIDChamberStatus));
1323 if (!cal) {
7754cd1f 1324 return -1;
6d50f529 1325 }
7754cd1f 1326
1327 return cal->GetStatus(det);
2745a409 1328
7754cd1f 1329}
1330
a7ac01d2 1331//_____________________________________________________________________________
1332AliTRDrecoParam* AliTRDcalibDB::GetRecoParam(Int_t */*eventtype*/)
1333{
be28462a 1334 //
1335 // Returns the TRD reconstruction parameters from the OCDB
1336 //
1337
a7ac01d2 1338 const TClonesArray *recos = dynamic_cast<const TClonesArray*>(GetCachedCDBObject(kIDRecoParam));
be28462a 1339 if (!recos) return 0x0;
a7ac01d2 1340
1341 // calculate entry based on event type info
1342 Int_t n = 0; //f(eventtype[0], eventtype[1], ....)
be28462a 1343
ea3eaa08 1344 return (AliTRDrecoParam *) recos->UncheckedAt(n);
a7ac01d2 1345
be28462a 1346}
a7ac01d2 1347
7754cd1f 1348//_____________________________________________________________________________
1349Bool_t AliTRDcalibDB::IsPadMasked(Int_t det, Int_t col, Int_t row)
1350{
1351 //
1352 // Returns status, see name of functions for details ;-)
1353 //
1354
95867fd1 1355 const AliTRDCalPadStatus *cal = dynamic_cast<const AliTRDCalPadStatus *>
1356 (GetCachedCDBObject(kIDPadStatus));
6d50f529 1357 if (!cal) {
7754cd1f 1358 return -1;
6d50f529 1359 }
7754cd1f 1360
95867fd1 1361 return cal->IsMasked(det,col,row);
2745a409 1362
7754cd1f 1363}
1364
1365//_____________________________________________________________________________
1366Bool_t AliTRDcalibDB::IsPadBridgedLeft(Int_t det, Int_t col, Int_t row)
1367{
1368 //
1369 // Returns status, see name of functions for details ;-)
1370 //
1371
95867fd1 1372 const AliTRDCalPadStatus *cal = dynamic_cast<const AliTRDCalPadStatus *>
1373 (GetCachedCDBObject(kIDPadStatus));
6d50f529 1374 if (!cal) {
7754cd1f 1375 return -1;
6d50f529 1376 }
7754cd1f 1377
95867fd1 1378 return cal->IsBridgedLeft(det,col,row);
2745a409 1379
7754cd1f 1380}
1381
1382//_____________________________________________________________________________
1383Bool_t AliTRDcalibDB::IsPadBridgedRight(Int_t det, Int_t col, Int_t row)
1384{
1385 //
1386 // Returns status, see name of functions for details ;-)
1387 //
1388
95867fd1 1389 const AliTRDCalPadStatus * cal = dynamic_cast<const AliTRDCalPadStatus *>
1390 (GetCachedCDBObject(kIDPadStatus));
6d50f529 1391 if (!cal) {
7754cd1f 1392 return -1;
6d50f529 1393 }
7754cd1f 1394
95867fd1 1395 return cal->IsBridgedRight(det,col,row);
2745a409 1396
7754cd1f 1397}
1398
fa7427d0 1399//_____________________________________________________________________________
1400Bool_t AliTRDcalibDB::IsPadNotConnected(Int_t det, Int_t col, Int_t row)
1401{
1402 //
1403 // Returns status, see name of functions for details ;-)
1404 //
f1dcad37 1405
fa7427d0 1406 const AliTRDCalPadStatus * cal = dynamic_cast<const AliTRDCalPadStatus *>
1407 (GetCachedCDBObject(kIDPadStatus));
1408 if (!cal) {
1409 return -1;
1410 }
1411
1412 return cal->IsNotConnected(det,col,row);
1413
1414}
1415
7754cd1f 1416//_____________________________________________________________________________
eda91732 1417Bool_t AliTRDcalibDB::IsChamberGood(Int_t det)
1418{
1419 //
1420 // Returns status, see name of functions for details ;-)
1421 //
1422
1423 const AliTRDCalChamberStatus * cal = dynamic_cast<const AliTRDCalChamberStatus *>
1424 (GetCachedCDBObject(kIDChamberStatus));
1425 if (!cal) {
1426 return -1;
1427 }
1428
1429 return cal->IsGood(det);
1430
1431}
1432
1433//_____________________________________________________________________________
1434Bool_t AliTRDcalibDB::IsChamberNoData(Int_t det)
7754cd1f 1435{
1436 //
1437 // Returns status, see name of functions for details ;-)
1438 //
1439
95867fd1 1440 const AliTRDCalChamberStatus * cal = dynamic_cast<const AliTRDCalChamberStatus *>
1441 (GetCachedCDBObject(kIDChamberStatus));
6d50f529 1442 if (!cal) {
7754cd1f 1443 return -1;
6d50f529 1444 }
7754cd1f 1445
eda91732 1446 return cal->IsNoData(det);
2745a409 1447
7754cd1f 1448}
1449
1450//_____________________________________________________________________________
eda91732 1451Bool_t AliTRDcalibDB::IsHalfChamberNoData(Int_t det, Int_t side)
7754cd1f 1452{
1453 //
1454 // Returns status, see name of functions for details ;-)
1455 //
1456
95867fd1 1457 const AliTRDCalChamberStatus * cal = dynamic_cast<const AliTRDCalChamberStatus *>
1458 (GetCachedCDBObject(kIDChamberStatus));
6d50f529 1459 if (!cal) {
7754cd1f 1460 return -1;
6d50f529 1461 }
7754cd1f 1462
eda91732 1463 return side > 0 ? cal->IsNoDataSideB(det) : cal->IsNoDataSideA(det);
2745a409 1464
7754cd1f 1465}
ea3eaa08 1466
7a6352a6 1467//_____________________________________________________________________________
eda91732 1468Bool_t AliTRDcalibDB::IsChamberBadCalibrated(Int_t det)
ea3eaa08 1469{
7a6352a6 1470 //
1471 // Returns status, see name of functions for details ;-)
1472 //
7754cd1f 1473
7a6352a6 1474 const AliTRDCalChamberStatus * cal = dynamic_cast<const AliTRDCalChamberStatus *>
1475 (GetCachedCDBObject(kIDChamberStatus));
1476 if (!cal) {
1477 return -1;
1478 }
1479
eda91732 1480 return cal->IsBadCalibrated(det);
7a6352a6 1481
1482}
ea3eaa08 1483
cc7cef99 1484//_____________________________________________________________________________
0d83b3a5 1485const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDpidUtil::ETRDPIDMethod method)
cc7cef99 1486{
1487 //
1488 // Returns the object storing the distributions for PID with likelihood
1489 //
9a96f175 1490
1491 switch(method) {
0d83b3a5 1492 case AliTRDpidUtil::kLQ:
4ba1d6ae 1493 return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDLQ));
0d83b3a5 1494 case AliTRDpidUtil::kNN:
4ba1d6ae 1495 return dynamic_cast<const AliTRDCalPID *>(GetCachedCDBObject(kIDPIDNN));
543691ea 1496 case AliTRDpidUtil::kESD:
1497 return 0x0; // To avoid compiler warnings
9a96f175 1498 }
1499
10f75631 1500 return 0x0;
2745a409 1501
cc7cef99 1502}
1503
9dcc64cc 1504//_____________________________________________________________________________
ea3eaa08 1505AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMethod method)
1506{
1507 //
1508 // Returns the PID response object for 1D-LQ
1509 //
1510
1511 if (!fPIDResponse) {
1512
9dcc64cc 1513 fPIDResponse = new AliTRDPIDResponse;
ea3eaa08 1514
9dcc64cc 1515 // Load Reference Histos from OCDB
1516 fPIDResponse->SetPIDmethod(method);
51a0ce25 1517 const TObjArray *references = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDPIDLQ1D));
ea3eaa08 1518
51a0ce25 1519 TIter refs(references);
1520 TObject *obj = NULL;
1521 AliTRDPIDReference *ref = NULL;
1522 Bool_t hasReference = kFALSE;
ea3eaa08 1523 while ((obj = refs())){
1524 if ((ref = dynamic_cast<AliTRDPIDReference *>(obj))){
1525 fPIDResponse->Load(ref);
1526 hasReference = kTRUE;
1527 break;
1528 }
51a0ce25 1529 }
ea3eaa08 1530
1531 if (!hasReference) {
1532 AliError("Reference histograms not found in the OCDB");
1533 }
1534
9dcc64cc 1535 }
ea3eaa08 1536
9dcc64cc 1537 return fPIDResponse;
ea3eaa08 1538
9dcc64cc 1539}
1540
1541//_____________________________________________________________________________
1542const AliTRDCalTrkAttach* AliTRDcalibDB::GetAttachObject()
1543{
1544 //
1545 // Returns the object storing likelihood distributions for cluster to track attachment
1546 //
ea3eaa08 1547
9dcc64cc 1548 return dynamic_cast<const AliTRDCalTrkAttach*>(GetCachedCDBObject(kIDAttach));
9dcc64cc 1549
ea3eaa08 1550}
9dcc64cc 1551
7754cd1f 1552//_____________________________________________________________________________
d4c6453d 1553const AliTRDCalMonitoring *AliTRDcalibDB::GetMonitoringObject()
7754cd1f 1554{
1555 //
1556 // Returns the object storing the monitoring data
1557 //
1558
d4c6453d 1559 return dynamic_cast<const AliTRDCalMonitoring *>
1560 (GetCachedCDBObject(kIDMonitoringData));
1b95a37b 1561
7754cd1f 1562}
1563
6a739e92 1564//_____________________________________________________________________________
1565void AliTRDcalibDB::SamplePRF()
1566{
1567 //
2745a409 1568 // Samples the pad response function (should maybe go somewhere else ...)
6a739e92 1569 //
1570
1571 const Int_t kPRFbin = 61;
1572
053767a4 1573 Float_t prf[kNlayer][kPRFbin] = {
6d50f529 1574 {2.9037e-02, 3.3608e-02, 3.9020e-02, 4.5292e-02,
6a739e92 1575 5.2694e-02, 6.1362e-02, 7.1461e-02, 8.3362e-02,
1576 9.7063e-02, 1.1307e-01, 1.3140e-01, 1.5235e-01,
1577 1.7623e-01, 2.0290e-01, 2.3294e-01, 2.6586e-01,
1578 3.0177e-01, 3.4028e-01, 3.8077e-01, 4.2267e-01,
1579 4.6493e-01, 5.0657e-01, 5.4655e-01, 5.8397e-01,
1580 6.1767e-01, 6.4744e-01, 6.7212e-01, 6.9188e-01,
1581 7.0627e-01, 7.1499e-01, 7.1851e-01, 7.1499e-01,
1582 7.0627e-01, 6.9188e-01, 6.7212e-01, 6.4744e-01,
1583 6.1767e-01, 5.8397e-01, 5.4655e-01, 5.0657e-01,
1584 4.6493e-01, 4.2267e-01, 3.8077e-01, 3.4028e-01,
1585 3.0177e-01, 2.6586e-01, 2.3294e-01, 2.0290e-01,
1586 1.7623e-01, 1.5235e-01, 1.3140e-01, 1.1307e-01,
1587 9.7063e-02, 8.3362e-02, 7.1461e-02, 6.1362e-02,
1588 5.2694e-02, 4.5292e-02, 3.9020e-02, 3.3608e-02,
1589 2.9037e-02},
1590 {2.5478e-02, 2.9695e-02, 3.4655e-02, 4.0454e-02,
1591 4.7342e-02, 5.5487e-02, 6.5038e-02, 7.6378e-02,
1592 8.9696e-02, 1.0516e-01, 1.2327e-01, 1.4415e-01,
1593 1.6794e-01, 1.9516e-01, 2.2573e-01, 2.5959e-01,
1594 2.9694e-01, 3.3719e-01, 3.7978e-01, 4.2407e-01,
1595 4.6889e-01, 5.1322e-01, 5.5569e-01, 5.9535e-01,
1596 6.3141e-01, 6.6259e-01, 6.8882e-01, 7.0983e-01,
1597 7.2471e-01, 7.3398e-01, 7.3761e-01, 7.3398e-01,
1598 7.2471e-01, 7.0983e-01, 6.8882e-01, 6.6259e-01,
1599 6.3141e-01, 5.9535e-01, 5.5569e-01, 5.1322e-01,
1600 4.6889e-01, 4.2407e-01, 3.7978e-01, 3.3719e-01,
1601 2.9694e-01, 2.5959e-01, 2.2573e-01, 1.9516e-01,
1602 1.6794e-01, 1.4415e-01, 1.2327e-01, 1.0516e-01,
1603 8.9696e-02, 7.6378e-02, 6.5038e-02, 5.5487e-02,
1604 4.7342e-02, 4.0454e-02, 3.4655e-02, 2.9695e-02,
1605 2.5478e-02},
1606 {2.2363e-02, 2.6233e-02, 3.0782e-02, 3.6140e-02,
1607 4.2535e-02, 5.0157e-02, 5.9197e-02, 6.9900e-02,
1608 8.2707e-02, 9.7811e-02, 1.1548e-01, 1.3601e-01,
1609 1.5998e-01, 1.8739e-01, 2.1840e-01, 2.5318e-01,
1610 2.9182e-01, 3.3373e-01, 3.7837e-01, 4.2498e-01,
1611 4.7235e-01, 5.1918e-01, 5.6426e-01, 6.0621e-01,
1612 6.4399e-01, 6.7700e-01, 7.0472e-01, 7.2637e-01,
1613 7.4206e-01, 7.5179e-01, 7.5551e-01, 7.5179e-01,
1614 7.4206e-01, 7.2637e-01, 7.0472e-01, 6.7700e-01,
1615 6.4399e-01, 6.0621e-01, 5.6426e-01, 5.1918e-01,
1616 4.7235e-01, 4.2498e-01, 3.7837e-01, 3.3373e-01,
1617 2.9182e-01, 2.5318e-01, 2.1840e-01, 1.8739e-01,
1618 1.5998e-01, 1.3601e-01, 1.1548e-01, 9.7811e-02,
1619 8.2707e-02, 6.9900e-02, 5.9197e-02, 5.0157e-02,
1620 4.2535e-02, 3.6140e-02, 3.0782e-02, 2.6233e-02,
1621 2.2363e-02},
1622 {1.9635e-02, 2.3167e-02, 2.7343e-02, 3.2293e-02,
1623 3.8224e-02, 4.5335e-02, 5.3849e-02, 6.4039e-02,
1624 7.6210e-02, 9.0739e-02, 1.0805e-01, 1.2841e-01,
1625 1.5216e-01, 1.7960e-01, 2.1099e-01, 2.4671e-01,
1626 2.8647e-01, 3.2996e-01, 3.7660e-01, 4.2547e-01,
1627 4.7536e-01, 5.2473e-01, 5.7215e-01, 6.1632e-01,
1628 6.5616e-01, 6.9075e-01, 7.1939e-01, 7.4199e-01,
1629 7.5838e-01, 7.6848e-01, 7.7227e-01, 7.6848e-01,
1630 7.5838e-01, 7.4199e-01, 7.1939e-01, 6.9075e-01,
1631 6.5616e-01, 6.1632e-01, 5.7215e-01, 5.2473e-01,
1632 4.7536e-01, 4.2547e-01, 3.7660e-01, 3.2996e-01,
1633 2.8647e-01, 2.4671e-01, 2.1099e-01, 1.7960e-01,
1634 1.5216e-01, 1.2841e-01, 1.0805e-01, 9.0739e-02,
1635 7.6210e-02, 6.4039e-02, 5.3849e-02, 4.5335e-02,
1636 3.8224e-02, 3.2293e-02, 2.7343e-02, 2.3167e-02,
1637 1.9635e-02},
1638 {1.7224e-02, 2.0450e-02, 2.4286e-02, 2.8860e-02,
1639 3.4357e-02, 4.0979e-02, 4.8966e-02, 5.8612e-02,
1640 7.0253e-02, 8.4257e-02, 1.0102e-01, 1.2094e-01,
1641 1.4442e-01, 1.7196e-01, 2.0381e-01, 2.4013e-01,
1642 2.8093e-01, 3.2594e-01, 3.7450e-01, 4.2563e-01,
1643 4.7796e-01, 5.2991e-01, 5.7974e-01, 6.2599e-01,
1644 6.6750e-01, 7.0344e-01, 7.3329e-01, 7.5676e-01,
1645 7.7371e-01, 7.8410e-01, 7.8793e-01, 7.8410e-01,
1646 7.7371e-01, 7.5676e-01, 7.3329e-01, 7.0344e-01,
1647 6.6750e-01, 6.2599e-01, 5.7974e-01, 5.2991e-01,
1648 4.7796e-01, 4.2563e-01, 3.7450e-01, 3.2594e-01,
1649 2.8093e-01, 2.4013e-01, 2.0381e-01, 1.7196e-01,
1650 1.4442e-01, 1.2094e-01, 1.0102e-01, 8.4257e-02,
1651 7.0253e-02, 5.8612e-02, 4.8966e-02, 4.0979e-02,
1652 3.4357e-02, 2.8860e-02, 2.4286e-02, 2.0450e-02,
1653 1.7224e-02},
1654 {1.5096e-02, 1.8041e-02, 2.1566e-02, 2.5793e-02,
1655 3.0886e-02, 3.7044e-02, 4.4515e-02, 5.3604e-02,
1656 6.4668e-02, 7.8109e-02, 9.4364e-02, 1.1389e-01,
1657 1.3716e-01, 1.6461e-01, 1.9663e-01, 2.3350e-01,
1658 2.7527e-01, 3.2170e-01, 3.7214e-01, 4.2549e-01,
1659 4.8024e-01, 5.3460e-01, 5.8677e-01, 6.3512e-01,
1660 6.7838e-01, 7.1569e-01, 7.4655e-01, 7.7071e-01,
1661 7.8810e-01, 7.9871e-01, 8.0255e-01, 7.9871e-01,
1662 7.8810e-01, 7.7071e-01, 7.4655e-01, 7.1569e-01,
1663 6.7838e-01, 6.3512e-01, 5.8677e-01, 5.3460e-01,
1664 4.8024e-01, 4.2549e-01, 3.7214e-01, 3.2170e-01,
1665 2.7527e-01, 2.3350e-01, 1.9663e-01, 1.6461e-01,
1666 1.3716e-01, 1.1389e-01, 9.4364e-02, 7.8109e-02,
1667 6.4668e-02, 5.3604e-02, 4.4515e-02, 3.7044e-02,
1668 3.0886e-02, 2.5793e-02, 2.1566e-02, 1.8041e-02,
1669 1.5096e-02}};
1670
1671 // More sampling precision with linear interpolation
6d50f529 1672 fPRFlo = -1.5;
1673 fPRFhi = 1.5;
6a739e92 1674 Float_t pad[kPRFbin];
1675 Int_t sPRFbin = kPRFbin;
6d50f529 1676 Float_t sPRFwid = (fPRFhi - fPRFlo) / ((Float_t) sPRFbin);
6a739e92 1677 for (Int_t iPad = 0; iPad < sPRFbin; iPad++) {
6d50f529 1678 pad[iPad] = ((Float_t) iPad + 0.5) * sPRFwid + fPRFlo;
6a739e92 1679 }
6d50f529 1680 fPRFbin = 500;
1681 fPRFwid = (fPRFhi - fPRFlo) / ((Float_t) fPRFbin);
1682 fPRFpad = ((Int_t) (1.0 / fPRFwid));
6a739e92 1683
6d50f529 1684 if (fPRFsmp) delete [] fPRFsmp;
053767a4 1685 fPRFsmp = new Float_t[kNlayer*fPRFbin];
6a739e92 1686
1687 Int_t ipos1;
1688 Int_t ipos2;
1689 Float_t diff;
1690
053767a4 1691 for (Int_t iLayer = 0; iLayer < kNlayer; iLayer++) {
6a739e92 1692
6d50f529 1693 for (Int_t iBin = 0; iBin < fPRFbin; iBin++) {
6a739e92 1694
6d50f529 1695 Float_t bin = (((Float_t) iBin) + 0.5) * fPRFwid + fPRFlo;
6a739e92 1696 ipos1 = ipos2 = 0;
1697 diff = 0;
1698 do {
1699 diff = bin - pad[ipos2++];
1700 } while ((diff > 0) && (ipos2 < kPRFbin));
1701 if (ipos2 == kPRFbin) {
053767a4 1702 fPRFsmp[iLayer*fPRFbin+iBin] = prf[iLayer][ipos2-1];
6a739e92 1703 }
1704 else if (ipos2 == 1) {
053767a4 1705 fPRFsmp[iLayer*fPRFbin+iBin] = prf[iLayer][ipos2-1];
6a739e92 1706 }
1707 else {
1708 ipos2--;
1709 if (ipos2 >= kPRFbin) ipos2 = kPRFbin - 1;
1710 ipos1 = ipos2 - 1;
053767a4 1711 fPRFsmp[iLayer*fPRFbin+iBin] = prf[iLayer][ipos2]
1712 + diff * (prf[iLayer][ipos2] - prf[iLayer][ipos1])
1713 / sPRFwid;
6a739e92 1714 }
1715
1716 }
1717 }
1718
1719}
1720
1721//_____________________________________________________________________________
1722Int_t AliTRDcalibDB::PadResponse(Double_t signal, Double_t dist
eb52b657 1723 , Int_t layer, Double_t *pad) const
6a739e92 1724{
1725 //
1726 // Applies the pad response
56178ff4 1727 // So far this is the fixed parametrization and should be replaced by
1728 // something dependent on calibration values
6a739e92 1729 //
1730
eb52b657 1731 Int_t iBin = ((Int_t) ((-dist - fPRFlo) / fPRFwid));
053767a4 1732 Int_t iOff = layer * fPRFbin;
6a739e92 1733
6d50f529 1734 Int_t iBin0 = iBin - fPRFpad + iOff;
6a739e92 1735 Int_t iBin1 = iBin + iOff;
6d50f529 1736 Int_t iBin2 = iBin + fPRFpad + iOff;
6a739e92 1737
1738 pad[0] = 0.0;
1739 pad[1] = 0.0;
1740 pad[2] = 0.0;
053767a4 1741 if ((iBin1 >= 0) && (iBin1 < (fPRFbin*kNlayer))) {
6a739e92 1742
1743 if (iBin0 >= 0) {
6d50f529 1744 pad[0] = signal * fPRFsmp[iBin0];
6a739e92 1745 }
6d50f529 1746 pad[1] = signal * fPRFsmp[iBin1];
053767a4 1747 if (iBin2 < (fPRFbin*kNlayer)) {
6d50f529 1748 pad[2] = signal * fPRFsmp[iBin2];
6a739e92 1749 }
1750
1751 return 1;
1752
1753 }
1754 else {
1755
1756 return 0;
1757
1758 }
ab0a4106 1759
3551db50 1760}
7a6352a6 1761
1762
a9e22145 1763AliTRDtrapConfig* AliTRDcalibDB::GetTrapConfig()
1764{
1765 // return an existing TRAPconfig or load it from the OCDB
1766 // in case of failure, a default TRAPconfig is created
1767
1768 if (fTrapConfig)
1769 return fTrapConfig;
1770 else {
1b1afdd3 1771 if ((fTrapConfigName.Length() <= 0) || (fTrapConfigVersion.Length() <= 0)) {
1772 // query the configuration to be used
1773 TString configName;
1774 this->GetGlobalConfiguration(configName);
1775 TString configVersion;
1776 this->GetGlobalConfigurationVersion(configVersion);
1777 }
a9e22145 1778
1b1afdd3 1779 // try to load the requested configuration
1780 this->LoadTrapConfig(fTrapConfigName, fTrapConfigVersion);
a9e22145 1781
1782 // if we still don't have a valid TRAPconfig, create a default one
1783 if (!fTrapConfig) {
1784 AliWarning("Falling back to default configuration");
1785 fTrapConfig = new AliTRDtrapConfig("default", "default TRAP configuration");
1786 AliTRDtrapConfigHandler cfgHandler(fTrapConfig);
1787 cfgHandler.Init();
1788 cfgHandler.LoadConfig();
1789 }
1790
1791 AliInfo(Form("using TRAPconfig \"%s\"", fTrapConfig->GetTitle()));
1792
1793 return fTrapConfig;
1794 }
1795}
1796
1797
1798AliTRDtrapConfig* AliTRDcalibDB::LoadTrapConfig(const TString &name, const TString &version)
1799{
1800 // try to load the specified configuration from the OCDB
1801 // if it fails, or it does not exist, return null
1802
1803 AliInfo(Form("looking for TRAPconfig \"%s.%s\"", name.Data(), version.Data()));
1804
1805 const AliTRDCalTrapConfig *caltrap = dynamic_cast<const AliTRDCalTrapConfig*> (GetCachedCDBObject(kIDTrapConfig));
1806
1807 if (caltrap) {
1808 TString configName(name);
1809 configName.Append(".");
1810 configName.Append(version);
1811 fTrapConfig = caltrap->Get(configName);
1812 }
1813 else {
1814 fTrapConfig = 0x0;
1815 AliError("No TRAPconfig entry found!");
1816 }
1817
1818 return fTrapConfig;
1819}