MC_JETS_BASE.hh

// -*- C++ -*-
#ifndef RIVET_MC_JETS_BASE_HH
#define RIVET_MC_JETS_BASE_HH
 
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
 
namespace Rivet {
 
 
  class MC_JETS_BASE : public Analysis {
  public:
 
    MC_JETS_BASE(const string& name,
                 size_t njet,
                 const string& jetpro_name,
                 double jetptcut=20*GeV);
 
 
    virtual void init();
    virtual void analyze(const Event& event);
    virtual void finalize();
 
 
  protected:
 
    size_t _njet;
 
    const std::string _jetpro_name;
 
    double _jetptcut;
 
 
    std::vector<Histo1DPtr> _h_pT_jet;
    std::vector<Histo1DPtr> _h_eta_jet;
    std::vector<Histo1DPtr> _h_eta_jet_plus, _h_eta_jet_minus;
    std::vector<Histo1DPtr> _h_rap_jet;
    std::vector<Histo1DPtr> _h_rap_jet_plus, _h_rap_jet_minus;
    std::vector<Histo1DPtr> _h_mass_jet;
    std::vector<Estimate1DPtr> tmpeta, tmprap;
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_deta_jets;
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_dphi_jets;
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_dR_jets;
    BinnedHistoPtr<int> _h_jet_multi_exclusive;
    BinnedHistoPtr<int> _h_jet_multi_inclusive;
    BinnedEstimatePtr<std::string> _h_jet_multi_ratio;
    Histo1DPtr _h_jet_HT, _h_mjj_jets;
 
  };
 
}
 
#endif