RivetAnalysis.hh

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
#include <string>
#include <vector>
 
#include "HEJ/Analysis.hh"
#include "HEJ/HepMCInterface.hh"
 
namespace Rivet {
  class AnalysisHandler;
}
namespace YAML {
  class Node;
}
 
namespace HEJ {
  class RivetAnalysis: public HEJ::Analysis {
  public:
    static std::unique_ptr<Analysis> create(YAML::Node const & config);
 
 
    RivetAnalysis(YAML::Node const & config);
    void fill(HEJ::Event const & event, HEJ::Event const &) override;
    bool pass_cuts(HEJ::Event const &, HEJ::Event const &) override
       {return true;} //< no additional cuts are applied
    void finalise() override;
  private:
    std::vector<std::string> analyses_names_;
    const std::string output_name_;
 
    struct rivet_info {
      std::unique_ptr<Rivet::AnalysisHandler> handler;
      std::string name;
      HEJ::HepMCInterface hepmc;
    };
    std::vector<rivet_info> rivet_runs_;
 
    void init(HEJ::Event const & event);
    bool first_event_;
  };
}