#include <Pythia6PartonEGun.h>

Inheritance diagram for gen::Pythia6PartonEGun:

Public Member Functions
	Pythia6PartonEGun (const edm::ParameterSet &)

virtual	~Pythia6PartonEGun ()

Public Member Functions inherited from gen::Pythia6PartonGun
	Pythia6PartonGun (const edm::ParameterSet &)

virtual	~Pythia6PartonGun ()

Public Member Functions inherited from gen::Pythia6Gun
void	beginJob ()

void	beginRun (edm::Run &, edm::EventSetup const &)

void	endJob ()

void	endRun (edm::Run &, edm::EventSetup const &)

void	produce (edm::Event &, const edm::EventSetup &)

	Pythia6Gun (const edm::ParameterSet &)

virtual	~Pythia6Gun ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Protected Member Functions
void	generateEvent ()

Protected Member Functions inherited from gen::Pythia6PartonGun
void	joinPartons (double qmax)

Protected Member Functions inherited from gen::Pythia6Gun
HepMC::GenParticle *	addAntiParticle (int &, int &, double &, double &, double &)

void	attachPy6DecaysToGenEvent ()

void	loadEvent (edm::Event &)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Private Attributes
double	fMaxE

double	fMaxEta

double	fMinE

double	fMinEta

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from gen::Pythia6PartonGun
int	fPartonID

Protected Attributes inherited from gen::Pythia6Gun
HepMC::GenEvent *	fEvt

bool	fHepMCVerbosity

int	fMaxEventsToPrint

double	fMaxPhi

double	fMinPhi

std::vector< int >	fPartIDs

Pythia6Service *	fPy6Service

int	fPylistVerbosity

Detailed Description

Definition at line 8 of file Pythia6PartonEGun.h.

Constructor & Destructor Documentation

Pythia6PartonEGun::Pythia6PartonEGun ( const edm::ParameterSet & pset )

Definition at line 18 of file Pythia6PartonEGun.cc.

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

                                                                :
    Pythia6PartonGun(pset)
 {
    
    // ParameterSet defpset ;
    ParameterSet pgun_params = 
       pset.getParameter<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"); // ,  20.);
    fMaxE       = pgun_params.getParameter<double>("MaxE"); // , 420.);
 
 }

Pythia6PartonEGun::~Pythia6PartonEGun ( )

virtual

Definition at line 32 of file Pythia6PartonEGun.cc.

33 {

34 }

Member Function Documentation

void Pythia6PartonEGun::generateEvent ( )

protectedvirtual

Implements gen::Pythia6Gun.

Definition at line 36 of file Pythia6PartonEGun.cc.

References gen::Pythia6Gun::addAntiParticle(), eta(), funct::exp(), gen::Pythia6Gun::fEvt, fMaxE, fMaxEta, gen::Pythia6Gun::fMaxPhi, fMinE, fMinEta, gen::Pythia6Gun::fMinPhi, gen::Pythia6PartonGun::fPartonID, gen::Pythia6Gun::fPy6Service, configurableAnalysis::GenParticle, gen::Pythia6PartonGun::joinPartons(), gen::p, phi, gen::py1ent_(), gen::pyexec_(), gen::pymass_(), and gen::pyr_().

 {
    
    Pythia6Service::InstanceWrapper guard(fPy6Service);  // grab Py6 instance
 
    // now actualy, start cooking up the event gun 
    //
 
    // 1st, primary vertex
    //
    HepMC::GenVertex* Vtx = new HepMC::GenVertex( HepMC::FourVector(0.,0.,0.));
 
    // here re-create fEvt (memory)
    //
    fEvt = new HepMC::GenEvent() ;
      
    int ip=1;
    
    int py6PID = HepPID::translatePDTtoPythia( fPartonID );
    int dum = 0;
    double ee=0, the=0, eta=0;
    double mass = pymass_(py6PID);
      
    // fill p(ip,5) (in PYJETS) with mass value right now, 
    // because the (hardcoded) mstu(10)=1 will make py1ent
    // pick the mass from there
    pyjets.p[4][ip-1]=mass; 
          
    double phi = (fMaxPhi-fMinPhi)*pyr_(&dum)+fMinPhi;   
    ee   = (fMaxE-fMinE)*pyr_(&dum)+fMinE;
    eta  = (fMaxEta-fMinEta)*pyr_(&dum)+fMinEta;                                                      
    the  = 2.*atan(exp(-eta));                                                                          
      
    py1ent_(ip, py6PID, ee, the, phi);
      
    double px     = pyjets.p[0][ip-1]; // pt*cos(phi) ;
    double py     = pyjets.p[1][ip-1]; // pt*sin(phi) ;
    double pz     = pyjets.p[2][ip-1]; // mom*cos(the) ;
          
    HepMC::FourVector p(px,py,pz,ee) ;
    HepMC::GenParticle* Part = new HepMC::GenParticle(p,fPartonID,1);
    Part->suggest_barcode( ip ) ;
    Vtx->add_particle_out(Part);
      
    // now add anti-quark
    ip = ip + 1;
    HepMC::GenParticle* APart = addAntiParticle( ip, fPartonID, ee, eta, phi );   
    if ( APart ) 
    {
       Vtx->add_particle_out(APart) ;
    }
    else
    {
       // otherwise it should throw !
    }        
 
    // this should probably be configurable...
    //
    double qmax = 2.*ee;
    
    joinPartons( qmax );
          
    fEvt->add_vertex(Vtx);
      
    // run pythia
    pyexec_();
    
    return;
 }