bug fixed
[u/mrichter/AliRoot.git] / RAW / AliCaloRawStream.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
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
42ClassImp(AliCaloRawStream)
43
44
45//_____________________________________________________________________________
467957dc 46AliCaloRawStream::AliCaloRawStream(AliRawReader* rawReader, TString calo, AliAltroMapping **mapping) :
d84933b0 47 AliAltroRawStream(rawReader),
48 fModule(-1),
49 fPrevModule(-1),
50 fRow(-1),
51 fPrevRow(-1),
52 fColumn(-1),
53 fPrevColumn(-1),
c7ecf7e8 54 fCaloFlag(0),
55 fFilter(0),
467957dc 56 fNRCU(0),
cee89af9 57 fNSides(0),
58 fCalo(calo),
467957dc 59 fExternalMapping(kFALSE)
d84933b0 60{
cee89af9 61 // create an object to read PHOS/EMCAL raw digits
d84933b0 62 SelectRawData(calo);
63
64 // PHOS and EMCAL have differen number of RCU per module
cee89af9 65 //For PHOS
d84933b0 66 fNRCU = 4;
cee89af9 67 fNSides = 1;
68 //For EMCAL
69 TString sides[]={"A","C"};
70 if(fCalo == "EMCAL") {
71 fNRCU = 2;
72 fNSides = 2;
73 }
d84933b0 74
467957dc 75 if (mapping == NULL) {
76 TString path = gSystem->Getenv("ALICE_ROOT");
cee89af9 77 path += "/"+fCalo+"/mapping/RCU";
467957dc 78 TString path2;
cee89af9 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());
cf19d394 86 fMapping[j*fNRCU+ i] = new AliCaloAltroMapping(path2.Data());
cee89af9 87 }
467957dc 88 }
89 }
90 else {
91 fExternalMapping = kTRUE;
cee89af9 92 //printf("AliCaloRawStream::External mapping N: RCU %d, sides %d \n", fNRCU,fNSides);
93 for(Int_t i = 0; i < fNRCU*fNSides; i++)
467957dc 94 fMapping[i] = mapping[i];
cee89af9 95
d84933b0 96 }
d84933b0 97}
98
99//_____________________________________________________________________________
100AliCaloRawStream::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),
c7ecf7e8 108 fCaloFlag(0),
109 fFilter(0),
467957dc 110 fNRCU(0),
cee89af9 111 fNSides(0),
112 fCalo(""),
467957dc 113 fExternalMapping(kFALSE)
d84933b0 114{
115 Fatal("AliCaloRawStream", "copy constructor not implemented");
116}
117
118//_____________________________________________________________________________
119AliCaloRawStream& AliCaloRawStream::operator = (const AliCaloRawStream&
120 /* stream */)
121{
122 Fatal("operator =", "assignment operator not implemented");
123 return *this;
124}
125
126//_____________________________________________________________________________
127AliCaloRawStream::~AliCaloRawStream()
128{
129// destructor
130
467957dc 131 if (!fExternalMapping)
cee89af9 132 for(Int_t i = 0; i < fNRCU*fNSides; i++)
467957dc 133 delete fMapping[i];
d84933b0 134}
135
136//_____________________________________________________________________________
137void AliCaloRawStream::Reset()
138{
139 // reset phos/emcal raw stream params
140 AliAltroRawStream::Reset();
141 fModule = fPrevModule = fRow = fPrevRow = fColumn = fPrevColumn = -1;
c7ecf7e8 142 fFilter = fCaloFlag = 0;
cee89af9 143 fCalo="";
d84933b0 144}
145
146//_____________________________________________________________________________
147Bool_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()) {
c7ecf7e8 156 if (IsNewHWAddress()) {
d84933b0 157 ApplyAltroMapping();
c7ecf7e8 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 }
d84933b0 165 return kTRUE;
166 }
167 else
168 return kFALSE;
169}
170
171//_____________________________________________________________________________
172void 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
cee89af9 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
d84933b0 186 Short_t hwAddress = GetHWAddress();
187 fRow = fMapping[rcuIndex]->GetPadRow(hwAddress);
188 fColumn = fMapping[rcuIndex]->GetPad(hwAddress);
c7ecf7e8 189 fCaloFlag = fMapping[rcuIndex]->GetSector(hwAddress);
d84933b0 190
191}