MC_PARTICLES_BASE.hh

// -*- C++ -*-
#ifndef RIVET_MC_PARTICLES_BASE_HH
#define RIVET_MC_PARTICLES_BASE_HH
 
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
 
namespace Rivet {
 

  class MC_PARTICLES_BASE : public Analysis {
  public:

    MC_PARTICLES_BASE(const string& name,
                      size_t nparticles,
                      const string& particle_name);
 


    virtual void init();

    virtual void analyze(const Event& event) = 0;

    virtual void finalize();

    virtual void _analyze(const Event& event, const Particles& particles);
 
 
  protected:

    size_t _nparts;

    std::string _pname;

    std::vector<Histo1DPtr> _h_pt;
    std::vector<Histo1DPtr> _h_eta;
    std::vector<Histo1DPtr> _h_eta_plus, _h_eta_minus;
    std::vector<Histo1DPtr> _h_rap;
    std::vector<Histo1DPtr> _h_rap_plus, _h_rap_minus;
    std::vector<Estimate1DPtr> tmpeta, tmprap;
 
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_deta;
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_dphi;
    std::map<std::pair<size_t, size_t>, Histo1DPtr> _h_dR;
    BinnedHistoPtr<int> _h_multi_exclusive, _h_multi_inclusive;
    BinnedHistoPtr<int> _h_multi_exclusive_prompt, _h_multi_inclusive_prompt;
    BinnedEstimatePtr<std::string> _h_multi_ratio, _h_multi_ratio_prompt;
 
  };
 
}
 
#endif