HarpoDccDecode.cxx

Go to the documentation of this file.

//
// File HarpoDccDecode.cxx 
//
#include "HarpoDccDecode.h"
#include "HarpoDebug.h"
#include <cstdlib>
#include <cstring>
#include <iostream>
using namespace std; 

ClassImp(HarpoDccDecode)
// code from subroutine decodag
HarpoDccDecode::HarpoDccDecode()
{
  Int_t raw, plane,iafter, strip, io, iv;
  Int_t isignal = 0;
  Int_t ibruit = 0;
  for (Int_t i=0; i<NSTRIPS; i++) 
    {
      iafter = i/76+1;
      strip = i + 1 - (iafter-1)*76 ;
      iafter = 5 - iafter;
      
      if(strip == 1)
        { 
          ibruit = 0;
          isignal = 0;
        }
      
      if(strip==13||strip==26||strip==51||strip==64) 
        {
          ibruit = ibruit +1;
          raw = ibruit + 288 + (iafter-1)*4;
        }
      else 
        {
          isignal = isignal + 1;
          raw = isignal + (iafter-1)*72;
        }

      for (plane= 0; plane < 2; plane ++)
        {
          if (plane == 0) 
            {
              io = raw%2;
              iv = -1 + 2*io;
              this->Row2Strip[plane][i] = raw + iv;
            }
          else 
            {
              this->Row2Strip[plane][i]  = raw; 
            }
        }
      // this->Strip2Raw[plane-1, raw] = i;
    }
}

// May be made inline
Int_t HarpoDccDecode::decode(Int_t raw, Int_t plane)
{ 
  if(gHarpoDebug>2) 
    Info("decode","row = %d, plane = %d",raw,plane);
  Int_t tmp = Row2Strip[plane][raw]-1;
  return tmp;
}

// May be made inline
Int_t HarpoDccDecode::invdecode(Int_t strip, Int_t /* plane */ )
{
  return strip;
}




void HarpoDccDecode::print()
{
  for(Int_t i=0; i<NSTRIPS; i++) {
    for(Int_t j=0; j<NPLANES; j++) {
      cout << this->Row2Strip[j][i]  << endl;
    }
  }
}