fesimdaq.cxx File Reference

#include <vector>
#include <stdio.h>
#include <algorithm>
#include <stdlib.h>
#include "midas.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)
	Polling routine for events.
INT	interrupt_configure (INT cmd, INT source, POINTER_T adr)
	Interrupt configuration (not implemented)
double	sampleNormal ()

Variables
char *	frontend_name = "fesimdaq"
	The frontend name (client name) as seen by other MIDAS clients.
char *	frontend_file_name = __FILE__
	The frontend file name, don't change it.
BOOL	frontend_call_loop = TRUE
	frontend_loop is called periodically if this variable is TRUE
INT	display_period = 0
	a frontend status page is displayed with this frequency in ms
INT	max_event_size = 3 * 1024 * 1024
	maximum event size produced by this frontend
INT	max_event_size_frag = 2 * 1024 * 1024
	maximum event size for fragmented events (EQ_FRAGMENTED)
INT	event_buffer_size = 20 * 1000000
	buffer size to hold events
void **	info
char	strin [256]
HNDLE	hDB
	main ODB handle
HNDLE	hSet
EQUIPMENT	equipment []
	Main structure for midas 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 145 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 186 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 196 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 178 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 157 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 217 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
	)

Interrupt configuration (not implemented)

Routine for interrupt configuration if equipment is set in EQ_INTERRUPT mode. Not implemented right now, returns SUCCESS.

Parameters

[in]	cmd	Command for interrupt events (see midas.h)
[in]	source	Equipment index number
[in]	adr	Interrupt routine (see mfe.c)

Returns

Midas status code

Definition at line 257 of file fesimdaq.cxx.

           { INT interrupt_configure(INT cmd, INT source, POINTER_T adr)
{
   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 call 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 203 of file fesimdaq.cxx.

{
   return SUCCESS;
}

poll_event()

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

Polling routine for events.

Parameters

[in]	source	Event source (LAM/IRQ)
[in]	count	Loop count for event polling timeout
[in]	test	flag used to time the polling

Returns

1 if event is available, 0 if done polling (no event). If test equals TRUE, don't return.

Definition at line 235 of file fesimdaq.cxx.

           { INT poll_event(INT source, INT count, BOOL test)
/* Polling routine for events. Returns TRUE if event
   is available. If test equals TRUE, don't return. The test
   flag is used to time the polling */
{
   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:

read_scaler_event()

INT read_scaler_event	(	char *	pevent,
		INT	off
	)

Scaler event readout.

Placeholder.

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

Here is the call graph for this function:

read_trigger_event()

INT read_trigger_event	(	char *	pevent,
		INT	off
	)

Definition at line 284 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 210 of file fesimdaq.cxx.

{
   return SUCCESS;
}

sampleNormal()

double sampleNormal

(

)

Definition at line 274 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:

SIMDAQSETTINGS_STR()

SIMDAQSETTINGS_STR

(

simdaqsettings_str

)

wait_end_cycle()

BOOL wait_end_cycle	(	int	transition,
		BOOL	first
	)

Definition at line 130 of file fesimdaq.cxx.

 {
   if(first){
     std::cout << "Starting deferred transition" << std::endl;
     nremaining = 5;
   }
 
   if (finished_readout) {
     return TRUE;
   }
   
   return FALSE;
 }

Here is the call graph for this function:

Variable Documentation

display_period

INT display_period = 0

a frontend status page is displayed with this frequency in ms

Definition at line 35 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,      
    },
 
   {""}
}

Main structure for midas equipment.

Definition at line 72 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

buffer size to hold events

Definition at line 44 of file fesimdaq.cxx.

finished_readout

BOOL finished_readout =FALSE

Definition at line 126 of file fesimdaq.cxx.

frontend_call_loop

BOOL frontend_call_loop = TRUE

frontend_loop is called periodically if this variable is TRUE

Definition at line 32 of file fesimdaq.cxx.

frontend_file_name

char* frontend_file_name = __FILE__

The frontend file name, don't change it.

Definition at line 29 of file fesimdaq.cxx.

frontend_name

char* frontend_name = "fesimdaq"

The frontend name (client name) as seen by other MIDAS clients.

Definition at line 27 of file fesimdaq.cxx.

hDB

HNDLE hDB

main ODB handle

Definition at line 47 of file fesimdaq.cxx.

hSet

HNDLE hSet

Definition at line 47 of file fesimdaq.cxx.

info

void** info

Definition at line 45 of file fesimdaq.cxx.

max_event_size

INT max_event_size = 3 * 1024 * 1024

maximum event size produced by this frontend

Definition at line 38 of file fesimdaq.cxx.

max_event_size_frag

INT max_event_size_frag = 2 * 1024 * 1024

maximum event size for fragmented events (EQ_FRAGMENTED)

Definition at line 41 of file fesimdaq.cxx.

nremaining

int nremaining = -1

Definition at line 128 of file fesimdaq.cxx.

strin

char strin[256]

Definition at line 46 of file fesimdaq.cxx.