CbmMcbm2018MonitorStsSync.h

Go to the documentation of this file.

/* Copyright (C) 2018-2019 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pierre-Alain Loizeau [committer] */
 
// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                   CbmMcbm2018MonitorStsSync                        -----
// -----                Created 11/05/18  by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#ifndef CBMMCBM2018MONITORSTSSYNC_H
#define CBMMCBM2018MONITORSTSSYNC_H
 
#include "Timeslice.hpp"
 
// Data
#include "StsXyterFinalHit.h"
#include "StsXyterMessage.h"
 
// CbmRoot
#include "CbmHistManager.h"
#include "CbmMcbmUnpack.h"
 
// C++11
#include <chrono>
 
// C/C++
#include <map>
#include <set>
#include <vector>
 
class CbmCern2017UnpackParHodo;
class CbmCern2017UnpackParSts;
 
class CbmMcbm2018MonitorStsSync : public CbmMcbmUnpack {
public:
  CbmMcbm2018MonitorStsSync();
  virtual ~CbmMcbm2018MonitorStsSync();
 
  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();
 
  virtual void AddMsComponentToList(size_t component, UShort_t usDetectorId);
  virtual void SetNbMsInTs(size_t uCoreMsNb, size_t uOverlapMsNb);
  void SetMsOverlap(size_t uOverlapMsNb = 1) { fuNbOverMsPerTs = uOverlapMsNb; }
  size_t GetMsOverlap() { return fuNbOverMsPerTs; }
 
  void SetRunStart(Int_t dateIn, Int_t timeIn, Int_t iBinSize = 5);
 
  void ResetAllHistos();
  void SaveAllHistos(TString sFileName = "");
  void SavePulserHistos(TString sFileName = "");
  void SetHistoFileName(TString sFileName = "data/SetupHistos.root") { fsHistoFileFullname = sFileName; }
 
  void SetPrintMessage(Bool_t bPrintMessOn             = kTRUE,
                       stsxyter::MessagePrintMask ctrl = stsxyter::MessagePrintMask::msg_print_Hex
                                                         | stsxyter::MessagePrintMask::msg_print_Human)
  {
    fbPrintMessages = bPrintMessOn;
    fPrintMessCtrl  = ctrl;
  }
  void SetLongDurationLimits(UInt_t uDurationSeconds = 3600, UInt_t uBinSize = 1);
  void SetCoincidenceBorder(Double_t dCenterPos, Double_t dBorderVal);
 
  void UseDaqBuffer(Bool_t) {};
 
private:
  // Parameters
  // FLES containers
  std::vector<size_t> fvMsComponentsList;    
  size_t fuNbCoreMsPerTs;                    
  size_t fuNbOverMsPerTs;                    
  Bool_t fbIgnoreOverlapMs;                  
                                             // Unpacking and mapping
  CbmCern2017UnpackParHodo* fUnpackParHodo;  
  UInt_t fuNrOfDpbs;                         
  std::map<UInt_t, UInt_t> fDpbIdIndexMap;   
  UInt_t fuNbElinksPerDpb;                   
  UInt_t fuNbStsXyters;                      
  UInt_t fuNbChanPerAsic;  
  std::vector<std::vector<UInt_t>>
    fvuElinkToAsic;  
 
  // Constants
  static const Int_t kiMaxNbFlibLinks   = 16;
  static const UInt_t kuBytesPerMessage = 4;
 
  // Internal Control/status of monitor
  // Histo File name and path
  TString fsHistoFileFullname;
  // Task configuration values
  Bool_t fbPrintMessages;
  stsxyter::MessagePrintMask fPrintMessCtrl;
  // TS/MS info
  ULong64_t fulCurrentTsIdx;
  ULong64_t fulCurrentMsIdx;
  // Current data properties
  std::map<stsxyter::MessType, UInt_t> fmMsgCounter;
  UInt_t fuCurrentEquipmentId;  
  UInt_t fuCurrDpbId;           
  UInt_t fuCurrDpbIdx;          
  Int_t fiRunStartDateTimeSec;  
  Int_t fiBinSizeDatePlots;     
 
  // Data format control
  std::vector<ULong64_t> fvulCurrentTsMsb;                  
  std::vector<UInt_t> fvuCurrentTsMsbCycle;                 
  std::vector<UInt_t> fvuInitialHeaderDone;                 
  std::vector<UInt_t> fvuInitialTsMsbCycleHeader;           
  std::vector<UInt_t> fvuElinkLastTsHit;                    
                                                            // Hits comparison
  std::vector<std::vector<ULong64_t>> fvulChanLastHitTime;  
  std::vector<std::vector<Double_t>> fvdChanLastHitTime;    
  std::vector<Double_t> fvdMsTime;                          
  std::vector<std::vector<std::vector<UInt_t>>> fvuChanNbHitsInMs;  
  std::vector<std::vector<std::vector<Double_t>>>
    fvdChanLastHitTimeInMs;  
  std::vector<std::vector<std::vector<UShort_t>>>
    fvusChanLastHitAdcInMs;  
  //   std::vector< std::vector< std::multiset< stsxyter::FinalHit > > > fvmChanHitsInTs; //! All hits (time & ADC) in bins in last TS for each Channel
  // Starting state book-keeping
  Double_t fdStartTime;     
  Double_t fdStartTimeMsSz; 
  std::chrono::steady_clock::time_point
    ftStartTimeUnix; 
  // Hits time-sorting
  std::vector<stsxyter::FinalHit>
    fvmHitsInMs;  
  std::vector<std::vector<stsxyter::FinalHit>>
    fvmAsicHitsInMs;  
                      // Coincidence histos
  UInt_t fuMaxNbMicroslices;
  // Rate evolution histos
  Bool_t fbLongHistoEnable;
  UInt_t fuLongHistoNbSeconds;
  UInt_t fuLongHistoBinSizeSec;
  UInt_t fuLongHistoBinNb;
 
  // Coincidences in sorted hits
  Double_t fdCoincCenter;  // ns
  Double_t fdCoincBorder;  // ns, +/-
  Double_t fdCoincMin;     // ns
  Double_t fdCoincMax;     // ns
 
  // Histograms
  CbmHistManager* fHM;  
  TH1* fhPulserMessType;
  TH1* fhPulserSysMessType;
  TH2* fhPulserMessTypePerDpb;
  TH2* fhPulserSysMessTypePerDpb;
  /*
   TH2* fhStsDpbRawTsMsb;
   TH2* fhStsDpbRawTsMsbSx;
   TH2* fhStsDpbRawTsMsbDpb;
*/
  TH2* fhPulserMessTypePerElink;
  TH2* fhPulserSysMessTypePerElink;
  TH2* fhPulserStatusMessType;
  TH2* fhPulserMsStatusFieldType;
 
  std::vector<TH1*> fhPulserChanCntRaw;           
  std::vector<TH1*> fhPulserChanCntRawGood;       
  std::vector<TH2*> fhPulserChanAdcRaw;           
  std::vector<TProfile*> fhPulserChanAdcRawProf;  
  std::vector<TH2*> fhPulserChanRawTs;            
  std::vector<TH2*> fhPulserChanMissEvt;          
  std::vector<TH2*> fhPulserChanMissEvtEvo;       
  std::vector<TH2*> fhPulserChanHitRateEvo;       
  std::vector<TH1*> fhPulserFebRateEvo;           
  std::vector<TH1*> fhPulserFebMissEvtEvo;        
  std::vector<TH2*> fhPulserChanHitRateEvoLong;   
  std::vector<TH1*> fhPulserFebRateEvoLong;       
 
  TCanvas* fcMsSizeAll;
  TH1* fhMsSz[kiMaxNbFlibLinks];
  TProfile* fhMsSzTime[kiMaxNbFlibLinks];
 
  Double_t fdStartTs;                                                         
  std::vector<stsxyter::FinalHit> fvmLastHitAsic;                             
  std::vector<TH2*> fhPulserTimeDiffPerAsic;                                  
  std::vector<std::vector<TH1*>> fhPulserTimeDiffPerAsicPair;                 
  std::vector<std::vector<TH1*>> fhPulserTimeDiffClkPerAsicPair;              
  std::vector<std::vector<TH2*>> fhPulserTimeDiffEvoPerAsicPair;              
  std::vector<std::vector<TProfile*>> fhPulserTimeDiffEvoPerAsicPairProf;     
  std::vector<std::vector<TProfile*>> fhPulserRawTimeDiffEvoPerAsicPairProf;  
  std::vector<std::vector<TH2*>> fhPulserTsLsbMatchPerAsicPair;               
  std::vector<std::vector<TH2*>> fhPulserTsMsbMatchPerAsicPair;               
  std::vector<std::vector<TProfile*>> fhPulserTsLsbDiffEvoPerAsicPairProf;    
  std::vector<std::vector<TProfile*>> fhPulserTsMsbDiffEvoPerAsicPairProf;    
  std::vector<TH1*> fhPulserIntervalAsic;                                     
  std::vector<TH1*> fhPulserIntervalLongAsic;                                 
 
  static const UInt_t kuNbValuesForTimeDiffMean        = 100;
  constexpr static const double kdPulserPeriod         = 154000;  
  constexpr static const double kdTimeDiffToMeanMargin = 50;
  std::vector<std::vector<std::vector<Double_t>>> fvdLastTimeDiffValuesAsicPair;  
  std::vector<std::vector<UInt_t>> fvuLastTimeDiffSlotAsicPair;                   
  std::vector<std::vector<Double_t>> fvdMeanTimeDiffAsicPair;                     
 
  void UpdatePairMeanValue(UInt_t uAsicA, UInt_t uAsicB, Double_t dNewValue);
 
  void CreateHistograms();
 
  void FillHitInfo(stsxyter::Message mess, const UShort_t& usElinkIdx, const UInt_t& uAsicIdx, const UInt_t& uMsIdx);
  void FillTsMsbInfo(stsxyter::Message mess, UInt_t uMessIdx = 0, UInt_t uMsIdx = 0);
  void FillEpochInfo(stsxyter::Message mess);
 
  CbmMcbm2018MonitorStsSync(const CbmMcbm2018MonitorStsSync&);
  CbmMcbm2018MonitorStsSync operator=(const CbmMcbm2018MonitorStsSync&);
 
  ClassDef(CbmMcbm2018MonitorStsSync, 1)
};
 
#endif  // CBMMCBM2018MONITORSTSSYNC_H