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