) args = parser.parse_args() cuts_jet_tag = [ ("%dj%dt"%(n,m), Cuts.n_jets(n)*Cuts.n_tags(m)) for n in [2,3] for m in [0,1,2] ] # cuts_jet_tag = [ # ("2j1t", Cuts.n_jets(2)*Cuts.n_tags(1)) # ] cuts = [] for cutname, cbline in cuts_jet_tag: for lep in ["mu", "ele"]: cn = "%s_%s" % (lep, cutname) baseline = Cuts.lepton(lep) * Cuts.hlt(lep) * Cuts.metmt(lep) * Cuts.rms_lj cuts += [ #Without the MET cut ("%s_nomet" % cn, Cuts.lepton(lep) * Cuts.hlt(lep) * Cuts.rms_lj * cbline), #Baseline for fit ("%s_baseline" % cn, baseline * cbline), #Cut-based check ("%s_cutbased_final" % cn, baseline * cbline * Cuts.top_mass_sig * Cuts.eta_lj ), #MVA-based selection ("%s_mva_loose" % cn, baseline * cbline * Cuts.mva_wp(lep)
for fi in sys.argv[1:]: print fi samp = Sample.fromFile(fi) print "Lumi scale factor: ", samp.lumiScaleFactor(20000) for lep in ["mu", "ele"]: if "/%s/"%lep not in fi: continue print lep cut = None for cutname, _cut in [ ("hlt", Cuts.hlt(lep)), ("lep", Cuts.single_lepton(lep)), ("2J", Cuts.n_jets(2)), ("1T", Cuts.n_tags(1)), ("MET/MtW", Cuts.metmt(lep)), ("rms", Cuts.rms_lj), ("Mtop", Cuts.top_mass_sig), ("etalj", Cuts.eta_lj) ]: if not cut: cut = _cut else: cut *= _cut try: hi = samp.drawHistogram("eta_lj", str(cut), binning=[50, -5, 5]) hi.Scale(samp.lumiScaleFactor(20000)) print cutname print hi.GetEntries(), hi.Integral() except: print "-1 -1"