#include <EgammaHLTPFPhotonIsolationProducer.h>

Inheritance diagram for EgammaHLTPFPhotonIsolationProducer:

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

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

	~EgammaHLTPFPhotonIsolationProducer ()

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

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

bool	doRhoCorrection_

double	drMax_

double	drVetoBarrel_

double	drVetoEndcap_

double	energyBarrel_

double	energyEndcap_

double	etaStripBarrel_

double	etaStripEndcap_

edm::InputTag	pfCandidates_

int	pfToUse_

edm::InputTag	recoEcalCandidateProducer_

float	rhoMax_

float	rhoScale_

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

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

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

Definition at line 24 of file EgammaHLTPFPhotonIsolationProducer.h.

Constructor & Destructor Documentation

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

explicit

Author: Matteo Sani (UCSD)

$Id:

Definition at line 33 of file EgammaHLTPFPhotonIsolationProducer.cc.

References conf_, drMax_, drVetoBarrel_, drVetoEndcap_, energyBarrel_, energyEndcap_, etaStripBarrel_, etaStripEndcap_, edm::ParameterSet::getParameter(), pfCandidates_, pfToUse_, and recoEcalCandidateProducer_.

                                                                                                     : conf_(config) {
 
   // use configuration file to setup input/output collection names
   pfCandidates_           = conf_.getParameter<edm::InputTag>("pfCandidatesProducer");
   recoEcalCandidateProducer_    = conf_.getParameter<edm::InputTag>("recoEcalCandidateProducer");
 
   drMax_          = conf_.getParameter<double>("drMax");
   drVetoBarrel_   = conf_.getParameter<double>("drVetoBarrel");
   drVetoEndcap_   = conf_.getParameter<double>("drVetoEndcap");
   etaStripBarrel_ = conf_.getParameter<double>("etaStripBarrel");
   etaStripEndcap_ = conf_.getParameter<double>("etaStripEndcap");
   energyBarrel_   = conf_.getParameter<double>("energyBarrel");
   energyEndcap_   = conf_.getParameter<double>("energyEndcap");
   pfToUse_        = conf_.getParameter<int>("pfCandidateType");
   
   produces < reco::RecoEcalCandidateIsolationMap >();
 }

EgammaHLTPFPhotonIsolationProducer::~EgammaHLTPFPhotonIsolationProducer ( )

Definition at line 51 of file EgammaHLTPFPhotonIsolationProducer.cc.

52 {}

Member Function Documentation

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

static

Definition at line 54 of file EgammaHLTPFPhotonIsolationProducer.cc.

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

                                                                                                     {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("pfCandidatesProducer", edm::InputTag("hltParticleFlowReg"));
   desc.add<edm::InputTag>("recoEcalCandidateProducer", edm::InputTag("hltL1SeededRecoEcalCandidatePF"));
   desc.add<double>("drMax", 0.3);
   desc.add<double>("drVetoBarrel", 0.0);
   desc.add<double>("drVetoEndcap", 0.070);
   desc.add<double>("etaStripBarrel", 0.015);
   desc.add<double>("etaStripEndcap", 0.0);
   desc.add<double>("energyBarrel", 0.0);
   desc.add<double>("energyEndcap", 0.0);
   desc.add<int>("pfCandidateType", 4);
   descriptions.add(("hltEgammaHLTPFPhotonIsolationProducer"), desc);
 }

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

virtual

Implements edm::EDProducer.

Definition at line 69 of file EgammaHLTPFPhotonIsolationProducer.cc.

References deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, drMax_, drVetoBarrel_, drVetoEndcap_, reco::LeafCandidate::energy(), energyBarrel_, energyEndcap_, etaStripBarrel_, etaStripEndcap_, edm::Event::getByLabel(), i, edm::AssociationMap< Tag >::insert(), reco::LeafCandidate::momentum(), reco::PFCandidate::particleId(), pfCandidates_, pfToUse_, edm::Handle< T >::product(), reco::LeafCandidate::pt(), edm::Event::put(), recoEcalCandidateProducer_, and reco::PFCandidate::vertex().

                                                                                              {
 
   edm::Handle<reco::PFCandidateCollection> pfHandle;
   iEvent.getByLabel(pfCandidates_, pfHandle);
 
   edm::Handle<reco::RecoEcalCandidateCollection> recoecalcandHandle;
   iEvent.getByLabel(recoEcalCandidateProducer_,recoecalcandHandle);
 
   reco::RecoEcalCandidateIsolationMap isoMap;
 
   float dRVeto = -1.;
   float etaStrip = -1;
   
   for (unsigned int iReco = 0; iReco < recoecalcandHandle->size(); iReco++) {
 
     reco::RecoEcalCandidateRef candRef(recoecalcandHandle, iReco);
     
     if (fabs(candRef->eta()) < 1.479) {
       dRVeto = drVetoBarrel_;
       etaStrip = etaStripBarrel_;
     } else {
       dRVeto = drVetoEndcap_;
       etaStrip = etaStripEndcap_;
     }
       
     const reco::PFCandidateCollection* forIsolation = pfHandle.product();
 
     float sum = 0;
     for(unsigned i=0; i<forIsolation->size(); i++) {
     
       const reco::PFCandidate& pfc = (*forIsolation)[i];
       
       if (pfc.particleId() ==  pfToUse_) {
     
     // FIXME
     // Do not include the PFCandidate associated by SC Ref to the reco::Photon
     //if(pfc.superClusterRef().isNonnull() && localPho->superCluster().isNonnull()) {
     //  if (pfc.superClusterRef() == localPho->superCluster()) 
     //    continue;
     //}
     
     if (fabs(candRef->eta()) < 1.479) {
       if (fabs(pfc.pt()) < energyBarrel_)
         continue;
     } else {
       if (fabs(pfc.energy()) < energyEndcap_)
         continue;
     }
     
     // Shift the photon direction vector according to the PF vertex
     math::XYZPoint pfvtx = pfc.vertex();
     math::XYZVector photon_directionWrtVtx(candRef->superCluster()->x() - pfvtx.x(),
                            candRef->superCluster()->y() - pfvtx.y(),
                            candRef->superCluster()->z() - pfvtx.z());
     
     float dEta = fabs(photon_directionWrtVtx.Eta() - pfc.momentum().Eta());
     float dR = deltaR(photon_directionWrtVtx.Eta(), photon_directionWrtVtx.Phi(), pfc.momentum().Eta(), pfc.momentum().Phi());
     
     if (dEta < etaStrip)
       continue;
     
     if(dR > drMax_ || dR < dRVeto)
       continue;
     
     sum += pfc.pt();
       }
     }
 
     isoMap.insert(candRef, sum);
   }
 
   std::auto_ptr<reco::RecoEcalCandidateIsolationMap> isolMap(new reco::RecoEcalCandidateIsolationMap(isoMap));
   iEvent.put(isolMap);
 }