if i == 4: _log.debug("-"*5+"EMTF+")
if i == 8: _log.debug("-"*5+"OMTF-")
if i == 12: _log.debug("-"*5+"OMTF+")
if i == 16: _log.debug("-"*5+"BMTF")
_log.debug("imd rank (pT+Q): {r} ({pt}+{q}={calc})".format(r=ranks[event*24+i], pt=imd_mu.ptBits, q=imd_mu.qualityBits, calc=imd_mu.ptBits+imd_mu.qualityBits))
cntr += 8
hist_rnk = create_and_fill_rank_hist(ranks, pattern)
plot_modifier(hist_rnk, "rank", "N", ROOT.kBlack)
hist_rnk.Draw()
canvas.Print("{f}/figures/hw_rank.pdf".format(f=fnames['base']))
hists_input = create_and_fill_muon_hists(hist_parameters, in_muons, pattern+"in")
hists_output = create_and_fill_muon_hists(hist_parameters, out_muons, pattern+"out")
hists_imd = create_and_fill_muon_hists(hist_parameters, intermediate_muons, pattern+"imd")
for var in hist_parameters:
hw_leg = TLegend(0.4, 0.7, 0.7, 0.85)
set_legend_style(hw_leg)
a = hists_input[var].GetBinContent(hists_input[var].GetMaximumBin())
b = hists_output[var].GetBinContent(hists_output[var].GetMaximumBin())
c = max(a,b)
plot_modifier(hists_input[var], hist_parameters[var][0], "N", ROOT.kAzure-4)
plot_modifier(hists_output[var], hist_parameters[var][0], "N", ROOT.kBlack, 20)
plot_modifier(hists_imd[var], hist_parameters[var][0], "N", ROOT.kBlack)
hists_imd[var].SetFillStyle(0)
hists_imd[var].SetLineStyle(2)
if i == 12: _log.debug("-" * 5 + "OMTF+")
if i == 16: _log.debug("-" * 5 + "BMTF")
_log.debug("imd rank (pT+Q): {r} ({pt}+{q}={calc})".format(
r=ranks[event * 24 + i],
pt=imd_mu.ptBits,
q=imd_mu.qualityBits,
calc=imd_mu.ptBits + imd_mu.qualityBits))
cntr += 8
hist_rnk = create_and_fill_rank_hist(ranks, pattern)
plot_modifier(hist_rnk, "rank", "N", ROOT.kBlack)
hist_rnk.Draw()
canvas.Print("{f}/figures/hw_rank.pdf".format(f=fnames['base']))
hists_input = create_and_fill_muon_hists(hist_parameters, in_muons,
pattern + "in")
hists_output = create_and_fill_muon_hists(hist_parameters, out_muons,
pattern + "out")
hists_imd = create_and_fill_muon_hists(hist_parameters,
intermediate_muons,
pattern + "imd")
for var in hist_parameters:
hw_leg = TLegend(0.4, 0.7, 0.7, 0.85)
set_legend_style(hw_leg)
a = hists_input[var].GetBinContent(
hists_input[var].GetMaximumBin())
b = hists_output[var].GetBinContent(
hists_output[var].GetMaximumBin())
c = max(a, b)
if pattern in phys_patterns:
del emu_imd_list[len(intermediate_muons):]
hist_rnk = create_and_fill_rank_hist(ranks, pattern)
plot_modifier(hist_rnk, "rank", "N", ROOT.kBlack)
hist_rnk.Draw()
canvas.Print("{f}/figures/hw_rank.pdf".format(f=fnames['base']))
print "{fn}_in_muons :".format(fn=pattern), non_zero(in_muons), "/", len(in_muons)
print "{fn}_num of final non-zero Output-Muons: ".format(fn=pattern), non_zero(out_muons), "/", len(out_muons)#,"), corresponds to ", len(out_muons)/8," Events"
if not options.nodebug:
print "{fn}_num of intermediate non-zero Output-Muons: ".format(fn=pattern), non_zero(intermediate_muons), "/" , len(intermediate_muons)#, "), corresponds to ", len(intermediate_muons)/24," Events"
print "{fn}_n_ranks".format(fn=pattern), rank_num_of_non_zeros, "/", len(ranks)
hists_input = create_and_fill_muon_hists(hist_parameters, in_muons, pattern+"in")
hists_output = create_and_fill_muon_hists(hist_parameters, out_muons, pattern+"out")
if gt_muons:
hists_gt = create_and_fill_muon_hists(hist_parameters, gt_muons, pattern+"gt")
if not options.nodebug:
hists_imd = create_and_fill_muon_hists(hist_parameters, intermediate_muons, pattern+"imd")
hists_emu_imd = create_and_fill_muon_hists(hist_parameters, emu_imd_list, pattern+"emuimd")
hists_emu_output = create_and_fill_muon_hists(hist_parameters, emu_out_list, pattern+"emuout")
hists_emu_input = create_and_fill_muon_hists(hist_parameters, emu_in_list, pattern+"emuin")
for var in hist_parameters:
hw_leg = TLegend(0.4, 0.7, 0.7, 0.85)
set_legend_style(hw_leg)
a = hists_input[var].GetBinContent(hists_input[var].GetMaximumBin())
print "{fn}_in_muons :".format(
fn=pattern), non_zero(in_muons), "/", len(in_muons)
print "{fn}_num of final non-zero Output-Muons: ".format(
fn=pattern), non_zero(out_muons), "/", len(
out_muons) #,"), corresponds to ", len(out_muons)/8," Events"
if not options.nodebug:
print "{fn}_num of intermediate non-zero Output-Muons: ".format(
fn=pattern
), non_zero(intermediate_muons), "/", len(
intermediate_muons
) #, "), corresponds to ", len(intermediate_muons)/24," Events"
print "{fn}_n_ranks".format(
fn=pattern), rank_num_of_non_zeros, "/", len(ranks)
hists_input = create_and_fill_muon_hists(hist_parameters, in_muons,
pattern + "in")
hists_output = create_and_fill_muon_hists(hist_parameters, out_muons,
pattern + "out")
if gt_muons:
hists_gt = create_and_fill_muon_hists(hist_parameters, gt_muons,
pattern + "gt")
if not options.nodebug:
hists_imd = create_and_fill_muon_hists(hist_parameters,
intermediate_muons,
pattern + "imd")
hists_emu_imd = create_and_fill_muon_hists(hist_parameters,
emu_imd_list,
pattern + "emuimd")
hists_emu_output = create_and_fill_muon_hists(hist_parameters,