FSUDAQ_Qt6/Aux/AggSeparator.h

#ifndef AGGSEPARATOR_H
#define AGGSEPARATOR_H

#include <stdio.h>
#include <string>
#include <sstream>
#include <cmath>
#include <cstring>  ///memset
#include <iostream> ///cout
#include <sstream>
#include <iomanip> // for setw
#include <algorithm>
#include <bitset>
#include <vector>
#include <sys/stat.h>

#define NumCoupledChannel 8 // the max numnber of Coupled/RegChannel is 8 for PHA, PSD, QDC


std::vector<std::string> AggSeperator(std::string inFileName, std::string saveFolder = "./", short verbose = false){

  printf("================ AggSeperator \n");

  std::vector<std::string> outputFileList;
  outputFileList.clear();

  FILE * file = fopen(inFileName.c_str(), "r");

  if( file == NULL ) {
    printf("file : %s cannot be open. exit program.\n", inFileName.c_str());
    return outputFileList;
  }

  std::string folder = "";
  size_t found = inFileName.find_last_of('/');
  std::string fileName = inFileName;
  if( found != std::string::npos ){
    folder = inFileName.substr(0, found + 1);
    fileName = inFileName.substr(found +1 );
  }

  if( saveFolder.empty() ) saveFolder = "./";
  if( saveFolder.back() != '/') saveFolder += '/';

  printf("    fileName : %s \n", fileName.c_str());
  printf("      folder : %s \n", folder.c_str());
  printf(" save folder : %s\n", saveFolder.c_str());

  char * buffer = nullptr;
  unsigned int word; // 4 bytes = 32 bits

  bool newFileFlag[NumCoupledChannel];

  for( int i = 0; i < NumCoupledChannel; i++){
    newFileFlag[i] = true;
    outputFileList.push_back( saveFolder + fileName + "." + std::to_string(i));
  }

  do{

    size_t dummy = fread(&word, 4, 1, file);
    if( dummy != 1 ){
      printf("=====> End of File.\n");
      break;
    }

    short header = ((word >> 28 ) & 0xF);
    if( header != 0xA ) {
      printf("header error. abort.\n");
      break;
    }
    //unsigned int aggSize = (word & 0x0FFFFFFF) * 4; ///byte

    dummy = fread(&word, 4, 1, file);
    //unsigned int BoardID = ((word >> 27) & 0x1F);
    //unsigned short pattern = ((word >> 8 ) & 0x7FFF );
    //bool BoardFailFlag =  ((word >> 26) & 0x1 );
    unsigned int ChannelMask = ( word & 0xFF ) ;

    dummy = fread(&word, 4, 1, file);
    unsigned int bdAggCounter = word;

    if( verbose ) printf("Agg counter : %u\n", bdAggCounter);

    dummy = fread(&word, 4, 1, file);
    //unsigned int aggTimeTag = word;

    for( int chMask = 0; chMask < NumCoupledChannel ; chMask ++ ){ 
      if( ((ChannelMask >> chMask) & 0x1 ) == 0 ) continue;
      if( verbose ) printf("==================== Dual/Group Channel Block, ch Mask : 0x%X (%d)\n", chMask *2, chMask );

      dummy = fread(&word, 4, 1, file);
      unsigned int dualChannelBlockSize = ( word & 0x7FFFFFFF ) * 4 ;

      if( verbose ) printf("dual channel size : %d words\n", dualChannelBlockSize / 4);

      buffer = new char[dualChannelBlockSize];
      fseek(file, -4, SEEK_CUR);
      dummy = fread(buffer, dualChannelBlockSize, 1, file);

      FILE * haha = nullptr;
      
      if( newFileFlag[chMask] ) {
        haha = fopen( outputFileList[chMask].c_str(), "wb");
        newFileFlag[chMask] = false;
      }else{
        haha = fopen( outputFileList[chMask].c_str(), "a+");
      }

      fwrite(buffer, dualChannelBlockSize, 1, haha);

      fclose(haha); 
      
    }

  }while( !feof(file));

  fclose(file);

  printf("======================= Duel channels seperated \n");

  for( int i = NumCoupledChannel -1 ; i >= 0 ; i--){
    if( newFileFlag[i] == true ) outputFileList.erase(outputFileList.begin() + i );
  }

  return outputFileList;

}

#endif