CbmAnaDielectronTaskDrawAll.h

Go to the documentation of this file.

 
#ifndef CBM_ANA_DIELECTRON_TASK_DRAW_ALL
#define CBM_ANA_DIELECTRON_TASK_DRAW_ALL
 
#include <map>
#include <string>
#include <vector>
 
#include "CbmLmvmHist.h"
 
#include "TObject.h"
 
class TH1;
class TH2D;
class TH1D;
class TFile;
class TCanvas;
class CbmHistManager;
 
enum SignalType { kInmed = 0, kQgp = 1, kOmega = 2, kPhi = 3, kOmegaD = 4 };
 
class CbmAnaDielectronTaskDrawAll : public TObject {
 
public:
  CbmAnaDielectronTaskDrawAll()
    : TObject()
    , fUseMvd(false)
    , fDrawQgp(false)
    , fHM()
    , fh_mean_bg_minv()
    , fh_mean_eta_minv()
    , fh_mean_pi0_minv()
    , fh_mean_eta_minv_pt()
    , fh_mean_pi0_minv_pt()
    , fh_mean_sbg_vs_minv()
    , fh_sum_s_minv()
    , fOutputDir("") {
    ;
  }
 
  virtual ~CbmAnaDielectronTaskDrawAll() { ; }
 
  void DrawHistosFromFile(const std::string& fileNameInmed,
                          const std::string& fileNameQgp,
                          const std::string& fileNameOmega,
                          const std::string& fileNamePhi,
                          const std::string& fileNameOmegaDalitz,
                          const std::string& outputDir = "",
                          Bool_t useMvd                = false);
 
private:
  static const int fNofSignals = 5;
 
  Bool_t
    fUseMvd;  // do you want to draw histograms related to the MVD detector?
  Bool_t fDrawQgp;  // do you wan to draq QGP signal
 
  //[0]=rho0, [1]=omega, [2]=phi, [3]=omegaDalitz
  std::vector<CbmHistManager*> fHM;
 
  // index: AnalysisSteps
  std::vector<TH1D*> fh_mean_bg_minv;  //mean histograms from all files
  std::vector<TH1D*> fh_mean_eta_minv;
  std::vector<TH1D*> fh_mean_pi0_minv;
  std::vector<TH2D*> fh_mean_eta_minv_pt;
  std::vector<TH2D*> fh_mean_pi0_minv_pt;
  std::vector<TH1D*>
    fh_mean_sbg_vs_minv;  //Coctail/BG vs. invariant mass for different analysis steps
 
 
  // index: AnalysisSteps
  std::vector<TH1D*> fh_sum_s_minv;  // sum of all signals
 
  std::string fOutputDir;  // output directory for figures
 
  TH1D* H1(int signalType, const std::string& name);
 
  TH2D* H2(int signalType, const std::string& name);
 
  TH1D* GetCoctailMinv(CbmLmvmAnalysisSteps step);
 
  void DrawSBgVsMinv();
 
  void DrawMinvAll();
 
  void DrawMinv(CbmLmvmAnalysisSteps step);
 
  void DrawMinvPtAll();
 
  void DrawMinvPt(CbmLmvmAnalysisSteps step);
 
  void FillMeanHist();
 
  void SaveHist();
 
  void FillSumSignalsHist();
 
  void CalcCutEffRange(Double_t minMinv, Double_t maxMinv);
 
  TH1D* SBgRange(Double_t minMinv, Double_t maxMinv);
 
  void SBgRangeAll();
 
  void DrawSBgSignals();
 
  void SaveCanvasToImage();
 
 
  CbmAnaDielectronTaskDrawAll(const CbmAnaDielectronTaskDrawAll&);
  CbmAnaDielectronTaskDrawAll operator=(const CbmAnaDielectronTaskDrawAll&);
 
  ClassDef(CbmAnaDielectronTaskDrawAll, 1);
};
 
#endif