onlinemonitor/rocmonitor/TRocRec.cxx (r4864/r3423)

Go to the documentation of this file.

#include "TRocRec.h"

#include "TLatex.h"

TRocRec::TRocRec()
{
   fRocId = 0;

   fCurrEpoch = 0;

   fAUXch = 0;
   fMsgTypes = 0;
   fSysTypes = 0;
   fSysUserTypes = 0;
   fTrigger_AllNX = 0;
   fTrigger_AllNX_100 = 0;
   fTriggerWind = 0; // time window to assign hits to event
   fAUXWind = 0; // time window to assign AUX to event
   fEvntMultipl = 0; // event multiplicity in time window after trigger signal

   fLastTm = 0;  //
   fLastTriggerTm = 0;

   for (int n=0;n<MAX_AUX;n++) fLastAuxTm[n] = 0;  
   for (int n=0;n<MAX_SYNC;n++) {
      fLastSyncTm[n] = 0;
      fLastSyncId[n] = 0;
   }
   fStartSyncTm = 0;
   fStopSyncTm = 0;

   for (int nx=0;nx<MAX_NX;nx++) {
      fADCs[nx] = 0; // individual ADC histogram for each ROC/nxyter
      fTrigADCs[nx] = 0; // individual ADC histogram for each ROC/nxyter
      fBaseline[nx] = 0; // individual baseline histogram for each ROC/nxyter
      fADCs_wo_baseline[nx] = 0; // individual ADC wo baseline for each ROC/nxyter
      fChs[nx] = 0; // individual ADC histogram for each ROC/nxyter
      fTrigger_NX[nx] = 0; // time shift between last trigger and nXYTER hit

      for (int ch=0; ch<NUM_NXCHANNELS; ch++) fNxTm[nx][ch] = 0;
      for (int ndiff = 0; ndiff < NUM_DIFFS; ndiff++)
         fTmDiff[nx][ndiff] = 0;
   }

   for (int aux=0;aux<MAX_AUX;aux++)
      fAUXt[aux] = 0; // time distribution of aux signals

   fIsEventComplete = kFALSE;
   fHasNewData = kFALSE;
   fHasNewTrigger=kFALSE;
   fHasEndOfBuffer=kFALSE;
   fTriggersPerBuffer=0;
   fLoopsPerBuffer=0;
   fDabcSeparator=-1;

   bCollectBaseline = kFALSE;
   bIgnoreData = kFALSE;
   fIgnoreCnt = 0;
   fIgnoreTime = 0;       // time when ignore message was generated
   fCollectTime = 0;      // time when collect message was generated

   /*************************************/

   for (int n=0;n<MAX_GET4;n++) {
      fLastEpoch2[n] = 0; // last epoch2 time
      fLastEpoch2Sync[n] = 0; // last epoch2 time with sync=1
      fGet4Channels[n] = 0;
      fGet4Tm[n] = 0; // distribution of GET4 time stamps
      for(int k=0;k<MAX_GET4_CH;k++) {
         fGet4EdgeCnt[n][k] = 0;
         fGet4ErrCnt[n][k] = 0;
         fGet4ChangeCnt[n][k] = 0;
         fGet4SuspiciousCnt[n][k] = 0;
         for (int fl=0;fl<2;fl++) {
            fGet4TmCh[n][k][fl] = 0;
            fGet4ChCnt[n][k][fl] = 0;
            fLastGet4Tm[n][k][fl] = 0;
         }
      }

      fGet4DiffSync[n] = 0;
      fGet4Info[n] = 0;
      fGet4Trigger[n] = 0;
      fGet4Trigger100[n] = 0;
      fLastEpoch2Info[n] = 0;
      fEpoch2cnt[n] = 0;
   }
   fGet4SysTypes = 0; // Get4 System message type
   fGet4Chips = 0;

}

void TRocRec::ProcessEpoch2(unsigned g4id, uint32_t epoch2, uint32_t sync, Bool_t debug)
{
   if (g4id>=MAX_GET4) return;

   uint32_t diff = epoch2 - fLastEpoch2[g4id];

   if ((fLastEpoch2[g4id]!=0) && (diff!=1) && (diff<10) && debug)
      printf("EPOCH2 JUMP ROC:%d epoch2:%u Get4:%u diff:%u\n", fRocId, epoch2, g4id, diff);

   if (sync && (epoch2 % 25 != 0) && debug)
      printf("ROC:%d Get4:%u EPOCH2 %u with sync mod 25 = %u \n", fRocId, g4id, epoch2, epoch2 % 25);

   // check about every 600 ms
   if (sync && (epoch2 % 25000 == 0)) {
      for (int nch=0;nch<MAX_GET4_CH;nch++) {
         if ((fGet4ChangeCnt[g4id][nch] > 10) &&
             (fGet4SuspiciousCnt[g4id][nch] > 0.2 * fGet4ChangeCnt[g4id][nch]) && debug)
                printf("ROC:%d Get4:%u Channel:%u Has suspicious difference between edges fullcnt:%d suspiciouscnt:%d\n",fRocId, g4id, nch, fGet4ChangeCnt[g4id][nch], fGet4SuspiciousCnt[g4id][nch]);
         if ((fGet4ErrCnt[g4id][nch]>4) && debug)
             printf("ROC:%d Get4:%u Channel:%u has too many contigious edges errors cnt:%d\n",fRocId, g4id, nch, fGet4ErrCnt[g4id][nch]);

         fGet4ChangeCnt[g4id][nch]=0;
         fGet4SuspiciousCnt[g4id][nch]=0;
         fGet4ErrCnt[g4id][nch]=0;
      }
   }

   fLastEpoch2[g4id] = epoch2;
   fEpoch2cnt[g4id]++;

   if (sync) fLastEpoch2Sync[g4id] = epoch2;

   if (fGet4Info[g4id]==0) return;

   if (fLastEpoch2Info[g4id]!=0) {

//      printf("Diff ROC:%d G4:%u %u %g \n", fRocId, g4id, fLastEpoch2[g4id] - fLastEpoch2Info[g4id], roc::Message::FullTimeStamp2(fLastEpoch2[g4id] - fLastEpoch2Info[g4id], 0) / 2e10);

      // 40000 GET4 epochs is about 1 sec
      if ((fLastEpoch2[g4id] > fLastEpoch2Info[g4id]) && (fLastEpoch2[g4id] < (fLastEpoch2Info[g4id]+40000))) return;

      if ((fLastEpoch2[g4id] >= (fLastEpoch2Info[g4id]+40000)) && (fEpoch2cnt[g4id]>0)) {

         // we take into account only effective time - excluding parts, skipped by streamserver
         double tm = roc::Message::FullTimeStamp2(fEpoch2cnt[g4id], 0) / 2e10;

         TString txt = "#splitline";

         txt+=Form("{#scale[1.0]{#color[2]{ROC%d GET4:%u}}}", fRocId, g4id);

         for(int k=0;k<MAX_GET4_CH;k++) {
            if ((fGet4ChCnt[g4id][k][0]==0) && (fGet4ChCnt[g4id][k][1]==0)) continue;

            double f0 = fGet4ChCnt[g4id][k][0]/tm;
            double f1 = fGet4ChCnt[g4id][k][1]/tm;

            // printf("Roc:%d Get4:%u Ch:%d %g %g\n", fRocId, g4id, k, f0,f1);

            txt+=Form("{#scale[0.8]{#color[4]{Ch%d fall:%6.1f Hz rise:%6.1f Hz diff:%4.1f}}}", k, f0, f1, f0-f1);
         }

         fGet4Info[g4id]->SetText(0.1, 0.7, txt.Data());
      }
   }

   fLastEpoch2Info[g4id] = fLastEpoch2[g4id];

   fEpoch2cnt[g4id] = 0;
   for(int k=0;k<MAX_GET4_CH;k++)
      for (int fl=0;fl<2;fl++)
         fGet4ChCnt[g4id][k][fl] = 0;
}