#include <EgammaHLTRemoveDuplicatedSC.h>

Inheritance diagram for EgammaHLTRemoveDuplicatedSC:

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

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

	~EgammaHLTRemoveDuplicatedSC ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
edm::InputTag	alreadyExistingSC_

std::string	outputCollection_

edm::InputTag	sCInputProducer_

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

typedef WorkerT< EDProducer >	WorkerType

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::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 25 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_, and sCInputProducer_.

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

EgammaHLTRemoveDuplicatedSC::~EgammaHLTRemoveDuplicatedSC ( )

Definition at line 29 of file EgammaHLTRemoveDuplicatedSC.cc.

 {
   ;
 }

Member Function Documentation

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

virtual

Implements edm::EDProducer.

Definition at line 35 of file EgammaHLTRemoveDuplicatedSC.cc.

References alreadyExistingSC_, edm::Event::getByLabel(), outputCollection_, edm::Event::put(), and sCInputProducer_.

 {
   using namespace edm;
 
  
   // Get raw SuperClusters from the event    
   Handle<reco::SuperClusterCollection> UnskimmedSuperCluster;
   evt.getByLabel(sCInputProducer_, UnskimmedSuperCluster);
   Handle<reco::SuperClusterCollection> L1IsoSuperCluster;
   evt.getByLabel(alreadyExistingSC_, L1IsoSuperCluster);
   /*
     edm::LogError("EgammaHLTRemoveDuplicatedSCError") 
       << "Error! can't get the rawSuperClusters " 
       << sCInputProducer_.label() ;
   */
   
 
   // Create a pointer to the existing SuperClusters
     const reco::SuperClusterCollection *UnskimmedL1NonIsoSC = UnskimmedSuperCluster.product();
     const reco::SuperClusterCollection *L1IsoSC = L1IsoSuperCluster.product();
 
     /*
     for(reco::SuperClusterCollection::const_iterator it = L1IsoSC->begin(); it != L1IsoSC->end(); it++){
       std::cout<<"L1 iso  E, eta, phi: "<<it->energy()<<" "<<it->eta()<<" "<<it->phi()<<std::endl;
     }
     std::cout<<std::endl;
     for(reco::SuperClusterCollection::const_iterator it = UnskimmedL1NonIsoSC->begin(); it != UnskimmedL1NonIsoSC->end(); it++){
       std::cout<<"L1 Non iso (not skimmed) E, eta, phi: "<<it->energy()<<" "<<it->eta()<<" "<<it->phi()<<std::endl;
     }
     std::cout<<std::endl;
     */
 
   // Define a collection of corrected SuperClusters to put back into the event
   std::auto_ptr<reco::SuperClusterCollection> corrClusters(new 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(corrClusters, outputCollection_);   
   
 }