---

PythiaProducer.h

Go to the documentation of this file.

#ifndef Pythia_Producer_h
#define Pythia_Producer_h

#define PYCOMP pycomp_

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

#include "GeneratorInterface/CommonInterface/interface/TauolaInterface.h"
#include "GeneratorInterface/Pythia6Interface/interface/PtYDistributor.h"

class Run;
namespace CLHEP {
class HepRandomEngine;
class RandFlat;
}

namespace edm
{
  class PythiaProducer : public EDProducer {
  public:

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

    void endRun( Run& r, const EventSetup& es);

  private:

    bool call_pygive(const std::string& iParm );
    bool call_txgive(const std::string& iParm );
    bool call_txgive_init();
    bool call_slhagive(const std::string& iParm );
    bool call_slha_init();
  
  private:
    
    virtual void produce(Event & e, const EventSetup& es);
    void clear();
    
    HepMC::GenEvent  *evt;

    // for pyinit
    std::string fFrame ;
    std::string fBeam1;
    std::string fBeam2;
    double      fCOMEnergy;
        
    unsigned int pythiaPylistVerbosity_;
    bool pythiaHepMCVerbosity_;
    bool imposeProperTimes_;
    unsigned int maxEventsToPrint_;    
   
    // external cross section and filter efficiency
    double extCrossSect;
    double extFilterEff;    
 
    // for single particle generation in pythia
    int    particleID;
    bool   doubleParticle;
    std::string kinedata;
    double ptmin, ptmax;
    double etamin, etamax;
    double phimin, phimax;
    double emin, emax;
    double ymin, ymax;
    bool flatEnergy;

    // external generators (tauola,...)
    bool useExternalGenerators_ ;
    bool useTauola_ ;
    bool useTauolaPolarization_ ;
    TauolaInterface tauola_ ;
    
    bool stopHadronsEnabled;
    bool gluinoHadronsEnabled;
    
    CLHEP::HepRandomEngine& fRandomEngine;
    CLHEP::RandFlat*        fRandomGenerator; 
    PtYDistributor*         fPtYGenerator;
    int                     eventNumber_;
  };
}

#endif

---