Added the most important Calo files, changed dependent files. Compiles. Calo library...
[u/mrichter/AliRoot.git] / HLT / CALO / AliHLTCaloMapper.cxx
1 // $Id: AliHLTCalorimeterMapper.cxx 34622 2009-09-04 13:22:01Z odjuvsla $
2
3 /**************************************************************************
4  * This file is property of and copyright by the Experimental Nuclear     *
5  * Physics Group, Dep. of Physics                                         *
6  * University of Oslo, Norway, 2006                                       *
7  *                                                                        * 
8  * Author: Per Thomas Hille perthi@fys.uio.no for the ALICE DCS Project.  *
9  * Contributors are mentioned in the code where appropriate.              *
10  * Please report bugs to perthi@fys.uio.no                                * 
11  *                                                                        *
12  * Permission to use, copy, modify and distribute this software and its   *
13  * documentation strictly for non-commercial purposes is hereby granted   *
14  * without fee, provided that the above copyright notice appears in all   *
15  * copies and that both the copyright notice and this permission notice   *
16  * appear in the supporting documentation. The authors make no claims     *
17  * about the suitability of this software for any purpose. It is          *
18  * provided "as is" without express or implied warranty.                  *
19  **************************************************************************/
20
21 //
22 // Mapping class fro mapping
23 // from hardware address to geometrical address
24 //
25 //
26
27 #include "AliHLTCaloMapper.h"
28 #include "AliHLTLogging.h"
29 #include "Rtypes.h"
30 #include "unistd.h"
31 #include <iostream>
32 #include "AliHLTCaloCoordinate.h"
33 #include "AliHLTCaloConstantsHandler.h"
34 #include "AliHLTCaloConstants.h"
35
36 //ClassImp(AliHLTCaloMapper)
37
38 AliHLTCaloMapper::AliHLTCaloMapper( const unsigned long  specification , TString det) :  
39   AliHLTLogging(), 
40   AliHLTCaloConstantsHandler(det),
41   fHw2geomapPtr(0),
42   fCellSize(0),
43   fSpecification(specification),
44   fIsInitializedMapping(false),
45   fSpecificationMapPtr(0)
46 {  
47   //see header file for class documentation
48 }
49
50
51 AliHLTCaloMapper::~AliHLTCaloMapper()
52 {
53   delete []  fHw2geomapPtr;
54   fHw2geomapPtr = 0;
55 }
56
57
58 //void AliHLTCaloMapper::InitAltroMapping(){
59   //Virtual base class
60 //}
61
62 //void AliHLTCaloMapper::InitDDLSpecificationMapping() {
63   //Virtual base class
64 //}
65
66
67
68 bool 
69 AliHLTCaloMapper::GetIsInitializedMapping()
70 {
71   return  fIsInitializedMapping;
72 }
73
74
75 char* 
76 AliHLTCaloMapper::GetFilePath()
77 {
78   return  fFilepath;
79 }
80
81
82 int 
83 AliHLTCaloMapper::GetChannelID(const AliHLTUInt32_t spec, const Int_t hadd)
84 {
85   Short_t index = GetDDLFromSpec(spec);
86   if( index < 0 )
87     {
88       //      HLTError("Specification 0x%X not consistent with single DDL in PHOS", spec);
89       return index;
90     }
91   else
92     {
93       return ((fHw2geomapPtr[hadd].fXCol   ) |
94               ((fHw2geomapPtr[hadd].fZRow  ) << 6) |
95               (fHw2geomapPtr[hadd].fGain << 12) |
96               fSpecificationMapPtr[index].fModId << 13);
97     }
98 }
99
100
101
102 // UShort_t
103 // AliHLTCaloMapper::GetChannelID(Int_t specification, Int_t hwAddress)
104 // {
105   
106 //   Short_t index = 0;
107
108 //   if(specification == 0x00001) index = 0;
109 //   else if(specification == 0x00002) index = 1;
110 //   else if(specification == 0x00004) index = 2;
111 //   else if(specification == 0x00008) index = 3;
112
113 //   else if(specification == 0x00010) index = 4;
114 //   else if(specification == 0x00020) index = 5;
115 //   else if(specification == 0x00040) index = 6;
116 //   else if(specification == 0x00080) index = 7;
117
118 //   else if(specification == 0x00100) index = 8;
119 //   else if(specification == 0x00200) index = 9;
120 //   else if(specification == 0x00400) index = 10;
121 //   else if(specification == 0x00800) index = 11;
122
123 //   else if(specification == 0x01000) index = 12;
124 //   else if(specification == 0x02000) index = 13;
125 //   else if(specification == 0x04000) index = 14;
126 //   else if(specification == 0x08000) index = 15;
127
128 //   else if(specification == 0x10000) index = 16;
129 //   else if(specification == 0x20000) index = 17;
130 //   else if(specification == 0x40000) index = 18;
131 //   else if(specification == 0x80000) index = 19;
132
133 //   else HLTError("Specification 0x%X not consistent with single DDL in PHOS", specification);
134 //   //  HLTError("Channel ID: 0x%X Coordinates: x = %d, z = %d, gain = %d", ((fHw2geomapPtr[hwAddress].fXCol + fSpecificationMapPtr[index].fRcuXOffset) |((fHw2geomapPtr[hwAddress].fZRow + fSpecificationMapPtr[index].fRcuZOffset) << 6) | (fHw2geomapPtr[hwAddress].fGain << 12) | fSpecificationMapPtr[index].fModId << 13),
135 //   //    fHw2geomapPtr[hwAddress].fXCol,
136 //   //    fHw2geomapPtr[hwAddress].fZRow, 
137 //   //    fHw2geomapPtr[hwAddress].fGain);
138 //   return ((fHw2geomapPtr[hwAddress].fXCol + fSpecificationMapPtr[index].fRcuXOffset) |
139 //        ((fHw2geomapPtr[hwAddress].fZRow + fSpecificationMapPtr[index].fRcuZOffset) << 6) |
140 //        (fHw2geomapPtr[hwAddress].fGain << 12) |
141 //        fSpecificationMapPtr[index].fModId << 13);
142 // }
143
144
145 void
146 AliHLTCaloMapper::GetChannelCoord(const UShort_t channelId, UShort_t* channelCoord)
147 {
148   channelCoord[0] = channelId&0x3f;
149   channelCoord[1] = (channelId >> 6)&0x3f;
150   channelCoord[2] = (channelId >> 12)&0x1;
151   channelCoord[3] = (channelId >> 13)&0x1f;
152   //  printf("Channel ID: 0x%X Coordinates: x = %d, z = %d, gain = %d\n", channelId, channelCoord[0], channelCoord[1], channelCoord[2]);
153 }
154
155 void
156 AliHLTCaloMapper::ChannelId2Coordinate(const int channelId,    AliHLTCaloCoordinate &channelCoord)
157 {
158   channelCoord.fX = channelId&0x3f;
159   channelCoord.fZ = (channelId >> 6)&0x3f;
160   channelCoord.fGain = (channelId >> 12)&0x1;
161   channelCoord.fModuleId  = (channelId >> 13)&0x1f;
162   //  printf("Channel ID: 0x%X Coordinates: x = %d, z = %d, gain = %d\n", channelId, channelCoord[0], channelCoord[1], channelCoord[2]);
163 }
164
165
166 void
167 AliHLTCaloMapper::GetLocalCoord(const int channelId, Float_t* localCoord) const
168 {
169   localCoord[0] = (static_cast<Float_t>(channelId&0x3f) - fCaloConstants->GetNXCOLUMNSMOD()/2)* fCaloConstants->GetCELLSTEP();
170   localCoord[1] = (static_cast<Float_t>((channelId >> 6)&0x3f) - fCaloConstants->GetNZROWSMOD()/2) * fCaloConstants->GetCELLSTEP();
171   //  printf("Local coordinates: x = %f, z = %f\n", channelCoord[0], channelCoord[1]);
172 }
173
174
175 int  
176 AliHLTCaloMapper::GetDDLFromSpec( const AliHLTUInt32_t spec )
177 {
178   int tmpIndex = -1;
179   for(int i=0; i < 32; i++ )
180     {
181       if (spec >> i ==1)
182         {
183           tmpIndex = i;
184           break;
185         }
186     }
187    if(  tmpIndex  < 0)
188     {
189       //   HLTError("Specification %d, not consistent with any DDL in PHOS or EMCAL", spec  );
190     }
191   return tmpIndex;
192 }
193
194
195 Int_t 
196 AliHLTCaloMapper::GetModuleFromSpec(Int_t specification)
197 {
198   Int_t module = -1;
199       
200   if(specification & 0xf) module = 0;
201   else if((specification >> 4) & 0xf) module = 1;
202   else if((specification >> 8) & 0xf) module = 2;
203   else if((specification >> 12) & 0xf) module = 3;
204   else if((specification >> 16) & 0xf) module = 4;
205  
206   else HLTError("Specification 0x%X not consistent with single module in PHOS", specification);
207
208   return module;
209 }
210
211
212 unsigned long 
213 AliHLTCaloMapper::GetSpecFromDDLIndex( const int ddlindex )
214 {
215   int iret = (unsigned long)1  <<  ddlindex;  
216
217   //  return  ((unsigned long)1)  <<  ddlindex ) ;
218
219   return iret;
220
221 }
222