CbmMcbm2018MonitorAlgoTofPulser.h

Go to the documentation of this file.

/* Copyright (C) 2019-2020 Facility for Antiproton and Ion Research in Europe, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pierre-Alain Loizeau [committer] */
 
// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                  CbmMcbm2018MonitorAlgoTofPulser                     -----
// -----               Created 12.10.2019 by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#ifndef CbmMcbm2018MonitorAlgoTofPulser_H
#define CbmMcbm2018MonitorAlgoTofPulser_H
 
#include "CbmStar2019Algo.h"
 
// Data
#include "CbmTofDigi.h"
 
#include "gDpbMessv100.h"
 
// CbmRoot
 
// C++11
#include <chrono>
 
// C/C++
#include <map>
#include <vector>
 
class CbmMcbm2018TofPar;
class TH1;
class TH2;
class TProfile;
 
class CbmMcbm2018MonitorAlgoTofPulser : public CbmStar2019Algo<CbmTofDigi> {
public:
  CbmMcbm2018MonitorAlgoTofPulser();
  ~CbmMcbm2018MonitorAlgoTofPulser();
 
  virtual Bool_t Init();
  virtual void Reset();
  virtual void Finish();
 
  Bool_t InitContainers();
  Bool_t ReInitContainers();
  TList* GetParList();
 
  Bool_t InitParameters();
 
  Bool_t ProcessTs(const fles::Timeslice& ts);
  Bool_t ProcessTs(const fles::Timeslice& ts, size_t /*component*/) { return ProcessTs(ts); }
  Bool_t ProcessMs(const fles::Timeslice& ts, size_t uMsCompIdx, size_t uMsIdx);
 
  void AddMsComponentToList(size_t component, UShort_t usDetectorId);
 
  Bool_t CreateHistograms();
  Bool_t FillHistograms();
  Bool_t UpdateStats();
  Bool_t ResetHistograms();
 
  inline void SetGdpbIndex(Int_t iGdpb = -1) { fiGdpbIndex = iGdpb; }
  inline void SetUpdateFreqTs(UInt_t uFreq = 100) { fuUpdateFreqTs = uFreq; }
  inline void SetPulserTotLimits(UInt_t uMin, UInt_t uMax)
  {
    fuPulserMinTot = uMin;
    fuPulserMaxTot = uMax;
  }
  inline void SetPulserChannel(UInt_t uChan) { fuPulserChannel = uChan; }
  inline void SetHistoryHistoSize(UInt_t inHistorySizeSec = 1800) { fuHistoryHistoSize = inHistorySizeSec; }
 
private:
  std::vector<Bool_t> fvbMaskedComponents;
  Int_t fiGdpbIndex;
  UInt_t fuUpdateFreqTs;
 
  CbmMcbm2018TofPar* fUnpackPar;             
  UInt_t fuNrOfGdpbs;                        
  std::map<UInt_t, UInt_t> fGdpbIdIndexMap;  
  UInt_t fuNrOfFeePerGdpb;                   
  UInt_t fuNrOfGet4PerFee;                   
  UInt_t fuNrOfChannelsPerGet4;              
  UInt_t fuNrOfChannelsPerFee;               
  UInt_t fuNrOfGet4;                         
  UInt_t fuNrOfGet4PerGdpb;                  
  UInt_t fuNrOfChannelsPerGdpb;              
 
  UInt_t fuPulserMinTot;
  UInt_t fuPulserMaxTot;
  UInt_t fuPulserChannel;
 
  static const Int_t kiMaxNbFlibLinks   = 32;
  static const UInt_t kuBytesPerMessage = 8;
 
  ULong64_t fulCurrentTsIdx;
  ULong64_t fulCurrentMsIdx;
  uint32_t fuCurrentMsSysId;    
  Double_t fdTsStartTime;       
  Double_t fdTsStopTimeCore;    
  Double_t fdMsTime;            
  UInt_t fuMsIndex;             
  UInt_t fuCurrentEquipmentId;  
  UInt_t fuCurrDpbId;           
  UInt_t fuCurrDpbIdx;          
  UInt_t fuGet4Id;  
  UInt_t fuGet4Nr;  
  std::vector<ULong64_t> fvulCurrentEpoch;       
  std::vector<ULong64_t> fvulCurrentEpochCycle;  
  std::vector<ULong64_t> fvulCurrentEpochFull;   
 
  std::vector<gdpbv100::Message> fvmEpSupprBuffer;
 
  std::vector<std::vector<Bool_t>> fvvbFeeHitFound;  
  std::vector<std::vector<Double_t>> fvvdFeeHits;    
 
  const Double_t kdMaxDtPulserPs = 300e3;
  const UInt_t kuNbBinsDt        = 2000;
  Double_t dMinDt;
  Double_t dMaxDt;
  const Double_t kdFitZoomWidthPs = 500.0;
  Double_t fdStartTime;       
  UInt_t fuHistoryHistoSize;  
 
  std::vector<std::vector<TH1*>> fvvhFeePairPulserTimeDiff;  
  TH2* fhPulserTimeDiffMean;                                 
  TH2* fhPulserTimeDiffRms;                                  
  TH2* fhPulserTimeDiffRmsZoom;                              
  TH2* fhPulserRmsGdpbToRefEvo;                              
  TH2* fhPulserRmsGbtxToRefEvo;                              
 
  TCanvas* fcSummary;  
 
  void ProcessEpochCycle(uint64_t ulCycleData);
  void ProcessEpoch(gdpbv100::Message mess);
 
  void ProcessEpSupprBuffer();
 
  void ProcessHit(gdpbv100::FullMessage mess);
 
  CbmMcbm2018MonitorAlgoTofPulser(const CbmMcbm2018MonitorAlgoTofPulser&);
  CbmMcbm2018MonitorAlgoTofPulser operator=(const CbmMcbm2018MonitorAlgoTofPulser&);
 
  ClassDef(CbmMcbm2018MonitorAlgoTofPulser, 1)
};
 
#endif  // CbmMcbm2018MonitorAlgoTofPulser_H