ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/src/ElectronTools.cc
(Generate patch)

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.22 by ceballos, Tue Mar 15 08:33:42 2011 UTC vs.
Revision 1.35 by sixie, Sun Apr 15 12:06:26 2012 UTC

# Line 3 | Line 3
3   #include "MitPhysics/Utils/interface/ElectronTools.h"
4   #include "MitAna/DataTree/interface/StableData.h"
5   #include <TFile.h>
6 + #include "MitPhysics/ElectronLikelihood/interface/ElectronLikelihood.h"
7 + #include "MitPhysics/ElectronLikelihood/interface/LikelihoodSwitches.h"
8 + #include "MitPhysics/ElectronLikelihood/interface/LikelihoodMeasurements.h"
9  
10   ClassImp(mithep::ElectronTools)
11  
# Line 35 | Line 38 | Bool_t ElectronTools::PassCustomID(const
38      {0.3, 0.92, 0.211, 0.0, 0.42, 0.88, 0.68, 0.0},             //eoverp
39      {0.8,0.2,0,0,0,0,0,0}};                                     //extra cuts fbrem and E_Over_P
40  
41 +  Double_t VBTFWorkingPointFakeable[6][8] = {
42 +    {0.12,  0.12,  0.12,  0.12,  0.10,   0.10,   0.10,   0.10  }, //hovere
43 +    {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03  }, //sigmaetaeta
44 +    {0.15,  0.15,  0.15,  0.15,  0.10,   0.10,   0.10,   0.10  }, //deltaphiin
45 +    {0.007, 0.007, 0.007, 0.007, 0.009,  0.009,  0.009,  0.009 }, //deltaetain
46 +    {0.0,   0.0,   0.0,   0.0,   0.0,    0.0,    0.0,    0.0   }, //eoverp
47 +    {0.0,   0.0,   0,     0,     0,      0,      0,      0     }  //extra cuts fbrem and E_Over_P
48 +  };            
49 +
50    Double_t VBTFWorkingPoint95[6][8] = {
51      {0.15,  0.15,  0.15,  0.15,  0.07,   0.07,   0.07,   0.07  }, //hovere
52      {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03  }, //sigmaetaeta
# Line 63 | Line 75 | Bool_t ElectronTools::PassCustomID(const
75    };            
76  
77    Double_t VBTFWorkingPoint80[6][8] = {
78 <    {0.04,  0.04,  0.04,  0.04,  0.025,  0.025,  0.025,  0.025}, //hovere
78 >    {0.04,  0.04,  0.04,  0.04,  0.10,   0.10,   0.10,   0.10 }, //hovere
79      {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
80      {0.06,  0.06,  0.06,  0.06,  0.03,   0.03,   0.03,   0.03 }, //deltaphiin
81      {0.004, 0.004, 0.004, 0.004, 0.007,  0.007,  0.007,  0.007}, //deltaetain
# Line 80 | Line 92 | Bool_t ElectronTools::PassCustomID(const
92      {0.0,   0.0,   0,     0,     0,      0,      0,      0    }  //extra cuts fbrem and E_Over_P
93    };            
94  
95 +  Double_t VBTFWorkingPoint80NoHOverEE[6][8] = {
96 +    {0.04,  0.04,  0.04,  0.04,  0.10,   0.10,   0.10,   0.10 }, //hovere
97 +    {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
98 +    {0.06,  0.06,  0.06,  0.06,  0.03,   0.03,   0.03,   0.03 }, //deltaphiin
99 +    {0.004, 0.004, 0.004, 0.004, 0.007,  0.007,  0.007,  0.007}, //deltaetain
100 +    {0.0,   0.0,   0.0,   0.0,   0.0,    0.0,    0.0,    0.0  }, //eoverp
101 +    {0.0,   0.0,   0,     0,     0,      0,      0,      0    }  //extra cuts fbrem and E_Over_P
102 +  };            
103 +
104 +  Double_t VBTFWorkingPoint70NoHOverEE[6][8] = {
105 +    {0.025, 0.025, 0.025, 0.025, 0.10,   0.10,   0.10,   0.10 }, //hovere
106 +    {0.01,  0.01,  0.01,  0.01,  0.03,   0.03,   0.03,   0.03 }, //sigmaetaeta
107 +    {0.03,  0.03,  0.03,  0.03,  0.02,   0.02,   0.02,   0.02 }, //deltaphiin
108 +    {0.004, 0.004, 0.004, 0.004, 0.005,  0.005,  0.005,  0.005}, //deltaetain
109 +    {0.0,   0.0,   0.0,   0.0,   0.0,    0.0,    0.0,    0.0  }, //eoverp
110 +    {0.0,   0.0,   0,     0,     0,      0,      0,      0    }  //extra cuts fbrem and E_Over_P
111 +  };            
112 +
113    switch (idType) {
114      case kCustomIdTight:    
115        memcpy(fCuts,tightcuts,sizeof(fCuts));
# Line 87 | Line 117 | Bool_t ElectronTools::PassCustomID(const
117      case kCustomIdLoose:
118        memcpy(fCuts,loosecuts,sizeof(fCuts));
119        break;
120 +    case kVBTFWorkingPointFakeableId:
121 +      memcpy(fCuts,VBTFWorkingPointFakeable,sizeof(fCuts));
122 +      break;
123      case kVBTFWorkingPoint95Id:
124        memcpy(fCuts,VBTFWorkingPoint95,sizeof(fCuts));
125        break;
# Line 99 | Line 132 | Bool_t ElectronTools::PassCustomID(const
132      case kVBTFWorkingPoint80Id:
133        memcpy(fCuts,VBTFWorkingPoint80,sizeof(fCuts));
134        break;
135 <    case kVBTFWorkingPoint80LowPtId:
136 <      memcpy(fCuts,VBTFWorkingPoint80,sizeof(fCuts));
135 >    case kVBTFWorkingPointLowPtId:
136 >      if(ele->Pt() < 20)
137 >        memcpy(fCuts,VBTFWorkingPoint70NoHOverEE,sizeof(fCuts));
138 >      else
139 >        memcpy(fCuts,VBTFWorkingPoint80NoHOverEE,sizeof(fCuts));
140        break;
141      case kVBTFWorkingPoint70Id:
142        memcpy(fCuts,VBTFWorkingPoint70,sizeof(fCuts));
# Line 110 | Line 146 | Bool_t ElectronTools::PassCustomID(const
146        break;
147    }
148  
113 –
149    // Based on RecoEgamma/ElectronIdentification/src/CutBasedElectronID.cc.
150    Double_t eOverP = ele->ESuperClusterOverP();
151    Double_t fBrem  = ele->FBrem();
# Line 138 | Line 173 | Bool_t ElectronTools::PassCustomID(const
173    Int_t eb = 1;
174    if (ele->IsEB())
175      eb = 0;
176 <
176 >
177    if (hOverE>fCuts[0][cat+4*eb])
178      return kFALSE;
179  
# Line 154 | Line 189 | Bool_t ElectronTools::PassCustomID(const
189    if(eSeedOverPin<fCuts[4][cat+4*eb])
190      return kFALSE;
191  
192 <  // Cuts only for pt<20 region and kVBTFWorkingPoint80LowPtId
193 <  if(ele->Pt() < 20 && idType == kVBTFWorkingPoint80LowPtId) {
192 >  // Cuts only for pt<20 region and kVBTFWorkingPointLowPtId
193 >  if(ele->Pt() < 20 && idType == kVBTFWorkingPointLowPtId) {
194      Bool_t isGoodLowPtEl = fBrem > 0.15 ||
195 <                          (ele->AbsEta() < 1.0 && eOverP > 0.95 &&
161 <                           TMath::Abs(ele->Charge()*ele->DeltaPhiSuperClusterTrackAtVtx()) < 0.006);
195 >                          (ele->SCluster()->AbsEta() < 1.0 && eOverP > 0.95);
196      if(!isGoodLowPtEl) return kFALSE;
197    }
198    
# Line 263 | Line 297 | Bool_t ElectronTools::PassConversionFilt
297                                             Double_t probMin,
298                                             Double_t lxyMin,
299                                             Bool_t matchCkf,
300 <                                           Bool_t requireArbitratedMerged)
300 >                                           Bool_t requireArbitratedMerged,
301 >                                           Double_t trkptMin)
302   {
303    Bool_t isGoodConversion = kFALSE;
304  
# Line 289 | Line 324 | Bool_t ElectronTools::PassConversionFilt
324            const Track *trk = dynamic_cast<const ChargedParticle*>
325              (conversions->At(ifc)->Daughter(d))->Trk();
326            if (trk) {
327 +            if (trk->Pt()<trkptMin) isGoodConversion = kFALSE;
328              const StableData *sd = dynamic_cast<const StableData*>
329                (conversions->At(ifc)->DaughterDat(d));
330              if (sd->NWrongHits() > nWrongHitsMax)
# Line 308 | Line 344 | Bool_t ElectronTools::PassConversionFilt
344   }
345  
346   //--------------------------------------------------------------------------------------------------
347 < Bool_t ElectronTools::PassD0Cut(const Electron *ele, const VertexCol *vertices, Double_t fD0Cut)
347 > Bool_t ElectronTools::PassD0Cut(const Electron *ele, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
348   {
349    Bool_t d0cut = kFALSE;
350 <  // d0 cut
351 <  Double_t d0_real = 99999;
352 <  for(UInt_t i0 = 0; i0 < vertices->GetEntries(); i0++) {
353 <    if(vertices->At(i0)->NTracks() > 0){
354 <      Double_t pD0 = ele->GsfTrk()->D0Corrected(*vertices->At(i0));
355 <      d0_real = TMath::Abs(pD0);
356 <      break;
350 >
351 >  Double_t d0_real = 1e30;
352 >  if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
353 >  else            {
354 >    Double_t distVtx = 999.0;
355 >    Int_t closestVtx = 0;
356 >    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
357 >      double dz = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nv)));
358 >      if(dz < distVtx) {
359 >        distVtx    = dz;
360 >        closestVtx = nv;
361 >      }
362      }
363 +    d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(closestVtx)));
364    }
365    if(d0_real < fD0Cut) d0cut = kTRUE;
366    
# Line 341 | Line 383 | Bool_t ElectronTools::PassD0Cut(const El
383   }
384  
385   //--------------------------------------------------------------------------------------------------
386 + Bool_t ElectronTools::PassDZCut(const Electron *ele, const VertexCol *vertices, Double_t fDZCut, Int_t nVertex)
387 + {
388 +  Bool_t dzcut = kFALSE;
389 +
390 +  Double_t distVtx = 999.0;
391 +  if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
392 +  else {
393 +    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
394 +      double dz = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nv)));
395 +      if(dz < distVtx) {
396 +        distVtx  = dz;
397 +      }
398 +    }
399 +  }
400 +
401 +  if(distVtx < fDZCut) dzcut = kTRUE;
402 +  
403 +  return dzcut;
404 + }
405 +
406 + //--------------------------------------------------------------------------------------------------
407   Bool_t ElectronTools::PassChargeFilter(const Electron *ele)
408   {
409    Bool_t passChargeFilter = kTRUE;
# Line 950 | Line 1013 | bool ElectronTools::compute_cut(double x
1013  
1014    return accept;
1015   }
1016 +
1017 + //--------------------------------------------------------------------------------------------------
1018 + Double_t ElectronTools::Likelihood(ElectronLikelihood *LH, const Electron *ele)
1019 + {
1020 +  if (!LH) {
1021 +    std::cout << "Error: Likelihood not properly initialized\n";
1022 +    return -9999;
1023 +  }
1024 +
1025 +  LikelihoodMeasurements measurements;
1026 +  measurements.pt = ele->Pt();
1027 +  if (ele->IsEB() && ele->AbsEta()<1.0) measurements.subdet = 0;
1028 +  else if (ele->IsEB())                 measurements.subdet = 1;
1029 +  else                                  measurements.subdet = 2;
1030 +  measurements.deltaPhi = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1031 +  measurements.deltaEta = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1032 +  measurements.eSeedClusterOverPout = ele->ESeedClusterOverPout();
1033 +  measurements.eSuperClusterOverP = ele->ESuperClusterOverP();
1034 +  measurements.hadronicOverEm = ele->HadronicOverEm();
1035 +  measurements.sigmaIEtaIEta = ele->CoviEtaiEta();
1036 +  measurements.sigmaIPhiIPhi = TMath::Sqrt(ele->SCluster()->Seed()->CoviPhiiPhi());
1037 +  measurements.fBrem = ele->FBrem();
1038 +  measurements.nBremClusters = ele->NumberOfClusters() - 1;
1039 +  //measurements.OneOverEMinusOneOverP = (1.0 / ele->SCluster()->Energy()) - (1.0 / ele->BestTrk()->P());
1040 +  measurements.OneOverEMinusOneOverP = (1.0 / ele->ESuperClusterOverP() / ele->BestTrk()->P()) - (1.0 / ele->BestTrk()->P());
1041 +  double likelihood = LH->result(measurements);
1042 +
1043 +  double newLik = 0.0;
1044 +  if     (likelihood<=0) newLik = -20.0;
1045 +  else if(likelihood>=1) newLik =  20.0;
1046 +  else                   newLik = log(likelihood/(1.0-likelihood));
1047 +
1048 +  Bool_t isDebug = kFALSE;
1049 +  if(isDebug == kTRUE){
1050 +    printf("LIKELIHOOD: %f %d %f %f %f %f %f %f %f %f %d %f %f %f - %f %f\n",measurements.pt,measurements.subdet,
1051 +    measurements.deltaPhi          ,measurements.deltaEta      ,measurements.eSeedClusterOverPout,
1052 +    measurements.eSuperClusterOverP,measurements.hadronicOverEm,measurements.sigmaIEtaIEta,
1053 +    measurements.sigmaIPhiIPhi     ,measurements.fBrem         ,measurements.nBremClusters,
1054 +    measurements.OneOverEMinusOneOverP,ele->SCluster()->Energy(),ele->BestTrk()->P(),
1055 +    likelihood,newLik);
1056 +  }
1057 +
1058 +  return newLik;
1059 +
1060 + }
1061 +
1062 + Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
1063 +
1064 +  Double_t EffectiveArea = 0;
1065 +  
1066 +  if (fabs(SCEta) < 1.0) {
1067 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1068 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
1069 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.045;
1070 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
1071 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1072 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.034;
1073 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.079;
1074 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
1075 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1076 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00016;
1077 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
1078 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1079 +  } else if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) {
1080 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1081 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.025;
1082 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.052;
1083 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
1084 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1085 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.050;
1086 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.073;
1087 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
1088 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1089 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00022;
1090 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
1091 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1092 +  } else if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) {
1093 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1094 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.030;
1095 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.170;
1096 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
1097 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1098 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.060;
1099 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.187;
1100 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
1101 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1102 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00030;
1103 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
1104 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;        
1105 +  } else if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.25 ) {
1106 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1107 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.022;
1108 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.623;
1109 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
1110 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1111 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.055;
1112 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.659;
1113 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
1114 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1115 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00054;
1116 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
1117 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
1118 +  } else if (fabs(SCEta) >= 2.25 && fabs(SCEta) < 2.5 ) {
1119 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1120 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.018;
1121 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 1.198;
1122 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
1123 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1124 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.073;
1125 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 1.258;
1126 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
1127 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1128 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00082;
1129 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
1130 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
1131 +  }
1132 +    
1133 +  //NoCorrections
1134 +  if (EffectiveAreaTarget == kEleEANoCorr) {
1135 +    return 0.0;
1136 +  }
1137 +  //2011 Data Effective Areas
1138 +  else if (EffectiveAreaTarget == kEleEAData2011) {
1139 +    if (type == kEleGammaIsoDR0p0To0p1) {
1140 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1141 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
1142 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
1143 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1144 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.004;
1145 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1146 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1147 +    }
1148 +    if (type == kEleGammaIsoDR0p1To0p2) {
1149 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1150 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1151 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
1152 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1153 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
1154 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.035;
1155 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.041;
1156 +    }
1157 +    if (type == kEleGammaIsoDR0p2To0p3) {
1158 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1159 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1160 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.014;
1161 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1162 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.039;
1163 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.042;
1164 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.048;
1165 +    }
1166 +    if (type == kEleGammaIsoDR0p3To0p4) {
1167 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1168 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.029;
1169 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1170 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1171 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
1172 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.047;
1173 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.054;
1174 +    }
1175 +    if (type == kEleGammaIsoDR0p4To0p5) {
1176 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1177 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
1178 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
1179 +     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1180 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.047;
1181 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.057;
1182 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1183 +    }
1184 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1185 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1186 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
1187 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
1188 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1189 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1190 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1191 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1192 +    }
1193 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1194 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
1195 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1196 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.008;
1197 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1198 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1199 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.001;
1200 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.003;
1201 +    }
1202 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1203 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1204 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1205 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1206 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1207 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1208 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1209 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.019;
1210 +    }
1211 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1212 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1213 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1214 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.025;
1215 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.030;
1216 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.036;
1217 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.038;
1218 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.084;
1219 +    }
1220 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1221 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1222 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1223 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
1224 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1225 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1226 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.107;
1227 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.228;
1228 +    }
1229 +  }
1230 +
1231 +  //Summer11 MC Effective Areas
1232 +  else if (EffectiveAreaTarget == kEleEASummer11MC) {
1233 +    if (type == kEleGammaIsoDR0p0To0p1) {
1234 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1235 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1236 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1237 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1238 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1239 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.024;
1240 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.023;
1241 +    }
1242 +    if (type == kEleGammaIsoDR0p1To0p2) {
1243 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1244 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1245 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1246 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.037;
1247 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.046;
1248 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.055;
1249 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.046;
1250 +    }
1251 +    if (type == kEleGammaIsoDR0p2To0p3) {
1252 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.021;
1253 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
1254 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1255 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.026;
1256 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1257 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.045;
1258 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1259 +    }
1260 +    if (type == kEleGammaIsoDR0p3To0p4) {
1261 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1262 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1263 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1264 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1265 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.058;
1266 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.073;
1267 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.083;
1268 +    }
1269 +    if (type == kEleGammaIsoDR0p4To0p5) {
1270 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.053;
1271 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
1272 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
1273 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.048;
1274 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.062;
1275 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.085;
1276 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.118;
1277 +    }
1278 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1279 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
1280 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
1281 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1282 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1283 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1284 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1285 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1286 +    }
1287 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1288 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
1289 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
1290 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1291 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
1292 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1293 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1294 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.004;
1295 +    }
1296 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1297 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.008;
1298 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
1299 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1300 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.013;
1301 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1302 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1303 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1304 +    }
1305 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1306 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1307 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1308 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1309 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1310 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.040;
1311 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.036;
1312 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.086;
1313 +    }
1314 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1315 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
1316 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
1317 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
1318 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1319 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1320 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.105;
1321 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.169;
1322 +    }
1323 +  }
1324 +  
1325 +  //Fall11 MC Effective Areas
1326 +  else if (EffectiveAreaTarget == kEleEAFall11MC) {
1327 +    if (type == kEleGammaIsoDR0p0To0p1) {
1328 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.014;
1329 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1330 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1331 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
1332 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1333 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.021;
1334 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.012;
1335 +    }
1336 +    if (type == kEleGammaIsoDR0p1To0p2) {
1337 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1338 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1339 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
1340 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1341 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.055;
1342 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.068;
1343 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.067;
1344 +    }
1345 +    if (type == kEleGammaIsoDR0p2To0p3) {
1346 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.024;
1347 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1348 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1349 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1350 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1351 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.066;
1352 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.080;
1353 +    }
1354 +    if (type == kEleGammaIsoDR0p3To0p4) {
1355 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.040;
1356 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1357 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1358 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.047;
1359 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1360 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.083;
1361 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.123;
1362 +    }
1363 +    if (type == kEleGammaIsoDR0p4To0p5) {
1364 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1365 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.041;
1366 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.037;
1367 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.057;
1368 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.095;
1369 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1370 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.133;
1371 +    }
1372 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1373 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1374 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.003;
1375 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1376 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1377 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1378 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1379 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1380 +    }
1381 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1382 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.006;
1383 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1384 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
1385 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1386 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1387 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.002;
1388 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.007;
1389 +    }
1390 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1391 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
1392 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1393 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.018;
1394 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1395 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
1396 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1397 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1398 +    }
1399 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1400 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1401 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
1402 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1403 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.035;
1404 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
1405 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1406 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.110;
1407 +    }
1408 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1409 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1410 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1411 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
1412 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1413 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.057;
1414 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1415 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.220;
1416 +    }
1417 +  }
1418 +
1419 +
1420 +
1421 +  return EffectiveArea;  
1422 + }
1423 +
1424 +

Diff Legend

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