CbmMcbm2018UnpackerAlgoSts.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] */
 
// -----------------------------------------------------------------------------
// -----                                                                   -----
// -----                  CbmMcbm2018UnpackerAlgoSts                       -----
// -----               Created 26.01.2019 by P.-A. Loizeau                 -----
// -----                                                                   -----
// -----------------------------------------------------------------------------
 
#ifndef CbmMcbm2018UnpackerAlgoSts_H
#define CbmMcbm2018UnpackerAlgoSts_H
 
#include "CbmStar2019Algo.h"
 
// Data
#include "CbmStsDigi.h"
 
#include "StsXyterFinalHit.h"
#include "StsXyterMessage.h"
 
// CbmRoot
 
// C++11
#include <chrono>
 
// C/C++
#include <map>
#include <vector>
 
class CbmMcbm2018StsPar;
/*
class TCanvas;
class THttpServer;
*/
class TH1;
class TH2;
class TProfile;
class TProfile2D;
 
class CbmMcbm2018UnpackerAlgoSts : public CbmStar2019Algo<CbmStsDigi> {
public:
  CbmMcbm2018UnpackerAlgoSts();
  ~CbmMcbm2018UnpackerAlgoSts();
 
  virtual Bool_t Init();
  virtual void Reset();
  virtual void Finish();
 
  Bool_t InitContainers();
  Bool_t ReInitContainers();
  TList* GetParList();
 
  Bool_t InitParameters();
 
  std::vector<CbmStsDigi>& GetPulserVector() { return fPulserDigiVect; }
  void ClearPulserVector() { fPulserDigiVect.clear(); }
 
  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 ResetHistograms();
 
  inline void SetMonitorMode(Bool_t bFlagIn = kTRUE) { fbMonitorMode = bFlagIn; }
  inline void SetTimeOffsetNs(Double_t dOffsetIn = 0.0) { fdTimeOffsetNs = dOffsetIn; }
  void SetTimeOffsetNsAsic(UInt_t uAsicIdx, Double_t dOffsetIn = 0.0);
  void MaskNoisyChannel(UInt_t uFeb, UInt_t uChan, Bool_t bMasked = kTRUE);
  void SetAdcCut(UInt_t uAdc) { fdAdcCut = uAdc; }
 
  void SetBinningFwFlag(Bool_t bEnable = kTRUE) { fbBinningFw = bEnable; }
 
  inline void SetVectCapInc(Double_t dIncFact) { fdCapacityIncFactor = dIncFact; }
 
  void SeparatePulserOutput(Bool_t bFlagIn) { fbPulserOutput = bFlagIn; }
 
private:
  Bool_t fbMonitorMode;       
  Bool_t fbDebugMonitorMode;  
  std::vector<Bool_t> fvbMaskedComponents;
  Bool_t fbBinningFw = kFALSE;
  Bool_t fbPulserOutput = kTRUE;  
 
  CbmMcbm2018StsPar* fUnpackPar;  
  UInt_t fuNbModules;             
  std::vector<Int_t>
    fviModuleType;  
  std::vector<Int_t> fviModAddress;         
  UInt_t fuNrOfDpbs;                        
  std::map<UInt_t, UInt_t> fDpbIdIndexMap;  
  std::vector<std::vector<Bool_t>>
    fvbCrobActiveFlag;   
  UInt_t fuNbFebs;       
  UInt_t fuNbStsXyters;  
  std::vector<std::vector<std::vector<Int_t>>>
    fviFebModuleIdx;  
  std::vector<std::vector<std::vector<Int_t>>>
    fviFebModuleSide;  
  std::vector<std::vector<std::vector<Int_t>>>
    fviFebType;  
  std::vector<bool> fvbFebPulser = {};  
  std::vector<Int_t> fviFebAddress;     
  std::vector<Int_t> fviFebSide;        
  std::vector<Double_t> fvdFebAdcGain;  
  std::vector<Double_t> fvdFebAdcOffs;  
 
  Double_t fdTimeOffsetNs;
  std::vector<Double_t> fvdTimeOffsetNsAsics;
  Bool_t fbUseChannelMask;
  std::vector<std::vector<bool>>
    fvvbMaskedChannels;  
  UInt_t fdAdcCut;
 
  static const Int_t kiMaxNbFlibLinks = 32;
 
  ULong64_t fulCurrentTsIdx;
  ULong64_t fulCurrentMsIdx;
  Double_t fdTsStartTime;     
  Double_t fdTsStopTimeCore;  
  Double_t fdMsTime;          
  UInt_t fuMsIndex;           
  std::map<stsxyter::MessType, UInt_t> fmMsgCounter;
  UInt_t fuCurrentEquipmentId;              
  UInt_t fuCurrDpbId;                       
  UInt_t fuCurrDpbIdx;                      
  Int_t fiRunStartDateTimeSec;              
  Int_t fiBinSizeDatePlots;                 
  std::vector<ULong64_t> fvulCurrentTsMsb;  
  std::vector<UInt_t> fvuCurrentTsMsbCycle;  
  Double_t fdStartTime;     
  Double_t fdStartTimeMsSz; 
  std::chrono::steady_clock::time_point
    ftStartTimeUnix; 
  std::vector<stsxyter::FinalHit>
    fvmHitsInMs;  
 
  std::vector<CbmStsDigi> fPulserDigiVect = {};
 
  static const UInt_t kuMaxTsMsbDiffDuplicates = 8;
  std::vector<std::vector<UShort_t>> fvvusLastTsChan;  
  std::vector<std::vector<UShort_t>>
    fvvusLastAdcChan;  
  std::vector<std::vector<UShort_t>>
    fvvusLastTsMsbChan;  
  std::vector<std::vector<UShort_t>>
    fvvusLastTsMsbCycleChan;  
 
  TH1* fhDigisTimeInRun;  
                          /*
      std::vector< TH1* > fvhHitsTimeToTriggerRaw;       //! [sector]
      std::vector< TH1* > fvhMessDistributionInMs;       //! [sector], extra monitor for debug
      TH1 *               fhEventNbPerTs;                //!
      TCanvas *           fcTimeToTrigRaw;               //! All sectors
*/
  TH1* fhVectorSize            = nullptr;
  TH1* fhVectorCapacity        = nullptr;
  size_t fuTsMaxVectorSize     = 0;
  Double_t fdCapacityIncFactor = 1.1;
 
 
  TH1* fhMsCntEvo    = nullptr;
  TH2* fhMsErrorsEvo = nullptr;
 
  /*
        TProfile2D* fhPulserVsTsAB   = nullptr;
      UShort_t fusTsA = 0.0;
      Double_t fdTimeA = 0.0;
*/
  void ProcessHitInfo(const stsxyter::Message& mess, const UShort_t& usElinkIdx, const UInt_t& uAsicIdx,
                      const UInt_t& uMsIdx);
  void ProcessTsMsbInfo(const stsxyter::Message& mess, UInt_t uMessIdx = 0, UInt_t uMsIdx = 0);
  void ProcessEpochInfo(const stsxyter::Message& mess);
  void ProcessStatusInfo(const stsxyter::Message& mess, const UInt_t& uAsicIdx);
 
  CbmMcbm2018UnpackerAlgoSts(const CbmMcbm2018UnpackerAlgoSts&);
  CbmMcbm2018UnpackerAlgoSts operator=(const CbmMcbm2018UnpackerAlgoSts&);
 
  ClassDef(CbmMcbm2018UnpackerAlgoSts, 1)
};
 
#endif