CandPtrProjector.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/Common/interface/PtrVector.h"
#include "DataFormats/Common/interface/RefToBaseVector.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/Math/interface/deltaR.h"

class CandPtrProjector : public edm::global::EDProducer<> {
public:

  explicit CandPtrProjector(edm::ParameterSet const& iConfig);
  void produce(edm::StreamID, edm::Event&, edm::EventSetup const&) const override;
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

private:
  edm::EDGetTokenT<edm::View<reco::Candidate>> candSrcToken_;
  edm::EDGetTokenT<edm::View<reco::Candidate>> vetoSrcToken_;
  bool useDeltaRforFootprint_;
};

CandPtrProjector::CandPtrProjector(edm::ParameterSet const& iConfig):
  candSrcToken_{consumes<edm::View<reco::Candidate>>(iConfig.getParameter<edm::InputTag>("src"))},
  vetoSrcToken_{consumes<edm::View<reco::Candidate>>(iConfig.getParameter<edm::InputTag>("veto"))},
  useDeltaRforFootprint_(iConfig.getParameter<bool>("useDeltaRforFootprint"))
{
  produces<edm::PtrVector<reco::Candidate>>();
}

void
CandPtrProjector::produce(edm::StreamID, edm::Event& iEvent, edm::EventSetup const&) const
{
  using namespace edm;
  Handle<View<reco::Candidate>> vetoes;
  iEvent.getByToken(vetoSrcToken_, vetoes);

  auto result = std::make_unique<PtrVector<reco::Candidate>>();
  std::set<reco::CandidatePtr> vetoedPtrs;
  for (auto const& veto : *vetoes) {
    auto const n = veto.numberOfSourceCandidatePtrs();
    for (size_t j {}; j<n; ++j) {
      vetoedPtrs.insert(veto.sourceCandidatePtr(j));
    }
  }

  Handle<View<reco::Candidate>> cands;
  iEvent.getByToken(candSrcToken_, cands);
  for (size_t i {}; i<cands->size(); ++i) {
    auto const c = cands->ptrAt(i);
    if (vetoedPtrs.find(c)==vetoedPtrs.cend()) {
      bool addcand = true;
      if (useDeltaRforFootprint_)
        for (const auto& it : vetoedPtrs)
          if (it.isNonnull() && it.isAvailable() && reco::deltaR2(*it, *c) < 0.00000025) {
        addcand = false;
        break;
      }
      if (addcand)
        result->push_back(c);
    }
  }
  iEvent.put(std::move(result));
}

void CandPtrProjector::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("src");
  desc.add<edm::InputTag>("veto");
  desc.add<bool>("useDeltaRforFootprint", false);
  descriptions.addWithDefaultLabel(desc);
}

DEFINE_FWK_MODULE(CandPtrProjector);