RAW/AliCaloAltroMapping.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 /* $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():
  34   AliAltroMapping(),
  35   fMinRow(0),
  36   fMaxRow(0),
  37   fMinCol(0),
  38   fMaxCol(0),
  39   fInvMappingLow(NULL),
  40   fInvMappingHigh(NULL)
  41 {
  42   // Default constructor
  43 }
  44
  45 //_____________________________________________________________________________
  46 AliCaloAltroMapping::AliCaloAltroMapping(const char *mappingFile):
  47   AliAltroMapping(mappingFile),
  48   fMinRow(0),
  49   fMaxRow(0),
  50   fMinCol(0),
  51   fMaxCol(0),
  52   fInvMappingLow(NULL),
  53   fInvMappingHigh(NULL)
  54 {
  55   // Constructor
  56   ReadMapping();
  57   CloseMappingFile();
  58 }
  59
  60 //_____________________________________________________________________________
  61 AliCaloAltroMapping::~AliCaloAltroMapping()
  62 {
  63   // destructor
  64   // Deletes the arrays which have been
  65   // allocated during the reading of the
  66   // mapping file
  67   if (fInvMappingLow) delete [] fInvMappingLow;
  68
  69   if (fInvMappingHigh) delete [] fInvMappingHigh;
  70 }
  71
  72 //_____________________________________________________________________________
  73 Bool_t AliCaloAltroMapping::ReadMapping()
  74 {
  75   // Initalizes the ALTRO mapping from a file
  76   // Look at the Calo module for the format of
  77   // the mapping file
  78   if (!fIn) {
  79     AliFatal("Mapping file has not been opened !");
  80     return kFALSE;
  81   }
  82
  83   fMinRow = 0x7fffffff;
  84   fMaxRow = 0;
  85   fMinCol = 0x7fffffff;
  86   fMaxCol = 0;
  87   fMappingSize = 3*(fMaxHWAddress+1);
  88   fMapping = new Short_t[fMappingSize];
  89   for (Int_t i = 0; i <= fMaxHWAddress; i++) {
  90     fMapping[3*i] = fMapping[3*i+1] = fMapping[3*i+2] = -1;
  91   }
  92
  93   for(Int_t i = 0; i < fNumberOfChannels ; i++) { // 1792 = 2*896 channels connected to each RCU
  94     Int_t hwAddress;
  95     if (!(*fIn >> hwAddress)) {
  96       AliFatal("Syntax of the mapping file is wrong !");
  97       return kFALSE;
  98     }
  99     if (hwAddress > fMaxHWAddress) {
 100       AliFatal(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
 101       return kFALSE;
 102     }
 103     Int_t row,col,caloFlag;
 104     if (!(*fIn >> row >> col >> caloFlag)) {
 105       AliFatal("Syntax of the mapping file is wrong !");
 106       return kFALSE;
 107     }
 108
 109     if (caloFlag < 0 || caloFlag > 3) {
 110       AliFatal(Form("Wrong CaloFlag value found (%d)! Should be 0 ,1, 2 or 3 !",caloFlag));
 111       return kFALSE;
 112     }
 113
 114     fMapping[3*hwAddress] = row;
 115     fMapping[3*hwAddress+1] = col;
 116     fMapping[3*hwAddress+2] = caloFlag;
 117
 118     if (row > fMaxRow) fMaxRow = row;
 119     if (row < fMinRow) fMinRow = row;
 120     if (col > fMaxCol) fMaxCol = col;
 121     if (col < fMinCol) fMinCol = col;
 122
 123   }
 124
 125   return kTRUE;
 126 }
 127
 128 //_____________________________________________________________________________
 129 Bool_t AliCaloAltroMapping::CreateInvMapping()
 130 {
 131   // Create the inverse mapping
 132   // needed for the simulation of
 133   // raw data
 134   if (fInvMappingLow) return kTRUE;
 135
 136   if (!fMapping) {
 137     AliWarning("Mapping array was not initalized correctly ! Impossible to create the inverse mapping !");
 138     return kFALSE;
 139   }
 140
 141   Int_t nRows = fMaxRow - fMinRow + 1;
 142   Int_t nCols = fMaxCol - fMinCol + 1;
 143   Int_t invMappingSize = nRows*nCols;
 144
 145   fInvMappingLow  = new Short_t[invMappingSize];
 146   fInvMappingHigh = new Short_t[invMappingSize];
 147   for (Int_t i = 0; i < nRows; i++) {
 148     for (Int_t j = 0; j < nCols; j++) {
 149       fInvMappingLow[nCols*i+j]  = -1;
 150       fInvMappingHigh[nCols*i+j] = -1;
 151     }
 152   }
 153   for(Int_t i = 0; i <= fMaxHWAddress; i++) {
 154     Int_t row = fMapping[3*i];
 155     Int_t col = fMapping[3*i+1];
 156     Int_t caloFlag = fMapping[3*i+2];
 157     if(row != -1 && col != -1) {
 158       if (caloFlag == 0)
 159         fInvMappingLow[nCols*(row-fMinRow)+(col-fMinCol)] = i;
 160       if (caloFlag == 1)
 161         fInvMappingHigh[nCols*(row-fMinRow)+(col-fMinCol)] = i;
 162     }
 163   }
 164
 165   return kTRUE;
 166 }
 167
 168 //_____________________________________________________________________________
 169 Int_t AliCaloAltroMapping::GetHWAddress(Int_t row, Int_t col, Int_t caloFlag)
 170 {
 171   // Get the content of the mapping array
 172   // return -1 in case there is no hardware
 173   // adress defined for these row-column-caloFlag
 174   if (!fInvMappingLow || !fInvMappingHigh) {
 175     if (!CreateInvMapping()) return -1;
 176   }
 177   if (row < fMinRow || row > fMaxRow) {
 178     AliWarning(Form("Index of row (%d) outside the range (%d -> %d) !",row,fMinRow,fMaxRow));
 179     return -1;
 180   }
 181   if (col < fMinCol || col > fMaxCol) {
 182     AliWarning(Form("Index of column (%d) outside the range (0 -> %d) !",col,fMinCol,fMaxCol));
 183     return -1;
 184   }
 185   if (caloFlag < 0 || caloFlag > 3) {
 186     AliWarning(Form("Invalid caloFlag (%d)! Should be 0, 1, 2 or 3 !",caloFlag));
 187     return -1;
 188   }
 189   Int_t hwAddress = -1;
 190   if (caloFlag == 0)
 191     hwAddress = fInvMappingLow[(fMaxCol - fMinCol + 1)*(row-fMinRow)+(col-fMinCol)];
 192   if (caloFlag == 1)
 193     hwAddress = fInvMappingHigh[(fMaxCol - fMinCol + 1)*(row-fMinRow)+(col-fMinCol)];
 194
 195   if (hwAddress == -1)
 196     AliWarning(Form("Hardware (ALTRO) adress is not defined for these row (%d), column (%d) and caloFlag (%d) !",row,col,caloFlag));
 197
 198   return hwAddress;
 199 }
 200
 201 //_____________________________________________________________________________
 202 Int_t AliCaloAltroMapping::GetPadRow(Int_t hwAddress) const
 203 {
 204   // Return the row index
 205   // Note the difference w.r.t to the base class notation
 206   if (!fMapping) {
 207     AliWarning("Mapping array was not initalized correctly !");
 208     return -1;
 209   }
 210   if (hwAddress > fMaxHWAddress) {
 211     AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
 212     return -1;
 213   }
 214   Int_t row = fMapping[3*hwAddress];
 215   if (row == -1)
 216     AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
 217
 218   return row;
 219 }
 220
 221 //_____________________________________________________________________________
 222 Int_t AliCaloAltroMapping::GetPad(Int_t hwAddress) const
 223 {
 224   // Return the column index
 225   // Note the difference w.r.t to the base class notation
 226   if (!fMapping) {
 227     AliWarning("Mapping array was not initalized correctly !");
 228     return -1;
 229   }
 230   if (hwAddress > fMaxHWAddress) {
 231     AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
 232     return -1;
 233   }
 234   Int_t col = fMapping[3*hwAddress+1];
 235   if (col == -1)
 236     AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
 237
 238   return col;
 239 }
 240
 241 //_____________________________________________________________________________
 242 Int_t AliCaloAltroMapping::GetSector(Int_t hwAddress) const
 243 {
 244   // Return the caloFlag factor (0/1)
 245   // Note the difference w.r.t to the base class notation
 246   if (!fMapping) {
 247     AliWarning("Mapping array was not initalized correctly !");
 248     return -1;
 249   }
 250   if (hwAddress > fMaxHWAddress) {
 251     AliWarning(Form("Hardware (ALTRO) adress (%d) outside the range (0 -> %d) !",hwAddress,fMaxHWAddress));
 252     return -1;
 253   }
 254   Int_t caloFlag = fMapping[3*hwAddress+2];
 255   if (caloFlag == -1)
 256     AliWarning(Form("Hardware (ALTRO) adress (%d) is not defined !",hwAddress));
 257
 258   return caloFlag;
 259 }