myReadOutDesy.cpp

#include <cstdio>
#include <iostream>
#include <fstream>
#include <cstdint>
#include <csignal>
#include <unistd.h>
#include <sys/time.h>

#include "myModules.h"                         // VME Module Id.s
#include "myV2718.h"                           // CAEN V2718 VME Bridge class
#include "myV513.h"                            // CAEN V513 I/O REG class
#include "myV792AC.h"                          // CAEN V792AC QDC class


using namespace std;

#define TRIGGER_IN_CHAN 7                    // channel 7
#define TRIGGER_MASK 7                       // channels 0, 1, 2
#define TRIGGER_DISA_ALL 15                  // channel 15
#define TRIGGER_UNLK_CHAN 14                 // channel 14
#define SCALER_RESET 13                      // channel 13
#define PEDESTAL_VETO 12                     // channel 12

#define NSLEEPT 1000                         // for sleeping 1000 ns

void dlwait( char* msg )
 {
  cout << "Press any key to continue ... " << msg;
  getchar();
  cout << endl;
 }

void nsleep ( uint32_t ns )
 {
  struct timespec ndelay = { 0, ns };
  nanosleep ( &ndelay, NULL );
 }

inline void enableTriggers( v513& ioreg ) { ioreg.clearOutputBit( TRIGGER_DISA_ALL ); }
inline void disableTriggers( v513& ioreg ) { ioreg.setOutputBit( TRIGGER_DISA_ALL ); }

inline void enablePedestal( v513& ioreg ) { ioreg.clearOutputBit( PEDESTAL_VETO ); }
inline void disablePedestal( v513& ioreg ) { ioreg.setOutputBit( PEDESTAL_VETO ); }

inline uint32_t ioregValue( v513& ioreg ) { return ioreg.readInputRegister(); }

inline bool isTrig( v513& ioreg ) { return (ioreg.readInputRegister() & 1<<TRIGGER_IN_CHAN); }
inline uint32_t getTrigMask( v513& ioreg ) { return (ioreg.readInputRegister() & TRIGGER_MASK); }

inline void resetScaler( v513& ioreg )
 {
  nsleep(NSLEEPT);
  ioreg.setOutputBit( SCALER_RESET );
  nsleep(NSLEEPT);
  ioreg.clearOutputBit( SCALER_RESET );
 }

inline void unlockTrigger( v513& ioreg )
 {
  nsleep(NSLEEPT);
  ioreg.setOutputBit( TRIGGER_UNLK_CHAN );
  nsleep(NSLEEPT);
  ioreg.clearOutputBit( TRIGGER_UNLK_CHAN );
 }

int initV513( v513& ioreg )
 {
  uint16_t w;
  ioreg.read16phys(0xFE, &w);
  cout << hex << " V513 FE " << w << endl;
  nsleep(NSLEEPT);
  ioreg.read16phys(0xFC, &w);
  cout << hex << " V513 FC " << w << endl;
  nsleep(NSLEEPT);
  ioreg.read16phys(0xFA, &w);
  cout << hex << " V513 FA " << w << endl;
  nsleep(NSLEEPT);
  w=0;
  ioreg.write16phys(0x48, w);
  nsleep(NSLEEPT);
  ioreg.write16phys(0x46, w);
  nsleep(NSLEEPT);
  ioreg.write16phys(0x42, w);
  nsleep(NSLEEPT);

  ioreg.reset();
  nsleep(NSLEEPT);

  ioreg.read16phys(0x04, &w);
  cout << hex << " V513 0x4 " << w << endl;
  nsleep(NSLEEPT);

  for (int i=0; i<8; i++)
   {
    int reg = 0x10+i*2;
    ioreg.setChannelInput(i);
    nsleep(NSLEEPT);
    ioreg.read16phys(reg, &w);
    cout << hex << " V513 input ch " << i << " reg " << reg << " " << w << endl;
    nsleep(NSLEEPT);
   }

  for (int i=8; i<16; i++)
   {
    int reg = 0x10+i*2;
    ioreg.setChannelOutput(i);
    nsleep(NSLEEPT);
    ioreg.read16phys(reg, &w);
    cout << hex << " V513 output ch " << i << " reg " << reg << " " << w << endl;
    nsleep(NSLEEPT);
   }

  disableTriggers( ioreg );

  return 0;
 }

void initV792( v792ac& v792adc )
 {
  v792adc.print();
  v792adc.reset();
  v792adc.disableSlide();
  v792adc.disableOverflowSupp();
  v792adc.disableZeroSupp();
  v792adc.clearEventCounter();
  v792adc.clearData();
  uint16_t ped;
  ped = v792adc.getI1();
  uint32_t bid = v792adc.id();
  cout << "v792ac addr 0x" << hex << bid << dec
       << ": default ped I1 value is " << ped << endl;

  ped = 255;                  // >~ minimum possible Iped (see manual par. 2.1.2.1)
  v792adc.setI1(ped);       // set I1 current to define pedestal position
  ped =v792adc.getI1();     // in the range [0x00-0xFF] (see manual par. 4.34.2)
  uint16_t thr = 0x00;
  v792adc.setChannelThreshold(thr);
  cout << "v792ac addr 0x" << hex << bid << dec << ": now ped  I1 value is " << ped << " thr " << thr << endl;
  
  v792adc.enableChannels();

  uint16_t reg1 = v792adc.getStatus1();
  uint16_t reg2 = v792adc.getControl1();
  uint16_t reg3 = v792adc.getStatus2();
  uint16_t reg4 = v792adc.getBit1();
  uint16_t reg5 = v792adc.getBit2();
  uint32_t evc = v792adc.eventCounter();
  cout << "v792ac addr 0x" << hex << bid << " status 1 0x" << reg1
       << " control 1 0x" << reg2 << " status 2 0x" << reg3 << " bit 1 0x" << reg4
       << " bit 2 0x" << reg5 << dec << endl;
  cout << "v792ac event counter " << evc << endl;
 }

volatile bool abort_run(false);
volatile bool pause_run(false);

void cntrl_c_handler ( int32_t sig )
 {
  time_t timestr = time(NULL);
  char * stime = ctime(&timestr);
  stime[24] = 0;
  fprintf(stderr,"%s cntrl_c_handler: sig%d\n\n", stime, sig);
  fprintf(stderr,"aborting run\n");
  abort_run = true;
 }

void sigusr1_handler ( int32_t sig )
 {
  time_t timestr = time(NULL);
  char * stime = ctime(&timestr);
  stime[24] = 0;
  fprintf(stderr,"%s sigusr1_handler: sig%d  pause_run is %d\n\n", stime, sig, pause_run);
  pause_run = not pause_run;
  pause_run ?  fprintf(stderr,"pausing run\n") : fprintf(stderr,"resuming run\n");
 }

uint32_t readEvent ( v792ac& v792adc, uint32_t phEvents, uint32_t* buffer )
 {
  nsleep(NSLEEPT);
  // v792adc.setEvents( phEvents, 0 );   // old version < 15.08.2021
  v792adc.setEvents( phEvents );
  uint32_t size = v792adc.readEvent( buffer );
  return size;
 }

int evt_realign ( v513& ion, v792ac& qdc, uint32_t* Tcts, ofstream& ofs, uint32_t& phytr )
 {
  uint32_t myBuffer[256];
  uint32_t tmask = getTrigMask( ion );
  cout << " RRR " << dec << phytr << " " << true << " tmask " << hex << tmask << " ioreg " << ioregValue( ion ) << endl;
  if (tmask & 1) Tcts[0] ++;    // physics triggers
  if (tmask & 2) Tcts[1] ++;    // pedestal triggers
  if (tmask & 4) Tcts[2] ++;
  (Tcts[1]*10 < Tcts[0]) ? enablePedestal( ion ) : disablePedestal( ion );
  uint32_t size = readEvent ( qdc, phytr, myBuffer );
  uint32_t head = myBuffer[0];
  uint32_t trail = myBuffer[size-1] & 0xff000000;
  uint32_t evnum = myBuffer[size-1] & 0xffffff;
  if (head != 0xfa002000) cout << dec << phytr << " Error in header " << hex << head << endl;
  if (trail != 0xfc000000) cout << dec << phytr << " Error in trailer " << hex << trail << endl;
  if (size != 34) cout << dec << phytr << " Size error " << size << endl;
  if (evnum != phytr)
   {
    cout << dec << phytr << " Error in event number " << evnum << endl;
    tmask |= 0xffff0000;
   }
  ofs << dec << " ev # " << phytr << " tow cts " << Tcts[0] << " " << Tcts[1] << " " << Tcts[2] << " trigger mask " << tmask << " values:";
  for (uint32_t j=1; j<size-1; j++)
   {
    uint32_t chan = (myBuffer[j] >> 16) & 0x1f;
    ofs << " " << dec << chan << " " << hex << (myBuffer[j] & 0x3fff);
   }
  ofs << endl;
  phytr ++;
  return 0;
 }

int main( int argc, char** argv )
 {
  ofstream ofs;
  ifstream ifs;

  uint32_t phytr(0);
  uint32_t Tcts[3] = { 0, 0, 0 };
  uint32_t myBuffer[256];

  signal(SIGINT, cntrl_c_handler);      // Control-C handler
  signal(SIGUSR1, sigusr1_handler);      // Control-C handler

  v2718 v2718_1("/V2718/cvA24_U_DATA/0");            // VME interface
  v513 ion(0xa00000,"/V2718/cvA24_U_DATA/0");          // I/O register
  v792ac qdc(0x05000000,"/V2718/cvA32_U_DATA/0");      // V792 QDC

  // Init I/O Register and QDC
  initV513 ( ion );

  initV792 ( qdc );
  sched_yield();

  disableTriggers( ion );

  resetScaler( ion );
  unlockTrigger( ion );
  // dlwait(" ");
  enableTriggers( ion );

  time_t stime = time(NULL);
  cout << " time is " << stime << endl;
  struct tm mytime;
  localtime_r( &stime, &mytime );

  cout << mytime.tm_year+1900 << "." << mytime.tm_mon+1 << "." << mytime.tm_mday << "." << mytime.tm_hour << ":" << mytime.tm_min << ":" << mytime.tm_sec << endl;
  cout << asctime(&mytime) << endl;

  uint32_t runnbr;
  ifs.open ("/home/dreamtest/working/run.number",std::ofstream::in);
  ifs >> runnbr;
  ifs.close();
  ofs.open ("/home/dreamtest/working/run.number",std::ofstream::out);
  ofs << ++runnbr;
  ofs.close();
  

  stime -= 1600000000;
  string fname = "desydata.run" + to_string(runnbr) + ".txt";

  if (argc == 2) fname = argv[1];

  cout << " filename is " << fname << endl;

  ofs.open (fname.c_str(), std::ofstream::out);

  time_t tr0 = time(NULL);
  bool running(true);

  while (1)
   {
    if (running and pause_run)
     {
      disableTriggers(ion);
      running = false;
  time_t tr = time(NULL) - tr0;
  cout << dec << tr << "(sec) " << phytr << "(evt) ioreg " << ioregValue( ion ) << " running was true (now false) and pause_run true" << endl;
     }
    else if (not running and not pause_run)
     {
      enableTriggers(ion);
      running = true;
  time_t tr = time(NULL) - tr0;
  cout << dec << tr << "(sec) " << phytr << "(evt) ioreg " << ioregValue( ion ) << " running was false (now true) and pause_run false" << endl;
     }
    bool isPh = isTrig( ion );
    if (isPh)
     {
      uint32_t tmask = getTrigMask( ion );
if ((phytr % 10000) == 0)
 {
  time_t tr = time(NULL) - tr0;
  double rate = phytr/double(tr);
  cout << dec << tr << "(sec) " << phytr << "(evt) " << " rate " << rate << "(evt/sec) isPh " << isPh << " tmask " << hex << tmask << " ioreg " << ioregValue( ion ) << endl;
 }
      if (tmask & 1) Tcts[0] ++;    // physics triggers
      if (tmask & 2) Tcts[1] ++;    // pedestal triggers
      if (tmask & 4) Tcts[2] ++;
      (Tcts[1]*10 < Tcts[0]) ? enablePedestal( ion ) : disablePedestal( ion );
      uint32_t size = readEvent ( qdc, phytr, myBuffer );
      uint32_t head = myBuffer[0];
      uint32_t trail = myBuffer[size-1] & 0xff000000;
      uint32_t evnum = myBuffer[size-1] & 0xffffff;
      if (head != 0xfa002000) cout << dec << phytr << " Error in header " << hex << head << endl;
      if (trail != 0xfc000000) cout << dec << phytr << " Error in trailer " << hex << trail << endl;
      if (size != 34) cout << dec << phytr << " Size error " << size << endl;
      if (evnum != phytr)
       {
        cout << dec << phytr << " Error in event number " << evnum << endl;
        tmask |= 0xffff0000;
       }
      ofs << dec << " ev # " << phytr << " tow cts " << Tcts[0] << " " << Tcts[1] << " " << Tcts[2] << " trigger mask " << hex << tmask << " values:";
      for (uint32_t j=1; j<size-1; j++)
	 {
	  uint32_t chan = (myBuffer[j] >> 16) & 0x1f;
	  ofs << " " << dec << chan << " " << hex << (myBuffer[j] & 0x3fff);
	 }
      ofs << endl;
      phytr ++;
      // while (phytr <= evnum) phytr ++; // evt_realign ( ion, qdc, Tcts, ofs, phytr );
      if (abort_run) disableTriggers( ion );
      unlockTrigger( ion );
     }
    if (abort_run) break;
    nsleep( NSLEEPT );
   }
  ofs.close();

  time_t etime = time(NULL);
  etime -= 1600000000;
  uint32_t tdiff = etime-stime;
  double rate = phytr/double(tdiff);
  cout << dec << " time is " << etime << " diff(sec) " << (etime-stime) << " events " << phytr
	<< " rate(Hz) " << rate << endl;

  return 0;
 }