IsolatedPFCandidateSelectorDefinition.h

Go to the documentation of this file.

#ifndef CommonTools_ParticleFlow_IsolatedPFCandidateSelectorDefinition
#define CommonTools_ParticleFlow_IsolatedPFCandidateSelectorDefinition

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "CommonTools/ParticleFlow/interface/PFCandidateSelectorDefinition.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Utilities/interface/transform.h"

namespace pf2pat {

  class IsolatedPFCandidateSelectorDefinition : public PFCandidateSelectorDefinition {
  public:
    typedef edm::ValueMap<double> IsoMap;

    IsolatedPFCandidateSelectorDefinition(const edm::ParameterSet& cfg, edm::ConsumesCollector&& iC)
        : isolationValueMapChargedTokens_(
              edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("isolationValueMapsCharged"),
                                    [&](edm::InputTag const& tag) { return iC.consumes<IsoMap>(tag); })),
          isolationValueMapNeutralTokens_(
              edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("isolationValueMapsNeutral"),
                                    [&](edm::InputTag const& tag) { return iC.consumes<IsoMap>(tag); })),
          doDeltaBetaCorrection_(cfg.getParameter<bool>("doDeltaBetaCorrection")),
          deltaBetaIsolationValueMapToken_(
              iC.mayConsume<IsoMap>(cfg.getParameter<edm::InputTag>("deltaBetaIsolationValueMap"))),
          deltaBetaFactor_(cfg.getParameter<double>("deltaBetaFactor")),
          isRelative_(cfg.getParameter<bool>("isRelative")),
          isolationCut_(cfg.getParameter<double>("isolationCut")) {}

    void select(const HandleToCollection& hc, const edm::Event& e, const edm::EventSetup& s) {
      selected_.clear();

      // read all charged isolation value maps
      std::vector<edm::Handle<IsoMap> > isoMapsCharged(isolationValueMapChargedTokens_.size());
      for (unsigned iMap = 0; iMap < isolationValueMapChargedTokens_.size(); ++iMap) {
        e.getByToken(isolationValueMapChargedTokens_[iMap], isoMapsCharged[iMap]);
      }

      // read all neutral isolation value maps
      std::vector<edm::Handle<IsoMap> > isoMapsNeutral(isolationValueMapNeutralTokens_.size());
      for (unsigned iMap = 0; iMap < isolationValueMapNeutralTokens_.size(); ++iMap) {
        e.getByToken(isolationValueMapNeutralTokens_[iMap], isoMapsNeutral[iMap]);
      }

      edm::Handle<IsoMap> dBetaH;
      if (doDeltaBetaCorrection_) {
        e.getByToken(deltaBetaIsolationValueMapToken_, dBetaH);
      }

      unsigned key = 0;
      for (collection::const_iterator pfc = hc->begin(); pfc != hc->end(); ++pfc, ++key) {
        reco::PFCandidateRef candidate(hc, key);

        bool passed = true;
        double isoSumCharged = 0.0;
        double isoSumNeutral = 0.0;

        for (unsigned iMap = 0; iMap < isoMapsCharged.size(); ++iMap) {
          const IsoMap& isoMap = *(isoMapsCharged[iMap]);
          double val = isoMap[candidate];
          isoSumCharged += val;
        }

        for (unsigned iMap = 0; iMap < isoMapsNeutral.size(); ++iMap) {
          const IsoMap& isoMap = *(isoMapsNeutral[iMap]);
          double val = isoMap[candidate];
          isoSumNeutral += val;
        }

        if (doDeltaBetaCorrection_) {
          const IsoMap& isoMap = *dBetaH;
          double dBetaVal = isoMap[candidate];
          double dBetaCorIsoSumNeutral = isoSumNeutral + deltaBetaFactor_ * dBetaVal;
          isoSumNeutral = dBetaCorIsoSumNeutral > 0 ? dBetaCorIsoSumNeutral : 0;  //follow muon POG definition in 2012
        }

        double isoSum = isoSumCharged + isoSumNeutral;

        if (isRelative_) {
          isoSum /= candidate->pt();
        }

        if (isoSum > isolationCut_) {
          passed = false;
        }

        if (passed) {
          // passed all cuts, selected
          selected_.push_back(reco::PFCandidate(*pfc));
          reco::PFCandidatePtr ptrToMother(hc, key);
          selected_.back().setSourceCandidatePtr(ptrToMother);
        }
      }
    }

  private:
    std::vector<edm::EDGetTokenT<IsoMap> > isolationValueMapChargedTokens_;
    std::vector<edm::EDGetTokenT<IsoMap> > isolationValueMapNeutralTokens_;
    bool doDeltaBetaCorrection_;
    edm::EDGetTokenT<IsoMap> deltaBetaIsolationValueMapToken_;
    double deltaBetaFactor_;
    bool isRelative_;
    double isolationCut_;
  };

}  // namespace pf2pat

#endif