--- /dev/null
+CC=g++
+CFLAGS=-c -g -Wall -std=c++11 `root-config --cflags`
+LDFLAGS=`root-config --glibs`
+CLDFLAGS = -std=c++11 `root-config --cflags --glibs` -lboost_system -lboost_filesystem
+SOURCES= main.C unpack.C unpack.h help.h header.h
+OBJECTS=$(subst .cc,.o,$(SOURCES))
+EXECUTABLE=Mi26NM
+
+all: $(SOURCES) $(EXECUTABLE)
+
+$(EXECUTABLE): $(OBJECTS)
+ $(CC) $(CLDFLAGS) $(OBJECTS) -o $@
+
+
+clean:
+ rm -f $(OBJS)
\ No newline at end of file
--- /dev/null
+#include <stdlib.h>
+#include <TClass.h>
+#include <fstream>
+#include <iostream>
+#include "Riostream.h"
+#include <string>
+#include "TString.h"
+#include "TMath.h"
+#include <TFile.h>
+#include <TNtuple.h>
+#include <TH1F.h>
+#include <TH2F.h>
+#include "TH1.h"
+#include "TF1.h"
+#include "TROOT.h"
+#include <TApplication.h>
+#include "TCanvas.h"
+#include "TObject.h"
+#include "TStyle.h"
+#include "TEnv.h"
+#include "TGraph.h"
+#include "TGaxis.h"
+#include "TStyle.h"
+#include "TFrame.h"
+#include "TVector.h"
+#include "TVectorD.h"
+#include "Foption.h"
+#include "TRandom.h"
+#include "TRandom3.h"
+#include "TSpline.h"
+#include "TPaveStats.h"
+#include "TVirtualFitter.h"
+#include "TVirtualPad.h"
+//#include "TVirtualGraphPainter.h"
+#include "TBrowser.h"
+#include "TSystem.h"
+#include "TPluginManager.h"
+#include "TPad.h"
+#include <TQObject.h>
+#include <TMultiGraph.h>
+#include <vector>
+#include "TControlBar.h"
+#include <string>
+#include <stdio.h>
+#include <stdlib.h>
+#include "TGraphErrors.h"
+#include "TMinuit.h"
+#include <iomanip>
+#include <sys/time.h>
+#include <TLegend.h>
+using namespace std;
\ No newline at end of file
--- /dev/null
+#ifndef HELP_H
+#define HELP_H
+
+#ifndef MAXSTATES
+#define MAXSTATES 1140
+#endif
+
+#include <iostream>
+#include <cstdlib>
+#include <unistd.h>
+#include <stdlib.h>
+#include <stddef.h>
+#include <stdio.h>
+#include <unistd.h>
+#include "TStopwatch.h"
+
+//####################################################################
+
+// Int_t print_progress(Float_t ProgressInPercent) {
+//
+// Int_t nCharacter = printf("\r [");
+// Int_t iBlock;
+//
+// for( iBlock=0; iBlock<ProgressInPercent/2; ++iBlock )
+// {
+// putchar( '=' );
+// nCharacter++;
+// }
+//
+// for( ; iBlock<100/2; ++iBlock )
+// {
+// putchar( ' ' );
+// nCharacter++;
+// }
+//
+// // Int_t sec = TimeInSec%60;
+// // Int_t min = TimeInSec/60;
+// // Int_t hours = TimeInSec/3600;
+//
+// nCharacter += printf( "] %4.0f%%", ProgressInPercent);
+// fflush( stdout );
+//
+// if((Int_t)ProgressInPercent==100)
+// {
+// printf("\r");
+// for(Int_t i=0; i<nCharacter; i++)
+// {
+// putchar( ' ' );
+// }
+// printf("\r ... done! \n");
+// }
+//
+// return nCharacter;
+// }
+
+//####################################################################
+
+// Int_t print_progress(Float_t ProgressInPercent, Int_t TimeInSec) {
+//
+// Int_t nCharacter = printf("\r [");
+// Int_t iBlock;
+//
+// for( iBlock=0; iBlock<ProgressInPercent/2; ++iBlock )
+// {
+// putchar( '=' );
+// nCharacter++;
+// }
+//
+// for( ; iBlock<100/2; ++iBlock )
+// {
+// putchar( ' ' );
+// nCharacter++;
+// }
+//
+// Int_t sec = TimeInSec%60;
+// Int_t min = TimeInSec/60;
+// Int_t hours = TimeInSec/3600;
+//
+// nCharacter += printf( "] %6.2f%% \t %3i:%02i:%02i ", ProgressInPercent, hours, min, sec );
+// fflush( stdout );
+//
+// if(ProgressInPercent==100)
+// {
+// printf("\r");
+// for(Int_t i=0; i<nCharacter; i++)
+// {
+// putchar( ' ' );
+// }
+// printf("\r ...done! \t ( %i:%02i:%02i ) \n", hours, min, sec );
+// }
+//
+// return nCharacter;
+// }
+
+//####################################################################
+
+struct frameInfo{
+// to be processed
+ UInt_t posInEvent;
+
+// to be stored
+ UInt_t event;
+ UInt_t framenr;
+ UInt_t sensorid;
+ UInt_t status;
+ UInt_t datalengthcount;
+ UInt_t datalength;
+ ULong_t timestamp;
+ ULong_t debug;
+
+ UInt_t states;
+ UInt_t row [MAXSTATES];
+ UInt_t column [MAXSTATES];
+ UInt_t nextn [MAXSTATES];
+ UInt_t overflow [MAXSTATES];
+
+ UInt_t clusters;
+ UInt_t pixelsPerCluster [MAXSTATES];
+ UInt_t statesPerCluster [MAXSTATES];
+ Float_t row_cg [MAXSTATES];
+ Float_t column_cg [MAXSTATES];
+ Int_t cerror;
+ Int_t rerror;
+ Int_t colerror;
+ UInt_t shapecode [MAXSTATES];
+ UInt_t xmin [MAXSTATES];
+ UInt_t ymin [MAXSTATES];
+ UInt_t dx [MAXSTATES];
+ UInt_t dy [MAXSTATES];
+};
+
+//####################################################################
+#endif
\ No newline at end of file
--- /dev/null
+// #include "config.h"
+
+#include "TApplication.h"
+#include <sstream>
+#include <TString.h>
+#include "unpack.h"
+#include "TROOT.h"
+
+//####################################################################
+
+
+int main(int argc, char** argv) {
+
+ gROOT->Reset();
+
+ TString Dir;
+ TString File;
+ UInt_t Start;
+ UInt_t Entries;
+
+ stringstream strValue;
+
+ if (argc>1) { Dir = argv[1]; }
+ if (argc>2) { File = argv[2]; }
+ if (argc>3) { strValue<<argv[3]; strValue >> Start; strValue.clear(); } else { Start =0; }
+ if (argc>4) { strValue<<argv[4]; strValue >> Entries; strValue.clear(); } else { Entries =0; }
+
+// TApplication theApp("Analysis", &argc, argv);
+ unpacker* a=new unpacker(Dir, File);
+ a->scanFrames(Start,Entries);
+ delete a;
+// theApp.Run();
+
+// TString Dir = "/local/samir/data/normalmode/";
+// TString File = "data.hld";
+// TString mode = "0";
+// UInt_t Start = 0;
+// UInt_t Entries = 0;
+}
\ No newline at end of file
--- /dev/null
+#ifndef UNPACK_H
+#define UNPACK_H
+
+#include "TMath.h"
+#include "unpack.h"
+// #include "help.h"
+#include "header.h"
+
+//####################################################################
+
+unpacker::unpacker(TString Dir, TString File) {
+
+ cout<<"--------------------------------"<<endl;
+ if( !checkFiles(Dir, File) ) { cerr << "ERROR: No Files found!" << endl << "-- Exit --" << endl; exit(-1); }
+}
+
+//####################################################################
+
+unpacker::unpacker() {
+}
+
+//####################################################################
+
+unpacker::~unpacker() {
+ if( fNrFiles>0 )
+ {
+ for(Int_t i=0; i<fNrFiles;i++)
+ {
+ fInn[i].close();
+ }
+ delete[] fFileLength;
+ delete[] fInn;
+ delete[] fRawData;
+#if SAVE > 0
+ cout<<"-------------"<<endl;
+ printf("\r Write data for Sensor:\n");
+ for(Int_t i=0; i<MAXSENSORS; i++)
+ {
+ if(fSensorpresent[i])
+ {
+ fFrameTree[i] ->Write("",TObject::kOverwrite);
+ printf(" ->%3i\n", i);
+ }
+ }
+ fOutputFile->Write();
+ fOutputFile->Save();
+ fOutputFile->Close();
+
+ cout<<"-------------"<<endl;
+ cout<<fOutFileName <<" created!"<<endl;
+ cout<<"-------------"<<endl;
+
+ TFile::Open(fOutFileName);
+// gROOT->ProcessLine("new TBrowser;");
+// gSystem->Exec(". ./runread.sh&");
+#endif
+ }
+
+// cout<<"--------------------------------"<<endl;
+// cout<<"' unpacker ' finished! "<<endl;
+// cout<<"--------------------------------"<<endl;
+
+}
+
+//####################################################################
+
+void unpacker::scanFrames(UInt_t Start, UInt_t Entries) {
+
+// cout<<"-----------"<<endl;
+ gSystem->Exec("date +%H:%M:%S");
+ cout<<"-----------"<<endl;
+
+// Progress bar
+ ULong_t progress_tmp=-1;
+ Float_t progress;
+
+// Variables
+ UInt_t eventSize;
+ UInt_t nrSensors;
+ UInt_t eventnrcount = 0;
+
+//---------------------------------- Tree branches -----------------------------------------
+ frameInfo *FrameInfo = new frameInfo[MAXSENSORS];
+ for(Int_t ii=0; ii<MAXSENSORS; ii++)
+ {
+ fFrameTree[ii]->Branch("event" , &FrameInfo[ii].event , "event/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("framenr" , &FrameInfo[ii].framenr , "framenr/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("sensorid" , &FrameInfo[ii].sensorid , "sensorid/s" , BUFSIZE);
+ fFrameTree[ii]->Branch("status" , &FrameInfo[ii].status , "status/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("datalength" , &FrameInfo[ii].datalength , "datalength/s" , BUFSIZE);
+ fFrameTree[ii]->Branch("datalengthcount" , &FrameInfo[ii].datalengthcount, "datalengthcount/s" , BUFSIZE);
+ fFrameTree[ii]->Branch("states" , &FrameInfo[ii].states , "states/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("clusters" , &FrameInfo[ii].clusters , "clusters/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("timestamp" , &FrameInfo[ii].timestamp , "timestamp/l" , BUFSIZE);
+ fFrameTree[ii]->Branch("debug" , &FrameInfo[ii].debug , "debug/l" , BUFSIZE);
+ fFrameTree[ii]->Branch("row" , FrameInfo[ii].row , "row [states]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("column" , FrameInfo[ii].column , "column [states]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("nextn" , FrameInfo[ii].nextn , "nextn [states]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("overflow" , FrameInfo[ii].overflow , "overflow [states]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("pixelsPerCluster" , FrameInfo[ii].pixelsPerCluster, "pixelsPerCluster [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("statesPerCluster" , FrameInfo[ii].statesPerCluster, "statesPerCluster [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("row_cg" , FrameInfo[ii].row_cg , "row_cg [clusters]/F" , BUFSIZE);
+ fFrameTree[ii]->Branch("column_cg" , FrameInfo[ii].column_cg , "column_cg [clusters]/F" , BUFSIZE);
+ fFrameTree[ii]->Branch("cerror" , &FrameInfo[ii].cerror , "cerror/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("rerror" , &FrameInfo[ii].rerror , "rerror/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("colerror" , &FrameInfo[ii].colerror , "colerror/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("shapecode" , FrameInfo[ii].shapecode , "shapecode [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("xmin" , FrameInfo[ii].xmin , "xmin [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("ymin" , FrameInfo[ii].ymin , "ymin [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("dx" , FrameInfo[ii].dx , "dx [clusters]/i" , BUFSIZE);
+ fFrameTree[ii]->Branch("dy" , FrameInfo[ii].dy , "dy [clusters]/i" , BUFSIZE);
+ }
+//------------------------------------------------------------------------------------------
+// Limit data loop:
+ ULong64_t start = 0;
+ ULong64_t entries = 0;
+
+ if( Entries > fEvents ) { entries = fEvents; printf("Changed 'Number Events' to: %lld\n", entries); }
+ else { entries = Entries; }
+ if( Start+entries > fEvents ) { start = fEvents-entries; printf("Changed 'First Event' to: %lld\n", start ); }
+ else { start = Start; }
+ if( Start==0 && Entries==0 ) { start = 0; entries = fEvents-1; }
+ eventnrcount = (UInt_t) start;
+//------------------------------------------------------------------------------------------
+// Data loop:
+ if( fNrFiles>0 )
+ {
+ ULong64_t i;
+ for(i=0;i<start;i++)
+ {
+ eventSize = getNextEvent();
+ }
+ for(;i<start+entries;i++)
+ {
+ eventSize = getNextEvent();
+ nrSensors = getFrameInfo(eventSize, eventnrcount, FrameInfo);
+ if( nrSensors!=fNrSensors )
+ {
+ printf("\r Data inconsistent -> Wrong number of sensors (%3u --> %3u )! Frame number %10u skipped! File Number: %10u\n", fNrSensors, nrSensors, eventnrcount, fCurrFile);
+ }
+ else
+ {
+// getStates(fNrSensors, FrameInfo, fFrameTree); // Note: nrSensors may sometimes != fNrSensors
+ getClusters(fNrSensors, FrameInfo, fFrameTree); // Note: nrSensors may sometimes != fNrSensors
+ }
+ eventnrcount++;
+//------------------------------------------------------------------------------------------
+// Progress bar
+ progress = (Int_t)(((i-start)*100)/(entries-1)*10);
+ if(progress!=progress_tmp) { print_progress( (((i-start)*100.)/(entries-1)) ); progress_tmp=progress;}
+//------------------------------------------------------------------------------------------
+ }
+ }
+
+ cout<<"-----------"<<endl;
+ gSystem->Exec("date +%H:%M:%S");
+// cout<<"-----------"<<endl;
+}
+
+//####################################################################
+//####################################################################
+
+bool unpacker::checkFiles(TString Dir, TString File) {
+
+ // Check whether files are there
+ // Open all files
+ // Check files: Size
+
+ fNrFiles = 0;
+
+ char *dir = (char*) Dir.Data();
+ char *file = (char*) File.Data();
+ char ending[] = ".hld";
+ bool has_ending = true;
+ bool has_ex = false;
+ uint len;
+
+ ifstream inn;
+ TString FileName;
+ TString FileNames[100];
+ ULong_t countEvents;
+ fEvents = 0;
+
+ ULong_t End;
+
+// Check naming of 'Dir' and 'File'
+ if( dir[strlen(dir)-1] != '/' )
+ {
+ strcat(dir, "/" );
+ Dir = dir;
+ }
+
+ len = strlen(file);
+
+ for(uint i=0;i<4;i++)
+ {
+ if( file [len-4+i] != ending[i] )
+ {
+ has_ending = false;
+ }
+ }
+
+ for(uint i=0;i<len;i++)
+ {
+ if( file [i] == '.' )
+ {
+ has_ex = true;
+ }
+ }
+
+ if(has_ex && !has_ending)
+ {
+ cerr<<"WARNING: Try to read file(s) with unknown filename extension: "<<file<<endl;
+ }
+
+ if( has_ending )
+ {
+ file[len-4] = '\0';
+ File = file;
+ }
+
+// Try to find 'Dir' and 'File'
+ DIR* dpdf = opendir(Dir);
+ if (dpdf == NULL)
+ {
+ cerr<<"ERROR: Directory not found: "<<Dir<<endl;
+ return 0;
+ }
+
+ do
+ {
+ FileName = Dir + File + Form("_%i",fNrFiles) + (TString)ending;
+ inn.open(FileName, ios::binary);
+
+ if(!inn.good())
+ {
+ inn.close();
+ break;
+ }
+
+ inn.close();
+ FileNames[fNrFiles] = FileName;
+ fNrFiles++;
+
+ }while(inn.good());
+
+ if( fNrFiles==0 )
+ {
+ FileName = Dir + File + (TString)ending;
+ inn.open(FileName, ios::binary);
+
+ if(!inn.good())
+ {
+ inn.close();
+ return 0;
+ }
+
+ inn.close();
+ FileNames[fNrFiles] = FileName;
+ fNrFiles++;
+ }
+
+// Link File(s)
+ if( fNrFiles>0 )
+ {
+ len = strlen(FileNames[fNrFiles-1]);
+
+ fRawData = new char[CONTENTIN];
+ fInn = new ifstream [fNrFiles];
+ fFileLength = new ULong_t [fNrFiles];
+ fFileEvents = new ULong_t [fNrFiles];
+ fLength = 0;
+
+// cout<<"--------"<<endl;
+ cout<<"Linked file(s):"<<endl;
+
+ for(int i=0; i<fNrFiles;i++)
+ {
+ fInn[i].open(FileNames[i], ios::binary);
+ fInn[i].seekg(0, ios::end);
+ End = fInn[i].tellg();
+ fFileLength[i] = End;
+ fLength+=End;
+ fInn[i].seekg(0, ios::beg);
+
+
+ cout<<" "<<setw(len+1)<<left<< FileNames[i];
+
+ if (End/(TMath::Power(2,40)) >= 1) { printf(" - %5.2f TBytes", (Float_t)(End/(TMath::Power(2,40))) ); }
+ else if (End/(TMath::Power(2,30)) >= 1) { printf(" - %5.2f GBytes", (Float_t)(End/(TMath::Power(2,30))) ); }
+ else if (End/(TMath::Power(2,20)) >= 1) { printf(" - %5.2f MBytes", (Float_t)(End/(TMath::Power(2,20))) ); }
+ else if (End/(TMath::Power(2,10)) >= 1) { printf(" - %5.2f kBytes", (Float_t)(End/(TMath::Power(2,10))) ); }
+ else { printf(" - %5.2f Bytes ", (Float_t)(End)); }
+
+// Count events
+ countEvents = 0;
+ for(ULong_t j=0; j<End; )
+ {
+ fInn[i].seekg(j, ios::beg);
+ fInn[i].read(fRawData, 4);
+
+ j += littleEndian32(fRawData,0);
+ if(j%8!=0) { j=j+8-j%8; }
+ countEvents++;
+ }
+
+ fFileEvents[i] = countEvents;
+ fEvents+=countEvents;
+
+// Output events:
+ printf(" --> %8ld Events\n", countEvents);
+ }
+// cout<<"--------"<<endl;
+ }
+
+ fCurrFile = 0;
+ fCurrPosInFile = 0;
+
+// Determine number of sensors:
+ for(Int_t i=0; i<fNrFiles;i++)
+ {
+ fInn[i].seekg(0, ios::beg);
+ }
+ Float_t NrSensors_tmp = 0;
+ UInt_t eventSize;
+ UInt_t eventnrcount = 0;
+ frameInfo *FrameInfo = new frameInfo[MAXSENSORS];
+ for(Int_t i=0; i<10;i++)
+ {
+ eventSize = getNextEvent();
+ NrSensors_tmp += getFrameInfo(eventSize, eventnrcount, FrameInfo);
+ }
+ NrSensors_tmp = NrSensors_tmp/10.+0.5;
+ fNrSensors = (UInt_t)NrSensors_tmp;
+ cout<<"--------------------------------"<<endl;
+ printf(" --> %u sensors found! <-- \n",fNrSensors);
+ cout<<"--------------------------------"<<endl;
+
+ for(Int_t i=0; i<fNrFiles;i++)
+ {
+ fInn[i].seekg(0, ios::beg);
+ }
+
+// ----------------
+ fOutFileName = Dir + File + Form("_DST.root");
+
+ fCurrFile = 0;
+ fCurrPosInFile = 0;
+// ----------------
+#if SAVE > 0
+ fOutputFile = new TFile (fOutFileName,"RECREATE");
+#endif
+ for(Int_t i=0; i<MAXSENSORS; i++)
+ {
+ fFrameTree[i] = new TTree(Form("sensor%i",i), Form("sensor%i",i));
+ fSensorpresent[i] = false;
+ }
+// ----------------
+
+ return fNrFiles;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextPseudoEvent() {
+
+ fCurrPosInFile = fInn[fCurrFile].tellg();
+
+ if( fCurrPosInFile >= fFileLength[fCurrFile] )
+ {
+ if( fCurrFile+1<fNrFiles )
+ {
+ fCurrFile++;
+ fCurrPosInFile = 0;
+ }
+ else
+ {
+ return 0;
+ }
+ }
+
+ fInn[fCurrFile].read(fRawData, 4);
+
+ UInt_t eventSize = littleEndian32(fRawData,0);
+ UInt_t padding;
+
+ if (eventSize%8!=0) { padding = 8 - eventSize%8; }
+ else { padding = 0; }
+
+ UInt_t eventRead = eventSize + padding - 4;
+ fInn[fCurrFile].read(fRawData, eventRead);
+
+ return eventSize;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextEvent() {
+
+ UInt_t eventSize = getNextPseudoEvent();
+
+ if(eventSize>32)
+ {
+ return eventSize;
+ }
+ else
+ {
+ return getNextPseudoEvent();
+ }
+}
+
+//####################################################################
+
+UInt_t unpacker::getFrameInfo( UInt_t EventSize, UInt_t FrameNrCount, frameInfo*& FrameInfo ) {
+
+ UInt_t eventSize = EventSize;
+ UInt_t subEventSize;
+ UInt_t subSubEventLength;
+
+ UInt_t posInEvent = 28; // Skip event header info - Length event header (32 bytes - 4 bytes (event size))
+ UInt_t posInSubEvent;
+ UInt_t posInSubSubEvent;
+
+ UInt_t count_subsubevents = 0; // Number of different sensors
+
+// UInt_t hheader;
+ UInt_t sensorid;
+ UInt_t status;
+ ULong_t debug=0;
+ ULong_t timestamp;
+ UInt_t header;
+ UInt_t framenr;
+ UInt_t datalength;
+ UInt_t framelength;
+ UInt_t trailer;
+ UInt_t padding=0;
+// UInt_t HDRver;
+ UInt_t DATAver;
+// UInt_t HDRsize;
+// UInt_t mtimestamp;
+//------------------------------------------------------------------------------------------
+//---------------------------------- Event loop --------------------------------------------
+//------------------------------------------------------------------------------------------
+ while( posInEvent+4 < eventSize ) // 4 bytes were skipped
+ {
+ subEventSize = bigEndian32(fRawData,posInEvent); // including itself
+ posInEvent += 16; // Skip sub-event header info - Length sub-event header (16 bytes) (trailer: (8 bytes))
+ posInSubEvent = 16; // -- " --
+//---------------------------------- Sub-event loop ------------------------------------
+//------------------------------------------------------------------------------------------
+ while( posInSubEvent+8 < subEventSize && posInEvent+4 < eventSize ) // 8 bytes -> trailer size
+ {
+ subSubEventLength = 4*bigEndian16(fRawData,posInEvent); // Number of 8 bit words in subsubevent excluding itself
+ posInEvent+=4; // Skip sub-sub-event header info
+ posInSubEvent+=4;
+ posInSubSubEvent=0;
+//---------------------------------- Individual sub-sub-events -----------------------------
+ while( posInSubSubEvent < subSubEventLength && posInSubEvent+8 < subEventSize && posInEvent+4 < eventSize )
+ {
+//---------------------------------- Check individual sub-sub-events -----------------------
+ /*HDRver = bigEndian8 (fRawData,posInEvent);*/ posInEvent+=1; posInSubEvent+=1; posInSubSubEvent+=1;
+ DATAver = bigEndian8 (fRawData,posInEvent); posInEvent+=1; posInSubEvent+=1; posInSubSubEvent+=1;
+ posInEvent+=1; posInSubEvent+=1; posInSubSubEvent+=1;
+ /*HDRsize = bigEndian8 (fRawData,posInEvent);*/ posInEvent+=1; posInSubEvent+=1; posInSubSubEvent+=1;
+ /*mtimestamp = bigEndian32(fRawData,posInEvent);*/ posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4;
+ framelength = bigEndian16(fRawData,posInEvent)*4; posInEvent+=2; posInSubEvent+=2; posInSubSubEvent+=2;
+ sensorid = bigEndian16(fRawData,posInEvent)%256; posInEvent+=2; posInSubEvent+=2; posInSubSubEvent+=2; // ID 00XX
+ status = bigEndian32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // status XXXX XXXX
+
+ if( status==0xF000000F ) // --> Frame ok!
+ {
+ if( DATAver!=0x01 ) { printf("\rData Version not identified: %#04x \n", DATAver ); }
+
+ debug = bigEndian32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Debug XXXX XXXX
+ timestamp = bigEndian32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Timestamp XXXX XXXX
+
+ header = bigEndian32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Header 5555 5555
+ framenr = bigEndian32Switch16(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Framenumber XXXX XXXX
+ datalength = datalength32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Datalength XXXX XXXX (0+1)
+
+ FrameInfo[count_subsubevents].event = FrameNrCount;
+ FrameInfo[count_subsubevents].posInEvent = posInEvent;
+ FrameInfo[count_subsubevents].sensorid = sensorid;
+ FrameInfo[count_subsubevents].status = status;
+ FrameInfo[count_subsubevents].debug = debug;
+ FrameInfo[count_subsubevents].timestamp = timestamp;
+ FrameInfo[count_subsubevents].framenr = framenr;
+ FrameInfo[count_subsubevents].datalength = datalength;
+
+ posInEvent += 2*datalength;
+ posInSubEvent += 2*datalength;
+ posInSubSubEvent += 2*datalength;
+
+ trailer = bigEndian32(fRawData,posInEvent); posInEvent+=4; posInSubEvent+=4; posInSubSubEvent+=4; // Trailer 8001 8001
+ if( header != 0x55555555 || trailer != 0x80018001 )
+ {
+ printf("\rHeader/Trailer: 0x55555555/0x55555555 != %#010x/%#010x ( File Number %3i - Event Number %10i ) \n", header, trailer, fCurrFile, FrameNrCount );
+
+ posInEvent -= (2*datalength+4);
+ posInSubEvent -= (2*datalength+4);
+ posInSubSubEvent -= (2*datalength+4);
+
+ posInEvent += (framelength+12-36);
+ posInSubEvent += (framelength+12-36);
+ posInSubSubEvent += (framelength+12-36);
+
+ datalength = 0;
+ }
+ FrameInfo[count_subsubevents].datalengthcount = datalength;
+ }
+ else // --> Frame not recorded / Header missing / Error
+ {
+// printf("\rStatus: 0xF000000F != %#010x ( File Number %3i - Event Number %10i ) \n", status, fCurrFile, FrameNrCount );
+
+ FrameInfo[count_subsubevents].event = FrameNrCount;
+ FrameInfo[count_subsubevents].posInEvent = posInEvent;
+ FrameInfo[count_subsubevents].sensorid = sensorid;
+ FrameInfo[count_subsubevents].status = status;
+ FrameInfo[count_subsubevents].debug = debug;
+ FrameInfo[count_subsubevents].timestamp = 0;
+ FrameInfo[count_subsubevents].framenr = 0;
+ FrameInfo[count_subsubevents].datalength = 0;
+ FrameInfo[count_subsubevents].datalengthcount = 0;
+
+ if( status!=0 )
+ {
+ posInSubSubEvent += (framelength+12-16);
+ posInSubEvent += (framelength+12-16);
+ posInEvent += (framelength+12-16);
+ }
+ }
+ count_subsubevents++;
+ }
+//---------------------------------- END: Sub-sub-event loop -------------------------------
+ }
+// printf(" TRAILER: %#010x ", bigEndian32(fRawData,posInEvent) );
+ posInEvent+=8; // Trailer of Sub-event
+ if (subEventSize%8!=0) { padding = 8 - subEventSize%8;}
+ else { padding = 0; }
+ posInEvent+=padding;
+//---------------------------------- END: Sub-event loop -----------------------------------
+//------------------------------------------------------------------------------------------
+ if( bigEndian32(fRawData,eventSize-12) != 0x00015555 )
+ {
+ printf("\rDATA ERROR ( File Number %3i - Event Number %10i ) -> TRAILER: 0x00015555 != %#010x\n",fCurrFile, FrameNrCount, bigEndian32(fRawData,eventSize-12));
+ }
+ }
+//------------------------------------------------------------------------------------------
+//---------------------------------- END: Event loop ---------------------------------------
+//------------------------------------------------------------------------------------------
+ return count_subsubevents;
+}
+
+//####################################################################
+
+void unpacker::getStates( UInt_t NrSensors, frameInfo*& FrameInfo, TTree** FrameTree ) {
+
+ UInt_t posInEvent;
+ UInt_t count_datalength = 0;
+ UInt_t count_statesPerLine = 0;
+ UInt_t addressLine;
+ UInt_t statesPerLine;
+ UInt_t overflow;
+
+ for(UInt_t i=0; i<NrSensors;i++)
+ {
+
+ posInEvent = FrameInfo[i].posInEvent;
+ FrameInfo[i].states = 0;
+ count_datalength = 0;
+
+ while( count_datalength < FrameInfo[i].datalengthcount )
+ {
+ addressLine = addressLine16 (fRawData,posInEvent);
+ statesPerLine = numberStatesInLine16 (fRawData,posInEvent);
+ overflow = overflowInLine16 (fRawData,posInEvent);
+
+ count_datalength ++;
+ posInEvent +=2;
+ count_statesPerLine = 0;
+
+ while( count_statesPerLine < statesPerLine && count_datalength < FrameInfo[i].datalengthcount )
+ {
+ FrameInfo[i].row [ FrameInfo[i].states ] = addressLine;
+ FrameInfo[i].column [ FrameInfo[i].states ] = addressColumn16 (fRawData,posInEvent);
+ FrameInfo[i].nextn [ FrameInfo[i].states ] = numberFiredPixelsInColumn16 (fRawData,posInEvent);
+ FrameInfo[i].overflow [ FrameInfo[i].states ] = overflow;
+
+ FrameInfo[i].states ++;
+ count_statesPerLine ++;
+ count_datalength ++;
+ posInEvent +=2;
+ }
+ }
+ FrameTree[i]->Fill();
+ fSensorpresent[i] = true;
+ }
+}
+
+//####################################################################
+
+void unpacker::getClusters( UInt_t NrSensors, frameInfo*& FrameInfo, TTree** FrameTree ) {
+
+ UInt_t posInEvent;
+ UInt_t count_datalength = 0;
+ UInt_t count_statesPerLine = 0;
+ UInt_t addressLine;
+ UInt_t statesPerLine;
+ UInt_t overflow;
+
+ vector< vector<UInt_t> > fClusterVector; //store the cluster information in fClusterVector ,including the index of state
+ vector< UInt_t > fCluster;
+ vector< UInt_t > fTempCluster;
+ vector< UInt_t > compareNumOfClusterVector;//store the num of state needed to compare in each cluster..
+ vector< UInt_t > newCompareNumOfClusterVector;
+
+
+ UInt_t priorRow = 0;
+ UInt_t newRow = 0;
+ UInt_t priorColumn = 0;
+ UInt_t priorNumPixel = 0;
+ UInt_t newColumn = 0;
+ UInt_t newNumPixel = 0;
+ UInt_t indexState = 0;
+ UInt_t clusterNum = 0;
+ UInt_t clusterSize = 0;
+ UInt_t compareNumOfCluster = 0;
+ UInt_t registerFlag =0;
+
+ UInt_t numActivePixels = 0;
+
+ Float_t sumX = 0.0;
+ Float_t sumY = 0.0;
+
+ for(UInt_t i=0; i<NrSensors;i++)
+ {
+
+
+ posInEvent = FrameInfo[i].posInEvent;
+ FrameInfo[i].states = 0;
+ FrameInfo[i].colerror = 0;
+ FrameInfo[i].rerror = 0;
+ FrameInfo[i].clusters = 0;
+ fClusterVector.clear();
+ newCompareNumOfClusterVector.clear();
+ compareNumOfClusterVector.clear();
+ priorColumn = 0;
+ newColumn = 0;
+ newNumPixel = 0;
+ priorNumPixel = 0;
+ priorRow = 0;
+ newRow = 0;
+ count_datalength = 0;
+ while( count_datalength < FrameInfo[i].datalengthcount )
+ {
+ addressLine = addressLine16 (fRawData,posInEvent);
+ statesPerLine = numberStatesInLine16 (fRawData,posInEvent);
+ overflow = overflowInLine16 (fRawData,posInEvent);
+ count_datalength++;
+ posInEvent += 2;
+ count_statesPerLine = 0;
+
+ priorRow = newRow;
+ registerFlag = 0;
+ compareNumOfClusterVector.clear();
+ compareNumOfClusterVector = newCompareNumOfClusterVector;
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ newCompareNumOfClusterVector.clear();
+ newRow = addressLine;
+ if( ((Int_t)(newRow - priorRow) < 0) && count_datalength!= FrameInfo[i].datalengthcount-1)
+ {
+ FrameInfo[i].rerror = 1;
+ }
+ else
+ {
+ if( statesPerLine>9 )
+ FrameInfo[i].colerror = 1;
+ else
+ {
+ if(((newRow - priorRow) > 1)||(compareNumOfCluster == 0))
+ {
+ while( count_statesPerLine<statesPerLine && count_datalength<FrameInfo[i].datalengthcount )
+ {
+
+ FrameInfo[i].row [ FrameInfo[i].states ] = addressLine;
+ FrameInfo[i].column [ FrameInfo[i].states ] = addressColumn16 (fRawData,posInEvent);
+ FrameInfo[i].nextn [ FrameInfo[i].states ] = numberFiredPixelsInColumn16 (fRawData,posInEvent);
+ FrameInfo[i].overflow [ FrameInfo[i].states ] = overflow;
+
+ clusterNum = fClusterVector.size();
+ priorColumn = newColumn;
+ newColumn = FrameInfo[i].column[ FrameInfo[i].states ];
+ newNumPixel = FrameInfo[i].nextn[FrameInfo[i].states];
+
+
+ if( count_statesPerLine == 0 )
+ priorNumPixel = 0;
+ else
+ priorNumPixel = FrameInfo[i].nextn[FrameInfo[i].states -1];
+ if (((newColumn - (priorColumn + priorNumPixel)) > 1)||(count_statesPerLine == 0))
+ {
+ fClusterVector.push_back(fCluster);
+ fClusterVector[clusterNum].push_back(clusterNum);
+ fClusterVector[clusterNum].push_back(FrameInfo[i].states);
+ newCompareNumOfClusterVector.push_back(1);
+
+ }
+ else
+ if((newColumn - (priorColumn + priorNumPixel)) == 1)
+ {
+ fClusterVector[clusterNum-1].push_back(FrameInfo[i].states);
+ Int_t iNum = newCompareNumOfClusterVector.size();
+ newCompareNumOfClusterVector[iNum-1]++;
+ }
+ else if(count_datalength != FrameInfo[i].datalengthcount -1)
+ {
+ FrameInfo[i].colerror = 1;
+ }
+ FrameInfo[i].states++;
+ count_statesPerLine++;
+ count_datalength++;
+ posInEvent+=2;
+
+ }
+ }
+ else if((newRow - priorRow) == 1)
+ {
+ while( count_statesPerLine<statesPerLine && count_datalength<FrameInfo[i].datalengthcount )
+ {
+
+ FrameInfo[i].row [ FrameInfo[i].states ] = addressLine;
+ FrameInfo[i].column [ FrameInfo[i].states ] = addressColumn16 (fRawData,posInEvent);
+ FrameInfo[i].nextn [ FrameInfo[i].states ] = numberFiredPixelsInColumn16 (fRawData,posInEvent);
+ FrameInfo[i].overflow [ FrameInfo[i].states ] = overflow;
+
+ newColumn = FrameInfo[i].column[FrameInfo[i].states];
+ newNumPixel = FrameInfo[i].nextn[FrameInfo[i].states];
+ clusterNum = fClusterVector.size();
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ Int_t stop = 0;
+ registerFlag = 0;
+
+ for(Int_t j=compareNumOfCluster;j>stop;j--)
+ {
+ clusterSize = fClusterVector[clusterNum - j].size();
+ for(Int_t m=compareNumOfClusterVector[compareNumOfCluster - j];m>0;m--)
+ {
+ indexState = fClusterVector[clusterNum - j][clusterSize - m];
+ priorColumn = FrameInfo[i].column[indexState];
+ priorNumPixel = FrameInfo[i].nextn[indexState];
+
+ if (((newColumn + newNumPixel + 1) >= priorColumn) && (newColumn <= (priorColumn + priorNumPixel +1)))
+ {
+ stop++;
+ m=0;
+ Int_t sub = 0;
+ if(registerFlag == 1)
+ {
+ indexState = fClusterVector[clusterNum - j][clusterSize-1];
+ Int_t iNum = newCompareNumOfClusterVector.size();
+ Int_t temp1 = newCompareNumOfClusterVector[iNum-1];
+ Int_t temp2 = 0;
+ if(FrameInfo[i].row[indexState] == addressLine)
+ {
+ Int_t sum=0;
+ sub = FrameInfo[i].states - indexState;
+ if(sub == 1)
+ {
+ temp2 = newCompareNumOfClusterVector[iNum-2];
+ newCompareNumOfClusterVector[iNum-2] = temp2 + temp1;
+ newCompareNumOfClusterVector.erase(newCompareNumOfClusterVector.begin()+iNum-1);
+ }
+ else
+ {
+ for(Int_t jj=iNum-2;jj>0;jj--)
+ {
+ sum = sum + newCompareNumOfClusterVector[jj];
+ if(sum == sub -1)
+ {
+ newCompareNumOfClusterVector[jj-1]++;
+ Int_t temp5 = newCompareNumOfClusterVector[jj-1];
+ for(Int_t ll = jj;ll<iNum-1;ll++)
+ newCompareNumOfClusterVector[ll-1] = newCompareNumOfClusterVector[ll];
+ newCompareNumOfClusterVector[iNum-2] = temp5;
+ newCompareNumOfClusterVector.erase(newCompareNumOfClusterVector.begin()+iNum-1);
+ jj = 0;
+ }
+ }
+ }
+ }
+ vector< Int_t > bb;
+ bb.clear();
+ Int_t cc = fClusterVector[clusterNum-1].size();
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ if(FrameInfo[i].row[indexState] != addressLine)
+ {
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ Int_t temp3 = compareNumOfClusterVector[compareNumOfCluster - 1];
+ for(UInt_t aa = temp3;aa>0;aa--)
+ {
+ bb.push_back(fClusterVector[clusterNum-1][cc-aa]);
+ }
+ for(Int_t aa = temp3;aa>0;aa--)
+ fClusterVector[clusterNum-1].pop_back();
+ for(UInt_t aa = 1;aa<fClusterVector[clusterNum - j].size();aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(fClusterVector[clusterNum - j][aa]);
+ }
+ for(UInt_t aa = 0;aa<bb.size();aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(bb[aa]);
+ }
+ }
+ else
+ {
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ Int_t temp4 = compareNumOfClusterVector[compareNumOfCluster - 1];
+ for(UInt_t aa = temp4;aa>0;aa--)
+ {
+ bb.push_back(fClusterVector[clusterNum-1][cc-aa]);
+ }
+ for(UInt_t aa = temp4;aa>0;aa--)
+ fClusterVector[clusterNum-1].pop_back();
+ for(UInt_t aa = 1;aa<(fClusterVector[clusterNum - j].size()-temp2);aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(fClusterVector[clusterNum - j][aa]);
+ }
+ for(UInt_t aa = 0;aa<(bb.size()-temp1);aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(bb[aa]);
+ }
+ for(UInt_t aa = (fClusterVector[clusterNum - j].size()-temp2);aa<fClusterVector[clusterNum - j].size();aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(fClusterVector[clusterNum - j][aa]);
+ }
+ for(UInt_t aa = (bb.size()-temp1);aa<bb.size();aa++)
+ {
+ fClusterVector[clusterNum-1].push_back(bb[aa]);
+ }
+ }
+ bb.clear();
+ Int_t ii = compareNumOfClusterVector.size();
+ compareNumOfClusterVector[ii-1] = compareNumOfClusterVector[ii-1] + compareNumOfClusterVector[compareNumOfCluster - j];
+ fClusterVector.erase(fClusterVector.begin()+clusterNum - j);
+ compareNumOfClusterVector.erase(compareNumOfClusterVector.begin()+compareNumOfCluster - j);
+ clusterNum = fClusterVector.size();
+ }
+ else
+ {
+ registerFlag = 1;
+ indexState = fClusterVector[clusterNum - j][clusterSize-1];
+ if(FrameInfo[i].row[indexState] == addressLine)
+ {
+ Int_t iNum = newCompareNumOfClusterVector.size();
+ Int_t sum=0;
+ sub = FrameInfo[i].states - indexState;
+ if(sub == 1)
+ {
+ newCompareNumOfClusterVector[iNum-1]++;
+ }
+ else
+ {
+ for(Int_t jj=iNum-1;jj>0;jj--)
+ {
+ sum = sum + newCompareNumOfClusterVector[jj];
+ if(sum == sub -1)
+ {
+ newCompareNumOfClusterVector[jj-1]++;
+ Int_t temp5 = newCompareNumOfClusterVector[jj-1];
+ for(Int_t ll = jj;ll<iNum;ll++)
+ newCompareNumOfClusterVector[ll-1] = newCompareNumOfClusterVector[ll];
+ newCompareNumOfClusterVector[iNum-1] = temp5;
+ jj = 0;
+ }
+ }
+ }
+ }
+ else
+ {
+ newCompareNumOfClusterVector.push_back(1);
+ }
+ fClusterVector[clusterNum - j].push_back(FrameInfo[i].states);
+ compareNumOfClusterVector[compareNumOfCluster - j]++;
+ fTempCluster.clear();
+ fTempCluster = fClusterVector[clusterNum - j];
+ Int_t iNumOfLoop = fClusterVector.size();
+ for(Int_t t = clusterNum - j;t<(iNumOfLoop-1);t++)
+ {
+ fClusterVector[t].clear();
+ fClusterVector[t] = fClusterVector[t+1];
+ }
+ fClusterVector[iNumOfLoop -1] = fTempCluster;
+ Int_t iTempNum = compareNumOfClusterVector[compareNumOfCluster - j];
+ iNumOfLoop = compareNumOfClusterVector.size();
+ for(Int_t t = compareNumOfCluster - j;t<(iNumOfLoop-1);t++)
+ {
+ compareNumOfClusterVector[t] = compareNumOfClusterVector[t+1];
+ }
+ compareNumOfClusterVector[iNumOfLoop -1] = iTempNum;
+ j++;
+ }
+ }
+ else if((m == 1) && (j == 1)&&(registerFlag == 0))
+ {
+ fClusterVector.push_back(fCluster);
+ fClusterVector[clusterNum].push_back(clusterNum);
+ fClusterVector[clusterNum].push_back(FrameInfo[i].states);
+ newCompareNumOfClusterVector.push_back(1);
+ compareNumOfClusterVector.push_back(1);
+ clusterNum = fClusterVector.size();
+ }
+ compareNumOfCluster = compareNumOfClusterVector.size();
+ }
+ }
+ FrameInfo[i].states++;
+ count_statesPerLine++;
+ count_datalength++;
+ posInEvent+=2;
+ }
+ }
+ }
+ }
+ }
+
+
+ clusterNum = fClusterVector.size();
+ FrameInfo[i].clusters = clusterNum;
+ vector < UInt_t > xVector;
+ vector < UInt_t > yVector;
+ vector < UInt_t > nPixel;
+ UInt_t code = 0;
+ sort(&fClusterVector[0],&fClusterVector[clusterNum]);
+ for(UInt_t num1=0;num1<clusterNum;num1++)
+ {
+ yVector.clear();
+ xVector.clear();
+ sumX = 0;
+ sumY = 0;
+ code = 0;
+ numActivePixels = 0;
+
+ clusterSize = fClusterVector[num1].size();
+ if(clusterSize > 10)
+ FrameInfo[i].cerror = 1;
+ else
+ FrameInfo[i].cerror = 0;
+ for(UInt_t num2=1;num2<clusterSize;num2++)
+ {
+ indexState = fClusterVector[num1][num2];
+ for(UInt_t num3 = 0;num3<FrameInfo[i].nextn[indexState]+1;num3++)
+ {
+ xVector.push_back(FrameInfo[i].column[indexState]+num3) ;
+ yVector.push_back(FrameInfo[i].row[indexState]);
+ }
+ }
+
+ numActivePixels = xVector.size();
+ UInt_t minX = 1152;
+ UInt_t minY = 576;
+ UInt_t maxX = 0;
+ UInt_t maxY = 0;
+
+ for(UInt_t m=0;m<numActivePixels;m++)
+ {
+ if(minX > xVector[m])
+ minX = xVector[m];
+ if(minY > yVector[m])
+ minY = yVector[m];
+ if(maxX < xVector[m])
+ maxX = xVector[m];
+ if(maxY < yVector[m])
+ maxY = yVector[m];
+ }
+ UInt_t dX = maxX - minX + 1;
+ UInt_t dY = maxY - minY + 1;
+ //coding
+ if((dX<16) && (dY<16) && (dX*dY<28))
+ {
+ for(UInt_t mm=0;mm<numActivePixels;mm++)
+ {
+ sumX = sumX + xVector[mm];
+ sumY = sumY + yVector[mm];
+ code = code + pow(2.0,dX*(yVector[mm]-minY))*pow(2.0,(xVector[mm]-minX));
+ }
+ code = code*16+dX;
+ }
+ else
+ code = 0;
+ FrameInfo[i].column_cg[num1] = (Float_t)sumX/numActivePixels;
+ FrameInfo[i].row_cg[num1] = (Float_t)sumY/numActivePixels;
+ FrameInfo[i].statesPerCluster[num1] = (clusterSize-1);
+ FrameInfo[i].pixelsPerCluster[num1] = numActivePixels;
+ FrameInfo[i].shapecode[num1] = code;//shapeCoder(code);
+ FrameInfo[i].xmin[num1] = minX;
+ FrameInfo[i].ymin[num1] = minY;
+ FrameInfo[i].dx[num1] = dX;
+ FrameInfo[i].dy[num1] = dY;
+ }
+ compareNumOfClusterVector.clear();
+ newCompareNumOfClusterVector.clear();
+ //plotClusters(10,FrameInfo[i].sensorid,12,FrameInfo[i].event,fClusterVector,FrameInfo[i].row,FrameInfo[i].column,FrameInfo[i].nextn);//plot all the cluster shape of one sensor one frame
+
+ FrameTree[i]->Fill();
+ fSensorpresent[i] = true;
+ fClusterVector.clear();
+ }
+}
+
+//####################################################################
+
+UInt_t unpacker::littleEndian32(char* RawData, Int_t Pos) {
+
+ return (Int_t)(UChar_t)RawData[3+Pos]*TMath::Power(2,24) + (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,16) + (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,8) + (Int_t)(UChar_t)RawData[0+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian32(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,24) +
+ (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,16) +
+ (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[3+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian32Switch16(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,24) +
+ (Int_t)(UChar_t)RawData[3+Pos]*TMath::Power(2,16) +
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[1+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian16(char* RawData, Int_t Pos) {
+
+ return (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,8) + (Int_t)(UChar_t)RawData[1+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian8(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::datalength32(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,0) +
+ (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[3+Pos]*TMath::Power(2,0);
+}
+
+//####################################################################
+
+UInt_t unpacker::numberStatesInLine16(char* RawData, Int_t Pos) {
+
+ return (( (Int_t)(UChar_t)RawData[1+Pos] >> 0) & 15);
+}
+
+//####################################################################
+
+UInt_t unpacker::addressLine16(char* RawData, Int_t Pos) {
+
+ return (( (Int_t)(UChar_t)RawData[0+Pos] >> 0) & 127)* 16 + (( (Int_t)(UChar_t)RawData[1+Pos] >> 4) & 15);
+}
+
+//####################################################################
+
+UInt_t unpacker::overflowInLine16(char* RawData, Int_t Pos) {
+
+ return (( (Int_t)(UChar_t)RawData[0+Pos] >> 7) & 1);
+}
+
+//####################################################################
+
+UInt_t unpacker::addressColumn16(char* RawData, Int_t Pos) {
+
+ return (( (Int_t)(UChar_t)RawData[0+Pos] >> 0) & 31)* 64 + (( (Int_t)(UChar_t)RawData[1+Pos] >> 2) & 63);
+}
+
+//####################################################################
+
+UInt_t unpacker::numberFiredPixelsInColumn16(char* RawData, Int_t Pos) {
+
+ return (( (Int_t)(UChar_t)RawData[1+Pos] >> 0) & 3);
+}
+
+//####################################################################
+
+Int_t unpacker::print_progress(Float_t ProgressInPercent) {
+
+ Int_t nCharacter = printf("\r [");
+ Int_t iBlock;
+
+ for( iBlock=0; iBlock<ProgressInPercent/2; ++iBlock )
+ {
+ putchar( '=' );
+ nCharacter++;
+ }
+
+ for( ; iBlock<100/2; ++iBlock )
+ {
+ putchar( ' ' );
+ nCharacter++;
+ }
+
+// Int_t sec = TimeInSec%60;
+// Int_t min = TimeInSec/60;
+// Int_t hours = TimeInSec/3600;
+
+ nCharacter += printf( "] %4.0f%%", ProgressInPercent);
+ fflush( stdout );
+
+ if((Int_t)ProgressInPercent==100)
+ {
+ printf("\r");
+ for(Int_t i=0; i<nCharacter; i++)
+ {
+ putchar( ' ' );
+ }
+ printf("\r ... done! \n");
+ }
+
+ return nCharacter;
+}
+
+//####################################################################
+#endif
\ No newline at end of file
--- /dev/null
+//-----------------
+#ifndef CONFIG_H
+#define CONFIG_H
+
+#define SAVE 4
+#define ENDING "hld"
+#define CONTENTIN 1048576 // --> 1 MByte
+#define BUFSIZE 32000
+#define MAXSENSORS 64
+#define MAXSTATES 1140
+#define MAXCLUSTERS 1140
+#endif
+//-----------------
+
+// #include <header.h>
+#include <dirent.h>
+#include <TFile.h>
+#include <TTree.h>
+#include <vector>
+#include "help.h"
+
+using namespace std;
+
+class unpacker {
+
+private:
+// File info
+ Int_t fNrFiles; // Number of files found (Unsigned integer 4 bytes)
+ ULong_t* fFileLength; // Length of each respective file (Unsigned long integer 8 bytes)
+ ULong_t fLength; // Length integrated
+ ULong_t* fFileEvents; // Number of events in each respective file
+ ULong_t fEvents; // Events integrated
+ UInt_t fNrSensors; // Number of sensors found according to the 10 first events
+
+// File access
+ ifstream* fInn; // Link to each file
+ char* fRawData; // in stream buffer
+
+// Coordinates in file(s)
+ Int_t fCurrFile; // Current file accessed
+ ULong_t fCurrPosInFile; // Current position in accessed file
+
+// Output file:
+ TFile* fOutputFile;
+ TString fOutFileName;
+
+// TTree
+ TTree* fFrameTree[MAXSENSORS];
+ bool fSensorpresent[MAXSENSORS];
+
+// Methods:
+ bool checkFiles ( TString Dir, TString File );
+ char* getEvent ( UInt_t EventNr );
+ UInt_t getNextPseudoEvent ();
+ UInt_t getNextEvent (); // returns eventSize and provides event-data
+ UInt_t getFrameInfo ( UInt_t EventSize, UInt_t FrameNrCount, frameInfo*& FrameInfo ); // returns number of Frames along with its info and position/length
+ void getStates ( UInt_t NrSensors, frameInfo*& FrameInfo, TTree** FrameTree ); // browses event data for frames and the corresponding states
+ void getClusters ( UInt_t NrSensors, frameInfo*& FrameInfo, TTree** FrameTree );
+
+ UInt_t littleEndian32 (char* RawData, Int_t Pos);
+ UInt_t bigEndian32 (char* RawData, Int_t Pos);
+ UInt_t bigEndian32Switch16 (char* RawData, Int_t Pos);
+ UInt_t bigEndian16 (char* RawData, Int_t Pos);
+ UInt_t bigEndian8 (char* RawData, Int_t Pos);
+ UInt_t datalength32 (char* RawData, Int_t Pos);
+ UInt_t numberStatesInLine16 (char* RawData, Int_t Pos);
+ UInt_t addressLine16 (char* RawData, Int_t Pos);
+ UInt_t overflowInLine16 (char* RawData, Int_t Pos);
+ UInt_t addressColumn16 (char* RawData, Int_t Pos);
+ UInt_t numberFiredPixelsInColumn16 (char* RawData, Int_t Pos);
+ Int_t print_progress (Float_t ProgressInPercent);
+
+public:
+// Methods
+ unpacker(TString Dir, TString FileName); // ready
+ unpacker(); // ready
+ ~unpacker(); // ready
+ void scanFrames(UInt_t Start, UInt_t Entries); // ready
+};
\ No newline at end of file
int main(int argc, char** argv) {
-// gROOT->Reset();
+ gROOT->Reset();
const char* InFile;
const char* PxFile;
--- /dev/null
+/****************************************************************************
+** Meta object code from reading C++ file 'scurves.h'
+**
+** Created by: The Qt Meta Object Compiler version 63 (Qt 4.8.5)
+**
+** WARNING! All changes made in this file will be lost!
+*****************************************************************************/
+
+#include "scurves.h"
+#if !defined(Q_MOC_OUTPUT_REVISION)
+#error "The header file 'scurves.h' doesn't include <QObject>."
+#elif Q_MOC_OUTPUT_REVISION != 63
+#error "This file was generated using the moc from 4.8.5. It"
+#error "cannot be used with the include files from this version of Qt."
+#error "(The moc has changed too much.)"
+#endif
+
+QT_BEGIN_MOC_NAMESPACE
+static const uint qt_meta_data_scurves[] = {
+
+ // content:
+ 6, // revision
+ 0, // classname
+ 0, 0, // classinfo
+ 13, 14, // methods
+ 0, 0, // properties
+ 0, 0, // enums/sets
+ 0, 0, // constructors
+ 0, // flags
+ 0, // signalCount
+
+ // slots: signature, parameters, type, tag, flags
+ 9, 8, 8, 8, 0x0a,
+ 17, 8, 8, 8, 0x0a,
+ 26, 8, 8, 8, 0x0a,
+ 40, 8, 8, 8, 0x0a,
+ 47, 8, 8, 8, 0x0a,
+ 56, 8, 8, 8, 0x0a,
+ 66, 8, 8, 8, 0x0a,
+ 74, 8, 8, 8, 0x0a,
+ 85, 8, 8, 8, 0x0a,
+ 97, 8, 8, 8, 0x0a,
+ 110, 8, 8, 8, 0x0a,
+ 122, 8, 8, 8, 0x0a,
+ 134, 8, 8, 8, 0x0a,
+
+ 0 // eod
+};
+
+static const char qt_meta_stringdata_scurves[] = {
+ "scurves\0\0Error()\0Unpack()\0OpenBrowser()\0"
+ "Save()\0Newest()\0RandPix()\0Close()\0"
+ "clickDir()\0clickFile()\0calcPix(int)\0"
+ "calcRC(int)\0RunAna(int)\0NextFile(int)\0"
+};
+
+void scurves::qt_static_metacall(QObject *_o, QMetaObject::Call _c, int _id, void **_a)
+{
+ if (_c == QMetaObject::InvokeMetaMethod) {
+ Q_ASSERT(staticMetaObject.cast(_o));
+ scurves *_t = static_cast<scurves *>(_o);
+ switch (_id) {
+ case 0: _t->Error(); break;
+ case 1: _t->Unpack(); break;
+ case 2: _t->OpenBrowser(); break;
+ case 3: _t->Save(); break;
+ case 4: _t->Newest(); break;
+ case 5: _t->RandPix(); break;
+ case 6: _t->Close(); break;
+ case 7: _t->clickDir(); break;
+ case 8: _t->clickFile(); break;
+ case 9: _t->calcPix((*reinterpret_cast< int(*)>(_a[1]))); break;
+ case 10: _t->calcRC((*reinterpret_cast< int(*)>(_a[1]))); break;
+ case 11: _t->RunAna((*reinterpret_cast< int(*)>(_a[1]))); break;
+ case 12: _t->NextFile((*reinterpret_cast< int(*)>(_a[1]))); break;
+ default: ;
+ }
+ }
+}
+
+const QMetaObjectExtraData scurves::staticMetaObjectExtraData = {
+ 0, qt_static_metacall
+};
+
+const QMetaObject scurves::staticMetaObject = {
+ { &QMainWindow::staticMetaObject, qt_meta_stringdata_scurves,
+ qt_meta_data_scurves, &staticMetaObjectExtraData }
+};
+
+#ifdef Q_NO_DATA_RELOCATION
+const QMetaObject &scurves::getStaticMetaObject() { return staticMetaObject; }
+#endif //Q_NO_DATA_RELOCATION
+
+const QMetaObject *scurves::metaObject() const
+{
+ return QObject::d_ptr->metaObject ? QObject::d_ptr->metaObject : &staticMetaObject;
+}
+
+void *scurves::qt_metacast(const char *_clname)
+{
+ if (!_clname) return 0;
+ if (!strcmp(_clname, qt_meta_stringdata_scurves))
+ return static_cast<void*>(const_cast< scurves*>(this));
+ if (!strcmp(_clname, "Ui::MainWindow"))
+ return static_cast< Ui::MainWindow*>(const_cast< scurves*>(this));
+ return QMainWindow::qt_metacast(_clname);
+}
+
+int scurves::qt_metacall(QMetaObject::Call _c, int _id, void **_a)
+{
+ _id = QMainWindow::qt_metacall(_c, _id, _a);
+ if (_id < 0)
+ return _id;
+ if (_c == QMetaObject::InvokeMetaMethod) {
+ if (_id < 13)
+ qt_static_metacall(this, _c, _id, _a);
+ _id -= 13;
+ }
+ return _id;
+}
+QT_END_MOC_NAMESPACE
--- /dev/null
+10000
+ 136386 33775 105102 2670 145304 114832 76719 120045 7319 74235 79021 22625 109765 132041 124878 137778 65704 74491 15394 27050 40117 106782 64306 55120 74263 57032 45689 25300 120210 3289 161686 153985 108392 100869 69913 159736 114100 126471 22447 39604 126192 132212 108158 85182 40614 7974 159673 56008 137424 158519 96126 78318 56938 151246 89093 50482 133447 114393 4804 136737 46703 95301 15753 84084 101727 9601 32296 62310 113168 54743 101914 73472 123467 146584 158654 100593 154558 152439 156602 62606 81582 23352 77436 75032 8710 125514 78669 51546 66830 151918 98250 162132 1783 16446 97971 11384 151143 96793 124552 142398 32819 134536 36490 115916 127302 137083 41098 50365 127797 103705 68460 151149 15253 143492 96371 118295 39631 111553 3953 106461 97583 38715 39217 99367 157562 73700 47263 142817 4605 5928 119327 139825 76976 92329 26365 140791 63525 67463 25268 127834 5280 93728 49608 122934 7845 82491 75341 149876 28156 79294 90449 62252 118010 66179 161619 109684 76391 42994 23125 17509 151322 78964 93846 5406 120211 139714 68931 21786 30877 27067 29335 16997 150001 139580 99489 59454 123568 127645 75260 150529 126409 129782 153220 122140 73578 63724 101647 33215 17745 23593 112180 111591 22516 117586 65914 98742 123029 87700 34349 90418 51279 63684 107416 137792 139776 41017 133758 33968 105174 43131 18609 65696 7025 108342 124348 17116 108578 162507 50331 126323 20213 99023 72026 145129 50721 137940 77984 110262 162152 48845 34793 149944 49041 142209 58360 22929 17338 26231 56897 59024 5874 12018 124720 115299 120360 83181 132415 63050 79800 119259 23485 36525 154794 32023 118167 142028 106475 132663 86402 39252 15620 121195 23308 64661 34028 81668 87590 51366 44411 80999 46903 50285 93017 5735 149890 25428 68624 149452 105229 124395 7050 78266 113302 39073 30545 25954 82061 99720 112356 121313 115340 4176 144621 116513 38204 162801 140615 26083 41325 55726 72986 91610 85256 78721 27819 5770 104150 96444 155222 145891 157351 98784 160669 41277 74370 127727 67231 156431 61559 13700 111856 113412 120276 90589 64037 158480 24014 38765 18675 65339 31003 28173 93462 52771 43407 121281 58541 84069 154237 47876 64072 82213 146660 161253 123490 157542 123092 24834 148085 18764 140934 94053 132176 95322 121154 132725 24426 145169 108002 145502 147020 75518 7787 74594 128289 51194 132388 123343 135263 57249 107731 135847 139462 88503 131213 97065 16670 88417 58411 43693 135857 92921 112381 1803 150587 79219 128629 66380 123733 108243 47513 33363 116031 58619 161653 103737 25119 119108 9625 82369 60951 145472 158343 85966 47309 26032 102636 72239 84443 38016 52444 156812 130937 95127 52148 14669 57869 118529 138402 108278 52767 58773 104043 83893 159663 102642 102774 157126 18738 95229 13716 2560 121262 116353 11311 142217 154369 63755 69654 55930 62693 108078 13874 159162 60719 88789 95899 163273 31179 148666 56159 135222 12307 76564 65509 51461 115850 56747 6712 45191 70463 9272 102965 123328 20583 15807 48321 84338 85461 40763 83544 20866 148842 97418 14141 146073 122719 110040 143459 153898 29330 136130 59744 144038 149206 125253 29611 99168 16112 36323 80871 23088 147995 17949 82928 2690 33756 131250 23541 55729 6125 107085 76595 91479 141015 27248 71665 97847 73800 151636 85857 103131 58390 145602 17793 41708 41479 47404 77388 159992 20240 158259 17192 2347 10320 100120 5038 146476 1994 130979 36317 8120 72176 49425 36111 149703 76673 44288 18174 150474 30036 40544 24229 24938 122658 144422 3158 164137 25938 80546 158241 148578 72918 111945 150926 19750 92476 102739 150572 57567 75568 95204 66255 124993 131316 152470 138179 9716 4757 59277 142153 45301 20018 103603 2071 164440 106762 126890 21420 95074 109581 30850 41131 31131 50601 23821 123607 89852 110906 117686 40222 120453 61038 108050 107035 33329 117767 111792 92606 30544 157093 49136 134147 159164 47688 75021 95997 111090 32954 25183 54783 63804 2826 22426 50917 129048 146033 140769 74066 34343 140302 50800 154796 137852 158851 95944 5293 47242 41848 34411 14298 135454 83547 148445 65242 67747 159979 161239 12949 129445 4245 129761 125761 26671 14791 88921 109217 92072 99499 143560 2998 150299 132469 140850 79774 82655 63528 143285 117066 77826 112851 137125 162784 12206 38984 93387 109957 154334 159344 67004 158579 123217 26877 19362 138008 115798 128579 151809 106252 3703 136220 9345 81066 152507 110782 163721 50147 88179 114900 64486 35143 22649 163782 149749 61634 91281 93818 50080 84737 160822 42771 42066 21811 116587 74121 127225 53804 60042 4101 159907 132775 115846 75085 119394 60740 9431 106053 148919 124331 4651 120574 83492 2545 104435 145126 30338 32366 131718 51587 129700 111001 30166 151512 109647 83265 162145 7496 137069 158700 11597 131088 125587 127443 142686 15605 22295 152117 121658 107726 47072 62822 62413 130564 1879 46315 32218 135726 12145 20317 99539 123147 50483 21675 3418 133748 120332 10914 104929 113144 124911 6642 9355 86466 149328 127360 45273 135557 19643 152999 119141 82465 151924 83817 20856 152885 130132 155474 59235 78790 9904 95286 36049 60387 116961 141867 28248 7918 152781 69689 121062 111804 76331 66930 134782 59771 28402 116567 131840 48045 103678 85093 67022 89715 105423 87879 13224 6179 13977 8971 84969 23881 104258 57530 84269 157731 33509 49029 122802 118718 57336 5230 131562 124266 140012 127845 89180 90691 93798 73738 130882 115403 140760 157109 157338 106845 163518 13241 52328 82599 139522 156586 140130 160303 84942 7751 43444 84703 67066 98675 78551 72296 64349 139329 148821 26306 62622 10136 56616 136360 77530 108532 47744 68751 99982 149408 112108 34124 99161 164437 116724 9307 91647 27478 3723 10701 137629 47167 31917 38807 145842 46980 111104 146703 20422 30549 109522 19556 40685 92428 118216 45294 140172 123479 81789 60204 6212 115913 95877 4761 3261 105185 32920 30739 108908 146022 2481 92587 12051 143688 9054 161431 88904 155757 15965 55965 99391 35521 96651 36154 127949 151379 81448 38746 45482 163237 35462 51694 113263 131340 56455 116524 70637 89376 83776 116057 69510 86257 42756 18073 64057 51810 116016 89474 41680 131982 145439 77583 104015 12714 113737 2589 100605 29298 143735 146088 129047 13309 31894 12934 81161 24862 65971 88310 114238 149747 38479 120260 70116 17748 74845 70685 69558 24973 96671 47750 93467 76223 61846 133995 25449 9695 73096 62567 141393 50943 42767 104553 74661 53999 81926 99523 119970 4348 150273 103829 145228 41157 110457 162976 52514 15255 66646 14000 111926 50909 107467 124661 112755 12086 86623 122450 149190 34468 136125 26069 75533 136890 37242 129532 155328 73278 83615 159676 57663 21556 139016 35333 132014 72616 87847 83781 139263 101847 29819 126684 145827 90993 73551 157913 11728 132513 13720 160918 149845 21099 13138 57359 160741 142671 149199 68131 162798 142988 125795 120866 52628 97640 23504 125245 19599 107285 35132 57959 73617 161816 37898 5991 132323 10450 138504 146043 107880 76110 66513 65491 89248 60384 60344 2543 43696 128476 101853 123196 24895 159232 112336 122535 16848 8205 78646 60646 43337 136605 134263 141665 111015 69497 147656 13963 16459 56785 96518 124339 132895 99543 23942 158655 159928 84286 97711 140136 149274 33676 33956 8281 129420 146292 67328 146269 154498 145975 41027 134347 53204 111802 110125 100732 117811 28405 114695 134270 85190 147725 92721 154597 81381 53175 147365 11933 73973 15700 88581 57360 151776 122537 2153 115309 102941 69482 95690 28063 151969 73229 98923 141685 121543 43160 76529 73466 8077 127736 144248 29780 109574 7593 18489 25067 60768 102366 139400 134741 118066 62093 128613 40467 18742 67279 155776 58195 136761 22090 86259 59354 31831 19294 35151 153374 111681 163352 109443 73529 141712 139223 17215 149305 157713 144682 146585 94191 118194 51950 148770 14399 14676 23349 135541 81955 115637 27849 155228 74239 114108 151094 106070 69914 20357 30068 5392 68550 27532 114835 142080 105756 88171 95807 25685 16508 11114 6382 110699 65820 160732 30093 16732 9520 155842 152273 27987 105591 14234 119727 13942 64854 104933 56524 134768 61803 86592 76672 130353 50636 25620 43057 92904 50303 75377 118589 66811 86491 61483 114022 88823 56328 144116 105555 2360 70582 91941 30348 112686 106175 86587 126628 5142 25633 119665 76422 87436 40369 153095 154301 27518 115227 31471 120422 106848 9636 2965 129851 71119 116987 52786 127447 31727 158342 66320 102310 84395 33180 49108 24682 119767 112248 132224 81912 66025 145159 105860 42907 68878 94274 6937 18217 62257 127359 120259 3217 136995 59736 69580 144627 10835 122366 42698 144963 114820 109018 81385 135727 142198 67005 96922 32590 13365 114502 15903 144929 156875 161210 47919 85261 2259 2648 147518 66130 59419 87247 139638 119155 156827 118377 66503 49817 97587 45578 101150 40718 126963 70989 119428 130480 104423 152018 80357 104194 37145 96260 83235 28132 28094 67667 49905 132753 70315 133935 32996 66247 157694 6746 121914 148633 61635 22529 134962 159222 68107 70224 136452 131582 77726 89993 96174 16261 12635 113044 120455 49780 43416 140203 14424 8023 41982 140776 48809 134776 110284 115056 126582 53542 7595 45839 115177 30124 117314 45024 34744 21650 117988 99376 42093 33125 115638 157129 146169 70205 143421 126097 44520 157846 70632 23014 158687 45521 71824 64088 92317 123392 127182 145860 130987 109534 31661 161112 60960 76685 132368 82610 131186 132806 118499 7391 102443 68249 101032 19236 74966 78566 81846 55999 7036 152478 79013 134511 87349 60941 147254 127618 143313 112353 7776 9086 44089 68736 22284 10569 87858 153470 79888 40469 97373 118843 45230 32518 138080 56709 111084 54038 112708 118120 40628 128233 54467 111652 152095 157303 6705 133461 119033 150018 79926 126809 95616 60528 132057 117900 7609 54027 105482 87497 31009 139368 40453 76239 5998 115045 132948 53594 105595 16280 108226 82735 144514 162693 28499 130721 154108 35204 34806 43765 121734 114732 107086 51463 111772 9767 105875 119382 63795 45470 40991 94804 18950 17956 107555 127348 69513 74616 15054 9220 90896 123280 91956 6034 56597 73267 147218 92172 108073 127495 150418 159318 68694 35993 105202 78461 78381 58696 78768 60363 36200 110084 79313 54156 51752 143173 60182 62880 158227 69402 90288 52131 161358 96323 45240 115926 106102 26570 144610 48288 154066 129140 41718 159272 146920 8357 77627 142129 23638 137990 12441 133722 153815 3109 121986 67612 63291 121378 162351 132694 45779 48594 64676 78614 30346 14714 18828 56917 159324 67116 45095 59089 45346 140879 58335 128779 149236 72474 105020 6986 146976 53973 77221 71415 57082 135719 139027 56886 91210 72002 126092 136989 57108 24880 152115 87454 141995 5055 144371 71943 8684 125978 131032 156430 100969 125879 119321 20830 134865 58453 27816 115953 112426 41549 123880 106021 11381 33531 162907 102591 105533 123111 10204 99153 84503 162319 20720 60610 103886 132554 49082 61694 97698 39625 162663 160090 158946 120005 129067 153912 147822 15645 36962 125883 139525 142983 137264 109569 140002 73967 49214 33737 84171 148368 118241 17114 105600 115363 57513 104803 82029 43107 103009 116240 82732 36297 110442 12303 156302 10133 102727 74748 25778 139689 116785 6120 108985 27411 80088 94711 163549 100771 13703 115902 54398 119303 1889 48423 160619 83919 91530 97740 136671 8375 134037 17737 123078 60964 27870 162317 135712 53649 136118 4568 155465 23527 10689 98562 50939 27289 129785 151000 64572 143489 101014 55482 33416 102903 103905 28147 123334 29548 62400 30629 140323 48366 34025 28025 12749 30454 163738 66398 2418 55975 126612 115507 91049 11663 79308 54946 80393 134947 34301 135875 2476 73716 73893 133023 133563 39953 29535 164192 14388 162485 149071 48413 24622 161820 15379 22472 16063 127291 157217 142675 13422 18890 90850 29243 73836 24137 109636 145296 160838 16135 84284 5179 90028 51419 138742 66493 80955 73558 80881 77552 56741 65806 38686 155073 81185 163559 90427 33761 124962 81756 112948 74896 100646 140311 104139 110995 100960 150287 26915 95911 111199 101090 90563 162618 10456 93569 14197 84014 8562 28261 77268 10881 66948 28578 64619 93393 62339 126093 111661 111800 137501 148820 22735 75753 30439 123695 60152 57354 156118 60687 2665 91320 87763 101795 38288 15444 52505 58815 126406 80766 136083 137287 84226 75560 102378 85357 105465 147962 51239 47141 119576 34171 69876 29441 64610 27684 26105 58476 17914 23305 163541 45747 111068 99448 20547 63024 151953 79362 125942 3344 151957 33854 24082 61630 136232 109440 71573 91514 154846 118715 45202 125529 125103 11155 126651 89299 37261 19239 43726 60566 16892 89473 108146 52853 110020 107682 141318 125895 67736 81174 111964 101590 105257 110106 71934 111314 80020 76835 36784 32847 58550 162314 157950 69705 59589 81362 43478 78829 125088 40556 32233 151073 85214 21598 31717 27008 99429 157612 31257 14715 40201 69359 56484 150307 77806 41805 32245 157826 118641 69030 24785 11303 1968 119247 17520 61557 137121 60999 140386 32833 101555 109132 120418 20882 67242 152136 150290 80372 15659 117899 57085 85019 8495 41505 162825 50301 73750 154763 105454 79292 116060 116757 81260 5931 134277 142818 143053 131788 53828 112398 3968 162960 66929 127250 155577 1610 6573 63824 119509 63659 85355 64517 105164 18804 114818 115426 110079 54384 28831 162651 107653 110091 2694 12554 23533 82259 144343 77362 28770 84823 10946 32211 46185 11773 21900 94349 13384 130873 158173 69405 28644 14152 133922 70320 32956 19364 122259 143035 10260 151090 139798 117913 31805 79005 130468 55339 161264 45435 132701 24146 130258 143647 158757 10555 91933 117169 41416 41829 82155 33702 47746 47311 47854 15781 117632 17323 137545 74003 160358 147806 161605 70781 99831 27522 149786 19373 145162 46358 152074 105821 113128 66346 98690 60195 94791 49972 101612 136620 68639 71826 18478 115950 119680 136659 4206 73515 108317 78209 67986 26747 10438 75279 63090 140361 59177 64014 159734 140851 46884 82433 17296 96525 85291 52499 156720 14194 102471 92444 87326 107622 164270 105804 57684 54575 13088 61891 128090 57917 76612 132588 84664 23563 144379 147754 163924 140068 148280 94282 115032 131677 10827 68840 62314 32630 121339 53146 149224 160322 82103 70662 102056 112979 159741 8084 126067 158144 72687 120496 5380 141787 39272 28943 120279 123538 26979 94459 42443 146003 8232 68601 48956 7058 37744 106807 21080 101242 78819 91743 37410 95817 141234 133663 103901 37925 62431 10700 158421 67812 152488 134205 96755 43391 28367 123735 74362 7322 15621 54478 15554 20734 39946 22612 58478 146753 19329 79559 18619 98148 107814 158430 28077 83160 62717 68491 57597 125149 15703 50130 27073 2303 120847 60340 148094 149214 120587 56569 156537 72720 111047 108603 93455 87505 67728 151933 4882 87057 2116 23502 121717 109930 16044 86307 129602 78761 154798 21311 140422 4613 7953 1607 109317 65312 164348 28035 48639 119047 84604 141688 25880 132163 20915 119335 53780 88643 41892 161063 112212 44009 18677 4554 90451 90861 156566 49994 16283 114389 24529 123296 58855 128536 3237 124167 63508 31273 109318 16668 115877 21630 144948 18665 42546 34907 8957 67701 76799 4132 116426 57320 125209 120980 96442 148353 74962 82949 125863 43990 122044 18830 109038 125282 142998 6659 93067 22940 23327 145456 44571 104787 164121 23629 139694 7191 91330 153005 113723 144268 146838 73045 99360 65234 3599 81825 76708 86548 17085 36683 67050 139130 157914 10201 35036 135024 16860 128103 94476 142587 107671 139047 81486 42417 162676 157692 152008 24631 144809 99843 3011 125759 7000 38884 127505 10600 57221 38325 33660 74307 11521 37223 149949 3547 47424 19097 75083 64284 83712 3671 143383 127895 79231 161381 106824 76019 153185 92944 37162 68618 26900 40174 130890 33900 15570 29019 146900 72791 3857 117073 83610 15378 154296 67671 121325 35832 23280 30520 36628 43504 136591 116523 5512 49934 112016 48848 62466 35825 141793 99628 40955 2805 76314 108357 139105 91884 137377 120118 7815 82398 156612 162111 86582 48561 134455 46374 79081 107595 89879 152184 58230 31903 36231 106758 80751 98697 79095 56656 32437 120050 59461 108752 62520 135079 137148 34009 89309 72448 111755 33636 91359 38991 132259 12053 87552 37338 58427 144933 84818 89441 139675 116721 125672 17057 133985 160881 96152 24753 129831 152714 20727 9207 49346 155806 146355 19867 15739 52916 131622 49375 144275 4725 118146 156328 92277 155484 151267 29534 134529 70198 118976 44828 123431 15272 164285 91528 10266 94549 52794 140097 17887 73521 149304 67234 66283 23613 118090 119199 155236 103977 97586 96473 56235 24538 125263 45831 112318 91309 116872 119028 146797 161700 13083 162070 96609 104612 108848 127670 93918 19569 145557 103951 105385 149303 103902 5780 7029 56104 124980 98777 160081 159078 29362 50428 120129 91137 32771 66559 16559 149643 122099 99868 81967 135182 32562 12688 10418 141410 140358 104336 160979 120027 42399 100476 103443 82813 106257 46984 75429 1861 145761 6134 160939 9235 158962 51692 36885 25845 54763 155844 9600 10974 89824 91567 146157 122387 40767 93087 97909 15237 31536 29513 71777 10447 66501 111732 93261 6870 158716 2802 111131 138589 8937 42695 84336 2011 94387 121221 130257 149151 47689 139857 96637 137514 2049 13418 94013 42816 106506 128434 160454 138042 157947 2855 85001 160961 114587 12374 1943 43927 117577 49587 16628 126514 92282 37476 65037 20781 95210 29406 106444 142899 105776 37194 114525 44337 50612 145050 87153 157118 44109 18231 65784 36168 123486 150786 31241 8697 99672 135585 52624 51361 121684 5764 114387 48078 43241 13537 5371 74963 145343 111816 51974 85231 85522 103012 129568 136134 18686 153234 63877 62795 107977 66173 35476 65576 153471 3229 10785 87256 138814 63410 138617 94610 69174 87117 79200 48927 37166 84572 123890 119021 132900 112377 140765 52534 151901 40957 125180 4699 28303 125569 67495 72793 25855 102971 74881 13438 106200 85666 37206 15639 149076 46761 135965 125962 140202 7243 147046 98205 62776 114058 147094 37653 103104 133107 78611 137806 43426 124478 141813 52731 86845 78896 127612 36795 121609 149791 74002 137248 69491 20450 18121 58366 156415 80595 32681 100170 164153 67398 162946 112323 151004 137112 151939 118223 49835 150848 26653 45950 104979 82505 69307 17987 44229 139596 116616 47468 97624 137066 2102 92502 64105 19209 61695 100787 17475 65605 34358 66310 5582 52362 105456 55417 139722 68622 122184 7713 1807 25603 128100 148437 131705 38321 113081 82434 85789 147217 156012 24403 73832 156630 43613 72039 91529 61088 137645 125887 127398 131469 116272 123398 26618 6093 145101 148802 13806 83420 110917 78418 65969 76735 116739 13163 95681 139041 160380 22317 163444 4836 13059 143569 13388 41101 38769 87545 102680 46535 69082 155464 149894 26592 16195 155987 108205 106305 25738 110027 18836 91707 123274 135575 41382 155467 138275 11896 145197 79623 127356 122878 93011 2569 98160 14668 106069 34952 163604 9263 126928 147610 138256 143123 137710 80573 78745 78127 42823 22884 33475 134531 82670 3163 112425 8761 150803 21324 123057 130112 37460 84525 87103 130471 23606 19375 81652 129676 156727 15880 138939 163490 47819 31515 64905 110260 86464 107728 69656 56451 12883 152326 59614 125309 161087 146930 83145 118256 111154 120605 139294 32369 21701 162900 154144 39865 63200 81495 55745 138652 135775 155747 164084 85488 111662 21172 27329 117830 77623 40212 104268 102033 35979 54792 19291 90747 76408 134828 98109 100078 123084 74486 163278 141092 66743 136042 110979 56603 93138 157199 12738 160555 12483 142467 112497 16792 153277 100368 55337 125768 155160 11140 153027 91731 87549 90205 60671 122170 24395 120267 30769 21785 95471 97512 94340 142962 90627 21590 25967 81938 136728 20634 30933 113308 69643 57552 66612 159432 94432 121949 119312 20217 133090 42964 111948 54751 133169 109131 113433 157564 144202 115861 32006 12339 44313 9080 39478 2415 35048 57929 139144 158082 88862 86564 164238 82926 153176 157782 11471 109237 47719 134088 12951 90683 80148 67702 160364 125791 117648 152040 62325 95962 38525 94331 44813 19351 39924 20804 124166 11484 78733 97422 106078 104107 120498 104428 21146 107786 135017 153536 144042 39729 102999 5349 119877 107214 16292 58974 88377 15129 91448 63414 109461 136262 82765 85897 93578 143444 33893 6423 74978 139971 110530 29589 78512 68188 73887 111347 139717 163935 26013 13558 101047 133762 69947 42373 133587 22751 37859 158476 37881 129307 56002 83854 36193 75280 106263 129771 52836 140156 136194 64326 114239 17349 30427 129263 22049 104315 11234 98279 38874 139647 48349 76433 107521 54809 118806 11732 14072 156665 4320 51953 162723 72319 92790 8627 12694 56674 163863 152850 129380 37714 83241 29241 103489 41803 70068 51236 76594 108942 24995 124943 19488 132517 116264 74806 80761 130337 2096 21594 118802 58693 18429 25234 151483 129456 37928 144669 63943 44674 62756 101517 127915 91997 141518 106230 98577 26866 16936 41632 51862 78392 163520 120891 28768 72438 35764 95617 74534 57358 48532 69739 12299 73766 157735 78267 48206 136516 142210 112009 15302 141479 47638 79730 4100 23269 20072 102678 50135 37009 144310 38509 115401 141942 95912 80681 150892 68189 10411 161939 125547 161343 65790 74359 69221 160037 89138 53939 67178 65461 82480 143452 150099 98722 147552 7480 118795 20854 160016 92316 101676 135037 144229 14242 65062 161422 133251 108345 161186 29422 103800 163488 40293 109533 94150 129432 163473 161328 131405 100046 14432 45481 84257 113155 27145 28250 2574 150400 124778 31402 86188 9743 100431 32013 5277 99678 101776 113623 31488 131199 153935 131488 108004 97581 59750 71548 95166 55190 139465 131724 69623 19058 152493 119290 148604 14855 121864 133116 76145 153266 53416 106585 99521 90359 138598 41310 126549 74487 91445 158037 79493 60150 20826 11186 119900 28887 42864 111603 2464 111100 117738 123923 97705 71140 106639 49073 129516 10379 125218 53406 2427 89439 90667 77538 130749 51328 152025 158707 145878 64847 72312 40140 85673 20010 96552 114560 42267 53537 110486 96517 11572 42303 1769 118211 91376 67797 65102 153107 121203 65409 155534 147154 156076 67184 143917 155721 41347 123907 157080 113659 100559 76866 70181 31223 127938 69567 10003 15587 14165 106520 27159 158868 44802 81882 84357 49111 146984 8088 4427 46506 100134 88093 139094 103831 34221 53635 93664 11654 21369 125739 48725 98235 30032 16461 162685 128864 12385 50276 6008 39544 43257 153210 57939 64126 36434 39035 72214 143261 22053 6460 65466 161148 46803 151295 140468 111768 162950 98349 8132 148187 30696 140564 27493 10788 64136 142278 61064 70145 15935 40833 57467 10386 41471 93901 49421 113685 7786 71475 56658 73253 3247 103461 109453 154542 14553 157733 88116 112902 6928 80110 77054 108088 44116 87842 6337 20506 85418 12994 138841 62764 70461 149227 104235 100875 135161 52033 108661 143148 45203 118426 82907 85176 61991 71561 36242 53837 96190 85656 156214 103118 102279 67380 45318 91734 154055 39925 11265 12879 74029 8135 98618 12376 45522 4403 90695 147551 148524 43234 64570 4325 41737 72644 142967 95574 2946 62735 22413 42576 101526 26305 24406 18545 118040 12574 58471 65817 116135 7862 76358 60782 106480 25246 106304 110884 128568 133512 92547 47716 10858 93630 52041 52595 29120 84682 105732 28368 43607 148308 129894 69912 6826 148440 124464 121800 41023 24393 8560 48885 37263 5854 91877 62510 112159 36873 127590 16295 65933 9418 27153 96075 61460 16260 96461 90580 37454 36305 118948 81061 121125 82955 87486 127951 2019 46062 20376 43042 6968 131336 28439 44231 137190 120316 43253 19973 93702 4955 36268 96147 14374 26334 75834 118594 59307 102926 92560 95612 158387 110134 50849 11966 31732 115312 13985 77794 72200 159427 21274 140048 21978 65506 111351 78806 108759 67836 6620 113715 40617 39279 128089 40550 65613 140435 95656 61432 77473 22329 157044 6484 68975 144405 18450 100707 30342 32435 115013 102542 128374 136288 76703 150352 138306 124566 63271 81177 128914 6403 131404 3643 45683 30117 146594 47808 107064 12874 45753 121050 137603 36909 127534 40690 117827 145985 141397 148169 114932 27035 77419 163323 98038 164283 135741 159116 61666 153430 58655 4582 55459 62298 152665 22088 145404 136985 129153 158279 16850 84315 129994 156160 45961 4797 108099 128458 82706 90380 77503 109741 111715 91434 43688 146266 89829 115941 139494 88008 39996 32261 92590 95455 31072 15879 54055 112988 152864 17320 41891 3827 101635 159987 84109 113195 102198 149079 30013 129090 163094 76267 74917 88640 119955 55295 12582 70009 131302 100590 110005 100075 129692 141972 67659 145571 30139 14760 132547 149860 56651 72886 22119 3497 106228 42207 25932 78882 5409 23138 91661 16839 111779 45728 8646 60873 52249 139948 161463 98766 10648 61779 11362 78307 41462 41502 110521 25474 86231 119920 47593 85392 123417 90334 134261 116266 149655 107546 139404 75428 124385 21807 57668 69544 82680 46430 43604 14767 145196 54252 76546 156559 69072 54520 134573 162139 101554 96559 18994 55586 144152 104387 115515 5110 115252 57888 158927 159311 133805 68467 54320 155612 62647 123864 72405 109077 1581 23684 24898 158233 100231 15569 61417 91263 86654 160069 26929 17325 13175 19027 97989 54074 134543 103100 63346 20419 97500 56385 13842 65417 61364 4675 124012 128539 3713 66632 24255 28611 161378 60998 146580 56907 152262 67346 153488 115703 84671 103176 134731 16773 157250 103386 54709 60317 153885 47606 10672 155815 108971 15347 83981 3607 143886 84552 7320 147031 45319 35932 142521 106318 16624 135940 29204 20483 123541 144907 105154 60829 113750 121927 154591 153648 12424 43412 150478 102822 91019 161150 29261 136502 10609 113242 140109 91007 134306 83941 72150 13737 48783 120055 73010 149259 93493 78889 64791 135159 76230 12653 91199 64933 102814 40135 44858 87404 79469 72389 19178 108730 145403 29787 158484 119624 120794 126902 37677 27057 140640 94063 12352 31319 103585 31188 14691 133590 46589 79482 39373 122819 28647 67084 124265 67973 107220 105635 91889 123201 114536 111067 2556 94051 140854 161040 47787 32273 122055 21976 161730 33319 52551 8194 64638 92648 141783 15841 162750 22484 31835 36236 81816 60483 103320 142593 64968 147052 18852 156858 40878 133388 38549 43434 61551 38586 45850 148189 97153 67826 144031 130472 56890 88737 29223 86050 64632 147464 82913 13412 55661 41957 10407 95493 75375 13170 139914 2857 54048 107414 143807 97482 105477 93835 72580 151327 12648 106246 53265 93191 70830 46667 16040 36565 132718 17185 18142 152143 143214 133954 144316 121683 144361 73921 76857 53848 87091 153283 159106 77651 31321 137025 111645 136798 1484 18338 122237 14132 124584 112014 107323 131926 158682 123363 105004 125512 77060 59658 48279 54386 27724 26707 112581 6197 100628 125950 162445 124232 113345 92175 138395 81178 84153 52088 39003 8437 15440 100099 56964 59275 66137 49758 119151 111782 30323 160061 21222 131547 20880 133803 137744 58020 93866 70814 118764 143723 162989 91272 59014 99326 111937 111102 100634 150940 56052 52587 85151 113016 111862 87800 99286 65125 131966 45180 95449 133389 141753 53183 99049 99145 21098 7417 93677 114964 14743 148954 92800 114245 74338 88326 47683 122787 33540 148317 107839 26104 137416 129502 139120 83391 153814 9030 85028 119892 156610 116989 23906 132475 4284 59467 2244 25383 3396 95922 76859 18140 78988 106171 132385 89838 28609 14180 46737 62150 99009 91088 24766 157102 163887 26953 145028 109429 111981 35545 100152 63083 59451 3251 67367 55430 5496 92750 58826 37930 106122 76966 53430 46405 145863 79780 75015 96555 63029 73677 132077 154117 98443 139127 81843 163383 142384 112173 33040 13160 29003 16411 13248 84433 124307 42510 79772 162237 148632 156738 152179 131550 136714 66071 143077 3893 129100 50866 135970 53841 85821 109209 135684 55288 45913 133180 134184 94599 102732 147344 40479 131736 163755 53727 152681 122175 32750 66565 55036 117894 159816 41328 83556 67154 107399 60745 71047 7124 48123 41130 60965 133945 150339 30762 23345 132765 163942 157529 61476 37298 75497 101955 3146 9876 92195 92340 68563 124945 95417 123600 76951 89345 101440 97020 156499 145351 157765 61659 152475 142401 102789 47553 110458 23752 14827 156517 115281 164356 154505 152579 73965 90573 92238 83841 119280 18690 152405 78337 114107 46629 91800 139965 84581 125332 130576 64044 53722 26347 153032 30235 65648 34697 77205 25913 84032 16542 101317 82500 5160 92977 14358 113330 151525 135448 36359 20180 19500 28798 160145 104081 90643 61345 2238 144365 24205 155270 8712 89853 126479 22429 115766 112515 68821 47944 125473 10493 72464 52562 145138 23335 102404 130775 158783 75276 67773 115475 94776 33083 109732 32970 123726 107590 137608 102203 131795 126990 47427 158160 24093 69856 108039 73120 112829 113372 57576 72414 123865 130040 61488 103115 153375 163893 4514 82783 73281 8799 32370 104569 41882 78615 137539 109259 38436 88624 6873 22376 80896 133366 28744 23047 40598 141574 72931 98174 48100 133308 62327 109589 70535 152214 44106 11561 69109 117387 20360 37992 158468 53119 130120 100876 9948 75824 98572 119277 98200 13580 86755 126945 139028 127354 102631 148471 59640 150731 115892 58479 30944 122939 44806 75050 134501 50427 128949 91373 24931 121530 90128 78050 22274 127516 87998 32277 37453 123082 151554 135653 73175 72422 33222 135853 28898 30040 57209 81302 25032 88153 38354 6350 99716 6967 56777 34852 18221 120819 124981 96271 143093 23121 120782 111883 60574 77976 97549 132740 87663 106483 133978 76883 135928 99389 43436 129649 14803 68468 154314 53157 74818 88142 162524 23944 31489 86328 144763 92982 119112 58481 116103 74006 4476 113190 88494 38537 80042 10270 145021 80116 80760 158416 152556 116661 35964 116317 67977 40944 104744 57646 115762 129398 156683 120382 153342 40822 68730 51778 159934 127211 104393 4564 68199 51695 29571 43248 131737 39841 22381 45966 57113 117310 32634 7887 89786 85464 130731 24320 70022 17117 90230 60817 137499 19613 114834 82927 71391 45392 146650 9897 49957 151361 163992 79528 28721 129842 55881 153503 112320 112994 104925 144954 57393 28823 64530 69769 96066 25362 76304 113184 115593 21307 66302 93247 136141 85741 15646 2115 154376 109152 18155 17209 9618 74036 107224 121938 123542 46261 37516 15048 102047 21329 44175 63921 97633 157359 116026 5379 12779 16440 98626 85432 102181 33888 37590 105196 39705 93684 127112 160260 47406 136730 68408 154630 29292 128463 137404 66808 80023 85513 2967 101352 129688 3401 135498 57672 119427 77389 70451 135868 112527 92395 8673 146415 129986 113870 164375 106203 41666 125599 100576 89072 32953 105496 14327 62245 68071 151731 129054 71356 68533 20071 35156 92081 92828 3582 99791 97854 116109 26299 106527 33148 92797 156909 31636 135512 32688 93747 58272 63213 112209 72599 125458 14392 160842 25136 98999 30182 119070 101467 38628 45263 130807 64102 151245 64711 98468 101466 141508 71126 56831 132529 39274 26455 133022 135568 160002 132747 81889 3225 28829 32793 164068 53966 131792 64490 84148 21486 102469 59288 66749 3901 123391 154506 68612 158371 90084 32646 133991 97723 25989 100633 136997 154844 36474 40643 124524 132988 7502 40525 136214 36332 73319 70906 26810 141623 71908 47470 163110 110890 114791 154990 119915 159145 15698 89073 63760 113421 115062 164393 84531 104018 137380 61686 62655 104480 69189 39692 11318 42033 49523 18736 68843 25259 90645 116313 22481 138047 159583 125364 86950 160357 114466 143377 153614 130165 66562 153886 77698 15736 152392 98741 56266 123884 160428 118921 95126 112707 42274 144649 131443 47629 4020 158600 100454 128901 130759 94149 88378 154221 88618 139356 131710 12856 103633 32384 117844 150520 89759 153097 40899 150155 147637 96332 85655 84390 25570 32474 63176 106732 163918 110805 110752 93142 147771 10278 58014 76033 98656 46347 101163 72124 114570 50532 112270 83466 153787 68099 14170 153835 45510 164325 78354 20604 156487 103925 53079 53776 147169 153509 101093 92033 80763 82977 38823 138777 95522 73991 121637 133197 82628 6831 120241 29010 90297 108140 29748 94908 58822 120095 140418 159659 26304 52885 14376 16904 93322 67455 70680 74603 157476 108285 103148 72351 127774 147641 57408 50075 39902 24718 69215 46733 144959 34737 137030 23724 2563 82546 121093 142981 76318 147397 29978 27206 164301 59812 31173 5605 134415 22761 113891 71676 31624 75777 150160 115777 69698 136747 92191 30928 40075 138924 9999 74812 46579 136123 75810 49142 155182 31015 128635 2124 12524 113338 154938 162903 118943 123466 122176 3458 131654 153800 15748 52438 40202 85446 23297 132393 116374 105430 62885 11209 88521 135521 23531 74175 124815 54546 36922 63451 132048 29293 53432 121099 26308 6488 15189 148484 112346 83355 72908 128094 135793 113110 150052 159090 79616 37050 95487 121558 99936 106696 44191 6081 130227 118366 130896 121285 155288 130859 61379 121449 160152 114812 13172 20572 57812 28361 105568 4270 111716 12588 68877 18133 125699 155441 11335 39427 26604 43334 160985 126540 150030 69133 50881 94166 34141 6278 20078 101512 4170 78039 95776 118982 91211 52860 113306 56084 158428 54088 104312 107528 59477 122445 67339 49031 70293 43278 75635 113627 140775 138687 34282 116575 41932 85163 147253 12585 27954 1444 114097 32124 15995 146385 87618 107207 33357 99803 128297 89124 38228 69937 85114 97185 137277 134145 1590 14667 43892 115218 155443 119091 149500 42642 161023 5287 24008 7720 33241 127852 58329 1877 143847 38826 89495 85166 72183 124531 121482 34592 150168 159710 104529 69394 27519 12430 140051 29110 129498 18055 80840 55565 73658 64452 98207 68793 6251 122215 13025 39493 84179 71354 143770 46692 67378 118875 128421 141811 89979 49213 72145 131020 55119 99665 143451 29282 128775 43573 149737 43727 99138 57507 44691 133857 62812 50942 26697 75837 26947 47388 83703 4830 46039 130395 8720 129857 19894 137141 105780 109873 20467 114438 75006 12098 48215 154969 41381 11102 32654 127630 54829 68304 19250 36032 36273 82062 23486 62970 94412 50434 110359 12227 157664 156398 79135 56879 99029 74149 99172 43015 31128 47722 54533 43227 95937 146014 84608 43551 115180 46350 98380 17596 2112 70924 156269 20687 94410 53352 115099 81356 100223 63838 9644 90733 79485 10140 84125 12627 20802 17409 158042 51930 1643 149087 95157 97580 129213 116277 141131 15017 99140 10135 135013 101252 17571 125394 58451 115258 110062 129850 49593 110413 76839 24010 86147 160964 36637 106949 114885 28791 158879 53040 11990 24661 87132 77715 140938 92732 10702 9022 61857 48467 128993 123764 106918 11598 73134 51093 141448 59240 161506 115055 136079 122028 137714 158666 78775 16664 123969 8278 69704 72472 135339 156836 86699 110390 155723 115944 57604 101257 114313 106071 72189 47102 12472 81836 34707 153921 141076 39600 111267 88865 11426 13046 26713 132110 78637 34991 138326 87621 4443 65786 8432 51345 158021 124376 108949 29902 9314 151533 30776 81503 32747 163339 3966 133681 138527 73203 7393 20418 162068 121219 135864 116736 147932 102086 131885 17036 11036 53618 123879 76822 9336 68955 54797 35369 68765 40442 66145 86780 9701 45906 84232 116067 13699 159271 23382 21093 116202 121963 142312 86178 72811 60869 124777 38808 77905 135813 28938 35896 149148 27500 147632 154615 86951 138941 24097 155716 13496 90242 76609 125597 136148 97353 86360 35671 107453 143450 157634 20389 64557 81258 125030 39877 95675 31467 5327 68083 14727 32828 49827 3455 119779 125281 129952 46120 54306 122729 34998 126967 156594 47287 47439 17954 82959 91404 97917 74705 48305 34682 75775 66076 96224 52164 98263 127691 57492 78931 90320 152686 18898 146611 112079 148850 26843 21480 139668 86084 56479 37259 76790 40278 84698 123237 112614 65686 97432 44342 163508 30704 33018 32395 94907 90510 32854 7950 117342 19653 26848 98066 68244 112210 61421 89725 85991 84018 82716 123250 160808 122994 144461 128577 80344 27699 28375 48911 125131 128743 29765 59263 95959 62783 91659 24979 89806 124513 32929 41260 80678 162178 75838 148923 108500 73772 72760 28603 157790 91988 88366 152710 49094 66939 51912 65950 31150 80287 51373 156282 145543 81138 49657 75614 80434 141316 4352 36454 70035 45612 53644 66325 57963 36679 111337 131735 45951 76453 123637 74451 46971 123546 2382 98883 23608 33532 115683 74982 23926 95338 92632 10096 107464 7178 87924 144570 113930 60890 151117 159543 114535 51554 154018 87726 99403 119865 70190 9968 14126 144641 111299 61097 38811 113681 96493 83726 46288 10426 44164 78138 103058 54260 122114 46749 78697 100796 160679 139587 86025 90846 24746 74091 78976 48985 7607 135353 119175 119975 149479 34440 65387 147089 73252 13180 77694 8696 33418 60494 19122 75144 122180 68355 31370 105441 132167 36745 57264 154704 127591 18522 62908 40680 67507 7027 112545 123194 127002 96137 157635 26501 77338 1511 142082 91544 10207 112012 152038 29329 23707 61294 88021 28575 92664 129975 112139 3515 128423 124438 149580 105615 26057 93119 52273 149252 54233 84922 77511 17247 162260 79022 95841 24428 89229 41965 10578 55070 2185 8384 79603 30760 37560 43690 135416 44522 56200 109458 41477 39892 151585 67534 133011 37970 153298 123756 59404 64921 77515 55776 143943 7468 80204 3796 151834 27294 161266 154019 138078 121291 9751 55184 98302 145167 36219 15341 35479 145677 159218 11883 67886 60865 81406 42369 48275 101773 113197 116790 94062 27764 60770 8378 133961 149116 138073 129339 73759 110263 38433 29162 56526 30130 127465 35805 66349 79318 7796 148538 72649 19679 50536 70026 101085 92905 118301 76872 2122 130174 122853 132287 25056 67743 100360 110685 142328 63811 18556 86704 38757 143230 68887 111696 115548 71748 25126 123344 34287 143024 41446 104313 14733 134352 159127 91605 99655 97761 55891 159020 64160 17460 60875 128145 139716 146701 60532 140113 130932 140274 43112 76740 26446 51372 38379 105770 72666 63439 83728 113492 78172 154592 106731 106290 88359 141004 162181 81491 141677 13753 78879 76821 141898 52707 59224 49751 50002 24268 126537 29626 101009 89496 80998 75900 73399 148566 73350 143521 96170 88034 132225 139413 28436 54697 51042 24730 72700 26831 140883 88091 103652 77310 115942 112630 63573 102456 136898 24223 68595 8531 113719 86105 84431 123630 145898 3622 133492 123531 36304 55638 89868 112230 20587 81077 163272 147717 90290 24215 122712 12493 64379 89804 116833 59372 153377 53402 132783 114112 121997 141314 164343 42214 162258 122085 124624 26201 82267 36296 18352 108648 148526 38939 23837 82422 123168 114127 43149 16504 63133 44041 129135 89449 160874 22619 76938 48388 155402 127563 106897 67341 126018 85624 18728 44360 44929 138911 69388 58062 38732 29737 140484 98412 80377 120146 114917 143510 164187 78164 67071 159173 37295 80521 144074 129210 42196 85083 30663 104727 107219 30886 59967 151580 133501 41408 48832 106524 41202 54732 101098 20981 75695 114859 101358 29953 63888 15492 28252 78564 19075 123937 115859 99597 102123 15693 78305 123719 46356 119544 65050 13754 13623 153142 83767 55032 138486 126804 96234 129730 62014 112680 150712 137709 164051 22694 104174 140587 68938 77127 159662 26988 129499 29883 65623 145192 44701 23454 128061 100757 88505 78327 114381 12271 162095 3525 150758 123011 36271 51112 19137 85463 138336 93359 83627 97543 134045 82190 72242 139496 95830 2528 102996 161841 134812 2731 77657 13625 128586 39830 50894 153603 54670 105312 110400 78095 138369 66049 47329 95920 17810 144872 36512 153626 140086 132342 156155 77194 64807 125079 16437 142465 75216 145023 118807 132738 7589 62010 69237 160866 103834 156493 41890 16042 65486 107940 116113 49328 140372 133923 28312 140962 4547 16051 115160 73402 108718 128866 138209 4421 114786 16149 163351 4218 78771 130202 11807 77293 33551 109186 15240 126556 151076 31282 128555 93128 147395 11995 4125 51942 142707 145087 56489 95270 30872 66403 38100 96250 141125 144921 12178 26535 161070 9642 30753 10466 76356 42561 87759 46419 151747 109487 136935 70793 72154 154700 20661 107213 40754 163369 86412 92751 53796 100159 67364 150047 75396 46397 98737 101932 143980 44891 69197 154446 121247 111758 12829 1778 34129 52341 111266 107577 123134 119932 44777 111947 140594 151990 89213 74587 9026 122969 1450 62823 159641 5326 149382 5661 51724 82231 44105 132216 127123 113303 120774 161573 70115 20773 29815 122456 68551 73904 79703 124995 55193 128162 99701 41295 51487 8400 152721 110969 112251 152056 104722 135550 46895 105813 51894 91001 72141 115529 38416 129427 134523 34101 33655 91808 156111 96871 74332 6438 55979 129525 134600 92192 107332 122600 37105 30678 67681 149356 16846 108915 37557 152397 155810 79883 140803 17435 152024 26956 55851 115564 161479 26465 85731 23912 26893 75954 57295 37738 82393 113274 103775 153505 141979 147620 110217 13196 12410 114410 99064 29256 159837 73133 15765 86272 153016 93080 103707 75665 120036 96071 127741 52140 122536 47584 76052 85941 60050 133347 123679 142443 17246 163967 66573 95737 145699 10902 108933 158109 61825 144509 21477 55774 51754 37243 78558 141283 130323 16378 51060 84472 48961 12913 73124 5609 162897 149176 91550 57059 53147 151742 136015 6905 149821 36700 129632 47602 148087 121853 45939 126708 143330 38226 114975 117085 116784 26882 18033 69674 77942 39017 118635 27367 112141 60756 24376 31941 88819 81435 85088 74673 51562 28506 161006 88262 131148 61262 72377 113348 17227 54828 10680 97069 93054 125655 150667 43951 152537 105212 50137 2885 28470 26994 155246 58935 88972 7306 80535 100157 161268 45553 95275 20155 10814 92532 60674 110276 83872 143858 11039 133522 65149 157390 120171 4473 41639 57787 146924 22769 112208 147078 12128 107655 70162 121834 124702 4432 117214 4368 137369 117582 129269 12368 10169 84398 83722 154027 95437 153756 155688 23451 108040 96673 65090 102339 15999 46127 125108 64719 129717 73748 108886 136391 29694 4213 140823 146909 8581 113555 54902 62675 76936 55428 75043 87105 139826 95277 11757 5887 83146 1557 29338 127698 98231 30940 64149 50742 13579 125769 115462 143296 136030 160860 113800 25247 83257 110166 138802 74672 6953 152250 130100 81996 9980 113786 21737 148837 31044 125694 12287 95254 160437 95915 45292 159207 55796 97266 136615 21039 128175 27936 46287 147944 138102 121601 56728 145056 44476 123340 163564 54456 163890 111462 12705 146839 142506 74911 95638 69461 75378 104299 51265 68697 160096 148531 141824 122546 11393 162864 21345 141729 145663 3402 113944 37888 162530 93112 82364 119983 90788 73332 117985 36363 86037 35449 115381 160949 131087 123766 40577 62177 52299 45786 158785 137342 21723 51956 85247 73301 61089 100874 9145 138762 9858 55239 129841 66069 128571 38944 151121 53899 154326 82694 19099 48716 162076 47405 72539 144375 142703 146393 28771 65752 121309 102072 63353 158695 115288 72498 131569 125146 47779 123320 89099 113849 107549 89305 7749 161449 14255 26955 117060 62971 110683 49760 46889 119734 28247 23704 100239 130857 154875 22790 27571 63457 15597 142859 72467 83679 38629 56759 41111 64240 107120 63627 108302 30537 7888 40375 44793 137243 93947 44276 18550 143707 91165 138284 6067 51381 72635 73436 142768 9251 96227 4452 111824 147311 18199 29615 122453 74958 7239 123205 118590 7378 65619 149128 117667 42506 130433 89022 136453 8821 107573 114273 36499 79969 56852 24392 89117 130288 34880 60627 5725 44100 108964 89548 75091 148513 10282 18842 105831 128872 26221 107962 48624 143888 150469 13169 3534 121034 124391 111107 97402 127589 121794 50818 129584 59579 22210 126723 1816 6311 91365 1633 81403 10502 11915 36757 116333 77300 62978 125924 143378 145389 146913 37047 34109 42979 131511 156233 23929 143563 108370 83509 69206 85325 2707 75517 113202 4341 93432 123705 118656 130190 10662 129680 155993 43683 87483 65711 127108 99820 155752 104022 1955 82808 25884 60285 25291 45905 162571 94497 106526 117457 106131 136471 75274 58900 100773 85937 88968 28031 15585 147160 81296 108380 10670 117628 34757 114692 119583 117367 31612 81980 114164 159303 64397 47359 24424 132140 46750 130951 83710 89393 38046 95496 148293 138819 117945 71374 145977 23471 57681 95604 104768 102573 106274 56508 137330 157478 112604 25321 125602 131096 75998 119018 29605 123357 79954 161746 106619 147417 16080 30125 19575 111576 114930 94906 122816 5427 146123 146288 12351 85168 102193 55137 78188 10148 149127 127304 35469 108841 92512 47979 164371 122118 107849 14950 54488 48580 162367 7080 15217 16055 55168 130148 110961 55314 87076 52900 137949 67476 52520 86812 89156 91226 141949 103857 101374 61700 65273 73355 4654 157786 121335 105537 50528 63296 120487 41528 53479 48608 63606 69534 103776 27866 117007 95603 51454 106420 67664 55443 158940 154477 81111 20790 19080 122164 128751 84354 132032 69917 12764 87479 9566 87287 66566 41332 135675 89940 33393 126091 130228 163659 13722 59943 151626 120142 127608 41181 49707 52709 122292 70497 119759 141373 129174 82622 162239 95318 152539 9115 119309 111319 40708 152651 10495 55853 13215 146289 118456 143443 80572 68690 37499 22945 101619 44003 9164 154328 102808 16173 108199 78293 145347 127334 11156 11289 113985 122671 130598 46715 68102 5418 148411 54865 118314 68080 35227 29856 148035 115799 98547 122046 16222 58004 57777 162625 67168 148617 99545 83341 27441 114350 154775 125506 10602 102872 18801 77713 86100 86903 5135 78280 137957 145912 82872 7296 112280 65020 59608 44939 21178 12342 102943 15468 9079 106623 164085 108625 126477 128038 57087 125790 116925 119106 72904 156310 137907 150617 76522 161323 4373 18169 10227 78842 164081 93100 22651 110474 158120 91925 115810 94601 131381 131278 40192 72116 65988 148817 32705 130538 142417 95007 81576 32147 2024 8510 106566 152641 21544 102001 93526 39713 112229 6481 37907 141841 29132 84893 70585 47903 20507 142504 78823 88298 119208 87452 90798 102138 120157 55448 15179 151677 137024 149726 153701 82046 90404 103590 128918 143304 75259 11086 117723 153612 142618 36728 58309 24633 150058 78816 103649 56969 86958 28536 95619 143627 83985 47310 129416 157521 31148 117229 73680 121553 28308 11382 84583 135313 91198 96354 82911 43033 84078 59641 16273 142387 20786 157715 60947 102350 95454 117916 23420 123990 150048 144487 31470 136121 62619 43913 120684 149137 157695 39379 118562 83006 72245 35585 126039 92835 31738 142313 5846 52524 79269 163561 49983 15731 29639 2012 141551 153630 88572 142710 132229 120042 43910 38974 119173 82887 73969 147989 11207 49860 37176 94213 122105 72761 156764 49052 41012 69701 54898 30048 85482 154972 80032 101213 18723 108865 107128 156099 11719 61282 136521 50693 116968 135373 70093 25049 117475 17812 11422 91163 48537 70039 39413 119092 45627 110502 131109 36098 22153 2947 54822 40709 16335 100199 147837 108946 48430 145631 79580 99124 96711 151465 105729 58273 103052 60053 69695 28327 108590 139734 128764 84515 29450 40836 130142 76465 46445 95364 112563 98311 103897 45819 114646 140609 130168 57705 23151 109912 73797 58787 143135 59374 137920 98939 58681 144127 63778 103783 54486 22411 20448 4003 152553 96913 152848 18541 145989 157959 116852 20510 37890 2123 161119 104571 162228 18383 150995 70137 77170 64754 129511 12311 36787 62562 7504 117426 126341 111288 6024 89507 70211 128873 30022 56876 162298 16983 75418 78911 128782 99422 149344 130905 94653 88027 127245 113036 9646 133894 126719 74401 97518 139030 47700 96592 83046 57045 28446 107663 83065 98657 70648 113087 92046 3570 66582 103976 82482 75636 66870 16016 161493 134288 47181 20144 32157 96730 29791 57561 104192 97682 133103 88404 130786 50261 87642 21944 15220 131817 41521 50389 138977 91120 142435 142548 157703 17035 59142 3963 83906 11670 154706 161251 19713 123375 92093 151904 123539 86166 26720 53381 115124 122632 40155 139278 81088 157775 105208 122609 42276 78297 150241 121224 54957 142056 138259 50611 146020 158677 145588 147495 160045 101813 104983 86250 24341 6529 51061 93392 162310 102697 152536 76088 130864 15408 126857 7253 30217 111210 128477 85175 87379 100848 72298 4023 93638 71092 149611 65249 53252 88011 46390 154799 31052 33964 82303 129875 136661 68951 102964 149261 136687 67941 101181 34168 11706 131399 145378 140183 50686 3381 11655 59496 7404 130588 93528 155963 132349 115577 43327 54473 98479 32238 85525 114541 49512 120004 88989 125502 27315 156930 60796 2748 26307 143374 79443 13505 146755 27610 73001 90672 37702 18312 161891 106563 70401 39330 161036 105392 8081 17186 8460 59134 3210 150589 13251 92199 40566 19753 107511 38562 124901 70330 16048 140856 20347 162803 93349 87587 110907 67563 42411 106910 8238 112812 82752 105787 52317 90833 59485 163177 149968 147878 99731 27919 131105 76809 47672 72729 115371 109086 79571 131419 20566 99919 125545 129780 88946 70564 31455 131358 11586 142093 14794 94339 81992 3623 77989 105764 77197 85302 142007 105116 50520 52929 89300 59761 2412 32498 139332 70344 9875 71703 115122 139655 160649 19799 62631 133785 83828 13938 102332 17561 120433 14434 120874 162709 28143 40288 133259 27320 128270 56671 67194 25682 162926 137538 15259 109313 43353 66893 83081 140515 23204 127215 37258 93502 45155 51196 92250 84000 5270 149195 34946 62656 146017 63089 102944 157870 132860 66777 120252 126043 90370 21558 88237 23920 95609 103496 69746 75474 106902 152827 50101 130106 50666 23872 57720 95822 11580 86483 116334 69790 151280 79506 152319 150881 118963 144301 117853 122252 98665 14520 149134 56736 102758 7167 88857 142766 76913 164331 83780 48545 47999 51030 72417 42231 83364 20509 128714 136210 37360 121602 53378 121448 108853 99862 97472 65217 135039 111993 48464 128287 48863 158031 51256 25741 5568 49700 5932 94033 159362 43686 72776 13351 64195 35603 149561 38067 4732 41788 91446 62692 84895 136811 99066 16480 36141 68218 128473 21117 133017 113848 115660 18386 139589 121228 4598 19735 63672 39355 113768 57146 83041 147236 95172 156571 19416 99904 32471 47374 99108 53878 18297 32287 70358 156838 37017 135343 114467 4146 83303 64239 22532 159405 121979 129530 13252 19763 2997 63533 76910 22550 20702 3899 115874 125715 152290 43585 109601 22022 33976 127898 54309 104334 55361 91326 10302 106340 31985 4692 156917 23634 126671 120560 139287 146435 123557 159857 146108 96546 14400 86519 148932 25515 148746 135334 5612 92459 157356 39588 156869 148178 80435 46342 10128 27249 89195 42113 57366 93887 33143 17513 157070 153703 93312 74129 47884 66756 4610 28104 163302 19011 51135 146346 146926 33993 52304 152538 126452 146172 128639 53946 64974 145586 36800 75103 6947 125995 117216 156394 150359 147577 74686 48162 158218 122571 114918 162829 87187 48844 118352 138323 131702 99390 108828 120518 22552 5905 37315 87703 59851 102289 67401 11504 74348 159159 65233 75174 149665 152104 30024 67866 60903 78187 60197 119986 129617 57138 41285 112002 116120 80788 45504 59061 35419 68056 64966 72734 155760 61329 111535 157980 123040 68134 124785 79820 71541 109844 139407 108416 124543 136116 62514 88273 129766 37359 75880 54032 118147 121384 49605 153566 125953 114571 162812 52337 10012 44972 46122 104505 50768 128780 63421 34433 8534 11057 110353 53463 135600 17093 115977 160385 146860 156289 31856 159252 146834 150004 51261 30551 74194 11326 81635 7631 28159 52603 148697 69176 157127 135978 32069 54660 37190 3014 48247 113367 53169 66973 3259 114857 47945 96060 146714 143709 77006 67342 29082 107558 141536 142808 125705 85679 153864 138282 62305 57153 129521 89222 120693 6097 92236 120400 156745 142116 7681 111217 43201 113340 162586 27658 145913 79924 5479 57031 147266 34562 101101 59427 113882 60918 145106 90978 48895 117501 89795 106048 17646 58702 131782 138340 64799 58130 92852 55656 136758 37045 103385 116471 86898 100084 80641 66923 116520 22633 123954 97899 159595 161568 93838 107589 56598 9568 32679 105493 127069 58986 148053 81228 117688 50459 53680 16600 83044 8768 15972 56601 112154 132444 143499 148750 47197 146934 99382 6342 41513 133793 37193 61743 107639 93791 71312 140318 135797 134893 135817 54474 50233 87617 104934 40425 40729 150036 123469 49498 102520 14183 161652 69076 94194 81026 75241 59128 116754 148314 59178 90459 44169 18362 51993 77759 90671 123399 81658 110144 71818 20704 150570 112547 107252 44663 98557 43884 161246 30833 49473 89553 111859 38771 101306 126380 52172 108806 93589 142631 89487 31030 97505 77993 45302 124689 35504 29502 40459 82160 37832 163563 66842 150380 104927 48018 19561 85323 43376 50395 134796 69441 162254 7679 4859 122746 42091 159431 2176 72192 72686 28176 103223 106704 106169 85037 2017 78185 51051 42476 160345 88884 142551 163700 9888 81590 45830 29449 25718 79844 72334 95160 76211 80013 36531 135469 122104 137646 130809 141673 4656 82489 42615 89693 84506 120800 140744 63494 115257 40157 113069 58259 95411 39590 161685 121130 119434 152802 29757 84656 125845 101739 143273 92432 73497 108194 4730 9943 112850 23731 52558 36655 44750 109870 113911 108244 59239 84914 108821 124304 38344 163894 75482 95986 53953 143613 19412 20222 62382 145258 121961 39767 8314 131970 84473 13044 141913 31435 36776 130983 4602 18038 74965 118513 126282 70717 3301 45308 116050 127606 124994 90906 62124 34588 21005 116077 114714 40417 72812 113608 122187 28885 130502 160856 8472 80058 73393 39907 53346 140889 44509 71384 152366 99534 31779 102836 13599 109757 75105 34776 46206 158181 52292 160920 32710 125104 108640 91410 153989 32639 56024 85469 136082 95375 81018 125941 6888 62039 31437 159254 35033 115728 33343 36099 24945 4631 35802 71151 88094 2695 132035 149710 111335 57557 137811 80486 50093 159793 121597 20287 89280 82740 32680 137678 114178 26046 133364 113906 115520 2509 63746 85376 141009 25204 90007 113323 32867 23444 137929 35563 155479 121751 83410 47148 30186 163897 97241 24091 56118 54040 49883 131758 136781 103548 21583 45122 7536 72001 97154 112446 135748 119042 87567 97464 145562 137402 130331 3118 102406 95109 1818 19929 78769 32004 120338 112522 158496 112968 144891 78838 16499 32062 145601 126435 104064 13379 72993 10436 68934 97072 107900 68586 72343 51726 8862 146835 116359 131191 162116 84876 108750 13996 135121 45936 95399 48072 56416 29990 88478 112103 127518 26654 125483 139490 28529 68208 81502 77137 136795 90358 65375 21960 35732 148722 138320 144950 159708 125588 87812 110216 94821 39727 79405 24276 6230 80488 49267 118333 144519 12433 77928 52167 42969 120375 67538 56618 27794 94408 121993 49755 130140 41339 124587 67687 20401 54919 27387 44725 122209 114986 141374 35726 75774 147604 52727 125041 100050 137475 114490 20037 123511 157460 140412 127561 48190 104719 56082 90986 122734 84547 49685 24534 41461 41116 154321 86186 40363 47154 137684 118249 19393 47571 36478 72120 6724 136528 142390 80711 21642 162427 38334 115614 73463 102408 100317 114694 158490 143525 39792 51848 62184 76382 103645 67105 107468 48482 98140 162230 67875 145711 32820 76507 152436 155409 3771 127502 151948 144506 13741 161924 91695 87035 61988 153879 113024 74883 28149 154493 58222 61205 86745 54564 65593 66569 23896 78612 155517 129756 159288 27883 137906 41624 89142 155444 92194 136385 18001 57532 4668 73080 73053 121310 85803 61658 13645 147008 148404 4721 149113 151485 28618 61838 141114 158374 129972 71026 20369 57190 145445 146333 135001 90699 116830 71210 45212 110874 69407 153605 52054 142487 7477 58941 123522 143857 128243 65660 129454 156861 127498 41192 85860 48730 42741 21796 28287 2241 99801 68458 55583 107523 113671 11042 37900 90041 35182 162501 88955 94123 120135 5501 72092 19002 71161 138058 112376 135171 13362 32348 99777 11601 121574 26892 155248 60327 63205 57819 115911 4840 116480 15882 145902 105616 42435 117596 102229 131390 148232 158876 73404 11991 129007 124367 152761 78881 93227 23162 77486 104828 144736 104379 94189 141576 1696 88520 91599 108937 133034 144591 61331 49560 20832 103767 59573 5781 15697 76961 64553 96775 40080 125432 20160 55473 43322 132960 157599 124571 71451 22412 100259 9659 110932 128370 118596 78078 116440 127638 64417 54319 65419 60503 162500 151431 113065 12310 62504 114130 45597 39096 73674 65757 94570 126445 45591 61642 118156 106674 69605 77081 41045 79264 124525 105927 134373 139116 149513 84925 37378 150442 75756 102797 45057 72368 88341 94634 21190 87357 42876 3299 62966 116550 157536 13619 51160 155690 68288 157834 59407 145369 135392 75183 104006 75431 146068 77234 161456 1617 114613 82523 77373 51522 64092 86254 76375 158727 43956 100245 35715 47256 163211 88778 52824 91371 102397 40496 81173 107197 32443 77092 86678 104347 88787 24797 116290 68968 38543 153156 130894 84471 141191 29098 107237 154078 21937 87705 88435 160053 71473 22270 19455 124298 113641 58364 131326 133749 42530 86352 131318 47661 125011 86205 7494 143895 9985 138388 87688 1598 106227 75879 125936 130445 100826 56613 36030 123097 12580 96840 7362 70944 32259 138689 7130 117731 29995 111241 83161 77656 70364 54071 100373 77858 134479 152758 133469 32159 90868 10321 44550 153316 140766 145377 146441 10908 39098 159021 44261 148860 76520 121661 13152 73505 37715 60905 93178 51884 131269 83762 152257 145639 52353 33240 69021 19934 1912 159890 132655 46462 83830 107533 128351 64384 54954 103961 57517 99215 86934 58107 112247 145107 2340 125399 155124 142455 20417 82415 28451 88198 102689 117221 4462 91554 86973 73483 48000 88885 3997 14768 71860 87828 122301 136723 88724 113767 74797 146241 149494 161731 140861 95853 140950 79713 130187 56280 149114 21244 37004 22427 74977 143866 113981 161950 153862 161981 21460 157859 113261 93320 16311 6187 64155 105035 119954 75464 21901 103561 71307 99274 135926 115499 127803 49581 5891 135207 129535 50126 38624 42239 34687 16078 140732 56147 110450 88106 85979 63273 94293 150134 104821 48359 59711 126722 88432 67530 162508 160871 152812 112119 122786 138905 54522 135908 44736 18170 20146 36482 56794 150650 15233 27993 103241 155966 84140 47803 78184 4231 47588 6589 90878 152409 157348 87101 113243 16405 91143 46375 11388 78068 158494 134174 153485 47129 32334 1811 60852 5328 58605 148015 20562 86599 85368 10640 107251 69683 88824 47995 117271 31925 138873 103793 23385 60086 153548 142190 87741 34036 153578 102321 129042 121865 89919 112683 122253 114494 159936 64093 109612 142063 84655 61543 95295 74086 131226 120631 122081 19121 89068 31578 122914 112454 91664 110575 88756 115918 81123 12959 52351 44277 71336 142270 93473 4531 35147 42079 99241 151691 77154 18008 118526 138697 49816 26725 40547 106959 85318 59669 30140 116897 119095 79106 145073 104374 95103 81417 117333 147455 62207 22781 123837 155680 27313 95497 31871 128292 28850 120075 39559 149258 72713 114768 33186 35950 91828 140146 57781 151497 106798 8790 41216 20016 153863 41511 124390 19591 122928 12348 121647 35129 61507 111439 157004 79823 127247 34010 105737 43235 52198 75435 144026 49693 73479 66147 93003 18781 107364 113019 6756 148875 8034 42427 20382 129905 164075 157911 25525 46138 156866 119041 62473 118225 75520 32995 55655 117769 12743 107853 27316 28441 85992 77010 38432 152139 106525 57213 30127 53657 63969 13114 164091 26829 157942 18585 93246 92641 10608 118771 138779 1586 137765 56323 20469 4872 48490 138238 17615 92856 46056 115360 13188 21000 38123 119714 78213 68251 7483 78695 17877 5686 105524 9931 32882 102572 137279 88166 11976 33115 149741 68213 53584 154613 53216 25934 108741 146072 25600 95544 38789 112310 133667 158503 24635 138430 102498 103330 156308 44696 145366 107879 114761 39436 15782 139439 51412 91160 6666 137665 95886 60250 62902 149102 86184 5755 129286 48296 161342 23597 43487 65633 118612 4634 38176 55222 44477 28596 36430 23955 131347 144309 138716 4895 160092 48779 21876 55445 66996 54274 52207 129899 37488 74903 72166 103286 123200 3343 35252 149230 100886 99522 90377 75574 91256 134854 104170 64199 158809 6141 145020 68150 113437 75736 116929 106256 97613 108887 151887 161094 73776 125888 70110 145942 27422 85798 35051 110731 5652 72449 146765 96029 148023 64995 22817 72845 160317 28958 154377 142395 64226 116020 19276 98351 59019 129041 20862 156626 83262 60848 55895 83061 24782 141693 118112 135513 83858 24673 52903 13999 109208 61549 15507 132025 9936 160983 59395 111427 60133 94996 57091 95344 47829 68785 33233 45189 124680 116294 69971 36998 5030 141996 57368 132103 29011 71367 75423 27072 23386 143960 34490 135722 139056 30398 81261 149571 90531 79333 19640 72936 45330 72011 120765 114115 105244 102466 9419 158051 99593 85057 40297 2321 50581 48177 77744 71563 55817 48656 143798 131385 79054 59171 51580 106097 75016 1566 62249 147952 46896 134260 39341 97523 10129 141807 43454 104692 84868 128784 140797 106438 3193 143118 93531 153770 157375 107696 161846 149704 92864 139756 115201 108430 135439 103293 48640 44568 104859 47401 129032 88267 118174 2486 122302 64815 144293 3619 128652 144416 106224 68358 58159 33867 56240 49646 78075 154598 33462 5052 128466 113482 34530 22580 98634 15610 63951 146036 144642 152218 34834 147128 108632 99649 62046 45012 103268 61832 110177 18308 104568 12647 76467 138435 5399 126113 50622 159998 96087 158074 59088 15374 105669 37954 140815 148140 38417 120963 126895 24747 155797 44647 69891 26070 106693 51416 129338 105037 161593 84159 43717 97138 118664 52663 57364 105799 149173 89963 97985 42374 105338 140166 9016 79804 115094 157157 118222 6681 118164 79481 98991 99323 149373 125061 40129 137301 25024 81678 69518 109183 125396 14382 40433 14684 67045 34309 120483 124273 154981 29217 129259 140633 80039 78465 131902 32373 21659 20690 111855 120650 120014 31852 96655 3265 104847 12445 72783 150542 74353 87165 25088 89038 90722 161797 77565 56694 135126 43294 56929 98498 18040 136969 113475 149942 3454 135134 107145 26408 61271 147161 106232 94438 86938 147591 43395 159721 132246 117749 17510 108233 89755 21044 19910 82962 92683 139891 127693 17757 47484 146927 45468 102760 109847 106739 20546 50191 137689 107484 31895 15197 101318 164141 132946 118828 92099 110357 163573 67913 119996 35440 97108 119713 111844 58913 137470 159328 70706 75793 44379 154145 83652 98104 104198 148296 136315 45794 116703 88024 147112 51468 55082 143567 101951 34250 32557 106376 154246 4510 37597 108071 116354 33022 79653 46306 103728 91959 90686 91985 47333 110850 24202 45958 49160 109010 118785 94954 59825 40921 12691 111293 32515 49256 25484 134466 83506 160441 11467 71865 151464 116448 51929 18598 132614 98236 122327 58685 125434 48424 106018 70396 9138 151976 119556 118148 104874 151022 12085 82307 163713 59890 114823 149481 21886 19913 67100 31380 75477 15503 119356 26037 67433 137955 158651 30906 51448 61727 79330 93978 132123 24981 16579 22303 143129 57965 7437 91727 140272 88129 89207 110016 109121 94868 126456 77013 4457 78471 30482 82416 104949 62170 18709 135855 50131 80436 151697 144109 10790 160688 107886 90432 155165 51836 118671 129550 159499 40912 155269 21379 150928 98502 116248 111496 112028 57217 126480 24224 43008 18083 22907 61718 90450 73038 78666 76259 51259 164250 23562 148460 106248 113994 137737 158084 37911 88207 44201 27293 109587 29242 162347 77250 8447 53676 37842 87469 14413 17363 42064 37320 132514 110358 94259 42886 98129 92621 2960 80701 135382 53466 154951 127578 56755 26974 49898 100956 54267 66710 92456 143961 125022 146132 118315 149004 160545 135678 25180 31977 151271 94207 78847 16155 73605 11089 45288 76565 28302 14782 130031 17365 142361 20898 146740 128771 35119 161280 125077 88206 87848 39662 4604 94489 132857 101662 64280 133639 152063 22868 46599 2330 32985 160088 47619 109550 22502 10205 39868 141274 133503 40669 25982 158239 36061 151059 80558 60421 127233 64411 65026 55834 133780 56626 70951 145202 124222 54062 147532 157207 48262 131663 37382 70765 130577 47587 47145 42475 15203 4377 19657 41185 162616 158118 128756 13798 52652 26613 14721 54190 82447 85013 54990 75586 20479 62453 54900 144701 116516 36544 72532 101290 104720 109914 108567 5921 157502 92224 48396 96601 4565 150402 29842 162684 49782 146040 49448 76395 160762 103638 95354 16399 95140 36878 157593 59952 118091 44733 33009 24736 82536 74241 71162 25215 80162 53952 65070 108505 150553 69636 93019 116907 66432 97060 52392 92220 92542 124087 147245 21794 65651 120636 115899 125604 72839 97145 95125 97575 13793 105878 2850 141408 122552 2138 29472 21734 110644 116538 27882
\ No newline at end of file
echo "Input: $1"
echo "Output: $4"
# ./../../cpp_based-on-dabc-lib/unpacker $1 $2 $3 $4
-./../../cpp_based-on-dabc-lib/unpacker $1 $4
+./../../cpp_based-on-dabc-lib/unpacker $1 $4 > /local/samir/unpacker.log
echo "Finished!"
--- /dev/null
+/local/samir/data/
+te15021111151_DST.root
+1 1 1 1 0 1
+10000 0
--- /dev/null
+/********************************************************************************
+** Form generated from reading UI file 'scurves.ui'
+**
+** Created by: Qt User Interface Compiler version 4.8.5
+**
+** WARNING! All changes made in this file will be lost when recompiling UI file!
+********************************************************************************/
+
+#ifndef UI_SCURVES_H
+#define UI_SCURVES_H
+
+#include <QtCore/QVariant>
+#include <QtGui/QAction>
+#include <QtGui/QApplication>
+#include <QtGui/QButtonGroup>
+#include <QtGui/QCheckBox>
+#include <QtGui/QFrame>
+#include <QtGui/QHeaderView>
+#include <QtGui/QLabel>
+#include <QtGui/QLineEdit>
+#include <QtGui/QMainWindow>
+#include <QtGui/QMenuBar>
+#include <QtGui/QPushButton>
+#include <QtGui/QRadioButton>
+#include <QtGui/QSpinBox>
+#include <QtGui/QStatusBar>
+#include <QtGui/QWidget>
+
+QT_BEGIN_NAMESPACE
+
+class Ui_MainWindow
+{
+public:
+ QWidget *centralwidget;
+ QLineEdit *lineDir;
+ QLineEdit *lineFile;
+ QLabel *label;
+ QLabel *label_2;
+ QFrame *frame;
+ QSpinBox *spinRandPix;
+ QPushButton *pushIdentRandPix;
+ QCheckBox *checkRandPix;
+ QFrame *frame_2;
+ QLabel *label_3;
+ QCheckBox *checkBank1;
+ QCheckBox *checkBank2;
+ QCheckBox *checkBank3;
+ QCheckBox *checkBank4;
+ QPushButton *pushPlotSc;
+ QPushButton *pushFitSc;
+ QPushButton *pushOpenBrowser;
+ QPushButton *pushSave;
+ QPushButton *pushEmpty1;
+ QFrame *frame_3;
+ QSpinBox *spinPix;
+ QPushButton *pushFitScSingle;
+ QPushButton *pushPlotScSingle;
+ QSpinBox *spinRow;
+ QLabel *label_4;
+ QSpinBox *spinCol;
+ QLabel *label_5;
+ QLabel *label_6;
+ QPushButton *pushClose;
+ QFrame *frame_4;
+ QRadioButton *radioHld;
+ QRadioButton *radioRoot;
+ QPushButton *pushFindNew;
+ QPushButton *pushFindBack;
+ QPushButton *pushFindForth;
+ QPushButton *pushUnpack;
+ QPushButton *pushRewrite;
+ QPushButton *pushEmpty2;
+ QPushButton *pushSum;
+ QPushButton *pushFakes;
+ QMenuBar *menubar;
+ QStatusBar *statusbar;
+
+ void setupUi(QMainWindow *MainWindow)
+ {
+ if (MainWindow->objectName().isEmpty())
+ MainWindow->setObjectName(QString::fromUtf8("MainWindow"));
+ MainWindow->resize(595, 336);
+ centralwidget = new QWidget(MainWindow);
+ centralwidget->setObjectName(QString::fromUtf8("centralwidget"));
+ lineDir = new QLineEdit(centralwidget);
+ lineDir->setObjectName(QString::fromUtf8("lineDir"));
+ lineDir->setGeometry(QRect(59, 12, 391, 26));
+ lineFile = new QLineEdit(centralwidget);
+ lineFile->setObjectName(QString::fromUtf8("lineFile"));
+ lineFile->setGeometry(QRect(59, 42, 391, 26));
+ label = new QLabel(centralwidget);
+ label->setObjectName(QString::fromUtf8("label"));
+ label->setGeometry(QRect(19, 18, 41, 16));
+ label_2 = new QLabel(centralwidget);
+ label_2->setObjectName(QString::fromUtf8("label_2"));
+ label_2->setGeometry(QRect(19, 48, 41, 16));
+ frame = new QFrame(centralwidget);
+ frame->setObjectName(QString::fromUtf8("frame"));
+ frame->setGeometry(QRect(198, 80, 251, 31));
+ frame->setFrameShape(QFrame::StyledPanel);
+ frame->setFrameShadow(QFrame::Raised);
+ spinRandPix = new QSpinBox(frame);
+ spinRandPix->setObjectName(QString::fromUtf8("spinRandPix"));
+ spinRandPix->setGeometry(QRect(30, 3, 101, 26));
+ spinRandPix->setAlignment(Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter);
+ spinRandPix->setMaximum(165888);
+ pushIdentRandPix = new QPushButton(frame);
+ pushIdentRandPix->setObjectName(QString::fromUtf8("pushIdentRandPix"));
+ pushIdentRandPix->setGeometry(QRect(140, 3, 101, 25));
+ checkRandPix = new QCheckBox(frame);
+ checkRandPix->setObjectName(QString::fromUtf8("checkRandPix"));
+ checkRandPix->setGeometry(QRect(5, 5, 21, 23));
+ frame_2 = new QFrame(centralwidget);
+ frame_2->setObjectName(QString::fromUtf8("frame_2"));
+ frame_2->setGeometry(QRect(198, 110, 251, 31));
+ frame_2->setFrameShape(QFrame::StyledPanel);
+ frame_2->setFrameShadow(QFrame::Raised);
+ label_3 = new QLabel(frame_2);
+ label_3->setObjectName(QString::fromUtf8("label_3"));
+ label_3->setGeometry(QRect(10, 6, 54, 16));
+ checkBank1 = new QCheckBox(frame_2);
+ checkBank1->setObjectName(QString::fromUtf8("checkBank1"));
+ checkBank1->setGeometry(QRect(65, 5, 41, 23));
+ checkBank2 = new QCheckBox(frame_2);
+ checkBank2->setObjectName(QString::fromUtf8("checkBank2"));
+ checkBank2->setGeometry(QRect(113, 5, 41, 23));
+ checkBank3 = new QCheckBox(frame_2);
+ checkBank3->setObjectName(QString::fromUtf8("checkBank3"));
+ checkBank3->setGeometry(QRect(161, 5, 41, 23));
+ checkBank4 = new QCheckBox(frame_2);
+ checkBank4->setObjectName(QString::fromUtf8("checkBank4"));
+ checkBank4->setGeometry(QRect(209, 5, 41, 23));
+ pushPlotSc = new QPushButton(centralwidget);
+ pushPlotSc->setObjectName(QString::fromUtf8("pushPlotSc"));
+ pushPlotSc->setGeometry(QRect(100, 150, 95, 51));
+ pushFitSc = new QPushButton(centralwidget);
+ pushFitSc->setObjectName(QString::fromUtf8("pushFitSc"));
+ pushFitSc->setGeometry(QRect(200, 150, 95, 51));
+ pushOpenBrowser = new QPushButton(centralwidget);
+ pushOpenBrowser->setObjectName(QString::fromUtf8("pushOpenBrowser"));
+ pushOpenBrowser->setGeometry(QRect(400, 150, 95, 51));
+ pushSave = new QPushButton(centralwidget);
+ pushSave->setObjectName(QString::fromUtf8("pushSave"));
+ pushSave->setGeometry(QRect(470, 15, 95, 51));
+ pushEmpty1 = new QPushButton(centralwidget);
+ pushEmpty1->setObjectName(QString::fromUtf8("pushEmpty1"));
+ pushEmpty1->setGeometry(QRect(500, 150, 95, 51));
+ frame_3 = new QFrame(centralwidget);
+ frame_3->setObjectName(QString::fromUtf8("frame_3"));
+ frame_3->setGeometry(QRect(97, 201, 201, 111));
+ frame_3->setFrameShape(QFrame::StyledPanel);
+ frame_3->setFrameShadow(QFrame::Raised);
+ spinPix = new QSpinBox(frame_3);
+ spinPix->setObjectName(QString::fromUtf8("spinPix"));
+ spinPix->setGeometry(QRect(114, 80, 81, 26));
+ spinPix->setAlignment(Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter);
+ spinPix->setMaximum(165887);
+ spinPix->setValue(0);
+ pushFitScSingle = new QPushButton(frame_3);
+ pushFitScSingle->setObjectName(QString::fromUtf8("pushFitScSingle"));
+ pushFitScSingle->setGeometry(QRect(103, 3, 95, 51));
+ pushPlotScSingle = new QPushButton(frame_3);
+ pushPlotScSingle->setObjectName(QString::fromUtf8("pushPlotScSingle"));
+ pushPlotScSingle->setGeometry(QRect(3, 3, 95, 51));
+ spinRow = new QSpinBox(frame_3);
+ spinRow->setObjectName(QString::fromUtf8("spinRow"));
+ spinRow->setGeometry(QRect(34, 56, 61, 26));
+ spinRow->setAlignment(Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter);
+ spinRow->setMaximum(575);
+ spinRow->setValue(0);
+ label_4 = new QLabel(frame_3);
+ label_4->setObjectName(QString::fromUtf8("label_4"));
+ label_4->setGeometry(QRect(5, 62, 31, 16));
+ spinCol = new QSpinBox(frame_3);
+ spinCol->setObjectName(QString::fromUtf8("spinCol"));
+ spinCol->setGeometry(QRect(34, 80, 61, 26));
+ spinCol->setAlignment(Qt::AlignRight|Qt::AlignTrailing|Qt::AlignVCenter);
+ spinCol->setMaximum(287);
+ spinCol->setValue(0);
+ label_5 = new QLabel(frame_3);
+ label_5->setObjectName(QString::fromUtf8("label_5"));
+ label_5->setGeometry(QRect(5, 86, 31, 16));
+ label_6 = new QLabel(frame_3);
+ label_6->setObjectName(QString::fromUtf8("label_6"));
+ label_6->setGeometry(QRect(119, 67, 31, 16));
+ pushClose = new QPushButton(centralwidget);
+ pushClose->setObjectName(QString::fromUtf8("pushClose"));
+ pushClose->setGeometry(QRect(470, 71, 95, 51));
+ frame_4 = new QFrame(centralwidget);
+ frame_4->setObjectName(QString::fromUtf8("frame_4"));
+ frame_4->setGeometry(QRect(59, 80, 135, 61));
+ frame_4->setFrameShape(QFrame::StyledPanel);
+ frame_4->setFrameShadow(QFrame::Raised);
+ radioHld = new QRadioButton(frame_4);
+ radioHld->setObjectName(QString::fromUtf8("radioHld"));
+ radioHld->setGeometry(QRect(10, 6, 51, 21));
+ radioHld->setLayoutDirection(Qt::RightToLeft);
+ radioRoot = new QRadioButton(frame_4);
+ radioRoot->setObjectName(QString::fromUtf8("radioRoot"));
+ radioRoot->setGeometry(QRect(62, 6, 51, 21));
+ pushFindNew = new QPushButton(frame_4);
+ pushFindNew->setObjectName(QString::fromUtf8("pushFindNew"));
+ pushFindNew->setGeometry(QRect(27, 26, 81, 31));
+ pushFindBack = new QPushButton(frame_4);
+ pushFindBack->setObjectName(QString::fromUtf8("pushFindBack"));
+ pushFindBack->setGeometry(QRect(3, 26, 21, 31));
+ pushFindForth = new QPushButton(frame_4);
+ pushFindForth->setObjectName(QString::fromUtf8("pushFindForth"));
+ pushFindForth->setGeometry(QRect(111, 26, 21, 31));
+ pushUnpack = new QPushButton(centralwidget);
+ pushUnpack->setObjectName(QString::fromUtf8("pushUnpack"));
+ pushUnpack->setGeometry(QRect(0, 150, 95, 51));
+ pushRewrite = new QPushButton(centralwidget);
+ pushRewrite->setObjectName(QString::fromUtf8("pushRewrite"));
+ pushRewrite->setGeometry(QRect(0, 205, 95, 51));
+ pushEmpty2 = new QPushButton(centralwidget);
+ pushEmpty2->setObjectName(QString::fromUtf8("pushEmpty2"));
+ pushEmpty2->setGeometry(QRect(500, 205, 95, 51));
+ pushSum = new QPushButton(centralwidget);
+ pushSum->setObjectName(QString::fromUtf8("pushSum"));
+ pushSum->setGeometry(QRect(300, 150, 95, 51));
+ pushFakes = new QPushButton(centralwidget);
+ pushFakes->setObjectName(QString::fromUtf8("pushFakes"));
+ pushFakes->setGeometry(QRect(300, 205, 95, 51));
+ MainWindow->setCentralWidget(centralwidget);
+ menubar = new QMenuBar(MainWindow);
+ menubar->setObjectName(QString::fromUtf8("menubar"));
+ menubar->setGeometry(QRect(0, 0, 595, 20));
+ MainWindow->setMenuBar(menubar);
+ statusbar = new QStatusBar(MainWindow);
+ statusbar->setObjectName(QString::fromUtf8("statusbar"));
+ MainWindow->setStatusBar(statusbar);
+
+ retranslateUi(MainWindow);
+
+ QMetaObject::connectSlotsByName(MainWindow);
+ } // setupUi
+
+ void retranslateUi(QMainWindow *MainWindow)
+ {
+ MainWindow->setWindowTitle(QApplication::translate("MainWindow", "MainWindow", 0, QApplication::UnicodeUTF8));
+ label->setText(QApplication::translate("MainWindow", "DIR", 0, QApplication::UnicodeUTF8));
+ label_2->setText(QApplication::translate("MainWindow", "FILE", 0, QApplication::UnicodeUTF8));
+ pushIdentRandPix->setText(QApplication::translate("MainWindow", "Rand. Pixels", 0, QApplication::UnicodeUTF8));
+ checkRandPix->setText(QApplication::translate("MainWindow", "CheckBox", 0, QApplication::UnicodeUTF8));
+ label_3->setText(QApplication::translate("MainWindow", "Bank(s)", 0, QApplication::UnicodeUTF8));
+ checkBank1->setText(QApplication::translate("MainWindow", "1", 0, QApplication::UnicodeUTF8));
+ checkBank2->setText(QApplication::translate("MainWindow", "2", 0, QApplication::UnicodeUTF8));
+ checkBank3->setText(QApplication::translate("MainWindow", "3", 0, QApplication::UnicodeUTF8));
+ checkBank4->setText(QApplication::translate("MainWindow", "4", 0, QApplication::UnicodeUTF8));
+ pushPlotSc->setText(QApplication::translate("MainWindow", "Plot \n"
+" S-Curves", 0, QApplication::UnicodeUTF8));
+ pushFitSc->setText(QApplication::translate("MainWindow", "Fit \n"
+" S-Curves", 0, QApplication::UnicodeUTF8));
+ pushOpenBrowser->setText(QApplication::translate("MainWindow", "Open \n"
+" Browser", 0, QApplication::UnicodeUTF8));
+ pushSave->setText(QApplication::translate("MainWindow", "Save", 0, QApplication::UnicodeUTF8));
+ pushEmpty1->setText(QApplication::translate("MainWindow", "Empty 1", 0, QApplication::UnicodeUTF8));
+ pushFitScSingle->setText(QApplication::translate("MainWindow", "Fit \n"
+" S-Curve", 0, QApplication::UnicodeUTF8));
+ pushPlotScSingle->setText(QApplication::translate("MainWindow", "Plot \n"
+" S-Curve", 0, QApplication::UnicodeUTF8));
+ label_4->setText(QApplication::translate("MainWindow", "Row", 0, QApplication::UnicodeUTF8));
+ label_5->setText(QApplication::translate("MainWindow", "Col", 0, QApplication::UnicodeUTF8));
+ label_6->setText(QApplication::translate("MainWindow", "Pix", 0, QApplication::UnicodeUTF8));
+ pushClose->setText(QApplication::translate("MainWindow", "Close \n"
+" windows", 0, QApplication::UnicodeUTF8));
+ radioHld->setText(QApplication::translate("MainWindow", "hld", 0, QApplication::UnicodeUTF8));
+ radioRoot->setText(QApplication::translate("MainWindow", "root", 0, QApplication::UnicodeUTF8));
+ pushFindNew->setText(QApplication::translate("MainWindow", "Find Newest", 0, QApplication::UnicodeUTF8));
+ pushFindBack->setText(QApplication::translate("MainWindow", "<", 0, QApplication::UnicodeUTF8));
+ pushFindForth->setText(QApplication::translate("MainWindow", ">", 0, QApplication::UnicodeUTF8));
+ pushUnpack->setText(QApplication::translate("MainWindow", "Unpack", 0, QApplication::UnicodeUTF8));
+ pushRewrite->setText(QApplication::translate("MainWindow", "Rewrite", 0, QApplication::UnicodeUTF8));
+ pushEmpty2->setText(QApplication::translate("MainWindow", "Empty 2", 0, QApplication::UnicodeUTF8));
+ pushSum->setText(QApplication::translate("MainWindow", "Sum", 0, QApplication::UnicodeUTF8));
+ pushFakes->setText(QApplication::translate("MainWindow", "RMS", 0, QApplication::UnicodeUTF8));
+ } // retranslateUi
+
+};
+
+namespace Ui {
+ class MainWindow: public Ui_MainWindow {};
+} // namespace Ui
+
+QT_END_NAMESPACE
+
+#endif // UI_SCURVES_H
--- /dev/null
+#ifndef CONFIG_H
+#define CONFIG_H
+
+#define SAVE 4
+#define ENDING "hld"
+#define CONTENTIN 1048576
+#define BUFSIZE 32000
+#define MAXSENSORS 64
+#define MAXSTATES 1140
+
+#endif
--- /dev/null
+#include "unpacker.C"
+
+void run()
+{
+// unpacker* a=new unpacker("/local/samir/data/normalmode/","60s_blank.hld");
+
+
+ unpacker* a=new unpacker("/local/samir/data","te15036113837");
+// unpacker* a=new unpacker("/local/samir/data","te15021111151");
+
+
+// unpacker* a=new unpacker("/local/samir/data","te15021111151.hld");
+// unpacker* a=new unpacker("/local/samir/data","test");
+// unpacker* a=new unpacker("/local/samir/data","test.hld");
+ delete a;
+}
\ No newline at end of file
--- /dev/null
+#include "config.h"
+#include "unpack.C"
+
+
+//####################################################################
+
+void rungui()
+{
+ unpacker* a=new unpacker("/local/samir/data/te15021111151.hld");
+ a->scanFrames(0,10000000);
+ delete a;
+
+}
\ No newline at end of file
--- /dev/null
+#ifndef UNPACKER_H
+#define UNPACKER_H
+
+#define CONTENTIN 1048576 // --> 1 MByte
+
+//-----------------
+#include "unpacker.h"
+
+#include <dirent.h>
+#include <help.h>
+// #include <time.h>
+// #include <sys/stat.h>
+//-----------------
+
+//####################################################################
+
+unpacker::unpacker(TString Dir, TString File) {
+
+ fPixelProb = new Float_t [1152*576];
+ fRowCount = new UInt_t [576];
+ for(int i=0;i<1152*576;i++) { fPixelProb [i]=0; }
+ for(int i=0;i<576;i++) { fRowCount [i]=0; }
+
+// fWrongRowNumber = 0;
+ fWrongRunNumber = 0;
+ fWrongThreshold = 0;
+ fHeaderVer = 0;
+
+ fRunNumber_tmp = 0;
+// fRow_tmp = 0;
+ fThreshold_tmp = 0;
+ fRunNumber_count= 0;
+ fRow_count = 0;
+
+ fDataTree = 0;
+ fOutFile = 0;
+ fFirstOut = 0;
+
+// Init
+ if( !openFiles(Dir, File) ) { cerr << "ERROR: No Files found!" << endl << "-- Exit --" << endl; exit(-1); }
+ scanFiles();
+
+}
+
+//####################################################################
+
+unpacker::unpacker() {
+}
+
+//####################################################################
+
+unpacker::~unpacker() {
+
+ if( fDataTree!=0 )
+ {
+ fOutFile->Write();
+ fOutFile->Save();
+ fOutFile->Close();
+ }
+
+ if( fNrFiles>0 )
+ {
+ for(uint i=0; i<fNrFiles;i++) { fInn[i].close(); }
+ delete[] fFileLength;
+ delete[] fInn;
+ delete[] fRawData;
+ }
+
+ delete[] fPixelProb;
+ delete[] fRowCount;
+}
+
+//####################################################################
+
+void unpacker::scanFiles() {
+
+ int eventNr = 0;
+ int subEventNr = 0;
+ int subSubEventNr = 0;
+
+ UInt_t HDRver;
+ UInt_t DATAver = 0;
+// UInt_t HDRsize;
+
+// Progress bar
+// ULong_t progress_tmp=-1;
+// Float_t progress;
+
+ while( getNextEvent() !=0 ) { eventNr++; subEventNr = 0;
+ while( getNextSubEvent() !=0 ) { subEventNr++; subSubEventNr = 0;
+ while( getNextSubSubEvent() !=0 ) { subSubEventNr++;
+// Progress bar
+// progress = (Int_t)(((fPosGlo)*100.)/(fLength-1)*1);
+// if(progress!=progress_tmp) { print_progress( (((fPosGlo)*100.)/(fLength-1)) ); progress_tmp=progress; }
+// --------------------
+// fRawData : allocates the event data
+// fPosInData : points to the position of the subsubevent in the data, beginning with the subsubevent hld-header
+// fSubSubEventSize : denotes the number of 8 bits words in the subsubevent including the hld-header
+// --------------------
+ if(fSubSubEventSize>4)
+ {
+ HDRver = bigEndian8(fRawData,fPosInData+4);
+ DATAver = bigEndian8(fRawData,fPosInData+5);
+// HDRsize = bigEndian8(fRawData,fPosInData+7);
+
+ if ( DATAver==0x01)
+ {
+
+ }
+ else if ( DATAver==0x02 )
+ {
+
+ }
+ else if ( DATAver==0xC0 )
+ {
+ if( !fHeaderVer && fSubSubEventSize==164 ) { fHeaderVer=1; cout<<"\r 'Frame Timestamp' missing in Header! Sub-Sub-Event Length: "<<fSubSubEventSize/4<<" "<<endl; }
+ anaScurves();
+ }
+ }
+
+// }}if(eventNr==600){exit(-1);}}
+ }}}
+
+ if( DATAver==0xC0 )
+ {
+ saveScurveData();
+ }
+
+ cout<<"----------"<<endl;
+ cout<<"Fertig!"<<endl;
+ cout<<"----------"<<endl;
+}
+
+//####################################################################
+
+void unpacker::anaScurves() {
+
+ UInt_t Threshold = bigEndian8 (fRawData,fPosInData+16);
+ UInt_t Bank = (fRawData[fPosInData+17]>>4) & 0xF;
+ UInt_t Row = ((Int_t)(UChar_t)fRawData[fPosInData+17]& 0xF)*TMath::Power(2,8) + (Int_t)(UChar_t)fRawData[fPosInData+18];
+ UInt_t RunNumber = bigEndian8 (fRawData,fPosInData+19);
+ UInt_t Data;
+
+// cout<<Threshold<<"\t"<<Bank<<"\t"<<Row<<"\t"<<RunNumber<<"\t"<<endl;
+// -------------------------------
+// Check: hardware row counting
+ if(Row>574)
+ {
+ if( !fWrongRowNumber )
+ {
+ fWrongRowNumber=1; cout<<"\r Wrong Row Number(s) found! "<<endl;
+ }
+
+ Row = Row%575;
+ }
+// -------------------------------
+// Save data
+ if( fThreshold_tmp!=Threshold && fFirstOut==1 )
+ {
+ saveScurveData();
+
+ for(int row=0;row<576;row++)
+ {
+ for(int col=0;col<1152;col++)
+ {
+ fPixelProb[1152*row+col] = 0;
+ }
+ fRowCount[row] = 0;
+ }
+ }
+
+ fFirstOut = 1;
+// -------------------------------
+// Check: whether next threshold starts after one whole frame
+ if( fRow_count!=0 && fThreshold_tmp!=Threshold )
+ {
+ if( !fWrongThreshold )
+ {
+ fWrongThreshold=1; cout<<"\r Wrong Threshold changing(s) found! "<<endl;
+ }
+ }
+// -------------------------------
+// Fill: fPixelProb array
+ for(UInt_t dataword=0; dataword<36; dataword++)
+ {
+ Data = bigEndian32 (fRawData,fPosInData+20+dataword);
+
+ if (Data == 0)
+ {
+ }
+ else if (Data == 0xFFFFFFFF)
+ {
+ for(int pixel=0;pixel<32;pixel++)
+ {
+ fPixelProb[1152*(fRow_count+1)+32*dataword+pixel ] ++;
+ }
+ }
+ else
+ {
+ for(int pixel=0;pixel<32;pixel++)
+ {
+ fPixelProb[1152*(fRow_count+1)+32*dataword+pixel ] += (Data & 0x1);
+ Data >>= 1;
+ }
+ }
+ }
+
+ fRowCount[fRow_count]++;
+
+ fRunNumber_tmp = RunNumber;
+// fRow_tmp = Row;
+ fThreshold_tmp = Threshold;
+ if( fRow_count==0 ) { fRunNumber_count++; }
+
+ fRow_count++;
+ fRow_count = fRow_count%575;
+
+// cout<<setw(10)<<right<< Bank;
+// cout<<setw(10)<<right<< Threshold;
+// cout<<setw(10)<<right<< RunNumber;
+// cout<<setw(10)<<right<< Row;
+// cout<<endl;
+}
+
+//####################################################################
+
+void unpacker::saveScurveData() {
+
+ if( fDataTree==0 )
+ {
+ cout<<"\r Tree started: "<<fOutFileName<<" "<<endl;
+ fOutFile = new TFile(fOutFileName,"RECREATE");
+ fDataTree = new TTree("scurvesall", "scurvesall");
+ fDataTree->Branch("threshold" , &fThreshold_tmp , "threshold/s" , 32000);
+ fDataTree->Branch("pixelprob" , fPixelProb , "pixelprob[663552]/F" , 32000);
+ fDataTree->Branch("frames" , fRowCount , "frames[576]/s" , 32000);
+ }
+
+ for(int row=0;row<576;row++)
+ {
+ for(int col=0;col<1152;col++)
+ {
+ if( fRowCount[row]!=0 && fPixelProb[1152*row+col]!=0 )
+ {
+// cout<<setw(25)<<setprecision(5)<<right<< fPixelProb[1152*row+col];
+ fPixelProb[1152*row+col] = 1.*fPixelProb[1152*row+col]/fRowCount[row];
+// cout<<setw(25)<<setprecision(5)<<right<< fPixelProb[1152*row+col];
+// cout<<setw(25)<<setprecision(5)<<right<< fRowCount[row];
+// cout<<endl;
+ }
+ }
+ }
+
+ fDataTree->Fill();
+}
+
+//####################################################################
+//####################################################################
+
+bool unpacker::openFiles(TString Dir, TString File) {
+
+ // Check whether files are there
+ // Open all files
+ // Check files: Size
+
+ fNrFiles = 0;
+
+ char *dir = (char*) Dir.Data();
+ char *file = (char*) File.Data();
+ char ending[] = ".hld";
+ bool has_ending = true;
+ bool has_ex = false;
+ uint len;
+
+ ifstream inn;
+ TString FileName;
+ TString FileNames[100];
+
+ ULong_t End;
+
+// Check naming of 'Dir' and 'File'
+ if( dir[strlen(dir)-1] != '/' )
+ {
+ strcat(dir, "/" );
+ Dir = dir;
+ }
+
+ len = strlen(file);
+
+ for(uint i=0;i<4;i++)
+ {
+ if( file [len-4+i] != ending[i] )
+ {
+ has_ending = false;
+ }
+ }
+
+ for(uint i=0;i<len;i++)
+ {
+ if( file [i] == '.' )
+ {
+ has_ex = true;
+ }
+ }
+
+ if(has_ex && !has_ending)
+ {
+ cerr<<"WARNING: Try to read file(s) with unknown filename extension: "<<file<<endl;
+ }
+
+ if( has_ending )
+ {
+ file[len-4] = '\0';
+ File = file;
+ }
+
+// Try to find 'Dir' and 'File'
+ DIR* dpdf = opendir(Dir);
+ if (dpdf == NULL)
+ {
+ cerr<<"ERROR: Directory not found: "<<Dir<<endl;
+ return 0;
+ }
+
+ do
+ {
+ FileName = Dir + File + Form("_%i",fNrFiles) + (TString)ending;
+ inn.open(FileName, ios::binary);
+
+ if(!inn.good())
+ {
+ inn.close();
+ break;
+ }
+
+ inn.close();
+ FileNames[fNrFiles] = FileName;
+ fNrFiles++;
+
+ }while(inn.good());
+
+ if( fNrFiles==0 )
+ {
+ FileName = Dir + File + (TString)ending;
+ inn.open(FileName, ios::binary);
+
+ if(!inn.good())
+ {
+ inn.close();
+ return 0;
+ }
+
+ inn.close();
+ FileNames[fNrFiles] = FileName;
+ fNrFiles++;
+ }
+
+// Link File(s)
+ if( fNrFiles>0 )
+ {
+ len = strlen(FileNames[fNrFiles-1]);
+
+ fRawData = new char[CONTENTIN];
+ fInn = new ifstream [fNrFiles];
+ fFileLength = new ULong_t [fNrFiles];
+ fLength = 0;
+
+ cout<<"--------"<<endl;
+ cout<<"Linked file(s):"<<endl;
+
+ for(uint i=0; i<fNrFiles;i++)
+ {
+ fInn[i].open(FileNames[i], ios::binary);
+ fInn[i].seekg(0, ios::end);
+ End = fInn[i].tellg();
+ fFileLength[i] = End;
+ fLength+=End;
+ fInn[i].seekg(0, ios::beg);
+
+
+ cout<<" "<<setw(len+1)<<left<< FileNames[i];
+
+ if (End/(TMath::Power(2,40)) >= 1) { printf(" - %5.2f TBytes", (Float_t)(End/(TMath::Power(2,40))) ); }
+ else if (End/(TMath::Power(2,30)) >= 1) { printf(" - %5.2f GBytes", (Float_t)(End/(TMath::Power(2,30))) ); }
+ else if (End/(TMath::Power(2,20)) >= 1) { printf(" - %5.2f MBytes", (Float_t)(End/(TMath::Power(2,20))) ); }
+ else if (End/(TMath::Power(2,10)) >= 1) { printf(" - %5.2f kBytes", (Float_t)(End/(TMath::Power(2,10))) ); }
+ else { printf(" - %5.2f Bytes ", (Float_t)(End)); }
+
+ cout<<endl;
+ }
+ cout<<"--------"<<endl;
+ }
+
+// ----------------
+ fOutFileName = Dir + File + Form("_DST.root");
+
+ fCurrFile = 0;
+ fCurrPosInFile = 0;
+ fPosGlo = 0;
+// ----------------
+
+ return fNrFiles;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextEvent() {
+
+ do
+ {
+ fEventSize = getNextPseudoEvent();
+ fPosGlo += fEventSize;
+
+ }while( fEventSize==32 );
+
+ fPosInEvent = 32;
+ fSubEventSize = 0;
+
+ return fEventSize;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextPseudoEvent() {
+
+ fCurrPosInFile = fInn[fCurrFile].tellg();
+
+ if( fCurrPosInFile >= fFileLength[fCurrFile] )
+ {
+ if( fCurrFile+1<fNrFiles )
+ {
+ fCurrFile++;
+ fCurrPosInFile = 0;
+ }
+ else
+ {
+ return 0;
+ }
+ }
+
+ fInn[fCurrFile].read (fRawData, 4);
+ fEventSize = littleEndian32(fRawData,0);
+
+ UInt_t padding;
+ if (fEventSize%8!=0) { padding = 8 - fEventSize%8; }
+ else { padding = 0; }
+
+ UInt_t eventRead = fEventSize + padding;
+ fInn[fCurrFile].seekg (fCurrPosInFile, ios::beg);
+ fInn[fCurrFile].read (fRawData, eventRead);
+
+ return fEventSize;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextSubEvent() {
+
+ fPosInEvent += fSubEventSize;
+
+ if( fPosInEvent>=fEventSize ) {return 0;} // There are no more Sub-Events
+
+ fSubEventSize = bigEndian32(fRawData,fPosInEvent); // including itself
+// fPosInEvent += 16;
+
+ if (fSubEventSize%8!=0) { fPaddingSub = 8 - fSubEventSize%8;}
+ else { fPaddingSub = 0; }
+ fSubEventSize +=fPaddingSub;
+ fSubSubEventSize = 0;
+ fPosInData = fPosInEvent+16;
+
+ return fSubEventSize;
+}
+
+//####################################################################
+
+UInt_t unpacker::getNextSubSubEvent() {
+
+ fPosInData += fSubSubEventSize;
+
+ if( fPosInData+fPaddingSub>=fPosInEvent+fSubEventSize-8 ) { return 0; }
+
+ fSubSubEventSize = 4*bigEndian16(fRawData,fPosInData) + 4; // Number of 8 bit words in subsubevent including itself
+
+ return fSubSubEventSize;
+}
+
+//####################################################################
+
+UInt_t unpacker::littleEndian32(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[3+Pos]*TMath::Power(2,24) +
+ (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,16) +
+ (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[0+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian32(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,24) +
+ (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,16) +
+ (Int_t)(UChar_t)RawData[2+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[3+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian24(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,16) +
+ (Int_t)(UChar_t)RawData[1+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[2+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian16(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos]*TMath::Power(2,8) +
+ (Int_t)(UChar_t)RawData[1+Pos];
+}
+
+//####################################################################
+
+UInt_t unpacker::bigEndian8(char* RawData, Int_t Pos) {
+
+ return
+ (Int_t)(UChar_t)RawData[0+Pos];
+}
+
+//####################################################################
+#endif
\ No newline at end of file
--- /dev/null
+#include <header.h>
+
+class unpacker {
+
+private:
+// File info
+ UInt_t fNrFiles; // Number of files found (Unsigned integer 4 bytes)
+ ULong_t* fFileLength; // Length of each respective file (Unsigned long integer 8 bytes)
+ ULong_t fLength; // Length integrated
+// File access
+ ifstream* fInn; // Link to each file
+ char* fRawData; // in stream buffer
+
+// Output file:
+// TFile* fOutputFile;
+ TString fOutFileName;
+
+// Coordinates in file(s)
+ UInt_t fCurrFile; // Current file accessed
+ ULong_t fCurrPosInFile; // Current position in accessed file
+
+ UInt_t fEventSize;
+ UInt_t fSubEventSize;
+ UInt_t fSubSubEventSize;
+ UInt_t fPosGlo;
+ UInt_t fPosInEvent;
+ UInt_t fPosInData;
+ UInt_t fPaddingSub;
+
+ Float_t* fPixelProb;
+ UInt_t* fRowCount;
+ UInt_t fRunNumber_count;
+ UInt_t fRunNumber_tmp;
+ UInt_t fRow_count;
+ UInt_t fRow_tmp;
+ UInt_t fThreshold_tmp;
+
+ bool fWrongRowNumber;
+ bool fWrongThreshold;
+ bool fWrongRunNumber;
+ bool fFirstOut;
+ bool fHeaderVer;
+
+ TTree* fDataTree;
+ TFile* fOutFile;
+
+// Methods
+ bool openFiles (TString Dir, TString File);
+ UInt_t getNextEvent ();
+ UInt_t getNextPseudoEvent ();
+ UInt_t getNextSubEvent ();
+ UInt_t getNextSubSubEvent ();
+
+ UInt_t littleEndian32 (char* RawData, Int_t Pos);
+ UInt_t bigEndian32 (char* RawData, Int_t Pos);
+ UInt_t bigEndian24 (char* RawData, Int_t Pos);
+ UInt_t bigEndian16 (char* RawData, Int_t Pos);
+ UInt_t bigEndian8 (char* RawData, Int_t Pos);
+
+public:
+// Methods
+ unpacker (TString Dir, TString File);
+ unpacker ();
+ ~unpacker ();
+ void scanFiles ();
+ void anaScurves ();
+ void saveScurveData ();
+
+};
\ No newline at end of file
jspc29 / mvd_soft.git / commitdiff