_log.debug("-- final muons: --")
for out_mu in out_muons[cntr:cntr+8]:
_log.debug("out rank (pT+Q): {r}".format(r=out_mu.ptBits+out_mu.qualityBits))
_log.debug("-- inter muons: --")
for i, imd_mu in enumerate(intermediate_muons[event*24:(event+1)*24]):
if i == 0: _log.debug("-"*5+"EMTF-")
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)
intermediate_muons[event * 24:(event + 1) * 24]):
if i == 0: _log.debug("-" * 5 + "EMTF-")
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(
if not options.nodebug:
intermediate_muons = output_parser.get_intermediate_muons()
imd_mu_non_zero = [imd_mu for imd_mu in intermediate_muons if imd_mu.bitword != 0 ]
ranks = output_parser.get_ranks()
rank_num_of_non_zeros = 0
for i in xrange(len(ranks)):
if ranks[i]!=0:
rank_num_of_non_zeros = rank_num_of_non_zeros+1
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")
intermediate_muons = output_parser.get_intermediate_muons()
imd_mu_non_zero = [
imd_mu for imd_mu in intermediate_muons if imd_mu.bitword != 0
]
ranks = output_parser.get_ranks()
rank_num_of_non_zeros = 0
for i in xrange(len(ranks)):
if ranks[i] != 0:
rank_num_of_non_zeros = rank_num_of_non_zeros + 1
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"
