Inheritance diagram for Pythia8Hadronizer:

Public Member Functions
const char *	classname () const

bool	decay ()

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

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

void	finalizeEvent ()

bool	generatePartonsAndHadronize ()

bool	hadronize ()

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons ()

	Pythia8Hadronizer (const edm::ParameterSet &params)

bool	residualDecay ()

void	statistics ()

	~Pythia8Hadronizer ()

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

edm::Event &	getEDMEvent () const

HepMC::GenEvent *	getGenEvent ()

GenEventInfoProduct *	getGenEventInfo ()

GenRunInfoProduct &	getGenRunInfo ()

const boost::shared_ptr < lhef::LHERunInfo > &	getLHERunInfo () const

void	resetEvent (HepMC::GenEvent *event)

void	resetEventInfo (GenEventInfoProduct *eventInfo)

void	setEDMEvent (edm::Event &event)

void	setLHEEvent (lhef::LHEEvent *event)

void	setLHERunInfo (lhef::LHERunInfo *runInfo)

	~BaseHadronizer ()

Private Attributes
double	comEnergy
	Center-of-Mass energy. More...

std::auto_ptr< LHAupLesHouches >	lhaUP

string	LHEInputFileName

unsigned int	maxEventsToPrint
	Events to print if verbosity. More...

ParameterCollector	parameters

std::auto_ptr< Pythia >	pythia

Event *	pythiaEvent

bool	pythiaHepMCVerbosity
	HepMC verbosity flag. More...

unsigned int	pythiaPylistVerbosity
	Pythia PYLIST Verbosity flag. More...

HepMC::I_Pythia8	toHepMC

UserHooks *	UserHookPointer

bool	useUserHook
	Switch User Hook flag. More...

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

std::auto_ptr < GenEventInfoProduct > &	eventInfo ()

lhef::LHEEvent *	lheEvent ()

lhef::LHERunInfo *	lheRunInfo ()

GenRunInfoProduct &	runInfo ()

Detailed Description

Definition at line 39 of file Pythia8Hadronizer.cc.

Constructor & Destructor Documentation

Pythia8Hadronizer::Pythia8Hadronizer ( const edm::ParameterSet & params )

Definition at line 87 of file Pythia8Hadronizer.cc.

References edm::ParameterSet::exists(), gen::getEngineReference(), edm::ParameterSet::getParameter(), randomEngine, and useUserHook.

                                                                   :
         BaseHadronizer(params),
     parameters(params.getParameter<edm::ParameterSet>("PythiaParameters")),
     comEnergy(params.getParameter<double>("comEnergy")),
     pythiaPylistVerbosity(params.getUntrackedParameter<int>("pythiaPylistVerbosity", 0)),
     pythiaHepMCVerbosity(params.getUntrackedParameter<bool>("pythiaHepMCVerbosity", false)),
     maxEventsToPrint(params.getUntrackedParameter<int>("maxEventsToPrint", 0)),
     LHEInputFileName(params.getUntrackedParameter<string>("LHEInputFileName","")),
     useUserHook(false),
     UserHookPointer(0)
 {
     if( params.exists( "useUserHook" ) )
       useUserHook = params.getParameter<bool>("useUserHook");
     randomEngine = &getEngineReference();
 }

Pythia8Hadronizer::~Pythia8Hadronizer ( )

Definition at line 103 of file Pythia8Hadronizer.cc.

104 {

105 }

Member Function Documentation

const char* Pythia8Hadronizer::classname ( ) const

inline

Definition at line 58 of file Pythia8Hadronizer.cc.

58 { return "Pythia8Hadronizer"; }

bool Pythia8Hadronizer::decay ( )

Definition at line 292 of file Pythia8Hadronizer.cc.

 {
     return true;
 }

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

Definition at line 245 of file Pythia8Hadronizer.cc.

 {
    return true;
 }

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

Definition at line 224 of file Pythia8Hadronizer.cc.

References i, and pythia.

 {
 
    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** 
       // 
       int PyID = HepPID::translatePDTtoPythia( pdgIds[i] ); 
       std::ostringstream pyCard ;
       pyCard << PyID <<":mayDecay=false";
       pythia->readString( pyCard.str() );
       // alternative:
       // set the 2nd input argument warn=false 
       // - this way Py8 will NOT print warnings about unknown particle code(s)
       // pythia->readString( pyCard.str(), false )
    }
    
    return true;
 
 }

void Pythia8Hadronizer::finalizeEvent ( )

Definition at line 302 of file Pythia8Hadronizer.cc.

References abs, gather_cfg::cout, gen::BaseHadronizer::event(), gen::BaseHadronizer::eventInfo(), gen::BaseHadronizer::lheEvent(), gen::BaseHadronizer::lheRunInfo(), maxEventsToPrint, pythia, pythiaHepMCVerbosity, pythiaPylistVerbosity, UserHookPointer, and useUserHook.

 {
   bool lhe = lheEvent() != 0;
 
   event()->set_signal_process_id(pythia->info.code());
   event()->set_event_scale(pythia->info.pTHat());   //FIXME
 
   // Putting pdf info into the HepMC record
   // There is the overloaded pythia8 HepMCInterface method fill_next_event
   // that does this, but CMSSW GeneratorInterface does not fill HepMC
   // record according to the HepMC convention (stores f(x) instead of x*f(x)
   // and converts gluon PDG ID to zero). For this reason we use the
   // method fill_next_event (above) that does NOT this and fill pdf info here
   //
   int id1 = pythia->info.id1();
   int id2 = pythia->info.id2();
   if (id1 == 21) id1 = 0;
   if (id2 == 21) id2 = 0;
   double x1 = pythia->info.x1();
   double x2 = pythia->info.x2();
   double Q = pythia->info.QRen();
   double pdf1 = pythia->info.pdf1() / pythia->info.x1();
   double pdf2 = pythia->info.pdf2() / pythia->info.x2();
   event()->set_pdf_info(HepMC::PdfInfo(id1,id2,x1,x2,Q,pdf1,pdf2));
 
   // Storing weights. Will be moved to pythia8 HepMCInterface
   //
   if (lhe && std::abs(lheRunInfo()->getHEPRUP()->IDWTUP) == 4)
     // translate mb to pb (CMS/Gen "convention" as of May 2009)
     event()->weights().push_back( pythia->info.weight() * 1.0e9 );
   else
     event()->weights().push_back( pythia->info.weight() );
   //
   // Below is the weight to be used in the case of reweighting by UserHook
   // This is a temporary solution. It requires the explicit conversion of
   // the pointer. In future versions (>150) the additional factor introduced
   // in the UserHook will be stored by pythia. Also, putting the weights
   // into the HepMC record will be made in the pythia8 HepMCInterface
   //
   if(useUserHook)
     event()->weights().push_back(1./((PtHatReweightUserHook*)UserHookPointer)->getFactor());
 
   // now create the GenEventInfo product from the GenEvent and fill
   // the missing pieces
   eventInfo().reset( new GenEventInfoProduct( event().get() ) );
 
   // in pythia pthat is used to subdivide samples into different bins
   // in LHE mode the binning is done by the external ME generator
   // which is likely not pthat, so only filling it for Py6 internal mode
   if (!lhe) {
     eventInfo()->setBinningValues(std::vector<double>(1, pythia->info.pTHat()));
   }
 
   //******** Verbosity ********
 
   if (maxEventsToPrint > 0 &&
       (pythiaPylistVerbosity || pythiaHepMCVerbosity)) {
     maxEventsToPrint--;
     if (pythiaPylistVerbosity) {
       pythia->info.list(std::cout); 
       pythia->event.list(std::cout);
     } 
 
     if (pythiaHepMCVerbosity) {
             std::cout << "Event process = "
                       << pythia->info.code() << "\n"
                       << "----------------------" << std::endl;
             event()->print();
         }
     }
 }

bool Pythia8Hadronizer::generatePartonsAndHadronize ( )

Definition at line 260 of file Pythia8Hadronizer.cc.

References gen::BaseHadronizer::event(), pythia, pythiaEvent, and toHepMC.

 {
     if (!pythia->next())
         return false;
 
     event().reset(new HepMC::GenEvent);
     toHepMC.fill_next_event(*pythiaEvent, event().get());
 
     return true;
 }

bool Pythia8Hadronizer::hadronize ( )

Definition at line 271 of file Pythia8Hadronizer.cc.

References lhef::LHEEvent::count(), gen::BaseHadronizer::event(), lhef::LHERunInfo::kAccepted, lhaUP, gen::BaseHadronizer::lheEvent(), LHEInputFileName, pythia, pythiaEvent, and toHepMC.

 {
     if(LHEInputFileName == string()) {
       //cout << "start loading event" << endl;
       lhaUP->loadEvent(lheEvent());
       //cout << "finish loading event" << endl;
     }
 
     if (!pythia->next())
         return false;
 
     // update LHE matching statistics
     //
     lheEvent()->count( lhef::LHERunInfo::kAccepted );
 
     event().reset(new HepMC::GenEvent);
     toHepMC.fill_next_event(*pythiaEvent, event().get());
 
     return true;
 }

bool Pythia8Hadronizer::initializeForExternalPartons ( )

Definition at line 156 of file Pythia8Hadronizer.cc.

References gen::ParameterCollector::begin(), gather_cfg::cout, gen::ParameterCollector::end(), lhaUP, LHEInputFileName, gen::BaseHadronizer::lheRunInfo(), geometryCSVtoXML::line, parameters, pythia, pythiaEvent, and pythiaPylistVerbosity.

 {
 
     std::cout << "Initializing for external partons" << std::endl;
 
     RandomP8* RP8 = new RandomP8();
 
     pythia.reset(new Pythia);
 
     pythia->setRndmEnginePtr(RP8);
 
     pythiaEvent = &pythia->event;
 
     for(ParameterCollector::const_iterator line = parameters.begin();
         line != parameters.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 (!pythia->readString(*line))
             throw cms::Exception("PythiaError")
                 << "Pythia 8 did not accept \""
                 << *line << "\"." << std::endl;
     }
 
     if(pythiaPylistVerbosity > 10) {
       if(pythiaPylistVerbosity == 11 || pythiaPylistVerbosity == 13)
         pythia->settings.listAll();
       if(pythiaPylistVerbosity == 12 || pythiaPylistVerbosity == 13)
         pythia->particleData.listAll();
     }
 
     if(LHEInputFileName != string()) {
 
       cout << endl;
       cout << "LHE Input File Name = " << LHEInputFileName << endl;
       cout << endl;
       pythia->init(LHEInputFileName);
 
     } else {
 
       lhaUP.reset(new LHAupLesHouches());
       lhaUP->loadRunInfo(lheRunInfo());
       pythia->init(lhaUP.get());
 
     }
 
     return true;
 }

bool Pythia8Hadronizer::initializeForInternalPartons ( )

Definition at line 107 of file Pythia8Hadronizer.cc.

References gen::ParameterCollector::begin(), comEnergy, gen::ParameterCollector::end(), geometryCSVtoXML::line, parameters, pythia, pythiaEvent, pythiaPylistVerbosity, UserHookPointer, and useUserHook.

 {
     //Old code that used Pythia8 own random engine
     //edm::Service<edm::RandomNumberGenerator> rng;
     //uint32_t seed = rng->mySeed();
     //Pythia8::Rndm::init(seed);
 
     RandomP8* RP8 = new RandomP8();
 
     pythia.reset(new Pythia);
 
     pythia->setRndmEnginePtr(RP8);
     if(useUserHook){
       UserHookPointer = new PtHatReweightUserHook();
       pythia->setUserHooksPtr(UserHookPointer);
     }
 
     pythiaEvent = &pythia->event;
 
     for(ParameterCollector::const_iterator line = parameters.begin();
         line != parameters.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 (!pythia->readString(*line))
             throw cms::Exception("PythiaError")
                 << "Pythia 8 did not accept \""
                 << *line << "\"." << std::endl;
     }
 
     if(pythiaPylistVerbosity > 10) {
       if(pythiaPylistVerbosity == 11 || pythiaPylistVerbosity == 13)
         pythia->settings.listAll();
       if(pythiaPylistVerbosity == 12 || pythiaPylistVerbosity == 13)
         pythia->particleData.listAll();
     }
 
     pythia->init(2212, 2212, comEnergy);
 
     pythia->settings.listChanged();
 
     return true;
 }

bool Pythia8Hadronizer::residualDecay ( )

Definition at line 297 of file Pythia8Hadronizer.cc.

 {
     return true;
 }

void Pythia8Hadronizer::statistics ( )

Definition at line 251 of file Pythia8Hadronizer.cc.

References pythia, gen::BaseHadronizer::runInfo(), and GenRunInfoProduct::setInternalXSec().

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