Forgot to commit this : added pt of the generator-level top

rappoccio · Dec 3, 2013 · ae7602a · ae7602a
1 parent c7d1ba3
commit ae7602a
Showing 1 changed file with 8 additions and 3 deletions.
diff --git a/IHeartNY/test/iheartny_topxs_fwlite.py b/IHeartNY/test/iheartny_topxs_fwlite.py
@@ -312,7 +312,9 @@ def match( self, jets ) :
 if options.makeResponse == True : 
     response = ROOT.RooUnfoldResponse(10, 300., 1300., 10, 300., 1300.)
     response.SetName('response_pt')
-    ptRecoTop = ROOT.TH1F("ptRecoTop", "Reconstructed top p_{T};p_{T} (GeV/c);Number", 10, 300., 1300.)
+    ptGenTop = ROOT.TH1F("ptGenTop", "Generated top p_{T};p_{T} (GeV/c);Number", 10, 300., 1300.)
+
+ptRecoTop = ROOT.TH1F("ptRecoTop", "Reconstructed top p_{T};p_{T} (GeV/c);Number", 10, 300., 1300.)
 
 events = Events (files)
 
@@ -593,6 +595,9 @@ def match( self, jets ) :
         else :
             hadTop = topQuarks[1]
             lepTop = topQuarks[0]
+	ptGenTop.Fill( hadTop.p4.Perp(), weight )
+    # endif (making response matrix)
+
 
 
     lepType = 0 # Let 0 = muon, 1 = electron    
@@ -1075,7 +1080,7 @@ def match( self, jets ) :
 
             passSelection = False
             if (jet.DeltaR( lepP4) > ROOT.TMath.Pi() / 2.0) :
-                if (topTagPt[ijet] > 250.):
+                if (topTagPt[ijet] > 200.):
 		    topTagptHistprecuts.Fill(topTagPt[ijet],t1weight)
 		    htLep3t1kin.Fill(htLepVal,t1weight)
 		    nsj.Fill(topTagNSub[ijet],t1weight)
@@ -1151,7 +1156,7 @@ def match( self, jets ) :
 							topTagMassHistPostBDMax.Fill(topTagMass[ijet],t1weight)
 	    if passSelection == True :
 		ptRecoTop.Fill( topTagPt[ijet], t1weight )
-            if options.makeResponse == True : 
+            if options.makeResponse == True :		
                 if passSelection == True :
                     response.Fill(  topTagPt[ijet], hadTop.p4.Perp(), t1weight )
                 else :