OverlapExclusionSelector.h

Go to the documentation of this file.

#ifndef CommonTools_UtilAlgos_OverlapExclusionSelector_h
#define CommonTools_UtilAlgos_OverlapExclusionSelector_h
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "DataFormats/Common/interface/Handle.h"

namespace edm { class EventSetup; }

template<typename C, typename T, typename O>
class OverlapExclusionSelector {
public:
  OverlapExclusionSelector(const edm::ParameterSet& cfg, edm::ConsumesCollector && iC) :
    OverlapExclusionSelector(cfg, iC) {};
  OverlapExclusionSelector(const edm::ParameterSet& cfg, edm::ConsumesCollector & iC);
  void newEvent(const edm::Event&, const edm::EventSetup&) const;
  bool operator()(const T&) const;
private:
  edm::EDGetTokenT<C> srcToken_;
  mutable typename C::const_iterator begin_, end_;
  O overlap_;
};

template<typename C, typename T, typename O>
OverlapExclusionSelector<C, T, O>::OverlapExclusionSelector(const edm::ParameterSet& cfg, edm::ConsumesCollector & iC) :
  srcToken_(iC.consumes<C>(cfg.template getParameter<edm::InputTag>("overlap"))),
  overlap_(cfg) {
}

template<typename C, typename T, typename O>
       void OverlapExclusionSelector<C, T, O>::newEvent(const edm::Event& evt, const edm::EventSetup&) const {
  edm::Handle<C> h;
  evt.getByToken(srcToken_, h);
  begin_ = h->begin();
  end_ = h->end();
}

template<typename C, typename T, typename O>
bool OverlapExclusionSelector<C, T, O>::operator()(const T& t) const {
  bool noOverlap = true;
  for(typename C::const_iterator i = begin_; i != end_; ++i) {
    if(overlap_(*i, t)) {
      noOverlap = false;
      break;
    }
  }
  return noOverlap;
}

#include "CommonTools/UtilAlgos/interface/EventSetupInitTrait.h"

EVENTSETUP_STD_INIT_T3(OverlapExclusionSelector);

#endif