MCatNLOSource.h

Go to the documentation of this file.

#ifndef GeneratorInterface_MCatNLOInterface_MCatNLOSource_h
#define GeneratorInterface_MCatNLOInterface_MCatNLOSource_h

#include <memory>

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Sources/interface/ProducerSourceFromFiles.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "SimDataFormats/GeneratorProducts/interface/LesHouches.h"

#include "GeneratorInterface/Herwig6Interface/interface/Herwig6Instance.h"

#include <fstream>

// Common Block for HERWIG parameters set by UPINIT
extern "C" {
extern struct MCPARS_ {
  double emmin;
  double emmax;
  double gammax;
  double rmass[1000];
  double gamw;
  double gamz;
  int emmins;
  int emmaxs;
  int gammaxs;
  int rmasss[1000];
  int gamws;
  int gamzs;
} mcpars_;
}

namespace lhef {
  class LHERunInfo;
  class LHEEvent;
}  // namespace lhef

class MCatNLOSource : public edm::ProducerSourceFromFiles, public gen::Herwig6Instance {
public:
  explicit MCatNLOSource(const edm::ParameterSet &params, const edm::InputSourceDescription &desc);
  ~MCatNLOSource() override;

private:
  void endJob() override;
  void beginRun(edm::Run &run) override;
  bool setRunAndEventInfo(edm::EventID &, edm::TimeValue_t &, edm::EventAuxiliary::ExperimentType &) override;
  void produce(edm::Event &event) override;

  void nextEvent();

  bool hwwarn(const std::string &fn, int code) override;

  std::string fileName;

  std::unique_ptr<std::ifstream> inputFile;

  unsigned int skipEvents;

  unsigned int nEvents;

  int ihpro;

  int processCode;

  std::unique_ptr<std::ifstream> reader;

  std::shared_ptr<lhef::LHERunInfo> runInfo;
  std::shared_ptr<lhef::LHEEvent> event;
};

#endif  // GeneratorInterface_MCatNLOInterface_MCatNLOSource_h