CbmRichUrqmdTest.h

Go to the documentation of this file.

/* Copyright (C) 2012-2024 UGiessen/JINR-LIT, Giessen/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Andrey Lebedev, Semen Lebedev [committer], Martin Beyer */
 
#ifndef CBM_RICH_URQMD_TEST
#define CBM_RICH_URQMD_TEST
 
#include "FairTask.h"
class TClonesArray;
class CbmRichRing;
class TCanvas;
class CbmHistManager;
class CbmMCDataArray;
class CbmMCEventList;
class CbmDigiManager;
class CbmLink;
 
#include <map>
#include <vector>
 
 
class CbmRichUrqmdTest : public FairTask {
 
 public:
  CbmRichUrqmdTest();
 
  virtual ~CbmRichUrqmdTest();
 
  virtual InitStatus Init();
 
  virtual void Exec(Option_t* option);
 
  virtual void Finish();
 
 
  void SetOutputDir(const std::string& dir) { fOutputDir = dir; }
 
 
 private:
  CbmHistManager* fHM = nullptr;
 
  std::string fOutputDir = "";  // output dir for results
 
  TClonesArray* fRichHits        = nullptr;
  TClonesArray* fRichRings       = nullptr;
  CbmMCDataArray* fRichPoints    = nullptr;
  CbmMCDataArray* fMcTracks      = nullptr;
  TClonesArray* fRichRingMatches = nullptr;
  TClonesArray* fRichProjections = nullptr;
  CbmDigiManager* fDigiMan       = nullptr;
  CbmMCEventList* fEventList     = nullptr;
 
  int fEventNum   = 0;
  int fMinNofHits = 7;  // Min number of hits in ring for detector acceptance calculation.
 
  std::map<CbmLink, int> fNofHitsInRingMap;  // Number of hits in the MC RICH ring
 
  std::vector<std::pair<int, int>> fVertexZStsSlices;  // Slices (Zmin - Zmax) inside STS
 
  void InitHistograms();
 
  void FillRichRingNofHits();
 
  void NofRings();
 
  void NofHitsAndPoints();
 
  void NofProjections();
 
  void Vertex();
 
 
  void DrawHist();
 
  void PmtXYSource();
 
  CbmRichUrqmdTest(const CbmRichUrqmdTest&);
 
  CbmRichUrqmdTest& operator=(const CbmRichUrqmdTest&);
 
  ClassDef(CbmRichUrqmdTest, 1)
};
 
#endif