#ifndef __GFDSTARFMNR_H
#define __GFDSTARFMNR_H

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// GFDstarFmnr                                                          //
//                                                                      //
//                                                                      //
//   Author:      Gero Flucke                                           //
//   Date:        Nov 3rd, 2004                                         //
//   last update: $Date: 2005/11/07 18:47:03 $
//   by:          $Author: flucke $
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include "TObject.h"
#include "TString.h"
#include "TArrayD.h"

#include "GFDstarVirtualTheory.h"

class TH1;
class GFHistArray;
class GFHistManager;

class GFDstarFmnr : public GFDstarVirtualTheory {
public:
//   GFDstarFmnr(const char *fileDirIncl = "q", const char *fileDirDsjet = "q2d",
  GFDstarFmnr(const char *fileDirIncl = "", const char *fileDirDsjet = "GFcyeq6mSlice032",
	      const char *mainDirectory = "/user/flucke/dstar/fmnrresults");
  virtual ~GFDstarFmnr();

  virtual void DrawHists(const char *var, Int_t flag, Bool_t rel = kFALSE);//1D,2D,PtDs;0=final, 1=m_c, 2=mu_f, 3=m_r, 4=e_pet
  void DrawHistsOrigBin(const char *var, const char *note = "");

  virtual GFHistArray* CreateHistsFinal(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsFinalOld(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  virtual GFHistArray* CreateHistsFinalNorm(const char *var, Int_t dirResFlag = 0); // < 0 res, > 0 dir
  virtual TArrayD  TotalCrossSec(ECrossSecType dstarJet, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsScaleR(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsScaleF(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsScaleOld(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsMass(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir
  GFHistArray* CreateHistsEpet(const char *var, Int_t dirResFlag = 0);// < 0 res, > 0 dir

  virtual TH1* CreateHist(const char *var, const char *note = "");
  virtual TH1* CreateHistDir(const char *var, const char *note = "");
  virtual TH1* CreateHistRes(const char *var, const char *note = "");

  TH1* CreateHistOrigBin(const char *var, const char *note = "");
  TH1* CreateHistOrigBinDir(const char *var, const char *note = "");
  TH1* CreateHistOrigBinRes(const char *var, const char *note = "");
  
//   TH1* CreateHist1D(TH2 *h2D, const char *var);

  virtual const char* GetName() const;
  void SetNameAdd(const char *add);
  void SetFileStartDirNlo(const char *fileStart) {fFileStartDirNlo = fileStart;}
  void SetFileStartResNlo(const char *fileStart) {fFileStartResNlo = fileStart;}

  Bool_t SetUse2D(Bool_t flag) 
  {const Bool_t old = fUse2D; fUse2D = flag; return old;}

private:
  const char* GetHistName(const char *var, const char *&fileDir) const;
  TH1* GetHistResW(const char *note);
  Bool_t Scale(TH1 *hist, const char *name) const;

  Bool_t       fUse2D;// if true always take plots from double diffmode

  TString      fName;
  TString      fFileDirIncl;
  TString      fFileDirDsjet;
  TString      fFileStartDirNlo;
  TString      fFileStartResNlo;

  Double_t     fCharmToDstar; // c->D* branching
  Double_t     fNbFactor;     // factor to convert FMNR into nanobarn

  TArrayD      fFluxesW;      // integrated photon fluxes in W-bins
  TArrayD      fFluxesXgam;   // integrated photon fluxes in y-bins to add up xGamma

  ClassDef(GFDstarFmnr,0)  //interface to NLO hists
};

typedef GFDstarFmnr GetDstarNloFmnr;
#endif