ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/MitPhysics/Utils/src/RecoilTools.cc (file contents):
Revision 1.3 by pharris, Sat Apr 7 09:36:33 2012 UTC vs.
Revision 1.4 by pharris, Fri Apr 13 14:24:00 2012 UTC

#	Line 80 | Line 80 | void RecoilTools::addNeut(const PFJet *i
80		if(iSign > 0) iVec -= lVec;
81		if(iSign < 0) iVec += lVec;
82		iSumEt += lPt;
83	+	//=== Above was a bug in the training
84	+	//if(iSign > 0) iSumEt += lPt;
85	+	//if(iSign < 0) iSumEt -= lPt;
86		}
87
88		//--------------------------------------------------------------------------------------------------
#	Line 89 | Line 92 | void RecoilTools::addNeut(const PFJet *i
92		double lPt = iJet->Pt();
93		lPt *= (iJet->NeutralEmEnergy()/iJet->E() + iJet->NeutralHadronEnergy()/iJet->E());
94		lVec.SetPt(lPt); lVec.SetEta(iJet->Eta()); lVec.SetPhi(iJet->Phi()); lVec.SetM(iJet->Mass());
95	<	if(iSign > 0) iVec -= lVec;
96	<	if(iSign < 0) iVec += lVec;
95	>	if(iSign > 0) iVec   -= lVec;
96	>	if(iSign < 0) iVec   += lVec;
97		iSumEt += lPt;
98	+	//=== Above was a bug in the training
99	+	//if(iSign > 0) iSumEt += lPt;
100	+	//if(iSign < 0) iSumEt -= lPt;
101		}
102
103		//--------------------------------------------------------------------------------------------------
#	Line 111 | Line 117 | Met RecoilTools::NoPUMet( const PFJetCol
117		lVec     -= pPF->Mom();
118		lSumEt   += pPF->Pt();
119		}
120	+	int lNPass = 0;
121		for(UInt_t i0 = 0; i0 < iJets->GetEntries(); i0++) {
122		const PFJet *pJet = iJets->At(i0);
116	–	if(!fJetIDMVA->pass(pJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity)) continue;
123		if(!filter(pJet,iPhi1,iEta1,iPhi2,iEta2))                                       continue; //Quick cleaning==> if not done already
124	+	if(!fJetIDMVA->pass(pJet,iVertex,iVertices,iJetCorrector,iPileupEnergyDensity)) continue;
125		addNeut(pJet,lVec,lSumEt,iJetCorrector,iPileupEnergyDensity);
126	+	lNPass++;
127		}
128		Met lMet(lVec.Px(),lVec.Py());
129	<	lMet.SetSumEt(lSumEt);
129	>	lMet.SetSumEt( lSumEt);
130		return lMet;
131		}
132		//--------------------------------------------------------------------------------------------------

Diff Legend

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