#ifndef __GFDSTARPYTHIA_H
#define __GFDSTARPYTHIA_H

//////////////////////////////////////////////////////////////////////////
//                                                                      //
// GFDstarPythia                                                      //
//                                                                      //
//                                                                      //
//   Author:      Gero Flucke                                           //
//   Date:        Jan 6th, 2005                                         //
//   last update: $Date: 2005/11/07 18:47:03 $
//   by:          $Author: flucke $
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include <TArrayD.h>
#include <TObject.h>
#include <TString.h>

#include "GFDstarVirtualTheory.h"


class TH1;
class GFHistArray;
class GFDstarHistsAnalysisMc;

// well PYTHIA is not NLO, but...
class GFDstarPythia : public GFDstarVirtualTheory {
public:
  explicit GFDstarPythia(const char *histMcLable);
  virtual ~GFDstarPythia();

  //  virtual void DrawHists(const char *var, Int_t flag, Bool_t rel = kFALSE);
  virtual TArrayD  TotalCrossSec(ECrossSecType dstarJet, Int_t dirResFlag = 0);// < 0 res, > 0 dir

  virtual GFHistArray* CreateHistsFinal(const char *var, Int_t dirResFlag=0);//flag obso.
  virtual GFHistArray* CreateHistsFinalNorm(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 = "");
  virtual TH1* CreateHistResMassive(const char *var, const char *note = "");
  virtual TH1* CreateHistExcitation(const char *var, const char *note = "");
  virtual const char* GetName() const;
  virtual Bool_t ShouldApplyHadCor() const;

  GFDstarHistsAnalysisMc* GetHistsAnalysisMc() {return fMc;}
  Bool_t DecodeDoubleDiff(const char *var, TString& var1, TString& var2,
			  Int_t& bin1, Int_t& bin2) const;
private:
  TH1* CreateHistFrom(const char *var, const char *note, GFDstarHistsAnalysisMc *mc) const;
  GFDstarHistsAnalysisMc*  fMc;
  GFDstarHistsAnalysisMc*  fMcDir;
  GFDstarHistsAnalysisMc*  fMcRes;
  GFDstarHistsAnalysisMc*  fMcExc;

  ClassDef(GFDstarPythia,0)  //interface to pythia hists
};

#endif