void start() {
while (1) {
state_check();
make_events();
move();
my_send("stepForward");
currentTime += dt;
counter++;
}
}
bool XiaEventReader::hasMoreEvents() {
if(events_index<events_cnt) {
return 1;
}
reset();
bool found = 0;
if(bufferset) {
treat_buffer(lastbuffer);
found = 1;
}
bufferset = 1;
int len = DGF_H_NWORDS + 1;
int status;
while((status = fread(lastbuffer, sizeof(uint16_t), len, fbuffer))>0) {
found = 1;
// read the DGF_H_NWORDS and find the size of the buffer
if(status!=len) {
long pos = ftell(fbuffer);
fprintf(stderr, "warning: filesize too short for buffer header, position %ld, %d words required but only %d available\n", pos, len, status);
continue;
//exit(1);
}
uint16_t nwords = lastbuffer[DGF_H_NWORDS];
int datlen = nwords-len;
status = fread(&lastbuffer[len], sizeof(uint16_t), datlen, fbuffer);
if(status!=datlen) {
long pos = ftell(fbuffer);
fprintf(stderr, "warning: filesize too short for buffer data, position %ld, %d words required but only %d available\n", pos, datlen, status);
continue;
//exit(1);
}
// buffer does not have the same time as previous buffers
if(!treat_buffer(lastbuffer)) {
return 1;
}
}
bufferset = 0;
// end of file reached --> make last events
make_events();
return found;
}
int XiaEventReader::treat_buffer(uint16_t *buffer) {
//DGF_MAXBUF
uint16_t nwords = buffer[DGF_H_NWORDS];
uint16_t module = buffer[DGF_H_MODULE];
uint16_t fmt = buffer[DGF_H_FORMAT] & DGF_FORMAT_MASK;
uint16_t timehi = buffer[DGF_H_TIMESTAMP];
uint16_t timemi = buffer[DGF_H_TIMESTAMP+1];
uint16_t timelo = buffer[DGF_H_TIMESTAMP+2];
// calc time
uint64_t buftime = get48bittime(timehi, timemi, timelo);
if(fmt==DGF_FORMAT_UNIX) {
uint16_t *evbuf = &buffer[DGF_BUFHEADLEN];
if(evbuf[0]!=1) {
fprintf(stderr, "incorrect start of eventbuffer for unixtimestamp, 1 expected but %u found\n", evbuf[0]);
exit(1);
}
uint16_t *utimebuf = &evbuf[DGF_EVENTHEADLEN];
unix_timestamp = (utimebuf[1] << 16) | utimebuf[0];
unix_timestamp_found = 1;
return 1;
}
// check if current buffer corresponds to previous buffers
int64_t diff = buftime - currenttime;
if((currenttime!=-1)&&(abs_int64(diff)>20)) {
// make events
make_events();
return 0;
}
currenttime = (int64_t)buftime;
// get index of empty buffer
int bufindex = get_emptybuffer();
// set module nr
dgf_buffers[bufindex].module = module;
uint16_t *evbuf = &buffer[DGF_BUFHEADLEN];
while((evbuf-buffer)<nwords) {
uint16_t pat = evbuf[DGF_E_PATTERN];
uint16_t evtimemi = evbuf[DGF_E_TIMESTAMP];
uint16_t evtimelo = evbuf[DGF_E_TIMESTAMP+1];
evbuf = &evbuf[DGF_EVENTHEADLEN];
uint64_t evtime = get48bittime(timehi, evtimemi, evtimelo);
// check if evtime is smaller than buffertime
if(evtime<buftime) {
timehi++;
evtime = get48bittime(timehi, evtimemi, evtimelo);
}
reset_current_dgfevent(bufindex);
set_dgfeventtime(bufindex, evtime);
for(int i=0; i<DGF_NCHAN; i++) {
// check if channel i had a hit
if(!(pat & (1 << i))) {
//printf("pat %d\n", i);
continue;
}
uint16_t chtimelo=0;
uint16_t en=0;
uint16_t ndat=0;
switch(fmt) {
// long format
case 0x0100: case 0x1100: case 0x0101: case 0x1101:
ndat = evbuf[DGF_C_LNWORDS];
chtimelo = evbuf[DGF_C_LTRIGGER_TIME];
en = evbuf[DGF_C_LENERGY];
evbuf = evbuf + ndat;
break;
// PSA format
case 0x0102: case 0x1102:
chtimelo = evbuf[DGF_C_PTRIGGER_TIME];
en = evbuf[DGF_C_PENERGY];
evbuf = evbuf + DGF_PCHANHEADLEN;
break;
// short format
case 0x0103: case 0x1103:
chtimelo = evbuf[DGF_C_STRIGGER_TIME];
en = evbuf[DGF_C_SENERGY];
evbuf = evbuf + DGF_SCHANHEADLEN;
break;
default:
fprintf(stderr, "unknown format: 0x%04x\n", fmt);
exit(1);
}
// calculate
uint64_t chtime = get48bittime(timehi, evtimemi, chtimelo);
if(chtime<evtime) {
uint64_t buf1 = get32bittime(timehi, evtimemi);
uint64_t buf2 = chtimelo;
buf1++;
chtime = buf1 << 16;
chtime = chtime | buf2;
}
// advance to next event
evbuf = &evbuf[ndat];
// add to event array
int veto = (pat & (1 << (i+12))) ? 1 : 0;
add_hit2dgfevent(bufindex, i, en, chtime, veto);
}
dgf_buffers[bufindex].evcnt++;
}
return 1;
}
