L1TEGMultiMerger.cc

Go to the documentation of this file.

#include "DataFormats/L1TCorrelator/interface/TkElectron.h"
#include "DataFormats/L1TCorrelator/interface/TkElectronFwd.h"
#include "DataFormats/L1Trigger/interface/EGamma.h"
#include "DataFormats/L1TCorrelator/interface/TkEm.h"
#include "DataFormats/L1TCorrelator/interface/TkEmFwd.h"

#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/transform.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include <vector>
#include <iostream>

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

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

  struct RefRemapper {
    BXVector<edm::Ref<BXVector<l1t::EGamma>>> oldRefs;
    std::map<edm::Ref<BXVector<l1t::EGamma>>, edm::Ref<BXVector<l1t::EGamma>>> old2newRefMap;
  };

  template <class T>
  class InstanceMerger {
  public:
    InstanceMerger(L1TEGMultiMerger* prod, const edm::ParameterSet& conf)
        : instanceLabel_(conf.getParameter<std::string>("instance")) {
      for (const auto& producer_tag : conf.getParameter<std::vector<edm::InputTag>>("pfProducers")) {
        tokens_.push_back(
            prod->consumes<T>(edm::InputTag(producer_tag.label(), producer_tag.instance(), producer_tag.process())));
      }
      // FIXME: move this outside
      prod->produces<T>(instanceLabel_);
    }

    void produce(edm::Event& iEvent, RefRemapper& refRemapper) const {
      edm::Handle<T> handle;
      auto out = std::make_unique<T>();
      for (const auto& token : tokens_) {
        iEvent.getByToken(token, handle);
        populate(out, handle, refRemapper);
      }
      remapRefs(iEvent, out, refRemapper);
      iEvent.put(std::move(out), instanceLabel_);
    }

  private:
    template <class TT>
    void remapRefs(edm::Event& iEvent, std::unique_ptr<TT>& out, RefRemapper& refRemapper) const {
      for (auto& egobj : *out) {
        auto newref = refRemapper.old2newRefMap.find(egobj.EGRef());
        if (newref != refRemapper.old2newRefMap.end()) {
          egobj.setEGRef(newref->second);
        }
      }
    }

    void remapRefs(edm::Event& iEvent, std::unique_ptr<BXVector<l1t::EGamma>>& out, RefRemapper& refRemapper) const {
      edm::RefProd<BXVector<l1t::EGamma>> ref_egs = iEvent.getRefBeforePut<BXVector<l1t::EGamma>>(instanceLabel_);
      edm::Ref<BXVector<l1t::EGamma>>::key_type idx = 0;
      for (std::size_t ix = 0; ix < out->size(); ix++) {
        refRemapper.old2newRefMap[refRemapper.oldRefs[ix]] = edm::Ref<BXVector<l1t::EGamma>>(ref_egs, idx++);
      }
    }

    template <class TT>
    void populate(std::unique_ptr<TT>& out, const edm::Handle<TT>& in, RefRemapper& refRemapper) const {
      out->insert(out->end(), in->begin(), in->end());
    }

    void populate(std::unique_ptr<BXVector<l1t::EGamma>>& out,
                  const edm::Handle<BXVector<l1t::EGamma>>& in,
                  RefRemapper& refRemapper) const {
      edm::Ref<BXVector<l1t::EGamma>>::key_type idx = 0;
      for (int bx = in->getFirstBX(); bx <= in->getLastBX(); bx++) {
        for (auto egee_itr = in->begin(bx); egee_itr != in->end(bx); egee_itr++) {
          out->push_back(bx, *egee_itr);
          // this to ensure that the old ref and the new object have the same index in the BXVector collection so that we can still match them
          // no matter which BX we will insert next
          refRemapper.oldRefs.push_back(bx, edm::Ref<BXVector<l1t::EGamma>>(in, idx++));
        }
      }
    }

    std::vector<edm::EDGetTokenT<T>> tokens_;
    std::string instanceLabel_;
  };

  std::vector<InstanceMerger<l1t::TkElectronCollection>> tkEleMerger;
  std::vector<InstanceMerger<l1t::TkEmCollection>> tkEmMerger;
  std::vector<InstanceMerger<BXVector<l1t::EGamma>>> tkEGMerger;
};

L1TEGMultiMerger::L1TEGMultiMerger(const edm::ParameterSet& conf) {
  for (const auto& config : conf.getParameter<std::vector<edm::ParameterSet>>("tkEgs")) {
    tkEGMerger.push_back(InstanceMerger<BXVector<l1t::EGamma>>(this, config));
  }
  for (const auto& config : conf.getParameter<std::vector<edm::ParameterSet>>("tkElectrons")) {
    tkEleMerger.push_back(InstanceMerger<l1t::TkElectronCollection>(this, config));
  }
  for (const auto& config : conf.getParameter<std::vector<edm::ParameterSet>>("tkEms")) {
    tkEmMerger.push_back(InstanceMerger<l1t::TkEmCollection>(this, config));
  }
}

L1TEGMultiMerger::~L1TEGMultiMerger() {}

void L1TEGMultiMerger::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup&) const {
  RefRemapper refmapper;
  for (const auto& egMerger : tkEGMerger)
    egMerger.produce(iEvent, refmapper);
  for (const auto& eleMerger : tkEleMerger)
    eleMerger.produce(iEvent, refmapper);
  for (const auto& emMerger : tkEmMerger)
    emMerger.produce(iEvent, refmapper);
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(L1TEGMultiMerger);