#include <ReducedEGProducer.h>

Inheritance diagram for ReducedEGProducer:

Public Member Functions
virtual void	produce (edm::Event &evt, const edm::EventSetup &es)

	ReducedEGProducer (const edm::ParameterSet &ps)

	~ReducedEGProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::EDGetTokenT < EcalRecHitCollection >	barrelEcalHits_

edm::EDGetTokenT < reco::ConversionCollection >	conversionT_

edm::EDGetTokenT < EcalRecHitCollection >	endcapEcalHits_

std::vector< edm::EDGetTokenT < edm::ValueMap< float > > >	gsfElectronIdTs_

edm::EDGetTokenT < edm::ValueMap< std::vector < reco::PFCandidateRef > > >	gsfElectronPfCandMapT_

edm::EDGetTokenT < reco::GsfElectronCollection >	gsfElectronT_

StringCutObjectSelector < reco::GsfElectron >	keepGsfElectronSel_

StringCutObjectSelector < reco::Photon >	keepPhotonSel_

std::string	outConversions_

std::string	outEBEEClusters_

std::string	outEBRecHits_

std::string	outEERecHits_

std::string	outESClusters_

std::string	outESRecHits_

std::string	outGsfElectronCores_

std::vector< std::string >	outGsfElectronIds_

std::string	outGsfElectronPfCandMap_

std::string	outGsfElectrons_

std::string	outPhotonCores_

std::vector< std::string >	outPhotonIds_

std::string	outPhotonPfCandMap_

std::string	outPhotons_

std::string	outSingleConversions_

std::string	outSuperClusters_

std::vector< edm::EDGetTokenT < edm::ValueMap< bool > > >	photonIdTs_

edm::EDGetTokenT < edm::ValueMap< std::vector < reco::PFCandidateRef > > >	photonPfCandMapT_

edm::EDGetTokenT < reco::PhotonCollection >	photonT_

edm::EDGetTokenT < EcalRecHitCollection >	preshowerEcalHits_

StringCutObjectSelector < reco::GsfElectron >	relinkGsfElectronSel_

StringCutObjectSelector < reco::Photon >	relinkPhotonSel_

edm::EDGetTokenT < reco::ConversionCollection >	singleConversionT_

StringCutObjectSelector < reco::GsfElectron >	slimRelinkGsfElectronSel_

StringCutObjectSelector < reco::Photon >	slimRelinkPhotonSel_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

Select subset of electrons and photons from input collections and produced consistently relinked output collections including associated SuperClusters, CaloClusters and ecal RecHits

Author: J.Bendavid (CERN)

Definition at line 49 of file ReducedEGProducer.h.

Constructor & Destructor Documentation

ReducedEGProducer::ReducedEGProducer ( const edm::ParameterSet & ps )

Definition at line 46 of file ReducedEGProducer.cc.

References barrelEcalHits_, edm::EDConsumerBase::consumes(), conversionT_, endcapEcalHits_, edm::ParameterSet::getParameter(), gsfElectronIdTs_, gsfElectronPfCandMapT_, gsfElectronT_, outConversions_, outEBEEClusters_, outEBRecHits_, outEERecHits_, outESClusters_, outESRecHits_, outGsfElectronCores_, outGsfElectronIds_, outGsfElectronPfCandMap_, outGsfElectrons_, outPhotonCores_, outPhotonIds_, outPhotonPfCandMap_, outPhotons_, outSingleConversions_, outSuperClusters_, photonIdTs_, photonPfCandMapT_, photonT_, preshowerEcalHits_, singleConversionT_, AlCaHLTBitMon_QueryRunRegistry::string, and GlobalPosition_Frontier_DevDB_cff::tag.

                                                                   :
   keepPhotonSel_(config.getParameter<std::string>("keepPhotons")),
   slimRelinkPhotonSel_(config.getParameter<std::string>("slimRelinkPhotons")),
   relinkPhotonSel_(config.getParameter<std::string>("relinkPhotons")),
   keepGsfElectronSel_(config.getParameter<std::string>("keepGsfElectrons")),
   slimRelinkGsfElectronSel_(config.getParameter<std::string>("slimRelinkGsfElectrons")),
   relinkGsfElectronSel_(config.getParameter<std::string>("relinkGsfElectrons"))
 {
 
 
   photonT_ = consumes<reco::PhotonCollection>(config.getParameter<edm::InputTag>("photons"));
   gsfElectronT_ = consumes<reco::GsfElectronCollection>(config.getParameter<edm::InputTag>("gsfElectrons"));
   conversionT_ = consumes<reco::ConversionCollection>(config.getParameter<edm::InputTag>("conversions"));
   singleConversionT_ = consumes<reco::ConversionCollection>(config.getParameter<edm::InputTag>("singleConversions"));
   
   barrelEcalHits_   = 
     consumes<EcalRecHitCollection>(config.getParameter<edm::InputTag>("barrelEcalHits"));
   endcapEcalHits_   = 
     consumes<EcalRecHitCollection>(config.getParameter<edm::InputTag>("endcapEcalHits"));
   preshowerEcalHits_   = 
     consumes<EcalRecHitCollection>(config.getParameter<edm::InputTag>("preshowerEcalHits"));
 
   photonPfCandMapT_ = consumes<edm::ValueMap<std::vector<reco::PFCandidateRef> > >(config.getParameter<edm::InputTag>("photonsPFValMap"));
   gsfElectronPfCandMapT_ = consumes<edm::ValueMap<std::vector<reco::PFCandidateRef> > >(config.getParameter<edm::InputTag>("gsfElectronsPFValMap"));
 
   std::vector<edm::InputTag> photonidinputs(config.getParameter<std::vector<edm::InputTag> >("photonIDSources"));
   for (edm::InputTag &tag : photonidinputs) {
     photonIdTs_.emplace_back(consumes<edm::ValueMap<bool> >(tag));
   }
   
   std::vector<edm::InputTag> gsfelectronidinputs(config.getParameter<std::vector<edm::InputTag> >("gsfElectronIDSources"));
   for (edm::InputTag &tag : gsfelectronidinputs) {
     gsfElectronIdTs_.emplace_back(consumes<edm::ValueMap<float> >(tag));
   }  
   
   //output collections    
   outPhotons_ = "reducedGedPhotons";
   outPhotonCores_ = "reducedGedPhotonCores";
   outGsfElectrons_ = "reducedGedGsfElectrons";
   outGsfElectronCores_ = "reducedGedGsfElectronCores";
   outConversions_ = "reducedConversions";
   outSingleConversions_ = "reducedSingleLegConversions";
   outSuperClusters_ = "reducedSuperClusters";
   outEBEEClusters_ = "reducedEBEEClusters";
   outESClusters_ = "reducedESClusters";
   outEBRecHits_ = "reducedEBRecHits";
   outEERecHits_ = "reducedEERecHits";
   outESRecHits_ = "reducedESRecHits";
   outPhotonPfCandMap_ = "reducedPhotonPfCandMap";
   outGsfElectronPfCandMap_ = "reducedGsfElectronPfCandMap";
   outPhotonIds_ = config.getParameter<std::vector<std::string> >("photonIDOutput");
   outGsfElectronIds_ = config.getParameter<std::vector<std::string> >("gsfElectronIDOutput");
   
   
   produces< reco::PhotonCollection >(outPhotons_);
   produces< reco::PhotonCoreCollection >(outPhotonCores_);
   produces< reco::GsfElectronCollection >(outGsfElectrons_);
   produces< reco::GsfElectronCoreCollection >(outGsfElectronCores_);
   produces< reco::ConversionCollection >(outConversions_);
   produces< reco::ConversionCollection >(outSingleConversions_);
   produces< reco::SuperClusterCollection >(outSuperClusters_);  
   produces< reco::CaloClusterCollection >(outEBEEClusters_);
   produces< reco::CaloClusterCollection >(outESClusters_);
   produces< EcalRecHitCollection >(outEBRecHits_);
   produces< EcalRecHitCollection >(outEERecHits_);
   produces< EcalRecHitCollection >(outESRecHits_);    
   produces< edm::ValueMap<std::vector<reco::PFCandidateRef> > >(outPhotonPfCandMap_);    
   produces< edm::ValueMap<std::vector<reco::PFCandidateRef> > >(outGsfElectronPfCandMap_);   
   for (const std::string &outid : outPhotonIds_) {
     produces< edm::ValueMap<bool> >(outid);   
   }
   for (const std::string &outid : outGsfElectronIds_) {
     produces< edm::ValueMap<float> >(outid);   
   }  
 
 }

ReducedEGProducer::~ReducedEGProducer ( )

Definition at line 123 of file ReducedEGProducer.cc.

124 {

125 }

Member Function Documentation

void ReducedEGProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::EDProducer.

Definition at line 130 of file ReducedEGProducer.cc.

References Reference_intrackfit_cff::barrel, barrelEcalHits_, begin, edm::PtrVectorBase::clear(), edm::RefVector< C, T, F >::clear(), reco::SuperCluster::clearHitsAndFractions(), reco::SuperCluster::clusters(), edm::conversion(), printConversionInfo::conversionHandle, conversions_cfi::conversions, reco::Photon::conversions(), reco::Photon::conversionsOneLeg(), conversionT_, reco::GsfElectron::core(), cond::rpcobgas::detid, DetId::Ecal, EcalBarrel, EcalEndcap, end, endcapEcalHits_, edm::helper::Filler< Map >::fill(), edm::EventSetup::get(), edm::Event::getByToken(), CaloTopology::getSubdetectorTopology(), CaloSubdetectorTopology::getWindow(), gsfElectronCores_cfi::gsfElectronCores, gsfElectronIdTs_, gsfElectronPfCandMapT_, gsfElectrons_cfi::gsfElectrons, gsfElectronT_, edm::helper::Filler< Map >::insert(), keep, keepGsfElectronSel_, keepPhotonSel_, ConversionTools::matchesConversion(), outConversions_, outEBEEClusters_, outEBRecHits_, outEERecHits_, outESClusters_, outESRecHits_, outGsfElectronCores_, outGsfElectronIds_, outGsfElectronPfCandMap_, outGsfElectrons_, outPhotonCores_, outPhotonIds_, outPhotonPfCandMap_, outPhotons_, outSingleConversions_, outSuperClusters_, photonCore_cfi::photonCore, reco::Photon::photonCore(), photonIdTs_, photonPfCandMapT_, interactiveExample::photons, photonT_, reco::SuperCluster::preshowerClusters(), preshowerEcalHits_, edm::Handle< T >::product(), edm::RefVector< C, T, F >::push_back(), edm::PtrVector< T >::push_back(), edm::Event::put(), relinkGsfElectronSel_, relinkPhotonSel_, reco::SuperCluster::seed(), singleConversionT_, edm::PtrVectorBase::size(), edm::RefVector< C, T, F >::size(), reco::CaloCluster::size(), slimRelinkGsfElectronSel_, slimRelinkPhotonSel_, reco::Photon::superCluster(), and reco::GsfElectron::superCluster().

                                                                                       {
 
   //get input collections
   
   
   edm::Handle<reco::PhotonCollection> photonHandle;
   theEvent.getByToken(photonT_, photonHandle);
 
   edm::Handle<reco::GsfElectronCollection> gsfElectronHandle;
   theEvent.getByToken(gsfElectronT_, gsfElectronHandle);  
   
   edm::Handle<reco::ConversionCollection> conversionHandle;
   theEvent.getByToken(conversionT_, conversionHandle);  
 
   edm::Handle<reco::ConversionCollection> singleConversionHandle;
   theEvent.getByToken(singleConversionT_, singleConversionHandle);    
   
   edm::Handle<EcalRecHitCollection> barrelHitHandle;
   theEvent.getByToken(barrelEcalHits_, barrelHitHandle);  
   
   edm::Handle<EcalRecHitCollection> endcapHitHandle;
   theEvent.getByToken(endcapEcalHits_, endcapHitHandle);
 
   edm::Handle<EcalRecHitCollection> preshowerHitHandle;
   theEvent.getByToken(preshowerEcalHits_, preshowerHitHandle);
   
   edm::Handle<edm::ValueMap<std::vector<reco::PFCandidateRef> > > photonPfCandMapHandle;
   theEvent.getByToken(photonPfCandMapT_, photonPfCandMapHandle);  
 
   edm::Handle<edm::ValueMap<std::vector<reco::PFCandidateRef> > > gsfElectronPfCandMapHandle;
   theEvent.getByToken(gsfElectronPfCandMapT_, gsfElectronPfCandMapHandle);    
   
   std::vector<edm::Handle<edm::ValueMap<bool> > > photonIdHandles(photonIdTs_.size());
   for (unsigned int itok=0; itok<photonIdTs_.size(); ++itok) {
     theEvent.getByToken(photonIdTs_[itok],photonIdHandles[itok]);
   }
   
   std::vector<edm::Handle<edm::ValueMap<float> > > gsfElectronIdHandles(gsfElectronIdTs_.size());
   for (unsigned int itok=0; itok<gsfElectronIdTs_.size(); ++itok) {
     theEvent.getByToken(gsfElectronIdTs_[itok],gsfElectronIdHandles[itok]);
   }  
   
   edm::ESHandle<CaloTopology> theCaloTopology;
   theEventSetup.get<CaloTopologyRecord>().get(theCaloTopology);  
   const CaloTopology *caloTopology = & (*theCaloTopology);  
   
   //initialize output collections
   std::auto_ptr<reco::PhotonCollection> photons(new reco::PhotonCollection);
   std::auto_ptr<reco::PhotonCoreCollection> photonCores(new reco::PhotonCoreCollection);
   std::auto_ptr<reco::GsfElectronCollection> gsfElectrons(new reco::GsfElectronCollection);
   std::auto_ptr<reco::GsfElectronCoreCollection> gsfElectronCores(new reco::GsfElectronCoreCollection);
   std::auto_ptr<reco::ConversionCollection> conversions(new reco::ConversionCollection);
   std::auto_ptr<reco::ConversionCollection> singleConversions(new reco::ConversionCollection);
   std::auto_ptr<reco::SuperClusterCollection> superClusters(new reco::SuperClusterCollection);
   std::auto_ptr<reco::CaloClusterCollection> ebeeClusters(new reco::CaloClusterCollection);
   std::auto_ptr<reco::CaloClusterCollection> esClusters(new reco::CaloClusterCollection);
   std::auto_ptr<EcalRecHitCollection> ebRecHits(new EcalRecHitCollection);
   std::auto_ptr<EcalRecHitCollection> eeRecHits(new EcalRecHitCollection);
   std::auto_ptr<EcalRecHitCollection> esRecHits(new EcalRecHitCollection);
   std::auto_ptr<edm::ValueMap<std::vector<reco::PFCandidateRef> > > photonPfCandMap(new edm::ValueMap<std::vector<reco::PFCandidateRef> >);
   std::auto_ptr<edm::ValueMap<std::vector<reco::PFCandidateRef> > > gsfElectronPfCandMap(new edm::ValueMap<std::vector<reco::PFCandidateRef> >);
   
   std::vector<std::auto_ptr<edm::ValueMap<bool> > > photonIds;
   for (unsigned int iid=0; iid<photonIdHandles.size(); ++iid) {
     photonIds.emplace_back(new edm::ValueMap<bool>);
   }
     
   std::vector<std::auto_ptr<edm::ValueMap<float> > > gsfElectronIds;
   for (unsigned int iid=0; iid<gsfElectronIdHandles.size(); ++iid) {
     gsfElectronIds.emplace_back(new edm::ValueMap<float>);
   }
   
   //maps to collection indices of output objects
   std::map<reco::PhotonCoreRef, unsigned int> photonCoreMap;
   std::map<reco::GsfElectronCoreRef, unsigned int> gsfElectronCoreMap;
   std::map<reco::ConversionRef, unsigned int> conversionMap;
   std::map<reco::ConversionRef, unsigned int> singleConversionMap;
   std::map<reco::SuperClusterRef, unsigned int> superClusterMap;
   std::map<reco::CaloClusterPtr, unsigned int> ebeeClusterMap;
   std::map<reco::CaloClusterPtr, unsigned int> esClusterMap;
   std::set<DetId> rechitMap;
   
   std::set<unsigned int> superClusterFullRelinkMap;
   
   //vectors for pfcandidate valuemaps
   std::vector<std::vector<reco::PFCandidateRef>> pfCandIsoPairVecPho;  
   std::vector<std::vector<reco::PFCandidateRef>> pfCandIsoPairVecEle;
   
   //vectors for id valuemaps
   std::vector<std::vector<bool> > photonIdVals(photonIds.size());
   std::vector<std::vector<float> > gsfElectronIdVals(gsfElectronIds.size());
   
   //loop over photons and fill maps
   for (unsigned int ipho=0; ipho<photonHandle->size(); ++ipho) {
     const reco::Photon &photon = (*photonHandle)[ipho];
     
     bool keep = keepPhotonSel_(photon);
     if (!keep) continue;
     
     reco::PhotonRef photonref(photonHandle,ipho);
     
     photons->push_back(photon);
     
     //fill pf candidate value map vector
     pfCandIsoPairVecPho.push_back((*photonPfCandMapHandle)[photonref]);
 
     //fill photon id valuemap vectors
     for (unsigned int iid=0; iid<photonIds.size(); ++iid) {
       photonIdVals[iid].push_back( (*photonIdHandles[iid])[photonref] );
     }    
     
     const reco::PhotonCoreRef &photonCore = photon.photonCore();
     if (!photonCoreMap.count(photonCore)) {
       photonCores->push_back(*photonCore);
       photonCoreMap[photonCore] = photonCores->size() - 1;
     }
     
     bool slimRelink = slimRelinkPhotonSel_(photon);
     //no supercluster relinking unless slimRelink selection is satisfied
     if (!slimRelink) continue;
     
     bool relink = relinkPhotonSel_(photon);
     
     const reco::SuperClusterRef &superCluster = photon.superCluster();
     const auto &mappedsc = superClusterMap.find(superCluster);
     //get index in output collection in order to keep track whether superCluster
     //will be subject to full relinking
     unsigned int mappedscidx = 0;
     if (mappedsc==superClusterMap.end()) {
       superClusters->push_back(*superCluster);
       mappedscidx = superClusters->size() - 1;
       superClusterMap[superCluster] = mappedscidx;
     }
     else {
       mappedscidx = mappedsc->second;
     }
     
     //additionally mark supercluster for full relinking
     if (relink) superClusterFullRelinkMap.insert(mappedscidx);
     
     //conversions only for full relinking
     if (!relink) continue;
     
     const reco::ConversionRefVector &convrefs = photon.conversions();
     for (const reco::ConversionRef &convref : convrefs) {
       if (!conversionMap.count(convref)) {
         conversions->push_back(*convref);
         conversionMap[convref] = conversions->size() - 1;
       }
     }
     
     //explicitly references conversions
     const reco::ConversionRefVector &singleconvrefs = photon.conversionsOneLeg();
     for (const reco::ConversionRef &convref : singleconvrefs) {
       if (!singleConversionMap.count(convref)) {
         singleConversions->push_back(*convref);
         singleConversionMap[convref] = singleConversions->size() - 1;
       }
     }    
     
   }
   
   //loop over electrons and fill maps
   for (unsigned int iele = 0; iele<gsfElectronHandle->size(); ++iele) {
     const reco::GsfElectron &gsfElectron = (*gsfElectronHandle)[iele];
     
     bool keep = keepGsfElectronSel_(gsfElectron);    
     if (!keep) continue;
     
     reco::GsfElectronRef gsfElectronref(gsfElectronHandle,iele);
     
     gsfElectrons->push_back(gsfElectron);
     pfCandIsoPairVecEle.push_back((*gsfElectronPfCandMapHandle)[gsfElectronref]);
     
     //fill electron id valuemap vectors
     for (unsigned int iid=0; iid<gsfElectronIds.size(); ++iid) {
       gsfElectronIdVals[iid].push_back( (*gsfElectronIdHandles[iid])[gsfElectronref] );
     }    
 
     const reco::GsfElectronCoreRef &gsfElectronCore = gsfElectron.core();
     if (!gsfElectronCoreMap.count(gsfElectronCore)) {
       gsfElectronCores->push_back(*gsfElectronCore);
       gsfElectronCoreMap[gsfElectronCore] = gsfElectronCores->size() - 1;
     }    
     
     bool slimRelink = slimRelinkGsfElectronSel_(gsfElectron);
     //no supercluster relinking unless slimRelink selection is satisfied
     if (!slimRelink) continue;
     
     bool relink = relinkGsfElectronSel_(gsfElectron);
     
     const reco::SuperClusterRef &superCluster = gsfElectron.superCluster();
     const auto &mappedsc = superClusterMap.find(superCluster);
     //get index in output collection in order to keep track whether superCluster
     //will be subject to full relinking
     unsigned int mappedscidx = 0;
     if (mappedsc==superClusterMap.end()) {
       superClusters->push_back(*superCluster);
       mappedscidx = superClusters->size() - 1;
       superClusterMap[superCluster] = mappedscidx;
     }
     else {
       mappedscidx = mappedsc->second;
     }
     
     //additionally mark supercluster for full relinking
     if (relink) superClusterFullRelinkMap.insert(mappedscidx);
     
     //conversions only for full relinking
     if (!relink) continue;
     
     //conversions matched by trackrefs
     for (unsigned int iconv = 0; iconv<conversionHandle->size(); ++iconv) {
       const reco::Conversion &conversion = (*conversionHandle)[iconv];
       reco::ConversionRef convref(conversionHandle,iconv);
       
       bool matched = ConversionTools::matchesConversion(gsfElectron,conversion,true,true);
       if (!matched) continue;
       
       if (!conversionMap.count(convref)) {
         conversions->push_back(conversion);
         conversionMap[convref] = conversions->size() - 1;
       }
       
     }
     
     //single leg conversions matched by trackrefs
     for (unsigned int iconv = 0; iconv<singleConversionHandle->size(); ++iconv) {
       const reco::Conversion &conversion = (*singleConversionHandle)[iconv];
       reco::ConversionRef convref(singleConversionHandle,iconv);
       
       bool matched = ConversionTools::matchesConversion(gsfElectron,conversion,true,true);
       if (!matched) continue;
       
       if (!singleConversionMap.count(convref)) {
         singleConversions->push_back(conversion);
         singleConversionMap[convref] = singleConversions->size() - 1;
       }
       
     }    
     
   }
   
   //loop over output SuperClusters and fill maps
   for (unsigned int isc = 0; isc<superClusters->size(); ++isc) {
     reco::SuperCluster &superCluster = (*superClusters)[isc];
     
     //link seed cluster no matter what
     if (!ebeeClusterMap.count(superCluster.seed())) {
       ebeeClusters->push_back(*superCluster.seed());
       ebeeClusterMap[superCluster.seed()] = ebeeClusters->size() - 1;
     }
         
     //only proceed if superCluster is marked for full relinking
     bool fullrelink = superClusterFullRelinkMap.count(isc);
     if (!fullrelink) {
       //zero detid vector which is anyways not useful without stored rechits
       superCluster.clearHitsAndFractions();
       continue; 
     }
     
     for (const reco::CaloClusterPtr &cluster : superCluster.clusters()) {
       if (!ebeeClusterMap.count(cluster)) {
         ebeeClusters->push_back(*cluster);
         ebeeClusterMap[cluster] = ebeeClusters->size() - 1;
       }
       for (std::pair<DetId,float> hitfrac : cluster->hitsAndFractions()) {
         rechitMap.insert(hitfrac.first);
       }
       //make sure to also take all hits in the 5x5 around the max energy xtal
       bool barrel = cluster->hitsAndFractions().front().first.subdetId()==EcalBarrel;
       const EcalRecHitCollection *rhcol = barrel ? barrelHitHandle.product() : endcapHitHandle.product();
       DetId seed = EcalClusterTools::getMaximum(*cluster, rhcol).first;
       
       std::vector<DetId> dets5x5 = (barrel) ? caloTopology->getSubdetectorTopology(DetId::Ecal,EcalBarrel)->getWindow(seed,5,5) : caloTopology->getSubdetectorTopology(DetId::Ecal,EcalEndcap)->getWindow(seed,5,5);
       for (const DetId &detid : dets5x5) {
         rechitMap.insert(detid);
       }
     }
     for (const reco::CaloClusterPtr &cluster : superCluster.preshowerClusters()) {
       if (!esClusterMap.count(cluster)) {
         esClusters->push_back(*cluster);
         esClusterMap[cluster] = esClusters->size() - 1;
       }
       for (std::pair<DetId,float> hitfrac : cluster->hitsAndFractions()) {
         rechitMap.insert(hitfrac.first);
       }      
     }
     
     //conversions matched geometrically
     for (unsigned int iconv = 0; iconv<conversionHandle->size(); ++iconv) {
       const reco::Conversion &conversion = (*conversionHandle)[iconv];
       reco::ConversionRef convref(conversionHandle,iconv);
       
       bool matched = ConversionTools::matchesConversion(superCluster,conversion,0.2);
       if (!matched) continue;
       
       if (!conversionMap.count(convref)) {
         conversions->push_back(conversion);
         conversionMap[convref] = conversions->size() - 1;
       }
       
     }
     
     //single leg conversions matched by trackrefs
     for (unsigned int iconv = 0; iconv<singleConversionHandle->size(); ++iconv) {
       const reco::Conversion &conversion = (*singleConversionHandle)[iconv];
       reco::ConversionRef convref(singleConversionHandle,iconv);
       
       bool matched = ConversionTools::matchesConversion(superCluster,conversion,0.2);
       if (!matched) continue;
       
       if (!singleConversionMap.count(convref)) {
         singleConversions->push_back(conversion);
         singleConversionMap[convref] = singleConversions->size() - 1;
       }
       
     }
     
   }
   
   //now finalize and add to the event collections in "reverse" order
   
   //rechits (fill output collections of rechits to be stored)
   for (const EcalRecHit &rechit : *barrelHitHandle) {
     if (rechitMap.count(rechit.detid())) {
       ebRecHits->push_back(rechit);
     }
   }
   
   for (const EcalRecHit &rechit : *endcapHitHandle) {
     if (rechitMap.count(rechit.detid())) {
       eeRecHits->push_back(rechit);
     }
   }
   
   for (const EcalRecHit &rechit : *preshowerHitHandle) {
     if (rechitMap.count(rechit.detid())) {
       esRecHits->push_back(rechit);
     }
   }
   
   theEvent.put(ebRecHits,outEBRecHits_);
   theEvent.put(eeRecHits,outEERecHits_);
   theEvent.put(esRecHits,outESRecHits_);  
   
   
   //CaloClusters
   //put calocluster output collections in event and get orphan handles to create ptrs
   const edm::OrphanHandle<reco::CaloClusterCollection> &outEBEEClusterHandle = theEvent.put(ebeeClusters,outEBEEClusters_);
   const edm::OrphanHandle<reco::CaloClusterCollection> &outESClusterHandle = theEvent.put(esClusters,outESClusters_);;  
   
   //loop over output superclusters and relink to output caloclusters
   for (reco::SuperCluster &superCluster : *superClusters) {
     //remap seed cluster
     const auto &seedmapped = ebeeClusterMap.find(superCluster.seed());
     if (seedmapped != ebeeClusterMap.end()) {
       //make new ptr
       reco::CaloClusterPtr clusptr(outEBEEClusterHandle,seedmapped->second);
       superCluster.setSeed(clusptr);
     }
     
     //remap all clusters
     reco::CaloClusterPtrVector clusters;
     for (const reco::CaloClusterPtr &cluster : superCluster.clusters()) {
       const auto &clustermapped = ebeeClusterMap.find(cluster);
       if (clustermapped != ebeeClusterMap.end()) {
         //make new ptr
         reco::CaloClusterPtr clusptr(outEBEEClusterHandle,clustermapped->second);
         clusters.push_back(clusptr);
       }
       else {
         //can only relink if all clusters are being relinked, so if one is missing, then skip the relinking completely
         clusters.clear();
         break;
       }
     }
     if (clusters.size()) {
       superCluster.setClusters(clusters);
     }
     
     //remap preshower clusters
     reco::CaloClusterPtrVector esclusters;
     for (const reco::CaloClusterPtr &cluster : superCluster.preshowerClusters()) {
       const auto &clustermapped = esClusterMap.find(cluster);
       if (clustermapped != esClusterMap.end()) {
         //make new ptr
         reco::CaloClusterPtr clusptr(outESClusterHandle,clustermapped->second);
         esclusters.push_back(clusptr);
       }
       else {
         //can only relink if all clusters are being relinked, so if one is missing, then skip the relinking completely
         esclusters.clear();
         break;
       }
     }
     if (esclusters.size()) {
       superCluster.setPreshowerClusters(esclusters);
     }
     
   }
   
   //put superclusters and conversions in the event
   const edm::OrphanHandle<reco::SuperClusterCollection> &outSuperClusterHandle = theEvent.put(superClusters,outSuperClusters_);
   const edm::OrphanHandle<reco::ConversionCollection> &outConversionHandle = theEvent.put(conversions,outConversions_);
   const edm::OrphanHandle<reco::ConversionCollection> &outSingleConversionHandle = theEvent.put(singleConversions,outSingleConversions_);
   
   //loop over photoncores and relink superclusters (and conversions)
   for (reco::PhotonCore &photonCore : *photonCores) {
     const auto &scmapped = superClusterMap.find(photonCore.superCluster());
     if (scmapped != superClusterMap.end()) {
       //make new ref
       reco::SuperClusterRef scref(outSuperClusterHandle,scmapped->second);
       photonCore.setSuperCluster(scref);
     }
     
     //conversions
     const reco::ConversionRefVector &convrefs = photonCore.conversions();
     reco::ConversionRefVector outconvrefs;
     for (const reco::ConversionRef &convref : convrefs) {
       const auto &convmapped = conversionMap.find(convref);
       if (convmapped != conversionMap.end()) {
         //make new ref
         reco::ConversionRef outref(outConversionHandle,convmapped->second);
       }
       else {
         //can only relink if all conversions are being relinked, so if one is missing, then skip the relinking completely
         outconvrefs.clear();
         break;
       }
     }
     if (outconvrefs.size()) {
       photonCore.setConversions(outconvrefs);
     }
     
     //single leg conversions
     const reco::ConversionRefVector &singleconvrefs = photonCore.conversionsOneLeg();
     reco::ConversionRefVector outsingleconvrefs;
     for (const reco::ConversionRef &convref : singleconvrefs) {
       const auto &convmapped = singleConversionMap.find(convref);
       if (convmapped != singleConversionMap.end()) {
         //make new ref
         reco::ConversionRef outref(outSingleConversionHandle,convmapped->second);
       }
       else {
         //can only relink if all conversions are being relinked, so if one is missing, then skip the relinking completely
         outsingleconvrefs.clear();
         break;
       }
     }
     if (outsingleconvrefs.size()) {
       photonCore.setConversionsOneLeg(outsingleconvrefs);
     }    
     
   }
   
   //loop over gsfelectroncores and relink superclusters
   for (reco::GsfElectronCore &gsfElectronCore : *gsfElectronCores) {
     const auto &scmapped = superClusterMap.find(gsfElectronCore.superCluster());
     if (scmapped != superClusterMap.end()) {
       //make new ref
       reco::SuperClusterRef scref(outSuperClusterHandle,scmapped->second);
       gsfElectronCore.setSuperCluster(scref);
     }
   }
   
   //put photon and gsfelectroncores into the event
   const edm::OrphanHandle<reco::PhotonCoreCollection> &outPhotonCoreHandle = theEvent.put(photonCores,outPhotonCores_);
   const edm::OrphanHandle<reco::GsfElectronCoreCollection> &outgsfElectronCoreHandle = theEvent.put(gsfElectronCores,outGsfElectronCores_);
   
   //loop over photons and electrons and relink the cores
   for (reco::Photon &photon : *photons) {
     const auto &coremapped = photonCoreMap.find(photon.photonCore());
     if (coremapped != photonCoreMap.end()) {
       //make new ref
       reco::PhotonCoreRef coreref(outPhotonCoreHandle,coremapped->second);
       photon.setPhotonCore(coreref);
     }
   }
 
   for (reco::GsfElectron &gsfElectron : *gsfElectrons) {
     const auto &coremapped = gsfElectronCoreMap.find(gsfElectron.core());
     if (coremapped != gsfElectronCoreMap.end()) {
       //make new ref
       reco::GsfElectronCoreRef coreref(outgsfElectronCoreHandle,coremapped->second);
       gsfElectron.setCore(coreref);
     }
   }
   
   //(finally) store the output photon and electron collections
   const edm::OrphanHandle<reco::PhotonCollection> &outPhotonHandle = theEvent.put(photons,outPhotons_);  
   const edm::OrphanHandle<reco::GsfElectronCollection> &outGsfElectronHandle = theEvent.put(gsfElectrons,outGsfElectrons_);
   
   //still need to output relinked valuemaps
   
   //photon pfcand isolation valuemap
   edm::ValueMap<std::vector<reco::PFCandidateRef>>::Filler fillerPhotons(*photonPfCandMap);
   fillerPhotons.insert(outPhotonHandle,pfCandIsoPairVecPho.begin(),pfCandIsoPairVecPho.end());
   fillerPhotons.fill();   
   
   //electron pfcand isolation valuemap
   edm::ValueMap<std::vector<reco::PFCandidateRef>>::Filler fillerGsfElectrons(*gsfElectronPfCandMap);
   fillerGsfElectrons.insert(outGsfElectronHandle,pfCandIsoPairVecEle.begin(),pfCandIsoPairVecEle.end());
   fillerGsfElectrons.fill();
   
   theEvent.put(photonPfCandMap,outPhotonPfCandMap_);
   theEvent.put(gsfElectronPfCandMap,outGsfElectronPfCandMap_);
   
   //photon id value maps
   for (unsigned int iid=0; iid<photonIds.size(); ++iid) {
     edm::ValueMap<bool>::Filler fillerPhotonId(*photonIds[iid]);
     fillerPhotonId.insert(outPhotonHandle,photonIdVals[iid].begin(),photonIdVals[iid].end());
     fillerPhotonId.fill();
     theEvent.put(photonIds[iid],outPhotonIds_[iid]);
   }
   
   //electron id value maps
   for (unsigned int iid=0; iid<gsfElectronIds.size(); ++iid) {
     edm::ValueMap<float>::Filler fillerGsfElectronId(*gsfElectronIds[iid]);
     fillerGsfElectronId.insert(outGsfElectronHandle,gsfElectronIdVals[iid].begin(),gsfElectronIdVals[iid].end());
     fillerGsfElectronId.fill();
     theEvent.put(gsfElectronIds[iid],outGsfElectronIds_[iid]);
   }  
   
 }