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