fesimdaq.cxx File Reference

#include <vector>
#include <stdio.h>
#include <algorithm>
#include <stdlib.h>
#include "midas.h"
#include "mfe.h"
#include <stdint.h>
#include <iostream>
#include <sstream>
#include <unistd.h>
#include "math.h"
#include <sys/time.h>

Include dependency graph for fesimdaq.cxx:

Go to the source code of this file.

Macros
#define	SIMDAQSETTINGS_STR(_name)

Functions
INT	frontend_init ()
	Frontend initialization.
INT	frontend_exit ()
	Frontend exit.
INT	begin_of_run (INT run_number, char *error)
	Begin of Run.
INT	end_of_run (INT run_number, char *error)
	End of Run.
INT	pause_run (INT run_number, char *error)
	Pause Run.
INT	resume_run (INT run_number, char *error)
	Resume Run.
INT	frontend_loop ()
	Frontend loop.
INT	read_trigger_event (char *pevent, INT off)
INT	read_scaler_event (char *pevent, INT off)
	Scaler event readout.
BOOL	wait_end_cycle (int transition, BOOL first)
	SIMDAQSETTINGS_STR (simdaqsettings_str)
INT	poll_event (INT source, INT count, BOOL test)
INT	interrupt_configure (INT cmd, INT source, POINTER_T adr)
double	sampleNormal ()

Variables
const char *	frontend_name = "fesimdaq"
const char *	frontend_file_name = __FILE__
BOOL	frontend_call_loop = TRUE
INT	display_period = 0
INT	max_event_size = 3 * 1024 * 1024
INT	max_event_size_frag = 2 * 1024 * 1024
INT	event_buffer_size = 20 * 1000000
void **	info
char	strin [256]
HNDLE	hSet
EQUIPMENT	equipment []
BOOL	finished_readout =FALSE
int	nremaining = -1

Macro Definition Documentation

SIMDAQSETTINGS_STR

#define SIMDAQSETTINGS_STR

(

_name

)

Value:

const char *_name[] = {\
"[.]",\
"Descrip = STRING : [256] ",\
"Parameter AAA = INT : 0",\
"Parameter BB2 = INT : 0",       \
"Parameter CC22 = INT : 0",\
"",\
NULL }

Definition at line 140 of file fesimdaq.cxx.

                                                        {\
"[.]",\
"Descrip = STRING : [256] ",\
"Parameter AAA = INT : 0",\
"Parameter BB2 = INT : 0",       \
"Parameter CC22 = INT : 0",\
"",\
NULL }

Function Documentation

begin_of_run()

INT begin_of_run	(	INT	run_number,
		char *	error
	)

Begin of Run.

Called every run start transition. Set equipment status in ODB, start acquisition on the modules.

Parameters

[in]	run_number	Number of the run being started
[out]	error	Can be used to write a message string to midas.log

Definition at line 181 of file fesimdaq.cxx.

{
 
  nremaining = -1;
  finished_readout=FALSE;
   return SUCCESS;
}

end_of_run()

INT end_of_run	(	INT	run_number,
		char *	error
	)

End of Run.

Called every stop run transition. Set equipment status in ODB, stop acquisition on the modules.

Parameters

[in]	run_number	Number of the run being ended
[out]	error	Can be used to write a message string to midas.log

Definition at line 191 of file fesimdaq.cxx.

{
   return SUCCESS;
}

frontend_exit()

INT frontend_exit

(

void

)

Frontend exit.

Runs at frontend shutdown. Disconnect hardware and set equipment status in ODB

Returns

Midas status code

Definition at line 173 of file fesimdaq.cxx.

{
  printf("Exiting fedcrc!\n");
   return SUCCESS;
}

Here is the call graph for this function:

frontend_init()

INT frontend_init

(

void

)

Frontend initialization.

Runs once at application startup. We initialize the hardware and optical interfaces and set the equipment status in ODB. We also lock the frontend to once physical cpu core.

Returns

Midas status code

Definition at line 152 of file fesimdaq.cxx.

{
 
  
  int status = db_check_record(hDB, 0, "/Equipment/SIMDAQ/Settings", strcomb1(simdaqsettings_str).c_str(), TRUE);
  printf("Status %i\n",status);
  if (status == DB_STRUCT_MISMATCH) {
    cm_msg(MERROR, "init_simdaqsettings", "Aborting on mismatching /Equipment/SIMDAQ/Settings");
    cm_disconnect_experiment();
    abort();
  }
 
  cm_register_deferred_transition(TR_STOP, wait_end_cycle);
   return SUCCESS;
}

Here is the call graph for this function:

frontend_loop()

INT frontend_loop

(

void

)

Frontend loop.

If frontend_call_loop is true, this routine gets called when the frontend is idle or once between every event.

Returns

Midas status code

Definition at line 212 of file fesimdaq.cxx.

{
   /* if frontend_call_loop is true, this routine gets called when
      the frontend is idle or once between every event */
  usleep(10000);
  return SUCCESS;
}

Here is the call graph for this function:

interrupt_configure()

INT interrupt_configure	(	INT	cmd,
		INT	source,
		POINTER_T	adr
	)

Definition at line 251 of file fesimdaq.cxx.

{
   switch (cmd) {
   case CMD_INTERRUPT_ENABLE:
      break;
   case CMD_INTERRUPT_DISABLE:
      break;
   case CMD_INTERRUPT_ATTACH:
      break;
   case CMD_INTERRUPT_DETACH:
      break;
   }
   return SUCCESS;
}

Here is the caller graph for this function:

pause_run()

INT pause_run	(	INT	run_number,
		char *	error
	)

Pause Run.

Called every pause run transition.

Parameters

[in]	run_number	Number of the run being ended
[out]	error	Can be used to write a message string to midas.log

Returns

Midas status code

Definition at line 198 of file fesimdaq.cxx.

{
   return SUCCESS;
}

poll_event()

INT poll_event	(	INT	source,
		INT	count,
		BOOL	test
	)

Definition at line 230 of file fesimdaq.cxx.

{
   int i;
 
   for (i = 0; i < count; i++) {
//      cam_lam_read(LAM_SOURCE_CRATE(source), &lam);
 
//      if (lam & LAM_SOURCE_STATION(source))
         if (!test)
            return 1;
   }
 
   usleep(1000);
   return 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

read_scaler_event()

INT read_scaler_event	(	char *	pevent,
		INT	off
	)

Scaler event readout.

Scaler event readout routine. Not used in this example.

Parameters

[in]	pevent	Pointer to event buffer
[in]	off	Caller info (unused here), see mfe.c

Returns

Size of the event

Definition at line 295 of file ebfe.cxx.

{
   DWORD *pdata, a;
 
   /* init bank structure */
   bk_init(pevent);
 
   /* create SCLR bank */
   bk_create(pevent, Sclx, TID_DWORD, (void**)&pdata);
 
   /* read scaler bank */
   for (a = 0; a < N_SCLR; a++)
      cam24i(CRATE, SLOT_SCLR, a, 0, pdata++);
 
   bk_close(pevent, pdata);
 
   return bk_size(pevent);
}

read_trigger_event()

INT read_trigger_event	(	char *	pevent,
		INT	off
	)

Definition at line 277 of file fesimdaq.cxx.

{
 
  /* init bank structure */
  bk_init32(pevent);
  
  uint32_t *pdata32;
  /* create structured ADC0 bank of double words (i.e. 4-byte words)  */
  bk_create(pevent, "ADC0", TID_DWORD, (void **)&pdata32);
 
  // Add a header, with number of words in event.
  // Use the top two bits to indicate different control words.
  // 11 -> 0xcXXXXXXX  : overall header
  // 01 -> 0x4XXXXXXX  : trigger header
  // 00 -> 0x0XXXXXXX  : channel header
  // 10 -> 0x8XXXXXXX  : trailer
  
  // Write the bank header to the bank, containing the number of triggers
  *pdata32++ = 0xfa000200;
 
  int sample = floor(sampleNormal()*80)+1200;
  if((rand() % 100) > 80) sample += 200;
  *pdata32++ = 0xf800406b + sample;
 
  int sample2 = floor(sampleNormal()*120)+1000;
  *pdata32++ = 0xf801405e + sample2;
 
  int sample3 = floor(sampleNormal()*100)+1000;
  *pdata32++ = 0xf802405e + sample3;
 
  int sample4 = floor(sampleNormal()*20)+950;
  *pdata32++ = 0xf803405e + sample4;
 
  int sample5 = floor(sampleNormal()*50)+1300;
  *pdata32++ = 0xf804405e + sample5;
 
  int sample6 = floor(sampleNormal()*200)+2000;
  *pdata32++ = 0xf805405e + sample6;
 
 
 
  *pdata32++ = 0xfcfd57e8;
 
  
  int size2 = bk_close(pevent, pdata32);    
  
 
  uint32_t *pdata322;
  bk_create(pevent, "ADC1", TID_DWORD, (void **)&pdata322);
 
 
  //for(int i = 0; i < 100000; i++)
    //*pdata322++ = i;
  
  size2 = bk_close(pevent, pdata322);    
 
  uint32_t *pdata720;
  bk_create(pevent, "W200", TID_DWORD, (void **)&pdata720);
 
 
  *pdata720++ = 0xa00001f8;
  *pdata720++ = 0xfff703;
  *pdata720++ = 0x0;
  *pdata720++ = 0x13132f;
  
  int pulse_position = floor(sampleNormal()*80) + 200;
  for(int j = 0; j < 2; j++){
 
    if (j==1) pulse_position+=  5;
    int pulse_height = floor(sampleNormal()*30) +200;  
    pulse_height = pulse_height * (1.0 - ((double)pulse_position)/600.0);
 
    for(int i = 0; i < 504; i++){
      uint32_t word = 0;
      uint32_t sample = 400 + floor(sampleNormal()*2);
      if(i == pulse_position) sample += 0.2 * pulse_height;
      if(i == pulse_position+1) sample += 0.8 * pulse_height;
      if(i == pulse_position+2) sample += 1.0 * pulse_height;
      if(i == pulse_position+3) sample += 0.9 * pulse_height;
      if(i == pulse_position+4) sample += 0.5 * pulse_height;
      if(i == pulse_position+5) sample += 0.3 * pulse_height;
      if(i == pulse_position+6) sample += 0.1 * pulse_height;    
      //std::cout << i << " " << pulse_position << " " << pulse_height << " " << sample << std::endl;
      word = sample;
      i++;
      sample = 400 + floor(sampleNormal()*2);
      if(i == pulse_position) sample += 0.2 * pulse_height;
      if(i == pulse_position+1) sample += 0.8 * pulse_height;
      if(i == pulse_position+2) sample += 1.0 * pulse_height;
      if(i == pulse_position+3) sample += 0.9 * pulse_height;
      if(i == pulse_position+4) sample += 0.5 * pulse_height;
      if(i == pulse_position+5) sample += 0.3 * pulse_height;
      if(i == pulse_position+6) sample += 0.1 * pulse_height;    
      //std::cout << j << " " << i << " " << pulse_position << " " << pulse_height << " " << sample << std::endl;
      word += (sample << 16);
      *pdata720++ = word;
    }  
  }
  
  
 
  size2 = bk_close(pevent, pdata720);    
 
 
 
  struct timeval start; //,end;
  gettimeofday(&start,NULL);
  if(0)printf ("simdaq request: %f\n",start.tv_sec
      + 0.000001*start.tv_usec); 
  // close the bank
 
  // handle the deferred transition
  if(nremaining > 0){ std::cout << "deferring " << nremaining << std::endl; nremaining--;
  };
 
  if(nremaining == 0){
    finished_readout=TRUE;
    std::cout << "Finished readout... " << std::endl;
  }
 
    /*  
  printf("Start sleep\n");
  bool done = false;
  double time;
  while(!done){
    sleep(5);
    gettimeofday(&end,NULL);
    time = end.tv_sec - start.tv_sec + (end.tv_usec - start.tv_usec)*0.000001;
    if(time > 10) done = true;
  }
  
      //}
  printf("End sleep. time %f\n", time);
  */
  
  usleep(10000);
  return bk_size(pevent);
}

Here is the call graph for this function:

resume_run()

INT resume_run	(	INT	run_number,
		char *	error
	)

Resume Run.

Called every resume run transition.

Parameters

[in]	run_number	Number of the run being ended
[out]	error	Can be used to write a message string to midas.log

Returns

Midas status code

Definition at line 205 of file fesimdaq.cxx.

{
   return SUCCESS;
}

sampleNormal()

double sampleNormal

(

)

Definition at line 267 of file fesimdaq.cxx.

                      {
    double u = ((double) rand() / (RAND_MAX)) * 2 - 1;
    double v = ((double) rand() / (RAND_MAX)) * 2 - 1;
    double r = u * u + v * v;
    if (r == 0 || r > 1) return sampleNormal();
    double c = sqrt(-2 * log(r) / r);
    return u * c;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SIMDAQSETTINGS_STR()

SIMDAQSETTINGS_STR

(

simdaqsettings_str

)

wait_end_cycle()

BOOL wait_end_cycle	(	int	transition,
		BOOL	first
	)

Definition at line 123 of file fesimdaq.cxx.

 {
   //printf("wait_end_cycle: transition %d, first %d\n", transition, first);
   
   if(first){
     std::cout << "Starting deferred transition" << std::endl;
     nremaining = 10;
   }
 
   if (finished_readout) {
     return TRUE;
   }
 
   return FALSE;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Variable Documentation

display_period

INT display_period = 0

Definition at line 31 of file fesimdaq.cxx.

equipment

EQUIPMENT equipment[]

Initial value:

= {
 
   {"SIMDAQ",               
    {1, 0,                   
     "SYSTEM",               
 
 
 
     EQ_PERIODIC,              
 
     LAM_SOURCE(0, 0xFFFFFF),        
     "MIDAS",                
     TRUE,                   
     RO_RUNNING,           
     1000,                    
     0,                      
     0,                      
     0,                      
     "", "", "",},
    read_trigger_event,      
    },
 
   {""}
}

Definition at line 69 of file fesimdaq.cxx.

                        {
 
   {"SIMDAQ",               /* equipment name */
    {1, 0,                   /* event ID, trigger mask */
     "SYSTEM",               /* event buffer */
#ifdef USE_INT
     EQ_INTERRUPT,           /* equipment type */
#else
     EQ_PERIODIC,              /* equipment type */
#endif
     LAM_SOURCE(0, 0xFFFFFF),        /* event source crate 0, all stations */
     "MIDAS",                /* format */
     TRUE,                   /* enabled */
     RO_RUNNING,           /* read only when running */
     1000,                    /* poll for 1000ms */
     0,                      /* stop run after this event limit */
     0,                      /* number of sub events */
     0,                      /* don't log history */
     "", "", "",},
    read_trigger_event,      /* readout routine */
    },
 
   {""}
};

event_buffer_size

INT event_buffer_size = 20 * 1000000

Definition at line 40 of file fesimdaq.cxx.

finished_readout

BOOL finished_readout =FALSE

Definition at line 119 of file fesimdaq.cxx.

frontend_call_loop

BOOL frontend_call_loop = TRUE

Definition at line 28 of file fesimdaq.cxx.

frontend_file_name

const char* frontend_file_name = __FILE__

Definition at line 25 of file fesimdaq.cxx.

frontend_name

const char* frontend_name = "fesimdaq"

Definition at line 23 of file fesimdaq.cxx.

hSet

HNDLE hSet

Definition at line 44 of file fesimdaq.cxx.

info

void** info

Definition at line 41 of file fesimdaq.cxx.

max_event_size

INT max_event_size = 3 * 1024 * 1024

Definition at line 34 of file fesimdaq.cxx.

max_event_size_frag

INT max_event_size_frag = 2 * 1024 * 1024

Definition at line 37 of file fesimdaq.cxx.

nremaining

int nremaining = -1

Definition at line 121 of file fesimdaq.cxx.

strin

char strin[256]

Definition at line 42 of file fesimdaq.cxx.