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