Increasing histo clu. lay.1 upper lim.
[u/mrichter/AliRoot.git] / 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 }