def do_all(files):
for f in files:
ph_cands = make_cut_histogram(f.photon_counts)
events = make_cut_histogram(f.event_counts)
#("oq", "grl", "pv", "pt_gt100").project_out().
print filename
print "\n".join(" %20s %12s %12s" % (a.title, comma_num(a.true), comma_num(a.total)) for a in ph_cands.axes_objs)
def do_cutflow_one_file(f, cuts, options):
ph_cands = make_cut_histogram(f.photon_counts)
#ph_cands = ph_cands.isConv.project_out(True)
if options.apply_cut:
cut, value = options.apply_cut.split(":")
ph_cands = getattr(ph_cands, cut).project_out(int(value))
rows = []
for axis in [None] + ordered_axes(ph_cands.axes_objs, cuts):
rows.append(do_cutflow_one_cut(ph_cands, axis, cuts, options))
return rows
