| 0,0 → 1,492 |
|---|
| #define USE_DAQ |
| |
| // Izberi ustrezni interface v meniju projektnega okna |
| // Options->CompilerDefines (dodaj /DSISVME ali /DWIENVME) |
| #ifdef USE_DAQ |
| # ifdef SISVME |
| # include "sisvme_dll.h" |
| # endif |
| # ifdef WIENVME |
| # include "wienvme_dll.h" |
| # endif |
| # ifdef WUSBVME |
| # include "wusbvme_dll.h" |
| # endif |
| # include "CAENV673A.h" |
| #endif |
| |
| #include <utility.h> |
| #include <userint.h> |
| #include <cvirte.h> /* Needed if linking in external compiler; harmless otherwise */ |
| #include "rich_ui.h" |
| |
| #define V673A_BASE_A24 0x130000 |
| #define USBVME_NAME "VM0120" |
| |
| #define MAPFILE "icfa2013.txt" |
| |
| #define MAXDATA 2000 |
| unsigned long dataa[3078]; |
| unsigned long data0[MAXDATA]; |
| unsigned long datac[1000]; |
| |
| |
| int p1h, p2h, p3h ; |
| int p1h,rph,rID; |
| |
| |
| int pmtmaparr[8][8] ; |
| int pmtmap_i[64] ; |
| int pmtmap_j[64] ; |
| int data[64][64] ; |
| int nndata[64] ; |
| unsigned short acum_allch[65536] ; |
| unsigned short acum_snglch[64][512] ; |
| //int offset=512 ; |
| int offset=0 ; |
| int min_tdc, max_tdc ; |
| int min_tdcR, max_tdcR ; |
| |
| int single_event ; |
| |
| unsigned long pmtbitmap[2] ; |
| unsigned long pmtbitmap_win[2] ; |
| |
| int minhit ; |
| |
| //***************************************************************************** |
| |
| int read_pmtmap ( ) { |
| FILE *fp=fopen(MAPFILE,"r"); |
| int ch, i, j ; |
| |
| while ( EOF != fscanf(fp,"%d %d %d\n",&ch,&i,&j) ) { |
| if ( i>7 || j>7 ) { |
| printf ( "Index out of range!!!! (%d,%d,%d)\n",ch,i,j ) ; |
| continue ; |
| } |
| pmtmaparr[i][j]=ch ; |
| pmtmap_i[ch] = i ; |
| pmtmap_j[ch] = j ; |
| } |
| |
| return 0 ; |
| } |
| |
| //***************************************************************************** |
| |
| int draw_histos ( ) { |
| int i, ch ; |
| int harr[4]={HISTOS_H_1,HISTOS_H_2,HISTOS_H_3,HISTOS_H_4} ; |
| |
| for ( ch=0; ch<4; ch++ ) { |
| DeleteGraphPlot (p2h, harr[ch], -1, VAL_IMMEDIATE_DRAW); |
| |
| for ( i=0 ; i<nndata[ch] ; i++ ) { |
| if ( (data[ch][i]>>16)&(0x1) ) { |
| //printf ( "falling\n" ) ; |
| PlotPoint (p2h, harr[ch], data[ch][i]&0xFFFF, 1.0, VAL_SOLID_SQUARE, VAL_RED); |
| PlotLine (p2h, harr[ch], data[ch][i]&0xFFFF, 1.0, data[ch][i]&0xFFFF, 0.0, VAL_RED); |
| } else { |
| //printf ( "rising %d %X\n", ch, data[ch][i]&0xFFFF ) ; |
| PlotPoint (p2h, harr[ch], data[ch][i]&0xFFFF, 1.0, VAL_EMPTY_SQUARE, VAL_BLUE); |
| PlotLine (p2h, harr[ch], data[ch][i]&0xFFFF, 1.0, data[ch][i]&0xFFFF, 0.0, VAL_BLUE); |
| } |
| } |
| |
| } |
| |
| |
| return 0 ; |
| } |
| |
| //***************************************************************************** |
| |
| int draw_acum_histos() { |
| int hsts[64]={P3_ACUM_1,P3_ACUM_2,P3_ACUM_3,P3_ACUM_4,P3_ACUM_5,P3_ACUM_6,P3_ACUM_7,P3_ACUM_8,P3_ACUM_9,P3_ACUM_10, |
| P3_ACUM_11,P3_ACUM_12,P3_ACUM_13,P3_ACUM_14,P3_ACUM_15,P3_ACUM_16,P3_ACUM_17,P3_ACUM_18,P3_ACUM_19,P3_ACUM_20, |
| P3_ACUM_21,P3_ACUM_22,P3_ACUM_23,P3_ACUM_24,P3_ACUM_25,P3_ACUM_26,P3_ACUM_27,P3_ACUM_28,P3_ACUM_29,P3_ACUM_30, |
| P3_ACUM_31,P3_ACUM_32,P3_ACUM_33,P3_ACUM_34,P3_ACUM_35,P3_ACUM_36,P3_ACUM_37,P3_ACUM_38,P3_ACUM_39,P3_ACUM_40, |
| P3_ACUM_41,P3_ACUM_42,P3_ACUM_43,P3_ACUM_44,P3_ACUM_45,P3_ACUM_46,P3_ACUM_47,P3_ACUM_48,P3_ACUM_49,P3_ACUM_50, |
| P3_ACUM_51,P3_ACUM_52,P3_ACUM_53,P3_ACUM_54,P3_ACUM_55,P3_ACUM_56,P3_ACUM_57,P3_ACUM_58,P3_ACUM_59,P3_ACUM_60, |
| P3_ACUM_61,P3_ACUM_62,P3_ACUM_63,P3_ACUM_64} ; |
| int i; |
| double grmin,grmax; |
| |
| GetCtrlVal (p1h, P1_TDC_MAX, &max_tdc ) ; |
| GetCtrlVal (p1h, P1_TDC_MIN, &min_tdc ) ; |
| GetCtrlVal (p1h, P1_TDC_MAX_2, &max_tdcR ) ; |
| GetCtrlVal (p1h, P1_TDC_MIN_2, &min_tdcR ) ; |
| //SetAxisScalingMode (p1h, P1_ACUM_WIN, VAL_XAXIS, VAL_MANUAL, min_tdc, max_tdc); |
| |
| DeleteGraphPlot (p1h, P1_ACUM, -1, VAL_IMMEDIATE_DRAW); |
| DeleteGraphPlot (p1h, P1_ACUM_WIN, -1, VAL_IMMEDIATE_DRAW); |
| |
| PlotY (p1h, P1_ACUM, acum_allch, /*65536*/ 65536, VAL_UNSIGNED_SHORT_INTEGER, |
| VAL_VERTICAL_BAR, VAL_NO_POINT, VAL_SOLID, 1, VAL_DK_BLUE); |
| PlotY (p1h, P1_ACUM_WIN, acum_allch, /*65536*/ 65536, VAL_UNSIGNED_SHORT_INTEGER, |
| VAL_VERTICAL_BAR, VAL_NO_POINT, VAL_SOLID, 1, VAL_DK_BLUE); |
| |
| GetAxisScalingMode (p1h, P1_ACUM_WIN, VAL_LEFT_YAXIS, 0, &grmin, &grmax ); |
| |
| PlotLine (p1h, P1_ACUM_WIN, min_tdcR, 0, min_tdcR, grmax, VAL_RED); |
| PlotLine (p1h, P1_ACUM_WIN, max_tdcR, 0, max_tdcR, grmax, VAL_RED); |
| PlotLine (p1h, P1_ACUM_WIN, min_tdc, 0, min_tdc, grmax, VAL_GREEN); |
| PlotLine (p1h, P1_ACUM_WIN, max_tdc, 0, max_tdc, grmax, VAL_GREEN); |
| |
| for ( i=0; i<64; i++ ) { |
| DeleteGraphPlot (p3h, hsts[i], -1, VAL_IMMEDIATE_DRAW); |
| |
| PlotY (p3h, hsts[i], acum_snglch[i], 256, VAL_UNSIGNED_SHORT_INTEGER, |
| VAL_VERTICAL_BAR, VAL_NO_POINT, VAL_SOLID, 1, VAL_DK_BLUE); |
| |
| } |
| |
| return 0 ; |
| } |
| |
| //***************************************************************************** |
| |
| int draw_evnt() { |
| int i,j,color,frame; |
| int iijj[8]={1,3, 7,9, 13,15, 19,21} ; |
| // int iijj[8]={1,3, 7,9, 13,15, 19,21} ; |
| |
| DeleteGraphPlot (p1h, P1_EV_DI, -1, VAL_IMMEDIATE_DRAW); |
| for ( i=1; i<20; i+=6 ) { |
| for ( j=1; j<20; j+=6 ) { |
| PlotRectangle (p1h,P1_EV_DI,i,j,i+4,j+4,VAL_DK_GRAY,VAL_TRANSPARENT); |
| } |
| } |
| for ( i=0; i<64; i++ ) { |
| if ( (pmtbitmap_win[i>>5]>>(i&0x1F))&0x1 ) { |
| color = VAL_DK_GREEN; |
| frame = VAL_DK_GREEN; |
| } else if ( (pmtbitmap[i>>5]>>(i&0x1F))&0x1 ) { |
| color = VAL_TRANSPARENT; |
| frame = VAL_RED; |
| } else { |
| color = -1 ; |
| } |
| |
| if ( color!=-1 ) { |
| PlotRectangle (p1h, P1_EV_DI, |
| iijj[pmtmap_i[i]]+0.1,iijj[pmtmap_j[i]]+0.1, |
| iijj[pmtmap_i[i]]+1.9,iijj[pmtmap_j[i]]+1.9, |
| frame,color); |
| } |
| } |
| return 0 ; |
| } |
| |
| //***************************************************************************** |
| |
| int num_gr_photons() { |
| |
| int rc=0; |
| long bitmap[2]; |
| int i; |
| |
| bitmap[0]=pmtbitmap_win[0]; |
| bitmap[1]=pmtbitmap_win[1]; |
| for (i=0;i<32;i++) { |
| if (bitmap[0]&0x1) rc++; |
| if (bitmap[1]&0x1) rc++; |
| bitmap[0]>>=1; |
| bitmap[1]>>=1; |
| } |
| printf("%08X\t%08X\n",pmtbitmap_win[0],pmtbitmap_win[1]); |
| printf("num_gr_photons %d\n",rc); |
| return rc; |
| } |
| |
| //***************************************************************************** |
| |
| int initvar(unsigned long xxx[],int n) { |
| for (;--n>=0;xxx[n]=0); |
| return 0; |
| } |
| |
| //***************************************************************************** |
| |
| int daq(){ |
| |
| static int cdaq=0; |
| int dum; |
| int error,i,j,fcount,fcountold=-1,nauto,fauto; |
| int nfirst[3],ntrig[3],ndata[3],ndatac,neve,ntrold; |
| unsigned short *j16=(unsigned short *) &j; |
| |
| unsigned long ext_c; |
| |
| char filename[256],path[256], fname[256]; |
| int maxeve,print,neve_sel=0,nueve=0; |
| |
| unsigned int bsr ; |
| unsigned short dataxxx ; |
| unsigned short edge_type, overflow, channel, ev_dat, tdc_num, last_dat, nval_dat; |
| unsigned short gerror ; |
| int hits=0 ; |
| //int last[2] ; |
| int last_ev_no = -1 ; |
| |
| |
| if (cdaq==1) return -1; |
| else cdaq=1; |
| |
| // initialize VME interface |
| #ifdef USE_DAQ |
| error = V673A_map (0,V673A_BASE_A24, 1); |
| error = V673A_init(0); |
| V673A_setreg(0,0x1a,0x003e); |
| V673A_setreg(0,0x1e,0x003e); |
| error = V673A_clallb(0); |
| #endif |
| |
| |
| SetCtrlVal (p1h, P1_LED_DAQ,1); |
| SetCtrlVal (p1h, P1_TRIG, 0 ) ; |
| GetCtrlVal (p1h, P1_DEBUG, &print ) ; |
| GetCtrlVal (p1h, P1_NAUTO, &nauto ) ; |
| GetCtrlVal (p1h, P1_AUTOSAVE, &fauto ) ; |
| GetCtrlVal (p1h, P1_NEVNT, &maxeve ) ; |
| GetCtrlVal (p1h, P1_FNAME, fname ); |
| GetCtrlVal (p1h, P1_PATH, path ); |
| GetCtrlVal (p1h, P1_TDC_MAX, &max_tdc ) ; |
| GetCtrlVal (p1h, P1_TDC_MIN, &min_tdc ) ; |
| GetCtrlVal (p1h, P1_TDC_MAX_2, &max_tdcR ) ; |
| GetCtrlVal (p1h, P1_TDC_MIN_2, &min_tdcR ) ; |
| GetCtrlVal (p1h, P1_SINGLE, &single_event ) ; |
| //PlotLine (p1h, P1_ACUM_WIN, min_tdcR, 0, min_tdcR, 1000, VAL_RED); |
| //PlotLine (p1h, P1_ACUM_WIN, max_tdcR, 0, max_tdcR, 1000, VAL_RED); |
| //PlotLine (p1h, P1_ACUM_WIN, min_tdc, 0, min_tdc, 1000, VAL_GREEN); |
| //PlotLine (p1h, P1_ACUM_WIN, max_tdc, 0, max_tdc, 1000, VAL_GREEN); |
| //SetAxisScalingMode (p1h, P1_ACUM_WIN, VAL_XAXIS, VAL_MANUAL, min_tdc, max_tdc); |
| // start measurement |
| neve=0; |
| neve_sel=0; |
| ntrold=0; |
| rID=0; |
| |
| do { |
| /* |
| if (( neve%nauto==0 && fauto==1 ) || (neve==0 && fauto==1)){ |
| fcount=neve/nauto; |
| if ( fcountold < fcount ) { |
| sprintf(filename,"%s\\%s_file%d.dat",path,fname,fcount); |
| printf("%s\n",filename); |
| fcountold=fcount ; |
| } |
| } |
| */ |
| |
| if ( rID == P1_STOP ) break ; |
| #ifdef USE_DAQ |
| ntrig[0] = V673A_ntrig (0); |
| #endif |
| //printf ( "%d\t", ntrig[0] ) ; |
| //if (ntrig[0]!=-1) { |
| if (ntrig[0]>ntrold) { |
| ntrold=ntrig[0]-1; |
| initvar ( data0, 1024 ) ; |
| |
| #ifdef USE_DAQ |
| ndata[0] = V673A_getevnt (0,ntrold,data0,MAXDATA); |
| #endif |
| //printf ( "%d\t", ndata[0] ) ; |
| |
| for ( i=0;i<ndata[0];i++ ) { |
| dataa[i]=0xc3ffffff & data0[i]; |
| // printf ( "%d\t", data0[i] ) ; |
| } |
| ntrold=(ntrold+1)&0xffff; |
| |
| if ( ndata[0] >0 ) { |
| pmtbitmap[0] = 0 ; |
| pmtbitmap[1] = 0 ; |
| pmtbitmap_win[0] = 0 ; |
| pmtbitmap_win[1] = 0 ; |
| |
| for (i=0;i<ndata[0];i++) { |
| dataxxx = data0[i] & 0xFFFF ; |
| edge_type = (data0[i]>>16)&0x1 ; |
| overflow = (data0[i]>>17)&0x1 ; |
| tdc_num = (data0[i]>>25)&0x1 ; |
| channel = (data0[i]>>18)&0x1F|tdc_num<<5 ; |
| ev_dat = (data0[i]>>23)&0x1 ; |
| last_dat = (data0[i]>>30)&0x1 ; |
| nval_dat = (data0[i]>>31)&0x1 ; |
| if (ev_dat == 0) { // event data |
| hits++ ; |
| if ( nndata[channel] < 63 ) { |
| data[channel][nndata[channel]++] = dataxxx|edge_type<<16 ; |
| |
| } |
| if ((edge_type & 0x1)==0 ) { // rising edge |
| acum_allch[dataxxx]++ ; |
| if ((0<=(dataxxx-offset))&&((dataxxx-offset)<512) ) { // siroko okno |
| acum_snglch[channel][dataxxx-offset]++ ; |
| } |
| if ((min_tdcR<=(dataxxx-offset)) && ((dataxxx-offset)<max_tdcR)) { // rdece okno |
| if ( channel>31 ) |
| pmtbitmap[1] |= 1<<(channel-32) ; |
| else |
| pmtbitmap[0] |= 1<<(channel) ; |
| } |
| if ( (min_tdc<=(dataxxx-offset)) && ((dataxxx-offset)<max_tdc) ) { // zeleno okno |
| if ( channel>31 ) |
| pmtbitmap_win[1] |= 1<<(channel-32); |
| else |
| pmtbitmap_win[0] |= 1<<(channel); |
| } |
| } |
| } else if (last_dat==0) { //not last datum |
| //pmtbitmap[0] = 0; |
| //pmtbitmap[1] = 0; |
| last_ev_no = dataxxx; |
| } |
| } |
| |
| if (num_gr_photons() >= minhit) { |
| draw_evnt(); |
| //draw_histos(); |
| draw_acum_histos(); |
| neve_sel++; |
| } |
| |
| } |
| neve++; |
| } else { |
| error = V673A_clallb(0); |
| ntrold=0; |
| } |
| SetCtrlVal (p1h, P1_PROC, (neve*100)/maxeve ) ; |
| SetCtrlVal (p1h, P1_TRIG, neve ) ; |
| do { |
| nueve=GetUserEvent (0, &rph, &rID); |
| } while ( rID != P1_STOP && nueve>0 ); |
| if ( rID == P1_STOP ) break ; |
| } while ( (neve < maxeve) && (neve_sel<1 || !single_event) ) ; |
| |
| SetCtrlVal (p1h, P1_LED_DAQ,0); |
| cdaq=0; |
| return 0; |
| } |
| |
| //***************************************************************************** |
| //*************************** MAIN ******************************************** |
| //***************************************************************************** |
| |
| int __stdcall WinMain (HINSTANCE hInstance, HINSTANCE hPrevInstance, |
| LPSTR lpszCmdLine, int nCmdShow) |
| { |
| int error; |
| int i,j; |
| int nevnt,cnt; |
| int pmtmap[32][32] ; |
| int gray[16] = {0x000000, 0x111111, 0x222222, 0x333333, 0x444444, 0x555555, 0x666666, |
| 0x777777, 0x888888, 0x999999, 0xAAAAAA, 0xBBBBBB, 0xCCCCCC, 0xDDDDDD, |
| 0xEEEEEE, 0xFFFFFF } ; |
| double x, y, min_r, max_r, r1, r2 ; |
| |
| ColorMapEntry grCol[16] ; |
| |
| grCol[0].dataValue.valInt = 0 ; |
| grCol[0].color = gray[0] ; |
| grCol[1].dataValue.valInt = 100000 ; |
| grCol[1].color = gray[15] ; |
| |
| |
| read_pmtmap ( ) ; |
| |
| for ( i=0 ; i<32; i++ ) { |
| for ( j=0 ; j<32; j++ ) { |
| if ( ( (((i+2)/4)%2) != 0 ) && ( (((j+2)/4)%2) != 0 ) ) |
| pmtmap[i][j]=1000000 ; |
| else pmtmap[i][j]=0 ; |
| //pmtmap[i][j]=(i+1)*(i+1)*(i+1)*(j+1) ; |
| } |
| } |
| |
| p1h = LoadPanel (0, "rich_ui.uir", P1); |
| // p2h = LoadPanel (0, "laprich.uir", HISTOS); |
| p3h = LoadPanel (0, "rich_ui.uir", P3); |
| DisplayPanel (p1h); |
| // DisplayPanel (p2h); |
| DisplayPanel (p3h); |
| for ( i=1; i<20; i+=6 ) { |
| for ( j=1; j<20; j+=6 ) { |
| PlotRectangle (p1h,P1_EV_DI,i,j,i+4,j+4,VAL_DK_GRAY,VAL_TRANSPARENT); |
| } |
| } |
| |
| #ifdef USE_DAQ |
| # ifdef WUSBVME |
| VME_START(USBVME_NAME); |
| # else |
| VME_START(NULL); |
| # endif |
| #endif |
| |
| do { |
| GetUserEvent (1, &rph, &rID); |
| switch (rID){ |
| case P1_EV_PL: |
| // PlotScaledIntensity (p1h, P1_EV_DI, pmtmap, 32, 32, VAL_INTEGER, |
| // 1.0, 0.0, 1.0, 0.0, grCol, 0xff0000, 2, 1, 0); |
| draw_evnt(); |
| break; |
| case P1_READ_TDC: |
| daq(); //read_tdc ( ) ; |
| break; |
| case P1_DRAWH: |
| // draw_histos ( ) ; |
| draw_acum_histos ( ) ; |
| break; |
| case P1_SINGLE: |
| GetCtrlVal (p1h, P1_SINGLE, &single_event ) ; |
| break; |
| case P1_CONT_DAQ: |
| nevnt=0; |
| cnt=0; |
| minhit=0; |
| GetCtrlVal (p1h, P1_NEVNT, &nevnt ); |
| GetCtrlVal (p1h, P1_MINHIT, &minhit ); |
| // printf("%d %d\n", nevnt, minhit); |
| // while((read_tdc()<minhit)&&(cnt++<nevnt)) ; |
| daq(); |
| // draw_evnt(); |
| break; |
| case P1_EV_DI: |
| GetCtrlVal (p1h, P1_R_MIN, &min_r ) ; |
| GetCtrlVal (p1h, P1_R_MAX, &max_r ) ; |
| r1=min_r*2.0 ; |
| r2=max_r*2.0 ; |
| GetGraphCursor (p1h, P1_EV_DI, 1, &x, &y); |
| printf ( "x:%5.2f\ty:%5.2f\n", x, y ) ; |
| draw_evnt(); |
| PlotOval (p1h, P1_EV_DI, x-r1, y-r1, x+r1, y+r1, VAL_RED, VAL_TRANSPARENT); |
| PlotOval (p1h, P1_EV_DI, x-r2, y-r2, x+r2, y+r2, VAL_RED, VAL_TRANSPARENT); |
| break; |
| case P1_PRINT: |
| SetCtrlVal (p1h, P1_PL,1); |
| PrintCtrl (p1h, P1_EV_DI, "event.tif", 1, VAL_FULL_PANEL); |
| SetCtrlVal (p1h, P1_PL,0); |
| break; |
| case P1_RESET: |
| break; |
| // case P1_PRINT: |
| // break; |
| default: |
| break; |
| } |
| // if ( rID==P1_I_GATE ) init_gate ( ) ; |
| } while ( (rID!=P1_EXIT) || (rph!=p1h) ) ; |
| |
| // CloseHandle(nHandle); |
| #ifdef USE_DAQ |
| VME_STOP(); |
| #endif |
| // printf ( "XXXX\n" ) ; |
| return 0; |
| } |