[u/mrichter/AliRoot.git] / RAW / AliCaloAltroMapping.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
/* $Id$ */

// This class handles the mapping of the Altro channels in the PHOS/EMCAL
// The mapping is read from an external mapping files
// Author: C.Cheshkov

/// Exported from PHOS to be used also by EMCAL
/// November 2006 Gustavo Conesa Balbastre

#include "AliCaloAltroMapping.h"
#include "AliLog.h"
#include <Riostream.h>
//#include <stdlib.h>


ClassImp(AliCaloAltroMapping)

//_____________________________________________________________________________
AliCaloAltroMapping::AliCaloAltroMapping(const char *mappingFile):
  AliAltroMapping(mappingFile),
  fMinRow(0),
  fMaxRow(0),
  fMinCol(0),
  fMaxCol(0),
  fMapping(NULL),
  fInvMappingLow(NULL),
  fInvMappingHigh(NULL)
{
  // Constructor
  ReadMapping();
  CloseMappingFile();
}

//_____________________________________________________________________________
AliCaloAltroMapping::~AliCaloAltroMapping()
{
  // destructor
  DeleteMappingArrays();
}

//_____________________________________________________________________________
AliCaloAltroMapping::AliCaloAltroMapping(const AliCaloAltroMapping& mapping):
  AliAltroMapping(mapping),
  fMinRow(mapping.fMinRow),
  fMaxRow(mapping.fMaxRow),
  fMinCol(mapping.fMinCol),
  fMaxCol(mapping.fMaxCol),
  fMapping(mapping.fMapping),
  fInvMappingLow(mapping.fInvMappingLow),
  fInvMappingHigh(mapping.fInvMappingHigh)
{
// Copy Constructor

  Fatal("AliCaloAltroMapping", "copy constructor not implemented");
}

//_____________________________________________________________________________
AliCaloAltroMapping& AliCaloAltroMapping::operator = (const AliCaloAltroMapping& /*mapping*/)
{
//Assigment operator

  Fatal("operator =", "assignment operator not implemented");
  return *this;
}

//_____________________________________________________________________________
Bool_t AliCaloAltroMapping::ReadMapping()
{
  // Initalizes the ALTRO mapping from a file
  // Look at the Calo module for the format of
  // the mapping file
  if (!fIn) {
    AliFatal("Mapping file has not been opened !");
    return kFALSE;
  }

  fMinRow = 0x7fffffff;
  fMaxRow = 0;
  fMinCol = 0x7fffffff;
  fMaxCol = 0;
  fMapping = new Short_t*[fMaxHWAddress+1];
  for (Int_t i = 0; i <= fMaxHWAddress; i++) {
    fMapping[i] = new Short_t[3];
    fMapping[i][0] = fMapping[i][1] = fMapping[i][2] = -1;
  }
 
  for(Int_t i = 0; i < fNumberOfChannels ; i++) { // 1792 = 2*896 channels connected to each RCU
    Int_t hwAddress;
    if (!(*fIn >> hwAddress)) {
      AliFatal("Syntax of the mapping file is wrong !");
      return kFALSE;
    }
    if (hwAddress > fMaxHWAddress) {
      AliFatal(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
      return kFALSE;
    }
    Int_t row,col,gain;
    if (!(*fIn >> row >> col >> gain)) {
      AliFatal("Syntax of the mapping file is wrong !");
      return kFALSE;
    }

    if (gain < 0 || gain > 1) {
      AliFatal(Form("Wrong gain value found (%d)! Should be 0 or 1 !",gain));
      return kFALSE;
    }
 
    fMapping[hwAddress][0] = row;
    fMapping[hwAddress][1] = col;
    fMapping[hwAddress][2] = gain;

    if (row > fMaxRow) fMaxRow = row;
    if (row < fMinRow) fMinRow = row;
    if (col > fMaxCol) fMaxCol = col;
    if (col < fMinCol) fMinCol = col;

  }

  fInvMappingLow  = new Short_t*[fMaxRow - fMinRow + 1];
  fInvMappingHigh = new Short_t*[fMaxRow - fMinRow + 1];
  for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++) {
    fInvMappingLow[i]  = new Short_t[fMaxCol - fMinCol + 1];
    fInvMappingHigh[i] = new Short_t[fMaxCol - fMinCol + 1];
    for (Int_t j = 0; j <= (fMaxCol - fMinCol); j++) {
      fInvMappingLow[i][j]  = -1;
      fInvMappingHigh[i][j] = -1;
    }
  }

  for(Int_t i = 0; i <= fMaxHWAddress; i++) {
    Int_t row = fMapping[i][0];
    Int_t col = fMapping[i][1];
    Int_t gain = fMapping[i][2];
    if(row != -1 && col != -1) {
      if (gain == 0)
	fInvMappingLow[row-fMinRow][col-fMinCol] = i;
      if (gain == 1)
	fInvMappingHigh[row-fMinRow][col-fMinCol] = i;
    }
  }

  return kTRUE;
}

//_____________________________________________________________________________
Int_t AliCaloAltroMapping::GetHWAddress(Int_t row, Int_t col, Int_t gain) const
{
  // Get the content of the mapping array
  // return -1 in case there is no hardware
  // adress defined for these row-column-gain
  if (!fInvMappingLow || !fInvMappingHigh) {
    AliWarning("Mapping array was not initalized correctly !");
    return -1;
  }
  if (row < fMinRow || row > fMaxRow) {
    AliWarning(Form("Index of row (%d) outside the range (%d -> %d) !",row,fMinRow,fMaxRow));
    return -1;
  }
  if (col < fMinCol || col > fMaxCol) {
    AliWarning(Form("Index of column (%d) outside the range (0 -> %d) !",col,fMinCol,fMaxCol));
    return -1;
  }
  if (gain < 0 || gain > 1) {
    AliWarning(Form("Invalid gain (%d)! Should be 0 or 1 !",gain));
    return -1;
  }
  Int_t hwAddress = -1;
  if (gain == 0)
    hwAddress = fInvMappingLow[row-fMinRow][col-fMinCol];
  if (gain == 1)
    hwAddress = fInvMappingHigh[row-fMinRow][col-fMinCol];

  if (hwAddress == -1)
    AliWarning(Form("Hardware (ALTRO) adress is not defined for these row (%d), column (%d) and gain (%d) !",row,col,gain));

  return hwAddress;
}

//_____________________________________________________________________________
Int_t AliCaloAltroMapping::GetPadRow(Int_t hwAddress) const
{
  // Return the row index
  // Note the difference w.r.t to the base class notation
  if (!fMapping) {
    AliWarning("Mapping array was not initalized correctly !");
    return -1;
  }
  if (hwAddress > fMaxHWAddress) {
    AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
    return -1;
  }
  Int_t row = fMapping[hwAddress][0];
  if (row == -1)
    AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));

  return row;
}

//_____________________________________________________________________________
Int_t AliCaloAltroMapping::GetPad(Int_t hwAddress) const
{
  // Return the column index
  // Note the difference w.r.t to the base class notation
  if (!fMapping) {
    AliWarning("Mapping array was not initalized correctly !");
    return -1;
  }
  if (hwAddress > fMaxHWAddress) {
    AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
    return -1;
  }
  Int_t col = fMapping[hwAddress][1];
  if (col == -1)
    AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));

  return col;
}

//_____________________________________________________________________________
Int_t AliCaloAltroMapping::GetSector(Int_t hwAddress) const
{
  // Return the gain factor (0/1)
  // Note the difference w.r.t to the base class notation
  if (!fMapping) {
    AliWarning("Mapping array was not initalized correctly !");
    return -1;
  }
  if (hwAddress > fMaxHWAddress) {
    AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
    return -1;
  }
  Int_t gain = fMapping[hwAddress][2];
  if (gain == -1)
    AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));

  return gain;
}

//_____________________________________________________________________________
void AliCaloAltroMapping::DeleteMappingArrays()
{
  // Deletes the arrays which have been
  // allocated during the reading of the
  // mapping file
  if (fMapping) {
    for (Int_t i = 0; i <= fMaxHWAddress; i++) delete [] fMapping[i];
    delete [] fMapping;
  }

  if (fInvMappingLow) {
    for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++)
      delete [] fInvMappingLow[i];
    delete [] fInvMappingLow;
  }

  if (fInvMappingHigh) {
    for (Int_t i = 0; i <= (fMaxRow - fMinRow); i++)
      delete [] fInvMappingHigh[i];
    delete [] fInvMappingHigh;
  }
}
Commit	Line	Data
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
	30	ClassImp(AliCaloAltroMapping)
	31
	32	//_____________________________________________________________________________
	33	AliCaloAltroMapping::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	//_____________________________________________________________________________
	49	AliCaloAltroMapping::~AliCaloAltroMapping()
	50	{
	51	// destructor
	52	DeleteMappingArrays();
	53	}
	54
	55	//_____________________________________________________________________________
	56	AliCaloAltroMapping::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	//_____________________________________________________________________________
72	AliCaloAltroMapping& AliCaloAltroMapping::operator = (const AliCaloAltroMapping& /mapping/)
73	{
74	//Assigment operator
75
76	Fatal("operator =", "assignment operator not implemented");
77	return *this;
78	}
79
80	//_____________________________________________________________________________
81	Bool_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	//_____________________________________________________________________________
160	Int_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	//_____________________________________________________________________________
194	Int_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	//_____________________________________________________________________________
214	Int_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	//_____________________________________________________________________________
234	Int_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	//_____________________________________________________________________________
254	void 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	}