LHAupLesHouches Class Reference

#include <LHAupLesHouches.h>

Inheritance diagram for LHAupLesHouches:

Public Member Functions
	LHAupLesHouches ()
void	loadEvent (lhef::LHEEvent *event)
void	loadRunInfo (lhef::LHERunInfo *runInfo)
void	setScalesFromLHEF (bool b)
	~LHAupLesHouches () override

Private Member Functions
bool	setEvent (int idProcIn) override
bool	setInit () override

Private Attributes
lhef::LHEEvent *	event
std::map< std::string, std::string > *	fEvAttributes
lhef::LHERunInfo *	runInfo
bool	setScalesFromLHEF_

Detailed Description

Definition at line 20 of file LHAupLesHouches.h.

Constructor & Destructor Documentation

LHAupLesHouches::LHAupLesHouches ( )

inline

Definition at line 22 of file LHAupLesHouches.h.

: setScalesFromLHEF_(false), fEvAttributes(nullptr) { ; }

LHAupLesHouches::~LHAupLesHouches ( )

inlineoverride

Definition at line 32 of file LHAupLesHouches.h.

References fEvAttributes, setEvent(), and setInit().

                              {
    if (fEvAttributes)
      delete fEvAttributes;
  }

Member Function Documentation

void LHAupLesHouches::loadEvent ( lhef::LHEEvent *  event )

inline

Definition at line 28 of file LHAupLesHouches.h.

References event.

{ this->event = event; }

void LHAupLesHouches::loadRunInfo ( lhef::LHERunInfo *  runInfo )

inline

Definition at line 25 of file LHAupLesHouches.h.

References runInfo.

{ this->runInfo = runInfo; }

bool LHAupLesHouches::setEvent ( int  idProcIn )

overrideprivate

Definition at line 50 of file LHAupLesHouches.cc.

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, edmScanValgrind::buffer, indexGen::comments, mps_fire::i, lhef::HEPEUP::ICOLUP, gen::PdfInfo::id, lhef::HEPEUP::IDPRUP, lhef::HEPEUP::IDUP, lhef::HEPEUP::ISTUP, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, lhef::HEPEUP::PUP, gen::PdfInfo::scalePDF, lhef::HEPEUP::SCALUP, lhef::HEPEUP::SPINUP, lhef::HEPEUP::VTIMUP, gen::PdfInfo::x, gen::PdfInfo::xPDF, and lhef::HEPEUP::XWGTUP.

Referenced by ~LHAupLesHouches().

                                           {
  if (!event)
    return false;

  if (event->getReadAttempts() > 0)
    return false;  // record already tried

  const lhef::HEPEUP &hepeup = *event->getHEPEUP();

  if (!hepeup.NUP)
    return false;

  setProcess(hepeup.IDPRUP, hepeup.XWGTUP, hepeup.SCALUP, hepeup.AQEDUP, hepeup.AQCDUP);

  const std::vector<float> &scales = event->scales();

  unsigned int iscale = 0;
  for (int i = 0; i < hepeup.NUP; i++) {
    //retrieve scale corresponding to each particle
    double scalein = -1.;

    //handle clustering scales if present,
    //applies to outgoing partons only
    if (setScalesFromLHEF_ && !scales.empty() && hepeup.ISTUP[i] == 1) {
      if (iscale >= scales.size()) {
        edm::LogError("InvalidLHEInput") << "Pythia8 requires"
                                         << "cluster scales for all outgoing partons or for none" << std::endl;
      }
      scalein = scales[iscale];
      ++iscale;
    }

    addParticle(hepeup.IDUP[i],
                hepeup.ISTUP[i],
                hepeup.MOTHUP[i].first,
                hepeup.MOTHUP[i].second,
                hepeup.ICOLUP[i].first,
                hepeup.ICOLUP[i].second,
                hepeup.PUP[i][0],
                hepeup.PUP[i][1],
                hepeup.PUP[i][2],
                hepeup.PUP[i][3],
                hepeup.PUP[i][4],
                hepeup.VTIMUP[i],
                hepeup.SPINUP[i],
                scalein);
  }

  if (!infoPtr->eventAttributes) {
    fEvAttributes->clear();
    infoPtr->eventAttributes = fEvAttributes;
  } else {
    infoPtr->eventAttributes = fEvAttributes;  // make sure still there
    infoPtr->eventAttributes->clear();
  }

  //fill parton multiplicities if available
  int npLO = event->npLO();
  int npNLO = event->npNLO();

  //default value of -99 indicates tags were not present in the original LHE file
  //don't pass to pythia in this case to emulate pythia internal lhe reader behaviour
  if (npLO != -99) {
    char buffer[100];
    snprintf(buffer, 100, "%i", npLO);
    (*infoPtr->eventAttributes)["npLO"] = buffer;
  }
  if (npNLO != -99) {
    char buffer[100];
    snprintf(buffer, 100, "%i", npNLO);
    (*infoPtr->eventAttributes)["npNLO"] = buffer;
  }

  //add #rwgt info from comments
  const std::vector<std::string> &comments = event->getComments();
  for (unsigned i = 0; i < comments.size(); i++) {
    if (comments[i].rfind("#rwgt", 0) == 0)
      (*infoPtr->eventAttributes)["#rwgt"] = comments[i];
  }

  const lhef::LHEEvent::PDF *pdf = event->getPDF();
  if (pdf) {
    this->setPdf(pdf->id.first,
                 pdf->id.second,
                 pdf->x.first,
                 pdf->x.second,
                 pdf->scalePDF,
                 pdf->xPDF.first,
                 pdf->xPDF.second,
                 true);
  } else {
    this->setPdf(hepeup.IDUP[0],
                 hepeup.IDUP[1],
                 hepeup.PUP[0][3] / runInfo->getHEPRUP()->EBMUP.first,
                 hepeup.PUP[1][3] / runInfo->getHEPRUP()->EBMUP.second,
                 0.,
                 0.,
                 0.,
                 false);
  }

  event->attempted();

  return true;
}

bool LHAupLesHouches::setInit ( )

overrideprivate

Definition at line 14 of file LHAupLesHouches.cc.

References lhef::HEPRUP::EBMUP, mps_fire::i, lhef::HEPRUP::IDBMUP, lhef::HEPRUP::IDWTUP, lhef::HEPRUP::LPRUP, lhef::HEPRUP::NPRUP, lhef::HEPRUP::PDFGUP, lhef::HEPRUP::PDFSUP, AlCaHLTBitMon_QueryRunRegistry::string, lhef::HEPRUP::XERRUP, lhef::HEPRUP::XMAXUP, and lhef::HEPRUP::XSECUP.

Referenced by ~LHAupLesHouches().

                              {
  if (!runInfo)
    return false;
  const lhef::HEPRUP &heprup = *runInfo->getHEPRUP();

  setBeamA(heprup.IDBMUP.first, heprup.EBMUP.first, heprup.PDFGUP.first, heprup.PDFSUP.first);
  setBeamB(heprup.IDBMUP.second, heprup.EBMUP.second, heprup.PDFGUP.second, heprup.PDFSUP.second);
  setStrategy(heprup.IDWTUP);

  for (int i = 0; i < heprup.NPRUP; i++)
    addProcess(heprup.LPRUP[i], heprup.XSECUP[i], heprup.XERRUP[i], heprup.XMAXUP[i]);

  //hadronisation->onInit().emit();

  //runInfo.reset();

  //fill SLHA header information if available
  std::vector<std::string> slha = runInfo->findHeader("slha");
  if (!slha.empty()) {
    std::string slhaheader;
    for (std::vector<std::string>::const_iterator iter = slha.begin(); iter != slha.end(); ++iter) {
      slhaheader.append(*iter);
    }
    infoPtr->setHeader("slha", slhaheader);
  }

  //will be used to work around missing initialization inside pythia8
  if (!fEvAttributes) {
    fEvAttributes = new std::map<std::string, std::string>;
  } else {
    fEvAttributes->clear();
  }

  return true;
}

void LHAupLesHouches::setScalesFromLHEF ( bool  b )

inline

Definition at line 30 of file LHAupLesHouches.h.

References b, and setScalesFromLHEF_.

{ setScalesFromLHEF_ = b; }

Member Data Documentation

lhef::LHEEvent* LHAupLesHouches::event

private

Definition at line 44 of file LHAupLesHouches.h.

Referenced by Types.EventID::cppID(), loadEvent(), looper.Looper::process(), and core.AutoHandle.AutoHandle::product().

std::map<std::string, std::string>* LHAupLesHouches::fEvAttributes

private

Definition at line 49 of file LHAupLesHouches.h.

Referenced by ~LHAupLesHouches().

lhef::LHERunInfo* LHAupLesHouches::runInfo

private

Definition at line 42 of file LHAupLesHouches.h.

Referenced by loadRunInfo().

bool LHAupLesHouches::setScalesFromLHEF_

private

Definition at line 47 of file LHAupLesHouches.h.

Referenced by setScalesFromLHEF().