ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Utils/src/ElectronTools.cc

Comparing UserCode/MitPhysics/Utils/src/ElectronTools.cc (file contents):
Revision 1.33 by sixie, Sat Dec 31 23:17:40 2011 UTC vs.
Revision 1.37 by ceballos, Fri Apr 20 20:25:58 2012 UTC

#	Line 189 | 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 ||
#	Line 1059 | Line 1069 | Double_t ElectronTools::Likelihood(Elect
1069
1070		}
1071
1072	<	Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t Eta) {
1072	>	Double_t ElectronTools::ElectronEffectiveArea(EElectronEffectiveAreaType type, Double_t SCEta, EElectronEffectiveAreaTarget EffectiveAreaTarget) {
1073
1074		Double_t EffectiveArea = 0;
1075	<
1076	<	if (fabs(Eta) < 1.0) {
1077	<	if (type == kEleChargedIso03) EffectiveArea = 0.000;
1078	<	if (type == kEleNeutralHadronIso03) EffectiveArea = 0.017;
1079	<	if (type == kEleGammaIso03) EffectiveArea = 0.045;
1080	<	if (type == kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.014;
1081	<	if (type == kEleChargedIso04) EffectiveArea = 0.000;
1082	<	if (type == kEleNeutralHadronIso04) EffectiveArea = 0.034;
1083	<	if (type == kEleGammaIso04) EffectiveArea = 0.079;
1084	<	if (type == kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.014;
1085	<	if (type == kEleNeutralHadronIso007) EffectiveArea = 0.000;
1086	<	if (type == kEleHoverE) EffectiveArea = 0.00016;
1087	<	if (type == kEleHcalDepth1OverEcal) EffectiveArea = 0.00016;
1088	<	if (type == kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
1089	<	} else if (fabs(Eta) >= 1.0 && fabs(Eta) < 1.479 ) {
1090	<	if (type == kEleChargedIso03) EffectiveArea = 0.000;
1091	<	if (type == kEleNeutralHadronIso03) EffectiveArea = 0.025;
1092	<	if (type == kEleGammaIso03) EffectiveArea = 0.052;
1093	<	if (type == kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.030;
1094	<	if (type == kEleChargedIso04) EffectiveArea = 0.000;
1095	<	if (type == kEleNeutralHadronIso04) EffectiveArea = 0.050;
1096	<	if (type == kEleGammaIso04) EffectiveArea = 0.073;
1097	<	if (type == kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.030;
1098	<	if (type == kEleNeutralHadronIso007) EffectiveArea = 0.000;
1099	<	if (type == kEleHoverE) EffectiveArea = 0.00022;
1100	<	if (type == kEleHcalDepth1OverEcal) EffectiveArea = 0.00022;
1101	<	if (type == kEleHcalDepth2OverEcal) EffectiveArea = 0.00000;
1102	<	} else if (fabs(Eta) >= 1.479 && fabs(Eta) < 2.0 ) {
1103	<	if (type == kEleChargedIso03) EffectiveArea = 0.000;
1104	<	if (type == kEleNeutralHadronIso03) EffectiveArea = 0.030;
1105	<	if (type == kEleGammaIso03) EffectiveArea = 0.170;
1106	<	if (type == kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.134;
1107	<	if (type == kEleChargedIso04) EffectiveArea = 0.000;
1108	<	if (type == kEleNeutralHadronIso04) EffectiveArea = 0.060;
1109	<	if (type == kEleGammaIso04) EffectiveArea = 0.187;
1110	<	if (type == kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.134;
1111	<	if (type == kEleNeutralHadronIso007) EffectiveArea = 0.000;
1112	<	if (type == kEleHoverE) EffectiveArea = 0.00030;
1113	<	if (type == kEleHcalDepth1OverEcal) EffectiveArea = 0.00026;
1114	<	if (type == kEleHcalDepth2OverEcal) EffectiveArea = 0.00002;
1115	<	} else if (fabs(Eta) >= 2.0 && fabs(Eta) < 2.25 ) {
1116	<	if (type == kEleChargedIso03) EffectiveArea = 0.000;
1117	<	if (type == kEleNeutralHadronIso03) EffectiveArea = 0.022;
1118	<	if (type == kEleGammaIso03) EffectiveArea = 0.623;
1119	<	if (type == kEleGammaIsoVetoEtaStrip03) EffectiveArea = 0.516;
1120	<	if (type == kEleChargedIso04) EffectiveArea = 0.000;
1121	<	if (type == kEleNeutralHadronIso04) EffectiveArea = 0.055;
1122	<	if (type == kEleGammaIso04) EffectiveArea = 0.659;
1123	<	if (type == kEleGammaIsoVetoEtaStrip04) EffectiveArea = 0.517;
1124	<	if (type == kEleNeutralHadronIso007) EffectiveArea = 0.000;
1125	<	if (type == kEleHoverE) EffectiveArea = 0.00054;
1126	<	if (type == kEleHcalDepth1OverEcal) EffectiveArea = 0.00045;
1127	<	if (type == kEleHcalDepth2OverEcal) EffectiveArea = 0.00003;
1128	<	} else if (fabs(Eta) >= 2.25 && fabs(Eta) < 2.5 ) {
1129	<	if (type == kEleChargedIso03) EffectiveArea = 0.000;
1130	<	if (type == kEleNeutralHadronIso03) EffectiveArea = 0.018;
1131	<	if (type == kEleGammaIso03) EffectiveArea = 1.198;
1132	<	if (type == kEleGammaIsoVetoEtaStrip03) EffectiveArea = 1.049;
1133	<	if (type == kEleChargedIso04) EffectiveArea = 0.000;
1134	<	if (type == kEleNeutralHadronIso04) EffectiveArea = 0.073;
1135	<	if (type == kEleGammaIso04) EffectiveArea = 1.258;
1136	<	if (type == kEleGammaIsoVetoEtaStrip04) EffectiveArea = 1.051;
1137	<	if (type == kEleNeutralHadronIso007) EffectiveArea = 0.000;
1138	<	if (type == kEleHoverE) EffectiveArea = 0.00082;
1139	<	if (type == kEleHcalDepth1OverEcal) EffectiveArea = 0.00066;
1140	<	if (type == kEleHcalDepth2OverEcal) EffectiveArea = 0.00004;
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

Diff Legend

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