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 CVICALLBACK cb_timer (int panel, int control, int event, void *callbackData,
                            int eventData1, int eventData2)
+{
   QueueUserEvent (9000, p1h, P1_TIMER);
   return (0);
 }
 int update_plots (void)
+{
   int irange, ch;
   GetCtrlVal (p1h, P1_PLCH, &ch);
   if (ph_adc>0) DeleteGraphPlot (p1h, P1_ADC, ph_adc, VAL_DELAYED_DRAW);
   GetCtrlVal (p1h, P1_ADCHL, &irange);
   ph_adc = PlotY (p1h, P1_ADC, &something[ch][irange], MAXCH, VAL_INTEGER,
                   VAL_VERTICAL_BAR, VAL_EMPTY_SQUARE, VAL_SOLID, 1, VAL_BLUE);
   return (0);
 }
 int module_header(int recid,unsigned long *data,int len){
    data[0] = recid;
    data[1] = (len >0)? len : 0 ;
   int timer_out = 0;
   //time_t t,told, tstart;
   int ntrig;
   int chip, len;
   unsigned long *pdata;
   FILE *fp;
   #define RUNREC_ID 1
   #define ENDREC_ID 2
   #define POSREC_ID 3
   #define EVTREC_ID 4
   typedef struct {
     unsigned long id,len;
     unsigned long fver,time;
     unsigned long nev,nch,ped,xy;
     long nx,x0,dx,ny,y0,dy;
     } RUNREC;
         RUNREC runrec;
   typedef struct {
     unsigned long id,len;
     unsigned long time;
     } ENDREC;
         ENDREC endrec;
   typedef struct {
     unsigned long id,len;
         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_YMIN, &runrec.y0);
         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);
         }
-  for (posrec.ix=0; posrec.ix<runrec.nx; posrec.ix++) {
-                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_GetPosition(MIKRO_X, &mikroX);
-        if (print) printf("%d\n", abs(posrec.xset - mikroX));
-        } while (abs(posrec.xset - mikroX) > STEP_TOLERANCE);
-#endif
-    posrec.xset = mikroX; // true value
+          // vertical moving stage loop
-          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;
+        posrec.yset = runrec.y0 + posrec.iy*runrec.dy;
 #ifdef USE_MIKRO_Y
-      if (print) printf("MIKRO_MoveTo (2, y);%d\n",posrec.yset);
+        if (print) printf("MIKRO_MoveTo (2, y);%d\n",posrec.yset);
-      do {
+        do {
-        MIKRO_MoveTo (MIKRO_Y, posrec.yset);
+                MIKRO_MoveTo (MIKRO_Y, offsetY + posrec.yset);
-        MIKRO_GetPosition(MIKRO_Y, &mikroY);
+                MIKRO_GetPosition(MIKRO_Y, &mikroY);
+                        mikroY -= offsetY;
-        if (print) printf("%d\n", abs(posrec.yset - mikroY));
+                if (print) printf("%d\n", abs(posrec.yset - mikroY));
-        } while (abs(posrec.yset - mikroY) > STEP_TOLERANCE);
+                } while (abs(posrec.yset - mikroY) > STEP_TOLERANCE);
-          if (print) printf("->MIKRO_MoveTo (2, y);%d\n",posrec.yset);
+                if (print) printf("->MIKRO_MoveTo (2, y);%d\n",posrec.yset);
 #endif
-      posrec.yset = mikroY;
+        posrec.yset = mikroY;
-      SetCtrlVal (p1h, P1_Y, posrec.yset);
+        SetCtrlVal (p1h, P1_Y, posrec.yset);
-      SetCtrlVal (p1h, P1_IY, posrec.iy);
+        SetCtrlVal (p1h, P1_IY, posrec.iy);
+        // horizontal moving stage loop
+    for (posrec.ix=0; posrec.ix<runrec.nx; posrec.ix++) {
+          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, 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
+          SetCtrlVal (p1h, P1_X, posrec.xset);
+          SetCtrlVal (p1h, P1_IX, posrec.ix);
-      //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);
                   //}
+                  }
                   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--;
           continue;
         }
         V673A_clallb(0);
 #endif
 #endif
         GetCtrlVal(p1h, P1_ADCHLSAVE, &runrec.xy);
         evtrec.len = count*sizeof(unsigned long) + sizeof(evtrec);
         status = fwrite (&evtrec, 1, sizeof(evtrec), fp);
         status = fwrite (&word, 1, count*sizeof(unsigned long), fp);
             if (!(evtrec.nev%1000)) SetCtrlVal (p1h, P1_CEVE, evtrec.nev);
             if (!((evtrec.nev+1)%1000)) {
             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");
             }
         }
         if(step_minutes > 0) {
                                   GetSystemTime(&start_hours, &start_minutes, &start_seconds);
                                   //cur_time_s = start_hours*3600 + start_minutes*60 + start_seconds;
                                   time(&cur_time_s);
         //told=t;
         GetCtrlVal(p1h, P1_DAQ, &daq_on);
             } while (evtrec.nev++<runrec.nev && daq_on);
-            if (!daq_on) break;
-          } // x loop
-          if (!daq_on) break;
+            if (!daq_on) break;
-        } // y loop
+          } // y loop
+          if (!daq_on) break;
+        } // 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);
     fclose(fp);
     daq_on=0;
         SetCtrlVal (p1h, P1_DAQ, daq_on);
         SetCtrlVal (p1h, P1_CEVE, evtrec.nev);
         if (print) printf("RUN END. \n\n");
         return 0;
 }
 int SetDac(int ch, double val){
         const unsigned int maxdac=0xFFF;
         const double RANGE = +9.9976; //V
         unsigned int dac;
         dac = (val/RANGE)*maxdac;
         CSSA_W(NDAC,ch,16,&dac);
         printf("DAC ch %d set to %f V dac=%d\n", ch, val, dac);
         return 0;
+}
 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;
   } RUNREC;
   RUNREC runrec;
   if (InitCVIRTE (hInstance, 0, 0) == 0)
         return -1;      /* out of memory */
   SetSleepPolicy(VAL_SLEEP_MORE);
   CmtNewThreadPool (MAX_THREADS, &poolHandle);
   SetStdioWindowOptions (1000, 0, 0);
   SetStdioWindowSize (200, 560);
 #ifdef USE_MIKRO
   MIKRO_Open (MIKRO_COM);
   MIKRO_Init (MIKRO_X,0);
 #ifdef USE_MIKRO_Y
   MIKRO_Init (MIKRO_Y,0);
 #endif
 #ifdef USE_MIKRO_Z
   MIKRO_Init (MIKRO_Z,0);
 #endif
 #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);
         GetCtrlVal (p1h, P1_REFON, &refon);
                     if (refon) update_plots();
             break;
       case P1_REFRESH:
                 update_plots();
             break;
           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 @ 39 – Rev 39 → 90