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)

            hists_input[var].GetYaxis().SetRangeUser(0, 1.1*c)
            hists_input[var].GetYaxis().SetTitle("N")
            
            txt = TLatex(0.48, 0.87, "N input events: {n}".format(n=input_parser.get_n_valid()/6))
            set_text_style(txt)
            
            hw_leg.AddEntry(hists_input[var], "RX content (N={tot})".format(tot=hists_input[var].Integral()), "f")
            hw_leg.AddEntry(hists_imd[var], "TX: intermediate (N={tot})".format(tot=hists_imd[var].Integral()), "f")
            hw_leg.AddEntry(hists_output[var], "TX: final (N={tot})".format(tot=hists_output[var].Integral()), "P")

            hists_input[var].Draw()
            hists_output[var].Draw("textPsame")
            hw_leg.Draw()
            hists_imd[var].Draw("same")

            txt.Draw()
            canvas.Print("{f}/figures/hw_{name}.pdf".format(f=fnames['base'], name=hist_parameters[var][0]))
            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)

            hists_input[var].GetYaxis().SetRangeUser(0, 1.1 * c)
            hists_input[var].GetYaxis().SetTitle("N")

            txt = TLatex(
                0.48, 0.87,
                "N input events: {n}".format(n=input_parser.get_n_valid() / 6))
            set_text_style(txt)

            hw_leg.AddEntry(
                hists_input[var],
                "RX content (N={tot})".format(tot=hists_input[var].Integral()),
                "f")
            hw_leg.AddEntry(
                hists_imd[var], "TX: intermediate (N={tot})".format(
                    tot=hists_imd[var].Integral()), "f")
            hw_leg.AddEntry(
                hists_output[var],
                "TX: final (N={tot})".format(tot=hists_output[var].Integral()),
                "P")

            hists_input[var].Draw()