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.24 by ceballos, Mon Mar 21 12:25:36 2011 UTC vs.
Revision 1.49 by mingyang, Mon Oct 8 17:22:59 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 >  // Apply detector isolation at high pt only
193 >  Bool_t isoCut = kTRUE;
194 >  if(idType == kVBTFWorkingPointFakeableId){
195 >    double isoEcal = ele->EcalRecHitIsoDr03();
196 >    if(ele->IsEB()) isoEcal = isoEcal - 1.0;
197 >    isoCut = (ele->TrackIsolationDr03() < ele->Pt()*0.2) &&
198 >             (isoEcal                   < ele->Pt()*0.2) &&
199 >             (ele->HcalTowerSumEtDr03() < ele->Pt()*0.2);
200 >  }
201 >  if(isoCut == kFALSE) return kFALSE;
202 >
203 >  // Cuts only for pt<20 region and kVBTFWorkingPointLowPtId
204 >  if(ele->Pt() < 20 && idType == kVBTFWorkingPointLowPtId) {
205      Bool_t isGoodLowPtEl = fBrem > 0.15 ||
206 <                          (ele->AbsEta() < 1.0 && eOverP > 0.95);
206 >                          (ele->SCluster()->AbsEta() < 1.0 && eOverP > 0.95);
207      if(!isGoodLowPtEl) return kFALSE;
208    }
209    
# Line 262 | Line 308 | Bool_t ElectronTools::PassConversionFilt
308                                             Double_t probMin,
309                                             Double_t lxyMin,
310                                             Bool_t matchCkf,
311 <                                           Bool_t requireArbitratedMerged)
311 >                                           Bool_t requireArbitratedMerged,
312 >                                           Double_t trkptMin)
313   {
314    Bool_t isGoodConversion = kFALSE;
315  
# Line 288 | Line 335 | Bool_t ElectronTools::PassConversionFilt
335            const Track *trk = dynamic_cast<const ChargedParticle*>
336              (conversions->At(ifc)->Daughter(d))->Trk();
337            if (trk) {
338 +            if (trk->Pt()<trkptMin) isGoodConversion = kFALSE;
339              const StableData *sd = dynamic_cast<const StableData*>
340                (conversions->At(ifc)->DaughterDat(d));
341              if (sd->NWrongHits() > nWrongHitsMax)
# Line 312 | Line 360 | Bool_t ElectronTools::PassD0Cut(const El
360    Bool_t d0cut = kFALSE;
361  
362    Double_t d0_real = 1e30;
363 +
364 +  if( nVertex >= (int) vertices->GetEntries() )
365 +    nVertex = vertices->GetEntries() - 1;
366 +  
367    if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
368    else            {
369      Double_t distVtx = 999.0;
# Line 346 | Line 398 | Bool_t ElectronTools::PassD0Cut(const El
398   }
399  
400   //--------------------------------------------------------------------------------------------------
401 < Bool_t ElectronTools::PassDZCut(const Electron *ele, const VertexCol *vertices, Double_t fDZCut)
401 > Bool_t ElectronTools::PassDZCut(const Electron *ele, const VertexCol *vertices, Double_t fDZCut, Int_t nVertex)
402   {
403    Bool_t dzcut = kFALSE;
404  
405    Double_t distVtx = 999.0;
406 <  for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
407 <    double dz = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nv)));
408 <    if(dz < distVtx) {
409 <      distVtx    = dz;
406 >
407 >  if( nVertex >= (int) vertices->GetEntries() )
408 >    nVertex = vertices->GetEntries()-1;
409 >  
410 >  if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
411 >  else {
412 >    for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
413 >      double dz = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nv)));
414 >      if(dz < distVtx) {
415 >        distVtx  = dz;
416 >      }
417      }
418    }
419 +  
420    if(distVtx < fDZCut) dzcut = kTRUE;
421    
422    return dzcut;
# Line 972 | Line 1032 | bool ElectronTools::compute_cut(double x
1032  
1033    return accept;
1034   }
1035 +
1036 + //--------------------------------------------------------------------------------------------------
1037 + Double_t ElectronTools::Likelihood(ElectronLikelihood *LH, const Electron *ele)
1038 + {
1039 +  if (!LH) {
1040 +    std::cout << "Error: Likelihood not properly initialized\n";
1041 +    return -9999;
1042 +  }
1043 +
1044 +  LikelihoodMeasurements measurements;
1045 +  measurements.pt = ele->Pt();
1046 +  if (ele->IsEB() && ele->AbsEta()<1.0) measurements.subdet = 0;
1047 +  else if (ele->IsEB())                 measurements.subdet = 1;
1048 +  else                                  measurements.subdet = 2;
1049 +  measurements.deltaPhi = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1050 +  measurements.deltaEta = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1051 +  measurements.eSeedClusterOverPout = ele->ESeedClusterOverPout();
1052 +  measurements.eSuperClusterOverP = ele->ESuperClusterOverP();
1053 +  measurements.hadronicOverEm = ele->HadronicOverEm();
1054 +  measurements.sigmaIEtaIEta = ele->CoviEtaiEta();
1055 +  measurements.sigmaIPhiIPhi = TMath::Sqrt(ele->SCluster()->Seed()->CoviPhiiPhi());
1056 +  measurements.fBrem = ele->FBrem();
1057 +  measurements.nBremClusters = ele->NumberOfClusters() - 1;
1058 +  //measurements.OneOverEMinusOneOverP = (1.0 / ele->SCluster()->Energy()) - (1.0 / ele->BestTrk()->P());
1059 +  measurements.OneOverEMinusOneOverP = (1.0 / ele->ESuperClusterOverP() / ele->BestTrk()->P()) - (1.0 / ele->BestTrk()->P());
1060 +  double likelihood = LH->result(measurements);
1061 +
1062 +  double newLik = 0.0;
1063 +  if     (likelihood<=0) newLik = -20.0;
1064 +  else if(likelihood>=1) newLik =  20.0;
1065 +  else                   newLik = log(likelihood/(1.0-likelihood));
1066 +
1067 +  Bool_t isDebug = kFALSE;
1068 +  if(isDebug == kTRUE){
1069 +    printf("LIKELIHOOD: %f %d %f %f %f %f %f %f %f %f %d %f %f %f - %f %f\n",measurements.pt,measurements.subdet,
1070 +    measurements.deltaPhi          ,measurements.deltaEta      ,measurements.eSeedClusterOverPout,
1071 +    measurements.eSuperClusterOverP,measurements.hadronicOverEm,measurements.sigmaIEtaIEta,
1072 +    measurements.sigmaIPhiIPhi     ,measurements.fBrem         ,measurements.nBremClusters,
1073 +    measurements.OneOverEMinusOneOverP,ele->SCluster()->Energy(),ele->BestTrk()->P(),
1074 +    likelihood,newLik);
1075 +  }
1076 +
1077 +  return newLik;
1078 +
1079 + }
1080 +
1081 + Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
1082 +
1083 +  Double_t EffectiveArea = 0;
1084 +  
1085 +  if(type == ElectronTools::kEleNeutralIso04){
1086 +    if     (fabs(SCEta) < 1.0                           ) EffectiveArea = 0.19;
1087 +    else if(fabs(SCEta) >= 1.0   && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1088 +    else if(fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.000 ) EffectiveArea = 0.12;
1089 +    else if(fabs(SCEta) >= 2.000 && fabs(SCEta) < 2.200 ) EffectiveArea = 0.21;
1090 +    else if(fabs(SCEta) >= 2.200 && fabs(SCEta) < 2.300 ) EffectiveArea = 0.27;
1091 +    else if(fabs(SCEta) >= 2.300 && fabs(SCEta) < 2.400 ) EffectiveArea = 0.44;
1092 +    else                                                  EffectiveArea = 0.52;
1093 +    return EffectiveArea;
1094 +  }
1095 +
1096 +  if (fabs(SCEta) < 1.0) {
1097 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1098 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
1099 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.045;
1100 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
1101 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1102 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.034;
1103 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.079;
1104 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
1105 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1106 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00016;
1107 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
1108 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1109 +  } else if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) {
1110 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1111 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.025;
1112 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.052;
1113 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
1114 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1115 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.050;
1116 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.073;
1117 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
1118 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1119 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00022;
1120 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
1121 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1122 +  } else if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) {
1123 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1124 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.030;
1125 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.170;
1126 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
1127 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1128 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.060;
1129 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.187;
1130 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
1131 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1132 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00030;
1133 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
1134 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;        
1135 +  } else if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.25 ) {
1136 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1137 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.022;
1138 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.623;
1139 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
1140 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1141 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.055;
1142 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.659;
1143 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
1144 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1145 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00054;
1146 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
1147 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
1148 +  } else if (fabs(SCEta) >= 2.25 && fabs(SCEta) < 2.5 ) {
1149 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1150 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.018;
1151 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 1.198;
1152 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
1153 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1154 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.073;
1155 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 1.258;
1156 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
1157 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1158 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00082;
1159 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
1160 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
1161 +  }
1162 +    
1163 +  //NoCorrections
1164 +  if (EffectiveAreaTarget == kEleEANoCorr) {
1165 +    return 0.0;
1166 +  }
1167 +  else if (EffectiveAreaTarget == kEleEAData2012) {
1168 +    if (type == kEleGammaIso03) {
1169 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.122;
1170 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.147;
1171 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.055;
1172 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.106;
1173 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.138;
1174 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.221;
1175 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.211;
1176 +    }
1177 +    if (type == kEleGammaIso04) {
1178 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.176;
1179 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.206;
1180 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.094;
1181 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.172;
1182 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.244;
1183 +      if (fabs(SCEta) >= 2.4 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.333;
1184 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.348;
1185 +    }
1186 +    if (type == kEleNeutralHadronIso03) {
1187 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1188 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1189 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1190 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1191 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1192 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1193 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.019;
1194 +    }
1195 +    if (type == kEleNeutralHadronIso04) {
1196 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.022;
1197 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.036;
1198 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1199 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
1200 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.052;
1201 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.063;
1202 +      if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.028;
1203 +    }
1204 +    if (type == kEleGammaAndNeutralHadronIso04) {
1205 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
1206 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1207 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.12;
1208 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.21;
1209 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.27;
1210 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.44;
1211 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.52;
1212 +    }                          
1213 +    if (type == kEleGammaIsoDR0p0To0p1) {
1214 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1215 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1216 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
1217 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1218 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1219 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.013;
1220 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.013;
1221 +    }
1222 +    if (type == kEleGammaIsoDR0p1To0p2) {
1223 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1224 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
1225 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1226 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
1227 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1228 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1229 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.102;
1230 +    }
1231 +    if (type == kEleGammaIsoDR0p2To0p3) {
1232 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
1233 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1234 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
1235 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
1236 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
1237 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.034;
1238 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.042;
1239 +    }
1240 +    if (type == kEleGammaIsoDR0p3To0p4) {
1241 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.039;
1242 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1243 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1244 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1245 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.053;
1246 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.059;
1247 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.069;
1248 +    }
1249 +    if (type == kEleGammaIsoDR0p4To0p5) {
1250 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1251 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.045;
1252 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
1253 +     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.043;
1254 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.056;
1255 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.065;
1256 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.074;
1257 +    }
1258 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1259 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1260 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.006;
1261 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
1262 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1263 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.001;
1264 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.003;
1265 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.008;
1266 +    }
1267 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1268 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1269 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
1270 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1271 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.003;
1272 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1273 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.006;
1274 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.010;
1275 +    }
1276 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1277 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.007;
1278 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1279 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1280 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
1281 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.007;
1282 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.018;
1283 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.028;
1284 +    }
1285 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1286 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1287 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1288 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
1289 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
1290 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.018;
1291 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.026;
1292 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.063;
1293 +    }
1294 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1295 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
1296 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.012;
1297 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1298 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
1299 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1300 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.051;
1301 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.143;
1302 +    }
1303 +  }
1304 +
1305 +  //2011 Data Effective Areas
1306 +  else if (EffectiveAreaTarget == kEleEAData2011) {
1307 +    if (type == kEleGammaAndNeutralHadronIso03) {
1308 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.100;
1309 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.120;
1310 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.085;
1311 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.110;
1312 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.120;
1313 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.120;
1314 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.130;
1315 +    }
1316 +    if (type == kEleGammaAndNeutralHadronIso04) {
1317 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.180;
1318 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.200;
1319 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.150;
1320 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.190;
1321 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.210;
1322 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.220;
1323 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.290;
1324 +    }
1325 +    if (type == kEleGammaIsoDR0p0To0p1) {
1326 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1327 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
1328 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
1329 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1330 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.004;
1331 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1332 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1333 +    }
1334 +    if (type == kEleGammaIsoDR0p1To0p2) {
1335 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1336 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1337 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
1338 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1339 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
1340 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.035;
1341 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.041;
1342 +    }
1343 +    if (type == kEleGammaIsoDR0p2To0p3) {
1344 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1345 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1346 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.014;
1347 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1348 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.039;
1349 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.042;
1350 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.048;
1351 +    }
1352 +    if (type == kEleGammaIsoDR0p3To0p4) {
1353 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1354 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.029;
1355 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1356 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1357 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
1358 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.047;
1359 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.054;
1360 +    }
1361 +    if (type == kEleGammaIsoDR0p4To0p5) {
1362 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1363 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
1364 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
1365 +     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1366 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.047;
1367 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.057;
1368 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1369 +    }
1370 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1371 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1372 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
1373 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
1374 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1375 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1376 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1377 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1378 +    }
1379 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1380 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
1381 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1382 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.008;
1383 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1384 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1385 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.001;
1386 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.003;
1387 +    }
1388 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1389 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1390 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1391 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1392 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1393 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1394 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1395 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.019;
1396 +    }
1397 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1398 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1399 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1400 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.025;
1401 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.030;
1402 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.036;
1403 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.038;
1404 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.084;
1405 +    }
1406 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1407 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1408 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1409 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
1410 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1411 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1412 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.107;
1413 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.228;
1414 +    }
1415 +  }
1416 +
1417 +  //Summer11 MC Effective Areas
1418 +  else if (EffectiveAreaTarget == kEleEASummer11MC) {
1419 +    if (type == kEleGammaIsoDR0p0To0p1) {
1420 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1421 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1422 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1423 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1424 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1425 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.024;
1426 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.023;
1427 +    }
1428 +    if (type == kEleGammaIsoDR0p1To0p2) {
1429 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1430 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1431 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1432 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.037;
1433 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.046;
1434 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.055;
1435 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.046;
1436 +    }
1437 +    if (type == kEleGammaIsoDR0p2To0p3) {
1438 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.021;
1439 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
1440 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1441 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.026;
1442 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1443 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.045;
1444 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1445 +    }
1446 +    if (type == kEleGammaIsoDR0p3To0p4) {
1447 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1448 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1449 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1450 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1451 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.058;
1452 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.073;
1453 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.083;
1454 +    }
1455 +    if (type == kEleGammaIsoDR0p4To0p5) {
1456 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.053;
1457 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
1458 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
1459 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.048;
1460 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.062;
1461 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.085;
1462 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.118;
1463 +    }
1464 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1465 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
1466 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
1467 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1468 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1469 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1470 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1471 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1472 +    }
1473 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1474 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
1475 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
1476 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1477 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
1478 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1479 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1480 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.004;
1481 +    }
1482 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1483 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.008;
1484 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
1485 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1486 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.013;
1487 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1488 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1489 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1490 +    }
1491 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1492 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1493 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1494 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1495 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1496 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.040;
1497 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.036;
1498 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.086;
1499 +    }
1500 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1501 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
1502 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
1503 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
1504 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1505 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1506 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.105;
1507 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.169;
1508 +    }
1509 +  }
1510 +  
1511 +  //Fall11 MC Effective Areas
1512 +  else if (EffectiveAreaTarget == kEleEAFall11MC) {
1513 +    if (type == kEleGammaIsoDR0p0To0p1) {
1514 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.014;
1515 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1516 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1517 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
1518 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1519 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.021;
1520 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.012;
1521 +    }
1522 +    if (type == kEleGammaIsoDR0p1To0p2) {
1523 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1524 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1525 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
1526 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1527 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.055;
1528 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.068;
1529 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.067;
1530 +    }
1531 +    if (type == kEleGammaIsoDR0p2To0p3) {
1532 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.024;
1533 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1534 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1535 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1536 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1537 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.066;
1538 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.080;
1539 +    }
1540 +    if (type == kEleGammaIsoDR0p3To0p4) {
1541 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.040;
1542 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1543 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1544 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.047;
1545 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1546 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.083;
1547 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.123;
1548 +    }
1549 +    if (type == kEleGammaIsoDR0p4To0p5) {
1550 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1551 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.041;
1552 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.037;
1553 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.057;
1554 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.095;
1555 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1556 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.133;
1557 +    }
1558 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1559 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1560 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.003;
1561 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1562 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1563 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1564 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1565 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1566 +    }
1567 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1568 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.006;
1569 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1570 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
1571 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1572 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1573 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.002;
1574 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.007;
1575 +    }
1576 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1577 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
1578 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1579 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.018;
1580 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1581 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
1582 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1583 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1584 +    }
1585 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1586 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1587 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
1588 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1589 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.035;
1590 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
1591 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1592 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.110;
1593 +    }
1594 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1595 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1596 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1597 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
1598 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1599 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.057;
1600 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1601 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.220;
1602 +    }
1603 +  }
1604 +
1605 +  return EffectiveArea;  
1606 + }
1607 +
1608 +
1609 + Bool_t ElectronTools::PassHggLeptonTagID(const Electron* ele) {
1610 +  
1611 +  float dist = ( ele->ConvPartnerDist()      == -9999.? 9999:TMath::Abs(ele->ConvPartnerDist()));
1612 +  float dcot = ( ele->ConvPartnerDCotTheta() == -9999.? 9999:TMath::Abs(ele->ConvPartnerDCotTheta()));  
1613 +  
1614 +  if (dist < 0.02) return false;
1615 +  if (dcot < 0.02) return false;
1616 +  
1617 +  int numInnerHits = ele->Trk()->NExpectedHitsInner();
1618 +  if( numInnerHits > 1 ) return false;
1619 +
1620 +  float coviEtaiEta = ele->CoviEtaiEta();
1621 +  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1622 +  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.031 ) return false;
1623 +
1624 +  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1625 +  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1626 +
1627 +  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.039 || deltaEtaIn > 0.005 ) ) return false;
1628 +  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.028 || deltaEtaIn > 0.007 ) ) return false;
1629 +
1630 +  return true;
1631 + }
1632 +
1633 + Bool_t ElectronTools::PassHggLeptonTagID2012(const Electron* ele) {
1634 +  
1635 +  if (TMath::Abs(1./ele->E()-1./ele->Pt())>0.05) return false;
1636 +
1637 +  int numInnerHits = ele->Trk()->NExpectedHitsInner();
1638 +  if( numInnerHits > 1 ) return false;
1639 +
1640 +  float coviEtaiEta = ele->CoviEtaiEta();
1641 +  if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1642 +  else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.03 ) return false; // h
1643 +
1644 +  Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1645 +  Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1646 +
1647 +  if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.15 || deltaEtaIn > 0.007 || ele->HadronicOverEm()>0.12) ) return false;   // h
1648 +  else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.10 || deltaEtaIn > 0.009 || ele->HadronicOverEm()>0.10 ) ) return false;   // h
1649 +
1650 +  return true;
1651 + }
1652 +

Diff Legend

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