From ae7602a900099e228d76cdb21ad81901d6ddcad0 Mon Sep 17 00:00:00 2001 From: Salvatore Rappoccio Date: Tue, 3 Dec 2013 11:52:44 -0600 Subject: [PATCH] Forgot to commit this : added pt of the generator-level top --- IHeartNY/test/iheartny_topxs_fwlite.py | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) diff --git a/IHeartNY/test/iheartny_topxs_fwlite.py b/IHeartNY/test/iheartny_topxs_fwlite.py index b831de4f..f6c9f2bd 100644 --- 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 :