CbmCheckDigisNbCorr.h

Go to the documentation of this file.

/* Copyright (C) 2019-2021 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Pierre-Alain Loizeau [committer] */
 
#ifndef CBMCHECKDIGISNBCORR_H
#define CBMCHECKDIGISNBCORR_H
 
#include "CbmTofDigi.h"
 
#include "FairTask.h"
 
#include "TString.h"
 
#include <vector>
 
class TClonesArray;
class TH1;
class TH2;
class TProfile;
class CbmDigiManager;
 
class CbmCheckDigisNbCorr : public FairTask {
public:
  CbmCheckDigisNbCorr();
 
  CbmCheckDigisNbCorr(const CbmCheckDigisNbCorr&) = delete;
  CbmCheckDigisNbCorr operator=(const CbmCheckDigisNbCorr&) = delete;
 
  //  CbmCheckDigisNbCorr(Int_t verbose);
 
 
  ~CbmCheckDigisNbCorr();
 
 
  virtual InitStatus Init();
 
  virtual InitStatus ReInit();
 
 
  virtual void Exec(Option_t*);
 
  virtual void SetParContainers();
 
  virtual void Finish();
 
  void SetStsOffset(Double_t dval = 1600) { fdStsOffset = dval; }
 
  void SetMuchOffset(Double_t dval = 1600) { fdMuchOffset = dval; }
 
  void SetTofOffset(Double_t dval = 50) { fdTofOffset = dval; }
 
  void SetRichOffset(Double_t dval = 1200) { fdRichOffset = dval; }
 
  inline void SetTimeInterval(Double_t dInterNs, Double_t dTsInNs)
  {
    fdBinWidthNs = dInterNs;
    fdTsLengthNs = dTsInNs;
  }
 
  inline void SetBmonPulserTotLimits(UInt_t uMin, UInt_t uMax)
  {
    fuMinTotPulserBmon = uMin;
    fuMaxTotPulserBmon = uMax;
  }
 
  inline void SetOutFilename(TString sNameIn) { fOutFileName = sNameIn; }
 
private:
  void CalcNrBins();
  void CreateHistos();
  void WriteHistos();
 
 
  TClonesArray* fTsMetaData = nullptr;
  CbmDigiManager* fDigiMan  = nullptr;  
  const std::vector<CbmTofDigi>* fBmonDigiVec = nullptr;  
  TClonesArray* fBmonDigiArr {nullptr};                   
 
  UInt_t fuMinTotPulserBmon;
  UInt_t fuMaxTotPulserBmon;
 
  //
  Int_t fNrTs = 0;
 
  Int_t fNrOfBmonErrors = 0;
  Int_t fNrOfBmonDigis  = 0;
  Int_t fNrOfStsErrors  = 0;
  Int_t fNrOfStsDigis   = 0;
  Int_t fNrOfMuchErrors = 0;
  Int_t fNrOfMuchDigis  = 0;
  Int_t fNrOfTofErrors  = 0;
  Int_t fNrOfTofDigis   = 0;
  Int_t fNrOfRichErrors = 0;
  Int_t fNrOfRichDigis  = 0;
 
  Double_t fdStsOffset  = 0.0;
  Double_t fdMuchOffset = 0.0;
  Double_t fdTofOffset  = 0.0;
  Double_t fdRichOffset = 0.0;
 
  Double_t fdTsLengthNs    = 100 * 102e3;
  Double_t fdBinWidthNs    = 10e3;
  Int_t fiBinNb            = fdTsLengthNs / fdBinWidthNs;
  Double_t fdTsStartOffset = 0;
 
  std::vector<UInt_t> fvuNbDigisPerBinBmon;
  std::vector<UInt_t> fvuNbDigisPerBinSts;
  std::vector<UInt_t> fvuNbDigisPerBinMuch;
  std::vector<UInt_t> fvuNbDigisPerBinTof;
  std::vector<UInt_t> fvuNbDigisPerBinRich;
 
  TH2* fBmonStsCorr  = nullptr;
  TH2* fBmonMuchCorr = nullptr;
  TH2* fBmonTofCorr  = nullptr;
  TH2* fBmonRichCorr = nullptr;
 
  TH2* fStsMuchCorr = nullptr;
  TH2* fStsTofCorr  = nullptr;
  TH2* fStsRichCorr = nullptr;
 
  TH2* fMuchTofCorr  = nullptr;
  TH2* fMuchRichCorr = nullptr;
 
  TH2* fTofRichCorr = nullptr;
 
  TProfile* fBmonStsCorrProf  = nullptr;
  TProfile* fBmonMuchCorrProf = nullptr;
  TProfile* fBmonTofCorrProf  = nullptr;
  TProfile* fBmonRichCorrProf = nullptr;
 
  TProfile* fStsMuchCorrProf = nullptr;
  TProfile* fStsTofCorrProf  = nullptr;
  TProfile* fStsRichCorrProf = nullptr;
 
  TProfile* fMuchTofCorrProf  = nullptr;
  TProfile* fMuchRichCorrProf = nullptr;
 
  TProfile* fTofRichCorrProf = nullptr;
 
  static const UInt_t kuMaxNbStsDpbs = 2;
  std::vector<std::vector<UInt_t>> fvuNbDigisPerBinStsDpb;
  TH2* fBmonStsDpbCorr[kuMaxNbStsDpbs];
  TH2* fStsMuchDpbCorr[kuMaxNbStsDpbs];
  TH2* fStsTofDpbCorr[kuMaxNbStsDpbs];
  TH2* fStsRichDpbCorr[kuMaxNbStsDpbs];
  TProfile* fBmonStsDpbCorrProf[kuMaxNbStsDpbs];
  TProfile* fStsMuchDpbCorrProf[kuMaxNbStsDpbs];
  TProfile* fStsTofDpbCorrProf[kuMaxNbStsDpbs];
  TProfile* fStsRichDpbCorrProf[kuMaxNbStsDpbs];
  /*
    static const UInt_t kuMaxNbMuchAsics = 24;
    TH2*      fBmonMuchAsicCorr[kuMaxNbMuchAsics];
    TProfile* fBmonMuchAsicCorrProf[kuMaxNbMuchAsics];
*/
  TString fOutFileName {"HistosCheckDigisNbCorr.root"};
 
  ClassDef(CbmCheckDigisNbCorr, 1);
};
 
#endif