Herwig7Interface.h

Go to the documentation of this file.

#ifndef GeneratorInterface_Herwig7Interface_Herwig7Interface_h
#define GeneratorInterface_Herwig7Interface_Herwig7Interface_h

#include <memory>
#include <string>

#include <HepMC/GenEvent.h>
#include <HepMC/PdfInfo.h>
#include <HepMC/IO_BaseClass.h>

#include <ThePEG/Repository/EventGenerator.h>
#include <ThePEG/EventRecord/Event.h>
#include <ThePEG/Vectors/HepMCTraits.h>

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "GeneratorInterface/Herwig7Interface/interface/RandomEngineGlue.h"
#include "GeneratorInterface/Herwig7Interface/interface/HerwigUIProvider.h"

namespace ThePEG {

  template <>
  struct HepMCTraits<HepMC::GenEvent> : public HepMCTraitsBase<HepMC::GenEvent,
                                                               HepMC::GenParticle,
                                                               HepMC::GenParticle *,
                                                               HepMC::GenVertex,
                                                               HepMC::GenVertex *,
                                                               HepMC::Polarization,
                                                               HepMC::PdfInfo> {};

}  // namespace ThePEG

namespace CLHEP {
  class HepRandomEngine;
}

class Herwig7Interface {
public:
  Herwig7Interface(const edm::ParameterSet &params);
  ~Herwig7Interface() noexcept;

  void setPEGRandomEngine(CLHEP::HepRandomEngine *);

  ThePEG::EGPtr eg_;

protected:
  void initRepository(const edm::ParameterSet &params);
  bool initGenerator();
  void flushRandomNumberGenerator();

  static std::auto_ptr<HepMC::GenEvent> convert(const ThePEG::EventPtr &event);

  static double pthat(const ThePEG::EventPtr &event);

  std::auto_ptr<HepMC::IO_BaseClass> iobc_;

  // HerwigUi contains settings piped to Herwig7
  std::shared_ptr<Herwig::HerwigUIProvider> HwUI_;

  void callHerwigGenerator();

  // The Inputfile ist created according to the parameter set
  void createInputFile(const edm::ParameterSet &params);

private:
  std::shared_ptr<ThePEG::RandomEngineGlue::Proxy> randomEngineGlueProxy_;

  const std::string dataLocation_;
  const std::string generator_;
  const std::string run_;
  // File name containing Herwig input config
  std::string dumpConfig_;
  const unsigned int skipEvents_;
  CLHEP::HepRandomEngine *randomEngine;
};

#endif  // GeneratorInterface_Herwig7Interface_Herwig7Interface_h