Inheritance diagram for gen::Py8JetGun:

Public Member Functions
const char *	classname () const override

bool	generatePartonsAndHadronize () override

	Py8JetGun (edm::ParameterSet const &)

	~Py8JetGun () 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
double	fMaxE

double	fMaxEta

double	fMaxP

double	fMinE

double	fMinEta

double	fMinP

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 Py8JetGun.cc.

Constructor & Destructor Documentation

◆ Py8JetGun()

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

Definition at line 31 of file Py8JetGun.cc.

                                                 : Py8GunBase(ps) {
     // ParameterSet defpset ;
     edm::ParameterSet pgun_params = ps.getParameter<edm::ParameterSet>("PGunParameters");  // , defpset ) ;
     fMinEta = pgun_params.getParameter<double>("MinEta");                                  // ,-2.2);
     fMaxEta = pgun_params.getParameter<double>("MaxEta");                                  // , 2.2);
     fMinP = pgun_params.getParameter<double>("MinP");                                      // ,  0.);
     fMaxP = pgun_params.getParameter<double>("MaxP");                                      // ,  0.);
     fMinE = pgun_params.getParameter<double>("MinE");                                      // ,  0.);
     fMaxE = pgun_params.getParameter<double>("MaxE");                                      // ,  0.);
   }

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

◆ ~Py8JetGun()

gen::Py8JetGun::~Py8JetGun ( )

inlineoverride

Definition at line 14 of file Py8JetGun.cc.

14 {}

Member Function Documentation

◆ classname()

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

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 115 of file Py8JetGun.cc.

115 { return "Py8JetGun"; }

◆ generatePartonsAndHadronize()

bool gen::Py8JetGun::generatePartonsAndHadronize ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 42 of file Py8JetGun.cc.

                                               {
     fMasterGen->event.reset();
  
     double totPx = 0.;
     double totPy = 0.;
     double totPz = 0.;
     double totE = 0.;
     double totM = 0.;
     double phi, eta, the, ee, pp;
  
     for (size_t i = 0; i < fPartIDs.size(); i++) {
       int particleID = fPartIDs[i];  // this is PDG - need to convert to Py8 ???
  
       phi = 2. * M_PI * randomEngine().flat();
       the = acos(-1. + 2. * randomEngine().flat());
  
       // from input
       //
       ee = (fMaxE - fMinE) * randomEngine().flat() + fMinE;
  
       double mass = (fMasterGen->particleData).m0(particleID);
  
       pp = sqrt(ee * ee - mass * mass);
  
       double px = pp * sin(the) * cos(phi);
       double py = pp * sin(the) * sin(phi);
       double pz = pp * cos(the);
  
       if (!((fMasterGen->particleData).isParticle(particleID))) {
         particleID = std::fabs(particleID);
       }
       (fMasterGen->event).append(particleID, 1, 0, 0, px, py, pz, ee, mass);
       int eventSize = (fMasterGen->event).size() - 1;
       // -log(flat) = exponential distribution
       double tauTmp = -(fMasterGen->event)[eventSize].tau0() * log(randomEngine().flat());
       (fMasterGen->event)[eventSize].tau(tauTmp);
  
       // values for computing total mass
       //
       totPx += px;
       totPy += py;
       totPz += pz;
       totE += ee;
     }
  
     totM = sqrt(totE * totE - (totPx * totPx + totPy * totPy + totPz * totPz));
  
     //now the boost (from input params)
     //
     pp = (fMaxP - fMinP) * randomEngine().flat() + fMinP;
     ee = sqrt(totM * totM + pp * pp);
  
     //the boost direction (from input params)
     //
     phi = (fMaxPhi - fMinPhi) * randomEngine().flat() + fMinPhi;
     eta = (fMaxEta - fMinEta) * randomEngine().flat() + fMinEta;
     the = 2. * atan(exp(-eta));
  
     double betaX = pp / ee * std::sin(the) * std::cos(phi);
     double betaY = pp / ee * std::sin(the) * std::sin(phi);
     double betaZ = pp / ee * std::cos(the);
  
     // boost all particles
     //
     (fMasterGen->event).bst(betaX, betaY, betaZ);
  
     if (!fMasterGen->next())
       return false;
  
     event() = std::make_unique<HepMC::GenEvent>();
     return toHepMC.fill_next_event(fMasterGen->event, event().get());
   }