---

edm::Pythia8Source Class Reference

#include <GeneratorInterface/Pythia8Interface/interface/Pythia8Source.h>

Inheritance diagram for edm::Pythia8Source:

List of all members.

Public Member Functions
void	endRun (Run &r)
	Pythia8Source (const ParameterSet &, const InputSourceDescription &)
	Constructor.
virtual	~Pythia8Source ()
	Destructor.
Private Member Functions
bool	call_txgive (const std::string &iParm)
	Interface to the PYGIVE pythia routine, with add'l protections.
bool	call_txgive_init ()
void	clear ()
virtual bool	produce (Event &e)
Private Attributes
double	comenergy
bool	doubleParticle
double	etamax
double	etamin
CLHEP::HepRandomEngine *	fRandomEngine
CLHEP::RandFlat *	fRandomGenerator
unsigned int	maxEventsToPrint_
	Events to print if verbosity.
int	particleID
double	phimax
double	phimin
double	ptmax
double	ptmin
Pythia8::Pythia *	pythia
Pythia8::Event *	pythia8event
bool	pythiaHepMCVerbosity_
	HepMC verbosity flag.
unsigned int	pythiaPylistVerbosity_
	Pythia PYLIST Verbosity flag.
HepMC::I_Pythia8 *	ToHepMC

---

Detailed Description

Definition at line 34 of file Pythia8Source.h.

---

Constructor & Destructor Documentation

Pythia8Source::Pythia8Source	(	const ParameterSet &	pset,
		const InputSourceDescription &	desc
	)

Constructor.

Definition at line 39 of file Pythia8Source.cc.

References comenergy, edm::errors::Configuration, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), cmsRelvalreport::exit, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init, maxEventsToPrint_, pars, pythia, pythia8event, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, and ToHepMC.

                                                                 :
  GeneratedInputSource(pset, desc),
  pythiaPylistVerbosity_ (pset.getUntrackedParameter<int>("pythiaPylistVerbosity",0)),
  pythiaHepMCVerbosity_ (pset.getUntrackedParameter<bool>("pythiaHepMCVerbosity",false)),
  maxEventsToPrint_ (pset.getUntrackedParameter<int>("maxEventsToPrint",1)),
  comenergy(pset.getUntrackedParameter<double>("comEnergy",14000.))
  
{
  
  cout << "Pythia8Source: initializing Pythia. " << endl;
  
  
  // PYLIST Verbosity Level
  // Valid PYLIST arguments are: 1, 2, 3, 5, 7, 11, 12, 13
  pythiaPylistVerbosity_ = pset.getUntrackedParameter<int>("pythiaPylistVerbosity",0);
  cout << "Pythia PYLIST verbosity level = " << pythiaPylistVerbosity_ << endl;
  
  // HepMC event verbosity Level
  pythiaHepMCVerbosity_ = pset.getUntrackedParameter<bool>("pythiaHepMCVerbosity",false);
  cout << "Pythia HepMC verbosity = " << pythiaHepMCVerbosity_ << endl; 

  //Max number of events printed on verbosity level 
  maxEventsToPrint_ = pset.getUntrackedParameter<int>("maxEventsToPrint",0);
  cout << "Number of events to be printed = " << maxEventsToPrint_ << endl;


  //In the future, we will get the random number seed on each event and tell 
  // pythia to use that new seed
    cout << "----------------------------------------------" << endl;
    cout << "Setting Pythia8 random number seed " << endl;
    cout << "----------------------------------------------" << endl;
  edm::Service<RandomNumberGenerator> rng;
  uint32_t seed = rng->mySeed();
  Pythia8::Rndm::init(seed);

    // Generator. Process selection. LHC initialization. Histogram.
  pythia = new Pythia8::Pythia;
  pythia8event = &(pythia->event);

  // Set PYTHIA parameters in a single ParameterSet
  ParameterSet pythia_params =
    pset.getParameter<ParameterSet>("PythiaParameters") ;

  // The parameter sets to be read (default, min bias, user ...) in the
  // proper order.
  vector<string> setNames =
    pythia_params.getParameter<vector<string> >("parameterSets");

  // Loop over the sets
  for ( unsigned i=0; i<setNames.size(); ++i ) {

    string mySet = setNames[i];

    // Read the PYTHIA parameters for each set of parameters
    vector<string> pars =
      pythia_params.getParameter<vector<string> >(mySet);

    if (mySet != "CSAParameters") {
      cout << "----------------------------------------------" << endl;
      cout << "Read PYTHIA parameter set " << mySet << endl;
      cout << "----------------------------------------------" << endl;

      // Loop over all parameters and stop in case of mistake
      for( vector<string>::const_iterator
           itPar = pars.begin(); itPar != pars.end(); ++itPar ) {

        if ( ! pythia->readString(*itPar) ) {
          throw edm::Exception(edm::errors::Configuration,"PythiaError") 
          <<" pythia8 did not accept the following \""<<*itPar<<"\"";
        }

      }

    } else {
      cout << " mySet == CSAParameters, not valid for Pythia8 " << endl;
      exit(1);
    }

  }

  pythia->init( 2212, 2212, comenergy);

  pythia->settings.listChanged();

  ToHepMC = new HepMC::I_Pythia8;
//  ToHepMC->set_crash_on_problem();

  cout << endl; // Stetically add for the output
  //********                                      
  
  produces<HepMCProduct>();
  produces<GenInfoProduct, edm::InRun>();

  cout << "Pythia8Source: starting event generation ... " << endl;

}

Pythia8Source::~Pythia8Source

(

)

[virtual]

Destructor.

Definition at line 138 of file Pythia8Source.cc.

References clear(), GenMuonPlsPt100GeV_cfg::cout, lat::endl(), pythia, and ToHepMC.

                             {
  cout << "PythiaSource: event generation done. " << endl;
  pythia->statistics();
  delete pythia;
  delete ToHepMC;
  clear();
}

---

Member Function Documentation

bool Pythia8Source::call_txgive

(

const std::string &

iParm

)

[private]

Interface to the PYGIVE pythia routine, with add'l protections.

Definition at line 239 of file Pythia8Source.cc.

                                                  {
  
//   TXGIVE( iParm.c_str(), iParm.length() );
//   cout << "     " <<  iParm.c_str() << endl; 

        return 1;  
}

bool Pythia8Source::call_txgive_init

(

)

[private]

Definition at line 248 of file Pythia8Source.cc.

                                {
  
//   TXGIVE_INIT();
//   cout << "  Setting CSA defaults.   "   << endl; 

        return 1;  
}

void Pythia8Source::clear

(

void

)

[private]

Definition at line 146 of file Pythia8Source.cc.

Referenced by ~Pythia8Source().

                          {
 
}

void Pythia8Source::endRun

(

Run &

r

)

[virtual]

Reimplemented from edm::ConfigurableInputSource.

Definition at line 151 of file Pythia8Source.cc.

References edm::Run::put(), and pythia.

                                  {

 double cs = pythia->info.sigmaGen(); // cross section in mb
 auto_ptr<GenInfoProduct> giprod (new GenInfoProduct());
 giprod->set_cross_section(cs);
// giprod->set_external_cross_section(extCrossSect);
// giprod->set_filter_efficiency(extFilterEff);
 r.put(giprod);

}

bool Pythia8Source::produce

(

Event &

e

)

[private, virtual]

Implements edm::ConfigurableInputSource.

Definition at line 163 of file Pythia8Source.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), edm::ConfigurableInputSource::event(), maxEventsToPrint_, edm::ConfigurableInputSource::numberEventsInRun(), edm::Event::put(), pythia, pythia8event, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, edm::InputSource::remainingEvents(), and ToHepMC.

                                     {

    auto_ptr<HepMCProduct> bare_product(new HepMCProduct());  
    //cout << "
    //********                                         
    //
    
    if (!pythia->next()) return false;  // generate one event with Pythia

    HepMC::GenEvent* hepmcevt = new HepMC::GenEvent();
    ToHepMC->fill_next_event( *pythia8event, hepmcevt );

    
//    hepmcevt->set_signal_process_id(pypars.msti[0]);
    hepmcevt->set_signal_process_id(pythia->info.code());
//    hepmcevt->set_event_scale(pypars.pari[16]);
    hepmcevt->set_event_scale(pythia->info.pTHat());
    hepmcevt->set_event_number(numberEventsInRun() - remainingEvents() - 1);

    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();
    hepmcevt->set_pdf_info( HepMC::PdfInfo(id1,id2,x1,x2,Q,pdf1,pdf2) ) ;

    hepmcevt->weights().push_back( pythia->info.weight() );

    //******** Verbosity ********
    
    if(event() <= maxEventsToPrint_ &&
       (pythiaPylistVerbosity_ || pythiaHepMCVerbosity_) ) {

      // Prints PYLIST info
      if(pythiaPylistVerbosity_) {
//      call_pylist(pythiaPylistVerbosity_);
        pythia->info.list();
        pythia->event.list();
    }
      
      // Prints HepMC event
      if(pythiaHepMCVerbosity_) {
//  cout << "Event process = " << pypars.msti[0] << endl 
    cout << "Event process = " << pythia->info.code() << endl
        << "----------------------" << endl;
        hepmcevt->print();
      }
    }
    

    if(hepmcevt)  bare_product->addHepMCData(hepmcevt);

    e.put(bare_product);

    return true;
}

---

Member Data Documentation

double edm::Pythia8Source::comenergy [private]

Definition at line 73 of file Pythia8Source.h.

Referenced by Pythia8Source().

bool edm::Pythia8Source::doubleParticle [private]

Definition at line 69 of file Pythia8Source.h.

double edm::Pythia8Source::etamax [private]

Definition at line 71 of file Pythia8Source.h.

double edm::Pythia8Source::etamin [private]

Definition at line 71 of file Pythia8Source.h.

CLHEP::HepRandomEngine* edm::Pythia8Source::fRandomEngine [private]

Definition at line 75 of file Pythia8Source.h.

CLHEP::RandFlat* edm::Pythia8Source::fRandomGenerator [private]

Definition at line 76 of file Pythia8Source.h.

unsigned int edm::Pythia8Source::maxEventsToPrint_ [private]

Events to print if verbosity.

Definition at line 61 of file Pythia8Source.h.

Referenced by produce(), and Pythia8Source().

int edm::Pythia8Source::particleID [private]

Definition at line 68 of file Pythia8Source.h.

double edm::Pythia8Source::phimax [private]

Definition at line 72 of file Pythia8Source.h.

double edm::Pythia8Source::phimin [private]

Definition at line 72 of file Pythia8Source.h.

double edm::Pythia8Source::ptmax [private]

Definition at line 70 of file Pythia8Source.h.

double edm::Pythia8Source::ptmin [private]

Definition at line 70 of file Pythia8Source.h.

Pythia8::Pythia* edm::Pythia8Source::pythia [private]

Definition at line 63 of file Pythia8Source.h.

Referenced by endRun(), produce(), Pythia8Source(), and ~Pythia8Source().

Pythia8::Event* edm::Pythia8Source::pythia8event [private]

Definition at line 64 of file Pythia8Source.h.

Referenced by produce(), and Pythia8Source().

bool edm::Pythia8Source::pythiaHepMCVerbosity_ [private]

HepMC verbosity flag.

Definition at line 59 of file Pythia8Source.h.

Referenced by produce(), and Pythia8Source().

unsigned int edm::Pythia8Source::pythiaPylistVerbosity_ [private]

Pythia PYLIST Verbosity flag.

Definition at line 57 of file Pythia8Source.h.

Referenced by produce(), and Pythia8Source().

HepMC::I_Pythia8* edm::Pythia8Source::ToHepMC [private]

Definition at line 65 of file Pythia8Source.h.

Referenced by produce(), Pythia8Source(), and ~Pythia8Source().

---

The documentation for this class was generated from the following files:

• GeneratorInterface/Pythia8Interface/interface/Pythia8Source.h
• GeneratorInterface/Pythia8Interface/src/Pythia8Source.cc

---