gen::Py8InterfaceBase Class Reference

#include <Py8InterfaceBase.h>

Inheritance diagram for gen::Py8InterfaceBase:

Public Member Functions
virtual const char *	classname () const =0
bool	decay ()
bool	declareSpecialSettings (const std::vector< std::string > &)
bool	declareStableParticles (const std::vector< int > &)
virtual void	finalizeEvent ()=0
virtual bool	generatePartonsAndHadronize ()=0
virtual bool	initializeForInternalPartons ()=0
void	p8SetRandomEngine (CLHEP::HepRandomEngine *v)
	Py8InterfaceBase (edm::ParameterSet const &ps)
P8RndmEngine &	randomEngine ()
bool	readSettings (int)
virtual bool	residualDecay ()
virtual void	statistics ()
	~Py8InterfaceBase ()

Protected Attributes
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

Private Attributes
P8RndmEngine	p8RndmEngine_

Detailed Description

Definition at line 19 of file Py8InterfaceBase.h.

Constructor & Destructor Documentation

gen::Py8InterfaceBase::Py8InterfaceBase

(

edm::ParameterSet const &

ps

)

Definition at line 15 of file Py8InterfaceBase.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), step1_ZMM_7Tev::maxEventsToPrint, step1_ZMM_7Tev::pythiaHepMCVerbosity, and step1_ZMM_7Tev::pythiaPylistVerbosity.

{
  fMasterGen.reset(new Pythia);
  fDecayer.reset(new Pythia);

  fMasterGen->readString("Next:numberShowEvent = 0");
  fDecayer->readString("Next:numberShowEvent = 0");

  fMasterGen->setRndmEnginePtr( &p8RndmEngine_ );
  fDecayer->setRndmEnginePtr( &p8RndmEngine_ );
  
  fParameters = ps.getParameter<edm::ParameterSet>("PythiaParameters");
  
  pythiaPylistVerbosity = ps.getUntrackedParameter<int>("pythiaPylistVerbosity", 0);
  pythiaHepMCVerbosity  = ps.getUntrackedParameter<bool>("pythiaHepMCVerbosity", false);
  maxEventsToPrint      = ps.getUntrackedParameter<int>("maxEventsToPrint", 0);

}

gen::Py8InterfaceBase::~Py8InterfaceBase ( )

inline

Definition at line 24 of file Py8InterfaceBase.h.

{}

Member Function Documentation

virtual const char* gen::Py8InterfaceBase::classname ( ) const

pure virtual

Implemented in Pythia8Hadronizer, gen::Py8EGun, gen::Py8JetGun, and gen::Py8PtGun.

bool gen::Py8InterfaceBase::decay ( )

inline

Definition at line 27 of file Py8InterfaceBase.h.

{ return true; } // NOT used - let's call it "design imperfection"

bool gen::Py8InterfaceBase::declareSpecialSettings

(

const std::vector< std::string > &

settings

)

Definition at line 79 of file Py8InterfaceBase.cc.

References spr::find().

{

   for ( unsigned int iss=0; iss<settings.size(); iss++ )
   {
      if ( settings[iss].find("QED-brem-off") == std::string::npos ) continue;
      fMasterGen->readString( "TimeShower:QEDshowerByL=off" );
   }

   return true;
}

bool gen::Py8InterfaceBase::declareStableParticles

(

const std::vector< int > &

pdgIds

)

Definition at line 54 of file Py8InterfaceBase.cc.

References i.

{
  
  for ( size_t i=0; i<pdgIds.size(); i++ )
  {
    // FIXME: need to double check if PID's are the same in Py6 & Py8,
    //        because the HepPDT translation tool is actually for **Py6** 
    // 
    // well, actually it looks like Py8 operates in PDT id's rather than Py6's
    //
    // int PyID = HepPID::translatePDTtoPythia( pdgIds[i] ); 
    int PyID = pdgIds[i]; 
    std::ostringstream pyCard ;
    pyCard << PyID <<":mayDecay=false";
    fMasterGen->readString( pyCard.str() );
    // alternative:
    // set the 2nd input argument warn=false 
    // - this way Py8 will NOT print warnings about unknown particle code(s)
    // fMasterPtr->readString( pyCard.str(), false )
  }
        
   return true;

}

virtual void gen::Py8InterfaceBase::finalizeEvent ( )

pure virtual

Implemented in Pythia8Hadronizer, and gen::Py8GunBase.

virtual bool gen::Py8InterfaceBase::generatePartonsAndHadronize ( )

pure virtual

Implemented in Pythia8Hadronizer, gen::Py8EGun, gen::Py8JetGun, and gen::Py8PtGun.

virtual bool gen::Py8InterfaceBase::initializeForInternalPartons ( )

pure virtual

Implemented in Pythia8Hadronizer, and gen::Py8GunBase.

void gen::Py8InterfaceBase::p8SetRandomEngine ( CLHEP::HepRandomEngine *  v )

inline

Definition at line 37 of file Py8InterfaceBase.h.

References p8RndmEngine_, and gen::P8RndmEngine::setRandomEngine().

Referenced by gen::Py8GunBase::setRandomEngine().

{ p8RndmEngine_.setRandomEngine(v); }

P8RndmEngine& gen::Py8InterfaceBase::randomEngine ( )

inline

Definition at line 38 of file Py8InterfaceBase.h.

References p8RndmEngine_.

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

{ return p8RndmEngine_; }

bool gen::Py8InterfaceBase::readSettings

(

int

)

Definition at line 34 of file Py8InterfaceBase.cc.

References geometryCSVtoXML::line.

{

   for ( ParameterCollector::const_iterator line = fParameters.begin();
         line != fParameters.end(); ++line ) 
   {
      if (line->find("Random:") != std::string::npos)
         throw cms::Exception("PythiaError") << "Attempted to set random number "
         "using Pythia commands. Please use " "the RandomNumberGeneratorService."
         << std::endl;

      if (!fMasterGen->readString(*line)) throw cms::Exception("PythiaError")
                          << "Pythia 8 did not accept \""
                      << *line << "\"." << std::endl;
   }

   return true;

}

bool gen::Py8InterfaceBase::residualDecay ( )

virtual

Reimplemented in gen::Py8GunBase.

Definition at line 91 of file Py8InterfaceBase.cc.

{
  
/*
  Event* pythiaEvent = &(fMasterPtr->event);
  
  assert(fCurrentEventState);
  
  int NPartsBeforeDecays = pythiaEvent->size();
  // int NPartsAfterDecays = event().get()->particles_size();
  int NPartsAfterDecays = fCurrentEventState->particles_size();
  int NewBarcode = NPartsAfterDecays;
   
  for ( int ipart=NPartsAfterDecays; ipart>NPartsBeforeDecays; ipart-- )
  {

    // HepMC::GenParticle* part = event().get()->barcode_to_particle( ipart );
    HepMC::GenParticle* part = fCurrentEventState->barcode_to_particle( ipart );

    if ( part->status() == 1 )
    {
      fDecayerPtr->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( (fDecayerPtr->particleData).tau0( part->pdg_id() ) );
      fDecayerPtr->event.append( py8part );
      int nentries = fDecayerPtr->event.size();
      if ( !fDecayerPtr->event[nentries-1].mayDecay() ) continue;
      fDecayerPtr->next();
      int nentries1 = fDecayerPtr->event.size();
      // --> fDecayerPtr->event.list(std::cout);
      if ( nentries1 <= nentries ) continue; //same number of particles, no decays...
        
      part->set_status(2);
        
      Particle& py8daughter = fDecayerPtr->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 ( ipart=nentries+1; ipart<nentries1; ipart++ )
      {
        py8daughter = fDecayerPtr->event[ipart];
        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 );
      fCurrentEventState->add_vertex( DecVtx );

    }
 } 
*/   
 return true;
 
}

void gen::Py8InterfaceBase::statistics ( )

virtual

Reimplemented in Pythia8Hadronizer, and gen::Py8GunBase.

Definition at line 173 of file Py8InterfaceBase.cc.

{
  
   fMasterGen->statistics();
   return;
   
}

Member Data Documentation

std::auto_ptr<Pythia8::Pythia> gen::Py8InterfaceBase::fDecayer

protected

Definition at line 43 of file Py8InterfaceBase.h.

Referenced by gen::Py8GunBase::initializeForInternalPartons(), and gen::Py8GunBase::residualDecay().

std::auto_ptr<Pythia8::Pythia> gen::Py8InterfaceBase::fMasterGen

protected

Definition at line 42 of file Py8InterfaceBase.h.

Referenced by gen::Py8GunBase::finalizeEvent(), Pythia8Hadronizer::finalizeEvent(), gen::Py8EGun::generatePartonsAndHadronize(), gen::Py8JetGun::generatePartonsAndHadronize(), gen::Py8PtGun::generatePartonsAndHadronize(), Pythia8Hadronizer::generatePartonsAndHadronize(), Pythia8Hadronizer::hadronize(), Pythia8Hadronizer::initializeForExternalPartons(), gen::Py8GunBase::initializeForInternalPartons(), Pythia8Hadronizer::initializeForInternalPartons(), Pythia8Hadronizer::Pythia8Hadronizer(), gen::Py8GunBase::residualDecay(), gen::Py8GunBase::statistics(), and Pythia8Hadronizer::statistics().

ParameterCollector gen::Py8InterfaceBase::fParameters

protected

Definition at line 45 of file Py8InterfaceBase.h.

Referenced by Pythia8Hadronizer::Pythia8Hadronizer().

unsigned int gen::Py8InterfaceBase::maxEventsToPrint

protected

Definition at line 49 of file Py8InterfaceBase.h.

Referenced by gen::Py8GunBase::finalizeEvent(), and Pythia8Hadronizer::finalizeEvent().

P8RndmEngine gen::Py8InterfaceBase::p8RndmEngine_

private

Definition at line 53 of file Py8InterfaceBase.h.

Referenced by p8SetRandomEngine(), and randomEngine().

bool gen::Py8InterfaceBase::pythiaHepMCVerbosity

protected

Definition at line 48 of file Py8InterfaceBase.h.

Referenced by gen::Py8GunBase::finalizeEvent(), and Pythia8Hadronizer::finalizeEvent().

unsigned int gen::Py8InterfaceBase::pythiaPylistVerbosity

protected

Definition at line 47 of file Py8InterfaceBase.h.

Referenced by gen::Py8GunBase::finalizeEvent(), Pythia8Hadronizer::finalizeEvent(), Pythia8Hadronizer::initializeForExternalPartons(), and Pythia8Hadronizer::initializeForInternalPartons().

HepMC::I_Pythia8 gen::Py8InterfaceBase::toHepMC

protected

Definition at line 44 of file Py8InterfaceBase.h.

Referenced by gen::Py8EGun::generatePartonsAndHadronize(), gen::Py8JetGun::generatePartonsAndHadronize(), gen::Py8PtGun::generatePartonsAndHadronize(), Pythia8Hadronizer::generatePartonsAndHadronize(), and Pythia8Hadronizer::hadronize().