#ifndef __GFDSTARZMVFNS_H
#define __GFDSTARZMVFNS_H

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// GFDstarZmvfns                                                      //
//                                                                      //
//                                                                      //
//   Author:      Gero Flucke                                           //
//   Date:        Dec 20th, 2004                                         //
//   last update: $Date: 2005/11/07 18:47:03 $
//   by:          $Author: flucke $
//                                                                      //
//////////////////////////////////////////////////////////////////////////

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

#include "GFDstarVirtualTheory.h"


class TH1;
class GFHistArray;

class GFDstarZmvfns : public GFDstarVirtualTheory {
public:
GFDstarZmvfns(const char *fileAdd = "_mrst04_grv", 
	      const char *filePath = "/afs/desy.de/user/f/flucke/h1/nlo/kniehl/dsjet/mrst04/withdoublediff",
	      const char *fileStartDir = "jgdfdstar",
	      const char *fileStartRes = "jgrfdstar");
  virtual ~GFDstarZmvfns();

  void DrawHists(const char *var, Int_t flag, Bool_t rel=kFALSE);//1D,2D,PtDs;0,1=final=scale
  virtual TArrayD TotalCrossSec(ECrossSecType dstarJet, Int_t dirResFlag = 0);

  virtual GFHistArray* CreateHistsFinal(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
  GFHistArray* CreateHistsScale(const char *var, Int_t dirResFlag);
  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 = "");

  virtual const char* GetName() const;

  void SetFileAdd(const char *add);

private:
  const char* GetHistName(const char *var, const char *& fileExt) const;
  //  TH1* GetHist(const char *hName, const char *fileName);
  Bool_t Scale(TH1 *hist, const char *var) const;

  void SetUpperHistEgde(TH1 *&hist, Double_t upper, Bool_t sw2 = kFALSE) const;

  //  TString      fFilePath;
  TString      fFileStartDirNlo;
  TString      fFileStartResNlo;
  TString      fFileAdd;
  TString      fName;

  //  TObjArray   *fOpenFiles;
  Double_t     fNbFactor;     // factor to convert massless into nanobarn

  ClassDef(GFDstarZmvfns,0)  //interface to massless NLO hists
};
typedef class GFDstarZmvfns GFDstarNloMassless;
#endif