ThePEG::powhegLesHouchesFileReader Class Reference

powhegLesHouchesFileReader is an base class to be used for objects which reads event files from matrix element generators. More...

#include <PowhegLesHouchesFileReader.h>

Inheritance diagram for ThePEG::powhegLesHouchesFileReader:

Collaboration diagram for ThePEG::powhegLesHouchesFileReader:

Public Member Functions
Standard constructors and destructors.
	powhegLesHouchesFileReader ()
	Default constructor. More...
	powhegLesHouchesFileReader (const powhegLesHouchesFileReader &)
	Copy-constructor. More...
virtual	~powhegLesHouchesFileReader ()
	Destructor. More...
Virtual functions specified by the LesHouchesReader base class.
virtual void	initialize (LesHouchesEventHandler &eh)
	Initialize. More...
virtual void	open ()
	Open a file with events. More...
virtual void	close ()
	Close the file from which events have been read. More...
virtual bool	doReadEvent ()
	Read the next event from the file or stream into the corresponding protected variables. More...
string	filename () const
	Return the name of the file from where to read events. More...
virtual vector< string >	optWeightsNamesFunc ()
	Return the optional weights information string ("Names") More...

Protected Member Functions
Clone Methods.
virtual IBPtr	clone () const
	Make a simple clone of this object. More...
virtual IBPtr	fullclone () const
	Make a clone of this object, possibly modifying the cloned object to make it sane. More...

Functions used by the persistent I/O system.
void	persistentOutput (PersistentOStream &os) const
	Function used to write out object persistently. More...
void	persistentInput (PersistentIStream &is, int version)
	Function used to read in object persistently. More...
void	erase_substr (std::string &subject, const std::string &search)
	Erases all occurences of a substring from a string. More...
static void Init	ATLAS_NOT_THREAD_SAFE ()
	Standard Init function used to initialize the interfaces. More...

Standard (and non-standard) Interfaced functions.
CFileLineReader	m_cfile
	The wrapper around the C FILE stream from which to read. More...
long	m_neve
	The number of events in this file. More...
long	m_ieve
	The current event number. More...
string	m_LHFVersion
	If the file is a standard Les Houches formatted file (LHF) this is its version number. More...
string	m_outsideBlock
	If LHF. More...
string	m_headerBlock
	If LHF. More...
string	m_initComments
	If LHF. More...
map< string, string >	m_initAttributes
	If LHF. More...
string	m_eventComments
	If LHF. More...
map< string, string >	m_eventAttributes
	If LHF. More...
string	m_theFileName
	The name of the file from where to read events. More...
bool	m_theQNumbers
	Whether or not to search for QNUMBERS stuff. More...
bool	m_theIncludeFxFxTags
	Include/Read FxFx tags. More...
bool	m_theIncludeCentral
	Include central weight (for backup use) More...
DecayerPtr	m_theDecayer
	Decayer for any decay modes read from the file. More...
map< string, string >	m_scalemap
	Further information on the weights. More...
map< string, double >	m_optionalWeightsTemp
	Temporary holder for optional weights. More...
map< string, string >	m_optionalWeightsLabel
static const ClassDescription< powhegLesHouchesFileReader >	m_initpowhegLesHouchesFileReader
	Describe an abstract base class with persistent data. More...
virtual void	doinit ()
	Initialize this object after the setup phase before saving an EventGenerator to disk. More...
virtual bool	preInitialize () const
	Return true if this object needs to be initialized before all other objects because it needs to extract PDFs from the event file. More...
powhegLesHouchesFileReader &	operator= (const powhegLesHouchesFileReader &)
	Private and non-existent assignment operator. More...

Detailed Description

powhegLesHouchesFileReader is an base class to be used for objects which reads event files from matrix element generators.

It inherits from LesHouchesReader and extends it by defining a file handle to be read from, which is opened and closed by the open() and close() functions. Note that the file handle is a standard C filehandle and not a C++ stream. This is because there is no standard way in C++ to connect a pipe to a stream for reading eg. gzipped files. This class is able to read plain event files conforming to the Les Houches Event File accord.

See also

Th1e interfaces defined for powhegLesHouchesFileReader.

Event

LesHouchesReader

Definition at line 43 of file PowhegLesHouchesFileReader.h.

Constructor & Destructor Documentation

powhegLesHouchesFileReader() [1/2]

ThePEG::powhegLesHouchesFileReader::powhegLesHouchesFileReader ( )

inline

Default constructor.

Definition at line 52 of file PowhegLesHouchesFileReader.h.

                               : m_neve(0), m_ieve(0), m_theQNumbers(false),
               m_theIncludeFxFxTags(false),
               m_theIncludeCentral(false) {}

powhegLesHouchesFileReader() [2/2]

powhegLesHouchesFileReader::powhegLesHouchesFileReader

(

const powhegLesHouchesFileReader &

x

)

Copy-constructor.

Note that a file which is opened in the object copied from will have to be reopened in this.

Definition at line 42 of file PowhegLesHouchesFileReader.cxx.

  : LesHouchesReader(x), m_neve(x.m_neve), m_ieve(0),
    m_LHFVersion(x.m_LHFVersion), m_outsideBlock(x.m_outsideBlock),
    m_headerBlock(x.m_headerBlock), m_initComments(x.m_initComments),
    m_initAttributes(x.m_initAttributes), m_eventComments(x.m_eventComments),
    m_eventAttributes(x.m_eventAttributes),
    m_theFileName(x.m_theFileName), m_theQNumbers(x.m_theQNumbers),
    m_theIncludeFxFxTags(x.m_theIncludeFxFxTags),
    m_theIncludeCentral(x.m_theIncludeCentral),
    m_theDecayer(x.m_theDecayer) {}

~powhegLesHouchesFileReader()

powhegLesHouchesFileReader::~powhegLesHouchesFileReader ( )

virtual

Destructor.

Definition at line 54 of file PowhegLesHouchesFileReader.cxx.

{}

Member Function Documentation

ATLAS_NOT_THREAD_SAFE()

static void Init ThePEG::powhegLesHouchesFileReader::ATLAS_NOT_THREAD_SAFE ( )

static

Standard Init function used to initialize the interfaces.

clone()

IBPtr powhegLesHouchesFileReader::clone ( ) const

protectedvirtual

Make a simple clone of this object.

Returns

a pointer to the new object.

Definition at line 56 of file PowhegLesHouchesFileReader.cxx.

                                              {
  return new_ptr(*this);
}

close()

void powhegLesHouchesFileReader::close ( )

virtual

Close the file from which events have been read.

Definition at line 815 of file PowhegLesHouchesFileReader.cxx.

                                       {
  m_cfile.close();
}

doinit()

void powhegLesHouchesFileReader::doinit ( )

protectedvirtual

Initialize this object after the setup phase before saving an EventGenerator to disk.

Exceptions

InitException

if object could not be initialized properly.

Definition at line 68 of file PowhegLesHouchesFileReader.cxx.

                                        {
  LesHouchesReader::doinit();
 
  // are we using QNUMBERS
  if(!m_theQNumbers) return;
  // parse the header block and create 
  // any new particles needed in QNUMBERS blocks
  string block = m_headerBlock;
  string line  = "";
  bool readingSLHA = false;
  int (*pf)(int) = tolower;
  unsigned int newNumber(0);
  do {
    line  = StringUtils::car(block,"\r\n");
    block = StringUtils::cdr(block,"\r\n");
    if(line[0]=='#') continue;
    // are we reading the SLHA block
    if(readingSLHA) {
      // reached the end of slha block 
      if(line.find("</slha") != string::npos) {
    readingSLHA = false;
    break;
      }
      // remove trailing comment from line
      vector<string> split = StringUtils::split(line,"#");
      // check for a qnumbers block
      transform(split[0].begin(), split[0].end(), split[0].begin(), pf);
      // if not contine
      if(split[0].find("block qnumbers")==string::npos)
    continue;
      // get name from comment
      string name;
      if(split.size()>=2) {
    name = StringUtils::stripws(split[1]);
      }
      else {
    ++newNumber;
    ostringstream tname;
    tname << "NP" << newNumber;
    name = tname.str();
      }
      // extract the PDG code
      split = StringUtils::split(split[0]," ");
      istringstream is(split[2]);
      long PDGCode(0);
      is >> PDGCode;
      // get the charge, spin, colour and whether an antiparticle
      int charge(0),spin(0),colour(0),anti(0);
      for(unsigned int ix=0;ix<4;++ix) {
    line = StringUtils::car(block,"\r\n");
    block = StringUtils::cdr(block,"\r\n");
    int dummy[2];
    istringstream is(line);
    is >> dummy[0] >> dummy[1];
    switch (dummy[0]) {
    case 1:
      charge = dummy[1];
      break;
    case 2:
      spin   = dummy[1];
      break;
    case 3:
      colour = dummy[1];
      break;
    case 4:
      anti   = dummy[1];
      break;
    default:
      assert(false);
    }
      }
      // check if particles already exist
      PDPair newParticle;
      newParticle.first = getParticleData(PDGCode);
      if(newParticle.first) Throw<SetupException>()
                  << "Particle with PDG code " << PDGCode 
                  << " whose creation was requested in a QNUMBERS Block"
                  << " already exists. Retaining the original particle" 
                  << Exception::warning;
      if(anti) {
    newParticle.second = getParticleData(-PDGCode);
    if(newParticle.second) Throw<SetupException>()
                 << "Anti-particle with PDG code " << -PDGCode 
                 << " whose creation was requested in a QNUMBERS Block"
                 << " already exists. Retaining the original particle" 
                 << Exception::warning;
    if(( newParticle.first  && !newParticle.second ) ||
       ( newParticle.second && !newParticle.first  ) )
      Throw<SetupException>()
        << "Either particle or anti-particle with PDG code " << PDGCode 
        << " whose creation was requested in a QNUMBERS Block"
        << " already exists, but not both the particle and antiparticle. "
        << " Something dodgy here stopping" 
        << Exception::runerror;
      }
      // already exists continue
      if(newParticle.first) continue;
      // create the particles
      // particle with no anti particle
      if( anti == 0 ) {
    // construct the name
    if(name=="") {
      ostringstream temp;
      temp << PDGCode;
      name = temp.str(); 
    }
    // create the ParticleData object
    newParticle.first = ParticleData::Create(PDGCode,name);
      }
      // particle anti-particle pair
      else {
    // construct the names
    string nameAnti;
    if(name=="") {
      ostringstream temp;
      temp << PDGCode;
      name = temp.str(); 
      ostringstream temp2;
      temp << -PDGCode;
      nameAnti = temp2.str(); 
    }
    else {
      nameAnti=name;
          auto swapSign=[](char & t){
            if (t=='-') t= '+';
            else if (t=='+') t= '-';
          };
          std::for_each(nameAnti.begin(),nameAnti.end(), swapSign);
      if(nameAnti==name) nameAnti += "bar";
    }
    // create the ParticleData objects
    newParticle = ParticleData::Create(PDGCode,name,nameAnti);
      }
      // set the particle properties
      if(colour==1) colour = 0;
      newParticle.first->iColour(PDT::Colour(colour));
      newParticle.first->iSpin  (PDT::Spin  (spin  ));
      newParticle.first->iCharge(PDT::Charge(charge));
      // register it
      generator()->preinitRegister(newParticle.first,
                   "/Herwig/Particles/"+newParticle.first->PDGName());
      // set the antiparticle properties
      if(newParticle.second) {
    if(colour==3||colour==6) colour *= -1;
    charge = -charge;
    newParticle.second->iColour(PDT::Colour(colour));
    newParticle.second->iSpin  (PDT::Spin  (spin  ));
    newParticle.second->iCharge(PDT::Charge(charge));
    // register it
    generator()->preinitRegister(newParticle.second,
                     "/Herwig/Particles/"+newParticle.second->PDGName());
      }
    }
    // start of SLHA block ?
    else if(line.find("<slha") != string::npos) {
      readingSLHA = true;
    }
  }
  while(line!="");
  // now set any masses/decay modes
  block = m_headerBlock;
  line="";
  readingSLHA=false;
  bool ok=true;
  do {
    line = StringUtils::car(block,"\r\n");
    block = StringUtils::cdr(block,"\r\n");
    // are we reading the SLHA block
    if(readingSLHA) {
      // reached the end?
      if(line.find("</slha") == 0 ) {
    readingSLHA = false;
    break;
      }
      // make lower case
      transform(line.begin(),line.end(),line.begin(), pf);
      // found the mass block ?
      if(line.find("block mass")!=string::npos) {
    // read it
    line = StringUtils::car(block,"\r\n");
    // check not at end
    while(line[0] != 'D' && line[0] != 'B' &&
          line[0] != 'd' && line[0] != 'b' &&
          line    != "") {
      // skip comment lines
      if(line[0] == '#') {
        block = StringUtils::cdr(block,"\r\n");
        line = StringUtils::car(block,"\r\n");
        continue;
      }
      // get the mass and PGD code
      istringstream temp(line);
      long id;
      double mass;
      temp >> id >> mass;
      // skip resetting masses on SM particles
      // as it can cause problems later on in event generation
      if(std::abs(id) <= 6 || (std::abs(id) >= 11 && std::abs(id) <= 16) ||
         std::abs(id) == 23 || std::abs(id) == 24) {
        block = StringUtils::cdr(block,"\r\n");
        line = StringUtils::car(block,"\r\n");
        continue;
      }
      // magnitude of mass for susy models
      mass = abs(mass);
      // set the mass
      tPDPtr particle = getParticleData(id);
      if(!particle) throw SetupException() 
              << "powhegLesHouchesFileReader::doinit() - Particle with PDG code not"
              << id << " not found." << Exception::runerror;
      const InterfaceBase * ifb = BaseRepository::FindInterface(particle,
                                    "NominalMass");
      ostringstream os;
      os << mass;
      ifb->exec(*particle, "set", os.str());
      // read the next line
      block = StringUtils::cdr(block,"\r\n");
      line = StringUtils::car(block,"\r\n");
    }; 
      }
      // found a decay block
      else if(line.find("decay") == 0) {
    // get PGD code and width
    istringstream iss(line);
    string dummy;
    long parent(0);
    Energy width(ZERO);
    iss >> dummy >> parent >> iunit(width, GeV);
    // get the ParticleData object
    PDPtr inpart = getParticleData(parent);
    if(!inpart)  {
      throw SetupException() 
        << "powhegLesHouchesFileReader::doinit() - A ParticleData object with the PDG code "
        << parent << " does not exist. " << Exception::runerror;
      return;
    }
    if ( std::abs(inpart->id()) == 6 || 
         std::abs(inpart->id()) == 15 || 
         std::abs(inpart->id()) == 23 || 
         std::abs(inpart->id()) == 24 || 
         std::abs(inpart->id()) == 25 ) {
      Throw<SetupException>() << "\n"
        "************************************************************************\n"
        "* Your LHE file changes the width of " << inpart->PDGName() << ".\n"
        "* This can cause serious problems in the event generation!\n"
        "************************************************************************\n"
        "\n" << Exception::warning;
    }
    else if (inpart->width() > ZERO && width <= ZERO) {
      Throw<SetupException>() << "\n"
        "************************************************************************\n"
        "* Your LHE file zeroes the non-zero width of " << inpart->PDGName() << ".\n"
        "* If " << inpart->PDGName() << " is a decaying SM particle,\n"
        "*     this can cause serious problems in the event generation!\n"
        "************************************************************************\n"
                  "\n" << Exception::warning;
    }
    // set the width
    inpart->width(width);
    if( width > ZERO ) {
      inpart->cTau(hbarc/width);
      inpart->widthCut(5.*width);
      inpart->stable(false);
    }
    // construct prefix for DecayModes
    string prefix(inpart->name() + "->"), tag(prefix),line("");
    unsigned int nmode(0);
    // read any decay modes
    line = StringUtils::car(block,"\r\n");
    while(line[0] != 'D' && line[0] != 'B' &&
          line[0] != 'd' && line[0] != 'b' && 
          line[0] != '<' && line    != "") {
      // skip comments
      if(line[0] == '#') {
        block = StringUtils::cdr(block,"\r\n");
        line = StringUtils::car(block,"\r\n");
        continue;
      }
      // read decay mode and construct the tag
      istringstream is(line);
      double brat(0.);
      unsigned int nda(0),npr(0);
      is >> brat >> nda;
      while( true ) {
        long t;
        is >> t;
        if( is.fail() ) break; 
        if( t == std::abs(parent) )
          throw SetupException() 
        << "An error occurred while read a decay of the " 
        << inpart->PDGName() << ". One of its products has the same PDG code "
        << "as the parent particle in powhegLesHouchesFileReader::doinit()."
        << " Please check the Les Houches file.\n"
        << Exception::runerror;
        tcPDPtr p = getParticleData(t);
        if( !p )
          throw SetupException()
        << "powhegLesHouchesFileReader::doinit() -"
        << " An unknown PDG code has been encounterd "
        << "while reading a decay mode. ID: " << t
        << Exception::runerror;
        ++npr;
        tag += p->name() + ",";
      }
      if( npr != nda )
        throw SetupException()
          << "powhegLesHouchesFileReader::doinit() - While reading a decay of the " 
          << inpart->PDGName() << " from an SLHA file, an inconsistency "
          << "between the number of decay products and the value in "
          << "the 'NDA' column was found. Please check if the spectrum "
          << "file is correct.\n"
          << Exception::warning;
      // create the DecayMode
      if( npr > 1 ) {
        if( nmode==0 ) {
          generator()->preinitInterface(inpart, "VariableRatio" , "set","false");
          if(inpart->massGenerator()) {
        ok = false;
        Throw<SetupException>()
          << inpart->PDGName() << " already has a MassGenerator set"
          << " this is incompatible with using QNUMBERS "
          << "Use\n"
          << "set " << inpart->fullName() << ":Mass_generator NULL\n"
          << "to fix this." << Exception::warning;
          }
          if(inpart->widthGenerator()) {
        ok = false;
        Throw<SetupException>()
          << inpart->PDGName() << " already has a WidthGenerator set"
          << " this is incompatible with using QNUMBERS "
          << "Use\n"
          << "set " << inpart->fullName() << ":Width_generator NULL\n"
          << "to fix this." << Exception::warning;
          }
          unsigned int ntemp=0;
          for(DecaySet::const_iterator dit = inpart->decayModes().begin();
          dit != inpart->decayModes().end(); ++dit ) {
        if((**dit).on()) ++ntemp;
          }
          if(ntemp!=0) {
        ok = false;
        Throw<SetupException>()
          << inpart->PDGName() << " already has DecayModes"
          << " this is incompatible with using QNUMBERS "
          << "Use\n"
          << "do " << inpart->fullName() << ":SelectDecayModes none\n" 
          << " to fix this." << Exception::warning;
          }
        }
        inpart->stable(false);
        tag.replace(tag.size() - 1, 1, ";");
        DMPtr dm = generator()->findDecayMode(tag);
        if(!m_theDecayer) Throw<SetupException>()
                  << "powhegLesHouchesFileReader::doinit() Decayer must be set using the "
                  << "powhegLesHouchesFileReader:Decayer"
                  << " must be set to allow the creation of new"
                  << " decay modes."
                  << Exception::runerror;
        if(!dm) {
          dm = generator()->preinitCreateDecayMode(tag);
          if(!dm)
        Throw<SetupException>()  
          << "powhegLesHouchesFileReader::doinit() - Needed to create "
          << "new decaymode but one could not be created for the tag " 
          << tag << Exception::warning;
        }
        generator()->preinitInterface(dm, "Decayer", "set",
                      m_theDecayer->fullName());
        ostringstream br;
        br << setprecision(13) << brat;
        generator()->preinitInterface(dm, "BranchingRatio", "set", br.str());
        generator()->preinitInterface(dm, "Active", "set", "Yes");
        if(dm->CC()) {
          generator()->preinitInterface(dm->CC(), "BranchingRatio", "set", br.str());
          generator()->preinitInterface(dm->CC(), "Active", "set", "Yes");
        }
        ++nmode;
      }
      tag=prefix;
      // read the next line
      block = StringUtils::cdr(block,"\r\n");
      line = StringUtils::car(block,"\r\n");
    };
    if(nmode>0) {
      inpart->update();
      if(inpart->CC())
        inpart->CC()->update();
    }
      }
    }
    // start of SLHA block ?
    else if(line.find("<slha") != string::npos) {
      readingSLHA = true;
    }
  }
  while(line!="");
  if(!ok)
    throw SetupException() << "Problem reading QNUMBERS blocks in powhegLesHouchesFileReader::doinit()"
               << Exception::runerror;
}

doReadEvent()

bool powhegLesHouchesFileReader::doReadEvent ( )

virtual

Read the next event from the file or stream into the corresponding protected variables.

Return false if there is no more events or if this was not a LHF event file.

Definition at line 624 of file PowhegLesHouchesFileReader.cxx.

                                             {
  if ( !m_cfile ) return false;
  if ( m_LHFVersion.empty() ) return false;
  if ( heprup.NPRUP < 0 ) return false;
  m_eventComments = "";
  m_outsideBlock = "";
  hepeup.NUP = 0;
  hepeup.XPDWUP.first = hepeup.XPDWUP.second = 0.0;
  optionalWeights.clear();
  m_optionalWeightsTemp.clear();
  // Keep reading lines until we hit the next event or the end of
  // the event block. Save any inbetween lines. Exit if we didn't
  // find an event.
  while ( m_cfile.readline() && !m_cfile.find("<event") )
    m_outsideBlock += m_cfile.getline() + "\n";
 
  // We found an event. First scan for attributes.
  m_eventAttributes = StringUtils::xmlAttributes("event", m_cfile.getline());
 
  /* information necessary for FxFx merging:
   * the npLO and npNLO tags
   */
  istringstream ievat(m_cfile.getline());
  int we(0), npLO(-99), npNLO(-99);
  do {
    string sub; ievat >> sub;
    if(we==2) { npLO = atoi(sub.c_str()); }
    if(we==5) { npNLO = atoi(sub.c_str()); }
    ++we;
  } while (ievat);
  optionalnpLO = npLO;
  optionalnpNLO = npNLO;
  std::stringstream npstringstream;
  npstringstream << "np " << npLO << " " << npNLO;
  std::string npstrings = npstringstream.str();
  /* the FxFx merging information 
   * becomes part of the optionalWeights, labelled -999 
   * for future reference
   */
 
  if(m_theIncludeFxFxTags) optionalWeights[npstrings.c_str()] = -999;
 
  if ( !m_cfile.readline()  ) return false;
 
  // The first line determines how many subsequent particle lines we
  // have.
  if ( !( m_cfile >> hepeup.NUP >> hepeup.IDPRUP >> hepeup.XWGTUP
            >> hepeup.SCALUP >> hepeup.AQEDUP >> hepeup.AQCDUP ) )
    return false;
  hepeup.resize();
  // Read all particle lines.
  for ( int i = 0; i < hepeup.NUP; ++i ) {
    if ( !m_cfile.readline() ) return false;
    if ( !( m_cfile >> 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] ) )
      return false;
    //print momenta for debugging
    if(std::isnan(hepeup.PUP[i][0])||std::isnan(hepeup.PUP[i][1])||
       std::isnan(hepeup.PUP[i][2])||std::isnan(hepeup.PUP[i][3])||
       std::isnan(hepeup.PUP[i][4])) 
      throw Exception() 
    << "nan's as momenta in Les Houches file "
    << Exception::eventerror;
    if(hepeup.MOTHUP[i].first -1==i || 
       hepeup.MOTHUP[i].second-1==i) {
      throw Exception()
    << "Particle has itself as a mother in Les Houches "
    << "file, this is not allowed\n"
    << Exception::eventerror;
    } 
  }
 
  // Now read any additional comments and named weights.
  // read until the end of rwgt is found
  bool readingWeights = false, readingaMCFast = false, readingMG5ClusInfo = false;
  while ( m_cfile.readline() && !m_cfile.find("</event>")) {
    
    if(m_cfile.find("</applgrid")) { readingaMCFast=false; } //aMCFast weights end
    if(m_cfile.find("</rwgt")) { readingWeights=false; } //optional weights end
    if(m_cfile.find("</clustering")) { readingMG5ClusInfo=false; } // mg5 mclustering scale info end
    
    /* reading of optional weights
     */
 
    /* BEGIN MOD*/
    if(readingWeights) { 
      if(!m_cfile.find("<wgt")) { continue; }
      istringstream iss(m_cfile.getline());
      double weightValue(0);
      std::string hs=m_cfile.getline();
      hs.erase(std::remove(hs.begin(), hs.end(), ' '), hs.end());
      std::string IdLabel = hs.substr(0, hs.find(">", 0));
      IdLabel = std::regex_replace(IdLabel, std::regex(R"([\D])"), "");
      std::string weightName = m_optionalWeightsLabel[IdLabel];
      std::string value=hs;
      erase_substr(value, "<wgtid='"+IdLabel+"'>");
      erase_substr(value, "</wgt>");
      erase_substr(value, "\n");
      weightValue=std::stod(value);
      m_optionalWeightsTemp[weightName] = weightValue; 
      optionalWeights[weightName] = weightValue; 
    }
    /* END MOD */
    
    /* reading of aMCFast weights
     */
    if(readingaMCFast) {
      std::stringstream amcfstringstream;
      amcfstringstream << "aMCFast " << m_cfile.getline();
      std::string amcfstrings = amcfstringstream.str();
      string str_newline = "\n";
      erase_substr(amcfstrings,str_newline);
      optionalWeights[amcfstrings.c_str()] = -111; //for the aMCFast we give them a weight -111 for future identification
    }
 
    /* read additional MG5 Clustering information 
     * used in LO merging
     */
    if(readingMG5ClusInfo) {
      string hs = m_cfile.getline();
      string startDEL = "<clus scale="; //starting delimiter
      string stopDEL = "</clus>"; //end delimiter
      unsigned firstLim = hs.find(startDEL); //find start of delimiter
      string mg5clusinfo = hs.substr(firstLim); //define the information for the scale
      erase_substr(mg5clusinfo,stopDEL);
      erase_substr(mg5clusinfo,startDEL);
      string str_arrow = ">";
      erase_substr(mg5clusinfo,str_arrow);
      string str_quotation = "\"";
      erase_substr(mg5clusinfo,str_quotation);
      optionalWeights[mg5clusinfo.c_str()] = -222; //for the mg5 scale info weights we give them a weight -222 for future identification
    }
    
    //store MG5 clustering information 
    if(m_cfile.find("<scales")) {
      string hs = m_cfile.getline();
      string startDEL = "<scales"; //starting delimiter
      string stopDEL = "</scales>"; //end delimiter
      unsigned firstLim = hs.find(startDEL); //find start of delimiter
      string mg5scinfo = hs.substr(firstLim); //define the information for the scale
      erase_substr(mg5scinfo,stopDEL);
      erase_substr(mg5scinfo,startDEL);
      string str_arrow = ">";
      erase_substr(mg5scinfo,str_arrow);
      optionalWeights[mg5scinfo.c_str()] = -333; //for the mg5 scale info weights we give them a weight -333 for future identification
    }
 
    //determine start of aMCFast weights
    if(m_cfile.find("<applgrid")) { readingaMCFast = true;}
    //determine start of optional weights
    if(m_cfile.find("<rwgt")) { readingWeights = true; }
    //determine start of MG5 clustering scale information
    if(m_cfile.find("<clustering")) { readingMG5ClusInfo = true;}
  }
  string str_quote = "'";
  string str_doublequote = "\"";
  // loop over the optional weights and add the extra information as found in the init
  for (map<string,double>::const_iterator it=m_optionalWeightsTemp.begin(); it!=m_optionalWeightsTemp.end(); ++it){
    //to avoid issues with inconsistencies of defining the weight ids, remove "" and ''
    string id_1 = it->first;
    erase_substr(id_1, str_quote);
    erase_substr(id_1, str_doublequote);
    for (map<string,string>::const_iterator it2=m_scalemap.begin(); it2!=m_scalemap.end(); ++it2){
      //find the scale id in the scale information and add this information
      string id_2 = it2->first;
      erase_substr(id_2, str_quote);
      erase_substr(id_2, str_doublequote);
      if(id_1==id_2) {
        string info = it2->second;
    string str_newline = "\n";
    erase_substr(info, str_newline);
    //set the optional weights
    optionalWeights[info] = it->second;
      }
    }
  }
  /* additionally, we set the "central" scale
   * this is actually the default event weight 
   */
  string central = "central";
  if (m_theIncludeCentral) optionalWeights[central] = hepeup.XWGTUP;
  
  if ( !m_cfile ) return false;
  return true;
  
}

erase_substr()

void powhegLesHouchesFileReader::erase_substr	(	std::string &	subject,
		const std::string &	search
	)

Erases all occurences of a substring from a string.

Definition at line 912 of file PowhegLesHouchesFileReader.cxx.

                                                                                         {
    size_t pos = 0;
    while((pos = subject.find(search, pos)) != std::string::npos) {
      subject.erase( pos, search.length() );
    }
}

filename()

string ThePEG::powhegLesHouchesFileReader::filename ( ) const

inline

Return the name of the file from where to read events.

Definition at line 102 of file PowhegLesHouchesFileReader.h.

{ return m_theFileName; }

fullclone()

IBPtr powhegLesHouchesFileReader::fullclone ( ) const

protectedvirtual

Make a clone of this object, possibly modifying the cloned object to make it sane.

Returns

a pointer to the new object.

Definition at line 60 of file PowhegLesHouchesFileReader.cxx.

                                                  {
  return new_ptr(*this);
}

initialize()

void powhegLesHouchesFileReader::initialize ( LesHouchesEventHandler &  eh )

virtual

Initialize.

This function is called by the LesHouchesEventHandler to which this object is assigned.

Definition at line 469 of file PowhegLesHouchesFileReader.cxx.

                                                                       {
  LesHouchesReader::initialize(eh);
  if ( m_LHFVersion.empty() )
    Throw<powhegLesHouchesFileError>()
      << "The file associated with '" << name() << "' does not contain a "
      << "proper formatted Les Houches event file. The events may not be "
      << "properly sampled." << Exception::warning;
}

open()

void powhegLesHouchesFileReader::open ( )

virtual

Open a file with events.

Derived classes should overwrite it and first calling it before reading in the run information into the corresponding protected variables.

Definition at line 480 of file PowhegLesHouchesFileReader.cxx.

                                      {
  if ( filename().empty() )
    throw powhegLesHouchesFileError()
      << "No Les Houches file name. "
      << "Use 'set " << name() << ":FileName'."
      << Exception::runerror;
  m_cfile.open(filename());
  if ( !m_cfile )
    throw powhegLesHouchesFileError()
      << "The powhegLesHouchesFileReader '" << name() << "' could not open the "
      << "event file called '" << m_theFileName << "'."
      << Exception::runerror;
 
  
  m_cfile.readline();
  if ( !m_cfile.find("<LesHouchesEvents") ) return;
  map<string,string> attributes =
    StringUtils::xmlAttributes("LesHouchesEvents", m_cfile.getline());
  m_LHFVersion = attributes["version"];
  if ( m_LHFVersion.empty() ) return;
 
  bool readingHeader = false;
  bool readingInit = false;
  m_headerBlock = "";
  string hs;
 
  // Loop over all lines until we hit the </init> tag.
  bool readingInitWeights = false, readingInitWeights_sc = false;
  string weightinfo;
  while ( m_cfile.readline() ) {
    // found the init block for multiple weights
    if(m_cfile.find("<initrwgt")) readingInitWeights = true;
    
    // found end of init block for multiple weights: end the while loop
    if(m_cfile.find("</initrwgt")) {
      readingInitWeights = false;
      readingInitWeights_sc = false;
      continue;
    }
 
    // found the end of init block
    if(m_cfile.find("</init")) {
      readingInit = false;
      break;
    } 
 
    /* read the weight information block 
     * optionalWeightNames will contain information about the weights
     * this will be appended to the weight information
     */ 
    if(readingInitWeights) {
      string scalename = "";
      if(m_cfile.find("<weightgroup")) {
    /* we found a weight group
     * start reading the information 
     * within it
     */
    readingInitWeights_sc = true;
    weightinfo = m_cfile.getline();
    /* to make it shorter, erase some stuff
     */
    string str_weightgroup = "<weightgroup";
    string str_arrow = ">";
    string str_newline = "\n";
    erase_substr(weightinfo, str_weightgroup);
    erase_substr(weightinfo, str_arrow);
    erase_substr(weightinfo, str_newline);
    continue;
      }
      /* if we are reading a new weightgroup, go on 
       * until we find the end of it
       */
      /* BEGIN MOD*/
      if(readingInitWeights_sc && !m_cfile.find("</weightgroup")) {
    hs = m_cfile.getline();
    hs.erase(std::remove(hs.begin(), hs.end(), ' '), hs.end());
    std::string IdLabel = hs.substr(0, hs.find(">", 0));
    IdLabel = std::regex_replace(IdLabel, std::regex(R"([\D])"), "");
    std::string name = hs;
    erase_substr(name, "<weightid='"+IdLabel+"'>");
    name.erase(name.find('<'));
    m_optionalWeightsLabel[IdLabel]=name;   
    optionalWeightsNames.push_back(name);
      }
      /*END MOD*/
    }
   
    if ( m_cfile.find("<header") ) {
      // We have hit the header block, so we should dump this and all
      // following lines to m_headerBlock until we hit the end of it.
      readingHeader = true;
      m_headerBlock = m_cfile.getline() + "\n";
    }
    if ( (m_cfile.find("<init ") && !m_cfile.find("<initrwgt")) || m_cfile.find("<init>") ) {
      // We have hit the init block, so we should expect to find the
      // standard information in the following. But first check for
      // attributes.
      m_initAttributes = StringUtils::xmlAttributes("init", m_cfile.getline());
      readingInit = true;
      m_cfile.readline();
      if ( !( m_cfile >> heprup.IDBMUP.first >> heprup.IDBMUP.second
            >> heprup.EBMUP.first >> heprup.EBMUP.second
                >> heprup.PDFGUP.first >> heprup.PDFGUP.second
                >> heprup.PDFSUP.first >> heprup.PDFSUP.second
            >> heprup.IDWTUP >> heprup.NPRUP ) ) {
    heprup.NPRUP = -42;
    m_LHFVersion = "";
    return;
      }
      heprup.resize();
      for ( int i = 0; i < heprup.NPRUP; ++i ) {
    m_cfile.readline();
    if ( !( m_cfile >> heprup.XSECUP[i] >> heprup.XERRUP[i]
                  >> heprup.XMAXUP[i] >> heprup.LPRUP[i] ) ) {
      heprup.NPRUP = -42;
      m_LHFVersion = "";
      return;
    }
      }
    }
    if ( m_cfile.find("</header") ) {
      readingHeader = false;
      m_headerBlock += m_cfile.getline() + "\n";
    }
    if ( readingHeader ) {
      /* We are in the process of reading the header block. Dump the
     line to m_headerBlock.*/
      m_headerBlock += m_cfile.getline() + "\n";
    }
    if ( readingInit ) {
      // Here we found a comment line. Dump it to m_initComments.
      m_initComments += m_cfile.getline() + "\n";
    }
  }
  string central = "central";
  optionalWeightsNames.push_back(central);
 
  if ( !m_cfile ) {
    heprup.NPRUP = -42;
    m_LHFVersion = "";
    return;
  }
}

operator=()

powhegLesHouchesFileReader& ThePEG::powhegLesHouchesFileReader::operator= ( const powhegLesHouchesFileReader &  )

private

Private and non-existent assignment operator.

optWeightsNamesFunc()

vector< string > powhegLesHouchesFileReader::optWeightsNamesFunc ( )

virtual

Return the optional weights information string ("Names")

Definition at line 478 of file PowhegLesHouchesFileReader.cxx.

{ return optionalWeightsNames; }

persistentInput()

void powhegLesHouchesFileReader::persistentInput	(	PersistentIStream &	is,
		int	version
	)

Function used to read in object persistently.

Parameters

is	the persistent input stream read from.
version	the version number of the object when written.

Definition at line 825 of file PowhegLesHouchesFileReader.cxx.

                                                                            {
  is >> m_neve >> m_LHFVersion >> m_outsideBlock >> m_headerBlock >> m_initComments
     >> m_initAttributes >> m_eventComments >> m_eventAttributes >> m_theFileName
     >> m_theQNumbers >> m_theIncludeFxFxTags >> m_theIncludeCentral >> m_theDecayer;
  m_ieve = 0;
}

persistentOutput()

void powhegLesHouchesFileReader::persistentOutput

(

PersistentOStream &

os

)

const

Function used to write out object persistently.

Parameters

os	the persistent output stream written to.

Definition at line 819 of file PowhegLesHouchesFileReader.cxx.

                                                                              {
  os << m_neve << m_LHFVersion << m_outsideBlock << m_headerBlock << m_initComments
     << m_initAttributes << m_eventComments << m_eventAttributes << m_theFileName
     << m_theQNumbers << m_theIncludeFxFxTags << m_theIncludeCentral << m_theDecayer ;
}

preInitialize()

bool powhegLesHouchesFileReader::preInitialize ( ) const

protectedvirtual

Return true if this object needs to be initialized before all other objects because it needs to extract PDFs from the event file.

Definition at line 64 of file PowhegLesHouchesFileReader.cxx.

                                                     {
  return true;
}

Member Data Documentation

m_cfile

CFileLineReader ThePEG::powhegLesHouchesFileReader::m_cfile

protected

The wrapper around the C FILE stream from which to read.

Definition at line 180 of file PowhegLesHouchesFileReader.h.

m_eventAttributes

map<string,string> ThePEG::powhegLesHouchesFileReader::m_eventAttributes

protected

If LHF.

Map of attributes (name-value pairs) found in the last event tag.

Definition at line 232 of file PowhegLesHouchesFileReader.h.

m_eventComments

string ThePEG::powhegLesHouchesFileReader::m_eventComments

protected

If LHF.

Additional comments found with the last read event.

Definition at line 226 of file PowhegLesHouchesFileReader.h.

m_headerBlock

string ThePEG::powhegLesHouchesFileReader::m_headerBlock

protected

If LHF.

All lines from the header block.

Definition at line 210 of file PowhegLesHouchesFileReader.h.

m_ieve

long ThePEG::powhegLesHouchesFileReader::m_ieve

protected

The current event number.

Definition at line 192 of file PowhegLesHouchesFileReader.h.

m_initAttributes

map<string,string> ThePEG::powhegLesHouchesFileReader::m_initAttributes

protected

If LHF.

Map of attributes (name-value pairs) found in the init tag.

Definition at line 221 of file PowhegLesHouchesFileReader.h.

m_initComments

string ThePEG::powhegLesHouchesFileReader::m_initComments

protected

If LHF.

Additional comments found in the init block.

Definition at line 215 of file PowhegLesHouchesFileReader.h.

m_initpowhegLesHouchesFileReader

const ClassDescription< powhegLesHouchesFileReader > powhegLesHouchesFileReader::m_initpowhegLesHouchesFileReader

staticprivate

Describe an abstract base class with persistent data.

Definition at line 280 of file PowhegLesHouchesFileReader.h.

m_LHFVersion

string ThePEG::powhegLesHouchesFileReader::m_LHFVersion

protected

If the file is a standard Les Houches formatted file (LHF) this is its version number.

If empty, this is not a Les Houches formatted file

Definition at line 199 of file PowhegLesHouchesFileReader.h.

m_neve

long ThePEG::powhegLesHouchesFileReader::m_neve

protected

The number of events in this file.

Definition at line 187 of file PowhegLesHouchesFileReader.h.

m_optionalWeightsLabel

map<string,string> ThePEG::powhegLesHouchesFileReader::m_optionalWeightsLabel

private

Definition at line 272 of file PowhegLesHouchesFileReader.h.

m_optionalWeightsTemp

map<string,double> ThePEG::powhegLesHouchesFileReader::m_optionalWeightsTemp

private

Temporary holder for optional weights.

Definition at line 270 of file PowhegLesHouchesFileReader.h.

m_outsideBlock

string ThePEG::powhegLesHouchesFileReader::m_outsideBlock

protected

If LHF.

All lines (since the last open() or readEvent()) outside the header, init and event tags.

Definition at line 205 of file PowhegLesHouchesFileReader.h.

m_scalemap

map<string,string> ThePEG::powhegLesHouchesFileReader::m_scalemap

private

Further information on the weights.

Definition at line 264 of file PowhegLesHouchesFileReader.h.

m_theDecayer

DecayerPtr ThePEG::powhegLesHouchesFileReader::m_theDecayer

private

Decayer for any decay modes read from the file.

Definition at line 259 of file PowhegLesHouchesFileReader.h.

m_theFileName

string ThePEG::powhegLesHouchesFileReader::m_theFileName

private

The name of the file from where to read events.

Definition at line 239 of file PowhegLesHouchesFileReader.h.

m_theIncludeCentral

bool ThePEG::powhegLesHouchesFileReader::m_theIncludeCentral

private

Include central weight (for backup use)

Definition at line 254 of file PowhegLesHouchesFileReader.h.

m_theIncludeFxFxTags

bool ThePEG::powhegLesHouchesFileReader::m_theIncludeFxFxTags

private

Include/Read FxFx tags.

Definition at line 249 of file PowhegLesHouchesFileReader.h.

m_theQNumbers

bool ThePEG::powhegLesHouchesFileReader::m_theQNumbers

private

Whether or not to search for QNUMBERS stuff.

Definition at line 244 of file PowhegLesHouchesFileReader.h.

---

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

• PowhegLesHouchesFileReader.h
• PowhegLesHouchesFileReader.cxx