GenEventInfoProduct.cc

Go to the documentation of this file.

#include <functional>
#include <numeric>
using std::ptrdiff_t;

#include <HepMC/GenEvent.h>
#include <HepMC/WeightContainer.h>
#include <HepMC/PdfInfo.h>

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"

using namespace edm;
using namespace std;

GenEventInfoProduct::GenEventInfoProduct() :
    signalProcessID_(0), qScale_(-1.), alphaQCD_(-1.), alphaQED_(-1.),
    nMEPartons_(-1), nMEPartonsFiltered_(-1)
{
}

GenEventInfoProduct::GenEventInfoProduct(const HepMC::GenEvent *evt) :
    weights_(evt->weights().begin(), evt->weights().end()),
    signalProcessID_(evt->signal_process_id()),
    qScale_(evt->event_scale()),
    alphaQCD_(evt->alphaQCD()),
    alphaQED_(evt->alphaQED()),
    nMEPartons_(-1), nMEPartonsFiltered_(-1)
{
    const HepMC::PdfInfo *hepPDF = evt->pdf_info();    
    if (hepPDF) {
        PDF pdf;

        pdf.id = std::make_pair(hepPDF->id1(), hepPDF->id2());
        pdf.x = std::make_pair(hepPDF->x1(), hepPDF->x2());
        pdf.xPDF = std::make_pair(hepPDF->pdf1(), hepPDF->pdf2());   
        pdf.scalePDF = hepPDF->scalePDF();

        setPDF(&pdf);
    }
}

GenEventInfoProduct::GenEventInfoProduct(GenEventInfoProduct const &other) :
    weights_(other.weights_),
    signalProcessID_(other.signalProcessID_),
    qScale_(other.qScale_),
    alphaQCD_(other.alphaQCD_),
    alphaQED_(other.alphaQED_),
    binningValues_(other.binningValues_),
    DJRValues_(other.DJRValues_),
    nMEPartons_(other.nMEPartons_), nMEPartonsFiltered_(other.nMEPartons_)
{
    setPDF(other.pdf());
}

GenEventInfoProduct::~GenEventInfoProduct()
{
}

GenEventInfoProduct &GenEventInfoProduct::operator = (GenEventInfoProduct const &other)
{
    weights_ = other.weights_;
    signalProcessID_ = other.signalProcessID_;
    qScale_ = other.qScale_;
    alphaQCD_ = other.alphaQCD_;
    alphaQED_ = other.alphaQED_;
    binningValues_ = other.binningValues_;
    DJRValues_ = other.DJRValues_;
    nMEPartons_=other.nMEPartons_;
        nMEPartonsFiltered_=other.nMEPartonsFiltered_;

    setPDF(other.pdf());

    return *this;
}

double GenEventInfoProduct::weight() const
{
    return std::accumulate(weights_.begin(), weights_.end(),
                           1., std::multiplies<double>());
}