130 lines
3.5 KiB
C++
130 lines
3.5 KiB
C++
#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 |