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 &&other)
    : weights_(std::move(other.weights_)),
      signalProcessID_(other.signalProcessID_),
      qScale_(other.qScale_),
      alphaQCD_(other.alphaQCD_),
      alphaQED_(other.alphaQED_),
      pdf_(other.pdf_.release()),
      binningValues_(std::move(other.binningValues_)),
      DJRValues_(std::move(other.DJRValues_)),
      nMEPartons_(other.nMEPartons_),
      nMEPartonsFiltered_(other.nMEPartons_) {}

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;
}

GenEventInfoProduct &GenEventInfoProduct::operator=(GenEventInfoProduct &&other) {
  weights_ = std::move(other.weights_);
  signalProcessID_ = other.signalProcessID_;
  qScale_ = other.qScale_;
  alphaQCD_ = other.alphaQCD_;
  alphaQED_ = other.alphaQED_;
  binningValues_ = std::move(other.binningValues_);
  DJRValues_ = std::move(other.DJRValues_);
  nMEPartons_ = other.nMEPartons_;
  nMEPartonsFiltered_ = other.nMEPartonsFiltered_;
  pdf_.reset(other.pdf_.release());

  return *this;
}

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