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.35 by sixie, Sun Apr 15 12:06:26 2012 UTC vs.
Revision 1.52 by pharris, Sat Jan 12 11:49:50 2013 UTC

#	Line 189 | Line 189 | Bool_t ElectronTools::PassCustomID(const
189		if(eSeedOverPin<fCuts[4][cat+4*eb])
190		return kFALSE;
191
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 ||
#	Line 344 | Line 355 | Bool_t ElectronTools::PassConversionFilt
355		}
356
357		//--------------------------------------------------------------------------------------------------
358	+	Bool_t ElectronTools::PassConversionFilterPFAOD(const Electron *ele,
359	+	const DecayParticleCol *conversions,
360	+	const BaseVertex *vtx,
361	+	UInt_t nWrongHitsMax,
362	+	Double_t probMin,
363	+	Double_t lxyMin,
364	+	Bool_t matchCkf,
365	+	Bool_t requireArbitratedMerged,
366	+	Double_t trkptMin)
367	+	{
368	+
369	+	Bool_t isGoodConversion = kFALSE;
370	+
371	+	for (UInt_t ifc=0; ifc<conversions->GetEntries(); ifc++) {
372	+	Bool_t ConversionMatchFound = kFALSE;
373	+	for (UInt_t d=0; d<conversions->At(ifc)->NDaughters(); d++) {
374	+	const ChargedParticle *pParticle = 0;
375	+	pParticle = dynamic_cast<const ChargedParticle*>(conversions->At(ifc)->Daughter(d));
376	+	if(pParticle == 0) continue;
377	+	const Track* trk = 0;
378	+	trk = pParticle->Trk();
379	+	if(trk == 0) continue;
380	+	if (ele->GsfTrk() == trk || (matchCkf && ele->TrackerTrk()==trk) ) {
381	+	ConversionMatchFound = kTRUE;
382	+	break;
383	+	}
384	+	}
385	+
386	+	// if match between the e-track and one of the conversion legs
387	+	if (ConversionMatchFound == kTRUE){
388	+	isGoodConversion =  (conversions->At(ifc)->Prob() > probMin) &&
389	+	(!requireArbitratedMerged || conversions->At(ifc)->Quality().Quality(ConversionQuality::arbitratedMerged)) &&
390	+	(conversions->At(ifc)->LxyCorrected(vtx) > lxyMin);
391	+
392	+	if (isGoodConversion == kTRUE) {
393	+	for (UInt_t d=0; d<conversions->At(ifc)->NDaughters(); d++) {
394	+	const ChargedParticle *pParticle = 0;
395	+	pParticle = dynamic_cast<const ChargedParticle*>(conversions->At(ifc)->Daughter(d));
396	+	if(pParticle == 0) continue;
397	+	const Track* trk = 0;
398	+	trk = pParticle->Trk();
399	+	if(trk == 0) continue;
400	+	if (trk) {
401	+	if (trk->Pt()<trkptMin) isGoodConversion = kFALSE;
402	+	const StableData *sd = dynamic_cast<const StableData*>
403	+	(conversions->At(ifc)->DaughterDat(d));
404	+	if (sd->NWrongHits() > nWrongHitsMax)
405	+	isGoodConversion = kFALSE;
406	+	} else {
407	+	isGoodConversion = kFALSE;
408	+	}
409	+	}
410	+	}
411	+	}
412	+
413	+	if (isGoodConversion == kTRUE) break;
414	+
415	+	} // loop over all conversions
416	+
417	+	return !isGoodConversion;
418	+	}
419	+	//--------------------------------------------------------------------------------------------------
420		Bool_t ElectronTools::PassD0Cut(const Electron *ele, const VertexCol *vertices, Double_t fD0Cut, Int_t nVertex)
421		{
422		Bool_t d0cut = kFALSE;
423
424		Double_t d0_real = 1e30;
425	+
426	+	if( nVertex >= (int) vertices->GetEntries() )
427	+	nVertex = vertices->GetEntries() - 1;
428	+
429		if(nVertex >= 0) d0_real = TMath::Abs(ele->GsfTrk()->D0Corrected(*vertices->At(nVertex)));
430		else            {
431		Double_t distVtx = 999.0;
#	Line 388 | Line 465 | Bool_t ElectronTools::PassDZCut(const El
465		Bool_t dzcut = kFALSE;
466
467		Double_t distVtx = 999.0;
468	+
469	+	if( nVertex >= (int) vertices->GetEntries() )
470	+	nVertex = vertices->GetEntries()-1;
471	+
472		if(nVertex >= 0) distVtx = TMath::Abs(ele->GsfTrk()->DzCorrected(*vertices->At(nVertex)));
473		else {
474		for(UInt_t nv=0; nv<vertices->GetEntries(); nv++){
#	Line 397 | Line 478 | Bool_t ElectronTools::PassDZCut(const El
478		}
479		}
480		}
481	<
481	>
482		if(distVtx < fDZCut) dzcut = kTRUE;
483
484		return dzcut;
#	Line 1063 | Line 1144 | Double_t ElectronTools::ElectronEffectiv
1144
1145		Double_t EffectiveArea = 0;
1146
1147	+	if(type == ElectronTools::kEleNeutralIso04){
1148	+	if     (fabs(SCEta) < 1.0                           ) EffectiveArea = 0.19;
1149	+	else if(fabs(SCEta) >= 1.0   && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1150	+	else if(fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.000 ) EffectiveArea = 0.12;
1151	+	else if(fabs(SCEta) >= 2.000 && fabs(SCEta) < 2.200 ) EffectiveArea = 0.21;
1152	+	else if(fabs(SCEta) >= 2.200 && fabs(SCEta) < 2.300 ) EffectiveArea = 0.27;
1153	+	else if(fabs(SCEta) >= 2.300 && fabs(SCEta) < 2.400 ) EffectiveArea = 0.44;
1154	+	else                                                  EffectiveArea = 0.52;
1155	+	return EffectiveArea;
1156	+	}
1157	+
1158		if (fabs(SCEta) < 1.0) {
1159		if (type == ElectronTools::kEleChargedIso03) EffectiveArea = 0.000;
1160		if (type == ElectronTools::kEleNeutralHadronIso03) EffectiveArea = 0.017;
#	Line 1134 | Line 1226 | Double_t ElectronTools::ElectronEffectiv
1226		if (EffectiveAreaTarget == kEleEANoCorr) {
1227		return 0.0;
1228		}
1229	+	else if (EffectiveAreaTarget == kEleEAData2012) {
1230	+	if (type == kEleGammaIso03) {
1231	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.122;
1232	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.147;
1233	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.055;
1234	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.106;
1235	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.138;
1236	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.221;
1237	+	if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.211;
1238	+	}
1239	+	if (type == kEleGammaIso04) {
1240	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.176;
1241	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.206;
1242	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.094;
1243	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.172;
1244	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.244;
1245	+	if (fabs(SCEta) >= 2.4 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.333;
1246	+	if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.348;
1247	+	}
1248	+	if (type == kEleNeutralHadronIso03) {
1249	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1250	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.021;
1251	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.013;
1252	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.010;
1253	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1254	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.020;
1255	+	if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.019;
1256	+	}
1257	+	if (type == kEleNeutralHadronIso04) {
1258	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.022;
1259	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.036;
1260	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.027;
1261	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
1262	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.052;
1263	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.063;
1264	+	if (fabs(SCEta) >= 2.4 ) EffectiveArea = 0.028;
1265	+	}
1266	+	if (type == kEleGammaAndNeutralHadronIso04) {
1267	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.19;
1268	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.25;
1269	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.12;
1270	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.21;
1271	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.27;
1272	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.44;
1273	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.52;
1274	+	}
1275	+	if (type == kEleGammaIsoDR0p0To0p1) {
1276	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.051;
1277	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1278	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.006;
1279	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.007;
1280	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.024;
1281	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.013;
1282	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.013;
1283	+	}
1284	+	if (type == kEleGammaIsoDR0p1To0p2) {
1285	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.013;
1286	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.013;
1287	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.021;
1288	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.052;
1289	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.066;
1290	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.043;
1291	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.102;
1292	+	}
1293	+	if (type == kEleGammaIsoDR0p2To0p3) {
1294	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.026;
1295	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.017;
1296	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
1297	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.028;
1298	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.041;
1299	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.034;
1300	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.042;
1301	+	}
1302	+	if (type == kEleGammaIsoDR0p3To0p4) {
1303	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.039;
1304	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.032;
1305	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.017;
1306	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.024;
1307	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.053;
1308	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.059;
1309	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.069;
1310	+	}
1311	+	if (type == kEleGammaIsoDR0p4To0p5) {
1312	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.059;
1313	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.045;
1314	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.033;
1315	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.043;
1316	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.056;
1317	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.065;
1318	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.074;
1319	+	}
1320	+	if (type == kEleNeutralHadronIsoDR0p0To0p1) {
1321	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.001;
1322	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.006;
1323	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.001;
1324	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.000;
1325	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.001;
1326	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.003;
1327	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.008;
1328	+	}
1329	+	if (type == kEleNeutralHadronIsoDR0p1To0p2) {
1330	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.002;
1331	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.001;
1332	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.004;
1333	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.003;
1334	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.005;
1335	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.006;
1336	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.010;
1337	+	}
1338	+	if (type == kEleNeutralHadronIsoDR0p2To0p3) {
1339	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.007;
1340	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.010;
1341	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.009;
1342	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.009;
1343	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.007;
1344	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.018;
1345	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.028;
1346	+	}
1347	+	if (type == kEleNeutralHadronIsoDR0p3To0p4) {
1348	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.010;
1349	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.011;
1350	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.012;
1351	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.008;
1352	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.018;
1353	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.026;
1354	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.063;
1355	+	}
1356	+	if (type == kEleNeutralHadronIsoDR0p4To0p5) {
1357	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.011;
1358	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.012;
1359	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.016;
1360	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.023;
1361	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.038;
1362	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.051;
1363	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.143;
1364	+	}
1365	+	}
1366	+
1367		//2011 Data Effective Areas
1368		else if (EffectiveAreaTarget == kEleEAData2011) {
1369	+	if (type == kEleGammaAndNeutralHadronIso03) {
1370	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.100;
1371	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.120;
1372	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.085;
1373	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.110;
1374	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.120;
1375	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.120;
1376	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.130;
1377	+	}
1378	+	if (type == kEleGammaAndNeutralHadronIso04) {
1379	+	if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.180;
1380	+	if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.200;
1381	+	if (fabs(SCEta) >= 1.479 && fabs(SCEta) < 2.0 ) EffectiveArea = 0.150;
1382	+	if (fabs(SCEta) >= 2.0 && fabs(SCEta) < 2.2 ) EffectiveArea = 0.190;
1383	+	if (fabs(SCEta) >= 2.2 && fabs(SCEta) < 2.3 ) EffectiveArea = 0.210;
1384	+	if (fabs(SCEta) >= 2.3 && fabs(SCEta) < 2.4 ) EffectiveArea = 0.220;
1385	+	if (fabs(SCEta) >= 2.4) EffectiveArea = 0.290;
1386	+	}
1387		if (type == kEleGammaIsoDR0p0To0p1) {
1388		if (fabs(SCEta) >= 0.0 && fabs(SCEta) < 1.0 ) EffectiveArea = 0.017;
1389		if (fabs(SCEta) >= 1.0 && fabs(SCEta) < 1.479 ) EffectiveArea = 0.033;
#	Line 1416 | Line 1664 | Double_t ElectronTools::ElectronEffectiv
1664		}
1665		}
1666
1667	+	return EffectiveArea;
1668	+	}
1669
1670
1671	<	return EffectiveArea;
1671	>	Bool_t ElectronTools::PassHggLeptonTagID(const Electron* ele) {
1672	>
1673	>	float dist = ( ele->ConvPartnerDist()      == -9999.? 9999:TMath::Abs(ele->ConvPartnerDist()));
1674	>	float dcot = ( ele->ConvPartnerDCotTheta() == -9999.? 9999:TMath::Abs(ele->ConvPartnerDCotTheta()));
1675	>
1676	>	if (dist < 0.02) return false;
1677	>	if (dcot < 0.02) return false;
1678	>
1679	>	int numInnerHits = ele->Trk()->NExpectedHitsInner();
1680	>	if( numInnerHits > 1 ) return false;
1681	>
1682	>	float coviEtaiEta = ele->CoviEtaiEta();
1683	>	if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1684	>	else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.031 ) return false;
1685	>
1686	>	Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1687	>	Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1688	>
1689	>	if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.039 || deltaEtaIn > 0.005 ) ) return false;
1690	>	else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.028 || deltaEtaIn > 0.007 ) ) return false;
1691	>
1692	>	return true;
1693		}
1694
1695	+	Bool_t ElectronTools::PassHggLeptonTagID2012(const Electron* ele) {
1696	+
1697	+	if (TMath::Abs(1./ele->E()-1./ele->Pt())>0.05) return false;
1698	+
1699	+	int numInnerHits = ele->Trk()->NExpectedHitsInner();
1700	+	if( numInnerHits > 1 ) return false;
1701	+
1702	+	float coviEtaiEta = ele->CoviEtaiEta();
1703	+	if( ele->SCluster()->AbsEta() < 1.5 && coviEtaiEta > 0.01 ) return false;
1704	+	else if( ele->SCluster()->AbsEta() > 1.5 && coviEtaiEta > 0.03 ) return false; // h
1705	+
1706	+	Double_t deltaPhiIn   = TMath::Abs(ele->DeltaPhiSuperClusterTrackAtVtx());
1707	+	Double_t deltaEtaIn   = TMath::Abs(ele->DeltaEtaSuperClusterTrackAtVtx());
1708	+
1709	+	if( ele->SCluster()->AbsEta() < 1.5 && ( deltaPhiIn > 0.15 || deltaEtaIn > 0.007 || ele->HadronicOverEm()>0.12) ) return false;   // h
1710	+	else if ( ele->SCluster()->AbsEta() > 1.5 && ( deltaPhiIn > 0.10 || deltaEtaIn > 0.009 || ele->HadronicOverEm()>0.10 ) ) return false;   // h
1711
1712	+	return true;
1713	+	}
1714	+
1715	+	std::pair<Double_t,Double_t> ElectronTools::ComputeEPCombination( const Electron * ele, const float regression_energy,
1716	+	const float regression_energy_error) {
1717	+
1718	+	enum Classification { GSF_ELECTRON_UNKNOWN=-1,
1719	+	GSF_ELECTRON_GOLDEN=0,
1720	+	GSF_ELECTRON_BIGBREM=1,
1721	+	GSF_ELECTRON_BADTRACK=2,
1722	+	GSF_ELECTRON_SHOWERING=3,
1723	+	GSF_ELECTRON_GAP=4 } ;
1724	+
1725	+	float newEnergyError_ = regression_energy_error;
1726	+	float scEnergy = regression_energy;
1727	+
1728	+	int elClass = ele->Classification();
1729	+
1730	+	float trackMomentum  = ele->PIn() ;
1731	+	float errorTrackMomentum_ = 999. ;
1732	+
1733	+	// the electron's track momentum error was not available in bambu versions less than 029
1734	+	if(ele->TrackMomentumError() > 0)
1735	+	errorTrackMomentum_ = ele->TrackMomentumError();
1736	+	else if ( ele->GsfTrk()->PtErr() > 0)
1737	+	errorTrackMomentum_ = ele->GsfTrk()->PtErr()*cosh(ele->GsfTrk()->Eta());
1738	+	else
1739	+	assert(0);
1740	+
1741	+	float finalMomentum = ele->E(); // initial
1742	+	float finalMomentumError = 999.;
1743	+
1744	+	// first check for large errors
1745	+
1746	+	if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy <= 0.5) {
1747	+	finalMomentum = regression_energy;    finalMomentumError = regression_energy_error;
1748	+	}
1749	+	else if (errorTrackMomentum_/trackMomentum <= 0.5 && regression_energy_error/regression_energy > 0.5){
1750	+	finalMomentum = trackMomentum;  finalMomentumError = errorTrackMomentum_;
1751	+	}
1752	+	else if (errorTrackMomentum_/trackMomentum > 0.5 && regression_energy_error/regression_energy > 0.5){
1753	+	if (errorTrackMomentum_/trackMomentum < regression_energy_error/regression_energy) {
1754	+	finalMomentum = trackMomentum; finalMomentumError = errorTrackMomentum_;
1755	+	}
1756	+	else{
1757	+	finalMomentum = regression_energy; finalMomentumError = regression_energy_error;
1758	+	}
1759	+	}
1760	+	// then apply the combination algorithm
1761	+	else {
1762	+	// calculate E/p and corresponding error
1763	+	float eOverP = regression_energy / trackMomentum;
1764	+	float errorEOverP = sqrt(
1765	+	(regression_energy_error/trackMomentum)*(regression_energy_error/trackMomentum) +
1766	+	(regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum)*
1767	+	(regression_energy*errorTrackMomentum_/trackMomentum/trackMomentum));
1768	+
1769	+
1770	+	bool eleIsNotInCombination = false ;
1771	+	if ( (eOverP  > 1 + 2.5*errorEOverP) || (eOverP  < 1 - 2.5*errorEOverP) || (eOverP < 0.8) || (eOverP > 1.3) )
1772	+	{ eleIsNotInCombination = true ; }
1773	+	if (eleIsNotInCombination)
1774	+	{
1775	+	if (eOverP > 1)
1776	+	{ finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
1777	+	else
1778	+	{
1779	+	if (elClass == GSF_ELECTRON_GOLDEN)
1780	+	{ finalMomentum = regression_energy; finalMomentumError = regression_energy_error; }
1781	+	if (elClass == GSF_ELECTRON_BIGBREM)
1782	+	{
1783	+	if (regression_energy<36)
1784	+	{ finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
1785	+	else
1786	+	{ finalMomentum = regression_energy ; finalMomentumError = regression_energy_error ; }
1787	+	}
1788	+	if (elClass == GSF_ELECTRON_BADTRACK)
1789	+	{ finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
1790	+	if (elClass == GSF_ELECTRON_SHOWERING)
1791	+	{
1792	+	if (regression_energy<30)
1793	+	{ finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_; }
1794	+	else
1795	+	{ finalMomentum = regression_energy; finalMomentumError = regression_energy_error;}
1796	+	}
1797	+	if (elClass == GSF_ELECTRON_GAP)
1798	+	{
1799	+	if (regression_energy<60)
1800	+	{ finalMomentum = trackMomentum ; finalMomentumError = errorTrackMomentum_ ; }
1801	+	else
1802	+	{ finalMomentum = regression_energy; finalMomentumError = regression_energy_error ; }
1803	+	}
1804	+	}
1805	+	}
1806	+	else
1807	+	{
1808	+	// combination
1809	+	finalMomentum = (regression_energy/regression_energy_error/regression_energy_error + trackMomentum/errorTrackMomentum_/errorTrackMomentum_) /
1810	+	(1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
1811	+	float finalMomentumVariance = 1 / (1/regression_energy_error/regression_energy_error + 1/errorTrackMomentum_/errorTrackMomentum_);
1812	+	finalMomentumError = sqrt(finalMomentumVariance);
1813	+	}
1814	+	}
1815	+
1816	+	std::pair<Double_t,Double_t> value;
1817	+	value.first = finalMomentum;
1818	+	value.second = finalMomentumError;
1819	+	return value;
1820	+
1821	+	}

Diff Legend

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