GenMETExtractor.cc

Go to the documentation of this file.

#include "DataFormats/METReco/interface/GenMET.h"
#include "DataFormats/PatCandidates/interface/MET.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include <memory>

namespace pat {

  class GenMETExtractor : public edm::global::EDProducer<> {
  public:
    explicit GenMETExtractor(const edm::ParameterSet &iConfig);
    ~GenMETExtractor() override;

    void produce(edm::StreamID streamID, edm::Event &iEvent, const edm::EventSetup &iSetup) const override;

  private:
    edm::EDGetTokenT<std::vector<pat::MET> > metSrcToken_;
  };

}  // namespace pat

using namespace pat;

GenMETExtractor::GenMETExtractor(const edm::ParameterSet &iConfig) {
  edm::InputTag metIT = iConfig.getParameter<edm::InputTag>("metSource");
  metSrcToken_ = consumes<pat::METCollection>(metIT);

  // produces vector of genMet
  produces<std::vector<reco::GenMET> >();
}

GenMETExtractor::~GenMETExtractor() {}

void GenMETExtractor::produce(edm::StreamID streamID, edm::Event &iEvent, const edm::EventSetup &iSetup) const {
  edm::Handle<std::vector<pat::MET> > src;
  iEvent.getByToken(metSrcToken_, src);
  if (src->empty())
    edm::LogError("GenMETExtractor::produce") << "input genMET collection is empty";

  const reco::GenMET *genMet = src->front().genMET();

  std::vector<reco::GenMET> *genMetCol = new std::vector<reco::GenMET>();
  genMetCol->push_back((*genMet));

  std::unique_ptr<std::vector<reco::GenMET> > genMETs(genMetCol);
  iEvent.put(std::move(genMETs));
}

#include "FWCore/Framework/interface/MakerMacros.h"

DEFINE_FWK_MODULE(GenMETExtractor);