Pythia6Hadronizer.h

Go to the documentation of this file.

#ifndef gen_Pythia6Hadronizer_h
#define gen_Pythia6Hadronizer_h

// -*- C++ -*-

// class Pythia6Hadronizer is an example of a class that models the
// Hadronizer concept.

#include <memory>
#include <string>
#include <vector>

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

#include "GeneratorInterface/Core/interface/BaseHadronizer.h"

namespace lhef {
  class LHERunInfo;
  class LHEEvent;
}  // namespace lhef

class LHEEventProduct;

namespace HepMC {
  class GenEvent;
}

namespace CLHEP {
  class HepRandomEngine;
}

namespace gen {

  class Pythia6Service;
  class JetMatching;

  class Pythia6Hadronizer : public BaseHadronizer {
  public:
    Pythia6Hadronizer(edm::ParameterSet const& ps);
    ~Pythia6Hadronizer() override;

    // bool generatePartons();
    bool generatePartonsAndHadronize();
    bool hadronize();
    bool decay();
    bool residualDecay();
    bool readSettings(int);
    bool initializeForExternalPartons();
    bool initializeForInternalPartons();
    bool declareStableParticles(const std::vector<int>&);
    bool declareSpecialSettings(const std::vector<std::string>&);

    static JetMatching* getJetMatching() { return fJetMatching; }

    void finalizeEvent();

    void statistics();

    const char* classname() const;

  private:
    // methods
    //

    void doSetRandomEngine(CLHEP::HepRandomEngine* v) override;
    std::vector<std::string> const& doSharedResources() const override { return theSharedResources; }

    void flushTmpStorage();
    void fillTmpStorage();

    void imposeProperTime();

    // data members
    //

    Pythia6Service* fPy6Service;

    // some of the following params are common for all generators(interfaces)
    // probably better to wrap them up in a class and reuse ?
    // (the event/run pointers are already moved to BaseHadronizer)
    //
    enum { PP, PPbar, ElectronPositron, ElectronProton, PositronProton };

    int fInitialState;  // pp, ppbar, e-e+, or e-p
    double fCOMEnergy;  // this one is irrelevant for setting py6 as hadronizer
                        // or if anything, it should be picked up from LHERunInfoProduct !
    double fBeam1PZ;
    double fBeam2PZ;

    static JetMatching* fJetMatching;

    bool fHepMCVerbosity;
    unsigned int fMaxEventsToPrint;

    // this is the only one specific to Pythia6
    //
    unsigned int fPythiaListVerbosity;
    bool fDisplayPythiaBanner;
    bool fDisplayPythiaCards;

    // these two params control stop- and r-handron features,
    // that are "custom" add-ons to Py6;
    // I doubt they should drag along Py6Int main library...
    //
    bool fStopHadronsEnabled;
    bool fGluinoHadronsEnabled;

    // this is a "trick" to generate enriched mu-samples and the likes
    bool fImposeProperTime;

    // and final touch - conversion of Py6 PID's into PDG convension
    bool fConvertToPDG;

    // tmp stuff, to deal with EvtGen corrupting pyjets
    // int NPartsBeforeDecays;

    static const std::vector<std::string> theSharedResources;
  };
}  // namespace gen

#endif