---

edm::ToprexSource Class Reference

#include <GeneratorInterface/TopRexInterface/interface/ToprexSource.h>

Inheritance diagram for edm::ToprexSource:

List of all members.

Public Member Functions
	ToprexSource (const ParameterSet &, const InputSourceDescription &)
	Constructor.
virtual	~ToprexSource ()
	Destructor.
Private Member Functions
bool	call_pygive (const std::string &iParm)
	Interface to the PYGIVE/TXGIVE pythia routine, with add'l protections.
bool	call_txgive (const std::string &iParm)
bool	call_txgive_init ()
void	clear ()
virtual bool	produce (Event &e)
Private Attributes
HepMC::GenEvent *	evt
CLHEP::HepRandomEngine *	fRandomEngine
CLHEP::RandFlat *	fRandomGenerator
unsigned int	maxEventsToPrint_
	Events to print if verbosity.
bool	pythiaHepMCVerbosity_
	HepMC verbosity flag.
unsigned int	pythiaPylistVerbosity_
	Pythia PYLIST Verbosity flag.

---

Detailed Description

Definition at line 34 of file ToprexSource.h.

---

Constructor & Destructor Documentation

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

Constructor.

Definition at line 49 of file ToprexSource.cc.

References call_pygive(), call_toprex(), call_txgive(), call_txgive_init(), edm::errors::Configuration, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), fRandomEngine, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, maxEventsToPrint_, pars, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, and randomEngine.

                                                                 :
  GeneratedInputSource(pset, desc), evt(0), 
  pythiaPylistVerbosity_ (pset.getUntrackedParameter<int>("pythiaPylistVerbosity",0)),
  pythiaHepMCVerbosity_ (pset.getUntrackedParameter<bool>("pythiaHepMCVerbosity",false)),
  maxEventsToPrint_ (pset.getUntrackedParameter<int>("maxEventsToPrint",1))
  
{
  
  cout << "ToprexSource: initializing TopReX " << 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;

  // 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 ) {
      static string sRandomValueSetting("MRPY(1)");
      if( 0 == itPar->compare(0,sRandomValueSetting.size(),sRandomValueSetting) ) {
        throw edm::Exception(edm::errors::Configuration,"PythiaError")
          <<" attempted to set random number using pythia command 'MRPY(1)' this is not allowed.\n  Please use the RandomNumberGeneratorService to set the random number seed.";
      }
      if( ! call_pygive(*itPar) ) {
        throw edm::Exception(edm::errors::Configuration,"PythiaError") 
          <<" pythia did not accept the following \""<<*itPar<<"\"";
      }
    }
    }else if(mySet == "CSAParameters"){   

   // Read CSA parameter
  
   pars = pythia_params.getParameter<vector<string> >("CSAParameters");

   cout << "----------------------------------------------" << endl; 
   cout << "Reading CSA parameter settings. " << endl;
   cout << "----------------------------------------------" << endl;                                                                           

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

}


// Read the TopReX parameters


  //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 Pythia random number seed " << endl;
    cout << "----------------------------------------------" << endl;
  edm::Service<RandomNumberGenerator> rng;
  randomEngine = fRandomEngine = &(rng->getEngine());
  uint32_t seed = rng->mySeed();
  ostringstream sRandomSet;
  sRandomSet <<"MRPY(1)="<<seed;
  call_pygive(sRandomSet.str());

 // srs insertion
  //  cout << "TopReX start called OK" << endl;
  call_toprex( "USER", "p", "p", 14000. );
  // cout << "TopReX was called OK" << endl;
  call_pyinit( "USER", "p", "p", 14000. );
  //
  //     call_pretauola(-1);     // TAUOLA initialization
  cout << endl; // Stetically add for the output

  produces<HepMCProduct>();
  cout << "ToprexSource: starting event generation ... " << endl;
}

ToprexSource::~ToprexSource

(

)

[virtual]

Destructor.

Definition at line 168 of file ToprexSource.cc.

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

                           {
  cout << "ToprexSource: event generation done. " << endl;
  call_pystat(1);
  //    call_pretauola(1);  // output from TAUOLA 
  clear(); 
}

---

Member Function Documentation

bool ToprexSource::call_pygive

(

const std::string &

iParm

)

[private]

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

Definition at line 230 of file ToprexSource.cc.

References pydat1, and PYGIVE.

Referenced by ToprexSource().

{
   int numWarn = pydat1.mstu[26]; //# warnings
   int numErr = pydat1.mstu[22];// # errors
//call the fortran routine pygive with a fortran string
  PYGIVE( iParm.c_str(), iParm.length() );  
//if an error or warning happens it is problem
   return pydat1.mstu[26] == numWarn && pydat1.mstu[22] == numErr;   
}

bool ToprexSource::call_txgive

(

const std::string &

iParm

)

[private]

Definition at line 242 of file ToprexSource.cc.

References TXGIVE.

Referenced by ToprexSource().

   {
   TXGIVE( iParm.c_str(), iParm.length() );  
        return 1;  
   }

bool ToprexSource::call_txgive_init

(

)

[private]

Definition at line 249 of file ToprexSource.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), and TXGIVE_INIT.

Referenced by ToprexSource().

{
   TXGIVE_INIT();
   cout << "  Setting CSA defaults.   "   << endl;
   return 1;
}

void ToprexSource::clear

(

void

)

[private]

Definition at line 175 of file ToprexSource.cc.

Referenced by ~ToprexSource().

{ }

bool ToprexSource::produce

(

Event &

e

)

[private, virtual]

Implements edm::ConfigurableInputSource.

Definition at line 178 of file ToprexSource.cc.

References call_findjets(), call_pylist(), conv, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), edm::ConfigurableInputSource::event(), evt, maxEventsToPrint_, edm::ConfigurableInputSource::numberEventsInRun(), edm::Event::put(), pypars, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, and edm::InputSource::remainingEvents().

                                    {

    auto_ptr<HepMCProduct> bare_product(new HepMCProduct());  
     //

    call_pyevnt();      // generate one event with Pythia

    //        call_pretauola(0);       // tau-lepton decays with TAUOLA 

    call_findjets();      // find jets 
    call_pyhepc( 1 );
    
    HepMC::GenEvent* evt = conv.read_next_event();
    evt->set_signal_process_id(pypars.msti[0]);
    evt->set_event_scale(pypars.pari[16]);
    evt->set_event_number(numberEventsInRun() - remainingEvents() - 1);

    //    Nest (*evt);




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

      // Prints PYLIST info
      if(pythiaPylistVerbosity_) {
        call_pylist(pythiaPylistVerbosity_);
      }
      
      // Prints HepMC event
      if(pythiaHepMCVerbosity_) {
        cout << "Event process = " << pypars.msti[0] << endl 
        << "----------------------" << endl;
        //      evt->print();
      }
    }
    
    
    //evt = reader_->fillCurrentEventData(); 
    //********                                      

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

    e.put(bare_product);

    return true;
}

---

Member Data Documentation

HepMC::GenEvent* edm::ToprexSource::evt [private]

Definition at line 56 of file ToprexSource.h.

Referenced by produce().

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

Definition at line 75 of file ToprexSource.h.

Referenced by ToprexSource().

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

Definition at line 76 of file ToprexSource.h.

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

Events to print if verbosity.

Definition at line 65 of file ToprexSource.h.

Referenced by produce(), and ToprexSource().

bool edm::ToprexSource::pythiaHepMCVerbosity_ [private]

HepMC verbosity flag.

Definition at line 63 of file ToprexSource.h.

Referenced by produce(), and ToprexSource().

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

Pythia PYLIST Verbosity flag.

Definition at line 61 of file ToprexSource.h.

Referenced by produce(), and ToprexSource().

---

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

• GeneratorInterface/TopRexInterface/interface/ToprexSource.h
• GeneratorInterface/TopRexInterface/src/ToprexSource.cc

---