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 ///////////////////////////////////////////////////////////////////////////////
26 #include "AliPHOSCalibData.h"
27 #include "AliCDBManager.h"
28 #include "AliCDBStorage.h"
30 #include "AliCDBEntry.h"
32 ClassImp(AliPHOSCalibData)
34 //________________________________________________________________
35 AliPHOSCalibData::AliPHOSCalibData():
36 TNamed(), fCalibDataEmc(0x0), fCalibDataCpv(0x0)
38 // Default constructor
40 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
41 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
45 //________________________________________________________________
46 AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
47 TNamed("phosCalib","PHOS Calibration Data Manager"),
48 fCalibDataEmc(0x0), fCalibDataCpv(0x0)
52 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
53 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
55 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
57 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
59 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
61 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
65 //________________________________________________________________
66 AliPHOSCalibData::~AliPHOSCalibData()
72 //________________________________________________________________
73 void AliPHOSCalibData::Reset()
75 // Set all pedestals to 0 and all ADC channels to 1
77 fCalibDataEmc->Reset();
78 fCalibDataCpv->Reset();
81 //________________________________________________________________
82 void AliPHOSCalibData::Print(Option_t *option) const
84 if (fCalibDataEmc) fCalibDataEmc->Print(option);
85 if (fCalibDataCpv) fCalibDataCpv->Print(option);
88 //________________________________________________________________
89 void AliPHOSCalibData::CreateNew()
91 if(fCalibDataEmc) delete fCalibDataEmc;
92 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
94 if(fCalibDataCpv) delete fCalibDataCpv;
95 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
99 //________________________________________________________________
100 Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
103 if(!fCalibDataEmc) return kFALSE;
105 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
107 storage = AliCDBManager::Instance()->GetDefaultStorage();
110 AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
111 storage->Put(fCalibDataEmc,id, md);
119 //________________________________________________________________
120 Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
123 if(!fCalibDataCpv) return kFALSE;
125 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
127 storage = AliCDBManager::Instance()->GetDefaultStorage();
130 AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
131 storage->Put(fCalibDataCpv,id, md);
139 //________________________________________________________________
140 Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
142 //module, column,raw should follow the internal PHOS convention:
143 //module 1:5, column 1:56, row 1:64
146 return fCalibDataEmc->GetADCchannelEmc(module,column,row);
148 return 0.0015; // default width of one EMC ADC channel in GeV
151 Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
154 return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
156 return 0.005; // default EMC ADC pedestal
159 void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
162 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
164 fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
167 void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
170 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
172 fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
175 //________________________________________________________________
176 Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
178 //module, column,raw should follow the internal CPV convention:
179 //module 1:5, column 1:64, row 1:128
182 return fCalibDataCpv->GetADCchannelCpv(module,column,row);
184 return 0.0012; // default width of one ADC channel in CPV 'popugais'
187 Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
190 return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
192 return 0.012; // default CPV ADC pedestal
195 void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
198 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
200 fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
203 void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
206 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
208 fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
211 //________________________________________________________________
212 void AliPHOSCalibData::RandomEmc()
215 if(fCalibDataEmc) delete fCalibDataEmc;
216 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
219 rn.SetSeed(0); //the seed is set to the current machine clock
221 Float_t ADCchanelEmc,ADCpedestalEmc;
223 for(Int_t module=1; module<6; module++) {
224 for(Int_t column=1; column<57; column++) {
225 for(Int_t row=1; row<65; row++) {
226 ADCchanelEmc=rn.Uniform(0.00075,0.00375); // Cmax/Cmin = 5, (Cmax-Cmin)/2 = 0.0015
227 ADCpedestalEmc=rn.Uniform(0.0045,0.0055); //+-10% spread of pedestals from 0.005
228 fCalibDataEmc->SetADCchannelEmc(module,column,row,ADCchanelEmc);
229 fCalibDataEmc->SetADCpedestalEmc(module,column,row,ADCpedestalEmc);
236 //________________________________________________________________
237 void AliPHOSCalibData::RandomCpv()
240 if(fCalibDataCpv) delete fCalibDataCpv;
241 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
244 rn.SetSeed(0); //the seed is set to the current machine clock
246 Float_t ADCchanelCpv,ADCpedestalCpv;
248 for(Int_t module=1; module<6; module++) {
249 for(Int_t column=1; column<65; column++) {
250 for(Int_t row=1; row<129; row++) {
251 ADCchanelCpv=TMath::Abs(rn.Uniform(0.0009,0.0015)); // 0.0012 +- 25%
252 ADCpedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
253 fCalibDataCpv->SetADCchannelCpv(module,column,row,ADCchanelCpv);
254 fCalibDataCpv->SetADCpedestalCpv(module,column,row,ADCpedestalCpv);