def eta(obj,stepsize,eta_low,eta_high): # reads the etaBits in twos_complement.
# The discommented lines transform etaBits to the physical values for eta and print a warning if eta is out of physical senseful boundaries
obj.etaBits = twos_complement_sign(obj.etaBits,9)
# print "len(hw_list) = ", len(hw_list), ", corresponds to ", len(hw_list)/8, " events"
# print "len(emu_out_list) = ", len(emu_out_list), ", corresponds to ", len(emu_out_list)/8, " events"
hist_input_dict = {}
for var in hist_parameters:
hpv = hist_parameters[var]
hist_input_dict[var] = ROOT.TH1D("input_"+var+"_{f}".format(f=filename), "hw_input_"+var+" [{f}]".format(f=filename), int(hpv[1]), hpv[2], hpv[3])
hist_input_dict[var].SetXTitle(hpv[0])
hist_input_dict[var].SetYTitle("N")
for i in xrange(len(in_muons)):
hist_input_dict["qualityBits"].Fill(in_muons[i].input_qualityBits)
hist_input_dict["ptBits"].Fill(in_muons[i].input_ptBits)
hist_input_dict["phiBits"].Fill(in_muons[i].input_phiBits)
hist_input_dict["etaBits"].Fill(twos_complement_sign(in_muons[i].input_etaBits))
#######
hist_dict = {}
for var in hist_parameters:
hpv = hist_parameters[var]
hist_dict[var] = TH1D("output_"+var+"_{f}".format(f=filename), "hw_output_"+var+" [{f}]".format(f=filename), int(hpv[1]), hpv[2], hpv[3])
hist_dict[var].SetXTitle(hpv[0])
hist_dict[var].SetYTitle("N")
for i in xrange(len(hw_list)):
hist_dict["qualityBits"].Fill(hw_list[i].qualityBits)
hist_dict["ptBits"].Fill(hw_list[i].ptBits)
hist_dict["phiBits"].Fill(hw_list[i].phiBits)
hist_dict["etaBits"].Fill(hw_list[i].etaBits)
