LHEEventProduct.cc

Go to the documentation of this file.

#include <iostream>
#include <iomanip>
#include <sstream>
#include <string>

#include "SimDataFormats/GeneratorProducts/interface/LesHouches.h"
#include "SimDataFormats/GeneratorProducts/interface/LHEEventProduct.h"

void LHEEventProduct::const_iterator::next()
{
    const lhef::HEPEUP &hepeup = event->hepeup();
    int line = this->line++;

    if (!line) {
        std::ostringstream ss;
        ss << std::setprecision(7)
           << std::scientific
           << std::uppercase
           << "    " << hepeup.NUP
           << "  " << hepeup.IDPRUP
           << "  " << event->originalXWGTUP()
           << "  " << hepeup.SCALUP
           << "  " << hepeup.AQEDUP
           << "  " << hepeup.AQCDUP << std::endl;
        tmp = ss.str();
        return;
    }
    line--;

    if (line < hepeup.NUP) {
        std::ostringstream ss;
        ss << std::setprecision(10)
           << std::scientific
           << std::uppercase
           << "\t" << hepeup.IDUP[line]
           << "\t" << hepeup.ISTUP[line]
           << "\t" << hepeup.MOTHUP[line].first
           << "\t" << hepeup.MOTHUP[line].second
           << "\t" << hepeup.ICOLUP[line].first
           << "\t" << hepeup.ICOLUP[line].second
           << "\t" << hepeup.PUP[line][0]
           << "\t" << hepeup.PUP[line][1]
           << "\t" << hepeup.PUP[line][2]
           << "\t" << hepeup.PUP[line][3]
           << "\t" << hepeup.PUP[line][4]
           << std::setprecision(3)
           << "\t" << hepeup.VTIMUP[line]
           << std::setprecision(1)
           << std::fixed
           << "\t" << hepeup.SPINUP[line] << std::endl;
        tmp = ss.str();
        return;
    }
    line -= hepeup.NUP;

    if (event->pdf()) {
        if (!line) {
            const PDF &pdf = *event->pdf();
            std::ostringstream ss;
            ss << std::setprecision(7)
               << std::scientific
               << std::uppercase
               << "#pdf  " << pdf.id.first
               << "  " << pdf.id.second
               << "  " << pdf.x.first
               << "  " << pdf.x.second
               << "  " << pdf.scalePDF
               << "  " << pdf.xPDF.first
               << "  " << pdf.xPDF.second << std::endl;
            tmp = ss.str();
            return;
        }
        line--;
    }

    if (line < (int)event->comments_size()) {
        tmp = *(event->comments_begin() + line);
        return;
    }
    line -= event->comments_size();

    if (!line) {
        tmp = "</event>\n";
        return;
    }

    tmp.clear();
    this->line = npos;
}

LHEEventProduct::const_iterator LHEEventProduct::begin() const
{
    const_iterator result;
    result.event = this;
    result.line = 0;
    result.tmp = "<event>\n";
    return result;
}