// created    : May 23rd, 2003
// by         : Gero Flucke
// last update: $Date: 2005/12/29 16:02:26 $
// by         : $Author: flucke $

#ifndef __GFDSTARCOMPARE_H
#define __GFDSTARCOMPARE_H

#include <vector>

#include "TObject.h"

class TH1;

class GFDstarHistsAnalysisData;
class GFDstarHistsAnalysisMc;

class GFHistManager;

class GFDstarCompare : public TObject {
public:
  GFDstarCompare();
  GFDstarCompare(const char* filePrefix);
  GFDstarCompare(const char* data, const char* mc);
  virtual ~GFDstarCompare();
  void AdoptData(GFDstarHistsAnalysisData *data);
  void AdoptMc(GFDstarHistsAnalysisMc *m);
  TH1* L4Eff(const char *var, Int_t dstarL4 = 1);// L4 D*-finder, 0: L4-class
  void RecEff(const char* var);
  void DetAcceptance(const char* var);
  void HadCorr(const char *var, const char *forwBack="");
  void Accept(const char* var);
  void DstarTracks();
  void TrigEff(Int_t TE = -1);//19, 31 (default: -1, i.e. both)
  void JetProfile(Int_t flagBins = 1, Int_t etaPhi = 0, Bool_t add2nd = kFALSE);//0: overall, 1: overall + pt/eta bins, 2: pt/eta_bins
  void DrawDmShape(const char* var);//
  void JetComposition(const char* var = NULL, Int_t flagBins=1);//Track,EmCl,HadCl; 0: overall, 1: overall + pt bins, 2: pt_bins
  void JetCompSRMean();
//   void JetPurityStability(const char* var = NULL);// n, eta, pt, phi, default: all
  void DsJetPurityStability(const char* var = NULL, Bool_t hadCorr = kFALSE);//Dphi Deta xGam PtJet PtDs,def: all

//   TH1* JetResEta(Int_t etaBin = 0); // relative pt-resolution
  void Ds1Jet(const char* forwBack);// "Forw", "Back", ""
  void TrackErr(Option_t * opt = "rel"); 
  void SetTrigger(Int_t subTr);
  GFHistManager* GetHistManager() {return fHistManager;}
  GFDstarHistsAnalysisData* GetHistsAnalysisData(Int_t i);
  GFDstarHistsAnalysisMc* GetHistsAnalysisMc(Int_t i);
protected:
  GFHistManager*                  fHistManager;

  //  GFDstarHistsAnalysisData*       fData;
  std::vector<GFDstarHistsAnalysisData*> fDatas;
  std::vector<GFDstarHistsAnalysisMc*> fMonteCarlo;

  

  ClassDef(GFDstarCompare,0) // compare several MC / Data
};
#endif