Subversion Repositories f9daq

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