---

AlpgenProducer.h

Go to the documentation of this file.

#ifndef Alpgen_Producer_h
#define Alpgen_Producer_h

#define PYCOMP pycomp_

#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include <map>
#include <vector>
#include <string>
#include <iostream>
#include <fstream>
#include "HepMC/GenEvent.h"

#include "SimDataFormats/GeneratorProducts/interface/LHERunInfoProduct.h"

namespace CLHEP
{
  class RandFlat ;
  class HepRandomEngine;
}

namespace edm
{
  class AlpgenProducer : public EDProducer {
  public:

    AlpgenProducer(const ParameterSet &);
    virtual ~AlpgenProducer();


  private:

    bool call_pygive(const std::string& iParm );
    bool call_txgive(const std::string& iParm );
    int Nev_; // number of events in the input file

    // .unw file with infos for AlpgenInfoProduct
//    std::ifstream* unwfile;

  private:
    
    virtual void produce(Event & e, const EventSetup& es);
    void beginRun(Run& r);
    void clear();    
    HepMC::GenEvent  *evt;
    
    unsigned int pythiaPylistVerbosity_;
    bool pythiaHepMCVerbosity_;
    unsigned int maxEventsToPrint_;    
    
    std::vector<std::string> fileNames_;
    std::string fileName_;
    unsigned int eventsRead_;
    
    // for single particle generation in pythia
    int    particleID;
    bool   doubleParticle;
    double ptmin, ptmax;
    double etamin, etamax;
    double phimin, phimax;
    
    CLHEP::HepRandomEngine* fRandomEngine;
    CLHEP::RandFlat*        fRandomGenerator; 

    LHERunInfoProduct::Header lheAlpgenUnwParHeader;
  };
} 

// 
#define alpgen_end alpgen_end_
    extern "C" {
        void alpgen_end();
      }

#endif

---