HepMC2Interface.hh

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <cstddef>
 
#include "HEJ/PDG_codes.hh"
 
namespace HepMC {
  class GenCrossSection;
  class GenEvent;
}
 
namespace LHEF {
  class HEPRUP;
}
 
namespace HEJ {
  class Event;
 
  class HepMC2Interface{
  public:
    HepMC2Interface(LHEF::HEPRUP const & /*heprup*/);
    HepMC::GenEvent operator()(Event const & event,
                               int weight_index = -1);
    void set_central(HepMC::GenEvent & out_ev, Event const & event,
                     int weight_index = -1);
 
    void set_xs_scale(double scale);
 
  protected:
    HepMC::GenEvent init_event(Event const & event) const;
 
  private:
    std::array<ParticleID,2> beam_particle_{};
    std::array<double,2> beam_energy_{};
    std::size_t event_count_;
    double tot_weight_;
    double tot_weight2_;
    double xs_scale_ = 1.;
    HepMC::GenCrossSection cross_section() const;
  };
} // namespace HEJ