#include <EgammaHLTRemoveDuplicatedSC.h>

Inheritance diagram for EgammaHLTRemoveDuplicatedSC:

Public Member Functions
	EgammaHLTRemoveDuplicatedSC (const edm::ParameterSet &)

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

	~EgammaHLTRemoveDuplicatedSC ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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 ()

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > 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::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< reco::SuperClusterCollection >	alreadyExistingSC_

std::string	outputCollection_

edm::EDGetTokenT< reco::SuperClusterCollection >	sCInputProducer_

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

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

Definition at line 31 of file EgammaHLTRemoveDuplicatedSC.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file EgammaHLTRemoveDuplicatedSC.cc.

References alreadyExistingSC_, edm::ParameterSet::getParameter(), outputCollection_, sCInputProducer_, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
  
  
   // the input producers
   sCInputProducer_ = consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("L1NonIsoUskimmedSC"));
   alreadyExistingSC_= consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("L1IsoSC"));
 
   // set the producer parameters
   outputCollection_ = ps.getParameter<std::string>("L1NonIsoSkimmedCollection");
   produces<reco::SuperClusterCollection>(outputCollection_);
  
 }

EgammaHLTRemoveDuplicatedSC::~EgammaHLTRemoveDuplicatedSC ( )

Definition at line 28 of file EgammaHLTRemoveDuplicatedSC.cc.

29 {}

Member Function Documentation

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

static

Definition at line 31 of file EgammaHLTRemoveDuplicatedSC.cc.

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

                                                                                              {
 
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>(("L1NonIsoUskimmedSC"), edm::InputTag("hltCorrectedMulti5x5EndcapSuperClustersWithPreshowerL1NonIsolatedTemp"));
   desc.add<edm::InputTag>(("L1IsoSC"), edm::InputTag("hltCorrectedMulti5x5EndcapSuperClustersWithPreshowerL1Isolated"));
   desc.add<std::string>(("L1NonIsoSkimmedCollection"), ""); 
   descriptions.add(("hltEgammaHLTRemoveDuplicatedSC"), desc);  
 }

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

overridevirtual

Definition at line 42 of file EgammaHLTRemoveDuplicatedSC.cc.

References alreadyExistingSC_, edm::Event::getByToken(), eostools::move(), outputCollection_, edm::Handle< T >::product(), edm::Event::put(), and sCInputProducer_.

                                                                            {
  
   using namespace edm;
   
   // Get raw SuperClusters from the event    
   Handle<reco::SuperClusterCollection> UnskimmedSuperCluster;
   evt.getByToken(sCInputProducer_, UnskimmedSuperCluster);
   Handle<reco::SuperClusterCollection> L1IsoSuperCluster;
   evt.getByToken(alreadyExistingSC_, L1IsoSuperCluster);
 
   // Create a pointer to the existing SuperClusters
     const reco::SuperClusterCollection *UnskimmedL1NonIsoSC = UnskimmedSuperCluster.product();
     const reco::SuperClusterCollection *L1IsoSC = L1IsoSuperCluster.product();
 
 
   // Define a collection of corrected SuperClusters to put back into the event
   auto corrClusters = std::make_unique<reco::SuperClusterCollection>();
   
   //  Loop over raw clusters and make corrected ones
   reco::SuperClusterCollection::const_iterator aClus;
   reco::SuperClusterCollection::const_iterator cit;
   for(aClus = UnskimmedL1NonIsoSC->begin(); aClus != UnskimmedL1NonIsoSC->end(); aClus++)
     {
       bool AlreadyThere = false;
       //reco::SuperCluster newClus;
       for(cit = L1IsoSC->begin(); cit != L1IsoSC->end(); cit++){
     if( fabs(aClus->energy()- cit->energy()) < 0.5 && fabs(aClus->eta()- cit->eta()) < 0.0175 ){
       float deltaphi=fabs( aClus->phi() - cit->phi() );
       if(deltaphi>6.283185308) deltaphi -= 6.283185308;
       if(deltaphi>3.141592654) deltaphi = 6.283185308-deltaphi;
 
       if( deltaphi < 0.035 ){AlreadyThere = true; break;}
     }
       }
       // if(AlreadyThere){std::cout<<"AAAA: SC removed: "<<aClus->energy()<<" "<<aClus->eta()<<" "<<aClus->phi()<<std::endl;}
       if(!AlreadyThere){corrClusters->push_back(*aClus);}
     }
 
   // Put collection of corrected SuperClusters into the event
   evt.put(std::move(corrClusters), outputCollection_);
   
 }