#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;
}