TOctalFADCBankReader.cpp

Go to the documentation of this file.

#include "TOctalFADCBankReader.h"
#include "midas.h"
#include <stdio.h>
#include <algorithm>


using std::string;
using std::vector;

TOctalFADCBankReader::TOctalFADCBankReader(string bankname)
  : fBankName(bankname)
{
}

void TOctalFADCBankReader::ProcessEvent(EVENT_HEADER* pheader, void *pevent)
{
  ClearEvent();

  unsigned char* raw; // Points at the raw data, one byte at a time.
  int bankSize = bk_locate(pevent,fBankName.c_str(),&raw);
  int nBlocks = bankSize / kDataBlockSize;

  vector<int> islandSamples;
  islandSamples.reserve(kSamplesArrayInitialSize);

  int islandTimestamp = 0; // time of first 4-sample block in stitched island
  int lastTimestamp = 0; // To see the end of the last island.
  bool firstIsland = true; // Beginning of first island is not end of an island
  for(int i=0; i < nBlocks; i++){
    // data format:
    //
    //  bits 
    //  79          FADC 0 / 1
    //  78-52       timestamp
    //  51-48       overflow
    //  47-36       sampleB0
    //  35-24       sampleA0
    //  23-12       sampleB1
    //  11-0        sampleA1
    //  
    //  Island samples in the order A1 B1 A0 B0

    int timestamp = ((raw[i*10+0] & 0x7f) << 20) |
                     (raw[i*10+1] << 12) |
                     (raw[i*10+2] << 4) |
                     (raw[i*10+3] >> 4);
    bool overflowB0 = ((raw[i*10+3] & 0x08) != 0);
    bool overflowA0 = ((raw[i*10+3] & 0x04) != 0);
    bool overflowB1 = ((raw[i*10+3] & 0x02) != 0);
    bool overflowA1 = ((raw[i*10+3] & 0x01) != 0);
    int sampleB0 = (overflowB0 << 12) |
                   (raw[i*10+4] << 4) |
                   (raw[i*10+5] >> 4);
    int sampleA0 = (overflowA0 << 12) |
                   ((raw[i*10+5] & 0xf) << 8) |
                   (raw[i*10+6]);
    int sampleB1 = (overflowB1 << 12) |
                   (raw[i*10+7] << 4) |
                   (raw[i*10+8] >> 4);
    int sampleA1 = (overflowA1 << 12) |
                   ((raw[i*10+8] & 0xf) << 8) |
                   (raw[i*10+9]);

    if(timestamp != lastTimestamp + 1) {
      if(!firstIsland) {
        if(islandSamples.size() > 4) {
          // This is a new island, so put the old 
          // island on fData and start again.
          // Remember, time stamps are actually 4 samples long
          int channel = fBankName.at(1) - 'a';
          fData.push_back(TOctalFADCIsland(islandTimestamp*4,islandSamples,channel));


#if 0
          printf("Time = %d, nsamples = %d, char = %c , channel = %d\n",
            islandTimestamp, islandSamples.size(),fBankName.at(1),channel);
          for(int k=0; k<islandSamples.size(); k++){
            printf("%d ",islandSamples[k]);
          }
          printf("\n");
#endif
          islandSamples.clear();
        }
        // Ignore islands with 4 or less samples.
      }
      firstIsland = false;
      islandTimestamp = timestamp;
    }
    lastTimestamp = timestamp;

#if 0
    for(int k=0; k<10; k++){
      printf("%x ",raw[i*10 + k]);
    }
    printf("\n");
    printf("fadc=%x t=%2x %2x %2x %2x %2x\n", i, timestamp, 
           sampleA1, sampleB1, sampleA0,sampleB0);
#endif

    islandSamples.push_back(sampleA1);
    islandSamples.push_back(sampleB1);
    islandSamples.push_back(sampleA0);
    islandSamples.push_back(sampleB0);
  }

  if(fData.size()>1)
  {
    //Time sort islands
    std::sort(fData.begin(), fData.end(),TOctalFADCIsland::TimeSortFADCIslands());
    StitchIslands();
  } // Some Pulses are split over two separate islands. Put them together again
 // printf(".............. FADC islands sorted .............\n");
}

void TOctalFADCBankReader::StitchIslands()
{
  
  vector<TOctalFADCIsland> fData_temp;
  //fData_temp.clear();
  
  TOctalFADCIsland last = fData.at(0);
  //if(fData.size()>1){printf("Size first Island: %d\n",fData.at(1).GetNSamples());}
  fData_temp.push_back(last);
  int counter_temp = 0; int counter = 0;
  
  for(int i = 1; i<fData.size(); i++)
  {
    //if(fData.at(i).GetFADCChannel()==2){printf("Time new = %d, End time last = %d, Difference = %d, Size new = %d, Size old =%d\n",fData.at(i).GetTime(),last.GetTime() + last.GetNSamples(),fData.at(i).GetTime() - (last.GetTime() + last.GetNSamples()),fData.at(i).GetNSamples(),last.GetNSamples());}
    if(fData.at(i).GetTime() - (last.GetTime() + last.GetNSamples()) >2 )
    {
      last = fData.at(i);
      fData_temp.push_back(last);
      counter_temp++;
      //printf("OK\n");
    }
    else
    {
      vector<int>& samples = fData_temp.at(counter_temp).GetSampleVector();
      //int timestamp = fData_temp.at(counter_temp).GetTime();
      //int channel = fData_temp.at(counter_temp).GetFADCChannel();
      //if(fData.at(i).GetFADCChannel()==2){printf("Stitch!!! : Size samples before = %d",samples.size());}

      for(int j=0;j<fData.at(i).GetNSamples();j++)
      {
        samples.push_back(fData.at(i).GetSample(j));
      }
      //if(fData.at(i).GetFADCChannel()==2){printf("   Size samples after = %d \n",samples.size());}  
      //fData_temp.erase(fData_temp.end());
      //fData_temp.push_back(TOctalFADCIsland(timestamp,samples,channel));
      counter ++;
      last = fData_temp.at(counter_temp);
    }
  }

//printf("Size original: %d, size stitched: %d, counter %d\n",fData.size(),fData_temp.size(),counter);

fData.clear();
fData = fData_temp;
fData_temp.clear();

//printf("Size new: %d, size temp: %d, counter %d , counter_temp %d\n",fData.size(),fData_temp.size(),counter,counter_temp);

}