WebSVN – f9daq – Diff – /cvi/instr/DRS/drs4.c

+#include "H1D.h"
+#include "H2D.h"
 #include <formatio.h>
 #include <utility.h>
 #include <ansi_c.h>
 #include <cvirte.h>
 #include <userint.h>
 #include "MIKRO.h"
 #endif
 static int daq_on;
 static int ph, p2,plothandle[4]= {0,0,0,0};
+static int scanplothandle;
 static int tfID;
 static int controlID;
 #define MAX_THREADS 10
 static CmtThreadPoolHandle poolHandle = 0;
+static float gSum[4]={0,0,0,0};
 int main (int argc, char *argv[]) {
         short port;
         if (InitCVIRTE (0, argv, 0) == 0)
                 return -1;  /* out of memory */
         if ((ph = LoadPanel (0, "drs4.uir", PANEL)) < 0)
         return ;
+}
+static float xs[4][1024];
 int CVICALLBACK daq(void *functionData) {
         int neve;
         char filename[0xff];
         int pfreq;
         int enabledoutput;
         int neveold = 0;
         double rate;
         int *args =  (int *) functionData;
+        double twin[2];
-        GetCtrlVal(ph, PANEL_FILENAME, filename );
         GetCtrlVal(ph, PANEL_CH0, &imask[0] );
         GetCtrlVal(ph, PANEL_CH1, &imask[1]  );
         GetCtrlVal(ph, PANEL_CH2, &imask[2]  );
         GetCtrlVal(ph, PANEL_CH3, &imask[3]  );
         mask = 0;
         GetCtrlVal(ph,PANEL_TRGTYPE, &trgtype);
         GetCtrlVal(ph,PANEL_TRGLEVEL, &trglevel);
         GetCtrlVal(ph,PANEL_TRGPOLARITY, &trgpolarity);
         GetCtrlVal(ph,PANEL_RANGE, &range);
+        GetCtrlVal(ph,PANEL_TWIN0, &twin[0]);
+  GetCtrlVal(ph,PANEL_TWIN1, &twin[1]);
         //printf("mask=0x%x\n",mask);
         DRSSetMask( (unsigned char)( mask & 0xF ) );
         DRSSetFrequency( frequency  );
         DRSSetTriggerType( trgtype );
         DRSSetTriggerLevel(trglevel);
         DRSSetTriggerPolarity(trgpolarity);
         DRSSetRange (  range );
         FILE *fp= NULL;
         if (enabledoutput) {
+                if (args[0]) {
+                        GetCtrlVal(ph, PANEL_FILENAME, filename );
-                if (args[0]) fp =  fopen(filename,"wb");
+                        fp =  fopen(filename,"wb");
+          } else {
+                        GetCtrlVal(p2, SCAN_FILENAME, filename );
-          else                     fp =  fopen(filename,"ab");
+                        fp =  fopen(filename,"ab");
+                }
         }
         static unsigned char *buffer;
         int buffer_size = 0;
         const int nBoards=1;
         const int waveDepth=1024;
         time_t t=0,told=0, tstart=0;
         if (!DRSInit()) {
                 time(&tstart);
                 told=tstart;
-                int i=0;
+                int nev=0;
+                for (int k = 0;k<4;k++){
+                         gSum[k]=0;
+                }
-                for (i=0; i<neve; i++) {
+                for (int i=0; i<neve; i++) {
                         start_timer(1);// 1 s timeout
                         int retval = DRSRead(0);
                         stop_timer();
                         int nb =  ( retval == 0 && fp ) ? DRSToBuffer( buffer , i ) : 0;
                         SetCtrlVal(ph,PANEL_CEVE,i);
                         if (retval) i--;
+                        nev++;
                         if (!daq_on) break;
                         time(&t);
                         if (t!=told ) {
                                 rate =  (i-neveold);
                                 printf("%d events in %2.2f min (%d s) Rate %f Hz %s   ",i+1, (double)(t-tstart)/60.,(t-tstart), rate , ctime(&t));
+                        }
                         told=t;
 // Save data
                         if (nb>0 && fp) fwrite(buffer, 1,nb ,fp);
 // Plot Data
-                        if (i % pfreq == 0) for (int k=0; k<4; k++) {
+                         for (int k=0; k<4; k++) {
-                                        if (! (mask & ( 0x1<<k ))  ) continue;
+                                        if ( (mask & ( 0x1<<k ))  ){
-                                        float *t=DRSGetTime(k);
+                                          float *t=DRSGetTime(k);
-                                        float *x=DRSGetWave(k);
+                                          float *x=DRSGetWave(k);
+                                                                 ;
+                                          for (int j=0 ; j<1024 ; j++) {
+                                                  xs[k][j]= x[j]*1e-3;
+                                                  if (t[j]> twin[0] && t[j] < twin[1]) gSum[k]+= fabs(xs[k][j]);
+                                                  if (verbose) printf("[%d] %d. x= %3.2f  y=%3.2f\n", k, i, t[j], x[j] );
+                                                  //h[k]->Fill( t[i], x[i]*1e-3);
+                                          }
+                                          if (i % pfreq == 0) {
-                                        const int col[4]= {VAL_WHITE,VAL_RED,VAL_GREEN,VAL_BLUE};
+                                            const int col[4]= {VAL_WHITE,VAL_RED,VAL_GREEN,VAL_BLUE};
-                                        if (plothandle[k])      DeleteGraphPlot (ph, PANEL_GRAPH, plothandle[k], VAL_IMMEDIATE_DRAW);
+                                          if (plothandle[k])    DeleteGraphPlot (ph, PANEL_GRAPH, plothandle[k], VAL_IMMEDIATE_DRAW);
-                                        plothandle[k] = PlotXY (ph, PANEL_GRAPH, t, x, 1024, VAL_FLOAT, VAL_FLOAT, VAL_THIN_LINE, VAL_NO_POINT, VAL_SOLID, 1, col[k]);
+                                            plothandle[k] = PlotXY (ph, PANEL_GRAPH, t, xs[k], 1024, VAL_FLOAT, VAL_FLOAT, VAL_THIN_LINE, VAL_NO_POINT, VAL_SOLID, 1, col[k]);
-                                        for (int i=0 ; i<1024 ; i++) {
+                                          }
-                                                if (verbose) printf("[%d] %d. x= %3.2f  y=%3.2f\n", k, i, t[i], x[i] );
-                                                //h[k]->Fill( t[i], x[i]*1e-3);
+                                        }
                                 }
                 }
                 time(&t);
-                printf("%d events in %2.2f min (%d s) %s",i+1, (double)(t-tstart)/60.,t-tstart, ctime(&t));
+                printf("%d events in %2.2f min (%d s) %s",nev, (double)(t-tstart)/60.,t-tstart, ctime(&t));
                 DRSEnd();
+        }
         if (fp) fclose(fp);
         free(buffer);
         return 0;
+}
         int  n[3];
         char filename[0xFF];
         int enabledoutput;
         FILE *fp;
-        GetCtrlVal(ph, PANEL_FILENAME, filename );
+        GetCtrlVal(p2, SCAN_FILENAME, filename );
         GetCtrlVal(ph,PANEL_ENABLEDOUTPUT, &enabledoutput);
         if ( GetFileInfo(filename,&size) ) {
-                 MessagePopup ("Warning", "File exist. Remove it first or choose another file");
+                MessagePopup ("Warning","File exist. Remove it first or choose another file");
-           return 0;
+          return 0;
+        }
         GetCtrlVal(p2, SCAN_STEPX, &dx[0]);
         GetCtrlVal(p2, SCAN_STEPY, &dx[1]);
         GetCtrlVal(p2, SCAN_NSTEPSX, &nx[0]);
         GetCtrlVal(p2, SCAN_NSTEPSY, &nx[1]);
         GetCtrlVal(p2, SCAN_STARTX, &x0[0]);
         GetCtrlVal(p2, SCAN_STARTY, &x0[1]);
+        for (int k=0;k<4;k++)  {
+                H2D_Init(k, "charge","Induced charge",
+                                                nx[0], x0[0] - dx[0]*0.5 ,x0[0] + dx[0] * ( nx[0] - 0.5) ,
+                                                nx[1], x0[1] - dx[1]*0.5 ,x0[1] + dx[1] * ( nx[1] - 0.5));
+          H1D_Init(k, "charge","Induced charge projection x",
+                                                nx[0], x0[0] - dx[0]*0.5 ,x0[0] + dx[0] * ( nx[0] - 0.5) );
+          H1D_Init(100 + k, "charge","Induced charge projection y",
+                                                nx[1], x0[1] - dx[1]*0.5 ,x0[1] + dx[1] * ( nx[1] - 0.5) );
+        }
         if (enabledoutput) {
                                 fp =  fopen(filename,"ab");
                                 if (fp) {
                                          size=36;
                                          fwrite(runbuf, 1,4 ,fp);
+                                }
+                        }
 #endif
                         daq(functionData);
+                        for (int k=0;k<4;k++) {
+                                H2D_Fill(k, ix[0], ix[1] , gSum[k] );
+                                H1D_Fill(k, ix[0] , gSum[k] );
+                                H1D_Fill(100+ k, ix[1] , gSum[k] );
+                        }
+                        PlotScanHistogramCB(0,0,EVENT_COMMIT,NULL, 0,0);
                         if (!daq_on) break;
+                }
                 if (!daq_on) break;
+        }
                                                                                                  void *callbackData, int eventData1, int eventData2) {
         ThreadFunctionPtr mythread = NULL;
         switch (event) {
                 case EVENT_COMMIT:
                         controlID=0;
                         if (panel == ph && control == PANEL_START) {
                                 mythread = daq;
                                 controlID= control;
                         }
                         if (panel == p2 && control == SCAN_SCAN)   mythread = scan;
                         if (mythread!=NULL) {
                                 printf("New Thread panel=%d button=%d\n", panel, control);
                                 // SetDimming(1);
                                 daq_on=1;
                                 CmtScheduleThreadPoolFunctionAdv (poolHandle, mythread, &controlID,
                                                                                                                                                                         DEFAULT_THREAD_PRIORITY,
                                                                                                                                                                         EndOfThread,
                                                                                                                                                                         EVENT_TP_THREAD_FUNCTION_END,
                                                                                                                                                                         NULL, RUN_IN_SCHEDULED_THREAD,
                                                                                                                                                                         &tfID);
                         }
                         break;
         }
         return 0;
+}
 int CVICALLBACK StopCB (int panel, int control, int event,
                                                                                                 void *callbackData, int eventData1, int eventData2) {
         switch (event) {
                 case EVENT_COMMIT:
                         daq_on=0;
                         break;
         }
         return 0;
+}
 int CVICALLBACK ExitCB (int panel, int control, int event,
                                                                                                 void *callbackData, int eventData1, int eventData2) {
         switch (event) {
                 case EVENT_COMMIT:
                         QuitUserInterface (0);
                         break;
+        }
         return 0;
+}
                 case EVENT_COMMIT:
                         if (panel == p2) {
                                 switch (control) {
                                         case SCAN_BR :
                                                 axis = 1;
                                                 direction = 1;
                                                 GetCtrlVal(p2, SCAN_STEPX, &step);
                                                 break;
                                         case SCAN_BL :
                                                 axis = 1;
                                                 direction = -1;
                                                 GetCtrlVal(p2, SCAN_STEPX, &step);
                                                 break;
                                         case SCAN_BU :
                                                 axis = 2;
                                                 direction = 1;
                                                 GetCtrlVal(p2, SCAN_STEPY, &step);
                                                 break;
                                         case SCAN_BD :
                                                 axis = 2;
                                 MIKRO_GetPosition(axis,&n);
                                 if (axis == 1) SetCtrlVal (p2, SCAN_XP, n);
                                 if (axis == 2) SetCtrlVal (p2, SCAN_YP, n);
 #endif // MIKRO
+                        }
                         break;
         }
         return 0;
 }
 int CVICALLBACK GoXCB (int panel, int control, int event,
                                                                                          void *callbackData, int eventData1, int eventData2) {
         int n2;
                         MIKRO_MoveTo(1,n2);
                         MIKRO_GetPosition(1,&n2);
 #endif
                         SetCtrlVal (p2, SCAN_XP, n2);
                         break;
         }
         return 0;
 }
 int CVICALLBACK GoYCB (int panel, int control, int event,
                                                                                          void *callbackData, int eventData1, int eventData2) {
         int n2;
         switch (event) {
                 case EVENT_COMMIT:
                         MIKRO_MoveTo(2,n2);
                         MIKRO_GetPosition(2,&n2);
 #endif
                         SetCtrlVal (p2, SCAN_YP, n2);
                         break;
         }
         return 0;
 }
 int CVICALLBACK GetCurrentPositionCB (int panel, int control, int event,
                                                                                                                                         void *callbackData, int eventData1, int eventData2) {
         int n[2];
         switch (event) {
                 case EVENT_COMMIT:
 #ifdef MIKRO
                         MIKRO_GetPosition(1,&n[0]);
                         SetCtrlVal (p2, SCAN_XP, n[0]);
                         MIKRO_GetPosition(2,&n[1]);
                         break;
+        }
         return 0;
+}
 int CVICALLBACK HomeCB (int panel, int control, int event,
                                                                                                 void *callbackData, int eventData1, int eventData2) {
+        switch (event) {
+                case EVENT_COMMIT:
+#ifdef MIKRO
+                        MIKRO_ReferenceMove(1);
+                        MIKRO_ReferenceMove(2);
+                        GetCurrentPositionCB(panel, control, event, NULL, 0, 0);
+#endif
+                        break;
+        }
+        return 0;
+}
+int CVICALLBACK PlotScanHistogramCB (int panel, int control, int event,
+                                                                                                                                                 void *callbackData, int eventData1, int eventData2) {
+        int hid=0;
+        int nx;
+        int ny;
+        switch (event) {
+                case EVENT_COMMIT:
+                        GetCtrlVal (p2, SCAN_CHANNEL, &hid);
+                        GetCtrlVal(p2, SCAN_NSTEPSX, &nx);
+            GetCtrlVal(p2, SCAN_NSTEPSY, &ny);
+                        if (nx>1 && ny>1) {
+                                H2D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
+                                printf("redraw 2d\n");
+                        } else {
+                                if (nx>1)  {
+                                        printf("redraw 1d x\n");
+                                        H1D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
+                                }
+                          if (ny>1) {
+                                        printf("redraw 1d x\n");
+                                        H1D_Draw(100+hid,p2,SCAN_GRAPH,&scanplothandle);
+                                }
+                        }
+                        break;
+        }
+        return 0;
+}
+int CVICALLBACK OpenGuiCB (int panel, int control, int event,
+                                                                                                         void *callbackData, int eventData1, int eventData2) {
         switch (event) {
                 case EVENT_COMMIT:
-#ifdef MIKRO
-                        MIKRO_ReferenceMove(1);
+                        DisplayPanel (ph);
-                        MIKRO_ReferenceMove(2);
+            DisplayPanel (p2);
-                        GetCurrentPositionCB(panel, control, event, NULL, 0, 0);
-#endif
                         break;
+        }
         return 0;
+}

Subversion Repositories f9daq

(root)/cvi/instr/DRS/drs4.c – Rev 210 → 217