CbmRichDigiQa.h

Go to the documentation of this file.

/* Copyright (C) 2023-2024 UGiessen, Giessen
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Martin Beyer [committer] */
 
#ifndef CBMRICHDIGIQA_H
#define CBMRICHDIGIQA_H
 
#include "FairTask.h"
 
#include <map>
#include <vector>
 
class CbmDigiManager;
class CbmRichDigi;
class CbmHistManager;
class TClonesArray;
 
class CbmRichDigiQa : public FairTask {
 public:
  CbmRichDigiQa() : FairTask("CbmRichDigiQa"){};
 
  ~CbmRichDigiQa() = default;
 
  CbmRichDigiQa(const CbmRichDigiQa&) = delete;
 
  CbmRichDigiQa operator=(const CbmRichDigiQa&) = delete;
 
  InitStatus Init();
 
  void Exec(Option_t* option);
 
  void Finish();
 
  void InitHistograms();
 
  void SetToTLimits(Double_t low, Double_t high)
  {
    fToTLimitLow  = low;
    fToTLimitHigh = high;
  }
 
  void SetNeighbourTimeLimit(Double_t limit) { fNeighbourTimeLimit = limit; }
 
 private:
  int fEventNum{};
 
  CbmDigiManager* fDigiMan{nullptr};
  CbmHistManager* fHM{nullptr};
 
  TClonesArray* fRichPoints{nullptr};
 
  std::map<Int_t, std::vector<Double_t>> fFiredTimes{};  // Store times of fired digis for each pixel in Ts
  std::map<Int_t, std::vector<std::pair<Double_t, Int_t>>>
    fPmtDigisTimeAddress{};  // Store times and addresses of fired digis for each pmt in Ts
 
  Double_t fNeighbourTimeLimit{5.};  // Time in which digi is considered a neighbour
  Double_t fToTLimitLow{-1.};        // Lower limit for ToT cut
  Double_t fToTLimitHigh{-1.};       // Upper limit for ToT cut
 
  ClassDef(CbmRichDigiQa, 1)
};
 
#endif