RAW/AliCaloRawStream.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 ///
  20 /// This class provides access to PHOS/EMCAL digits in raw data.
  21 ///
  22 /// It loops over all PHOS/EMCAL digits in the raw data given by the AliRawReader.
  23 /// The Next method goes to the next digit. If there are no digits left
  24 /// it returns kFALSE.
  25 /// Several getters provide information about the current digit.
  26 /// usage:
  27 /// root > AliRawReaderFile rawReader ;
  28 /// root > AliCaloRawStream input(&rawReader) ;
  29 /// root > while (input.Next()) .....
  30 ///
  31 ///Modification: Class exported from PHOS to be used by EMCAL and PHOS
  32 ///November 2006 Gustavo Conesa Balbastre
  33 ///////////////////////////////////////////////////////////////////////////////
  34
  35 #include <TString.h>
  36 #include <TSystem.h>
  37
  38 #include "AliCaloRawStream.h"
  39 #include "AliRawReader.h"
  40 #include "AliCaloAltroMapping.h"
  41
  42 ClassImp(AliCaloRawStream)
  43
  44
  45 //_____________________________________________________________________________
  46 AliCaloRawStream::AliCaloRawStream(AliRawReader* rawReader, TString calo, AliAltroMapping **mapping) :
  47   AliAltroRawStream(rawReader),
  48   fModule(-1),
  49   fPrevModule(-1),
  50   fRow(-1),
  51   fPrevRow(-1),
  52   fColumn(-1),
  53   fPrevColumn(-1),
  54   fCaloFlag(0),
  55   fFilter(0),
  56   fNRCU(0),
  57   fNSides(0),
  58   fCalo(calo),
  59   fExternalMapping(kFALSE)
  60 {
  61   // create an object to read PHOS/EMCAL raw digits
  62   SelectRawData(calo);
  63
  64   // PHOS and EMCAL have differen number of RCU per module
  65   //For PHOS
  66   fNRCU = 4;
  67   fNSides = 1;
  68   //For EMCAL
  69   TString sides[]={"A","C"};
  70   if(fCalo == "EMCAL")  {
  71     fNRCU = 2;
  72     fNSides = 2;
  73   }
  74
  75   if (mapping == NULL) {
  76     TString path = gSystem->Getenv("ALICE_ROOT");
  77     path += "/"+fCalo+"/mapping/RCU";
  78     TString path2;
  79     for(Int_t j = 0; j < fNSides; j++){
  80       for(Int_t i = 0; i < fNRCU; i++) {
  81         path2 = path;
  82         path2 += i;
  83         if(fCalo == "EMCAL") path2 += sides[j];
  84         path2 += ".data";
  85         //printf("AliCaloRawStream::RCU:  %s\n",path2.Data());
  86         fMapping[j*fNRCU+ i] = new AliCaloAltroMapping(path2.Data());
  87       }
  88     }
  89   }
  90   else {
  91     fExternalMapping = kTRUE;
  92     //printf("AliCaloRawStream::External mapping N: RCU %d, sides %d \n",  fNRCU,fNSides);
  93     for(Int_t i = 0; i < fNRCU*fNSides; i++)
  94       fMapping[i] = mapping[i];
  95
  96   }
  97 }
  98
  99 //_____________________________________________________________________________
 100 AliCaloRawStream::AliCaloRawStream(const AliCaloRawStream& stream) :
 101   AliAltroRawStream(stream),
 102   fModule(-1),
 103   fPrevModule(-1),
 104   fRow(-1),
 105   fPrevRow(-1),
 106   fColumn(-1),
 107   fPrevColumn(-1),
 108   fCaloFlag(0),
 109   fFilter(0),
 110   fNRCU(0),
 111   fNSides(0),
 112   fCalo(""),
 113   fExternalMapping(kFALSE)
 114 {
 115   Fatal("AliCaloRawStream", "copy constructor not implemented");
 116 }
 117
 118 //_____________________________________________________________________________
 119 AliCaloRawStream& AliCaloRawStream::operator = (const AliCaloRawStream&
 120                                               /* stream */)
 121 {
 122   Fatal("operator =", "assignment operator not implemented");
 123   return *this;
 124 }
 125
 126 //_____________________________________________________________________________
 127 AliCaloRawStream::~AliCaloRawStream()
 128 {
 129 // destructor
 130
 131   if (!fExternalMapping)
 132     for(Int_t i = 0; i < fNRCU*fNSides; i++)
 133       delete fMapping[i];
 134 }
 135
 136 //_____________________________________________________________________________
 137 void AliCaloRawStream::Reset()
 138 {
 139   // reset phos/emcal raw stream params
 140   AliAltroRawStream::Reset();
 141   fModule = fPrevModule = fRow = fPrevRow = fColumn = fPrevColumn = -1;
 142   fFilter = fCaloFlag = 0;
 143   fCalo="";
 144 }
 145
 146 //_____________________________________________________________________________
 147 Bool_t AliCaloRawStream::Next()
 148 {
 149   // Read next PHOS/EMCAL signal
 150   // Apply the PHOS/EMCAL altro mapping to get
 151   // the module,row and column indeces
 152   fPrevModule = fModule;
 153   fPrevRow = fRow;
 154   fPrevColumn = fColumn;
 155   if (AliAltroRawStream::Next()) {
 156     if (IsNewHWAddress()) {
 157       ApplyAltroMapping();
 158       if ( fFilter > 0 ) { // some data should be filtered out
 159         if ( (fFilter & (1<<fCaloFlag)) != 0) {
 160           // this particular data should be filtered out
 161           Next(); // go to the next address instead
 162         }
 163       }
 164     }
 165     return kTRUE;
 166   }
 167   else
 168     return kFALSE;
 169 }
 170
 171 //_____________________________________________________________________________
 172 void AliCaloRawStream::ApplyAltroMapping()
 173 {
 174   // Take the DDL index, load
 175   // the corresponding altro mapping
 176   // object and fill the sector,row and pad indeces
 177   Int_t ddlNumber = GetDDLNumber();
 178   fModule = ddlNumber / fNRCU;
 179
 180   Int_t rcuIndex = ddlNumber % fNRCU;
 181
 182   if(fCalo=="EMCAL"){ // EMCAL may need to increase RCU index for the maps
 183     if (fModule%2 == 1) { rcuIndex += 2; } // other='C' side maps
 184   }
 185
 186   Short_t hwAddress = GetHWAddress();
 187   fRow = fMapping[rcuIndex]->GetPadRow(hwAddress);
 188   fColumn = fMapping[rcuIndex]->GetPad(hwAddress);
 189   fCaloFlag = fMapping[rcuIndex]->GetSector(hwAddress);
 190
 191 }