ATLAS Offline Software
Functions | Variables
maketauturnon Namespace Reference

Functions

def rebin (binning, data)
 
def rebin2 (binning, data1, data2)
 
def turn_on_curve (all_data, selected_data, bins)
 
def main (path)
 

Variables

 inputfile
 

Function Documentation

◆ main()

def maketauturnon.main (   path)

Definition at line 52 of file maketauturnon.py.

52 def main(path):
53  myntuple = eFEXNTupleLoader(path)
54  all_l_ET = []
55  selected_l_ET_TLV = []
56  selected_l_eta_TLV = []
57  selected_l_ET_Oregon = []
58  selected_l_eta_Oregon = []
59  all_l_eta = []
60  for each in myntuple.entries():
61  each.set_class_type("tau")
62  if each.leading_l_id() is None:
63  continue
64  # if len(each.truth_l_ET)==0:
65  # continue
66  if each.is_in_crack(each.leading_l_id()):
67  continue
68  if abs(each.truth_tauvisible_eta[each.leading_l_id()]) > 2.3:
69  continue
70  # if abs(each.truth_tauvisible_ET[each.leading_l_id()] <= 40):
71  # continue
72  all_l_ET.append(each.truth_tauvisible_ET[each.leading_l_id()])
73  all_l_eta.append(each.truth_tauvisible_eta[each.leading_l_id()])
74  matched_tobs = each.get_eg_matchedtob_id()[each.leading_l_id()]
75  tlv_selected = False
76  oregon_selected = False
77  for i in matched_tobs:
78  if each.eg_tauTLV_Et[i] <= 12:
79  continue
80  if each.eg_tauTLV_Iso[i] < 0.66 and each.eg_tauTLV_Et[i] <= 15:
81  continue
82  tlv_selected = True
83  break
84  for i in matched_tobs:
85  if each.eg_tauOregon_Et[i] <= 12:
86  continue
87  if each.eg_tauOregon_Iso[i] < 0.66 and each.eg_tauOregon_Et[i] <= 15:
88  continue
89  oregon_selected = True
90  break
91  if tlv_selected:
92  selected_l_ET_TLV.append(each.truth_tauvisible_ET[each.leading_l_id()])
93  #selected_l_eta_TLV.append(each.truth_tauvisible_eta[each.leading_l_id()])
94  if oregon_selected:
95  selected_l_ET_Oregon.append(each.truth_tauvisible_ET[each.leading_l_id()])
96  #selected_l_eta_Oregon.append(each.truth_tauvisible_eta[each.leading_l_id()])
97  hist_height, bin_centre, hist_error = turn_on_curve(all_l_ET, selected_l_ET_TLV, range(0,100,5))
98  histplot([bin_centre], [hist_height], [hist_error], errorbar_limit=1, filename="tautlvtrunon")
99  hist_height, bin_centre, hist_error = turn_on_curve(all_l_ET, selected_l_ET_Oregon, range(0,100,5))
100  histplot([bin_centre], [hist_height], [hist_error], errorbar_limit=1, filename="tauoregontrunon")
101  # hist_height, bin_centre, hist_error = turn_on_curve(all_l_eta, selected_l_eta, np.linspace(-2.3,2.3,10))
102  # histplot([bin_centre], [hist_height], [hist_error], errorbar_limit=1, filename="eta", xlabel="eta")
103 
104 

◆ rebin()

def maketauturnon.rebin (   binning,
  data 
)

Definition at line 15 of file maketauturnon.py.

15 def rebin(binning, data):
16  hist, edge = np.histogram(data, binning)
17  newbin = [edge[0]]
18  total_i = 0
19  for i in range(len(hist)):
20  total_i += hist[i]
21  if total_i > 0:
22  newbin.append(edge[i+1])
23  total_i = 0
24  newbin[-1] = binning[-1]
25  return newbin
26 

◆ rebin2()

def maketauturnon.rebin2 (   binning,
  data1,
  data2 
)

Definition at line 27 of file maketauturnon.py.

27 def rebin2(binning, data1, data2):
28  hist1, edge = np.histogram(data1, binning)
29  hist2, edge = np.histogram(data2, binning)
30  newbin = [edge[0]]
31  total_1 = 0
32  total_2 = 0
33  for i in range(len(hist1)):
34  total_1 += hist1[i]
35  total_2 += hist2[i]
36  if total_1 > 0 and total_2 > 0:
37  newbin.append(edge[i+1])
38  total_1 = 0
39  total_2 = 0
40  newbin[-1] = binning[-1]
41  return newbin
42 

◆ turn_on_curve()

def maketauturnon.turn_on_curve (   all_data,
  selected_data,
  bins 
)

Definition at line 43 of file maketauturnon.py.

43 def turn_on_curve(all_data, selected_data, bins):
44  bins = rebin(bins, all_data)
45  hist_all, edges = np.histogram(all_data, bins)
46  hist_selected, edges = np.histogram(selected_data, bins)
47  bin_centre = [(edges[i] + edges[i+1])/2. for i in range(len(edges)-1)]
48  hist_height = hist_selected / hist_all
49  hist_error = ((hist_selected**0.5 / hist_all)**2. + (hist_all**0.5 * hist_selected / hist_all**2.)**2.)**0.5
50  return (hist_height, bin_centre, hist_error)
51 

Variable Documentation

◆ inputfile

maketauturnon.inputfile

Definition at line 114 of file maketauturnon.py.

maketauturnon.rebin
def rebin(binning, data)
Definition: maketauturnon.py:15
maketauturnon.rebin2
def rebin2(binning, data1, data2)
Definition: maketauturnon.py:27
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:194
maketauturnon.turn_on_curve
def turn_on_curve(all_data, selected_data, bins)
Definition: maketauturnon.py:43
plotlib.histplot
def histplot(data_x, data_y, error=None, label=None, **kwargs)
Definition: plotlib.py:18
maketauturnon.main
def main(path)
Definition: maketauturnon.py:52