]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - FMD/AliFMDRawReader.cxx
Improved ALTRO mapping: Now does the full conversion from
[u/mrichter/AliRoot.git] / FMD / AliFMDRawReader.cxx
index fba982b38ccbe345ce9e4edb2221b6563e06cfa5..bc4c14255e4c2397f30dac16bb8f110c8b6a224e 100644 (file)
@@ -77,7 +77,8 @@ AliFMDRawReader::AliFMDRawReader(AliRawReader* reader, TTree* tree)
     fReader(reader), 
     fSampleRate(1),
     fData(0),
-    fNbytes(0)
+    fNbytes(0), 
+    fSeen()
 {
   // Default CTOR
 }
@@ -114,6 +115,9 @@ AliFMDRawReader::ReadAdcs(TClonesArray* array)
     AliError("No TClonesArray passed");
     return kFALSE;
   }
+  const UShort_t kUShortMax = (1 << 16) - 1;
+  fSeen.Reset(kUShortMax);
+  
   //  if (!fReader->ReadHeader()) {
   //    AliError("Couldn't read header");
   //    return kFALSE;
@@ -125,8 +129,8 @@ AliFMDRawReader::ReadAdcs(TClonesArray* array)
   input.SetShortDataHeader(!pars->HasCompleteHeader());
 
   UShort_t stripMin = 0;
-  UShort_t stripMax = 127;
-  UShort_t preSamp  = 14+5;
+  UShort_t stripMax = 0; // 127;
+  UShort_t preSamp  = 0; // 14+5;
   
   UInt_t ddl    = 0;
   UInt_t rate   = 0;
@@ -146,44 +150,56 @@ AliFMDRawReader::ReadAdcs(TClonesArray* array)
     }
 
     AliFMDDebug(5, ("Read channel 0x%x of size %d", hwaddr, last));
-    UShort_t det, sec, str;
-    Char_t   ring;
-    if (!pars->Hardware2Detector(ddl, hwaddr, det, ring, sec, str)) {
-      AliError(Form("Failed to get detector id from DDL %d "
-                   "and hardware address 0x%x", ddl, hwaddr));
-      continue;
-    }
+
     
-    rate     = pars->GetSampleRate(det, ring, sec, str);
-    stripMin = pars->GetMinStrip(det, ring, sec, str);
-    stripMax = pars->GetMaxStrip(det, ring, sec, str);
-    preSamp  = pars->GetPreSamples(det, ring, sec, str);
-    AliFMDDebug(5, ("DDL 0x%04x, address 0x%03x maps to FMD%d%c[%2d,%3d]", 
-                      ddl, hwaddr, det, ring, sec, str));
     
     // Loop over the `timebins', and make the digits
     for (size_t i = 0; i < last; i++) {
-      if (i < preSamp) continue;
-      Int_t    n      = array->GetEntriesFast();
-      Short_t  curStr = str + stripMin + (i-preSamp) / rate;
-      if ((curStr-str) > stripMax) {
-       // AliInfo(Form("timebin %4d -> (%3d+%3d+(%4d-%d)/%d) -> %d", 
-       //              i, str, stripMin, i, preSamp, rate, curStr));
-       // AliError(Form("Current strip is %3d (0x%04x,0x%03x,%4d) "
-       //               "but DB says max is %3d (rate=%d)", 
-       //               curStr, ddl, hwaddr, i, str+stripMax, rate));
-       // Garbage timebins - ignore 
+      // if (i < preSamp) continue;
+
+      UShort_t det, sec, samp;
+      Short_t str;
+      Char_t   ring;
+      if (!pars->Hardware2Detector(ddl, hwaddr, i, det, ring, sec, str, samp)) {
+       AliError(Form("Failed to get detector id from DDL %d, "
+                     "hardware address 0x%03x, timebin %d", ddl, hwaddr, i));
        continue;
       }
-      AliFMDDebug(5, ("making digit for FMD%d%c[%2d,%3d] from sample %4d", 
-                      det, ring, sec, curStr, i));
-      new ((*array)[n]) AliFMDDigit(det, ring, sec, curStr, data[i], 
-                                   (rate >= 2 ? data[i+1] : 0),
-                                   (rate >= 3 ? data[i+2] : 0),
-                                   (rate >= 4 ? data[i+3] : 0));
-      if (rate >= 2) i++;
-      if (rate >= 3) i++;
-      if (rate >= 4) i++;
+      AliFMDDebug(10, ("0x%04x/0x%03x/%04d maps to FMD%d%c[%2d,%3d]-%d", 
+                     ddl, hwaddr, i, det, ring, sec, str, samp));
+      if (str < 0) { 
+       AliFMDDebug(8, ("Got presamples at timebin %d", i));
+       continue;
+      }
+      
+      stripMin = pars->GetMinStrip(det, ring, sec, str);
+      stripMax = pars->GetMaxStrip(det, ring, sec, str);
+      preSamp  = pars->GetPreSamples(det, ring, sec, str);
+      rate     = pars->GetSampleRate(det, ring, sec, str);
+      Short_t lstrip = (i - preSamp) / rate + stripMin;
+      
+      AliFMDDebug(15, ("Checking if strip %d (%d) in range [%d,%d]", 
+                     lstrip, str, stripMin, stripMax));
+      if (lstrip < stripMin || lstrip > stripMax) {
+       AliFMDDebug(5, ("FMD%d%c[%02d,%03d]-%d out of range (%3d->%3d)", 
+                       det, ring, sec, samp, str, stripMin, stripMax));
+       continue;
+      }
+      
+      // Check the cache of indicies
+      Int_t idx = fSeen(det, ring, sec, str);
+      if (idx == kUShortMax) { 
+       // We haven't seen this strip yet. 
+       fSeen(det, ring, sec, str) = idx = array->GetEntriesFast();
+       AliFMDDebug(7,("making digit for FMD%d%c[%2d,%3d]-%d from timebin %4d", 
+                      det, ring, sec, str, samp, i));
+       new ((*array)[idx]) AliFMDDigit(det, ring, sec, str);
+      }
+      AliFMDDigit* digit = static_cast<AliFMDDigit*>(array->At(idx));
+      AliFMDDebug(10,
+                 ("Setting from FMD%d%c[%2d,%3d]-%d from timebin %4d = %4d", 
+                  det, ring, sec, str, samp, i, data[i]));
+      digit->SetCount(samp, data[i]);
     }
   }
   return kTRUE;
@@ -388,7 +404,8 @@ AliFMDRawReader::Exec(Option_t*)
 Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate, 
                                     AliFMDCalibStripRange* stripRange, 
                                     TArrayS &pulseSize, 
-                                    TArrayS &pulseLength) {
+                                    TArrayS &pulseLength) 
+{
 
   AliFMDDebug(0, ("Start of SOD/EOD"));
   
@@ -398,6 +415,7 @@ Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate,
   UInt_t strip_high[18];
   UInt_t pulse_size[18];
   UInt_t pulse_length[18];  
+  AliFMDParameters* param = AliFMDParameters::Instance();
   
   while(fReader->ReadNextData(fData)) {
     
@@ -410,249 +428,251 @@ Bool_t AliFMDRawReader::ReadSODevent(AliFMDCalibSampleRate* sampleRate,
     ULong_t nWords      = GetNwords();
     UInt_t  trailerLast = Get32bitWord(nWords);
     AliFMDDebug(20, (" # Bytes: %d, # Words: %d, Last word: 0x%08x", 
-                   fNbytes, nWords, trailerLast));
+                    fNbytes, nWords, trailerLast));
     if ((trailerLast & 0xFFFF0000) != 0xAAAA0000) {
       AliWarning(Form("Last word 0x%08x does not match RCU II trailer",
                      trailerLast));
       return kFALSE;
     }
     
-    ULong_t nTrailerWords = trailerLast & 0x7f; // 7 last bits is size of trailer
-    
+    // 7 last bits is size of trailer
+    ULong_t nTrailerWords = trailerLast & 0x7f; 
     ULong_t nPayloadWords = Get32bitWord(nWords - nTrailerWords+1);
+
     AliFMDDebug(20, (" # trailer words: %d, # payload words: %d", 
-                   nTrailerWords, nPayloadWords));
+                    nTrailerWords, nPayloadWords));
     
     
     for (ULong_t i = 1; i <= nPayloadWords ; i++) {
       UInt_t payloadWord = Get32bitWord(i);
       
-       // address is only 24 bit
-       UInt_t address       = (0xffffff & payloadWord);
-       UInt_t type          = ((address >> 21) & 0xf);
-       UInt_t error         = ((address >> 20) & 0x1);
-       UInt_t bcast         = ((address >> 18) & 0x1);
-       UInt_t bc_not_altro  = ((address >> 17) & 0x1);
-       UInt_t board         = ((address >> 12) & 0x1f);
-       UInt_t instruction   = 0;
-       UInt_t chip          = 0;
-       UInt_t channel       = 0;
-       if(bc_not_altro)
-         instruction        = address & 0xfff;
-       else {
-         chip               = ((address >> 9) & 0x7);
-         channel            = ((address >> 5) & 0x5);
-         instruction        = (address & 0x1f);
-       }
-       
-       Bool_t readDataWord = kFALSE;
-       switch(type) {
-       case 0x0: // Fec read
-         readDataWord = kTRUE;  
-       case 0x1: // Fec cmd
-       case 0x2: // Fec write
-         i++;  
-         break;
-       case 0x4: // Loop
-       case 0x5: // Wait
-         break;
-       case 0x6: // End sequence
-       case 0x7: // End Mem
-         i = nPayloadWords + 1;
-         break;
-       default:    
-         break;
-       }
-       
+      // address is only 24 bit
+      UInt_t address       = (0xffffff & payloadWord);
+      UInt_t type          = ((address >> 21) & 0xf);
+      UInt_t error         = ((address >> 20) & 0x1);
+      UInt_t bcast         = ((address >> 18) & 0x1);
+      UInt_t bc_not_altro  = ((address >> 17) & 0x1);
+      UInt_t board         = ((address >> 12) & 0x1f);
+      UInt_t instruction   = 0;
+      UInt_t chip          = 0;
+      UInt_t channel       = 0;
+      if(bc_not_altro)
+       instruction        = address & 0xfff;
+      else {
+       chip               = ((address >> 9) & 0x7);
+       channel            = ((address >> 5) & 0x5);
+       instruction        = (address & 0x1f);
+      }
        
+      Bool_t readDataWord = kFALSE;
+      switch(type) {
+      case 0x0: // Fec read
+       readDataWord = kTRUE;  
+      case 0x1: // Fec cmd
+      case 0x2: // Fec write
+       i++;  
+       break;
+      case 0x4: // Loop
+      case 0x5: // Wait
+       break;
+      case 0x6: // End sequence
+      case 0x7: // End Mem
+       i = nPayloadWords + 1;
+       break;
+      default:    
+       break;
+      }
        
-       if(!readDataWord)  //Don't read unless we have a FEC_RD
-         continue;
-       UInt_t dataWord      = Get32bitWord(i);
-       UInt_t data          = (0xFFFFF & dataWord) ;
-       //UInt_t data          = (0xFFFF & dataWord) ;
+      //Don't read unless we have a FEC_RD
+      if(!readDataWord)  continue;
+
+      UInt_t dataWord      = Get32bitWord(i);
+      UInt_t data          = (0xFFFFF & dataWord) ;
+      //UInt_t data          = (0xFFFF & dataWord) ;
        
-       if(error) {
-         AliWarning(Form("error bit detected at Word 0x%06x; "
-                         "error % d, type %d, bc_not_altro %d, "
-                         "bcast %d, board 0x%02x, chip 0x%x, "
-                         "channel 0x%02x, instruction 0x%03x",
-                         address, error, type, bc_not_altro, 
-                         bcast,board,chip,channel,instruction));
-         
-         
-         //process error
-         continue;
-       }
+      if(error) {
+       AliWarning(Form("error bit detected at Word 0x%06x; "
+                       "error % d, type %d, bc_not_altro %d, "
+                       "bcast %d, board 0x%02x, chip 0x%x, "
+                       "channel 0x%02x, instruction 0x%03x",
+                       address, error, type, bc_not_altro, 
+                       bcast,board,chip,channel,instruction));
+       //process error
+       continue;
+      }
        
        
-       switch(instruction) {
+      switch(instruction) {
          
-       case 0x01: break;  // First ADC T           
-       case 0x02: break; // I  3.3 V              
-       case 0x03: break; // I  2.5 V altro digital
-       case 0x04: break; // I  2.5 V altro analog 
-       case 0x05: break; // I  2.5 V VA           
-       case 0x06: break; // First ADC T           
-       case 0x07: break; // I  3.3 V              
-       case 0x08: break; // I  2.5 V altro digital
-       case 0x09: break; // I  2.5 V altro analog 
-       case 0x0A: break; // I  2.5 V VA           
-       case 0x2D: break; // Second ADC T           
-       case 0x2E: break; // I  1.5 V VA            
-       case 0x2F: break; // I -2.0 V               
-       case 0x30: break; // I -2.0 V VA            
-       case 0x31: break; //    2.5 V Digital driver
-       case 0x32: break; // Second ADC T           
-       case 0x33: break; // I  1.5 V VA            
-       case 0x34: break; // I -2.0 V               
-       case 0x35: break; // I -2.0 V VA            
-       case 0x36: break; //    2.5 V Digital driver
-       case 0x37: break; // Third ADC T             
-       case 0x38: break; // Temperature sens. 1     
-       case 0x39: break; // Temperature sens. 2     
-       case 0x3A: break; // U  2.5 altro digital (m)
-       case 0x3B: break; // U  2.5 altro analog (m) 
-       case 0x3C: break; // Third ADC T             
-       case 0x3D: break; // Temperature sens. 1     
-       case 0x3E: break; // Temperature sens. 2     
-       case 0x3F: break; // U  2.5 altro digital (m)
-       case 0x40: break; // U  2.5 altro analog (m) 
-       case 0x41: break; // Forth ADC T  
-       case 0x42: break; // U  2.5 VA (m)
-       case 0x43: break; // U  1.5 VA (m)
-       case 0x44: break; // U -2.0 VA (m)
-       case 0x45: break; // U -2.0 (m)   
-       case 0x46: break; // Forth ADC T  
-       case 0x47: break; // U  2.5 VA (m)
-       case 0x48: break; // U  1.5 VA (m)
-       case 0x49: break; // U -2.0 VA (m)
-       case 0x4A: break;  // U -2.0 (m)   
-         // Counters 
-       case 0x0B: break; // L1 trigger CouNTer
-       case 0x0C: break; // L2 trigger CouNTer
-       case 0x0D: break; // Sampling CLK CouNTer
-       case 0x0E: break; // DSTB CouNTer
-         // Test mode 
-       case 0x0F: break; // Test mode word
-       case 0x10: break; // Undersampling ratio.
-         // Configuration and status 
-       case 0x11: break; // Config/Status Register 0
-       case 0x12: break; // Config/Status Register 1
-       case 0x13: break; // Config/Status Register 2
-       case 0x14: break; // Config/Status Register 3
-       case 0x15: break; // Free
-     // Comands:
-       case 0x16: break; // Latch L1, L2, SCLK Counters
-       case 0x17: break; // Clear counters
-       case 0x18: break; // Clear CSR1
-       case 0x19: break; // rstb ALTROs
-       case 0x1A: break; // rstb BC
-       case 0x1B: break; // Start conversion
-       case 0x1C: break; // Scan event length
-       case 0x1D: break; // Read event length
-       case 0x1E: break; // Start test mode
-       case 0x1F: break; // Read acquisition memory
-         // FMD
-       case 0x20: break; // FMDD status
-       case 0x21: break; // L0 counters
-       case 0x22: break; // FMD: Wait to hold
-       case 0x23: break; // FMD: L1 timeout
-       case 0x24: break; // FMD: L2 timeout
-       case 0x25: // FMD: Shift clk 
-         shift_clk[board] = ((data >> 8 ) & 0xFF); 
-         break; 
-       case 0x26: // FMD: Strips 
-         strip_low[board]  = ((data >> 0 ) & 0xFF); 
-         strip_high[board] = ((data >> 8 ) & 0xFF);  
-         break; 
-       case 0x27: // FMD: Cal pulse 
-         pulse_size[board]  =  ((data >> 8 ) & 0xFF);
-         break; 
-       case 0x28: break; // FMD: Shape bias
-       case 0x29: break; // FMD: Shape ref
-       case 0x2A: break; // FMD: Preamp ref
-       case 0x2B: // FMD: Sample clk 
-         sample_clk[board] = ((data >> 8 ) & 0xFF); 
-         break; 
-       case 0x2C: break; // FMD: Commands
-       case 0x4B: // FMD: Cal events 
-         pulse_length[board] = ((data >> 0 ) & 0xFF);
-         break; 
-       default: break;
+      case 0x01: break;  // First ADC T           
+      case 0x02: break; // I  3.3 V              
+      case 0x03: break; // I  2.5 V altro digital
+      case 0x04: break; // I  2.5 V altro analog 
+      case 0x05: break; // I  2.5 V VA           
+      case 0x06: break; // First ADC T           
+      case 0x07: break; // I  3.3 V              
+      case 0x08: break; // I  2.5 V altro digital
+      case 0x09: break; // I  2.5 V altro analog 
+      case 0x0A: break; // I  2.5 V VA           
+      case 0x2D: break; // Second ADC T           
+      case 0x2E: break; // I  1.5 V VA            
+      case 0x2F: break; // I -2.0 V               
+      case 0x30: break; // I -2.0 V VA            
+      case 0x31: break; //    2.5 V Digital driver
+      case 0x32: break; // Second ADC T           
+      case 0x33: break; // I  1.5 V VA            
+      case 0x34: break; // I -2.0 V               
+      case 0x35: break; // I -2.0 V VA            
+      case 0x36: break; //    2.5 V Digital driver
+      case 0x37: break; // Third ADC T             
+      case 0x38: break; // Temperature sens. 1     
+      case 0x39: break; // Temperature sens. 2     
+      case 0x3A: break; // U  2.5 altro digital (m)
+      case 0x3B: break; // U  2.5 altro analog (m) 
+      case 0x3C: break; // Third ADC T             
+      case 0x3D: break; // Temperature sens. 1     
+      case 0x3E: break; // Temperature sens. 2     
+      case 0x3F: break; // U  2.5 altro digital (m)
+      case 0x40: break; // U  2.5 altro analog (m) 
+      case 0x41: break; // Forth ADC T  
+      case 0x42: break; // U  2.5 VA (m)
+      case 0x43: break; // U  1.5 VA (m)
+      case 0x44: break; // U -2.0 VA (m)
+      case 0x45: break; // U -2.0 (m)   
+      case 0x46: break; // Forth ADC T  
+      case 0x47: break; // U  2.5 VA (m)
+      case 0x48: break; // U  1.5 VA (m)
+      case 0x49: break; // U -2.0 VA (m)
+      case 0x4A: break;  // U -2.0 (m)   
+       // Counters 
+      case 0x0B: break; // L1 trigger CouNTer
+      case 0x0C: break; // L2 trigger CouNTer
+      case 0x0D: break; // Sampling CLK CouNTer
+      case 0x0E: break; // DSTB CouNTer
+       // Test mode 
+      case 0x0F: break; // Test mode word
+      case 0x10: break; // Undersampling ratio.
+       // Configuration and status 
+      case 0x11: break; // Config/Status Register 0
+      case 0x12: break; // Config/Status Register 1
+      case 0x13: break; // Config/Status Register 2
+      case 0x14: break; // Config/Status Register 3
+      case 0x15: break; // Free
+       // Comands:
+      case 0x16: break; // Latch L1, L2, SCLK Counters
+      case 0x17: break; // Clear counters
+      case 0x18: break; // Clear CSR1
+      case 0x19: break; // rstb ALTROs
+      case 0x1A: break; // rstb BC
+      case 0x1B: break; // Start conversion
+      case 0x1C: break; // Scan event length
+      case 0x1D: break; // Read event length
+      case 0x1E: break; // Start test mode
+      case 0x1F: break; // Read acquisition memory
+       // FMD
+      case 0x20: break; // FMDD status
+      case 0x21: break; // L0 counters
+      case 0x22: break; // FMD: Wait to hold
+      case 0x23: break; // FMD: L1 timeout
+      case 0x24: break; // FMD: L2 timeout
+      case 0x25: // FMD: Shift clk 
+       shift_clk[board] = ((data >> 8 ) & 0xFF); 
+       break; 
+      case 0x26: // FMD: Strips 
+       strip_low[board]  = ((data >> 0 ) & 0xFF); 
+       strip_high[board] = ((data >> 8 ) & 0xFF);  
+       break; 
+      case 0x27: // FMD: Cal pulse 
+       pulse_size[board]  =  ((data >> 8 ) & 0xFF);
+       break; 
+      case 0x28: break; // FMD: Shape bias
+      case 0x29: break; // FMD: Shape ref
+      case 0x2A: break; // FMD: Preamp ref
+      case 0x2B: // FMD: Sample clk 
+       sample_clk[board] = ((data >> 8 ) & 0xFF); 
+       break; 
+      case 0x2C: break; // FMD: Commands
+      case 0x4B: // FMD: Cal events 
+       pulse_length[board] = ((data >> 0 ) & 0xFF);
+       break; 
+      default: break;
        
-       }
-       AliFMDDebug(50, ("instruction 0x%x, dataword 0x%x",instruction,dataWord));
+      }
+      AliFMDDebug(50, ("instruction 0x%x, dataword 0x%x",
+                      instruction,dataWord));
     }
     
-    UShort_t det,sector,strip;
+    UShort_t det,sector;
+    Short_t strip;
     Char_t ring;
    
     const UInt_t boards[4] = {0,1,16,17};
     for(Int_t i=0;i<4;i++) {
-      if(ddl==0 && (i==1 || i==3))
-       continue;
+      if(ddl==0 && (i==1 || i==3)) continue;
+
       UInt_t chip =0, channel=0;
-      
-      AliFMDParameters::Instance()->Hardware2Detector(ddl,boards[i],chip,channel,det,ring,sector,strip);
+      param->Hardware2Detector(ddl,boards[i],chip,channel,
+                              det,ring,sector,strip);
      
       UInt_t samplerate = 1;
-      
       if(sample_clk[boards[i]] == 0) {
-         if(ddl == 0) {
-           Int_t sample1 = sample_clk[boards[0]];
-           Int_t sample2 = sample_clk[boards[2]];          
-           if(sample1) sample_clk[boards[i]] = sample1;
-           else sample_clk[boards[i]] = sample2;
-         }
-         
-         if(ddl!=0) {
-           Int_t sample1 = sample_clk[boards[0]];
-           Int_t sample2 = sample_clk[boards[1]];
-           Int_t sample3 = sample_clk[boards[2]];
-           Int_t sample4 = sample_clk[boards[3]];
-           Int_t agreement = 0;
-           if(sample1 == sample2) agreement++;
-           if(sample1 == sample3) agreement++;
-           if(sample1 == sample4) agreement++;
-           if(sample2 == sample3) agreement++;
-           if(sample2 == sample4) agreement++;
-           if(sample3 == sample4) agreement++;
+       if(ddl == 0) {
+         Int_t sample1 = sample_clk[boards[0]];
+         Int_t sample2 = sample_clk[boards[2]];            
+         if(sample1) sample_clk[boards[i]] = sample1;
+         else sample_clk[boards[i]] = sample2;
+       }
+       else {
+         Int_t sample1 = sample_clk[boards[0]];
+         Int_t sample2 = sample_clk[boards[1]];
+         Int_t sample3 = sample_clk[boards[2]];
+         Int_t sample4 = sample_clk[boards[3]];
+         Int_t agreement = 0;
+         if(sample1 == sample2) agreement++;
+         if(sample1 == sample3) agreement++;
+         if(sample1 == sample4) agreement++;
+         if(sample2 == sample3) agreement++;
+         if(sample2 == sample4) agreement++;
+         if(sample3 == sample4) agreement++;
            
-           Int_t idx = 0;
-           if(i<3) idx = i+1;
-           else  idx = i-1;
-           if(agreement == 3) {
-             sample_clk[boards[i]] = sample_clk[boards[idx]];
-             shift_clk[boards[i]] = shift_clk[boards[idx]];
-             strip_low[boards[i]] = strip_low[boards[idx]];
-             strip_high[boards[i]] = strip_high[boards[idx]];
-             pulse_length[boards[i]] = pulse_length[boards[idx]];
-             pulse_size[boards[i]] = pulse_size[boards[idx]];
-             AliFMDDebug(0, ("Vote taken for ddl %d, board 0x%x",ddl,boards[i]));
-           }
+         Int_t idx = 0;
+         if(i<3) idx = i+1;
+         else  idx = i-1;
+         if(agreement == 3) {
+           sample_clk[boards[i]] = sample_clk[boards[idx]];
+           shift_clk[boards[i]] = shift_clk[boards[idx]];
+           strip_low[boards[i]] = strip_low[boards[idx]];
+           strip_high[boards[i]] = strip_high[boards[idx]];
+           pulse_length[boards[i]] = pulse_length[boards[idx]];
+           pulse_size[boards[i]] = pulse_size[boards[idx]];
+           AliFMDDebug(0, ("Vote taken for ddl %d, board 0x%x",
+                           ddl,boards[i]));
          }
+       }
       } 
       
       if(sample_clk[boards[i]])
-      samplerate = shift_clk[boards[i]]/sample_clk[boards[i]];
+       samplerate = shift_clk[boards[i]]/sample_clk[boards[i]];
       sampleRate->Set(det,ring,sector,0,samplerate);
-      stripRange->Set(det,ring,sector,0,strip_low[boards[i]],strip_high[boards[i]]);
+      stripRange->Set(det,ring,sector,0,
+                     strip_low[boards[i]],strip_high[boards[i]]);
       
       AliFMDDebug(20, ("det %d, ring %c, ",det,ring));
-      pulseLength.AddAt(pulse_length[boards[i]],GetHalfringIndex(det,ring,boards[i]/16));
-      pulseSize.AddAt(pulse_size[boards[i]],GetHalfringIndex(det,ring,boards[i]/16));
+      pulseLength.AddAt(pulse_length[boards[i]],
+                       GetHalfringIndex(det,ring,boards[i]/16));
+      pulseSize.AddAt(pulse_size[boards[i]],
+                     GetHalfringIndex(det,ring,boards[i]/16));
       
       
       
       AliFMDDebug(20, (": Board: 0x%02x\n"
-                     "\tstrip_low  %3d, strip_high   %3d\n"
-                     "\tshift_clk  %3d, sample_clk   %3d\n"
-                     "\tpulse_size %3d, pulse_length %3d",
-                     boards[i], 
-                     strip_low[boards[i]], strip_high[boards[i]],
-                     shift_clk[boards[i]], sample_clk[boards[i]],
-                     pulse_size[boards[i]],pulse_length[boards[i]]));
+                      "\tstrip_low  %3d, strip_high   %3d\n"
+                      "\tshift_clk  %3d, sample_clk   %3d\n"
+                      "\tpulse_size %3d, pulse_length %3d",
+                      boards[i], 
+                      strip_low[boards[i]], strip_high[boards[i]],
+                      shift_clk[boards[i]], sample_clk[boards[i]],
+                      pulse_size[boards[i]],pulse_length[boards[i]]));
     }
     
   }
@@ -695,7 +715,8 @@ UInt_t AliFMDRawReader::Get32bitWord(Int_t idx)
   return word;
 }
 //_____________________________________________________________________ 
-Int_t AliFMDRawReader::GetHalfringIndex(UShort_t det, Char_t ring, UShort_t board) {
+Int_t AliFMDRawReader::GetHalfringIndex(UShort_t det, Char_t ring, 
+                                       UShort_t board) {
 
   UShort_t iring  =  (ring == 'I' ? 1 : 0);