Inheritance diagram for Pythia8Hadronizer:

Public Member Functions
const char *	classname () const override

void	finalizeEvent () override

bool	generatePartonsAndHadronize () override

std::unique_ptr< GenLumiInfoHeader >	getGenLumiInfoHeader () const override

bool	hadronize ()

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons () override

	Pythia8Hadronizer (const edm::ParameterSet &params)

virtual bool	residualDecay ()

void	statistics () override

	~Pythia8Hadronizer () override

Public Member Functions inherited from gen::Py8InterfaceBase
bool	decay ()

bool	declareSpecialSettings (const std::vector< std::string > &)

bool	declareStableParticles (const std::vector< int > &)

void	p8SetRandomEngine (CLHEP::HepRandomEngine *v)

	Py8InterfaceBase (edm::ParameterSet const &ps)

P8RndmEngine &	randomEngine ()

bool	readSettings (int)

	~Py8InterfaceBase () override

Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)

void	cleanLHE ()

void	generateLHE (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine, unsigned int ncpu)

edm::Event &	getEDMEvent () const

std::unique_ptr< HepMC::GenEvent >	getGenEvent ()

std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()

GenRunInfoProduct &	getGenRunInfo ()

std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()

const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const

const std::string &	gridpackPath () const

int	randomIndex () const

const std::string &	randomInitConfigDescription () const

void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)

void	resetEvent (std::unique_ptr< HepMC::GenEvent > event)

void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)

virtual bool	select (HepMC::GenEvent *) const

void	setEDMEvent (edm::Event &event)

void	setLHEEvent (std::unique_ptr< lhef::LHEEvent > event)

void	setLHERunInfo (std::unique_ptr< lhef::LHERunInfo > runInfo)

void	setRandomEngine (CLHEP::HepRandomEngine *v)

std::vector< std::string > const &	sharedResources () const

virtual	~BaseHadronizer () noexcept(false)

Private Types
enum	{ PP, PPbar, ElectronPositron }

Private Member Functions
void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

const std::vector< std::string > &	doSharedResources () const override

Private Attributes
double	comEnergy
	Center-of-Mass energy. More...

vector< float >	DJR

int	EV1_emittedMode

int	EV1_maxVetoCount

bool	EV1_MPIvetoOn

int	EV1_nFinal

int	EV1_nFinalMode

int	EV1_pTdefMode

int	EV1_pTempMode

int	EV1_pThardMode

int	EV1_QEDvetoMode

bool	EV1_vetoOn

double	fBeam1PZ

double	fBeam2PZ

std::unique_ptr< BiasedTauDecayer >	fBiasedTauDecayer

std::unique_ptr< Pythia8::Dire >	fDire

std::unique_ptr< PowhegHooks >	fEmissionVetoHook

std::unique_ptr< EmissionVetoHook1 >	fEmissionVetoHook1

int	fInitialState

std::unique_ptr< JetMatchingHook >	fJetMatchingHook

std::unique_ptr< Pythia8::JetMatchingMadgraph >	fJetMatchingPy8InternalHook

std::unique_ptr< Pythia8::amcnlo_unitarised_interface >	fMergingHook

std::unique_ptr< MultiUserHook >	fMultiUserHook

std::unique_ptr< PowhegHooksBB4L >	fPowhegHooksBB4L

std::unique_ptr< PowhegResHook >	fPowhegResHook

std::unique_ptr< PTFilterHook >	fPTFilterHook

std::unique_ptr< ResonanceDecayFilterHook >	fResonanceDecayFilterHook

std::unique_ptr< UserHooks >	fReweightEmpUserHook

std::unique_ptr< UserHooks >	fReweightPtHatRapUserHook

std::unique_ptr< UserHooks >	fReweightRapUserHook

std::unique_ptr< UserHooks >	fReweightUserHook

std::unique_ptr< Pythia8::VinciaPlugin >	fvincia

std::unique_ptr< LHAupLesHouches >	lhaUP

std::string	LHEInputFileName

int	nFSRveto

int	nISRveto

int	nME

int	nMEFiltered

Static Private Attributes
static const std::vector< std::string >	p8SharedResources = {edm::SharedResourceNames::kPythia8}

Additional Inherited Members
Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > &	event ()

std::unique_ptr< GenEventInfoProduct > &	eventInfo ()

lhef::LHEEvent *	lheEvent ()

lhef::LHERunInfo *	lheRunInfo ()

GenRunInfoProduct &	runInfo ()

Protected Attributes inherited from gen::Py8InterfaceBase
HepMC::IO_AsciiParticles *	ascii_io

std::shared_ptr< EvtGenDecays >	evtgenDecays

std::string	evtgenDecFile

std::string	evtgenPdlFile

std::vector< std::string >	evtgenUserFiles

std::unique_ptr< Pythia8::Pythia >	fDecayer

std::unique_ptr< Pythia8::Pythia >	fMasterGen

edm::ParameterSet	fParameters

unsigned int	maxEventsToPrint

bool	pythiaHepMCVerbosity

bool	pythiaHepMCVerbosityParticles

unsigned int	pythiaPylistVerbosity

std::string	slhafile_

HepMC::Pythia8ToHepMC	toHepMC

bool	useEvtGen

Protected Attributes inherited from gen::BaseHadronizer
std::string	lheFile_

int	randomIndex_

Detailed Description

Definition at line 79 of file Pythia8Hadronizer.cc.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

private

Enumerator
PP
PPbar
ElectronPositron

Definition at line 113 of file Pythia8Hadronizer.cc.

113 { PP, PPbar, ElectronPositron };

Constructor & Destructor Documentation

◆ Pythia8Hadronizer()

Pythia8Hadronizer::Pythia8Hadronizer ( const edm::ParameterSet & params )

Definition at line 176 of file Pythia8Hadronizer.cc.

     : Py8InterfaceBase(params),
       comEnergy(params.getParameter<double>("comEnergy")),
       LHEInputFileName(params.getUntrackedParameter<std::string>("LHEInputFileName", "")),
       fInitialState(PP),
       nME(-1),
       nMEFiltered(-1),
       nISRveto(0),
       nFSRveto(0) {
   // J.Y.: the following 3 parameters are hacked "for a reason"
   //
   if (params.exists("PPbarInitialState")) {
     if (fInitialState == PP) {
       fInitialState = PPbar;
       edm::LogImportant("GeneratorInterface|Pythia8Interface")
           << "Pythia8 will be initialized for PROTON-ANTIPROTON INITIAL STATE. "
           << "This is a user-request change from the DEFAULT PROTON-PROTON initial state.";
     } else {
       // probably need to throw on attempt to override ?
     }
   } else if (params.exists("ElectronPositronInitialState")) {
     if (fInitialState == PP) {
       fInitialState = ElectronPositron;
       edm::LogInfo("GeneratorInterface|Pythia8Interface")
           << "Pythia8 will be initialized for ELECTRON-POSITRON INITIAL STATE. "
           << "This is a user-request change from the DEFAULT PROTON-PROTON initial state.";
     } else {
       // probably need to throw on attempt to override ?
     }
   } else if (params.exists("ElectronProtonInitialState") || params.exists("PositronProtonInitialState")) {
     // throw on unknown initial state !
     throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
         << " UNKNOWN INITIAL STATE. \n The allowed initial states are: PP, PPbar, ElectronPositron \n";
   }
  
   // Reweight user hook
   //
   if (params.exists("reweightGen")) {
     edm::LogInfo("Pythia8Interface") << "Start setup for reweightGen";
     edm::ParameterSet rgParams = params.getParameter<edm::ParameterSet>("reweightGen");
     fReweightUserHook.reset(
         new PtHatReweightUserHook(rgParams.getParameter<double>("pTRef"), rgParams.getParameter<double>("power")));
     edm::LogInfo("Pythia8Interface") << "End setup for reweightGen";
   }
   if (params.exists("reweightGenEmp")) {
     edm::LogInfo("Pythia8Interface") << "Start setup for reweightGenEmp";
     edm::ParameterSet rgeParams = params.getParameter<edm::ParameterSet>("reweightGenEmp");
  
     std::string tuneName = "";
     if (rgeParams.exists("tune"))
       tuneName = rgeParams.getParameter<std::string>("tune");
     fReweightEmpUserHook.reset(new PtHatEmpReweightUserHook(tuneName));
     edm::LogInfo("Pythia8Interface") << "End setup for reweightGenEmp";
   }
   if (params.exists("reweightGenRap")) {
     edm::LogInfo("Pythia8Interface") << "Start setup for reweightGenRap";
     edm::ParameterSet rgrParams = params.getParameter<edm::ParameterSet>("reweightGenRap");
     fReweightRapUserHook.reset(new RapReweightUserHook(rgrParams.getParameter<std::string>("yLabSigmaFunc"),
                                                        rgrParams.getParameter<double>("yLabPower"),
                                                        rgrParams.getParameter<std::string>("yCMSigmaFunc"),
                                                        rgrParams.getParameter<double>("yCMPower"),
                                                        rgrParams.getParameter<double>("pTHatMin"),
                                                        rgrParams.getParameter<double>("pTHatMax")));
     edm::LogInfo("Pythia8Interface") << "End setup for reweightGenRap";
   }
   if (params.exists("reweightGenPtHatRap")) {
     edm::LogInfo("Pythia8Interface") << "Start setup for reweightGenPtHatRap";
     edm::ParameterSet rgrParams = params.getParameter<edm::ParameterSet>("reweightGenPtHatRap");
     fReweightPtHatRapUserHook.reset(new PtHatRapReweightUserHook(rgrParams.getParameter<std::string>("yLabSigmaFunc"),
                                                                  rgrParams.getParameter<double>("yLabPower"),
                                                                  rgrParams.getParameter<std::string>("yCMSigmaFunc"),
                                                                  rgrParams.getParameter<double>("yCMPower"),
                                                                  rgrParams.getParameter<double>("pTHatMin"),
                                                                  rgrParams.getParameter<double>("pTHatMax")));
     edm::LogInfo("Pythia8Interface") << "End setup for reweightGenPtHatRap";
   }
  
   if (params.exists("useUserHook"))
     throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
         << " Obsolete parameter: useUserHook \n Please use the actual one instead \n";
  
   // PS matching prototype
   //
   if (params.exists("jetMatching")) {
     edm::ParameterSet jmParams = params.getUntrackedParameter<edm::ParameterSet>("jetMatching");
     std::string scheme = jmParams.getParameter<std::string>("scheme");
     if (scheme == "Madgraph" || scheme == "MadgraphFastJet") {
       fJetMatchingHook.reset(new JetMatchingHook(jmParams, &fMasterGen->info));
     }
   }
  
   // Pythia8Interface emission veto
   //
   if (params.exists("emissionVeto1")) {
     EV1_nFinal = -1;
     if (params.exists("EV1_nFinal"))
       EV1_nFinal = params.getParameter<int>("EV1_nFinal");
     EV1_vetoOn = true;
     if (params.exists("EV1_vetoOn"))
       EV1_vetoOn = params.getParameter<bool>("EV1_vetoOn");
     EV1_maxVetoCount = 10;
     if (params.exists("EV1_maxVetoCount"))
       EV1_maxVetoCount = params.getParameter<int>("EV1_maxVetoCount");
     EV1_pThardMode = 1;
     if (params.exists("EV1_pThardMode"))
       EV1_pThardMode = params.getParameter<int>("EV1_pThardMode");
     EV1_pTempMode = 0;
     if (params.exists("EV1_pTempMode"))
       EV1_pTempMode = params.getParameter<int>("EV1_pTempMode");
     if (EV1_pTempMode > 2 || EV1_pTempMode < 0)
       throw edm::Exception(edm::errors::Configuration, "Pythia8Interface") << " Wrong value for EV1_pTempMode code\n";
     EV1_emittedMode = 0;
     if (params.exists("EV1_emittedMode"))
       EV1_emittedMode = params.getParameter<int>("EV1_emittedMode");
     EV1_pTdefMode = 1;
     if (params.exists("EV1_pTdefMode"))
       EV1_pTdefMode = params.getParameter<int>("EV1_pTdefMode");
     EV1_MPIvetoOn = false;
     if (params.exists("EV1_MPIvetoOn"))
       EV1_MPIvetoOn = params.getParameter<bool>("EV1_MPIvetoOn");
     EV1_QEDvetoMode = 0;
     if (params.exists("EV1_QEDvetoMode"))
       EV1_QEDvetoMode = params.getParameter<int>("EV1_QEDvetoMode");
     EV1_nFinalMode = 0;
     if (params.exists("EV1_nFinalMode"))
       EV1_nFinalMode = params.getParameter<int>("EV1_nFinalMode");
     fEmissionVetoHook1.reset(new EmissionVetoHook1(EV1_nFinal,
                                                    EV1_vetoOn,
                                                    EV1_maxVetoCount,
                                                    EV1_pThardMode,
                                                    EV1_pTempMode,
                                                    EV1_emittedMode,
                                                    EV1_pTdefMode,
                                                    EV1_MPIvetoOn,
                                                    EV1_QEDvetoMode,
                                                    EV1_nFinalMode,
                                                    0));
   }
  
   if (params.exists("VinciaPlugin")) {
     fMasterGen.reset(new Pythia);
     fvincia.reset(new Pythia8::VinciaPlugin(fMasterGen.get()));
   }
   if (params.exists("DirePlugin")) {
     fMasterGen.reset(new Pythia);
     fDire.reset(new Pythia8::Dire());
     fDire->initSettings(*fMasterGen.get());
     fDire->initShowersAndWeights(*fMasterGen.get(), nullptr, nullptr);
   }
 }

References edm::errors::Configuration, ElectronPositron, EV1_emittedMode, EV1_maxVetoCount, EV1_MPIvetoOn, EV1_nFinal, EV1_nFinalMode, EV1_pTdefMode, EV1_pTempMode, EV1_pThardMode, EV1_QEDvetoMode, EV1_vetoOn, Exception, edm::ParameterSet::exists(), fDire, fEmissionVetoHook1, fInitialState, fJetMatchingHook, gen::Py8InterfaceBase::fMasterGen, fReweightEmpUserHook, fReweightPtHatRapUserHook, fReweightRapUserHook, fReweightUserHook, fvincia, edm::ParameterSet::getParameter(), CalibrationSummaryClient_cfi::params, PP, PPbar, generator_cfi::scheme, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~Pythia8Hadronizer()

Pythia8Hadronizer::~Pythia8Hadronizer ( )

override

Definition at line 327 of file Pythia8Hadronizer.cc.

327 {}

Member Function Documentation

◆ classname()

const char* Pythia8Hadronizer::classname ( ) const

inlineoverridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 96 of file Pythia8Hadronizer.cc.

96 { return "Pythia8Hadronizer"; }

◆ doSetRandomEngine()

void Pythia8Hadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine * v )

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 104 of file Pythia8Hadronizer.cc.

104 { p8SetRandomEngine(v); }

References findQualityFiles::v.

◆ doSharedResources()

const std::vector<std::string>& Pythia8Hadronizer::doSharedResources ( ) const

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 105 of file Pythia8Hadronizer.cc.

105 { return p8SharedResources; }

◆ finalizeEvent()

void Pythia8Hadronizer::finalizeEvent ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 892 of file Pythia8Hadronizer.cc.

                                       {
   bool lhe = lheEvent() != nullptr;
  
   // now create the GenEventInfo product from the GenEvent and fill
   // the missing pieces
   eventInfo().reset(new GenEventInfoProduct(event().get()));
  
   // in pythia pthat is used to subdivide samples into different bins
   // in LHE mode the binning is done by the external ME generator
   // which is likely not pthat, so only filling it for Py6 internal mode
   if (!lhe) {
     eventInfo()->setBinningValues(std::vector<double>(1, fMasterGen->info.pTHat()));
   }
  
   eventInfo()->setDJR(DJR);
   eventInfo()->setNMEPartons(nME);
   eventInfo()->setNMEPartonsFiltered(nMEFiltered);
  
   //******** Verbosity ********
  
   if (maxEventsToPrint > 0 && (pythiaPylistVerbosity || pythiaHepMCVerbosity || pythiaHepMCVerbosityParticles)) {
     maxEventsToPrint--;
     if (pythiaPylistVerbosity) {
       fMasterGen->info.list();
       fMasterGen->event.list();
     }
  
     if (pythiaHepMCVerbosity) {
       std::cout << "Event process = " << fMasterGen->info.code() << "\n"
                 << "----------------------" << std::endl;
       event()->print();
     }
     if (pythiaHepMCVerbosityParticles) {
       std::cout << "Event process = " << fMasterGen->info.code() << "\n"
                 << "----------------------" << std::endl;
       ascii_io->write_event(event().get());
     }
   }
 }

References gen::Py8InterfaceBase::ascii_io, gather_cfg::cout, DJR, gen::BaseHadronizer::event(), gen::BaseHadronizer::eventInfo(), gen::Py8InterfaceBase::fMasterGen, get, gen::BaseHadronizer::lheEvent(), gen::Py8InterfaceBase::maxEventsToPrint, nME, nMEFiltered, gen::Py8InterfaceBase::pythiaHepMCVerbosity, gen::Py8InterfaceBase::pythiaHepMCVerbosityParticles, and gen::Py8InterfaceBase::pythiaPylistVerbosity.

◆ generatePartonsAndHadronize()

bool Pythia8Hadronizer::generatePartonsAndHadronize ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 662 of file Pythia8Hadronizer.cc.

                                                     {
   DJR.resize(0);
   nME = -1;
   nMEFiltered = -1;
  
   if (fJetMatchingHook.get()) {
     fJetMatchingHook->resetMatchingStatus();
     fJetMatchingHook->beforeHadronization(lheEvent());
   }
  
   if (!fMasterGen->next())
     return false;
  
   double mergeweight = fMasterGen.get()->info.mergingWeightNLO();
   if (fMergingHook.get()) {
     mergeweight *= fMergingHook->getNormFactor();
   }
  
   //protect against 0-weight from ckkw or similar
   if (std::abs(mergeweight) == 0.) {
     event().reset();
     return false;
   }
  
   if (fJetMatchingPy8InternalHook.get()) {
     const std::vector<double> djrmatch = fJetMatchingPy8InternalHook->getDJR();
     //cap size of djr vector to save storage space (keep only up to first 6 elements)
     unsigned int ndjr = std::min(djrmatch.size(), std::vector<double>::size_type(6));
     for (unsigned int idjr = 0; idjr < ndjr; ++idjr) {
       DJR.push_back(djrmatch[idjr]);
     }
  
     nME = fJetMatchingPy8InternalHook->nMEpartons().first;
     nMEFiltered = fJetMatchingPy8InternalHook->nMEpartons().second;
   }
  
   if (evtgenDecays.get())
     evtgenDecays->decay();
  
   event().reset(new HepMC::GenEvent);
   bool py8hepmc = toHepMC.fill_next_event(*(fMasterGen.get()), event().get());
  
   if (!py8hepmc) {
     return false;
   }
  
   //add ckkw/umeps/unlops merging weight
   if (mergeweight != 1.) {
     event()->weights()[0] *= mergeweight;
   }
  
   if (fEmissionVetoHook.get()) {
     nISRveto += fEmissionVetoHook->getNISRveto();
     nFSRveto += fEmissionVetoHook->getNFSRveto();
   }
  
   //fill additional weights for systematic uncertainties
   if (fMasterGen->info.getWeightsDetailedSize() > 0) {
     for (const string &key : fMasterGen->info.initrwgt->weightsKeys) {
       double wgt = (*fMasterGen->info.weights_detailed)[key];
       event()->weights().push_back(wgt);
     }
   } else if (fMasterGen->info.getWeightsCompressedSize() > 0) {
     for (unsigned int i = 0; i < fMasterGen->info.getWeightsCompressedSize(); i++) {
       double wgt = fMasterGen->info.getWeightsCompressedValue(i);
       event()->weights().push_back(wgt);
     }
   }
  
   // fill shower weights
   // http://home.thep.lu.se/~torbjorn/pythia82html/Variations.html
   if (fMasterGen->info.nWeights() > 1) {
     for (int i = 0; i < fMasterGen->info.nWeights(); ++i) {
       double wgt = fMasterGen->info.weight(i);
       event()->weights().push_back(wgt);
     }
   }
  
   // VINCIA shower weights
   // http://vincia.hepforge.org/current/share/Vincia/htmldoc/VinciaUncertainties.html
   if (fvincia.get()) {
     event()->weights()[0] *= fvincia->weight(0);
     for (int iVar = 1; iVar < fvincia->nWeights(); iVar++) {
       event()->weights().push_back(fvincia->weight(iVar));
     }
   }
  
   // Retrieve Dire shower weights
   if (fDire.get()) {
     fDire->weightsPtr->calcWeight(0.);
     fDire->weightsPtr->reset();
  
     //Make sure the base weight comes first
     event()->weights()[0] *= fDire->weightsPtr->getShowerWeight("base");
  
     unordered_map<string, double>::iterator it;
     for (it = fDire->weightsPtr->getShowerWeights()->begin(); it != fDire->weightsPtr->getShowerWeights()->end();
          it++) {
       if (it->first == "base")
         continue;
       event()->weights().push_back(it->second);
     }
   }
  
   return true;
 }

References funct::abs(), DJR, gen::BaseHadronizer::event(), gen::Py8InterfaceBase::evtgenDecays, fDire, fEmissionVetoHook, fJetMatchingHook, fJetMatchingPy8InternalHook, gen::Py8InterfaceBase::fMasterGen, fMergingHook, fvincia, mps_fire::i, crabWrapper::key, gen::BaseHadronizer::lheEvent(), min(), nFSRveto, nISRveto, nME, nMEFiltered, and gen::Py8InterfaceBase::toHepMC.

◆ getGenLumiInfoHeader()

std::unique_ptr< GenLumiInfoHeader > Pythia8Hadronizer::getGenLumiInfoHeader ( ) const

overridevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 932 of file Pythia8Hadronizer.cc.

                                                                                {
   auto genLumiInfoHeader = BaseHadronizer::getGenLumiInfoHeader();
  
   //fill lhe headers
   //*FIXME* initrwgt header is corrupt due to pythia bug
   for (const std::string &key : fMasterGen->info.headerKeys()) {
     genLumiInfoHeader->lheHeaders().emplace_back(key, fMasterGen->info.header(key));
   }
  
   //check, if it is not only nominal weight
   int weights_number = fMasterGen->info.nWeights();
   if (fMasterGen->info.initrwgt)
     weights_number += fMasterGen->info.initrwgt->weightsKeys.size();
   if (weights_number > 1) {
     genLumiInfoHeader->weightNames().reserve(weights_number + 1);
     genLumiInfoHeader->weightNames().push_back("nominal");
   }
  
   //fill weight names
   if (fMasterGen->info.initrwgt) {
     for (const std::string &key : fMasterGen->info.initrwgt->weightsKeys) {
       std::string weightgroupname;
       for (const auto &wgtgrp : fMasterGen->info.initrwgt->weightgroups) {
         const auto &wgtgrpwgt = wgtgrp.second.weights.find(key);
         if (wgtgrpwgt != wgtgrp.second.weights.end()) {
           weightgroupname = wgtgrp.first;
         }
       }
  
       std::ostringstream weightname;
       weightname << "LHE, id = " << key << ", ";
       if (!weightgroupname.empty()) {
         weightname << "group = " << weightgroupname << ", ";
       }
       weightname << fMasterGen->info.initrwgt->weights[key].contents;
       genLumiInfoHeader->weightNames().push_back(weightname.str());
     }
   }
  
   //fill shower labels
   // http://home.thep.lu.se/~torbjorn/pythia82html/Variations.html
   // http://home.thep.lu.se/~torbjorn/doxygen/classPythia8_1_1Info.html
   if (fMasterGen->info.nWeights() > 1) {
     for (int i = 0; i < fMasterGen->info.nWeights(); ++i) {
       genLumiInfoHeader->weightNames().push_back(fMasterGen->info.weightLabel(i));
     }
   }
  
   // VINCIA shower weights
   // http://vincia.hepforge.org/current/share/Vincia/htmldoc/VinciaUncertainties.html
   if (fvincia.get()) {
     for (int iVar = 0; iVar < fvincia->nWeights(); iVar++) {
       genLumiInfoHeader->weightNames().push_back(fvincia->weightLabel(iVar));
     }
   }
  
   if (fDire.get()) {
     //Make sure the base weight comes first
     genLumiInfoHeader->weightNames().push_back("base");
  
     unordered_map<string, double>::iterator it;
     for (it = fDire->weightsPtr->getShowerWeights()->begin(); it != fDire->weightsPtr->getShowerWeights()->end();
          it++) {
       if (it->first == "base")
         continue;
       genLumiInfoHeader->weightNames().push_back(it->first);
     }
   }
  
   return genLumiInfoHeader;
 }

References fDire, gen::Py8InterfaceBase::fMasterGen, fvincia, nano_cff::genLumiInfoHeader, mps_fire::i, crabWrapper::key, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ hadronize()

bool Pythia8Hadronizer::hadronize ( )

Definition at line 769 of file Pythia8Hadronizer.cc.

                                   {
   DJR.resize(0);
   nME = -1;
   nMEFiltered = -1;
   if (LHEInputFileName.empty())
     lhaUP->loadEvent(lheEvent());
  
   if (fJetMatchingHook.get()) {
     fJetMatchingHook->resetMatchingStatus();
     fJetMatchingHook->beforeHadronization(lheEvent());
   }
  
   bool py8next = fMasterGen->next();
  
   double mergeweight = fMasterGen.get()->info.mergingWeightNLO();
   if (fMergingHook.get()) {
     mergeweight *= fMergingHook->getNormFactor();
   }
  
   //protect against 0-weight from ckkw or similar
   if (!py8next || std::abs(mergeweight) == 0.) {
     lheEvent()->count(lhef::LHERunInfo::kSelected, 1.0, mergeweight);
     event().reset();
     return false;
   }
  
   if (fJetMatchingPy8InternalHook.get()) {
     const std::vector<double> djrmatch = fJetMatchingPy8InternalHook->getDJR();
     //cap size of djr vector to save storage space (keep only up to first 6 elements)
     unsigned int ndjr = std::min(djrmatch.size(), std::vector<double>::size_type(6));
     for (unsigned int idjr = 0; idjr < ndjr; ++idjr) {
       DJR.push_back(djrmatch[idjr]);
     }
  
     nME = fJetMatchingPy8InternalHook->nMEpartons().first;
     nMEFiltered = fJetMatchingPy8InternalHook->nMEpartons().second;
   }
  
   // update LHE matching statistics
   //
   lheEvent()->count(lhef::LHERunInfo::kAccepted, 1.0, mergeweight);
  
   if (evtgenDecays.get())
     evtgenDecays->decay();
  
   event().reset(new HepMC::GenEvent);
   bool py8hepmc = toHepMC.fill_next_event(*(fMasterGen.get()), event().get());
   if (!py8hepmc) {
     return false;
   }
  
   //add ckkw/umeps/unlops merging weight
   if (mergeweight != 1.) {
     event()->weights()[0] *= mergeweight;
   }
  
   if (fEmissionVetoHook.get()) {
     nISRveto += fEmissionVetoHook->getNISRveto();
     nFSRveto += fEmissionVetoHook->getNFSRveto();
   }
  
   // fill shower weights
   // http://home.thep.lu.se/~torbjorn/pythia82html/Variations.html
   if (fMasterGen->info.nWeights() > 1) {
     for (int i = 0; i < fMasterGen->info.nWeights(); ++i) {
       double wgt = fMasterGen->info.weight(i);
       event()->weights().push_back(wgt);
     }
   }
  
   return true;
 }

References funct::abs(), lhef::LHEEvent::count(), DJR, gen::BaseHadronizer::event(), gen::Py8InterfaceBase::evtgenDecays, fEmissionVetoHook, fJetMatchingHook, fJetMatchingPy8InternalHook, gen::Py8InterfaceBase::fMasterGen, fMergingHook, mps_fire::i, lhef::LHERunInfo::kAccepted, lhef::LHERunInfo::kSelected, lhaUP, gen::BaseHadronizer::lheEvent(), LHEInputFileName, min(), nFSRveto, nISRveto, nME, nMEFiltered, and gen::Py8InterfaceBase::toHepMC.

◆ initializeForExternalPartons()

bool Pythia8Hadronizer::initializeForExternalPartons ( )

Definition at line 493 of file Pythia8Hadronizer.cc.

                                                      {
   edm::LogInfo("Pythia8Interface") << "Initializing for external partons";
  
   bool status = false, status1 = false;
  
   fMultiUserHook.reset(new MultiUserHook);
  
   if (fReweightUserHook.get())
     fMultiUserHook->addHook(fReweightUserHook.get());
   if (fReweightEmpUserHook.get())
     fMultiUserHook->addHook(fReweightEmpUserHook.get());
   if (fReweightRapUserHook.get())
     fMultiUserHook->addHook(fReweightRapUserHook.get());
   if (fReweightPtHatRapUserHook.get())
     fMultiUserHook->addHook(fReweightPtHatRapUserHook.get());
   if (fJetMatchingHook.get())
     fMultiUserHook->addHook(fJetMatchingHook.get());
   if (fEmissionVetoHook1.get()) {
     edm::LogInfo("Pythia8Interface") << "Turning on Emission Veto Hook 1 from CMSSW Pythia8Interface";
     fMultiUserHook->addHook(fEmissionVetoHook1.get());
   }
  
   if (fMasterGen->settings.mode("POWHEG:veto") > 0 || fMasterGen->settings.mode("POWHEG:MPIveto") > 0) {
     if (fJetMatchingHook.get() || fEmissionVetoHook1.get())
       throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
           << " Attempt to turn on PowhegHooks by pythia8 settings but there are incompatible hooks on \n Incompatible "
              "are : jetMatching, emissionVeto1 \n";
  
     fEmissionVetoHook.reset(new PowhegHooks());
  
     edm::LogInfo("Pythia8Interface") << "Turning on Emission Veto Hook from pythia8 code";
     fMultiUserHook->addHook(fEmissionVetoHook.get());
   }
  
   bool PowhegRes = fMasterGen->settings.flag("POWHEGres:calcScales");
   if (PowhegRes) {
     edm::LogInfo("Pythia8Interface") << "Turning on resonance scale setting from CMSSW Pythia8Interface";
     fPowhegResHook.reset(new PowhegResHook());
     fMultiUserHook->addHook(fPowhegResHook.get());
   }
  
   bool PowhegBB4L = fMasterGen->settings.flag("POWHEG:bb4l");
   if (PowhegBB4L) {
     edm::LogInfo("Pythia8Interface") << "Turning on BB4l hook from CMSSW Pythia8Interface";
     fPowhegHooksBB4L.reset(new PowhegHooksBB4L());
     fMultiUserHook->addHook(fPowhegHooksBB4L.get());
   }
  
   //adapted from main89.cc in pythia8 examples
   bool internalMatching = fMasterGen->settings.flag("JetMatching:merge");
   bool internalMerging = !(fMasterGen->settings.word("Merging:Process") == "void");
  
   if (internalMatching && internalMerging) {
     throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
         << " Only one jet matching/merging scheme can be used at a time. \n";
   }
  
   if (internalMatching) {
     fJetMatchingPy8InternalHook.reset(new Pythia8::JetMatchingMadgraph);
     fMultiUserHook->addHook(fJetMatchingPy8InternalHook.get());
   }
  
   if (internalMerging) {
     int scheme = (fMasterGen->settings.flag("Merging:doUMEPSTree") || fMasterGen->settings.flag("Merging:doUMEPSSubt"))
                      ? 1
                      : ((fMasterGen->settings.flag("Merging:doUNLOPSTree") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSSubt") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSLoop") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSSubtNLO"))
                             ? 2
                             : 0);
     fMergingHook.reset(new Pythia8::amcnlo_unitarised_interface(scheme));
     fMultiUserHook->addHook(fMergingHook.get());
   }
  
   bool biasedTauDecayer = fMasterGen->settings.flag("BiasedTauDecayer:filter");
   if (biasedTauDecayer) {
     fBiasedTauDecayer.reset(new BiasedTauDecayer(&(fMasterGen->info),
                                                  &(fMasterGen->settings),
                                                  &(fMasterGen->particleData),
                                                  &(fMasterGen->rndm),
                                                  &(fMasterGen->couplings)));
     std::vector<int> handledParticles;
     handledParticles.push_back(15);
     fMasterGen->setDecayPtr(fBiasedTauDecayer.get(), handledParticles);
   }
  
   bool resonanceDecayFilter = fMasterGen->settings.flag("ResonanceDecayFilter:filter");
   if (resonanceDecayFilter) {
     fResonanceDecayFilterHook.reset(new ResonanceDecayFilterHook);
     fMultiUserHook->addHook(fResonanceDecayFilterHook.get());
   }
  
   bool PTFilter = fMasterGen->settings.flag("PTFilter:filter");
   if (PTFilter) {
     fPTFilterHook.reset(new PTFilterHook);
     fMultiUserHook->addHook(fPTFilterHook.get());
   }
  
   if (fMultiUserHook->nHooks() > 0) {
     fMasterGen->setUserHooksPtr(fMultiUserHook.get());
   }
  
   if (!LHEInputFileName.empty()) {
     edm::LogInfo("Pythia8Interface") << "Initialize direct pythia8 reading from LHE file " << LHEInputFileName;
     edm::LogInfo("Pythia8Interface") << "Some LHE information can be not stored";
     fMasterGen->settings.mode("Beams:frameType", 4);
     fMasterGen->settings.word("Beams:LHEF", LHEInputFileName);
     status = fMasterGen->init();
  
   } else {
     lhaUP.reset(new LHAupLesHouches());
     lhaUP->setScalesFromLHEF(fMasterGen->settings.flag("Beams:setProductionScalesFromLHEF"));
     lhaUP->loadRunInfo(lheRunInfo());
  
     if (fJetMatchingHook.get()) {
       fJetMatchingHook->init(lheRunInfo());
     }
  
     fMasterGen->settings.mode("Beams:frameType", 5);
     fMasterGen->setLHAupPtr(lhaUP.get());
     edm::LogInfo("Pythia8Interface") << "Initializing MasterGen";
     status = fMasterGen->init();
   }
  
   //clean up temp file
   if (!slhafile_.empty()) {
     std::remove(slhafile_.c_str());
   }
  
   if (pythiaPylistVerbosity > 10) {
     if (pythiaPylistVerbosity == 11 || pythiaPylistVerbosity == 13)
       fMasterGen->settings.listAll();
     if (pythiaPylistVerbosity == 12 || pythiaPylistVerbosity == 13)
       fMasterGen->particleData.listAll();
   }
  
   // init decayer
   fDecayer->settings.flag("ProcessLevel:all", false);  // trick
   fDecayer->settings.flag("ProcessLevel:resonanceDecays", true);
   edm::LogInfo("Pythia8Interface") << "Initializing Decayer";
   status1 = fDecayer->init();
  
   if (useEvtGen) {
     edm::LogInfo("Pythia8Hadronizer") << "Creating and initializing pythia8 EvtGen plugin";
     evtgenDecays.reset(new EvtGenDecays(fMasterGen.get(), evtgenDecFile, evtgenPdlFile));
     for (unsigned int i = 0; i < evtgenUserFiles.size(); i++)
       evtgenDecays->readDecayFile(evtgenUserFiles.at(i));
   }
  
   return (status && status1);
 }

◆ initializeForInternalPartons()

bool Pythia8Hadronizer::initializeForInternalPartons ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 329 of file Pythia8Hadronizer.cc.

                                                      {
   bool status = false, status1 = false;
  
   if (lheFile_.empty()) {
     if (fInitialState == PP)  // default
     {
       fMasterGen->settings.mode("Beams:idA", 2212);
       fMasterGen->settings.mode("Beams:idB", 2212);
     } else if (fInitialState == PPbar) {
       fMasterGen->settings.mode("Beams:idA", 2212);
       fMasterGen->settings.mode("Beams:idB", -2212);
     } else if (fInitialState == ElectronPositron) {
       fMasterGen->settings.mode("Beams:idA", 11);
       fMasterGen->settings.mode("Beams:idB", -11);
     } else {
       // throw on unknown initial state !
       throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
           << " UNKNOWN INITIAL STATE. \n The allowed initial states are: PP, PPbar, ElectronPositron \n";
     }
     fMasterGen->settings.parm("Beams:eCM", comEnergy);
   } else {
     fMasterGen->settings.mode("Beams:frameType", 4);
     fMasterGen->settings.word("Beams:LHEF", lheFile_);
   }
  
   fMultiUserHook.reset(new MultiUserHook);
  
   if (fReweightUserHook.get())
     fMultiUserHook->addHook(fReweightUserHook.get());
   if (fReweightEmpUserHook.get())
     fMultiUserHook->addHook(fReweightEmpUserHook.get());
   if (fReweightRapUserHook.get())
     fMultiUserHook->addHook(fReweightRapUserHook.get());
   if (fReweightPtHatRapUserHook.get())
     fMultiUserHook->addHook(fReweightPtHatRapUserHook.get());
   if (fJetMatchingHook.get())
     fMultiUserHook->addHook(fJetMatchingHook.get());
   if (fEmissionVetoHook1.get()) {
     edm::LogInfo("Pythia8Interface") << "Turning on Emission Veto Hook 1 from CMSSW Pythia8Interface";
     fMultiUserHook->addHook(fEmissionVetoHook1.get());
   }
  
   if (fMasterGen->settings.mode("POWHEG:veto") > 0 || fMasterGen->settings.mode("POWHEG:MPIveto") > 0) {
     if (fJetMatchingHook.get() || fEmissionVetoHook1.get())
       throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
           << " Attempt to turn on PowhegHooks by pythia8 settings but there are incompatible hooks on \n Incompatible "
              "are : jetMatching, emissionVeto1 \n";
  
     fEmissionVetoHook.reset(new PowhegHooks());
  
     edm::LogInfo("Pythia8Interface") << "Turning on Emission Veto Hook from pythia8 code";
     fMultiUserHook->addHook(fEmissionVetoHook.get());
   }
  
   bool PowhegRes = fMasterGen->settings.flag("POWHEGres:calcScales");
   if (PowhegRes) {
     edm::LogInfo("Pythia8Interface") << "Turning on resonance scale setting from CMSSW Pythia8Interface";
     fPowhegResHook.reset(new PowhegResHook());
     fMultiUserHook->addHook(fPowhegResHook.get());
   }
  
   bool PowhegBB4L = fMasterGen->settings.flag("POWHEG:bb4l");
   if (PowhegBB4L) {
     edm::LogInfo("Pythia8Interface") << "Turning on BB4l hook from CMSSW Pythia8Interface";
     fPowhegHooksBB4L.reset(new PowhegHooksBB4L());
     fMultiUserHook->addHook(fPowhegHooksBB4L.get());
   }
  
   //adapted from main89.cc in pythia8 examples
   bool internalMatching = fMasterGen->settings.flag("JetMatching:merge");
   bool internalMerging = !(fMasterGen->settings.word("Merging:Process") == "void");
  
   if (internalMatching && internalMerging) {
     throw edm::Exception(edm::errors::Configuration, "Pythia8Interface")
         << " Only one jet matching/merging scheme can be used at a time. \n";
   }
  
   if (internalMatching) {
     fJetMatchingPy8InternalHook.reset(new Pythia8::JetMatchingMadgraph);
     fMultiUserHook->addHook(fJetMatchingPy8InternalHook.get());
   }
  
   if (internalMerging) {
     int scheme = (fMasterGen->settings.flag("Merging:doUMEPSTree") || fMasterGen->settings.flag("Merging:doUMEPSSubt"))
                      ? 1
                      : ((fMasterGen->settings.flag("Merging:doUNLOPSTree") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSSubt") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSLoop") ||
                          fMasterGen->settings.flag("Merging:doUNLOPSSubtNLO"))
                             ? 2
                             : 0);
     fMergingHook.reset(new Pythia8::amcnlo_unitarised_interface(scheme));
     fMultiUserHook->addHook(fMergingHook.get());
   }
   bool biasedTauDecayer = fMasterGen->settings.flag("BiasedTauDecayer:filter");
   if (biasedTauDecayer) {
     fBiasedTauDecayer.reset(new BiasedTauDecayer(&(fMasterGen->info),
                                                  &(fMasterGen->settings),
                                                  &(fMasterGen->particleData),
                                                  &(fMasterGen->rndm),
                                                  &(fMasterGen->couplings)));
     std::vector<int> handledParticles;
     handledParticles.push_back(15);
     fMasterGen->setDecayPtr(fBiasedTauDecayer.get(), handledParticles);
   }
  
   bool resonanceDecayFilter = fMasterGen->settings.flag("ResonanceDecayFilter:filter");
   if (resonanceDecayFilter) {
     fResonanceDecayFilterHook.reset(new ResonanceDecayFilterHook);
     fMultiUserHook->addHook(fResonanceDecayFilterHook.get());
   }
  
   bool PTFilter = fMasterGen->settings.flag("PTFilter:filter");
   if (PTFilter) {
     fPTFilterHook.reset(new PTFilterHook);
     fMultiUserHook->addHook(fPTFilterHook.get());
   }
  
   if (fMultiUserHook->nHooks() > 0) {
     fMasterGen->setUserHooksPtr(fMultiUserHook.get());
   }
  
   edm::LogInfo("Pythia8Interface") << "Initializing MasterGen";
   if (fvincia.get()) {
     fvincia->init();
     status = true;
   } else if (fDire.get()) {
     //fDire->initTune(*fMasterGen.get());
     fDire->weightsPtr->setup();
     fMasterGen->init();
     fDire->setup(*fMasterGen.get());
     status = true;
   } else {
     status = fMasterGen->init();
   }
  
   //clean up temp file
   if (!slhafile_.empty()) {
     std::remove(slhafile_.c_str());
   }
  
   if (pythiaPylistVerbosity > 10) {
     if (pythiaPylistVerbosity == 11 || pythiaPylistVerbosity == 13)
       fMasterGen->settings.listAll();
     if (pythiaPylistVerbosity == 12 || pythiaPylistVerbosity == 13)
       fMasterGen->particleData.listAll();
   }
  
   // init decayer
   fDecayer->settings.flag("ProcessLevel:all", false);  // trick
   fDecayer->settings.flag("ProcessLevel:resonanceDecays", true);
   edm::LogInfo("Pythia8Interface") << "Initializing Decayer";
   status1 = fDecayer->init();
  
   if (useEvtGen) {
     edm::LogInfo("Pythia8Hadronizer") << "Creating and initializing pythia8 EvtGen plugin";
     evtgenDecays.reset(new EvtGenDecays(fMasterGen.get(), evtgenDecFile, evtgenPdlFile));
     for (unsigned int i = 0; i < evtgenUserFiles.size(); i++)
       evtgenDecays->readDecayFile(evtgenUserFiles.at(i));
   }
  
   return (status && status1);
 }

◆ residualDecay()

bool Pythia8Hadronizer::residualDecay ( )

virtual

Definition at line 842 of file Pythia8Hadronizer.cc.

                                       {
   Event *pythiaEvent = &(fMasterGen->event);
  
   int NPartsBeforeDecays = pythiaEvent->size();
   int NPartsAfterDecays = event().get()->particles_size();
  
   if (NPartsAfterDecays == NPartsBeforeDecays)
     return true;
  
   bool result = true;
  
   for (int ipart = NPartsAfterDecays; ipart > NPartsBeforeDecays; ipart--) {
     HepMC::GenParticle *part = event().get()->barcode_to_particle(ipart);
  
     if (part->status() == 1 && (fDecayer->particleData).canDecay(part->pdg_id())) {
       fDecayer->event.reset();
       Particle py8part(part->pdg_id(),
                        93,
                        0,
                        0,
                        0,
                        0,
                        0,
                        0,
                        part->momentum().x(),
                        part->momentum().y(),
                        part->momentum().z(),
                        part->momentum().t(),
                        part->generated_mass());
       HepMC::GenVertex *ProdVtx = part->production_vertex();
       py8part.vProd(ProdVtx->position().x(), ProdVtx->position().y(), ProdVtx->position().z(), ProdVtx->position().t());
       py8part.tau((fDecayer->particleData).tau0(part->pdg_id()));
       fDecayer->event.append(py8part);
       int nentries = fDecayer->event.size();
       if (!fDecayer->event[nentries - 1].mayDecay())
         continue;
       fDecayer->next();
       int nentries1 = fDecayer->event.size();
       if (nentries1 <= nentries)
         continue;  //same number of particles, no decays...
  
       part->set_status(2);
  
       result = toHepMC.fill_next_event(*(fDecayer.get()), event().get(), -1, true, part);
     }
   }
  
   return result;
 }

References gen::BaseHadronizer::event(), gen::Py8InterfaceBase::fDecayer, gen::Py8InterfaceBase::fMasterGen, GenParticle::GenParticle, mps_fire::result, and gen::Py8InterfaceBase::toHepMC.

◆ statistics()

void Pythia8Hadronizer::statistics ( )

overridevirtual

Reimplemented from gen::Py8InterfaceBase.

Definition at line 646 of file Pythia8Hadronizer.cc.

                                    {
   fMasterGen->stat();
  
   if (fEmissionVetoHook.get()) {
     edm::LogPrint("Pythia8Interface") << "\n"
                                       << "Number of ISR vetoed = " << nISRveto;
     edm::LogPrint("Pythia8Interface") << "Number of FSR vetoed = " << nFSRveto;
   }
  
   double xsec = fMasterGen->info.sigmaGen();  // cross section in mb
   xsec *= 1.0e9;                              // translate to pb (CMS/Gen "convention" as of May 2009)
   double err = fMasterGen->info.sigmaErr();   // cross section err in mb
   err *= 1.0e9;                               // translate to pb (CMS/Gen "convention" as of May 2009)
   runInfo().setInternalXSec(GenRunInfoProduct::XSec(xsec, err));
 }

References runTheMatrix::err, fEmissionVetoHook, gen::Py8InterfaceBase::fMasterGen, nFSRveto, nISRveto, gen::BaseHadronizer::runInfo(), and GenRunInfoProduct::setInternalXSec().

Member Data Documentation

◆ comEnergy

double Pythia8Hadronizer::comEnergy

private

Center-of-Mass energy.

Definition at line 108 of file Pythia8Hadronizer.cc.

Referenced by initializeForInternalPartons().

◆ DJR

vector<float> Pythia8Hadronizer::DJR

private

Definition at line 166 of file Pythia8Hadronizer.cc.

Referenced by finalizeEvent(), generatePartonsAndHadronize(), and hadronize().

◆ EV1_emittedMode

int Pythia8Hadronizer::EV1_emittedMode

private

Definition at line 158 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_maxVetoCount

int Pythia8Hadronizer::EV1_maxVetoCount

private

Definition at line 155 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_MPIvetoOn

bool Pythia8Hadronizer::EV1_MPIvetoOn

private

Definition at line 160 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_nFinal

int Pythia8Hadronizer::EV1_nFinal

private

Definition at line 153 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_nFinalMode

int Pythia8Hadronizer::EV1_nFinalMode

private

Definition at line 162 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_pTdefMode

int Pythia8Hadronizer::EV1_pTdefMode

private

Definition at line 159 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_pTempMode

int Pythia8Hadronizer::EV1_pTempMode

private

Definition at line 157 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_pThardMode

int Pythia8Hadronizer::EV1_pThardMode

private

Definition at line 156 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_QEDvetoMode

int Pythia8Hadronizer::EV1_QEDvetoMode

private

Definition at line 161 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ EV1_vetoOn

bool Pythia8Hadronizer::EV1_vetoOn

private

Definition at line 154 of file Pythia8Hadronizer.cc.

Referenced by Pythia8Hadronizer().

◆ fBeam1PZ

double Pythia8Hadronizer::fBeam1PZ

private

Definition at line 116 of file Pythia8Hadronizer.cc.

◆ fBeam2PZ

double Pythia8Hadronizer::fBeam2PZ

private

Definition at line 117 of file Pythia8Hadronizer.cc.

◆ fBiasedTauDecayer

std::unique_ptr<BiasedTauDecayer> Pythia8Hadronizer::fBiasedTauDecayer

private

Definition at line 145 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fDire

std::unique_ptr<Pythia8::Dire> Pythia8Hadronizer::fDire

private

Definition at line 102 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), getGenLumiInfoHeader(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fEmissionVetoHook

std::unique_ptr<PowhegHooks> Pythia8Hadronizer::fEmissionVetoHook

private

Definition at line 137 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), initializeForExternalPartons(), initializeForInternalPartons(), and statistics().

◆ fEmissionVetoHook1

std::unique_ptr<EmissionVetoHook1> Pythia8Hadronizer::fEmissionVetoHook1

private

Definition at line 138 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fInitialState

int Pythia8Hadronizer::fInitialState

private

Definition at line 114 of file Pythia8Hadronizer.cc.

Referenced by initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fJetMatchingHook

std::unique_ptr<JetMatchingHook> Pythia8Hadronizer::fJetMatchingHook

private

Definition at line 131 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fJetMatchingPy8InternalHook

std::unique_ptr<Pythia8::JetMatchingMadgraph> Pythia8Hadronizer::fJetMatchingPy8InternalHook

private

Definition at line 132 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), initializeForExternalPartons(), and initializeForInternalPartons().

◆ fMergingHook

std::unique_ptr<Pythia8::amcnlo_unitarised_interface> Pythia8Hadronizer::fMergingHook

private

Definition at line 133 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), initializeForExternalPartons(), and initializeForInternalPartons().

◆ fMultiUserHook

std::unique_ptr<MultiUserHook> Pythia8Hadronizer::fMultiUserHook

private

Definition at line 120 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fPowhegHooksBB4L

std::unique_ptr<PowhegHooksBB4L> Pythia8Hadronizer::fPowhegHooksBB4L

private

Definition at line 142 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fPowhegResHook

std::unique_ptr<PowhegResHook> Pythia8Hadronizer::fPowhegResHook

private

Definition at line 141 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fPTFilterHook

std::unique_ptr<PTFilterHook> Pythia8Hadronizer::fPTFilterHook

private

Definition at line 151 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fResonanceDecayFilterHook

std::unique_ptr<ResonanceDecayFilterHook> Pythia8Hadronizer::fResonanceDecayFilterHook

private

Definition at line 148 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), and initializeForInternalPartons().

◆ fReweightEmpUserHook

std::unique_ptr<UserHooks> Pythia8Hadronizer::fReweightEmpUserHook

private

Definition at line 125 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fReweightPtHatRapUserHook

std::unique_ptr<UserHooks> Pythia8Hadronizer::fReweightPtHatRapUserHook

private

Definition at line 127 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fReweightRapUserHook

std::unique_ptr<UserHooks> Pythia8Hadronizer::fReweightRapUserHook

private

Definition at line 126 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fReweightUserHook

std::unique_ptr<UserHooks> Pythia8Hadronizer::fReweightUserHook

private

Definition at line 124 of file Pythia8Hadronizer.cc.

Referenced by initializeForExternalPartons(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ fvincia

std::unique_ptr<Pythia8::VinciaPlugin> Pythia8Hadronizer::fvincia

private

Definition at line 101 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), getGenLumiInfoHeader(), initializeForInternalPartons(), and Pythia8Hadronizer().

◆ lhaUP

std::unique_ptr<LHAupLesHouches> Pythia8Hadronizer::lhaUP

private

Definition at line 111 of file Pythia8Hadronizer.cc.

Referenced by hadronize(), and initializeForExternalPartons().

◆ LHEInputFileName

std::string Pythia8Hadronizer::LHEInputFileName

private

Definition at line 110 of file Pythia8Hadronizer.cc.

Referenced by hadronize(), and initializeForExternalPartons().

◆ nFSRveto

int Pythia8Hadronizer::nFSRveto

private

Definition at line 171 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), and statistics().

◆ nISRveto

int Pythia8Hadronizer::nISRveto

private

Definition at line 170 of file Pythia8Hadronizer.cc.

Referenced by generatePartonsAndHadronize(), hadronize(), and statistics().

◆ nME

int Pythia8Hadronizer::nME

private

Definition at line 167 of file Pythia8Hadronizer.cc.

Referenced by finalizeEvent(), generatePartonsAndHadronize(), and hadronize().

◆ nMEFiltered

int Pythia8Hadronizer::nMEFiltered

private

Definition at line 168 of file Pythia8Hadronizer.cc.

Referenced by finalizeEvent(), generatePartonsAndHadronize(), and hadronize().

◆ p8SharedResources

const std::vector< std::string > Pythia8Hadronizer::p8SharedResources = {edm::SharedResourceNames::kPythia8}

staticprivate

Definition at line 164 of file Pythia8Hadronizer.cc.

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ Pythia8Hadronizer()

◆ ~Pythia8Hadronizer()

Member Function Documentation

◆ classname()

◆ doSetRandomEngine()

◆ doSharedResources()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ getGenLumiInfoHeader()

◆ hadronize()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ residualDecay()

◆ statistics()

Member Data Documentation

◆ comEnergy

◆ DJR

◆ EV1_emittedMode

◆ EV1_maxVetoCount

◆ EV1_MPIvetoOn

◆ EV1_nFinal

◆ EV1_nFinalMode

◆ EV1_pTdefMode

◆ EV1_pTempMode

◆ EV1_pThardMode

◆ EV1_QEDvetoMode

◆ EV1_vetoOn

◆ fBeam1PZ

◆ fBeam2PZ

◆ fBiasedTauDecayer

◆ fDire

◆ fEmissionVetoHook

◆ fEmissionVetoHook1

◆ fInitialState

◆ fJetMatchingHook

◆ fJetMatchingPy8InternalHook

◆ fMergingHook

◆ fMultiUserHook

◆ fPowhegHooksBB4L

◆ fPowhegResHook

◆ fPTFilterHook

◆ fResonanceDecayFilterHook

◆ fReweightEmpUserHook

◆ fReweightPtHatRapUserHook

◆ fReweightRapUserHook

◆ fReweightUserHook

◆ fvincia

◆ lhaUP

◆ LHEInputFileName

◆ nFSRveto

◆ nISRveto

◆ nME

◆ nMEFiltered

◆ p8SharedResources