Herwig7.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
// -*- C++ -*-

 
#ifndef GENERATORMODULES_HERWIG7_H
#define GENERATORMODULES_HERWIG7_H
 
#include "GeneratorModules/GenModule.h"
 
#include "ThePEG/Repository/Repository.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "ThePEG/EventRecord/Event.h"
#include "Herwig/API/HerwigUI.h"
 
#include "xAODEventInfo/EventInfo.h"

class Herwig7API : public Herwig::HerwigUI {
 
public:

  Herwig7API(){};
  //   ATH_MSG_DEBUG("HERWIG7_PATH is "+getenv("HERWIG7_PATH"));
  //   ThePEG::Repository::load(getenv("HERWIG7_PATH")+"/share/Herwig/HerwigDefaults.rpo");
  // };

  ~Herwig7API(){};
 

  Herwig::RunMode::Mode runMode() const { return(m_runMode); }

  virtual bool resume() const { return false; }

  virtual bool tics() const { return false; }
 

  virtual std::string repository() const { return(m_repository); }

  virtual std::string inputfile() const { return(m_inputfile); }
  void inputfile(const std::string& inputfile) { m_inputfile = inputfile; }

  virtual std::string setupfile() const { return(m_setupfile); }

  virtual std::string tag() const { return(m_tag); }

  virtual std::string integrationList() const { return(m_integrationList); }
 

  virtual const std::vector<std::string> & prependReadDirectories() const { return(m_prependReadDirectories); }
  void prependReadDirectory(const std::string& prependReadDirectory) { m_prependReadDirectories.push_back(prependReadDirectory); }

  virtual const std::vector<std::string> & appendReadDirectories() const { return(m_appendReadDirectories); }
 
 
  virtual long N() const { return 0; }      
 
  virtual int seed() const { return(m_seed); } 
  void seed(int seed) { m_seed = seed; }
 
  virtual int jobs() const { return(m_jobs); } 
  virtual unsigned int jobSize() const { return(m_jobsize); } 
  virtual unsigned int maxJobs() const { return(m_maxjobs); } 

  virtual void quitWithHelp() const {
    std::cerr << "=================================\n";
    std::cerr << "HERE I CAN PUT SOME ERROR MESSAGE\n";
    std::cerr << "=================================\n";
    quit();
  }

  virtual void quit() const {
    ThePEG::Repository::cleanup();
    std::abort();
  }

  virtual std::ostream& outStream() const { return(std::cout); }

  virtual std::ostream& errStream() const { return(std::cerr); }

  virtual std::istream& inStream() const { return(std::cin); }
 
 
private:
 
  Herwig::RunMode::Mode m_runMode{};
 
  std::string m_repository;
  std::string m_inputfile;
  std::string m_setupfile;
  std::string m_tag;
 
  std::string m_integrationList;
 
  std::vector<std::string> m_prependReadDirectories;
  std::vector<std::string> m_appendReadDirectories;
 
  int m_seed{};
  int m_jobs{};
  unsigned int m_jobsize{};
  unsigned int m_maxjobs{};
 
};
 
 

class Herwig7 : public GenModule {
 
public:

  Herwig7(const std::string& name, ISvcLocator* pSvcLocator);

  StatusCode genInitialize();

  StatusCode callGenerator();

  StatusCode fillEvt(HepMC::GenEvent* evt);

  StatusCode genFinalize();

  StatusCode writeRunFileFromText(const std::string& share_path);
 
 
private:

  Herwig7API m_api;

  ThePEG::EGPtr m_gen;

  ThePEG::EventPtr m_event;
 



  std::string m_runname;

  std::string m_runfile;

  std::string m_runSettings;

  std::string m_setupfile;
 

  bool m_use_seed_from_generatetf;

  int m_seed_from_generatetf;

  std::string m_pdfname_me;

  // std::string m_pdfname_ps;

  std::string m_pdfname_mpi;

  bool m_cleanup_herwig_scratch;

  double m_xsscale;

  int m_dsid;
 
  SG::ReadHandleKey<xAOD::EventInfo>  m_evtInfoKey{this
      , "EventInfo"
      ,  "TMPEvtInfo"
      , "ReadHandleKey for xAOD::EventInfo"};

 
};
 
#endif