CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
gen::Py8PtGun Class Reference
Inheritance diagram for gen::Py8PtGun:
gen::Py8GunBase gen::Py8InterfaceBase gen::BaseHadronizer

Public Member Functions

const char * classname () const override
 
bool generatePartonsAndHadronize () override
 
 Py8PtGun (edm::ParameterSet const &)
 
 ~Py8PtGun () 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)
 
P8RndmEnginerandomEngine ()
 
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::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)
 

Private Attributes

bool fAddAntiParticle
 
double fMaxEta
 
double fMaxPt
 
double fMinEta
 
double fMinPt
 

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 ()
 
- 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
unsigned int ivhepmc = 2
 
std::string lheFile_
 
int randomIndex_
 

Detailed Description

Definition at line 10 of file Py8PtGun.cc.

Constructor & Destructor Documentation

◆ Py8PtGun()

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

Definition at line 29 of file Py8PtGun.cc.

References fAddAntiParticle, fMaxEta, fMaxPt, fMinEta, fMinPt, and edm::ParameterSet::getParameter().

29  : Py8GunBase(ps) {
30  // ParameterSet defpset ;
31  edm::ParameterSet pgun_params = ps.getParameter<edm::ParameterSet>("PGunParameters"); // , defpset ) ;
32  fMinEta = pgun_params.getParameter<double>("MinEta"); // ,-2.2);
33  fMaxEta = pgun_params.getParameter<double>("MaxEta"); // , 2.2);
34  fMinPt = pgun_params.getParameter<double>("MinPt"); // , 0.);
35  fMaxPt = pgun_params.getParameter<double>("MaxPt"); // , 0.);
36  fAddAntiParticle = pgun_params.getParameter<bool>("AddAntiParticle"); //, false) ;
37  }
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
bool fAddAntiParticle
Definition: Py8PtGun.cc:24
double fMinPt
Definition: Py8PtGun.cc:22
double fMinEta
Definition: Py8PtGun.cc:20
Py8GunBase(edm::ParameterSet const &ps)
Definition: Py8GunBase.cc:15
double fMaxPt
Definition: Py8PtGun.cc:23
double fMaxEta
Definition: Py8PtGun.cc:21

◆ ~Py8PtGun()

gen::Py8PtGun::~Py8PtGun ( )
inlineoverride

Definition at line 13 of file Py8PtGun.cc.

13 {}

Member Function Documentation

◆ classname()

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

Implements gen::Py8InterfaceBase.

Definition at line 114 of file Py8PtGun.cc.

114 { return "Py8PtGun"; }

◆ generatePartonsAndHadronize()

bool gen::Py8PtGun::generatePartonsAndHadronize ( )
overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 39 of file Py8PtGun.cc.

References funct::abs(), mps_setup::append, funct::cos(), PVValHelper::eta, gen::BaseHadronizer::event(), gen::Py8GunBase::evtGenDecay(), JetChargeProducer_cfi::exp, fAddAntiParticle, gen::P8RndmEngine::flat(), gen::Py8InterfaceBase::fMasterGen, fMaxEta, gen::Py8GunBase::fMaxPhi, fMaxPt, fMinEta, gen::Py8GunBase::fMinPhi, fMinPt, gen::Py8GunBase::fPartIDs, mps_fire::i, CrabHelper::log, EgHLTOffHistBins_cfi::mass, EgammaObjectsElectrons_cfi::particleID, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, gen::Py8InterfaceBase::randomEngine(), funct::sin(), mathSSE::sqrt(), metsig::tau, and gen::Py8InterfaceBase::toHepMC.

39  {
40  fMasterGen->event.reset();
41 
42  int NTotParticles = fPartIDs.size();
43  if (fAddAntiParticle)
44  NTotParticles *= 2;
45 
46  // energy below is dummy, it is not used
47  (fMasterGen->event).append(990, -11, 0, 0, 2, 1 + NTotParticles, 0, 0, 0., 0., 0., 15000., 15000.);
48 
49  for (size_t i = 0; i < fPartIDs.size(); i++) {
50  int particleID = fPartIDs[i]; // this is PDG - need to convert to Py8 ???
51 
52  double phi = (fMaxPhi - fMinPhi) * randomEngine().flat() + fMinPhi;
53  double eta = (fMaxEta - fMinEta) * randomEngine().flat() + fMinEta;
54  double the = 2. * atan(exp(-eta));
55 
56  double pt = (fMaxPt - fMinPt) * randomEngine().flat() + fMinPt;
57 
58  double mass = (fMasterGen->particleData).m0(particleID);
59 
60  double pp = pt / sin(the); // sqrt( ee*ee - mass*mass );
61  double ee = sqrt(pp * pp + mass * mass);
62 
63  double px = pt * cos(phi);
64  double py = pt * sin(phi);
65  double pz = pp * cos(the);
66 
67  if (!((fMasterGen->particleData).isParticle(particleID))) {
69  }
70  if (1 <= std::abs(particleID) && std::abs(particleID) <= 6) // quarks
71  (fMasterGen->event).append(particleID, 23, 1, 0, 0, 0, 101, 0, px, py, pz, ee, mass);
72  else if (std::abs(particleID) == 21) // gluons
73  (fMasterGen->event).append(21, 23, 1, 0, 0, 0, 101, 102, px, py, pz, ee, mass);
74  // other
75  else {
76  (fMasterGen->event).append(particleID, 1, 1, 0, 0, 0, 0, 0, px, py, pz, ee, mass);
77  int eventSize = (fMasterGen->event).size() - 1;
78  // -log(flat) = exponential distribution
79  double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
80  (fMasterGen->event)[eventSize].tau(tauTmp);
81  }
82 
83  // Here also need to add anti-particle (if any)
84  // otherwise just add a 2nd particle of the same type
85  // (for example, gamma)
86  //
87  if (fAddAntiParticle) {
88  if (1 <= std::abs(particleID) && std::abs(particleID) <= 6) { // quarks
89  (fMasterGen->event).append(-particleID, 23, 1, 0, 0, 0, 0, 101, -px, -py, -pz, ee, mass);
90  } else if (std::abs(particleID) == 21) { // gluons
91  (fMasterGen->event).append(21, 23, 1, 0, 0, 0, 102, 101, -px, -py, -pz, ee, mass);
92  } else {
93  if ((fMasterGen->particleData).isParticle(-particleID)) {
94  (fMasterGen->event).append(-particleID, 1, 1, 0, 0, 0, 0, 0, -px, -py, -pz, ee, mass);
95  } else {
96  (fMasterGen->event).append(particleID, 1, 1, 0, 0, 0, 0, 0, -px, -py, -pz, ee, mass);
97  }
98  int eventSize = (fMasterGen->event).size() - 1;
99  // -log(flat) = exponential distribution
100  double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
101  (fMasterGen->event)[eventSize].tau(tauTmp);
102  }
103  }
104  }
105 
106  if (!fMasterGen->next())
107  return false;
108  evtGenDecay();
109 
110  event() = std::make_unique<HepMC::GenEvent>();
111  return toHepMC.fill_next_event(fMasterGen->event, event().get());
112  }
size
Write out results.
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
double fMinPhi
Definition: Py8GunBase.h:58
double flat() override
Definition: P8RndmEngine.cc:7
bool fAddAntiParticle
Definition: Py8PtGun.cc:24
void evtGenDecay()
Definition: Py8GunBase.cc:144
T sqrt(T t)
Definition: SSEVec.h:23
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::unique_ptr< HepMC::GenEvent > & event()
std::vector< int > fPartIDs
Definition: Py8GunBase.h:57
P8RndmEngine & randomEngine()
double fMinPt
Definition: Py8PtGun.cc:22
double fMinEta
Definition: Py8PtGun.cc:20
double fMaxPhi
Definition: Py8GunBase.h:59
HepMC::Pythia8ToHepMC toHepMC
std::unique_ptr< Pythia8::Pythia > fMasterGen
double fMaxPt
Definition: Py8PtGun.cc:23
double fMaxEta
Definition: Py8PtGun.cc:21
Definition: event.py:1

Member Data Documentation

◆ fAddAntiParticle

bool gen::Py8PtGun::fAddAntiParticle
private

Definition at line 24 of file Py8PtGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8PtGun().

◆ fMaxEta

double gen::Py8PtGun::fMaxEta
private

Definition at line 21 of file Py8PtGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8PtGun().

◆ fMaxPt

double gen::Py8PtGun::fMaxPt
private

Definition at line 23 of file Py8PtGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8PtGun().

◆ fMinEta

double gen::Py8PtGun::fMinEta
private

Definition at line 20 of file Py8PtGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8PtGun().

◆ fMinPt

double gen::Py8PtGun::fMinPt
private

Definition at line 22 of file Py8PtGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8PtGun().