HepMC3Interface.hh

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cstddef>
#include <memory>
 
#include "HEJ/PDG_codes.hh"
 
namespace HepMC3 {
  class GenCrossSection;
  class GenEvent;
  class GenRunInfo;
}
 
namespace LHEF {
  class HEPRUP;
}
 
namespace HEJ {
  class Event;
 
  class HepMC3Interface{
  public:
    HepMC3Interface(LHEF::HEPRUP heprup);
    ~HepMC3Interface();
    HepMC3Interface & operator=(HepMC3Interface const & other) = default;
    HepMC3Interface(HepMC3Interface const & other) = default;
    HepMC3Interface & operator=(HepMC3Interface && other) = default;
    HepMC3Interface(HepMC3Interface && other) = default;
    HepMC3::GenEvent operator()(Event const & event, int weight_index = -1);
    void set_central(HepMC3::GenEvent & out_ev, Event const & event,
      int weight_index = -1);
 
    std::shared_ptr<HepMC3::GenRunInfo> run_info(){
      return run_info_;
    }
 
    void set_xs_scale(double scale);
 
  protected:
    HepMC3::GenEvent init_event(Event const & event) const;
 
  private:
    std::array<ParticleID,2> beam_particle_{};
    std::array<double,2> beam_energy_{};
    std::shared_ptr<HepMC3::GenRunInfo> run_info_;
    std::size_t event_count_;
    double tot_weight_;
    double tot_weight2_;
    double xs_scale_ = 1.;
    std::shared_ptr<HepMC3::GenCrossSection> xs_;
  };
} // namespace HEJ