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

 #  ifdef WIENVME
 #    include "wienvme_dll.h"
 #  endif
 # include "CAENV673A.h"
 # include "CAENV462.h"
 #endif
 #ifdef USE_MIKRO
 #  include "MIKRO.h"
 #endif
 #include <userint.h>
 #include <ansi_c.h>
 #include <utility.h>
 #include <analysis.h>
 #include "c:\HOME\dino\sipmScan\include\sipmBias_ui.h"
 #ifdef USE_DAQ
 #  define USE_CAMAC
 #  include "camac.h"
+#  define NDAC 1 // CAEN C221
+#    define ASD8 0 // channels
+#    define AMPDISSH 1
-#  define NDAC 1
+#    define SIPM 2
 #  define CAEN_V673A   0x22220000 // IJS
 //#define CAEN_V673A   0x10110000 // FMF
 #  define CAEN_V462    0x100300
 #endif
 #define MAX_THREADS 10
 #define IWAIT 200
 #define NCH 64
+#define POWERSUPPLY 66 // Power supply voltage for SiPM bias
 static int p1h, pID, rID, tfID;
 static int ph_tdc, ph_adc;
 static int dtdc[NCH][2][MAXCH];
 static int something[NCH][2][MAXCH];
-  #define RUNREC_ID 1
+  #define  RUNREC_ID 1
-  #define ENDREC_ID 2
+  #define  ENDREC_ID 2
-  #define POSREC_ID 3
+  #define  POSREC_ID 3
-  #define EVTREC_ID 4
+  #define  EVTREC_ID 4
-  #define THRREC_ID 5
+  #define  THRREC_ID 5
+  #define BIASREC_ID 6
   typedef struct {
     unsigned long id,len;
     unsigned long fver,time;
+    unsigned long thLow, thUp, thStep;
+    unsigned long biasLow, biasUp, biasStep;
     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;
     unsigned long time;
     long ix,x,xset,iy,y,yset;
     } POSREC;
         POSREC posrec;
   typedef struct {
     unsigned long id,len;
     unsigned long nev;
     //unsigned short data[NCH];
     } EVTREC;
         EVTREC evtrec;
   typedef struct {
         unsigned long id;
         unsigned long len;
-        unsigned long thLow, thUp, thStep;
         unsigned long threshold;
         } THRREC;
   THRREC thrrec;
+  typedef struct {
+        unsigned long id;
+        unsigned long len;
+        unsigned long bias;
+        } BIASREC;
+  BIASREC biasrec;
         double thresholdLow;
         double thresholdUp;
         double thresholdStep;
+        double biasLow;
+        double biasUp;
+        double biasStep;
 #define BSIZE 20000
         int maxn = BSIZE-4; // 2 words x 2 headers
         int tdcdata, edge_type, overflow, tdc_num, channel, ev_dat, last_dat, nval_dat;
         posrec.len = sizeof(posrec);
         evtrec.id = EVTREC_ID;
         evtrec.len = sizeof(evtrec);
         thrrec.id = THRREC_ID;
         thrrec.len = sizeof(thrrec);
+        biasrec.id = BIASREC_ID;
+        biasrec.len = sizeof(biasrec);
         cres = 0;
         GetCtrlVal (p1h, P1_NEVE, &runrec.nev);
         if (supr0) {
           GetCtrlVal (p1h, P1_TDCMIN, &tdcmin);
           GetCtrlVal (p1h, P1_FRAC, &frac);
+        }
-        GetCtrlVal (p1h, P1_STHRLOW, &thresholdLow);
+        GetCtrlVal (p1h, P1_SASD8THR, &thresholdLow);
         runrec.thLow = thresholdLow*1000; // mV
         GetCtrlVal (p1h, P1_STHRUP, &thresholdUp);
         runrec.thUp = thresholdUp*1000; // mV
         GetCtrlVal (p1h, P1_STEP, &thresholdStep);
         runrec.thStep = thresholdStep*1000; // mV
+        GetCtrlVal (p1h, P1_SBIAS, &biasLow);
+        GetCtrlVal (p1h, P1_SBIASUP, &biasUp);
+        GetCtrlVal (p1h, P1_SBIASSTEP, &biasStep);
+        runrec.biasLow  = biasLow*1000;
+        runrec.biasUp   = biasUp*1000;
+        runrec.biasStep = biasStep*1000;
 #ifdef USE_DAQ
         V673A_map(0,CAEN_V673A,0);
         V673A_init(0);
         if (supr0) {
           fseed = runrec.time & 0x7fffffff;
           Uniform (1, fseed, &fracg);
+        }
+  // Set SiPM bias voltage
+  //SetDac(SIPM, 1.0);
+  biasrec.bias = runrec.biasLow;
+  while (biasrec.bias <= runrec.biasUp) {
+        SetDac(SIPM, biasrec.bias/1000.0 - POWERSUPPLY);
+        SetCtrlVal(p1h, P1_SBIAS, biasrec.bias/1000.0);
+        status = fwrite(&biasrec, 1, biasrec.len, fp);
+        if (print) printf("BIASREC status %d len %d bias %d\n", status, biasrec.len, biasrec.bias);
+  // Set ASD threshold voltage
   thrrec.threshold = (unsigned long) ceil(thresholdLow*1000); // threshold in mV
   while (thrrec.threshold <= thresholdUp*1000) {
-    SetDac(0, thrrec.threshold/1000.0);
+    SetDac(ASD8, thrrec.threshold/1000.0);
-        SetCtrlVal(p1h, P1_STHRLOW, thrrec.threshold/1000.0);
+        SetCtrlVal(p1h, P1_SASD8THR, thrrec.threshold/1000.0);
-    //wait_loop(100000);
     status = fwrite(&thrrec, 1, thrrec.len, fp);
     if (print) printf("THRREC status %d len %d threshold %d\n", status, thrrec.len, thrrec.threshold);
   for (posrec.ix=0; posrec.ix<runrec.nx; posrec.ix++) {
                 posrec.xset = runrec.x0 + posrec.ix*runrec.dx;
       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);
                   //}
       do {
         unsigned long word[BSIZE];
             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
         //if (timer_out) {
         //   fprintf(stderr,"-----------------> V673a timeout !!!\n");
         //   V673A_clallb(0);
         //   continue;
         //}
         if (ntics>=2) {
            fprintf(stderr,"-----------------> V673a timeout !!!\n");
            printf("Timer tics %d\n", ntics);
            ntics = 0;
            evtrec.nev--;
             else          len = V673A_read0 (0, &word[count+2],maxn-count-2);
             if (len>16*32*4){
               if (print) printf("Length > 2048 per chip!\n");
               chip=2;
               V673A_clallb(0);
               evtrec.nev--;
               continue;
+            }
             pdata=&word[count+2];
             for(i=0;i< len;i++){
           evtrec.nev--;
           continue;
+        }
         V673A_clallb(0);
 #endif
 #endif
         GetCtrlVal(p1h, P1_ADCHLSAVE, &runrec.xy);
         evtrec.len = count*sizeof(unsigned long) + sizeof(evtrec);
         if (print) printf("EVTREC.len = %d bytes.\n", evtrec.len);
         status = fwrite (&evtrec, 1, sizeof(evtrec), fp);
         status = fwrite (&word, 1, count*sizeof(unsigned long), fp);
                       status = fwrite (&endrec, 1, endrec.len, fp);
               fcount+=1;
                   sprintf(dfile,"%s_file%02d.dat",dfile0,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);
                                   if(cur_time_s >= end_time_s) {
                                           end_time_s = cur_time_s + step_minutes*60;
                                           printf("STEP (nev):%2d-%2d-%2d @ %d\n", start_hours, start_minutes, start_seconds, posrec.xset);
                         break;
                 }
         }
         //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;
         } // y loop
         if (!daq_on) break;
         thrrec.threshold += thresholdStep*1000;
         } // threhsold loop
+        if (!daq_on) break;
+        biasrec.bias += runrec.biasStep;
+        } // bias loop
+        SetDac(SIPM, 0); // safety
   //if (dsave) {
     time (&endrec.time);
     status = fwrite (&endrec, 1, endrec.len, fp);
     if (print) printf("ENDREC status %d len %d\n", status, endrec.len);
     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
         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;
   SetStdioWindowSize (200, 560);
   SetStdioWindowPosition (825, 250);
 #ifdef USE_DAQ
   VME_START(NULL);
 #endif
 #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/sipmBias_ui.uir", P1)) < 0) return -1;
   DisplayPanel (p1h);
   SetCtrlAttribute (p1h, P1_PLCH, ATTR_MAX_VALUE, NCH-1);
   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);
   do {
     GetUserEvent (1, &pID, &rID);
         switch (rID) {
       case P1_TIMER:
         ntics+=1;
           case P1_ZSET:
             if (!daq_on) {
               GetCtrlVal (p1h, P1_ZSET, &zpos);
 #ifdef USE_MIKRO_Z
           MIKRO_MoveTo (MIKRO_Z, zpos);
 #endif
         }
             break;
           case P1_REREAD:
             if (!daq_on) {
           status = FileSelectPopup ("", "*.dat", ".dat",
                                     "Choose the previous file",
                           SetCtrlVal (p1h, P1_NY, runrec.ny);
                           SetCtrlVal (p1h, P1_YSTEP, runrec.dy);
                           SetCtrlVal (p1h, P1_YMIN, runrec.y0);
                           SetCtrlVal (p1h, P1_NEVE, runrec.nev);
+            }
               }
         }
             break;
           case P1_MGET:
           #ifdef USE_MIKRO
                 MIKRO_GetPosition(MIKRO_X,&xpos);
                 Delay(0.01);
                 SetCtrlVal (p1h, P1_X, xpos);
           #ifdef USE_MIKRO_Y
                 MIKRO_GetPosition(MIKRO_Y,&ypos);
                 Delay(0.01);
                 SetCtrlVal (p1h, P1_Y, ypos);
           #endif
           #ifdef USE_MIKRO_Z
                 MIKRO_GetPosition(MIKRO_Z,&zpos);
                 Delay(0.01);
                   MIKRO_ReferenceMove (MIKRO_Y);
                 #endif
                 #ifdef USE_MIKRO_Z
                   MIKRO_ReferenceMove (MIKRO_Z);
                 #endif
 #endif
                   SetWaitCursor (0);
+            }
             break;
           case P1_RESET:
         for (j=0;j<NCH;j++) {
           for (i=0;i<MAXCH;i++){
             dtdc[j][0][i]=0;
             dtdc[j][1][i]=0;
+          }
+        }
                 update_plots();
             break;
           case P1_TDCLOG:
           case P1_ADCLOG:
             GetCtrlVal (p1h, P1_ADCLOG, &status);
         SetCtrlAttribute (p1h, P1_ADC, ATTR_YMAP_MODE, status);
                 update_plots();
             break;
-          /*case P1_ASD8THR:
+          case P1_ASD8THR:
+            {
               double value;
               GetCtrlVal (p1h, P1_SASD8THR, &value);
-              SetDac(0, value);
+              SetDac(ASD8, value);
+            }
-            break;*/
+            break;
+          case P1_BIAS:
+                {
+                  double value;
+                  GetCtrlVal (p1h, P1_SBIAS, &value);
+                  value -= POWERSUPPLY;
+                  SetDac(SIPM, value);
+                }
+                break;
           case P1_AMPDISSHTHR:
+            {
               double value;
               GetCtrlVal (p1h, P1_SAMPDISSHTHR, &value);
-              SetDac(1, value);
+              SetDac(AMPDISSH, value);
+            }
             break;
+        }

Subversion Repositories f9daq

(root)/sipmScan/src/sipmBias.c – Rev 33 → 34