ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc

Comparing UserCode/MitPhysics/Mods/src/PhotonTreeWriter.cc (file contents):
Revision 1.77 by bendavid, Fri Dec 13 01:06:15 2013 UTC vs.
Revision 1.78 by veverka, Fri Dec 13 02:58:46 2013 UTC

#	Line 2301 | Line 2301 | void PhotonTreeWriter::ApplyLeptonTag(co
2301		// electron tagged
2302		fDiphotonEvent->leptonTag = 1;
2303
2304	<	fDiphotonEvent-> elePt = fLeptonTagElectrons->At(0)->Pt();
2305	<	fDiphotonEvent-> eleEta = fLeptonTagElectrons->At(0)->Eta();
2306	<	fDiphotonEvent-> elePhi = fLeptonTagElectrons->At(0)->Phi();
2307	<	fDiphotonEvent-> eleSCEta = fLeptonTagElectrons->At(0)->SCluster()->Eta();
2308	<	fDiphotonEvent-> eleIso1 = (fLeptonTagElectrons->At(0)->TrackIsolationDr03() + fLeptonTagElectrons->At(0)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(0)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoRandomLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2304	>	fDiphotonEvent-> elePt = electron->Pt();
2305	>	fDiphotonEvent-> eleEta = electron->Eta();
2306	>	fDiphotonEvent-> elePhi = electron->Phi();
2307	>	fDiphotonEvent-> eleSCEta = electron->SCluster()->Eta();
2308	>	fDiphotonEvent-> eleIso1 = (electron->TrackIsolationDr03() + electron->EcalRecHitIsoDr03() + electron->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoRandomLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2309
2310		fDiphotonEvent-> eleIso2 = -99.;
2311
2312		if ( fDoSynching ) {
2313		Double_t distVtx = 999.0;
2314		for(UInt_t nv=0; nv<fPV->GetEntries(); nv++){
2315	<	double dz = TMath::Abs(fLeptonTagElectrons->At(0)->GsfTrk()->DzCorrected(*fPV->At(nv)));
2315	>	double dz = TMath::Abs(electron->GsfTrk()->DzCorrected(*fPV->At(nv)));
2316		if(dz < distVtx) {
2317		distVtx    = dz;
2318		closestVtx = nv;
2319		}
2320		}
2321	<	fDiphotonEvent-> eleIdMva = fElectronIDMVA->MVAValue(fLeptonTagElectrons->At(0), fPV->At(closestVtx));
2321	>	fDiphotonEvent-> eleIdMva = fElectronIDMVA->MVAValue(electron, fPV->At(closestVtx));
2322		}
2323
2324
2325	<	fDiphotonEvent-> eleIso3 = (fLeptonTagElectrons->At(0)->TrackIsolationDr03() + fLeptonTagElectrons->At(0)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(0)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2326	<	fDiphotonEvent-> eleIso4 = (fLeptonTagElectrons->At(0)->TrackIsolationDr03() + fLeptonTagElectrons->At(0)->EcalRecHitIsoDr03() + fLeptonTagElectrons->At(0)->HcalTowerSumEtDr03() - fPileUpDen->At(0)->Rho() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2327	<	fDiphotonEvent-> eleDist = fLeptonTagElectrons->At(0)->ConvPartnerDist();
2328	<	fDiphotonEvent-> eleDcot = fLeptonTagElectrons->At(0)->ConvPartnerDCotTheta();
2329	<	fDiphotonEvent-> eleCoviee = fLeptonTagElectrons->At(0)->CoviEtaiEta();
2330	<	fDiphotonEvent-> eleDphiin = TMath::Abs(fLeptonTagElectrons->At(0)->DeltaPhiSuperClusterTrackAtVtx());
2331	<	fDiphotonEvent-> eleDetain = TMath::Abs(fLeptonTagElectrons->At(0)->DeltaEtaSuperClusterTrackAtVtx());
2332	<	fDiphotonEvent-> eleDR1 = MathUtils::DeltaR(fLeptonTagElectrons->At(0),phHard);
2333	<	fDiphotonEvent-> eleDR2 = MathUtils::DeltaR(fLeptonTagElectrons->At(0),phSoft);
2334	<	fDiphotonEvent-> eleMass1 = (phHard->Mom()+fLeptonTagElectrons->At(0)->Mom()).M();
2335	<	fDiphotonEvent-> eleMass2 = (phSoft->Mom()+fLeptonTagElectrons->At(0)->Mom()).M();
2336	<	fDiphotonEvent-> eleNinnerHits =      fLeptonTagElectrons->At(0)->Trk()->NExpectedHitsInner();
2325	>	fDiphotonEvent-> eleIso3 = (electron->TrackIsolationDr03() + electron->EcalRecHitIsoDr03() + electron->HcalTowerSumEtDr03() - fPileUpDen->At(0)->RhoLowEta() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2326	>	fDiphotonEvent-> eleIso4 = (electron->TrackIsolationDr03() + electron->EcalRecHitIsoDr03() + electron->HcalTowerSumEtDr03() - fPileUpDen->At(0)->Rho() * TMath::Pi() * 0.3 * 0.3)/fDiphotonEvent-> elePt;
2327	>	fDiphotonEvent-> eleDist = electron->ConvPartnerDist();
2328	>	fDiphotonEvent-> eleDcot = electron->ConvPartnerDCotTheta();
2329	>	fDiphotonEvent-> eleCoviee = electron->CoviEtaiEta();
2330	>	fDiphotonEvent-> eleDphiin = TMath::Abs(electron->DeltaPhiSuperClusterTrackAtVtx());
2331	>	fDiphotonEvent-> eleDetain = TMath::Abs(electron->DeltaEtaSuperClusterTrackAtVtx());
2332	>	fDiphotonEvent-> eleDR1 = MathUtils::DeltaR(electron,phHard);
2333	>	fDiphotonEvent-> eleDR2 = MathUtils::DeltaR(electron,phSoft);
2334	>	fDiphotonEvent-> eleMass1 = (phHard->Mom()+electron->Mom()).M();
2335	>	fDiphotonEvent-> eleMass2 = (phSoft->Mom()+electron->Mom()).M();
2336	>	fDiphotonEvent-> eleNinnerHits =      electron->Trk()->NExpectedHitsInner();
2337		} // electron tagged
2338
2339		if(false){
#	Line 2547 | Line 2547 | void PhotonTreeWriter::ApplyVHLepTag(con
2547		!MassOfPairIsWithinWindowAroundMZ(phSoft, electron, 10)) {
2548		isVHLepLoose = true;   // low MET event
2549		} // Low MET event.
2550	+	fDiphotonEvent->elePt  = electron->Pt ();
2551	+	fDiphotonEvent->eleEta = electron->Eta();
2552	+	fDiphotonEvent->elePhi = electron->Phi();
2553		} // Found a good VH(lep) tag electron.
2554
2555		//printf("check electron, tight  = %i, loose = %i\n",int(isVHLepTight),int(isVHLepLoose));

Diff Legend

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