CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
gen::Py8EGun Class Reference
Inheritance diagram for gen::Py8EGun:
gen::Py8GunBase gen::Py8InterfaceBase gen::BaseHadronizer

Public Member Functions

const char * classname () const override
 
bool generatePartonsAndHadronize () override
 
 Py8EGun (edm::ParameterSet const &)
 
 ~Py8EGun () override
 
- Public Member Functions inherited from gen::Py8GunBase
void evtGenDecay ()
 
void finalizeEvent () override
 
bool initializeForInternalPartons () override
 
 Py8GunBase (edm::ParameterSet const &ps)
 
virtual bool residualDecay ()
 
void setRandomEngine (CLHEP::HepRandomEngine *v)
 
std::vector< std::string > const & sharedResources () const
 
void statistics () override
 
 ~Py8GunBase () override
 
- Public Member Functions inherited from gen::Py8InterfaceBase
bool decay ()
 
bool declareSpecialSettings (const std::vector< std::string > &)
 
bool declareStableParticles (const std::vector< int > &)
 
void makeTmpSLHA (const std::string &)
 
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 ()
 
virtual std::unique_ptr< GenLumiInfoHeader > getGenLumiInfoHeader () const
 
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 Attributes

bool fAddAntiParticle
 
double fMaxE
 
double fMaxEta
 
double fMinE
 
double fMinEta
 

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::Py8GunBase
double fMaxPhi
 
double fMinPhi
 
std::vector< int > fPartIDs
 
- Protected Attributes inherited from gen::Py8InterfaceBase
HepMC::IO_AsciiParticles * ascii_io
 
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_
 
HepMC::Pythia8ToHepMC toHepMC
 
bool useEvtGen
 
- Protected Attributes inherited from gen::BaseHadronizer
std::string lheFile_
 
int randomIndex_
 

Detailed Description

Definition at line 11 of file Py8EGun.cc.

Constructor & Destructor Documentation

◆ Py8EGun()

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

Definition at line 30 of file Py8EGun.cc.

30  : Py8GunBase(ps) {
31  // ParameterSet defpset ;
32  edm::ParameterSet pgun_params = ps.getParameter<edm::ParameterSet>("PGunParameters"); // , defpset ) ;
33  fMinEta = pgun_params.getParameter<double>("MinEta"); // ,-2.2);
34  fMaxEta = pgun_params.getParameter<double>("MaxEta"); // , 2.2);
35  fMinE = pgun_params.getParameter<double>("MinE"); // , 0.);
36  fMaxE = pgun_params.getParameter<double>("MaxE"); // , 0.);
37  fAddAntiParticle = pgun_params.getParameter<bool>("AddAntiParticle"); //, false) ;
38  }

References fAddAntiParticle, fMaxE, fMaxEta, fMinE, fMinEta, and edm::ParameterSet::getParameter().

◆ ~Py8EGun()

gen::Py8EGun::~Py8EGun ( )
inlineoverride

Definition at line 14 of file Py8EGun.cc.

14 {}

Member Function Documentation

◆ classname()

const char * gen::Py8EGun::classname ( ) const
overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 120 of file Py8EGun.cc.

120 { return "Py8EGun"; }

◆ generatePartonsAndHadronize()

bool gen::Py8EGun::generatePartonsAndHadronize ( )
overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 40 of file Py8EGun.cc.

40  {
41  fMasterGen->event.reset();
42 
43  int colorindex = 101;
44 
45  for (size_t i = 0; i < fPartIDs.size(); i++) {
46  int particleID = fPartIDs[i]; // this is PDG - need to convert to Py8 ???
47  if ((std::abs(particleID) <= 6 || particleID == 21) && !(fAddAntiParticle)) {
48  throw cms::Exception("PythiaError") << "Attempting to generate quarks or gluons without setting "
49  "AddAntiParticle to true. This will not handle color properly."
50  << std::endl;
51  }
52 
53  double phi = (fMaxPhi - fMinPhi) * randomEngine().flat() + fMinPhi;
54  double ee = (fMaxE - fMinE) * randomEngine().flat() + fMinE;
55  double eta = (fMaxEta - fMinEta) * randomEngine().flat() + fMinEta;
56  double the = 2. * atan(exp(-eta));
57 
58  double mass = (fMasterGen->particleData).m0(particleID);
59 
60  double pp = sqrt(ee * ee - mass * mass);
61  double px = pp * sin(the) * cos(phi);
62  double py = pp * sin(the) * sin(phi);
63  double pz = pp * cos(the);
64 
65  if (!((fMasterGen->particleData).isParticle(particleID))) {
66  particleID = std::fabs(particleID);
67  }
68  if (1 <= std::abs(particleID) && std::abs(particleID) <= 6) { // quarks
69  (fMasterGen->event).append(particleID, 23, colorindex, 0, px, py, pz, ee, mass);
70  if (!fAddAntiParticle)
71  colorindex += 1;
72  } else if (std::abs(particleID) == 21) { // gluons
73  (fMasterGen->event).append(21, 23, colorindex, colorindex + 1, px, py, pz, ee, mass);
74  if (!fAddAntiParticle) {
75  colorindex += 2;
76  }
77  }
78  // other
79  else {
80  (fMasterGen->event).append(particleID, 1, 0, 0, px, py, pz, ee, mass);
81  int eventSize = (fMasterGen->event).size() - 1;
82  // -log(flat) = exponential distribution
83  double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
84  (fMasterGen->event)[eventSize].tau(tauTmp);
85  }
86 
87  // Here also need to add anti-particle (if any)
88  // otherwise just add a 2nd particle of the same type
89  // (for example, gamma)
90  //
91  if (fAddAntiParticle) {
92  if (1 <= std::abs(particleID) && std::abs(particleID) <= 6) { // quarks
93  (fMasterGen->event).append(-particleID, 23, 0, colorindex, -px, -py, -pz, ee, mass);
94  colorindex += 1;
95  } else if (std::abs(particleID) == 21) { // gluons
96  (fMasterGen->event).append(21, 23, colorindex + 1, colorindex, -px, -py, -pz, ee, mass);
97  colorindex += 2;
98  } else {
99  if ((fMasterGen->particleData).isParticle(-particleID)) {
100  (fMasterGen->event).append(-particleID, 1, 0, 0, -px, -py, -pz, ee, mass);
101  } else {
102  (fMasterGen->event).append(particleID, 1, 0, 0, -px, -py, -pz, ee, mass);
103  }
104  int eventSize = (fMasterGen->event).size() - 1;
105  // -log(flat) = exponential distribution
106  double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
107  (fMasterGen->event)[eventSize].tau(tauTmp);
108  }
109  }
110  }
111 
112  if (!fMasterGen->next())
113  return false;
114  evtGenDecay();
115 
116  event() = std::make_unique<HepMC::GenEvent>();
117  return toHepMC.fill_next_event(fMasterGen->event, event().get());
118  }

References funct::abs(), mps_setup::append, funct::cos(), PVValHelper::eta, gen::BaseHadronizer::event(), gen::Py8GunBase::evtGenDecay(), Exception, JetChargeProducer_cfi::exp, fAddAntiParticle, gen::P8RndmEngine::flat(), gen::Py8InterfaceBase::fMasterGen, fMaxE, fMaxEta, gen::Py8GunBase::fMaxPhi, fMinE, fMinEta, gen::Py8GunBase::fMinPhi, gen::Py8GunBase::fPartIDs, get, mps_fire::i, dqm-mbProfile::log, EgHLTOffHistBins_cfi::mass, EgammaObjectsElectrons_cfi::particleID, createTree::pp, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, gen::Py8InterfaceBase::randomEngine(), funct::sin(), findQualityFiles::size, mathSSE::sqrt(), metsig::tau, and gen::Py8InterfaceBase::toHepMC.

Member Data Documentation

◆ fAddAntiParticle

bool gen::Py8EGun::fAddAntiParticle
private

Definition at line 25 of file Py8EGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8EGun().

◆ fMaxE

double gen::Py8EGun::fMaxE
private

Definition at line 24 of file Py8EGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8EGun().

◆ fMaxEta

double gen::Py8EGun::fMaxEta
private

Definition at line 22 of file Py8EGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8EGun().

◆ fMinE

double gen::Py8EGun::fMinE
private

Definition at line 23 of file Py8EGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8EGun().

◆ fMinEta

double gen::Py8EGun::fMinEta
private

Definition at line 21 of file Py8EGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8EGun().

gen::Py8InterfaceBase::fMasterGen
std::unique_ptr< Pythia8::Pythia > fMasterGen
Definition: Py8InterfaceBase.h:47
gen::Py8GunBase::Py8GunBase
Py8GunBase(edm::ParameterSet const &ps)
Definition: Py8GunBase.cc:15
mps_fire.i
i
Definition: mps_fire.py:428
gen::Py8GunBase::fMinPhi
double fMinPhi
Definition: Py8GunBase.h:58
gen::Py8EGun::fMinEta
double fMinEta
Definition: Py8EGun.cc:21
metsig::tau
Definition: SignAlgoResolutions.h:49
gen::P8RndmEngine::flat
double flat() override
Definition: P8RndmEngine.cc:7
multPhiCorr_741_25nsDY_cfi.py
py
Definition: multPhiCorr_741_25nsDY_cfi.py:12
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
EgammaObjectsElectrons_cfi.particleID
particleID
Definition: EgammaObjectsElectrons_cfi.py:4
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
gen::Py8GunBase::evtGenDecay
void evtGenDecay()
Definition: Py8GunBase.cc:144
PVValHelper::eta
Definition: PVValidationHelpers.h:70
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
gen::Py8InterfaceBase::randomEngine
P8RndmEngine & randomEngine()
Definition: Py8InterfaceBase.h:44
edm::ParameterSet
Definition: ParameterSet.h:47
mps_setup.append
append
Definition: mps_setup.py:85
gen::Py8GunBase::fPartIDs
std::vector< int > fPartIDs
Definition: Py8GunBase.h:57
gen::Py8EGun::fMinE
double fMinE
Definition: Py8EGun.cc:23
get
#define get
DDAxes::phi
multPhiCorr_741_25nsDY_cfi.px
px
Definition: multPhiCorr_741_25nsDY_cfi.py:10
gen::Py8EGun::fAddAntiParticle
bool fAddAntiParticle
Definition: Py8EGun.cc:25
gen::BaseHadronizer::event
std::unique_ptr< HepMC::GenEvent > & event()
Definition: BaseHadronizer.h:86
Exception
Definition: hltDiff.cc:245
gen::Py8InterfaceBase::toHepMC
HepMC::Pythia8ToHepMC toHepMC
Definition: Py8InterfaceBase.h:49
EgHLTOffHistBins_cfi.mass
mass
Definition: EgHLTOffHistBins_cfi.py:34
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
createTree.pp
pp
Definition: createTree.py:17
gen::Py8EGun::fMaxEta
double fMaxEta
Definition: Py8EGun.cc:22
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
gen::Py8GunBase::fMaxPhi
double fMaxPhi
Definition: Py8GunBase.h:59
event
Definition: event.py:1
gen::Py8EGun::fMaxE
double fMaxE
Definition: Py8EGun.cc:24
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
Generated for CMSSW Reference Manual by doxygen 1.8.16