MakeImage is now a method of AliCheckerBase, was AliQADataMaker before. This will...
[u/mrichter/AliRoot.git] / FMD / AliFMDRawReader.cxx
CommitLineData
e802be3e 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 **************************************************************************/
e802be3e 15/* $Id$ */
c2fc1258 16/** @file AliFMDRawReader.cxx
17 @author Christian Holm Christensen <cholm@nbi.dk>
18 @date Mon Mar 27 12:45:23 2006
19 @brief Class to read raw data
02a27b50 20 @ingroup FMD_rec
c2fc1258 21*/
e802be3e 22//____________________________________________________________________
23//
24// Class to read ADC values from a AliRawReader object.
25//
7684b53c 26// This class uses the AliFMDRawStreamer class to read the ALTRO
27// formatted data.
28//
29// +-------+
30// | TTask |
31// +-------+
32// ^
33// |
34// +-----------------+ <<references>> +--------------+
35// | AliFMDRawReader |<>----------------| AliRawReader |
36// +-----------------+ +--------------+
37// | ^
38// | <<uses>> |
39// V |
40// +-----------------+ <<uses>> |
41// | AliFMDRawStream |------------------------+
42// +-----------------+
43// |
44// V
45// +----------------+
46// | AliAltroStream |
47// +----------------+
48//
f95a63c4 49// #include <AliLog.h> // ALILOG_H
50#include "AliFMDDebug.h" // Better debug macros
1a1fdef7 51#include "AliFMDParameters.h" // ALIFMDPARAMETERS_H
e802be3e 52#include "AliFMDDigit.h" // ALIFMDDIGIT_H
faf80567 53#include "AliFMDSDigit.h" // ALIFMDSDIGIT_H
e802be3e 54#include "AliFMDRawStream.h" // ALIFMDRAWSTREAM_H
3effc6e7 55#include "AliRawReader.h" // ALIRAWREADER_H
e802be3e 56#include "AliFMDRawReader.h" // ALIFMDRAWREADER_H
3effc6e7 57#include "AliFMDDebug.h"
58#include "AliFMDCalibSampleRate.h"
59#include "AliFMDCalibStripRange.h"
b995fc28 60#include "AliFMDAltroMapping.h"
c2fc1258 61// #include "AliFMDAltroIO.h" // ALIFMDALTROIO_H
3effc6e7 62#include <TArrayS.h> // ROOT_TArrayS
1a1fdef7 63#include <TTree.h> // ROOT_TTree
64#include <TClonesArray.h> // ROOT_TClonesArray
3effc6e7 65#include <TString.h>
7af3df7f 66#include <iostream>
67#include <climits>
02a27b50 68// #include <iomanip>
e802be3e 69
70//____________________________________________________________________
925e6570 71ClassImp(AliFMDRawReader)
1a1fdef7 72#if 0
73 ; // This is here to keep Emacs for indenting the next line
74#endif
e802be3e 75
76//____________________________________________________________________
1a1fdef7 77AliFMDRawReader::AliFMDRawReader(AliRawReader* reader, TTree* tree)
e802be3e 78 : TTask("FMDRawReader", "Reader of Raw ADC values from the FMD"),
1a1fdef7 79 fTree(tree),
80 fReader(reader),
faf80567 81 // fSampleRate(1),
3effc6e7 82 fData(0),
f38b1653 83 fNbytes(0),
84 fSeen()
e802be3e 85{
56b1929b 86 // Default CTOR
faf80567 87 for (Int_t i = 0; i < 3; i++) {
88 fSampleRate[i] = 0;
89 fZeroSuppress[i] = kFALSE;
90 fNoiseFactor[i] = 1;
91 }
e802be3e 92}
93
e802be3e 94//____________________________________________________________________
c2fc1258 95void
96AliFMDRawReader::Exec(Option_t*)
97{
02a27b50 98 // Read the data
c2fc1258 99 TClonesArray* array = new TClonesArray("AliFMDDigit");
100 if (!fTree) {
101 AliError("No tree");
102 return;
103 }
104 fTree->Branch("FMD", &array);
3effc6e7 105
106
c2fc1258 107 ReadAdcs(array);
108 Int_t nWrite = fTree->Fill();
f95a63c4 109 AliFMDDebug(1, ("Got a grand total of %d digits, wrote %d bytes to tree",
bf37515f 110 array->GetEntriesFast(), nWrite));
c2fc1258 111}
112
b5014544 113//____________________________________________________________________
114Bool_t
faf80567 115AliFMDRawReader::NextSample(UShort_t& det, Char_t& rng, UShort_t& sec,
116 UShort_t& str, UShort_t& sam, UShort_t& rat,
117 Short_t& adc, Bool_t& zs, UShort_t& fac)
b5014544 118{
119 // Scan current event for next signal. It returns kFALSE when
120 // there's no more data in the event.
121 static AliAltroRawStream stream(fReader); // = 0;
122 static AliFMDParameters* pars = 0;
123 static AliFMDAltroMapping* map = 0;
124 static Int_t ddl = -1;
faf80567 125 // static UInt_t rate = 0;
b5014544 126 static UShort_t tdet = 0;
127 static Char_t trng = '\0';
128 static UShort_t tsec = 0;
129 static Short_t tstr = 0;
130 static Short_t bstr = -1;
faf80567 131 static Short_t tsam = -1;
132 static UInt_t trate = 0;
b5014544 133 static Int_t hwaddr = -1;
134 static UShort_t stripMin = 0;
135 static UShort_t stripMax = 0; // 127;
136 static UShort_t preSamp = 0; // 14+5;
137 if (stream.GetDDLNumber() < 0) {
138 fReader->Select("FMD");
139
140 // Reset "seen" array
141 // const UShort_t kUShortMax = (1 << 16) - 1;
142 // fSeen.Reset(kUShortMax);
143
144 pars = AliFMDParameters::Instance();
145 map = pars->GetAltroMap();
146 // stream = new AliAltroRawStream(fReader);
147
148 AliFMDDebug(5, ("Setting %d word headers",
149 pars->HasCompleteHeader() ? 8 : 7));
150 stream.SetShortDataHeader(!pars->HasCompleteHeader());
b5014544 151
152 // Reset variables
153 ddl = -1;
faf80567 154 trate = 0;
b5014544 155 tdet = 0;
156 trng = '\0';
157 tsec = 0;
158 tstr = 0;
faf80567 159 tsam = -1;
b5014544 160 hwaddr = -1;
161 }
162 do {
163 Bool_t next = stream.Next();
164 if (!next) {
165 // if (stream) delete stream;
166 // stream = 0;
167 return kFALSE;
168 }
169 Int_t thisDDL = stream.GetDDLNumber();
170 AliFMDDebug(10, ("RCU @ DDL %d", thisDDL));
171 if (thisDDL != ddl) {
faf80567 172 ddl = thisDDL;
173 fZeroSuppress[ddl] = zs = stream.GetZeroSupp();
174 fNoiseFactor[ddl] = fac = stream.GetNPostsamples();
175 fSampleRate[ddl] = trate = 0; // stream.GetNPresamples();
176 tdet = map->DDL2Detector(ddl);
b5014544 177 AliFMDDebug(10, ("RCU @ DDL %d zero suppression: %s",ddl, zs?"yes":"no"));
178 AliFMDDebug(10, ("RCU @ DDL %d noise factor: %d", ddl,fac));
faf80567 179 AliFMDDebug(10, ("RCU @ DDL %d sample rate: %d", ddl, trate));
b5014544 180 }
181 Int_t thisAddr = stream.GetHWAddress();
182 AliFMDDebug(10, ("RCU @ DDL %d, Address 0x%03x", ddl, thisAddr));
183 if (thisAddr != hwaddr) {
184 hwaddr = thisAddr;
185 UShort_t board, chip, channel;
186 map->ChannelAddress(hwaddr, board, chip, channel);
187 map->Channel2StripBase(board, chip, channel, trng, tsec, bstr);
188 AliFMDDebug(10, ("0x%04x/0x%03x maps to FMD%d%c[%2d]-%3d",
189 ddl, hwaddr, tdet, trng, tsec, bstr));
190
191 stripMin = pars->GetMinStrip(tdet, trng, tsec, bstr);
192 stripMax = pars->GetMaxStrip(tdet, trng, tsec, bstr);
193 preSamp = pars->GetPreSamples(tdet, trng, tsec, bstr);
faf80567 194 if (trate == 0)
195 fSampleRate[ddl] = trate = pars->GetSampleRate(tdet, trng, tsec, bstr);
b5014544 196 AliFMDDebug(10, ("RCU @ DDL %d, Address 0x%03x sample rate: %d",
faf80567 197 ddl, hwaddr, trate));
b5014544 198
199 Int_t nChAddrMismatch = stream.GetNChAddrMismatch();
200 Int_t nChLenMismatch = stream.GetNChLengthMismatch();
201 if (nChAddrMismatch != 0)
202 AliWarning(Form("Got %d channels with address mis-matches for 0x%03x",
203 nChAddrMismatch, hwaddr));
204 if (nChLenMismatch != 0)
205 AliWarning(Form("Got %d channels with length mis-matches for 0x%03x",
206 nChLenMismatch, hwaddr));
207 }
208 // Get the signal
209 adc = stream.GetSignal();
210
211 // Sanity check - if the total bunch length is less than 1, then
212 // read until we get the next bunch.
213 Int_t b = stream.GetTimeLength();
214 if (b < 1) {
215 AliWarning(Form("Bunch length %0d is less than 0 for "
216 "DDL %4d address 0x%03x", b, ddl, hwaddr));
217 continue;
218 }
219 // Sanity check - if the current time is less than 0, then read
220 // until we get a new bunch.
221 Int_t t = stream.GetTime();
222 if (t < 0) {
223 AliWarning(Form("Time %0d is less than 0 for DDL %4d address 0x%03x",
224 t, ddl, hwaddr));
225 continue;
226 }
227
228 Short_t strOff = 0;
229 UShort_t samp = 0;
faf80567 230 map->Timebin2Strip(tsec, t, preSamp, trate, strOff, samp);
b5014544 231 tstr = bstr + strOff;
faf80567 232 tsam = samp;
b5014544 233 AliFMDDebug(20, ("0x%04x/0x%03x/%04d maps to FMD%d%c[%2d,%3d]-%d",
234 ddl, hwaddr, t, tdet, trng, tsec, tstr, samp));
235
b5014544 236 // Local strip number to channel
237 Short_t l = (tstr > bstr ? tstr - bstr : bstr - tstr);
238 AliFMDDebug(10, ("Checking if strip %d in range [%d,%d]",
239 l, stripMin, stripMax));
240 if (l < stripMin || l > stripMax) {
241 AliFMDDebug(5, ("VA channel %3d (t: %4d) of DDL %4d address 0x%03x "
242 "is out of range (%3d->%3d)",
243 l, t, ddl, hwaddr, stripMin, stripMax));
244 continue;
245 }
246
247
248 det = tdet;
249 rng = trng;
250 sec = tsec;
251 str = tstr;
faf80567 252 sam = tsam;
253 rat = trate;
b5014544 254 // adc = stream.GetSignal();
255
256 break;
257 } while (true);
258 return kTRUE;
259}
260
faf80567 261//____________________________________________________________________
262Bool_t
263AliFMDRawReader::NextSignal(UShort_t& det, Char_t& rng,
264 UShort_t& sec, UShort_t& str,
265 Short_t& adc, Bool_t& zs,
266 UShort_t& fac)
267{
268
269 do {
270 UShort_t samp, rate;
271 if (!NextSample(det, rng, sec, str, samp, rate, adc, zs, fac))
272 return kFALSE;
273
274 Bool_t take = kFALSE;
275 switch (rate) {
276 case 1: take = kTRUE; break;
277 case 2: if (samp == 1) take = kTRUE; break;
278 case 3: if (samp == 1) take = kTRUE; break;
279 case 4: if (samp == 2) take = kTRUE; break;
280 default: if (samp == rate-2) take = kTRUE; break;
281 }
282 if (!take) continue;
283 break;
284 } while (true);
285 return kTRUE;
286}
287
288//____________________________________________________________________
289Bool_t
290AliFMDRawReader::SelectSample(UShort_t samp, UShort_t rate)
291{
292 Bool_t take = kFALSE;
293 switch (rate) {
294 case 1: take = kTRUE; break;
295 case 2: if (samp == 1) take = kTRUE; break;
296 case 3: if (samp == 1) take = kTRUE; break;
297 case 4: if (samp == 2) take = kTRUE; break;
298 default: if (samp == rate-2) take = kTRUE; break;
299 }
300
301 return take;
302}
303
c2fc1258 304#if 1
305//____________________________________________________________________
306Bool_t
307AliFMDRawReader::ReadAdcs(TClonesArray* array)
308{
fd8d0298 309
310 fReader->Select("FMD");
c2fc1258 311 // Read raw data into the digits array, using AliFMDAltroReader.
312 if (!array) {
313 AliError("No TClonesArray passed");
314 return kFALSE;
315 }
f38b1653 316 const UShort_t kUShortMax = (1 << 16) - 1;
317 fSeen.Reset(kUShortMax);
318
d56f2b53 319 // if (!fReader->ReadHeader()) {
320 // AliError("Couldn't read header");
321 // return kFALSE;
322 // }
c2fc1258 323 // Get sample rate
b995fc28 324 AliFMDParameters* pars = AliFMDParameters::Instance();
325 AliFMDAltroMapping* map = pars->GetAltroMap();
326 AliFMDRawStream input(fReader);
9a090ccd 327 AliFMDDebug(5, ("Setting 7 word headers"));
625d7886 328 input.SetShortDataHeader(!pars->HasCompleteHeader());
c2fc1258 329
330 UShort_t stripMin = 0;
f38b1653 331 UShort_t stripMax = 0; // 127;
332 UShort_t preSamp = 0; // 14+5;
5cf05dbb 333
334 Int_t oldddl = -1;
c2fc1258 335 UInt_t ddl = 0;
faf80567 336 // UInt_t rate = 0;
c2fc1258 337 UInt_t last = 0;
338 UInt_t hwaddr = 0;
339 // Data array is approx twice the size needed.
340 UShort_t data[2048];
5cf05dbb 341 for (size_t i = 0; i < 2048; i++) data[i] = 0; // kUShortMax;
2b893216 342
625d7886 343 Bool_t isGood = kTRUE;
344 while (isGood) {
345 isGood = input.ReadChannel(ddl, hwaddr, last, data);
7868d9e2 346 // if (!isGood) break;
7af3df7f 347 if (ddl >= UINT_MAX /* UInt_t(-1) */) {
fd8d0298 348 AliFMDDebug(5, ("At end of event with %d digits",
349 array->GetEntriesFast()));
350 break;
351 }
7af3df7f 352 if (UInt_t(oldddl) != ddl) {
5cf05dbb 353 fZeroSuppress[ddl] = input.GetZeroSupp();
354 AliFMDDebug(20, ("RCU @ DDL %d zero suppression: %s",
355 ddl, (fZeroSuppress[ddl] ? "yes" : "no")));
2b893216 356
5cf05dbb 357 // WARNING: We store the noise factor in the 2nd baseline
358 // filters excluded post samples, since we'll never use that
359 // mode.
360 fNoiseFactor[ddl] = input.GetNPostsamples();
361 AliFMDDebug(20, ("RCU @ DDL %d noise factor: %d", ddl,fNoiseFactor[ddl]));
362
faf80567 363 // WARNING: We store the noise factor in the 2nd baseline
364 // filters excluded post samples, since we'll never use that
365 // mode.
366 fSampleRate[ddl] = input.GetNPretriggerSamples();
367 AliFMDDebug(20, ("RCU @ DDL %d Sample rate: %d", ddl,fNoiseFactor[ddl]));
368
5cf05dbb 369 Int_t nChAddrMismatch = input.GetNChAddrMismatch();
370 Int_t nChLenMismatch = input.GetNChLengthMismatch();
371 if (nChAddrMismatch != 0)
372 AliWarning(Form("Got %d channels with address mis-matches for 0x%03x",
373 nChAddrMismatch, hwaddr));
374 if (nChLenMismatch != 0)
375 AliWarning(Form("Got %d channels with length mis-matches for 0x%03x",
376 nChLenMismatch, hwaddr));
377 oldddl = ddl;
378 }
379 // AliFMDDebug(5, ("Read channel 0x%x of size %d", hwaddr, last));
f38b1653 380
b995fc28 381 UShort_t det, sec, samp, board, chip, channel;
1f12303f 382 Short_t strbase;
383 Char_t ring;
fd8d0298 384
b995fc28 385
386 if (map->DDL2Detector(ddl) < 0) break;
387 det = map->DDL2Detector(ddl);
388 map->ChannelAddress(hwaddr, board, chip, channel);
389 if (!map->Channel2StripBase(board, chip, channel, ring, sec, strbase)) {
1f12303f 390 AliError(Form("Failed to get detector id from DDL %d, "
391 "hardware address 0x%03x", ddl, hwaddr));
392 continue;
393 }
faf80567 394 AliFMDDebug(5, ("Board: 0x%02x, Altro: 0x%x, Channel: 0x%x, Length: %4d",
7af3df7f 395 board, chip, channel, last));
1f12303f 396
397 stripMin = pars->GetMinStrip(det, ring, sec, strbase);
398 stripMax = pars->GetMaxStrip(det, ring, sec, strbase);
399 preSamp = pars->GetPreSamples(det, ring, sec, strbase);
faf80567 400 if (fSampleRate[ddl] == 0)
401 fSampleRate[ddl] = pars->GetSampleRate(det, ring, sec, strbase);
c2fc1258 402
403 // Loop over the `timebins', and make the digits
404 for (size_t i = 0; i < last; i++) {
f38b1653 405 // if (i < preSamp) continue;
5cf05dbb 406 AliFMDDebug(15, ("0x%04x/0x%03x/%04d %4d", ddl, hwaddr, i, data[i]));
f38b1653 407
b995fc28 408 Short_t stroff = 0;
faf80567 409 map->Timebin2Strip(sec, i, preSamp, fSampleRate[ddl], stroff, samp);
b995fc28 410 Short_t str = strbase + stroff;
1f12303f 411
412 AliFMDDebug(10, ("0x%04x/0x%03x/%04d maps to FMD%d%c[%2d,%3d]-%d",
413 ddl, hwaddr, i, det, ring, sec, str, samp));
f38b1653 414 if (str < 0) {
415 AliFMDDebug(8, ("Got presamples at timebin %d", i));
5cf05dbb 416 data[i] = 0; // Reset cache
f38b1653 417 continue;
418 }
419
faf80567 420 Short_t lstrip = (i - preSamp) / fSampleRate[ddl] + stripMin;
f38b1653 421
422 AliFMDDebug(15, ("Checking if strip %d (%d) in range [%d,%d]",
423 lstrip, str, stripMin, stripMax));
424 if (lstrip < stripMin || lstrip > stripMax) {
425 AliFMDDebug(5, ("FMD%d%c[%02d,%03d]-%d out of range (%3d->%3d)",
426 det, ring, sec, samp, str, stripMin, stripMax));
5cf05dbb 427 data[i] = 0; // Reset cache
f38b1653 428 continue;
429 }
faf80567 430 // Possibly do pedestal subtraction of signal
431 Int_t counts = data[i];
432
f38b1653 433
434 // Check the cache of indicies
435 Int_t idx = fSeen(det, ring, sec, str);
436 if (idx == kUShortMax) {
437 // We haven't seen this strip yet.
438 fSeen(det, ring, sec, str) = idx = array->GetEntriesFast();
439 AliFMDDebug(7,("making digit for FMD%d%c[%2d,%3d]-%d from timebin %4d",
440 det, ring, sec, str, samp, i));
441 new ((*array)[idx]) AliFMDDigit(det, ring, sec, str);
442 }
faf80567 443 AliFMDBaseDigit* digit = static_cast<AliFMDBaseDigit*>(array->At(idx));
f38b1653 444 AliFMDDebug(10,
faf80567 445 ("Setting FMD%d%c[%2d,%3d]-%d from timebin %4d=%4d (%4d)",
446 det, ring, sec, str, samp, i, counts, data[i]));
447 digit->SetCount(samp, counts);
5cf05dbb 448 data[i] = 0; // Reset cache
c2fc1258 449 }
450 }
451 return kTRUE;
452}
453#else
454//____________________________________________________________________
d760ea03 455Bool_t
456AliFMDRawReader::ReadAdcs(TClonesArray* array)
e802be3e 457{
1e8f773e 458 // Read raw data into the digits array, using AliFMDAltroReader.
d760ea03 459 if (!array) {
460 AliError("No TClonesArray passed");
461 return kFALSE;
462 }
d56f2b53 463 // if (!fReader->ReadHeader()) {
464 // AliError("Couldn't read header");
465 // return kFALSE;
466 // }
1e8f773e 467 // Get sample rate
468 AliFMDParameters* pars = AliFMDParameters::Instance();
469
470 // Select FMD DDL's
362c9d61 471 fReader->Select("FMD");
1e8f773e 472
473 UShort_t stripMin = 0;
474 UShort_t stripMax = 127;
475 UShort_t preSamp = 0;
476
477 do {
478 UChar_t* cdata;
479 if (!fReader->ReadNextData(cdata)) break;
480 size_t nchar = fReader->GetDataSize();
362c9d61 481 UShort_t ddl = fReader->GetDDLID();
c2fc1258 482 UShort_t rate = 0;
f95a63c4 483 AliFMDDebug(1, ("Reading %d bytes (%d 10bit words) from %d",
1e8f773e 484 nchar, nchar * 8 / 10, ddl));
485 // Make a stream to read from
486 std::string str((char*)(cdata), nchar);
487 std::istringstream s(str);
488 // Prep the reader class.
489 AliFMDAltroReader r(s);
490 // Data array is approx twice the size needed.
491 UShort_t data[2048], hwaddr, last;
492 while (r.ReadChannel(hwaddr, last, data) > 0) {
f95a63c4 493 AliFMDDebug(5, ("Read channel 0x%x of size %d", hwaddr, last));
1e8f773e 494 UShort_t det, sec, str;
495 Char_t ring;
496 if (!pars->Hardware2Detector(ddl, hwaddr, det, ring, sec, str)) {
497 AliError(Form("Failed to detector id from DDL %d "
498 "and hardware address 0x%x", ddl, hwaddr));
499 continue;
500 }
c2fc1258 501 rate = pars->GetSampleRate(det, ring, sec, str);
502 stripMin = pars->GetMinStrip(det, ring, sec, str);
503 stripMax = pars->GetMaxStrip(det, ring, sec, str);
f95a63c4 504 AliFMDDebug(5, ("DDL 0x%04x, address 0x%03x maps to FMD%d%c[%2d,%3d]",
1e8f773e 505 ddl, hwaddr, det, ring, sec, str));
506
507 // Loop over the `timebins', and make the digits
508 for (size_t i = 0; i < last; i++) {
509 if (i < preSamp) continue;
d760ea03 510 Int_t n = array->GetEntries();
1e8f773e 511 UShort_t curStr = str + stripMin + i / rate;
512 if ((curStr-str) > stripMax) {
513 AliError(Form("Current strip is %d but DB says max is %d",
514 curStr, stripMax));
515 }
f95a63c4 516 AliFMDDebug(5, ("making digit for FMD%d%c[%2d,%3d] from sample %4d",
1e8f773e 517 det, ring, sec, curStr, i));
518 new ((*array)[n]) AliFMDDigit(det, ring, sec, curStr, data[i],
519 (rate >= 2 ? data[i+1] : 0),
520 (rate >= 3 ? data[i+2] : 0));
521 if (rate >= 2) i++;
522 if (rate >= 3) i++;
523 }
524 if (r.IsBof()) break;
525 }
526 } while (true);
d760ea03 527 return kTRUE;
528}
529
530
531
1e8f773e 532// This is the old method, for comparison. It's really ugly, and far
533// too convoluted.
534//____________________________________________________________________
535void
536AliFMDRawReader::Exec(Option_t*)
537{
56b1929b 538 // Read raw data into the digits array
d56f2b53 539 // if (!fReader->ReadHeader()) {
540 // Error("ReadAdcs", "Couldn't read header");
541 // return;
542 // }
e802be3e 543
1a1fdef7 544 Int_t n = 0;
545 TClonesArray* array = new TClonesArray("AliFMDDigit");
546 fTree->Branch("FMD", &array);
547
8f6ee336 548 // Get sample rate
549 AliFMDParameters* pars = AliFMDParameters::Instance();
362c9d61 550 fSampleRate = pars->GetSampleRate(0);
8f6ee336 551
e802be3e 552 // Use AliAltroRawStream to read the ALTRO format. No need to
553 // reinvent the wheel :-)
554 AliFMDRawStream input(fReader, fSampleRate);
555 // Select FMD DDL's
362c9d61 556 fReader->Select("FMD");
e802be3e 557
558 Int_t oldDDL = -1;
559 Int_t count = 0;
560 UShort_t detector = 1; // Must be one here
561 UShort_t oldDetector = 0;
562 Bool_t next = kTRUE;
563
564 // local Cache
565 TArrayI counts(10);
566 counts.Reset(-1);
567
568 // Loop over data in file
569 while (next) {
570 next = input.Next();
571
572 count++;
573 Int_t ddl = fReader->GetDDLID();
f95a63c4 574 AliFMDDebug(10, ("Current DDL is %d", ddl));
e802be3e 575 if (ddl != oldDDL || input.IsNewStrip() || !next) {
576 // Make a new digit, if we have some data (oldDetector == 0,
577 // means that we haven't really read anything yet - that is,
578 // it's the first time we get here).
579 if (oldDetector > 0) {
580 // Got a new strip.
f95a63c4 581 AliFMDDebug(10, ("Add a new strip: FMD%d%c[%2d,%3d] "
e802be3e 582 "(current: FMD%d%c[%2d,%3d])",
583 oldDetector, input.PrevRing(),
584 input.PrevSector() , input.PrevStrip(),
585 detector , input.Ring(), input.Sector(),
586 input.Strip()));
1a1fdef7 587 new ((*array)[n]) AliFMDDigit(oldDetector,
588 input.PrevRing(),
589 input.PrevSector(),
590 input.PrevStrip(),
591 counts[0], counts[1], counts[2]);
592 n++;
e802be3e 593#if 0
594 AliFMDDigit* digit =
595 static_cast<AliFMDDigit*>(fFMD->Digits()->
596 UncheckedAt(fFMD->GetNdigits()-1));
597#endif
598 }
599
600 if (!next) {
f95a63c4 601 AliFMDDebug(10, ("Read %d channels for FMD%d",
e802be3e 602 count + 1, detector));
603 break;
604 }
605
606
607 // If we got a new DDL, it means we have a new detector.
608 if (ddl != oldDDL) {
609 if (detector != 0)
f95a63c4 610 AliFMDDebug(10, ("Read %d channels for FMD%d", count + 1, detector));
e802be3e 611 // Reset counts, and update the DDL cache
612 count = 0;
613 oldDDL = ddl;
614 // Check that we're processing a FMD detector
615 Int_t detId = fReader->GetDetectorID();
362c9d61 616 if (detId != (AliDAQ::DetectorID("FMD"))) {
1e8f773e 617 AliError(Form("Detector ID %d != %d",
362c9d61 618 detId, (AliDAQ::DetectorID("FMD"))));
e802be3e 619 break;
620 }
621 // Figure out what detector we're deling with
622 oldDetector = detector;
623 switch (ddl) {
624 case 0: detector = 1; break;
625 case 1: detector = 2; break;
626 case 2: detector = 3; break;
627 default:
1e8f773e 628 AliError(Form("Unknown DDL 0x%x for FMD", ddl));
e802be3e 629 return;
630 }
f95a63c4 631 AliFMDDebug(10, ("Reading ADCs for 0x%x - That is FMD%d",
e802be3e 632 fReader->GetEquipmentId(), detector));
633 }
634 counts.Reset(-1);
635 }
636
637 counts[input.Sample()] = input.Count();
638
f95a63c4 639 AliFMDDebug(10, ("ADC of FMD%d%c[%2d,%3d] += %d",
e802be3e 640 detector, input.Ring(), input.Sector(),
641 input.Strip(), input.Count()));
642 oldDetector = detector;
643 }
1a1fdef7 644 fTree->Fill();
e802be3e 645 return;
646
647}
1e8f773e 648#endif
3effc6e7 649//____________________________________________________________________
650Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate,
651 AliFMDCalibStripRange* stripRange,
652 TArrayS &pulseSize,
f38b1653 653 TArrayS &pulseLength)
654{
3effc6e7 655
656 AliFMDDebug(0, ("Start of SOD/EOD"));
657
658 UInt_t shift_clk[18];
659 UInt_t sample_clk[18];
660 UInt_t strip_low[18];
661 UInt_t strip_high[18];
662 UInt_t pulse_size[18];
663 UInt_t pulse_length[18];
b995fc28 664 AliFMDParameters* param = AliFMDParameters::Instance();
665 AliFMDAltroMapping* map = param->GetAltroMap();
3effc6e7 666
667 while(fReader->ReadNextData(fData)) {
668
669 Int_t ddl = fReader->GetDDLID();
670 Int_t detID = fReader->GetDetectorID();
671 AliFMDDebug(0, (" From reader: DDL number is %d , det ID is %d",ddl,detID));
672
673 fNbytes = fReader->GetDataSize();
674
675 ULong_t nWords = GetNwords();
676 UInt_t trailerLast = Get32bitWord(nWords);
677 AliFMDDebug(20, (" # Bytes: %d, # Words: %d, Last word: 0x%08x",
f38b1653 678 fNbytes, nWords, trailerLast));
3effc6e7 679 if ((trailerLast & 0xFFFF0000) != 0xAAAA0000) {
680 AliWarning(Form("Last word 0x%08x does not match RCU II trailer",
681 trailerLast));
682 return kFALSE;
683 }
684
f38b1653 685 // 7 last bits is size of trailer
686 ULong_t nTrailerWords = trailerLast & 0x7f;
3effc6e7 687 ULong_t nPayloadWords = Get32bitWord(nWords - nTrailerWords+1);
f38b1653 688
3effc6e7 689 AliFMDDebug(20, (" # trailer words: %d, # payload words: %d",
f38b1653 690 nTrailerWords, nPayloadWords));
3effc6e7 691
692
693 for (ULong_t i = 1; i <= nPayloadWords ; i++) {
694 UInt_t payloadWord = Get32bitWord(i);
76d06ffa 695
f38b1653 696 // address is only 24 bit
697 UInt_t address = (0xffffff & payloadWord);
698 UInt_t type = ((address >> 21) & 0xf);
699 UInt_t error = ((address >> 20) & 0x1);
700 UInt_t bcast = ((address >> 18) & 0x1);
701 UInt_t bc_not_altro = ((address >> 17) & 0x1);
702 UInt_t board = ((address >> 12) & 0x1f);
703 UInt_t instruction = 0;
704 UInt_t chip = 0;
705 UInt_t channel = 0;
706 if(bc_not_altro)
707 instruction = address & 0xfff;
708 else {
709 chip = ((address >> 9) & 0x7);
710 channel = ((address >> 5) & 0x5);
711 instruction = (address & 0x1f);
712 }
76d06ffa 713
f38b1653 714 Bool_t readDataWord = kFALSE;
715 switch(type) {
716 case 0x0: // Fec read
717 readDataWord = kTRUE;
718 case 0x1: // Fec cmd
719 case 0x2: // Fec write
720 i++;
721 break;
722 case 0x4: // Loop
723 case 0x5: // Wait
724 break;
725 case 0x6: // End sequence
726 case 0x7: // End Mem
727 i = nPayloadWords + 1;
728 break;
729 default:
730 break;
731 }
76d06ffa 732
f38b1653 733 //Don't read unless we have a FEC_RD
734 if(!readDataWord) continue;
735
736 UInt_t dataWord = Get32bitWord(i);
737 UInt_t data = (0xFFFFF & dataWord) ;
738 //UInt_t data = (0xFFFF & dataWord) ;
3effc6e7 739
f38b1653 740 if(error) {
741 AliWarning(Form("error bit detected at Word 0x%06x; "
742 "error % d, type %d, bc_not_altro %d, "
743 "bcast %d, board 0x%02x, chip 0x%x, "
744 "channel 0x%02x, instruction 0x%03x",
745 address, error, type, bc_not_altro,
746 bcast,board,chip,channel,instruction));
747 //process error
748 continue;
749 }
3effc6e7 750
751
f38b1653 752 switch(instruction) {
3effc6e7 753
f38b1653 754 case 0x01: break; // First ADC T
755 case 0x02: break; // I 3.3 V
756 case 0x03: break; // I 2.5 V altro digital
757 case 0x04: break; // I 2.5 V altro analog
758 case 0x05: break; // I 2.5 V VA
759 case 0x06: break; // First ADC T
760 case 0x07: break; // I 3.3 V
761 case 0x08: break; // I 2.5 V altro digital
762 case 0x09: break; // I 2.5 V altro analog
763 case 0x0A: break; // I 2.5 V VA
764 case 0x2D: break; // Second ADC T
765 case 0x2E: break; // I 1.5 V VA
766 case 0x2F: break; // I -2.0 V
767 case 0x30: break; // I -2.0 V VA
768 case 0x31: break; // 2.5 V Digital driver
769 case 0x32: break; // Second ADC T
770 case 0x33: break; // I 1.5 V VA
771 case 0x34: break; // I -2.0 V
772 case 0x35: break; // I -2.0 V VA
773 case 0x36: break; // 2.5 V Digital driver
774 case 0x37: break; // Third ADC T
775 case 0x38: break; // Temperature sens. 1
776 case 0x39: break; // Temperature sens. 2
777 case 0x3A: break; // U 2.5 altro digital (m)
778 case 0x3B: break; // U 2.5 altro analog (m)
779 case 0x3C: break; // Third ADC T
780 case 0x3D: break; // Temperature sens. 1
781 case 0x3E: break; // Temperature sens. 2
782 case 0x3F: break; // U 2.5 altro digital (m)
783 case 0x40: break; // U 2.5 altro analog (m)
784 case 0x41: break; // Forth ADC T
785 case 0x42: break; // U 2.5 VA (m)
786 case 0x43: break; // U 1.5 VA (m)
787 case 0x44: break; // U -2.0 VA (m)
788 case 0x45: break; // U -2.0 (m)
789 case 0x46: break; // Forth ADC T
790 case 0x47: break; // U 2.5 VA (m)
791 case 0x48: break; // U 1.5 VA (m)
792 case 0x49: break; // U -2.0 VA (m)
793 case 0x4A: break; // U -2.0 (m)
794 // Counters
795 case 0x0B: break; // L1 trigger CouNTer
796 case 0x0C: break; // L2 trigger CouNTer
797 case 0x0D: break; // Sampling CLK CouNTer
798 case 0x0E: break; // DSTB CouNTer
799 // Test mode
800 case 0x0F: break; // Test mode word
801 case 0x10: break; // Undersampling ratio.
802 // Configuration and status
803 case 0x11: break; // Config/Status Register 0
804 case 0x12: break; // Config/Status Register 1
805 case 0x13: break; // Config/Status Register 2
806 case 0x14: break; // Config/Status Register 3
807 case 0x15: break; // Free
808 // Comands:
809 case 0x16: break; // Latch L1, L2, SCLK Counters
810 case 0x17: break; // Clear counters
811 case 0x18: break; // Clear CSR1
812 case 0x19: break; // rstb ALTROs
813 case 0x1A: break; // rstb BC
814 case 0x1B: break; // Start conversion
815 case 0x1C: break; // Scan event length
816 case 0x1D: break; // Read event length
817 case 0x1E: break; // Start test mode
818 case 0x1F: break; // Read acquisition memory
819 // FMD
820 case 0x20: break; // FMDD status
821 case 0x21: break; // L0 counters
822 case 0x22: break; // FMD: Wait to hold
823 case 0x23: break; // FMD: L1 timeout
824 case 0x24: break; // FMD: L2 timeout
825 case 0x25: // FMD: Shift clk
826 shift_clk[board] = ((data >> 8 ) & 0xFF);
827 break;
828 case 0x26: // FMD: Strips
829 strip_low[board] = ((data >> 0 ) & 0xFF);
830 strip_high[board] = ((data >> 8 ) & 0xFF);
831 break;
832 case 0x27: // FMD: Cal pulse
833 pulse_size[board] = ((data >> 8 ) & 0xFF);
834 break;
835 case 0x28: break; // FMD: Shape bias
836 case 0x29: break; // FMD: Shape ref
837 case 0x2A: break; // FMD: Preamp ref
838 case 0x2B: // FMD: Sample clk
839 sample_clk[board] = ((data >> 8 ) & 0xFF);
840 break;
841 case 0x2C: break; // FMD: Commands
842 case 0x4B: // FMD: Cal events
843 pulse_length[board] = ((data >> 0 ) & 0xFF);
844 break;
845 default: break;
3effc6e7 846
f38b1653 847 }
848 AliFMDDebug(50, ("instruction 0x%x, dataword 0x%x",
849 instruction,dataWord));
3effc6e7 850 }
851
f38b1653 852 UShort_t det,sector;
853 Short_t strip;
3effc6e7 854 Char_t ring;
76d06ffa 855
3effc6e7 856 const UInt_t boards[4] = {0,1,16,17};
857 for(Int_t i=0;i<4;i++) {
f38b1653 858 if(ddl==0 && (i==1 || i==3)) continue;
859
3effc6e7 860 UInt_t chip =0, channel=0;
b995fc28 861 det = map->DDL2Detector(ddl);
862 map->Channel2StripBase(boards[i], chip, channel, ring, sector, strip);
76d06ffa 863
864 UInt_t samplerate = 1;
76d06ffa 865 if(sample_clk[boards[i]] == 0) {
f38b1653 866 if(ddl == 0) {
867 Int_t sample1 = sample_clk[boards[0]];
868 Int_t sample2 = sample_clk[boards[2]];
869 if(sample1) sample_clk[boards[i]] = sample1;
870 else sample_clk[boards[i]] = sample2;
871 }
872 else {
873 Int_t sample1 = sample_clk[boards[0]];
874 Int_t sample2 = sample_clk[boards[1]];
875 Int_t sample3 = sample_clk[boards[2]];
876 Int_t sample4 = sample_clk[boards[3]];
877 Int_t agreement = 0;
878 if(sample1 == sample2) agreement++;
879 if(sample1 == sample3) agreement++;
880 if(sample1 == sample4) agreement++;
881 if(sample2 == sample3) agreement++;
882 if(sample2 == sample4) agreement++;
883 if(sample3 == sample4) agreement++;
76d06ffa 884
f38b1653 885 Int_t idx = 0;
886 if(i<3) idx = i+1;
887 else idx = i-1;
888 if(agreement == 3) {
889 sample_clk[boards[i]] = sample_clk[boards[idx]];
890 shift_clk[boards[i]] = shift_clk[boards[idx]];
891 strip_low[boards[i]] = strip_low[boards[idx]];
892 strip_high[boards[i]] = strip_high[boards[idx]];
893 pulse_length[boards[i]] = pulse_length[boards[idx]];
894 pulse_size[boards[i]] = pulse_size[boards[idx]];
895 AliFMDDebug(0, ("Vote taken for ddl %d, board 0x%x",
896 ddl,boards[i]));
76d06ffa 897 }
f38b1653 898 }
76d06ffa 899 }
900
3effc6e7 901 if(sample_clk[boards[i]])
f38b1653 902 samplerate = shift_clk[boards[i]]/sample_clk[boards[i]];
3effc6e7 903 sampleRate->Set(det,ring,sector,0,samplerate);
f38b1653 904 stripRange->Set(det,ring,sector,0,
905 strip_low[boards[i]],strip_high[boards[i]]);
3effc6e7 906
76d06ffa 907 AliFMDDebug(20, ("det %d, ring %c, ",det,ring));
f38b1653 908 pulseLength.AddAt(pulse_length[boards[i]],
909 GetHalfringIndex(det,ring,boards[i]/16));
910 pulseSize.AddAt(pulse_size[boards[i]],
911 GetHalfringIndex(det,ring,boards[i]/16));
76d06ffa 912
3effc6e7 913
76d06ffa 914
915 AliFMDDebug(20, (": Board: 0x%02x\n"
f38b1653 916 "\tstrip_low %3d, strip_high %3d\n"
917 "\tshift_clk %3d, sample_clk %3d\n"
918 "\tpulse_size %3d, pulse_length %3d",
919 boards[i],
920 strip_low[boards[i]], strip_high[boards[i]],
921 shift_clk[boards[i]], sample_clk[boards[i]],
922 pulse_size[boards[i]],pulse_length[boards[i]]));
3effc6e7 923 }
924
925 }
926
927 AliFMDParameters::Instance()->SetSampleRate(sampleRate);
928 AliFMDParameters::Instance()->SetStripRange(stripRange);
929
930 AliFMDDebug(0, ("End of SOD/EOD"));
931
7af3df7f 932 return kTRUE;
3effc6e7 933}
934//____________________________________________________________________
935
936UInt_t AliFMDRawReader::Get32bitWord(Int_t idx)
937{
938 // This method returns the 32 bit word at a given
939 // position inside the raw data payload.
940 // The 'index' points to the beginning of the next word.
941 // The method is supposed to be endian (platform)
942 // independent.
943 if (!fData) {
944 AliFatal("Raw data paylod buffer is not yet initialized !");
945 }
946
947 Int_t index = 4*idx;
948
949 if (index < 4) {
950 // fRawReader->AddFatalErrorLog(k32bitWordReadErr,Form("pos = %d",index));
951 // PrintDebug();
952 AliWarning(Form("Invalid raw data payload position (%d) !",index));
953 }
954
955 UInt_t word = 0;
76d06ffa 956
3effc6e7 957 word = fData[--index] << 24;
958 word |= fData[--index] << 16;
959 word |= fData[--index] << 8;
960 word |= fData[--index] << 0 ;
e802be3e 961
3effc6e7 962 return word;
963}
964//_____________________________________________________________________
f38b1653 965Int_t AliFMDRawReader::GetHalfringIndex(UShort_t det, Char_t ring,
966 UShort_t board) {
3effc6e7 967
968 UShort_t iring = (ring == 'I' ? 1 : 0);
969
970 Int_t index = (((det-1) << 2) | (iring << 1) | (board << 0));
971
76d06ffa 972 return index-2;
3effc6e7 973
974}
e802be3e 975//____________________________________________________________________
976//
977// EOF
978//