#ifndef __GFDSTARHISTSANALYSISDATA_H
#define __GFDSTARHISTSANALYSISDATA_H


#include "GFDstarHistsAnalysis.h"

class DstarDmFitter;
class GFHistArray;

class GFDstarHistsAnalysisData : public GFDstarHistsAnalysis {
public:
  GFDstarHistsAnalysisData(const char* file);
  GFDstarHistsAnalysisData(const char* prefix, const char* period);
  virtual ~GFDstarHistsAnalysisData();

  virtual void SetTrigger(Int_t trigger);
  //  void SetFitMode(Int_t mode) {fFitMode = mode;}
//   void SetFitOrNot(Short_t mode) {fFitOrNotFlag = (mode >=0 && mode <= 3) ? mode : fFitOrNotFlag;}// 0: fit, 1: # signal WITH or, 2: WITHOUT subtr. backgr. estimation, 3: count all delta m

  void DrawDmWc(Option_t* drawOpt = "");
  TArrayD DrawDmFitWc();

  TH1*  DrawDmFitWc(const char* var); //pt,eta,phi,wGammaP

  TH1*  DrawAccept(const char* var);// Pt,Phi,Eta,WgammaP
  TH1* CreateDedxEff(const char* var);//pt,eta,phi,ptSumEt,wGammaP,y,nTrack
  TH1* CreateDedxEff2(const char* var);//Theta,Length,Nhit,PK,P + K,Pi
  TArrayD TotalAccept();
  TArrayD L4Efficiency(Int_t mode = 0);// only 0 so far: inclusive D*

  void DrawTrigRefVsDm(const char* varName);
  TH1*  CreateTrigEff(const char* var, Int_t TE = -1);//pt eta phi wGammaP y; 19 31 both
  TArrayD TotalTrigEff(Int_t TE = -1);// 19 31 (both is default)

  TH1*   GetHistPt1CutWc();
  TH1*   GetHistPt2CutWc();
  TH1*   GetHistPtCutWc();
  TH2*   GetHist2DWc(const char* variable);
  TH1*   GetHistAccept(const char* variable);
  TH1*   GetHistAcceptNom(const char* variable);

  void DrawYresolution(const char* var); // Pt, Eta, Phi, WgammaP, Yjb, Ytag

protected:
  void CheckTrigEff(TH1* hEff, Int_t TE);

  //  virtual TArrayD FitDstar(TH1* dmHist, Int_t mode = 0) const;


  TArrayD  fL4Eff;    // store eff(L4) (inclusive) to avoid recalculations

//   Short_t  fFitOrNotFlag;// 0: fit, 1: # signal WITH or, 2: WITHOUT subtr. backgr. estimation, 3: count all delta m

  ClassDef(GFDstarHistsAnalysisData,0) // manipulation of histograms (real data)
};
#endif // __GFDSTARHISTsANALYSISDATA_H