// -*- C++ -*-
#ifndef RIVET_MC_JetMatrix_HH
#define RIVET_MC_JetMatrix_HH

//#include "/afs/cern.ch/sw/lcg/external/MCGenerators/rivet/1.2.1/share/include/Rivet/Analysis.hh"
#include "Rivet/Analysis.hh"
//#include "/afs/cern.ch/sw/lcg/external/MCGenerators/rivet/1.2.1/share/include/Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/FinalState.hh"
#include <iostream>
namespace Rivet {


  class MC_JetMatrix : public Analysis {

  public:

    /// Default constructor.
    MC_JetMatrix(const string& name,
                   const size_t& njet, 
                   const string& jetpro_name);


    /// @name Analysis methods
    //@{
    virtual void init();
    virtual void analyze(const Event& event);
    virtual void finalize();
    //@}


  protected:
 
    /// The number of jets for which histograms are to be initialised
    size_t m_njet;
 
    /// The name of the jet projection to be used for this analysis
    /// (this projection has to be registered by the derived analysis!)
    const std::string m_jetpro_name;

    /// @name Histograms
    //@{

   std::ofstream JetFile;
    std::vector<AIDA::IHistogram1D *> _h_log10_d;
    std::vector<AIDA::IDataPointSet *> _h_log10_R;
    std::vector<AIDA::IHistogram1D *> _h_pT_jet;
    std::vector<AIDA::IHistogram1D *> _h_eta_jet;
    std::map<std::pair<size_t, size_t>, AIDA::IHistogram1D*> _h_deta_jets;
    std::map<std::pair<size_t, size_t>, AIDA::IHistogram1D*> _h_dR_jets;
    AIDA::IHistogram1D * _h_jet_multi_exclusive;
    AIDA::IHistogram1D * _h_jet_multi_inclusive;
    AIDA::IDataPointSet * _h_jet_multi_ratio;
    AIDA::IHistogram1D * _h_jet1_etaVBF;
    AIDA::IHistogram1D * _h_jet2_etaVBF;
    AIDA::IHistogram1D * _h_jet1_ptVBF;
    AIDA::IHistogram1D * _h_jet2_ptVBF;
    AIDA::IHistogram1D * _h_jet_mass;
    AIDA::IHistogram1D * _h_jet_massVBF;
    AIDA::IHistogram1D * _h_jet_DetaVBF;
    AIDA::IHistogram1D * _h_jet_DRVBF;

    AIDA::IHistogram1D * _h_jet1_etaVBFex;
    AIDA::IHistogram1D * _h_jet2_etaVBFex;
    AIDA::IHistogram1D * _h_jet1_ptVBFex;
    AIDA::IHistogram1D * _h_jet2_ptVBFex;
    AIDA::IHistogram1D * _h_jet_massex;
    AIDA::IHistogram1D * _h_jet_massVBFex;
    AIDA::IHistogram1D * _h_jet_DetaVBFex;
    AIDA::IHistogram1D * _h_jet_DRVBFex;
    //@}

  };

}

#endif