[ViewVC] Diff of: cvsroot/UserCode/MitPhysics/Utils/src/ElectronTools.cc

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs.
Revision 1.52 by pharris, Sat Jan 12 11:49:50 2013 UTC

+#include "MitPhysics/Utils/interface/ElectronTools.h"
+#include "MitAna/DataTree/interface/StableData.h"
+#include <TFile.h>
-+
+#include "MitPhysics/ElectronLikelihood/interface/ElectronLikelihood.h"
-+
+#include "MitPhysics/ElectronLikelihood/interface/LikelihoodSwitches.h"
-+
+#include "MitPhysics/ElectronLikelihood/interface/LikelihoodMeasurements.h"
+ClassImp(mithep::ElectronTools)
+  };
+  Double_t VBTFWorkingPoint80[6][8] = {
-<
+    {0.04,  0.04,  0.04,  0.04,  0.025,  0.025,  0.025,  0.025}, //hovere
->
+    {0.04,  0.04,  0.04,  0.04,  0.10,   0.10,   0.10,   0.10 }, //hovere
+    {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
+    {0.06,  0.06,  0.06,  0.06,  0.03,   0.03,   0.03,   0.03 }, //deltaphiin
+    {0.004, 0.004, 0.004, 0.004, 0.007,  0.007,  0.007,  0.007}, //deltaetain
+  if(eSeedOverPin<fCuts[4][cat+4*eb])
+    return kFALSE;
-+
+  // Apply detector isolation at high pt only
-+
+  Bool_t isoCut = kTRUE;
-+
+  if(idType == kVBTFWorkingPointFakeableId){
-+
+    double isoEcal = ele->EcalRecHitIsoDr03();
-+
+    if(ele->IsEB()) isoEcal = isoEcal - 1.0;
-+
+    isoCut = (ele->TrackIsolationDr03() < ele->Pt()*0.2) &&
-+
+             (isoEcal                   < ele->Pt()*0.2) &&
-+
+             (ele->HcalTowerSumEtDr03() < ele->Pt()*0.2);
-+
+  }
-+
+  if(isoCut == kFALSE) return kFALSE;
-+
+  // Cuts only for pt<20 region and kVBTFWorkingPointLowPtId
+  if(ele->Pt() < 20 && idType == kVBTFWorkingPointLowPtId) {
+    Bool_t isGoodLowPtEl = fBrem > 0.15 ||
+//--------------------------------------------------------------------------------------------------
-+
+Bool_t ElectronTools::PassConversionFilterPFAOD(const Electron *ele,
-+
+                                           const DecayParticleCol *conversions,
-+
+                                           const BaseVertex *vtx,
-+
+                                           UInt_t nWrongHitsMax,
-+
+                                           Double_t probMin,
-+
+                                           Double_t lxyMin,
-+
+                                           Bool_t matchCkf,
-+
+                                           Bool_t requireArbitratedMerged,
-+
+                                           Double_t trkptMin)
-+
+{
-+
-+
+  Bool_t isGoodConversion = kFALSE;
-+
-+
+  for (UInt_t ifc=0; ifc<conversions->GetEntries(); ifc++) {
-+
+    Bool_t ConversionMatchFound = kFALSE;
-+
+    for (UInt_t d=0; d<conversions->At(ifc)->NDaughters(); d++) {
-+
+      const ChargedParticle *pParticle = 0;
-+
+      pParticle = dynamic_cast<const ChargedParticle*>(conversions->At(ifc)->Daughter(d));
-+
+      if(pParticle == 0) continue;
-+
+      const Track* trk = 0;
-+
+      trk = pParticle->Trk();
-+
+      if(trk == 0) continue;
-+
+      if (ele->GsfTrk() == trk || (matchCkf && ele->TrackerTrk()==trk) ) {
-+
+        ConversionMatchFound = kTRUE;
-+
+        break;
-+
+      }
-+
+    }
-+
-+
+    // if match between the e-track and one of the conversion legs
-+
+    if (ConversionMatchFound == kTRUE){
-+
+      isGoodConversion =  (conversions->At(ifc)->Prob() > probMin) &&
-+
+        (!requireArbitratedMerged || conversions->At(ifc)->Quality().Quality(ConversionQuality::arbitratedMerged)) &&
-+
+        (conversions->At(ifc)->LxyCorrected(vtx) > lxyMin);
-+
-+
+      if (isGoodConversion == kTRUE) {
-+
+        for (UInt_t d=0; d<conversions->At(ifc)->NDaughters(); d++) {
-+
+          const ChargedParticle *pParticle = 0;
-+
+          pParticle = dynamic_cast<const ChargedParticle*>(conversions->At(ifc)->Daughter(d));
-+
+          if(pParticle == 0) continue;
-+
+          const Track* trk = 0;
-+
+          trk = pParticle->Trk();
-+
+          if(trk == 0) continue;
-+
+          if (trk) {
-+
+            if (trk->Pt()<trkptMin) isGoodConversion = kFALSE;
-+
+            const StableData *sd = dynamic_cast<const StableData*>
-+
+              (conversions->At(ifc)->DaughterDat(d));
-+
+            if (sd->NWrongHits() > nWrongHitsMax)
-+
+              isGoodConversion = kFALSE;
-+
+          } else {
-+
+            isGoodConversion = kFALSE;
-+
+          }
-+
+        }
-+
+      }
-+
+    }
-+
-+
+    if (isGoodConversion == kTRUE) break;
-+
-+
+  } // loop over all conversions
-+
-+
+  return !isGoodConversion;
-+
+}
-+
+//--------------------------------------------------------------------------------------------------
+Bool_t ElectronTools::PassD0Cut(const Electron *ele, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
+  Bool_t d0cut = kFALSE;
+  Double_t d0_real = 1e30;
-+
-+
+  if( nVertex >= (int) vertices->GetEntries() )
-+
+    nVertex = vertices->GetEntries() - 1;
-+
+  if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
+  else            {
+    Double_t distVtx = 999.0;
+  Bool_t dzcut = kFALSE;
+  Double_t distVtx = 999.0;
-+
-+
+  if( nVertex >= (int) vertices->GetEntries() )
-+
+    nVertex = vertices->GetEntries()-1;
-+
+  if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
+  else {
+    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
-<
->
+  if(distVtx < fDZCut) dzcut = kTRUE;
+  return dzcut;
+  return accept;
-+
-+
+//--------------------------------------------------------------------------------------------------
-+
+Double_t ElectronTools::Likelihood(ElectronLikelihood *LH, const Electron *ele)
-+
+{
-+
+  if (!LH) {
-+
+    std::cout << "Error: Likelihood not properly initialized\n";
-+
+    return -9999;
-+
+  }
-+
-+
+  LikelihoodMeasurements measurements;
-+
+  measurements.pt = ele->Pt();
-+
+  if (ele->IsEB() && ele->AbsEta()<1.0) measurements.subdet = 0;
-+
+  else if (ele->IsEB())                 measurements.subdet = 1;
-+
+  else                                  measurements.subdet = 2;
-+
+  measurements.deltaPhi = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
-+
+  measurements.deltaEta = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
-+
+  measurements.eSeedClusterOverPout = ele->ESeedClusterOverPout();
-+
+  measurements.eSuperClusterOverP = ele->ESuperClusterOverP();
-+
+  measurements.hadronicOverEm = ele->HadronicOverEm();
-+
+  measurements.sigmaIEtaIEta = ele->CoviEtaiEta();
-+
+  measurements.sigmaIPhiIPhi = TMath::Sqrt(ele->SCluster()->Seed()->CoviPhiiPhi());
-+
+  measurements.fBrem = ele->FBrem();
-+
+  measurements.nBremClusters = ele->NumberOfClusters() - 1;
-+
+  //measurements.OneOverEMinusOneOverP = (1.0 / ele->SCluster()->Energy()) - (1.0 / ele->BestTrk()->P());
-+
+  measurements.OneOverEMinusOneOverP = (1.0 / ele->ESuperClusterOverP() / ele->BestTrk()->P()) - (1.0 / ele->BestTrk()->P());
-+
+  double likelihood = LH->result(measurements);
-+
-+
+  double newLik = 0.0;
-+
+  if     (likelihood<=0) newLik = -20.0;
-+
+  else if(likelihood>=1) newLik =  20.0;
-+
+  else                   newLik = log(likelihood/(1.0-likelihood));
-+
-+
+  Bool_t isDebug = kFALSE;
-+
+  if(isDebug == kTRUE){
-+
+    printf("LIKELIHOOD: %f %d %f %f %f %f %f %f %f %f %d %f %f %f - %f %f\n",measurements.pt,measurements.subdet,
-+
+    measurements.deltaPhi          ,measurements.deltaEta      ,measurements.eSeedClusterOverPout,
-+
+    measurements.eSuperClusterOverP,measurements.hadronicOverEm,measurements.sigmaIEtaIEta,
-+
+    measurements.sigmaIPhiIPhi     ,measurements.fBrem         ,measurements.nBremClusters,
-+
+    measurements.OneOverEMinusOneOverP,ele->SCluster()->Energy(),ele->BestTrk()->P(),
-+
+    likelihood,newLik);
-+
+  }
-+
-+
+  return newLik;
-+
-+
+}
-+
-+
+Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
-+
-+
+  Double_t EffectiveArea = 0;
-+
-+
+  if(type == ElectronTools::kEleNeutralIso04){
-+
+    if     (fabs(SCEta) < 1.0                           ) EffectiveArea = 0.19;
-+
+    else if(fabs(SCEta) >= 1.0   && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
-+
+    else if(fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.000 ) EffectiveArea = 0.12;
-+
+    else if(fabs(SCEta) >= 2.000 && fabs(SCEta) < 2.200 ) EffectiveArea = 0.21;
-+
+    else if(fabs(SCEta) >= 2.200 && fabs(SCEta) < 2.300 ) EffectiveArea = 0.27;
-+
+    else if(fabs(SCEta) >= 2.300 && fabs(SCEta) < 2.400 ) EffectiveArea = 0.44;
-+
+    else                                                  EffectiveArea = 0.52;
-+
+    return EffectiveArea;
-+
+  }
-+
-+
+  if (fabs(SCEta) < 1.0) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.045;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.034;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.079;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00016;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
-+
+  } else if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.025;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.052;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.050;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.073;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00022;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
-+
+  } else if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.030;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.170;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.060;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.187;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00030;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;
-+
+  } else if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.25 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.022;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.623;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.055;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.659;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00054;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
-+
+  } else if (fabs(SCEta) >= 2.25 && fabs(SCEta) < 2.5 ) {
-+
+    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.018;
-+
+    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 1.198;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
-+
+    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.073;
-+
+    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 1.258;
-+
+    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
-+
+    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
-+
+    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00082;
-+
+    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
-+
+    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
-+
+  }
-+
-+
+  //NoCorrections
-+
+  if (EffectiveAreaTarget == kEleEANoCorr) {
-+
+    return 0.0;
-+
+  }
-+
+  else if (EffectiveAreaTarget == kEleEAData2012) {
-+
+    if (type == kEleGammaIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.122;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.147;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.106;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.138;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.221;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.211;
-+
+    }
-+
+    if (type == kEleGammaIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.176;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.206;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.094;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.172;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.244;
-+
+      if (fabs(SCEta) >= 2.4 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.333;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.348;
-+
+    }
-+
+    if (type == kEleNeutralHadronIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.019;
-+
+    }
-+
+    if (type == kEleNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.022;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.052;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.063;
-+
+      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.028;
-+
+    }
-+
+    if (type == kEleGammaAndNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.12;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.21;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.27;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.44;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.52;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.013;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.102;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.034;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.042;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.039;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.053;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.059;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.069;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
-+
+     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.056;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.065;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.074;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.008;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.010;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.028;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.063;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.143;
-+
+    }
-+
+  }
-+
-+
+  //2011 Data Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEAData2011) {
-+
+    if (type == kEleGammaAndNeutralHadronIso03) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.100;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.085;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.110;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.120;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.130;
-+
+    }
-+
+    if (type == kEleGammaAndNeutralHadronIso04) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.180;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.200;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.150;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.190;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.210;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.220;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.290;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.041;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.039;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.048;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.054;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
-+
+     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.001;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.003;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.019;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.025;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.084;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.107;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.228;
-+
+    }
-+
+  }
-+
-+
+  //Summer11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEASummer11MC) {
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.023;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.046;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.046;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.058;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.073;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.083;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.053;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.048;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.062;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.085;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.118;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.004;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.040;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.086;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.105;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.169;
-+
+    }
-+
+  }
-+
-+
+  //Fall11 MC Effective Areas
-+
+  else if (EffectiveAreaTarget == kEleEAFall11MC) {
-+
+    if (type == kEleGammaIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.012;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.055;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.068;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.067;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.024;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.080;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.040;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.047;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.083;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.123;
-+
+    }
-+
+    if (type == kEleGammaIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.041;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.095;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.133;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.003;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.002;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.007;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.018;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.035;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.110;
-+
+    }
-+
+    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
-+
+      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
-+
+      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
-+
+      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
-+
+      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
-+
+      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.057;
-+
+      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
-+
+      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.220;
-+
+    }
-+
+  }
-+
-+
+  return EffectiveArea;
-+
+}
-+
-+
-+
+Bool_t ElectronTools::PassHggLeptonTagID(const Electron* ele) {
-+
-+
+  float dist = ( ele->ConvPartnerDist()      == -9999.? 9999:TMath::Abs(ele->ConvPartnerDist()));
-+
+  float dcot = ( ele->ConvPartnerDCotTheta() == -9999.? 9999:TMath::Abs(ele->ConvPartnerDCotTheta()));
-+
-+
+  if (dist < 0.02) return false;
-+
+  if (dcot < 0.02) return false;
-+
-+
+  int numInnerHits = ele->Trk()->NExpectedHitsInner();
-+
+  if( numInnerHits > 1 ) return false;
-+
-+
+  float coviEtaiEta = ele->CoviEtaiEta();
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
-+
+  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.031 ) return false;
-+
-+
+  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
-+
+  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
-+
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.039 || deltaEtaIn > 0.005 ) ) return false;
-+
+  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.028 || deltaEtaIn > 0.007 ) ) return false;
-+
-+
+  return true;
-+
+}
-+
-+
+Bool_t ElectronTools::PassHggLeptonTagID2012(const Electron* ele) {
-+
-+
+  if (TMath::Abs(1./ele->E()-1./ele->Pt())>0.05) return false;
-+
-+
+  int numInnerHits = ele->Trk()->NExpectedHitsInner();
-+
+  if( numInnerHits > 1 ) return false;
-+
-+
+  float coviEtaiEta = ele->CoviEtaiEta();
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
-+
+  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.03 ) return false; // h
-+
-+
+  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
-+
+  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
-+
-+
+  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.15 || deltaEtaIn > 0.007 || ele->HadronicOverEm()>0.12) ) return false;   // h
-+
+  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.10 || deltaEtaIn > 0.009 || ele->HadronicOverEm()>0.10 ) ) return false;   // h
-+
-+
+  return true;
-+
+}
-+
-+
+std::pair<Double_t,Double_t> ElectronTools::ComputeEPCombination( const Electron * ele, const float regression_energy,
-+
+                                                  const float regression_energy_error) {
-+
-+
+  enum Classification { GSF_ELECTRON_UNKNOWN=-1,
-+
+                        GSF_ELECTRON_GOLDEN=0,
-+
+                        GSF_ELECTRON_BIGBREM=1,
-+
+                        GSF_ELECTRON_BADTRACK=2,
-+
+                        GSF_ELECTRON_SHOWERING=3,
-+
+                        GSF_ELECTRON_GAP=4 } ;
-+
-+
+  float newEnergyError_ = regression_energy_error;
-+
+  float scEnergy = regression_energy;
-+
-+
+  int elClass = ele->Classification();
-+
-+
+  float trackMomentum  = ele->PIn() ;
-+
+  float errorTrackMomentum_ = 999. ;
-+
-+
+  // the electron's track momentum error was not available in bambu versions less than 029
-+
+  if(ele->TrackMomentumError() > 0)
-+
+    errorTrackMomentum_ = ele->TrackMomentumError();
-+
+  else if ( ele->GsfTrk()->PtErr() > 0)
-+
+    errorTrackMomentum_ = ele->GsfTrk()->PtErr()*cosh(ele->GsfTrk()->Eta());
-+
+  else
-+
+    assert(0);
-+
-+
+  float finalMomentum = ele->E(); // initial
-+
+  float finalMomentumError = 999.;
-+
-+
+  // first check for large errors
-+
-+
+  if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy <= 0.5) {
-+
+    finalMomentum = regression_energy;    finalMomentumError = regression_energy_error;
-+
+  }
-+
+  else if (errorTrackMomentum_/trackMomentum <= 0.5 && regression_energy_error/regression_energy > 0.5){
-+
+    finalMomentum = trackMomentum;  finalMomentumError = errorTrackMomentum_;
-+
+  }
-+
+  else if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy > 0.5){
-+
+    if (errorTrackMomentum_/trackMomentum < regression_energy_error/regression_energy) {
-+
+      finalMomentum = trackMomentum; finalMomentumError = errorTrackMomentum_;
-+
+    }
-+
+    else{
-+
+      finalMomentum = regression_energy; finalMomentumError = regression_energy_error;
-+
+    }
-+
+  }
-+
+  // then apply the combination algorithm
-+
+  else {
-+
+     // calculate E/p and corresponding error
-+
+    float eOverP = regression_energy / trackMomentum;
-+
+    float errorEOverP = sqrt(
-+
+                             (regression_energy_error/trackMomentum)*(regression_energy_error/trackMomentum) +
-+
+                             (regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum)*
-+
+                             (regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum));
-+
-+
-+
+    bool eleIsNotInCombination = false ;
-+
+    if ( (eOverP  > 1 + 2.5*errorEOverP) || (eOverP  < 1 - 2.5*errorEOverP) || (eOverP < 0.8) || (eOverP > 1.3) )
-+
+      { eleIsNotInCombination = true ; }
-+
+     if (eleIsNotInCombination)
-+
+       {
-+
+         if (eOverP > 1)
-+
+           { finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
-+
+         else
-+
+           {
-+
+             if (elClass == GSF_ELECTRON_GOLDEN)
-+
+               { finalMomentum = regression_energy; finalMomentumError = regression_energy_error; }
-+
+             if (elClass == GSF_ELECTRON_BIGBREM)
-+
+               {
-+
+                 if (regression_energy<36)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
-+
+               }
-+
+             if (elClass == GSF_ELECTRON_BADTRACK)
-+
+               { finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
-+
+             if (elClass == GSF_ELECTRON_SHOWERING)
-+
+               {
-+
+                 if (regression_energy<30)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy; finalMomentumError = regression_energy_error;}
-+
+               }
-+
+             if (elClass == GSF_ELECTRON_GAP)
-+
+               {
-+
+                 if (regression_energy<60)
-+
+                   { finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
-+
+                 else
-+
+                   { finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
-+
+               }
-+
+           }
-+
+       }
-+
+     else
-+
+       {
-+
+         // combination
-+
+         finalMomentum = (regression_energy/regression_energy_error/regression_energy_error + trackMomentum/errorTrackMomentum_/errorTrackMomentum_) /
-+
+           (1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
-+
+         float finalMomentumVariance = 1 / (1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
-+
+         finalMomentumError = sqrt(finalMomentumVariance);
-+
+       }
-+
+  }
-+
-+
+  std::pair<Double_t,Double_t> value;
-+
+  value.first = finalMomentum;
-+
+  value.second = finalMomentumError;
-+
+  return value;
-+
-+
+}

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents): Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs. Revision 1.52 by pharris, Sat Jan 12 11:49:50 2013 UTC

Diff Legend

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.29 by ceballos, Fri Jun 10 10:42:44 2011 UTC vs.
Revision 1.52 by pharris, Sat Jan 12 11:49:50 2013 UTC