Inheritance diagram for ParticleBasedIsoProducer:

Public Member Functions
void	beginRun (edm::Run const &r, edm::EventSetup const &es) override

	ParticleBasedIsoProducer (const edm::ParameterSet &conf)

void	produce (edm::Event &e, const edm::EventSetup &c) 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

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

Private Attributes
edm::ParameterSet	conf_

std::string	electronCollection_

edm::InputTag	electronProducer_

edm::EDGetTokenT< reco::GsfElectronCollection >	electronProducerT_

edm::InputTag	electronTmpProducer_

edm::EDGetTokenT< reco::GsfElectronCollection >	electronTmpProducerT_

edm::EDGetTokenT< reco::PFCandidateCollection >	pfCandidates_

edm::EDGetTokenT< reco::PFCandidateCollection >	pfEgammaCandidates_

std::string	photonCollection_

edm::InputTag	photonProducer_

edm::EDGetTokenT< reco::PhotonCollection >	photonProducerT_

edm::InputTag	photonTmpProducer_

edm::EDGetTokenT< reco::PhotonCollection >	photonTmpProducerT_

std::unique_ptr< PFBlockBasedIsolation >	thePFBlockBasedIsolation_

edm::EDGetTokenT< edm::ValueMap< reco::GsfElectronRef > >	valMapPFCandToEle_

edm::EDGetTokenT< edm::ValueMap< reco::PhotonRef > >	valMapPFCandToPhoton_

std::string	valueMapElePFCandIso_

std::string	valueMapPFCandEle_

std::string	valueMapPFCandPhoton_

std::string	valueMapPhoPFCandIso_

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 16 of file ParticleBasedIsoProducer.cc.

Constructor & Destructor Documentation

◆ ParticleBasedIsoProducer()

ParticleBasedIsoProducer::ParticleBasedIsoProducer ( const edm::ParameterSet & conf )

Definition at line 56 of file ParticleBasedIsoProducer.cc.

References conf_, electronProducer_, electronProducerT_, electronTmpProducer_, electronTmpProducerT_, edm::ParameterSet::getParameter(), pfCandidates_, pfEgammaCandidates_, photonProducer_, photonProducerT_, photonTmpProducer_, photonTmpProducerT_, AlCaHLTBitMon_QueryRunRegistry::string, valMapPFCandToEle_, valMapPFCandToPhoton_, valueMapElePFCandIso_, valueMapPFCandEle_, valueMapPFCandPhoton_, and valueMapPhoPFCandIso_.

                                                                               : conf_(conf) {
   photonTmpProducer_ = conf_.getParameter<edm::InputTag>("photonTmpProducer");
   photonProducer_ = conf_.getParameter<edm::InputTag>("photonProducer");
   electronProducer_ = conf_.getParameter<edm::InputTag>("electronProducer");
   electronTmpProducer_ = conf_.getParameter<edm::InputTag>("electronTmpProducer");
 
   photonProducerT_ = consumes<reco::PhotonCollection>(photonProducer_);
 
   photonTmpProducerT_ = consumes<reco::PhotonCollection>(photonTmpProducer_);
 
   electronProducerT_ = consumes<reco::GsfElectronCollection>(electronProducer_);
 
   electronTmpProducerT_ = consumes<reco::GsfElectronCollection>(electronTmpProducer_);
 
   pfCandidates_ = consumes<reco::PFCandidateCollection>(conf_.getParameter<edm::InputTag>("pfCandidates"));
 
   pfEgammaCandidates_ = consumes<reco::PFCandidateCollection>(conf_.getParameter<edm::InputTag>("pfEgammaCandidates"));
 
   valueMapPFCandPhoton_ = conf_.getParameter<std::string>("valueMapPhoToEG");
   valueMapPFCandEle_ = conf_.getParameter<std::string>("valueMapEleToEG");
 
   valMapPFCandToPhoton_ = consumes<edm::ValueMap<reco::PhotonRef>>({"gedPhotonsTmp", valueMapPFCandPhoton_});
 
   valMapPFCandToEle_ =
       consumes<edm::ValueMap<reco::GsfElectronRef>>({"gedGsfElectronValueMapsTmp", valueMapPFCandEle_});
 
   valueMapPhoPFCandIso_ = conf_.getParameter<std::string>("valueMapPhoPFblockIso");
   valueMapElePFCandIso_ = conf_.getParameter<std::string>("valueMapElePFblockIso");
 
   produces<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(valueMapPhoPFCandIso_);
   produces<edm::ValueMap<std::vector<reco::PFCandidateRef>>>(valueMapElePFCandIso_);
 }

Member Function Documentation

◆ beginRun()

void ParticleBasedIsoProducer::beginRun	(	edm::Run const &	r,
		edm::EventSetup const &	es
	)

override

Definition at line 89 of file ParticleBasedIsoProducer.cc.

References conf_, edm::ParameterSet::getParameter(), and thePFBlockBasedIsolation_.

                                                                                  {
   thePFBlockBasedIsolation_ = std::make_unique<PFBlockBasedIsolation>();
   edm::ParameterSet pfBlockBasedIsolationSetUp = conf_.getParameter<edm::ParameterSet>("pfBlockBasedIsolationSetUp");
   thePFBlockBasedIsolation_->setup(pfBlockBasedIsolationSetUp);
 }

◆ fillDescriptions()

void ParticleBasedIsoProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 192 of file ParticleBasedIsoProducer.cc.

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

                                                                                           {
   // particleBasedIsolation
   edm::ParameterSetDescription desc;
   desc.add<std::string>("valueMapEleToEG", "");
   desc.add<std::string>("valueMapPhoToEG", "valMapPFEgammaCandToPhoton");
   desc.add<edm::InputTag>("electronTmpProducer", {"gedGsfElectronsTmp"});
   desc.add<edm::InputTag>("pfCandidates", {"particleFlow"});
   desc.add<std::string>("valueMapElePFblockIso", "gedGsfElectrons");
   desc.add<edm::InputTag>("electronProducer", {"gedGsfElectrons"});
   desc.add<edm::InputTag>("photonTmpProducer", {"gedPhotonsTmp"});
   desc.add<edm::InputTag>("pfEgammaCandidates", {"particleFlowEGamma"});
   {
     edm::ParameterSetDescription psd0;
     psd0.add<std::string>("ComponentName", "pfBlockBasedIsolation");
     psd0.add<double>("coneSize", 9999999999.);
     desc.add<edm::ParameterSetDescription>("pfBlockBasedIsolationSetUp", psd0);
   }
   desc.add<edm::InputTag>("photonProducer", {"gedPhotons"});
   desc.add<std::string>("valueMapPhoPFblockIso", "gedPhotons");
   descriptions.add("particleBasedIsolation", desc);
   // or use the following to generate the label from the module's C++ type
   //descriptions.addWithDefaultLabel(desc);
 }

◆ produce()

void ParticleBasedIsoProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)

override

Definition at line 95 of file ParticleBasedIsoProducer.cc.

References electronProducerT_, electronTmpProducerT_, edm::Event::get(), edm::Event::getHandle(), edm::Ref< C, T, F >::isNonnull(), eostools::move(), pfCandidates_, pfEgammaCandidates_, photonProducerT_, photonTmpProducerT_, edm::Event::put(), thePFBlockBasedIsolation_, valMapPFCandToEle_, valMapPFCandToPhoton_, valueMapElePFCandIso_, and valueMapPhoPFCandIso_.

                                                                                  {
   auto photonHandle = theEvent.getHandle(photonProducerT_);
   auto photonTmpHandle = theEvent.getHandle(photonTmpProducerT_);
   auto electronTmpHandle = theEvent.getHandle(electronTmpProducerT_);
   auto electronHandle = theEvent.getHandle(electronProducerT_);
   auto pfEGCandidateHandle = theEvent.getHandle(pfEgammaCandidates_);
   auto pfCandidateHandle = theEvent.getHandle(pfCandidates_);
   auto const& pfEGCandToPhotonMap = theEvent.get(valMapPFCandToPhoton_);
   auto const& pfEGCandToElectronMap = theEvent.get(valMapPFCandToEle_);
 
   std::vector<std::vector<reco::PFCandidateRef>> pfCandIsoPairVecPho;
 
   //  std::cout << " ParticleBasedIsoProducer  photonHandle size " << photonHandle->size() << std::endl;
   for (unsigned int lSC = 0; lSC < photonTmpHandle->size(); lSC++) {
     reco::PhotonRef phoRef(reco::PhotonRef(photonTmpHandle, lSC));
 
     // loop over the unbiased candidates to retrieve the ref to the unbiased candidate corresponding to this photon
     unsigned nObj = pfEGCandidateHandle->size();
     reco::PFCandidateRef pfEGCandRef;
 
     std::vector<reco::PFCandidateRef> pfCandIsoPairPho;
     for (unsigned int lCand = 0; lCand < nObj; lCand++) {
       pfEGCandRef = reco::PFCandidateRef(pfEGCandidateHandle, lCand);
       const reco::PhotonRef& myPho = (pfEGCandToPhotonMap)[pfEGCandRef];
 
       if (myPho.isNonnull()) {
         //std::cout << "ParticleBasedIsoProducer photons PF SC " << pfEGCandRef->superClusterRef()->energy() << " Photon SC " << myPho->superCluster()->energy() << std::endl;
         if (myPho != phoRef)
           continue;
         //  std::cout << " ParticleBasedIsoProducer photons This is my egammaunbiased guy energy " <<  pfEGCandRef->superClusterRef()->energy() << std::endl;
         pfCandIsoPairPho = thePFBlockBasedIsolation_->calculate(myPho->p4(), pfEGCandRef, pfCandidateHandle);
 
         //  for ( std::vector<reco::PFCandidateRef>::const_iterator iPair=pfCandIsoPairPho.begin(); iPair<pfCandIsoPairPho.end(); iPair++) {
         // float dR= deltaR(myPho->eta(),  myPho->phi(), (*iPair)->eta(),  (*iPair)->phi() );
         // std::cout << " ParticleBasedIsoProducer photons  checking the pfCand bool pair " << (*iPair)->particleId() << " dR " << dR << " pt " <<  (*iPair)->pt() << std::endl;
         //  }
       }
     }
 
     pfCandIsoPairVecPho.push_back(pfCandIsoPairPho);
   }
 
   std::vector<std::vector<reco::PFCandidateRef>> pfCandIsoPairVecEle;
   //  std::cout << " ParticleBasedIsoProducer  electronHandle size " << electronHandle->size() << std::endl;
   for (unsigned int lSC = 0; lSC < electronTmpHandle->size(); lSC++) {
     reco::GsfElectronRef eleRef(reco::GsfElectronRef(electronTmpHandle, lSC));
 
     // loop over the unbiased candidates to retrieve the ref to the unbiased candidate corresponding to this electron
     unsigned nObj = pfEGCandidateHandle->size();
     reco::PFCandidateRef pfEGCandRef;
 
     std::vector<reco::PFCandidateRef> pfCandIsoPairEle;
     for (unsigned int lCand = 0; lCand < nObj; lCand++) {
       pfEGCandRef = reco::PFCandidateRef(pfEGCandidateHandle, lCand);
       const reco::GsfElectronRef& myEle = (pfEGCandToElectronMap)[pfEGCandRef];
 
       if (myEle.isNonnull()) {
         //  std::cout << "ParticleBasedIsoProducer Electorns PF SC " << pfEGCandRef->superClusterRef()->energy() << " Electron SC " << myEle->superCluster()->energy() << std::endl;
         if (myEle != eleRef)
           continue;
 
         //math::XYZVector candidateMomentum(myEle->p4().px(),myEle->p4().py(),myEle->p4().pz());
         //math::XYZVector myDir=candidateMomentum.Unit();
         //  std::cout << " ParticleBasedIsoProducer  Electrons This is my egammaunbiased guy energy " <<  pfEGCandRef->superClusterRef()->energy()  << std::endl;
         //  std::cout << " Ele  direction " << myDir << " eta " << myEle->eta() << " phi " << myEle->phi() << std::endl;
         pfCandIsoPairEle = thePFBlockBasedIsolation_->calculate(myEle->p4(), pfEGCandRef, pfCandidateHandle);
         //for ( std::vector<reco::PFCandidateRef>::const_iterator iPair=pfCandIsoPairEle.begin(); iPair<pfCandIsoPairEle.end(); iPair++) {
         // float dR= deltaR(myEle->eta(),  myEle->phi(), (*iPair)->eta(),  (*iPair)->phi() );
         // std::cout << " ParticleBasedIsoProducer Electron  checking the pfCand bool pair " << (*iPair)->particleId() << " dR " << dR << " pt " <<  (*iPair)->pt() << " eta " << (*iPair)->eta() << " phi " << (*iPair)->phi() <<  std::endl;
         //  }
       }
     }
 
     pfCandIsoPairVecEle.push_back(pfCandIsoPairEle);
   }
 
   auto phoToPFCandIsoMap_p = std::make_unique<edm::ValueMap<std::vector<reco::PFCandidateRef>>>();
   edm::ValueMap<std::vector<reco::PFCandidateRef>>::Filler fillerPhotons(*phoToPFCandIsoMap_p);
 
   fillerPhotons.insert(photonHandle, pfCandIsoPairVecPho.begin(), pfCandIsoPairVecPho.end());
   fillerPhotons.fill();
   theEvent.put(std::move(phoToPFCandIsoMap_p), valueMapPhoPFCandIso_);
 
   auto eleToPFCandIsoMap_p = std::make_unique<edm::ValueMap<std::vector<reco::PFCandidateRef>>>();
   edm::ValueMap<std::vector<reco::PFCandidateRef>>::Filler fillerElectrons(*eleToPFCandIsoMap_p);
 
   fillerElectrons.insert(electronHandle, pfCandIsoPairVecEle.begin(), pfCandIsoPairVecEle.end());
   fillerElectrons.fill();
   theEvent.put(std::move(eleToPFCandIsoMap_p), valueMapElePFCandIso_);
 }