GenericBoostedTauSeedsProducer< JetType, CandType > Class Template Reference

Inheritance diagram for GenericBoostedTauSeedsProducer< JetType, CandType >:

Public Member Functions
	GenericBoostedTauSeedsProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~GenericBoostedTauSeedsProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Types
typedef std::vector< CandType >	CandTypeCollection
typedef edm::AssociationMap < edm::OneToMany< std::vector < JetType >, std::vector < CandType >, unsigned int > >	JetToPFCandidateAssociation
typedef std::vector< JetType >	JetTypeCollection
typedef edm::View< reco::Jet >	JetView

Private Attributes
std::string	moduleLabel_
edm::EDGetTokenT < CandTypeCollection >	srcPFCandidates_
edm::EDGetTokenT< JetView >	srcSubjets_
int	verbosity_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

template<class JetType, class CandType>
class GenericBoostedTauSeedsProducer< JetType, CandType >

Definition at line 46 of file BoostedTauSeedsProducer.cc.

Member Typedef Documentation

template<class JetType , class CandType >

typedef std::vector<CandType> GenericBoostedTauSeedsProducer< JetType, CandType >::CandTypeCollection

private

Definition at line 56 of file BoostedTauSeedsProducer.cc.

template<class JetType , class CandType >

typedef edm::AssociationMap<edm::OneToMany<std::vector<JetType>, std::vector<CandType>, unsigned int> > GenericBoostedTauSeedsProducer< JetType, CandType >::JetToPFCandidateAssociation

private

Definition at line 54 of file BoostedTauSeedsProducer.cc.

template<class JetType , class CandType >

typedef std::vector<JetType> GenericBoostedTauSeedsProducer< JetType, CandType >::JetTypeCollection

private

Definition at line 55 of file BoostedTauSeedsProducer.cc.

template<class JetType , class CandType >

typedef edm::View<reco::Jet> GenericBoostedTauSeedsProducer< JetType, CandType >::JetView

private

Definition at line 60 of file BoostedTauSeedsProducer.cc.

Constructor & Destructor Documentation

template<class JetType , class CandType >

GenericBoostedTauSeedsProducer< JetType, CandType >::GenericBoostedTauSeedsProducer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 68 of file BoostedTauSeedsProducer.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), GenericBoostedTauSeedsProducer< JetType, CandType >::srcPFCandidates_, GenericBoostedTauSeedsProducer< JetType, CandType >::srcSubjets_, and GenericBoostedTauSeedsProducer< JetType, CandType >::verbosity_.

    : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
  srcSubjets_ = consumes<JetView>(cfg.getParameter<edm::InputTag>("subjetSrc"));
  srcPFCandidates_ = consumes<CandTypeCollection>(cfg.getParameter<edm::InputTag>("pfCandidateSrc"));

  verbosity_ = (cfg.exists("verbosity")) ? cfg.getParameter<int>("verbosity") : 0;

  produces<JetTypeCollection>();
  produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsolation");
  //produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsoDepositVetos");
}

template<class JetType , class CandType >

GenericBoostedTauSeedsProducer< JetType, CandType >::~GenericBoostedTauSeedsProducer ( )

inlineoverride

Definition at line 49 of file BoostedTauSeedsProducer.cc.

{}

Member Function Documentation

template<class JetType , class CandType >

void GenericBoostedTauSeedsProducer< JetType, CandType >::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 195 of file BoostedTauSeedsProducer.cc.

References cms::cuda::assert(), gather_cfg::cout, reco::LeafCandidate::eta(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), mps_fire::i, dqmiolumiharvest::j, reco::Jet::jetArea(), submitPVResolutionJobs::key, reco::LeafCandidate::mass(), moduleLabel_(), eostools::move(), reco::Jet::nConstituents(), slimmedMuons_cfi::pfCandidates, reco::LeafCandidate::phi(), edm::Event::productGetter(), reco::LeafCandidate::pt(), and edm::Event::put().

                                                                                                      {
  if (verbosity_ >= 1) {
    std::cout << "<BoostedTauSeedsProducer::produce (moduleLabel = " << moduleLabel_ << ")>:" << std::endl;
  }

  edm::Handle<JetView> subjets;
  evt.getByToken(srcSubjets_, subjets);
  if (verbosity_ >= 1) {
    std::cout << "#subjets = " << subjets->size() << std::endl;
  }
  assert((subjets->size() % 2) == 0);  // CV: ensure that subjets come in pairs

  edm::Handle<CandTypeCollection> pfCandidates;
  evt.getByToken(srcPFCandidates_, pfCandidates);
  if (verbosity_ >= 1) {
    std::cout << "#pfCandidates = " << pfCandidates->size() << std::endl;
  }

  auto selectedSubjets = std::make_unique<JetTypeCollection>();
  edm::RefProd<JetTypeCollection> selectedSubjetRefProd = evt.getRefBeforePut<JetTypeCollection>();

  auto selectedSubjetPFCandidateAssociationForIsolation =
      std::make_unique<JetToPFCandidateAssociation>(&evt.productGetter());
  //auto selectedSubjetPFCandidateAssociationForIsoDepositVetos = std::make_unique<JetToPFCandidateAssociation>(&evt.productGetter());

  // cache for jet->pfcandidate
  JetToConstitMap constitmap(subjets->size());

  // fill constituents map
  const auto& thesubjets = *subjets;
  for (unsigned i = 0; i < thesubjets.size(); ++i) {
    for (unsigned j = 0; j < thesubjets[i].numberOfDaughters(); ++j) {
      constitmap[i].emplace(thesubjets[i].daughterPtr(j).key());
    }
  }

  for (size_t idx = 0; idx < (subjets->size() / 2); ++idx) {
    const reco::Jet* subjet1 = &subjets->at(2 * idx);
    const reco::Jet* subjet2 = &subjets->at(2 * idx + 1);
    assert(subjet1 && subjet2);
    if (verbosity_ >= 1) {
      std::cout << "processing jet #" << idx << ":" << std::endl;
      std::cout << " subjet1: Pt = " << subjet1->pt() << ", eta = " << subjet1->eta() << ", phi = " << subjet1->phi()
                << ", mass = " << subjet1->mass() << " (#constituents = " << subjet1->nConstituents()
                << ", area = " << subjet1->jetArea() << ")" << std::endl;
      std::cout << " subjet2: Pt = " << subjet2->pt() << ", eta = " << subjet2->eta() << ", phi = " << subjet2->phi()
                << ", mass = " << subjet2->mass() << " (#constituents = " << subjet2->nConstituents()
                << ", area = " << subjet2->jetArea() << ")" << std::endl;
    }

    if (!(subjet1->nConstituents() >= 1 && subjet1->pt() > 1. && subjet2->nConstituents() >= 1 && subjet2->pt() > 1.))
      continue;  // CV: skip pathological cases

    // find PFCandidate constituents of each subjet
    reco::Jet::Constituents subjetConstituents1;
    getJetConstituents(*subjet1, subjetConstituents1);
    reco::Jet::Constituents subjetConstituents2;
    getJetConstituents(*subjet2, subjetConstituents2);
    if (verbosity_ >= 1) {
      printJetConstituents("subjetConstituents1", subjetConstituents1);
      printJetConstituents("subjetConstituents2", subjetConstituents2);
    }

    selectedSubjets->push_back(convertToPFJet<JetType, CandType>(*subjet1, subjetConstituents1));
    edm::Ref<JetTypeCollection> subjetRef1(selectedSubjetRefProd, selectedSubjets->size() - 1);
    selectedSubjets->push_back(convertToPFJet<JetType, CandType>(*subjet2, subjetConstituents2));
    edm::Ref<JetTypeCollection> subjetRef2(selectedSubjetRefProd, selectedSubjets->size() - 1);

    // find all PFCandidates that are not constituents of the **other** subjet
    std::vector<edm::Ref<std::vector<CandType>>> pfCandidatesNotInSubjet1 =
        getPFCandidates_exclJetConstituents<CandType>(*subjet1, pfCandidates, constitmap[2 * idx + 1], false);
    std::vector<edm::Ref<std::vector<CandType>>> pfCandidatesNotInSubjet2 =
        getPFCandidates_exclJetConstituents<CandType>(*subjet2, pfCandidates, constitmap[2 * idx], false);
    if (verbosity_ >= 1) {
      std::cout << "#pfCandidatesNotInSubjet1 = " << pfCandidatesNotInSubjet1.size() << std::endl;
      std::cout << "#pfCandidatesNotInSubjet2 = " << pfCandidatesNotInSubjet2.size() << std::endl;
    }

    // build JetToPFCandidateAssociation
    // (key = subjet, value = collection of PFCandidates that are not constituents of subjet)
    for (auto const& pfCandidate : pfCandidatesNotInSubjet1) {
      selectedSubjetPFCandidateAssociationForIsolation->insert(subjetRef1, pfCandidate);
    }
    for (auto const& pfCandidate : pfCandidatesNotInSubjet2) {
      selectedSubjetPFCandidateAssociationForIsolation->insert(subjetRef2, pfCandidate);
    }
  }

  evt.put(std::move(selectedSubjets));
  evt.put(std::move(selectedSubjetPFCandidateAssociationForIsolation), "pfCandAssocMapForIsolation");
}

Member Data Documentation

template<class JetType , class CandType >

std::string GenericBoostedTauSeedsProducer< JetType, CandType >::moduleLabel_

private

Definition at line 58 of file BoostedTauSeedsProducer.cc.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

template<class JetType , class CandType >

edm::EDGetTokenT<CandTypeCollection> GenericBoostedTauSeedsProducer< JetType, CandType >::srcPFCandidates_

private

Definition at line 62 of file BoostedTauSeedsProducer.cc.

Referenced by GenericBoostedTauSeedsProducer< JetType, CandType >::GenericBoostedTauSeedsProducer().

template<class JetType , class CandType >

edm::EDGetTokenT<JetView> GenericBoostedTauSeedsProducer< JetType, CandType >::srcSubjets_

private

Definition at line 61 of file BoostedTauSeedsProducer.cc.

Referenced by GenericBoostedTauSeedsProducer< JetType, CandType >::GenericBoostedTauSeedsProducer().

template<class JetType , class CandType >

int GenericBoostedTauSeedsProducer< JetType, CandType >::verbosity_

private

Definition at line 64 of file BoostedTauSeedsProducer.cc.

Referenced by GenericBoostedTauSeedsProducer< JetType, CandType >::GenericBoostedTauSeedsProducer().