Inheritance diagram for gen::Py8PtGun:

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	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 ()(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< 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< 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 8 of file Py8PtGun.cc.

Constructor & Destructor Documentation

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

Definition at line 31 of file Py8PtGun.cc.

References fAddAntiParticle, fMaxEta, fMaxPt, fMinEta, fMinPt, 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);
    fMinPt      = pgun_params.getParameter<double>("MinPt"); // ,  0.);
    fMaxPt      = pgun_params.getParameter<double>("MaxPt"); // ,  0.);
    fAddAntiParticle = pgun_params.getParameter<bool>("AddAntiParticle"); //, false) ;  
 
 }

gen::Py8PtGun::~Py8PtGun ( )

inlineoverride

Definition at line 13 of file Py8PtGun.cc.

References classname(), and generatePartonsAndHadronize().

13 {}

Member Function Documentation

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

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 116 of file Py8PtGun.cc.

Referenced by ~Py8PtGun().

 {
    return "Py8PtGun"; 
 }

bool gen::Py8PtGun::generatePartonsAndHadronize ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 45 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, cmsBatch::log, ResonanceBuilder::mass, EgammaObjectsElectrons_cfi::particleID, createTree::pp, EnergyCorrector::pt, gen::Py8InterfaceBase::randomEngine(), funct::sin(), findQualityFiles::size, mathSSE::sqrt(), metsig::tau, and gen::Py8InterfaceBase::toHepMC.

Referenced by ~Py8PtGun().

 {
 
    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 pt   = (fMaxPt-fMinPt) * randomEngine().flat() + fMinPt;
       
       double mass = (fMasterGen->particleData).m0( particleID );
 
       double pp = pt / sin(the); // sqrt( ee*ee - mass*mass );
       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().reset(new HepMC::GenEvent);
    return toHepMC.fill_next_event( fMasterGen->event, event().get() );
   
 }