1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // class for PHOS calibration //
22 ///////////////////////////////////////////////////////////////////////////////
25 #include "AliPHOSCalibData.h"
26 #include "AliCDBManager.h"
27 #include "AliCDBStorage.h"
29 #include "AliCDBEntry.h"
30 #include "AliPHOSEmcCalibData.h"
31 #include "AliPHOSCpvCalibData.h"
32 #include "AliCDBMetaData.h"
34 ClassImp(AliPHOSCalibData)
36 //________________________________________________________________
37 AliPHOSCalibData::AliPHOSCalibData():
38 TNamed(), fCalibDataEmc(0x0), fCalibDataCpv(0x0)
40 // Default constructor
42 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
43 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
45 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
47 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
49 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
51 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
55 //________________________________________________________________
56 AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
57 TNamed("phosCalib","PHOS Calibration Data Manager"),
58 fCalibDataEmc(0x0), fCalibDataCpv(0x0)
62 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
63 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
65 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
67 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
69 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
71 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
75 //________________________________________________________________
76 AliPHOSCalibData::AliPHOSCalibData(AliPHOSCalibData & phosCDB) :
81 fCalibDataEmc = phosCDB.fCalibDataEmc;
82 fCalibDataCpv = phosCDB.fCalibDataCpv;
84 fEmcDataPath = phosCDB.fEmcDataPath;
85 fCpvDataPath = phosCDB.fCpvDataPath;
87 //________________________________________________________________
88 AliPHOSCalibData::~AliPHOSCalibData()
94 //________________________________________________________________
95 void AliPHOSCalibData::Reset()
97 // Set all pedestals to 0 and all ADC channels to 1
99 fCalibDataEmc->Reset();
100 fCalibDataCpv->Reset();
103 //________________________________________________________________
104 void AliPHOSCalibData::Print(Option_t *option) const
106 // Print EMC and CPV calibration containers
107 // Input: option="ped" to print pedestals
108 // option="gain" to print calibration coefficients
109 if (fCalibDataEmc) fCalibDataEmc->Print(option);
110 if (fCalibDataCpv) fCalibDataCpv->Print(option);
113 //________________________________________________________________
114 void AliPHOSCalibData::CreateNew()
116 // Create new EMC and CPV calibration containers with ideal coefficients
118 if(fCalibDataEmc) delete fCalibDataEmc;
119 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
121 if(fCalibDataCpv) delete fCalibDataCpv;
122 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
126 //________________________________________________________________
127 Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
129 // Write EMC calibration container to CDB
131 if(!fCalibDataEmc) return kFALSE;
133 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
135 storage = AliCDBManager::Instance()->GetDefaultStorage();
138 AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
139 storage->Put(fCalibDataEmc,id, md);
147 //________________________________________________________________
148 Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
150 // Write CPV calibration container to CDB
152 if(!fCalibDataCpv) return kFALSE;
154 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
156 storage = AliCDBManager::Instance()->GetDefaultStorage();
159 AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
160 storage->Put(fCalibDataCpv,id, md);
168 //________________________________________________________________
169 Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
171 // Return EMC calibration coefficient
172 // for channel defined by (module,column,row)
173 // module, column,raw should follow the internal PHOS convention:
174 // module 1:5, column 1:56, row 1:64
175 // if CBD instance exists, the value is taken from CDB.
176 // Otherwise it is an ideal one
179 return fCalibDataEmc->GetADCchannelEmc(module,column,row);
181 return 1.0; // default width of one EMC ADC channel in GeV
184 //________________________________________________________________
185 Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
187 // Return EMC pedestal for channel defined by (module,column,row)
188 // module, column,raw should follow the internal PHOS convention:
189 // module 1:5, column 1:56, row 1:64
190 // if CBD instance exists, the value is taken from CDB.
191 // Otherwise it is an ideal one
194 return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
196 return 0.0; // default EMC ADC pedestal
199 //________________________________________________________________
200 void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
202 // Set EMC calibration coefficient for (module,column,row)
205 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
207 fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
210 //________________________________________________________________
211 void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
213 // Set EMC pedestal for (module,column,row)
216 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
218 fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
221 //________________________________________________________________
222 Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
224 // Return CPV calibration coefficient
225 // for channel defined by (module,column,row)
226 // module, column,raw should follow the internal CPV convention:
227 // module 1:5, column 1:56, row 1:128
228 // if CBD instance exists, the value is taken from CDB.
229 // Otherwise it is an ideal one
232 return fCalibDataCpv->GetADCchannelCpv(module,column,row);
234 return 0.0012; // default width of one ADC channel in CPV arbitrary units
237 //________________________________________________________________
238 Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
240 // Return CPV pedestal
241 // for channel defined by (module,column,row)
242 // module, column,raw should follow the internal CPV convention:
243 // module 1:5, column 1:56, row 1:128
244 // if CBD instance exists, the value is taken from CDB.
245 // Otherwise it is an ideal one
248 return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
250 return 0.012; // default CPV ADC pedestal
253 //________________________________________________________________
254 void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
256 // Set CPV calibration coefficient for (module,column,row)
259 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
261 fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
264 //________________________________________________________________
265 void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
267 // Set CPV pedestal for (module,column,row)
270 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
272 fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
275 //________________________________________________________________
276 void AliPHOSCalibData::RandomEmc(Float_t ccMin, Float_t ccMax)
278 // Create decalibrated EMC with calibration coefficients and pedestals
279 // randomly distributed within hard-coded limits
280 // Default spread of calibration parameters is Cmax/Cmin = 4, (Cmax-Cmin)/2 = 1
282 if(fCalibDataEmc) delete fCalibDataEmc;
283 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
286 rn.SetSeed(0); //the seed is set to the current machine clock
288 Float_t adcChannelEmc,adcPedestalEmc;
290 for(Int_t module=1; module<6; module++) {
291 for(Int_t column=1; column<57; column++) {
292 for(Int_t row=1; row<65; row++) {
293 adcChannelEmc =rn.Uniform(ccMin,ccMax);
294 adcPedestalEmc=rn.Uniform(0.0,0.0); // 0 spread of pedestals
295 fCalibDataEmc->SetADCchannelEmc(module,column,row,adcChannelEmc);
296 fCalibDataEmc->SetADCpedestalEmc(module,column,row,adcPedestalEmc);
303 //________________________________________________________________
304 void AliPHOSCalibData::RandomCpv(Float_t ccMin, Float_t ccMax)
306 // Create decalibrated CPV with calibration coefficients and pedestals
307 // randomly distributed within hard-coded limits
308 // Default spread of calibration parameters is 0.0012 +- 25%
310 if(fCalibDataCpv) delete fCalibDataCpv;
311 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
314 rn.SetSeed(0); //the seed is set to the current machine clock
316 Float_t adcChannelCpv,adcPedestalCpv;
318 for(Int_t module=1; module<6; module++) {
319 for(Int_t column=1; column<57; column++) {
320 for(Int_t row=1; row<129; row++) {
321 adcChannelCpv =rn.Uniform(ccMin,ccMax);
322 adcPedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
323 fCalibDataCpv->SetADCchannelCpv(module,column,row,adcChannelCpv);
324 fCalibDataCpv->SetADCpedestalCpv(module,column,row,adcPedestalCpv);