CbmStar2019MonitorTof.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                    CbmStar2019MonitorTof                          -----
// -----               Created 10.07.2018 by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#ifndef CbmStar2019MonitorTof_H
#define CbmStar2019MonitorTof_H
 
#include "Timeslice.hpp"
#include "gDpbMessv100.h"
//#include "CbmTofStarData.h"
//#include "CbmTofStarData2018.h"
 
#include "CbmMcbmUnpack.h"
 
#include "Rtypes.h"
#include "TClonesArray.h"
 
#include <chrono>
#include <map>
#include <vector>
 
class CbmStar2019TofPar;
 
class TCanvas;
class TH1;
class TH2;
class TProfile;
class TProfile2D;
 
class CbmStar2019MonitorTof : public CbmMcbmUnpack {
public:
  CbmStar2019MonitorTof();
  virtual ~CbmStar2019MonitorTof();
 
  virtual Bool_t Init();
 
  virtual Bool_t DoUnpack(const fles::Timeslice& ts, size_t component);
 
  virtual void Reset();
 
  virtual void Finish();
 
  void SetParContainers();
 
  Bool_t InitContainers();
 
  Bool_t ReInitContainers();
 
  void SetMsLimitLevel(size_t uAcceptBoundaryPct = 100) {
    fuMsAcceptsPercent = uAcceptBoundaryPct;
  }
  size_t GetMsLimitLevel() { return fuMsAcceptsPercent; }
 
  virtual void AddMsComponentToList(size_t component, UShort_t usDetectorId);
  virtual void SetNbMsInTs(size_t uCoreMsNb, size_t uOverlapMsNb);
  void SetIgnoreMsOverlap(Bool_t bEnaFlag = kTRUE) {
    fbIgnoreOverlapMs = bEnaFlag;
  }
  void SetMsOverlap(size_t uOverlapMsNb = 1) { fuOverlapMsNb = uOverlapMsNb; }
  size_t GetMsOverlap() { return fuOverlapMsNb; }
 
  inline void SetFitZoomWidthPs(Double_t inZoomWidth = 1000.0) {
    fdFitZoomWidthPs = inZoomWidth;
  }
  inline void SetHistoryHistoSize(UInt_t inHistorySizeSec = 1800) {
    fuHistoryHistoSize = inHistorySizeSec;
  }
  inline void SetHistoryHistoSizeLong(UInt_t inHistorySizeMin = 1800) {
    fuHistoryHistoSizeLong = inHistorySizeMin;
  }
 
  inline void EnablePulserMode(Bool_t bEnaFlag = kTRUE) {
    fbPulserModeEnable = bEnaFlag;
  }
  inline void EnableCoincidenceMaps(Bool_t bEnaFlag = kTRUE) {
    fbCoincMapsEnable = bEnaFlag;
  }
 
  void SaveAllHistos(TString sFileName = "");
  void ResetAllHistos();
  void ResetEvolutionHistograms();
  void ResetLongEvolutionHistograms();
  void UpdateNormedFt();
  void UpdateZoomedFit();
 
private:
  std::vector<size_t> fvMsComponentsList;  
  size_t fuNbCoreMsPerTs;                  
  size_t fuNbOverMsPerTs;                  
  Bool_t
    fbIgnoreOverlapMs;  
 
  size_t
    fuMsAcceptsPercent; 
  size_t fuTotalMsNb; 
  size_t
    fuOverlapMsNb; 
  size_t fuCoreMs; 
  Double_t fdMsSizeInNs;
  Double_t fdTsCoreSizeInNs;
  UInt_t fuMinNbGdpb;
  UInt_t fuCurrNbGdpb;
 
  CbmStar2019TofPar* fUnpackPar;  
  UInt_t fuNrOfGdpbs;             // Total number of GDPBs in the system
  UInt_t fuNrOfFeePerGdpb;        // Number of FEBs per GDPB
  UInt_t fuNrOfGet4PerFee;        // Number of GET4s per FEE
  UInt_t fuNrOfChannelsPerGet4;   // Number of channels in each GET4
 
  UInt_t fuNrOfChannelsPerFee;   // Number of channels in each FEE
  UInt_t fuNrOfGet4;             // Total number of Get4 chips in the system
  UInt_t fuNrOfGet4PerGdpb;      // Number of GET4s per GDPB
  UInt_t fuNrOfChannelsPerGdpb;  // Number of channels per GDPB
 
  UInt_t fuNrOfGbtx;
  UInt_t fuNrOfModules;
  std::vector<Int_t> fviNrOfRpc;
  std::vector<Int_t> fviRpcType;
  std::vector<Int_t> fviRpcSide;
  std::vector<Int_t> fviModuleId;
 
  const UInt_t kuNbFeePerGbtx  = 5;
  const UInt_t kuNbGbtxPerGdpb = 6;
 
  UInt_t fuRawDataPrintMsgNb;
  UInt_t fuRawDataPrintMsgIdx;
  Bool_t fbPrintAllHitsEnable;
  Bool_t fbPrintAllEpochsEnable;
  Bool_t fbPulserModeEnable;
  Bool_t fbCoincMapsEnable;
 
  uint64_t fulCurrentTsIndex;  // Idx of the current TS
  size_t fuCurrentMs;          // Idx of the current MS in TS (0 to fuTotalMsNb)
  Double_t fdMsIndex;          // Time in ns of current MS from its index
  UInt_t fuGdpbId;  // Id (hex number) of the GDPB for current message
  UInt_t
    fuGdpbNr;  // running number (0 to fuNrOfGdpbs) of the GDPB for current message
  UInt_t
    fuGet4Id;  // running number (0 to fuNrOfGet4PerGdpb) of the Get4 chip of a unique GDPB for current message
  UInt_t
    fuGet4Nr;  // running number (0 to fuNrOfGet4) of the Get4 chip in the system for current message
  Int_t fiEquipmentId;
  std::vector<int> fviMsgCounter;
 
  std::vector<ULong64_t> fvulGdpbTsMsb;
  std::vector<ULong64_t> fvulGdpbTsLsb;
  std::vector<ULong64_t> fvulStarTsMsb;
  std::vector<ULong64_t> fvulStarTsMid;
  std::vector<ULong64_t> fvulGdpbTsFullLast;
  std::vector<ULong64_t> fvulStarTsFullLast;
  std::vector<UInt_t> fvuStarTokenLast;
  std::vector<UInt_t> fvuStarDaqCmdLast;
  std::vector<UInt_t> fvuStarTrigCmdLast;
 
  std::vector<ULong64_t> fvulCurrentEpoch;  
  std::vector<Bool_t> fvbFirstEpochSeen;    
  std::vector<ULong64_t>
    fvulCurrentEpochCycle;  
  std::vector<ULong64_t> fvulCurrentEpochFull;  
 
  ULong64_t fulCurrentEpochTime; 
  std::map<UInt_t, UInt_t> fGdpbIdIndexMap;
 
  std::vector<std::vector<gdpbv100::Message>> fvmEpSupprBuffer;
 
  std::vector<UInt_t> fvuFeeNbHitsLastMs;   
  std::vector<Double_t> fdTsLastPulserHit;  
 
  std::vector<std::vector<UInt_t>>
    fvuCoincNbHitsLastMs;  
  std::vector<std::vector<Double_t>>
    fvdCoincTsLastHit;  
 
  // Default value for nb bins in Pulser time difference histos
  const UInt_t kuNbBinsDt = 2000;
  Double_t dMinDt;
  Double_t dMaxDt;
  // Default number of FEE per channels histograms
  UInt_t fuNbFeePlot;
  UInt_t fuNbFeePlotsPerGdpb;
  // Evolution plots control
  Double_t
    fdStartTime; 
  Double_t
    fdStartTimeLong; 
  Double_t
    fdStartTimeMsSz; 
  UInt_t fuHistoryHistoSize; 
  UInt_t
    fuHistoryHistoSizeLong; 
                            // Pulser plots
  Double_t fdLastRmsUpdateTime;
  Double_t fdFitZoomWidthPs;
  // Flesnet
  TCanvas* fcMsSizeAll;
  std::vector<TH1*> fvhMsSzPerLink;
  std::vector<TProfile*> fvhMsSzTimePerLink;
 
  // Messages types and flags
  TH1* fhMessType;
  TH1* fhSysMessType;
  TH2* fhGet4MessType;
  TH2* fhGet4ChanScm;
  TH2* fhGet4ChanErrors;
  TH2* fhGet4EpochFlags;
  TH2* fhGdpbMessType;
  TH2* fhGdpbSysMessType;
  TH2* fhGdpbSysMessPattType;
  TH2* fhGdpbEpochFlags;
  TH2* fhGdpbEpochSyncEvo;
  TH2* fhGdpbEpochMissEvo;
  std::vector<TH2*> fvhGdpbGet4MessType;
  std::vector<TH2*> fvhGdpbGet4ChanScm;
  std::vector<TH2*> fvhGdpbGet4ChanErrors;
  TH2* fhScmScalerCounters;
  TH2* fhScmDeadtimeCounters;
  TH2* fhScmSeuCounters;
  TH2* fhScmSeuCountersEvo;
  TH2* fhPatternMissmatch;
  TH2* fhPatternEnable;
  TH2* fhPatternResync;
  std::vector<TH2*> fvhGdpbPatternMissmatchEvo;
  std::vector<TH2*> fvhGdpbPatternEnableEvo;
  std::vector<TH2*> fvhGdpbPatternResyncEvo;
 
  // Hit messages
  std::vector<TH2*> fvhRawFt_gDPB;
  std::vector<TH2*> fvhRawTot_gDPB;
  std::vector<TH1*> fvhChCount_gDPB;
  std::vector<TH2*> fvhChannelRate_gDPB;
  std::vector<TH2*> fvhRemapTot_gDPB;
  std::vector<TH1*> fvhRemapChCount_gDPB;
  std::vector<TH2*> fvhRemapChRate_gDPB;
  std::vector<TH1*> fvhFeeRate_gDPB;
  std::vector<TH1*> fvhFeeErrorRate_gDPB;
  std::vector<TProfile*> fvhFeeErrorRatio_gDPB;
  std::vector<TH1*> fvhFeeRateLong_gDPB;
  std::vector<TH1*> fvhFeeErrorRateLong_gDPB;
  std::vector<TProfile*> fvhFeeErrorRatioLong_gDPB;
 
  const UInt_t kuNbFeeSide = 5;
  std::vector<TH2*> fvhRemapTotSideA_mod;
  std::vector<TH2*> fvhRemapTotSideB_mod;
  std::vector<TH1*> fvhModRate;
  std::vector<TH1*> fvhModErrorRate;
  std::vector<TProfile*> fvhModErrorRatio;
 
  std::vector<TH1*> fvhTokenMsgType;
  std::vector<TH1*> fvhTriggerRate;
  std::vector<TH2*> fvhCmdDaqVsTrig;
  std::vector<TH2*> fvhStarTokenEvo;
  std::vector<TProfile*> fvhStarTrigGdpbTsEvo;
  std::vector<TProfile*> fvhStarTrigStarTsEvo;
 
  const Double_t kdMaxDtPulserPs = 300e3;
  std::vector<std::vector<TH1*>> fvhTimeDiffPulser;
  TH2* fhTimeMeanPulser;
  TH2* fhTimeRmsPulser;
  TH2* fhTimeRmsZoomFitPuls;
  TH2* fhTimeResFitPuls;
  std::vector<TProfile*> fvhPulserTimeDiffEvoGbtxGbtx;
  std::vector<std::vector<TProfile*>> fvvhPulserTimeDiffEvoGdpbGdpb;
 
  std::vector<TH2*> fvhCoincMapAllChanGdpb;
  std::vector<TProfile2D*> fvhCoincMeanAllChanGdpb;
 
  void CreateHistograms();
 
  void ProcessEpochCycle(uint64_t ulCycleData);
  void FillHitInfo(gdpbv100::Message);
  void FillEpochInfo(gdpbv100::Message);
  void PrintSlcInfo(gdpbv100::Message);
  void PrintSysInfo(gdpbv100::Message);
  void FillPattInfo(gdpbv100::Message);
  void PrintGenInfo(gdpbv100::Message);
  void FillStarTrigInfo(gdpbv100::Message);
 
  inline Int_t GetArrayIndex(Int_t gdpbId, Int_t get4Id) {
    return gdpbId * fuNrOfGet4PerGdpb + get4Id;
  }
 
  std::vector<UInt_t> fvuPadiToGet4;
  std::vector<UInt_t> fvuGet4ToPadi;
 
  static const UInt_t kuNbGet4PerGbtx = 5 * 8;  
  std::vector<UInt_t> fvuElinkToGet4;
  std::vector<UInt_t> fvuGet4ToElink;
  inline UInt_t ConvertElinkToGet4(UInt_t uElinkIdx) {
    return fvuElinkToGet4[uElinkIdx % kuNbGet4PerGbtx]
           + kuNbGet4PerGbtx * (uElinkIdx / kuNbGet4PerGbtx);
  }
  inline UInt_t ConvertGet4ToElink(UInt_t uGet4Idx) {
    return fvuGet4ToElink[uGet4Idx % kuNbGet4PerGbtx]
           + kuNbGet4PerGbtx * (uGet4Idx / kuNbGet4PerGbtx);
  }
 
  std::chrono::time_point<std::chrono::system_clock> fTimeLastHistoSaving;
 
  CbmStar2019MonitorTof(const CbmStar2019MonitorTof&);
  CbmStar2019MonitorTof operator=(const CbmStar2019MonitorTof&);
 
  ClassDef(CbmStar2019MonitorTof, 1)
};
 
#endif