MultiIsolator.cc

Go to the documentation of this file.

#include "PhysicsTools/PatAlgos/interface/MultiIsolator.h"
#include "PhysicsTools/PatAlgos/interface/SimpleIsolator.h"
#include "PhysicsTools/PatAlgos/interface/IsoDepositIsolator.h"
#include "DataFormats/PatCandidates/interface/Flags.h"
#include <sstream>

using namespace pat::helper;

MultiIsolator::MultiIsolator(const edm::ParameterSet &conf, edm::ConsumesCollector &&iC, bool cuts) {
  using pat::Flags;
  if (conf.exists("tracker"))
    addIsolator(conf.getParameter<edm::ParameterSet>("tracker"), iC, cuts, Flags::Isolation::Tracker, pat::TrackIso);
  if (conf.exists("ecal"))
    addIsolator(conf.getParameter<edm::ParameterSet>("ecal"), iC, cuts, Flags::Isolation::ECal, pat::EcalIso);
  if (conf.exists("hcal"))
    addIsolator(conf.getParameter<edm::ParameterSet>("hcal"), iC, cuts, Flags::Isolation::HCal, pat::HcalIso);
  if (conf.exists("calo"))
    addIsolator(conf.getParameter<edm::ParameterSet>("calo"), iC, cuts, Flags::Isolation::Calo, pat::CaloIso);
  if (conf.exists("calo") && (conf.exists("ecal") || conf.exists("hcal"))) {
    throw cms::Exception("Configuration")
        << "MultiIsolator: you can't specify both 'calo' isolation and 'ecal'/'hcal', "
        << "as the 'calo' isolation flag is just the logical OR of 'ecal' and 'hcal'.\n";
  }
  if (conf.exists("pfAllParticles"))
    addIsolator(
        conf.getParameter<edm::ParameterSet>("pfAllParticles"), iC, cuts, Flags::Isolation::Calo, pat::PfAllParticleIso);
  if (conf.exists("pfChargedHadron"))
    addIsolator(conf.getParameter<edm::ParameterSet>("pfChargedHadron"),
                iC,
                cuts,
                Flags::Isolation::Calo,
                pat::PfChargedHadronIso);
  if (conf.exists("pfNeutralHadron"))
    addIsolator(conf.getParameter<edm::ParameterSet>("pfNeutralHadron"),
                iC,
                cuts,
                Flags::Isolation::Calo,
                pat::PfNeutralHadronIso);
  if (conf.exists("pfGamma"))
    addIsolator(conf.getParameter<edm::ParameterSet>("pfGamma"), iC, cuts, Flags::Isolation::Calo, pat::PfGammaIso);
  if (conf.exists("user")) {
    std::vector<edm::ParameterSet> psets = conf.getParameter<std::vector<edm::ParameterSet> >("user");
    if (psets.size() > 5) {
      throw cms::Exception("Configuration") << "MultiIsolator: you can specify at most 5 user isolation collections.\n";
    }
    uint32_t bit = Flags::Isolation::User1;
    for (std::vector<edm::ParameterSet>::const_iterator it = psets.begin(), ed = psets.end(); it != ed;
         ++it, bit <<= 1) {
      addIsolator(*it, iC, cuts, bit, pat::IsolationKeys(pat::UserBaseIso + (it - psets.begin())));
    }
  }
}

void MultiIsolator::addIsolator(BaseIsolator *iso, uint32_t mask, pat::IsolationKeys key) {
  isolators_.emplace_back(iso);
  masks_.push_back(mask);
  keys_.push_back(key);
}

BaseIsolator *MultiIsolator::make(const edm::ParameterSet &conf, edm::ConsumesCollector &iC, bool withCut) {
  if (conf.empty())
    return nullptr;
  if (conf.exists("placeholder") && conf.getParameter<bool>("placeholder"))
    return nullptr;
  if (conf.exists("deltaR")) {
    return new IsoDepositIsolator(conf, iC, withCut);
  } else {
    return new SimpleIsolator(conf, iC, withCut);
  }
}

void MultiIsolator::addIsolator(
    const edm::ParameterSet &conf, edm::ConsumesCollector &iC, bool withCut, uint32_t mask, pat::IsolationKeys key) {
  BaseIsolator *iso = make(conf, iC, withCut);
  if (iso)
    addIsolator(iso, mask, key);
}

void MultiIsolator::beginEvent(const edm::Event &event, const edm::EventSetup &eventSetup) {
  for (auto it = isolators_.begin(), ed = isolators_.end(); it != ed; ++it) {
    (*it)->beginEvent(event, eventSetup);
  }
}

void MultiIsolator::endEvent() {
  for (auto it = isolators_.begin(), ed = isolators_.end(); it != ed; ++it) {
    (*it)->endEvent();
  }
}

void MultiIsolator::print(std::ostream &out) const {
  for (auto it = isolators_.cbegin(), ed = isolators_.end(); it != ed; ++it) {
    out << " * ";
    (*it)->print(out);
    out << ": Flag " << pat::Flags::bitToString(masks_[it - isolators_.begin()]) << "\n";
  }
  out << "\n";
}

std::string MultiIsolator::printSummary() const {
  std::ostringstream isoSumm;
  print(isoSumm);
  return isoSumm.str();
}