CbmRichRecoTbQa.h

Go to the documentation of this file.

/* Copyright (C) 2018-2020 GSI/JINR-LIT, Darmstadt/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Andrey Lebedev [committer] */
 
#ifndef CBM_RICH_RECO_TB_QA
#define CBM_RICH_RECO_TB_QA
 
#include "FairTask.h"
class CbmMCDataArray;
class TClonesArray;
class CbmHistManager;
class CbmRichPoint;
class CbmMCTrack;
class CbmMCEventList;
class CbmDigiManager;
class CbmLink;
 
#include <map>
#include <vector>
 
using namespace std;
 
class CbmRichRecoTbQa : public FairTask {
 
 public:
  CbmRichRecoTbQa();
 
  virtual ~CbmRichRecoTbQa() {}
 
  virtual InitStatus Init();
 
  virtual void Exec(Option_t* option);
 
  virtual void Finish();
 
 
  void SetOutputDir(const string& dir) { fOutputDir = dir; }
 
  void DrawFromFile(const string& fileName, const string& outputDir);
 
 private:
  void InitHistograms();
 
  void DrawHist();
 
  void Process();
 
  void ProcessMc();
 
  void RingRecoEfficiency();
 
  Int_t GetNofPrimaryMcTracks(Int_t iEv);
 
  Bool_t IsCherenkovPhoton(const CbmRichPoint* point, Int_t fileId, Int_t eventId);
 
  Bool_t IsCherenkovPhotonFromPrimaryElectron(const CbmRichPoint* point, Int_t fileId, Int_t eventId);
 
  Bool_t IsCherenkovPhotonFromSecondaryElectron(const CbmRichPoint* point, Int_t fileId, Int_t eventId);
 
  Bool_t IsMcPrimaryElectron(const CbmMCTrack* mctrack);
 
  Bool_t IsCherenkovPhotonFromPion(const CbmRichPoint* point, Int_t fileId, Int_t eventId);
 
  Bool_t IsMcPion(const CbmMCTrack* mctrack);
 
  void DrawTimeLog(const string& hMainName, Int_t nofLogEvents, bool withNoise = false);
 
  CbmRichRecoTbQa(const CbmRichRecoTbQa&);
 
  CbmRichRecoTbQa& operator=(const CbmRichRecoTbQa&);
 
 
  CbmHistManager* fHM;
 
  Int_t fTimeSliceNum;
 
  Int_t fNofLogEvents;
  string fOutputDir;  // output dir for results
 
  CbmMCDataArray* fMCTracks;
  CbmMCDataArray* fRichPoints;
  CbmMCDataArray* fStsPoints;
  CbmDigiManager* fDigiMan;
  TClonesArray* fRichHits;
  TClonesArray* fRichRings;
  TClonesArray* fRichRingMatches;
  CbmMCEventList* fEventList;
 
  vector<CbmLink> fRecRings;
 
  ClassDef(CbmRichRecoTbQa, 1)
};
 
#endif