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)