ThePEGInterface.h

Go to the documentation of this file.

#ifndef GeneratorInterface_ThePEGInterface_ThePEGInterface_h
#define GeneratorInterface_ThePEGInterface_ThePEGInterface_h

#include <memory>
#include <string>

#include <boost/shared_ptr.hpp>

#include <HepMC/GenEvent.h>
#include <HepMC/PdfInfo.h>
#include <HepMC/IO_BaseClass.h>

#include <ThePEG/Repository/EventGenerator.h>
#include <ThePEG/EventRecord/Event.h>

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

#include "GeneratorInterface/ThePEGInterface/interface/RandomEngineGlue.h"
#include "GeneratorInterface/ThePEGInterface/interface/HepMCTemplate.h"

namespace CLHEP {
  class HepRandomEngine;
}

class ThePEGInterface {
    public:
    ThePEGInterface(const edm::ParameterSet &params);
    virtual ~ThePEGInterface();

        void setPEGRandomEngine(CLHEP::HepRandomEngine*);

    protected:
    void initRepository(const edm::ParameterSet &params) const;
    void initGenerator();
    void flushRandomNumberGenerator();

    static std::auto_ptr<HepMC::GenEvent>
                convert(const ThePEG::EventPtr &event);
    static void clearAuxiliary(HepMC::GenEvent *hepmc,
                               HepMC::PdfInfo *pdf);
    static void fillAuxiliary(HepMC::GenEvent *hepmc,
                              HepMC::PdfInfo *pdf,
                              const ThePEG::EventPtr &event);

    static double pthat(const ThePEG::EventPtr &event);

    std::string dataFile(const std::string &fileName) const;
    std::string dataFile(const edm::ParameterSet &pset,
                         const std::string &paramName) const;

    ThePEG::EGPtr               eg_;
    std::auto_ptr<HepMC::IO_BaseClass>  iobc_;

    private:
    boost::shared_ptr<ThePEG::RandomEngineGlue::Proxy>
                        randomEngineGlueProxy_;

    const std::string           dataLocation_;
    const std::string           generator_;
    const std::string           run_;
    const std::string           dumpConfig_;
    const unsigned int          skipEvents_;
};







#endif // GeneratorInterface_ThePEGInterface_ThePEGInterface_h