Inheritance diagram for LeptonFSRProducer:

Public Member Functions
	LeptonFSRProducer (const edm::ParameterSet &iConfig)

	~LeptonFSRProducer () override=default

Public Member Functions inherited from edm::global::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

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
double	computeRelativeIsolation (const pat::PackedCandidate &photon, const pat::PackedCandidateCollection &pfcands, const double &isoConeMax2, const double &isoConeMin2) const

bool	electronFootprintVeto (pat::PackedCandidateRef &pfcandRef, edm::Handle< pat::ElectronCollection > electronsForVeto) const

void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

Private Attributes
const double	drEtCut

const double	drSafe

const edm::EDGetTokenT < edm::View< reco::GsfElectron > >	electrons_

const edm::EDGetTokenT < pat::ElectronCollection >	electronsForVeto_

const double	etaCutE

const double	etaCutMu

const double	isoCut

const edm::EDGetTokenT < edm::View< reco::Muon > >	muons_

const edm::EDGetTokenT < pat::PackedCandidateCollection >	pfcands_

const double	photonPtCut

const double	ptCutE

const double	ptCutMu

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Search for FSR photons for muons and electrons.

Photon candidates are searched among the "packedPFCandidates" collection with the specified cuts, and are required to be isolated (relIso, with a cone of 0.3) and not to be in the footprint of all electrons in the "electrons" collection. Each photon is matched by DeltaR to the closest among all muons and electrons and stored if passing dR/ET^2<deltaROverEt2Max. In addition ValueMaps are stored, with links to one photon per muon/electron. For this purpose, if more than a photon is matched to a lepton, the lowest-DR/ET^2 is chosen.

Definition at line 31 of file LeptonFSRProducer.cc.

Constructor & Destructor Documentation

LeptonFSRProducer::LeptonFSRProducer ( const edm::ParameterSet & iConfig )

inlineexplicit

Definition at line 33 of file LeptonFSRProducer.cc.

References edm::ParameterSet::getParameter().

       : pfcands_{consumes<pat::PackedCandidateCollection>(iConfig.getParameter<edm::InputTag>("packedPFCandidates"))},
         electronsForVeto_{consumes<pat::ElectronCollection>(iConfig.getParameter<edm::InputTag>("slimmedElectrons"))},
         muons_{consumes<edm::View<reco::Muon>>(iConfig.getParameter<edm::InputTag>("muons"))},
         electrons_{consumes<edm::View<reco::GsfElectron>>(iConfig.getParameter<edm::InputTag>("electrons"))},
         ptCutMu(iConfig.getParameter<double>("muonPtMin")),
         etaCutMu(iConfig.getParameter<double>("muonEtaMax")),
         ptCutE(iConfig.getParameter<double>("elePtMin")),
         etaCutE(iConfig.getParameter<double>("eleEtaMax")),
         photonPtCut(iConfig.getParameter<double>("photonPtMin")),
         drEtCut(iConfig.getParameter<double>("deltaROverEt2Max")),
         isoCut(iConfig.getParameter<double>("isolation")),
         drSafe(0.0001) {
     produces<std::vector<pat::GenericParticle>>();
     produces<edm::ValueMap<int>>("muFsrIndex");
     produces<edm::ValueMap<int>>("eleFsrIndex");
   }

LeptonFSRProducer::~LeptonFSRProducer ( )

overridedefault

Member Function Documentation

double LeptonFSRProducer::computeRelativeIsolation	(	const pat::PackedCandidate &	photon,
		const pat::PackedCandidateCollection &	pfcands,
		const double &	isoConeMax2,
		const double &	isoConeMin2
	)		const

private

Definition at line 239 of file LeptonFSRProducer.cc.

References funct::abs(), reco::deltaR2(), drSafe, pat::PackedCandidate::eta(), pat::PackedCandidate::phi(), and pat::PackedCandidate::pt().

Referenced by produce().

                                                                                     {
   double ptsum = 0;
 
   for (const auto& pfcand : pfcands) {
     // Isolation cone
     double dRIsoCone2 = deltaR2(photon.eta(), photon.phi(), pfcand.eta(), pfcand.phi());
     if (dRIsoCone2 > isoConeMax2 || dRIsoCone2 < isoConeMin2)
       continue;
 
     // Charged hadrons
     if (pfcand.charge() != 0 && abs(pfcand.pdgId()) == 211 && pfcand.pt() > 0.2 && dRIsoCone2 > drSafe * drSafe) {
       ptsum += pfcand.pt();
       // Neutral hadrons + photons
     } else if (pfcand.charge() == 0 && (abs(pfcand.pdgId()) == 22 || abs(pfcand.pdgId()) == 130) && pfcand.pt() > 0.5 &&
                dRIsoCone2 > 0.01 * 0.01) {
       ptsum += pfcand.pt();
     }
   }
 
   return ptsum / photon.pt();
 }

bool LeptonFSRProducer::electronFootprintVeto	(	pat::PackedCandidateRef &	pfcandRef,
		edm::Handle< pat::ElectronCollection >	electronsForVeto
	)		const

private

Definition at line 264 of file LeptonFSRProducer.cc.

References edm::Ref< C, T, F >::id(), and edm::Ref< C, T, F >::key().

Referenced by produce().

                                                                                                        {
   bool skipPhoton = false;
   for (auto electrons_iter = electronsForVeto->begin(); electrons_iter != electronsForVeto->end(); ++electrons_iter) {
     for (auto const& cand : electrons_iter->associatedPackedPFCandidates()) {
       if (!cand.isAvailable())
         continue;
       if (cand.id() != pfcandRef.id())
         throw cms::Exception("Configuration")
             << "The electron associatedPackedPFCandidates item does not have "
             << "the same ID of packed candidate collection used for cleaning the electron footprint: " << cand.id()
             << " (" << pfcandRef.id() << ")\n";
       if (cand.key() == pfcandRef.key()) {
         skipPhoton = true;
         break;
       }
     }
   }
   return skipPhoton;
 }

static void LeptonFSRProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

inlinestatic

Definition at line 50 of file LeptonFSRProducer.cc.

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and edm::ParameterDescriptionNode::setComment().

                                                                            {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag>("packedPFCandidates", edm::InputTag("packedPFCandidates"))
         ->setComment("packed pf candidates where to look for photons");
     desc.add<edm::InputTag>("slimmedElectrons", edm::InputTag("slimmedElectrons"))
         ->setComment(
             "electrons to check for footprint, the electron collection must have proper linking with the "
             "packedCandidate collection");
     desc.add<edm::InputTag>("muons", edm::InputTag("slimmedMuons"))
         ->setComment("collection of muons to match with FSR ");
     desc.add<edm::InputTag>("electrons", edm::InputTag("slimmedElectrons"))
         ->setComment("collection of electrons to match with FSR ");
     desc.add<double>("muonPtMin", 3.)->setComment("minimum pt of the muon to look for a near photon");
     desc.add<double>("muonEtaMax", 2.4)->setComment("max eta of the muon to look for a near photon");
     desc.add<double>("elePtMin", 5.)->setComment("minimum pt of the electron to look for a near photon");
     desc.add<double>("eleEtaMax", 2.5)->setComment("max eta of the electron to look for a near photon");
     desc.add<double>("photonPtMin", 2.0)->setComment("minimum photon Pt");
     desc.add<double>("deltaROverEt2Max", 0.05)->setComment("max ratio of deltaR(lep,photon) over et2 of the photon");
     desc.add<double>("isolation", 2.0)->setComment("photon relative isolation cut");
 
     descriptions.addWithDefaultLabel(desc);
   }

void LeptonFSRProducer::produce	(	edm::StreamID	streamID,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overrideprivatevirtual

Implements edm::global::EDProducerBase.

Definition at line 100 of file LeptonFSRProducer.cc.

References funct::abs(), computeRelativeIsolation(), HLT_FULL_cff::deltaR, HLT_FULL_cff::dRMin, drSafe, electronFootprintVeto(), HI_PhotonSkim_cff::electrons, electrons_, electronsForVeto_, etaCutE, etaCutMu, edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), edm::helper::Filler< Map >::insert(), isoCut, eostools::move(), HLT_FULL_cff::muon, patZpeak::muons, muons_, HLT_FULL_cff::pfcands, pfcands_, photonPtCut, DiDispStaMuonMonitor_cfi::pt, ptCutE, ptCutMu, and edm::Event::put().

                                                                                                          {
   using namespace std;
 
   edm::Handle<pat::PackedCandidateCollection> pfcands;
   iEvent.getByToken(pfcands_, pfcands);
   edm::Handle<edm::View<reco::Muon>> muons;
   iEvent.getByToken(muons_, muons);
   edm::Handle<edm::View<reco::GsfElectron>> electrons;
   iEvent.getByToken(electrons_, electrons);
   edm::Handle<pat::ElectronCollection> electronsForVeto;
   iEvent.getByToken(electronsForVeto_, electronsForVeto);
 
   // The output collection of FSR photons
   auto fsrPhotons = std::make_unique<std::vector<pat::GenericParticle>>();
 
   std::vector<int> muPhotonIdxs(muons->size(), -1);
   std::vector<double> muPhotonDRET2(muons->size(), 1e9);
 
   std::vector<int> elePhotonIdxs(electrons->size(), -1);
   std::vector<double> elePhotonDRET2(electrons->size(), 1e9);
 
   //----------------------
   // Loop on photon candidates
   //----------------------
 
   for (auto pc = pfcands->begin(); pc != pfcands->end(); pc++) {
     // consider only photons, with pT and eta cuts
     if (abs(pc->pdgId()) != 22 || pc->pt() < photonPtCut || abs(pc->eta()) > 2.5)
       continue;
 
     //------------------------------------------------------
     // Get the closest lepton
     //------------------------------------------------------
     double dRMin(0.5);
     int closestMu = -1;
     int closestEle = -1;
     double photon_relIso03 = 1e9;  // computed only if necessary
     bool skipPhoton = false;
 
     for (auto muon = muons->begin(); muon != muons->end(); ++muon) {
       if (muon->pt() < ptCutMu || std::abs(muon->eta()) > etaCutMu)
         continue;
 
       int muonIdx = muon - muons->begin();
       double dR = deltaR(muon->eta(), muon->phi(), pc->eta(), pc->phi());
       if (dR < dRMin && dR > drSafe && dR < drEtCut * pc->pt() * pc->pt()) {
         // Check if photon is isolated
         photon_relIso03 = computeRelativeIsolation(*pc, *pfcands, 0.3 * 0.3, drSafe * drSafe);
         if (photon_relIso03 > isoCut) {
           skipPhoton = true;
           break;  // break loop on muons -> photon will be skipped
         }
         // Check that photon is not in footprint of an electron
         pat::PackedCandidateRef pfcandRef = pat::PackedCandidateRef(pfcands, pc - pfcands->begin());
         skipPhoton = electronFootprintVeto(pfcandRef, electronsForVeto);
         if (skipPhoton)
           break;  // break loop on muons -> photon will be skipped
 
         // Candidate matching
         dRMin = dR;
         closestMu = muonIdx;
       }
     }  // end of loop on muons
 
     if (skipPhoton)
       continue;  // photon does not pass iso or ele footprint veto; do not look for electrons
 
     for (auto ele = electrons->begin(); ele != electrons->end(); ++ele) {
       if (ele->pt() < ptCutE || std::abs(ele->eta()) > etaCutE)
         continue;
 
       int eleIdx = ele - electrons->begin();
       double dR = deltaR(ele->eta(), ele->phi(), pc->eta(), pc->phi());
       if (dR < dRMin && dR > drSafe && dR < drEtCut * pc->pt() * pc->pt()) {
         // Check if photon is isolated (no need to recompute iso if already done for muons above)
         if (photon_relIso03 > 1e8) {
           photon_relIso03 = computeRelativeIsolation(*pc, *pfcands, 0.3 * 0.3, drSafe * drSafe);
         }
         if (photon_relIso03 > isoCut) {
           break;  // break loop on electrons -> photon will be skipped
         }
         // Check that photon is not in footprint of an electron
         pat::PackedCandidateRef pfcandRef = pat::PackedCandidateRef(pfcands, pc - pfcands->begin());
         if (electronFootprintVeto(pfcandRef, electronsForVeto)) {
           break;  // break loop on electrons -> photon will be skipped
         }
 
         // Candidate matching
         dRMin = dR;
         closestEle = eleIdx;
         closestMu = -1;  // reset match to muons
       }
     }  // end loop on electrons
 
     if (closestMu >= 0 || closestEle >= 0) {
       // Add FSR photon to the output collection
       double dRET2 = dRMin / pc->pt() / pc->pt();
       int iPhoton = fsrPhotons->size();
       fsrPhotons->push_back(pat::GenericParticle(*pc));
       fsrPhotons->back().addUserFloat("relIso03", photon_relIso03);
       fsrPhotons->back().addUserFloat("dROverEt2", dRET2);
 
       if (closestMu >= 0) {
         fsrPhotons->back().addUserCand("associatedMuon", reco::CandidatePtr(muons, closestMu));
         // Store the backlink to the photon: choose the lowest-dRET2 photon for each mu...
         if (dRET2 < muPhotonDRET2[closestMu]) {
           muPhotonDRET2[closestMu] = dRET2;
           muPhotonIdxs[closestMu] = iPhoton;
         }
       } else if (closestEle >= 0) {
         // ...and same for eles
         fsrPhotons->back().addUserCand("associatedElectron", reco::CandidatePtr(electrons, closestEle));
         if (dRET2 < elePhotonDRET2[closestEle]) {
           elePhotonDRET2[closestEle] = dRET2;
           elePhotonIdxs[closestEle] = iPhoton;
         }
       }
     }
   }  // end of loop over pfCands
 
   edm::OrphanHandle<std::vector<pat::GenericParticle>> oh = iEvent.put(std::move(fsrPhotons));
 
   {
     std::unique_ptr<edm::ValueMap<int>> bareIdx(new edm::ValueMap<int>());
     edm::ValueMap<int>::Filler fillerBareIdx(*bareIdx);
     fillerBareIdx.insert(muons, muPhotonIdxs.begin(), muPhotonIdxs.end());
     fillerBareIdx.fill();
     iEvent.put(std::move(bareIdx), "muFsrIndex");
   }
 
   {
     std::unique_ptr<edm::ValueMap<int>> bareIdx(new edm::ValueMap<int>());
     edm::ValueMap<int>::Filler fillerBareIdx(*bareIdx);
     fillerBareIdx.insert(electrons, elePhotonIdxs.begin(), elePhotonIdxs.end());
     fillerBareIdx.fill();
     iEvent.put(std::move(bareIdx), "eleFsrIndex");
   }
 }

Member Data Documentation

const double LeptonFSRProducer::drEtCut

private

Definition at line 95 of file LeptonFSRProducer.cc.

const double LeptonFSRProducer::drSafe

private

Definition at line 97 of file LeptonFSRProducer.cc.

Referenced by computeRelativeIsolation(), and produce().

const edm::EDGetTokenT<edm::View<reco::GsfElectron> > LeptonFSRProducer::electrons_

private

Definition at line 89 of file LeptonFSRProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<pat::ElectronCollection> LeptonFSRProducer::electronsForVeto_

private

Definition at line 87 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::etaCutE

private

Definition at line 93 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::etaCutMu

private

Definition at line 91 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::isoCut

private

Definition at line 96 of file LeptonFSRProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<edm::View<reco::Muon> > LeptonFSRProducer::muons_

private

Definition at line 88 of file LeptonFSRProducer.cc.

Referenced by produce().

const edm::EDGetTokenT<pat::PackedCandidateCollection> LeptonFSRProducer::pfcands_

private

Definition at line 86 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::photonPtCut

private

Definition at line 94 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::ptCutE

private

Definition at line 92 of file LeptonFSRProducer.cc.

Referenced by produce().

const double LeptonFSRProducer::ptCutMu

private

Definition at line 90 of file LeptonFSRProducer.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation