Pythia6Service.h

Go to the documentation of this file.

#ifndef gen_Pythia6Service_h
#define gen_Pythia6Service_h

#include <vector>

#include "FWCore/ParameterSet/interface/ParameterSetfwd.h"
#include "FWCore/Utilities/interface/Exception.h"

#include "GeneratorInterface/Core/interface/FortranInstance.h"

// #include "HepMC/PythiaWrapper6_2.h"

namespace CLHEP {
  class HepRandomEngine;
}

namespace gen {

  // the callbacks from Pythia which are passed on to the Pythia6Service
  extern "C" {
  double pyr_(int* idummy);
  }

  class Pythia6Service : public FortranInstance {
  public:
    // ctor & dtor
    Pythia6Service();
    Pythia6Service(edm::ParameterSet const&);
    ~Pythia6Service() override;

    void setGeneralParams();
    void setCSAParams();
    void setSLHAParams();
    void setPYUPDAParams(bool afterPyinit);
    void setSLHAFromHeader(const std::vector<std::string>& lines);

    void openSLHA(const char*);
    void closeSLHA();
    void openPYUPDA(const char*, bool write_file);
    void closePYUPDA();

    // initialise Pythia on first call from "dummy" instance
    void enter() override;

    CLHEP::HepRandomEngine* randomEngine() const { return fRandomEngine; }
    void setRandomEngine(CLHEP::HepRandomEngine* v) { fRandomEngine = v; }

  private:
    friend double gen::pyr_(int*);

    bool fInitialising;

    CLHEP::HepRandomEngine* fRandomEngine;

    std::vector<std::string> fParamGeneral;
    std::vector<std::string> fParamCSA;
    std::vector<std::string> fParamSLHA;
    std::vector<std::string> fParamPYUPDA;
    int fUnitSLHA;
    int fUnitPYUPDA;

    static Pythia6Service* fPythia6Owner;
  };

}  // namespace gen

#endif