EcalPFClusterIsolation< T1 > Class Template Reference

#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.

{}

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(), and edm::Ref< C, T, F >::isNonnull().

                                                                                    {

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

                                                                                                          {

  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 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;
}

Member Data Documentation

template<typename T1>

const double EcalPFClusterIsolation< T1 >::drMax_

private

Definition at line 53 of file EcalPFClusterIsolation.h.

template<typename T1>

double EcalPFClusterIsolation< T1 >::drVeto2_

private

Definition at line 52 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::drVetoBarrel_

private

Definition at line 54 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::drVetoEndcap_

private

Definition at line 55 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::energyBarrel_

private

Definition at line 58 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::energyEndcap_

private

Definition at line 59 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::etaStripBarrel_

private

Definition at line 56 of file EcalPFClusterIsolation.h.

template<typename T1>

const double EcalPFClusterIsolation< T1 >::etaStripEndcap_

private

Definition at line 57 of file EcalPFClusterIsolation.h.