---

EgammaEcalIsolation Class Reference

#include <RecoEgamma/EgammaIsolationAlgos/interface/EgammaEcalIsolation.h>

List of all members.

Public Member Functions
	EgammaEcalIsolation (double extRadius, double etLow, const reco::BasicClusterCollection *, const reco::SuperClusterCollection *)
double	getEcalEtSum (const reco::Candidate *)
	~EgammaEcalIsolation ()
Private Attributes
const reco::BasicClusterCollection *	basicClusterCollection_
double	conesize
double	etMin
const reco::SuperClusterCollection *	superClusterCollection_

---

Detailed Description

Definition at line 26 of file EgammaEcalIsolation.h.

---

Constructor & Destructor Documentation

EgammaEcalIsolation::EgammaEcalIsolation	(	double	extRadius,
		double	etLow,
		const reco::BasicClusterCollection *	basicClusterCollection,
		const reco::SuperClusterCollection *	superClusterCollection
	)

Definition at line 23 of file EgammaEcalIsolation.cc.

                                                                                                   :
  etMin(etLow),
  conesize(extRadius),
  basicClusterCollection_(basicClusterCollection),
  superClusterCollection_(superClusterCollection)
{
}

EgammaEcalIsolation::~EgammaEcalIsolation

(

)

Definition at line 34 of file EgammaEcalIsolation.cc.

{}

---

Member Function Documentation

double EgammaEcalIsolation::getEcalEtSum

(

const reco::Candidate *

candidate

)

Definition at line 36 of file EgammaEcalIsolation.cc.

References reco::BasicCluster::algo(), basicClusterCollection_, reco::BasicCluster::chi2(), reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), conesize, reco::CaloCluster::energy(), reco::CaloCluster::eta(), etMin, funct::exp(), reco::Candidate::get(), edm::Ref< C, T, F >::get(), reco::CaloCluster::position(), reco::SuperCluster::seed(), funct::sin(), and superClusterCollection_.

Referenced by PhotonIDAlgo::calculateBasicClusterIso(), and EgammaEcalIsolationProducer::produce().

                                                                      {

  
  double ecalIsol=0.;
  reco::SuperClusterRef sc = candidate->get<reco::SuperClusterRef>();
  math::XYZVector position(sc.get()->position().x(),
                   sc.get()->position().y(),
                   sc.get()->position().z());
  
  // match the photon hybrid supercluster with those with Algo==0 (island)
  double delta1=1000.;
  double deltacur=1000.;
  const reco::SuperCluster *matchedsupercluster=0;
  bool MATCHEDSC = false;
  
  for(reco::SuperClusterCollection::const_iterator scItr = superClusterCollection_->begin(); scItr != superClusterCollection_->end(); ++scItr){
    
    const reco::SuperCluster *supercluster = &(*scItr);
 
    math::XYZVector currentPosition(supercluster->position().x(),
                     supercluster->position().y(),
                     supercluster->position().z());
 
    
    if(supercluster->seed()->algo() == 0){
      deltacur = DeltaR(currentPosition, position); 
      
      if (deltacur < delta1) {
        delta1=deltacur;
        matchedsupercluster = supercluster;
        MATCHEDSC = true;
      }
    }
  }

  const reco::BasicCluster *cluster= 0;
  
  //loop over basic clusters
  for(reco::BasicClusterCollection::const_iterator cItr = basicClusterCollection_->begin(); cItr != basicClusterCollection_->end(); ++cItr){
 
    cluster = &(*cItr);
    double ebc_bcchi2 = cluster->chi2();
    int   ebc_bcalgo = cluster->algo();
    double ebc_bce    = cluster->energy();
    double ebc_bceta  = cluster->eta();
    double ebc_bcet   = ebc_bce*sin(2*atan(exp(ebc_bceta)));
    double newDelta = 0.;


    if (ebc_bcet > etMin && ebc_bcalgo == 0) {
      if (ebc_bcchi2 < 30.) {
        
        if(MATCHEDSC){
          bool inSuperCluster = false;

          reco::basicCluster_iterator theEclust = matchedsupercluster->clustersBegin();
          // loop over the basic clusters of the matched supercluster
          for(;theEclust != matchedsupercluster->clustersEnd();
              theEclust++) {
            if ((**theEclust) ==  (*cluster) ) inSuperCluster = true;
          }
          if (!inSuperCluster) {

            math::XYZVector basicClusterPosition(cluster->position().x(),
                                                  cluster->position().y(),
                                                  cluster->position().z());
            newDelta=DeltaR(basicClusterPosition,position);
            if(newDelta < conesize) {
              ecalIsol+=ebc_bcet;
            }
          }
        }
      } // matches ebc_bcchi2
    } // matches ebc_bcet && ebc_bcalgo

  }
  
  //  std::cout << "Will return ecalIsol = " << ecalIsol << std::endl; 
  return ecalIsol;
  
}

---

Member Data Documentation

const reco::BasicClusterCollection* EgammaEcalIsolation::basicClusterCollection_ [private]

Definition at line 49 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::conesize [private]

Definition at line 47 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::etMin [private]

Definition at line 46 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

const reco::SuperClusterCollection* EgammaEcalIsolation::superClusterCollection_ [private]

Definition at line 50 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

---

The documentation for this class was generated from the following files:

• RecoEgamma/EgammaIsolationAlgos/interface/EgammaEcalIsolation.h
• RecoEgamma/EgammaIsolationAlgos/src/EgammaEcalIsolation.cc

---