Fixed coding convention issues as given by the automatic
[u/mrichter/AliRoot.git] / FMD / AliFMDRawReader.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15 /* $Id$ */
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 
20     @ingroup FMD_rec
21 */
22 //____________________________________________________________________
23 //
24 // Class to read ADC values from a AliRawReader object. 
25 //
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 //
49 // #include <AliLog.h>          // ALILOG_H
50 #include "AliFMDDebug.h" // Better debug macros
51 #include "AliFMDParameters.h"   // ALIFMDPARAMETERS_H
52 #include "AliFMDDigit.h"        // ALIFMDDIGIT_H
53 #include "AliFMDSDigit.h"       // ALIFMDSDIGIT_H
54 // #include "AliFMDRawStream.h" // ALIFMDRAWSTREAM_H 
55 #include "AliRawReader.h"       // ALIRAWREADER_H 
56 #include "AliFMDRawReader.h"    // ALIFMDRAWREADER_H 
57 #include "AliFMDDebug.h"
58 #include "AliFMDCalibSampleRate.h"
59 #include "AliFMDCalibStripRange.h"
60 #include "AliFMDAltroMapping.h"
61 // #include "AliFMDAltroIO.h"   // ALIFMDALTROIO_H 
62 #include "AliAltroRawStreamV3.h"
63 #include <TArrayS.h>            // ROOT_TArrayS
64 #include <TTree.h>              // ROOT_TTree
65 #include <TClonesArray.h>       // ROOT_TClonesArray
66 #include <TString.h>
67 #include <iostream>
68 #include <climits>
69 // #include <iomanip>
70
71 //____________________________________________________________________
72 ClassImp(AliFMDRawReader)
73 #if 0
74   ; // This is here to keep Emacs for indenting the next line
75 #endif
76
77 //____________________________________________________________________
78 AliFMDRawReader::AliFMDRawReader(AliRawReader* reader, TTree* tree) 
79   : TTask("FMDRawReader", "Reader of Raw ADC values from the FMD"),
80     fTree(tree),
81     fReader(reader), 
82     // fSampleRate(1),
83     fData(0),
84     fNbytes(0), 
85     fMinStrip(0),
86     fMaxStrip(127), 
87     fPreSamp(14+5),
88     fSeen(0)
89 {
90   // Default CTOR
91   for (Int_t i = 0; i < 3; i++) { 
92     fSampleRate[i]   = 0;
93     fZeroSuppress[i] = kFALSE;
94     fNoiseFactor[i]  = 1;
95   }
96 }
97
98 //____________________________________________________________________
99 void
100 AliFMDRawReader::Exec(Option_t*) 
101 {
102   // Read the data 
103   TClonesArray* array = new TClonesArray("AliFMDDigit");
104   if (!fTree) {
105     AliError("No tree");
106     return;
107   }
108   fTree->Branch("FMD", &array);
109   
110   
111   ReadAdcs(array);
112   Int_t nWrite = fTree->Fill();
113   AliFMDDebug(1, ("Got a grand total of %d digits, wrote %d bytes to tree", 
114                    array->GetEntriesFast(), nWrite));
115   delete array;
116 }
117
118 //____________________________________________________________________
119 Int_t
120 AliFMDRawReader::NewDDL(AliAltroRawStreamV3& input, UShort_t& det)
121 {
122   // Process a new DDL.  Sets the internal data members fZeroSuppress, 
123   // fSampleRate, and fNoiseFactor based on information in the RCU trailer. 
124   // 
125   // Parameters:
126   //   input Input stream
127   //   det   On return, the detector number
128   // 
129   // Return value:
130   //   negative value in case of problems, the DDL number otherwise
131
132   // Get the DDL number
133   UInt_t ddl = input.GetDDLNumber();
134   AliFMDDebug(2, ("DDL number %d", ddl));
135
136   // Note, previously, the ALTROCFG1 register was interpreted as 
137   // 
138   // Bits    Value    Description
139   //   0- 3      0/1   1st Baseline filter, mode 
140   //   4- 5   Over-1   2nd baseline filter, # of pre-samples
141   //   6- 9   factor   2nd baseline filter, # of post-samples 
142   //  10-          0   2nd baseline filter, enable
143   //  11-12       00   Zero suppression, glitch filter mode
144   //  13-15      001   Zero suppression, # of post samples
145   //  16-17       01   Zero suppression, # of pre  samples
146   //  18         0/1   Zero suppression, enable
147   //
148   // The interpretation used in AliAltroRawStreamerV3 - which
149   // corresponds directly to ALTRO DPCFG register - is 
150   //
151   // Bits    Value  Description
152   //   0- 3    0/1   1st Baseline filter, mode 
153   //   4         0   Polarity (if '1', then "1's inverse")
154   //   5- 6     01   Zero suppression, # of pre samples
155   //   7-10   0001   Zero suppression, # of post samples
156   //  11         0   2nd baseline filter, enable
157   //  12-13     00   Zero suppression, glitch filter mode
158   //  14-16 factor   2nd baseline filter, # of post-samples
159   //  17-18     01   2nd baseline filter, # of pre-samples 
160   //  19       0/1   Zero suppression, enable
161   //
162   //  Writing 'x' for variable values, that means we have the
163   //  following patterns for the 2 cases 
164   //
165   //    bit #  20   16   12    8    4    0
166   //     old    |0x01|0010|00xx|xxxx|xxxx|
167   //     new    |x01x|xx00|0000|1010|xxxx|
168   //
169   //  That means that we can check if bits 10-13 are '1000' or
170   //  '0000', which will tell us if the value was written with the
171   //  new or the old interpretation.    That is, we can check that 
172   //
173   //    if (((altrocfg1 >> 10) & 0x8) == 0x8) { 
174   //      // old interpretation 
175   //    }
176   //    else { 
177   //      // New interpretation 
178   //    }
179   //
180   // That means, that we should never 
181   //
182   //  - change the # of zero suppression post samples 
183   //  - Turn on 2nd baseline correction 
184   //  - Change the zero-suppression glitch filter mode
185   //
186   // This change as introduced in version 1.2 of Rcu++
187   //
188   UInt_t cfg1 = input.GetAltroCFG1();
189   if (((cfg1 >> 10) & 0x8) == 0x8) {
190     UInt_t cfg2 = input.GetAltroCFG2();
191     AliFMDDebug(3, ("We have data from older MiniConf 0x%x cfg2=0x%08x", 
192                     ((cfg1 >> 10) & 0x8), cfg2));
193     fZeroSuppress[ddl] = (cfg1 >>  0) & 0x1;
194     fNoiseFactor[ddl]  = (cfg1 >>  6) & 0xF;
195     fSampleRate[ddl]   = (cfg2 >> 20) & 0xF;
196   }
197   else {
198     AliFMDDebug(3, ("We have data from newer MiniConf 0x%x", 
199                     ((cfg1 >> 10) & 0x8)));
200     fZeroSuppress[ddl] = input.GetZeroSupp();
201     // WARNING: We store the noise factor in the 2nd baseline
202     // filters excluded post samples, since we'll never use that
203     // mode. 
204     fNoiseFactor[ddl]  = input.GetNPostsamples();
205     // WARNING: We store the sample rate in the number of pre-trigger
206     // samples, since we'll never use that mode.
207     fSampleRate[ddl]     = input.GetNPretriggerSamples();
208   }
209   AliFMDDebug(3, ("RCU @ DDL %d zero suppression: %s", 
210                    ddl, (fZeroSuppress[ddl] ? "yes" : "no")));
211   AliFMDDebug(3, ("RCU @ DDL %d noise factor: %d", ddl,fNoiseFactor[ddl]));    
212   AliFMDDebug(3, ("RCU @ DDL %d sample rate: %d", ddl,fSampleRate[ddl]));
213
214
215   // Get Errors seen 
216   Int_t nChAddrMismatch = input.GetNChAddrMismatch();
217   Int_t nChLenMismatch  = input.GetNChLengthMismatch();
218   if (nChAddrMismatch != 0) 
219     AliWarning(Form("Got %d channels with address mis-matches for 0x%03x",
220                     nChAddrMismatch, ddl));
221   if (nChLenMismatch != 0) 
222     AliWarning(Form("Got %d channels with length mis-matches for 0x%03x",
223                     nChLenMismatch, ddl));
224
225   // Map DDL number to the detector number 
226   AliFMDParameters*    pars   = AliFMDParameters::Instance();
227   AliFMDAltroMapping*  map    = pars->GetAltroMap();
228   if (map->DDL2Detector(ddl) < 0) return -1;
229   det = map->DDL2Detector(ddl);
230
231   if (AliLog::GetDebugLevel("FMD", 0) > 5) 
232     input.PrintRCUTrailer();
233   return ddl;
234 }
235
236 //____________________________________________________________________
237 Int_t
238 AliFMDRawReader::NewChannel(const AliAltroRawStreamV3& input,  UShort_t det,
239                             Char_t&  ring, UShort_t& sec, Short_t& strbase)
240 {
241   // Processs a new channel.  Sets the internal data members
242   // fMinStrip, fMaxStrip, and fPreSamp. 
243   //
244   // Parameter:
245   //   input   Input stream
246   //   ring    On return, the ring identifier 
247   //   sec     On return, the sector number
248   //   strbase On return, the strip base
249   // 
250   // Return value
251   //   negative value in case of problems, hardware address otherwise
252
253   // Get the hardware address, and map that to detector coordinates 
254   UShort_t board, chip, channel;
255   Int_t    ddl    = input.GetDDLNumber();
256   Int_t    hwaddr = input.GetHWAddress();
257   if (input.IsChannelBad()) { 
258     AliError(Form("Channel 0x%03x is marked as bad!", hwaddr));
259   }
260   
261   AliFMDParameters*    pars   = AliFMDParameters::Instance();
262   AliFMDAltroMapping*  map    = pars->GetAltroMap();
263   // Map to hardware stuff 
264   map->ChannelAddress(hwaddr, board, chip, channel);
265   // Then try to map to detector address
266   if (!map->Channel2StripBase(board, chip, channel, ring, sec, strbase)) {
267     AliError(Form("Failed to get detector id from DDL %d, "
268                   "hardware address 0x%03x", ddl, hwaddr));
269     return -1;
270   }
271   AliFMDDebug(4, ("Board: 0x%02x, Altro: 0x%x, Channel: 0x%x", 
272                   board, chip, channel));
273
274   // Get the 'conditions'
275   fMinStrip = pars->GetMinStrip(det, ring, sec, strbase);
276   fMaxStrip = pars->GetMaxStrip(det, ring, sec, strbase);
277   fPreSamp  = pars->GetPreSamples(det, ring, sec, strbase);
278   if (fSampleRate[ddl] == 0) 
279     fSampleRate[ddl] = pars->GetSampleRate(det, ring, sec, strbase);
280
281   return hwaddr;
282 }
283
284 //____________________________________________________________________
285 Int_t
286 AliFMDRawReader::NewSample(const AliAltroRawStreamV3& input, 
287                            Int_t i, UShort_t t, UShort_t sec,
288                            UShort_t  strbase, Short_t&  str, UShort_t& samp)
289 {
290   // Process a new timebin
291   // 
292   // Parameters:
293   //   input   Input stream
294   //   i       Index into bunch data
295   //   t       Time
296   //   strbase Base of strip numbers for this channel
297   //   str     On return, the strip number
298   //   samp    On return, the sample number
299   // 
300   // Return value
301   //   negative value in case of problems, ADC value otherwise
302   if (t < fPreSamp) return -1;
303
304   Int_t           ddl  = input.GetDDLNumber();
305   Int_t           hwa  = input.GetHWAddress();
306   const UShort_t* data = input.GetSignals();
307   Short_t         adc  = data[i];
308   AliFMDDebug(10, ("0x%04x/0x%03x/%04d %4d", ddl, hwa, t, adc));
309
310   AliFMDParameters*    pars   = AliFMDParameters::Instance();
311   AliFMDAltroMapping*  map    = pars->GetAltroMap();
312
313   samp            = 0;
314   Short_t  stroff = 0;
315   map->Timebin2Strip(sec, t, fPreSamp, fSampleRate[ddl], stroff, samp);
316   str             = strbase + stroff;
317       
318   AliFMDDebug(20, ("0x%04x/0x%03x/%04d=%4d maps to strip %3d sample %d " 
319                    "(pre: %d, min: %d, max: %d, rate: %d)",
320                   ddl, hwa, t, adc, str, samp, fPreSamp, 
321                   fMinStrip, fMaxStrip, fSampleRate[ddl]));
322   if (str < 0) { 
323     AliFMDDebug(10, ("Got presamples at timebin %d", i));
324     return -1;
325   }
326           
327   // VA1 Local strip number 
328   Short_t lstrip = (t - fPreSamp) / fSampleRate[ddl] + fMinStrip;
329       
330   AliFMDDebug(15, ("Checking if strip %d (%d) in range [%d,%d]", 
331                    lstrip, str, fMinStrip, fMaxStrip));
332   if (lstrip < fMinStrip || lstrip > fMaxStrip) {
333     AliFMDDebug(10, ("Strip %03d-%d (%d,%d) from t=%d out of range (%3d->%3d)", 
334                     str, samp, lstrip, stroff, t, fMinStrip, fMaxStrip));
335     adc = -1;
336   }
337   // Possibly do pedestal subtraction of signal 
338   if (adc > 1023) 
339     AliWarning(Form("ADC value out of range: %4d", adc));
340   return adc;
341 }
342
343 //____________________________________________________________________
344 Bool_t
345 AliFMDRawReader::NextSample(UShort_t& det, Char_t&   rng, UShort_t& sec, 
346                             UShort_t& str, UShort_t& sam, UShort_t& rat, 
347                             Short_t&  adc, Bool_t&   zs,  UShort_t& fac)
348 {
349   // Scan current event for next signal.   It returns kFALSE when
350   // there's no more data in the event. 
351   static AliAltroRawStreamV3 stream(fReader); //    = 0;
352   static Int_t               ddl      = -1;
353   static UShort_t            tdet     = 0;
354   static Char_t              trng     = '\0';
355   static UShort_t            tsec     = 0;
356   static Short_t             tstr     = 0;   
357   static Short_t             bstr     = -1;
358   static UShort_t            tsam     = 0;   
359   static UInt_t              trate    = 0;
360   static Int_t               hwaddr   = -1;
361   static UShort_t            start    = 0;
362   static UShort_t            length   = 0;
363   static Short_t             t        = -1;
364   static Int_t               i        = 0; 
365   // First entry!
366   if (stream.GetDDLNumber() < 0) { 
367     fReader->Select("FMD");
368
369     // Reset variables
370     ddl    = -1;  
371     trate  = 0;   
372     tdet   = 0;   
373     trng   = '\0';
374     tsec   = 0;   
375     tstr   = 0;  
376     tsam   = -1;
377     hwaddr = -1;
378   }
379
380   do { 
381     if (t < start - length + 1) { 
382       if (!stream.NextBunch()) { 
383         if (!stream.NextChannel()) { 
384           if (!stream.NextDDL()) {
385             stream.Reset();
386             return kFALSE;
387           }
388           ddl = NewDDL(stream, tdet);
389         }
390         hwaddr = NewChannel(stream, tdet, trng, tsec, bstr);
391       }
392       start  = stream.GetStartTimeBin();
393       length = stream.GetBunchLength();
394       t      = start;
395       i      = 0;
396     }
397     Int_t tadc = NewSample(stream, i, t, tsec, bstr, tstr, tsam);
398     if (tadc >= 0) { 
399       det = tdet;
400       rng = trng;
401       sec = tsec;
402       str = tstr;
403       sam = tsam;
404       adc = tadc;
405       rat = fSampleRate[ddl];
406       zs  = fZeroSuppress[ddl];
407       fac = fNoiseFactor[ddl];
408       t--;
409       i++;
410       break;
411     }
412   } while (true);
413
414   return kTRUE;
415 }
416
417
418 //____________________________________________________________________
419 Bool_t
420 AliFMDRawReader::NextSignal(UShort_t& det, Char_t&   rng, 
421                             UShort_t& sec, UShort_t& str, 
422                             Short_t&  adc, Bool_t&   zs, 
423                             UShort_t& fac)
424 {
425   // 
426   // Get the next signal
427   // 
428   // Parameters:
429   //    det  On return, the detector
430   //    rng  On return, the ring
431   //    sec  On return, the sector
432   //    str  On return, the strip
433   //    adc  On return, the ADC value
434   //    zs   On return, whether zero-supp. is enabled
435   //    fac  On return, the usd noise factor
436   // 
437   // Return:
438   //    true if valid data is returned
439   //
440   
441   do { 
442     UShort_t samp, rate;
443     if (!NextSample(det, rng, sec, str, samp, rate, adc, zs, fac)) 
444       return kFALSE;
445
446     Bool_t take = kFALSE;
447     switch (rate) { 
448     case 1:                      take = kTRUE; break;
449     case 2:  if (samp == 1)      take = kTRUE; break;
450     case 3:  if (samp == 1)      take = kTRUE; break; 
451     case 4:  if (samp == 2)      take = kTRUE; break;
452     default: if (samp == rate-2) take = kTRUE; break;
453     }
454     if (!take) continue;
455     break;
456   } while (true);
457   return kTRUE;
458 }
459
460 //____________________________________________________________________
461 Bool_t
462 AliFMDRawReader::SelectSample(UShort_t samp, UShort_t rate) 
463 {
464   // Check if the passed sample is the one we need
465   Bool_t take = kFALSE;
466   switch (rate) { 
467   case 1:                      take = kTRUE; break;
468   case 2:  if (samp == 1)      take = kTRUE; break;
469   case 3:  if (samp == 1)      take = kTRUE; break; 
470   case 4:  if (samp == 2)      take = kTRUE; break;
471   default: if (samp == rate-2) take = kTRUE; break;
472   }
473   
474   return take;
475 }
476   
477 //____________________________________________________________________
478 Bool_t
479 AliFMDRawReader::ReadAdcs(TClonesArray* array) 
480 {
481   // Read ADC values from raw input into passed TClonesArray of AliFMDDigit
482   // objects. 
483   AliFMDDebug(2, ("Reading ADC values into a TClonesArray"));
484
485   // Read raw data into the digits array, using AliFMDAltroReader. 
486   if (!array) {
487     AliError("No TClonesArray passed");
488     return kFALSE;
489   }
490   const UShort_t kUShortMax = (1 << 16) - 1;
491   fSeen.Reset(kUShortMax);
492
493   AliAltroRawStreamV3  input(fReader);
494   input.Reset();
495   input.SelectRawData("FMD");
496   
497   // Loop over input RORCs
498   while (input.NextDDL()) { 
499     UShort_t det = 0;
500     Int_t    ddl = NewDDL(input, det);
501     if (ddl < 0) break;
502
503     while (input.NextChannel()) { 
504       // Get the hardware address, and map that to detector coordinates 
505       Char_t   ring;
506       UShort_t sec;
507       Short_t  strbase;
508       Int_t    hwaddr = NewChannel(input, det, ring, sec, strbase);
509       if (hwaddr < 0) break;
510
511       // Loop over bunches 
512       while (input.NextBunch()) { 
513         // Get Lenght of bunch, and pointer to the data 
514         const UShort_t* data   = input.GetSignals();
515         UShort_t        start  = input.GetStartTimeBin();
516         UShort_t        length = input.GetBunchLength();
517       
518         // Loop over the data and store it. 
519         for (Int_t i = 0; i < length; i++) { 
520           // Time 
521           Short_t  str;
522           UShort_t samp;
523           Int_t    t    = start - i;
524           Int_t    adc  = NewSample(input, i, t, sec, strbase, str, samp);
525           if (adc < 0) continue;
526           UShort_t counts = adc;
527       
528           AliFMDDebug(10, ("FMD%d%c[%02d,%03d]-%d: %4d", 
529                            det, ring, sec, str, samp, counts));
530           // Check the cache of indicies
531           Int_t idx = fSeen(det, ring, sec, str);
532           AliFMDDigit* digit = 0;
533           if (idx == kUShortMax) { 
534             // We haven't seen this strip yet. 
535             fSeen(det, ring, sec, str) = idx = array->GetEntriesFast();
536             AliFMDDebug(7,("making digit for FMD%d%c[%2d,%3d]-%d "
537                            "from timebin %4d", 
538                            det, ring, sec, str, samp, t));
539             digit = new ((*array)[idx]) AliFMDDigit(det, ring, sec, str);
540             digit->SetDefaultCounts(fSampleRate[ddl]);
541           }
542           else 
543             digit = static_cast<AliFMDDigit*>(array->At(idx));
544           AliFMDDebug(10, ("Setting FMD%d%c[%2d,%3d]-%d "
545                            "from timebin %4d=%4d (%4d)", 
546                            det, ring, sec, str, samp, t, counts, data[i]));
547           digit->SetCount(samp, counts);
548         } // for (i)
549       } // while (bunch)
550     } // while (channel)
551   } // while (ddl)
552   return kTRUE;
553 }
554
555 //____________________________________________________________________
556 Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate, 
557                                      AliFMDCalibStripRange* stripRange, 
558                                      TArrayS &pulseSize, 
559                                      TArrayS &pulseLength, 
560                                      Bool_t* detectors) 
561 {
562   // 
563   // Read SOD event into passed objects.
564   // 
565   // Parameters:
566   //    samplerate   The sample rate object to fill
567   //    striprange   The strip range object to fill
568   //    pulseSize    The pulse size object to fill
569   //    pulseLength  The pulse length (in events) object to fill
570   // 
571   // Return:
572   //    @c true on success
573   //  
574   AliFMDDebug(0, ("Start of SOD/EOD"));
575   
576   UInt_t shift_clk[18];
577   UInt_t sample_clk[18];
578   UInt_t strip_low[18];
579   UInt_t strip_high[18];
580   UInt_t pulse_size[18];
581   UInt_t pulse_length[18];  
582   for (size_t i = 0; i < 18; i++) { 
583     shift_clk[i]    = 0;
584     sample_clk[i]   = 0;
585     strip_low[i]    = 0;
586     strip_high[i]   = 0;
587     pulse_size[i]   = 0;
588     pulse_length[i] = 0;
589   }
590   AliFMDParameters*   param = AliFMDParameters::Instance();
591   AliFMDAltroMapping* map   = param->GetAltroMap();
592   
593   AliAltroRawStreamV3  streamer(fReader);
594   streamer.Reset();
595   streamer.SelectRawData("FMD");
596   //fReader->GetDDLID();
597   while (streamer.NextDDL()) {
598     Int_t ddl   = streamer.GetDDLNumber();
599     Int_t detID = fReader->GetDetectorID();
600     if (detectors) detectors[map->DDL2Detector(ddl)-1] = kTRUE;
601     AliFMDDebug(0, (" From reader: DDL number is %d , det ID is %d",ddl,detID));
602     
603     ULong_t  nPayloadWords = streamer.GetRCUPayloadSizeInSOD();
604     UChar_t* payloadData   = streamer.GetRCUPayloadInSOD();
605     UInt_t*  payloadWords  = reinterpret_cast<UInt_t*>(payloadData);
606     //UInt_t*   payloadWords  = streamer.GetRCUPayloadInSOD();
607
608     //std::cout<<nPayloadWords<<"    "<<ddl<<std::endl;
609     for (ULong_t i = 1; i <= nPayloadWords ; i++, payloadWords++) {
610       UInt_t payloadWord = *payloadWords; // Get32bitWord(i);
611     
612       //std::cout<<i<<Form("  word: 0x%x",payloadWord)<<std::endl;
613       // address is only 24 bit
614       UInt_t address       = (0xffffff & payloadWord);
615       UInt_t type          = ((address >> 21) & 0xf);
616       UInt_t error         = ((address >> 20) & 0x1);
617       UInt_t bcast         = ((address >> 18) & 0x1);
618       UInt_t bc_not_altro  = ((address >> 17) & 0x1);
619       UInt_t board         = ((address >> 12) & 0x1f);
620       UInt_t instruction   = 0;
621       UInt_t chip          = 0;
622       UInt_t channel       = 0;
623       if(bc_not_altro)
624         instruction        = address & 0xfff;
625       else {
626         chip               = ((address >> 9) & 0x7);
627         channel            = ((address >> 5) & 0x5);
628         instruction        = (address & 0x1f);
629       }
630         
631       Bool_t readDataWord = kFALSE;
632       switch(type) {
633       case 0x0: // Fec read
634         readDataWord = kTRUE;  
635       case 0x1: // Fec cmd
636       case 0x2: // Fec write
637         i++;  
638         payloadWords++;
639         break;
640       case 0x4: // Loop
641       case 0x5: // Wait
642         break;
643       case 0x6: // End sequence
644       case 0x7: // End Mem
645         i = nPayloadWords + 1;
646         break;
647       default:    
648         break;
649       }
650         
651       //Don't read unless we have a FEC_RD
652       if(!readDataWord)  continue;
653
654       UInt_t dataWord      = *payloadWords;//Get32bitWord(i);
655       UInt_t data          = (0xFFFFF & dataWord) ;
656       //UInt_t data          = (0xFFFF & dataWord) ;
657         
658       if(error) {
659         AliWarning(Form("error bit detected at Word 0x%06x; "
660                         "error % d, type %d, bc_not_altro %d, "
661                         "bcast %d, board 0x%02x, chip 0x%x, "
662                         "channel 0x%02x, instruction 0x%03x",
663                         address, error, type, bc_not_altro, 
664                         bcast,board,chip,channel,instruction));
665         //process error
666         continue;
667       }
668         
669         
670       switch(instruction) {
671           
672       case 0x01: break;  // First ADC T           
673       case 0x02: break; // I  3.3 V              
674       case 0x03: break; // I  2.5 V altro digital
675       case 0x04: break; // I  2.5 V altro analog 
676       case 0x05: break; // I  2.5 V VA           
677       case 0x06: break; // First ADC T           
678       case 0x07: break; // I  3.3 V              
679       case 0x08: break; // I  2.5 V altro digital
680       case 0x09: break; // I  2.5 V altro analog 
681       case 0x0A: break; // I  2.5 V VA           
682       case 0x2D: break; // Second ADC T           
683       case 0x2E: break; // I  1.5 V VA            
684       case 0x2F: break; // I -2.0 V               
685       case 0x30: break; // I -2.0 V VA            
686       case 0x31: break; //    2.5 V Digital driver
687       case 0x32: break; // Second ADC T           
688       case 0x33: break; // I  1.5 V VA            
689       case 0x34: break; // I -2.0 V               
690       case 0x35: break; // I -2.0 V VA            
691       case 0x36: break; //    2.5 V Digital driver
692       case 0x37: break; // Third ADC T             
693       case 0x38: break; // Temperature sens. 1     
694       case 0x39: break; // Temperature sens. 2     
695       case 0x3A: break; // U  2.5 altro digital (m)
696       case 0x3B: break; // U  2.5 altro analog (m) 
697       case 0x3C: break; // Third ADC T             
698       case 0x3D: break; // Temperature sens. 1     
699       case 0x3E: break; // Temperature sens. 2     
700       case 0x3F: break; // U  2.5 altro digital (m)
701       case 0x40: break; // U  2.5 altro analog (m) 
702       case 0x41: break; // Forth ADC T  
703       case 0x42: break; // U  2.5 VA (m)
704       case 0x43: break; // U  1.5 VA (m)
705       case 0x44: break; // U -2.0 VA (m)
706       case 0x45: break; // U -2.0 (m)   
707       case 0x46: break; // Forth ADC T  
708       case 0x47: break; // U  2.5 VA (m)
709       case 0x48: break; // U  1.5 VA (m)
710       case 0x49: break; // U -2.0 VA (m)
711       case 0x4A: break;  // U -2.0 (m)   
712         // Counters 
713       case 0x0B: break; // L1 trigger CouNTer
714       case 0x0C: break; // L2 trigger CouNTer
715       case 0x0D: break; // Sampling CLK CouNTer
716       case 0x0E: break; // DSTB CouNTer
717         // Test mode 
718       case 0x0F: break; // Test mode word
719       case 0x10: break; // Undersampling ratio.
720         // Configuration and status 
721       case 0x11: break; // Config/Status Register 0
722       case 0x12: break; // Config/Status Register 1
723       case 0x13: break; // Config/Status Register 2
724       case 0x14: break; // Config/Status Register 3
725       case 0x15: break; // Free
726         // Comands:
727       case 0x16: break; // Latch L1, L2, SCLK Counters
728       case 0x17: break; // Clear counters
729       case 0x18: break; // Clear CSR1
730       case 0x19: break; // rstb ALTROs
731       case 0x1A: break; // rstb BC
732       case 0x1B: break; // Start conversion
733       case 0x1C: break; // Scan event length
734       case 0x1D: break; // Read event length
735       case 0x1E: break; // Start test mode
736       case 0x1F: break; // Read acquisition memory
737         // FMD
738       case 0x20: break; // FMDD status
739       case 0x21: break; // L0 counters
740       case 0x22: break; // FMD: Wait to hold
741       case 0x23: break; // FMD: L1 timeout
742       case 0x24: break; // FMD: L2 timeout
743       case 0x25: // FMD: Shift clk 
744         shift_clk[board] = ((data >> 8 ) & 0xFF); 
745         AliFMDDebug(30, ("Read shift_clk=%d for board 0x%02x", 
746                          shift_clk[board], board));
747         break; 
748       case 0x26: // FMD: Strips 
749         strip_low[board]  = ((data >> 0 ) & 0xFF); 
750         strip_high[board] = ((data >> 8 ) & 0xFF);  
751         break; 
752       case 0x27: // FMD: Cal pulse 
753         pulse_size[board]  =  ((data >> 8 ) & 0xFF);
754         break; 
755       case 0x28: break; // FMD: Shape bias
756       case 0x29: break; // FMD: Shape ref
757       case 0x2A: break; // FMD: Preamp ref
758       case 0x2B: // FMD: Sample clk 
759         sample_clk[board] = ((data >> 8 ) & 0xFF); 
760         AliFMDDebug(30, ("Read sample_clk=%d for board 0x%02x", 
761                          sample_clk[board], board));
762         break; 
763       case 0x2C: break; // FMD: Commands
764       case 0x4B: // FMD: Cal events 
765         pulse_length[board] = ((data >> 0 ) & 0xFF);
766         break; 
767       default: break;
768         
769       }
770       AliFMDDebug(50, ("instruction 0x%x, dataword 0x%x",
771                        instruction,dataWord));
772     } // End of loop over Result memory event
773     
774     UShort_t det    = 0;
775     UShort_t sector = 0;
776     Short_t  strip  = -1;
777     Char_t   ring   = '\0';
778    
779     const UInt_t boards[4] = {0,1,16,17};
780     for(Int_t i=0;i<4;i++) {
781       if(ddl==0 && (i==1 || i==3)) continue;
782
783       UInt_t chip =0, channel=0;
784       det = map->DDL2Detector(ddl);
785       map->Channel2StripBase(boards[i], chip, channel, ring, sector, strip);
786      
787       UInt_t samplerate = 0;
788 #if USE_VOTE
789       if(sample_clk[boards[i]] == 0) {
790         if(ddl == 0) {
791           Int_t sample1 = sample_clk[boards[0]];
792           Int_t sample2 = sample_clk[boards[2]];            
793           if(sample1) sample_clk[boards[i]] = sample1;
794           else sample_clk[boards[i]] = sample2;
795         }
796         else {
797           Int_t sample1 = sample_clk[boards[0]];
798           Int_t sample2 = sample_clk[boards[1]];
799           Int_t sample3 = sample_clk[boards[2]];
800           Int_t sample4 = sample_clk[boards[3]];
801           Int_t agreement = 0;
802           if(sample1 == sample2) agreement++;
803           if(sample1 == sample3) agreement++;
804           if(sample1 == sample4) agreement++;
805           if(sample2 == sample3) agreement++;
806           if(sample2 == sample4) agreement++;
807           if(sample3 == sample4) agreement++;
808             
809           Int_t idx = 0;
810           if(i<3) idx = i+1;
811           else  idx = i-1;
812           if(agreement == 3) {
813             sample_clk[boards[i]] = sample_clk[boards[idx]];
814             shift_clk[boards[i]] = shift_clk[boards[idx]];
815             strip_low[boards[i]] = strip_low[boards[idx]];
816             strip_high[boards[i]] = strip_high[boards[idx]];
817             pulse_length[boards[i]] = pulse_length[boards[idx]];
818             pulse_size[boards[i]] = pulse_size[boards[idx]];
819             AliFMDDebug(3, ("Vote taken for ddl %d, board 0x%x",
820                             ddl,boards[i]));
821           }
822         }
823       } 
824 #endif
825       
826       if(sample_clk[boards[i]])
827         samplerate = shift_clk[boards[i]]/sample_clk[boards[i]];
828       AliFMDDebug(10, ("Sample rate for board 0x%02x is %d", 
829                       boards[i], samplerate));
830       sampleRate->Set(det,ring,sector,0,samplerate);
831       stripRange->Set(det,ring,sector,0,
832                       strip_low[boards[i]],strip_high[boards[i]]);
833       
834       AliFMDDebug(20, ("det %d, ring %c, ",det,ring));
835       pulseLength.AddAt(pulse_length[boards[i]],
836                         GetHalfringIndex(det,ring,boards[i]/16));
837       pulseSize.AddAt(pulse_size[boards[i]],
838                       GetHalfringIndex(det,ring,boards[i]/16));
839       
840       
841       
842       AliFMDDebug(20, (": Board: 0x%02x\n"
843                        "\tstrip_low  %3d, strip_high   %3d\n"
844                        "\tshift_clk  %3d, sample_clk   %3d\n"
845                        "\tpulse_size %3d, pulse_length %3d",
846                        boards[i], 
847                        strip_low[boards[i]], strip_high[boards[i]],
848                        shift_clk[boards[i]], sample_clk[boards[i]],
849                        pulse_size[boards[i]],pulse_length[boards[i]]));
850     }
851     
852   }
853   
854   AliFMDParameters::Instance()->SetSampleRate(sampleRate);
855   AliFMDParameters::Instance()->SetStripRange(stripRange);
856   
857   AliFMDDebug(0, ("End of SOD/EOD"));
858   
859   return kTRUE;
860 }
861 //____________________________________________________________________
862
863 UInt_t AliFMDRawReader::Get32bitWord(Int_t idx)
864 {
865   // This method returns the 32 bit word at a given
866   // position inside the raw data payload.
867   // The 'index' points to the beginning of the next word.
868   // The method is supposed to be endian (platform)
869   // independent.
870   if (!fData) {
871     AliFatal("Raw data paylod buffer is not yet initialized !");
872   }
873
874   Int_t index = 4*idx;
875   
876   if (index < 4) {
877     //  fRawReader->AddFatalErrorLog(k32bitWordReadErr,Form("pos = %d",index));
878     //    PrintDebug();
879     AliWarning(Form("Invalid raw data payload position (%d) !",index));
880   }
881
882   UInt_t word = 0;
883    
884   word  = fData[--index] << 24;
885   word |= fData[--index] << 16;
886   word |= fData[--index] << 8;
887   word |= fData[--index] << 0 ;
888
889   return word;
890 }
891 //_____________________________________________________________________ 
892 Int_t AliFMDRawReader::GetHalfringIndex(UShort_t det, Char_t ring, 
893                                         UShort_t board) const
894 {
895   // 
896   // Get short index for a given half-ring
897   // 
898   // Parameters:
899   //    det   Detector number
900   //    ring  Ring identifer
901   //    board Board number
902   // 
903   // Return:
904   //    
905   //
906   UShort_t iring  =  (ring == 'I' ? 1 : 0);
907   
908   Int_t index = (((det-1) << 2) | (iring << 1) | (board << 0));
909   
910   return index-2;
911   
912 }
913
914 //____________________________________________________________________
915 // 
916 // EOF
917 //
918