#include <PhysicsTools/NanoAOD/plugins/LeptonJetVarProducer.cc>

Inheritance diagram for LeptonJetVarProducer< T >:

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

	~LeptonJetVarProducer () override

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

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)

virtual	~ProducerBase () noexcept(false)

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

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::vector< ModuleDescription const * > &modules, 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	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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
std::tuple< float, float, float >	calculatePtRatioRel (edm::Ptr< reco::Candidate > lep, edm::Ptr< pat::Jet > jet, const reco::Vertex &vtx) const

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

Private Attributes
edm::EDGetTokenT< edm::View< pat::Jet > >	srcJet_

edm::EDGetTokenT< edm::View< T > >	srcLep_

edm::EDGetTokenT< std::vector< reco::Vertex > >	srcVtx_

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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 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)

Detailed Description

template<typename T>
class LeptonJetVarProducer< T >

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 48 of file LeptonJetVarProducer.cc.

Constructor & Destructor Documentation

template<typename T >

LeptonJetVarProducer< T >::LeptonJetVarProducer ( const edm::ParameterSet & iConfig )

inlineexplicit

Definition at line 50 of file LeptonJetVarProducer.cc.

                                                                :
     srcJet_(consumes<edm::View<pat::Jet>>(iConfig.getParameter<edm::InputTag>("srcJet"))),
     srcLep_(consumes<edm::View<T>>(iConfig.getParameter<edm::InputTag>("srcLep"))),
     srcVtx_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("srcVtx")))
   {
     produces<edm::ValueMap<float>>("ptRatio");
     produces<edm::ValueMap<float>>("ptRel");
     produces<edm::ValueMap<float>>("jetNDauChargedMVASel");
     produces<edm::ValueMap<reco::CandidatePtr>>("jetForLepJetVar");
   }

template<typename T >

LeptonJetVarProducer< T >::~LeptonJetVarProducer ( )

inlineoverride

Definition at line 60 of file LeptonJetVarProducer.cc.

References LeptonJetVarProducer< T >::calculatePtRatioRel(), LeptonJetVarProducer< T >::fillDescriptions(), metsig::jet, LeptonJetVarProducer< T >::produce(), and badGlobalMuonTaggersAOD_cff::vtx.

60 {};

Member Function Documentation

template<typename T >

std::tuple< float, float, float > LeptonJetVarProducer< T >::calculatePtRatioRel	(	edm::Ptr< reco::Candidate >	lep,
		edm::Ptr< pat::Jet >	jet,
		const reco::Vertex &	vtx
	)		const

private

Definition at line 156 of file LeptonJetVarProducer.cc.

References pat::Jet::correctedP4(), edmIntegrityCheck::d, pat::Jet::daughterPtrVector(), deltaR(), MillePedeFileConverter_cfg::e, objects.autophobj::float, pat::Jet::jecFactor(), reco::Candidate::p4(), reco::Vertex::position(), reco::LeafCandidate::pt(), and dttmaxenums::R.

Referenced by LeptonJetVarProducer< T >::produce(), and LeptonJetVarProducer< T >::~LeptonJetVarProducer().

                                                                                                                        {
  
   auto rawp4 = jet->correctedP4("Uncorrected");
   auto lepp4 = lep->p4();
 
   if ((rawp4-lepp4).R()<1e-4) return std::tuple<float,float,float>(1.0,0.0,0.0);
 
   auto l1corrFactor = jet->jecFactor("L1FastJet")/jet->jecFactor("Uncorrected");
 
   auto jetp4 = (rawp4 - lepp4*(1.0/l1corrFactor))*(jet->pt()/rawp4.pt())+lepp4;
   auto ptratio = lepp4.pt()/jetp4.pt();
   auto ptrel = lepp4.Vect().Cross((jetp4-lepp4).Vect().Unit()).R();
 
   unsigned int jndau = 0;
   for(const auto _d : jet->daughterPtrVector()) {
     const auto d = dynamic_cast<const pat::PackedCandidate*>(_d.get());
     if (d->charge()==0) continue;
     if (d->fromPV()<=1) continue;
     if (deltaR(*d,*lep)>0.4) continue;
     if (!(d->hasTrackDetails())) continue;
     auto tk = d->pseudoTrack();
     if(tk.pt()>1 &&
        tk.hitPattern().numberOfValidHits()>=8 &&
        tk.hitPattern().numberOfValidPixelHits()>=2 &&
        tk.normalizedChi2()<5 &&
        fabs(tk.dxy(vtx.position()))<0.2 &&
        fabs(tk.dz(vtx.position()))<17
        ) jndau++;
   }
 
   return std::tuple<float,float,float>(ptratio,ptrel,float(jndau));
 }

template<typename T >

void LeptonJetVarProducer< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 192 of file LeptonJetVarProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by LeptonJetVarProducer< T >::~LeptonJetVarProducer().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("srcJet")->setComment("jet input collection");
   desc.add<edm::InputTag>("srcLep")->setComment("lepton input collection");
   desc.add<edm::InputTag>("srcVtx")->setComment("primary vertex input collection");
   std::string modname;
   if (typeid(T) == typeid(pat::Muon)) modname+="Muon";
   else if (typeid(T) == typeid(pat::Electron)) modname+="Electron";
   modname+="JetVarProducer";
   descriptions.add(modname,desc);
 }

template<typename T >

void LeptonJetVarProducer< T >::produce	(	edm::StreamID	streamID,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overrideprivate

Definition at line 92 of file LeptonJetVarProducer.cc.

References LeptonJetVarProducer< T >::calculatePtRatioRel(), edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), edm::helper::Filler< Map >::insert(), metsig::jet, electrons_cff::jetForLepJetVar, electrons_cff::jetNDauChargedMVASel, matchByCommonSourceCandidatePtr(), eostools::move(), electrons_cff::ptRatio, ptRel(), edm::Event::put(), MetAnalyzer::pv(), electrons_cff::srcJet, LeptonJetVarProducer< T >::srcJet_, electrons_cff::srcLep, LeptonJetVarProducer< T >::srcLep_, electrons_cff::srcVtx, and LeptonJetVarProducer< T >::srcVtx_.

Referenced by LeptonJetVarProducer< T >::~LeptonJetVarProducer().

 {
 
   edm::Handle<edm::View<pat::Jet>> srcJet;
   iEvent.getByToken(srcJet_, srcJet);
   edm::Handle<edm::View<T>> srcLep;
   iEvent.getByToken(srcLep_, srcLep);
   edm::Handle<std::vector<reco::Vertex>> srcVtx;
   iEvent.getByToken(srcVtx_, srcVtx);
 
   unsigned int nJet = srcJet->size();
   unsigned int nLep = srcLep->size();
 
   std::vector<float> ptRatio(nLep,-1);
   std::vector<float> ptRel(nLep,-1);
   std::vector<float> jetNDauChargedMVASel(nLep,0);
   std::vector<reco::CandidatePtr> jetForLepJetVar(nLep,reco::CandidatePtr());
 
   const auto & pv = (*srcVtx)[0];
 
   for (unsigned int il = 0; il<nLep; il++){
     for (unsigned int ij = 0; ij<nJet; ij++){
       auto lep = srcLep->ptrAt(il);
       auto jet = srcJet->ptrAt(ij);
       if(matchByCommonSourceCandidatePtr(*lep,*jet)){
       auto res = calculatePtRatioRel(lep,jet,pv);
       ptRatio[il] = std::get<0>(res);
       ptRel[il] = std::get<1>(res);
       jetNDauChargedMVASel[il] = std::get<2>(res);
       jetForLepJetVar[il] = jet;
       break; // take leading jet with shared source candidates
     }
     }
   }
 
   std::unique_ptr<edm::ValueMap<float>> ptRatioV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerRatio(*ptRatioV);
   fillerRatio.insert(srcLep,ptRatio.begin(),ptRatio.end());
   fillerRatio.fill();
   iEvent.put(std::move(ptRatioV),"ptRatio");
 
   std::unique_ptr<edm::ValueMap<float>> ptRelV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerRel(*ptRelV);
   fillerRel.insert(srcLep,ptRel.begin(),ptRel.end());
   fillerRel.fill();
   iEvent.put(std::move(ptRelV),"ptRel");
 
   std::unique_ptr<edm::ValueMap<float>> jetNDauChargedMVASelV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerNDau(*jetNDauChargedMVASelV);
   fillerNDau.insert(srcLep,jetNDauChargedMVASel.begin(),jetNDauChargedMVASel.end());
   fillerNDau.fill();
   iEvent.put(std::move(jetNDauChargedMVASelV),"jetNDauChargedMVASel");
 
   std::unique_ptr<edm::ValueMap<reco::CandidatePtr>> jetForLepJetVarV(new edm::ValueMap<reco::CandidatePtr>());
   edm::ValueMap<reco::CandidatePtr>::Filler fillerjetForLepJetVar(*jetForLepJetVarV);
   fillerjetForLepJetVar.insert(srcLep,jetForLepJetVar.begin(),jetForLepJetVar.end());
   fillerjetForLepJetVar.fill();
   iEvent.put(std::move(jetForLepJetVarV),"jetForLepJetVar");
 
 
 }