WebSVN – f9daq – Diff – /sipmScan/src/analysisScan.cpp

 #include "TF1.h"
 #include "TMath.h"
 #include "TStyle.h"
 #include "TCanvas.h"
 #include "TLine.h"
-#include "zlib.h"
+//#include "zlib.h"
 // ------------------------------------------------------------------------------
 #define POSMARG 1000
 // data format
 #define MAXDATA 16
 #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 OFFSETX 0  // Right edge of SiPM+Lightguide
-#define OFFSETY 5860 // Lower edge of SiPM+Lightguide
+#define OFFSETY 0 // Lower edge of SiPM+Lightguide
 #define RUNREC_ID 1
 #define ENDREC_ID 2
 #define POSREC_ID 3
 #define EVTREC_ID 4
 THRREC *thrrec;
 THRREC thr;
 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);
   char fullname[256];
   char sbuff[256];
 /*    }
     fclose(fp);
   }*/
   //for(int i=0; i<NCH; i++) printf("%.2lf %.2lf %.2lf\n", corpar[i][0], corpar[i][1], corpar[i][2]);
         //histograms
         char hname[256];
         //double tdc;
         TH2F *htdc;
-  TH2F* h_correctedTDC;
+        TH2F* h_correctedTDC;
         TH1F *hnhitsx[NCH], *hnhitsy[NCH];
         TH2F *h2d[NCH];
         TNtuple *nt;
         //data buffer
         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;
         do {
           switch(ftype++) {
               break;
             default:
         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);
         //opens ROOT file
         //TFile *rootfile;
         int nint;
         int nb;
         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];
                 if(dbg) printf("-----------------------------------------------\n");
                 if(dbg) printf("[%d] rec_id = %d | rec_len = %u\n", ceve, rec_id, rec_len);
                 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;
                                 if(dbg) {
                                         printf("RUNREC_ID\n");
+                                }
                                 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;
                                 if(dbg) {
                                         printf("POSREC_ID\n");
                                 } else  printf("        [%d,%d] %d, %d\n", pos.ix, pos.iy, pos.xset, pos.yset);
                                 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) {
                           if (rec_len > READBUFFERLENGTH) {
                             printf("Len %u\n", rec_len);
                             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");
                                         printf("id = %d, len = %d, nev = %d\n", evtrec->id, evtrec->len, evtrec->nev);
                                         //for(int datai = 0; datai < NCH; datai++) printf("%u ", evtrec->data[datai]);
           int recid = buf[ii++];
           int len   = buf[ii++];
           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  ;
               int overflow  = (dbuf[i]>>17)&0x1  ;
               int tdc_num   = (dbuf[i]>>25)&0x1  ;
+          }
         } */
                                 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",
                                     thrrec->id, thrrec->len, thrrec->threshold);
                     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) {
                                         printf("ENDREC_ID\n");
                                         printf("id = %d, len = %d, time = %d\n", endrec->id, endrec->len, endrec->time);
             case 2:
               sprintf(dfile, "./data/%s_file%02d.gz", dfile0, fcount);
               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;
         } else {
           printf("      Opened data file: %s\n", dfile);
                 //if(dbg) if( ceve>dbg ) break;
                 if(end_of_file) break;
+        }
         if(dfp) {
-          gzclose(dfp);
+          //gzclose(dfp);
+          fclose(dfp);
-          delete dfp;
+          //delete dfp;
+        }
         if(dbg) return 1;
         if(rootfile) {
           nt->Write();
                 rootfile->Write();
                 printf("Saved to %s\n", fnameroot);
                 rootfile->Close();
-                delete rootfile;
+                //delete rootfile;
+        }
         return 1;
+}

Subversion Repositories f9daq

(root)/sipmScan/src/analysisScan.cpp – Rev 50 → 90