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 }
 116
 117 //____________________________________________________________________
 118 Int_t
 119 AliFMDRawReader::NewDDL(AliAltroRawStreamV3& input, UShort_t& det)
 120 {
 121   // Process a new DDL.  Sets the internal data members fZeroSuppress,
 122   // fSampleRate, and fNoiseFactor based on information in the RCU trailer.
 123   //
 124   // Parameters:
 125   //   input Input stream
 126   //   det   On return, the detector number
 127   //
 128   // Return value:
 129   //   negative value in case of problems, the DDL number otherwise
 130
 131   // Get the DDL number
 132   UInt_t ddl = input.GetDDLNumber();
 133   AliFMDDebug(2, ("DDL number %d", ddl));
 134
 135   // Get zero suppression flag
 136   fZeroSuppress[ddl] = input.GetZeroSupp();
 137   AliFMDDebug(3, ("RCU @ DDL %d zero suppression: %s",
 138                    ddl, (fZeroSuppress[ddl] ? "yes" : "no")));
 139
 140   // WARNING: We store the noise factor in the 2nd baseline
 141   // filters excluded post samples, since we'll never use that
 142   // mode.
 143   fNoiseFactor[ddl]  = input.GetNPostsamples();
 144   AliFMDDebug(3, ("RCU @ DDL %d noise factor: %d", ddl,fNoiseFactor[ddl]));
 145
 146   // WARNING: We store the sample rate in the number of pre-trigger
 147   // samples, since we'll never use that mode.
 148   fSampleRate[ddl]     = input.GetNPretriggerSamples();
 149   AliFMDDebug(3, ("RCU @ DDL %d sample rate: %d", ddl,fSampleRate[ddl]));
 150
 151   // Get Errors seen
 152   Int_t nChAddrMismatch = input.GetNChAddrMismatch();
 153   Int_t nChLenMismatch  = input.GetNChLengthMismatch();
 154   if (nChAddrMismatch != 0)
 155     AliWarning(Form("Got %d channels with address mis-matches for 0x%03x",
 156                     nChAddrMismatch, ddl));
 157   if (nChLenMismatch != 0)
 158     AliWarning(Form("Got %d channels with length mis-matches for 0x%03x",
 159                     nChLenMismatch, ddl));
 160
 161   // Map DDL number to the detector number
 162   AliFMDParameters*    pars   = AliFMDParameters::Instance();
 163   AliFMDAltroMapping*  map    = pars->GetAltroMap();
 164   if (map->DDL2Detector(ddl) < 0) return -1;
 165   det = map->DDL2Detector(ddl);
 166
 167   if (AliLog::GetDebugLevel("FMD", 0) > 5)
 168     input.PrintRCUTrailer();
 169   return ddl;
 170 }
 171
 172 //____________________________________________________________________
 173 Int_t
 174 AliFMDRawReader::NewChannel(AliAltroRawStreamV3& input,  UShort_t det,
 175                             Char_t&  ring, UShort_t& sec, Short_t& strbase)
 176 {
 177   // Processs a new channel.  Sets the internal data members
 178   // fMinStrip, fMaxStrip, and fPreSamp.
 179   //
 180   // Parameter:
 181   //   input   Input stream
 182   //   ring    On return, the ring identifier
 183   //   sec     On return, the sector number
 184   //   strbase On return, the strip base
 185   //
 186   // Return value
 187   //   negative value in case of problems, hardware address otherwise
 188
 189   // Get the hardware address, and map that to detector coordinates
 190   UShort_t board, chip, channel;
 191   Int_t    ddl    = input.GetDDLNumber();
 192   Int_t    hwaddr = input.GetHWAddress();
 193   if (input.IsChannelBad()) {
 194     AliError(Form("Channel 0x%03x is marked as bad!", hwaddr));
 195   }
 196
 197   AliFMDParameters*    pars   = AliFMDParameters::Instance();
 198   AliFMDAltroMapping*  map    = pars->GetAltroMap();
 199   // Map to hardware stuff
 200   map->ChannelAddress(hwaddr, board, chip, channel);
 201   // Then try to map to detector address
 202   if (!map->Channel2StripBase(board, chip, channel, ring, sec, strbase)) {
 203     AliError(Form("Failed to get detector id from DDL %d, "
 204                   "hardware address 0x%03x", ddl, hwaddr));
 205     return -1;
 206   }
 207   AliFMDDebug(3, ("Board: 0x%02x, Altro: 0x%x, Channel: 0x%x",
 208                   board, chip, channel));
 209
 210   // Get the 'conditions'
 211   fMinStrip = pars->GetMinStrip(det, ring, sec, strbase);
 212   fMaxStrip = pars->GetMaxStrip(det, ring, sec, strbase);
 213   fPreSamp  = pars->GetPreSamples(det, ring, sec, strbase);
 214   if (fSampleRate[ddl] == 0)
 215     fSampleRate[ddl] = pars->GetSampleRate(det, ring, sec, strbase);
 216
 217   return hwaddr;
 218 }
 219
 220 //____________________________________________________________________
 221 Int_t
 222 AliFMDRawReader::NewSample(AliAltroRawStreamV3& input,
 223                            Int_t i, UShort_t t, UShort_t sec,
 224                            UShort_t  strbase, Short_t&  str, UShort_t& samp)
 225 {
 226   // Process a new timebin
 227   //
 228   // Parameters:
 229   //   input   Input stream
 230   //   i       Index into bunch data
 231   //   t       Time
 232   //   strbase Base of strip numbers for this channel
 233   //   str     On return, the strip number
 234   //   samp    On return, the sample number
 235   //
 236   // Return value
 237   //   negative value in case of problems, ADC value otherwise
 238   if (t < fPreSamp) return -1;
 239
 240   Int_t           ddl  = input.GetDDLNumber();
 241   Int_t           hwa  = input.GetHWAddress();
 242   const UShort_t* data = input.GetSignals();
 243   Short_t         adc  = data[i];
 244   AliFMDDebug(10, ("0x%04x/0x%03x/%04d %4d", ddl, hwa, t, adc));
 245
 246   AliFMDParameters*    pars   = AliFMDParameters::Instance();
 247   AliFMDAltroMapping*  map    = pars->GetAltroMap();
 248
 249   samp            = 0;
 250   Short_t  stroff = 0;
 251   map->Timebin2Strip(sec, t, fPreSamp, fSampleRate[ddl], stroff, samp);
 252   str             = strbase + stroff;
 253
 254   AliFMDDebug(20, ("0x%04x/0x%03x/%04d maps to strip %3d sample %d "
 255                    "(pre: %d, min: %d, max: %d, rate: %d)",
 256                   ddl, hwa, t, str, samp, fPreSamp,
 257                   fMinStrip, fMaxStrip, fSampleRate[ddl]));
 258   if (str < 0) {
 259     AliFMDDebug(10, ("Got presamples at timebin %d", i));
 260     return -1;
 261   }
 262
 263   // VA1 Local strip number
 264   Short_t lstrip = (t - fPreSamp) / fSampleRate[ddl] + fMinStrip;
 265
 266   AliFMDDebug(15, ("Checking if strip %d (%d) in range [%d,%d]",
 267                    lstrip, str, fMinStrip, fMaxStrip));
 268   if (lstrip < fMinStrip || lstrip > fMaxStrip) {
 269     AliFMDDebug(20, ("Strip %03d-%d (%d,%d) from t=%d out of range (%3d->%3d)",
 270                     str, samp, lstrip, stroff, t, fMinStrip, fMaxStrip));
 271     return -1;
 272   }
 273   // Possibly do pedestal subtraction of signal
 274   if (adc > 1023)
 275     AliWarning(Form("ADC value out of range: %4d", adc));
 276   return adc;
 277 }
 278
 279 //____________________________________________________________________
 280 Bool_t
 281 AliFMDRawReader::NextSample(UShort_t& det, Char_t&   rng, UShort_t& sec,
 282                             UShort_t& str, UShort_t& sam, UShort_t& rat,
 283                             Short_t&  adc, Bool_t&   zs,  UShort_t& fac)
 284 {
 285   // Scan current event for next signal.   It returns kFALSE when
 286   // there's no more data in the event.
 287   static AliAltroRawStreamV3 stream(fReader); //    = 0;
 288   static Int_t               ddl      = -1;
 289   static UShort_t            tdet     = 0;
 290   static Char_t              trng     = '\0';
 291   static UShort_t            tsec     = 0;
 292   static Short_t             tstr     = 0;
 293   static Short_t             bstr     = -1;
 294   static UShort_t            tsam     = 0;
 295   static UInt_t              trate    = 0;
 296   static Int_t               hwaddr   = -1;
 297   static UShort_t            start    = 0;
 298   static UShort_t            length   = 0;
 299   static Short_t             t        = -1;
 300   static Int_t               i        = 0;
 301   // First entry!
 302   if (stream.GetDDLNumber() < 0) {
 303     fReader->Select("FMD");
 304
 305     // Reset variables
 306     ddl    = -1;
 307     trate  = 0;
 308     tdet   = 0;
 309     trng   = '\0';
 310     tsec   = 0;
 311     tstr   = 0;
 312     tsam   = -1;
 313     hwaddr = -1;
 314   }
 315
 316   do {
 317     if (t < start - length + 1) {
 318       if (!stream.NextBunch()) {
 319         if (!stream.NextChannel()) {
 320           if (!stream.NextDDL()) {
 321             stream.Reset();
 322             return kFALSE;
 323           }
 324           ddl = NewDDL(stream, tdet);
 325         }
 326         hwaddr = NewChannel(stream, tdet, trng, tsec, bstr);
 327       }
 328       start  = stream.GetStartTimeBin();
 329       length = stream.GetBunchLength();
 330       t      = start;
 331       i      = 0;
 332     }
 333     Int_t tadc = NewSample(stream, i, t, tsec, bstr, tstr, tsam);
 334     if (tadc >= 0) {
 335       det = tdet;
 336       rng = trng;
 337       sec = tsec;
 338       str = tstr;
 339       sam = tsam;
 340       adc = tadc;
 341       rat = fSampleRate[ddl];
 342       zs  = fZeroSuppress[ddl];
 343       fac = fNoiseFactor[ddl];
 344       t--;
 345       i++;
 346       break;
 347     }
 348   } while (true);
 349
 350   return kTRUE;
 351 }
 352
 353
 354 //____________________________________________________________________
 355 Bool_t
 356 AliFMDRawReader::NextSignal(UShort_t& det, Char_t&   rng,
 357                             UShort_t& sec, UShort_t& str,
 358                             Short_t&  adc, Bool_t&   zs,
 359                             UShort_t& fac)
 360 {
 361
 362   do {
 363     UShort_t samp, rate;
 364     if (!NextSample(det, rng, sec, str, samp, rate, adc, zs, fac))
 365       return kFALSE;
 366
 367     Bool_t take = kFALSE;
 368     switch (rate) {
 369     case 1:                      take = kTRUE; break;
 370     case 2:  if (samp == 1)      take = kTRUE; break;
 371     case 3:  if (samp == 1)      take = kTRUE; break;
 372     case 4:  if (samp == 2)      take = kTRUE; break;
 373     default: if (samp == rate-2) take = kTRUE; break;
 374     }
 375     if (!take) continue;
 376     break;
 377   } while (true);
 378   return kTRUE;
 379 }
 380
 381 //____________________________________________________________________
 382 Bool_t
 383 AliFMDRawReader::SelectSample(UShort_t samp, UShort_t rate)
 384 {
 385   // Check if the passed sample is the one we need
 386   Bool_t take = kFALSE;
 387   switch (rate) {
 388   case 1:                      take = kTRUE; break;
 389   case 2:  if (samp == 1)      take = kTRUE; break;
 390   case 3:  if (samp == 1)      take = kTRUE; break;
 391   case 4:  if (samp == 2)      take = kTRUE; break;
 392   default: if (samp == rate-2) take = kTRUE; break;
 393   }
 394
 395   return take;
 396 }
 397
 398 //____________________________________________________________________
 399 Bool_t
 400 AliFMDRawReader::ReadAdcs(TClonesArray* array)
 401 {
 402   // Read ADC values from raw input into passed TClonesArray of AliFMDDigit
 403   // objects.
 404   AliFMDDebug(2, ("Reading ADC values into a TClonesArray"));
 405
 406   // Read raw data into the digits array, using AliFMDAltroReader.
 407   if (!array) {
 408     AliError("No TClonesArray passed");
 409     return kFALSE;
 410   }
 411   const UShort_t kUShortMax = (1 << 16) - 1;
 412   fSeen.Reset(kUShortMax);
 413
 414   AliAltroRawStreamV3  input(fReader);
 415   input.Reset();
 416   input.SelectRawData("FMD");
 417
 418   // Loop over input RORCs
 419   while (input.NextDDL()) {
 420     UShort_t det = 0;
 421     Int_t    ddl = NewDDL(input, det);
 422     if (ddl < 0) break;
 423
 424     while (input.NextChannel()) {
 425       // Get the hardware address, and map that to detector coordinates
 426       Char_t   ring;
 427       UShort_t sec;
 428       Short_t  strbase;
 429       Int_t    hwaddr = NewChannel(input, det, ring, sec, strbase);
 430       if (hwaddr < 0) break;
 431
 432       // Loop over bunches
 433       while (input.NextBunch()) {
 434         // Get Lenght of bunch, and pointer to the data
 435         const UShort_t* data   = input.GetSignals();
 436         UShort_t        start  = input.GetStartTimeBin();
 437         UShort_t        length = input.GetBunchLength();
 438
 439         // Loop over the data and store it.
 440         for (Int_t i = 0; i < length; i++) {
 441           // Time
 442           Short_t  str;
 443           UShort_t samp;
 444           Int_t    t    = start - i;
 445           Int_t    adc  = NewSample(input, i, t, sec, strbase, str, samp);
 446           if (adc < 0) continue;
 447           UShort_t counts = adc;
 448
 449           AliFMDDebug(10, ("FMD%d%c[%02d,%03d]-%d: %4d",
 450                            det, ring, sec, str, samp, counts));
 451           // Check the cache of indicies
 452           Int_t idx = fSeen(det, ring, sec, str);
 453           if (idx == kUShortMax) {
 454             // We haven't seen this strip yet.
 455             fSeen(det, ring, sec, str) = idx = array->GetEntriesFast();
 456             AliFMDDebug(7,("making digit for FMD%d%c[%2d,%3d]-%d "
 457                            "from timebin %4d",
 458                            det, ring, sec, str, samp, t));
 459             new ((*array)[idx]) AliFMDDigit(det, ring, sec, str);
 460           }
 461           AliFMDBaseDigit* digit =
 462             static_cast<AliFMDBaseDigit*>(array->At(idx));
 463           AliFMDDebug(10, ("Setting FMD%d%c[%2d,%3d]-%d "
 464                            "from timebin %4d=%4d (%4d)",
 465                            det, ring, sec, str, samp, t, counts, data[i]));
 466           digit->SetCount(samp, counts);
 467         } // for (i)
 468       } // while (bunch)
 469     } // while (channel)
 470   } // while (ddl)
 471   return kTRUE;
 472 }
 473
 474 //____________________________________________________________________
 475 Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate,
 476                                      AliFMDCalibStripRange* stripRange,
 477                                      TArrayS &pulseSize,
 478                                      TArrayS &pulseLength,
 479                                      Bool_t* detectors)
 480 {
 481
 482   AliFMDDebug(0, ("Start of SOD/EOD"));
 483
 484   UInt_t shift_clk[18];
 485   UInt_t sample_clk[18];
 486   UInt_t strip_low[18];
 487   UInt_t strip_high[18];
 488   UInt_t pulse_size[18];
 489   UInt_t pulse_length[18];
 490   AliFMDParameters*   param = AliFMDParameters::Instance();
 491   AliFMDAltroMapping* map   = param->GetAltroMap();
 492
 493   while(fReader->ReadNextData(fData)) {
 494
 495     Int_t ddl   = fReader->GetDDLID();
 496     Int_t detID = fReader->GetDetectorID();
 497     if (detectors) detectors[map->DDL2Detector(ddl)-1] = kTRUE;
 498     AliFMDDebug(0, (" From reader: DDL number is %d , det ID is %d",ddl,detID));
 499
 500     fNbytes = fReader->GetDataSize();
 501
 502     ULong_t nWords      = GetNwords();
 503     UInt_t  trailerLast = Get32bitWord(nWords);
 504     AliFMDDebug(20, (" # Bytes: %d, # Words: %d, Last word: 0x%08x",
 505                      fNbytes, nWords, trailerLast));
 506     if ((trailerLast & 0xFFFF0000) != 0xAAAA0000) {
 507       AliWarning(Form("Last word 0x%08x does not match RCU II trailer",
 508                       trailerLast));
 509       //return kFALSE;
 510     }
 511
 512     // 7 last bits is size of trailer
 513     ULong_t nTrailerWords = trailerLast & 0x7f;
 514     ULong_t nPayloadWords = Get32bitWord(nWords - nTrailerWords+1);
 515
 516     AliFMDDebug(20, (" # trailer words: %d, # payload words: %d",
 517                      nTrailerWords, nPayloadWords));
 518
 519
 520     for (ULong_t i = 1; i <= nPayloadWords ; i++) {
 521       UInt_t payloadWord = Get32bitWord(i);
 522
 523       // address is only 24 bit
 524       UInt_t address       = (0xffffff & payloadWord);
 525       UInt_t type          = ((address >> 21) & 0xf);
 526       UInt_t error         = ((address >> 20) & 0x1);
 527       UInt_t bcast         = ((address >> 18) & 0x1);
 528       UInt_t bc_not_altro  = ((address >> 17) & 0x1);
 529       UInt_t board         = ((address >> 12) & 0x1f);
 530       UInt_t instruction   = 0;
 531       UInt_t chip          = 0;
 532       UInt_t channel       = 0;
 533       if(bc_not_altro)
 534         instruction        = address & 0xfff;
 535       else {
 536         chip               = ((address >> 9) & 0x7);
 537         channel            = ((address >> 5) & 0x5);
 538         instruction        = (address & 0x1f);
 539       }
 540
 541       Bool_t readDataWord = kFALSE;
 542       switch(type) {
 543       case 0x0: // Fec read
 544         readDataWord = kTRUE;
 545       case 0x1: // Fec cmd
 546       case 0x2: // Fec write
 547         i++;
 548         break;
 549       case 0x4: // Loop
 550       case 0x5: // Wait
 551         break;
 552       case 0x6: // End sequence
 553       case 0x7: // End Mem
 554         i = nPayloadWords + 1;
 555         break;
 556       default:
 557         break;
 558       }
 559
 560       //Don't read unless we have a FEC_RD
 561       if(!readDataWord)  continue;
 562
 563       UInt_t dataWord      = Get32bitWord(i);
 564       UInt_t data          = (0xFFFFF & dataWord) ;
 565       //UInt_t data          = (0xFFFF & dataWord) ;
 566
 567       if(error) {
 568         AliWarning(Form("error bit detected at Word 0x%06x; "
 569                         "error % d, type %d, bc_not_altro %d, "
 570                         "bcast %d, board 0x%02x, chip 0x%x, "
 571                         "channel 0x%02x, instruction 0x%03x",
 572                         address, error, type, bc_not_altro,
 573                         bcast,board,chip,channel,instruction));
 574         //process error
 575         continue;
 576       }
 577
 578
 579       switch(instruction) {
 580
 581       case 0x01: break;  // First ADC T
 582       case 0x02: break; // I  3.3 V
 583       case 0x03: break; // I  2.5 V altro digital
 584       case 0x04: break; // I  2.5 V altro analog
 585       case 0x05: break; // I  2.5 V VA
 586       case 0x06: break; // First ADC T
 587       case 0x07: break; // I  3.3 V
 588       case 0x08: break; // I  2.5 V altro digital
 589       case 0x09: break; // I  2.5 V altro analog
 590       case 0x0A: break; // I  2.5 V VA
 591       case 0x2D: break; // Second ADC T
 592       case 0x2E: break; // I  1.5 V VA
 593       case 0x2F: break; // I -2.0 V
 594       case 0x30: break; // I -2.0 V VA
 595       case 0x31: break; //    2.5 V Digital driver
 596       case 0x32: break; // Second ADC T
 597       case 0x33: break; // I  1.5 V VA
 598       case 0x34: break; // I -2.0 V
 599       case 0x35: break; // I -2.0 V VA
 600       case 0x36: break; //    2.5 V Digital driver
 601       case 0x37: break; // Third ADC T
 602       case 0x38: break; // Temperature sens. 1
 603       case 0x39: break; // Temperature sens. 2
 604       case 0x3A: break; // U  2.5 altro digital (m)
 605       case 0x3B: break; // U  2.5 altro analog (m)
 606       case 0x3C: break; // Third ADC T
 607       case 0x3D: break; // Temperature sens. 1
 608       case 0x3E: break; // Temperature sens. 2
 609       case 0x3F: break; // U  2.5 altro digital (m)
 610       case 0x40: break; // U  2.5 altro analog (m)
 611       case 0x41: break; // Forth ADC T
 612       case 0x42: break; // U  2.5 VA (m)
 613       case 0x43: break; // U  1.5 VA (m)
 614       case 0x44: break; // U -2.0 VA (m)
 615       case 0x45: break; // U -2.0 (m)
 616       case 0x46: break; // Forth ADC T
 617       case 0x47: break; // U  2.5 VA (m)
 618       case 0x48: break; // U  1.5 VA (m)
 619       case 0x49: break; // U -2.0 VA (m)
 620       case 0x4A: break;  // U -2.0 (m)
 621         // Counters
 622       case 0x0B: break; // L1 trigger CouNTer
 623       case 0x0C: break; // L2 trigger CouNTer
 624       case 0x0D: break; // Sampling CLK CouNTer
 625       case 0x0E: break; // DSTB CouNTer
 626         // Test mode
 627       case 0x0F: break; // Test mode word
 628       case 0x10: break; // Undersampling ratio.
 629         // Configuration and status
 630       case 0x11: break; // Config/Status Register 0
 631       case 0x12: break; // Config/Status Register 1
 632       case 0x13: break; // Config/Status Register 2
 633       case 0x14: break; // Config/Status Register 3
 634       case 0x15: break; // Free
 635         // Comands:
 636       case 0x16: break; // Latch L1, L2, SCLK Counters
 637       case 0x17: break; // Clear counters
 638       case 0x18: break; // Clear CSR1
 639       case 0x19: break; // rstb ALTROs
 640       case 0x1A: break; // rstb BC
 641       case 0x1B: break; // Start conversion
 642       case 0x1C: break; // Scan event length
 643       case 0x1D: break; // Read event length
 644       case 0x1E: break; // Start test mode
 645       case 0x1F: break; // Read acquisition memory
 646         // FMD
 647       case 0x20: break; // FMDD status
 648       case 0x21: break; // L0 counters
 649       case 0x22: break; // FMD: Wait to hold
 650       case 0x23: break; // FMD: L1 timeout
 651       case 0x24: break; // FMD: L2 timeout
 652       case 0x25: // FMD: Shift clk
 653         shift_clk[board] = ((data >> 8 ) & 0xFF);
 654         break;
 655       case 0x26: // FMD: Strips
 656         strip_low[board]  = ((data >> 0 ) & 0xFF);
 657         strip_high[board] = ((data >> 8 ) & 0xFF);
 658         break;
 659       case 0x27: // FMD: Cal pulse
 660         pulse_size[board]  =  ((data >> 8 ) & 0xFF);
 661         break;
 662       case 0x28: break; // FMD: Shape bias
 663       case 0x29: break; // FMD: Shape ref
 664       case 0x2A: break; // FMD: Preamp ref
 665       case 0x2B: // FMD: Sample clk
 666         sample_clk[board] = ((data >> 8 ) & 0xFF);
 667         break;
 668       case 0x2C: break; // FMD: Commands
 669       case 0x4B: // FMD: Cal events
 670         pulse_length[board] = ((data >> 0 ) & 0xFF);
 671         break;
 672       default: break;
 673
 674       }
 675       AliFMDDebug(50, ("instruction 0x%x, dataword 0x%x",
 676                        instruction,dataWord));
 677     }
 678
 679     UShort_t det,sector;
 680     Short_t strip;
 681     Char_t ring;
 682
 683     const UInt_t boards[4] = {0,1,16,17};
 684     for(Int_t i=0;i<4;i++) {
 685       if(ddl==0 && (i==1 || i==3)) continue;
 686
 687       UInt_t chip =0, channel=0;
 688       det = map->DDL2Detector(ddl);
 689       map->Channel2StripBase(boards[i], chip, channel, ring, sector, strip);
 690
 691       UInt_t samplerate = 1;
 692       if(sample_clk[boards[i]] == 0) {
 693         if(ddl == 0) {
 694           Int_t sample1 = sample_clk[boards[0]];
 695           Int_t sample2 = sample_clk[boards[2]];
 696           if(sample1) sample_clk[boards[i]] = sample1;
 697           else sample_clk[boards[i]] = sample2;
 698         }
 699         else {
 700           Int_t sample1 = sample_clk[boards[0]];
 701           Int_t sample2 = sample_clk[boards[1]];
 702           Int_t sample3 = sample_clk[boards[2]];
 703           Int_t sample4 = sample_clk[boards[3]];
 704           Int_t agreement = 0;
 705           if(sample1 == sample2) agreement++;
 706           if(sample1 == sample3) agreement++;
 707           if(sample1 == sample4) agreement++;
 708           if(sample2 == sample3) agreement++;
 709           if(sample2 == sample4) agreement++;
 710           if(sample3 == sample4) agreement++;
 711
 712           Int_t idx = 0;
 713           if(i<3) idx = i+1;
 714           else  idx = i-1;
 715           if(agreement == 3) {
 716             sample_clk[boards[i]] = sample_clk[boards[idx]];
 717             shift_clk[boards[i]] = shift_clk[boards[idx]];
 718             strip_low[boards[i]] = strip_low[boards[idx]];
 719             strip_high[boards[i]] = strip_high[boards[idx]];
 720             pulse_length[boards[i]] = pulse_length[boards[idx]];
 721             pulse_size[boards[i]] = pulse_size[boards[idx]];
 722             AliFMDDebug(0, ("Vote taken for ddl %d, board 0x%x",
 723                             ddl,boards[i]));
 724           }
 725         }
 726       }
 727
 728       if(sample_clk[boards[i]])
 729         samplerate = shift_clk[boards[i]]/sample_clk[boards[i]];
 730       sampleRate->Set(det,ring,sector,0,samplerate);
 731       stripRange->Set(det,ring,sector,0,
 732                       strip_low[boards[i]],strip_high[boards[i]]);
 733
 734       AliFMDDebug(20, ("det %d, ring %c, ",det,ring));
 735       pulseLength.AddAt(pulse_length[boards[i]],
 736                         GetHalfringIndex(det,ring,boards[i]/16));
 737       pulseSize.AddAt(pulse_size[boards[i]],
 738                       GetHalfringIndex(det,ring,boards[i]/16));
 739
 740
 741
 742       AliFMDDebug(20, (": Board: 0x%02x\n"
 743                        "\tstrip_low  %3d, strip_high   %3d\n"
 744                        "\tshift_clk  %3d, sample_clk   %3d\n"
 745                        "\tpulse_size %3d, pulse_length %3d",
 746                        boards[i],
 747                        strip_low[boards[i]], strip_high[boards[i]],
 748                        shift_clk[boards[i]], sample_clk[boards[i]],
 749                        pulse_size[boards[i]],pulse_length[boards[i]]));
 750     }
 751
 752   }
 753
 754   AliFMDParameters::Instance()->SetSampleRate(sampleRate);
 755   AliFMDParameters::Instance()->SetStripRange(stripRange);
 756
 757   AliFMDDebug(0, ("End of SOD/EOD"));
 758
 759   return kTRUE;
 760 }
 761 //____________________________________________________________________
 762
 763 UInt_t AliFMDRawReader::Get32bitWord(Int_t idx)
 764 {
 765   // This method returns the 32 bit word at a given
 766   // position inside the raw data payload.
 767   // The 'index' points to the beginning of the next word.
 768   // The method is supposed to be endian (platform)
 769   // independent.
 770   if (!fData) {
 771     AliFatal("Raw data paylod buffer is not yet initialized !");
 772   }
 773
 774   Int_t index = 4*idx;
 775
 776   if (index < 4) {
 777     //  fRawReader->AddFatalErrorLog(k32bitWordReadErr,Form("pos = %d",index));
 778     //    PrintDebug();
 779     AliWarning(Form("Invalid raw data payload position (%d) !",index));
 780   }
 781
 782   UInt_t word = 0;
 783
 784   word  = fData[--index] << 24;
 785   word |= fData[--index] << 16;
 786   word |= fData[--index] << 8;
 787   word |= fData[--index] << 0 ;
 788
 789   return word;
 790 }
 791 //_____________________________________________________________________
 792 Int_t AliFMDRawReader::GetHalfringIndex(UShort_t det, Char_t ring,
 793                                         UShort_t board) {
 794
 795   UShort_t iring  =  (ring == 'I' ? 1 : 0);
 796
 797   Int_t index = (((det-1) << 2) | (iring << 1) | (board << 0));
 798
 799   return index-2;
 800
 801 }
 802
 803 //____________________________________________________________________
 804 //
 805 // EOF
 806 //
 807