#include <EcalPFClusterIsolation.h>

Public Types
typedef std::vector< T1 >	T1Collection

typedef edm::Ref< T1Collection >	T1Ref

Public Member Functions
	EcalPFClusterIsolation (double drMax, double drVetoBarrel, double drVetoEndcap, double etaStripBarrel, double etaStripEndcap, double energyBarrel, double energyEndcap)

double	getSum (T1Ref, edm::Handle< std::vector< reco::PFCluster > >)

	~EcalPFClusterIsolation ()

Private Member Functions
bool	computedRVeto (T1Ref candRef, reco::PFClusterRef pfclu)

template<>
bool	computedRVeto (T1Ref candRef, reco::PFClusterRef pfclu)

Private Attributes
const double	drMax_

double	drVeto2_

const double	drVetoBarrel_

const double	drVetoEndcap_

const double	energyBarrel_

const double	energyEndcap_

const double	etaStripBarrel_

const double	etaStripEndcap_

Detailed Description

template<typename T1>
class EcalPFClusterIsolation< T1 >

Definition at line 32 of file EcalPFClusterIsolation.h.

Member Typedef Documentation

template<typename T1>

typedef std::vector<T1> EcalPFClusterIsolation< T1 >::T1Collection

Definition at line 35 of file EcalPFClusterIsolation.h.

template<typename T1>

typedef edm::Ref<T1Collection> EcalPFClusterIsolation< T1 >::T1Ref

Definition at line 36 of file EcalPFClusterIsolation.h.

Constructor & Destructor Documentation

template<typename T1 >

EcalPFClusterIsolation< T1 >::EcalPFClusterIsolation	(	double	drMax,
		double	drVetoBarrel,
		double	drVetoEndcap,
		double	etaStripBarrel,
		double	etaStripEndcap,
		double	energyBarrel,
		double	energyEndcap
	)

Definition at line 19 of file EcalPFClusterIsolation.cc.

                                                :
   drMax_(drMax),
   drVetoBarrel_(drVetoBarrel),
   drVetoEndcap_(drVetoEndcap),
   etaStripBarrel_(etaStripBarrel),
   etaStripEndcap_(etaStripEndcap),
   energyBarrel_(energyBarrel),
   energyEndcap_(energyEndcap)
 {}

template<typename T1 >

EcalPFClusterIsolation< T1 >::~EcalPFClusterIsolation ( )

Definition at line 36 of file EcalPFClusterIsolation.cc.

37 {}

Member Function Documentation

template<typename T1 >

bool EcalPFClusterIsolation< T1 >::computedRVeto	(	T1Ref	candRef,
		reco::PFClusterRef	pfclu
	)

private

Definition at line 76 of file EcalPFClusterIsolation.cc.

References reco::deltaR2(), EcalPFClusterIsolation< T1 >::drMax_, and edm::Ref< C, T, F >::isNonnull().

Referenced by EcalPFClusterIsolation< T1 >::getSum().

                                                                                     {
 
   float dR2 = deltaR2(candRef->eta(), candRef->phi(), pfclu->eta(), pfclu->phi());
   if(dR2 > (drMax_*drMax_))
     return false;
 
   if (candRef->superCluster().isNonnull()) {
     // Exclude clusters that are part of the candidate
     for (reco::CaloCluster_iterator it = candRef->superCluster()->clustersBegin(); it != candRef->superCluster()->clustersEnd(); ++it) {
       if ((*it)->seed() == pfclu->seed()) {
     return false;
       }
     }
   }
 
   return true;
 }

template<>

bool EcalPFClusterIsolation< reco::RecoChargedCandidate >::computedRVeto	(	T1Ref	candRef,
		reco::PFClusterRef	pfclu
	)

private

Definition at line 95 of file EcalPFClusterIsolation.cc.

References reco::deltaR2(), EcalPFClusterIsolation< T1 >::drMax_, and EcalPFClusterIsolation< T1 >::drVeto2_.

                                                                                                           {
 
   float dR2 = deltaR2(candRef->eta(), candRef->phi(), pfclu->eta(), pfclu->phi());
   if(dR2 > (drMax_*drMax_) || dR2 < drVeto2_)
     return false;
   else
     return true;
 }

template<typename T1>

double EcalPFClusterIsolation< T1 >::getSum	(	T1Ref	,
		edm::Handle< std::vector< reco::PFCluster > >
	)

Definition at line 40 of file EcalPFClusterIsolation.cc.

References EcalPFClusterIsolation< T1 >::computedRVeto(), particleFlow_cfi::dEta, EcalPFClusterIsolation< T1 >::drVeto2_, EcalPFClusterIsolation< T1 >::drVetoBarrel_, EcalPFClusterIsolation< T1 >::drVetoEndcap_, EcalPFClusterIsolation< T1 >::energyBarrel_, EcalPFClusterIsolation< T1 >::energyEndcap_, EcalPFClusterIsolation< T1 >::etaStripBarrel_, EcalPFClusterIsolation< T1 >::etaStripEndcap_, and mps_fire::i.

Referenced by EgammaEcalPFClusterIsolationProducer< T1 >::produce(), and HLTEcalPFClusterIsolationProducer< T1 >::produce().

                                                                                                                {
   
   drVeto2_ = -1.;
   float etaStrip = -1;
  
   if (fabs(candRef->eta()) < 1.479) {
     drVeto2_ = drVetoBarrel_*drVetoBarrel_;
     etaStrip = etaStripBarrel_;
   } else {
     drVeto2_ = drVetoEndcap_*drVetoEndcap_;
     etaStrip = etaStripEndcap_;
   }
    
   float etSum = 0;
   for (size_t i=0; i<clusterHandle->size(); i++) {
     reco::PFClusterRef pfclu(clusterHandle, i);
 
     if (fabs(candRef->eta()) < 1.479) {
       if (fabs(pfclu->pt()) < energyBarrel_)
     continue;
     } else {
       if (fabs(pfclu->energy()) < energyEndcap_)
     continue;
     }
     
     float dEta = fabs(candRef->eta() - pfclu->eta());
     if(dEta < etaStrip) continue;
     if (not computedRVeto(candRef, pfclu)) continue;
 
     etSum += pfclu->pt();
   }
 
   return etSum;
 }