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