10 from __future__
import division
16 output = abs(phi1 - phi2)
18 output = math.pi * 2. - output
22 return (
delta_phi(phi1, phi2)**2. + (eta1 - eta2)**2.)**0.5
26 self.
tfile = ROOT.TFile(filename)
30 for each_entry
in self.
tree:
32 entry_dict[
"truth_e_ET"] = np.asarray(each_entry.truth_e_ET)/1000.
33 entry_dict[
"truth_e_phi"] = np.asarray(each_entry.truth_e_phi)
34 entry_dict[
"truth_e_eta"] = np.asarray(each_entry.truth_e_eta)
35 entry_dict[
"eg_eta"] = np.asarray(each_entry.eg_eta) * 0.1
36 entry_dict[
"eg_phi"] = np.asarray(each_entry.eg_phi) * math.pi/32.
37 entry_dict[
"eg_ET"] = np.asarray(each_entry.eg_ET)/1000.
38 entry_dict[
"eg_wstotnum"] = np.asarray(each_entry.eg_wstotnum)
39 entry_dict[
"eg_wstotden"] = np.asarray(each_entry.eg_wstotden)
40 entry_dict[
"eg_rhadnum"] = np.asarray(each_entry.eg_wstotnum)
41 entry_dict[
"eg_rhadden"] = np.asarray(each_entry.eg_wstotden)
42 entry_dict[
"eg_retanum"] = np.asarray(each_entry.eg_wstotnum)
43 entry_dict[
"eg_retaden"] = np.asarray(each_entry.eg_wstotden)
44 entry_dict[
"eg_haveseed"] = np.asarray(each_entry.eg_haveseed).astype(dtype=bool)
45 entry_dict[
"eg_nTOBs"] =
int(each_entry.eg_nTOBs)
47 entry_dict[
"truth_tauvisible_ET"] = np.asarray(each_entry.truth_tauvisible_ET)/1000.
48 entry_dict[
"truth_tauvisible_phi"] = np.asarray(each_entry.truth_tauvisible_phi)
49 entry_dict[
"truth_tauvisible_eta"] = np.asarray(each_entry.truth_tauvisible_eta)
51 entry_dict[
"eg_tauOregon_Et"] = np.asarray(each_entry.eg_tauOregon_Et)/1000.
52 entry_dict[
"eg_tauOregon_Iso"] = np.asarray(each_entry.eg_tauOregon_Iso)
53 entry_dict[
"eg_tauTLV_Et"] = np.asarray(each_entry.eg_tauTLV_Et)/1000.
54 entry_dict[
"eg_tauTLV_Iso"] = np.asarray(each_entry.eg_tauTLV_Iso)
66 for key, value
in inputs.items():
67 setattr(self, key, value)
75 for i_e
in range(len(self.truth_e_ET)):
77 for i_tob
in range(self.eg_nTOBs):
79 if self.eg_ET[i_tob] < 15.:
81 if delta_R(self.truth_e_eta[i_e], self.truth_e_phi[i_e], self.eg_eta[i_tob], self.eg_phi[i_tob]) > 0.12:
83 tob_id_tem.append(i_tob)
84 output.append(tob_id_tem)
86 for i_e
in range(len(self.truth_tauvisible_ET)):
88 for i_tob
in range(self.eg_nTOBs):
92 if delta_R(self.truth_tauvisible_eta[i_e], self.truth_tauvisible_phi[i_e], self.eg_eta[i_tob], self.eg_phi[i_tob]) > 0.12:
94 tob_id_tem.append(i_tob)
95 output.append(tob_id_tem)
100 return np.argmax(self.truth_e_ET)
102 if len(self.truth_tauvisible_ET) == 0:
104 return np.argmax(self.truth_tauvisible_ET)
109 return (abs(self.truth_e_eta[i_e]) > 1.375
and abs(self.truth_e_eta[i_e]) < 1.52)
111 return (abs(self.truth_tauvisible_eta[i_e]) > 1.375
and abs(self.truth_tauvisible_eta[i_e]) < 1.52)