Subversion Repositories f9daq

Rev

Rev 210 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
217 f9daq 1
#include "H1D.h"
2
#include "H2D.h"
203 f9daq 3
#include <formatio.h>
195 f9daq 4
#include <utility.h>
5
#include <ansi_c.h>
203 f9daq 6
#include <cvirte.h>
195 f9daq 7
#include <userint.h>
8
#include "drs4.h"
203 f9daq 9
#include "drsread.h"
195 f9daq 10
 
203 f9daq 11
#define MIKRO 
12
 
13
#ifdef MIKRO
14
#include "MIKRO.h"
15
#endif
16
 
195 f9daq 17
static int daq_on;
203 f9daq 18
static int ph, p2,plothandle[4]= {0,0,0,0};
217 f9daq 19
static int scanplothandle;
195 f9daq 20
static int tfID;
21
static int controlID;
22
 
203 f9daq 23
#define MAX_THREADS 10
195 f9daq 24
 
203 f9daq 25
static CmtThreadPoolHandle poolHandle = 0;
195 f9daq 26
 
217 f9daq 27
static float gSum[4]={0,0,0,0};
195 f9daq 28
int main (int argc, char *argv[]) {
203 f9daq 29
        short port;
30
        if (InitCVIRTE (0, argv, 0) == 0)
31
                return -1;  /* out of memory */
32
        if ((ph = LoadPanel (0, "drs4.uir", PANEL)) < 0)
33
                return -1;
34
        if ((p2 = LoadPanel (0, "drs4.uir", SCAN)) < 0)
35
                return -1;
36
        SetStdioPort (CVI_STDIO_WINDOW);
37
        SetSleepPolicy(VAL_SLEEP_MORE);
38
        CmtNewThreadPool (MAX_THREADS,  &poolHandle);
39
 
40
        DisplayPanel (ph);
41
        DisplayPanel (p2);
42
 
43
#ifdef MIKRO
44
        GetCtrlVal(p2, SCAN_PORT, &port);
45
        if (MIKRO_Open (port)) MessagePopup ("Error", "Mikro Port Not found !\n Change in the GUI") ;
46
        MIKRO_Init(1,0);
47
        MIKRO_Init(2,0);
48
#endif
49
 
50
        RunUserInterface ();
51
        DiscardPanel (ph);
52
        DiscardPanel (p2);
53
        CmtDiscardThreadPool (poolHandle);
54
        MIKRO_Close ();
55
        return 0;
195 f9daq 56
}
57
 
58
 
210 f9daq 59
char strbuf[0xFF];
195 f9daq 60
 
210 f9daq 61
int gLog=0;
195 f9daq 62
 
210 f9daq 63
int printf(const char *format, ...) {
64
  va_list aptr;
65
  int ret;
66
  FILE *flog;
195 f9daq 67
 
210 f9daq 68
  va_start(aptr, format);
69
  ret = vsprintf(strbuf, format, aptr);
70
  va_end(aptr);
71
  SetCtrlVal(ph,PANEL_STDIO,strbuf);
72
 
73
  if (gLog) {
74
    flog = fopen ("stdio.log", "a");
75
    fprintf (flog, "%s", strbuf);
76
    fclose (flog);
77
  }
78
  return(ret);
79
}
80
 
81
 
82
 
195 f9daq 83
static void start_timer (double tout) {
203 f9daq 84
        SetCtrlAttribute (ph, PANEL_TIMER, ATTR_INTERVAL, tout);
85
        SetCtrlAttribute (ph, PANEL_TIMER, ATTR_ENABLED, 1);
195 f9daq 86
}
87
 
88
static void stop_timer ( void ) {
203 f9daq 89
        SetCtrlAttribute (ph, PANEL_TIMER, ATTR_ENABLED, 0);
90
        DRSSetTimeout();
195 f9daq 91
}
92
 
93
 
94
 
95
 
96
void CVICALLBACK EndOfThread ( CmtThreadPoolHandle poolhandle,
203 f9daq 97
                                                                                                                         CmtThreadFunctionID functionID, unsigned int event,
98
                                                                                                                         int value, void *callbackData  ) {
195 f9daq 99
 
203 f9daq 100
        daq_on=0;
101
        //SetDimming(0);
102
        printf("End of Thread \n");
103
        return ;
195 f9daq 104
 
105
}
106
 
107
 
217 f9daq 108
static float xs[4][1024];      
203 f9daq 109
int CVICALLBACK daq(void *functionData) {
195 f9daq 110
 
111
 
203 f9daq 112
        int neve;
113
        char filename[0xff];
114
        int imask[4];
115
        unsigned long mask;
116
        int frequency;
117
        double trgdelay;
118
        double trglevel;
119
        int trgtype;
120
        int trgchannel;
121
        int trgpolarity;
122
        int verbose;
123
        double range;
124
        int pfreq;
125
        int enabledoutput;
126
        int neveold = 0;
127
        double rate;
128
        int *args =  (int *) functionData;
217 f9daq 129
        double twin[2];
201 f9daq 130
 
217 f9daq 131
 
203 f9daq 132
        GetCtrlVal(ph, PANEL_CH0, &imask[0] );
133
        GetCtrlVal(ph, PANEL_CH1, &imask[1]  );
134
        GetCtrlVal(ph, PANEL_CH2, &imask[2]  );
135
        GetCtrlVal(ph, PANEL_CH3, &imask[3]  );
136
        mask = 0;
137
        for (int i=0; i<4; i++) {
138
                if (imask[i]) mask |= (1<<i);
201 f9daq 139
        }
203 f9daq 140
 
141
        GetCtrlVal(ph,PANEL_NEVE, &neve);
142
        GetCtrlVal(ph,PANEL_DEBUG, &verbose);
143
        GetCtrlVal(ph,PANEL_PFREQ, &pfreq);
144
        GetCtrlVal(ph,PANEL_ENABLEDOUTPUT, &enabledoutput);
145
 
146
        GetCtrlVal(ph,PANEL_FREQUENCY, &frequency);
147
        GetCtrlVal(ph,PANEL_TRGDELAY, &trgdelay);
148
        GetCtrlVal(ph,PANEL_TRGCHANNEL, &trgchannel);
149
 
150
        GetCtrlVal(ph,PANEL_TRGTYPE, &trgtype);
151
        GetCtrlVal(ph,PANEL_TRGLEVEL, &trglevel);
152
        GetCtrlVal(ph,PANEL_TRGPOLARITY, &trgpolarity);
153
        GetCtrlVal(ph,PANEL_RANGE, &range);
154
 
217 f9daq 155
        GetCtrlVal(ph,PANEL_TWIN0, &twin[0]);
156
  GetCtrlVal(ph,PANEL_TWIN1, &twin[1]);
203 f9daq 157
        //printf("mask=0x%x\n",mask);
158
 
159
        DRSSetMask( (unsigned char)( mask & 0xF ) );
160
 
161
        DRSSetFrequency( frequency  );
162
        DRSSetTriggerDelay(trgdelay );
163
        DRSSetTriggerChannel(trgchannel );
164
        DRSSetTriggerType( trgtype );
165
        DRSSetTriggerLevel(trglevel);
166
        DRSSetTriggerPolarity(trgpolarity);
167
 
168
        DRSSetRange (  range );
169
 
170
 
171
 
172
        FILE *fp= NULL;
173
 
174
        if (enabledoutput) {
217 f9daq 175
                if (args[0]) {
176
                        GetCtrlVal(ph, PANEL_FILENAME, filename );  
177
                        fp =  fopen(filename,"wb");
178
          } else {
179
                        GetCtrlVal(p2, SCAN_FILENAME, filename );  
180
                        fp =  fopen(filename,"ab");
181
                }      
203 f9daq 182
        }
183
 
184
        static unsigned char *buffer;
185
 
186
        int buffer_size = 0;
187
        const int nBoards=1;
188
        const int waveDepth=1024;
207 f9daq 189
 
203 f9daq 190
        if (buffer_size == 0) {
207 f9daq 191
         buffer_size =   8; // file header + time header
192
         buffer_size +=  nBoards * (4 + 4*(4+waveDepth*4)); // bin widths
193
         buffer_size += 24 + nBoards * (8 + 4*(8+waveDepth*2));
194
         buffer = (unsigned char *)malloc(buffer_size);
195
  }
196
 
203 f9daq 197
        time_t t=0,told=0, tstart=0;
198
 
199
        if (!DRSInit()) {
200
                time(&tstart);
201
                told=tstart;
217 f9daq 202
                int nev=0;
203
                for (int k = 0;k<4;k++){
204
                         gSum[k]=0;
205
                }
206
                for (int i=0; i<neve; i++) {
203 f9daq 207
                        start_timer(1);// 1 s timeout
208
                        int retval = DRSRead(0);
209
                        stop_timer();
210
                        int nb =  ( retval == 0 && fp ) ? DRSToBuffer( buffer , i ) : 0;
211
                        SetCtrlVal(ph,PANEL_CEVE,i);
212
                        if (retval) i--;
217 f9daq 213
                        nev++;
203 f9daq 214
                        if (!daq_on) break;
215
                        time(&t);
216
                        if (t!=told ) {
217
                                rate =  (i-neveold);
218
                                printf("%d events in %2.2f min (%d s) Rate %f Hz %s   ",i+1, (double)(t-tstart)/60.,(t-tstart), rate , ctime(&t));
219
                                GetCtrlVal(ph,PANEL_PFREQ, &pfreq);
220
                                neveold = i;
221
                        }
222
                        told=t;
223
// Save data
224
                        if (nb>0 && fp) fwrite(buffer, 1,nb ,fp);
225
// Plot Data
217 f9daq 226
                         for (int k=0; k<4; k++) {
227
                                        if ( (mask & ( 0x1<<k ))  ){
228
                                          float *t=DRSGetTime(k);
229
                                          float *x=DRSGetWave(k);
230
                                                                 ;
231
                                          for (int j=0 ; j<1024 ; j++) {
232
                                                  xs[k][j]= x[j]*1e-3;
233
 
234
                                                  if (t[j]> twin[0] && t[j] < twin[1]) gSum[k]+= fabs(xs[k][j]);
235
                                                  if (verbose) printf("[%d] %d. x= %3.2f  y=%3.2f\n", k, i, t[j], x[j] );
236
                                                  //h[k]->Fill( t[i], x[i]*1e-3);
237
                                          }
203 f9daq 238
 
217 f9daq 239
                                          if (i % pfreq == 0) {
240
 
241
                                            const int col[4]= {VAL_WHITE,VAL_RED,VAL_GREEN,VAL_BLUE};
242
                                          if (plothandle[k])    DeleteGraphPlot (ph, PANEL_GRAPH, plothandle[k], VAL_IMMEDIATE_DRAW);
243
                                            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]);
244
                                          }    
203 f9daq 245
                                        }
217 f9daq 246
 
203 f9daq 247
                                }
248
 
249
 
250
                }
251
                time(&t);
217 f9daq 252
                printf("%d events in %2.2f min (%d s) %s",nev, (double)(t-tstart)/60.,t-tstart, ctime(&t));
203 f9daq 253
                DRSEnd();
254
        }
255
 
256
        if (fp) fclose(fp);
257
 
258
        free(buffer);
259
 
260
        return 0;
261
 
195 f9daq 262
}
263
 
264
 
265
 
203 f9daq 266
int CVICALLBACK scan(void *functionData) {
195 f9daq 267
 
207 f9daq 268
        int dx[3]={0,0,0};
269
        int nx[3]={0,0,0};
270
        int x0[3]={0,0,0};
271
        int ix[3]={0,0,0};
272
        int idx[3]={0,0,0};  
203 f9daq 273
        int size;
207 f9daq 274
        char posrec[4]="POSR"; 
275
        char runbuf[4]="PRUN";
203 f9daq 276
 
207 f9daq 277
        int  n[3];
203 f9daq 278
        char filename[0xFF];
279
        int enabledoutput;
280
 
281
        FILE *fp;
217 f9daq 282
        GetCtrlVal(p2, SCAN_FILENAME, filename );
203 f9daq 283
 
284
        GetCtrlVal(ph,PANEL_ENABLEDOUTPUT, &enabledoutput);
285
 
286
        if ( GetFileInfo(filename,&size) ) {
217 f9daq 287
                MessagePopup ("Warning","File exist. Remove it first or choose another file");
288
          return 0;
203 f9daq 289
        }        
290
        GetCtrlVal(p2, SCAN_STEPX, &dx[0]);
291
        GetCtrlVal(p2, SCAN_STEPY, &dx[1]);
292
        GetCtrlVal(p2, SCAN_NSTEPSX, &nx[0]);
293
        GetCtrlVal(p2, SCAN_NSTEPSY, &nx[1]);
294
        GetCtrlVal(p2, SCAN_STARTX, &x0[0]);
295
        GetCtrlVal(p2, SCAN_STARTY, &x0[1]);
217 f9daq 296
        for (int k=0;k<4;k++)  {
297
                H2D_Init(k, "charge","Induced charge",
298
                                                nx[0], x0[0] - dx[0]*0.5 ,x0[0] + dx[0] * ( nx[0] - 0.5) ,
299
                                                nx[1], x0[1] - dx[1]*0.5 ,x0[1] + dx[1] * ( nx[1] - 0.5));  
300
          H1D_Init(k, "charge","Induced charge projection x",
301
                                                nx[0], x0[0] - dx[0]*0.5 ,x0[0] + dx[0] * ( nx[0] - 0.5) );  
302
          H1D_Init(100 + k, "charge","Induced charge projection y",
303
                                                nx[1], x0[1] - dx[1]*0.5 ,x0[1] + dx[1] * ( nx[1] - 0.5) );  
304
        }
207 f9daq 305
        if (enabledoutput) {
306
                                fp =  fopen(filename,"ab");
307
                                if (fp) {
308
                                         size=36;
309
                                         fwrite(runbuf, 1,4 ,fp);
310
                                         fwrite(&size  , 1,4 ,fp);
311
                                         fwrite(x0    , 1,4*3 ,fp);
312
                                         fwrite(dx    , 1,4*3 ,fp);
313
                                         fwrite(nx    , 1,4*3 ,fp);
314
                                         fclose(fp);
315
                                }    
316
        }      
317
 
203 f9daq 318
        for (int i=0; i<nx[0]; i++) {
319
 
320
                ix[0]= x0[0]+i*dx[0];
321
#ifdef MIKRO
322
                MIKRO_MoveTo(1,ix[0]);
323
#endif
324
                SetCtrlVal (p2, SCAN_IX, i);
325
                for (int j=0; j<nx[1]; j++) {
326
 
327
                        SetCtrlVal (p2, SCAN_IY, j);
328
 
329
                        ix[1]= x0[1]+j*dx[1];
330
#ifdef MIKRO
331
                        MIKRO_MoveTo(2,ix[1]);
332
 
333
                        MIKRO_GetPosition(1,&n[0]);
334
                        SetCtrlVal (p2, SCAN_XP, n[0]);
335
                        MIKRO_GetPosition(2,&n[1]);
336
                        SetCtrlVal (p2, SCAN_YP, n[1]);
207 f9daq 337
 
338
                        if (enabledoutput) {
203 f9daq 339
                                fp =  fopen(filename,"ab");
340
                                if (fp) {
207 f9daq 341
                                         idx[0]=i;
342
                                         idx[1]=j;
343
                                         size=24;
344
                                         fwrite(posrec, 1,4 ,fp);
345
                                         fwrite(&size , 1,4 ,fp);
346
                                         fwrite(n     , 1,4*3 ,fp);
347
                                         fwrite(idx    , 1,4*3 ,fp);
203 f9daq 348
                                         fclose(fp);
349
                                }    
350
                        }      
351
 
352
#endif
353
                        daq(functionData);
217 f9daq 354
 
355
                        for (int k=0;k<4;k++) {
356
                                H2D_Fill(k, ix[0], ix[1] , gSum[k] );
357
                                H1D_Fill(k, ix[0] , gSum[k] );
358
                                H1D_Fill(100+ k, ix[1] , gSum[k] );  
359
                        }      
360
                        PlotScanHistogramCB(0,0,EVENT_COMMIT,NULL, 0,0);
203 f9daq 361
                        if (!daq_on) break;
362
                }
363
                if (!daq_on) break;
364
        }
365
 
366
        return 0;
195 f9daq 367
}
368
 
203 f9daq 369
 
195 f9daq 370
int CVICALLBACK StartCB (int panel, int control, int event,
203 f9daq 371
                                                                                                 void *callbackData, int eventData1, int eventData2) {
372
        ThreadFunctionPtr mythread = NULL;
373
        switch (event) {
374
 
375
                case EVENT_COMMIT:
376
 
377
                        controlID=0;
378
                        if (panel == ph && control == PANEL_START) {
379
                                mythread = daq;
380
                                controlID= control;
381
                        }
382
                        if (panel == p2 && control == SCAN_SCAN)   mythread = scan;
383
                        if (mythread!=NULL) {
384
                                printf("New Thread panel=%d button=%d\n", panel, control);
385
 
386
                                // SetDimming(1);
387
 
388
                                daq_on=1;
389
                                CmtScheduleThreadPoolFunctionAdv (poolHandle, mythread, &controlID,
390
                                                                                                                                                                        DEFAULT_THREAD_PRIORITY,
391
                                                                                                                                                                        EndOfThread,
392
                                                                                                                                                                        EVENT_TP_THREAD_FUNCTION_END,
393
                                                                                                                                                                        NULL, RUN_IN_SCHEDULED_THREAD,
394
                                                                                                                                                                        &tfID);
395
                        }
396
                        break;
397
        }
398
        return 0;
195 f9daq 399
}
400
 
401
int CVICALLBACK StopCB (int panel, int control, int event,
203 f9daq 402
                                                                                                void *callbackData, int eventData1, int eventData2) {
403
        switch (event) {
404
                case EVENT_COMMIT:
405
                        daq_on=0;
406
                        break;
407
        }
408
        return 0;
195 f9daq 409
}
410
 
411
int CVICALLBACK ExitCB (int panel, int control, int event,
203 f9daq 412
                                                                                                void *callbackData, int eventData1, int eventData2) {
413
        switch (event) {
414
                case EVENT_COMMIT:
415
                        QuitUserInterface (0);
416
                        break;
417
        }
418
        return 0;
195 f9daq 419
}
203 f9daq 420
 
421
int CVICALLBACK MoveStageCB (int panel, int control, int event,
422
                                                                                                                 void *callbackData, int eventData1, int eventData2) {
423
        int axis=0, step=1000, direction=1, n;
424
        switch (event) {
425
                case EVENT_COMMIT:
426
 
427
                        if (panel == p2) {
428
                                switch (control) {
429
                                        case SCAN_BR :
430
                                                axis = 1;
431
                                                direction = 1;
432
                                                GetCtrlVal(p2, SCAN_STEPX, &step);
433
                                                break;
434
                                        case SCAN_BL :
435
                                                axis = 1;
436
                                                direction = -1;
437
                                                GetCtrlVal(p2, SCAN_STEPX, &step);
438
                                                break;
439
                                        case SCAN_BU :
440
                                                axis = 2;
441
                                                direction = 1;
442
                                                GetCtrlVal(p2, SCAN_STEPY, &step);
443
                                                break;
444
                                        case SCAN_BD :
445
                                                axis = 2;
446
                                                direction = -1;
447
                                                GetCtrlVal(p2, SCAN_STEPY, &step);
448
                                                break;
449
                                }
450
#ifdef MIKRO
451
                                MIKRO_MoveFor(axis,  direction*step );
452
                                MIKRO_GetPosition(axis,&n);
453
                                if (axis == 1) SetCtrlVal (p2, SCAN_XP, n);
454
                                if (axis == 2) SetCtrlVal (p2, SCAN_YP, n);
455
#endif // MIKRO                                 
456
                        }
457
 
458
                        break;
459
        }
460
        return 0;
461
}
462
 
463
 
464
 
465
 
466
 
467
int CVICALLBACK GoXCB (int panel, int control, int event,
468
                                                                                         void *callbackData, int eventData1, int eventData2) {
469
        int n2;
470
        switch (event) {
471
                case EVENT_COMMIT:
472
                        GetCtrlVal (p2, SCAN_XG, &n2);
473
#ifdef MIKRO
474
                        MIKRO_MoveTo(1,n2);
475
                        MIKRO_GetPosition(1,&n2);
476
 
477
#endif
478
                        SetCtrlVal (p2, SCAN_XP, n2);
479
                        break;
480
        }
481
        return 0;
482
}
483
 
484
int CVICALLBACK GoYCB (int panel, int control, int event,
485
                                                                                         void *callbackData, int eventData1, int eventData2) {
486
        int n2;
487
        switch (event) {
488
                case EVENT_COMMIT:
489
                        GetCtrlVal (p2, SCAN_YG, &n2);
490
#ifdef MIKRO
491
                        MIKRO_MoveTo(2,n2);
492
                        MIKRO_GetPosition(2,&n2);
493
 
494
#endif
495
                        SetCtrlVal (p2, SCAN_YP, n2);
496
                        break;
497
        }
498
        return 0;
499
}
500
 
501
int CVICALLBACK GetCurrentPositionCB (int panel, int control, int event,
502
                                                                                                                                        void *callbackData, int eventData1, int eventData2) {
503
 
504
        int n[2];      
505
        switch (event) {
506
                case EVENT_COMMIT:
507
#ifdef MIKRO                    
508
                        MIKRO_GetPosition(1,&n[0]);
509
                        SetCtrlVal (p2, SCAN_XP, n[0]);
510
                        MIKRO_GetPosition(2,&n[1]);
511
                        SetCtrlVal (p2, SCAN_YP, n[1]);
512
#endif                  
513
                        break;
514
        }
515
        return 0;
516
}
517
 
518
int CVICALLBACK HomeCB (int panel, int control, int event,
519
                                                                                                void *callbackData, int eventData1, int eventData2) {
520
        switch (event) {
521
                case EVENT_COMMIT:
522
#ifdef MIKRO                    
523
                        MIKRO_ReferenceMove(1);
524
                        MIKRO_ReferenceMove(2);
525
                        GetCurrentPositionCB(panel, control, event, NULL, 0, 0);
526
#endif                  
527
                        break;
528
        }
529
        return 0;
530
}
217 f9daq 531
 
532
int CVICALLBACK PlotScanHistogramCB (int panel, int control, int event,
533
                                                                                                                                                 void *callbackData, int eventData1, int eventData2) {
534
        int hid=0;
535
        int nx;
536
        int ny;
537
        switch (event) {
538
                case EVENT_COMMIT:
539
                        GetCtrlVal (p2, SCAN_CHANNEL, &hid);
540
                        GetCtrlVal(p2, SCAN_NSTEPSX, &nx);
541
            GetCtrlVal(p2, SCAN_NSTEPSY, &ny);
542
                        if (nx>1 && ny>1) {
543
                                H2D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
544
                                printf("redraw 2d\n");
545
                        } else {
546
                                if (nx>1)  {
547
                                        printf("redraw 1d x\n");
548
                                        H1D_Draw(hid,p2,SCAN_GRAPH,&scanplothandle);
549
                                }      
550
                          if (ny>1) {
551
                                        printf("redraw 1d x\n");
552
                                        H1D_Draw(100+hid,p2,SCAN_GRAPH,&scanplothandle);
553
                                }
554
                        }      
555
                        break;
556
        }
557
        return 0;
558
}
559
 
560
int CVICALLBACK OpenGuiCB (int panel, int control, int event,
561
                                                                                                         void *callbackData, int eventData1, int eventData2) {
562
        switch (event) {
563
                case EVENT_COMMIT:
564
                        DisplayPanel (ph);
565
            DisplayPanel (p2);
566
                        break;
567
        }
568
        return 0;
569
}