JetMatching.cc

Go to the documentation of this file.

#include <memory>

#include <memory>
#include <string>

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

#include "GeneratorInterface/LHEInterface/interface/LHERunInfo.h"
#include "GeneratorInterface/LHEInterface/interface/LHEEvent.h"

#include "GeneratorInterface/PartonShowerVeto/interface/JetMatchingMadgraph.h"
#include "GeneratorInterface/PartonShowerVeto/interface/JetMatchingAlpgen.h"

namespace gen {

  JetMatching::JetMatching(const edm::ParameterSet& params) { fMatchingStatus = false; }

  JetMatching::~JetMatching() {}

  void JetMatching::init(const lhef::LHERunInfo* runInfo) {}

  void JetMatching::beforeHadronisation(const lhef::LHEEvent* event) {}

  void JetMatching::beforeHadronisationExec() {}

  std::set<std::string> JetMatching::capabilities() const {
    std::set<std::string> result;
    result.insert("psFinalState");
    result.insert("hepmc");
    return result;
  }

  std::unique_ptr<JetMatching> JetMatching::create(const edm::ParameterSet& params) {
    std::string scheme = params.getParameter<std::string>("scheme");

    std::unique_ptr<JetMatching> matching;

    if (scheme == "Madgraph") {
      matching = std::make_unique<JetMatchingMadgraph>(params);
    } else if (scheme == "Alpgen") {
      matching = std::make_unique<JetMatchingAlpgen>(params);
    } else if (scheme == "MLM") {
      matching.reset();
    } else
      throw cms::Exception("InvalidJetMatching") << "Unknown scheme \"" << scheme
                                                 << "\""
                                                    " specified for parton-shower matching."
                                                 << std::endl;

    if (!matching.get())
      throw cms::Exception("InvalidJetMatching") << "Port of " << scheme << "scheme \""
                                                 << "\""
                                                    " for parton-shower matching is still in progress."
                                                 << std::endl;

    return matching;
  }

}  // namespace gen