Inheritance diagram for gen::Py8EGun:

Public Member Functions
const char *	classname () const

bool	generatePartonsAndHadronize ()

	Py8EGun (edm::ParameterSet const &)

	~Py8EGun ()

Public Member Functions inherited from gen::Py8GunBase
void	finalizeEvent ()

edm::Event &	getEDMEvent () const

HepMC::GenEvent *	getGenEvent ()

GenEventInfoProduct *	getGenEventInfo ()

GenRunInfoProduct &	getGenRunInfo ()

bool	initializeForInternalPartons ()

	Py8GunBase (edm::ParameterSet const &ps)

void	resetEvent (HepMC::GenEvent *event)

void	resetEventInfo (GenEventInfoProduct *eventInfo)

virtual bool	residualDecay ()

virtual bool	select (HepMC::GenEvent *) const

void	setEDMEvent (edm::Event &event)

void	statistics ()

	~Py8GunBase ()

Public Member Functions inherited from gen::Py8InterfaceBase
bool	decay ()

bool	declareSpecialSettings (const std::vector< std::string > &)

bool	declareStableParticles (const std::vector< int > &)

	Py8InterfaceBase (edm::ParameterSet const &ps)

bool	readSettings (int)

	~Py8InterfaceBase ()

Private Attributes
bool	fAddAntiParticle

double	fMaxE

double	fMaxEta

double	fMinE

double	fMinEta

Additional Inherited Members
Protected Member Functions inherited from gen::Py8GunBase
std::auto_ptr< HepMC::GenEvent > &	event ()

std::auto_ptr < GenEventInfoProduct > &	eventInfo ()

GenRunInfoProduct &	runInfo ()

Protected Attributes inherited from gen::Py8GunBase
double	fMaxPhi

double	fMinPhi

std::vector< int >	fPartIDs

Protected Attributes inherited from gen::Py8InterfaceBase
std::auto_ptr< Pythia8::Pythia >	fDecayer

std::auto_ptr< Pythia8::Pythia >	fMasterGen

ParameterCollector	fParameters

unsigned int	maxEventsToPrint

bool	pythiaHepMCVerbosity

unsigned int	pythiaPylistVerbosity

HepMC::I_Pythia8	toHepMC

Detailed Description

Definition at line 10 of file Py8EGun.cc.

Constructor & Destructor Documentation

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

Definition at line 33 of file Py8EGun.cc.

References fAddAntiParticle, fMaxE, fMaxEta, fMinE, fMinEta, 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);
    fMinE       = pgun_params.getParameter<double>("MinE"); // ,  0.);
    fMaxE       = pgun_params.getParameter<double>("MaxE"); // ,  0.);
    fAddAntiParticle = pgun_params.getParameter<bool>("AddAntiParticle"); //, false) ;  
 
 }

gen::Py8EGun::~Py8EGun ( )

inline

Definition at line 15 of file Py8EGun.cc.

15 {}

Member Function Documentation

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

virtual

Implements gen::Py8InterfaceBase.

Definition at line 108 of file Py8EGun.cc.

 {
    return "Py8EGun"; 
 }

bool gen::Py8EGun::generatePartonsAndHadronize ( )

virtual

Implements gen::Py8InterfaceBase.

Definition at line 48 of file Py8EGun.cc.

References python.multivaluedict::append(), funct::cos(), eta(), gen::Py8GunBase::event(), create_public_lumi_plots::exp, fAddAntiParticle, gen::Py8InterfaceBase::fMasterGen, fMaxE, fMaxEta, gen::Py8GunBase::fMaxPhi, fMinE, fMinEta, gen::Py8GunBase::fMinPhi, gen::Py8GunBase::fPartIDs, i, phi, createTree::pp, randomEngine, funct::sin(), mathSSE::sqrt(), 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 ???
 
       // FIXME !!!
       // Ouch, it's using bare randomEngine pointer - that's NOT safe.
       // Need to hold a pointer somewhere properly !!!
       //
       double phi = (fMaxPhi-fMinPhi) * randomEngine->flat() + fMinPhi;
       double ee   = (fMaxE-fMinE) * randomEngine->flat() + fMinE;
       double eta  = (fMaxEta-fMinEta) * randomEngine->flat() + fMinEta;                                                      
       double the  = 2.*atan(exp(-eta));                                                                          
       
       double mass = (fMasterGen->particleData).mass( particleID );
 //      double mass = (pythia->particleData).m0( particleID );
 
       double 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 ); 
 
 // 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 ( (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 );
      }
       }
 
    }
    
    if ( !fMasterGen->next() ) return false;
    
    event().reset(new HepMC::GenEvent);
    toHepMC.fill_next_event( fMasterGen->event, event().get() );
       
    return true;   
   
 }