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 **************************************************************************/
16 ///////////////////////////////////////////////////////////////////////////////
18 // class for PHOS EmCal calibration //
20 ///////////////////////////////////////////////////////////////////////////////
22 #include "AliPHOSEmcCalibData.h"
24 ClassImp(AliPHOSEmcCalibData)
26 //________________________________________________________________
27 AliPHOSEmcCalibData::AliPHOSEmcCalibData():
31 // Default constructor
35 //________________________________________________________________
36 AliPHOSEmcCalibData::AliPHOSEmcCalibData(const char* name):
41 TString namst = "Calib_";
43 SetName(namst.Data());
44 SetTitle(namst.Data());
48 //________________________________________________________________
49 AliPHOSEmcCalibData::AliPHOSEmcCalibData(const AliPHOSEmcCalibData& calibda) :
51 fSampleTimeStep(calibda.fSampleTimeStep)
54 SetName(calibda.GetName());
55 SetTitle(calibda.GetName());
57 for(Int_t module=0; module<5; module++) {
58 for(Int_t column=0; column<56; column++) {
59 for(Int_t row=0; row<64; row++) {
60 fADCchannelEmc[module][column][row] = calibda.fADCchannelEmc[module][column][row];
61 fADCpedestalEmc[module][column][row] = calibda.fADCpedestalEmc[module][column][row];
62 fHighLowRatioEmc[module][column][row] = calibda.fHighLowRatioEmc[module][column][row];
63 fTimeShiftEmc[module][column][row] = calibda.fTimeShiftEmc[module][column][row];
64 fAltroOffsets[module][column][row] = calibda.fAltroOffsets[module][column][row];
70 //________________________________________________________________
71 AliPHOSEmcCalibData &AliPHOSEmcCalibData::operator =(const AliPHOSEmcCalibData& calibda)
73 // assignment operator
75 if(this != &calibda) {
77 SetName(calibda.GetName());
78 SetTitle(calibda.GetName());
80 for(Int_t module=0; module<5; module++) {
81 for(Int_t column=0; column<56; column++) {
82 for(Int_t row=0; row<64; row++) {
83 fADCchannelEmc[module][column][row] = calibda.fADCchannelEmc[module][column][row];
84 fADCpedestalEmc[module][column][row] = calibda.fADCpedestalEmc[module][column][row];
85 fHighLowRatioEmc[module][column][row] = calibda.fHighLowRatioEmc[module][column][row];
86 fTimeShiftEmc[module][column][row] = calibda.fTimeShiftEmc[module][column][row];
87 fAltroOffsets[module][column][row] = calibda.fAltroOffsets[module][column][row];
96 //________________________________________________________________
97 AliPHOSEmcCalibData::~AliPHOSEmcCalibData()
102 //________________________________________________________________
103 void AliPHOSEmcCalibData::Reset()
105 // Set all pedestals and all ADC channels to its ideal values = 5. (MeV/ADC)
107 for (Int_t module=0; module<5; module++){
108 for (Int_t column=0; column<56; column++){
109 for (Int_t row=0; row<64; row++){
110 fADCpedestalEmc[module][column][row] = 0.;
111 fADCchannelEmc[module][column][row] = 0.005;
112 fHighLowRatioEmc[module][column][row] = 16. ;
113 fTimeShiftEmc[module][column][row] = 0. ;
114 fAltroOffsets[module][column][row] = 0 ;
118 fSampleTimeStep=100.e-9 ; //100 ns
121 //________________________________________________________________
122 void AliPHOSEmcCalibData::Print(Option_t *option) const
124 // Print tables of pedestals and ADC channels
126 if (strstr(option,"ped")) {
127 printf("\n ---- EMC Pedestal values ----\n\n");
128 for (Int_t module=0; module<5; module++){
129 printf("============== Module %d\n",module+1);
130 for (Int_t column=0; column<56; column++){
131 for (Int_t row=0; row<64; row++){
132 printf("%4.1f",fADCpedestalEmc[module][column][row]);
139 if (strstr(option,"gain")) {
140 printf("\n ---- EMC ADC channel values ----\n\n");
141 for (Int_t module=0; module<5; module++){
142 printf("============== Module %d\n",module+1);
143 for (Int_t column=0; column<56; column++){
144 for (Int_t row=0; row<64; row++){
145 printf("%4.1f",fADCchannelEmc[module][column][row]);
152 if (strstr(option,"hilo")) {
153 printf("\n ---- EMC High/Low ratio ----\n\n");
154 for (Int_t module=0; module<5; module++){
155 printf("============== Module %d\n",module+1);
156 for (Int_t column=0; column<56; column++){
157 for (Int_t row=0; row<64; row++){
158 printf("%4.1f",fHighLowRatioEmc[module][column][row]);
164 if (strstr(option,"time")) {
165 printf("\n ---- EMC t0 shifts ----\n\n");
166 for (Int_t module=0; module<5; module++){
167 printf("============== Module %d\n",module+1);
168 for (Int_t column=0; column<56; column++){
169 for (Int_t row=0; row<64; row++){
170 printf("%6.3e",fTimeShiftEmc[module][column][row]);
176 if (strstr(option,"altro")) {
177 printf("\n ---- EMC altro offsets ----\n\n");
178 for (Int_t module=0; module<5; module++){
179 printf("============== Module %d\n",module+1);
180 for (Int_t column=0; column<56; column++){
181 for (Int_t row=0; row<64; row++){
182 printf("%5d",fAltroOffsets[module][column][row]);
190 //________________________________________________________________
191 Float_t AliPHOSEmcCalibData::GetADCchannelEmc(Int_t module, Int_t column, Int_t row) const
193 //Return EMC calibration coefficient
194 //module, column,raw should follow the internal PHOS convention:
195 //module 1:5, column 1:56, row 1:64
197 return fADCchannelEmc[module-1][column-1][row-1];
200 //________________________________________________________________
201 Float_t AliPHOSEmcCalibData::GetADCpedestalEmc(Int_t module, Int_t column, Int_t row) const
203 //Return EMC pedestal
204 //module, column,raw should follow the internal PHOS convention:
205 //module 1:5, column 1:56, row 1:64
207 return fADCpedestalEmc[module-1][column-1][row-1];
210 //________________________________________________________________
211 Float_t AliPHOSEmcCalibData::GetHighLowRatioEmc(Int_t module, Int_t column, Int_t row) const
213 //Return EMC pedestal
214 //module, column,raw should follow the internal PHOS convention:
215 //module 1:5, column 1:56, row 1:64
217 return fHighLowRatioEmc[module-1][column-1][row-1];
220 //________________________________________________________________
221 Float_t AliPHOSEmcCalibData::GetTimeShiftEmc(Int_t module, Int_t column, Int_t row) const
223 //Return EMC pedestal
224 //module, column,raw should follow the internal PHOS convention:
225 //module 1:5, column 1:56, row 1:64
227 return fTimeShiftEmc[module-1][column-1][row-1];
229 //________________________________________________________________
230 Float_t AliPHOSEmcCalibData::GetSampleTimeStep()const
232 //Returns conversion coefficient from ALTRO smaple time step and secods
233 //Negative value not used in reconstruction (conversion taken from TRU trailer)
234 //and only in raw simulation
235 return fSampleTimeStep ;
237 //________________________________________________________________
238 Int_t AliPHOSEmcCalibData::GetAltroOffsetEmc(Int_t module, Int_t column, Int_t row) const
240 //Return EMC altro offsets
241 //module, column,raw should follow the internal PHOS convention:
242 //module 1:5, column 1:56, row 1:64
244 return fAltroOffsets[module-1][column-1][row-1];
246 //________________________________________________________________
247 void AliPHOSEmcCalibData::SetADCchannelEmc(Int_t module, Int_t column, Int_t row, Float_t value)
249 //Set EMC calibration coefficient
250 //module, column,raw should follow the internal PHOS convention:
251 //module 1:5, column 1:56, row 1:64
253 fADCchannelEmc[module-1][column-1][row-1] = value;
256 //________________________________________________________________
257 void AliPHOSEmcCalibData::SetADCpedestalEmc(Int_t module, Int_t column, Int_t row, Float_t value)
260 //module, column,raw should follow the internal PHOS convention:
261 //module 1:5, column 1:56, row 1:64
262 fADCpedestalEmc[module-1][column-1][row-1] = value;
265 //________________________________________________________________
266 void AliPHOSEmcCalibData::SetHighLowRatioEmc(Int_t module, Int_t column, Int_t row, Float_t value)
269 //module, column,raw should follow the internal PHOS convention:
270 //module 1:5, column 1:56, row 1:64
271 fHighLowRatioEmc[module-1][column-1][row-1] = value;
273 //________________________________________________________________
274 void AliPHOSEmcCalibData::SetTimeShiftEmc(Int_t module, Int_t column, Int_t row, Float_t value)
277 //module, column,raw should follow the internal PHOS convention:
278 //module 1:5, column 1:56, row 1:64
279 fTimeShiftEmc[module-1][column-1][row-1] = value;
281 //________________________________________________________________
282 void AliPHOSEmcCalibData::SetAltroOffsetEmc(Int_t module, Int_t column, Int_t row, Int_t value)
285 //module, column,raw should follow the internal PHOS convention:
286 //module 1:5, column 1:56, row 1:64
287 fAltroOffsets[module-1][column-1][row-1] = value;
289 //________________________________________________________________
290 void AliPHOSEmcCalibData::SetSampleTimeStep(Float_t step)
292 //Sets conversion coefficient from ALTRO smaple time step and secods
293 //Negative value not used in reconstruction (conversion taken from TRU trailer)
294 //and only in raw simulation
295 fSampleTimeStep = step ;