add getter for the array with digits and initialization inside, just for analysis...
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALCalibData.cxx
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 ///*-- Author: Yves Schutz (SUBATECH)
20 //           : Aleksei Pavlinov (WSU); Jun 30, 2006 - ALICE numbering scheme
21 //           : Add decalibration and time calibration arrays: Jul 21, 2011 (GCB)
22 ///////////////////////////////////////////////////////////////////////////////
23 //                                                                           //
24 // class for EMCAL calibration                                               //
25 //                                                                           //
26 ///////////////////////////////////////////////////////////////////////////////
27
28 #include <TMath.h>
29
30 #include "AliEMCALCalibData.h"
31
32 ClassImp(AliEMCALCalibData)
33
34 //________________________________________________________________
35 AliEMCALCalibData::AliEMCALCalibData()
36 {
37   // Default constructor
38   Reset();
39 }
40
41 //________________________________________________________________
42 AliEMCALCalibData::AliEMCALCalibData(const char* name)
43 {
44   // Constructor
45   TString namst = "Calib_";
46   namst += name;
47   SetName(namst.Data());
48   SetTitle(namst.Data());
49   Reset();
50 }
51
52 //________________________________________________________________
53 AliEMCALCalibData::AliEMCALCalibData(const AliEMCALCalibData& calibda) :
54   TNamed(calibda)
55 {
56   // copy constructor
57   SetName(calibda.GetName());
58   SetTitle(calibda.GetName());
59   Reset();
60
61   Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
62   Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
63   Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
64   Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows;    //12 - Modules 11 and 12 are half modules
65
66   for(Int_t supermodule = 0; supermodule < nSMod; supermodule++) {
67     
68     if(supermodule >= 10)
69       nRow = nRow2;
70     
71     for(Int_t column = 0; column<nCol; column++) {
72       
73       for(Int_t row = 0; row<nRow; row++) {
74         
75         fADCchannel[supermodule][column][row] = 
76         calibda.GetADCchannel(supermodule,column,row);
77         
78         fADCchannelDecal[supermodule][column][row] = 
79         calibda.GetADCchannelDecal(supermodule,column,row);
80         
81         fADCpedestal[supermodule][column][row] = 
82         calibda.GetADCpedestal(supermodule,column,row);
83         
84         fTimeChannel[supermodule][column][row] = 
85         calibda.GetTimeChannel(supermodule,column,row);
86         
87         fTimeChannelDecal[supermodule][column][row] = 
88         calibda.GetTimeChannelDecal(supermodule,column,row);
89         
90         
91       }
92     }
93   }
94 }
95
96
97 //________________________________________________________________
98 AliEMCALCalibData &AliEMCALCalibData::operator =(const AliEMCALCalibData& calibda)
99 {
100   // assignment operator
101   SetName(calibda.GetName());
102   SetTitle(calibda.GetName());
103   Reset();
104
105   Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
106   Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
107   Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
108   Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //12 - Modules 11 and 12 are half modules
109
110   for(Int_t supermodule = 0; supermodule < nSMod; supermodule++) {
111     
112     if(supermodule >= 10)
113       nRow = nRow2;
114     
115     for(Int_t column = 0; column<nCol; column++) {
116       
117       for(Int_t row = 0; row<nRow; row++) {
118         
119         fADCchannel[supermodule][column][row] = 
120         calibda.GetADCchannel(supermodule,column,row);
121         
122         fADCchannelDecal[supermodule][column][row] = 
123         calibda.GetADCchannelDecal(supermodule,column,row);
124         
125         fADCpedestal[supermodule][column][row] = 
126         calibda.GetADCpedestal(supermodule,column,row);
127         
128         fTimeChannel[supermodule][column][row] = 
129         calibda.GetTimeChannel(supermodule,column,row);
130
131         fTimeChannelDecal[supermodule][column][row] = 
132         calibda.GetTimeChannelDecal(supermodule,column,row);
133         
134       }
135     }
136   }
137   return *this;
138 }
139
140 //________________________________________________________________
141 AliEMCALCalibData::~AliEMCALCalibData()
142 {
143   // Destructor
144 }
145
146 //________________________________________________________________
147 void AliEMCALCalibData::Reset()
148 {
149   // Set all pedestals to 0 and all ADC channels widths to 1
150   //memset(fADCchannel ,1,12*48*24*sizeof(Float_t));
151   //memset(fADCpedestal,0,12*48*24*sizeof(Float_t));
152
153   Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
154   Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
155   Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
156   Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //12 - Modules 11 and 12 are half modules
157   for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
158     if(supermodule >= 10)
159       nRow = nRow2;
160     for (Int_t column=0; column < nCol; column++){
161       
162       for (Int_t row = 0; row < nRow; row++){
163         
164         fADCpedestal     [supermodule][column][row]=0.;
165         
166         fADCchannelDecal [supermodule][column][row]=1.;
167         fADCchannel      [supermodule][column][row]=1.;
168         
169         fTimeChannel     [supermodule][column][row]=0.;
170         fTimeChannelDecal[supermodule][column][row]=0.;
171
172       }
173     }
174   } 
175 }
176
177 //________________________________________________________________
178 void  AliEMCALCalibData::Print(Option_t *option) const
179 {
180   // Print tables of pedestals and ADC channels widths
181   // options are: "gain", "ped", "decal", "time", "all"
182   
183   Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
184   Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
185   Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
186   Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //12 - Modules 11 and 12 are half modules
187   if (strstr(option,"ped") || strstr(option,"all")) {
188     printf("\n  ----    Pedestal values ----\n\n");
189     for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
190       if(supermodule >= 10)
191         nRow = nRow2;
192       printf("============== Supermodule %d\n",supermodule+1);
193       for (Int_t column=0; column<nCol; column++){
194         for (Int_t row=0; row<nRow; row++){
195           printf(" %2.4f ",fADCpedestal[supermodule][column][row]);
196         }
197         printf("\n");
198       }
199     } 
200   }
201   
202   if (strstr(option,"gain") || strstr(option,"all")) {
203     printf("\n  ----    ADC channel values      ----\n\n");
204     for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
205       if(supermodule >= 10) 
206         nRow = nRow2;
207       printf("============== Supermodule %d\n",supermodule+1);
208       for (Int_t column=0; column<nCol; column++){
209         for (Int_t row=0; row<nRow; row++){
210           printf(" %2.4f ",fADCchannel[supermodule][column][row]);
211         }
212         printf("\n");
213       }
214     }   
215   }
216   
217   if (strstr(option,"adcdecal") || strstr(option,"all")) {
218     printf("\n  ----    ADC decalibration channel values        ----\n\n");
219     for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
220       if(supermodule >= 10) 
221         nRow = nRow2;
222       printf("============== Supermodule %d\n",supermodule+1);
223       for (Int_t column=0; column<nCol; column++){
224         for (Int_t row=0; row<nRow; row++){
225           printf(" %2.4f ",fADCchannelDecal[supermodule][column][row]);
226         }
227         printf("\n");
228       }
229     }   
230   }
231   
232   if (strstr(option,"time") || strstr(option,"all")) {
233     printf("\n  ----    time channel values     ----\n\n");
234     for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
235       if(supermodule >= 10) 
236         nRow = nRow2;
237       printf("============== Supermodule %d\n",supermodule+1);
238       for (Int_t column=0; column<nCol; column++){
239         for (Int_t row=0; row<nRow; row++){
240           printf(" %2.4f ",fTimeChannel[supermodule][column][row]);
241         }
242         printf("\n");
243       }
244     }   
245   }
246   
247   
248   if (strstr(option,"time") || strstr(option,"all")) {
249     printf("\n  ----    time decalibration channel values       ----\n\n");
250     for (Int_t supermodule=0; supermodule<nSMod; supermodule++){
251       if(supermodule >= 10) 
252         nRow = nRow2;
253       printf("============== Supermodule %d\n",supermodule+1);
254       for (Int_t column=0; column<nCol; column++){
255         for (Int_t row=0; row<nRow; row++){
256           printf(" %2.4f ",fTimeChannelDecal[supermodule][column][row]);
257         }
258         printf("\n");
259       }
260     }   
261   }
262   
263   
264 }
265
266 //________________________________________________________________
267 Float_t AliEMCALCalibData::GetADCchannel(Int_t supermodule, Int_t column, Int_t row) const
268 {
269   // Set ADC channel witdth values
270   // All indexes start from 0!
271   // Supermodule, column,raw should follow the ALICE convention:
272   // supermodule 0:11, column 0:47, row 0:23
273
274   return fADCchannel[supermodule][column][row];
275 }
276
277 //________________________________________________________________
278 Float_t AliEMCALCalibData::GetADCchannelDecal(Int_t supermodule, Int_t column, Int_t row) const
279 {
280   // Set ADC channel decalibration witdth values
281   // All indexes start from 0!
282   // Supermodule, column,raw should follow the ALICE convention:
283   // supermodule 0:11, column 0:47, row 0:23
284   
285   return fADCchannelDecal[supermodule][column][row];
286 }
287
288 //________________________________________________________________
289 Float_t AliEMCALCalibData::GetADCpedestal(Int_t supermodule, Int_t column, Int_t row) const
290 {
291   // Get ADC pedestal values
292  return fADCpedestal[supermodule][column][row];
293 }
294
295 //________________________________________________________________
296 Float_t AliEMCALCalibData::GetTimeChannel(Int_t supermodule, Int_t column, Int_t row) const
297 {
298   // Set channel time witdth values
299   return fTimeChannel[supermodule][column][row];
300 }
301
302 //________________________________________________________________
303 Float_t AliEMCALCalibData::GetTimeChannelDecal(Int_t supermodule, Int_t column, Int_t row) const
304 {
305   // Set channel time witdth values
306   return fTimeChannelDecal[supermodule][column][row];
307 }
308
309 //________________________________________________________________
310 void AliEMCALCalibData::SetADCchannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
311
312   // Set ADC channel width values
313   fADCchannel[supermodule][column][row] = value;
314 }
315
316 //________________________________________________________________
317 void AliEMCALCalibData::SetADCchannelDecal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
318
319   // Set ADC channel width values
320   fADCchannelDecal[supermodule][column][row] = value;
321 }
322
323 //________________________________________________________________
324 void AliEMCALCalibData::SetADCpedestal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
325 {
326   // Set ADC pedestal values
327   fADCpedestal[supermodule][column][row] = value;
328 }
329
330 //________________________________________________________________
331 void AliEMCALCalibData::SetTimeChannel(Int_t supermodule, Int_t column, Int_t row, Float_t value)
332 {
333   // Set ADC pedestal values
334   fTimeChannel[supermodule][column][row] = value;
335 }
336
337 //________________________________________________________________
338 void AliEMCALCalibData::SetTimeChannelDecal(Int_t supermodule, Int_t column, Int_t row, Float_t value)
339 {
340   // Set ADC pedestal values
341   fTimeChannelDecal[supermodule][column][row] = value;
342 }
343
344
345 //________________________________________________________________
346 void AliEMCALCalibData::Fill(const AliEMCALCalibData *cd1, const AliEMCALCalibData *cd2, Bool_t print)
347 {
348   // Fill for (relative) recalibration to undo calibration 
349   // from 1 and apply calibration from 2.
350
351   Int_t nSMod = AliEMCALGeoParams::fgkEMCALModules; //12
352   Int_t nCol  = AliEMCALGeoParams::fgkEMCALCols;    //48
353   Int_t nRow  = AliEMCALGeoParams::fgkEMCALRows;    //24
354   Int_t nRow2 = AliEMCALGeoParams::fgkEMCALRows/2;  //12 - Modules 11 and 12 are half modules
355   for (Int_t supermodule=0; supermodule<nSMod; supermodule++) {
356     if(supermodule >= 10)
357       nRow = nRow2;
358     for (Int_t column=0; column<nCol; column++) {
359       for (Int_t row=0; row<nRow; row++){
360         Double_t adc1 = cd1->GetADCchannel(supermodule, column, row);
361         Double_t ped1 = cd1->GetADCpedestal(supermodule, column, row);
362         Double_t adc2 = cd2->GetADCchannel(supermodule, column, row);
363         Double_t ped2 = cd2->GetADCpedestal(supermodule, column, row);
364         Double_t adc = adc2/adc1;
365         Double_t ped = ped2-adc*ped1;
366         if (print)
367           printf("%d %d %d: %f %f %f %f -> %f %f\n",supermodule,column,row,adc1,ped1,adc2,ped2,adc,ped);
368         fADCchannel [supermodule][column][row]=adc;
369         fADCpedestal[supermodule][column][row]=ped;
370       }
371     }
372   } 
373 }
374