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); |
|
156,7 → 160,7 |
char hname[256]; |
//double tdc; |
TH2F *htdc; |
TH2F* h_correctedTDC; |
TH2F* h_correctedTDC; |
TH1F *hnhitsx[NCH], *hnhitsy[NCH]; |
TH2F *h2d[NCH]; |
|
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; |
|
gzread(dfp, (voidp)&readbuf, 2*ulsize); |
//if(gzeof(dfp)) end_of_file = 1; |
if(feof(dfp)) end_of_file = 1; |
|
//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"); |
342,10 → 350,10 |
if (dbg) printf("Buffer pointer %d\n", ii); |
unsigned int *dbuf = (unsigned int *)&buf[ii]; |
//if (n%1000==0) |
if (dbg) printf("%d 0x%03x Len=%d\n",evtrec->nev,recid,len); |
if (dbg) printf("%d 0x%03x Len=%d\n",evtrec->nev,recid,len); |
//unsigned short edge; |
//int nhits; |
if (recid==0x140 || recid==0x141) { |
if (recid==0x140 || recid==0x141) { |
for (int i=0; i<len; i++) { |
int data = dbuf[i] & 0xFFFF ; |
int edge_type = (dbuf[i]>>16)&0x1 ; |
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; |