CbmRichRecoTbQa.h

Go to the documentation of this file.

 
#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;
 
#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<pair<Int_t, Int_t>> fRecRings;
 
  ClassDef(CbmRichRecoTbQa, 1)
};
 
#endif