from project_histos import Cuts, filter_alnum, Sample from make_plots import plot_hists, canvas_margin, legend import argparse if __name__=="__main__": parser = argparse.ArgumentParser() args = parser.parse_args() #fname = "/Users/joosep/Documents/stpol/data/step3_trees_Mar28/T_t.root" fnames = [ "/Users/joosep/Documents/stpol/data/T_t.root", "/Users/joosep/Documents/stpol/data/step3_trees_Apr04/iso/mc/T_t.root", "/Users/joosep/Documents/stpol/data/step3_trees_Apr05_noPuClean/iso/mc/T_t.root", ] samples = [Sample.fromFile(f) for f in fnames] sample_names = ["Apr16 (fixed norm)", "Apr04", "Apr05 (no PU cl.)"] samps = zip(samples, sample_names) cut = Cuts.mt_mu*Cuts.n_jets(2)*Cuts.n_tags(1)*Cuts.eta_lj*Cuts.top_mass_sig hists = [] for sample, name in samps: hist = sample.drawHistogram("cos_theta", cut.cut_str, plot_range=[40, -1, 1]) #normalize to 20/fb hist.normalize_lumi(20000) hist.hist.SetName(name) hist.hist.SetTitle("%s: %.3E" % (name, hist.hist.Integral())) hist.update()
from project_histos import Cuts, filter_alnum, Sample from make_plots import plot_hists, canvas_margin, legend import argparse from collections import OrderedDict as dict if __name__=="__main__": parser = argparse.ArgumentParser() parser.add_argument('-i', '--infile', type=str, required=True) parser.add_argument('--normalize_lumi', action="store_true", default=True) args = parser.parse_args() #fname = "/Users/joosep/Documents/stpol/data/step3_trees_Mar28/T_t.root" sample = Sample.fromFile(args.infile) # ev_ids = dict() # for ev in sample.tree: # evid = int(ev.event_id) # if evid not in ev_ids.keys(): # ev_ids[evid] = 1 # else: # ev_ids[evid] += 1 cuts = dict() cuts["MT_mu"] = Cuts.mt_mu cuts["2J"] = Cuts.n_jets(2) cuts["1T"] = Cuts.n_tags(1) #cuts["eta_lj"] = Cuts.eta_lj cuts["rms_lj"] = Cuts.rms_lj cuts["eta_jet"] = Cuts.eta_jet cuts["top_mass_sig"] = Cuts.top_mass_sig
from project_histos import Cuts, filter_alnum, Sample from make_plots import plot_hists, canvas_margin, legend import argparse from collections import OrderedDict as dict if __name__ == "__main__": parser = argparse.ArgumentParser() parser.add_argument("-i", "--infile", type=str, required=True) parser.add_argument("--normalize_lumi", action="store_true", default=True) args = parser.parse_args() # fname = "/Users/joosep/Documents/stpol/data/step3_trees_Mar28/T_t.root" sample = Sample.fromFile(args.infile) # ev_ids = dict() # for ev in sample.tree: # evid = int(ev.event_id) # if evid not in ev_ids.keys(): # ev_ids[evid] = 1 # else: # ev_ids[evid] += 1 cuts = dict() cuts["MT_mu"] = Cuts.mt_mu cuts["2J"] = Cuts.n_jets(2) cuts["1T"] = Cuts.n_tags(1) # cuts["eta_lj"] = Cuts.eta_lj cuts["rms_lj"] = Cuts.rms_lj cuts["eta_jet"] = Cuts.eta_jet cuts["top_mass_sig"] = Cuts.top_mass_sig