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