| 17,7 → 17,7 |
|---|
| #include "TStyle.h" |
| #include "TCanvas.h" |
| #include "TLine.h" |
| #include "zlib.h" |
| //#include "zlib.h" |
| |
| // ------------------------------------------------------------------------------ |
| |
| 30,8 → 30,8 |
|---|
| #define NCH 64 |
| #define TDC_BIN 1.0416 // 1 TDC bin in ns |
| #define MIKRO_BIN 0.49609/1000. //1 mikro step in mm; stage MM3MF |
| #define OFFSETX 5220 // Right edge of SiPM+Lightguide |
| #define OFFSETY 5860 // Lower edge of SiPM+Lightguide |
| #define OFFSETX 0 // Right edge of SiPM+Lightguide |
| #define OFFSETY 0 // Lower edge of SiPM+Lightguide |
| |
| #define RUNREC_ID 1 |
| #define ENDREC_ID 2 |
| 84,14 → 84,18 |
|---|
| |
| double padCenter[NCH][2]; |
| |
| // Test if the hit(x,y) lies inside the coordinates of the channel |
| // file with pad centers (without offset) should be provided |
| int position(int, int, int); |
| |
| // ------------------------------------------------------------------------------ |
| |
| int analysisScan(char* dfile0="test", int dbg=0, double tdcCut=5.0) |
| int analysisScan(char* dfile0="test", double c_tdcOffset=97, double tdcCut=5.0, int dbg=0) |
| { |
| |
| const double c_tdcOffset = +97*TDC_BIN; // ns, SiPM4 |
| //const double c_tdcOffset = +98*TDC_BIN; // ns |
| c_tdcOffset *= TDC_BIN; //ns |
| |
| printf(" Data to root conversion program\nUsage:\nd2r(input file name <without .dat>, debug on/off, TDC cut +-[ns])\n\n"); |
| printf(" OFFSETS: \n x: %d\n y: %d\n t:%f\n", OFFSETX, OFFSETY, c_tdcOffset); |
| |
| 166,12 → 170,9 |
|---|
| unsigned int readbuf[READBUFFERLENGTH]; |
| unsigned int buf[READBUFFERLENGTH]; |
| |
| |
| |
| |
| |
| //data file |
| gzFile dfp; |
| //gzFile dfp; |
| FILE *dfp; |
| char dfile[256]; |
| int ftype=0; |
| int fcount=1; |
| 190,7 → 191,8 |
|---|
| printf(" Cannot find data file for %s !!!\n", dfile0); |
| return -1; |
| } |
| dfp=gzopen(dfile,"rb"); |
| //dfp=gzopen(dfile,"rb"); |
| dfp=fopen(dfile,"rb"); |
| } while(!dfp); |
| printf("Opened data file %s\n", dfile); |
| |
| 219,9 → 221,11 |
|---|
| int status; |
| |
| while(1) { |
| if(gzeof(dfp)) end_of_file = 1; |
| //if(gzeof(dfp)) end_of_file = 1; |
| if(feof(dfp)) end_of_file = 1; |
| |
| gzread(dfp, (voidp)&readbuf, 2*ulsize); |
| //gzread(dfp, (voidp)&readbuf, 2*ulsize); |
| fread((void*)&readbuf, 2*ulsize, 1, dfp); |
| rec_id=readbuf[0]; |
| rec_len=readbuf[1]; |
| |
| 232,7 → 236,8 |
|---|
| switch(rec_id) |
| { |
| case RUNREC_ID: |
| gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| //gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp); |
| runrec = (RUNREC*) readbuf; |
| run = *runrec; |
| |
| 298,7 → 303,8 |
|---|
| break; |
| |
| case POSREC_ID: |
| gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| //gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp); |
| posrec = (POSREC *) readbuf; |
| pos=*posrec; |
| |
| 312,7 → 318,8 |
|---|
| break; |
| |
| case EVTREC_ID: |
| gzread(dfp, (voidp)&readbuf[2], ulsize); // last field of event record |
| //gzread(dfp, (voidp)&readbuf[2], ulsize); // last field of event record |
| fread((void*)&readbuf[2], ulsize, 1, dfp); |
| evtrec = (EVTREC *) readbuf; |
| //evtrec->nev = buf[0]; |
| //if (rec_len < 0 || rec_len > 10000) { |
| 321,7 → 328,8 |
|---|
| return(0); |
| } |
| nb = rec_len - 3*ulsize; // no. of bytes to read |
| gzread(dfp, (voidp)&buf, nb); |
| //gzread(dfp, (voidp)&buf, nb); |
| fread((void*)&buf, nb, 1, dfp); |
| |
| if(dbg) { |
| printf("EVTREC_ID\n"); |
| 423,7 → 431,8 |
|---|
| break; |
| |
| case THRREC_ID: |
| status = gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| //status = gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| status = fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp); |
| thrrec = (THRREC*) readbuf; |
| thr = *thrrec; |
| if (dbg) printf("THRREC id = %d len = %d threshold %d\n", |
| 431,7 → 440,8 |
|---|
| break; |
| |
| case ENDREC_ID: |
| gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| //gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize)); |
| fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp); |
| endrec = (ENDREC *) readbuf; |
| |
| if(dbg) { |
| 452,9 → 462,11 |
|---|
| break; |
| } |
| |
| if(dfp) gzclose(dfp); |
| //if(dfp) gzclose(dfp); |
| if(dfp) fclose(dfp); |
| |
| dfp=gzopen(dfile,"rb"); |
| //dfp=gzopen(dfile,"rb"); |
| dfp=fopen(dfile,"rb"); |
| if(!dfp) { |
| printf(" Cannot open data file: %s ---> Exiting\n", dfile); |
| end_of_file = 1; |
| 478,8 → 490,9 |
|---|
| } |
| |
| if(dfp) { |
| gzclose(dfp); |
| delete dfp; |
| //gzclose(dfp); |
| fclose(dfp); |
| //delete dfp; |
| } |
| if(dbg) return 1; |
| if(rootfile) { |
| 487,7 → 500,7 |
|---|
| rootfile->Write(); |
| printf("Saved to %s\n", fnameroot); |
| rootfile->Close(); |
| delete rootfile; |
| //delete rootfile; |
| } |
| |
| return 1; |