EgammaHLTHcalIsolationDoubleCone Class Reference

#include <RecoEgamma/EgammaHLTAlgos/interface/EgammaHLTHcalIsolationDoubleCone.h>

Public Member Functions
	EgammaHLTHcalIsolationDoubleCone (double egHcalIso_PtMin, double egHcalIso_ConeSize, double egHcalIso_Exclusion)
float	getConeSize () const
	Get isolation cone size. More...
float	getExclusion () const
	Get exclusion region. More...
float	getptMin () const
	Get pt cut for hcal hits. More...
float	isolPtSum (const reco::RecoCandidate *recocandidate, const HBHERecHitCollection *hbhe, const HFRecHitCollection *hf, const CaloGeometry *geometry) const

Private Attributes
const float	conesize
const float	exclusion
const float	ptMin

Detailed Description

Description: sum pt hcal hits in cone around egamma candidate but exlude central cone mostly identical to EgammaHLTHcalIsolation, but with an inner exclusion cone

Usage: <usage>

Definition at line 37 of file EgammaHLTHcalIsolationDoubleCone.h.

Constructor & Destructor Documentation

EgammaHLTHcalIsolationDoubleCone::EgammaHLTHcalIsolationDoubleCone ( double  egHcalIso_PtMin, double  egHcalIso_ConeSize, double  egHcalIso_Exclusion  )

inline

Definition at line 42 of file EgammaHLTHcalIsolationDoubleCone.h.

References photonIsolationHIProducer_cfi::hbhe, photonIsolationHIProducer_cfi::hf, and isolPtSum().

                                                                                                                  :
    ptMin(egHcalIso_PtMin),conesize(egHcalIso_ConeSize),exclusion(egHcalIso_Exclusion){
      /* 
     std::cout << "EgammaHLTHcalIsolation instance:"
     << " ptMin=" << ptMin << "|" << ptMinG
     << " conesize="<< conesize << "|" << conesizeG
     << std::endl;
      */
    }

Member Function Documentation

float EgammaHLTHcalIsolationDoubleCone::getConeSize ( ) const

inline

Get isolation cone size.

Definition at line 59 of file EgammaHLTHcalIsolationDoubleCone.h.

References conesize.

{ return conesize; }

float EgammaHLTHcalIsolationDoubleCone::getExclusion ( ) const

inline

Get exclusion region.

Definition at line 61 of file EgammaHLTHcalIsolationDoubleCone.h.

References exclusion.

{ return exclusion; }

float EgammaHLTHcalIsolationDoubleCone::getptMin ( ) const

inline

Get pt cut for hcal hits.

Definition at line 57 of file EgammaHLTHcalIsolationDoubleCone.h.

References ptMin.

{ return ptMin; }

float EgammaHLTHcalIsolationDoubleCone::isolPtSum	(	const reco::RecoCandidate *	recocandidate,
		const HBHERecHitCollection *	hbhe,
		const HFRecHitCollection *	hf,
		const CaloGeometry *	geometry
	)		const

Definition at line 24 of file EgammaHLTHcalIsolationDoubleCone.cc.

References edm::SortedCollection< T, SORT >::begin(), conesize, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, exclusion, JetChargeProducer_cfi::exp, CaloGeometry::getPosition(), phi, PI, ptMin, funct::sin(), reco::RecoCandidate::superCluster(), and TWOPI.

Referenced by EgammaHLTHcalIsolationDoubleCone(), and EgammaHLTHcalIsolationDoubleConeProducers::produce().

                                                                                                                                                                                            {

  float hcalIsol=0.;

  float candSCphi = recocandidate->superCluster()->phi();
  float candSCeta = recocandidate->superCluster()->eta();
  if(candSCphi<0) candSCphi+=TWOPI;
  float conesizeSquared=conesize*conesize;
  float exclusionSquared= exclusion*exclusion;

  for(HBHERecHitCollection::const_iterator hbheItr = hbhe->begin(); hbheItr != hbhe->end(); ++hbheItr){
    double HcalHit_eta=geometry->getPosition(hbheItr->id()).eta(); //Attention getpos
    if(fabs(HcalHit_eta-candSCeta)<conesize) {
      float HcalHit_pth=hbheItr->energy()*sin(2*atan(exp(-HcalHit_eta)));
      if(HcalHit_pth>ptMin) {
    double HcalHit_phi=geometry->getPosition(hbheItr->id()).phi();
    float deltaeta=fabs(HcalHit_eta-candSCeta);
    if(HcalHit_phi<0) HcalHit_phi+=TWOPI;
    float deltaphi=fabs(HcalHit_phi-candSCphi);
    if(deltaphi>TWOPI) deltaphi-=TWOPI;
    if(deltaphi>PI) deltaphi=TWOPI-deltaphi;
    float newDelta= (deltaphi*deltaphi+ deltaeta*deltaeta);
    if(newDelta<conesizeSquared && newDelta>exclusionSquared ) hcalIsol+=HcalHit_pth;
      }
    }      
  }

  for(HFRecHitCollection::const_iterator hfItr = hf->begin(); hfItr != hf->end(); ++hfItr){
    double HcalHit_eta=geometry->getPosition(hfItr->id()).eta(); //Attention getpos
    if(fabs(HcalHit_eta-candSCeta)<conesize) {
      float HcalHit_pth=hfItr->energy()*sin(2*atan(exp(-HcalHit_eta)));
      if(HcalHit_pth>ptMin) {
    double HcalHit_phi=geometry->getPosition(hfItr->id()).phi();
    float deltaeta=fabs(HcalHit_eta-candSCeta);
    float deltaphi;
    if(HcalHit_phi<0) HcalHit_phi+=TWOPI;
    if(candSCphi<0) candSCphi+=TWOPI;
    deltaphi=fabs(HcalHit_phi-candSCphi);
    if(deltaphi>TWOPI) deltaphi-=TWOPI;
    if(deltaphi>PI) deltaphi=TWOPI-deltaphi;
    float newDelta= (deltaphi*deltaphi+ deltaeta*deltaeta);
    if(newDelta<conesizeSquared && newDelta>exclusionSquared ) hcalIsol+=HcalHit_pth;
      }
    }      
  }


  return hcalIsol;
  
}

Member Data Documentation

const float EgammaHLTHcalIsolationDoubleCone::conesize

private

Definition at line 68 of file EgammaHLTHcalIsolationDoubleCone.h.

Referenced by getConeSize(), and isolPtSum().

const float EgammaHLTHcalIsolationDoubleCone::exclusion

private

Definition at line 69 of file EgammaHLTHcalIsolationDoubleCone.h.

Referenced by getExclusion(), and isolPtSum().

const float EgammaHLTHcalIsolationDoubleCone::ptMin

private

Definition at line 67 of file EgammaHLTHcalIsolationDoubleCone.h.

Referenced by getptMin(), and isolPtSum().