add set and getter for neutral energy fraction
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALTriggerElectronics.cxx
CommitLineData
916f1e76 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/*
17
18
19EMCal trigger electronics manager L0/L1
20can handle both simulated digits and raw data
21Author: R. GUERNANE LPSC Grenoble CNRS/IN2P3
22*/
23
24#include "AliEMCALTriggerElectronics.h"
25#include "AliEMCALTriggerTRU.h"
26#include "AliEMCALTriggerSTU.h"
27#include "AliEMCALGeometry.h"
28#include "AliRunLoader.h"
29#include "AliEMCAL.h"
30#include "AliRun.h"
fff39dd1 31#include "AliEMCALTriggerDCSConfig.h"
916f1e76 32#include "AliEMCALTriggerData.h"
33#include "AliEMCALDigit.h"
34#include "AliCaloRawStreamV3.h"
916f1e76 35#include "AliEMCALTriggerSTURawStream.h"
36#include "AliEMCALDigit.h"
de39a0ff 37#include "AliEMCALTriggerRawDigit.h"
38#include "AliEMCALTriggerPatch.h"
916f1e76 39
40#include <TVector2.h>
41#include <TClonesArray.h>
916f1e76 42
43namespace
44{
8cf986c8 45 const Int_t kNTRU = 30;
916f1e76 46}
47
48ClassImp(AliEMCALTriggerElectronics)
49
50//__________________
fff39dd1 51AliEMCALTriggerElectronics::AliEMCALTriggerElectronics(const AliEMCALTriggerDCSConfig *dcsConf) : TObject(),
916f1e76 52fTRU(new TClonesArray("AliEMCALTriggerTRU",32)),
53fSTU(0x0)
54{
55 TVector2 rSize;
56
57 rSize.Set( 24., 4. );
58
59 // 32 TRUs
fff39dd1 60 for (Int_t i=0;i<kNTRU;i++)
61 {
62 AliEMCALTriggerTRUDCSConfig* truConf = dcsConf->GetTRUDCSConfig(i);
de39a0ff 63 new ((*fTRU)[i]) AliEMCALTriggerTRU(truConf, rSize, i % 2);
fff39dd1 64 }
65
916f1e76 66 rSize.Set( 48., 64. );
67
68 // 1 STU
fff39dd1 69 AliEMCALTriggerSTUDCSConfig* stuConf = dcsConf->GetSTUDCSConfig();
70 fSTU = new AliEMCALTriggerSTU(stuConf, rSize);
916f1e76 71
de39a0ff 72 TString str = "map";
73 for (Int_t i=0;i<kNTRU;i++) fSTU->Build(str,
74 i,
75 (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Map(),
76 (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize()
77 );
916f1e76 78}
79
80//________________
81AliEMCALTriggerElectronics::~AliEMCALTriggerElectronics()
82{
83 //
84 fTRU->Delete();
85 delete fSTU;
86}
87
88//__________________
de39a0ff 89void AliEMCALTriggerElectronics::Digits2Trigger(TClonesArray* digits, const Int_t V0M[], AliEMCALTriggerData* data)
916f1e76 90{
91 //
92 AliEMCALGeometry* geom = 0x0;
93
94 AliRunLoader *rl = AliRunLoader::Instance();
a51e676d 95 if (rl->GetAliRun() && rl->GetAliRun()->GetDetector("EMCAL")){
96 AliEMCAL* emcal = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"));
97 if(emcal)geom = emcal->GetGeometry();
98 }
99
100 if(!geom) geom = AliEMCALGeometry::GetInstance(AliEMCALGeometry::GetDefaultGeometryName());
101
102 if(!geom) AliError("Cannot access geometry!");
916f1e76 103
de39a0ff 104 // digits->Sort();
105
106 Int_t region[48][64], posMap[48][64];
107 for (Int_t i = 0; i < 48; i++) for (Int_t j = 0; j < 64; j++)
108 {
109 region[i][j] = 0;
110 posMap[i][j] = -1;
111 }
112
113 for (Int_t i = 0; i < digits->GetEntriesFast(); i++)
916f1e76 114 {
de39a0ff 115 AliEMCALTriggerRawDigit* digit = (AliEMCALTriggerRawDigit*)digits->At(i);
116
117 Int_t id = digit->GetId();
118
119 Int_t iTRU, iADC;
120
121 Bool_t isOK1 = geom->GetTRUFromAbsFastORIndex(id, iTRU, iADC);
122
c33e980f 123 if ((isOK1 && iTRU >= kNTRU) || !isOK1) continue;
124
de39a0ff 125 for (Int_t j = 0; j < digit->GetNSamples(); j++)
916f1e76 126 {
de39a0ff 127 Int_t time, amp;
128 Bool_t isOK2 = digit->GetTimeSample(j, time, amp);
916f1e76 129
de39a0ff 130 if (isOK1 && isOK2 && amp) (static_cast<AliEMCALTriggerTRU*>(fTRU->At(iTRU)))->SetADC(iADC, time, amp);
131 }
132
133 Int_t px, py;
134 if (geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py))
135 {
136 posMap[px][py] = i;
916f1e76 137
de39a0ff 138 if (fSTU->GetRawData() && digit->GetL1TimeSum() >= 0)
916f1e76 139 {
de39a0ff 140 region[px][py] = digit->GetL1TimeSum();
916f1e76 141 }
142 }
143 }
916f1e76 144
de39a0ff 145 Int_t iL0 = 0;
916f1e76 146
147 for (Int_t i=0; i<kNTRU; i++)
148 {
de39a0ff 149 AliDebug(999, Form("===========< TRU %2d >============\n", i));
916f1e76 150
de39a0ff 151 AliEMCALTriggerTRU* iTRU = static_cast<AliEMCALTriggerTRU*>(fTRU->At(i));
916f1e76 152
de39a0ff 153 // L0 is always computed from F-ALTRO
154 if (iTRU->L0())
916f1e76 155 {
de39a0ff 156 iL0 += iTRU->L0();
157
ddd131cc 158 Int_t sizeX = (Int_t) ((iTRU->PatchSize())->X() * (iTRU->SubRegionSize())->X());
de39a0ff 159
ddd131cc 160 Int_t sizeY = (Int_t) ((iTRU->PatchSize())->Y() * (iTRU->SubRegionSize())->Y());
de39a0ff 161
162 // transform local to global
163 TIter Next(&iTRU->Patches());
164 while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)Next())
165 {
166 Int_t px, py, id; p->Position(px, py);
167
168 if (geom->GetAbsFastORIndexFromPositionInTRU(i, px, py, id)
169 &&
170 geom->GetPositionInEMCALFromAbsFastORIndex(id, px, py)) p->SetPosition(px, py);
171
a61738e1 172 if (!data->GetMode()) // Simulation
de39a0ff 173 {
174 Int_t peaks = p->Peaks();
175
a61738e1 176 Int_t pos;
177 AliEMCALTriggerRawDigit* dig = 0x0;
178
de39a0ff 179 for (Int_t j = 0; j < sizeX * sizeY; j++)
180 {
181 if (peaks & (1 << j))
182 {
a61738e1 183 pos = posMap[px + j % sizeX][py + j / sizeX];
de39a0ff 184
185 if (pos == -1)
186 {
a61738e1 187 // Add a new digit
de39a0ff 188 new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
189
190 dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
191 }
192 else
193 {
194 dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
195 }
196
197 dig->SetL0Time(p->Time());
198 }
a61738e1 199 }
200
201 pos = posMap[px][py];
202
203 if (pos == -1)
204 {
205 // Add a new digit
206 new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
207
208 dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
209 }
210 else
211 {
212 dig = (AliEMCALTriggerRawDigit*)digits->At(pos);
213 }
214
215 dig->SetTriggerBit(kL0,0);
de39a0ff 216 }
217 }
218
219 data->SetL0Trigger(0, i, 1);
de39a0ff 220 }
221 else
222 data->SetL0Trigger(0, i, 0);
916f1e76 223 }
224
225 // A L0 has been issued, run L1
de39a0ff 226 // Depending on raw data enabled or not in STU data: L1 computation
227 // should be done from F-ALTRO or directly on TRU time sums in STU raw data
228 if (iL0)
916f1e76 229 {
de39a0ff 230 // Use L1 threshold from raw data when reconstructing raw data
231 if (data->GetMode())
232 {
233 fSTU->SetThreshold(kL1Gamma, data->GetL1GammaThreshold());
234 fSTU->SetThreshold(kL1Jet, data->GetL1JetThreshold() );
235 }
236 else
237 {
238 fSTU->ComputeThFromV0(V0M); // C/A
239 data->SetL1GammaThreshold(fSTU->GetThreshold(kL1Gamma));
240 data->SetL1JetThreshold( fSTU->GetThreshold(kL1Jet) );
241 }
916f1e76 242
de39a0ff 243 if (fSTU->GetRawData())
244 {
245 // Compute L1 from STU raw data
246 fSTU->SetRegion(region);
247 }
248 else
249 {
250 // Build STU raw data from F-ALTRO
251 TString str = "region";
252 for (Int_t i = 0; i < kNTRU; i++) fSTU->Build(str,
253 i,
254 (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->Region(),
255 (static_cast<AliEMCALTriggerTRU*>(fTRU->At(i)))->RegionSize());
256 }
916f1e76 257
de39a0ff 258 fSTU->L1(kL1Gamma);
916f1e76 259
a61738e1 260 Int_t id, px, py;
261 AliEMCALTriggerRawDigit* dig = 0x0;
262
263 TIterator* nP = 0x0;
264
265 nP = (fSTU->Patches()).MakeIterator();
266
267 while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next())
268 {
269 p->Position(px, py);
270
271 if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
272 {
273 if (posMap[px][py] == -1)
274 {
275 // Add a new digit
276 new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
277
278 dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
279 }
280 else
281 {
282 dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);
283 }
284
285 dig->SetTriggerBit(kL1Gamma,0);
286 }
287 }
916f1e76 288
289 fSTU->Reset();
290
de39a0ff 291 fSTU->L1(kL1Jet);
916f1e76 292
a61738e1 293 nP = (fSTU->Patches()).MakeIterator();
294
295 while (AliEMCALTriggerPatch* p = (AliEMCALTriggerPatch*)nP->Next())
296 {
297 p->Position(px, py);
298
299 px *= (Int_t)((fSTU->SubRegionSize())->X());
300
301 py *= (Int_t)((fSTU->SubRegionSize())->Y());
302
303 if (geom->GetAbsFastORIndexFromPositionInEMCAL(px, py, id))
304 {
305 if (posMap[px][py] == -1)
306 {
307 // Add a new digit
308 new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
309
310 dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
311 }
312 else
313 {
314 dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[px][py]);
315 }
316
317 dig->SetTriggerBit(kL1Jet,0);
318 }
319 }
de39a0ff 320
321 Int_t** reg = fSTU->Region();
322
323 if (!fSTU->GetRawData())
324 {
325 // Update digits w/ L1 time sum
326 // Done in raw digit maker when raw data enabled
327 for (Int_t i = 0; i < 48; i++)
328 {
329 for (Int_t j = 0; j < 64; j++)
330 {
331 if (reg[i][j])
332 {
de39a0ff 333 if (geom->GetAbsFastORIndexFromPositionInEMCAL(i, j, id))
334 {
de39a0ff 335 if (posMap[i][j] == -1)
336 {
337 // Add a new digit with L1 time sum
338 new((*digits)[digits->GetEntriesFast()]) AliEMCALTriggerRawDigit(id, 0x0, 0);
339
340 dig = (AliEMCALTriggerRawDigit*)digits->At(digits->GetEntriesFast() - 1);
341 }
342 else
343 {
344 dig = (AliEMCALTriggerRawDigit*)digits->At(posMap[i][j]);
345 }
346
347 dig->SetL1TimeSum(reg[i][j]);
348 }
349 }
350 }
351 }
352 }
916f1e76 353 }
354
de39a0ff 355 if (AliDebugLevel() >= 999) data->Scan();
39c05eac 356
916f1e76 357 // Now reset the electronics for a fresh start with next event
358 Reset();
359}
360
361//__________________
362void AliEMCALTriggerElectronics::Reset()
363{
364 //
365 TIter NextTRU(fTRU);
366 while ( AliEMCALTriggerTRU *TRU = (AliEMCALTriggerTRU*)NextTRU() ) TRU->Reset();
367
368 fSTU->Reset();
369}