CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Attributes | Private Attributes
gen::Py8HMC3InterfaceBase Class Referenceabstract

#include <Py8HMC3InterfaceBase.h>

Inheritance diagram for gen::Py8HMC3InterfaceBase:
gen::BaseHadronizer Pythia8HepMC3Hadronizer

Public Member Functions

virtual const char * classname () const =0
 
bool decay ()
 
bool declareSpecialSettings (const std::vector< std::string > &)
 
bool declareStableParticles (const std::vector< int > &)
 
virtual void finalizeEvent ()=0
 
virtual bool generatePartonsAndHadronize ()=0
 
virtual bool initializeForInternalPartons ()=0
 
void makeTmpSLHA (const std::string &)
 
void p8SetRandomEngine (CLHEP::HepRandomEngine *v)
 
 Py8HMC3InterfaceBase (edm::ParameterSet const &ps)
 
P8RndmEnginerandomEngine ()
 
bool readSettings (int)
 
virtual void statistics ()
 
 ~Py8HMC3InterfaceBase () override=default
 
- 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::EventgetEDMEvent () const
 
std::unique_ptr< HepMC::GenEventgetGenEvent ()
 
std::unique_ptr< HepMC3::GenEvent > getGenEvent3 ()
 
std::unique_ptr< GenEventInfoProductgetGenEventInfo ()
 
std::unique_ptr< GenEventInfoProduct3getGenEventInfo3 ()
 
virtual std::unique_ptr< GenLumiInfoHeadergetGenLumiInfoHeader () const
 
GenRunInfoProductgetGenRunInfo ()
 
std::unique_ptr< lhef::LHEEventgetLHEEvent ()
 
const std::shared_ptr< lhef::LHERunInfo > & getLHERunInfo () const
 
unsigned int getVHepMC ()
 
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 resetEvent3 (std::unique_ptr< HepMC3::GenEvent > event3)
 
void resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)
 
void resetEventInfo3 (std::unique_ptr< GenEventInfoProduct3 > 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)
 

Protected Attributes

std::shared_ptr< Pythia8::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_
 
HepMC3::Pythia8ToHepMC3 toHepMC
 
bool useEvtGen
 
- Protected Attributes inherited from gen::BaseHadronizer
unsigned int ivhepmc = 2
 
std::string lheFile_
 
int randomIndex_
 

Private Attributes

std::shared_ptr< P8RndmEnginep8RndmEngine_
 

Additional Inherited Members

- Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > & event ()
 
std::unique_ptr< HepMC3::GenEvent > & event3 ()
 
std::unique_ptr< GenEventInfoProduct > & eventInfo ()
 
std::unique_ptr< GenEventInfoProduct3 > & eventInfo3 ()
 
lhef::LHEEventlheEvent ()
 
lhef::LHERunInfolheRunInfo ()
 
GenRunInfoProductrunInfo ()
 

Detailed Description

Definition at line 24 of file Py8HMC3InterfaceBase.h.

Constructor & Destructor Documentation

◆ Py8HMC3InterfaceBase()

gen::Py8HMC3InterfaceBase::Py8HMC3InterfaceBase ( edm::ParameterSet const &  ps)

Definition at line 15 of file Py8HMC3InterfaceBase.cc.

References web.browse_db::env, evtgenDecFile, evtgenPdlFile, evtgenUserFiles, Exception, edm::ParameterSet::exists(), ztee::fd, fParameters, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, maxEventsToPrint, p8RndmEngine_, pythiaHepMCVerbosity, pythiaHepMCVerbosityParticles, pythiaPylistVerbosity, AlCaHLTBitMon_QueryRunRegistry::string, useEvtGen, and writeEcalDQMStatus::write.

16  : BaseHadronizer(ps), useEvtGen(false), evtgenDecays(nullptr) {
17  fParameters = ps;
18 
19  pythiaPylistVerbosity = ps.getUntrackedParameter<int>("pythiaPylistVerbosity", 0);
20  pythiaHepMCVerbosity = ps.getUntrackedParameter<bool>("pythiaHepMCVerbosity", false);
21  pythiaHepMCVerbosityParticles = ps.getUntrackedParameter<bool>("pythiaHepMCVerbosityParticles", false);
22  maxEventsToPrint = ps.getUntrackedParameter<int>("maxEventsToPrint", 0);
23  p8RndmEngine_ = std::make_shared<P8RndmEngine>();
24 
25  //if (pythiaHepMCVerbosityParticles)
26  // ascii_io = new HepMC::IO_AsciiParticles("cout", std::ios::out);
27 
28  if (ps.exists("useEvtGenPlugin")) {
29  useEvtGen = true;
30  auto env = std::getenv("EVTGENDATA");
31  if (not env) {
32  throw cms::Exception("EvtGenMissingEnv") << "The environment variable EVTGENDATA must be defined";
33  }
34  string evtgenpath(env);
35  evtgenDecFile = evtgenpath + string("/DECAY_2010.DEC");
36  evtgenPdlFile = evtgenpath + string("/evt.pdl");
37 
38  if (ps.exists("evtgenDecFile")) {
39  edm::FileInPath decay_table(ps.getParameter<std::string>("evtgenDecFile"));
40  evtgenDecFile = decay_table.fullPath();
41  }
42 
43  if (ps.exists("evtgenPdlFile")) {
44  edm::FileInPath pdt(ps.getParameter<std::string>("evtgenPdlFile"));
45  evtgenPdlFile = pdt.fullPath();
46  }
47 
48  if (ps.exists("evtgenUserFile")) {
49  std::vector<std::string> user_decays = ps.getParameter<std::vector<std::string> >("evtgenUserFile");
50  for (unsigned int i = 0; i < user_decays.size(); i++) {
51  edm::FileInPath user_decay(user_decays.at(i));
52  evtgenUserFiles.push_back(user_decay.fullPath());
53  }
54  }
55 
56  if (ps.exists("evtgenUserFileEmbedded")) {
57  std::vector<std::string> user_decay_lines =
58  ps.getParameter<std::vector<std::string> >("evtgenUserFileEmbedded");
59  char tempslhaname[] = "pythia8evtgenXXXXXX";
60  int fd = mkstemp(tempslhaname);
61 
62  for (unsigned int i = 0; i < user_decay_lines.size(); i++) {
63  user_decay_lines.at(i) += "\n";
64  write(fd, user_decay_lines.at(i).c_str(), user_decay_lines.at(i).size());
65  }
66  close(fd);
67  evtgenUserFiles.push_back(std::string(tempslhaname));
68  }
69  }
70  }
BaseHadronizer(edm::ParameterSet const &ps)
std::shared_ptr< P8RndmEngine > p8RndmEngine_
std::shared_ptr< Pythia8::EvtGenDecays > evtgenDecays
std::vector< std::string > evtgenUserFiles
fd
Definition: ztee.py:136

◆ ~Py8HMC3InterfaceBase()

gen::Py8HMC3InterfaceBase::~Py8HMC3InterfaceBase ( )
overridedefault

Member Function Documentation

◆ classname()

virtual const char* gen::Py8HMC3InterfaceBase::classname ( ) const
pure virtual

Implemented in Pythia8HepMC3Hadronizer.

◆ decay()

bool gen::Py8HMC3InterfaceBase::decay ( )
inline

Definition at line 30 of file Py8HMC3InterfaceBase.h.

30 { return true; } // NOT used - let's call it "design imperfection"

◆ declareSpecialSettings()

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

Definition at line 206 of file Py8HMC3InterfaceBase.cc.

References fDecayer, spr::find(), fMasterGen, and AlCaHLTBitMon_QueryRunRegistry::string.

206  {
207  for (unsigned int iss = 0; iss < settings.size(); iss++) {
208  if (settings[iss].find("QED-brem-off") != std::string::npos) {
209  fMasterGen->readString("TimeShower:QEDshowerByL=off");
210  } else {
211  size_t fnd1 = settings[iss].find("Pythia8:");
212  if (fnd1 != std::string::npos) {
213  std::string value = settings[iss].substr(fnd1 + 8);
214  fDecayer->readString(value);
215  }
216  }
217  }
218  return true;
219  }
std::unique_ptr< Pythia8::Pythia > fDecayer
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
Definition: value.py:1
std::unique_ptr< Pythia8::Pythia > fMasterGen

◆ declareStableParticles()

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

Definition at line 185 of file Py8HMC3InterfaceBase.cc.

References fMasterGen, mps_fire::i, and CosmicGenFilterHelix_cfi::pdgIds.

185  {
186  for (size_t i = 0; i < pdgIds.size(); i++) {
187  int PyID = pdgIds[i];
188  std::ostringstream pyCard;
189  pyCard << PyID << ":mayDecay=false";
190 
191  if (fMasterGen->particleData.isParticle(PyID)) {
192  fMasterGen->readString(pyCard.str());
193  } else {
194  edm::LogWarning("DataNotUnderstood") << "Pythia8 does not "
195  << "recognize particle id = " << PyID << std::endl;
196  }
197  // alternative:
198  // set the 2nd input argument warn=false
199  // - this way Py8 will NOT print warnings about unknown particle code(s)
200  // fMasterPtr->readString( pyCard.str(), false )
201  }
202 
203  return true;
204  }
std::unique_ptr< Pythia8::Pythia > fMasterGen
Log< level::Warning, false > LogWarning

◆ finalizeEvent()

virtual void gen::Py8HMC3InterfaceBase::finalizeEvent ( )
pure virtual

Implemented in Pythia8HepMC3Hadronizer.

◆ generatePartonsAndHadronize()

virtual bool gen::Py8HMC3InterfaceBase::generatePartonsAndHadronize ( )
pure virtual

Implemented in Pythia8HepMC3Hadronizer.

◆ initializeForInternalPartons()

virtual bool gen::Py8HMC3InterfaceBase::initializeForInternalPartons ( )
pure virtual

Implemented in Pythia8HepMC3Hadronizer.

◆ makeTmpSLHA()

void gen::Py8HMC3InterfaceBase::makeTmpSLHA ( const std::string &  slhatable)

Definition at line 173 of file Py8HMC3InterfaceBase.cc.

References ztee::fd, fMasterGen, slhafile_, and writeEcalDQMStatus::write.

Referenced by readSettings().

173  {
174  char tempslhaname[] = "pythia8SLHAtableXXXXXX";
175  int fd = mkstemp(tempslhaname);
176  write(fd, slhatable.c_str(), slhatable.size());
177  close(fd);
178 
179  slhafile_ = tempslhaname;
180 
181  fMasterGen->settings.mode("SLHA:readFrom", 2);
182  fMasterGen->settings.word("SLHA:file", slhafile_);
183  }
std::unique_ptr< Pythia8::Pythia > fMasterGen
fd
Definition: ztee.py:136

◆ p8SetRandomEngine()

void gen::Py8HMC3InterfaceBase::p8SetRandomEngine ( CLHEP::HepRandomEngine *  v)
inline

Definition at line 40 of file Py8HMC3InterfaceBase.h.

References p8RndmEngine_, and gen::v.

40 { p8RndmEngine_->setRandomEngine(v); }
double v[5][pyjets_maxn]
std::shared_ptr< P8RndmEngine > p8RndmEngine_

◆ randomEngine()

P8RndmEngine& gen::Py8HMC3InterfaceBase::randomEngine ( )
inline

Definition at line 41 of file Py8HMC3InterfaceBase.h.

References p8RndmEngine_.

41 { return *p8RndmEngine_; }
std::shared_ptr< P8RndmEngine > p8RndmEngine_

◆ readSettings()

bool gen::Py8HMC3InterfaceBase::readSettings ( int  )

Definition at line 72 of file Py8HMC3InterfaceBase.cc.

References gen::ParameterCollector::begin(), gen::ParameterCollector::end(), edm::ParameterSet::exists(), validate-o2o-wbm::f1, fDecayer, fMasterGen, fParameters, get, edm::ParameterSet::getParameter(), mps_splice::line, makeTmpSLHA(), p8RndmEngine_, gen::BaseHadronizer::randomIndex(), SQLiteCheck_cfg::reader, slhafile_, and AlCaHLTBitMon_QueryRunRegistry::string.

72  {
73  if (!fMasterGen.get())
74  fMasterGen = std::make_unique<Pythia>();
75  fDecayer = std::make_unique<Pythia>();
76 
77  //add settings for resonance decay filter
78  fMasterGen->settings.addFlag("BiasedTauDecayer:filter", false);
79  fMasterGen->settings.addFlag("BiasedTauDecayer:eDecays", true);
80  fMasterGen->settings.addFlag("BiasedTauDecayer:muDecays", true);
81 
82  //add settings for resonance decay filter
83  fMasterGen->settings.addFlag("ResonanceDecayFilter:filter", false);
84  fMasterGen->settings.addFlag("ResonanceDecayFilter:exclusive", false);
85  fMasterGen->settings.addFlag("ResonanceDecayFilter:eMuAsEquivalent", false);
86  fMasterGen->settings.addFlag("ResonanceDecayFilter:eMuTauAsEquivalent", false);
87  fMasterGen->settings.addFlag("ResonanceDecayFilter:allNuAsEquivalent", false);
88  fMasterGen->settings.addFlag("ResonanceDecayFilter:udscAsEquivalent", false);
89  fMasterGen->settings.addFlag("ResonanceDecayFilter:udscbAsEquivalent", false);
90  fMasterGen->settings.addFlag("ResonanceDecayFilter:wzAsEquivalent", false);
91  fMasterGen->settings.addMVec("ResonanceDecayFilter:mothers", std::vector<int>(), false, false, 0, 0);
92  fMasterGen->settings.addMVec("ResonanceDecayFilter:daughters", std::vector<int>(), false, false, 0, 0);
93 
94  //add settings for PT filter
95  fMasterGen->settings.addFlag("PTFilter:filter", false);
96  fMasterGen->settings.addMode("PTFilter:quarkToFilter", 5, true, true, 3, 6);
97  fMasterGen->settings.addParm("PTFilter:scaleToFilter", 0.4, true, true, 0.0, 10.);
98  fMasterGen->settings.addParm("PTFilter:quarkRapidity", 10.0, true, true, 0.0, 10.);
99  fMasterGen->settings.addParm("PTFilter:quarkPt", -.1, true, true, -.1, 100.);
100 
101  //add settings for powheg resonance scale calculation
102  fMasterGen->settings.addFlag("POWHEGres:calcScales", false);
103  fMasterGen->settings.addFlag("POWHEG:bb4l", false);
104  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:onlyDistance1", false);
105  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:veto", false);
106  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:dryRun", false);
107  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:vetoAtPL", false);
108  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:vetoQED", false);
109  fMasterGen->settings.addFlag("POWHEG:bb4l:PartonLevel:veto", false);
110  fMasterGen->settings.addFlag("POWHEG:bb4l:PartonLevel:excludeFSRConflicting", false);
111  fMasterGen->settings.addFlag("POWHEG:bb4l:DEBUG", false);
112  fMasterGen->settings.addFlag("POWHEG:bb4l:ScaleResonance:veto", false);
113  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:vetoDipoleFrame", false);
114  fMasterGen->settings.addFlag("POWHEG:bb4l:FSREmission:pTpythiaVeto", false);
115  fMasterGen->settings.addParm("POWHEG:bb4l:pTminVeto", 10.0, true, true, 0.0, 10.);
116 
117  fMasterGen->setRndmEnginePtr(p8RndmEngine_);
118  fDecayer->setRndmEnginePtr(p8RndmEngine_);
119 
120  fMasterGen->readString("Next:numberShowEvent = 0");
121  fDecayer->readString("Next:numberShowEvent = 0");
122 
123  edm::ParameterSet currentParameters;
124  if (randomIndex() >= 0) {
125  std::vector<edm::ParameterSet> randomizedParameters =
126  fParameters.getParameter<std::vector<edm::ParameterSet> >("RandomizedParameters");
127  currentParameters = randomizedParameters[randomIndex()];
128  } else {
129  currentParameters = fParameters;
130  }
131 
132  ParameterCollector pCollector = currentParameters.getParameter<edm::ParameterSet>("PythiaParameters");
133 
134  for (ParameterCollector::const_iterator line = pCollector.begin(); line != pCollector.end(); ++line) {
135  if (line->find("Random:") != std::string::npos)
136  throw cms::Exception("PythiaError") << "Attempted to set random number "
137  "using Pythia commands. Please use "
138  "the RandomNumberGeneratorService."
139  << std::endl;
140 
141  if (!fMasterGen->readString(*line))
142  throw cms::Exception("PythiaError") << "Pythia 8 did not accept \"" << *line << "\"." << std::endl;
143 
144  if (line->find("ParticleDecays:") != std::string::npos) {
145  if (!fDecayer->readString(*line))
146  throw cms::Exception("PythiaError") << "Pythia 8 Decayer did not accept \"" << *line << "\"." << std::endl;
147  }
148  }
149 
150  slhafile_.clear();
151 
152  if (currentParameters.exists("SLHAFileForPythia8")) {
153  std::string slhafilenameshort = currentParameters.getParameter<std::string>("SLHAFileForPythia8");
154  edm::FileInPath f1(slhafilenameshort);
155 
156  fMasterGen->settings.mode("SLHA:readFrom", 2);
157  fMasterGen->settings.word("SLHA:file", f1.fullPath());
158  } else if (currentParameters.exists("SLHATableForPythia8")) {
159  std::string slhatable = currentParameters.getParameter<std::string>("SLHATableForPythia8");
160 
161  makeTmpSLHA(slhatable);
162  } else if (currentParameters.exists("SLHATreeForPythia8")) {
163  auto slhaReaderParams = currentParameters.getParameter<edm::ParameterSet>("SLHATreeForPythia8");
164  std::unique_ptr<SLHAReaderBase> reader =
165  SLHAReaderFactory::get()->create(slhaReaderParams.getParameter<std::string>("name"), slhaReaderParams);
166 
167  makeTmpSLHA(reader->getSLHA(currentParameters.getParameter<std::string>("ConfigDescription")));
168  }
169 
170  return true;
171  }
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
bool exists(std::string const &parameterName) const
checks if a parameter exists
std::unique_ptr< Pythia8::Pythia > fDecayer
void makeTmpSLHA(const std::string &)
int randomIndex() const
std::unique_ptr< Pythia8::Pythia > fMasterGen
std::shared_ptr< P8RndmEngine > p8RndmEngine_
#define get

◆ statistics()

void gen::Py8HMC3InterfaceBase::statistics ( )
virtual

Reimplemented in Pythia8HepMC3Hadronizer.

Definition at line 221 of file Py8HMC3InterfaceBase.cc.

References fMasterGen.

221  {
222  fMasterGen->stat();
223  return;
224  }
std::unique_ptr< Pythia8::Pythia > fMasterGen

Member Data Documentation

◆ evtgenDecays

std::shared_ptr<Pythia8::EvtGenDecays> gen::Py8HMC3InterfaceBase::evtgenDecays
protected

◆ evtgenDecFile

std::string gen::Py8HMC3InterfaceBase::evtgenDecFile
protected

◆ evtgenPdlFile

std::string gen::Py8HMC3InterfaceBase::evtgenPdlFile
protected

◆ evtgenUserFiles

std::vector<std::string> gen::Py8HMC3InterfaceBase::evtgenUserFiles
protected

◆ fDecayer

std::unique_ptr<Pythia8::Pythia> gen::Py8HMC3InterfaceBase::fDecayer
protected

◆ fMasterGen

std::unique_ptr<Pythia8::Pythia> gen::Py8HMC3InterfaceBase::fMasterGen
protected

◆ fParameters

edm::ParameterSet gen::Py8HMC3InterfaceBase::fParameters
protected

Definition at line 47 of file Py8HMC3InterfaceBase.h.

Referenced by Py8HMC3InterfaceBase(), and readSettings().

◆ maxEventsToPrint

unsigned int gen::Py8HMC3InterfaceBase::maxEventsToPrint
protected

◆ p8RndmEngine_

std::shared_ptr<P8RndmEngine> gen::Py8HMC3InterfaceBase::p8RndmEngine_
private

◆ pythiaHepMCVerbosity

bool gen::Py8HMC3InterfaceBase::pythiaHepMCVerbosity
protected

◆ pythiaHepMCVerbosityParticles

bool gen::Py8HMC3InterfaceBase::pythiaHepMCVerbosityParticles
protected

◆ pythiaPylistVerbosity

unsigned int gen::Py8HMC3InterfaceBase::pythiaPylistVerbosity
protected

◆ slhafile_

std::string gen::Py8HMC3InterfaceBase::slhafile_
protected

◆ toHepMC

HepMC3::Pythia8ToHepMC3 gen::Py8HMC3InterfaceBase::toHepMC
protected

◆ useEvtGen

bool gen::Py8HMC3InterfaceBase::useEvtGen
protected