Adding a trigger selection combo-box to be used by the online display/reco. The combo...
[u/mrichter/AliRoot.git] / PHOS / AliPHOSCalibData.cxx
CommitLineData
52783dbc 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 for PHOS calibration //
21// //
22///////////////////////////////////////////////////////////////////////////////
23
fc6706cb 24#include "TRandom.h"
76d78b7a 25#include "AliLog.h"
52783dbc 26#include "AliPHOSCalibData.h"
fc6706cb 27#include "AliCDBManager.h"
28#include "AliCDBStorage.h"
29#include "AliCDBId.h"
30#include "AliCDBEntry.h"
1fdb1c5e 31#include "AliPHOSEmcCalibData.h"
32#include "AliPHOSCpvCalibData.h"
81701069 33#include "AliPHOSEmcBadChannelsMap.h"
1fdb1c5e 34#include "AliCDBMetaData.h"
52783dbc 35
36ClassImp(AliPHOSCalibData)
37
38//________________________________________________________________
fc6706cb 39 AliPHOSCalibData::AliPHOSCalibData():
3f7dbdb7 40 TNamed(),
41 fCalibDataEmc(0x0),
42 fCalibDataCpv(0x0),
81701069 43 fEmcBadChannelsMap(0x0),
3f7dbdb7 44 fEmcDataPath("PHOS/Calib/EmcGainPedestals"),
81701069 45 fCpvDataPath("PHOS/Calib/CpvGainPedestals"),
46 fEmcBadChannelsMapPath("PHOS/Calib/EmcBadChannels")
52783dbc 47{
76d78b7a 48 // Default constructor.
49 // Open CDB entry, get EMC and CPV calibration data and bad channel map.
50 // If EMC or CPV calibration data does not exist, stop the run
3f7dbdb7 51
81701069 52 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
abcff5fc 53 if(entryEmc)
54 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
76d78b7a 55
56 if(!fCalibDataEmc)
f36dfa9a 57 AliWarning("Calibration parameters for PHOS EMC not found. Create a new set.\n");
76d78b7a 58
abcff5fc 59 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
60 if(entryCpv)
61 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
62
76d78b7a 63 if(!fCalibDataCpv)
f36dfa9a 64 AliWarning("Calibration parameters for PHOS CPV not found. Create a new set.\n");
76d78b7a 65
81701069 66 AliCDBEntry* entryEmcBadMap = AliCDBManager::Instance()->Get(fEmcBadChannelsMapPath.Data());
67 if(entryEmcBadMap)
68 fEmcBadChannelsMap = (AliPHOSEmcBadChannelsMap*)entryEmcBadMap->GetObject();
69
52783dbc 70}
71
72//________________________________________________________________
fc6706cb 73AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
74 TNamed("phosCalib","PHOS Calibration Data Manager"),
81701069 75 fCalibDataEmc(0x0), fCalibDataCpv(0x0), fEmcBadChannelsMap(0x0),
3f7dbdb7 76 fEmcDataPath("PHOS/Calib/EmcGainPedestals"),
81701069 77 fCpvDataPath("PHOS/Calib/CpvGainPedestals"),
78 fEmcBadChannelsMapPath("PHOS/Calib/EmcBadChannels")
52783dbc 79{
80 // Constructor
76d78b7a 81 // Open CDB entry, get EMC and CPV calibration data and bad channel map.
82 // If EMC or CPV calibration data does not exist, stop the run
83
fc6706cb 84 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
85 if(entryEmc)
86 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
76d78b7a 87
88 if(!fCalibDataEmc)
89 AliFatal("Calibration parameters for PHOS EMC not found. Stop reconstruction!\n");
fc6706cb 90
91 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
92 if(entryCpv)
93 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
52783dbc 94
76d78b7a 95 if(!fCalibDataCpv)
96 AliFatal("Calibration parameters for PHOS CPV not found. Stop reconstruction!\n");
97
81701069 98 AliCDBEntry* entryEmcBadMap = AliCDBManager::Instance()->
99 Get(fEmcBadChannelsMapPath.Data(),runNumber);
100 if(entryEmcBadMap)
101 fEmcBadChannelsMap = (AliPHOSEmcBadChannelsMap*)entryEmcBadMap->GetObject();
102
52783dbc 103}
104
105//________________________________________________________________
1fdb1c5e 106AliPHOSCalibData::AliPHOSCalibData(AliPHOSCalibData & phosCDB) :
3f7dbdb7 107 TNamed(phosCDB),
108 fCalibDataEmc(phosCDB.fCalibDataEmc),
109 fCalibDataCpv(phosCDB.fCalibDataCpv),
81701069 110 fEmcBadChannelsMap(phosCDB.fEmcBadChannelsMap),
3f7dbdb7 111 fEmcDataPath(phosCDB.fEmcDataPath),
81701069 112 fCpvDataPath(phosCDB.fCpvDataPath),
113 fEmcBadChannelsMapPath(phosCDB.fEmcBadChannelsMapPath)
1fdb1c5e 114{
115 // Copy constructor
1fdb1c5e 116}
117//________________________________________________________________
52783dbc 118AliPHOSCalibData::~AliPHOSCalibData()
119{
120 // Destructor
fc6706cb 121
52783dbc 122}
123
3f7dbdb7 124AliPHOSCalibData & AliPHOSCalibData::operator = (const AliPHOSCalibData & rhs)
125{
126 //Copy-assignment. Does not delete anything (see destructor)
127 //compiler generated is ok, but ... because -Weffc++ and pointer
128 //members we have to define it explicitly.
129 TNamed::operator=(rhs);
130 fCalibDataEmc = rhs.fCalibDataEmc;
131 fCalibDataCpv = rhs.fCalibDataCpv;
81701069 132 fEmcBadChannelsMap = rhs.fEmcBadChannelsMap;
3f7dbdb7 133 fEmcDataPath = rhs.fEmcDataPath;
134 fCpvDataPath = rhs.fCpvDataPath;
81701069 135 fEmcBadChannelsMapPath = rhs.fEmcBadChannelsMapPath;
3f7dbdb7 136
137 return *this;
138}
139
52783dbc 140//________________________________________________________________
141void AliPHOSCalibData::Reset()
142{
b8628054 143 // Set all pedestals to 0 and all ADC channels to 1,
144 // and all channels are good (alive)
fc6706cb 145
b8628054 146 fCalibDataEmc ->Reset();
147 fCalibDataCpv ->Reset();
148 fEmcBadChannelsMap->Reset();
52783dbc 149}
150
151//________________________________________________________________
152void AliPHOSCalibData::Print(Option_t *option) const
153{
1fdb1c5e 154 // Print EMC and CPV calibration containers
155 // Input: option="ped" to print pedestals
156 // option="gain" to print calibration coefficients
df0d5123 157 if (fCalibDataEmc) fCalibDataEmc->Print(option);
158 if (fCalibDataCpv) fCalibDataCpv->Print(option);
fc6706cb 159}
160
161//________________________________________________________________
162void AliPHOSCalibData::CreateNew()
163{
1fdb1c5e 164 // Create new EMC and CPV calibration containers with ideal coefficients
165
fc6706cb 166 if(fCalibDataEmc) delete fCalibDataEmc;
167 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
168
169 if(fCalibDataCpv) delete fCalibDataCpv;
170 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
171
81701069 172 if(fEmcBadChannelsMap) delete fEmcBadChannelsMap;
173 fEmcBadChannelsMap = new AliPHOSEmcBadChannelsMap();
174
fc6706cb 175}
176
177//________________________________________________________________
178Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
179{
1fdb1c5e 180 // Write EMC calibration container to CDB
fc6706cb 181
182 if(!fCalibDataEmc) return kFALSE;
183
18669c5c 184 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS/*");
df0d5123 185 if(!storage)
186 storage = AliCDBManager::Instance()->GetDefaultStorage();
187
fc6706cb 188 if(storage) {
189 AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
190 storage->Put(fCalibDataEmc,id, md);
191 return kTRUE;
52783dbc 192 }
fc6706cb 193 else
194 return kFALSE;
52783dbc 195
fc6706cb 196}
197
198//________________________________________________________________
199Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
200{
1fdb1c5e 201 // Write CPV calibration container to CDB
fc6706cb 202
203 if(!fCalibDataCpv) return kFALSE;
204
18669c5c 205 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS/*");
df0d5123 206 if(!storage)
207 storage = AliCDBManager::Instance()->GetDefaultStorage();
208
fc6706cb 209 if(storage) {
210 AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
211 storage->Put(fCalibDataCpv,id, md);
212 return kTRUE;
52783dbc 213 }
fc6706cb 214 else
215 return kFALSE;
216
52783dbc 217}
072de3a8 218
81701069 219
220//________________________________________________________________
221Bool_t AliPHOSCalibData::WriteEmcBadChannelsMap(Int_t firstRun,Int_t lastRun,AliCDBMetaData *md)
222{
223 //Write EMC bad channels map into CDB.
224
225 if(!fEmcBadChannelsMap) return kFALSE;
226
227 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS/*");
228 if(!storage)
229 storage = AliCDBManager::Instance()->GetDefaultStorage();
230
231 if(storage) {
232 AliCDBId id(fEmcBadChannelsMapPath.Data(),firstRun,lastRun);
233 storage->Put(fEmcBadChannelsMap,id, md);
234 return kTRUE;
235 }
236 else
237 return kFALSE;
238}
239
fc6706cb 240//________________________________________________________________
072de3a8 241Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
242{
1fdb1c5e 243 // Return EMC calibration coefficient
244 // for channel defined by (module,column,row)
245 // module, column,raw should follow the internal PHOS convention:
246 // module 1:5, column 1:56, row 1:64
247 // if CBD instance exists, the value is taken from CDB.
248 // Otherwise it is an ideal one
072de3a8 249
fc6706cb 250 if(fCalibDataEmc)
251 return fCalibDataEmc->GetADCchannelEmc(module,column,row);
252 else
5f8bb860 253 return 1.0; // default width of one EMC ADC channel in GeV
072de3a8 254}
255
1fdb1c5e 256//________________________________________________________________
e3310625 257void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
258{
259 // Set EMC calibration coefficient for (module,column,row)
260
261 if(!fCalibDataEmc)
262 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
263
264 fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
265}
266
267//________________________________________________________________
072de3a8 268Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
269{
1fdb1c5e 270 // Return EMC pedestal for channel defined by (module,column,row)
271 // module, column,raw should follow the internal PHOS convention:
272 // module 1:5, column 1:56, row 1:64
273 // if CBD instance exists, the value is taken from CDB.
274 // Otherwise it is an ideal one
275
fc6706cb 276 if(fCalibDataEmc)
277 return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
278 else
5f8bb860 279 return 0.0; // default EMC ADC pedestal
072de3a8 280}
281
1fdb1c5e 282//________________________________________________________________
e3310625 283void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
072de3a8 284{
e3310625 285 // Set EMC pedestal for (module,column,row)
1fdb1c5e 286
fc6706cb 287 if(!fCalibDataEmc)
288 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
289
e3310625 290 fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
072de3a8 291}
292
1fdb1c5e 293//________________________________________________________________
e3310625 294Float_t AliPHOSCalibData::GetHighLowRatioEmc(Int_t module, Int_t column, Int_t row) const
072de3a8 295{
e3310625 296 // Return EMC calibration coefficient
297 // for channel defined by (module,column,row)
298 // module, column,raw should follow the internal PHOS convention:
299 // module 1:5, column 1:56, row 1:64
300 // if CBD instance exists, the value is taken from CDB.
301 // Otherwise it is an ideal one
302
303 if(fCalibDataEmc)
304 return fCalibDataEmc->GetHighLowRatioEmc(module,column,row);
305 else
306 return 1.0; // default width of one EMC ADC channel in GeV
307}
308
309//________________________________________________________________
310void AliPHOSCalibData::SetHighLowRatioEmc(Int_t module, Int_t column, Int_t row, Float_t value)
311{
312 // Set EMC calibration coefficient for (module,column,row)
313
fc6706cb 314 if(!fCalibDataEmc)
315 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
e3310625 316
317 fCalibDataEmc->SetHighLowRatioEmc(module,column,row,value);
fc6706cb 318}
e3310625 319
320//________________________________________________________________
321Float_t AliPHOSCalibData::GetTimeShiftEmc(Int_t module, Int_t column, Int_t row) const
322{
323 // Return EMC calibration coefficient
324 // for channel defined by (module,column,row)
325 // module, column,raw should follow the internal PHOS convention:
326 // module 1:5, column 1:56, row 1:64
327 // if CBD instance exists, the value is taken from CDB.
328 // Otherwise it is an ideal one
329
330 if(fCalibDataEmc)
331 return fCalibDataEmc->GetTimeShiftEmc(module,column,row);
332 else
333 return 1.0; // default width of one EMC ADC channel in GeV
334}
335
336//________________________________________________________________
337void AliPHOSCalibData::SetTimeShiftEmc(Int_t module, Int_t column, Int_t row, Float_t value)
338{
339 // Set EMC calibration coefficient for (module,column,row)
340
341 if(!fCalibDataEmc)
342 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
343
344 fCalibDataEmc->SetTimeShiftEmc(module,column,row,value);
345}
bafc1087 346//________________________________________________________________
347Int_t AliPHOSCalibData::GetAltroOffsetEmc(Int_t module, Int_t column, Int_t row) const
348{
349 // Return ALTRO pedestal coefficient
350 // for channel defined by (module,column,row)
351 // module, column,raw should follow the internal PHOS convention:
352 // module 1:5, column 1:56, row 1:64
353 // if CBD instance exists, the value is taken from CDB.
354 // Otherwise it is an ideal one
355
356 if(fCalibDataEmc)
357 return fCalibDataEmc->GetAltroOffsetEmc(module,column,row);
358 else
359 return 0; // default width of one EMC ADC channel in GeV
360}
361
362//________________________________________________________________
363void AliPHOSCalibData::SetAltroOffsetEmc(Int_t module, Int_t column, Int_t row, Int_t value)
364{
365 // Set altro offset for (module,column,row)
366
367 if(!fCalibDataEmc)
368 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
369
370 fCalibDataEmc->SetAltroOffsetEmc(module,column,row,value);
371}
372
e3310625 373
fc6706cb 374//________________________________________________________________
375Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
376{
1fdb1c5e 377 // Return CPV calibration coefficient
378 // for channel defined by (module,column,row)
379 // module, column,raw should follow the internal CPV convention:
380 // module 1:5, column 1:56, row 1:128
381 // if CBD instance exists, the value is taken from CDB.
382 // Otherwise it is an ideal one
fc6706cb 383
384 if(fCalibDataCpv)
385 return fCalibDataCpv->GetADCchannelCpv(module,column,row);
386 else
5f8bb860 387 return 0.0012; // default width of one ADC channel in CPV arbitrary units
fc6706cb 388}
389
1fdb1c5e 390//________________________________________________________________
fc6706cb 391Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
392{
1fdb1c5e 393 // Return CPV pedestal
394 // for channel defined by (module,column,row)
395 // module, column,raw should follow the internal CPV convention:
396 // module 1:5, column 1:56, row 1:128
397 // if CBD instance exists, the value is taken from CDB.
398 // Otherwise it is an ideal one
399
fc6706cb 400 if(fCalibDataCpv)
401 return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
402 else
403 return 0.012; // default CPV ADC pedestal
404}
405
1fdb1c5e 406//________________________________________________________________
fc6706cb 407void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
408{
1fdb1c5e 409 // Set CPV calibration coefficient for (module,column,row)
410
fc6706cb 411 if(!fCalibDataCpv)
412 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
413
414 fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
415}
416
1fdb1c5e 417//________________________________________________________________
fc6706cb 418void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
419{
1fdb1c5e 420 // Set CPV pedestal for (module,column,row)
421
fc6706cb 422 if(!fCalibDataCpv)
423 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
424
425 fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
426}
427
428//________________________________________________________________
e2e977c3 429void AliPHOSCalibData::RandomEmc(Float_t ccMin, Float_t ccMax)
fc6706cb 430{
1fdb1c5e 431 // Create decalibrated EMC with calibration coefficients and pedestals
432 // randomly distributed within hard-coded limits
e2e977c3 433 // Default spread of calibration parameters is Cmax/Cmin = 4, (Cmax-Cmin)/2 = 1
fc6706cb 434
435 if(fCalibDataEmc) delete fCalibDataEmc;
436 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
437
438 TRandom rn;
439 rn.SetSeed(0); //the seed is set to the current machine clock
440
1fdb1c5e 441 Float_t adcChannelEmc,adcPedestalEmc;
fc6706cb 442
443 for(Int_t module=1; module<6; module++) {
444 for(Int_t column=1; column<57; column++) {
445 for(Int_t row=1; row<65; row++) {
e2e977c3 446 adcChannelEmc =rn.Uniform(ccMin,ccMax);
5f8bb860 447 adcPedestalEmc=rn.Uniform(0.0,0.0); // 0 spread of pedestals
1fdb1c5e 448 fCalibDataEmc->SetADCchannelEmc(module,column,row,adcChannelEmc);
449 fCalibDataEmc->SetADCpedestalEmc(module,column,row,adcPedestalEmc);
fc6706cb 450 }
451 }
452 }
453
454}
455
456//________________________________________________________________
e2e977c3 457void AliPHOSCalibData::RandomCpv(Float_t ccMin, Float_t ccMax)
fc6706cb 458{
1fdb1c5e 459 // Create decalibrated CPV with calibration coefficients and pedestals
460 // randomly distributed within hard-coded limits
e2e977c3 461 // Default spread of calibration parameters is 0.0012 +- 25%
fc6706cb 462
463 if(fCalibDataCpv) delete fCalibDataCpv;
464 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
465
466 TRandom rn;
467 rn.SetSeed(0); //the seed is set to the current machine clock
468
1fdb1c5e 469 Float_t adcChannelCpv,adcPedestalCpv;
fc6706cb 470
471 for(Int_t module=1; module<6; module++) {
5f8bb860 472 for(Int_t column=1; column<57; column++) {
fc6706cb 473 for(Int_t row=1; row<129; row++) {
e2e977c3 474 adcChannelCpv =rn.Uniform(ccMin,ccMax);
1fdb1c5e 475 adcPedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
476 fCalibDataCpv->SetADCchannelCpv(module,column,row,adcChannelCpv);
477 fCalibDataCpv->SetADCpedestalCpv(module,column,row,adcPedestalCpv);
fc6706cb 478 }
479 }
480 }
072de3a8 481}
81701069 482
483//________________________________________________________________
484Bool_t AliPHOSCalibData::IsBadChannelEmc(Int_t module, Int_t col, Int_t row) const
485{
486 //If no bad channels map found, channel considered good
487
488 if(fEmcBadChannelsMap)
489 return fEmcBadChannelsMap->IsBadChannel(module,col,row);
490 else
491 return kFALSE;
492}
493
494//________________________________________________________________
495Int_t AliPHOSCalibData::GetNumOfEmcBadChannels() const
496{
497 if(fEmcBadChannelsMap)
498 return fEmcBadChannelsMap->GetNumOfBadChannels();
499 else
500 return 0;
501}
d940eb29 502
503//________________________________________________________________
504void AliPHOSCalibData::EmcBadChannelIds(Int_t *badIds)
505{
506 //Fill array badIds by the Ids of EMC bad channels.
507 //Array badIds of length GetNumOfBadChannels() should be prepared in advance.
508
509 if(fEmcBadChannelsMap)
510 fEmcBadChannelsMap->BadChannelIds(badIds);
511}