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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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
bool	correctlyExcludeOverlap_

std::string	moduleLabel_

edm::EDGetTokenT< reco::PFCandidateCollection >	srcPFCandidates_

edm::EDGetTokenT< JetView >	srcSubjets_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 43 of file BoostedTauSeedsProducer.cc.

Member Typedef Documentation

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.

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

private

Definition at line 55 of file BoostedTauSeedsProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 63 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"));
 
   correctlyExcludeOverlap_ = (cfg.exists("correctlyExcludeOverlap")) ? cfg.getParameter<bool>("correctlyExcludeOverlap") : false;
   if (!correctlyExcludeOverlap_) {
     edm::LogWarning("") << "The \"correctlyExcludeOverlap\" flag set to \"False\":\n"
       << "=> module is working in a buggy backward compatibility mode; is it intended?";
   }
 
   verbosity_ = (cfg.exists("verbosity")) ? cfg.getParameter<int>("verbosity") : 0;
 
   produces<reco::PFJetCollection>();
   produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsolation");
   //produces<JetToPFCandidateAssociation>("pfCandAssocMapForIsoDepositVetos");
 }

BoostedTauSeedsProducer::~BoostedTauSeedsProducer ( )

inlineoverride

Definition at line 46 of file BoostedTauSeedsProducer.cc.

References produce().

46 {}

Member Function Documentation

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

override

Definition at line 194 of file BoostedTauSeedsProducer.cc.

References edm::View< T >::at(), correctlyExcludeOverlap_, gather_cfg::cout, DEFINE_FWK_MODULE, reco::LeafCandidate::eta(), edm::Event::getByToken(), edm::Event::getRefBeforePut(), mps_fire::i, training_settings::idx, reco::Jet::jetArea(), crabWrapper::key, reco::LeafCandidate::mass(), moduleLabel_, eostools::move(), reco::Jet::nConstituents(), patMuonMerger_cfi::pfCandidates, reco::LeafCandidate::phi(), edm::Event::productGetter(), reco::LeafCandidate::pt(), edm::Event::put(), edm::View< T >::size(), srcPFCandidates_, srcSubjets_, and verbosity_.

Referenced by ~BoostedTauSeedsProducer().

                                                                             {
   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, pfCandidatesNotInSubjet2;
     //Overlapping constituents can be incorrectly selected as in previous
     //buggy version of this module to fulfill non-changing policy in
     //10_6_X release series.
     if (correctlyExcludeOverlap_) {
       pfCandidatesNotInSubjet1 = getPFCandidates_exclJetConstituents(*subjet1, pfCandidates, constitmap[2 * idx + 1], false);
       pfCandidatesNotInSubjet2 = getPFCandidates_exclJetConstituents(*subjet2, pfCandidates, constitmap[2 * idx], false);
     }
     else {
       pfCandidatesNotInSubjet1 = getPFCandidates_exclJetConstituents(*subjet1, pfCandidates, constitmap[2 * idx], false);
       pfCandidatesNotInSubjet2 = getPFCandidates_exclJetConstituents(*subjet2, pfCandidates, constitmap[2 * idx + 1], 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");
 }