EgammaEcalIsolation Class Reference

#include <EgammaEcalIsolation.h>

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 28 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 basicClusterCollection_, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), conesize, etMin, create_public_lumi_plots::exp, reco::Candidate::get(), edm::Ref< C, T, F >::get(), reco::CaloCluster::position(), position, reco::SuperCluster::seed(), funct::sin(), and superClusterCollection_.

                                                                      {

  
  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::CaloCluster_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 51 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::conesize

private

Definition at line 49 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::etMin

private

Definition at line 48 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

const reco::SuperClusterCollection* EgammaEcalIsolation::superClusterCollection_

private

Definition at line 52 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().