79 vec_comb_err =
array(
'd')
82 grl = pt.get_grl(year)
85 livetime, zlumi, zerr, olumi, timestamp, dfz_small = pt.get_dfz(args.indir, year, run, channel)
88 if livetime < pt.runlivetimecut:
89 if livetime >= 0.:
print(f
"Skip Run {run} because of live time {livetime/60:.1f} min")
92 dfz_small[
'CombEff'] = dfz_small[channel +
'EffComb']
93 dfz_small[
'CombErr'] = dfz_small[channel +
'ErrComb']
95 dfz_small.loc[dfz_small[
'CombErr'] == 0.,
'CombErr'] = 1.
99 dfz_small[
'CombEff'] *= dq_cf.correction(dfz_small[
'OffMu'], channel, campaign,
int(run))
100 dfz_small[
'CombErr'] *= dq_cf.correction(dfz_small[
'OffMu'], channel, campaign,
int(run))
102 for index, event
in dfz_small.iterrows():
103 pileup =
int(event.OffMu)
105 weight_comb = 1/
pow(event.CombErr, 2)
106 if pileup
not in dict_mu:
107 dict_mu[pileup] = pileup
108 dict_comb[pileup] = weight_comb * event.CombEff
109 dict_comb_err[pileup] = weight_comb
111 dict_comb[pileup] += weight_comb * event.CombEff
112 dict_comb_err[pileup] += weight_comb
115 print(
"File has no filled lumi blocks!")
118 for pileup
in dict_mu:
119 comb_weighted_average = dict_comb[pileup]/dict_comb_err[pileup]
120 comb_error = sqrt(1/dict_comb_err[pileup])
121 vec_comb.append(comb_weighted_average)
122 vec_comb_err.append(comb_error)
124 vec_mu.append(pileup)
126 all_graphs.append((year, R.TGraphErrors(len(vec_comb), vec_mu, vec_comb, R.nullptr, vec_comb_err)))
130 leg = R.TLegend(0.645, 0.7, 0.805, 0.9)
132 ymin, ymax = 0.52, 0.74
133 elif channel ==
"Zmumu":
134 ymin, ymax = 0.74, 0.84
137 for igraph
in range(len(all_graphs)):
138 comb_graph = all_graphs[igraph][1]
139 comb_graph.SetMarkerSize(1)
140 comb_graph.SetMarkerColor(R.kBlack+igraph)
141 comb_graph.SetMarkerStyle(R.kFullCircle+igraph)
143 comb_graph.Draw(
"ap")
144 comb_graph.GetHistogram().SetYTitle(
"#varepsilon_{event}^{"+pt.plotlabel[channel]+
"}#times F^{MC}")
145 comb_graph.GetHistogram().GetYaxis().SetRangeUser(ymin, ymax)
146 comb_graph.GetHistogram().GetXaxis().SetLimits(xmin, xmax)
147 comb_graph.GetHistogram().SetXTitle(
"Pileup (#mu)")
149 comb_graph.Draw(
"p same")
151 leg.SetTextSize(0.07)
152 leg.AddEntry(comb_graph,
"Data 20"+all_graphs[igraph][0],
"ep")
156 pt.drawAtlasLabel(0.6, ymax-0.46,
"Internal")
157 pt.drawText(0.2, ymax-0.46, date_tag)
158 pt.drawText(0.2, ymax-0.52, pt.plotlabel[channel] +
" counting")
159 elif channel ==
"Zmumu":
160 pt.drawAtlasLabel(0.6, ymax-0.56,
"Internal")
161 pt.drawText(0.2, ymax-0.62, pt.plotlabel[channel] +
" counting")
163 c1.SaveAs(outdir + channel +
"_eventeff_vs_mu_"+out_tag+
".pdf")