EgammaHLTHcalIsolation Class Reference

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

Public Member Functions
	EgammaHLTHcalIsolation (float eMinHB, float eMinHE, float etMinHB, float etMinHE, float innerCone, float outerCone, int depth)
float	getESum (float candEta, float candPhi, const HBHERecHitCollection *hbhe, const CaloGeometry *geometry) const
float	getESum (float candEta, float candPhi, const HBHERecHitCollection *hbhe, const CaloGeometry *geometry, const HcalSeverityLevelComputer *hcalSevLvlAlgo, const HcalChannelQuality *dbHcalChStatus) const
float	getEtSum (float candEta, float candPhi, const HBHERecHitCollection *hbhe, const CaloGeometry *geometry) const
float	getEtSum (float candEta, float candPhi, const HBHERecHitCollection *hbhe, const CaloGeometry *geometry, const HcalSeverityLevelComputer *hcalSevLvlAlgo, const HcalChannelQuality *dbHcalChStatus) const
std::pair< float, float >	getSum (float candEta, float candPhi, const HBHERecHitCollection *hbhe, const CaloGeometry *geometry, const HcalSeverityLevelComputer *hcalSevLvlAlgo=0, const HcalChannelQuality *dbHcalChStatus=0) const

Static Public Member Functions
static int	getEffectiveDepth (const HcalDetId id)

Private Member Functions
bool	acceptHit_ (const HcalDetId id, const GlobalPoint &pos, const float hitEnergy, const float candEta, const float candPhi) const
bool	passCleaning_ (const CaloRecHit *hit, const HcalSeverityLevelComputer *hcalSevLvlComp, const HcalChannelQuality *hcalChanStatus) const
bool	passDepth_ (const HcalDetId id) const
bool	passMinE_ (float energy, const HcalDetId id) const
bool	passMinEt_ (float et, const HcalDetId id) const

Private Attributes
int	depth_
float	eMinHB_
float	eMinHE_
float	etMinHB_
float	etMinHE_
int	hcalAcceptSeverityLevel_
float	innerCone_
float	outerCone_
bool	useRecoveredHcalHits_

Detailed Description

Description: sum pt hcal hits in cone around egamma candidate

Usage: <usage>

Definition at line 43 of file EgammaHLTHcalIsolation.h.

Constructor & Destructor Documentation

EgammaHLTHcalIsolation::EgammaHLTHcalIsolation ( float  eMinHB, float  eMinHE, float  etMinHB, float  etMinHE, float  innerCone, float  outerCone, int  depth  )

inline

Definition at line 48 of file EgammaHLTHcalIsolation.h.

References getSum(), and photonIsolationHIProducer_cfi::hbhe.

                                                        :
    eMinHB_(eMinHB),eMinHE_(eMinHE),etMinHB_(etMinHB),etMinHE_(etMinHE),
    innerCone_(innerCone),outerCone_(outerCone),depth_(depth),
    useRecoveredHcalHits_(true),hcalAcceptSeverityLevel_(9){}//{std::cout <<"innerCone "<<innerCone<<" outerCone "<<outerCone<<std::endl;}

Member Function Documentation

bool EgammaHLTHcalIsolation::acceptHit_ ( const HcalDetId  id, const GlobalPoint &  pos, const float  hitEnergy, const float  candEta, const float  candPhi  ) const

private

Definition at line 65 of file EgammaHLTHcalIsolation.cc.

References reco::deltaR2(), PV3DBase< T, PVType, FrameType >::eta(), JetChargeProducer_cfi::exp, innerCone_, outerCone_, passDepth_(), passMinE_(), passMinEt_(), PV3DBase< T, PVType, FrameType >::phi(), and funct::sin().

Referenced by getEtSum().

{
  if(passMinE_(hitEnergy,id) && passDepth_(id)){ //doing the energy and depth cuts first will avoid potentially slow eta calc
    float innerConeSq = innerCone_*innerCone_; 
    float outerConeSq = outerCone_*outerCone_;

    float hitEta = pos.eta();
    float hitPhi = pos.phi();
    
    float dR2 = reco::deltaR2(candEta,candPhi,hitEta,hitPhi);
    
    if(dR2>=innerConeSq && dR2<outerConeSq) { //pass inner and outer cone cuts
      float hitEt = hitEnergy*sin(2*atan(exp(-hitEta)));
      if(passMinEt_(hitEt,id)) return true; //and we've passed the last requirement
    }//end dR check
  }//end min energy + depth check
  
  return false;
}

int EgammaHLTHcalIsolation::getEffectiveDepth ( const HcalDetId  id )

static

Definition at line 131 of file EgammaHLTHcalIsolation.cc.

References particleFlowClusterECALTimeSelected_cfi::depth.

Referenced by getEtSum(), and passDepth_().

{
  int iEtaAbs = id.ietaAbs();
  int depth = id.depth();
  if(iEtaAbs<=17 ||
     (iEtaAbs<=29 && depth==1) ||
     (iEtaAbs>=27 && iEtaAbs<=29 && depth==2)){   
    return 1;
  }else return 2;
  
}

float EgammaHLTHcalIsolation::getESum ( float  candEta, float  candPhi, const HBHERecHitCollection *  hbhe, const CaloGeometry *  geometry  ) const

inline

Definition at line 59 of file EgammaHLTHcalIsolation.h.

References getSum().

Referenced by EgammaHLTHcalIsolationProducersRegional::produce().

                                            {return getSum(candEta,candPhi,hbhe,geometry).first;}

float EgammaHLTHcalIsolation::getESum ( float  candEta, float  candPhi, const HBHERecHitCollection *  hbhe, const CaloGeometry *  geometry, const HcalSeverityLevelComputer *  hcalSevLvlAlgo, const HcalChannelQuality *  dbHcalChStatus  ) const

inline

Definition at line 65 of file EgammaHLTHcalIsolation.h.

References getSum().

                                                        {return getSum(candEta,candPhi,hbhe,geometry,
                                           hcalSevLvlAlgo,dbHcalChStatus).first;}

float EgammaHLTHcalIsolation::getEtSum ( float  candEta, float  candPhi, const HBHERecHitCollection *  hbhe, const CaloGeometry *  geometry  ) const

inline

Definition at line 62 of file EgammaHLTHcalIsolation.h.

References getSum().

Referenced by EgammaHLTHcalIsolationProducersRegional::produce().

                                             {return getSum(candEta,candPhi,hbhe,geometry).second;}

float EgammaHLTHcalIsolation::getEtSum ( float  candEta, float  candPhi, const HBHERecHitCollection *  hbhe, const CaloGeometry *  geometry, const HcalSeverityLevelComputer *  hcalSevLvlAlgo, const HcalChannelQuality *  dbHcalChStatus  ) const

inline

Definition at line 71 of file EgammaHLTHcalIsolation.h.

References acceptHit_(), stringResolutionProvider_cfi::et, getEffectiveDepth(), getSum(), passCleaning_(), passDepth_(), passMinE_(), and passMinEt_().

                                                         {return getSum(candEta,candPhi,hbhe,geometry,
                                        hcalSevLvlAlgo,dbHcalChStatus).second;}

std::pair< float, float > EgammaHLTHcalIsolation::getSum	(	float	candEta,
		float	candPhi,
		const HBHERecHitCollection *	hbhe,
		const CaloGeometry *	geometry,
		const HcalSeverityLevelComputer *	hcalSevLvlAlgo = 0,
		const HcalChannelQuality *	dbHcalChStatus = 0
	)		const

Definition at line 24 of file EgammaHLTHcalIsolation.cc.

References particleFlowClusterECALTimeSelected_cfi::depth, photonIsolationHIProducer_cfi::hbhe, funct::sin(), objects.METAnalyzer::sumEt, PV3DBase< T, PVType, FrameType >::theta(), and HcalDigiParam_cfi::zside.

Referenced by EgammaHLTHcalIsolation(), getESum(), and getEtSum().

{
  float sumE=0.;
  float sumEt=0.;
  
  for(HBHERecHitCollection::const_iterator hbheItr = hbhe->begin(); hbheItr != hbhe->end(); ++hbheItr){
    if(passCleaning_(&(*hbheItr),hcalSevLvlAlgo,dbHcalChStatus)){
      //      if(hbheItr->id().ietaAbs()==29) continue;

      HcalDetId id = hbheItr->id();
      if(!(id.ietaAbs()==28 && id.depth()==3)){ //default normal case
    float energy = hbheItr->energy();
    const GlobalPoint& pos = geometry->getPosition(id);
    if(acceptHit_(id,pos,energy,candEta,candPhi)){
      sumE+=energy;
      sumEt+=energy*sin(pos.theta());
    }
      }else{
    //the special case, tower 28 depth 3 is split between tower 28 and 29 when using calo towers so we have to emulate it. To do this we need to divide energy by 2 and then check seperately if 28 and 29 are accepted
    float energy = hbheItr->energy()/2.;
    HcalDetId tower28Id(id.subdet(),28*id.zside(),id.iphi(),2);
    const GlobalPoint& tower28Pos = geometry->getPosition(tower28Id);
    if(acceptHit_(id,tower28Pos,energy,candEta,candPhi)){
      sumE+=energy;
      sumEt+=energy*sin(tower28Pos.theta());
    }
    HcalDetId tower29Id(id.subdet(),29*id.zside(),id.iphi(),2);
    const GlobalPoint& tower29Pos = geometry->getPosition(tower29Id);
    if(acceptHit_(id,tower29Pos,energy,candEta,candPhi)){
      sumE+=energy;
      sumEt+=energy*sin(tower29Pos.theta());
    }
      }//end of the special case for tower 28 depth 3
    }//end cleaning check
  }//end of loop over all rec hits
  return std::make_pair(sumE,sumEt);
  
}

bool EgammaHLTHcalIsolation::passCleaning_ ( const CaloRecHit *  hit, const HcalSeverityLevelComputer *  hcalSevLvlComp, const HcalChannelQuality *  hcalChanStatus  ) const

private

Definition at line 109 of file EgammaHLTHcalIsolation.cc.

References CaloRecHit::detid(), CaloRecHit::flags(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), hcalAcceptSeverityLevel_, HcalSeverityLevelComputer::recoveredRecHit(), interestingDetIdCollectionProducer_cfi::severityLevel, and useRecoveredHcalHits_.

Referenced by getEtSum().

{
  if(hcalSevLvlComp==nullptr || hcalChanStatus==nullptr) return true; //return true if we dont have valid pointers

  const DetId id = hit->detid();
  
  const uint32_t recHitFlag = hit->flags();
  const uint32_t dbStatusFlag = hcalChanStatus->getValues(id)->getValue();
  
  int severityLevel = hcalSevLvlComp->getSeverityLevel(id,recHitFlag,dbStatusFlag);
  bool isRecovered = hcalSevLvlComp->recoveredRecHit(id,recHitFlag);
  
  if(severityLevel == 0) return true;
  else if(isRecovered) return useRecoveredHcalHits_;
  else if(severityLevel <= hcalAcceptSeverityLevel_) return true;
  else return false;
}

bool EgammaHLTHcalIsolation::passDepth_ ( const HcalDetId  id ) const

private

Definition at line 100 of file EgammaHLTHcalIsolation.cc.

References depth_, and getEffectiveDepth().

Referenced by acceptHit_(), and getEtSum().

{
  if(depth_==-1) return true; //I wish we had chosen 0 as all depths but EgammaTowerIsolation chose -1 and 0 as invalid
  else if(getEffectiveDepth(id)==depth_) return true;
  else return false;
    
}

bool EgammaHLTHcalIsolation::passMinE_ ( float  energy, const HcalDetId  id  ) const

private

Definition at line 85 of file EgammaHLTHcalIsolation.cc.

References eMinHB_, eMinHE_, HcalBarrel, and HcalEndcap.

Referenced by acceptHit_(), and getEtSum().

{
  if(id.subdet()==HcalBarrel && energy>=eMinHB_) return true;
  else if(id.subdet()==HcalEndcap && energy>=eMinHE_) return true;
  else return false;
}

bool EgammaHLTHcalIsolation::passMinEt_ ( float  et, const HcalDetId  id  ) const

private

Definition at line 92 of file EgammaHLTHcalIsolation.cc.

References etMinHB_, etMinHE_, HcalBarrel, and HcalEndcap.

Referenced by acceptHit_(), and getEtSum().

{
  if(id.subdet()==HcalBarrel && et>=etMinHB_) return true;
  else if(id.subdet()==HcalEndcap && et>=etMinHE_) return true;
  else return false;

}

Member Data Documentation

int EgammaHLTHcalIsolation::depth_

private

Definition at line 102 of file EgammaHLTHcalIsolation.h.

Referenced by passDepth_().

float EgammaHLTHcalIsolation::eMinHB_

private

Definition at line 96 of file EgammaHLTHcalIsolation.h.

Referenced by passMinE_().

float EgammaHLTHcalIsolation::eMinHE_

private

Definition at line 97 of file EgammaHLTHcalIsolation.h.

Referenced by passMinE_().

float EgammaHLTHcalIsolation::etMinHB_

private

Definition at line 98 of file EgammaHLTHcalIsolation.h.

Referenced by passMinEt_().

float EgammaHLTHcalIsolation::etMinHE_

private

Definition at line 99 of file EgammaHLTHcalIsolation.h.

Referenced by passMinEt_().

int EgammaHLTHcalIsolation::hcalAcceptSeverityLevel_

private

Definition at line 106 of file EgammaHLTHcalIsolation.h.

Referenced by passCleaning_().

float EgammaHLTHcalIsolation::innerCone_

private

Definition at line 100 of file EgammaHLTHcalIsolation.h.

Referenced by acceptHit_().

float EgammaHLTHcalIsolation::outerCone_

private

Definition at line 101 of file EgammaHLTHcalIsolation.h.

Referenced by acceptHit_().

bool EgammaHLTHcalIsolation::useRecoveredHcalHits_

private

Definition at line 105 of file EgammaHLTHcalIsolation.h.

Referenced by passCleaning_().