/************************************** * * root script for using the BinReader.h * * ***************************************/ #include "BinReader.h" #include "TH1.h" #include "TStyle.h" #include "TCanvas.h" #define NChannel 16 TH1I * hBoard; TH1I * hChannel; TH1I * he[NChannel]; // BinReader * reader = new BinReader("/home/ryan/ExpData/fsu_testing/testing_014_single.bin"); BinReader * reader = new BinReader("data/run_123/UNFILTERED/Data_CH9@V1725S_19555_run_123.BIN"); void test(){ // reader->ReadBlock(); // reader->data.Print(); reader->ScanNumHit(); printf("################## Number of data Block : %llu \n", reader->GetNumHit()); unsigned int count = 0; do{ reader->ReadBlock(); reader->data.Print(); count ++; }while( !reader->IsEndOfFile() ); printf(">>>>>>>>>>>>>>. count : %ld\n", count); reader->PrintStatus(1); printf("\n\n\n"); // hBoard = new TH1I ("hBoard", "BoardID", 100, 0, 10); // hChannel = new TH1I ("hChannel", "Channel", NChannel, 0, NChannel); // for( int i = 0 ; i < NChannel ; i++){ // he[i] = new TH1I(Form("he%0d", i), Form("he%0d", i), 500, 0, 10000); // } // for( int i = 0; i < reader->GetNumberOfBlock(); i++){ // reader->ReadBlock(); // hBoard->Fill( reader->data.BoardID ); // unsigned short ch = reader->data.Channel; // hChannel->Fill( ch); // he[ch]->Fill( reader->data.Energy); // if( i < 10 ) reader->data.Print(); // if( i % 5000000 == 0 ) { // printf("------- %d \n", i); // reader->data.Print(); // } // } // gStyle->SetOptStat("neiou"); // TCanvas * haha = new TCanvas("haha", "haha", 4*400, 5*400); // haha->Divide(4,5); // haha->cd(1); hBoard->Draw(); // haha->cd(2); hChannel->Draw(); // for( int i = 0; i < 16; i++){ // haha->cd(i+5); // he[i]->Draw(); // } }