Inheritance diagram for gen::Py8PtotGun:

Public Member Functions
const char *	classname () const override

bool	generatePartonsAndHadronize () override

	Py8PtotGun (edm::ParameterSet const &)

	~Py8PtotGun () 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	fMaxEta

double	fMaxPtot

double	fMinEta

double	fMinPtot

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 10 of file Py8PtotGun.cc.

Constructor & Destructor Documentation

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

Definition at line 29 of file Py8PtotGun.cc.

References fAddAntiParticle, fMaxEta, fMaxPtot, fMinEta, fMinPtot, and edm::ParameterSet::getParameter().

                                                   : 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);
     fMinPtot = pgun_params.getParameter<double>("MinPtot");                                // ,  0.);
     fMaxPtot = pgun_params.getParameter<double>("MaxPtot");                                // ,  0.);
     fAddAntiParticle = pgun_params.getParameter<bool>("AddAntiParticle");                  //, false) ;
   }

gen::Py8PtotGun::~Py8PtotGun ( )

inlineoverride

Definition at line 13 of file Py8PtotGun.cc.

13 {}

Member Function Documentation

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

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 106 of file Py8PtotGun.cc.

106 { return "Py8PtotGun"; }

bool gen::Py8PtotGun::generatePartonsAndHadronize ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 39 of file Py8PtotGun.cc.

References funct::abs(), bitset_utilities::append(), funct::cos(), PVValHelper::eta, gen::BaseHadronizer::event(), gen::Py8GunBase::evtGenDecay(), funct::exp(), fAddAntiParticle, gen::P8RndmEngine::flat(), gen::Py8InterfaceBase::fMasterGen, fMaxEta, gen::Py8GunBase::fMaxPhi, fMaxPtot, fMinEta, gen::Py8GunBase::fMinPhi, fMinPtot, gen::Py8GunBase::fPartIDs, mps_fire::i, log, ResonanceBuilder::mass, createTree::pp, DiDispStaMuonMonitor_cfi::pt, gen::Py8InterfaceBase::randomEngine(), funct::sin(), findQualityFiles::size, mathSSE::sqrt(), metsig::tau, and gen::Py8InterfaceBase::toHepMC.

                                                {
     fMasterGen->event.reset();
 
     for (size_t i = 0; i < fPartIDs.size(); i++) {
       int particleID = fPartIDs[i];  // this is PDG - need to convert to Py8 ???
 
       double phi = (fMaxPhi - fMinPhi) * randomEngine().flat() + fMinPhi;
       double eta = (fMaxEta - fMinEta) * randomEngine().flat() + fMinEta;
       double the = 2. * atan(exp(-eta));
 
       double pp = (fMaxPtot - fMinPtot) * randomEngine().flat() + fMinPtot;
 
       double mass = (fMasterGen->particleData).m0(particleID);
 
       double pt = pp * sin(the);
       double ee = sqrt(pp * pp + mass * mass);
       double px = pt * cos(phi);
       double py = pt * sin(phi);
       double pz = pp * cos(the);
 
       if (!((fMasterGen->particleData).isParticle(particleID))) {
         particleID = std::abs(particleID);
       }
       if (1 <= std::abs(particleID) && std::abs(particleID) <= 6)  // quarks
         (fMasterGen->event).append(particleID, 23, 101, 0, px, py, pz, ee, mass);
       else if (std::abs(particleID) == 21)  // gluons
         (fMasterGen->event).append(21, 23, 101, 102, px, py, pz, ee, mass);
       // other
       else {
         (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);
       }
 
       // Here also need to add anti-particle (if any)
       // otherwise just add a 2nd particle of the same type
       // (for example, gamma)
       //
       if (fAddAntiParticle) {
         if (1 <= std::abs(particleID) && std::abs(particleID) <= 6) {  // quarks
           (fMasterGen->event).append(-particleID, 23, 0, 101, -px, -py, -pz, ee, mass);
         } else if (std::abs(particleID) == 21) {  // gluons
           (fMasterGen->event).append(21, 23, 102, 101, -px, -py, -pz, ee, mass);
         } else {
           if ((fMasterGen->particleData).isParticle(-particleID)) {
             (fMasterGen->event).append(-particleID, 1, 0, 0, -px, -py, -pz, ee, mass);
           } else {
             (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);
         }
       }
     }
 
     if (!fMasterGen->next())
       return false;
     evtGenDecay();
 
     event() = std::make_unique<HepMC::GenEvent>();
     return toHepMC.fill_next_event(fMasterGen->event, event().get());
   }