#include <Py8GunBase.h>

Inheritance diagram for gen::Py8GunBase:

Public Member Functions
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	setRandomEngine (CLHEP::HepRandomEngine *v)

std::vector< std::string > const &	sharedResources () const

void	statistics ()

	~Py8GunBase ()

Public Member Functions inherited from gen::Py8InterfaceBase
virtual const char *	classname () const =0

bool	decay ()

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

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

virtual bool	generatePartonsAndHadronize ()=0

void	p8SetRandomEngine (CLHEP::HepRandomEngine *v)

	Py8InterfaceBase (edm::ParameterSet const &ps)

P8RndmEngine &	randomEngine ()

bool	readSettings (int)

	~Py8InterfaceBase ()

Protected Member Functions
std::auto_ptr< HepMC::GenEvent > &	event ()

std::auto_ptr < GenEventInfoProduct > &	eventInfo ()

GenRunInfoProduct &	runInfo ()

Protected Attributes
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::Pythia8ToHepMC	toHepMC

Private Attributes
edm::Event *	edmEvent_

std::auto_ptr< HepMC::GenEvent >	genEvent_

std::auto_ptr < GenEventInfoProduct >	genEventInfo_

GenRunInfoProduct	genRunInfo_

Static Private Attributes
static const std::vector < std::string >	p8SharedResources = { edm::SharedResourceNames::kPythia8 }

Detailed Description

Definition at line 41 of file Py8GunBase.h.

Constructor & Destructor Documentation

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

Definition at line 10 of file Py8GunBase.cc.

References fMaxPhi, fMinPhi, fPartIDs, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), runInfo(), GenRunInfoProduct::setExternalXSecLO(), GenRunInfoProduct::setExternalXSecNLO(), and GenRunInfoProduct::setFilterEfficiency().

    : Py8InterfaceBase(ps)
 {
 
    runInfo().setFilterEfficiency(
       ps.getUntrackedParameter<double>("filterEfficiency", -1.) );
    runInfo().setExternalXSecLO(
       GenRunInfoProduct::XSec(ps.getUntrackedParameter<double>("crossSection", -1.)) );
    runInfo().setExternalXSecNLO(
        GenRunInfoProduct::XSec(ps.getUntrackedParameter<double>("crossSectionNLO", -1.)) );
 
   
   // PGun specs
   //
    edm::ParameterSet pgun_params = 
       ps.getParameter<edm::ParameterSet>("PGunParameters");
       
    // although there's the method ParameterSet::empty(),  
    // it looks like it's NOT even necessary to check if it is,
    // before trying to extract parameters - if it is empty,
    // the default values seem to be taken
    //
    fPartIDs    = pgun_params.getParameter< std::vector<int> >("ParticleID");
    fMinPhi     = pgun_params.getParameter<double>("MinPhi"); // ,-3.14159265358979323846);
    fMaxPhi     = pgun_params.getParameter<double>("MaxPhi"); // , 3.14159265358979323846);
    
 }

gen::Py8GunBase::~Py8GunBase ( )

inline

Definition at line 44 of file Py8GunBase.h.

44 {}

Member Function Documentation

std::auto_ptr<HepMC::GenEvent>& gen::Py8GunBase::event ( )

inlineprotected

Definition at line 69 of file Py8GunBase.h.

References genEvent_.

Referenced by finalizeEvent(), gen::Py8JetGun::generatePartonsAndHadronize(), gen::Py8PtGun::generatePartonsAndHadronize(), gen::Py8EGun::generatePartonsAndHadronize(), residualDecay(), and setEDMEvent().

69 { return genEvent_; }

gen::Py8GunBase::genEvent_

std::auto_ptr< HepMC::GenEvent > genEvent_

Definition: Py8GunBase.h:80

std::auto_ptr<GenEventInfoProduct>& gen::Py8GunBase::eventInfo ( )

inlineprotected

Definition at line 70 of file Py8GunBase.h.

References genEventInfo_.

70 { return genEventInfo_; }

gen::Py8GunBase::genEventInfo_

std::auto_ptr< GenEventInfoProduct > genEventInfo_

Definition: Py8GunBase.h:81

void gen::Py8GunBase::finalizeEvent ( )

virtual

Implements gen::Py8InterfaceBase.

Definition at line 137 of file Py8GunBase.cc.

References gather_cfg::cout, event(), gen::Py8InterfaceBase::fMasterGen, gen::Py8InterfaceBase::maxEventsToPrint, gen::Py8InterfaceBase::pythiaHepMCVerbosity, and gen::Py8InterfaceBase::pythiaPylistVerbosity.

 {
    
 
 // FIXME !!!
 
   //******** Verbosity ********
 
    if (maxEventsToPrint > 0 &&
       (pythiaPylistVerbosity || pythiaHepMCVerbosity)) 
    {
       maxEventsToPrint--;
       if (pythiaPylistVerbosity) 
       {
          fMasterGen->info.list(std::cout); 
          fMasterGen->event.list(std::cout);
       } 
 
       if (pythiaHepMCVerbosity) 
       {
          std::cout << "Event process = "
                    << fMasterGen->info.code() << "\n"
                    << "----------------------" << std::endl;
          event()->print();
       }
    }
       
    return;
 }

edm::Event& gen::Py8GunBase::getEDMEvent ( ) const

inline

Definition at line 56 of file Py8GunBase.h.

References edmEvent_.

56 { return *edmEvent_; }

gen::Py8GunBase::edmEvent_

edm::Event * edmEvent_

Definition: Py8GunBase.h:83

HepMC::GenEvent* gen::Py8GunBase::getGenEvent ( )

inline

Definition at line 48 of file Py8GunBase.h.

References genEvent_.

48 { return genEvent_.release(); }

gen::Py8GunBase::genEvent_

std::auto_ptr< HepMC::GenEvent > genEvent_

Definition: Py8GunBase.h:80

GenEventInfoProduct* gen::Py8GunBase::getGenEventInfo ( )

inline

Definition at line 49 of file Py8GunBase.h.

References genEventInfo_.

49 { return genEventInfo_.release(); }

gen::Py8GunBase::genEventInfo_

std::auto_ptr< GenEventInfoProduct > genEventInfo_

Definition: Py8GunBase.h:81

GenRunInfoProduct& gen::Py8GunBase::getGenRunInfo ( )

inline

Definition at line 47 of file Py8GunBase.h.

References genRunInfo_.

47 { return genRunInfo_; }

gen::Py8GunBase::genRunInfo_

GenRunInfoProduct genRunInfo_

Definition: Py8GunBase.h:79

bool gen::Py8GunBase::initializeForInternalPartons ( )

virtual

Implements gen::Py8InterfaceBase.

Definition at line 40 of file Py8GunBase.cc.

References gen::Py8InterfaceBase::fDecayer, and gen::Py8InterfaceBase::fMasterGen.

 {
 
    // NO MATTER what was this setting below, override it before init 
    // - this is essencial for the PGun mode 
    
    // Key requirement: switch off ProcessLevel, and thereby also PartonLevel.
    fMasterGen->readString("ProcessLevel:all = off");
    fMasterGen->readString("ProcessLevel::resonanceDecays=on");
    fMasterGen->init();
    
    // init decayer
    fDecayer->readString("ProcessLevel:all = off"); // Same trick!
    fDecayer->readString("ProcessLevel::resonanceDecays=on");
    fDecayer->init();
   
    return true;
 
 }

void gen::Py8GunBase::resetEvent ( HepMC::GenEvent * event )

inline

Definition at line 51 of file Py8GunBase.h.

References genEvent_.

51 { genEvent_.reset(event); }

gen::Py8GunBase::genEvent_

std::auto_ptr< HepMC::GenEvent > genEvent_

Definition: Py8GunBase.h:80

gen::Py8GunBase::event

std::auto_ptr< HepMC::GenEvent > & event()

Definition: Py8GunBase.h:69

void gen::Py8GunBase::resetEventInfo ( GenEventInfoProduct * eventInfo )

inline

Definition at line 52 of file Py8GunBase.h.

References genEventInfo_.

52 { genEventInfo_.reset(eventInfo); }

gen::Py8GunBase::genEventInfo_

std::auto_ptr< GenEventInfoProduct > genEventInfo_

Definition: Py8GunBase.h:81

bool gen::Py8GunBase::residualDecay ( )

virtual

Definition at line 60 of file Py8GunBase.cc.

References event(), gen::Py8InterfaceBase::fDecayer, gen::Py8InterfaceBase::fMasterGen, and configurableAnalysis::GenParticle.

 {
 
    Event* pythiaEvent = &(fMasterGen->event);
   
    int NPartsBeforeDecays = pythiaEvent->size()-1; // do NOT count the very 1st "system" particle 
                                                    // in Pythia8::Event record; it does NOT even
                            // get translated by the HepMCInterface to the
                            // HepMC::GenEvent record !!!
    int NPartsAfterDecays = event().get()->particles_size();
    int NewBarcode = NPartsAfterDecays;
    
    for ( int ipart=NPartsAfterDecays; ipart>NPartsBeforeDecays; ipart-- )
    {
 
       HepMC::GenParticle* part = event().get()->barcode_to_particle( ipart );
 
       if ( part->status() == 1 )
       {
          fDecayer->event.reset();
          Particle py8part(  part->pdg_id(), 93, 0, 0, 0, 0, 0, 0,
                          part->momentum().x(),
                          part->momentum().y(),
                          part->momentum().z(),
                          part->momentum().t(),
                          part->generated_mass() );
          HepMC::GenVertex* ProdVtx = part->production_vertex();
          py8part.vProd( ProdVtx->position().x(), ProdVtx->position().y(), 
                      ProdVtx->position().z(), ProdVtx->position().t() );
          py8part.tau( (fDecayer->particleData).tau0( part->pdg_id() ) );
          fDecayer->event.append( py8part );
          int nentries = fDecayer->event.size();
          if ( !fDecayer->event[nentries-1].mayDecay() ) continue;
          fDecayer->next();
          int nentries1 = fDecayer->event.size();
          if ( nentries1 <= nentries ) continue; //same number of particles, no decays...
         
          part->set_status(2);
         
          Particle& py8daughter = fDecayer->event[nentries]; // the 1st daughter
          HepMC::GenVertex* DecVtx = new HepMC::GenVertex( HepMC::FourVector(py8daughter.xProd(),
                                                           py8daughter.yProd(),
                                                           py8daughter.zProd(),
                                                           py8daughter.tProd()) );
 
          DecVtx->add_particle_in( part ); // this will cleanup end_vertex if exists, replace with the new one
                                           // I presume (vtx) barcode will be given automatically
         
          HepMC::FourVector pmom( py8daughter.px(), py8daughter.py(), py8daughter.pz(), py8daughter.e() );
         
          HepMC::GenParticle* daughter =
                              new HepMC::GenParticle( pmom, py8daughter.id(), 1 );
         
          NewBarcode++;
          daughter->suggest_barcode( NewBarcode );
          DecVtx->add_particle_out( daughter );
                 
          for ( int ipart1=nentries+1; ipart1<nentries1; ipart1++ )
          {
             py8daughter = fDecayer->event[ipart1];
             HepMC::FourVector pmomN( py8daughter.px(), py8daughter.py(), py8daughter.pz(), py8daughter.e() );       
             HepMC::GenParticle* daughterN =
                                 new HepMC::GenParticle( pmomN, py8daughter.id(), 1 );
             NewBarcode++;
             daughterN->suggest_barcode( NewBarcode );
             DecVtx->add_particle_out( daughterN );
          }
         
          event().get()->add_vertex( DecVtx );
 
       }
    } 
 
    return true;
 
 }

GenRunInfoProduct& gen::Py8GunBase::runInfo ( )

inlineprotected

Definition at line 68 of file Py8GunBase.h.

References genRunInfo_.

Referenced by Py8GunBase(), and statistics().

68 { return genRunInfo_; }

gen::Py8GunBase::genRunInfo_

GenRunInfoProduct genRunInfo_

Definition: Py8GunBase.h:79

virtual bool gen::Py8GunBase::select ( HepMC::GenEvent * ) const

inlinevirtual

Definition at line 57 of file Py8GunBase.h.

57 { return true;}

void gen::Py8GunBase::setEDMEvent ( edm::Event & event )

inline

Definition at line 55 of file Py8GunBase.h.

References edmEvent_, and event().

55 { edmEvent_ = &event; }

gen::Py8GunBase::edmEvent_

edm::Event * edmEvent_

Definition: Py8GunBase.h:83

gen::Py8GunBase::event

std::auto_ptr< HepMC::GenEvent > & event()

Definition: Py8GunBase.h:69

void gen::Py8GunBase::setRandomEngine ( CLHEP::HepRandomEngine * v )

inline

Definition at line 64 of file Py8GunBase.h.

References gen::Py8InterfaceBase::p8SetRandomEngine().

64 { p8SetRandomEngine(v); }

gen::Py8InterfaceBase::p8SetRandomEngine

void p8SetRandomEngine(CLHEP::HepRandomEngine *v)

Definition: Py8InterfaceBase.h:36

gen::v

double v[5][pyjets_maxn]

Definition: Cascade2Hadronizer.cc:80

std::vector<std::string> const& gen::Py8GunBase::sharedResources ( ) const

inline

Definition at line 65 of file Py8GunBase.h.

References p8SharedResources.

65 { return p8SharedResources; }

gen::Py8GunBase::p8SharedResources

static const std::vector< std::string > p8SharedResources

Definition: Py8GunBase.h:85

void gen::Py8GunBase::statistics ( )

virtual

Reimplemented from gen::Py8InterfaceBase.

Definition at line 167 of file Py8GunBase.cc.

References gen::Py8InterfaceBase::fMasterGen, runInfo(), and GenRunInfoProduct::setInternalXSec().

 {
 
    fMasterGen->statistics();
 
    double xsec = fMasterGen->info.sigmaGen(); // cross section in mb
    xsec *= 1.0e9; // translate to pb (CMS/Gen "convention" as of May 2009)
    runInfo().setInternalXSec(xsec);
    return;
 
 }