#include <RecoEgamma/EgammaHLTProducers/interface/EgammaHLTClusterShapeProducer.h>

Inheritance diagram for EgammaHLTClusterShapeProducer:

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

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

	~EgammaHLTClusterShapeProducer ()

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::ParameterSet	conf_

edm::InputTag	ecalRechitEBTag_

edm::InputTag	ecalRechitEETag_

bool	EtaOrIeta_

edm::InputTag	recoEcalCandidateProducer_

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)

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

Author: Roberto Covarelli (CERN)

$Id:

Definition at line 32 of file EgammaHLTClusterShapeProducer.h.

Constructor & Destructor Documentation

EgammaHLTClusterShapeProducer::EgammaHLTClusterShapeProducer ( const edm::ParameterSet & config )

explicit

Definition at line 25 of file EgammaHLTClusterShapeProducer.cc.

References conf_, ecalRechitEBTag_, ecalRechitEETag_, EtaOrIeta_, edm::ParameterSet::getParameter(), and recoEcalCandidateProducer_.

                                                                                           : conf_(config)
 {
  // use configuration file to setup input/output collection names
   recoEcalCandidateProducer_ = conf_.getParameter<edm::InputTag>("recoEcalCandidateProducer");
 
   ecalRechitEBTag_ = conf_.getParameter< edm::InputTag > ("ecalRechitEB");
   ecalRechitEETag_ = conf_.getParameter< edm::InputTag > ("ecalRechitEE");
   EtaOrIeta_ = conf_.getParameter< bool > ("isIeta");
 
   //register your products
   produces < reco::RecoEcalCandidateIsolationMap >();
 }

EgammaHLTClusterShapeProducer::~EgammaHLTClusterShapeProducer ( )

Definition at line 39 of file EgammaHLTClusterShapeProducer.cc.

39 {}

Member Function Documentation

void EgammaHLTClusterShapeProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDProducer.

Definition at line 48 of file EgammaHLTClusterShapeProducer.cc.

References EcalClusterLazyTools::covariances(), ecalRechitEBTag_, ecalRechitEETag_, EtaOrIeta_, edm::Event::getByLabel(), edm::AssociationMap< Tag >::insert(), EcalClusterLazyTools::localCovariances(), edm::Event::put(), recoEcalCandidateProducer_, and mathSSE::sqrt().

 {
   
   // Get the HLT filtered objects
   edm::Handle<reco::RecoEcalCandidateCollection> recoecalcandHandle;
   iEvent.getByLabel(recoEcalCandidateProducer_,recoecalcandHandle);
 
   EcalClusterLazyTools lazyTools( iEvent, iSetup, ecalRechitEBTag_, ecalRechitEETag_ );
   
   reco::RecoEcalCandidateIsolationMap clshMap;
    
   for(reco::RecoEcalCandidateCollection::const_iterator iRecoEcalCand = recoecalcandHandle->begin(); iRecoEcalCand != recoecalcandHandle->end(); iRecoEcalCand++){
     
     reco::RecoEcalCandidateRef recoecalcandref(recoecalcandHandle,iRecoEcalCand-recoecalcandHandle->begin());
     
     std::vector<float> vCov ; 
     double sigmaee;
     if (EtaOrIeta_) {
       vCov = lazyTools.localCovariances( *(recoecalcandref->superCluster()->seed()) );
       sigmaee = sqrt(vCov[0]);
     } else {
       vCov = lazyTools.covariances( *(recoecalcandref->superCluster()->seed()) );
       sigmaee = sqrt(vCov[0]);
       double EtaSC = recoecalcandref->eta();
       if (EtaSC > 1.479) sigmaee = sigmaee - 0.02*(EtaSC - 2.3); 
     }
 
     clshMap.insert(recoecalcandref, sigmaee);
     
   }
 
   std::auto_ptr<reco::RecoEcalCandidateIsolationMap> clushMap(new reco::RecoEcalCandidateIsolationMap(clshMap));
   iEvent.put(clushMap);
 
 }