PATVertexSlimmer.cc

Go to the documentation of this file.

#include <string>
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/ValueMap.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/Math/interface/libminifloat.h"

namespace pat {
  class PATVertexSlimmer : public edm::global::EDProducer<> {
  public:
    explicit PATVertexSlimmer(const edm::ParameterSet&);
    ~PATVertexSlimmer() override;

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

  private:
    const edm::EDGetTokenT<std::vector<reco::Vertex>> src_;
    const edm::EDGetTokenT<edm::ValueMap<float>> score_;
    const bool rekeyScores_;
  };
}  // namespace pat

pat::PATVertexSlimmer::PATVertexSlimmer(const edm::ParameterSet& iConfig)
    : src_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("src"))),
      score_(mayConsume<edm::ValueMap<float>>(
          iConfig.existsAs<edm::InputTag>("score") ? iConfig.getParameter<edm::InputTag>("score") : edm::InputTag())),
      rekeyScores_(iConfig.existsAs<edm::InputTag>("score")) {
  produces<std::vector<reco::Vertex>>();
  if (rekeyScores_)
    produces<edm::ValueMap<float>>();
}

pat::PATVertexSlimmer::~PATVertexSlimmer() {}

void pat::PATVertexSlimmer::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup& iSetup) const {
  edm::Handle<std::vector<reco::Vertex>> vertices;
  iEvent.getByToken(src_, vertices);
  auto outPtr = std::make_unique<std::vector<reco::Vertex>>();

  outPtr->reserve(vertices->size());
  for (unsigned int i = 0, n = vertices->size(); i < n; ++i) {
    const reco::Vertex& v = (*vertices)[i];
    auto co = v.covariance4D();
    if (i > 0) {
      for (size_t j = 0; j < 4; j++) {
        for (size_t k = j; k < 4; k++) {
          co(j, k) = MiniFloatConverter::reduceMantissaToNbits<10>(co(j, k));
        }
      }
    }
    outPtr->push_back(reco::Vertex(v.position(), co, v.t(), v.chi2(), v.ndof(), 0));
  }

  auto oh = iEvent.put(std::move(outPtr));
  if (rekeyScores_) {
    edm::Handle<edm::ValueMap<float>> scores;
    iEvent.getByToken(score_, scores);
    auto vertexScoreOutput = std::make_unique<edm::ValueMap<float>>();
    edm::ValueMap<float>::const_iterator idIt = scores->begin();
    for (; idIt != scores->end(); idIt++) {
      if (idIt.id() == vertices.id())
        break;
    }
    // std::find_if(scores->begin(), scores->end(), [vertices] (const edm::ValueMap<float>::const_iterator& s) { return s.id() == vertices.id(); } );
    edm::ValueMap<float>::Filler vertexScoreFiller(*vertexScoreOutput);
    vertexScoreFiller.insert(oh, idIt.begin(), idIt.end());
    vertexScoreFiller.fill();
    iEvent.put(std::move(vertexScoreOutput));
  }
}

using pat::PATVertexSlimmer;
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(PATVertexSlimmer);