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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Constructor & Destructor Documentation

◆ ParticleBasedIsoProducer()

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

Definition at line 56 of file ParticleBasedIsoProducer.cc.

                                                                               : 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_);
 }

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_.

Member Function Documentation

◆ beginRun()

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

override

Definition at line 89 of file ParticleBasedIsoProducer.cc.

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

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

◆ fillDescriptions()

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

static

Definition at line 192 of file ParticleBasedIsoProducer.cc.

                                                                                           {
   // 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);
 }

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

◆ produce()

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

override

Definition at line 95 of file ParticleBasedIsoProducer.cc.

                                                                                  {
   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);
       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);
       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_);
 }