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

 #  endif
 #  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];
   int ntrig;
   int chip, len;
   unsigned long *pdata;
   FILE *fp;
-  #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;
         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;
           do {
         MIKRO_MoveTo (MIKRO_X, posrec.xset);
         MIKRO_GetPosition(MIKRO_X, &mikroX);
       } while (abs(posrec.xset - mikroX) > STEP_TOLERANCE);
           //printf("->MIKRO_MoveTo (1, x);%d\n",posrec.xset);
 #endif
           //SetCtrlVal (p1h, P1_X, posrec.xset);
           posrec.xset = mikroX; // set the true value
           SetCtrlVal(p1h, P1_X, mikroX);
           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
       //printf("MIKRO_MoveTo (2, y);%d\n",y);
       do {
         MIKRO_MoveTo (MIKRO_Y, posrec.yset);
         MIKRO_GetPosition(MIKRO_Y, &mikroY);
         } while (abs(posrec.yset - mikroY) > STEP_TOLERANCE);
           //printf("->MIKRO_MoveTo (2, y);%d\n",posrec.yset);
 #endif
             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 (!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);
   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",
                                     VAL_LOAD_BUTTON, 0, 0, 1, 0, dfile);
                           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);
           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