CMS 3D CMS Logo

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

Public Member Functions

const char * classname () const override
 
bool generatePartonsAndHadronize () override
 
 Py8JetGun (edm::ParameterSet const &)
 
 ~Py8JetGun ()
 
- Public Member Functions inherited from gen::Py8GunBase
void finalizeEvent ()
 
bool initializeForInternalPartons ()
 
 Py8GunBase (edm::ParameterSet const &ps)
 
virtual bool residualDecay ()
 
void setRandomEngine (CLHEP::HepRandomEngine *v)
 
std::vector< std::string > const & sharedResources () const
 
void statistics ()
 
virtual ~Py8GunBase ()
 
- 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)
 
P8RndmEnginerandomEngine ()
 
bool readSettings (int)
 
virtual ~Py8InterfaceBase ()
 
- 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
 
HepMC::GenEvent * getGenEvent ()
 
GenEventInfoProductgetGenEventInfo ()
 
virtual GenLumiInfoHeadergetGenLumiInfoHeader () const
 
GenRunInfoProductgetGenRunInfo ()
 
const boost::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 (HepMC::GenEvent *event)
 
void resetEventInfo (GenEventInfoProduct *eventInfo)
 
virtual bool select (HepMC::GenEvent *) const
 
void setEDMEvent (edm::Event &event)
 
void setLHEEvent (lhef::LHEEvent *event)
 
void setLHERunInfo (lhef::LHERunInfo *runInfo)
 
void setRandomEngine (CLHEP::HepRandomEngine *v)
 
std::vector< std::string > const & sharedResources () const
 
virtual ~BaseHadronizer ()(false)
 

Private Attributes

double fMaxE
 
double fMaxEta
 
double fMaxP
 
double fMinE
 
double fMinEta
 
double fMinP
 

Additional Inherited Members

- Protected Member Functions inherited from gen::BaseHadronizer
std::auto_ptr< HepMC::GenEvent > & event ()
 
std::auto_ptr< GenEventInfoProduct > & eventInfo ()
 
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::auto_ptr< EvtGenDecays > evtgenDecays
 
std::string evtgenDecFile
 
std::string evtgenPdlFile
 
std::vector< std::string > evtgenUserFiles
 
std::auto_ptr< Pythia8::Pythia > fDecayer
 
std::auto_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 9 of file Py8JetGun.cc.

Constructor & Destructor Documentation

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

Definition at line 33 of file Py8JetGun.cc.

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

34  : Py8GunBase(ps)
35 {
36 
37  // ParameterSet defpset ;
38  edm::ParameterSet pgun_params =
39  ps.getParameter<edm::ParameterSet>("PGunParameters"); // , defpset ) ;
40  fMinEta = pgun_params.getParameter<double>("MinEta"); // ,-2.2);
41  fMaxEta = pgun_params.getParameter<double>("MaxEta"); // , 2.2);
42  fMinP = pgun_params.getParameter<double>("MinP"); // , 0.);
43  fMaxP = pgun_params.getParameter<double>("MaxP"); // , 0.);
44  fMinE = pgun_params.getParameter<double>("MinE"); // , 0.);
45  fMaxE = pgun_params.getParameter<double>("MaxE"); // , 0.);
46 
47 }
T getParameter(std::string const &) const
double fMinEta
Definition: Py8JetGun.cc:22
double fMinE
Definition: Py8JetGun.cc:26
double fMaxP
Definition: Py8JetGun.cc:25
double fMaxE
Definition: Py8JetGun.cc:27
Py8GunBase(edm::ParameterSet const &ps)
Definition: Py8GunBase.cc:15
double fMinP
Definition: Py8JetGun.cc:24
double fMaxEta
Definition: Py8JetGun.cc:23
gen::Py8JetGun::~Py8JetGun ( )
inline

Definition at line 14 of file Py8JetGun.cc.

References classname(), and generatePartonsAndHadronize().

14 {}

Member Function Documentation

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

Implements gen::Py8InterfaceBase.

Definition at line 128 of file Py8JetGun.cc.

Referenced by ~Py8JetGun().

129 {
130  return "Py8JetGun";
131 }
bool gen::Py8JetGun::generatePartonsAndHadronize ( )
overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 49 of file Py8JetGun.cc.

References mps_setup::append, funct::cos(), PVValHelper::eta, gen::BaseHadronizer::event(), JetChargeProducer_cfi::exp, gen::P8RndmEngine::flat(), gen::Py8InterfaceBase::fMasterGen, fMaxE, fMaxEta, fMaxP, gen::Py8GunBase::fMaxPhi, fMinE, fMinEta, fMinP, gen::Py8GunBase::fMinPhi, gen::Py8GunBase::fPartIDs, mps_fire::i, cmsBatch::log, M_PI, ResonanceBuilder::mass, EgammaObjectsElectrons_cfi::particleID, createTree::pp, gen::Py8InterfaceBase::randomEngine(), funct::sin(), findQualityFiles::size, mathSSE::sqrt(), metsig::tau, and gen::Py8InterfaceBase::toHepMC.

Referenced by ~Py8JetGun().

50 {
51 
52  fMasterGen->event.reset();
53 
54  double totPx = 0.;
55  double totPy = 0.;
56  double totPz = 0.;
57  double totE = 0.;
58  double totM = 0.;
59  double phi, eta, the, ee, pp;
60 
61  for ( size_t i=0; i<fPartIDs.size(); i++ )
62  {
63 
64  int particleID = fPartIDs[i]; // this is PDG - need to convert to Py8 ???
65 
66  phi = 2. * M_PI * randomEngine().flat() ;
67  the = acos( -1. + 2.*randomEngine().flat() );
68 
69  // from input
70  //
71  ee = (fMaxE-fMinE)*randomEngine().flat() + fMinE;
72 
73  double mass = (fMasterGen->particleData).m0( particleID );
74 
75  pp = sqrt( ee*ee - mass*mass );
76 
77  double px = pp * sin(the) * cos(phi);
78  double py = pp * sin(the) * sin(phi);
79  double pz = pp * cos(the);
80 
81  if ( !((fMasterGen->particleData).isParticle( particleID )) )
82  {
83  particleID = std::fabs(particleID) ;
84  }
85  (fMasterGen->event).append( particleID, 1, 0, 0, px, py, pz, ee, mass );
86  int eventSize = (fMasterGen->event).size()-1;
87  // -log(flat) = exponential distribution
88  double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
89  (fMasterGen->event)[eventSize].tau( tauTmp );
90 
91  // values for computing total mass
92  //
93  totPx += px;
94  totPy += py;
95  totPz += pz;
96  totE += ee;
97 
98  }
99 
100  totM = sqrt( totE*totE - (totPx*totPx+totPy*totPy+totPz*totPz) );
101 
102  //now the boost (from input params)
103  //
104  pp = (fMaxP-fMinP)*randomEngine().flat() + fMinP;
105  ee = sqrt( totM*totM + pp*pp );
106 
107  //the boost direction (from input params)
108  //
109  phi = (fMaxPhi-fMinPhi)*randomEngine().flat() + fMinPhi;
110  eta = (fMaxEta-fMinEta)*randomEngine().flat() + fMinEta;
111  the = 2.*atan(exp(-eta));
112 
113  double betaX = pp/ee * std::sin(the) * std::cos(phi);
114  double betaY = pp/ee * std::sin(the) * std::sin(phi);
115  double betaZ = pp/ee * std::cos(the);
116 
117  // boost all particles
118  //
119  (fMasterGen->event).bst( betaX, betaY, betaZ );
120 
121  if ( !fMasterGen->next() ) return false;
122 
123  event().reset(new HepMC::GenEvent);
124  return toHepMC.fill_next_event( fMasterGen->event, event().get() );
125 
126 }
size
Write out results.
std::auto_ptr< Pythia8::Pythia > fMasterGen
double fMinEta
Definition: Py8JetGun.cc:22
double fMinE
Definition: Py8JetGun.cc:26
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
double fMinPhi
Definition: Py8GunBase.h:59
double fMaxP
Definition: Py8JetGun.cc:25
virtual double flat() override
Definition: P8RndmEngine.cc:7
double fMaxE
Definition: Py8JetGun.cc:27
std::auto_ptr< HepMC::GenEvent > & event()
T sqrt(T t)
Definition: SSEVec.h:18
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
std::vector< int > fPartIDs
Definition: Py8GunBase.h:58
#define M_PI
P8RndmEngine & randomEngine()
double fMaxPhi
Definition: Py8GunBase.h:60
HepMC::Pythia8ToHepMC toHepMC
double fMinP
Definition: Py8JetGun.cc:24
double fMaxEta
Definition: Py8JetGun.cc:23
Definition: event.py:1

Member Data Documentation

double gen::Py8JetGun::fMaxE
private

Definition at line 27 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().

double gen::Py8JetGun::fMaxEta
private

Definition at line 23 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().

double gen::Py8JetGun::fMaxP
private

Definition at line 25 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().

double gen::Py8JetGun::fMinE
private

Definition at line 26 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().

double gen::Py8JetGun::fMinEta
private

Definition at line 22 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().

double gen::Py8JetGun::fMinP
private

Definition at line 24 of file Py8JetGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8JetGun().