Fix for 64-bit platforms
[u/mrichter/AliRoot.git] / RAW / AliCaloAltroMapping.cxx
CommitLineData
d84933b0 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/* $Id$ */
16
17// This class handles the mapping of the Altro channels in the PHOS/EMCAL
18// The mapping is read from an external mapping files
19// Author: C.Cheshkov
20
21/// Exported from PHOS to be used also by EMCAL
22/// November 2006 Gustavo Conesa Balbastre
23
24#include "AliCaloAltroMapping.h"
25#include "AliLog.h"
26#include <Riostream.h>
27//#include <stdlib.h>
28
29
30ClassImp(AliCaloAltroMapping)
31
32//_____________________________________________________________________________
33AliCaloAltroMapping::AliCaloAltroMapping(const char *mappingFile):
34 AliAltroMapping(mappingFile),
35 fMinRow(0),
36 fMaxRow(0),
37 fMinCol(0),
38 fMaxCol(0),
39 fMapping(NULL),
40 fInvMappingLow(NULL),
41 fInvMappingHigh(NULL)
42{
43 // Constructor
44 ReadMapping();
45 CloseMappingFile();
46}
47
48//_____________________________________________________________________________
49AliCaloAltroMapping::~AliCaloAltroMapping()
50{
51 // destructor
52 DeleteMappingArrays();
53}
54
55//_____________________________________________________________________________
56AliCaloAltroMapping::AliCaloAltroMapping(const AliCaloAltroMapping& mapping):
57 AliAltroMapping(mapping),
58 fMinRow(mapping.fMinRow),
59 fMaxRow(mapping.fMaxRow),
60 fMinCol(mapping.fMinCol),
61 fMaxCol(mapping.fMaxCol),
62 fMapping(mapping.fMapping),
63 fInvMappingLow(mapping.fInvMappingLow),
64 fInvMappingHigh(mapping.fInvMappingHigh)
65{
66// Copy Constructor
67
68 Fatal("AliCaloAltroMapping", "copy constructor not implemented");
69}
70
71//_____________________________________________________________________________
72AliCaloAltroMapping& AliCaloAltroMapping::operator = (const AliCaloAltroMapping& /*mapping*/)
73{
74//Assigment operator
75
76 Fatal("operator =", "assignment operator not implemented");
77 return *this;
78}
79
80//_____________________________________________________________________________
81Bool_t AliCaloAltroMapping::ReadMapping()
82{
83 // Initalizes the ALTRO mapping from a file
84 // Look at the Calo module for the format of
85 // the mapping file
86 if (!fIn) {
87 AliFatal("Mapping file has not been opened !");
88 return kFALSE;
89 }
90
91 fMinRow = 0x7fffffff;
92 fMaxRow = 0;
93 fMinCol = 0x7fffffff;
94 fMaxCol = 0;
95 fMapping = new Short_t*[fMaxHWAddress+1];
96 for (Int_t i = 0; i <= fMaxHWAddress; i++) {
97 fMapping[i] = new Short_t[3];
98 fMapping[i][0] = fMapping[i][1] = fMapping[i][2] = -1;
99 }
100
101 for(Int_t i = 0; i < fNumberOfChannels ; i++) { // 1792 = 2*896 channels connected to each RCU
102 Int_t hwAddress;
103 if (!(*fIn >> hwAddress)) {
104 AliFatal("Syntax of the mapping file is wrong !");
105 return kFALSE;
106 }
107 if (hwAddress > fMaxHWAddress) {
108 AliFatal(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
109 return kFALSE;
110 }
111 Int_t row,col,gain;
112 if (!(*fIn >> row >> col >> gain)) {
113 AliFatal("Syntax of the mapping file is wrong !");
114 return kFALSE;
115 }
116
117 if (gain < 0 || gain > 1) {
118 AliFatal(Form("Wrong gain value found (%d)! Should be 0 or 1 !",gain));
119 return kFALSE;
120 }
121
122 fMapping[hwAddress][0] = row;
123 fMapping[hwAddress][1] = col;
124 fMapping[hwAddress][2] = gain;
125
126 if (row > fMaxRow) fMaxRow = row;
127 if (row < fMinRow) fMinRow = row;
128 if (col > fMaxCol) fMaxCol = col;
129 if (col < fMinCol) fMinCol = col;
130
131 }
132
133 fInvMappingLow = new Short_t*[fMaxRow - fMinRow + 1];
134 fInvMappingHigh = new Short_t*[fMaxRow - fMinRow + 1];
135 for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++) {
136 fInvMappingLow[i] = new Short_t[fMaxCol - fMinCol + 1];
137 fInvMappingHigh[i] = new Short_t[fMaxCol - fMinCol + 1];
138 for (Int_t j = 0; j <= (fMaxCol - fMinCol); j++) {
139 fInvMappingLow[i][j] = -1;
140 fInvMappingHigh[i][j] = -1;
141 }
142 }
143
144 for(Int_t i = 0; i <= fMaxHWAddress; i++) {
145 Int_t row = fMapping[i][0];
146 Int_t col = fMapping[i][1];
147 Int_t gain = fMapping[i][2];
148 if(row != -1 && col != -1) {
149 if (gain == 0)
150 fInvMappingLow[row-fMinRow][col-fMinCol] = i;
151 if (gain == 1)
152 fInvMappingHigh[row-fMinRow][col-fMinCol] = i;
153 }
154 }
155
156 return kTRUE;
157}
158
159//_____________________________________________________________________________
160Int_t AliCaloAltroMapping::GetHWAddress(Int_t row, Int_t col, Int_t gain) const
161{
162 // Get the content of the mapping array
163 // return -1 in case there is no hardware
164 // adress defined for these row-column-gain
165 if (!fInvMappingLow || !fInvMappingHigh) {
166 AliWarning("Mapping array was not initalized correctly !");
167 return -1;
168 }
169 if (row < fMinRow || row > fMaxRow) {
170 AliWarning(Form("Index of row (%d) outside the range (%d -> %d) !",row,fMinRow,fMaxRow));
171 return -1;
172 }
173 if (col < fMinCol || col > fMaxCol) {
174 AliWarning(Form("Index of column (%d) outside the range (0 -> %d) !",col,fMinCol,fMaxCol));
175 return -1;
176 }
177 if (gain < 0 || gain > 1) {
178 AliWarning(Form("Invalid gain (%d)! Should be 0 or 1 !",gain));
179 return -1;
180 }
181 Int_t hwAddress = -1;
182 if (gain == 0)
183 hwAddress = fInvMappingLow[row-fMinRow][col-fMinCol];
184 if (gain == 1)
185 hwAddress = fInvMappingHigh[row-fMinRow][col-fMinCol];
186
187 if (hwAddress == -1)
188 AliWarning(Form("Hardware (ALTRO) adress is not defined for these row (%d), column (%d) and gain (%d) !",row,col,gain));
189
190 return hwAddress;
191}
192
193//_____________________________________________________________________________
194Int_t AliCaloAltroMapping::GetPadRow(Int_t hwAddress) const
195{
196 // Return the row index
197 // Note the difference w.r.t to the base class notation
198 if (!fMapping) {
199 AliWarning("Mapping array was not initalized correctly !");
200 return -1;
201 }
202 if (hwAddress > fMaxHWAddress) {
203 AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
204 return -1;
205 }
206 Int_t row = fMapping[hwAddress][0];
207 if (row == -1)
208 AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
209
210 return row;
211}
212
213//_____________________________________________________________________________
214Int_t AliCaloAltroMapping::GetPad(Int_t hwAddress) const
215{
216 // Return the column index
217 // Note the difference w.r.t to the base class notation
218 if (!fMapping) {
219 AliWarning("Mapping array was not initalized correctly !");
220 return -1;
221 }
222 if (hwAddress > fMaxHWAddress) {
223 AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
224 return -1;
225 }
226 Int_t col = fMapping[hwAddress][1];
227 if (col == -1)
228 AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
229
230 return col;
231}
232
233//_____________________________________________________________________________
234Int_t AliCaloAltroMapping::GetSector(Int_t hwAddress) const
235{
236 // Return the gain factor (0/1)
237 // Note the difference w.r.t to the base class notation
238 if (!fMapping) {
239 AliWarning("Mapping array was not initalized correctly !");
240 return -1;
241 }
242 if (hwAddress > fMaxHWAddress) {
243 AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
244 return -1;
245 }
246 Int_t gain = fMapping[hwAddress][2];
247 if (gain == -1)
248 AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
249
250 return gain;
251}
252
253//_____________________________________________________________________________
254void AliCaloAltroMapping::DeleteMappingArrays()
255{
256 // Deletes the arrays which have been
257 // allocated during the reading of the
258 // mapping file
259 if (fMapping) {
260 for (Int_t i = 0; i <= fMaxHWAddress; i++) delete [] fMapping[i];
261 delete [] fMapping;
262 }
263
264 if (fInvMappingLow) {
265 for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++)
266 delete [] fInvMappingLow[i];
267 delete [] fInvMappingLow;
268 }
269
270 if (fInvMappingHigh) {
271 for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++)
272 delete [] fInvMappingHigh[i];
273 delete [] fInvMappingHigh;
274 }
275}