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