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.26 by ceballos, Tue Apr 5 05:40:06 2011 UTC vs.
Revision 1.37 by ceballos, Fri Apr 20 20:25:58 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 81 | Line 93 | Bool_t ElectronTools::PassCustomID(const
93    };            
94  
95    Double_t VBTFWorkingPoint80NoHOverEE[6][8] = {
96 <    {0.04,  0.04,  0.04,  0.04,  1.000,  1.000,  1.000,  1.000}, //hovere
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
# Line 90 | Line 102 | Bool_t ElectronTools::PassCustomID(const
102    };            
103  
104    Double_t VBTFWorkingPoint70NoHOverEE[6][8] = {
105 <    {0.025, 0.025, 0.025, 0.025, 1.000,  1.000,  1.000,  1.000}, //hovere
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
# Line 105 | 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 174 | Line 189 | Bool_t ElectronTools::PassCustomID(const
189    if(eSeedOverPin<fCuts[4][cat+4*eb])
190      return kFALSE;
191  
192 +  Bool_t isoCut = kTRUE;
193 +  if(idType == kVBTFWorkingPointFakeableId){
194 +    double isoEcal = ele->EcalRecHitIsoDr03();
195 +    if(ele->IsEB()) isoEcal = isoEcal - 1.0;
196 +    isoCut = (ele->TrackIsolationDr03() < ele->Pt()*0.2) &&
197 +             (isoEcal                   < ele->Pt()*0.2) &&
198 +             (ele->HcalTowerSumEtDr03() < ele->Pt()*0.2);
199 +  }
200 +  if(isoCut == kFALSE) return kFALSE;
201 +
202    // Cuts only for pt<20 region and kVBTFWorkingPointLowPtId
203    if(ele->Pt() < 20 && idType == kVBTFWorkingPointLowPtId) {
204      Bool_t isGoodLowPtEl = fBrem > 0.15 ||
205 <                          (ele->AbsEta() < 1.0 && eOverP > 0.95);
205 >                          (ele->SCluster()->AbsEta() < 1.0 && eOverP > 0.95);
206      if(!isGoodLowPtEl) return kFALSE;
207    }
208    
# Line 282 | Line 307 | Bool_t ElectronTools::PassConversionFilt
307                                             Double_t probMin,
308                                             Double_t lxyMin,
309                                             Bool_t matchCkf,
310 <                                           Bool_t requireArbitratedMerged)
310 >                                           Bool_t requireArbitratedMerged,
311 >                                           Double_t trkptMin)
312   {
313    Bool_t isGoodConversion = kFALSE;
314  
# Line 308 | Line 334 | Bool_t ElectronTools::PassConversionFilt
334            const Track *trk = dynamic_cast<const ChargedParticle*>
335              (conversions->At(ifc)->Daughter(d))->Trk();
336            if (trk) {
337 +            if (trk->Pt()<trkptMin) isGoodConversion = kFALSE;
338              const StableData *sd = dynamic_cast<const StableData*>
339                (conversions->At(ifc)->DaughterDat(d));
340              if (sd->NWrongHits() > nWrongHitsMax)
# Line 996 | Line 1023 | bool ElectronTools::compute_cut(double x
1023  
1024    return accept;
1025   }
1026 +
1027 + //--------------------------------------------------------------------------------------------------
1028 + Double_t ElectronTools::Likelihood(ElectronLikelihood *LH, const Electron *ele)
1029 + {
1030 +  if (!LH) {
1031 +    std::cout << "Error: Likelihood not properly initialized\n";
1032 +    return -9999;
1033 +  }
1034 +
1035 +  LikelihoodMeasurements measurements;
1036 +  measurements.pt = ele->Pt();
1037 +  if (ele->IsEB() && ele->AbsEta()<1.0) measurements.subdet = 0;
1038 +  else if (ele->IsEB())                 measurements.subdet = 1;
1039 +  else                                  measurements.subdet = 2;
1040 +  measurements.deltaPhi = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1041 +  measurements.deltaEta = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1042 +  measurements.eSeedClusterOverPout = ele->ESeedClusterOverPout();
1043 +  measurements.eSuperClusterOverP = ele->ESuperClusterOverP();
1044 +  measurements.hadronicOverEm = ele->HadronicOverEm();
1045 +  measurements.sigmaIEtaIEta = ele->CoviEtaiEta();
1046 +  measurements.sigmaIPhiIPhi = TMath::Sqrt(ele->SCluster()->Seed()->CoviPhiiPhi());
1047 +  measurements.fBrem = ele->FBrem();
1048 +  measurements.nBremClusters = ele->NumberOfClusters() - 1;
1049 +  //measurements.OneOverEMinusOneOverP = (1.0 / ele->SCluster()->Energy()) - (1.0 / ele->BestTrk()->P());
1050 +  measurements.OneOverEMinusOneOverP = (1.0 / ele->ESuperClusterOverP() / ele->BestTrk()->P()) - (1.0 / ele->BestTrk()->P());
1051 +  double likelihood = LH->result(measurements);
1052 +
1053 +  double newLik = 0.0;
1054 +  if     (likelihood<=0) newLik = -20.0;
1055 +  else if(likelihood>=1) newLik =  20.0;
1056 +  else                   newLik = log(likelihood/(1.0-likelihood));
1057 +
1058 +  Bool_t isDebug = kFALSE;
1059 +  if(isDebug == kTRUE){
1060 +    printf("LIKELIHOOD: %f %d %f %f %f %f %f %f %f %f %d %f %f %f - %f %f\n",measurements.pt,measurements.subdet,
1061 +    measurements.deltaPhi          ,measurements.deltaEta      ,measurements.eSeedClusterOverPout,
1062 +    measurements.eSuperClusterOverP,measurements.hadronicOverEm,measurements.sigmaIEtaIEta,
1063 +    measurements.sigmaIPhiIPhi     ,measurements.fBrem         ,measurements.nBremClusters,
1064 +    measurements.OneOverEMinusOneOverP,ele->SCluster()->Energy(),ele->BestTrk()->P(),
1065 +    likelihood,newLik);
1066 +  }
1067 +
1068 +  return newLik;
1069 +
1070 + }
1071 +
1072 + Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
1073 +
1074 +  Double_t EffectiveArea = 0;
1075 +  
1076 +  if (fabs(SCEta) < 1.0) {
1077 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1078 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
1079 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.045;
1080 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
1081 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1082 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.034;
1083 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.079;
1084 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
1085 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1086 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00016;
1087 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
1088 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1089 +  } else if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) {
1090 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1091 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.025;
1092 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.052;
1093 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
1094 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1095 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.050;
1096 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.073;
1097 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
1098 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1099 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00022;
1100 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
1101 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;    
1102 +  } else if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) {
1103 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1104 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.030;
1105 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.170;
1106 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
1107 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1108 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.060;
1109 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.187;
1110 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
1111 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1112 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00030;
1113 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
1114 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;        
1115 +  } else if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.25 ) {
1116 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1117 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.022;
1118 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 0.623;
1119 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
1120 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1121 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.055;
1122 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 0.659;
1123 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
1124 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1125 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00054;
1126 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
1127 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
1128 +  } else if (fabs(SCEta) >= 2.25 && fabs(SCEta) < 2.5 ) {
1129 +    if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1130 +    if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.018;
1131 +    if (type == ElectronTools::kEleGammaIso03) EffectiveArea = 1.198;
1132 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
1133 +    if (type == ElectronTools::kEleChargedIso04) EffectiveArea = 0.000;
1134 +    if (type == ElectronTools::kEleNeutralHadronIso04) EffectiveArea = 0.073;
1135 +    if (type == ElectronTools::kEleGammaIso04) EffectiveArea = 1.258;
1136 +    if (type == ElectronTools::kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
1137 +    if (type == ElectronTools::kEleNeutralHadronIso007) EffectiveArea = 0.000;
1138 +    if (type == ElectronTools::kEleHoverE) EffectiveArea = 0.00082;
1139 +    if (type == ElectronTools::kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
1140 +    if (type == ElectronTools::kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
1141 +  }
1142 +    
1143 +  //NoCorrections
1144 +  if (EffectiveAreaTarget == kEleEANoCorr) {
1145 +    return 0.0;
1146 +  }
1147 +  //2011 Data Effective Areas
1148 +  else if (EffectiveAreaTarget == kEleEAData2011) {
1149 +    if (type == kEleGammaIsoDR0p0To0p1) {
1150 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1151 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
1152 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.005;
1153 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1154 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.004;
1155 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1156 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1157 +    }
1158 +    if (type == kEleGammaIsoDR0p1To0p2) {
1159 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1160 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1161 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.019;
1162 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1163 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
1164 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.035;
1165 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.041;
1166 +    }
1167 +    if (type == kEleGammaIsoDR0p2To0p3) {
1168 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1169 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1170 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.014;
1171 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1172 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.039;
1173 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.042;
1174 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.048;
1175 +    }
1176 +    if (type == kEleGammaIsoDR0p3To0p4) {
1177 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1178 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.029;
1179 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1180 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.029;
1181 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.042;
1182 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.047;
1183 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.054;
1184 +    }
1185 +    if (type == kEleGammaIsoDR0p4To0p5) {
1186 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1187 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.038;
1188 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.028;
1189 +     if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1190 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.047;
1191 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.057;
1192 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1193 +    }
1194 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1195 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1196 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.002;
1197 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.002;
1198 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1199 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1200 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1201 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1202 +    }
1203 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1204 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.005;
1205 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1206 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.008;
1207 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1208 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1209 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.001;
1210 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.003;
1211 +    }
1212 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1213 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1214 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1215 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1216 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1217 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1218 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1219 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.019;
1220 +    }
1221 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1222 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1223 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1224 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.025;
1225 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.030;
1226 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.036;
1227 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.038;
1228 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.084;
1229 +    }
1230 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1231 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.020;
1232 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1233 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.035;
1234 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1235 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1236 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.107;
1237 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.228;
1238 +    }
1239 +  }
1240 +
1241 +  //Summer11 MC Effective Areas
1242 +  else if (EffectiveAreaTarget == kEleEASummer11MC) {
1243 +    if (type == kEleGammaIsoDR0p0To0p1) {
1244 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.015;
1245 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1246 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1247 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1248 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1249 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.024;
1250 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.023;
1251 +    }
1252 +    if (type == kEleGammaIsoDR0p1To0p2) {
1253 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1254 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1255 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1256 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.037;
1257 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.046;
1258 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.055;
1259 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.046;
1260 +    }
1261 +    if (type == kEleGammaIsoDR0p2To0p3) {
1262 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.021;
1263 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.018;
1264 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1265 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.026;
1266 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1267 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.045;
1268 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.059;
1269 +    }
1270 +    if (type == kEleGammaIsoDR0p3To0p4) {
1271 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.036;
1272 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.030;
1273 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1274 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.036;
1275 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.058;
1276 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.073;
1277 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.083;
1278 +    }
1279 +    if (type == kEleGammaIsoDR0p4To0p5) {
1280 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.053;
1281 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.037;
1282 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.032;
1283 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.048;
1284 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.062;
1285 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.085;
1286 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.118;
1287 +    }
1288 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1289 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.000;
1290 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.000;
1291 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1292 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1293 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1294 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1295 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1296 +    }
1297 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1298 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.004;
1299 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.007;
1300 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1301 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.004;
1302 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.003;
1303 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1304 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.004;
1305 +    }
1306 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1307 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.008;
1308 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
1309 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1310 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.013;
1311 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.014;
1312 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.016;
1313 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1314 +    }
1315 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1316 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1317 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1318 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.020;
1319 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1320 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.040;
1321 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.036;
1322 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.086;
1323 +    }
1324 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1325 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.016;
1326 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.026;
1327 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.030;
1328 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1329 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1330 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.105;
1331 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.169;
1332 +    }
1333 +  }
1334 +  
1335 +  //Fall11 MC Effective Areas
1336 +  else if (EffectiveAreaTarget == kEleEAFall11MC) {
1337 +    if (type == kEleGammaIsoDR0p0To0p1) {
1338 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.014;
1339 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1340 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1341 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.012;
1342 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.016;
1343 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.021;
1344 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.012;
1345 +    }
1346 +    if (type == kEleGammaIsoDR0p1To0p2) {
1347 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.012;
1348 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1349 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.015;
1350 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.042;
1351 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.055;
1352 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.068;
1353 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.067;
1354 +    }
1355 +    if (type == kEleGammaIsoDR0p2To0p3) {
1356 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.024;
1357 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.020;
1358 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1359 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.038;
1360 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.051;
1361 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.066;
1362 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.080;
1363 +    }
1364 +    if (type == kEleGammaIsoDR0p3To0p4) {
1365 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.040;
1366 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1367 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1368 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.047;
1369 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1370 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.083;
1371 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.123;
1372 +    }
1373 +    if (type == kEleGammaIsoDR0p4To0p5) {
1374 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1375 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.041;
1376 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.037;
1377 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.057;
1378 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.095;
1379 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1380 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.133;
1381 +    }
1382 +    if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1383 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1384 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.003;
1385 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.000;
1386 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1387 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.000;
1388 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.000;
1389 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.000;
1390 +    }
1391 +    if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1392 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.006;
1393 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.008;
1394 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.010;
1395 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.006;
1396 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1397 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.002;
1398 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.007;
1399 +    }
1400 +    if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1401 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.009;
1402 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.014;
1403 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.018;
1404 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.016;
1405 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.017;
1406 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1407 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.021;
1408 +    }
1409 +    if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1410 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1411 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.019;
1412 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1413 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.035;
1414 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.037;
1415 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1416 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.110;
1417 +    }
1418 +    if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1419 +      if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1420 +      if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.027;
1421 +      if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.036;
1422 +      if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.045;
1423 +      if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.057;
1424 +      if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.123;
1425 +      if (fabs(SCEta) >= 2.4) EffectiveArea = 0.220;
1426 +    }
1427 +  }
1428 +
1429 +
1430 +
1431 +  return EffectiveArea;  
1432 + }
1433 +
1434 +

Diff Legend

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