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 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 22 of file EgammaEcalIsolation.cc.

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

EgammaEcalIsolation::~EgammaEcalIsolation

(

)

Definition at line 31 of file EgammaEcalIsolation.cc.

{}

Member Function Documentation

double EgammaEcalIsolation::getEcalEtSum

(

const reco::Candidate *

candidate

)

Definition at line 33 of file EgammaEcalIsolation.cc.

References basicClusterCollection_, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), conesize, electronAnalyzer_cfi::DeltaR, etMin, JetChargeProducer_cfi::exp, reco::Candidate::get(), edm::Ref< C, T, F >::get(), reco::CaloCluster::position(), position, SimDataFormats::CaloAnalysis::sc, 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 = nullptr;
  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 = nullptr;

  //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 45 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::conesize

private

Definition at line 43 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

double EgammaEcalIsolation::etMin

private

Definition at line 42 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().

const reco::SuperClusterCollection* EgammaEcalIsolation::superClusterCollection_

private

Definition at line 46 of file EgammaEcalIsolation.h.

Referenced by getEcalEtSum().