Subversion Repositories f9daq

Compare Revisions

Ignore whitespace Rev 90 → Rev 89

/sipmScan/src/analysisScan.cpp
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 0 // Right edge of SiPM+Lightguide
#define OFFSETY 0 // Lower edge of SiPM+Lightguide
#define OFFSETX 5220 // Right edge of SiPM+Lightguide
#define OFFSETY 5860 // Lower edge of SiPM+Lightguide
 
#define RUNREC_ID 1
#define ENDREC_ID 2
84,18 → 84,14
 
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", double c_tdcOffset=97, double tdcCut=5.0, int dbg=0)
int analysisScan(char* dfile0="test", int dbg=0, double tdcCut=5.0)
{
 
//const double c_tdcOffset = +98*TDC_BIN; // ns
c_tdcOffset *= TDC_BIN; //ns
const double c_tdcOffset = +97*TDC_BIN; // ns, SiPM4
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);
160,7 → 156,7
char hname[256];
//double tdc;
TH2F *htdc;
TH2F* h_correctedTDC;
TH2F* h_correctedTDC;
TH1F *hnhitsx[NCH], *hnhitsy[NCH];
TH2F *h2d[NCH];
170,9 → 166,12
unsigned int readbuf[READBUFFERLENGTH];
unsigned int buf[READBUFFERLENGTH];
//data file
//gzFile dfp;
FILE *dfp;
gzFile dfp;
char dfile[256];
int ftype=0;
int fcount=1;
191,8 → 190,7
printf(" Cannot find data file for %s !!!\n", dfile0);
return -1;
}
//dfp=gzopen(dfile,"rb");
dfp=fopen(dfile,"rb");
dfp=gzopen(dfile,"rb");
} while(!dfp);
printf("Opened data file %s\n", dfile);
221,11 → 219,9
int status;
while(1) {
//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);
if(gzeof(dfp)) end_of_file = 1;
 
gzread(dfp, (voidp)&readbuf, 2*ulsize);
rec_id=readbuf[0];
rec_len=readbuf[1];
236,8 → 232,7
switch(rec_id)
{
case RUNREC_ID:
//gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp);
gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
runrec = (RUNREC*) readbuf;
run = *runrec;
303,8 → 298,7
break;
case POSREC_ID:
//gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp);
gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
posrec = (POSREC *) readbuf;
pos=*posrec;
318,8 → 312,7
break;
case EVTREC_ID:
//gzread(dfp, (voidp)&readbuf[2], ulsize); // last field of event record
fread((void*)&readbuf[2], ulsize, 1, dfp);
gzread(dfp, (voidp)&readbuf[2], ulsize); // last field of event record
evtrec = (EVTREC *) readbuf;
//evtrec->nev = buf[0];
//if (rec_len < 0 || rec_len > 10000) {
328,8 → 321,7
return(0);
}
nb = rec_len - 3*ulsize; // no. of bytes to read
//gzread(dfp, (voidp)&buf, nb);
fread((void*)&buf, nb, 1, dfp);
gzread(dfp, (voidp)&buf, nb);
if(dbg) {
printf("EVTREC_ID\n");
350,10 → 342,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 ;
431,8 → 423,7
break;
case THRREC_ID:
//status = gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
status = fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp);
status = gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
thrrec = (THRREC*) readbuf;
thr = *thrrec;
if (dbg) printf("THRREC id = %d len = %d threshold %d\n",
440,8 → 431,7
break;
case ENDREC_ID:
//gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
fread((void*)&readbuf[2], (rec_len-2*ulsize), 1, dfp);
gzread(dfp, (voidp)&readbuf[2], (rec_len-2*ulsize));
endrec = (ENDREC *) readbuf;
if(dbg) {
462,11 → 452,9
break;
}
//if(dfp) gzclose(dfp);
if(dfp) fclose(dfp);
if(dfp) gzclose(dfp);
//dfp=gzopen(dfile,"rb");
dfp=fopen(dfile,"rb");
dfp=gzopen(dfile,"rb");
if(!dfp) {
printf(" Cannot open data file: %s ---> Exiting\n", dfile);
end_of_file = 1;
490,9 → 478,8
}
 
if(dfp) {
//gzclose(dfp);
fclose(dfp);
//delete dfp;
gzclose(dfp);
delete dfp;
}
if(dbg) return 1;
if(rootfile) {
500,7 → 487,7
rootfile->Write();
printf("Saved to %s\n", fnameroot);
rootfile->Close();
//delete rootfile;
delete rootfile;
}
return 1;