#include <EgammaHLTPFPhotonIsolationProducer.h>

Inheritance diagram for EgammaHLTPFPhotonIsolationProducer:

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

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

	~EgammaHLTPFPhotonIsolationProducer ()

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)

	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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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
bool	doRhoCorrection_

double	drMax_

double	drVetoBarrel_

double	drVetoEndcap_

float	effectiveAreaBarrel_

float	effectiveAreaEndcap_

edm::EDGetTokenT < reco::ElectronCollection >	electronProducer_

double	energyBarrel_

double	energyEndcap_

double	etaStripBarrel_

double	etaStripEndcap_

edm::EDGetTokenT < reco::PFCandidateCollection >	pfCandidateProducer_

int	pfToUse_

edm::EDGetTokenT < reco::RecoEcalCandidateCollection >	recoEcalCandidateProducer_

float	rhoMax_

edm::EDGetTokenT< double >	rhoProducer_

float	rhoScale_

bool	useSCRefs_

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 32 of file EgammaHLTPFPhotonIsolationProducer.h.

Constructor & Destructor Documentation

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

explicit

Author: Matteo Sani (UCSD)

$Id:

Definition at line 32 of file EgammaHLTPFPhotonIsolationProducer.cc.

References doRhoCorrection_, drMax_, drVetoBarrel_, drVetoEndcap_, effectiveAreaBarrel_, effectiveAreaEndcap_, electronProducer_, energyBarrel_, energyEndcap_, etaStripBarrel_, etaStripEndcap_, edm::ParameterSet::getParameter(), pfCandidateProducer_, pfToUse_, recoEcalCandidateProducer_, rhoMax_, rhoProducer_, rhoScale_, and useSCRefs_.

                                                                                                     {
 
   pfCandidateProducer_       = consumes<reco::PFCandidateCollection>(config.getParameter<edm::InputTag>("pfCandidatesProducer"));
 
   useSCRefs_ = config.getParameter<bool>("useSCRefs");
 
   drMax_          = config.getParameter<double>("drMax");
   drVetoBarrel_   = config.getParameter<double>("drVetoBarrel");
   drVetoEndcap_   = config.getParameter<double>("drVetoEndcap");
   etaStripBarrel_ = config.getParameter<double>("etaStripBarrel");
   etaStripEndcap_ = config.getParameter<double>("etaStripEndcap");
   energyBarrel_   = config.getParameter<double>("energyBarrel");
   energyEndcap_   = config.getParameter<double>("energyEndcap");
   pfToUse_        = config.getParameter<int>("pfCandidateType");
 
   doRhoCorrection_                = config.getParameter<bool>("doRhoCorrection");
   if (doRhoCorrection_)
     rhoProducer_                    = consumes<double>(config.getParameter<edm::InputTag>("rhoProducer"));
   
   rhoMax_                         = config.getParameter<double>("rhoMax"); 
   rhoScale_                       = config.getParameter<double>("rhoScale"); 
   effectiveAreaBarrel_            = config.getParameter<double>("effectiveAreaBarrel");
   effectiveAreaEndcap_            = config.getParameter<double>("effectiveAreaEndcap");
 
   if(useSCRefs_) {
     produces < reco::RecoEcalCandidateIsolationMap >(); 
     recoEcalCandidateProducer_ = consumes<reco::RecoEcalCandidateCollection>(config.getParameter<edm::InputTag>("recoEcalCandidateProducer"));
   } else {
     produces < reco::ElectronIsolationMap >();
     electronProducer_          = consumes<reco::ElectronCollection>(config.getParameter<edm::InputTag>("electronProducer"));
   }
 }

EgammaHLTPFPhotonIsolationProducer::~EgammaHLTPFPhotonIsolationProducer ( )

inline

Definition at line 35 of file EgammaHLTPFPhotonIsolationProducer.h.

35 {};

Member Function Documentation

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

static

Definition at line 65 of file EgammaHLTPFPhotonIsolationProducer.cc.

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

                                                                                                     {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("electronProducer", edm::InputTag("hltEle27WP80PixelMatchElectronsL1SeededPF"));
   desc.add<edm::InputTag>("recoEcalCandidateProducer", edm::InputTag("hltL1SeededRecoEcalCandidatePF"));
   desc.add<edm::InputTag>("pfCandidatesProducer",  edm::InputTag("hltParticleFlowReg"));
   desc.add<edm::InputTag>("rhoProducer", edm::InputTag("fixedGridRhoFastjetAllCalo"));
   desc.add<bool>("doRhoCorrection", false);
   desc.add<double>("rhoMax", 9.9999999E7); 
   desc.add<double>("rhoScale", 1.0); 
   desc.add<double>("effectiveAreaBarrel", 0.101);
   desc.add<double>("effectiveAreaEndcap", 0.046);
   desc.add<bool>("useSCRefs", false);
   desc.add<double>("drMax", 0.3);
   desc.add<double>("drVetoBarrel", 0.0);
   desc.add<double>("drVetoEndcap", 0.0);
   desc.add<double>("etaStripBarrel", 0.0);
   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
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 88 of file EgammaHLTPFPhotonIsolationProducer.cc.

References b, reco::PFBlockElement::clusterRef(), deltaR(), doRhoCorrection_, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, drMax_, drVetoBarrel_, drVetoEndcap_, reco::PFBlockElement::ECAL, effectiveAreaBarrel_, effectiveAreaEndcap_, electronProducer_, bookConverter::elements, energyBarrel_, energyEndcap_, etaStripBarrel_, etaStripEndcap_, edm::Event::getByToken(), i, edm::AssociationMap< Tag >::insert(), edm::Ref< C, T, F >::isNull(), pfCandidateProducer_, pfToUse_, edm::Handle< T >::product(), edm::Event::put(), recoEcalCandidateProducer_, rho, rhoMax_, rhoProducer_, rhoScale_, reco::PFBlockElement::type(), and useSCRefs_.

                                                                                              {
 
   edm::Handle<double> rhoHandle;
   double rho = 0.0;
   if (doRhoCorrection_) {
     iEvent.getByToken(rhoProducer_, rhoHandle);
     rho = *(rhoHandle.product());
   }
   
   if (rho > rhoMax_)
     rho = rhoMax_;
   
   rho = rho*rhoScale_;
 
   edm::Handle<reco::ElectronCollection> electronHandle;
   edm::Handle<reco::RecoEcalCandidateCollection> recoecalcandHandle;
   edm::Handle<reco::PFCandidateCollection> pfHandle;
 
   iEvent.getByToken(pfCandidateProducer_, pfHandle);
 
   if(useSCRefs_) {
 
     iEvent.getByToken(recoEcalCandidateProducer_,recoecalcandHandle);
     reco::RecoEcalCandidateIsolationMap recoEcalCandMap(recoecalcandHandle);
 
     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_;
       }
       
       float sum = 0;
 
       // Loop over the PFCandidates
       for(unsigned i=0; i<pfHandle->size(); i++) {
     reco::PFCandidateRef pfc(pfHandle, i);
     
     //require that the PFCandidate is a photon
     if (pfc->particleId() == pfToUse_) {
       
       if (fabs(candRef->eta()) < 1.479) {
         if (fabs(pfc->pt()) < energyBarrel_)
           continue;
       } else {
         if (fabs(pfc->energy()) < energyEndcap_)
           continue;
       }
       
       // Shift the RecoEcalCandidate direction vector according to the PF vertex
       math::XYZPoint pfvtx = pfc->vertex();
       math::XYZVector candDirectionWrtVtx(candRef->superCluster()->x() - pfvtx.x(),
                           candRef->superCluster()->y() - pfvtx.y(),
                           candRef->superCluster()->z() - pfvtx.z());
       
       float dEta = fabs(candDirectionWrtVtx.Eta() - pfc->momentum().Eta());
       if(dEta < etaStrip) continue;
       
       float dR = deltaR(candDirectionWrtVtx.Eta(), candDirectionWrtVtx.Phi(), pfc->momentum().Eta(), pfc->momentum().Phi());
       if(dR > drMax_ || dR < dRVeto) continue;
 
       // Exclude PF photons which clusters are part of the candidate 
       bool clusterOverlap = false;
       for(unsigned b=0; b<pfc->elementsInBlocks().size(); b++){
         reco::PFBlockRef blockRef = pfc->elementsInBlocks()[b].first;
         unsigned elementIndex = pfc->elementsInBlocks()[b].second;
         if(blockRef.isNull()) continue;
         const edm::OwnVector< reco::PFBlockElement >& elements = blockRef->elements();
         const reco::PFBlockElement& pfbe(elements[elementIndex]); 
         if( pfbe.type() == reco::PFBlockElement::ECAL ){
           reco::PFClusterRef myPFClusterRef = pfbe.clusterRef();
           if(myPFClusterRef.isNull()) continue;
           for(reco::CaloCluster_iterator it = candRef->superCluster()->clustersBegin(); it != candRef->superCluster()->clustersEnd(); ++it){
         if( myPFClusterRef->seed() == (*it)->seed() ){
           clusterOverlap = true;
           break;
         }
           }
         }
         if(clusterOverlap) break;
       }
       if(clusterOverlap) continue;
       
       sum += pfc->pt();
     }
       }
 
       if (doRhoCorrection_) {
       if (fabs(candRef->eta()) < 1.479) 
     sum = sum - rho*effectiveAreaBarrel_;
       else
     sum = sum - rho*effectiveAreaEndcap_;
       }
       
       recoEcalCandMap.insert(candRef, sum);
     }
     std::auto_ptr<reco::RecoEcalCandidateIsolationMap> mapForEvent(new reco::RecoEcalCandidateIsolationMap(recoEcalCandMap));
     iEvent.put(mapForEvent);
     
   } else {
 
     iEvent.getByToken(electronProducer_,electronHandle);
     reco::ElectronIsolationMap eleMap(electronHandle);
     
     float dRVeto = -1.;
     float etaStrip = -1;
 
     for(unsigned int iEl=0; iEl<electronHandle->size(); iEl++) {
       reco::ElectronRef eleRef(electronHandle, iEl);
 
       if (fabs(eleRef->eta()) < 1.479) {
     dRVeto = drVetoBarrel_;
     etaStrip = etaStripBarrel_;
       } else {
     dRVeto = drVetoEndcap_;
     etaStrip = etaStripEndcap_;
       }
       
       float sum = 0;
 
       // Loop over the PFCandidates
       for(unsigned i=0; i<pfHandle->size(); i++) {
     reco::PFCandidateRef pfc(pfHandle, i);
     
     //require that the PFCandidate is a photon
     if (pfc->particleId() == pfToUse_) {
       
       if (fabs(eleRef->eta()) < 1.479) {
         if (fabs(pfc->pt()) < energyBarrel_)
           continue;
       } else {
         if (fabs(pfc->energy()) < energyEndcap_)
           continue;
       }
 
       float dEta = fabs(eleRef->eta() - pfc->momentum().Eta());
       if(dEta < etaStrip) 
         continue;
 
       float dR = deltaR(eleRef->eta(), eleRef->phi(), pfc->momentum().Eta(), pfc->momentum().Phi());
       if(dR > drMax_ || dR < dRVeto) 
         continue;
 
       // Exclude PF photons which clusters are part of the electron supercluster
       bool clusterOverlap = false;
       for(unsigned b=0; b<pfc->elementsInBlocks().size(); b++){
         reco::PFBlockRef blockRef = pfc->elementsInBlocks()[b].first;
         unsigned elementIndex = pfc->elementsInBlocks()[b].second;
         if(blockRef.isNull()) continue;
         const edm::OwnVector< reco::PFBlockElement >& elements = blockRef->elements();
         const reco::PFBlockElement& pfbe(elements[elementIndex]); 
         if( pfbe.type() == reco::PFBlockElement::ECAL ){
           reco::PFClusterRef myPFClusterRef = pfbe.clusterRef();
           if(myPFClusterRef.isNull()) continue;
           for(reco::CaloCluster_iterator it = eleRef->superCluster()->clustersBegin(); it != eleRef->superCluster()->clustersEnd(); ++it){
         if( myPFClusterRef->seed() == (*it)->seed() ){
           clusterOverlap = true;
           break;
         }
           }
         }
         if(clusterOverlap) break;
       }
       if(clusterOverlap) continue;
 
       sum += pfc->pt();
     }
       }
 
       if (doRhoCorrection_) {
     if (fabs(eleRef->eta()) < 1.479) 
       sum = sum - rho*effectiveAreaBarrel_;
     else
       sum = sum - rho*effectiveAreaEndcap_;
       }
 
       eleMap.insert(eleRef, sum);
     }   
     std::auto_ptr<reco::ElectronIsolationMap> mapForEvent(new reco::ElectronIsolationMap(eleMap));
     iEvent.put(mapForEvent);
   }
 }