BoostedTauSeedsProducer Class Reference

Inheritance diagram for BoostedTauSeedsProducer:

Public Member Functions
	BoostedTauSeedsProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~BoostedTauSeedsProducer () 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 edm::AssociationMap< edm::OneToMany< std::vector< reco::PFJet >, std::vector< reco::PFCandidate >, unsigned int > >	JetToPFCandidateAssociation
typedef edm::View< reco::Jet >	JetView

Private Attributes
std::string	moduleLabel_
edm::EDGetTokenT< reco::PFCandidateCollection >	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

Definition at line 43 of file BoostedTauSeedsProducer.cc.

Member Typedef Documentation

JetToPFCandidateAssociation

typedef edm::AssociationMap<edm::OneToMany<std::vector<reco::PFJet>, std::vector<reco::PFCandidate>, unsigned int> > BoostedTauSeedsProducer::JetToPFCandidateAssociation

private

Definition at line 51 of file BoostedTauSeedsProducer.cc.

JetView

typedef edm::View<reco::Jet> BoostedTauSeedsProducer::JetView

private

Definition at line 55 of file BoostedTauSeedsProducer.cc.

Constructor & Destructor Documentation

BoostedTauSeedsProducer()

BoostedTauSeedsProducer::BoostedTauSeedsProducer ( const edm::ParameterSet &  cfg )

explicit

Definition at line 62 of file BoostedTauSeedsProducer.cc.

    : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
  srcSubjets_ = consumes<JetView>(cfg.getParameter<edm::InputTag>("subjetSrc"));
  srcPFCandidates_ = consumes<reco::PFCandidateCollection>(cfg.getParameter<edm::InputTag>("pfCandidateSrc"));
 
  verbosity_ = (cfg.exists("verbosity")) ? cfg.getParameter<int>("verbosity") : 0;
 
  produces<reco::PFJetCollection>();
  produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsolation");
  //produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsoDepositVetos");
}

References looper::cfg, srcPFCandidates_, srcSubjets_, and verbosity_.

~BoostedTauSeedsProducer()

BoostedTauSeedsProducer::~BoostedTauSeedsProducer ( )

inlineoverride

Definition at line 46 of file BoostedTauSeedsProducer.cc.

{}

Member Function Documentation

produce()

void BoostedTauSeedsProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 187 of file BoostedTauSeedsProducer.cc.

                                                                            {
  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<reco::PFCandidateCollection> pfCandidates;
  evt.getByToken(srcPFCandidates_, pfCandidates);
  if (verbosity_ >= 1) {
    std::cout << "#pfCandidates = " << pfCandidates->size() << std::endl;
  }
 
  auto selectedSubjets = std::make_unique<reco::PFJetCollection>();
  edm::RefProd<reco::PFJetCollection> selectedSubjetRefProd = evt.getRefBeforePut<reco::PFJetCollection>();
 
  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(*subjet1, subjetConstituents1));
    edm::Ref<reco::PFJetCollection> subjetRef1(selectedSubjetRefProd, selectedSubjets->size() - 1);
    selectedSubjets->push_back(convertToPFJet(*subjet2, subjetConstituents2));
    edm::Ref<reco::PFJetCollection> subjetRef2(selectedSubjetRefProd, selectedSubjets->size() - 1);
 
    // find all PFCandidates that are not constituents of the **other** subjet
    std::vector<reco::PFCandidateRef> pfCandidatesNotInSubjet1 =
        getPFCandidates_exclJetConstituents(*subjet1, pfCandidates, constitmap[2 * idx + 1], false);
    std::vector<reco::PFCandidateRef> pfCandidatesNotInSubjet2 =
        getPFCandidates_exclJetConstituents(*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");
}

References cms::cuda::assert(), edm::View< T >::at(), gather_cfg::cout, reco::LeafCandidate::eta(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), mps_fire::i, heavyIonCSV_trainingSettings::idx, dqmiolumiharvest::j, reco::Jet::jetArea(), crabWrapper::key, reco::LeafCandidate::mass(), moduleLabel_, eostools::move(), reco::Jet::nConstituents(), zmumugammaAnalyzer_cfi::pfCandidates, reco::LeafCandidate::phi(), edm::Event::productGetter(), reco::LeafCandidate::pt(), edm::Event::put(), edm::View< T >::size(), srcPFCandidates_, srcSubjets_, and verbosity_.

Member Data Documentation

moduleLabel_

std::string BoostedTauSeedsProducer::moduleLabel_

private

Definition at line 53 of file BoostedTauSeedsProducer.cc.

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

srcPFCandidates_

edm::EDGetTokenT<reco::PFCandidateCollection> BoostedTauSeedsProducer::srcPFCandidates_

private

Definition at line 57 of file BoostedTauSeedsProducer.cc.

Referenced by BoostedTauSeedsProducer(), and produce().

srcSubjets_

edm::EDGetTokenT<JetView> BoostedTauSeedsProducer::srcSubjets_

private

Definition at line 56 of file BoostedTauSeedsProducer.cc.

Referenced by BoostedTauSeedsProducer(), and produce().

verbosity_

int BoostedTauSeedsProducer::verbosity_

private

Definition at line 59 of file BoostedTauSeedsProducer.cc.

Referenced by BoostedTauSeedsProducer(), and produce().