Pythia6Gun.h

Go to the documentation of this file.

#ifndef gen_Pythia6Gun_h
#define gen_Pythia6Gun_h

#include <string>
#include <vector>

//#include "HepPDT/defs.h"
//#include "HepPDT/TableBuilder.hh"
//#include "HepPDT/ParticleDataTable.hh"

#include "HepMC/GenEvent.h"

// #include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "GeneratorInterface/Pythia6Interface/interface/Pythia6Service.h"
#include "GeneratorInterface/Pythia6Interface/interface/Pythia6Declarations.h"

#include "HepPID/ParticleIDTranslations.hh"

namespace CLHEP {
  class HepRandomEngine;
}

namespace edm {
  class EventSetup;
  class LuminosityBlock;
}  // namespace edm

namespace gen {

  // class Pythia6Service;

  class Pythia6Gun : public edm::EDProducer {
  public:
    Pythia6Gun(const edm::ParameterSet&);
    ~Pythia6Gun() override;
    void beginJob() override;
    void beginLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) override;
    void beginRun(edm::Run const&, edm::EventSetup const&) override;
    void endRun(edm::Run const&, edm::EventSetup const&) override;
    void produce(edm::Event&, const edm::EventSetup&) override;

  protected:
    void attachPy6DecaysToGenEvent();
    void loadEvent(edm::Event&);
    virtual void generateEvent(CLHEP::HepRandomEngine*) = 0;
    HepMC::GenParticle* addAntiParticle(int&, int&, double&, double&, double&);

    Pythia6Service* fPy6Service;

    // gun particle(s) characteristics
    //
    std::vector<int> fPartIDs;
    double fMinPhi;
    double fMaxPhi;

    // the event format itself
    //
    HepMC::GenEvent* fEvt;

    // HepMC/HepPDT related things
    // (for particle/event construction)
    //
    // edm::ESHandle<HepPDT::ParticleDataTable> fPDGTable ;

    bool fHepMCVerbosity;
    int fPylistVerbosity;
    int fMaxEventsToPrint;
  };

}  // namespace gen

#endif