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

Inheritance diagram for IsoValueMapProducer< T >:

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

	~IsoValueMapProducer () 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
void	doMiniIso (edm::Event &) const

template<>
void	doMiniIso (edm::Event &iEvent) const

void	doPFIsoEle (edm::Event &) const

template<>
void	doPFIsoEle (edm::Event &iEvent) const

void	doPFIsoPho (edm::Event &) const

template<>
void	doPFIsoPho (edm::Event &iEvent) const

float	getEtaForEA (const T *) const

template<>
float	getEtaForEA (const pat::Electron *el) const

template<>
float	getEtaForEA (const pat::Photon *ph) const

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

Private Attributes
std::unique_ptr< EffectiveAreas >	ea_miniiso_

std::unique_ptr< EffectiveAreas >	ea_pfiso_

std::unique_ptr< EffectiveAreas >	ea_pfiso_chg_

std::unique_ptr< EffectiveAreas >	ea_pfiso_neu_

std::unique_ptr< EffectiveAreas >	ea_pfiso_pho_

edm::EDGetTokenT< edm::ValueMap< float > >	mapIsoChg_

edm::EDGetTokenT< edm::ValueMap< float > >	mapIsoNeu_

edm::EDGetTokenT< edm::ValueMap< float > >	mapIsoPho_

bool	relative_

edm::EDGetTokenT< double >	rho_miniiso_

edm::EDGetTokenT< double >	rho_pfiso_

edm::EDGetTokenT< edm::View< T > >	src_

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 IsoValueMapProducer< T >

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 45 of file IsoValueMapProducer.cc.

Constructor & Destructor Documentation

template<typename T >

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

inlineexplicit

Definition at line 47 of file IsoValueMapProducer.cc.

References IsoValueMapProducer< T >::ea_miniiso_, IsoValueMapProducer< T >::ea_pfiso_, IsoValueMapProducer< T >::ea_pfiso_chg_, IsoValueMapProducer< T >::ea_pfiso_neu_, IsoValueMapProducer< T >::ea_pfiso_pho_, edm::ParameterSet::getParameter(), IsoValueMapProducer< T >::mapIsoChg_, IsoValueMapProducer< T >::mapIsoNeu_, IsoValueMapProducer< T >::mapIsoPho_, IsoValueMapProducer< T >::rho_miniiso_, and IsoValueMapProducer< T >::rho_pfiso_.

                                                               :
     src_(consumes<edm::View<T>>(iConfig.getParameter<edm::InputTag>("src"))),
     relative_(iConfig.getParameter<bool>("relative"))
   {
     if ((typeid(T) == typeid(pat::Muon)) || (typeid(T) == typeid(pat::Electron)) || typeid(T) == typeid(pat::IsolatedTrack)) {
       produces<edm::ValueMap<float>>("miniIsoChg");
       produces<edm::ValueMap<float>>("miniIsoAll");
       ea_miniiso_.reset(new EffectiveAreas((iConfig.getParameter<edm::FileInPath>("EAFile_MiniIso")).fullPath()));
       rho_miniiso_ = consumes<double>(iConfig.getParameter<edm::InputTag>("rho_MiniIso"));
     }
     if ((typeid(T) == typeid(pat::Electron))) {
       produces<edm::ValueMap<float>>("PFIsoChg");
       produces<edm::ValueMap<float>>("PFIsoAll");
       produces<edm::ValueMap<float>>("PFIsoAll04");
       ea_pfiso_.reset(new EffectiveAreas((iConfig.getParameter<edm::FileInPath>("EAFile_PFIso")).fullPath()));
       rho_pfiso_ = consumes<double>(iConfig.getParameter<edm::InputTag>("rho_PFIso"));
     }
     else if ((typeid(T) == typeid(pat::Photon))) {
       produces<edm::ValueMap<float>>("PFIsoChg");
       produces<edm::ValueMap<float>>("PFIsoAll");
       mapIsoChg_ = consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("mapIsoChg"));
       mapIsoNeu_ = consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("mapIsoNeu"));
       mapIsoPho_ = consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("mapIsoPho"));
       ea_pfiso_chg_.reset(new EffectiveAreas((iConfig.getParameter<edm::FileInPath>("EAFile_PFIso_Chg")).fullPath()));
       ea_pfiso_neu_.reset(new EffectiveAreas((iConfig.getParameter<edm::FileInPath>("EAFile_PFIso_Neu")).fullPath()));
       ea_pfiso_pho_.reset(new EffectiveAreas((iConfig.getParameter<edm::FileInPath>("EAFile_PFIso_Pho")).fullPath()));
       rho_pfiso_ = consumes<double>(iConfig.getParameter<edm::InputTag>("rho_PFIso"));
     }
   }

template<typename T >

IsoValueMapProducer< T >::~IsoValueMapProducer ( )

inlineoverride

Definition at line 76 of file IsoValueMapProducer.cc.

References IsoValueMapProducer< T >::fillDescriptions(), and IsoValueMapProducer< T >::produce().

76 {}

Member Function Documentation

template<typename T >

void IsoValueMapProducer< T >::doMiniIso ( edm::Event & iEvent ) const

private

Definition at line 137 of file IsoValueMapProducer.cc.

References chg, IsoValueMapProducer< T >::ea_miniiso_, edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), IsoValueMapProducer< T >::getEtaForEA(), edm::helper::Filler< Map >::insert(), SiStripPI::max, min(), electrons_cff::miniIsoAll, electrons_cff::miniIsoChg, eostools::move(), NewMuonAssociatorByHits_cfi::obj, funct::pow(), edm::Event::put(), dttmaxenums::R, IsoValueMapProducer< T >::relative_, rho, IsoValueMapProducer< T >::rho_miniiso_, Scenarios_cff::scale, TrackRefitter_38T_cff::src, and IsoValueMapProducer< T >::src_.

Referenced by IsoValueMapProducer< T >::produce().

                                                        {
 
   edm::Handle<edm::View<T>> src;
   iEvent.getByToken(src_, src);
   edm::Handle<double> rho;
   iEvent.getByToken(rho_miniiso_,rho);
 
   unsigned int nInput = src->size();
   
   std::vector<float> miniIsoChg, miniIsoAll;
   miniIsoChg.reserve(nInput);
   miniIsoAll.reserve(nInput);
  
   for (const auto & obj : *src) { 
     auto iso = obj.miniPFIsolation();
     auto chg = iso.chargedHadronIso();
     auto neu = iso.neutralHadronIso();
     auto pho = iso.photonIso();
     auto ea = ea_miniiso_->getEffectiveArea(fabs(getEtaForEA(&obj)));
     float R = 10.0/std::min(std::max(obj.pt(), 50.0),200.0);
     ea *= std::pow(R/0.3,2);
     float scale = relative_ ? 1.0/obj.pt() : 1;
     miniIsoChg.push_back(scale*chg);
     miniIsoAll.push_back(scale*(chg+std::max(0.0,neu+pho-(*rho)*ea)));
   }
   
   std::unique_ptr<edm::ValueMap<float>> miniIsoChgV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerChg(*miniIsoChgV);
   fillerChg.insert(src,miniIsoChg.begin(),miniIsoChg.end());
   fillerChg.fill();
   std::unique_ptr<edm::ValueMap<float>> miniIsoAllV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerAll(*miniIsoAllV);
   fillerAll.insert(src,miniIsoAll.begin(),miniIsoAll.end());
   fillerAll.fill();
 
   iEvent.put(std::move(miniIsoChgV),"miniIsoChg");
   iEvent.put(std::move(miniIsoAllV),"miniIsoAll");
 }

template<>

void IsoValueMapProducer< pat::Photon >::doMiniIso ( edm::Event & iEvent ) const

private

Definition at line 178 of file IsoValueMapProducer.cc.

178 {}

template<typename T >

void IsoValueMapProducer< T >::doPFIsoEle ( edm::Event & iEvent ) const

private

Definition at line 183 of file IsoValueMapProducer.cc.

Referenced by IsoValueMapProducer< T >::produce().

183 {}

template<>

void IsoValueMapProducer< pat::Electron >::doPFIsoEle ( edm::Event & iEvent ) const

private

Definition at line 187 of file IsoValueMapProducer.cc.

References chg, IsoValueMapProducer< T >::ea_pfiso_, edm::helper::Filler< Map >::fill(), edm::Event::getByToken(), IsoValueMapProducer< T >::getEtaForEA(), edm::helper::Filler< Map >::insert(), SiStripPI::max, eostools::move(), NewMuonAssociatorByHits_cfi::obj, electrons_cff::PFIsoAll, electrons_cff::PFIsoAll04, electrons_cff::PFIsoChg, edm::Event::put(), IsoValueMapProducer< T >::relative_, rho, IsoValueMapProducer< T >::rho_pfiso_, Scenarios_cff::scale, TrackRefitter_38T_cff::src, and IsoValueMapProducer< T >::src_.

                                                                   {
 
   edm::Handle<edm::View<pat::Electron>> src;
   iEvent.getByToken(src_, src);
   edm::Handle<double> rho;
   iEvent.getByToken(rho_pfiso_,rho);
   
   unsigned int nInput = src->size();
 
   std::vector<float> PFIsoChg, PFIsoAll, PFIsoAll04;
   PFIsoChg.reserve(nInput);
   PFIsoAll.reserve(nInput);
   PFIsoAll04.reserve(nInput);
   
   for (const auto & obj : *src) { 
     auto iso = obj.pfIsolationVariables();
     auto chg = iso.sumChargedHadronPt;
     auto neu = iso.sumNeutralHadronEt;
     auto pho = iso.sumPhotonEt;
     auto ea = ea_pfiso_->getEffectiveArea(fabs(getEtaForEA(&obj)));
     float scale = relative_ ? 1.0/obj.pt() : 1;
     PFIsoChg.push_back(scale*chg);
     PFIsoAll.push_back(scale*(chg+std::max(0.0,neu+pho-(*rho)*ea)));
     PFIsoAll04.push_back(scale*(obj.chargedHadronIso()+std::max(0.0,obj.neutralHadronIso()+obj.photonIso()-(*rho)*ea*16./9.)));
   }
   
   std::unique_ptr<edm::ValueMap<float>> PFIsoChgV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerChg(*PFIsoChgV);
   fillerChg.insert(src,PFIsoChg.begin(),PFIsoChg.end());
   fillerChg.fill();
   std::unique_ptr<edm::ValueMap<float>> PFIsoAllV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerAll(*PFIsoAllV);
   fillerAll.insert(src,PFIsoAll.begin(),PFIsoAll.end());
   fillerAll.fill();
   std::unique_ptr<edm::ValueMap<float>> PFIsoAll04V(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerAll04(*PFIsoAll04V);
   fillerAll04.insert(src,PFIsoAll04.begin(),PFIsoAll04.end());
   fillerAll04.fill();
 
   iEvent.put(std::move(PFIsoChgV),"PFIsoChg");
   iEvent.put(std::move(PFIsoAllV),"PFIsoAll");
   iEvent.put(std::move(PFIsoAll04V),"PFIsoAll04");
 
 }

template<typename T >

void IsoValueMapProducer< T >::doPFIsoPho ( edm::Event & iEvent ) const

private

Definition at line 234 of file IsoValueMapProducer.cc.

Referenced by IsoValueMapProducer< T >::produce().

234 {}

template<>

void IsoValueMapProducer< pat::Photon >::doPFIsoPho ( edm::Event & iEvent ) const

private

Definition at line 238 of file IsoValueMapProducer.cc.

                                                                  {
 
   edm::Handle<edm::View<pat::Photon>> src;
   iEvent.getByToken(src_, src);
   edm::Handle<double> rho;
   iEvent.getByToken(rho_pfiso_,rho);
   edm::Handle<edm::ValueMap<float> > mapIsoChg;
   iEvent.getByToken(mapIsoChg_, mapIsoChg);
   edm::Handle<edm::ValueMap<float> > mapIsoNeu;
   iEvent.getByToken(mapIsoNeu_, mapIsoNeu);
   edm::Handle<edm::ValueMap<float> > mapIsoPho;
   iEvent.getByToken(mapIsoPho_, mapIsoPho);
   
   unsigned int nInput = src->size();
 
   std::vector<float> PFIsoChg, PFIsoAll;
   PFIsoChg.reserve(nInput);
   PFIsoAll.reserve(nInput);
   
   for (unsigned int i=0; i<nInput; i++){
     auto obj = src->ptrAt(i);
     auto chg = (*mapIsoChg)[obj];
     auto neu = (*mapIsoNeu)[obj];
     auto pho = (*mapIsoPho)[obj];
     auto ea_chg = ea_pfiso_chg_->getEffectiveArea(fabs(getEtaForEA(obj.get())));
     auto ea_neu = ea_pfiso_neu_->getEffectiveArea(fabs(getEtaForEA(obj.get())));
     auto ea_pho = ea_pfiso_pho_->getEffectiveArea(fabs(getEtaForEA(obj.get())));
     float scale = relative_ ? 1.0/obj->pt() : 1;
     PFIsoChg.push_back(scale*std::max(0.0,chg-(*rho)*ea_chg));
     PFIsoAll.push_back(PFIsoChg.back()+scale*(std::max(0.0,neu-(*rho)*ea_neu)+std::max(0.0,pho-(*rho)*ea_pho)));
   }
   
   std::unique_ptr<edm::ValueMap<float>> PFIsoChgV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerChg(*PFIsoChgV);
   fillerChg.insert(src,PFIsoChg.begin(),PFIsoChg.end());
   fillerChg.fill();
   std::unique_ptr<edm::ValueMap<float>> PFIsoAllV(new edm::ValueMap<float>());
   edm::ValueMap<float>::Filler fillerAll(*PFIsoAllV);
   fillerAll.insert(src,PFIsoAll.begin(),PFIsoAll.end());
   fillerAll.fill();
 
   iEvent.put(std::move(PFIsoChgV),"PFIsoChg");
   iEvent.put(std::move(PFIsoAllV),"PFIsoAll");
 
 }

template<typename T >

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

static

Definition at line 289 of file IsoValueMapProducer.cc.

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

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

                                                                                    {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("src")->setComment("input physics object collection");
   desc.add<bool>("relative")->setComment("compute relative isolation instead of absolute one");
   if ((typeid(T) == typeid(pat::Muon)) || (typeid(T) == typeid(pat::Electron)) || typeid(T) == typeid(pat::IsolatedTrack)) {
     desc.add<edm::FileInPath>("EAFile_MiniIso")->setComment("txt file containing effective areas to be used for mini-isolation pileup subtraction");
     desc.add<edm::InputTag>("rho_MiniIso")->setComment("rho to be used for effective-area based mini-isolation pileup subtraction");
   }
   if ((typeid(T) == typeid(pat::Electron))) {
     desc.add<edm::FileInPath>("EAFile_PFIso")->setComment("txt file containing effective areas to be used for PF-isolation pileup subtraction for electrons");
     desc.add<edm::InputTag>("rho_PFIso")->setComment("rho to be used for effective-area based PF-isolation pileup subtraction for electrons");
   }
   if ((typeid(T) == typeid(pat::Photon))) {
     desc.add<edm::InputTag>("mapIsoChg")->setComment("input charged PF isolation calculated in VID for photons");
     desc.add<edm::InputTag>("mapIsoNeu")->setComment("input neutral PF isolation calculated in VID for photons");
     desc.add<edm::InputTag>("mapIsoPho")->setComment("input photon PF isolation calculated in VID for photons");
     desc.add<edm::FileInPath>("EAFile_PFIso_Chg")->setComment("txt file containing effective areas to be used for charged PF-isolation pileup subtraction for photons");
     desc.add<edm::FileInPath>("EAFile_PFIso_Neu")->setComment("txt file containing effective areas to be used for neutral PF-isolation pileup subtraction for photons");
     desc.add<edm::FileInPath>("EAFile_PFIso_Pho")->setComment("txt file containing effective areas to be used for photon PF-isolation pileup subtraction for photons");
     desc.add<edm::InputTag>("rho_PFIso")->setComment("rho to be used for effective-area based PF-isolation pileup subtraction for photons");
   }
   std::string modname;
   if (typeid(T) == typeid(pat::Muon)) modname+="Muon";
   else if (typeid(T) == typeid(pat::Electron)) modname+="Ele";
   else if (typeid(T) == typeid(pat::Photon)) modname+="Pho";
   else if (typeid(T) == typeid(pat::IsolatedTrack)) modname+="IsoTrack";
   modname+="IsoValueMapProducer";
   descriptions.add(modname,desc);
 }

template<typename T >

float IsoValueMapProducer< T >::getEtaForEA ( const T * obj ) const

private

Definition at line 114 of file IsoValueMapProducer.cc.

Referenced by IsoValueMapProducer< T >::doMiniIso(), IsoValueMapProducer< T >::doPFIsoEle(), and IsoValueMapProducer< T >::doPFIsoPho().

                                                                                 {
   return obj->eta();
 }

template<>

float IsoValueMapProducer< pat::Electron >::getEtaForEA ( const pat::Electron * el ) const

private

Definition at line 117 of file IsoValueMapProducer.cc.

References pat::Electron::superCluster().

                                                                                          {
   return el->superCluster()->eta();
 }

template<>

float IsoValueMapProducer< pat::Photon >::getEtaForEA ( const pat::Photon * ph ) const

private

Definition at line 120 of file IsoValueMapProducer.cc.

References pat::Photon::superCluster().

                                                                                      {
   return ph->superCluster()->eta();
 }

template<typename T >

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

overrideprivate

Definition at line 126 of file IsoValueMapProducer.cc.

References IsoValueMapProducer< T >::doMiniIso(), IsoValueMapProducer< T >::doPFIsoEle(), and IsoValueMapProducer< T >::doPFIsoPho().

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

 {
 
   if ((typeid(T) == typeid(pat::Muon)) || (typeid(T) == typeid(pat::Electron)) || typeid(T) == typeid(pat::IsolatedTrack)) { doMiniIso(iEvent); };
   if ((typeid(T) == typeid(pat::Electron))) { doPFIsoEle(iEvent); }
   if ((typeid(T) == typeid(pat::Photon))) { doPFIsoPho(iEvent); }
 
 }