git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Additional protection
[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 "TMath.h"
25#include "TRandom.h"
52783dbc 26#include "AliPHOSCalibData.h"
fc6706cb 27#include "AliCDBManager.h"
28#include "AliCDBStorage.h"
29#include "AliCDBId.h"
30#include "AliCDBEntry.h"
52783dbc 31
32ClassImp(AliPHOSCalibData)
33
34//________________________________________________________________
fc6706cb 35 AliPHOSCalibData::AliPHOSCalibData():
36 TNamed(), fCalibDataEmc(0x0), fCalibDataCpv(0x0)
52783dbc 37{
fc6706cb 38 // Default constructor
39
df0d5123 40 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
41 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
fc6706cb 42
abcff5fc 43 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data());
44 if(entryEmc)
45 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
46
47 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data());
48 if(entryCpv)
49 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
50
52783dbc 51}
52
53//________________________________________________________________
fc6706cb 54AliPHOSCalibData::AliPHOSCalibData(Int_t runNumber) :
55 TNamed("phosCalib","PHOS Calibration Data Manager"),
56 fCalibDataEmc(0x0), fCalibDataCpv(0x0)
52783dbc 57{
58 // Constructor
fc6706cb 59
df0d5123 60 fEmcDataPath="PHOS/Calib/EmcGainPedestals";
61 fCpvDataPath="PHOS/Calib/CpvGainPedestals";
fc6706cb 62
63 AliCDBEntry* entryEmc = AliCDBManager::Instance()->Get(fEmcDataPath.Data(),runNumber);
64 if(entryEmc)
65 fCalibDataEmc = (AliPHOSEmcCalibData*)entryEmc->GetObject();
66
67 AliCDBEntry* entryCpv = AliCDBManager::Instance()->Get(fCpvDataPath.Data(),runNumber);
68 if(entryCpv)
69 fCalibDataCpv = (AliPHOSCpvCalibData*)entryCpv->GetObject();
52783dbc 70
52783dbc 71}
72
73//________________________________________________________________
74AliPHOSCalibData::~AliPHOSCalibData()
75{
76 // Destructor
fc6706cb 77
52783dbc 78}
79
80//________________________________________________________________
81void AliPHOSCalibData::Reset()
82{
83 // Set all pedestals to 0 and all ADC channels to 1
fc6706cb 84
85 fCalibDataEmc->Reset();
86 fCalibDataCpv->Reset();
52783dbc 87}
88
89//________________________________________________________________
90void AliPHOSCalibData::Print(Option_t *option) const
91{
df0d5123 92 if (fCalibDataEmc) fCalibDataEmc->Print(option);
93 if (fCalibDataCpv) fCalibDataCpv->Print(option);
fc6706cb 94}
95
96//________________________________________________________________
97void AliPHOSCalibData::CreateNew()
98{
99 if(fCalibDataEmc) delete fCalibDataEmc;
100 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
101
102 if(fCalibDataCpv) delete fCalibDataCpv;
103 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
104
105}
106
107//________________________________________________________________
108Bool_t AliPHOSCalibData::WriteEmc(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
109{
110
111 if(!fCalibDataEmc) return kFALSE;
112
113 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
df0d5123 114 if(!storage)
115 storage = AliCDBManager::Instance()->GetDefaultStorage();
116
fc6706cb 117 if(storage) {
118 AliCDBId id(fEmcDataPath.Data(),firstRun,lastRun);
119 storage->Put(fCalibDataEmc,id, md);
120 return kTRUE;
52783dbc 121 }
fc6706cb 122 else
123 return kFALSE;
52783dbc 124
fc6706cb 125}
126
127//________________________________________________________________
128Bool_t AliPHOSCalibData::WriteCpv(Int_t firstRun, Int_t lastRun, AliCDBMetaData *md)
129{
130
131 if(!fCalibDataCpv) return kFALSE;
132
133 AliCDBStorage* storage = AliCDBManager::Instance()->GetSpecificStorage("PHOS");
df0d5123 134 if(!storage)
135 storage = AliCDBManager::Instance()->GetDefaultStorage();
136
fc6706cb 137 if(storage) {
138 AliCDBId id(fCpvDataPath.Data(),firstRun,lastRun);
139 storage->Put(fCalibDataCpv,id, md);
140 return kTRUE;
52783dbc 141 }
fc6706cb 142 else
143 return kFALSE;
144
52783dbc 145}
072de3a8 146
fc6706cb 147//________________________________________________________________
072de3a8 148Float_t AliPHOSCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
149{
150 //module, column,raw should follow the internal PHOS convention:
151 //module 1:5, column 1:56, row 1:64
152
fc6706cb 153 if(fCalibDataEmc)
154 return fCalibDataEmc->GetADCchannelEmc(module,column,row);
155 else
156 return 0.0015; // default width of one EMC ADC channel in GeV
072de3a8 157}
158
159Float_t AliPHOSCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
160{
fc6706cb 161 if(fCalibDataEmc)
162 return fCalibDataEmc->GetADCpedestalEmc(module,column,row);
163 else
164 return 0.005; // default EMC ADC pedestal
072de3a8 165}
166
167void AliPHOSCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
168{
fc6706cb 169 if(!fCalibDataEmc)
170 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
171
172 fCalibDataEmc->SetADCchannelEmc(module,column,row,value);
072de3a8 173}
174
175void AliPHOSCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
176{
fc6706cb 177 if(!fCalibDataEmc)
178 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
179
180 fCalibDataEmc->SetADCpedestalEmc(module,column,row,value);
181}
182
183//________________________________________________________________
184Float_t AliPHOSCalibData::GetADCchannelCpv(Int_t module, Int_t column, Int_t row) const
185{
186 //module, column,raw should follow the internal CPV convention:
187 //module 1:5, column 1:64, row 1:128
188
189 if(fCalibDataCpv)
190 return fCalibDataCpv->GetADCchannelCpv(module,column,row);
191 else
192 return 0.0012; // default width of one ADC channel in CPV 'popugais'
193}
194
195Float_t AliPHOSCalibData::GetADCpedestalCpv(Int_t module, Int_t column, Int_t row) const
196{
197 if(fCalibDataCpv)
198 return fCalibDataCpv->GetADCpedestalCpv(module,column,row);
199 else
200 return 0.012; // default CPV ADC pedestal
201}
202
203void AliPHOSCalibData::SetADCchannelCpv(Int_t module, Int_t column, Int_t row, Float_t value)
204{
205 if(!fCalibDataCpv)
206 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
207
208 fCalibDataCpv->SetADCchannelCpv(module,column,row,value);
209}
210
211void AliPHOSCalibData::SetADCpedestalCpv(Int_t module, Int_t column, Int_t row, Float_t value)
212{
213 if(!fCalibDataCpv)
214 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
215
216 fCalibDataCpv->SetADCpedestalCpv(module,column,row,value);
217}
218
219//________________________________________________________________
220void AliPHOSCalibData::RandomEmc()
221{
222
223 if(fCalibDataEmc) delete fCalibDataEmc;
224 fCalibDataEmc = new AliPHOSEmcCalibData("PHOS-EMC");
225
226 TRandom rn;
227 rn.SetSeed(0); //the seed is set to the current machine clock
228
229 Float_t ADCchanelEmc,ADCpedestalEmc;
230
231 for(Int_t module=1; module<6; module++) {
232 for(Int_t column=1; column<57; column++) {
233 for(Int_t row=1; row<65; row++) {
234 ADCchanelEmc=rn.Uniform(0.00075,0.00375); // Cmax/Cmin = 5, (Cmax-Cmin)/2 = 0.0015
235 ADCpedestalEmc=rn.Uniform(0.0045,0.0055); //+-10% spread of pedestals from 0.005
236 fCalibDataEmc->SetADCchannelEmc(module,column,row,ADCchanelEmc);
237 fCalibDataEmc->SetADCpedestalEmc(module,column,row,ADCpedestalEmc);
238 }
239 }
240 }
241
242}
243
244//________________________________________________________________
245void AliPHOSCalibData::RandomCpv()
246{
247
248 if(fCalibDataCpv) delete fCalibDataCpv;
249 fCalibDataCpv = new AliPHOSCpvCalibData("PHOS-CPV");
250
251 TRandom rn;
252 rn.SetSeed(0); //the seed is set to the current machine clock
253
254 Float_t ADCchanelCpv,ADCpedestalCpv;
255
256 for(Int_t module=1; module<6; module++) {
257 for(Int_t column=1; column<65; column++) {
258 for(Int_t row=1; row<129; row++) {
259 ADCchanelCpv=TMath::Abs(rn.Uniform(0.0009,0.0015)); // 0.0012 +- 25%
260 ADCpedestalCpv=rn.Uniform(0.0048,0.0192); // Ped[max]/Ped[min] = 4, <Ped> = 0.012
261 fCalibDataCpv->SetADCchannelCpv(module,column,row,ADCchanelCpv);
262 fCalibDataCpv->SetADCpedestalCpv(module,column,row,ADCpedestalCpv);
263 }
264 }
265 }
072de3a8 266}