CbmMcbmCheckTimingAlgo.h

Go to the documentation of this file.

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *              GNU Lesser General Public Licence (LGPL) version 3,             *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
#ifndef CBMMCBMCHECKTIMINGALGO_H
#define CBMMCBMCHECKTIMINGALGO_H
 
#include "CbmModuleList.h"
#include "CbmTofDigi.h"
#include "TString.h"
#include <vector>
 
class TH1;
class TH2;
class CbmDigiManager;
 
class CheckTimingDetector {
public:
  CheckTimingDetector() { ; }
  CheckTimingDetector(ECbmModuleId detIdIn, std::string sNameIn) {
    detId = detIdIn;
    sName = sNameIn;
  }
 
  ECbmModuleId detId     = ECbmModuleId::kNotExist;
  std::string sName      = "Invalid";
  Double_t dTimeRangeBeg = -1000.0;
  Double_t dTimeRangeEnd = 1000.0;
  UInt_t uRangeNbBins    = 320;
  UInt_t uChargeCutMin =
    0;  
  UInt_t uChargeCutMax =
    0;  
 
  Double_t dPrevTime      = 0.;
  Int_t iPrevRefFirstDigi = 0;
};
 
class CbmMcbmCheckTimingAlgo {
public:
  CbmMcbmCheckTimingAlgo();
 
  CbmMcbmCheckTimingAlgo(const CbmMcbmCheckTimingAlgo&) = delete;
  CbmMcbmCheckTimingAlgo operator=(const CbmMcbmCheckTimingAlgo&) = delete;
 
  //  CbmMcbmCheckTimingAlgo(Int_t verbose);
 
 
  ~CbmMcbmCheckTimingAlgo();
 
 
  Bool_t Init();
 
  Bool_t ReInit();
 
 
  void ProcessTs();
 
  void SetParContainers();
 
  void Finish();
 
  inline void SetOutFilename(TString sNameIn) { fOutFileName = sNameIn; }
  void WriteHistos();
 
  void SetReferenceDetector(ECbmModuleId refDetIn,
                            std::string sNameIn,
                            Double_t dTimeRangeBegIn = -1000.0,
                            Double_t dTimeRangeEndIn = 1000.0,
                            UInt_t uRangeNbBinsIn    = 320,
                            UInt_t uChargeCutMinIn   = 0,
                            UInt_t uChargeCutMaxIn   = 0);
  void AddCheckDetector(ECbmModuleId detIn,
                        std::string sNameIn,
                        Double_t dTimeRangeBegIn = -1000.0,
                        Double_t dTimeRangeEndIn = 1000.0,
                        UInt_t uRangeNbBinsIn    = 320,
                        UInt_t uChargeCutMinIn   = 0,
                        UInt_t uChargeCutMaxIn   = 0);
  void RemoveCheckDetector(ECbmModuleId detIn);
 
private:
  void CheckDataPresence(CheckTimingDetector detToCheck);
  void CreateHistos();
 
  template<class DigiRef>
  void CheckInterSystemOffset();
  template<class Digi>
  void FillTimeOffsetHistos(const Double_t dRefTime,
                            const Double_t dRefCharge,
                            UInt_t uDetIdx);
 
 
  CbmDigiManager* fDigiMan = nullptr;  
 
  const std::vector<CbmTofDigi>* fpT0DigiVec = nullptr;  
 
  //
  UInt_t fuNbTs = 0;
 
  CheckTimingDetector fRefDet {CheckTimingDetector(ECbmModuleId::kT0, "T0")};
  std::vector<CheckTimingDetector> fvDets {
    CheckTimingDetector(ECbmModuleId::kSts, "Sts"),
    CheckTimingDetector(ECbmModuleId::kMuch, "Much"),
    CheckTimingDetector(ECbmModuleId::kTrd, "Trd"),
    CheckTimingDetector(ECbmModuleId::kTof, "Tof"),
    CheckTimingDetector(ECbmModuleId::kRich, "Rich"),
    CheckTimingDetector(ECbmModuleId::kPsd, "Psd")};
 
  std::vector<TH1*> fvhDetSelfDiff           = {};
  std::vector<TH1*> fvhDetToRefDiff          = {};
  std::vector<TH2*> fvhDetToRefDiffRefCharge = {};
  std::vector<TH2*> fvhDetToRefDiffDetCharge = {};
  std::vector<TH2*> fvhDetToRefDiffEvo       = {};
  std::vector<TH2*> fvhDetToRefDiffEvoLong   = {};
 
  TString fOutFileName = "data/HistosCheckTiming.root";
 
  ClassDefNV(CbmMcbmCheckTimingAlgo, 1);
};
 
#endif