#include <PFGSFixLinker.h>

Inheritance diagram for PFGSFixLinker:

Public Member Functions
template<typename T >
edm::Handle< T >	getHandle (edm::Event const &_event, edm::EDGetTokenT< T > const &token, std::string const &name)

template<typename T >
void	getToken (edm::EDGetTokenT< T > &token, edm::ParameterSet const &pset, std::string const &label)

	PFGSFixLinker (const edm::ParameterSet &)

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

	~PFGSFixLinker ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef edm::ValueMap < reco::GsfElectronRef >	ElectronRefMap

typedef edm::View < reco::PFCandidate >	PFCandidateView

typedef edm::ValueMap < reco::PhotonRef >	PhotonRefMap

Private Attributes
edm::EDGetTokenT< ElectronRefMap >	electronMapToken_

std::string	electronsMapName_

edm::EDGetTokenT< PFCandidateView >	inputCandidatesToken_

edm::EDGetTokenT < reco::GsfElectronCollection >	inputElectronsToken_

edm::EDGetTokenT < reco::PhotonCollection >	inputPhotonsToken_

edm::EDGetTokenT< PhotonRefMap >	photonMapToken_

std::string	photonsMapName_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Relink GS-fixed e/gamma objects

Author: Y. Iiyama (MIT)

Definition at line 22 of file PFGSFixLinker.h.

Member Typedef Documentation

typedef edm::ValueMap<reco::GsfElectronRef> PFGSFixLinker::ElectronRefMap

private

Definition at line 48 of file PFGSFixLinker.h.

typedef edm::View<reco::PFCandidate> PFGSFixLinker::PFCandidateView

private

Definition at line 47 of file PFGSFixLinker.h.

typedef edm::ValueMap<reco::PhotonRef> PFGSFixLinker::PhotonRefMap

private

Definition at line 49 of file PFGSFixLinker.h.

Constructor & Destructor Documentation

PFGSFixLinker::PFGSFixLinker ( const edm::ParameterSet & _config )

explicit

Definition at line 18 of file PFGSFixLinker.cc.

References electronMapToken_, electronsMapName_, getToken(), inputCandidatesToken_, inputElectronsToken_, inputPhotonsToken_, photonMapToken_, and photonsMapName_.

                                                            :
   electronsMapName_(_config.getParameter<std::string>("ValueMapElectrons")),
   photonsMapName_(_config.getParameter<std::string>("ValueMapPhotons"))
 {
   getToken(inputCandidatesToken_, _config, "PFCandidate");
   getToken(inputElectronsToken_, _config, "GsfElectrons");
   getToken(inputPhotonsToken_, _config, "Photons");
   getToken(electronMapToken_, _config, "GsfElectrons");
   getToken(photonMapToken_, _config, "Photons");
   
   produces<reco::PFCandidateCollection>();
   // new to old
   produces<PFPtrMap>();
 
   // e/g collections to PF
   produces<PFPtrMap>(electronsMapName_);
   produces<PFPtrMap>(photonsMapName_);
 }

PFGSFixLinker::~PFGSFixLinker ( )

Definition at line 37 of file PFGSFixLinker.cc.

38 {

39 }

Member Function Documentation

template<typename T >

edm::Handle<T> PFGSFixLinker::getHandle	(	edm::Event const &	_event,
		edm::EDGetTokenT< T > const &	token,
		std::string const &	name
	)

inline

Definition at line 37 of file PFGSFixLinker.h.

References edm::Event::getByToken(), patZpeak::handle, and mergeVDriftHistosByStation::name.

Referenced by produce().

   {
     edm::Handle<T> handle;
     if (!_event.getByToken(token, handle))
       throw cms::Exception("ProductNotFound") << name;
 
     return handle;
   }

template<typename T >

void PFGSFixLinker::getToken	(	edm::EDGetTokenT< T > &	token,
		edm::ParameterSet const &	pset,
		std::string const &	label
	)

inline

Definition at line 31 of file PFGSFixLinker.h.

References edm::ParameterSet::getParameter(), and diffTwoXMLs::label.

Referenced by PFGSFixLinker().

   {
     token = consumes<T>(pset.getParameter<edm::InputTag>(label));
   }

void PFGSFixLinker::produce	(	edm::Event &	_event,
		const edm::EventSetup &
	)

overridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 42 of file PFGSFixLinker.cc.

References electronMapToken_, electronsMapName_, edm::helper::Filler< Map >::fill(), GainSwitchTools::findNewRef(), getHandle(), customizeTrackingMonitorSeedNumber::idx, inputCandidatesToken_, inputElectronsToken_, inputPhotonsToken_, edm::helper::Filler< Map >::insert(), reco::GsfElectron::P4_COMBINATION, pfElectrons_cff::pfElectrons, electronProducer_cff::pfPhotons, photonMapToken_, photonsMapName_, edm::Event::put(), and reco::Photon::regression2.

 {
   std::auto_ptr<reco::PFCandidateCollection> pOutput(new reco::PFCandidateCollection);
 
   auto& inCandidates(*getHandle(_event, inputCandidatesToken_, "PFCandidate"));
   auto&& inElectronsHandle(getHandle(_event, inputElectronsToken_, "GsfElectrons"));
   auto& electronMap(*getHandle(_event, electronMapToken_, "GsfElectronsMap"));
   auto&& inPhotonsHandle(getHandle(_event, inputPhotonsToken_, "Photons"));
   auto& photonMap(*getHandle(_event, photonMapToken_, "PhotonsMap"));
 
   std::vector<reco::PFCandidatePtr> oldPtrs;
   std::vector<unsigned> pfElectronIndices(inElectronsHandle->size(), -1);
   std::vector<unsigned> pfPhotonIndices(inPhotonsHandle->size(), -1);
 
   unsigned iP(0);
   for (auto& inCand : inCandidates) {
     oldPtrs.emplace_back(inCandidates.ptrAt(iP++));
 
     pOutput->emplace_back(inCand);
     auto& outCand(pOutput->back());
     
     auto&& eRef(inCand.gsfElectronRef());
     if (eRef.isNonnull()) {
       auto&& newERef(GainSwitchTools::findNewRef(eRef, inElectronsHandle, electronMap));
       auto& newE(*newERef);
       auto& newSC(*newE.superCluster());
 
       pfElectronIndices[newERef.key()] = pOutput->size() - 1;
 
       outCand.setGsfElectronRef(newERef);
       outCand.setSuperClusterRef(newE.superCluster());
       outCand.setEcalEnergy(newSC.rawEnergy(), newE.ecalEnergy());
       outCand.setDeltaP(newE.p4Error(reco::GsfElectron::P4_COMBINATION));
       outCand.setP4(newE.p4(reco::GsfElectron::P4_COMBINATION));
     }
 
     auto&& phRef(inCand.photonRef());
     if (phRef.isNonnull()) {
       auto&& newPhRef(GainSwitchTools::findNewRef(phRef, inPhotonsHandle, photonMap));
       auto& newPh(*newPhRef);
       auto& newSC(*newPh.superCluster());
 
       pfPhotonIndices[newPhRef.key()] = pOutput->size() - 1;
     
       outCand.setPhotonRef(newPhRef);
       outCand.setSuperClusterRef(newPh.superCluster());
       outCand.setEcalEnergy(newSC.rawEnergy(), newPh.getCorrectedEnergy(reco::Photon::regression2));
       outCand.setDeltaP(newPh.getCorrectedEnergyError(reco::Photon::regression2));
       outCand.setP4(newPh.p4(reco::Photon::regression2));
     }
   }
 
   auto&& outCandsHandle(_event.put(pOutput));
 
   std::auto_ptr<PFPtrMap> pPFMap(new PFPtrMap);
   PFPtrMap::Filler pfMapFiller(*pPFMap);
   pfMapFiller.insert(outCandsHandle, oldPtrs.begin(), oldPtrs.end());
   pfMapFiller.fill();
   _event.put(pPFMap);
 
   std::vector<reco::PFCandidatePtr> pfElectrons;
   for (unsigned idx : pfElectronIndices) {
     if (idx == unsigned(-1))
       pfElectrons.emplace_back();
     else
       pfElectrons.emplace_back(outCandsHandle, idx);
   }
 
   std::auto_ptr<PFPtrMap> pEleMap(new PFPtrMap);
   PFPtrMap::Filler eleMapFiller(*pEleMap);
   eleMapFiller.insert(inElectronsHandle, pfElectrons.begin(), pfElectrons.end());
   eleMapFiller.fill();
   _event.put(pEleMap, electronsMapName_);
 
   std::vector<reco::PFCandidatePtr> pfPhotons;
   for (unsigned idx : pfPhotonIndices) {
     if (idx == unsigned(-1))
       pfPhotons.emplace_back();
     else
       pfPhotons.emplace_back(outCandsHandle, idx);
   }
 
   std::auto_ptr<PFPtrMap> pPhoMap(new PFPtrMap);
   PFPtrMap::Filler phoMapFiller(*pPhoMap);
   phoMapFiller.insert(inPhotonsHandle, pfPhotons.begin(), pfPhotons.end());
   phoMapFiller.fill();
   _event.put(pPhoMap, photonsMapName_);
 }