WeightManager.cc

Go to the documentation of this file.

#include "Validation/EventGenerator/interface/WeightManager.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/Common/interface/Handle.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"

using namespace edm;

WeightManager::WeightManager(const ParameterSet& iConfig, edm::ConsumesCollector iC)
    : _useHepMC(iConfig.getParameter<bool>("UseWeightFromHepMC")) {
  if (_useHepMC) {
    _hepmcCollection = iConfig.getParameter<InputTag>("hepmcCollection");
    hepmcCollectionToken_ = iC.consumes<HepMCProduct>(_hepmcCollection);
  } else {
    _genEventInfos = iConfig.getParameter<std::vector<InputTag>>("genEventInfos");
    for (unsigned int i = 0; i < _genEventInfos.size(); i++)
      genEventInfosTokens_.push_back(iC.consumes<GenEventInfoProduct>(_genEventInfos[i]));
  }
}

double WeightManager::weight(const Event& iEvent) {
  if (_useHepMC) {
    edm::Handle<HepMCProduct> evt;
    iEvent.getByToken(hepmcCollectionToken_, evt);
    const HepMC::GenEvent* myGenEvent = evt->GetEvent();

    double weight = 1.;
    if (!myGenEvent->weights().empty())
      weight = myGenEvent->weights()[0];
    return weight;
  } else {
    double weight = 1.;
    for (unsigned int i = 0; i < genEventInfosTokens_.size(); ++i) {
      edm::Handle<GenEventInfoProduct> info;
      iEvent.getByToken(genEventInfosTokens_[i], info);
      weight *= info->weight();
    }
    return weight;
  }
}

std::vector<std::vector<double>> WeightManager::weightsCollection(const Event& iEvent) {
  std::vector<std::vector<double>> weightsCollection;
  for (unsigned int i = 0; i < genEventInfosTokens_.size(); ++i) {
    edm::Handle<GenEventInfoProduct> info;
    iEvent.getByToken(genEventInfosTokens_[i], info);
    weightsCollection.push_back(info->weights());
  }

  return weightsCollection;
}