WebSVN – f9daq – Diff – /sipmScan/src/sipmScan.c

 #include <userint.h>
 #include <ansi_c.h>
 #include <utility.h>
 #include <analysis.h>
+#include <time.h>
+//#include <zlib.h>
 #include "c:\HOME\dino\sipmScan\include\sipmScan_ui.h"
 #ifdef USE_DAQ
 #  define USE_CAMAC
 #  include "camac.h"
 #  define NDAC 1
 #    define ASD8 0 // channels
 #    define AMPDISSH 1
-#    define SIPM 2
+#    define SIPM 15
 //#  define CAEN_V673A   0x10110000 // FMF
 #  define CAEN_V673A   0x22220000 // IJS
 #  define CAEN_V462    0x100300
 #endif
 #ifdef USE_MIKRO
-#  define MIKRO_COM 5
+#  define MIKRO_COM 3
 #  define MIKRO_X 1
 #  define USE_MIKRO_Y
 #  define MIKRO_Y 2
 #  define USE_MIKRO_Z
 #  define MIKRO_Z 3
 #  define STEP_TOLERANCE 1
 #endif
-#define MAXCH 512
+#define MAXCH 128
 #define MAX_THREADS 10
 #define IWAIT 200
 #define POWERSUPPLY 66 // Power supply voltage for SiPM bias
 static int something[NCH][2][MAXCH];
 static int daq_on;
 static int poolHandle = 0;
 static int ntics,dummy;
+static int offsetX = 5600;
+static int offsetY = 5400;
 /************************************************************/
 void wait_loop(unsigned long iloop)
+{
         int maxn = BSIZE-4; // 2 words x 2 headers
         int tdcdata, edge_type, overflow, tdc_num, channel, ev_dat, last_dat, nval_dat;
         // end of declarations
         // ***
         runrec.id = RUNREC_ID;
         runrec.len = sizeof(runrec);
         //runrec.fver = 0x10000;
         runrec.nch = NCH;
         GetCtrlVal(p1h, P1_ADCHLSAVE, &runrec.xy);//runrec.xy = 1;
         endrec.len = sizeof(endrec);
         posrec.id = POSREC_ID;
         posrec.len = sizeof(posrec);
         evtrec.id = EVTREC_ID;
         evtrec.len = sizeof(evtrec);
         cres = 0;
         GetCtrlVal (p1h, P1_NEVE, &runrec.nev);
         GetCtrlVal (p1h, P1_PEDESTAL, &runrec.ped);
         GetCtrlVal (p1h, P1_NX, &runrec.nx);
         GetCtrlVal (p1h, P1_XSTEP, &runrec.dx);
         GetCtrlVal (p1h, P1_XMIN, &runrec.x0);
         GetCtrlVal (p1h, P1_NY, &runrec.ny);
         GetCtrlVal (p1h, P1_YSTEP, &runrec.dy);
         GetCtrlVal (p1h, P1_NMIN, &step_minutes);
         GetSystemTime(&start_hours, &start_minutes, &start_seconds);
         //cur_time_s = start_hours*3600 + start_minutes*60 + start_seconds;
         time(&cur_time_s);
-        end_time_s = cur_time_s + step_minutes*60;
+        end_time_s = cur_time_s + (int) (runrec.nx * runrec.ny * runrec.nev/10000.0 * 4.5); // expected seconds
         printf("START:%2d-%2d-%2d (cur_time = %u s, end_time = %u s)\n", start_hours, start_minutes, start_seconds, cur_time_s, end_time_s);
+        SetCtrlAttribute(p1h, P1_EXPTIME, ATTR_DIMMED, 0);
+        SetCtrlVal(p1h, P1_EXPTIME, ctime(&end_time_s));
         //GetCtrlVal (p1h, P1_DSAVE, &dsave);
   //if (dsave) {
           GetCtrlVal (p1h, P1_DFILE, dfile0);
           fev=0;
           fcount=1;
           GetCtrlVal (p1h, P1_NEWF, &fmax);
         GetCtrlVal (p1h, P1_SUPR, &supr0);
         if (supr0) {
           GetCtrlVal (p1h, P1_TDCMIN, &tdcmin);
           GetCtrlVal (p1h, P1_FRAC, &frac);
+        }
 #ifdef USE_DAQ
         V673A_map(0,CAEN_V673A,0);
         V673A_init(0);
         V462_map(0,CAEN_V462,0);
         V462_set0(0,1);
 #endif
   //if (dsave) {
     sprintf(dfile,"%s_file%02d.dat",dfile0,fcount);
+        SetCtrlVal(p1h, P1_FCOUNT,fcount);
     fp = fopen (dfile, "wb");
     time (&runrec.time);
     status = fwrite (&runrec, 1, runrec.len, fp);
         //}
         if (supr0) {
           fseed = runrec.time & 0x7fffffff;
           Uniform (1, fseed, &fracg);
+        }
+          // vertical moving stage loop
+          for (posrec.iy=0; posrec.iy<runrec.ny; posrec.iy++) {
+        posrec.yset = runrec.y0 + posrec.iy*runrec.dy;
+#ifdef USE_MIKRO_Y
+        if (print) printf("MIKRO_MoveTo (2, y);%d\n",posrec.yset);
+        do {
+                MIKRO_MoveTo (MIKRO_Y, offsetY + posrec.yset);
+                MIKRO_GetPosition(MIKRO_Y, &mikroY);
+                        mikroY -= offsetY;
+                if (print) printf("%d\n", abs(posrec.yset - mikroY));
+                } while (abs(posrec.yset - mikroY) > STEP_TOLERANCE);
+                if (print) printf("->MIKRO_MoveTo (2, y);%d\n",posrec.yset);
+#endif
+        posrec.yset = mikroY;
+        SetCtrlVal (p1h, P1_Y, posrec.yset);
+        SetCtrlVal (p1h, P1_IY, posrec.iy);
+        // horizontal moving stage loop
-  for (posrec.ix=0; posrec.ix<runrec.nx; posrec.ix++) {
+    for (posrec.ix=0; posrec.ix<runrec.nx; posrec.ix++) {
-                posrec.xset = runrec.x0 + posrec.ix*runrec.dx;
+          posrec.xset = runrec.x0 + posrec.ix*runrec.dx;
 #ifdef USE_MIKRO
       do {
         if (print) printf("MIKRO_MoveTo (1, x);%d\n",posrec.xset);
-        MIKRO_MoveTo (MIKRO_X, posrec.xset);
+        MIKRO_MoveTo (MIKRO_X, offsetX + posrec.xset);
         MIKRO_GetPosition(MIKRO_X, &mikroX);
+                mikroX -= offsetX;
         if (print) printf("%d\n", abs(posrec.xset - mikroX));
         } while (abs(posrec.xset - mikroX) > STEP_TOLERANCE);
 #endif
-    posrec.xset = mikroX; // true value
+      posrec.xset = mikroX; // true value
           SetCtrlVal (p1h, P1_X, posrec.xset);
           SetCtrlVal (p1h, P1_IX, posrec.ix);
-          for (posrec.iy=0; posrec.iy<runrec.ny; posrec.iy++) {
-      posrec.yset = runrec.y0 + posrec.iy*runrec.dy;
-#ifdef USE_MIKRO_Y
-      if (print) printf("MIKRO_MoveTo (2, y);%d\n",posrec.yset);
-      do {
-        MIKRO_MoveTo (MIKRO_Y, posrec.yset);
-        MIKRO_GetPosition(MIKRO_Y, &mikroY);
-        if (print) printf("%d\n", abs(posrec.yset - mikroY));
-        } while (abs(posrec.yset - mikroY) > STEP_TOLERANCE);
-          if (print) printf("->MIKRO_MoveTo (2, y);%d\n",posrec.yset);
-#endif
-      posrec.yset = mikroY;
-      SetCtrlVal (p1h, P1_Y, posrec.yset);
-      SetCtrlVal (p1h, P1_IY, posrec.iy);
-      //if (dsave) {
+        //if (dsave) {
-      if (fmax && (ftell(fp) > fmax)) {
+        if (fmax && (ftell(fp) > fmax)) {
         fcount+=1;
         sprintf(dfile,"%s_file%02d.dat",dfile0,fcount);
         fclose(fp);
         fp = fopen (dfile, "wb");
       }
       time (&posrec.time);
             status = fwrite (&posrec, 1, posrec.len, fp);
             if (print) printf("POSREC status %d len %d\n", status, posrec.len);
                   //}
             count=0;
             ntics=0;
             if (print) printf("Event counter: %d\n", evtrec.nev);
                 //if (t!=told ) printf("%d in %2.2f min daq::event() %s\n",i, (double)(t-tstart)/60., ctime(&t));
 #ifdef USE_DAQ
 #ifdef CAEN_V673A
         V462_start0(0);
                 //tmlnk (tout);
+                  }
                   if (print) printf("Received triggers: %d \n", ntrig);
         } else {
           if (ntrig==0) {
                 fprintf(stderr,"-----------------> no trigger");
-                break;
+                //break;
           } else {
                 if (print) fprintf(stderr,"-----------------> wrong number of triggers in chips V673a  !!!\n");
+          }
           V673A_clallb(0);
           evtrec.nev--;
             if (fmax && (ftell(fp) > fmax)) {
               time (&endrec.time);
                       status = fwrite (&endrec, 1, endrec.len, fp);
               fcount+=1;
                   sprintf(dfile,"%s_file%02d.dat",dfile0,fcount);
+                          SetCtrlVal(p1h, P1_FCOUNT,fcount);
                   fclose(fp);
                   fp = fopen (dfile, "wb");
+            }
+        }
         GetCtrlVal(p1h, P1_DAQ, &daq_on);
             } while (evtrec.nev++<runrec.nev && daq_on);
             if (!daq_on) break;
-          } // x loop
+          } // y loop
           if (!daq_on) break;
-        } // y loop
+        } // x loop
   //if (dsave) {
     time (&endrec.time);
     status = fwrite (&endrec, 1, endrec.len, fp);
     if (print) printf("ENDREC status %d len %d\n", status, endrec.len);
 int __stdcall WinMain (HINSTANCE hInstance, HINSTANCE hPrevInstance,
                        LPSTR lpszCmdLine, int nCmdShow)
+{
   int i,j,status,refon;
   long int xpos, ypos, zpos;
-  char dfile[300];
+  char dfile[512];
+  char message[256];
   FILE *fp;
   typedef struct {
   unsigned long id,len;
   unsigned long fver,time;
 #endif
   if ((p1h = LoadPanel (0, "include/sipmScan_ui.uir", P1)) < 0) return -1;
   DisplayPanel (p1h);
   SetCtrlAttribute (p1h, P1_PLCH, ATTR_MAX_VALUE, NCH-1);
+  SetCtrlVal (p1h, P1_OFFSETX, offsetX);
+  SetCtrlVal (p1h, P1_OFFSETY, offsetY);
   GetCtrlVal(p1h, P1_ADCHLSAVE, &runrec.xy);
   SetCtrlVal(p1h, P1_ADCHL, runrec.xy-2);
   QueueUserEvent (1000, p1h, P1_RESET);
   QueueUserEvent (1000, p1h, P1_ASD8THR);
   QueueUserEvent (1000, p1h, P1_AMPDISSHTHR);
   QueueUserEvent (1000, p1h, P1_BIAS);
+  QueueUserEvent (1000, p1h, P1_OFFSETX);
+  QueueUserEvent (1000, p1h, P1_OFFSETY);
   printf("BIAS offset set to: %d\n", POWERSUPPLY);
   do {
     GetUserEvent (1, &pID, &rID);
         switch (rID) {
           case P1_DAQ:
                 GetCtrlVal (p1h, P1_DAQ, &daq_on);
                 if (daq_on) {
                   CmtScheduleThreadPoolFunction (poolHandle, daq_run, (void *)&dummy, &tfID);
                 } else {
-          CmtWaitForThreadPoolFunctionCompletion (poolHandle, tfID,
+          CmtWaitForThreadPoolFunctionCompletion (poolHandle, tfID, OPT_TP_PROCESS_EVENTS_WHILE_WAITING);
-                                                                        OPT_TP_PROCESS_EVENTS_WHILE_WAITING);
               CmtReleaseThreadPoolFunctionID (poolHandle, tfID);
+            }
             break;
+          case P1_XMIN:
+                  if (!daq_on) {
+                          GetCtrlVal (p1h, P1_XMIN, &xpos);
+                          //backlash
+                          MIKRO_MoveTo(MIKRO_X, xpos+offsetX-1000);
+                          Delay(0.01);
+                          MIKRO_MoveTo(MIKRO_X, xpos+offsetX);
+                  }
+                  break;
+          case P1_YMIN:
+                  if (!daq_on) {
+                          GetCtrlVal(p1h, P1_YMIN, &ypos);
+                          MIKRO_MoveTo(MIKRO_Y, ypos+offsetY-1000);
+                          Delay(0.01);
+                          MIKRO_MoveTo(MIKRO_Y, ypos+offsetY);
+                  }
+                  break;
           case P1_ZSET:
             if (!daq_on) {
               GetCtrlVal (p1h, P1_ZSET, &zpos);
 #ifdef USE_MIKRO_Z
           MIKRO_MoveTo (MIKRO_Z, zpos);
             break;
           case P1_MGET:
           #ifdef USE_MIKRO
                 MIKRO_GetPosition(MIKRO_X,&xpos);
                 Delay(0.01);
-                SetCtrlVal (p1h, P1_X, xpos);
+                SetCtrlVal (p1h, P1_X, xpos-offsetX);
           #ifdef USE_MIKRO_Y
                 MIKRO_GetPosition(MIKRO_Y,&ypos);
                 Delay(0.01);
-                SetCtrlVal (p1h, P1_Y, ypos);
+                SetCtrlVal (p1h, P1_Y, ypos-offsetY);
           #endif
           #ifdef USE_MIKRO_Z
                 MIKRO_GetPosition(MIKRO_Z,&zpos);
                 Delay(0.01);
                 SetCtrlVal (p1h, P1_Z, zpos);
                   GetCtrlVal (p1h, P1_SBIAS, &value);
                   value -= POWERSUPPLY;
                   SetDac(SIPM, value);
+                }
                 break;
+                // reset SiPM bias to initial power supply voltage
+                case P1_BIASOFF:
+                        {
+                                SetDac(SIPM, 0);
+                                SetCtrlVal(p1h, P1_SBIAS, (double) POWERSUPPLY);
+                        }
+                case P1_OFFSETX:
+                        {
+                                GetCtrlVal(p1h, P1_OFFSETX, &offsetX);
+                        }
+                case P1_OFFSETY:
+                        {
+                                GetCtrlVal(p1h, P1_OFFSETY, &offsetY);
+                        }
+        }
   } while ((rID != P1_EXIT)||daq_on);

Subversion Repositories f9daq

(root)/sipmScan/src/sipmScan.c – Rev 39 → 90