newmonitor/get4/Processor.cxx (r3879/r3878)

Go to the documentation of this file.

#include "get4/Processor.h"

#include <stdlib.h>
#include <stdio.h>
#include <math.h>

#include <algorithm>

#include "base/commons.h"
#include "base/ProcMgr.h"

#include "get4/SubEvent.h"

get4::Get4Rec::Get4Rec() :
   used(false),
   fChannels(0)
{
   for(unsigned n=0;n<NumChannels;n++) {
      fRisTm[n] = 0;
      fFalTm[n] = 0;
      lastRisingEdge[0] = 0;
   }
}



get4::Processor::Processor(unsigned rocid, unsigned get4mask) :
   base::SysCoreProc("ROC", rocid),
   fIter(),
   fIter2()
{
   mgr()->RegisterProc(this, roc::proc_RocEvent, rocid);

//   printf("Start histo creation\n");

   fMsgsKind = MakeH1("MsgKind", "kind of messages", 8, 0, 8, "xbin:NOP,-,EPOCH,SYNC,AUX,EPOCH2,GET4,SYS;kind");
   fSysTypes = MakeH1("SysTypes", "Distribution of system messages", 16, 0, 16, "systype");

   CreateBasicHistograms();

//   printf("Histo creation done\n");

   for (unsigned get4=0; (get4<16) && (get4mask!=0); get4++) {
      GET4.push_back(Get4Rec());
      GET4[get4].used = (get4mask & 1) == 1;
      get4mask = get4mask >> 1;
      if (!GET4[get4].used) continue;

      SetSubPrefix("GET4_", get4);

      GET4[get4].fChannels = MakeH1("Channels", "GET4 channels", 8, 0, 4., "ch");

      for(unsigned n=0;n<NumChannels;n++) {
         SetSubPrefix("GET4_", get4, "Ch", n);

         GET4[get4].fRisTm[n] = MakeH1("RisingTm", "rising edge time", 1024*8, 0, 1024.*512, "bin");

         GET4[get4].fRisFineTm[n] = MakeH1("RisingFine", "rising fine time", 128, 0, 128, "bin");

         GET4[get4].fFalTm[n] = MakeH1("FallingTm", "falling edge time", 1024*8, 0, 1024.*512, "bin");

         GET4[get4].fFalFineTm[n] = MakeH1("FallingFine", "falling fine time", 128, 0, 128, "bin");

         GET4[get4].fWidth[n] = MakeH1("Width", "time-over-threshold", 1024, 0, 1024, "bin");

         GET4[get4].fRisRef[n] = 0;
         GET4[get4].fFalRef[n] = 0;
      }
   }

   SetSubPrefix("");

//   printf("Histo creation done\n");

   SetNoSyncSource(); // use not SYNC at all
   SetNoTriggerSignal();

   fRefGet4 = 9999;
   fRefChannel = NumChannels;

   fNumHits = 0;
   fNumBadHits = 0;
}

get4::Processor::~Processor()
{
   //printf("ROC%u   NumHits %7d  NumBad %7d rel = %5.3f\n",
   //      fRocId, fNumHits, fNumBadHits, fNumHits > 0 ? 1.* fNumBadHits / fNumHits : 0.);

   printf("get4::Processor::~Processor\n");
}

void get4::Processor::setRefChannel(unsigned ref_get4, unsigned ref_ch)
{
   fRefGet4 = ref_get4;
   fRefChannel = ref_ch;
   if (!isRefChannel()) return;

   for (unsigned get4=0; get4<GET4.size(); get4++) {
      if (!GET4[get4].used) continue;

      for(unsigned ch=0;ch<NumChannels;ch++) {
         SetSubPrefix("GET4_", get4, "Ch", ch);

         GET4[get4].fRisRef[ch] = MakeH1("Rising_to_Ref", "difference to ref signal, rising edge", 1024*8, -4096., 4096., "bin");
         GET4[get4].fFalRef[ch] = MakeH1("Falling_to_Ref", "difference to ref signal, falling edge", 1024*8, -4096., 4096., "bin");
      }
   }

}


void get4::Processor::AssignBufferTo(roc::Iterator& iter, const base::Buffer& buf)
{
   if (buf.null()) return;

   iter.setFormat(buf().format);
   iter.setRocNumber(buf().boardid);

   unsigned msglen = roc::Message::RawSize(buf().format);

   // we exclude last message which is duplication of SYNC
   iter.assign(buf().buf, buf().datalen - msglen);
}



bool get4::Processor::FirstBufferScan(const base::Buffer& buf)
{
   if (buf.null()) return false;

   AssignBufferTo(fIter, buf);

   roc::Message& msg = fIter.msg();

   unsigned cnt(0), msgcnt(0);

   bool first = true;

//   printf("======== MESSAGES for board %u ==========\n", fRocId);


//   static int doprint = 1;
//   static double mymarker = 4311410447414.;

   while (fIter.next()) {

      cnt++;

      // ignore epoch message at the end of the buffer
      if (fIter.islast() && msg.isEpochMsg()) continue;

      unsigned rocid = msg.getRocNumber();

      if (rocid != GetBoardId()) {
         printf("Message from wrong ROCID %u, expected %u\n", rocid, GetBoardId());
         continue;
      }

      msgcnt++;

      FillH1(fMsgsKind, msg.getMessageType());


//      fIter.printMessage(roc::msg_print_Human);

      // keep time of first non-epoch message as start point of the buffer
      if (first && !msg.isEpochMsg()) {
         first = false;
         buf().local_tm = fIter.getMsgFullTimeD();
      }

      // this time used for histograming and print selection
      double msgtm = (fIter.getMsgFullTime() % 1000000000000LLU)*1e-9;

      FillH1(fALLt, msgtm);


      switch (msg.getMessageType()) {

         case roc::MSG_NOP:
            break;

         case roc::MSG_HIT: {
            printf("FAILURE - no any nXYTER hits in GET4 processor!!!\n");
            exit(5);
            break;
         }

         case roc::MSG_EPOCH: {
            break;
         }

         case roc::MSG_GET4: {
            // must be ignored

            unsigned get4 = msg.getGet4Number();
            if (!get4_in_use(get4)) break;

            unsigned ch = msg.getGet4ChNum();
            unsigned edge = msg.getGet4Edge();
            unsigned ts = msg.getGet4Ts();
            uint64_t fulltm = fIter.getMsgFullTime();

            // fill rising and falling edges together
            FillH1(GET4[get4].fChannels, ch + edge*0.5);

            if ((get4==fRefGet4) && (ch == fRefChannel) && (edge==0)) {
               base::LocalTriggerMarker marker;
               marker.localid = 1000 + get4*16 + ch*2 + edge;
               marker.localtm = fIter.getMsgFullTimeD();

               // printf("Select TRIGGER: %10.9f\n", marker.localtm*1e-9);

               AddTriggerMarker(marker);

            }


            if ((get4==fRefGet4) && (ch == fRefChannel)) {
               // remember time of reference signal
               if (edge==0)
                  fLastRefRising = fulltm;
               else
                  fLastRefFalling = fulltm;
               // recalculate all differences

               for (unsigned ii=0; ii<GET4.size(); ii++) {
                  if (!GET4[ii].used) continue;

                  for(unsigned jj=0;jj<NumChannels;jj++) {
                     if (edge==0)
                        FillH1(GET4[ii].fRisRef[jj], Get4TimeDiff(fLastRefRising, GET4[ii].lastRisingEdge[jj]));
                     else
                        FillH1(GET4[ii].fFalRef[jj], Get4TimeDiff(fLastRefFalling, GET4[ii].lastFallingEdge[jj]));
                  }
               }

            }

            if (edge==0) {
               FillH1(GET4[get4].fRisTm[ch], ts);
               FillH1(GET4[get4].fRisFineTm[ch], ts % 128);
               GET4[get4].lastRisingEdge[ch] = fulltm;

               if (isRefChannel())
                  FillH1(GET4[get4].fRisRef[ch], Get4TimeDiff(fLastRefRising, fulltm));

            } else {
               FillH1(GET4[get4].fFalTm[ch], ts);
               FillH1(GET4[get4].fFalFineTm[ch], ts % 128);
               GET4[get4].lastFallingEdge[ch] = fulltm;

               FillH1(GET4[get4].fWidth[ch], GET4[get4].lastFallingEdge[ch] - GET4[get4].lastRisingEdge[ch]);

               if (isRefChannel())
                  FillH1(GET4[get4].fFalRef[ch], Get4TimeDiff(fLastRefFalling, fulltm));
            }


            break;
         }

         case roc::MSG_EPOCH2: {
            break;
         }

         case roc::MSG_SYNC: {
            unsigned sync_ch = msg.getSyncChNum();
            // unsigned sync_id = msg.getSyncData();

            FillH1(fSYNCt[sync_ch], msgtm);

            // for the moment use DABC format where SYNC is always second message
            if (sync_ch == fSyncSource) {

               base::SyncMarker marker;
               marker.uniqueid = msg.getSyncData();
               marker.localid = msg.getSyncChNum();
               marker.local_stamp = fIter.getMsgFullTime();
               marker.localtm = fIter.getMsgFullTimeD(); // nx gave us time in ns

               // printf("ROC%u Find sync %u tm %6.3f\n", rocid, marker.uniqueid, marker.localtm*1e-9);
               // if (marker.uniqueid > 11798000) exit(11);

               // marker.globaltm = 0.; // will be filled by the StreamProc
               // marker.bufid = 0; // will be filled by the StreamProc

               AddSyncMarker(marker);
            }

            if (fTriggerSignal == (10 + sync_ch)) {

               base::LocalTriggerMarker marker;
               marker.localid = 10 + sync_ch;
               marker.localtm = fIter.getMsgFullTimeD();

               AddTriggerMarker(marker);
            }

            break;
         }

         case roc::MSG_AUX: {
            unsigned auxid = msg.getAuxChNum();

            FillH1(fAUXt[auxid], msgtm);

            if (fTriggerSignal == auxid) {

               base::LocalTriggerMarker marker;
               marker.localid = auxid;
               marker.localtm = fIter.getMsgFullTimeD();

               AddTriggerMarker(marker);
            }

            break;
         }

         case roc::MSG_SYS: {
            FillH1(fSysTypes, msg.getSysMesType());
            break;
         }
      } // switch

   }

   FillMsgPerBrdHist(msgcnt);

   return true;
}

unsigned get4::Processor::GetTriggerMultipl(unsigned indx)
{
   get4::SubEvent* ev = (get4::SubEvent*) fGlobalTrig[indx].subev;

   return ev ? ev->fExtMessages.size() : 0;
}


bool get4::Processor::SecondBufferScan(const base::Buffer& buf)
{
   if (buf.null()) return false;

//   printf("Start second buffer scan left %u  right %u  size %u\n", lefttrig, righttrig, fGlobalTrig.size());
//   for (unsigned n=0;n<fGlobalTrig.size();n++)
//      printf("TRIG %u %12.9f flush:%u\n", n, fGlobalTrig[n].globaltm*1e-9, fGlobalTrig[n].isflush);

   AssignBufferTo(fIter2, buf);

   roc::Message& msg = fIter2.msg();

   unsigned cnt(0);

   unsigned help_index = 0; // special index to simplify search of respective syncs

   while (fIter2.next()) {

      cnt++;

      // ignore epoch message at the end of the buffer
      if (fIter2.islast() && msg.isEpochMsg()) continue;

      unsigned rocid = msg.getRocNumber();

      if (rocid != GetBoardId()) continue;

//      printf("Scan message %u type %u\n", cnt, msg.getMessageType());

      switch (msg.getMessageType()) {

         case roc::MSG_NOP:
            break;

         case roc::MSG_SYNC:
         case roc::MSG_AUX:
         case roc::MSG_HIT:
            break;


         case roc::MSG_GET4: {
            base::GlobalTime_t localtm = fIter2.getMsgFullTimeD();

            base::GlobalTime_t globaltm = LocalToGlobalTime(localtm, &help_index);

            unsigned indx = TestHitTime(globaltm, get4_in_use(msg.getGet4Number()));

            if ((indx < fGlobalTrig.size()) && !fGlobalTrig[indx].isflush) {
               get4::SubEvent* ev = (get4::SubEvent*) fGlobalTrig[indx].subev;

               if (ev==0) {
                  ev = new get4::SubEvent;
                  fGlobalTrig[indx].subev = ev;
               }

               ev->fExtMessages.push_back(get4::MessageExtended(msg, globaltm));
            }
            break;
         }

         case roc::MSG_EPOCH:
            break;

         case roc::MSG_EPOCH2:
            break;

         case roc::MSG_SYS:
            break;
      } // switch

   }

   return true;
}

void get4::Processor::SortDataInSubEvent(base::SubEvent* subev)
{
   get4::SubEvent* nxsub = (get4::SubEvent*) subev;

   std::sort(nxsub->fExtMessages.begin(), nxsub->fExtMessages.end());
}