CbmStsTimeCal.h

Go to the documentation of this file.

/* Copyright (C) 2023 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Dario Ramirez [committer] */
 
#ifndef CBMSTSTIMECAL_H
#define CBMSTSTIMECAL_H
 
#include "CbmDefs.h"
#include "CbmDigiManager.h"
#include "CbmStsAnaBase.h"
#include "CbmStsUtils.h"
#include "CbmYaml.h"
#include "algo/detectors/sts/WalkMap.h"
 
#include <FairTask.h>
 
#include <TClonesArray.h>
 
#include <cstring>
#include <map>
#include <unordered_set>


class CbmStsTimeCal : public FairTask, public CbmStsAnaBase {
 public:
  CbmStsTimeCal()  = default;
  ~CbmStsTimeCal() = default;

  CbmStsTimeCal(ECbmModuleId system, double lower_lim = -60, double upper_lim = +60);

  CbmStsTimeCal(int run_id, ECbmModuleId system, double lower_lim, double upper_lim);

  void SetWalkFile(std::string par_file);
 
  InitStatus Init();
  void Exec(Option_t*);
  void Finish();
 
 
 private:
  double fGlobalTimeOffset{0};
  double fTimeWindowMin{-60};
  double fTimeWindowMax{+60};
  const double kStsClock{3.125};
  std::string fParFile{""};
  std::string fOutputPath{""};
 
  std::map<int /* Module */, std::map<int /* ASIC */, std::vector<double> /* ADC offset */>> fWalkMapRaw;
  cbm::algo::sts::WalkMap fWalkMap;
 
  std::unique_ptr<TFile> fReportOffset;
  std::unique_ptr<TFile> fReportFit;
 
  ECbmModuleId fRefSystem{ECbmModuleId::kBmon};
 
  CbmDigiManager* fDigiManager{nullptr};
 
  TClonesArray* fCbmEvtArray{nullptr};

  void BookHistograms(int32_t address);

  void CheckTimeWalk();

  double FindModuleOffset(int32_t);

  double FindGlobalOffset();

  template<class Digi>
  void CheckTiming();

  void InitTimeWalkMap(int32_t address);

  void LoadWalkFromFile();

  void WriteYAML();
 
  void DrawResults();
 
  ClassDef(CbmStsTimeCal, 1);
};
#endif