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
