LesHouchesReader.hh

Go to the documentation of this file.

 
#pragma once
 
#include <string>
 
#include "LHEF/LHEF.h"
 
#include "HEJ/EventReader.hh"
#include "HEJ/stream.hh"
 
namespace HEJ {
 
  class LesHouchesReader : public EventReader{
  public:
    explicit LesHouchesReader(std::string const & filename):
      stream_{filename},
      reader_{stream_}
    {
      // always use the newest LHE version
      reader_.heprup.version = LHEF::HEPRUP().version;
    }
 
    bool read_event() override {
      return reader_.readEvent();
    }
 
    std::string const & header() const override {
      return reader_.headerBlock;
    }
 
    LHEF::HEPRUP const & heprup() const override {
      return reader_.heprup;
    }
 
    LHEF::HEPEUP const & hepeup() const override {
      return reader_.hepeup;
    }
 
  private:
    istream stream_;
    LHEF::Reader reader_;
 
  protected:
    LHEF::Reader & reader(){
      return reader_;
    }
  };
 
  class SherpaLHEReader : public LesHouchesReader{
  public:
    explicit SherpaLHEReader(std::string const & filename);
 
    bool read_event() override;
 
    optional<size_t> number_events() const override {
      return num_events_;
    }
 
  private:
    double num_trials_;
    size_t num_events_;
  };
 
} // namespace HEJ