---

edm::AlpgenProducer Class Reference

#include <GeneratorInterface/AlpgenInterface/interface/AlpgenProducer.h>

Inheritance diagram for edm::AlpgenProducer:

List of all members.

Public Member Functions
	AlpgenProducer (const ParameterSet &)
	Constructor.
virtual	~AlpgenProducer ()
	Destructor.
Private Member Functions
void	beginRun (Run &r)
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)
void	clear ()
virtual void	produce (Event &e, const EventSetup &es)
Private Attributes
bool	doubleParticle
double	etamax
double	etamin
unsigned int	eventsRead_
HepMC::GenEvent *	evt
std::string	fileName_
std::vector< std::string >	fileNames_
CLHEP::HepRandomEngine *	fRandomEngine
CLHEP::RandFlat *	fRandomGenerator
LHERunInfoProduct::Header	lheAlpgenUnwParHeader
unsigned int	maxEventsToPrint_
	Events to print if verbosity.
int	Nev_
int	particleID
double	phimax
double	phimin
double	ptmax
double	ptmin
bool	pythiaHepMCVerbosity_
	HepMC verbosity flag.
unsigned int	pythiaPylistVerbosity_
	Pythia PYLIST Verbosity flag.

---

Detailed Description

Definition at line 36 of file AlpgenProducer.h.

---

Constructor & Destructor Documentation

AlpgenProducer::AlpgenProducer

(

const ParameterSet &

pset

)

Constructor.

Definition at line 51 of file AlpgenProducer.cc.

References LHERunInfoProduct::Header::addLine(), call_pygive(), call_txgive(), edm::errors::Configuration, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), fileName_, fileNames_, fRandomEngine, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), lheAlpgenUnwParHeader, GenMuonPlsPt100GeV_cfg::maxEvents, maxEventsToPrint_, Nev_, pars, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, and randomEngine.

                                                         :
  EDProducer(), evt(0), 
  pythiaPylistVerbosity_ (pset.getUntrackedParameter<int>("pythiaPylistVerbosity",0)),
  pythiaHepMCVerbosity_ (pset.getUntrackedParameter<bool>("pythiaHepMCVerbosity",false)),
  maxEventsToPrint_ (pset.getUntrackedParameter<int>("maxEventsToPrint",1)),
// JMM experimenting
  fileNames_ (pset.getUntrackedParameter<std::vector<std::string> >("fileNames")),
  eventsRead_(0),
  lheAlpgenUnwParHeader("AlpgenUnwParFile")  
// end JMM experimenting
{
  
  fileName_ = fileNames_[0];
  // strip the file: 
  if ( fileName_.find("file:") || fileName_.find("rfio:")){
    fileName_.erase(0,5);
  }   

  // open the .unw file to store additional 
  // informations in the AlpgenInfoProduct
//  unwfile = new ifstream((fileName_+".unw").c_str());
  // get the number of input events from  _unw.par files
  char buffer[256];
  ifstream reader((fileName_+"_unw.par").c_str());
  char sNev[80];
  lheAlpgenUnwParHeader.addLine("\n");
  while ( reader.getline (buffer,256) ) {
    istringstream is(buffer);
    lheAlpgenUnwParHeader.addLine(std::string(buffer) + "\n");
    is >> sNev;
    Nev_ = atoi(sNev);
  }

// JMM experimenting
#ifdef NEVER
  //check that N(asked events) <= N(input events)
  if(maxEvents()>Nev_) {
    cout << "ALPGEN warning: Number of events requested > Number of unweighted events." << endl;
    cout << "                Execution will stop after processing the last unweighted event" << endl;  
  }

  if(maxEvents() != -1 && maxEvents() < Nev_) // stop at N(asked events) if N(asked events)<N(input events)
    Nev_ = maxEvents();
#endif
// end JMM experimenting
  
  // PYLIST Verbosity Level
  // Valid PYLIST arguments are: 1, 2, 3, 5, 7, 11, 12, 13
  pythiaPylistVerbosity_ = pset.getUntrackedParameter<int>("pythiaPylistVerbosity",0);
  
  // HepMC event verbosity Level
  pythiaHepMCVerbosity_ = pset.getUntrackedParameter<bool>("pythiaHepMCVerbosity",false);
  
  //Max number of events printed on verbosity level 
  maxEventsToPrint_ = pset.getUntrackedParameter<int>("maxEventsToPrint",0);
  
  // Set PYTHIA parameters in a single ParameterSet
  {
    ParameterSet pythia_params = 
      pset.getParameter<ParameterSet>("PythiaParameters") ;
    
    // Read the PYTHIA parameters for each set of parameters
    vector<string> pars = 
      pythia_params.getParameter<vector<string> >("pythia");
    
    // 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 seed.";
      }
      if( ! call_pygive(*itPar) ) {
        throw edm::Exception(edm::errors::Configuration,"PythiaError") 
          <<" pythia did not accept the following \""<<*itPar<<"\"";
      }
    }
  }

  // Read the Alpgen parameters

  // read External Generator parameters
  {   ParameterSet generator_params = 
        pset.getParameter<ParameterSet>("GeneratorParameters") ;
  vector<string> pars = 
    generator_params.getParameter<vector<string> >("generator");
  for( vector<string>::const_iterator  
         itPar = pars.begin(); itPar != pars.end(); ++itPar ) 
    {
      call_txgive(*itPar);          
    }
  // giving to txgive a string with the alpgen filename
  string tmpstring = "UNWFILE = " + fileName_;
  call_txgive(tmpstring);
  }
  
  
  //In the future, we will get the random number seed on each event and tell 
  // pythia to use that new seed
  edm::Service<RandomNumberGenerator> rng;
  randomEngine = fRandomEngine = &(rng->getEngine());
  uint32_t seed = rng->mySeed();
  ostringstream sRandomSet;
  sRandomSet <<"MRPY(1)="<<seed;
  call_pygive(sRandomSet.str());
  
  //  call_pretauola(-1);     // TAUOLA initialization
  call_pyinit( "USER", "p", "p", 14000. );
  
  cout << endl; // Stetically add for the output
  //********                                      
  
  produces<HepMCProduct>();
  produces<LHEEventProduct>();

//  produces<AlpWgtFileInfoProduct, edm::InRun>();
  produces<LHERunInfoProduct, edm::InRun>();
}

AlpgenProducer::~AlpgenProducer

(

)

[virtual]

Destructor.

Definition at line 172 of file AlpgenProducer.cc.

References alpgen_end, and clear().

                               {
  call_pystat(1);
  //  call_pretauola(1);  // output from TAUOLA 
  alpgen_end();
  clear(); 
}

---

Member Function Documentation

void AlpgenProducer::beginRun

(

Run &

r

)

[private]

Definition at line 183 of file AlpgenProducer.cc.

References LHERunInfoProduct::Header::addLine(), indexGen::comments, fileName_, lheAlpgenUnwParHeader, edm::Run::put(), lhef::CommonBlocks::readHEPRUP(), and edm::runInfo.

                                     {
  // get the LHE init information from Fortran code
  lhef::HEPRUP heprup;
  lhef::CommonBlocks::readHEPRUP(&heprup);
  auto_ptr<LHERunInfoProduct> runInfo(new LHERunInfoProduct(heprup));

  // information on weighted events
  LHERunInfoProduct::Header lheAlpgenWgtHeader("AlpgenWgtFile");
  lheAlpgenWgtHeader.addLine("\n");
  ifstream wgtascii((fileName_+".wgt").c_str());
  char buffer[512];
  while(wgtascii.getline(buffer,512)) {
    lheAlpgenWgtHeader.addLine(std::string(buffer) + "\n");
  }

  // comments on top
  LHERunInfoProduct::Header comments;
  comments.addLine("\n");
  comments.addLine("\tExtracted by AlpgenInterface\n");

  // build the final Run info object
  runInfo->addHeader(comments);
  runInfo->addHeader(lheAlpgenUnwParHeader);
  runInfo->addHeader(lheAlpgenWgtHeader);
  r.put(runInfo);
}

bool AlpgenProducer::call_pygive

(

const std::string &

iParm

)

[private]

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

Definition at line 321 of file AlpgenProducer.cc.

References pydat1, and PYGIVE.

Referenced by AlpgenProducer().

                                                   {

  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() );  
  //  PYGIVE( iParm );  
//if an error or warning happens it is problem
  return pydat1.mstu[26] == numWarn && pydat1.mstu[22] == numErr;   
}

bool AlpgenProducer::call_txgive

(

const std::string &

iParm

)

[private]

Definition at line 334 of file AlpgenProducer.cc.

References TXGIVE.

Referenced by AlpgenProducer().

   {
    //call the fortran routine txgive with a fortran string
    TXGIVE( iParm.c_str(), iParm.length() );  
    return 1;  
   }

void AlpgenProducer::clear

(

void

)

[private]

Definition at line 179 of file AlpgenProducer.cc.

Referenced by ~AlpgenProducer().

                           {
  
}

void AlpgenProducer::produce	(	Event &	e,
		const EventSetup &	es
	)			[private, virtual]

Implements edm::EDProducer.

Definition at line 210 of file AlpgenProducer.cc.

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, call_pylist(), conv2, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), edm::EventID::event(), eventsRead_, evt, Exception, i, edm::Event::id(), maxEventsToPrint_, Nev_, lhef::HEPEUP::NUP, edm::Event::put(), pyint1, pypars, pythiaHepMCVerbosity_, pythiaPylistVerbosity_, lhef::CommonBlocks::readHEPEUP(), and lhef::HEPEUP::SPINUP.

                                                            {
  
  // exit if N(events asked) has been exceeded
  if(e.id().event()> Nev_) {
    throw cms::Exception("Generator") << "Can't produce event because _unw.par file is over.";
  } else {
    
    auto_ptr<HepMCProduct> bare_product(new HepMCProduct());  
    
//    // Additional information from unweighted file
//    auto_ptr<AlpgenInfoProduct> alp_product(new AlpgenInfoProduct());

    // Extract from .unw file the info for AlpgenInfoProduct
    
//    char buffer[512];
//    if(unwfile->getline(buffer,512)) {
//      alp_product->EventInfo(buffer);
//    }
//    if(unwfile->getline(buffer,512)) 
//      alp_product->InPartonInfo(buffer);
//    if(unwfile->getline(buffer,512)) 
//      alp_product->InPartonInfo(buffer);
//    for(int i_out = 0; i_out <  alp_product->nTot()-2; i_out++) {
//      if(unwfile->getline(buffer,512)) 
//      alp_product->OutPartonInfo(buffer);
//    }

    call_pyevnt();      // generate one event with Pythia
    //        call_pretauola(0);  // tau-lepton decays with TAUOLA 

    // fill the parton level LHE event information
    lhef::HEPEUP hepeup;
    lhef::CommonBlocks::readHEPEUP(&hepeup);
    hepeup.AQEDUP = hepeup.AQCDUP = -1.0; // alphas are not saved by Alpgen
    for(int i = 0; i < hepeup.NUP; i++)
      hepeup.SPINUP[i] = -9;    // Alpgen does not store spin information
    auto_ptr<LHEEventProduct> lheEvent(new LHEEventProduct(hepeup));
    // Moved to later, first we have to make sure that the event generated by Pythia is OK.
    // e.put(lheEvent);

    // The way this is implemented, call_pyevnt() will keep reading events from the .unw until
    // one of them is able to pass the matching veto. If it can't read a last event, it will
    // simply continue with whatever event was in memory. In this case, there are two possibilities:
    // A) The event passed the matching veto. There is meaningful information in HEPEUP. Continue.
    // (Will fail next Event with HEPEUP.NUP == 0. Harmless.)
    // B) The event didn't pass. HEPEUP.NUP == 0. This should never happen, and we should abort here.
    if(hepeup.NUP==0) {
      edm::LogInfo("Generator|TooLittleData") << "ALPGEN warning: last unweighted event reached.\n"
                                              << "                (hepeup.NUP == 0)\n"
                                              << "                The event number " << e.id().event() << " will not be written to disk.";  
      throw cms::Exception("Generator") << "Can't produce event because _unw.par file is over.";
    }
    
    call_pyhepc( 1 );
    //    HepMC::GenEvent* evt = conv2.getGenEventfromHEPEVT();
    HepMC::GenEvent* evt = conv2.read_next_event();
    
    evt->set_signal_process_id(pypars.msti[0]);
    ++eventsRead_;
    evt->set_event_number(eventsRead_);
    
    int id1 = pyint1.mint[14];
    int id2 = pyint1.mint[15];
    if ( id1 == 21 ) id1 = 0;
    if ( id2 == 21 ) id2 = 0; 
    double x1 = pyint1.vint[40];
    double x2 = pyint1.vint[41];  
    double Q  = pyint1.vint[50];
    double pdf1 = pyint1.vint[38];
    pdf1 /= x1 ;
    double pdf2 = pyint1.vint[39];
    pdf2 /= x2 ;
    evt->set_pdf_info( HepMC::PdfInfo(id1,id2,x1,x2,Q,pdf1,pdf2) ) ;
   
    evt->weights().push_back( pyint1.vint[96] );

    //******** Verbosity ********
    
    if(e.id().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();
      }
    }

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

    // Last check to make sure there are no empty events. This DOES NOT
    // supersedes the AlpgenEmptyEventFilter, because this is a Producer,
    // not a source. So, it can only return, instead of returning false.
    // Should we take extra care here?

    if(!(evt->particles_empty())) {
      e.put(lheEvent);
      e.put(bare_product);
    }
    
    return;
  }
}

---

Member Data Documentation

bool edm::AlpgenProducer::doubleParticle [private]

Definition at line 75 of file AlpgenProducer.h.

double edm::AlpgenProducer::etamax [private]

Definition at line 77 of file AlpgenProducer.h.

double edm::AlpgenProducer::etamin [private]

Definition at line 77 of file AlpgenProducer.h.

unsigned int edm::AlpgenProducer::eventsRead_ [private]

Definition at line 71 of file AlpgenProducer.h.

Referenced by produce().

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

Definition at line 60 of file AlpgenProducer.h.

Referenced by produce().

std::string edm::AlpgenProducer::fileName_ [private]

Definition at line 70 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and beginRun().

std::vector<std::string> edm::AlpgenProducer::fileNames_ [private]

Definition at line 69 of file AlpgenProducer.h.

Referenced by AlpgenProducer().

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

Definition at line 80 of file AlpgenProducer.h.

Referenced by AlpgenProducer().

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

Definition at line 81 of file AlpgenProducer.h.

LHERunInfoProduct::Header edm::AlpgenProducer::lheAlpgenUnwParHeader [private]

Definition at line 83 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and beginRun().

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

Events to print if verbosity.

Definition at line 67 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and produce().

int edm::AlpgenProducer::Nev_ [private]

Definition at line 50 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and produce().

int edm::AlpgenProducer::particleID [private]

Definition at line 74 of file AlpgenProducer.h.

double edm::AlpgenProducer::phimax [private]

Definition at line 78 of file AlpgenProducer.h.

double edm::AlpgenProducer::phimin [private]

Definition at line 78 of file AlpgenProducer.h.

double edm::AlpgenProducer::ptmax [private]

Definition at line 76 of file AlpgenProducer.h.

double edm::AlpgenProducer::ptmin [private]

Definition at line 76 of file AlpgenProducer.h.

bool edm::AlpgenProducer::pythiaHepMCVerbosity_ [private]

HepMC verbosity flag.

Definition at line 65 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and produce().

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

Pythia PYLIST Verbosity flag.

Definition at line 63 of file AlpgenProducer.h.

Referenced by AlpgenProducer(), and produce().

---

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

• GeneratorInterface/AlpgenInterface/interface/AlpgenProducer.h
• GeneratorInterface/AlpgenInterface/src/AlpgenProducer.cc

---