/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/GeneratorInterface/PartonShowerVeto/src/JetMatching.cc

Go to the documentation of this file.

#include <string>
#include <memory>

#include <boost/shared_ptr.hpp>

#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::auto_ptr<JetMatching> JetMatching::create(const edm::ParameterSet &params)
{
        std::string scheme = params.getParameter<std::string>("scheme");

        std::auto_ptr<JetMatching> matching;

        if (scheme == "Madgraph")
        {
                matching.reset(new JetMatchingMadgraph(params));
        }
        else if (scheme == "Alpgen")
        {
                matching.reset(new 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