CbmMcbm2018MonitorMcbmSync.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] */
 
// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                   CbmMcbm2018MonitorMcbmSync                      -----
// -----                Created 11/05/18  by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#ifndef CBMMCBM2018MONITORMCBMSYNC_H
#define CBMMCBM2018MONITORMCBMSYNC_H
 
#include "Timeslice.hpp"
 
// Data
#include "StsXyterFinalHit.h"
#include "StsXyterMessage.h"
#include "gDpbMessv100.h"
 
// CbmRoot
#include "CbmHistManager.h"
#include "CbmMcbmUnpack.h"
 
// C++11
#include <chrono>
 
// C/C++
#include <map>
#include <set>
#include <vector>
 
class CbmMcbm2018StsPar;
 
class CbmMcbm2018TofPar;
 
class CbmMcbm2018MonitorMcbmSync : public CbmMcbmUnpack {
public:
  CbmMcbm2018MonitorMcbmSync();
  virtual ~CbmMcbm2018MonitorMcbmSync();
 
  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);
  void AddMsComponentToListSts(size_t component);
  void AddMsComponentToListTof(size_t component);
 
  virtual void SetNbMsInTs(size_t uCoreMsNb, size_t uOverlapMsNb);
  void SetMsOverlap(size_t uOverlapMsNb = 1) { fuNbOverMsPerTs = uOverlapMsNb; }
  size_t GetMsOverlap() { return fuNbOverMsPerTs; }
  void SetIgnoreMsOverlap(Bool_t bIgnoreOver = kTRUE) { fbIgnoreOverlapMs = bIgnoreOver; }
 
  void ResetAllHistos();
  void SaveAllHistos(TString sFileName = "");
  void SavePulserHistos(TString sFileName = "");
  void SetHistoFileName(TString sFileName = "data/SetupHistos.root") { fsHistoFileFullname = sFileName; }
 
  inline void SetDiamondDpbIdx(UInt_t uIdx = 2) { fuDiamondDpbIdx = uIdx; }
  inline void SetMuchDpbIdx(UInt_t uIdx = 1) { fuMuchDpbIdx = uIdx; }
  inline void SetStsTofOffsetNs(Double_t dOffsetIn = 0.0) { fdStsTofOffsetNs = dOffsetIn; }
  inline void SetMuchTofOffsetNs(Double_t dOffsetIn = 0.0) { fdMuchTofOffsetNs = dOffsetIn; }
  inline void SetUseBestPair(Bool_t bInFLag = kTRUE) { fbUseBestPair = bInFLag; }
  inline void SetTsLevelAna(Bool_t bInFLag = kTRUE) { fbTsLevelAna = bInFLag; }
  inline void SetSpillLimits(Double_t dSpillA, Double_t dSpillB, Double_t dSpillC)
  {
    fdSpillStartA = dSpillA;
    fdSpillStartB = dSpillB;
    fdSpillStartC = dSpillC;
  }
 
  void UseDaqBuffer(Bool_t) {};
 
private:
  std::vector<size_t> fvMsComponentsListSts;  
  std::vector<size_t> fvMsComponentsListTof;  
  size_t fuNbCoreMsPerTs;                     
  size_t fuNbOverMsPerTs;                     
  Bool_t fbIgnoreOverlapMs;          
  CbmMcbm2018StsPar* fUnpackParSts;  
  UInt_t fuStsNrOfDpbs;              
  std::map<UInt_t, UInt_t> fmStsDpbIdIndexMap;  
  UInt_t fuMuchDpbIdx;
  CbmMcbm2018TofPar* fUnpackParTof;  
  UInt_t fuTofNrOfDpbs;              
  std::map<UInt_t, UInt_t> fmTofDpbIdIndexMap;
  UInt_t fuDiamondDpbIdx;
  UInt_t fuTotalNrOfDpb;
  Double_t fdStsTofOffsetNs;
  Double_t fdMuchTofOffsetNs;
  Bool_t fbUseBestPair;
  Bool_t fbTsLevelAna;
 
  // Constants
  static const UInt_t kuStsBytesPerMessage = 4;  
  static const UInt_t kuTofBytesPerMessage = 8;
 
  // Internal Control/status of monitor
  // Histo File name and path
  TString fsHistoFileFullname;
  // Task configuration values
  Bool_t fbPrintMessages;
  stsxyter::MessagePrintMask fPrintMessCtrlSts;
  // 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;     
 
  /****************** STS Sync ******************************************/
  // Data format control
  std::vector<ULong64_t> fvulStsCurrentTsMsb;   
  std::vector<UInt_t> fvuStsCurrentTsMsbCycle;  
 
  // Hits time-sorting
  std::vector<std::vector<stsxyter::FinalHit>>
    fvmStsSdpbHitsInMs;  
  std::vector<std::vector<stsxyter::FinalHit>>
    fvmStsSdpbHitsInTs;  
 
  void CreateStsHistograms();
 
  Bool_t ProcessStsMs(const fles::Timeslice& ts, size_t uMsComp, UInt_t uMsIdx);
 
  void FillStsHitInfo(stsxyter::Message mess, const UInt_t& uMsIdx);
  void FillStsTsMsbInfo(stsxyter::Message mess, UInt_t uMessIdx = 0, UInt_t uMsIdx = 0);
  void FillStsEpochInfo(stsxyter::Message mess);
  /****************** STS Sync ******************************************/
 
  /****************** TOF Sync ******************************************/
 
  uint64_t fulTofCurrentTsIndex;  // Idx of the current TS
  size_t fuTofCurrentMs;          // Idx of the current MS in TS (0 to fuTotalMsNb)
  Double_t fdTofMsIndex;          // Time in ns of current MS from its index
  UInt_t fuTofGdpbId;             // Id (hex number) of the GDPB for current message
  UInt_t fuTofGdpbNr;             // running number (0 to fuNrOfGdpbs) of the GDPB for current message
  Int_t fiTofEquipmentId;
  std::vector<int> fviTofMsgCounter;
 
  std::vector<ULong64_t> fvulTofCurrentEpoch;       
  std::vector<ULong64_t> fvulTofCurrentEpochCycle;  
  std::vector<ULong64_t> fvulTofCurrentEpochFull;   
 
  ULong64_t fulTofCurrentEpochTime; 
  std::vector<std::vector<gdpbv100::Message>> fvmTofEpSupprBuffer;
 
  std::vector<std::vector<gdpbv100::FullMessage>> fvmTofGdpbHitsInMs;
  std::vector<std::vector<gdpbv100::FullMessage>> fvmTofGdpbHitsInTs;
 
  void CreateTofHistograms();
 
  Bool_t ProcessTofMs(const fles::Timeslice& ts, size_t uMsComp, UInt_t uMsIdx);
 
  void FillTofEpochCycle(uint64_t ulCycleData);
  void FillTofHitInfo(gdpbv100::Message);
  void FillTofEpochInfo(gdpbv100::Message);
  /****************** TOF Sync ******************************************/
 
  /****************** mCBM Sync *****************************************/
  TH2* fhMcbmHitsNbPerMs;
 
  TH2* fhMcbmTimeDiffToDiamond;
  TH2* fhMcbmTimeDiffToDiamondWide;
  TH2* fhMcbmTimeDiffToDiamondTs;
 
  TH2* fhMcbmTimeDiffToMuch;
  TH2* fhMcbmTimeDiffToMuchWide;
  TH2* fhMcbmTimeDiffToMuchTs;
 
  TH2* fhMcbmStsTimeDiffToMuchVsAdc;
  TH2* fhMcbmStsTimeDiffToMuchWideVsAdc;
  TH2* fhMcbmStsTimeDiffToMuchTsVsAdc;
 
  std::vector<TH2*> fvhMcbmTimeDiffToDiamondEvoDpb;
  std::vector<TH2*> fvhMcbmTimeDiffToDiamondWideEvoDpb;
  std::vector<TH2*> fvhMcbmTimeDiffToDiamondTsEvoDpb;
 
  Double_t fdSpillStartA;
  Double_t fdSpillStartB;
  Double_t fdSpillStartC;
  std::vector<TH1*> fvhHitsTimeEvoSpillA;
  std::vector<TH1*> fvhHitsTimeEvoSpillB;
 
  std::vector<TH2*> fvhMcbmTimeDiffToDiamondEvoSpillA;
  std::vector<TH2*> fvhMcbmTimeDiffToDiamondEvoSpillB;
 
  std::vector<TH2*> fvhMcbmTimeDiffToMuchEvoSpillA;
  std::vector<TH2*> fvhMcbmTimeDiffToMuchEvoSpillB;
 
  void CreateMcbmHistograms();
  /****************** mCBM Sync *****************************************/
 
 
  CbmMcbm2018MonitorMcbmSync(const CbmMcbm2018MonitorMcbmSync&);
  CbmMcbm2018MonitorMcbmSync operator=(const CbmMcbm2018MonitorMcbmSync&);
 
  ClassDef(CbmMcbm2018MonitorMcbmSync, 1)
};
 
#endif  // CBMMCBM2018MONITORMCBMSYNC_H